乌干达恶性疟原虫临床分离株的全基因组测序
系统鉴定与青蒿素部分抗性和苯芴醇敏感性下降相关的近期方向性选择信号。
对190份低复杂性感染样本进行选择性全基因组扩增和Illumina测序，158份样本通过质量控制，最终157份用于下游分析。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
与乌干达疟原虫对青蒿素类联合疗法药物敏感性下降相关的恶性疟原虫PX1多态性的出现与传播
Emergence and spread of Plasmodium falciparum PX1 polymorphisms associated with decreased susceptibility to antimalarials in Uganda
实验完整度
高全基因组测序、离体药物敏感性分析及px1敲除虫株体外功能验证等多个独立证据层级。
主要模型
乌干达恶性疟原虫临床分离株 3D7恶性疟原虫株 PX1-KO-1和PX1-KO-2敲除虫株 PX1-GFP整合对照虫株
重点核对
全基因组测序样本纳入标准：低复杂性感染、按年份和地点匹配的不同K13基因型或药物敏感性表型 离体药物敏感性分析：仅纳入测序深度≥25×且等位基因支持率≥90%的样本 PIN单倍型定义：L1222P、M1701I、D1705N突变以及delV1680–N1685和delN811–Y822缺失 体外72小时药物敏感性测定：三种生物学重复，每种药物设稀释系列 混合效应模型协变量：采样点、寄生虫血症、复杂性感染、采样年份和K13基因型
摘要
青蒿素类联合疗法是疟疾治疗和控制的基石。在非洲，蒿甲醚-苯芴醇是使用最广泛的一线青蒿素类联合疗法，但其在乌干达的疗效日益受到青蒿素部分抗性出现和苯芴醇敏感性降低的威胁。为识别导致这种敏感性下降的位点，本研究评估了来自乌干达的157份恶性疟原虫全基因组序列中的选择信号。尽管在Kelch13 C469Y和A675V突变周围观察到延长的单倍型，但近期选择的最强信号集中在编码磷酸肌醇结合蛋白（PX1，PF3D7_0720700）的7号染色体片段上。一个由三个PX1突变（L1222P、M1701I和D1705N）和两个缺失（命名为PIN）代表的单倍型于2008年首次出现并迅速增加，到2016年在乌干达北部达到>50%的流行率，到2023年在乌干达东部达到>50%。携带PIN的虫株对苯芴醇、甲氟喹和双氢青蒿素（蒿甲醚的活性代谢物）的离体敏感性显著降低。在体外破坏px1的虫株对这三种药物的敏感性增加。因此，PX1多态性似乎影响了非洲疟原虫对关键药物的敏感性。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 乌干达恶性疟原虫中对青蒿素部分抗性和苯芴醇敏感性下降的基因组选择信号及遗传决定因素是什么？
- 核心机制
- px1基因的PIN单倍型（L1222P、M1701I、D1705N及两个缺失）受到近期方向性选择，与恶性疟原虫对苯芴醇、甲氟喹和双氢青蒿素敏感性下降相关；px1敲除导致对这些药物的敏感性增加，提示PX1调控药物反应。
- 主要证据
- 157份乌干达恶性疟原虫全基因组测序显示px1区域最强的iR/iHS选择信号；PIN单倍型携带虫株的离体IC50显著升高；px1敲除虫株在72小时体外药物敏感性试验中对苯芴醇、甲氟喹和DHA的IC50显著降低。
- 研究意义
- PX1多态性可能影响非洲疟原虫对关键抗疟药的敏感性，PIN单倍型的快速传播提示有必要将其纳入分子监测panel并开展临床研究以评估其对治疗效果的影响。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
基因组选择信号扫描
识别全基因组范围内近期正选择的候选区域。
使用isoRelate iR统计量和rehh iHS分析扫描全基因组，并对样本按地理区域和K13基因型分层分析。
px1候选区域单倍型解析与PIN定义
确定候选选择区域内的共享单倍型及其关键变异。
对候选区域进行单倍型可视化和聚类，定义PIN单倍型为L1222P、M1701I、D1705N及两个框内缺失。
PIN单倍型时空流行率分析
评估PIN单倍型在乌干达随时间及地理区域的流行率变化。
对2004、2008、2012年东部样本及2016–2024年东部和北部样本进行ONT扩增子测序和单倍型分型。
PIN单倍型与离体药物敏感性关联分析
评估px1单倍型与恶性疟原虫离体药物敏感性之间的关系。
对465份测序深度≥25×且等位基因支持率≥90%的样本进行离体IC50和RSA分析，并按K13基因型分层。
px1敲除虫株的体外功能验证
验证PX1对疟原虫药物敏感性的功能贡献。
使用两个独立来源的3D7 px1敲除克隆（PX1-KO-1、PX1-KO-2）及PX1-GFP对照，在标准化72小时体外药物敏感性试验中检测多种药物IC50。
全球PX1等位基因分布分析
评估PIN单倍型在全球恶性疟原虫群体中的历史分布。
重新分析MalariaGEN Pf6数据集（6388份样本，2001–2015年）以检测px1单倍型频率和地理分布。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 全基因组测序文库制备 | Watchmaker DNA文库构建试剂盒（含片段化） | Watchmaker Genomics | -- |
| px1长读长测序 | GoTaq预混液 | Promega | -- |
| ONT测序文库制备 | 连接法测序试剂盒SQK-LSK114 | Oxford Nanopore Technologies | SQK-LSK114 |
| 体外药物敏感性试验 | SYBR Green | Invitrogen | -- |
| WR99210 | Jacobus Pharmaceuticals | -- | |
| 细胞培养 | Albumax II | Gibco | -- |
| 含L-谷氨酰胺的RPMI 1640 | Gibco | -- | |
| 庆大霉素 | Gibco | -- | |
| 次黄嘌呤 | Calbiochem | -- | |
| HEPES | Corning | -- | |
| qPCR板读数 | Flexstation 3 | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 乌干达临床分离株全基因组测序样本选择 | 样本来源：2016–2024年乌干达北部和东部卫生机构监测样本；纳入标准：低复杂性感染，并按年份和地点匹配K13突变型与野生型或不同药物敏感性表型；测序深度≥50× 阅读提示：Methods中Study design and sample collection和Library preparation, whole-genome sequencing and variant calling小节；Supplementary Table 1和2 |
| 基因组选择信号分析参数 | iR分析使用MAF≥2%的SNP（n=157）；iHS分析使用单基因组样本（n=118）并排除Tajima's D>1或Pf6中MAF≥5%的SNP；FDR校正；IBD片段长度≥50 kb且至少20个SNP支持 阅读提示：Methods中Selection analysis小节；Fig. 1和Extended Data Fig. 1–4 |
| px1 PIN单倍型定义和分型 | PIN单倍型：L1222P、M1701I、D1705N及delV1680–N1685和delN811–Y822；ONT扩增区域894851–897457；测序深度≥25×；突变支持率≥50% 阅读提示：Methods中px1 genotyping using Oxford Nanopore Technologies小节；Supplementary Table 4和5，Supplementary Fig. 5 |
| 离体药物敏感性分析样本纳入 | 测序深度≥25×；等位基因支持率≥90%；离体IC50通过72小时生长抑制试验测定；RSA使用700 nM DHA处理6小时；混合效应模型控制采样点、寄生虫血症、复杂性感染、采样年份和K13基因型 阅读提示：Methods中Ex vivo drug assays和Statistical analysis小节；Fig. 5a和Extended Data Table 1 |
| px1敲除虫株体外药物敏感性验证 | 虫株：3D7亲本、PX1-GFP对照、PX1-KO-1和PX1-KO-2；药物：氯喹、环氯胍、甲氟喹、苯芴醇、DHA；72小时试验，三次生物学重复；单因素方差分析加Tukey多重比较 阅读提示：Methods中In vitro 72-hour drug susceptibility assay of genetically manipulated parasites小节；Fig. 5b |