逆转录病毒整合位点的扩增，新一代测序和基因组DNA制图

摘要

我们描述了一种方案，用于从感染细胞的基因组 DNA 中扩增逆转录病毒整合位点，对扩增的病毒-宿主连接进行测序，然后将这些序列映射到参考基因组。我们还描述了使用 BEDTools 量化整合位点相对于各种基因组注释分布的技术。

摘要

在本地和全球尺度上表现出逆转录病毒整合签名的偏好。这里，我们提出了（1）的生成用连接介导的PCR（LM-PCR）扩增和新一代测序（NGS）逆转录病毒整合位点的不同的库的详细的协议;（2）绘制的每个病毒 - 的基因组位置使用BEDTools交界处，和（3）分析统计相关数据主机。从感染的细胞中提取的基因组DNA是通过消化用限制性酶或通过超声处理分段。合适的DNA末端修复后，双链连接子连接到DNA末端，和半巢式PCR是使用这两个病毒的长末端重复（LTR）端部和连接接头DNA互补的引物进行的。 PCR引物携带NGS期间用于DNA聚类所需的序列，否定为单独的适配器连接的要求。质量控制（QC）下进行，以评估DNA片段大小分布和适应呃之前NGS DNA结合。序列输出文件过滤含有LTR-读取，并确定LTR和连接器的序列裁剪掉。修剪宿主细胞序列被映射到使用BLAT参照基因组和被过滤为最小97％的同一性的参照基因组中的独特点。独特的整合位点是审查相邻核苷酸（nt）的顺序和分配相对于不同的基因特性。使用该协议，高复杂度的整合位点文库可从基因组DNA在三天构成。该包括易感组织培养细胞至整合位点分析的外源病毒感染的整个协议可以因此在大约一到两周进行。这项技术的应用近期涉及到艾滋病毒感染患者的整合位点的纵向分析。

引言

病毒DNA（VDNA）到宿主细胞基因组中的整合是在逆转录病毒的生命周期的一个重要步骤。积分是由病毒酶整合酶（IN），其执行，导致建立稳定地插入原病毒¹的两个不同的催化过程来完成的。在亚基接合是通过逆转录生成的线性VDNA的端部，形成具有VDNA高阶intasome端通过一个IN多聚体^2-4保持在一起。在切割中称作3'-处理的处理顺序的3'从不变5'-CA-3下游VDNA的端部"，留下凹入3'在每个VDNA末端^5-8活性羟基结束。该intasome随后导入细胞核作为大型组件的主机的一部分，并称为preintegration复合物^{（PIC），9-11}的病毒蛋白。遇到细胞靶DNA（T-DNA）后，使用VDNA 3'-羟基GRO起坐裂解TDNA顶部和底部股线以交错的方式，同时通过链转移^12,13的过程中加入VDNA到TDNA 5'磷酸基团。

在本地和全球尺度上表现出逆转录病毒整合位点的偏好。在当地，共识的整合位点包括了从VDNA插入位点^14,15大约五到十个基点的上游和下游跨越弱保守的回文T-DNA序列。全球范围内，逆转录病毒针对具体的染色质注释^16。有七种不同的逆转录病毒属 - 通过ε，伦蒂和spuma阿尔法。慢病毒，包括HIV-1，有利于积极转录基因¹⁷机构内部整合，而gammaretroviruses优先整合到转录起始位点（起始位置跟）和有源增强地区^18-20。与之形成鲜明对比，泡沫病毒强烈向heterochrom偏见ATIC区域，如基因贫乏椎板相关领域^21。当地TDNA基地偏好由IN和T-DNA ^13,22,23之间核蛋白的接触特定网络决定的很大一部分。对于慢病毒和gammaretroviruses，相对于基因组注释集成是通过中和同源细胞因子^24-27之间的相互作用支配很大一部分。改变IN-T-DNA相互作用网络^13,22,23,28的细节，扰乱或重新工程宿主因子交互^25-27,29-32是行之有效的战略分别重新定位在本地和全球层面的整合。

用来编目逆转录病毒整合位点的DNA测序方法的功率在过去几十年的极大增加。整合位点用费力净化和手工克隆技术产生只是根据研究^33,34独特的网站屈指可数回收创举。LTR宿主DNA结的LM-PCR扩增与映射单个整合位点，以转化领域人类和小鼠草案基因组，与来自外源组织培养细胞的感染中回收网的增加至几百到上千¹⁷数量的能力的结合^18。 LM-PCR检测方法NGS的更近的组合已派出图书馆深度跃然而出。具体而言，焦磷酸测序，得到的独特的集成网站^30,35-38数以万计的量级，而库测序通过使用DNA聚类可以产生数以百万计的独特序列^19-21,39的。在这里，我们描述了扩增和测序使用DNA集群NGS逆转录病毒整合位点的优化LM-PCR方法。该方法结合了所需的适配器序列引入PCR引物，因此直接进入扩增的DNA分子，从而排除了要求用于向sequen之前一个额外的适配器连接步骤庆安^40。在生物信息学分析管道中，从LTR宿主DNA结到的独特的整合位点的映射到相关的基因组特征的原始测序数据的解析，也一般地描述。按照这一领域^36,38,41-43从以前的方法确定的协议的优先级，自定义脚本可以开发，以帮助在生物信息学管道的具体步骤的完成。该实用程序和协议的灵敏度通过扩增，测序和映射从在感染的1.0（MOI），以及一个滴定系列该DNA的近似多重感染的组织培养细胞中的HIV-1的整合位点与代表性数据示出通过未感染的细胞的DNA中稀释5倍的步骤为1的最大稀释度:15625，得到的6.4×10 ^-5的近似等效惯性矩。

方案

1.生成病毒原液

注意:此协议的湿式清洗台方面的流程图在图1中描绘的病毒原料生产和组织培养细胞的随后感染的细节通常将适用于不同类型的逆转录病毒。对于一些实验，靶细胞可能不表达内源性病毒受体（多个），并且在这种情况下，假型逆转录病毒颗粒携带异源病毒包膜糖蛋白的结构， 例如 ，从水泡性口炎病毒（VSV-G）对G糖蛋白，会需要感染^44,45。

注:与HIV-1工作时，须采取预防措施。虽然具体准则会有所不同，从机构的机构，所有的基于病毒的工作应在一个专门的，操作者限制生物安全柜（通常被称为组织培养罩）来进行。适当的个人防护设备包括面部防护，鞋套，手套双层和全身套装西装应随时佩戴。从病毒相关的实验产生的所有废液应与漂白剂（10％终浓度）灭活，以及包括固体废弃物都应该在处置之前进行高压灭菌。

1. 一个转染前一天，补充有10％（体积/体积）胎牛血清和1％的板3.3×10 ^6个 HEK293T细胞于10ml的Dulbecco改良的Eagle培养基（DMEM）的（体积/体积）青霉素/链霉素（10,000ù / ml的股票）在每一个100 mm培养皿。
   注意:辅以的DMEM上称为DMEM-FPS从这一点。
2. 在随后的日子，转染使用市售的转染试剂或磷酸钙细胞用质粒携带全长逆转录病毒分子克隆的10微克或9微克信封缺失单轮载体用VSV-G表达构建体的1微克。
   1. 孵育的C在5％的CO ₂的潮湿的细胞培养箱中，在37℃厄尔（在此条件下称"组织培养孵化"）。后约48小时，使用一个容积移液管收获含有病毒的细胞培养基，并通过它通过0.45μm的过滤器在重力作用下流动。
   2. 集中超速离心病毒为20万XG在4℃1小时。重悬在500微升DMEM-FPS含20单位DNA酶病毒颗粒，并在37℃孵育1小时。
      注意:DNA酶步骤有助于通过消除来自转染过程仍然存在的质粒DNA的冲击，以减少不必要的质粒序列的恢复。
3. 确定使用HIV-1 p24抗原捕获试剂盒，按照制造商的说明P24浓度^46。
   注:病毒浓度还可以通过逆转录酶活性测定^47,48确定。备选地，功能性病毒的水平可以通过测量MOI确定。这是通过使用荧光激活细胞与表达荧光报告基因，如绿色荧光蛋白病毒排序最容易完成的。与可能不支持感染的相同的水平进行优化的细胞系的初级细胞时MOI确定可以是特别有用的。

2.感染细胞与病毒

1. 板，每孔3.0×10 ^5个 HEK293T细胞于在2.5ml的DMEM-FPS一个6孔板，并在组织培养箱孵育过夜。
   注意:此协议中回收独特的整合位点的数目成正比在感染中使用的细胞和活性病毒的数量的数量。
2. 传染为500纳克/毫升的终病毒的p24浓度的细胞在500微升新鲜DMEM-FPS终体积用于在组织培养箱2小时，再加入2毫升的DMEM-FPS预热至每孔在37℃和继续孵化。
3. 在48小时后感染，删除媒体并用2ml磷酸盐缓冲盐水（PBS）洗涤细胞。加0.5 ml胰蛋白酶 - EDTA预加热到37℃，和几秒钟后目视检查细胞移位的孔中。
4. 加2ml的预温热的DMEM-FPS和悬浮轻轻向上/向下吹打使细胞与一个容积移液管〜10倍。将该溶液转移至含有18毫升预热DMEM-FPS一个75cm ²的组织培养烧瓶中，并孵育所述细胞在组织培养培养箱中。
5. 后，从感染的起始微创五天，收集细胞通过除去介质，用5ml PBS洗涤，加入2毫升预热胰蛋白酶-EDTA，并通过移液用5ml预热的DMEM-FPS悬浮。离心在室温下5分钟将溶液在2500 xg离心，并弃去上清液。
   注意:虽然在约48小时后感染^49,50在这些条件下的高原整合，需要培养的附加​​3天sufficientlý稀释，从基于细胞的DNA重组或病毒介导的自动积分导致未集成的DNA分子的浓度。
6. 提取使用市售的试剂盒（ 例如 ，见^51）将细胞沉淀基因组DNA。洗脱从用200μl的10mM Tris-HCl中，pH值8.5所提供的离子交换柱将DNA。
   注:细胞的等分试样应在48小时后的感染（步骤2.3）用于感染性测定，以确保适当的病毒感染前NGS分摊。

3.片段基因组DNA通过超声处理或限制性内切酶消化

注意:超声处理的碎片在一个几乎序列无关的方式的基因组DNA，因此碎片的优选模式时用低的预期回 ​​收率测序样品（ 例如 ，在相对 ​​低的MOI启动感染的患者的细胞或感染）。此外，超声允许一个区分的党重复PCR从独特的集成丘拉尔整合位点序列在同一地点，这是至关重要的，以区分含有原病毒细胞的感染患者（见下文第11步^）39,52-54的克隆扩增。
注意:将DNA应立即下游从上游的LTR被切割的LM-PCR过程中，以减少内部的病毒序列的扩增。限制性内切酶BglⅡ即位于43碱基下游从上游U5序列，这是用于与MseI-产生的DNA随后连接不相容结束作品以及与许多HIV-1株（ 图1B）。当通过超声准备DNA，内部裂解限制性内切酶应连接结扎后应用（ 见图1C - E和下面的步骤4.3）。

1. 用于超声处理，混合在不含核酸酶的水的基因组DNA的10微克至120微升的最终体积。超声处理使用的参数为500 bp的平均规模突破（两轮以下的段米:占空比:5％;强度:3;每一次脉冲周期:200;时间:80秒）。
2. 纯化使用PCR纯化试剂盒声振的DNA。修复DNA末端用DNA末端修复包，净化使用PCR纯化试剂盒的DNA。 A-尾用Klenow外的DNA ^-酶，净化使用PCR纯化试剂盒的A-DNA尾。请参阅^51,52的配套使用的其他详细信息。
3. 用于限制性内切酶消化，在100μl的缓冲液由生产和产生5'- TA突出酶的混合物（100单位每）供给的体积切割基因组DNA的10微克，37℃过夜，以及一个不相容酶如BglII位劈开从上游病毒LTR下游。使用PCR纯化试剂盒纯化该DNA的第二天。
   注意:限制性内切酶的无应的终端内切〜由该LM-PCR方案扩增病毒DNA端的30碱基对。该协议专门放大U5HIV-1 DNA的末端。

4.退火接头寡核苷酸和结扎支离破碎的基因组DNA

注:准备含有悬垂是与上述DNA片段（ 见表1中这个协议所使用的寡核苷酸序列）兼容的非对称接头。链接器超声处理的DNA用于必须包含一个兼容的T-3'突出，而对于MseI消化的DNA的连接体必须包含一个兼容5'-TA突出（ 图1）。短连接链必须还含有不可延伸的化学修饰，如3'-胺，来约束朝向感兴趣的DNA随后的扩增反应。
注意:当制备多个并联不同整合位点的库和/或当在相同的测序运行多路复用独特样本，它建议使用独特的接头对每个样品，以限制对样品交叉沾污的可能性PCR过程中通货膨胀。此附加意味着在半巢式PCR（以下描述）的使用独特接头的引物对每个样品。独特接头链和接头引物可以通过加扰表1中列出，同时保持相似的总体％GC含量和适用突出位置上的连接体寡核苷酸序列来设计。

1. 退火短和长接头链在35微升10毫摩尔Tris盐酸，pH值8.0-0.1毫摩尔EDTA（10μM的每种寡核苷酸的最终浓度）通过加热至90℃并缓慢冷却到室温以1°步每分钟℃。
2. 制备每基因组DNA样品的至少四个并联连接反应，其中含有1.5μM的连接接头，1微克片段化DNA，并在50μl800单位T4 DNA连接酶。在12℃结扎过夜。净化用PCR纯化试剂盒的第二天。
3. 用于通过超声处理制备的样品，消化有100U一个restric的纯化连接反应化的酶，根据制造商的下游，从上游LTR（ 例如 ，BglⅡ位为HIV-1）的切割推荐的条件过夜。纯化使用PCR纯化试剂盒将DNA。

5.扩增出病毒LTR-主机的基因组DNA结由半巢式PCR

注意:为了确保最佳文库多样性，至少4-8平行的PCR，这取决于所回收的连接反应的DNA浓度，应该对每个样品为两种PCR回合来制备。 DNA模板浓度应通过分光光度法进行量化。在这个协议中的PCR的第一轮和第二轮采用嵌套LTR特异性引物，但相同的特定连接体引物被用于两轮（ 表1）。第二轮LTR特异引物和用于DNA聚类以及测序引物结合位点的特异性连接子引物编码接头序列。嵌套LTR特异性引物还编码6核苷酸索引序列，WHICH可以为同一测序运行中多路复用库不同的引物之间变化。

1. 制备含有每管的成分如表2中列出的第一轮的PCR。
   注意:特定接头引物怀有互补的22个核苷酸的连接体，53℃，45％的GC含量，并且其3的熔融温度'末端位于15-16 bp的上游的3'末端不同的连接子长链（ 表1）。第一轮27 nt的LTR引物具有的59℃，48％的GC含量的熔融温度，和它的3'末端位于34碱基上游的HIV-1的U5末端。第二轮26 nt的LTR的引物的区域即HIV-1的LTR互补为60℃，50％的GC含量，并且其3'末端的熔融温度位于18碱基上游从病毒U5总站。建议在寡核苷酸的熔化温度和GC含量应如果用户模仿这些参数设计PCR引物用（包括与其他逆转录病毒使用^）21改变的序列。
2. 运行一次PCR轮在以下热循环参数:一循环:94℃2分钟; 30个循环:94℃15秒，55℃，30秒，68℃45秒;一个循环:68℃10分钟。
3. 池反应和使用PCR纯化试剂盒纯化。制备含有每管中的成分按照表3中的第二轮的PCR。使用在步骤5.2中所述的热循环仪参数运行所述第二轮PCR。凝聚反应和使用市售的PCR纯化试剂盒按照制造商的说明纯化DNA。
   注:各种与DNA集群NGS兼容推荐指数序列可用^71。

6.执行QC和NGS（通常由测序完成设施）

1. （QC测定＃1），使用氟确认步骤5.3文库DNA浓度米^55。简要地说，在200微升不含核酸酶的水的最终体积制备标准和实验样品。涡流管2-3秒，室温下培养2分钟，然后在荧光计读出的样品。
   注:样品应包含在15微升体积最小的为2nM文库DNA的最低浓度。
2. （QC测定＃2）确认用的基于磁带的测定法⁵⁶的DNA片段大小分布。
   注:理想的分布是比较宽泛的DNA峰中心在长约500基点。如果材料的显著量大于1 kb的大，则建议掺入大小选过程，以消除更长的DNA的物种，这将聚类过程中妨碍桥扩增。相反，如果一个显著峰是明显的约100至200碱基，引物二聚体可能PCR过程中形成的。在这种情况下，该过程应当被优化，以尽量减少引物二聚体的形成。
3. （QC检测＃3）有限公司通过定量PCR ⁵⁷并将通过适配器插入DNA文库的适当结合。
4. 制造商的应用程序文学下列情况进行NGS。利用10％的尖峰在（W / W）ΦX174的DNA，这将通过对测序运行提供均衡基础组合物优化实时质量度量。
   注:整合位点测序实验通常受到单端150基点（SE150）或配对末端150基点（PE150）测序。 PE150是特别有用的以捕获在每个DNA分子的连接体的连接点（ 例如 ，审议整合位时为宿主细胞克隆扩增的证据）。

7.使用自定义的Python或Perl脚本解析测序数据包含LTR-序列，作物离开LTR和接头序列，并映射到参考基因组与BLAT

1. 对于含LTR序列读取，作物LTR和接头序列从宿主基因组DNA序列路程，扫描FASTA文件这些序列导出到新FASTA文件。图裁剪读取到两个参照基因组（ 例如 ，人类基因组的版本hg19或GRCh38），并使用BLAT ⁵⁸的病毒基因组，与输出整合位点坐标输出到一个单独的.txt文件，使用以下设置:
   stepSize的= 6，minIdentity = 97，和maxIntron = 0
2. 解析BLAT输出.txt文件，删除autointegrations（ 即证据表明，LTR月底已整合到病毒DNA基因组的内部区域）等的序列映射到HIV-1的基因组，并创建.txt文件一个单独的输出中，所有重复的整合位点已经被浓缩成单一的，唯一的坐标命中。

8.创建一个包含15-NT间隔周围集成.bed文件，这些转换为FASTA文件，构建序列标志来展示基地周边首整合位点

1. 创建列出的基地间隔.bed文件每个整合位点。至少有15个碱基（5上游和下游10）的建议对序列的标志产生。通过使用BEDTools ⁵⁹ fastaFromBed功能，该命令从这些.bed文件的FASTA文件:
   fastaFromBed -fi /目录/至/参考/基因组/ -name -s -bed 15_base_pair_file.bed -fo output_file.fasta
   注意:不变病毒5'-CA-3'二核苷酸接合集成期间主办的DNA，和验证LTR末端到细胞DNA的交界处是一个重要的初始滤波器，以确定真正的整合位点。我们还从该主机DNA序列人口编译序列标识验证实验结果。作为逆转录病毒展示周围的整合位点^14,15签名碱的喜 ​​好，该序列标识用来验证该映射基因组位点产生通过中介导的整合，相对于其他的重组机制，如非同源DNA末端连接^60,61。
2. 使用WebLogo 3（http://weblogo.threeplusone.com/create.cgi）创建从FASTA文件序列标识。点击"选择文件"上传FASTA文件，并使用以下设置输出格式，PDF（矢量）;标志的大小，大的;第一个位置号，-5;标志齐全，-5〜5; Y轴的规模，0.1，Y轴间距抽动，0.5，颜色方案，经典（NA）。

9.创建中央碱基对.bed文件，检查样品交叉污染，以及地图独特的集成点的相对分配到相关的功能基因组

1. 由于逆转录病毒整合整个T-DNA链交错的方式出现时，调整整合位点的精确坐标，以反映目标站点复制相对于基因组功能基因组分布的正确映射的核心基点。
   1. 因此，5基点复制病毒如HIV-1，建立与中央基点从我偏移.bed文件通过向正链的整合绘图下游两个基地，两基地上游集成映射负链ntegration网站。
2. 来检查样品的交叉污染，通过使用BEDTools 相交函数相交中央碱基.bed文件两个不同的样品，并按照该命令计算不同的库之间的共同的整合位点的数目:
   bedtools相交-a central_basepair_1.bed -b central_basepair_2.bed -f 1.00 -r -s> overlap1v2.txt
3. 算，以量化的两个库通过使用以下命令中常见的部位的确切数目的输出overlap1v2.txt文件内的行数:
   厕所-l overlap1v2.txt
4. 下载的RefSeq注释.bed文件被用于从UCSC基因组注释数据库整合位映射的参照基因组的版本（ 例如 http://hgdownload.cse.ucsc.edu/goldenPath/hg38/da^{tabase）62。}
   1. 通过使用BEDTools功能相交相交，那是对样品产生的RefSeq的.bed文件下面这个命令中央碱基对.bed文件计算整合位点下跌的RefSeq基因中的数量:
      bedtools相交-a central_basepair_1.bed -b RefSeq_hg38.bed -u> RefSeq_sample1.bed
5. 算，以量化通过使用以下命令中的RefSeq基因落位点的确切数目的输出RefSeq_sample1.bed文件内的行数:
   厕所-l RefSeq_sample1.bed
6. 映射整合位点与感兴趣的任何其他注释的量的间隔.bed文件是可重复步骤9.3和9.4。下载来自UCSC基因组注释数据库感兴趣的参考基因组的最新CpG岛标注.bed文件作为指导步骤9.4。
   1. 计算整合位点的落下一定二内的数姿态（在本实施例中所示是一个5 kb的窗口）CpG岛通过使用BEDTools 窗口函数和跟随此指令的:
      bedtools窗口-w 2500 central_basepair_1.bed -b CpG_hg38.bed -u> CpG_sample1.bed
7. 算，以量化落入2.5 kb以内的上游或下游的CpG岛的通过使用下面的命令位点的确切数目的输出CpG_sample1.bed文件内的行数:
   厕所-l CpG_sample1.bed
8. 重复步骤9.6和9.7绘图整合位点附近的起始位置跟。产生RefSeq.bed文件，其中的基因组坐标映射到多个基因的候选版本进行了调整，以反映只有一个单一的基因存在于该位置。这可以防止周围的整合位点基因密度高估。计算通过使用BEDTools 窗口函数和该命令之后围绕每个整合位点的1兆区的基因密度:
   bedtools窗口-w 500000 central_basepair_1.bed -b RefSeq_hg38_NonRedundant.bed -u> GeneDensity_sample1.bed
9. 通过下面这个命令计算数据集中所有集成的平均密度基因:
   awk的'（总和+ = $ 7）END（打印"平均="总和/ NR）"GeneDensity_sample1.bed

10.统计比较整合位点分布样品采用双尾中的R Fisher精确检验和双尾秩和检验中

注意:使用Fisher精确测试的RefSeq基因内或CpG岛或起始位置跟的窗口内比​​较整合位点的比例，但使用Wilcoxon秩和检验的分布在周围的整合位点基因密度比较。 R程序可在http://www.r-project.org/。
双尾Fisher精确检验:

1. 使用如步骤9.4和9.7指示计算出的数字，铬eate矩阵中观察到的事件与按照该命令余下的网站（注释内或周围的注释窗口内集成）的R各自的比较:
   （annotation_of_interest < - 矩阵（C（SampleA＃中，SampleA＃剩余SampleB＃中，SampleB＃剩余），nrow = 2，dimnames =名单（C（'中心'，'余'），C（'SampleA'， 'SampleB'））））
2. 计算双尾费舍尔的使用以下命令精确检验比较P值:
   fisher.test（annotation_of_interest，另类='two.sided'）$ p.value
   双尾秩和检验:
3. 创建一个制表符分隔的.txt文件，其中每列包含在顶部电池的样品名称，由该库（从9.9步产生的.bed文件中获得的）所有整合位点基因密度值低于紧随其后。使用以下命令导入这个制表符分隔的.txt文件到R和NAvigating到正确的文件目录:
   FILENAME < - as.data.frame（read.delim（file.choose（），标题= T，check.names = FALSE，填写= TRUE，九月=' t'））
4. 计算通过使用以下命令双尾秩和检验比较P值:
   wilcox.test（文件名$ SampleA，FILENAME $ SampleB，另类='two.sided'，配对= F，准确= T）$ p.value
   注:P值只能下降到R中一定（非常低）极限来计算，在此之后零将由程序返回。对于大量不同的样品能产生在R A P = 0，估计P值作为<2.2×10 ^-308。

11.检查原始测序数据对于包含集成病毒DNA细胞的克隆膨胀的迹象

注意:存在在参考基因组完全相同的核苷酸多个结合一个小的潜力。可替代地，一个单一的tegration事件可能成为在测序数据冗余地存在由于文库制备期间和/或之前，通过DNA制备细胞复制的使用PCR。从HIV感染患者的基因组DNA的最近的分析已经通过识别独特的超声波剪切点/连接器的连接点包含相同的整合位点^52-54的DNA序列内（可前PCR才会出现）来区分这些可能性。目前，以克隆扩增的细胞内窝藏原病毒是否向潜病毒储，因而它是特别令人感兴趣的研究在人类患者中的整合位点时以表征其膨胀水​​平的辩论。

1. 类似于在步骤8.1中所列的程序，生成.bed文件列出的延伸碱基的间隔，在这种情况下，25个核苷酸的下游从每个唯一整合位点（上游的碱是不必要这里）。产生这些.bed文件FASTA文件（如指示步骤8.1），通过使用来自BEDTools的fastaFromBed功能，下面这个命令:
   fastaFromBed -fi /目录/至/参考/基因组/ -name -s -bed 25_base_pair_file.bed -fo output_file.fasta
   注意:为了提高每个的特异性查询，建议从各积分位点提取至少25个核苷酸的下游用于克隆扩增的分析。
2. 最好使用自定义脚本，搜索包含完全匹配的25个核苷酸下游每一个独特的整合位点的所有字符串原始序列数据FASTA文件，并存放这些序列到一个新文件。从原始的字符串修剪LTR和接头序列。合并的PE序列通过转换读取读取到的反向互补，修整LTR和接头序列，然后分配READ2字符串其READ1对如果字符串共享至少20个核苷酸的重叠。
3. 扫描每个整合位点块的连接器连接点。每个分类整合为"克隆扩增＆＃34;如果接头连接点≥3沸点分开。
   注意:不合并序列克隆扩增分析的协议内容已被描述^52。
   注意:在通过超声处理的完全相同位置的基因组中的碎片导致克隆扩增的程度的低估，和方法来纠正所产生的实验偏压已经描述^63,64。

结果

表4列出一个代表性实验的结果，以示出用于从感染细胞的培养物中回收整合位点的NGS的灵敏度。未感染的细胞的DNA被用于串联从其中在平均每个细胞包含一个积分⁴⁰的感染稀释的基因组DNA。 15625:在五个步骤的1的最大稀释制备稀释液。在滴定系列的基因组DNA，然后通过超声处理片段化或通过用限制性内切核酸酶消化MseI和BglⅡ，随后的LM-PCR的独特整合位点的数目，以及网站地图近端到选定的基因组注释的数目，是根据上述方案计算。数据分析显示，数十个独特的集成网站（量的1-2％，从整齐的基因组DNA恢复）从哪里理论上只有15625人被感染的细胞制备库恢复。

当分析整合位数据集，这是至关重要的数据进行比较，以匹配组随机基因组位点，这被称为匹配的随机控制或MRC。作为代表性的结果通过限制性内切酶消化或超声剪切基因组DNA，两个不同的数据集MRC构建。 MRC _恩茨包含从hg19随机选择地点邻近MseI和BglⅡ酶切的部位产生50000独特的基因组位点，而MRC _随机窝藏1...

讨论

一种用于逆转录病毒整合位点的分析协议，从通过基因组分布模式映射的初始病毒感染步骤，进行说明。这个协议是适用于任何反转录病毒和任何可感染的细胞类型。此外，该测定管道是相当敏感，以回收的从基因组DNA等效的系列稀释独特的整合位点令人满意编号以与6.4×10 ^-5的MOI启动的感染的可能性。当从可能含有低病毒载量，其中只有细胞的一小部分会窝藏集成原病​​毒感染的患者施加到样本这种敏感性使得协议特别有用。在这一领域^36,38,41-43之前方法论的论文一致，在这个协议的生物信息学部分的多个步骤将定制的脚本开发受益顺序处理大量数据文件。虽然BLAT ⁵⁸是马pping在本协议中所述程序，用户可能会发现蝴蝶结^{67（http://bowtie-bio.sourceforge.net/index.shtml）}是一个合适的选择。

一个替代的生物信息学流水线最近报道的测定莫洛尼鼠白血病病毒（的MoMLV）整合位点^19。在它被开发成独立的软件，是可公开获得的，并且是相当强大的，因为它最初是用于映...

材料

本文使用的材料清单

姓名	公司	目录编号	评论
DMEM	Gibco	11965-084	标准细胞培养基，与HEK293T细胞兼容
胎牛血清	Thermo Scientific	SH 30088.03	可能需要预先筛选不同批次的血清以实现最佳病毒生产
青霉素/链霉素	Corning	30-002-Cl	抗生素将添加到 DMEM
磷酸盐缓冲盐水中	Mediatech	21-040-CV	用于洗涤细胞
胰蛋白酶 EDTA	康宁	25-053-CI	用于从组织培养板中分离贴壁细胞
PolyJet	SignaGen Laboratories	SL100688	DNA 转染试剂
0.45 &微量;m 过滤器	Thermo Scientific	09-740-35B	用于过滤含有病毒颗粒的细胞培养基
Turbo DNase	Ambion	AM2239	用于降解病毒原液中的残留质粒
DNA HIV-1 p24 抗原捕获测定	ABL Inc.	5447	用于定量病毒产量
DNeasy Blood 和组织试剂盒	Qiagen	69506	用于从细胞中纯化基因组DNA
超声仪	Covaris	S2	使用这种型号的超声仪执行两轮占空比，5%;强度，3;每次爆发循环，200;时间，80 秒
无核酸酶水	GeneMate	G-3250-125	建议使用市售水，以减少样品交叉污染的可能性
QIAQuick PCR 纯化试剂盒	Qiagen	28106	用于在文库构建过程中纯化
DNA End-It DNA End-Repair Kit	Epicentre	ER81050	用于修复超声处理的 DNA 样品的 DNA 末端
Klenow 片段 （3'-5' exo–）	New England Biolabs （NEB）	M0212S	与 dATP 一起用于 A 尾修复的 DNA 片段
dATP	Thermo Scientific	R0141	三磷酸脱氧腺苷
MseI	NEB	R0525L	用于基因组 DNA 切割的限制性核酸内切酶
BglII	NEB	R0144L	内切酶抑制上游 HIV-1 U5 序列
T4 DNA 连接酶的扩增	NEB	M0202L/6218	用于相容 DNA 末端共价连接的酶
DNA 寡核苷酸	集成 DNA 技术	定制	让公司纯化寡核苷酸。HPLC 纯化足以检测 DNA <30 个核苷酸;PAGE 纯化较长的 DNA
优势 2 聚合酶混合物	Clontech	639202	含有 DNA 聚合酶的商业混合物，用于 PCR
dNTP（100 mM 溶液）Thermo	Scientific	R0181	用无菌水在冰上稀释四种化学品，以达到每个 dNTP 2.5 mM 的中等浓度
NanoDrop Thermo	Scientific	NanoDrop 2000	分光光度计，用于测定 DNA 浓度
Qubit 荧光计	Life Technologies	Qubit® 3.0	荧光计，用于确认整合位点文库 DNA 浓度
2200 TapeStation 系统	Agilent	G2964AA	基于磁带的分析，用于确认整合位点文库 DNA 大小分布
MiSeq	Illumina	SY-410-1003	用于 NGS