- Research article
The complete mitochondrial genomes of three parasitic nematodes of birds: a unique gene order and insights into nematode phylogeny
BMC Genomicsvolume 14, Article number: 414 (2013)
Analyses of mitochondrial (mt) genome sequences in recent years challenge the current working hypothesis of Nematoda phylogeny proposed from morphology, ecology and nuclear small subunit rRNA gene sequences, and raise the need to sequence additional mt genomes for a broad range of nematode lineages.
We sequenced the complete mt genomes of three Ascaridia species (family Ascaridiidae) that infest chickens, pigeons and parrots, respectively. These three Ascaridia species have an identical arrangement of mt genes to each other but differ substantially from other nematodes. Phylogenetic analyses of the mt genome sequences of the Ascaridia species, together with 62 other nematode species, support the monophylies of seven high-level taxa of the phylum Nematoda: 1) the subclass Dorylaimia; 2) the orders Rhabditida, Trichinellida and Mermithida; 3) the suborder Rhabditina; and 4) the infraorders Spiruromorpha and Oxyuridomorpha. Analyses of mt genome sequences, however, reject the monophylies of the suborders Spirurina and Tylenchina, and the infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. Monophyly of the infraorder Ascaridomorpha varies depending on the methods of phylogenetic analysis. The Ascaridomorpha was more closely related to the infraorders Rhabditomorpha and Diplogasteromorpha (suborder Rhabditina) than they were to the other two infraorders of the Spirurina: Oxyuridorpha and Spiruromorpha. The closer relationship among Ascaridomorpha, Rhabditomorpha and Diplogasteromorpha was also supported by a shared common pattern of mitochondrial gene arrangement.
Analyses of mitochondrial genome sequences and gene arrangement has provided novel insights into the phylogenetic relationships among several major lineages of nematodes. Many lineages of nematodes, however, are underrepresented or not represented in these analyses. Expanding taxon sampling is necessary for future phylogenetic studies of nematodes with mt genome sequences.
Nematodes (also called roundworms) are an extremely diverse group of bilateral animals with an estimate of 1–10 million species although only ~25,000 species have been described . Nematode taxonomy has traditionally been formed from morphology and ecology, e.g. the early systems proposed by Schneider (1866)  on somatic musculature and by Cobb (1919)  on stoma armature, and the system proposed by Filipjev (1929)  on the presence and absence of zooparasitism. Chitwood (1937)  proposed initially the bipartite phylogenetic and taxonomic system based on the presence and absence of phasmids; this system was further elaborated later by Maggenti (1963, 1983) [6, 7] based on pharyngeal structure and excretory systems. However, due to the lack of homologous characters and informative fossil records, and the extensive convergent evolution, it has been extremely difficult to derive a consistent phylogenetic and taxonomic framework for the phylum Nematoda from morphological and ecological characters [8–12]. In the past two decades, nematode phylogeny and taxonomy have been revised with analyses of the nuclear small subunit (SSU) rRNA gene sequences [12–14]. The current working hypothesis of Nematoda phylogeny incorporated evidence from morphology, ecology and SSU rRNA gene sequence analyses [13, 14].
The mitochondrial (mt) genomes of bilateral animals usually contain 37 genes and a control region on a circular chromosome, ~16 kb in size [15–17]. Sequences of individual mt genes and whole mt genomes are widely used to infer phylogenetic relationships among animals at different taxonomic levels [18, 19]. In a number of cases, arrangement of mt genes has also been used to resolve long-standing phylogenetic relationships that could not be resolved by other means [20–22].
Recently, mt genome sequences have also been analyzed to understand the phylogenetic relationships among nematodes. Kang et al.  inferred the Nematoda phylogeny with the mt genome sequences of 25 species; Park et al.  and Sultana et al.  expanded further the analysis to include 36 and 41 species respectively. A major difference between mt genome phylogenies and the current working hypothesis of the Nematoda is on the relationship among four infraorders: Ascaridomorpha, Rhabditomorpha, Spiruromorpha and Oxyuridomorpha. In the current working hypothesis, Ascaridomorpha, Spiruromorpha and Oxyuridomorpha are in the suborder Spirurina whereas Rhabditomorpha is in the suborder Rhabditina [13, 14]. Furthermore, Ascaridomorpha is most closely related to Spiruromorpha, whereas Oxyuridomorpha is sister to Ascaridomorpha + Spiruromorpha + Rhigonematomorph. In the phylogenies inferred from mt genome sequences, however, Ascaridomorpha is most closely related to Rhabditomorpha; in most analyses, Oxyuridomorpha is sister to Ascaridomorpha + Rhabditomorpha, whereas Spiruromorpha is sister to the group that contains Ascaridomorpha, Rhabditomorpha and Oxyuridomorpha [23, 24]. Due to limited taxon sampling, however, the novel phylogenetic relationships inferred from mt genome sequences in these studies need to be interpreted with caution and further test with more taxa from a wide range of lineages is necessary.
In this study, we sequenced the complete mt genomes of three Ascaridia species in the family Ascaridiidae of the infraorder Ascaridomorpha. Ascaridia species are among the most prevalent and pathogenic parasitic nematodes found in domestic and wild birds and have a worldwide distribution . These Ascaridia species share an identical arrangement of mt genes to each other but differ substantially from those of other nematodes. We inferred the phylogenetic relationships with the complete mt genome sequences of the Ascaridia species and 62 other nematode species that have been sequenced to date. Our analyses support the close relationship between Ascaridomorpha and Rhabditomorpha, and provide novel insights into other phylogenetic relationships in the Nematoda.
Mitochondrial genomes of three Ascaridiaspecies
The complete mt genomes of A. galli, A. columbae and Ascaridia sp. (GHL-2012) were 13,977 bp, 13,931 bp and 13,862 bp long, respectively (GenBank accession numbers JX624728, JX624729 and JX624730). Each mt genome contains 12 protein-coding genes, 22 tRNA genes and two rRNA genes on a circular chromosome (Figure 1). As in most other nematodes, atp8 gene is not present in the mt genome of A. galli, A. columbae and Ascaridia sp. (GHL-2012); all mt genes are transcribed from the same direction. TTT is used as the initiation codon in nad6 gene in the mt genome of A. galli. TTT was previously reported to be the initiation codon for cox2, cytb, nad4, nad3, nad2 and cox3 genes in the mt genome of Strongyloides stercoralis (Strongyloididae), but not for any other Ascaridida nematodes. All of the 12 protein-coding genes have complete termination codons, TAA or TAG. The tRNA genes in mt genomes of A. galli, A. columbae and Ascaridia sp. (GHL-2012) range from 51 to 76 bp. The tRNA-Ser(AGN) and tRNA-Ser(UCN) have truncated secondary structures; these two tRNAs lack a DHU stem but possess a TΨC loop. In all other 20 tRNA genes, the TΨC arm and variable loop are replaced with a TV replacement loop. There are two non-coding regions (NCR) in the mt genomes of A. galli, A. columbae and Ascaridia sp. (GHL-2012) (Table 1). The longer NCR (NCRL) is between trnC gene and trnN gene, and is 610 bp (A. galli), 563 bp (A. columbae) and 566 bp [Ascaridia sp. (GHL-2012)] respectively. The shorter NCR (NCRS) is between nad4 and trnM gene, and is 157 bp (A. galli), 91 bp (A. columbae) and 101 bp [Ascaridia sp. (GHL-2012)] respectively. A NCR at this location, i.e. between nad4 and trnM, is also present in the mt genomes of most other nematodes [27, 28].
Gene rearrangement in the mitochondrial genomes of Ascaridiaspecies
A. galli, A. columbae and Ascaridia sp. (GHL-2012) have an identical gene arrangement in their mt genomes. The arrangement of mt genes in these three species, however, is substantially different from those of other nematodes (Figure 2). The mt genomes of the 62 species of nematodes that have been sequenced to date showed 25 types of gene arrangements (GA1–GA25, Figure 2), of which GA3 is the most common type and has been found in 32 species of nematodes. In comparison to the GA3 type, at least three rearrangement events occurred in the Ascaridia species: trnM was translocated, and a block of 11 genes (trnN, trnY, nad1, atp6, trnK, trnL 2 , trnS 1 , nad2, trnI, trnR and trnQ) swapped position with another block of 12 genes (trnG, cox2, trnH, rrnL, nad3, nad5, trnA, trnP, trnV, nad6, nad4L and trnW) (Figure 3).
Phylogeny of nematodes inferred from mitochondrial genome sequences
Monophylies of the order Rhabditida and the subclass Dorylaimia
Of the 65 species of nematodes included in the phylogenetic analyses in this study, 54 species were from the order Rhabditida of the class Chromadorea, and 11 species were from the subclass Dorylaimia of the class Enoplea. Both Rhabditida and Dorylaimia were monophyletic in all of the trees inferred with Bayesian, ML and MP methods from concatenated amino acid sequences deduced from the sequences of the 12 mt protein-coding genes. The monophyly of the order Rhabditida was strongly supported with a posterior probability (PP) of 1 in Bayesian analysis (Figure 4), a bootstrapping frequency (BF) of 100% in ML analysis (Figure 5), and a BF of 94% in MP analysis (Figure 6). The monophyly of the subclass Dorylaimia was strongly supported in Bayesian and ML analyses (PP = 1, Figure 4; BF = 97, Figure 5), and was moderately supported in MP analysis (BF = 67%, Figure 6).
Phylogenetic relationships within the order Rhabditida
The 54 species of nematodes in the order Rhabditida included in this study were from three suborders: Rhabditina (24 species), Spirurina (26 species), and Tylenchina (4 species). The monophyly of the suborder Rhabditina was strongly supported in Bayesian and ML analyses (PP = 1, Figure 4; BF = 99%, Figure 5), and was moderately supported in MP analysis (BF = 79%, Figure 6). The Spirurina and the Tylenchina, however, were not monophyletic in all of the three phylogenetic analyses in this study.
Of the 24 species from the suborder Rhabditina, 23 species were from the infraorder Rhabditomorpha and one species from the infraorder Diplogasteromorpha. Species of these two infraorders were most closely related in all of the phylogenetic analyses in this study. The Rhabditomorpha, however, was paraphyletic with respect to the Diplogasteromorpha. Three species from the superfamily Rhabitoidea of the Rhabditomorpha were more closely related to Pristionchus pacificus (Diplogasteromorpha) than they were to the other 20 species from the Rhabditomorpha. The close relationship between the species of the superfamily Rhabitoidea and P. pacificus was strongly supported in Bayesian analysis (PP = 1, Figure 4), and was moderately supported in ML and MP analyses (BF = 83%, Figure 5; BF = 82%, Figure 6). In addition to the Rhabitoidea, three other superfamilies of the infraorder Rhabditomorpha were also represented in our analyses: Ancylostomatoidea (4 species), Strongyloidea (10 species), and Trichostrongyloidea (6 species). The Trichostrongyloidea was monophyletic with strong support in all of the three phylogenetic analyses in this study (PP = 1, Figure 4; BF = 88%, Figure 5; BF = 88%, Figure 6). The Rhabitoidea was monophyletic with strong support in Bayesian analysis (PP = 0.99) and weak support in MP analysis (BF = 52%, Figure 6), but was paraphyletic in ML analysis with weak support (BF = 51%, Figure 5). The Ancylostomatoidea and the Strongyloidea were paraphyletic in all of the three phylogenetic analyses in this study; there was strong support for the paraphylies of these two superfamilies in Bayesian and MP analyses (PP > 0.96, Figure 4; BF = 99%, Figure 5).
The 26 species of the suborder Spirurina included in this study were from three infraorders: Ascaridomorpha (14 species), Oxyuridomorpha (2 species) and Spiruromorpha (10 species). The Oxyuridomorpha and the Spiruromorpha were both monophyletic with strong support in all of the three phylogenetic analyses (PP = 1, Figure 4; BF = 100%, Figure 5; BF = 100%, Figure 6). The Ascaridomorpha, however, was monophyletic only in ML analysis with moderate support (BF = 84%, Figure 5), and was paraphyletic in Bayesian analysis with strong support (PP = 0.99, Figure 4) and MP analysis with weak support (BF = 52%, Figure 6). In both Bayesian and MP analyses, the 10 species of the superfamily Ascarodoidea of the infraorder Ascaridomorpha were more closely related to those of the Rhabditina and Steinernema carpocapsae (Tylenchina) than they were to the three Ascaridia species and Cucullanus robustus, which were also from the infraorder Ascaridomorpha. The three Ascaridia species, for which mt genomes were sequenced in this study, form a monoplyletic group with strong support in all of the three phylogenetic analyses (Figures 4, 5, 6). Cucullanus robustus was sister to Ascaridomorpha + Rhabditina + Steinernema carpocapsae in Bayesian analysis with strong support (PP = 1, Figure 4). Cucullanus robustus, however, was sister the Ascaridia species with moderate support in ML analysis (BF = 83%, Figure 5) and weak support in MP analysis (55%, Figure 6).
The four species of the suborder Tylenchina included in this study were from two infraorders: Panagrolaimorpha (2 species), and Tylenchomorpha (2 species). Both of these infraorders were paraphyletic in all of the three analyses (PP > 0.99, Figure 4; BF > 57%, Figure 5; BF > 51%, Figure 6).
Phylogenetic relationships within the subclass Dorylaimia
The 11 species of the subclass Dorylaimia included in this study were from three orders: Mermithida (7 species), Dorylaimida (1 species), and Trichinellida (3 species). The Mermitida and the Trichinellida were both monophyletic with strong support in all of the three analyses (PP = 1, Figure 4; BF = 100%, Figure 5; BF = 100%, Figure 6). Among the three orders, the Mermithida and the Dorylaimida were more closely related than either of them to the Trichocephalida, with strong support in Bayesian analysis (PP = 1, Figure 4) and moderate support in ML analysis (BF = 73%, Figure 5). The Dorylaimida was grouped with the Trichocephalida in MP analysis but the support to this grouping was weak (BF = 40%, Figure 6).
Phylogenetic relationships among nematodes have been revised in the past two decades using sequences of nuclear small subunit (SSU) rRNA gene [12–14]. The current working hypothesis of Nematoda phylogeny incorporated evidence from both morphology, ecology and nuclear SSU rRNA gene sequence analyses [14, 15]. Recently, Kang et al.  tested the hypothesis of the Nematoda phylogeny with the mt genome sequences of 25 species. Kang et al.  showed strong support to the monophyly of the class Chromadorea, represented by 16 species from the order Rhabditida; this study, however, could not establish the monophyly of the class Enoplea, represented by nine species from the subclass Dorylaimia. Park et al.  and Sultana et al.  expanded the taxon sampling and inferred phylogenies with mt genome sequences of 36 and 41 species of nematodes respectively. Park et al.  showed strong support for the monophylies of both Chromadorea and Enoplea. Sultana et al. , however, could not establish the monophyly of the class Enoplea in most of their analyses. A major difference between the mt genome phylogenies and the current working hypothesis of the Nematode phylogeny is on the monophyly of the suborder Spirurina, which includes Clade III nematodes proposed in Blaxter et al. . In the current working hypothesis, the suborder Spirurina is monophyletic, which includes the infraorders Ascaridomorpha, Spiruromorpha and Oxyuridomorpha, whereas the infraorder Rhabditomorpha is in the suborder Rhabditina [13, 14]. Further, Ascaridomorpha is most closely related to Spiruromorpha, whereas Oxyuridomorpha is sister to the group that includes Ascaridomorpha + Spiruromorpha. In the phylogenies inferred from mt genome sequences, however, Ascaridomorpha is most closely related to Rhabditomorpha; in most analyses, Oxyuridomorpha is sister to Ascaridomorpha + Rhabditomorpha, whereas Spiruromorpha is sister to the group that contains Ascaridomorpha, Rhabditomorpha and Oxyuridomorpha. Kang et al. , Park et al.  and Sultana et al.  contained much less taxa compared to previous phylogenetic studies on nematodes with nuclear SSU rRNA gene sequences (e.g. >200 taxa in Meldal et al. ). The novel phylogenetic relationship inferred from mt genome sequences, thus, needs to be tested further with more taxa from a wide range of nematode lineages.
In the present study, we sequenced the complete mt genomes of three Ascaridia species from the family Ascaridiidae of the superfamily Heterakoidea, which was not represented in previous studies (Kang et al. ; Park et al. ; Sultana et al. ). Furthermore, we inferred the phylogenetic relationships among 65 nematode species, for which complete mt genomes have been sequenced to date. Our analyses support the division of nematodes into two classes, Chromadorea and Enoplea, represented by the order Rhabditida and the subclass Dorylaimia, respectively, in the present study. Both the Rhabditida and the Dorylaimia were monophyletic with strong support, regardless the methods of phylogenetic analysis used. Within the order Rhabditida, the suborder Rhabditina was monophyletic; the other two suborders, Spirurina and Tylenchina, however, were both paraphyletic in all of the three phylogenetic analyses in this study. Monophyly of the suborder Spirurina was well supported in nuclear SSU rRNA gene sequence analyses [12, 14], but was rejected in mt genome sequence analyses in Kang et al. , Park et al. , Sultana et al.  and the present study. Phylogenetic analyses of mt genome sequences in all of these four studies support strongly a close relationship between Ascaridomorpha and Rhabditomorpha to the exclusion of Oxyuridomorpha and Spiruromorpha. Apparently, additional markers for phylogenetic inference are required to resolve the controversy on the Spirurina between mt genome sequence phylogeny and nuclear SSU rRNA phylogeny.
Gene arrangement in mt genomes provides a source of information for phylogenetic inference, independent from mt genome sequences . It is noteworthy that 10 species of Ascaridomorpha and 20 species of Rhabditomorpha included in our phylogenetic analyses share a common pattern of gene arrangement, GA3, in their mt genomes (Figure 2). GA3 is also present in Pristionchus pacificus, the only species from the infraorder Diplogasteromorpha, which is closely related to the species from the superfamily Rhabitoidea of the Rhabditomorpha (Figures 4, 5, 6). Indeed, GA3 is the most common pattern of mt gene arrangement observed in nematodes to date; nevertheless, this pattern is present only in species of the Ascaridomorpha, the Rhabditomorpha and the Diplogasteromorpha. The possibility that GA3 is ancestral to the phylum Nematoda can be rejected with confidence thanks to the discovery of mt gene arrangement pattern GA23 in three Trichinellida species [29, 30]. GA23 shares obvious similarity with that of many arthropods and would resemble the ancestral pattern of mt gene arrangement of the phylum Nematoda more than any other patterns found in the nematodes. There is no evidence either that GA3 is ancestral to the class Chromadorea. Given the strong support to the close relationship between Ascaridomorpha and Rhabditomorpha + Diplogasteromorpha indicated by mt genome sequence analyses in the present study, the most parsimonious explanation is that GA3 is a shared derived character (i.e. synapomorphy) of the clade Ascaridomorpha + Rhabditomorpha + Diplogasteromorpha. The other six patterns, GA1, GA2, GA4, GA5, GA6 and GA25, observed in either Ascaridomorpha and Rhabditomorpha or in species closely related to these two infraorders, can be derived from GA3 by only a few rearrangement events (Figure 2), e.g. GA1, observed in the three Ascaridia species, can be converted from GA3 by three rearrangement events (see above).
Monophyly of the suborder Tylenchina could not be established previously in nuclear SSU RNA gene sequence analyses . No species from this suborder were included in two previous analyses of mt genome sequences [23, 24]. The present study included four species from the Tylenchina and indicated the paraphyly of this suborder. The positions of the four species of the Tylenchina change depending on the methods of phylogenetic inference used in this study. For instance, Steinernema carpocapsae is sister to the suborder Rhabditina in Bayesian analysis (Figure 4) and MP analysis (Figure 6); it is, however, sister to Ascaridomorpha + Rhabditomorpha + Diplogasteromorpha in ML analysis (Figure 5). Further test with more species from this suborder apparently is needed in order to clarify the phylogenetic position of the suborder Tylenchina. Indeed, a very recent study by Sultana et al.  included six species from the Tylenchina and showed the paraphyly of this suborder with strong support. The paraphyly of the infraorder Rhabditomorpha with respect to the Diplogasteromorpha will also need further test with expanded taxon sampling as only one species, Pristionchus pacificus, from the Diplogasteromorpha was included in the present study. No species from the subclass Enoplia, which has two orders and eight suborders , was included in our analyses. For the subclass Chromadoria, all of the species included in this study were from the order Rhabditida whereas the other five orders, Plectida, Araeolaimida, Monhysterida, Desmodorida, and Chromadorida, were not represented. Expanding taxon sampling from these lineages of nematodes is clearly the next step for phylogenetic studies of nematodes with mt genome sequences.
We sequenced the complete mt genomes of three Ascaridia species (family Ascaridiidae) that infest chickens, pigeons and parrots, respectively. These species have an identical arrangement of mt genes to each other but differ substantially from other nematodes. Phylogenetic analyses of the mt genome sequences of the Ascaridia species, in combination with 62 other nematode species, support the monophylies of the subclass Dorylaimia; the orders Rhabditida, Trichinellida and Mermithida; the suborder Rhabditina; and the infraorders Spiruromorpha and Oxyuridomorpha. Analyses of mt genome sequences, however, reject the monophylies of the suborders Spirurina and Tylenchina, and the infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. Monophyly of the infraorder Ascaridomorpha varies depending on the methods of phylogenetic analysis. The Ascaridomorpha was more closely related to the infraorders Rhabditomorpha and Diplogasteromorpha (suborder Rhabditina) than they were to the other two infraorders of the Spirurina: Oxyuridorpha and Spiruromorpha. The closer relationship among Ascaridomorpha, Rhabditomorpha and Diplogasteromorpha was also supported by a shared common pattern of mitochondrial gene arrangement. Analyses of mitochondrial genome sequences and gene arrangement in the current study and several previous studies have provided novel insights into the phylogenetic relationships among several major lineages of nematodes. Many lineages of nematodes, however, are underrepresented or not represented in these analyses. Expanding taxon sampling is necessary for future phylogenetic studies of nematodes with mt genome sequences.
This study was approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Approval No. LVRIAEC2011-009). The animals, from which the nematodes were collected post-mortem, were handled in strict accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People's Republic of China.
Sample collection and DNA extraction
Adult Ascaridia galli were obtained from naturally infected chickens in Hunan province, China. Adult A. columbae and an undescribed Ascaridia sp. [hereafter Ascaridia sp. (GHL-2012)] were obtained from naturally infected pigeons and parrots in Guangdong province, China. These nematodes were washed in physiological saline, identified preliminarily to species based on host preference, morphological characters and predilection sites [31–33], fixed in 70% ethanol, and stored at −20°C.
Total genomic DNA was isolated from individual nematodes using sodium dodecyl sulphate/proteinase K treatment, followed by spin-column purification (Wizard® SV Genomic DNA Purification System, Promega). In order to independently verify the identity of the specimen, the region spanning ITS-1, 5.8S rRNA gene and ITS-2 was amplified from each individual nematode by PCR using previously reported primers  and sequenced directly.
PCR amplification of the mitochondrial genomes of Ascaridiaspecies
We amplified fragments (400–500 bp) of cox1, nad1 and nad4 genes by PCR using conserved primers (Table 2). The amplicons were sequenced from both directions using BigDye terminator v3.1 on ABI PRISM 3730 platform. We then designed specific primers from these fragments for long PCR amplification (Table 3). We amplified the entire mt genome of individual specimens of A. galli, A. columbae and Ascaridia sp. (GHL-2012) by long PCR in three or four overlapping fragments, respectively. The three overlapping long-PCR fragments for A. galli were between cox1 and nad1 (~2.5 kb), between nad1 and nad4 (~7.5 kb), and between nad4 and cox1 (~5.0 kb). The four overlapping long-PCR fragments for A. columbae and Ascaridia sp. (GHL-2012) were between cox1 and nad1 (~2.5 kb), between nad1 and nad3 (~5 kb), between nad3 and cox2 (~2.5 kb) and between cox2 and cox1 (~5.0 kb). Each long-PCR reaction (25 μl) contains 2.5 μL 2 mM MgCl2, 4.0 μL 0.2 mM each of dNTPs, 2.5 μL 2.5 μl 10× rTaq buffer, 0.25 μL 2.5 μM of each primer, 0.25 μL 1.25 U rTaq polymerase (Takara), and 2 μl of total genomic DNA. The PCR conditions were: 92°C for 2 min (initial denaturation), then 92°C for 10 sec (denaturation), 48–51°C for 30 sec (annealing) and 60°C for 10 min (extension) for 10 cycles, followed by 92°C for 10 sec, 48–51°C for 30 sec, and 60°C for 10 min for 20 cycles and a final extension at 60°C for 10 min. PCR amplicons were column-purified and then sequenced using a primer walking strategy .
Sequence assembly and mitochondrial genome annotation
Sequence reads of A. galli, A. columbae and Ascaridia sp. (GHL-2012) were assembled with ContigExpress program of the Vector NTI software package version 6.0 (Invitrogen, Carlsbad, CA). The mt genome sequences of the Ascaridia species were aligned with those of other nematode species available in GenBank using Clustal X 1.83  to infer gene boundaries. The open-reading frames and codon usage profiles of protein-coding genes were analysed by the Open Reading Frame Finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html) using the invertebrate mitochondrial genetic code, and subsequently compared with those of Ascaris suum and Contracaecum rudolphii B . Sequences of protein-coding genes were translated into amino acid sequences using the invertebrate mitochondrial genetic code in MEGA 5.0 . Translation initiation and termination codons were identified by comparison with those of the nematodes reported previously [27, 28]. The secondary structures of 22 tRNA genes were predicted using tRNAscan-SE  with manual adjustment . Tandem repeats in the non-coding regions were found using Tandem Repeat Finder program (http://tandem.bu.edu/trf/trf.html) .
Considering the high degree of intraspecific variation in nucleotide sequences of mt genes of nematodes , we used the deduced amino acid sequences of mt proteins for phylogenetic analyses. Amino acid sequences were deduced from the sequences of the 12 mt protein-coding genes common for all of the 65 nematode species and two outgroup species included in this study. The deduced amino acid sequences of each of the 12 mt proteins of these species were aligned individually and then were concatenated into a single alignment. Ambiguous sites and regions in the alignment were excluded using Gblocks (http://molevol.cmima.csic.es/castresana/Gblocks_server.html)  with the default parameters. Two Taenia cestodes (T. multiceps and T. hydatigena, GenBank accession number FJ495086 and FJ518620, respectively) were used as the outgroup.
Phylogenetic analyses were conducted using three methods: Bayesian, maximum likelihood (ML), and maximum parsimony (MP). The MtArt 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). Bayesian analysis was with MrBayes 3.1.1 . As MtArt model is not implemented in the current version of MrBayes, an alternative model, MtREV, was used in Bayesian analysis. Four independent Markov chains were run for 1,000,000 metropolis-coupled MCMC generations, sampling a tree every 100 generations. The first 2,500 trees represented burn-in, and the remaining trees were used to calculate Bayesian posterior probabilities (PP). Maximum likelihood (ML) analysis was performed using PhyML 3.0  with the MtArt model. 100 bootstrap replicates were run and bootstrapping frequencies (BF) were calculated. MP analysis was with PAUP* 4.0b10 ; indels were treated as missing characters. A total of 1,000 random addition searches using tree bisection-reconnection (TBR) branch swapping were performed for each MP analysis. Bootstrapping frequency (BF) was calculated from 1,000 bootstrap replicates with 10 random additions per replicate in PAUP. Phylograms were drawn using FigTree v.1.31 (http://tree.bio.ed.ac.uk/software/figtree/).
Dieterich C, Sommer RJ: How to become a parasite - lessons from the genomes of nematodes. Trends Genet. 2009, 25: 203-209. 10.1016/j.tig.2009.03.006.
Schneider A: Monoden. 1866, Berlin: Georg Reimer, 357-
Cobb NA: The orders and classes of nemas. Contrib Sci Nematol. 1919, 8: 213-216.
Filipjev IN: Classification of freeliving Nematoda and their relations to parasitic forms. J Parasitol. 1929, 15: 281-282.
Chitwood BG: “A revised classification of the Nematoda”. Papers on helminthology, 30 year jubileum. Edited by: Skrjabin KJ. 1937, Moscow: All-Union Lenin Academy of Agricultural Sciences, 67-79.
Maggenti AR: Comparative morphology in nemic phylogeny. The Lower Metazoa: Comparative Biology and Phylogeny. Edited by: Dougherty EC, Brown ZN, Hanson ED, Hartman WD. 1963, Berkeley: University of California Press, 273-282.
Maggenti AR: Nematode Higher Classification as Influenced by Species and Family Concepts. 1983, In concepts in nematodesystematics: Academic Press, London
Inglis WG: An outline classification on the Phylum Nematoda. Aust J Zool. 1983, 31: 243-255. 10.1071/ZO9830243.
Anderson RC: Nematode Parasites of Vertebrates: Their Development and Transmission. 1992, Wallingford: CAB International
Lorenzen S: The Phylogenetic Systematics of Free living Nematodes. 1994, London: The Ray Society
Adamson M: Phylogenetic analysis of the higher classification of the Nematoda. Can J Zool. 1987, 65: 1478-1482. 10.1139/z87-230.
Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, Vanfleteren JR, Mackey LY, Dorris M, Frisse LM, Vida JT, Thomas WK: A molecular evolutionary framework for the phylum Nematoda. Nature. 1998, 392: 71-75. 10.1038/32160.
De Ley P, Blaxter M: Systematic position and phylogeny. The Biology of Nematodes. Edited by: Lee DL. 2002, London and New York: Taylor & Francis, 1-30.
Meldal BH, Debenham NJ, De Ley P, De Ley IT, Vanfleteren JR, Vierstraete AR, Bert W, Borgonie G, Moens T, Tyler PA, Austen MC, Blaxter ML, Rogers AD, Lambshead PJ: An improved molecular phylogeny of the Nematoda with special emphasis on marine taxa. Mol Phylogenet Evol. 2007, 42: 622-636. 10.1016/j.ympev.2006.08.025.
Wolstenholme DR: Animal mitochondrial DNA, structure and evolution. Int Rev Cytol. 1992, 141: 173-216.
Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.
Lavrov DV: Key transitions in animal evolution: a mitochondrial DNA perspective. Integr Comp Biol. 2007, 47: 734-743. 10.1093/icb/icm045.
Boore JL, Brown WM: Big trees from little genomes: mitochondrial gene order as a phylogenetic tool. Curr Opin Genet Dev. 1998, 8: 668-674. 10.1016/S0959-437X(98)80035-X.
Shao R, Barker SC: Mitochondrial genomes of parasitic arthropods: implications for studies of population genetics and evolution. Parasitology. 2007, 134: 153-167. 10.1017/S0031182006001429.
Boore JL, Lavrov DV, Brown WM: Gene translocation links insects and crustaceans. Nature. 1998, 392: 667-668. 10.1038/33577.
Boore JL, Brown WM: Mitochondrial genomes of Galathealinum, Helobdella, and Platynereis: sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa. Mol Biol Evol. 2000, 17: 87-106. 10.1093/oxfordjournals.molbev.a026241.
Lavrov DV, Brown WM, Boore JL: Phylogenetic position of the Pentastomida and (pan) crustacean relationships. Proc Biol Sci. 2004, 271: 537-544. 10.1098/rspb.2003.2631.
Kang S, Sultana T, Eom KS, Park YC, Soonthornpong N, Nadler SA, Park JK: The mitochondrial genome sequence of Enterobius vermicularis (Nematoda: Oxyurida)–an idiosyncratic gene order and phylogenetic information for chromadorean nematodes. Gene. 2009, 429: 87-97. 10.1016/j.gene.2008.09.011.
Park JK, Sultana T, Lee SH, Kang S, Kim HK, Min GS, Eom KS, Nadler SA: Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics. 2011, 12: 392-10.1186/1471-2164-12-392.
Sultana T, Kim J, Lee SH, Han H, Kim S, Min GS, Nadler SA, Park JK: Comparative analysis of complete mitochondrial genome sequences confirms independent origins of plant-parasitic nematodes. BMC Evol Biol. 2013, 13: 12-10.1186/1471-2148-13-12.
Abdelqader A, Gauly M, Wollny CB, Abo-Shehada MN: Prevalence and burden of gastrointestinal helminthes among local chickens, in northern Jordan. Prevent Vet Med. 2008, 85: 17-22. 10.1016/j.prevetmed.2008.01.009.
Liu GH, Wu CY, Song HQ, Wei SJ, Xu MJ, Lin RQ, Zhao GH, Huang SY, Zhu XQ: Comparative analyses of the complete mitochondrial genomes of Ascaris lumbricoides and Ascaris suum from humans and pigs. Gene. 2012, 492: 110-116. 10.1016/j.gene.2011.10.043.
Lin RQ, Liu GH, Zhang Y, D’Amelio S, Zhou DH, Yuan ZG, Zou FC, Song HQ, Zhu XQ: Contracaecum rudolphii B: gene content, arrangement and composition of its complete mitochondrial genome compared with Anisakis simplex s.l. Exp Parasitol. 2012, 130: 135-140. 10.1016/j.exppara.2011.11.003.
Liu GH, Gasser RB, Su A, Nejsum P, Peng L, Lin RQ, Li MW, Xu MJ, Zhu XQ: Clear genetic distinctiveness between human- and pig-derived Trichuris based on analyses of mitochondrial datasets. PLoS Negl Trop Dis. 2012, 6: e1539-10.1371/journal.pntd.0001539.
Lavrov DV, Brown WM: 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. 2001, 157: 621-637.
Ramadan HH, Abouznada NY: Morphology and life history of Ascaridia galli in the Domestic Fowl that are raised in Jeddah. J K A U Science. 1992, 4: 87-99.
Wehr EE, Hwang JC: The life cycle and morphology of Ascaridia columbae (Gmelin, 1790) Travassos, 1913 (Nematoda:Ascarididae) in the Domestic Pigeon (Columba livia domestica). J Parasitol. 1964, 50: 131-137. 10.2307/3276047.
Kajerova V, Barus V, Literak I: New records of Ascaridia platycerium (Nematodes) in parrots (Psittaciformes). Vet Med-Czech. 2004, 49: 237-241.
Zhu X, Chilton NB, Jacobs DE, Boes J, Gasser RB: Characterisation of Ascaris from human and pig hosts by nuclear ribosomal DNA sequences. Int J Parasitol. 1999, 29: 469-478. 10.1016/S0020-7519(98)00226-4.
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-2344. 10.1038/nprot.2007.358.
Bowles J, McManus DP: NADH dehydrogenase 1 gene sequences compared for species and strains of the genus Echinococcus. Int J Parasitol. 1993, 23: 969-972. 10.1016/0020-7519(93)90065-7.
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: 4876-4882. 10.1093/nar/25.24.4876.
Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR: The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992, 130: 471-498.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.
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-158. 10.1016/S0020-7519(01)00316-2.
Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.
Hu M, Gasser RB: Mitochondrial genomes of parasitic nematodes progress and perspectives. Trends Parasitol. 2006, 22: 78-84. 10.1016/j.pt.2005.12.003.
Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577. 10.1080/10635150701472164.
Abascal F, Posada D, Zardoya R: MtArt: a new model of amino acid replacement for Arthropoda. Mol Biol Evol. 2007, 24: 1-5.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Swofford DL: Paup*: Phylogenetic analysis using parsimony, version 4.0b10. 2002, Sunderland, MA: Sinauer Associates
This work was supported, in part, by 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. RS acknowledges the funding support from the Australian Research Council (Discovery Project DP120100240) and the Australian Government (Australia-China Science & Research Fund, ACSRF00980).
The authors declare no competing interests.
XQZ, RS and GHL designed the research. GHL and JYL performed the research. DHZ contributed reagents/materials/analyses. GHL, RS and HL analyzed the data. RS, GHL and XQZ wrote the manuscript. All authors have read and approved the final manuscript.
Guo-Hua Liu, Renfu Shao contributed equally to this work.