Abstract
Gyrodactylus gurleyi, was inhabited on the fins and gills of goldfish (Carassius auratus), which belonged to the family Gyrodactylidae. In this study, we sequenced the complete mitochondrial genome of G. gurleyi with the total length of 14 771 bp. The mitogenome contained 12 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and two major non-coding regions (NC1 and NC2). The overall AT content was 72.1%. In phylogenetic analysis, G. gurleyi and G. kobayashii clustered together and then united with the clade of other three Gyrodactylus species (G. salaris, G. thymalli and G. derjavinoides) with high nodal support.
Gyrodactylus gurleyi was collected on the fins and gills of goldfish (Carassius auratus) from Wuhan (30°31′23′′N, 114°23′01′′E), China. It was identified by morphology and ITS molecular marker (Li et al. Citation2013). The specimen (accession no. IHB20150315006) was stored in the Museum of Aquatic Organisms, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China. The sequences of total mt genomic DNA of G. gurleyi were retrieved using long PCR and Sanger method of DNA sequencing.
The complete mt genome of G. gurleyi was circular and 14 771 bp in size (GenBank accession no. KU659806). It contained 12 protein-coding genes (PCGs, lacking Atp8), 22 tRNA genes, two rRNA genes and two major non-coding regions (NC1 and NC2) (). All the genes were transcribed from the same strand. The nucleotide composition was computed at 29.2%A, 16.9% G, 42.9% T and 11.1% C. The AT content of 72.1% was lower than that of Paragyrodactylus variegatus (76.3%), but higher than that of other four Gyrodactylus species (G. salaris, 62.3%; G. thymalli, 62.8%; G. derjavinoides, 68.2%; G. kobayashii, 71.6%) (Huyse et al. Citation2007; Plaisance et al. Citation2007; Huyse et al. Citation2008; Ye et al. Citation2014; Zhang et al. Citation2016).
Table 1. Organization of the mitochondrial genome of Gyrodactylus gurleyi.
The length of 12 PCGs was 9942 bp, accounting for 67.3% of the full length of the genome. As the four Gyrodactylus species (G. salaris, G. thymalli, G. derjavinoides and G. kobayashii), ATG was the unique start codon. The stop codon TAG was only found in three PCGs (Nad3, Nad5 and Nad6), whereas TAA in the rest of the PCGs. The size of the 22 tRNA genes was 1374 bp, varying from 58 bp (tRNASer(AGN)) to 71 bp (tRNAGlu). Nineteen of them had conventional secondary structure, while tRNASer(AGN), tRNASer(UCN) and tRNACys lacked DHU arms, which was similar to the other known Gyrodactylidae species. The rrnL and rrnS were 954 bp and 710 bp in size, respectively. They were adjacent to tRNAThr (upstream) and Cox2 (downstream), and separated by tRNACys, as described with other reported monopisthocotyleans (Huyse et al. Citation2007; Plaisance et al. Citation2007; Huyse et al. Citation2008; Perkins et al. Citation2010; Ye et al. Citation2014; Zhang et al. Citation2014a, Citation2014b, Citation2016). In addition, there were 37 bp overlapping sequences and 1763 space sequences, among which the biggest were NC1 (783 bp) and NC2 (783 bp).
Phylogenetic relationships between G. gurleyi and other nine monopisthocotyleans were inferred by using concatenated amino acid sequences of the 12 PCGs. The same tree topology was obtained by two different computational algorithms: Bayesian inference (BI) and maximum likelihood (ML), in which G. gurleyi and G. kobayashii clustered together and then united with the clade of 3 Gyrodactylus species (G. salaris, G. thymalli and G. derjavinoides) with high nodal support (). In addition, although G. gurleyi and G. kobayashii generally parasitized on the same host, sequence alignments showed that 99% (14 623 bp) mitogenome sequence of G. gurleyi was covered by that of G. kobayashii with only 80% identity. In contrast, G. salaris from Atlantic salmon and G. thymalli from grayling share 98% identity with 100% mitogenome alignment coverage.
Figure 1. Phylogenetic relationships between Gyrodactylus gurleyi and other 9 monopisthocotyleans based on 3044 concatenated amino acid sequences representing 12 mitochondrial protein-coding genes, with Polylabris halichoeres used as an outgroup. The MtZoa model for maximum-likelihood analysis and MtREV model for Bayes analysis are selected. Scale bar corresponds to the estimated number of substitutions per site. Bootstrap support values in percent units above nodes are displayed as follows: maximum likelihood bootstrap/Bayesian posterior probabilities.
Disclosure statement
The authors report no conflict of interest. The authors alone are responsible for the content and writing of the manuscript.
Funding information
This work was supported by the Earmarked Fund for China Agriculture Research System (CARS-46-08), the National Natural Science Foundation of China (31272695, 31572658) and the major scientific and technological innovation project of Hubei Province (2015ABA045).
Related Research Data
References
- Huyse T, Buchmann K, Littlewood DTJ. 2008. The mitochondrial genome of Gyrodactylus derjavinoides (Platyhelminthes: Monogenea) – a mitogenomic approach for Gyrodactylus species and strain identification. Gene. 417:27–34.
- Huyse T, Plaisance L, Webster BL, Mo TA, Bakke TA, Bachmann L, Littlewood DTJ. 2007. The mitochondrial genome of Gyrodactylus salaris (Platyhelminthes: Monogenea), a pathogen of Atlantic salmon (Salmo salar). Parasitology. 134:739–747.
- Li RR, Li WX, Wu XD, Wang GT. 2013. Identification of Gyrodactylus species in goldfish (Carassius auratus) through morphological study and the analysis of the rDNA ITS sequence. Acta Hydrobiol Sin. 38:903–909.
- Perkins EM, Donnellan SC, Bertozzi T, Whittington ID. 2010. Closing the mitochondrial circle on paraphyly of the Monogenea (Platyhelminthes) infers evolution in the diet of parasitic flatworms. Int J Parasitol. 40:1237–1245.
- Plaisance L, Huyse T, Littlewood DTJ, Bakke TA, Bachmann L. 2007. The complete mitochondrial DNA sequence of the monogenean Gyrodactylus thymalli (Platyhelminthes: Monogenea), a parasite of grayling (Thymallus thymallus). Mol Biochem Parasitol. 154:190–194.
- Ye F, King SD, Cone DK, You P. 2014. The mitochondrial genome of Paragyrodactylus variegatus (Platyhelminthes: Monogenea): differences in major non-coding region and gene order compared to Gyrodactylus. Parasit Vectors. 7:377.
- Zhang D, Zou H, Zhou S, Wu SG, Li WX, Wang GT. 2016. The complete mitochondrial genome of Gyrodactylus kobayashii (Platyhelminthes: Monogenea). Mitochondrial DNA Part B: Resoureces. 1:154–155.
- Zhang J, Wu X, Li Y, Xie M, Li A. 2014a. The complete mitochondrial genome of Tetrancistrum nebulosi (Monogenea: Ancyrocephalidae). Mitochondrial DNA. 1736:1–2.
- Zhang J, Wu X, Li Y, Zhao M, Xie M, Li A. 2014b. The complete mitochondrial genome of Neobenedenia melleni (Platyhelminthes: Monogenea): mitochondrial gene content, arrangement and composition compared with two Benedenia species. Mol Biol Rep. 41:6583–6589.