Next Article in Journal
Ketogenic Diet and Microbiota: Friends or Enemies?
Previous Article in Journal
Maternally Inherited Differences within Mitochondrial Complex I Control Murine Healthspan
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 7
10.3390/genes10070533
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The First Complete Mitochondrial Genome of the Flathead Cociella crocodilus (Scorpaeniformes: Platycephalidae) and the Phylogenetic Relationships within Scorpaeniformes Based on Whole Mitogenomes

by
Lei Cui
1,†,
Rongbo Cao
1,†,
Yuelei Dong
1,
Xingchen Gao
2,
Jingyi Cen
1,* and
Songhui Lu
1,*
1
Key Laboratory of Eutrophication and Red Tide Prevention, Research Center for Harmful Algae and Marine Biology, Jinan University, Guangzhou 510632, China
2
Chinese Sturgeon Research Institute, Three Gorges Corporation, Yichang 443100, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2019, 10(7), 533; https://0-doi-org.brum.beds.ac.uk/10.3390/genes10070533
Submission received: 21 May 2019 / Revised: 12 July 2019 / Accepted: 12 July 2019 / Published: 15 July 2019
(This article belongs to the Section Population and Evolutionary Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Complete mitochondrial genomes (mitogenomes) are important molecular markers for understanding the phylogenetics of various species. Although recent studies on the mitogenomes of the Scorpaeniformes species have been greatly advanced, information regarding molecular studies and the taxonomic localization of Platycephalidae is still sparse. To further analyze the phylogeny of Platycephalidae, we sequenced the complete mitogenome of Cociella crocodilus of the Platycephalidae family within Scorpaeniformes for the first time. The mitogenome was 17,314 bp in length, contained two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and two typical noncoding control regions (the control region (CR) and origin of the light strand (OL)). All PCGs used standard initiation codons ATG, apart from cox1. The majority of the tRNA genes could be folded into cloverleaf secondary structures, whereas the secondary structure of tRNASer (AGN) lacked a dihydrouridine (DHU) stem. The CR contained several conserved sequence blocks (CSBs) and eight tandem repeats. In addition, the phylogenetic relationship based on the concatenated nucleotides sequences of 13 PCGs indicated that the Platycephalidae species are relatively basal in the phylogenetic relationships of Scorpaeniformes. Our results may not only advance the origin and the evolution of Scorpaeniformes, but also provide information for the genetic evolution and taxonomy of the teleostean species.

1. Introduction

Within the order of Scorpaeniformes, the Platycephalidae family is usually distributed in tropical and temperate coastal or estuarine waters in the Indian and Pacific Oceans. Most of the species live in typical sandy and muddy benthic habitats [1,2]. Platycephalidae, most of which have important commercial and medicinal value [3], consists of approximately 70 valid species of 17 genera [4]. The identification of flatheads usually depends on morphological characteristics such as an elongated body with a spiny depressed head and a large mouth [5]. With the development of molecular biology and molecular genetics, DNA markers (e.g., nuclear and mitochondrial genes) have been widely applied in studies of classification, phylogenetics and genetic diversity of fish.
The typical mitochondrial genome (mitogenome) is a circular, double-stranded molecule, that ranges in size from 15 kb to 19 kb, and generally contains 37 genes (13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA genes (tRNAs) and two noncoding regions (the control region (CR) and origin of the light strand (OL)), and has significant functions in the regulation and initiation of mitogenome transcription and replication [6,7]. In recent years, mitogenomes have been widely used as effective molecular markers for molecular evolution, phylogenetic studies and population genetics [8,9] due to their unique characteristics of maternal inheritance, relatively conserved genome structures, coding content conservation, high mutation rate and low intermolecular genetic recombination level [6]. Based on these advantages, the complete mitogenome of genetic code alteration, gene rearrangement, secondary structures of tRNAs and regions of transcription and replication are also widely utilized for taxonomic phylogenetic analysis [10,11].
To date, complete mitogenomes have been sequenced for approximately 39 Scorpaeniformes species. However, only one complete mitogenome from Platycephalidae, Platycephalus indicus, has been sequenced. We previously used the mitochondrial genome to study the genetic relationships of Scorpaeniformes [12], which does not include the Platycephalidae family, although phylogenetic trees based on nuclear genes and partial mitochondrial gene sequences were constructed to analyze interspecies relationships with Scorpaeniformes in previous research. Nevertheless, the difference between whether flatheads or scorpionfishes are ‘the most basal clade’ of Scorpaeniformes was studied [13,14,15]. The markers of small fragments may not provide sufficient evidence to explain phylogenetic relationships. To research higher-level relationships among Scorpaeniformes, the complete mitogenome of one Platycephalidae species, Cociella crocodilus, was sequenced in our study and we analyzed the gene content and organization compared with other Scorpaeniformes species. In addition, the phylogenetic tree based on the thirteen PCG sequences was reconstructed by Maximum Likelihood (ML) and Bayesian inference (BI) methods to understand the higher phylogeny of Scorpaeniformes.

2. Methods

2.1. Sampling and DNA Extraction

A specimen of C. crocodilus was captured in the Pearl River estuary (N 21°45′, E 133°36′), China, in July 2016.The field study did not involve endangered or protected species, in accordance with the IUCN Red List. The research was approved by the Institutional Animal Care and Use Committee at Jinan University, with no ethical code associated. The experiments were conducted in accordance with international guidelines for the care and treatment of laboratory animals. The samples were identified according to morphological characteristics and were then conserved in 95% ethanol and stored at −80 °C until processing. Dorsal muscle tissue was removed to extract the total genomic DNA using the Animal Tissue Genomic DNA Extraction Kit (SangonBiotech, Shanghai, China) according to the manufacturer’s instructions. The extracted DNA was used to amplify the complete mitogenome of C. crocodilus by polymerase chain reaction (PCR).

2.2. Polymerase Chain Reaction Amplification and Sequencing

To obtain the complete C. crocodilus mitogenome sequence, PCR amplification was performed with several primer pairs designed based on aligned mitogenome sequences of Chelidonichthys kumu (GenBank: KY379222) and Lepidotrigla microptera (GenBank: KY012348.1) (Table S1) [12]. All PCR reactions were conducted with LA Taq DNA polymerase using Premix LA Taq (Takara, Dalian, China) according to the following cycling conditions: beginning with an initial denaturation step at 95 °C for 1 min, followed by 35 cycles of denaturation at 95 °C for 20 s, annealing at 55 °C for 45 s and elongation at 72 °C for 1–5 min depending on the length of the segment. The sequences of the PCR-amplified products were determined on a 3730XL DNA Analyzer (Beijing Genomics Institute, Shenzhen China).

2.3. Sequence Analysis

In order to obtain the final complete sequence, the obtained sequenced fragments were assembled through the program Seqman within the Lasergene software and then manually checked [16]. Sequence annotation was performed using NCBI BLAST (http://blast.ncbi.nlm.nih.gov/Blast) and the DNASTAR package (DNASTAR Inc, Madison, WI, USA). The tRNA genes were identified using the default search mode of the tRNAscan-SE program [17]. Secondary structures of tRNA genes and OL were inferred using the software RNAstructure [18] with vertebrate mitochondrial predictors. The Dual Organellar GenoMe Annotator (DOGMA) automated the annotation of two rRNA and 13 PCGs [19]. The codon usage of the protein-coding genes and the nucleotide composition were determined using MEGA 6.0 [20]. The tandem repeats in the control region were predicted using the Tandem Repeats Finder available online [21]. The formulas AT-skew [(A − T)/(A + T)] and GC-skew [(G − C)/(G + C)] were used to measure the nucleotide composition skewness [22]. The final mitogenome sequence of the C. crocodilus was submitted to the GenBank database under the accession number MH521260.1.

2.4. Phylogenetic Analysis

Based on 39 Scorpaeniformes mitogenomes available in the GenBank, phylogenetic analysis was performed (Table 1). Three outgroups taxa were assigned, i.e., Perciformes fishes (Percalates novemaculeata (NC_024850.1) and Siniperca knerii (NC_015987.1) [23,24]. The concatenated nucleotide sequence alignment was used in IQ-TREE and MrBayes v 3.2.4 to determine maximum likelihood (ML) and Bayesian inference (BI) for phylogenetic analysis [25,26]. All the sequences of the genes were aligned individually by ClustalX [27]. GTR+F+R6 was selected as the appropriate model for the nucleotide sequences by ModelFinder in IQ-TREE based on Akaike’s information criterion (AIC) [25,28]. Four independent Markov chains were simultaneously used at 1,000,000 generations with sampling every 1000 generations for BI analyses. The first 25% were discarded as burn-in. The BI Tree was considered to be reached since the average standard deviation of the split frequencies was below 0.01. ML analysis was carried out with 100 bootstrap replicates based on the default parameters.

3. Results and Discussion

3.1. Mitogenome Organization, Structure and Skewness

The complete mitogenome of C. crocodilus formed a closed circular molecule with a total length of 17,314 bp in size, which was the largest of all the sequenced Scorpaeniformes species available in the GenBank to date. The mitogenome contained 37 genes identical to other Scorpaeniformes species: two rRNA genes (12S rRNA and 16S rRNA), 22 tRNAs (one for each amino acid, two for Leucine and Serine), 13 PCGs (cox1-3, nad1-6, nad4L, atp6, atp8 and cytb) and two noncoding regions (control region (CR) and OL) (Figure S1 and Table 2). Except ND6 and eight tRNA genes (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro), which were distributed on the light strand (L-strand), most of the other mitochondrial genes were distributed on the heavy strand (H-strand). The mitogenome of C. crocodilus had a total overlap of 29 bp (ranging from 1 to 10 bp) between genes at nine locations, among which atp8- atp6 (10 bp) and nad4L-nad4 (7 bp) often occurred in the Metazoa. A total of 162 bp nucleotides ranged from 1 to 79 bp in 15 spacers. The longest 79 bp spacer was between tRNAL1 and nad1 (Table 2). The nucleotide composition of the C. crocodilus is as follows: A = 29.14%, T = 27.19%, G = 15.91% and C = 27.76%, which shows a high A+T content of 56.3%. The A+T content is higher than most sequenced Scorpaeniformes species, such as Sebates schlegeli (53.8%), Scorpaenopsis cirrosa (54.27%), and Hexagrammos lagocephalus (53.27%). The exception is Synanceia verrucosa (59.35%). The fact that the A + T content was higher than the G+C content is in accordance with the characteristics of the mitogenome of metazoan animals [29]. For the complete mitogenome sequence of Cociella crocodilus, the AT skew was positive (0.0348) and the GC skew was negative (−0.2713). Generally, the teleostean mitogenomes tend to have identical positive AT-skew and negative GC-skew. All of the GC skew was negative among Scorpaeniformes species and most of the AT skew was positive, apart from Satyrichthys amiscus (−0.0058), Pleurogrammus monopterygius (−0.0012), Cottus hangiongensis (−0.0011), and Cottus poecilopus (−0.0075) [30,31].

3.2. Protein-Coding Genes

The complete mitogenome of C. crocodilus had 13 PCGs with a total length of 11,445 bp, 12 of which utilize ATG for the start codon, whereas cox1 utilizes GTG. In teleost mitogenomes, initiator codons are often used selectively (e.g., cox1 of the whiting, Merlangius merlangus; cox1 of the haddock, Melanogrammus aeglefinus; and cox1 and nad4l of the Globose Head Whiptail, Cetonurus globiceps) [32,33]. Eight of the 13 PCGs (nad1, nad2, cox1, atp8, atp6, cox3, nad4l, and nad5) used the typical stop codon TAA. The two genes nad3 and nad6 were terminated with TAG. The incomplete termination codon T was used by cox2, nad4, and cytb (Table 2). Such incomplete termination codons were due to the 3’ end of the mRNA appearing to have posttranscriptional polyadenylation [34]. Except for the AGA and AGG termination codons, the total number of codons in the PCGs of C. crocodilus was 3,804, among which three amino acids (Leu1, 513, Ala, 353, and Thr, 302) were the most frequent (Figure 1 and Table S2). In addition to the three genes of C. crocodilus, cox2, nad4 and nad5, most of the AT skew was negative. Most of the GC skew values were slightly negative, except for nad6, indicating that there were more Ts and Cs in most PCGs, which is consistent with previous research (Figure 2).

3.3. Transfer RNA Genes and Ribosomal RNA Genes

There were 22 tRNA genes ranging from 67 bp (trnP) to 76 bp (trnL1) identified in the C. crocodilus mitogenome. Among them, 14 tRNA genes were encoded on the H-strand, and the other eight tRNA genes were coded by the L-strand (Table 2). Most of the tRNAs were folded into an ordinary cloverleaf secondary structure (Figure S2), except trnS2 (AUN), which lacked a DHU stem. This phenomenon is identical to other mitogenomes of teleost fishes, including the Scorpaeniformes species [12,35]. The standard structure of the amino acid acceptor stem was 7 bp, while trnaV, trnF, and trnaL1 (UUR) were 9 bp. Further, there were 24 base-pair mismatches present in the tRNA secondary structures, all of which formed a weak bond (G-U pairs). The compositions of the complete sequences of all 22 tRNAs showed both slightly positive AT skew (0.0136) and GC skew (0.0473), indicating that tRNAs biased toward As and Gs. The 12S and 16S rRNA genes were 1,021 bp and 1,713 bp in length, respectively, and similarly, the AT skew (0.0763 and 0.2227) was slightly positive and the GC skew (−0.1952 and −0.1064) was negative. They were situated between trnF and trnL1 and there is a trnV between rRNAs similar to other vertebrate mitogenomes to date [36,37,38].

3.4. Noncoding Regions

There were two long noncoding regions in the mitogenome of C. crocodilus, OL and CR, which were closely related to the initiation of the replication and transcription of the mitogenome [39,40]. Located between trnN and trnC with a length of 30 bp, the relatively short gene OL was folded into a hairpin secondary structure (Figure S3). The CR was the longest noncoding region and extended over 1447 bp between trnP and trnF, with a 68.76% A + T content, which is the highest among Scorpaeniformes species to date. Similarly to other teleost fish, there were several conserved sequence blocks (CSBs) identified in CR, which were shown to play indispensable roles in mitochondrial metabolism [41]. The central conserved domain contained blocks CSB-D, -E and -F, while the typical conserved sequence block region involved blocks CSB-1, CSB-2 and CSB-3 (Figure S4). Consistent with other Scorpaeniformes species, the conserved motifs ATGTA and its complement TACAT as the recognition sites are commonly found in the CR of teleosts. In addition, eight concatenated repeats in C. crocodilus were confirmed with the Tandem Repeats Finder program [21]. This phenomenon occurs only in Enophrys diceraus and Sebastiscus marmoratus, but never in other Scorpaeniformes species [42,43]. These tandem repeats may be one reason why the C. crocodilus mitogenome was bigger than in other sequenced Scorpaeniformes species.

3.5. Phylogenetic Relationships

To determine the phylogenetic position of C. crocodilus, 39 Scorpaeniformes species’ mitogenomes were obtained from the GenBank and the mitogenomes of Perciformes fishes (P. novemaculeata (NC_024850.1) and S. knerii (NC_015987.1) were used for outgroups (Figure 3). Phylogenetic trees inferred from the ML and BI methods were constructed for the nucleotide sequences of 13 PCGs. The results clearly indicate that flatheads, scorpionfishes and rockfishes together comprised Scorpaeniformes as a paraphyletic group. The phylogenetic trees revealed three strongly supported clades: (I) Synanceiidae and Platycephalidae, (II) Scorpaenidae and Sebastidae, and (III) Anoplopomatidae, Agonidae, Triglidae, Peristediidae, Hexagrammidae and Cottidae. The best supported phylogenetic relationships of Scorpaeniformes are as follows: (Synanceiidae + (Platycephalidae + ((Scorpaenidae + Sebastidae) + ((Triglidae + Peristediidae) + (Agonidae + (Hexagrammidae + (Anoplopomatidae + Cottidae))))))). There was a conflict with our work based on complete mitogenomes [12]. In a previous study, Synanceiidae was as a sister to Scorpaenidae + Sebastidae. The phylogenetic analysis results were similar to earlier studies (although only Synanceiidae, Scorpaenidae and Sebastidae were available) [44]. Based on the results of the current study, Synanceiidae was proposed to be one of the most basal taxa in the order of Scorpaeniformes after adding new taxa (Anoplopomatidae, Agonidae and Platycephalidae) to the phylogenetic analysis. As the only sequenced mitogenome to date, the Synanceiidae family requires further phylogenetic research. C. crocodilus, within Platycephalidae, which branches into an independent sub-branch, shares a close ancestry with P. indicus. The topology data obtained in this study indicates that the flatheads were more basal than scorpionfishes among Scorpaeniformes, which is not consistent with previous studies involving phylogenetic analysis on the basis of segmental mitogenomes (16S rRNA) and nuclear genes (rag1, rag2, and rh) [13]. Whether the difference in the phylogenetic analysis is due to introgression, hybridization and lineage sorting is unknown. However, some phylogenetic results also showed that flatheads were relatively basal species based on a modified time-calibrated model [14,15]. Additional research should be conducted to discover other un-sequenced families and employ other molecular datasets to elucidate phylogenetic relationships among Scorpaeniformes species.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/10/7/533/s1, Figure S1: Organization of the complete mitochondrial genome of C. crocodilus, Figure S2: Predicted secondary structures of 22 tRNA genes in the C. crocodilus mitogenome. Figure S3: The putative secondary hairpin structural features of the OL in the mitogenome of C. crocodilus, Figure S4: The structure of the control region in the mitogenome of C. crocodilus, Table S1: Primer pairs used for PCR amplification of C. crocodilus mitogenome, Table S2: Codon number and RSCU in C. crocodilus mitochondrial PCGs.

Author Contributions

Conceptualization, L.C., J.C. and S.L.; data curation, L.C., R.C. and Y.D.; formal analysis, L.C. and R.C.; funding acquisition, J.C. and S.L.; investigation, X.G. and J.C.; methodology, X.G.; resources, X.G. and J.C.; supervision, S.L.; validation, Y.D.; writing—original draft, L.C., R.C. and Y.D.; writing—review & editing, J.C. and S.L.

Funding

This work was supported by the National Natural Science Foundation of China (No. 41806127) and the Natural Science Foundation of Guangdong Province (No. 2018A030313956).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keenan, C.P. Phylogeny of Australian species of flatheads (Teleostei, Platycephalidae) as determined by allozyme electrophoresis. J. Fish Biol. 2010, 39, 237–249. [Google Scholar] [CrossRef]
  2. Dulčić, J.; Pallaoro, A.; Dragičević, B.; Stagličić-Radica, N. First record of dwarf flathead Elates ransonnetii (Platycephalidae) in the Adriatic Sea. Cybium 2010, 34, 222–223. [Google Scholar]
  3. Imamura, H. Phylogeny of the family platycephalidae and related taxa (Pisces: Scorpaeniformes). Species Divers. 1996, 1, 123–233. [Google Scholar] [CrossRef]
  4. Sriramachandra Murty, V.; Manikyam, Y. Taxonomic Revision of the Flatheads (Platycephalidae: Pisces) of India; Zoological Survey of India: Kolkata, India, 2007. [Google Scholar]
  5. Imamura, H.; Sakashita, M. Redescription of a rare flathead, Onigocia grandisquama (Regan, 1908) (Scorpaeniformes: Platycephalidae). Ichthyol. Res. 1997, 44, 119–124. [Google Scholar] [CrossRef]
  6. Boore, J.L. Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27, 1767–1780. [Google Scholar] [CrossRef] [Green Version]
  7. Cameron, S.L. Insect mitochondrial genomics: Implications for evolution and phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef]
  8. Simon, C.; Buckley, T.R.; Frati, F.; Stewart, J.B.; Beckenbach, A.T. Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 2006, 37, 545–579. [Google Scholar] [CrossRef]
  9. Miya, M.; Kawaguchi, A.; Nishida, M. Mitogenomic exploration of higher teleostean phylogenies: A case study for moderate-scale evolutionary genomics with 38 newly determined complete mitochondrial DNA sequences. Mol. Biol. Evol. 2001, 18, 1993–2009. [Google Scholar] [CrossRef]
  10. Boore, J.L.; Lavrov, D.V.; Brown, W.M. Gene translocation links insects and crustaceans. Nature 1998, 392, 667–668. [Google Scholar] [CrossRef]
  11. Dowton, M.; Castro, L.R.; Austin, A.D. Mitochondrial gene rearrangements as phylogenetic characters in the invertebrates: The examination of genome ‘morphology’. Invertebr. Syst. 2002, 16, 345–356. [Google Scholar] [CrossRef]
  12. Cui, L.; Dong, Y.; Liu, F.; Gao, X.; Zhang, H.; Li, L.; Cen, J.; Lu, S. The first two complete mitochondrial genomes for the family Triglidae and implications for the higher phylogeny of Scorpaeniformes. Sci. Rep. 2017, 7, 1553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Betancur, R.R.; Broughton, R.E.; Wiley, E.O.; Carpenter, K.; Ortí, G. The tree of life and a new classification of bony fishes. PLoS Curr. 2013, 5, 1001550. [Google Scholar] [CrossRef] [PubMed]
  14. Rabosky, D.L.; Chang, J.; Title, P.O.; Cowman, P.F.; Sallan, L.; Friedman, M.; Kaschner, K.; Garilao, C.; Near, T.J.; Coll, M.; et al. An inverse latitudinal gradient in speciation rate for marine fishes. Nature 2018, 559, 392–395. [Google Scholar] [CrossRef] [PubMed]
  15. Betancur, R.R.; Wiley, E.O.; Arratia, G.; Acero, A.; Bailly, N.; Miya, M.; Lecointre, G.; Ortí, G. Phylogenetic classification of bony fishes. BMC Evol. Biol. 2017, 17, 162. [Google Scholar] [CrossRef] [PubMed]
  16. Swindell, S.R.; Plasterer, T.N. SEQMAN. Contig assembly. Methods Mol Biol. 1997, 70, 75–89. [Google Scholar] [CrossRef] [PubMed]
  17. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  18. Mathews, D.H. Using the RNAstructure software package to predict conserved RNA structures. Curr. Protoc. Bioinform. 2014, 46, 12.14.11–12.14.22. [Google Scholar] [CrossRef]
  19. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [Green Version]
  20. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  21. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef]
  22. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 41, 353–358. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, D.X.; Chu, W.Y.; Liu, X.L.; Nong, X.X.; Li, Y.L.; Du, S.J.; Zhang, J.S. Phylogenetic studies of three sinipercid fishes (Perciformes: Sinipercidae) based on complete mitochondrial DNA sequences. Mitochondrial Dna 2012, 23, 70. [Google Scholar] [CrossRef] [PubMed]
  24. Lavoué, S.; Nakayama, K.; Jerry, D.R.; Yamanoue, Y.; Yagishita, N.; Suzuki, N.; Nishida, M.; Miya, M. Mitogenomic phylogeny of the Percichthyidae and Centrarchiformes (Percomorphaceae): Comparison with recent nuclear gene-based studies and simultaneous analysis. Gene 2014, 549, 46–57. [Google Scholar] [CrossRef] [PubMed]
  25. Lam-Tung, N.; Schmidt, H.A.; Arndt, V.H.; Bui Quang, M. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar]
  26. Ronquist, F.; Teslenko, M.; Mark, P.V.D.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed]
  27. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinform. 2002. [Google Scholar] [CrossRef] [PubMed]
  28. MrModeltest, V2, program distributed by the author; Nylander, J.A.A., Evolutionary Biology Centre at Uppsala University; Science Open; Berlin. 2004.
  29. Asakawa, S.; Kumazawa, Y.; Araki, T.; Himeno, H.; Miura, K.; Watanabe, K. Strand-specific nucleotide composition bias in echinoderm and vertebrate mitochondrial genomes. J. Mol. Evol. 1991, 32, 511–520. [Google Scholar] [CrossRef]
  30. Miya, M.; Takeshima, H.; Endo, H.; Ishiguro, N.B.; Inoue, J.G.; Mukai, T.; Shirai, S.M. Major patterns of higher teleostean phylogenies: A new perspective based on 100 complete mitochondrial DNA sequences. Mol. Phylogenetics Evol. 2003, 26, 121–138. [Google Scholar] [CrossRef]
  31. Hwang, D.S.; Byeon, H.K.; Lee, J.S. Complete mitochondrial genome of the freshwater sculpin Cottus hangiongensis (Scorpaeniformes, Cottidae). Mitochondrial Dna 2013, 24, 490–491. [Google Scholar] [CrossRef]
  32. Shi, X.; Peng, T.; Lin, R.; Huang, D.; Wang, J. Characterization of the complete mitochondrial genome sequence of the globose head whiptail Cetonurus globiceps (Gadiformes: Macrouridae) and its phylogenetic analysis. PloS ONE 2016, 11, e0153666. [Google Scholar] [CrossRef]
  33. Roques, S.; Fox, C.J.; Villasana, M.I.; Rico, C. The complete mitochondrial genome of the whiting, Merlangius merlangus and the haddock, Melanogrammus aeglefinus: A detailed genomic comparison among closely related species of the Gadidae family. Gene 2006, 383, 12–23. [Google Scholar] [CrossRef] [PubMed]
  34. Ojala, D.; Montoya, J.; Attardi, G. tRNA punctuation model of RNA processing in human mitochondrial. Nature 1981, 290, 470–474. [Google Scholar] [CrossRef] [PubMed]
  35. Lavrov, D.V.; Brown, W.M.; Boore, J.L. A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc. Natl. Acad. Sci. USA 2000, 97, 13738–13742. [Google Scholar] [CrossRef] [PubMed]
  36. Niu, S.; Liu, Y.; Qin, C.; Wang, X.; Wu, R. The complete mitochondrial genome and phylogenetic analysis of (Perciformes, Moronidae). Mitochondrial Dna 2015, 69, 1–3. [Google Scholar] [CrossRef] [PubMed]
  37. Chou, C.E.; Wang, F.Y.; Chang, S.K.; Chang, H.W. Complete mitogenome of three flyingfishes Cheilopogon unicolor, Cheilopogon arcticeps and Cheilopogon atrisignis (Teleostei: Exocoetidae). Mitochondrial Dna Part B Resour. 2016, 1, 134–135. [Google Scholar] [CrossRef]
  38. Li, Y.; Lu, J.; Wang, Z. Complete mitochondrial genome of Lasiopodomys mandarinus mandarinus (Arvicolinae, Rodentia). Mitochondrial Dna A Dna Mapp. Seq. Anal. 2016, 27, 1459–1460. [Google Scholar] [CrossRef] [PubMed]
  39. Zhang, D.X.; Hewitt, G.M. Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem. Syst. Ecol. 1997, 25, 99–120. [Google Scholar] [CrossRef]
  40. Fernã, N.-S.P.; Enriquez, J.A.; Montoya, J. Replication and transcription of mammalian mitochondrial DNA. Exp. Physiol. 2010, 88, 41–56. [Google Scholar] [CrossRef]
  41. Lee, W.J.; Conroy, J.; Howell, W.H.; Kocher, T.D. Structure and evolution of teleost mitochondrial control regions. J. Mol. Evol. 1995, 41, 54–66. [Google Scholar] [CrossRef]
  42. Li, A.; Munehara, H. Complete mitochondrial genome of the Antlered sculpin Enophrys diceraus (Scorpaeniformes, Cottidae). Mitochondrial Dna 2013, 26, 125–126. [Google Scholar] [CrossRef]
  43. Xu, T.J.; Cheng, Y.Z.; Liu, X.Z.; Shi, G.; Wang, R.X. The complete mitochondrial genome of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae): Genome characterization and phylogenetic considerations. Mol. Biol. 2011, 45, 392–403. [Google Scholar] [CrossRef]
  44. Wang, Q.; Wang, J.; Luo, J.; Chen, G. The complete mitochondrial genome of the Synanceia verrucosa (Scorpaeniformes: Synanceiidae). Mitochondrial Dna A Dna Mapp. Seq. Anal. 2016, 27, 4234–4235. [Google Scholar] [CrossRef] [PubMed]
Genes 10 00533 g001 550
Figure 1. Relative synonymous codon usage (RSCU) of the C. crocodilus mitogenome.
Figure 1. Relative synonymous codon usage (RSCU) of the C. crocodilus mitogenome.
Genes 10 00533 g001
Genes 10 00533 g002 550
Figure 2. The AT skew and GC skew in the PCGs of the C. crocodilus mitogenome.
Figure 2. The AT skew and GC skew in the PCGs of the C. crocodilus mitogenome.
Genes 10 00533 g002
Genes 10 00533 g003 550
Figure 3. Phylogenetic tree based on the nucleotide sequences of the 13 PCGs, using Bayesian inference (BI) and maximum likelihood (ML) methods, among Scorpaeniformes species. Perciformes fishes (P. novemaculeata (NC_024850.1) and S. knerii (NC_015987.1)) were used as outgroups. Numbers above the branches refer to Bayesian posterior probability values (blue numbers) and ML bootstrap values (red numbers).
Figure 3. Phylogenetic tree based on the nucleotide sequences of the 13 PCGs, using Bayesian inference (BI) and maximum likelihood (ML) methods, among Scorpaeniformes species. Perciformes fishes (P. novemaculeata (NC_024850.1) and S. knerii (NC_015987.1)) were used as outgroups. Numbers above the branches refer to Bayesian posterior probability values (blue numbers) and ML bootstrap values (red numbers).
Genes 10 00533 g003
Table
Table 1. Summary of the base composition and skewness of whole mitochondrial genomes and 13 protein-coding genes (PCGs) among 39 Scorpaeniformes species.
Table 1. Summary of the base composition and skewness of whole mitochondrial genomes and 13 protein-coding genes (PCGs) among 39 Scorpaeniformes species.
FamilySpeciesAccession NumberSize (bp)Whole Genome CompositionPCGs
A%G%T%C%A + T%AT skewGC skewAT skewGC skew
AnoplopomatidaeErilepis zoniferAP017436.116,500 26.7017.7929.0326.4853.180.0041 −0.2400 −0.0965 −0.2715
AnoplopomatidaeAnoplopoma fimbriaKP777542.116,50726.0718.3429.5826.0152.080.0010−0.2346−0.0994−0.2701
AgonidaeAspidophoroides olrikiiKT004434.117,20027.3016.9027.2828.4955.80−0.0214−0.2350−0.1153−0.2662
CottidaeIcelus spatulaKT004432.116,384 26.4317.4330.0426.0352.460.0078 −0.2656 −0.0987 −0.2935
CottidaeTrachidermus fasciatusJX017305.116,536 26.3318.1330.0725.4751.800.0167 −0.2478 −0.0815 −0.2750
CottidaeCottus poecilopusEU332750.116,560 25.6918.1830.3925.7451.43−0.0011 −0.2513 −0.1040 −0.2822
CottidaeCottus hangiongensisNC_014851.116,598 25.4818.2230.4025.8751.36−0.0075 −0.2506 −0.1135 −0.2840
CottidaeEnophrys dicerausNC_022147.116,976 27.5316.6528.6427.1954.710.0062 −0.2648 −0.0869 −0.2987
CottidaeCottus reiniiAP004442.116,56126.3017.6330.2825.7852.090.0100−0.2640−0.0919−0.2999
CottidaeMesocottus haitejNC_022181.116,527 26.6417.3529.8826.1252.760.0099 −0.2653 −0.0853 −0.2995
CottidaeCottus szanagaKX762050.116,518 26.5117.3929.9426.1652.670.0067 −0.2650 −0.0853 −0.3000
CottidaeCottus czerskiiNC_025242.116,53426.4017.4730.0726.0752.470.0063−0.2651−0.0850−0.3130
CottidaeCottus dzungaricusNC_024739.116,527 26.9317.0729.7026.3053.230.0119 −0.2701 −0.0822 −0.3135
HexagrammidaePleurogrammus monopterygiusNC_023475.116,575 27.0517.1528.6827.1254.17−0.0012 −0.2517 −0.0956 −0.2913
HexagrammidaeHexagrammos lagocephalusKP682334.116,505 26.9817.2629.4826.2953.270.0130 −0.2615 −0.0836 −0.3002
HexagrammidaeHexagrammos otakiiAP017443.116,513 26.9017.3329.8725.9052.800.0189 −0.2656 −0.0727 −0.3067
HexagrammidaeHexagrammos agrammusAB763992.116,514 26.8817.2329.7226.1553.030.0137 −0.2659 −0.0839 −0.3009
PeristediidaeSatyrichthys amiscusAP004441.116,526 27.0617.0028.5627.3854.44−0.0058 −0.2537 −0.1030 −0.2883
PlatycephalidaePlatycephalus indicusAP006783.116,542 27.3216.9930.1225.5752.880.0332 −0.2787 −0.0537 −0.3143
PlatycephalidaeCociella crocodilaMH521260.117,314 29.1415.9127.7627.1956.330.0348 −0.2713 −0.0375 −0.2911
ScorpaenidaePterois milesNC_024746.116,497 27.4018.2628.7425.6053.000.0340 −0.2229 −0.0509 −0.2420
ScorpaenidaeScorpaenopsis cirrosaNC_027735.116,966 27.9117.7128.0226.3554.270.0288 −0.2254 −0.0601 −0.2448
SebastidaeSebastes oblongusNC_024549.116,396 27.9116.9628.7226.4154.320.0276 −0.2574 −0.0595 −0.2842
SebastidaeSebastes thompsoniNC_027447.116,405 27.9817.0428.2126.7754.750.0222 −0.2468 −0.0677 −0.2781
SebastidaeSebastes minorNC_027444.116,408 27.7917.3327.6327.2555.040.0099 −0.2292 −0.0815 −0.2609
SebastidaeSebastes trivittatusNC_027446.116,409 27.8617.0928.3726.6754.540.0218 −0.2480 −0.0690 −0.2789
SebastidaeSebastes taczanowskiiKJ525744.116,452 27.7117.2928.5326.4754.180.0229 −0.2452 −0.0709 −0.2748
SebastidaeSebastes owstoniKJ834063.116,465 27.7117.3028.4126.5754.280.0210 −0.2430 −0.0748 −0.2723
SebastidaeSebates schlegeliAY491978.1 16,525 27.4717.4528.7426.3453.800.0210 −0.2444 −0.0680 −0.2742
SebastidaeSebastiscus marmoratusAP005989.117,301 28.6916.5128.1426.6755.360.0364 −0.2605 −0.0592 −0.2855
SebastidaeSebastes pachycephalusKF836442.116,415 27.7717.2428.6626.3354.100.0266 −0.2488 −0.0635 −0.2786
SebastidaeSebastes longispinisKJ834061.116,445 27.9117.1228.3126.6654.570.0230 −0.2462 −0.0699 −0.2713
SebastidaeSebastes steindachneriNC_027445.116,450 27.3617.5428.0127.0954.460.0049 −0.2298 −0.0873 −0.2587
SebastidaeSebastes hubbsiKJ525745.116,453 27.8617.2028.2826.6654.520.0221 −0.2436 −0.0671 −0.2714
SebastidaeSebastes vulpesNC_027438.116,462 27.7117.1428.6126.5554.250.0214 −0.2506 −0.0706 −0.2808
SebastidaeSebastes inermisKF725093.116,504 27.7717.1328.3426.7654.530.0186 −0.2466 −0.0714 −0.2735
SynanceiidaeSynanceia verrucosaNC_026989.116,506 31.0115.0625.6028.3459.350.0451 −0.2593 −0.0286 −0.2883
TriglidaeChelidonichthys kumuKY379222.116,495 26.6317.0431.1325.2051.830.0277 −0.2925 −0.0640 −0.3351
TriglidaeLepidotrigla micropteraKY012348.116,610 26.5317.1531.2825.0551.580.0286 −0.2918 −0.0666 −0.3326
Table
Table 2. Organization constituents of the mitochondrial genome of Cociella crocodilus.
Table 2. Organization constituents of the mitochondrial genome of Cociella crocodilus.
FeatureStrand *PositionSpacer (+)/Overlap (−)Start/Stop codon
tRNA-Phe (F)H1–670
12S rRNAH68–10880
tRNA-Val (V)H1089–1160−2
16S rRNAH1161–287379
tRNA-Leu (L1)H2872–29473
nad1H3027–4001−1ATG/TAA
tRNA-Ile (I)H4005–407423
tRNA-Gln (Q)L4074–41441
tRNA-Met (M)H4168–42376
nad2H4239–52880ATG/TAA
tRNA-Trp (W)H5295–53661
tRNA-Ala (A)L5367–543530
tRNA-Asn (N)L5437–55090
tRNA-Cys (C)L5540–56071
tRNA-Tyr (Y)L5608–56750
cox1H5677–72273GTG/TAA
tRNA-Ser (S1)L7228–72985
tRNA-Asp (D)H7302–73740
cox2H7380–80701ATG/T
tRNA-Lys (K)H8071–8145−10
atp 8H8147–8311−1ATG/TAA
atp 6H8302–8985−1ATG/TAA
cox3H8985–97700ATG/TAA
tRNA-Gly (G)H9770–9841−2
nad3H9842–10,1920ATG/TAG
tRNA-Arg (R)H10,191–10,259−7
nad4lH10,260–10,5560ATG/TAA
nad4H10,550–11,930−7ATG/T
tRNA-His (H)H11,931–11,9990
tRNA-Ser (S1)H12,002–12,0692
tRNA-Leu (L1)H12,074–12,1464
nad5H12,147–13,9850ATG/TAA
nad6L13,982–14,503−4ATG/TAG
tRNAGlu (E)L14,505–14,5731
cytbH14,576–15,7252ATG/T
tRNA-Thr (T)H15,726–15,7970
tRNA-Pro (P)L15,797–15,867−1
Control regionH15,868–17,3140
* H and L refer to the heavy and light strand, respectively.

Share and Cite

MDPI and ACS Style

Cui, L.; Cao, R.; Dong, Y.; Gao, X.; Cen, J.; Lu, S. The First Complete Mitochondrial Genome of the Flathead Cociella crocodilus (Scorpaeniformes: Platycephalidae) and the Phylogenetic Relationships within Scorpaeniformes Based on Whole Mitogenomes. Genes 2019, 10, 533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10070533

AMA Style

Cui L, Cao R, Dong Y, Gao X, Cen J, Lu S. The First Complete Mitochondrial Genome of the Flathead Cociella crocodilus (Scorpaeniformes: Platycephalidae) and the Phylogenetic Relationships within Scorpaeniformes Based on Whole Mitogenomes. Genes. 2019; 10(7):533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10070533

Chicago/Turabian Style

Cui, Lei, Rongbo Cao, Yuelei Dong, Xingchen Gao, Jingyi Cen, and Songhui Lu. 2019. "The First Complete Mitochondrial Genome of the Flathead Cociella crocodilus (Scorpaeniformes: Platycephalidae) and the Phylogenetic Relationships within Scorpaeniformes Based on Whole Mitogenomes" Genes 10, no. 7: 533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes10070533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop