Production losses in plantain and banana by Helicotylenchus spp. depend on two main factors. On the one side, the susceptibility of planting material and nematode density, because differences in the population level of H. multicinctus have been registered among Musa spp. cultivars with a high correlation between its density and the necrotic and dead roots percentage (Speijer and Ssango, 1999). On the other hand, losses depend on environmental conditions, since Helicotylenchus spp. are predominant in crops cultivated between 1,000 and 1,350 meters above sea level in andisols and vertisols soils with high contents of Ca, P, Mn and reduced organic matter (Speijer and Ssango, 1999; Karakas, 2007; Araya et al., 2011).
Although Helicotylenchus has been reported in the banana and plantain crops of Colombia and Brazil, morphological, morphometric, and molecular data are scarce for these geographical regions (Zuñiga et al., 1979; Villegas, 1989; Guzmán-Piedrahita and Cataño-Zapata, 2004; Torrado-Jaime and Castaño-Zapata, 2009). In order to contribute to the knowledge of the taxonomic identity of the species of Helicotylenchus associated with Musaceae in Colombia and Brazil, the present study has the following objectives: (i) to identify the species of Helicotylenchus that are associated with banana and plantain crops, using a combination of morphological, morphometric, and molecular analysis, and (ii) to elucidate the phylogenetic relationships of Helicotylenchus species that are associated with banana and plantain crops in the departments of Quindío, Risaralda y Valle del Cauca in Colombia and Minas Gerais in Brazil.
Materials and methods
Area of study, sampling, and nematodes extraction
Root and soil samples were collected from banana and plantain rhizosphere crops in the departments of Quindío (Calarcá and Córdoba municipalities), Risaralda (La Celia municipality), and Valle del Cauca (Palmira and Buenaventura municipalities) in Colombia and in Janauba, Minas Gerais state, Brazil during 2015 to 2018. In each crop, compound samples were collected, comprising 15 to 20 plants/ha. Nematodes were extracted from roots and soil using a modification of Cobb’s method (Ravichandra, 2014). Afterwards, representative populations of Helicotylenchus were selected for posterior analysis in the microbiology and molecular biology laboratories of the National University of Colombia (Palmira, Valle del Cauca, Colombia) and the Laboratory of Nematology of the Biological Institute (Campinas, São Paulo, Brazil).
Morphologic and morphometric analysis
Nematodes extracted and identified as Helicotylenchus were killed with heat at 60°C for 4 min and fixated in 2% formalin. Semipermanent preparations were made and morphometric data were registered following Boag and Shamim Jairajpuri (1985) and Uzma et al. (2015). Microphotographs were taken using a compound microscope equipped with differential interference contrast–DIC (DM2500, Leica, Germany).
Statistical analysis:
Morphometric data were analyzed using the Community Analysis Package (PISCES Conservation Ltd, Lymington UK, 1995) with principal components analysis (PCA) and conglomerate analysis to determine groupings and evaluate those characters that could discriminate species.
Extraction, amplification, and DNA sequencing
For molecular analysis, DNA was extracted by the Proteinase K method of Múnera et al. (2009) with modifications. A single nematode was crushed with a sterile scalpel and transferred to an Eppendorf tube with 15 µl of worm lysis buffer (50 mM KCl, 10 mM Tris pH 8.0, 15 mM MgCl2, 0.5% Triton x–100, 4.5% Tween–20, 0.09% Proteinase K). The tubes were incubated at −80°C (15 min), 65°C (1 h), and 95°C (15 min), centrifuged to 16,000 g (1 min) and stored at −20°C. Amplification of D2 to D3 expansion segment of the large subunit – LSU of ribosomal DNA (28S) was done using forward primer D2A (5′–ACAAGTACCGTGAGGGAAAGTTG–3′) and reverse D3B (5′–TCCTCGGAAGGAACCAGCTACTA–3′) (De Ley et al., 1999). The PCR conditions were initial denaturation during 2 min at 94°C followed by 40 cycles of 45 s at 94°C, 45 s at 55°C, and 1 min at 72°C and final extension of 10 min at 72°C. PCR products were sequenced in both directions at BIONNER Korea.
Phylogenetic analysis
Consensus sequences were edited using the software Geneious (Kearse et al., 2012) and BLAST (Basic Local Alignment Search Tool) at NCBI (National Center for Biotechnology Information) was used to confirm the species’ identity of the sequences. To estimate the phylogenetic history of Helicotylenchus, D2 to D3 sequences of others specimens of the species were downloaded from GenBank, including a sequence of Rotylenchus magnus, which was used as outgroup. Sequence alignment was performed using MAFFT v7 (Katoh et al., 2002) (protocol Q-INS-i), and jModelTest v2.1.7 (Posada, 2008) to find the best nucleotide substitution model, based on the Akaike information criterion corrected for small sample sizes. Afterwards maximum likelihood (ML) was used to estimate a tree with 250 bootstraps and the general time reversible model with allowance for gamma distribution of rate variation (GTR+Γ) in RAxML v8 (Stamatakis, 2014). The phylogeny of Helicotylenchus was inferred using MrBayes v3.2.6 (Ronquist et al., 2012) with the GTR + Γ model. Two independent Metropolis-coupled Markov chain Monte Carlo (MCMCMC) searches were performed for 2 million generations sampling every 2,000 steps. Convergence was assessed using Tracer v1.5 (burn-in = 20% of the samples), and by examining the average standard deviation of split frequencies among parallel chains. The consensus tree was calculated from the posterior distribution of 1,600 phylogenies. The Bayesian analysis was performed in the CIPRES Science Gateway (Miller et al., 2010).
Results
Morphological and morphometric identification of nematodes
Four species of Helicotylenchus were identified in the study areas: H. multicinctus, H. dihystera, H. erythrinae, and H. californicus. Morphological and morphometric data from each species closely resembled type and reference populations species (Tables 1-4).
Morphometric data of H. multicinctus and reference populations.
Note: Measurements in µm; mean ± s.d. (range), except for L in mm.
Helicotylenchus multicinctus populations were identified in plantain crops of Colombia (Córdoba, Quindío and Rozo, Valle del Cauca) and banana crops of Brazil (Minas Gerais). Morphologically, these populations show a habitus post-mortem open C form, a hemispherical cephalic region, the shape of the stylet knobs flattened anteriorly and rounded posteriorly, a rounded tail, functional spermatheca and males present (Fig. 1, Table 1).
Populations of H. dihystera were founded in plantain and banana crops from Colombia (Córdoba, Quindío). Morphologically, the populations of this species show a habitus post-mortem with a spiral shape, a hemispherical cephalic region with the shape of stylet knobs indented or flattened anteriorly and rounded posteriorly, a conoide tail and males absent (Fig. 2, Table 2).
Morphometric data of H. dihystera and reference populations.
Note: Measurements in µm; mean ± s.d. (range), except for L in mm.
The populations identified as H. erythrinae occurred in Colombian plantain and banana crops (Calarcá, Quindío and La Celia, Risaralda). Morphologically, the populations of this species show a habitus post-mortem of a loose spiral, a hemispherical cephalic region with the shape of stylet knobs indented or flattened anteriorly and rounded posteriorly, a tail with a long ventral projection, females with functional spermatheca and males present (Fig. 2, Table 3).
Morphometric data of H. erythrinae and reference populations.
Note: Measurements in µm; mean ± s.d. (range), except for L in mm.
The species H. californicus was identified in a plantain crop from Colombia (Delfina–Buenaventura, Colombia) and a banana crop from Brazil (Minas Gerais). Morphologically, these populations presented a habitus post-mortem in a spiral, a hemispherical cephalic region with the shape of stylet knobs flattened anteriorly and rounded posteriorly, an irregular tail with a short sharply pointed or blunt ventral projection and males present (Fig. 2; Table 4).
PCA and conglomerate analysis separated the species into four groups (H. multicinctus, H. dihystera, H. erythrinae, and H. californicus) (Figs. 3, 4). The principal components 1 to 4 have eigenvalues greater than or equal to 1 and explain 91.39% of variance (Table 5). However, PC1 and PC2 axes better separated the species. In PC1, the most discriminating variables were maximum body diameter, the number of tail annuli, and ratio c′, while in PC2 they were vulva position, anal body diameter, and tail length (Table 6).
Eigenvalues and percent total variance accounted for each principal component.
Principal component
Eigenvalues
Cumulative total
% of Total variance
Cum. % of total variance
1
4.92
4.92
37.88
37.88
2
3.75
8.68
28.88
66.76
3
2.16
10.84
16.63
83.39
4
1.04
11.88
7.99
91.39
5
0.35
12.23
2.70
94.09
6
0.29
12.52
2.22
96.31
7
0.22
12.74
1.68
97.99
8
0.16
12.89
1.23
99.22
9
0.06
12.95
0.43
99.65
10
0.04
12.99
0.29
99.94
11
0.01
13
0.06
100
12
1.5128E–7
13
1.16369E–6
100
13
–4.65844E–10
13
–3.58342E–9
100
Correlations between the first four component principals and the morphometric parameters of females in Helicotylenchus spp.
Vector
Diagnostic character
1
2
3
4
1
L
–0.25
0.21
0.02
0.66
2
a
0.09
–0.01
–0.57
0.35
3
c
–0.34
–0.33
–0.01
0.02
4
c′
0.36
0.28
–0.12
–0.11
5
V %
–0.06
–0.48
–0.13
0.01
6
DGO
–0.22
0.35
–0.25
–0.27
7
Stylet length
–0.31
–0.05
–0.39
–0.03
8
Number of head annuli
0.01
0.05
0.61
0.07
9
Number of tail annuli
0.37
0.07
–0.09
0.33
10
Position of phasmid
–0.29
0.31
–0.09
–0.40
11
Maximum body diameter
–0.39
0.22
0.02
0.15
12
Tail length
0.32
0.35
–0.10
–0.05
13
Anal body diameter
–0.24
0.37
0.17
0.22
Note: Key diagnostics for discriminating Helicotylenchus species are denoted in italic.
Molecular identification of nematodes
In relation to the D2 to D3 expansion segment of ribosomal DNA, four consensus sequences were obtained for H. multicinctus with similarity of 99% with different sequences of this species deposited in the NCBI (KF443214, DQ328745, DQ328746, HM014290, HM014291, and HM014292). On the other hand, two consensus sequences were obtained for H. dihystera, which presented similarity of 99% with other sequences of this species already deposited in NCBI (HM014251, KF486503, HM014250, HM014245, HM014246, HM014247, HM014248, HM014249, KM506834, KM506835, and KM506836). Molecular identity of the sequence obtained from the individual morphologically identified as H. erythrinae was not confirmed, because there are no reference sequences of H. erythrinae in NCBI or any other molecular database. When BLAST was performed, no sequence deposited in NCBI presented similarity greater than or equal to 99% with the sequence obtained in this research for H. erythrinae and closer sequences were of H. labiodiscinus with a similarity of 90.85 to 91.03% (HM014293; HM014294; HM014295). For individuals identified morphologically as H. californicus were obtained four sequences that showed similarity of 99% with an isolate of Helicotylenchus sp. labeled as CD761, with accession numbers KM506844 and KM506845. The sequences obtained in this study for individuals identified morphologically as H. erythrinae and H. californicus are proposed as standard and reference populations until topotype specimens become available and were molecularly characterized. All sequences obtained in the present study were deposited in NCBI with accession numbers MT321729-MT321739 (Table 7).
Information of sequences D2 to D3 of ribosomal DNA downloaded from GenBank and obtained in the present study for Helicotylenchus.
Note: Information of new sequences obtained in the present study are marked in italic.
Phylogenetic relationships of nematodes
The evolutionary relationships of 51 sequences of the D2 to D3 expansion segment of ribosomal DNA for Helicotylenchus, including those obtained in the present study (Table 7), are depicted in Figures 5 and 6. The sequences of H. multicinctus obtained in this research clustered in the same clade with other sequences of this species, including some populations from Musa sp. plantations, such as DQ328745 and DQ328746 from Sudan, HM014290 and HM014291 from Lambani, South Africa, and KF443214 from Fujian, China. Monophyly of this species was recovered both with maximum likelihood (bootstrap support or BS = 96%) and Bayesian inference (posterior probability of PP = 1). Likewise, the sequences of H. dihystera from Colombia obtained in this research, grouped in the same clade with other sequences of that species, including HM014248 and HM014250 from grasses in Hawaii and KF486503 from Musa sp. in Fujian, China (BS = 95%, PP = 1). The sequence of H. erythrinae grouped with the clade of H. labiodiscinus, although the support for this group was weak (BS = 63%, PP = 0.75). Finally, the sequences of H. californicus formed a clade with a sequence of Helicotylenchus sp. isolate CD761 from Calathea in USA (KM506844) (BS = 98%, PP = 1). These results support the presence of four species of Helicotylenchus in the Musa crops that were sampled.
Discussion
A total of 10 Helicotylenchus populations associated with Musa spp. crops were identified to the species level in the present study: four H. multicinctus, two H. dihystera, two H. californicus, and two H. erythrinae. All species identified in this study occur in Colombia, while only two species were registered in Brazil (H. multicinctus and H. californicus). Morphometric measurements recorded for these species closely resemble the type and reference populations (Golden, 1956; Sher, 1961; Van Den Berg and Heyns, 1975; Krall, 1990; Mizukubo et al., 1992; Wount and Yeates, 1994; Uzma et al., 2015).
Although it is reported that morphological and morphometric identification of species in Helicotylenchus is a difficult task because many species share very similar diagnostic characters and overlapping morphometrics, the populations of Helicotylenchus studied associated with Musa spp. were satisfactorily identified to the species level through morphological and morphometric evidence. PCA and conglomerate analysis clearly separated the four species identified, confirming the utility of morphological and morphometric data but also of multivariate statistical analysis to discriminate among species of the genus (Fortuner and Maggenti, 1991; Subbotin et al., 2015; Uzma et al., 2015). In accordance with PCA, the variables associated with the tail are diagnostic characters powerful enough to discriminate between or to separate among H. multicinctus, H. dihystera, H. californicus, and H. erythrinae, which could be useful for preparing a dichotomous key for identification of Helicotylenchus species in Colombia (Mizukubo et al., 1992; Uzma et al., 2015).
The diagnostic characters: number of tail annuli, ratio c′, maximum body diameter, tail length, anal body diameter, and vulva position separated the species identified in this study in accordance with PCA and conglomerate analysis. In various publications, tail length and DGO have been suggested by their discriminate values among Helicotylenchus species (Perry et al., 1959; Fortuner et al., 1984). Additionally, intraspecific variability at the morphometric level was observed in the four species identified in this study, due to the various morphological and morphometric unconstant characters with a coefficient of high variability within of Helicotylenchus species (Fortuner et al., 1981; Fortuner, 1984).
Different morphological and morphometric characters were constants within the species identified, as habitus, ratio a and V as reported for various species of Helicotylenchus (Fortuner, 1984). Variation in the shape of the tail was registered in females of the Colombian population of H. californicus. This species is characterized by showing individuals with short ventral projections and sharply pointed or blunt tails (Van Den Berg and Heyns, 1975; Krall, 1990). Variability in the shape of the tail to the intraspecific level has been recorded for various species of Helicotylenchus (Fortuner, 1979).
Molecular analysis supported the presence of H. multicinctus, H. dihystera, H. californicus, and H. erythrinae in Musa spp. crops from Colombia. Interestingly, in the phylogenetic tree, Colombian sequences of H. multicinctus and H. dihystera clustered in the same clade with sequences of nematodes isolated from banana crops in Africa and China, respectively (Subbotin et al., 2011; Xiao et al., 2014). The sequences obtained for H. californicus grouped with an isolate of Helicotylenchus from Calathea not identified to the species level and labeled as CD761 in Subbotin et al. (2015) (Table 7). Recently, Fortuner et al. (2018) identified isolate CD761 as H. pseudorobustus using morphological and morphometric analysis. However, our populations are similar to preliminary reports of H. californicus (Van Den Berg and Heyns, 1975; Krall, 1990). Unfortunately, there are no sequences of reference for H. erythrinae deposited in any of the databases of genes in the world for comparison. Therefore, the present study is reporting the first sequences of H. californicus and H. erythrinae associated to Musa spp. in a public database.
In relation to H. multicinctus, it is considered the most limiting species of Helicotylenchus in plantain and banana. It was present in the Musa crop production of both Colombia and Brazil, confirming the wide distribution of this plant-parasitic nematode, which also has been previously registered in Puerto Rico, Vietnam, India, Mexico, and Greece (Ngoc Chau et al., 1997; Ravichandra, 2014; Lara et al., 2016; Tzortzakakis et al., 2017). The presence of H. multicinctus in the different production zones, Quindío and Valle del Cauca, suggests a wide spread of this nematode in Colombia, directly related with the dispersion of contaminated seedlings between production zones, but also to environmental conditions that have favored the establishment of this species (altitude, soil type, and nutrient availability for the plant) (Nath et al., 1998; Speijer and Ssango, 1999; Araya et al., 2011; Godefroid et al., 2017).