Comparative analysis of complete mitochondrial genome sequences confirms independent origins of plant-parasitic nematodes

Background: The nematode infraorder Tylenchomorpha (Class Chromadorea) includes plant parasites that are of agricultural and economic importance, as well as insect-associates and fungal feeding species. Among tylenchomorph plant parasites, members of the superfamily Tylenchoidea, such as root-knot nematodes, have great impact on agriculture. Of the five superfamilies within Tylenchomorpha, one (Aphelenchoidea) includes mainly fungal-feeding species, but also some damaging plant pathogens, including certain Bursaphelenchus spp. The evolutionary relationships of tylenchoid and aphelenchoid nematodes have been disputed based on classical morphological features and molecular data. For example, similarities in the structure of the stomatostylet suggested a common evolutionary origin. In contrast, phylogenetic hypotheses based on nuclear SSU ribosomal DNA sequences have revealed paraphyly of Aphelenchoidea, with, for example, fungal-feeding Aphelenchus spp. within Tylenchomorpha, but Bursaphelenchus and Aphelenchoides spp. more closely related to infraorder Panagrolaimomorpha. We investigated phylogenetic relationships of plant-parasitic tylenchoid and aphelenchoid species in the context of other chromadorean nematodes based on comparative analysis of complete mitochondrial genome data, including two newly sequenced genomes from Bursaphelenchus xylophilus (Aphelenchoidea) and Pratylenchus vulnus (Tylenchoidea). Results: The complete mitochondrial genomes of B. xylophilus and P. vulnus are 14,778 bp and 21,656 bp, respectively, and identical to all other chromadorean nematode mtDNAs in that they contain 36 genes (lacking atp8) encoded in the same direction. Their mitochondrial protein-coding genes are biased toward use of amino acids encoded by T-rich codons, resulting in high A+T richness. Phylogenetic analyses of both nucleotide and amino acid sequence datasets using maximum likelihood and Bayesian methods did not support B. xylophilus as most closely related to Tylenchomorpha (Tylenchoidea). Instead, B. xylophilus, was nested within a strongly supported clade consisting of species from infraorders Rhabditomorpha, Panagrolaimomorpha, Diplogasteromorpha, and Ascaridomorpha. The clade containing sampled Tylenchoidea (P. vulnus, H. glycines, and R. similis) was sister to all analyzed chromadoreans. Comparison of gene arrangement data was also consistent with the phylogenetic relationships as inferred from sequence data. Alternative tree topologies depicting a monophyletic grouping of B. xylophilus (Aphelenchoidea) plus Tylenchoidea, Tylenchoidea plus Diplogasteromorpha (Pristionchus pacificus), or B. xylophilus plus Diplogasteromorpha were significantly worse interpretations of the mtDNA data. Conclusions: Phylogenetic trees inferred from nucleotide and amino acid sequences of mtDNA coding genes are in agreement that B. xylophilus (the single representative of Aphelenchoidea) is not closely related to Tylenchoidea, indicating that these two groups of plant parasites do not share an exclusive most recent common ancestor, and that certain morphological similarities between these stylet-bearing nematodes must result from convergent evolution. In addition, the exceptionally large mtDNA genome size of P. vulnus, which is the largest among chromadorean nematode mtDNAs sequenced to date, results from lengthy repeated segments in non-coding regions.