在草原植被中,植物种类和花卉特征对节肢动物群落的影响比新烟碱类农药污染更为显著

时间:2026年5月18日
来源:Ecosphere

编辑推荐:

摘要 农业实践是生物多样性丧失的重要驱动因素,导致各地区生态功能和服务下降。为了缓解这些影响,农业生态系统内的栖息地恢复已成为提高生态稳定性和增加生物多样性的策略。然而,新烟碱类杀虫剂在农业生态系统中的广泛使用引发了人们对它们对受恢复栖息地影响的担忧。本研究调查了新烟碱类污染对

广告
   X   

摘要

农业实践是生物多样性丧失的重要驱动因素,导致各地区生态功能和服务下降。为了缓解这些影响,农业生态系统内的栖息地恢复已成为提高生态稳定性和增加生物多样性的策略。然而,新烟碱类杀虫剂在农业生态系统中的广泛使用引发了人们对它们对受恢复栖息地影响的担忧。本研究调查了新烟碱类污染对早期草原恢复生态系统中节肢动物群落结构和植物适应性的多营养级影响。通过一个受控的野外实验和结构方程建模,我们评估了新烟碱类暴露如何影响不同摄食类的节肢动物数量、摄食类之间的数量相互作用以及通过种子产量和地上生物量衡量的植物适应性。我们的研究结果表明,新烟碱类污染对节肢动物摄食类的影响并不均匀;相反,植物特征是驱动群落组成和相互作用的主要因素。虽然在新烟碱类暴露的情况下,一个模型中传粉者的访问次数显著增加,而杂食者的数量显著减少,但植物种类和花序数量对节肢动物摄食类和群落结构的影响更大。这些结果支持了先前的研究,证明了植物群落组成对草原恢复稳定性的重要性,并表明扩大和增强恢复的草原栖息地可能比新烟碱类污染对生态系统结构和稳定性的影响更为显著。

引言

农业实践是生物多样性丧失的最大贡献者之一,会大幅减少一个地区内的生态位数量(Dudley & Alexander, 2017; Wagner et al., 2021)。为了应对这一下降趋势,草原保护工作,特别是在农业生态系统中,已成为通过增加生物多样性和发展有韧性的生态系统来提高生态稳定性和功能的普遍策略(Biondini, 2007; Raven & Wagner, 2021; Tilman et al., 1998; Outhwaite et al., 2022)。草原保护既包括在先前耕种的土壤上进行重建,也包括恢复退化但未耕种的地区——这里称为恢复。然而,这两种保护策略都可能受到遗留土壤残留物中的新烟碱类污染(Aseperi et al., 2020),或者由于在相邻农田中播种时种子涂层磨损产生的叶面喷雾或灰尘的漂移(Cheng, 2021; Nuyttens et al., 2013)的影响。这种污染可能会显著影响恢复区域改善生态系统稳定性的能力(Grman et al., 2010; Ruiz-Jaen & Aide, 2005)。鉴于节肢动物占据的生态位数量众多(Longcore, 2003),农药污染可能会通过不同地影响节肢动物的多样性并改变它们的生态相互作用,从而威胁到节肢动物提供的众多关键生态系统服务。全球使用最广泛的杀虫剂类别是新烟碱类杀虫剂(Hall et al., 2022; Humann-Guilleminot et al., 2019)。新烟碱类经常出现在恢复的栖息地中(Hall et al., 2022; Main et al., 2014),其环境持久性比其他传统杀虫剂更长(Khan, 1981; Rexrode et al., 2003),并且已被确定为昆虫数量下降的主要驱动因素(Forister et al., 2016, 2019; Goulson, 2013; Pisa et al., 2021; Wagner et al., 2021)。尽管如此,我们仍缺乏从系统角度全面了解节肢动物群落与恢复区内开花植物群落如何响应新烟碱类污染的理解。大多数关于节肢动物与新烟碱类的研究仅限于个别物种或狭窄的分类群,因此缺乏生态相关的野外背景(Disque et al., 2019; Main et al., 2018)。虽然孤立的研究发现了节肢动物对新烟碱类的特定反应(Kessler et al., 2015; Morales-Rodriguez & Peck, 2009; Nauen & Elbert, 1997),但这些孤立的结果未能描述新烟碱类对更广泛的生态系统层面的影响,如摄食类数量(即食草动物、捕食者、传粉者等)和多营养级相互作用。例如,Harmon等人(2023)观察到由于群落水平的新烟碱类污染,早期季节的食草动物生物量减少,捕食者生物量长期下降。这表明新烟碱类污染有可能通过影响某些营养级而破坏节肢动物群落的稳定性,进而改变恢复区提供的生态系统和农业服务。此外,鉴于节肢动物行为、生物量和多样性对新烟碱类暴露的广泛变化,可能会影响草原植物的适应性和生长。这些植物适应性的变化在农业商品中已有充分记录(Johnson et al., 2009; Krupke et al., 2017; Orlowski et al., 2016);然而,在自然系统中这些变化了解较少,可能会对恢复生态系统的稳定性产生巨大影响(Burghardt & Tallamy, 2013; Harvey et al., 2008; Schuldt et al., 2014)。为了更全面和系统地了解新烟碱类污染如何影响节肢动物群落及其提供的生态系统服务,我们使用了一个受控的野外实验来研究新烟碱类污染如何影响生态系统稳定性(以种子产量和植物地上生物量衡量)。这种设计允许更好地控制混杂变量,同时利用已建立的自然群落,从而保持其生物学相关性。在这里,我们探讨了(1)不同摄食类的节肢动物数量是否受到新烟碱类污染的不同影响,(2)新烟碱类污染是否改变了摄食类之间的相互作用,以及(3)这些因素是否影响植物适应性——在本研究中通过植物种子产量和地上生物量来表征。我们使用结构方程建模(SEM)来解析植物、节肢动物摄食类数量和新烟碱类之间的相互关系。这种方法使我们能够利用假设的物种间关系结构来研究新烟碱类在整个系统中的直接和间接影响,并评估每个连接在系统中的相对重要性。我们构建了两个SEM,分别研究植物地上生物量和种子产量,以更好地评估这些生态系统中草原植物的短期和长期建立情况(Biondini, 2007)。为了进一步研究趋势并在不同尺度上验证发现,我们使用非参数多元方差分析(PERMANOVA)来比较新烟碱类污染和植物种类对节肢动物家族群落结构的影响。这种方法评估了节肢动物群落结构的分类差异,而不是摄食类数量。这种实验设计结合这些分析技术,提供了新烟碱类污染对早期草原恢复群落中植物-昆虫相互作用稳定性的更全面的影响表示。

方法

实验设计

为了测试新烟碱类对昆虫访问植物的影响,我们使用了四种常见于高草草原恢复区的草原植物(Kurtz, 2013),这些植物被认为对传粉者有吸引力(Angelella & O'Rourke, 2017):Chamaecrista fasciculata(豆科)、Coreopsis tinctoria(菊科)、Rudbeckia hirta(菊科)和Monarda fistulosa(唇形科)。实验地块由56.78升的布质种植盆(247 Garden, Montebello, CA)组成,填充了伊利诺伊大学农业、消费者和环境科学学院(ACES)植物护理设施(Urbana, IL)准备的消毒通用土壤混合物(1:2:2,土壤:泥炭:珍珠岩),并在六月初种植了这四种植物。多年生植物物种(M. fistulosa和R. hirta)的一岁幼苗是从Prairie Nursery(Westfield, WI)购买的,而一年生植物则是从Prairie Moon Nursery(Winona, MN)购买的种子中培育的。两种处理方式,一种含有clothianidin,另一种作为对照,在Urbana, IL的Phillips Tract Natural Area重复了20次。采样设计和实际地点的照片见附录S1:图S1。包含实验地块的草原面积约为2.1公顷,于1976年从农业用途中恢复,这意味着它从未接受过新烟碱类的处理,新烟碱类是在20世纪90年代初引入的。该草原周围有大约105公顷的自然区域,包括原始森林、恢复的草原和管理林地。clothianidin污染的盆中的土壤在种植时用Arena 0.25G(Nufarm Americas Inc., Alsip, IL)处理,处理剂量是根据玉米施用率(每公顷80克活性成分)计算得出的(Badua et al., 2021)。选择这种施用率是为了代表从农业用地或受到大量新烟碱类灰尘漂移影响的栖息地重建的最坏情况(Nuyttens et al., 2013; Nuyttens & Verboven, 2015)。我们使用颗粒状杀虫剂作为农业种子涂层的代理,这是新烟碱类最主要的施用技术(Hladik et al., 2018)。这种实验设计捕捉了植物生活史特征(一年生和多年生)和开花时间的变化。它还允许更好地控制新烟碱类浓度、土壤组成和土壤微生物组的异质性,以减少数据中的混杂变量。通过控制这些因素和空间范围,我们可以观察到在污染广泛且异质的情况下可能不太明显的处理效果。

昆虫访问

从六月底开始开花时起,每隔一周记录一次昆虫发生情况,持续到2022年九月的第一周,共进行了七个采样期。每个采样期包括1分钟的观察者适应期(之前确定这是昆虫访问恢复的足够时间),然后是5分钟的观察期,记录每株植物上的活跃花序数量、每株植物上的所有节肢动物以及发现节肢动物的植物。节肢动物被鉴定到最低的分类群,并拍照辅助鉴定(Borror et al., 1998; Thomas et al., 2002; Whitfield et al., 2013)。为了限制植物损伤,没有收集昆虫凭证标本。

植物适应性采样

2022年九月末,在植物衰老后收集了地上生物量(包括土壤以上的所有植物组织)和平均种子产量。夏季有几株植物的盆被鹿损坏,因此被排除在所有分析之外。每个植物和处理的重复次数见附录S1:表S1。所有剩余植物的地上生物量都被采样。样本在60°C下干燥并称重。从随机选取的盆中收集花序(附录S1:表S1)。只收集了未受鹿损坏且完全衰老的植物的花序。如果可能,每株植物取三个花序,用于分析的平均值。花序被干燥、单独脱粒,并手动计数每个花序中的种子。M. fistulosa在整个实验期间没有开花;因此,没有评估该物种的种子产量。

统计分析—SEM

为了分析节肢动物摄食类数量、植物适应性结果和clothianidin土壤污染之间的群落级响应,我们采用了SEM。由于之前提到的地上生物量和种子产量之间的采样组差异,构建了两个SEM。此外,这些植物适应性指标可能受到不同摄食类的不同影响——这些摄食类利用植物的不同部分(Nemec et al., 2014)——并最终可能影响草原恢复的长期建立和繁衍。为每个SEM创建了单独的数据集,以便每个SEM分析只包括完整的地上生物量和种子产量样本。传统的SEM(使用如lavaan之类的R包)仅依赖于变量之间的方差-协方差矩阵,不考虑变量间误差方差的非正态性,无法解释样本的聚类(随机效应),并且复杂模型需要大量的自由度。因为我们的感兴趣的变量来自多个植物物种并且包含计数数据,这些数据需要不同的误差分布,所以我们使用了piecewiseSEM R包(Lefcheck, 2016)。Piecewise SEM将单独的广义线性混合效应模型结合成一个因果网络,可以适应随机效应结构和非正态分布(Shipley, 2009)。我们根据假设的变量关系构建了一个假设的SEM(附录S1:图S2)。基于这些假设关系,构建了广义线性混合效应模型(GLMMs),并在piecewise SEM分析中使用了它们。广义线性混合模型(GLMMs)是使用lme4包中的glmer函数(版本1.1-2;Bates等人,2015年)拟合的,旨在建立用于选择 feeding guilds 丰富度、地上植物生物量和每朵花平均种子产量的全局模型。feeding guilds的丰富度和每个植物及样地的活跃花序数量被汇总到一年中,以便在单一分析中同时考虑花卉适应性和节肢动物发生率。模型选择是针对每个数据子集分别进行的。我们原始结构方程模型(SEM)中假设的所有feeding guilds、clothianidin处理和花序数量之间的交互作用都被包含在每个全局模型中。植物物种也被作为随机效应包含在每个模型中。这使我们能够在考虑植物物种间种子产量和地上生物量差异的同时,研究这个模拟系统中的更广泛趋势。然后使用MuMin包中的dredge函数(Bartoń,2023年)来确定我们SEM中每个因变量的最佳模型。接着使用psem函数(Lefcheck,2016年)对GLMMs进行分段SEM分析。使用fisherC函数(Lefcheck,2016年)评估了整体模型的拟合优度。只有p值大于0.05的模型被认为是对数据的良好拟合。通过计算方差膨胀因子(VIFs)来评估自变量之间的多重共线性。所有模型的VIF分数显示预测变量之间的相关性较低,表明我们的协变量足够独立。所有模型的残差都通过视觉检查和Shapiro–Wilk正态性检验进行了评估。异方差性也是通过视觉检查和Breusch–Pagan检验来评估的。

为了评估家族级别的节肢动物群落差异,我们使用了vegan包中的adonis2函数(Oksanen等人,2024年)进行了非PERMANOVA分析,进行了999次排列。比较了clothianidin处理组和对照组、不同植物物种之间以及植物物种与处理组之间的群落结构。排列使用日期作为分层项来考虑季节性变化。对于所有比较,都使用了成对Bray–Curtis相似性指数。对所有相应的PERMANOVA分析进行了方差同质性评估。使用非度量多维缩放(NMDS)排序可视化了clothianidin处理组和对照组以及不同植物物种之间的差异。所有统计建模都是在R v4.3.2(R核心团队,2023年)中进行的。

在研究的七个采样期间,共记录了1533次节肢动物访问,代表了60个家族。观察到的最常见家族包括Formicidae(对照组:347次,clothianidin处理组:243次)、Chrysomelidae(对照组:86次,clothianidin处理组:116次)、Halictidae(对照组:47次,clothianidin处理组:58次)、Carabidae(对照组:39次,clothianidin处理组:24次)、Cicadellidae(对照组:27次,clothianidin处理组:20次)和Reduviidae(对照组:23次,clothianidin处理组:14次)。观察到的所有标本的完整列表可以在伊利诺伊数据银行找到(https://doi.org/10.13012/B2IDB-8830381_V1)。我们的SEM显示,clothianidin土壤污染并没有均匀影响节肢动物feeding guilds或植物适应性指标(图1和图2)。相反,植物特征中的花序丰富度是影响节肢动物feeding guilds丰富度的更重要因素(图1和图2)。

结构方程模型展示了clothianidin土壤污染、节肢动物feeding guilds丰富度以及生长季节末每朵花平均种子产量之间的关系。箭头表示变量之间的方向性效应。实线表示统计上显著的关系(p<0.05),而虚线表示不显著的关系(p>0.05)。箭头宽度根据标准化系数进行缩放(表1),并根据其效应的方向进行着色(蓝色表示正面效应;红色表示负面效应)。还报告了每个响应变量的边际(仅固定效应)和条件(固定和随机效应)R2值;条件值用括号表示。

传粉者丰富度与clothianidin污染呈正相关(t44=2.656,p=0.008)(表1,图1)。此外,clothianidin对平均种子产量有正面影响,这种影响独立于节肢动物的交互作用(t44=9.373,p=0.000)。除了clothianidin污染外,花序丰富度对feeding guilds丰富度的影响最大。传粉者、捕食者和杂食者的丰富度都与花序丰富度呈正相关(t44=2.037,p=0.042;t44=2.037,p=0.0003;t44=3.384,p=0.0007)。植食者的丰富度与花序丰富度呈正相关,但相关性不显著(t44=1.302,p=0.193)。节肢动物内部群落之间的关联没有显著的因果关系;然而,捕食者的丰富度对杂食者的丰富度有边际显著的负面影响(t44=-1.675,p=0.094)。种子产量与传粉者丰富度、捕食者丰富度和clothianidin污染呈正相关(t44=5.845,p=0.000;t44=5.674,p=0.000;t44=9.373,p=0.000),并与植食者丰富度和花序丰富度呈负相关(t44=-7.140,p=0.000;t44=-3.419,p=0.0006)。

在研究的七个采样期间,共记录了1533次节肢动物访问,代表了60个家族。观察到的最常见家族包括Formicidae(对照组:347次,clothianidin处理组:243次)、Chrysomelidae(对照组:86次,clothianidin处理组:116次)、Halictidae(对照组:47次,clothianidin处理组:58次)、Carabidae(对照组:39次,clothianidin处理组:24次)、Cicadellidae(对照组:27次,clothianidin处理组:20次)和Reduviidae(对照组:23次,clothianidin处理组:14次)。观察到的所有标本的完整列表可以在伊利诺伊数据银行找到(https://doi.org/10.13012/B2IDB-8830381_V1)。我们的SEM显示,clothianidin土壤污染并没有均匀影响节肢动物feeding guilds或植物适应性指标(图1和图2)。相反,植物特征中的花序丰富度是影响节肢动物feeding guilds丰富度的更重要因素(图1和图2)。

结构方程模型展示了clothianidin土壤污染、节肢动物feeding guilds丰富度以及生长季节末每朵花平均种子产量之间的关系。箭头表示变量之间的方向性效应。实线表示统计上显著的关系(p<0.05),而虚线表示不显著的关系(p>0.05)。箭头宽度根据标准化系数进行缩放(表1),并根据其效应的方向进行着色(蓝色表示正面效应;红色表示负面效应)。还报告了每个响应变量的边际(仅固定效应)和条件(固定和随机效应)R2值;条件值用括号表示。

传粉者丰富度与clothianidin污染呈正相关(t44=2.656,p=0.008)(表1,图1)。此外,clothianidin对平均种子产量有正面影响,这种影响独立于节肢动物的交互作用(t44=9.373,p=0.000)。除了clothianidin污染外,花序丰富度对feeding guilds丰富度的影响最大。传粉者、捕食者和杂食者的丰富度都与花序丰富度呈正相关(t44=2.037,p=0.042;t44=2.037,p=0.0003;t44=3.384,p=0.0007)。植食者的丰富度与花序丰富度呈正相关,但相关性不显著(t44=1.302,p=0.193)。节肢动物内部群落之间的关联没有显著的因果关系;然而,捕食者的丰富度对杂食者的丰富度有边际显著的负面影响(t44=-1.675,p=0.094)。种子产量与传粉者丰富度、捕食者丰富度和clothianidin污染呈正相关(t44=5.845,p=0.000;t44=5.674,p=0.000;t44=9.373,p=0.000),并与植食者丰富度和花序丰富度呈负相关(t44=-7.140,p=0.000;t44=-3.419,p=0.0006)。

在研究的七个采样期间,共记录了1533次节肢动物访问,代表了60个家族。观察到的最常见家族包括Formicidae(对照组:347次,clothianidin处理组:243次)、Chrysomelidae(对照组:86次,clothianidin处理组:116次)、Halictidae(对照组:47次,clothianidin处理组:58次)、Carabidae(对照组:39次,clothianidin处理组:24次)、Cicadellidae(对照组:27次,clothianidin处理组:20次)和Reduviidae(对照组:23次,clothianidin处理组:14次)。观察到的所有标本的完整列表可以在伊利诺伊数据银行找到(https://doi.org/10.13012/B2IDB-8830381_V1)。我们的SEM显示,clothianidin土壤污染并没有均匀影响节肢动物feeding guilds或植物适应性指标(图1和图2)。相反,植物特征中的花序丰富度是影响节肢动物feeding guilds丰富度的更重要因素(图1和图2)。

结构方程模型展示了clothianidin土壤污染、节肢动物feeding guilds丰富度以及生长季节末每朵花平均种子产量之间的关系。箭头表示变量之间的方向性效应。实线表示统计上显著的关系(p<0.05),而虚线表示不显著的关系(p>0.05)。箭头宽度根据标准化系数进行缩放(表1),并根据其效应的方向进行着色(蓝色表示正面效应;红色表示负面效应)。还报告了每个响应变量的边际(仅固定效应)和条件(固定和随机效应)R2值;条件值用括号表示。

传粉者丰富度与clothianidin污染呈正相关(t44=2.656,p=0.008)(表1,图1)。此外,clothianidin对平均种子产量有正面影响,这种影响独立于节肢动物的交互作用(t44=9.373,p=0.000)。除了clothianidin污染外,花序丰富度对feeding guilds丰富度的影响最大。传粉者、捕食者和杂食者的丰富度都与花序丰富度呈正相关(t44=2.037,p=0.042;t44=2.037,p=0.0003;t44=3.384,p=0.0007)。植食者的丰富度与花序丰富度呈正相关,但相关性不显著(t44=1.302,p=0.193)。节肢动物内部群落之间的关联没有显著的因果关系;然而,捕食者的丰富度对杂食者的丰富度有边际显著的负面影响(t44=-1.675,p=0.094)。种子产量与传粉者丰富度、捕食者丰富度和clothianidin污染呈正相关(t44=5.845,p=0.000;t44=5.674,p=0.000;t44=9.373,p=0.000),并与植食者丰富度和花序丰富度呈负相关(t44=-7.140,p=0.000;t44=-3.419,p=0.0006)。

为了评估家族级别的节肢动物群落差异,我们使用了vegan包中的adonis2函数(Oksanen等人,2024年)进行了非PERMANOVA分析,进行了999次排列。比较了clothianidin处理组和对照组、不同植物物种之间以及植物物种与处理组之间的群落结构。排列使用日期作为分层项来考虑季节性变化。对于所有比较,都使用了成对Bray–Curtis相似性指数。对所有相应的PERMANOVA分析进行了方差同质性评估。clothianidin处理组和对照组之间以及不同植物物种之间的差异通过非度量多维缩放(NMDS)排序进行了可视化。所有统计建模都是在R v4.3.2(R核心团队,2023年)中进行的。

在研究的七个采样期间,共记录了1533次节肢动物访问,代表了60个家族。观察到的最常见家族包括Formicidae(对照组:347次,clothianidin处理组:243次)、Chrysomelidae(对照组:86次,clothianidin处理组:116次)、Halictidae(对照组:47次,clothianidin处理组:58次)、Carabidae(对照组:39次,clothianidin处理组:24次)、Cicadellidae(对照组:27次,clothianidin处理组:20次)和Reduviidae(对照组:23次,clothianidin处理组:14次)。观察到的所有标本的完整列表可以在伊利诺伊数据银行找到(https://doi.org/10.13012/B2IDB-8830381_V1)。我们的SEM显示,clothianidin土壤污染并没有均匀影响节肢动物feeding guilds或植物适应性指标(图1和图2)。相反,植物特征中的花序丰富度是影响节肢动物feeding guilds丰富度的更重要因素(图1和图2)。

结构方程模型展示了clothianidin土壤污染、节肢动物feeding guilds丰富度以及生长季节末每朵花平均种子产量之间的关系。箭头表示变量之间的方向性效应。实线表示统计上显著的关系(p<0.05),而虚线表示不显著的关系(p>0.05)。箭头宽度根据标准化系数进行缩放(表1),并根据其效应的方向进行着色(蓝色表示正面效应;红色表示负面效应)。还报告了每个响应变量的边际(仅固定效应)和条件(固定和随机效应)R2值;条件值用括号表示。

传粉者丰富度与clothianidin污染呈正相关(t44=2.656,p=0.008)(表1,图1)。此外,clothianidin对平均种子产量有正面影响,这种影响独立于节肢动物的交互作用(t44=9.373,p=0.000)。除了clothianidin污染外,花序丰富度对feeding guilds丰富度的影响最大。传粉者、捕食者和杂食者的丰富度都与花序丰富度呈正相关(t44=2.037,p=0.042;t44=2.037,p=0.0003;t44=3.384,p=0.0007)。植食者的丰富度与花序丰富度呈正相关,但相关性不显著(t44=1.302,p=0.193)。节肢动物内部群落之间的关联没有显著的因果关系;然而,捕食者的丰富度对杂食者的丰富度有边际显著的负面影响(t44=-1.675,p=0.094)。种子产量与传粉者丰富度、捕食者丰富度和clothianidin污染呈正相关(t44=5.845,p=0.000;t44=5.674,p=0.000;t44=9.373,p=0.000),并与植食者丰富度和花序丰富度呈负相关(t44=-7.140,p=0.000;t44=-3.419,p=0.0006)。

表1. 来自平均种子产量结构方程模型的回归系数估计值及其标准误差(图1)。

预测变量 估计值 标准误差 自由度 t值 p值 标准化估计值

| 传粉者丰富度 | 0.6687 | 0.2517 | 44 | 2.6564 | 0.0079 |
| 花序丰富度 | 0.0156 | 0.0077 | 44 | 2.0372 | 0.0416 |
| 捕食者丰富度 | 0.0320 | 0.0613 | 44 | 0.5223 | 0.6014 |
| 杂食者丰富度 | 0.0040 | 0.0115 | 44 | 0.3491 | 0.7270 |
| 植食者丰富度 | 0.1958 | 0.1920 | 44 | 1.0200 | 0.3077 |
| 捕食者丰富度 | 0.0695 | 0.0447 | 44 | 1.5539 | 0.1202 |
| 花序丰富度 | 0.0084 | 0.0065 | 44 | 1.3024 | 0.1928 |
| 杂食者丰富度 | 0.0775 | 0.2036 | 44 | 0.3804 | 0.7036 |
| 捕食者丰富度 | -0.1351 | 0.0807 | 44 | -1.6751 | 0.0939 |
| 花序丰富度 | 0.0248 | 0.0073 | 44 | 3.3838 | 0.0007 |
| 捕食者丰富度 | 0.1054 | 0.2073 | 44 | 0.5083 | 0.6112 |
| 花序丰富度 | 0.0208 | 0.0058 | 44 | 3.6131 | 0.0003 |
| 平均种子产量 | 0.0345 | 0.0059 | 44 | 5.8449 | 0.0000 |
| 植食者丰富度 | -0.0290 | 0.0041 | 44 | -7.1404 | 0.0000 |
| clothianidin | 0.2044 | 0.0218 | 44 | 9.3732 | 0.0000 |
| 杂食者丰富度 | 0.0096 | 0.0056 | 44 | 1.7054 | 0.0881 |
| 捕食者丰富度 | 0.0319 | 0.0056 | 44 | 5.6741 | 0.0000 |
| 花序丰富度 | -0.0035 | 0.0010 | 44 | -3.4185 | 0.0006 |

注:显著的标准化系数以粗体显示(p<0.05)。除了上述直接效应外,clothianidin污染和花序丰富度还通过传粉者丰富度间接对平均种子产量产生了正面影响(图1)。花序丰富度通过传粉者和捕食者丰富度对平均种子产量产生了正面间接影响。

杂食者丰富度是唯一受到clothianidin污染显著影响的feeding guild(t80=-2.007,p=0.045)(表2,图2)。花序丰富度与捕食者丰富度、杂食者丰富度和地上生物量呈正相关(t80=3.325,p=0.001;t80=3.263,p=0.001;t80=3.182,p=0.002;t80=3.840,p=0.0003)。除了直接效应外,花序丰富度还通过植食者丰富度对地上生物量产生了正面间接影响。没有发现节肢动物内部群落之间的显著关联。植食者丰富度与地上生物量呈正相关(t80=3.884,p=0.0002)。

表2. 来自地上生物量结构方程模型的回归系数估计值及其标准误差(图2)。

注:显著的标准化系数以粗体显示(p<0.05)。

节肢动物群落组成——通过我们的PERMANOVA分析确定——受到植物物种的显著影响(F1,59=8.584,p=0.001),但未受到clothianidin污染的显著影响(F1,59=0.823,p=0.552),也未受到植物物种与clothianidin污染交互作用的显著影响(F1,59=0.821,p=0.691;表3)。方差同质性分析显示,clothianidin处理组的离散度没有显著差异(F1,59=0.059,p=0.771)(附录S1:表S2)。然而,植物PERMANOVA的方差同质性分析显示分散度不均匀(F1,59=7.342,p<0.001),这可能会影响估计值和显著性值。尽管如此,NMDS排序(图3)支持了一个聚类模式,证实了植物物种对节肢动物群落结构的重要性。

表3. 参数化多元方差分析的估计值、拟合优度和p值(图3)。

注:显著的标准化系数以粗体显示(p<0.05)。

在本次原位操纵的草原保护生态系统中进行的实验中,我们发现,总体而言,花卉特征而非新烟碱类农药污染是影响节肢动物feeding guilds丰富度的主要因素。同样,花卉物种显著影响了节肢动物家族结构,而clothianidin对节肢动物家族结构没有显著影响。然而,我们发现clothianidin污染通过直接增加种子产量以及通过传粉者丰富度间接增加种子产量,对草原草本植物的种子产量产生了正面影响。这些结果对于草原恢复和重建工作具有广泛的相关性,并为描述农业生态系统中早期草原建立期间生态重要群体之间群落级相互作用的相对重要性提供了有价值的框架。

恢复生态系统的长期栖息地稳定性在很大程度上依赖于草原植物建立和繁殖的能力(Benson & Hartnett,2006)。我们发现,clothianidin污染在草原植物早期建立阶段可以直接和间接通过增强传粉者的招募来影响种子产量。这两种直接和间接效应在图1的实线中都有体现。直接效应由clothianidin和种子产量之间的实线表示,间接效应由clothianidin、传粉者和种子产量之间的实线表示。这表明传粉者被吸引到受污染的花朵上,从而对种子产量产生了后续的正面影响。这些结果表明,当存在新烟碱类农药时,开花植物在新建立的恢复区中的繁殖能力得到了偶然的增强。蜜蜂是研究中观察到的主要传粉者,占总访问次数的82%以上。蜜蜂与新烟碱类农药之间的这种正面关联证实了早期研究的结果,这些研究表明蜜蜂倾向于消耗受污染的资源(Arce等人,2018;Kessler等人,2015)并在受污染的区域筑巢(Tetlie & Harmon-Threatt,2024)。尽管新烟碱类农药污染可能改善了觅食栖息地,但长期暴露于这些化学物质可能对蜜蜂群落健康有害。新烟碱类农药的暴露会降低非蜜蜂物种的繁殖力(Sandrock等人,2014)、幼虫发育速率(Anderson & Harmon-Threatt,2019)、成虫寿命(Anderson & Harmon-Threatt,2019)和种群增长率(Stuligross & Williams,2021)。因此,尽管clothianidin污染可能有利于传粉者对草原花朵的招募和随后的花卉适应性,但长期暴露于新烟碱类农药污染的栖息地可能会对传粉者群落产生负面影响。虽然捕食者和植食者在这项研究中没有受到新烟碱类农药的影响,但它们分别对种子产量有显著的正面和负面影响。这些发现与其他研究结果一致,表明新烟碱类农药污染并不会严重破坏这些生态群体,特别是在恢复栖息地中发现的水平(Alford & Krupke,2017;Ding等人,2018)(附录S1:表S3)。虽然捕食者不受clothianidin的影响,但它们受到花序丰富度的正面影响。这种正面关联可能是由常见于花朵上的节肢动物捕食者驱动的,如蟹蛛(Thomisidae)、跳蛛(Salticidae)和埋伏虫(Phymata sp.)。这些捕食者经常捕食传粉者(Gil-Santana等人,2015;Morse,1984;Trindade-Santos等人,2023)——它们也受到花序丰富度的正面影响。然而,我们在扫描电子显微镜(SEM)分析中并未观察到传粉者与捕食者之间存在任何显著的相互作用。捕食者与传粉者之间无显著差异,加上捕食者与食草动物之间也无显著差异,这可能是由于将与花朵相关的捕食者与非花朵相关的捕食者混合在一起所导致的稀释效应。尽管这些节肢动物食物群的丰度趋势具有参考价值,但重要的是要认识到,这种类型的分析无法完全涵盖这些群体的分类范围。很可能这些食物群中的不同物种受到噻虫胺污染的影响程度不同,这可能会导致节肢动物群落的长期变化以及随之而来的生态系统服务的变化。因此,需要进一步的野外研究来更好地了解这些恢复系统中与新烟碱类农药污染相关的节肢动物种群趋势。虽然我们的研究结果表明噻虫胺污染可能会提高植物的适应性,但成功的植物建立并不能保证生态功能的恢复(Hansen & Gibson, 2014)。噻虫胺与某些食物群丰度之间缺乏因果关系,可能未能捕捉到节肢动物可能经历的潜在亚致死效应,这些效应可能会对节肢动物群落结构和功能产生长期影响。

地上生物量指标对于评估草原植物的短期和长期生产力和适应性非常有用(Biondini, 2007)。已有研究表明新烟碱类农药可以增加某些农作物的株高(Preetha & Stanley, 2012),但这一点在草原植物中尚未得到验证。我们的结果显示,噻虫胺污染对所研究的草原草本植物的地上生物量既没有显著的直接影响,也没有间接影响。尽管我们预期食草动物在受污染植物上的数量会减少,但我们发现它们并未受到噻虫胺污染的显著影响。这一发现与一些行作物研究的结果一致(Krupke et al., 2017);然而,这些结果并不一致,因为其他研究显示新烟碱类农药对害虫具有有效的控制作用(Kundoo et al., 2018; Lee & Davis, 2023)。此外,需要指出的是,我们划分的食草动物群体实际上包含了多种可能对新烟碱类农药反应不同的物种。因此,需要不同的实验设计和统计分析来直接在物种层面分析这些趋势。虽然我们植物中的噻虫胺含量(附录S1:表S3)在恢复生境中可能被认为是较高的(Cheng, 2021),但在行作物农业中则不足以构成显著影响(Alford & Krupke, 2017)。事实上,某些作物害虫,如黑切根虫(Agrotis ipsilon),在受到亚致死剂量噻虫胺处理时,其繁殖率会增加(Ding et al., 2018)。虽然这不能代表所有食草动物,但某些叶食性昆虫,如A. ipsilon,对新烟碱类农药的敏感性远低于其他节肢动物群体,如蜜蜂(Beadle et al., 2019; Costa et al., 2015; Mayack & Boff, 2019),这表明节肢动物食物群可能受到新烟碱类农药污染的不同影响。这些差异可能由多种因素造成,包括基于进食习惯的不同暴露途径(Bonmatin et al., 2015)以及解毒新烟碱类农药的基因组能力差异(Ahn et al., 2012; Beadle et al., 2019; Berenbaum & Johnson, 2015; Hayward et al., 2019)。虽然需要更高的分类分辨率来确定噻虫胺污染如何影响长期食草动物结构,但我们的研究表明,总体而言,食草动物的访问频率并未受到我们研究中使用的噻虫胺水平的影响。捕食者的数量同样没有受到噻虫胺污染的显著影响。关于新烟碱类农药污染对捕食者的影响,目前的研究结果尚无定论,因此需要进一步研究(Prabhaker et al., 2011, 2017; Rix & Cutler, 2020)。新烟碱类农药对捕食者的亚致死效应包括对 Pentatomid 捕食者繁殖的促进作用(Rix & Cutler, 2020),导致捕食者总体生物量的显著减少(Harmon et al., 2023),以及严重影响某些捕食性甲虫的移动和摄食能力(Tooming et al., 2017)。尽管许多这些发现是特定于某些类群的,并且常常缺乏野外环境的生物学背景,但这些结果的变化可能预示着这些受污染生态系统中的节肢动物食物群结构的重组或不稳定。因此,需要进行长期研究以评估这些变化的趋势。虽然在我们的种子设定模型中这一效应不显著(见上文),但在地上生物量模型中,杂食动物的数量确实因噻虫胺而显著减少。这一趋势可能是由于蚂蚁数量的增加所致,特别是那些在 C. fasciculata 上发现的蚂蚁,它们占本研究中观察到的所有杂食动物的66%。蚂蚁通常会在 C. fasciculata 的花外蜜腺(EFNs)上活动,并经常提供对食草动物的防御(Fehling, 2022)。除了总体数量减少外,我们没有观察到由于蚂蚁数量增加而导致的其他食草动物减少的情况。在这些 EFNs 上觅食的蚂蚁可能比不直接摄取花蜜的杂食动物暴露于更多的新烟碱类农药中。观察到的杂食动物减少与先前研究表明的蚂蚁暴露于新烟碱类农药后觅食减少的结果一致(Thiel & Köhler, 2016)。此外,与其他真社会性膜翅目昆虫一样,长期暴露于新烟碱类农药会降低蚁后的繁殖能力(Schläppi et al., 2020; Wu-Smart & Spivak, 2016)。蚂蚁的负面影响可能特别严重,因为蚂蚁在草原生态系统中占据了许多生态位并提供了多种生态系统服务(Nemec, 2014; Wills & Landis, 2018)。尽管新烟碱类农药暴露的影响不显著,但地上生物量仍受到花序数量的直接影响和间接影响。花序数量对地上生物量的间接影响是通过食草动物数量介导的,这与早期研究结果一致,即花序数量的增加与食草活动的增加相关(Sletvold & Grindeland, 2008)。此外,食草活动可以在多种植物中引发补偿性生长,在某些情况下甚至过度补偿(Garcia & Eubanks, 2019)。花序数量与地上生物量之间的正相关关系也可能仅仅与较大植物尺寸带来的更多资源有关(Younginger et al., 2017)。

与我们的 SEM 结果一致,这些结果强调了植物对节肢动物的重要性。我们的 PERMANOVA 分析表明,节肢动物群落在科水平上受到植物种类的显著影响,但未受到噻虫胺污染或植物种类与噻虫胺之间相互作用的影响。此外,这表明节肢动物的植物资源和潜在生态位在不同植物科之间存在差异,从而导致不同植物之间的节肢动物群落关联存在显著差异。这一点可以通过比较菊科植物(其节肢动物群落结构最为相似,见图3)与其他研究中的植物来看出。先前的研究支持了植物种类组成在塑造节肢动物群落结构中的重要性(Schaffers et al., 2008),尽管我们在研究中添加了新烟碱类农药,但仍然得到了类似的结果。我们认为,改善植物群落组成是促进早期恢复区和退化景观中节肢动物群落恢复力和稳定性的最重要因素。尽管由于节肢动物的分组方式不同,PERMANOVA 和 SEM 的结果不能直接比较,但两项分析都表明,植物种类和花特征比新烟碱类农药污染对访问节肢动物群落结构和食物群丰度的影响更大。

新烟碱类农药对新的和现有的自然栖息地的污染可能会持续,并且在未来可能还会增加,因为美国目前种植的大部分玉米和大豆田都使用了某种形式的新烟碱类农药(Douglas & Tooker, 2015),并且每粒种子的施用率有增加的趋势(Tooker et al., 2017)。我们的研究表明,在农业生态系统中建立多样化的植物群落对于节肢动物群落结构和食物群丰度的稳定性更为关键,而不是新烟碱类农药污染可能带来的负面影响。然而,由于亚致死剂量新烟碱类农药暴露的潜在代际影响,需要进一步的纵向研究来更准确地评估新烟碱类农药污染对恢复生态系统的潜在影响。虽然本研究中没有进行相关研究,但新烟碱类农药在这些系统中影响较小的一个可能解释是,经过多年的栖息地退化、自然区域丧失和持续使用农药后,这些地区的昆虫群落已经减少为耐干扰的物种。

**作者贡献**
Jonathan Tetlie:概念化;数据管理;正式分析;资金获取;调查;方法论;项目管理;资源;软件;验证;可视化;初稿撰写;审稿和编辑。
Alexandra Harmon-Threatt:概念化;资金获取;项目管理;资源;监督;验证;审稿和编辑。

**致谢**
作者感谢 Adam Dolezal、May Berenbaum 和 Anthony Yannarell 在项目建立过程中的贡献,感谢 Jamie Ellis 和 Nate Hudson 在选址和建立研究地点方面的帮助,以及 Nick Anderson、Scott Clem、Marissa Chase、Annaliese Wargin、Virginia Roberts、Kristine Schoenecker、Adrien Seabloom、Timo Wayman、Amauri Middleton、Jacob Chow 和 Jenaya Wilder 对手稿提供的反馈。作者还要感谢以下资金来源:本研究得到了 USDA-NIFA 奖项的支持:2018-67013-27537 和 2022-67013-36273(授予 Alexandra Harmon-Threatt)以及 2023-67011-40487(授予 Jonathan Tetlie);以及伊利诺伊大学厄巴纳-香槟分校综合生物学学院的支持(Fred H. Schmidt 夏季学者奖、Francis M. 和 Harlie M. Clark 研究支持基金以及 Harley J. Van Cleave 研究奖)。

**利益冲突声明**
作者声明没有利益冲突。

**数据可用性声明**
数据和代码(Tetlie & Harmon-Threatt, 2026)可从伊利诺伊数据银行获取:https://doi.org/10.13012/B2IDB-8830381_V1。

生物通微信公众号
微信
新浪微博


生物通 版权所有