淡水蝽 Aphelocheirus aestivalis（昆虫纲：异翅亚目）中沃尔巴克氏体（Wolbachia）感染的多样性与空间分布及其与线粒体DNA变异的关联

《Hydrobiologia》：Diversity and spatial distribution of Wolbachia infection and its association with mitochondrial DNA variation in the freshwater bug Aphelocheirus aestivalis (Insecta: Heteroptera)

　　

沃尔巴克氏体（Wolbachia）是一种广泛分布的经母系遗传的节肢动物内共生菌，可改变宿主的繁殖、种群结构与遗传多样性。然而其在水生昆虫中的流行率与进化效应仍知之甚少。研究人员调查了感染欧洲淡水蝽 Aphelocheirus aestivalis 的沃尔巴克氏体的流行率、空间分布与遗传多样性，并评估其对宿主线粒体与核遗传变异的潜在影响。利用 wsp 标记筛选了63个种群的个体。76%的种群检出感染，个体感染率介于0–100%，平均42.1%。感染频率在生命周期阶段间存在差异，成虫高于若虫，但与性别和地理无关联。系统发育分析在超群A（supergroup A）内揭示两个分歧亚群，分布重叠且在中欧存在共感染。宿主多样性通过线粒体和核标记评估：线粒体数据呈现异质性（heteroplasmy）、低核苷酸多样性与优势单倍型，核变异有限。感染个体线粒体核苷酸多样性降低，但单倍型多样性中等。中性检验指示稀有单倍型过剩，宿主单倍型与沃尔巴克氏体群组无显著关联。结果表明沃尔巴克氏体在 A. aestivalis 中广泛但异质分布，并可能通过种群或选择过程影响线粒体多样性。

研究背景与立项依据

沃尔巴克氏体（Wolbachia）作为一类经卵母细胞母系垂直传播、并可借助水平转移扩散的革兰氏阴性内共生菌，在陆生节肢动物中感染率可达66%，而在水生昆虫中报道偏少但估算可达52%。该菌通过胞质不相容（cytoplasmic incompatibility, CI）、孤雌生殖诱导、雄性致死（male killing）及雌性化等繁殖操纵机制，深刻塑造宿主的核与线粒体遗传结构；母系共传递还可导致线粒体单倍型搭车（selective sweep）及异质性（heteroplasmy）。淡水底栖蝽 Aphelocheirus aestivalis 在欧洲广布但受栖息地不稳定威胁，既有研究仅限波兰西北部小尺度，发现低遗传变异与高沃尔巴克氏体感染，并报道线粒体异质性，但大尺度流行规律与菌-线粒体耦合关系未明。研究人员因此在《Hydrobiologia》发表本工作，首次以欧洲63个种群的大样本，联用 wsp 与多线粒体/核标记，量化沃尔巴克氏体超群A内两个亚群的时空分布，并检验其是否驱动宿主线粒体多样性变化。

主要技术方法概括

研究人员于2008–2022年从欧洲63个河流采样点采集 A. aestivalis 成体与若虫标本（部分样本既往用于异质性研究），腹组织提取基因组DNA后，采用常规PCR扩增沃尔巴克氏体 wsp 片段及宿主线粒体 COI、CYTB、16S rRNA 与核 EF-1α 片段；测序后利用 Geneious Prime 识别异质性位点，Clustal Omega/MAFFT 比对，PHASE 重构二倍体核单倍型，DnaSP 抽提变异，MrBayes 建贝叶斯系统树，PopART 构中位连接网络，Arlequin 计算 Tajima's D、Fu's Fs 与 AMOVA，R 包 ape 中 PACo 与 ParaFitGlobal 检验菌-线粒体协同系统发育信号，QGIS 局部 Moran's I 做空间自相关（kNN, BH-FDR 校正）。

研究结果

沃尔巴克氏体流行率：63个种群中48个（76%）检出，个体感染率42.1%（678/1610），成虫显著高于若虫（Mann–Whitney U=875.0, P=0.01），雌雄无差异，与经纬度无相关；局部 Moran's I（k=8, BH-FDR q<0.05）在西北波兰与匈牙利 Karasica 呈高-高聚集，德国 Ostpeene 为低-高离群，k=10时信号消失，说明聚类具邻域尺度依赖。

沃尔巴克氏体遗传多样性：从感染标本获128条 wsp 序列、25个变异型；比对 PubMLST 显示短型（463 bp, H_W1–14）近缘性低，长型（507–508 bp, H_W15–25）近缘 allele 111/70。贝叶斯树归为超群A下两新亚群 Aesti1（Hd=0.286, π=0.00093）与 Aesti2（Hd=0.768, π=0.00291），前者偏北、后者偏西，中欧13个种群共感染且多>50%感染。

Aphelocheirus aestivalis 遗传多样性：COI（719 bp, 136单倍型, Hd=0.7627, π=0.001518）以 C_N1 占49%；CYTB（361 bp）仅异质性位点区分4单倍型，π=0；16S rRNA（439 bp, 11单倍型, Hd=0.041, π=0.00046）以 S_N1 占92%；EF-1α（731 bp）杂合位IUPAC编码得7变异型（π=0.000029），PHASE 重构等位基因后59变异型（Hd=0.6901, π=0.001832）。串联线粒体贝叶斯树浅分两支，98.8%归主支。

选择性清除与菌-线粒体演化关联：感染组 COI 核苷酸多样性比未感染组降约71%，Tajima's D（?2.5814）与 Fu's Fs（?30.33987）更负，指示稀有单倍型过剩；AMOVA 群内变异为主（F_ST>0.15）；PACo（P=0.110）与 ParaFitGlobal（P=0.205）均不支持 COI 单倍型与 wsp 变异型协同共进化，C_N1 可携任一亚群，H_W1 跨19个宿主单倍型。

讨论与结论翻译

异翅亚目（Heteroptera）中沃尔巴克氏体研究多停留检测层面，水生类群更薄弱；本研究打破既往波兰小样本“固定感染”印象，揭示欧洲尺度下0–100%的种群异质流行，与“菌种在物种间常见、种内未必”假说吻合。成虫偏高感染可能源于若虫菌体载量低致PCR漏检，非必然水平传播差异。超群A两亚群空间部分重叠与共感染提示历史水平转移或复合定殖。宿主线粒体低多样性、优势单倍型与异质性延续既往发现，感染组 π 降低伴中性检验负值符合菌群经CI或搭车引发的选择性清除或瓶颈效应，但菌-线粒体无显著系统发育耦合，说明 Aesti1/Aesti2 未长期锁定特定宿主单倍型。研究人员结论：沃尔巴克氏体在 A. aestivalis 中广泛但异质分布，可能经种群动态或选择过程影响线粒体多样性，而非简单母系共分化；未来需组织特异性定量、多基因座分型与功能表型验证以区分瞬态肠道菌与系统感染。