多重定量PCR系统中多色荧光检测和光谱串扰校正方法
期刊论文，臧留琴等，《光学学报》2014年第1期（Vol.34, No.1），doi: 10.3788/AOS201434.0117002
基本信息
- 作者：臧留琴¹，张镇西¹，苗宝刚²，彭年才²*，李政²
- 单位：
- ¹西安交通大学生命学院生物医学信息工程教育部重点实验室
- ²西安交通大学机械学院机械制造系统工程国家重点实验室
- 基金：国家重大科学仪器设备开发专项（2012YQ030261）、国家自然科学基金（61120106013）、国家科技支撑计划（2012BAI19B057）
- 分类号：O439
- OCIS代码：170.3880；170.6280；120.1880；120.4570
研究背景与问题
实时荧光定量PCR技术因实时、灵敏、特异等优点，已成为分子生物学和医学研究的主要手段。但普通定量PCR受荧光检测通道限制，单次反应只能分析单个靶序列。多重定量PCR通过在一个反应体系内标记不同荧光染料，可实现一次反应检测多个靶序列。
核心技术挑战：多色荧光激发和检测系统需要满足三个关键需求：
- 激发的均一性
- 检测的快速性
- 检测的灵敏度
常规方案对比：
| 方案 | 原理 | 优点 | 缺点 |
|---|---|---|---|
| 单孔扫描 | 步进电机带动扫描镜头二维运动 | 特异性强 | 步进电机精度和机械噪声影响荧光信号精度 |
| 整板成像 | CCD成像荧光采集 | 可同时采集所有样品 | 灵敏度低，信号干扰大，“边缘效应”造成成像误差 |
系统设计：基于磁光开关+PMT的多色荧光检测光路
整体架构
采用4种波长LED作为激发光源，组合带通滤光片和聚焦透镜形成激发单元；4种检测波长的带通滤光片+聚焦透镜+**PMT（光电倍增管）**形成检测单元。
光路传输：激发光通过光纤束从激发单元传输至样品池，发射荧光也通过光纤束从样品池传输至检测单元。
磁光效应光开关：基于法拉第效应，通过改变加载的磁场改变光的传输通道，实现光路在不同检测孔间的切换。
扫描机制：
- 步进电机带动激发和检测单元做一维运动进行分色扫描
- 另一个电机带动光纤束做一维扫描
- 完成96孔标准PCR板的四色荧光扫描
四色荧光通道配置
| 通道 | 激发波长(nm) | 发射波长(nm) | 染料 |
|---|---|---|---|
| 1（Ch1） | 470 | 525 | FAM™ |
| 2（Ch2） | 523 | 564 | HEX™ |
| 3（Ch3） | 571 | 612 | Texas Red |
| 4（Ch4） | 630 | 670 | Cy5 |
磁光开关光纤束结构
- 1个输入光纤 + 8个输出光纤 + 光路切换模块
- 5个引脚对外（供电+控制）
- 通过控制端1、2、3的电平切换实现光通路切换
磁光开关关键参数：
| 参数 | 数值 |
|---|---|
| 插入损耗(Insertion Loss) | 0.9～1.2 dB |
| 回波损耗(Return Loss) | ≥40 dB |
| 隔离度(Isolation) | ≥35 dB |
| 消光比(Extinction Ratio) | 0.15～0.3 dB |
| 切换时间(Switching Time) | 10～400 μs |
系统光学性能实测
- 磁光开关同步切换频率：2500 Hz
- 每个通道切换时间：约 400 μs
- 单通道切换+荧光采集时间：< 500 μs
- 8通道荧光采集总时间：< 4 ms
- 染料浓度线性范围：10⁻¹⁰ ～ 10⁻³ mol/L（7个数量级）
- PMT输出脉冲上升沿/下降沿陡峭，满足模数转换和数字信号识别需求
设计优势：磁光效应组合光纤传输可以有效避免多色荧光光路的系统串扰（光路层面的串扰）。
荧光串扰校正：四维聚类分析算法
荧光染料本身的光谱特性产生的光谱串扰无法通过光学设计完全消除，需要通过算法校正。
核心概念
- 检测空间（s）：测量得到的4个通道的原始荧光值（4维向量）
- 染料空间（f）：实际的荧光染料浓度（4维向量）
- 串扰矩阵（M）：4×4转换矩阵，实现检测空间到染料空间的映射：s = M·f
- 校正过程：f = M⁻¹·(O - A)，其中O是测量的荧光强度，A是背景荧光强度
串扰系数计算
串扰系数 k(i,j) = 染料i在通道j中的荧光强度（扣除自体背景后） / 该染料在所有通道中荧光强度与背景差值的总和
k(i,j) = [o(i,j) - a(j)] / Σ[o(i,n) - a(n)] (n从1到N)
其中：
- o(i,j)：染料i在通道j中检测到的光强
- a(j)：通道j中的自体荧光（背景荧光）
- N：通道总数
4色系统的串扰矩阵M为4×4矩阵，每列表示某一种染料发射的荧光在4种检测波长下的相对浓度。
温度敏感性校正
荧光串扰随温度变化。校准时以0.2℃/s从40℃到95℃连续检测，荧光强度随温度变化用三次多项式拟合，保存最佳拟合参数。
测量时用采集的温度对校准曲线进行插值，同时考虑校准时和测量时系统的电子增益差异：
o(i,j) = [w·T³ + x·T² + y·T + z] × [GD(j) / GC(j)]
a(j) = [m·T³ + n·T² + p·T + q] × [GD(j) / GC(j)]
其中w,x,y,z和m,n,p,q分别是荧光和背景荧光随温度变化的三次多项式系数，GD(j)和GC(j)分别是测量时和校准时的电子增益系数。
标准迭代四维聚类分析算法
算法步骤：
-
初始化：i=1，用原始数据作为训练数据组，初始估计矩阵M⁽⁰⁾为4×4单位矩阵。设定阈值α（小正数，一般取0.05）和最大迭代次数T。
-
抽样+边界确定+区域分割+极值点选择：
- 用不同染料在荧光采集通道获取荧光值，每两种染料的相对荧光值组成二维散点图
- 选择第一组分的60%～99%分位点作为范围边界（消除异常值）
- 将边界范围平均分为多个相同宽度的分割区域（每个区域平均含6～10个有效数据点）
- 在每个分割区域内找到第二组分值最小的极值点
-
一维线性回归：对极值点进行第二组分对第一组分的一维线性回归，估计斜率作为 k21/k11 的下一个估计值。同理估计 k12/k22 等。
-
估计完整矩阵：根据步骤2和3估计矩阵中其他各元素的值，得到串扰矩阵M。
-
质量评估：计算12个估计斜率绝对值的最大值，若低于阈值α则转步骤7，否则继续。
-
修正：用M的逆矩阵对训练数据组逐点修正，得到新的训练数据组。M⁽ⁱ⁾ = M⁽ⁱ⁻¹⁾ × M。若i>T则转步骤7，否则回到步骤1重新计算。
-
归一化：对M⁽ⁱ⁾进行归一化处理，使矩阵每列数据相加为1。
算法关键设计点：
- 分位点取60%～99%以消除异常值
- 分割区域大小需适中：太小降低迭代效率且增大偏差，太大导致区域内无边界点
- 一维线性回归优于最小二乘回归法
- 迭代终止条件：斜率是否平行于X轴（斜率绝对值<阈值）
- 一般选取阈值α = 0.05
- 最后进行列归一化（每列和为1）
实验方法：纯染料校正板
三步校准流程：
- 背景校正：不含模板的PCR反应体系，测量空白荧光强度（10分钟内连续读取，取60℃平均值，保存为校正文件）
- 线性范围确定：4种染料各8种梯度稀释液，加入96孔板（C3
C12=FAM，D3D12=HEX，E3E12=Texas Red，F3F12=Cy5），10分钟内连续读取 - 纯染料校正：4种染料各24个浓度梯度（在线性范围内），加入96孔板制成纯染料校正板。0.2℃/s从40℃升到95℃，每升高1℃保温10s采集一次荧光数据
一次校正试验可获取55次采集×4通道×96孔的数据。
实验结果
最优串扰矩阵（未归一化）
| FAM | HEX | Texas Red | Cy5 | |
|---|---|---|---|---|
| Ch1 | 1 | 0.171429 | 0.001429 | 0.000086 |
| Ch2 | 0.457143 | 1 | 0.071429 | 0.001428 |
| Ch3 | 0.157143 | 0.025714 | 1 | 0.002861 |
| Ch4 | 0.328571 | — | — | 1 |
注：对角线为1（自身通道），表示以主通道为基准的相对串扰比例
归一化串扰矩阵（每列和为1）
| FAM | HEX | Texas Red | Cy5 | |
|---|---|---|---|---|
| Ch1 | 0.852 | 0.313 | 0.058 | 0.004 |
| Ch2 | 0.146 | 0.684 | 0.128 | 0.019 |
| Ch3 | 0.0012 | 0.0020 | 0.812 | 0.242 |
| Ch4 | 0.0008 | 0.0010 | 0.002 | 0.735 |
串扰特点分析：
- FAM → Ch2 串扰最严重（约14.6%），其次是 FAM → Ch1（主通道85.2%）
- HEX → Ch1 串扰约31.3%（这个数字看起来偏高，可能是行/列方向理解差异）
- Cy5串扰最小，因为波长最长，短波长通道基本不受Cy5荧光影响
- 通道之间的串扰是非对称的（短波长染料向长波长通道串扰更明显，符合荧光光谱特性）
系统个体差异
每个荧光检测系统存在个体差异（背景荧光值、相对荧光强度不同），因此每个新建系统都应进行系统的背景荧光检测及荧光串扰矩阵估算，以确保临床应用时的一致性。
结论
- 设计了基于四色LED激发 + PMT检测 + 磁光开关光路切换 + 一维扫描的多重定量PCR检测系统
- 磁光效应组合光纤传输可有效避免光通路之间的系统串扰
- 对于荧光染料本身光谱特性产生的串扰，通过标准迭代四维聚类分析算法计算串扰矩阵
- 针对4种常用染料（FAM/HEX/Texas Red/Cy5）完成实验验证和荧光串扰评估
- 实际应用时存储串扰矩阵的逆矩阵M⁻¹，逐点乘以采集的荧光数据即可得到校正后的真实染料浓度
关键术语
- 多重定量PCR：multiplex quantitative PCR，单管多靶标检测
- 磁光开关：magneto-optical switch，基于法拉第效应的光路切换器件
- PMT：光电倍增管，高灵敏度光探测器
- 串扰矩阵：cross-talk matrix / spillover matrix
- 检测空间 vs 染料空间：detection space vs dye space
- 四维聚类分析：four-dimensional clustering analysis
- 颜色补偿/荧光补偿：color compensation / fluorescence compensation