Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Halibut mitochondrial genomes contain extensive heteroplasmic tandem repeat arrays involved in DNA recombination

Kenneth A Mjelle1, Bård O Karlsen1, Tor E Jørgensen1, Truls Moum1 and Steinar D Johansen12*

Author Affiliations

1 Department of Fisheries and Natural Sciences, Bodø University College, N-8049 Bodø, Norway

2 Department of Molecular Biotechnology, Institute of Medical Biology, University of Tromsø, N-9037 Tromsø, Norway

For all author emails, please log on.

BMC Genomics 2008, 9:10  doi:10.1186/1471-2164-9-10


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/9/10


Received:6 July 2007
Accepted:11 January 2008
Published:11 January 2008

© 2008 Mjelle et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Halibuts are commercially important flatfish species confined to the North Pacific and North Atlantic Oceans. We have determined the complete mitochondrial genome sequences of four specimens each of Atlantic halibut (Hippoglossus hippoglossus), Pacific halibut (Hippoglossus stenolepis) and Greenland halibut (Reinhardtius hippoglossoides), and assessed the nucleotide variability within and between species.

Results

About 100 variable positions were identified within the four specimens in each halibut species, with the control regions as the most variable parts of the genomes (10 times that of the mitochondrial ribosomal DNA). Due to tandem repeat arrays, the control regions have unusually large sizes compared to most vertebrate mtDNAs. The arrays are highly heteroplasmic in size and consist mainly of different variants of a 61-bp motif. Halibut mitochondrial genomes lacking arrays were also detected.

Conclusion

The complexity, distribution, and biological role of the heteroplasmic tandem repeat arrays in halibut mitochondrial control regions are discussed. We conclude that the most plausible explanation for array maintenance includes both the slipped-strand mispairing and DNA recombination mechanisms.

Background

Halibuts (family Pleuronectidae) represent the largest of the flatfish species. Whereas Atlantic halibut (Hippoglossus hippoglossus) and Pacific halibut (Hippoglossus stenolepis) are endemic species confined to the North Atlantic and North Pacific Oceans, respectively, the Greenland halibut (Reinhardtius hippoglossoides) has an Arctic-boreal distribution in both the Atlantic and Pacific Oceans. All three species are commercially important flatfishes with extensive annual catch volumes, and the Atlantic halibut has further become increasingly popular in North European aquaculture [1]. Phylogenetic analysis based on partial mitochondrial DNA (mtDNA) sequences supports a sister taxa affiliation of the Hippoglossus and Reinhardtius halibuts among the Pleuronectidae [2].

Genetic markers have been developed to investigate and assess genetic issues within e.g. taxonomy, systematics, conservation biology, population structuring, or breeding programs. Mitochondrial DNA (mtDNA) has become one of the most popular genetic markers [3,4] due to its small size and stable organization, its simple inheritance pattern (maternal without apparent DNA recombination), high copy number, and elevated mutation rate compared to single-copy nuclear DNA. Vertebrate mtDNA is usually less than 17 kb in size with a plasmid-like organization, encoding only 37 gene products (13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes) as well as a main control region (CR) containing transcriptional promoters, at least one of the replication origins as well as the displacement loop (D-loop) [5].

Mitogenomics has been developed to increase the resolution of mtDNA markers by including the complete mitochondrial genome sequence in the analyses. Recently, several genetic issues in bony fishes have been successfully investigated and resolved by mitogenomic analyses, e.g. higher-order taxonomy [6,7], within-family taxonomy [8,9], within-genus taxonomy [10,11], and intraspecific variability among geographically separated populations [11,12]. However, vertebrate mtDNA has some limitations and possible shortcomings as a molecular marker that are important to be aware of, and to further investigate [13,14]. Occasional biparental inheritance has been reported, which challenges the clonal maternal nature of vertebrate mtDNA, and some of the best-known examples are found in human and mice [15,16]. Mitochondrial DNA recombination, sometimes recognizable as a consequence of biparental inheritance, appears more frequently than originally assumed but is still a rare event in vertebrates [17]. Here, heteroplasmic tandem repeat (HTR) arrays in the CR may change due to DNA recombination, with some notable examples reported from bony fishes [18-20]. Finally, mtDNA is not always a strictly neutral marker, and both direct and indirect selection has been noted [13]. ATPase6 gene variation in humans [21] and the inherited bacterial symbionts in arthropods [22] represent fascinating examples of mtDNA selection.

In the present study we have assessed the nucleotide variability in halibut mitochondrial genomes within and between species. The complete mitochondrial genome sequences from four individuals each of the Atlantic-, Pacific-, and Greenland halibuts were determined and analysed. A complex organized HTR array in the mitochondrial CR was discovered and investigated in further detail. These composite arrays provide new evidence of DNA recombination in vertebrate mitochondria.

Results

Gene content and organization of halibut mitochondrial genomes

The complete mitochondrial genome sequences were determined for four individuals each of Atlantic halibut (H. hippoglossus), Pacific halibut (H. stenolepsis), and Greenland halibut (R. hippoglossoides) (Table 1). The circular mtDNAs were identical in gene content (13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes) and organization compared to most vertebrates, but varied in size between 17.546 kb and 18.139 kb (Figure 1; Table 1), which are about 1 kb larger than most bony fish mitochondrial genomes [23]. However, these sizes are not absolute since HTRs are observed within the mitochondrial CR (see below). The GC contents of the mitochondrial genomes were 46.1%, 45.7%, and 45.1% for Atlantic-, Pacific, and Greenland halibuts, respectively. These values are similar to most sequenced bony fish mitochondrial genomes. Furthermore, the codon usage was found to be very similar among the three halibut species investigated and with general discrimination against G at the third nucleotide position.

thumbnailFigure 1. Gene content and organization of halibut mitochondrial genomes. Circular gene map representing the mtDNAs of Atlantic-, Pacific-, and Greenland halibuts. All genes, except ND6 and eight of the transfer RNA genes (indicated by the standard one-letter symbols for amino acids), are encoded by the H-strand. Protein genes and ribosomal RNA genes are indicated by blue and yellow boxes, respectively. The tRNA genes are indicated by red bars, and the control regions in grey boxes. Abbreviations: SSU and LSU, mitochondrial small- and large-subunit ribosomal RNA genes; ND1-6, NADH dehydrogenase subunit 1 to 6; COI-III, cytochrome c oxidase subunit I to III; A6 and A8, ATPase subunit 6 and 8; Cyt b, cytochrome b; oriH and oriL, origin of H-strand and L-strand replication; CR, control region containing the D-loop. The mitochondrial genome sizes of the 12 sequenced specimens, representing three different halibut species, are indicated below the map.

Table 1. Key features of halibut specimens and complete mtDNA sequences

The mitochondrial control regions contain heteroplasmic tandem repeat arrays

Intergenic regions are practically lacking in halibut mtDNAs, except for the short spacer between the tRNA-Asn and tRNA-Cys genes that contains the origin of light-strand (oriL) replication, and the major CR (contains the D-loop) located between the tRNA-Pro and tRNA-Phe genes (Figure 1). Whereas the former region is completely conserved in sequence between the 12 analysed halibut specimens, the latter is more variable and contains control elements like the origin of heavy-strand (oriH) replication and the transcriptional promoters.

A schematic presentation of the CR organization is shown in Figure 2A. All specimens contain extensive direct repeat arrays located between the conserved sequence box (CSB) 3 and the tRNA-Phe gene. These arrays were similar among the halibut species investigated, and consist of a free-standing 11-bp motif flanking each side of the array as well as variable numbers of a 61-bp motif in between. In Atlantic halibut the investigated specimens Hh-1, Hh-2, Hh-3, and Hh-4 contain 12, 13, 19, and 15 copies, respectively, of the 61-bp motif in a plasmid cloned representative of the CRs (Figure 2B). Interestingly, eight different variants of the 61-bp motifs were identified (HTR motifs I to VIII) with a distinct, but scattered distribution pattern within and between specimens (Figure 2B). Similar patterns and distributions were observed both for the Pacific halibut (17–19 HTR copies; Figure 2C) and Greenland halibut (17–21 HTR copies; Figure 2D), but with less complexity compared to the Atlantic halibut.

thumbnailFigure 2. Organization of tandem repeat arrays located within the control region of halibut mtDNAs. (A) Schematic organization of the control region (CR) representing all analysed halibut species and specimens. CR is located between the tRNA genes Pro (P) and Phe (F), and contains the highly conserved termination association sequence (TAS: TACATGTATAA), conserved sequence box D (CSB-D: CCTGGCATTTGGTTCC), pyrimidine sequence run (Py-run: TTCTCTTTTTTTTTTTCCTTTC), and the two conserved sequence boxes associated with oriH (CSB-2: AAACCCCCCTACCCCCC, and CSB-2: TGAAAACCCCCCGGAAACA). The heteroplasmic tandem repeat (HTR) array is located between CSB-3 and tRNA-F gene. (B) Detailed view of the HTR array in cloned and sequenced CR from the Atlantic halibut (H. hippoglossus) specimens. The 61-bp HTR motif was found in 12, 13, 19, and 15 copies in Hh-1, Hh-2, Hh-3, and Hh-4, respectively. The HTR motif starts with an 11-bp submotif (boxed), also found freestanding flanking the 61-bp array motif. Eight different variants of the HTR motif (I-VIII) were noted in Atlantic halibut. Hh-2 has one additional freestanding 11-bp motifs compared to the other specimens. (C) HTR array in Pacific halibut (H. stenolepsis) specimens Hs-1 to Hs-4. Two different variants of the motif were found. (D) HTR array in Greenland halibut (R. hippoglossoides) specimens Rh-1 to Rh-4. Five different variants of the motif were found.

A PCR-amplification and DNA sequencing approach was included to assess the heteroplasmic patterns of the HTR-motif arrays. The HTR region and some flanking sequences were amplified from DNA isolated from all 12 specimens (Figure 3A) and subsequently separated by agarose electrophoresis (Figure 3B). Extensive heteroplasmic features were observed for all 12 specimens, with the most common copy numbers of the 61-bp repeat of about 15–20. The four smallest amplified fragments (named

    a-d
in Figure 3B) from one specimen of Atlantic halibut (Hh-2) were eluted from the gel, cloned into a plasmid vector, and subsequently DNA sequenced. The results are summarized in Figure 3C and confirm that the fragments present in the agarose gel as ladder patterns differs in size by one 61-bp motif. Surprisingly, the smallest fragment (fragment
    d
) lacks a complete 61-bp motif and only contains a few copies of the free-standing 11-bp motif.

thumbnailFigure 3. Heteoplasmic features of HTR arrays in halibut mtDNA. (A) Schematic map of the analysed HTR and flanking regions. The HTR arrays were amplified using a forward primer (FP) and a reverse primer (RP). Here, primer sets L16293/H17347 and L16376/H17369 were used on Atlantic/Pacific halibuts and Greenland halibut, respectively. (B) Separation of amplified products in a 2.5% agarose gel. M, the size marker 1 Kb Plus DNA Ladder from Invitrogen (right). Number of repeat motifs is indicated (left). The four smallest amplified fragments from Hh-2 (marked a, b, c, and d) were eluted from the gel and further sequence analysed. (C) Summary of sequence analyses of fragment

    a-d
from Hh-2. Note that fragment
    a
contains three repeats of the 61-bp motif (I*, II, II), fragment
    b
contains two repeat motifs (II, II), fragment
    c
contains a single motif (II), and fragment
    d
lacks a complete motif. Motif I* is not present in the plasmid cloned array of Hh-2 (Figure 2) and may represent site heteroplasmy at array motif position 33.

Distribution of sequence variation within the halibut mitochondrial genomes

Nucleotide substitutions and deletions were assessed by comparing the complete mtDNA sequence of the 4 specimens of each halibut species. The total numbers of variable sites identified were 105, 103, and 119 in Atlantic-, Pacific-, and Greenland halibuts, respectively. The variable sites include all protein coding and ribosomal RNA genes, the CR, and nine of the 22 transfer RNA genes (Figure 4). Transition substitutions at third codon positions of protein coding genes were the most common changes, and nucleotide deletions were only observed at one site in the Atlantic halibut CR.

thumbnailFigure 4. Variability of halibut mitochondrial genomes. (A) Distribution of variable sites in Atlantic halibut mtDNA numbered according to the sequence of specimen Hh-1 (AM749122; Table 1). The variable sites were aligned to that of Hh-1. Identical sites are indicated by dots and deletions by dashes. Abbreviations are according to the legend to Figure 1. The first and last three 61-bp motifs in the CR heteroplasmic array were included in the alignment. (B) Variable sites in Pacific halibut mtDNA numbered according to the sequence of specimen Hs-1 (AM749126; Table 1). (C) Variable sites in Greenland halibut mtDNA numbered according to the sequence of specimen Rh-1 (AM749130; Table 1).

In order to further evaluate the distribution of sequence variation among the halibut mitochondrial genes, both within and between species, substitution versus nucleotide position was estimated for the different gene regions. This ratio was then divided by the ratio obtained for the mitochondrial ribosomal DNA (mt-rDNA; mtSSU + mtLSU; 2665 bp) in order to estimate a relative number of the mitochondrial gene region variation that could be compared between different species and datasets. Despite the fact that numbers of specimens are too low to perform statistics, some general trends were seen (Table 2). First, protein coding genes were about 3–4 times more variable than mt-rDNA both within and between species. The cytochrome c oxydase (CO) subunits were slightly more conserved than the other mitochondrial encoded subunits. Second, the tRNA gene pool (ca 1550 bp) has a similar sequence variation rate as the mt-rDNA. Finally, whereas the CR possesses a between-species sequence variation similar to that of most of the protein genes, it was clearly elevated within the halibut species. In fact, the CRs in Atlantic-, Pacific, and Greenland halibuts were about 10 times more variable per site than the corresponding mt-rDNAs.

Table 2. Estimates of gene specific variation related to the mitochondrial ribosomal DNA gene region

Based on the knowledge about gene specific variability within halibut mtDNA we selected three regions with moderate to high sequence variation (parts of ND1, COI, and CR) for an extended study of Atlantic halibut. Thirty additional specimens (15 wild caught and 15 farmed progenies) from a halibut hatchery were subjected to mtDNA PCR analysis (Figure 5A) and subsequently DNA sequencing. Here, a 1770 bp region was analysed and compared, and the results are summarized in Figure 5B. When including the 4 completely sequenced specimens (Hh-1 to Hh-4), 24 variable sites were detected resulting in 13 distinct haplotypes (Haplotypes a-m). Interestingly, 15 specimens possess an identical haplotype (Haplotype g; Figure 5B), all corresponding to the farmed progenies that probably reflect siblings, or half siblings. This result supports the potential of mtDNA as molecular marker in breeding programs, stock assessments, or population studies.

thumbnailFigure 5. Mitochondrial haplotype in Atlantic halibut. (A) Amplification strategy of partial ND2, COI, and CR mitochondrial sequences. PCR primers and sequencing primers (Table 2) are indicated below and above the schematic line-drawings, respectively. (B) Summary of haplotypes detected within the ca 1770 bp sequence. Haplotype g refers to sibling or half sibling progenies from halibut hatchery.

Discussion

We have sequenced and compared the complete mitochondrial genome sequences of 4 individuals each of the flatfish species Atlantic halibut, Pacific halibut, and Greenland halibut, all related members of the family Pleuronectidae. The mitochondrial genomes were similar to most other bony fish species, except for an unusual large and complex CR located between the tRNA-Pro and tRNA-Phe genes. Halibut CR contains a HTR array of a 61-bp motif, most frequently present in 15–20 copies of each individual.

The within-species variation in mtDNA includes only about 20–100 sequence positions between the individuals. These numbers correlates well with those observed among 12 individuals of Theragra pollocks [11,12]. The variable sites are not equally distributed along the mitochondrial genome sequence, with the structural RNA genes as the most conserved sequence regions (Table 2). The latter observation is best explained by the complex structural constrains of their corresponding tRNA and rRNAs due to secondary and tertiary RNA:RNA interactions, as well as RNA:protein interactions. Interestingly, the structure determination of the vertebrate mitochondrial ribosome explains some of the dramatic reduction in size of the mitochondrial rRNAs, which leaves almost exclusively the highly conserved regions involved in ribosome function and ribosomal protein binding [24].

The elevated within-species sequence variation in CR observed in all three halibut species (Table 2) appears unique compared to other investigated fish mtDNA genomes. The only fish species were complete mtDNA sequences have been recovered from multiple specimens (12 individuals) is the Theragra pollocks [11,12]. Intraspecific sequence variability estimates were similar to that of the halibut species, but with a notable exception of the CR. The Theragra CR showed variability similar to that of the protein coding genes, an observation significantly different from that of halibuts (Table 2). The variable sites in halibuts are almost exclusively located in the extended termination associated sequence (ETAS) and CSB regions located at the 5' end and 3' end of the CR, respectively. What molecular processes that causes this elevated sequence variability is currently not known, but DNA recombination events at HTR arrays (see below) are likely to be involved.

HTRs in mitochondrial CR are widespread, but scattered among vertebrates [18]. Five different locations within the mitochondrial CR have been noted to harbour HTRs [25]. Whereas the RS1 and RS2 sites are located at the CR 5' end in proximity to the termination association sequence, RS3 to RS5 are located close to the oriH replication at the 3' end of the CR. Thus, the presence of HTR in CR is probably associated with the DNA replication processes in vertebrate mtDNA [13]. The complexity of HTR motifs vary greatly among different vertebrates, from simple di- and tetra-nucleotide microsatellites to motifs more than 150 bp in length and at high copy numbers [26-28].

The halibut 61-bp motif HTR array is located at site RS5 between CSB-3 and the 3' end of CR (Figure 2A), and differs from the RS1 HTR arrays seen in e.g. Atlantic cod and Asian arowana that consist of only 2–6 copies of approximately 40-bp motifs [28-30]. The RS5 HTR is a conserved feature among the Pleuronectidae where a ca 60-bp motif array is present in e.g. Spotted halibut (Verasper variegatus; DQ403797), Barfin flounder (V. moseri; EF025506), Winter flounder (Pseudopleuronectes americanus), Yellowtail flounder (Limanda ferruginea), and American plaice (Hippologlossoides platessoides) [31], in addition to the three halibut species investigated in this study. Interestingly, RS5 direct repeats are also noted in Soleidae, but these appear unrelated in sequence to the Pleuronectidae HTRs and do not create heteroplasmy in mtDNA [32]. Partial sequencing of the mitochondrial CR in European flounder (Platichthys flesus) identified a different repeat motif at RS1 [19,33]. This 19-bp motif was involved in extensive heteroplasmy identified in a study including 168 individuals [19]. Interestingly, two different types of repeat motifs were noted among the 18 individuals studied in more detail, and one of these contains a compound array consisting of both motif types. Our finding of multiple types of HTR motifs in Atlantic-, Pacific-, and Greenland halibuts represents an extended support of the observation in European flounder. Errors during mtDNA replication (e.g. slipped strand mispairing) [34] cannot fully explain the halibut length heteroplasmy since repeat motifs in arrays of most individuals are not identical (Figure 2). Furthermore, technically generated mutations in the sequences, as well as the possibility that an ancestral sequence variant that contained all motif variants, are both highly unlikely explanations since the same type of motifs appear in more than one species and that eight different types were present among the four individuals of Atlantic halibut. Thus, we strongly favour DNA recombination as the most plausible mechanism, a conclusion supporting the findings of Hoarau and co-workers in European flounder mtDNA [19].

Is there a biological role of the mitochondrial HTR arrays in halibuts? The facts that repeat motifs are highly conserved in sequence both between individuals and between Pleuronectidae species (Figure 2) indicate a functional role in the mitochondria. However, mitochondrial genomes lacking the motif, or with only a single copy present (Figure 3), favour no essential role of the array or the motif sequences. The deletion variant (fragment d in Figure 3C) may represent a dead-end of array heteroplasmy unless the HTR motif is reintroduced by DNA recombination. Interestingly, the deleted region is flanked by identical copies of the 11-bp motif and thus probably is generated by a slipped strand mispairing-like process [34], similar to that reported in mitochondria associated with some human diseases [35]. The HTR arrays in halibuts are located between the putative promoter region (3' end of CR) and oriH, and HTRs in RS5 have been functionally linked to the initiation of mtDNA replication [27]. A role of stable secondary structures of nucleotide repeats nucleotides has been suggested. Such putative structures might act at the RNA or DNA levels [36], but at present no experimental biochemical evidence has been provided to support this notion in mitochondria. To further elucidate the molecular evolution and biological roles of HTR arrays in halibut mitochondrial genomes, investigations of the distribution and variation of arrays among different tissues and at different developmental stages should be performed. Studying array variability of mother and progeny would be of particular interest in order to identify possible DNA recombination events. The well studied example of similar RS5 arrays in mitochondria of European rabbits provides an interesting model system for such analyses [27,37-39].

Conclusion

Unusual molecular features of halibut mitochondrial genomes are located in the control region. Extensive size heteroplasmy was detected in Atlantic-, Pacific-, and Greenland halibut mitochondrial control regions. Heteroplasmic tandem repeat arrays contain different variants of a 61-bp motif in compound organization. We conclude that the most plausible explanation for array maintenance includes both slipped-strand mispairing and DNA recombination mechanisms.

Methods

Fish samples and DNA extraction

Key-features of fish samples and mitochondrial DNA sequences used in this study are listed in Table 1. Of the 30 additional Atlantic halibut specimens obtained from a halibut hatchery at Bodø University College, 15 were wild caught (Northern Norway) and 15 were farmed progenies from the hatchery. DNA was extracted from muscle tissue and fin clip by using the High Pure PCR Template Preparation Kit (Roche).

PCR amplification, cloning, DNA sequencing, and data analysis

Specific primer sets consisting of one heavy (H) and one light (L) strand primer (Additional file 1) were used to amplify the complete halibut mitochondrial genomes in five overlapping fragments (L466/H3978, L3851/H7461, L7109/H10004, L9620/H13706, and L12991/H530). In general, the PCR reactions were performed with the following cycling parameters: 94°C initial denaturation for 3 min, 15 cycles with 94°C denaturation for 60 sec, 48°C annealing for 60 sec, 72°C elongation for 4 min. Then, 15 cycles with 94°C denaturation for 60 sec, 53°C annealing for 60 sec, 72°C elongation for 4 min and finally 72°C for 10 min. Products were run on agarose gels containing ethidium bromide, and bands were excised and purified essentially as previously described [11]. When appropriate, PCR products were inserted into the pCR4-TOPO vector (Invitrogen) and transformed in E. coli competent cells. PCR products were sequenced on both strands by using the BigDye version 3.1 kit (Applied Biosystems) with the same primers as in the PCR and internal primers (Additional file 1). The sequencing products were analysed on an ABI genetic analyser (Applied Biosystems). In general, computer analyses of DNA sequences were performed using software package programs from DNASTAR Inc.

Additional file 1. PCR and DNA sequencing primers. Sequence information and location of DNA primers used in PCR and sequencing reactions

Format: DOC Size: 26KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Authors' contributions

KAM, TEJ, and BOK organized the sequencing of the mitochondrial genomes. SDJ, BOK, KAM, TM, and TEJ contributed to mtDNA sequence analyses. SDJ directed the research in collaboration with BOK and TM. SDJ wrote the paper in collaboration with BOK. All authors read and approved the final manuscript version.

Acknowledgements

We thank Dag H. Coucheron for valuable discussions, Anita Ursvik for technical advices, Arild Eeg (Risør Aquarium, Norway) for providing the Hh-4 specimen of Atlantic halibut, and Michael E. Cobb (Pacific Northwest National Laboratory, Washingthon, USA) and Ruth Withler (Pacific Biological Station, Nanaimo, Canada) for providing the Pacific halibut samples Hs-1 and Hs-2/Hs-3/Hs-4, respectively. This work was supported by grants from Bodø University College and University of Tromsø.

References

  1. Gibson RN: Flatfishes. Biology and Exploitation. In Fish and Aquatic Resources, Series 9. Blackwell Publishing; 2005.

    ISBN 0-632-05926-5

    OpenURL

  2. Sotelo CG, Calo-Mata P, Chalela MJ, Perez-Martin RI, Rehbein H, Hold GL, Russell VJ, Pryde S, Quinteiro J, Izquierdo M, Rey-Mendez M, Rosa C, Santos AT: Identification of flatfish (Pleuronectiforme) species using DNA-based technoques.

    J Agric Food Chem 2001, 49:4562-4569. PubMed Abstract | Publisher Full Text OpenURL

  3. Bruford MW, Bradley DG, Luikart G: DNA markers reveal the complexity of livestock domestication.

    Nature Rev Genet 2003, 4:900-910. Publisher Full Text OpenURL

  4. Waugh J: DNA barcoding in animal species: progress, potential and pitfalls.

    BioEssays 2007, 29:188-197. PubMed Abstract | Publisher Full Text OpenURL

  5. Boore JL: Animal mitochondrial genomes.

    Nucleic Acids Res 1999, 27:1767-1780. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Miya M, Takeshima H, Endo H, Ishiguro NB, Inoue JG, Mukai T, Satoh TP, Yamaguchi M, Kawaguchi A, Mabuchi K, Shirai SM, Nishida M: Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences.

    Mol Phylogenet Evol 2003, 26:121-38. PubMed Abstract | Publisher Full Text OpenURL

  7. Yamanoue Y, Miya M, Matsuura K, Yagishita N, Mabuchi K, Sakai H, Katoh M, Nishida M: Phylogenetic position of tetraodontiform fishes within the higher teleosts: Bayesian inferences based on 44 whole mitochondrial genome sequences.

    Mol Phylogenet Evol 2007, 45:89-101. PubMed Abstract | Publisher Full Text OpenURL

  8. Minegishi Y, Aoyama J, Inoue JG, Miya M, Nishida M, Tsukamoto K: Molecular phylogeny and evolution of the freshwater eels genus Anguilla based on the whole mitochondrial genome sequences.

    Mol Phylogenet Evol 2005, 34:134-146. PubMed Abstract | Publisher Full Text OpenURL

  9. Lavoue S, Miya M, Saitoh K, Ishiguro NB, Nishida M: Phylogenetic relationships among anchovies, sardines, herrings and their relatives (Clupeiformes), inferred from whole mitogenome sequences.

    Mol Phylogenet Evol 2007, 43:1096-1105. PubMed Abstract | Publisher Full Text OpenURL

  10. Doiron S, Bernatchez L, Blier PU: A comparative mitogenomic analysis of the potential adaptive value of Arctic charr mtDNA introgression in brook charr populations (Salvelinus fontinalis Mitchill).

    Mol Biol Evol 2002, 19:1902-1909. PubMed Abstract | Publisher Full Text OpenURL

  11. Ursvik A, Breines R, Christiansen JS, Fevolden S-E, Coucheron DH, Johansen SD: A mitogenomic approach to the taxonomy of pollocks: Theragra chalcogramma and T. finnmarchica represent one single species.

    BMC Evol Biol 2007, 7:87. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  12. Yanagimoto T, Kitamura T, Kobayashi T: Complete nucleotide sequence and variation of mitochondrial DNA from 10 individuals of walleye Pollock, Theragra chalcogramma.

    Fisheries Science 2004, 70:885-895. Publisher Full Text OpenURL

  13. Ballard WO, Whitlock MC: The incomplete natural history of mitochondria.

    Mol Ecol 2004, 13:729-744. PubMed Abstract | Publisher Full Text OpenURL

  14. Barr CM, Neiman M, Taylor DR: Inheritance and recombination of mitochondrial genomes in plants, fungi and animals.

    New Phytologist 2005, 168:39-50. PubMed Abstract | Publisher Full Text OpenURL

  15. Kraytsberg Y, Schwarts M, Brown A, Ebralidse K, Kunz WS, Clayton DA, Vissing J, Khrapko K: Recombination of human mitochondrial DNA.

    Science 2004, 304:981. PubMed Abstract | Publisher Full Text OpenURL

  16. Gyllenstein U, Wharton D, Josefsson A, Wilson AC: Paternal inheritance of mitochondrial DNA in mice.

    Nature 1991, 352:255-257. PubMed Abstract | Publisher Full Text OpenURL

  17. Piganeau G, Gardner M, Eyre-Walker A: A broad survey of recombination in animal mitochondria.

    Mol Biol Evol 2004, 21:2329-2325. Publisher Full Text OpenURL

  18. Lunt DH, Whipple LE, Hyman BC: Mitochondrial DNA variable number tandem repeats (VNTRs): utility and problems in molecular ecology.

    Mol Ecol 1998, 7:1441-1455. PubMed Abstract | Publisher Full Text OpenURL

  19. Hoarau G, Holla S, Lescasse R, Stam WT, Olsen JL: Heteroplasmy and evidence for recombination in the mitochondrial control region of the flatfish Platichthys flesus.

    Mol Biol Evol 2002, 19:2261-2264. PubMed Abstract | Publisher Full Text OpenURL

  20. Tatarenkov A, Avise JC: Rapid concerted evolution in animal mitochondrial DNA.

    Proc R Soc B 2007, 274:1795-1798. PubMed Abstract | Publisher Full Text OpenURL

  21. Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, Brandon M, Easley K, Chen E, Brown MD, Sukernik RI, Olckers A, Wallace DC: Natural selection shaped regional mtDNA variation in human.

    Proc Natl Acad Sci USA 2003, 100:171-176. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Hurst GDD, Jiggins FM: Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effect of inherited symbionts.

    Proc R Soc R 2005, 272:1525-1534. Publisher Full Text OpenURL

  23. National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov/] webcite

  24. Sharma MR, Koc EC, Datta PP, Booth TM, Spremulli LL, Agrawal RK: Structure of the mammalian mitochondrial ribosome reveals an expanded functional role for its component proteins.

    Cell 2003, 115:97-108. PubMed Abstract | Publisher Full Text OpenURL

  25. Savolainen P, Arvestad L, Lundeberg J: mtDNA tandem repeats in domestic dogs and wolves: mutation mechanism studied by analysis of the sequence of imperfect repeats.

    Mol Biol Evol 2000, 17:474-488. PubMed Abstract | Publisher Full Text OpenURL

  26. Berg T, Moum T, Johansen S: Variable numbers of simple tandem repeats make birds of the order Ciconiiformes heteroplasmic in their mitochondrial genomes.

    Curr Genet 1995, 27:257-262. PubMed Abstract | Publisher Full Text OpenURL

  27. Dufresne C, Mignotte F, Gueride M: The presence of tandem repeats and the initiation of replication in rabbit mitochondrial DNA.

    Eur J Biochem 1996, 235:593-600. PubMed Abstract | Publisher Full Text OpenURL

  28. Yue GH, Liew WC, Orban L: The complete mitochondrial genome of a basal teleost, the Asian arowana (Sclerophages formosus, Osteoglossidae).

    BMC Genomics 2006, 7:242. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Arnason E, Rand DM: Heteroplasmy of short tandem repeats in mitochondrial DNA of Atlantic cod, Gadus morhua.

    Genetics 1992, 132:211-220. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Johansen S, Bakke I: The complete mitochondrial DNA sequence of Atlantic cod (Gadus morhua): relevance to taxonomic studies among codfishes.

    Mol Mar Biol Biotechnol 1996, 5:203-214. PubMed Abstract OpenURL

  31. Lee W-J, Conroy J, Howell WH, Kocher TD: Structure and evolution of teleost mitochondrial control region.

    J Mol Evol 1995, 41:54-66. PubMed Abstract | Publisher Full Text OpenURL

  32. Manchado M, Catanese G, Ponce M, Funes V, Infante C: The complete mitochondrial genome of the Senegal sole, Solea senegalensis Kaup. Comparative analysis of tandem repeats in the control region among soles.

    DNA Sequence 2007, 18:169-175. PubMed Abstract | Publisher Full Text OpenURL

  33. Tinti F, Colombari A, Vallisneri M, Piccinetti C, Stagni AM: Comparative analysis of a mitochondrial DNA control region amplified from three Adriatic flatfish species and molecular phylogenesis of Pleuronectiformes.

    Mar Biotechnol 1999, 1:20-24. PubMed Abstract | Publisher Full Text OpenURL

  34. Levinson G, Gutman GA: Slipped-strand mispairing: a major mechanism for DNA sequence evolution.

    Mol Biol Evol 1987, 4:203-221. PubMed Abstract | Publisher Full Text OpenURL

  35. Bianchi NO, Bianchi MS, Richard SM: Mitochondrial genome instability in human cancers.

    Mutat Res 2001, 488(1):9-23. PubMed Abstract | Publisher Full Text OpenURL

  36. Mirkin SM: Expandable DNA repeats and human disease.

    Nature 2007, 447:932-940. PubMed Abstract | Publisher Full Text OpenURL

  37. Casane D, Dennebouy N, deRochambeau H, Mounolou JC, Monnerot M: Genetic analysis of systematic mitochondrial heteroplasmy in rabbits.

    Genetics 1994, 138:471-480. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Casane D, Dennebouy N, deRochambeau H, Mounolou JC, Monnerot M: Nonneutral evolution of tandem repeats in the mitochondrial DNA control region of Lagomorphs.

    Mol Biol Evol 1997, 14:779-789. PubMed Abstract | Publisher Full Text OpenURL

  39. Casane D, Gueride M: Evolution of heteroplasmy at a mitochondrial tandem repeat locus in cultured rabbit cells.

    Curr Genet 2002, 42:66-72. PubMed Abstract | Publisher Full Text OpenURL