Society of Nematologists
Subject: Life Sciences
ISSN: 0022-300X
eISSN: 2640-396X
SEARCH WITHIN CONTENT
Wenhao Li / Huixia Li * / Chunhui Ni / Deliang Peng / Yonggang Liu / Ning Luo / Xuefen Xu
Keywords : Goettingiana group, Heterodera , Morphology, New species, Phylogeny, Taxonomy
Citation Information : Journal of Nematology. Volume 52, Pages 1-16, DOI: https://doi.org/10.21307/jofnem-2020-097
License : (CC-BY-4.0)
Received Date : 05-July-2020 / Published Online: 24-November-2020
A new cyst-forming nematode,
Cyst-forming nematodes are the economical pests of cultivated crops and known to be reported from all the continents (Jones et al., 2013). The genus Heterodera was erected by Schmidt (1871) and currently contains about 80 species (Subbotin et al., 2010). Literature studies have indicated the presence of 14 Heterodera species from China mainland, including H. avenae (Chen et al., 1991), H. glycines (Liu et al., 1994), H. sinensis (Chen and Zheng, 1994), H. filipjevi (Li et al., 2010), H. koreana (Wang et al., 2012; Wang et al., 2012b), H. elachista (Ding et al., 2012), H. ripae (Wang et al., 2012a; Wang et al., 2012b), H. hainanensis (Zhuo et al., 2013), H. fengi (Wang et al., 2013), H. guangdongensis (Zhuo et al., 2014), H. zeae (Wu et al., 2017), H. sojae (Zhen et al., 2018), H. schachtii, and H. vallicola (Peng et al., 2020).
Due to overlapping morphological characters and phenotypic plasticity, it is difficult to distinguish closely related Heterodera species; therefore, sequence-based diagnosis is gaining more reliability for precise and accurate identification of cyst-forming nematodes (Peng et al., 2003). The internal transcribed spacer region of the ribosomal DNA (ITS-rDNA), the D2 and D3 expansion fragments of the 28S ribosomal DNA genes (D2-D3 of 28S-rDNA), and mitochondrial DNA (COI gene) units are good candidate genes for molecular taxonomic and phylogenetic studies (Subbotin et al., 2001; Subbotin et al., 2006; Madani et al., 2004; Vovlas et al., 2017). Based on morphomolecular characterizations, Handoo and Subbotin (2018) divided Heterodera into nine distinct groups such as Afenestrata, Avenae, Bifenestra, Cardiolata, Cyperi, Goettingiana, Humuli, Sacchari, and Schachtii group. Sequence analysis of these groups is significant to study the phylogenetic relationship and identifying the Heterodera species.
During 2018 and 2019, a population of cyst nematode was collected from the rhizosphere of Microula sikkimensis in Tianzhu county of Gansu Province, China. Considering the economic value of the cyst nematode, morphomolecular studies were performed; the preliminary studies indicated that the population belongs to Goettingiana group of Heterodera. The species characters were then compared with all the related species and concluded that this population possess unique characters and it is described herein as Heterodera microulae sp. n.
The nematodes were extracted from root and soil samples of Microula sikkimensis in Tianzhu county, Gansu Province, China. Cysts and white females were collected using sieving-decanting method, while second-stage juveniles (J2s) were recovered from hatched eggs and kept in water suspension until further use (Hooper, 1970; Golden, 1990). Males were not found. For morphometric studies, second-stage juveniles were killed by gentle heating, fixed in TAF solution (formalin: triethanolamine: water = 7:2:91), and processed to ethanol-glycerin dehydration according to Seinhorst (1959) as modified by De Grisse (1969) and mounted on permanent slides. Vulval cones were mostly mounted in glycerin jelly. Measurements were made on mounted specimens using a Nikon Eclipse E100 Microscope (Nikon, Tokyo, Japan). Light micrographs and illustrations were produced using a Zeiss Axio Scope A1 microscope (Zeiss, Jena, Germany) equipped with an AxioCam 105 color camera and Nikon YS 100 with a drawing tube (Nikon, Tokyo, Japan), respectively.
DNA samples were prepared according to Maria et al. (2018). Three sets of primers (synthesized by Tsingke Biotech Co. Ltd., Xi’an, China) were used in the PCR analyses to amplify sequences of the ITS, D2-D3 expansion segments of 28S, and COI gene. The ITS region was amplified with TW81 (5′-GTTTCCGTAGGTGAACCTGC-3′) and AB28 (5′-ATATGCTTAAGTTCAGCGGGT-3′) (Maafi et al., 2003). The 28S D2-D3 region was amplified with the D2A (5′-ACAAGTACCGTGAGGGAAAGTTG-3′) and D3B (5′-TCGGAAGGAACCAGCTACTA-3′) (De Ley et al., 2005; Ye et al., 2007). Finally, the partial COI gene was amplified using primers Het-coxiF (5′-TAGTTGATCGTAATTTTAATGG-3′) and Het-coxiR (5′-CCTAAAACATAATGAAAATGWGC-3′) (Subbotin, 2015). PCR conditions were as described by Ye et al. (2007), De Ley et al. (2005), and Subbotin (2015). PCR products were separated on 1% agarose gels and visualized by staining with ethidium bromide. PCR products of sufficiently high quality were purified for cloning and sequencing by Tsingke Biotech Co. Ltd., Xi’an, China. The PCR products were purified by the Tiangen Gel Extraction Kit (Tiangen Biotech Co. Ltd., Beijing, China), cloned into pMD18-T vectors and transformed into DH5α-competent cells, and then sequenced by Tsingke Biotech Co. Ltd (Xi’an, China).
The newly obtained sequences for each gene (ITS-rDNA, D2-D3 region of 28S-rDNA, and COI gene) were compared with known sequences of Heterodera using BLASTn homology search program. Outgroup taxa for phylogenetic analyses were selected based on the previously published studies (Subbotin et al., 2001; Maafi et al., 2003; Mundo-Ocampo et al., 2008; Kang et al., 2016; Madani et al., 2018; Vovlas et al., 2017). The selected sequences were aligned by MAFFT (Kazutaka and Standley, 2013) with default parameters and edited using Gblock (Castresana, 2000). Phylogenetic analyses were based on Bayesian inference (BI) using MrBayes 3.1.2 (Huelsenbeck and Ronquist, 2001). The GTR + I + G model was selected as the best-fit model of DNA evolution for both 28S D2-D3, ITS, and COI regions using MrModeltest version 2.3 (Nylander, 2004), according to the Akaike information criterion (AIC). BI analysis for each gene was initiated with a random starting tree and run with four Markov chains for 1,000,000 generations. The Markov chains were sampled at intervals of 100 generations and the burn-in value was 25%. Two runs were performed for each analysis. After discarding burn-in samples, the remaining samples were used to generate a 50% majority-rule consensus tree. Posterior probabilities (PP) were given on appropriate clades. The phylogenetic consensus trees were visualized using FigTree v.1.4.3 software (http://tree.bio.ed.ac.uk/software/figtree/) (Rambaut, 2016). The species in Goettingiana group and their localities, hosts, and GenBank accession numbers used in this study were presented in Table S1.
It is lemon-shaped with an obtuse vulval cone, neck extruding, and cuticle thick with an irregular zig-zag pattern. The color was white to pale to medium brown; remnants of the subcrystalline layer were rarely present. The egg sac was usually absent (Figs. 1G, 3B, C). The vulval cone was ambifenestrate-like waning crescent moon and separated by a well-developed vulval bridge. The anus area was distinct, bullae were absent (Figs. 1F, 3D, E). The vulval slit was longer than fenestral width (39.00 vs 37.75 µm); the underbridge was weak and often lost during cone preparation.
The female was lemon-shaped, pearl white, or pale yellow in color. It was rarely rounded with a protruding neck and vulva, the subcrystalline layer was present, and the egg sac absent (Figs. 2A, B, 3A). There was a labial region with two annuli. Labial sclerotization was weak, the stylet was strong, and basal knobs were rounded and anteriorly flattened. The excretory pore was indistinct, median bulb was rounded and massive, and other parts of the pharynx were not clearly discernable. There was vulval slit in a cleft on the cone terminus (Fig. 2C, D).
The body was straight or slightly curved ventrally after heat treatment (Fig. 4A). The lip region was offset and rounded, measuring 3.90 to 5.50 (4.63) µm in height and 9.65 to 12.75 (11.01)-µm wide. The cephalic framework was strongly sclerotized (Figs. 1B, 4C). The stylet was strong; knobs were well developed, rounded and flat, or slightly concave anteriorly (Figs. 1C, 4C). The dorsal esophageal gland orifice measured from 5.32 to 6.32 (5.61) µm posterior to the stylet knob. Median bulb was rounded with a strong valvular apparatus. The pharyngeal glands were well developed, overlapping the intestine dorsoventrally (Figs. 1A, 4B). The hemizoind was distinct from one to three annuli long (Fig. 4H), the excretory pore was situated 102.46 to 130.79 (114.40) µm from the anterior end, and one to two annules were posterior to the hemizonid (Fig. 4G). There was a lateral field with four incisures (Figs. 1D, 4G). The dorsal gland nucleus and subventral gland nuclei were distinct (Fig. 4E, F). Genital primordium situated at 59 to 62% of body length behind the anterior end, with two distinct nucleate cells (Fig. 4I). The tail was conoid, gradually tapering to a finely rounded terminus. The hyaline portion was irregularly annulated occupying 50% of tail length. Phasmid was absent (Figs. 1E, 4D).
Light micrographs of second-stage juvenile of H. microulae sp. n. A: Entire body; B: Anterior region of; C: Head region; D: Tail region; E: Posterior pharyngeal region arrow showing the position of dorsal gland nucleus; F: Posterior pharyngeal region arrow showing the position of subventral gland nuclei; G: Lateral field; H: Hemizonid; I: Genital primordium; J: Excretory pore (scale bar: A = 100 μm, B, H, G = 50 μm, C, D, E, F, I, J = 20 μm).
Holotype and paratype material (20 cysts, 20 females, and 20 second-stage juveniles) were deposited in the nematode collection of the Department of Plant Protection, Biocontrol Engineering Laboratory of Crop Diseases and Pests of Gansu Province, Lanzhou, China.
Heterodera microulae sp. n. was collected from the roots and rhizosphere soil of Microula sikkimensis Hemsl. (Boraginaceae, Tubiflorae, Metachlamydeae) in Tianzhu county of Gansu Province, China. The geographical position is N 37°11′46″; E 102°47′6″. This site was located in continental highland with the vegetation type of meadow grassland and the soil is composed of chernozems. The climatic parameters of the site include 450 mm of average rainfall and an approximate −2 air temperature.
Heterodera microulae sp. n. is characterized by having lemon-shaped cysts that have protruding necks and obtuse vulval cones. The cysts are 414 to 543-µm long and 305 to 456-µm wide having ambifenestrate vulval cone and bullae are absent. Females are white in color with a subcrystalline layer. Second-stage juveniles are straight or slightly curved ventrally with four incisures in the lateral field. The juveniles are 506 to 628-µm long having strong stylets with well-developed rounded stylet knobs, genital primordium situated at 59 to 62% of body length, and tail 49 to 61-µm long with a hyaline portion of 24 to 33 µm. Eggs are hyaline without any markings; juveniles inside the eggs form sixfold.
The new species belongs to the Goettingiana group of Heterodera; currently, the group contains seven valid species, viz, Heterodera goettingiana (Liebscher, 1892), H. carotae (Jones, 1950), H. cruciferae (Franklin, 1945), H. circeae (Subbotin and Turhan, 2004), H. scutellariae (Subbotin and Turhan, 2004), H. urticae (Cooper, 1955), and H. persica (Maafi et al., 2006).
The new species differs from H. goettingiana by having a shorter fenestral length (31 µm vs 35 µm), absence of bullae (vs few), weak underbridge (vs 117 µm), longer J2s body length (568 µm vs 486 µm), stylet knobs rounded and flat or slightly concave anteriorly vs smoothly rounded to slightly hook-shaped with a recurved anterior surface, longer distance of median bulb from the anterior end (MB) (86 µm vs 70 µm), shorter excretory pore distance from the anterior end (114 µm vs 158 µm), and shorter length of hyaline tail portion (29 µm vs 37 µm).
The new species is differentiated from H. carotae by having a bigger size of cysts (495 × 384 µm vs 408 × 309 µm), shorter vulval slit length (39 µm vs 47 µm), longer J2s body length (568 µm vs 422 µm), stylet knobs rounded and flat or slightly concave anteriorly vs concave anterior face, higher MB value (86 µm vs 66 µm), longer excretory pore distance from the anterior end (114 µm vs 99 µm), and longer tail length (57 µm vs 52 µm).
The new species differs from H. cruciferae by having a bigger size of cysts (495 × 384 µm vs 429 × 333 µm), slightly shorter fenestral length (31 µm vs 34 µm), shorter vulval length (39 µm vs 45 µm), longer J2s body length (568 µm vs 431 µm), higher MB value (86 µm vs 68 µm), longer excretory pore distance from the anterior end (114 µm vs 101 µm), longer tail length (57 µm vs 50 µm), and longer length of hyaline tail portion (29 µm vs 25 µm).
The new species differs from H. persica by a shorter fenestral length (31 µm vs 47 µm), absence of bullae (vs present), shorter vulval slit length (39 µm vs 49 µm), longer J2s body length (568 µm vs 440 µm), stylet knobs (flat or concave anteriorly vs projecting slightly anteriorly, convex posteriorly), longer stylet (26 µm vs 23 µm), higher MB value (86 µm vs 70 µm), longer excretory pore distance from the anterior end (114 µm vs 103 µm), longer tail length (57 µm vs 47 µm), and longer length of hyaline tail portion (29 µm vs 24 µm).
Compared with H. urticae, the new species has a smaller size of cysts (495 × 384 µm vs 492 × 435 µm), vulval cone obtrusive (vs unobtrusive) and absence of egg sac (vs presence), shorter fenestral length (31 µm vs 38 µm), shorter vulval slit length (39 µm vs 46 µm), longer J2s body length (568 µm vs 541 µm), shorter DGO (8 µm vs 5 µm), and shorter excretory pore distance from the anterior end (114 µm vs 130 µm).
The new species differs from H. circeae having a smaller size of cysts (495 × 384 µm vs 555 × 397 µm), a shorter fenestral length (31 µm vs 43 µm), vulval slit length (39 µm vs 48 µm), longer J2s body length (568 µm vs 434 µm), stylet knobs (rounded and slightly sloping posteriorly vs rounded and flat or slightly concave anteriorly), higher MB value (86 µm vs 70 µm), longer excretory pore distance from the anterior end (114 µm vs 101 µm), longer tail length (57 µm vs 52 µm), and longer length of hyaline tail portion (29 µm vs 26 µm).
The new species differs from H. scutellariae, having smaller cysts (495 × 384 µm vs 560 × 424 µm), by a shorter fenestral length (31 µm vs 35 µm), vulval slit length (39 µm vs 43 µm), longer J2s body length (568 µm vs 408 µm), higher MB value (86 µm vs 62 µm), longer excretory pore distance from the anterior end (114 µm vs 89 µm), longer tail length (57 µm vs 47 µm), and longer length of hyaline tail portion (29 µm vs 25 µm).
Additionally, comparative morphological and morphometric characters of H. microulae sp. n. with other valid species of Goettingiana group are listed in Table 2.
The H. microulae sp. n. sequences of D2-D3 region of 28 S (734 bp), ITS (993 bp), and COI (415 bp) gene were obtained and submitted to the GenBank.
The D2-D3 of 28S-rRNA sequence (accession no. MT573436) of H. microulae sp. n. showed 97.09% (19-bp difference), 97.66 to 98.49% (11-17-bp difference), 98.38% (9-bp difference), 98.62% (9-bp difference), 98.45% (11-bp difference), and 99.86 to 100% (0-1-bp difference) sequence identities with H. goettingiana (DQ328697), H. carotae (KX463292 and KX463293), H. cruciferae (KP114546), H. urticae (DQ328696), Heterodera sp. RH-2010 (GU456692) from Iran, and Heterodera sp. DP-2010 (HM560856 and HM560855) from Qinghai, China, respectively. The Bayesian phylogenetic tree of the D2-D3 of 28S gene (Fig. 5) represented a highly supported (posterior probability PP = 100) clade of Heterodera species, where Goettingiana group species occupied a basal position. It is noted that H. microulae sp. n. clustered together with Heterodera sp. DP-2010 (HM560855, HM560856) from Qinghai, China and forms a 100% supported clade.
Molecular phylogenetic tree of H. microulae sp. n. (highlighted in bold) inferred from 28S D2/D3 extension region under GTR + I + G model. The posterior probability values exceeding 50% are given on appropriate clades. *Identified as Heterodera sp. by Ye et al. (unpublished) and Peng et al. (unpublished) in the GenBank.
The ITS-rDNA sequence (accession no. MT573437) divergence of H. microulae sp. n. with other Goettingiana group species is as follows: 0.20% (2-bp difference), 0.4 to 0.5% (4-bp difference), 3.02% (29-bp difference), 5.01% (48-bp difference), 5.11% (49-bp difference), 7.45% (72-bp difference), 6.77 to 6.95% (67-68-bp difference), 6.29 to 7.25% (66-70-bp difference), and 7.41 to 8% (74-77-bp difference) for Heterodera sp. DP-2010 (HM560791), H. goettingiana (HM370423, HM370425), H. persica (AF498377), H. scutellariae (AY368995), H. circeae (AY368994), H. urticae (AF274412), H. carotae (AF274413; MG976790), H. cruciferae (AF274411; GU126668), and H. goettingiana (KY129827; AF274411; AF498374), respectively. The Bayesian phylogenetic tree of the ITS gene (Fig. 6) represented a highly supported (posterior probability PP = 100) clade of Heterodera species. As in the 28S tree, the ITS tree also positioned the Goettingiana group species. H. microulae sp. n. (MT573437) clustered with H. persica (AF498377), H. scutellariae (AY368994), H. circeae (AY368995), Heterodera sp. DP-2010 (HM560791), and H. goettingiana (HM370423, HM370425) from Qinghai, China with high-probability support (pp = 91%). It is also noted that sequences of H. goettingiana (HM370423, HM370425) from Qinghai, China, clustered outside with other H. goettingiana (KY129827, AF274411, and AF498374) subclades and should be considered a misidentification. However, H. microulae sp. n. (MT573437) is clustered with H. sp. DP-2010 (HM560791) and H. goettingiana (HM370423, HM370425) from Qinghai, China, with 100% support. It is also noted that H. microulae sp. n. (MT573437) clustered with two Chinese populations of Heterodera species (HM560791; HM370425) with 100% support.
Molecular phylogenetic tree of H. microulae sp. n. (highlighted in bold) inferred from ITS region under GTR + I + G model. The posterior probability values exceeding 50% are given on appropriate clades. *Identified as Heterodera goettingiana by Peng et al. (unpublished); **Identified as Heterodera sp. by Peng et al. (unpublished); ***Identified as Heterodera goettingiana by Huang et al. (unpublished) in the GenBank.
The COI gene sequence of H. microulae sp. n. showed 87.21 to 89.53% (differing from 36 to 44 bp), 88.19% (differing from 43 bp), 88.67 to 88.92% (differing from 46 to 47 bp), and 88.67 to 89.40% (differing from 44 to 47 bp), sequence identities with H. goettingiana (KY129829-KY129831), H. urticae (MK093155 and MK093156), H. cruciferae (MG563230 and MG563234), and H. carotae (KX463299-KX463306, MG563227, MG563229, MG563231-MG563233, and MN820659), respectively. The Bayesian phylogenetic tree of the COI gene (Fig. 7) represented a highly supported (posterior probability PP = 100) clade of Heterodera species. In this tree, H. microulae sp. n. clustered with H. goettingiana, H. urticae, H. cruciferae, and H. caratae with 98% support; however, H. microulae sp. n. formed a separate clade from those sequences.
Molecular phylogenetic tree of H. microulae sp. n. (highlighted in bold) inferred from COI gene under GTR + I + G model. The posterior probability values exceeding 50% are given on appropriate clades. *KC172916 identified as H. pratensis by Toumi et al. (2013) and later corrected to Heterodera carotae by Madani et al. (2018).
Taxonomy of Heterodera species has been revised extensively in the past; Baldwin and Mundo-Ocampo (1991) placed 23 Heterodera species into Goettingiana group. However, Sturhan (1998) and Subbotin et al. (2001) used J2’s lateral field characters and host preferences to separate Heterodera species into different groups (such as Bifenestra, Cyperi, and Humuli groups). The key morphological characters of the Goettingiana group include lemon-shaped cysts having a protruding neck, ambifenestration, and absence of bullae (small bullae occasionally present); some species may have an egg sac, vulval slit length > 35 µm, a thin vulval bridge, fenestral length (30-45 µm), and a weak underbridge. There were second-stage juveniles with body length > 400 µm, stylet length > 20 µm, tail length > 45 µm, hyaline tail portion > 20 µm, and lateral field with four lines (Subbotin and Turhan, 2004). The new species also belong to the Goettingiana group and morphologically very close to H. urticae; however, morphometrics of J2s body lengths, DGO and excretory pore position, fenestral length, vulval slit length, and cyst width can be used to differentiate both species.
Phylogenetically, it is evident that H. microulae sp. n. is a member of Goettingiana group. In our analyses, it is also noted that Heterodera sp. DP-2010 (HM560791, HM560855, and HM560856) and H. goettingiana (HM370423 and HM370425) from Qinghai, China, formed a well-supported molecular clade with the H. microulae sp. n. Moreover, the nucleotide differences of these sequences with our new species sequences are also very low (2-4-bp difference for ITS and 0-1 bp for 28S). Previously, Escobar-Avila et al. (2018) indicated that the sequences of H. goettingiana (HM370423 and HM370425) from Qinghai, China, might be a case of misidentification. Based on our phylogenetic and sequence analysis results, we regard Heterodera sp. DP-2010 (HM560791, HM560855, and HM560856) and H. goettingiana (HM370423 and HM370425) as H. microulae sp. n.
Heterodera microulae sp. n. is isolated from Microula sikkimensis, it is a biennial herbaceous plant that grows in forests, meadows, and forest edges at altitudes of 2,200 to 4,700 m, and it is widely distributed in South and East Asian countries (Pi et al., 2014). H. microulae sp. n. was found in Gansu and Qinghai Provinces, but we speculate that it is likely to be found in some localities that are characterized by low temperature, high rainfall, and high altitude.
The present study described a new species found in the rhizosphere of M. sikkimensis; further research is needed to understand the distribution and biology of the new species. In addition, plenty of leguminous crops (pea, kidney bean, pole bean, etc.) are growing in the same locality. Therefore, host-suitability tests of H. microulae sp. n. are an open research field to investigate the damage potential of this species.
This research was supported by the National Natural Science Foundation of China No. 31760507 and the Special R & D Fund for Public Benefit Agricultural, No. 201503114. The authors thank the assistance of the Institute of Plant Protection of China for the light micrographs and Dr. Maria Munawar from the Department of Biological Sciences, University of Lethbridge, Canada, for English polishing in this paper.
Figure 4:
Light micrographs of second-stage juvenile of
Figure 5:
Molecular phylogenetic tree of
Figure 6:
Molecular phylogenetic tree of
Figure 7:
Molecular phylogenetic tree of