疟疾，尤其是由恶性疟原虫（Plasmodium falciparum）引起的恶性疟，至今仍是全球最重要的公共卫生挑战之一。在包括东南亚在内的许多地区，疟疾的传播强度已显著降低，进入了消除阶段。然而，在通向消除的道路上，新的科学问题也随之浮现：当我们采取极为强力的控制措施——例如在社区广泛设立“疟疾哨点”（Malaria Posts, MPs）以实现快速诊断和治疗，并在高传播“热点”村庄实施“大规模药物投放”（Mass Drug Administration, MDA）——来近乎“清剿”寄生虫时，会对疟原虫种群本身产生什么样的深层影响？这些措施是会像我们希望的那样，安全、有效地将寄生虫推向灭绝边缘，还是可能“按下葫芦浮起瓢”，意外地加速抗药性寄生虫的筛选和扩散，为最终消除埋下隐患？
为了回答这些至关重要的问题，一个国际研究团队在缅甸东部的克耶邦（Kayin State）展开了一项为期近5年（2015-2020年）的大规模基因组流行病学研究。该地区自2014年起实施了由疟疾消除特遣部队（Malaria Elimination Task Force, METF）主导的消除计划，核心策略正是上述MPs网络和目标性MDA。这项雄心勃勃的计划取得了显著的流行病学效果，将恶性疟的发病率降低了97%。但在此背景下，寄生虫种群在基因组层面经历了怎样的演变，却是一个未知的“黑箱”。研究人员意识到，解开这个黑箱，不仅能够从分子层面验证控制措施的效果，更能揭示在极低传播强度下寄生虫种群动态的基本规律，以及抗药性进化的新模式，从而为全球范围内的疟疾消除行动提供关键的决策依据。他们的研究成果最终发表在了顶级期刊《自然·微生物学》（Nature Microbiology）上。
为了开展这项研究，作者团队运用了几个关键的技术方法。首先，他们建立了大规模的样本队列，从克耶邦283个疟疾哨点收集了2270份恶性疟原虫感染者的干血斑样本，时间跨度达58个月。其次，他们利用选择性全基因组扩增结合高通量测序技术，成功获得了这些临床样本的全基因组数据，并经过严格质控，最终基于1927个高质量样本的25461个单核苷酸多态性进行分析。最后，研究采用了先进的群体基因组学分析手段，包括基于身份同源（Identity by Descent, IBD）的遗传相关性分析、有效种群大小估计、时空相关性分析等，以多维度解析寄生虫的种群结构、传播动态和演化历程。
极端克隆扩张与近交
研究表明，克耶邦的疟原虫种群经历了极端的克隆扩张和严重的近交。在1781个单基因型感染样本中，仅鉴定出166个独特的基因组，遗传丰富度极低。研究发现了93个IBD簇（指基因组相似度≥90%的样本组），其中14个为大型克隆簇。前十大的克隆簇就占据了总样本量的51.6%，表明种群多样性严重匮乏。更有甚者，在166个独特基因组中，有110个（占感染样本的66.9%）归属于9个密切相关的“家族”，成员间遗传相关性≥0.45，其中两个家族甚至能鉴定出亲本和子代，这直接证明了在极小的有效种群中，近交和克隆传播是主要繁殖模式。
局部化传播与单倍型的区域稳定性
寄生虫的传播表现出强烈的空间局限性。不同的IBD克隆簇在地理分布上呈现清晰的聚集性。例如，携带kelch13-R561H的最大克隆簇主要分布在Hpapun镇北部，而携带F446I和P441L的克隆簇则分别集中在其中部和西部。空间相关图分析证实，在≤20公里范围内，寄生虫间的遗传相关性呈正相关，而在27.5-90公里范围内则呈负相关，这清晰地划定了不同遗传簇或传播灶的边界。同时，许多IBD克隆簇的传播时间可持续长达数年，表明在低传播环境下，独特的寄生虫基因型可以长期、稳定地在局部地区循环。
东南亚地区寄生虫种群的长距离连通性
通过比较东南亚东西部地区的寄生虫种群，研究发现西部东南亚（包括克耶邦、泰缅边境诊所及其他缅甸地区）和东部东南亚（柬埔寨、越南、老挝）的疟原虫形成了两个遗传结构清晰的亚群。没有证据表明在2015-2020年间，东西部种群之间存在克隆传播或近期基因交流。这一结论得到了耐药性标记分布的支持：在东部广泛存在的、与哌喹耐药相关的pfcrt基因突变，在克耶邦及西部泰国样本中完全缺失。
评估寄生虫种群大小的基因组学指标
研究人员评估了多种遗传指标在监测传播强度变化中的应用。在克耶邦，多基因型感染的比例一直维持在较低水平（约7.6%），且随时间变化不显著。遗传丰富度同样维持在很低水平（0.10-0.30）。然而，基于连锁不平衡计算的有效种群大小在克耶邦和泰缅边境诊所的低传播时期都显示出显著的下降趋势。这表明，在传播强度极低（接近消除）时，有效种群大小可能是一个比多基因型感染比例或遗传丰富度更敏感的监测指标。
kelch13等位基因的演化
在青蒿素抗性方面，克耶邦的优势等位基因与邻近的泰缅边境诊所截然不同。后者以C580Y为主导，而克耶邦则以P441L、F446I、R561H和G449A为主，C580Y频率很低（2.28%）。值得注意的是，在2016至2019年间，主要抗性等位基因的频率保持相对稳定。然而，到了2020年，携带kelch13-R561H的IBD克隆簇1的频率急剧上升至74.2%，这并非源于强药物选择下的选择性清除，而是由于疟疾在其他地区被基本消除，唯独在该克隆簇占主导的北部Hpapun地区仍有残留，从而在种群瓶颈和遗传漂变的作用下，其相对频率出现快速、剧烈的变化。
大规模药物投放对寄生虫种群结构的影响
研究为MDA的有效性提供了基因组学证据。在三个实施了MDA的区域内，比较干预前6个月和干预后6个月收集的寄生虫，发现它们的遗传相关性显著低于未实施MDA的对照区域在相同时期内的寄生虫相关性。这表明MDA有效地清除了当地的寄生虫库，之后重新定植的寄生虫与之前的种群遗传关系较远，产生了“奠基者效应”。全基因组扫描未发现MDA导致明显的耐药性选择信号。
结论与意义
这项研究系统揭示了在接近消除的背景下，高强度疟疾控制措施对恶性疟原虫种群基因组产生的深刻影响。主要结论包括：强化控制（MPs联合目标性MDA）能显著降低寄生虫的有效种群大小，增强遗传漂变的作用；导致寄生虫种群呈现高度克隆化、近交化和强烈的空间遗传结构；在极小的种群中，遗传漂变而非正选择，可能成为驱动抗性等位基因频率发生快速、剧烈变化的主要力量，如2020年观察到的kelch13-R561H克隆株的频率爆发；MDA能有效切断局部传播链，产生可检测的基因组瓶颈效应；东西部东南亚的疟原虫种群间存在显著的遗传隔离，限制了多重耐药菌株的跨区域传播。
这些发现具有多重重要意义。首先，它从基因组层面证实了以社区为基础的综合消除策略（特别是MDA）能够有效清除寄生虫库，为这类干预措施的有效性提供了分子证据。其次，研究强调了在消除后期，监测重点需要从传统的流行病学指标，部分转向能反映种群规模和遗传结构的基因组指标（如有效种群大小）。再者，研究揭示了一种新的抗性演化模式：在接近消除的微小种群中，遗传漂变可能超越药物选择，成为等位基因频率动态的主导因素，这对预测和应对消除过程中的“最后阶段”疫情具有警示和指导价值。最后，研究明确了东南亚地区疟原虫种群的地理隔离现状，这为分区制定和调整消除与抗药性遏制策略提供了科学依据。总之，这项工作将基因组学深度融入消除科学，为我们理解在人类强力干预下病原体的演化轨迹，以及最终实现疟疾消除的目标，提供了至关重要的新见解。
打赏
