The genus Aphelenchoides Fischer, 1894 (family Aphelenchoididae Skarbilovich, 1947), has around 175 nominal species (Mobasseri et al., 2018). This could also be an underestimate, and several sequences deposited at the GenBank database for example, are identified only at the genus level. This species-rich genus is the type genus of the family Aphelenchoididae, and is well known by prevalence of species lacking conspicuous apomorphies, helpful for its species delimitation. The term “foliar nematodes” is the common name of plant-parasitic forms commonly used by plant pathologists. Actually, some Aphelenchoides species are important plant parasites (Kanzaki and Giblin-Davis, 2012). They could be recovered from soil, mosses, mushrooms, decaying organic materials or in some cases, from plant tissues (Khusainov, 2013). Historically, several authors, e.g., Hunt (1993, 2008), Shahina (1996) and Andrássy (2007) have provided list of valid species for the genus. The book of Baranovskaya (1981), however, is a useful resource, especially for some species with inaccessible descriptions. The Index to Organism Names website (www.organismnames.com), however, includes updated data of most of the newly described species. According to Shahina (1996), the species could be grouped based on the number of their lateral lines, i.e., the species usually have obscure, unknown, 2, 3 and 4 lateral lines, while species with 4 to 6, and 6 lateral lines do also rarely occur. Andrássy (2007), however, pointed out that six lines in lateral fields rarely occur. As far as our knowledge, until 1995 that the key of Shahina (1996) was published, A. shamimi, Khera 1970 was the only, having five lines in the lateral fields. Since 1995 until 2008-2009 that the checklist of Hunt (2008) was published, it was only A. paramonovi Eroshenko and Kruglik, 2004 that had five lines in the lateral fields. No species having five lines in lateral fields has been described since 2008 to 2009 till date. During recent years, revisions have been performed on taxonomic status and placement of some Aphelenchoides species. In a recent study, the species Laimaphelenchus heidelbergi Zhao et al., 2007 was transferred to the genus Aphelenchoides (Carta et al., 2016). In another study, the species Tylaphelenchus christinae Lieutier and Laumond, 1978 was transferred to the genus Aphelenchoides (Pedram et al., 2018a, 2018b) and finally, two species A. subtenuis (Cobb, 1926) Steiner and Buhrer, 1932 and A. arachidis Bos, 1977 have been transferred to the genus Robustodorus Andrássy, 2007 (Kanzaki et al., 2018). The latter action, in revising the taxonomic placement of those two species, was supported by both traditional and molecular criteria.
The history of the reported or described species of Aphelenchoides from Iran was given by Mobasseri et al. 2018. During our nematological surveys conducted in northern provinces of the country, two populations of the genus were recovered from Asalem forests in Gilan province, and a forest in Semnan province in north of Iran. Further morphological, morphometric and phylogenetic studies, especially detailed SEM studies, revealed the first species belongs to an unknown species, described herein as Aphelenchoides giblindavisi n. sp. The second species was identified as A. helicus Heyns, 1964, its detailed morphological studies using fresh females yielded new observations, meriting its transferring to the genus Robustodorus, an action that was supported by the molecular phylogenetic analyses using two genomic markers too.
Present study aims to (i) describe Aphelenchoides giblindavisi n. sp. and characterize it using morphological and molecular data and (ii) to revise the current taxonomic status of A. helicus using newly observed morphological traits and molecular phylogenetic criteria.
Materials and methods
Sampling, nematode extraction, and morphological observation
Several soil, moss, rotten wood, bark and insect cadavers were collected from the natural forests in Gilan and Semnan provinces, northwestern and northern Iran. The samples were placed in plastic bags, transferred to nematology laboratory of Tarbiat Modares University and maintained at 4 °C. Nematodes were extracted using the tray method (Whitehead and Hemming, 1965), heat killed by adding hot 4% formalin solution, transferred to anhydrous glycerin according to De Grisse (1969), and mounted on permanent slides. Both species were studied in detail in temporary slides using fresh females in water. Photographs were taken using an Olympus DP72 digital camera attached to an Olympus BX51 microscope powered with differential interference contrast (DIC). For SEM studies, after their examination and identification, a few specimens preserved in glycerin were selected for observation under SEM following the protocol of Álvarez-Ortega and Peña-Santiago (2016). The nematodes were hydrated in distilled water, dehydrated in a graded ethanol and acetone series, critical point-dried, coated with gold, and observed with a Zeiss Merlin Scanning Electron Microscope.
DNA Extraction, PCR and sequencing
DNA of two recovered species was extracted from one single female nematode. Each specimen was picked out, studied onto a temporary slide, transferred to a small drop of TE buffer (10 mM Tris-Cl, 0.5 mM EDTA; pH 9.0, 100 QIAGEN Inc., Valencia CA, USA) on a clean slide and squashed using a clean slide cover glass. The suspension was collected by adding 15 μl of the aforementioned buffer. The DNA samples were stored at −20 °C until using 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 two genomic fragments (SSU and LSU rDNA D2–D3) was as follows: denaturation at 95 °C for 4 min, followed by 32 cycles of denaturation at 94 °C for 30 sec, annealing at 52 °C for 40 sec, and extension at 72 °C for 80 sec. A final extension was performed at 72 °C for 10 min (Pedram, 2017; Mobasseri et al., 2017). Primers for 28S rDNA D2–D3 amplification were forward primer D2A (5′-ACAAGTACCGTGAGGGAAAGT-3′) and reverse primer D3B (5′-TGCGAAGGAACCAGCTACTA-3′) (Nunn, 1992). Primers for amplification of 18S rDNA were forward primer SSU 1813F (5′-CTGCGTGAGAGGTGAAAT-3′) and reverse primer SSU 2646R (5′-GCTACCTTGTTACGACTTTT-3′) as used by Holterman et al. (2006).
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 were deposited into the GenBank database (accession numbers: MG545999 for the partial SSU and MG546000 for the partial LSU rDNA D2–D3 of the new species, and KP264116 for the SSU and KP264117 for the partial LSU rDNA D2–D3 of A. helicus).
Phylogenetic analyses
The newly obtained SSU and LSU rDNA D2–D3 sequences were compared with those of other nematode species available in GenBank using the BLAST homology search program. The selected sequences for reconstructing of each phylogenetic trees were aligned using the Q-INS-i algorithm of online version of MAFFT version 7 (http://mafft.cbrc.jp/alignment/server/) (Katoh and Standley, 2013). The Gblocks program (version 0.91b) with all 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 the alignments, i.e., to eliminate poorly aligned regions or divergent positions. The most appropriate model of nucleotide substitution was selected using the Akaike information criterion in MrModeltest 2 (Nylander, 2004). 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 both datasets. Bayesian inference (BI) was performed using MrBayes v3.1.2 (Ronquist and Huelsenbeck, 2003) running the chains for five million generations. 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 Bayesian posterior probabilities (BPP) of the phylogenetic trees (Larget and Simon, 1999) using the 50% majority rule. Adequacy of the posterior sample size was evaluated using autocorrelation statistics as implemented in TRACER v.1.5 (Drummond and Rambaut, 2007). A maximum likelihood (ML) tree was reconstructed by using RaxmlGUI 1.1 (Silvestro and Michalak, 2012) software using the same nucleotide substitution model as in the BI including 1000 bootstrap (BS) pseudoreplicates. For SSU phylogeny, the species, Panagrolaimus detritophagus Fuchs, 1930, Plectonchus sp. and Brevibucca saprophaga Goodey, 1935 (accession numbers EU543176, AY593920 and EU196018 respectively) and for LSU rDNA D2–D3 dataset, the species Panagrellus redivivus (Linnaeus, 1767) Goodey, 1945 and Poikilolaimus piniperdae Fuchs, 1930 (accession numbers DQ408249 and DQ059060 respectively) were used as outgroup taxa (according to previous studies, e.g., van Megen et al., 2009; Ryss et al., 2013; Pedram, 2017). The inferred trees were visualized using Dendroscope V.3.2.8 (Huson and Scornavacca, 2012) and re-drawn in CorelDRAW software version 16.
Results
Systematics
Aphelenchoides giblindavisi n. sp.
Table 1
Morphometrics of Aphelenchoides giblindavisi n. sp. from Iran. All measurements are in μm and in the form: mean ± S.D. (range).
| Female | Male | ||
|---|---|---|---|
| Character | Holotype | Paratypes | Paratypes |
| n | – | 12 | 12 |
| L | 653 | 671 ± 72.4 (546–795) | 632 ± 41.4 (523–679) |
| a | 39.6 | 34.9 ± 2.2 (31.4–39.6) | 37.4 ± 2.2 (32.4–41.0) |
| b | 8.6 | 8.9 ± 0.8 (7.4–10.1) | 8.7 ± 0.4 (8.0–9.4) |
| b′ | 4.5 | 4.5 ± 0.4 (3.9–5.0) | 4.5 ± 0.2 (4.2–5.0) |
| c | 18.7 | 18.9 ± 1.0 (17.2–20.4) | 18.9 ± 1.7 (15.8–21.5) |
| c′ | 3.5 | 3.4 ± 0.2 (3.1–3.7) | 2.8 ± 0.2 (2.6–3.1) |
| V or T | 70.1 | 70.6 ± 1.6 (68.8–73.5) | 52.8 ± 3.4 (47.4–59.2) |
| Head height | 2.5 | 2.5 ± 0.0 (2.0–2.5) | 2.5 ± 0.0 (2.0–2.5) |
| Head diam. | 5.5 | 5.5 ± 0.4 (5–6) | 5.5 ± 0.3 (5–6) |
| Stylet | 10.5 | 10.5 ± 0.5 (10–11) | 10.0 ± 0.4 (10–11) |
| Stylet conus | 4.5 | 4.5 ± 0.2 (4.0–4.5) | 4.0 ± 0.2 (3.5–4.5) |
| m | 42.9 | 38.9 ± 2.3 (36.4–42.9) | 38.4 ± 2.1 (35–40) |
| Median bulb | 58 | 58.0 ± 3.2 (54–65) | 56.0 ± 1.4 (54–58) |
| Excretory pore | 69 | 69.0 ± 6.2 (60–84) | 66.5 ± 1.8 (64–69) |
| Hemizonid | 80 | 89.0 ± 10.6 (77–110) | 81.0 ± 2.6 (78–85) |
| Pharynx | 76 | 75.5 ± 2.8 (71–81) | 72.0 ± 3.4 (65–78) |
| Nerve ring | 75 | 75.0 ± 2.8 (71–81) | 73.0 ± 2.3 (70–78) |
| Median bulb length | 13.5 | 14 ± 1 (12–15) | 13.0 ± 0.9 (11–14) |
| Median bulb diam. | 10.5 | 11.0 ± 1.1 (9.0–12.5) | 10.0 ± 0.9 (8.0–11.5) |
| Median bulb length/diam. | 1.3 | 1.3 ± 0.1 (1.1–1.4) | 1.3 ± 0.1 (1.2–1.6) |
| Pharyngeal overlapping | 70 | 73.0 ± 6.7 (63–87) | 67.0 ± 6.4 (60–80) |
| Maximum body diam. | 16.5 | 19.0 ± 2.4 (15–23) | 17.0 ± 1.4 (14–19) |
| Vulval body diam. (VBD) | 16 | 18.0 ± 1.8 (14.5–20.0) | – |
| Body diam. at median bulb | 13.5 | 14.0 ± 1.4 (12.0–15.5) | 13.0 ± 0.8 (11–14) |
| Postvulval uterine sac (PUS) | 65 | 62.4 ± 6.0 (52–69) | – |
| PUS/VBD | 4.1 | 3.5 ± 0.4 (2.7–4.1) | – |
| Vulva–anus | 160 | 161 ± 21.6 (114–190) | – |
| Ovary or testis length | 295 | 257 ± 58.7 (183–356) | 334 ± 35 (248–394) |
| Anal (cloacal) body width | 10 | 10 ± 1 (9.0–12.5) | 12.0 ± 0.6 (10.5–13.0) |
| Tail length | 35 | 35.0 ± 2.5 (30–39) | 34.0 ± 1.8 (31–37) |
| Spicules length (arc line) | – | – | 17.0 ± 0.7 (16–18) |
| Capitulum | – | – | 7.0 ± 0.5 (6–8) |
| Females | |
|---|---|
| n | 15 |
| L | 464 ± 37.3 (395–547) |
| a | 29.2 ± 1.8 (26.7–33.5) |
| b | 7.8 ± 0.6 (6.7–8.8) |
| b′ | 4.6 ± 0.4 (3.8–5.4) |
| c | 19.7 ± 1.7 (16.8–22.8) |
| c′ | 2.8 ± 0.3 (2.4–3.1) |
| V | 69.3 ± 1.6 (64.4–71.1) |
| Head height | 2.0 ± 0.3 (2.0–2.5) |
| Head diam. | 5.0 ± 0.4 (4.5–5.5) |
| Stylet | 10.0 ± 0.4 (10–11) |
| Stylet conus | 4.0 ± 0.3 (4.0–4.5) |
| m | 41.3 ± 2.8 (36.4–45.0) |
| Median bulb | 47.5 ± 2.9 (39–51) |
| Excretory pore | 53.0 ± 4.8 (46–62) |
| Hemizonid | 63.0 ± 5.6 (53–69) |
| Pharynx | 60 ± 4 (49–65) |
| Nerve ring | 58 ± 4 (45–62) |
| Median bulb length | 11.5 ± 0.9 (10.0–13.5) |
| Median bulb diam. | 10.0 ± 0.8 (8–11) |
| Median bulb length/diam. | 1.2 ± 0.1 (1.1–1.3) |
| Pharyngeal overlapping | 40 ± 5 (31–49) |
| Maximum body diam. | 16.0 ± 1.6 (13–20) |
| Vulval body diam. (VBD) | 15.0 ± 1.4 (12–19) |
| Body diam. at median bulb | 13.0 ±1.1 (11.5–15.0) |
| Postvulval uterine sac (PUS) | 38 ± 8 (28–55) |
| PUS/VBD | 2.5 ± 0.5 (2.0–3.7) |
| Vulva–anus | 119 ± 13 (101–147) |
| Ovary length | 161 ± 20 (130–193) |
| Anal body diam. | 8.0 ± 0.9 (7–10) |
| Tail length | 24.0 ± 1.4 (21–26) |








