The complete mitochondrial genome of the gullet worm Gongylonema pulchrum: gene content, arrangement, composition and phylogenetic implications
- Guo-Hua Liu†1,
- Yan-Qing Jia†2,
- Ya-Nan Wang2,
- Guang-Hui Zhao2Email author and
- Xing-Quan Zhu1Email author
© Liu et al.; licensee BioMed Central. 2015
Received: 28 October 2014
Accepted: 26 January 2015
Published: 15 February 2015
Gongylonema pulchrum (Nematoda: Gongylonematidae), a thread-like spirurid gullet worm, infects a range of mammalian definitive hosts, including cattle, pigs, equines, goats, primates and humans, and can cause gongylonemiasis.
In the present study, the complete mitochondrial (mt) genome of G. pulchrum was obtained using Long-range PCR and subsequent primer walking. The phylogenetic position of G. pulchrum within the Spiruromorpha was established using Bayesian analyses of the protein-coding genes at the amino acid level.
The length of this AT-rich (75.94%) mt genome is 13,798 bp. It contains 12 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and one non-coding region. The gene arrangement is the same as those of Thelazia callipaeda (Thelaziidae) and Setaria digitata (Onchocercidae), but distinct from that of Heliconema longissimum (Physalopteridae). Phylogenetic analyses, based on the concatenated amino acid sequence data for all 12 protein-coding genes using Bayesian inference (BI) method, showed that G. pulchrum (Gongylonematidae) was more closely related to Spirocerca lupi (Spiruroidea) than other members of the infraorder Spiruromorpha.
The present study represents the first mt genome sequence for the family Gongylonematidae, which provides the opportunity to develop novel genetic markers for studies of epidemiology, population genetics and systematics of this nematode of human and animal health significance.
Gongylonema pulchrum Molin, 1857, known as the ‘gullet worm’ because of its location in the upper digestive system of its definitive hosts, has a worldwide distribution and is responsible for gongylonemiasis of humans, sometimes causing serious complaints, e.g., expectoration of blood, numbness of tongue, vomiting, pharyngitis and stomatitis . G. pulchrum also infects domestic and wild animal hosts, including bears, camels, cattle, cervids, donkeys, equines, goats, non-human primates, pigs, and sheep, causing major economic losses . Fortunately, gongylonemiasis can be treated effectively using anthelmintics, such as levamisole, mebendazole and ivermectin .
Human gongylonemiasis, a neglected parasitic disease, has been reported from many countries (e.g., Australia, Bulgaria, Ceylon, China, France, Germany, Hungary, Iran, Japan, Laos, Morocco, New Zealand, Soviet Union, Spain, Sri Lanka, Turkey and USA). Recently, human gongylonemiasis has frequently been reported from China , France [5,6] and USA . Clearly, the increased frequency and more widespread occurrence of clinical cases of human gongylonemiasis mean that knowledge on enhanced diagnosis, treatment and control is needed. Although molecular markers in portion of mitochondrial (mt) DNA and the internal transcribed spacer (ITS) regions of nuclear ribosomal DNA (rDNA), have found utility for taxonomic and epidemiological studies of G. pulchrum , there is still a paucity of information on G. pulchrum in different hosts and countries around the world.
Due to its maternal inheritance, fast rate of evolutionary change, lack of recombination and relatively conserved genome structures , mt genomes have been widely used as genetic markers for population genetic structure and the study of phylogenetic relationships among nematodes and trematodes . In particular, concatenated amino acid sequences, derived from the protein-coding genes, lend themselves for assessing systematic relationships of parasitic nematodes [11-19]. The objectives of the present study were to characterize the mt genome of G. pulchrum, the first representative of the family Gongylonematidae, and to assess the phylogenetic position of this zoonotic nematode in relation to other nematodes of the infraorder Spiruromorpha.
The performance of this study was strictly according to the recommendations of the Guide for the Care and Use of Laboratory Animals of the Ministry of Health, China, and our protocol was reviewed and approved by the Research Ethics Committee of Northwest A&F University.
Adult specimens of G. pulchrum were collected from the oesophagus of a naturally infected goat in Shenmu, Shaanxi province, China, with no specific permits being required by the authority for the sample collection.
Genomic DNA extraction
The gullet worms were washed extensively in physiological saline, identified morphologically to species according to existing keys and descriptions , fixed in ethanol and then stored at −20°C until use. The mid-body section of each worm was used for the isolation of total genomic DNA using proteinase K digestion and mini-column purification (TIANamp Genomic DNA Purification System, TIANGEN). The identity of the specimen was verified by sequencing regions of the ITS-1 and ITS-2 rDNA using an established method ; the two regions were 99.7% and 100% identical to previously published sequences for G. pulchrum from Bos taurus in Iran (GenBank accession nos. AB495392 and AB513721, respectively).
Long-PCR and sequencing
Primers used to amplify the mitochondrial genome of Gongylonema pulchrum in three overlapping long PCR fragments
Name of primer
Sequence (5′ to 3′)
Sequences were assembled manually and aligned against the complete mt genome sequences of Spirocerca lupi  using the computer program MAFFT 7.122  to identify gene boundaries. Each gene was translated into its amino acid sequence using the invertebrate mitochondrial genetic code in MEGA 5 . The translation initiation and termination codons were identified to avoid gene overlap and to optimize the similarity between the gene lengths of closely related species of the infraorder Spiruromorpha. The program tRNAscan-SE  was used to find tRNA and infer their secondary structure, putative secondary structures of 19 tRNA genes were identified, and the remaining three tRNA genes (tRNA-Arg, tRNA-Ser AGN and tRNA-Ser UCN) were inferred by recognizing potential secondary structures and anticodon sequences by eye . Two rRNA genes were predicted by comparison with those of closely related nematodes of the infraorder Spiruromorpha [15,22].
Amino acid sequences inferred from published mt genomes representing 12 species of the infraorder Spiruromorpha (Brugia malayi: GenBank accession no. NC_004298; Wuchereria bancrofti: JN367461; Chandlerella quiscali: NC_014486; Loa loa: NC_016199; Acanthocheilonema viteae: NC_016197; Onchocerca flexuosa: NC_016172; Onchocerca volvulus: AF015193; Dirofilaria immitis: NC_005305; Setaria digitata: NC_014282; S. lupi: KC305876; Thelazia callipaeda: JX069968; Heliconema longissimum: NC_016127) were included in the present analysis, using Toxascaris leonina (NC_023504) as the outgroup . 12 amino acid sequences were separately aligned using MAFFT 7.122 and then concatenated, with ambiguously aligned regions excluded using Gblocks 0.91b (doc)  with the default parameters using the options for a less stringent selection. Phylogenetic analyses were conducted using Bayesian inference (BI) method. The JTT+G+F model of amino acid evolution was selected as the most suitable model of evolution by ProtTest 2.4  based on the Akaike information criterion (AIC). As the JTT model is not implemented in the current version of MrBayes, an alternative model, CpREV, was used in BI and four chains (three heated and one cold) were run simultaneously for the Monte Carlo Markov Chain. Two independent runs for 1,000,000 metropolis-coupled MCMC generations were used, sampling a tree every 100 generation in MrBayes 3.1.1 . At the end of each run, the average standard deviation of split frequencies was less than 0.01. In addition, the potential scale reduction factor (~1) was examined to ensure that the convergence had been achieved. A 50% majority rule consensus tree was obtained from BI. Of 10,000 trees, the first 2,500 trees represented burn-in and the remaining trees were used to calculate Bayesian posterior probabilities (Bpp). Phylograms were drawn using the program FigTree v.1.4 .
Results and discussion
Mitochondrial genome organization of Gongylonema pulchrum
Number of aa a
Comparison of A+T content (%) of gene and region of the mt genomes of spirurid nematodes sequenced to date (alphabetical order), including Gongylonema pulchrum (in bold)
Nucleotide composition of the mt genomes of spirurid nematodes, including that of Gongylonema pulchrum (in bold)
Nucleotide frequency (%)
Whole genome sequence
The most common initiation codon for G. pulchrum is TGG, followed by ATG, ATT, GTT and GTG (4, 3, 2, 2, 1 genes, respectively; Table 2). The most frequent complete termination codon is TAA (4 genes); nad4 is terminated with the codon TAG. Of the remaining genes, nad1, nad3, nad5, atp6 and cox2 are terminated with the abbreviated stop codon T, whereas cytb and nad2 are terminated with the abbreviated stop codon TA. This is consistent with the arrangement in the mt genomes of other nematodes [27,33-35].
Twenty-two tRNA genes were predicted from the mt genome of G. pulchrum and varied from 52 to 59 bp in length. Twenty of the 22 tRNA genes (excluding two tRNA-Ser) have a predicted secondary structure with a 3–5 bp DHU arm and a DHU loop of 7–9 bases, in which the variable TψC arm and loop are replaced by a “TV-replacement loop” of 8–10 bases. As seen in almost all other nematode mtDNAs , the tRNA-Ser gene of G. pulchrum mt genome is equipped with a TψC arm and loop but lacks the DHU arm and loop, consisting of a 6–8 bp TψC arm, TψC loop of 4–6 bases and a variable loop of 4 bases. The majority of nematode mtDNA sequences usually contain two non-coding regions with significant size difference , but there is only one non-coding region (AT-rich region) in the mt genome of G. pulchrum which is located between cox3 and tRNA-Ala (Figure 1 and Table 2), with 81.11% of A+T content (Table 3). Furthermore, this region is devoid of consecutive sequences of [A] and [T], and there are no AT dinucleotide repeat sequences; these repeat regions have been reported in the mt genome of A. simplex s.l. and S. digitata [32,34].
Identification and differentiation of G. pulchrum has traditionally been based on morphological features. However, these criteria are often insufficient for specific identification and differentiation, particularly at the larval and/or egg stages. Molecular tools, using genetic markers in mt cox1 and ITS-1 region of nuclear rDNA, have been used to support clinical diagnosis and to assist in undertaking molecular identification and epidemiological investigations of G. pulchrum . Because sequence polymorphism (heterogeneity) in ITS rDNA sequences occurs within individual spirurid specimens , mt genome sequences appear to be better solutions for such studies. Additionally, the cox1 sequences of G. pulchrum were further divided into multiple haplotypes and two groups of haplotypes (i.e. those from a majority of sika deer, wild boars and Japanese macaques and those from cattle and zoo animals, were clearly differentiated) . Nonetheless, the cox1 is a relatively conserved mt gene in nematodes [13,16,39], and, to date, there is no genetic information for G. pulchrum from other mt genes.
Many studies have indicated that the mtDNA sequence is a valuable genetic marker for phylogenetic studies [11,12,14,17]. The mt genome sequence of G. pulchrum could promote to reassess the systematic relationships within the spirurid nematodes using mt genomic datasets. Over the last decades, there have been considerable debate concerning the systematics of members of the spirurid nematodes (including superfamilies Filarioidea, Physalopteroidea, Spiruroidea and Thelazoidea) . Some studies using nuclear small subunit (SSU) rDNA and mtDNA sequences have indicated that S. lupi (Spirocercidae) is the sister taxon of T. callipaeda (Thelaziidae), suggesting that S. lupi belongs to the superfamily Thelazioidea [41,42], but this finding was contradicted by other studies that used the same markers [43,44]. The results of the present study support that the Spirocercidae (represented by S. lupi) was more closely related to the family Gongylonematidae (represented by G. pulchrum) than to other families within Spiruromorpha, indicating that both families Spirocercidae and Gongylonematidae belong to superfamily Spiruroidea, consistent with conclusions of a previous study . Given this utility of mt genomic datasets, further work should include sequencing of mt genomes of other spirurid nematodes in order to reconstruct the phylogenetic relationships of spirurid nematodes.
The present study determined the complete mt genome sequence of G. pulchrum, and ascertained its phylogenetic position within the infraorder Spiruromorpha. The complete mt genome represents the first sequenced mt genome of any member of the family Gongylonematidae. It will provide an important resource for the design of novel primers for the study of epidemiology, population genetics and systematics of Gongylonematidae.
This work was supported, in part, by grants from Funds of Basic Research Key Program (ZD2012010) to GHZ, the International Science & Technology Cooperation Program of China (Grant No. 2013DFA31840) and the Science Fund for Creative Research Groups of Gansu Province (Grant No. 1210RJIA006) to XQZ.
- Cappucci Jr DT, Augsburg JK, Klinck PC. Gongylonemiasis. In: Steele JH, editor. Handbook Series in Zoonoses, Section C: Parasitic Zoonoses, vol. II. Boca Raton, Florida: CRC Press; 1982. p. 181–92.Google Scholar
- Anderson RC: Nematode parasites of vertebrates; their development and transmission. C.A.B. International, Oxon, 1992.Google Scholar
- Kudo N, Kubota H, Gotoh H, Ishida H, Ikadai H, Oyamada T. Efficacy of thiabendazole, mebendazole, levamisole and ivermectin against gullet worm, Gongylonema pulchrum: In vitro and in vivo studies. Vet Parasitol. 2008;151:46–52.View ArticlePubMedGoogle Scholar
- Zhang ZW, Zhang HF, Gao CY. The endoscopic diagnosis and treatment of esophageal gongylonemiasis. Chin Foreign Med Treat. 2011;6:31 (in Chinese).Google Scholar
- Battistelli-Lux C. Buccal infection with Gongylonema pulchrum: an indigenous case in France. Ann Dermatol Venereol. 2013;140:623–7 (In France).View ArticlePubMedGoogle Scholar
- Pesson B, Hersant C, Biehler JF, Abou-Bacar A, Brunet J, Pfaff AW, et al. First case of human gongylonemosis in France. Parasite. 2013;20:5.View ArticlePubMed CentralPubMedGoogle Scholar
- Allen JD, Esquela-Kerscher A. Gongylonema pulchrum Infection in a Resident of Williamsburg, Virginia, Verified by Genetic Analysis. Am J Trop Med Hyg. 2013;89:755–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Makouloutou P, Setsuda A, Yokoyama M, Tsuji T, Saita E, Torii H, et al. Genetic variation of Gongylonema pulchrum from wild animals and cattle in Japan based on ribosomal RNA and mitochondrial cytochrome c oxidase subunit I genes. J Helminthol. 2013;87:326–35.View ArticlePubMedGoogle Scholar
- Wolstenholme DR. Animal mitochondrial DNA, structure and evolution. Int Rev Cytol. 1992;141:173–216.View ArticlePubMedGoogle Scholar
- Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80.View ArticlePubMed CentralPubMedGoogle Scholar
- Kang S, Sultana T, Eom KS, Park YC, Soonthornpong N, Nadler SA, et al. The mitochondrial genome sequence of Enterobius vermicularis (Nematoda: Oxyurida)–an idiosyncratic gene order and phylogenetic information for chromadorean nematodes. Gene. 2009;429:87–97.View ArticlePubMedGoogle Scholar
- Park JK, Sultana T, Lee SH, Kang S, Kim HK, Min GS, et al. Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics. 2011;12:392.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu GH, Gasser RB, Su A, Nejsum P, Peng L, Lin RQ, et al. Clear genetic distinctiveness between human- and pig-derived Trichuris based on analyses of mitochondrial datasets. PLoS Negl Trop Dis. 2012;6:e1539.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu GH, Shao R, Li JY, Zhou DH, Li H, Zhu XQ. The complete mitochondrial genomes of three parasitic nematodes of birds: a unique gene order and insights into nematode phylogeny. BMC Genomics. 2013;14:414.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu GH, Gasser RB, Otranto D, Xu MJ, Shen JL, Mohandas N, et al. Mitochondrial genome of the eyeworm, Thelazia callipaeda (nematoda: spirurida), as the first representative from the family Thelaziidae. PLoS Negl Trop Dis. 2013;7:e2029.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu GH, Zhao L, Song HQ, Zhao GH, Cai JZ, Zhao Q, et al. Chabertia erschowi (Nematoda) is a distinct species based on nuclear ribosomal DNA sequences and mitochondrial DNA sequences. Parasit Vectors. 2014;7:44.View ArticlePubMed CentralPubMedGoogle Scholar
- Sultana T, Kim J, Lee SH, Han H, Kim S, Min GS, et al. Comparative analysis of complete mitochondrial genome sequences confirms independent origins of plant-parasitic nematodes. BMC Evol Biol. 2013;13:12.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhao GH, Hu B, Cheng WY, Jia YQ, Li HM, Yu SK, et al. The complete mitochondrial genomes of Oesophagostomum asperum and Oesophagostomum columbianum in small ruminants. Infect Genet Evol. 2013;19:205–11.View ArticlePubMedGoogle Scholar
- Gao JF, Zhao Q, Liu GH, Zhang Y, Zhang Y, Wang WT, et al. Comparative analyses of the complete mitochondrial genomes of the two ruminant hookworms Bunostomum trigonocephalum and Bunostomum phlebotomum. Gene. 2014;541:92–100.View ArticlePubMedGoogle Scholar
- Sato H, Une Y, Takada M. High incidence of the gullet worm, Gongylonema pulchrum, in a squirrel monkey colony in a zoological garden in Japan. Vet Parasitol. 2005;127:131–7.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. Nature Protoc. 2007;2:2339–44.View ArticleGoogle 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
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.View ArticlePubMed CentralPubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA 5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.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:145–58.View ArticlePubMedGoogle Scholar
- Liu GH, Zhou DH, Zhao L, Xiong RC, Liang JY, Zhu XQ. The complete mitochondrial genome of Toxascaris leonina: comparison with other closely related species and phylogenetic implications. Infect Genet Evol. 2014;21:329–33.View ArticlePubMedGoogle Scholar
- Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56:564–77.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21:2104–5.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.View ArticlePubMedGoogle Scholar
- Rambaut A: FigTree, a graphical viewer of phylogenetic trees. 2012. Available from http://tree.bio.ed.ac.uk/software/figtree.
- Yatawara L, Wickramasinghe S, Rajapakse RP, Agatsuma T. The complete mitochondrial genome of Setaria digitata (Nematoda: Filarioidea): Mitochondrial gene content, arrangement and composition compared with other nematodes. Mol Biochem Parasitol. 2010;173:32–8.View ArticlePubMedGoogle Scholar
- Li MW, Lin RQ, Song HQ, Wu XY, Zhu XQ. The complete mitochondrial genomes for three Toxocara species of human and animal health significance. BMC Genomics. 2008;9:224.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim KH, Eom KS, Park JK. The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications. Int J Parasitol. 2006;36:319–28.View ArticlePubMedGoogle Scholar
- Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR. The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992;130:471–98.PubMed CentralPubMedGoogle Scholar
- Jex AR, Hall RS, Littlewood DT, Gasser RB. An integrated pipeline for next-generation sequencing and annotation of mitochondrial genomes. Nucleic Acids Res. 2010;38:522–33.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu M, Gasser RB. Mitochondrial genomes of parasitic nematodes–progress and perspectives. Trends Parasitol. 2006;22:78–84.View ArticlePubMedGoogle Scholar
- Gasser RB, LeGoff L, Petit G, Bain O. Rapid delineation of closely-related filarial parasites using genetic markers in spacer rDNA. Acta Trop. 1996;62:143–50.View ArticlePubMedGoogle Scholar
- Liu GH, Gasser RB, Nejsum P, Wang Y, Chen Q, Song HQ, et al. Mitochondrial and nuclear ribosomal DNA evidence supports the existence of a new Trichuris species in the endangered françois’ leaf-monkey. PLoS One. 2013;8:e66249.View ArticlePubMed CentralPubMedGoogle Scholar
- De Ley P, Blaxter M. Systematic position and phylogeny. In: Lee DL, editor. The Biology of Nematodes. London and New York.: Taylor & Francis; 2002. p. 1–30.View ArticleGoogle Scholar
- Iorio R, Slapeta J, Otranto D, Paoletti B, Giangaspero A, Traversa D. Phylogenetic relationships of Habronema microstoma and Habronema muscae (Spirurida: Habronematidae) within the order Spirurida inferred using mitochondrial cytochrome c oxidase subunit 1 (cox1) gene analysis. Parasitol Res. 2009;104:979–84.View ArticlePubMedGoogle Scholar
- Černotíková E, Horák A, Moravec F. Phylogenetic relationships of some spirurine nematodes (Nematoda: Chromadorea: Rhabditida: Spirurina) parasitic in fishes inferred from SSU rRNA gene sequences. Folia Parasitol. 2011;58:135–48.View ArticlePubMedGoogle Scholar
- Traversa D, Costanzo F, Iorio R, Aroch I, Lavy E. Mitochondrial cytochrome c oxidase subunit 1 (cox1) gene sequence of Spirocerca lupi (Nematoda, Spirurida): avenues for potential implications. Vet Parasitol. 2007;146:263–70.View ArticlePubMedGoogle Scholar
- Nadler SA, Carreno RA, Mejía-Madrid H, Ullberg J, Pagan C, Houston R, et al. Molecular phylogeny of clade III nematodes reveals multiple origins of tissue parasitism. Parasitology. 2007;134:1421–42.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.