Year-Long Assessment of Soil Nematode Diversity and Root Inhibition-Indicator Nematode Genera in Rice Fields

Simple Summary Nematodes are the soil microbes that function for nutrient regulation and biological degradation. Nematode diversity changes during different seasons of a year due to changes in ecological factors. Soil characteristics change significantly due to lowland rice cultivation. This research investigated soil nematode diversity and seasonal changes in rice fields during three months of the year. Soil nematode abundance was different during the three seasons, and plant-parasite nematodes were more abundant during summer than in spring or winter. Soil characteristics, such as soil moisture, carbon content, and nitrogen content, were more common during the summer season than in the spring and winter seasons, while soil pH was low in the summer season. The plant-parasitic nematodes showed a stronger correlation with the soil characteristics during the summer season than in the spring and winter seasons. In addition, the abundance of some free-living nematode genera functioned as ecological indicators. The community and diversity indices of this study will help farmers and microbiologists in nematode management in crop fields. Abstract Soil nematodes contribute to nutrient cycling. This year-long study aimed to investigate the changes in the diversity of soil nematodes during the spring, summer, and winter seasons in rice fields at 24 sites and to determine the indicator nematode genera that inhibit the roots of rice plants. A total of 216 soil samples were collected during three seasons, and the collection of 72 root samples was carried out during rice cropping. Forty-four soil nematode genera were identified. They exhibited significant changers in their abundance, which were dependent on the seasons and on soil characteristics. In particular, the abundance of plant-parasitic nematodes (PPNs) and free-living soil nematodes was 49% and 15% higher during the summer than during the spring and winter seasons, respectively. Soil characteristics, such as soil nitrogen (N) contents, carbon (C) contents, and soil moisture were significantly higher during the summer than in the spring and winter seasons, but soil pH was significantly lower during the summer than in the spring and winter seasons. Moreover, Hirschmanniella, Meloidogyne, and Heterodera emerged as good indicators for rice root inhibition, corroborating the frequency, density, and prominence value of PPNs of the sampled soil and rice roots. This study also indicated that free-living nematode genera, such as Rhabdolaimus, Diplogaster, and Rhabditis, might function as ecological indicators for soil health.


Introduction
More than half of the world's population depends on rice as the primary source of food, most of which (>90%) is cultivated and consumed in Asia [1]. Nematodes parasitize most crop plants and cause over USD 150 billion of crop loss per year globally [2,3]. In latitude 33 • 38 24 N and longitude 75 • 01 12 E (Figure 1). These selected sites were used for rice crops for many years. The sites were selected based on rice planting times, irrigation systems, and management of fertilization by nitrogen (N), phosphorus (P), and potassium (K) (NPK) (N = 120, P = 40, and K = 40 kg hac −1 ). A total of 216 independent soil samples were collected during three seasons of the year: spring, summer, and winter. The soil sampling in the spring season was carried out before rice plantation, on 12 April 2018 to 14 April 2018; the summer sampling was carried out when the rice plants were at the maximum tiller stage, on 20 July 2018 to 22 July 2018; the winter sampling was carried out after the rice crop was harvested, on 27 November 2018 to 29 November 2018. Of the soil samples collected, 72 were from each season. At each of 24 selected sites, three random soil samples were collected during each sampling time. For each sample, 500 g of soil were taken, with the help of a soil probe, from 5 cm to 10 cm below ground surface and kept separately in properly labeled plastic bags. The collected soil samples were stored at 4 • C during transfer to the lab for analysis, to maintain moisture. The rice root samples were collected during the summer sampling from the same 24 selected sites of the rice fields. A total of 72 independent rice root samples were collected from the 24 sites, and at each site, three random root samples were obtained. Fifty grams of the roots were collected from several plants and kept in labeled plastic bags for each sample.

Measurement of Soil pH and Moisture Content
A sub-sample of 20 g of soil from each sample was dissolved in 100 mL of distilled water separately and stirred with a glass rod till the formation of uniform soil suspension. An automatic pH meter was used to determine the pH of the soil suspension. The gravimetric method was employed to determine the soil moisture content. A sub-sample of 10 g of soil was used to measure fresh weight; then, the soil samples were dried in an oven for 24 h at 105 • C before measuring dry weight. The soil moisture contents are repre-sented in percentages using a formula: fresh weight of the sample-the dry weight of the sample × 100/dry weight of the sample.

Determination of Soil Nitrogen and Carbon Content
A CHNS analyzer (EuroVector 3000, Pavia, Lombardy, Italy) measured the soil N and C content based a simultaneous determination by a gas chromatography (GC) system. The instrument was calibrated, and a standard curve was prepared by using certified standards. A sub-sample of 2 g of soil from each sample was oven-dried at 105 • C for 24 h, then finely powdered with the help of a ball mill machine. After that, 10-12 mg of powdered soil was used for analysis, which was completed within 10 min after loading to a combustion reactor [33].

Nematode Extraction from the Soil Samples
For nematode extraction, a sub-sample of 200 g of soil was taken from each sample, mixed well to obtain a uniform distribution of nematodes, and mixed thoroughly with water to obtain soil suspension. The soil suspension was passed through a coarse sieve to remove debris and lumps. The suspension was also passed onto the fine sieve to obtain minor soil particles. Finally, 50 mL volume of soil suspension was used for nematode extraction. Nematodes were isolated sequentially using modified Cobb's sieving and decantation (with mesh sieves of 833, 74, and 38 µm apertures) and Baermann's funnel technique [34].

Nematode Extraction from the Root Samples
The nematodes of the rice root samples were extracted by a maceration sieving method [35]. The root samples (50 g) were washed five times, sheared into 2 cm pieces, then air-dried. The root samples were softened for 10 s, crushed in a blender, then transferred to Baermann's funnel. The succeeding steps were similar to those of nematode extraction from the soil samples.

Identification of Nematode Genera
Glass slides were prepared for the identification of nematodes. Using an Olympus microscope, BX51 nematodes were identified up to the genus level, based on the morphological characters of nematode genera at 100× and 40× magnifications. Trophic groups were allocated according to Yeates [36] and arranged in colonizer-persister (c-p) groups according to Bongers [27]. The counting of nematodes was performed using a Syracuse counting dish. Identification of nematodes up to a generic level was carried out on the basis of appropriate literature [37][38][39]. Hunt [40] taxonomic keys for plant-feeding dorylaimids (Trichodorids and Longidorids) and Aphelenchids were also used.

Nematode Community and Diversity Indices
The relationship between nematodes and rice fields was determined by community analysis. Nine nematode indices were calculated, showing different aspects of the community. The Shannon-Wiener index (H') and the Simpson index (D) were calculated, which indicated the nematode species diversity in a community [41,42]. Both the Shannon-Wiener index and the Simpson index accounted for the abundance and evenness of the species present. The species evenness was illustrated by Pielou's evenness index (J'), which is closely related to species dominance [43]. Species richness was measured by the Margalef index (MgI) and was calculated as follows: MgI = (G−1)/ln (n), where G is the total genera number and n is the total number of individuals [44].
Environmental disturbances in the soil were signified by the maturity index (MI) [45,46], based on the relative contribution of nematodes of different c-p values. The MI considered all nematode groups (c-p 1 to 5), and the high MI revealed stable soil conditions. The functional structure of the community was referred to by the Wasilewska index (WI), the channel Index (CI), the plant-parasitic index (PPI), and food web complexity (FWC). The WI represents the ratio of bacterial feeders (BF) plus fungal feeders (FF) to plant-parasites (PP) as WI = (BF + FF)/PP [47]. The CI represents the fungal participation in decomposition channels of soil food webs. The higher value signifies dominated fungal feeding decomposition, whereas the low values indicate a dominated bacterial decomposition pathway [14]. The PPI is similar to MI, but only for plant-parasitic nematodes, i.e., 1/N (c − p)i × ni, where (c − p)i is the c-p value for plant-parasitic nematodes and ni is the total number of individuals of a plant-feeding nematode i [48]. The FWC is the ratio between predatory nematodes (PR) and plant-feeding nematodes. It represents the top-down control in which predatory nematodes control herbivores, i.e., plant-feeding nematodes (PR/PF) [49]. The functional indices were calculated using the online program "NINJA: An automated calculation system for nematode-based biological monitoring" [50], https://sieriebriennikov.shinyapps.io/ninja/) (accessed on 5 June 2019).

Statistical Analysis
The datasets were analyzed by a one-way ANOVA model where pH, moisture, nitrogen content of the soil, carbon content of the soil, and soil nematode measures function as responses and seasons function as factors. The absolute frequency (AF), relative frequency (RF), mean density (MD), and relative density (RD) were calculated from the mean data of nematode soil samples. The N denotes the number of samples in which the genus was present, and the AF denotes the frequency of the genus by the total number of soil/root samples. MD is the ratio between the number of nematode specimens counted in all samples and the total number of samples collected. RD represents the percentage of the mean density of the genus × 100/sum of the mean density of all nematode genera. The prominence value (PV), which represents the dominant genus, was calculated as PV = population density √ relative frequency, and the relative prominence value (RPV) was calculated as the ratio between the PV of the nematode genus and the sum of the PVs of all nematode genera ×100. The PV and RPV indicate the dominant genus and a relationship of population density and frequency of the identified nematode genera. Analysis of variance was performed by Minitab 19.0, and Duncan's multiple range test was applied separately for each season. A principal component analysis (PCA) was used on the nematode-community data for the three-season samplings, with explanatory variables of soil pH and moisture. Correlation between nematodes found in root and soil samples was performed to determine indicator plant-parasitic genera.

Results
The variation of soil factors (pH, moisture, nitrogen content, and carbon content) and nematode indices (density, frequency, and abundance of soil nematodes) were significant between the three seasons of the year in the rice fields ( Table 1). The highest variation was depicted by soil pH (66.12%) and moisture (58.23%), followed by soil nitrogen content (48.34%). Among the soil nematodes that were measured, abundance showed the highest variation of 55.25% and frequency showed the lowest variation of 44.61%. In addition, the root nematode frequency depicted significant interaction with the soil nematode frequency of the three seasons in the rice fields. Table 1. One-way ANOVA analysis of soil ecology-pH, moisture, nitrogen content, and carbon content; soil nematode traits-density, frequency, and abundance; and interaction between soil and root nematodes.

Soil pH, Moisture, Nitrogen, and Carbon Content in Different Seasons
There was a significant change in soil pH, moisture content, total nitrogen content, and soil carbon content between the three seasons of the year in the rice fields ( Figure 2). The highest soil pH (6.80) was observed in spring, followed by the winter and summer seasons, whereas the highest soil moisture content was during summer (87.50%), followed by the winter (50.10%) and spring (31.20%) seasons. The soil nitrogen content (24.21 mg g −1 soil) and carbon content (34.56 mg g −1 soil) were highest in summer, compared with the winter and spring seasons.

Abundance of Soil Nematodes in Different Seasons
A significant change was observed in nematode abundance during the three seasons of the year ( Figure 3). The nematode abundance was highest in summer, followed by spring and winter. The contribution of PPNs was observed to be highest in the summer season and lowest in the spring season. In contrast, the highest contribution of free-living nematodes was found in the spring and the lowest was found in the summer season in the rice fields.

Trophic Nematode Structure in Different Seasons
A diversified nematode population was found during the different seasons of the year ( Figure 4). Bacterivorous nematodes and PPNs dominated over other nematode groups during the spring and summer seasons, respectively. Nevertheless, they were equally dominant in the rice fields during the winter season. A minor abundance of nematodes was observed in the predatory group during each of the year's three seasons.

Variation of Soil Nematode Frequency and Density in Different Seasons
The frequency and density of nematode genera varied frequently in three seasons of the year ( Figure 5). During spring, the genus Diplogaster had the highest frequency (17) and density (112). During the summer season, the genus Hirschmanniella had the highest frequency (23) and the genus Meloidogyne had the highest density (130). During the winter season, the genus Dorylaimellus had a higher frequency (18) and density (116) among the 44 nematode genera found in the soil samples. Similar trends were observed in other nematode diversity measures, such as absolute frequency, relative frequency, mean density, and relative density (see the Supplementary Materials). A correlation between the density and the frequency of the nematode genera of soil samples changed with the seasons in the rice fields. During the winter season, the density and frequency of the nematode genera predicted a strong correlation (R 2 > 0.9117; p = 0.021), which decreased in spring, but little correlation (R 2 = 0.8213; p = 0.035) was observed during the summer season in the fields.

Shift of PPNs during Different Seasons in the Rice Fields
The PPNs showed dynamic variation with the seasonal changes of the year in the rice fields ( Figure 6). Genera such as Longidorus, Meloidogyne, Hirschmanniella, Tylenchus, and Paratylenchus were persistent during the three seasons of the year. These genera displayed more than a two-fold increase in frequency and density during the summer season, and more than a one-fold increase during the winter season, compared with the spring season. PPN genera such as Rotylenchus, Hoplolaimus, Hexatylenchus, and Xiphinema appeared only in the summer season in the rice fields. A correlation between the density and the frequency of PPN was strong during the summer season (R 2 = 0.9632) and weak during the winter season (R 2 = 0.7318). The changes in plant-parasitic nematode diversity of the soil during the spring, summer, and winter seasons in the rice fields. Linear trendlines were obtained from mean data of the frequency and the density of the soil PPN genera; R 2 represents the correlation strength in the three seasons.

PV and RPV of Nematode Genera in the Rice Fields
The soil nematode genera prominence and relative prominence were diverse during the three seasons of the year in the rice fields ( Table 2). The highest prominence value (PV = 459) was observed during spring in the genus Rhabdolaimus, while the least PV was found in the genus Clarkus (PV = 11.95). Among the PPNs, the genus Meloidogyne had the maximum PV (351.5), followed by the genus Ditylenchus (PV = 325.5), while the least PV was noted in the genus Xiphinema (PV = 27.8). * Data represent the mean value ± standard error. The nematode genera with higher persistent prominence values are denoted with an asterisk (*) and lost persistence during the summer season is denoted with a dollar sign ($).
During the spring season, the genus Rhabdolaimus had the highest relative prominence value (RPV = 6.4), followed by the genus Diploscapter (RPV = 5.7), while the least RPV was found in the genera Acrobelus and Protorhabditis (RPV = 0.2). During the summer season, the highest RPV (87.8) was found in the genus Meloidogyne, followed by the genus Ditylenchus (RPV = 81.2), and the lowest (RPV = 2.8) was found in the genus Miconchulus. The genus Dorylaimoides displayed the highest relative prominence value (RPV = 82.5) during the winter season, while the lowest (PRV = 3.4) was found in the genus Longidorus.

Principal Component Analysis of Soil Nematodes
The nematode-community data was used for the PCA, which determined the favorable position of nematode diversity during the spring followed by the summer and winter seasons of the year (Figure 7). Four groups of nematodes were determined in the PCA plot, based on PC 1 and PC 2 contributions. Group-I mainly consisted of bacterivores; group II mainly consisted of omnivores, bacterivores, and PPNs; group III mainly consisted of PPNs; and group IV mainly consisted of predatory and fungivores. According to the position of the groups, PPNs had a strong correlation with bacterivores during the change of seasons durng the year. At the same time, the predatory and the fungivores had a strong correlation with the omnivores and some genera of the bacterivores and PPNs during the different seasons of the year.

Principal Component Analysis of Soil Nematodes
The nematode-community data was used for the PCA, which determined the favorable position of nematode diversity during the spring followed by the summer and winter seasons of the year (Figure 7). Four groups of nematodes were determined in the PCA plot, based on PC 1 and PC 2 contributions. Group-I mainly consisted of bacterivores; group II mainly consisted of omnivores, bacterivores, and PPNs; group III mainly consisted of PPNs; and group IV mainly consisted of predatory and fungivores. According to the position of the groups, PPNs had a strong correlation with bacterivores during the change of seasons during the year. At the same time, the predatory and the fungivores had a strong correlation with the omnivores and some genera of the bacterivores and PPNs during the different seasons of the year.

Nematode Genera Inhibited Rice Roots
Seven parasitic nematode genera were found in the root samples of the rice plants ( Figure 8). The nematode genera Hirschmanniella showed the highest abundance in root samples, followed by Meloidogyne; the abundant minor genus was Xiphinema. Three rootnematode genera with high frequency and density, such as Hirschmanniella (F = 21), Meloidogyne (F = 19), and Heterodera (F = 15), were indicators of the nematode community. The effects caused by the dominant root nematodes on the rice root were root-knots, cysts, and lesion formations (Figure 8c-e). Among the root nematodes, Hirschmanniella had the highest (4.87) relative frequency, followed by Meloidogyne (4.41) and Heterodera (3.21), while Xiphinema (0.58) had the lowest. Similar trends were observed in absolute frequency, mean density, relative density, prominence value, and relative prominence value of the root-nematode genera in the rice root samples (Table 3).

Nematode Genera Inhibited Rice Roots
Seven parasitic nematode genera were found in the root samples of the rice plants ( Figure 8). The nematode genera Hirschmanniella showed the highest abundance in root samples, followed by Meloidogyne; the abundant minor genus was Xiphinema. Three root-nematode genera with high frequency and density, such as Hirschmanniella (F = 21), Meloidogyne (F = 19), and Heterodera (F = 15), were indicators of the nematode community. The effects caused by the dominant root nematodes on the rice root were root-knots, cysts, and lesion formations (Figure 8c-e). Among the root nematodes, Hirschmanniella had the highest (4.87) relative frequency, followed by Meloidogyne (4.41) and Heterodera (3.21), while Xiphinema (0.58) had the lowest. Similar trends were observed in absolute frequency, mean density, relative density, prominence value, and relative prominence value of the root-nematode genera in the rice root samples (Table 3).  Table 3. The absolute and relative frequency, mean and relative density and prominence value and relative prominence value of plant-parasitic nematode genera found in rice root samples.

Community and Diversity Indices of Soil Nematodes
The community and diversity indices of the soil nematodes in the spring, summer, and winter seasons explained the soil structure and stability of the rice fields (Table 4). The highest H' was observed in summer (3.60) and the lowest in winter (3.40). The J' (1.00), MgI (6.21), and CI (9.10) mainly dominated in the summer season in the rice fields. The highest WI (1.50) was observed during the winter season in the rice fields. In addition, the D (0.97), FWC (0.60), and MI (2.40) of the soil nematodes were most significant during the spring season in the rice fields. The highest (2.97) PPI of the soil nematodes was found in summer, while the lowest (2.47) PPI was found in winter.

Discussion
In the present study, nematode diversity was examined in Indian rice fields during the spring, summer, and winter seasons. It was found that the bacterivorous nematodes dominated at higher pH. In contrast, the plant-parasitic nematodes dominated in acidic pH of the soil and decreased the bacterial feeders during all three seasons in the rice fields (Figures 2 and 4). However, omnivores were not much affected by the change of soil pH during the seasons. A previous study found that low pH decreases bacterivorous nematodes but increases fungivorous nematodes [51]. Factors such as soil pH and soil moisture significantly impact the nematode community, and their changes result in a change of the nematode community structure [13,52,53]. The habitat of nematodes [54] and their density, frequency, and diversity change with changing edaphic and ecological factors [55,56]. The pH change in the rice fields and nematode abundance may occur due to inorganic fertilizers [57][58][59][60]. In addition, soil moisture plays an essential role in the distribution of nematodes. The PPNs were highest with higher soil moisture content, followed by predatory and bacterivorous nematodes. Thus, soil moisture is a regulating factor for nematode abundance [61]. In some instances, the low soil moisture causes the reproduction rate in nematodes to decline [62], as it functions as a medium during interchange [63]. Based on the above observations and statements, nematodes are indicators for community analysis, soil health, and water conditions for soil microbes in changing environmental habitats [13,30,64].
During three seasons, the nematode dominance varied, and free-living nematodes were dominant in the spring and winter seasons; however, PPNs were dominant in the rice fields during the summer season. Nematode frequency and density were lower during summer (r 2 = 0.8213) than in the spring (r 2 = 0.9117) and winter (r 2 = 0.9269) seasons in the soil samples ( Figure 5). The decrease in free-living nematode populations during summer may be due to the low adaptation of these nematodes under flood irrigation conditions (anaerobic conditions) in the rice fields. Among the five trophic groups of 44 identified soil nematode genera, PPNs constituted the highest genera in terms of abundance in the rice fields. These results coincide with those of previous studies [65]. The higher population of plant-parasitic nematodes during the summer rice crop may be due to the growing crop, manures, and fertilizers [58,59]. Among the PPN genera, Hirschmanniella, Meloidogyne, Heterodera, Tylenchus, Paratylenchus, and Longidorus had high frequency and density during all three seasons of the year. However, regarding the frequency of the genera, Hirschmanniella, Meloidogyne, and Heterodera increased more than two-fold during the summer rice crop than the spring season (with no crop) ( Figure 6). The dominance of these three nematode genera was previously reported in rice fields worldwide [66,67]. Among the free-living nematode genera, Diplogaster, Rhabditis, Mesodorylaimus, Dorylaimoides, and Aphelenchus may function as indicators for community structure and soil stability in rice fields, due to their persistence during all three seasons ( Table 2).
The three nematode infections (root-knots, cysts, and lesion formations) associated with seven PPN genera were found in the roots of rice plants in this experiment. Rootknot infections in rice was mostly caused by Meloidogyne graminicola; cysts were caused by the Heterodera species; and lesion formations were caused by the Pratylenchus species. Of seven identified root-nematode genera, Hirschmanniella, Meloidogyne, and Heterodera dominate rice roots. The frequency of the three root nematode genera approached the frequency of nematodes in the soil during summer; in the spring, it was was < 0.5-fold (Figures 6 and 8). These results suggest that the frequency and density of Hirschmanniella, Meloidogyne, and Heterodera function as the best indicators for the nematode inhibition rate in the rice roots and the position of PPNs in the soil community structure. In previous reports, Meloidogyne was the most prevalent and abundant nematode in flooded and rainy lowland areas [68][69][70][71]. It has worldwide dispersion [67] and its spread potential to different topographical places [72] creates an in-time alarm of possible crop destruction in the future. Similarly, Hirschmanniella and Heterodera have great potential to spread to other ranges of topography, especially in rice fields, and to harm agriculture [6,73].
The H' was highest during summer (H = 3.6), which suggests high diversity. The H' is used to characterize species diversity in a community and helps find the disturbances of the habitat [74]. Meanwhile, the J was also highest during summer, which signifies greater evenness because of resource distribution [43]. Other indices, such as the MgI, PPI, and CI of soil nematodes, were significant in the rice fields during summer. The value of MgI was 6.12 during summer, which explains the high species richness in a community [75], which may be due to fertilizers in the rice fields. The PPI was highest (2.97) in the soil nematodes during summer; this refers to dominance and better adaptation under flood irrigation. The CI was a good indicator of decomposition [58,76]. It was highest (9.1) during summer and lowest during spring; the lower values explain the dominant bacterial pathway in the spring season of the rice fields. The FWC and MI of the soil nematodes were more significant during spring in the rice fields. The decrease of FWC during summer may be due to the abundance of PPNs in the soil. The MI value shows disturbance in the soil [48], and it was found to be lowest (2.2) during the winter season in the rice fields, which is due to the dominance of bacterial feeders. The decreased MI suggests the nematode community's decreasing structure, because the FWC decreased with increasing N deposition in the soil during summer. The soil mineralization process was indicated by WI and explained the relative balance of positive-to-negative impacts of nematodes on primary productivity [77,78], which was highest (1.50) during the winter season. Due to agricultural practices, immense trophic diversity was mainly associated with an increased frequency of less abundant groups, such as predator groups, fungivorous groups, and omnivore groups [49].
Rice cultivation continuously alters the diversity and abundances with different feeding habits, and thus alters the overall complexity and architecture of the detritus food web, which coincides with the findings of Korobushkin [79]. The above-mentioned community and diversity indices of soil nematodes are helpful in investigating soil stability and efficiency, explaining life strategies to exist as colonizer-persisters, and feeding pathways. Some factors regulate nematode diversity indices; for example, the application of manures and fertilizers [59,80].

Conclusions
The composition and structure of nematodes exhibited apparent diversity in rice fields during three seasons. The pH of the soil was not the same during the three seasons, and higher nematode abundance was observed, at slightly acidic pH, during the summer. The acidic nature of the soil during the summer decreased free-living nematodes and, in turn, increased PPNs and possibly affected the growing crop. A disruption of nematode trophic communities during the summer occurred because of the amplification of the PPNs. Hirschmanniella, Meloidogyne, and Heterodera were identified as the best indicators for rice root inhibition and soil health. These PPN genera could also function as indicators for crop management in rice fields, because they exhibited the same frequency rate in the rice roots and the soil samples during the summer season. In addition, Diplogaster, Rhabditis, and Aphelenchus were important ecological balancing indicators for soil health, due to their persistence in the soil during all three seasons. The community and diversity indices of this study will assist farmers and microbiologists in nematode management in rice fields.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/biology11111572/s1, Table S1: Datasets of frequency and density and their relative data of nematode genera during the three seasons in the soil of the rice fields.