寄生性蝇类引入后蟋蟀景观中适应性的基因组特征

《SCIENCE ADVANCES》：Genomic signatures of adaptation across a landscape of crickets following the introduction of a parasitoid fly

　　

新型物种互作可为评估遗传适应的最早阶段提供机会。研究人员结合群体基因组学与野外选择调查，探究夏威夷小型、碎片化海岛蟋蟀种群在引入一种以鸣叫雄蟋蟀为寄主目标的寄生性蝇类（parasitoid fly）后的基因组与地理适应景观。多种保护性雄性静音（male-silencing）蟋蟀翅型在该新型选择压力下近期扩散，尽管鸣声在配偶吸引中具有重要作用。研究人员发现蝇类引入后蟋蟀有效种群大小（effective population size，Ne）急剧下降，并识别出由侵染风险介导的空间变异选择区域，这些区域与一个适应性翅型及其他参与抵抗侵染的基因相关。尽管经历反复瓶颈，遗传变异仍以看似受平衡选择（balancing selection）维持的结构变异（structural variant，SV）多态性为主导。该研究揭示了突变选择变化对小种群的广泛后果。蝇介导的选择在适应性雄性静音表型扩散后依然强烈，而这些表型在配偶吸引背景下也会降低雄蟋蟀适合度。

研究背景方面，人类活动通过改变物种分布造成突发选择压力，其中新捕食者或寄生者引入可形成宿主—寄生生物新型共进化互作。夏威夷大洋蟋蟀Teleogryllus oceanicus与引入的窃听型寄生性蝇类Ormia ochracea仅在该地共存，雌蝇利用蟋蟀鸣声进行趋声定位并产内寄生幼虫，导致鸣叫雄蟋蟀高致死率。过去25年出现flatwing、curly-wing、small-wing等雄性静音翅型，虽降低被寄生风险，但鸣声也是雌蟋蟀择偶关键信号，因此这些翅型很少固定。已有研究多聚焦大效应翅型位点，但本地寄生压力是否驱动静音表型频率与分布尚缺全基因组证据。该研究即在此背景下，由研究人员在夏威夷三岛11个地点采样，用300个雄蟋蟀基因组与489次鸣声回放捕蝇试验，分析寄生蝇引入后群体基因组与本地适应格局，论文发表于《SCIENCE ADVANCES》。

主要关键技术方法包括：在夏威夷Kauai、Oahu、Hawaii三岛11个地点采集野生雄蟋蟀后腿样本形成300基因组队列，并结合已发表澳大利亚群体RAD测序数据；以约15×深度全基因组测序进行变异检出与过滤；用Lostruct窗口主成分分析（PCA）和连锁不平衡（LD）识别大片段结构变异（SV）；在无非大SV的染色体5、9、14上用FastTree系统发育、F-branch统计、Phlash基于溯祖的N_e历史和GONE2基于连锁不平衡的近期N_e推断；用鸣声回放吸引雌蝇次数量化本地选择压；用GEMMA做flatwing与curly-wing全基因组关联分析（GWAS），用潜因子混合模型（LFMM）检验等位基因频率与蝇攻击率关系；用DESeq2分析未侵染与侵染4天蟋蟀体组织转录组差异表达，并用topGO做基因本体（GO）富集。

研究结果如下。

结构变异主导遗传与转录组变异：研究人员通过窗口PCA、LD和染色体级PCA发现14条染色体中11条存在大于10 Mb疑似倒位的大SV；与澳大利亚比较显示部分SV为种分布广布型、部分仅存夏威夷。神经组织RNA测序方差分配与差异表达显示SV对局部基因表达有强效应，提示其受平衡选择维持。

寄生性蝇类与蟋蟀近期种群衰退相关：在排除大SV的染色体上做群体结构分析，群体主要按三岛聚类。贝叶斯溯祖与连锁不平衡推断均显示，自澳大利亚种群分流后N_e先降后升，蝇类1989年引入后N_e骤降；Oahu.AC与Hawaii.CL后期采不到蟋蟀，说明雄性静音表型扩散期间瓶颈仍持续。

寄生性蝇类的当代选择存在地点差异：489次10分钟回放中鸣声陷阱引到415只雌蝇、静音陷阱仅1只；地点间攻击率可重复且受降水正影响。研究人员据寿命、侵染致死率等估算减鸣雄蟋蟀相对适合度地点间可变，净相对适合度跨地点方向不一，说明鸣声收益与寄生代价存在空间权衡。

选择影响近期演化的雄性静音翅型频率：跨岛GWAS在Chr1的prospero邻近doublesex区域找到10个两岛共有flatwing关联变异；curly-wing位点位于Chr2的ABDA与UBX同源基因区。LFMM显示40个显著选择关联变异中82.5%靠近curly-wing位点；还检出SIFAR、U520、GBRB3、GLT35等基因，GLT35涉及气管发育。侵染4天转录组有1099个差异表达基因，富集于丝氨酸内肽酶活动、整合素结合与背干气管发育，支持封装反应与气管共优化是寄生互作靶标。

讨论部分总结：研究人员指出寄生蝇引入造成强选择并引发严重N_e下降，小种群快速适应伴随瓶颈、漂变与未来适应力受损；静音翅型虽避寄生却削弱性选择信号，可能放大繁殖输出下降。部分无鸣唱雄蟋蟀地点仍有蝇攻击，可能因O. ochracea对其他共存蟋蟀种鸣声保持趋声性。除静音翅型外，歌唱行为调制与气管发育相关基因也受选择，显示共进化不只走“消去鸣声”单一路线。大SV多态在多次瓶颈小种群中保留，提示平衡选择来源与抑制重组间的适应潜力权衡，是检验结构多态功能的理想系统。

研究结论部分翻译：理解促进遗传适应的条件既是进化生物学目标，也具有现实紧迫性。Ormia ochracea引入夏威夷后对Teleogryllus oceanicus施加新型极端选择，该研究罕见揭示快速遗传适应期间及之后的种群情境。蝇类致N_e锐减，田野试验表明即便静音表型扩散，多数地点当代选择仍强。选择强度空间差异与curly-wing区域等位基因频率及其他共进化相关基因相连。小种群强选择适应必然伴遗传瓶颈，侵蚀后续选择可利用变异并增加有害突变风险。雄性静音表型虽避寄生，但因鸣声是重要性信号，也可能通过降低繁殖输出削弱种群韧性。部分地点虽多年无鸣唱雄蟋蟀仍见非可忽略蝇攻击，可能因蝇对其他共存蟋蟀种鸣声具phonotaxis，使宿主—寄生者共进化复杂化。除male-silencing表型，SIFAR等歌唱行为基因与GLT35等气管发育基因也呈选择信号，封装反应与幼虫气管共优化构成潜在选择靶标。倒位或多态SV占夏威夷蟋蟀大量遗传变异并对表达有强效应，其受平衡选择维持，但是何种平衡机制及抑制重组如何权衡适应潜力仍待解；该系统可检验结构多态的平衡选择来源与适应重要性。