Figures
Abstract
Human thelaziosis is an underestimated parasitic disease caused by Thelazia species (Spirurida: Thelaziidae). The oriental eyeworm, Thelazia callipaeda, infects a range of mammalian definitive hosts, including canids, felids and humans. Although this zoonotic parasite is of socio-economic significance in Asian countries, its genetics, epidemiology and biology are poorly understood. Mitochondrial (mt) DNA is known to provide useful genetic markers to underpin fundamental investigations, but no mt genome had been characterized for any members of the family Thelaziidae. In the present study, we sequenced and characterized the mt genome of T. callipaeda. This AT-rich (74.6%) mt genome (13,668 bp) is circular and contains 12 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes, but lacks an atp8 gene. All protein-coding genes are transcribed in the same direction; the gene order is the same as those of Dirofilaria immitis and Setaria digitata (Onchocercidae), but distinct from Dracunculus medinensis (Dracunculidae) and Heliconema longissimum (Physalopteridae). Phylogenetic analyses of the concatenated amino acid sequence data for all 12 protein-coding genes by Bayesian inference (BI) showed that T. callipaeda (Thelaziidae) is related to the family Onchocercidae. This is the first mt genome of any member of the family Thelaziidae and should represent a new source of genetic markers for studying the epidemiology, ecology, population genetics and systematics of this parasite of humans and other mammals.
Author Summary
Human thelaziosis is an underestimated parasitic disease caused by the eyeworm Thelazia callipaeda (Spirurida: Thelaziidae). Although this parasite is of significance in humans in many Asian countries, its genetics, epidemiology and biology are poorly understood. Mitochondrial (mt) DNA can provide useful genetic markers for fundamental investigations, but no mt genome had been characterized for any members of the family Thelaziidae. In this study, we sequenced and characterized the mt genome of T. callipaeda. This circular mt genome is 13,668 bp long and contains 12 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes, but lacks an atp8 gene. Phylogenetic analyses of the concatenated amino acid sequence data for all 12 protein-coding genes by Bayesian inference showed that T. callipaeda is closely related to the family Onchocercidae, consistent with previous study. This is the first mt genome of any member of the family Thelaziidae, and represents a new source of genetic markers for studies of the epidemiology, ecology, population genetics and systematics of this parasite of human and animal health significance.
Citation: Liu G-H, Gasser RB, Otranto D, Xu M-J, Shen J-L, Mohandas N, et al. (2013) Mitochondrial Genome of the Eyeworm, Thelazia callipaeda (Nematoda: Spirurida), as the First Representative from the Family Thelaziidae. PLoS Negl Trop Dis 7(1): e2029. https://doi.org/10.1371/journal.pntd.0002029
Editor: Peter U. Fischer, Washington University School of Medicine, United States of America
Received: September 18, 2012; Accepted: December 7, 2012; Published: January 31, 2013
Copyright: © 2013 Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported in part by the National S & T Major Program (Grant No. 2012ZX10004220) and the Program for Outstanding Scientists in Agricultural Research. The present study was also supported by the Australian Research Council (RBG). Current research in the Gasser Lab is funded mainly by the Australian Research Council (ARC), the National Health and Medical Research Council (NHMRC) and Melbourne Water Corporation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Thelazia callipaeda Railliet and Henry, 1910, known as the ‘oriental eye-worm’, because of its geographical distribution in Asian countries (including China, India, Japan, Korea and Thailand), is frequently reported as being responsible for thelaziosis of humans, carnivores (dogs, foxes and cats) and rabbits, causing mild to severe clinical signs (including lacrimation, epiphora, conjunctivitis, keratitis and/or sometimes corneal ulcers) [1]. Fortunately, thelaziosis can be treated effectively using anthelminthics, such as milbemycin oxime or macrocyclic lactones (e.g., moxidectin), and anti-inflammatory compounds [2]–[4]. Although T. callipaeda may seem to be of minor importance to some clinicians and scientists, human thelaziosis is highly endemic in some under-developed communities in Asia, particularly in China [3]. Clearly, scant attention has been paid to human thelaziosis, and there are difficulties in its clinical diagnosis and differentiation from allergic conjunctivitis, particularly when small numbers of adult or larval stages of T. callipaeda are present in the eyes of infected patients.
The transmission of human thelaziosis occurs when the intermediate host, a drosophilid fly of the genus Phortica, feeds on lacrimal secretions from humans and other animals, and ingests first-stage larvae (L1s) produced by adult females of T. callipaeda, which live together with males in the conjunctival sac. After being ingested by the fly, the T. callipaeda larvae migrate in the vector's body (i.e. testis of the male) and undergo development from the L1 to the infective, third-stage larvae (L3) within 14–21 days. Following this migration, the L3s of Thelazia emerge from the labella of the infected fly, are deposited on the eye, as the vector feeds on lacrimal secretions, and then develop into the dioecious adult stages in the ocular cavity within ∼35 days [5].
In spite of the significance of human thelaziasis, little is known about the biology and epidemiology of T. callipaeda and its close relatives [1], [3]. This relates mainly to a lack of reliable morphological characters for their specific identification and for comparative study. Although molecular tools, employing genetic markers in short regions of nuclear ribosomal and mitochondrial (mt) DNA, have found utility for taxonomic and epidemiological studies of some species, such as T. gulosa, T. rhodesi, T. skrjabini, T. lacrymalis and T. callipaeda [1], there is still a paucity of information on T. callipaeda in different human populations and countries around the world.
Mt genomes can provide markers for genetic and epidemiological investigations of spirurid nematodes (e.g., [6], [7]), and provide the potential to discover population variants or cryptic species and investigate transmission patterns linked to particular haplotypes [8]–[11]. In addition, mt proteomic datasets could be used for reassessing systematic relationships of Thelazia and other spirurids. Recent studies have shown that concatenated mt proteomic datasets can be used effectively to retest hypotheses regarding the systematic relationships of different groups of nematodes (e.g., [12], [13]). Such amino acid sequence datasets are relatively large and informative, usually achieving excellent phylogenetic signal and mostly attaining nodal support values of 98–100% in tree reconstructions [12], [13]. Long-range PCR-coupled sequencing and bioinformatic methods [14], [15] have underpinned these advances, which now allow mt proteomic barcodes to be defined rapidly for Thelazia spp. and their relatives from a range of mammalian and invertebrate hosts. Here, as a first step, we (i) determined the sequence and structure of the mt genome for T. callipaeda, (ii) assessed the phylogenetic position of this nematode in relation to other Spirurida for which whole mt sequence datasets are available, and (iii) discussed the implications of the new dataset as a new resource for future genetic studies of T. callipaeda populations.
Materials and Methods
Ethics statement
This study was approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Approval No. LVRIAEC2010-008). The dog from which the adult specimens of T. callipaeda were collected was handled in accordance with good animal practice (GAP) required by the Animal Ethics Procedures and Guidelines of the People's Republic of China.
Parasites and total genomic DNA isolation
Adult specimens of T. callipaeda were collected from a conjunctival sac of an infected dog at a veterinary hospital in Zhanjiang, Guangdong Province, China. The worms were washed extensively in physiological saline, fixed in ethanol and then stored at −20°C until use. Upon thawing, the anterior and posterior ends of each nematode were cut off and cleared in lactophenol for subsequent morphological identification [16]. The mid-body section of each worm was used for the isolation of total genomic DNA by small-scale sodium dodecyl-sulphate (SDS)/proteinase K digestion [17] and mini-column purification (TIANamp Genomic DNA kit). The molecular identity of each specimen was then verified by PCR-based sequencing of regions in the cox1 and rrnS genes using an established method [18], [19], and both regions had 99% identity to previously published sequences for T. callipaeda from China and Italy (GenBank accession nos. AM042555 and AJ544858, respectively).
Long-PCR, sequencing and annotation
Using primers (Table S1) designed to relatively conserved regions within the cox1 and rrnS regions (see Figure 1), the complete mt genome was amplified by long-PCR as two overlapping amplicons (∼5 kb and ∼9 kb) from the genomic DNA from the mid-body section of a single female specimen of T. callipaeda. PCR was conducted in 25 µl using 2 mM MgCl2, 0.2 mM each of dNTPs, 2.5 µl 10× Taq buffer, 2.5 µM of each primer and 0.5 µl LA Taq DNA polymerase (5 U/µl, Takara) in a thermocycler (Biometra) under the following conditions: 92°C for 2 min (initial denaturation), then 92°C for 10 s (denaturation), 58°C (5 kb) or 42°C (9 kb) for 30 s (annealing), and 60°C for 10 min (extension) for 10 cycles, followed by 92°C for 10 s, 58°C (∼5 kb) or 42°C (∼9 kb) for 30 s (annealing), and 60°C for 10 min for 20 cycles, with a cycle elongation of 10 s for each cycle and a final extension at 60°C for 10 min. Genomic DNA (80 ng in 2 µl) was added to PCR, and no-template and known-positive controls were included in each run. Amplicons were column-purified (Wizard PCR Preps, Promega). Subsequently, the amount of DNA in each purified amplicon was estimated spectrophotometrically (ND-1000 UV-VIS spectrophotometer, v.3.2.1, NanoDrop Technologies). Following an electrophoretic analysis of quality, purified amplicons were sequenced using a primer walking strategy [20]. The whole mt genome sequence (GenBank accession no. JX069968) was then assembled using the ContigExpress program of the Vector NTI software package v.6.0 (Invitrogen, Carlsbad, CA).
Scale is approximate. All genes have standard nomenclature except for the 22 tRNA genes, which are designated by the one-letter code for the corresponding amino acid, with numerals differentiating each of the two leucine- and serine-specifying tRNAs (L1 and L2 for codon families CUN and UUR, respectively; S1 and S2 for codon families UCN and AGN, respectively). All genes are transcribed in the clockwise direction. ‘AT’ indicates the non-coding region.
The mt genome was annotated using an approach similar to that of Yatawara et al. [21]. In brief, each protein-encoding mt gene was identified by local alignment comparison using amino acid sequences conceptually translated from corresponding genes from the mt genome of a reference species (i.e. Setaria digitata; accession number: NC_014282) [21]. The tRNA (trn) genes were identified using the program tRNAscan-SE [22] or by visual inspection [23]; rRNA (rrn) genes were predicted by comparison with those of S. digitata [21].
Phylogenetic analysis of concatenated amino acid sequence data
The amino acid sequences conceptually translated from individual genes of the mt genome of T. callipaeda were concatenated. Selected for comparison were concatenated amino acid sequences predicted from published mt genomes from key nematodes representing the order Spirurida, including the superfamily Filarioidea (Acanthocheilonema viteae [24], Brugia malayi [25], Chandlerella quiscali [24], Dirofilaria immitis [26], Loa loa [24], Onchocerca flexuosa [24], O. volvulus [27], S. digitata [21] and Wuchereria bancrofti [7]), the superfamily Dracunculoidea (Dracunculus medinensis) and the superfamily Physalopteroidea (Heliconema longissimum) [13] (GenBank accession numbers NC_016197, NC_004298, NC_014486, NC_005305, NC_016199, NC_016172, AF015193, NC_014282, JN367461, NC_016019 and NC_016127, respectively), using Ascaris suum (GenBank accession number HQ704901) as an outgroup [28]. All amino acid sequences (considering all homologous characters) were aligned using MUSCLE [29] and then subjected to phylogenetic analysis using BI as described previously [30], [31]. Phylograms were drawn using the Tree View program v.1.65 [32].
Results and Discussion
General features of the mt genome of T. callipaeda
The complete mt genome sequence of T. callipaeda (GenBank accession no. JX069968) was 13,668 bp in length (Figure 1). This genome contains 12 protein-coding genes (cox1–3, nad1–6, nad4L, atp6 and cytb), 22 trn genes, two rrn genes (rrnL and rrnS) and a non-coding (control or AT-rich) region, but lacks an atp8 gene (Table 1). The gene content and arrangement are the same as those of D. immitis and S. digitata [21], [26], but distinct from those of D. medinensis (rearrangement markedly) and H. longissimum (tRNA-Met and tRNA-Val change) [13]. All genes are transcribed in the same direction. In addition, the mt genes of T. callipaeda overlap by 102 nt in 14 locations (1 to 32 nt per location) (Table 1). The mt genome of T. callipaeda has 14 intergenic regions, which range from 1 to 62 nt in length. The longest region is between tRNA-Pro and tRNA-Asp (Table 1).
The nucleotide content of the entire mt genome sequence of T. callipaeda is biased toward A+T (74.6%), in accordance with mt genomes of other nematodes of the order Spirurida (e.g., [21], [26]) (Table 2). One non-coding region (AT-loop) (328 bp), located between cox3 and tRNA-Ala, has the highest A+T content of 79.6% (Table 2). AT- and GC-skews of the whole mt genome were calculated for T. callipaeda and other spirurid nematodes studied to date (see Table 3). This composition of the mt genome sequence of T. callipaeda was strongly skewed away from A, in favour of T (AT skew = −0.40), and the GC skew was 0.449 (Table 3). All spirurid nematodes reported to date and in the present study show strand asymmetry (GC skew between 0.354 and 0.521) (Table 3).
Protein-encoding genes
The boundaries between protein-coding genes of the mt genome of T. callipaeda were determined by aligning its sequence and by identifying translation initiation and termination codons with those of H. longissimum and S. digitata [13], [21]. In this mt genome, all protein-coding genes had ATT, ATA and TTG as their initiation codons, and TAA or TAG as their termination codon. Incomplete termination codons (T or TA) were not identified, which is inconsistent with studies of some other nematodes, including Anisakis simplex (s. l.), A. suum, Caenorhabditis elegans, S. digitata, Toxocara spp. and Trichinella spiralis [21], [33]–[36]. Codons composed of A and T were more frequently used in protein-coding genes, reflecting the high A+T content in the mt genome of T. callipaeda. The most frequently used amino acid was Phe (19.3%), followed by Leu (13.4%), Val (7.6%), Gly (7.1%) and IIe (6.2%) (Table 4).
Other genes
In the mt genome of T. callipaeda, the rrnL was located between tRNA-His and nad3, and rrnS was between nad4L and tRNA-Tyr (Table 1). The sizes of the rrnL and rrnS genes of T. callipaeda were 966 bp and 666 bp, respectively (Table 1). The 22 trn genes ranged from 52 to 66 bp in size. The secondary structures predicted for the latter genes were similar to those of S. digitata [21].
Substitution ratios
As synonymous and non-synonymous substitution rates assist in predicting evolutionary processes [37], the rate of non-synonymous substitutions (Ka), the rate of synonymous substitutions (Ks) and the Ka/Ks ratios were calculated for all 12 protein-coding genes encoded in the mt genomes of T. callipaeda and 11 other spirurid nematodes, including A. viteae, B. malayi, C. quiscali, L. loa, S. digitata and W. bancrofti (Table 3). The Ka/Ks ratio is a measure of selective pressures acting on gene that indicates neutral mutation (ka/ks = 1), negative or purifying selection (Ka/Ks of <1), and positive or diversifying selection (Ka/Ks of >1) [38], [39]. Here, nad2 showed the highest ratio, followed by nad3, while cox1 appeared to have the lowest ratio (Figure 2). Notably, the Ka/Ks ratio of eight protein-coding genes was <1 (range: 0.346 to 0.873), indicating that these genes are evolving under negative or purifying selection [40], [41]. The Ka/Ks ratio of 4 protein-coding genes (nad2, nad3, nad5 and nad6) was >1 (range: 1.105 to 1.331), suggesting that these genes have evolved under positive or diversifying selection [42].
The rate of non-synonymous (Ka), the rate of synonymous (Ks) substitutions, and the respective ratios (Ka/Ks) for individual protein-coding genes are shown.
Sequence comparisons and phylogenetic relationships of T. callipaeda with selected members of the Spirurida
The amino acid sequences predicted from individual protein-coding mt genes of T. callipaeda were compared with those of 11 other spirurid nematodes (see Table 5). Pairwise comparisons of the concatenated amino acid sequences revealed identities of 40.3–91.8% among them. Based on identity, COX1 was the most conserved protein, whereas nad4L and nad3 were the least conserved (see Table 5). Phylogenetic analysis of the concatenated amino acid sequence data for all 12 mt proteins showed that T. callipaeda (Thelaziidae) was a sister taxon to a clade containing S. digitata (Setariidae) and other members of the Onchocercidae, including B. malayi and D. immitis, consistent with results of a previous study [18]. Basal to these taxa were H. longissimum (Physalopteridae) and D. medinensis (Dracunculidae) (posterior probability = 1.00) (Figure 3).
The concatenated amino acid sequences of 12 protein-coding genes were subjected to analysis by Bayesian inference (BI) using Ascaris suum as an outgroup. Posterior probability (pp) values are indicated.
Fundamental and applied implications
Although much attention has been paid to soil-transmitted helminths as pathogens because of their major socioeconomic impact on human populations [43], [44], parasitic nematodes that cause relatively subtle, but chronic disease, such as members of the genus Thelazia, have been seriously neglected [3], [45]. The main reservoirs for human thelaziosis seem to be dogs, since they often live in areas populated by a large entomo-fauna [1], [46]. T. callipaeda is usually prevalent in dogs, cats and humans in disadvantaged, rural areas of the former Soviet Union [47] and the Asian continent, including China [48], India [49], Indonesia [50], Japan [51], Korea [52], Taiwan [53] and Thailand [54], [55]. More recently, T. callipaeda has also been reported in Europe, with a high prevalence (60%) in dogs being recorded in some areas of Southern Italy [56]. Autochthonous cases of canine thelaziosis have also been recorded in France [57], [58], Portugal [59], Spain [60] and Switzerland [61], suggesting that the latitude range of endemicity of canine thelaziosis in Europe (between 39° and 46°N) is similar to that of Asia (between 10° and 45°N for India and Japan) [56]. Interestingly, in spite of the high prevalence of canine thelaziosis reported for southern parts of Europe [56], only a small number of human cases have yet been reported in this geographical region [45].
In the present study, the characterization of the mt genome of T. callipaeda provides a foundation for the improved diagnosis of human thelaziosis using molecular methods as well as future, detailed studies of the population genetics and epidemiology/ecology of this parasite in Asia. As adult and larval stages of T. callipaeda from the eyes of patients cannot be identified reliably by morphology to species, molecular tools, using genetic markers in the first internal transcribed spacer (ITS-1) region of nuclear rDNA and cox1, have been used to support clinical diagnosis and to assist in undertaking molecular epidemiological investigations of T. callipaeda [1], [19]. Because sequence heterogeneity in ITS rDNA can be high in individual spirurid specimens (e.g., [62]), sometimes complicating sequence analyses, protein-coding mt genes appear to be better suited for such studies [19].
Having available the mt genome of T. callipaeda now sets the scene to develop combined DNA-based analytical and diagnostic tools, whereby mt genetic regions with differing levels of within-species divergence [19] might be used to explore haplotypic variation of individuals within and among T. callipaeda populations infecting humans and other definitive hosts as well as fly intermediate hosts, such as Phortica variegata. This could be done effectively using PCR-coupled mutation scanning and selective sequencing [17], already effectively applied, on a small scale, to T. callipaeda [19]. A previous investigation, employing cox1 alone, showed that, despite a relatively high degree of genetic variability among specimens isolated from Asia (i.e. China and Korea), no genetic variation was detected among individual specimens from different host species (i.e. dogs, cats and foxes) and localities within Europe (i.e. France, Germany, Italy, Netherlands, and Spain) [19]. These data were supported by additional studies [58]–[60], suggesting a genetically homogenized population in Europe, a tighter affiliation of this nematode to intermediate hosts than to the definitive hosts, and, thus, that the distribution of the parasite might be expected to resemble that of the vector [19]. Nonetheless, cox1 is a relatively conserved mt gene [63], and, to date, there is no genetic information for T. callipaeda from humans.
In the future, it would be interesting to assess whether various haplotypes or genotypes of T. callipaeda might relate to different clinical symptoms of thelaziosis in humans, and whether particular subpopulations of T. callipaeda undergo arrested development (hypobiosis, hibernation or aestivation) and survive for long periods of time in their intermediate hosts, as indicated in southern Europe [5]. Ecological aspects would be interesting to study in T. callipaeda in different countries and hosts, particularly dogs and cats, given the apparent complexity of the parasite's life cycle and biology. In addition, there is an applied human disease management imperative for improved molecular tools in being able to track the transmission of T. callipaeda and address the question as to whether infected domestic and wild canids and felids represent reservoir hosts for infection to humans in the same natural environment.
The characterization of the mt genome of T. callipaeda also stimulates a reassessment of the systematic relationships of nematodes within the order Spirurida using mt genomic/proteomic datasets. For decades, there have been controversies surrounding the systematics of members of the Spirurida (including superfamilies Acuarioidea, Aproctoidea, Diplotriaenoidea, Filarioidea, Gnathostomatoidea, Habronematoidea, Physalopteroidea, Rictularioidea and Spiruroidea) [64], [65]. Given the demonstrated utility of mt proteomic datasets, high phylogenetic signal and strong statistical support in trees [12], [13], there is now an opportunity to test the phylogenetic relationships of a wide range of spirurid nematodes using expanded mt datasets.
Supporting Information
Table S1.
Sequences of oligonucleotide primers for amplifying regions of the mitochondrial genome of Thelazia callipaeda.
https://doi.org/10.1371/journal.pntd.0002029.s001
(DOCX)
Author Contributions
Conceived and designed the experiments: XQZ RBG GHL. Performed the experiments: GHL MJX. Analyzed the data: GHL RBG XQZ. Contributed reagents/materials/analysis tools: MJX JLS DHZ NM. Wrote the paper: GHL RBG DO XQZ.
References
- 1. Otranto D, Traversa D (2004) Molecular characterization of the first internal transcribed spacer of ribosomal DNA of the most common species of eyeworms (Thelazioidea: Thelazia). J Parasitol 90: 185–188.
- 2. Bianciardi P, Otranto D (2005) Treatment of dog thelaziosis caused by Thelazia callipaeda (Spirurida, Thelaziidae) using a topical formulation of imidacloprid 10% and moxidectin 2.5%. Vet Parasitol 129: 89–93.
- 3. Shen J, Gasser RB, Chu D, Wang Z, Yuan X, et al. (2006) Human thelaziosis–a neglected parasitic disease of the eye. J Parasitol 92: 872–875.
- 4. Motta B, Schnyder M, Basano FS, Nägeli F, Nägeli C, et al. (2012) Therapeutic efficacy of milbemycin oxime/praziquantel oral formulation (Milbemax) against Thelazia callipaeda in naturally infested dogs and cats. Parasit Vectors 5: 85.
- 5. Otranto D, Lia RP, Cantacessi C, Testini G, Troccoli A, et al. (2005) Nematode biology and larval development of Thelazia callipaeda (Spirurida, Thelaziidae) in the drosophilid intermediate host in Europe and China. Parasitology 131: 847–855.
- 6. Unnasch TR, Williams SA (2000) The genomes of Onchocerca volvulus. Int J Parasitol 30: 543–552.
- 7. Ramesh A, Small ST, Kloos ZA, Kazura JW, Nutman TB, et al. (2012) The complete mitochondrial genome sequence of the filarial nematode Wuchereria bancrofti from three geographic isolates provides evidence of complex demographic history. Mol Biochem Parasitol 183: 32–41.
- 8. Traversa D, Costanzo F, Iorio R, Aroch I, Lavy E (2007) Mitochondrial cytochrome c oxidase subunit 1 (cox1) gene sequence of Spirocerca lupi (Nematoda, Spirurida): avenues for potential implications. Vet Parasitol 146: 263–270.
- 9. Iorio R, Slapeta J, Otranto D, Paoletti B, Giangaspero A, et al. (2009) 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 104: 979–984.
- 10. Wu SG, Wang GT, Xi BW, Xiong F, Liu T, et al. (2009) Population genetic structure of the parasitic nematode Camallanus cotti inferred from DNA sequences of ITS1 rDNA and the mitochondrial COI gene. Vet Parasitol 164: 248–256.
- 11. Otranto D, Latrofa MS, Brianti E, Annoscia G, Parisi A, et al. (2012) An assessment of genetic variability in the mitochondrial cytochrome c oxidase subunit 1 gene of Cercopithifilaria sp. (Spirurida, Onchocercidae) from dog and Rhipicephalus sanguineus populations. Mol Cell Probes 26: 81–89.
- 12. Jex AR, Hall RS, Littlewood DT, Gasser RB (2010) An integrated pipeline for next-generation sequencing and annotation of mitochondrial genomes. Nucleic Acids Res 38: 522–533.
- 13. Park JK, Sultana T, Lee SH, Kang S, Kim HK, et al. (2011) Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics 12: 392.
- 14. Hu M, Gasser RB (2006) Mitochondrial genomes of parasitic nematodes–progress and perspectives. Trends Parasitol 22: 78–84.
- 15. Liu GH, Gasser RB, Su A, Nejsum P, Peng L, et al. (2012) Clear genetic distinctiveness between human- and pig-derived Trichuris based on analyses of mitochondrial datasets. PLoS Negl Trop Dis 6: e1539.
- 16. Otranto D, Lia RP, Traversa D, Giannetto S (2003) Thelazia callipaeda (Spirurida, Thelaziidae) of carnivores and humans: morphological study by light and scanning electron microscopy. Parassitology 45: 125–133.
- 17. Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AJ, et al. (2006) Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc 1: 3121–3128.
- 18. Casiraghi M, Bain O, Guerrero R, Martin C, Pocacqua V, et al. (2004) Mapping the presence of Wolbachia pipientis on the phylogeny of filarial nematodes: evidence for symbiont loss during evolution. Int J Parasitol 34: 191–203.
- 19. Otranto D, Testini G, De Luca F, Hu M, Shamsi S, et al. (2005) Analysis of genetic variability within Thelazia callipaeda (Nematoda: Thelazioidea) from Europe and Asia by sequencing and mutation scanning of the mitochondrial cytochrome c oxidase subunit 1 gene. Mol Cell Probes 19: 306–313.
- 20. Hu M, Jex AR, Campbell BE, Gasser RB (2007) Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing. Nature Protoc 2: 2339–2344.
- 21. Yatawara L, Wickramasinghe S, Rajapakse RP, Agatsuma T (2010) The complete mitochondrial genome of Setaria digitata (Nematoda: Filarioidea): Mitochondrial gene content, arrangement and composition compared with other nematodes. Mol Biochem Parasitol 173: 32–38.
- 22. Lowe TM, Eddy SR (1997) tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 25: 955–964.
- 23. Hu M, Chilton NB, Gasser RB (2002) The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea). Int J Parasitol 32: 145–158.
- 24. McNulty SN, Mullin AS, Vaughan JA, Tkach VV, Weil GJ, et al. (2012) Comparing the mitochondrial genomes of Wolbachia-dependent and independent filarial nematode species. BMC Genomics 13: 145.
- 25. Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, et al. (2007) Draft genome of the filarial nematode parasite Brugia malayi. Science 317: 1756–1760.
- 26. Hu M, Gasser RB, Abs EL-Osta YG, Chilton NB (2003) Structure and organization of the mitochondrial genome of the canine heartworm, Dirofilaria immitis. Parasitology 127: 37–51.
- 27. Keddie EM, Higazi T, Unnasch TR (1998) The mitochondrial genome of Onchocerca volvulus: sequence, structure and phylogenetic analysis. Mol Biochem Parasitol 95: 111–127.
- 28. Liu GH, Wu CY, Song HQ, Wei SJ, Xu MJ, et al. (2012) Comparative analyses of the complete mitochondrial genomes of Ascaris lumbricoides and Ascaris suum from humans and pigs. Gene 492: 110–116.
- 29. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.
- 30. Jex AR, Waeschenbach A, Hu M, van Wyk JA, Beveridge I, et al. (2009) The mitochondrial genomes of Ancylostoma caninum and Bunostomum phlebotomum-two hookworms of animal health and zoonotic importance. BMC Genomics 10: 79.
- 31. Liu GH, Wang SY, Huang WY, Zhao GH, Wei SJ, et al. (2012) The complete mitochondrial genome of Galba pervia (Gastropoda: Mollusca), an intermediate host snail of Fasciola spp. PLoS One 7: e42172.
- 32. Page RD (1996) TREEVIEW: an application to display phylogenetic trees on personal computers. Comput Appl Biosci 12: 357–358.
- 33. Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR (1992) The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics 130: 471–498.
- 34. Lavrov DV, Brown WM (2001) Trichinella spiralis mtDNA: A nematode mitochondrial genome that encodes a putative atp8 and normally structured tRNAs and has a gene arrangement relatable to those of coelomate metazoans. Genetics 157: 621–637.
- 35. Kim KH, Eom KS, Park JK (2006) The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications. Int J Parasitol 36: 319–328.
- 36. Li MW, Lin RQ, Song HQ, Wu XY, Zhu XQ (2008) The complete mitochondrial genomes for three Toxocara species of human and animal health significance. BMC Genomics 9: 224.
- 37. Yang Z, Nielsen R (2000) Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol 17: 32–43.
- 38. Meganathan PR, Dubey B, Batzer MA, Ray DA, Haque I (2011) Complete mitochondrial genome sequences of three Crocodylus species and their comparison within the order Crocodylia. Gene 478: 35–41.
- 39. Li H, Liu H, Shi A, Stys P, Zhou XG, et al. (2012) The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae). PLoS One 7: e29419.
- 40. Roques S, Fox CJ, Villasana MI, Rico C (2006) The complete mitochondrial genome of the whiting, Merlangius merlangus and the haddock, Melanogrammus aeglefinus: a detailed genomic comparison among closely related species of the Gadidae family. Gene 383: 12–23.
- 41. Yuan ML, Wei DD, Wang BJ, Dou W, Wang JJ (2010) The complete mitochondrial genome of the citrus red mite Panonychus citri (Acari: Tetranychidae): high genome rearrangement and extremely truncated tRNAs. BMC Genomics 11: 597.
- 42. Jordan IK, Rogozin IB, Wolf YI, Koonin EV (2002) Essential genes are more evolutionarily conserved than are nonessential genes in bacteria. Genome Res 12: 962–968.
- 43. Bethony JM, Cole RN, Guo X, Kamhawi S, Lightowlers MW, et al. (2011) Vaccines to combat the neglected tropical diseases. Immunol Rev 239: 237–270.
- 44. Hotez PJ, Savioli L, Fenwick A (2012) Neglected tropical diseases of the Middle East and North Africa: review of their prevalence, distribution, and opportunities for control. PLoS Negl Trop Dis 6: e1475.
- 45. Otranto D, Dutto M (2008) Human thelaziasis, Europe. Emerg Infect Dis 14: 647–649.
- 46. Seo M, Yu JR, Park HY, Huh S, Kim SK, et al. (2002) Enzooticity of the dogs, the reservoir host of Thelazia callipaeda, in Korea. Korean J Parasitol 40: 101–103.
- 47. Miroshnichenko VA, Desiaterik MP, Novik AP, Gorbach AP, Papernova NIu (1988) [A case of ocular thelaziasis in a 3-year-old child]. Vestn Oftalmol 104: 64 [Article in Russian].
- 48. Chen W, Zheng J, Hou P, Li L, Hu Y (2010) A case of intraocular thelaziasis with rhegmatogenous retinal detachment. Clin Exp Optom 93: 360–362.
- 49. Singh TS, Singh KN (1993) Thelaziasis: report of two cases. Br J Ophthalmol 77: 528–529.
- 50. Kosin E, Kosman ML, Depary AA (1989) First case of human Thelaziasis in Indonesia. Southeast Asian J Trop Med Public Health 20: 233–236.
- 51. Koyama Y, Ohira A, Kono T, Yoneyama T, Shiwaku K (2000) Five cases of thelaziasis. Br J Ophthalmol 84: 441.
- 52. Sohn WM, Na BK, Yoo JM (2011) Two cases of human thelaziasis and brief review of Korean cases. Korean J Parasitol 49: 265–271.
- 53. Cheung WK, Lu HJ, Liang CH, Peng ML, Lee HH (1998) Conjunctivitis caused by Thelazia callipaeda infestation in a woman. J Formos Med Assoc 97: 425–427.
- 54. Bhaibulaya M, Prasertsilpa S, Vajrasthira S (1970) Thelazia callipaeda Railliet and Henry, 1910, in man and dog in Thailand. Am J Trop Med Hyg 19: 476–479.
- 55. Yospaiboon Y, Sithithavorn P, Maleewong V, Ukosanakarn U, Bhaibulaya M (1989) Ocular thelaziasis in Thailand: a case report. J Med Assoc Thai 72: 469–473.
- 56. Otranto D, Ferroglio E, Lia RP, Traversa D, Rossi L (2003) Current status and epidemiological observation of Thelazia callipaeda (Spirurida, Thelaziidae) in dogs, cats and foxes in Italy: a “coincidence” or a parasitic disease of the Old Continent? Vet Parasitol 116: 315–325.
- 57. Dorchies P, Chaudieu G, Siméon LA, Cazalot G, Cantacessi C, et al. (2007) First reports of autochthonous eyeworm infection by Thelazia callipaeda (Spirurida, Thelaziidae) in dogs and cat from France. Vet Parasitol 149: 294–297.
- 58. Ruytoor P, Déan E, Pennant O, Dorchies P, Chermette R, et al. (2010) Ocular thelaziosis in dogs, France. Emerg Infect Dis 16: 1943–1945.
- 59. Vieira L, Rodrigues FT, Costa A, Diz-Lopes D, Machado J, et al. (2012) First report of canine ocular thelaziosis by Thelazia callipaeda in Portugal. Parasit Vectors 5: 124.
- 60. Miró G, Montoya A, Hernández L, Dado D, Vázquez MV (2011) Thelazia callipaeda: infection in dogs: a new parasite for Spain. Parasit Vectors 4: 148.
- 61. Malacrida F, Hegglin D, Bacciarini L, Otranto D, Nägeli F, et al. (2008) Emergence of canine ocular Thelaziosis caused by Thelazia callipaeda in southern Switzerland. Vet Parasitol 157: 321–327.
- 62. Gasser RB, LeGoff L, Petit G, Bain O (1996) Rapid delineation of closely-related filarial parasites using genetic markers in spacer rDNA. Acta Trop 62: 143–150.
- 63. Liu GH, Wang Y, Xu MJ, Zhou DH, Ye YG, et al. (2012) Characterization of the complete mitochondrial genomes of two whipworms Trichuris ovis and Trichuris discolor (Nematoda: Trichuridae). Infect Genet Evol 12: 1635–1641.
- 64.
De Ley P, Blaxter M (2002) Systematic position and phylogeny. In: Lee DL (ed). The Biology of Nematodes. Taylor & Francis, London and New York. pp. 1–30.
- 65. Nadler SA, Carreno RA, Mejía-Madrid H, Ullberg J, Pagan C, et al. (2007) Molecular phylogeny of clade III nematodes reveals multiple origins of tissue parasitism. Parasitology 134: 1421–1442.