Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

New evidence of nematode-endosymbiont bacteria coevolution based on one new and one known dagger nematode species of Xiphinema americanum-group (Nematoda, Longidoridae)

  • Mahyar Mobasseri,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Plant Pathology, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran

  • Matthew C. Hutchinson,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Ecology and Evolutionary Biology, Princeton University, Princeton, New Jersey, United States of America

  • Farahnaz Jahanshahi Afshar,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Iranian Research Institute of Plant Protection, Agricultural Research, Education and Extension Organization (AREEO), Tehran, Iran

  • Majid Pedram

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    majid.pedram@modares.ac.ir

    Affiliation Department of Plant Pathology, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran

Abstract

Three populations of Xiphinema primum n. sp. and two populations of X. pachtaicum were recovered from natural forests and cultural regions of northern Iran. Both species belong to the X. americanum-group and were characterized by their morphological, morphometric and molecular data. The new species, which was recovered in three locations, belongs to the X. brevicolle-complex and is characterized by 2124–2981 μm long females with a widely rounded lip region separated from the rest of the body by a depression, 103–125 μm long odontostyle, two equally developed genital branches with endosymbiont bacteria inside the ovary, which are visible under light microscope (LM), vulva located at 51.8–58.0%, the tail is 26–37 μm long with a bluntly rounded end and four juvenile developmental stages. It was morphologically compared with nine similar species viz. X. brevicolle, X. diffusum, X. incognitum, X. himalayense, X. luci, X. parabrevicolle, X. paramonovi, X. parataylori and X. taylori. The second species, X. pachtaicum, was recovered in two geographically distant points close to city of Amol. Molecular phylogenetic studies of the new species were performed using partial sequences of the D2-D3 expansion segments of the large subunit ribosomal RNA gene (LSU rDNA D2-D3), the internal-transcribed spacer rDNA (ITS = ITS1+5.8S+ITS2), and the mitochondrial cytochrome c oxidase I gene (COI mtDNA) regions. The Iranian population of X. pachtaicum was also phylogenetically studied based upon its LSU rDNA D2-D3 sequences. Both species were also inspected for their putative endosymbiont bacteria. Candidatus Xiphinematobacter sp. was detected from two examined populations of the new species, whereas the second endosymbiont bacterium, detected from three examined isolates of X. pachtaicum, was related to the plant and fungal endosymbionts of the family Burkholderiaceae. The phylogenetic analyses of the two endosymbiont bacteria were performed using partial sequences of 16S rDNA. In cophylogenetic analyses, significant levels of cophylogenetic signal were observed using both LSU rDNA D2-D3 and COI mtDNA markers of the host nematodes and 16S rDNA marker of the endosymbiont bacteria.

Introduction

The nematode species of the genus Xiphinema Cobb, 1913 [1] are migratory ectoparasites of roots in a wide range of plants. They cause damage through direct feeding on the root cells of host plants and also through the transmission of some pathogenic plant nepoviruses [2]. Species of this genus are traditionally and historically divided into two groups, those of the americanum-group, and those of non-americanum-group [36].

The X. americanum-group species are generally small nematodes (compared to non-americanum species; usually, 1–3 mm long females) that are unified by the morphological similarities in the lip region, a spiral body after relaxation, and in genital tracts characters (for the revised grouping characters, see Lamberti et al. [5]). Historically, Lima [7] and Tarjan [8] noted the species with these morphological similarities as the X. americanum-complex. Species delimitation has been an issue due to their close morphology and usually overlapping morphometric ranges [9] and thus yielded a “complexity” of their identification, and designation of the X. americanum-complex. Mean values of morphometric data ranges have recently been proposed for primary clustering of populations [10]. The aforementioned “complex”, is nowadays known as X. americanum-group, mainly after exploiting high-throughput molecular techniques for species delimitation [11], which are now routinely applied [1218]. On the other hand, available molecular data is lacking for most species of the group. Lamberti et al. [5, 19] listed 49 species in this group, including some species synonymized by Luc [4]. Since 2004, 10 extra species have been added to that group (Table 1).

thumbnail
Table 1. Xiphinema americanum-group spp. described after 2004 (in chronological order).

https://doi.org/10.1371/journal.pone.0217506.t001

From an agricultural vantage point, some species of the group are regarded as quarantine pests. Indeed, in Europe, the European Economic Union has ratified funding for projects dealing with the identification of virus vector species [5]. Some species are also rated as A1 and A2 quarantine nematodes by the European and Mediterranean Plant Protection Organizations [10, 20]. The species X. pachtaicum (Tulaganov, 1938 [22]) Kirjanova, 1951 [23] seemingly occurs worldwide, as it is common in in Europe [2427] and North America [10]. According to Pedram [28], it is common in Iran.

The ability of these nematodes to retain viable plant viruses in absence of host plants adds to their importance as agriculture pests. Indeed, nine species of genus Xiphinema have been proven to transmit viruses (genus Nepovirus, family Secoviridae [27]), but this could be an underestimate as the capacity of several described Xiphinema species to transmit viruses has not examined. By 2000, seven species of the americanum-group Xiphinemas were known as capable of transmitting plant pathogenic nepoviruses in natural or laboratory conditions [5] and the aforementioned research measures and funding allocation, emphasize the need for reliable species identifications in this group.

Besides the biological ability of X. americanum-group members to transmit viruses, their relationships with endosymbiont Candidatus Xiphinematobacter spp. belonging to subdivision 2 of the Verrucomicrobia is another interesting area of research, already documented and discussed [2933]. The nematode species in the group have prokaryotic endosymbionts and reproduce by thelytokous parthenogenesis; therefore, endosymbiont bacteria are maternally inherited [31]. The bacteria primarily occupy the gut of juveniles, while their final destination isfemale’s ovaries for transovarial vertical transmission through thelytokous parthenogenesis [30, 3233]. Long-term association of the nematode with its endosymbiont and species-specific associations are already hypothesized [34]. Recent studies have emphasized identification [16] and cophylogenetic analyses of these endosymbiont bacteria with their host nematodes [10, 34], to corroborate the long-term association and significant cophylogeny between nematode host and endosymbiont. In one of the aforementioned studies [34], a new endosymbiont, belonging to the plant and fungal endosymbionts of the family Burkholderiaceae Garrity, Bell & Lilburn, 2005 [35] was discovered. The cophylogenetic analyses of this group and the host nematodes also demonstrated significant association between phylogenies including species-specific associations (for a few exceptions in species specificity of the endosymbionts, see [34]). In a more recent study, Howe et al. [36] also showed the congruence between the host and endosymbiont phylogeny using genetic markers of varying evolutionary rates.

Several members of the Xiphinema americanum-group have a global distribution [5, 6, 11, 14]. Currently nine species of this group have been reported from Iran [28]. During our surveys in two Mazandaran and Tehran provinces conducted in 2015–2017, which followed our previous studies on Longidoridae Thorne, 1935 [37] in the country (the history within [38, 39]), several populations of X. americanum-group were recovered. Primary morphological studies revealed three of them could be conspecific and do not fit known species of the group thus representing a new species, described herein as Xiphinema primum n. sp. The other populations belonged to X. pachtaicum. Study with light microscopes revealed endosymbiont bacteria occupying ovaries of several females of both species. Accordingly, the present study aims to (i) describe X. primum n. sp. and characterize the recovered populations of X. pachtaicum using morphological and molecular data, (ii) characterize the endosymbiont bacteria of the two nematode species, and (iii) perform cophylogenetic analyses using selected LSU rDNA D2-D3 and COI mtDNA sequences of host nematodes and 16S rDNA sequences of their endosymbiont bacteria.

Material and methods

Ethics statement

Specific permissions were not needed for collecting the nematodes studied. All soil samples were collected from public regions and natural forests, that were not under protection and no endangered animal species were under protection in those regions. The garden samples were collected after taking the permission from the landowner. There was no need for an animal research ethics committee to approve this research.

Nematode materials

The X. americanum-group species/populations studied in the present research were recovered from soil samples collected from cultivated fields and natural forests in two Mazandaran and Tehran provinces using the flotation method [40] and direct extraction of the nematodes using a series of 20, 60 and 270 mesh sieves (USA standard mesh numbers), equal to 850, 250, 53 μm openings size, respectively.

Morphological observations

The hand-picked specimens were killed with 4% formalin solution, transferred to anhydrous glycerin according to De Grisse [41], and mounted on permanent slides. Photographs were taken using an Olympus DP72 digital camera attached to an Olympus BX51 microscope powered with differential interference contrast (DIC). Drawings of the new species were made using a drawing tube attached to the microscope and were redrawn using CorelDRAW software version 16.

Molecular profiles

For DNA extraction, each nematode was washed with distilled water, surface sterilized using sodium hypochlorite (0.5% v/v) for 5 min, or 5% (v/v) solution of commercial Dettol for 1 min., followed by washing in distilled water three times, transferred to an Eppendorf tube containing 1 μl proteinase K (CinnaGen, Tehran) (10 mg/ml) and 49 μl of extraction buffer (worm lysis buffer [42]; containing: 50 mM KCL, 10 mM Tris-Cl pH 8.3, 2.5 mM MgCl2, 0.45% NP40, and 0.45% Tween 20), frozen at −80°C (20 min), incubated at 65°C (2 h), and consequently treated at 95°C (10 min) for inaction of the proteinase K enzyme. DNA samples were stored at −20°C until their use as PCR templates. PCR was carried out in a total volume of 30 μl [19.2 μl distilled water, 3 μl 10× PCR buffer, 0.6 μl 10 mM dNTP mixture, 1.2 μl 50 mM MgCl2, 1.2 μl of each primer (10 pmol/μl), 0.6 μl of Taq DNA polymerase (5 unit/μl, CinnaGen, Tehran, Iran) and 3 μl of DNA template]. The thermal cycling program for amplifying three nematode genomic fragments (SSU, LSU D2-D3 and ITS rDNA) was as follows: denaturation at 95°C for 4 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 52°C for 40 s, and extension at 72°C for 80 s. A final extension was performed at 72°C for 10 min. The COI mtDNA was amplified according to aforementioned PCR program, except the annealing temperature was set to 46°C. The PCR mixture rations and the thermal cycling program for amplifying bacterial 16S rDNA was as already described for nematode genomic fragments amplifications.

The information of the primers used for amplification of aforementioned fragments are given in Table 2.

thumbnail
Table 2. Primers used to amplify and sequence nematode and bacterial ribosomal, and nematode mitochondrial DNA.

https://doi.org/10.1371/journal.pone.0217506.t002

The PCR products were sequenced in both directions using the same primers with an ABI 3730XL sequencer (Applied Biosystems) at Macrogen (Seoul, South Korea). Newly obtained sequences of the studied nematode species and their endosymbiont bacteria were deposited into the GenBank database under the accession numbers given in Table 3.

thumbnail
Table 3. The species and isolate names, newly generated accession numbers and the sequenced regions.

https://doi.org/10.1371/journal.pone.0217506.t003

Phylogenetic analyses

The newly obtained SSU, LSU D2-D3, ITS, COI mtDNA and 16S rDNA sequences were compared with those of other nematodes and bacteria sequences available in GenBank database using the BLAST homology search program.

The DNA sequences for reconstructing the phylogenetic trees were selected according to the recent studies [10, 21, 34, 48]. The selected sequences were updated after BLAST search.

The LSU rDNA D2-D3 and COI mt DNA sequences (94 and 75 sequences respectively, including outgroup sequences) were aligned using the MUSCLE algorithm implemented in MEGA 6 [50]; and the resultant alignment was manually edited.

The nematode ITS sequences (53 sequences including outgroup sequences) and the bacterial 16S rDNA (42 sequences for Candidatus Xiphinamtobacter spp. and 25 sequences for Burkholderiaceae spp. phylogenies respectively, including outgroup sequences) were aligned using the Q-INS-i algorithm of online version of MAFFT version 7 (http://mafft.cbrc.jp/alignment/server/) [51]. The Gblocks program (version 0.91b) with the three less stringent parameters, a server tool at the Castresana Lab (http://molevol.cmima.csic.es/castresana/ Gblocks_server.html), was used for post-editing of this alignment to eliminate poorly aligned regions or divergent positions.

The most appropriate model of nucleotide substitution was selected using the Akaike information criterion (AIC) in MrModeltest 2 [52]. The general time-reversible model, including a gamma distribution for rates across sites and a proportion of invariant sites (GTR+G+I) was selected for all phylogenetic analyses in this study. Bayesian inference (BI) was performed using MrBayes v3.1.2 [53], running the chains for 5 × 106 generations (for all datasets). After discarding burn-in samples, the remaining samples were retained for further analyses. The Markov-chain Monte Carlo method within a Bayesian framework was used to estimate the posterior probabilities of the phylogenetic trees [54], using the 50% majority rule. Adequacy of the posterior sample size was evaluated using autocorrelation statistics as implemented in TRACER v.1.5 [55]. A maximum likelihood (ML) tree was reconstructed with RaxmlGUI 1.1 [56] software and the same nucleotide substitution model as that used for BI including 1000 bootstrap (BS) pseudoreplicates. For the LSU rDNA dataset, the species Tylencholaimus teres Thorne, 1939 [57] and Ecumenicus monohystera De Man, 1880 [58], (accession number EF207243, AY593013 respectively), for the ITS rDNA phylogeny, the species Heterodorus veletensis Guerrero, Leibanas & Pena-Sañtiago, 2007 [59] (accession number EU477380) and for COI phylogeny, Xiphinema diversicaudatum (Micoletzky, 1922 [60]) Thorne, 1939 [57] (KF292292) and X. index Thorne & Allen, 1950 [61] (HM921380) (according to Orlando et al. [10]) were used as outgroup taxa. For the bacterial 16S rDNA phylogeny of Candidatus Xiphinematobacter spp., two species Prosthecobacter fusiformis Staley, Bont & Jonge, 1976 [62] (U60015) and Verrucomicrobium spinosum Schlesner, 1987 [63] (X90515), and for the Burkholderiaceae 16S rDNA phylogeny, the species Candidatus Glomeribacter gigasporarum Bianciotto, Lumini, Bonfante & Vandamme, 2003 [64] (AJ251634) were used as outgroup species in accordance with previous studies [10, 34].

The resultant phylogenetic trees were visualized using Dendroscope V.3.2.8 [65] and re-drawn in CorelDRAW software version 16.

Co-phylogenetic analyses

In order to understand how the new species identified in our study, and its endosymbiont, fit into the broader evolutionary picture of the Xiphinema americanum-group, we undertook a cophylogenetic analysis of the nematodes and their endosymbionts. Our analysis examined cophylogenetic signal between the endosymbiont 16S phylogenies and both nematode phylogenies (one based on the LSU sequences, the other based on COI sequences). Genetic data was available for nine pairs of host and endosymbiont for the LSU region and seven pairs of species for the COI region. Sequences for the host nematode species were selected when sequences for their endosymbiont bacteria were also available. Cophylogenetic signal in the associations of host and endosymbiont can allude to a parallel evolutionary history between the two as is expected from results elsewhere in similar systems [6668].

We used the cophylogenetic tool, PACo [69, 70], to perform our analyses. PACo uses a Procrustes approach to superimpose one phylogeny on the other based on the observed host-endosymbiont associations and estimates cophylogenetic signal as the sum of squared residuals (ss) from this superimposition [69]. Observed cophylogenetic signal is considered significant if the observed test statistic is significantly smaller than the test statistic observed when the matrix of host-endosymbiont associations is randomized and the analysis re-run [70]. In our analysis, we performed PACo analysis with the asymmetric statistic which assesses the goodness-of-fit of the endosymbiont phylogeny onto the host phylogeny (since we would expect host evolution to drive endosymbiont evolution). We compared our observed cophylogenetic signal to the same from 10,000 randomizations of the association matrix and considered cophylogenetic signal significant if the observed signal was smaller than the random expectation at α = 0.05. We also examined the results from the perspective of individual host-endosymbiont associations to examine how cophylogenetic signal of our newly identified species compared to that between the other host-endosymbiont pairs.

Nomenclatural acts

The electronic vision of this paper meets the requirements of the amended international code of zoological nomenclature (ICZN), and hence the new name contained herein is available under that code from the electronic vision of this paper. This published work and the contained nomenclatural acts have been registered in the online registration system for the ICZN in ZooBank. The ZooBank LSID (life science identifiers) for this publication is: urn:lsid:zoobank.org:pub:96D34294-16C8-45C2-BEFD-5A0A5C5AA148. The related LSID information can be viewed through any standard web browser by appending the LSID to the prefix "http://www.zoobank.org/References/". The electronic edition of this work was published in a journal with an ISSN, and has been archived and is available from the following digital repositories: PubMed Central, LOCKSS.

Results

Xiphinema primum Mobasseri, Hutchinson, Jahanshahi Afshar & Pedram n. sp. urn:lsid:zoobank.org:act:6CC653DD-05BE-4F9A-A3BB-990E73B0D2BF

(Figs 1, 2 and 3)

thumbnail
Fig 1. Line drawings of Xiphinema primum n. sp. (Amol population).

(A) Pharynx. (B) Part of female reproductive system, showing endosymbiont bacteria inside the ovary. (C) Anterior region. (D and E) Entire female. (F) Pharyngeal bulb. (G-J) Tail of juvenile developmental stages from J1-J4, respectively. (K-N) Tail of female.

https://doi.org/10.1371/journal.pone.0217506.g001

thumbnail
Fig 2. Light microphotographs of Xiphinema primum n. sp. (Amol population).

(A) Entire female. (B-D) Anterior end. (E, H and K) Pharyngeal bulb. (F and J) Ovary and endosymbiont bacteria. (G) Odontophore. (I) Vulval region. (L-P) Tail of female. (All scale bars = 10 μm, except A = 100 μm).

https://doi.org/10.1371/journal.pone.0217506.g002

thumbnail
Fig 3.

Light microphotographs of Xiphinema primum n. sp. (A-C, E, G, H-J: Tehran population, D, F, K: Ramsar population, paratype female). (A-D) Anterior region. (E) Pharyngeal bulb. (F and G) Ovary and endosymbiont bacteria. (H) Part of female reproductive system, showing less developed uterus and oviduct. (I-K) Tail of female. (All scale bars = 10 μm).

https://doi.org/10.1371/journal.pone.0217506.g003

Description.

Female. Body cylindrical, very gradually narrowing towards both extremities, assuming a loose spiral after fixation (Fig 1D and 1E). Cuticle two layered, smooth under light microscopy, fine transverse striae visible on outer layer mostly in tail. The hyaline part of tail 11–15 μm thick, occupying about 1/3 of tail length. Lateral cord 12–28 μm wide, occupying about 1/4 of corresponding body diam. Lip region widely rounded, separated from the rest body by a shallow depression (Fig 1C and Fig 2B–2D and Fig 3A–3D). Amphidial fovea cup-shaped; aperture 12–13 μm wide, located slightly anterior to lip region–body junction region. Odontostyle ca. 1.7 times longer than odontophore, the latter with well-developed flanges. Guiding ring double and guiding sheath 5–10 μm long depending on the degree of protraction/retraction of the stylet. Pharyngeal bulb 78–85 μm long and 19.5–26.0 μm wide. The larger dorsal gland nucleus (DN) located at 15.3–24.7% of pharyngeal bulb length, the two smaller ventrosublateral nuclei (S1N) located at about the same level, at 50.2–60.3% of terminal bulb length. Cardia 9–19 × 5–13 μm in size. Intestine simple, prerectum often not well seen, and rectum 0.6–0.9 times the anal body diameter. Reproductive system didelphic-amphidelphic, the branches about equal in size, the endosymbiont bacteria within the ovaries visible in most examined individuals (Fig 1B and Fig 2F and 2J and Fig 3F–3H), uterus short ca 50 μm long, and the sphincter not clearly seen. Vagina perpendicular to body axis, about 1/2 to 1/3 of corresponding body diameter, the pars distalis vaginae 15–20 μm long, the pars proximalis vaginae about 1.5 times wider than tall. Tail short, conical, dorsally more convex, with bluntly rounded end (Fig 1K–1N and Fig 2M–2P and Fig 3I–3K). Measurements are listed in Table 4.

thumbnail
Table 4. Morphometric data for females of Xiphinema primum n. sp. from Iran.

All measurements in μm and in the format: mean ± s.d. (range).

https://doi.org/10.1371/journal.pone.0217506.t004

Male. Not found.

Juveniles. Recovered in one of the three populations of the new species (the type population with code D, from Ramsar). All four juvenile developmental stages were identified and separated from each other according to Robbins et al. [71]. The correlation between body length of juveniles and females, replacement odontostyle of juveniles and functional odontostyle of females and juveniles is given in Fig 4. The first juvenile (J1) was characterized by the replacement odontostyle tip close to base of functional odontostyle, located at the level of the odontophore, and three other stages had replacement odontostyle posterior to odontophore flanges in resting position of the spear. All stages had dorsally convex conical tails, with a bluntly rounded terminus, however tail in J1and J2 was narrower (Fig 1G–1J). Measurements are listed in Table 5.

thumbnail
Fig 4. Graph of correlation of functional and replacement odontostyle to body length in juveniles and females of X. primum n. sp. of the type population with the code D.

The replacement odontostyle in each juvenile stage is equal to the functional odontostyle in the next stage, and in J4, is equal to the functional odontostyle in females. (Ods. = odontostyle).

https://doi.org/10.1371/journal.pone.0217506.g004

thumbnail
Table 5. Morphometric data of paratype juveniles of Xiphinema primum n. sp. from Iran.

All measurements in μm and in the format: mean ± s.d. (range).

https://doi.org/10.1371/journal.pone.0217506.t005

Type host and locality

Rhizospheric soils of an unidentified forest tree in Ramsar, Mazandaran province (the population code D) (GPS coordinates: 36°55′57.22″N, 50°36′47.00″E).

Other localities

Two other populations of the new species were recovered from the rhizospheric soil of Fagus sp. collected from Amol, Mazandaran province (the population code Dkond), GPS coordinates: 36°24′7.89″N, 52°18′13.054″E; and rhizospheric soil of unidentified trees in village of Āhār, Shemiranat county, Tehran province (the population code Dbana), GPS coordinates: 35°56′13.518″N, 51°27′ 6.615″E.

Type materials

Holotype, six paratype females and 14 juvenile paratypes (J1-J4) were deposited in the Nematode Collection at the Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran. Seven females of Tehran population (the Dbana population) were deposited in the USDA Nematode Collection, Beltsville, MD, USA and 10 females of Amol population (the Dkond population) were deposited in Ghent University Museum, Zoology Collections, Ghent, Belgium. Specific sequences are deposited in GenBank with accession numbers given in Table 3.

Etymology

The specific epithet shows this is the first species of X. americanum-group described from Iran.

Diagnosis and relationships

Xiphinema primum n. sp. belongs to the brevicolle-complex of the Xiphinema americanum-group, mainly characterized by the lip region being offset from the rest of the body by a slight depression and conical tail, with a bluntly rounded tip. It is further characterized by 2124–2981 μm long females with a widely rounded lip region, separated from the rest of the body by a shallow depression, a 103–125 μm long odontostyle, two equally developed genital branches with visible endosymbiont bacteria under light microscopy, vulva located at 51.8–58.0%, 26–37 μm long conical tail with a more convex dorsal side, bluntly rounded tip and four juvenile developmental stages. The alpha-numeric identification codes of the new species according to Lamberti et al. [19] are: A56, B23, C1, D23, E3, F-, G1, H2, I123.

By having a widely rounded lip region, and conical tail with bluntly rounded tip, the new species comes close to nine species of the brevicolle-complex sensu Orlando et al. [10] and Palomares-Rius et al. [34] viz. X. brevicolle Lordello & Da Costa, 1961 [72], X. diffusum Lamberti & Bleve-Zacheo, 1979 [9], X. incognitum Lamberti & Bleve-Zacheo, 1979 [9], X. himalayense Ahmad, Lamberti, Rawat, Agostinelli & Srivastava, 1998 [73], X. luci Lamberti & Bleve-Zacheo, 1979 [9], X. parabrevicolle Gutiérrez-Gutiérrez, Navarrete, Decraemer, Vovlas, Prior, Palomares-Rius & Castillo, 2012 [14], X. paramonovi Romanenko, 1981 [74], X. parataylori Lazarova, Oliveira, Prior, Peneva & Kumari, 2019 [21] and X. taylori Lamberti, Ciancio, Agostinelli & Coiro, 1991 [75]. The comparisons with aforementioned species are as follows:

Differentiated from X. brevicolle, besides distant placement in LSU, ITS and COI mtDNA trees, by a longer body (2497.7 (2124–2981) vs (1800–2200) μm), posteriorly located guiding ring (100 (91–112) vs 72.3–92.3 μm) from anterior end and slightly longer tail (31 (26–37) vs 23.0–31.2) μm) (data of X. brevicolle based on topotypes given by Lamberti et al. [9] and Lazarova et al. [21]

Differentiated from X. diffusum, besides distant placement in LSU, ITS and COI mtDNA trees, by longer body (2497.7 (2124–2981) vs 1630 (1400–1800) μm), odontostyle and odontophore (113.5 (103–125) vs 85.5 (73–95) μm; and 68.3 (60–74) vs 49.7 (44–57) μm, respectively), more posteriorly located guiding ring (100 (91–112) vs 70.2 (60–77) μm) and greater body width at anus (35 (30–42) vs 22.8 (20–30) μm) [9].

Differentiated from X. parabrevicolle, the sister taxon in COI mtDNA tree by distant placement in LSU and ITS trees, longer odontophore (68.3 (60–74) vs 58.3 (54.5–61.0) μm), greater c′ value (0.9 (0.8–1.0) vs 0.7 (0.7–0.8)) and longer tail (31 (26–37) vs 26.6 (24.5–30) μm) [14].

Differentiated from X. himalayense, the tentative cryptic species of the new species, besides distant placement in both LSU and ITS trees, by lip region separated from the body by a constriction (vs continuous), conical tail, dorsally more convex (vs conoid, ventrally straight or slightly concave with rounded terminus) and wider body at anus (35 (30–42) vs 30.5 (28.9–32.4) [73].

Differentiated from X. paramonovi by a longer body (2497.7 (2124–2981) vs 2100 (2000–2300) μm), narrower lip region (12.5 (11–14) vs 14.6 (13.5–15) μm), longer odontophore (68.3 (60–74) vs 56.7 (53–60) μm), greater c value (81 (66.2–95.6) vs 60.5 (49.1–68)), smaller c′ value (0.9 (0.8–1.0) vs 1.1 (0.9–1.2)), greater body width at vulva (52.5 (41–65) vs 43.4 (39.0–47.2) μm) and shorter tail (31 (26–37) vs 36.1 (33–47) μm) [74].

Differentiated from X. parataylori, besides distant placement in LSU, ITS and COI mtDNA trees, by posteriorly located vulva (V = 54.7 (51.8–58.0) vs 49.7 (47–50)), longer odontostyle (113.5 (103–125) vs 91 (85–98) μm) and odontophore (68.3 (60–74) vs 58.6 (54–66) μm), posteriorly located guiding ring (100 (91–112) vs 78.5 (73–89) μm from anterior end) [21].

Differentiated from X. taylori, besides distant placement in LSU and COI mtDNA trees, by lip region separated from the rest body by a depression (vs constriction), longer body (2497.7 (2124–2981) vs 2300 (1800–2300) μm), odontostyle (113.5 (103–125) vs 95 (83.0–100.6) μm) and odontophore (68.3 (60–74) vs 58.8 (54.0–61.8) μm) and greater V value (54.7 (51.8–58.0) vs 50.4 (48–52)%), greater body width at anus (35 (30–42) vs 30.1 (25–34.1) μm) [75].

Differentiated from X. incognitum, besides distant placement in LSU, ITS and COI mtDNA trees, by longer body (2497.7 (2124–2981) vs 1900 (1700–2100) μm), greater a and c values (48 (39.2–61.7) vs 45 (41–49) and (81 (66.2–95.6) vs 62 (47–75), respectively), posteriorly located vulva (V = 54.7 (51.8–58.0) vs 51 (48–53%), longer odontostyle and odontophore (113.5 (103–125) vs 87 (82–93) μm, and (68.3 (60–74) vs 52 (46–56) μm, respectively), wider anal body region (35 (30–42) vs 28 (24–33) μm), and narrower tail tip (vs wider) [9].

Differentiated from X. luci, besides distant placement in LSU and ITS trees, by longer body (2497.7 (2124–2981) vs 1800 (1700–1900) μm), longer odontostyle and odontophore (113.5 (103–125) vs 95 (93–99) μm, and (68.3 (60–74) vs 50 (47–53) μm, respectively), posteriorly located guiding ring (100 (91–112) vs 76 (68–85) μm from anterior end), greater body width at anus (35 (30–42) vs 24 (20–30) μm), smaller c′ value (0.9 (0.8–1.0) vs 1.2 (1.0–1.4)) and posteriorly located vulva (V = 54.7 (51.8–58.0) vs 51 (49–52)%) [9].

Iranian population of X. pachtaicum from Amol

(Fig 5)

thumbnail
Fig 5. Light microphotographs of Iranian population of Xiphinema pachtaicum (Tulaganov, 1938) Kirjanova, 1951.

(A-D) Anterior region. (E) Pharyngeal bulb. (F) Vulval region. (G and H) Endosymbiont bacteria inside the ovary. (I-L) Female tail. (All scale bars = 10 μm).

https://doi.org/10.1371/journal.pone.0217506.g005

Two populations of the species were recovered from the rhizospheric soil samples of Citrus sp. and Alnus sp. trees, collected respectively in a garden and a forest close to the city of Amol in Mazandaran province, northern Iran (GPS coordinates 36°26′39″N, 52°19′8.8″E and 36°24′7.8″N, 52°18′13.05″E, respectively). The recovered populations were morphologically compared with some other populations of the species from different parts of the world. Compared with the Spanish populations reported by Gutiérrez-Gutiérrez et al. [13], Hungarian populations reported by Repasi et al. [76], and Yugoslavian populations reported by Barsi & Lamberti [20], no remarkable differences in the presently studied population were observed. Compared to the American populations described by Orlando et al. [10], differences in vulva location (V = 55.9–60.5% in Iranian populations vs 48–57%) were observed. As mentioned by Orlando et al. [10], American populations show differences with other populations of the species from other parts of the world for this index. Compared to the Iranian populations of the species reported by Fadaei et al. [77], only minor intraspecies variation was observed for body width at anus (17–24 vs 13.5–19 μm). Measurements are listed in Table 6.

thumbnail
Table 6. Morphometric data of Xiphinema pachtaicum (Tulaganov, 1938) Kirjanova, 1951 from Iran.

All measurements in μm and in the format: mean ± s.d. (range).

https://doi.org/10.1371/journal.pone.0217506.t006

Molecular phylogenetic relationships.

Molecular characterization of Xiphinema primum n. sp. and Iranian population of X. pachtaicum

The basic local alignment search tool (BLAST) results using newly obtained sequences given in Table 3

The BLAST search using partial LSU rDNA D2-D3 sequences of four isolates of the new species (MF372941, MF372942, MF372943, MF372947) revealed they have 99% identity with three accession numbers MH248814, MH248815 and AY601602 recorded in the GenBank as Xiphinema sp. Kra-BG, Xiphinema sp. Per-BG and X. taylori, respectively.

The BLAST search using ITS sequences of four isolates of the new species (MF372951, MF372952, MF372953, MF372954) revealed their identity with available sequences having a high coverage were never higher than 91%.

The BLAST search using COI mtDNA sequences of four isolates of the new species (MK202795, MK202796) revealed their identity with currently available sequences is at maximum 90%.

The BLAST search using partial sequences of LSU rDNA D2-D3 expansion segments of the Iranian isolates of X. pachtaicum (accession numbers MF372944, MF372945, MF372946) revealed they have 99–100% identity with sequences of other isolates of the species (KU250155, HM921389, KJ802888, HM921397, HM921365, HM921353, JQ990033, KP268968 and some other sequences).

The BLAST search using the 16S rDNA of the endosymbiont bacterium from three isolates of X. pachtaicum, revealed they have 99–100% identity with several sequences of Burkholderiacea bacterium (KT735078, KT735074, KT735068, KU899555, KU899554, KU899551, KT735076, KT735082 and some others) already isolated from the same species.

The BLAST search using the 16S rDNA of the endosymbiont bacterium from two isolates of X. primum n. sp. revealed they have at maximum, 95% identity with several sequences of Candidatus Xiphinematobacter spp. (KJ614453, KJ614452, KJ614450, AF217462, KJ614451, KT735101, KT735096).

The BLAST search using partial SSU rDNA of an isolate of X. pachtaicum sequenced during the present study (MF372940), revealed it has 99% identity with other sequences of the species (KU250139, KP407873, AM086682, MH484520) and an isolate of X. penevi Lazarova, Peneva & Kumari, 2016 [16] (KU250141) available in GenBank. No phylogenetic analyses were performed in the present study using this sequence.

The D2-D3 expansion segments of LSU rDNA phylogeny based on Xiphinema primum n. sp. and Iranian population of X. pachtaicum sequences.

This dataset was composed of 727 characters of which 314 characters were variable. Fig 6 represents the Bayesian phylogenetic tree reconstructed using this dataset. The tree is divided to two major clades I and II, according to the nomenclature used in previous studies [10, 18]. In this tree, the new species appeared as an independent lineage in a clade containing brevicolle-complex spp. However, the relationships between most species in this clade were not resolved due to polytomy. The Iranian populations of X. pachtaicum formed a clade with two other selected isolates of the species (KU250155, HM921390). The five isolates of X. pachtaicum clade had almost identical LSU rDNA D2-D3 sequences (only one indel in overlapping region was observed), however, the length of their sequences differed, yielding intraspecies cladogenesis and varied branch lengths inside this species clade. The topotype population of X. brevicolle is marked by asterisk in this tree, corroborating its restricted distribution [21].

thumbnail
Fig 6. Bayesian 50% majority rule consensus tree inferred using LSU rDNA D2-D3 sequences of X. primum n. sp. and Iranian populations of X. pachtaicum under the GTR+G+I model.

Bayesian posterior probabilities and maximum likelihood bootstrap values are given for appropriate clades in the shape BPP/ML BS. The newly generated sequences are in bold. According to Lazarova et al. [21], the accession number assigned to X. brevicolle marked by asterisk belongs to the topotype population of the species.

https://doi.org/10.1371/journal.pone.0217506.g006

The ITS phylogeny based on Xiphinema primum n. sp. and Iranian population of X. pachtaicum sequences.

The ITS dataset was composed of 1077 characters of which 452 characters were variable. Fig 7 represents the Bayesian phylogenetic tree reconstructed using this dataset. In this tree, a similar pattern to that observed in LSU phylogeny was observed, i.e., the sequences in this tree were divided into two major clades I and II. However, the relationships between the clades inside group II were not resolved due to polytomy. The new species occupies a basal placement to the clade of the X. brevicolle-complex, mostly having similar morphology, and appeared as an independent lineage. Iranian isolates of X. pachtaicum are not included in this tree as attempts to get their ITS sequences failed. Again, the topotype population of X. brevicolle is marked by asterisk, corroborating its restricted distribution [21].

thumbnail
Fig 7. Bayesian 50% majority rule consensus tree inferred using ITS sequences of X. primum n. sp. under the GTR+G model.

Bayesian posterior probabilities and maximum likelihood bootstrap values are given for appropriate clades in the shape BPP/ML BS. The newly generated sequences are in bold. According to Lazarova et al. [21], the accession number assigned to X. brevicolle marked by asterisk belongs to the topotype population of the species.

https://doi.org/10.1371/journal.pone.0217506.g007

The COI mtDNA phylogeny based on Xiphinema primum n. sp. and Iranian population of X. pachtaicum sequences.

The COI dataset was composed of 348 characters of which 201 characters were variable. Fig 8 represents the Bayesian phylogenetic tree reconstructed using this dataset. In this tree, the relationship between the new species and other species and several members of the X. brevicolle-complex is not resolved due to polytomy, and it is sister to X. parabrevicolle and X. peruvianum Lamberti & Bleve-Zacheo, 1979 [9] with poor clade support. Similar to two previous phylogenies, the topotype population of X. brevicolle is marked by asterisk, corroborating its restricted distribution [21].

thumbnail
Fig 8. Bayesian 50% majority rule consensus tree inferred using COI mtDNA sequences of X. primum n. sp. under the GTR+G model.

Bayesian posterior probabilities and maximum likelihood bootstrap values are given for appropriate clades in the shape BPP/ML BS. The newly generated sequences are in bold. According to Lazarova et al. [21], the accession number assigned to X. brevicolle marked by asterisk belongs to the topotype population of the species.

https://doi.org/10.1371/journal.pone.0217506.g008

Molecular characterization of endosymbiont bacterial species

Ca. Xiphinematobacter spp. phylogeny.

The 16S rDNA dataset of Ca. Xiphinematobacter spp., including two newly amplified sequences of the bacterial 16S rDNA from the new species, included 526 total characters of which 158 characters were variable. Fig 9 represents the Bayesian phylogenetic tree reconstructed using this dataset. In this tree, the relation of the major clades was not resolved due to polytomy, and the newly generated sequences (MF372955, MF372956) were placed inside the clade including sequences of endosymbionts from X. brevicolle-complex (AF217462, KJ614448, KX263241, KJ614453). However, the clade received poor support (0.66/-) and different alignment, alignment post-editing, and reconstructing methods need to be applied to test any putative scenario for the relationship between these sequences.

thumbnail
Fig 9. Bayesian 50% majority rule consensus tree inferred using 16S rDNA sequences of the endosymbiont bacteria, isolated from Xiphinema primum n. sp. under the GTR+G+I model.

Bayesian posterior probabilities and maximum likelihood bootstrap values are given for appropriate clades in the shape BPP/ML BS. The newly generated sequences are in bold. The host nematodes are indicated in right.

https://doi.org/10.1371/journal.pone.0217506.g009

The Burkholderiaceae spp. phylogeny.

The 16S rDNA dataset of Burkholderiaceae spp. including three newly generated sequences of endosymbionts from X. pachtaicum had 1228 characters of which 214 characters were variable. In this tree (Fig 10), the newly generated sequences occupied a placement inside the clade of several sequences of Burkholderiaceae from X. pachtaicum.

thumbnail
Fig 10. Bayesian 50% majority rule consensus tree inferred using 16S rDNA sequences of the endosymbiont bacteria, isolated from Iranian populations of Xiphinema pachtaicum under the GTR+G+I model.

Bayesian posterior probabilities and maximum likelihood bootstrap values are given for appropriate clades in the shape BPP/ML BS. The newly generated sequences are in bold. The host nematodes are indicated in right.

https://doi.org/10.1371/journal.pone.0217506.g010

Cophylogenetic analyses

We examined the cophylogenetic patterns between selected sequences of Xiphinema americanum-group spp. and Ca. Xiphinematobacter spp. using both COI and LSU sequences for nematodes. In both cases, we observed strong, significant cophylogenetic signal between the two groups (LSU: ss = 0.076, n = 10,000, p < 0.01, Fig 11A; COI: ss = 0.168, n = 10,000, p < 0.01, Fig 12A). In Figs 11A and 12A, we plot the phylogenies of endsymbionts and hosts (LSU and COI, respectively) and blue lines connect host with endosymbiont where line thickness represents cophylogenetic signal strength (thinner lines represent smaller discrepancy between host and endosymbiont and therefore stronger cophylogenetic signal, dashed lines represent the species identified in this study). Only 0.5% and 0.1% of randomized host-endosymbiont associations, respectively, generated a stronger cophylogenetic signal than the observed associations. When comparing host-endosymbiont associations to each other, we saw that the endosymbiont-host pairs identified in this study display a medium level of cophylogenetic association (Figs 11B and 12B) with some associations (e.g. pairs 4, 6, 7 in Fig 11A and pair 6 in Fig 12A) showing stronger signal and others showing a weaker signal (e.g. pairs 3 & 9 in Fig 11A and pairs 1–5 in Fig 12A). In Figs 11B and 12B, we plot the distribution of cophylogenetic signal between pairs of hosts and endosymbionts and mark the species/populations described in this study with vertical lines that demonstrate their medium level of cophylogenetic signal relative to the rest of the group for the LSU marker and relatively strong signal with the COI marker. Overall, these results confirm the established pattern of strong cophylogenetic signal between these two groups, identify average levels of cophylogenetic signal between the host and endosymbiont introduced in this study, and provide another example of parallel cladogenesis between endosymbionts and their hosts.

thumbnail
Fig 11. Cophylogenetic analysis of Xiphinema americanum spp. and their endosymbionts of Ca. Xiphinematobacter spp. using nematode LSU rDNA D2-D3 and bacteria 16S rDNA.

(A) The phylogenies of endosymbiont bacteria alongside that of nematode host (Xiphinema americanum spp.). Blue lines connect host and endosymbiont and the thickness of each line represents the strength of cophylogenetic signal between host and endosymbiont. Thinner lines reflect stronger cophylogenetic signal (i.e. line thickness is the level of discordance between host and endosymbiont). Dashed blue lines and red numbers distinguish the species described in this study from the others used in the analysis. (B) The cophylogenetic signal of individual host-endosymbiont associations. Cophylogenetic signal is displayed on the x-axis where smaller values represent stronger cophylogenetic signal. The density profile (grey) describes the cophylogenetic signal between the all pairs of species and the two pairs identified in this study are also represented by the blue vertical lines.

https://doi.org/10.1371/journal.pone.0217506.g011

thumbnail
Fig 12. Cophylogenetic analysis of Xiphinema americanum spp. and their endosymbionts of Ca. Xiphinematobacter spp. using nematode COI mtDNA and bacteria 16S rDNA.

(A) The phylogenies of endosymbiont bacteria alongside that of nematode host (Xiphinema americanum spp.) Blue lines connect host and endosymbiont and the thickness of each line represents the strength of cophylogenetic signal between host and endosymbiont. Thinner lines reflect stronger cophylogenetic signal (i.e. line thickness is the level of discordance between host and endosymbiont). Dashed blue lines and red numbers distinguish the species described in this study from the others used in the analysis. (B) The cophylogenetic signal of individual host-endosymbiont associations. Cophylogenetic signal is displayed on the x-axis where smaller values represent stronger cophylogenetic signal. The density profile (grey) describes the cophylogenetic signal between the all pairs of species and the pair identified in this study are also represented by the blue vertical line.

https://doi.org/10.1371/journal.pone.0217506.g012

Discussion

The nematodes of the Xiphinema americanum-group are infamous for their difficulty in species delimitation using traditional criteria. The conserved morphology of even phylogenetically distant species or the close morphometric data ranges further complicate their identification. Using the mean values of the morphometric data as proposed by Orlando et al. [10] and applied in this study (see Differential and Diagnosis), could help primary species comparisons.

In this study, a new species in the Xiphinema americanum-group was described from Iran, as the first Iranian representative of the group, using an integrative approach. It belongs to X. brevicolle-complex, mainly characterized by a bluntly conoid tail, and the lip region not separated from the body by a sharp constriction. The species X. himalayense was regarded as its tentative cryptic species that besides morphological traits, and especially the tail characters, could be separated from it based on molecular data.

The molecular phylogenetic studies of the group were performed in some previous studies [10, 21, 34]. Using LSU rDNA, He et al. [12] found that the species of this group are divided into two subclades: the X. americanum and the X. brevicolle subgroups [12]. That latter was recently considered as “X. brevicolle species complex” [10]. The former phylogenetic studies based upon LSU rDNA D2-D3 and ITS markers, usually yielded congruent phylogenies [14], dividing the species into two main subclades I and II [10, 14, present study]. Two aforementioned studies, the present study, and the most recent study by Lazarova et al. [21] show that the COI mtDNA marker has acceptable interspecific variation to be used as an alternative approach that is helpful in overcoming the low divergence of LSU rDNA D2-D3 [10]. X. primum n. sp. was however separated from its morphologically close species and its cryptic species by both LSU and ITS markers in the corresponding phylogenetic trees.

The cryptic nature of some species and the misidentification of some species in GenBank remains an issue. The observed intraspecific divergence of D2-D3 marker in X. brevicolle s. l. [10, 14] documents the cryptic nature of nematodes of this group well. Sequencing of the topotype population of X. brevicolle s. str. by Lazarova et al. [21] revealed several populations are indeed misidentified and are erroneously assigned to this species. As a result, any observed paraphyletic or polyphyletic status of a given species, should be checked by the correct identification of that species first, and more preferably, topotype individuals of the studied species should be included in sequencing experiments whenever possible.

The recovered Iranian populations of X. pachtaicum and the new species were inspected for the presence of endosymbiont bacteria inside their ovaries. The rod-shaped endosymbiont bacteria were observed inside the ovaries of both recovered species under light microscopy (Fig 5G and 5H). The amplification of 16S rDNA and sequencing of the amplified fragments, revealed that the endosymbiont bacterium of the new species belongs to Ca. Xiphinematobacter sp., forming an independent lineage in the corresponding phylogenetic tree, and was not identical to other available sequences of this genus. The endosymbiont relation of Ca. Xiphinematobacter spp. with X. americanum-group members has recently been studied by several authors [34, 48, 10]. The endosymbiont of X. pachtaicum had an affinity with plant- and fungi-associated species of Burkholderiaceae, and was placed in the clade of species already isolated from X. pachtaicum. According to Palomares-Rius et al. [34], X. pachtaicum could harbor both endosymbionts; either a Ca. Xiphinematobacter sp. or a Burkholderiaceae sp.

Significant cophylogenetic signal has already been observed between nematodes of this group and both groups of the endosymbiont bacteria phylogenies, albeit with variation in signal strength [10, 34]. In this study, significant cophylogenetic signal was observed between both recovered nematode species and their endosymbionts phylogeny. The pair of species introduced in this study showed similar cophylogenetic signal to the rest of the group and do not disrupt the overall signal within the group. These results confirm the generally held expectation that host-endosymbiont evolution often occurs in parallel and indeed the evolution in the dependent group, parallels that of the host [6668] as previously documented for X. americanum-group members and their endosymbionts [10, 34]. Indeed, the vertical transmission of endosymbionts via parthenogenesis means a significant cophylogenetic signal is expected.

Conclusion

Hereby, we studied two Iranian populations of Xiphinema pachtaicum and described a new species, both of which belonged to the X. americanum-group, using an integrated approach. Significant cophylogeny between the newly generated sequences of nematode LSU rDNA D2-D3 and COI mtDNA sequences and the 16S rDNA sequences of the endosymbiont bacteria were also observed.

Acknowledgments

The technical helps of Dr. Yousef Panahandeh are appreciated.

References

  1. 1. Cobb NA. New nematode genera found inhabiting fresh water and non-brackish soils. Journal of the Washington Academy of Sciences. 1913; 3(16):432–44.
  2. 2. Taylor CE, Brown DJ. Nematode vectors of plant viruses. Wallingford, UK. CAB International; 1997. 286. p.
  3. 3. Loof PAA, Luc M. A revised polytomous key for the identification of species of the genus Xiphinema Cobb, 1913 (Nematoda: Longidoridae) with exclusion of the X. americanum-group. Systematic Parasitology. 1990; 16:35–66.
  4. 4. Luc M, Coomans A, Loof PA, Baujard P. The Xiphinema americanum-group (Nematoda: Longidoridae). 2. Observations on Xiphinema brevicollum Lordello & da Costa, 1961 and comments on the group. Fundamental and Applied Nematology. 1998; 21(5):475–90.
  5. 5. Lamberti F, Molinari S, Moens M, Brown DJ. The Xiphinema americanum group. I. Putative species, their geographical occurrence and distribution, and regional polytomous identification keys for the group. Russian Journal of Nematology. 2000; 8(1):65–84.
  6. 6. Coomans A, Huys R, Heyns J, Luc M. Character analysis, phylogeny and biogeography of the genus Xiphinema Cobb, 1913 (Nematoda: Longidoridae). Annales Musée Royal de L'Afrique Centrale (Zoologie), Tervuren. Belgique. 2001; 287:1–289.
  7. 7. Lima MB. Studies on species of the genus Xiphinema and other nematodes. Ph.D. Thesis, University of London, UK. 1965; 165 pp.
  8. 8. Tarjan AC. Variation within the Xiphinema americanum group (Nematoda: Longidoridae). Nematologica. 1969; 15(2):241–52.
  9. 9. Lamberti F, Bleve-Zacheo T. Studies on Xiphinema americanum sensu lato with descriptions of fifteen new species (Nematoda, Longidoridae). Nematologia Mediterranea. 1979; 7(1):51–106.
  10. 10. Orlando V, Chitambar JJ, Dong K, Chizhov VN, Mollov D, Bert W, Subbotin SA. Molecular and morphological characterisation of Xiphinema americanum-group species (Nematoda: Dorylaimida) from California, USA, and other regions, and co-evolution of bacteria from the genus Candidatus Xiphinematobacter with nematodes. Nematology. 2016; 18(9):1015–43. https://doi.org/10.1163/15685411-00003012
  11. 11. Lamberti F, Molinari S, Moens M, Brown DJ. The Xiphinema americanum-group. 11. Morphometric relationships. Russian Journal of Nematology. 2002; 10:99–112.
  12. 12. He Y, Subbotin SA, Rubtsova TV, Lamberti F, Brown DJ, Moens M. A molecular phylogenetic approach to Longidoridae (Nematoda: Dorylaimida). Nematology. 2005; 7(1):111–24. https://doi.org/10.1163/1568541054192108
  13. 13. Gutiérrez-Gutiérrez C, Rius JE, Cantalapiedra-Navarrete C, Landa BB, Castillo P. Prevalence, polyphasic identification, and molecular phylogeny of dagger and needle nematodes infesting vineyards in southern Spain. European Journal of Plant Pathology. 2011; 129(3):427–53. https://doi.org/10.1007/s10658-010-9705-y
  14. 14. Gutiérrez-Gutiérrez C, Cantalapiedra-Navarrete C, Decraemer W, Vovlas N, Prior T, Rius JE, Castillo P. Phylogeny, diversity, and species delimitation in some species of the Xiphinema americanum-group complex (Nematoda: Longidoridae), as inferred from nuclear and mitochondrial DNA sequences and morphology. European Journal of Plant Pathology. 2012; 134(3):561–97. https://doi.org/10.1007/s10658-012-0039-9
  15. 15. Lazarova SS, Malloch G, Oliveira CM, Hübschen J, Neilson R. Ribosomal and mitochondrial DNA analyses of Xiphinema americanum-group populations. Journal of Nematology. 2006; 38(4):404–410. pmid:19259456
  16. 16. Lazarova S, Peneva V, Kumari S. Morphological and molecular characterisation, and phylogenetic position of X. browni sp. n., X. penevi sp. n. and two known species of Xiphinema americanum-group (Nematoda, Longidoridae). ZooKeys. 2016; (574):1–42. pmid:27110175
  17. 17. Kumari S, Decraemer W, De Luca F, Tiefenbrunner W. Cytochrome c oxidase subunit 1 analysis of Xiphinema diversicaudatum, X. pachtaicum, X. simile and X. vuittenezi (Nematoda, Dorylaimida). European Journal of Plant Pathology. 2010; 127(4):493–9. https://doi.org/10.1007/s10658-010-9614-0
  18. 18. Archidona-Yuste A, Navas-Cortés JA, Cantalapiedra-Navarrete C, Palomares-Rius JE, Castillo P. Cryptic diversity and species delimitation in the Xiphinema americanum-group complex (Nematoda: Longidoridae) as inferred from morphometrics and molecular markers. Zoological Journal of the Linnean Society. 2016; 176(2):231–65. https://doi.org/10.1111/zoj.12316
  19. 19. Lamberti F, Hockland S, Agostinelli A, Moens M, Brown DJ. The Xiphinema americanum group. III. Keys to species identification. Nematologia Mediterranea. 2004; 32(1): 53–6.
  20. 20. Barsi L, Lamberti F. Xiphinema parasimile sp. n. from Serbia and X. simile, first record from Bosnia and Herzegovina (Nematoda, Dorylaimida). Nematologia Mediterranea. 2004; 32(1):101–9.
  21. 21. Lazarova S, Oliveira CM, Prior T, Peneva V, Kumari S. An integrative approach to the study of Xiphinema brevicolle Lordello and Da Costa 1961, supports its limited distribution worldwide (Nematoda: Longidoridae). European Journal of Plant Pathology. 2019; 153(2), 441–464. https://doi.org/10.1007/s10658-018-1571-z
  22. 22. Tulaganov AT. The fauna of nematodes of cotton and surrounding soil in Katta-Kurgan district of the Uzbek SSR. Trudy Uzbekskogo Gosudarstvennogo Universiteta. 1938; 12:1–25.
  23. 23. Kirjanova ES. Nematodypoehvykholkovogopolja i theliny v Golodnoj step (Uzbekistan). Trudy Zoologicheskogo Instituta Akademiya Nauk SSSR. 1951; 9:625–57.
  24. 24. EPPO. EPPO Standards. EPPO A1 and A2 lists of pests recommended for regulation as quarantine pests. Available online at http://archives.eppo.int/EPPOStandards/PM1_GENERAL/pm1-02(24)_A1A2_2015;pdf.
  25. 25. Siddiqi MR. Xiphinema mediterraneum Martelli et Lamberti a junior synonym of X. pachtaicum (Tulaganov) Kijanova. Nematologia Mediterranea. 1977; 5(1):133–135.
  26. 26. Roca F, Lamberti F, Agostinelli A. I Longidoridae (Nematoda, Dorylaimida) delle regioni italiane. IX. La sicilia. Nematologia Mediterranea. 1989; 17(2):151–165.
  27. 27. Decraemer W, Robbins RT. The who, what and where of Longidoridae and Trichodoridae. Journal of Nematology. 2007; 39(4):295–297. pmid:19259501
  28. 28. Pedram M. Nematodes of the family Longidoridae. In: Ghaderi, R., Kashi, L. & Karegar, A. (Eds). Plant-parasitic nematodes in Iran. Science Reference in collaboration with the Iranian Society of Nematology. 2018; 627–67. (In Persian).
  29. 29. Adams RE, Eichenmuller JJ. A bacterial infection of Xiphinema americanum. In Phytopathology 1963; 53(7):745
  30. 30. Coomans A, Claeys M. Structure of the female reproductive system of Xiphinema americanum (Nematoda: Longidoridae). Fundamental and Applied Nematology. 1998; 21(5):569–80.
  31. 31. Coomans A, Vandekerckhove TT, Claeys M. Transovarial transmission of symbionts in Xiphinema brevicollum (Nematoda: Longidoridae). Nematology. 2000; 2(4):443–9. https://doi.org/10.1163/156854100509303
  32. 32. Vandekerckhove TT, Willems A, Gillis M, Coomans A. Occurrence of novel verrucomicrobial species, endosymbiotic and associated with parthenogenesis in Xiphinema americanum-group species (Nematoda, Longidoridae). International Journal of Systematic and Evolutionary Microbiology. 2000; 50(6):2197–205. https://doi.org/10.1099/00207713-50-6-2197
  33. 33. Vandekerckhove TT, Coomans A, Cornelis K, Baert P, Gillis M. Use of the Verrucomicrobia-specific probe EUB338-III and fluorescent in situ hybridization for detection of “Candidatus Xiphinematobacter” cells in nematode hosts. Applied and Environmental Microbiology. 2002; 68(6):3121–5. pmid:12039775
  34. 34. Palomares‐Rius JE, Archidona‐Yuste A, Cantalapiedra‐Navarrete C, Prieto P, Castillo P. Molecular diversity of bacterial endosymbionts associated with dagger nematodes of the genus Xiphinema (Nematoda: Longidoridae) reveals a high degree of phylogenetic congruence with their host. Molecular Ecology. 2016; 25(24):6225–47. pmid:27813204
  35. 35. Garrity GM, Bell JA, Lilburn T. Class I. Alphaproteobacteria class. nov. InBergey’s Manual of Systematic Bacteriology. Springer, Boston, MA. 2005. 574. p.
  36. 36. Howe DK, Smith M, Tom DM, Brown AM, Peetz AB, Zasada IA, Denver DR. Analysis of nematode-endosymbiont coevolution in the Xiphinema americanum species complex using molecular markers of variable evolutionary rates. Nematology. 2019; 21(5):533–46. https://doi.org/10.1163/15685411-00003233
  37. 37. Thorne G. Notes on free-living and plant-parasitic nematodes, II. Proceedings of the Helminthological Society of Washington. 1935; 2(2):96–8.
  38. 38. Pedram M, Pourjam E, Namjou S, Atighi MR, Cantalapiedra-Navarrete C, Liébanas G, Palomares-Rius JE, Castillo P. Molecular and morphological characterisation of Paralongidorus iranicus n. sp. and P. bikanerensis (Lal & Mathur, 1987) Siddiqi, Baujard & Mounport, 1993 (Nematoda: Longidoridae) from Iran. Nematology. 2012; 14(4):427–43. https://doi.org/10.1163/156854111X610272
  39. 39. Gharibzadeh F, Pourjam E, Pedram M. Description of Longidorus azarbaijanensis n. sp. (Dorylaimida: Longidoridae) from Iran. Journal of Nematology. 2018; 50:207–18. pmid:30451438
  40. 40. Jenkins W. A rapid centrifugal-flotation technique for separating nematodes from soil. Plant Disease Reporter. 1964; 48(9): 692 p.
  41. 41. De Grisse AT. Redescription ou modifications de quelques technique utilis [a] es dan l'etude des n [a] ematodes phytoparasitaires. Mededelingen van de Rijksfaculteit der Landbouwwetenschappen Gent. 1969; 34:351–69.
  42. 42. Williams BD, Schrank B, Huynh C, Shownkeen R, Waterston RH. A genetic mapping system in Caenorhabditis elegans based on polymorphic sequence-tagged sites. Genetics. 1992; 131(3):609–24. PMCID: PMC1205034 pmid:1321065
  43. 43. Holterman M, van der Wurff A, van den Elsen S, van Megen H, Bongers T, Holovachov O, Bakker J, Helder J. Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown clades. Molecular Biology and Evolution. 2006; 23(9):1792–800. pmid:16790472
  44. 44. Nunn GB. Nematode molecular evolution. Ph.D. dissertation. Nottingham: University of Nottingham; 1992.
  45. 45. Kiontke K, Gavin NP, Raynes Y, Roehrig C, Piano F, Fitch DH. Caenorhabditis phylogeny predicts convergence of hermaphroditism and extensive intron loss. Proceedings of the National Academy of Sciences. 2004; 101(24):9003–8. https://doi.org/10.1073/pnas.0403094101
  46. 46. Subbotin SA, Halford PD, Warry A, Perry RN. Variations in ribosomal DNA sequences and phylogeny of Globodera parasitising solanaceous plants. Nematology. 2000; 2(6):591–604. https://doi.org/10.1163/156854100509484
  47. 47. He Y, Jones J, Armstrong M, Lamberti F, Moens M. The mitochondrial genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): considerable economization in the length and structural features of encoded genes. Journal of Molecular Evolution. 2005; 61(6):819–33. pmid:16315110
  48. 48. Lazarova SS, Brown DJ, Oliveira CM, Fenton B, MacKenzie K, Wright F, Malloch G, Neilson R. Diversity of endosymbiont bacteria associated with a non-filarial nematode group. Nematology. 2016; 18(5):615–23. https://doi.org/10.1163/15685411-00002982
  49. 49. Brosius J, Palmer ML, Kennedy PJ, Noller HF. Complete nucleotide sequence of a 16S ribosomal RNA gene from Escherichia coli. Proceedings of the National Academy of Sciences. 1978; 75(10):4801–5. https://doi.org/10.1073/pnas.75.10.4801
  50. 50. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution. 2013; 30(12):2725–9. pmid:24132122
  51. 51. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution. 2013; 30(4):772–80. pmid:23329690
  52. 52. Nylander JA. MrModeltest v2. Program distributed by the author. Evolutionary Biology Center, Uppsala University. 2004.
  53. 53. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003; 19(12):1572–4. pmid:12912839
  54. 54. Larget B, Simon DL. Markov chain Monte Carlo algorithms for the Bayesian analysis of phylogenetic trees. Molecular Biology and Evolution. 1999; 16(6):750–9.
  55. 55. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007; 7(1):214. https://doi.org/10.1186/1471-2148-7-214
  56. 56. Silvestro D, Michalak I. raxmlGUI: a graphical front-end for RAxML. Organisms Diversity and Evolution. 2012; 12(4):335–7. https://doi.org/10.1007/s13127-011-0056-0
  57. 57. Thorne G. A monograph of the nematodes of the superfamily Dorylaimoidea. Capita Zoologica. 1939; 8(5): 261 pp.
  58. 58. De Man JG. Die einheimischen, frei in der reinen Erde und im süssen Wasser lebende Nematoden, monographischen bearbeitet. Tijdschrift Nederlandse Dierkundige vereniging. 1880. 104. p.
  59. 59. Guerrero P, Liébanas G, Peña-Santiago R. Nematodes of the order Dorylaimida from Andalucía Oriental, Spain. The genus Enchodelus Thorne, 1939. 1. Description of two new and two known species with conical tail. Nematology. 2007; 9(4):515–36. https://doi.org/10.1163/156854107781487251
  60. 60. Micoletzky H. Die freilebenden Erd-Nematoden: mit besonderer Berücksichtigung der Steiermark und der Bukowina, zugleich mit einer Revision sämtlicher nicht mariner, freilebender Nematoden in Form von Genus-Beschreibungen und Bestimmungsschlüsseln. Archiv für naturgeschichte, Abteilung A. 1922; 7:59–79.
  61. 61. Thorne G, Allen MW. Paratylenchus hamatus n. sp, and Xiphinema index n. sp., two nematodes associated with fig roots, with a note on Paratylenchus anceps Cobb. Proceedings of the Helminthological Society of Washington. 1950; 17(1):27–35.
  62. 62. Staley JT, De Bont JA, de Jonge K. Prosthecobacter fusiformis nov. gen. et sp., the fusiform caulobacter. Antonie Van Leeuwenhoek. 1976; 42(3):333–42. https://doi.org/10.1007/BF00394132 pmid:1086646
  63. 63. Schlesner H. Verrucomicrobium spinosum gen. nov., sp. nov.: a fimbriated prosthecate bacterium. Systematic and Applied Microbiology. 1987; 10(1):54–6. https://doi.org/10.1016/S0723-2020(87)80010-3
  64. 64. Bianciotto V, Lumini E, Bonfante P, Vandamme P. ‘Candidatus Glomeribacter gigasporarum’gen. nov., sp. nov., an endosymbiont of arbuscular mycorrhizal fungi. International Journal of Systematic and Evolutionary Microbiology. 2003; 53(1):121–4. https://doi.org/10.1099/ijs.0.02382-0
  65. 65. Huson DH, Scornavacca C. Dendroscope 3: an interactive tool for rooted phylogenetic trees and networks. Systematic Biology. 2012; 61(6):1061–7. pmid:22780991
  66. 66. Hafner MS, Page RD. Molecular phylogenies and host-parasite cospeciation: gophers and lice as a model system. Philosophical Transactions of The Royal Society of London Series B-Biologic.1995; 349(1327):77–83. https://doi.org/10.1098/rstb.1995.0093
  67. 67. Maneesakorn P, An R, Daneshvar H, Taylor K, Bai X, Adams BJ, Grewal PS, Chandrapatya A. Phylogenetic and cophylogenetic relationships of entomopathogenic nematodes (Heterorhabditis: Rhabditida) and their symbiotic bacteria (Photorhabdus: Enterobacteriaceae). Molecular Phylogenetics and Evolution. 2011; 59(2):271–80. pmid:21335093
  68. 68. Rosenblueth M, Sayavedra L, Sámano‐Sánchez H, Roth A, Martínez‐Romero E. Evolutionary relationships of flavobacterial and enterobacterial endosymbionts with their scale insect hosts (Hemiptera: Coccoidea). Journal of Evolutionary Biology. 2012; 25(11):2357–68. pmid:22994649
  69. 69. Balbuena JA, Míguez-Lozano R, Blasco-Costa I. PACo: a novel procrustes application to cophylogenetic analysis. PloS One. 2013; 8(4):e61048. pmid:23580325
  70. 70. Hutchinson MC, Cagua EF, Balbuena JA, Stouffer DB, Poisot T. paco: implementing Procrustean Approach to Cophylogeny in R. Methods in Ecology and Evolution. 2017; 8(8):932–40.
  71. 71. Robbins RT, Brown DJ, Halbrendt JM, Vrain TC. Compendium of juvenile stages of Xiphinema species (Nematoda: Longidoridae). Russian Journal of Nematology. 1996; 4(2):163–71. ISSN: 0869–6918
  72. 72. Lordello LG, Costa CD. A new nematode parasite of coffee roots in Brazil. Revista Brasileira de Biologia. 1961; 21(4):363–6. ISSN: 0034–7108
  73. 73. Ahmad M, Lamberti F, Rawat VS, Agostinelli A, Srivastava N. Two new species within the Xiphinema americanum-group (Nematoda, Dorylaimida) from Garhwal Himalayas, India. Nematologia Mediterranea. 1998; 26(2):131–8.
  74. 74. Romanenko ND. Discovery of a new species of nematode Xiphinema paramonovi n. sp. (Nematoda: Longidoridae) in the territory of the Soviet Union. Tezisy Dokladov, Pervaya Konferentsiya po Nematodam Rastenii, Nasekomykh, Pochvy i Vod. 1981; 68–9.
  75. 75. Lamberti FD, Ciancio A, Agostinelli A, Coiro MI. Relationship between Xiphinema brevicolle and X. diffusum whit a redescription of C. brevicolle and descriptions of three new species of Xiphinema (Nematoda: Dorylaimida). Nematologia Mediterranea. 1991; 19(2):311–26.
  76. 76. Repasi V, Agostinelli A, Nagy P, Coiro MI, Hecker K, Lamberti F. Distribution and morphometrical characterization of Xiphinema pachtaicum, X. simile and X. brevicollum from Hungary. Helminthologia. 2008; 45(2):96–102. https://doi.org/10.2478/s11687-008-0018-z
  77. 77. Fadaei AA, Coomans A, Kheiri A. Three species of the Xiphinema americanum lineage (Nematoda: Longidoridae) from Iran. Nematology. 2003; 5(3):453–61. https://doi.org/10.1163/156854103769224430