CN110804655A

CN110804655A - 一种宏基因组绝对定量的方法

Info

Publication number: CN110804655A
Application number: CN201911116970.7A
Authority: CN
Inventors: 胡悦; 许冬瑾; 蒙洁妙
Original assignee: Kang Meihua Big Gene Technology Co Ltd
Current assignee: Kang Meihua Big Gene Technology Co Ltd
Priority date: 2019-11-15
Filing date: 2019-11-15
Publication date: 2020-02-18
Anticipated expiration: 2039-11-15
Also published as: CN110804655B

Abstract

本发明涉及宏基因组绝对定量的方法。方法包括：在样品中添加第一定量DNA片段；对所述样品进行基因组DNA的提取；在所述提取的基因组DNA中添加第二定量DNA片段，得到总基因组；对所述总基因组进行宏基因组文库构建和测序；利用得到的测序结果的数据计算单位质量的所述样品中的细菌数量。通过本发明的方法获得细菌的绝对定量，不用采用高成本的流式细胞计数，也不用对特定的DNA进行定量PCR，简化了过程并且降低了成本。

Description

一种宏基因组绝对定量的方法

技术领域

本发明涉及宏基因组绝对定量的方法。

背景技术

目前16S和宏基因组测序研究中，菌群的数量和结构主要通过相对丰度来表示，然而基于相对值的菌群分析方法局限性较大。因此，实现对单个菌株或者目标基因的绝对定量，对于研究菌群功能具有十分重要的意义。在绝对定量的过程中，不同样品(如粪便或土壤样本)之间、甚至同一样品的多个重复之间DNA提取得率的差异同样也会显著影响目标菌株的绝对定量结果。

发明内容

本发明通过分别在粪便样本基因组DNA(genomic DNA,gDNA)提取前和提取后加入两段已知长度、序列和准确质量的标准品DNA片段(sDNA1和sDNA2)；通过一系列计算方法，结合提取gDNA所用的粪便质量，最终得出单位质量粪便中细菌的分子数。本发明的方法无需进行qPCR定量或者流式细菌计数，操作简单易行且实验成本低。

本发明提供了一种细菌定量的方法，包括：在样品中添加第一定量DNA片段；对所述样品进行基因组DNA的提取；在所述提取的基因组DNA中添加第二定量DNA片段，得到总基因组；对所述总基因组进行宏基因组文库构建和测序；利用得到的测序结果的数据计算单位质量的所述样品中的细菌数量。

在上述方法中，其中，所述样品为粪便样品。

在上述方法中，其中，所述第一定量DNA片段的序列为序列1。

在上述方法中，其中，所述第二定量DNA片段的序列为序列2。

在上述方法中，其中，所述样品中的所述第一定量DNA片段的分子数为：测序所得的所述第一定量DNA片段的拷贝数与单位质量的所述样品中加入的所述第二定量DNA片段的分子数的乘积再除以测序所得的所述第二定量DNA片段的分子数。

在上述方法中，其中，所述样品中的基因组DNA的得率为：单位质量的所述样品中的所述总基因组中的所述第一定量DNA片段的分子数除以加入的所述第一定量DNA片段的分子数，之后再乘以所述样品中的基因组DNA的总量。

在上述方法中，其中，单位质量的所述样品的细菌的数量为：单位质量的所述样品中加入的所述第二定量DNA片段的分子数、测序所得的细菌的总拷贝数和所述样品的基因组DNA的总量的乘积，再除以测序所得的所述第二定量DNA片段的拷贝数、所述样品中的基因组DNA的得率和所述样品的质量的乘积。

本发明不用针对粪便样本进行流式细菌计数实验，极大节省了该方面的实验成本及时间。通过在提取后的粪便样本gDNA后加入sDNA2(已知序列、相对分子质量以及加入量，可得出其准确分子数)，不用针对该片段进行专门的qPCR实验进行定量，同样极大节省了该方面的实验成本及时间。另外，sDNA1和sDNA2片段均来自植物基因组，与细菌的参考基因组序列无同源性；此外，本发明的方法中采用的是宏基因组测序，获得的菌群数据信息相比16S测序更全面准确。

附图说明

图1示出了本发明的方法的流程的示意图。

图2示出了本发明的结果的计算方法的示意图。

具体实施方式

下面的实施例可以使本领域技术人员更全面地理解本发明，但不以任何方式限制本发明。

为了提高例如粪便样本中的目标菌株的绝对定量结果的准确性，一种方法是利用流式细胞仪检测粪便样本中细菌细胞的个数，之后结合16S测序获得的肠道菌群丰度数据，实现对菌群的绝对定量；另一种方法则是通过向样品DNA中添加一定量外标序列，进行16S扩增子文库构建、测序，再根据外标序列的16S扩增子读数及其绝对拷贝数绘制标准曲线，计算出样品中操作分类单元(OTU)代表序列对应物种16S rRNA基因绝对拷贝数。

本发明通过分别在粪便样本gDNA提取前和提取后加入两段已知长度、序列和准确质量的标准品DNA片段(sDNA1和sDNA2)；其次，通过计算后加入sDNA2片段的准确分子数，利用sDNA2的分子数与其测序拷贝数的比值和sDNA1分子数与其宏基因组测序数据分析得出的拷贝数比值相等：(1)计算出提取后sDNA1的分子数，再与提取前加入的sDNA1分子数相比，由sDNA1的得率来代表粪便样本gDNA的提取率；(2)用同样的方法，计算得出单位质量gDNA样本中细菌分子数；最后，通过换算提取gDNA所用的粪便质量，最终得出单位质量粪便中细菌的分子数。本发明的方法的大致流程如图1所示，本发明的结果的计算过程如图2所示。

下面详细描述本发明的方法的具体操作流程，以更好地理解本发明。

1.粪便样本gDNA的提取

本实验中选取了3例粪便样本进行测试，gDNA提取所用的是Perkin Elmer粪便gDNA提取试剂盒(CMG-1076 chemagic DNA Stool Kit)，样本前处理完成后在PerkinElmer自动核酸提取仪(PE chemagic 360)上完成粪便样本gDNA的提取。

(1)粪便样本前处理：

a.从粪便样本品保存管(含DNA保护液(主要成分是三羟甲基氨基甲烷(Trisbase)；乙二胺四乙酸(EDTA)；氯化钠(NaCl)～常温饱和；二甲基亚砜(DMSO)；乙醇；甘油，具体请见申请号为201910327561.5的专利))中用阔口枪头吸取400μL粪便混悬液，装入2mLEP管中，随后依次加入1000μL裂解液(包含50mmol/L Tris-HCl；1.0mmol/L EDTA；150mmol/L NaCl；0.1％SDS)、30μL蛋白酶K、适量镐珠，10ng sDNA1；

b.放入组织研磨器60Hz研磨60s；

c.70℃金属浴恒温加热10mins，95℃金属浴恒温加热5mins；

d.离心机最大转速(>15,000g)离心5mins；

(2)样品盒加样

a.样品盒2中加入600μL步骤(1)d中上清液、600μL Binding Buffer、600μL无水乙醇；

b.样品盒3中每孔加入75μL磁珠；

c.样品盒8中每孔加入200μL洗脱缓冲液(Elution Buffer，10mM Tris-HCl和1mMEDTA溶液，pH8.0)；

(3)自动核酸提取机器设定

a.打开自动核酸提取机，打开紫外灯(UV Lamp)提前消毒5-10mins；

b.打开桌面Chemical360.exe，待2-3mins后程序就绪完毕；

c.打开隔离板，按照仪器指示位置依次放入1-8号样品盒；

d.选择程序Chemagic DNA Stool 360 WB Prefilling Drying VD180907.che选择孔位置。

(4)仪器结束工作后，打开隔离板，取出样品盒8，吸出提取好的gDNA至洁净的1.5ml EP管中，即为提取好的粪便样本gDNA；妥善处理好使用过的样品盒后，检查机器无异常即关闭程序，关闭仪器。

2.宏基因组文库的构建以及测序

本实验中所使用的宏基因组建库试剂盒为华大智造生产的MGIEasy DNA文库制备试剂盒，测序试剂盒为BGISEQ-500RS高通量测序试剂套装(PE100)，测序在BGISEQ-500RS测序仪平台上完成。

(1)粪便样本gDNA的质控、定量及sDNA2的添加

a.使用酶标仪对所提取的gDNA进行浓度和质量(OD_260/280)检测，其中浓度小于10ng/μL，或OD_260/280值异常(OD_260/280>2.5,or OD_260/280<1.5)的样本则判定不合格，不进入标准品添加和建库实验环节；

b.取2μL gDNA进行琼脂糖凝胶电泳，15K DNA marker位置左右无明显主带或降解严重的样本则判定不合格，不进入标准品添加和建库实验环节；

c.经质控合格的gDNA样本，取2μg(根据a中的浓度计算体积)加入0.2ng的DNA2(万分之一比例)，混匀；

(2)宏基因组文库的构建

a.取1μg gDNA样本，加入Covaris microTUBE(打断管)中，用TE缓冲液将体积补至130μL，开启Covaris打断仪，选择350bp片段的打断程序，进行超声打断；

b.磁珠双选：提前30分钟取出AMPure XP磁珠置于室温，使用前充分震荡混匀；先后加入0.6倍和0.2倍体积的磁珠，纯化打断后的gDNA，获得峰值在350bp的片段；

c.末端修复:在PCR管中配制反应体系(DNA 40μL；ERAT Buffer 5μL；ERAT dNTPMix 0.6μL；ERAT Enzyme I 0.6μL；ERAT Enzyme II 0.1μL；ERAT Enzyme III 0.2μL；ERATEnzyme IV 2μL；NF H2O 1.5μL；总体系50μL)，混匀后在PCR仪上进行反应，程序如下：热盖On；37℃30mins；65℃15min；

d.接头连接及反应产物纯化：向上一步反应液的PCR管中加入5μL Adapter Mix，吹打混匀；反应混合液按组分比例添加(含adapter的DNA 55μL；Ligation Buffer I 3μL；Ligation Buffer II 16μL；Ligation Enhancer 0.8μL；Ligation Enzyme 1.6μL；NF水3.6μL；总体系80μL)，混匀后在PCR仪上反应，程序如下：热盖On；23℃60mins；反应结束后，取出PCR管，加入20μL TE到PCR管中至体积为100μL，再加入50μL AMPure XP磁珠(0.5倍)至100μL连接产物中，获得44μL的纯化产物；

e.PCR扩增及产物纯化：取44μL纯化后的连接产物，按组分配制反应体系(DNA 44μL；PCR Enzyme Mix 100μL；PCR Primer Mix 12μL；NF水44μL；总体系200μL)，PCR反应体系如下：热盖On；95℃3mins；7个循环(98℃20s；60℃15s；72℃30s)；72℃10mins；反应完成后加入等体积的AMPure XP磁珠进行纯化，获得30μL的纯化产物；

f.文库质控：PCR产物纯化后可选用Agilent 2100 Bioanalyzer或者LabChip检测PCR产物的长度和分布范围。要求PCR产物片段主带在450bp左右，无dimer、无其他污染；

g.PCR产物环化及消化处理：使用

dsDNA HS Assay Kit试剂盒对纯化的PCR产物进行定量，等量混合所有的待测文库(总量为330ng)后，补充NF H₂O到PCR管中至总体积为60μL，按照组分要求配制热变性反应体系(均一化后PCR产物60μL；Splint Oligo 10μL；总体积70μL)，在PCR仪上95℃孵育3mins后立即冰上放置3mins；之后按组分要求配制环化反应混合液(热变性后PCR产物70μL；Splint Buffer 12μL；Ligation Enhancer 1.2μL；Ligation Enzyme 0.4μL；NF水36.4μL；总体积120μL)，PCR仪上37℃孵育60mins；反应完成后，按照组分要求配制消化反应体系(环化产物120μL；Digestion Buffer 0.8μL；Digestion Enzyme I 3.9μL；Digestion Enzyme III 1.3μL；NF水2μL；总体积128μL)，PCR仪上37℃孵育30mins；反应完成后加入15μL Digestion Stop Buffer震荡混匀5s终止反应，之后加入170μL试剂盒自带的Library Purification Beads对产物进行纯化，最终获得40μL产物；

(3)文库上机测序

a.DNB(DNA Nanoballs)的制备：取出上一步骤中环化好的文库，按照组分要求在0.2ml的PCR管中配制反应液(文库40fmol；DNB制备缓冲液20μL；补NF H₂O至总体积40μL)，震荡5s混匀后置于PCR仪上进行反应，体系如下：热盖On；95℃1min；65℃1min；40℃1min；当PCR仪温度达到4℃时取出PCR管，迷你离心机离心5s后，在冰上将如下组分加入管中(DNB聚合酶混合液40μL；DNB聚合酶混合液II 4μL)；反应混合液用漩涡振荡器震荡混匀，迷你离心机离心5s后即刻置于PCR仪中开始反应，反应条件如下：热盖On(80℃)；30℃20mins；当PCR仪温度达到4℃后立即取出PCR管(切勿离心)，加入20μL DNB终止缓冲液，用移液器和阔口吸头缓慢地吹打混匀5次，切勿震荡及剧烈吹打，放于4℃保存备用(48小时内使用)；

b.DNB浓度的测定(质控)：DNB制备完成后，用

ssDNA Assay Kit和 Fluorometer仪器进行浓度检测。浓度8ng/μL以上为合格，对浓度过低的DNB需重新制备；

c.DNB的加载及处理：提前取出样本加载试剂板于室温融化后并放置好；检查全自动样本加载系统是否正常，并按仪器说明书放置好芯片，将DNB转移至适配的PCR管中，并加入DNB加载缓冲液II 32μL；DNB聚合酶混合液II 1μL；用阔口吸头温和混匀5-8次，将以上混合液放于加载系统指定DNB放置区，选择DNB加载程序Sample load 2.0，点击开始，等待加载完成；加载完成后，保持芯片在室温孵育30mins待用；建议制备好的芯片应尽快安排上机测序，如果当天不能使用，请置于干净的PE手套内并平放于2-8℃冰箱中备用，有效期48小时；

d.测序：取出PE100测序试剂槽I和PE100试剂槽II室温融化后(约3h)，置于4℃冰箱备用；提前1h取出PE100 dNTPs混合液和PE100 dNTPs混合液II，室温融化后置于4℃冰箱或者冰盒上备用；使用前取出PE100测序酶混合液，置于冰盒上备用；移取PE100测序酶混合液(3600μL)、PE100 dNTPs混合液(1800μL)和PE100 dNTPs混合液II(1800μL)前，短暂离心后用移液器吹打混匀5次，加入到对应孔位中；依照BGISEQ500-RS使用说明书，启动测序仪对应的内置控制程序；首先用清洗芯片进行仪器的流道清洗，以上预处理完成的PE100测序试剂槽I和PE100测序试剂槽II按照孔位一一对应，并使用试剂盒IV的试剂槽盖板组装在一起，然后放入测序仪冰箱，参照BGISEQ500-RS使用说明书进行试剂预载，之后安装芯片开始进行测序；MDA试剂加载如进行单芯片测序，则待PE100测序进行在60-100cycles间加入试剂，如进行双芯片测序，则待PE100测序进行在80-100cycles间加入，仪器上界面会出现相应提示；测序完成以后，取出芯片和试剂槽，参照BGISEQ500-RS使用说明书清洗测序仪。

3.下机数据质控及分析

(1)在此次测试实验中，使用了3个样本的数据，经过初步的处理，数据见下表1。

表1

(2)测序样本中标准品片段及细菌的拷贝数计算

经过数据进一步比对到参考基因组上，通过软件计算得出菌群的拷贝数值见下表2。

表2

(3)测序样本细菌分子数的定量计算：

测序样本单位质量粪便样本分子数计算结果见表3。

表3

项目	KMS1	KMS2	KMS3
				测序所得sDNA2的拷贝数(copies)	15.6849602	18.15123586	19.48387097
测序所得sDNA1的拷贝数(copies)	50.09312053	174.1325261	27.07253715
				实际1μgDNA样本中sDNA1的分子数	55391665.8	166388559	24099247.05
DNA提取率	0.004547156	0.127477542	0.0025409
				细菌拷贝数(copies)	3.866859231	4.467251411	4.548159557
1μgDNA测序样本中细菌分子数	100727582.1	4268585.206	4048649.751
				400μl样本的DNA提取总量(μg)	6.91257	0.74067	5.38209
400μl粪便提取的DNA所含分子数	1.53126E+11	24801333.32	8575780792
				原采样管粪便质量(g)	1.3	0.6	0.36
原采样管粪便体积(ml)	2.8	3	2
				400μl粪便样本的质量(g)	0.185714286	0.08	0.072
每克粪便的细菌分子数(个/g)	8.25E+11	3.10E+08	1.19E+11

其中，sDNA2加入1μg DNA中的分子数为：17344019.78；sDNA1在gDNA提取时加入的分子数为1762239883，那换算到每提取1μg DNA所加入的sDNA1的分子数＝sDNA1在gDNA提取时加入的分子数/400μl样本的DNA提取总量。

从KMS1和KMS3这两个样本实际计算所得的单位质量的细菌数来看，与常规文献里所得的结果数量级相符(文献中报道的单位质量细菌分子数数量级在1010-1011，参考文献：Vandeputte,D.,Kathagen,G.,D'Hoe,K.,Vieira-Silva,S.,Valles-Colomer,M.,Sabino,J.,Wang,J.,Tito,R.Y.,De Commer,L.,Darzi,Y.,et al.(2017).Quantitativemicrobiome profiling links gut community variationto microbial load.Nature551,507-511.)。

如上所述，通过本发明的方法获得细菌的绝对定量，不用采用高成本的流式细胞计数，也不用对特定的DNA进行定量PCR，简化了过程并且降低了成本。

本领域技术人员应理解，以上实施例仅是示例性实施例，在不背离本申请的精神和范围的情况下，可以进行多种变化、替换以及改变。

序列表

<110> 康美华大基因技术有限公司

<120> 一种宏基因组绝对定量的方法

<130> HP190001LZ

<141> 2019-11-13

<160> 2

<170> SIPOSequenceListing 1.0

<210> 1

<211> 5530

<212> DNA

<213> 人工序列()

<400> 1

aataaccgaa agcctctgga actcaatacc aagcctaatc aaagtcgttg agacttgcac 60

ttgcagatca ccgtgcttgc ttatcgtgac attaagaaca tcaccaacat gtctcatatg 120

atcaacaaag ttctgtaact gattcgagtc ttgaacttga acaagagtcg gaggaaggtc 180

ttgagccgag tcaagagcag tctgcagctc gataacttga gaccgagata acggtttcgt 240

aatcggaaca tcttgaataa ctgcagactt gtatcctttc gtttcaaacg tgaggaaagg 300

catcggttgt gtacagttag gagggagctt tttgacaagc ttgatctgca gtctgtttga 360

agctaatcca ccactgaact cagtacaaat actaacactg ctcttcaccg ctcgatacag 420

aagagcaaca tctaaggaga aagcgatacg gtcctcgttc tggctcgaga tacggtaatc 480

atcaaagaga acatctttgc gaaactgagc aatgcattga actccatcac cattgagaag 540

gttatggagg aagaacaagt gttctttagt gagaagaaga tgacaattct tccccatttt 600

gtcaaacgct ggaagaaacc tcttctctaa gagattcaca ccattttcag tcagaaacgc 660

cttgaacttc atttcactgt tacaaagaaa caagacgatg agtgaatctt tagtagaaat 720

cactctacga actcaaattt tggcacagat ttcggagaaa aagcgattct tcaaaactaa 780

aagaccagat tcaatcaaac acaagactat tcaaatttag gttaaattcg aatgcaatca 840

aaacatcgga ttgtatgaaa ttgaaataat taaacgattg gagaatcgtg aattgttaga 900

atcgtttcaa gtatacagtt cataaactca ccggaatcaa tagtttgtaa tattctgacc 960

ggagatcacc ggaatcgtct gctggatgaa ggaaactgta gctgcgttga agaaatcaaa 1020

caatctatat caaaattcat gagcacagta gagaaagatg acagattcag agaatttagg 1080

aattcaccaa ctaagaaatt ttctgaattt ttgcgccaaa aacaactctg agattgaagg 1140

agagggaaac tggtagattt ttatttcttc tgaaattttg tgtaccaagc ccattaccat 1200

tagtgtcata tgcatacatg ctcaaaagcc cataggggta aaatgggaaa ctggttataa 1260

ctttttaaag aaaaaaaact taagtttgaa gaacttttga atgatcaaat ttatttacca 1320

agctgcctat ttgcattatt tttggcttca ttaatctgtt attacagctc tatgtcaaga 1380

tataaaattt tacagcaaac aaaaattaaa ataagctcct tttctttata tacatacaga 1440

tatagtttca gaaattcttc tataaacact agtcttgttt ttcttctaca ctatgaactt 1500

tatccaaaga ctcttgaaga acttgcattg cgactttata atgatccaaa gataccgctc 1560

cttcttcaac tacctgaaga gagtttttat acaactggtc aaaccattga acccggtccg 1620

taccatcgac aaaaccagtc aagccgttat cagcgaagtg gaggcgtttg tagtcttttc 1680

tccaccgagg tagaatatac cgtaaaggga tttcttcaac gccgttgaaa ttgagaacac 1740

atagagcatg tctgcataag tatccgtaga agttgaagca gctacaaatg caacgtactt 1800

cgccaactga tctgttgtac agaacttcga aatcgcggat ttctcttctg ctggattcgc 1860

ctcggacccg ttctttcaca aggaagatca cgaaaggtcc atcgacatgg acctgcgttg 1920

tgctgaaaca tgagtacatt tcttctacct ctatttggaa tttcttgaac atgtctctcg 1980

tgtagattct cgagagctgt gtctcgaaag aacactttgt tttaagctcc gcggtgttca 2040

aggtttgaga ctctatgtcg ctgagagttt cttccctgtg cttcttctga agagctagct 2100

cgtacttatc gagaaactct ttcaatggtg tttgcttgtg aacatacctc tcgaaaaatg 2160

gtttcagggt ttctccagga tgagcagcgg cgattcctgc aaaaaaagtg tctttgagat 2220

aaaccggagc ccatttagct ctctcttcat acaaggaacg gagccattca ttctcgataa 2280

ccccaaaatt gtgaaccatg aacccccaag ctgcttcaaa ctcaacaact tttaaggttt 2340

cgtaaactgc tttcgtaaat gccttcctca cggcatcata attatgcagc ccccctaatt 2400

tctctggaat tttcctcatt atatgtgtca gactgaatct ctgatgagat cttgggaaca 2460

cttgagagat ggcagcctct aaaggcttgc atctgtctgt gacaatagtt tgaggagaac 2520

gcttcatgac actgagccat actttaagca accaatggta agactccatt gtctctccag 2580

caagaaaacc acaactaagc aatgtagttt tcccatggtg attcactccg gtaaatgtca 2640

ccagtggaat ctcgaatttt cctgatatat aagaactgtc gatgaagatg acatcgccaa 2700

agtaactgca ggaaacttta gagaaagcat ctgcccagaa tacattccgt agttgacctt 2760

catcattcac atccatcaag taaaaaaagt ttggatttgt caattgcata cggcaaaagt 2820

agttgtaaat ggctgcagaa tctcctctct taaggttcag caaatcaggt gacccagtag 2880

aattctgaaa cttcttgtta agagtcgagt ttgggttcac attgctccca ttatccacca 2940

cacaagctcg gtacagctta atcgtcttag cgtctgagac tggactagat acacacttcc 3000

tcttcctttt aaccgattta tataacttac aaccgagtaa atggttatga tcaagtgtaa 3060

cttcaactac ccgccacctc ttagaatcaa cctgtctcat tcttatcata gcaggacaac 3120

cagttcgagt ttcctttcta acacgattca catcgttaat cctcttaaaa ccttgactac 3180

tacaacaaag caccgcacca tacttctctt tactacgccg cttaaaccac gaattcttca 3240

ctcttacacg aaacccaact tcactagcat aacagttata atagttataa gcatcatcat 3300

acgactcaaa ctccattcca accgcaggtg catcgaattc ttttctctca tcaacaagta 3360

aaccactctg ttcatcaaac cctcccattg tttcttcttc actctttctt tcttcgtttt 3420

gacaaaatcc aacatcaatc tcatgctcac tactacaaac tagagcttct tctcctatcc 3480

gtaaagtatc tgacattgga tcatctcttt cacgaacctc gtcgttctct aaatatgcag 3540

atgcttggac ctgtaataac aaaaagctat aaactttgat tctagcaaac cctatttact 3600

caattgaggc atcaaacttg tgaaggagct gaacaaaaat gataatcaac ccaatgaaat 3660

caaatgtaac atcatcacaa cccttgtgcg aatcgattca taatcttcac aaaatccaga 3720

attgaaaaga ggagttttag agttaaacac atacatcttc atcaacggac tcacttctct 3780

ccattggagg tttggttcct cttccttctt caattcgcgg gtaaattttg atgaaagcac 3840

ttgactttat tcgattgaag gagaacactc gtgaaattga tcggtccggg atagatccgg 3900

acgggtctta cacggaaaga tccgacccga gagtgatctt gcgcgggtat tgacctttga 3960

aagcggttag gtagtttaac ttcttcttct gacggtgagg tggtcttttg gggctcgact 4020

ggtgactcta gcattttcac caacatagag acaccttagg cttagggttt tcacacacag 4080

gttcaaagca cgtgctgcaa acgtcgtcgt tcgcacttaa acatatcatt tttgtttttc 4140

cttcacttaa gacacttaaa cgacaattac actgaaacgg agggaatgct atgaaatgaa 4200

atccaaaatt gcattacata gattccctga aatgaaatct atgagtctat gacattagta 4260

aaacgcacaa acttgtgaat gtgatattta tgagcttata aaacagcata tataaaaaaa 4320

aaatatatca atttctagat tttacccaca atcttaattt atgtcatcat taaccactaa 4380

tttttcaata ttttctctct ctttcgaggc gcatatataa atacaggaaa tacgcttttt 4440

tacttaccaa tcaaaaatgt ctttgctctc tcctcttcct ctgctccact cattctcctc 4500

caccgtcgcc acaaaatcca ccgcctctag aatcactgcc acgccgtcga aaatccgctt 4560

ctctgtcata aacgccacgt cggagaatgg aaacagcggt ggtagtaaaa acgatcgtga 4620

tgaggatcca tcttttaacc catttggttt cgttacagac aacccgtcta gccggagtgc 4680

aattcagtta ccagagtctc cggcggaaga cggcaatgtc ggccaaatgc tctacgtaaa 4740

tttcacatat ttactcagtt caattttgaa ctgtgattga aattgaaaag atgcatggtt 4800

catatttgct tctattgatt tcaaattcca ttttttatca tagtgttata tttttgttga 4860

ttgttatgtt gtgggtttag aggactgaag ataaaggtaa agagtatggc tcaactatca 4920

aatctgggaa gcttagatgg tttgtgagag agactggtaa agacttaagt tttgtgtttg 4980

ttttttggat tcaagaagag tttttgatca ctgattaagc attttcctat tttccggtag 5040

gatcaaagga gagtcgtcga ggaaccattg tatttgttca tggagctcct actcagtctt 5100

ttagctacag aactgtcatg tctgaggtgc tacactgaat cttgaatcat gtgtctatat 5160

gggaaacaat tgtttctaat gaaagactct tatgatttcc aatgaatcat gtgtatgtat 5220

atgaagcttg agtggatctt ttgtgtgttc ccaacaatgt aatagttgtc agatgctgga 5280

tttcattgct ttgcacctga ctggatagga ttcggattca gtgacaaacc gcagcctgga 5340

tacggtttta attacacagg ttcgttgatg gtacatttaa tagtttgtga actgcaaaat 5400

cagtgtacat ggttctttgt tctaataagt tttaaatctt ctcagaaaaa gagtaccatg 5460

aggcgtttga taaactgctt gaagtgctag aggtcaaatc tcctttcttt cttgttgttc 5520

aggttcggat 5530

<210> 2

<211> 5619

<212> DNA

<213> 人工序列()

<400> 2

gtgcttatct tctcggtggt tgcagccgct gtggctgcac cgcctgttgt ggcggctata 60

gctggtgcat tggctgttcc tgtaggttct gttgggaaat ggtgcaacac tttgtggacg 120

aagtatgaga aagtggttag agggcagaag gagattatta catcgattag aatcgggact 180

tatatctcgg ttaaggagat ggataatata agtattcttg tgcgtaaagt ggaagtagag 240

attgaatctt tgttgaagaa agctgagttt gctattactg aggagaaaga agtgaggctt 300

gctatagatg agatcaagaa gaagctcgat gtgtttactg aaaccattga agaacttggg 360

gaacatgcgg gtaagtattg tagcgatgtg acaaaggcga gaaccgtgat tctgcagagg 420

ataatccgat acccggctgg ttcaccgaag gatgaagctc cttggtaaga aatcatttgc 480

ccctaaagct aactctgaat agtaaactta tgcttcatgt ttttgctata gaattctaga 540

gaacaaactt gtgaacatga ttcatcacta ggatttatag atgatttgga ttgagttttt 600

ttttcgtgtt ttcgattgat atcaggacgg agatgatgta gtgagaatgt ggattcggat 660

tcagatttgt gtttggcttt ggaaaattgt tggtttgtca tctttagcac aagtattatg 720

aacagctaaa ggataatctt tacctaatgt tattagttgc cctgaggaaa cctttgttta 780

atacatttaa ttaaagagac acttgtaagg tttctatgtc atttattttg ttggctttct 840

ttaatcagtg ttgtatgttt ataaatgaat gtctatgtgc ttagtcaagg cttggagaaa 900

ttataagctg gatccttaaa gaatcaatca attttaaaaa taggacgtga agagttagtt 960

gatgagacga gggagggatg tttctcatag acaataggga cataatacta cttcttaagc 1020

aacatttgcg atgttttaat gttttaatgt tctttaacgt taatttggaa acatgtttta 1080

cgtgacagac gtggtctatg aaagcccatt ttgacagaag catgatataa aacccatgtt 1140

tcaagaatat ggtgttgaat gtaaaaaggt ttaaatgaat tgatttagct acctatgtca 1200

tcaaagtcca cttatctttt cggtcaaaag tctttagaga atttcatagt taaaattatt 1260

attcagtttt aaagtttgtt ctctttcaat tttaaaatta ttgtgttaga atttgagtcc 1320

atgattaaat cttatgatat ctttatgtta gttcgtatgt aaagtcaatc ggtttgtgtt 1380

aatcgtgact tgttaagaaa gggcgtgagt aagagttgtg tcttgttgga taagtgcaag 1440

tgttttccta atgtaaagta tcgtcgagta aataatcaaa tttcagattc aataaagagt 1500

tatccttgtt cttgtcttct ttacaagaga agatacaaat gtctcttgat ctgtttttct 1560

tctttctcta aacattccaa cataacccag gaggatactt gaattcttga agcatataaa 1620

gagatgggaa caacagacca tgaccatcaa tgtattaata gtgtgaaaga tgaaattaga 1680

aactccgtat agtttcatta caaagttggt atccttcaac aaaaataaag gacatatcaa 1740

acacgaaatc attttggcta ctaagcacaa tggaaactac aagctctcaa taaagaaact 1800

tctatagttt gtgctggtca cattcatatt ccactgttac tatatatttt gacaacatta 1860

tcacatctta acgtaacttg tttagcaaac aacaattttc agtggtttct gttggcaatg 1920

aagttagtaa gagctatcaa aggtttgact ttgtcagaat caaacccttg aagatgttct 1980

ctcgcatctc tcttcaactt ctccacaaac tccttagatt tctcaacacc aatcactttt 2040

ggatacgtca gctttccggc gaccttatct tttccggcgg tttttcccaa ttcctccgac 2100

gacttggtct catccaagat atcatcaacc acctgaaaca acaatccgat acacctcgca 2160

aagtttctca ctttctctat ctctttctct gttccaccgc ctataaccgc tccgataacc 2220

gcagatgctt caagcaacga accggtttta tgaacgtgaa tatattcaag atcctctaat 2280

ccgacatcgt tttgttccaa cccttcactg cttaaatcct tcgcctgtcc cgctacgagt 2340

cctttcgttc caatagactt agccaattct ttaacggctc gaaccatttt cttggatgac 2400

acgtcagctt ccgtcaaatg ctcgaaagct agagctaaaa gtgcgccgcc ggagagtatg 2460

gctacgcttt ctccgaagac tttatgtgtt gttggtttcc ctcgtcggag atcgtcgttg 2520

tccatacaag gaagatcgtc tttgatgagt gacattgtgt gaatcatttc aactgcacaa 2580

gccgcttgaa tcgctaaacg ctcttccccg ccgactagct cgcaagcggc taggcagagt 2640

attggcctga cacgttttcc gcctgcgaga atcgcgtaac gcattgcctt gtggatgttg 2700

agtggctcaa taagtggaat tgcttcatct aatgctctgt ttattgtctc tgctttgctg 2760

atcatgtacg gtttgaaatc gaaactgaaa ccgaaatcgt cttcttcgtt tactaatttg 2820

attgatgaat caagaaactc acgggtgtct tttctagcga aaatagcagc ggcgaccgga 2880

gatttgatcg ggatatggag aaggcggata agccgccgtg tattcgcggg actcagccgt 2940

cctttcagag ttgtaaacac gaagtgtaga gatagaaaga cgatcaaaag atagagaaag 3000

atattttgag cttccatttt gatcaacaat ataggtgatt gattgatctt tgtttcttct 3060

tctttgtaat aactttttgt gttctgtttt tttttaaggg gttttaggaa aaggtcaatg 3120

gcagggcacg tcgtgggatc atatcattct tggtaagatc ggttaaaact gaaattattt 3180

gaagcgtggt ttggatccaa gccgtttact ttgaagttcg gtgtggagtg aagtgagcca 3240

acgagtgcat agtttttttt ttttccttct aaattcattc ataccaaatc aacattacag 3300

aagacaaaaa ctttcgaaaa caaatgagta cataatacca tattgagaaa aatatacatt 3360

aagtgaattt aacgatactt tccacgtata agcttttgag tcatagtgta aatactgaac 3420

aagtgagtaa gacttgtttg tttatcgact tttacaacca tatatataag aaaaggttta 3480

caaagaaaac caattcgaac cgggagaaaa gataaatcaa accggagttt tgcatgaggc 3540

tgatggtatc ttcattgcaa atgaagcgga gaaacacgac ctccctagtt ttcatcaact 3600

tttcacatct ggttttgttt tcatgtggtt tattttttat aatcgtttct caagtttcag 3660

tggaggattt aaagatttat gatattggat atcaacgtgc tcgacgcatg gtgtcctaac 3720

gttcaacaat tagtgtcgta acgttgaacc tgtctttttt atattatagt aaatagtaac 3780

attttggatg ctttctttta tgaaacaaag tcaacgaatg caagtctctt cgtggaatct 3840

ggtcatatat cggttttatg ttatcggaat atgttttaaa ccatacatgc gttgattaat 3900

catactgtac ttgtttcacg gcaagatggg agttactagt cccttattct ccgacccgaa 3960

ttttagtttg gttaattttg ttaccttagc ccttaaccga aataatggaa aaccggaata 4020

cctgaaccga gttgaaccaa acttggtcct tagttctttc tagtccatat cagaggttat 4080

actatttata cattagagcc tcatcaatgt gttttcgatt gattttgata cacataagcc 4140

ctccataaac tacaaactat ggaggttttc agcacactca tttgtaacat ttaattttta 4200

ttttcagtaa catttaagat tttattagat tagatatttt attcaaggat tttgtgcaag 4260

aaaacccttc aactcgttaa tttttcgatt ttaccatcta tttccgaaat gtaaattaat 4320

tttatgaagt caaaaataat acgcaatata acacttatac aaagacgaat actaatatgc 4380

ctaaagtaaa aaacaagttc taatttacag aaaacatata ttaagatata agagttaaga 4440

aaaatctatg taaaatttta gtaattatta cgttaagcat aaggcgaatg ctaatctgcg 4500

taaggtaaaa aaacaatttc taatttacag attttttatt aagatctgcg taaggtagaa 4560

atgaattata atttgtggaa cgaattctag gatgcgttaa gtatctaatt tgtattgttc 4620

tcgtgcgaaa aaacaaaaaa aaaaactaca ttgcgtatta ttgttgactt tttgtgatta 4680

atgtaactat gatgaattag aatgaaagga tccatccaaa ttcttacaag acaaagtcat 4740

aaattcataa tataaagact ctggtaatta acatgattgc gatacattaa ctttgtttaa 4800

actgctcaag taaggagatt agaacaatgg ctaaaatgtt tgcagaccga ccaaattttc 4860

aatgggccta agaagtacga accctacatt tctaattaat ctcccactat tcatctactt 4920

tctctctaat gaaaacacac aaaatttcca aacacaattg taaagaaaaa aataaaaatt 4980

taataagcaa aagaaaataa aaacctagac attggatgta gaaggctgaa caatcttgtt 5040

ttgtctccct gcaagccaat gcatggtcct cctagcaaaa gatggcgttt tcgtctccaa 5100

agattccgaa ttttcgctct taagactctc ttcgggagga atccgatggt ttggaaaccg 5160

gagtgaagag gaagacacga aagaacaatc cgaatctcta gcactagtcc ccttttgctt 5220

ctccttttta ccccttagcc aaattttcga taaagaaaag ctctccctct ccacaacact 5280

cttctcaatt cccttaactg ttggatcaaa accgcattcg gacataacgg tcctgtggct 5340

cggcacccgg tcttttccgg atagtttctt ggtcgaaaca tgaaccttaa gtgtagcttc 5400

ttgatccata atatactcat atgatcccat tgacaaacac ctccttccat ccacattttt 5460

actattacca cttatgctat ttttcttctg gtcacttcct tcaccaacat gatctatgtt 5520

cctaaacttc cctaacttaa cttctaaagg aacaaccttt tcatctaaac catccatgtc 5580

accatccggg tcacatgatt tacgacccga cctaatccg 5619

Claims

1.一种细菌定量的方法，包括：

在样品中添加第一定量DNA片段；

对所述样品进行基因组DNA的提取；

在所述提取的基因组DNA中添加第二定量DNA片段，得到总基因组；

对所述总基因组进行宏基因组文库构建和测序；

利用得到的测序结果的数据计算单位质量的所述样品中的细菌数量。

2.根据权利要求1所述的方法，其中，所述样品为粪便样品。

3.根据权利要求1所述的方法，其中，所述第一定量DNA片段的序列为序列1。

4.根据权利要求1所述的方法，其中，所述第二定量DNA片段的序列为序列2。

5.根据权利要求1所述的方法，其中，所述样品中的所述第一定量DNA片段的分子数为：测序所得的所述第一定量DNA片段的拷贝数与单位质量的所述样品中加入的所述第二定量DNA片段的分子数的乘积再除以测序所得的所述第二定量DNA片段的分子数。

6.根据权利要求5所述的方法，其中，所述样品中的基因组DNA的得率为：单位质量的所述样品中的所述总基因组中的所述第一定量DNA片段的分子数除以加入的所述第一定量DNA片段的分子数，之后再乘以所述样品中的基因组DNA的总量。

7.根据权利要求6所述的方法，其中，单位质量的所述样品的细菌的数量为：单位质量的所述样品中加入的所述第二定量DNA片段的分子数、测序所得的细菌的总拷贝数和所述样品的基因组DNA的总量的乘积，再除以测序所得的所述第二定量DNA片段的拷贝数、所述样品中的基因组DNA的得率和所述样品的质量的乘积。