https://orcid.org/0000-0002-9433-5092
Abstract
The complete mitochondrial genome sequence of Omobranchus fasciolatoceps was firstly described in this article. The total length of mitogenome was 16,569 bp. It contains 13 protein-coding genes, 22 tRNA genes, and two ribosomal RNA genes. The overall base composition of H-strand was 29.04% A, 27.14% C, 27.89% T, and 15.93% G, with an A＋T bias of 56.93%. The phylogenetic analysis results showed that the O. fasciolatoceps was most closely related to O. elegans.
Omobranchus fasciolatoceps (Richardson 1846) belongs to genus Omobranchus in family Blenniidae of Blenniiformes (Kawaguchi et al. Citation1999). It is a native species distributed along southern China coast, and mainly appears in estuaries and the ocean shallows near shore. It commonly grows to an average length of 8 cm and can adapt to a wide range of temperatures. Its fleshy crest and zebra-striped face make this substrate dweller appealing.
The complete mitochondrial genome of O. fasciolatoceps was firstly reported in this paper. From that, we expected to provide useful information on the population genetics of O. fasciolatoceps and promote further molecular phylogenetic studies. The sample of O. fasciolatoceps in this article was collected from the south China coast (122°1′20″ E, 29°2′33″ N). The specimen was deposited at the Key Laboratory of Fisher Equipment and Engineering, Ministry of Agriculture (contacts, Liu Yimeng, [email protected]), with the voucher number Ofacsi20190710001. DNA libraries were prepared using Nextera DNA Flex Library Prep (Illumina, San Diego, CA, USA). High-throughput sequencing was performed on an Illumina HiSeq X-10 sequencer in 150-bp paired-end mode. Subsequently, the clean paired-end reads were assembled by the Perl script Novoplasty v 4.3.1(Dierckxsens et al. Citation2017).
The whole length of O. fasciolatoceps mitogenome was 16,569 bp. The nucleotide composition of the heavy strand was 29.04% A, 27.14% C, 27.89% T, and 15.93% G, with an A＋T bias of 56.93%. It contained 13 protein-coding genes (PCGs), 22 tRNAs and 2 rRNAs. Most genes were located on the heavy strand, except for ND6 and 8 tRNA genes (tRNAGln, RNAAla, tRNAAsn, tRNACys,tRNATyr, tRNASer, tRNAGlu, tRNAPro). The control region (913 bp) was identified between tRNAPro and tRNAPhe. Eleven PCGs shared the start codon ATG, except for COXI and ND6 (GTG as start codon). Six PCGs (ND1, COX1, ATPase 8, ND4L, ND6) were terminated with the typical stop codons (TAA), while the incomplete termination codons (T––/TA–) were identified among seven PCGs (ND2, COX2, ATPase 6, COX3, ND3, ND4, and Cyt b). In addition, ND5 ended with AGA.
For phylogenetic analysis, 15 taxa were selected: 11 species from order Blenniiformes (including Omobranchus fasciolatoceps), two species from order Cichliformes, and two species from order Lutjaniformes. Nucleotide sequences of the other mitochondrial genomes were downloaded from GenBank. The individual nucleotide sequences of 13 PCGs, two rRNAs genes, and 22 tRNAs genes were aligned by the ClustalW program (Larkin et al. Citation2007). A Maximum Likelihood phylogenetic tree (shown in ) was constructed using MEGA 6 (Stamatakis Citation2006; Tamura et al. Citation2013) with GTR + I + G model. It showed O. fasciolatoceps was situated near O. elegans, and had a closer relationship with Petroscirtes breviceps.
Figure 1. A phylogenetic tree for Omobranchus fasciolatoceps and other 14 species representing orders Blenniiformes, Cichliformes, and Lutjaniformes based on assembled nucleotide sequences of 13 protein-coding, two rRNAs genes, and 22 tRNAs genes. The number on each node indicates the values of ultrafast bootstrap (UFB) of 1000 replications.
Ethical approval
Experiments were performed in accordance with the recommendations of the Ethics Committee of Marine Fisheries Research Institute of Zhejiang.
These policies were enacted according to the Chinese Association for the Laboratory Animal Sciences and the Institutional Animal Care and Use Committee (IACUC) protocols.
Author contributions
Tao Zhang: Methodology, Data curation, Formal analysis, Writing- manuscript. Shiwei Wu: Methodology, Writing- manuscript. Lei Yang: Methodology, Software, Data curation. Hongling Ping: Methodology. Bin Lu: Investigation. Jun Liang: Investigation. Xuejun Yu: Investigation. Yanming Sui: Investigation. Yimeng Liu: Methodology, Investigation. Huilai Shi: Methodology, Funding acquisition. All authors agreed to the published version of the manuscript.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at [https://www.ncbi.nlm.nih.gov] under the accession no. OK501311. The associated **BioProject**, **SRA**, and **Bio-Sample** numbers are PRJNA770635, SRR16298008, and SAMN22225431 respectively.
Additional information
Funding
References
- Dierckxsens N, Mardulyn P, Smits G. 2017. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 45(4):e18.
- Kawaguchi T, Kohno H, Fujita K, Taki Y. 1999. Early morphological development of Omobranchus fasciolatoceps and O. Punctatus (Blenniidae: Omobranchini) reared in an aquarium. Ichthyological Research. 46(2):163–170.
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. 2007. Clustal W and Clustal X version 2.0. Bioinformatics. 23(21):2947–2948.
- Stamatakis A. 2006. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 22(21):2688–2690.
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 30(12):2725–2729.