CN114262743B - 中风标志微生物及其应用 - Google Patents

中风标志微生物及其应用 Download PDF

Info

Publication number
CN114262743B
CN114262743B CN202111672791.9A CN202111672791A CN114262743B CN 114262743 B CN114262743 B CN 114262743B CN 202111672791 A CN202111672791 A CN 202111672791A CN 114262743 B CN114262743 B CN 114262743B
Authority
CN
China
Prior art keywords
genus
bacteria
clostridium
species
microorganisms
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.)
Active
Application number
CN202111672791.9A
Other languages
English (en)
Other versions
CN114262743A (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.)
Qingdao Ruiyi Precision Medical Test Co ltd
Original Assignee
Qingdao Ruiyi Precision Medical Test 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 Qingdao Ruiyi Precision Medical Test Co ltd filed Critical Qingdao Ruiyi Precision Medical Test Co ltd
Priority to CN202111672791.9A priority Critical patent/CN114262743B/zh
Publication of CN114262743A publication Critical patent/CN114262743A/zh
Application granted granted Critical
Publication of CN114262743B publication Critical patent/CN114262743B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明提出了中风标志微生物及其应用,该中风标志微生物包括第一微生物集,因此,进一步提出了一种试剂盒,包括适于检测第一微生物集中的至少一种菌种的试剂,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。本发明提出的标志微生物在健康人群和中风患者中的丰度具有显著差异,可以作为检测和/或治疗中风的标志物。

Description

中风标志微生物及其应用
技术领域
本发明涉及生物技术领域,具体的,本发明涉及中风标志微生物及其用途,更具体的,本申请涉及一种试剂盒、试剂在制备试剂盒中的用途、用于预防或者治疗中风的药物组合物、确定个体是否患有中风方法、确定个体是否患有中风的装置、一种装置、一种筛选药物的方法。
背景技术
中风(Apoplexy),又称为脑卒中,是在气血失调、阴阳失衡的基础上遇劳倦内伤、忧思恼怒、气候变化、嗜食烟酒或病久失治等诱因,进而引起脏腑阴阳失调,气血逆乱直冲犯脑,导致脑泳痹阻或血溢脑脉之外而发病,临床以半身不遂、口眼喎斜、神志昏蒙、舌强言謇或不语、偏身麻木为主症。该病具有起病急、变化快、病情危重易损害神志的特点,多见于体衰的中老年人。据现代流行病学研究调查所得,自20世纪80年代以来,我国多数地区脑卒中发病率呈上升趋势,我国平均发病年龄约为60.9岁,脑卒中存活患者中有70%~85%不同程度地丧失活动能力和工作能力,严重影响生活质量。脑血管病已成为城乡居民疾病死亡原因的前列,给社会及家庭带来沉重的负担。因此,脑卒中的预防、治疗及康复已成为国内外医家的研究热点,积极探索中风相关标志物对治疗及预防中风具有重要意义。
发明内容
本申请是基于发明人对以下事实和问题的发现和认识作出的:
本申请的申请人通过前期大量研究,本发明意外地发现一种微生物组合可以作为检测中风的微生物物种标志物,为早期发现中风提供一种非侵入性方法;合理有效地应用微生物物种标志物,扶持肠道有益菌的生长,抑制肠道潜在致病菌,可以治疗或减轻中风的临床症状。
为此,在本发明的第一方面,本发明提出了一种试剂盒。根据本发明的实施例,包括适于检测第一微生物集中至少之一的菌种的试剂,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。根据本发明具体实施例的试剂盒,包含检测第一微生物集中的至少一种菌种的试剂,因此,可以准确检测生物样品中是否含有第一微生物集中的至少一种菌种,更进一步地,也可以准确区分或诊断中风患者和健康个体。
在本发明的第二方面,本发明提出了试剂在制备试剂盒中的用途,所述试剂适于检测第一微生物集中的至少一种菌种。根据本发明的实施例,所述试剂盒用于诊断中风或检测中风的治疗效果,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。根据本发明具体实施例的试剂制备的试剂盒,可以准确的检测生物样品中是否包含所述第一微生物集中的至少一种菌种,更进一步地,也可以极准确的区分中风患者和健康个体,由此,可以有效的在早期进行中风诊断,或用于检测治疗过程中中风的变化。
在本发明的第三方面,本发明提出了一种用于预防或者治疗中风的药物组合物。根据本发明的实施例,含有第一微生物集中的至少一种菌种,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。根据本发明实施例的中风标志微生物中第一微生物集的菌种可以非侵入性的在早期发现或辅助检测中风,确定个体患有中风的概率高低或者个体处于健康状态的概率高低;同时,提高中风高风险人群或已中风患者肠道内的所述第一微生物集中的各种菌种,可以降低患中风的概率或减缓、治愈中风,因此,所述包含所述第一微生物集中的至少一种菌种的药物能够用于平衡肠道菌群,有效预防或治疗中风。
在本发明的第四方面,本发明提出了一种确定个体是否患有中风的方法。根据本发明的实施例,包括:(1)确定所述个体的粪便样本中标志微生物的丰度,所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种;(2)将步骤(1)中得到的所述丰度与预定的阈值进行比较,以便确定所述个体是否患有中风;其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。根据本发明实施例的方法可以依据个体的粪便样本中的所述标志微生物的丰度确定个体是否患有中风,所述标志微生物是发明人对大量已知状态的粪便样本进行验证,通过差异比较分析各种肠道微生物在中风组和健康组粪便样本中的丰度,而确定下来的。
在本发明的第五方面,本发明提出了一种确定个体是否患有中风的装置。根据本发明的实施例,包括:丰度确定单元,用于确定所述个体的粪便样本中标志微生物的丰度,所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种;比较单元,用于将所得到的所述丰度与预定的阈值进行比较,以便确定所述个体是否患有中风;其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridiumsensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。所述标志微生物是发明人通过差异比较分析各种肠道微生物在中风患者和健康人群的粪便样本中的丰度并经过分析和大量已知状态的粪便样本的验证而确定下来的,根据本发明实施例的装置可以准确确定个体是否为中风的高风险人群或中风患者。
在本发明的第六方面,本发明提出了一种装置。根据本发明的实施例,包括:计算机可读存储介质,其上存储有计算机程序,所述程序用于执行第四方面所述的方法;以及一个或者多个处理器,用于执行所述计算机可读存储介质中的程序。根据本发明实施例的装置可以准确确定个体是否为中风的高风险人群或中风患者。
在本发明的第七方面,本发明提出了一种筛选药物的方法。根据本发明的实施例,所述药物用于治疗或者预防中风,所述方法包括:将候选药物施用于受试者,检测施用前后,所述受试者粪便中标志微生物的丰度,所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种,其中,满足下列条件至少之一的候选药物适于用于治疗或者预防中风:(1)进行所述施用后,所述第一微生物集中的至少一种菌种的所述丰度升高;和(2)进行所述施用后,所述第二微生物集中的至少一种菌种的所述丰度降低;其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensustricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。根据本发明实施例的方法可以生产或筛选出促进所述标志微生物中第一微生物集中各种菌种生长和/或抑制肠道标志微生物中第二微生物集中的各种菌种生长的药物,对于辅助减轻中风的临床症状具有重要意义。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是根据本发明实施例的筛选中风标志微生物的实验分析流程示意图;以及
图2是根据本发明实施例的23个标志微生物综合指标AUC评价结果示意图,其中,Specificity表示特异度,即预测为阳性且实际为阳性,真阳性,纵坐标Sensitivity表示敏感度,即真阴性,Confidence interval表示置信区间:
2-A为一期163个样品数据ROC曲线下AUC值和置信区间结果图;
2-B为二期70个样品数据ROC曲线下AUC值和置信区间结果图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
术语“任选地”仅用于描述目的,而不能理解为指示或暗示相对重要性。由此,限定有“任选地”的特征可以明示或者隐含地包括或不包括该特征。
生物学标志物是从生物学介质中可以检测到的细胞/生物化学或分子改变。生物学介质包括各种体液、组织、细胞、粪便、头发、呼气等。
微生物的丰度是指在某一微生物群体中该种微生物的丰富程度,例如在肠道微生物群体中的该种微生物程度,可表示为该种微生物在该群体中的含量。
根据本发明提供的一种试剂盒,包括适于检测第一微生物集中的至少一种菌种的试剂,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。
根据本发明的具体实施方案,所述试剂盒进一步包括适于检测第二微生物集中的至少一种菌种的试剂,所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。
根据本发明的具体实施方案,所述试剂盒包括适于检测所述第一微生物集中全部所述菌种的试剂。
根据本发明的具体实施方案,所述试剂盒包括适于检测所述第二微生物集中全部所述菌种的试剂。
根据本发明的具体实施方案,所述标志微生物是发明人通过对大量患中风个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证,而确定下来的,明确了肠道微生物中中风相关的标志微生物。利用包含检测所述标志微生物的试剂的试剂盒能够准确检测生物样品中是否含有所述微生物菌种,更进一步地,也可以准确确定个体处于患有中风状态的概率高低或者处于健康状态的概率高低,能够用于非侵入性的早期发现或辅助检测中风。
根据本发明的具体实施方案,所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制,任何可以检测所述微生物菌种的试剂均包含在本发明的范围内,如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂,具体的,如抗体、酶、核酸分子等。
本文中,所述微生物形态学特征指:在显微镜下观察微生物外形大小、形状、排列等,细胞构造,革兰氏染色反应,能否运动、鞭毛着生部位和数目,有无芽孢和荚膜、芽孢的大小和位置,放线菌和真菌的繁殖器官的形状、构造,孢子的数目、形状、大小、颜色和表面特征等。
本文中,所述微生物生理生化反应特征指:微生物利用物质的能力、代谢产物的特殊性,如是否产生H2S、吲哚、CO2、醇、有机酸,能否还原硝酸盐,能否使牛奶凝固、冻化等、生长环境(温度、湿度、氧气及二氧化碳等气体的浓度、PH、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如:共生、寄生、宿主范围及致病情况)等。
本文中,所述微生物血清学反应指:利用抗原与抗体的高度敏感特异性反应,来鉴别相似的菌种,或对同种微生物分型,如用已知菌种、型或菌株制成的抗血清,与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
本文中,分子生物学方法检测微生物主要包括:利用PCR技术、高通量测序等方法。
根据本发明提供的试剂在制备试剂盒中的用途,所述试剂适于检测第一微生物集中的至少一种菌种,所述试剂盒用于诊断中风或者检测中风的治疗效果,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensustricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。
根据本发明的具体实施例,所述标志微生物是发明人通过对大量患中风个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证,而确定下来的,明确了肠道微生物中中风相关的微生物标志物。利用检测所述标志微生物的试剂能够确定个体患有中风的概率高低或者处于健康状态的概率高低,能够用于非侵入性的早期发现或辅助检测中风。
根据本发明一些具体的实施例,所述试剂进一步适于检测第二微生物集中的至少一种菌种,所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。
根据本发明的具体实施例,所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制,可以检测所述微生物菌种的试剂均包含在本发明的范围内,如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂等,具体的,如抗体、酶、核酸分子。
本文中,所述微生物形态学特征指:在显微镜下观察微生物外形大小、形状、排列等,细胞构造,革兰氏染色反应,能否运动、鞭毛着生部位和数目,有无芽孢和荚膜、芽孢的大小和位置,放线菌和真菌的繁殖器官的形状、构造,孢子的数目、形状、大小、颜色和表面特征等。
本文中,所述微生物生理生化反应特征指:微生物利用物质的能力、代谢产物的特殊性,如是否产生H2S、吲哚、CO2、醇、有机酸,能否还原硝酸盐,能否使牛奶凝固、冻化等、生长环境(适宜生长的温度、湿度、氧气及二氧化碳等气体的浓度、PH、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如:共生、寄生、宿主范围及致病情况)等。
本文中,所述微生物血清学反应指:利用抗原与抗体的高度敏感特异性反应,来鉴别相似的菌种,或对同种微生物分型,如用已知菌种、型或菌株制成的抗血清,与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
本文中,分子生物学方法检测微生物主要包括:利用PCR技术、高通量测序等方法。
根据本发明提供的一种用于预防或者治疗中风的药物组合物,含有第一微生物集中的至少一种菌种,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属。
上述标志微生物是发明人通过差异比较分析各种肠道微生物在中风疾病组和健康组的粪便样本中的丰度,以及经过大量已知状态的粪便样本的验证,而确定下来的。所述标志微生物中的菌种相较于中风患者群体,在健康群体组中显著富集,所述显著富集是指与在中风患者组中的丰度相比,上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在中风患者组中的丰度;能够使该部分菌种丰度提高的物质能够用于治疗中风或者益于中风患者服用,能够使其丰度提高的物质不限于治疗中风的药物。该实施方式提供的标志微生物能够用于制备治疗中风的药物,所述药物可有效治疗或缓解中风。
根据本发明提供的一种确定个体是否患有中风的方法,包括步骤(1)和(2)。
(1)确定所述个体的粪便样本中的标志微生物的丰度。
所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种。其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridiumsensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属、Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。
根据本发明的一些具体的实施例,步骤(1)进一步包括:获得所述个体的粪便样本中的核酸测序数据;将所述测序数据与参考基因组进行比对;基于所述比对的结果,确定所述标志微生物的丰度。
根据本发明的具体实施例,在步骤(1)中,按照下列公式确定所述标志微生物的丰度:Ab(S)=Ab(US)+Ab(MS),其中,S表示所述标志微生物的编号,Ab(S)表示所述标志微生物S的丰度,Ab(US)=US/lS,US为所述测序数据中与所述标志微生物S的参考基因组唯一比对的读段数目,lS为所述标志微生物S的参考基因组的总长度,MS为所述测序数据中与所述标志微生物S的参考基因组非唯一比对的读段的数目,i表示所述非唯一比对读段的编号,Coi为所述第i读段对应的丰度系数,/>Coi,s表示针对所述标志微生物S,所述非唯一比对的读段i的丰度系数,N为所述非唯一比对的读段i能够比对的微生物的总数,j表示所述非唯一比对的读段i能够比对的微生物的编号。
比对可以利用已知比对软件进行,例如SOAP、BWA和TeraMap等,在比对过程中,一般对比对参数进行设置,设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch),例如设置s≤2,若reads中有超过s个碱基发生错配,则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况,包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。
所述菌种/微生物的参考基因组指预先确定的该微生物物种的序列,可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板,例如,目标是待测样本中的微生物,参考序列可选择NCBI数据库中的各种微生物的参考基因组和/或HMP、MetaHIT项目公开的DACC肠道参考基因组,进一步地,也可以预先配置包含更多参考序列的资源库,例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例,各种微生物的参考基因组从公开的数据库中获得,通常的,一种微生物的参考基因组有多个版本,即一种微生物有多个公开的参考基因组。
reads与物种的参考基因组比对,比对上的可以被分为两部分:a)Unique reads(U):唯一比对上一个物种的参考基因组;称这些reads为unique reads。即,如果reads比对上的参考基因组均来自同一物种,定义这些reads为unique reads;b)Multiple reads(M):比对上一个以上物种的参考基因组,定义为multiple reads。即,如果reads比对上的参考基因组来自至少两种物种,定义这些reads为multiple reads。
(2)丰度比较,以确定个体是否患有中风。
根据本发明的一个实施例,将步骤(1)中得到的所述丰度与预定的阈值进行比较,以便确定所述个体是否患有中风。
根据本发明的一些具体实施例,所述阈值为预先设定的,包括患中风丰度阈值和不患中风丰度阈值。将标志微生物在待测个体样本中的丰度与所述的阈值比较,进行待测个体的状态的确定。所述阈值可以为一数值或者数值范围,例如,基于已知患病或健康状态个体中的标志微生物的丰度均值,该微生物对应的阈值可以设为该丰度均值的95%的置信区间。
所述的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中,一个概率样本的置信区间(Confidence interval)是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度,即前面所要求的“一定概率”,这个概率被称为置信水平。
根据本发明的一些具体实施例,当步骤(1)中确定的标志微生物的丰度达到所述患中风丰度阈值,未达到所述不患中风丰度阈值时,确定所述个体患有中风,当(1)中确定的标志微生物的丰度达到不患中风丰度阈值,未达到患中风丰度阈值时,确定所述个体不患中风。
需要说明的是,根据目的或要求不同,可能对确定个体状态结果的可信程度有不同的要求,本领域技术人员可以选择不同的显著性水平或阈值。
该方法基于检测个体的粪便样本中的标志微生物中的各种菌种的丰度,分别将检测确定的各种菌种的丰度与其阈值进行比较,依据获得的比较结果能够确定个体为中风个体或者为健康个体的概率。为早期发现中风提供一种非侵入性的辅助检测或者辅助干预治疗的方法。
以上任一实施例中的利用标志微生物确定个体是否患有中风的方法的全部或部分步骤,可以利用包含可拆分的相应单元功能模块的装置/系统来施行,或者将方法程序化、存储于机器可读介质,利用机器运行该可读介质来实现。
根据本发明提供的一种确定个体是否患有中风的装置,该装置包括:丰度确定单元,用于确定所述个体的粪便样本中标志微生物的丰度,所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种;比较单元,用于将所得到的所述丰度与预定的阈值进行比较,以便确定所述个体是否患有中风;其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。上述对本发明任一实施例的利用标志微生物确定个体是否患有中风的方法的技术特征和优点的描述,同样适用本发明这一方面的装置,在此不再赘述。
根据本发明的实施例,所述丰度确定单元适于通过下列步骤确定所述丰度:获得所述个体的粪便样本中的核酸测序数据;将所述测序数据与参考基因组进行比对;基于所述比对的结果,确定所述标志微生物的丰度。
比对可以利用已知比对软件进行,例如SOAP、BWA和TeraMap等,在比对过程中,一般对比对参数进行设置,设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch),例如设置s≤2,若reads中有超过s个碱基发生错配,则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况,包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例,利用SOAPalign 2.21进行比对,设置参数为-r 2-m 100-x 1000。
所述微生物的参考基因组指预先确定的该微生物物种的序列,可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板,例如,目标是待测样本中的微生物,参考序列可选择NCBI数据库中的各种微生物的参考基因组和/或HMP、MetaHIT项目公开的DACC肠道基因组,进一步地,也可以预先配置包含更多参考序列的资源库,例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例,各种微生物的参考基因组从公开的数据库中获得,通常的,一种微生物的参考基因组有多个版本,即一种微生物有多个公开的参考基因组。
reads与物种的参考基因组比对,比对上的可以被分为两部分:a)Unique reads(U):唯一比对上一个物种的基因组;称这些reads为unique reads。即,如果reads比对上的参考基因组均来自同一物种,定义这些reads为unique reads。b)Multiple reads(M):比对上一个以上物种的参考基因组,定义为multiple reads。即,如果reads比对上的参考基因组来自至少两种物种,定义这些reads为multiple reads。
根据本发明的一个实施例,按照下列公式确定所述标志微生物的丰度:Ab(S)=Ab(US)+Ab(MS),其中其中,S表示所述标志微生物的编号,Ab(S)表示所述标志微生物S的丰度,Ab(US)=US/lS,US为所述测序数据中与所述标志微生物S的参考基因组唯一比对的读段数目,lS为所述标志微生物S的参考基因组的总长度,MS为所述测序数据中与所述标志微生物S的参考基因组非唯一比对的读段的数目,i表示所述非唯一比对读段的编号,Coi为所述第i读段对应的丰度系数,/>Coi,s表示针对所述标志微生物S,所述非唯一比对的读段i的丰度系数,N为所述非唯一比对的读段i能够比对的微生物的总数,j表示所述非唯一比对的读段i能够比对的微生物的编号。上述对本发明任一实施例的利用标志微生物确定个体是否患有中风的方法的技术特征和优点的描述,同样适用本发明这一方面的装置,在此不再赘述。
根据本发明的又一个实施方式提供的一种装置,包括:计算机可读存储介质,其上存储有计算机程序,所述程序用于执行前面所述的一种确定个体是否患有中风方法;以及一个或者多个处理器,用于执行所述计算机可读存储介质中的程序。
根据本发明的又一个实施方式提供的一种筛选药物的方法,所述药物用于治疗或者预防中风,所述方法包括:将候选药物施用于受试者,检测施用前后,所述受试者粪便中标志微生物的丰度,所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种,其中,满足下列条件至少之一的候选药物适于用于治疗或者预防中风:(1)进行所述施用后,所述第一微生物集中的至少一种菌种的所述丰度升高;和(2)进行所述施用后,所述第二微生物集中的至少一种菌种的所述丰度降低;其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)、Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属、Butyricicoccus菌属;所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XlVb、巨球菌属(Megasphaera)、大肠杆菌属(Escherichia/Shigella)、卟啉单胞菌属(Porphyromonas)、Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属、Alistipes菌属。
利用本发明这一方面的生产或筛选治疗中风的药物的方法,通过合理有效地应用确定的中风生物标志物进行筛选,能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物,对于辅助减轻中风的临床症状具有重要意义。
下面将对实施例作具体介绍。以下实施例中涉及的未特别交待的试剂、序列、软件及仪器,都是常规市售产品。
实施例1生物标志物的鉴定
1、样本收集和DNA提取
中风患者来自南方医科大学,实验人员共采集141例中风患者和92例健康对照的粪便样品,其中每个个体的新鲜粪便样品分成200mg/份,共5份,立即-80℃冰箱冷冻保存。采用苯酚三氯甲烷处理提取DNA方法提取141例中风患者和92例健康对照的粪便样品总DNA。
2、构建DNA文库和测序
DNA建库按仪器制造商(Illumina)的操作指南进行。对文库进行PE150 bp测序。Illumina HiSeq2000(Illumina,San Diego,CA)平台对141例中风患者和92例健康对照的文库16S rDNA高变区即V4区的序列进行测序。每个样本平均产生6.05Mb(sd.±1.18Mb)高质量测序结果,总计1.411Gb测序数据量。其中99例中风患者和64例健康对照样品作为一期,42例中风患者和28例健康对照样品作为二期。
参照图1的实验流程,鉴定中风的相关生物标志物,其中省略的步骤或者细节为本领域技术人员所熟知,几个重要步骤介绍如下面所述。
3、微生物物种丰度分析
3.1序列优化统计
所有样本的数据每个通过Pandaseq(V2.9)软件利用重叠关系将双末端测序得到的成对Reads拼接成一条序列,得到高变区的长Reads。然后使用内部撰写的程序对拼接后的Reads进行如下处理,以获取CleanReads:1)去除平均质量值低于20的Reads;2)去除Reads含N的碱基数超过3个的Reads;3)Reads长度范围为150~500nt,以统计Clean Reads的长度分布和数量。
3.2OTU分析
为便于下游物种多样性分析,将长Reads聚类为OTUs(Operational TaxonomicUnits)。首先把拼接的长Reads中的singletons(对应reads只有一条的序列)过滤掉,因为singletons可能由于测序错误造成,故将这部分序列去除,不加入聚类分析,利用Usearch(V7.0.1090)在0.97相似度下进行聚类,对聚类后的序列进行嵌合体过滤后,得到用于物种分类的OTU,每个OTU被认为可代表一个物种。
3.3物种注释和OTU丰度表分析
从各个OTU中挑选出一条序列,作为该OTU的代表序列。将该代表序列,用RDP方法,与已知物种的16S数据库(RDP,http://rdp.cme.msu.edu)进行比对,从而对每个OTU进行物种归类。归类后,根据每个OTU中序列的条数,从而得到OTU丰度表。这里得到的OTU丰度表就是16S测序技术对应的物种丰度表。
对于物种S,其丰度为Ab(S),与特有的U reads和共享的M reads相关,丰度的计算方式如下:
Ab(S)=Ab(US)+Ab(MS),
其中,S表示所述标志微生物的编号,
Ab(S)表示所述标志微生物S的丰度,
Ab(US)=US/lS
US为所述测序数据中与所述标志微生物S的参考基因组唯一比对的读段数目,
lS为所述标志微生物S的参考基因组的总长度,
MS为所述测序数据中与所述标志微生物S的参考基因组非唯一比对的读段的数目,
Coi为所述第i读段对应的丰度系数,
i表示所述非唯一比对读段的编号,
Coi,s表示针对所述标志微生物S,
所述非唯一比对的读段i的丰度系数,
N为所述非唯一比对的读段i能够比对的微生物的总数,
j表示所述非唯一比对的读段i能够比对的微生物的编号。
3.4物种丰度分析
在门、纲、目、科、属水平,将每个注释上的物种或OTU在不同样品中的序列数整理在一张表格中,形成profiling柱状图、星形图及统计表。
实施例2筛选微生物物种标记物
为了获得与中风疾病密切相关的肠道微生物物种标记物,发明人利用一期中风患者组(99例)与正常人组(64例)两组肠道微生物物种丰度数据,在物种级别做了一个与疾病相关性的研究。基于实施例1中得到的物种丰度表(即OTU丰度表),发明人设置标准如下:(1)通过结合Benjamini Hochberg的多重检验的Wilcoxon秩和检验进行检验,得到每个物种和中风疾病的相关性p值和q值;(2)使用一个具有统计意义的的阈值(p值<0.05),利用上述参数进行筛选。发明人得到23个与中风疾病密切相关的肠道微生物物种,数据分析如下所示:
99例中风患者和64例健康对照样品来自于实施例1所述的一期数据,发明人得到23个与中风疾病密切相关的肠道微生物物种,数据如表1所示,其中,在中风(stroke)患者肠道内中富集的微生物有12个物种,在正常的健康人群(HD)肠道中富集的微生物有11个物种。
表1:
实施例3微生物物种标记物的验证
为了证实实施例2中的发现,发明人参照实施例2所述的方法确定二期验证群中的28个健康人及42个中风患者的粪便样本中的表1所示23种菌属的丰度,并判断各样本种的这23种菌属丰度是否落入实施例2确定出的疾病组或健康组的95%的置信区间,判定23种菌种的丰度落入疾病组对应区间的样本所对应的个体的状态为中风患者,判定23种菌种的丰度落入健康组的对应区间的样本所对应的个体的状态为非中风患者。根据分析结果对表1中所示微生物物种标记物做出删选,验证群体的测序数据的获得与处理参照实施例2进行。
验证结果如下:验证富集在中风患者中的12个物种,1个在验证集找到高质量的验证(p值<0.05),中风患者富集的微生物物种标志物验证的p值和q值情况如表2所示。
表2:
分类 p值 q值 来源
Megasphaera 0.042122301 0.298639012 stroke
发明人将中风患者富集的1个微生物物种标记物作为中风疾病患病的正向指标,特别用于中风疾病非创伤式的检测和诊断。
物种标志物可以选自中风患者富集的物种标记物或者健康人群中富集的物种标记物中的一种或者多种。优选的,对于中风患者或易感人群,应该检测到表2中的物种标志物被富集。
治疗方案中,优选的,使表2中的物种标志物生长得到抑制或者清除,使表1中来自于健康个体的物种标记物被富集。
其他22个标记物虽然没有在验证过程中满足p值<0.05的条件,但是在以下ROC分类判断的一期数据和二期数据评价中,都显示出一定的应用价值。
发明人使用23个微生物物种标记物,构建一个综合指标,估计ROC(receiver-operating characteristic)曲线下面积AUC(AUC越大,表示诊断能力越高),评价综合得分对中风的诊断能力。通过一期163个样品和二期70个样品进行评测,在一期得到AUC=79.34%,二期得到AUC=63.99%,都表现出了一定的诊断能力,详细结果如图2a、图2b所示。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (2)

1.试剂在制备试剂盒中的用途,所述试剂适于检测第一微生物集和第二微生物集中的菌种,所述试剂盒用于诊断中风,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;
所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium)XIVb、巨球菌属(Megasphaera)Escherichia/Shigella、卟啉单胞菌属(Porphyromonas)Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属。
2.一种装置,其特征在于,包括:
计算机可读存储介质,其上存储有计算机程序,所述程序用于执行一种确定个体是否患有中风的方法,包括:
(1)确定待检测个体的粪便样本中标志微生物的丰度,所述标志微生物由第一微生物集和第二微生物集中的菌种组成;
(2)将步骤(1)中得到的所述丰度与预定的阈值进行比较,以便确定所述个体是否患有中风;
其中,所述第一微生物集由以下菌种组成:拟杆菌属(Bacteroides)、梭状芽胞杆菌属(Clostridium sensu stricto)、瘤胃球菌属(Ruminococcus)、柔嫩梭菌属(Faecalibacterium)、普氏菌属(Prevotella)Paraprevotella菌属、Gemmiger菌属、Fusicatenibacter菌属、Dorea菌属、Coprococcus菌属和Butyricicoccus菌属;
所述第二微生物集由以下菌种组成:梭菌属(Clostridium)XVIII、梭菌属(Clostridium) XIVb、巨球菌属(Megasphaera)Escherichia/Shigella、卟啉单胞菌属(Porphyromonas)Oscillibacter菌属、Intestinimonas菌属、Flavonifractor菌属、Erysipelotrichaceae incertae sedis菌属、Cloacibacillus菌属、Butyricimonas菌属和Alistipes菌属;
以及一个或者多个处理器,用于执行所述计算机可读存储介质中的程序。
CN202111672791.9A 2021-12-31 2021-12-31 中风标志微生物及其应用 Active CN114262743B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111672791.9A CN114262743B (zh) 2021-12-31 2021-12-31 中风标志微生物及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111672791.9A CN114262743B (zh) 2021-12-31 2021-12-31 中风标志微生物及其应用

Publications (2)

Publication Number Publication Date
CN114262743A CN114262743A (zh) 2022-04-01
CN114262743B true CN114262743B (zh) 2024-02-13

Family

ID=80832195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111672791.9A Active CN114262743B (zh) 2021-12-31 2021-12-31 中风标志微生物及其应用

Country Status (1)

Country Link
CN (1) CN114262743B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107217089A (zh) * 2016-03-22 2017-09-29 浙江中医药大学 确定个体状态的方法及装置
TWI698187B (zh) * 2019-05-28 2020-07-11 景岳生物科技股份有限公司 一種預防中風及改善中風嚴重度之組合物及其用途
CN111727263A (zh) * 2018-02-14 2020-09-29 Md保健株式会社 通过细菌性宏基因组分析来诊断中风的方法
CN112007051A (zh) * 2019-05-28 2020-12-01 景岳生物科技股份有限公司 一种预防中风及改善中风严重度的组合物及其用途
CN112011607A (zh) * 2020-09-19 2020-12-01 河北医科大学第二医院 源于肠道与脑卒中诊断和治疗效果相关的菌群及其应用
CN112251507A (zh) * 2020-09-19 2021-01-22 河北医科大学第二医院 一种脑卒中诊疗效果评价相关微生物
CN112553358A (zh) * 2020-09-19 2021-03-26 河北医科大学第二医院 微生物菌群在脑卒中中的应用

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019090353A1 (en) * 2017-11-06 2019-05-09 uBiome, Inc. Control processes for microorganism-related characterization processes

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107217089A (zh) * 2016-03-22 2017-09-29 浙江中医药大学 确定个体状态的方法及装置
CN111727263A (zh) * 2018-02-14 2020-09-29 Md保健株式会社 通过细菌性宏基因组分析来诊断中风的方法
TWI698187B (zh) * 2019-05-28 2020-07-11 景岳生物科技股份有限公司 一種預防中風及改善中風嚴重度之組合物及其用途
CN112007051A (zh) * 2019-05-28 2020-12-01 景岳生物科技股份有限公司 一种预防中风及改善中风严重度的组合物及其用途
CN112011607A (zh) * 2020-09-19 2020-12-01 河北医科大学第二医院 源于肠道与脑卒中诊断和治疗效果相关的菌群及其应用
CN112251507A (zh) * 2020-09-19 2021-01-22 河北医科大学第二医院 一种脑卒中诊疗效果评价相关微生物
CN112553358A (zh) * 2020-09-19 2021-03-26 河北医科大学第二医院 微生物菌群在脑卒中中的应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
脑卒中与肠道菌群及其代谢产物TMAO相关性的初步研究;廖硕希;中国优秀硕士学位论文全文数据库 医药卫生科技辑;摘要、第29页,图3 *

Also Published As

Publication number Publication date
CN114262743A (zh) 2022-04-01

Similar Documents

Publication Publication Date Title
CN107217089B (zh) 确定个体状态的方法及装置
CN107034279A (zh) 结核病微生物标志物在制备诊断结核病的试剂中的应用
CN107217088B (zh) 强直性脊柱炎微生物标志物
CN110541026A (zh) 一种检测溃疡性结肠炎的生物标志物及应用
CN114182007B (zh) 白塞病标志基因及其应用
CN113999922B (zh) 急性腹泻标志微生物及其应用
CN114262743B (zh) 中风标志微生物及其应用
CN114836508A (zh) 慢性阻塞性肺病标志微生物及其应用
CN107217086B (zh) 疾病标志物及应用
CN114317674B (zh) 类风湿性关节炎标志微生物及其应用
CN113913490B (zh) 非酒精性脂肪肝标志微生物及其应用
CN105671177B (zh) 强直性脊柱炎标志物及应用
CN113930479B (zh) 系统性红斑狼疮标志微生物及其应用
CN114381493A (zh) 炎症性肠病标志微生物及其应用
CN109161590A (zh) 整合素β4基因DNA甲基化位点在制备哮喘和或COPD早期诊断的生物标志物的应用
CN114317784B (zh) 白塞病标志微生物及其应用
CN116042866A (zh) 用于评估二型糖尿病患者粪菌移植疗效的微生物标志物及其应用
CN114381507B (zh) 格雷夫斯病标志微生物及其应用
CN114085886B (zh) 儿童克罗恩标志微生物及其应用
CN113999923B (zh) 小柳-原田综合症标志微生物及其应用
CN113913490A (zh) 非酒精性脂肪肝标志微生物及其应用
CN114891901A (zh) 溃疡性结肠炎标志微生物及其应用
CN114410809A (zh) 慢传输便秘标志微生物及其应用
CN114381492B (zh) 肌痛性脑脊髓炎标志微生物及其应用
CN114807392A (zh) 精神分裂症标志微生物及其应用

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