CN113061648A - 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 - Google Patents
一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 Download PDFInfo
- Publication number
- CN113061648A CN113061648A CN202110314939.5A CN202110314939A CN113061648A CN 113061648 A CN113061648 A CN 113061648A CN 202110314939 A CN202110314939 A CN 202110314939A CN 113061648 A CN113061648 A CN 113061648A
- Authority
- CN
- China
- Prior art keywords
- library
- sample
- sequencing
- transposase
- dna
- 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.)
- Granted
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6806—Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B50/00—Methods of creating libraries, e.g. combinatorial synthesis
- C40B50/06—Biochemical methods, e.g. using enzymes or whole viable microorganisms
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Biochemistry (AREA)
- Zoology (AREA)
- Engineering & Computer Science (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Medicinal Chemistry (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于RNA修饰位点检测技术领域,具体涉及一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用,本发明在文库的构建中,通过Tn5转座酶在mRNA/DNA杂合链上添加测序接头,通过该方法构建得到的m6A修饰检测文库可以极大地降低了高通量测序样品的初始投入量,使用本发明方法可以检测现有MeRIP‑seq方法不能应用的微量样品,使高通量测序样品的初始总RNA投入量降低至0.06微克。
Description
技术领域
本发明属于RNA修饰位点检测技术领域,具体涉及一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用。
背景技术
目前已知在RNA上存在超过100种修饰,其中包括:N6-腺苷酸甲基化(m6A)、N1-腺苷酸甲基化(m1A)、胞嘧啶羟基化(m5C)等,而m6A是其中含量最多的修饰,影响RNA 的代谢,肿瘤的发生、发展以及细胞的生长、发育等。
目前已知的m6A检测方法主要有二大类:一类是质谱检测法,可以检测样本m6A修饰的总体水平,但是无法得知m6A修饰具体的修饰位置以及是哪些基因上存在修饰。另一类是高通量测序法,如应用最广泛的Methylated RNA immunoprecipitation andsequencing (MeRIP-seq or m6A-seq)技术,可以精确的知道基因的某些区域上存在m6A修饰。但是,基因的功能体现在mRNA上,而mRNA只占总RNA的1-2%,因此MeRIP-seq需要从大量的总RNA中提取mRNA,另外,MeRIP-seq是在抗体免疫沉淀后获得RNA样品上开始构建测序文库,受限于连接法构建文库的策略,因此MeRIP-seq的初始样品通常需要大于100微克的总RNA作为实验材料,从而极大的限制了该方法在一些难以获得样品的材料以及病人的组织样品上的应用。虽然最新的研究方案通过改进实验流程和策略,将初始样品的用量降低至0.5微克的总RNA,也可以获得较好的m6A高通量数据,但是仍然不能满足微量样品的需要,特别是单细胞的需要。因此,有必要创建一种新的高通量测序方法,以便于研究微量样品的m6A修饰。
发明内容
为了克服上述现有技术的不足,本发明提出了一种采用Tn5转座酶辅助构建微量样品 m6A修饰检测文库的方法,通过Tn5转座酶在mRNA/DNA杂合链上添加测序接头得到的m6A修饰检测文库可以极大地降低了高通量测序样品的初始投入量。
为了实现上述目的,本发明所采用的技术方案是:
本发明提供了一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,即通过 Tn5转座酶在mRNA/DNA杂合链上添加测序接头,构建得到Input测序文库和IP测序文库,所述测序接头的核苷酸序列如SEQ ID NO:1所示。
作为本发明的一个优选的实施方式,上述的一种采用Tn5转座酶辅助构建微量样品m6A 修饰检测文库的方法,包括以下步骤:
S1、提取待测样品的总RNA或含有总RNA的细胞裂解物;
S2、通过Tn5转座酶建库试剂盒去除总RNA中的基因组DNA;
S3、以去除基因组DNA的总RNA作为模板,Tn5转座酶建库试剂盒中的Oligo(dT)20VN引物作为反转录的引物,通过转录获得mRNA/DNA的杂合链产物;
S4、通过Tn5转座酶建库试剂盒使mRNA/DNA的杂合链片段化,并且加上测序接头,所述测序接头的核苷酸序列如SEQ ID NO:1所示;
S5、通过Tn5转座酶建库试剂盒将片段化产物上产生的缺口补平;
S6、将补平缺口的产物分为两份:1/9体积作为Input样品,直接用于构建常规的Tn5测序文库;剩下的8/9体积作为IP样品,用于构建IP测序文库;
S7、首先对IP样品进行变性处理,然后往变性后的IP样品中加入m6A抗体和免疫磁珠,对含有m6A甲基化的mRNA/DNA片段进行富集;
S8、对IP样品和Input样品进行第一轮测序文库的扩增;
S9、使用磁珠纯化第一轮扩增的产物,并洗脱磁珠上的DNA;
S10、对IP样品和Input样品的洗脱产物进行第二轮测序文库的扩增;
S11、使用磁珠纯化第二轮扩增的产物,并洗脱磁珠上的DNA,选取符合文库浓度和文库大小的扩增产物,即可构建得到m6A修饰检测文库,包括Input测序文库和IP测序文库。
优选地,所述总RNA的量不小于60ng或细胞裂解物中的细胞数不少于2000个(可保证细胞裂解物中的总RNA的量不小于60ng)。
优选地,在步骤S8中,可以将富集后的IP样品直接进行第一轮测序文库的扩增,也可以先将富集后的IP样品转变为cDNA,然后再进行第一轮测序文库的扩增。
优选地,步骤S2为使用RNA建库试剂盒中的5×gDNAwiper Mix试剂将总RNA中的基因组DNA去除。
进一步地,具体的DNA去除方法如下:
将总RNA置于70℃下加热3min,迅速置于冰上骤冷,并在冰上静置2min(预变性),然后按照下列组分和成分配制体系:
用移液器轻轻吹打混匀,42℃孵育2min。
优选地,步骤S3中的转录为:首先按照下列组分和成分配制体系:
去除基因组DNA的总RNA | 10μL |
10×RT Mix | 2μL |
HiScript III Enzyme Mix | 2μL |
Oligo(dT)<sub>20</sub>VN | 1μL |
Nuclease-free ddH<sub>2</sub>O | 5μL |
然后混匀反应体系,在PCR仪内,按照以下参数设置反应:42℃60min,10cycles of[50℃2min,42℃2min],70℃10min。
优选地,步骤S4为使用Tn5转座酶建库试剂盒中的Tn5 VR100试剂使mRNA/DNA的杂合链片段化,并且加上测序接头。
进一步地,具体的片段化方法如下:
首先按照下列组分和成分配制体系:
杂合链产物 | 20μL |
Tagment buffer | 10μL |
Tn5 VR100 | 2μL |
然后混匀反应体系,在PCR仪内,按照以下参数设置反应:55℃30min,反应完成后,加入2μL TStop Solution;充分混匀,室温孵育5min。
优选地,步骤S5为使用Tn5转座酶建库试剂盒中的TSE试剂将片段化产物上产生的缺口补平。
进一步地,具体的缺口补平方法如下:
首先按照下列组分和成分配制体系:
片段化产物 | 34μL |
VAHTS HiFi Amplification Mix | 50μL |
TSE | 1μL |
Nuclease-free ddH<sub>2</sub>O | 5μL |
然后混匀反应体系,在PCR仪内,按照以下参数设置反应:55℃15min,60℃15min,85℃5min。
优选地,所述免疫磁珠在使用前需要进行预处理,即用免疫磁珠的预先处理,用反应缓冲液洗2次免疫磁珠(Protein G beads),然后用反应缓冲液重悬免疫磁珠。同时,步骤S7的免疫沉淀反应(IP)后需要依次用反应缓冲液,低盐缓冲液,高盐缓冲液冲洗免疫磁珠2次,降低非特异性结合。
反应缓冲液、低盐缓冲液和高盐缓冲液的组成如下表所示:
反应缓冲液组分 | 最终浓度 |
NaCl | 150mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
低盐缓冲液组分 | 最终浓度 |
NaCl | 50mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
高盐缓冲液组分 | 最终浓度 |
NaCl | 500mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
优选地,步骤S7中,IP样品在富集前需进行变性处理。进一步地,所述变性处理的方法为:置于PCR仪内,95℃10min或者95℃20min。
通过在抗体免疫沉淀反应前,将样品进行95℃10min或者95℃20min的变性处理,从而将mRNA/DNA杂合链部分或者完全打开,更有利于后续对含有m6A甲基化的mRNA/DNA 片段进行富集。
本发明还提供了通过上述的方法构建得到的m6A修饰检测文库在m6A修饰检测中的应用。
优选地,上述应用的领域为非疾病诊断治疗领域。
本发明还提供了一种Tn5转座酶辅助的检测微量样品m6A修饰的高通量测序方法(此方法命名为tMeRIP-seq),具体为:将通过上述的方法构建得到的m6A修饰检测文库采用高通量测序平台进行高通量测序,下载测序结果经数据分析后即可得到m6A修饰检测结果。
优选地,所述高通量测序平台包括但不局限于illumina的HiseqX、MiSeq等测序平台, ABI的Series Genetic Analyzer和Ion Torrent PGM测序平台。
与现有技术相比,本发明的有益效果是:
本发明提出了一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,在文库的构建中,通过Tn5转座酶在mRNA/DNA杂合链上添加测序接头,通过该方法构建得到的m6A修饰检测文库可以极大地降低了高通量测序样品的初始投入量。
总体而言,本发明具有以下优点:
(1)通过使用Oligo(dT)20VN作为反转录引物,从而获得mRNA/DNA杂合链产物,不仅节省了初始样品rRNA去除的步骤,而且可以有效的避免rRNA污染测序数据,并降低成本。
(2)通过Tn5转座酶将mRNA/DNA杂合链片段化,并且加上测序接头,从而可以极大的降低高通量测序样品的初始投入量:降低至0.06微克。
(3)使用本发明方法可以检测现有MeRIP-seq方法不能应用的微量样品:例如微量病理样品。
附图说明
图1为片段化流程图;
图2为用于m6A的IP文库构建流程图;
图3为测序文库大小分布图;
图4为实施例1的m6A修饰检测的数据分析图;
图5为实施例2的m6A修饰检测的数据分析图;
图6为血液样品的m6A修饰检测流程图;
图7为实施例2的m6A修饰检测的数据分析图。
具体实施方式
下面对本发明的具体实施方式作进一步说明。在此需要说明的是,对于这些实施方式的说明用于帮助理解本发明,但并不构成对本发明的限定。此外,下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
下述实施例中的实验方法,如无特殊说明,均为常规方法,下述实施例中所用的试验材料,如无特殊说明,均为可通过常规的商业途径购买得到的。
本实施例中所采用的免疫磁珠在使用前需要进行预处理,即用免疫磁珠的预先处理,用反应缓冲液洗2次免疫磁珠(Protein G beads),然后用250μL反应缓冲液重悬免疫磁珠。
同时,免疫沉淀反应(IP)后需要依次用反应缓冲液,低盐缓冲液,高盐缓冲液冲洗免疫磁珠2次,降低非特异性结合。
反应缓冲液、低盐缓冲液和高盐缓冲液的组成如下表所示:
反应缓冲液组分 | 最终浓度 |
NaCl | 150mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
低盐缓冲液组分 | 最终浓度 |
NaCl | 50mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
高盐缓冲液组分 | 最终浓度 |
NaCl | 500mM |
Tris-HCl(PH7.4) | 10mM |
Igepal CA-630 | 0.01% |
例1一种Tn5转座酶辅助的检测微量样品m6A修饰的高通量测序方法
本实施例需要用到的引物(Vazyme,cat.no.TD202)如下(N5-index包括N501-508、N7-index包括N701-712):
本实施例需要的总RNA为0.06微克,步骤(1)-(7)的流程如图1所示,步骤(8)-(16)的流程如图2所示:
(1)使用二氧化碳培养箱,将293T细胞(购买于中国科学院细胞库)置于改良DMEM培养基(DMEM;273 Corning,cat.no.10-017-CV),并添加10%(v/v)胎牛血清(FBS;Gibco),5%CO2,37℃培养至细胞数达到十万级。
(2)将细胞培养基弃掉,用PBS清洗两次,使用胰酶消化细胞,对细胞悬浊液进行离心(1000rpm/5min),弃上清,保留底层细胞。
(3)使用经典的Trizol(Thermo Fisher Scientific,catalog no.15596026)法提取总RNA,按照产品说明书进行操作,获得总RNA(取0.06微克作为检测样品)。
(4)使用Tn5转座酶建库试剂盒(Vazyme,TR501-01)中的5×gDNA wiper Mix试剂,将总RNA中的基因组DNA去除,具体如下:
将60ng总RNA置于70℃下加热3min,迅速置于冰上骤冷,并在冰上静置2min(预变性),然后按照下列组分和成分配制体系:
预变性RNA | 60ng |
5×gDNA wiper Mix(Vazyme,TR501-01) | 2μL |
Nuclease-free ddH<sub>2</sub>O | up to 10μL |
用移液器轻轻吹打混匀反应体系,42℃孵育2min。
(5)使用Tn5转座酶建库试剂盒(Vazyme,TR501-01)中的试剂合成mRNA/DNA杂合链,具体如下:
首先按照下列组分和成分配制体系:
然后用移液器轻轻吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:42℃60min, 10cycles of[50℃2min,42℃2min],70℃10min。
(6)使用Tn5转座酶建库试剂盒(Vazyme,TR501-01)中的Tn5 VR100试剂将mRNA/DNA 杂合链片段化,并且在片段的两端分别随机加上测序接头(接头存在于Tn5 VR100酶上,Tn5 VR100酶在片段化模板的时候,会将两端的接头序列随机插入模板的两端): 5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3′(接头1)和 5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3′(接头2),如SEQ ID NO:1所示。
具体操作如下:
首先按照下列组分和成分配制体系:
上一步的反应液(杂合链产物) | 20μL |
Tagment buffer | 10μL |
Tn5 VR100(Vazyme,TR501-01) | 2μL |
然后用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:55℃30min,反应完成后,加入2μL TStop Solution;充分混匀,室温孵育5min。
(7)使用Tn5转座酶建库试剂盒(Vazyme,TR501-01)中的TSE试剂将片段化后产生的缺口补平,具体如下:
首先按照下列组分和成分配制体系:
上一步的产物(片段化产物) | 34μL |
VAHTS HiFi Amplification Mix | 50μL |
TSE(Vazyme,TR501-01) | 1μL |
Nuclease-free ddH<sub>2</sub>O | 5μL |
然后用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:55℃15min,60℃15min,85℃5min。
(8)将上一步补平缺口的产物分成两份:1/9体积(10微升)作为Input样品,直接用于构建常规的Tn5测序文库;剩下的8/9体积(80微升)作为IP样品,用于构建IP测序文库;
(9)首先对IP样品进行变性,即在PCR仪内,按照以下参数设置反应:95℃10min(方法一)或者95℃20min(方法二),然后迅速置于冰上骤冷,并在冰上静置2min;然后在变性后的IP样品中加入m6A抗体和免疫磁珠(需要预先处理),对含有m6A甲基化的 mRNA/DNA片段进行富集,具体如下:
首先按照下列组分和成分配制体系:
等分后的变性IP样品(加入了m6A抗体和免疫磁珠) | 40μL |
m6A antibody(Cell Signaling Technology,56593) | 2μL |
Protein G beads(Thermo Fisher Scientific,1004D) | 250μL |
Recombinant ribonuclease inhibitor(Takara,2313) | 4μL |
然后用移液器充分吹打混匀反应体系,在4℃冰箱内,将样品置于旋转仪上颠倒混匀过夜。
(10)去除免疫磁珠的非特异性结合,即依次用反应缓冲液,低盐缓冲液,高盐缓冲液冲洗免疫磁珠2次。
(11)第一轮扩增测序文库:
1)IP样品第一轮扩增测序文库:方法一,将免疫磁珠中的液体吸干净,随后加入PCR 反应体系开始第一轮扩增测序文库(参考文献:McEvoy,C.R.,Semple,T.,Yellapu,B.,Choong, D.Y.,Xu,H.,Mir Arnau,G.,Fellowes,A.P.,and Fox,S.B.(2020).Improvednext-generation sequencing pre-capture library yields and sequencingparameters using on-bead PCR.Biotechniques 68,48-51),以免疫磁珠为样品,加入以下反应体系:
2xKAPA HiFi Hot Start Ready Mix(KAPA Biosystems,KK2601) | 25μL |
Indexed common primers N5引物序列 | 1.5μL |
Indexed common primers N7引物序列 | 1.5μL |
P5-primer | 1.5μL |
P7-primer | 1.5μL |
Nuclease-free ddH<sub>2</sub>O | 19μL |
用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:98℃3min,10 cycles of[98℃20s,60℃20s,72℃1min],72℃3min。
方法二,将免疫磁珠中的液体吸干净,随后加入反转录体系,将RNA反转录成cDNA,以免疫磁珠为样品,加入以下反应体系:
步骤一 | |
2xKAPA HiFi Hot Start Ready Mix(KAPA Biosystems,KK2601) | 10μL |
Indexed common primers(Vazyme,cat.no.TD202) | 0.5μL |
Nuclease-free ddH<sub>2</sub>O | 9.5μL |
步骤二 | |
RiboLock RNase Inhibitor(Thermo Scientific,K1622) | 0.5μL |
RevertAid M-MuLV RT(Thermo Scientific,K1622) | 1μL |
首先按照步骤一加入相应的试剂,用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:85℃30s,60℃2min,缓慢冷却到4℃;然后加入步骤二的试剂,用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:42℃30min,反应结束后,加入PCR反应体系开始第一轮扩增测序文库,具体按照以下反应体系进行:
2xKAPA HiFi Hot Start Ready Mix(KAPA Biosystems,KK2601) | 25μL |
Indexed common primers N5引物序列 | 1.0μL |
Indexed common primers N7引物序列 | 1.0μL |
P5-primer | 1.0μL |
P7-primer | 1.0μL |
上一步反应产物(步骤二的PCR产物) | 21μL |
用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:98℃3min,10 cycles of[98℃20s,60℃20s,72℃1min],72℃3min。
2)Input样品第一轮扩增测序文库:将Input样品加入PCR反应体系开始第一轮扩增测序文库,反应体系具体如下:
(12)按1:1的体积比使用商业化的PCR产物纯化磁珠(Vazyme,N411)纯化PCR产物(IP样品和Input样品的PCR产物),最后用10μL无核酶的纯水洗脱磁珠上的DNA。
(13)第二轮扩增测序文库,具体按照以下反应体系进行:
上一步的洗脱产物(IP样品和Input样品的DNA) | 10μL |
KAPA HiFi Hot Start Ready Mix(KAPA Biosystems,KK2601) | 25μL |
P5-primer | 1.5μL |
P7-primer | 1.5μL |
Nuclease-free ddH<sub>2</sub>O | 12μL |
用移液器充分吹打混匀反应体系,在PCR仪内,按照以下参数设置反应:
Input样品:98℃3min,8 cycles of[98℃20s,60℃20s,72℃1min],72℃3min。
IP样品:98℃3min,12 cycles of[98℃20s,60℃20s,72℃1min],72℃3min。
(14)按1:1的体积比使用商业化的PCR产物纯化磁珠(Vazyme,N411)纯化PCR产物(IP样品和Input样品的PCR产物),最后用15μL洗脱磁珠上的DNA。一方面用商业化试剂Qubit dsDNA HS Assay kit(Thermo Fisher Scientific,Q32851)测定文库浓度(文库浓度大于 4ng/μL,符合测序要求[即不低于2ng/μL]),另一方面检测文库大小的分布情况(如图3所示,文库大小为200-500bp,符合要求)。
由此构建得到Input测序文库和IP测序文库。
(15)将Input测序文库样品和IP测序文库样品送测序公司,在HiSeqX测序平台进行高通量测序。
(16)下载测序结果,并进行分析。
首先,使用FastQC(http://www.bioinformatics.babraham.ac.uk/projects/fastqc/)软件对测序结果进行质量评估,并使用Cutadapt(https://cutadapt.readthedocs.io/en/stable/)进行质量控制,去除接头序列、低质量的碱基。然后,使用HISAT2(https://daehwankimlab.github.io/hisat2/) 将测序片段比对到参考基因组,使用Samtools(https://github.com/samtools/samtools)软件对比对结果进行排序,建索引等处理。最后,使用R语言包exomePeak2 (https://github.com/ZW-xjtlu/exomePeak2)进行call peak分析,以鉴定出m6A修饰的富集区域。
目前,大量的研究都表明m6A修饰趋向于富集在终止密码子附近,并且许多研究都以 m6A修饰在mRNA上的分布模式作为判断其结果好坏的依据。根据方法一和方法二鉴定得到的m6APeak,展示了其在mRNA的分布情况,结果如图4-A所示,样品的m6Apeak主要分布在终止编码区域附近,Motif分析结果所示如图4-B和4-C,Motif具有明显的“RRACH”(R =Gor A;H=A,C,or U)富集,与已知的结论相符合(参考文献:Meyer,K.D.,Saletore,Y.,Zumbo,P.,Elemento,O.,Mason,C.E.,and Jaffrey,S.R.(2012).Comprehensive analysisof mRNA methylation reveals enrichment in 3'UTRs and near stop codons.Cell149,1635-1646)。
实施例2一种Tn5转座酶辅助的检测微量样品m6A修饰的高通量测序方法
本实施例的检测样品为293T细胞(2000个细胞),即以细胞裂解物作为检测样品。
本实施例需要用到的引物同实施例1。
步骤(1)-(7)的流程如图1所示,步骤(8)-(16)的流程如图2所示:
(1)使用二氧化碳培养箱,将293T细胞置于改良DMEM培养基(DMEM;273Corning,cat.no.10-017-CV)并添加10%(v/v)胎牛血清(FBS;Gibco)。
(2)将细胞培养基弃掉,用PBS清洗两次,使用胰酶消化细胞,对细胞悬浊液进行离心(1000rpm/5min),弃上清,用PBS重选细胞,随后用流式细胞分选仪,分选出2000个细胞。
(3)将3μL细胞裂解液加入步骤2中分选好的细胞中并且充分混匀,然后将样品置于 PCR仪器中,72℃3min,随后立即插入冰中,以备使用(裂解细胞膜,获得含有总RNA的细胞裂解物)。
细胞裂解液的组成如下:
细胞裂解液 | 50μL |
TritonX-100(Sigma-Aldrich,cat.no.T9284) | 0.1μL |
Recombinant ribonuclease inhibitor(Takara,cat.no.2313) | 2μL |
Nuclease-free ddH<sub>2</sub>O | 47.8μL |
(4)以下的m6A修饰检测方法与同实施例1。
根据方法一和方法二鉴定得到的m6A Peak,展示了其在mRNA的分布情况,结果如图5 所示,样品的m6A peak主要分布在终止编码区域附近,Motif分析结果所示如图5-A和5-B, Motif具有明显的“RRACH”(R=G orA;H=A,C,or U)富集,与已知的结论相符合(参考文献:Meyer,K.D.,Saletore,Y.,Zumbo,P.,Elemento,O.,Mason,C.E.,and Jaffrey,S.R.(2012). Comprehensive analysis of mRNA methylation reveals enrichment in 3'UTRs and near stop codons. Cell 149,1635-1646)。
实施例3一种Tn5转座酶辅助的检测微量样品m6A修饰的高通量测序方法
本实施例的检测样品为一滴血样品(5微升)。
本实施例需要用到的引物同实施例1。
本实施例的操作流程如图6所示:
(1)使用采血针从健康的志愿者获得指尖外周血5微升/样品。
(2)使用红细胞裂解液(Roche,11814389001),将外周血中的红细胞去除,从而获得大量的白细胞。
(3)使用经典的Trizol(Thermo Fisher Scientific,catalog no.15596026)法提取总RNA,按照说明书进行操作,获得总RNA。
(4)以下的m6A修饰检测方法与同实施例1,其中步骤(11)采用方法二。
根据方法二鉴定得到的m6APeak,展示了其在mRNA的分布情况,结果如图7-A所示,样品的m6Apeak主要分布在终止编码区域附近,与已知的结论相符合。Motif分析结果所示如图7-B,Motif具有明显的“RRACH”(R=G orA;H=A,C,or U)富集,与已知的结论相符合(参考文献:Meyer,K.D.,Saletore,Y.,Zumbo,P.,Elemento,O.,Mason,C.E.,and Jaffrey,S.R. (2012).Comprehensive analysis of mRNA methylation reveals enrichment in3'UTRs and near stop codons.Cell 149,1635-1646)。
以上对本发明的实施方式作了详细说明,但本发明不限于所描述的实施方式。对于本领域的技术人员而言,在不脱离本发明原理和精神的情况下,对这些实施方式进行多种变化、修改、替换和变型,仍落入本发明的保护范围内。
序列表
<110> 中山大学
<120> 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用
<160> 27
<170> SIPOSequenceListing 1.0
<210> 1
<211> 33
<212> DNA
<213> 接头1(人工序列)
<400> 1
tcgtcggcag cgtcagatgt gtataagaga cag 33
<210> 2
<211> 34
<212> DNA
<213> 接头2(人工序列)
<400> 2
gtctcgtggg ctcggagatg tgtataagag acag 34
<210> 3
<211> 8
<212> DNA
<213> N501(人工序列)
<400> 3
tagatcgc 8
<210> 4
<211> 8
<212> DNA
<213> N502(人工序列)
<400> 4
ctctctat 8
<210> 5
<211> 8
<212> DNA
<213> N503(人工序列)
<400> 5
tatcctct 8
<210> 6
<211> 8
<212> DNA
<213> N504(人工序列)
<400> 6
agagtaga 8
<210> 7
<211> 8
<212> DNA
<213> N505(人工序列)
<400> 7
gtaaggag 8
<210> 8
<211> 8
<212> DNA
<213> N506(人工序列)
<400> 8
actgcata 8
<210> 9
<211> 8
<212> DNA
<213> N507(人工序列)
<400> 9
aaggagta 8
<210> 10
<211> 8
<212> DNA
<213> N508(人工序列)
<400> 10
ctaagcct 8
<210> 11
<211> 8
<212> DNA
<213> N701(人工序列)
<400> 11
taaggcga 8
<210> 12
<211> 8
<212> DNA
<213> N702(人工序列)
<400> 12
cgtactag 8
<210> 13
<211> 8
<212> DNA
<213> N703(人工序列)
<400> 13
aggcagaa 8
<210> 14
<211> 8
<212> DNA
<213> N704(人工序列)
<400> 14
tcctgagc 8
<210> 15
<211> 8
<212> DNA
<213> N705(人工序列)
<400> 15
ggactcct 8
<210> 16
<211> 8
<212> DNA
<213> N706(人工序列)
<400> 16
taggcatg 8
<210> 17
<211> 8
<212> DNA
<213> N707(人工序列)
<400> 17
ctctctac 8
<210> 18
<211> 8
<212> DNA
<213> N708(人工序列)
<400> 18
cagagagg 8
<210> 19
<211> 8
<212> DNA
<213> N709(人工序列)
<400> 19
gctacgct 8
<210> 20
<211> 8
<212> DNA
<213> N710(人工序列)
<400> 20
cgaggctg 8
<210> 21
<211> 8
<212> DNA
<213> N711(人工序列)
<400> 21
aagaggca 8
<210> 22
<211> 8
<212> DNA
<213> N712(人工序列)
<400> 22
gtagagga 8
<210> 23
<211> 20
<212> DNA
<213> P5-Primer(人工序列)
<400> 23
aatgatacgg cgaccaccga 20
<210> 24
<211> 21
<212> DNA
<213> P7-Primer(人工序列)
<400> 24
caagcagaag acggcatacg a 21
<210> 25
<211> 22
<212> DNA
<213> Oligo (dT)20VN(人工序列)
<220>
<223> V=dA, dC, or dG;N= dA, dC, dG or dT
<220>
<221> misc_feature
<222> (22)..(22)
<223> n is a, c, g, t or u
<400> 25
tttttttttt tttttttttt vn 22
<210> 26
<211> 43
<212> DNA
<213> N5引物(人工序列)
<220>
<223> 5′-aatgatacggcgaccaccgagatctacac-[N5-index]-tcgtcggcagcgtc-3′
<400> 26
aatgatacgg cgaccaccga gatctacact cgtcggcagc gtc 43
<210> 27
<211> 39
<212> DNA
<213> N7引物(人工序列)
<220>
<223> 5′-caagcagaagacggcatacgagat-[N7-index]-gtctcgtgggctcgg-3′
<400> 27
caagcagaag acggcatacg agatgtctcg tgggctcgg 39
Claims (10)
1.一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,通过Tn5转座酶在mRNA/DNA杂合链上添加测序接头,构建得到Input测序文库和IP测序文库,所述测序接头的核苷酸序列如SEQ ID NO:1所示。
2.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,包括以下步骤:
S1、提取待测样品的总RNA或含有总RNA的细胞裂解物;
S2、通过Tn5转座酶建库试剂盒去除总RNA中的基因组DNA;
S3、以去除基因组DNA的总RNA作为模板,Tn5转座酶建库试剂盒中的Oligo(dT)20VN引物作为反转录的引物,通过转录获得mRNA/DNA的杂合链产物;
S4、通过Tn5转座酶建库试剂盒使mRNA/DNA的杂合链片段化,并且加上测序接头,所述测序接头的核苷酸序列如SEQ ID NO:1所示;
S5、通过Tn5转座酶建库试剂盒将片段化产物上产生的缺口补平;
S6、将补平缺口的产物分为两份:1/9体积作为Input样品,直接用于构建常规的Tn5测序文库;剩下的8/9体积作为IP样品,用于构建IP测序文库;
S7、首先对IP样品进行变性处理,然后往变性后的IP样品中加入m6A抗体和免疫磁珠,对含有m6A甲基化的mRNA/DNA片段进行富集;
S8、对IP样品和Input样品进行第一轮测序文库的扩增;
S9、使用磁珠纯化第一轮扩增的产物,并洗脱磁珠上的DNA;
S10、对IP样品和Input样品的洗脱产物进行第二轮测序文库的扩增;
S11、使用磁珠纯化第二轮扩增的产物,并洗脱磁珠上的DNA,选取符合文库浓度和文库大小的扩增产物,即可构建得到m6A修饰检测文库,包括Input测序文库和IP测序文库。
3.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,所述总RNA的量不小于60ng或细胞裂解物中的细胞数不少于2000个。
4.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,在步骤S8中,可以将富集后的IP样品直接进行第一轮测序文库的扩增,也可以先将富集后的IP样品转变为cDNA,然后再进行第一轮测序文库的扩增。
5.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,步骤S2为使用RNA建库试剂盒中的5×gDNAwiper Mix试剂将总RNA中的基因组DNA去除。
6.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,步骤S4为使用Tn5转座酶建库试剂盒中的Tn5 VR100试剂使mRNA/DNA的杂合链片段化,并且加上测序接头。
7.根据权利要求1所述的一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法,其特征在于,步骤S5为使用Tn5转座酶建库试剂盒中的TSE试剂将片段化产物上产生的缺口补平。
8.通过权利要求1-7任一项所述的方法构建得到的m6A修饰检测文库在m6A修饰检测中的应用。
9.根据权利要求8所述的应用,其特征在于,所述应用的领域为非疾病诊断治疗领域。
10.一种Tn5转座酶辅助的检测微量样品m6A修饰的高通量测序方法,其特征在于,将通过权利要求1-7任一项所述的方法构建得到的m6A修饰检测文库采用高通量测序平台进行高通量测序,下载测序结果经数据分析后即可得到m6A修饰检测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110314939.5A CN113061648B (zh) | 2021-03-24 | 2021-03-24 | 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110314939.5A CN113061648B (zh) | 2021-03-24 | 2021-03-24 | 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113061648A true CN113061648A (zh) | 2021-07-02 |
CN113061648B CN113061648B (zh) | 2022-04-19 |
Family
ID=76561638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110314939.5A Active CN113061648B (zh) | 2021-03-24 | 2021-03-24 | 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113061648B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113584135A (zh) * | 2021-07-26 | 2021-11-02 | 中山大学 | 一种混样检测rna修饰并实现精准定量的方法 |
CN113718343A (zh) * | 2021-09-03 | 2021-11-30 | 翌圣生物科技(上海)股份有限公司 | 极速rna建库方法及试剂盒 |
CN114645074A (zh) * | 2022-03-07 | 2022-06-21 | 华南农业大学 | 一种单花蕾微量rna建库方法 |
CN114758723A (zh) * | 2022-03-31 | 2022-07-15 | 广州华银医学检验中心有限公司 | 一种基于MeRIP测序技术检测肿瘤治疗靶点的方法和系统 |
CN115386966A (zh) * | 2022-10-26 | 2022-11-25 | 北京寻因生物科技有限公司 | Dna表观修饰的建库方法、测序方法及其建库试剂盒 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105018617A (zh) * | 2015-07-20 | 2015-11-04 | 浙江大学 | 单个基因mRNA甲基化水平检测方法 |
CN106047997A (zh) * | 2016-05-27 | 2016-10-26 | 浙江大学 | mRNA甲基化高通量检测方法 |
CN108588200A (zh) * | 2018-05-06 | 2018-09-28 | 湖南大地同年生物科技有限公司 | 一种R-Loop高通量测序文库构建方法 |
CN110904192A (zh) * | 2018-12-28 | 2020-03-24 | 广州表观生物科技有限公司 | 一种超微量RNA甲基化m6A检测方法及其应用 |
CN111154837A (zh) * | 2019-09-02 | 2020-05-15 | 浙江大学 | 一种全转录组范围单碱基分辨率检测rna n6-甲基腺嘌呤修饰的方法及试剂盒 |
CN111434777A (zh) * | 2019-01-11 | 2020-07-21 | 中国人民解放军军事科学院军事医学研究院 | 一种m6A修饰RNA检测方法及用途 |
CN111560423A (zh) * | 2020-06-05 | 2020-08-21 | 中山大学孙逸仙纪念医院 | 一种高通量高灵敏度单碱基分辨率检测RNA m6A的方法及其应用 |
WO2020169826A1 (en) * | 2019-02-21 | 2020-08-27 | Universität Heidelberg | Biomarker panel for diagnosis and prognosis of cancer |
-
2021
- 2021-03-24 CN CN202110314939.5A patent/CN113061648B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105018617A (zh) * | 2015-07-20 | 2015-11-04 | 浙江大学 | 单个基因mRNA甲基化水平检测方法 |
CN106047997A (zh) * | 2016-05-27 | 2016-10-26 | 浙江大学 | mRNA甲基化高通量检测方法 |
CN108588200A (zh) * | 2018-05-06 | 2018-09-28 | 湖南大地同年生物科技有限公司 | 一种R-Loop高通量测序文库构建方法 |
CN110904192A (zh) * | 2018-12-28 | 2020-03-24 | 广州表观生物科技有限公司 | 一种超微量RNA甲基化m6A检测方法及其应用 |
CN111434777A (zh) * | 2019-01-11 | 2020-07-21 | 中国人民解放军军事科学院军事医学研究院 | 一种m6A修饰RNA检测方法及用途 |
WO2020169826A1 (en) * | 2019-02-21 | 2020-08-27 | Universität Heidelberg | Biomarker panel for diagnosis and prognosis of cancer |
CN111154837A (zh) * | 2019-09-02 | 2020-05-15 | 浙江大学 | 一种全转录组范围单碱基分辨率检测rna n6-甲基腺嘌呤修饰的方法及试剂盒 |
CN111560423A (zh) * | 2020-06-05 | 2020-08-21 | 中山大学孙逸仙纪念医院 | 一种高通量高灵敏度单碱基分辨率检测RNA m6A的方法及其应用 |
Non-Patent Citations (3)
Title |
---|
BO LU等: "Transposase-assisted tagmentation of RNA/DNA hybrid duplexes", 《ELIFE.》 * |
LIN DI等: "RNA sequencing by direct tagmentation of RNA/DNA hybrids", 《PNAS》 * |
YONG ZENG等: "Refined RIP-seq protocol for epitranscriptome analysis with low input materials", 《PLOS BIOL.》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113584135A (zh) * | 2021-07-26 | 2021-11-02 | 中山大学 | 一种混样检测rna修饰并实现精准定量的方法 |
CN113584135B (zh) * | 2021-07-26 | 2022-05-27 | 中山大学 | 一种混样检测rna修饰并实现精准定量的方法 |
CN113718343A (zh) * | 2021-09-03 | 2021-11-30 | 翌圣生物科技(上海)股份有限公司 | 极速rna建库方法及试剂盒 |
CN114645074A (zh) * | 2022-03-07 | 2022-06-21 | 华南农业大学 | 一种单花蕾微量rna建库方法 |
CN114758723A (zh) * | 2022-03-31 | 2022-07-15 | 广州华银医学检验中心有限公司 | 一种基于MeRIP测序技术检测肿瘤治疗靶点的方法和系统 |
CN115386966A (zh) * | 2022-10-26 | 2022-11-25 | 北京寻因生物科技有限公司 | Dna表观修饰的建库方法、测序方法及其建库试剂盒 |
CN115386966B (zh) * | 2022-10-26 | 2023-03-21 | 北京寻因生物科技有限公司 | Dna表观修饰的建库方法、测序方法及其建库试剂盒 |
Also Published As
Publication number | Publication date |
---|---|
CN113061648B (zh) | 2022-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113061648B (zh) | 一种采用Tn5转座酶辅助构建微量样品m6A修饰检测文库的方法及其应用 | |
CN109996892B (zh) | 单细胞测序文库的构建方法及其应用 | |
JP6324962B2 (ja) | 標的rna枯渇化組成物を調製するための方法およびキット | |
CN111139532B (zh) | 一种dna和rna同时构建测序文库的方法及试剂盒 | |
EP2725125B1 (en) | High throughput methylation detection method | |
WO2017084023A1 (zh) | 一种高通量的单细胞转录组建库方法 | |
CN110886021B (zh) | 一种单细胞dna文库的构建方法 | |
CN111808854A (zh) | 带有分子条码的平衡接头及快速构建转录组文库的方法 | |
CN113668068A (zh) | 基因组甲基化文库及其制备方法和应用 | |
CN113308514A (zh) | 微量m6A的检测文库构建方法和试剂盒、高通量检测方法 | |
CN114774527A (zh) | 一种高通量单细胞转录组测序方法及其应用 | |
CN114507711A (zh) | 一种单细胞转录组测序方法及其应用 | |
CN113638055B (zh) | 一种制备双链rna测序文库的方法 | |
CN110951827B (zh) | 一种转录组测序文库快速构建方法及其应用 | |
CN115715323A (zh) | 一种高兼容性的PCR-free建库和测序方法 | |
US20240182962A1 (en) | Ultra-high-throughput single cell sequencing method | |
CN113584135B (zh) | 一种混样检测rna修饰并实现精准定量的方法 | |
CN114875118A (zh) | 确定细胞谱系的方法、试剂盒和装置 | |
CN106222164B (zh) | 利用转座子酶进行体外核酸单向扩增的方法、组合物和试剂盒 | |
CN111705119A (zh) | 一种检测毛竹环状rna表观遗传修饰的方法 | |
CN112011625B (zh) | 一种用于评估猪组蛋白修饰的富集结果的检测方法 | |
WO2023159416A1 (zh) | 一种单细胞转录组测序方法及其应用 | |
Tian | mRNA 3’End Processing and Metabolism | |
CN117625595A (zh) | 染色质可及性/转录组共检测试剂盒及其应用 | |
CN112410406A (zh) | 一种确定文库扩增循环数的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |