《Ecology and Evolution》：Population Genomics Approaches Identify a Cryptic, Emerging Generalist Pest Complex
编辑推荐：
摘要：关于害虫管理所必需的生物学性状信息，如物种身份、食性和扩散能力，通常需要在实验室和田间对害虫自然史进行费时且高强度的研究。然而，新的农业害虫威胁不断出现，往往需要在信息有限的情况下迅速作出反应。研究人员结合分子肠道内容物分析和RAD-seq技术，研究了美
摘要：关于害虫管理所必需的生物学性状信息，如物种身份、食性和扩散能力，通常需要在实验室和田间对害虫自然史进行费时且高强度的研究。然而，新的农业害虫威胁不断出现，往往需要在信息有限的情况下迅速作出反应。研究人员结合分子肠道内容物分析和RAD-seq技术，研究了美国一种新兴的重要农业害虫——花生掘土蝽（Pangaeus bilineatus Say，半翅目Hemiptera，土蝽科Cydnidae）的物种身份、植物食性组成和种群遗传结构。研究人员发现，两种形态相似的掘土蝽科昆虫（包括P. bilineatus）常在花生田附近的灯光诱捕器中被捕获，其中一种（Dallasiellus lugubris）此前未被视为花生的害虫。分子肠道内容物分析显示，这两种蝽的食性范围广泛但略有差异。令人惊讶的是，花生在两种蝽的食谱中均占比极少。RAD-seq分析揭示了两种蝽均存在与弱距离隔离相符的证据以及适度的空间遗传分化。综上，这些结果表明可能存在一个害虫复合体，而此前仅关注其中一种物种。此外，它们广泛的食性和受限的空间种群扩散模式也可能解释了该潜在害虫复合体危害的零星发生特征。对新兴害虫挑战的应对措施可受益于种群基因组学技术所提供的见解，这将为研究开辟新途径，并支持快速制定针对新型害虫的管理策略。
论文解读：种群基因组学揭示花生掘土蝽的隐秘害虫复合体
一、研究背景与科学问题
农业生态系统中，新害虫的持续出现对粮食生产构成重大威胁。有效应对新害虫往往需要快速获取其物种身份、寄主范围、扩散能力等基础生物学信息，然而传统研究方法如田间观察、温室试验和标记-回捕法等通常耗时费力，难以满足应急管理需求。在此背景下，利用分子肠道内容物分析和简化基因组测序（如RAD-seq）等现代分子工具，可在较短时间内评估害虫的食性宽度和种群结构，为制定管理策略提供关键信息。
花生掘土蝽（Pangaeus bilineatus Say）是美国东南部新近兴起的重要花生病虫，其危害具有零星爆发的特点，给防治带来极大困难。尽管该虫被长期认定为花生田的优势土蝽物种，但关于其寄主专一性、田间扩散格局以及是否存在同域共存的其他近缘物种等问题尚不清楚。尤其值得关注的是，花生田中还捕获了形态相似的另一种土蝽Dallasiellus lugubris，该物种的生物学信息极为匮乏，是否取食花生亦属未知。因此，研究人员以美国佐治亚州南部花生主产区为研究区域，运用分子肠道内容物分析和RAD-seq技术，系统研究了这两种土蝽的物种鉴定、植物食性组成及种群遗传结构，旨在阐明实际害虫复合体的组成及其生态学机制。
二、主要技术方法
研究样本采自佐治亚州兰道夫县（Randolph County）和沃斯县（Worth County）花生田边的灯光诱捕器，并补充了2019年从伊曼纽尔县（Emanuel County）、杰斐逊县（Jefferson County）和沃斯县手采的个体，共获得127头P. bilineatus和123头D. lugubris。提取全基因组DNA后，一部分用于RAD-seq测序（nextRAD方法，Illumina NovaSeq平台），另一部分用于扩增叶绿体trnF基因并进行PacBio三代测序，以解析肠道内容物中的植物DNA组成。RAD-seq数据经Stacks 2软件denovo组装鉴定SNP，分析种群遗传结构；肠道内容物数据经DADA2降噪、基于自建质体数据库进行物种注释，并通过多样性和多元统计方法比较两种蝽的食性差异。此外，利用线粒体COI基因序列构建系统发育树验证物种鉴定结果。
三、研究结果
（一）样本鉴定与物种组成
分子鉴定结果显示，在118个成功鉴定的样本中，67个为D. lugubris，34个为P. bilineatus，另有17个无法确定物种归属。COI条形码分析进一步验证了鉴定结果。时间动态上，P. bilineatus主要在夏季出现（2020年6月达高峰），而D. lugubris主要在秋季出现（2021年9月达高峰），两者的发生高峰期存在明显错位。
（二）食性分析
研究人员从160个样本中获得120,041条原始reads，经质量控制后保留138个田间样本和10个人工饲养样本。共鉴定出188个植物科和186个植物属。两种蝽均表现出极广的杂食性，各自食谱中均含有43–50个植物属。常见植物包括繁缕属（Stellaria）、萝卜属（Raphanus）、雀稗属（Paspalum）、野芝麻属（Lamium）、山核桃属（Carya）和松属（Pinus）等野生及木本植物，也检测到芭蕉属（Musa）、菊苣属（Cichorium）、芝麻属（Sesamum）等经济作物。然而，花生（Arachis）的检出率极低，仅在4头D. lugubris和2头P. bilineatus个体的肠道中检测到少量花生DNA。相比之下，10头人工饲喂花生的P. bilineatus均检测到大量花生DNA（占全部检出序列的99.78%以上），验证了检测方法的可靠性。
食性多样性分析显示，两种蝽的食谱组成无显著差异（F5,60=1.099, R2=0.085, p=0.330），但采样月份对食谱组成有边缘显著影响（F7,60=1.415, R2=0.154, p=0.045），这可能与两种蝽发生高峰期的植物物候差异有关。进一步仅对比8月份同时捕获的两种蝽，仍无显著食性差异（Bray-Curtis: F1,19=1.110, R2=0.06, p=0.335；Jaccard: F1,19=1.204, R2=0.06, p=0.222）。未发现任何一种植物在所有个体中稳定出现，繁缕属的检出率最高，也仅在12头D. lugubris中的5头中检出。
（三）种群遗传结构
RAD-seq分析在过滤后保留P. bilineatus 5774个SNP位点（113个样本）和D. lugubris 8845个SNP位点（105个样本）。主成分分析（PCA）显示，P. bilineatus的样本沿PC1（解释1.5%变异）和PC2（解释1.4%变异）有轻微的地理分离；D. lugubris的样本则没有表现出按地点的分离，但沿PC3和PC4按年份有轻微分离。空间主成分分析（sPCA）显示两种蝽均存在显著的全局空间遗传结构（P. bilineatus: 标准化统计量=5.81, p=0.005；D. lugubris: 标准化统计量=5.21, p=0.005），但无显著的局部空间结构，表明研究区域内存在弱但明显的空间遗传分化。成对FST值方面，P. bilineatus范围为0.017–0.15，D. lugubris范围为0.043–0.119。Mantel检验表明两种蝽均符合弱距离隔离模式（P. bilineatus: Mantel r=0.33, p=0.034；D. lugubris: Mantel r=0.41, p=0.038），即地理距离越远的种群遗传分化越大。P. bilineatus的FST和PCA结果略高于D. lugubris，提示后者可能在农田景观中具有更强的扩散能力。
四、结论与讨论
本研究的核心发现是：在花生田灯光诱捕中，D. lugubris与P. bilineatus同域共存且数量相当，并确认D. lugubris也能取食花生，尽管两者对花生的取食频率均很低。这表明花生土蝽的经济危害可能并非由单一物种造成，而是由至少两个形态相似的物种组成的害虫复合体共同导致。研究结果也提示，以往仅针对P. bilineatus的防治策略可能需要调整为覆盖整个复合体的综合治理方案。
两种蝽广谱的杂食性和忠实的年度田间滞留性，解释了为何花生受害具有零星爆发的特点——当野生寄主或替代作物不足时，它们可能转向取食花生，从而造成偶发性经济损失。此外，两种蝽均表现出的弱距离隔离模式意味着成虫主要通过夜间飞行在农田景观中进行有限距离的扩散，种群间基因流较为有限，这可能影响抗性基因的空间传播，也提示局部管理措施可能更为有效。
研究同时存在一定局限性：由于测序的是叶绿体基因trnF，对花粉污染不敏感，但松树和禾本科植物的花粉可携带叶绿体DNA，因此对这些类群的检出需进一步实验验证；样本仅来自灯光诱捕的成虫，无法反映若虫阶段（可能正是花生荚果期危害的主要虫态）的食性；某些植物属（如芭蕉属、芝麻属）的检出可能受到样品运输和加工过程的环境污染影响。未来研究应结合室内饲养实验、若虫期食性分析和景观尺度种群基因组学，进一步明确该害虫复合体中各物种的生态位分化及对花生产量的实际影响。
本研究展示了分子肠道内容物分析与RAD-seq联用在快速评估新兴农业害虫食性宽度和种群结构方面的巨大潜力，为害虫复合体的识别和针对性管理策略的制定提供了有力支撑，相关方法和结论对其它新兴害虫的研究与防控具有重要参考价值。
