Skip to main content
Have a personal or library account? Click to login
Nematodes and the effect of seasonality in grassland habitats of South Africa Cover

Nematodes and the effect of seasonality in grassland habitats of South Africa

Journal of Nematology
Volume 52 (2020): Issue 1 (January 2020)
By: ,  ,  ,   and    
Open Access
|Dec 2020
Article

Full Article

Nematodes occupy most terrestrial habitats on earth (Liu et al., 2019), even the furthest reaches of caves (Du Preez et al., 2017) and deep underground mines (Borgonie et al., 2011). Estimates also indicate that nematodes represent 80% of all multicellular organisms (Eisenhauer and Guerra, 2019; Van Den Hoogen et al., 2019). However, despite their omnipresent distribution and dominating abundance, many nematode communities are poorly studied with the majority of species remaining undescribed (Archidona-Yuste et al., 2016). This is especially true for many habitats in South Africa since most nematode-related studies are focused on agricultural systems and nematodes of economic importance (i.e. crop pests) (Procter, 2012; Fourie et al., 2017).

Grasslands, for example, host diverse soil ecosystems, represent more than 40% of Earth’s terrestrial surface area and provide essential ecosystem services including the provision of food, storage of carbon, mitigating droughts and floods, and erosion control (Hewins et al., 2018; Bengtsson et al., 2019; Gibson and Newman, 2019). However, although the structure and ecological status of grassland nematode communities have been investigated in some parts of the world (Ekschmitt et al., 2001; Biederman and Boutton, 2009; Kergunteuil et al., 2016), a literature survey revealed no ecological reports on the nematode communities of South African grasslands.

Studying the nematode communities of grasslands is not only relevant from an ecological perspective, but also from a conservation perspective. Nematodes are valuable indicators of soil ecosystem disturbance affected by, for example, agricultural activities and climate change (Ferris and Bongers, 2009; Zhong et al., 2017; Zhang et al., 2020). Therefore, studying the nematode communities of natural, undisturbed grasslands will generate valuable baseline data for monitoring disturbance and implementing conservation policies. This would be especially valuable for South Africa as the country’s grasslands are under continuous threat from overgrazing, cultivation, and urban expansion (Haddad and Butler, 2018). Woody (or bush) encroachment, driven by rising temperatures and increasing atmospheric carbon dioxide levels, also presents a major threat (Skowno et al., 2017). This while only 2.0 to 2.8% of South Africa’s grasslands are protected within conserved areas (Carbutt et al., 2011; Haddad and Butler, 2018).

When studying terrestrial ecosystems, it is also important to consider seasonal variation as water availability, temperature, and vegetative growth can have an important effect on the structure and functioning of soil ecosystems (Bongers and Ferris, 1999; Cardoso et al., 2016; Song et al., 2016; Siebert et al., 2020). Song et al. (2016), for example, showed that higher precipitation rates significantly increased nematode abundance and richness in a grassland ecosystem. Bakonyi et al. (2007), in turn, found that the effects of soil drying and warming on nematode communities were dependant on the presence of vegetation with the smallest effect recorded in grasses. Ultimately, understanding how seasonal variation affects soil ecosystems is necessary in order to accurately assess and monitor the potential threats posed by anthropogenic activities and climate change.

For these reasons, this study was undertaken and aimed at (i) studying the nematode community structure in the open grassland, shrubland with rocky outcrops and riparian zone habitats of the Telperion Nature Reserve and (ii) investigating the effect of seasonal-induced changes on the condition, food web status and functioning of the soil ecosystems.

Materials and methods

Site description

The study was conducted in the Telperion Nature Reserve (Mpumalanga province, South Africa) (Fig. 1) that has a total surface area of 9,061.3 ha. This reserve is located in the Rand Highveld Grassland (of which only 1% is conserved) and is dominated by grassy plains interspersed with rocky outcrops and woody plant species (Grobler, 1999). The greater region is characterized by strong summer rainfall (570-730 mm per annum) and dry and cold winters (frost may occur frequently). This is a fire prone ecosystem with high lightning flash densities making lightning-induced fires common. The lithology of Telperion Nature Reserve is dominated by Arenite-Conglomerate, producing dystrophic or mesotrophic soils with some red soils, as well as rocky areas with miscellaneous soils (Grobler, 1999). Three streams flow through the reserve that originates from higher lying wetlands and sponge areas.

Figure 1:

Satellite image of Telperion Nature Reserve (South Africa) indicating the sampling locations associated with the open grassland, shrubland with rocky outcrops and riparian zone habitats.

Rainfall and temperature data

Climatic data (rainfall and temperature) from January 2014 to December 2016 were obtained from the South African Weather Service (www.weathersa.co.za). The total amount of rain for 90 days prior to each sampling event (as listed in the following section) was calculated, while the mean, minimum, and maximum daily temperatures were extracted for the same period. Since sampling was based on climatic factors (as listed in the following section) and thus not evenly dated over 1 year, this 90-day period was considered before each sampling event.

Site selection and sample collections

Within Telperion Nature Reserve three main terrestrial habitats, namely, (i) open grassland (OG), (ii) shrubland with rocky outcrops (SRO), and (iii) riparian zone (RZ) were identified for investigation. Sites within each habitat were selected based on accessibility as some areas were restricted and/or not accessible by vehicle or foot. The elevation of the sites in the each habitat ranged as follows: OG: 1333 to 1482 m, SRO: 1361 to 1480 m, and RZ: 1304 to 1383 m above sea level. Ultimately, 26 sampling sites (Fig. 1) were selected and included: 11 OG, 10 SRO sites, and 5 RZ sites.

In total, 104 soil samples (one per site per season) were collected from the three terrestrial habitats during four consecutive seasons [Winter (11 June 2015), Spring (30 November 2015), Summer (24 February 2016), and Autumn (19 April 2016)]. These sampling dates were selected based on climatic patterns through consultation with the reserve manager. Winter samples were the first to be collected when frost was at its peak. Spring sampling followed the first rains, which due to a drought in South Africa came later than normally expected. Summer sampling was undertaken when grasses reached maturity and produced seeds, while Autumn samples were collected before the first frost occurred.

Multiple discrete soil samples (500 cm3 each) for the analysis of nematode community structure were collected from each site up to a depth of 20 cm using a garden trowel, which was cleaned between sampling sites to prevent contamination. Samples were collected beneath plants (if present). For the analysis of selected soil properties (soil texture and total organic carbon), a sub-sample was collected from each site. The latter was performed only during the summer (February 2016) sampling interval as the selected soil properties were unlikely to change significantly during the sampling period. All the collected samples were labeled and transported in cool boxes to the Nematology Unit of the Agricultural Research Council – Plant Health and Protection (Roodeplaat, South Africa). Samples were stored at 10°C until further analysis.

Nematode extraction, preservation, and identification

Nematodes were extracted from 250 cm3 soil samples using an adapted decanting and sieving method, followed by a sugar flotation method (Marais et al., 2017). The nematodes were then fixed in a heated 4% formaldehyde and 1% propionic acid (FPG) solution (Netscher and Seinhorst, 1969), dehydrated in a glycerin solution (Seinhorst, 1962), and mounted in glycerin on permanent glass slides using a wax ring (Marais et al., 2017). Nematodes were counted using a De Grisse counting dish (De Grisse, 1963) and identified to genus level using an Olympus BX53F microscope. Taxonomic classification was based on Maggenti et al. (1988) and Geraert (2019) for Tylenchina; Kleynhans et al. (1996) for Heteroderidae; Hunt (1993) for Aphelenchida; Decraemer (1995) and Duarte et al. (2010) for Trichodoridae; Andrássy (2009) for Dorylaimidae; and Andrássy (2005, 2007) for other free-living nematodes. Feeding groups were assigned, based on Yeates et al. (1993), as herbivores (He), bacterivores (Ba), fungivores (Fu), omnivores (Om), or predators (Pr). Nematodes were also assigned a cp-value based on Bongers and Bongers (1998).

Physical and biological properties of soil samples

The selected physical and biological properties were analyzed by Eco-Analytica (North-West University, South Africa) as follows: total organic carbon (C) content of the soil samples were determined using the loss-on-ignition method (Donkin, 1991) and the soil texture as described by Laker and Dupreez (1982).

Statistical analysis

Rarefaction curves were used to compare taxa richness between habitats. This was achieved by calculating the sample-based Mao Tau estimator (with inter- and extrapolation) using Estimate S 9.1 Software Package. Furthermore, differences in the abundance of nematode trophic groups between habitats were investigated by pooling samples from across seasons. Abundance values were log10(x + 1) transformed and visualized, while statistical significance between trophic groups and habitats was inferred using two-way analysis of variance (ANOVA) and Tukey’s multiple comparisons tests. These analyses were performed using Graphpad Prism 6 Software Package.

Selected nematode-based indices (i.e. maturity index 2-5, enrichment index, structure index, channel index, basal index, and metabolic footprints) were used to quantify the condition, food web status, and functioning of soil ecosystems (Ferris et al., 2001; Ferris and Bongers, 2009; Ferris, 2010). These indices were calculated using the Nematode Indicator Joint Analyses (NINJA) online tool (Sieriebriennikov et al., 2014). The Graphpad Prism 6 Software Package was used to create a violin plot of the maturity index 2-5, a measure of soil ecosystem condition. Also, values of the structure and enrichment indices were used to plot a faunal analysis and characterize the food web status of the studied soil ecosystems (Ferris et al., 2001). Furthermore, linear mixed models with pairwise comparisons were used to determine if the independent variables, i.e. season (Winter, Spring, Summer, and Autumn) and habitat (OG, SRO, and RZ) significantly affected (singularly and interactively) the calculated nematode-based indices. Prior to this analysis the nematode-based indices were log transformed (where needed). All linear mixed models were performed using SPSS Statistics 25 Software Package.

Lastly, a redundancy analysis (RDA) was used to investigate the relationship between the response (nematode-based indices) and explanatory variables [(factors (seasons and habitats) and numeric variables (soil texture and total organic carbon)]. Predictor effects of the listed explanatory variables were calculated using a Monte-Carlo permutation test (Šmilauer and Lepš, 2014). These multivariate analyses were performed and illustrated on a biplot using Canoco 5 Software Package. Significance for all univariate and multivariate analyses was regarded at p < 0.05.

Results

Environmental conditions

The total rainfall for 90 days prior to each sampling event is illustrated in Figure 2A. This shows a general increase over time with the lowest and highest rainfall recorded prior to the Winter and Autumn sampling events, respectively. However, it is important to note that South Africa experienced a drought in 2015. According to the South African Weather Service (www.weathersa.co.za), the Telperion Nature Reserve region received only 388 mm from January to December 2015. By contrast, the recorded annual rainfall in 2014 and 2016 were 794 and 748 mm, respectively.

Figure 2:

The (A) total rainfall and (B) temperature (mean, minimum, and maximum) for the 90 days prior to each sampling event at Telperion Nature Reserve (South Africa).

Temperature (mean, minimum, and maximum) values are illustrated in Figure 2B with the lowest and highest mean temperatures recorded, as expected, prior to the Winter and Summer sampling events, respectively. The largest range between minimum (−2.6°C) and maximum (31.5°C) values were recorded before Winter sampling. Finally, the soil texture and organic carbon content of the collected and analyzed soil samples are provided as Supplementary material (Table S1).

Table S1.

Investigated soil properties measured at the open grassland, shrubland with rocky outcrops and riparian zone habitats in the Telperion Nature Reserve (Mpumalanga province, South Africa).

Note: Mean ± standard deviation values are reported.

Nematode community structure

Nematode communities associated with the studied habitats were investigated in terms of their taxonomic composition and richness, as well as their trophic abundance. In total, 93 genera, representing 50 families, were recorded in samples collected from Telperion Nature Reserve (Supplementary material Tables S2–S6) and were representative of the five major nematode trophic groups, namely, herbivores (Table S2), bacterivores (Table S3), fungivores (Table S4), omnivores (Table S5), and predators (Table S6). The largest number of taxa were found in the shrubland with rocky outcrops (SRO) habitat (77 genera), followed by the riparian zone (RZ) (74 genera) and open grasslands (OG) (68 genera) habitats. In terms of the number of recorded genera per trophic group, bacterivores dominated with 30 genera recorded in both the SRO and RZ habitats, while 24 bacterivore genera were found in the OG habitat. This was followed by herbivores (OG: 23; SRO: 23; RZ: 20), predators (OG: 11; SRO: 14; RZ: 13), fungivores (OG: 8; SRO: 9; RZ: 8), and omnivores (OG: 2; SRO: 2; RZ: 3). The results also showed that some nematode genera were present in all three habitats during all four seasons. This included Helicotylenchus, Rotylenchulus, Scutellonema, Tylenchus, Xiphinema (herbivores, Table S2); Acrobeles, Acrobeloides, Panagrolaimus, Prismatolaimus, Zeldia (bacterivores; Table S3); Aphelenchoides, Diphtherophora (fungivores; Table S4), and Dorylaimidae (omnivores, Table S5).

Table S2.

Occurrence and abundance of herbivores in 250 cm3 soil samples obtained from open grassland, shrubland with rocky outcrops and riparian zone sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

HabitatSand (%)Silt (%)Clay (%)Organic carbon (%)
Open grassland91.77  ±  3.033.75  ±  2.334.47  ±  0.940.78  ±  0.24
Shrubland with rocky outcrops86.45  ±  4.446.99  ±  3.136.54  ±  2.081.65  ±  0.84
Riparian zone86.48  ±  7.188.18  ±  5.165.33  ±  2.241.59  ±  1.26

Notes: Mean ± standard deviation values are reported. aGenus could not be identified.

Table S3.

Occurrence and abundance of bacterivores in 250 cm3 soil samples obtained from open grassland, shrubland with rocky outcrops and riparian zone sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

Open grasslandShrubland with rocky outcropsRiparian zone
FamilyGenus or FamilyWinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Anguinidae (He1) Subanguina 1.36 ± 4.5202.73 ± 9.050.45 ± 1.51000001.00 ± 2.2401.00 ± 2.24
Belondiridae (He5) Axonchium 0.91  ±  3.023.64 ± 4.5200.91 ± 3.0205.00 ± 12.693.00 ± 7.890.50 ± 1.58001.00 ± 2.240
Belonolaimidae (He3) Tylenchorhynchus 1.82 ± 6.030030.91 ± 65.9118.00 ± 56.920.50 ± 1.5825.00 ± 36.828.00 ± 13.980000
Criconematidae (He3)Criconematidaea 1.82 ± 4.625.00 ± 13.420.91 ± 3.020.91 ± 3.023.00 ± 3.503.00 ± 6.321.00 ± 3.164.50 ± 9.561.00 ± 2.2401.00 ± 2.240
Criconema 0.91 ± 3.0200000000002.00 ± 4.47
Criconemoides 6.36 ± 15.6710.45 ± 17.247.27 ± 11.918.18 ± 22.390.50 ± 1.589.00 ± 21.195.00 ± 14.1405.00 ± 11.185.00 ± 11.181.00 ± 2.240
Hemicriconemoides 01.82 ± 6.030001.00 ± 3.160001.00 ± 2.2400
Hemicycliophora 5.91 ± 19.602.27 ± 7.540.45 ± 1.511.82 ± 6.031.50 ± 4.740.50 ± 1.5801.00 ± 3.16165.00 ± 301.5813.00 ± 29.0769.00 ± 154.291.00 ± 2.24
Dorylaimellidae (He5) Dorylaimellus 58.18 ± 98.826.36 ± 13.808.64 ± 18.0411.36 ± 31.2321.50 ± 25.2812.50 ± 20.858.50 ± 16.5112.00 ± 25.1905.00 ± 11.181.00 ± 2.2455.00 ± 122.98
Heteroderidae (He3) Meloidogyne 4.09 ± 10.682.27 ± 7.541.82 ± 4.0525.45 ± 71.40154.50 ± 350.5513.00 ± 19.3211.50 ± 16.678.00 ± 17.030162.00 ± 342.9045.00 ± 81.7054.00 ± 120.98
Hoplolaimidae (He3)Hoplolaimidaea 6.36 ± 10.517.27 ± 16.9402.27 ± 4.1015.50 ± 26.2914.00 ± 26.0110.50 ± 31.4902.00 ± 4.4701.00 ± 2.240
Helicotylenchus 38.64 ± 43.25139.09 ± 246.4142.27 ± 112.3740.45 ± 64.6752.00 ± 82.2767.50 ± 102.2931.50 ± 57.4061.50 ± 12.836.00 ± 13.4283.00 ± 185.593.00 ± 6.7173.00 ± 139.45
Rotylenchulus 34.55 ± 71.5750.91 ± 128.7831.82 ± 47.2955.91 ± 125.8934.50 ± 92.8746.00 ± 143.7253.00 ± 139.9296.00 ± 275.772.00 ± 4.4728.00 ± 59.853.00 ± 6.711.00 ± 2.24
Rotylenchus 29.55 ± 42.5121.36 ± 38.61016.36 ± 22.5918.50 ± 35.045.50 ± 11.650.50 ± 1.587.00 ± 10.5915.00 ± 28.28008.00 ± 7.89
Scutellonema 110.45 ± 131.25172.73 ± 231.3838.64 ± 71.9135.45 ± 36.7744.00 ± 51.4764.00 ± 92.04131.00 ± 343.2421.00 ± 23.7819.00 ± 26.5540.00 ± 58.7435.00 ± 57.77137.00 ± 275.83
Leptonchidae (He4) Xiphinemella 6.82 ± 22.6101.36 ± 3.235.45 ± 9.344.50 ± 14.232.50 ± 7.910.50 ± 1.585.50 ± 12.57001.00 ± 2.242.00 ± 4.47
Longidoridae (He5) Longidorus 0000000003.00 ± 6.7100
Xiphinema 6.82 ± 13.652.27 ± 4.107.27 ± 10.815.91 ± 10.6820.22 ± 31.093.00 ± 5.3711.50 ± 16.1711.00 ± 16.963.00 ± 4.472.00 ± 2.746.00 ± 13.421.00 ± 2.24
Pratylenchidae (He3) Pratylenchus 10.45 ± 19.4222.73 ± 44.742.27 ± 3.445.91 ± 13.3827.00 ± 34.5016.00 ± 35.8115.50 ± 29.675.50 ± 10.660000
Psilenchidae (He2) Psilenchus 0.91 ± 2.02000003.00 ± 9.49031.00 ± 69.327.00 ± 15.6515.00 ± 33.540
Trichodoridae (He4) Nanidorus 02.73 ± 9.050002.50 ± 5.401.50 ± 4.740.50 ± 1.585.00 ± 7.0718.00 ± 28.42313.00 ± 683.1728.00 ± 49.32
Paratrichodorus 00002.00 ± 6.3200006.00 ± 13.4200
Trichodoridaea 00.45 ± 1.510000000000
Trichodorus 000000.50 ± 1.58000000
Tylenchidae (He2) Filenchus 00000000041.00 ± 80.9600
Coslenchus 02.73 ± 9.050.45 ± 1.51000.50 ± 1.58000000
Tylenchidaea 0.45 ± 1.511.36 ± 3.230002.50 ± 7.9100.50 ± 1.581.00 ± 2.2416.00 ± 30.5000
Tylenchus 77.73 ± 73.3333.64 ± 40.3815.91 ± 36.5321.82 ± 24.21121.50 ± 65.7518.00 ± 18.5914.00 ± 11.0126.00 ± 27.37139.00 ± 144.8085.00 ± 100.93246.00 ± 222.64189.00 ± 263.71
Tylencholaimidae (He4) Chitwoodius 0000.45 ± 1.510.50 ± 1.580000000
Tylenchulidae (He3) Meloidoderita 3.64 ± 12.0600000000000
Tylenchulus 000.45 ± 1.51002.00 ± 6.32000000
Total herbivores407.73 ± 194.43489.09 ± 476.05162.27 ± 214.97270.00 ± 206.95539.00 ± 383.70289.00 ± 240.71326.50 ± 363.48268.50 ± 289.86394.00 ± 353.77516.00 ± 496.61741.00 ± 922.38552.00 ± 734.72

Notes: Mean ± standard deviation values are reported. aGenus could not be identified.

Table S4.

Occurrence and abundance of fungivores in 250 cm3 soil samples obtained from open grassland, shrubland with rocky outcrops and riparian zone sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

Open grasslandShrubland with rocky outcropsRiparian zone
Family (functional guild)Genus or FamilyWinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Alaimidae (Ba4) Alaimus 98.64 ± 91.3317.73 ± 19.152.27 ± 11.0712.73 ± 19.0292.00 ± 57.0721.50 ± 22.7418.50 ± 22.6157.50 ± 63.5610.00 ± 14.1422.00 ± 49.191.00 ± 2.240
Amphidelus 000000.51 ± 1.585.50 ± 7.625.00 ± 11.550001.00 ± 2.24
Aphanolaimidae (Ba3) Aphanolaimus 000000001.00 ± 2.24000
Cephalobidae (Ba2)Cephalobidaea 013.18 ± 18.0787.73 ± 192.33113.64 ± 198.92019.00 ± 45.2632.00 ± 38.8236.50 ± 45.03003.00 ± 2.7410.00 ± 14.14
Acrobeles 122.27 ± 112.4869.09 ± 103.4636.82 ± 55.51160.45 ± 192.5751.50 ± 69.7652.50 ± 123.5218.00 ± 39.2424.50 ± 36.403.00 ± 6.7147.00 ± 64.5832.00 ± 54.1513.00 ± 16.43
Acrobeloides 60.45 ± 41.4420.45 ± 21.6225.00 ± 53.0623.64 ± 23.6724.00 ± 43.8366.50 ± 152.1943.00 ± 103.9821.50 ± 39.8625.00 ± 30.0052.00 ± 52.5162.00 ± 107.6263.00 ± 105.69
Cephalobus 3.18 ± 9.0242.27 ± 59.2254.09 ± 47.2197.73 ± 171.3718.00 ± 50.0164.50 ± 118.75116.00 ± 272.0163.50 ± 76.52075.00 ± 148.6658.00 ± 57.0726.00 ± 28.59
Cervidellus 08.64 ± 15.517.27 ± 15.557.73 ± 10.0908.00 ± 15.673.50 ± 5.304.50 ± 6.43050.00 ± 111.801.00 ± 2.245.00 ± 8.66
Chiloplacus 20.45 ± 22.524.09 ± 6.253.18 ± 10.550.91 ± 3.0213.00 ± 15.4913.50 ± 17.492.00 ± 6.3216.00 ± 44.0260.00 ± 128.654.00 ± 8.941.00 ± 2.240
Eucephalobus 00.91 ± 3.021.82 ± 4.6200121.00 ± 380.888.50 ± 16.345.50 ± 11.1709.00 ± 20.12167.00 ± 365.1040.00 ± 61.64
Heterocephalo-bellus 0000.45 ± 1.510000.50 ± 1.580000
Seleborca 53.64 ± 90.1157.73 ± 84.4289.55 ± 186.8357.73 ± 71.398.00 ± 22.0118.50 ± 46.197.50 ± 6.3514.50 ± 16.24017.00 ± 30.336.00 ± 10.841.00 ± 2.24
Zeldia 28.18 ± 32.9635.45 ± 77.052.73 ± 5.1813.18 ± 16.176.00 ± 7.388.00 ± 12.748.50 ± 13.754.00 ± 6.583.00 ± 6.7122.00 ± 49.192.00 ± 4.474.00 ± 8.94
Chromadoridae (Ba3) Punctodora 000001.00 ± 3.16009.00 ± 12.45000
Chronogasteridae (Ba3) Chronogaster 0000.45 ± 1.510.50 ± 1.5802.00 ± 4.830.50 ± 1.58279.00 ± 509.2012.00 ± 14.4095.00 ± 212.4328.00 ± 57.18
Cyatholaimidae (Ba4) Achromadora 004.55 ± 15.08000009.00 ± 15.172.00 ± 2.7418.00 ± 40.250
Diplogasteridae (Ba1) Butlerius 00001.50 ± 3.374.50 ± 14.23003.00 ± 6.712.00 ± 4.4700
Diplogasteroididae (Ba1)Diplogasteroididaea 000000000014.00 ± 31.301.00 ± 2.24
Elaphonematidae (Ba3) Elaphonema 8.64 ± 17.761.82 ± 4.628.64 ± 15.5112.73 ± 42.216.50 ± 11.5607.00 ± 22.141.50 ± 3.370000
Mesorhabditidae (Ba1) Mesorhabditis 03.18 ± 10.5508.18 ± 27.14012.50 ± 32.681.50 ± 3.372.00 ± 6.32014.00 ± 21.912.00 ± 4.470
Metateratocepha-lidae (Ba3) Euteratocephalus 000000005.00 ± 8.66000
Monhysteridae (Ba2) Monhystera 1.82 ± 6.030002.50 ± 4.862.00 ± 6.32008.00 ± 13.044.00 ± 6.523.00 ± 4.476.00 ± 10.84
Monhystrella 000000000002.00 ± 4.47
Neodiplogastridae (Ba1) Fictor 00000000.50 ± 1.582.00 ± 4.47000
Koerneria 00002.00 ± 6.3200001.00 ± 2.243.00 ± 6.710
Odontolaimidae (Ba3) Odontolaimus 000000002.00 ± 4.47000
Osstellidae (Ba2) Drilocephalobus 0018.18 ± 50.7137.73 ± 71.46003.50 ± 4.7412.00 ± 14.570001.00 ± 2.24
Panagrolaimidae (Ba1) Panagrolaimus 26.36 ± 25.3115.00 ± 23.020.91 ± 3.0217.27 ± 35.3117.50 ± 14.1917.50 ± 39.3212.50 ± 23.834.00 ± 7.3839.00 ± 39.1215.00 ± 21.2119.00 ± 21.336.00 ± 6.52
Plectidae (Ba2) Anaplectus 00000000.50 ± 1.580000
Ceratoplectus 0000.91 ± 2.0210.50 ± 33.204.50 ± 11.172.50 ± 7.9111.00 ± 25.0302.00 ± 4.4700
Plectus 5.45 ± 12.340004.00 ± 12.6501.50 ± 2.423.50 ± 11.074.00 ± 4.181.00 ± 2.242.00 ± 4.471.00 ± 2.24
Wilsonema 3.18 ± 7.830.45 ± 1.5101.36 ± 2.346.50 ± 10.551.50 ± 2.422.50 ± 5.406.50 ± 17.1706.00 ± 13.421.00 ± 2.241.00 ± 2.24
Prismatolaimidae (Ba3) Prismatolaimus 0.91 ± 2.020.91 ± 2.020.91 ± 3.023.18 ± 10.556.50 ± 11.807.00 ± 11.357.50 ± 13.997.00 ± 12.523.00 ± 2.745.00 ± 11.186.00 ± 13.423.00 ± 4.47
Rhabditidae (Ba1) Cruznema 00000000.50 ± 1.580000
Rhabditidaea 0.45 ± 1.5100000000000
Rhabditis 1.36 ± 3.232.27 ± 4.6703.18 ± 5.1318.00 ± 31.5512.50 ± 34.422.00 ± 4.225.00 ± 11.557.00 ± 13.046.00 ± 8.9400
Rhabdolaimidae (Ba3) Rhabdolaimus 000.45 ± 1.510000.50 ± 1.5800000
Teratocephalidae (Ba3) Teratocephalus 0.45 ± 1.510.45 ± 1.510004.00 ± 11.01004.00 ± 6.52000
Total bacterivores435.45 ± 275.30293.64 ± 291.33349.09 ± 507.99573.18 ± 517.51288.50 ± 147.82460.50 ± 719.66306.00 ± 435.38308.00 ± 126.34476.00 ± 506.13368.00 ± 364.02496.00 ± 603.15212.00 ± 129.79

Notes: Mean ± standard deviation values are reported. aGenus could not be identified.

Table S5.

Occurrence and abundance of omnivores in 250 cm3 soil samples obtained from open grassland, shrubland with rocky outcrops and riparian zone sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

Open grasslandShrubland with rocky outcropsRiparian zone
FamilyGenus or FamilyWinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Anguinidae (Fu2) Ditylenchus 8.64 ± 13.9829.55 ± 40.2828.18 ± 33.8610.91 ± 13.751.50 ± 4.7420.50 ± 38.1819.00 ± 19.9723.50 ± 30.56025.00 ± 34.64130.00 ± 144.27110.00 ± 186.38
Aphelenchidae (Fu2) Aphelenchus 15.91 ± 16.71133.18 ± 190.5686.36 ± 246.8024.55 ± 35.6712.50 ± 13.5951.50 ± 63.2080.00 ± 161.579.00 ± 9.66019.00 ± 26.0872.00 ± 102.0268.00 ± 141.01
Aphelenchoididae (Fu2) Aphelenchoides 8.64 ± 15.3410.45 ± 15.7221.36 ± 61.085.91 ± 7.0113.00 ± 23.12284.50 ± 887.4016.00 ± 30.264.00 ± 6.5815.00 ± 23.9813.00 ± 13.0421.00 ± 20.7412.00 ± 21.39
Diphtherophoridae (Fu3) Diphtherophora 23.18 ± 22.288.64 ± 13.625.00 ± 7.420.91 ± 2.0224.00 ± 19.553.00 ± 6.3223.50 ± 67.3732.50 ± 86.359.00 ± 15.176.00 ± 13.4234.00 ± 73.2632.00 ± 35.81
Leptonchidae (Fu4) Leptonchus 01.36 ± 3.230.91 ± 3.020.91 ± 3.0201.00 ± 2.1100.50 ± 1.5801.00 ± 2.2400
Proleptonchus 003.18 ± 10.5500000.50  ± 1.580000
Neotylenchidae (Fu2)Neotylenchidaea 5.91 ± 10.203.64 ± 5.058.64 ± 15.0210.00 ± 13.045.00 ± 15.8111.50 ± 20.1512.50 ± 18.4518.50 ± 19.87017.00 ± 28.2023.00 ± 51.4310.00 ± 17.32
Paraphelenchidae (Fu2) Paraphelenchus 00002.00 ± 6.320001.00 ± 2.2409.00 ± 20.120
Tylencholaimellidae (Fu4) Tylencholaimellus 6.36 ± 11.8500.45 ± 1.510.45 ± 1.153.50 ± 8.180.50 ± 1.582.00 ± 6.3200001.00 ± 2.24
Tylencholaimidae (Fu4)Tylencholaimidaea 2.73  ± 9.0500000001.00 ± 2.24000
Tylencholaimus 7.27 ± 24.121.82 ± 4.0524.09 ± 65.831.36 ± 4.527.00 ± 14.7609.00 ± 28.4602.00 ± 4.4708.00 ± 15.250
Total fungivores78.64 ± 54.13188.64 ± 229.54178.18 ± 342.9955.00 ± 41.2368.50 ± 35.83372.50 ± 863.64162.00 ± 179.2088.50 ± 120.2128.00 ± 22.5381.00 ± 93.10297.00 ± 341.31233.00 ± 360.60

Notes: Mean ± standard deviation values are reported. aGenus not identified.

Table S6.

Occurrence and abundance of predators in 250 cm3 soil samples obtained from open grassland, shrubland with rocky outcrops and riparian zone sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

Open grasslandShrubland with rocky outcropsRiparian zone
FamilyGenus or FamilyWinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Dorylaimidae (Om4)Dorylaimidaea 55.91 ± 22.8920.45 ± 23.0711.36 ± 19.2517.27 ± 21.1449.00 ± 29.2314.00 ± 14.1017.50 ± 29.6516.00 ± 16.9615.00 ± 12.7515.00 ± 14.146.00 ± 6.5218.00 ± 40.25
Dorylaimus 0.45 ± 1.510.45 ± 1.51000.50 ± 1.5801.00 ± 3.1603.00 ± 6.71000
Mesodorylaimus 00003.00 ± 9.491.50 ± 4.7402.50 ± 7.911.00 ± 2.243.00 ± 4.471.00 ± 2.242.00 ± 4.47
Tylencholaimidae (Om4) Discomyctus 000000000006.00 ± 13.42
Mydonomidae (Om4) Dorylaimoides 0.91 ± 3.0200000000000
Total omnivores57.27 ± 23.6020.91 ± 23.0011.36 ± 19.2517.27 ± 21.1452.50 ± 31.3815.50 ± 14.99018.50 ± 29.16018.50 ± 20.01019.00 ± 9.6218.00 ± 16.817.00 ± 5.7026.00 ± 52.73

Notes: Mean ± standard deviation values are reported. aGenus could not be identified.

A comparison of taxa richness between habitats was made using rarefaction curves (Fig. 3). These curves illustrate the sample-based observed (solid line) and extrapolated (dotted line) Mau Tau richness, while 95% confidence intervals are indicated as shaded bands. When comparing richness at the 20 samples mark (total number of samples collected at the RZ habitat), it is clear that RZ presented the greatest richness with 84 genera. This was followed by SRO and OG, which presented 73 and 65 genera, respectively. Nonetheless, it is important to note that these curves did not reach an asymptote (leveling), suggesting that additional sampling might be required to record the complete nematode community. To this end, extrapolation revealed that more than 50 samples per habitat are likely required at which point an estimated 98, 91, and 83 genera would be predicted to be collected at the RZ, SRO, and OG habitats, respectively.

Figure 3:

Rarefaction curves illustrate the sample-based observed (solid line) and extrapolated (dotted line) Mau Tau richness of nematode genera at the open grassland (OG), shrubland with rocky outcrops (SRO), and riparian zone (RZ) habitats in Telperion Nature Reserve (South Africa). The shaded bands represent the 95% confidence intervals.

Lastly, the nematode trophic abundances per habitat were considered. Figure 4 shows that herbivores and bacterivores were the most abundant trophic groups, followed by fungivores, omnivores, and predators. The multiple comparisons tests revealed that although no significant differences occurred between herbivores and bacterivores, both these trophic groups differed significantly from fungivores, omnivores, and predators at all three habitats. Also, significant differences between fungivore, omnivore, and predator abundances were recorded at all three habitats with the exception of no significant difference between fungivores and omnivores/predators at the RZ habitat. No significant differences in nematode trophic abundances were recorded between habitats.

Figure 4:

Log10(x + 1)-transformed nematode abundances (mean, minimum, and maximum) per trophic group at the open grassland (OG), shrubland with rocky outcrops (SRO), and riparian zone (RZ) habitats in Telperion Nature Reserve (South Africa). Bars with common superscript do not differ significantly (p > 0.05).

Soil ecosystem condition and food web status

The maturity index 2-5 was used to investigate the soil ecosystem condition of the studied grassland habitats. A violin plot (Fig. 5) illustrates the median, quartiles, minimum, maximum, and data distribution (shape outline) of this index per habitat per season. Minimum and maximum values ranged from 2 (RZ in Spring) to 3.7 (SRO in Winter), respectively, while the highest and lowest median values for all three habitats were recorded during the Winter and Spring, respectively. Also worth noting is that the SRO habitat had the highest median values during all four seasons. This habitat was thus classified as having the most stable soil ecosystem. The results from the linear mixed models (Table 1) showed that habitat and season had a significant effect on the maturity index 2-5. Yet, pairwise comparisons revealed that these effects were only evident between seasons, which showed significantly higher maturity index 2-5 values in Winter compared to the rest of the seasons. Differences between Spring, Summer, and Autumn were not significant.

Table 1

Mixed linear models reveal the effect (singularly and interactively) of the independent variables, season and habitat, on the calculated nematode-specific indices.

Open grasslandShrubland with rocky outcropsRiparian zone
FamilyGenus or FamilyWinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Actinolaimidae (Pr5) Neoactinolaimus 00000000002.00 ± 4.470
Anatonchidae (Pr4) Iotonchus 4.09 ± 7.010.45 ± 1.510.45 ± 1.5109.50 ± 12.791.50 ± 2.423.00 ± 4.221.00 ± 3.162.00 ± 4.47014.00 ± 31.304.00 ± 8.94
Aphelenchoididae (Pr2) Seinura 0000000000044.00 ± 98.39
Aporcelaimidae (Pr5)Aporcelaimidaea 1.36 ± 4.52001.82 ± 4.6201.00 ± 3.1600.50 ± 1.580003.00 ± 6.71
Aporcelaimellus 16.36 ± 16.750.91 ± 3.025.91 ± 9.951.36 ± 4.523.50 ± 11.073.50 ± 4.129.50 ± 13.834.50 ± 9.562.00 ± 4.47000
Aporcelaimus 1.36 ± 3.235.00 ± 9.495.00 ± 9.224.09 ± 12.0003.50 ± 5.803.50 ± 11.0701.00 ± 2.242.00 ± 2.746.00 ± 13.4219.00 ± 39.75
Makatinus 00.91 ± 3.0200001.50 ± 4.740.50 ± 1.580000
Ironidae (Pr4) Ironus 000000004.00 ± 8.94000
Mononchidae (Pr4)Mononchidaea 0.45 ± 1.510.45 ± 1.5100001.00 ± 3.16002.00 ± 2.7400
Cobbonchus 000000003.00 ± 4.47000
Granonchulus 000001.50 ± 4.74000000
Mononchus 0.45 ± 1.510.45 ± 1.5101.82 ± 6.0314.00 ± 29.040.50 ± 1.58002.00 ± 2.74000
Mylonchulidae (Pr4) Mylonchulus 000.91 ± 3.0202.00 ± 4.830000000
Paraxonchiidae (Pr5) Paraxhonchium 4.55 ± 10.601.36 ± 3.233.64 ± 12.0601.50 ± 4.7403.00 ± 7.896.00 ± 14.490000
Qudsianematidae (Pr4)Discolaiminaea 1.36 ± 4.520002.00 ± 4.83004.00 ± 8.760000
Discolaimoides 9.09 ± 14.802.73 ± 5.1812.73 ± 25.4320.91 ± 27.281.50 ± 3.3701.00 ± 3.169.00 ± 12.872.00 ± 2.74003.00 ± 6.71
Discolaimus 0.91 ± 3.020.45 ± 1.5101.36 ± 4.520.50 ± 1.582.50 ± 5.400.50 ± 1.580.50 ± 1.5801.00 ± 2.241.00 ± 2.242.00 ± 2.74
Eudorylaimus 000000.50 ± 1.5800002.00 ± 4.470
Labronema 01.82 ± 6.03000.50 ± 1.580.50 ± 1.581.00 ± 3.160.50 ± 1.58006.00 ± 8.940
Tobrilidae (Pr4) Tobrilus 00001.00 ± 3.160002.00 ± 2.741.00 ± 2.241.00 ± 2.241.00 ± 2.24
Tripylidae (Pr3) Tripyla 0000.45 ± 1.511.50 ± 4.740.50 ± 1.581.00 ± 3.1600000
Total predators40.00 ± 28.8114.55 ± 20.3028.64 ± 24.4031.82 ± 27.0437.5 ± 34.5615.5 ± 11.2024.5 ± 21.6726.5 ± 24.2318.00 ± 13.516.00 ± 4.1832.00 ± 39.1576.00 ± 131.79
Figure 5:

Maturity index 2-5 violin plot of values recorded during the studied seasons at the open grassland (OG), shrubland with rocky outcrops (SRO), and riparian zone (RZ) habitats in Telperion Nature Reserve (South Africa). The bold line within each plot indicates the median value, dashed lines are quartiles, minimum, and maximum are the bounds of the plot and data distribution is shown by plot shape outline.

The status of the soil food webs were investigated using the nematode faunal analysis (Fig. 6). This analysis showed that all three habitats had the highest structure (complexity) in the Winter with the SRO and RZ habitats classified as maturing and enriched with balanced (bacterial vs. fungal) decomposition pathways. The OG habitat was classified as mature and fertile characterized by fungal dominated decomposition. In Spring, however, nematode community composition in all three habitats represented degraded and depleted food webs. During Summer and Autumn, the communities represented either degraded and depleted or mature and fertile food webs. The faunal analysis also revealed the clustering of habitats (i.e. low variability between habitats) especially during Winter and Spring. The effect of season and habitat on the two indices used to determine the food web status, namely the structure index and enrichment index (Fig. 6), were further investigated using linear mixed models. This revealed that while both season and habitat had a significant effect on the enrichment index, only season had a significant effect on the structure index. Pairwise comparisons showed, as with the maturity index 2-5, that communities in Winter had significantly higher enrichment and structure compared to the other seasons, while the OG habitat had significantly lower enrichment compared to the SRO and RZ habitats.

Figure 6:

Nematode faunal analysis indicating the food web status during the studied seasons at the open grassland (green), shrubland with rocky outcrops (orange), and riparian zone (blue) habitats in Telperion Nature Reserve (South Africa). Winter and Spring seasonal clusters are circled.

Finally, a redundancy analysis illustrated on a biplot (Fig. 7) was used to investigate the relationship between the explanatory and response variables. The explanatory variables accounted for 42% (p < 0.001) of the observed variation in the response variables with axes 1 and 2 representing 35.3% (p < 0.001) and 5.5% (p < 0.001), respectively. This analysis revealed that the three habitats could be differentiated based on the predominant soil fraction, i.e. the OG, RZ, and SRO habitats consisted of more sandy, silty, or clayey soils, respectively. Nonetheless, only clay presented a significant effect of 4.3% on the response variables. In addition, organic carbon levels were positively correlated with the RZ habitat and enrichment index, but had only a minor effect (< 4%) on the response variables. Furthermore, the redundancy analysis biplot showed that a positive correlation existed between the OG habitat and the channel and basal indices. The SRO habitat, in turn, presented a positive correlation to the maturity 2-5 and structure indices. Finally, the Monte-Carlo permutation test revealed that 30% of the variation in response variables can be attributed to seasonality.

Figure 7:

Redundancy analysis of the relationship between the response (nematode-based indices) and explanatory variables (factors: seasons and habitats; numeric variables: soil texture and total organic carbon). Measurements were taken during four consecutive seasons at open grassland (OG), shrubland with rocky outcrops (SRO), and riparian zone (RZ) habitats in Telperion Nature Reserve (South Africa).

Soil ecosystem functioning

Metabolic footprints (Table 2) were used to assess the soil ecosystem functioning, i.e. the magnitude of functions and services delivered by the soil ecosystems. Generally, relatively high herbivore footprints were recorded, while especially high herbivore footprints were observed in Winter and Spring at the SRO and RZ habitats, respectively. However, even though clear differences were recorded in the herbivore footprint between seasons and habitats, the linear mixed models (Table 1) revealed that none of these differences were significant. Of the remaining footprints, significant seasonal effects were only recorded for the predator and structure footprints. Pairwise comparisons revealed that the predator footprints were significantly higher in the Winter compared to Spring and Summer, while the structure footprints were significantly higher in the Winter compared to Spring and Autumn. No significant habitat effects were recorded for any of the footprints.

Table 2.

Metabolic footprints calculated from 250 cm3 soil samples obtained from riparian zone, shrubland with rocky outcrops and open grassland sites in the Telperion Nature Reserve (Mpumalanga province, South Africa) as part of an ecological study that was conducted over four seasons from 2015 to 2016.

VariableSource F ratio p valueVariableSource F ratio p value
Maturity index 2–5Habitat3.66 0.04 Bacterivore footprintHabitat0.040.96
Season14.73 <0.001 Season2.460.07
Habitat × Season0.420.86Habitat × Season0.580.75
Channel indexHabitat2.490.11Fungivore footprintHabitat0.350.71
Season11.67 <0.001 Season1.520.22
Habitat × Season1.380.24Habitat × Season1.010.42
Basal indexHabitat4.96 0.01 Omnivore footprintHabitat0.760.48
Season11.08 <0.001 Season2.540.06
Habitat × Season0.350.91Habitat × Season0.290.94
Enrichment indexHabitat8.77 <0.001 Predator footprintHabitat1.450.26
Season7.65 <0.001 Season4.53 <0.01
Habitat × Season1.080.38Habitat × Season0.710.65
Structure indexHabitat2.000.16Enrichment footprintHabitat1.490.25
Season13.33 <0.001 Season0.290.83
Habitat × Season0.640.70Habitat × Season0.840.54
Herbivore footprintHabitat0.560.58Structure footprintHabitat0.110.90
Season1.600.20Season4.73 <0.01
Habitat × Season1.340.25Habitat × Season0.520.79

Note: Mean ± standard deviation values are reported.

Discussion

The nematode community structure

Grasslands typically support a high richness of nematode taxa (Song et al., 2017). Yeates and Lee (1997) and Bell et al. (2005) reported 32 and 70 genera, respectively, from tussock grasslands in New Zeeland. Popovici and Ciobanu (2000), in turn, listed between 33 and 67 genera from different grassland regions in Romania, while Čerevková (2006) reported 64 genera in grasslands from the Slovakia Republic. The present study, however, reported a substantially higher number of nematode genera of 93 in total. Rarefaction curves suggest that increased sampling efforts at these sites would result in further taxa being added to that total. At the RZ habitat, for example, extrapolation of the Mau Tao richness indicator suggests that close to a 100 nematode genera could be recovered with increased sampling efforts. The high number of recorded and estimated genera is important as it is commonly accepted that a high richness in nematode taxa are indicative of healthy and functioning soil ecosystems (Ekschmitt et al., 2001; Ferris, 2010; Sánchez-Moreno and Ferris, 2018). A high richness of herbivore taxa in natural environments are often the result of a high diversity in plant species (Hodda et al., 2009; Cortois et al., 2017). According to Hodda et al. (2009), some evidence suggests that a potential reason for this is the lack of competition resulting from root systems creating very defined niches.

Another important consideration when looking at nematode community structure is the abundance distribution between trophic groups (Ferris et al., 2001; Ferris, 2010; Van Den Hoogen et al., 2019). A recent study by Van Den Hoogen et al. (2019) investigated the abundance of nematode trophic groups across all major terrestrial biomes and showed that, on average, bacterivores dominated, followed by herbivores, fungivores, omnivores, and predators. The same trend was observed for temperate grasslands, which is the classification under which the South African grassland biome is listed (Van Den Hoogen et al., 2019). However, even though bacterivores are reportedly the most abundant trophic group on a global and biome scale, studies have reported regional grassland nematode communities with large and even dominating herbivore populations (Porazinska et al., 2003; Ayres et al., 2008). Similarly, during the present study, no significant difference was recorded between the two most abundant groups, namely herbivores and bacterivores, at any of the habitats.

It is worth noting that herbivores of economic importance (e.g. Meloidogyne) occurred in relatively high numbers at some sampling sites. Other herbivores (e.g. Helicotylenchus, Rotylenchulus, and Scutellonema) were present in all the habitats during all four seasons. The presence of these herbivore taxa, sometimes in high numbers, poses a threat to crop production in this region. According to various studies conducted in production areas where major grain crops (e.g. maize and soybean) are produced, Meloidogyne is generally the most predominant nematode pest. Substantial yield losses in maize (up to 60%) (Riekert and Henshaw, 1998) and soybean (up to 100%) (Smit and De Beer, 1998; Fourie et al., 2013) due to infection by this genus demonstrate its adverse impact on crop production. Moreover, Helicotylenchus, Rotylenchulus, and Scutellonema are also abundant in local crop fields, indicating their potential impact on crop production (Fourie et al., 2001; Bekker et al., 2016; Mbatyoti et al., 2018). Therefore, grasslands may act as reservoirs for herbivores that pose a threat to crop production should such areas be converted into agricultural fields.

Seasonal and habitat induced changes in grassland ecosystems

Seasonality had a large and significant effect on the condition, food web status and functioning of the studied grassland soil ecosystems. This effect was most evident progressing from Winter to Spring. During Winter the habitats were enriched with stable ecosystems and structured food webs, while depletion with reduced stability and structure were evident in Spring. Furthermore, during the Winter increased ecosystem functioning were evident with significantly higher predator and structure footprints. This even though low temperatures and reduced precipitation (as recorded during the 90 days period prior to Winter sampling) typically exert a negative influence on nematode communities (Verschoor et al., 2001; Song et al., 2017; Andriuzzi et al., 2020). According to Song et al. (2017), the optimal temperature for the survival and reproduction of soil nematodes ranges from 20 to 25°C, while extremes below 5 and above 30°C significantly inhibit development. Similarly, water availability influences primary production and thus energy flow into the soil food web, while also regulating organic carbon decomposition (Andriuzzi et al., 2020). It should be noted that the especially high herbivore footprints can be attributed to the large number of Meloidogyne that were present in some of the samples. Metabolic footprints are calculated using the biomass of adult females (Ferris, 2010), which can substantially increase footprint values if large numbers of, for example, some sedentary endoparasitic nematodes are recorded.

In contrast to seasonal effects, habitats and the measured soil properties had minimal influence on the studied soil ecosystems. This despite habitat type, soil texture and organic carbon levels being known to substantially influence soil faunal communities (Du Preez et al., 2018; Sánchez-Moreno and Ferris, 2018; Sprunger et al., 2019). Nonetheless, the positive correlation (as evidenced by the redundancy analysis) between organic carbon and the enrichment index is in line with the general understanding that higher organic carbon levels support opportunistic bacterial-feeding nematodes (Sánchez-Moreno and Ferris, 2018). Increased levels of soil organic carbon at the RZ habitat is therefore likely the reason for this habitat also having the highest nematode community enrichment levels during all four seasons in the soil food web analysis. Conversely, the basal and channel indices were positively correlated with the OG habitat, indicating that soil ecosystems associated with this habitat were likely more degraded, but with greater fungal decomposition (Ferris et al., 2001; Ferris, 2010; Sánchez-Moreno and Ferris, 2018).

Interestingly, the abundance of nematode taxa (Supplementary material: Tables S2–S6) made it clear that the dramatic shift in nematode community in Winter can largely be attributed to fluctuations in the populations of only a few predacious nematodes (namely Aporcelaimellus, Iotonchus, and Mononchus), as well as one omnivorous nematode group, the Dorylaimidae. These taxa are assigned a colonizer-persister (c-p) value of 4 (Dorylaimidae and Iotonchus) and 5 (Aporcelaimellus and Mononchus) (Sieriebriennikov et al., 2014), which means that they are regarded as persisters (or k-strategists). Their presence in large numbers therefore infers greater ecosystem stability and food web structure (Ferris et al., 2001; Sánchez-Moreno and Ferris, 2018), as was recorded during Winter at the SRO and OG habitats. The reason for the increased number of members of these nematode groups recorded during Winter is not explicitly clear, however, findings from previous studies provide some potential answers. According to De Ley (1992), desert and dune sands often contain a remarkably high proportion of dorylaims, while Arpin (1969) repeatedly found dorylaims (predators and omnivores) in desiccation experiments. De Ley and Mundo-Ocampo (2004), in turn, stated that the ecological sensitivity of these nematodes may have been overestimated and that certain species may be capable of surviving desiccation and freezing. Omnivorous nematodes typically also have versatile feeding habits and can probably interact at various levels of the soil food web (Hanel, 2003). It is therefore possible that dorylaims in the studied Rand Highveld Grasslands of South Africa are not adversely affected by the conditions present in Winter and may even increase in numbers following reduced competition for resources from other faunal groups. The possibility that certain nematode families and genera are favored by Winter conditions (low temperature and rainfall), infers that nematode community responses to altered conditions are taxa dependent (see also Papatheodorou et al., 2004).

Conclusion

The grassland habitats of Telperion Nature Reserve host a diversity of nematode taxa, many of which are likely still to be observed or described. Although the studied nematode communities were dominated by herbivores and bacterivores, it was the less abundant nematode groups, i.e. omnivores and predators, which had a significant influence on the condition, food web status, and functioning of the soil ecosystems. Furthermore, evidence suggests that only a few taxa are largely responsible for seasonal shifts in measurable soil ecosystem parameters.

This study also recorded a high abundance in nematodes of economic importance, which has important implications for crop production in this grassland biome.

Acknowledgements

The authors thank Oppenheimer Generations, South Africa for access to their property, Telperion Nature Reserve, and the opportunity to collect samples. The following people are also thanked: Duncan MacFadyen, manager of research and conservation at Oppenheimer Generations, Elsabe Bosch, manager at Telperion Nature Reserve, Adoration Shubane and Elsa van Niekerk at ARC-PHP for technical assistance, Christel du Preez for constructive comments, Pieter Holtzhausen for help with creating maps, the biometrician Suria Ellis of the NWU for her assistance with the analyses of the data and Robert Colwell for help with rarefaction curves. We thank the Professional Development Program of the Agricultural Research Council for funding this project.

Open grassland (OG)Shrubland with rocky outcrops (SRO)Riparian zone (RZ)
WinterSpringSummerAutumnWinterSpringSummerAutumnWinterSpringSummerAutumn
Herbivore footprint120 ± 70118 ± 11953 ± 63168 ± 265682 ± 1,326116 ± 111124 ± 117116.11 ± 13497 ± 81683 ± 1,300331 ± 326289 ± 434
Bacterivore footprint84 ± 4956 ± 6250 ± 7499 ± 8785 ± 6391 ± 15647 ± 5257.21 ± 3089 ± 6266 ± 5865 ± 6729 ± 12
Fungivore footprint14 ± 1123 ± 2625 ± 399 ± 711 ± 634 ± 6023 ± 2317.56 ± 224 ± 416 ± 1849 ± 5837 ± 55
Omnivore footprint179 ± 97103 ± 11089 ± 9491 ± 155123 ± 6778 ± 69102 ± 12458 ± 5063 ± 3454 ± 5389 ± 150241 ± 477
Predator footprint12 ± 113 ± 46 ± 1011 ± 1131 ± 267 ± 97 ± 88 ± 815 ± 144 ± 322 ± 4215 ± 13
Enrichment footprint16 ± 830 ± 4019 ± 3821 ± 1846 ± 6568 ± 14724 ± 2123.46 ± 2528 ± 2635 ± 3451 ± 4433 ± 55
Structure footprint217 ± 109114 ± 112103 ± 94107 ± 153178 ± 6291 ± 78122 ± 14086 ± 62123 ± 4766 ± 61132 ± 140266 ± 481
DOI: https://doi.org/10.21307/jofnem-2020-118 | Journal eISSN: 2640-396X | Journal ISSN: 0022-300X
Journal RSS Feed
Language: English
Page range: 1 - 22
Published online: Dec 14, 2020
Published in print: Jan 1, 2020
Published by: Society of Nematologists, Inc.
In partnership with: Paradigm Publishing Services
Publication frequency: 1 issue per year
Keywords:
Ecosystem functioning,
Food web status,
Grassland habitats,
Seasonal variation,
Soil ecology
Related subjects:
Life sciences,
Life sciences, other

© 2020 Chantelle Girgan, Gerhard du Preez, Mariette Marais, Antoinette Swart, Hendrika Fourie, published by Society of Nematologists, Inc.
This work is licensed under the Creative Commons Attribution 4.0 License.

Volume 52 (2020): Issue 1 (January 2020)
Previous article
Next article