CN112609024A - 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用 - Google Patents

一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用 Download PDF

Info

Publication number
CN112609024A
CN112609024A CN202011628621.6A CN202011628621A CN112609024A CN 112609024 A CN112609024 A CN 112609024A CN 202011628621 A CN202011628621 A CN 202011628621A CN 112609024 A CN112609024 A CN 112609024A
Authority
CN
China
Prior art keywords
patients
microbial
infected
model
microorganism
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011628621.6A
Other languages
English (en)
Other versions
CN112609024B (zh
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.)
First Affiliated Hospital of Zhengzhou University
Original Assignee
First Affiliated Hospital of Zhengzhou University
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 First Affiliated Hospital of Zhengzhou University filed Critical First Affiliated Hospital of Zhengzhou University
Priority to CN202011628621.6A priority Critical patent/CN112609024B/zh
Publication of CN112609024A publication Critical patent/CN112609024A/zh
Application granted granted Critical
Publication of CN112609024B publication Critical patent/CN112609024B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/70—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving virus or bacteriophage
    • C12Q1/701—Specific hybridization probes
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869—Methods for sequencing
    • C—CHEMISTRY; METALLURGY
    • 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/158—Expression markers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/30—Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change

Landscapes

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

Abstract

本发明属于生物医药技术领域，具体涉及一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用。肠道微生物基因标志物由SEQ ID NO:1‑7所示的7种微生物基因组成，所述微生物基因在人体肠道中富集。本发明的微生物基因区别模型在新冠肺炎感染患者和非感染者中实现了良好的区别能力，也验证了该区别模型能够弥补RT‑PCR检测的不足，将感染新冠病毒的疑似患者从健康人群的鉴定出来。

Description

一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因 标志物及其应用
技术领域
本发明属于生物医药技术领域，具体涉及一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用。
背景技术
新型冠状病毒肺炎(Coronavirus Disease 2019,COVID-19)在世界范围内迅速传播，截止到2020年12月，已有超过5900万例感染病例和100万例死亡病例。新冠肺炎的广泛传播给全世界的卫生健康事业带来了极大的负担。通过逆转录聚合酶链反应(reversetranscription-polymerase chain reaction,RT-PCR)进行核酸检测来判断是否感染新型冠状病毒(Severe Acute Respiratory Syndrome Coronavirus 2，SARS-CoV-2)，是目前公认的诊断“金标准”。但是，由于病毒窗口期等方面的原因，感染新冠病毒后，部分患者经RT-PCR检测出现假阴性的结果，且假阴性率高达20％，这极大的增加了疾病的传染源和传播范围。因此，建立一种新的区别健康人和新型冠状病毒感染者的诊断模型，以弥补现有检测方法的不足，对于新冠肺炎防治工作具有重大意义。
人体肠道微生态系统与健康和疾病密切相关。人体肠道内定植着数目庞大(1014)、结构复杂(超过1000种细菌)的微生物群落(约1.5公斤)。其细胞总量几乎是人体自身细胞数量的10倍，基因总量是人类自身基因的150倍。肠道菌群和宿主共生并共进化过程中，在调节宿主的消化吸收、免疫反应、代谢、炎症等方面发挥重要作用。肠道微生态失调促进慢性疾病的发生发展，包括肝硬化、心血管疾病、糖尿病、结肠癌、慢性肾脏病等。
肠道微生物中关健功能菌可成为人类疾病的新型生物标志物。肠道微生态的特征或者基于肠道微生物建立的区别模型作为特定疾病或肿瘤的区分工具正越来越多的被广泛报道和认可。肠道菌群的“肠型(Enterotypes)”可反映人对疾病的易感性，提示肠道菌群具有潜在的预警和诊断作用。Qin J等率先使用宏基因组学报道了肠道微生态与2型糖尿病的相关性，指出23个菌种可能成为区分2型糖尿病的肠道微生物标志物。Qin N,Li A等(2014年Nature杂志)通过宏基因组测序技术解析了中国人群肝硬化患者肠道微生态结构，在微生物基因和功能水平上，鉴定了肝硬化特异性的15个标记物，创建了一个高准确度的肝硬化患者区分指数。Ren Z等(Advanced Science,2020)指出，肠道微生态对慢性肾脏疾病更为敏感，5个菌种可成为区分慢性肾脏疾病的肠道微生物标志物，并且肠道微生物也可反应慢性肾脏疾病进展的不同阶段。Wang Z等(Nature,2011)指出依赖于肠道微生态的卵磷脂代谢促进了心血管疾病的进展，这一发现为心血管疾病新型的诊断方法和治疗策略提供了理论基础。Oh PL等(AJT,2012)发现肠道微生态的变化与小肠移植排斥紧密相关，可成为一个潜在的移植排斥的诊断标记物。Yu J等在不同民族患者中揭示并验证了结直肠癌的微生物标志物，指出微生物标志物作为一种可支付的、无创的结直肠癌的早期诊断标志物。因而，肠道微生物可能是不同疾病诊断的有力工具。然而，用于区别新型冠状病毒感染和非感染者的肠道微生物模型还未曾报道过。
发明内容
本发明提供了一种用于区别新型冠状病毒感染和非感染者的肠道微生物基因标志物，由SEQ ID NO:1-7所示的7种微生物基因组成，所述微生物基因在人体肠道中富集。
SEQ ID NO:1-7所示的7种微生物基因的序列如下：
>OTU470
CCTACGGggggCTGCAGTGGGGAATATTGCACAATGGgggAAACCCTGATGCAGCGACGCCGCGTGAGTGAAGAAGTATCTCGGTATGTAAAGCTCTATCAGCAGGGAAGAAAATGACGGTACCTGACTAAGAAGCCCCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGgggCAAGCGTTATCCGGATTTACTGGGTGTAAAGGGAGCGTAGACGGTGTGGCAAGTCTGATGTGAAAGGCATGGGCTCAACCTGTGGACTGCATTGGAAACTGTCATACTTGAGTGCCGGAGGGGTAAGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGCTTACTGGACAGTAACTGACGTTGAGGCTCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTTGTAGTC
>OTU1000
CCTACGGGTGGCTGCAGTGGGGAATATTGCACAATGGgggAAACCCTGATGCAGCGACGCCGCGTGAAGGAAGAAGTATCTCGGTATGTAAACTTCTATCAGCAGGGAAGATAGTGACGGTACCTGACTAAGAAGCCCCGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGgggCAAGCGTTATCCGGATTTACTGGGTGTAAAGGGAGCGTAGACGGTGTGGCAAGTCTGATGTGAAAGGCATGGGCTCAACCTGTGGACTGCATTGGAAACTGTCATACTTGAGTGCCGGAGGGGTAAGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGCTTACTGGACGGTAACTGACGTTGAGGCTCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTAGTAGTC
>OTU1092
CCTACGGGTGGCTGCAGTGGGGAATATTGCACAATGGgggAAACCCTGATGCAGCGACGCCGCGTGGAGGAAGAAGGTCTTCGGATTGTAAACTCCTGTTGTTGAGGAAGATAATGACGGTACTCAACAAGGAAGTGACGGCTAACTACGTGCCAGCAGCCGCGGTAAAACGTAGGTCACAAGCGTTGTCCGGAATTACTGGGTGTAAAGGGAGCGTAGACGGTGTGGCAAGTCTGATGTGAAAGGCATGGGCTCAACCTGTGGACTGCATTGGAAACTGTCATACTTGAGTGCCGGAGGGGTAAGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGCTTACTGGACGGTAACTGACGTTGAGGCTCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCGAGTAGTC
>OTU1741
CCTACGGGTGGCAGCAGTGGGGAATATTGGACAATGGggggAACCCTGATCCAGCCATGCCGCGTgtgtgAAGAAGGCCCTCGGGTTGTAAAGCACTTTCAGTGAGGAAGAACGCCTGGTGGTTAATACCCATCAGGAAAGACATCACTCACAGAAGAAGCACCGGCTAACTCCGTGCCAGCAGCCGCGGTAATACGGAGGGTGCGAGCGTTAATCGGAATTACTGGGCGTAAAGCgcgcgTAGGTGGCTTGATAAGCCGGTTGTGAAAGCCCCGGGCTCAACCTGGGAACGGCATCCGGAACTGTCAAGCTAGAGTGCAGGAGAGGAAGGTAGAATTCCCGGTGTAGCGGTGAAATGCGTAGAGATCGGGAGGAATACCAGTGGCGAAGGCGGCCTTCTGGACTGACACTGACACTGAGGTGCGAAAGCGTGGGTAGCAAACAGGATTAGATACCCTAGTAGTC
>OTU1314
CCTACGGGTGGCAGCAGTGGGGAATATTGGACAATGGGCGCAAGCCTGATCCAGCCATGCCGCGTGAGTGATGAAGGCCTTAGGGTTGTAAAGCTCTTTCACCGGAGAAGATAATGACGGTATCCGGAGAAGAAGCCCCGGCTAACTTCGTGCCAGCAGCCGCGGTAATACGAAGGgggCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGACTATTAAGTCAGGGGTGAAATCCCGGAGCTCAACTCCGGAACTGCCTTTGATACTGGTAGTCTAGAGACCGAgagagGTAAGTGGAACTCCTAGTGTAGAGGTGGAATTCGTAGATATTAGGAAGAACACCAGTGGCGAAGGCGGCTTACTGGCTCGGAACTGACGCTGAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCCTGTAGTC
>OTU1473
CCTACGGGTGGCAGCAGTGGGGAATATTGCACAATGGgggAAACCCTGATGCAGCGACGCCGCGTGAGCGAAGAAGTATTTCGGTATGTAAAGCTCTATCAGCAGGGAAGAAGAAATGACGGTACCTGACTAAGAAGCACCGGCTAAATACGTGCCAGCAGCCGCGGTAAAACGTAGGTCACAAGCGTTGTCCGGAATTACTGGGTGTAAAGGGAGCGCAGGCGGGAGAACAAGTTGGAAGTGAAATCCATGGGCTCAACCCATGAACTGCTTTCAAAACTGTTtttCTTGAGTAGTGCAGAGGTAGGCGGAATTCCCGGTGTAGCGGTGGAATGCGTAGATATCGGGAGGAACACCAGTGGCGAAGGCGGCCTACTGGGCACCAACTGACGCTGAGGCTCGAAAGTGTGGGTAGCAAACAGGATTAGATACCCCTGTAGTC
>OTU239
CCTACGGGCGGCAGCAGTGGGGAATATTGCACAATGGgggAAACCCTGATGCAGCGACGCCGCGTGAGCGAAGAAGTATTTCGGTATGTAAAGCTCTATCAGCAGGGAAGAAAATGACGGTACCTGACTAAGAAGCACCGGCTAAATACGTGCCAGCAGCCGCGGTAATACGTATGGTGCAAGCGTTATCCGGATTTACTGGGTGTAAAGGGAGCGTAGACGGCAAGGCAAGTCTGATGTGAAAACCCAGGGCTTAACCCTGGGACTGCATTGGAAACTGTCTGGCTCGAGTGCCGGAGAGGTAAGCGGAATTCCTAGTGTAGCGGTGAAATGCGTAGATATTAGGAGGAACACCAGTGGCGAAGGCGGCTTACTGGACGGTAACTGACGTTGAGGCTCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTAGTAGTC
另外本发明还提供了一种用于检测试剂，包括用于检测SEQ ID NO:1-7所示的7种微生物基因的引物。
所述引物序列为SEQ ID NO:8-9。
引物Primers
测序区域V3+V4：338F-805R
上游引物：338F ACTCCTACGGGAGGCAGCA
下游引物：805R GGACTACHVGGGTWTCTAAT
检测试剂在制备新型冠状病毒感染和非感染的区分检测试剂盒中的应用，所述检测试剂适用于检测SEQ ID NO:1-7所示的7种微生物基因。
适用于区别新型冠状病毒感染和非感染者，新型冠状病毒感染者包括RT-PCR检测阳性的确诊患者及RT-PCR检测阴性但IgG抗体阳性的临床诊断患者(疑似患者)。
对所述对象的粪便进行检测，以便确定该样本是否包含所述的微生物基因，是否可以建立区别新型冠状病毒感染和非感染者的微生物基因模型。
通过收集入组对象的粪便样本，抽提微生物总DNA，完成微生物DNA的16S rDNAMiseq测序，检测是否存在权利要求1所述的7种微生物基因。
进一步的，通过收集入组对象的粪便样本，抽提微生物总DNA，进行肠道菌群的16SrDNA Miseq测序。基于高通量测序数据，在训练集中建立新型冠状病毒感染和非感染者的微生物区别模型，建立新冠肺炎患病率(probability of disease，POD)指数；POD指数在验证集中计算其区别能力，进行验证；进一步在新冠肺炎疑似患者中进行验证，证实微生物基因区别模型弥补RT-PCR检测不足的可能性。
本发明还提供了一种区别新型冠状病毒感染和非感染者的肠道微生物基因区别模型，包括SEQ ID NO:1-7所示的7种微生物基因。
肠道微生物基因区别模型的建立方法具体包括：
(1)收集入组对象(新型冠状病毒感染和非感染者)的粪便样本，按照DNA的标准抽提方法完成粪便样本中微生物总DNA的抽提，在Illumina MiSeq平台完成肠道菌群的16SrDNA高通量测序工作；
(2)基于高通量测序数据，在微生物区别模型的训练集中，在24例新冠病毒感染患者和48例健康对照之间，基于一个随机森林模型，通过一个五倍交叉验证的算法，鉴定了用于该模型的最佳的7个微生物基因标志物。
(3)基于7个微生物基因标志物，通过使用随机生成的决策树的比率来计算新冠肺炎患病率的患病率(Probability of disease,POD)指数。
(4)该微生物区别模型在24例新冠肺炎确诊患者和48例健康对照之间的区别能力达到99.74％，POD指数在新冠病毒感染确诊患者中明显升高，两组之间有显著性差异(p＝1.2x10-13)。
(5)在验证集中，该微生物区别模型在12例新冠肺炎确诊患者和24例健康对照之间的区别能力达到99.31％；POD指数在确诊患者中明显升高，两组之间有显著性差异(p＝2.0x10-6)。
(6)在验证集中，该微生物区别模型在23例新冠肺炎疑似患者和24例健康对照之间的区别能力达到98.01％；POD指数在疑似患者中明显升高，两组之间有显著性差异(p＝1.8x10-8)。
因此，本发明的微生物基因区别模型在新冠病毒感染和非感染患者中实现了良好的区别能力，也验证了该区别模型能够将核酸检测阴性而IgG抗体阳性的疑似患者从健康人群中鉴别出来。
另外，还提供了一种用于区别新冠病毒感染和非感染的肠道微生物模型的试剂盒，包括用于检测权利要求1所述的SEQ ID NO:1-7所示的7种微生物基因的引物。
本发明的具体操作步骤如下：
(1)按照前瞻性临床试验的设计原则，本发明的研究设计如图1所示。该研究方案得到了郑州大学第一附属医院和光山县人民医院伦理委员会的批准。所有入组的患者签署研究方案知情同意书和临床样本收集知情同意书。
(2)每一个入组的新冠病毒感染和非感染患者提供一份新鲜的粪便样本，研究实验人员将样本在70℃灭活1小时，然后按照每份200mg的重量分装成10份，并立即冻存于-80℃冰箱。粪便细菌总DNA的抽提方法按照试剂盒的说明书进行。
(3)完成粪便细菌总DNA样本的扩增和DNA文库构建，在IlluminaMiseq测序平台完成16S rDNA测序。所有的输出序列完成基本的预处理和基本的生物信息学分析。
(4)从所有样本中随机抽选等量的序列数，按照UPARSE传递途径拼接成对应的16SrDNA基因序列分类单元(Operational Taxonomy Units,OTUs)。按照微生物基因标志物训练集和微生物基因标志物验证集，将产生的所有样本的OTUs基因序列收集整理。基于微生物基因序列，使用RDP分类器2.6版本注释。
(5)基于高通量测序数据产生的代表性序列，计算出微生物基因标志物发现集的OTUs频率文件和验证集的OTUs频率文件。这些OTUs用于一个相关性研究来鉴定在新冠病毒感染和非感染患者之间差异明显的OTUs丰度。使用Wilcoxon检验方法统计分析新冠病毒感染和非感染患者之间差异的微生物基因标志物。选择出的7个OTUs微生物基因标志物做进一步的分析。
(6)在微生物区别模型的训练集中，包括24例新冠病毒感染确诊患者和48例健康对照，使用筛选出的OTUs丰度文件，在一个随机森林模型(R软件3.4.1和随机森林软件包4.6–12)中采用五倍交叉验证的算法(除了设置“importance＝TRUE”之外，软件参数默认)进行微生物基因标志物的筛选。采用五倍交叉验证的5次试验，获得了交叉验证错误曲线，其中最小的交叉验证错误点作为cut-off值使用。最小的交叉验证错误值加上对应值的标准差为cut-off值。筛选出小于cut-off值的错误率的7个以下的OTUs标志物的集合，选择最小数目OTUs的集合作为最佳的微生物基因标志物的集合，最终鉴定了用于该模型的最佳的7个微生物基因标志物(图2)。选择出的7个微生物OTUs标志物的基因序列见SEQ ID NO:1-7。
(7)通过使用随机生成的决策树的比率来计算患病率(Probability of disease,POD)指数。决策树预测样本为“CP”，设置的参数预测为：proximity＝T,norm.votes＝T,predict.all＝TRUE。在LOO模式中构建的随机森林模型用于预测验证集中每一个样本的POD指数，最终计算每一个样本的平均预测的POD指数。
(8)使用R 3.3.0程序包中的pROC工具计算受试者工作曲线(ROC)，用来评估微生物区别模型，曲线下面积(AUC)用于指定ROC的效应值。
(9)该微生物区别模型在24例新冠病毒确诊患者和48例健康对照之间的区别能力达到99.74％(图4)，POD指数在新冠病毒确诊患者中明显升高，两组之间有显著性差异(p＝1.2x10-13)(图3)。
(10)在验证集中，POD指数在12例新冠肺炎确诊患者中明显升高，两组之间有显著性差异(p＝2.0x10-6)(图5)，该微生物区别模型在12例新冠肺炎确诊患者和24例健康对照之间的区别能力达到99.31％(图6)。
(11)在验证集中，POD指数在23例中新冠疑似患者中明显升高，两组之间有显著性差异(p＝1.8x10-8)(图7)，该微生物区别模型在在23例新冠疑似患者和24例健康对照之间的区别能力达到98.01％(图8)。
因此，本发明的微生物基因区别模型在新冠肺炎感染患者和非感染者中实现了良好的区别能力，也验证了该区别模型能够弥补RT-PCR检测的不足，将感染新冠病毒的疑似患者从健康人群的鉴定出来。
附图说明
图1.一种用于区别新冠肺炎感染和非感染的肠道微生物模型的研究设计和临床应用。
图2.基于随机森林模型采用五倍交叉验证法鉴定的最佳的肠道微生物基因标志物。
图3.在24例新冠肺炎确诊患者和48例健康对照的训练集中，患病率(POD)指数在两组之间的表达差异；
图4.在24例新冠肺炎确诊患者和48例健康对照的训练集中，微生物基因区别模型实现的区别效能；
图5.在验证集中，与24例健康对照人群相比，患病率(POD)指数在12例新冠肺炎确诊患者的表达差异；
图6.在验证集中，患病率(POD)指数在24例健康对照和12例新冠肺炎确诊患者之间的区别能力；
图7.在验证集中，与24例健康对照人群相比，患病率(POD)指数在23例新冠肺炎疑似患者的表达差异；
图8.在验证集中，患病率(POD)指数在24例健康对照和23例新冠肺炎疑似患者之间的区别能力；
具体实施方式
下面结合实施例对本发明作进一步的阐述，但本发明的保护内容不仅限于这些实施例。
下列实施例中所用方法如无特别说明，均为常规方法。下列实施例中所需要的材料或试剂，如无特殊说明均为公开商业途径获得。
本发明通过收集入组对象的粪便样本，抽提微生物总DNA，进行肠道菌群的16SrDNA Miseq测序。基于高通量测序数据，在训练集中建立新型冠状病毒感染和非感染者的微生物区别模型，建立新冠肺炎患病率(probability of disease，POD)指数；POD指数在验证集中计算其区别能力，进行验证；进一步在新冠肺炎疑似患者中进行验证，以证明微生物基因区别模型弥补RT-PCR的可能性。
其操作步骤如下：
(1)按照前瞻性临床试验的设计原则，本发明的研究设计如图1所示。该研究方案得到了郑州大学第一附属医院和光山县人民医院伦理委员会的批准。所有入组的患者签署研究方案知情同意书和临床样本收集知情同意书。
(2)每一个入组的新型冠状病毒感染患者和非感染患者提供一份新鲜的粪便样本，研究实验人员将样本在70℃灭活1小时，然后按照每份200mg的重量分装成10份，并立即冻存于-80℃冰箱。粪便细菌总DNA的抽提方法按照试剂盒的说明书进行。
(3)完成粪便细菌总DNA样本的扩增和DNA文库构建，在IlluminaMiseq测序平台完成16S rDNA测序。所有的输出序列完成基本的预处理和基本的生物信息学分析。
(4)从所有样本中随机抽选等量的序列数，按照UPARSE传递途径拼接成对应的16SrDNA基因序列分类单元(Operational Taxonomy Units,OTUs)。按照微生物基因标志物训练集和微生物基因标志物验证集，将产生的所有样本的OTUs基因序列收集整理。基于微生物基因序列，使用RDP分类器2.6版本注释。
(5)基于高通量测序数据产生的代表性序列，计算出微生物基因标志物发现集的OTUs频率文件和验证集的OTUs频率文件。这些OTUs用于一个相关性研究来鉴定在新型冠状病毒感染和非感染者之间差异明显的OTUs丰度。使用Wilcoxon检验方法统计分析新型冠状病毒感染和非感染者之间差异的微生物基因标志物。选择出的OTUs微生物基因标志物做进一步的分析。
(6)在微生物区别模型的训练集中，包括24例新冠肺炎确诊患者和48例健康对照，使用筛选出的OTUs丰度文件，在一个随机森林模型(R软件3.4.1和随机森林软件包4.6–12)中采用五倍交叉验证的算法(除了设置“importance＝TRUE”之外，软件参数默认)进行微生物基因标志物的筛选。采用五倍交叉验证的5次试验，获得了交叉验证错误曲线，其中最小的交叉验证错误点作为cut-off值使用。最小的交叉验证错误值加上对应值的标准差为cut-off值。筛选出小于cut-off值的错误率的7个以下的OTUs标志物的集合，选择最小数目OTUs的集合作为最佳的微生物基因标志物的集合，最终鉴定了用于该模型的最佳的7个微生物基因标志物(图2)。选择出的7个微生物OTUs标志物的基因序列见SEQ ID NO:1-7。
(7)通过使用随机生成的决策树的比率来计算患病率(Probability of disease,POD)指数。决策树预测样本为“CP”，设置的参数预测为：proximity＝T,norm.votes＝T,predict.all＝TRUE。在LOO模式中构建的随机森林模型用于预测验证集中每一个样本的POD指数，最终计算每一个样本的平均预测的POD指数。
(8)使用R 3.3.0程序包中的pROC工具计算受试者工作曲线(ROC)，用来评估微生物区别模型，曲线下面积(AUC)用于指定ROC的效应值。
(9)该微生物区别模型在24例新冠肺炎确诊患者和48例健康对照之间的区别能力达到99.74％(图4)，POD指数在新冠肺炎确诊患者中明显升高，两组之间有显著性差异(p＝1.2x10-13)(图3)。
(10)在验证集中，POD指数在12例新冠肺炎确诊患者中明显升高，两组之间有显著性差异(p＝2.0x10-6)(图5)，该微生物区别模型在12例新冠肺炎确诊患者和24例健康对照之间的区别能力达到99.31％(图6)。
(11)在验证集中，POD指数在23例中新冠疑似患者中明显升高，两组之间有显著性差异(p＝1.8x10-8)(图7)，该微生物区别模型在在23例新冠疑似患者和24例健康对照之间的区别能力达到98.01％(图8)。
因此，本发明的微生物基因区别模型在新冠肺炎感染患者和非感染者中实现了良好的区别能力，也验证了该区别模型能够弥补RT-PCR检测的不足，将核酸检测阴性的新冠肺炎感染者从健康人群中鉴别出来。
序列表
<110> 郑州大学第一附属医院
<120> 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用
<160> 9
<170> SIPOSequenceListing 1.0
<210> 1
<211> 440
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 1
cctacggggg gctgcagtgg ggaatattgc acaatggggg aaaccctgat gcagcgacgc 60
cgcgtgagtg aagaagtatc tcggtatgta aagctctatc agcagggaag aaaatgacgg 120
tacctgacta agaagccccg gctaactacg tgccagcagc cgcggtaata cgtagggggc 180
aagcgttatc cggatttact gggtgtaaag ggagcgtaga cggtgtggca agtctgatgt 240
gaaaggcatg ggctcaacct gtggactgca ttggaaactg tcatacttga gtgccggagg 300
ggtaagcgga attcctagtg tagcggtgaa atgcgtagat attaggagga acaccagtgg 360
cgaaggcggc ttactggaca gtaactgacg ttgaggctcg aaagcgtggg gagcaaacag 420
gattagatac ccttgtagtc 440
<210> 2
<211> 440
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 2
cctacgggtg gctgcagtgg ggaatattgc acaatggggg aaaccctgat gcagcgacgc 60
cgcgtgaagg aagaagtatc tcggtatgta aacttctatc agcagggaag atagtgacgg 120
tacctgacta agaagccccg gctaactacg tgccagcagc cgcggtaata cgtagggggc 180
aagcgttatc cggatttact gggtgtaaag ggagcgtaga cggtgtggca agtctgatgt 240
gaaaggcatg ggctcaacct gtggactgca ttggaaactg tcatacttga gtgccggagg 300
ggtaagcgga attcctagtg tagcggtgaa atgcgtagat attaggagga acaccagtgg 360
cgaaggcggc ttactggacg gtaactgacg ttgaggctcg aaagcgtggg gagcaaacag 420
gattagatac cctagtagtc 440
<210> 3
<211> 440
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 3
cctacgggtg gctgcagtgg ggaatattgc acaatggggg aaaccctgat gcagcgacgc 60
cgcgtggagg aagaaggtct tcggattgta aactcctgtt gttgaggaag ataatgacgg 120
tactcaacaa ggaagtgacg gctaactacg tgccagcagc cgcggtaaaa cgtaggtcac 180
aagcgttgtc cggaattact gggtgtaaag ggagcgtaga cggtgtggca agtctgatgt 240
gaaaggcatg ggctcaacct gtggactgca ttggaaactg tcatacttga gtgccggagg 300
ggtaagcgga attcctagtg tagcggtgaa atgcgtagat attaggagga acaccagtgg 360
cgaaggcggc ttactggacg gtaactgacg ttgaggctcg aaagcgtggg gagcaaacag 420
gattagatac ccgagtagtc 440
<210> 4
<211> 465
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 4
cctacgggtg gcagcagtgg ggaatattgg acaatggggg gaaccctgat ccagccatgc 60
cgcgtgtgtg aagaaggccc tcgggttgta aagcactttc agtgaggaag aacgcctggt 120
ggttaatacc catcaggaaa gacatcactc acagaagaag caccggctaa ctccgtgcca 180
gcagccgcgg taatacggag ggtgcgagcg ttaatcggaa ttactgggcg taaagcgcgc 240
gtaggtggct tgataagccg gttgtgaaag ccccgggctc aacctgggaa cggcatccgg 300
aactgtcaag ctagagtgca ggagaggaag gtagaattcc cggtgtagcg gtgaaatgcg 360
tagagatcgg gaggaatacc agtggcgaag gcggccttct ggactgacac tgacactgag 420
gtgcgaaagc gtgggtagca aacaggatta gataccctag tagtc 465
<210> 5
<211> 440
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 5
cctacgggtg gcagcagtgg ggaatattgg acaatgggcg caagcctgat ccagccatgc 60
cgcgtgagtg atgaaggcct tagggttgta aagctctttc accggagaag ataatgacgg 120
tatccggaga agaagccccg gctaacttcg tgccagcagc cgcggtaata cgaagggggc 180
tagcgttgtt cggaattact gggcgtaaag cgcacgtagg cggactatta agtcaggggt 240
gaaatcccgg agctcaactc cggaactgcc tttgatactg gtagtctaga gaccgagaga 300
ggtaagtgga actcctagtg tagaggtgga attcgtagat attaggaaga acaccagtgg 360
cgaaggcggc ttactggctc ggaactgacg ctgaggtgcg aaagcgtggg gagcaaacag 420
gattagatac ccctgtagtc 440
<210> 6
<211> 442
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 6
cctacgggtg gcagcagtgg ggaatattgc acaatggggg aaaccctgat gcagcgacgc 60
cgcgtgagcg aagaagtatt tcggtatgta aagctctatc agcagggaag aagaaatgac 120
ggtacctgac taagaagcac cggctaaata cgtgccagca gccgcggtaa aacgtaggtc 180
acaagcgttg tccggaatta ctgggtgtaa agggagcgca ggcgggagaa caagttggaa 240
gtgaaatcca tgggctcaac ccatgaactg ctttcaaaac tgtttttctt gagtagtgca 300
gaggtaggcg gaattcccgg tgtagcggtg gaatgcgtag atatcgggag gaacaccagt 360
ggcgaaggcg gcctactggg caccaactga cgctgaggct cgaaagtgtg ggtagcaaac 420
aggattagat acccctgtag tc 442
<210> 7
<211> 440
<212> DNA
<213> 肠道细菌(Microorganism)
<400> 7
cctacgggcg gcagcagtgg ggaatattgc acaatggggg aaaccctgat gcagcgacgc 60
cgcgtgagcg aagaagtatt tcggtatgta aagctctatc agcagggaag aaaatgacgg 120
tacctgacta agaagcaccg gctaaatacg tgccagcagc cgcggtaata cgtatggtgc 180
aagcgttatc cggatttact gggtgtaaag ggagcgtaga cggcaaggca agtctgatgt 240
gaaaacccag ggcttaaccc tgggactgca ttggaaactg tctggctcga gtgccggaga 300
ggtaagcgga attcctagtg tagcggtgaa atgcgtagat attaggagga acaccagtgg 360
cgaaggcggc ttactggacg gtaactgacg ttgaggctcg aaagcgtggg gagcaaacag 420
gattagatac cctagtagtc 440
<210> 8
<211> 19
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 8
actcctacgg gaggcagca 19
<210> 9
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 9
ggactachvg ggtwtctaat 20

Claims (10)

1.一种用于区别新型冠状病毒感染和非感染者的肠道微生物基因标志物，其特征在于：由SEQ ID NO:1-7所示的7种微生物基因组成，所述微生物在肠道中富集。
2.一种用于检测权利要求1所述肠道微生物基因标志物的检测试剂，包括用于检测权利要求1所述的SEQ ID NO:1-7所示的7种微生物基因的引物，所述引物序列为SEQ ID NO:8-9。
3.权利要求2所述检测试剂在制备新型冠状病毒感染和非感染的区分检测试剂盒中的应用，所述检测试剂适用于检测权利要求1所述的肠道微生物基因。
4.根据权利要求3所述的应用，其特征在于：新型冠状病毒感染者包括RT-PCR检测阳性的确诊患者及RT-PCR检测阴性但IgG抗体阳性的临床诊断患者。
5.根据权利要求3所述的应用，其特征在于：对所述对象的粪便进行检测，以便确定该样本是否包含所述的微生物基因，是否可以建立区别新型冠状病毒感染和非感染的微生物基因模型。
6.根据权利要求5所述的应用，其特征在于：通过收集入组对象的粪便样本，抽提微生物总DNA，完成微生物DNA的16S rDNA Miseq测序，检测是否存在权利要求1所述的7种微生物基因。
7.根据权利要求5所述的应用，其特征在于：通过收集入组对象的粪便样本，抽提微生物总DNA，进行肠道菌群的16S rDNA Miseq测序；基于高通量测序数据，在训练集中建立新型冠状病毒感染和非感染者的微生物区别模型，建立新冠肺炎患病率(probability ofdisease，POD)指数；POD指数在验证集中计算其区别能力，进行验证；进一步在新冠肺炎疑似患者中进行验证，以证明微生物基因区别模型弥补RT-PCR的可能性。
8.一种用于区别新型冠状病毒感染和非感染者的肠道微生物模型的试剂盒，包括用于检测权利要求1所述的SEQ ID NO:1-7所示的7种微生物基因的引物。
9.一种用于区别新型冠状病毒感染和非感染者的肠道微生物基因区别模型，所述肠道微生物模型包含权利要求1所述的肠道微生物基因标志物，适用于区分新型冠状病毒感染和非感染者，其中新型冠状病毒感染者包括RT-PCR检测阳性的确诊患者及RT-PCR检测阴性但IgG抗体阳性的临床诊断患者。
10.权利要求9所述的微生物基因区别模型的建构方法，包括如下步骤：
(1)收集入组对象的粪便样本，入组对象包括24例新冠肺炎确诊患者和48例健康对照，按照DNA的标准抽提方法完成粪便样本中微生物总DNA的抽提，在Illumina MiSeq平台完成肠道菌群的16S rDNA高通量测序工作；
(2)基于高通量测序数据，在微生物区别模型的训练集中，在24例新冠肺炎确诊患者和48例健康对照者之间，基于一个随机森林模型，通过一个五倍交叉验证的算法，鉴定了用于该模型的最佳的7个微生物基因标志物；
(3)基于7个微生物基因标志物，通过使用随机生成的决策树的比率来计算新冠肺炎的患病率(Probability of disease,POD)指数；
(4)该微生物区别模型在24例新冠肺炎确诊患者和48例健康对照者之间的区别能力达到99.74％，POD指数在新冠肺炎确诊患者中明显升高，两组之间有显著性差异(p＝1.2x10-13)；
(5)在验证集中，该微生物区别模型在12例新冠肺炎确诊患者和24例健康对照之间的区别能力达到99.31％；POD指数在早期肝癌患者中明显升高，两组之间有显著性差异(p＝2.0x10-6)；
(6)在验证集中，POD指数在23例中新冠疑似患者中明显升高，两组之间有显著性差异(p＝1.8x10-8)，该微生物区别模型在在23例新冠疑似患者和24例健康对照之间的区别能力达到98.01％。
CN202011628621.6A 2020-12-31 2020-12-31 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用 Active CN112609024B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011628621.6A CN112609024B (zh) 2020-12-31 2020-12-31 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011628621.6A CN112609024B (zh) 2020-12-31 2020-12-31 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用

Publications (2)

Publication Number Publication Date
CN112609024A true CN112609024A (zh) 2021-04-06
CN112609024B CN112609024B (zh) 2023-08-15

Family

ID=75253130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011628621.6A Active CN112609024B (zh) 2020-12-31 2020-12-31 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用

Country Status (1)

Country Link
CN (1) CN112609024B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109371147A (zh) * 2018-11-01 2019-02-22 任志刚 一种用于区别肝癌和非肝癌的肠道微生物基因标志物及其应用
CN111041089A (zh) * 2020-03-13 2020-04-21 广州微远基因科技有限公司 Covid-19感染的宿主标志物的应用
CN111378788A (zh) * 2020-04-08 2020-07-07 广州微远基因科技有限公司 辅助covid-19诊断的菌种标志物及其应用

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109371147A (zh) * 2018-11-01 2019-02-22 任志刚 一种用于区别肝癌和非肝癌的肠道微生物基因标志物及其应用
CN111041089A (zh) * 2020-03-13 2020-04-21 广州微远基因科技有限公司 Covid-19感染的宿主标志物的应用
CN111378788A (zh) * 2020-04-08 2020-07-07 广州微远基因科技有限公司 辅助covid-19诊断的菌种标志物及其应用

Also Published As

Publication number Publication date
CN112609024B (zh) 2023-08-15

Similar Documents

Publication Publication Date Title
CN109852714B (zh) 一种肠癌早期诊断和腺瘤诊断标志物及用途
CN105296590B (zh) 大肠癌标志物及其应用
CN111254207A (zh) 用于区别自身免疫性肝炎和健康人群的肠道微生物标志物及其应用
CN112899368B (zh) 一种用于原发性肝细胞癌早期诊断的生物标志物及其检测试剂和应用
CN109371147A (zh) 一种用于区别肝癌和非肝癌的肠道微生物基因标志物及其应用
CN110838365A (zh) 肠易激综合症相关菌群标志物及其试剂盒
CN107058521A (zh) 一种检测人体免疫状态的检测系统
CN115873956A (zh) 用于预测受试者患有结直肠癌风险的试剂盒、系统、应用和预测模型的建模方法
CN112680549A (zh) 一种用于区别自身免疫性肝炎和健康人群的口腔微生物基因标志物及其应用
CN106156541A (zh) 分析个体两类状态的免疫差异的方法和装置
KR101940446B1 (ko) 미생물 메타게놈 분석을 통한 난소암 진단방법
CN113913490A (zh) 非酒精性脂肪肝标志微生物及其应用
CN115927588A (zh) 一种用于新型冠状病毒感染患者无创诊断的口腔真菌微生物基因标志物及其应用
CN111630189B (zh) 通过细菌性宏基因组分析来诊断胆管癌的方法
CN112609024A (zh) 一种用于新型冠状病毒感染患者无创诊断的肠道微生物基因标志物及其应用
CN114317725B (zh) 克罗恩病生物标志物、试剂盒及生物标志物的筛选方法
CN113308560B (zh) 一种用于新型冠状病毒感染患者无创诊断的口腔微生物基因标志物及其应用
CN114214438B (zh) 胆道菌群检测试剂在制备预测胆道结石术后早期复发的试剂中的应用
CN114480636B (zh) 胆汁细菌作为肝门部胆管癌诊断及预后标志物的用途
CN113862382B (zh) 肠道菌群的生物标志物在制备诊断成人免疫性血小板减少症的产品中的应用
CN117106946A (zh) 用于新型冠状病毒感染患者无创诊断的口腔真菌微生物基因标志物及其应用
CN116445599A (zh) 一种用于系统性红斑狼疮患者无创诊断的口腔微生物基因标志物及其临床应用
CN114717342B (zh) 用于预测新冠肺炎患者一年后中和抗体水平的肠道微生物基因标志物及其应用
CN114703305B (zh) 用于预测新冠肺炎患者一年后中和抗体水平的口腔微生物基因标志物及其应用
WO2016049927A1 (en) Biomarkers for obesity related diseases

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant