Research Article |
Corresponding author: Rui Song ( ryain1983@163.com ) Academic editor: David Gibson
© 2020 Yuan-An Wu, Jin-Wei Gao, Xiao-Fei Cheng, Min Xie, Xi-Ping Yuan, Dong Liu, Rui Song.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Wu Y-A, Gao J-W, Cheng X-F, Xie M, Yuan X-P, Liu D, Song R (2020) Characterization and comparative analysis of the complete mitochondrial genome of Azygia hwangtsiyui Tsin, 1933 (Digenea), the first for a member of the family Azygiidae. ZooKeys 945: 1-16. https://doi.org/10.3897/zookeys.945.49681
|
Azygia hwangtsiyui (Trematoda, Azygiidae), a neglected parasite of predatory fishes, is little-known in terms of its molecular epidemiology, population ecology and phylogenetic study. In the present study, the complete mitochondrial genome of A. hwangtsiyui was sequenced and characterized: it is a 13,973 bp circular DNA molecule and encodes 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as two non-coding regions. The A+T content of the A. hwangtsiyui mitogenome is 59.6% and displays a remarkable bias in nucleotide composition with a negative AT skew (–0.437) and a positive GC skew (0.408). Phylogenetic analysis based on concatenated amino acid sequences of twelve protein-coding genes reveals that A. hwangtsiyui is placed in a separate clade, suggesting that it has no close relationship with any other trematode family. This is the first characterization of the A. hwangtsiyui mitogenome, and the first reported mitogenome of the family Azygiidae. These novel datasets of the A. hwangtsiyui mt genome represent a meaningful resource for the development of mitochondrial markers for the identification, diagnostics, taxonomy, homology and phylogenetic relationships of trematodes.
gene arrangement, phylogenetic analysis, Trematoda
The genus Azygia Looss, 1899 is an endoparasitic helminth found in the stomach and intestine of freshwater feral carnivorous fish (
Azygia hwangtsiyui Tsin, 1933 is a member of the family Azygiidae Odhner, 1911 and is often overlooked; it is parasitic in the gastrointestinal tract of species of the family Channidae Fowler, 1934 but caused only slight clinical signs, including malnutrition and weight loss. In China, Azygia hwangtsiyui-infected freshwater predatory fishes have been described from Shandong, Heilongjiang, Jiangsu, Fujian, Sichuan and Hunan Provinces (
Morphology is the most commonly used method for species identification and differentiation of metazoans and is widely adopted globally by parasitologists and taxonomists. A huge disadvantage of using morphological criteria, however, is that it is difficult to identify and distinguish closely related and cryptic species. Although the family Azygiidae was erected more than a century ago, its situation, and that of several species of Azygia, is still controversial and uncertain.
Mitochondrial (mt) genome and nuclear ribosomal DNA sequences are effective molecular tools for taxonomic identification, phylogeny and biogeographical research (
The specimens of flatworms were isolated from the stomach of their definitive host, in this case snakehead fish (Ophiocephalus argus (Cantor, 1842)) obtained from east Dongting Lake in Yueyang, Hunan province, China (29°22'N, 113°06'E). Azygia hwangtsiyui was morphologically identified according to the original and other descriptions (
According to conserved regions of mitochondrial genes in other available digenea mitogenomes, six partial gene fragments for cytb, nad4, nad1, 16S, 12S and cox2 were amplified using six generic primers sets HWF1/HWR1 (for cytb), HWF3/HWR3 (for nad4), HWF5/HWR5 (for nad1), HWF7/HWR7 (for 16S), HWF9/HWR9 (for 12S), and HWF11/HWR11 (for cox2), respectively. On the basis of these obtained nucleotide sequences, A. hwangtsiyui-specific primers were designed for amplification and sequencing of the remaining mitogenome (Suppl. material
According to sequence chromatograms, all raw fragments were quality-proofed using CHROMAS (https://www.technelysium.com.au) to remove ambiguity codes and low-quality bases. Whenever the quality was sub-optimal, sequencing was repeated until the amplicon is the consensus sequence. Before manual assembly of the entire mitochondrial genomic sequence, identification of all amplicons was performed by BLASTN check (
For phylogenetic analyses, we utilized translated and concatenated amino acid sequences of twelve protein-coding genes for 49 Platyhelminthes including A. hwangtsiyui mitogenome determined in this study. Two tapeworm species, Cloacotaenia megalops (Nitzsch in Creplin, 1829) (NC_032295.1) and Dibothriocephalus latus (Linnaeus, 1758) (NC_008945.1) were included as outgroup taxa representing two different families. Species information including systematic positions and GenBank accession numbers is provided in Suppl. material
The entire A. hwangtsiyui mtDNA is 13,973 bp in length (GenBank accession number: MN844889) and comprised of 12 protein-coding genes (cox1-3, nad1-6, nad4L, cytb, and atp6), 22 tRNA genes, two rRNA genes (rrnL and rrnS), and two non-coding regions. The 12 protein-coding gene order arrangement is cox3-cytb-nad4L-nad4-atp6-nad2-nad1-nad3-cox1-cox2-nad6-nad5 (Fig.
Gene | Position | Size | Intergenic nucleotides | Codon | Anti-codon | Strand | ||
---|---|---|---|---|---|---|---|---|
From | To | Start | Stop | |||||
cox3 | 1 | 660 | 660 | – | ATG | TAG | – | H |
trnH | 666 | 729 | 64 | +5 | – | – | GTG | H |
cytb | 732 | 1841 | 1110 | +2 | ATG | TAG | – | H |
nad4L | 1848 | 2108 | 261 | +6 | ATG | TAG | – | H |
nad4 | 2069 | 3340 | 1272 | –40 | ATG | TAG | – | H |
trnQ | 3345 | 3409 | 65 | +4 | – | – | TTG | H |
trnF | 3423 | 3488 | 66 | +13 | – | – | GAA | H |
trnM | 3490 | 3555 | 66 | +1 | – | – | CAT | H |
atp6 | 3556 | 4068 | 513 | – | ATG | TAG | – | H |
nad2 | 4072 | 4932 | 861 | +3 | GTG | TAG | – | H |
trnV | 4946 | 5009 | 64 | +13 | – | – | TAC | H |
trnA | 5013 | 5076 | 64 | +3 | – | – | TGC | H |
trnD | 5081 | 5146 | 66 | +4 | – | – | GTC | H |
nad1 | 5149 | 6054 | 906 | +2 | GTG | TAG | – | H |
trnN | 6070 | 6134 | 65 | +15 | – | – | GTT | H |
trnP | 6148 | 6212 | 65 | +13 | – | – | TGG | H |
trnI | 6216 | 6279 | 64 | +3 | – | – | GAT | H |
trnK | 6280 | 6348 | 69 | – | – | – | CTT | H |
nad3 | 6349 | 6708 | 360 | – | ATG | TAA | – | H |
trnS1 | 6712 | 6770 | 59 | +3 | – | – | GCT | H |
trnW | 6781 | 6842 | 62 | +10 | – | – | TCA | H |
cox1 | 6843 | 8396 | 1554 | – | TTG | TAG | – | H |
trnT | 8410 | 8474 | 65 | +13 | – | – | TGT | H |
rrnL | 8475 | 9449 | 975 | – | – | – | – | H |
trnC | 9450 | 9506 | 57 | – | – | – | GCA | H |
rrnS | 9507 | 10246 | 740 | – | – | – | H | |
cox2 | 10247 | 10828 | 582 | – | GTG | TAA | – | H |
nad6 | 10834 | 11277 | 444 | +5 | GTG | TAG | – | H |
trnY | 11284 | 11352 | 69 | +6 | – | – | GTA | H |
trnL1 | 11352 | 11416 | 65 | –1 | – | – | TAG | H |
trnS2 | 11421 | 11490 | 70 | +4 | – | – | TGA | H |
trnL2 | 11491 | 11555 | 65 | – | – | – | TAA | H |
trnR | 11558 | 11617 | 60 | +2 | – | – | TCG | H |
nad5 | 11626 | 13225 | 1600 | +8 | GTG | T | – | H |
trnE | 13226 | 13288 | 63 | – | – | – | TTC | H |
trnG | 13604 | 13669 | 66 | – | – | – | TCC | H |
Nucleotide contents of genes and the non–coding region within the mitochondrial genome of Azygia hwangtsiyui.
Regions | Size (bp) | T | C | A | G | AT (%) | GC (%) | AT skew | GC skew |
---|---|---|---|---|---|---|---|---|---|
atp6 | 513 | 50.5 | 11.3 | 12.9 | 25.3 | 63.4 | 36.6 | –0.594 | 0.383 |
cox1 | 1554 | 44.0 | 13.0 | 15.9 | 27.1 | 59.9 | 40.1 | –0.469 | 0.352 |
cox2 | 582 | 39.7 | 13.1 | 18.0 | 29.2 | 57.7 | 42.3 | –0.375 | 0.382 |
cox3 | 660 | 44.8 | 11.4 | 15.9 | 27.9 | 60.7 | 39.3 | –0.476 | 0.421 |
cytb | 1110 | 44.4 | 13.5 | 15.8 | 26.3 | 60.2 | 39.8 | –0.476 | 0.321 |
nad1 | 906 | 44.0 | 10.0 | 16.3 | 29.6 | 60.3 | 39.6 | –0.459 | 0.493 |
nad2 | 861 | 46.2 | 11.3 | 11.8 | 30.7 | 58.0 | 42.0 | –0.592 | 0.463 |
nad3 | 360 | 45.8 | 8.1 | 14.2 | 31.9 | 60.0 | 40.0 | –0.528 | 0.597 |
nad4 | 1272 | 45.7 | 12.8 | 12.3 | 29.2 | 58.0 | 42.0 | –0.577 | 0.391 |
nad4L | 264 | 47.0 | 8.0 | 17.8 | 27.3 | 64.8 | 35.3 | –0.450 | 0.548 |
nad5 | 1600 | 46.4 | 9.5 | 14.1 | 30.0 | 60.5 | 39.5 | –0.535 | 0.519 |
nad6 | 444 | 45.9 | 11.9 | 14.2 | 27.9 | 60.1 | 39.8 | –0.528 | 0.401 |
rrnL | 975 | 37.4 | 14.1 | 22.5 | 26.1 | 59.9 | 40.2 | –0.250 | 0.299 |
rrnS | 740 | 36.6 | 14.5 | 21.6 | 27.3 | 58.2 | 41.8 | –0.258 | 0.307 |
tRNAs | 1349 | 35.3 | 14.5 | 21.7 | 28.5 | 57.0 | 43.0 | –0.238 | 0.328 |
rRNAs | 1715 | 37.1 | 14.2 | 22.1 | 26.6 | 59.2 | 40.8 | –0.253 | 0.303 |
PCGs | 10126 | 45.2 | 11.5 | 14.7 | 28.6 | 59.9 | 40.1 | –0.509 | 0.425 |
Genome | 13973 | 42.8 | 12.0 | 16.8 | 28.5 | 59.6 | 40.5 | –0.437 | 0.408 |
A total of 3364 amino acids was encoded by the A. hwangtsiyui mtDNA. The full scale of 12 concatenated protein-coding genes was 10126 bp, composed of 45.2% T, 11.5% C, 14.7% A, and 28.6% G. Average A+T content of concatenated 12 protein-coding genes was 59.9%, varying from 57.7% (cox2) to 64.8% (nad4L) (Table
For the A. hwangtsiyui mitogenome, codon ends in G or T were more continual than those ending in A or C. The most frequently used start codon in protein-coding genes was ATG (for six PCGs), secondly was GTG (for five PCGs), which resembles that of the most frequent extrapolated start codons for mitogenome protein-encoding genes of digenean species (
Comparative analysis of gene arrangement among 47 selected digenean taxa, two gene blocks (cox1-trnT-rrnL-trnC-rrnS-cox2-nad6 and cytb-nad4L-nad4-trnQ) are shared by all selected taxa (Fig.
Phylogenetic relationships and gene arrangement of Azygia hwangtsiyui with other selected digeneas based on translated mitochondrial proteins. The concatenated amino-acid sequence datasets of the 12 protein-coding genes were analyzed by Bayesian Inference (BI) and Maximum Likelihood (ML), utilizing Cloacotaenia megalops (NC_032295.1) and Dibothriocephalus latus (NC_008945.1) as the outgroups. Both ML and BI analyses constructed identical tree topologies.
To assess phylogenetic relationships among available flatworms, we utilized concatenated amino acid sequence dataset representing 12 protein-coding genes of A. hwangtsiyui, 46 other digenean representatives, and two tapeworm species (C. megalops and D. latus) for analyzing molecular-based phylogeny. In this study, the topological structure is divided into two large clades: one consists of seven members of the family Schistosomatidae; and the other clade comprises 40 members from 16 families including the family Azygiidae (A. hwangtsiyui) (Fig.
This work was supported by the Earmarked Fund for China Agriculture Research System (CARS-45) and Important Research Project of Hunan Provincial Science and Technology Department (2016NK2176).
Table S1
Data type: molecular data
Explanation note: Primers for amplification and sequencing mitochondrial genome of Azygia hwangtsiyui.
Table S2
Data type: molecular data
Explanation note: Information of the Digenea and the outgroups for which complete mitogenomes are available in GenBank.
Table S3
Data type: molecular data
Explanation note: A+T content (%) for 12 protein-coding genes in the available 49 Platyhelminthes mitogenomes.
Figure S1
Data type: molecular data
Explanation note: Secondary structure of tandem repeats in non-coding region 2 (NCR2) of Azygia hwangtsiyui mitochondrial genome.
Figure S2
Data type: molecular data
Explanation note: Relative synonymous codon usage (RSCU) of Azygia hwangtsiyui mitochondrial genome.