CN112646866A - 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 - Google Patents
一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 Download PDFInfo
- Publication number
- CN112646866A CN112646866A CN202011575802.7A CN202011575802A CN112646866A CN 112646866 A CN112646866 A CN 112646866A CN 202011575802 A CN202011575802 A CN 202011575802A CN 112646866 A CN112646866 A CN 112646866A
- Authority
- CN
- China
- Prior art keywords
- probe
- gene
- ddpcr
- mutation
- specific
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000010362 genome editing Methods 0.000 title claims abstract description 49
- 238000010453 CRISPR/Cas method Methods 0.000 title claims abstract description 23
- 206010064571 Gene mutation Diseases 0.000 title claims abstract description 19
- 231100000089 gene mutation induction Toxicity 0.000 title claims abstract description 12
- 239000000523 sample Substances 0.000 claims abstract description 87
- 230000035772 mutation Effects 0.000 claims abstract description 68
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 53
- 238000012408 PCR amplification Methods 0.000 claims abstract description 20
- 238000010791 quenching Methods 0.000 claims abstract description 8
- 230000000171 quenching effect Effects 0.000 claims abstract description 8
- 108020004414 DNA Proteins 0.000 claims description 42
- 238000006243 chemical reaction Methods 0.000 claims description 19
- 235000007164 Oryza sativa Nutrition 0.000 claims description 15
- 235000009566 rice Nutrition 0.000 claims description 13
- 241000196324 Embryophyta Species 0.000 claims description 12
- 230000036438 mutation frequency Effects 0.000 claims description 9
- 239000002773 nucleotide Substances 0.000 claims description 7
- 125000003729 nucleotide group Chemical group 0.000 claims description 7
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 239000011541 reaction mixture Substances 0.000 claims description 4
- 240000007594 Oryza sativa Species 0.000 claims description 3
- 238000000137 annealing Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 19
- 230000035945 sensitivity Effects 0.000 abstract description 7
- 241000209094 Oryza Species 0.000 description 12
- 238000003752 polymerase chain reaction Methods 0.000 description 12
- 108091033409 CRISPR Proteins 0.000 description 11
- 238000011529 RT qPCR Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 235000021067 refined food Nutrition 0.000 description 6
- 238000010354 CRISPR gene editing Methods 0.000 description 5
- 208000020584 Polyploidy Diseases 0.000 description 5
- 230000003321 amplification Effects 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 240000002791 Brassica napus Species 0.000 description 4
- 235000011293 Brassica napus Nutrition 0.000 description 4
- 238000007847 digital PCR Methods 0.000 description 4
- 239000003550 marker Substances 0.000 description 4
- 108700001094 Plant Genes Proteins 0.000 description 3
- 238000012217 deletion Methods 0.000 description 3
- 230000037430 deletion Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000011304 droplet digital PCR Methods 0.000 description 3
- 231100000350 mutagenesis Toxicity 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 108091093088 Amplicon Proteins 0.000 description 2
- 238000010356 CRISPR-Cas9 genome editing Methods 0.000 description 2
- 238000007400 DNA extraction Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005782 double-strand break Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000006780 non-homologous end joining Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000003753 real-time PCR Methods 0.000 description 2
- 238000007789 sealing Methods 0.000 description 2
- 239000008223 sterile water Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241000589158 Agrobacterium Species 0.000 description 1
- 108700028369 Alleles Proteins 0.000 description 1
- 108010008532 Deoxyribonuclease I Proteins 0.000 description 1
- 102000007260 Deoxyribonuclease I Human genes 0.000 description 1
- AHCYMLUZIRLXAA-SHYZEUOFSA-N Deoxyuridine 5'-triphosphate Chemical compound O1[C@H](COP(O)(=O)OP(O)(=O)OP(O)(O)=O)[C@@H](O)C[C@@H]1N1C(=O)NC(=O)C=C1 AHCYMLUZIRLXAA-SHYZEUOFSA-N 0.000 description 1
- 101710163270 Nuclease Proteins 0.000 description 1
- 108020005120 Plant DNA Proteins 0.000 description 1
- 108091027544 Subgenomic mRNA Proteins 0.000 description 1
- 208000035199 Tetraploidy Diseases 0.000 description 1
- 108700019146 Transgenes Proteins 0.000 description 1
- 238000000246 agarose gel electrophoresis Methods 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 210000000349 chromosome Anatomy 0.000 description 1
- 238000003776 cleavage reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000007865 diluting Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000012165 high-throughput sequencing Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001404 mediated effect Effects 0.000 description 1
- 238000011880 melting curve analysis Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000007857 nested PCR Methods 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000007480 sanger sequencing Methods 0.000 description 1
- 230000007017 scission Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
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/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
-
- 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/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
- C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
-
- 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Biochemistry (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Genetics & Genomics (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Botany (AREA)
- Mycology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供了一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法。所述方法包括根据CRISPR/Cas靶序列的突变位置，设计一对变异位点特异性PCR扩增引物和变异位点特异性探针，以及一对内源参考基因特异性PCR扩增引物和内源参考基因特异性探针；其中，所述变异位点特异性探针设置在突变位点区域，并且探针的5’端位于PAM区且连接有第一荧光基团，3’端连接有第一淬灭基团；所述内源参考基因特异性探针的5’端连接有第二荧光基团，3’端连接有第二淬灭基团。所述方法对于基因突变检测灵敏度高，对基因编辑频率的检测下限(LOD)更低，更接近预期的编辑频率值。
Description
技术领域
本发明涉及生物医学技术领域，尤其是涉及一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法。
背景技术
CRISPR/Cas系统已被广泛应用于生物体基因组编辑和其他靶向修饰。与传统的农杆菌介导的T-DNA转基因方法相比，农杆菌介导法依赖于植物的随机重组或整合，由Cas核酸酶产生的位置特异性双链断裂(double-strand break,DSB)，主要通过容易出错的非同源末端连接(NHEJ)机械进行修复。修复的结果是可能导致各种替换、插入或缺失(indels)突变，最常见的是只有1bp变异的缺失突变。
目前已报道了几种检测CRISPR/Cas系统诱导的靶基因突变的方法，包括基于聚合酶链反应(PCR)的检测、T7核酸内切酶I识别切割、高分辨率熔融曲线分析(high-resolution melting，HRM)和基于NGS的方法。但这些方法都是半定量的，很难检测到含有低初始浓度DNA的加工食品样品中的变异，而且也难以精确量化基因编辑频率。此外，对于复杂的多倍体植物基因组中非常低频率的突变，现有的方法也难以很好地检测。
微滴式数字PCR(ddPCR)是一项突破性的技术，它依赖于将单个扩增分为不同的隔间，并检测其末端扩增产物。它提供超灵敏和绝对的核酸定量，没有标准曲线。到目前为止，几乎很少有关于ddPCR技术检测基因编辑植物及其衍生的加工食品相关报道。
鉴于，提出本发明。
发明内容
本发明的目的在于提供一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法。解决了现有技术难以定量检测CRISPR/Cas诱导的植物基因突变以及难以准确量化植物基因编辑频率的技术问题。
本发明提供的技术方案如下：
一种基于ddPCR检测CRISPR/Cas诱导的基因突变的方法，所述方法包括：
根据CRISPR/Cas靶序列的突变位置，设计一对变异位点特异性PCR扩增引物和变异位点特异性探针，以及一对内源参考基因特异性PCR扩增引物和内源参考基因特异性探针；
其中，所述变异位点特异性探针设置在突变位点区域，并且探针的5’端位于PAM区且连接有第一荧光基团，3’端连接有第一淬灭基团；所述内源参考基因特异性探针的5’端连接有第二荧光基团，3’端连接有第二淬灭基团。所述第一荧光基团和所述第二荧光基团不同。
变异位点特异性PCR扩增引物设计跨越突变位点。使用普通PCR筛选引物，确定可以扩增单一的，片段大小正确的产物。CRISPR/Cas9诱导的突变是可预测的，大多数突变发生在原间隔基序邻接基序(PAM)5’端上游的3个碱基处。变异位点特异性探针应位于PAM区。为了保持探针对突变的敏感性，最好将PAM区标记为探针的5’端。
内源性参考基因是已验证的参考基因，内源性参考基因的检测区域不含基因编辑的突变位点，所选择的内源性参考基因拷贝数需要与所编辑的目标基因拷贝数在该植物样本中相同。基因编辑引起的突变将会导致变异位点特异性探针无法与目标序列结合，但不会影响内源性参考基因探针的结合，变异位点特异性探针主要用来评估突变基因的数量，而内源性参考基因特异性探针主要用来分析样本中的等位基因总数量。
在一个实施方案中，所述方法还包括：采用所述的变异位点特异性PCR扩增引物和变异位点特异性探针以及内源参考基因特异性PCR扩增引物和内源参考基因特异性探针进行ddPCR反应；
依据荧光微滴的荧光信号得出结果：第一荧光和第二荧光两种荧光信号均为阳性的微滴为野生型微滴，第二荧光信号阳性但第一荧光信号阴性的微滴为突变型微滴。
野生型微滴和纯合突变型微滴可以通过ddPCR分析的二维视图明确区分。
一种基于ddPCR检测CRISPR/Cas诱导的基因编辑频率的方法，所述方法在根据前述的方法获得荧光微滴的荧光信号结果的基础上，计算突变型微滴拷贝数与野生型微滴拷贝数的比值来计算基因编辑突变频率。
在一个实施方案中，在上述基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法中，所述第一淬灭基团选自BHQ或MGB；所述第一荧光基团和所述第二荧光基团选自FAM、HEX和TEX中的任一种；优选地，所述第一荧光基团为FAM，所述第二荧光基团为HEX。
内源参考基因探针为5’端HEX标记，突变位点特异性探针为5’端FAM标记。
在一个实施方案中，所述CRISPR/Cas靶序列为水稻TGW6基因(Os06g0623700)的一段序列，其核苷酸序列如SEQ ID NO.17所示。
在一个实施方案中，所述变异位点特异性PCR扩增引物的上下游引物序列分别如SEQ ID NO.1和SEQ ID NO.2所示，所述变异位点特异性探针的序列如SEQ ID NO.3所示。
在一个实施方案中，所述内源参考基因特异性PCR扩增引物的上下游引物序列分别如SEQ ID NO.4和SEQ ID NO.5所示，所述内源参考基因特异性探针的序列如SEQ IDNO.6所示。
用于ddPCR技术检测CRISPR/Cas诱导的水稻TGW6基因突变和基因编辑频率的引物对和探针，所述引物对由上述变异位点特异性PCR扩增引物和变异位点特异性探针以及内源参考基因特异性PCR扩增引物和内源参考基因特异性探针的序列组成。
在一个实施方案中，所述方法还包括提取植物模板DNA，然后将所述模板DNA与所述的变异位点特异性PCR扩增引物和变异位点特异性探针以及内源参考基因特异性PCR扩增引物和内源参考基因特异性探针一起制备反应混合物，使用所述反应混合物进行ddPCR反应。
在一个实施方案中，所述ddPCR反应体系中，植物模板DNA的量为0.116-10ng/μL。
在一个实施方案中，所述ddPCR反应中，引物的浓度为400-600nM，优选地，引物的浓度为450-500nM；探针的浓度为200-400nM；优选地，探针的浓度为250-300nM。
在一个具体的实施方案中，所述ddPCR反应中，引物的浓度为450nM；探针的浓度为250nM。
在一个实施方案中，所述ddPCR反应体系中，还包含ddPCR SuperMix和无菌水。将PCR反应液放入微滴生成卡中，加入微滴生成油在微滴发生器中生成微滴。将微滴转移到96孔PCR板上进行PCR反应。
在一个实施方案中，所述ddPCR反应的程序为93℃-96℃5-15min，93℃-95℃5-15s，58℃或68℃退火和延伸50-70s，35-40个循环，95-100℃下8-12min终止反应。
在一个具体的实施方案中，所述ddPCR反应的程序为95℃10min，94℃10s，58℃或68℃退火和延伸60s，40个循环，98℃下10min终止反应，冷却至4℃。
在一个实施方案中，所述方法具体包括以下步骤：
(1)提取基因编辑植物DNA，获得DNA模板；
(2)将杂合子基因编辑DNA样本稀释为不同的浓度梯度，优选地，DNA浓度分别为10、5、0.4、0.08、0.016ng/μL；
(3)根据突变位点，设计变异位点特异性PCR扩增引物和变异位点特异性探针；探针的5’端位于PAM区且连接有第一荧光基团(优选FAM)，3’端连接有第一淬灭基团；选择内源参考基因的特异性引物和探针，探针的5’端连接有第二淬灭基团(优选HEX)；
(4)进行微滴式数字PCR并进行数据分析，分析了突变型微滴和野生型微滴的浓度并确定阈值以区分阳性微滴和阴性微滴；
(5)通过突变型微滴(如只有HEX阳性微滴)与野生型微滴(如HEX/FAM双阳性微滴)的比值确定基因编辑的突变频率。
有益效果：
本发明提供的基于微滴式数字PCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法，对于基因编辑技术诱导的不同类型基因突变，该方法灵敏性较高。与基于qPCR和NGS的方法相比，ddPCR方法对编辑频率的检测下限(LOD)更低，更接近预期的编辑频率值。
基于微滴式数字PCR检测CRISPR/Cas诱导的基因突变具有更高的准确度和精密度，可以准确检测基因编辑突变体的突变频率，将成为检测和评价植物基因编辑频率的有力工具。
本发明提供的方法可以检测到含有低初始浓度DNA的加工食品样品中的基因突变，同时适用于检测复杂的多倍体植物基因组中低频率的基因突变，具有更低的错误率。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的引物和探针设计原理图以及通过ddPCR检测野生型水稻样本、纯合型和杂合型基因突变的液滴二维分布图；
图2为本发明实施例提供的利用ddPCR检测经Sanger测序验证的不同突变类型纯合子水稻样品的结果图；
图3为本发明实施例提供的ddPCR法检测相同基因编辑样本不同初始DNA浓度的灵敏性测定结果图；
图4为本发明实施例提供的不同方法(ddPCR、qPCR和NGS)分析基因编辑频率的检测结果比较；
图5为ddPCR法验证多倍体生物(甘蓝型油菜)中的基因编辑频率的结果图。
具体实施方式
下面将结合实施例对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
1.植物材料
CRISPR/Cas9基因编辑水稻突变体，来源于发明人团队利用CRISPR/Cas9基因编辑，自主研发创制产生的突变株系编号为1-8。
利用网络工具CRISPR-P设计sgRNA，其靶位点区域核苷酸序列如SEQ ID NO.17所示
(CCTCGTGCTGCTCTCGCCGTCCCCTACTGCCGCCGCCACAGCCACAACGAGAATGTTCAA，其中TCC为PAM位点)。
基因编辑油菜株系(S1-14，S1-18，S1-24，S1-53，S1-104)为CRISPR/Cas9基因编辑作物，其靶位点区域核苷酸序列如SEQ ID NO.18所示(GGCCACTTGCCCATGCTGGCTAAGTACGGCCCTGAC，其中CCA为PAM位点)。由中国农科院油料所郑明博士提供。
2.DNA提取
利用天根试剂盒提取植物DNA，利用NanoDrop 1000测定DNA的浓度，琼脂糖凝胶电泳进一步鉴定DNA提取的完整性。所有DNA模板均稀释至10ng/μL，所有的DNA模板均在-20℃冰箱里保存。
3.样品制备
为了分析低初始浓度的DNA，杂合子基因编辑DNA样本用水稀释，并使用
2.0(美国生命技术公司)测量，DNA浓度分别为10、5、0.4、0.08、0.016ng/μL。
4.引物和探针设计
CRISPR/Cas9诱导的大多数突变发生在原间隔基序邻接基序(PAM)5'端上游的3个碱基处。在微滴式数字PCR系统中，设计两对引物(针对内源性参考基因的引物对和针对突变位点特异性基因的引物对)和各自的检测探针(内源性参考基因的探针和突变位点特异性探针)(见图1中A图)：其中一对引物，扩增区域包含突变位点，探针位于突变位点区域并由突变位点特异探针(FAM修饰)进行检测；另一对引物和探针可以选择已被验证的内源性参考基因(HEX修饰)，该检测区域不含基因编辑的突变位点。所选择的内源性参考基因拷贝数需要与所编辑的目标基因拷贝数在该植物样本中相同。利用突变型微滴(只有FAM阳性液滴)拷贝数与野生型微滴(HEX/FAM双阳性微滴)拷贝数的比值来计算基因编辑的突变频率。
使用Primer Express软件3.0设计突变位置的引物和探针。设计还要遵循以下原则：(1)引物跨越突变位点，(2)探针应位于PAM区，(3)为了保持探针对突变的敏感性，PAM区域标记为探针的5'端。使用普通PCR筛选引物，确定可以扩增单一的，片段大小正确的产物。
参考基因探针为5′端HEX标记，突变位点特异性探针为5′端FAM标记。
下表1列出了本研究中使用的引物和探针。
表1 ddPCR中所用的引物及荧光探针
a用MGB修饰的探针，*作为内参基因
5.微滴式数字PCR
液滴式数字PCR检测包括以下成分：
20μL总反应体系中：10μL ddPCR SuperMix(无dUTP)(Bio-Rad)，每个引物对450nM(引物对包括：针对内源性参考基因的引物对和针对突变位点特异性基因的引物对)，每个探针250nM(探针包括：内源性参考基因的探针和突变位点特异性探针)，加无菌水调节最终体积至19μL，然后加入1μL DNA模板。
将20μL混合物放入BioRad DG8的微滴生成卡中，并加入70μL微滴生成油。将微滴生成卡放入QX200微滴发生器(Bio-Rad)中以生成微滴。将微滴转移到96孔PCR板上。在用铝箔密封热封后，将PCR板置于7500实时PCR系统中，并在以下循环条件下进行扩增：95℃10min，40个循环94℃10s，58或68℃退火和延伸60s，在98℃下10min终止反应，并冷却至4℃。扩增后，将96孔PCR板放入QX200微滴阅读器(Bio-Rad)中进行数据分析。
野生型和纯合突变型可以通过ddPCR分析的二维视图明确区分。
利用Bio-Rad QuantaSoftTM分析了突变型微滴和野生型微滴的浓度并确定阈值以区分阳性微滴和阴性微滴。基因编辑的突变频率通过突变微滴(只有HEX阳性微滴)与野生型微滴(HEX/FAM双阳性微滴)的比值确定。
6.实时定量PCR
根据说明，使用FastStart Universal Probe Master(Roche)和ROX参考染料，在7500实时PCR仪(Life Technologies AB)上进行定量实时PCR。
7.基于NGS的测序
设计PCR引物(表2)用于使用套式PCR方法扩增目标突变位点两侧的产物，并在PCR产物上附加独特的样本特异性条形码。高通量测序由中国北京诺和基因生物信息研究所使用Illumina HiSeq平台(美国Illumina)进行。最初使用
2.0(Life Technologies，美国)测量库的浓度。利用hitom程序对高通量突变序列解码进行基因编辑突变频率进行分析(http://www.hi-tom.net/hi-tom/)。
表2本发明中NGS所用的引物
结果分析：
1、基于ddPCR检测基因编辑水稻
用水稻样本(水稻TGW6野生型基因、CRISPR-Cas9基因编辑诱导的纯合子突变型和CRISPR-Cas9基因编辑诱导的杂合子突变型)对基于ddPCR的基因编辑植物突变检测方法的性能进行了评价。
使用表1中所设计的针对水稻(Oryza sativa)TGW6基因(Os06g0623700)突变位点设计的突变位点特异性基因的引物对(SEQ ID NO.1和SEQ ID NO.2)以及突变位点特异性探针(SEQ ID NO.3)以及表1中所设计的针对水稻(Oryza sativa)内源性参考基因(LOC_Os2g42314)的引物对(SEQ ID NO.4和SEQ ID NO.5)以及针对内源性参考基因的探针(SEQID NO.6)进行了ddPCR检测，结果显示，三个水稻样本能够通过ddPCR分析的二维视图清晰地区分(图1中B图所示)。
在ddPCR分析的二维视图中，含有两种荧光信号的微滴是野生型扩增子，而含有HEX阳性但FAM阴性信号的微滴是突变型扩增子。杂合子突变的扩增包含突变微滴或野生型微滴。
2、ddPCR平台能有效地识别基因编辑诱导的不同突变类型
鉴于Cas产生的突变通常是只有1bp变异，因此新开发的方法能够有效地检测这些微小的基因突变显得尤为非常重要。
为了验证ddPCR方法是否能对不同类型的突变(如单核苷酸指数)敏感，我们利用ddPCR平台检测了经测序证实的不同类型纯合子突变的CRISPR/Cas9基因编辑的TGW6水稻样品。这些突变不仅包含单核苷酸的缺失，而且包含单核苷酸突变(图2中A图)。
ddPCR结果表明，在所有这些样品中，微滴含有HEX阳性信号，而FAM阴性信号，并且当有数千个HEX阳性微滴时，FAM阳性微滴的数量几乎为零(图2中B、C图)。因此，基于ddPCR的检测方法可以有效地识别基因编辑引起的各种突变。
3、基因编辑样本的灵敏性测定
为了评价ddPCR检测基因编辑样本的灵敏性，我们使用低初始浓度的DNA样本，甚至是加工过的食品样本来评估ddPCR的性能。
制备了一系列DNA浓度的杂合子样品(10、5、0.4、0.08、0.016ng/μL)使用ddPCR进行检测。
结果表明，ddPCR法可准确检测低至0.08ng/μL的样本，而0.016ng/μL的DNA样本由于背景噪音太大而无法识别。对于加工食品样品，ddPCR方法也表现出良好的检测性能(图3中A图)，说明ddPCR具有广泛的应用前景。
为了确定ddPCR方法突变频率的检测限(LOD)，我们将突变模板(纯合突变DNA)和野生型模板以不同的比例混合(突变DNA在50％-0.1％之间)。随着突变模板含量的减少，混合样品中突变微滴的浓度逐渐降低，而野生型微滴的浓度基本保持不变(图3中B图)。
结果表明，基于ddPCR检测方法可以在0.1％的突变频率下检测到突变模板。因此，ddPCR是一种超灵敏的基因编辑检测和评价方法，可用于检测低初始浓度DNA的样品，如加工食品样品。
4、ddPCR与qPCR和NGS方法的比较
目前，基于qPCR和NGS的方法已经被证明可以用于定量检测由基因编辑诱导的植物突变频率。我们进一步使用qPCR和NGS方法检测和分析相同的混合样品。
与ddPCR方法相比，基于qPCR和NGS的方法的结果显示，当突变模板率低于5％(图4中A-C，表3)时，没有任何阳性信号，表明qPCR和基于NGS的方法的LOD约为5％(图4中B-C)。ddPCR方法检测到的结果与预期编辑频率呈线性关系(R2>0.999)，而基于qPCR和基于NGS的方法分别为R2＝0.982和0.998。此外，通过使用基于NGS的方法，我们发现一些混合样本(表3)的编辑目标区域存在意外的不匹配(T->C)，这可能是由NGS样本制备引起的。总之，与其他方法相比，ddPCR为基因编辑频率的检测和评估提供了一种更准确、更定量的方法。
表3.NGS法检测基因编辑频率
5、利用ddPCR平台检测同源四倍体油菜基因编辑突变
为了进一步证实ddPCR法是否能准确检测多倍体生物中的基因编辑频率，选择了甘蓝型油菜(Brassica napus L.，AACC，2n＝38)样品进行分析。我们根据先前已建立的基因编辑油菜，在不同染色体上有基因编辑靶点，其中一个靶点是杂合子(A基因组或C基因组)，而另一个是Sanger测序的纯合突变或野生型(图5中A图)，理论上这些样本，一个基因组纯合突变，另一个杂合子突变的编辑频率应为75％或一个基因组纯合野生型，另一个杂合子突变的编辑频率为25％。当使用ddPCR分析，正如预期的那样，S1-14(Aacc)、S1-18(aaCc)和S1-24(Aacc)的编辑频率为75％，S1-53(AaCC)和S1-104(AACc)的编辑频率为25％，这与之前的结果一致(图5中B图)，这表明我们的ddPCR方法在检测多倍体生物突变方面的稳定性。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
SEQUENCE LISTING
<110> 浙江省农业科学院
<120> 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法
<130> PA20034568
<160> 18
<170> PatentIn version 3.3
<210> 1
<211> 19
<212> DNA
<213> 人工序列
<400> 1
tcgcgctcat tgtcttcct 19
<210> 2
<211> 23
<212> DNA
<213> 人工序列
<400> 2
tggtcttgaa cattctcgtt gtg 23
<210> 3
<211> 14
<212> DNA
<213> 人工序列
<400> 3
ctactgccgc cgcc 14
<210> 4
<211> 20
<212> DNA
<213> 人工序列
<400> 4
ccttcggaga caccttttga 20
<210> 5
<211> 20
<212> DNA
<213> 人工序列
<400> 5
ttgaaatgca cattcgggtg 20
<210> 6
<211> 21
<212> DNA
<213> 人工序列
<400> 6
ctccttcctc cgcaagttcg c 21
<210> 7
<211> 21
<212> DNA
<213> 人工序列
<400> 7
cgaccttcct ggtccgtact c 21
<210> 8
<211> 20
<212> DNA
<213> 人工序列
<400> 8
gcttggcaag aacggagaag 20
<210> 9
<211> 15
<212> DNA
<213> 人工序列
<400> 9
ttgcccatgc tggct 15
<210> 10
<211> 18
<212> DNA
<213> 人工序列
<400> 10
ggccagggct tccgtgat 18
<210> 11
<211> 21
<212> DNA
<213> 人工序列
<400> 11
ccgtcgttgt agaaccattg g 21
<210> 12
<211> 29
<212> DNA
<213> 人工序列
<400> 12
agtccttatg tgctccactt tctggtgca 29
<210> 13
<211> 40
<212> DNA
<213> 人工序列
<400> 13
ggagtgagta cggtgtgcag cgtacgtgca gtgcagccaa 40
<210> 14
<211> 40
<212> DNA
<213> 人工序列
<400> 14
gagttggatg ctggatggag tgatccgccg aggtccagat 40
<210> 15
<211> 40
<212> DNA
<213> 人工序列
<400> 15
ggagtgagta cggtgtgcta ctcaatgtcc ccggtccaac 40
<210> 16
<211> 43
<212> DNA
<213> 人工序列
<400> 16
gagttggatg ctggatggaa ctgaaccaaa atccaaatat ccc 43
<210> 17
<211> 60
<212> DNA
<213> 人工序列
<400> 17
cctcgtgctg ctctcgccgt cccctactgc cgccgccaca gccacaacga gaatgttcaa 60
<210> 18
<211> 38
<212> DNA
<213> 人工序列
<400> 18
ggccacttgc ccatgctggc taagtacggc cctgac 36
Claims (10)
1.一种基于ddPCR检测CRISPR/Cas诱导的基因突变的方法，其特征在于，所述方法包括：
根据CRISPR/Cas靶序列的突变位置，设计一对变异位点特异性PCR扩增引物和变异位点特异性探针，以及一对内源参考基因特异性PCR扩增引物和内源参考基因特异性探针；
其中，所述变异位点特异性探针设置在突变位点区域，并且探针的5’端位于PAM区且连接有第一荧光基团，3’端连接有第一淬灭基团；所述内源参考基因特异性探针的5’端连接有第二荧光基团，3’端连接有第二淬灭基团，所述第一荧光基团和所述第二荧光基团不同。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
采用权利要求1的引物和探针进行ddPCR反应；
依据荧光微滴的荧光信号得出结果：第一荧光和第二荧光两种荧光信号均为阳性的微滴为野生型微滴，第二荧光信号阳性但第一荧光信号阴性的微滴为突变型微滴。
3.一种基于ddPCR检测CRISPR/Cas诱导的基因编辑频率的方法，其特征在于，所述方法在根据权利要求2所述的方法获得荧光微滴的荧光信号结果的基础上，计算突变型微滴拷贝数与野生型微滴拷贝数的比值来计算基因编辑突变频率。
4.根据权利要求1或2所述的方法，其特征在于，所述第一淬灭基团选自BHQ或MGB；所述第一荧光基团和所述第二荧光基团选自FAM、HEX和TEX中的任一种；优选地，所述第一荧光基团为FAM，所述第二荧光基团为HEX。
5.根据权利要求1-3任一项所述的方法，其特征在于，所述CRISPR/Cas靶序列为水稻TGW6基因的一段序列，其核苷酸序列如SEQ ID NO.17所示。
6.根据权利要求5所述的方法，其特征在于，所述变异位点特异性PCR扩增引物的上下游引物序列分别如SEQ ID NO.1和SEQ ID NO.2所示，所述变异位点特异性探针的序列如SEQ ID NO.3所示。
7.根据权利要求5所述的方法，其特征在于，所述内源参考基因特异性PCR扩增引物的上下游引物序列分别如SEQ ID NO.4和SEQ ID NO.5所示，所述内源参考基因特异性探针的序列如SEQ ID NO.6所示。
8.根据权利要求7所述的方法，其特征在于，所述方法还包括提取植物模板DNA，然后将所述模板DNA与变异位点特异性PCR扩增引物和变异位点特异性探针以及内源参考基因特异性PCR扩增引物和内源参考基因特异性探针一起制备反应混合物，使用所述反应混合物进行ddPCR反应。
9.根据权利要求8所述的方法，其特征在于，所述ddPCR反应中，引物的浓度为400-600nM，优选地，引物的浓度为450-500nM；探针的浓度为200-400nM；优选地，探针的浓度为250-300nM。
10.根据权利要求8所述的方法，其特征在于，所述ddPCR反应的程序为93℃-96℃5-15min，93℃-95℃5-15s，55℃-60℃退火和延伸50-70s，35-40个循环，95-100℃下8-12min终止反应。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011575802.7A CN112646866A (zh) | 2020-12-28 | 2020-12-28 | 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011575802.7A CN112646866A (zh) | 2020-12-28 | 2020-12-28 | 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112646866A true CN112646866A (zh) | 2021-04-13 |
Family
ID=75363233
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011575802.7A Pending CN112646866A (zh) | 2020-12-28 | 2020-12-28 | 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112646866A (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106868137A (zh) * | 2017-03-01 | 2017-06-20 | 浙江省农业科学院 | 转基因水稻多重数字pcr定量检测方法 |
| CN112029838A (zh) * | 2020-07-23 | 2020-12-04 | 东南大学 | 一种用于DNA均相检测的CRISPR/Cas9分型PCR方法及其应用 |
-
2020
- 2020-12-28 CN CN202011575802.7A patent/CN112646866A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106868137A (zh) * | 2017-03-01 | 2017-06-20 | 浙江省农业科学院 | 转基因水稻多重数字pcr定量检测方法 |
| CN112029838A (zh) * | 2020-07-23 | 2020-12-04 | 东南大学 | 一种用于DNA均相检测的CRISPR/Cas9分型PCR方法及其应用 |
Non-Patent Citations (1)
| Title |
|---|
| CHENG PENG等: "Accurate Detection and Evaluation of the Gene-Editing Frequency in Plants Using Droplet Digital PCR", 《 FRONT. PLANT SCI.》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104894271B (zh) | 一种检测基因融合的方法及装置 | |
| Peng et al. | Accurate detection and evaluation of the gene-editing frequency in plants using droplet digital PCR | |
| CN111052249B (zh) | 确定预定染色体保守区域的方法、确定样本基因组中是否存在拷贝数变异的方法、系统和计算机可读介质 | |
| Loftus et al. | Development and validation of InnoQuant® HY, a system for quantitation and quality assessment of total human and male DNA using high copy targets | |
| Kim et al. | A real-time PCR-based amelogenin Y allele dropout assessment model in gender typing of degraded DNA samples | |
| CN111647649B (zh) | 一种基于ccdc39基因cnv检测辅助选择黄牛生长性状的方法 | |
| CN105316418B (zh) | 用于检测肉制品中鸭源性成分的特异性引物、探针、试剂盒及其检测方法 | |
| CN108866155B (zh) | 一种下一代测序文库的制备方法 | |
| CN108359720A (zh) | 突变体的检测方法及其应用 | |
| CN116863997A (zh) | 一种基于动物模型三黄鸡选育系统及方法 | |
| CN105695581B (zh) | 一种基于二代测试平台的中通量基因表达分析方法 | |
| CN112646866A (zh) | 一种基于ddPCR检测CRISPR/Cas诱导的基因突变和基因编辑频率的方法 | |
| Mitrečić et al. | An improved method for determination of gene copy numbers in transgenic mice by serial dilution curves obtained by real-time quantitative PCR assay | |
| CN111206076B (zh) | 一种基于片段分析技术的拷贝数变异通用验证方法 | |
| CN108384841A (zh) | 一种高通量区分转基因植物纯合体和杂合体的方法 | |
| CN103589797A (zh) | 一种单核苷酸多态性分型的方法及其应用 | |
| He et al. | An efficient and accurate droplet digital PCR method for rapid transgene copy number detection and homozygous identification in cotton (Gossypium hirsutum) | |
| Barbisin et al. | Quantifiler® Duo DNA Quantification Kit: a guiding tool for short tandem repeat genotyping of forensic samples | |
| EP4289966A1 (en) | Genetic information analysis system and genetic information analysis method | |
| Wu et al. | CTREP-finder: A web service for quick identification and visualization of clean transgenic and genome-edited plants | |
| STRES | The first decade of terminal restriction fragment length polymorphism (T-RFLP) in microbial ecology | |
| RU2817377C1 (ru) | Способ выявления событий редактирования гена gbssi у зерновых культур с помощью набора олигонуклеотидных последовательностей | |
| CN114807407B (zh) | 一种检测大豆转基因品系的引物对组合、试剂盒及检测方法 | |
| CN114369678B (zh) | 一种检测甜菜转基因成分和转基因品系的引物组合、试剂盒、检测方法及应用 | |
| CN116732227A (zh) | 转基因作物目的基因gat-ms拷贝数检测方法及试剂盒 |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210413 |
|
| RJ01 | Rejection of invention patent application after publication |