CN113791018A - 基于电阻抗信号的单细胞机械本征参数测量系统和方法 - Google Patents

基于电阻抗信号的单细胞机械本征参数测量系统和方法 Download PDF

Info

Publication number
CN113791018A
CN113791018A CN202111063359.XA CN202111063359A CN113791018A CN 113791018 A CN113791018 A CN 113791018A CN 202111063359 A CN202111063359 A CN 202111063359A CN 113791018 A CN113791018 A CN 113791018A
Authority
CN
China
Prior art keywords
cell
flow channel
impedance
mechanical
cells
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
Application number
CN202111063359.XA
Other languages
English (en)
Other versions
CN113791018B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN202111063359.XA priority Critical patent/CN113791018B/zh
Publication of CN113791018A publication Critical patent/CN113791018A/zh
Application granted granted Critical
Publication of CN113791018B publication Critical patent/CN113791018B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/1031Investigating individual particles by measuring electrical or magnetic effects
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L3/00Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
    • B01L3/50Containers for the purpose of retaining a material to be analysed, e.g. test tubes
    • B01L3/502Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
    • B01L3/5027Containers 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L3/00Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
    • B01L3/50Containers for the purpose of retaining a material to be analysed, e.g. test tubes
    • B01L3/502Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
    • B01L3/5027Containers 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/50273Containers 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L3/00Containers or dishes for laboratory use, e.g. laboratory glassware; Droppers
    • B01L3/50Containers for the purpose of retaining a material to be analysed, e.g. test tubes
    • B01L3/502Containers for the purpose of retaining a material to be analysed, e.g. test tubes with fluid transport, e.g. in multi-compartment structures
    • B01L3/5027Containers 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/502761Containers 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/06Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L2200/00Solutions for specific problems relating to chemical or physical laboratory apparatus
    • B01L2200/10Integrating sample preparation and analysis in single entity, e.g. lab-on-a-chip concept
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01LCHEMICAL OR PHYSICAL LABORATORY APPARATUS FOR GENERAL USE
    • B01L2300/00Additional constructional details
    • B01L2300/08Geometry, shape and general structure
    • B01L2300/0887Laminated structure
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N2015/1006Investigating 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变化将能够反映细胞的机械特性;Ⅲ阶段,细胞完全变形并进入收缩流道并快速向前移动,该阶段的移动速度主要与细胞和流道的摩擦力有关。
所述由机械模型得到单细胞机械本征参数,在流道中产生细胞的受力变形,通过细胞变形过程和力学模型求取机械特性参数;为获得单细胞的机械本征参数,利用电阻抗测量的高时间分辨率特性,通过流道结构与电极形状的耦合实现对单细胞的形状与机械特性的测量,其关键在于电信号能够关联机械形变;对应于流道壁接触式挤压变形结构,将测量电极部署在收缩流道的入口,当细胞挤压变形进入收缩流道时将等体积置换测量区域内的溶液;此时细胞-培养基混合体系的阻抗可由下式计算
Figure BDA0003257281460000051
其中CDL为双电层电容,
Figure BDA0003257281460000052
为挤压细胞的阻抗,
Figure BDA0003257281460000053
为溶液阻抗,ω为激励信号的角频率。由于寄生电容Cs的存在,测量系统的整体阻抗为:
Figure BDA0003257281460000054
由于Cs的值一般在0.1-0.5pF之间,在激励信号频率为1MHz时,ω*Cs≈0,因此
Figure BDA0003257281460000055
而由电阻计算公式可知:
Figure BDA0003257281460000056
其中ρ为电阻率,L为长度,s为流道横截面,又由于Lcell+Lmed=Lchannel,可得:
Figure BDA0003257281460000057
由于Lchannel为设计常量,因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关,也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联。
所述计算求解细胞的机械本征参数,需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化,在Ⅱ阶段中,细胞的剪切蠕变J(t)由下式
Figure BDA0003257281460000061
其中Rp是流道的有效直径,可由
Figure BDA0003257281460000062
计算,Φp是一个常数,在矩形截面的微流道中通常取值为2.1,Lp(t)是细胞在t时刻的挤压伸长量,ΔP为施加压强;
根据幂律流变模型,瞬时剪切蠕变J(t)也可表示为:
Figure BDA0003257281460000063
其中AJ是反映细胞变形能力的剪切柔量常数,α是幂律指数也称为流度,t0一般取t0=1s;
细胞的杨氏模量AG可由下式计算:
Figure BDA0003257281460000064
其中Γ为伽马函数,因此由公式(5)-(7)可得
Figure BDA0003257281460000065
因此,细胞的机械本征参数AG和α可以通过最小化细胞的测量挤压伸长长度Lcell(t)及其模型估计值Lp(t)之间的方差来获得
Figure BDA0003257281460000066
其中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)所示,细胞-培养基混合体系的阻抗可由下式计算
Figure BDA0003257281460000091
其中CDL为双电层电容,
Figure BDA0003257281460000092
为挤压细胞的阻抗,
Figure BDA0003257281460000093
为溶液阻抗,ω为激励信号的角频率。由于寄生电容Cs的存在,测量系统的整体阻抗为:
Figure BDA0003257281460000094
由于Cs的值一般在0.1-0.5pF之间,在激励信号频率为1MHz时,ω*Cs≈0,因此
Figure BDA0003257281460000095
而由电阻计算公式可知:
Figure BDA0003257281460000096
其中ρ为电阻率,L为长度,s为流道横截面,又由于Lcell+Lmed=Lchannel,可得:
Figure BDA0003257281460000101
由于Lchannel为设计常量,因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关,也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联。
本发明中的机械本征参数测量原理是基于广泛使用的单细胞幂律流变模型,如图2中(a)与图3中(b)所示,细胞的挤压伸长长度Lp定义为细胞前沿与收缩入口的距离,当细胞挤压通过收缩流道时存在三个阶段:I阶段,细胞刚进入收缩区,Lp快速非线性增加;Ⅱ阶段,细胞开始收缩进入收缩流道并发生弹性变形,导致细胞进入收缩流道的速度降低并呈幂律变化,从而该阶段的Lp变化将能够反映细胞的机械特性;Ⅲ阶段,细胞完全变形并进入收缩流道并快速向前移动,该阶段的移动速度主要与细胞和流道的摩擦力有关。因此,为了计算求解细胞的机械本征参数,需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化。在Ⅱ阶段中,细胞的剪切蠕变J(t)由下式给出:
Figure BDA0003257281460000102
其中Rp是流道的有效直径,可由
Figure BDA0003257281460000103
计算,Φp是一个常数,在矩形截面的微流道中通常取值为2.1,Lp(t)是细胞在t时刻的挤压伸长量,ΔP为施加压强。
根据幂律流变模型,瞬时剪切蠕变J(t)也可表示为:
Figure BDA0003257281460000104
其中AJ是反映细胞变形能力的剪切柔量常数,α是幂律指数也称为流度,t0一般取t0=1s。
细胞的杨氏模量AG可由式(7)表示,
Figure BDA0003257281460000111
其中Γ为伽马函数,因此由式(5)-(7)可得
Figure BDA0003257281460000112
因此,细胞的机械本征参数AG和α可以通过最小化细胞的测量挤压伸长长度Lcell(t)及其模型估计值Lp(t)之间的方差来获得
Figure BDA0003257281460000113
其中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 (8)

1.一种基于电阻抗信号的单细胞机械本征参数测量系统,其特征在于,所述基于电阻抗信号的单细胞机械本征参数测量系统由微流控系统、多频阻抗测量系统和实时处理算法系统串联组成:
所述微流控系统主要由微流控器件和进样驱动装置组成,负责完成驱动大量单细胞聚束列队高速通过测量区域;
所述阻抗测量系统主要由直接数字频率合成(DDS)信号发生模块、锁定放大器模块和信号采集模块组成,负责测量并采集单细胞的阻抗值;
所述实时处理算法系统主要由数字滤波、事件检测与有效信号提取算法组成,负责实现对阻抗信号进行滤波、单细胞事件检测与细胞有效形变区域信号提取以及实时获取待测单细胞杨氏模量等机械本征参数;
所述微流控芯片结构采用收缩流道的尺寸设计,测量电极部署在收缩流道的入口与出口处,以使单细胞挤压通过;当细胞挤压通过收缩流道时,测量的阻抗幅度将与细胞的挤压伸长长度Lp成正比;因此,细胞进入收缩流道的形变过程能够通过阻抗信号来描述。
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)其弹性形变由阻抗测量系统采集;
5)实时求解本征参数,电阻抗信号实时解算,进而由机械模型得到单细胞机械本征参数。
6.根据权利要求5所述基于电阻抗信号的单细胞机械本征参数测量系统的测量方法,其特征在于,所述利用电阻抗测量信号的高时间分辨率特性来实现对单细胞的高通量单细胞机械本征参数测量,该单细胞机械本征参数测量广泛使用的单细胞幂律流变模型,单细胞的挤压伸长长度Lp定义为细胞前沿与收缩入口的距离,当细胞挤压通过收缩流道时存在三个阶段:I阶段,细胞刚进入收缩区,Lp快速非线性增加;Ⅱ阶段,细胞开始收缩进入收缩流道并发生弹性变形,导致细胞进入收缩流道的速度降低并呈幂律变化,从而该阶段的Lp变化将能够反映细胞的机械特性;Ⅲ阶段,细胞完全变形并进入收缩流道并快速向前移动,该阶段的移动速度主要与细胞和流道的摩擦力有关。
7.根据权利要求5所述基于电阻抗信号的单细胞机械本征参数测量系统的测量方法,其特征在于,所述由机械模型得到单细胞机械本征参数,在流道中产生细胞的受力变形,通过细胞变形过程和力学模型求取机械特性参数;为获得单细胞的机械本征参数,利用电阻抗测量的高时间分辨率特性,通过流道结构与电极形状的耦合实现对单细胞的形状与机械特性的测量,其关键在于电信号能够关联机械形变;对应于流道壁接触式挤压变形结构,将测量电极部署在收缩流道的入口,当细胞挤压变形进入收缩流道时将等体积置换测量区域内的溶液;此时细胞-培养基混合体系的阻抗可由下式计算
Figure FDA0003257281450000031
其中CDL为双电层电容,
Figure FDA0003257281450000032
为挤压细胞的阻抗,
Figure FDA0003257281450000033
为溶液阻抗,ω为激励信号的角频率。由于寄生电容Cs的存在,测量系统的整体阻抗为:
Figure FDA0003257281450000034
由于Cs的值一般在0.1-0.5pF之间,在激励信号频率为1MHz时,ω*Cs≈0,因此
Figure FDA0003257281450000041
而由电阻计算公式可知:
Figure FDA0003257281450000042
其中ρ为电阻率,L为长度,s为流道横截面,又由于Lcell+Lmed=Lchannel,可得:
Figure FDA0003257281450000043
由于Lchannel为设计常量,因此Roverall的大小与细胞的挤压伸长长度Lcell呈正相关,也就是细胞通过收缩流道过程中的挤压变形与阻抗大小形成关联。
8.根据权利要求5所述基于电阻抗信号的单细胞机械本征参数测量系统的测量方法,其特征在于,所述计算求解细胞的机械本征参数,需要提取细胞挤压通过收缩流道的Ⅱ阶段Lp动态变化,在Ⅱ阶段中,细胞的剪切蠕变J(t)由下式给出:
Figure FDA0003257281450000044
其中Rp是流道的有效直径,可由
Figure FDA0003257281450000045
计算,Φp是一个常数,在矩形截面的微流道中通常取值为2.1,Lp(t)是细胞在t时刻的挤压伸长量,ΔP为施加压强;
根据幂律流变模型,瞬时剪切蠕变J(t)也可表示为:
Figure FDA0003257281450000046
其中AJ是反映细胞变形能力的剪切柔量常数,α是幂律指数也称为流度,t0一般取t0=1s;
细胞的杨氏模量AG可由下式计算:
Figure FDA0003257281450000051
其中Γ为伽马函数,因此由公式(5)-(7)可得
Figure FDA0003257281450000052
因此,细胞的机械本征参数AG和α可以通过最小化细胞的测量挤压伸长长度Lcell(t)及其模型估计值Lp(t)之间的方差来获得
Figure FDA0003257281450000053
其中i表示连续的时间点,Lcell(t)为测量值,由公式(4)获得;Lp(t)为模型估计值,基于网格搜索与最小二乘拟合方法,从参数网格中搜索AG和α代入公式(8)中获得。
CN202111063359.XA 2021-09-10 2021-09-10 基于电阻抗信号的单细胞机械本征参数测量系统和方法 Active CN113791018B (zh)

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 true CN113791018A (zh) 2021-12-14
CN113791018B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115417585A (zh) * 2022-09-22 2022-12-02 重庆理工大学 一种多收缩孔径玻璃毛细管的加工装置及加工方法

Citations (8)

* Cited by examiner, † Cited by third party
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

Patent Citations (8)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
ANDREA ADAMO等: "Microfluidics-Based Assessment of Cell Deformability", 《ANALYTICAL CHEMISTRY》 *
JIAN CHEN等: "Classification of cell types using a microfluidic device for mechanical and electrical measurement on single cells", 《LAB ON A CHIP》 *
XIAONAN YANG等: "High-throughput and label-free parasitemia quantification and stage differentiation for malaria-infected red blood cells", 《BIOSENSORS AND BIOELECTRONICS》 *
YING ZHOU等: "Characterizing Deformability and Electrical Impedance of Cancer Cells in a Microfluidic Device", 《ANALYTICAL CHEMISTRY》 *
冯迪等: "微流控阻抗流式细胞仪在单细胞检测中的应用", 《化学进展》 *
唐文来等: "基于微流控技术的单细胞生物物理特性表征", 《化学进展》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115417585A (zh) * 2022-09-22 2022-12-02 重庆理工大学 一种多收缩孔径玻璃毛细管的加工装置及加工方法
CN115417585B (zh) * 2022-09-22 2023-08-29 重庆理工大学 一种多收缩孔径玻璃毛细管的加工装置及加工方法

Also Published As

Publication number Publication date
CN113791018B (zh) 2022-07-15

Similar Documents

Publication Publication Date Title
CN104251810B (zh) 同时表征单细胞杨氏模量和细胞膜比电容的系统
CN101726585A (zh) 一种基于微流控芯片的流式细胞仪
CN201548547U (zh) 一种基于微流控芯片的流式细胞分析装置
US20090280518A1 (en) System for high throughput measurement of mechanical properties of cells
CN103439241A (zh) 单细胞多参数表征的微流控芯片检测系统
Sobek Microfabricated fused silica flow chambers for flow cytometry
CN113567326A (zh) 一种高通量实时单细胞电学本征参数测量系统及方法
CN103923825A (zh) 一种集成细胞分选及检测的微流控芯片系统
CN113791018B (zh) 基于电阻抗信号的单细胞机械本征参数测量系统和方法
CN109107621A (zh) 基于细胞形变量和介电泳力的癌细胞分离装置及控制系统
CN113008766B (zh) 一种基于微流控多次挤压的癌细胞动态行为检测系统
CN112268934B (zh) 一种针对循环肿瘤细胞的检测芯片及其检测方法
CN111359687B (zh) 基于电学及荧光信号的微流控芯片和血细胞分析装置
CN113533171A (zh) 基于深度学习和微流控芯片的细胞形变检测方法及系统
Carminati et al. Design and characterization of a current sensing platform for silicon-based nanopores with integrated tunneling nanoelectrodes
Gajasinghe et al. Label and immobilization free detection and differentiation of tumor cells
Zhang et al. Tunable microfluidic chip for single-cell deformation study
US20240248074A1 (en) Method and system for single-cell biophysical profiling using a microfluidic device
WO2020097603A2 (en) Microfluidic antibody microarray with an electronic sensor array
CN108279202B (zh) 一种检测精度可调的颗粒计数装置和方法
CN101231262B (zh) 用于单细胞组学的方法
CN107515297B (zh) 一种自驱动自传感微悬臂传感器、制作方法及其应用
CN113413148A (zh) 一种提高阻抗生物传感器效果的信号处理方法及装置
Mi et al. High-precision high-speed nanopore ping-pong control system based on field programmable gate array
CN208155808U (zh) 一种利用电场聚焦提高颗粒检测精度和通量的装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant