Felicity V. Crotty
Sina M. Adl- Department of Soil Science, University of Saskatchewan, Saskatoon, SK, Canada
The soil food web is often described as having three main energy channels: root, bacterial and fungal. Here we provide quantitative data using a sensitive stable isotope ratio mass spectrometry procedure with microcosms on species interactions in the fungal pathway. We measured 15N and 13C enrichment in microarthropods through grazing rare isotope enriched fungal mycelia. Experimental treatments were various combinations of 1, 2, 3, 4 microarthropods species. We used three fungivores (the collembolan Lepidocyrtus curvicollis, the Astigmata Tyrophagus putrescentiae, the Oribatida Oribatula tibialis), and the Mesostigmata predator Hypoaspis acquilifer. We collected individuals of each species separately, as well as their feces, and molt where available. All three fungivorous microarthropods consumed significantly more than their own body weight per day. The three fungivores differed in their consumption of the mycelium as it was not equally palatable to each. The Mesostigmata predator Hypoaspis also differed in its microarthropod prey preference. In multiple species combinations microarthropod behavioral interactions modified consumption and predation rates. Our selection of mites of different sizes, with varied preference for the mycelium, combined with differing predation rates on each mite, demonstrate that even three trophic level interactions with only five interacting species are not predictably simple. The interpretation of the stable isotope results and consumed-excreted weights indicate that: (a) behavior and microscopic observations should not be ignored in competition-predation interactions, and (b) functional guilds can take advantage of more diverse food opportunities. The reality of mixed diets complicates functional guild assignments that are reflected in 15N and 13C isotope levels at natural abundances in the environment. Microcosm experiments with this sensitive technique can help decipher the interpretation of rare isotope natural abundance values, as well as providing measured consumption, growth, and excretion rate values for modeling soil food web interactions.
Introduction
Soil food web models typically recognize three main routes energy can flow through the soil – roots, the fungal or bacterial pathway (Moore et al., 2005). An early study of soil food web structure found that both top-down control and bottom-up control existed simultaneously and were necessary for ecosystem stability (de Ruiter et al., 1995). Increased predation rate on prey, or of grazing consumption (e.g., on bacteria lawn, a mycelium, or pasture) decreases the biomass and abundance of the species consumed, while increasing the abundance of the consumers over time. In food webs, the increased abundance or biomass of the consumer trophic level, relative to the consumed trophic level is called a trophic cascade. It makes allusion to the appearance in graphs of successive reduced and increased biomass in the trophic levels in the system. These increasing and decreasing consumption rates oscillate through time and ripple through generations and are one cause of population dynamics. An interesting and useful parameter for ecologists is to find out at what relative consumption rates trophic cascades occur, for they cannot occur under elevated food resources (Polis et al., 2000). It has been postulated that the microarthropods, particularly the Collembola and Oribatida, have a more significant role in the fungal pathway (Faber, 1991; Schneider and Maraun, 2005; Jonas et al., 2007). Several additional papers discussed fungivory by soil microarthropods (e.g., Maraun et al., 2003; Schneider et al., 2004; Pollierer et al., 2009; Thiele-Bruhn et al., 2012), but few tried to quantify rates of hyphae consumption except for individual species, such as Collembola (Jonas et al., 2007), while others focused on feeding preferences (e.g., Ruess et al., 2005; Koukol et al., 2009; Semenina and Tiunov, 2011). However, little is known regarding the quantities of fungal mycelium consumed and assimilated by microarthropods species when in monoculture, in competition, or with predation pressure. Likewise there is little insight into the microarthropods foraging behavior within the soil itself, or their response to the distribution and quality of food resources (Adl, 2003; Hassall et al., 2006). Considering there are 101–102 m of hyphae in one gram of fertile soil (Leake et al., 2003), this is a large biomass contribution to soil nutrient turnover that ought to be studied more. The difficulty has been in designing microcosms with compatible species, and to develop a technique with sufficient sensitivity to work with the small number of individuals in a microcosm experiment.
One of the main methods of differentiating between feeding preferences and measuring consumption rates is the use of stable isotopes (Tiunov, 2007; Semenina and Tiunov, 2011), as this can provide a time-integrated measure of the trophic position of soil animals (Pollierer et al., 2009). An organisms’ tissues exhibit a fixed isotopic enrichment in relation to their diet; therefore stable isotopes can be measured to assess the assimilation rates over the long-term (Peterson and Fry, 1987). Introducing an enriched food source as a pulse provides a mechanism to accurately monitor consumption in a controlled environment (Bradford et al., 2012). The addition of a substrate with a distinct isotopic signature can be traced into newly synthesized compounds and tissues (Elfstrand et al., 2008). The introduction of a food source which is enriched above natural abundance levels with stable isotopes into the soil food web can elucidate feeding interactions as they are happening. In ecology, shifts in 13C indicate change in diet, and the 13C natural abundance value indicates an equilibrium average of the various food sources. In contrast 15N values have been used to indicate position in trophic level in a food web – the higher the δ 15N value, the higher the trophic level. It also shows categorically that the food source has been consumed and assimilated.
However, Maraun et al. (2011) advised there was a need to investigate the individual feeding channels, but there is a difficulty differentiating basal resources using stable isotopes at natural abundance. The problem with the interpretation of natural abundance studies with the rare stable isotopes 13C and 15N is that it provides an equilibrium and average value of all the sources of food consumed (Moore et al., 2004; Schneider et al., 2004; Pollierer et al., 2009; Schneider and Maraun, 2009). These studies indicated there was sometimes poor resolution between some trophic compartments without more detailed studies. To address these concerns with natural abundance field data, it is necessary to turn toward microcosms to make specific measurements. In these microcosm studies, it is important to consider both competition and predation and not just single species consumption rates, as single species results might not be representative of field data (Wardle and Yeates, 1993). A number of studies employing stable isotopes focused on quantifying the bacterial energy pathway (e.g., Murray et al., 2009; Crotty et al., 2011b, 2012a), and here this study focuses on the fungal pathway. Tracking the energy flux through the fungal-feeding channel within soil food webs is more complicated than the bacterial channel (Crotty et al., 2012b). Fungi are easily shredded during sampling, and are not easily re-introduced in situ, highlighting a need to validate the use of enriched isotope tracers in microcosm settings to study fungal energy pathways.
Energy flux dynamics in multitrophic interactions food webs provide the data to understand biodiversity-ecosystem functioning. Yet, there are few studies portraying this, despite the crucial link these studies would provide in linking trophic interactions with ecosystem function (Barnes et al., 2018). The overall aim of the experiments presented here, were to study the effects of species interactions (competition and predation) on the rate of consumption and assimilation as determined by isotope ratio mass spectrometry (IRMS). To do this the rate of consumption of saprotrophic fungal biomass by soil microarthropods needed to be determined. We tested the hypothesis that mycelium assimilation by the fungivores, measured as body mass enrichment, will be the same for each of the species used. We tested the hypothesis that in the presence of a predator, there is evidence of a trophic cascade, so that consumption of the mycelium is reduced. We also tested the hypothesis that the δ13C and δ15N isotopic signatures of the predator would be affected by its prey’s isotopic composition.
Materials and Methods
General Methods
Soil cores were obtained from a wooded area of Point Pleasant Park in Halifax, Nova Scotia, Canada (latitude 44.621876277, longitude -63.5711053). The microarthropods were extracted from the soil using a Tullgren funnel heat extraction (5 mm mesh) into prepared microcosms. The microcosms were plastic vessels (54 mm diameter, 60 mm height) with a ∼10 mm layer of plaster of Paris mixed with activated charcoal (5:1 by weight, 2 parts powder to 1 part water). Numerous species were found but only some were successfully cultured on fungal mycelia or baker’s yeast in the necessary quantities for the microcosm experiments. Laboratory cultures were established of Lepidocyrtus curvicollis (Collembola: Entomobryidae), Oribatula tibialis (Acari: Oribatida: Oribatulidae), Tyrophagus putrescentiae (Acari: Astigmata: Acaridae), and Hypoaspis acquilifer (Acari: Mesostigmata: Laelapidae). Microarthropods were grown on a mixed diet of baker’s yeast (added every few days) and ground dried nettle as organic matter (added weekly), in darkness at room temperature (about 24–26°C).
Fungal cultures were obtained from the same soil as the microarthropods. Saprotrophic fungi were initially isolated by placing peds of soil on potato dextrose agar (1.5%) with Rose Bengal (0.1%) to suppress bacterial contamination. Sub-cultures of the hyphae, which grew out of the soil, were taken in a series to obtain single fungal species in culture that was identified as an ascomycete by microscopy. The minimal medium was composed of 7 g Na2HPO4, 3 g KH2PO, 0.12 g MgSO4, 0.011 g CaCl2, 0.5 g NaCl2 per liter deionised water, with 1.2% w/v technical agar (number 3), which were combined and sterilized (autoclaved at 121°C for 15 min), before the addition of 2.5 and 1 g, respectively, of filter sterilized 13C-glucose (C6H12O6) and 15N-ammonium chloride (NH4Cl) until the agar substrate was enriched to 99.9% (Crotty et al., 2011a). Reference control cultures were also grown prepared in the same way but the glucose and ammonium chloride were at natural abundance of C and N isotopes. Nitrocellulose disks (45 mm dia. Bio-Rad, Ontario, Canada) were placed on top of the agar and inoculated with the mycelium. Fungal hyphae grow without penetrating the agar to ease the transfer of mycelium to microcosms.
The feeding experiments were conducted in Petri dish (diameter 6 cm) with a <5 mm layer of plaster of Paris and activated charcoal at the bottom. In preparation for each experiment, microarthropods were transferred from culture to a new Petri dish without food, and starved for 48 h before the feeding experiments commenced. Fungal biomass was harvested by lifting the nitrocellulose disks from the agar, and each placed in the center of a fresh Petri dish microcosms. Each microcosm received about 20 mg wet weight of mycelium (weight measured and recorded for stable isotope analysis). We verified in controls the fungal mycelium transferred to the experiment Petri plate microcosms did not grow, as there was no substrate.
In preliminary trials to determine optimal experimental conditions, we established that 20 mg of mycelium in the Petri dish area (approximately a quarter of the nitrocellulose disk) could maintain a stable environment for about 20 Lepidocyrtus or Tyrophagus, with up to 10 Oribatula; more than 10 individuals of the predator Hypoaspis was too crowded in one dish, especially without prey, but 5–10 with prey provided stable microcosms. Oribatula was slow-growing on the cultured fungus, and this species clearly was not a preferred food source. We kept this species so as to provide a range of grazing rates for the competition and predation experiments. These preliminary growth experiments on labeled and un-labeled mycelium were conducted for up to 21 days to establish optimal durations. We established that optimal multispecies microcosm incubations were 5–7 days; shorter incubations (<5 days) showed more variability and lower enrichment; longer periods of up to 2–3 weeks allow recycling nutrients by coprophagy, losing too many prey to the predator, and over-grazing the mycelium. We calculated the amount of consumption using the Petersen and Luxton (1982) equation, Assimilation = consumption – feces.
The experiments were terminated by removing the microarthropods from each dish, grouped by species, killed by submersion in 100% ethanol for <1 min, and samples were processed immediately for analysis by IRMS. The microarthropod filled capsules were dried at 65°C for 48 h, re-weighed, and 13C/12C and 15N/14N ratios measured using Costech ECS4040 elemental analyzer coupled to a Delta V Advantage mass spectrometer (Thermo Fisher Scientific, Finnegan, Germany) with a Conflo IV interface (Thermo Fisher Scientific, Finnegan, Germany), following the methods described in Crotty et al. (2013) for low mass samples.
The fecal pellets and exoskeletons remaining in the microcosms were also collected, counted, added to pre-weighed tin capsules, dried at 65°C for 48 h, re-weighed, and 13C/12C and 15N/14N ratios measured as described above. The amount of fungi consumed by the microarthropods over the incubation period was calculated and expressed as dry weight per individual per day.
Competition Experiments With and Without a Predator
This experiment focused on the effect competition or predation had on feeding rates, using the level of isotope enrichment acquired by each species as a proxy for consumption. Stable isotope enriched mycelium was added to ten different treatments (each treatment replicated five times) in microcosms. The treatments were, four single species monocultures (Lepidocyrtus, Tyrophagus, Oribatula, and Hypoaspis); one with two fungivores in competition (Lepidocyrtus and Oribatula); one with a fungivore and a predator (Lepidocyrtus and Hypoaspis); one with three fungivores (Lepidocyrtus, Tyrophagus, and Oribatula); two microcosm treatments with two fungivores and a predator (either Lepidocyrtus, Oribatula and Hypoaspis or Lepidocyrtus, Tyrophagus, and Hypoaspis); and the final microcosm treatment which had all faunal species in competition / predation with each other (Lepidocyrtus, Tyrophagus, Oribatula, and Hypoaspis) Microcosms were incubated for 7 days, as established above. The number of individuals at the end of each treatment were counted. Following incubation, the Petri dishes were prepared for analysis by mass spectrometry as described above. Individuals of the same species in the same microcosm were bulked due to obtain sufficient biomass for analysis.
Statistical Analysis
Microarthropod effects on mycelium consumption, fecal pellet biomass and isotopic composition where assessed by a general analysis of variance. Where applicable multiple comparison within tables of means were made using the Student-Newman-Keuls test (Sokal and Rohlf, 1995), with the experiment-wise type-1 error rate set at 5%. In order for the mycelium samples to be within the enrichment thresholds of the mass spectrometer a companion reference sample was used to dilute the isotope signal, consequently the isotopic results were back-calculated using the equations of Hauck and Bremner (1976) to ascertain original enrichment levels. All data were analyzed using GenStat (14th Edition, Payne et al., 2011) and are presented as mean ± standard error, or SED (standard error of the difference), unless otherwise stated.
Results
Consumption of Mycelium by Microarthropods
The mean ingestion (consumption rate) and fecal pellet excretion rates over a 21 day single microarthropod incubation with unlabelled mycelium and with 13C and 15N enriched mycelium (except without Oribatula as they were too few) were compared (Table 1). There were no significant differences in microarthropod weight when comparing unlabelled mycelium with highly enriched mycelium between the two experiments after 21 days (P = 0.073). There were no significant differences between the consumption rate of Lepidocyrtus, Oribatula, and Tyrophagus (P = 0.352). The weight of each fecal pellet produced per organism was: Oribatula (0.04 ± 0.012 μg), Tyrophagus (0.33 ± 0.032 μg), and for Lepidocyrtus (0.10 ± 0.020 μg). There were no significant differences in the fecal pellet weight excreted per organism per day in the isotopically enriched cultures compared to the natural abundance measurements (P = 0.15). There was a significant difference between the fungivorous organisms, with Tyrophagus excreting significantly greater quantities than Lepidocyrtus and Oribatula (P = 0.011). The feces of both species showed a lower level of enrichment in 13C than the mycelium at the end of the experiment (Lepidocyrtus APE 3.23 ( ± 0.506) P = 0.033; Tyrophagus APE 2.15 ( ± 0.475) P = 0.010) compared to mycelium APE 5.08 ( ± 0.5735), although not in APE 15N [Lepidocyrtus APE 18.54 ( ± 2.780) P = 0.375; Tyrophagus APE 23.99 ( ± 3.054) P = 0.438] compared to mycelium APE 22.15 ( ± 3.449). These results taken together (Table 1) indicate that the three species of fungivores do not assimilate the same amount of the mycelium consumed, and that differences accumulate over time between the enrichment of fecal pellets and body mass among the fungivores.
Table 1. Weight of microarthropods and mycelium consumed in daily dry-weight per individual (n = 5).
The mesostigmatid predator Hypoaspis on its own with mycelium but without prey lost weight and did not excrete fecal pellets. The three fungivorous microarthropods consumed mycelium at about the same rate ∼5 μg individual-1 day-1, and the differences were not statistically significant. The oribatid Oribatula was the least efficient on this fungus mycelium, excreting ∼31% of the ingested fungus while consuming more than three times its body-weight daily. The astigmatid Tyrophagus was the most efficient, excreting only about 11½% of the ingested mycelium, but it also consumed about three times its body-weight. The collembolan ingested about 1½ times its body weight daily and excreted about 16% of the ingested mycelium. Lepidocyrtus was the only organism to produce exoskeletons in enough numbers to collect and weigh. Each collembolan exoskeleton weighed 1.10 μg (± 0.244) and equated to each individual collembolan shedding its exoskeleton 1.17 (± 0.096) times within a 21 day period.
Analyzing the organisms, their fecal pellets, and the mycelium that remained at the end of the incubations, there were no significant differences between the mycelium 13C and 15N atom% levels and those in the body mass of Lepidocyrtus, Oribatula, or Tyrophagus P > 0.1 for all fungivores.
Survival
In competition experiments with and without the predator, the number of surviving individuals at the end of the incubation in each microcosm is informative (Table 2). Individuals are lost to predation and poor adaptation to the incubation conditions. The predator Hypoaspis survival was about 75% in monoculture without prey. Over half of the Lepidocyrtus and Tyrophagus survived in monoculture on the mycelium. Surprisingly, even though the mycelium did not provide a good food source to Oribatula, there was only a 12% loss of individuals in monoculture (Table 2). The highest survival rates of predators was when there were two or more prey species, although only some combinations are statistically significant. The collembolan Lepidocyrtus did not ameliorate survival of the predator. The Oribatula losses were greatest when competing for the mycelium with Lepidocyrtus, with indiscernible predation effect. Tyrophagus losses were least with all three fungivores present and the predator. The Lepidocyrtus microcosms show a distinct increase in losses in the presence of the predator, but this observation is nuanced by the observation (above) that the collembolan did not improve survival of the predator. Visually, at the microscope, collembolan were better at escaping, and it was more effort for the predator to capture one.
Table 2. Loss (%) of individuals in each treatment.
Species Interactions Effect on Consumption
The results of the competition-predation microcosm incubation experiments with stable isotope enrichment are summarized in Figure 1. For all these graphs the microarthropod 13C and 15N value at the beginning of each experiment is at the 0 mark of both axis, and we plot the enrichment transferred to the microarthropods in atom percent enrichment (APE). It does not represent the natural abundance value in nature, as these organisms were cultivated in the laboratory for about 2 years.
Figure 1. Isotopic label of microarthropods through consumption of fungal mycelium for (A) Lepidocyrtus (Collembola), (B) Oribatula (Oribatida), (C) Tyrophagus (Astigmata), and (D) Hypoaspis (Mesostigmata) in microcosms with different combination of competitors with/out predators (n = 5). The Hypoaspis graph (D) is at a different scale from the others. The microarthropod 13C and 15N value at the beginning of each experiment is at the 0 mark of both axis, and we plot the enrichment transferred to the microarthropods in atom percent enrichment (APE). It does not represent the natural abundance value in nature, as these organisms were cultivated in the laboratory prior to the feeding experiment.
Lepidocyrtus (Figure 1A) consumption of mycelium differed significantly depending on competition/predation for both isotopes (APE 13C P = 0.010 and APE 15N P = 0.018). Consumption was numerically greatest when incubated as a single microcosm. When Lepidocyrtus was incubated solely with the predator Hypoaspis (black triangle, Figure 1A) there was very little change in 13C or 15N value, indicating consumption did not decrease significantly, although survival reduced to 33% (Table 2). When in competition with Oribatula (open circle, Figure 1A) or when the microcosms contained two competitors, the value is not statistically different than when it was on its own, indicating that Oribatula (and Tyrophagus) in fact did not pose a competition threat to Lepidocyrtus. When Lepidocyrtus is incubated with the predator Hypoaspis the 13C value is significantly reduced (less consumption) but not the 15N value (no change in trophic level, and as results above). When both competitors are present (clear triangle), the three fungivores interfere with each other’s grazing so as to reduce the 15N value and the 13C value compared to Lepidocyrtus alone or with Oribatula. In the presence of the predator and Oribatula as the only competitor, Lepidocyrtus (shaded square) values are not statistically different than with Oribatula alone, indicating it is relatively unbothered by predation in this combination. In the absence of the predator but both competitors (clear triangle), there is an intermediate amount of grazing. It clearly shows the difference in competition imposed on Lepidocyrtus by Tyrophagus but not by Oribatula. Without Oribatula (clear circle), but with competition from Tyrophagus and predation from Hypoaspis (clear square), the 13C value is reduced but not the 15N value when compared to Lepidocyrtus alone, indicating reduced mycelium consumption. Curiously, when the predator and both competitors are present with Lepidocyrtus (shaded diamond), the value overlaps that of competition with Oribatula alone (clear circle). This value is statistically the same as Lepidocyrtus with the predator alone for 15N, and for Lepidocyrtus alone for 13C. However, when the microcosm included Tyrophagus as a competitor and predators (open square, black diamond, Figure 1A), significantly less consumption occurred. At the microscope in these microcosms, Lepidocyrtus seemed unbothered or less-so by the predator, which preferred to chase after the two other fungivores, and thus alleviating competition from the two other fungivores. Feces enrichment levels were also measured for all the different treatments containing Lepidocyrtus and there were no significant differences in enrichment levels dependent on species interaction (P = 0.634 for 13C and P = 0.498 for 15N).
The amount of isotopic enrichment was relatively lower in Oribatula (Figure 1B), as this fungus was not a preferred food source for this species. However, there were no significant differences between treatments in enrichment levels for either 13C or 15N (P = 0.194 and P = 0.276, respectively) for Oribatula. Comparison of enrichment through grazing between Lepidocyrtus and Oribatula across the different treatments showed significant differences in enrichment levels in both 13C and 15N (P < 0.001 for both).
Tyrophagus showed a different pattern of response to grazing competition and predation (Figure 1C). The level of enrichment from the mycelium between Lepidocyrtus and Tyrophagus were the same (P < 0.001 for both isotopes). As a single species microcosm and in the presence of two competitors, the three species microcosm had significantly greater for both 13C and 15N when compared to Tyrophagus with predation (P < 0.001 for both). The APE 15N result was significantly highest when two competitors were present (open circle, Figure 1C), this could be through coprophagy on fecal pellets of one or both competitors. In the presence of a predator and either Lepidocyrtus or Lepidocyrtus and Oribatula as competitors, the amount of grazing significantly decreased and the isotope APE were greatly reduced for both isotopes. In this composition of the organisms (Figure 1C, black/open triangle) grazing on the mycelium is interfered with through effective grazing competition, combined with organisms trying to escape predation, and by the predation itself.
The level of enrichment obtained by the predator Hypoaspis was significantly lower than all other species (P < 0.001 for both 13C and 15N) by one order of magnitude. This is to be expected with the predator being one trophic level above the fungi consumers. When Hypoaspis was in a single species monoculture, negligible amounts of enrichment occurred for both 13C and 15N, however, significant quantities of enrichment were obtained when in microcosms with prey species (P = 0.007 for 13C and P < 0.001 for 15N) (Figure 1D). There is a trend for greater levels of 15N enrichment when in a microcosm with Tyrophagus in comparison to Lepidocyrtus or Oribatula (open triangle / black square).
Comparing the microcosms with all four organisms together (Figure 1) shows the variety of consumption preferences detected between the predator (Mesostigmata) on each of its prey, and of each of the fungivores (Collembola, Astigmata, Oribatida) on the mycelium. The predator (Table 2 and Figure 1D) far preferred the Astigmata Tyrophagus as prey, compared to the oribatid Oribatula, or the collembolan Lepidocyrtus. This is consistent with our results above, and with our visual observations. In the combination with all four species present, Lepidocyrtus consumption of mycelium (Figure 1A) was not affected, as the other fungivores were preferred alternative prey for the predator, and Oribatula was not an effective competitor for mycelium consumption. Tyrophagus was an effective competitor for Lepidocyrtus as both consumed the mycelium at the same rate, but it could barely eat (Figure 1C) because of interference from Lepidocyrtus, and preferential predation by Hypoaspis. In this combination, Oribatula was a secondary prey for Hypoaspis but non-etheless preyed upon. In addition, it was in competition with both Lepidocyrtus and Tyrophagus that both preferred the mycelium, which was less palatable to Oribatula. Thus, Oribatula and Tyrophagus were less enriched in this four species combination.
Discussion
We selected an ascomycete culture that would sustain three fungivorous microarthropods, and one predatory microarthropod that would prey on the fungivores. We quantified rates of hyphae consumption using 13C and 15N stable isotopes enriched mycelium. We compared microcosms with and without predation, and with or without competition for changes in hyphae consumption and predation rates.
Precision and Sensitivity
We showed that standard mass spectrometry instruments could perform at much higher sensitivity by modifying the chemistry of the combustion stage in the sample chamber (Crotty et al., 2013). The IRMS results in that paper showed that 1–10 individual microarthropods were sufficient for reliably accurate measurements. The microcosms in this study consisted of 10–60 individuals at the start of treatments. The resolution of the experiments (Figure 1) are the mean values with five replications. The exception was the Oribatula (Figure 1B) where low consumption showed the pattern but without statistical validity. However, the spread of the estimate of error in these experiments (Figure 1) are reasonably low for IRMS microcosms.
The objective was to obtain measurements of fungal consumption, assimilation and foraging behavior of each microarthropod when affected by competition and/or predation. We understand that in these microcosms we have a flat surface that does not mimic the soil three dimensional structure. Within the complexity of the soil matrix there would be a diversity of more or less palatable species to choose from, and behavior of escape and competition interference could be different. Conceivably, in the microcosm the Furca (jumping organ, unique to Collembola), may have reduced the time needed to avoid predation (Hopkin, 2002) and resume fungal consumption; and in contrast in the soil matrix, Collembola could be harder to find but less able to escape. In addition, the fungivore consumers were fed a single species of fungus continuously that does not represent a more mixed diet in the soil. Our 7 days incubation period (from trials of up to 21 days) was chosen to provide sufficient levels of consumption to occur for stable isotope enrichment in 2nd and 3rd trophic levels, while providing the number of individuals necessary for isotopic analysis.
Differential Consumption-Predation
In natural abundance stable isotope studies, due to isotope fractionation during metabolism, δ13C values indicate the source of food and the amount of consumption while the δ15N values indicate the amount consumed and the trophic level of the species (Tiunov, 2007; Heijboer et al., 2018). In environmental samples with natural abundances, the rare stable isotopes ratio to common isotopes (13C and 15N to 12C and 14N) reflect equilibrium values of the isotopes from diverse food sources. That is not necessarily the case in microcosms with species in culture and then in experimental treatments, even though the microarthropods were kept on the same mycelium for many months prior to the treatment. The natural abundance ratios of the three fungivores growing on the same mycelium for several months showed variations in isotope levels between the different microarthropods. If all three species had the same diet and consumed the same amount of the mycelium, then Lepidocyrtus, Oribatula and Tyrophagus would show the same isotope ratios, but that is not the case. The differences are due to probably an array of reasons, such as differences in development stage or exoskeleton - coprophagy / fecal pellet consumption, and selective metabolism (Semenina and Tiunov, 2011; Maraun et al., 2011; Hatch, 2012; Potapov et al., 2014).
The δ13C values of each fungivore in monoculture (Figures 1A–C, dark circles) reflects the equilibrium amount and value of δ13C when undisturbed. During isotopic enrichment treatment incubation, the rare stable isotope is accumulated and the graphs show enrichment for both isotopes compared to the starting point. However, with increasing disturbance from competition and predation in mixed cultures, the δ13C values show less enrichment (13C APE) moving closer to the zero value (the equilibrium value at the start).
There was a visible trend that Tyrophagus was a preferred prey for Hypoaspis, as those microcosms provided a greater amount of enrichment to the predator. It is often stated that oribatida live relatively sheltered from predation due to their size and sclerotisation (Peschel et al., 2006; Schneider and Maraun, 2009) which are reasonable assumptions based on the low predation rate on Oribatula (Figures 1B,D) although the greatest losses of Oribatula were in the microcosms that contained predators (Table 2). Moreover, Oribatula tibialis is known to have a defensive cyanogenic aromatic ester (mandelonitrile hexanoate) which would further explain the low predation rate (Brückner et al., 2017). These experiments provide evidence of differential consumption and differential prey preference, so that in the soil optimal foraging would contribute to the behavior of both the fungivores and predators (Hassall et al., 2006; Schneider and Maraun, 2009).
The natural abundance δ15N levels were similar between the predatory mite Hypoaspis and two fungivores Lepidocyrtus and Oribatula before commencement of the competition/predation experiment. However, when Hypoaspis was in microcosms with a prey species, enrichment from fungal mycelium was traceable to a third trophic level. This exemplifies the caution that is needed when assessing isotope levels in relation to trophic level. Based on natural abundance values from field data increases of 3.4% (Post, 2002) or 1–2% (McCutchan et al., 2003; Illig et al., 2005) have been promoted as indicators of trophic level differences in food webs across ecosystems (see also Moore et al., 2004). Thus, we cannot ignore that these organisms could have consumed eggs, exoskeletons and fecal pellets, engaged in cannibalism, or even predation on juveniles as alternative sources of food in the microcosms, thus appearing to be at a higher trophic level. Some studies have even found Mesostigmata to act as more than just predators by engaging is some opportunist omnivory (Beaulieu, 2012) which would lower its apparent trophic level in natural abundance data.
The literature is rich with examples of microarthropods having a more mixed diet in the field than typically construed, based on casual laboratory observations and from observations in microcosms (Chahartaghi et al., 2005; Erdmann et al., 2007). Other modes of feeding such as coprophagy using exoskeletons and fecal pellets, on eggs, by cannibalism or predation on juveniles have been observed elsewhere (Ponge, 1991; Briones et al., 1999; Ladygina et al., 2008; Endlweber et al., 2009; Fiera, 2014; Feng et al., 2019). For example, Lepidocyrtus has a hemimetabolous lifecycle compared to Oribatula and Tyrophagus which are holometabolous, so that different levels of fractionation may have occurred (Spence and Rosenheim, 2005). Tyrophagus is also known to be somewhat omnivorous (Behan-Pelletier, 1999) and therefore, could have potentially found an alternative food source within the microcosms, consuming eggs or immature hatchlings. Not implying that all species are opportunistically or a little omnivorous, we must nonetheless recognize that not all are as specialized as the literature might suggest. These additional sources of consumption that must be considered would affect trophic level interpretations in field data, as they do in microcosms. Microcosm observations then become useful in identifying these cases, and in quantifying energy fluxes in interactions food webs.
Top-Down Control and Competition Interference
We have shown the fungivorous microarthropods investigated had a voracious appetite, consuming more than their own body weight per day (Table 1). In our preliminary microcosm runs, with incubation of up to 21 days, to avoid the fungivores grazing down the mycelium we provided sufficient non-growing mycelium so as to avoid bottom-up resource limitation on the microcosm. We observed several parameters of optimal foraging theory to be at play in our resource competition interference and predation food webs. One was the interference by competition for the mycelium by other fungivores (Figure 1, and see section “Species Interactions Effect on Consumption”), another was a differential preference for, consumption and excretion of the consumed mycelium biomass. There were more complicated behavioral effects, such as the mycelium consumption of Tyrophagus (Figure 1C) reduced by competition from Lepidocyrtus especially in the presence of the predator Hypoaspis as Tyrophagus was a preferred prey. Last, we speculated that the detailed interpretation of the stable isotope results (Figure 1) suggest a more mixed diet that included fecal pellets, exoskeleton, or even juveniles. Excretion by Mesostigmata is mainly fluid with no fecal pellets being formed (Koehler, 1999), but Collembola and Oribatida fecal pellets are important in soil aggregate formation, as a food source for other soil invertebrates, and sometimes contribute significantly to soil structure (Behan-Pelletier, 2003; Maas et al., 2015).
Although the three-dimensional nature of the soil pore space is dissimilar to the Petri dish microcosms, our results convey the pivotal role fungivores must play in nutrient cycling, when considering the amount of fungal mycelium consumed and excreted over time, as both groups represent major trophic and functional components of the soil environment (Leake et al., 2003; Crowther and A’Bear, 2012). The interactive effects of biotic (microarthropods and mycelium) and abiotic (habitat) are central in resolving the functional and ecological responses of the soil food web and this study provides quantified values for use in global models. However, in the natural environment abiotic factors or resource accessibility may affect populations more than competition. Similar studies by others will provide additional comparative data to better understand the complex soil food web interactions, including quantifying the fungivory pathway in detrital food webs.
Here, we have shown methods to ascertain the levels of consumption, competition and predation within the soil food web, knowledge gaps that have been highlighted historically, (Moore et al., 1988; Polis, 1994) but until now have not been quantified. Our results suggest the importance of trophic cascades, these have been modeled, and observed in low productivity systems (de Ruiter et al., 1995; Moore and De Ruiter, 2000; Moore, 2018) but have not been shown empirically until now. Studies such as this set of microcosm studies described here with stable isotope tracers, provide both the required resolution and microscope observations that allow collection of the necessary data for predictive model simulations to be developed.
Data Availability
The datasets generated for this study are available on request to the corresponding author.
Author Contributions
This manuscript was carried out by FC during a postdoctoral year with SA. SA wrote the research proposal, obtained the grant, and the overall design of the experiments. FC obtained the measurements and conducted the analysis and calculations. FC wrote the first version of this manuscript. It has since been rewritten, with additional statistical analysis by SA.
Funding
This work was supported through NSERC grant 249889 to SA.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Myles Stocki for operating the mass spectrometer and maintaining its accuracy and precision to the extent necessary for this work with low biomass samples.
References
Barnes, A. D., Jochum, M., Lefcheck, J. S., Eisenchauer, N., Scherber, C., O’Connor, M. I., et al. (2018). Energy flux: the link between multitrophic biodiversity and ecosystem functioning. Trends Ecol. Evol. 33, 186–197. doi: 10.1016/j.tree.2017.12.007
Beaulieu, F. (2012). Saproxyly in predatory mites? Mesostigmata in decaying log habitats versus litter in a wet eucalypt forest, Tasmania, Australia. Int. J. Acarol. 38, 313–323. doi: 10.1080/01647954.2011.647072
Behan-Pelletier, V. M. (1999). Oribatid mite biodiversity in agroecosystems: role for bioindication. Agric. Ecosyst. Environ. 74, 411–423. doi: 10.1016/b978-0-444-50019-9.50023-6
Behan-Pelletier, V. M. (2003). Acari and Collembola biodiversity in Canadian agricultural soils. Can. J. Soil Sci. 83, 279–288. doi: 10.4141/s01-063
Bradford, M. A., Strickland, M. S., Devore, J. L., and Maerz, J. C. (2012). Root carbon flow from an invasive plant to belowground foodwebs. Plant Soil 359, 233–244. doi: 10.1007/s11104-012-1210-y
Briones, M. J. I., Ineson, P., and Sleep, D. (1999). Use of 13C to determine food selection in collembolan species. Soil Biol. Biochem. 31, 937–940. doi: 10.1016/s0038-0717(98)00179-5
Brückner, A., Raspotnig, G., Wehner, K., Meusinger, R., Norton, R. A., and Heethoff, M. (2017). Storage and release of hydrogen cyanide in a chelicerate (Oribatula tibialis). Proc. Natl. Acad. Sci. U.S.A. 114, 3469–3472. doi: 10.1073/pnas.1618327114
Chahartaghi, M., Langel, R., Scheu, S., and Ruess, L. (2005). Feeding guilds in Collembola based on nitrogen stable isotope ratios. Soil Biol. Biochem. 37, 1718–1725. doi: 10.1016/j.soilbio.2005.02.006
Crotty, F. V., Adl, S. M., Blackshaw, R. P., and Murray, P. J. (2012a). Protozoan pulses unveil their pivotal position within the soil food web. Microb. Ecol. 63, 905–918. doi: 10.1007/s00248-011-9956-y
Crotty, F. V., Adl, S. M., Blackshaw, R. P., and Murray, P. J. (2012b). Review: using stable isotopes to differentiate trophic feeding channels within soil food webs. J. Eukaryot. Microbiol. 59, 520–526. doi: 10.1111/j.1550-7408.2011.00608.x
Crotty, F. V., Blackshaw, R. P., and Murray, P. J. (2011a). Differential growth of the fungus Absidia cylindrospora on 13C/15N-labelled media. Rapid Commun. Mass Spectrom. 25, 1479–1484. doi: 10.1002/rcm.4888
Crotty, F. V., Blackshaw, R. P., and Murray, P. J. (2011b). Tracking the flow of bacterially derived 13C and 15N through soil faunal feeding channels. Rapid Commun. Mass Spectrom. 25, 1503–1513. doi: 10.1002/rcm.4945
Crotty, F. V., Stocki, M., Knight, J. D., and Adl, S. M. (2013). Improving accuracy and sensitivity of isotope ratio mass spectrometry for δ13C and δ15N values in very low mass samples for ecological studies. Soil Biol. Biochem. 65, 75–77. doi: 10.1016/j.soilbio.2013.04.020
Crowther, T. W., and A’Bear, A. D. (2012). Impacts of grazing soil fauna on decomposer fungi are species-specific and density-dependent. Fungal Ecol. 5, 277–281. doi: 10.1016/j.funeco.2011.07.006
de Ruiter, P. C., Neutel, A.-M., and Moore, J. C. (1995). Energetics, patterns of interaction strength, and stability in real ecosystems. Science 269, 1257–1260. doi: 10.1126/science.269.5228.1257
Elfstrand, S., Lagerlöf, J., Hedlund, K., and Mårtensson, A. (2008). Carbon routes from decomposing plant residues and living roots into soil food webs assessed with 13C labelling. Soil Biol. Biochem. 40, 2530–2539. doi: 10.1016/j.soilbio.2008.06.013
Endlweber, K., Ruess, L., and Scheu, S. (2009). Collembola switch diet in presence of plant roots thereby functioning as herbivores. Soil Biol. Biochem. 41, 1151–1154. doi: 10.1016/j.soilbio.2009.02.022
Erdmann, G., Otte, V., Langel, R., Scheu, S., and Maraun, M. (2007). The trophic structure of bark-living oribatid mite communities analysed with stable isotopes (N-15, C-13) indicates strong niche differentiation. Exp. Appl. Acarol. 41, 1–10. doi: 10.1007/s10493-007-9060-7
Feng, L., Chang, L., Zhang, S., Zhu, X., Adl, S., and Wu, D. (2019). What is the carcass-usage mode of the collembola? A case study of Entomobrya proxima in the laboratory. Insects 10:67. doi: 10.3390/insects10030067
Fiera, C. (2014). Application of stable isotopes and lipid analysis to understand trophic interactions in springtails. North West. J. Zool. 10, 227–235.
Hassall, M., Adl, S., Berg, M., Griffiths, B., and Scheu, S. (2006). Soil fauna-microbe interactions: towards a conceptual framework for research. Eur. J. Soil Biol. 42, S54–S60.
Hatch, K. A. (2012). “The use and application of stable isotope analysis to the study of starvation, fasting, and nutritional stress in animals,” in Comparative Physiology of Fasting, Starvation, and Food Limitation. ed. M. McCue (Berlin: Springer), 337–364. doi: 10.1007/978-3-642-29056-5_20
Hauck, R. D., and Bremner, J. M. (1976). Use of tracers for soil and fertilizer nitrogen research. Adv. Agron. 28, 219–266. doi: 10.1016/s0065-2113(08)60556-8
Heijboer, A., Ruess, L., Traugott, M., Jousset, A., and de Ruiter, P. C. (2018). “Empirical methods of identifying and quantifying trophic interactions for constructing soil food-web models,” in Adaptive Food Webs: Stability and Transitions of Real and Model Ecosystems. eds J. C. Moore, P. C. de Ruiter, K. S. McCann, and V. Wolters (Cambridge: Cambridge Univ Press), 257–280.
Hopkin, S. P. (2002). “Collembola,” in Encyclopedia of Soil Science. ed. R. Lal (New York, NY: Marcel Dekker), 207–210.
Illig, J., Langel, R., Norton, R. A., Scheu, S., and Maraun, M. (2005). Where are the decomposers? Uncovering the soil food web of a tropical montane rain forest in southern Ecuador using stable isotopes (N-15). J. Trop. Ecol. 21, 589–593. doi: 10.1017/s0266467405002646
Jonas, J. L., Wilson, G. W. T., White, P. M., and Joern, A. (2007). Consumption of mycorrhizal and saprophytic fungi by Collembola in grassland soils. Soil Biol. Biochem. 39, 2594–2602. doi: 10.1016/j.soilbio.2007.05.004
Koehler, H. H. (1999). Predatory mites (Gamasina, Mesostigmata). Agric. Ecosyst. Environ. 74, 395–410. doi: 10.1016/b978-0-444-50019-9.50022-4
Koukol, O., Mourek, J., Janovsky, Z., and Cerna, K. (2009). Do oribatid mites (Acari: Oribatida) show a higher preference for ubiquitous vs. specialized saprotrophic fungi from pine litter?. Soil Biol. Biochem. 41, 1124–1131. doi: 10.1016/j.soilbio.2009.02.018
Ladygina, N., Caruso, T., and Hedlund, K. (2008). Dietary switching of collembola in grassland soil food web. Soil Biol. Biochem. 40, 2898–2903. doi: 10.1016/j.soilbio.2008.08.012
Leake, J. R., Johnson, D., Donnelly, D. P., Muckle, G. E., Boddy, L., Read, D. J., et al. (2003). “Networks of power and influence: the role of mycorrhizal mycelium in controlling plant communities and agroecosystem functioning,” in Proceedimgs of the 4th International Conference on Mycorrizae (ICOM 2003). (Montreal), 1016–1045. doi: 10.1139/b04-060
Maas, S., Caruso, T., and Rillig, M. (2015). Functional role of microarthropods in soil aggregation. Pedobiologia 58, 59–63. doi: 10.1016/j.pedobi.2015.03.001
Maraun, M., Erdmann, G., Fischer, B. M., Pollierer, M. M., Norton, R. A., Schneider, K., et al. (2011). Stable isotopes revisited: their use and limits for oribatid mite trophic ecology. Soil Biol. Biochem. 43, 877–882. doi: 10.1016/j.soilbio.2011.01.003
Maraun, M., Martens, H., Migge, S., Theenhaus, A., and Scheu, S. (2003). Adding to ’the enigma of soil animal diversity’: fungal feeders and saprophagous soil invertebrates prefer similar food substrates. Eur. J. Soil Biol. 39, 85–95. doi: 10.1016/s1164-5563(03)00006-2
McCutchan, J. H., Lewis, W. M., Kendall, C., and Mcgrath, C. C. (2003). Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos 102, 378–390. doi: 10.1371/journal.pone.0140946
Moore, J. C. (2018). Predicting tipping-points in complex environmental systems. Proc. Natl. Acad. Sci. U.S.A. 115, 635–636. doi: 10.1073/pnas.1721206115
Moore, J. C., Berlow, E. L., Coleman, D. C., de Ruiter, P. C., Dong, Q., Hastings, A., et al. (2004). Detritus, trophic dynamics and biodiversity. Ecol. Lett. 7, 584–600. doi: 10.1111/j.1461-0248.2004.00606.x
Moore, J. C., and De Ruiter, P. C. (2000). “Invertebrates in detrital food web along gradients of productivity,” in Invertebrates as Webmasters in Ecosystems. eds D. C. Coleman and P. Hendrix (Wallingford: CABI), 161–184. doi: 10.1079/9780851993942.0161
Moore, J. C., McCann, K., and de Ruiter, P. C. (2005). Modeling trophic pathways, nutrient cycling, and dynamic stability in soils. Pedobiologia 49, 499–510. doi: 10.1016/j.pedobi.2005.05.008
Moore, J. C., Walter, D. E., and Hunt, H. W. (1988). Arthropod regulation of micro- and mesobiota in below-ground detrital food webs. Annu. Rev. Entomol. 33, 419–439. doi: 10.1146/annurev.ento.33.1.419
Murray, P. J., Clegg, C. D., Crotty, F. V., De La Fuente Martinez, N., Williams, J. K., and Blackshaw, R. P. (2009). Dissipation of bacterially derived C and N through the meso- and macrofauna of a grassland soil. Soil Biol. Biochem. 41, 1146–1150. doi: 10.1016/j.soilbio.2009.02.021
Payne, R. W., Murray, D. A., Harding, S. A., Baird, D. B., and Soutar, D. M. (2011). Introduction to GenStat for Windows. 14th Edn. Hemel Hempstead: VSN International.
Peschel, K., Norton, R. A., Scheu, S., and Maraun, M. (2006). Do oribatid mites live in enemy-free space? Evidence from feeding experiments with the predatory mite Pergamasus septentrionalis. Soil Biol. Biochem. 38, 2985–2989. doi: 10.1016/j.soilbio.2006.04.035
Petersen, H., and Luxton, M. (1982). A comparative-analysis of soil fauna populations and their role in decomposition processes. Oikos 39, 287–388.
Peterson, B. J., and Fry, B. (1987). Stable isotopes in ecosystem studies. Annu. Rev. Ecol. Syst. 18, 293–320. doi: 10.1146/annurev.es.18.110187.001453
Polis, G. A. (1994). Food webs, trophic cascades, and community structure. Aust. J. Ecol. 19, 121–136. doi: 10.1111/j.1442-9993.1994.tb00475.x
Polis, G. A., Sears, A. L., Huxel, G. R., Strong, D. R., and Maron, J. (2000). When is a trophic cascade a trophic cascade? Trends Ecol. Evol. 15, 473–475. doi: 10.1016/s0169-5347(00)01971-6
Pollierer, M. M., Langel, R., Scheu, S., and Maraun, M. (2009). Compartmentalization of the soil animal food web as indicated by dual analysis of stable isotope ratios (15N/14N and 13C/12C). Soil Biol. Biochem. 41, 1221–1226. doi: 10.1016/j.soilbio.2009.03.002
Ponge, J. F. (1991). Food Resources and diets of soil animals in a small area of Scots pine litter. Geoderma 49, 33–62. doi: 10.1016/0016-7061(91)90090-g
Post, D. M. (2002). Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology 83, 703–718. doi: 10.1890/0012-9658(2002)083%5B0703:usitet%5D2.0.co;2
Potapov, A. M., Semenyuk, I. I., and Tiunov, A. V. (2014). Seasonal and age-related changes in the stale isotope composition (15N/14N and 13C/12C) of millipedes and collembolans in a temperate forest soil. Pedobiologia 57, 215–222. doi: 10.1016/j.pedobi.2014.09.005
Ruess, L., Tiunov, A., Haubert, D., Richnow, H. H., Haggblom, M. M., and Scheu, S. (2005). Carbon stable isotope fractionation and trophic transfer of fatty acids in fungal based soil food chains. Soil Biol. Biochem. 37, 945–953. doi: 10.1016/j.soilbio.2004.09.015
Schneider, K., and Maraun, M. (2005). Feeding preferences among dark pigmented fungal taxa (”Dematiacea”) indicate limited trophic niche differentiation of oribatid mites (Oribatida, Acari). Pedobiologia 49, 61–67. doi: 10.1016/j.pedobi.2004.07.010
Schneider, K., and Maraun, M. (2009). Top-down control of soil microarthropods - evidence from a laboratory experiment. Soil Biol. Biochem. 41, 170–175. doi: 10.1016/j.soilbio.2008.10.013
Schneider, K., Migge, S., Norton, R. A., Scheu, S., Langel, R., Reneking, A., et al. (2004). Trophic niche differentiation in soil microarthropods (Oribatida, Acari): evidence from stabel isotope ratios (15N/14N). Soil Biol. Biochem. 36, 1769–1774. doi: 10.1016/j.soilbio.2004.04.033
Semenina, E. E., and Tiunov, A. V. (2011). Trophic fractionation ([Delta]15N) in Collembola depends on nutritional status: a laboratory experiment and mini-review. Pedobiologia 54, 101–109. doi: 10.1016/j.pedobi.2010.10.004
Sokal, R. R., and Rohlf, F. J. (1995). Biometry - The Principles and Practice of Statistics in Biological Research. New York, NY: Freeman and Company.
Spence, K. O., and Rosenheim, J. A. (2005). Isotopic enrichment in herbivorous insects: a comparative field-based study of variation. Oecologia 146, 89–97. doi: 10.1007/s00442-005-0170-9
Thiele-Bruhn, S., Bloem, J., De Vries, F. T., Kalbitz, K., and Wagg, C. (2012). Linking soil biodiversity and agricultural soil management. Curr. Opin. Environ. Sustain. 4, 523–528. doi: 10.1016/j.cosust.2012.06.004
Tiunov, A. V. (2007). Stable isotopes of carbon and nitrogen in soil ecological studies. Biol. Bull. 34, 395–407. doi: 10.1134/s1062359007040127
Keywords: fungivory, microarthropods, nutrient cycling, stable isotopes, trophic interactions
Citation: Crotty FV and Adl SM (2019) Competition and Predation in Soil Fungivorous Microarthropods Using Stable Isotope Ratio Mass Spectrometry. Front. Microbiol. 10:1274. doi: 10.3389/fmicb.2019.01274
Received: 18 March 2019; Accepted: 22 May 2019;
Published: 07 June 2019.
Edited by:
Stefan Geisen, Netherlands Institute of Ecology (NIOO-KNAW), NetherlandsReviewed by:
Mark Maraun, University of Göttingen, GermanyKerstin Heidemann, University of Göttingen, Germany
Copyright © 2019 Crotty and Adl. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sina M. Adl, c2luYS5hZGxAdXNhc2suY2E=
†Present address: Felicity V. Crotty, Royal Agricultural University, Cirencester, United Kingdom