Occurrence of Sheraphelenchus sucus (Nematoda: Aphelenchoidinae) and Panagrellus sp. (Rhabditida: Panagrolaimidae) Associated with Decaying Pomegranate Fruit in Italy

Fruit symptoms and nematode extraction: Pomegranate fruit showing rotting symptoms (Fig. 1) were collected at Metaponto (Matera province, Basilicata) and at Fasano (Brindisi province, Apulia), in Southern Italy on October 2015. Three out of twelve sampled localities revealed the occurrence of decaying fruit, still hanging on the tree, with simultaneous presence of fungi, bacterial- and fungal-feeder nematodes. Pomegranate fruit was fragmented into small pieces of 1 × 2 cm which were incubated in water in Petri dishes. After 2–3 hr of incubation, the emerged nematodes were used to conduct morphological and molecular observations.

FIG. 1.

Adult specimens of the mycetophagus population for microscopic observation were handpicked and killed by gentle heat, fixed in a solution of 4% of formaldehyde +1% propionic acid, and processed to glycerol according to Seinhorst’s method (Hooper, 1986). Glycerin-infiltrated specimens were examined by light microscopy for diagnosis. Photographs were taken with a Leica DFC 425 system mounted on a Leitz Wetzlar optical microscope, using specimens mounted on water agar temporary slides (Troccoli, 2002) (Fig. 2), whereas measurements were made at the camera lucida on glycerine infiltrated specimens. Abbreviations and terms used in the morphometric data presentation as defined in Siddiqi (2000).

FIG. 2.

DNA extraction, PCR amplification, and sequencing: Individual nematodes were crushed with a sterile microspatula under a stereomicroscope. Genomic DNA was extracted from fifteen and ten individual mycetophagus and bacterial nematodes, respectively, as described by De Luca et al. (2004). The crude DNA isolated from each individual nematode was directly amplified. The ITS1-5.8S-ITS2 regions were amplified using the forward primer TW81 (5′-GTTTCCGTAGGTGAACCTGC-3′) and the reverse primer AB28 (5′-ATATGCTTAAGTTCAGCGGGT-3′) (Joyce et al., 1994); the 18S rDNA was amplified using the 18SnF (5′-TGGATAACTGTGGTAATTCTAGAGC-3′) and 18SnR (5′-TTACGACTTTTGCCCGGTTC-3′) (Kanzaki and Futai, 2002); the D2A-D3B expansion segments of 28S rRNA gene amplified using the primers D2A (5′-ACAAGTACCGTGGGGAAAGTTG-3′) and the D3B (5′-TCGGAAGGAACCAGCTACTA-3′) (Nunn, 1992) and the portion of the mitochondrial cytochrome oxidase c subunit 1 (mtCOI) gene was amplified with this primer set: COI-F1 (5′-CCTACTATGATTGGTGGTTTTGGTAATTG-3′) and COI-R2 (5′-GTAGCAGCAGTAAAATAAGCACG-3′) (Kanzaki and Futai, 2002). PCR cycling conditions used for amplification were an initial denaturation at 94°C for 5 min, followed by 35 cycles of denaturation at 94°C for 50 sec, annealing at 55°C for 50 sec, extension at 72°C for 1 min, and a final step at 72°C for 7 min. The size of the amplification products was determined by comparison with the molecular weight marker ladder 100 (Fermentas, St. Leon-Rot, Germany) after electrophoresis of 10 μl on a 1% agarose gel. PCR products of the ITS containing region, the 18 S gene, the D2–D3 expansion domains, and the mitochondrial COI, from three individual nematodes, were purified using the protocol suggested by the manufacturer (High Pure PCR elution kit; Roche, Germany). Purified DNA fragments were cloned and sent for sequencing, in both directions, at MWG-Eurofins in Germany.

Phylogenetic analysis: A BLAST (Basic Local Alignment Search Tool) search at NCBI (National Center for Biotechnology Information) was performed to confirm their nematode origins and species (Altschul et al., 1997). The newly obtained sequences for the ITS, 18S rRNA gene, the D2–D3 expansion domains of 28S gene, and mtcoI were aligned using ClustalW with default parameters with the corresponding published gene sequences of Aphelenchoididae. Sequence alignments were manually edited using BioEdit to improve the multialignment. Outgroup taxa for each dataset were chosen based on the results of previously published data. Phylogenetic trees, obtained for the ITS, 18S rRNA gene, and the D2–D3 expansion domains were performed with Neighbour-Joining, Minimum Evolution, Maximum Likelihood (ML), and Maximum Parsimony methods using MEGA version 7 software (Kumar et al., 2016). No significant conflict in branching order and support level among methods is observed and, therefore, only ML tree is shown for each marker. The phylograms were bootstrapped 1,000 times to assess the degree of support for the phylogenetic branching indicated by the optimal tree for each method. The newly ITS and D2–D3 sequences of the bacterial nematodes Panagrellus sp. were aligned with the corresponding sequences of Rhabditida. Rhabditophanes sp. and Halicephalobus gingivalis were used as outgroup taxa. The newly obtained sequences were submitted to GenBank with the following accession numbers: LT903555–LT903556 for the ITS of S. sucus; LT837694 for the 18S gene of S. sucus; LT837695–LT837696 for the D2–D3 domains of 28S gene of S. sucus; LT-837697–LT837698 for the mitochondrial COI of S. sucus; LT837699–LT8837700 for the ITS of Panagrellus sp., and LT908055–LT908056 for the D2–D3 domains of 28S gene of Panagrellus sp.

Identification of Sheraphelenchus population: In Table 1, the comparative morphometrics of populations of S. entomophagus (type and Iranian populations) and of S. sucus from the present study and from literature are presented. The morphology of the two populations of S. sucus from Italy closely resembles that of the type population of S. sucus from Japan (Kanzaki and Tanaka, 2013) and that from South Korea (Fang et al., 2015). However, they are also morphologically close to that of the type population of S. entomophagus from Hawaii (Nickle, 1970) as well as to the Iranian S. entomophagus (Moosavi et al., 2014) although with little variation concerning the slightly smaller size of specimens of the latter population. The stylet length, which was originally considered the only diagnostic character useful to separate these species (sensu Kanzaki and Tanaka, 2013), was afterward amended as being underestimated (Kanzaki, 2015) resulting in the close similarity between the two species S. sucus and S. entomophagus.

Molecular identification: Amplification of the ITS, 18S rRNA gene, the D2–D3 expansion segments of the 28S gene, and the mitochondrial COI of Italian S. sucus yielded single fragments of 1,218, 1,610, 779, and 711 bp, respectively.

The ITS sequences of Italian S. sucus from pomegranate showed 99% similarity (8–9 nucleotides different) to S. sucus original population (KC875230 and KC875231). The partial 18S rRNA sequences of S. sucus from pomegranate were identical to S. sucus (AB808720 and KC875229) from the database (99% similarity and 5–8 different nucleotides, respectively) and S. parabrevigulonis (KC875226 and 98% similarity, 30 different nucleotides). The D2–D3 expansion domains of the Italian S. sucus showed a 99% similarity with the corresponding sequences of S. sucus (KC894755 and AB808721; 3 and 8 different nucleotides, respectively) and S. entomophagus from Iran (KF373753 and KF373754; 1 different nucleotide), present in the database. No corresponding sequences were present in the database for the partial mitochondrial COI of the Italian S. sucus population.

Phylogenetic relationships of the Italian population of S. sucus based on the ITS of a multiple edit alignment including 22 sequences, using ML method with Aphelenchus avenae as outgroup are reported in Fig. 3. All populations of S. sucus clustered together in the same clade with 99% of bootstrap support and revealed sister relationships with S. parabrevigulonis.

FIG. 3.

Phylogenetic relationships of the Italian population of S. sucus based on 18S rRNA gene of a multiple edited alignment including 31 sequences, using ML method with Aphelenchus avenae as an outgroup is reported in Fig. 4. This ML tree showed that all populations of S. sucus formed a well-supported group with 100% support and resulted sister species with S. parabrevigulonis.

FIG. 4.

Phylogenetic relationships of S. sucus from Italy based on D2–D3 expansion domains of a multiple edited alignment including 34 sequences were also determined (Fig. 5). Sheraphelenchus sucus populations formed a well-supported grouping with S. entomophagus from Iran, congruent with those obtained by using the 18S rRNA gene and ITS.

FIG. 5.

The ITS containing region of Panagrellus sp. produced a single fragment of 1,261 bp. No corresponding sequences were present in the database. The D2–D3 sequences of Panagrellus sp. from pomegranate showed 99% similarity (8 nucleotide different) with Panagrellus sp. isolated from aberrant adult specimen of R. ferrugineus (Camerota et al., 2016). Thus, selected Panagrellus and rhabditid sequences were aligned with those of D2–D3 of Panagrellus sp. from pomegranate.

In the present study, the molecular identification of the two nematode species occurring in the same rotting pomegranate fruit revealed that the ITS, the partial 18S rRNA gene, the D2–D3 expansion domains, and the partial mitochondrial COI clearly identified S.sucus, whereas the ITS and D2–D3 sequences of the bacterial feeder nematode from pomegranate were identified as Panagrellus sp. The results obtained by sequence and phylogenetic analyses conducted in the present study strongly suggest that S. sucus and S. entomophagus from Iran are conspecific. In addition, morphology and measurements of these two species, including the type population of S. entomophagus from Hawaii fit well with the intraspecific variability of species as reported for other nematode species (Gutiérrez-Gutiérrez et al., 2013; Groza et al., 2017) thus giving strong support for synonymization. However, no molecular data are available for the S. entomophagus type population, suggesting the necessity of reisolation of the S. entomophagus type population from the original locality to organize the species status using morphological characters and molecular barcoding.

Very little is known about the morphological and molecular characterization of Panagrellus species and their interaction with other organisms and corresponding habitats (Stock and Nadler, 2006). Panagrellus species have been already reported as insect associates (Massey, 1974) in bark beetles. In Italy, Smith et al. (1992) described a yeast–nematode–insect association with Drosophila flies in grape sour rot suggesting that Drosophila contributed to the dissemination of these insect associates.

In the present study, phylogenetic relationships of Panagrellus species based on the D2–D3 sequences supported the grouping of Panagrellus sp. from rotting pomegranate with Panagrellus sp. from R. ferrugineus, with 99% overlap (Fig. 6). This finding suggests that the Italian isolates of Panagrellus sp. share the same habitat preference, decaying fruits or rotting materials, and the phoretic relationships occurred casually when they occupy the same habitat of insects (Félix and Duveau, 2012; Camerota et al., 2016; Kanzaki, 2017). These observations strongly suggest that the Italian Panagrellus sp. isolates undergo host switching.

FIG. 6.

The coexistence of the mycetophagus nematode, S. sucus, and the bacterial feeder nematode, Panagrellus sp. in the same rotting pomegranate has never been reported so far. As both nematode species were found in hanging fruit, it is conceivable that they have been carried by insects with which they are phoretically or occasionally associated. More studies are needed to understand the role of phoretic insect–nematode associations in insect ecology and evolution and how they influence reproduction, survival, etc.