CN113791018B - 基于电阻抗信号的单细胞机械本征参数测量系统和方法 - Google Patents
基于电阻抗信号的单细胞机械本征参数测量系统和方法 Download PDFInfo
- Publication number
- CN113791018B CN113791018B CN202111063359.XA CN202111063359A CN113791018B CN 113791018 B CN113791018 B CN 113791018B CN 202111063359 A CN202111063359 A CN 202111063359A CN 113791018 B CN113791018 B CN 113791018B
- Authority
- CN
- China
- Prior art keywords
- cell
- impedance
- flow channel
- cells
- mechanical
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 51
- 238000002847 impedance measurement Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000008878 coupling Effects 0.000 claims abstract description 11
- 238000010168 coupling process Methods 0.000 claims abstract description 11
- 238000005859 coupling reaction Methods 0.000 claims abstract description 11
- 238000004070 electrodeposition Methods 0.000 claims abstract description 4
- 230000000877 morphologic effect Effects 0.000 claims abstract description 3
- 230000008602 contraction Effects 0.000 claims description 38
- 238000001125 extrusion Methods 0.000 claims description 25
- 230000008569 process Effects 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 14
- 238000013461 design Methods 0.000 claims description 9
- 239000000243 solution Substances 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 8
- 230000005284 excitation Effects 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 8
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000000875 corresponding effect Effects 0.000 claims description 6
- 230000005489 elastic deformation Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000003990 capacitor Substances 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000003071 parasitic effect Effects 0.000 claims description 3
- 229920001992 poloxamer 407 Polymers 0.000 claims description 3
- 239000004094 surface-active agent Substances 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 7
- 210000004027 cell Anatomy 0.000 description 169
- 238000010586 diagram Methods 0.000 description 17
- 238000012512 characterization method Methods 0.000 description 7
- 108010076504 Protein Sorting Signals Proteins 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000001413 cellular effect Effects 0.000 description 3
- 239000004205 dimethyl polysiloxane Substances 0.000 description 3
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 102000029749 Microtubule Human genes 0.000 description 2
- 108091022875 Microtubule Proteins 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 239000006285 cell suspension Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 239000007850 fluorescent dye Substances 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 210000004688 microtubule Anatomy 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- -1 Polydimethylsiloxane Polymers 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 239000006143 cell culture medium Substances 0.000 description 1
- 210000000170 cell membrane Anatomy 0.000 description 1
- 230000008614 cellular interaction Effects 0.000 description 1
- 238000012258 culturing Methods 0.000 description 1
- 230000003436 cytoskeletal effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 238000001215 fluorescent labelling Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 210000003963 intermediate filament Anatomy 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 201000004792 malaria Diseases 0.000 description 1
- 239000002609 medium Substances 0.000 description 1
- 210000003632 microfilament Anatomy 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000009832 plasma treatment Methods 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 208000007056 sickle cell anemia Diseases 0.000 description 1
- 238000002174 soft lithography Methods 0.000 description 1
- 238000005476 soldering Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/1031—Investigating individual particles by measuring electrical or magnetic effects
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/50273—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip characterised by the means or forces applied to move the fluids
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L3/00—Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
- B01L3/50—Containers for the purpose of retaining a material to be analysed, e.g. test tubes
- B01L3/502—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
- B01L3/5027—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip
- B01L3/502761—Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures by integrated microfluidic structures, i.e. dimensions of channels and chambers are such that surface tension forces are important, e.g. lab-on-a-chip specially adapted for handling suspended solids or molecules independently from the bulk fluid flow, e.g. for trapping or sorting beads, for physically stretching molecules
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L2200/00—Solutions for specific problems relating to chemical or physical laboratory apparatus
- B01L2200/10—Integrating sample preparation and analysis in single entity, e.g. lab-on-a-chip concept
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01L—CHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
- B01L2300/00—Additional constructional details
- B01L2300/08—Geometry, shape and general structure
- B01L2300/0887—Laminated structure
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Clinical Laboratory Science (AREA)
- Hematology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Electrochemistry (AREA)
- Fluid Mechanics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明公开了属于细胞分析、测量技术领域，特别涉及一种基于电阻抗信号的单细胞机械本征参数测量系统和方法。该基于电阻抗信号的单细胞机械本征参数测量系统由微流控系统、多频阻抗测量系统和实时处理算法系统串联组成：其单细胞机械本征参数测量系统的测量方法包括阻抗测量系统测量、采集单细胞的阻抗值是利用流道结构和电极位置的空间耦合，以及电阻抗测量电信号的高时间分辨率特性，得到单细胞的时空形态信息，通过拟合单细胞力学模型实现对单细胞的高通量实时本征机械参数测量。本发明无需传统的高速相机，即能求解细胞动态空间位置与形变，实现对单细胞机械本征参数的高通量实时测量。
Description
技术领域
本发明属于细胞分析、测量技术领域，特别涉及一种基于电阻抗信号的单细胞机械本征参数测量系统和方法。
背景技术
单细胞表征提供了细胞基本结构、功能信息和病理状态，在揭示细胞异质性中起着重要作用，对生命科学研究、疾病诊断和个性化医学意义重大。一直以来，荧光标记是单细胞分析和成像的主要工具，它通过对细胞标记具有分子特异性的检测标签(如荧光染料、量子点、磁珠和稳定同位素等)来识别细胞中的组分及其状态。对细胞进行荧光标记不但需要细胞特异性的先验知识，而且标记过程的侵入式操作会改变细胞状态，使分析过程复杂化，限制后续分析。相比之下，细胞的生物物理特性(如电学特性和机械特性)同样与细胞内的分子组成有关，已被证明是诊断疾病(癌症、疟疾、糖尿病和镰状细胞性贫血等)的有效生物标记，且生物物理特性表征无需标记细胞，表征过程中细胞状态基本不变，表征后细胞仍然可以进行下一步操作和分析，如分选、培养、组学分析等，因此广受关注。
单细胞生物物理特性表征已有不少经典方法，如：利用原子力显微镜、微管吸吮、光纤拉伸、流体挤压、力诱导变形等测量单细胞机械参数。这些方法大都在准确性上表现优异，在帮助人类对细胞的理解和认识上取得了巨大的进步，至今也在不断演进，但在细胞样品处理上存在通量低和系统不易操作等固有问题，实用性上大打折扣。近年来，微流控技术已发展成为单细胞生物物理特性表征的有力工具。微流控技术所需样品体积小、整合能力强、生物兼容性好、响应速度快，在通量和系统集成度方面具有很强的优势，促进了细胞特性表征的发展。
作为生物物理特性之一，细胞机械特性反映了细胞膜、细胞骨架网络(肌动蛋白丝、中间丝和微管)及细胞核间的联系与特征。机械特性表征的基本原理是在流道中产生细胞的受力变形，通过细胞变形过程和力学模型求取机械特性参数，往往要使用高速相机捕捉毫秒级的受力变形动态过程。美国加州大学伯克利分校Rosenbluth课题组提出具有多平行收缩微流道结构的微流控器件，利用结构诱导变形原理，以细胞穿越收缩流道的时间来表征细胞的机械特性。德国德累斯顿工业大学Guck课题组提出基于流体诱导变形的微流控器件，利用鞘液的作用力使细胞发生变形，避免了细胞与流道壁粘附作用对变形性的影响，采用变形量表征细胞的机械特性并实现分选。为保证高通量，这些方法需要使用复杂昂贵的高速成像系统得到细胞的动态受力变形过程，且需要对应的图像处理算法与高算力实现对细胞变形量的提取，因此只能实时得到对算力要求很低的变形量或通过时间等现象学参数来定性地表征单细胞，无法解算出细胞的机械特性本征参数(杨氏模量、流度)，导致测量结果高度依赖于测量平台，并且无法在不同平台之间进行比较。
在高通量测量场景下，只能以现象学参数(如通过时间、变形量和松弛指数)来定性地表征单个细胞，使得测量结果高度依赖于测量平台，并且无法在不同平台之间进行比较。对单细胞机械本征参数的测量则往往需要使用昂贵的高速成像装置和复杂的图像处理算法来捕捉毫秒级的细胞受力变形动态过程，产生的大量图像数据，耗费大量算力，只能离线分析，难以实现实时细胞检测。为了摒弃单细胞机械本征参数测量传统方法中必须使用的复杂庞大的高速成像装置，从而降低系统和处理算法的复杂性，为了提高高通量场景下对单细胞机械特性表征的能力，本发明利用阻抗测量的高时间分辨率特性，采用流道结构与电极形状耦合的框架，提出基于阻抗脉冲提取脉冲宽度和形状的方法，无需传统的高速相机，即能求解细胞动态空间位置与形变，实现对单细胞机械本征参数的高通量实时测量。
发明内容
本发明的目的是提出一种基于电阻抗信号的单细胞机械本征参数测量系统和方法。其特征在于，所述基于电阻抗信号的单细胞机械本征参数测量系统，包括微流控系统、多频阻抗测量系统和实时处理算法系统串联组成：
所述微流控系统主要由微流控器件和进样驱动装置组成，负责完成驱动大量单细胞聚束列队高速通过测量区域；
所述阻抗测量系统主要由直接数字频率合成(DDS)信号发生模块、锁定放大器模块和信号采集模块组成，负责测量并采集单细胞的阻抗值；
所述实时处理算法系统主要由数字滤波、事件检测与有效信号提取算法组成，负责实现对阻抗信号进行滤波、单细胞事件检测与细胞有效形变区域信号提取以及实时获取待测单细胞杨氏模量等机械本征参数；
所述微流控芯片结构采用收缩流道的尺寸设计，测量电极部署在收缩流道的入口与出口处，以使单细胞挤压通过；当细胞挤压通过收缩流道时，测量的阻抗幅度将与细胞的挤压伸长长度Lp成正比；因此，细胞进入收缩流道的形变过程能够通过阻抗信号来描述。
所述微流控芯片结构，考虑到制造和光学观察的便利性，微流控器件芯片设计为两层。
所述阻抗测量系统测量并采集单细胞的阻抗值是利用流道结构和电极位置的空间耦合，以及电阻抗测量电信号的高时间分辨率特性，得到单细胞的时空形态信息，通过拟合单细胞力学模型实现对单细胞的高通量实时本征机械参数测量。
所述收缩流道的尺寸设计为10μm宽和20μm长以使单细胞挤压通过，测量电极设计为30μm长，20μm宽和20μm间距。
一种基于电阻抗信号的单细胞机械本征参数测量系统的测量方法，其特征在于，该方法是一种结合微流控技术和阻抗检测技术的高通量实时机械参数测量方法；利用阻抗流式芯片设计中测量电极与收缩流道结构的空间耦合，使得测量的阻抗信号中包含空间信息，从而得到细胞通过收缩流道的挤压变形动态信息，并利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量，降低测量系统的复杂程度；具体包括如下步骤：
1)流道修饰，为了减少细胞与流道表面、侧壁粘附，在进样之前，将微流道用1wt％的Pluronic F-127表面活性剂预处理15分钟；
2)细胞进给，使用注射泵驱动细胞液，将细胞液样品以1μL/min的流速通入微流道中；
3)细胞收缩挤压，单细胞在流道中列队挤压通过收缩流道测量区域，
4)其弹性形变由阻抗测量系统采集；
5)实时求解本征参数，电阻抗信号实时解算，进而由机械模型得到单细胞机械本征参数。
所述利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量，该单细胞机械本征参数测量广泛使用的单细胞幂律流变模型，单细胞的挤压伸长长度Lp定义为细胞前沿与收缩入口的距离，当细胞挤压通过收缩流道时存在三个阶段：I阶段，细胞刚进入收缩区，Lp快速非线性增加；Ⅱ阶段，细胞开始收缩进入收缩流道并发生弹性变形，导致细胞进入收缩流道的速度降低并呈幂律变化，从而该阶段的Lp变化将能够反映细胞的机械特性；Ⅲ阶段，细胞完全变形并进入收缩流道并快速向前移动，该阶段的移动速度主要与细胞和流道的摩擦力有关。
所述由机械模型得到单细胞机械本征参数，在流道中产生细胞的受力变形，通过细胞变形过程和力学模型求取机械特性参数；为获得单细胞的机械本征参数，利用电阻抗测量的高时间分辨率特性，通过流道结构与电极形状的耦合实现对单细胞的形状与机械特性的测量，其关键在于电信号能够关联机械形变；对应于流道壁接触式挤压变形结构，将测量电极部署在收缩流道的入口，当细胞挤压变形进入收缩流道时将等体积置换测量区域内的溶液；此时细胞-培养基混合体系的阻抗可由下式计算
其中CDL为双电层电容，
为挤压细胞的阻抗，
为溶液阻抗，ω为激励信号的角频率。由于寄生电容Cs的存在，测量系统的整体阻抗为：
由于Cs的值一般在0.1-0.5pF之间，在激励信号频率为1MHz时，ω*Cs≈0，因此
而由电阻计算公式可知：
其中ρ为电阻率，L为长度，s为流道横截面，又由于Lcell+Lmed＝Lchannel，可得：
由于Lchannel为设计常量，因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关，也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联。
所述计算求解细胞的机械本征参数，需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化，在Ⅱ阶段中，细胞的剪切蠕变J(t)由下式
其中Rp是流道的有效直径，可由
计算，Φp是一个常数，在矩形截面的微流道中通常取值为2.1，Lp(t)是细胞在t时刻的挤压伸长量，ΔP为施加压强；
根据幂律流变模型，瞬时剪切蠕变J(t)也可表示为：
其中AJ是反映细胞变形能力的剪切柔量常数，α是幂律指数也称为流度，t0一般取t0＝1s；
细胞的杨氏模量AG可由下式计算：
其中Γ为伽马函数，因此由公式(5)-(7)可得
因此，细胞的机械本征参数AG和α可以通过最小化细胞的测量挤压伸长长度Lcell(t)及其模型估计值Lp(t)之间的方差来获得
其中i表示连续的时间点，Lcell(t)为测量值，由公式(4)获得；Lp(t)为模型估计值，基于网格搜索与最小二乘拟合方法，从参数网格中搜索AG和α代入公式(8)中获得。
本发明的有益效果是本发明在单细胞阻抗流式分析技术的基础上，利用流道结构和电极位置的空间耦合，以及电阻抗测量电信号的高时间分辨率特性，得到单细胞的时空形态信息，通过拟合单细胞力学模型实现对单细胞的高通量实时本征机械参数测量。
附图说明
图1为单细胞机械本征参数实时测量系统示意框图。
图2为单细胞测量系统示意图；其中，(a)微流控芯片结构示意图；(b)微流控器件加工示意图；(c)一种微流控器件实例图；(d)测量区域显微结构图。
图3为单细胞阻抗测量示意图；其中，(a)单细胞阻抗测量等效电路模型；(b)细胞进入过程中的阻抗信号序列。
图4为测量电路示意图。
图5为单细胞性能测试，其中(a)细胞弹性形变过程中的有效阻抗及其导数信号图；(b)细胞幂律流变模型拟合结果图。
图6为细胞事件的阻抗特性，其中，(a)连续细胞事件的阻抗信号序列；(b)挤压细胞的阻抗与对应的伸长长度；(c)细胞通过收缩流道的连续显微图像，其中椭圆框标示了细胞的位置与挤压时的前端。
图7为单细胞机械本征参数测量流程框图。
具体实施方式
本发明的目的是提出一种基于电阻抗信号的单细胞机械本征参数测量系统和方法。下面结合附图对本发明予以说明。
所述基于电阻抗信号的单细胞机械本征参数测量系统，
如图1所示的单细胞机械本征参数实时测量系统示意框图，其包括三个子系统：微流控系统、多频阻抗测量系统和实时处理算法系统串联组成：其中微流控系统主要由微流控器件和进样驱动装置组成，负责完成驱动大量单细胞聚束列队高速通过测量区域；阻抗测量系统主要由直接数字频率合成(DDS)信号发生模块、锁定放大器模块和信号采集模块组成，负责测量并采集单细胞的阻抗值；实时处理算法系统主要由数字滤波、事件检测与有效信号提取算法组成，负责实现对阻抗信号进行滤波、单细胞事件检测与细胞有效形变区域信号提取以及实时获取待测单细胞杨氏模量等机械本征参数。
如图2所示为单细胞测量系统示意图；其图2中，(a)微流控芯片结构示意图；(b)微流控器件加工示意图；(c)一种微流控器件实例图；(d)测量区域显微结构图。
在(a)所示微流控芯片结构示意图中，以直径为15-20μm的单细胞为实例，检测区域流道入口的尺寸设计为20μm以使单细胞聚束列队依次通过检测区域，检测区域收缩流道的尺寸设计为宽10μm、长20μm，以使单细胞挤压通过，测量电极设计为长30μm、宽20μm、间距20μm；部署在收缩流道的入口与出口处。当细胞挤压通过收缩流道时，测量的阻抗幅度将与细胞的挤压伸长长度Lp成正比，因此，细胞进入收缩流道的形变过程可以通过阻抗信号来描述。考虑到制造和光学观察的便利性，微流控器件芯片设计为两层，其中一个为，底层以玻璃作为基底并配置有一对平面电极以进行阻抗检测，顶层则配置有聚二甲基硅氧烷(PDMS)微流道。对应于这个设计，微流控器件的制造过程如图2(b)所示，微流道层通过软光刻技术制作，利用PDMS倒模形成，电极层通过lift-off技术图案化到玻璃基底上。在用氧等离子体处理并烘烤之后，将电极层和微流道层对准牢固地键合在一起。最终，为了方便施加激励信号和采集测量信号，电极与其图案匹配的定制印刷电路板通过焊锡焊接在一起。图2(c)与图2(d)分别为微流控器件的实例图与器件的显微结构图，其中虚线方框内为测量区域，细胞悬浮液进口出口通过塑料软管与微流泵相连，通过微流泵将细胞悬浮液推注入微流道中。
如图3所示为单细胞阻抗测量示意图；其中，(a)单细胞阻抗测量等效电路模型；(b)细胞进入过程中的阻抗信号序列。
单细胞阻抗测量中机械特性表征的基本原理是在流道中产生细胞的受力变形，通过细胞变形过程和力学模型求取机械特性参数。为获得单细胞的机械本征参数，利用电阻抗测量的高时间分辨率特性，通过流道结构与电极形状的耦合实现对单细胞的形状与机械特性的测量，其关键在于电信号能够关联机械形变。对应于流道壁接触式挤压变形结构，本发明将测量电极部署在收缩流道的入口，当细胞挤压变形进入收缩流道时将等体积置换测量区域内的溶液，此时电学等效模型如图3中(a)所示，细胞-培养基混合体系的阻抗可由下式计算
其中CDL为双电层电容，
为挤压细胞的阻抗，
为溶液阻抗，ω为激励信号的角频率。由于寄生电容Cs的存在，测量系统的整体阻抗为：
由于Cs的值一般在0.1-0.5pF之间，在激励信号频率为1MHz时，ω*Cs≈0，因此
而由电阻计算公式可知：
其中ρ为电阻率，L为长度，s为流道横截面，又由于Lcell+Lmed＝Lchannel，可得：
由于Lchannel为设计常量，因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关，也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联。
本发明中的机械本征参数测量原理是基于广泛使用的单细胞幂律流变模型，如图2中(a)与图3中(b)所示，细胞的挤压伸长长度Lp定义为细胞前沿与收缩入口的距离，当细胞挤压通过收缩流道时存在三个阶段：I阶段，细胞刚进入收缩区，Lp快速非线性增加；Ⅱ阶段，细胞开始收缩进入收缩流道并发生弹性变形，导致细胞进入收缩流道的速度降低并呈幂律变化，从而该阶段的Lp变化将能够反映细胞的机械特性；Ⅲ阶段，细胞完全变形并进入收缩流道并快速向前移动，该阶段的移动速度主要与细胞和流道的摩擦力有关。因此，为了计算求解细胞的机械本征参数，需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化。在Ⅱ阶段中，细胞的剪切蠕变J(t)由下式给出：
其中Rp是流道的有效直径，可由
计算，Φp是一个常数，在矩形截面的微流道中通常取值为2.1，Lp(t)是细胞在t时刻的挤压伸长量，ΔP为施加压强。
根据幂律流变模型，瞬时剪切蠕变J(t)也可表示为：
其中AJ是反映细胞变形能力的剪切柔量常数，α是幂律指数也称为流度，t0一般取t0＝1s。
细胞的杨氏模量AG可由式(7)表示，
其中Γ为伽马函数，因此由式(5)-(7)可得
因此，细胞的机械本征参数AG和α可以通过最小化细胞的测量挤压伸长长度Lcell(t)及其模型估计值Lp(t)之间的方差来获得
其中i表示连续的时间点，Lcell(t)为测量值，由公式(4)获得；Lp(t)为模型估计值，基于网格搜索与最小二乘拟合方法，从参数网格中搜索AG和α代入公式(8)中获得。
如图4所示的阻抗测量的基本原理，将电流信号通过跨阻放大器(TIA)转换为电压信号，并通过锁相放大器(LIA)检测所得的微弱响应信号。一个简单的实例为，信号发生器基于DDS芯片AD9958(ADI，美国)搭建，TIA基于高增益带宽、低偏置电流运算放大器芯片OPA657(TI，美国)搭建，解调器基于宽带四象限电压输出乘法器AD835芯片(ADI，美国)搭建和具有5KHz截止频率的低通滤波器和放大器由高精度，低噪声运算放大器OPA227(TI，美国)搭建。最终输出的直流信号由计算机通过数据采集卡(NI，PCI-6289)以100KHz的采样率进行采样获取，并在计算机上由软件系统进行后续处理。
如图5所示，其中(a)细胞弹性形变过程中的有效阻抗及其导数信号图；(b)细胞幂律流变模型拟合结果图；是对图4所示的阻抗测量系统采集的阻抗信号的实时处理，首先基于事件检测原理提取细胞通过收缩流道的有效阻抗信号，随后计算阻抗信号一阶导数，当细胞挤压通过收缩流道分为三个阶段：第I阶段，细胞进入收缩流道区域，细胞的前沿以较快的速度进入收缩流道测量区，其范围从导数的第一个零点到峰值点；在第Ⅱ阶段，由于细胞剪切蠕变特性，细胞挤压通过的速度降低，细胞发生弹性变形挤压通过收缩流道，并呈幂律变化，从而该阶段的Lp变化将能够反映细胞的机械特性；从峰值点到第二个零点；在第Ⅲ阶段，细胞完全挤入收缩区并开始挤出，该范围从第二个零点到第三个零点。因此，根据单细胞的机械模型，该算法利用寻峰原理定位峰值点及峰值点后的第一个零点，并从有效阻抗信号中提取Ⅱ阶段的弹性变形相关信号用于拟合幂律流变模型，拟合结果如图5中(b)所示的细胞幂律流变模型拟合结果图。该算法经验证，在普通PC机上能够在100ms级时间里完成单细胞的机械本征参数杨氏模量AG和流度α的计算，从而实现>100cells/min的高通量单细胞检测。
为了验证该测量系统和方法，将K562细胞系通入上述系统，图6中(a)显示了包含252个细胞事件的阻抗信号序列，图6中(b)放大显示了一个细胞经历挤压全过程的阻抗信号。为了确认阻抗信号对应于机械变形；图6中(c)所示基于高速相机同步记录了此细胞的连续显微图像(显示了具有代表性的8帧)，其中图像帧1中的圆圈标示出了收缩流道之前的整个细胞，图像帧2-7中的椭圆标示了细胞挤压通过收缩流道的细胞前端，图像帧8中的椭圆标示出了细胞最后离开收缩流道的细胞尾端，可以看到，连续图像帧1-7提取的细胞挤压伸长量Lp与阻抗匹配得很好，证明了建立的模型能够很好地描述细胞挤压通过收缩流道的动态变化，验证了基于阻抗信号的测量策略能够适用于单细胞机械本征参数求解。
如图7所示的单细胞机械本征参数测量系统；搭建的锁定放大器系统的输入与输出分别与微流控器件的测量电极相连，用于施加激励信号和读出响应信号。及描述了阻抗检测技术的高通量实时机械参数测量方法；利用阻抗流式芯片设计中测量电极与收缩流道结构的空间耦合，使得测量的阻抗信号中包含空间信息，从而得到细胞通过收缩流道的挤压变形动态信息，并利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量，降低测量系统的复杂程度；为了减少细胞与流道表面的相互作用如侧壁粘附，在进样之前，将微流道用1wt％(1x PBS)的Pluronic F-127表面活性剂预处理15分钟。随后使用合适的细胞液驱动方式如注射泵将样品以合适的流速(如1μL/min)通入微流道中，单细胞在流道中列队挤压通过收缩流道测量区域，其弹性形变由阻抗测量系统采集电阻抗信号实时解算，进而由机械模型得到单细胞机械本征参数。
Claims (6)
1.一种基于电阻抗信号的单细胞机械本征参数测量系统，其特征在于，所述基于电阻抗信号的单细胞机械本征参数测量系统由微流控系统、多频阻抗测量系统和实时处理算法系统串联组成：
所述微流控系统主要由微流控器件和进样驱动装置组成，负责完成驱动大量单细胞聚束列队高速通过测量区域；
所述阻抗测量系统主要由直接数字频率合成(DDS)信号发生模块、锁定放大器模块和信号采集模块组成，负责测量并采集单细胞的阻抗值；
所述实时处理算法系统主要由数字滤波、事件检测与有效信号提取算法组成，负责实现对阻抗信号进行滤波、单细胞事件检测与细胞有效形变区域信号提取以及实时获取待测单细胞杨氏模量机械本征参数，细胞有效形变区域信号为细胞挤压通过收缩流道的ⅠⅠ阶段区域时的阻抗信号；
微流控芯片结构采用收缩流道的尺寸设计，测量电极部署在收缩流道的入口与出口处，以使单细胞挤压通过；当细胞挤压通过收缩流道时，测量的阻抗幅度将与细胞的挤压伸长长度成正比；因此，细胞进入收缩流道的形变过程能够通过阻抗信号来描述。
2.根据权利要求1所述基于电阻抗信号的单细胞机械本征参数测量系统，其特征在于，所述微流控芯片结构，考虑到制造和光学观察的便利性，微流控器件芯片设计为两层。
3.根据权利要求1所述基于电阻抗信号的单细胞机械本征参数测量系统，其特征在于，所述阻抗测量系统测量并采集单细胞的阻抗值是利用流道结构和电极位置的空间耦合，以及电阻抗测量电信号的高时间分辨率特性，得到单细胞的时空形态信息，通过拟合单细胞力学模型实现对单细胞的高通量实时本征机械参数测量。
4.根据权利要求1所述基于电阻抗信号的单细胞机械本征参数测量系统，其特征在于，所述收缩流道的尺寸设计为10μm宽和20μm长以使单细胞挤压通过，测量电极设计为30μm长，20μm宽和20μm间距。
5.一种基于电阻抗信号的单细胞机械本征参数测量系统的测量方法，其特征在于，该方法是一种结合微流控技术和阻抗检测技术的高通量实时机械参数测量方法；利用阻抗流式芯片设计中测量电极与收缩流道结构的空间耦合，使得测量的阻抗信号中包含空间信息，从而得到细胞通过收缩流道的挤压变形动态信息，并利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量，降低测量系统的复杂程度；具体包括如下步骤：
1)流道修饰，为了减少细胞与流道表面、侧壁粘附，在进样之前，将微流道用1wt％的Pluronic F-127表面活性剂预处理15分钟；
2)细胞进给，使用注射泵驱动细胞液，将细胞液样品以1μL/min的流速通入微流道中；
3)细胞收缩挤压，单细胞在流道中列队挤压通过收缩流道测量区域；
4)其弹性形变由阻抗测量系统采集；为获得单细胞的机械本征参数，利用电阻抗测量的高时间分辨率特性，通过流道结构与电极形状的耦合实现对单细胞的形状与机械特性的测量，其关键在于电信号能够关联机械形变；对应于流道壁接触式挤压变形结构，将测量电极部署在收缩流道的入口，当细胞挤压变形进入收缩流道时将等体积置换测量区域内的溶液；此时细胞-培养基混合体系的阻抗可由下式计算，
其中CDL为双电层电容，
为挤压细胞的阻抗，
为溶液阻抗，ω为激励信号的角频率；由于寄生电容Cs的存在，测量系统的整体阻抗为：
由于Cs的值一般在0.1-0.5pF之间，在激励信号频率为1MHz时，ω*Cs≈0，因此
而由电阻计算公式可知：
其中ρcell为细胞电阻率，Lcell为细胞伸长长度，ρmed为溶液电阻率，Lmed为溶液长度，s为流道横截面，又由于Lcell+Lmed＝Lchannel，可得：
由于流道长度Lchannel为设计常量，因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关，也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联；
5)实时求解本征参数，电阻抗信号实时解算，进而由机械模型得到单细胞机械本征参数；计算求解细胞的机械本征参数，需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化，在Ⅱ阶段中，细胞的剪切蠕变J(t)由下式给出：
其中Rp是流道的有效直径，可由
计算，Φp是一个常数，在矩形截面的微流道中通常取值为2.1，Lp(t)是细胞在t时刻的挤压伸长量，ΔP为施加压强；
根据幂律流变模型，瞬时剪切蠕变J(t)也可表示为：
其中AJ是反映细胞变形能力的剪切柔量常数，α是幂律指数也称为流度，t0一般取t0＝1s；
细胞的杨氏模量AG可由下式计算：
其中Γ为伽马函数，因此由公式(5)-(7)可得
因此，细胞的机械本征参数AG和α可以通过最小化细胞的Lcell(t)及其Lp(t)之间的方差来获得
其中i表示连续的时间点，Lcell(t)由测量得到，由公式(4)获得；Lp(t)由模型估计得到，基于网格搜索与最小二乘拟合方法，从参数网格中搜索AG和α代入公式(8)中获得。
6.根据权利要求5所述基于电阻抗信号的单细胞机械本征参数测量系统的测量方法，其特征在于，所述利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量，该单细胞机械本征参数测量使用单细胞幂律流变模型，单细胞的挤压伸长长度定义为细胞前沿与收缩入口的距离，当细胞挤压通过收缩流道时存在三个阶段：Ⅰ阶段，细胞刚进入收缩区，单细胞的挤压伸长长度快速非线性增加；Ⅱ阶段，细胞开始收缩进入收缩流道并发生弹性变形，导致细胞进入收缩流道的速度降低并呈幂律变化，从而该阶段的单细胞的挤压伸长长度变化将能够反映细胞的机械特性；Ⅲ阶段，细胞完全变形并进入收缩流道并快速向前移动，该阶段的移动速度主要与细胞和流道的摩擦力有关。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111063359.XA CN113791018B (zh) | 2021-09-10 | 2021-09-10 | 基于电阻抗信号的单细胞机械本征参数测量系统和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111063359.XA CN113791018B (zh) | 2021-09-10 | 2021-09-10 | 基于电阻抗信号的单细胞机械本征参数测量系统和方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113791018A CN113791018A (zh) | 2021-12-14 |
| CN113791018B true CN113791018B (zh) | 2022-07-15 |
Family
ID=79183193
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111063359.XA Active CN113791018B (zh) | 2021-09-10 | 2021-09-10 | 基于电阻抗信号的单细胞机械本征参数测量系统和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113791018B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115417585B (zh) * | 2022-09-22 | 2023-08-29 | 重庆理工大学 | 一种多收缩孔径玻璃毛细管的加工装置及加工方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439241A (zh) * | 2013-08-23 | 2013-12-11 | 东南大学 | 单细胞多参数表征的微流控芯片检测系统 |
| CN103674813A (zh) * | 2013-09-22 | 2014-03-26 | 中国科学院电子学研究所 | 基于微流控技术测量单个细胞杨氏模量的方法 |
| CN104251810A (zh) * | 2013-12-18 | 2014-12-31 | 中国科学院电子学研究所 | 同时表征单细胞杨氏模量和细胞膜比电容的系统 |
| CN106124388A (zh) * | 2016-06-12 | 2016-11-16 | 中国科学院电子学研究所 | 毛细管进样系统及进样方法、单细胞电学特性检测系统 |
| CN107462512A (zh) * | 2017-08-18 | 2017-12-12 | 中国科学院电子学研究所 | 单细胞固有电学特性检测装置及方法 |
| CN108303364A (zh) * | 2018-02-02 | 2018-07-20 | 中国科学院电子学研究所 | 高通量检测细胞尺寸、尺寸无关单细胞电参数装置和方法 |
| CN110777077A (zh) * | 2019-10-22 | 2020-02-11 | 哈尔滨工业大学（深圳） | 检测细胞弹性模量的装置及方法 |
| WO2020180252A1 (en) * | 2019-03-05 | 2020-09-10 | Singapore University Of Technology And Design | Microfluidic device for single cell processing and method and system for single cell biophysical phenotyping using the microfludic device |
-
2021
- 2021-09-10 CN CN202111063359.XA patent/CN113791018B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439241A (zh) * | 2013-08-23 | 2013-12-11 | 东南大学 | 单细胞多参数表征的微流控芯片检测系统 |
| CN103674813A (zh) * | 2013-09-22 | 2014-03-26 | 中国科学院电子学研究所 | 基于微流控技术测量单个细胞杨氏模量的方法 |
| CN104251810A (zh) * | 2013-12-18 | 2014-12-31 | 中国科学院电子学研究所 | 同时表征单细胞杨氏模量和细胞膜比电容的系统 |
| CN106124388A (zh) * | 2016-06-12 | 2016-11-16 | 中国科学院电子学研究所 | 毛细管进样系统及进样方法、单细胞电学特性检测系统 |
| CN107462512A (zh) * | 2017-08-18 | 2017-12-12 | 中国科学院电子学研究所 | 单细胞固有电学特性检测装置及方法 |
| CN108303364A (zh) * | 2018-02-02 | 2018-07-20 | 中国科学院电子学研究所 | 高通量检测细胞尺寸、尺寸无关单细胞电参数装置和方法 |
| WO2020180252A1 (en) * | 2019-03-05 | 2020-09-10 | Singapore University Of Technology And Design | Microfluidic device for single cell processing and method and system for single cell biophysical phenotyping using the microfludic device |
| CN110777077A (zh) * | 2019-10-22 | 2020-02-11 | 哈尔滨工业大学（深圳） | 检测细胞弹性模量的装置及方法 |
Non-Patent Citations (6)
| Title |
|---|
| Characterizing Deformability and Electrical Impedance of Cancer Cells in a Microfluidic Device;Ying Zhou等;《analytical chemistry》;20171127;第90卷;第912-919页 * |
| Classification of cell types using a microfluidic device for mechanical and electrical measurement on single cells;Jian Chen等;《Lab on a Chip》;20111231;第11卷;第3174–3181页 * |
| High-throughput and label-free parasitemia quantification and stage differentiation for malaria-infected red blood cells;Xiaonan Yang等;《Biosensors and Bioelectronics》;20170708;第98卷;第408-411页 * |
| Microfluidics-Based Assessment of Cell Deformability;Andrea Adamo等;《analytical chemistry》;20120620;第84卷;第6438-6443页 * |
| 基于微流控技术的单细胞生物物理特性表征;唐文来等;《化学进展》;20140525;第26卷(第6期);第1050-1064页 * |
| 微流控阻抗流式细胞仪在单细胞检测中的应用;冯迪等;《化学进展》;20201130;第33卷(第4期);第555-567页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113791018A (zh) | 2021-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104251810B (zh) | 同时表征单细胞杨氏模量和细胞膜比电容的系统 | |
| CN101726585A (zh) | 一种基于微流控芯片的流式细胞仪 | |
| CN201548547U (zh) | 一种基于微流控芯片的流式细胞分析装置 | |
| CN113567326A (zh) | 一种高通量实时单细胞电学本征参数测量系统及方法 | |
| CN113791018B (zh) | 基于电阻抗信号的单细胞机械本征参数测量系统和方法 | |
| CN103439241A (zh) | 单细胞多参数表征的微流控芯片检测系统 | |
| CN103923825A (zh) | 一种集成细胞分选及检测的微流控芯片系统 | |
| Kang et al. | Differentiation between normal and cancerous cells at the single cell level using 3-D electrode electrical impedance spectroscopy | |
| CN109107621A (zh) | 基于细胞形变量和介电泳力的癌细胞分离装置及控制系统 | |
| CN113008766B (zh) | 一种基于微流控多次挤压的癌细胞动态行为检测系统 | |
| CN111359687B (zh) | 基于电学及荧光信号的微流控芯片和血细胞分析装置 | |
| CN109967143A (zh) | 一种基于微流控显微系统的细胞大小检测方法 | |
| CN104713932A (zh) | 一种交流模式的多参数纳米孔单分子分析仪 | |
| CN113533171A (zh) | 基于深度学习和微流控芯片的细胞形变检测方法及系统 | |
| CN112268934A (zh) | 一种针对循环肿瘤细胞的检测芯片及其检测方法 | |
| Li et al. | Intelligent droplet tracking with correlation filters for digital microfluidics | |
| Zhang et al. | Tunable microfluidic chip for single-cell deformation study | |
| EP4347128A2 (en) | Method and system for single-cell biophysical profiling using a microfluidic device | |
| CN101231262B (zh) | 用于单细胞组学的方法 | |
| Gajasinghe et al. | Miniaturized system for tumor cell detection and differentiation | |
| CN113413148A (zh) | 一种提高阻抗生物传感器效果的信号处理方法及装置 | |
| CN105181852A (zh) | 一种鉴别延胡索生熟药材的方法 | |
| Yan et al. | Preliminary study on the quick detection of acquired immure deficiency syndrome by saliva analysis using surface enhanced Raman spectroscopic technique | |
| CN113533178B (zh) | 多物理特性融合感知的细胞流式检测方法 | |
| Mi et al. | High-precision high-speed nanopore ping-pong control system based on field programmable gate array |
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 |