CN107904298A - 一种用于分析肠道微生物的试剂盒及其应用 - Google Patents

一种用于分析肠道微生物的试剂盒及其应用 Download PDF

Info

Publication number
CN107904298A
CN107904298A CN201711478123.6A CN201711478123A CN107904298A CN 107904298 A CN107904298 A CN 107904298A CN 201711478123 A CN201711478123 A CN 201711478123A CN 107904298 A CN107904298 A CN 107904298A
Authority
CN
China
Prior art keywords
sequence
primer
seq
sequencing
analysis
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
Application number
CN201711478123.6A
Other languages
English (en)
Inventor
朱永亮
穆延召
周燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Puruisen Gene Technology Co Ltd
Original Assignee
Suzhou Puruisen Gene Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Puruisen Gene Technology Co Ltd filed Critical Suzhou Puruisen Gene Technology Co Ltd
Priority to CN201711478123.6A priority Critical patent/CN107904298A/zh
Publication of CN107904298A publication Critical patent/CN107904298A/zh
Pending legal-status Critical Current

Links

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/6869—Methods for sequencing

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Immunology (AREA)
  • Biotechnology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明属于分子生物学领域，涉及一种用于分析肠道微生物的试剂盒及其应用，具体涉及通过两轮PCR构建文库，分析分类肠道微生物的试剂盒及应用，本发明引物组合物采用两步法扩增获得测序文库，其中包含了随机碱基序列，增加了文库序列的复杂度，降低了phix的比例，保证更多测序质量，提高了文库的覆盖度。

Description

一种用于分析肠道微生物的试剂盒及其应用
技术领域
本发明属于分子生物学领域，涉及一种用于分析肠道微生物的试剂盒及其应用，具体涉及通过两轮PCR构建文库，包括分析分类肠道微生物的引物组合物的试剂盒及应用。
背景技术
动物胃肠道内庞大多样的微生物群落与动物的食性、机体的免疫功能、疾病与健康等有着密切联系，越来越引起重视并成为研究热点。然而采用传统方法破环性的采集胃肠道样品对其进行研究受到很大制约，且对于野生保护动物并不可取，因此很多研究都以动物的粪便作为研究样品。
科学家们还发现，不同国家人群中的肠道微生物存在很多差异，差异最显著的是微生物的多样性程度。例如：印第安人和马拉维人的肠道微生物组比美国人的肠道微生物具有更大的多样性。有观点认为，微生物多样性程度越高，人体越健康。该研究还发现，尽管来自三种不同地域人群的肠道微生物组存在很多差异，但它们之间也存在惊人的相似性。如，三个不同国家的婴儿微生物组形成过程具有共同的模式，即婴儿需要6-9个月的时间来获得第一组6-700个细菌，然后再经过几年的时间才能获得成人的微生物组。该研究还发现了一个非常有趣的现象，即肠道微生物组的构成会随年龄增长而发生改变，而这一变化恰恰适应了不同年龄段人体的需求。
肠道微生物的研究大体经历了培养依赖的方法，非培养依赖的传统分子生物学方法，基于测序的高通量组学方法3个阶段。
通过比较3个阶段研究方法的特点与应用，人们可以发现培养依赖的方法在鉴定菌种的同时即可获得相应菌株，方便后续研究，但培养法本身费时费力，肠道微生物以厌氧菌和兼性厌氧菌为主，培养起来更加困难，并且在培养的过程中菌种比例会发生改变，使得其应用存在瓶颈。非培养依赖的传统分子生物学方法可以不经培养直接从样品中提取肠道微生物基因组，利用分子生物学手段进行分离、鉴定和定量，使其结果可以比较准确的反应肠道微生物中高丰度菌种的组成和真实比例。
但是传统分子生物学方法存在着通量低的缺陷，靶向方法如实时定量PCR每次只能研究一种或一类肠道微生物，而非靶向方法如最常用的梯度变形凝胶电泳(Denaturinggradient gel electrophoresis，DGGE)受限于灵敏度，往往只能研究肠道中高丰度的微生物。肠道微生物群落组成复杂，每种细菌都与其他细菌形成复杂的相互关系网络，低丰度的菌种同样扮演者重要的角色，所以传统分子生物学研究结果往往具有片面性。
16SrDNA因其序列在物种间的高度多样性，成为细菌分类学研究的“分子钟”，具9个可变区和10个保守区间隔排列的特征，可变区因细菌而异，且变异程度与细菌的***发育密切相关。因此通过分析可变区的序列即可得到各细菌的分类学特征。
CN 105937053 A公开了一种基于高通量基因测序建立粪便菌群基因文库的方法，该方法采用“巢式PCR”的方法对以16S rDNA进行富集扩增，最大限度减少了宿主和食物残渣基因组污染，同时将V3和V6组合，进行特异扩增并进行大规模并行测序，获得菌群的靶标基因序列库。CN 107058490 A公开了一种基于新一代高通量测序技术的肠道及口腔菌群多样性及差异性的分析方法，该方法利用illumina MiSeq测序***进行测序，从而分析菌群的组成与功能。但构建的文库覆盖率不够高，检测灵敏度较低。
基于上述问题，如何设计引物，能够尽可能全面覆盖所述样本的微生物，提高检测的特异性和灵敏度，就成为亟待解决的问题。
发明内容
针对现有技术的不足，本发明提供了一种分析肠道微生物的试剂盒及其应用，本发明引物组合物采用两步法扩增获得测序文库，其中包含了1-7bp的随机碱基序列，增加了文库序列的复杂度，降低了phix的比例，保证更多测序质量，提高了文库的覆盖度。
第一方面，本发明提供了一种用于分析肠道微生物的试剂盒，其包含用于构建肠道微生物文库的引物组合物，所述引物组合物包括两组引物，具体如下：
第一组引物，所述第一组引物包含第一正向引物和第一反向引物，其中，所述第一正向引物从5’到3’端包括依次相连的三个部分：第一测序通用引物、1-7bp的随机碱基序列和肠道微生物16S rRNA V3-V4区特异性正向引物；所述第一反向引物从5’到3’端包括依次相连的三个部分：第二测序通用引物、1-7bp的随机碱基序列和肠道微生物16S rRNA V3-V4区特异性反向引物；
第二组引物，所述第二组引物包含第二正向引物和第二反向引物，其中，所述第二正向引物从5’到3’端包括依次相连的三个部分：P5端接头序列、第一标签序列和第一测序通用引物的部分序列；所述第二反向引物从5’到3’端包括依次相连的三个部分：P7端接头序列、第二标签序列和第二测序通用引物的部分序列。
本发明中，发明人发现，通过采用两步法扩增获得测序文库，所述第一组引物包含V3-V4位点特异性简并引物，提高了扩增片段覆盖率、检测特异性及灵敏度，结合数据分析，可获得更丰富的微生物组信息，此外，所述第一组引物包含随机碱基，增加了文库序列的复杂度，降低了phix的比例，保证更多测序质量；第二组引物包含双index(标签)的设计，可以更大程度地混合更多的样本。
根据本发明，所述第一测序通用引物和所述第二测序通用引物独立地选自任意一种二代测序平台的测序通用引物，本发明采用illumina平台的测序通用引物。
根据本发明，所述第一测序通用引物的核酸序列如SEQ ID NO.1所示，所述第二测序通用引物的核酸序列如SEQ ID NO.2所示，所述核酸序列如下：
第一测序通用引物(SEQ ID NO.1)：TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG；
第二测序通用引物(SEQ ID NO.2)：GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG.
根据本发明，所述随机碱基序列为1-7个随机碱基，如(N)1-7，优选为3-5个随机碱基，如(N)3-5，其中，N表示为A\T\G\C中的任意一个碱基。
发明人发现，通过随机碱基序列可以进一步增加文库的复杂性，提升扩增后的文库的质量，随机碱基不会对所述引物扩增的文库造成影响，选取任意的随机碱基都可以构建文库。
根据本发明，所述肠道微生物16S rRNA V3-V4区特异性正向引物的核酸序列如SEQ ID NO.3所示，所述肠道微生物16S rRNA V3-V4区特异性反向引物的核酸序列如SEQID NO.4所示，具体核酸序列如下：
肠道微生物16S rRNA V3-V4区特异性正向引物(SEQ ID NO.3)：CCTACGGRRBGCASCAGKVRVGAAT；
其中，R选自A/G中的任意一个碱基，B选自T/C/G中的任意一个碱基，S选自C/G中的任意一个碱基，K选自T/G中的任意一个碱基，V选自A/C/G中的任意一个碱基；
肠道微生物16S rRNA V3-V4区特异性反向引物(SEQ ID NO.4)：GGACTACNVGGGTWTCTAATC；
其中，N选自A\G\C\T中的任意一个碱基，V选自A/C/G中的任意一个碱基，W选自A/T中的任意一个碱基。
本发明中，所述第一正向引物通过SEQ ID NO.1所示的第一测序通用引物+1-7bp的随机碱基+SEQ ID NO.3所示的肠道微生物16S rRNA V3-V4区特异性正向引物组合得到，为：5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-(N)1-7-CCTACGGRRBGCASCAGKVRVGAAT-3'，例如可以是5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-ATGCA-CCTACGGRRBGCASCAGKVRVGAAT-3'；所述第一反向引物通过SEQ ID NO.2所示的第二测序通用引物+1-7bp的随机碱基+SEQID NO.4所示的肠道微生物16S rRNA V3-V4区特异性反向引物组合得到，为5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-TCAGA-GGACTACNVGGGTWTCTAATC-3'，例如可以是5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-(N)1-7-CCTACGGRRBGCASCAGKVRVGAAT-3'.
根据本发明，所述P5端接头序列的核酸序列如SEQ ID NO.5所示，所述P7端接头序列的核酸序列如SEQ ID NO.6所示，具体核酸序列如下：
P5端接头序列(SEQ ID NO.5)：AATGATACGGCGACCACCGAGATCTACAC；
P7端接头序列(SEQ ID NO.6)：CAAGCAGAAGACGGCATACGAGAT.
根据本发明，所述第一标签序列的核酸序列如SEQ ID NO.7-16所示，所述核酸序列如下：
序列 接头名称
SEQ ID NO.7 CTCTCTAT S502
SEQ ID NO.8 TATCCTCT S503
SEQ ID NO.9 GTAAGGAG S505
SEQ ID NO.10 ACTGCATA S506
SEQ ID NO.11 AAGGAGTA S507
SEQ ID NO.12 CTAAGCCT S508
SEQ ID NO.13 CGTCTAAT S510
SEQ ID NO.14 TCTCTCCG S511
SEQ ID NO.15 TCGACTAG S513
SEQ ID NO.16 TTCTAGCT S515
根据本发明，所述第二标签序列的核酸序列如SEQ ID NO.17-26所示，所述核酸序列如下：
根据本发明，所述第一测序通用引物的部分序列如SEQ ID NO.27所示，所述核酸序列如下：
第一测序通用引物的部分序列(SEQ ID NO.27)：TCGTCGGCAGCGTC.
根据本发明，所述第二测序通用引物的部分序列如SEQ ID NO.28所示，所述核酸序列如下：
第二测序通用引物的部分序列(SEQ ID NO.28)：GTCTCGTGGGCTCGG.
本发明中，所述第二正向引物通过SEQ ID NO.5所示的P5接头+第一标签序列+SEQID NO.27所示的第一测序通用引物的部分序列组合得到，例如可以是：AATGATACGGCGACCACCGAGATCTACAC-CTCTCTAT-TCGTCGGCAGCGTC；所述第二反向引物通过SEQID NO.6所示的P7接头+第二标签序列+SEQ ID NO.28所示的第二测序通用引物的部分序列组合得到，例如可以是CAAGCAGAAGACGGCATACGAGAT-TCGCCTTA-GTCTCGTGGGCTCGG.
第二方面，本发明提供一种肠道微生物的分析方法，采用如第一方面所述的引物组合物，包括如下步骤：
(1)样本基因组DNA的提取；
(2)扩增子文库的构建：进行两轮PCR扩增，第一轮PCR扩增采用第一组引物组，第二轮PCR扩增采用第二组引物组，构建高通量测序文库；
(3)数据分析。
根据本发明，所述样本来源于粪便、土壤、尿液或唾液中的任意一种或至少两种的组合，优选为来源于粪便。
本发明，步骤(1)所述样本基因组DNA的提取的方法为常规的DNA提取方法，本领域技术人员可以根据需要进行选择，在此不作特殊限定。
根据本发明，步骤(2)所述第一轮PCR的条件为：93-96℃预变性1-5min；93-96℃预变性20-40s，50-55℃20-40s，70-75℃20-40s，共进行25-35个循环；70-75℃延伸3-6min。
优选地，步骤(2)所述第一轮PCR的条件为：95℃预变性3min；95℃预变性30s，53℃30s，72℃30s，共进行28个循环；72℃延伸5min。
根据本发明，步骤(2)所述第二轮PCR的条件为：93-96℃预变性1-5min；93-96℃预变性20-40s，54-58℃20-40s，70-75℃20-40s，共进行6-12个循环；70-75℃延伸3-6min。
优选地，步骤(2)所述第二轮PCR的条件为：95℃预变性3min；95℃预变性30s，55.5℃延伸30s，72℃30s，8个循环；72℃5min，1个循环。
根据本发明，步骤(2)所述第一轮PCR之后还包括纯化的步骤，优选采用磁珠进行纯化。
根据本发明，所述磁珠与扩增产物的体积比为(0.4-0.7):1，例如可以是0.4:1、0.5:1、0.6:1或0.7:1，优选为0.5:1。
根据本发明，步骤(2)所述第二轮PCR之后还包括纯化的步骤，优选采用磁珠进行纯化。
根据本发明，所述磁珠与扩增产物的体积比为(0.7-1.1):1，例如可以是0.7:1、0.8:1、0.9:1、1:1或1.1:1，优选为0.9:1。
根据本发明，步骤(3)所述数据分析为采用Illumina Miseq测序仪进行测序分析。
根据本发明，所述测序分析具体包括如下步骤：通过Illumina Miseq测序仪测序得到按照标签区分的测序读长，利用读长的重叠关系组装得到高可变区V3-V4的全长序列，再对全长序列进行分类分析，得到所述微生物群体的分类。
根据本发明，所述全长序列进行分类分析具体包括计算全长序列差异度，根据序列差异度执行操作分类学单元OTU的分类，将全长序列分配到OTU中，再将每一个OTU分类中的全长序列比对到16S rRNA的V6数据库中，将比对结果根据众数原则对OTU进行物种注释。
第三方面，本发明提供一种如第一方面所述的试剂盒或如第三方面所述的分析方法用于分析分类肠道微生物。
第四方面，本发明提供一种如第三方面所述的分析方法用于免疫组学分析。
第五方面，本发明提供一种如第一方面所述的试剂盒用于制备癌症的诊断试剂。
优选地，所述癌症选自但不限于肠癌、腺癌或膀胱癌中的任意一种或至少两种的组合。
与现有技术相比，本发明具有如下有益效果：
(1)本发明中，通过采用两步法扩增获得测序文库，所述试剂盒中的引物组合物可以更大程度地混合更多的样本，提高了扩增片段覆盖率、检测特异性及灵敏度，降低了PCR扩增偏好性，结合数据分析，可获得更丰富的微生物组信息，增加了文库序列的复杂度，降低了phix的比例，保证更多测序质量；
(2)本发明提供的数据处理分析方法，实现更多微生物物种在分类学属、种上的鉴别及多样性分析，并得到相关微生物群体的相对丰度值；
(3)本发明方法对肠道微生物进行测序及数据分析，实现了更多物种在分类学属甚至种上的鉴别，即获得更多肠道微生物生物学信息，可为免疫组学研究以及临床用药上提供一定指导意义；
(4)本发明方法可用于检测癌症，癌症的检出率可达85％以上，尤其针对肠癌的检出率可达96.92％，针对腺瘤的检出率可达85.25％。
附图说明
图1为本发明实施例中样本A数据分析后相关性分析的结果图；
图2为本发明实施例中样本B数据分析后相关性分析的结果图；
图3为本发明实施例中样本C数据分析后相关性分析的结果图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
材料与试剂：
缓冲液A：0.5M Tris-Cl缓冲液(pH8.0)，1.0M NaCl；
裂解液B：1M Tris-Cl缓冲液(pH 8.0)，0.2M EDTA，10％SDS，12％Tween-40；
悬浮液C：30mg/ml RNAaseA；
去除剂D：1.0M Al2(SO4)3，250g/L CTAB，1.5M NaCl；
结合液E：0.5M Tris-Cl缓冲液(pH 8.0)，5M异硫腈酸胍；
漂洗液F：5M Tris-Cl缓冲液(pH 8.0)，25％(体积比)无水乙醇；
洗脱液G：无菌去离子水。
实施例1 检测试剂盒的制备
(1)引物组合物的制备
第一组引物：
所述第一测序通用引物的核酸序列如SEQ ID NO.1所示，所述第二测序通用引物的核酸序列如SEQ ID NO.2所示，所述核酸序列如下：
第一测序通用引物(SEQ ID NO.1)：TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG；
第二测序通用引物(SEQ ID NO.2)：GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG；
所述随机碱基序列为1-7个随机碱基，如(N)1-7，优选为3-5个随机碱基，如(N)3-5，其中，N表示为A\T\G\C中的任意一个碱基；
所述肠道微生物16S rRNA V3-V4区特异性正向引物的核酸序列如SEQ ID NO.3所示，所述肠道微生物16S rRNA V3-V4区特异性反向引物的核酸序列如SEQ ID NO.4所示，具体核酸序列如下：
肠道微生物16S rRNA V3-V4区特异性正向引物(SEQ ID NO.3)：CCTACGGRRBGCASCAGKVRVGAAT；
肠道微生物16S rRNA V3-V4区特异性反向引物(SEQ ID NO.4)：GGACTACNVGGGTWTCTAATC.
第二组引物：
所述P5端接头序列的核酸序列如SEQ ID NO.5所示，所述P7端接头序列的核酸序列如SEQ ID NO.6所示，具体核酸序列如下：
P5端接头序列(SEQ ID NO.5)：AATGATACGGCGACCACCGAGATCTACAC；
P7端接头序列(SEQ ID NO.6)：CAAGCAGAAGACGGCATACGAGAT；
所述第一标签序列的核酸序列如SEQ ID NO.7-16所示，所述核酸序列如下：
序列 接头名称
SEQ ID NO.7 CTCTCTAT S502
SEQ ID NO.8 TATCCTCT S503
SEQ ID NO.9 GTAAGGAG S505
SEQ ID NO.10 ACTGCATA S506
SEQ ID NO.11 AAGGAGTA S507
SEQ ID NO.12 CTAAGCCT S508
SEQ ID NO.13 CGTCTAAT S510
SEQ ID NO.14 TCTCTCCG S511
SEQ ID NO.15 TCGACTAG S513
SEQ ID NO.16 TTCTAGCT S515
所述第二标签序列的核酸序列如SEQ ID NO.17-26所示，所述核酸序列如下：
序列 接头名称
SEQ ID NO.17 TCGCCTTA N701
SEQ ID NO.18 CTAGTACG N702
SEQ ID NO.19 TTCTGCCT N703
SEQ ID NO.20 GCTCAGGA N704
SEQ ID NO.21 AGGAGTCC N705
SEQ ID NO.22 CATGCCTA N706
SEQ ID NO.23 GTAGAGAG N707
SEQ ID NO.24 CAGCCTCG N710
SEQ ID NO.25 TGCCTCTT N711
SEQ ID NO.26 TCCTCTAC N712
根据本发明，所述第一测序通用引物的部分序列如SEQ ID NO.27所示，所述核酸序列如下：
第一测序通用引物的部分序列(SEQ ID NO.27)：TCGTCGGCAGCGTC；
所述第二测序通用引物的部分序列如SEQ ID NO.28所示，所述核酸序列如下：
第二测序通用引物的部分序列(SEQ ID NO.28)：GTCTCGTGGGCTCGG.
(2)试剂盒的组装
将所述引物组合物和检测探针与试剂盒相关的试剂进行组装，制备成所述检测试剂盒。
实施例2 样本DNA的提取
选取3个粪便样本：A、B、C，平均分成2份，标记为：A1/A2，B1/B2，C1/C2，A1、B1、C1三个样本用本专利说明的方法处理并获得数据，A2、B2、C2三个样本直接送去金唯智测序，具体的DNA提取方法包括如下步骤：
(1)取700～800μl缓冲液A和250mg玻璃珠至2ml离心管中；
(2)在上述2ml离心管中加入250mg样本，涡旋混匀10s；
(3)向样本中加入70μl裂解液B，涡旋振荡7min混匀样本，12,000rpm(～13,400×g)，离心60s。转移上清液(500μl)至新的2ml离心管；
(4)加入250μl悬浮液C，涡旋振荡10s，室温放置2min，12,000rpm(～13,400×g)，离心60s，沉淀样本颗粒；
(5)转移上清至新的2ml离心管，加入300μl去除剂D混匀，4℃放置2～5min；
(6)12,000rpm(～13,400×g)，离心60s。转移上清液至新的2ml离心管，加入1000μl结合液E颠倒混匀；
(7)取上一步所得溶液700μl加入到一个吸附柱中(吸附柱放入收集管中)，12,000rpm(～13,400×g)，离心30s，倒掉废液，将吸附柱放入收集管中；
(8)向吸附柱中，加入700μl漂洗液F，12,000rpm(～13,400×g)，离心30s，倒掉废液，将吸附柱放入收集管中；
(9)向吸附柱中，加入500μl 75％无水乙醇，12,000rpm(～13,400×g)，离心30s，倒掉废液，将吸附柱放入收集管中；
(10)将吸附柱放回空收集管中，12,000rpm(～13,400×g)，离心60s，尽量除去漂洗液，以免漂洗液中残留乙醇抑制下游反应；
(11)取出吸附柱，放入一个干净的离心管中，在吸附膜的中间部位加50～100μl洗脱液G(洗脱液可事先在60～70℃水浴中预热)，室温放置2min，12,000rpm(～13,400×g)，离心60s，将得到的溶液重新加入离心吸附柱中，室温放置2min，12,000rpm(～13,400×g)，离心1min；
(12)使用分光光度计或Qubit对样品进行精确定量，DNA可以存放在0～10℃，如果要长时间存放，可以放置在-20℃。
将提取后的DNA进行纯度测试，OD260/OD280的比值在1.7～1.9之间，可以直接用于后续的PCR建库。
实施例3 文库构建
(1)在冰上配制以下反应液，除DNA提取液和H2O以外，依照如下组份先按反应数+α的量配制扩增预混Mix，取配制好的预混Mix分装到PCR反应管中，具体引物序列如下：
上游引物(SEQ ID NO.29)：5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-ATGCA-CCTACGGA GTGCAGCAGTAGCGAAT-3'；
下游引物(SEQ ID NO.30)：5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-TCAGA-GGACTAC ACGGGTATCTAATC-3'；
反应体系如下：
设置一组以水为样本的阴性对照。
反应条件如下：
(2)第一轮PCR产物的纯化
使用Agencourt AMPure XP(Beckman Coulter公司)对第一轮PCR扩增产物进行纯化，具体步骤如下：
1)4℃冰箱取出，室温放置25min，将Agencourt AMPure XP磁珠充分打散，混匀；
2)扩增产物与Agencourt AMPure XP磁珠的结合，所述扩增产物与AgencourtAMPure XP磁珠的体积比为2:1；
3)充分混匀，室温孵育15min，瞬离，磁力架上静置5min，去除上清液；
4)第一次清洗：反应管放置在磁力架上不动，加入200μl新配的80％的乙醇室温静置(>30s)，去除80％的乙醇；
5)第二次清洗：重复一次上述4)第一次清洗操作，轻轻离心，去除残留80％的乙醇；
6)磁珠的干燥：室温干燥1min；
7)洗脱：加入30μl的洗脱液，轻轻地上下吸打，充分混匀，室温孵育15min，瞬离，磁力架上静置5min，回收上清液；
(3)片段大小确定
琼脂糖电泳确定扩增是否准确，扩增产物长度约550bp。
(4)第二轮PCR
将步骤(2)中纯化得到的PCR产物为模板，进行PCR导入测序接头，具体步骤如下：
1)在冰上配制以下反应液，除第一轮PCR纯化产物以外，依照如下组份先按反应数+α的量配制扩增预混Mix，取配制好的预混Mix分装到PCR反应管中，具体引物序列如下：
上游引物(SEQ ID NO.31)：AATGATACGGCGACCACCGAGATCTACAC-CTCTCTAT-TCGTCGGCAGCGTC；
下游引物(SEQ ID NO.32)：CAAGCAGAAGACGGCATACGAGAT-TCGCCTTA-GTCTCGTGGGCTCGG；
反应体系如下：
试剂 使用量
2×PCR混合液(Mg2+,dNTPs) 12.5μl
上游引物(1μM) 5μl
下游引物(1μM) 5μl
DNA(第一轮PCR纯化产物) 2.5μl
Total 25μl
设置一组以水为样本的阴性对照。
反应条件如下：
(5)第二轮PCR产物的纯化
具体步骤同第一轮PCR产物纯化，即步骤(2)，其中PCR产物与Agencourt AMPureXP磁珠的体积比为10:9；
(6)片段大小及浓度测定
琼脂糖电泳确定扩增是否准确，扩增产物长度约600bp；
(7)使用Qubit dsDNA HS Assay Kit(Thermo Fisher)等测定PCR产物的浓度；
(8)测序建好的文库用illumina公司的Miseq测序进行测序，测序侧率为PE300测序。
实施例4 数据分析
(1)对测序得到的原始数据进行质控：通过Barcode和引物等信息筛选出高质量的数据，筛选标准如下：含有完整的序列；具有完整前引物序列；去除Barcode、引物、接头，序列长度必须大于400bp，不含有模糊碱基；不能存在超过20个的单碱基重复，将得到的序列根据信息分别与各个样品相对应；
(2)通过数据库比对工具进行比对，将序列聚类生成OUT；
(3)将得到的OTU表，按照A1/A2，B1/B2，C1/C2进行相关性分析，结果如图1-3所示。
从图1可以看出，样本A，分成两份后，分别用本专利方法和金唯智公司方法建库测序，得到的数据分析后进行相关性分析发现，两种方法得到的OUT相关性达到0.9857；从图2可以看出，样本B，分成两份后，分别用本专利方法和金唯智公司方法建库测序，得到的数据分析后进行相关性分析发现，两种方法得到的OUT相关性达到0.9774；从图3可以看出，样本C，分成两份后，分别用本专利方法和金唯智公司方法建库测序，得到的数据分析后进行相关性分析发现，两种方法得到的OUT相关性达到0.9809。
实施例5 分析验证
选取65例肠癌患者粪便样本，61例腺瘤患者粪便样本和65例健康患者粪便样本，用本专利所述的方法提取建库测序获得DNA序列信息，根据序列信息预测肠癌风险，具体试验方法同实施例2-实施例4，结果如表1所示。
表1
肠癌 腺瘤 健康人
实际数目 65 61 65
肠癌检出数目 63 8 0
腺瘤检出数目 2 52 7
健康检出数目 0 1 58
准确率 96.92％ 85.25％ 89.23％
从表1可以看出，本发明检测中，肠癌的检出率为96.92％，腺瘤的检出率为85.25％，健康人的准确率为89.23％，可见，本发明方法用于检测癌症的准确率可达96％以上。
综上所述，本发明中，通过采用两步法扩增获得测序文库，所述引物组合物可以更大程度地混合更多的样本，提高了扩增片段覆盖率、检测特异性及灵敏度，结合数据分析，可获得更丰富的微生物组信息，增加了文库序列的复杂度，降低了phix的比例，保证更多测序质量，降低了PCR扩增偏好性；本发明提供的数据处理分析方法，实现更多微生物物种在分类学属、种上的鉴别及多样性分析，并得到相关微生物群体的相对丰度值。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
序列表
<110> 苏州普瑞森基因科技有限公司
<120> 一种用于分析肠道微生物的试剂盒及其应用
<130> 2017
<141> 2017-12-29
<160> 32
<170> SIPOSequenceListing 1.0
<210> 1
<211> 33
<212> DNA
<213> 人工合成序列()
<400> 1
tcgtcggcag cgtcagatgt gtataagaga cag 33
<210> 2
<211> 34
<212> DNA
<213> 人工合成序列()
<400> 2
gtctcgtggg ctcggagatg tgtataagag acag 34
<210> 3
<211> 25
<212> DNA
<213> 人工合成序列()
<400> 3
cctacggrrb gcascagkvr vgaat 25
<210> 4
<211> 21
<212> DNA
<213> 人工合成序列()
<220>
<221> misc_feature
<222> (8)..(8)
<223> n is a, c, g, or t
<400> 4
ggactacnvg ggtwtctaat c 21
<210> 5
<211> 29
<212> DNA
<213> 人工合成序列()
<400> 5
aatgatacgg cgaccaccga gatctacac 29
<210> 6
<211> 24
<212> DNA
<213> 人工合成序列()
<400> 6
caagcagaag acggcatacg agat 24
<210> 7
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 7
ctctctat 8
<210> 8
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 8
tatcctct 8
<210> 9
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 9
gtaaggag 8
<210> 10
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 10
actgcata 8
<210> 11
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 11
aaggagta 8
<210> 12
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 12
ctaagcct 8
<210> 13
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 13
cgtctaat 8
<210> 14
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 14
tctctccg 8
<210> 15
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 15
tcgactag 8
<210> 16
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 16
ttctagct 8
<210> 17
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 17
tcgcctta 8
<210> 18
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 18
ctagtacg 8
<210> 19
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 19
ttctgcct 8
<210> 20
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 20
gctcagga 8
<210> 21
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 21
aggagtcc 8
<210> 22
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 22
catgccta 8
<210> 23
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 23
gtagagag 8
<210> 24
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 24
cagcctcg 8
<210> 25
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 25
tgcctctt 8
<210> 26
<211> 8
<212> DNA
<213> 人工合成序列()
<400> 26
tcctctac 8
<210> 27
<211> 14
<212> DNA
<213> 人工合成序列()
<400> 27
tcgtcggcag cgtc 14
<210> 28
<211> 15
<212> DNA
<213> 人工合成序列()
<400> 28
gtctcgtggg ctcgg 15
<210> 29
<211> 63
<212> DNA
<213> 人工合成序列()
<400> 29
tcgtcggcag cgtcagatgt gtataagaga cagatgcacc tacggagtgc agcagtagcg 60
aat 63
<210> 30
<211> 60
<212> DNA
<213> 人工合成序列()
<400> 30
gtctcgtggg ctcggagatg tgtataagag acagtcagag gactacacgg gtatctaatc 60
<210> 31
<211> 51
<212> DNA
<213> 人工合成序列()
<400> 31
aatgatacgg cgaccaccga gatctacacc tctctattcg tcggcagcgt c 51
<210> 32
<211> 47
<212> DNA
<213> 人工合成序列()
<400> 32
caagcagaag acggcatacg agattcgcct tagtctcgtg ggctcgg 47

Claims (10)

1.一种用于分析肠道微生物的试剂盒，其包含用于构建肠道微生物文库的引物组合物，所述引物组合物包括两组引物，具体如下：
第一组引物，所述第一组引物包含第一正向引物和第一反向引物，其中，所述第一正向引物从5’到3’端包括依次相连的三个部分：第一测序通用引物、1-7bp的随机碱基序列和肠道微生物16S rRNA V3-V4区特异性正向引物；所述第一反向引物从5’到3’端包括依次相连的三个部分：第二测序通用引物、1-7bp的随机碱基序列和肠道微生物16S rRNA V3-V4区特异性反向引物；
第二组引物，所述第二组引物包含第二正向引物和第二反向引物，其中，所述第二正向引物从5’到3’端包括依次相连的三个部分：P5端接头序列、第一标签序列和第一测序通用引物的部分序列；所述第二反向引物从5’到3’端包括依次相连的三个部分：P7端接头序列、第二标签序列和第二测序通用引物的部分序列。
2.根据权利要求1所述的试剂盒，其特征在于，所述第一测序通用引物和所述第二测序通用引物独立地选自任意一种二代测序平台的测序通用引物；
优选地，所述第一测序通用引物的核酸序列如SEQ ID NO.1所示，所述第二测序通用引物的核酸序列如SEQ ID NO.2所示；
优选地，所述随机碱基的碱基个数为3-5个；
优选地，所述肠道微生物16S rRNA V3-V4区特异性正向引物的核酸序列如SEQ IDNO.3所示，所述肠道微生物16S rRNA V3-V4区特异性反向引物的核酸序列如SEQ ID NO.4所示。
3.根据权利要求1或2所述的试剂盒，其特征在于，所述P5端接头序列的核酸序列如SEQID NO.5所示，所述P7端接头序列的核酸序列如SEQ ID NO.6所示；
优选地，所述第一标签序列如SEQ ID NO.7-16所示，所述第二标签序列如SEQ IDNO.17-26所示；
优选地，所述第一测序通用引物的部分序列如SEQ ID NO.27所示，所述第二测序通用引物的部分序列如SEQ ID NO.28所示。
4.一种肠道微生物的分析方法，其特征在于，采用如权利要求1-3中任一项所述的试剂盒，包括如下步骤：
(1)样本基因组DNA的提取；
(2)扩增子文库的构建：进行两轮PCR扩增，第一轮PCR扩增采用第一组引物组，第二轮PCR扩增采用第二组引物组，构建高通量测序文库；
(3)数据分析。
5.根据权利要求4所述的方法，其特征在于，所述样本来源于粪便、土壤、尿液或唾液中的任意一种或至少两种的组合，优选为来源于粪便。
6.根据权利要求4或5所述的方法，其特征在于，步骤(2)所述第一轮PCR的条件为：93-96℃预变性1-5min；93-96℃预变性20-40s，50-55℃20-40s，70-75℃20-40s，共进行25-35个循环；70-75℃延伸3-6min；
优选地，步骤(2)所述第一轮PCR的条件为：95℃预变性3min；95℃预变性30s，53℃30s，72℃30s，共进行28个循环；72℃延伸5min；
优选地，步骤(2)所述第二轮PCR的条件为：93-96℃预变性1-5min；93-96℃预变性20-40s，54-58℃20-40s，70-75℃20-40s，共进行6-12个循环；70-75℃延伸3-6min；
优选地，步骤(2)所述第二轮PCR的条件为：95℃预变性3min；95℃预变性30s，55.5℃延伸30s，72℃30s，8个循环；72℃5min，1个循环。
7.根据权利要求4-6中任一项所述的方法，其特征在于，步骤(2)所述第一轮PCR之后还包括纯化的步骤，优选采用磁珠进行纯化；
优选地，所述磁珠与扩增产物的体积比为(0.4-0.7):1，优选为0.5:1；
优选地，步骤(2)所述第二轮PCR之后还包括纯化的步骤，优选采用磁珠进行纯化；
优选地，所述磁珠与扩增产物的体积比为(0.7-1.1):1，优选为0.9:1。
8.根据权利要求4-7中任一项所述的方法，其特征在于，步骤(3)所述数据分析为采用Illumina Miseq测序仪进行测序分析；
优选地，所述测序分析具体包括如下步骤：通过Illumina Miseq测序仪测序得到按照标签区分的测序读长，利用读长的重叠关系组装得到高可变区V3-V4的全长序列，再对全长序列进行分类分析，得到所述微生物群体的分类；
优选地，所述全长序列进行分类分析具体包括计算全长序列差异度，根据序列差异度执行操作分类学单元OTU的分类，将全长序列分配到OTU中，再将每一个OTU分类中的全长序列比对到16S rRNA的V6数据库中，将比对结果根据众数原则对OTU进行物种注释。
9.一种如权利要求1-3中任一项所述的试剂盒或如权利要求4-8中任一项所述的分析方法用于分析分类肠道微生物；
优选地，一种如权利要求4-10中任一项所述的分析方法用于免疫组学分析。
10.一种如权利要求1-3中任一项所述的试剂盒用于制备癌症的诊断试剂；
优选地，所述癌症为肠癌、腺癌或膀胱癌中的任意一种或至少两种的组合。
CN201711478123.6A 2017-12-29 2017-12-29 一种用于分析肠道微生物的试剂盒及其应用 Pending CN107904298A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711478123.6A CN107904298A (zh) 2017-12-29 2017-12-29 一种用于分析肠道微生物的试剂盒及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711478123.6A CN107904298A (zh) 2017-12-29 2017-12-29 一种用于分析肠道微生物的试剂盒及其应用

Publications (1)

Publication Number Publication Date
CN107904298A true CN107904298A (zh) 2018-04-13

Family

ID=61872026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711478123.6A Pending CN107904298A (zh) 2017-12-29 2017-12-29 一种用于分析肠道微生物的试剂盒及其应用

Country Status (1)

Country Link
CN (1) CN107904298A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109576386A (zh) * 2019-01-29 2019-04-05 苏州普瑞森基因科技有限公司 一种鉴定肠道微生态状态的引物组合物及其应用
CN109988819A (zh) * 2019-04-19 2019-07-09 南京诺唯赞生物科技有限公司 一种特异性接头分子标签及其在预测可变剪切中的应用
CN110029181A (zh) * 2019-04-19 2019-07-19 中山大学附属第三医院 一种***性红斑狼疮的病情评估***及其试剂盒
CN110176305A (zh) * 2019-05-27 2019-08-27 天益健康科学研究院（镇江）有限公司 一种利用高通量基因测序评估肠道菌群健康的方法
CN112029862A (zh) * 2020-09-09 2020-12-04 上海宝藤生物医药科技股份有限公司 一种基于无核酸提取的文库构建检测结直肠癌标志物的试剂盒、方法及其应用
CN113088561A (zh) * 2019-12-23 2021-07-09 深圳华大智造科技有限公司 构建侧翼已知序列的测序文库的引物组和方法
CN115058490A (zh) * 2022-06-28 2022-09-16 广州市金圻睿生物科技有限责任公司 一种用于构建微生物靶向测序文库的引物组合及其应用

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104153004A (zh) * 2014-08-11 2014-11-19 上海美吉生物医药科技有限公司 一种用于扩增子测序的建库方法
EP2876167A1 (en) * 2013-11-21 2015-05-27 Institut Gustave Roussy Microbiota composition, as a marker of responsiveness to chemotherapy, and use of microbial modulators (pre-,pro- or synbiotics) for improving the efficacy of a cancer treatment
CN106192023A (zh) * 2016-08-08 2016-12-07 中国科学院北京基因组研究所 一种基于多维Index的多重测序文库构建方法
CN106497926A (zh) * 2016-11-03 2017-03-15 承启医学（深圳）科技有限公司 一种用于构建微生物细菌16s　rDNA可变区测序文库的扩增子引物及构建方法
CN106555226A (zh) * 2016-04-14 2017-04-05 北京京诺玛特科技有限公司 一种构建高通量测序文库的方法和试剂盒
CN107292123A (zh) * 2016-03-31 2017-10-24 苏州普瑞森基因科技有限公司 一种基于高通量测序的微生物群落组成的方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2876167A1 (en) * 2013-11-21 2015-05-27 Institut Gustave Roussy Microbiota composition, as a marker of responsiveness to chemotherapy, and use of microbial modulators (pre-,pro- or synbiotics) for improving the efficacy of a cancer treatment
CN104153004A (zh) * 2014-08-11 2014-11-19 上海美吉生物医药科技有限公司 一种用于扩增子测序的建库方法
CN107292123A (zh) * 2016-03-31 2017-10-24 苏州普瑞森基因科技有限公司 一种基于高通量测序的微生物群落组成的方法和装置
CN106555226A (zh) * 2016-04-14 2017-04-05 北京京诺玛特科技有限公司 一种构建高通量测序文库的方法和试剂盒
CN106192023A (zh) * 2016-08-08 2016-12-07 中国科学院北京基因组研究所 一种基于多维Index的多重测序文库构建方法
CN106497926A (zh) * 2016-11-03 2017-03-15 承启医学（深圳）科技有限公司 一种用于构建微生物细菌16s　rDNA可变区测序文库的扩增子引物及构建方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
付玉荣等: "《临床微生物学检验实验》", 31 August 2013, 武汉：华中科技大学出版社 *
樊龙江等: "《生物信息学》", 30 September 2017, 杭州：浙江大学出版社 *
苏为华等: "《综合评价技术的扩展与集成问题研究》", 30 December 2007, 北京：中国统计出版社 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109576386A (zh) * 2019-01-29 2019-04-05 苏州普瑞森基因科技有限公司 一种鉴定肠道微生态状态的引物组合物及其应用
CN109988819A (zh) * 2019-04-19 2019-07-09 南京诺唯赞生物科技有限公司 一种特异性接头分子标签及其在预测可变剪切中的应用
CN110029181A (zh) * 2019-04-19 2019-07-19 中山大学附属第三医院 一种***性红斑狼疮的病情评估***及其试剂盒
CN110176305A (zh) * 2019-05-27 2019-08-27 天益健康科学研究院（镇江）有限公司 一种利用高通量基因测序评估肠道菌群健康的方法
CN113088561A (zh) * 2019-12-23 2021-07-09 深圳华大智造科技有限公司 构建侧翼已知序列的测序文库的引物组和方法
CN113088561B (zh) * 2019-12-23 2024-05-14 深圳华大智造科技股份有限公司 构建侧翼已知序列的测序文库的引物组和方法
CN112029862A (zh) * 2020-09-09 2020-12-04 上海宝藤生物医药科技股份有限公司 一种基于无核酸提取的文库构建检测结直肠癌标志物的试剂盒、方法及其应用
CN115058490A (zh) * 2022-06-28 2022-09-16 广州市金圻睿生物科技有限责任公司 一种用于构建微生物靶向测序文库的引物组合及其应用

Similar Documents

Publication Publication Date Title
CN107904298A (zh) 一种用于分析肠道微生物的试剂盒及其应用
CN107868837B (zh) 一种用于分析肠道微生物的引物组合物及其应用
Lesker et al. An integrated metagenome catalog reveals new insights into the murine gut microbiome
McCartney Application of molecular biological methods for studying probiotics and the gut flora
Boughner et al. Microbial ecology: where are we now?
CN109273053B (zh) 一种高通量测序的微生物数据处理方法
Brennan et al. Characterization of environmentally persistent Escherichia coli isolates leached from an Irish soil
CN107937582A (zh) 一种用于分析肠道微生物的引物组及其应用
JP2005526952A (ja) 個体群における多型マーカーを同定するための方法
CN111315884B (zh) 测序文库的归一化
IE840683L (en) Method for characterizing an unknown organism species by¹identifying evolutionary conserved sequences in genetic¹material.
CN110904250B (zh) 用于检测多种菌的多重荧光定量pcr引物、试剂盒及检测方法
CN109971871A (zh) 一种筛选和/或鉴定乳杆菌的方法及其应用
US20140011686A1 (en) Selection of single nucleic acids based on optical signature
Jarocki et al. Molecular routes to specific identification of the Lactobacillus casei group at the species, subspecies and strain level
CN110452974A (zh) 一种检测细菌16S rDNA全长的建库测序方法
Kashyap et al. Ribotyping: A tool for molecular taxonomy
Zhang et al. Rapid strain-specific identification of two Lactobacillus rhamnosus strains using PCR based on gene family analysis
Wang et al. Substrate utilization and competitive interactions among soil bacteria vary with life-history strategies
CN104560982B (zh) 不同种属微生物间种类和丰度比较的人工外源性参照分子
Lan et al. Massively parallel single-cell sequencing of diverse microbial populations
Thies Molecular methods for studying soil ecology
CN107674920A (zh) 嵌合体多重pcr引物组合物和检测方法
CN109652573B (zh) 用于鼠伤寒沙门菌或其单相菌变种分型检测的vntr位点、检测引物组及检测分析方法
Huang et al. Development of a high-resolution single-nucleotide polymorphism strain-typing assay using whole genome-based analyses for the Lactobacillus acidophilus probiotic strain

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180413