CN1717585A - 利用微分方程从时序基因表达数据推断基因调控网络 - Google Patents
利用微分方程从时序基因表达数据推断基因调控网络 Download PDFInfo
- Publication number
- CN1717585A CN1717585A CNA2003801040561A CN200380104056A CN1717585A CN 1717585 A CN1717585 A CN 1717585A CN A2003801040561 A CNA2003801040561 A CN A2003801040561A CN 200380104056 A CN200380104056 A CN 200380104056A CN 1717585 A CN1717585 A CN 1717585A
- Authority
- CN
- China
- Prior art keywords
- overbar
- gene
- sigma
- centerdot
- following formula
- 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.)
- Pending
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 101
- 230000014509 gene expression Effects 0.000 title claims abstract description 64
- 230000001105 regulatory effect Effects 0.000 title description 21
- 238000000034 method Methods 0.000 claims abstract description 61
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 238000005259 measurement Methods 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000002299 complementary DNA Substances 0.000 claims description 4
- 230000035484 reaction time Effects 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 abstract description 2
- 238000007476 Maximum Likelihood Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 13
- 244000063299 Bacillus subtilis Species 0.000 description 12
- 235000014469 Bacillus subtilis Nutrition 0.000 description 12
- 102000004169 proteins and genes Human genes 0.000 description 8
- 238000011161 development Methods 0.000 description 7
- 239000003814 drug Substances 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 230000003993 interaction Effects 0.000 description 6
- 108091008053 gene clusters Proteins 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 108020004999 messenger RNA Proteins 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 101150029062 15 gene Proteins 0.000 description 3
- 101150096316 5 gene Proteins 0.000 description 3
- 108020004635 Complementary DNA Proteins 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000010804 cDNA synthesis Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 101710134784 Agnoprotein Proteins 0.000 description 2
- 230000010558 Gene Alterations Effects 0.000 description 2
- 101710124584 Probable DNA-binding protein Proteins 0.000 description 2
- 240000004808 Saccharomyces cerevisiae Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008827 biological function Effects 0.000 description 2
- 230000001364 causal effect Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000002609 medium Substances 0.000 description 2
- 238000002493 microarray Methods 0.000 description 2
- 238000012775 microarray technology Methods 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- 238000000528 statistical test Methods 0.000 description 2
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 1
- 101150000874 11 gene Proteins 0.000 description 1
- 101150066838 12 gene Proteins 0.000 description 1
- 101150092328 22 gene Proteins 0.000 description 1
- 101150029857 23 gene Proteins 0.000 description 1
- 101150110188 30 gene Proteins 0.000 description 1
- 101150072006 33 gene Proteins 0.000 description 1
- 101150086149 39 gene Proteins 0.000 description 1
- 101150017816 40 gene Proteins 0.000 description 1
- 101150106774 9 gene Proteins 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- ZDXPYRJPNDTMRX-VKHMYHEASA-N L-glutamine Chemical compound OC(=O)[C@@H](N)CCC(N)=O ZDXPYRJPNDTMRX-VKHMYHEASA-N 0.000 description 1
- 235000014680 Saccharomyces cerevisiae Nutrition 0.000 description 1
- 238000000692 Student's t-test Methods 0.000 description 1
- 241000192593 Synechocystis sp. PCC 6803 Species 0.000 description 1
- 230000037354 amino acid metabolism Effects 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 1
- 230000003570 biosynthesizing effect Effects 0.000 description 1
- 230000023852 carbohydrate metabolic process Effects 0.000 description 1
- 235000021256 carbohydrate metabolism Nutrition 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000022131 cell cycle Effects 0.000 description 1
- 210000000170 cell membrane Anatomy 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000009089 cytolysis Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- ZDXPYRJPNDTMRX-UHFFFAOYSA-N glutamine Natural products OC(=O)C(N)CCC(N)=O ZDXPYRJPNDTMRX-UHFFFAOYSA-N 0.000 description 1
- 239000001963 growth medium Substances 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 150000002611 lead compounds Chemical class 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000010208 microarray analysis Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000004963 pathophysiological condition Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 229940081969 saccharomyces cerevisiae Drugs 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 108700026220 vif Genes Proteins 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- 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
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- 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
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- 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
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
- G16B5/20—Probabilistic models
-
- 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
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- 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
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
- G16B5/10—Boolean models
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Probability & Statistics with Applications (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Evolutionary Computation (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Bioethics (AREA)
- Artificial Intelligence (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
通过利用时间过程表达数据和一组线性微分方程,本方法的实施例能够被用来估计生物体基因之间的网络关系。通过确定哪些元素为0或者哪些元素在所研究的条件下不显著地改变,能够使用Aikaike信息准则和掩模工具来减少矩阵中元素的数目。使用最大似然估计和新的统计方法来估计所提出的网络关系的显著性。
Description
相关申请:
本申请要求在35 U.S.C§119(e)于2002年11月25日提交的美国临时专利申请No.60/428,827的优先权。其全部内容结合于此作为参考。
技术领域
本发明涉及确定生物体基因之间关系的方法。特别地,本发明包括利用微分方程线性***,从时间过程基因表达数据(time course gene expressiondata)推断基因调控网络的新方法。
背景技术
在生命科学、医学、药物发明和发展及制药行业中,当前研究和发展的最重要的方面中的一个为,对发展用于解释大量的第一手数据并基于这些数据引出结论的方法和器件的需求。生物信息学对***生物学的理解做出了重大贡献并保证产生对生命***成分之间的复杂关系更重大的理解。特别地,随着用于迅速探测被表达的基因和定量基因表达的新方法的出现,生物信息学甚至在不确定的情况下能够用来预测潜在的治疗目标,特定的基因在生物体生态学中可能起的确切作用。
遗传***的模拟是***生物学的中心主题。因为模拟能够基于生物学知识,通过预测或者推断以前未知的关系,网络评估方法能够支持生物模拟。
特别地,微阵列技术允许来自多种生物体的很多基因表达的研究。大量的第一手数据能够从来自生物体的许多基因得到,并且能够通过干扰或者通过突变、疾病或药物来研究基因表达。在特定的疾病中或者响应特定的干扰而特定的基因表达增加的发现,可以使人相信基因直接牵涉到疾病过程或者药物响应。然而,在生物有机体中,基因很少通过任何这样的干扰独立地调控,因为许多基因能够被一特定的干扰影响。因为很多不同的基因可以被如此影响,在这样的研究中了解基因之间的因果关系是非常困难的。因而正花费大量的努力来发展用于确定基因之间因果关系的方法,哪种基因对一种生物现象重要,哪种基因的表达对正在研究的生物过程不重要。尽管这样的不重要的基因表达作为生物或病理生理条件下的标记可能有用,如果这样的基因对于生理学或者病理生理学条件不重要,基于这样的基因发展药物可能不值得努力。相反,对于对一个过程识别为重要的基因,药物或者其它干扰的发展可能对发展用于与改变的基因表达关联的条件的治疗至关重要。
微阵列技术允许在同一时间测量很多基因的基因表达水平。使用互补DNA(cDNA),微阵列分析能够容易地实现,但是RNA微阵列也能够用来研究基因表达。随着可利用的基因表达数据的数量迅速地增加,用来分析这种数据的技术仍在发展中。数学方法被日益地用来确定表达基因之间的关系。然而,从基因表达数据精确地导出基因调控网络会很困难。
在时序基因表达测量中,能够通过在少量的时间点测量基因表达水平来研究基因表达的时效模式。例如,已经在酿酒酵母菌(Saccharomycescerevisiae)的细胞周期测量了周期变化的基因表达水平(见参考文献1)。基因对缓慢改变的环境的响应已经在同种酵母菌的二次转变期间进行了测量(见参考文献2)。其它的实验测量响应生物体环境突然变化的时效基因表达模式。作为例子,在外部光强突然移位之后,测量了藻青菌集胞藻6803(Synechocystis sp.PCC 6803)的基因表达响应。
已经提出了几种方法来从表达数据推断基因相互关系(见参考文献2、5和6)。在簇分析中,基因基于他们的基因表达谱的相似点集群在一起。从测量的基因表达数据来推断布尔(Boolean)或贝氏(Bayesian)网络(见参考文献7、8、9、10、11和美国专利申请No:10/259.723及申请题目“NonlinearModeling of Gene Networks From Time Series Gene Expression Data”,2003年11月18日提交;律师案号No:GENN 1008 US1DBB,两个申请都全部结合于此作为参考)和使用任意的微分方程***模拟基因表达数据(见参考文献12)先前已经公开了。然而,为了可靠地推断这样一个任意的微分方程***,需要长序列的时序基因表达数据,这在当前常常还不可获得。
发明内容
为了克服现有技术的缺点,在本发明的某些方面,我们发展了使用微分方程线性***及从基因表达数据导出的信息来推断基因网络的方法。这种途径保持了微分方程的定量及内在因果关系的优点,而又足够简单易于计算处理。我们也发展了用于检验涉及基因调控网络的假定的新方法。
附图说明
参照其具体的例子,描述本发明的各个方面。本发明的其它特征能够通过参照附图理解,其中:
图1描绘来自枯草芽孢杆菌(Bacillus subtilis)的5条基因簇的基因表达随时间的图表。
图2描绘利用本发明的方法导出的在图1中描绘的5条基因簇的基因网络。
具体实施方式
Chen从理论上考虑了使用线性微分方程模拟生物数据(见参考文献13)。在这个模型中,mRNA和蛋白质浓度都由线性微分方程***描述。这样的***能够描述为:
在cDNA微阵列实验中,当蛋白质浓度未知时,通过测量相应的mRNA浓度通常仅只确定基因表达水平。因此我们集中在仅描述基因相互作用的微分方程***上。矩阵元Λij于是代表基因j对基因i的作用,[Λij]-1为反应时间。
为了从所测量的数据推断微分方程***中的系数,先前已提出离散微分方程***(见参考文献13),代替所测量的mRNA和蛋白质浓度,并求解因而产生的方程线性***以得到线性微分方程***中的系数Λij。该方程***通常是不确定的。使用附加的必要条件,即基因调控网络应该是稀疏的,Chen表示模型能够以o(mh+1)次构造,这里m是基因的数目以及h是***中每个微分方程所允许的非零系数的数目(见参考文献13)。
参数h特别选定,这会有两个意外的结果。因为矩阵
中每行将精确地具有h个非零元素,网络中每个基因或者蛋白质具有h个亲本基因或蛋白质,并且因此没有基因或者蛋白质能够存在于网络的顶点。其次,每个基因不可避免地为一个反馈环的成员。虽然反馈环可能存在于基因调控网络中,他们的存在应该从测量的数据中来确定而不是人为地创造。
另一方面,贝式网络,不允许环的存在。贝氏网络依赖于所估计的网络的联合概率分布,以便可分解为条件概率分布的积。仅在不存在环时这种分解是可能的。我们进一步注意到贝氏网络往往包含许多参数,并且因此需要大量的数据以进行可靠的估计。
因此,我们旨在找到允许网络中存在环的方法,但是不一定需要他们存在。利用方程1,通过限制可能出现在***中的非零系数的数目,我们构造了一个稀疏矩阵。代替特别选定这个数目,我们通过利用Akaike信息准则(Akaike’s Information Criterion)(AIC)从数据中估计在相互作用矩阵中哪些系数为0,允许基因调控路径的数目对于每个基因不相同。
我们的方法的多个方面能够用于寻找单独基因之间的网络,也可用于寻找基因簇之间的调控网络。作为例子,我们能够利用枯草芽孢杆菌的时间过程数据推断基因簇之间的基因调控网络。簇能够利用k平均聚类演算法(k-means clustering algorithm)创造。簇的生物功能能够从属于每个簇的基因的功能范畴来确定。
在一些实施例中,我们按照微分方程(方程1)线性***考虑m个基因之间的调控网络,这里矢量
x(t)包含在t时刻m个基因的表达率。这个微分方程***能够求解为:
其中x0包含在零时刻的基因表达率。在这个方程中,矩阵指数按照泰勒(Taylor)展开定义为:
或
通过利用这个方程,我们能够按照多维线性马尔可夫(Markov)模型有效地描绘基因调控网络。
能够假定误差有一个不依赖于时间的正态分布,如下所示:
在所有的时刻对于所有的基因有一个相等的标准偏差σ。对于一系列在时刻ti,i∈{1,...,n),n个时间点的时序测量值
x i的对数似然函数于是为:
其中
为在时刻ti从测量数据估计的测量误差。
方差σ2的最大似然估计能够通过关于σ2取对数似然函数的最大值得到。这得出:
将这个式子代入对数似然函数(方程8)得到:
(12)
其中矩阵
和
定义为:
在某些实施例中,如果总平方差因此增加很小,如方程12所给出的,能够设置该矩阵元等于零。正式地,我们将使用Akaike信息准则(见参考文献15、16):
AIC=2·[估计的模型的对数似然值]+2·[估计的参数的数目], (16)来决定哪个矩阵元素应该设置为等于零。通过比较模型中使用的参数的数目和所估计的模型中的总误差,AIC能够用来避免模型到数据的过拟合(overfitting)。具有最低AIC的模型被认为是最佳的。AIC基于信息论并且广泛地用于统计模型识别,尤其用于时序模型拟合(见参考文献17)。
这里о表示哈达玛(Hadamard)(基于元素的(element-wise))乘积,并且掩模
是一个元素不为1就为0矩阵。相应的总平方误差
能够通过在方程12中用
替换
得到。给出掩模
能够通过求解一组方程18最小化总平方误差,
如果
如果
估计的参数为
和我们所允许的非零的矩阵
元素。从这个方程,能够看到当平方误差减少时,随着非零元素的数目增加,AIC可能增加。现在可以通过找到对于AIC来说产生最低值的掩模
根据基因表达数据推断基因调控网络。
对于任何但是最平凡的情形,可能的掩模
的数目非常大,使得进行彻底地搜寻以找到最佳掩模不可行。替代地,我们能够使用贪婪(greedy)搜寻方法。开始,以对于每个掩模元素为1和0相等的概率,我们能够随机地选择掩模。通过改变每个掩模元素Mij,能够减小AIC。这个过程能够继续直到找到最后的掩模,对于这个掩模不能实现进一步AIC的减小。这种算法能够从不同的(例如,随机的)初始掩模重复开始,并且能够用来确定具有最小的相应AIC的最终掩模
如果在几十次试验中都找到这个最佳的掩模,能够合理地断定没有更好的掩模存在。
我们已经描述并证明了从测量基因表达数据以微分方程线性***的形式来推断基因调控网络的方法。由于进行典型测量的时间点的数目有限,得到基因调控网络通常是一个不确定的难题。因为生物学上作为结果的基因调控网络预期为稀疏的,我们设置某些矩阵项等于零,并且仅利用非零项推断网络。非零项的数目,并且因此,网络的稀疏性,利用Akaike信息准则从数据确定,而不利用任何特别参数。
按照微分方程描述基因网络至少具有三个优点。第一,微分方程组描述基因之间的因果关系:系数矩阵的系数Λij确定基因j对基因i的作用。第二,它以清楚的数字形式描述基因相互作用。第三,由于微分方程***中呈现的大量信息,其它网络形式能够容易地从它导出。另外,我们能够将推断的网络连接到其它分析或者可视化工具,如GON(Genomic ObjectNet)。
在先前描述的方法中,或者不能够找到任何环(如在贝氏网络模型中)或者该方法人为地在网络中产生回路。尽管这里描述的方法允许网络中出现环,但并不需要他们的存在。仅只在有数据保证时可找到环。例如,当利用MMGE培养基中枯草芽孢杆菌的时间过程数据推断基因簇之间的调控网络时,我们发现一些簇是环的部分,而其它的不是(参看下面的例子和图2)。
如果基因的数目m等于或者大于试验的数目n,方程18中的矩阵
是奇异的。这个问题于是就不确定,并且能够找到具有零总误差
和为-∞的AIC的相互作用矩阵
通过对足够小数量的基因或者基因簇应用这种方法,或者通过限制网络中亲本的数目,能够避免我们方法的这种不足。
用于评估网络关系的统计显著性的方法
在本发明的其它实施例中,提供了用于确定网络关系分析的统计显著性的方法。在无效假设(null hypothesis)下,能够假定基因不被实验操作影响。在不同时间点的测量对数比(log-ratio)于是相等。我们能够进一步假设对数比具有零平均的正态分布。在一些情形,统计测试,如学生的t-检验(Student’s t-test),会在每个时间点进行以确定哪个对数比显著地不同于零。然而,学生的t-检验会不可靠,因为仅只少许测量的数据组。因此,在包括每个时间点仅两次测量的数据组的一些实施例中,我们设计了一种新的统计测试,将多个时间点的测量合并在一起。特别地,如例2中所示,我们将这种方法应用到来自所有的8个时间点的数据上。令人满意的是该方法能够用于其它类型的实验,并且将在下面描述。
实现本方法的步骤在下面描述。
步骤1:在每个时间点,计算平均对数比为
在无效假设下,
xj·(在一个时间点两个基因表达对数比的平均)是具有零平均正态分布并且估计的标准偏差为
的随机变量。
步骤2:接着从所有的测量估计标准偏差(如,对于如例1中包括的数据组8×2=16):
其中xji[k]表示对于基因j在时间点i测量k的数据值。
步骤3:就绝对值而言比测量值
xji大的
xj·的联合概率为:
其中,erf是误差函数。对于这个乘积中的单个因素Pi,我们通常会挑选一个显著性水平α,并且如果Pi<α则舍弃无效假设。
步骤4:采用准则:P<αn用于舍弃无效假设。这允许我们通过利用关于那个基因所有可利用的数据来确定在实验期间该基因的表达水平是否显著地改变。
步骤5:确定一个基因改变的表达水平是否显著。
用于确定基因之间网络关系的方法和新的统计方法能够用在研究,生物医学科学,包括诊断学中,以便发展新的诊断及用于在制药工业中选择先导化合物。
例子
下面的例子旨在阐明本发明的实施例,而不限制其范围。能够发展其它的实施例而不背离本发明的范围,并且本发明的方法及其变体能够在不存在不适当的实验下用来推断枯草芽孢杆菌或其它生物体中不同基因的调控网络。所有这样的实施例被认为是本发明的部分。
例1:枯草芽孢杆菌中的基因网络
最近在枯草芽孢杆菌的MMGE基因表达实验中测量了用于利用基因表达数据寻找基因调控网络的本发明实施例。MMGE是包含葡萄糖和谷氨酰胺(如碳和氮源)的合成的基本培养基。在这种介质中,诱导了小分子生物合成所需要的基因的表达,如氨基酸。在这个实验中,在一个小时的时间间隔中于八个时间点测量了4320 ORF的表达水平,在每个时间点进行两次测量。
数据准备和分析
为减少出现在数据中的测量噪声的影响,每个基因的表达水平与测量背景水平比较。无论在红色或是绿色通道中具有低于平均背景水平的平均基因表达水平的基因从分析中除去。
对3823个剩下的基因应用全局归一化,并且计算了基因表达率以2为底的对数。我们对测量的对数比应用统计测试来确定它们是否显著地不同于零。
以上描述的方法的流程图重现在下面的总结中。
步骤1:计算在每个时间点,每个基因表达的平均对数比;
步骤2:从所有的测量计算标准偏差;
步骤3:计算联合概率;
步骤4:采用用于统计显著性的准则;及
步骤5:确定一个基因改变的表达水平是否显著。
在这个例子中,我们选择一个显著性水平α=0.00025以便预期的假阳性数(0.00025×3823=1)可接受。通过应用这个准则到这3823个基因,我们发现684个基因显著地受到影响。
例2:枯草芽孢杆菌基因的聚类
使用k平均聚类,枯草芽孢杆菌的这684个基因顺序地群集成5组。使用欧几里得(Euclidean)距离测量基因之间的距离,而簇的质心(centroid)定义为簇中所有基因的中部(median)。挑选簇的数目以避免显著的重叠。k平均算法从不同的随机初始聚类开始而重复1000000次。找到最佳解81次。
完整的聚类结果可在下列网站得到:
http://bonsai.ims.u-tokyo.ac.jp/-mdehoon/publications/Subtilis/
clusters.html.
为了确定所创造的簇的生物学功能,我们考虑每个簇中所有基因在枯草芽孢杆菌数据库中的功能范畴。表1列出所形成的5个簇的主要功能范畴。
图1表示对于每个簇,基因表达的对数比作为时间的函数。尽管在该时间过程期间,簇I、II和V的表达水平相当大地改变,簇II和III具有相当稳定的表达水平。特别地,簇IV能够看作为一个包罗万象的簇,被分配到其上的基因不能适合其它簇。
表1
利用k平均聚类创造的5个簇的主要功能范畴。
簇 | 基因的数目 | 主要功能范畴 |
I | 42 | 2.2:11基因;1.1:9基因 |
II | 62 | 1.2:15基因;2.2:12基因 |
III | 187 | 5.1:30基因;6.0:23基因;1.2:22基因 |
IV | 343 | 5.1:40基因;5.2:39基因;1.2:33基因 |
V | 50 | 1.2:15基因;2.1.1:15基因 |
基因的功能范畴
参考法国巴斯德研究所(Institute Pasteur)枯草芽孢杆菌数据库的功能范畴
1.1: 细胞壁1.2: 运输/结合蛋白质和脂蛋白2.1.1: 碳水化合物代谢和相关分子----特定通路2.2: 氨基酸代谢和相关分子5.1: 与来自枯草芽孢杆菌的未知蛋白质类似的5.2: 与来自其它生物体的未知蛋白质类似的6.0: 无相似性 |
图1表示对于每个簇,作为时间的函数的基因表达的对数比,其从测量的基因表达数据确定。
分段网络构造
从那12个基因的测量对数比中,我们构造矩阵
和
并计算矩阵
从一个随机的初始掩模开始,计算掩模
的过程重复1000次。找到最佳的解55次。因此,不太可能存在其它具有较低AIC的掩模。注意到可能的掩模的总数目为225=33,554,432。
图2表示所找到的网络。在网络中簇的亲本的数目在0和5之间变化。簇III和IV作为网络的顶点出现,而簇I、II和V连接到一个环中。注意到这种网络既不能由先前提出的方法(见参考文件13)产生,也不能由贝氏网络模型产生。
网络中两个最强的相互作用分别是簇IV对簇V和簇II的正效应和负效应。簇II和V的基因表达水平的相反行为极大可能地由簇IV引起,而不是簇II和V之间的直接相互作用。
参考文献
1. P.T.Spellman,G.Sherlock,M.Q.Zhang,V.R.Iyer,K.Anders,M.B.Eisen,P.O.Brown,D.Botstein,and B.Futcher,″Comprehensive identification of cellcycle-regulated genes of the yeast Saccharomyces cerevisiae by microarrayhybridization″Mol.Biol.Cell 9(1998)3273-3297.
2. J.L.DeRisi,V.R.Iyer,and P.O.Brown,″Exploring the metabolic and geneticcontrol of gene expression on a genomic scale″Science 278(1997)680-686.
3. Y.Hihara,A.Kamei,M.Kanehisa,A.Kaplan,and M.Ikeuchi,″DNAmicroarray analysis of cyanobacterial gene expression during acclimation to high light″The Plant Cell 13(2001)793-806.
4. M.J.L.de Hoon,S.Imoto,and S.Miyano,″Statistical analysis of a small set oftime-ordered gene expression data using linear splines″Bioinformatics,in press.
5. M.B.Eisen,P.T.Spellman,P.O.Brown,and D.Botstein,″Cluster analysis anddisplay of genome-wide expression patterns″Proc.Natl.Acad.Sci.USA 95(1998)14863-14868.
6. P.Tamayo,D.Slonim,J.Mesirov,Q.Zhu,S.Kitareewan,E.Dmitrovsky,E.S.Lander,and T.R.Golub,″Interpreting patterns of gene expression with self-organizingmaps:Methods and application to hematopoietic differentiation″Proc.Natl.Acad.Sci.USA 96(1999)2907-02912.
7. S.Liang,S.Fuhrman,and R.Somogyi,″REVEAL,a general reverseengineering algorithm for inference of genetic network architectures″Proc.Pac.Symp.on Biocomputing 3(1998)18-29.
8. T.Akutsu,S.Miyano,and S.Kuhara,″Inferring qualitative relations in geneticnetworks and metabolic pathways″Bioinformatics 16(2000)727-734.
9. N.Friedman,M.Linial,I.Nachman,and D.Pe′er,″Using Bayesian networks toanalyze expression data″J.Comp.Biol.7(2000)601-620.
10. S.Imoto,T.Goro,and S.Miyano,″Estimation of genetic networks andfunctional structures between genes by using Bayesian networks and nonparametricregression″Proc.Pac.Symp.on Biocomputing 7(2002)175-186.
11. S.Imoto,S.-Y.Kim,T.Goto,S.Aburatani,K.Tashiro,S.Kuhara,and S.Miyano,″Bayesian network and nonparametric heteroscedastic regression for nonlinearmodeling of genetic network″Proc.IEEE Computer Society Bioinformatics Conference(2002)219-227.
12. E.Sakamoto and H.Iba,″Evolutionary inference of a biological network asdifferential equations by genetic programming″Genome Informatics 12(2001)276-277.
13. T.Chen,H.L.He,and G.M.Church,″Modeling gene expression withdifferential equations″Proc.Pac.Symp.on Biocomputing 4(1999)29-40.
14. R.A.Hom and C.R.Johnson,Matrix Analysis.Cambridge University Press,Cambridge,UK(1999).
15. H.Akaike,″Information theory and an extension of the maximum likelihoodprinciple″Research Memorandum No.46,Institute of Statistical Mathematics,Tokyo(1971).In B.N.Petrov and F.Csaki(editors),2nd Int.Symp.on Inf.Theory.AkadémiaiKiadó,Budapest (1973)267-281.
16. H.Akaike,″A new look at the statistical model identification″IEEE Trans.Automat.Contr.AC-19(1974)716-723.
17. M.B.Priestley,Spectral Analysis and Time Series,Academic Press,London(1994).
18. Microbial Advanced Database Organization (Micado).http://www-mig.versailles.inra.fr/bdsi/Micado/.
19. I.Moszer,P.Glaser,and A.Danchin,″SubtiList:a relational database for theBacillus subtilis genome″Microbiology 141(1995)261-268.
20. I.Moszer,″The complete genome of Bacillus subtilis:From sequence annotationto data management and analysis″FEBS Letters 430(1998)28-36
21. T.W.Anderson and J.D.Finn,The New Statistical Analysis of Data.SpringerVerlag,New York(1996).
22. H.Matsuno,A.Doi,Y.Hirata,and S.Miyano,″XML documentation ofbiopathways and their simulation in Genomic Object Net″Genome Informatics 12(2001)54-62.Genomic Object Net is available at http://www.GenomicObject.net.
Claims (28)
1.一种用于推断基因之间网络关系的方法,包括:
(a)提供一个用于生物体的一组基因的定量时间过程数据库,所述库包括基于在所述基因组中的每个基因表达的时间过程的表达结果,量化所述基因彼此间的平均作用和每个时间点的变异性的测量;
(b)从所述库创造一个稀疏矩阵,所述矩阵具有从其上除去的零系数;
(c)从所述矩阵产生一组线性微分方程;及
(d)求解所述方程组以产生所述网络关系。
2.如权利要求1所述的方法,其中利用Akaike信息准则(AIC)来识别所述零系数。
4.如权利要求1-3中任何一个所述的方法,其中所述矩阵包含元素Λij,Λij代表基因j对基因i的作用,并且其中[Λij]-1代表所述基因j对基因i的作用的反应时间。
5.如权利要求1-4中任何一个所述的方法,其中所述求解的微分方程为:
6.如权利要求1-5中任何一个所述的方法,其中所述指数Λt(exp(Λ))根据以下表达求解:
7.如权利要求1-6中任何一个所述的方法,其中所述微分方程通过求解以下差分方程估计:
8.如权利要求1-7中任何一个所述的方法,其中,根据以下公式,所述稀疏矩阵进一步包括一个估计的误差:
9.如权利要求1-8中任何一个所述的方法,其中,根据以下公式,所述误差具有不依赖于时间的正态分布:
其中标准偏差σ在所有时刻对于所述基因的每一个相等。
10.如权利要求1-9中任何一个所述的方法,其中,根据以下公式,方差σ2的最大似然估计通过关于σ2取对数似然函数的最大值确定:
11.如权利要求1-10中任何一个所述的方法,其中所述方差σ2根据以下公式确定:
12.如权利要求2-11中任何一个所述的方法,其中所述AIC根据以下公式取最小值:
AIC=2·[估计的模型的对数似然值]+2·[估计的参数的数目].
14.如权利要求13所述的方法,其中通过应用由最小化以下公式所产生的掩模
来设置矩阵元素为0:
如果
如果
从而得出最大似然估计
17.一种介质,其含有一个或多个利用上述权利要求所述的方法获得的基因之间网络关系的结果,所述结果存储在所述介质上。
18.一种用于确定网络关系的统计显著性的方法,包括:
(a)在每个时间点对于每个基因计算表达的平均对数比;
(b)从所有的测量中计算标准偏差;
(c)计算联合概率;及
(d)采用用于统计显著性的准则。
19.如权利要求18所述的方法,其中所述步骤(a)利用以下公式确定:
20.如权利要求18或者19所述的方法,其中所述步骤(b)利用以下公式确定:
其中xji[k]为对于基因j在时间点i测量k的数据值。
21.如权利要求18-20中任何一个所述的方法,其中就绝对值而言比测量值
xji大的
xj·的联合概率利用以下公式计算:
其中erf为误差函数。
22.如权利要求18-21中任何一个所述的方法,其中选择一个显著性水平α。
23.如权利要求18-22中任何一个所述的方法,其中如果Pi<α则舍弃无效假设。
24.如权利要求18-23中任何一个所述的方法,其中如果P<αn则舍弃无效假设,其中n为进行基因表达计算的时间点的数目。
25.一种用于确定网络关系的统计显著性的方法,包括:
(a)利用以下公式,在每个时间点对于每个基因计算表达测量的平均对数比:
(b)利用以下公式,计算所述测量的标准偏差:
其中
xji[k]为对于基因j在时间点i测量k的数据值。
(c)利用以下公式,计算就绝对值而言比测量值
xji大的
xj·的联合概率:
其中erf为误差函数;及
(d)应用用于统计显著性的准则以确定是否舍弃无效假设。
26.如权利要求25所述的方法,其中如果P<αn则舍弃无效假设,其中n为进行基因表达计算的时间点的数目。
27.一种用于推断如在这里充分描述的基因网络的方法。
28.一种用于确定如在这里充分描述的网络关系的统计显著性的方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US42882702P | 2002-11-25 | 2002-11-25 | |
US60/428,827 | 2002-11-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1717585A true CN1717585A (zh) | 2006-01-04 |
Family
ID=32393460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2003801040561A Pending CN1717585A (zh) | 2002-11-25 | 2003-11-25 | 利用微分方程从时序基因表达数据推断基因调控网络 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20040142362A1 (zh) |
EP (1) | EP1565741A4 (zh) |
JP (1) | JP2006507605A (zh) |
CN (1) | CN1717585A (zh) |
AU (1) | AU2003295842A1 (zh) |
CA (1) | CA2504856A1 (zh) |
WO (1) | WO2004048532A2 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646159A (zh) * | 2013-09-30 | 2014-03-19 | 温州大学 | 一种基于约束性布尔网络的最大评分预测方法 |
CN108491686A (zh) * | 2018-03-30 | 2018-09-04 | 中南大学 | 一种基于双向XGBoost的基因调控网络构建方法 |
CN109726352A (zh) * | 2018-12-12 | 2019-05-07 | 青岛大学 | 一种基于微分方程模型的基因调控网络的构建方法 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102004030296B4 (de) * | 2004-06-23 | 2008-03-06 | Siemens Ag | Verfahren zur Analyse eines regulatorischen genetischen Netzwerks einer Zelle |
JP2009169831A (ja) * | 2008-01-18 | 2009-07-30 | Mitsubishi Space Software Kk | 遺伝子相互作用データベース装置、遺伝子相互作用検索プログラムおよび遺伝子相互作用検索方法 |
AU2009307884B2 (en) | 2008-10-22 | 2014-07-31 | Merck Sharp & Dohme Corp. | Novel cyclic benzimidazole derivatives useful anti-diabetic agents |
US8329914B2 (en) | 2008-10-31 | 2012-12-11 | Merck Sharp & Dohme Corp | Cyclic benzimidazole derivatives useful as anti-diabetic agents |
WO2011106273A1 (en) | 2010-02-25 | 2011-09-01 | Merck Sharp & Dohme Corp. | Novel cyclic benzimidazole derivatives useful anti-diabetic agents |
US8796258B2 (en) | 2011-02-25 | 2014-08-05 | Merck Sharp & Dohme Corp. | Cyclic azabenzimidazole derivatives useful as anti-diabetic agents |
US20140045746A1 (en) | 2012-08-02 | 2014-02-13 | Merck Sharp & Dohme Corp. | Antidiabetic tricyclic compounds |
JP2016516004A (ja) | 2013-02-22 | 2016-06-02 | メルク・シャープ・アンド・ドーム・コーポレーションMerck Sharp & Dohme Corp. | 抗糖尿病二環式化合物 |
EP2970119B1 (en) | 2013-03-14 | 2021-11-03 | Merck Sharp & Dohme Corp. | Novel indole derivatives useful as anti-diabetic agents |
WO2015051496A1 (en) | 2013-10-08 | 2015-04-16 | Merck Sharp & Dohme Corp. | Antidiabetic tricyclic compounds |
EP3551176A4 (en) | 2016-12-06 | 2020-06-24 | Merck Sharp & Dohme Corp. | ANTIDIABETIC HETEROCYCLIC COMPOUNDS |
EP3558298A4 (en) | 2016-12-20 | 2020-08-05 | Merck Sharp & Dohme Corp. | ANTIDIABETIC SPIROCHROMAN COMPOUNDS |
EP3584727A4 (en) * | 2017-02-14 | 2020-03-04 | Fujifilm Corporation | METHOD AND DEVICE FOR ANALYZING BIOLOGICAL SUBSTANCE, AND PROGRAM |
CN113609652B (zh) * | 2021-07-14 | 2023-10-13 | 中国地质大学(武汉) | 分数阶环状基因调控网络的状态反馈控制方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030018457A1 (en) * | 2001-03-13 | 2003-01-23 | Lett Gregory Scott | Biological modeling utilizing image data |
CA2459570A1 (en) * | 2001-09-05 | 2003-03-13 | Genicon Sciences Corporation | Apparatus for reading signals generated from resonance light scattered particle labels |
US20030144823A1 (en) * | 2001-11-01 | 2003-07-31 | Fox Jeffrey J. | Scale-free network inference methods |
US7415359B2 (en) * | 2001-11-02 | 2008-08-19 | Gene Network Sciences, Inc. | Methods and systems for the identification of components of mammalian biochemical networks as targets for therapeutic agents |
-
2003
- 2003-11-25 JP JP2004555640A patent/JP2006507605A/ja active Pending
- 2003-11-25 US US10/722,033 patent/US20040142362A1/en not_active Abandoned
- 2003-11-25 CN CNA2003801040561A patent/CN1717585A/zh active Pending
- 2003-11-25 AU AU2003295842A patent/AU2003295842A1/en not_active Abandoned
- 2003-11-25 CA CA002504856A patent/CA2504856A1/en not_active Abandoned
- 2003-11-25 WO PCT/US2003/037475 patent/WO2004048532A2/en active Application Filing
- 2003-11-25 EP EP03787054A patent/EP1565741A4/en not_active Withdrawn
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646159A (zh) * | 2013-09-30 | 2014-03-19 | 温州大学 | 一种基于约束性布尔网络的最大评分预测方法 |
CN103646159B (zh) * | 2013-09-30 | 2016-07-06 | 温州大学 | 一种基于约束性布尔网络的最大评分预测方法 |
CN108491686A (zh) * | 2018-03-30 | 2018-09-04 | 中南大学 | 一种基于双向XGBoost的基因调控网络构建方法 |
CN108491686B (zh) * | 2018-03-30 | 2021-06-18 | 中南大学 | 一种基于双向XGBoost的基因调控网络构建方法 |
CN109726352A (zh) * | 2018-12-12 | 2019-05-07 | 青岛大学 | 一种基于微分方程模型的基因调控网络的构建方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2004048532A2 (en) | 2004-06-10 |
US20040142362A1 (en) | 2004-07-22 |
JP2006507605A (ja) | 2006-03-02 |
EP1565741A2 (en) | 2005-08-24 |
AU2003295842A1 (en) | 2004-06-18 |
EP1565741A4 (en) | 2008-04-02 |
WO2004048532A3 (en) | 2004-09-30 |
CA2504856A1 (en) | 2004-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1717585A (zh) | 利用微分方程从时序基因表达数据推断基因调控网络 | |
Celton et al. | Comparative analysis of missing value imputation methods to improve clustering and interpretation of microarray experiments | |
MacIsaac et al. | Practical strategies for discovering regulatory DNA sequence motifs | |
Quackenbush | Computational analysis of microarray data | |
Dondrup et al. | EMMA: a platform for consistent storage and efficient analysis of microarray data | |
Yu et al. | Incorporating nonlinear relationships in microarray missing value imputation | |
Kim et al. | Microbial forensics: predicting phenotypic characteristics and environmental conditions from large-scale gene expression profiles | |
Žitnik et al. | Gene prioritization by compressive data fusion and chaining | |
Wang et al. | AC-PCoA: Adjustment for confounding factors using principal coordinate analysis | |
Yoosefzadeh-Najafabadi et al. | Genome-wide association study statistical models: A review | |
Gerber et al. | Automated discovery of functional generality of human gene expression programs | |
Lopes et al. | Entropic biological score: a cell cycle investigation for GRNs inference | |
Mosleth et al. | Analysis of megavariate data in functional omics | |
Huang et al. | Computational strategies for the identification of a transcriptional biomarker panel to sense cellular growth states in Bacillus subtilis | |
Martínez | Time course gene expression experiments | |
Zeng et al. | A link-free sparse group variable selection method for single-index model | |
Srivastava et al. | Genome-wide functional annotation by integrating multiple microarray datasets using meta-analysis | |
Singh et al. | Accuracy of functional gene community detection in Saccharomyces cerevisiae by maximizing Generalized Modularity Density | |
CN1714371A (zh) | 根据时间序列基因表达数据的基因网络的非线性模拟 | |
Arani et al. | Genetic variant effect prediction by supervised nonnegative matrix tri-factorization | |
Samee et al. | K4. Gene network construction and pathways analysis for high throughput microarrays | |
Carey et al. | A Big Data Pipeline: Identifying Dynamic Gene Regulatory Networks from Time Course GEO Data with Applications to Influenza Infection | |
Zhao et al. | Tutorial: guidelines for survival analysis with omics data | |
Kamgnia Wonkap | Gene Regulatory Network Inference Using Machine Learning Techniques | |
Khachatryan et al. | Reference-free resolution of long-read metagenomic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1084149 Country of ref document: HK |
|
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: WD Ref document number: 1084149 Country of ref document: HK |