More stories

  • in

    Development of a treatment for water contaminated with Cr (VI) using cellulose xanthogenate from E. crassipes on a pilot scale

    Analysis of FTIRUnderstanding the functional groups involved in the biosorption of toxic metals is essential to elucidate the mechanism of this process. Groups such as carboxylic, hydroxyl and amine are among the main responsible for the absorption of metals by cellulose34 In the Fig. 1, show the FTIR of ECx.Figure 1FTIR of ECx before and after of adsorptions of Cr (VI).Full size imageAccording to13 the bandwidth at 3000–3600 cm−1 corresponds to bonds related to the -OH group. These hydrogen bonds are useful tools for cation exchange with heavy metals. This evidenced in the color spectrum (dark green) that represents an ECx sample with attached Cr (VI) after the adsorption process, where the stretching of the (OH) group lost part of its extension. The change observed in the peak from 3420 cm−1 of ECx to 3440 cm−1 in ECx-Cr indicates that these groups have a participation in the bond with the Cr (VI) ions. The variation of bands in the peak of the amines after adsorption confirms the participation of these groups in the adsorption process. This result confirmed by the ion exchange evaluation experiment discussed later in section SEM–EDX.The change in peak 3280, after Cr (VI) adsorption, indicates that EC removed Cr (VI) based on interaction with (OH), part of (OH) lost due to formation of vibrations of ascension O–Cr. Also, after Cr (VI) biosorption on ECx, the peak of the EC-S group is shifted to 590. This can be explained by surface complexation or ion exchange35.In general, comparable results reported in the literature for cellulose in the absorption of other toxic metals, as for other cellulose-derived biosorbentes in the removal of Cr (VI) ions36.One way to corroborate the information presented in the FTIR measurements is through SEM images since with these images it is possible to observe the distribution of the reagents in the ECx biomass treatment and subsequently the Cr (VI) adsorption process.SEM–EDXFigure 2 shows the micrographs obtained for the biomass before (a) the adsorption of Cr (VI), in addition to showing the distribution of the different biomass chemical modifications in (b) and in (c) it shows the distribution of chromium around all biomasses.Figure 2Biomass before (a) Cr (VI) adsorption, biomass chemical modifications in (b) and shows the distribution of chromium around the whole biomass (c).Full size imageFrom Fig. 2a, it can see that the biomass has a very irregular rough surface, with macropores and cracks. Many of these irregularities may associated with damage caused by the delignification process of E. crassipes cellulose with NaOH14. In Fig. 2b it is possible to visualize the components of the cellulose xanthogenate, coming from sodium, distributed throughout the biomass, a result like that reported in other studies35 The colored dots represent the elements in the samples, green dots represent carbon, red dots represent oxygen, and yellow dots represent the places where sodium lodged.Table 2 shows that, in addition to carbon and oxygen, the element with the greatest presence in the composition of pure waste is sodium and sulfur from the xanthogenate cellulose transformation process. Table 2 shows the physicochemical characterization of the ECx sample, through EDS.Table 2 Features of sample of ECx.Full size tableCellulose xanthogenate, is one of the cellulose transformations to improve the adsorption performance of heavy metals, this compound produced from dry and ground biomass, mixing with sodium hydroxide (NaOH) to remove lignin, creating alkaline biomass, then disulfide (CS2) added13,14. (CS2) reacts with hydratable hydroxycellulose, forming C-SNa complexes; these are responsible for the cation exchange with heavy metals. Metal ions enter the interior of E. crassipes with (CS2), exchanging with Na36,37.The SEM morphology of ECx and coupled with the high content of sulfides (7.3%) determined by the spectrum in Table 2, it further confirms that xanthate groups are successfully grafted onto the biomass of E. crassipes, and Fig. 3 represents this information based on13,36,37,38.Figure 3Prototype.Full size imageExchange biochemistry is usually identified as the main mechanism for the adsorption of metals in cellulose and its derivatives35 and through the evaluation of EDS this process could verify. Similar observations were made by36 where the adhesion of Cr (VI) in this biomass was observed. Also, in xanthogenate cellulose processes, the adhesion of Pb (II) to this type of biomass verified, concluding that this cellulose is important in the removal of heavy metals from water13.The SEM morphology of ECx with Cr (VI) coupled with the high content of sulfides determined by the spectrum in Table 3, was the determinate for the chemisorption’s of Cr (VI). The mechanism of Cr (VI) sorption by cellulose xanthate is:$$left[ {{4}left( {{text{C}}_{{6}} {text{H}}_{{{12}}} {text{O}}_{{6}} } right)} right]*{text{2CS}}_{{2}} {text{Na }} + {text{ Cr}}_{{2}} {text{O}}_{7}^{ – 2} to left{ {left[ {{4}left( {{text{C}}_{{6}} {text{H}}_{{5}} {text{O}}_{{6}} } right)} right] , *{text{2CS}}_{{2}} } right}*{mathbf{Cr}}_{{mathbf{2}}} + {text{Na}} + {text{7H}}_{{2}} {text{O}}$$where [4(C6H12O6)] *2CS2Na represents the xanthogenate biomass, and Cr2O7–2 represents the Cr (VI), that 4 parts of glucose xanthate react with the dichromate. In the Tables 3 and 4, the relationship between cellulose xanthogenate and Cr (VI), with related weights of 10.4 for Cr (VI).Table 3 Features of sample of ECx with Cr (VI).Full size tableTable 4 Researcher of process of the desorption.Full size tableMass balance in treatmentAdsorption is the phenomenon through which the removal of Cr (VI) achieved in the treatment systems; this quantified by means of the general balance equation of the treatment system as shown in Fig. 3.Adsorption is the phenomenon through which the removal of Cr (VI) achieved in treatment systems, this quantified by mass balance. Equation (1) shows the general balance of matter in the treatment system, together with the accumulation, inputs, and outputs of the system and the chemical process of adsorption.$${text{Acumulation }}upvarepsilon *frac{{partial {text{Cr}}left( {{text{VI}}} right)}}{{partial {text{t}}}} = {text{In}} frac{{partial {text{Cr}}left( {{text{VI}}} right)_{0} }}{{partial {text{t}}}} – {text{Out}}frac{{partial {text{Cr}}left( {{text{VI}}} right)}}{{partial {text{t}}}} – {text{Adsortion}},{rho b}frac{{partial {text{q}}}}{{partial {text{t}}}}$$
    (1)
    Accumulation represents by Eq. (1), where ∂C(VI) is the contaminant input to the treatment system, (ε) is the porosity of the bed, which calculated as the ratio between the density of the bed of treatment and the density of the microparticle of this biomass. This parameter must be above 0.548 achieved using particle diameters less than 0.212 mm, which favors contact between the contaminant and the particle49. The contaminant input to the treatment system represents by the design speed and the amount of contaminant that the system could treat. The output in the treatment system represents by the same input speed and the amount of contaminant that comes out. With these equations, the general material balance will be complete, summarized in Eq. (2), where it can see that the accumulation is equal to the input to the system, minus the output, and minus the adsorption.$$upvarepsilon *frac{{partial {text{Cr}}left( {{text{VI}}} right)}}{{partial {text{t}}}} = frac{{partial {text{Cr}} left( {{text{VI}}} right)}}{{partial {text{t}}}} – frac{{partial {text{Cr}} left( {{text{VI}}} right)}}{{partial {text{t}}}} – frac{{text{M}}}{{text{V}}}*frac{{partial {text{q}}}}{{partial {text{t}}}}$$
    (2)
    where V = System volume (ml), ε = Porosity, Co = Initial concentration of Cr (VI) (mg/ml), C = Final concentration Cr (VI) in the treated solution (mg/ml), Q = design flow (ml/min), Tb = Breaking time (Min), M = amount of biomass used (g), q = Adsorption capacity of the biomass used (mg/g).$${text{V}}*upvarepsilon *{text{Co}} = {text{Q}}*{text{Tb}}*{text{Co}} – {text{Q}}*{text{Tb}}*{text{C}} – {text{M}}*{text{q}}$$
    (3)
    Depending on the most important parameters when building a treatment system, Eq. (3) could use to model and validate the best form of treatment, for example, the necessary amount of biomass to use to treat a certain amount of contaminant, in the present investigation it used to establish the adsorption capacity in these initial treatment conditions. The remaining Eq. (4) determines the adsorption capacity.$${text{q}} = frac{{{text{QTbCo}}}}{{text{M}}} – frac{{{text{QTbCf}}}}{{text{M}}} – frac{{upvarepsilon {text{VCo}}}}{{text{M}}}$$
    (4)
    Adsorption capacity is generally taken through Eq. (5) for both batch and continuous experiments20,21But unlike Eqs. (5), (4) takes into account design variables such as flow rate (Q), rupture time (Tb), particle bed porosity ε, and vessel design volume (v).$${text{q}} = frac{{{text{v}}left( {{text{Co}} – {text{C}}} right)}}{{text{m }}}$$
    (5)
    where m: Mass used in the treatment, V: Volume, Co: Initial concentration, C: Final Concentration, Q: adsorption capacity.However, unlike Eqs. (5),  (4) considers the design variables such as flow rate (Q), rupture time (Tb), particle bed porosity ε and vessel design volume (v).When a desorption-elution process is involved for the reuse of biomass, Eq. (4) would be:$${text{q}}_{{text{T}}} = mathop sum limits_{j = 1}^{n} left[ {frac{{{text{QTbjCo}}}}{{text{M}}} – frac{{{text{QTbjCj}}}}{{text{M}}} – frac{{upvarepsilon {text{VCo}}}}{{text{M}}}} right]$$
    (6)
    where Q = design flow (ml/min), Tbj = Break time of use number j (Min), Co = Initial concentration of Cr (VI) (mg/ml), C = Final concentration Cr (VI) in the treated solution (mg/ml), V = System volume (ml), ε = Porosity, M = amount of biomass used (g), q_T = Total adsorption capacity of the biomass used (mg/g).This model (6) is design to determine the adsorption capacity when different elution processes have conducted, it will used to determine the new adsorption capacity and is one of the contributions of the present investigation.Result process of adsorptionsIn Fig. 4 shows the Cr (VI) adsorption process of the system.Figure 4Percentages of Cr (VI) removal the system for ECx.Full size imageVarious researchers have extensively studied the influence of factors such as bed height, flow rate and metal inlet concentration on rupture (Tb) curves. For example, the influence and similarity of the initial contaminant concentrations should be reflected as in the case of a tannery, with initial concentrations of 600 mg/l. Figure 4 shows the progress curves obtained for the study of Cr (VI) removal by the studied biomasses, reflecting the percentage of Cr (VI) removal in contrast to the treated volume, which is a very important parameter to time to scale the process.Regarding the effect of the input concentration, it can see in Fig. 5 that the breakpoint had a better performance in all the initial concentrations in the ECx biomass. comparing it with the EC-Na biomass (see Fig. 5), always obtaining breakpoints with more treated volume.Figure 5Percentages of Cr (VI) removal the system for EC-Na.Full size imageThe difference between the rupture curves between ECx and EC-Na indicates that the cellulose xanthate modification scheme should completed, although it can also elucidate that the EC-Na biomass has high yields compared to other biomass studied. for example, in Ref.34 investigate the biomass of E. crassipes without modifying, having removals below this alkaline cellulose.Adsorption capacitiesThrough Eq. (3), the adsorption capacity of ECx, using the initial concentration of 600 mg/l, since it was the maximum concentration used.The break point was around 1200 ml according to Fig. 6 and together with the flow rate of 15 ml/min; the break time obtained in 80 min.$${text{q}} = frac{{80{*}15{*}0.6}}{40} – frac{{80{*}15{*}0.04}}{40} – frac{{0.66{*}78{*}0.6}}{40}$$q: Adsorption capacity, Co: 0.6 mg/ml, C: 0.06 mg/ml, M: 40 g, Tb: rupture time 80 min, Q: 15 Flow rate ml/min, ε: 0.6649, V: Occupied volume: 70 ml.Figure 6Adsorption capacities in the different adsorption processes in the biomass ECx.Full size imageA result of 16 mg/g obtained in this continuous study for the biomass ECx. With this same equation it gives the capacity of the biomass EC-Na, with 11 mg/g.Desorption-Elution and reuseThrough Eq. (6), the sum of the Cr (VI) adsorption capacities established, after different biomass reuses due to EDTA elution. In the second treatment process, it yielded the following results under concentrations of 6 g/l of EDTA.$${text{q}}left( {text{T}} right) = frac{{60{*}15{*}0.6}}{40} – frac{{50{*}15{*}0.06}}{40} – frac{{0.66{*}68{*}0.6}}{40}$$Co: 0.6 mg/ml, C: 0.06 mg/ml, M: 45 g Biomass eluted with EDTA, Tb: rupture time: 60 min, Q: 15 Flow ml/min, ε: 0.6649, V: Occupied volume: 68 ml, q: 10 mg/g.Five Cr (VI) adsorption cycles performed using ECx and EC-Na cellulose in a continuous system to evaluate the regeneration and reuse potential. Between each biosorption cycle, a desorption cycle performed using three different concentrations of EDTA eluent.According to Figs. 6 and 7, although the adsorption capacity gradually decreases from the first adsorption process, it could consider that it is a satisfactory biomass recycling process and a design parameter for later stages of this treatment system.Figure 7Adsorption capacities in the different adsorption processes in the biomass EC-Na.Full size imageIn the experiments with concentrations of 6 g/l, five reuse processes obtained, obtaining a final sum of 52 mg/g. In concentrations of 3 g/l of EDTA, final capacities of 51 mg/g obtained lower than concentrations of 6 g/l but with half of this reagent. With concentrations of 1 g/l, final capacities of 33 mg/g obtained.The desorption processes of the EC-Na biomass with initial capacities of 11 mg/g were also evaluated and through desorption processes with EDTA of 3 g/l this biomass recycled on 5 occasions, reaching 32 mg/l in capacities of adsorption and like the EC-Na biomass, the ideal concentration in the process for desorption processes is 3 g/l, due to the considerable increase in reuse processes and low concentration compared to 6 g/l, which, although higher, does not this value is significant in the absorption capacity.Through Eq. (6) and with different bibliographic references, representative data obtained to feed this equation, determining the capacities of each of these biomasses together with the new capacities determining the desorption power of the different eluents shown and summarized in Table 4.For the EDTA eluent and with Eq. (6), satisfactory results evidenced by removing Al (II), reaching almost 150% of its adsorption capacity, corroborating what presented in the present investigation, also the EDTA reagent obtained interesting yields to recycle the cassava biomass increasing up to 40 mg/g. In Ref.39 used the biomass of Phanera vahlii to remove Cr (VI) obtaining results of 30 mg/g and with NaOH they reached capacities in the reuse process of this biomass up to 62 mg/g, reaching almost double of its total capacity41, also used NaOH for desorption processes with green synthesized nanocrystalline chlorapatite biomass, achieving results of 75% more. The eluent HCl is also a good chemical agent to use in desorption processes since it reached more than 100% in the reuse of biochar alginate for Cr (VI) but not so significant with biomass A. barbadensis Miller to remove Ni (II) and in40 significant results were also obtained to remove Pb (II) with pine cone Shell biomass. With the chemical agent HNO3, interesting contaminant recycling processes obtained, since more than 100% of the adsorption capacity of the biomasses used in this process used1,45.Mathematical models of adsorptionIn general, the models presented R2 greater than 0.95 for the adjustment of all the advance curves, which indicates a good adherence to the data, the model that best describes the behavior of the ECx system was the phenomenological model Thomas, which presented all the R2 values above 0.99.This model could use for the extension of the Cr (VI) ion biosorption system using cellulose xanthogenate, in the literature it is possible to observe that this model often tends to better adapt to the experimental data of the adsorption systems that use cellulose for the absorption of toxic metals28,30,31.With qt values remarkably close to the experimental values of Eq. (4) designed and presented in this investigation, indicating the validity of this equation where it reflects the maximum capacity obtained. Table 5 shows the adsorption constant of the Thomas model (Kt), which corresponds to the adsorption rate of Cr (VI) in the biomass49 This value was 0.048 (ml/mg*min) reflecting the speed with which Cr (VI) is chemisorbed in the biomass of ECx, in the EC-Na cellulose there was a Thomas model speed of 0.039 (ml/ mg*min) evidencing a lower adsorption rate than ECx. In the adsorption of Cr (VI) by rice biomass, the Thomas constant is 0.1 (ml/mg*min)47,50 also in the adsorption of Cr (VI) by biomass. Nanocrystalline chlorapatite biomass obtained at the Thomas constant 0.013 (ml/mg*min)49.Table 5 Summary of the experiments obtained with material ECx.Full size tableIn the Table 6, it presents summary of the experiments obtained with material EC-Na.Table 6 Summary of the experiments obtained with material EC-Na.Full size tableThe Cr (VI) adsorption process in the EC-Na biomass represented through the Bohart equation, since the sorption rate is proportional to the biomass capacity, obtaining an adsorption rate of 0.85(ml/mg*min). Having an alkalized biomass represents this model due to the homogeneity of this adsorbent.Mathematical models in desorption processesThe continuous desorption process with its fit to the Thomas model for biomass ECx always shows the fit of this model with significance, because this type of model fits representatively to desorption processes with good performance32,51 It can also verify that with values of qt it is close to the experimental values of Eq. (6) designed and presented in this research, indicating the validity of this equation again, where it reflects the maximum capacity obtained.In the Table 7. Show Summary of the experiments obtained with material ECx in process of desorption’s.Table 7 Summary of the experiments obtained with material ECx in process of desorption’s.Full size tableIn the Table 8 the EC-Na biomass had a different behavior and in its second and third cycle it adjusted to the Yoon model and later to the Bohart model.Table 8 Summary of the experiments obtained with material EC-Na in process of desorption’s.Full size tableThis behavior is due to the alkalinization of the biomass and this process makes the biomass a little more unstable. The values of qt, although a resemblance evidenced, were not so representative due to the little adjustment that there was with respect to the Thomas model. More

  • in

    Different effects of pesticides on transcripts of the endocrine regulation and energy metabolism in honeybee foragers from different colonies

    Eilers, E. J., Kremen, C., Smith Greenleaf, S., Garber, A. K. & Klein, A. M. Contribution of pollinator-mediated crops to nutrients in the human food supply. PLoS ONE 6, 21363 (2011).ADS 

    Google Scholar 
    Williams, P. H. The dependence of crop pollination within the European Union on pollination by honey bees. Agric. Zool. Rev. 6, 229–257 (1994).
    Google Scholar 
    Burd, M. Bateman’s principle and plant reproduction: the role of pollen limitation in fruit and seed set. Bot. Rev. 60, 83–139 (1994).MathSciNet 

    Google Scholar 
    Aguilar, R., Ashworth, L., Galetto, L. & Aizen, M. A. Plant reproductive susceptibility to habitat fragmentation: Review and synthesis through a meta-analysis. Ecol. Lett. 9, 968–980 (2006).
    Google Scholar 
    Potts, S. G. et al. Declines of managed honey bees and beekeepers in Europe. J. Apic. Res. 49, 15–22 (2010).
    Google Scholar 
    van Engelsdorp, D., Hayes, J., Underwood, R. M. & Pettis, J. A survey of honey bee colony losses in the U.S., fall 2007 to spring 2008. PLoS ONE 3, e4071 (2008).ADS 

    Google Scholar 
    Aizen, M. A. & Harder, L. D. The global stock of domesticated honey bees is growing slower than agricultural demand for pollination. Curr. Biol. 19, 915–918 (2009).CAS 

    Google Scholar 
    Van Engelsdorp, D. et al. Colony collapse disorder: A descriptive study. PLoS ONE 4, e6481 (2009).ADS 

    Google Scholar 
    Genersch, E. American foulbrood in honeybees and its causative agent, Paenibacillus larvae. J. Invertebr. Pathol. 103(suppl 1), 10–19 (2010).
    Google Scholar 
    Graystock, P., Yates, K., Darvill, B., Goulson, D. & Hughes, W. O. H. Emerging dangers: Deadly effects of an emergent parasite in a new pollinator host. J. Invertebr. Pathol. 114, 114–119 (2013).
    Google Scholar 
    Insolia, L. et al. Honey bee colony loss linked to parasites, pesticides and extreme weather across the United States. Sci. Rep. 12(1), 20787. https://doi.org/10.1038/s41598-022-24946-4 (2022).Article 
    ADS 
    CAS 

    Google Scholar 
    Sanchez-Bayo, F. & Goka, K. Pesticide residues and bees: A risk assessment. PLoS ONE 9(4), e94482 (2014).ADS 

    Google Scholar 
    Bolognesi, C. & Merlo, F. D. Pesticides: Human health effects. Encyclop. Environ. Health 1, 438–453 (2011).
    Google Scholar 
    Mullin, C. A. et al. High levels of miticides and agrochemicals in North American apiaries: Implications for honey bee health. PLoS ONE 1, e9754 (2015).
    Google Scholar 
    Calatayud-Vernich, P., Calatayud, F., Simó, E. & Picó, Y. Pesticide residues in honey bees, pollen and beeswax: Assessing beehive exposure. Environ. Pollut. 241, 106–114. https://doi.org/10.1016/j.envpol.2018.05.062 (2018).Article 
    CAS 

    Google Scholar 
    Woodcock, B. A. et al. Impacts of neonicotinoid use on long-term population changes in wild bees in England. Nat. Commun. 2016(7), 12459 (2016).ADS 

    Google Scholar 
    Zhao, H. et al. Review on effects of some insecticides on honey bee health. Pestic. Biochem. Physiol. 188, 105219. https://doi.org/10.1016/j.pestbp.2022.105219 (2022).Article 
    CAS 

    Google Scholar 
    Ludicke, J. C. & Nieh, J. C. Thiamethoxam impairs honey bee visual learning, alters decision times, and increases abnormal behaviors. Ecotoxicol. Environ. Saf. 193, 110367 (2020).CAS 

    Google Scholar 
    Tison, L., Duer, A., Púčiková, V., Greggers, U. & Menzel, R. Detrimental effects of clothianidin on foraging and dance communication in honey bees. PLoS ONE 15(10), e0241134 (2020).CAS 

    Google Scholar 
    Fent, K., Schmid, M. & Christen, V. Global transcriptome analysis reveals relevant effects at environmental concentrations of cypermethrin in honey bees (Apis mellifera). Environ. Pollut. 259, 113715 (2020).CAS 

    Google Scholar 
    Christen, V., Krebs, J., Bünter, I. & Fent, K. Biopesticide spinosad induces transcriptional alterations in genes associated with energy production in honey bees (Apis mellifera) at sublethal concentrations. J. Hazard. Mater. 378, 120736 (2019).CAS 

    Google Scholar 
    Christen, V., Krebs, J. & Fent, K. Fungicides chlorothanolin, azoxystrobin and folpet induce transcriptional alterations in genes encoding enzymes involved in oxidative phosphorylation and metabolism in honey bees (Apis mellifera) at sublethal concentrations. J. Hazard. Mater. 377, 215–226 (2019).CAS 

    Google Scholar 
    Fent, K., Haltiner, T., Kunz, P. & Christen, V. Insecticides cause transcriptional alterations of endocrine related genes in the brain of honey bee foragers. Chemosphere 260, 127542 (2020).ADS 
    CAS 

    Google Scholar 
    Christen, V., Grossar, D., Charrière, J. D., Eyer, M. & Jeker, L. Correlation between increased homing flight duration and altered gene expression in the brain of honey bee foragers after acute oral exposure to thiacloprid and thiamethoxam. Insect Sci. 1, 1–15 (2021).
    Google Scholar 
    Mao, W., Schuler, M. A. & Berenbaum, M. R. Disruption of quercetin metabolism by fungicide affects energy production in honey bees (Apis mellifera). Proc. Natl. Acad. Sci. USA 114(10), 2538–2543 (2017).ADS 
    CAS 

    Google Scholar 
    Christen, V., Kunz, P. Y. & Fent, K. Endocrine disruption and chronic effects of plant protection products in bees: Can we better protect our pollinators?. Environ. Pollut. 243(Pt B), 1588–1601 (2018).CAS 

    Google Scholar 
    Testai, E., Buratti, F. & Di Consiglio, E. Chlorpyrifos Hayes’ Handbook of Pesticide Toxicology 1505–1526 (Academic Press, 2010).
    Google Scholar 
    Eastmond, D. & Balakrishnan, S. Genotoxicity of Pesticides Hayes’ Handbook of Pesticide Toxicology 357–380 (Academic Press, 2010).
    Google Scholar 
    Urlacher, E. et al. Measurements of chlorpyrifos levels in forager bees and comparison with levels that disrupt honey bee odor-mediated learning under laboratory conditions. J. Chem. Ecol. 42(2), 127–138 (2016).CAS 

    Google Scholar 
    Li, Z. et al. Effects of sublethal concentrations of chlorpyrifos on olfactory learning and memory performances in two bee species, Apis mellifera and Apis cerana. Sociobiology 64, 174 (2017).
    Google Scholar 
    DeGrandi-Hoffman, G., Chen, Y. & Simonds, R. The effects of pesticides on queen rearing and virus titers in honey bees (Apis mellifera L.). Insects 4, 71–89 (2013).
    Google Scholar 
    Cutler, G. C., Purdy, J., Giesy, J. P. & Solomon, K. R. Risk to pollinators from the use of chlorpyrifos in the United States. In Ecological Risk Assessment for Chlorpyrifos in Terrestrial and Aquatic Systems in the United States Reviews of Environmental Contamination and Toxicology (eds Giesy, J. & Solomon, K.) (Springer, 2014).
    Google Scholar 
    Christen, V. & Fent, K. Exposure of honey bees (Apis mellifera) to different classes of insecticides exhibit distinct molecular effect patterns at concentrations that mimic environmental contamination. Environ. Pollut. 226, 48–59 (2017).CAS 

    Google Scholar 
    Stevenson, J. H. The acute toxicity of unformulated pesticides to worker honey bees (Apis mellifera L.). Plant Pathol. 27, 38–40 (1978).CAS 

    Google Scholar 
    Bartlett, D. W. et al. The strobilurin fungicides. Pest. Manag. Sci. 58, 649–662 (2002).CAS 

    Google Scholar 
    Ostiguy, N. et al. Honey bee exposure to pesticides: A four-year nationwide study. Insects. 10, 13 (2019).
    Google Scholar 
    Inoue, L. V. B., Domingues, C. E. C., Gregorc, A., Silva-Zacarin, E. C. M. & Malaspina, O. Harmful effects of pyraclostrobin on the fat body and pericardial cells of foragers of africanized honey bee. Toxics 10, 530. https://doi.org/10.3390/toxics10090530 (2022).Article 
    CAS 

    Google Scholar 
    Nicodemo, D. et al. Mitochondrial respiratory inhibition promoted by pyraclostrobin in fungi is also observed in honey bees. Environ. Toxicol. Chem. 39, 1267–1272 (2020).CAS 

    Google Scholar 
    Domingues, C. E. C., Inoue, L. V. B., Silva-Zacarin, E. C. M. & Malaspina, O. Foragers of Africanized honeybee are more sensitive to fungicide pyraclostrobin than newly emerged bees. Environ. Pollut. 266, 115267 (2020).
    Google Scholar 
    Tadei, R. et al. Late effect of larval co-exposure to the insecticide clothianidin and fungicide pyraclostrobin in Africanized Apis mellifera. Sci. Rep 9, 3277 (2019).ADS 

    Google Scholar 
    Zioga, E., Kelly, R., White, B. & Stout, J. C. Plant protection product residues in plant pollen and nectar: A review of current knowledge. Environ. Res. 189, 109873 (2020).CAS 

    Google Scholar 
    Corona, M. et al. Vitellogenin, juvenile hormone, insulin signaling, and queen honey bee longevity. Proc. Natl. Acad. Sci. USA 104, 7128–7133 (2007).ADS 
    CAS 

    Google Scholar 
    Winston, M. L. The Biology of the Honey Bee (Harvard University Press, 1987).
    Google Scholar 
    Ueno, T., Nakaoka, T., Takeuchi, H. & Kubo, T. Differential gene expression in the hypopharyngeal glands of worker honeybees (Apis mellifera L.) associated with an age-dependent role change. Zool. Sci. 8, 557–563 (2009).
    Google Scholar 
    Kubo, T. et al. Change in the expression of hypopharyngealgland proteins of the worker honeybees (Apis mellifera L.) with age and/or role. J. Biochem. 119, 291–295 (1996).CAS 

    Google Scholar 
    Ohashi, K., Sawata, M., Takeuchi, H., Natori, S. & Kubo, T. Molecular cloning of cDNA and analysis of expression of the gene for alpha-glucosidase from the hypopharyngeal gland of the honeybee Apis mellifera L. Biochem. Biophys. Res. Commun. 221, 380–385 (1996).CAS 

    Google Scholar 
    Ohashi, K., Natori, S. & Kubo, T. Expression of amylase and glucose oxidase in the hypopharyngeal gland with an age dependent role change of the worker honeybee (Apis mellifera L.). Eur. J. Biochem. 265, 127–133 (1999).CAS 

    Google Scholar 
    Chanchao, C., Padoongsupalai, R. & Sangvanich, P. Expression and characterization of α-glucosidase III in the dwarf honeybee, Apis florea (Hymenoptera: Apoidea: Apidae). Insect Sci. 14(4), 283–293 (2007).CAS 

    Google Scholar 
    Corby-Harris, V. & Snyder, L. A. Measuring hypopharyngeal gland acinus size in honey bee (Apis mellifera) Workers. J. Vis. Exp. 139, 58261 (2018).
    Google Scholar 
    Yamada, T. & Yamada, K. Comparison of long-term changes in size and longevity of bee colonies in mid-west Japan and Maui with and without exposure to pesticide, cold winters, and mites. PeerJ 8, e9505 (2020).
    Google Scholar 
    Rinkevich, F. D. et al. Genetics, synergists, and age affect insecticide sensitivity of the honey bee, Apis mellifera. PLoS ONE 10(10), e0139841 (2015).
    Google Scholar 
    Weidenmüller, A. The control of nest climate in bumblebee (Bombus terrestris) colonies: Interindividual variability and self reinforcement in fanning response. Behav. Ecol. 15(1), 120–128 (2004).MathSciNet 

    Google Scholar 
    Flatt, T., Tu, M. P. & Tatar, M. Hormonal pleiotropy and the juvenile hormone regulation of Drosophila development and life history. BioEssays 27, 999–1010 (2005).CAS 

    Google Scholar 
    Wu, M. C., Chang, Y. W., Lu, K. H. & Yang, E. C. Gene expression changes in honey bees induced by sublethal imidacloprid exposure during the larval stage. Insect. Biochem. Mol. Biol. 88, 12–20 (2017).CAS 

    Google Scholar 
    Ament, S. A., Corona, M., Pollock, H. S. & Robinson, G. E. Insulin signaling is involved in the regulation of worker division of labor in honey bee colonies. Proc. Natl. Acad. Sci. USA 105, 4226–4231 (2008).ADS 
    CAS 

    Google Scholar 
    Nicodemo, D. et al. Fipronil and imidacloprid reduce honeybee mitochondrial activity. Environ. Toxicol. Chem. 33(9), 2070–2075 (2014).CAS 

    Google Scholar 
    Syromyatnikov, M. Y., Lopatin, A. V., Starkov, A. A. & Popov, V. N. Isolation and properties of flight muscle mitochondria of the bumblebee Bombus terrestris (L.). Biochemistry 78(8), 909–914 (2013).CAS 

    Google Scholar 
    Dayer, F. C. Coadaptation of colony design and worker performance in honeybees. In Diversity in the Genus Apis (ed. Smith, D. R.) 2133–2245 (Westview Press, 1991).
    Google Scholar 
    Simon-Delso, N., Amaral-Rogers, V. & Belzunces, L. P. Systemic insecticides (neonicotinoids and fipronil): Trends, uses, mode of action and metabolites. Environ. Sci. Pollut. Res. 22, 5–34 (2015).CAS 

    Google Scholar 
    Evans, J. D. et al. Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect. Mol. Biol. 5, 645–656 (2006).
    Google Scholar 
    Pankiw, T. & Page, R. E. Response thresholds to sucrose predict foraging division of labor in honeybees. Behav. Ecol. Sociobiol. 47, 265–267 (2000).
    Google Scholar  More

  • in

    Salp blooms drive strong increases in passive carbon export in the Southern Ocean

    Roemmich, D. et al. Unabated planetary warming and its ocean structure since 2006. Nat. Clim. Change 5, 240–245 (2015).Article 
    ADS 

    Google Scholar 
    Frölicher, T. L. et al. Dominance of the Southern Ocean in anthropogenic carbon and heat uptake in CMIP5 models. J. Clim. 28, 862–886 (2015).Article 
    ADS 

    Google Scholar 
    Buesseler, K. O. & Boyd, P. W. Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean. Limnol. Oceanogr. 54, 1210–1232 (2009).Article 
    ADS 
    CAS 

    Google Scholar 
    Arteaga, L., Haentjens, N., Boss, E., Johnson, K. S. & Sarmiento, J. L. Assessment of export efficiency equations in the Southern Ocean applied to satellite-based net primary production. J. Geophys. Res.-Oceans 123, 2945–2964 (2018).Article 
    ADS 

    Google Scholar 
    Siegel, D. A. et al. Prediction of the export and fate of global ocean net primary production: the EXPORTS science plan. Front. Marine Sc. 3, 22 (2016).Perissinotto, R. & Pakhomov, E. A. The trophic role of the tunicate Salpa thompsoni in the Antarctic marine ecosystem. J. Mar. Syst. 17, 361–374 (1998).Article 

    Google Scholar 
    Atkinson, A., Siegel, V., Pakhomov, E. & Rothery, P. Long-term decline in krill stock and increase in salps within the Southern Ocean. Nature 432, 100–103 (2004).Article 
    ADS 
    CAS 

    Google Scholar 
    Perissinotto, R. & Pakhomov, E. A. Contribution of salps to carbon flux of marginal ice zone of the Lazarev sea, southern ocean. Mar. Biol. 131, 25–32 (1998).Article 
    CAS 

    Google Scholar 
    Phillips, B., Kremer, P. & Madin, L. P. Defecation by Salpa thompsoni and its contribution to vertical flux in the Southern Ocean. Mar. Biol. 156, 455–467 (2009).Article 

    Google Scholar 
    Stone, J. P. & Steinberg, D. K. Salp contributions to vertical carbon flux in the Sargasso Sea. Deep-Sea Res. Part I 113, 90–100 (2016).Article 
    CAS 

    Google Scholar 
    Ramaswamya, V., Sarin, M. M. & Rengarajan, R. Enhanced export of carbon by salps during the northeast monsoon period in the northern Arabian Sea. Deep-Sea Res. Part II 52, 1922–1929 (2005).Article 
    ADS 

    Google Scholar 
    Smith, K. L. et al. Large salp bloom export from the upper ocean and benthic community response in the abyssal northeast Pacific: day to week resolution. Limnol. Oceanogr. 59, 745–757 (2014).Article 
    ADS 
    CAS 

    Google Scholar 
    Madin, L. P. & Kremer, P. Determination of the filter feeding rates of salps (Tunicata, Thaliacea). ICES J. Mar. Sci. 52, 583–595 (1995).Article 

    Google Scholar 
    Wiebe, P. H., Madin, L. P., Haury, L. R., Harbison, G. R. & Philbin, L. M. Diel vertical migration by Salpa aspera and its potential for large-scale particulate organic matter transport to the deep-sea. Mar. Biol. 53, 249–255 (1979).Article 

    Google Scholar 
    Dadon-Pilosof, A., Lombard, F., Genin, A., Sutherland, K. R. & Yahel, G. Prey taxonomy rather than size determines salp diets. Limnol. Oceanogr. 64, 1996–2010 (2019).Article 
    ADS 

    Google Scholar 
    Stukel, M. R., Décima, M., Selph, K. E. & Gutiérrez-Rodríguez, A. Size-specific grazing and competitive interactions between large salps and protistan grazers. Limnol. Oceanogr. 66, 2521–2534 (2021).Madin, L. P. Production, composition and sedimentation of salp fecal pellets in oceanic waters. Mar. Biol. 67, 39–45 (1982).Article 

    Google Scholar 
    Michaels, A. F. & Silver, M. W. Primary production, sinking fluxes and the microbial food web. Deep-Sea Res. Part I 35, 473–490 (1988).Article 
    ADS 

    Google Scholar 
    Luo, J. Y. et al. Gelatinous zooplankton‐mediated carbon flows in the global oceans: a data‐driven modeling study. Glob. Biogeochem. Cycle 34, e2020GB006704 (2020).Kremer, P. & Madin, L. P. Particle retention efficiency of salps. J. Plankton Res. 14, 1009–1015 (1992).Article 

    Google Scholar 
    Harbison, G. R. & Gilmer, R. W. Feeding rates of pelagic tunicate Pegea confederata and 2 other salps. Limnol. Oceanogr. 21, 517–528 (1976).Article 
    ADS 
    CAS 

    Google Scholar 
    Harbison, G. R. & McAlister, V. L. The filter-feeding rates and particle retention efficiencies of 3 species of Cyclosalpa (Tunicata, Thaliacea). Limnol. Oceanogr. 24, 875–892 (1979).Article 
    ADS 

    Google Scholar 
    Mullin, M. M. In situ measurement of filtering rates of the salp Thalia democratica, on phytoplankton and bacteria. J. Plankton Res. 5, 279–288 (1983).Article 

    Google Scholar 
    Deibel, D. Clearance rates of the salp Thalia democratica fed naturally occurring particles. Mar. Biol. 86, 47–54 (1985).Article 

    Google Scholar 
    Fender, C. K. et al. Prey size spectra and predator:prey ratios of 7 species of New Zealand salps. Mar. Biol. (in press).Chiswell, S. M., Bostock, H. C., Sutton, P. J. H. & Williams, M. J. M. Physical oceanography of the deep seas around New Zealand: a review. N.Z. J. Mar. Freshw. Res. 49, 286–317 (2015).Article 

    Google Scholar 
    Henschke, N. et al. Salp-falls in the Tasman Sea: a major food input to deep-sea benthos. Mar. Ecol. Prog. Ser. 491, 165–175 (2013).Article 
    ADS 

    Google Scholar 
    Childerhouse, S., Dix, B. & Gales, N. Diet of new Zealand sea lions (Phocarctos hookeri) at the Auckland islands. Wildl. Res. 28, 291–298 (2001).Article 

    Google Scholar 
    Horn, P. L., Burrell, T., Connell, A. & Dunn, M. R. A comparison of the diets of silver (Seriolella punctata) and white (Seriolella caerulea) warehou. Mar. Biol. Res. 7, 576–591 (2011).Article 

    Google Scholar 
    Horn, P. L., Forman, J. S. & Dunn, M. R. Dietary partitioning by two sympatric fish species, red cod (Pseudophycis bachus) and sea perch (Helicolenus percoides), on Chatham Rise, New Zealand. Mar. Biol. Res. 8, 624–634 (2012).Article 

    Google Scholar 
    Forman, J. S., Horn, P. L. & Stevens, D. W. Diets of deepwater Oreos (Oreosomatidae) and orange roughy Hoplostethus atlanticus. J. Fish. Biol. 88, 2275–2302 (2016).Article 
    CAS 

    Google Scholar 
    Carroll, E. L. et al. Multi-locus DNA metabarcoding of zooplankton communities and scat reveal trophic interactions of a generalist predator. Sci. Rep. 9, 1–14 (2019).Savoye, N. et al. 234Th sorption and export models in the water column: a review. Mar. Chem. 100, 234–249 (2006).Article 
    CAS 

    Google Scholar 
    Sutton, P. J. H. The Southland Current: a subantarctic current. N.Z. J. Mar. Freshw. Res. 37, 645–652 (2003).Article 

    Google Scholar 
    Foxton, P. The distribution and life history of Salpa thompsoni Foxton with observations on a related species, Salpa gerlachei Foxton. Discovery Rep. 34, 1–116 (1966).Loeb, V. J. & Santora, J. A. Population dynamics of Salpa thompsoni near the Antarctic Peninsula: growth rates and interannual variations in reproductive activity (1993–2009). Prog. Oceanogr. 96, 93–107 (2012).Article 
    ADS 

    Google Scholar 
    Pakhomov, E. A. & Hunt, B. P. V. Trans-Atlantic variability in ecology of the pelagic tunicate Salpa thompsoni near the Antarctic Polar Front. Deep-Sea Res. Part II 138, 126–140 (2017).Article 

    Google Scholar 
    Lüskow, F., Pakhomov, E. A., Stukel, M. R. & Décima, M. Biology of Salpa thompsoni at the Chatham Rise, New Zealand: demography, growth, and diel vertical migration. Mar. Biol. 167, 1–18 (2020).Pakhomov, E. A. & Froneman, P. W. Zooplankton dynamics in the eastern Atlantic sector of the Southern Ocean during the austral summer 1997/1998—Part 2: grazing impact. Deep-Sea Res. Part II 51, 2617–2631 (2004).Article 
    ADS 

    Google Scholar 
    Iversen, M. H. et al. Sinkers or floaters? Contribution from salp pellets to the export flux during a large bloom event in the Southern. Ocean. Deep-Sea Res. Part II 138, 116–125 (2017).Article 
    CAS 

    Google Scholar 
    Buesseler, K. O., Boyd, P. W., Black, E. E. & Siegel, D. A. Metrics that matter for assessing the ocean biological carbon pump. Proc. Natl Acad. Sci. USA 117, 9679 (2020).Article 
    ADS 
    CAS 

    Google Scholar 
    Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).Article 

    Google Scholar 
    Hall, J., Safi, K. & Cumming, A. Role of microzooplankton grazers in the subtropical and subantarctic waters to the east of New Zealand. N.Z. J. Mar. Freshw. Res. 38, 91–101 (2004).Article 

    Google Scholar 
    Zeldis, J. R. & Décima, M. Mesozooplankton connect the microbial food web to higher trophic levels and vertical export in the New Zealand Subtropical Convergence Zone. Deep-Sea Res. Part I 155, 103146 (2020).Article 
    CAS 

    Google Scholar 
    Hall, J. A., James, M. R. & Bradford-Grieve, J. M. Structure and dynamics of the pelagic microbial food web of the Subtropical Convergence region east of New Zealand. Aquat. Micro. Ecol. 20, 95–105 (1999).Article 

    Google Scholar 
    Bradford-Grieve, J. M. et al. Pelagic ecosystem structure and functioning in the subtropical front region east of New Zealand in austral winter and spring 1993. J. Plankton Res. 21, 405–428 (1999).Article 

    Google Scholar 
    Nodder, S. & Gall, M. Pigment fluxes from the Subtropical Convergence region, east of New Zealand: relationships to planktonic community structure. N.Z. J. Mar. Freshw. Res. 32, 441–465 (1998).Article 
    CAS 

    Google Scholar 
    Nodder, S. D., Chiswell, S. M. & Northcote, L. C. Annual cycles of deep-ocean biogeochemical export fluxes in subtropical and subantarctic waters, southwest Pacific Ocean. J. Geophys. Res.: Oceans 121, 2405–2424 (2016).Article 
    ADS 
    CAS 

    Google Scholar 
    Kiko, R. et al. Biological and physical influences on marine snowfall at the equator. Nat. Geosci. 10, 852-+ (2017).Article 
    ADS 
    CAS 

    Google Scholar 
    Kelly, R. P., Shelton, A. O. & Gallego, R. Understanding PCR processes to draw meaningful conclusions from environmental DNA studies. Sci. Rep. 9, 12133 (2019).Article 
    ADS 

    Google Scholar 
    Caron, D. A., Madin, L. P. & Cole, J. J. Composition and degradation of salp fecal pellets: implications for vertical flux in oceanic environments. J. Mar. Res. 47, 829–850 (1989).Article 
    CAS 

    Google Scholar 
    Sempere, R., Yoro, S., Wambeke, F. V. & Charriere, B. Microbial decomposition of large organic particles in the northwestern Mediterranean Sea: an experimental approach. Mar. Ecol. Prog. Ser. 198, 61–72 (2000).Article 
    ADS 

    Google Scholar 
    Dell’Anno, A. & Corinaldesi, C. Degradation and turnover of extracellular DNA in marine sediments: ecological and methodological considerations. Appl. Environ. Microbiol. 70, 4384–4386 (2004).Article 
    ADS 

    Google Scholar 
    Torti, A., Lever, M. A. & Jørgensen, B. B. Origin, dynamics, and implications of extracellular DNA pools in marine sediments. Mar. Genomics 24, 185–196 (2015).Article 

    Google Scholar 
    Norris, R. Sediments of the Chatham Rise. N.Z. Dep. Sci. Ind. Res. Res. Bull. 159, 38 (1964).
    Google Scholar 
    Waite, A. M., Safi, K. A., Hall, J. A. & Nodder, S. D. Mass sedimentation of picoplankton embedded in organic aggregates. Limnol. Oceanogr. 45, 87–97 (2000).Article 
    ADS 

    Google Scholar 
    Gafar, N. A. & Schulz, K. G. A three-dimensional niche comparison of Emiliania huxleyi and Gephyrocapsa oceanica: reconciling observations with projections. Biogeosciences 15, 3541–3560 (2018).Article 
    ADS 
    CAS 

    Google Scholar 
    Eynaud, F., Giraudeau, J., Pichon, J. J. & Pudsey, C. J. Sea-surface distribution of coccolithophores, diatoms, silicoflagellates and dinoflagellates in the South Atlantic Ocean during the late austral summer 1995. Deep-Sea Res. Part I 46, 451–482 (1999).Article 

    Google Scholar 
    Hagino, K., Okada, H. & Matsuoka, H. Coccolithophore assemblages and morphotypes of Emiliania huxleyi in the boundary zone between the cold Oyashio and warm Kuroshio currents off the coast of Japan. Mar. Micropaleontol. 55, 19–47 (2005).Article 
    ADS 

    Google Scholar 
    Rhodes, L. L., Peake, B. M., Mackenzie, A. L. & Marwick, S. Coccolithophores Gephyrocapsa oceanica and Emiliana Huxleyi (Prymnesiophyceae = Haptophyceae) in New Zealand’s coastal waters: characteristics of blooms and growth in laboratory culture. N.Z. J. Mar. Freshw. Res. 29, 345–357 (1995).Article 

    Google Scholar 
    Ziveri, P., de Bernardi, B., Baumann, K.-H., Stoll, H. M. & Mortyn, P. G. Sinking of coccolith carbonate and potential contribution to organic carbon ballasting in the deep ocean. Deep-Sea Res. Part II 54, 659–675 (2007).Article 
    ADS 

    Google Scholar 
    Buesseler, K. O. et al. VERTIGO (VERtical Transport In the Global Ocean): a study of particle sources and flux attenuation in the North Pacific. Deep Sea Res. II 55, 1522–1539 (2008).Article 
    ADS 

    Google Scholar 
    Billett, D. S. M., Lampitt, R. S., Rice, A. L. & Mantoura, R. F. C. Seasonal sedimentation of phytoplankton to the deep-sea benthos. Nature 302, 520–522 (1983).Article 
    ADS 
    CAS 

    Google Scholar 
    Martin, J. H., Fitzwater, S. E., Gordon, R. M., Hunter, C. N. & Tanner, S. J. Iron, primary production and carbon nitrogen fluxes during the JGOFS North Atlantic Bloom Experiment. Deep-Sea Res. Part II 40, 115–134 (1993).Article 
    ADS 
    CAS 

    Google Scholar 
    Buesseler, K. O. et al. The effect of marginal ice-edge dynamics on production and export in the Southern Ocean along 170 degrees W. Deep-Sea Res. Part II 50, 579–603 (2003).Article 
    ADS 
    CAS 

    Google Scholar 
    Kiko, R. et al. Zooplankton-mediated fluxes in the eastern tropical north. Atl. Front. Mar. Sci. 7, 21 (2020).
    Google Scholar 
    Kelly, T. B. et al. The importance of mesozooplankton diel vertical migration for sustaining a mesopelagic food web. Front. Mar. Sci. 6, 508 (2019).Maiti, K., Charette, M. A., Buesseler, K. O. & Kahru, M. An inverse relationship between production and export efficiency in the Southern Ocean. Geophys. Res. Lett. 40, 1557–1561 (2013).Article 
    ADS 

    Google Scholar 
    Loeb, V. et al. Effects of sea-ice extent and krill or salp dominance on the Antarctic food web. Nature 387, 897–900 (1997).Article 
    ADS 
    CAS 

    Google Scholar 
    Steinberg, D. K. et al. Long-term (1993–2013) changes in macrozooplankton off the Western Antarctic Peninsula. Deep-Sea Res. Part I 101, 54–70 (2015).Article 

    Google Scholar 
    Cabanes, D. J. E. et al. First evaluation of the role of salp fecal pellets on iron biogeochemistry. Front. Mar. Sci. 3, 10 (2017).Article 

    Google Scholar 
    Belcher, A. et al. Krill faecal pellets drive hidden pulses of particulate organic carbon in the marginal ice zone. Nat. Commun. 10, 1–8 (2019).Manno, C. et al. Continuous moulting by Antarctic krill drives major pulses of carbon export in the north Scotia Sea, Southern Ocean. Nat. Commun. 11, 6051 (2020).Article 
    ADS 
    CAS 

    Google Scholar 
    Law, C. S. et al. Did dilution limit the phytoplankton response to iron addition in HNLCLSi sub-Antarctic waters during the SAGE experiment? Deep-Sea Res. Part II 58, 786–799 (2011).Article 
    ADS 
    CAS 

    Google Scholar 
    Gutiérrez‐Rodríguez, A. et al. Decoupling between phytoplankton growth and microzooplankton grazing enhances productivity in Subantarctic waters on Campbell Plateau, southeast of New Zealand. J. Geophys. Res.: Oceans 125, e2019JC015550 (2020).Sherman, J., Gorbunov, M. Y., Schofield, O. & Falkowski, P. G. Photosynthetic energy conversion efficiency in the West Antarctic Peninsula. Limnol. Oceanogr. 65, 1–14 (2020).Peterson, B. J. Aquatic primary productivity and the 14C-CO2 method: a history of the productivity problem. Annu Rev. Ecol. Syst. 11, 359–385 (1980).Article 

    Google Scholar 
    Landry, M. R. & Hassett, R. P. Estimating the grazing impact of marine microzooplankton. Mar. Biol. 67, 283–288 (1982).Article 

    Google Scholar 
    Landry, M. R., Haas, L. W. & Fagerness, V. L. Dynamics of microbial plankton communities—experiments in Kaneohe Bay, Hawaii. Mar. Ecol. Prog. Ser. 16, 127–133 (1984).Article 
    ADS 
    CAS 

    Google Scholar 
    Gutierrez-Rodriguez, A., Latasa, M., Estrada, M., Vidal, M. & Marrase, C. Carbon fluxes through major phytoplankton groups during the spring bloom and post-bloom in the Northwestern Mediterranean Sea. Deep Sea Res. Part I 57, 486–500 (2010).Article 
    CAS 

    Google Scholar 
    Gutierrez-Rodriguez, A. & Latasa, M. Pigment-based measurements of phytoplankton rates. in Phytoplankton Pigments Characterization, Chemotaxonomy and Applications in Oceanography (eds Roy, S. et al.) (Cambridge University Press, 2011) 472–495.Lorenzen, C. J. Determination of chlorophyll and pheo-pigments: spectrophotometric equations. Limnol. Oceanogr. 12, 343–346 (1967).Article 
    ADS 
    CAS 

    Google Scholar 
    Conover, R. J., Durvasula, R., Roy, S. & Wang, R. Probable loss of chlorophyll-derived pigments during passage through the gut of zooplankton, and some of the consequences. Limnol. Oceanogr. 31, 878–887 (1986).Article 
    ADS 
    CAS 

    Google Scholar 
    Latasa, M. A simple method to increase sensitivity for RP-HPLC phytoplankton pigment analysis. Limnol. Oceanogr. Meth 12, 46–53 (2014).Article 

    Google Scholar 
    Caporaso, J. G., Paszkiewicz, K., Field, D., Knight, R. & Gilbert, J. A. The Western English Channel contains a persistent microbial seed bank. ISME J. 6, 1089–1093 (2012).Article 
    CAS 

    Google Scholar 
    Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581-+ (2016).Article 
    CAS 

    Google Scholar 
    Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).Article 
    CAS 

    Google Scholar 
    McMurdie, P. J. & Holmes, S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217 (2013).Oksanen, J. et al. vegan: community ecology package. R package version 2.5-6. (2019).Piredda, R. et al. Diversity and temporal patterns of planktonic protist assemblages at a Mediterranean Long Term Ecological Research site. FEMS Microbiol. Ecol. 93, fiw200 (2017).Guillou, L. et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D597–D604 (2013).Article 
    CAS 

    Google Scholar 
    Pike, S. M., Buesseler, K. O., Andrews, J. & Savoye, N. Quantification of Th-234 recovery in small volume sea water samples by inductively coupled plasma-mass spectrometry. J. Radioanal. Nucl. Chem. 263, 355–360 (2005).Article 
    CAS 

    Google Scholar 
    Benitez-Nelson, C. R. et al. Testing a new small-volume technique for determining Th-234 in seawater. J. Radioanal. Nucl. Chem. 248, 795–799 (2001).Article 
    CAS 

    Google Scholar 
    Bone, Q. The Biology of Pelagic Tunicates (Oxford University Press, 1998).Foxton, P. An aid to the detailed examination of salps [Tunicata: Salpidae]. J. Mar. Biol. Assoc. UK 45, 679–681 (1965).Article 

    Google Scholar 
    Thompson, H. Pelagic Tunicates of Australia (Commonwealth Council for Scientific and Industrial Research, 1948).Iguchi, N. & Ikeda, T. Metabolism and elemental composition of aggregate and solitary forms of Salpa thompsoni (Tunicata: Thaliacea) in waters off the Antarctic Peninsula during austral summer 1999. J. Plankton Res. 26, 1025–1037 (2004).Article 
    CAS 

    Google Scholar 
    von Harbou, L. et al. Salps in the Lazarev Sea, Southern Ocean: I. feeding dynamics. Mar. Biol. 158, 2009–2026 (2011).Article 

    Google Scholar 
    Pakhomov, E. A., Dubischar, C. D., Strass, V., Brichta, M. & Bathmann, U. V. The tunicate Salpa thompsoni ecology in the Southern Ocean. I. Distribution, biomass, demography and feeding ecophysiology. Mar. Biol. 149, 609–623 (2006).Article 

    Google Scholar 
    Huntley, M. E., Sykes, P. F. & Marin, V. Biometry and trophodynamics of Salp thompsoni Foxton (Tunicata, Thaliacea) near the Antarctic peninsula in austral summer 1983–1984. Polar Biol. 10, 59–70 (1989).Article 

    Google Scholar 
    Knauer, G. A., Martin, J. H. & Bruland, K. W. Fluxes of particulate carbon, nitrogen, and phosphorus in the upper water column of the northeast Pacific. Deep Sea Res. Part I 26, 97–108 (1979).Article 
    ADS 
    CAS 

    Google Scholar 
    Karl, D. M. et al. Seasonal and interannual variability in primary production and particle flux at Station ALOHA. Deep-Sea Res. Part II 43, 539–568 (1996).Article 
    ADS 
    CAS 

    Google Scholar  More

  • in

    Analysis of Himalayan marmot distribution and plague risk in Qinghai province of China using the “3S” technology

    Himalayan marmot habitat analysisThe environmental factors including temperature, vegetation and elevation are the key drivers for the wildlife in alpine ecosystems32. Specific landform attributes such as slope and elevation and vegetation cover affect the population and burrowing of rodents33. For example, rodent burrows in the Western Usambara Mountains in Tanzania were only found at an elevation of above 1600 m33. However, the Himalayan marmot seems to prefer to inhabit areas with low elevation and high land surface temperature34. In this study, the data showed that 76.25% of the Himalayan marmots were found in areas with elevation values of 3400–4600 m. The majority of marmots were found in areas with slopes of 5–20° and vegetation cover higher than 60%. Most marmots were found in alpine meadows, a few were found in temperate grasslands and alpine grasslands, and none were found in other grassland types.Preliminary statistical analysis of vegetation cover, grass type, vegetation type, and Himalayan marmot distribution sample sites obtained using spatial geographic information technology revealed that the meadow grassland areas with lush grass growth, more dominant plants, and abundant food had more marmots. When the vegetation cover reached 0.60–1.00, the number of marmot distribution sample sites was the highest. Dense grass is an ideal habitat and provides concealment for Himalayan marmots, and the abundant plant types provide sufficient food for marmots. In contrast, no marmots were distributed in the alpine scrub, coniferous forest, and alpine snow/ice covered areas where vegetation growth was poor, vegetation cover was low, and food was relatively scarce. Moreover, 70.24% of Himalayan marmots were found in alpine meadows with a wide variety of plant species, including Poaceae, Cyperaceae, and grasses. This finding indicated that alpine meadows are more suitable for Himalayan marmots and have more advantageous habitat conditions compared with other grassland types. The elevation of alpine meadows is 3236–5126 m, and the vegetation is mainly meadows with simple vegetation structure, substantial vegetation cover and dense vegetation growth, and a wide variety of plants, rich food, soft grass, and good palatability. Therefore, alpine meadows provide good natural habitats and foraging sites for marmots.Habitat selection of large rodents is influenced by a combination of vegetation cover availability, food availability, and population density35. Vegetation cover is an important parameter that describes vegetation communities and ecosystems and is closely related to vegetation quantity and productivity. The quality of habitat vegetation is an important factor that affects the spatial distribution of plateau rodents. Both feeding and concealment depend on vegetation, and the height and cover of edible plants and vegetation suitable for concealment determine the choice of vegetation type by marmots. Thus, vegetation cover becomes an important factor for habitat selection by marmots. Different grassland types determine different plant conditions, and selection of different vegetation conditions can increase the chances of survival and improve the reproductive success of marmots; therefore, grassland type is an important ecological factor in habitat selection by marmots. A study showed that the ecological factors affecting habitat selection of Himalayan marmots are mainly topography, anthropogenic disturbance, and vegetation8. Another study concluded that habitat selection by Himalayan marmots is closely related to elements such as topography, landform, temperature, precipitation, and vegetation24.The functions of burrows’ physical parameters is to protect the Himalayan marmots from natural enemies and bad weather36. There is clearly influence of slope on habitat selection by marmots. When the slope is large, wind is strong, and burrows are not well hidden; this makes them difficult to defend against enemies, unsafe for survival, and not conducive to hibernation during winter. In addition, Himalayan marmots prefer to burrow on sunny aspect, because the temperature is suitable and the vegetation is lush, which is suitable for marmots to breed. Therefore, the number of marmot burrows gradually decreases with increasing slope and ubac. Although flat and low-lying areas with small slopes are good for marmots to create dens, rainwater will easily flow into the dens during summer rainfall, which will kill marmots. Therefore, a suitable slope and sunny aspect are also very important for habitat selection by marmots.Application of the predictive spatial distribution map of Himalayan marmots in Qinghai provincePlague surveillance is the main measure used for plague prevention and control in China. Although we have made many improvements in plague surveillance, the traditional method of dragnet surveillance still consumes a lot of human and material resources, is inefficient. The pasture area of Qinghai province is approximately 380,000 km2, and the identified natural plague focus is approximately 180,000 km2; therefore, there is still 200,000 km2 of pasture where the distribution of Himalayan marmots and plague have not been identified. Currently, RS technology is widely used in the fields of mapping and ecological surveillance18,19,21,22,37.Applications of RS technology in areas such as malaria, dengue, schistosomiasis and plague have been previously reported27,37. Using GIS combined with remotely sensed data, Proches Hieronimo et al. found that the presence of small mammals was positively influenced by elevation, whereas the presence of fleas was clearly influenced by land management features, and thus these observations have positive implications for plague surveillance27. In this study, RS technology combined with field validations were used to determine the distribution and areas of different types of grasslands in Qinghai province, and the average density of Himalayan marmot distribution in different types of grasslands. The high-, low-, and very low-density areas of Himalayan marmot distribution were identified. The soil map, vegetation map, administrative map, and marmot density statistics were merged to form the spatial data and attribute data basis for the information system to map the distribution of Himalayan marmot and determine the area of Himalayan marmot distribution. Generally speaking, the occurrence of human plague epidemic is closely related to the local animal plague epidemic2. However, a large part of the high-density distribution of Himalayan marmots is located in uninhabited areas and the areas are generally sparsely populated, which also indicates that we should reasonably allocate plague prevention and control resources to areas where human plague is most likely to occur to prevent the occurrence of human plague epidemics.Field validation for verificationThrough field validation and information from local farmers and herdsmen, we confirmed that Himalayan marmots inhabited 68 sample sites in Tongde, Zeku, Guinan, Xunhua, Haiyan, Ulan, Qilian, Hualong, and Huzhu counties. Among them, Tongde, Zeku, Guinan, Xunhua, Haiyan, Ulan, and Qilian counties have all historically experienced marmot plague outbreaks and can be considered as reliable natural plague foci38. The data from this field validation are consistent with the previous survey data and the epidemic history of the counties in Qinghai province39.MAE can better reflect the actual number of errors in prediction values; the smaller the MAE value, the higher the prediction accuracy. The MAE derived from the field validation data was 0.1331 and the prediction accuracy was 0.8669. The accuracy of the predicted Himalayan marmot spatial distribution reached 87%, which indicated that the predicted probability map of the Himalayan marmot spatial distribution can better predict the potential marmot distribution.The predicted spatial distribution map of Himalayan marmot in Qinghai province was then compared with environmental information such as elevation, vegetation, grass type, slope, and aspect of 352 field survey sites. The obtained RS data showed that the prediction results were excellent, and the predicted spatial distribution map of Himalayan marmot in Qinghai province was drawn with high accuracy. The prediction map visually reflects the different density distribution of Himalayan marmots; this allows us to optimize the settings and reasonable spatial layout of animal plague surveillance sites and improve surveillance efficiency.Application of marmot information collection system V3.0Marmot information collection system V3.0 was developed based on the “3S” technology standardizing the collection of surveillance data, and makes the management and analysis of information more convenient and faster. This study revolutionized the traditional method of considering plague-stricken counties as the plague foci, and effectively reduces the work intensity of operators and improves the data collection efficiency. In 2016 and 2017, we applied this system to the animal plague surveillance tasks in the plague-stricken counties of Haidong, Hainan, and Haibei in Qinghai province, and standardized the collection of provincial geographic location data of animal plague surveillance (data not shown). In 2018, we also applied this system in Wulan County, which frequently experiences plague, and achieved a good application effect (data not shown).In the next step, we will expand the pilot areas (mainly national and provincial plague surveillance sites), collect surveillance data from each surveillance site, continuously optimize and update the system, improve the efficiency of data analysis and utilization, detect the plague epidemic in marmot in a timely and accurate manner, correctly determine the epidemic trend of plague in marmots, and attempt to strictly prevent the plague from spreading to humans. We plan to use a new model of drone surveillance to create a multidimensional, three-dimensional, real-time big data plague surveillance information reporting system to enhance early plague warnings and prediction in Qinghai province and even in the country, which will be of positive practical significance to serve and guarantee the Belt and Road Initiative. These approaches are expected to provide new technical means for plague investigation and research, and to provide references for setting up plague surveillance programs and prediction for the natural Himalayan marmot plague focus in Qinghai province and the QTP. More

  • in

    Environmental factors driving the abundance of Philaenus spumarius in mesomediterranean habitats of Corsica (France)

    Saponari, M. et al. Infectivity and transmission of Xylella fastidiosa by Philaenus spumarius (Hemiptera: Aphrophoridae) in Apulia, Italy. J. Econ. Entomol. 107, 1316–1319. https://doi.org/10.1603/EC14142 (2014).Article 

    Google Scholar 
    Cornara, D., Bosco, D. & Fereres, A. Philaenus spumarius: When an old acquaintance becomes a new threat to European agriculture. J. Pest Sci. 91, 957–972. https://doi.org/10.1007/s10340-018-0966-0 (2018).Article 

    Google Scholar 
    Weaver, C. R. & King, D. Meadow spittlebug, Philaenus leucophthalmus (L.). Ohio Agric. Exp. Stn. Res. Bull. 741, 258 (1954).
    Google Scholar 
    Halkka, A., Halkka, L., Halkka, O., Roukka, K. & Pokki, J. Lagged effects of North Atlantic Oscillation on spittlebug Philaenus spumarius (Homoptera) abundance and survival. Glob. Change Biol. 12, 2250–2262. https://doi.org/10.1111/j.1365-2486.2006.01266.x (2006).Article 

    Google Scholar 
    Cruaud, A. et al. Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: A case study in Corsica. Sci. Rep. 8, 15628. https://doi.org/10.1038/s41598-018-33957-z (2018).Article 
    CAS 

    Google Scholar 
    Godefroid, M. et al. Climate tolerances of Philaenus spumarius should be considered in risk assessment of disease outbreaks related to Xylella fastidiosa. J. Pest Sci. 2021, 1–14. https://doi.org/10.1007/s10340-021-01413-z (2021).Article 

    Google Scholar 
    Farigoule, P. et al. Vectors as sentinels: Rising temperatures increase the risk of Xylella fastidiosa outbreaks. Biology 11, 1299. https://doi.org/10.3390/biology11091299 (2022).Article 

    Google Scholar 
    Drosopoulos, S. & Asche, M. Biosystematic studies on the spittlebug genus Philaenus with the description of a new species. Zool. J. Linn. Soc. 101, 169–177. https://doi.org/10.1111/j.1096-3642.1991.tb00891.x (1991).Article 

    Google Scholar 
    Godefroid, M. & Durán, J. M. Composition of landscape impacts the distribution of the main vectors of Xylella fastidiosa in southern Spain. J. Appl. Entomol. 146, 666–675. https://doi.org/10.1111/jen.13003 (2022).Article 

    Google Scholar 
    Karban, R. & Strauss, S. Y. Physiological tolerance, climate change, and a northward range shift in the spittlebug, Philaenus spumarius. Ecol. Entomol. 29, 251–254. https://doi.org/10.1111/j.1365-2311.2004.00576.x (2004).Article 

    Google Scholar 
    Chmiel, S. M. & Wilson, M. C. Estimation of the lower and upper developmental threshold temperatures and duration of the nymphal stages of the meadow Spittlebug, Philaenus spumarius. Environ. Entomol. 8, 682–685. https://doi.org/10.1093/ee/8.4.682 (1979).Article 

    Google Scholar 
    Yurtsever, S. On the polymorphic meadow spittlebug, Philaenus spumarius (L.) (Homoptera: Cercopidae). Turk. J. Zool. 24, 447–460 (2000).
    Google Scholar 
    Ahmed, D. D. & Davidson, R. H. Life history of the meadow spittlebug in Ohio. J. Econ. Entomol. 43, 905–908. https://doi.org/10.1093/jee/43.6.905 (1950).Article 

    Google Scholar 
    Whittaker, J. B. Cercopid spittle as a microhabitat. Oikos 21, 59–64. https://doi.org/10.2307/3543839 (1970).Article 

    Google Scholar 
    Drosopoulos, S. New data on the nature and origin of colour polymorphism in the spittlebug genus Philaenus (Hemiptera: Aphorophoridae). Ann. Soc. Entomol. Fr. NS 39, 31–42. https://doi.org/10.1080/00379271.2003.10697360 (2003).Article 

    Google Scholar 
    Bodino, N. et al. Phenology, seasonal abundance, and host-plant association of spittlebugs (Hemiptera: Aphrophoridae) in vineyards of Northwestern Italy. Insects 12, 1012. https://doi.org/10.3390/insects12111012 (2021).Article 

    Google Scholar 
    Cornara, D. et al. Natural areas as reservoir of candidate vectors of Xylella fastidiosa. Bull. Insectol. 74, 173–180 (2021).
    Google Scholar 
    Gargani, E. et al. A five-year survey in Tuscany (Italy) and detection of Xylella fastidiosa subspecies multiplex in potential insect vectors, collected in Monte Argentario. Redia 104, 75–88. https://doi.org/10.19263/REDIA-104.21.09 (2021).Article 

    Google Scholar 
    Morente, M. et al. Distribution and relative abundance of insect vectors of Xylella fastidiosa in olive groves of the iberian peninsula. Insects 9, 175. https://doi.org/10.3390/insects9040175 (2018).Article 

    Google Scholar 
    Delong, D. et al. Spittle-insect vectors of Pierce’s disease virus. I. Characters, distribution, and food plants. Hilgardia 19, 339–356 (1950).Article 

    Google Scholar 
    Bodino, N. et al. Phenology, seasonal abundance and stage-structure of spittlebug (Hemiptera: Aphrophoridae) populations in olive groves in Italy. Sci. Rep. 9, 1–17. https://doi.org/10.1038/s41598-019-54279-8 (2019).Article 
    CAS 

    Google Scholar 
    Wiegert, R. G. Population energetics of meadow spittlebugs (Philaenus spumarius L.) as affected by migration and habitat. Ecol. Monogr. 34, 217–241. https://doi.org/10.2307/1948501 (1964).Article 

    Google Scholar 
    Dongiovanni, C. et al. Plant selection and population trend of spittlebug immatures (Hemiptera: Aphrophoridae) in olive groves of the Apulia region of Italy. J. Econ. Entomol. 112, 67–74. https://doi.org/10.1093/jee/toy289 (2019).Article 

    Google Scholar 
    Bodino, N. et al. Spittlebugs of mediterranean olive groves: Host-plant exploitation throughout the year. Insects 11, 130. https://doi.org/10.3390/insects11020130 (2020).Article 

    Google Scholar 
    Villa, M., Rodrigues, I., Baptista, P., Fereres, A. & Pereira, J. A. Populations and host/non-host plants of spittlebugs nymphs in olive orchards from northeastern Portugal. Insects 11, 720. https://doi.org/10.3390/insects11100720 (2020).Article 

    Google Scholar 
    Antonatos, S. et al. Seasonal appearance, abundance, and host preference of Philaenus spumarius and Neophilaenus campestris (Hemiptera: Aphrophoridae) in olive groves in Greece. Environ. Entomol. 50, 1474–1482. https://doi.org/10.1093/ee/nvab093 (2021).Article 

    Google Scholar 
    Hasbroucq, S., Casarin, N., Ewelina, C., Bragard, C. & Grégoire, J.-C. Distribution, adult phenology and life history traits of potential insect vectors of Xylella fastidiosa in Belgium. Belg. J. Entomol. 92, 2569 (2020).
    Google Scholar 
    Mesmin, X. et al. Interaction networks between spittlebugs and vegetation types in and around olive and clementine groves of Corsica; implications for the spread of Xylella fastidiosa. Agric. Ecosyst. Environ. 334, 107979. https://doi.org/10.1016/j.agee.2022.107979 (2022).Article 

    Google Scholar 
    Albre, J., García-Carrasco, J. M. & Gibernau, M. Ecology of the meadow spittlebug Philaenus spumarius in the Ajaccio region (Corsica)—I: Spring. Bull. Entomol. Res. 111, 246–256. https://doi.org/10.1017/S0007485320000711 (2021).Article 

    Google Scholar 
    Andersson, P., Löfstedt, C. & Hambäck, P. A. Insect density–plant density relationships: A modified view of insect responses to resource concentrations. Oecologia 173, 1333–1344. https://doi.org/10.1007/s00442-013-2737-1 (2013).Article 

    Google Scholar 
    Hambäck, P. A., Inouye, B. D., Andersson, P. & Underwood, N. Effects of plant neighborhoods on plant–herbivore interactions: Resource dilution and associational effects. Ecology 95, 1370–1383. https://doi.org/10.1890/13-0793.1 (2014).Article 

    Google Scholar 
    Otway, S. J., Hector, A. & Lawton, J. H. Resource dilution effects on specialist insect herbivores in a grassland biodiversity experiment. J. Anim. Ecol. 74, 234–240 (2005).Article 

    Google Scholar 
    Lago, C. et al. Flight performance and the factors affecting the flight behaviour of Philaenus spumarius the main vector of Xylella fastidiosa in Europe. Sci. Rep. 11, 17608. https://doi.org/10.1038/s41598-021-96904-5 (2021).Article 
    CAS 

    Google Scholar 
    Casarin, N. et al. Investigating dispersal abilities of Aphrophoridae in European temperate regions to assess the threat of potential Xylella fastidiosa-based pathosystems. J. Pest Sci. https://doi.org/10.1007/s10340-022-01562-9 (2022).Article 

    Google Scholar 
    Bodino, N. et al. Dispersal of Philaenus spumarius (Hemiptera: Aphrophoridae), a vector of Xylella fastidiosa, in olive grove and meadow agroecosystems. Environ. Entomol. 50, 267–279. https://doi.org/10.1093/ee/nvaa140 (2020).Article 
    CAS 

    Google Scholar 
    Santoiemma, G., Tamburini, G., Sanna, F., Mori, N. & Marini, L. Landscape composition predicts the distribution of Philaenus spumarius, vector of Xylella fastidiosa, in olive groves. J. Pest. Sci. 92, 1101–1109. https://doi.org/10.1007/s10340-019-01095-8 (2019).Article 

    Google Scholar 
    Cappellari, A. et al. Spatio-temporal dynamics of vectors of Xylella fastidiosa subsp. pauca across heterogeneous landscapes. Entomol. Gen. 42, 515–521. https://doi.org/10.1127/entomologia/2022/1427 (2022).Article 

    Google Scholar 
    Avosani, S., Tattoni, C., Mazzoni, V. & Ciolli, M. Occupancy and detection of agricultural threats: The case of Philaenus spumarius, European vector of Xylella fastidiosa. Agric. Ecosyst. Environ. 324, 107707. https://doi.org/10.1016/j.agee.2021.107707 (2022).Article 

    Google Scholar 
    Allier, C. & Lacoste, A. Processus dynamiques de reconstitution dans la série du Quercus ilex en Corse. In Vegetation Dynamics in Grasslans, Healthlands and Mediterranean Ligneous Formations 83–91 (Springer, 1981).Delbosc, P., Bioret, F. & Panaïotis, C. Plant landscape of Corsica: Typology and mapping plant landscape of Cap Corse region and Biguglia Pond (Springer Nature, 2020).Book 

    Google Scholar 
    Chessel, D., Dufour, A.-B. & Thioulouse, J. The ade4 package—I: One-table methods. R. News 4, 5–10 (2004).
    Google Scholar 
    Biedermann, R. & Niedringhaus, R. The Plant-and Leafhoppers of Germany: Identification Key to All Species (Wabv Fründ, 2009).
    Google Scholar 
    Stöckmann, M., Biedermann, R., Nickel, H. & Niedringhaus, R. The Nymphs of the Planthoppers and Leafhoppers of Germany (WABV, 2013).
    Google Scholar 
    INRAE-CBGP. Arthemis DB@se – ARTHropod Ecology, Molecular Identification and Systematics. https://arthemisdb.supagro.inrae.frhttps://doi.org/10.15454/TBGRIB. Accessed 2021.Xu, T. & Hutchinson, M. ANUCLIM version 6.1 user guide. Aust. Natl. Univ. Fenner Sch. Environ. Soc. Canberra 2011, 256 (2011).
    Google Scholar 
    Quintana-Seguí, P. et al. Analysis of near-surface atmospheric variables: Validation of the SAFRAN analysis over France. J. Appl. Meteorol. Climatol. 47, 92–107. https://doi.org/10.1175/2007JAMC1636.1 (2008).Article 

    Google Scholar 
    Hijmans, R. J., Phillips, S., Leathwick, J. & Elith, J. Dismo: Species Distribution Modeling https://CRAN.R-project.org/package=dismo (2017).Faraway, J. J. Extending the Linear Model with R: Generalized Linear, Mixed Effects and Nonparametric Regression Models (Chapman and Hall/CRC, 2006).MATH 

    Google Scholar 
    Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R. J. 9, 378–400. https://doi.org/10.3929/ethz-b-000240890 (2017).Article 

    Google Scholar 
    R Core Team. R: A Language and Environment for Statistical Computing https://www.R-project.org/ (2019).Hardin, J. W. & Hilbe, J. M. Generalized Linear Models and Extensions 4th edn. (Stata Press, 2018).MATH 

    Google Scholar 
    Harrison, X. A. Using observation-level random effects to model overdispersion in count data in ecology and evolution. PeerJ 2, e616. https://doi.org/10.7717/peerj.616 (2014).Article 

    Google Scholar 
    Bolker, B. M. et al. Generalized linear mixed models: A practical guide for ecology and evolution. Trends Ecol. Evol. 24, 127–135. https://doi.org/10.1016/j.tree.2008.10.008 (2009).Article 

    Google Scholar 
    Hartig, F. DHARMa: Residual Diagnostics for Hierarchical (Multi-Level/Mixed) Regression Models https://CRAN.R-project.org/package=DHARMa (2020).Lüdecke, D., Ben-Shachar, M. S., Patil, I., Waggoner, P. & Makowski, D. performance: An R package for assessment, comparison and testing of statistical models. J. Open Sourc. Softw. 6, 3139. https://doi.org/10.21105/joss.03139 (2021).Article 

    Google Scholar 
    Fox, J. & Weisberg, S. An {R} Companion to Applied Regression Third edn, https://socialsciences.mcmaster.ca/jfox/Books/Companion/ (Sage, Thousand Oaks CA, 2019).Lenth, R. V. Emmeans: Estimated Marginal Means, Aka Least-Squares Means https://CRAN.R-project.org/package=emmeans (2021).Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biom. J. 50, 346–363. https://doi.org/10.1002/bimj.200810425 (2008).Article 
    MATH 

    Google Scholar 
    Fernández-Mazuecos, M. & Vargas, P. Ecological rather than geographical isolation dominates Quaternary formation of Mediterranean Cistus species. Mol. Ecol. 19, 1381–1395. https://doi.org/10.1111/j.1365-294X.2010.04549.x (2010).Article 
    CAS 

    Google Scholar 
    Berenbaum, M. R. & Feeny, P. P. 1. Chemical mediation of host-plant specialization: The papilionid paradigm. In Specialization, Speciation, and Radiation (ed. Tilmon, K.) 3–19 (University of California Press, 2008). https://doi.org/10.1525/california/9780520251328.003.0001.Kapantaidaki, D. E., Antonatos, S., Evangelou, V., Papachristos, D. P. & Milonas, P. Genetic and endosymbiotic diversity of Greek populations of Philaenus spumarius, Philaenus signatus and Neophilaenus campestris, vectors of Xylella fastidiosa. Sci. Rep. 11, 3752. https://doi.org/10.1038/s41598-021-83109-z (2021).Article 
    CAS 

    Google Scholar 
    Mesmin, X. et al. Ooctonus vulgatus (Hymenoptera, Mymaridae), a potential biocontrol agent to reduce populations of Philaenus spumarius (Hemiptera, Aphrophoridae) the main vector of Xylella fastidiosa in Europe. PeerJ 8, e8591. https://doi.org/10.7717/peerj.8591 (2020).Article 

    Google Scholar 
    Denancé, N. et al. Several subspecies and sequence types are associated with the emergence of Xylella fastidiosa in natural settings in France. Plant Pathol. 66, 1054–1064. https://doi.org/10.1111/ppa.12695 (2017).Article 
    CAS 

    Google Scholar 
    EFSA, Delbianco, A., Gibin, D., Pasinato, L. & Morelli, M. Update of the Xylella spp host plant database—systematic literature search up to 31 December 2020. EFSA J. 19, 6. https://doi.org/10.2903/j.efsa.2021.6674 (2021).Article 

    Google Scholar 
    Soubeyrand, S. et al. Inferring pathogen dynamics from temporal count data: The emergence of Xylella fastidiosa in France is probably not recent. New Phytol. 219, 824–836. https://doi.org/10.1111/nph.15177 (2018).Article 

    Google Scholar 
    Roy, J. & Sonié, L. Germination and population dynamics of Cistus species in relation to fire. J. Appl. Ecol. 29, 647–655. https://doi.org/10.2307/2404472 (1992).Article 

    Google Scholar 
    Whittaker, J. B. Density regulation in a population of Philaenus spumarius (L.) (Homoptera: Cercopidae). J. Anim. Ecol. 42, 163–172. https://doi.org/10.2307/3410 (1973).Article 

    Google Scholar 
    Chapman, D. et al. Improving knowledge of Xylella fastidiosa vector ecology: modelling vector occurrence and abundance in the wider landscape in Scotland. Project Final Report. PHC2020/04, Scotland’s Centre of Expertise for Plant Health (PHC) https://doi.org/10.5281/zenodo.6523478 (2022).Saponari, M., Giampetruzzi, A., Loconsole, G., Boscia, D. & Saldarelli, P. Xylella fastidiosa in olive in Apulia: Where we stand. Phytopathology 109, 175–186. https://doi.org/10.1094/PHYTO-08-18-0319-FI (2019).Article 
    CAS 

    Google Scholar 
    López-Mercadal, J. et al. Collection of data and information in Balearic Islands on biology of vectors and potential vectors of Xylella fastidiosa (GP/EFSA/ALPHA/017/01). EFSA Supp. Publ. 18, 10. https://doi.org/10.2903/sp.efsa.2021.EN-6925 (2021).Article 

    Google Scholar  More

  • in

    Short-term sedimentation dynamics in mesotidal marshes

    No plants were collected or harmed during this study, and all research involving plants followed relevant national, and international guidelines and legislation.Study areaThe study site encloses a wetland area bordering Ramalhete Channel, in the western part of the Ria Formosa lagoon, a mesotidal system located in southern Portugal (Fig. 1). Lunar tides are semi-diurnal, with a mean tidal range of about 2 m that can reach up to 3.5 m during spring tides. Offshore waves have no major propagation inside the lagoon33,34. Water circulation inside the lagoon is mostly driven by tides. The lagoon extends over 55 km along the coast and is connected to the ocean through six tidal inlets35. The three westmost inlets of the system (Ancão, Faro-Olhão, and Armona), which together capture ca. 90% of the total prism, are highly interconnected, with a strong residual circulation from Faro-Olhão Inlet directed towards Ancão and Armona inlets (located in Fig. 1), during both spring and neap tides36. The tidal currents in Ramalhete Channel, connecting the Faro-Olhão and Ancão Inlet, have high tidal asymmetry and shifts in tidal dominance, from flood to ebb. There are no significant fluvial inputs into the lagoon, with a yearly average terrestrial sediment influx of around 2 × 105 m3/yr37, reaching the system through small streams. The main sediment delivery to the system is through the inlets, though there are few studies assessing related fluxes. The net sediment entry through the stabilized Faro-Olhão Inlet is estimated at 1.4 × 105 m3/year38. Recent sedimentation rates in the marsh of the westmost edge of the lagoon were estimated at 1.1 ± 0.1 mm/yr39.The lagoon system is composed of large salt marsh patches, tidal flats and a complex net of natural, and partially dredged tidal channels. The tidal flats (vegetated and non-vegetated) and salt marshes represent more than 2/3 of the total lagoon area. The salt marshes comprise silt and fine sand40, while coarser (sand to shingle) shell-rich sediment, of marine provenance, is found on tidal channels and the lower domain of intertidal flats41. The dominant intertidal species are Spartina maritima and the seagrass Zostera noltei, the latter occupying an estimated area of 1304 ha, which represent 45% of the total intertidal area42.Figure 1Location of the field site in the Ria Formosa lagoon western sector over a satellite image collected in 2019 (South Portugal; upper panel); zoom to monitoring stations S1 to S4 (left lower panel); and field view of the studied site (right lower panel). Map generated with ArcGIS 10.8 (http://www.esri.com) and Adobe Illustrator 2022. Map data: Google Earth 7.3, image Landsat / Copernicus.Full size imageExperimental setup and data analysisAn experimental setup was deployed in the study area to assess dominant local topography, hydrodynamics (water levels and current velocities), Suspended Sediment Concentrations (SSCs), Deposition Rates (DRs), vegetation characteristics, and bed sediment grain size and organic matter content. Measurements were made during a full tide cycle, on a spring tide (tidal range = 3.2 m), and on a neap tide (tidal range = 1.8 m). Sampling was conducted in four wetland stations: S1 and S2 in a vegetated tidal flat comprising Zostera noltei; S3 in the low marsh comprising Spartina maritima; and S4 in the mid-upper marsh with the most abundant species of Sarcocornia perennis and Atriplex portucaloides (see S1 to S4, Fig. 1); the tidal flat is interrupted by a small oblique secondary tidal creek that flows near S2 station.Stations of sediment sampling and equipment deployment along the transect are illustrated in Fig. 2. During neap tide there was no data collection in S4, since the inundation time of the station was very short. The profile elevation was measured using Real Time Kinematic Differential Global Positioning System (RTK-DGPS, Trimble R6; vertical error in the order of few centimetres), and the slope of each habitat within a transect was calculated and expressed in percentage (%). Vegetation at each point was characterized by the canopy height, calculated as the average shoot length.Suspended Sediment Samplers (SSSs) were installed during low tide in the monitored stations using siphon samplers (Fig. 2) and recovered in the next low tide. These samplers consist of 0.5 L bottles with two holes on the cap, one for water intake and the other for air exhaust, according to the method described in13. Each intake tube is adjusted to form a siphon (i.e., inverse U), allowing to control the water level at which intake starts. Siphons were aligned at the same elevation along the transect for spring and neap tides, which means that all SSSs were collecting at the same time within the tidal cycle. During spring tide, in S1 and S2 at the tidal flat, SSSs were sampling at 0.1, 0.9, and 1.2 m from the bed, while at S3 SSSs were sampling at 0.7 and 1.0 m from the bed, and at S4 the SSS was sampling at 0.1 m from the bed (Fig. 2). During neap tide, in S1 and S2, SSSs were sampling at 0.1 and 0.9 m from the bed, while at S3 the SSS was sampling at 0.7 m from the bed.Surficial sediment samples were collected in each habitat to characterize the sediment grain size (d50) and content of organic matter (% OM). Sediment traps were installed in 3 replicates, during low tide, at each sampling point to measure the short-term sediment deposition rate (i.e., deposition over a tidal cycle, following procedures of43). Traps consisted of 3 cm diameter pre-labeled cylindrical tubes (Falcon® tubes, 50 ml). Traps and sediment samples were transported to the laboratory and maintained in a fridge. The sediment content was washed, and both the inorganic and organic weights were determined.The measured inorganic DR (g/m2/hr) was calculated as:$${text{DR}} = {raise0.7exhbox{${{text{W}}_{{{text{DS}}}} }$} !mathord{left/ {vphantom {{{text{W}}_{{{text{DS}}}} } {{text{A}} cdot {text{T}}}}}right.kern-0pt} !lower0.7exhbox{${{text{A}} cdot {text{T}}}$}}$$
    (1)
    where WDS is the weight of deposited sediment (in grams), A is the area of the sediment trap opening (m2), and T is in hours. Two different tide durations were considered to compute DRs, one assuming T equal to the hydroperiod in each station, and one assuming T equal to the entire tide duration (~ 12.4 h). These measured DRs are hereon mentioned as flood and tide DRs (DRflood and DRtide, respectively). The former is an expression of the actual deposition rate within the flood phase, during the period in which each station is inundated (and therefore active deposition can take place). The latter is the value used to compare with DRs in literature, which typically corresponds to values averaged over multiple tidal cycles (thus accounting for the entire tide duration).Tide levels were measured in the field using pressure sensors (PT, InSitu Inc. Level TROLL; ~ 2 cm from the bed), deployed from S2 towards S4 (Fig. 2). Velocity currents were measured at 20 cm from the bed, using an electromagnetic current meter (EMCM; Infinity Series JFE Advantech Co., Ltd; in S2 to S4; Fig. 2), and raw data (recording interval: 30 s) were filtered using a 10 min moving average for cross-shore and longshore components. To identify tidal asymmetry and assess the related phase dominance, tidal current skewness was calculated through the formula described in44 by which:
    $$Sk_{U} = frac{{frac{1}{N – 1}mathop sum nolimits_{t = 1}^{N} left( {U_{t} – overline{U}} right)^{3} }}{{left( {frac{1}{N – 1}mathop sum nolimits_{t = 1}^{N} left( {U_{t} – overline{U}} right)^{2} } right)^{{{raise0.7exhbox{$3$} !mathord{left/ {vphantom {3 2}}right.kern-0pt} !lower0.7exhbox{$2$}}}} }}$$
    (2)
    where N is the number of recordings, Ut is the input velocity signal and (overline{U}) is the mean velocity. Positive/negative skewness indicates flood/ebb dominance (assuming that flood currents are positive).Figure 2Deployment of the sediment traps, SSSs and devices (electromagnetic current meter EMCM; pressure transducer PT) in the stations (S1 to S4) during spring tide (sketch is exaggerated in the vertical).Full size imageComplementary to the measured DRs, theoretical DRs were also determined from the data, allowing us to link the sediment and flow data collected, and validate the deposition patterns from the traps. The theoretical deposition rate was determined based on45 formula:$${text{DR}} = left{ {begin{array}{*{20}c} {{text{C}}_{{text{b}}} cdot {text{w}}_{{text{s}}} cdot left( {1 – frac{{{uptau }_{{text{b}}} }}{{{uptau }_{{{text{cd}}}} }}} right)} & {{uptau }_{{text{b}}} < {uptau }_{{{text{cd}}}} } \ 0 & {{uptau }_{{text{b}}} ge {uptau }_{{{text{cd}}}} } \ end{array} } right.$$ (3) where Cb is the SSC at the bed, ws is the flock settling velocity, τb is the bed shear stress and τcd is the corresponding critical value for deposition.To determine the settling rate of the flocculates, the modified Stokes’ velocity for cohesive sediment was used, taking shape factors α and β (α = β = 1 for perfectly spherical particles):$${text{w}}_{{text{s}}} = frac{{upalpha }}{{upbeta }} cdot frac{{left( {{uprho }_{{text{s}}} - {uprho }_{{text{w}}} } right) cdot {text{g}} cdot {text{D}}_{50}^{2} }}{{{uprho }_{{text{w}}} cdot 18 cdot {upnu }}}$$ (4) where ρw and ρs are the densities of the water and sediment, respectively and ν is the kinematic viscosity of water (~ 106 m2/s).The bed shear stress τb was calculated from the measured current magnitude, |U| using the law of the wall:$$begin{array}{*{20}c} \ {{uptau }_{{text{b}}} = {uprho }_{{text{w}}} cdot {text{u}}_{*}^{2} , {text{u}}_{*} = frac{left| U right| cdot kappa }{{ln left( {{raise0.7exhbox{$z$} !mathord{left/ {vphantom {z {z_{0} }}}right.kern-0pt} !lower0.7exhbox{${z_{0} }$}}} right)}} } \ end{array} { }$$ (5) where κ is the von Kármán constant (~ 0.4) and z0 is the roughness length. For Zostera noltei, the roughness length was estimated at 5 mm46, value that was also used in the other stations, in lack of related estimate for marsh plants.The critical shear for deposition, τcd, was calculated using the formula47:$$sqrt {frac{{{uptau }_{{{text{cd}}}} }}{{{uprho }_{{text{w}}} }}} = left{ {begin{array}{*{20}c} {0.008} & {{text{w}}_{{text{s}}} le 5 cdot 10^{ - 5} {text{m}}/{text{s}}} \ {0.094 + 0.02 cdot {text{log}}_{10} left( {{text{w}}_{{text{s}}} } right)} & {3 cdot 10^{ - 4} le {text{w}}_{{text{s}}} le 5 cdot 10^{ - 5} {text{m}}/{text{s}}} \ {0.023} & {{text{w}}_{{text{s}}} ge 3 cdot 10^{ - 4} {text{m}}/{text{s}}} \ end{array} } right.$$ (6) Theoretical values of minimum SSCs needed for these DRs were also calculated, assuming that there is constant deposition (i.e., setting τb = 0), and compared with the field results. More

  • in

    The temperature dependence of microbial community respiration is amplified by changes in species interactions

    Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M. & Charnov, E. L. Effects of size and temperature on metabolic rate. Science 293, 2248–2251 (2001).Article 
    CAS 

    Google Scholar 
    Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173 (2006).Article 
    CAS 

    Google Scholar 
    Lopez-Urrutia, A., San Martin, E., Harris, R. P. & Irigoien, X. Scaling the metabolic balance of the oceans. Proc. Natl Acad. Sci. USA 103, 8739–8744 (2006).Article 
    CAS 

    Google Scholar 
    Yvon-Durocher, G. et al. Reconciling the temperature dependence of respiration across timescales and ecosystem types. Nature 487, 472–476 (2012).Article 
    CAS 

    Google Scholar 
    Crowther, T. W. et al. Quantifying global soil carbon losses in response to warming. Nature 540, 104–108 (2016).Article 
    CAS 

    Google Scholar 
    Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on Earth. Proc. Natl Acad. Sci. USA 115, 6506–6511 (2018).Article 
    CAS 

    Google Scholar 
    Rivkin, R. B. & Legendre, L. Biogenic carbon cycling in the upper ocean: effects of microbial respiration. Science 291, 2398–2400 (2001).Article 
    CAS 

    Google Scholar 
    Friedlingstein, P. et al. Uncertainties in CMIP5 climate projections due to carbon cycle feedbacks. J. Clim. 27, 511–526 (2014).Article 

    Google Scholar 
    Smith, T. P. et al. Community-level respiration of prokaryotic microbes may rise with global warming. Nat. Commun. 10, 5124 (2019).Article 

    Google Scholar 
    Antwis, R. E. et al. Fifty important research questions in microbial ecology. FEMS Microbiol. Ecol. 93, fix044 (2017).Bardgett, R. D., Freeman, C. & Ostle, N. J. Microbial contributions to climate change through carbon cycle feedbacks. ISME J. 2, 805–814 (2008).Article 
    CAS 

    Google Scholar 
    Enquist, B. J. et al. Scaling from traits to ecosystems: developing a general trait driver theory via integrating trait-based and metabolic scaling theories. Adv. Ecol. Res. 52, 249–318 (2015).Article 

    Google Scholar 
    Allen, A. P., Gillooly, J. F. & Brown, J. H. Linking the global carbon cycle to individual metabolism. Funct. Ecol. 19, 202–213 (2005).Article 

    Google Scholar 
    Schramski, J. R., Dell, A. I., Grady, J. M., Sibly, R. M. & Brown, J. H. Metabolic theory predicts whole-ecosystem properties. Proc. Natl Acad. Sci. USA 112, 2617–2622 (2015).Article 
    CAS 

    Google Scholar 
    Alster, C. J., Koyama, A., Johnson, N. G., Wallenstein, M. D. & von Fischer, J. C. Temperature sensitivity of soil microbial communities: an application of macromolecular rate theory to microbial respiration. J. Geophys. Res. Biogeosci. 121, 1420–1433 (2016).Article 

    Google Scholar 
    Yvon-Durocher, G. et al. Five years of experimental warming increases the biodiversity and productivity of phytoplankton. PLoS Biol. 13, e1002324 (2015).Article 

    Google Scholar 
    Garzke, J., Connor, S. J., Sommer, U. & O’Connor, M. I. Trophic interactions modify the temperature dependence of community biomass and ecosystem function. PLoS Biol. 17, e2006806 (2019).Foster, K. R. & Bell, T. Competition, not cooperation, dominates interactions among culturable microbial species. Curr. Biol. 22, 1845–1850 (2012).Article 
    CAS 

    Google Scholar 
    Coyte, K. Z., Schluter, J. & Foster, K. R. The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666 (2015).Article 
    CAS 

    Google Scholar 
    Machado, D. et al. Polarization of microbial communities between competitive and cooperative metabolism. Nat. Ecol. Evol. 5, 195–203 (2021).Article 

    Google Scholar 
    Bradford, M. A. et al. Cross-biome patterns in soil microbial respiration predictable from evolutionary theory on thermal adaptation. Nat. Ecol. Evol. 3, 223–231 (2019).Article 

    Google Scholar 
    Garcia-Martin, E. E., McNeill, S., Serret, P. & Leakey, R. J. G. Plankton metabolism and bacterial growth efficiency in offshore waters along a latitudinal transect between the UK and Svalbard. Deep Sea Res. I 92, 141–151 (2014).Article 
    CAS 

    Google Scholar 
    Davidson, E. A., Richardson, A. D., Savage, K. E. & Hollinger, D. Y. A distinct seasonal pattern of the ratio of soil respiration to total ecosystem respiration in a spruce-dominated forest. Glob. Change Biol. 12, 230–239 (2006).Article 

    Google Scholar 
    Dutkiewicz, S., Follows, M. J. & Bragg, J. G. Modeling the coupling of ocean ecology and biogeochemistry. Glob. Biogeochem. Cycles 23, GB4017 (2009).Article 

    Google Scholar 
    Follows, M. J., Dutkiewicz, S., Ward, B. & Follett, C. in Microbial Ecology of the Oceans 3rd edn (eds Gasol, J. & Kirchman, D.) Ch. 12 (John Wiley, 2018).Letten, A. D. & Stouffer, D. B. The mechanistic basis for higher-order interactions and non-additivity in competitive communities. Ecol. Lett. 22, 423–436 (2019).Article 

    Google Scholar 
    Grilli, J., Barabás, G., Michalska-Smith, M. J. & Allesina, S. Higher-order interactions stabilize dynamics in competitive network models. Nature 548, 210–213 (2017).Article 
    CAS 

    Google Scholar 
    Maynard, D. S., Crowther, T. W. & Bradford, M. A. Competitive network determines the direction of the diversity–function relationship. Proc. Natl Acad. Sci. USA 114, 11464–11469 (2017).Article 
    CAS 

    Google Scholar 
    Fiegna, F., Moreno-Letelier, A., Bell, T. & Barraclough, T. G. Evolution of species interactions determines microbial community productivity in new environments. ISME J. 9, 1235–1245 (2015).Article 

    Google Scholar 
    Lawrence, D. et al. Species interactions alter evolutionary responses to a novel environment. PLoS Biol. 10, e1001330 (2012).Article 
    CAS 

    Google Scholar 
    Harcombe, W. R., Chacón, J. M., Adamowicz, E. M., Chubiz, L. M. & Marx, C. J. Evolution of bidirectional costly mutualism from byproduct consumption. Proc. Natl Acad. Sci. USA 115, 12000–12004 (2018).Article 
    CAS 

    Google Scholar 
    Goldford, J. E. et al. Emergent simplicity in microbial community assembly. Science 361, 469–474 (2018).Article 
    CAS 

    Google Scholar 
    Yvon-Durocher, G. et al. Methane fluxes show consistent temperature dependence across microbial to ecosystem scales. Nature 507, 488–491 (2014).Article 
    CAS 

    Google Scholar 
    Fox, J. W. & Harpole, W. S. Revealing how species loss affects ecosystem function: the trait-based price equation partition. Ecology 89, 269–279 (2008).Article 

    Google Scholar 
    Kontopoulos, D., Smith, T. P., Barraclough, T. G. & Pawar, S. Adaptive evolution shapes the present-day distribution of the thermal sensitivity of population growth rate. PLoS Biol. 18, e3000894 (2020).Article 
    CAS 

    Google Scholar 
    Wilson, W. G. & Lundberg, P. Biodiversity and the Lotka–Volterra theory of species interactions: open systems and the distribution of logarithmic densities. Proc. R. Soc. Lond. B 271, 1977–1984 (2004).Article 

    Google Scholar 
    Rossberg, A. G. in Food Webs and Biodiversity 181–191 (John Wiley & Sons, 2013).Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).Article 
    CAS 

    Google Scholar 
    Garcia, F. C., Bestion, E., Warfield, R. & Yvon-Durocher, G. Changes in temperature alter the relationship between biodiversity and ecosystem functioning. Proc. Natl Acad. Sci. USA 115, 10989–10999 (2018).Article 
    CAS 

    Google Scholar 
    Padfield, D., O’Sullivan, H. & Pawar, S. rTPC and nls.multstart: a new pipeline to fit thermal performance curves in R. Methods Ecol. Evol. 12, 1138–1143 (2021).Article 

    Google Scholar  More

  • in

    Pathogen evasion of social immunity

    Ant hostWe used workers of the invasive Argentine ant, Linepithema humile, as host species. As typical for invasive ants, populations of this species lack territorial structuring and instead consist of interconnected nests forming a single supercolony with constant exchange of individuals between nests40. We collected L. humile queens, workers and brood in 2011, 2016 and 2022 from its main supercolony in Europe that extends more than 6,000 km along the coasts of Portugal, Spain and France40,41,42, from a field population close to Sant Feliu de Guíxols, Spain (41° 49’ N, 3° 03’ E). Field-collected ants were reared in large stock colonies in the laboratory. For the experiments, we sampled worker ants from outside the brood chambers and placed them into petri dishes with plastered ground (Alabastergips, Boesner, BAG), subjected to their respective treatments as detailed below. Experiments were carried out in a temperature- and humidity-controlled room at 23 °C, 65% relative humidity and a 12 h day/night light cycle. During experiments, ants were provided with ad libitum access to a sucrose-water solution (100 g l−1) and plaster was watered every 2–3 d to keep humidity high.Collection of this unprotected species from the field was in compliance with international regulations, such as the Convention on Biological Diversity and the Nagoya Protocol on Access and Benefit-Sharing (ABS, permit numbers ABSCH-IRCC-ES-260624-1 ESNC126 and SF0171/22). All experimental work followed European and Austrian law and institutional ethical guidelines.Fungal pathogensAs pathogen, we used the obligate-killing entomopathogenic fungus Metarhizium, whose infectious conidiospores naturally infect ants43,44,45 by penetrating their cuticles, killing them and growing out to produce highly infectious sporulating carcasses23,46. We used a total of six strains of the two species M. robertsii and M. brunneum, all isolated from the soil of the same natural population—an agricultural field at the Research Centre Årslev, Denmark27,47, which makes host co-infections with these sympatric strains in the field likely. As in ref. 24, we used three strains of M. robertsii (R1: KVL 12-36, R2: KVL 12-38, R3: KVL 12-35) and three of M. brunneum (B1: KVL 13-13, B2: KVL 12-37, B3: KVL 13-14), all obtained from the University of Copenhagen, Denmark (B. M. Steinwender, J. Eilenberg and N. V. Meyling).We started our selection experiment by exposing the ants to a mix of the six strains in equal proportions. To this end, each strain was grown separately from monospore cultivates from its long-term storage (43% glycerol (Sigma-Aldrich, G2025) in skimmed milk, −80 °C) on SDA plates (Sabouraud-4% dextrose agar, Sigma-Aldrich, 84088-500G) at 23 °C until sporulation. Conidiospores (abbreviated to ‘spores’) were collected by suspending them in sterile 0.05% Triton X-100 (Sigma-Aldrich, X-100; in milliQ water, autoclaved) and mixed in equal amounts to a total concentration of 1 × 106 spores ml−1. Before mixing, we confirmed that all strains had ≥98% germination.We exposed worker ants individually to the fungal pathogen by dipping them into the spore suspension using clean forceps (Gebrüder Martin; bioform, B32d). Afterwards, each ant was brieftly placed on filter paper (Whatman; VWR, 512-1027) to remove excess liquid before being placed into its experimental Petri dish.Serial passage experimentWe tested for the long-term effect of social immunity on pathogen selection, in which the pathogen was serially cycled through the host in the absence or presence of social immunity while the host population remained constant.Experimental design and procedureAfter exposure to the fungal spore mix, worker ants were either kept alone (individual host treatment, n = 10 replicate lines) or together with two untreated nestmates (social host treatment, n = 10 replicate lines; Fig. 1a). Individual ants could only protect themselves by individual immunity (selfgrooming behaviour and their physiological immune system), while the attended ants experienced both individual and social immunity due to the additional allogrooming by their caregiving nestmates. Thus, comparing the two host conditions revealed the effect of social immunity.As sanitary care by the nestmates reduces the pathogens’ success to kill the exposed individuals, we had to set up more experimental dishes of the social host treatment to obtain equal numbers of sporulating carcasses under both selection treatments, from which we then collected the spores for the next host infection cycle. For the individual treatment, we exposed an average of 23 workers per cycle, while an average of 40 workers per cycle were exposed in the social host treatment. The experiment was run for 10 host passages, that is, 27 weeks. In total, 6,312 workers (2,299 in the individual and 4,013 in the social host treatment) were exposed during the course of the experiment, and 8,026 nestmates were used. To obtain the spore suspensions for the next steps, we then collected and pooled the outgrowing spores of the first 8 carcasses produced per replicate line and cycle (that is, a total of n = 800 carcasses from the individual and n = 800 carcasses from the social host treatment, over the 10 host passages). Dead nestmates were not considered (see below).In detail, at each host cycle, the freshly exposed ants were placed into Petri dishes with plastered, humidified ground (Ø 3.5 cm for the individual and Ø 6 cm for the social host condition; both Bioswisstec AG, 10035 and 10060) in the absence (individual host treatment) or presence (social host treatment) of two untreated nestmates. We checked survival daily for 8 d. Ants that died within 24 h after exposure were excluded from the experiment as their mortality could not yet have resulted from infection, but rather from treatment procedures. Ants dying from days 2 to 8 were checked for internal Metarhizium infections by surface-sterilization (washing the carcass in 70% ethanol (Honeywell; Bartelt, 24194-2.5l; diluted with water) for a few seconds, rinsing it in distilled water, incubating in 3% bleach (Sigma-Aldrich, 1056142500) in sterile 0.05% Triton X-100 for 3 min and rinsing it again three times in water48), followed by incubation in a Petri dish on humidified filter paper at 23 °C until day 13, when they were checked for Metarhizium spore outgrowth. This timeline was chosen as preliminary work showed that the exposed ants die mostly on days 4 to 8 (median day 5, for both individual and social host treatments) after exposure and that sporulation required no longer than 5 d in our experimental conditions, so that a duration of 13 d per cycle also allowed for the later dying ants to complete sporulation. Preliminary work further revealed that in cases where nestmates contracted the disease, they died at a delayed timepoint and with spore outgrowth mostly around the mouthparts. These characteristics were used to distinguish between the directly exposed ants and infected nestmates in the experiment where ants were not colour-marked. The carcasses of sporulating nestmates were excluded from further procedures. An additional control experiment using 120 sham-treated ants showed no Metarhizium outgrowth, so that all Metarhizium outgrowth in our experiment could be attributed to our experimental infections. Carcasses with saprophytic outgrowth were not considered. For each host passage and each replicate line, we collected the spores of the first 8 ants dying after day 1 from their Metarhizium-sporulating carcasses at day 13 in 0.05% Triton X-100, pooled and counted them using an automated cell counter (Cellometer Auto M10, Nexcelom Bioscience). The concentration of each pool was then adjusted to 1 × 106 spores ml−1, and was used directly (that is, in the absence of any intermediate fungal growth step on agar plates) for exposing the ants in the next host infection cycle. The ants of each host passage were thus dipped in the same spore concentration. The remaining spore suspension was frozen at −80 °C in a long-term storage for further analysis.Pathogen diversity and strain compositionWe analysed which strains were present and in which proportion after 5 and 10 passages in each of the 10 individual and 10 social replicate lines. To this end, we first extracted total DNA from the respective spore pools (n = 40), which we analysed (1) quantitatively for the respective representation of M. robertsii vs M. brunneum (using species-specific real-time PCR targeting the PR1-gene sequence; detailed below) and (2) qualitatively for which of the 6 original strains were still present in the pool (using strain-specific microsatellite analysis; detailed below). We used this first estimate of remaining strain diversity and composition of each pool to determine how many spores we had to analyse separately for their strain identity after individualization by FACS sorting and growing them individually as colony forming units (c.f.u.s). This clone-level strain identification was again performed using microsatellite analysis (n = 1,347 individualized clones from the 40 spore mixes, in addition to n = 27 spores from the 6 ancestral strains; detailed below). Such clonal separation was needed since expansion of the spore mix by growth on SDA plates was not representative of the genetic composition of the strains in the pool, due to strong strain–strain growth inhibition when growing in a mix.In detail, we extracted the DNA of the 6 ancestral strains and the 40 spore mixes (10 each for individual and social lines at passages 5 and 10), as well as of 27 individualized clones of the ancestral strains and 1,374 clones from the 40 pools of passages 5 and 10, by centrifuging 100 µl of the spore suspensions in 1.5 ml tubes (Eppendorf, 0030120086) at full speed for 1 min and discarding the supernatant. Nuclease-free water (50 µl) was added and the spores were crushed in a bead mill (Qiagen TissueLyser II, 85300) at 30 Hz for 10 min using acid-washed glass beads (425–600 µm; Sigma-Aldrich, G8772). DNA was extracted using a DNeasy blood and tissue kit (Qiagen, 69506) following the manufacturer’s instructions, using a final elution volume of 50 µl buffer AE.For the quantitative species-level analysis of the pools, we performed quantitative real-time PCR (qPCR) using primers and differently labelled probes24 that we had developed on the basis of the sequence of the PR1 gene49 (forward: 5′ TCGATATTTTCGCTCCTG, reverse 5′-TTGTTAGAGCTGGTTCTGAAG, PR1 probe M. brunneum: 5′-(6-carboxyfluorescein (6FAM))TATTGTACCTACCTCGATAAGCTTAGAGAC(BHQ1), PR1 probe M. robertsii: 5′-(hexachloro-fluorescein (HEX))AGTATTGTACCTCGATAAGCTCGGAGAC(BHQ1)). Reactions were performed in 20 μl volumes using 10 μl iQ Multiplex Powermix (Bio-Rad, 1725849), with 600 nM of each primer (Sigma-Aldrich), 200 nM of each probe (Sigma-Aldrich) and 2 μl of extracted DNA. The amplification programme was initiated with a first step at 95 °C for 3 min, followed by 40 cycles of 10 s at 95 °C and 45 s at 60 °C. Primer efficiency was above 92% for both primer/probe combinations using standard curves of 10-fold dilutions of known input amounts. Data were analysed using Bio-Rad CFX Manager software.For the strain-specific analysis of both the pools and the individualized clones, we used two microsatellite loci, Ma30750 and Ma205451. Microsatellite locus Ma307 (forward: 5′-(6FAM)CATGCTCCGCCTTATTCCTC-3′, reverse: 5′-GGGTGGCGAAGAAGTAGACG-3′) allowed distinction of all strains except two of the M. brunneum strains (B1 and B3), which were distinguished by microsatellite locus Ma2054 (forward: 5′-(6FAM)GCCTGATCCAGACTCCCTCAGT-3′, reverse: 5′-GCTTTCGTACCGAGGGCG-3′). We analysed the microsatellites by E-Gel high-resolution 4% agarose gels (ILife Technologies, G501804) and fragment length analysis (done by Eurofins MWG) using Peak Scanner software 2.For clone individualization, we used flow cytometry to sort single spores out of the 40 spore pools (and the 6 ancestral strains for comparison) on 96-well plates (TPP; Biomedica, TP-92696) containing SDA (100 µl per well). The unstained spore population was detected using the FSC (forward scatter)/SSC (side scatter) in linear mode (70 μm nozzle, FACS ARIA III, BD Biosciences, as exemplified in Supplementary Fig. 1). Purity mode was set to ‘single cell’ and spore clones were obtained by sorting 1 particle event into each well. Sorting and data analysis were performed using Diva 6.2 software. The number of spores that we obtained for microsatellite analysis varied for each replicate, as it was adjusted to the remaining strain diversity estimate that we obtained from the quantitative and qualitative analysis of the pools. In total, we analysed 4–5 clones per ancestral strain (total n = 27) and a median of 5, but up to 101 different clones for the pools (total n = 1,347), as we intensified analysis for the strains that were revealed to be present at low frequency on the basis of previous analysis.Common garden experimentExperimental design and procedureWe then tested whether the successful lines at the end of the experiment (that is, after 10 host passages) differed in their virulence (induced host mortality) and investment into transmission stages (produced spore number) depending on their selection history (individual vs social), when current host social context either reflected the selection history or not. This common garden experiment thus led to 20 matched combinations of selection history and current condition (10 each of the individual lines in current individual host conditions (individual–individual) and the social lines in current social host conditions (social–social)) and 20 non-matched conditions (10 each of the individual lines in current social host conditions (individual–social) and the social lines in current individual host conditions (social–individual)).We obtained the lines for performance of the common garden experiment by the following procedure: (1) for the 16 out of the 20 replicate lines, where a single strain was the sole remaining representative at the end of the experiment (Fig. 1b), we expanded one of the c.f.u.s grown after FACS sorting (see above) by plating on SDA; (2) for the 4 remaining replicates in which two strains had remained (two individual and two social replicate lines), we expanded one c.f.u. of each of the remaining strains and mixed the spores in their representative proportion, as determined above.Virulence and transmissionFor the 10 individual and 10 social lines, we determined the induced host mortality as a measure of virulence and the outgrowing spore number as transmission stage production under their matched and non-matched current host conditions. We exposed the workers as in the selection treatment, kept them either alone or with two untreated nestmates, and monitored their mortality daily for 8 d. Again, ants dying in the first 24 h after treatment and dying nestmates were excluded from the analysis. In total, we obtained survival data of 797 ants (19–20 ants exposed for each of the 10 replicates from each of 4 combinations of selection history and current host condition). Dead ants were treated as above and their outgrowing spores collected by a needle dipped in sterile 0.05% Triton X-100 directly from the carcass, and resuspended in 100 µl of sterile 0.05% Triton X-100. The number of spores per carcass was counted individually using the automated cell counter, as described above (n = 215; median of 5 per replicate). We excluded one outlier carcass(from replicate I5) where we expected a counting error as this single carcass showed approx. 100-fold higher spore count than the other carcasses of this replicate. Exclusion of this outlier did not affect the statistical outcome. The proportion of ants dying per replicate line for each combination of selection history and current host condition and the number of spores produced by all carcasses per replicate were respectively used as measures of virulence and transmission (mean carcass spore load per replicate plotted in Fig. 2).Allogrooming elicitation by the fungal linesWe determined the allogrooming elicited by the individual and the social lines. To this end, we exposed workers as above and observed the allogrooming performed by two untreated nestmates towards the exposed ant. In detail, we performed 3 biological replicates for each of the 20 replicate lines (n = 10 individual and 10 social lines, resulting in a total of 60 videos), where the exposed ant was placed with two untreated nestmates within 10 min after exposure, and filmed with Ueye cameras for 30 min (whereby 4 cameras were used in parallel, each filming 3 replicates simultaneously, and using StreamPix 5 software (NorPix 2009-2001) for analysis). Videos were obtained in a randomized manner and labels did not contain treatment information so that the observer was blind to both the selection history and individual treatment during the behavioural annotations. For each ant, we observed both self- and allogrooming. Start and end times for each grooming event were determined, supported by use of the software BioLogic (Dimitri Missoh, 2010 (https://sourceforge.net/projects/biologic/)).As the ants in our serial passage and common garden experiments were not colour-marked, we also used unmarked ants for this behavioural experiment to keep conditions the same. This was possible as preliminary data with colour-coded nestmates (n = 18 videos) had shown that exposure alters the ant’s behaviour and that of its untreated nestmates in a predictable way that allows reliable classification of the pathogen-exposed individuals from the untreated nestmates; we used the following rules to classify an ant as the exposed individual: (1) the individual spent >5% more time (of the 30 min observation period) selfgrooming than the other individuals; (2) if the difference in selfgrooming time between the individuals was More