Mitochondrial genome of Hypoderaeum conoideum – comparison with selected trematodes
© Yang et al.; licensee BioMed Central. 2015
Received: 24 November 2014
Accepted: 6 February 2015
Published: 12 February 2015
Hypoderaeum conoideum is a neglected but important trematode. The life cycle of this parasite is complex: snails serve as the first intermediate hosts: bivalves, fishes or tadpoles serve as the second intermediate hosts, and poultry (such as chickens and ducks) act as definitive hosts. In recent years, H. conoideum has caused significant economic losses to the poultry industry in some Asian countries. Despite its importance, little is known about the molecular ecology and population genetics of this parasite. Knowledge of mitochondrial (mt) genome of H. conoideum can provide a foundation for phylogenetic studies as well as epidemiological investigations.
The entire mt genome of H. conoideum was amplified in five overlapping fragments by PCR and sequenced, annotated and compared with mt genomes of selected trematodes. A phylogenetic analysis of concatenated mt amino acid sequence data for H. conoideum, eight other digeneans (Clonorchis sinensis, Fasciola gigantica, F. hepatica, Opisthorchis felineus, Schistosoma haematobium, S. japonicum, S. mekongi and S. spindale) and one tapeworm (Taenia solium; outgroup) was conducted to assess their relationships.
The complete mt genome of H. conoideum is 14,180 bp in length, and contains 12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and one non-coding region (NCR). The gene arrangement is the same as in Fasciola spp, with all genes being transcribed in the same direction. The phylogenetic analysis showed that H. conoideum had a relatively close relationship with F. hepatica and other members of the Fasciolidae, followed by the Opisthorchiidae, and then the Schistosomatidae.
The mt genome of H. conoideum should be useful as a resource for comparative mt genomic studies of trematodes and for DNA markers for systematic, population genetic and epidemiological studies of H. conoideum and congeners.
Echinostomatid trematodes comprise a group of at least 60 species , some of which are of socioeconomic significance in animals. Hypoderaeum conoideum (Bloch, 1782) is an important member of the family. This echinostomatid was originally found in the intestines of birds and is known to infect chickens, ducks and geese in many countries around the world [2-4]. It has also been found to infect humans and cause echinostomiasis in Thailand [5,6]. Freshwater snails, Planorbis corneus, Indoplanorbis exustus, Lymnaea stagnalis, L. limosa, L. ovata and L. rubiginosa, act as first intermediate hosts and shed the cercariae; bivalves, fishes or tadpoles can act as second intermediate hosts [3,5].
The accurate identification of species and genetic variants of Hypoderaeum conoideum will be central to investigating its biology, epidemiology and ecology, and also has implications for the diagnosis of infections. Although morphological features are used to identify this and other trematodes, such characters are not always reliable . Due to these constraints, various molecular methods have been established for specific identification . For instance, PCR-based techniques using genetic markers in nuclear ribosomal (r) and mitochondrial (mt) DNA have been widely used . The sequences of the first and second internal transcribed spacers (ITS-1 and ITS-2 = ITS) of nuclear rDNA have been particularly useful for specific identification, based on consistent levels of sequence difference between species and little variation within individual species , while the mitochondrial gene cox1 has been used for studying genetic variation and relationships among different species [8-10]. As a basis for the development of molecular tools to study H. conoideum populations (irrespective of developmental stage), we have characterized the complete mt genome of this parasite, compared this genome with those of selected trematodes and undertaken a phylogenetic analysis of concatenated amino acid sequence data for 12 protein-coding genes to assess the genetic relationship of H. conoideum with these other trematodes.
Parasites and DNA isolation
H. conoideum adults were collected from the intestine of a naturally infected free-range duck in Hubei province, China, in accordance with the Animal Ethics Procedures and Guidelines of Huazhong Agricultural University. These worms were washed in physiological saline and identified morphologically according to existing morphological descriptions . A reference specimen was stained and mounted  and the remaining specimens were fixed in 70% (v/v) ethanol and stored at −20°C until use . Total genomic DNA was extracted from one specimen using E.Z.N.A.® Tissue DNA Kit. To provide further identification for this specimen, the ITS-2 region was amplified and sequenced , it was identical to a reference sequence available for H. conoideum (GenBank accession no. KJ 944311.1).
Amplification and sequencing of partial cox1, cox3, nad4, nad5 and rrnS
Sequences of primers used to amplify fragments from Hypoderaeum conoideum
Sequences(5′ to 3′)
rrnS - nad52
rrnS - nad52
Long-PCR amplification and sequencing
Ten additional primers (see Table 1) were then designed from the sequences obtained, and used to amplify genomic DNA (~40-80 ng) from five regions (see Table 1) by long-PCR; PCRs (25 μl) were performed in a reaction buffer containing 2 mM MgCl2, 1× LA Taq Buffer II, 0.4 mM dNTP mixture, 0.8 μM of each primer, 2.5 U LA Taq polymerase (Takara) and 2.5 μl of genomic DNA or H2O (no-DNA control) for 35 cycles of 94°C/30 s (denaturation), 50°C/30 s (annealing) and 72°C/1 min (extension) per kb. Amplicons were cloned into pGEM-T-Easy vector (Promega, USA) according to the manufacturer’s protocol; inserts were amplified by long-range PCR (employing vector primers M13 and M14) and then sequenced using a primer-walking strategy .
Sequences were assembled using the software ContigExpress program (Invitrogen, Carlsbad, CA), and aligned against the mt genome sequences of other available trematodes (including F. hepatica) using the programs Clustal X v.1.83  to infer gene boundaries. The open reading frames (ORFs) were identified using ORF Finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html) employing the flatworm mitochondrial genetic code. Translation initiation and termination codons were identified as described previously [14,17,18]. The secondary structures of the 22 tRNA genes were predicted using tRNAscan-SE and/or manual adjustment [9,19]. The two rRNA genes were identified by comparison with those from the mt genome of F. hepatica . Amino acid sequences of the protein-coding genes were obtained by using the flatworm mt code, and aligned using the program MUSCLE  employing default settings.
Sliding window analysis of nucleotide variation
Sequence variability between H. conoideum and F. hepatica was conducted by sliding window analysis using the software DnaSP v.5 . A sliding window analyses was implemented as described previously .
Amino acid sequences conceptually translated from individual genes of the mt genome of H. conoideum were concatenated and aligned with those from available mt genomes of trematodes, including Clonorchis sinensis (NC_012147) [14,23], Fasciola gigantica (NC_024025) , F. hepatica (NC_002546) , Opisthorchis felineus (NC_011127) , Schistosoma haematobium (NC_008074) , Schistosoma japonicum (AF215860) , Schistosoma mekongi (NC_002529) , Schistosoma spindale (NC_008067) , and the cestode Taenia solium (outgroup) (NC_004022.1) . The phylogenetic analysis was conducted using the neighbour-joining (NJ) method employing the Tamura-Nei model . Confidence limits were assessed using bootstrap procedure with 1000 pseudo-replicates for neighbour-joining tree, and other settings were obtained using the default values in MEGA v.6.0 . In addition, maximum parsimony (MP), Bayesian (MB) and maximum likelihood (ML) analyses were implemented as described previously by other workers [20,26,27].
Features of the mt genome of H. conoideum
The organization of the mitochondrial genome of Hypoderaeum conoideum
Number of aa 1
Ini/Ter codons 2
Non coding region
Nucleotide contents of genes and the non-coding region within the mitochondrial genome of Hypoderaeum conoideum
A + T (%)
Non coding region
Codon usage for 12 protein-coding genes in the mitochondrial genome of Hypoderaeum conoideum
Transfer RNA and ribosomal RNA genes, and non-coding regions
The H. conoideum mt genome encodes 22 tRNAs; all of them have a typical cloverleaf structure. The length of 22 tRNA genes ranges from 60 bp to 75 bp (Table 2). There are intergenic and overlapping nucleotides between adjacent tRNA genes (Table 2). The rrnS and rrnL are 751 bp and 979 bp in length, respectively (Table 2). The location of rrnS is between tRNA-Cys and cox2, and that of rrnL is between tRNA-Thr and tRNA-Cys, which is the same as other trematodes. In contrast to some other trematodes (two AT-rich regions), such as F. hepatica and F. gigantica [14,23], O. felineus  and S. haematobium , there is only one AT-rich region (348 bp) in the mt genome of H. conoideum, which is located between tRNA-Glu and cox3 (Figure 1 and Table 2), with an A + T content of 60.19% (Table 3).
A comparison of nucleotide variability between H. conoideum and F. hepatica
The present characterization of the mt genome of H. conoideum provides a basis for addressing questions regarding the biology, epidemiology and population genetics of Hypoderaeum spp. In addition, it will also assist in supporting taxonomic studies of Hypoderaeum spp. of other animals (e.g., chickens, ducks, geese and humans) as well as in tracking life cycles by identifying larval stages in different intermediate hosts using molecular tools.
Assisted by sliding window analysis, PCR primers could be selectively designed to regions conserved among different trematode species and flanking variable regions in the mt genome that are informative (based on sequencing from a small number of individuals from particular populations). PCR-coupled single-strand conformation polymorphism (SSCP) analysis  could then be employed to screen large numbers of individuals representing different populations and, based on such an analysis, samples representing all detectable genetic variability could be selected for subsequent sequencing and analyses. Such an approach has been applied to study the genetic make-up of the blood fluke S. japonicum from seven provinces in China [30,31].
Now that the H. conoideum mt genome is available, it would be interesting to undertake a comprehensive study of this morphospecies from various host species from different countries by integrating morphological data with PCR-based genetic analyses of adult worms and larval stages (from intermediate hosts) to begin to understand the epidemiology and ecology of H. conoideum. In addition to conducting targeted mt genetic analyses, it would also be useful to include analyses of sequence variability in the two internal transcribed spacers (ITS-1 and ITS-2), 18S and 28S of nuclear ribosomal DNA, because, for trematodes, these markers usually allow specific identification of trematodes. Importantly, although H. conoideum is recognized as a species, it is possible that cryptic species of this taxon might exist. This proposal could be tested using the mt markers defined here, together with ITS-1 and/or ITS-2.
Our analysis showed that H. conoideum is genetically closely related to F. hepatica comparing with other trematodes. The mt genome of H. conoideum should be useful as a resource for comparative mt genomic studies of trematodes and DNA markers for systematic, population genetic and epidemiological studies of H. conoideum and congeners.
Sincere thanks to Professor Ian Beveridge for comments on the manuscript and assistance with the staining of trematodes. Thanks to Ross Hall and Namitha Mohandas for bioinformatic assistance. This work was supported by the “Special Fund for Agro-scientific Research in the Public Interest” (Grant no. 201303037), “Fundamental Research Funds for the Central Universities” (Program. 2013PY059) and the Australian Research Council (RBG).
- Sorensen RE, Curtis J, Minchella DJ. Intraspecific variation in the rDNA its loci of 37-collar-spined echinostomes from North America: implications for sequence-based diagnoses and phylogenetics. J Parasitol. 1998;84(5):992–7.View ArticlePubMedGoogle Scholar
- Steele JH. CRC Handbook Series in Zoonoses. Section C: Parasitic Zoonoses (Vol. III). Boca Raton, Florida, USA: CRC Press; 1982.Google Scholar
- Yamaguti S. Systema helminthum (Vol. I): part I. The digenetic trematodes of vertebrates. New York, USA: Interscience Publishers Inc.; 1958.Google Scholar
- Rim HJ. Echinostomiasis. Boca Raton, Florida, USA: CRC Press; 1982.Google Scholar
- Harinasuta T, Bunnag D, Radomyos P. Intestinal fluke infections. In: Bailliere's clinical tropical medicine and communicable diseases, vol. 2. 1987. p. 695–721.Google Scholar
- Yokogawa M, Harinasuta C, Charoenlarp P. Hypoderaeum conoideum (Block, 1782) Diez, 1909, a common intestinal fluke of man in the north-east Thailand. Jpn J Parasitol. 1965;14:148–53.Google Scholar
- Nolan MJ, Cribb TH. The use and implications of ribosomal DNA sequencing for the discrimination of digenean species. Adv Parasitol. 2005;60:101–63.View ArticlePubMedGoogle Scholar
- Liu GH, Wang Y, Song HQ, Li MW, Ai L, Yu XL, et al. Characterization of the complete mitochondrial genome of Spirocerca lupi: sequence, gene organization and phylogenetic implications. Parasit Vectors. 2013;6:45.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu M, Chilton NB, Gasser RB. The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea). Int J Parasitol. 2002;32(2):145–58.View ArticlePubMedGoogle Scholar
- Saijuntha W, Sithithaworn P, Duenngai K, Kiatsopit N, Andrews RH, Petney TN. Genetic variation and relationships of four species of medically important echinostomes (Trematoda: Echinostomatidae) in South-East Asia. Infect Genet Evol. 2011;11(2):375–81.View ArticlePubMedGoogle Scholar
- Chai JY, Shin EH, Lee SH, Rim HJ. Foodborne intestinal flukes in Southeast Asia. Korean J Parasitol. 2009;47(Suppl):S69–102.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang GL. Laboratory diagnostic techniques of parasites. Feeding Livestock. 2013;3:39–43.Google Scholar
- Tantrawatpan C, Saijuntha W, Sithithaworn P, Andrews RH, Petney TN. Genetic differentiation of Artyfechinostomum malayanum and A. sufrartyfex (Trematoda: Echinostomatidae) based on internal transcribed spacer sequences. Parasitol Res. 2013;112(1):437–41.View ArticlePubMedGoogle Scholar
- Le TH, Blair D, McManus DP. Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda). Parasitology. 2001;123(6):609–21.View ArticlePubMedGoogle Scholar
- Hu M, Jex AR, Campbell BE, Gasser RB. Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing. Nat Protoc. 2007;2(10):2339–44.View ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Ai L, Weng YB, Elsheikha HM, Zhao GH, Alasaad S, Chen JX, et al. Genetic diversity and relatedness of Fasciola spp. isolates from different hosts and geographic regions revealed by analysis of mitochondrial DNA sequences. Vet Parasitol. 2011;181(2–4):329–34.View ArticlePubMedGoogle Scholar
- Le TH, Blair D, Agatsuma T, Humair PF, Campbell NJ, Iwagami M, et al. Phylogenies inferred from mitochondrial gene orders-a cautionary tale from the parasitic flatworms. Mol Biol Evol. 2000;17(7):1123–5.View ArticlePubMedGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.View ArticlePubMed CentralPubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Liu GH, Gasser RB, Young ND, Song HQ, Ai L, Zhu XQ. Complete mitochondrial genomes of the ‘intermediate form’ of Fasciola and Fasciola gigantica, and their comparison with F. hepatica. Parasit Vectors. 2014;7:150.View ArticlePubMed CentralPubMedGoogle Scholar
- Shekhovtsov SV, Katokhin AV, Kolchanov NA, Mordvinov VA. The complete mitochondrial genomes of the liver flukes Opisthorchis felineus and Clonorchis sinensis (Trematoda). Parasitol Int. 2010;59(1):100–3.View ArticlePubMedGoogle Scholar
- Littlewood DT, Lockyer AE, Webster BL, Johnston DA, Le TH. The complete mitochondrial genomes of Schistosoma haematobium and Schistosoma spindale and the evolutionary history of mitochondrial genome changes among parasitic flatworms. Mol Phylogenet Evol. 2006;39(2):452–67.View ArticlePubMedGoogle Scholar
- Nakao M, Sako Y, Ito A. The mitochondrial genome of the tapeworm Taenia solium: a finding of the abbreviated stop codon U. J Parasitol. 2003;89(3):633–5.View ArticlePubMedGoogle Scholar
- Yin F, Gasser RB, Li F, Bao M, Huang W, Zou F, et al. Genetic variability within and among Haemonchus contortus isolates from goats and sheep in China. Parasit Vectors. 2013;6(1):279.View ArticlePubMed CentralPubMedGoogle Scholar
- Mohandas N, Pozio E, La Rosa G, Korhonen PK, Young ND, Koehler AV, et al. Mitochondrial genomes of Trichinella species and genotypes-a basis for diagnosis, and systematic and epidemiological explorations. Int J Parasitol. 2014;44(14):1073–80.View ArticlePubMedGoogle Scholar
- Yan HB, Wang XY, Lou ZZ, Li L, Blair D, Yin H, et al. The mitochondrial genome of Paramphistomum cervi (Digenea), the first representative for the family Paramphistomidae. PLoS One. 2013;8(8):e71300.View ArticlePubMed CentralPubMedGoogle Scholar
- Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AJ, Otranto D, et al. Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc. 2006;1(6):3121–8.View ArticlePubMedGoogle Scholar
- Chilton NB, Bao-Zhen Q, Bøgh HO, Nansen P. An electrophoretic comparison of Schistosoma japonicum (Trematoda) from different provinces in the People’s Republic of China suggests the existence of cryptic species. Parasitology. 1999;119(Pt 4):375–83.View ArticlePubMedGoogle Scholar
- Zhu XQ, Bøgh HO, Gasser RB. Dideoxy fingerprinting of low-level nucleotide variation in mitochondrial DNA of the human blood fluke, Schistosoma japonicum. Electrophoresis. 1999;20(14):2830–3.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0) which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.