一种流式细胞系统荧光寿命高斯拟合测量方法技术方案

本发明专利技术提供了一种流式细胞系统荧光寿命高斯拟合测量方法，所述方法包括以下步骤：1)利用时间微元法求出瞬时散射光强度和荧光信号强度，2)利用卷积法对峰值位置的变化进行分析，3)利用高斯拟合对散射光信号与荧光信号进行表示，得到的脉冲信号进行高斯拟合。

【技术实现步骤摘要】

本专利技术涉及荧光检测及光谱分析领域，具体涉及流式细胞仪的荧光强度检测及荧光寿命时域估算领域。
技术介绍
流式细胞仪是一种集激光技术、电子物理技术、光电测量技术、电子计算机技术、细胞荧光化学技术和单克隆抗体技术为一体的新型高科技仪器。对悬液中处于高速、直线流动的单细胞或其他颗粒，通过检测散射光信号和(或)标记的荧光信号，实现高速逐一的多参数定量分析。在细胞生物学、细胞周期动力学、免疫学、血液学及肿瘤学等领域具有广泛的应用。流式细胞仪对悬液中处于高速、直线流动的单细胞或其他颗粒，通过检测散射光信号和(或)标记的荧光信号，实现高速逐一的多参数定量分析。在细胞生物学、细胞周期动力学、免疫学、血液学及肿瘤学等领域具有广泛的应用。流式细胞仪实现前向散射光(FS)、侧向散射光(SS)及各色荧光信号(FLn)的收集及光电转换，并对转换后的脉冲信号进行多参数提取，最终根据提取得到的多参数实现对细胞/微球的统计学分析。荧光寿命测量能够提供有机荧光素、荧光蛋白和其他可发荧光的无机分子相关的荧光弛豫、强度衰减等信息。并且荧光寿命不存在荧光强度测量过程中存在的非线性问题，所以荧光寿命被引入到流式细胞仪单细胞测量中作为一个多参数信息。应用到流式分析的荧光寿命测量方法主要有频域方法和时域方法两种。频域方法系统中，通过高频正弦信号对激发光源进行调制，从而使激发出的荧光信号与激发光具有相同的调制频率，并且两者之间存在一定的相移。时域方法系统中，利用飞秒激光器作为激发光源对样品或细胞进行激发，通过单光子计数器及指数衰减模型对荧光衰减时间进行观测。频域方法和时域方法都会大大增加传统流式细胞仪的成本及复杂程度。2014年C.Ruofan等人提出利用数字信号处理方法对传统流式细胞仪荧光信号时延进行分析，在不改变流式细胞仪结构的情况下采用250Mbps高速ADC芯片对荧光信号进行采样并利用(Directmethod、GaussianFitting、HalfArea)3种方法对荧光寿命进行计算，计算结果的时域分辨能力受限于ADC的采样频率。
技术实现思路
为了解决上述问题，本专利技术的目的在于提供一种流式细胞系统荧光寿命高斯拟合测量方法，所述方法包括以下步骤：1)利用时间微元法求出瞬时散射光强度和荧光信号强度，由以下公式得出t时刻散射光强度：其中，R为荧光微球的半径，dx为微元的宽度，x为微元距位置原点的距离；由以下公式得出t时刻荧光信号强度：其中，τ为衰减周期；2)利用卷积法对峰值位置的变化进行分析：3)利用高斯拟合对散射光信号与荧光信号进行表示，得到的脉冲信号进行高斯拟合。优选地，所述步骤3)中对拟合结果的评价指标为：均方分误差：RMSE(root-mean-square-error)确定系数：R-square预测数据与原始数据均值之差的平方和SSR。原始数据和均值之差的平方和SST。应当理解，前述大体的描述和后续详尽的描述均为示例性说明和解释，并不应当用作对本专利技术所要求保护内容的限制。附图说明参考随附的附图，本专利技术更多的目的、功能和优点将通过本专利技术实施方式的如下描述得以阐明，其中：图1为散射光脉冲产生原理示意图；图2为散射光产生机理示意图；其中图2(a)为第一个微元S1的散射光fs_S1产生过程示意图；图2(b)为最后一个微元Sn的散射光fs_Sn产生过程示意图；图3为荧光产生机理示意图；图4为利用卷积法与时间微元法进行分析的结果图；图5为前向散射光信号fs及高斯拟合结果图；图6为荧光信号fl1及高斯拟合结果图；图7为荧光信号fl2及高斯拟合结果图；图8为荧光信号fl3及高斯拟合结果图；图9为三路荧光信号荧光时延直方图。具体实施方式通过参考示范性实施例，本专利技术的目的和功能以及用于实现这些目的和功能的方法将得以阐明。然而，本专利技术并不受限于以下所公开的示范性实施例；可以通过不同形式来对其加以实现。说明书的实质仅仅是帮助相关领域技术人员综合理解本专利技术的具体细节。在下文中，将参考附图描述本专利技术的实施例。在附图中，相同的附图标记代表相同或类似的部件，或者相同或类似的步骤。本专利技术提供了一种对荧光脉冲信号进行高斯拟合的时延估计方法，利用荧光脉冲信号与散射光脉冲信号之间的时延，实现对荧光寿命的测量。高斯拟合的时延估计方法具有：(1)硬件电路系统结构简单，易于实现；(2)避免了光强调制及高速处理模块之间的同步等问题；(3)完全由数字信号处理芯片完成，计算速度快，灵活性好；(4)计算所得荧光寿命的时域分辨率不受ADC采样频率限制等优势。在流式细胞仪检测过程中，散射光信号强度与光斑位置关系的产生原理如图1所示，将荧光微球分成n个径向等宽度的微元，宽度为dx，令微元距位置原点(l＝0)的距离为x(0≤x≤2R)，则该微元对应的表面积为第一个微元S1的散射光fs_S1产生过程如图2(a)所示，在时间微元范围内散射光强度不存在衰减，故用矩形进行表示。fs_S1起始时刻为0，终止时刻为L/v。最后一个微元Sn的散射光fs_Sn产生过程如图2(b)所示。fs_Sn起始时刻为2R/v，终止时刻为(L+2R)/v。光斑范围内各个微元的散射光叠加可得对应时刻的瞬时散射光强度。t时刻散射光强度为：荧光信号具有衰减特性，假定为单指数衰减其中τ为衰减周期(荧光信号衰减为初始值的1/e)。散射光信号fs(t)包含微球位置、相应的激发光斑强度分布等信息，故利用fs(t)作为荧光信号时间微元分析的初始时刻光强信号。图3为荧光信号按照时间微元分析得到的效果图。t时刻荧光信号强度为：另外，利用卷积法对峰值位置的变化进行分析：fl(t)与fs(t)的相频差即为的相频即，对应的时域时延为τ。散射光信号：荧光信号：本专利技术的一个实施例激发光斑短轴长度L＝22μm，沿短轴的强度分布为荧光微球半径R＝10μm，微球流速v＝10.5m/s。FS、FLm的采样频率为100MHz，即脉冲序列相邻数据点之间的时间间隔为10ns。本专利技术利用卷积法及时间微元法计算得到的荧光信号强度曲线如图4所示。时间微元法曲线的包络与卷积法结果一致，两者峰值点相对于fs(t)峰值点时延均为衰减周期τ。散射光信号具有中心对称性，且对称中心为光斑中心位置l＝L/2，对应时刻为t＝(L+2R)/2。荧光信号相对于散射光信号的时延为衰减函数的衰减周期τ，同时衰减函数会引起荧光信号波形的非对称性，并且τ越大波形的非对称性越明显。由于荧光信号宽度约为τ的102～103倍，故引起的非对称性不是很显著。本专利技术利用高斯拟合对散射光信号与荧光信号进行表示，对100MSPS得到的脉冲信号进行高斯拟合，拟合结果如表1所示。具体波形如图5～8所示。拟合结果的评价指标为：均方分误差：RMSE(root-mean-square-error)确定系数：R-squareSSR:Sumofsquaresoftheregression.预测数据与原始数据均值之差的平方和。SST:Totalsumofsquares.原始数据和均值之差的平方和。拟合结果为：a1＝10.51；b1＝2.01；c1＝0.5；评价指标：RMSE＝2.21×10-5；R-square＝1.000；即，包含单指数衰减特性的荧光信号可以用高斯模型进行有效的表征。本专利技术分别对100本文档来自技高网...

【技术保护点】
一种流式细胞系统荧光寿命高斯拟合测量方法，所述方法包括以下步骤：1)利用时间微元法求出瞬时散射光强度和荧光信号强度，由以下公式得出t时刻散射光强度：fs(t)=∫02R2πR2-(R-x)2e-(vt-x-L/2)22σ2dx---(1)]]> 其中，R为荧光微球的半径，dx为微元的宽度，x为微元距位置原点的距离； 由以下公式得出t时刻荧光信号强度：fl(t)=Σfs(t)δ(t-nτ)e-t-nττ---(2)]]> 其中，τ为衰减周期； 2)利用卷积法对峰值位置的变化进行分析：fl(t)=fs(t)*e-tτ---(3)]]> 3)利用高斯拟合对散射光信号与荧光信号进行表示，得到的脉冲信号进行高斯拟合。

【技术特征摘要】
1.一种流式细胞系统荧光寿命高斯拟合测量方法，所述方法包括以下步骤：1)利用时间微元法求出瞬时散射光强度和荧光信号强度，由以下公式得出t时刻散射光强度：fs(t)=∫02R2πR2-(R-x)2e-(vt-x-L/2)22σ2dx---(1)]]>其中，R为荧光微球的半径，dx为微元的宽度，x为微元距位置原点的距离；由以下公式得出t时刻荧光信号强度：fl(t)=Σfs(t)δ(t-nτ)e-t-nττ---(2)]]>其中，τ为衰减周期；2)利用卷积法对峰值位置的变化进行分析：fl(t)=fs(t)*e-tτ---(3)]]>3)利用高斯拟合对散射光信号与荧光信...

【专利技术属性】
技术研发人员：祝连庆，张文昌，娄小平，刘超，潘志康，董明利，孟晓辰，
申请(专利权)人：北京信息科技大学，
类型：发明
国别省市：北京;11