《Austral Entomology》：Attack of the Clone(s): Genomic Analysis of Russian Wheat Aphid in Australia Suggests a Single Extant Clonal Lineage but Multiple Invasions
编辑推荐：
俄罗斯小麦蚜虫（Diuraphis noxia）于2016年首次在澳大利亚被检测到。尽管该虫被认定为高度适应性害虫，但对其在澳大利亚入侵范围内的遗传变异，尤其是相对于其他入侵地区的情况，仍知之甚少。研究人员利用在澳大利亚谷物种植区采集的D. noxia基因组数
俄罗斯小麦蚜虫（Diuraphis noxia）于2016年首次在澳大利亚被检测到。尽管该虫被认定为高度适应性害虫，但对其在澳大利亚入侵范围内的遗传变异，尤其是相对于其他入侵地区的情况，仍知之甚少。研究人员利用在澳大利亚谷物种植区采集的D. noxia基因组数据，表征其克隆多样性，并与国际生物型进行比较。研究结果支持澳大利亚存在一个单一现存克隆谱系（clonal lineage），该谱系与南非生物型亲缘关系最近。Delta D（ΔD）估计表明，除2016年澳大利亚首次检测后不久采集的单个种群外，澳大利亚种群之间无显著分化。澳大利亚D. noxia与国际生物型之间也检测到显著分化。这些发现提示，现存澳大利亚D. noxia代表一个单一克隆谱系，而另一次入侵可能已经发生但未建立种群。该结果为澳大利亚D. noxia的监测和管理提供了基线，并增进了对入侵遗传学及入侵蚜虫种群进化动态的理解。
研究背景方面，俄罗斯小麦蚜虫（Diuraphis noxia）是一种入侵性谷物害虫，分布于欧洲、亚洲、非洲、北美、南美和澳大利亚等地，通过取食及注入植物毒性唾液成分造成叶片褪绿条纹、卷叶等直接危害，卷叶还可为蚜虫提供躲避接触性杀虫剂和天敌的保护环境。在入侵范围内，该虫的有性生殖环节常缺失或不成功，主要营孤雌生殖（parthenogenesis），但可借显著表型可塑性适应新环境压力并形成生物型（biotype；原文以宿主抗性基因上的差异毒力定义，而非稳定遗传分化）。既往核DNA和线粒体DNA分析显示，遗传变异最高见于西亚和俄罗斯南部有性生殖的本地种群，而孤雌生殖入侵范围内变异有限；美国与南非生物型全基因组测序也发现尽管毒力不同，分化仍有限，唾液蛋白多态性与毒力功能关联不明确，因而其破坏性毒力生物型出现的机制仍存困惑。澳大利亚于2016年首次检测到D. noxia，随后扩散至整个谷物带；当地尚未观察到有性生殖，目前仅识别为单一生物型RWAau1，但尚不清楚当前分布是由单一克隆谱系还是多次不同谱系入侵所致。由于澳大利亚入侵范围内环境异质性明显，明确基因组变异程度有助于判断克隆适应潜力；同时，引入新型兼性内共生菌（facultative endosymbionts）进行生物防治的兴趣增加，而澳大利亚D. noxia不携带兼性内共生菌、仅携带专性内共生菌Buchnera aphidicola，使其成为转感染研究候选，评估克隆多样性对相关策略亦有必要。研究人员据此利用单核苷酸多态性（SNP）位点，对2016—2023年澳大利亚13个地点田间采集的D. noxia进行基因分型，并与已发表国际生物型同一位点数据比较，以刻画澳大利亚及全球范围内的基因组分化。研究结论支持现存澳大利亚种群主要为单一克隆谱系RWAau1，且与南非谱系最近；2016年采集的V02A种群显著分化，可能代表另一未定殖或已消失的入侵谱系。该研究为澳大利亚D. noxia监测、管理和入侵遗传学提供了基线。
研究人员采用的关键技术方法为：2016—2023年从澳大利亚谷物种植区13个地点采集D. noxia，样本来自栽培谷物及非作物禾草；国际比较样本为NCBI SRA中4个南非生物型（SA1、SA5、SAM、SAM2）和1个美国生物型（RWA2）。对澳大利亚种群主要采用DArTseq平台进行SNP基因分型，V02A采用Illumina双端全基因组测序；序列比对至D. noxia参考基因组（GCF_001186385.1），经变异检测与过滤后获得SNP；用polyRAD估计基因型概率，以ΔD统计量（基于Shannon熵、按等位基因频率加权、不假设繁殖方式的遗传分化指标）评估分化，并结合主成分分析（PCA）和邻接树（neighbour-joining tree）分析。
研究结果如下。3.1 测序与过滤：DArT-seq对182个D. noxia个体、80个文库测序，获得127 307 029条reads，平均每文库1 591 338条；其中92 369 759条比对至参考基因组，平均比对率72.56%。对13个个体、4个文库的Illumina双端全基因组测序获得420 270 554条reads，平均每文库105 067 639条；410 333 930条比对成功，平均比对率97.57%。两种测序策略共享覆盖660 869 bp。国际生物型数据112 131 872条reads中92 924 680条比对成功，平均比对率92.02%。全部过滤后仅保留308个SNP位点（占比对覆盖的0.004%），表明所采样D. noxia种群基因组变异有限。3.2 种群与生物型间的基因组分化：ΔD置换检验显示，除2016年6月采集的V02A外，所有澳大利亚种群两两比较均不显著（p>0.05）；V02A与所有其他澳大利亚种群显著不同（p<0.05，ΔD=0.0331–0.0341），包括2019年同地点采集的V02B。排除V02A后，澳大利亚种群ΔD近乎为零（0.0003–0.0021，平均0.0011），与这些种群代表单一克隆谱系RWAau1相一致；数据提示澳大利亚至少出现过两个不同谱系，但是否同时到达尚不清楚。澳大利亚种群与国际生物型比较中，V02A与SA1、SAM、SAM2和RWA2差异不显著（p>0.05），V02A与RWA2比较边缘不显著（p=0.051），且V02A与RWA2的ΔD（0.058）高于V02A与其他澳大利亚显著分化种群比较；SAM与V02B、S07、S22比较亦边缘不显著（p=0.051–0.054），但ΔD较高（0.0645–0.0659）。基于ΔD，这些种群与国际生物型仍分化；国际生物型重复不足可能造成置换检验虚假显著性，应谨慎解释。PCA结果与ΔD一致：仅澳大利亚种群时，前两轴解释58.2%和5.1%遗传变异，除V02A单独沿轴1明显分离外，其余种群紧密聚类；澳大利亚与国际生物型合并时，前两轴解释54.4%和21.1%，V02A与其余澳大利亚种群分离仍存在，SA1、SAM和SAM2形成单一簇，V02A与SA1组之间遗传分离小于其与澳大利亚其他种群，SA5与RWA2相对更近但仍存在变异。邻接树中，5个节点支持率>90%，其中3个达100%；除V02A外所有澳大利亚种群组成100%支持的最根部支系，V02A随后分支至下一个100%支持的衍生支系；国际生物型以96%节点支持分开，SA1、SAM、SAM2组成100%支持支系，SA5与RWA2组成93%支持支系；在100%支持支系内，ΔD为0.00129–0.0103，可作为区分谱系的ΔD上限估计。
讨论部分总结：系统发育分析和除V02A外澳大利亚种群间缺乏基因组分化，支持现存澳大利亚D. noxia为单一克隆谱系RWAau1；分析支持SAM和SAM2为SA1克隆后代，RWAau1与V02A最接近该南非谱系，并提供SA5与RWA2共享更近期祖先的新证据。V02A采集于2016年6月Horsham，距离澳大利亚首次检测地南澳大利亚Tarlee约400 km；其从未建立实验室群体，不能完全排除实验室维持期间遗传漂变影响，但其余种群采样和维持时间差异大却未检测到分化，故研究人员认为V02A可能代表2016年或更早进入澳大利亚的独立克隆谱系，现存RWAau1与V02A是否共同来源或独立引入仍不清楚。南非SA1克隆、RWAau1与V02A之间有限分化提示经南非入侵澳大利亚的可能路径；但既往 ancestry 报道及RWA2测序时间与本研究系统发育存在矛盾。COI等保守标记在孤雌生殖种群中多态性近无，难以解析近期祖先关系；RAPD-PCR、同工酶和微卫星亦受低变异限制，高覆盖全基因组分析可能更有信息。Zhang et al. 曾结合宿主基因组标记与Buchnera基因组比较，解析南非至美国及美洲后续入侵路径，Buchnera基因组亦可为宿主祖先研究提供变异。单一克隆谱系并不排除未来毒力生物型演化；基因组分化本身不能可靠区分D. noxia毒力生物型，SA1与SA5虽共享祖先但明显分化，而生物型可在入侵后于单一克隆内出现；DnDNMT3和DnTET差异表达驱动的甲基化梯度与南非生物型毒力相关，提示表观遗传适应可能重要并可监测。单克隆谱系对控制策略有意义：若存在多个克隆谱系，则需要多重转感染以验证表型效应是否跨遗传背景一致；除非Buchnera变异在澳大利亚单一克隆中分离，否则当前无需担忧，但仍需持续监测。研究结论翻译：总之，该初步基因组筛查提示，现存澳大利亚D. noxia种群以单一克隆谱系RWAau1为主，并可能与南非D. noxia共享近期祖先。然而，结论应结合若干局限解释：采样为机会性，未覆盖D. noxia现已定殖的全部澳大利亚地区，采集宿主限于相对有限的栽培谷物和非作物禾草；仅308个SNP位点保留用于分析，国际比较仅限少量可用生物型。因此，不能排除澳大利亚其他地区存在额外谱系，也不能自信识别入侵来源；需更广泛地理采样、更宽宿主范围以及来自澳大利亚和潜在源种群的更多全基因组数据，才能全面解析D. noxia在澳大利亚的入侵历史。
