CN113567326A - 一种高通量实时单细胞电学本征参数测量***及方法 - Google Patents
一种高通量实时单细胞电学本征参数测量***及方法 Download PDFInfo
- Publication number
- CN113567326A CN113567326A CN202110821829.8A CN202110821829A CN113567326A CN 113567326 A CN113567326 A CN 113567326A CN 202110821829 A CN202110821829 A CN 202110821829A CN 113567326 A CN113567326 A CN 113567326A
- Authority
- CN
- China
- Prior art keywords
- cell
- impedance
- real
- time
- measurement
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 34
- 210000004027 cell Anatomy 0.000 claims abstract description 144
- 238000012545 processing Methods 0.000 claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 238000002847 impedance measurement Methods 0.000 claims abstract description 17
- 239000012528 membrane Substances 0.000 claims abstract description 9
- 210000000805 cytoplasm Anatomy 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 210000000170 cell membrane Anatomy 0.000 claims abstract description 6
- 238000005516 engineering process Methods 0.000 claims description 22
- 239000000243 solution Substances 0.000 claims description 20
- 238000013528 artificial neural network Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 11
- 239000004205 dimethyl polysiloxane Substances 0.000 claims description 10
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 claims description 10
- 239000011521 glass Substances 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 239000006285 cell suspension Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000002156 mixing Methods 0.000 claims description 5
- 210000002569 neuron Anatomy 0.000 claims description 5
- -1 Polydimethylsiloxane Polymers 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 4
- 239000003990 capacitor Substances 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 230000001086 cytosolic effect Effects 0.000 claims description 3
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 210000004962 mammalian cell Anatomy 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 229920001992 poloxamer 407 Polymers 0.000 claims description 3
- 229920000642 polymer Polymers 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000002174 soft lithography Methods 0.000 claims description 3
- 239000004094 surface-active agent Substances 0.000 claims description 3
- 239000000725 suspension Substances 0.000 claims description 3
- 238000003491 array Methods 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 230000005611 electricity Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 9
- 238000012512 characterization method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 5
- 238000000157 electrochemical-induced impedance spectroscopy Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 230000005684 electric field Effects 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000005206 flow analysis Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 239000007850 fluorescent dye Substances 0.000 description 2
- 238000001453 impedance spectrum Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000026676 system process Effects 0.000 description 2
- 102000004310 Ion Channels Human genes 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000007853 buffer solution Substances 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000012258 culturing Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001215 fluorescent labelling Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 201000004792 malaria Diseases 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000009832 plasma treatment Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 230000008844 regulatory mechanism Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 208000007056 sickle cell anemia Diseases 0.000 description 1
- 238000001228 spectrum 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
- B01L2200/00—Solutions for specific problems relating to chemical or physical laboratory apparatus
- B01L2200/06—Fluid handling related problems
- B01L2200/0605—Metering of fluids
-
- 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/0809—Geometry, shape and general structure rectangular shaped
- B01L2300/0819—Microarrays; Biochips
-
- 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/01—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials specially adapted for biological cells, e.g. blood cells
-
- 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)
- Dispersion Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Hematology (AREA)
- Clinical Laboratory Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本发明公开了一种高通量实时单细胞电学本征参数测量***和方法,所述测量***包括微流控***、多频阻抗测量***和实时处理算法***共三个子***。其中,微流控***负责完成驱动大量单细胞聚束列队高速通过测量区域;多频阻抗测量***负责测量并采集单细胞在多个不同频点下的阻抗值;实时处理算法***负责实现对阻抗值进行滤波、脉冲检测提取以及实时获取待测单细胞的细胞半径、细胞质电导率和细胞膜单位膜电容等本征参数。应用本发明所述的测量***可以高通量实时单细胞电学本征参数,在生命科学仪器上具有广阔应用前景。
Description
技术领域
本发明单细胞电学参数测量领域,具体涉及一种高通量实时单细胞电学本征参数测量***及方法。
背景技术
单细胞表征提供了细胞基本结构,功能信息和病理状态,在揭示细胞异质性中起着重要作用,对生命科学研究、疾病诊断和个性化医学意义重大。一直以来,荧光标记是单细胞分析和成像的主要工具,它通过对细胞标记具有分子特异性的检测标签(如荧光染料、量子点、磁珠和稳定同位素等)来识别细胞中的组分及其状态。对细胞进行荧光标记不但需要细胞特异性的先验知识,而且标记过程的侵入式操作会改变细胞状态,使分析过程复杂化,限制后续分析。相比之下,细胞的生物物理特性(如电学特性和机械特性)同样与细胞内的分子组成有关,已被证明是诊断疾病(癌症、疟疾、糖尿病和镰状细胞性贫血等)的有效生物标记,且生物物理特性表征无需标记细胞,表征过程中细胞状态基本不变,表征后细胞仍然可以进行下一步操作和分析,如分选、培养、组学分析等,因此广受关注。
作为生物物理特性之一,细胞电学特性反映了膜形态、离子通道状态、核大小和细胞质的特征。电学特性表征的基本原理是利用电旋转、电阻抗谱和阻抗流式术,对细胞施加交流电激励,以细胞运动学或电学测量信号表征细胞的电学特性。电旋转技术中,单细胞被定位于微流控器件的多个电极中,受介电作用而发生旋转,通过拟合细胞的转速-电信号频率曲线,解算出单细胞的电学本征参数(如细胞大小,细胞质电导率和单位膜电容)。此技术对细胞的定位操作过程繁琐,测量时间较长,通量只有1个细胞/秒,且不太适用于高电导率的真实细胞缓冲液。电阻抗谱技术中,单细胞被捕获在电极中,对电极施加扫频信号并测量输出电流,得到单细胞的阻抗谱后从细胞电学模型中解算出细胞的电学本征参数。由于电信号响应快,该技术能够实现对捕获细胞的快速扫描,但是测量结果易受捕获结构的影响,通量因为捕获操作而无法提升。相比之下,阻抗流式技术中,单细胞连续高速通过电极检测区,直接以细胞的少数频点的绝对或相对阻抗值来表征其电学特性,通量高达1000个细胞/秒,但现有的基于模型拟合的本征参数解算方法耗时较长,测量***在线使用时只能给出现象学参数(如阻抗幅度、相位和阻抗透明度)来表征单个细胞,导致测量结果高度依赖于测量平台,并且无法在不同平台之间进行比较。
现有技术至少存在以下缺点:1)由于单细胞操作的局限性如对细胞的定位与捕获操作过程繁琐、测量时间较长,导致测量***复杂且单细胞测量通量较低;2)阻抗流式分析技术中单细胞连续高速通过电极检测区,单频点测量导致信息不足,只能以现象学参数(如阻抗幅度、相位和阻抗透明度)来表征单个细胞,使得测量结果高度依赖于测量平台,并且无法在不同平台之间进行比较;3)在高通量背景下,传统的基于梯度下降拟合或决策树等算法计算量大、所需时间长,无法满足实时处理的需求,导致现有技术无法实现对单细胞进行高通量实时在线的电学本征参数测量。
因此,需要寻找一种高通量实时单细胞电学本征参数测量***及方法,可以提高高通量场景下单细胞阻抗流式技术有效分析原始阻抗数据的能力,实现对单细胞的高通量实时在线电学本征参数测量。
发明目的
本发明的目的就在于解决现有技术中所存在的不足,提供一种高通量实时单细胞电学本征参数测量***及方法,一是利用微流控技术、阻抗流式分析技术与锁定放大技术,提高单细胞电学表征***的通量的同时降低测量***的复杂程度;二是利用多频锁定放大原理,得到单细胞的多频阻抗并解算出电学本征参数;三是利用神经网络技术快速分析复杂的输入数据并以高精度预测输出的优势,实现对单细胞进行高通量实时在线电学本征参数测量。
发明内容
根据本发明的一个方面,提供了一种高通量实时单细胞电学本征参数测量***,包括共三个子***:微流控***、多频阻抗测量***和实时处理算法***;
其中,所述微流控***由微流控器件和进样驱动装置组成,负责完成驱动大量单细胞聚束列队高速通过测量区域;
所述多频阻抗测量***由DDS信号发生模块、锁定放大器模块和信号采集模块组成,负责测量并采集单细胞在多个不同频点下的阻抗值;
所述实时处理算法***由数字滤波算法、事件检测与脉冲提取算法和神经网络回归算法组成,负责实现对阻抗值进行滤波、脉冲检测提取以及实时获取待测单细胞的细胞半径、细胞质电导率和细胞膜单位膜电容这些本征参数。
优选地,所述微流控***中的微流控器件被芯片化设计为微流控器件芯片,所述微流控器件芯片被设计为两层,其中,底层以玻璃作为基底并配置有一对平面电极以进行阻抗检测,顶层配置有聚二甲基硅氧烷PDMS微通道,以使单细胞聚束列队依次通过检测区域,所述检测区域流道的尺寸由细胞大小决定。
更优选地,对于直径为8-18μm的哺乳动物细胞,所述测量区域流道的宽度和高度均为20μm;在所述测量区域的入口配置有间隔为20μm的柱状阵列作为片上过滤器,并根据伯努利原理,将所述过滤器设计为漏斗状;所述一对平面电极被配置为30μm长,20μm宽和20μm间距。
优选地,所述微流控器件芯片的微流道层通过软光刻技术制作,利用PDMS倒模形成;所述微流控器件芯片的顶层的平面电极层是通过lift-off技术图案化到玻璃基底上,用氧等离子体处理并烘烤之后,将所述平面电极层和微通道层对准牢固地键合在一起,接着,所述平面电极与其图案匹配的定制印刷电路板焊接在一起。
优选地,所述微流控***中,细胞悬浮液进口出口通过塑料软管与流体驱动装置相连,通过流体驱动装置将细胞悬浮液推注入所述微流控器件芯片顶层的微流道。
优选地,所述DDS信号发生模块为基于DDS芯片AD9958搭建的四频信号发生器;所述锁定放大器模块包括基于电流运算放大器芯片OPA657搭建的跨阻放大器TIA、基于宽带四象限电压输出乘法器AD835芯片搭建的解调器、具有5KHz截止频率的低通滤波器和由运算放大器OPA227搭建的放大器;所述信号采集模块为连接计算机的数据采集卡;所述多频阻抗测量***最终输出的直流信号由计算机通过所述数据采集卡以20KHz的采样率进行采样获取,并在计算机上由时处理算法***进行后续处理;
优选地,所述实时处理算法***基于Python,Matlab和LabView混合编程,其处理过程为:首先通过数字低通滤波器对数据流进行数字滤波,以进一步提高信噪比SNR;随后在Matlab脚本中检测单细胞事件并提取相应的单细胞阻抗峰值;所提取的阻抗峰值通过TCP/IP协议发送到Python节点,以通过训练好的神经网络计算单细胞的本征参数,并将计算结果显示在图形界面中,整个处理过程在实时在线模式下进行。
根据本发明的另一个方面,提供了一种应用上述高通量实时单细胞电学本征参数测量***实时测量单细胞电学本征参数的方法,包括以下步骤:
步骤S1:流道修饰,具体是将所述微流控器件芯片的微流道用1wt%,即1x PBS的Pluronic F-127表面活性剂预处理15分钟;
步骤S2:细胞进给,具体是使用注射泵将单细胞样品以1μL/min的恒定流速通入微通道中,其中,在测量区域的等效速度为20mm/s,从而实现超过1000个细胞/分钟的高通量;
步骤S3:信号测量,具体是当单细胞在流道中列队通过测量区域时,所产生的电信号由阻抗测量***采集并经过实时处理算法***求解得到单细胞本征参数;
步骤S4:将步骤S3中所得到的单细胞本征参数实时显示在用户界面上。
优选地,步骤S3中阻抗测量***采集电信号并经过实时处理算法***求解得到单细胞本征参数的过程包括:
子步骤S31:用单细胞和悬浮介质的特性来表示细胞-溶液混合体系的电学特性,其复数阻抗由式(1)给出:
子步骤S32:构建当单细胞位于所述微流控器件芯片的一对平面测量电极之间时的等效电路模型ECM,其中,阻抗由溶液电阻Rmed、溶液电容Cmed、细胞质电阻Ri、细胞膜电容Cmem、双电层电容CDL和杂散电容Cs组成,将细胞-溶液混合体系的阻抗表征为如式(4)所示:
考虑到Cs和CDL,将模型总阻抗表示为如式(5)所示:
子步骤S33:分别以45KHz,250KHz,750KHz和1200KHz四个频率进行模型拟合和参数提取,其中45KHz和250KHz用于确定CDL和细胞大小,750KHz用于确定细胞膜介电常数,1200KHz用于确定细胞质电导率和Cs;利用最小二乘法与梯度下降法,通过最小化测量阻抗值与模型估计值,及模型总阻抗之间的方差确定单细胞的本征参数如式(6)所示:
步骤S34:基于混频原理,将步骤S33中的四个频率点的余弦信号通过加法器进行调制,表示为如式(7)所示:
S(t)=cosω1t+cosω2t+cosω3t+cosω4t (7),
再分别通过四个解调器将响应信号R(t)解调为实电压和虚电压,如式(8)所示:
利用跨阻放大器TIA中使用的已知参考反馈电阻RF和测量***的校准系数,从而计算出复数阻抗。
优选地,所述实时处理算法***的处理过程包括以下序列:
序列1:在开源框架Pytorch下搭建了一个深层全连接神经网络FCN用于处理阻抗数据流,从而实现对单细胞本征参数的实时求解;所述设计的FCN由五层组成,第一层为输入层,负责接收输入数据;中间三层为隐藏层,负责将输入数据映射到输出数据;最后一层是输出层,负责导出输出结果;所述FCN的每一层都由多个神经元组成,这些神经元获得前一层的加权总和,并利用线性整流ReLU激活函数将加权总和转换为输出;
序列2:将FCN的输入设为四频点的单细胞阻抗幅度,输出则为本征参数,选择均方误差函数作为损失函数计算标签值与输出值的偏差,并通过使用随机梯度下降SGD优化器反向传播误差对FCN进行监督训练;通过基于梯度下降算法模型拟合获得的本征参数估计值,用作训练数据的目标标签;完整的数据集分为占70%的训练集与占30%的测试集两个数据集,分别用于训练和测试。
附图说明
图1是本发明所述高通量实时单细胞电学本征参数实时测量***示意框图。
图2是测量***中微流控***器件的结构示意图:(a)微流控器件测结构示意图;(b)微流控器件测量尺寸图。
图3是测量***中微流控***器件加工示意图及示例图:(a)微流控器件加工示意图;(b)一种器件实例图。
图4是单细胞等效电路模型及电阻抗谱:(a)单细胞等效电路模型;(b)单细胞电阻抗谱
图5是单细胞通过测量电极产生的电信号示意图以及锁定放大器原理图:(a)单细胞通过测量电极产生的电信号示意图;(b)锁定放大器原理图。
图6是测量***中多频阻抗测量硬件***视为示意图及软件***处理流程图:(a)多频阻抗测量硬件***示意图;(b)软件***处理流程图。
图7是本发明实施例所应用的神经网络结构示意图。
图8是本发明所述高通量实时单细胞电学本征参数实时测量***的使用流程框图。
具体实施方式
以下结合附图详细阐述本发明的具体实施方式,本领域技术人员应当理解,此处的具体实施例仅用于示例性地描述本发明,而不应被视作对本发明保护范围的限定。
图1是本发明所述高通量实时单细胞电学本征参数实时测量***示意框图。由图可知,所述高通量实时单细胞电学本征参数实时测量***包括共三个子***:微流控***、多频阻抗测量***和实时处理算法***。以下详细描述各个子***的结构设计。
微流控***及器件结构设计
由图1可以看出,微流控***由微流控器件与进样装置组成。图2是测量***中微流控***器件的结构示意图,如图2所示,考虑到制造和光学观察的便利性,本实施例中微流控器件芯片被设计为两层,结构示意图如图2(a)所示。底层以玻璃作为基底并配置有一对平面电极以进行阻抗检测,顶层则配置有聚二甲基硅氧烷(PDMS)微通道以使单细胞聚束列队依次通过检测区域,检测区域流道的尺寸由细胞大小决定。特别地,针对常见的哺乳动物细胞,在典型直径为8-18μm的范围时,为了实现高灵敏度且准确的单细胞电学表征并避免细胞阻塞和电流泄漏问题,如图2(b)所示,测量区域流道的宽度和高度均为20μm。同时为了降低发生多个细胞同时通过测量区域的可能性,测量区域的入口配置有间隔为20μm的柱状阵列作为片上过滤器,并根据伯努利原理设计为漏斗状从而产生压力分离距离太近的细胞。此外,通过对单细胞电学模型的分析优化的电极尺寸为30μm长,20μm宽和20μm间距。需要指出,给出的这个结构和尺寸并不是唯一的,可以根据应用场合进行适当调整。
图3是测量***中微流控***器件加工示意图及示例图。如图3所示,微流控器件的一种制造过程如图3(a)所示。微流道层通过软光刻技术制作,利用PDMS倒模形成,电极层通过lift-off技术图案化到玻璃基底上。在用氧等离子体处理并烘烤之后,将电极层和微通道层对准牢固地键合在一起。最终,为便于施加激励信号和采集测量信号,电极与其图案匹配的定制印刷电路板焊接在一起。图3(b)为器件的实例图,细胞悬浮液进口出口通过塑料软管与流体驱动装置相连,通过流体驱动装置将细胞悬浮液推注入微流道中。
单细胞电学参数测量原理
本实施例中的电学本征参数测量原理可以球形单壳模型为一种实例,实际上,任何模型都可以用本***来实现。在球形单壳模型中,细胞可以表示为被电介质膜包裹的均匀导电细胞质。当悬浮的细胞处于交流电场中时,通过麦克斯韦的混合理论(MMT),可以用单细胞和悬浮介质的特性来表示细胞-溶液混合体系的电学特性,其复数阻抗由式(1)给出:
图4是单细胞等效电路模型及电阻抗谱,当细胞位于一对测量电极之间时,如图4(a)所示可以推导出等效电路模型ECM。在最简单的分析模型中,阻抗由溶液电阻Rmed,溶液电容Cmed,细胞质电阻Ri,细胞膜电容Cmem,双电层电容CDL和杂散电容Cs组成。利用这些参数,计算细胞-溶液混合体系的阻抗如式(4)所示:
考虑到Cs和CDL,模型总阻抗表示为如式(5)所示:
如图4(b)所示,计算出的电阻抗谱具有四个区域。在第一个区域(<0.1MHz)时,阻抗由CDL主导,导致细胞-溶液混合体系的阻抗被淹没在高阻抗中;随着频率的增加,CDL的影响在第二个区域(0.1MHz–0.5MHz)处呈指数下降,介电膜开始被电场穿透导致阻抗由膜电容调节,该膜的介电常数通常为10mF/m2;当频率在第三个区域(0.5MHz–5MHz)处连续增加时,介电膜被电场完全穿透,此时细胞质的电导率在阻抗谱中占主导地位;但是,第四个区域(>5MHz)的阻抗受Cs的影响,导致难以探索细胞核的性质。图4(b)展示了单细胞的电阻抗如何随本征参数的变化而变化,并为多频阻抗测量的频率选择提供了指导。根据上述讨论,本发明方法以其中四个频率(即45KHz,250KHz,750KHz和1200KHz)为例进行模型拟合和参数提取,其中45KHz和250KHz用于确定CDL和细胞大小,750KHz用于确定细胞膜介电常数,1200KHz用于确定细胞质电导率和Cs。利用最小二乘法与梯度下降法,通过最小化测量阻抗值与模型估计值(模型总阻抗)之间的方差确定单细胞的本征参数如式(6)所示:
多频阻抗测量***设计与原理
细胞穿过测量区域时将替换等体积的介质并扰动电场,从而导致电路中阻抗的变化。图5是单细胞通过测量电极产生的电信号示意图以及锁定放大器原理图。如图5(a)所示,该信号具有高斯形状,其中的峰值对应于特定频率下的细胞阻抗。如图5(b)所示,阻抗测量的基本原理是将电流信号通过跨阻放大器TIA转换为电压信号,并通过锁相放大器LIA检测所得的微弱响应信号,从而实现对单细胞阻抗的测量。本实施例基于混频原理单细胞测量原理中的四频点余弦信号通过加法器进行调制,表示为如式(7)所示:
S(t)=cosω1t+cosω2t+cosω3t+cosω4t (7),
随后分别通过四个解调器将响应信号R(t)解调为实电压和虚电压,如式(8)所示:
利用跨阻放大器TIA中使用的已知参考反馈电阻RF和测量***的校准系数,计算出复数阻抗。
图6是测量***中多频阻抗测量硬件***示意图及软件***处理流程图,其中,多频阻抗测量硬件***如图6(a)所示,其用于负责产生激励信号,解调响应信号并获取输出直流信号。作为一种实例,四频信号发生器基于DDS芯片AD9958(ADI,美国)搭建,跨阻放大器TIA基于高增益带宽、低偏置电流运算放大器芯片OPA657(TI,美国)搭建,解调器基于宽带四象限电压输出乘法器AD835芯片(ADI,美国)搭建和具有5KHz截止频率的低通滤波器和放大器由高精度,低噪声运算放大器OPA227(TI,美国)搭建。最终输出的直流信号由计算机通过数据采集卡(NI,PCI-6289)以20KHz的采样率进行采样获取,并在计算机上由软件***进行后续处理。软件***处理流程如图6(b)所示,用于负责实时处理,显示和记录数据。软件***基于Python,Matlab(Math Works,美国)和LabView(NI,美国)混合编程,首先通过数字低通滤波器对数据流进行数字滤波,以进一步提高信噪比(SNR),随后在Matlab脚本中检测单细胞事件并提取相应的单细胞阻抗峰值。提取的阻抗峰值通过TCP/IP协议发送到Python节点,以通过训练好的神经网络计算单细胞的本征参数,并将计算结果显示在图形界面中,整个工作流程都在实时在线模式下进行。
实时算法设计与原理
本实施例利用神经网络技术快速分析复杂的输入数据并以高精度预测输出的优势。图7是本实施例所应用的神经网络结构示意图,如图7所示,在开源框架Pytorch下搭建了一个深层全连接神经网络FCN用于处理阻抗数据流,从而实现对单细胞本征参数的实时求解。设计的FCN由五层组成,第一层为输入层(负责接收输入数据),中间三层为隐藏层(负责将输入数据映射到输出数据),最后一层是输出层(负责导出输出结果)。每一层都由多个神经元组成,这些神经元获得前一层的加权总和,并利用线性整流ReLU激活函数将加权总和转换为输出。ReLU激活函数将非线性特性引入了神经网络***。为了求解单细胞的本征参数,FCN的输入为四频点的单细胞阻抗幅度,输出则为本征参数,因此选择均方误差函数作为损失函数计算标签值与输出值的偏差,并通过使用随机梯度下降(SGD)优化器反向传播误差对FCN进行监督训练。通过传统的基于梯度下降算法模型拟合获得的本征参数估计值用作训练数据的目标标签。完整的数据集分为训练集(70%)与测试集(30%)两个数据集,分别用于训练和测试。在测试过程中,单细胞事件的推导时间约为0.3ms,而基于梯度下降的模型拟合方法则需要几秒钟才能解算得本征参数,从而使处理速度提高了数个数量级,从而实现了实时在线单个细胞的本征参数求解。
图8是本发明所述高通量实时单细胞电学本征参数实时测量***的使用流程框图。如图所示,使用所述测量***实时测量单细胞电学本证参数时,依次包括以下步骤:
步骤S1:流道修饰,具体是将微通道用1wt%,即1x PBS的Pluronic F-127表面活性剂预处理15分钟;
步骤S2:细胞进给,具体是使用注射泵将单细胞样品以1μL/min的恒定流速通入微通道中,其中,在测量区域的等效速度为20mm/s,从而实现超过1000个细胞/分钟的高通量;
步骤S3:信号测量,具体是当单细胞在流道中列队通过测量区域时,所产生的电信号由阻抗测量***采集并经过实时处理算法求解得到单细胞本征参数;
步骤S4:将步骤S3中所得到的单细胞本征参数实时显示在用户界面上。
在电学测量过程外,***使用过程中也可使用配备CCD相机的倒置显微镜监视和记录细胞的运动,并通过时间戳与电信号同步,从而可以离线分析以将电阻抗信号与微观观察相关联。
相较于现有技术,本发明具有如下优势:
1)利用微流控技术、阻抗流式分析技术与锁定放大技术,提高单细胞电学表征***的通量的同时降低测量***的复杂程度,构筑了一套样本入-结果出的单细胞表征方法和***;
2)利用多频锁定放大原理,通过测量四个频点的阻抗值基于梯度下降的模型拟合方法离线解算出电学本征参数;
3)利用神经网络技术快速分析复杂的输入数据并以高精度预测输出的优势,将解算时间由秒级下降至亚毫秒级,实现对单细胞进行高通量实时在线电学本征参数测量。
综上所述,本发明结合阻抗流式芯片结构设计、电学本征参数提取方法与人工智能与机器学***调控机制提供技术基础,推动单细胞分析向即时诊断和临床应用纵深发展。这种包含微流控技术、单细胞测量和人工智能和机器学***,对生命科学研究、疾病诊断和个性化医学具有重大意义。
Claims (10)
1.一种高通量实时单细胞电学本征参数测量***,其特征在于,包括三个子***:微流控***、多频阻抗测量***和实时处理算法***;
其中,所述微流控***由微流控器件和进样驱动装置组成,负责完成驱动大量单细胞聚束列队高速通过测量区域;
所述多频阻抗测量***由DDS信号发生模块、锁定放大器模块和信号采集模块组成,负责测量并采集单细胞在多个不同频点下的阻抗值;
所述实时处理算法***由数字滤波算法、事件检测与脉冲提取算法和神经网络回归算法组成,负责实现对阻抗值进行滤波、脉冲检测提取以及实时获取待测单细胞的细胞半径、细胞质电导率和细胞膜单位膜电容这些本征参数。
2.根据权利要求1所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,所述微流控***中的微流控器件被芯片化设计为微流控器件芯片,所述微流控器件芯片被设计为两层,其中,底层以玻璃作为基底并配置有一对平面电极以进行阻抗检测,顶层配置有聚二甲基硅氧烷PDMS微通道,以使单细胞聚束列队依次通过检测区域,所述检测区域流道的尺寸由细胞大小决定。
3.根据权利要求2所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,对于直径为8-18μm的哺乳动物细胞,所述测量区域流道的宽度和高度均为20μm;在所述测量区域的入口配置有间隔为20μm的柱状阵列作为片上过滤器,并根据伯努利原理,将所述过滤器设计为漏斗状;所述一对平面电极被配置为30μm长,20μm宽和20μm间距。
4.根据权利要求2所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,所述微流控器件芯片的微流道层通过软光刻技术制作,利用PDMS倒模形成;所述微流控器件芯片的顶层的平面电极层是通过lift-off技术图案化到玻璃基底上,用氧等离子体处理并烘烤之后,将所述平面电极层和微通道层对准牢固地键合在一起,接着,所述平面电极与其图案匹配的定制印刷电路板焊接在一起。
5.根据权利要求1-4任一所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,所述微流控***中,细胞悬浮液进口出口通过塑料软管与流体驱动装置相连,通过流体驱动装置将细胞悬浮液推注入所述微流控器件芯片顶层的微流道。
6.根据权利要求1所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,所述DDS信号发生模块为基于DDS芯片AD9958搭建的四频信号发生器;所述锁定放大器模块包括基于电流运算放大器芯片OPA657搭建的跨阻放大器TIA、基于宽带四象限电压输出乘法器AD835芯片搭建的解调器、具有5KHz截止频率的低通滤波器和由运算放大器OPA227搭建的放大器;所述信号采集模块为连接计算机的数据采集卡;所述多频阻抗测量***最终输出的直流信号由计算机通过所述数据采集卡以20KHz的采样率进行采样获取,并在计算机上由时处理算法***进行后续处理。
7.根据权利要求1所述的一种高通量实时单细胞电学本征参数测量***,其特征在于,所述实时处理算法***基于Python,Matlab和LabView混合编程,其处理过程为:首先通过数字低通滤波器对数据流进行数字滤波,以进一步提高信噪比SNR;随后在Matlab脚本中检测单细胞事件并提取相应的单细胞阻抗峰值;所提取的阻抗峰值通过TCP/IP协议发送到Python节点,以通过训练好的神经网络计算单细胞的本征参数,并将计算结果显示在图形界面中,整个处理过程。
8.一种应用根据权利要求1-7任一所述的一种高通量实时单细胞电学本征参数测量***实时测量单细胞电学本征参数的方法,其特征在于,包括以下步骤:
步骤S1:流道修饰,具体是将所述微流控器件芯片的微流道用1wt%,即1x PBS的Pluronic F-127表面活性剂预处理15分钟;
步骤S2:细胞进给,具体是使用注射泵将单细胞样品以1μL/min的恒定流速通入微通道中,其中,在测量区域的等效速度为20mm/s,从而实现超过1000个细胞/分钟的高通量;
步骤S3:信号测量,具体是当单细胞在流道中列队通过测量区域时,所产生的电信号由阻抗测量***采集并经过实时处理算法***求解得到单细胞本征参数;
步骤S4:将步骤S3中所得到的单细胞本征参数实时显示在用户界面上。
9.根据权利要求8所述的实时测量单细胞电学本征参数的方法,其特征在于,步骤S3中阻抗测量***采集电信号并经过实时处理算法***求解得到单细胞本征参数的过程包括:
子步骤S31:用单细胞和悬浮介质的特性来表示细胞-溶液混合体系的电学特性,其复数阻抗由式(1)给出:
子步骤S32:构建当单细胞位于所述微流控器件芯片的一对平面测量电极之间时的等效电路模型ECM,其中,阻抗由溶液电阻Rmed、溶液电容Cmed、细胞质电阻Ri、细胞膜电容Cmem、双电层电容CDL和杂散电容Cs组成,将细胞-溶液混合体系的阻抗表征为如式(4)所示:
考虑到Cs和CDL,将模型总阻抗表示为如式(5)所示:
子步骤S33:分别以45KHz,250KHz,750KHz和1200KHz四个频率进行模型拟合和参数提取,其中45KHz和250KHz用于确定CDL和细胞大小,750KHz用于确定细胞膜介电常数,1200KHz用于确定细胞质电导率和Cs;利用最小二乘法与梯度下降法,通过最小化测量阻抗值与模型估计值,及模型总阻抗之间的方差确定单细胞的本征参数如式(6)所示:
步骤S34:基于混频原理,将步骤S33中的四个频率点的余弦信号通过加法器进行调制,表示为如式(7)所示:
S(t)=cosω1t+cosω2t+cosω3t+cosω4t (7),
再分别通过四个解调器将响应信号R(t)解调为实电压和虚电压,如式(8)所示:
利用跨阻放大器TIA中使用的已知参考反馈电阻RF和测量***的校准系数,从而计算出复数阻抗。
10.根据权利要求8所述的实时测量单细胞电学本征参数的方法,其特征在于,所述实时处理算法***的处理过程包括以下序列:
序列1:在开源框架Pytorch下搭建了一个深层全连接神经网络FCN用于处理阻抗数据流,从而实现对单细胞本征参数的实时求解;所述设计的FCN由五层组成,第一层为输入层,负责接收输入数据;中间三层为隐藏层,负责将输入数据映射到输出数据;最后一层是输出层,负责导出输出结果;所述FCN的每一层都由多个神经元组成,这些神经元获得前一层的加权总和,并利用线性整流ReLU激活函数将加权总和转换为输出;
序列2:将FCN的输入设为四频点的单细胞阻抗幅度,输出则为本征参数,选择均方误差函数作为损失函数计算标签值与输出值的偏差,并通过使用随机梯度下降SGD优化器反向传播误差对FCN进行监督训练;通过基于梯度下降算法模型拟合获得的本征参数估计值,用作训练数据的目标标签;完整的数据集分为占70%的训练集与占30%的测试集两个数据集,分别用于训练和测试。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110821829.8A CN113567326A (zh) | 2021-07-19 | 2021-07-19 | 一种高通量实时单细胞电学本征参数测量***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110821829.8A CN113567326A (zh) | 2021-07-19 | 2021-07-19 | 一种高通量实时单细胞电学本征参数测量***及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113567326A true CN113567326A (zh) | 2021-10-29 |
Family
ID=78165834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110821829.8A Pending CN113567326A (zh) | 2021-07-19 | 2021-07-19 | 一种高通量实时单细胞电学本征参数测量***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113567326A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114112808A (zh) * | 2021-11-05 | 2022-03-01 | 国家纳米科学中心 | 细胞质力学性质的表征方法 |
CN114350490A (zh) * | 2021-12-30 | 2022-04-15 | 杭州电子科技大学 | 一种用于测量细胞s参数的检测平台及其检测方法 |
CN117487883A (zh) * | 2023-12-29 | 2024-02-02 | 嘉兴朝云帆生物科技有限公司 | 一种用于核酸分析的微流控检测方法和*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460114A (zh) * | 2009-06-05 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 用于鉴别和计数表达特定标志物的颗粒的多频阻抗方法和装置 |
CN103439241A (zh) * | 2013-08-23 | 2013-12-11 | 东南大学 | 单细胞多参数表征的微流控芯片检测*** |
CN104941704A (zh) * | 2015-05-27 | 2015-09-30 | 东南大学 | 一种集成细胞聚焦与检测的方法及其微型化*** |
CN109321456A (zh) * | 2018-10-12 | 2019-02-12 | 江苏大学 | 一种微流控芯片细胞培养控制装置与方法 |
CN109682745A (zh) * | 2019-01-07 | 2019-04-26 | 清华大学 | 一种单细胞参数测量方法及装置 |
CN111564183A (zh) * | 2020-04-24 | 2020-08-21 | 西北工业大学 | 融合基因本体和神经网络的单细胞测序数据降维方法 |
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 |
CN113029917A (zh) * | 2021-02-22 | 2021-06-25 | 中国科学院空天信息创新研究院 | 一种细胞及细胞核生物电学特性检测装置及方法 |
-
2021
- 2021-07-19 CN CN202110821829.8A patent/CN113567326A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460114A (zh) * | 2009-06-05 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 用于鉴别和计数表达特定标志物的颗粒的多频阻抗方法和装置 |
CN103439241A (zh) * | 2013-08-23 | 2013-12-11 | 东南大学 | 单细胞多参数表征的微流控芯片检测*** |
CN104941704A (zh) * | 2015-05-27 | 2015-09-30 | 东南大学 | 一种集成细胞聚焦与检测的方法及其微型化*** |
CN109321456A (zh) * | 2018-10-12 | 2019-02-12 | 江苏大学 | 一种微流控芯片细胞培养控制装置与方法 |
CN109682745A (zh) * | 2019-01-07 | 2019-04-26 | 清华大学 | 一种单细胞参数测量方法及装置 |
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 |
CN111564183A (zh) * | 2020-04-24 | 2020-08-21 | 西北工业大学 | 融合基因本体和神经网络的单细胞测序数据降维方法 |
CN113029917A (zh) * | 2021-02-22 | 2021-06-25 | 中国科学院空天信息创新研究院 | 一种细胞及细胞核生物电学特性检测装置及方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114112808A (zh) * | 2021-11-05 | 2022-03-01 | 国家纳米科学中心 | 细胞质力学性质的表征方法 |
CN114112808B (zh) * | 2021-11-05 | 2024-02-20 | 国家纳米科学中心 | 细胞质力学性质的表征方法 |
CN114350490A (zh) * | 2021-12-30 | 2022-04-15 | 杭州电子科技大学 | 一种用于测量细胞s参数的检测平台及其检测方法 |
CN114350490B (zh) * | 2021-12-30 | 2024-03-22 | 杭州电子科技大学 | 一种用于测量细胞s参数的检测平台及其检测方法 |
CN117487883A (zh) * | 2023-12-29 | 2024-02-02 | 嘉兴朝云帆生物科技有限公司 | 一种用于核酸分析的微流控检测方法和*** |
CN117487883B (zh) * | 2023-12-29 | 2024-03-29 | 嘉兴朝云帆生物科技有限公司 | 一种用于核酸分析的微流控检测方法和*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113567326A (zh) | 一种高通量实时单细胞电学本征参数测量***及方法 | |
CN103323383B (zh) | 一种基于电阻抗技术的微流控芯片微粒计数*** | |
Feng et al. | Neural network-enhanced real-time impedance flow cytometry for single-cell intrinsic characterization | |
US9915599B2 (en) | Microparticle analysis apparatus and microparticle analysis system | |
US8828212B2 (en) | Dielectric cytometric apparatus and dielectric-cytometric cell sorting method | |
CN106644900B (zh) | 一种基于非均匀电场的阻抗脉冲颗粒计数装置及其计数方法 | |
CN104713932B (zh) | 一种交流模式的多参数纳米孔单分子分析仪 | |
US11604133B2 (en) | Use of multi-frequency impedance cytometry in conjunction with machine learning for classification of biological particles | |
CN203337513U (zh) | 一种基于电阻抗技术的微流控芯片微粒计数*** | |
Caselli et al. | A Bayesian approach for coincidence resolution in microfluidic impedance cytometry | |
Sherif et al. | Integration of tri-polar microelectrodes for performance enhancement of an impedance biosensor | |
US20050114041A1 (en) | Method for discrimination of particles in a flow cytometer | |
Zhu et al. | An easy-fabricated and disposable polymer-film microfluidic impedance cytometer for cell sensing | |
Ramos et al. | LabVIEW 2010 computer vision platform based virtual instrument and its application for pitting corrosion study | |
Carminati et al. | Design and characterization of a current sensing platform for silicon-based nanopores with integrated tunneling nanoelectrodes | |
CN113791018B (zh) | 基于电阻抗信号的单细胞机械本征参数测量***和方法 | |
CN114870917B (zh) | 用于识别不同细胞的微流控芯片及其制备方法和检测平台 | |
CN114152557B (zh) | 基于图像分析的血细胞计数方法和*** | |
Wang et al. | Integrated sensor networks with error correction for multiplexed particle tracking in microfluidic chips | |
Ell et al. | Label-free identification of nonelectrogenic cancer cells using adhesion noise | |
CN115044468A (zh) | 一种用于单细胞电参数测量的流式电旋转微器件 | |
Tang et al. | An impedance microsensor with coplanar electrodes and vertical sensing apertures | |
US20230408396A1 (en) | Method and device for identifing the material type of particles in liquid | |
Baur et al. | Joint particle detection and analysis by a cnn and adaptive norm minimization approach | |
González-Murillo et al. | Electrical impedance spectroscopy microflow cytometer for cell viability tests |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211029 |
|
RJ01 | Rejection of invention patent application after publication |