CN115101129A - 一种基于宏基因组测序数据组装病原微生物基因组的方法 - Google Patents
一种基于宏基因组测序数据组装病原微生物基因组的方法 Download PDFInfo
- Publication number
- CN115101129A CN115101129A CN202210732531.4A CN202210732531A CN115101129A CN 115101129 A CN115101129 A CN 115101129A CN 202210732531 A CN202210732531 A CN 202210732531A CN 115101129 A CN115101129 A CN 115101129A
- Authority
- CN
- China
- Prior art keywords
- genome
- sequence
- contig
- host
- sequences
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biophysics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明涉及一种基于宏基因组测序数据组装病原微生物基因组的方法。该方法将原始数据过滤后与宿主数据库进行比对,以达到在数据层面进行去除宿主序列的目的;再使用soapdenovo对去宿主的reads进行组装,得到无参组装的contig序列,将病源数据库作为参考基因组,然后统计去宿主的reads比对情况中每个位点的测序深度,得到有参组装的contig,将无参组装的contig序列和有参组装的contig序列进行整合,得到合并后的contig序列,并进行病原的判别。本发明通过去宿主后,采取无参组装和有参组装相结合的方法进行病原微生物的组装,得到的宏基因组没有宿主的污染,并且在点突变以外加入了结构变异的信息,准确度更高。
Description
技术领域
本发明涉及宏基因组测序与组装领域,具体涉及一种基于宏基因组测序数据组装病原微生物基因组的方法。
背景技术
微生物群体几乎存在于这个世界每一个生态群落之中,从个体体表到肠道,从高原空气到深海海底淤泥,从冰川冻湖到火山岩浆都无处不在,并扮演着不可或缺的角色。对微生物的研究从列文虎克发明显微镜开始的数百年中,主要基于纯培养的研究方式。在数以万亿计的微生物种类中,仅0.1%~1%的物种可培养,极大地限制了对微生物多样性资源的研究和开发。
Metagenomics(翻译成元基因组学,或者翻译成宏基因组学),是一种直接对微生物群体中包含的全部基因组信息进行研究的手段,即绕过对微生物个体进行分离培养,应用基因组学技术对自然环境中的微生物群落进行研究的学科。它规避了对样品中的微生物进行分离培养,提供了一种对不可分离培养的微生物进行研究的途径,更真实的反应样本中微生物组成、互作情况,同时在分子水平对其代谢通路、基因功能进行研究。
近年来,随着测序技术和信息技术的快速发展,利用新一代测序技术(NextGeneration Sequencing),即二代测序技术研究Metagenomics,能快速准确的得到大量生物数据和丰富的微生物研究信息,从而成为研究微生物多样性和群落特征的重要手段。如致力于研究微生物与人类疾病健康关系的人体微生物组计划(HMP,Human MicrobiomeProject),研究全球微生物组成和分布的全球微生物组计划(EMP,Earth MicrobiomeProject)都主要利用高通量测序技术进行研究。
宏基因组的二代测序数据常用的组装软件为soapdenovo,将病原微生物组装到contig水平。样本中由于可能存在大量宿主细胞或宿主的游离核酸,病原微生物核酸在提取出的核酸样本中的丰度通常极低。因此宏基因组测序数据组装之后的序列会含有宿主的contig序列,导致组装准确度降低,影响宏基因组的后续分析。
发明内容
本发明所要解决的技术问题就是提供一种基于宏基因组测序,通过对测序数据质控、与宿主数据库比对、去除宿主序列、进行有参和无参组装,进而得到较准确的宏基因组组装结果的方法,即提供一种基于宏基因组测序数据组装病原微生物基因组的方法。采用的技术方案为:
一种基于宏基因组测序数据组装病原微生物基因组的方法,具体步骤如下:
(1)使用fastp对二代下机数据进行质控,得到质控后高质量的二代数据;
(2)使用bwa将质控后高质量的二代数据比对到宿主数据库,得到比对情况;
(3)使用samtools软件,根据步骤(2)中比对情况提取未比对上的数据,得到去宿主的reads,即Clean reads;
(4)使用soapdenovo对去宿主的reads进行组装,得到无参组装的contig序列;
(5)使用bwa将去宿主的reads比对到病原数据库,得到去宿主的reads比对情况;
(6)使用samtools软件统计去宿主的reads比对情况中每个位点的测序深度,得到有参组装的contig序列;
(7)将无参组装的contig序列和有参组装的contig序列进行整合,得到合并后的contig序列;
(8)使用去宿主的Clean reads,以合并后的contig序列为参考进行比对,并重复步骤6,继续进行有参组装,得到最终版的contig序列;
(9)使用blast软件,将最终的contig序列与已有的病原微生物数据库进行比对,根据比对结果中序列之间的相似性与一致性等参数,对contig序列进行病原属种和亚型鉴定。
优选的,所述步骤(1)的具体方法为,根据实验信息与数据的实际情况,使用fastp软件对原始下机数据进行过滤,除去接头序列,低质量序列以及长度不符合要求的序列,得到符合质控要求的下机数据。
优选的,所述步骤(2)的具体方法为,质控后的下机数据中同时存在宿主的DNA序列以及病原微生物的DNA序列,使用bwa软件,以宿主物种的基因组作为参考基因组,将质控后的下机数据与参考基因组进行比对,获取下机数据中所有序列的比对情况,即all_bam文件。
优选的,所述步骤(3)的具体方法为,根据步骤(2)中得到的比对情况,使用samtools软件对all_bam文件中所有序列的比对情况进行过滤,对比对到参考基因组也就是宿主基因组的序列进行过滤,提取未比对上的所有序列,得到最终的去宿主的病原微生物的序列集合,即Clean reads。
优选的,所述步骤(4)的具体方法为,使用soapdenovo软件对去宿主的reads进行组装,得到不依赖病原微生物参考基因组的contig序列。
优选的,所述步骤(5)的具体方法为,使用bwa软件,以病原微生物的基因组作为参考基因组,将去宿主的reads比对到参考基因组上,获得去宿主的reads在病原微生物参考基因组上的比对情况,即pathogen_bam文件。
优选的,所述步骤(6)的具体方法为,使用samtools软件统计步骤(5)中pathogen_bam文件每个位点的测序深度,将每个位置上占比最高的碱基确定为该位置的碱基序列,并依次进行延伸;对于比对中出现的indel,根据位点两端序列以及位置信息,保留reads中的indel,得到有参组装的contig序列。
优选的,所述步骤(7)的具体方法为,使用自主编写脚本,以有参组装的contig序列为骨架,根据无参组装contig序列中的大片段易位或者倒位信息对骨架序列进行调整或重排,从而将无参组装的contig序列和有参组装的contig序列进行合并,得到合并后的contig序列。
优选的,所述步骤(8)的具体方法为,使用bwa软件,以步骤(7)中得到的合并的contig序列作为参考基因组,将Clean reads比对到参考基因组上,得到最终的比对情况,即final_bam文件,并参照步骤(6)中的方法,对得到的final_bam文件进行二次有参组装,得到最终版的contig序列。
与现有技术相比,本发明的有益之处在于:
本发明将原始数据过滤后与宿主数据库进行比对,以达到在数据层面进行去除宿主序列的目的;
本发明使用去宿主的序列分别进行无参组装和有参组装,并将无参和有参的结果进行合并,得到整合的基因组,然后将该整合基因组作为参考基因组,然后进行第二轮有参组装,通过去宿主后,采取无参组装和有参组装相结合的方法进行病原微生物的组装,组装的宏基因组没有宿主的污染,并且在点突变以外加入了结构变异的信息,准确度更高。
附图说明
图1为本发明一种基于宏基因组测序数据组装病原微生物基因组的方法流程图。
具体实施方式
下面结合附图、实施例、表格对本发明进行详细说明,附图仅用于示例性说明;应当理解,下面所提到的案例仅仅用来解释本发明,是为了便于描述本发明和简化描述,因此,不能理解为对本发明的限制。
如图1所示,一种基于宏基因组测序数据组装病原微生物基因组的方法,具体步骤如下:
步骤一.根据实验信息与数据的实际情况,使用fastp软件对含有甲型流感病毒的鸡胚尿囊液样本二代下机数据(Raw data)进行过滤,除去接头序列,低质量序列以及长度不符合要求的序列,得到符合质控要求的下机数据,即质控后高质量的二代数据。
步骤二.质控后的下机数据中同时存在宿主的DNA序列以及病原微生物的DNA序列,使用bwa软件将质控后高质量的二代数据比对到宿主(鸡)基因组数据库,以宿主(鸡)的基因组作为参考基因组,将质控后的下机数据与参考基因组进行比对,获取下机数据中所有序列比对情况的all_bam文件。
步骤三.根据步骤二中得到的比对情况,使用samtools软件对all_bam文件中所有序列的比对情况进行过滤,对比对到参考基因组(鸡基因组)的序列进行过滤,提取未比对上的所有序列,得到最终的去宿主的病原微生物的序列集合,即Clean reads。
步骤四.使用soapdenovo软件对Clean reads进行组装,得到不依赖病原微生物参考基因组的重叠群(contig)序列。
步骤五.使用bwa软件,以病原微生物的基因组作为参考基因组,将去宿主的Cleanreads比对到参考基因组上,获得Clean reads在病原微生物参考基因组上的比对情况,即pathogen_bam文件。
步骤六.使用samtools软件统计步骤五中pathogen_bam文件每个位点的测序深度,将每个位置上占比最高的碱基确定为该位置的碱基序列,并依次进行延伸;对于比对中出现的indel,根据位点两端序列以及位置信息,保留reads中的indel,得到有参组装的contig序列(即图1中的微生物基因组V1)。
步骤七.使用自主编写脚本,以有参组装的contig序列为骨架,根据无参组装contig序列中的大片段易位或者倒位信息对骨架序列进行调整或重排,从而将无参组装的contig序列和有参组装的contig序列进行合并,得到合并后的contig序列(即图1中的微生物基因组V2)。
步骤八.使用bwa软件,以步骤七中得到的合并的contig序列作为参考基因组,将Clean reads比对到参考基因组上,得到最终的比对情况(final_bam),并参照步骤六中的方法,对得到的final_bam文件进行二次有参组装,得到最终版的contig序列(即图1中的微生物基因组V3)。
步骤九.使用blast软件,将最终的contig序列与已有的病原微生物数据库进行比对,根据比对结果中序列之间的相似性与一致性等参数,对contig序列进行病原属种和亚型鉴定。
为了设置生物学重复,避免实验中可能出现的错误并且减小误差,随机选取宿主鸡的三个样本进行甲型流感病毒数据统计,统计结果如表1所示。
表1甲型流感病毒数据统计
其中,HQ:高质量reads数量;
Host:宿主(鸡)的reads数量;
IAV:流感病毒的reads数量;
Others:除去宿主和流感病毒后剩余的其他reads数量;
Host%:宿主reads数占高质量reads数的百分比;
IAV%:流感病毒reads数占高质量reads数的百分比;
Others:其他reads数占高质量reads数的百分比。
在对宿主鸡进行数据统计,将最终的contig序列与已有的病原微生物数据库进行比对,每一个样本分别选择8个目标片段,对这8个基因片段再进行分别的统计,统计结果如表2所示,表2是甲型流感病毒病毒基因组组装结果。
表2甲型流感病毒病毒基因组组装结果
其中,Gene:流感病毒基因节段;
Accession:参考序列GenBank号;
Length:组装序列长度;
Coverage:组装序列长度占参考序列长度的百分比;
Depth:组装序列中每个碱基的平均测序深度;
≥10X:组装序列中测序深度≥10X的碱基百分比;
根据表2的统计结果表明,采用本发明的方法,三个随机样本中,每个样本中的所有8个目标片段均被完整(100%)组装出来,并且目标片段的所有位点的覆盖度均在10X以上,准确度很高。
本发明的组装方法适用于各类病原微生物,如细菌、真菌、病毒等的基因组组装。
本发明采用的实施例只是较佳的实施例,用到的技术软件或物品如果没有明确说明,皆采用现有技术中常规使用标准的软件或物品,其他未述及的方法采用现有技术。
本领域的普通技术人员将会意识到,这里所示的实施例是为了帮助读者理解本发明的原理,应理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其他各方面变换仍然在本发明的保护范围内。
Claims (9)
1.一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,具体步骤如下:
(1)使用fastp对二代下机数据进行质控,得到质控后高质量的二代数据;
(2)使用bwa将质控后高质量的二代数据比对到宿主数据库,得到比对情况;
(3)使用samtools软件,根据步骤(2)中比对情况提取未比对上的数据,得到去宿主的reads,即Clean reads;
(4)使用soapdenovo对去宿主的reads进行组装,得到无参组装的contig序列;
(5)使用bwa将去宿主的reads比对到病原数据库,得到去宿主的reads比对情况;
(6)使用samtools软件统计去宿主的reads比对情况中每个位点的测序深度,得到有参组装的contig序列;
(7)将无参组装的contig序列和有参组装的contig序列进行整合,得到合并后的contig序列;
(8)使用去宿主的reads,以合并后的contig序列为参考进行比对,并重复步骤6,继续进行有参组装,得到最终版的contig序列;
(9)使用blast软件,将最终的contig序列与已有的病原微生物数据库进行比对,根据比对结果中序列之间的相似性与一致性的参数,对contig序列进行病原属种和亚型鉴定。
2.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(1)的具体方法为,根据实验信息与数据的实际情况,使用fastp软件对原始下机数据进行过滤,除去接头序列,低质量序列以及长度不符合要求的序列,得到符合质控要求的下机数据。
3.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(2)的具体方法为,质控后的下机数据中同时存在宿主的DNA序列以及病原微生物的DNA序列,使用bwa软件,以宿主物种的基因组作为参考基因组,将质控后的下机数据与参考基因组进行比对,获取下机数据中所有序列的比对情况,即all_bam文件。
4.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(3)的具体方法为,根据步骤(2)中得到的比对情况,使用samtools软件对all_bam文件中所有序列的比对情况进行过滤,对比对到参考基因组也就是宿主基因组的序列进行过滤,提取未比对上的所有序列,得到最终的去宿主的病原微生物的序列集合,即Clean reads。
5.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(4)的具体方法为,使用soapdenovo软件对Clean reads进行组装,得到不依赖病原微生物参考基因组的contig序列。
6.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(5)的具体方法为,使用bwa软件,以病原微生物的基因组作为参考基因组,将Clean reads比对到参考基因组上,获得Clean reads在病原微生物参考基因组上的比对情况,即pathogen_bam文件。
7.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(6)的具体方法为,使用samtools软件统计步骤(5)中pathogen_bam文件每个位点的测序深度,将每个位置上占比最高的碱基确定为该位置的碱基序列,并依次进行延伸;对于比对中出现的indel,根据位点两端序列以及位置信息,保留reads中的indel,得到有参组装的contig序列。
8.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(7)的具体方法为,使用自主编写脚本,以有参组装的contig序列为骨架,根据无参组装contig序列中的大片段易位或者倒位信息对骨架序列进行调整或重排,从而将无参组装的contig序列和有参组装的contig序列进行合并,得到合并后的contig序列。
9.根据权利要求1所述的一种基于宏基因组测序数据组装病原微生物基因组的方法,其特征在于,所述步骤(8)的具体方法为,使用bwa软件,以步骤(7)中得到的合并的contig序列作为参考基因组,将Clean reads比对到参考基因组上,得到最终的比对情况,即final_bam文件,并参照步骤(6)中的方法,对得到的final_bam文件进行二次有参组装,得到最终版的contig序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210732531.4A CN115101129B (zh) | 2022-06-27 | 2022-06-27 | 一种基于宏基因组测序数据组装病原微生物基因组的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210732531.4A CN115101129B (zh) | 2022-06-27 | 2022-06-27 | 一种基于宏基因组测序数据组装病原微生物基因组的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115101129A true CN115101129A (zh) | 2022-09-23 |
CN115101129B CN115101129B (zh) | 2023-03-24 |
Family
ID=83292374
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210732531.4A Active CN115101129B (zh) | 2022-06-27 | 2022-06-27 | 一种基于宏基因组测序数据组装病原微生物基因组的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115101129B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190252042A1 (en) * | 2016-10-13 | 2019-08-15 | bioMérieux | Indentification and antibiotic characterization of pathogens in metagenomic sample |
CN110349629A (zh) * | 2019-06-20 | 2019-10-18 | 广州赛哲生物科技股份有限公司 | 一种利用宏基因组或宏转录组检测微生物的分析方法 |
CN111009286A (zh) * | 2018-10-08 | 2020-04-14 | 深圳华大因源医药科技有限公司 | 对宿主样本进行微生物分析的方法和装置 |
CN113539378A (zh) * | 2021-07-16 | 2021-10-22 | 明科生物技术(杭州)有限公司 | 一种病毒数据库的数据分析方法、系统、设备及存储介质 |
CN113689912A (zh) * | 2020-12-14 | 2021-11-23 | 广东美格基因科技有限公司 | 基于宏基因组测序的微生物对比结果校正的方法和系统 |
CN114121160A (zh) * | 2021-11-25 | 2022-03-01 | 广东美格基因科技有限公司 | 一种检测样本中宏病毒组的方法和系统 |
-
2022
- 2022-06-27 CN CN202210732531.4A patent/CN115101129B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190252042A1 (en) * | 2016-10-13 | 2019-08-15 | bioMérieux | Indentification and antibiotic characterization of pathogens in metagenomic sample |
CN111009286A (zh) * | 2018-10-08 | 2020-04-14 | 深圳华大因源医药科技有限公司 | 对宿主样本进行微生物分析的方法和装置 |
CN110349629A (zh) * | 2019-06-20 | 2019-10-18 | 广州赛哲生物科技股份有限公司 | 一种利用宏基因组或宏转录组检测微生物的分析方法 |
CN113689912A (zh) * | 2020-12-14 | 2021-11-23 | 广东美格基因科技有限公司 | 基于宏基因组测序的微生物对比结果校正的方法和系统 |
CN113539378A (zh) * | 2021-07-16 | 2021-10-22 | 明科生物技术(杭州)有限公司 | 一种病毒数据库的数据分析方法、系统、设备及存储介质 |
CN114121160A (zh) * | 2021-11-25 | 2022-03-01 | 广东美格基因科技有限公司 | 一种检测样本中宏病毒组的方法和系统 |
Non-Patent Citations (4)
Title |
---|
商商SUTKWKB7DV: "宏基因组专题 | 宏基因组分析流程,基于reads VS 基于组装", 《个人图书馆》 * |
孙新蕾 等: "宏基因组测序在中枢神经系统感染性疾病诊断中的应用及研究进展", 《微生物学报》 * |
梅楚刚等: "基于高通量测序的动物基因组研究进展", 《西北农林科技大学学报(自然科学版)》 * |
郑智俊等: "人体肠道宏基因组生物信息分析方法", 《微生物学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115101129B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110349629B (zh) | 一种利用宏基因组或宏转录组检测微生物的分析方法 | |
CN109273053B (zh) | 一种高通量测序的微生物数据处理方法 | |
CN107810277B (zh) | 用于高通量微生物学应用的高分辨率系统、试剂盒、设备和方法 | |
TW201713775A (zh) | 用於自複雜的異源群集中分析微生物品系、預測及識別其機能性關係及交互作用,且依據前述選擇並合成微生物系集之方法、裝置及系統 | |
CN105420375B (zh) | 一种环境微生物基因组草图的构建方法 | |
CN109923217A (zh) | 宏基因组样品中病原体的鉴定和抗生素表征 | |
WO2022028624A1 (zh) | 通过测序获取微生物物种及相关信息的方法、装置、计算机可读存储介质和电子设备 | |
CN111187813B (zh) | 全流程质控的病原微生物高通量测序检测方法 | |
CN109610008A (zh) | 基于高通量测序的中枢系统感染病原检测文库构建方法、检测方法和试剂盒 | |
AU2015263055B2 (en) | Sequencing process | |
WO2017129110A1 (zh) | 一种人体微生物定性与定量的检测方法 | |
US20140288844A1 (en) | Characterization of biological material in a sample or isolate using unassembled sequence information, probabilistic methods and trait-specific database catalogs | |
CN109913461B (zh) | 根尖牙乳头干细胞外泌体piRNA生物标志物及其筛选应用 | |
CN114974411A (zh) | 宏基因组病原微生物基因组数据库及其构建方法 | |
CN110875082A (zh) | 一种基于靶向扩增测序的微生物检测方法和装置 | |
CN115101129B (zh) | 一种基于宏基因组测序数据组装病原微生物基因组的方法 | |
CN112342284B (zh) | 一种分析微生物群落功能基因转录和翻译活性的方法 | |
Mayali et al. | Chip-SIP: stable isotope probing analyzed with rRNA-targeted microarrays and NanoSIMS | |
EP2978858B1 (de) | In vitro-verfahren zur prognostischen beurteilung der erfolgsaussichten einer implantation und/oder transplantation | |
CN113470752B (zh) | 一种基于纳米孔测序仪的细菌测序数据鉴定方法 | |
CN115691679A (zh) | 一种基于二代和三代测序技术的宏病毒组分析方法 | |
CN116504319A (zh) | 一种宏基因组磷元素循环注释绘图方法 | |
CN113293201A (zh) | 一种新型氯霉素抗性基因的快速定量检测方法及其应用 | |
CN110684830A (zh) | 一种石蜡切片组织rna分析方法 | |
CN117690485B (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 |