方法文章

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

DOI:

10.3791/53840

2016年3月22日

本文内容

摘要

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

摘要

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

在本地和全球尺度上表现出逆转录病毒整合签名的偏好。这里,我们提出了(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在三天构成。该包括易感组织培养细胞至整合位点分析的外源病毒感染的整个协议可以因此在大约一到两周进行。这项技术的应用近期涉及到艾滋病毒感染患者的整合位点的纵向分析。

引言

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

病毒DNA(VDNA)到宿主细胞基因组中的整合是在逆转录病毒的生命周期的一个重要步骤。积分是由病毒酶整合酶(IN),其执行,导致建立稳定地插入原病毒1的两个不同的催化过程来完成的。在亚基接合是通过逆转录生成的线性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,有利于积极转录基因17机构内部整合,而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扩增与映射单个整合位点,以转化领域人类和小鼠草案基因组,与来自外源组织培养细胞的感染中回收网的增加至几百到上千17数量的能力的结合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的近似等效惯性矩。

访问受限。请登录或开始试用以查看此内容。

方案

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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 2的潮湿的细胞培养箱中,在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 2的组织培养烧瓶中,并孵育所述细胞在组织培养培养箱中。
  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)确认用的基于磁带的测定法56的DNA片段大小分布。
    注:理想的分布是比较宽泛的DNA峰中心在长约500基点。如果材料的显著量大于1 kb的大,则建议掺入大小选过程,以消除更长的DNA的物种,这将聚类过程中妨碍桥扩增。相反,如果一个显著峰是明显的约100至200碱基,引物二聚体可能PCR过程中形成的。在这种情况下,该过程应当被优化,以尽量减少引物二聚体的形成。
  3. (QC检测#3)有限公司通过定量PCR 57并将通过适配器插入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 58的病毒基因组,与输出整合位点坐标输出到一个单独的.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 59 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/database)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。

访问受限。请登录或开始试用以查看此内容。

结果

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

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

访问受限。请登录或开始试用以查看此内容。

讨论

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

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

访问受限。请登录或开始试用以查看此内容。

披露

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

作者没有什么可透露的。

致谢

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

我们非常感谢我们的同事斯蒂芬·休斯和亨利·莱文的意见,这是建立逆转录病毒整合位点的测序NGS协议在恩格尔曼实验的关键。这项工作是由美国国立卫生研究院授予AI039394和​​AI052014(以ANE)和AI060354(哈佛大学艾滋病研究中心)的支持。

访问受限。请登录或开始试用以查看此内容。

材料

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

参考文献

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Craigie, R., Bushman, F. D. HIV DNA integration. Cold Spring Harb. Perspect. Med. 2, a006890(2012).
  2. Li, M., Mizuuchi, M., Burke, T. R. J., Craigie, R. Retroviral DNA integration: reaction pathway and critical intermediates. EMBO J. 25, 1295-1304 (2006).
  3. Hare, S., Gupta, S. S., Valkov, E., Engelman, A., Cherepanov, P. Retroviral intasome assembly and inhibition of DNA strand transfer. Nature. 464, 232-236 (2010).
  4. Hare, S., Maertens, G. N., Cherepanov, P. 3'-processing and strand transfer catalysed by retroviral integrase in crystallo. EMBO J. 31, 3020-3028 (2012).
  5. Fujiwara, T., Mizuuchi, K. Retroviral DNA integration: structure of an integration intermediate. Cell. 54, 497-504 (1988).
  6. Roth, M. J., Schwartzberg, P. L., Goff, S. P. Structure of the termini of DNA intermediates in the integration of retroviral DNA: dependence on IN function and terminal DNA sequence. Cell. 58, 47-54 (1989).
  7. Brown, P. O., Bowerman, B., Varmus, H. E., Bishop, J. M. Retroviral integration: structure of the initial covalent product and its precursor, and a role for the viral IN protein. Proc. Natl. Acad. Sci. USA. 86, 2525-2529 (1989).
  8. Pauza, C. D. Two bases are deleted from the termini of HIV-1 linear DNA during integrative recombination. Virology. 179, 886-889 (1990).
  9. Bowerman, B., Brown, P. O., Bishop, J. M., Varmus, H. E. A nucleoprotein complex mediates the integration of retroviral DNA. Genes Dev. 3, 469-478 (1989).
  10. Bukrinsky, M. I., et al. Active nuclear import of human immunodeficiency virus type 1 preintegration complexes. Proc. Natl. Acad. Sci. USA. 89, 6580-6584 (1992).
  11. Miller, M. D., Farnet, C. M., Bushman, F. D. Human immunodeficiency virus type 1 preintegration complexes: studies of organization and composition. J. Virol. 71, 5382-5390 (1997).
  12. Engelman, A., Mizuuchi, K., Craigie, R. HIV-1 DNA integration: mechanism of viral DNA cleavage and DNA strand transfer. Cell. 67, 1211-1221 (1991).
  13. Maertens, G. N., Hare, S., Cherepanov, P. The mechanism of retroviral integration from X-ray structures of its key intermediates. Nature. 468, 326-329 (2010).
  14. Holman, A. G., Coffin, J. M. Symmetrical base preferences surrounding HIV-1, avian sarcoma/leukosis virus, and murine leukemia virus integration sites. Proc. Natl. Acad. Sci. USA. 102, 6103-6107 (2005).
  15. Wu, X., Li, Y., Crise, B., Burgess, S. M., Munroe, D. J. Weak palindromic consensus sequences are a common feature found at the integration target sites of many retroviruses. J. Virol. 79, 5211-5214 (2005).
  16. Kvaratskhelia, M., Sharma, A., Larue, R. C., Serrao, E., Engelman, A. Molecular mechanisms of retroviral integration site selection. Nucleic Acids Res. 42, 10209-10225 (2014).
  17. Schroder, A. R., et al. HIV-1 integration in the human genome favors active genes and local hotspots. Cell. 110, 521-529 (2002).
  18. Wu, X., Li, Y., Crise, B., Burgess, S. M. Transcription start regions in the human genome are favored targets for MLV integration. Science. 300, 1749-1751 (2003).
  19. LaFave, M. C., et al. MLV integration site selection is driven by strong enhancers and active promoters. Nucleic Acids Res. 42, 4257-4269 (2014).
  20. De Ravin, S. S., et al. Enhancers are major targets for murine leukemia virus vector integration. J. Virol. 88, 4504-4513 (2014).
  21. Maskell, D. P., et al. Structural basis for retroviral integration into nucleosomes. Nature. 523, 366-369 (2015).
  22. Serrao, E., et al. Integrase residues that determine nucleotide preferences at sites of HIV-1 integration: implications for the mechanism of target DNA binding. Nucleic Acids Res. 42, 5164-5176 (2014).
  23. Aiyer, S., et al. Structural and sequencing analysis of local target DNA recognition by MLV integrase. Nucleic Acids Res. 43, 5647-5663 (2015).
  24. Ciuffi, A., et al. A role for LEDGF/p75 in targeting HIV DNA integration. Nat. Med. 11, 1287-1289 (2005).
  25. Sharma, A., et al. BET proteins promote efficient murine leukemia virus integration at transcription start sites. Proc. Natl. Acad. Sci. USA. 110, 12036-12041 (2013).
  26. Gupta, S. S., et al. Bromo- and extraterminal domain chromatin regulators serve as cofactors for murine leukemia virus integration. J. Virol. 87, 12721-12736 (2013).
  27. De Rijck, J., et al. The BET family of proteins targets moloney murine leukemia virus integration near transcription start sites. Cell Rep. 5, 886-894 (2013).
  28. Demeulemeester, J., et al. HIV-1 integrase variants retarget viral integration and are associated with disease progression in a chronic infection cohort. Cell Host Microbe. 16, 651-662 (2014).
  29. Meehan, A. M., et al. LEDGF/p75 proteins with alternative chromatin tethers are functional HIV-1 cofactors. PLoS Pathog. 5, e1000522(2009).
  30. Ferris, A. L., et al. Lens epithelium-derived growth factor fusion proteins redirect HIV-1 DNA integration. Proc. Natl. Acad. Sci. USA. 107, 3135-3140 (2010).
  31. Gijsbers, R., et al. LEDGF hybrids efficiently retarget lentiviral integration into heterochromatin. Mol. Ther. 18, 552-560 (2010).
  32. Aiyer, S., et al. Altering murine leukemia virus integration through disruption of the integrase and BET protein family interaction. Nucleic Acids Res. 42, 5917-5928 (2014).
  33. Jahner, D., Jaenisch, R. Integration of Moloney leukaemia virus into the germ line of mice: correlation between site of integration and virus activation. Nature. 287, 456-458 (1980).
  34. Stevens, S. W., Griffith, J. D. Human immunodeficiency virus type 1 may preferentially integrate into chromatin occupied by L1Hs repetitive elements. Proc. Natl. Acad. Sci. USA. 91, 5557-5561 (1994).
  35. Wang, G. P., Ciuffi, A., Leipzig, J., Berry, C. C., Bushman, F. D. HIV integration site selection: analysis by massively parallel pyrosequencing reveals association with epigenetic modifications. Genome Res. 17, 1186-1194 (2007).
  36. Wang, G. P., et al. DNA bar coding and pyrosequencing to analyze adverse events in therapeutic gene transfer. Nucleic Acids Res. 36, e49(2008).
  37. Roth, S. L., Malani, N., Bushman, F. D. Gammaretroviral integration into nucleosomal target DNA in vivo. J. Virol. 85, 7393-7401 (2011).
  38. Ciuffi, A., Barr, S. D. Identification of HIV integration sites in infected host genomic DNA. Methods. 53, 39-46 (2011).
  39. Gillet, N. A., et al. The host genomic environment of the provirus determines the abundance of HTLV-1-infected T-cell clones. Blood. 117, 3113-3122 (2011).
  40. Matreyek, K. A., et al. Host and viral determinants for MxB restriction of HIV-1 infection. Retrovirology. 11, 90(2014).
  41. Ciuffi, A., et al. Methods for integration site distribution analyses in animal cell genomes. Methods. 47, 261-268 (2009).
  42. Brady, T., et al. A method to sequence and quantify DNA integration for monitoring outcome in gene therapy. Nucleic Acids Res. 39, e72(2011).
  43. Beard, B. C., Adair, J. E., Trobridge, G. D., Kiem, H. P. High-throughput genomic mapping of vector integration sites in gene therapy studies. Methods Mol. Biol. 1185, 321-344 (2014).
  44. Page, K. A., Landau, N. R., Littman, D. R. Construction and use of a human immunodeficiency virus vector for analysis of virus infectivity. J. Virol. 64, 5270-5276 (1990).
  45. Emi, N., Friedmann, T., Yee, J. K. Pseudotype formation of murine leukemia virus with the G protein of vesicular stomatitis virus. J. Virol. 65, 1202-1207 (1991).
  46. Wehrly, K., Chesebro, B. p24 antigen capture assay for quantification of human immunodeficiency virus using readily available inexpensive reagents. Methods. 12, 288-293 (1997).
  47. Goff, S., Traktman, P., Baltimore, D. Isolation and properties of Moloney murine leukemia virus mutants: use of a rapid assay for release of virion reverse transcriptase. J. Virol. 38, 239-248 (1981).
  48. Willey, R. L., et al. In vitro mutagenesis identifies a region within the envelope gene of the human immunodeficiency virus that is critical for infectivity. J. Virol. 62, 139-147 (1988).
  49. Butler, S. L., Hansen, M. S., Bushman, F. D. A quantitative assay for HIV DNA integration in vivo. Nat. Med. 7, 631-634 (2001).
  50. Brussel, A., Sonigo, P. Analysis of early human immunodeficiency virus type 1 DNA synthesis by use of a new sensitive assay for quantifying integrated provirus. J. Virol. 77, 10119-10124 (2003).
  51. Serrao, E., Ballandras-Colas, A., Cherepanov, P., Maertens, G. N., Engelman, A. N. Key determinants of target DNA recognition by retroviral intasomes. Retrovirology. 12, 39(2015).
  52. Maldarelli, F., et al. HIV latency. Specific HIV integration sites are linked to clonal expansion and persistence of infected cells. Science. 345, 179-183 (2014).
  53. Wagner, T. A., et al. HIV latency. Proliferation of cells with HIV integrated into cancer genes contributes to persistent infection. Science. 345, 570-573 (2014).
  54. Cohn, L. B., et al. HIV-1 integration landscape during latent and active infection. Cell. 160, 420-432 (2015).
  55. Li, X., Ben-Dov, I. Z., Mauro, M., Williams, Z. Lowering the quantification limit of the QubitTM RNA HS assay using RNA spike-in. BMC Mol. Biol. 16, 9(2015).
  56. Padmanaban, A., Walker, D. M. Analysis of high molecular weight genomic DNA using the Agilent 2200 TapeStation and genomic DNA ScreenTape. Publication number 5991-1797EN Agilent Technologies. , Agilent Technologies, Inc. Santa Clara, CA. (2013).
  57. Kapa library quantification technical guide version v1.14. , KapaBiosystems. Boston, MA. (2014).
  58. Kent, W. J. BLAT--the BLAST-like alignment tool. Genome Res. 12, 656-664 (2002).
  59. Quinlan, A. R., Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 26, 841-842 (2010).
  60. Gaur, M., Leavitt, A. D. Mutations in the human immunodeficiency virus type 1 integrase D,D(35)E motif do not eliminate provirus formation. J. Virol. 72, 4678-4685 (1998).
  61. Varadarajan, J., McWilliams, M. J., Hughes, S. H. Treatment with suboptimal doses of raltegravir leads to aberrant HIV-1 integrations. Proc. Natl. Acad. Sci. USA. 110, 14747-14752 (2013).
  62. Kent, W. J., et al. The human genome browser at UCSC. Genome Res. 12, 996-1006 (2002).
  63. Berry, C. C., et al. Estimating abundances of retroviral insertion sites from DNA fragment length data. Bioinformatics. 28, 755-762 (2012).
  64. Firouzi, S., et al. Development and validation of a new high-throughput method to investigate the clonality of HTLV-1-infected cells based on provirus integration sites. Genomic Med. 6, 46(2014).
  65. Mitchell, R. S., et al. Retroviral DNA integration: ASLV, HIV, and MLV show distinct target site preferences. PLoS Biol. 2, E234(2004).
  66. Schneider, T. D., Stormo, G. D., Gold, L., Ehrenfeucht, A. Information content of binding sites on nucleotide sequences. J. Mol. Biol. 188, 415-431 (1986).
  67. Langmead, B. Chapter 11, Unit 11. 17, Aligning sort sequencing reads with Bowtie. Curr. Protoc. Bioinformatics. , (2010).
  68. LaFave, M. C., Varshney, G. K., Burgess, S. M. GeIST: a pipeline for mapping integrated DNA elements. Bioinformatics. 31, 3219-3221 (2015).
  69. Hocum, J. D., et al. VISA - Vector Integration Site Analysis server: a web-based server to rapidly identify retroviral integration sites from next-generation sequencing. BMC Bioinformatics. 16, 212(2015).
  70. Gabriel, R., et al. Comprehensive genomic access to vector integration in clinical gene therapy. Nat. Med. 15, 1431-1436 (2009).
  71. TruSeq Library Prep Pooling Guide. Guidelines for pooling TruSeq libraries for Illumina sequencing systems that require balanced index combinations. , Source: https://support.illumina.com/downloads/truseq-library-prep-pooling-guide-15042173.html (2015).

访问受限。请登录或开始试用以查看此内容。

重印与许可

申请许可以重复使用本 JoVE 文章的文本或图表

申请许可

标签

PCR DNA PCR BLAT BEDTools

相关文章