Comparative and Evolutionary Analyses of Meloidogyne spp. Based on Mitochondrial Genome Sequences

Molecular taxonomy and evolution of nematodes have been recently the focus of several studies. Mitochondrial sequences were proposed as an alternative for precise identification of Meloidogyne species, to study intraspecific variability and to follow maternal lineages. We characterized the mitochondrial genomes (mtDNAs) of the root knot nematodes M. floridensis, M. hapla and M. incognita. These were AT rich (81–83%) and highly compact, encoding 12 proteins, 2 rRNAs, and 22 tRNAs. Comparisons with published mtDNAs of M. chitwoodi, M. incognita (another strain) and M. graminicola revealed that they share protein and rRNA gene order but differ in the order of tRNAs. The mtDNAs of M. floridensis and M. incognita were strikingly similar (97–100% identity for all coding regions). In contrast, M. floridensis, M. chitwoodi, M. hapla and M. graminicola showed 65–84% nucleotide identity for coding regions. Variable mitochondrial sequences are potentially useful for evolutionary and taxonomic studies. We developed a molecular taxonomic marker by sequencing a highly-variable ~2 kb mitochondrial region, nad5-cox1, from 36 populations of root-knot nematodes to elucidate relationships within the genus Meloidogyne. Isolates of five species formed monophyletic groups and showed little intraspecific variability. We also present a thorough analysis of the mitochondrial region cox2-rrnS. Phylogenies based on either mitochondrial region had good discrimination power but could not discriminate between M. arenaria, M. incognita and M. floridensis.


Introduction
Mitochondria are subcellular organelles derived from an alpha-proteobacterial endosymbiont, wherein oxidative phosphorylation and other important biochemical functions take place [1].Most of the genes of the original bacterial chromosome were lost or transferred to the nucleus, with a handful of genes remaining in animal mitochondrial genomes (mtDNAs) [2].Animal mtDNAs share several features, including relatively constant gene content and order, maternal inheritance, reduced recombination rate, and high evolutionary rate.This molecule is central in disease, apoptosis and ageing studies, in genetic and genomic analyses, and in developing genetic markers for molecular systematics [3].Animal mtDNAs are generally small and circular molecules containing 12-13 protein-coding genes for enzymes of the oxidative phosphorylation pathway, 22 transfer RNA-coding genes (tRNAs), and 2 ribosomal RNA-coding genes [2].This gene content is variable only in a few groups of metazoans, such as some nematodes that lack the gene atp8 [3], a bivalve that also lacks atp8 and has an extra tRNA gene [4], and cnidarians that have lost nearly all tRNAs genes [5].Usually, no introns are present in animal mtDNAs and intergenic sequences are small, although one large non-coding region with control elements is found in some metazoan species [6].To date, more than 100 complete mtDNA sequences of the phylum Nematoda have been published.A phylum-wide comparison revealed that nematode mtDNAs are characterized by a very rapid rate of mitochondrial gene rearrangement and are smaller in size than those of other metazoan groups, varying in size from 13 kb to 15 kb with some exceptions as Romanomermis culicivorax (26 kb), Pratylenchus vulnus (21 kb) Hexamermis agrotis (24 kb) [2,3,7,8,9].
Root-knot nematodes (genus Meloidogyne, infraorder Tylenchomorpha, suborder Tylenchina, order Rhabditida, class Chromadorea) [10] are plant parasites, comprising more than 100 species, in addition to several described races [11] and rank first based on a list of scientifically and economically important plant parasitic nematodes [12].Their common name refers to typical galls on the roots of host plants that reduce the uptake of water and nutrients resulting in lower crop productivity [13].The genus Meloidogyne includes the most widespread and economically damaging nematodes worldwide: M. incognita, M. javanica, M. arenaria, M. hapla, M. chitwoodi and M. enterolobii [14].Accurate identification of Meloidogyne species is critical for effective crop management.Traditionally, species diagnosis was done based on morphological characters, host range and esterase isozyme electrophoresis and were later combined with DNA-based methods [15].These molecular techniques were initially based on analyses of DNA fragments and sequencing of the nuclear rDNA operon and became obsolete or less effective with the increasing description of new species [16].Mitochondrial sequences were proposed as an alternative for precise identification of Meloidogyne species, to study intraspecific variability and to follow maternal lineages [17].Mitochondrial markers have been used in the past, but with a different fragment than in this current study [17,18,19,20].
To gain insight into the evolutionary patterns, sequence diversity and potential taxonomic use of the mitochondrial genome of the genus Meloidogyne, we characterized the mtDNAs of three root-knot nematodes: a facultative meiotic parthenogenetic strain of M. hapla and M. floridensis and a strain of the obligate mitotic parthenogenetic M. incognita, taking advantage of genomic projects that focused on these species.We undertook detailed genomic comparisons among them and with mtDNAs of Meloidogyne spp. recently published [21,22,23].Based on genomic alignments of Meloidogyne, we identified a 2 kb variable mitochondrial marker and performed phylogenetic analyses to aid in Meloidogyne species identification and evolution and to elucidate relationships within the genus.

Mitochondrial genome annotation
The nuclear genomes of three species of root-knot nematodes, M. floridensis [24], M. hapla [25], and M. incognita [26] have been reported and mitochondrial sequences were made available in public sequence databases as part of these genome sequencing projects but have not been yet assembled or analyzed.Sequence similarity searches were performed with BLASTn [27] against the GenBank databases nr, WGS (whole-genome shotgun), and EST (expressed sequence tag), restricted to the genus Meloidogyne without filtering low complexity regions and using the mtDNA of several Nematoda as the query sequence.A total of 17 contigs of the rootknot nematodes M. floridensis, M. hapla and M. incognita had similarity to mitochondrial genes (S1 Table ).Mitochondrial fragments of each species were assembled based on overlapping sequences and by amplification and sequencing of one joining fragment.We designed primers to link two mitochondrial contigs of M. hapla ABLG01002800 and ABLG01002664 (S1 Table ): Mh32800R (5`-AGAGTGAATTGGTAAGAGG-3`) and Mh12664F (5`-CGGTAA CCAAAAACCTCCAAGC-3`).We amplified and sequenced a product of 1500 bp using DNA of M. hapla 7J2 (S2 Table ).It is not possible to rule out the possibility that some of these mitochondrial sequences belong to NUMTs in the nuclear genome.However, the lack of nuclear sequences in the fragments containing mitochondrial genes suggests that those DNA fragments reside in the mitochondria.
Open reading frames of the mitochondrial sequences of M. floridensis, M. hapla and M. incognita were recognized with Sequencher 5.2.2 (Genes Codes Corporation) using the invertebrate mitochondrial genetic code.Each putative protein-coding gene was used as query for BLASTx searches against GenBank protein databases to identify them.Ribosomal RNA genes (rRNAs) were recognized by comparison to other rRNAs of nematodes using BLASTn [27].The tRNA genes were identified using tRNAscan-SE search server [28] or by similarity to those reported in other Meloidogyne species.Secondary structures of tRNAs were predicted with tRNAscan-SE.
Gene boundaries, i.e. start and stop codons, are generally predicted based on the presence of a canonical start codon followed by a standard stop codon.When non-universal genetic codes or non-canonical stop codons are involved, as reported for animal mtDNAs [3], gene boundaries are less obvious and several start or stop codons may be available.A disparity in the criterion for establishing gene boundaries was observed for several genes in Meloidogyne mtDNA, even between two strains of the same species [21,23].In this study, start codons were set at the first start codon available (including non-universal ones) that had the minimum overlap with an upstream gene.Complete universal stop codons (TAA or TAG) were selected even when they overlapped with a downstream gene.

Analyses of mitochondrial sequences
Pairwise comparisons of protein, rRNA, and tRNA genes of M. chitwoodi, M. floridensis, M. graminicola, M. hapla and M. incognita were done in MEGA6 [29].Codon usage of proteincoding genes and nucleotide composition were calculated with the software SMSv2 (Sequence Manipulation Suite) [30].Whole-genome alignments of Meloidogyne spp.mitochondrial genomes were performed with the VISTA pipeline infrastructure [31].Tandem repetitive elements were identified using Tandem Repeat Finder software [32].

PCR amplification and sequencing of a mitochondrial region
DNA samples of 28 populations of 5 Meloidogyne species were kindly provided by Andrea Skantar (USDA) (S2 Table ).DNA of a M. arenaria isolate from Argentina was extracted as reported previously [33].We amplified and sequenced the mitochondrial region between the genes nad5 and cox1 from different Meloidogyne species and isolates using the following primers: Mmt5 (5`-GGTTTAATTGGTGGTTTTGG-3`), Mmt2rc (5`-TGTCCTCAAACTAAACAACC-3`), Mmt9 (5`-TTGGTTGATTGGTGAAAGC-3`), and Mmt4rc (5`-AAACCACCAATTAAAC CAGG-3`).Primers were designed with the software Primaclade [34] based on a sequence alignment of available Meloidogyne species.PCR conditions included an initial denaturing step at 94°C for 2 min, followed by 35 cycles of 94°C for 30 seconds, 45°C for 30 seconds and 72°C for 1 min; followed by a final extension step at 72°C for 8 min.Amplified fragments were sequenced by Sanger sequencing with Applied Biosystems 3730XL.Sequences were deposited in GenBank (accession numbers KM491188-KM491216).

Phylogenetic analyses
Three mitochondrial nucleotide data sets were constructed and analyzed.One consisted of 8 protein-coding genes (cox1, cox2, cox3, cob, nad2, nad3, nad4 and nad5) from 43 diverse nematode species that were retrieved from GenBank.The 8 protein-coding genes were concatenated in a single nucleotide data set of 4,920 bp.The second data set included sequences of the nad5-cox1 region from 36 populations of 8 species of Meloidogyne (S2 Table ).The third data set consisted of mitochondrial sequences of the cox2-rrnS region from 151 isolates of 20 Meloidogyne species that were retrieved from GenBank.Nucleotide sequences were manually aligned with MacClade 4.07 [35].Alignment regions for which positional homology could not be determined with confidence by visual inspection were excluded.Alignments and trees have been deposited in Treebase (study S17102).
Maximum Likelihood (ML) phylogenetic analyses were performed with Garli 0.951 [36] under the General Time Reversible model with parameters for invariable sites and gamma-distributed rate heterogeneity.The GTR substitution model was chosen by using Modeltest for each data set [37].Ten independent runs were conducted using either the automated stopping criterion or for up to 5,000,000 generations to ensure convergence to a similar topology and likelihood score.A hundred ML bootstrap replicates were performed.Maximum Parsimony (MP) analyses were conducted using PAUP Ã 4.0 [38].The search strategy involved tree-bisection reconnection branch swapping, MulTrees option in effect.To assess node support, MP bootstrap analyses were performed using 1,000-replicate heuristic search.

Results and Discussion
Characterization of the mitochondrial genomes of M. incognita, M. floridensis, and M. hapla The mtDNA of M. incognita assembled as a circular molecule of 17,985-18,332 bp, depending on the number of 63-bp repeats (shown linearized in Fig. 1).The incomplete mitochondrial genomes of Meloidogyne hapla and M. floridensis assembled as single molecules of 17,355 bp and 15,811 bp respectively (Fig. 1, Table 1).The total length of the latter genomes is unknown due to missing data in the long intergenic region between the genes nad4 and atp6 (Fig. 1).Complete mitochondrial genomes of three Meloidogyne spp.have been described with similar genome length (Table 1, Fig. 1).Also, the size of the mtDNA of M. javanica (20.5 kb) has been reported [39].A comparison to other animal mtDNAs revealed that Meloidogyne mtDNAs fall within the size range of Metazoa (8 kb to 48 kb) [6,40].
We identified 22 discrete sequences (ranging from 49 to 59 bp) predicted to fold into secondary structures of tRNAs in the mtDNAs of M. floridensis, M. hapla, and M. incognita, respectively.In each of the species, 20 tRNAs had the predicted T-arm-lacking structure [5] and two putative tRNA genes folded into different secondary structures, such as the genes tRNA-S 1 (UCN) and tRNA-S 2 (AGN) that had a TyC arm and loop but not the DHU arm.
All genes in Meloidogyne mtDNAs were encoded on the same strand, with very few and often short intergenic sequences (Fig. 1; S3 Table ).This asymmetry in transcriptional architecture is typical of chromadoreans, while enopleans and most animal species have circular mitochondrial genomes with coding genes on both strands [42,43].The nucleotide composition of Meloidogyne spp.mtDNAs showed a large bias towards AT (Table 1), similar to that of their nuclear genomes [24].In general, mitochondrial genomes of the phylum Nematoda showed a trend for AT richness, ranging from 66% to 85%.In addition, a bias towards T and G content S a , V a , S b , and V b correspond to tRNA genes annotated differently in two previous studies (S a and V a in GenBank KJ139963 [23]; S b and V b in GenBank HG529223 [21]).Tandem repeats are represented by boxes labeled by the repeat length in kb followed by "R".doi:10.1371/journal.pone.0121142.g001a Numbers in parentheses indicate the copy-number of each repeat.
b The genome size depends on the number of copies of the R63 repeat, which varies from 7 to 13 in the genomic assembly.
c Not determined because the total genome length is unknown.
doi:10.1371/journal.pone.0121142.t001 in the coding strand was observed in the mtDNAs of M. floridensis, M. hapla and M. incognita, as reported for other nematodes [23,44,45,46].Start and stop codons of protein-coding genes in Meloidogyne mtDNAs were often non-universal.Putative start codons for the 12 mitochondrial protein-coding genes in M. hapla, M. incognita and M. floridensis were ATA, ATG, ATT, TTA, TTG; stop codons were TAA and TAG (S3 Table ).Some of these start codons are non-universal and have been previously reported for nematode mitochondrial genomes [21,22,23,47,48,49].In addition, the termination codon TGA encoded tryptophan, among other deviations from the universal genetic code in nematode mitochondria [50].

Phylogeny of Nematoda based on mitochondrial sequences
We performed a Maximum Likelihood phylogenetic analysis based on nucleotide sequences of eight mitochondrial protein-coding genes (4,920 bp) from 43 species of the phylum Nematoda, including 5 Meloidogyne spp.and 8 taxa of the class Enoplea as outgroups (S1 Fig. ).Mitochondrial genes could resolve most phylogenetic relationships within the phylum Nematoda and these were consistent with previous studies based on mitochondrial [7,8,22] and nuclear genes [51,52].Meloidogyne spp.grouped together as sister to Pratylenchus vulnus within the infraorder Tylenchomorpha with high bootstrap support (S1 Fig. ).Within the clade of Meloidogyne, M. incognita + M. floridensis were sister to M. hapla with great bootstrap support (BS = 100%), and these were sister to M. chitwoodi + M. graminicola.These findings agree with previous studies based on the nuclear gene 18S rRNA [53].

Gene order in nematode mitochondrial genomes
Protein and rRNA gene order was identical within the genus Meloidogyne (Fig. 1), but different to other members of the infraorder Tylenchomorpha and any other nematode (S1 Fig. ).Furthermore, the location of the tRNA genes was unique for four species of Meloidogyne, except for M. incognita and M. floridensis that were identical in gene order (Fig. 1).Differences occurred mainly in two blocks of tRNAs; those located between rrnS and nad1 and the group found between the genes cox3 and nad6.
Identification of tRNA genes is controversial in nematode mitochondrial genomes given that tRNAs are short, divergent and have atypical structures [54].As a result, we noticed discrepancies in the prediction of tRNA genes in published Meloidogyne mtDNAs.For example, there was no consensus between the tRNA annotations of the highly similar (see below) mitochondrial genomes of two strains of M. graminicola (asterisks in Fig. 1) using two different tRNAs prediction programs [21,23].
Nematode mtDNAs, including those of root-knot nematodes, showed mitochondrial evolutionary dynamics that differ from most metazoan [55].For example, mitochondrial genomic rearrangements within the phylum Nematoda were significantly higher than those of other animal groups [3].A previous study found 25 different gene rearrangements (GA) in 62 mitochondrial genomes of nematodes, considering the position of tRNAs; of them, the most common was GA3, which was shared by 32 nematode species [8].Within the infraorder Tylenchomorpha, several mtDNAs have recently become available.Heterodera glycines and Pratylenchus vulnus [7,46] showed two different gene arrangements to those already described (indicated as GA26 and GA27 in S1 Fig .).The mtDNAs of Meloidogyne spp.added four novel gene rearrangements to the list of described GA in Nematoda (GA28-GA31; S1 Fig. ).The genus Meloidogyne (and the infraorder Tylenchomorpha as a whole) are particularly variable in gene order in comparison to other lineages of Nematoda, along with the highly re-arranged enopleans [43,56], and the genus Onchocerca [49,57].
The cox2-rrnS region has been largely studied and used in molecular systematics [17,18,19,58,59] because it has indels and polymorphisms useful to differentiate species of the genus Meloidogyne.Until recently, it was possible to identify the major species of Meloidogyne by PCR amplification and digestion of this fragment, or by DNA sequencing of this region.Lately, the analysis of the cox2-rrnS in additional Meloidogyne species uncovered limitations in the taxonomic use of this marker (see below).
The largest intergenic region in all Meloidogyne spp.mtDNAs was found between the protein-coding genes nad4 and atp6 and included 2-5 tRNA-coding genes (Fig. 1).This region also contained tandem repeats in most species, but the repeats were not homologous with a few exceptions (Fig. 1).M. graminicola had a 111-bp (111R) tandem repeat that was similar to the 102-bp (102R) tandem repeat of M. incognita, but different to the 111-bp (111R) tandem repeat of M. chitwoodi.Also, M. chitwoodi had a 48-bp (48R) tandem repeat that aligned with the 94-bp (94R) repeat of M. graminicola (Fig. 1).

Variability of mitochondrial genomes within a single Meloidogyne species
The mtDNAs of 2 strains of M. incognita (this study and [15]) and 2 strains of M. graminicola [21,23] are now available for intraspecific comparisons.Variation in the copy-number of tandem repeats and the presence of indels in intergenic regions are responsible for the differences in genome length within each species (Table 1).A pairwise comparison between the mtDNA of the 2 strains of M. incognita and M. graminicola, respectively, indicated that they shared the same gene content and order (not considering the discrepancy in tRNA annotation of the M. graminicola mtDNAs shown in Fig. 1 given that the putative tRNAs were present in both strains).
Mitochondrial sequences between the strains of M. incognita were extremely similar showing an identity of 99.2% (excluding gaps) with only a few substitutions and indels across coding and non-coding regions (S2-A Fig. ).A pairwise comparison between the mtDNAs of both strains of M. graminicola showed that they were almost identical (99.9% identity, excluding gaps) with very few polymorphisms in the coding regions and most substitutions clustered within the trnD-trnV non-coding region (S2-B Fig. ).The closer geographic distance between the M. graminicola strains (China and Philippines) in comparison to the two isolates of M. incognita (Mexico and France) could relate to the higher mtDNA similarity observed within M. graminicola.It is also possible that the mitochondrial sequence of M. incognita strain Morelos were contaminated with NUMTs, showing more substitutions than the actual mtDNA residing in the mitochondria [60,61].
Variability of mitochondrial genomes among Meloidogyne spp.
The mtDNA of M. floridensis was compared to that of other Meloidogyne spp. with the program VISTA (Fig. 2).M. floridensis and M. incognita mtDNAs were strikingly similar across the genome with an identity of 97-100% for all coding regions.The genes cox3, nad6 and nad3, along with 19 tRNAs, were identical in both species while the other genes showed 1-7 substitutions/gene.This result is in agreement with recent evidence showing that M. incognita is a double hybrid of M. floridensis and another species [24].Comparisons to M. hapla, M. graminicola and M. chitwoodi revealed the presence of conserved genic regions and highly divergent (including some non-homologous) intergenic sequences (Fig. 2).The same is true for pairwise comparisons between the other Meloidogyne species.

Analysis of mitochondrial regions for taxonomic use
Variable intergenic mitochondrial sequences are potentially useful for phylogenetic studies and species identification.Based on mitochondrial genomic comparisons (i.e.Fig. 2), we identified a highly variable ~2 kb mitochondrial segment starting on the gene nad5 through the middle of the gene cox1, including 6 tRNAs and intergenic regions.We selected this region because it showed interspecific variability and it was easy to amplify by PCR from a wide range of species given the presence of highly conserved flanking sequences.The nad5-cox1 region was sequenced and analyzed in 36 populations of 8 Meloidogyne spp.(S2 Table ).The overall topologies of the ML and MP phylogenetic trees were almost identical (Fig. 3).
Isolates of each species formed monophyletic groups with high bootstrap support (95-100%), except for populations of M. arenaria, M. incognita and M. floridensis that grouped in a single clade but were not separated by species (Fig. 3).The few polymorphisms observed in this heterogeneous group fell within the range of intraspecific variability.Within each of the monophyletic species in the tree, sequences were highly conserved with very few or no polymorphisms.This limited intraspecific variability (Fig. 3) has been reported for parthenogenetic nematodes [17,62] and it is now extended to sexual species (Fig. 3).The mitochondrial genomic comparisons between the 2 strains of M. incognita and M. graminicola (S2 Fig. ) also agreed with this statement.Phylogenetic relationships observed among the sampled species (Fig. 3) were congruent with previous phylogenies based on the nuclear 18S rRNA [63] or the mitochondrial region cox2-rrnS.
Besides the mitochondrial marker nad5-cox1, a single mitochondrial region, cox2-rrnS, has been extensively analyzed within the genus Meloidogyne [17,18,19,20].Phylogenetic trees based on this marker had good discriminating power but relationships among several species of root-knot nematodes were unresolved and several isolates could not be identified at the species level [58,64,65,66].A phylogeny including all cox2-rrnS sequences available in public databases for Meloidogyne species (S3 Fig. ) showed that M. graminis, M. marylandi, M. hapla, M. partityla, M. haplanaria, M. enterolobii, M. naasi and M. graminicola could be separated based on a number of shared substitutions that took place since each species diverged.In contrast, some species had highly similar cox2-rrnS sequences and could not be differentiated.These and previous studies based on nuclear or mitochondrial data [53,58,63,67] found that mitotic parthenogenetic Meloidogyne species (shown in red in Fig. 3  In addition to the results found with both mitochondrial markers, other studies also indicated that mitotic parthenogenetic species and M. floridensis are hard to distinguish using the nuclear regions ITS and 18S rDNA [33,53,65,68] and occasionally showed contradictory results between nuclear and mitochondrial phylogenies [17,58].Species boundaries, especially for mitotic parthenogenetic nematodes, have been discussed [17] and interspecific hybridizations prior to the recent establishment of parthenogenesis were proposed [17,69].A comparative analysis of gene families based on the complete nuclear genomes of M. hapla [25], M. incognita [26] and M. floridensis [24] suggested that M. incognita is a double-hybrid between M. floridensis and an unidentified species [24].
It remains challenging to distinguish most mitotic parthenogenetic Meloidogyne isolates and to identify the species, if they, indeed, belong to different phylogenetic species.Extensive studies that include a wide range of populations and species of Meloidogyne while comparing mitochondrial and nuclear sequences are needed to elucidate the complex evolution of this group.Single copy nuclear genes may aid in the identification of mitotic parthenogenetic species.

Fig 1 .
Fig 1. Linear maps of the mitochondrial genomes of five Meloidogyne species.The maps were drawn to scale.Dots interrupting the main genomic line indicate unknown sequences.All genes are encoded on the same strand of the mtDNA.Protein and rRNA genes have standard nomenclature.tRNA genes are designated by single-letter abbreviation.Two tRNAs for Leucine (L) and Serine (S) are present: L 1 (anticodon UAG), L 2 (UAA), S 1 (UCN) and S 2 (AGN).S a , V a , S b , and V b correspond to tRNA genes annotated differently in two previous studies (S a and V a in GenBank KJ139963[23]; S b and V b in GenBank HG529223[21]).Tandem repeats are represented by boxes labeled by the repeat length in kb followed by "R".

Fig 2 .
Fig 2. Pairwise comparisons between the mitochondrial genome of M. floridensis (this study) and that of M. incognita (KJ47615), M. hapla (this study), M. chitwoodi (KJ476150) and M. graminicola (KJ139963) using the program VISTA.The Y-axis represents the sequence identity (50-100%).Protein (in violet) and rRNA (in yellow) genes are indicated in the standard nomenclature.tRNA genes (in red) are designated by a single-letter abbreviation.The range of identity (%) of protein-coding genes among Meloidogyne spp. is shown at the top.doi:10.1371/journal.pone.0121142.g002

Fig 3 .
Fig 3. Maximum Likelihood phylogenetic tree of Meloidogyne spp.based on the nad5-cox1 mitochondrial region.The tree was rooted according to Holterman et al. [53].Numbers above and below branches represent ML and MP bootstrap support values >50%, respectively.Species names in red and green indicate mitotic and meiotic parthenogenetic species, respectively.doi:10.1371/journal.pone.0121142.g003 For example, M. fallax was paraphyletic with respect to M. chitwoodi with almost identical sequences (S3 Fig.).The greatest unresolved clade included the major and most economically damaging Meloidogyne spp, along with other minor root-knot nematodes (S3 Fig.).These 10 species formed a monophyletic clade with highly similar cox2-rrnS sequences and indels shared by different taxa, which prevented a clear distinction of species based on this marker.
and S3 Fig.) were evolutionary distinct from the meiotic parthenogenetic ones (shown in green in Fig. 3 and S3 Fig.), with a few exceptions such as M. floridensis that is meiotic parthenogenetic and clusters within the first group (Fig. 3; S3 Fig.).

Table 1 .
Comparisons between sequenced mitochondrial genomes of the genus Meloidogyne.