CN116790576A - 一种消解微生物组文库中rRNA的方法 - Google Patents
一种消解微生物组文库中rRNA的方法 Download PDFInfo
- Publication number
- CN116790576A CN116790576A CN202310383401.9A CN202310383401A CN116790576A CN 116790576 A CN116790576 A CN 116790576A CN 202310383401 A CN202310383401 A CN 202310383401A CN 116790576 A CN116790576 A CN 116790576A
- Authority
- CN
- China
- Prior art keywords
- sequences
- rrna
- sequence
- sgrna
- characteristic
- 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
- 244000005700 microbiome Species 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 15
- 108091027544 Subgenomic mRNA Proteins 0.000 claims abstract description 18
- 230000029087 digestion Effects 0.000 claims abstract description 18
- 239000011324 bead Substances 0.000 claims abstract description 15
- 238000012163 sequencing technique Methods 0.000 claims abstract description 13
- 238000013461 design Methods 0.000 claims abstract description 8
- 238000000246 agarose gel electrophoresis Methods 0.000 claims abstract description 7
- 238000012216 screening Methods 0.000 claims abstract description 7
- 238000002360 preparation method Methods 0.000 claims abstract description 5
- 238000007865 diluting Methods 0.000 claims abstract description 4
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 108091033409 CRISPR Proteins 0.000 claims description 21
- 239000002299 complementary DNA Substances 0.000 claims description 21
- 239000000047 product Substances 0.000 claims description 21
- 108091005804 Peptidases Proteins 0.000 claims description 20
- 239000004365 Protease Substances 0.000 claims description 17
- 102100037486 Reverse transcriptase/ribonuclease H Human genes 0.000 claims description 17
- 150000007523 nucleic acids Chemical group 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 9
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 6
- 238000012408 PCR amplification Methods 0.000 claims description 6
- 238000011529 RT qPCR Methods 0.000 claims description 6
- 239000012264 purified product Substances 0.000 claims description 6
- 108090000623 proteins and genes Proteins 0.000 claims description 5
- 238000000746 purification Methods 0.000 claims description 5
- 101710163270 Nuclease Proteins 0.000 claims description 4
- 238000001976 enzyme digestion Methods 0.000 claims description 4
- 230000000813 microbial effect Effects 0.000 claims description 4
- 238000010839 reverse transcription Methods 0.000 claims description 4
- 241000894007 species Species 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 238000010354 CRISPR gene editing Methods 0.000 claims description 3
- 108020004414 DNA Proteins 0.000 claims description 3
- 102100031780 Endonuclease Human genes 0.000 claims description 3
- 108010042407 Endonucleases Proteins 0.000 claims description 3
- 102000035195 Peptidases Human genes 0.000 claims description 3
- 241000193996 Streptococcus pyogenes Species 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000007853 buffer solution Substances 0.000 claims description 3
- 101150038500 cas9 gene Proteins 0.000 claims description 3
- 238000010790 dilution Methods 0.000 claims description 3
- 239000012895 dilution Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000002184 metal Substances 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000007481 next generation sequencing Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 108091032973 (ribonucleotides)n+m Proteins 0.000 claims 2
- 238000003776 cleavage reaction Methods 0.000 claims 1
- 230000007017 scission Effects 0.000 claims 1
- 241000894006 Bacteria Species 0.000 abstract description 10
- 230000001580 bacterial effect Effects 0.000 abstract description 8
- 241000736262 Microbiota Species 0.000 abstract description 4
- 238000003559 RNA-seq method Methods 0.000 abstract description 4
- 108020005004 Guide RNA Proteins 0.000 description 21
- 238000013459 approach Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010356 CRISPR-Cas9 genome editing Methods 0.000 description 2
- 108091034117 Oligonucleotide Proteins 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010804 cDNA synthesis Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 108020004999 messenger RNA Proteins 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 241000203069 Archaea Species 0.000 description 1
- 208000002064 Dental Plaque Diseases 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- JLCPHMBAVCMARE-UHFFFAOYSA-N [3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-[[3-[[3-[[3-[[3-[[3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-[[5-(2-amino-6-oxo-1H-purin-9-yl)-3-hydroxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxyoxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(5-methyl-2,4-dioxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(6-aminopurin-9-yl)oxolan-2-yl]methoxy-hydroxyphosphoryl]oxy-5-(4-amino-2-oxopyrimidin-1-yl)oxolan-2-yl]methyl [5-(6-aminopurin-9-yl)-2-(hydroxymethyl)oxolan-3-yl] hydrogen phosphate Polymers Cc1cn(C2CC(OP(O)(=O)OCC3OC(CC3OP(O)(=O)OCC3OC(CC3O)n3cnc4c3nc(N)[nH]c4=O)n3cnc4c3nc(N)[nH]c4=O)C(COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3COP(O)(=O)OC3CC(OC3CO)n3cnc4c(N)ncnc34)n3ccc(N)nc3=O)n3cnc4c(N)ncnc34)n3ccc(N)nc3=O)n3ccc(N)nc3=O)n3ccc(N)nc3=O)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cc(C)c(=O)[nH]c3=O)n3cc(C)c(=O)[nH]c3=O)n3ccc(N)nc3=O)n3cc(C)c(=O)[nH]c3=O)n3cnc4c3nc(N)[nH]c4=O)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)n3cnc4c(N)ncnc34)O2)c(=O)[nH]c1=O JLCPHMBAVCMARE-UHFFFAOYSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 238000012350 deep sequencing Methods 0.000 description 1
- 238000001962 electrophoresis Methods 0.000 description 1
- 210000003527 eukaryotic cell Anatomy 0.000 description 1
- 238000001502 gel electrophoresis Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011222 transcriptome analysis Methods 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1003—Extracting or separating nucleic acids from biological samples, e.g. pure separation or isolation methods; Conditions, buffers or apparatuses therefor
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/10—Processes for the isolation, preparation or purification of DNA or RNA
- C12N15/1003—Extracting or separating nucleic acids from biological samples, e.g. pure separation or isolation methods; Conditions, buffers or apparatuses therefor
- C12N15/1006—Extracting or separating nucleic acids from biological samples, e.g. pure separation or isolation methods; Conditions, buffers or apparatuses therefor by means of a solid support carrier, e.g. particles, polymers
- C12N15/1013—Extracting or separating nucleic acids from biological samples, e.g. pure separation or isolation methods; Conditions, buffers or apparatuses therefor by means of a solid support carrier, e.g. particles, polymers by using magnetic beads
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Biotechnology (AREA)
- General Engineering & Computer Science (AREA)
- Crystallography & Structural Chemistry (AREA)
- Plant Pathology (AREA)
- Molecular Biology (AREA)
- Microbiology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于分子生物学技术领域,具体公开了一种消解微生物组文库中rRNA的方法,包括以下步骤:sgRNA设计:汇总58种细菌rRNA序列形成细菌rRNA基因组,将多个含有微生物群的复杂样本的转录组测序序列映射到上述细菌rRNA基因组,挑选覆盖更多序列的特征序列,经过重重筛选最后得到300条特征序列;sgRNA制备:将由生物公司合成的sgDNA通过PCR扩增加上tracrDNA序列,之后通过琼脂糖凝胶电泳验证PCR产物条带大小在130bp,得到的产物经磁珠纯化并转录得到sgRNA,将得到的sgRNA测定浓度并稀释后待用;本发明参考常见的58种细菌基因组序列,筛选设计出300条sgRNA,该300条sgRNA可靶向常见的几乎所有细菌的rRNA,适用于含有多种细菌的复杂转录组文库中rRNA的消解,增加微生物组各种水平RNA测序的覆盖度。
Description
技术领域
本发明涉及分子生物学技术领域,具体为一种消解微生物组文库中rRNA的方法。
背景技术
微生物组(microbiome)是特定时间特定生境中微生物群所包含的基因序列(含同源序列)的总和,而微生物群(microbiota)是特定时间特定生境所有微生物有机体的总称,其组成包括非细胞结构的病毒(包括噬菌体)、原核生物中的真细菌和古细菌,以及真核细胞微生物。目前,包括转录组分析和核糖体图谱分析在内的原核生物高通量RNA测序方法都需要深度测序以获得足够数量的有效读取,从而进一步找到低表达的遗传元件。而在测序细菌总RNA时,80%-90%的读取片段都映射到rRNA,映射到mRNA的高质量读取片段受总RNA中大量rRNA的影响(Rosenow C et al.,2001)。因此,rRNA的消解可以增加测序覆盖的mRNA范围,从而降低测序成本。
在目前已经实现的技术中,主要包含三种方法来实现rRNA消解。第一种是基于寡核苷酸与rRNA的杂交以及随后用磁珠捕获和分离(Stewart FJ etal.,2010;Kim IV etal.,2019)。第二种是基于寡核苷酸退火之后的互补和特定核酸酶对双链杂交体的消解(Yiet al.,2011;Adiconis X et al.,2013;Lyden A et al.,2019)。第三种方法是将cDNA合成限制于非rRNA模板,在cDNA合成中使用非随机且与rRNA不互补的引物(Armour CD etal.,2009)。其中第一种方法的优势在于消解的无偏差性且消解效果较好,但处理单个样本的成本较高,远高于其他两种方法,不利于推广。而第三种方法又由于引物的局限性,不能保证cDNA文库的完整性。因此,我们在第二种方法中选择了基于CRISPR-Cas9技术的rRNA消解手段进行优化。
利用CRISPR-Cas9技术来消解rRNA的基本原理是通过人工设计的sgRNA(smallguide RNA)来识别经逆转录和扩增后形成的cDNA文库,并引导Cas9蛋白酶进行有效切割rDNA双链,使得rDNA都断裂成小片段,之后经过定量的磁珠纯化去掉rDNA小片段,从而达到消解rRNA的目的。这种方法之前运用在比较单一的物种上,针对不同物种样本需要单独设计一套sgRNA(100条),这也大大提高了成本(Loi DSC et al.,2021)。而且在很多微生物没有参考基因组或者参考基因组质量较差的情况下会导致sgRNA脱靶,从而使得rRNA消解效果较差。对于含有多种微生物的复杂样本,面向单一物种设计的sgRNA明显不能满足需求。我们通过对常见的58种细菌rRNA序列分析并设计筛选了一套可以靶向细菌微生物组的sgRNA,该套sgRNA不仅可以靶向设计时参考的58种细菌的rRNA,且基于细菌rRNA序列一定程度的保守性,该套sgRNA还能靶向其他细菌rRNA。因此,我们设计的这套sgRNA能靶向几乎所有常见的细菌rRNA,弥补了之前的成本弱势,适用性也更广,大大增加了微生物组各种水平RNA测序的覆盖度。
发明内容
本发明的目的在于提供一种消解微生物组文库中rRNA的方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:一种消解微生物组文库中rRNA的方法,包括以下步骤:
S1、sgRNA设计:
1、选取常见的58种进化距离较远的微生物作为代表微生物,覆盖8个门,48属;
2、特征序列识别:从58种微生物的rRNA序列中提取针对化脓性链球菌CRISPR相关核酸酶9(SpCas9)的特征核酸序列,即长度为20nt的5′特征序列-NGG 3′特征序列,额外要求gc含量在30%-80%之间;
3、筛选出保守性较高的特征序列:通过使用blastn(Version:2.13.0)将上一步提取的特征序列与58种微生物的rRNA序列进行比对,允许1个错配,若可以映射到的物种小于等于2,则为保守性较差,滤去;
4、剔除脱靶序列:再次使用blastn,将上一步保守性较高的特征序列与58种微生物的非rRNA序列进行比对,允许3个错配,若该序列能够比对上,则说明该序列是脱靶的,滤去;
5、剔除与微液滴体系不兼容的特征序列:通过使用blastn,将上一步剔除脱靶序列后的特征序列与微液滴合成序列和barcode进行比对,允许3个错配,若该序列能够比对上,则说明该序列与微液滴体系不兼容,滤去;
6、对特征序列的保守性进行评估:建立58种微生物rRNA序列的参考基因组,使用bwa aln(Version:0.7.17-r1188)将来自多种不同环境的单细胞微生物组NGS测序数据映射到该参考基因,为上一步中的特征序列按照覆盖reads数进行降序排序,产生优先级;
7、循环初始化,令上一步中剩余的特征核酸序列为S,A为空集;
8、若S为空集循环结束,输出A,否则选取S集中优先级最高的特征核酸序列a,令S=S-{a},A=A∪{a};
9、对于序列a上下游50bp内的序列b,执行S=S-{b};
10、重复进行第8步,算法结束后,合成300条特征序列,即300条sgRNA;
S2、sgRNA制备:将由生物公司合成的sgDNA通过PCR扩增加上tracrDNA序列,之后通过琼脂糖凝胶电泳验证PCR产物条带大小在130bp,得到的产物经磁珠纯化并转录得到sgRNA,将得到的sgRNA测定浓度并稀释后待用;
S3、cDNA文库制备:将微生物RNA逆转录成cDNA,之后通过PCR扩增形成DNA双链,PCR产物经磁珠纯化后得到干净的cDNA文库,该cDNA文库测定浓度并稀释后待用;
S4、rRNA消解:取Cas9蛋白酶、Cas9蛋白酶缓冲液、sgRNA及无酶水均匀混合,在37℃金属浴中运行15min以形成sgRNA-Cas9复合体,之后加cDNA文库,继续过夜运行16h以保证Cas9蛋白酶充分酶切rDNA,加入RNA内切酶以酶解多余的sgRNA,最后加入蛋白溶解酶使Cas9蛋白失活;
S5、未酶切和酶切之后的cDNA用等体积的磁珠纯化,取少许纯化产物进行qPCR2验证,剩余纯化产物进行PCR扩增,qPCR结果显示酶切和未酶切产物的循环数存在6个左右的差异,且两者的PCR产物经过琼脂糖凝胶电泳后也显示出带型差异;
S6、上述4中PCR产物经磁珠纯化后建库测序,测序结果显示未经过酶切样本中rRNA占比明显高于经过酶切的样本中的。
优选的,在所述S3步骤中,在Cas9蛋白酶酶切时设置一个不加Cas9蛋白酶的对照组,本发明涉及一种新颖的技术方案,通过利用58种具有代表性的微生物的基因组,设计CRISPR Cas9系统的guide RNA,以达到筛选出最佳的guide RNA组合的目的。
与现有技术相比,本发明的有益效果是:
本发明提供的sgRNA设计流程及序列,参考了常见的58种细菌基因组序列,从其中rRNA序列中筛选Cas9蛋白酶能够靶向的含有5′-NGG-3′特征的核酸序列,进一步比对去除保守性较差和容易脱靶的特征序列,筛选设计出300条sgRNA,该组合旨在高效覆盖尽可能多的细菌的rRNA,同时兼顾避免对编码基因和微液滴体系中的人工序列造成脱靶效应,从而达到节省成本,扩大适用范围的效果,也大大增加了微生物组各种水平RNA测序的覆盖度。
附图说明
图1为本发明sgRNA筛选流程图;
图2为本发明未进行rRNA消解的含微生物组的样本转录组序列占比情况图;
图3为本发明rRNA消解的含微生物组的样本转录组序列占比情况图;
图4为本发明rRNA消解的含微生物组样本cDNA文库凝胶电泳图,其中,左边泳道为未进行rRNA消解的含微生物组的样本cDNA文库,中间泳道为rRNA消解的含微生物组的样本cDNA文库,右边泳道为marker;
图5为本发明qPCR验证rRNA消解效果图,其中,包括未经过rRNA消解的微生物组样本cDNA扩增曲线,Ct值为8.35,rRNA消解的微生物组样本cDNA扩增曲线,Ct值为14.18;
图6为本发明sgRNA制备的电泳图,其中,左边泳道为sgRNA,右边泳道为marker。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供一种消解微生物组文库中rRNA的方法,包括以下步骤:
S1、sgRNA设计:
1、选取常见的58种进化距离较远的微生物作为代表微生物,覆盖8个门,48属,如下表1所示;
2、特征序列识别:从58种微生物的rRNA序列中提取针对化脓性链球菌CRISPR相关核酸酶9(SpCas9)的特征核酸序列,即长度为20nt的5′特征序列-NGG 3′特征序列,额外要求gc含量在30%-80%之间;
3、筛选出保守性较高的特征序列:通过使用blastn(Version:2.13.0)将上一步提取的特征序列与58种微生物的rRNA序列进行比对,允许1个错配,若可以映射到的物种小于等于2,则为保守性较差,滤去;
4、剔除脱靶序列:再次使用blastn,将上一步保守性较高的特征序列与58种微生物的非rRNA序列进行比对,允许3个错配,若该序列能够比对上,则说明该序列是脱靶的,滤去;
5、剔除与微液滴体系不兼容的特征序列:通过使用blastn,将上一步剔除脱靶序列后的特征序列与微液滴合成序列和barcode进行比对,允许3个错配,若该序列能够比对上,则说明该序列与微液滴体系不兼容,滤去;
6、对特征序列的保守性进行评估:建立58种微生物rRNA序列的参考基因组,使用bwa aln(Version:0.7.17-r1188)将来自多种不同环境的单细胞微生物组NGS测序数据映射到该参考基因,为上一步中的特征序列按照覆盖reads数进行降序排序,产生优先级;
7、循环初始化,令上一步中剩余的特征核酸序列为S,A为空集;
8、若S为空集循环结束,输出A,否则选取S集中优先级最高的特征核酸序列a,令S=S-{a},A=A∪{a};
9、对于序列a上下游50bp内的序列b,执行S=S-{b};
10、重复进行第8步,算法结束后,合成300条特征序列,即300条sgRNA;
表1
S2、sgRNA制备:将由生物公司合成的sgDNA通过PCR扩增加上tracrDNA序列,之后通过琼脂糖凝胶电泳验证PCR产物条带大小在130bp,得到的产物经磁珠纯化并转录得到sgRNA,将得到的sgRNA测定浓度并稀释后待用;
S3、cDNA文库制备:将牙菌斑中微生物RNA逆转录成cDNA,之后通过PCR扩增形成DNA双链,PCR产物经磁珠纯化后得到干净的cDNA文库,该cDNA文库测定浓度并稀释后待用;
S4、rRNA消解:取Cas9蛋白酶、Cas9蛋白酶缓冲液、sgRNA及无酶水均匀混合,在37℃金属浴中运行15min以形成sgRNA-Cas9复合体,之后加cDNA文库,继续过夜运行16h以保证Cas9蛋白酶充分酶切rDNA,加入RNA内切酶以酶解多余的sgRNA,最后加入蛋白溶解酶使Cas9蛋白失活,在Cas9蛋白酶酶切时设置一个不加Cas9蛋白酶的对照组;
S5、未酶切和酶切之后的cDNA用等体积的磁珠纯化,取少许纯化产物进行qPCR验证,剩余纯化产物进行PCR扩增,qPCR结果显示酶切和未酶切产物的循环数存在6个左右的差异,且两者的PCR产物经过琼脂糖凝胶电泳后也显示出带型差异;
S6、上述4中PCR产物经磁珠纯化后建库测序,测序结果显示未经过rRNA消解的样本中rRNA占比高达90%,而经过rRNA消解过的样本中rRNA占比不到20%,rRNA消解效率高达98%。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (2)
1.一种消解微生物组文库中rRNA的方法,其特征在于,包括以下步骤:
S1、sgRNA设计:
1、选取常见的58种进化距离较远的微生物作为代表微生物,覆盖8个门,48属;
2、特征序列识别:从58种微生物的rRNA序列中提取针对化脓性链球菌CRISPR相关核酸酶9(SpCas9)的特征核酸序列,即长度为20nt的5′特征序列-NGG 3′特征序列,额外要求gc含量在30%-80%之间;
3、筛选出保守性较高的特征序列:通过使用blastn(Version:2.13.0)将上一步提取的特征序列与58种微生物的rRNA序列进行比对,允许1个错配,若可以映射到的物种小于等于2,则为保守性较差,滤去;
4、剔除脱靶序列:再次使用blastn,将上一步保守性较高的特征序列与58种微生物的非rRNA序列进行比对,允许3个错配,若该序列能够比对上,则说明该序列是脱靶的,滤去;
5、剔除与微液滴体系不兼容的特征序列:通过使用blastn,将上一步剔除脱靶序列后的特征序列与微液滴合成序列和barcode进行比对,允许3个错配,若该序列能够比对上,则说明该序列与微液滴体系不兼容,滤去;
6、对特征序列的保守性进行评估:建立58种微生物rRNA序列的参考基因组,使用bwaaln(Version:0.7.17-r1188)将来自多种不同环境的单细胞微生物组NGS测序数据映射到该参考基因,为上一步中的特征序列按照覆盖reads数进行降序排序,产生优先级;
7、循环初始化,令上一步中剩余的特征核酸序列为S,A为空集;
8、若S为空集循环结束,输出A,否则选取S集中优先级最高的特征核酸序列a,令S=S-{a},A=A∪{a};
9、对于序列a上下游50bp内的序列b,执行S=S-{b};
10、重复进行第8步,算法结束后,合成300条特征序列,即300条sgRNA;
S2、sgRNA制备:将由生物公司合成的sgDNA通过PCR扩增加上tracrDNA序列,之后通过琼脂糖凝胶电泳验证PCR产物条带大小在130bp,得到的产物经磁珠纯化并转录得到sgRNA,将得到的sgRNA测定浓度并稀释后待用;
S3、cDNA文库制备:将微生物RNA逆转录成cDNA,之后通过PCR扩增形成DNA双链,PCR产物经磁珠纯化后得到干净的cDNA文库,该cDNA文库测定浓度并稀释后待用;
S4、rRNA消解:取Cas9蛋白酶、Cas9蛋白酶缓冲液、sgRNA及无酶水均匀混合,在37℃金属浴中运行15min以形成sgRNA-Cas9复合体,之后加cDNA文库,继续过夜运行16h以保证Cas9蛋白酶充分酶切rDNA,之后加入RNA内切酶以酶解多余的sgRNA,最后加入蛋白溶解酶使Cas9蛋白失活;
S5、未酶切和酶切之后的cDNA用等体积的磁珠纯化,取少许纯化产物进行qPCR验证,剩余纯化产物进行PCR扩增,qPCR结果显示酶切和未酶切产物的循环数存在差异,且两者的PCR产物经过琼脂糖凝胶电泳后也显示出带型差异;
S6、上述4中PCR产物经磁珠纯化后建库测序,测序结果显示未经过酶切样本中rRNA占比明显高于经过酶切的样本中的。
2.根据权利要求1所述的一种消解微生物组文库中rRNA的方法,其特征在于:在所述S4步骤中,在Cas9蛋白酶酶切时设置一个不加Cas9蛋白酶的对照组。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310383401.9A CN116790576A (zh) | 2023-04-12 | 2023-04-12 | 一种消解微生物组文库中rRNA的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310383401.9A CN116790576A (zh) | 2023-04-12 | 2023-04-12 | 一种消解微生物组文库中rRNA的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116790576A true CN116790576A (zh) | 2023-09-22 |
Family
ID=88038775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310383401.9A Pending CN116790576A (zh) | 2023-04-12 | 2023-04-12 | 一种消解微生物组文库中rRNA的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116790576A (zh) |
-
2023
- 2023-04-12 CN CN202310383401.9A patent/CN116790576A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220389416A1 (en) | COMPOSITIONS AND METHODS FOR CONSTRUCTING STRAND SPECIFIC cDNA LIBRARIES | |
US20230056763A1 (en) | Methods of targeted sequencing | |
EP2576780B1 (en) | Method for the preparation and amplification of representative and strand- specific libraries of cdna for high throughput sequencing, use thereof, kit and cartridges for automation kit | |
CN104805099A (zh) | 一种安全编码Cas9蛋白的核酸分子及其表达载体 | |
WO2019128743A1 (zh) | 一种猪全基因组sgRNA文库及其构建方法和应用 | |
CN112680797B (zh) | 一种去除高丰度rna的测序文库及其构建方法 | |
CN109023537B (zh) | 一种微量dna样品高通量测序文库的构建技术 | |
WO2018144216A1 (en) | Methods and compositions for enrichment of target polynucleotides | |
US20170152574A1 (en) | Artificial exogenous reference molecule for type and abundance comparison among different species of microorganisms | |
CN113462748A (zh) | Dna测序文库的制备方法及试剂盒 | |
CN110218811B (zh) | 一种筛选水稻突变体的方法 | |
CN116790576A (zh) | 一种消解微生物组文库中rRNA的方法 | |
AU2021105278A4 (en) | Whole Genome High-Efficiency Gene Region Enriching and Sequencing Method | |
CN110592183A (zh) | 一种检测超低频基因突变的高通量测序文库的制备方法 | |
CN110699428B (zh) | 一种寡核苷酸库的均一化方法 | |
CN113166756B (zh) | 用于三代测序建库的融合引物、建库方法、测序方法和建库试剂盒 | |
CN110144387B (zh) | 一种多重pcr方法 | |
US20230122979A1 (en) | Methods of sample normalization | |
US20240141325A1 (en) | Generation of novel crispr genome editing agents using combinatorial chemistry | |
US20210163922A1 (en) | Assembly and error reduction of synthetic genes from oligonucleotides | |
WO2024157194A1 (en) | Methods and assays for off-target analysis | |
Crawford et al. | Consistent production of cost-effective LongSAGE libraries | |
CN115772555A (zh) | 一种利用dCas13去除核糖体RNA序列的方法 | |
CN118745425A (zh) | 用于转录组测序过程中核糖体rna去除的dna探针组合及试剂盒 | |
CN110144383A (zh) | 利用多重pcr富集目标dna片段的方法 |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Huang Yingjuan Inventor after: Ni Chengyu Inventor after: Meng Hongen Inventor after: Wu Jun Inventor before: Huang Yingjuan Inventor before: Ni Chengyu Inventor before: Wang Yongcheng Inventor before: Wu Jun Inventor before: Meng Hongen |