CN108693099B - 用于流式细胞仪的全自动荧光补偿方法 - Google Patents

用于流式细胞仪的全自动荧光补偿方法 Download PDF

Info

Publication number
CN108693099B
CN108693099B CN201810220169.6A CN201810220169A CN108693099B CN 108693099 B CN108693099 B CN 108693099B CN 201810220169 A CN201810220169 A CN 201810220169A CN 108693099 B CN108693099 B CN 108693099B
Authority
CN
China
Prior art keywords
fluorescence
compensation
matrix
detection channel
detection
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
Application number
CN201810220169.6A
Other languages
English (en)
Other versions
CN108693099A (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.)
Suzhou Institute of Biomedical Engineering and Technology of CAS
Original Assignee
Suzhou Institute of Biomedical Engineering and Technology of CAS
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 Suzhou Institute of Biomedical Engineering and Technology of CAS filed Critical Suzhou Institute of Biomedical Engineering and Technology of CAS
Priority to CN201810220169.6A priority Critical patent/CN108693099B/zh
Publication of CN108693099A publication Critical patent/CN108693099A/zh
Application granted granted Critical
Publication of CN108693099B publication Critical patent/CN108693099B/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/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434Optical arrangements

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明公开了一种用于流式细胞仪的全自动荧光补偿方法,包括以下步骤:步骤1)制备多色样品上机测试,得到原始荧光信号结果;步骤2)保持电压增益与步骤1)中的相同,各单色小球和空白小球按比例混合上机测试,得到各单色小球和空白小球的子集在各个检测通道的荧光信号原始输出值;步骤3)计算荧光泄漏矩阵K和自发荧光矩阵A;步骤4)计算荧光补偿矩阵KC;步骤5)进行全自动荧光补偿,得到荧光补偿后的检测结果。本发明利用混合小球,只需测试一次即可得到补偿结果,而不需要多次测试阴性对照样本和各单阳性样本,简化了操作,缩短了时间,消除了实验间的误差;本发明还公开了荧光泄漏矩阵K的具体计算方法,本发明方法简便、效果显著。

Description

用于流式细胞仪的全自动荧光补偿方法
技术领域
本发明涉及流式数据分析领域,特别涉及一种用于流式细胞仪的全自动荧光补偿方法。
背景技术
在进行流式细胞分析的时候,待测样本通常携带两种或两种以上的荧光素,如藻红蛋白(PE)、异硫氰荧光素(FITC)、叶丝素蛋白(PerCP)等。荧光素经过激光激发以后发射出不同波长的荧光。理论上来讲,通过选择合适的分光组件可以将不同的荧光分开,使得每个探测器只接收一种荧光信号,而不会检测到其他的荧光信号。然而实际上,目前常用的荧光染料的激发波长或发射波长都是正态或者偏态分布,有很宽的范围,因此尽管他们的发射峰值各不一样,但是发射谱常有重叠,如图1所示,为FITC、PE的发射波长,可以看出两者的波谱有重叠的现象,因而PE探测器能检测到由FITC发射的525nm左右波长的光信号。同样,FITC探测器也能检测到PE的信号,只不过每种荧光探测器检测到的荧光信号以一种荧光素为主。这样就直接影响了检测结果的准确性,因此,多色分析时必须进行光谱重叠的校正,即进行荧光补偿。
现有流式细胞仪一般采用手动补偿的方法来实现。手动补偿对于两色和三色分析可以很好地完成补偿功能,但是对于四色及四色以上的分析操作难度太大而导致几乎不可能实现。而流式细胞仪提供的自动补偿方法需要制备阴性对照样本以及所有单阳性对照样本,对样品需求量大,样品制备繁琐,样品检测次数多,导致实施难度高。并且由于多次检测间的误差,导致补偿效果往往并不好。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种用于流式细胞仪的全自动荧光补偿方法。
本发明提供的方法,不需要制备大量的荧光补偿对照样本,且一次上样就可以得到补偿结果,方法简便有效。
为解决上述技术问题,本发明采用的技术方案是:一种用于流式细胞仪的全自动荧光补偿方法,包括以下步骤:
步骤1)制备多色样品上机测试,得到理想的原始荧光信号结果;
步骤2)保持电压增益与步骤1)中的相同,各单色小球和空白小球按比例混合上机测试,得到各单色小球和空白小球的子集在各个检测通道的荧光信号原始输出值;
步骤3)根据步骤2)的检测结果,计算荧光泄漏矩阵K和自发荧光矩阵A;
步骤4)计算荧光补偿矩阵Kc
步骤5)根据补偿公式S,进行全自动荧光补偿,得到荧光补偿后的检测结果。
优选的是,在所述步骤3)中,根据步骤2)所得的检测结果,设门得到空白小球以及各单色小球的子集;
通过计算各检测通道之间的泄漏系数kij,从而得到自发荧光矩阵A和泄漏矩阵K:
A=[a1 a2 a3 … an]T
Figure BDA0001599444380000021
其中,i=1,…n;j=1,…n;n为总检测通道数;
kij为理论上只能在第i个检测通道检测到的荧光泄漏到第j个检测通道的比例;
Figure BDA0001599444380000022
为第i个检测通道对应的荧光单色小球子集在第j个检测通道的荧光输出平均值;aj为空白小球子集在第j个检测通道的自发荧光或者非特异性荧光信号平均值;
Figure BDA0001599444380000023
为第i个检测通道对应的荧光单色小球子集在第m个检测通道的荧光输出平均值;ak为空白小球子集在第k个检测通道的自发荧光或者非特异性荧光信号平均值。
优选的是,所述步骤4)中,荧光补偿矩阵KC为步骤3)中所得的泄漏矩阵K的转置矩阵的逆矩阵,即KC=(KT)-1
优选的是,所述步骤5)中,根据补偿公式S=Kc×[O-A],对多色样品分析结果进行全自动荧光补偿,并最终得到荧光补偿后的检测结果,即为真实检测结果。
优选的是,自发荧光矩阵A由所述步骤3)所得,荧光补偿矩阵KC由所述步骤4)中所得,补偿公式S=Kc×[O-A]中,
O=[o1 o2 o3 … on]T,oj为第j个检测通道的原始输出信号;
S=[s1 s2 s3 … sn]T,sj为第j个检测通道的理论输出信号,即最终的荧光补偿后的真实检测结果。
优选的是,所述步骤1)中还包括调整包含电压增益在内的参数至合适的范围,获得理想的原始荧光信号结果。
优选的是,在所述步骤2)中,各单色小球和空白小球按照等比例混合,以获得更好的统计结果。
其中,单色小球和空白小球按照其他任何比例混合也可以得到补偿矩阵。
其中,步骤3)、步骤4)以及步骤5)涉及大量的运算,一般由计算机完成。
其中,优先推荐上述的步骤顺序。但是也可以选择其他组合方式,其中步骤1)和步骤2)的顺序可以互换,但是要确保电压增益保持先执行的步骤中的数值。步骤3)、4)必须在步骤2)之后,且在步骤5)之前。
本发明的用于流式细胞仪的全自动荧光补偿方法,对于5色以内(包含5色)的荧光检测可以取得较好的效果。对于现有的流式细胞仪,单激光所包含的荧光通道数不会超过5色。而不同激光之间的荧光通道不会同时检测,因此不需要补偿。所以5色补偿能满足现有的所有商业流式细胞仪的需求。
本发明的有益效果是:本发明的用于流式细胞仪的全自动荧光补偿方法,利用混合小球,只需测试一次即可得到补偿结果,而不需要多次测试阴性对照样本和各单阳性样本,简化了操作,缩短了时间,消除了实验间的误差;本发明还公开了荧光泄漏矩阵K的具体计算方法,本发明增加了自发荧光,改善了补偿效果。本发明可以解决流式细胞仪多色分析过程中因为荧光光谱重叠造成的设门困难,统计不准确等问题,且不需要制备大量的荧光补偿对照样本,一次上样就可以得到补偿结果,方法简便、效果显著。
附图说明
图1为荧光光谱重叠的一个示例的示意图;
图2为本发明的用于流式细胞仪全自动荧光补偿方法的一种实施例的步骤示意图;
图3为本发明的一种实施例中的三色荧光补偿时混合小球的检测结果;
图4为本发明的一种实施例中的三色荧光检测时荧光泄漏的示意图;
图5为本发明的一种实施例中的荧光补偿前的示意图;
图6为本发明的一种实施例中的荧光补偿后的示意图;
图7为本发明的另一种实施例中的荧光补偿前的示意图;
图8为本发明的另一种实施例中的不增加自发荧光a时的荧光补偿后的示意图;
图9为本发明的另一种实施例中的增加自发荧光a时的荧光补偿后的示意图。
具体实施方式
下面结合实施例对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
应当理解,本文所使用的诸如“具有”、“包含”以及“包括”术语并不排除一个或多个其它元件或其组合的存在或添加。
以下提供本发明的本实施例的用于流式细胞仪的全自动荧光补偿方法的具体实施例,以FITC、PE、PE-Cy5三色实验为例,参照图2,其包括以下步骤:
步骤S10:制备FITC、PE、PE-Cy5三色实验样品上样分析。调整好电压增益,以获得较理想的实验结果。
步骤S11:保持步骤S10中的电压增益不变,FITC、PE、PE-Cy5单色小球和空白小球等比例混合上机测试,得到各单色小球和空白小球的子集在各个检测通道的荧光信号原始输出值。获得如图3所示的测试结果。其中也可以选用商业的包含FITC、PE、PE-Cy5单色小球和空白小球的混合小球。
步骤S12:计算荧光泄漏矩阵K和自发荧光矩阵A。
根据荧光泄漏的原理,三色实验的荧光泄漏模型如图4所示。
kij为理论上应该只能在第i个通道检测到的荧光泄漏到第j个通道的比例;si为第i个通道的原始输入值,也是我们要求的真实荧光信号;oi为第i个通道检测到的荧光值;ai为第i个通道检测到的自发荧光或者非特异性荧光信号。
由图4所示的模型可以得到:
Figure BDA0001599444380000051
令O=[o1 o2 o3]T,S=[s1 s2 s3]T,A=[a1 a2 a3]T
Figure BDA0001599444380000052
则上式可以记为:O=KT×S+A
变换可得:S=(KT)-1×[O-A],将(KT)-1记为Kc,则S=Kc×[O-A]。
具体地,根据步骤S11所得的检测结果,设门得到空白小球以及各单色小球的子集,如图3所示。计算各通道之间的泄漏系数kij,从而得到泄漏矩阵K。
假设
Figure BDA0001599444380000054
为第i个通道对应的荧光单色小球子集在第j个通道的荧光输出平均值;aj为空白小球子集在第j个通道的自发荧光或者非特异性荧光信号平均值,则对于通道1对应的单色小球子集,s2=0,s3=0,因此可得:
Figure BDA0001599444380000053
且:k11+k12+k13=1,因此可得:
Figure BDA0001599444380000061
Figure BDA0001599444380000062
Figure BDA0001599444380000063
其他同理可得,因此:
Figure BDA0001599444380000064
其中:
Figure BDA0001599444380000065
其中,i=1,…n;j=1,…n;n为总检测通道数;
kij为理论上只能在第i个检测通道检测到的荧光泄漏到第j个检测通道的比例;
Figure BDA0001599444380000067
为第i个检测通道对应的荧光单色小球子集在第j个检测通道的荧光输出平均值;aj为空白小球子集在第j个检测通道的自发荧光或者非特异性荧光信号平均值;
Figure BDA0001599444380000068
为第i个检测通道对应的荧光单色小球子集在第m个检测通道的荧光输出平均值;ak为空白小球子集在第k个检测通道的自发荧光或者非特异性荧光信号平均值。
例如对于图3所示的实验:
Figure BDA0001599444380000066
步骤S13:计算荧光补偿矩阵Kc,Kc=(KT)-1;。
步骤S14:进行全自动荧光补偿。
具体地,根据补偿公式S=Kc×[O-A]对多色分析结果进行全自动荧光补偿。
其中,Kc、A为步骤S13、S12中所得。
O=[o1 o2 o3]T,oj为第j个通道的原始输出信号。
S=[s1 s2 s3]T,sj为第j个通道的理论输出信号,也就是我们需要的补偿后的结果。
图1是荧光光谱重叠的一个示例的示意图,其示出了FITC(异硫氰荧光素)与PE(藻红蛋白)发射波长的叠加情况。
图5、图6为本实验荧光补偿前后的对比图。可以看出,本发明能够有效解决流式细胞仪荧光检测中产生的荧光重叠问题,完成荧光补偿。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节。

Claims (3)

1.一种用于流式细胞仪的全自动荧光补偿方法,其特征在于,包括以下步骤:
步骤1)制备多色样品上机测试,得到原始荧光信号结果;
步骤2)保持电压增益与步骤1)中的相同,各单色小球和空白小球按比例混合上机测试,得到各单色小球和空白小球的子集在各个检测通道的荧光信号原始输出值;
步骤3)根据步骤2)的检测结果,计算荧光泄漏矩阵K和自发荧光矩阵A;
步骤4)计算荧光补偿矩阵Kc
步骤5)根据补偿公式S,进行全自动荧光补偿,得到荧光补偿后的真实检测结果;
在所述步骤3)中,根据步骤2)所得的检测结果,设门得到空白小球以及各单色小球的子集;
通过计算各检测通道之间的泄漏系数kij,从而得到自发荧光矩阵A和泄漏矩阵K:
A=[a1 a2 a3 L an]T
Figure FDA0003033047270000021
其中,i=1,...n;j=1,...n;n为总检测通道数;
kij为理论上只能在第i个检测通道检测到的荧光泄漏到第j个检测通道的比例;
Figure FDA0003033047270000022
为第i个检测通道对应的荧光单色小球子集在第j个检测通道的荧光输出平均值;aj为空白小球子集在第j个检测通道的自发荧光或者非特异性荧光信号平均值;
Figure FDA0003033047270000023
为第i个检测通道对应的荧光单色小球子集在第m个检测通道的荧光输出平均值;ak为空白小球子集在第k个检测通道的自发荧光或者非特异性荧光信号平均值;
所述步骤4)中,荧光补偿矩阵KC为步骤3)中所得的泄漏矩阵K的转置矩阵的逆矩阵,即KC=(KT)-1
所述步骤5)中,根据补偿公式S=Kc×[O-A],对多色样品分析结果进行全自动荧光补偿,并最终得到荧光补偿后的荧光信号值,即为真实检测结果;
自发荧光矩阵A由所述步骤3)所得,荧光补偿矩阵KC由所述步骤4)中所得,补偿公式S=Kc×[O-A]中,
O=[o1 o2 o3 L on]T,oj为第j个检测通道的原始输出信号;
S=[s1 s2 s3 L sn]T,sj为第j个检测通道的理论输出信号,即最终的荧光补偿后的真实检测结果。
2.如权利要求1所述的用于流式细胞仪的全自动荧光补偿方法,其特征在于,所述步骤1)中还包括调整包含电压增益在内的参数至合适的范围,获得理想的原始荧光信号结果。
3.如权利要求1所述的用于流式细胞仪的全自动荧光补偿方法,其特征在于,在所述步骤2)中,各单色小球和空白小球按照等比例混合。
CN201810220169.6A 2018-03-16 2018-03-16 用于流式细胞仪的全自动荧光补偿方法 Active CN108693099B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810220169.6A CN108693099B (zh) 2018-03-16 2018-03-16 用于流式细胞仪的全自动荧光补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810220169.6A CN108693099B (zh) 2018-03-16 2018-03-16 用于流式细胞仪的全自动荧光补偿方法

Publications (2)

Publication Number Publication Date
CN108693099A CN108693099A (zh) 2018-10-23
CN108693099B true CN108693099B (zh) 2021-07-09

Family

ID=63844550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810220169.6A Active CN108693099B (zh) 2018-03-16 2018-03-16 用于流式细胞仪的全自动荧光补偿方法

Country Status (1)

Country Link
CN (1) CN108693099B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109580458B (zh) * 2018-11-30 2020-02-14 苏州深析智能科技有限公司 流式细胞智能免疫分型方法、装置以及电子设备
CN109490176B (zh) * 2018-11-30 2019-11-22 苏州深析智能科技有限公司 流式细胞检测数据处理方法、装置及分型方法、装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106644901A (zh) * 2015-10-14 2017-05-10 北京信息科技大学 一种用于流式细胞仪的荧光补偿方法

Also Published As

Publication number Publication date
CN108693099A (zh) 2018-10-23

Similar Documents

Publication Publication Date Title
CN105718744B (zh) 一种基于深度学习的代谢质谱筛查方法及系统
Kettman et al. Classification and properties of 64 multiplexed microsphere sets
CN108693099B (zh) 用于流式细胞仪的全自动荧光补偿方法
US10371639B2 (en) Detecting fluorescent material in a stained particle by comparison with an unstained particle over a plurality of frequency bands and by estimating a linear combination of base vectors
US5084394A (en) Method for corrective calibration of a flow cytometry using a mixture of fluorescent microbeads and cells
US20230349819A1 (en) Fluorescent spectrum correcting method and fluorescent spectrum measuring device
Roederer Compensation in flow cytometry
EP2551665B1 (en) Information processing apparatus, method and program for fitting a measured spectral distribution.
EP2187199B1 (en) Instrument setup method for a fluorescence analyzer
US7613583B2 (en) Method and system for particle data processing
Solly et al. Comparable flow cytometry data can be obtained with two types of instruments, Canto II, and Navios. A GEIL study
JP6707637B2 (ja) 塗料表面の質感パラメータを決定するための方法
US8779387B2 (en) Method and system for detecting fluorochromes in a flow cytometer
Boglione et al. Assessing signals of TMD physics in SIDIS azimuthal asymmetries and in the extraction of the Sivers function
KR100914952B1 (ko) 컴퓨터 조색장치를 이용한 컬러의 처방예측 방법
CN109030321B (zh) 用于流式细胞仪的数据处理方法
CN103954713A (zh) 一种染料纯度的快速检测方法
US5314824A (en) Method of setting up a flow cytometer
CN109402924A (zh) 一种提高染色加工准确性的方法
CN214427452U (zh) 便携式压力与精神疾病客观指标评估试纸条
WO2021046502A3 (en) Kits and methods for testing for lung cancer risks
JPH08510051A (ja) 保存対照細胞を蛍光および光散乱測定の校正において用いる方法
Baghla et al. Rapid quantitative analysis of fermentation broth samples to assess efficiency of engineered yeast strain turnover
WO2023153316A1 (ja) 信号処理方法、信号処理装置、及び信号処理システム
WO2018216181A1 (ja) 測定方法

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