IF50.5｜非洲锥虫抗原表达层级由同源修复模板决定

单细胞测序揭示VSG基因双链断裂后，修复模板有无决定抗原切换路径。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

非洲锥虫通过周期性更换表面VSG蛋白逃避宿主免疫，其基因组中约2500个VSG基因大多沉默于亚端粒区，仅少数端粒近端表达位点（BES）可活跃转录。VSG表达呈半可预测的层级模式，但此前研究受限于群体水平分析，无法区分不同VSG的激活频率与表达后生长差异。已有研究表明，在活跃VSG基因中诱导双链断裂（DSB）可触发抗原切换，但断裂后具体修复机制及VSG选择规则仍不清楚。

科学问题

卡在哪里： 此前缺乏在单个细胞中同时追踪转录组变化和基因组重排的手段，无法在切换事件发生后立即确定新激活的VSG及其激活机制，因此难以判断VSG层级表达是由激活频率还是生长优势决定，以及同源序列在修复中的作用。

本文要回答： 本文旨在利用高灵敏度单细胞RNA测序，在单个锥虫细胞中解析VSG基因双链断裂后的修复路径，揭示决定抗原表达层级的基因组因素。

技术路线

作者主要用了：SL-Smart-seq3xpress、CRISPR-Cas9、BLISS、ATAC-seq、流式分选。

作者首先建立可诱导Cas9切割活跃VSG基因的细胞系，通过γH2A和BLISS验证DSB产生。随后开发了针对锥虫优化的SL-Smart-seq3xpress单细胞测序方法，利用SL序列引物替代模板转换，并加入index hopping过滤。在VSG-2或VSG-8中诱导DSB后，对单细胞进行转录组测序，通过分析ESAG转录判断切换类型，结合de novo组装检测嵌合VSG形成，并用ATAC-seq探究BES染色质状态。

图 1　Fig. 1。图1展示CRISPR-Cas9系统在VSG-2 CDS第1140位诱导DSB的验证。a为BES1示意图及切割位点；b的Western blot显示诱导后Cas9和γH2A表达增加，对照无变化；c的BLISS覆盖图显示切割位点附近信号富集；d的流式分析显示诱导后第4天超过94%的细胞失去VSG-2表达，证明该系统能高效触发VSG切换。

核心亮点

亮点 1｜SL-Smart-seq3xpress实现高灵敏单细胞VSG检测

作者开发的SL-Smart-seq3xpress在相同测序深度下检测到中位2876个基因和4640个UMI，显著高于10X Genomics的1052个基因和1552个UMI。在更高深度下可检测到16797个转录本，约占预测总mRNA的84%。混合表达不同VSG的细胞实验显示，该方法100%的细胞达到互斥表达阈值，而5′ Chromium有11%的细胞未达标。该平台的低背景和高灵敏度为后续单细胞水平解析VSG切换机制奠定了基础。

图 2　Fig. 2。图2展示SL-Smart-seq3xpress方法原理及性能。a为文库构建流程，省略模板转换，利用SL序列引物进行第二链合成；b比较不同方法检测基因数，SL-Smart-seq3xpress中位基因数（2876）远高于10X（1052）和Smart-seq2；c显示高深度下可检测16797个转录本；d为混合细胞实验，SL-Smart-seq3xpress中100%细胞达到互斥表达阈值，而5′ Chromium有11%未达标。

亮点 2｜VSG-2断裂后仅激活端粒邻近VSG

在VSG-2 CDS第1140位诱导DSB后，单细胞测序显示VSG-2表达迅速下降，总VSG转录本减少近10倍，提示BES转录停滞。至第2天，大多数细胞无主导VSG。随后出现一组明确的VSG激活，包括BES、亚端粒阵列、微染色体和metacyclic表达位点中的VSG，但所有新激活的VSG均紧邻端粒。两个生物学重复中激活的VSG集合高度相似，表明存在硬编码的激活偏好。第10天时，表达VSG-9和VSG-18的克隆开始占优，提示生长优势也参与塑造群体异质性。

图 3　Fig. 3。图3展示VSG-2切割后单细胞VSG表达动态。a为BES1及切割位点示意图；b为不同时间点单细胞VSG表达热图，显示VSG-2表达迅速下降，随后出现多种VSG激活，第10天VSG-9和VSG-18占优；c为激活VSG的基因组位置分布，全部位于端粒邻近；d为切换类型分析示意图；e显示转录切换比例随时间下降；f和g为ATAC-seq结果，显示BES1在切换后保持开放而BES17关闭。

亮点 3｜BES1转录赋予体外生长优势

分析切换类型发现，早期时间点（第3天）转录切换比例高达62.5%，但到第10天降至3%，而通过BIR重组到BES1的细胞占比升至95.4%。ATAC-seq显示，从BES1切换到BES17后，BES1染色质仍保持开放，VSG-2转录仅部分下降；反向切换时BES17迅速关闭。这表明BES1的转录在体外具有强选择优势，转录切换细胞可能被淘汰或作为中间状态等待BES1修复后重新激活。

图 4　Fig. 4。图4展示VSG-8切割后节段性基因转换产生嵌合VSG。a为重组BES1及切割位点；b为单细胞de novo组装的VSG转录本与VSG-8及同源基因比对，显示切割位点附近出现SNP变异；c和d为不同细胞中重组片段长度和供体基因的多样性，表明修复模板不固定，产生多种嵌合VSG。

亮点 4｜同源模板存在时DSB通过节段性基因转换修复

在表达VSG-8（与基因组中多个VSG有同源性）的细胞中诱导DSB后，超过99.4%的细胞未发生70-bp重复区重组或转录切换，而是在切割位点附近出现单核苷酸多态性（SNP）。单细胞de novo组装显示不同细胞中SNP位置和数量各异，且这些序列变异与基因组中同源VSG假基因匹配，表明发生了节段性基因转换，产生了新的嵌合VSG。在VSG-11中也观察到类似现象。这说明当存在合适同源模板时，锥虫优先通过节段性基因转换修复活跃VSG中的DSB。

图 5　Fig. 5。图5为模型图，总结DSB在活跃VSG中的两种修复路径：存在同源模板时通过节段性基因转换形成嵌合VSG；无同源模板时通过BIR将端粒邻近VSG复制到活跃BES。该图直观呈现了同源序列可用性如何决定抗原激活机制。

亮点 5｜修复模板可用性决定抗原激活路径

综合VSG-2（无同源模板）和VSG-8/VSG-11（有同源模板）的实验结果，作者提出模型：DSB发生后，若基因组中存在同源序列，则通过节段性基因转换修复，产生嵌合VSG；若无同源序列，则通过断裂诱导复制（BIR）将端粒邻近的VSG复制到活跃BES中。该模型解释了为何VSG-2切割后仅激活端粒邻近VSG，以及为何VSG-8切割后未发生BIR。这揭示了同源序列可用性是决定VSG激活层级的关键因素。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： SL-Smart-seq3xpress、CRISPR-Cas9、BLISS、ATAC-seq、流式分选
• 公开数据： 原文未给出公开组学登录号
• 代码： 公开仓库 https://github.com/colomemaria/scSwitchFilter
• 谁该点开原文： 研究抗原变异、DNA修复或单细胞基因组学的科研人员，尤其是关注锥虫或病原体免疫逃逸的学者。

带走一句

单细胞测序不仅能看转录组，还能推断基因组重排；设计实验时，切割位点的序列背景（同源性）会决定修复路径，影响表型解读。