027-62435310 | 
service@benagen.com | 
文献解读 | 人体抗体图谱革命：单细胞全长转录组测序+单倍型基因组组装揭秘血液未知抗体
英文标题：De novo antibody identification in human blood from full-length single B cell transcriptomics and matching haplotype-resolved germline assemblies
发表期刊：Genome Research
发表时间：2025.4.14
研究背景
免疫系统分为先天免疫和适应性免疫。适应性免疫中的 B 细胞通过 V(D)J 重排和体细胞超突变产生高度多样化的抗体。B 细胞成熟后分化为浆母细胞（plasmablasts）（短暂抗体分泌细胞）、浆细胞（plasma cells）（长期抗体分泌细胞）和记忆 B 细胞（MBCs）。当再次接触抗原时，记忆 B 细胞可迅速分化为浆母细胞，启动高效免疫应答。
长读长单细胞转录组测序能捕获完整转录本，明确可变区（V）与恒定区（C）的配对，克服短读长测序遗漏异构体和新转录本的问题。结合单倍型基因组组装，可精准注释等位基因特异性转录本，并过滤未在胚系中存在的无效克隆。然而，IGH/IGK/IGL 基因座的高度多态性和结构复杂性（如单倍型变异、重复序列）仍是解析抗体遗传机制的挑战。
本研究结合单倍型基因组测序与单细胞全长转录组测序，从同一供体的麻疹、腮腺炎、风疹（MMR）疫苗接种后血液样本中，首次实现了对活化 B 细胞抗体库的系统解析。通过纳米孔长读长测序技术，获得了高精度的重链与轻链转录本一致性序列，并合成功能性抗体验证其结合与中和能力。研究揭示了抗体糖基化模式、克隆多样性及单倍型特异性表达规律，为精准免疫医学提供了新工具。
图1 实验设计示意图
主要研究成果
1. 单倍型高精度组装与注释
研究人员先对供体全血提取DNA，进行纳米孔全基因组测序，组装鉴定到IGH区域位于79 Mb连续的contig中。通过手动解析杂合缺失和结构变异，首次实现了母源与父源单倍型的完全分型（图2A）。单倍型1缺失D3-3至D2-8基因簇，单倍型2缺失D6-6至D5-1。IGL J-C 基因簇包含9个J-C基因簇（3.5–5.6 kb），并发现新型IGLC3等位基因（与已知基因仅相差一个单核苷酸变异）。
图2 使用Flye和HapDup组装单倍型基因组
在注释的295个IG基因中，288个（97.6%）与 IMGT 或 OGRDB 数据库完全匹配。并鉴定到4个新的等位基因，包括2个IGLV、1个IGHV和1个IGLC。
2. 单细胞全长转录组揭示B细胞异质性
为解析疫苗接种后B细胞免疫应答机制，研究对接种MMR疫苗的受试者进行纵向采血，通过荧光激活细胞分选（FACS）捕获CD27++/CD38++抗体分泌细胞（ASCs）和CD27+/CD38−记忆B细胞（MBCs）。分选后分别获得约1000个ASCs和30,000个MBCs进行单细胞全长转录组测序。
图3 MBC的单细胞转录组数据分析
通过比较短读长与纳米孔测序平台的差异，UMI比较显示两种技术对ASCs和MBCs的检测结果高度一致。纳米孔测序凭借全长转录本覆盖能力，不仅支持基因表达分析，还可精确量化转录本表达差异。例如，MBCs中BCL2基因在cluster 3主要表达BCL2-201（占69%），而cluster 8中BCL2-202占比升至32%（图3D）。尽管两种异构体编码区相同，5'UTR存在差异，转录本的使用偏好表明了仅基于基因水平分析可能遗漏的重要生物学过程。
对ASCs的纳米孔测序数据聚类发现多个亚群：CD38++浆母细胞（增殖态，MKI67高/IRF4低）、浆细胞（非增殖态，MKI67低/IRF4高）及少量误分选的CD27+/CD38+激活MBCs（图4A,B）。浆母细胞与浆细胞均高表达免疫球蛋白重链（IGHC），而误分选MBCs则呈现CD20/HLA-DR高表达特征。
图4 ASC的单细胞转录组数据分析
研究特别解析了抗体同种型分析的准确性。短读长测序因恒定区序列相似性导致IGHG1和IGHG3亚型分配模糊，而纳米孔测序通过全长转录本覆盖实现精准判定，可明确区分IGHG1和IGHG3（图4C）。
3. 抗体糖基化表型
在抗体糖基化方面，研究检测了 ASC 和 MBC 中核心抗体糖基化基因的表达。在ASC中，IgA 分泌细胞的岩藻糖基转移酶（FUT8）表达显著高于B4GALT1，而IgG1分泌细胞中则呈现相反趋势（图 4E）。由于岩藻糖基化水平升高与ADCC功能呈负相关，结果提示分泌型IgA的ADCC活性可能被抑制，而IgG的ADCC活性则增强。
在 MBC 中，糖基化模式与特定同种型没有广泛关联。但MGAT3 和 SLC35C1 在部分MBCs中互斥表达，体现了它们在调节 ADCC 活性中的相反作用。
4. 抗体合成与功能验证
通过分析IGHC M外显子（编码跨膜区）的转录本比例，发现MBCs主要表达膜结合型抗体，而ASCs则主要表达分泌型抗体（图5A）。为重建高精度抗体序列，研究开发了定制流程，结合纳米孔测序与生物素探针富集技术，使MBCs的IG转录本富集130倍。最终从ASC和MBC中分别获得761和7653条功能性重轻链一致性序列（图5B）。序列准确性验证显示，CH1区与胚系序列一致性达96%（平均准确率99.984%）。
图5 单细胞IG序列和抗体合成的表征
CDR3长度分析显示，MBCs的CDR3中位数为14个氨基酸，与健康人群一致（图5C）。V区突变频率（平均7.7%）符合生发中心后B细胞特征，其中IgD型MBCs和浆细胞突变率超12%，且偏好使用JH6片段（图5D）。利用R10测序芯片优化后，功能性IGH序列比例从80%提升至95%。
基于 ASCs 中鉴定的抗体序列，研究人员合成 100 种单克隆抗体。经 ELISA 试剂盒和 FACS 筛选，发现 MMR-82 能特异性结合麻疹病毒（MeV）表面蛋白，能有效中和 MeV 感染（图 5F-G）。同时筛选出 MMR-8 可识别人类巨细胞病毒（HCMV），证明该方法在多病原抗体发现中的潜力。
5. 胚系特异性抗体谱分析
分析IGH/IGL基因家族的使用模式显示，MBCs偏好使用VH3/VH4基因家族及JH4片段（图6A）。IgD型MBCs显著富集JH6片段（63.89%，总MBCs的19.1%），与IgD类别转换B细胞特征一致。相比之下，浆母细胞和浆细胞的VH家族更为受限且偏倚。
图6 供体特异性V基因使用和突变率
基于个性化参考数据库分析，等位基因模糊性降低90%。单倍型来源追踪显示55%的IGH序列可分配至特定单倍型（图6B,C）。CDR区的替换/沉默突变比（R/S=3.7-3.9）显著高于FWR区以及随机突变的估算值（图6D）。该模式在ASCs及MBCs中均存在，印证抗体谱系的适应性演化特征。
总结
本研究通过单倍型基因组测序与单细胞全长转录组测序的联合，首次在个体化水平揭示了抗体多样性的遗传基础和功能特性。未来需更大规模的单细胞测序与胚系组装数据集，将有助于深入理解特定等位基因的单倍型特异性表达模式、等位基因沉默和失衡、体细胞超突变以及重轻链配对等机制，推动精准医学发展。
参考文献：
Beaulaurier J, Ly L, Duty JA, Tyer C, Stevens C, Hung CT, Sookdeo A, Drong AW, Kowdle S, Guzman-Solis A, Tortorella D, Turner DJ, Juul S, Hickey S, Lee B. De novo antibody identification in human blood from full-length single B cell transcriptomics and matching haplotype-resolved germline assemblies. Genome Res. 2025
