Genus Laimaphelenchus has been defined by the presence of pedunculate tubercles that expanded to finger-like projections on the tail terminus, although some species with this character have been demonstrated to be polyphyletic and transferred to Aphelenchoides (Zhao et al., 2006a, 2007; Asghari et al., 2012; Carta et al., 2016; Maleita et al., 2018). At present, the genus contains 15 species distributed across different climatic zones and environment (Hunt, 1993; Peneva and Chipev, 1999; Asghari et al., 2012; Fang et al., 2019). Members of this genus are known to exhibit a global distribution as they have been reported from six continents (Hunt, 1993; Swart, 1997; Peneva and Chipev, 1999; Zhao et al., 2007; Negi et al., 2009; Pedram et al., 2018). Due to the presence of potential pest species in family Aphelenchoididae, members of this family are diagnosed with caution. However, none of the Laimaphelenchus species were reported to cause potential damage to conifers (Raghavendra and Newcombe, 2013), although their possible association with oak decline syndrome was suggested by Pedram et al. (2018).
During the present study, a population of Laimaphelenchus species was isolated from declining Chinese pine, Pinus tabuliformis Carrière, in Beijing, China, in November, 2018. The population was examined carefully and preliminary studies reveal the status of this species as a new species. Therefore, the objectives of the study are: to provide morphological and molecular characterization of L. sinensis n. sp.; and to demonstrate phylogenetic relationships of the new species with related aphelenchids.
Materials and methods
Nematode isolation and morphological study
Several twigs collected from declining Chinese pine (Pinus tabuliformis) were sliced into small pieces approximately 1 cm wide. The nematodes were isolated by the modified Baermann funnel technique for 24 hr. Adults for observation and measurements were collected from the declining twig samples as the cultures were unsuccessful. Permanent slides were prepared by heat-killed nematodes fixed with FA 4:1 and ethanol-glycerin dehydration according to Seinhorst (1959) as modified by De Grisse (1969). Morphometrics, drawings and light micrographs of nematodes were made with the aid of a Zeiss microscope equipped with a Zeiss AxioCam MRm CCD camera (Carl Zeiss Shanghai Co. Ltd. Shanghai, China).
Molecular and phylogenetic analyses
DNA samples were prepared according to Li et al. (2008). Four sets of primers (synthesized by Majorbio, Shanghai, China) were used in the PCR analyses to amplify the near full length 18 S, full length ITS region and D2-D3 expansion segments of the 28 S ribosomal RNA genes (rDNA). The near full length 18 S region was amplified as two partially overlapping fragments; for the first fragment, 988 F (5-CTC AAA GAT TAA GCC ATG C-3) and 1912R (5-TTT ACG GTC AGA ACT AGG G-3) were used and for the second fragment 1813F (5-CTG CGT GAG AGG TGA AAT-3) and 2646 R (5-GCT ACC TTG TTA CGA CTT TT-3) (Holterman et al., 2006). The full length ITS region was amplified with the forward primer TW81 (5’-GTT TCC GTA GGT GAA CCT GC-3’) and the reverse primer AB28 (5’-ATA TGC TTA AGT TCA GCG GGT-3’) (Joyce et al., 1994). The 28 S D2-D3 region was amplified with the forward primer D2A (5’-ACA AGT ACC GTG AGG GAA AGT TG-3’) and the reverse primer D3B (5’-TCG GAA GGA ACC AGC TAC TA-3’) (De Ley et al., 1999). PCR conditions were as described by Li et al. (2008) and Ye et al. (2007). PCR products were separated on 1.5% agarose gels and visualized by staining with ethidium bromide. PCR products of sufficiently high quality were sent for sequencing by Invitrogen, Shanghai, China.
The newly generated near full length 18 S and 28 S D2-D3 rDNA sequences of L. sinensis n. sp. were compared with other aphelenchid sequences available in GenBank using the BLAST homology search program (Altschul et al., 1990). The alignments of selected sequences were conducted with MAFFT (Katoh and Standley, 2013) with the default parameters and edited with AliView (Larsson, 2014). The best-fitted model of DNA evolution and the base frequency, the proportion of invariable sites and the gamma distribution shape parameters and substitution rates were obtained using jModelTest2 (Darriba et al., 2012) with the Akaike information criterion. The phylogenetic tree for each gene was obtained separately using MrBayes 3.2.3 (Ronquist and Huelsenbeck, 2003) with four chains (three heated and one cold). The number of generations for the total analysis was set to 1 × 107, with the chain sampled every 1,000 generations and the burn-in value set at 25%. The Markov chain Monte Carlo method within a Bayesian framework was used to estimate the posterior probabilities of the phylogenetic trees using the 50% majority rule (Larget and Simon, 1999). The consensus trees were selected to represent the phylogenetic relationships as well as the branch length and support level, all visualized using TreeGraph 2 (Stöver and Müller, 2010).
Results
Systematics
Laimaphelenchus sinensis n. sp.

Figure 1:
Line drawings of Laimaphelenchus sinensis n. sp. A: Entire female; B: Entire male; C: Anterior region; D: Female posterior region showing vulva and post-uterine sac; E: Lateral lines F, G: Female tail terminus; H: Male tail; I: Spicule. (Scale bars = A, B = 20 μm; C-I = 10 μm).

Figure 2:
Light photomicrographs of Laimaphelenchus sinensis n. sp. A: Entire female; B: Entire male; C: Lateral lines; D: Anterior region; E: Female posterior region showing vulva and post-uterine sac; F, G: Vulval regions; H: Female tail; I-K: Female tail terminus; L-N: Male tails arrows showing position of caudal papillae (Scale bars = A, B = 20 μm; C-N = 10 μm; Abbreviations: ex, excretory pore).
Measurements
Measurements of the new species are given in Table 1.
Table 1.
Morphometrics data for Laimaphelenchus sinensis n. sp.
| Female | Male | ||
|---|---|---|---|
| Character | Holotype | Paratypes | Paratypes |
| n | – | 6 | 5 |
| L | 914 | 968±46.1 (914-1064) | 876±69.1 (750-956) |
| a | 42.1 | 41.5±1.9 (38.6-44.8) | 46.2±2.6 (42.4-50) |
| b | 11.2 | 11.7±0.6 (11.2-12.7) | 10.5±0.6 (9.4-11.1) |
| b’ | 4.3 | 4.8±0.4 (4.3-5.6) | 4.6±0.3 (4.1-4.9) |
| c | 25.3 | 28±3.1 (25-32.9) | 19.9±1.4 (17.9-22.3) |
| c’ | 2.6 | 2.6±0.3 (2.1-2.9) | 2.7±0.2 (2.4-2.9) |
| V or T | 70.5 | 69.1±0.9 (67.7-70.7) | 70.8±4.4 (63.4-76.5) |
| Lip region height | 2.3 | 2.4±0.1 (2.2-2.6) | 2.7±0.1 (2.6-2.8) |
| Lip region width | 6.9 | 7.1±0.3 (6.8-7.4) | 7.1±0.6 (6.4-7.9) |
| Stylet length | 12.2 | 12.3±0.3 (11.8-12.6) | 12.3±0.6 (11.1-12.9) |
| Body diam. | 21.7 | 23.4±1.7 (21.6-26.4) | 19±1.1 (17.7-20.9) |
| Median bulb width | 18.9 | 13.2±0.7 (12.4-14) | 12±0.7 (11-12.8) |
| Median bulb length | 12.6 | 18.5±0.5 (17.6-18.9) | 17.4±0.8 (16.2-18) |
| Median bulb length/diam. ratio | 1.5 | 1.4±0.1 (1.3-1.5) | 1.5±0.1 (1.4-1.5) |
| Excretory pore from anterior end | 92.2 | 96.3±4.7 (88-104) | 83.2±2.5 (80.1-87) |
| Ovary length or testis | 430 | 454±24 (420-480) | 622.6±76.6 (476-690) |
| Post-uterine sac | 121 | 130.3±6.1 (119-138) | – |
| Vulva to anus distance | 270 | 264.2±20.9 (233-300) | – |
| Post-uterine sac length/vulva to anus (%) | 44.8 | 49.8±5.5 (39.7-55.8) | – |
| Anal (cloacal) body diameter | 14.1 | 13.3±0.7 (12-14.1) | 16.3±1.0 (15.3-18.2) |
| Tail length | 36.1 | 34.9±3.2 (29-38.6) | 44.1±1.8 (41.8-46.4) |
| Spicule (curved median line) | – | – | 14.0±0.6 (13.2-15) |
| Spicule (Chord) | – | – | 15.6±0.9 (14.1-16.6) |

