No evidence of immunosurveillance in mutation-hotspot-driven clonal hematopoiesis
免疫系统真的在清场吗？克隆性造血未见免疫监视
克隆性造血、免疫监视、MHC结合、体细胞突变、UK Biobank
免疫监视理论认为，携带体细胞突变的细胞会产生正常细胞中不存在的突变肽，这些突变肽经主要组织相容性复合体，即MHC或HLA分子呈递后，可被T细胞识别并清除，从而阻止具有潜在致癌性的细胞克隆持续扩增。然而，这一理论的证据主要来自已经形成的恶性肿瘤，免疫系统是否会在克隆性造血等癌前阶段便对突变细胞施加负向选择，尚缺少大规模人群证据。
这篇文章通过整合UK Biobank的大规模全外显子组测序数据、个体MHC基因型及突变肽-MHC结合预测，探索了免疫识别能力是否影响克隆性造血驱动突变的发生概率和克隆扩增程度，这对于理解免疫系统从何时开始介入肿瘤演化、判断癌前克隆是否受到免疫选择以及开发针对早期血液肿瘤的免疫预防策略具有重要意义。研究者提出，如果某个个体携带的MHC分子能够高效结合并呈递某一突变肽，那么携带该突变的造血干细胞克隆理论上更容易被T细胞识别，其在人群中的出现率应当更低，或者即使出现，其克隆规模也应更小。
研究分析了40个发生于11个常见克隆性造血驱动基因中的热点错义或截短突变，包括DNMT3A、JAK2、SRSF2、SF3B1、IDH1、IDH2、TP53、KRAS和NRAS等。研究人员为每位参与者计算其MHC-I和MHC-II分子与对应突变肽之间的预测结合能力，并将该指标与突变阳性率、变异等位基因频率、年龄、MHC杂合度和HLA等位基因多样性进行比较。
结果显示，无论按照相对结合强弱还是绝对结合阈值划分人群，MHC对突变肽的预测结合能力均未显著影响任何热点突变所驱动的克隆性造血患病率。携带较强结合型MHC的个体并未表现出更低的克隆性造血风险，其克隆大小也未显著小于弱结合或不结合人群。即使针对已有实验支持、被认为具有免疫原性的KRAS G12D、IDH1 R132H、IDH2 R140Q、TP53 R175H和JAK2 V617F等突变，研究者也没有观察到相应呈递型HLA等位基因在突变携带者中明显减少。
因此，文章并不是证明免疫监视绝对不存在，而是表明：在UK Biobank能够检测到的、变异等位基因频率通常超过2%至3%的克隆性造血中，MHC基因型并未对哪些热点突变能够成功扩增产生可检测的强烈选择作用。该结论挑战了将肿瘤阶段免疫编辑机制直接外推至癌前克隆扩增阶段的做法，也提示克隆性造血可能通过低免疫原性、骨髓免疫抑制、MHC表达下降、早期免疫耐受或长期抗原暴露导致的T细胞耗竭等方式逃避免疫清除。
研究以UK Biobank约45万名中老年参与者的全外显子组测序数据为基础。经过MHC基因型质量控制后，MHC-I分析纳入约384,600人，MHC-II分析纳入约327,961人。研究者从外显子组CRAM文件中重新检测26个热点基因组位置上的常见克隆性造血突变，并补充TP53体细胞变异分析，最终保留40个在至少15名参与者中出现的热点突变，共识别出9,309名至少携带一种目标突变的个体。
MHC基因型通过HLA*IMP:02算法进行两字段水平的插补。研究利用NetMHCpan-4.1预测MHC-I结合，利用NetMHCIIpan-4.3预测MHC-II结合。对每个“个体—突变”组合，生成覆盖突变位点的候选肽段，并在该个体全部MHC等位基因中选取最强结合肽的最低洗脱百分位排名作为综合结合分数。
统计分析包括：比较强结合组与弱结合组的突变患病率；比较突变阳性者与阴性者的结合分数分布；按照绝对结合阈值比较结合者与不结合者；使用外周血变异等位基因频率评估克隆大小；分析年龄分层、MHC杂合度及等位基因数量；并针对具有实验性HLA限制证据的突变比较特定HLA等位基因频率。研究还通过克隆动力学模拟评估统计效能，以确定能够被检出的最小选择压力。
首先，在40个热点突变中，将人群按照预测MHC-I结合能力划分为强结合和弱结合两组后，任何突变的克隆性造血阳性率均未出现显著差异。突变阳性者的结合分数也没有系统性低于突变阴性者，说明成功扩增的克隆并没有优先出现在免疫呈递能力较弱的个体中。该分析根据突变频率和适应度差异，具有检测约0.5%至3.0%负向选择压力的统计能力。
其次，年龄分层分析没有发现强结合型个体在任何年龄段具有更低的克隆性造血患病率。MHC等位基因数量和杂合度增加也未与克隆性造血风险下降相关，说明更高的抗原呈递多样性并未转化为可检测的癌前克隆保护效应。
再次，以变异等位基因频率表示克隆规模时，无论将所有突变合并分析，还是分别分析DNMT3A R882C、SRSF2 P95H、IDH2 R140Q等具体突变，强结合组与弱结合组的克隆大小分布均无显著差异。根据突变携带者数量，该分析可检测约1%至4%的克隆适应度差异。采用绝对结合与完全不结合的分组方式后，结果仍然一致。
最后，MHC-II结合预测以及针对特定HLA限制性突变的分析也未支持免疫选择。整体而言，研究未发现MHC基因型能够决定哪些热点驱动突变在血液中成功扩增，也未发现克隆达到一定规模后才触发免疫清除的阈值效应。
样本规模和统计效能突出。 研究利用数十万人的基因组数据和9,000余名突变携带者，对较弱的免疫选择效应也具有较高检出能力。
将“是否发生”与“扩增多大”同时纳入。 研究不仅比较克隆性造血患病率，还利用变异等位基因频率检验免疫系统是否限制克隆后续扩张。
多层次验证阴性结论。 相对结合、绝对结合、MHC-I、MHC-II、年龄分层、HLA多样性及已知免疫原性突变均得到一致结果，使“未发现证据”的结论较为稳健。
纠正潜在聚合偏差。 文章指出，将不同突变直接合并可能产生由高频突变驱动的虚假免疫选择信号，对同类免疫遗传学研究具有方法学警示意义。
MHC结合由计算模型预测，强结合并不等同于突变肽被加工、呈递并诱导有效T细胞反应，无法完全替代免疫肽组学和T细胞功能实验。
外显子组测序主要能够可靠检测变异等位基因频率约2%至3%以上的克隆，无法排除免疫监视只发生于更小、更早期克隆的可能。
研究主要关注40个常见热点突变，结论未必适用于CALR移码突变等免疫原性更强的突变、非热点突变或融合新抗原。
克隆可能通过降低MHC表达、形成骨髓免疫抑制环境、诱导T细胞耗竭或在免疫耐受形成前获得突变而逃逸，这些机制无法仅凭基因型和外显子组数据排除。
Spatial transcriptomic analyses highlight distinct erythroid niches in mice and humans
红细胞岛的物种分歧：人类造血不靠中央巨噬细胞
空间转录组、红细胞岛、巨噬细胞、C1q、ICAM4
红细胞终末分化传统上被认为发生于红细胞母细胞岛中，其经典结构是一个中央巨噬细胞周围环绕多个不同成熟阶段的红系祖细胞。中央巨噬细胞被认为能够提供铁、细胞因子和黏附支持，并吞噬红细胞脱核后形成的核质小体。然而，这一经典模型主要来自小鼠研究和体外二维培养，红细胞岛在完整组织中的真实空间结构，尤其是人类红细胞岛是否同样依赖中央巨噬细胞，一直没有得到系统验证。
这篇文章通过整合Visium和Xenium空间转录组、公开单细胞RNA测序数据、组织学成像、流式细胞术、基因敲除小鼠、人源骨髓类器官和功能阻断实验，探索了小鼠与人类红系造血生态位的空间构成及调控机制，这对于正确解释动物模型、认识人类红细胞生成障碍以及寻找骨髓增生异常综合征等疾病的微环境治疗靶点具有重要意义。
研究首先在小鼠胚胎肝、出生后骨髓和成年造血组织中发现，表达补体成分C1q的巨噬细胞与红系祖细胞和成熟中红细胞高度共定位。C1q不仅可作为小鼠红细胞岛中央巨噬细胞的标志，还参与识别和吞噬红细胞脱核产生的核质小体。值得注意的是，C1q缺失并未完全破坏红细胞岛结构，说明其主要作用更偏向核质清除和稳态维持，而非红细胞岛组装所必需的结构支架。
与此形成鲜明对比的是，在人胚胎肝和成人骨髓中，红系细胞主要与其他红系祖细胞聚集，而不是围绕C1Q阳性巨噬细胞排列。人类红细胞岛因而更多表现为巨噬细胞非依赖性红系细胞簇。在人胚胎肝中，一个约80微米范围内平均可包含约50个红系细胞，明显不同于经典的小鼠中央巨噬细胞模式。研究进一步确定，红细胞特异性黏附分子ICAM4是维持人类红系细胞簇完整性的关键分子。阻断ICAM4会减少红系细胞数量、缩小细胞簇规模并抑制红系集落形成，而阻断RHAG或αv整合素则没有相同效应。
在骨髓增生异常综合征患者中，巨噬细胞非依赖性红细胞岛结构遭到破坏，红细胞之间的空间相关性和ICAM4表达下降。接受低甲基化药物等治疗后，部分患者的红系细胞簇和空间组织关系得到恢复。研究由此提出，人类红细胞生成并非简单复制小鼠模式，而是依赖红系细胞之间的黏附网络。该发现重新定义了红细胞岛的基本概念，也提醒研究者不能将小鼠中央巨噬细胞机制未经验证地直接外推至人类。
研究在多个发育阶段和生理状态下采集小鼠胚胎肝、出生后骨髓及成年造血组织，并分析12周和16周人胚胎肝、健康成人骨髓凝块切片以及骨髓增生异常综合征患者治疗前后的骨髓样本。为避免人骨髓脱钙造成RNA降解，研究采用临床常规制备且无需脱钙的骨髓凝块切片。
研究首先使用10x Genomics Visium V1进行全转录组空间测序，通过Leiden聚类和标志基因识别不同空间区域。由于每个Visium点位包含约10至20个造血细胞，研究进一步整合公开单细胞RNA测序数据，对各点位进行细胞类型富集和去卷积，并以Pearson相关、细胞间距离及邻域分析判断细胞共定位关系。
随后使用Xenium原位空间转录组达到单细胞和亚细胞分辨率，验证C1q阳性巨噬细胞、红系祖细胞和成熟红系细胞的直接空间关系。功能实验包括C1qa基因敲除及骨髓移植、苯肼诱导的应激性造血、巨噬细胞与红系细胞共培养、核质小体吞噬检测，以及外源C1q或C1qa过表达实验。
在人类体系中，研究利用配体—受体分析筛选红系细胞之间的黏附分子，并在人诱导多能干细胞来源的骨髓类器官和CD34阳性造血干祖细胞红系培养体系中，分别加入ICAM4、RHAG和αv整合素阻断抗体，通过流式细胞术、共聚焦成像和集落形成实验评价红系分化及细胞簇完整性。
在小鼠胚胎肝中，Visium分析显示红系细胞、红系祖细胞与C1q阳性巨噬细胞存在显著空间相关性；Xenium单细胞分析进一步证实，C1q阳性巨噬细胞位于红系细胞簇中央。典型小鼠红细胞岛直径约80微米，每个C1q阳性巨噬细胞周围平均分布约25个红系细胞。相同空间关系在出生后骨髓和应激性造血条件下仍然存在。
C1q阳性巨噬细胞同时高表达Hmox1、Slc40a1、Dnase2a、Timd4和Mertk等与铁代谢、凋亡细胞识别和核质降解有关的基因。功能实验表明，C1q可促进巨噬细胞吞噬红细胞脱核后的核质小体，并改善核质清除；但C1qa缺失并未完全阻止红细胞岛形成，说明C1q主要维持清除功能，而非决定红细胞岛是否组装。
在人类16周胚胎肝中，红系标志基因区域与巨噬细胞标志基因区域呈明显空间分离。红系细胞只与红系祖细胞保持稳定正相关，而与C1Q阳性或其他巨噬细胞无稳定关联。12周胚胎肝、健康成人骨髓凝块、Xenium分析和组织学成像均支持这一结果。
配体—受体分析将ICAM4相关黏附作用确定为人类红系细胞簇的重要候选机制。阻断ICAM4显著减少CD71和CD235a双阳性红系细胞数量，降低大于5个细胞的红系细胞簇数量和体积，并抑制纯红系培养中的集落形成；阻断RHAG或αv整合素未产生同等效应。MDS患者红系细胞的ICAM4表达下降，红系空间聚集遭到破坏，而治疗后部分患者的红系空间相关性重新增强，说明红细胞岛结构具有疾病相关性和一定可逆性。
直接推翻了红细胞岛结构在不同物种中高度保守的默认假设。 研究证明，小鼠以C1q阳性中央巨噬细胞为核心，而人类主要依赖红系细胞自身聚集。
实现多尺度空间验证。 Visium提供组织范围的无偏扫描，Xenium提供单细胞验证，组织学和功能实验进一步形成完整证据链。
从空间描述推进到功能机制。 研究不仅观察到人类红系细胞聚集，还通过抗体阻断和集落实验确定ICAM4是维持其结构和增殖的重要分子。
连接临床疾病和治疗反应。 MDS患者的红系细胞簇受损、ICAM4下降且治疗后部分恢复，使空间结构具有潜在生物标志物价值。
人类胚胎肝和健康骨髓样本数量相对有限，且多为特定孕周或中年个体，尚不足以覆盖完整年龄、性别和疾病谱。
Visium点位包含多个细胞，部分细胞共定位关系依赖单细胞数据去卷积；尽管Xenium进行了验证，仍可能遗漏未纳入探针的低表达分子。
ICAM4功能主要在类器官和体外红系培养中验证，尚缺少人体内直接干预证据，也未完全阐明ICAM4在红系细胞之间的具体结合受体。
MDS治疗前后样本数量较少，且治疗方式和疾病阶段可能存在差异，暂不能证明红细胞岛恢复是治疗有效的原因而非伴随结果。
Longitudinal changes in DNA methylation in IDH-mutant glioma fuel disease progression through altered cell state differentiation
胶质瘤如何变凶？DNA去甲基化锁住干性状态
IDH突变胶质瘤、DNA甲基化、单核多组学、细胞状态、肿瘤演化
IDH1或IDH2突变胶质瘤通常起病较慢，但最终会复发并向更高级别、更难治疗的恶性肿瘤演进。IDH突变会导致D-2-羟基戊二酸积累，抑制DNA和组蛋白去甲基化酶，形成以广泛高甲基化为特征的胶质瘤CpG岛甲基化表型，即G-CIMP。然而，部分复发肿瘤会逐渐失去这一高甲基化状态并转变为G-CIMP-low。此前虽然已经知道G-CIMP-low与不良预后相关，但DNA甲基化丢失究竟是恶性进展的驱动因素，还是干性细胞增加后的被动结果，并不清楚。
这篇文章通过对19名患者的36份原发和复发IDH突变胶质瘤样本进行纵向单核多组学、全外显子组测序和系统发育分析，探索了DNA甲基化丢失如何改变恶性细胞分化、状态遗传性和肿瘤进化轨迹，这对于建立IDH突变胶质瘤的进展标志物、理解表观遗传可塑性以及优化IDH抑制剂等分化治疗策略具有重要意义。
研究发现，胶质瘤复发时不仅获得新的拷贝数异常和驱动基因改变，还出现恶性细胞整体性的DNA低甲基化。重要的是，同一肿瘤内部不同恶性细胞的G-CIMP状态相对一致，而不是只有少数干样细胞发生去甲基化。这说明甲基化丢失更可能位于细胞状态改变的上游，为整个恶性细胞群提供重新进入干性状态的表观遗传基础。
伴随G-CIMP水平下降，复发肿瘤中的神经祖细胞样和少突胶质前体细胞样干性群体增加，而星形胶质细胞样和少突胶质细胞样分化群体减少。差异甲基化和转录分析提示，PRC2靶基因区域、胶质瘤干细胞相关基因及干性转录因子结合位点发生低甲基化，使SOX10等干性程序更容易被重新激活。
通过由单细胞甲基化信息重建的高分辨率系统发育树，研究进一步发现，G-CIMP-low肿瘤中干性状态的遗传性增强，细胞更倾向于在代际间维持干样表型，而向分化状态转换的概率下降。这种变化相当于将原本具有分化出口的细胞状态系统，重塑为更容易自我维持和增殖的“干性吸引态”。
研究还发现，低甲基化并非仅见于IDH突变型星形细胞瘤，部分IDH突变并伴1p/19q共缺失的少突胶质细胞瘤同样可转变为低甲基化状态，并具有更高级别和更差生存结局。因此，G-CIMP低甲基化可能是跨IDH胶质瘤亚型的恶性进展标志，而不只是某一病理亚型的特征。
研究收集来自美国、加拿大和日本三家医疗机构的36份冷冻IDH突变胶质瘤样本，涉及19名患者，其中15名患者具有匹配的原发和复发样本，包括11组IDH突变星形细胞瘤和4组IDH突变少突胶质细胞瘤。所有肿瘤均携带IDH1 R132H突变，并结合临床分级、治疗及复发时间信息进行纵向比较。
每份样本均进行全外显子组测序，用于鉴定体细胞单核苷酸变异、插入缺失、拷贝数异常、CDKN2A纯合缺失和替莫唑胺相关超突变。研究同时开展10x Genomics 3′端单核RNA测序，在32份样本中保留158,143个高质量细胞核，用于解析恶性和非恶性细胞组成及细胞状态。
核心技术为联合捕获单核转录组和DNA甲基化组：同一个细胞核中的mRNA采用Smart-seq2检测，基因组DNA采用扩展代表性亚硫酸氢盐测序，即snXRBS检测，共获得2,117个转录—甲基化配对细胞核。研究根据684个已知在胶质瘤进展中易发生去甲基化的CpG位点构建单细胞G-CIMP评分，并以混合模型区分G-CIMP-high和G-CIMP-low细胞。
研究采用非负矩阵分解识别星形胶质样、少突胶质样、干样和增殖相关表达程序；结合DNA甲基化、转录表达和转录因子基序富集寻找潜在调控机制。最后，根据单细胞甲基化变异构建高分辨率谱系树，将转录细胞状态映射到树上，量化干性、分化、细胞周期和拷贝数特征的遗传性及状态转换概率。
全外显子组分析表明，多数复发肿瘤获得了新的遗传改变，但不同患者的进化路径差异较大。12组可评估配对样本中有4组获得CDKN2A纯合缺失；部分替莫唑胺治疗后肿瘤出现错配修复基因异常及C>T突变显著增加，也存在与治疗无关的自发超突变病例。
单核甲基化分析显示，复发肿瘤的平均G-CIMP评分显著下降，而非恶性细胞在不同时间点间保持稳定，说明低甲基化并非技术误差或组织组成变化造成。同一肿瘤内恶性细胞的G-CIMP评分方差通常较低，提示甲基化丢失是肿瘤细胞群体层面的状态，而非少数亚克隆特征。
细胞状态分析发现，复发肿瘤中干样细胞比例增加、分化样细胞比例下降，配对比较达到统计学显著性。甲基化变化在不同恶性细胞状态之间相对一致，因此更支持“低甲基化促进干性增加”，而不是“干性细胞增加导致总体甲基化下降”。
机制分析识别出三类联系：第一，PRC2靶区域低甲基化，使发育调控基因更容易被重新激活；第二，SOX10等胶质瘤干细胞相关基因表达上升；第三，低甲基化区域富集干性转录因子结合基序，可能提高这些转录因子的染色质可及性。
系统发育分析显示，G-CIMP-low肿瘤具有更高的干性状态遗传性、更少的分化方向转换和更高的细胞周期活性。TCGA队列验证发现，低甲基化IDH突变少突胶质细胞瘤富集于高级别肿瘤，并具有显著较差的总生存期，说明G-CIMP状态可能成为跨IDH亚型的预后指标。
真正的纵向多组学设计。 同一患者原发与复发肿瘤配对，减少了横断面研究中个体差异对进展机制判断的干扰。
同一细胞核同步获得转录和甲基化信息。 使研究能够直接分析表观遗传状态与细胞表型的对应关系，而不是在不同细胞群之间进行间接推断。
将谱系关系引入细胞状态研究。 通过单细胞甲基化系统发育树量化细胞状态遗传性和转换概率，使肿瘤可塑性从描述性概念变成可计算指标。
扩展G-CIMP-low的临床范围。 证明低甲基化也可发生于IDH突变少突胶质细胞瘤，并与较差预后相关。
snXRBS虽然提高了非启动子区域覆盖率，但单细胞甲基化数据仍然稀疏，需要通过基因组窗口平均降低缺失值影响。
纵向临床样本珍贵但总体数量仍有限，每个样本获得的配对多组学细胞核数量也相对较少。
snXRBS适合推断大尺度拷贝数变化，但不足以在单个细胞中可靠识别点突变，因而难以精确连接特定基因亚克隆和细胞状态。
研究主要为回顾性观察，不能完全区分治疗、复发时间和自然进化对甲基化变化的贡献，也尚未通过实验直接证明去甲基化足以驱动干性增强。
A blended genome and exome sequencing method captures genetic variation in an unbiased and cost-effective manner
一次测序兼顾全基因组与外显子组：BGE破解成本难题
混合基因组外显子组测序、低深度全基因组、遗传变异、基因型插补、祖源多样性
大规模复杂疾病遗传学研究通常面临一个核心矛盾：高深度全基因组测序能够覆盖编码和非编码变异、结构变异以及不同频率的遗传变异，但成本较高；基因分型芯片价格较低，却受预设探针限制，对非欧洲祖源人群的遗传变异存在明显选择偏倚；高深度外显子组测序能够较好地识别罕见编码变异，却缺少全基因组常见变异信息。传统上，研究者往往需要分别进行芯片和外显子组测序，增加实验流程、批次效应和数据整合成本。
这篇文章通过开发混合基因组—外显子组测序方法BGE，在同一次测序中同时获得1至4倍低深度全基因组数据和30至40倍深度全外显子组数据，探索了如何在成本、覆盖范围和变异检测准确性之间取得平衡，这对于扩大大规模基因组研究、降低非欧洲祖源人群的技术偏倚以及促进全球基因组研究公平性具有重要意义。
BGE的核心思想是将PCR-free全基因组文库与经过外显子捕获的文库按一定比例混合后共同测序。经过六轮比例和数据量优化，研究者最终确定33%外显子捕获文库加67%全基因组文库的方案。该方案可以利用深度外显子数据直接检测罕见编码SNV、indel和编码区拷贝数变异，同时利用低深度全基因组数据和参考单倍型面板插补数千万个常见及低频基因组变异。
研究将BGE应用于PUMAS项目的53,000余份样本，参与者包括非洲、非裔美国人、拉丁美洲和混合祖源人群。结果显示，不同队列中外显子组平均深度超过30倍、编码区检出率超过0.99，绝大多数样本达到预定质量标准。低深度全基因组数据经GLIMPSE2和HGDP+1000 Genomes参考面板插补后，与Illumina全球筛查芯片基因型高度一致。对于次要等位基因频率不低于1%的变异，准确度通常达到95%以上；更低频变异也普遍达到90%左右，而且在非洲、欧洲和美洲局部祖源片段上保持较稳定表现。
在编码区拷贝数变异方面，跨越至少3个外显子的缺失和重复具有约90%的阳性预测值；跨越5个以上外显子的CNV召回率约87%，高覆盖BGE样本中可达到100%。其基本实验成本约为高深度全基因组测序的28%，即每个样本约99美元，而30倍全基因组测序约350美元。文章因此将BGE定位为介于芯片加外显子组和深度全基因组之间的可扩展方案，尤其适合需要数万至数十万样本的疾病关联研究。
BGE首先从基因组DNA制备PCR-free全基因组文库，并进行定量、归一化和混池。研究者从PCR-free文库中取一部分进行PCR扩增和外显子捕获，随后将捕获后的外显子文库与原始全基因组文库按纳摩尔浓度混合。经过六轮迭代实验，比较不同WES:WGS比例、每个测序泳道样本数、总数据量、外显子覆盖和基因型一致性，最终确定33% WES加67% WGS、每个NovaSeq S4泳道相当于64个样本的方案。
研究建议每份样本产生超过9.5 Gb通过质量过滤的数据，超过90%的外显子目标区域达到10倍覆盖，10%至65%的读段位于外显子探针及其邻近区域。大规模验证纳入53,446份PUMAS样本，并基于测序深度、外显子质量、祖源分层、性别一致性和异常质量指标实施QC。
在变异检测方面，外显子SNV和indel采用高深度直接分型；编码区CNV使用GATK-gCNV检测，并在Simons Simplex Collection的400名家系成员中以高深度WGS结果为金标准验证。全基因组结构变异分别采用VISTA和INSURVeyor检测。
低深度全基因组数据通过GLIMPSE2进行定相和基因型插补，参考面板由HGDP与1000 Genomes高深度数据组成。研究将插补结果与相同个体的Illumina GSA芯片数据比较，并按照等位基因频率、编码或非编码区域以及局部非洲、欧洲和美洲祖源进行分层评价。
六轮优化显示，33%外显子加67%全基因组的混合比例可获得约37倍外显子和约2.5倍全基因组平均深度，超过96%的外显子碱基达到10倍覆盖。对于次要等位基因频率高于5%的常见变异，与30倍全基因组测序相比的基因型一致性R²超过0.99，且优于GSA芯片。
在53,446份多祖源样本中，不到1%的样本因外显子覆盖不足而失败，仅6份样本因全基因组平均深度低于1倍被剔除。各队列外显子平均检出率超过0.99、平均深度超过30倍、平均基因型质量超过39，最终约47,264份样本通过完整QC。
CNV分析中，BGE对跨越5个以上外显子的CNV召回率为87%，跨越3个以上外显子的阳性预测值约90%。在高覆盖BGE样本中，5个以上外显子CNV的召回率达到100%；高深度WGS确认的11个新发CNV全部被BGE识别，且未产生额外的新发假阳性。
对于全基因组结构变异，BGE检测结果具有较高特异性但召回率较低：插入的平均阳性预测值约79.31%，召回率约23.04%；缺失的平均阳性预测值约99.14%，召回率约34.78%。因此，BGE适合获得高可信结构变异，但不能替代深度WGS进行全面结构变异发现。
研究共插补超过6,700万个双等位基因SNP，经INFO分数过滤后保留3,000余万个高质量位点。所有队列在各频率区间的聚合R²均超过0.90。常见变异在不同局部祖源背景中的准确性均超过90%；罕见变异存在小幅祖源差异，但总体R²差异通常小于10%。
一次实验同时获得两类互补数据。 既保留高深度外显子对罕见功能变异的检测能力，又获得可用于全基因组插补的低深度读段。
在超过5万份多祖源样本中完成真实规模验证。 不只是小样本技术展示，证明了方案可以进入大型队列生产流程。
降低芯片的祖源选择偏倚。 低深度全基因组测序不依赖预先选择的探针位点，对非欧洲祖源和混合祖源人群更具可扩展性。
成本与信息量平衡较好。 约以深度WGS 28%的实验成本，获得接近“外显子组加芯片插补”的分析能力，并额外支持编码CNV检测。
插补准确性高度依赖参考单倍型面板。对于参考面板中代表不足的祖源和超罕见变异，性能仍会下降。
研究主要以GSA芯片作为大规模基因型一致性参照，但芯片本身覆盖有限且可能包含分型误差，并非真正的全基因组金标准。
BGE对全基因组结构变异的召回率较低，尤其难以发现罕见、新发和复杂结构变异，不能完全替代深度WGS。
文章列出的每样本成本主要是实验室测序和基础处理成本，不包括DNA提取、项目管理、长期存储、计算分析及数据共享等费用。
与成熟芯片流程相比，BGE仍需要专门的低深度定相插补、质量控制和多祖源参考面板，计算流程及跨中心标准化门槛较高。
致谢橙子牛奶糖（陈文燕），请用参考模版：We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX
感谢小可爱们多年来的陪伴， 我与你们一起成长~
此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。