Round and pointed-head grenadier fishes (Actinopterygii: Gadiformes) represent a single sister group: Evidence from the complete mitochondrial genome sequences

Complete Mitochondrial Genome of Four Peristediidae Fish Species: Genome Characterization and Phylogenetic Analysis

The systematic revision of the family Peristediidae remains an unresolved issue due to their diverse and unique morphology. Despite the popularity of using mitochondrial genome research to comprehensively understand phylogenetic relationships in fish, genetic data for peristediid fish need to be included. Therefore, this study aims to investigate the mitochondrial genomic characteristics and intra-family phylogenetic relationships of Peristediidae by utilizing mitochondrial genome analysis. Therefore, this study aims to investigate the phylogenetic relationship of Peristediidae by utilizing mitochondrial genome analysis. The mitochondrial genome of four species of Peristediidae (Peristedion liorhynchus, Satyrichthys welchi, Satyrichthys rieffeli, and Scalicus amiscus) collected in the East China Sea was studied. The mitochondrial gene sequence lengths of four fish species were 16,533 bp, 16,526 bp, 16,527 bp, and 16,526 bp, respectively. They had the same mitochondrial structure and were all composed of 37 genes and one control region. Most PCGs used ATG as the start codon, and a few used GTG as the start codon. An incomplete stop codon (TA/T) occurred. The AT-skew and GC-skew values of 13 PCGs from four species were negative, and the GC-skew amplitude was greater than that of AT-skew. All cases of D-arm were found in tRNA-Ser (GCT). The Ka/Ks ratio analysis indicated that 13 PCGs were suffering purifying selection. Based on 12 PCGs (excluding ND6) sequences, a phylogenetic tree was constructed using Bayesian inference (BI) and maximum likelihood (ML) methods, providing a further supplement to the scientific classification of Peristediidae fish. According to the results of divergence time, the four species of fish had apparent divergence in the Early Cenozoic, which indicates that the geological events at that time caused the climax of species divergence and evolution.

Comparative mitochondrial genome brings insights to slight variation in gene proportion and large intergenic spacer and phylogenetic relationship of mudskipper species

Fish mitochondrial genome have been largely studied worldwide for evolutionary and other genetic purposes and the structure and gene organization are commonly conservative. However, several studies have demonstrated that this scenario may present variations in some taxa, showing differentiation on the gene rearrangement. In this study, the complete mitogenome of terrestrial fish Boleophthalmus dussumieri was generated and compared with other species of the Exudercidae fishes. The newly complete mitogenome generated is circular and 16,685 bp of length, and it contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA genes (tRNAs), and one control region (CR), with high conservative structure, like other Mudskippers. Most of the PCG showed similar codon usage bias. The gene length was found to be different specially for the CR, 12S rRNA gene and ND5 gene in some taxon. All the Boleophthalmus species showed a gene duplication in the CR, except for B. dussumieri, and they presented a long intergenic spacer specially on the tRNA-Pro/ OH Tandem duplication/random loss (TDRL) and dimer-mitogenome and nonrandom loss (DMNL) are suitable to explain the mitogenome rearrangement observed in this study. The phylogenetic analysis well supported the monophyly of all mudskipper species and the analysis positioned the Periophthalmus clade as the most basal of the terrestrial fishes. This finding provides basis and brings insights for gene variation, gene rearrangements and replications showing evidence for variety of mitochondrial structure diversity within mudskippers.

Comparative analyses of the complete mitochondrial genomes of two southern African endemic guitarfish, Acroteriobatus annulatus and A. blochii

Shark-like rays (order Rhinopristiformes) are among the most threatened cartilaginous fish globally. Despite this, unresolved taxonomic issues still exist within the group. To date, no studies have used complete mitochondrial genomes to assess the phylogenetic placement of Acroteriobatus within the non-monophyletic family Rhinobatidae. The current study reports the first complete mitochondrial genomes for Acroteriobatus annulatus and A. blochii. Similar to other rhinopristiforms, the complete sequences of A. annulatus (16,773 bp) and A. blochii (16,771 bp) were circular molecules with gene organisations identical to that of the typical vertebrate mitogenome. The A + T content was higher than the G + C content, with a bias towards A and C nucleotides observed in all complete mitogenomes. The stem-and-loop secondary structures of the putative origin of light-strand replication were found to have highly conserved synthesis and stem regions, with all substitutions and indels restricted to the loop structure. The ratios of non-synonymous to synonymous substitution rates indicated that purifying selection has been the dominant driver of evolution in rhinopristiform mitogenomes. Phylogenetic reconstructions placed Acroteriobatus as a sister-group to Rhinobatos, confirming its affiliation with the family Rhinobatidae. However, based on its apparent polyphyly with the aforementioned genera, the familial assignment of Pseudobatos is not fully resolved and requires further investigation.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

Comparative Analysis of the Mitochondrial Genomes of Nicotiana tabacum: Hints Toward the Key Factors Closely Related to the Cytoplasmic Male Sterility Mechanism

Background

Cytoplasmic male sterility (CMS) is a complex phenomenon of plant sterility that can produce non-functional pollen. It is caused by mutation, rearrangement or recombination in the mitochondrial genome. So far, the systematic structural characteristics of the changes in the mitochondrial genome from the maintainer lines to the CMS lines have not been reported in tobacco.

Results

The mitochondrial genomes of the flower buds from both CMS lines and maintainer lines of two Nicotiana tabacum cultivars (YY85, sYY85, ZY90, and sZY90) were sequenced using the PacBio and Illumina Hiseq technology, and several findings were produced by comparative analysis based on the de novo sequencing. (1) The genomes of the CMS lines were larger, and the different areas were mostly non-coding regions. (2) A large number of rearrangement regions were detected in the CMS lines, with many translocation regions. (3) Thirteen gene clusters were shared by the four mitochondrial genomes, among which two of the gene clusters, nad2-sdh3 and nad6-rps4, were far from each other in the CMS lines. (4) Thirty-three protein-coding genes were conserved in four mitochondrial genomes. However, nad3 was detected one additional copy in the maintainer lines, and sequence differences were revealed in the four candidate genes (atp6, cox2, nad2, and sdh3). Importantly, the evolutionary tree based on the four genes could be used to distinguish the CMS lines and the maintainer lines well for the sequenced mitochondrial genomes of the tobacco. (5) Sixteen CMS-specific open reading frames (ORFs) were found, three of which (orf91, orf115b, and orf100) were previously reported. (6) The differences in intensity of the protein–protein (PPI) interaction in ATP6 were further verified using the yeast two-hybrid analysis.

Conclusion

Although the majority of the sequences, genes and gene clusters were shared by the mitochondrial genomes of the maintainer and the CMS lines in tobacco, extensive structural variations identified with comprehensive analysis based on the mitochondrial genomes, including rearrangement, gene order, the mitochondrial genome expansion and shrinkage events, might be related to CMS. Additionally, the candidate protein-coding genes and CMS-specific ORFs were closely associated with the CMS mechanism. Verification experiments of one of the candidate genes were performed, and the validity of our research results was supported.

Mechanisms of gene rearrangement in 13 bothids based on comparison with a newly completed mitogenome of the threespot flounder, Grammatobothus polyophthalmus (Pleuronectiformes: Bothidae)

Background

The mitogenomes of 12 teleost fish of the bothid family (order Pleuronectiformes) indicated that the genomic-scale rearrangements characterized in previous work. A novel mechanism of genomic rearrangement called the Dimer-Mitogenome and Non-Random Loss (DMNL) model was used to account for the rearrangement found in one of these bothids, Crossorhombus azureus.

Results

The 18,170 bp mitogenome of G. polyophthalmus contains 37 genes, two control regions (CRs), and the origin of replication of the L-strand (O_L). This mitogenome is characterized by genomic-scale rearrangements: genes located on the L-strand are grouped in an 8-gene cluster (Q-A-C-Y-S₁-ND6-E-P) that does not include tRNA-N; genes found on the H-strand are grouped together (F-12S … CytB-T) except for tRNA-D that was translocated inside the 8-gene L-strand cluster. Compared to non-rearranged mitogenomes of teleost fishes, gene organization in the mitogenome of G. polyophthalmus and in that of the other 12 bothids characterized thus far is very similar. These rearrangements could be sorted into four types (Type I, II, III and IV), differing in the particular combination of the CR, tRNA-D gene and 8-gene cluster and the shuffling of tRNA-V. The DMNL model was used to account for all but one gene rearrangement found in all 13 bothid mitogenomes. Translocation of tRNA-D most likely occurred after the DMNL process in 10 bothid mitogenomes and could have occurred either before or after DMNL in the three other species. During the DMNL process, the tRNA-N gene was retained rather than the expected tRNA-N′ gene. tRNA-N appears to assist in or act as O_L function when the O_L secondary structure could not be formed from intergenic sequences. A striking finding was that each of the non-transcribed genes has degenerated to a shorter intergenic spacer during the DMNL process. These findings highlight a rare phenomenon in teleost fish.

Conclusions

This result provides significant evidence to support the existence of dynamic dimeric mitogenomes and the DMNL model as the mechanism of gene rearrangement in bothid mitogenomes, which not only promotes the understanding of mitogenome structural diversity, but also sheds light on mechanisms of mitochondrial genome rearrangement and replication.

Resolving deep-sea pelagic saccopharyngiform eel mysteries: Identification of Neocyema and Monognathidae leptocephali and establishment of a new fish family "Neocyematidae" based on larvae, adults and mitogenomic gene orders

Deep-sea midwater "saccopharyngiform" eels of the families Cyematidae, Monognathidae, Eurypharyngidae and Saccopharyngidae (order Anguilliformes) are extraordinary fishes having major skeletal reductions and modifications compared to the general anguilliform body structure. Little is known about most aspects of the systematics, phylogeny, and ecology of these families, and few of the approximately 30 species described from adult specimens have been matched with their leptotocephalus larvae. Based on mitogenomic sequence data from rare new specimens, we show that the long-speculated-about larval form referred to as "Leptocephalus holti", which was thought to possibly be the larva of the rare orange-colored eels of Neocyema (5 known specimens; speculated to belong to the Cyematidae) are actually the larvae of the one-jaw eels of the family Monognathidae. One of the 5 types of L. holti larvae that were collected in the Pacific is genetically matched with Monognathus jesperseni, but multiple species exist based on larval sequence data and the morphology of adult specimens. A rare leptocephalus from the Sargasso Sea, with unique morphological characteristics including many small orange spots on the gut, was found to be the larva of Neocyema, which is presently only known from the Atlantic Ocean. We demonstrate that Neocyema constitutes a separate family being most closely related to Eurypharyngidae and Saccopharyngidae based on mitogenomic DNA sequences and unique mitochondrial gene orders.

Characterization of the Complete Mitochondrial Genome Sequences of Three Croakers (Perciformes, Sciaenidae) and Novel Insights into the Phylogenetics

The three croakers (Nibea coibor, Protonibea diacanthus and Argyrosomus amoyensis, Perciformes, Sciaenidae) are important commercial species inhabiting the Eastern Indian Ocean and Western Pacific. Molecular data employed in previous research on phylogenetic reconstruction have not been adequate and complete, and systematic and comprehensive phylogenetic relationships for these fish are unresolved. We sequenced the complete mitochondrial genomes of the three croakers using next-generation sequencing for the first time. We analyzed the composition and phylogenies between 19 species in the family Sciaenidae using the mitochondrial protein coding sequences of 204 species in the Series Eupercaria. We present the characterization of the complete mitochondrial genome sequences of the three croakers. Gene arrangement and distribution of the three croakers are canonically identical and consistent with other vertebrates. We found that the family Sciaenidae is an independent branch that is isolated from the order Perciformes and does not belong to any extant classification. Therefore, this family is expected to belong to a new classification at the order level and needs further analysis. The evolution of Sciaenidae has lagged far behind the Perciformes differentiation. This study presents a novel insight into the phylogenetics of the family Sciaenidae from the order Perciformes and facilitates additional studies on the evolution and phylogeny of Series Eupercaria.

Characterization of the complete mitogenome of Indian Mouse Deer, Moschiola indica (Artiodactyla: Tragulidae) and its evolutionary significance

The mitochondrial genome of Indian mouse deer (Moschiola indica) was sequenced, assembled and characterized for the first time using 22 pairs of polymerase chain reaction (PCR) primers. The mitogenome of M. indica which is 16,444 bp in size was found very similar to most vertebrates in organisation that harbours 13 protein-coding genes, 22 transfer RNA, 2 ribosomal RNA and 1A + T-rich region. Its comparison with over 52 mitogenomes of the order Artiodactyla, showed a conserved nature of gene organisation, codon usage, gene orientation and evolutionary rates of proteins except that M. indica possesses an extra copy of trnF. The complete mitogenome and protein-coding genes of M. indica were found to be highly A + T biased. Rate of protein evolution was highest in atp8 and lowest in cox3. Further, a higher purifying selection pressure was found to be acting on family Tragulidae compared to Bovidae and Cervidae. The phylogenetic analysis of M. indica placed the Tragulidae as sister-group of all other ruminants, similar to previous analyses. Moschiola forms the sister-group to the other two tragulid genera Tragulus (from Asia) and Hyemoschus (from Africa), which is unexpected as usually the Asian species are thought to form a monophyletic group.

Sequencing and characterization of the complete mitochondrial genome of Japanese Swellshark (Cephalloscyllium umbratile)

To further comprehend the genome features of Cephalloscyllium umbratile (Carcharhiniformes), an endangered species, the complete mitochondrial DNA (mtDNA) was firstly sequenced and annotated. The full-length mtDNA of C. umbratile was 16,697 bp and contained ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, and a major non-coding control region. Each PCG was initiated by an authoritative ATN codon, except for COX1 initiated by a GTG codon. Seven of 13 PCGs had a typical TAA termination codon, while others terminated with a single T or TA. Moreover, the relative synonymous codon usage of the 13 PCGs was consistent with that of other published Carcharhiniformes. All tRNA genes had typical clover-leaf secondary structures, except for tRNA-Ser (GCT), which lacked the dihydrouridine 'DHU' arm. Furthermore, the analysis of the average Ka/Ks in the 13 PCGs of three Carcharhiniformes species indicated a strong purifying selection within this group. In addition, phylogenetic analysis revealed that C. umbratile was closely related to Glyphis glyphis and Glyphis garricki. Our data supply a useful resource for further studies on genetic diversity and population structure of C. umbratile.

Structure and variation of the mitochondrial genome of fishes

Background

The mitochondrial (mt) genome has been used as an effective tool for phylogenetic and population genetic analyses in vertebrates. However, the structure and variability of the vertebrate mt genome are not well understood. A potential strategy for improving our understanding is to conduct a comprehensive comparative study of large mt genome data. The aim of this study was to characterize the structure and variability of the fish mt genome through comparative analysis of large datasets.

Results

An analysis of the secondary structure of proteins for 250 fish species (248 ray-finned and 2 cartilaginous fishes) illustrated that cytochrome c oxidase subunits (COI, COII, and COIII) and a cytochrome bc1 complex subunit (Cyt b) had substantial amino acid conservation. Among the four proteins, COI was the most conserved, as more than half of all amino acid sites were invariable among the 250 species. Our models identified 43 and 58 stems within 12S rRNA and 16S rRNA, respectively, with larger numbers than proposed previously for vertebrates. The models also identified 149 and 319 invariable sites in 12S rRNA and 16S rRNA, respectively, in all fishes. In particular, the present result verified that a region corresponding to the peptidyl transferase center in prokaryotic 23S rRNA, which is homologous to mt 16S rRNA, is also conserved in fish mt 16S rRNA. Concerning the gene order, we found 35 variations (in 32 families) that deviated from the common gene order in vertebrates. These gene rearrangements were mostly observed in the area spanning the ND5 gene to the control region as well as two tRNA gene cluster regions (IQM and WANCY regions). Although many of such gene rearrangements were unique to a specific taxon, some were shared polyphyletically between distantly related species.

Conclusions

Through a large-scale comparative analysis of 250 fish species mt genomes, we elucidated various structural aspects of the fish mt genome and the encoded genes. The present results will be important for understanding functions of the mt genome and developing programs for nucleotide sequence analysis. This study demonstrated the significance of extensive comparisons for understanding the structure of the mt genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3054-y) contains supplementary material, which is available to authorized users.

Characterization of the Complete Mitochondrial Genome Sequence of the Globose Head Whiptail Cetonurus globiceps (Gadiformes: Macrouridae) and Its Phylogenetic Analysis

The particular environmental characteristics of deep water such as its immense scale and high pressure systems, presents technological problems that have prevented research to broaden our knowledge of deep-sea fish. Here, we described the mitogenome sequence of a deep-sea fish, Cetonurus globiceps. The genome is 17,137 bp in length, with a standard set of 22 transfer RNA genes (tRNAs), two ribosomal RNA genes, 13 protein-coding genes, and two typical non-coding control regions. Additionally, a 70 bp tRNA(Thr)-tRNA(Pro) intergenic spacer is present. The C. globiceps mitogenome exhibited strand-specific asymmetry in nucleotide composition. The AT-skew and GC-skew values in the whole genome of C. globiceps were 0 and -0.2877, respectively, revealing that the H-strand had equal amounts of A and T and that the overall nucleotide composition was C skewed. All of the tRNA genes could be folded into cloverleaf secondary structures, while the secondary structure of tRNA(Ser(AGY)) lacked a discernible dihydrouridine stem. By comparing this genome sequence with the recognition sites in teleost species, several conserved sequence blocks were identified in the control region. However, the GTGGG-box, the typical characteristic of conserved sequence block E (CSB-E), was absent. Notably, tandem repeats were identified in the 3' portion of the control region. No similar repetitive motifs are present in most of other gadiform species. Phylogenetic analysis based on 12 protein coding genes provided strong support that C. globiceps was the most derived in the clade. Some relationships however, are in contrast with those presented in previous studies. This study enriches our knowledge of mitogenomes of the genus Cetonurus and provides valuable information on the evolution of Macrouridae mtDNA and deep-sea fish.

A new classification of viviparous brotulas (Bythitidae) - with family status for Dinematichthyidae - based on molecular, morphological and fossil data

The order Ophidiiformes is a large but not very well known group of fishes, unique among teleosts for showing high diversity in both deep sea and shallow reef habitats. The current classification includes more than 500 species, 115 genera and four families, based primarily on mode of reproduction: viviparous Aphyonidae and Bythitidae vs oviparous Carapidae and Ophidiidae. Since 2004 we revised the bythitid tribe Dinematichthyini, described more than 100 new species and noticed that this group has unique morphological characters, perhaps supporting a higher level of classification than the current status. Here we study the viviparous families phylogenetically with partial mitochondrial (nd4, 16s) and nuclear (Rag1) DNA sequences (2194bp). We use a fossil calibration of otolith-based taxa to calibrate the age of the clade comprising bythitid and dinematicththyid representatives, together with fossil calibrations adopted from previous phylogenetic studies. The separation of the order into two major lineages, the viviparous Bythitoidei and the oviparous Ophidioidei is confirmed. At the familial level, however, a new classification is presented for the viviparous clades, placing Aphyonidae as a derived, pedomorphic member of Bythitidae (new diagnosis provided, 33 genera and 118 species). The current subfamily Brosmophycinae is considered polyphyletic and we propose family status for Dinematichthyidae (25 genera, 114 species), supported by unique, morphological synapomorphic characters in the male copulatory apparatus. Previous use of the caudal fin separation or fusion with vertical fins is ambiguous. Age estimates based on calibrated molecular phylogeny agrees with fossil data, giving an origin within the Cretaceous (between 84 and 104mya) for a common ancestor to Ophidiiformes.

The inversion of the Control Region in three mitogenomes provides further evidence for an asymmetric model of vertebrate mtDNA replication

Mitochondrial genomes are known to have a strong strand-specific compositional bias that is more pronounced at fourfold redundant sites of mtDNA protein-coding genes. This observation suggests that strand asymmetries, to a large extent, are caused by mutational asymmetric mechanisms. In vertebrate mitogenomes, replication and not transcription seems to play a major role in shaping compositional bias. Hence, one can better understand how mtDNA is replicated – a debated issue – through a detailed picture of mitochondrial genome evolution. Here, we analyzed the compositional bias (AT and GC skews) in protein-coding genes of almost 2,500 complete vertebrate mitogenomes. We were able to identify three fish mitogenomes with inverted AT/GC skew coupled with an inversion of the Control Region. These findings suggest that the vertebrate mitochondrial replication mechanism is asymmetric and may invert its polarity, with the leading-strand becoming the lagging-strand and vice-versa, without compromising mtDNA maintenance and expression. The inversion of the strand-specific compositional bias through the inversion of the Control Region is in agreement with the strand-displacement model but it is also compatible with the RITOLS model of mtDNA replication.

An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA

Background

Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNA^Thr and tRNA^Pro 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.

Results

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).

Conclusions

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.

Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei)

Background

A skewed assemblage of two epi-, meso- and bathypelagic fish families makes up the order Myctophiformes – the blackchins Neoscopelidae and the lanternfishes Myctophidae. The six rare neoscopelids show few morphological specializations whereas the divergent myctophids have evolved into about 250 species, of which many show massive abundances and wide distributions. In fact, Myctophidae is by far the most abundant fish family in the world, with plausible estimates of more than half of the oceans combined fish biomass. Myctophids possess a unique communication system of species-specific photophore patterns and traditional intrafamilial classification has been established to reflect arrangements of photophores. Myctophids present the most diverse array of larval body forms found in fishes although this attribute has both corroborated and confounded phylogenetic hypotheses based on adult morphology. No molecular phylogeny is available for Myctophiformes, despite their importance within all ocean trophic cycles, open-ocean speciation and as an important part of neoteleost divergence. This study attempts to resolve major myctophiform phylogenies from both mitogenomic sequences and corroborating evidence in the form of unique mitochondrial gene order rearrangements.

Results

Mitogenomic evidence from DNA sequences and unique gene orders are highly congruent concerning phylogenetic resolution on several myctophiform classification levels, corroborating evidence from osteology, larval ontogeny and photophore patterns, although the lack of larval morphological characters within the subfamily Lampanyctinae stands out. Neoscopelidae is resolved as the sister family to myctophids with Solivomer arenidens positioned as a sister taxon to the remaining neoscopelids. The enigmatic Notolychnus valdiviae is placed as a sister taxon to all other myctophids and exhibits an unusual second copy of the tRNA-Met gene – a gene order rearrangement reminiscent of that found in the tribe Diaphini although our analyses show it to be independently derived. Most tribes are resolved in accordance with adult morphology although Gonichthyini is found within a subclade of the tribe Myctophini consisting of ctenoid scaled species. Mitogenomic sequence data from this study recognize 10 reciprocally monophyletic lineages within Myctophidae, with five of these clades delimited from additional rearranged gene orders or intergenic non-coding sequences.

Conclusions

Mitogenomic results from DNA sequences and unique gene orders corroborate morphology in phylogeny reconstruction and provide a likely scenario for the phylogenetic history of Myctophiformes. The extent of gene order rearrangements found within the mitochondrial genomes of myctophids is unique for phylogenetic purposes.

The complete mitochondrial genomes of largemouth bass of the northern subspecies (Micropterus salmoides salmoides) and Florida subspecies (Micropterus salmoides floridanus) and their applications in the identification of largemouth bass species

The largemouth bass belongs to the family Centrarchidae, which includes two subspecies: the northern subspecies, Micropterus salmoides salmoides, and the Florida subspecies, Micropterus salmoides floridanus. In this study, the complete mitochondrial genomes of the two subspecies were sequenced, and their genetic differences were identified. The mitogenomes of M. s. salmoides and M. s. floridanus are 16,486 and 16,479 bp in length, respectively. The two subspecies consisted of 37 genes (13 protein-coding genes, 2 ribosomal RNA, and 22 transfer RNA), which are typical for vertebrate mtDNA. Phylogenetic analysis provided statistical support for the monophyly of the family Centrarchidae. Comparison of the two subspecies' mitogenomes revealed a relatively high number (450) of single nucleotide polymorphisms (SNPs) in protein-coding genes. We characterized SNPs in the partial cytochrome c oxidase subunit 1 gene of different individuals from three cultured populations, one wild northern subspecies population, and one wild Florida subspecies population. Twenty-eight SNPs were fixed with alternative nucleotides in the two subspecies, which could be used for differentiating them. Based on this gene, phylogenetic tree and genetic distance analyses supported that cultured largemouth bass in China belongs to the northern subspecies.

Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals

Expression and processing of mitochondrial gene transcripts are fundamental to mitochondrial function, but information from early vertebrates like teleost fishes is essentially lacking. We have analyzed mitogenome sequences of ten codfishes (family Gadidae), and provide complete sequences from three new species (Saithe, Pollack and Blue whiting). Characterization of the mitochondrial mRNAs in Saithe and Atlantic cod identified a set of ten poly(A) transcripts, and six UAA stop codons are generated by posttranscriptional polyadenylation. Structural assessment of poly(A) sites is consistent with an RNaseP cleavage activity 5' of tRNA acceptor-like stems. COI, ND5 and ND6 mRNAs were found to harbor 3' UTRs with antisense potential extending into neighboring gene regions. While the 3' UTR of COI mRNA is complementary to the tRNA(Ser UCN) and highly similar to that detected in human mitochondria, the ND5 and ND6 3' UTRs appear more heterogenic. Deep sequencing confirms expression of all mitochondrial mRNAs and rRNAs, and provides information about the precise 5' ends in mature transcripts. Our study supports an overall evolutionary conservation in mitochondrial RNA processing events among vertebrates, but reveals some unique 5' and 3' end characteristics in codfish mRNAs with implications to antisense regulation of gene expression.

The complete mitochondrial genome of the relict frog Leiopelma archeyi: insights into the root of the frog Tree of Life

Determining the root of the anuran Tree of Life is still a contentious and open question in frog systematics. Two genera with disjunct distributions have been traditionally considered the most basal among extant frogs: Leiopelma, which is endemic to New Zealand, and Ascaphus, which lives in North America. However, their specific phylogenetic position is rather elusive because each genus shows many autapomorphies, and together they retain many symplesiomorphic characters. Therefore, several alternative hypotheses have been proposed regarding the relative phylogenetic position of both Leiopelma and Ascaphus. In order to distinguish among these competing phylogenetic hypotheses, we sequenced the complete mitochondrial (mt) genome of Leiopelma archeyi and used it along with previously reported frog mt genomes (including that of Ascaphus truei) to infer a robust phylogeny of major anuran lineages. The reconstructed maximum likelihood and Bayesian inference phylogenies recovered identical topology, which supports the sister group relationship of Ascaphus and Leiopelma, and the placement of this clade at the base of the anuran tree. Interestingly, the mt genome of L. archeyi displays a novel gene arrangement in frog mt genomes affecting the relative position of cytochrome b, trnT, NADH dehydrogenase subunit 6, trnE, and trnP genes. The tandem duplication-random loss model of gene order change explains the origin of this novel frog mt genome arrangement, which is convergent with others reported in some fishes and salamanders. These results, together with comparative data for other available vertebrate mt genomes, provide evidence that the 5' end of the control region is a hot spot for gene order rearrangement.

Evolutionary history of anglerfishes (Teleostei: Lophiiformes): a mitogenomic perspective

BACKGROUND

The teleost order Lophiiformes, commonly known as the anglerfishes, contains a diverse array of marine fishes, ranging from benthic shallow-water dwellers to highly modified deep-sea midwater species. They comprise 321 living species placed in 68 genera, 18 families and 5 suborders, but approximately half of the species diversity is occupied by deep-sea ceratioids distributed among 11 families. The evolutionary origins of such remarkable habitat and species diversity, however, remain elusive because of the lack of fresh material for a majority of the deep-sea ceratioids and incompleteness of the fossil record across all of the Lophiiformes. To obtain a comprehensive picture of the phylogeny and evolutionary history of the anglerfishes, we assembled whole mitochondrial genome (mitogenome) sequences from 39 lophiiforms (33 newly determined during this study) representing all five suborders and 17 of the 18 families. Sequences of 77 higher teleosts including the 39 lophiiform sequences were unambiguously aligned and subjected to phylogenetic analysis and divergence time estimation.

RESULTS

Partitioned maximum likelihood analysis confidently recovered monophyly for all of the higher taxa (including the order itself) with the exception of the Thaumatichthyidae (Lasiognathus was deeply nested within the Oneirodidae). The mitogenomic trees strongly support the most basal and an apical position of the Lophioidei and a clade comprising Chaunacoidei + Ceratioidei, respectively, although alternative phylogenetic positions of the remaining two suborders (Antennarioidei and Ogcocephaloidei) with respect to the above two lineages are statistically indistinguishable. While morphology-based intra-subordinal relationships for relatively shallow, benthic dwellers (Lophioidei, Antennarioidei, Ogcocephaloidei, Chaunacoidei) are either congruent with or statistically indistinguishable from the present mitogenomic tree, those of the principally deep-sea midwater dwellers (Ceratioidei) cannot be reconciled with the molecular phylogeny. A relaxed molecular-clock Bayesian analysis of the divergence times suggests that all of the subordinal diversifications have occurred during a relatively short time period between 100 and 130 Myr ago (early to mid Cretaceous).

CONCLUSIONS

The mitogenomic analyses revealed previously unappreciated phylogenetic relationships among the lophiiform suborders and ceratioid familes. Although the latter relationships cannot be reconciled with the earlier hypotheses based on morphology, we found that simple exclusion of the reductive or simplified characters can alleviate some of the conflict. The acquisition of novel features, such as male dwarfism, bioluminescent lures, and unique reproductive modes allowed the deep-sea ceratioids to diversify rapidly in a largely unexploited, food-poor bathypelagic zone (200-2000 m depth) relative to the other lophiiforms occurring in shallow coastal areas.

A novel rearrangement in the mitochondrial genome of tongue sole, Cynoglossus semilaevis: control region translocation and a tRNA gene inversion

The organization of fish mitochondrial genomes (mitogenomes) is quite conserved, usually with the heavy strand encoding 12 of 13 protein-coding genes and 14 of 22 tRNA genes, and the light strand encoding ND6 and the remaining 8 tRNA genes. Currently, there are only a few reports on gene reorganization of fish mitogenomes, with only two types of rearrangements (shuffling and translocation) observed. No gene inversion has been detected in approximately 420 complete fish mitogenomes available so far. Here we report a novel rearrangement in the mitogenome of Cynoglossus semilaevis (Cynoglossinae, Cynoglossidae, Pleuronectiformes). The genome is 16 371 bp in length and contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 2 main noncoding regions, the putative control region and the light-strand replication origin. A striking finding of this study is that the tRNAGln gene is translocated from the light to the heavy strand (Q inversion). This is accompanied by shuffling of the tRNAIle gene and long-range translocation of the putative control region downstream to a site between ND1 and the tRNAGln gene. The remaining gene order is identical to that of typical fish mitogenomes. Additionally, unique characters of this mitogenome, including a high A+T content and length variations of 8 protein-coding genes, were found through comparison of the mitogenome sequence with those from other flatfishes. All the features detected and their relationships with the rearrangements, as well as a possible rearrangement pathway, are discussed. These data provide interesting information for better understanding the molecular mechanisms of gene reorganization in fish mitogenomes.

Phylogenetic relationships among families of Gadiformes (Teleostei, Paracanthopterygii) based on nuclear and mitochondrial data

Phylogenetic hypotheses among Gadiformes fishes at the suborder, family, and subfamily levels are controversial. To address this problem, we analyze nuclear and mitochondrial DNA (mtDNA) sequences for the most extensive taxonomic sampling compiled to date, representing all of the recognized families and subfamilies in the order (except the monotypic family Lyconidae). Our study sampled 117 species from 46 genera, comprising around 20% of the species described for the order (more than 60% of all genera in the order) and produced 2740 bp of DNA sequence data for each species. Our analysis was successful in confirming the monophyly of Gadiformes and most of the proposed families for the order, but alternative hypotheses of sister-group relationships among families were poorly resolved. Our results are consistent with dividing Gadiformes into 12 families in three suborders, Muraenolepidoidei, Macrouroidei, and Gadoidei. Muraenolepidoidei contains the single family Muraenolepididae. The suborder Macrouroidei includes at least three families: Macrouridae, Macruronidae and Steindachneriidae. Macrouridae is deeply divided into two well-supported subfamilies: Macrourinae and Bathygadinae, suggesting that Bathygadinae may be ranked at the family level. The suborder Gadoidei includes the families: Merlucciidae, Melanonidae, Euclichthyidae, Gadidae, Ranicipitidae, and Bregmacerotidae. Additionally, Trachyrincinae could be ranked at family level including two subfamilies: Trachyrincinae and Macrouroidinae within Gadoidei. Further taxonomic sampling and sequencing efforts are needed in order to corroborate these relationships.

The complete mitochondrial genome of rock carp Procypris rabaudi (Cypriniformes: Cyprinidae) and phylogenetic implications

Rock carp, Procypris rabaudi (Tchang), is an endemic fish species in China. We sequenced the complete mitochondrial genome of it by high-fidelity polymerase chain reaction with conserved primers and primer walking sequencing method. The complete mitochondrial genome of rock carp is 16595 bp in length and contains 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region, with an identical order to that of most other vertebrates. The origin of L-strand replication (OL) in rock carp mitochondrion is located in a cluster of five tRNA genes (WANCY region) with 35 nucleotides in length. The control region is located between the tRNA-Pro and tRNA-Phe genes and is 943 bp in length. Three conserved sequence blocks (CSB), an extended termination associated sequence (ETAS), an AT-repeat microsatellite sequence and a putative promoter sequence for H-strand transcription (HSP) were identified within this region. The microsatellite sequence has a very low variation, with only one repeat alteration in 50 checked individuals (from 12 to 13 repeats). The phylogenetic analysis for rock carp was performed with Bayesian and Maximum likelihood (ML) methods based on the concatenated nucleotide sequence of 12 protein-coding genes on the heavy strand. The result suggested that traditional taxonomic barbines possibly originated more early than cyprininaes; rock carp was placed at the position between barbines and cyprininaes, while has a closer relationship with cyprininaes than barbines.