Complete mitochondrial genome of the blackspot seabream, Pagellus bogaraveo (Perciformes: Sparidae), with high levels of length heteroplasmy in the WANCY region

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

Mitochondrial Analysis of Sparidae Species to Detect a New DNA Barcoding Marker for Dentex gibbosus to Utilize against Fraud

Dentex gibbosus (Pink dentex) is a fish species of increasing economic interest in the Mediterranean Sea that is consumed both whole and processed. The growing value of this sparid in European markets is responsible for its substitution with fraudulent species. The distinctive morphologic feature of D. gibbosus is the conspicuous hump on the forehead in the older and larger specimens. However, the head is regularly convex in young individuals, requiring high skills and competencies for correct identification. Authentication becomes even more challenging in the case of prepared and processed products. Therefore, the molecular characterization of Pink dentex plays a crucial role in preventing commercial fraud with species substitution. This paper proposes a comparative mitogenome analysis between 19 sparid species of commercial interest as a tool to accurately design species-specific primers targeting a fragment of the NAD2 gene for the identification of D. gibbosus. We successfully detected Pink dentex DNA both using endpoint and real-time PCR. The findings showed the high specificity of the designed primers, demonstrating this a suitable, fast, and cost-effective method that could be used for the unambiguous identification of Pink dentex. This innovative approach for sparid authentication is expected to contribute to seafood traceability, public health assurance, integrity, and the credibility of the seafood industry.

A Rapid Method for the Identification of Fresh and Processed Pagellus erythrinus Species against Frauds

The commercialization of porgies or seabreams of the family Sparidae has greatly increased in the last decade, and some valuable species have become subject to seafood substitution. DNA regions currently used for fish species identification in fresh and processed products belong to the mitochondrial (mt) genes cytochrome b (Cytb), cytochrome c oxidase I (COI), 16S and 12S. However, these markers amplify for fragments with lower divergence within and between some species, failing to provide informative barcodes. We adopted comparative mitogenomics, through the analysis of complete mtDNA sequences, as a compatible approach toward studying new barcoding markers. The intent is to develop a specific and rapid assay for the identification of the common pandora Pagellus erythrinus, a sparid species frequently subject to fraudulent replacement. The genetic diversity analysis (Hamming distance, p-genetic distance, gene-by-gene sequence variability) between 16 sparid mtDNA genomes highlighted the discriminating potential of a 291 bp NAD2 gene fragment. A pair of species-specific primers were successfully designed and tested by end-point and real-time PCR, achieving amplification only in P. erythrinus among several fish species. The use of the NAD2 barcoding marker provides a rapid presence/absence method for the identification of P. erythrinus.

One Hundred Mitochondrial Genomes of Cicadas

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

Rapid divergence, molecular evolution, and morphological diversification of coastal host-parasite systems from southern Brazil

This study assessed the role of historical processes on the geographic isolation, molecular evolution, and morphological diversification of host-parasite populations from the southern Brazilian coast. Adult specimens of Scleromystax barbatus and Scleromystax macropterus were collected from the sub-basin of the Nhundiaquara River and the sub-basin of the Paranaguá Bay, state of Paraná, Brazil. Four species of Gyrodactylus were recovered from the body surface of both host species. Morphometric analysis of Gyrodactylus spp. and Scleromystax spp. indicated that subpopulations of parasites and hosts could be distinguished from different sub-basins and locations, but the degree of morphological differentiation seems to be little related to geographic distance between subpopulations. Phylogenetic relationships based on DNA sequences of Gyrodactylus spp. and Scleromystax spp. allowed distinguishing lineages of parasites and hosts from different sub-basins. However, the level of genetic structuring of parasites was higher in comparison to host species. Evidence of positive selection in mtDNA sequences is likely associated with local adaptation of lineages of parasites and hosts. A historical demographic analysis revealed that populations of Gyrodactylus and Scleromystax have expanded in the last 250 000 years. The genetic variation of parasites and hosts is consistent with population-specific selection, population expansions, and recent evolutionary co-divergence.

Complete mitochondrial genome of Ophichthus brevicaudatus reveals novel gene order and phylogenetic relationships of Anguilliformes

Generally, a teleostean group possesses only one type or a set of similar mitochondrial gene arrangement. However, two types of gene arrangement have been identified in the mitochondrial genomes (mitogenomes) of Anguilliformes. Here, a newly sequenced mitogenome of Ophichthus brevicaudatus (Anguilliformes; Ophichthidae) was presented. The total length of the O. brevicaudatus mitogenome was 17,773 bp, and it contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and two identical control regions (CRs). The gene order differed from that of the typical vertebrate mitogenomes. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The duplication-random loss model was adopted to explain the gene rearrangement events in this mitogenome. The most comprehensive phylogenetic trees of Anguilliformes based on complete mitogenome was constructed. The non-monophyly of Congridae was well supported, whereas the non-monophyly of Derichthyidae and Chlopsidae was not supported. These results provide insight into gene arrangement features of anguilliform mitogenomes and lay the foundation for further phylogenetic studies on Anguilliformes.

Comparison of mitochondrial DNA enrichment and sequencing methods from fish tissue

Sparid fish species have different commercial values related to their organoleptic features. Mitochondrial (mt) DNA provides a potential tool to distinguish species, but the enrichment of high-quality mtDNA from total genomic DNA is critical to obtain entire mtDNA sequences. Conventional mtDNA isolation is relatively low-cost and proficient. However, high numbers of PCR cycles can lead to artefacts (10^-6 mutations/bp). We describe a rapid protocol for mtDNA extraction and enrichment from fish tissues, based on conventional miniprep columns and paramagnetic bead-based purification, without the need to employ PCR amplification. This newly described method generates a substrate for next-generation sequencing (NGS) analysis and is likely to have wider applications for mitochondrial studies in other fish families to help ensure traceability and differentiation of fish with high commercial values.

Phylogenetic analysis of sooty grunter and other major freshwater fishes in the suborder Percoidei based on mitochondrial DNA

Perciformes is the largest order of fishes and vertebrates. Sooty grunter (Hephaestus fuliginosus) is an economic fish species in the Terapontidae family of Percoidei, a suborder within Perciformes. To conduct molecular-level analysis of the phylogenetic relationships between sooty grunter and major freshwater fishes in Percoidei, we analysed the entire sooty grunter mitochondrial genome sequence and obtained the mitochondrial genome information of 19 fishes from Terapontidae, Serranidae, and Centrarchidae families in Percoidei from GenBank. The complete length of the sooty grunter mitochondrial genome was 16,770 bp; it encoded 13 proteins, 2 rRNAs, 22 tRNAs, and a displacement loop (D-loop). Other than ND6 and eight tRNA genes that are encoded by the light strand, the majority of genes are encoded by the heavy strand. The sequence and distribution of sooty grunter mitochondrial-encoded genes and non-coding segment were similar to those of most vertebrates. The results of neighbour joining, maximum parsimony, and Bayesian inference analyses of the complete mitochondrial genome and six genes, including cytochrome oxidase I, cytochrome B, 12S rRNA, ND2, ND4, and ND5, were consistent. In the phylogenetic trees, fishes in Terapontidae and Centrarchidae formed monophyletic clades, whereas those in Serranidae were divided into two clades, each containing Lateolabrax and Siniperca species. Among the three freshwater fish species in Terapontidae, the freshwater Terapontidae were more closely related to jade perch than with silver perch, suggesting that freshwater Terapontidae fishes originate from marine fishes. In addition, the phylogenetic results indicated that Micropterus salmoides salmoides and Micropterus salmoides floridanus in Centrarchidae should be designated as two independent species, and Siniperca in Serranidae should be considered an independent family. The sooty grunter mitochondrial genome sequence obtained in this study could be used to conduct population genetic diversity and germplasm resource studies. Furthermore, the phylogenetic analysis results of freshwater fishes in Percoidei could provide a molecular basis for cross-breeding.

Multiple independent structural dynamic events in the evolution of snake mitochondrial genomes

BACKGROUND

Mitochondrial DNA sequences have long been used in phylogenetic studies. However, little attention has been paid to the changes in gene arrangement patterns in the snake's mitogenome. Here, we analyzed the complete mitogenome sequences and structures of 65 snake species from 14 families and examined their structural patterns, organization and evolution. Our purpose was to further investigate the evolutionary implications and possible rearrangement mechanisms of the mitogenome within snakes.

RESULTS

In total, eleven types of mitochondrial gene arrangement patterns were detected (Type I, II, III, III-A, III-B, III-B1, III-C, III-D, III-E, III-F, III-G), with mitochondrial genome rearrangements being a major trend in snakes, especially in Alethinophidia. In snake mitogenomes, the rearrangements mainly involved three processes, gene loss, translocation and duplication. Within Scolecophidia, the O_L was lost several times in Typhlopidae and Leptotyphlopidae, but persisted as a plesiomorphy in the Alethinophidia. Duplication of the control region and translocation of the tRNA^Leu gene are two visible features in Alethinophidian mitochondrial genomes. Independently and stochastically, the duplication of pseudo-Pro (P*) emerged in seven different lineages of unequal size in three families, indicating that the presence of P* was a polytopic event in the mitogenome.

CONCLUSIONS

The WANCY tRNA gene cluster and the control regions and their adjacent segments were hotspots for mitogenome rearrangement. Maintenance of duplicate control regions may be the source for snake mitogenome structural diversity.

The complete mitochondrial genome of Gymnogobius Petschiliensis (Gobiiformes; Gobiidae; Gobionellinae) and its phylogenetic implications

Of the goby fishes, many Gymnogobius species have been poor recognized mainly because of the absence of enough molecular information and clear phylogenetic framework. In this study, the complete mitochondrial genome of Gymnogobius petschiliensis was determined and described. The mitogenome is 16,422 bp in length and consists of 22 tRNAs, 13 protein-coding genes, two rRNAs, one control region and a light strand replication origin (O_L). The arrangement of this mitogenome is identical to that of the typical teleost. The overall base composition is 27.5%, 29.5%, 26.1%, and 16.9%, for A, T, C, and G, respectively, with a slight bias on A + T content (57.0%). The 13 protein-coding genes use the initiation codon ATG except COI, which uses GTG. Most of them use TAA or TAG as the stop codon, while COII, COIII, and Cyt b use an incomplete T or TA and ND4 uses an unusual AGA. The maximum-likelihood phylogeny tree of 19 Gobionellinae species demonstrated that G. petschiliensis had a very closely relationship with the same genus G. urotaenia. This study is expected to contributing to the phylogenetic evolution of G. petschiliensis and further phylogenetic relationship of Gobionellinae and Gobiiformes.

The complete mitochondrial genome characterization of Tridentiger obscurus (Gobiiformes: Gobiidae) and phylogenetic analyses of Gobionellinae

The dusky tripletooth goby, Tridentiger obscurus, is a good model organism for the small in size and reaching maturity within a single year. Previous studies mainly focused on the annual reproductive cycle, social behavior and life history, but little information is available of the mitochondrial genome and phylogenetic evolution of this gobioid fish. In this article, we described the complete mitogenome of T. obscurus and reconstructed the phylogenetic relationship of the relative species of Gobionellinae. The genome is 16,501 bp in length including 13 protein-coding, two ribosomal RNA, 22 transfer RNA genes, as well as a putative control region and an L-strand replication origin. The overall base composition is 28.1%, 27.0%, 28.0% and 16.9% for A, T, C and G, respectively. This result is expected for better understanding the systematic evolution of the genus Tridentiger and further phylogenetic study of Gobiiformes.

Development of a Method to Extract and Amplify the Complete Mitogenome of Some Sparidae Species

Previous studies showed that fish mitochondrial DNA (mtDNA) is set up by a closed circular molecule of 16-17 kilobases (kb), comprising 2 ribosomal RNA genes (rRNA), 22 transfer RNA genes (tRNA), 13 protein-coding genes and 2 non-coding regions. The analysis of single mtDNA genes, such as Cytb, COI, 16S and 12S, or short segment of them, has been widely used against species substitution in both fresh and processed fish products. The analysis of the complete mitochondrial genome of fishery products allows to better study and characterise fish species. The aim of this research was to extract and amplify the complete mtDNA of some fish species of commercial interest belonging to the Sparidae family. The studied species were Dentex dentex, Dentex gibbosus, Dentex nufar, Pagellus acarne and Pagellus erythrinus. The entire mitogenome was obtained by gene amplification using long polymerase chain reactions. The analysis of the complete mitochondrial sequences will allow to gain further insights on these species and to find polymorphic sites that assess the degree of genetic variability of the species belonging to the family Sparidae.

Characterization of the complete mitochondrial genome of Cynoglossus gracilis and a comparative analysis with other Cynoglossinae fishes

Mitochondrial genomes can provide basic information for phylogenetic analysis and evolutionary studies. We present here the mitochondrial genome of Cynoglossus gracilis, which is 16,565bp in length. Numerous distinct regions were identified, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, a light-strand replication origin, and a control region. Interestingly, we detected rearrangement of genes in C. gracilis, including a control region translocation, tRNA(Gln) gene inversion, and tRNA(Ile) gene shuffling. Additionally, a phylogenetic analysis based on the nucleotide sequences of the 13 PCGs using maximum likelihood and Bayesian inference methods reveals that C. gracilis is closely related to Cynoglossus semilaevis. This study provides important mitogenomic data for analyzing phylogenetic relationships in the Cynoglossinae.

Tandem Duplication and Random Loss for mitogenome rearrangement in Symphurus (Teleost: Pleuronectiformes)

Background

The mitochondrial genomes (mitogenomes) of flatfishes (Pleuronectiformes) exhibit highly diversified types of large-scale gene rearrangements. We have reported that the mitogenomes of Crossorhombus azureus (Bothidae), Samariscus latus (Samaridae) and Cynoglossus fishes (Cynoglossidae) show different types of gene rearrangements.

Results

In the present study, the complete mitogenomes of two Symphurus species (Cynoglossidae), Symphurus plagiusa and Symphurus orientalis, were determined. The gene order in the S. plagiusa mitogenome is the same as that of a typical vertebrate (without any gene rearrangements). Surprisingly, large-scale gene rearrangements have occurred in S. orientalis. In the rearranged fragment from the control region (CR) to the WANCY tRNA cluster (tRNA cluster of tRNA-W, tRNA-A, tRNA-N, tRNA-C and tRNA-Y) in the S. orientalis mitogenome, tRNA-V and tRNA-M have been translocated to the 3’ end of the 16S rRNA gene, with six large intergenic spacers over 20 bp in length. In addition, an origin for light-strand replication (O_L) structure that is typically located in the WANCY region was absent in both the S. plagiusa and S. orientalis mitogenomes. It is generally recognized that a sequence in the WANCY region that encodes tRNAs forms a hairpin structure (O_L-like structure) and can act as the O_L when the typical locus is lost. Moreover, an additional O_L-like structure was identified near the control region in the S. plagiusa mitogenome.

Conclusions

The positions of the intergenic spacers and the rearranged genes of the S. orientalis mitogenome strongly indicate that the mechanism underlying the rearrangement of this mitogenome was Tandem Duplication and Random Loss. Additionally, two O_L-like regions substituting for the typical locus were found in the S. plagiusa mitogenome. We speculate that the ancestral mitogenomes of S. plagiusa and S. orientalis also had this characteristic, such that if both O_L-like structures functioned during mitochondrial replication, they could initiate duplicate replications of the light strand (L-strand), leading to duplication of the region between the two structures. We consider that this mechanism may account for the gene duplication that occurred during the gene rearrangement process in the evolution of the ancestral mitogenome to the S. orientalis mitogenome.

Electronic supplementary material

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

The complete mitochondrial genome sequence of Hemibagrus sp. (Siluriformes: Bagridae), molecular data for species identification

In this study, we determined and described the complete mitogenome sequence of Hemibagrus sp. for the first time, which is 16,513 bp in length, and contains 37 genes, including 13 protein-coding genes, 2 rRNAs, 22 tRNAs, 1 origin of replication on the light-strand (OL) and a putative control region. The overall base composition was 31.1% A, 26.9% T, 26.9% C, 15.1% G, with a slight AT bias (58.0%). All protein-coding genes shared the start codon ATG, except for COI, which began with GTG. The tRNA-Ser(UGC) couldn't be folded into the typical cloverleaf secondary structure because its dihydrouridine arm is replaced by a simple loop. These results are expected to provide useful molecular data for species identification and further phylogenetic studies of Bagridae and Siluriformes.

The complete mitochondrial genome of the gilthead seabream Sparus aurata L. (Sparidae)

The complete mitochondrial genome of the gilthead seabream Sparus aurata Linnaeus 1758, one of the world's most important mariculture species, was sequenced using next generation sequencing technology. The genome sequence is comprised of 16,652 bp exhibiting the canonical vertebrate mitochondria gene order. Regions of gene overlap, tRNA length, as well as start and stop codon were similar to those observed in other Sparidae. Phylogenetic reconstructions based on mitochondrial protein coding genes corroborate the view that Sparidae is paraphyletic and includes Centracanthidae.

A novel model of double replications and random loss accounts for rearrangements in the Mitogenome of Samariscus latus (Teleostei: Pleuronectiformes)

BACKGROUND

Although more than one thousand complete mitochondrial DNA (mtDNA) sequences have been determined in teleostean fishes, only a few gene rearrangements have been observed, and genome-scale rearrangements are even rarer. However, flatfishes (Pleuronectiformes) have been identified as having diverse types of mitochondrial gene rearrangements. It has been reported that tongue soles and the blue flounder mitogenomes exhibit different types of large-scale gene rearrangements.

RESULTS

In the present study, the complete mitochondrial genome of another flatfish, Samariscus latus, was sequenced, and genome-scale rearrangements were observed. The genomic features of this flounder are different from those of any other studied vertebrates, including flatfish species too. The mitogenome of S. latus is characterized by the duplication and translocation of the control region (CR). The genes located between the two CRs are divided into two clusters in which their relative orders are maintained.

CONCLUSIONS

We propose a "Double Replications and Random Loss" model to explain the rearrangement events in S. latus mitogenome. This model consists of the following steps. First, the CR was duplicated and translocated. Subsequently, double replications of the mitogenome were successively initiated from the two CRs, leading to the duplication of the genes between the two CRs. Finally, one of each pair of duplicated genes was lost in a random event.

The complete mitochondrial genome of Japanese Ammodytes personatus (Perciformes, Ammodytidae)

The mitochondrial genome of Japanese Ammodytes personatus is a circular molecule of 16,527 bp in length, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 non-coding regions (control region and L-strand replication origin OL). The mitogenome exhibits the canonical gene arrangement as those of most vertebrates, the nucleotide composition of which reveals an obvious anti-G bias. In the control region, we identified the termination-associated sequence domain, the central conserved sequence block domains (F, E, D, C, B and A), and CSB domains (1, 2 and 3).

Complete sequence and characterization of mitochondrial genome in great snakehead, Channa marulius (Hamilton, 1822)

The complete nucleotide sequence of mitogenome of the great snakehead, Channa marulius (Channidae), was determined and found to be 16,569 base pairs in length. The content and arrangement of different genes on the mitogenome was found similar to other typical teleosts. The overall base composition of the L-strand was found to be T (19.1%), C (31.5%), A (34.8%) and G (14.6%). The control region was 915 nt long and without any repetitive region. The mitogenome sequence data would be useful for studying phylogenetic relationship of C. marulius with other perciform species.

Complete mitogenome sequences of four flatfishes (Pleuronectiformes) reveal a novel gene arrangement of L-strand coding genes

Background

Few mitochondrial gene rearrangements are found in vertebrates and large-scale changes in these genomes occur even less frequently. It is difficult, therefore, to propose a mechanism to account for observed changes in mitogenome structure. Mitochondrial gene rearrangements are usually explained by the recombination model or tandem duplication and random loss model.

Results

In this study, the complete mitochondrial genomes of four flatfishes, Crossorhombus azureus (blue flounder), Grammatobothus krempfi, Pleuronichthys cornutus, and Platichthys stellatus were determined. A striking finding is that eight genes in the C. azureus mitogenome are located in a novel position, differing from that of available vertebrate mitogenomes. Specifically, the ND6 and seven tRNA genes (the Q, A, C, Y, S₁, E, P genes) encoded by the L-strand have been translocated to a position between tRNA-T and tRNA-F though the original order of the genes is maintained.

Conclusions

These special features are used to suggest a mechanism for C. azureus mitogenome rearrangement. First, a dimeric molecule was formed by two monomers linked head-to-tail, then one of the two sets of promoters lost function and the genes controlled by the disabled promoters became pseudogenes, non-coding sequences, and even were lost from the genome. This study provides a new gene-rearrangement model that accounts for the events of gene-rearrangement in a vertebrate mitogenome.

A comprehensive description and evolutionary analysis of 22 grouper (perciformes, epinephelidae) mitochondrial genomes with emphasis on two novel genome organizations

Groupers of the family Epinephelidae are a diverse and economically valuable group of reef fishes. To investigate the evolution of their mitochondrial genomes we characterized and compared these genomes among 22 species, 17 newly sequenced. Among these fishes we identified three distinct genome organizations, two of them never previously reported in vertebrates. In 19 of these species, mitochondrial genomes followed the typical vertebrate canonical organization with 13 protein-coding genes, 22 tRNAs, two rRNAs, and a non-coding control region. Differing from this, members of genus Variola have an extra tRNA-Ile between tRNA-Val and 16S rRNA. Evidence suggests that this evolved from tRNA-Val via a duplication event due to slipped strand mispairing during replication. Additionally, Cephalopholisargus has an extra tRNA-Asp in the midst of the control region, likely resulting from long-range duplication of the canonical tRNA-Asp through illicit priming of mitochondrial replication by tRNAs. Along with their gene contents, we characterized the regulatory elements of these mitochondrial genomes' control regions, including putative termination-associated sequences and conserved sequence blocks. Looking at the mitochondrial genomic constituents, rRNA and tRNA are the most conserved, followed by protein-coding genes, and non-coding regions are the most divergent. Divergence rates vary among the protein-coding genes, and the three cytochrome oxidase subunits (COI, II, III) are the most conserved, while NADH dehydrogenase subunit 6 (ND6) and the ATP synthase subunit 8 (ATP8) are the most divergent. We then tested the phylogenetic utility of this new mt genome data using 12 protein-coding genes of 48 species from the suborder Percoidei. From this, we provide further support for the elevation of the subfamily Epinephelinae to family Epinephelidae, the resurrection of the genus Hyporthodus, and the combination of the monotypic genera Anyperodon and Cromileptes to genus Epinephelus, and Aethaloperca to genus Cephalopholis.

Structure comparison of control region in Stromateoidei and complete mitochondrial genome of the bluefin driftfish Psenes pellucidus (Perciformes, Nomeidae)

In the present study, the complete mitochondrial genome of the bluefin driftfish Psenes pellucidus was determined first. The genome was 16,572 bp in length and consisted of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes and 2 main non-coding regions. The mitochondrial genome of bluefin driftfish had common features about gene arrangement, base composition and tRNA structures compared with those of most of the bony fishes. In the control region, a termination-associated sequence, the central conserved block (CSB-F, CSB-E and CSB-D) and the conserved sequence blocks (CSB-1, CSB-2 and CSB-3) were determined. Meanwhile, the conserved motif 5'-GCCGG-3' was identified in the origin of light strand replication of bluefin driftfish. This mitogenome sequence data would play key role in phylogenetic analysis of Stromateoidei.

The complete mitochondrial genome of black grunt Hapalogenys nigripinnis

The complete sequence of the mitochondrial genome of Hapalogenys nigripinnis was determined in this study. The genome was 16,478 bp in size with a base composition of 28.6% A, 15.8% G, 27.4% T, and 28.2% C, containing a typical structure of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 noncoding control region. The genomic composition, organization, and gene order of H. nigripinnis was similar to that obtained in most vertebrates. These results may provide molecular information on the future phylogenetic relationships of H. nigripinnis and its position within the suborder Percoidei.

Phylogenetic studies of three sinipercid fishes (Perciformes: Sinipercidae) based on complete mitochondrial DNA sequences

The sinipercids are a group of 12 species of freshwater percoid fish endemic to East Asia and their phylogenetic placements have perplexed generations of taxonomists. We cloned and sequenced the complete mitochondrial DNA (mtDNA) of three sinipercid fishes (Siniperca chuatsi, S. kneri, and S. scherzeri) to characterize and compare their mitochondrial genomes. The mitochondrial genomes of S. chuatsi, S. kneri, and S. scherzeri were 16,496, 17,002, and 16,585 bp in length, respectively. The organization of the three mitochondrial genomes is similar to those reported from other fish mitochondrial genomes, which contains 37 genes (13 protein-coding genes, 2 ribosomal RNAs, and 22 transfer RNAs) and a major non-coding control region. Among the 13 protein-coding genes of all the three sinipercid fishes, three reading-frame overlaps were found on the same strand. There is an 81-bp tandem repeat cluster at the end of CSB-3 in the S. scherzeri control region. The complete mitochondrial genomes of the three sinipercids should be useful for the evolutionary studies of sinipercids and other vertebrate species.

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus+Collichthys lucidus and Larimichthys crocea+Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.

Complete mitochondrial genome of three Branchiostegus (Perciformes, Malacanthidae) species: genome description and phylogenetic considerations

We cloned and sequenced the complete mitochondrial DNA (mtDNA) of three tilefishes (Branchiostegus albus, Branchiostegus argentatus, and Branchiostegus japonicus) to characterize and compare their mitochondrial genomes (mitogenomes). The mitogenomes of B. albus, B. argentatus, and B. japonicus were 16,532, 16,550, and 16,541 bp long, respectively, and all consisted of 37 genes (13 protein-coding genes, 2 ribosomal RNA, and 22 transfer RNA (tRNAs)), which are typical for vertebrate mtDNA. As in other bony fishes, most genes were encoded on the H-strand, except for the nad6 and eight tRNA genes that were encoded on the L-strand. Among the 13 protein-coding genes of all three tilefishes, 2 reading-frame overlaps were found on the same strand: atp8 and atp6 overlapped by 10 nucleotides, and nad4L and nad4 overlapped by 7 nucleotides. The identity of the nad4 gene between B. albus and B. argentatus was the lowest at 87%. Conversely, the identity of the nad6 gene between B. albus and B. japonicus was the highest at 99%. Most tRNA genes were similar in length among the three species, while the tRNA-Ser((AGY)) of B. japonicus was 9 bp longer than those of B. albus and B. argentatus. The control region of the mitogenome spanned 853, 862, and 856 bp in B. albus, B. argentatus, and B. japonicus, respectively. A maximum likelihood tree constructed using 11,035 sites contained five independent groups with bootstrap values of 100% in support of their divergence. All three tilefishes examined were clustered with the Pomacanthidae species in Group II.

The mitochondrial genome of the spinyhead croaker Collichthys lucida: genome organization and phylogenetic consideration

The complete mitochondrial genome of the spiny head croaker Collichthys lucida was determined in the present study. The mitochondrial DNA was 16,442 base pairs in length, and contained 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and one major non-coding control region, with the content and order of genes being similar to those in typical teleosts. Most of the genes of C. lucida were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro) genes were encoded on the L-strand. The reading frames of two pairs of genes overlapped: ATPase 8 and 6 and ND4L and ND4 by ten and seven nucleotides, respectively. The control region was unusually short at only 768bp, and absence of typical conserved blocks (CSB-D, CSB-E, and CSB-F). Phylogenetic analyses indicated that C. lucida was located in the cluster of fish species from the family Sciaenidae, supporting the traditional taxonomic classification of fish, and in the cluster of Serranidae, the divergence time in Plectropomus leopardus is longer than that among its coordinal species. On the other hand, phylogenetic analyses do not support the monophyletic of family Centracanthidae and genera Larimichthys and Collichthys, which is against the morphological results.

Evolutionary relatedness of mackerels of the genus Scomber based on complete mitochondrial genomes: strong support to the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as distinct species

Mackerels of the genus Scomber are commercially important species, but their taxonomic status is still controversial. Although previous phylogenetic data support the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as separate species, it is only based on the analysis of partial mitochondrial and nuclear DNA sequences. In an attempt to shed light on this relevant issue, we have determined the complete mitochondrial DNA sequence of S. colias, S. japonicus, and Scomber australasicus. The total length of the mitogenomes was 16,568 bp for S. colias and 16,570 bp for both S. japonicus and S. australasicus. All mitogenomes had a gene content (13 protein-coding, 2 rRNAs, and 22 tRNAs) and organization similar to that observed in Scomber scombrus and most other vertebrates. The major noncoding region (control region) ranged between 865 and 866 bp in length and showed the typical conserved blocks. Phylogenetic analyses revealed a monophyletic origin of Scomber species with regard to other scombrid fish. The major finding of this study is that S. colias and S. japonicus were significantly grouped in distinct lineages within Scomber cluster, which phylogenetically constitutes evidence that they may be considered as separate species. Additionally, molecular data here presented provide a useful tool for evolutionary as well as population genetic studies.

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.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae): Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796bp in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O(L)) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T. nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T. japonicus together with T. brevis, T. lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae.

The complete mitochondrial genome of the large yellow croaker, Larimichthys crocea (Perciformes, Sciaenidae): unusual features of its control region and the phylogenetic position of the Sciaenidae

To understand the systematic status of Larimichthys crocea in the Percoidei, we determined the complete mitochondrial (mt) genome sequence using 454 sequencing-by-synthesis technology. The complete mt genome is 16,466 bp in length including the typical structure of 22 tRNAs, 2 rRNAs, 13 protein-coding genes and the noncoding control region (CR). Further sequencing for the complete CR was performed using the primers Cyt b-F and 12S-R on six L. crocea individuals and two L. polyactis individuals. Interestingly, all seven CR sequences from L. crocea were identical while the three sequences from L. polyactis were distinct (including one from GenBank). Although the conserved blocks such as TAS and CSB-1, -2, and -3 are readily identifiable in the control regions of the two species, the typical central conserved blocks CSB-D, -E, and -F could not be detected, while they are found in Cynoscion acoupa of Sciaenidae and other Percoidei species. Phylogenetic analysis shows that L. crocea is a relatively recently emerged species in Sciaenidae and this family is closely related to family Pomacanthidae within the Percoidei. L. crocea, as the first species of Sciaenidae with complete mitochondrial genome available, will provide important information on the molecular evolution of the group. Moreover, the genus-specific pair of primers designed in this study for amplifying the complete mt control region will be very useful in studies on the population genetics and conservation biology of Larimichthys.