CN114609318A - 一种基于分子结构关联网络的规模化代谢组定性方法 - Google Patents
一种基于分子结构关联网络的规模化代谢组定性方法 Download PDFInfo
- Publication number
- CN114609318A CN114609318A CN202011407875.5A CN202011407875A CN114609318A CN 114609318 A CN114609318 A CN 114609318A CN 202011407875 A CN202011407875 A CN 202011407875A CN 114609318 A CN114609318 A CN 114609318A
- Authority
- CN
- China
- Prior art keywords
- metabolites
- metabolite
- metabolome
- molecular structure
- network
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000002207 metabolite Substances 0.000 claims abstract description 177
- 230000014759 maintenance of location Effects 0.000 claims abstract description 37
- 239000000284 extract Substances 0.000 claims abstract description 15
- 239000012472 biological sample Substances 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 150000002500 ions Chemical class 0.000 claims description 43
- 238000001819 mass spectrum Methods 0.000 claims description 35
- 239000000523 sample Substances 0.000 claims description 22
- 238000013051 Liquid chromatography–high-resolution mass spectrometry Methods 0.000 claims description 16
- 238000002474 experimental method Methods 0.000 claims description 7
- 238000012797 qualification Methods 0.000 claims description 7
- 238000000547 structure data Methods 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000004811 liquid chromatography Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- 238000007637 random forest analysis Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000013136 deep learning model Methods 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 230000002503 metabolic effect Effects 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 238000010801 machine learning Methods 0.000 claims 1
- 238000011173 large scale experimental method Methods 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 12
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 9
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 9
- 238000004949 mass spectrometry Methods 0.000 description 7
- 238000004885 tandem mass spectrometry Methods 0.000 description 7
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 210000004185 liver Anatomy 0.000 description 6
- 238000002705 metabolomic analysis Methods 0.000 description 6
- 230000001431 metabolomic effect Effects 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 241001465754 Metazoa Species 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 5
- XPFVYQJUAUNWIW-UHFFFAOYSA-N furfuryl alcohol Natural products OCC1=CC=CO1 XPFVYQJUAUNWIW-UHFFFAOYSA-N 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- BDAGIHXWWSANSR-UHFFFAOYSA-N methanoic acid Natural products OC=O BDAGIHXWWSANSR-UHFFFAOYSA-N 0.000 description 4
- 239000006228 supernatant Substances 0.000 description 4
- AYWGPQFATRWJMF-UHFFFAOYSA-N dihydroxy-beta-ionone Natural products CC(=O)C=CC1(O)C(C)(C)CCCC1(C)O AYWGPQFATRWJMF-UHFFFAOYSA-N 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004896 high resolution mass spectrometry Methods 0.000 description 3
- 210000005228 liver tissue Anatomy 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 239000003643 water by type Substances 0.000 description 3
- CHBOSHOWERDCMH-UHFFFAOYSA-N 1-chloro-2,2-bis(4-chlorophenyl)ethane Chemical compound C=1C=C(Cl)C=CC=1C(CCl)C1=CC=C(Cl)C=C1 CHBOSHOWERDCMH-UHFFFAOYSA-N 0.000 description 2
- NCYCYZXNIZJOKI-IOUUIBBYSA-N 11-cis-retinal Chemical compound O=C/C=C(\C)/C=C\C=C(/C)\C=C\C1=C(C)CCCC1(C)C NCYCYZXNIZJOKI-IOUUIBBYSA-N 0.000 description 2
- VHTFHZGAMYUZEP-UHFFFAOYSA-N 2,6,6-Trimethyl-1-cyclohexen-1-acetaldehyde Chemical compound CC1=C(CC=O)C(C)(C)CCC1 VHTFHZGAMYUZEP-UHFFFAOYSA-N 0.000 description 2
- OSWFIVFLDKOXQC-UHFFFAOYSA-N 4-(3-methoxyphenyl)aniline Chemical compound COC1=CC=CC(C=2C=CC(N)=CC=2)=C1 OSWFIVFLDKOXQC-UHFFFAOYSA-N 0.000 description 2
- KDCGOANMDULRCW-UHFFFAOYSA-N 7H-purine Chemical compound N1=CNC2=NC=NC2=C1 KDCGOANMDULRCW-UHFFFAOYSA-N 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- DRTQHJPVMGBUCF-XVFCMESISA-N Uridine Chemical compound O[C@@H]1[C@H](O)[C@@H](CO)O[C@H]1N1C(=O)NC(=O)C=C1 DRTQHJPVMGBUCF-XVFCMESISA-N 0.000 description 2
- UZFLPKAIBPNNCA-BQYQJAHWSA-N alpha-ionone Chemical compound CC(=O)\C=C\C1C(C)=CCCC1(C)C UZFLPKAIBPNNCA-BQYQJAHWSA-N 0.000 description 2
- UZFLPKAIBPNNCA-UHFFFAOYSA-N alpha-ionone Natural products CC(=O)C=CC1C(C)=CCCC1(C)C UZFLPKAIBPNNCA-UHFFFAOYSA-N 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- SQFRYZPEWOZAKJ-UHFFFAOYSA-N dihydrodehydro-beta-ionone Chemical compound CC(=O)CCC1=C(C)C=CCC1(C)C SQFRYZPEWOZAKJ-UHFFFAOYSA-N 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000010828 elution Methods 0.000 description 2
- 235000019253 formic acid Nutrition 0.000 description 2
- 238000004108 freeze drying Methods 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 150000003384 small molecules Chemical class 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000001195 ultra high performance liquid chromatography Methods 0.000 description 2
- PSQYTAPXSHCGMF-BQYQJAHWSA-N β-ionone Chemical compound CC(=O)\C=C\C1=C(C)CCCC1(C)C PSQYTAPXSHCGMF-BQYQJAHWSA-N 0.000 description 2
- SFEOKXHPFMOVRM-UHFFFAOYSA-N (+)-(S)-gamma-ionone Natural products CC(=O)C=CC1C(=C)CCCC1(C)C SFEOKXHPFMOVRM-UHFFFAOYSA-N 0.000 description 1
- PCSKKIUURRTAEM-UHFFFAOYSA-N 5-hydroxymethyl-2-furoic acid Chemical group OCC1=CC=C(C(O)=O)O1 PCSKKIUURRTAEM-UHFFFAOYSA-N 0.000 description 1
- -1 5-hydroxymethylfuran aldehyde Chemical class 0.000 description 1
- QGZKDVFQNNGYKY-UHFFFAOYSA-O Ammonium Chemical compound [NH4+] QGZKDVFQNNGYKY-UHFFFAOYSA-O 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 125000003158 alcohol group Chemical group 0.000 description 1
- 125000003172 aldehyde group Chemical group 0.000 description 1
- SHGAZHPCJJPHSC-YCNIQYBTSA-N all-trans-retinoic acid Chemical compound OC(=O)\C=C(/C)\C=C\C=C(/C)\C=C\C1=C(C)CCCC1(C)C SHGAZHPCJJPHSC-YCNIQYBTSA-N 0.000 description 1
- KKTCWAXMXADOBB-UHFFFAOYSA-N azanium;hydrogen carbonate;hydrate Chemical compound [NH4+].O.OC([O-])=O KKTCWAXMXADOBB-UHFFFAOYSA-N 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 238000010296 bead milling Methods 0.000 description 1
- DRTQHJPVMGBUCF-PSQAKQOGSA-N beta-L-uridine Natural products O[C@H]1[C@@H](O)[C@H](CO)O[C@@H]1N1C(=O)NC(=O)C=C1 DRTQHJPVMGBUCF-PSQAKQOGSA-N 0.000 description 1
- 230000008827 biological function Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 125000003178 carboxy group Chemical group [H]OC(*)=O 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000007876 drug discovery Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002546 full scan Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 102000004196 processed proteins & peptides Human genes 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 229930002330 retinoic acid Natural products 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 125000001424 substituent group Chemical group 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 229960001727 tretinoin Drugs 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 238000004704 ultra performance liquid chromatography Methods 0.000 description 1
- DRTQHJPVMGBUCF-UHFFFAOYSA-N uracil arabinoside Natural products OC1C(O)C(CO)OC1N1C(=O)NC(=O)C=C1 DRTQHJPVMGBUCF-UHFFFAOYSA-N 0.000 description 1
- 229940045145 uridine Drugs 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
- G01N30/8686—Fingerprinting, e.g. without prior knowledge of the sample components
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/90—Programming languages; Computing architectures; Database systems; Data warehousing
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Analytical Chemistry (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Library & Information Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于分子结构关联网络的规模化代谢组定性方法。首先,收集开源代谢组数据库中的内源性代谢物,基于代谢物分子结构相似性构建代谢组水平分子结构关联网络;其次,对生物样本提取物进行非靶向代谢组学分析,并构建保留时间预测模型;进一步通过少量结构确证的代谢物作为种子代谢物,基于网络相邻代谢物与种子代谢物有相似MS2这一前提,对代谢物进行定性。本发明方法不依赖于大规模实验MS2数据库,即可真正实现代谢组规模的快速定性,且定性结果更可靠。
Description
技术领域
本发明涉及分析化学和代谢组学领域,是一种基于分子结构关联网络的规模化代谢组定性方法。
研究背景
代谢组学(Nicholson,JK et al.,Metabonomics:a platform for studyingdrug toxicity and gene function.Nature Reviews Drug Discovery 2002,1,153-161.)是以生物体中分子量小于1000的小分子代谢物为研究对象,运用多种分析手段如高分辨质谱(HRMS)、核磁共振(NMR)、色谱-质谱联用等从整体水平上研究体内代谢物组成及其与生理病理相关的变化规律。代谢组学已被广泛应用于生命科学、疾病诊断、药物安全性评价、药理研究、营养科学等多个领域。
代谢物结构鉴定是代谢组学发展的瓶颈问题之一。生物体内代谢物成分及其复杂、理化性质各异、含量差别巨大,生物样品中的代谢物浓度分布范围达7~9个数量级。随着分析技术的不断进步,特别是高分辨质谱技术迅猛发展,生物样本非靶向代谢组学分析获得了海量的代谢组数据,但其中可被鉴定的代谢物不超过25%(Baker,M.,Metabolomics:from small molecules to big ideas.Nature Methods 2011,8,117-121.)。质谱数据库检索是代谢物鉴定的最主要方式,其收录情况直接决定了代谢物鉴定的数量与质量。Metlin数据库目前已收录超过1000,000种化合物,其中50万多种有二级质谱(MS/MS)数据;人类代谢组数据库(HMDB)(http://www.hmdb.ca/)收录114,260种化合物,仅有7,418多个化合物有实测MS/MS谱图数据,279,972种为预测MS/MS数据。预测二级谱图是解决二级短缺问题的一种手段,但大多预测谱图与实验谱图相差很大,仍然会给代谢物的定性带来困难。基于实验二级质谱相似性的分子网络如GNPS(https://gnps.ucsd.edu/)已被证明是一种有效的代谢物辅助定性工具。其基于具有相似二级谱图的代谢物其结构具有相似性的假设,但采用该方法对代谢物结构鉴定仍需依赖大规模的实验二级谱图数据库,因此其定性能力仍受到目前代谢组数据库实验二级谱图短缺的制约。
目前代谢组学研究主要基于其中的已知代谢物,大量稳定、可重复检测的“未知代谢物”虽可能具有重要生物功能,却因无法鉴定分子结构,未能对其利用或开展深入研究,导致质谱信息利用率低。因此,亟待发展功能代谢组学新技术、新方法,包括更高效可靠的算法、知识库及相关工具来实现对非靶向超高效液相色谱-高分辨质谱(UHPLC-HRMS)代谢谱数据的深度解析。因此,发展能够辅助代谢物结构解析,尤其是能够鉴定未知代谢物的分析新方法显得尤为迫切。
为此,本发明充分利用代谢组数据库中丰富的一级及其结构信息,发展一种不依赖大规模二级数据库的基于分子结构关联网络的规模化代谢组定性方法,突破保留时间和二级信息缺乏造成的定性瓶颈。
发明内容
本发明提供了一种基于分子结构关联网络的代谢组规模化定性方法。为了实现本发明目的,收集开源代谢组数据库中的内源性代谢物,基于所获取的代谢物分子指纹相似性构建代谢物分子结构关联网络;对生物样本提取物进行基于超高效液相色谱-高分辨质谱的非靶向代谢组学分析,获取生物样本提取物的代谢组相关色谱-质谱信息;基于标样的分子描述符和实验保留时间,构建保留时间预测模型;利用保留时间预测模型计算分子结构关联网络中每个节点的预测保留时间;利用非靶向超高效液相色谱-高分辨质谱代谢组实验数据,以分子结构关联网络作为背景网络,进行代谢组规模化定性。本发明采用的技术方案步骤如下:
第一步,构建代谢组的分子结构关联网络;根据开源代谢组数据库如人类代谢组HMDB数据库、京都基因与基因组百科全书KEGG数据库获取内源性代谢物,包括其名称,单同位素质量,分子结构数据文件;基于其分子结构数据文件获得每个代谢物的分子指纹,这里分子指纹可以是圆形指纹、MACCS指纹、原子对指纹和拓扑指纹中的任意一种;计算任意两个候选代谢物分子指纹间的相似度;设定相似性阈值,以代谢物为节点,分子指纹相似性为边,构建代谢组水平的分子结构关联网络;
第二步,采用超高效液相色谱-高分辨质谱对待测生物样本的提取物进行非靶向代谢组学分析;获取提取物代谢组的相关色谱-质谱信息,包括实验测得的代谢物峰的保留时间tR实测,一级质谱信息即一级离子质荷比m/z实测,和相应的二级质谱信息即二级离子的质荷比和强度;一级离子指化合物被离子化后直接采集到的离子;二级离子指一级离子在施加一定能量碰撞碎裂后采集到的离子;
第三步,构建保留时间预测模型;采用与生物提取物相同的超高效液相色谱-高分辨质谱数据采集条件,分析内源性代谢物的混合标样,获取其液相色谱实验保留时间;采用开源网站ChemDes(http://www.scbdd.com/chemdes)和Dragon7(https://chm.kode-solutions.net/index.php)获取标样的分子描述符;基于实验保留时间和分子描述符,通过定量结构-保留关系构建保留时间预测模型;保留时间预测模型采用多重线性回归模型、弹性网络回归模型、随机森林模型、支持向量机模型和深度学习模型中任意一种方法建立;
第四步,计算分子结构关联网络中代谢物的预测保留时间;基于其分子结构数据文件获得每个代谢物的分子描述符;利用保留时间预测模型计算分子结构关联网络中每个节点的预测保留时间;
第五步,基于分子结构关联网络对非靶向代谢组学数据进行规模化定性;以第一步构建的分子结构关联网络作为背景网络,利用标准样品从非靶向超高效液相色谱-高分辨质谱代谢组实验数据鉴定50~100个代谢物作为种子代谢物,将其映射到已建立的分子结构关联网络中,从网络中获取种子代谢物的相邻代谢物,这里相邻代谢物指分子结构关联网络中有直接边连接的代谢物;将种子代谢物的二级质谱作为相邻代谢代谢物二级质谱,设定搜索阈值,在实验数据中搜索与相邻代谢物m/z理论,tR预测,二级质谱匹配的代谢物峰,匹配成功则该代谢物峰鉴定完成;当有多个匹配结果时,对匹配结果进行打分;鉴定出的代谢物再作为新的种子,重复定性过程,直至没有新的代谢物被定性出来。
本发明依据结构相似的代谢物,其MS/MS具有相似性这一前提,建立了代谢组数据库导向的基于分子结构关联网络的规模化定性方法,实现对未知代谢物的结构鉴定。通过建立代谢物数据库及其代谢物分子结构关联网络,采用分子结构关联网络去鉴定没有标准MS/MS谱图的代谢物,使得代谢物的结构鉴定可不依赖大规模的标准MS/MS数据库。本发明是一种不依赖于大规模实验二级谱图数据库的代谢组深度注释方法,可实现代谢组规模化、可靠定性,显著扩大代谢组注释的覆盖度。
附图说明
图1正离子模式分子结构关联网络及其局部放大示意图;
图2基于分子结构关联网络的代谢物定性过程示意图;
图3正负离子模式下小鼠肝脏的鉴定结果(共鉴定518个代谢物);
图4分子结构关联网络定性方法与现有方法对小鼠肝脏组织代谢物的鉴定结果比较。
具体实施方式
下面结合附图对本发明实施作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例1
为了证实本发明的有效性和可行性,在组织提取物中添加由885个代谢物(包括尿嘧啶核苷和嘌呤等代谢物,添加终浓度为0.1~5微克/毫升)组成的混合标样,混合标样为从构建的分子结构网络中选取的,有市售标准品的代谢物组成。
组织代谢组的提取:采用动物代谢组学方法,提取小鼠肝脏中的代谢物。首先称取10毫克小鼠肝脏组织于1.5毫升离心管中,加入600微升的甲醇/水(体积比,4:1)提取剂(含由885个代谢物组成的混合标样),研磨珠打进行匀浆处理(25Hz×1min*2),在4℃,12000rpm条件下离心15min。取400μL上清液,放入真空离心浓缩仪冻干。在冻干样品粉末中加入50微升乙腈/水(体积比,1:3),涡旋1min,在高速离心机中采用4摄氏度,14000rpm转速,离心10分钟。
将加标后的生物样本提取物进行超高效液相色谱-高分辨质谱数据采集,以采集的非靶向代谢组学数据中混合标样的定性为例进行发明的原理性说明。
基于代谢组数据库的分子结构关联网络构建:分子结构关联网络中的代谢物来自开源代谢组数据库人类代谢组数据库(HMDB)(https://hmdb.ca/)和京都基因与基因组百科全书数据库(KEGG)(https://www.genome.jp/kegg/)。以正离子模式为例,初始代谢物源自HMDB数据库收录的5170个代谢物以及KEGG数据库收录的5102个代谢物,并从上述数据库中获取每个代谢物的名称,单同位素质量及分子结构数据文件。由分子结构数据文件计算每个代谢物的分子指纹,本例中采用圆形指纹,计算任意两个分子圆形指纹间的相似度(相似度计算方法来源于开源工具RDKit)。设置分子指纹相似度阈值为0.6,符合该阈值的代谢物共有8562个(4141个代谢物来自HMDB,4421个代谢物来自KEGG)。以代谢物为节点,分子间指纹相似性为边,构建分子结构关联网络(共有8,562个节点和264,095条边),正离子模式下的分子结构关联网络见图1(左图)。以代谢物HMDB0002432(HMDB数据库ID)为例,其名称为5-羟甲基糠酸(5-Hydroxymethyl-2-Furoate),其分子结构见图1(右下图);它的电喷雾质谱正离子模式下[M+H]+的m/z理论,tR预测分别为m/z 143.0344,1.36min。将代谢物HMDB0002432及其周边局部放大见图1(右上图),分子结构关联网络中有直接边连接的代谢物称之为相邻代谢物,HMDB0002432的两个相邻代谢物分别为5-羟甲基呋喃醛(5-Hydroxymethylfurfural,KEGG数据库ID:C11101)和5-羟甲基呋喃醇(5-Hydroxymethylfurfuralalcohol,KEGG数据库ID:C20443)。从图中三个代谢物的分子结构可知,它们均具有5-羟甲基呋喃子结构,仅在2位的取代基有差异,分别为醛基,醇基和羧基。因此,其分子结构均有较好的相似性,其中HMDB0002432与C11101和C20443的指纹相似度均为0.65,C11101和C20443的指纹相似度为0.64,均满足相似度阈值,故在网络中有直接连边。
实验非靶向色谱-质谱信息获取:数据采集于ACQUITY UHPLC超高效液相色谱系统(UPLC,Waters,Milford,MA,USA.)和Q Exactive HF高分辨质谱(Thermo FisherScientific,Rockford,IL,USA.)联用的分析仪器。
正离子模式下液相色谱条件为:A相和B相分别为0.1%的甲酸/水(体积比)和0.1%的甲酸/乙腈(体积比)。流速为0.35mL/min。总分析时间为30min。洗脱梯度以5%B开始,保持1min。然后在23min内线性增加到100%B,并保持4min。接着在0.1min内线性回到初始比例,并保持2.9min。样品采用ACQUITY BEH C8色谱柱(100mm x 2.1mm,1.7μm,Waters,Milford,MA,U.S.A.)进行分离。柱温为50℃。样品进样室温度设置为4℃,样品进样量为5μL。
负离子模式下液相色谱条件为:A相和B相分别为6.5mM碳酸氢铵-水和6.5mM碳酸氢铵-95%甲醇/水(体积比)。流速为0.35mL/min。总分析时间为25min。洗脱梯度以2%B开始,1min后开始线性增长,在18min时到100%B,并保持到22min。接着以线性方式在22.1min时回到初始比例,并保持到25min。样品采用ACQUITY HSS T3色谱柱(100mm x 2.1mm,1.8μm,Waters,Milford,MA,U.S.A.)进行分离。柱温为50℃。样品进样室温度设置为4℃,样品进样量为5μL。
Q Exactive HF质谱分析条件为:扫描模式为全扫加二级质谱采集模式(full MS/data-dependent(ddMS2))。在全扫质谱设置中,分辨率采用120,000。自动增益控制目标(AGC target)和最大注入时间(maximum IT)分别设置为3×106离子容量和100ms。全扫质量扫描范围为m/z 85-1250。在二级质谱设置中,自动增益控制目标(AGC target)和最大注入时间(maximum IT)分别设置为1×105离子容量和50ms。隔离窗口为m/z 1.0。碰撞能采用混合归一化能量15%,30%和45%。二级的采集由每个全扫描循环中响应最强的前10个离子触发。加入Inclusion list并设置为on。正、负离子模式下喷雾电压分别为3.5kV和3.0kV,离子传输管的温度为320℃,辅助气的加热器温度为350℃。鞘气和辅助气流速分别为45和10(in arbitrary units)。S-lens设置为50.0(in arbitrary units)。
在采集到的非靶向代谢组学数据中,提取出885个代谢物的实验保留时间tR实测,一级质谱信息即一级离子质荷比m/z实测和相应的二级质谱信息即二级离子的质荷比和强度。
保留时间预测模型构建:采用与动物提取物相同的超高效液相色谱-高分辨质谱数据采集条件,分析1551个标样,获取液相色谱实验保留时间。利用标样的SDF文件在开源网站ChemDes(http://www.scbdd.com/chemdes)中计算获得每个标样的1D&2D分子描述符,采用随机森林方法,构建保留时间预测模型。
基于分子结构关联网络定性:以构建的分子结构关联网络作为背景网络,对非靶向超高效液相色谱-高分辨质谱代谢组采集的加标代谢物进行鉴定。以正离子模式为例,定性的具体过程如下:
1)利用标准样品从加标的组织提取物非靶向超高效液相色谱-高分辨质谱代谢组实验数据中鉴定出60个代谢物作为种子代谢物,将其映射到已建立的分子结构关联网络中,从网络中获取种子代谢物相邻代谢物。如图2所示,种子代谢物HMDB0036565(β-紫罗兰酮,beta-Ionone)的相邻代谢物共有4个,分别为HMDB0032913(二羟基-β-紫罗兰酮,dihydro-β-ionone)、C02110(视黄酸,9-cis-Retinoic acid)、HMDB0059883(α-紫罗兰酮,alpha-Ionone)和HMDB0041011(β-环高柠檬醛,beta-Cyclocitral)。
2)将种子代谢物的二级质谱赋给相邻代谢代谢物,作为其拟二级质谱。设定搜索阈值,实验数据中代谢物峰的一级离子质荷比m/z实测,实验保留时间tR实测与相邻代谢物m/z理论,tR预测匹配的阈值同时满足|tR预测-tR实测|/tR预测*100%<30%与|m/z理论-m/z实测|/m/z理论*1000000<10ppm;且其实验二级质谱和相邻代谢物的拟二级质谱的相似度≥0.5。在实验数据中搜索与相邻代谢物m/z理论,tR预测以及拟二级质谱匹配的代谢物峰,匹配成功则该代谢物峰鉴定完成。如图2左下图所示,种子代谢物的二级谱图为下方红色谱图,将其作为4个相邻代谢物的二级谱图,每个相邻代谢物分别在实验数据中寻找与其m/z理论,tR预测,以及拟二级质谱匹配的代谢物峰。相邻代谢物HMDB0032913,C02110,HMDB0059883和HMDB0041011的m/z理论,tR预测分别为m/z 195.1749,12.69min;m/z 285.2218,17.11min;m/z 193.1592,13.38min和m/z 153.1279,9.42min。在实验数据中搜索与相邻代谢物满足阈值要求的代谢物峰,如图2左下图所示,在实验数据中搜索到保留时间14.32min,m/z 195.1745的代谢物峰,与HMDB0032913的|tR预测-tR实测|/tR预测*100%=12.84%,|m/z理论-m/z实测|/m/z理论*1000000=2.05ppm,且该峰的实验二级谱图(蓝色谱图)与HMDB0032913的拟二级谱图(红色谱图)的相似度为0.80,则该代谢物峰被定性为二羟基-β-紫罗兰酮(HMDB0032913)。为了进一步证实采用分子结构网络及相邻代谢物的“拟二级谱图”方法定性的可行性和可靠性,在图2右下图中对比了标准样品HMDB0032913的实验二级谱图(红色谱图)与混合标准样品中代谢物峰(14.32min,m/z 195.1745)的实验二级谱图(蓝色谱图),二者的二级谱图有很高的相似性,证明利用相邻代谢物拟二级谱图的分子结构关联网络定性方法是准确可靠的。采用上述定性方法,实验数据中2个代谢物峰(m/z实测,tR实测,二级相似度)m/z 193.1590,13.39min,0.96和m/z 153.1268,11.08min,0.54分别和α-紫罗兰酮(HMDB0059883)和β-环高柠檬醛(HMDB0041011)匹配,这2个代谢物峰也被成功鉴定。
3)当实验数据搜索出多个匹配结果时,对匹配结果进行打分,得分规则为:
得分=0.25*(1-(|m/z理论-m/z实测|/m/z理论*1000000)/10)+0.25*(1-(|tR预测-tR实测|/tR预测*100%)/30%)+0.5*二级谱图相似度
其中,m/z理论为利用代谢物的分子式所获得的理论一级离子质荷比;m/z实测为实验测得的一级离子质荷比;tR预测为利用前述构建的保留时间预测模型所预测出的保留时间;tR实测为实验测得的代谢物峰的保留时间。
如实验数据中搜索到3个代谢物峰满足搜索阈值,与HMDB0059883匹配,其m/z实测,tR实测,二级质谱相似度分别为m/z193.1581,13.38min,0.96;m/z 193.1583,13.49min,0.99和m/z 193.1590,13.39min,0.96,对3个结果进行打分,对应的得分为0.87,0.89和0.96,按照从大到小排序。得分越高的代谢物峰被鉴定的可靠性越高,但有多个搜索结果的代谢物峰不再作为新的种子参与下一轮定性。
4)鉴定出的代谢物再作为新的种子,重复定性过程,直至没有新的代谢物被定性出来。如实验数据中代谢物峰被成功鉴定为二羟基-β-紫罗兰酮(图2中HMDB0032913),代谢物峰(m/z 195.1745,14.32min)的相应二级谱图被赋予下一级相邻代谢物HMDB0037139(二氢脱氢β-紫罗兰酮,Dihydrodehydro-beta-ionone),作为其拟二级谱图。HMDB0037139的m/z理论,tR预测为193.1592,11.94min。在实验数据中寻找到满足阈值的代谢物峰m/z 193.1590,13.39min,其实验二级谱图与拟二级谱图相似度为0.72,匹配成功。代谢物峰(m/z193.1590,13.39min)被成功鉴定为二氢脱氢β-紫罗兰酮,其作为新的种子重复上述定性过程。
正离子模式下,885个代谢物混合标样的实验代谢物峰均符合一级质谱离子质荷比和保留时间的搜索阈值,从885个代谢物峰中选取60个代谢物作为种子代谢物映射到分子结构关联网络,采用上述鉴定过程,共可关联出578个代谢物峰,网络关联覆盖度为72.1%(638/885)。网络中未关联出的代谢物峰主要因为混合标样较少,在网络中缺少有效边连接。关联出的代谢物峰可满足全部鉴定阈值可成功鉴定的为461个代谢物峰(含种子),鉴定率为72.2%(461/638),其中392个代谢物峰鉴定结果正确,鉴定正确率为85.0%(392/461)。
上述结果显示,本发明无需依赖于大规模实验二级谱图数据库的代谢组定性方法,可实现代谢组规模化快速定性,且结果准确可靠。
实施例2
为了进一步证实本发明在生物样本的提取物中同样有效可行,我们对动物组织(小鼠肝脏)进行了代谢组的提取,并对动物组织提取物进行超高效液相色谱-高分辨质谱数据采集,以采集的非靶向代谢组学数据的定性为例进行说明。
组织代谢组的提取:采用动物代谢组学方法,提取小鼠肝脏中的代谢物。首先称取10毫克小鼠肝脏组织于1.5毫升离心管中,加入600微升的甲醇/水(体积比,4:1)提取剂,研磨珠打进行匀浆处理(25Hz×1min*2),在4℃,12000rpm条件下离心15min。取400μL上清液,放入真空离心浓缩仪冻干。在冻干样品粉末中加入50微升乙腈/水(体积比,1:3),涡旋1min,在高速离心机中采用4摄氏度,14000rpm转速,离心10分钟。
过程和条件同实施例1,与其不同之处在于:
非靶向代谢组学数据采集:同实施例1。
实验色谱-质谱信息获取:基于小鼠肝脏提取物的非靶向代谢组学数据,采用软件CompoundDiscovery3.1获得峰表,包括实验保留时间tR实测,一级质谱信息即一级离子质荷比m/z实测,导出Excel表。采用软件ProteoWizard将原始数据进行转化,获得.mgf的二级文件,包含相应的二级质谱信息即二级离子的质荷比和强度。
保留时间预测模型构建和分子结构关联网络构建过程同实例1所示,这里不再赘述。
基于分子结构关联网络定性:以构建的分子结构关联网络作为背景网络,对非靶向超高效液相色谱-高分辨质谱代谢组采集的实验数据进行鉴定,确定待测生物样本中的代谢物。正离子模式下,512个代谢物峰在分子结构关联网络中符合一级质谱离子质荷比和保留时间的搜索阈值。采用同实施例1相同的鉴定过程,利用60个代谢物作为种子代谢物,共关联到422个代谢物峰,正离子模式的关联覆盖率为94.1%,其中成功鉴定353个代谢物峰;负离子模式下,342个代谢物峰在分子结构关联网络中。采用同实施例1相同的鉴定过程,利用55个代谢物作为种子代谢物,共关联出258个代谢物峰,负离子模式的关联覆盖率为91.5%,其中成功鉴定230个代谢物峰。见图3,正负离子模式下共成功鉴定出518个代谢物。
将分子结构关联网络定性方法与传统的搜库鉴定方法比较,以小鼠肝脏为例,采用Compound Discoverer(CD)质谱数据库,基于一级离子质荷比和二级质谱信息共鉴定到139个代谢物。分子结构关联网络定性方法基于保留时间,一级离子质荷比和二级相似性共鉴定出518个代谢物。如图4,本方法鉴定到更多的代谢物,且对传统方法所鉴定的代谢物有较好的包容性。本方法充分利用代谢组数据库中丰富的代谢物一级离子质荷比及其结构信息,发展不依赖大规模MS/MS数据库的基于分子结构关联网络的规模化代谢组定性方法,突破保留时间和二级质谱信息缺乏造成的定性瓶颈。从而,达到大规模、准确定性的目的。
Claims (9)
1.一种基于分子结构关联网络的规模化的代谢组定性方法,其特征在于:
第一步,构建代谢组分子结构关联网络;从开源代谢组数据库获取内源性代谢物,包括其名称,单同位素质量,分子结构数据文件;基于分子结构数据文件获得每个代谢物的分子指纹,并计算任意两个代谢物分子指纹间相似度;设定分子指纹间相似度阈值0.6-0.9,以代谢物为节点,分子指纹相似性为边,将大于等于分子指纹间相似度域值的代谢物间连线,构建分子结构关联网络;
第二步,获取待测生物样本提取物中代谢组的色谱-质谱信息;采用超高效液相色谱-高分辨质谱对生物样本提取物进行非靶向代谢组学分析,获取包括实验测得的代谢物峰的保留时间tR实测,一级质谱信息即一级离子质荷比m/z实测,和相应的二级质谱信息即二级离子的质荷比和强度;
第三步,构建保留时间预测模型;采用与生物提取物相同的超高效液相色谱-高分辨质谱数据采集条件,分析已知内源性代谢物的混合标样,获取其液相色谱实验保留时间;进而基于标样的分子描述符和实验保留时间,采用机器学习算法,构建保留时间预测模型;
第四步,利用保留时间预测模型计算分子结构关联网络中每个节点的预测保留时间;
第五步,基于分子结构关联网络对非靶向代谢组学数据进行规模化定性;以第一步构建的分子结构关联网络作为背景网络,对非靶向超高效液相色谱-高分辨质谱代谢组采集的色谱峰进行鉴定,确定待测生物样本中的代谢物。
2.根据权利要求1所述方法,其特征在于:第一步所述开源代谢组数据库为人类代谢组HMDB数据库和/或京都基因与基因组百科全书KEGG数据库。
3.根据权利要求1所述方法,其特征在于:第一步所述内源性代谢物为人源、鼠源中一种或二种以上内源性代谢物。
4.根据权利要求1所述方法,其特征在于:第一步所述分子指纹可以是圆形指纹、MACCS指纹、原子对指纹和拓扑指纹中的任意一种。
5.根据权利要求1所述方法,其特征在于:第二步所述待测生物样本为人源、鼠源样本中一种或二种以上。
6.根据权利要求1所述方法,其特征在于:第二步所述一级质谱离子为化合物被质谱电离离子化后直接采集到的离子;二级质谱离子为一级质谱离子在施加一定能量碰撞碎裂后采集到的离子。
7.根据权利要求1所述方法,其特征在于:第三步所述的分子描述符是利用标样的结构信息在开源网站中计算所得;保留时间预测模型采用多重线性回归模型、弹性网络回归模型、随机森林模型、支持向量机模型和深度学习模型中任意一种方法建立。
8.根据权利要求1所述方法,其特征在于:第五步非靶向代谢组学数据色谱峰的鉴定方法为,以分子结构关联网络为参考,利用内源性代谢物的标准样品从非靶向超高效液相色谱-高分辨质谱代谢组实验数据鉴定50~100个代谢物作为种子代谢物,将其映射到第一步建立的分子结构关联网络中,从网络中获取与种子代谢物有直接连边的代谢物定义为相邻代谢物;将种子代谢物的二级质谱赋给相邻代谢代谢物,作为其拟二级质谱;设定搜索阈值,在实验数据中搜索与相邻代谢物m/z理论,tR预测,拟二级质谱匹配的代谢物峰,匹配成功则该代谢物峰鉴定完成;当有多个匹配结果时,对匹配结果进行打分,按得分由高到低排序,得分越高代谢物峰则其鉴定的准确度越高,但有多个搜索结果的代谢物峰不再作为新的种子参与下一轮定性;鉴定出的代谢物再作为新的种子,重复定性过程,直至没有新的代谢物被定性出来。
9.根据权利要求8所述方法,其特征在于:搜索阈值:
|tR预测-tR实测|/tR预测*100%<30%与|m/z理论-m/z实测|/m/z理论*1000000<10ppm,且实验二级质谱和相邻代谢物的拟二级质谱的相似度≥0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011407875.5A CN114609318B (zh) | 2020-12-03 | 2020-12-03 | 一种基于分子结构关联网络的规模化代谢组定性方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011407875.5A CN114609318B (zh) | 2020-12-03 | 2020-12-03 | 一种基于分子结构关联网络的规模化代谢组定性方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114609318A true CN114609318A (zh) | 2022-06-10 |
CN114609318B CN114609318B (zh) | 2024-03-12 |
Family
ID=81856899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011407875.5A Active CN114609318B (zh) | 2020-12-03 | 2020-12-03 | 一种基于分子结构关联网络的规模化代谢组定性方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114609318B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297355A (zh) * | 2013-07-17 | 2015-01-21 | 中国科学院大连化学物理研究所 | 一种基于液相色谱/质谱联用的拟靶标代谢组学分析方法 |
CN107729721A (zh) * | 2017-10-17 | 2018-02-23 | 中国科学院上海有机化学研究所 | 一种代谢物鉴定及紊乱通路分析方法 |
CN107817309A (zh) * | 2016-09-13 | 2018-03-20 | 中国科学院大连化学物理研究所 | 一种规模化代谢组定性方法 |
CN109920473A (zh) * | 2019-04-02 | 2019-06-21 | 中国科学院城市环境研究所 | 一种代谢组学标志物权重分析通用方法 |
-
2020
- 2020-12-03 CN CN202011407875.5A patent/CN114609318B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104297355A (zh) * | 2013-07-17 | 2015-01-21 | 中国科学院大连化学物理研究所 | 一种基于液相色谱/质谱联用的拟靶标代谢组学分析方法 |
CN107817309A (zh) * | 2016-09-13 | 2018-03-20 | 中国科学院大连化学物理研究所 | 一种规模化代谢组定性方法 |
CN107729721A (zh) * | 2017-10-17 | 2018-02-23 | 中国科学院上海有机化学研究所 | 一种代谢物鉴定及紊乱通路分析方法 |
CN109920473A (zh) * | 2019-04-02 | 2019-06-21 | 中国科学院城市环境研究所 | 一种代谢组学标志物权重分析通用方法 |
Non-Patent Citations (7)
Also Published As
Publication number | Publication date |
---|---|
CN114609318B (zh) | 2024-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Domingo-Almenara et al. | Annotation: a computational solution for streamlining metabolomics analysis | |
Kind et al. | Identification of small molecules using accurate mass MS/MS search | |
Draper et al. | Metabolite signal identification in accurate mass metabolomics data with MZedDB, an interactive m/z annotation tool utilising predicted ionisation behaviour'rules' | |
Hufsky et al. | Computational mass spectrometry for small-molecule fragmentation | |
Werner et al. | Mass spectrometry for the identification of the discriminating signals from metabolomics: current status and future trends | |
Vaniya et al. | Using fragmentation trees and mass spectral trees for identifying unknown compounds in metabolomics | |
Ridder et al. | Automatic compound annotation from mass spectrometry data using MAGMa | |
US8975577B2 (en) | System and method for grouping precursor and fragment ions using selected ion chromatograms | |
CN104237364B (zh) | 同位素模式识别 | |
Oberacher et al. | Evaluation of the sensitivity of the ‘Wiley registry of tandem mass spectral data, MSforID’with MS/MS data of the ‘NIST/NIH/EPA mass spectral library’ | |
WO2012104956A1 (ja) | 質量分析方法及び装置 | |
Zeng et al. | Alternative moving window factor analysis for comparison analysis between complex chromatographic data | |
JP6695556B2 (ja) | イメージング質量分析データ処理装置及び方法 | |
JP6149810B2 (ja) | 代謝物解析システム及び代謝物解析方法 | |
JP2011510300A (ja) | 複合混合物の化学成分の組成を決定するためのシステム、方法、及びコンピュータ読み取り可能な媒体 | |
JP2020532738A (ja) | 質量スペクトルデータの取得および解析方法 | |
JP5810983B2 (ja) | 質量分析を用いた化合物同定方法及び化合物同定システム | |
Godzien et al. | Metabolite annotation and identification | |
Yates III | Mass spectrometry as an emerging tool for systems biology | |
CN112634997A (zh) | 一种甾醇数据库建立及甾醇分析方法 | |
CN114609318B (zh) | 一种基于分子结构关联网络的规模化代谢组定性方法 | |
CN114594171B (zh) | 一种代谢组深度注释方法 | |
Samaraweera et al. | MolFind 2: A Protocol for Acquiring and Integrating MS 3 Data to Improve In Silico Chemical Structure Elucidation for Metabolomics | |
James | XLIM-MS Towards the Development of a Novel approach to Cross-linking Mass Spectrometry | |
Nash | Improving metabolite annotation and identification in untargeted UHPLC-MS metabolomics studies |
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 |