An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA
© Jørgensen et al.; licensee BioMed Central Ltd. 2014
Received: 3 May 2014
Accepted: 5 August 2014
Published: 22 August 2014
Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNAThr and tRNAPro genes (“T-P spacers”) have been observed in several taxa, including gadiform species, but information about their biological roles and putative functions is still lacking.
Sequence characterization of the complete European hake Merluccius merluccius mitogenome identified a complex T-P spacer ranging in size from 223–532 bp. Further analyses of 32 gadiform species, representing 8 families and 28 genera, revealed the evolutionary preserved presence of T-P spacers across all taxa. Molecular complexity of the T-P spacers was found to be coherent with the phylogenetic relationships, supporting a common ancestral origin and gain of function during codfish evolution. Intraspecific variation of T-P spacer sequences was assessed in 225 Atlantic cod specimens and revealed 26 haplotypes. Pyrosequencing data representing the mito-transcriptome poly (A) fraction in Atlantic cod identified an abundant H-strand specific long noncoding RNA of about 375 nt. The T-P spacer corresponded to the 5′ part of this transcript, which terminated within the control region in a tail-to-tail configuration with the L-strand specific transcript (the 7S RNA).
The T-P spacer is inferred to be evolutionary preserved in gadiform mitogenomes due to gain of function through a long noncoding RNA. We suggest that the T-P spacer adds stability to the H-strand specific long noncoding RNA by forming stable hairpin structures and additional protein binding sites.
KeywordsAtlantic cod European hake Heteroplasmy lncRNA Mitogenome
One characteristic feature of vertebrate mitochondrial genomes (mitogenomes) is the economical gene organization. Vertebrates usually lack mitochondrial intergenic sequences of any appreciable size other than the approximately 1-kb control region (CR) and the 30-bp origin of light (L) strand replication (OriL). The CR harbors the origin of heavy (H) strand replication (OriH) and initiation sites for H-strand and L-strand transcription, as well as the displacement loop (D-loop), but does not contain any canonical protein coding or RNA coding genes .
Among the 2600 completely sequenced vertebrate mitogenomes, a great majority shares an identical arrangement of 37 genes. However, alternative gene orders due to minor rearrangements and sequence duplications have been noted ,. A common feature of most rearrangements is DNA duplications that involve tRNA genes and intergenic sequences. Mitogenomes of several snake species contain two CR copies ,, and the four gene orders known in birds differ in CR copy number and CR location -. Related rearrangements have also been noted in bony fishes, and include both CR, OriL and tRNA gene duplications -.
A mitogenome intergenic spacer between the tRNAThr and tRNAPro genes (the T-P spacer), located close to the CR, has been reported in a few species representing distantly related vertebrate taxa. Short T-P spacers have been found in ostrich birds , in some toad linages including Xenopus,, and in the mole lizard . Larger T-P spacers have been reported in various species of the salamander family Ambystomatidae ,. These complex spacers vary in size from 240 bp to 680 bp and tandem repeat sequence motifs appear common, but no functional role has yet been assigned. A T-P spacer in a fish species was first discovered in the Atlantic cod (Gadus morhua) mitogenome -, and later noted in the related Walleye pollock (Theragra chalcogramma)  and six additional gadiform species representing three families (Gadidae, Lotidae and Phycidae) . These T-P spacers vary in size from 25 bp in fourbeard rockling (Enchelyopus cimbrius) to 99 bp in haddock (Melanogrammus aeglefinus), and contain one or two copies of a conserved 17-bp motif (Box-motif) . T-P spacers demonstrate interspecific sequence variation and have been applied as genetic markers in population studies of both Walleye pollock  and Atlantic cod . More recently T-P spacers were reported in two additional gadid species . While Greenland cod (Gadus ogac) contains two size variants of the spacer (73 bp and 102 bp), highly similar in sequence to that of Atlantic cod, the Arctic cod (Arctogadus glacialis) T-P spacer appears more complex and variable in size due to short heteroplasmic tandem repeat motifs .
The vertebrate mitochondrial transcriptome (mito-transcriptome) has been investigated in human cells and tissues ,. The H and L strand polycistronic precursor transcripts are processed into the 22 tRNAs, 2 rRNAs, and 11 mRNAs . Here, the L-strand specific promoter (LSP) generates a transcript that gives rise to the ND6 mRNA and 8 tRNAs. The H-strand specific transcription, however, is more complex since two promoters (HSP1 and HSP2) are involved. Whereas the HSP1-specific transcript is short and highly abundant, and is the main source of mitochondrial rRNAs, the HSP2-specific transcript is processed into 8 monocistronic and two bicistronic mRNAs ,. This organization appears highly conserved among vertebrates since Atlantic cod and saithe (Pollachius virens) mito-transcriptomes are very similar to that of humans .
Long noncoding RNAs (lncRNAs), typically longer than 200 nt, are highly abundant in vertebrates where they possess key roles in gene regulation linked to tissue specificity, development and disease ,. The human mitogenome also code for lncRNAs in addition to their 37 canonical genes. The first mitochondrial lncRNA to be discovered was the 7S RNA, transcribed from LSP within the CR ,. Additional mitochondrial lncRNAs have since then been identified in human cells, and include lncND5, lncND6, lncCytB , as well as the prognostic LIPCAR lncRNAs antisense to CytB and COII mRNAs . The lncND5 RNA is of particular interest due to its antisense organization to ND5 mRNA, and the fact that ND5 is the only tightly regulated protein gene in vertebrate mitochondria ,. Mitochondrial lncRNAs appear highly conserved among vertebrates, and we recently reported noncoding RNA corresponding to lncND5 in the mito-transcriptomes of Atlantic cod and saithe .
Here we report two non-overlapping mitochondrial lncRNAs (lncCR-H and lncCR-L) transcribed from opposite strands within the Atlantic cod CR, both terminating at the termination-associated sequence (TAS). While the 500 nt lncCR-L RNA apparently corresponds to the human 7S RNA , the 375 nt lncCR-H RNA has not been reported previously. The 5′ end of lncCR-H corresponds to the T-P spacer that apparently adds stability to the noncoding RNA. The T-P spacer was found to be present in all 32 gadiform species investigated (representing 8 families and 28 genera), and feature analysis indicates directional evolution of spacers from a simple organization in early branching gadiform families (e.g. Moridae or Macrouridae) to complex structures including heteroplasmy in Merluciidae and Gadidae.
Complete mitogenome sequence of the European hake
The phylogenetic position of European hake among Gadiformes was assessed from concatenated mitochondrial protein sequences (3814 amino acid positions) derived from complete mitogenome sequences (Additional file 3: Table S1), and a representative neighbor-joining tree is presented in Additional file 4: Figure S3. Eight gadiform families were included in the analysis and European hake (family Merlucciidae) was found to occupy a more basal position compared to Gadidae and Lotidae, but apparenty more recent than Macrouridae, Bregmacerotidae and Moridae. These observations are in general agreement with previous works on gadiform phylogeny -.
T-P spacer in European hake contains heteroplasmic tandem repeats
The T-P spacer included in the complete European hake mitogenome sequence (Mm1) was determined to 532 bp. This sequence consists of a 34 bp conserved direct repeat containing the 17-bp Box-motif , an optional direct repeat of about 60 bp, and 20 copies of an 8-bp heteroplasmic tandem repeat (HTR) (Figure 1B; Additional file 5: Figure S4A). To investigate individual sequence variations, a T-P spacer region from another specimen (Mm2) was analyzed by PCR amplification, plasmid cloning, and subsequent DNA sequencing. The Mm2 T-P spacer showed extensive length heteroplasmy, varying in size from 258 bp to 430 bp (Figure 1B; Additional file 5: Figure S4B). Size variation was partly caused by the 8-bp HTR motif, which varied from 7 to 23 copies among six sequenced clones. Furthermore, four clones contained an additional 76 bp sequence (Figure 1B). All plasmid clones investigated gave identical sequences, except the HTR copy number and two single nucleotide polymorphic sites (Additional file 5: Figure S4B).
Distribution and variability of T-P spacers among gadiform species
Key features of gadiform T-P spacers
FR751402; This work
Intraspecific variation of the Atlantic cod T-P spacer
Analysis of mito-transcriptome data from Atlantic cod liver tissue
We report a complex intergenic spacer between the tRNAThr and tRNAPro genes in the European hake mitochondrial genome. The T-P spacer contains HTRs, and significant size variation was observed among specimens. Furthermore, we found T-P spacers to be present in mitogenomes of all the 32 gadiform species investigated, representing 8 families and 28 genera, and the majority of spacers harbor a conserved 17-bp Box-motif. Analysis of the Atlantic cod liver mito-transcriptome identified a long noncoding RNA (lncCR-H) that contains the T-P spacer within its 5′ end.
Higher order phylogenetic relationships within Gadiformes are still weakly supported by molecular data -. Other branching patterns than those depicted here would imply instances of secondary loss of the Box-motif, which can of course not be precluded. However, the Box-motif apparently represents a selective feature, of which the exact functional role is still unknown. It is interesting to note that complex T-P spacers in salamanders , contains strikingly similar Box-motifs to those reported in codfishes (Additional file 7: Table S2), suggesting that an evolutionary conserved mitochondrial factor may interact with the motif.
Analysis of the Atlantic cod liver mito-transcriptome identified a relatively abundant lncRNA (lncCR-H) of about 375 nt, that contains the T-P spacer sequence within its 5′ end. The first nucleotide of lncCR-H corresponds to 3 nt downstream the tRNAThr gene, and the 5′ end of the lncRNA was probably generated after RNase Z processing of the tRNA precursor . The 3′ end of lncCR-H appears heterogenous and polyadenylated, and corresponds to sequences in proximity to the TAS box element within the CR. LncCR-H harbours secondary structure elements, including a mirror tRNAPro and tetraloop helices, probably adding stability to the RNA. Interestingly, recent RNA-Seq experiments in breast cell lines identified an expressed human homolog to lncCR-H (our unpublished results), and HeLa cell RNA profiles supported the existence of this mitochondrial lncRNA (see Figure 1B in ). Thus, lncCR-H appears conserved at a wide taxonomic range, which implies a functional role in vertebrate mitochondria.
Mitochondrial lncRNAs corresponding to antisense transcripts of mRNAs have been characterized in human cells ,, as well as in Atlantic cod . In that study we identified a second mitochondrial lncRNA within CR, named lncCR-L. The L-strand specific lncCR-L, which corresponds to the human 7S RNA ,, is apparently initiated at LSP and terminated at TAS. Our finding implicates that TAS is a transcription termination site both for H-strand and L-strand specific lncRNAs. TAS in Atlantic cod consists of a perfect inverted repeat motif (UUAUACAUAUGUAUAA) and thus forms identical tetraloop hairpin structures at the H- and L-strand transcripts. TAS-binding proteins have been reported in human, rat and bovine mitochondria ,. How lncCR-H and lncCR-L are terminated and regulated in Atlantic cod, as well as their biological functions, are currently unknown. However, we speculate that the biological role of the T-P spacer is to add stability to an evolutionary conserved lncRNA where the Box-motif forms an RNA structure signal to be recognized by mitochondrial proteins. lncCR-H may have a role in antisense regulation by premature termination of L-strand transcription at TAS, resulting in a coordination of H and L transcript levels. Anyway, these possibilities need further investigations in different cell and tissue types, and in different species.
Two long noncoding RNA candidates from the mitochondrial control region were described in gadiform fishes. These RNAs were transcribed from different DNA strands, and found to terminate at the TAS region in a tail-to-tail configuration. Whereas the L-strand specific lncRNA (lncCR-L) corresponded to the previously described mitochondrial 7S RNA, lncCR-H was found to be unique by containing a 5′ domain of an intergenic spacer. Analyses of 32 gadiform species revealed that a T-P spacer was present in all investigated species, representing 8 families and 28 genera. Some species contained complex spacer sequences due to conserved motifs, duplications, and heteroplasmy. We propose that the T-P spacer adds stability to lncCR-H by forming stable hairpin structures and additional protein binding sites, but a biological role of lncCR-H is currently unknown.
DNA extraction and PCR amplification
DNA was extracted from muscle tissue samples as previously described ,,,. DNA for complete mitogenome sequencing of European hake was extracted from frozen muscle tissue by using the mtDNA Extractor CT Kit from Wako . The PCR and sequencing primers (Additional file 8: Table S2) applied in the European hake sequencing were designed from the published Atlantic cod mitogenome sequence  and other available gadiform sequences (Additional file 3: Table S1). These primers were used to amplify the mitogenomes in 1 to 4 kb fragments essentially as described by . The primer pair L15498/H15666 was used for amplification of the T-P spacer from most species. For European hake, a specific primer set (L15498/H15642) was designed for T-P spacer amplifications. The PCR products were treated with exonuclease and shrimp alkaline phosphatase (USB) prior to plasmid cloning and sequencing. The amplified fragments were subcloned into the SmaI-site of pUC18 vector using the Sure-Clone Ligation kit (Pharmacia Biotech).
Sanger sequencing of the European hake mitogenome and Atlantic cod T-P spacers was performed directly on PCR products as previously described  using the BigDye kit (Applied Biosystems). Sequencing primers are found in Additional file 8: Table S2. Manual sequencing reactions were performed on most T-P spacers using the Thermo Sequencing Terminator Cycle Sequencing kit (Amersham) and [α-33P] ddNTP (Amersham) as the label. Sequencing reactions were run on 7 M urea/6% polyacrylamide gels after denaturation at 85°C for 2 minutes. Roche 454 pyrosequencing of the Atlantic cod transcriptome was performed as a service given by Eurofins MWG Operon (Germany). A poly(A) enriched normalized liver-specific cDNA library (based on total cellular RNA) was generated by random hexamer first strand synthesis. About 1.2 million reads were obtained from pooling equal amounts of total RNA from 10 Atlantic cod specimens .
A sequence alignment representing complete mitogenomes of 20 gadiform species and one out-group was generated using T-Coffee v/9  with manual refinements. The alignment was build from the complete set of protein codons (except stop codons) creating a concatenated sequence of 11442 nt positions (3814 codons) corresponding to the 13 protein genes. The tree-building methods of neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML) in MEGA version 5  were used to reconstruct molecular phylogenies. NJ trees were build using Jones-Taylor-Thornton model. MP trees were reconstructed using the Close-Neighbor-Interchange search option. ML trees were built from best-fit models of protein evolution generated by MEGA 5 and ProtTest 3 . The topologies of the ME, MP and ML trees were evaluated by bootstrap analyses (2000 replications). Analyses of mito-transcriptome 454 pyrosequencing reads were as previously described .
TEJ and IB participated in DNA sequencing and design of the study, performed analysis, and contributed in the discussions of the results. AU and MA participated in mitogenome sequencing and transcriptome analysis, respectively. TM participated in the design of this study and in discussions of results. SDJ directed the research, participating in the analysis, and wrote the manuscript in collaboration with all the authors. All authors read and approved the final manuscript.
We thank Willy Richardsen, Kees O. Ekeli, Odd M. Smedstad, Wim Vader, Stefan Hugel, Jon Einar Hansen, Bjørn Steinar Sætre and Åse Emblem for kindly providing fish samples. We also thank Ragna Breines, Dag H. Coucheron, Marianne Nymark, Bård O. Karlsen, Åse Emblem and Jarle T. Nordeide for discussions. This work was supported by grants to SDJ from the Norwegian Research Council and the University of Tromsø.
- Wallace DC: Why do we still have a maternally inherited mitochondrial DNA? Insights from evolutionary medicine. Annu Rev Biochem. 2007, 76: 781-821.PubMedView ArticleGoogle Scholar
- Brown C, Dowling TE, Brown WM: Evolution of animal mitochondrial DNA: relevance for population biology and systematics. Annu Rev Ecol Syst. 1987, 18: 269-292.View ArticleGoogle Scholar
- San Mauro D, Gower DJ, Zardoya R, Wilkinson M: A hotspot of gene order rearrangement by tandem duplication and random loss in the vertebrate mitochondrial genome. Mol Biol Evol. 2006, 23: 227-234.PubMedView ArticleGoogle Scholar
- Kumazawa Y, Ota H, Nishida M, Ozawa T: Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster. Mol Biol Evol. 1996, 13: 1242-1254.PubMedView ArticleGoogle Scholar
- Kumazawa Y, Ota H, Nishida M, Ozawa T: The complete nucleotide sequence of a snake (Dinodon semicarinatus) mitochondrial genome with two identical control regions. Genetics. 1998, 150: 313-329.PubMedPubMed CentralGoogle Scholar
- Mindell DP, Sorenson MD, Dimcheff DE: Multiple independent origins of mitochondrial gene order in birds. Proc Natl Acad Sci U S A. 1998, 95: 10693-10697.PubMedPubMed CentralView ArticleGoogle Scholar
- Eberhard JR, Wright TF, Bermingham E: Duplication and concerted evolution of the mitochondrial control region in the parrot genus Amazona. Mol Biol Evol. 2001, 18: 1330-1342.PubMedView ArticleGoogle Scholar
- Abbott CL, Double MC, Trueman JW, Robinson A, Cockburn A: An unusual source of apparent mitochondrial heteroplasmy: duplicate mitochondrial control regions in Thalassarche albatrosses. Mol Ecol. 2005, 14: 3605-3613.PubMedView ArticleGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of Conger myriaster (Teleostei: Anguilliformes): novel gene order for vertebrate mitochondrial genomes and the phylogenetic implications for Auguilliform families. J Mol Evol. 2001, 52: 311-320.PubMedGoogle Scholar
- Satoh TP, Miya M, Endo H, Nishida M: Round and pointed-head grenadier fishes (Actinopterygii: Gadiformes) represent a single sister group: evidence from complete mitochondrial genome sequences. Mol Phylogen Evol. 2006, 40: 129-138.View ArticleGoogle Scholar
- Breines R, Ursvik A, Nymark M, Johansen SD, Coucheron DH: Complete mitochondrial genome sequences of the Arctic Ocean codfishes Arctogadus glacialis and Boreogadus saida reveal oriL and tRNA gene duplications. Polar Biol. 2008, 31: 1245-1252.View ArticleGoogle Scholar
- Poulsen JY, Byrkjedal I, Willassen E, Rees D, Takeshima H, Satoh TP, Shinohara G, Nishida M, Miya M: Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei). BMC Evol Biol. 2013, 13: 111-PubMedPubMed CentralView ArticleGoogle Scholar
- Härlid A, Janke A, Arnason U: The mtDNA sequence of the ostrich and the divergence between Paleognathous and Neognathous birds. Mol Biol Evol. 1997, 14: 754-761.PubMedView ArticleGoogle Scholar
- Roe BA, Ma DP, Wilson RK, Wong JF: The complete nucleotide sequence of the Xenopus laevis mitochondrial genome. J Biol Chem. 1985, 260: 9759-9774.PubMedGoogle Scholar
- Lloyd RE, Foster PG, Guille M, Littlewood DTJ: Next generation sequencing and comparative analyses of Xenopus mitogenomes. BMC Genomics. 2012, 13: 496-PubMedPubMed CentralView ArticleGoogle Scholar
- Macay JR, Schulte JA, Larson A, Papenfuss TJ: Tandem duplications via light-strand synthesis may provide a precursor for mitochondrial genomic rearrangement. Mol Biol Evol. 1998, 15: 71-75.View ArticleGoogle Scholar
- McKnight ML, Shaffer HB: Large, rapidly evolving intergenic spacers in the mitochondrial DNA of the salamander family Ambystomatidae (Amphibia: Caudata). Mol Biol Evol. 1997, 14: 1167-1176.PubMedView ArticleGoogle Scholar
- Samuels AK, Weisrock DW, Smith JJ, France KJ, Walker JA, Putta S, Voss SR: Transcriptional and phylogenetic analysis of five complete ambystomatid salamander mitochondrial genomes. Gene. 2005, 349: 43-53.PubMedView ArticleGoogle Scholar
- Johansen S, Guddal PH, Johansen T: Organization of the mitochondrial genome of Atlantic cod, Gadus morhua. Nucleic Acids Res. 1990, 18: 411-419.PubMedPubMed CentralView ArticleGoogle Scholar
- Johansen S, Johansen T: Sequence analysis of 12 structural genes and a novel non-coding region from mitochondrial DNA of Atlantic cod, Gadus morhua. Biochem Biophys Acta. 1994, 1218: 213-217.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Shields GF, Gust JR: Lack of geographic structure in mitochondrial DNA sequences of Bering Sea walleye Pollock, Theragra chalcogramma. Mol Mar Biol Biotechnol. 1995, 4: 69-82.PubMedGoogle Scholar
- Bakke I, Shields GF, Johansen S: Sequence characterization of a unique intergenic spacer in gadiformes mitochondrial DNA. Mar Biotechnol. 1999, 1: 411-415.PubMedView ArticleGoogle Scholar
- Sigurgíslason H, Árnason E: Extent of mitochondrial DNA sequence variation in Atlantic cod from the Faroe Islands: a resolution of gene genealogy. Heredity. 2003, 91: 557-564.PubMedView ArticleGoogle Scholar
- Pálsson S, Paulsen J, Árnason E: Rapid evolution of the intergenic T-P spacer in the mtDNA of Arctic cod Arctogadus glacialis. Mar Biotechnol. 2008, 10: 270-277.PubMedView ArticleGoogle Scholar
- Mercer TR, Neph S, Dinger ME, Crawford J, Smith MA, Shearwood AMJ, Haugen E, Bracken CP, Rackham O, Stamatoyannopoulos JA, Filipovska A, Mattick JS: The human mitochondrial transcriptome. Cell. 2011, 146: 645-658.PubMedPubMed CentralView ArticleGoogle Scholar
- Rackham O, Shearwood AMJ, Mercer TR, Davies SMK, Mattick JS, Filipovska A: Long noncoding RNAs are generated from the mitochondrial genome and regulated by nuclear-encoded proteins. RNA. 2011, 17: 2085-2093.PubMedPubMed CentralView ArticleGoogle Scholar
- Temperley RJ, Wydro M, Lightowlers RN, Chrzanowska-Lightowlers ZM: Human mitochondrial mRNAs – like members of all families, similar but different. Biochim Biophys Acta. 2010, 1797: 1081-1085.PubMedPubMed CentralView ArticleGoogle Scholar
- Coucheron DH, Nymark M, Breines R, Karlsen BO, Andreassen M, Jørgensen TE, Moum T, Johansen SD: Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals. Curr Genet. 2011, 57: 213-222.PubMedPubMed CentralView ArticleGoogle Scholar
- Fatica A, Bozzoni I: Long non-coding RNAs: new players in cell differentiation and development. Nature Rev Genet. 2014, 15: 7-21.PubMedView ArticleGoogle Scholar
- Necsulea A, Soumillon M, Warnefors M, Liechti A, Daish T, Zeller U, Baker JC, Grutzner F, Kaessmann H: The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature. 2014, 505: 635-640.PubMedView ArticleGoogle Scholar
- Ojala D, Crews S, Montoya J, Gelfand R, Attardi G: A small polyadenylated RNA (7S RNA), containing a putative ribosome attachment site, maps near the origin of human mitochondrial DNA replication. J Mol Biol. 1981, 150: 303-314.PubMedView ArticleGoogle Scholar
- Chang DD, Clayton DA: Precise identification of individual promoter for transcription of each strand of human mitochondrial DNA. Cell. 1984, 36: 635-643.PubMedView ArticleGoogle Scholar
- Kumarswamy R, Bauters C, Volkmann I, Maury F, Fetisch J, Holzmann A, Lemesle G, de Groote P, Pinet F, Thum T: Circulating long noncoding RNA, LIPCAR, predicts survival in patients with heart failure. Circ Res. 2014, 114: 1569-1575.PubMedView ArticleGoogle Scholar
- Bai Y, Shakeley RM, Attardi G: Tight control of respiration by NADH dehydrogenase ND5 subunit gene expression in mouse mitochondria. Mol Cell Biol. 2000, 20: 805-815.PubMedPubMed CentralView ArticleGoogle Scholar
- Chomyn A: Mitochondrial genetic control of assembly and function of complex I in mammalian cells. J Bioenerg Biomembr. 2001, 33: 251-257.PubMedView ArticleGoogle Scholar
- Moum T, Willassen NP, Johansen S: Intragenic rearrangements in the mitochondrial NADH dehydrogenase subunit 6 gene of vertebrates. Curr Genet. 1994, 25: 554-557.PubMedView ArticleGoogle Scholar
- Bakke I, Johansen SD: Molecular phylogenetics of gadidae and related gadiformes based on mitochondrial DNA sequences. Mar Biotechnol. 2005, 7: 61-69.PubMedView ArticleGoogle Scholar
- Teletchea F, Laudet V, Hanni C: Phylogeny of the Gadidae (sensu Svetovidov, 1948) based on their morphology and two mitochondrial genes. Mol Phylogen Evol. 2006, 38: 189-199.View ArticleGoogle Scholar
- Roa-Varón A, Ortí G: Phylogenetic relationships among families of Gadiformes (Teleostei, Paracanthopterygii) based on nuclear and mitochondrial data. Mol Phylogenet Evol. 2009, 52: 688-704.PubMedView ArticleGoogle Scholar
- Johansen SD, Coucheron DH, Andreassen M, Karlsen BO, Furmanek T, Jørgensen TE, Emblem Å, Breines R, Nordeide JT, Moum T, Nederbragt AJ, Stenseth NC, Jakobsen KS: Large-scale sequence analyses of Atlantic cod. N Biotechnol. 2009, 25: 263-271.PubMedView ArticleGoogle Scholar
- 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-View ArticleGoogle Scholar
- Karlsen BO, Emblem A, Jørgensen TE, Klingan KA, Nordeide JT, Moum T, Johansen SD: Mitogenome sequence variation in migratory and stationary ecotypes of North-east Atlantic cod. Mar Genomics. 2014, 15: 103-108.PubMedView ArticleGoogle Scholar
- Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrøm M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, Wetten OF, Lanzén A, Winer R, Knight J, Vogel JH, Aken B, Andersen O, Lagesen K, Tooming-Klunderud A, Edvardsen RB, Tina KG, Espelund M, Nepal C, Previti C, Karlsen BO, Moum T, Skage M, Berg PR, Gjøen T, Kuhl H, et al: The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011, 477: 207-210.PubMedPubMed CentralView ArticleGoogle Scholar
- Rossmanith W: Of P and Z: Mitochondrial tRNA processing enzymes. Biochim Biophys Acta. 2012, 1810: 1017-1026.View ArticleGoogle Scholar
- Madsen CS, Ghivizzani SC, Hauswirth WW: Protein binding to a single termination-associated sequence in the mitochondrial DNA D-loop region. Mol Cell Biol. 1993, 13: 2162-2171.PubMedPubMed CentralView ArticleGoogle Scholar
- Roberti M, Musicco C, Polosa PL, Milella F, Gadaleta MN, Cantatore P: Multiple protein-binding sites in the TAS-region of human and rat mitochondrial DNA. Biochem Biophys Res Com. 1998, 243: 36-40.PubMedView ArticleGoogle Scholar
- Bakke I, Johansen S: Characterization of mitochondrial ribosomal RNA genes in gadiformes: sequence variations, secondary structural features, and phylogenetic implications. Mol Phylogen Evol. 2002, 25: 87-100.View ArticleGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2001, 302: 205-217.View ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739.PubMedPubMed CentralView ArticleGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinform. 2005, 21: 2104-2105.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.