Fitness cost analysis of resistance to dimefluthrin in Aedes albopictus

敏感品系白纹伊蚊（SS）由广东省中山市疾病预防控制中心提供。抗药品系白纹伊蚊（RR）来自广东省中山市野外采集，在华南农业大学饲养，将幼虫暴露于四氟甲醚菊酯溶液中，筛选出具有抗性的白纹伊蚊（2014年至今）。在此过程中，由本实验室人工饲养多代，浓度对其幼虫保持在半致死剂量，通过连续筛选获得稳定的遗传抗药种群（抗性倍数为27.02），其饲养环境保持在光照周期为白天∶夜晚=14 h∶10 h，温度为27±1 ℃，相对湿度为（75±5）%的室内环境。实验室采用标准化昆虫饲养^[6]。幼虫饲养于22 cm×14 cm×6.3 cm的饲养盆中，加入盆容量3/5的去氯自来水，将小鼠粮切成小块，每日投喂定量小鼠粮，每2天更换一次去氯自来水，保证水质干净清澈。幼虫化蛹后，使用一次性吸管将蛹吸出，置于装有去氯自来水的100 mL小烧杯中，将烧杯置于23 cm×23 cm×23 cm的养虫笼中，蚊蛹羽化为成蚊后，喂食10%的葡萄糖水（将棉花芯浸渍在葡萄糖水中）。成蚊在羽化后的1～2天后，每日16:00～18:00喂食小白鼠血供雌蚊产卵。

参考赵紫华（2020）的生命表方法^[7] 组建SS和RR品系白纹伊蚊的试验种群生命表，以羽化为发育期起始点，计算种群净增殖率（R₀）、世代平均周期（T）、内禀增长率（r_m）、周限增长率（λ）和种群加倍时间（D）等参数，计算公式如下：

式中x为以d为单位的时间间隔，平均后代数m_x是指成虫羽化后第x天,平均每头雌虫后代数^[8]。

将两种品系卵块分别泡在去氯水中，孵化出幼虫后，用一次性吸管各取50只SS组和RR组幼虫，相同环境下采用标准化昆虫饲养。从一龄幼虫期开始每天记录，内容包括幼虫的数量，蛹的数量，幼虫的发育时长，蛹的发育时长。至到所有幼虫羽化成成蚊，停止记录。

SS品系和RR品系的白纹伊蚊幼虫在相同的条件下饲养，羽化为成蚊后，用吸蚊器各取50头24 h内羽化出的成蚊，放入养虫笼中观察。养虫笼中每天添加干净的10%的葡萄糖水，并记录两种品系死亡的成蚊数量，直至所有成蚊死亡。

用Excel记录数据，SPSS 19.0分析数据，GraphPad Prism 9绘制图表。一龄期幼虫发育至蛹的平均时间、蛹羽化的平均时间、一龄期幼虫发育至羽化的平均时间用t检验进行判定，p<0.05被认为有统计学差异；幼虫化蛹率和羽化率之间的差异用t检验进行判定，p<0.05被认为有统计学差异；成蚊生存时间之间的差异采用Kaplan-Meier生存分析和log-rank检验进行判定，生存时间中位数之间的差异用Mann-Whitney U检验判定，p<0.05被认为有统计学差异^[9]。

敏感品系种群净增长率（R₀）、世代平均周期（T）、内禀增长率（r_m）、周限增长率（λ）、种群加倍时间（D）分别为41.13、11.77 d、0.32 d⁻¹、1.38 d⁻¹、2.17 d；抗性种群分别为26.48、14.27 d、0.23 d⁻¹、1.26 d⁻¹、3.01 d。敏感品系白纹伊蚊较抗性品系的种群加倍时间短，表明其繁殖速度较快。

SS组和RR组的化蛹率和羽化率均为100%，两个品系的化蛹率和羽化率均无显著差异（p>0.05）。SS和RR从一龄幼虫发育至蛹的平均时间分别9.66±0.17 d和11.31±0.17 d，蛹的平均羽化时间分别为2.15±0.04 d和2.90±0.03 d（图1），结果显示RR组一龄幼虫发育至蛹所需时间比SS组长，差异有统计学意义（p<0.005），且蛹发育至成蚊所需时间也具有统计学差异（p<0.05）。总的来说，RR组从一龄幼虫发育至成蚊的时间比SS组长，分别为14.27±0.14 d和11.77±0.16 d，说明四氟甲醚菊酯抗性的产生对白纹伊蚊的发育产生影响，并延长了白纹伊蚊幼虫发育至成蚊的时间。

Figure 1.  Life table analysis for the susceptible and resistant lines of Aedes albopictus from larvae to adults

从SS组和RR组雌蚊的生存曲线来看，敏感品系和抗性品系雌蚊的平均生存时间分别为30.03±1.60 d和21.83±1.12 d，生存时间中位数分别为32.50±0.69 d和22.00±1.09 d（表2, 图2-b），有统计学差异（c2=21.954, df=1, p<0.001, Mann-Whitney U test p<0.001）。结合成蚊的生存曲线，说明四氟甲醚菊酯抗性的产生对白纹伊蚊雌性成蚊的寿命产生了影响，缩短了雌性成蚊的寿命。

长期、广泛使用杀虫剂会导致白纹伊蚊产生相应的抗药性，抗药性的产生往往伴随着适合度代价的出现^[10]。昆虫产生的适合度代价通常与它们对某种环境压力或杀虫剂的抗性发展相关^{[11, 12]}。这种代价可能体现在昆虫滞育、寿命缩短、雌性或者雄性的生殖力的下降等，即昆虫在受到外界某种压力时，其生存和繁殖能力受到一定程度的影响^[13]。昆虫产生适合度代价的现象和背后的机制是非常复杂的，不同品种的昆虫对不同杀虫剂产生抗性后所表现出的适合度代价也各不相同。例如小菜蛾对阿维菌素和多杀菌素产生抗性后，分别表现出繁殖力降低和耐低温能力下降等适合度代价^[15]。烟粉虱对新烟碱类药剂产生抗性后，出现卵巢发育异常的生殖代价^[15]。同时也存在现象相反情况，例如褐飞虱对吡虫啉产生抗性后，表现出了适合度下降的趋势^[16]。不同的适合度代价意味着需要采取不同的治理策略来应对昆虫抗药性问题。目前已有很多关于对拟除虫菊酯类杀虫剂产生抗药性的媒介蚊虫的研究，而四氟甲醚菊酯作为一种常用的拟除虫菊酯类杀虫剂，对其产生抗药性的白纹伊蚊所伴随的适合度代价还尚不清晰^[17]。了解白纹伊蚊对四氟甲醚菊酯产生抗药性机理和伴随的适合度代价对今后抗药性蚊虫治理具有非常重要的意义。

有研究表明溴氰菊酯抗性品系白纹伊蚊发育时间明显延长，成蚊的平均寿命明显缩短^[19]，证明了白纹伊蚊溴氰菊酯抗性的产生增加了其适合度代价。而白纹伊蚊对四氟甲醚菊酯产生抗药性所伴随的适合度代价尚不明确，本研究将实验室筛选的白纹伊蚊四氟甲醚菊酯抗性品系与敏感品系进行对比，结果显示，抗性品系从幼虫发育至成蚊的时间比敏感品系长，但二者的化蛹率和羽化率无显著差异，说明抗性品系明显滞育，但并不会导致它们在生长期的死亡。结合成蚊的生存曲线来看，雌蚊敏感品系平均寿命要显著长于抗性品系，因此四氟甲醚菊酯抗性的产生缩短了白纹伊蚊成蚊的寿命，导致成蚊过早死亡，表明抗药性种群产生了生长发育和生存的适合度代价，这与溴氰菊酯抗性品系白纹伊蚊适合度代价的报道结果相一致。此外本研究还通过分析白纹伊蚊四氟甲醚菊酯抗性品系的种群数量增长，发现抗性品系白纹伊蚊较敏感品系的种群加倍时间长，证明了抗性种群数量增长比敏感种群慢，表明抗药性种群产生了生殖适合度代价。我们推测可能是由于靶标基因位点的改变导致耗能增加，如昆虫的抗药性是由于靶标基因位点敏感性降低引起的，那么昆虫可能投入额外的能量用来维持这种改变后的靶标状态。但是昆虫的生长同样需要大量能量维持细胞分裂、组织构建和器官发育等生理过程。充足的能量供应能促使昆虫快速生长，使其体型增大、体重增加，而当昆虫调动过多能量投入到维持改变后的靶标状态，就会导致能量匮乏，出现昆虫生长缓慢，生长停滞，发育不良甚至生殖障碍的情况^[19]。所以靶标位点改变耗能对昆虫的生存和繁殖会产生多方面的影响。长期使用传统杀虫剂导致许多昆虫产生了抗药性，使得杀虫效果逐渐降低，而通过作用于新的靶标位点使白纹伊蚊体内能量失衡，以此研发新型杀虫剂，为解决日益严重的昆虫抗药性问题提供了新的思路和方法。

本研究通过种群生命表比较法分析了抗药性白纹伊蚊适合度代价，发现抗性白纹伊蚊幼虫期和蛹期发育时间延长、成蚊雌蚊寿命缩短和繁殖力下降，表明白纹伊蚊对四氟甲醚菊酯产生抗性后存在一定适合度代价。适合度代价的研究对于探讨抗性发展规律以及完善抗性治理策略具有重要意义。通过了解白纹伊蚊种群生命适合度代价的产生机制和影响因素，我们可以更有效地预测和管理害虫的抗性进化，从而延长杀虫剂的使用寿命，减少化学农药的使用量，保护生态环境和人类健康。