VDJ-SEQ:重排免疫球蛋白重链基因的深度测序分析，以揭示的B细胞淋巴瘤克隆演变模式

摘要

This protocol describes an approach to interrogate the recombined immunoglobulin heavy chain VDJ regions of lymphomas by deep-sequencing and retrieve VDJ rearrangement and somatic hypermutation status to delineate clonal architecture of individual tumor. Comparing clonal architecture between paired diagnosis and relapse samples reveals lymphoma relapse clonal evolution modes.

摘要

了解肿瘤的克隆是了解参与肿瘤的发生和疾病进展的机制是至关重要的。此外，理解内肿瘤发生响应于特定微环境或治疗方法，所述组合物克隆变化可能会导致更复杂的和有效的方法的设计消除肿瘤细胞。然而，跟踪肿瘤克隆亚种群已具有挑战性，由于缺乏可区分的标记。为了解决这个问题，一个VDJ-SEQ协议创建通过利用VDJ重组和体细胞突变（SHM），B细胞淋巴瘤两个独特的功能，以跟踪弥漫性大B细胞淋巴瘤（DLBCL）复发的克隆演变图案。

在这个协议中，新一代测序（NGS）库，索引潜在的扩增重排的免疫球蛋白重链（IgH重）VDJ区，构建由对初步诊断和复发DLBCL SAMP的LES。每个样品平均超过五十万的VDJ序列测序，同时含有的VDJ重排和SHM信息后获得。此外，自定义生物信息学管道开发，以充分利用序列信息FR3区​​域，VDJ剧目这些样品中的表征。此外，该管道允许重建和个体肿瘤的克隆构，这使得肿瘤的诊断和演绎的诊断和复发性肿瘤对之间克隆演变模式内的克隆异质性的检查的比较。当应用这一分析的几个诊断复发对，我们发现关键证据多个独特的肿瘤进化模式可能会导致DLBCL复发。此外，这种方法可以扩展为其它临床方面，如识别微小残留病，监测复发进展和治疗反应和免疫repertoir的调查ES在非淋巴瘤上下文。

引言

癌症是一种克隆性疾病。自从三十年前，当彼得C.诺维尔提出了癌症克隆演化模型^1，许多研究都试图肿瘤样本中剖析克隆群体和重建克隆扩增和发展模式，背后的肿瘤发生过程中^2。最近，全基因组测序，使调查人员深吸看克隆异质性和演化^3,4。然而，由于缺乏在许多细胞类型易于处理的标记物，它是很难推断精确的克隆构和进化路径。幸运的是在成熟B细胞，许多淋巴组织的恶性肿瘤，包括DLBCL，源于自然的克隆性标志。响应于抗原刺激，每个B细胞可以从这些区段的大池通过加入为V _H（可变的），一个D（多样性）形成单个生产性的IgH VDJ序列和_歼轰（接合）片段一起。在这个公关ocess，原始序列的小部分可以被删除和附加非模板的核苷酸可以添加到创建一个唯一的VDJ重排。这个特定的VDJ重排可在该B细胞的所有子代继承，因此标注个体成熟的B细胞和其后代^5。此外，SHM发生在随后的生发中心（GC）的反应重组的VDJ序列引入另外的突变的抗体池的膨胀和抗体亲和⁶的提高。因此，通过比较和对比已经经过这些过程淋巴瘤样品的VDJ和SHM图案，肿瘤内的异质性，可以划定及疾病的克隆演变路径可以被推断出来。

先前，VDJ重排和SHM可以通过PCR鉴定扩增重组的区域，克隆PCR产物，并随后Sanger测序来获得序列信息。这种方法我的低通量和低产量，检索整个重组的VDJ剧目只有一个非常小的部分，并阻碍给定样品中的克隆群的整体表现的表征。一种改进的方法是由来自VDJ PCR产物生成索引NGS测序文库并进行PE 2x150碱基序列，获得半数以上一万元左右的样品重组VDJ序列创建。此外，定制的管道的开发是为了进行质量控制（QC），对齐，滤波器的VDJ测序读数，以确定每次读取的重排和SHMS，以及执行系统发育分析每个样品的克隆构。此外，一种新的方法已经建立，以进一步表征克隆演变模式在不同疾病阶段收集的样本。

我们应用这种技术来DLBCL患者样本。 DLBCL是频繁复发的非霍奇金淋巴瘤的最多一个第一种侵袭性IRD的患者^7。 DLBCL复发通常发生早，在2至3年的初步诊断，但也有一些确实发生5年⁸之后。愈后复发的患者差，只有10％实现3年无进展生存期，由于有限的治疗选择。这是基础，迫切需要新的方法来治疗DLBCL复发^9,10。然而，DLBCL复发相关的分子机制仍知之甚少。具体地，克隆的异质性，在诊断和克隆演变的DLBCL复发发育过程中的作用是当前未表征，因此很难确定准确的和有用的生物标志物来预测复发。为了解决这些问题，我们运用我们的VDJ测序的方法对多对匹配的初步诊断复发DLBCL样本对。复发的两个不同的克隆的进化方案从克隆构的诊断和复发SAMP之间的比较出现莱这表明多个分子机制可能参与了DLBCL复发。

研究方案

1. VDJ放大

1.1）的肿瘤样本DNA提取

1. 提取冻结十月嵌入式正常或恶性组织的薄片（10-20微米）的DNA。
   1. 消化十月三十零日薄于4毫升核酸裂解缓冲液（0.0075摩尔Tris盐酸，pH值8.2; 0.3M NaCl的; 0.002，1M ​​的Na ₂ EDTA）中削减了一个低温恒温器切片包埋组织切片用蛋白酶K（0.5毫克/毫升，终浓度）和0.625％SDS中在37℃水浴中过夜15ml离心管中。
   2. 加入1毫升饱和NaCl（5米）的消化混合物，剧烈摇晃15秒。
   3. 离心在1100×g离心15分钟，在室温下。
   4. 上清转移到一个新的15ml离心管中，加两卷100％的乙醇。
   5. 通过倒转试管6-8次混匀。离心机在5000×g离心60分钟，在4℃下以收集沉淀的DNA。
   6. 用70％乙醇洗两次DNA沉淀。 CENTRifuge在5000 XG为每次收集沉淀15分钟。
   7. 溶于100-400微升TE缓冲液中的DNA在室温下在振荡器上过夜。的DNA产量是5至200微克（50-500纳克/微升之间最终浓度）取决于组织的大小
2. 提取福尔马林固定，石蜡包埋（FFPE）正常或恶性组织薄切片（10-20微米）的DNA。
   1. 孵化石蜡切片在1毫升二甲苯在1.5ml微量离心管中室温下每次10分钟两次，以解除paraffinize。通过在13000×g离心纺丝5分钟在室温下收集的组织切片。
   2. 孵育在1ml 100％乙醇部分两次在室温下，每次以除去残余的二甲苯10分钟。通过在13000×g离心纺丝5分钟在室温下收集的组织切片。石蜡是完全在这一点上溶解并只将组织部分的剩余量。
   3. 空气干燥，在室温TE的章节mperature 10-15分钟。
   4. 使0.5毫克/毫升蛋白酶K溶液用1×PCR缓冲液（10×从PCR缓冲液稀释无核酸酶的水）。添加蛋白酶K溶液的样品中的50-100微升最终体积孵育过夜，在37℃。
   5. 加热在95℃下对样品进行10分钟以灭活蛋白酶K.
      注意:在这个步骤中，样本DNA溶解在PCR缓冲液，并且可以在步骤直接1.2和1.3被使用。的DNA产量是0.5至20微克（10-200纳克/微升之间最终浓度）取决于组织的大小。

1.2）DNA的质量评价

1. 混合0.25微升Taq DNA聚合酶用45微升主混合物的从在一个PCR管商业梯试剂盒。
2. 添加从1.1.1.7或1.1.2.5制备5微升DNA导入PCR管并通过上下抽吸拌匀至少5倍。
3. 使用以下条件来扩增该DNA:95℃7分钟;接着45秒的35个循环的95℃，45秒60℃，和90秒72℃;然后72℃10分钟，并保持在15℃。
4. 准备在TBE液（Tris /硼酸盐/ EDTA）在2％琼脂糖凝胶。
5. 混合20微升与4μl的6×上样染料PCR反应，并装载到2％琼脂糖凝胶。
6. 染色用0.5微克/毫升溴化乙锭溶液中的琼脂糖凝胶和检测PCR产物用凝胶想象系统。注意:得到5- PCR产物，在100，200，300，400，和600碱基对大小的样品将被继续产生的VDJ扩增子。
   小心:溴化乙锭是一种强效的诱变剂。戴防护装备， 如手套非常谨慎，请处理和处置成每个机构的指导方针特定的容器。

1.3）VDJ PCR

1.3.1）的框架区1中扩增重组的IgH VDJ段（IgVHFR1）

1. 从管标签混合45μl反应混合液编为"混合2"商业体细胞超变异IGH检测方法的凝胶检测试剂盒的，0.25微升Taq DNA聚合酶，和从1.1.1.7或1.1.2.5在PCR管制备5μl的样品DNA。
2. 使用以下条件来扩增该DNA:95℃7分钟;接着45秒的35个循环的95℃，45秒60℃，和90秒72℃;然后72℃10分钟，并保持在15℃。
3. 解决整个PCR产物通过电泳在2％琼脂糖凝胶。
4. 染色用0.5微克/毫升溴化乙锭溶液中的琼脂糖凝胶和检测PCR产物用凝胶成像系统。单克隆扩增，预计用310-380基点的大小范围。
5. 切下含有310-380碱基之间的单克隆扩增的凝胶部分。
6. 使用根据制造商的协议标准的凝胶提取试剂盒纯化从切下的凝胶中的DNA。注意:对于FR1片段可以得到的样品，没有必要充分IgV结构HFR2和IgVHFR3。

1.3.2）的框架区2中扩增重组的IgH VDJ段（IgVHFR2）

1. 从标记为从1.1.1.7或1.1.2.5在PCR管准备凝胶检测试剂盒"B管"商业IgH基因克隆测定的，0.25微升Taq DNA聚合酶，和5μl样本DNA的管混合45微升主混合物。
2. 使用以下条件来扩增该DNA:95℃7分钟;接着45秒的35个循环的95℃，45秒60℃，和90秒72℃;然后72℃10分钟，并保留在15℃。
3. 解决整个PCR产物通过电泳在2％琼脂糖凝胶。
4. 0.5微克/毫升溴化乙锭染色进行可视化的PCR产物（多个）。单克隆扩增子可以由250-295基点的尺寸范围内观察到。
5. 切下含有250-295碱基之间的单克隆扩增的凝胶部分。
6. 使用标准的凝胶中提取纯化从凝胶切除DNA根据制造商的协议离子试剂盒。

1.4）优化VDJ PCR

1. 使用下面的修饰的PCR条件IgVHFR1和IgVHFR2使用次优的DNA扩增:95℃7分钟，40个循环的95℃60秒，60℃60秒，72℃90秒，在最后的延伸72℃下进行10分钟。

2. VDJ扩增子文库制备和测序

2.1）文库制备

2.1.1）末端修复

1. 从1.3传输的VDJ PCR产物成PCR管并从DNA样品制备试剂盒中添加悬浮缓冲液，使体积达60微升。
2. 加入40微升末端修复混合的DNA样品制备试剂盒调匀。
3. 孵育在30℃下进行30分钟的反应中预热的热循环仪（30℃）与预先加热盖在100℃下。
4. 混合136微升磁珠和24微升的PCR克RADE水首先在1.5ml管中，然后传送整个最终修复反应从2.1.1.3至1.5毫升管，并与小珠溶液充分混合。
5. 在磁性支架2分钟地方管，以允许从溶液中珠的分离。吸出上清液，并用80％的新鲜制备的乙醇的珠两次，而该管上的磁性支架上。
6. 完全吸出乙醇溶液，并允许珠粒空气干燥15分钟，在室温下搅拌15分钟。
7. 以管从磁力架，重悬珠在17.5微升再悬浮缓冲液。
8. 地方管重新安装到磁静置2分钟，以从所述重悬缓冲液单独的珠。除去重悬浮缓冲液（现在最终修复产物再悬浮于重悬浮缓冲液）到一个干净的PCR管。

2.1.2）A尾

1. 添加12.5微升拖尾混合到含有末端修复产品中的PCR管调匀。
2. 孵育将该PCR管在37℃下30分钟在预先加热的热循环仪（37℃）与预先加热盖在100℃下。

2.1.3）接头连接

1. 加入2.5微升再悬浮缓冲液，2.5微升式反应，以及2.5微升DNA适配器指数进入尾反应。
2. 孵育在30℃进行30 MINL反应在预加热的热循环仪（30℃）与预先加热盖在100℃下。
3. 添加5微升站连接缓冲液到每个反应调匀。
4. 拌42.5微升充分混合磁珠清洁以下步骤2.1.1.5至2.1.1.8反应。加入50微升重悬浮缓冲液以洗脱衔接子 - 连接产物。
5. 通过使用50微升充分混合AMPure XP珠粒清洗衔接子连接的产物进行第二次，和洗脱在25μl的重悬缓冲液的产品到一个干净的PCR管。

2.1.4）扩增DNA片段

1. 添加5微升的PCR引物鸡尾酒和25微升的PCR主混合物含有衔接子连接的产物，将该PCR管并充分混合。
2. 在一个预编程热循环具有以下条件进行扩增:98℃30秒; 10个循环的98℃，10秒，60℃30秒和72℃30秒;然后72℃5分钟，并保持在10℃。
3. 清理反应用50μl的充分混合的磁珠，并且洗脱最终产品在30微升缓冲液重悬。

2.1.5）库验证

1. 量化，使用制造商的协议荧光计的最终产品。最终库浓度为2.5〜20纳克/微升之间。
2. 通过使用具有根据制造商的协议高灵敏度的DNA芯片的分析仪器评价最终产品的质量。期望与任一IGVHFR1，IGVHFR2或IGVHFR3预期大小的单条带。电子库产物的xpected大小等于原始的VDJ扩增子的大小加上适配器（〜125 bp）的大小。

2.2）的VDJ-PCR库池和测序

1. 使用下列公式计算每个文库级分的摩尔浓度:纳米= [（毫微克/ 1,000）/ bp的* ^660]×10 9，其中n是在步骤2.1.5.1测定的VDJ-PCR文库的浓度和基点是峰在VDJ-PCR库的大小测定步骤2.1.5.2。
2. 稀释库为2nM（10微升总）与无DNA酶的水。
3. 结合将10μl稀释为2nM库一起放入1.5ml试管。
4. 添加的PhiX尖峰控制等体积为1.5 ml管。
5. 加载游泳池晚上7点集中到流动池。
6. 测序使用根据制造商的协议的配对末端150循环音序器的库。

3.数据分析

注意:SUMM在本节中使用的生物信息学脚本进制可以发现作为辅助代码文件。

3.1）校准和QC

1. 地图配对末端测序读针对人IGH V，D和J区数据库使用基于"BLASTN'可得自NCBI与缺口开放罚分2，1缺口延伸罚一个适于核苷酸BLAST算法从IMGT网站¹¹下载的， 7字长，和¹⁰ -4 e值的阈值。丢弃所读一对不映射到既是IGH V及J区。
2. 算剩余的读对具有IGH V和J的映射，以获得匹配的VJ组合在所有读取从样品上的测序运行得到的对的频率。计数读一对，如果它被映射到一个IGH V和J基因，然后将其等位基因增加相应的计用于特定VJ组合。
3. 排名从3.1.2高得到的每个重组VDJ区的数美国东部时间到最低。该VDJ地区具有最高的读取次数被定义为主要的重排组合。
4. 丢弃，覆盖在3.1.3中标识的域大的重排的小于35％比对的序列。

3.2）SHM档案鉴定

1. 计算所有SHM模式在整个读取步骤3.1.3。定义每个SHM图案作为亚克隆。
2. 计数分别对应于不同的唯一SHM模式，将亚克隆的数量，和对读取映射到各个亚克隆该号码。

结果3.3）图形表示

1. 执行使用基于邻域加入与R包"猿"的方法，根据制造商的协议及其对应的SHM模式的亚克隆系统发育分析。计算从各不同路线的距离矩阵来重新创建亚克隆系统发育植根于V的种系序列Ĵ组合为每个配对样本集的问题。
2. 使用每个亚克隆作为字符载体的核苷酸序列，并计算在主要的VJ重排为诊断和使用Levenshtein距离量度相应复发样品，其中数学2路线之间的距离为d _的x，y给出的所有亚克隆之间的近似串距离（I，J）其中d _{的x，y（I，J）=} 1，如果I≠j和否则为0，然后求和这个函数在字符串的长度相对应的核苷酸序列i和j。
3. 图形方式显示所产生的亚克隆的距离及其相应的亚克隆计数以下列两种方式中的矩阵:
   1. 使用ř包"。MASS"根据制造商的协议来多维尺度应用于亚克隆距离矩阵和生成的二维坐标原理图，然后将其与对数作图沿的亚克隆计为圆的半径。
   2. 构造基于所述Levenshtein距离，其中该顶点对应于从主要VJ重排所产生的每一个不同的亚克隆无向图。
      注意:如果两个不同的克隆之间的距离等于一个则存在这两个顶点之间的边。如果距离大于一，那么就没有边缘连接它们。
   3. 使用tkplot函数在"IGRAPH"R包具有根据制造商的协议的镰河合布局的绘制3.3.3.2的曲线图。顶点的半径成正比的克隆映射到它的总数量的立方根和遮光使用红色和蓝色的范围来表示的克隆的，要么属于诊断（蓝色）或复发（红色）的样品的比例。

3.4）异质性（熵）测量

1. 检查的数量和亚克隆计数体细胞突变的主要VJ重排的每个个体出现的图案配对设置的样品。
2. 计算经验熵和残余经验性熵调整不同的克隆在每个样品中的使用标准直方图基于熵估计的数量。

结果

VDJ测序（VDJ-SEQ），包括DNA提取的整体过程，重组的VDJ区域扩增和纯化，测序文库构建，读取处理，和系统发育分析，表示在图1中。按常规5-200微克的DNA可以从被检索冷冻固体组织切片或从福尔马林固定的石蜡包埋组织切片0.5-20微克的DNA。视质量，重排的模式，和个体样本SHM程度，可以得到各种的VDJ PCR产物来自不同样品（图2），包括IGVHFR1（310-380碱基），FR2（250-295碱基?...

讨论

因为几乎无限数量的通过的VDJ重排和SHM在人B细胞的IGH轨迹编码序列信息迭代的，检查整个IGH剧目由高通量深测序证明是划定克隆一种有效和全面的方式和子克隆B细胞群。此外，该策略可用于通过比较沿该疾病的不同阶段收集的患者样本的克隆和子克隆构研究B细胞肿瘤的发展，缓解，复发的克隆演变路径。虽然很容易在实验室环境中使用可商购的引物进行从初级样品DNA重排的VDJ序列的扩增，扩增?...

材料

Name	Company	Catalog Number	Comments
Ethanol	VWR	89125-170	200 proof, for molecular biology
TE buffer	Life Technologies	12090-015	10 mM Tris·Cl, pH 8.0; 10mM EDTA
Xylenes	VWR	EM-XX0055-6	500 ml
Proteinase K	Life Technonogies	25530-015	100 mg
Deoxynucleotide triphosphate (dNTP) Solution Mix	Promega	U1515	10 mM each nucleotide
10x PCR Buffer	Roche	11699105001	Without MgCl2
AmpliTaq Gold DNA Polymerase with Gold Buffer and MgCl2	Life Technonogies	4311806	50 µl at 5 U/µl
Specimen Control Size Ladder	Invivoscribe Technologies	2-096-0020	33 reactions
IGH Somatic Hypermutation Assay v2.0 - Gel Detection	Invivoscribe Technologies	5-101-0030	33 reactions, Mix 2 for IGVHFR1 detection
IGH Gene Clonality Assay - Gel Detection	Invivoscribe Technologies	1-101-0020	33 reactions, Tube B for IGVHFR2 detection
GoTaq Flexi DNA Polymerase	Promega	M8291	20 µl at 5 U/µl
10X Tris-Borate-EDTA (TBE) buffer	Corning (cellgro)	46-011-CM	6x1 L
50X TAE BUFFER	VWR	101414-298	1 L
Ethidium bromide solution	Sigma-Aldrich	E1510	10 mg/ml
25 bp DNA Ladder	Life Technologies	10597011	1 µg/µl
100 bp DNA Ladder	Life Technologies	15628-019	1 µg/µl
UltraPure Agarose	Life Technologies	16500-500	500 g
40% Acrylamide/Bis Solution	Bio-Rad Laboratories	1610144	500 ml
QIAquick Gel Extraction Kit	Qiagen	28704	50 columns
Agencourt AMPure XP	Beckman Coulter	A63881
Qubit dsDNA High Sensitivity Assay Kit	Life Technologies	Q32854
High Sensitivity DNA Kit	Agilent Technologies	5067-4626
2100 Bioanalyzer	Agilent Technologies
PhiX Control v3	Illumina	FC-110-3001
MiSeq	Illumina
Qubit 2.0 Fluorometer	Life Technologies	Q32872
Resuspension buffer (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
End repair mix (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
A-tailing mix (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
Ligation mix (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
DNA adaptor index (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
Stop ligation buffer (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
PCR primer cocktail (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
PCR master mix (Illumina TruSeq DNA Sample Preparation Kit v2)	Illumina	FC-121-2001
Magnetic stand	Life Technologies	4457858
Gel imaging system	Bio-Rad Laboratories	170-8370