Mitochondrial DNA differentiation between two forms of trout Salmo letnica, endemic to the Balkan Lake Ohrid, reflects their reproductive isolation

Mitochondrial haplotype diversity in sympatric populations of Ohrid trout, Salmo letnica was investigated by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of the mtDNA control region and ND1, ND3/4, ND5/6 segments. A 310 bp fragment at the 5' end, and a 340-572 bp fragment at the 3' end of the control region were sequenced from representatives of the populations studied. Based on pairwise comparison of the sequences, five new haplotypes were identified plus one identical with the brown trout Andalusian haplotype from the southern Iberian Peninsula. The combination of both RFLP and sequence data sets yielded a total of 10 composite haplotypes. A high degree of genetic subdivision between S. letnica typicus and S. letnica aestivalis populations was observed. The notion of a sympatric origin for the two morphs is discussed. Length variation of the mtDNA control region due to the presence of an 82 bp unit, tandemly repeated one to four times, in the region between the conserved sequence block-3 (CSB-3) and the gene for phenylalanine tRNA is reported. Further, we demonstrate that a single duplication of the approximately 82 bp repeat unit is a common element of the salmonid mitochondrial control region. The unique genetic structure of Ohrid trout represents a highly valuable genetic resource that deserves appropriate management and conservation.     

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.     

Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.     

Characterization of the red knot (Calidris canutus) mitochondrial control region.

We sequenced the complete mitochondrial control regions of 11 red knots (Calidris canutus). The control region is 1168 bp in length and is flanked by tRNA glutamate (glu) and the gene ND6 at its 5' end and tRNA phenylalanine (phe) and the gene 12S on its 3' end. The sequence possesses conserved sequence blocks F, E, D, C, CSB-1, and the bird similarity box (BSB), as expected for a mitochondrial copy. Flanking tRNA regions show correct secondary structure, and a relative rate test indicated no significant difference between substitution rates in the sequence we obtained versus the known mitochondrial sequence of turnstones (Charadriiformes: Scolopacidae). These characteristics indicate that the sequence is mitochondrial in origin. To confirm this, we sequenced the control region of a single individual using both purified mitochondrial DNA and genomic DNA. The sequences were identical using both methods. The sequence and methods presented in this paper may now serve as a reference for future studies using knot and other avian control regions. Furthermore, the discovery of five variable sites in 11 knots towards the 3' end of the control region, and the variability of this region in contrast to the more conserved central domain in the alignment between knots and other Charadriiformes, highlights the importance of this area as a source of variation for future studies in knots and other birds.     

Structural conservation and variation in the mitochondrial control region of fringilline finches (Fringilla spp.) and the greenfinch (Carduelis chloris)

We sequenced the entire control region and portions of flanking genes (tRNA(Phe), tRNA(Glu), and ND6) in the common chaffinch (Fringilla coelebs), blue chaffinch (F. teydea), brambling (F. montifringilla), and greenfinch (Carduelis chloris). In these finches the control region is similar in length (1,223-1,237 bp) and has the same flanking gene order as in other birds, and contains a putative TAS element and the highly conserved CSB-1 and F, D, and C boxes recognizable in most vertebrates. Cloverleaf-like structures associated with the TAS element at the 5' end and CSB-1 at the 3' end of the control region may be involved with the stop and start of D-loop synthesis, respectively. The pattern of nucleotide and substitution bias is similar to that in other vertebrates, and consequently the finch control region can be subdivided into a central, conserved G-rich domain (domain II) flanked by hypervariable 5'-C-rich (domain I) and 3'-AT-rich (domain III) segments. In pairwise comparisons among finch species, the central domain has unusually low transition/transversion ratios, which suggests that increased G + T content is a functional constraint, possibly for DNA primase efficiency. In finches the relative rates of evolution vary among domains according to a ratio of 4.2 (domain III) to 2.2 (domain I) to 1 (domain II), and extensively among sites within domains I and II. Domain I and III sequences are extremely useful in recovering intraspecific phylogeographic splits between populations in Africa and Europe, Madeira, and a basal lineage in Nefza, Tunisia. Domain II sequences are highly conserved, and are therefore only useful in conjunction with sequences from domains I and III in phylogenetic studies of closely related species.      

Complete sequence of the mitochondrial genome of Odontamblyopus rubicundus (Perciformes: Gobiidae): genome characterization and phylogenetic analysis

Odontamblyopus rubicundus is a species of gobiid fishes, inhabits muddy-bottomed coastal waters. In this paper, the first complete mitochondrial genome sequence of O. rubicundus is reported. The complete mitochondrial genome sequence is 17119 bp in length and contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a control region and an L-strand origin as in other teleosts. Most mitochondrial genes are encoded on H-strand except for ND6 and seven tRNA genes. Some overlaps occur in protein-coding genes and tRNAs ranging from 1 to 7 bp. The possibly nonfunctional L-strand 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. The TAS, CSB-2 and CSB-3 could be detected in the control region. However, in contrast to most of other fishes, the central conserved sequence block domain and the CSB-1 could not be recognized in O. rubicundus, which is consistent with Acanthogobius hasta (Gobiidae). In addition, phylogenetic analyses based on different sequences of species of Gobiidae and different methods showed that the classification of O. rubicundus into Odontamblyopus due to morphology is debatable.     

The complete mitochondrial DNA sequence of the Atlantic salmon, Salmo salar.

The complete sequence of the Atlantic salmon (Salmo salar) mitochondrial genome has been determined. The entire sequence is 16665 base pairs (bp) in length, with a gene content (13 protein-coding, two ribosomal RNA [rRNA] and 22 transfer RNA [tRNA] genes) and order conforming to that observed in most other vertebrates. Base composition and codon usage have been detailed. Nucleotide and derived amino acid sequences of the 13 protein-coding genes from Atlantic salmon have been compared with their counterparts in rainbow trout. A putative structure for the origin of L-strand replication (O(L)) is proposed, and sequence features of the control region (D-loop) are described.     

Molecular cloning and characterization of cDNAs coding for apopolysialoglycoproteins in cherry salmon (Oncorhynchus masou) eggs

Recently we have cloned the cDNAs and genomic DNAs for apopolysialoglycoproteins (apoPSGPs) of Salmo gairdneri (rainbow trout) [Sorimachi, H., Emori, Y., Kawasaki, H., Kitajima, K., Inoue, S., Suzuki, K., & Inoue, Y. (1988) J. Biol. Chem. 262, 17678-17684], and the sequence analyses have indicated that the mRNAs for apoPSGPs vary in length and contain different numbers of identical 39-bp repeating units encoding the tridecapeptide (Asp-Asp-Ala-Thr-Ser-Glu-Ala-Ala-Thr-Gly-Pro-Ser-Gly) as well as highly conserved sequences encoding pre-, pro-, and telo-peptide regions. In this study we isolated cDNA clones for yamame (cherry salmon, river resident form; Oncorhynchus masou ishikawai) apoPSGP using a genomic DNA fragment for rainbow trout apoPSGP as a probe. The nucleotide sequence analyses revealed that the structures of mRNAs for yamame apoPSGP including the noncoding regions are essentially identical to those for rainbow trout, showing 90% sequence identity. Within the repeating region, 4 bp out of the 39 were replaced, producing a different tridecapeptide, Asp-Asp-Ala-Thr-Ser-Glu-Ala-Ala-Thr-Gly-Pro-Ser-Ser. This tridecapeptide is unique to yamame and common among all cDNAs obtained from yamame. Genomic Southern blot analysis showed that the yamame apoPSGP genes constituted a multiple gene family with a similar gene organization to that of rainbow trout. Oligodeoxynucleotide probes (18 bases) synthesized based on specific sequences for the yamame repeating unit hybridized only to the yamame DNA and not to the rainbow trout DNA, and vice versa.(ABSTRACT TRUNCATED AT 250 WORDS)     

The mitochondrial genome structure of the clouded leopard (Neofelis nebulosa).

The complete 16 844 bp mitochondrial genome of Neofelis nebulosa has been sequenced and compared with the complete mitochondrial genomes of Felis catus and the Acinonyx jubatus. The base composition of the mitochondrial genome of N. nebulosa is as follows: A, 5343 bp (31.7%); C, 4441 bp (26.4%); G, 2491 bp (14.8%); T, 4569 bp (27.1%). The genome complement and the gene order of this mitochondrial genome was found to be typical of those reported for other mammals. Several unusual features of this genome, however, were found. First, in protein-coding regions, AT bias in the genome was not prevalent in the third position of codons, as it is in most other mammals, but was found in the second position of codons. Second, in tRNA regions, tRNASer (AGY), which lacked the "DHU" arm, could not be folded into the typical cloverleaf-shaped structure. Third, in the control region, no repetitive sequences (RS)-2 were found. However, RS-2 repetitive motifs usually occurr in the control regions of most great cats. In addition, 4 variable sites were found in CSB-3 of the control region. Fourth, AT content in the control region of the mtDNA from the clouded leopard was lower than it is in other regions.     

The complete nucleotide sequence of the mitochondrial DNA genome of the rainbow trout,Oncorhynchus mykiss

The complete nucleotide sequence of the mitochondrial DNA of the rainbow trout, Onchorynchus mykiss, has been determined. The total length of the molecule is 16,660 bp. The rainbow trout mitochondrial DNA has the same organization described in eutherian mammals, the clawed frog (Xenopus laevis), and the two fish species, Oriental stream loach (Crossotoma lacustre) and carp (Cyprinus carpio). Alignment and comparison of the deduced amino acid sequences of the 13 proteins encoded by rainbow trout and other vertebrate mitochondrial genomes allowed us to estimate that COI is the most conserved mitochondrial subunit (amino acid identity ranging from 85.6% to 94.8%) whereas ATPase 8 is the most variable one (amino acid identity ranging from 30.8% to 70.4%). Putative secondary structures for the 22 tRNAs found in the molecule are given along with an extensive comparison of tRNA sequences among representative species of each major group of vertebrates. In this sense, an unusual cloverleaf structure for the tRNA^Ser(AGY) is proposed. A stem-loop structure inferred for the origin of the L-strand replication (O_L) and the presence of a large polycytidine tract in the O_L loop is described. The existence of this stretch instead of the usual T-rich sequence reported so far in mammal mtDNAs is explained in terms of a less-strict template dependence of the RNA primase involved in the initiation of L-strand replication. Correspondence to: J.M. Bautista     

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.     

Assembly of the mitochondrial membrane system. Sequences of yeast mitochondrial tRNA genes

Two cytoplasmic "petite" (rho-) clones of Saccharomyces cerevisiae have been selected for the retention of the aspartic acid tRNA gene. The two clones, designated DS200/A102 and DS200/A5, have tandemly repeated segments of mitochondrial DNA (mtDNA) with unit lengths of 1,000 and 6,400 base pairs, respectively. The DS200/A102 genome has a single tRNA gene with a 3'-CUG-5' anticodon capable of recognizing the 5'-GAC-3' and 5'-GAU-3' codons for aspartic acid. The mtDNA segment of DS200/A102 has been determined to represent the wild type sequence from 5.3 to 6.8 map units. The genome of DS200/A5 is more complex encompassing the region of wild type mtDNA from 3.5 to 12.7 units. A continuous sequence has been obtained from 3.5 to 8.6 units. In addition to the aspartic acid tRNA, this region codes for the tRNAUGCAla,tRNAUCUArg, tRNAACGArg, tRNAGCUSer,tRNAUCCGly and tRNAUUULys. The DNA sequence of the DS200/A5 genome has allowed us to deduce the secondary structures of the seven tRNAs and to assign precise map positions for their genes. All the tRNAs except tRNA GUCAsp exhibit most of the invariant features of prokaryotic and eukaryotic tRNAs. The aspartic acid tRNA has unusual D and T psi C loops. The structure of this tRNA is similar to the mitochondrial initiator tRNA of Neurospora crassa (Heckman, J.E., Hecker, L.I., Shwartzbach, S.D., Barnett, W.E., Baumstark, B., and RajBhandary, U.L. Cell 13, 83-95).     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Complete sequences of mitochondrial genomes from the Baltic mussel Mytilus trossulus

Marine mussels Mytilus possess two mitochondrial (mt) genomes, which undergo doubly uniparental inheritance (DUI). Female (F) and male (M) genomes are usually highly diverged at the sequence level. Both genomes contain the same set of metazoan genes (for 12 proteins, 2 rRNAs and 23 tRNAs), both lack the atp8 gene and have two tRNAs for methionine. However, recently recombination between those variants has been reported. Both original F and M mt genomes of M. trossulus were replaced by M. edulis mtDNA in the Baltic populations. Highly diverged M genome occurs rarely in the Baltic mussels. Full sequences of the M genome identified in males (sperm) and F genome in females (eggs) were obtained. Both genomes were diverged by 24% in nucleotide sequence, but had similar nucleotide composition and codon usage bias. Constant domain (CD) of the control region (CR), the tRNA and rRNA genes were the most conserved. The most diverged was the variable domain 1 (VD1) of the control region. The F genome was longer than M by 147 bp. and the main difference was localised in the VD1 region. No recombination was observed in whole mtDNA of both studied variants. Nuclear mitochondrial pseudogenes (numts) have not been found by hybridisation with probes complementary to several fragments of the Baltic M. trossulus mtDNA.     

The mitochondrial DNA control region comparison studies of four hinged turtles and its phylogentic significance of the genusCuora sensu lato (Testudinata: Geoemydidae)

The complete sequences of the mitochondrial DNA (mtDNA) control region (CR) ofCistoclemmys flavomarginata, Cistoclemmys galbinifrons, Cuora aurocapitata andCyclemys atripons were amplified by long-polymerase chain reaction (Long-PCR). The lengths were 1207 bp, 1722 bp, 1379 bp and 980 bp, respectively. Combining with the CR sequence ofPyxidea mouhotii (DQ659152), we compared the CR structure, and identified three functional domains (TAS, CD and CSB) in which the conservation sequences (TAS, CSB-F, CSB-1, CSB-2 and CSB-3) were also successfully identified according to their homology to those of other turtles. These 5 turtles have the identical CSB-2 and CSB-3 sequences, and 4 of them have the same CSB-1 sequence while there is one base transversion (T → A) inCy. atripons. We analyzed the variable number of tandem repeat (VNTR) sequences or microsatellites at the 3′ end of CR. The motifs of tandem repeats (7 types) are made up of 2–8 nucleotides, and the copy numbers are from 4 to 48. All of the 5 turtles exceptCy. atripons have the “TATTATAT” repeats and are ended by TA. The results of CR structure analysis displayed that theCuora, Cistoclemmys, andPyxidea have many similarities, but differ fromCyclemys. WithIndotestudo elongate (DQ080043) andIndotestudo forstenii (DQ080044) as outgroups, using the CR sequences (1123bp) excluded the microsatellites at the 3′ end of CR, we constructed the molecular phylogenetic trees using the MP, ML and BI methods. The results showed that there was a strong support to the monophyly of theCuora group consisting ofCuora,Cistoclemmys andPyxidea, which has a close relationship withMauremys andChinemys but far fromCyclemys, which are consistent with the analysis of the CR structure of the 5 turtles.     

Tandemly repeated sequences in mtDNA control region of whitefish, Coregonus lavaretus.

Length variation of the mitochondrial DNA control region was observed with PCR amplification of a sample of 138 whitefish (Coregonus lavaretus). Nucleotide sequences of representative PCR products showed that the variation was due to the presence of an approximately 100-bp motif tandemly repeated two, three, or five times in the region between the conserved sequence block-3 (CSB-3) and the gene for phenylalanine tRNA. This is the first report on the tandem array composed of long repeat units in mitochondrial DNA of salmonids.