CN110804645B - 一种评价香蕉镰刀菌枯萎病抗病性的方法 - Google Patents

一种评价香蕉镰刀菌枯萎病抗病性的方法 Download PDF

Info

Publication number
CN110804645B
CN110804645B CN201911079667.4A CN201911079667A CN110804645B CN 110804645 B CN110804645 B CN 110804645B CN 201911079667 A CN201911079667 A CN 201911079667A CN 110804645 B CN110804645 B CN 110804645B
Authority
CN
China
Prior art keywords
disease
grade
equal
banana
resistance
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
CN201911079667.4A
Other languages
English (en)
Other versions
CN110804645A (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.)
Pomology Research Institute Guangdong Academy of Agricultural Sciences
Original Assignee
Pomology Research Institute Guangdong Academy of Agricultural Sciences
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 Pomology Research Institute Guangdong Academy of Agricultural Sciences filed Critical Pomology Research Institute Guangdong Academy of Agricultural Sciences
Priority to CN201911079667.4A priority Critical patent/CN110804645B/zh
Publication of CN110804645A publication Critical patent/CN110804645A/zh
Priority to ZA2020/06929A priority patent/ZA202006929B/en
Application granted granted Critical
Publication of CN110804645B publication Critical patent/CN110804645B/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/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/18—Testing for antimicrobial activity of a material
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Microbiology (AREA)
  • Toxicology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Analytical Chemistry (AREA)
  • Biophysics (AREA)
  • Software Systems (AREA)
  • Biotechnology (AREA)
  • Immunology (AREA)
  • Evolutionary Biology (AREA)
  • Molecular Biology (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Breeding Of Plants And Reproduction By Means Of Culturing (AREA)

Abstract

本发明公开了一种对香蕉镰刀菌枯萎病的抗病性进行评价的方法，属于农业—植物保护技术领域。与田间鉴定法和苗期人工接种鉴定法不同，生根试管苗离体接种鉴定法是在无菌条件下将Foc接种到香蕉试管苗的基部，然后按照1‑6级的病害评价等级进行抗病性鉴定。鉴于病害级数在病害严重度上是非线性的，利用Logistic回归分析进行模型构建和发病等级概率预测，再根据预测结果进一步将香蕉品种的抗病性划分为高抗、抗病、中抗、感病和高感等5个级别。本发明提供的方法实现了生根试管苗离体接种鉴定法和其它鉴定法的整合，通过对采用不同鉴定方法获得的结果进行比对，并在此基础上调整实验方案和技术路线，推动抗病育种工作的进程。

Description

一种评价香蕉镰刀菌枯萎病抗病性的方法
技术领域
本发明属于农业-植物保护技术领域，特别涉及一种评价香蕉镰刀菌枯萎病抗病性的方法。
背景技术
香蕉镰刀菌枯萎病抗病性的鉴定是抗病育种工作的重要环节，主要分为田间鉴定法、苗期人工接种鉴定法和生根试管苗离体接种鉴定法等。
虽然田间的病害观察结果仍然是对抗性水平进行评价的最终依据，但在田间开展香蕉镰刀菌枯萎病的抗性鉴定，需要有大面积发病均匀的地块，并且成本较高。抗枯萎病香蕉新品种引种评价(刘文清等)以及几个香蕉品种对枯萎病的抗性与主要性状比较(黄素梅等)均报道了，研究学者在田间发病率的基础上，将香蕉种质的抗病性划分为高抗、抗病、中抗、感病和高感五个级别。
苗期人工接种鉴定法：在温室大棚中建立盆栽系统或水培系统，接着将病原菌(Fusarium oxysporum f.sp.cubense,Foc)接种到香蕉苗的根部并记录发病情况。根据调查结果，按照下列公式计算各个品种的病情指数：
Figure BDA0002263545860000011
最后依据苗期病情指数将供试香蕉品种的抗病性划分为高抗、抗病、中抗、感病和高感等5个级别。相关文献报道如：18份广东香蕉种质对枯萎病的抗性评价(宋晓兵等)，香蕉对尖孢镰刀菌热带4号小种的抗性评价方法的建立(左存武等)。
生根试管苗离体接种鉴定法是在无菌条件下将Foc接种到香蕉生根试管苗的基部，接种后将培养容器放置在组培室中观察发病情况，然后按照1-6级的病害评价等级对单株生根试管苗进行病害等级鉴定。由于生根试管苗离体接种鉴定法的病害级数在病害严重度上是非线性的，代入前述病情指数计算公式所得的数值并不能反映真正的发病程度。也即在病害分级没有数量界限的情况下，病情指数计算公式并不适用。
综上所述，根据病害等级数据的特点构建合理的数学模型，并在此基础上进一步科学地划分抗病性级别成为亟需解决的问题。
发明内容
本发明的目的在于提供一种评价香蕉镰刀菌枯萎病抗病性的方法。
为达到上述目的，本发明所采取的技术方案是：
一种评价香蕉镰刀菌枯萎病抗病性的方法，包括以下步骤：
1)数据收集：观察获得香蕉生根试管苗接种Foc后的病害等级数；
2)基于累积Logistic回归的模型构建：
Logistic回归模型对应包含5个logit函数：
Figure BDA0002263545860000021
Figure BDA0002263545860000022
Figure BDA0002263545860000023
Figure BDA0002263545860000024
Figure BDA0002263545860000025
其中：p1，p2，p3，p4，p5，p6为事件概率，分别代表1-6级病害评价等级，用于对比的基础水平是第6级；xk(k＝1，2，...，K)代表第k个香蕉品种；β0j(j＝1,2,…,5)代表回归的截距项；βk(k＝1，2，...，K)代表回归系数。每个logit函数具有相同的系数项和不同的截距项，各累积logit的回归线是相互平行的；
Logistic回归模型使用的估计方法为极大似然法，根据前述为预测香蕉生根试管苗发病程度设计的Logistic模型函数，得到的累积Logistic回归模型描述如下：
Figure BDA0002263545860000026
其中y′代表香蕉生根试管苗的发病情况；α代表截距项；βk(k＝1，2，...，K)代表回归系数；xk(k＝1，2，...，K)代表第k个香蕉品种；ε是误差项；
3)计算累积概率：
将生根试管苗离体接种鉴定法设定的6个病害等级赋予相应的取值y＝1，y＝2，...，y＝6，
各y值间关系为(y＝1)＜(y＝2)＜...＜(y＝6)，共有5个分界线分开相邻类别：
如果y′≤μ1，则y＝1；
如果μ1＜y′≤μ2，则y＝2；
如果μ2＜y′≤μ3，则y＝3；
如果μ3＜y′≤μ4，则y＝4；
如果μ4＜y′≤μ5，则y＝5；
如果μ5＜y′，则y＝6；
μj是分开类别的分界点，且μ1＜μ2＜μ3＜μ4＜μ5；
则累积概率的计算公式如下：
Figure BDA0002263545860000031
由此可以得到某香蕉品种的生根试管苗处于某病害等级的概率，计算方法如下：
P(y＝1)＝P(y≤1)
P(y＝2)＝P(y≤2)-P(y≤1)
P(y＝3)＝P(y≤3)-P(y≤2)
P(y＝4)＝P(y≤4)-P(y≤3)
P(y＝5)＝P(y≤5)-P(y≤4)
P(y＝6)＝1-P(y≤5)
且各等级概率值和为1，即P(y＝1)+P(y＝2)+…+P(y＝6)＝1；
4)根据病害等级的预测概率进一步划分抗病性级别。
根据本发明的实施例，步骤1)所述病害等级分为1-6级，1级-假茎基部的小叶片枯萎，但假茎本身的颜色并未发生变化；2级-假茎颜色变暗的区域小于或等于整个假茎高度的1/2；3级-假茎颜色变暗的区域超过整个假茎高度的1/2；4级-变黄或枯萎的上部叶片数小于或等于小植株上部总叶片数的1/2；5级-变黄或枯萎的上部叶片数超过小植株上部总叶片数的1/2；6级-整株试管苗枯萎死亡。
根据本发明的实施例，步骤4)抗病性级别为：当P(y＝2)≥50％，P(y＝1)≤50％，为高抗。
根据本发明的实施例，步骤4)抗病性级别为：当P(y＝2)≥50％，P(y＝3)≤50％，为抗病。
根据本发明的实施例，步骤4)抗病性级别为：当P(y＝3)≥50％，P(y＝2)≤50％或P(y＝4)≤50％，为中抗。
根据本发明的实施例，步骤4)抗病性级别为：当P(y＝4)≥50％，P(y＝3)≤50％或P(y＝5)≤50％，为感病。
根据本发明的实施例，步骤4)抗病性级别为：当P(y＝5)≥50％，P(y＝4)≤50％或P(y＝6)≤50％，为高感。
本发明的有益效果是：
Logistic回归模型属于概率型非线性回归模型，它对数据的正态性、方差齐性以及自变量类型等不做要求，并且具有系数的可解释性等优点，采用Logistic回归分析病害等级数据将大大提高生根试管苗离体接种鉴定法的准确性；根据发病等级概率的预测结果，进一步将供试香蕉品种的抗病性划分为高抗、抗病、中抗、感病和高感等5个级别，实现了生根试管苗离体接种鉴定法和其它鉴定法的整合，通过对采用不同鉴定方法获得的结果进行比对，并在此基础上调整实验方案和技术路线，可以有效地推动香蕉抗枯萎病育种的进程。
附图说明
图1为利用本发明的方法对香蕉生根试管苗的病害等级数据进行分析的结果。
具体实施方式
下面结合实施例对本发明中的技术方案进行清楚、完整的说明，但并不局限于此。
实施例1
实验材料：巴西香蕉、抗枯5号香蕉、宝岛蕉、农科1号香蕉、广粉1号粉蕉和东莞大蕉等6个品种。
方法：按照一种对香蕉镰刀菌枯萎病的抗病性进行快速鉴定的方法(专利号为：ZL200910192176.0)，建立香蕉—Foc离体互作体系。接种Foc后24天，观察记录生根试管苗的发病情况，按照1-6级的病害评价等级(表1)，对单株生根试管苗进行病害等级鉴定。
表1香蕉镰刀菌枯萎病生根试管苗离体接种鉴定法的病害等级
病害等级 分级标准
1级 假茎基部的小叶片枯萎，但假茎本身的颜色并未发生变化
2级 假茎颜色变暗的区域小于或等于整个假茎高度的1/2；
3级 假茎颜色变暗的区域超过整个假茎高度的1/2；
4级 变黄或枯萎的上部叶片数小于或等于小植株上部总叶片数的1/2；
5级 变黄或枯萎的上部叶片数超过小植株上部总叶片数的1/2；
6级 整株试管苗枯萎死亡。
将各品种的病害等级录入SPSS 19.0，执行Logistic回归分析，结果见图1。
从图1可以看出：巴西香蕉5级病害的发生概率为70％、6级病害的发生概率为30％；广粉1号粉蕉4级病害的发生概率为70％、5级病害的发生概率为30％；宝岛蕉3级病害的发生概率为50％、4级病害的发生概率为50％；农科1号香蕉2级病害的发生概率为80％、3级病害的发生概率为20％；抗枯5号香蕉1级病害的发生概率为30％、2级病害的发生概率为70％；东莞大蕉1级病害的发生概率为40％、2级病害的发生概率为60％。
根据图1病害等级概率的预测结果，对照本发明制定的抗病性级别划分标准(表2)，6个供试香蕉品种的抗病性级别划分为：高感(巴西香蕉)、感病(广粉1号粉蕉)、中抗(宝岛蕉)、抗病(农科1号香蕉)和高抗(抗枯5号香蕉和东莞大蕉)。
上述抗病性级别划分结果与实际情况相符。
表2病害等级的预测概率和抗病级别的关系
病害等级的预测概率 抗病级别
2级≥50％，1级≤50％ 高抗
2级≥50％，3级≤50％ 抗病
3级≥50％，2级≤50％或4级≤50％ 中抗
4级≥50％，3级≤50％或5级≤50％ 感病
5级≥50％，4级≤50％或6级≤50％ 高感
对比例：
苗期人工接种鉴定法：在温室大棚中建立盆栽系统或水培系统，接着将病原菌(Foc)接种到香蕉苗的根部并记录发病情况。按照下列公式计算各个品种的病情指数：
Figure BDA0002263545860000051
依据苗期病情指数，将供试香蕉品种的抗病性划分为高抗、抗病、中抗、感病和高感等5个级别。相关文献报道可见18份广东香蕉种质对枯萎病的抗性评价(宋晓兵等)，香蕉对尖孢镰刀菌热带4号小种的抗性评价方法的建立(左存武等)。
对比例的生根试管苗离体接种鉴定法的病害级数在病害严重度上是非线性的，代入前述病情指数计算公式所得的数值并不能反映真正的发病程度。也即在病害分级没有数量界限的情况下，病情指数计算公式并不适用。根据病情指数计算公式所得的数值划分抗病性级别不完全符合实际。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (7)

1.一种评价香蕉镰刀菌枯萎病抗病性的方法，其特征在于，包括以下步骤：
1)数据收集：观察获得香蕉生根试管苗接种Foc后的病害等级数；
2)基于累积Logistic回归的模型构建：
Figure FDA0002263545850000011
其中y′代表香蕉生根试管苗的发病情况；α代表截距项；βk(k＝1，2，...，K)代表回归系数；xk(k＝1，2，...，K)代表第k个香蕉品种；ε是误差项；
3)计算累积概率：
将生根试管苗离体接种鉴定法设定的6个病害等级赋予相应的取值y＝1，y＝2，...，y＝6，各y值间关系为(y＝1)＜(y＝2)＜...＜(y＝6)，共有5个分界线分开相邻类别：
如果y′≤μ1，则y＝1；
如果μ1＜y′≤μ2，则y＝2；
如果μ2＜y′≤μ3，则y＝3；
如果μ3＜y′≤μ4，则y＝4；
如果μ4＜y′≤μ5，则y＝5；
如果μ5＜y′，则y＝6；
μj是分开类别的分界点，且μ1＜μ2＜μ3＜μ4＜μ5；
则累积概率的计算公式如下：
Figure FDA0002263545850000012
由此得到某香蕉品种的生根试管苗处于某病害等级的概率，且各等级概率值和为1，即P(y＝1)+P(y＝2)+…+P(y＝6)＝1；
4)根据病害等级的预测概率进一步划分抗病性级别。
2.根据权利要求1所述的方法，其特征在于，步骤1)所述病害等级分为1-6级，1级-假茎基部的小叶片枯萎，但假茎本身的颜色并未发生变化；2级-假茎颜色变暗的区域小于或等于整个假茎高度的1/2；3级-假茎颜色变暗的区域超过整个假茎高度的1/2；4级-变黄或枯萎的上部叶片数小于或等于小植株上部总叶片数的1/2；5级-变黄或枯萎的上部叶片数超过小植株上部总叶片数的1/2；6级-整株试管苗枯萎死亡。
3.根据权利要求1所述的方法，其特征在于，步骤4)抗病性级别为：当P(y＝2)≥50％，P(y＝1)≤50％，为高抗。
4.根据权利要求1所述的方法，其特征在于，步骤4)抗病性级别为：当P(y＝2)≥50％，P(y＝3)≤50％，为抗病。
5.根据权利要求1所述的方法，其特征在于，步骤4)抗病性级别为：当P(y＝3)≥50％，P(y＝2)≤50％或P(y＝4)≤50％，为中抗。
6.根据权利要求1所述的方法，其特征在于，步骤4)抗病性级别为：当P(y＝4)≥50％，P(y＝3)≤50％或P(y＝5)≤50％，为感病。
7.根据权利要求1所述的方法，其特征在于，步骤4)抗病性级别为：当P(y＝5)≥50％，P(y＝4)≤50％或P(y＝6)≤50％，为高感。
CN201911079667.4A 2019-11-07 2019-11-07 一种评价香蕉镰刀菌枯萎病抗病性的方法 Active CN110804645B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911079667.4A CN110804645B (zh) 2019-11-07 2019-11-07 一种评价香蕉镰刀菌枯萎病抗病性的方法
ZA2020/06929A ZA202006929B (en) 2019-11-07 2020-11-06 Method for evaluating fusarium wilt resistance in banana

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911079667.4A CN110804645B (zh) 2019-11-07 2019-11-07 一种评价香蕉镰刀菌枯萎病抗病性的方法

Publications (2)

Publication Number Publication Date
CN110804645A CN110804645A (zh) 2020-02-18
CN110804645B true CN110804645B (zh) 2022-06-28

Family

ID=69501442

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911079667.4A Active CN110804645B (zh) 2019-11-07 2019-11-07 一种评价香蕉镰刀菌枯萎病抗病性的方法

Country Status (2)

Country Link
CN (1) CN110804645B (zh)
ZA (1) ZA202006929B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101643772A (zh) * 2009-09-09 2010-02-10 广东省农业科学院果树研究所 一种对香蕉镰刀菌枯萎病的抗病性进行快速鉴定的方法
WO2012069768A1 (fr) * 2010-11-24 2012-05-31 Assistance Publique - Hôpitaux De Marseille Procede d'identification d'un champignon, notamment d'un champignon filamenteux, et utilisations

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101643772A (zh) * 2009-09-09 2010-02-10 广东省农业科学院果树研究所 一种对香蕉镰刀菌枯萎病的抗病性进行快速鉴定的方法
WO2012069768A1 (fr) * 2010-11-24 2012-05-31 Assistance Publique - Hôpitaux De Marseille Procede d'identification d'un champignon, notamment d'un champignon filamenteux, et utilisations

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Logistic 模型在植病流行学中应用的有关问题的探讨;檀根甲等;《安徽农业大学学报(自然科学版)》;19980415(第02期);138-141 *
Serial passage through resistant and susceptible cucumber cultivars affects the virulence of Fusarium oxysporum f. sp. cucumerinum;Xiaoqing Huang 等;《Microbiology Open》;20181231(第8期);e641 *
小麦品种对赤霉病的慢病性研究;胡小平等;《西北农业学报》;20030322(第01期);121-123 *
植物病害I-S关系――logistic衍生模型的研究;张连翔等;《沈阳农业大学学报》;20011025(第04期);270-273 *
香蕉枯萎病抗性离体接种鉴定方法的优化;吴元立等;《园艺学报》;20200825(第08期);154-161 *

Also Published As

Publication number Publication date
CN110804645A (zh) 2020-02-18
ZA202006929B (en) 2021-09-29

Similar Documents

Publication Publication Date Title
CN108346142B (zh) 一种基于植物光照图像的植物生长状态识别方法
Ozrenk et al. Levels of fire blight (Erwinia amylovora) susceptibility of native apple, pear and quince germplasm from Lake Van Basin, Turkey
Khadivi-Khub et al. Analysis of grape germplasm from Iran based on fruit characteristics
Faruq et al. Morphological charecterization of kenaf (Hibiscus cannabinus L.) in Malaysian tropical environment using multivariate analysis
CN109146119B (zh) 小麦耐荫性综合鉴定模型、方法及应用
Alizadeh-Moghaddam et al. Genomic markers analysis associated with resistance to Alternaria alternata (fr.) keissler—tomato pathotype, Solanum lycopersicum L.
CN108486214B (zh) 一种荞麦茎枯病抗性鉴定方法
CN110804645B (zh) 一种评价香蕉镰刀菌枯萎病抗病性的方法
CN106171898A (zh) 一种苗期鉴定花生黑腐病抗性品种的方法
CN112111551A (zh) 一种马铃薯黑胫病抗性的鉴定方法
Khan et al. Influence of inoculum density of Verticillium dahliae on root discoloration of horseradish
Viana et al. Genetic diversity in yellow passion fruit populations
CN112710661B (zh) 一种马铃薯育种监测分析方法及系统
Jafari et al. Evaluation of herbage yield, agronomic traits and powdery mildew disease in 35 populations of sainfoin (Onobrychis sativa) across 5 environments of Iran
Kabi et al. Phenotyping of sesame (Sesamum indicum L.) genotypes for powdery mildew resistance
Ibáñez et al. Characterization of the reproductive performance of a collection of grapevine cultivars
Prashantha et al. Identification of resistant sources against bakanae disease in short grained aromatic rice (Oryza sativa)
CN111524022A (zh) 一种植物品种dus测试方法
Russell The whole embryo method of testing barley for loose smut as a routine test
CN111367967A (zh) 一种植物品种一致性测试方法
da Costa Filho et al. Reaction of watermelon accessions to Meloidogyne enterolobii.
CN112163769B (zh) 一种利用蛛网图的茶树品种综合评价方法
CN114051853B (zh) 一种离体春梢快速鉴定高感松材线虫病马尾松的方法及其应用
De et al. Numerical Taxonomic analysis for the estimation of Genetic Diversity among some traditional rice (Oryza sativa L.) varieties of West Bengal
Singh et al. Evaluation of genetic divergence in Asiatic cotton

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