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. (Continued on next page) * Correspondence: jkpyou@chungbuk.ac.kr Graduate Program in Cell Biology and Genetics, College of Medicine, Chungbuk National University, Cheongju 361-763, Republic of Korea Graduate Program in Cell Biology and Genetics and Department of Parasitology, College of Medicine, Chungbuk National University, Cheongju 361-763, Republic of Korea Full list of author information is available at the end of the article © 2013 Sultana et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Sultana et al. BMC Evolutionary Biology 2013, 13:12 http://www.biomedcentral.com/1471-2148/13/12