多重定量PCR系统中多色荧光检测和光谱串扰校正方法.docx
- 文档编号:30578724
- 上传时间:2023-08-17
- 格式:DOCX
- 页数:7
- 大小:426.45KB
多重定量PCR系统中多色荧光检测和光谱串扰校正方法.docx
《多重定量PCR系统中多色荧光检测和光谱串扰校正方法.docx》由会员分享,可在线阅读,更多相关《多重定量PCR系统中多色荧光检测和光谱串扰校正方法.docx(7页珍藏版)》请在冰豆网上搜索。
多重定量PCR系统中多色荧光检测和光谱串扰校正方法
臧留琴1 张镇西1 苗宝刚2 彭年才2* 李 政2
( )
1西安交通大学生命学院生物医学信息工程教育部重点实验室,陕西西安710049
2西安交通大学机械学院机械制造系统工程国家重点实验室,陕西西安710049
摘要多色荧光的激发检测光路是多重定量聚合酶链式反应(PCR)系统的核心组成部分。
根据系统对荧光激发的均一性、荧光检测的时间和灵敏度需求,提出一种基于磁光开关和光电倍增管(PMT)的多色荧光激发检测光路。
通过四色LED单独激发和PMT检测来提高灵敏度;利用磁光效应实现光路的电控切换,配合一维扫描机构和滤光片切换装置,完成96孔标准PCR板的四色荧光扫描。
该设计避免了多色荧光光路的系统串扰。
由于荧光染料本身光谱特性产生的荧光光谱串扰是通过标准迭代的四维聚类分析算法计算串扰矩阵,并通过4种常用染料的实验研究,对所建荧光检测系统进行荧光光谱串扰评估。
关键词光谱学;多重定量;多色荧光检测;光谱串扰校正;磁光开关
中图分类号O439 文献标识码A doi:
10.3788/AOS201434.0117002
MulticolorFluorescenceDetectionintheMultiplexQuantitativePCRSystemandSpectraCrosstalkCorrectionMethod
ZangLiuqin1ZhangZhenxi1 MiaoBaogang2PengNiancai2 LiZheng2
1KeyLaboratoryofBiomedicalInformationEngineering,MinistryofEducation,SchoolofLifeScience
烄 烌
andTechnology,Xi′anJiaotongUniversity,Xi′an,Shaanxi710049,China
2StateKeyLaboratoryforManufacturingSystemEngineering,SchoolofMechanicalEngineering,
烆 Xi′anJiaotongUniversity,Xi′an,Shaanxi710049,China 烎
AbstractOpticaldesignformulticolorfluorescencedetectionisthecorepartofmultiplexquantitativepolymerasechainreaction(PCR)system.Accordingtothesystemrequirementsforexcitationuniform,detectiontimeandsensitivity,anopticalsystemispresented,whichisbasedonmagneto-opticalswitchandphotomultipliertube(PMT)detector.FourLEDsexcitationunitandPMTdetectionunitareadoptedtoimprovesystemsensitivity.Electronicopticalswitchcombinedwiththeone-dimensionalscanningmechanismandfiltertoimplementfour-colorfluorescencescanningona96-wellstandardPCRplate.Thedesignavoidsthecrosstalkamongthefourcoloropticalpaths.Whereasthecrosstalkgeneratedbythespectralcharacteristicsofthefluorescentdyeitselfcanbeestimatedbystandarditerativefour-dimensionalclusteringanalysisalgorithm.Wefinallyevaluatethesystematicalfluorescentspectracrosstalkbytheexperimentalstudyoffourcommonlyuseddyes.
Keywordsspectroscopy;multiplexquantitative;multicolorfluorescencedetection;spectracrosstalkcorrection;
magneto-opticalswitch
OCIScodes170.3880;170.6280;120.1880;120.4570
收稿日期:
2013-07-12;收到修改稿日期:
2013-09-13
基金项目:
国家重大科学仪器设备开发专项(2012YQ030261)、国家自然科学基金(61120106013)、国家科技支撑计划
(2012BAI19B057)
作者简介:
臧留琴(1982—),女,博士研究生,主要从事生命科学检测仪器及诊断技术等方面的研究。
E-mail:
lqzang@stu.xjtu.edu.cn
导师简介:
张镇西(1951—),男,博士,教授,博士生导师,主要从事生物医学光子学影像与光谱分析技术、纳米光生物物理与生物光量子调控技术等方面的研究。
E-mail:
zxzhang@mail.xjtu.edu.cn
*通信联系人。
E-mail:
pnc@mail.xjtu.edu.cn
1 引 言
光 学 学 报
代数法将测量的荧光值转化为实际的荧光值。
设计
实时荧光定量聚合酶链式反应(PCR)技术由于具有实时、灵敏、特异等优点,已经成为分子研究的主要手段之一,广泛应用于分子生物学和医学研究及临床。
但是普通的实时定量PCR受荧光检测通道限制,单次反应只能分析单个靶序列,使其在科学研究中的使用受到限制。
基于多色荧光检测和分析技术的多重定量PCR方法,通过在一个反应体系内标记不同荧光染料,可实现一次反应检测和多个靶序列鉴别[1]。
基于多色荧光检测的多重PCR方法中,关键技术是多色荧光激发和检测系统的实现,而系统对激发的均一性、检测的快速性和灵敏度有很高的要求。
常规的多色荧光激发和检测系统主要有单孔扫描和整板成像两大类,其中,单孔扫描方式一般以步进电机带动扫描镜头进行二维运动以实现所有荧光信号的采集,这种方式可以获取特异性较强的荧光信号,但是步进电机控制的精确性及机械运动噪声会影响荧光信号的精度[2-3];而整板成像是基于电荷耦合器件(CCD)成像的荧光采集技术,其优点在于可以同时采集所有样品的荧光信号,但其灵敏度低,荧光信号之间有较大干扰,存在明显的“边缘效应”而造成较大的成像误差[4]。
本文提出一种多色荧光激发
和检测装置,并在该装置下对实际获取的实时PCR荧光数据进行分析、处理及验证。
的多重PCR采用4个波长段的荧光报告基团,理论
上在一次实验中,可以用4种染料作为探针的荧光报告基团,结合4种不同的靶基因,进行4重PCR反应。
2实验方法
2.1多色荧光激发检测系统的设计
采用4种波长冷光源发光二极管(LED)作为激发光源,组合相应的带通滤光片及聚焦透镜,形成激发单元;相应的4种检测波长的带通滤光和聚焦透镜,组合光电倍增管(PMT)形成检测单元;激发光通过光纤束从激发单元传输至样品池,而发射荧光也通过光纤束从样品池传输至检测单元;激发光纤束和发射光纤束利用基于磁光效应的光开关[5]实现光路在不同检测孔间的切换;步进电机带动激发和检测单元做一维运动进行分色扫描,另一个电机带动光纤束做一维扫描,完成96孔的扫描检测,光路示意图如图1所示。
系统中设计的激发和荧光发射波长如表1所示,使用的荧光染料分别为羧基荧光素(FAM),亚磷酰胺(HEX),德克萨斯红(TexasRed)和菁染料琥珀亚酰胺酯(Cy5)。
表1系统所采用激发和荧光发射波长
Table1Excitationandemissionwavelengthinthesystem
Channel Excitation
Emission
Dye
针对所设计的光学检测系统,需要解决其由于
荧光染料光谱重叠引起的色彩补偿问题。
要分析每一个检测通道获得的荧光数据首先必须得到纯荧光染料在每个检测通道的校准数据和自动荧光控制数据,利用这些数据计算出信号交叉常数,并通过矩阵
wavelength/nmwavelength/nm
1 470 525 FAMTM
2 523 564 HEXTM
3 571 612 Texasred
4 630 670 Cy5
图1光学检测系统
Fig.1Opticaldetectionsystem
臧留琴等:
多重定量PCR系统中多色荧光检测和光谱串扰校正方法
荧光激发和检测光纤束采用一种磁光效应,利用法拉第效应,通过改变所加载的磁场改变光的传输通道[6],如图2所示。
在4个检测通道的荧光值,计算串扰影响矩阵M[8]。
确定荧光串扰矩阵,必须明确各个荧光染料在各检测通道的荧光串扰系数k(i,j),由该系数形
成多通道的荧光串扰矩阵M。
荧光染料在检测通
道的荧光串扰系数是该染料在该检测通道读取荧光强度去除荧光染料自体背景荧光值后,与该染料在所有通道中读取的荧光强度值与背景值之差的比值,可以表示为
N
k(i,j)= o(i,j)-a(j)
, (1)
图2基于磁光效应的光通路切换
Fig.2Lightpathswitchbasedonmagneto-opticeffect
激发光检测光纤束示意图如图3所示,由一个
[o(i,n)-a(n)]
∑
n=1
式中k(i,j)是染料i在通道j中的串扰常数,o(i,j)是检测到的染料i在通道j中的光强,a(j)是通道j中的自体荧光(即背景荧光),n是通道序号,N是通
道总数。
荧光串扰矩阵可表示为
熿k11k12k13k14燄
输入光纤,8个输出光纤及光路切换模块组成,其中光路切换模块有5个引脚对外,实现光路切换模块
k21k22k23k24
M=
k31k32k33k34
. (2)
的供电和控制。
图3检测光纤束
燀k41k42k43k44燅
大多数荧光染料和荧光染料之间的串扰都是随温度变化的[9]。
校准时以0.2℃/s的升温速度从40℃~95℃进行连续荧光检测。
荧光强度随温度的变化曲线可以用3次多项式来拟合,保存最佳拟
Fig.3Opticalfiberbundlefordetection
通过控制端1、2、3的电平切换,实现光通路切换,所选用的光通路切换装置具体参数如表2所示。
表2光通路切换装置参数
Table2Parametersofopticalpathswitchingdevice
Parameter Defaultvalue
合参数用于之后的插值。
为了实现荧光数据的颜色补偿,每次采集的温度用来对校准曲线温度敏感的荧光值进行插值。
荧光强度还要考虑校准时和测量时系统的电子增益。
o(i,j)=[w(i,j)T3+x(i,j)T2+y(i,j)T+z(i,j)][GD(j)/GC(j)], (3)
Insertionloss/dB 0.9~1.2
a(j)=[m(j)T3
+n(j)T2
+p(j)T+q(j)]×
Returnloss/dB ≥40
Isolation/dB ≥35
Extinctionratio/dB 0.15~0.3
[GD(j)/GC(j)], (4)
式中w(i,j),x(i,j),y(i,j),z(i,j)是荧光随温度变化的三次多项式系数,m(i,j),n(i,j),p(i,j),q(i,j)是通道自体荧光随温度变化的三
Switchingtime/μs 10~400
次多项式系数。
GD(j)和
GC(j)分别是测量时和校
2.2荧光串扰的校正方法
在基于4色荧光检测的PCR系统中,测得的原始荧光数据包含4个波长的荧光浓度[7],组成“检测
空间”向量s,而实际需要得到的模版浓度由荧光染料浓度表征,形成“染料空间”向量f,此处需要通过一个4×4的转换矩阵M实现检测空间到染料空间的转换,此转换矩阵即为荧光串扰矩阵,s=Mf。
M的一列表示某一种染料发射的荧光在4种检测波长
下的相对浓度值。
通过测定多个浓度的某种纯染料
准时系统的电子增益系数。
测量时每种染料的实际荧光信号强度f为
f=M-1[O-A], (5)式中f是每种染料的实际荧光信号强度,M是各种染料对不同通道的串扰矩阵,O是每个通道测量得到的荧光强度,A是荧光染料的背景荧光强度。
计算荧光串扰矩阵的最简单的方法是:
首先从原始数据中移除基线,其次需要选择4个与各荧光染料相关的典型峰值,然后计算出这4个峰值的相
0117002-3
光 学 学 报
对荧光强度。
这种方法虽然简便,但是由于只有4个数据,精确性较差,一般在实际应用中不会采用。
熿1k12k22
k13k33
k14燄k44
以这种方法为基础,对其进行优化提出一种通过多次扫描得到多组数据形成多个矩阵,然后对多个矩
k21k11
M=
1k23
k33
k24k44
. (6)
阵进行平均的方法。
这种方法是在每次扫描时,每种染料在每个通道只取最大值,而其他值都设为0,然后利用最小二乘法对串扰矩阵进行优化。
一种通过荧光检测数据直接计算荧光串扰矩阵
的算法核心为四维聚类分析[10],包括数据搜集、基
k31
k11
k41
燀k11
具体算法如下:
k32
k22
k42k22
1k34
k44
k43 1
k33 燅
线校准、峰值识别以及自动归一化。
估计出的串扰矩阵的质量可以通过对训练数组中转换峰值点的三维散点分析进行评估。
采用一系列标准迭代聚类技术在模式识别中进行特征选择来确定串扰矩阵。
这种方法所需的处理步骤较少,但是需要移除基线且需要一组初始值来进行迭代。
对该方法可通过观察色彩分离后的数据轨迹进行评估。
在此基础上,提出了另一种对串扰矩阵进行估计的方法[11],其串扰矩阵每一列数据可通过一种特定的实验进行计算,该实验方法是在每次荧光扫描和读值时,只采用单一的荧光标记物,而数据采集则在4个通道同时进行,以此评估出每一种染料在各个通道的荧光串扰。
该方法利用迭代算法和奇异值分解算法使两次测量差异、一维散度和均方误差和达到最小值。
计算荧光串扰矩阵,需要考虑以下因素。
1)必须确定所读取的荧光值的准确性,以此作为评价各种串扰矩阵估计质量的基础,并优选出最佳矩阵估计算法。
2)在估计串扰矩阵时,需在训练过程(迭代)中采用尽量多的扫描数据组。
3)如何确定采集的荧光值的基线是一个关键问题,此处提出以下几种方法:
寻求一种不需要预先确定基线轨迹的算法;然后用预先给定的串扰矩阵来近似估计基线并对荧光数据进行校正。
4)峰值的确定通常可以通过低通滤波对数据进行多点平滑处理。
提出一种实验方法确定4种荧光染料在4个检测通道的荧光值,通过确定不同荧光染料在不同通道的相互关系,确定各染料在各通道的干扰,以此确定荧光串扰矩阵。
由于荧光定量PCR的检测值为相对荧光强度,在荧光串扰矩阵中,可以认为每一种荧光染料在其发射荧光波长的检测通道下,其影响因子为1,这样,实际荧光串扰矩阵只需要确定12
个元素。
因此串扰矩阵M可以改写为如下表达式:
1)初始化:
设i=1,用原始数据作为训练数据
组,初始估计矩阵M(0)是一个4×4单位矩阵。
同时需要设一个阈值α和一个正整数T,其中α是一个小正数。
2)抽样:
用不同染料在荧光采集通道获取荧光值,每两种染料的相对荧光值组成多个二维散点图。
以第一组分为X坐标,第二组分为Y坐标组成的二
维散点图为例。
确定边界:
选择第一组分的两个分位点作为范围的边界,在本次迭代中,当前训练数据组中处于边界范围内的这些点可用于估计出串扰矩阵中的元素
k21/k11。
区域分割:
将两边界所包含的范围平均分为多个相同的宽度。
选择极值点:
找到所有扫描数据点中,第一组分值落在一个给定分割的区域内而其第二组分值最小的散点。
3)一维线性回归:
对从最后一个步骤中获得的扫描数据进行第二组分对第一组分的一维线性回归计算,将此处估计出的斜率值作为k21/k11的下一个
估计值。
同理,该图中估计出的另一臂的斜率可作
为k12/k22的下一个估计值。
4)对其他参数进行估计:
根据上述步骤2和3,可估计出矩阵中其他各元素的值,并估计出串扰矩阵M。
5)评估颜色校正的质量:
计算步骤2、3和4中
得到的12个估计斜率的绝对值的最大值,如果该值低于阈值α,则进行步骤7,否则即按步骤6继续计算。
6)修正:
用M的逆矩阵以逐点相乘的方法对训练数据组进行修正,得到新的训练数据组。
设M(i)=M(i-1)×M,其中,i为迭代次数,每次增加1,
如果i>T,则进行步骤7,否则就从第一步骤重新计
算。
0117002-4
臧留琴等:
多重定量PCR系统中多色荧光检测和光谱串扰校正方法
7)归一化:
对M(i)进行归一化处理,使矩阵每列的数据相加为1。
该算法的输出结果为估算出的串扰矩阵。
首先,步骤2中选择的区域的分位点一般为训练数据组中每个组分值的60%~99%,选取99%的分位点来消除异常值。
其次,分割的区域不能太小也不能
太大,太小容易使迭代效率下降,且容易导致更大的偏差,而太大则会引起许多分割区域内不存在一个边界点,一般确定分割区域的数量时要求每个分割区域内平均包含6~10个有效数据点。
第三,数值实验结果表明一维线性回归要优于最小二乘回归法。
第四,结束迭代的基本原则是检查斜率是否平
行于X轴,此处通过为斜率设定阈值实现该原则,一旦估计的斜率值小于设定的阈值,迭代则终止。
阈值的选择依赖于精度要求和噪声水平。
根据多组数据分析和判定,一般选取0.05作为阈值。
最后,还对估计出的串扰矩阵进行了归一化处理,使矩阵的每一列相加之和为1。
2.3估算荧光串扰校正矩阵的实验方法
为了评估4色PCR光学系统的荧光串扰矩阵,设计以下试验进行原始数据采集。
4个检测通道的滤光片中心波长及带宽与4种染料发射谱相对关系如图4所示。
图4检测波长及带宽与染料发射谱关系
Fig.4Relationshipbetweendyeemissionspectrumanddetectionwavelengthandbandwidth
确定系统的荧光串扰矩阵,需首先对系统进行背景校正,背景荧光是PCR仪的环境荧光强度,主要由背景电信号、样本块中的污染物、塑料耗材中的水以及塑料耗材自身引起,背景校正就是通过实验获取系统的背景荧光值,生成校正文件用于后续实验数据处理消除背景荧光影响。
然后确定此光学系统的荧光染料浓度线性范围,最后进行纯染料校正,通过对获取的数据进行分析处理得到系统的荧光串扰矩阵。
背景校正方法:
在反应管中加入不含模版的PCR反应体系,制成背景校正板,测量其空白荧光强度。
在运行校正程序期间,定量PCR系统在10min内连续读取背景校正板的荧光强度,信号收集的温度为60℃(一般PCR反应中进行荧光采集的温度)。
随后,用软件计算所收集到的荧光强度的平均值,保存到校正文件中。
软件在今后的分析中将自动调用此校正文件,从实验数据中扣除背景信号。
梯度稀释液25μL;在D3~D12孔中各加入HEX的8种梯度稀释液25μL;在E3~E12孔中各加入TexasRed的8种梯度稀释液25μL;在F3~F12孔中各加入Cy5的8种梯度稀释液25μL。
将反应
板放入检测仪器中进行检测,在10min内连续读取的反应板的荧光强度,对读取的荧光数据进行分析,确定荧光染料的浓度线性范围。
纯染料校正方法:
将4种染料在其浓度线性范围内各稀释成24个浓度梯度,并按浓度大小顺序加入到96孔标准反应板中,制成纯染料校正板。
将纯染料校正反应板放置入仪器,进行光学扫描和荧光数据采集。
设置仪器以0.2℃/s的升温速率从40℃升到95℃,从40℃开始,每升高1℃,保温10s,进行一次荧光数据采集,记录采集的荧光数值用于后续计算。
3 结果与讨论
3.1光学系统性能
随机选择经定量PCR扩增的DNA试剂标本
确定系统荧光染料浓度线性范围方法:
在96孔
装入96
孔载台,设置激发端和检测端的磁光开关同
标准PCR板的C3~C12孔中各加入FAM的8种
步切换频率为
2500Hz,相当于每个通道的切换时
0117002-5
光 学 学 报
图5荧光检测速度测试图
Fig.5Speedtestinfluorecencedetection
本对应的输出幅度各异的脉冲输出上升沿和下降沿很陡,能满足之后模数转换及数字和软件系统对各路信号的识别和采集,实现多路荧光检测通道的快速扫描检测。
实测的单通道切换加上荧光采集时间小于500μs,8通道的荧光采集不到4ms。
3.2四色PCR光学系统荧光串
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多重 定量 PCR 系统 多色 荧光 检测 光谱 校正 方法
![提示](https://static.bdocx.com/images/bang_tan.gif)