CN104561013A - 基于高通量测序技术进行核酸适配体序列优化的方法 - Google Patents
基于高通量测序技术进行核酸适配体序列优化的方法 Download PDFInfo
- Publication number
- CN104561013A CN104561013A CN201510002930.5A CN201510002930A CN104561013A CN 104561013 A CN104561013 A CN 104561013A CN 201510002930 A CN201510002930 A CN 201510002930A CN 104561013 A CN104561013 A CN 104561013A
- Authority
- CN
- China
- Prior art keywords
- sequence
- parsing
- dna
- aptamer
- dna fragmentation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于高通量测序技术进行核酸适配体序列优化的方法,包括将核酸适配体以单个碱基顺移的方式解析后,合成带有接头DNA片段的解析DNA片段;将解析DNA片段与靶标孵育,获得其中与靶标结合的部分;以结合解析DNA片段为模板,采用逆转录酶催化获得双链DNA文库;采用高通量测序技术对获得的双链DNA文库进行测序,然后根据测序结果,选择百分比最高的解析DNA片段作为核酸适配体的优化序列。本发明比普通的核酸适配体方法更准确、快速,可以实现多种核酸适配体、多种靶标类型的同步、高通量分析,提高了核酸适配体优化研究的效率。
Description
技术领域
本发明属于分子生物学领域,更具体地说,涉及一种基于高通量测序技术进行核酸适配体序列优化的方法。
背景技术
核酸适配体是能通过特定三维结构与多种靶标紧密结合的寡核苷酸链(DNA或RNA)。核酸适配体具有亲和力高,特异性好,分子量小,免疫原性低,结构稳定,易合成等优点,因此核酸适配体替代传统抗体在多种疾病的基础研究、药物筛选、临床诊断及疾病治疗中具有广泛的应用前景。
通常,核酸适配体是通过指数富集配体的系统进化技术(systematic evolution of ligands by exponential enrichment, SELEX) 从含有大量序列的随机文库中筛选得到的。随机文库的基本序列结构包括:中间长度为35-60mer的随机序列,两端各为约20mer的引物结合序列。筛选到的核酸适配体同样具有上述序列结构,总长度约60-100mer。在已报道研究中发现核酸适配体长度过长不利于形成稳定的结构,还具有潜在的免疫激活作用。而通过去除核酸适配体与靶标结合的非必需序列,仅保留必需序列可以使核酸适配体长度缩短,并提高核酸适配体与靶标结合的亲和力和特异性,便于后续研究,同时也节省了成本,提高了安全性。目前,应用广泛的核酸适配体,例如凝血酶适配体、ATP适配体、IgE适配体等,均经过序列优化,长度在15-25mer水平。对于筛选到的核酸适配体进行序列优化,已成为其进一步应用前所必须的。
对核酸适配体进行序列优化的方法包括生物信息学预测法、酶切保护实验、芯片分析法等。生物信息学预测法是基于热力学最小自由能等算法预测核酸适配体的二级结构,但由于各算法可预测到的核酸二级结构十分有限,其可靠性难以保证;酶切保护实验是利用核酶水解核酸适配体与靶标结合的二聚体,这样可以得知核酸适配体与靶标结合的部分序列,但该方法无法得知核酸适配体中对稳定性起关键作用的非结合序列;芯片技术是将核酸适配体的截短片段制成芯片,利用荧光修饰靶标与芯片的结合信号,优选亲和力更好的优化序列,该方法准确性高,但成本过于昂贵。由于目前对于核酸适配体进行序列优化尚缺乏快速、准确、令人满意的方法,大量筛选到的核酸适配体因为缺乏有效优化而序列过长,从而限制了其广泛应用。
发明内容
本发明所要解决的技术问题是克服现有技术中种种不足之处,提供一种基于高通量测序技术进行核酸适配体序列优化的方法。
本发明提供的技术方案是:
一种基于高通量测序技术进行核酸适配体序列优化的方法,其特征在于所属方法包括以下步骤:
(1) 将候选核酸适配体以单个碱基顺移的方式分解为解析序列;
(2) 合成解析DNA片段,所述解析DNA片段包含解析序列,一端与DNA接头片段相连,所述一端为解析序列的3’端或5’端;
(3) 将解析DNA片段按等摩尔混合后制备成解析DNA文库,将解析DNA文库与靶标的混合物,在一定条件下孵育。所述靶标为筛选核酸适配体所采用的物质;所述一定条件为核酸适配体筛选时的孵育条件,包括孵育的温度、时间以及缓冲液的类型、pH值;
(4) 将结合解析DNA片段与游离解析DNA片段分离,所述结合解析DNA片段为与靶标结合的解析DNA片段;
(5) 将结合解析DNA片段与靶标分离,制备结合解析DNA文库;
(6) 以结合解析DNA文库为模板,利用引物和逆转录酶合成结合解析DNA片段的互补链,制备双链DNA文库。所述的引物与所述解析DNA片段中DNA接头序列互补;
(7) 采用高通量测序技术对获得的双链DNA文库进行高通量测序,然后对测序结果进行分析,双链DNA文库中序列最多的解析序列为核酸适配体的优化序列。
优选的,所述方法中所述候选核酸适配体为一种序列或多种序列;
优选的,所述方法中所述解析序列长度范围为15-30碱基,长度范围也可根据候选核酸适配体的不同而进行调整;
优选的,所述方法中所述DNA接头片段为不同长度的多聚A序列、多聚T序列;当候选核酸适配体为多种序列时,DNA接头片段的序列和长度可对应不同。
优选的,所述方法中所述靶标为小分子、大分子、细胞及病原体。
优选的,所述方法中所述逆转录酶为M-MLV逆转录酶、AMV逆转录酶。
优选的,所述步骤(1)中进行候选核酸适配体的解析是利用计算机软件(Sequence Split Tool)来完成。
与现有技术相比,本发明首次利用高通量测序原理完成适配体核酸的序列优化。该方法以核酸适配体一级序列为基础,避免了生物信息学分析的不确定性,序列优化的准确性提高显著;通芯片方法相比,成本显著降低,无需标价靶标,可用于任何类型靶标的核酸适配体。针对不同靶标或不同核酸适配体,设计不同的DNA接头片段加以区分,这样使得一次高通量测序,可同时对多种候选核酸适配体进行优化。
附图说明
下面结合附图和实施例对本发明作进一步描述:
图1为本发明进行核酸适配体序列优化的方法流程示意图;
图2为本发明解析DNA片段的结构示意图;在序列结构组成上,所属解析DNA片段由接头DNA片段+解析序列组成。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图激具体实施例对上述方案做进一步说明。应理解,该具体实施例是用于说明本发明而不限于限制本发明的范围。实施例中采用的实施条件可以根据候选核酸适配体及靶标的不同做进一步调整,未注明的实施条件通常为常规实验中的条件。
实施例半乳甘露聚糖(GM)核酸适配体的序列优化方法示例
请参阅图1,为根据本发明的优选事实例的基于高通量测序技术进行核酸适配体序列优化的方法流程示意图。主要包括以下步骤:
一.GM核酸适配体的解析
在步骤①中,以GM为靶标筛选到特异性核酸适配体AT74,长度为62个碱基,序列如下:AT74: TCCCTACGGCGCTAACTGGAGTGAGGAACCCGCCATACCGATGTCTGCC- ACCGTGCTACAAC;
利用Sequence Split Tool软件平台(http://genestamp.sinica.edu.tw/marray/fun/ seqSplit.htm)将AT74序列解析为15-30mer的解析序列。在Minimum length栏中输入15,在Maximum length栏中输入30,在Input sequence栏中输入所述AT74的碱基序列。点击“提交查询内容”,AT74被解析为648条解析序列。
二.核酸适配体解析DNA片段的合成
请参阅图2,解析DNA片段的序列为在每条AT74解析序列的3’端加入polyT(n=15)序列。所有解析DNA片段委托生工生物工程(上海)股份有限公司合成并通过HPLC纯化。
三.解析DNA文库与GM结合
1. 将合成的解析DNA片段按照4 μmol/L的浓度溶解于Binding buffer(0.05 mol/L Tris-HCl, 0.1 mol/L NaCl, 1mmol/LMgCl2, 5mmol/L KCl,pH=7.4)。
2. 取等份的解析DNA片段混合成解析DNA文库,将解析DNA文库在95℃加热5min,然后放到冰上迅速冷却。
3. 在步骤②中,取解析DNA文库500 μl与GM纯品混合,混合溶液在37℃条件下,震动混合孵育30min,以使DNA与GM能够充分结合。
四.结合解析DNA片段的分离
在步骤③中,将混合溶液转移至10 kDa的超滤离心柱中,12000g 离心5min。超滤离心柱中GM由于分子量大于10 kDa,结合解析DNA片段随着GM被留下来,并于游离解析DNA片段分离。用500 μl Binding Buffer超滤洗涤3次。
五.结合解析DNA文库的制备
1.将结合了结合解析DNA片段的GM中加入200 μl灭菌蒸馏水,100℃加热5min;
2.在步骤④中,将去离子水转移至10 kDa的超滤离心柱中,12000g 离心5min,留取上清。
六.双链DNA文库的制备
在步骤⑤中,以结合解析DNA文库为模板,oligo (dA)15为引物,42℃反应30min,合成双链DNA文库,反应体系如下:
Temple | 100 μl |
oligo (dA)15 | 5 μl (100pmol) |
10 mM dNTP Mix | 5 μl |
5ⅹRT Buffer | 30 μl |
AMV | 5 μl |
灭菌蒸馏水 | Up to 150 μl |
七.高通量测序分析
在步骤⑥中,高通量测序分析双链DNA文库中的序列组成,采用ABI SOLiD连接法测序,委托商业公司完成。测序分析表明, 648个解析片段中的7个与GM具有较高的结合率,最高可达36%;提示核心序列“CGCTAACTGGAGTGAGGAACCCG”可能是核酸适配体AT74的作用位点,而两边侧翼序列起到稳定结构的作用,统计结果见下表:
编号 | 序列 | 序列组成(%) |
30.7 | ACGGCGCTAACTGGAGTGAGGAACCCGCCA- | 36 |
25.10 | --GGCGCTAACTGGAGTGAGGAACCCG--- | 27 |
26.10 | ---GCGCTAACTGGAGTGAGGAACCCGCC-- | 18 |
26.7 | ----CGCTAACTGGAGTGAGGAACCCGCCA- | 9 |
26.9 | -CGGCGCTAACTGGAGTGAGGAACCCG--- | 7 |
29.8 | -CGGCGCTAACTGGAGTGAGGAACCCGCCA- | 1.5 |
28.6 | ---GCGCTAACTGGAGTGAGGAACCCGCCAT | 0.6 |
八.GM亲和力验证
进一步比较核酸适配体核心序列(AT74C)、碱基突变适配体(AT74MT)和完整DNA适配体(AT74)与GM的亲和力,发现AT74CAT74和与GM的亲和力相当,均显著高于AT74MT,经优化,核酸适配体在保持亲和力不变的情况下,长度被显著缩短。
Claims (6)
1.一种基于高通量测序技术进行核酸适配体序列优化的方法,其特征在于所述方法包括以下步骤:
(1) 将候选核酸适配体以单个碱基顺移的方式分解为解析序列;
(2) 合成解析DNA片段,所述解析DNA片段包含解析序列,一端与DNA接头相连,所述一端为解析序列的3’端或5’端;
(3) 将解析DNA片段按等摩尔混合后制备成解析DNA文库,将解析DNA文库与靶标的混合物,在一定条件下孵育,所述靶标为筛选适配体所采用的物质;所述一定条件为核酸适配体筛选时的孵育条件,包括孵育的温度、时间以及缓冲液的类型、pH值;
(4) 将结合解析DNA片段与游离解析DNA片段分离,所述结合解析DNA片段为与靶标结合的解析DNA片段;
(5) 将结合解析DNA片段与靶标分离,制备结合解析DNA文库;
(6) 以结合解析DNA文库为模板,利用引物和逆转录酶合成结合解析DNA片段的互补链,制备双链DNA文库,所述的引物与所述解析DNA片段中DNA接头序列互补;
(7) 采用高通量测序技术对获得的双链DNA文库进行高通量测序,然后对测序结果进行分析,双链DNA文库中序列最多的解析序列为核酸适配体的优化序列。
2.根据权利要求1所述的方法,其特征在于所述方法中所述解析序列长度范围为15-30碱基。
3.根据权利要求1所述的方法,其特征在于所述方法中所述DNA接头为多聚A序列、多聚T序列。
4.根据权利要求1所述的方法,其特征在于所述方法中所述靶标为小分子、大分子、细胞及病原体。
5.根据权利要求1所述的方法,其特征在于所述方法中所述逆转录酶为M-MLV逆转录酶、AMV逆转录酶。
6.根据权利要求1所述的方法,其特征在于所述方法步骤(1)中进行候选核酸适配体的解析是利用计算机软件(Sequence Split Tool)来完成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510002930.5A CN104561013A (zh) | 2015-01-05 | 2015-01-05 | 基于高通量测序技术进行核酸适配体序列优化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510002930.5A CN104561013A (zh) | 2015-01-05 | 2015-01-05 | 基于高通量测序技术进行核酸适配体序列优化的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104561013A true CN104561013A (zh) | 2015-04-29 |
Family
ID=53078092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510002930.5A Pending CN104561013A (zh) | 2015-01-05 | 2015-01-05 | 基于高通量测序技术进行核酸适配体序列优化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104561013A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017133159A1 (zh) * | 2016-02-03 | 2017-08-10 | 中国农业科学院北京畜牧兽医研究所 | 一种计算机辅助筛选小分子化合物靶标适配体的实现方法 |
CN108387561A (zh) * | 2018-01-18 | 2018-08-10 | 大连理工大学 | 一种基于碱基淬灭荧光原理实现低成本高通量核酸适配体最优序列的优化方法 |
CN108823211A (zh) * | 2018-05-04 | 2018-11-16 | 中国医学科学院基础医学研究所 | 一种cd19核酸适配体及其应用 |
CN109868268A (zh) * | 2017-12-04 | 2019-06-11 | 上海药明康德新药开发有限公司 | 一种优化dna编码化合物文库起始片段的方法 |
CN111445951A (zh) * | 2020-03-17 | 2020-07-24 | 大连理工大学 | 一种基于免标记杂交探针竞争法实现低成本高通量核酸适配体最优序列的优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104069512A (zh) * | 2010-04-12 | 2014-10-01 | 私募蛋白质体公司 | 针对β-NGF的适体及其在治疗β-NGF介导的疾病和失调中的用途 |
WO2014159669A2 (en) * | 2013-03-14 | 2014-10-02 | Somalogic, Inc. | Aptamers that bind to il-6 and their use in treating or diagnosing il-6 mediated conditions |
-
2015
- 2015-01-05 CN CN201510002930.5A patent/CN104561013A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104069512A (zh) * | 2010-04-12 | 2014-10-01 | 私募蛋白质体公司 | 针对β-NGF的适体及其在治疗β-NGF介导的疾病和失调中的用途 |
WO2014159669A2 (en) * | 2013-03-14 | 2014-10-02 | Somalogic, Inc. | Aptamers that bind to il-6 and their use in treating or diagnosing il-6 mediated conditions |
Non-Patent Citations (2)
Title |
---|
ALEXEY BEREZHNOY ET AL.: "Isolation and Optimization of Murine IL-10 Receptor Blocking Oligonucleotide Aptamers Using High-throughput Sequencing", 《MOLECULAR THERAPY》 * |
王成龙 等: "特异识别已分化 PC12 细胞的 ssDNA 适体的结构和功能研究", 《军事医学科学院院刊》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017133159A1 (zh) * | 2016-02-03 | 2017-08-10 | 中国农业科学院北京畜牧兽医研究所 | 一种计算机辅助筛选小分子化合物靶标适配体的实现方法 |
CN109868268A (zh) * | 2017-12-04 | 2019-06-11 | 上海药明康德新药开发有限公司 | 一种优化dna编码化合物文库起始片段的方法 |
CN108387561A (zh) * | 2018-01-18 | 2018-08-10 | 大连理工大学 | 一种基于碱基淬灭荧光原理实现低成本高通量核酸适配体最优序列的优化方法 |
CN108387561B (zh) * | 2018-01-18 | 2019-08-09 | 大连理工大学 | 一种基于碱基淬灭荧光原理实现低成本高通量核酸适配体最优序列的优化方法 |
CN108823211A (zh) * | 2018-05-04 | 2018-11-16 | 中国医学科学院基础医学研究所 | 一种cd19核酸适配体及其应用 |
CN111445951A (zh) * | 2020-03-17 | 2020-07-24 | 大连理工大学 | 一种基于免标记杂交探针竞争法实现低成本高通量核酸适配体最优序列的优化方法 |
CN111445951B (zh) * | 2020-03-17 | 2022-09-20 | 大连理工大学 | 一种基于免标记杂交探针竞争法实现低成本高通量核酸适配体最优序列的优化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Komarova et al. | Inside the black box: what makes SELEX better? | |
Zhu et al. | RNA polymerase II activity revealed by GRO-seq and pNET-seq in Arabidopsis | |
Glenn et al. | Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext) | |
KR102640255B1 (ko) | 감소된 증폭 편향을 갖는 고속-대량 단일 세포 서열분석 | |
Glenn et al. | Adapterama I: universal stubs and primers for thousands of dual-indexed Illumina libraries (iTru & iNext) | |
Gamble et al. | Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing | |
Meyer et al. | Fast evolving 18S rRNA sequences from Solenogastres (Mollusca) resist standard PCR amplification and give new insights into mollusk substitution rate heterogeneity | |
CN104561013A (zh) | 基于高通量测序技术进行核酸适配体序列优化的方法 | |
BR112020012696A2 (pt) | diagnóstico multiplex com base em sistema efetor crispr | |
CA3075303A1 (en) | Multi-effector crispr based diagnostic systems | |
Maier et al. | Essential requirements for the detection and degradation of invaders by the Haloferax volcanii CRISPR/Cas system IB | |
Billis et al. | Comparative transcriptomics between Synechococcus PCC 7942 and Synechocystis PCC 6803 provide insights into mechanisms of stress acclimation | |
AU2014409073B2 (en) | Linker element and method of using same to construct sequencing library | |
JP7332733B2 (ja) | 次世代シークエンシングのための高分子量dnaサンプル追跡タグ | |
Komarova et al. | Implementation of high-throughput sequencing (HTS) in aptamer selection technology | |
Kumazawa et al. | Gene rearrangements in gekkonid mitochondrial genomes with shuffling, loss, and reassignment of tRNA genes | |
CN102181943B (zh) | 一种配对双末端文库构建方法及用该文库进行基因组测序的方法 | |
US20220033811A1 (en) | Method and kit for preparing complementary dna | |
US20220033805A1 (en) | High-throughput single-nuclei and single-cell libraries and methods of making and of using | |
ES2947437T3 (es) | Creación de Códigos de barra compartidos en el ADN, en un solo tubo con perlas, para la secuenciación, haplotipado y ensamblaje preciso y rentable | |
CA3128098A1 (en) | Haplotagging - haplotype phasing and single-tube combinatorial barcoding of nucleic acid molecules using bead-immobilized tn5 transposase | |
CN106192019A (zh) | 用于制备测序文库的组合物和方法 | |
CN105349675A (zh) | 基于双酶切的大黄鱼全基因组SNP和InDel分子标记方法 | |
ES2776202T3 (es) | Captura de conformación cromosómica dirigida | |
Jørgensen et al. | An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150429 |