CN117995279A - 不同pcr仪器间定量或定性检测结果同质化换算方法 - Google Patents
不同pcr仪器间定量或定性检测结果同质化换算方法 Download PDFInfo
- Publication number
- CN117995279A CN117995279A CN202410397223.XA CN202410397223A CN117995279A CN 117995279 A CN117995279 A CN 117995279A CN 202410397223 A CN202410397223 A CN 202410397223A CN 117995279 A CN117995279 A CN 117995279A
- Authority
- CN
- China
- Prior art keywords
- value
- instrument
- standard curve
- test
- line segment
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 58
- 238000000265 homogenisation Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 117
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000003321 amplification Effects 0.000 claims description 35
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 35
- 230000008569 process Effects 0.000 claims description 4
- 239000012895 dilution Substances 0.000 claims description 3
- 238000010790 dilution Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 239000000523 sample Substances 0.000 description 57
- 238000003752 polymerase chain reaction Methods 0.000 description 46
- 238000004422 calculation algorithm Methods 0.000 description 12
- 108020004414 DNA Proteins 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000007405 data analysis Methods 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 108091028043 Nucleic acid sequence Proteins 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000011282 treatment Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000010222 PCR analysis Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000011088 calibration curve Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003753 real-time PCR Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000011425 standardization method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6851—Quantitative amplification
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/20—Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/10—Signal processing, e.g. from mass spectrometry [MS] or from PCR
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Organic Chemistry (AREA)
- Biotechnology (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Genetics & Genomics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Theoretical Computer Science (AREA)
- Zoology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Wood Science & Technology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Software Systems (AREA)
- Analytical Chemistry (AREA)
- Data Mining & Analysis (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioethics (AREA)
- Signal Processing (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
不同PCR仪器间定量或定性检测结果同质化换算方法，涉及医学检测技术领域，对于定量检测，利用仪器标准曲线及等比例计算方法进行换算。对于定性检测，步骤为：通过选择N个样本分别在参比仪器和测试仪器检测，获得N‑1对一一对应子线段，建立参比仪器标准曲线和测试仪器标准曲线，该标准曲线各为一根横线，横坐标为循环次数Cq值，纵坐标为仪器信号值；再对待测样本，获得其在测试仪器上荧光信号达到阈值线时对应Cq值的实际测试值，找到该实际测试值处于横线上的具体子直线段；再利用等比例计算方法计算出实际测试值对应到参比仪器上的同质化结果值。本发明可使不同仪器检测结果具有一致性和可比性，实现PCR检测结果互认。
Description
技术领域
本发明涉及医学检测技术领域，特别涉及不同PCR仪器间定量或定性检测结果同质化换算方法。
背景技术
聚合酶链反应（PCR）是一种在分子生物学和医学研究中广泛应用的技术。它将特定基因或序列的拷贝数通过测量PCR过程中DNA模板产生的荧光信号进行量化。PCR具有灵敏度高、动态范围广等优点，因此被广泛应用于许多领域。PCR的检测原理是根据每个PCR循环中荧光的增加情况进行实时监测。扩增反应的定量循环Cq值（又称Ct值）定义为荧光达到定量阈值所需循环的次数。目标复制越多，达到定量阈值所需要的扩展周期就越小。这种简单的关系是原始PCR数据分析的基础，也是大部分PCR分析方法的出发点。PCR定性检测通常用于确定样品中是否存在特定的基因或DNA序列，通过比较已知阳性样本和未知样本的扩增曲线，可以确定未知样本中是否存在相同的基因或DNA序列。如果未知样本的扩增曲线与阳性样本的扩增曲线相似或一致，则可以认为未知样本是阳性的，即存在目标基因或DNA序列。但不同医院使用的仪器和试剂不同，扩增效率和不同平台的阈值设置也有差异，这是导致检测结果在不同医院之间存在差异，从而无法互相比对和认可的主要原因。
因此，2009年Stephen A等人发表文章规定了PCR结果想要发表需要提供的最少信息，旨在更好的解读PCR检测结果从而达到实验PCR结果的标准化。与此同时，为了实现PCR检测结果同质化这一目标，大部分国内外的研究学者主要从以下几个传统方面出发：1、标准化方法：研究人员通过制定和推广一致的实验操作方法和分析流程，以确保不同仪器在样品处理和数据分析等方面的一致性。这包括标准化反应体系、引物和探针设计、PCR条件等。2、校准样品：研究人员使用特定的校准样品来评估不同仪器之间的差异，并建立校准曲线或校准模型。这样可以通过校准样品的测量结果来对不同仪器的数据进行校正和比较。3、多中心研究：通过多个独立实验室之间的合作研究，针对不同仪器进行大规模的样本测试和数据比较，以评估不同仪器之间的一致性和可比性。4、质量控制和认证：研究人员对PCR仪器进行质量控制和认证，确保其性能和准确性。这包括仪器的校准和维护，以及参与相关质量认证项目。5、数据分析和统计方法：研究人员开发和应用不同的数据分析和统计方法，以评估不同仪器之间的结果差异，并提供一致性和可比性的指标。
近年来，有部分学者开始尝试通过数学算法来解决不同仪器检测间结果差异这一问题。相比于传统的方法，数学算法只需通过数学计算的过程，从而就可解决不同仪器之间的结果差异问题，无需对仪器进行仪器校正、标准曲线校标等工作，这对于实验室、仪器制造商以及国家医疗系统可以节省大量的人力和财力成本。如Aldo Clerico小组在2020年的研究中采用百分位转换和重新校准数学算法实现TSH检测结果的协调化，这一研究为实现不同仪器检测结果之间具有可比性和一致性提供了新的方向和策略。Aldo Clerico通过收集125，419个来自7个不同实验室的3种检测仪器的TSH检测结果，对这些TSH数值进行百分位数转换和自举分布从而推导出不同的校正方程，以实现不同仪器检测结果的协调化。
结果同质化发展需要多方面的努力和推动，包括标准化和规范化、实验室认可和认证、信息化和数字化、人工智能和大数据分析以及跨学科合作等。目前还没有相关文献和研究报道通过数学算法的方法实现不同仪器PCR检测结果同质化。
发明内容
基于以上问题，本发明的目的在于提供不同PCR仪器间定量或定性检测结果同质化换算方法，实现不同仪器检测结果间具有一致性和可比性，进而实现PCR检测结果相互认可。
对于定量检测，本发明实现其发明目的所采用的技术方案是，不同PCR仪器间定量检测结果同质化换算方法，包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过梯度稀释标准品N个梯度得到的N个样本分别在两台仪器上检测，分别获得参比仪器标准曲线和测试仪器标准曲线，其纵坐标为循环次数Cq值，横坐标为初始样本浓度对数值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
所述标准曲线是仪器根据检测结果输出的曲线；
S2、对待测样本，获得其在测试仪器上的初始样本浓度对数值的实际测试值，找到所述初始样本浓度对数值的实际测试值处于测试仪器标准曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述初始样本浓度对数值的测试值对应到参比仪器标准曲线上的同质化结果值。
等比例计算是指：测试值到测试仪器曲线上对应子直线段第一端点的距离与该子直线段长度的比例，等于同质化结果值到参比仪器曲线上对应子直线段第一端点的距离与该子直线段长度的比例。
进一步，所述比例值，其中XN为待测样本在测试仪器上的初始样本浓度对数值的实际测试值，Xi-1为测试仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，Xi为测试仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
更进一步，所述参比仪器标准曲线上的初始样本浓度对数值的同质化结果值为：，其中/>为参比仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，/>为参比仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
对于定性检测，本发明实现其发明目的所采用的技术方案是，不同PCR仪器间定性检测结果同质化换算方法，包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过Cq值为均匀分布于参比仪器检测范围的N个样本分别在两台仪器上检测，首先分别获得参比仪器扩增曲线和测试仪器扩增曲线，其横坐标为循环次数Cq值，纵坐标为应过程中实时荧光强度为纵坐标，再利用所述参比仪器扩增曲线和测试仪器扩增曲线分别获得参比仪器标准曲线和测试仪器标准曲线，所述参比仪器标准曲线和测试仪器标准曲线各为一根横线，所述横线的横坐标为循环次数Cq值，纵坐标分别等于参比仪器和测试仪器的仪器阈值，所述仪器阈值为提前设定的反应过程中荧光强度阈值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
所述扩增曲线是仪器根据检测结果输出的曲线；而参比仪器标准曲线和测试仪器标准曲线各为一根横线，这两根横线是本发明根据N个样本检测后获得的扩增曲线及仪器阈值而建立的标准曲线，用于后续转换计算；
仪器阈值的默认设置通常是仪器3～15个循环的荧光信号的标准差的10倍；
S2、对待测样本，获得其在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，找到所述实际测试值处于测试仪器曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述荧光信号达到阈值线时对应的循环次数Cq值的实际测试值对应到参比仪器标准曲线上的同质化结果值。
同样的，等比例计算是指：测试值到测试仪器标准曲线上对应子直线段第一端点的距离与该子直线段长度的比例，等于同质化结果值到参比仪器标准曲线上对应子直线段第一端点的距离与该子直线段长度的比例。
进一步，所述比例值，其中XM为待测样本在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，Xj-1为测试仪器标准曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，Xj为测试仪器标准曲线上对应子直线段第二端点的荧光信号达到阈值线时的循环次数Cq值。
更进一步，所述参比仪器曲线上的循环次数Cq同质化结果值为:，其中/>为参比仪器标准曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，/>为参比仪器标准曲线上对应子直线段第二端点的荧光信号达到阈值线时循环次数Cq值。
本发明的有益效果为：
（一）算法创新：通过数学算法解决PCR仪器间检测结果差异导致的结果互认困难的问题。此算法不仅可以用于不同医疗机构间检测结果的比对与相互认可，也可用于科室新使用的仪器的校正工作。本算法将数学算法和检验临床结合应用，不依赖大量的数据和计算工作即可实现转化，具有简单方便的特点，应用场景多种多样，推广实用性强。
（二）临床效益：通过本发明的换算方法实现不同PCR仪器的检测结果互认，可以提高工作效率，对于已经获得的检测、诊断或学习成果，可以不必重复进行相同或类似的评估或测试。这样可以大大提高工作效率，减少不必要的资源和时间浪费。此外，由于避免了不必要的重复评估或测试，结果互认可以显著降低成本。这不仅包括经济成本，还包括患者或个体所需要的时间和精力等非经济成本。
（三）医院效益：结果互认有助于促进机构间的合作。这种互认机制也有助于建立一个更加开放的医疗或教育服务环境，使得服务提供者更加关注质量，而不是重复的评估或测试。
（四）社会效益：结果互认有助于医疗机构间资源共享，为疾病的早期预防、筛查、诊疗提供更加可靠的检测结果，减少重复检验，从而改善人民的就医体验。近年来，为有效利用医疗资源，简化患者就医环节，减少患者不必要的重复性检查，降低患者医疗费用，提高患者满意度。
附图说明
图1为本发明实施例1的PCR定量参比仪器标准曲线和测试仪器标准曲线图（其中，横坐标起始lg DNA量代表样本浓度的对数值，纵坐标Cq代表扩增的循环次数）；
图2为本发明实施例1的PCR定量检测结果转换计算原理示意图（其中，横坐标起始lg DNA量代表样本浓度的对数值，纵坐标Cq代表扩增的循环次数）；
图3为本发明实施例2的PCR定性参比仪器扩增曲线和测试仪器扩增曲线（其中，横坐标Cq代表扩增的循环次数，纵坐标F代表荧光信号值）；
图4为本发明实施例2的PCR定性参比仪器标准曲线和测试仪器标准曲线（其中，横坐标Cq代表扩增的循环次数，纵坐标F代表荧光信号值）；
图5为本发明实施例2的 PCR定性检测结果转换计算原理示意图（其中，横坐标Cq代表扩增的循环次数，纵坐标F代表荧光信号值）。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。
实施例1
图1～2示出本发明的第一种具体实施方式，不同PCR仪器间定量检测结果同质化换算方法，其中我们选取6个比对样本作为示例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述地实施例式本发明一部分实施例，而不是全部地实施例。
此外，图1是我们通过6例样本分别在两台仪器检测所获得的标准曲线示意图，该图表明了两台仪器之间检测相同样本之间的对应关系。图2是我们选取第一段子直线段为例，实施不同PCR仪器间定量检测结果同质化换算。
包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过梯度稀释标准品5个梯度得到的6个样本分别在两台仪器上检测，分别获得参比仪器标准曲线和测试仪器标准曲线，其纵坐标为循环次数Cq值，横坐标为初始样本浓度对数值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
S2、对待测样本，获得其在测试仪器上的初始样本浓度对数值的实际测试值，找到所述初始样本浓度对数值的实际测试值处于测试仪器标准曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述初始样本浓度对数值的测试值对应到参比仪器标准曲线上的同质化结果值。
等比例计算是指：测试值到测试仪器曲线上对应子直线段第一端点的距离与该子直线段长度的比例，等于同质化结果值到参比仪器曲线上对应子直线段第一端点的距离与该子直线段长度的比例。
本实施例中，所述比例值，其中XN为待测样本在测试仪器上的初始样本浓度对数值的实际测试值，Xi-1为测试仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，Xi为测试仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
所述参比仪器标准曲线上的初始样本浓度对数值的同质化结果值为：，其中/>为参比仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，/>为参比仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
如图2所示，假设某一样本N在测试仪器检测初始样本DNA浓度对数值为XN，落在子直线段AB的区间上，现需要使用同质化算法计算出参比仪器上该同一样本初始样本DNA浓度对数值。
第一段AB区间在测试仪器的区间检测结果为：7.11~9.7
第一段A'B'区间在参比仪器的区间检测结果为：7.09~9.64
样本N使用测试仪器检测结果为：XN=8.84
由于A和A'、B和B'是相同的样本在两台仪器分别检测的结果，因此理论上AB区间和A'B'区间上的点之间是一一对应的关系。
故N点在直线AB上的位置的比例和N'点在直线上的位置的比例是相同的。
因为∆ABG和∆NHG是相似三角形，∆A'B'G'和∆N'H'G'也是相似三角形；
则H点在直线GB上的位置的比例等于N点在直线AB上的位置的比例；
H'点在直线G'B'上的位置的比例等于N'点在直线A'B'上的位置的比例；
即可以得到：
代入即可得到：
因此样本N使用等比例计算方法转换至参比仪器的结果为=8.83。
即换算到参比仪器，该同一样本初始样本DNA浓度对数值为8.83。
上述第一种具体实施方式，是不同仪器PCR定量检测结果同质化的算法，用于对待测样本的初始样本DNA浓度对数值进行换算。
实施例2
图3～5示出本发明的第二种具体实施方式，不同PCR仪器间定性检测结果同质化换算方法，其中我们选取5个比对样本作为示例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述地实施例式本发明一部分实施例，而不是全部地实施例。
此外，图3是我们通过5例样本分别在两台仪器检测所获得的扩增曲线示意图。图4是我们通过仪器的扩增曲线的特点，建立了两台仪器的标准曲线的示意图，该标准曲线是一条以y=仪器阈值线的直线。图5是我们选取第一段子直线段为例，实施不同PCR仪器间定性检测结果同质化换算。
包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过Cq值为均匀分布于参比仪器检测范围的N个样本分别在两台仪器上检测，首先分别获得参比仪器扩增曲线和测试仪器扩增曲线，其横坐标为循环次数Cq值，纵坐标为应过程中实时荧光强度为纵坐标，再利用所述参比仪器扩增曲线和测试仪器扩增曲线分别获得参比仪器标准曲线和测试仪器标准曲线，所述参比仪器标准曲线和测试仪器标准曲线各为一根横线，所述横线的横坐标为循环次数Cq值，纵坐标分别等于参比仪器和测试仪器的仪器阈值，所述仪器阈值为提前设定的反应过程中荧光强度阈值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
仪器阈值的默认设置通常是仪器3～15个循环的荧光信号的标准差的10倍；
S2、对待测样本，获得其在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，找到所述实际测试值处于测试仪器标准曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述荧光信号达到阈值线时对应的循环次数Cq值的实际测试值对应到参比仪器标准曲线上的同质化结果值。
同样的，等比例计算是指：测试值到测试仪器标准曲线上对应子直线段第一端点的距离与该子直线段长度的比例，等于同质化结果值到参比仪器标准曲线上对应子直线段第一端点的距离与该子直线段长度的比例。
本实施例中，所述比例值，其中XM为待测样本在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，Xj-1为测试仪器标准曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，Xj为测试仪器标准曲线上对应子直线段第二端点的荧光信号达到阈值线时的循环次数Cq值。
所述参比仪器标准曲线上的循环次数Cq同质化结果值为：，其中/>为参比仪器标准曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，/>为参比仪器标准曲线上对应子直线段第二端点的荧光信号达到阈值线时循环次数Cq值。
如图5所示，假设某一样本M在测试仪器的检测结果为XM，落在子直线段AB区间上，现需要计算出参比仪器上的检测结果。
第一段AB区间在测试仪器的Cq值为：15~20.68
第一段A'B'区间在参比仪器的Cq值为：14~19.82
样本M使用测试仪器检测的Cq值为：XM=16.64
由于A和A'、B和B'是相同的样本在两台仪器分别检测的结果，因此理论上AB区间和A'B'区间上的点之间是一一对应的关系。
故M点在直线AB上的位置的比例和M'点在直线A'B'上的位置的比例是相同的。
则将已知的点的结果代入上式可得：
则样本M使用等比例计算方法转换至参比仪器的结果为=15.68。
即换算到参比仪器，荧光信号达到阈值线时的循环次数Cq值为15.68。
上述第二种具体实施方式，是不同仪器PCR定性检测结果同质化的算法，用于对待测样本达到荧光强度阈值时的循环次数Cq值进行换算。
本发明的上述实施例仅仅是为说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其他不同形式的变化和变动。这里无法对所有的实施方式予以穷举。凡是属于本发明的技术方案所引申出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (6)
1.不同PCR仪器间定量检测结果同质化换算方法，其特征在于，包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过梯度稀释标准品得到的N个样本分别在两台仪器上检测，分别获得参比仪器标准曲线和测试仪器标准曲线，其纵坐标为循环次数Cq值，横坐标为初始样本浓度对数值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
S2、对待测样本，获得其在测试仪器上的初始样本浓度对数值的实际测试值，找到所述初始样本浓度对数值的实际测试值处于测试仪器标准曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述初始样本浓度对数值的测试值对应到参比仪器标准曲线上的同质化结果值。
2.根据权利要求1所述的不同PCR仪器间定量检测结果同质化换算方法，其特征在于：所述比例值其中XN为待测样本在测试仪器上的初始样本浓度对数值的实际测试值，Xi-1为测试仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，Xi为测试仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
3.根据权利要求2所述的不同PCR仪器间定量检测结果同质化换算方法，其特征在于：所述参比仪器标准曲线上的初始样本浓度对数值的同质化结果值为：其中/>为参比仪器标准曲线上对应子直线段第一端点的初始样本浓度对数值，/>为参比仪器标准曲线上对应子直线段第二端点的初始样本浓度对数值。
4.不同PCR仪器间定性检测结果同质化换算方法，其特征在于，包括下述步骤：
S1、选取两台PCR仪器，一台作为参比仪器，另一台作为测试仪器，通过Cq值为均匀分布于参比仪器检测范围的N个样本分别在两台仪器上检测，首先分别获得参比仪器扩增曲线和测试仪器扩增曲线，其横坐标为循环次数Cq值，纵坐标为应过程中实时荧光强度为纵坐标，再利用所述参比仪器扩增曲线和测试仪器扩增曲线分别获得参比仪器标准曲线和测试仪器标准曲线，所述参比仪器标准曲线和测试仪器标准曲线各为一根横线，所述横线的横坐标为循环次数Cq值，纵坐标分别等于参比仪器和测试仪器的仪器阈值，所述仪器阈值为提前设定的反应过程中荧光强度阈值；所述参比仪器标准曲线和测试仪器标准曲线均由多段子直线段组成，每一段子直线段的第一端点对应一个样本，第二端点与下一段子直线段的第一端点重合；
S2、对待测样本，获得其在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，找到所述实际测试值处于测试仪器标准曲线上的具体子直线段，同时找到对应的参比仪器标准曲线上的具体子直线段；
S3、在步骤S2中确定的两段具体子直线段上，利用等比例计算方法，先求出比例值，再计算出所述荧光信号达到阈值线时对应的循环次数Cq值的实际测试值对应到参比仪器标准曲线上的同质化结果值。
5.根据权利要求4所述的不同PCR仪器间定性检测结果同质化换算方法，其特征在于：所述比例值，其中XM为待测样本在测试仪器上荧光信号达到阈值线时对应的循环次数Cq值的实际测试值，Xj-1为测试仪器标准曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，Xj为测试仪器曲线上对应子直线段第二端点的荧光信号达到阈值线时的循环次数Cq值。
6.根据权利要求5所述的不同PCR仪器间定性检测结果同质化换算方法，其特征在于：所述参比仪器标准曲线上的循环次数Cq同质化结果值为：，其中为参比仪器曲线上对应子直线段第一端点的荧光信号达到阈值线时的循环次数Cq值，为参比仪器曲线上对应子直线段第二端点的荧光信号达到阈值线时循环次数Cq值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410397223.XA CN117995279B (zh) | 2024-04-03 | 2024-04-03 | 不同pcr仪器间定量或定性检测结果同质化换算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410397223.XA CN117995279B (zh) | 2024-04-03 | 2024-04-03 | 不同pcr仪器间定量或定性检测结果同质化换算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117995279A true CN117995279A (zh) | 2024-05-07 |
| CN117995279B CN117995279B (zh) | 2024-06-18 |
Family
ID=90887774
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410397223.XA Active CN117995279B (zh) | 2024-04-03 | 2024-04-03 | 不同pcr仪器间定量或定性检测结果同质化换算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117995279B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1588002A (zh) * | 2004-09-16 | 2005-03-02 | 邓平建 | 降低实时荧光pcr仪器定量分析系统误差的分析方法及应用 |
| JP2009103913A (ja) * | 2007-10-23 | 2009-05-14 | Kyocera Mita Corp | 画像形成装置，画像形成装置のパラメータ設定方法 |
| KR20120107717A (ko) * | 2011-03-22 | 2012-10-04 | 삼성테크윈 주식회사 | 실시간 중합효소 연쇄 반응의 증폭 프로파일 곡선으로부터 임계값을 선택하는 방법 |
| WO2016060536A1 (fr) * | 2014-10-17 | 2016-04-21 | Morrocan Foundation For Advanced Science Innovation & Research (Mascir) | Kit de diagnostic du cancer de sein en reaction qpcr multiplex et son application dans l'evaluation du statut de la proteine her2 et le choix du traitement |
| US20170283859A1 (en) * | 2014-10-24 | 2017-10-05 | The Trustees Of Columbia University In The City Of New York | MICROFLUIDIC qRT-PCR ANALYSIS OF SINGLE CELLS |
| JP2019126193A (ja) * | 2018-01-17 | 2019-07-25 | 一般財団法人電力中央研究所 | 情報評価プログラム、情報評価装置及び情報評価方法 |
| CN114166768A (zh) * | 2022-02-14 | 2022-03-11 | 四川大学华西医院 | 不同设备检测同一指标同质化换算方法、装置、电子设备 |
| CN114780905A (zh) * | 2022-06-21 | 2022-07-22 | 四川大学华西医院 | 一种比对样本的确定方法及装置、存储介质、电子设备 |
| CN115048613A (zh) * | 2022-08-16 | 2022-09-13 | 四川大学华西医院 | 一种指标同质化换算方法、装置、电子设备及存储介质 |
| CN116189784A (zh) * | 2023-01-18 | 2023-05-30 | 杭州奥盛仪器有限公司 | 核酸扩增曲线分析方法、装置、设备及存储介质 |
| CN116340331A (zh) * | 2023-03-29 | 2023-06-27 | 杭州电子科技大学 | 一种基于区块链的大型仪器实验结果存证方法和系统 |
| US20230273973A1 (en) * | 2022-02-14 | 2023-08-31 | West China Hospital, Sichuan University | Method and device for homogenization conversion of same index detected by different equipment, and electronic equipment |
-
2024
- 2024-04-03 CN CN202410397223.XA patent/CN117995279B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1588002A (zh) * | 2004-09-16 | 2005-03-02 | 邓平建 | 降低实时荧光pcr仪器定量分析系统误差的分析方法及应用 |
| US20060058970A1 (en) * | 2004-09-16 | 2006-03-16 | Pingjian Deng | Analytical method for reducing quantitative analysis systematic error of real-time fluorescence thermal cycler and the applications of the same |
| JP2009103913A (ja) * | 2007-10-23 | 2009-05-14 | Kyocera Mita Corp | 画像形成装置，画像形成装置のパラメータ設定方法 |
| KR20120107717A (ko) * | 2011-03-22 | 2012-10-04 | 삼성테크윈 주식회사 | 실시간 중합효소 연쇄 반응의 증폭 프로파일 곡선으로부터 임계값을 선택하는 방법 |
| WO2016060536A1 (fr) * | 2014-10-17 | 2016-04-21 | Morrocan Foundation For Advanced Science Innovation & Research (Mascir) | Kit de diagnostic du cancer de sein en reaction qpcr multiplex et son application dans l'evaluation du statut de la proteine her2 et le choix du traitement |
| US20170283859A1 (en) * | 2014-10-24 | 2017-10-05 | The Trustees Of Columbia University In The City Of New York | MICROFLUIDIC qRT-PCR ANALYSIS OF SINGLE CELLS |
| JP2019126193A (ja) * | 2018-01-17 | 2019-07-25 | 一般財団法人電力中央研究所 | 情報評価プログラム、情報評価装置及び情報評価方法 |
| CN114166768A (zh) * | 2022-02-14 | 2022-03-11 | 四川大学华西医院 | 不同设备检测同一指标同质化换算方法、装置、电子设备 |
| US20230273973A1 (en) * | 2022-02-14 | 2023-08-31 | West China Hospital, Sichuan University | Method and device for homogenization conversion of same index detected by different equipment, and electronic equipment |
| CN114780905A (zh) * | 2022-06-21 | 2022-07-22 | 四川大学华西医院 | 一种比对样本的确定方法及装置、存储介质、电子设备 |
| CN115048613A (zh) * | 2022-08-16 | 2022-09-13 | 四川大学华西医院 | 一种指标同质化换算方法、装置、电子设备及存储介质 |
| CN116189784A (zh) * | 2023-01-18 | 2023-05-30 | 杭州奥盛仪器有限公司 | 核酸扩增曲线分析方法、装置、设备及存储介质 |
| CN116340331A (zh) * | 2023-03-29 | 2023-06-27 | 杭州电子科技大学 | 一种基于区块链的大型仪器实验结果存证方法和系统 |
Non-Patent Citations (2)
| Title |
|---|
| ZHOU, X等: "A mathematical algorithm to harmonize measurements for thyroid-stimulating hormone between instruments", CLINICA CHIMICA ACTA, vol. 548, 1 August 2023 (2023-08-01) * |
| 刘在栓等: "基于互联网医院实现线上检验项目同质化的管理构想", 国际检验医学杂志, vol. 43, no. 15, 15 August 2022 (2022-08-15), pages 1911 - 1914 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117995279B (zh) | 2024-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110349629A (zh) | 一种利用宏基因组或宏转录组检测微生物的分析方法 | |
| JP6695899B2 (ja) | 標準物質として外部生体分子を用いた生体分子分析方法およびそのキット | |
| Barnes et al. | Whole tumor section quantitative image analysis maximizes between-pathologists' reproducibility for clinical immunohistochemistry-based biomarkers | |
| Scott Jr et al. | Assessing the reproductive competence of individual embryos: a proposal for the validation of new “-omics” technologies | |
| CN116052778A (zh) | 实时监测生物反应器中细胞培养液组分浓度的方法 | |
| CN117995279B (zh) | 不同pcr仪器间定量或定性检测结果同质化换算方法 | |
| CN102866258B (zh) | 一种泌尿系结石病因诊断系统及其应用方法 | |
| CN101825558B (zh) | 一种研究聚丙烯腈基碳纤维热稳定化进程的方法 | |
| Rajasekharan et al. | From bench to MS bedside: challenges translating biomarker discovery to clinical practice | |
| Jing et al. | Reference materials for molecular diagnostics: Current achievements and future strategies | |
| CN107630097A (zh) | 一种快速检测泌尿生殖道感染常见病原体的引物组及试剂盒 | |
| CN117079723A (zh) | 一种与肌萎缩侧索硬化症相关的生物标志物、诊断模型及其应用 | |
| CN111896497A (zh) | 一种基于预测值的光谱数据修正方法 | |
| CN104087672A (zh) | 一种多重实时荧光定量pcr技术快速检测人21号染色体数目的试剂盒 | |
| CN105603081A (zh) | 一种肠道微生物定性与定量的检测方法 | |
| CN102367494B (zh) | 用于鼠源rna病毒多重实时荧光pcr检测的引物、探针及其试剂盒 | |
| CN113990485A (zh) | 一种基于常规检验数据与知识图谱联合诊断疾病的方法 | |
| CN117012274B (zh) | 基于高通量测序识别基因缺失的装置 | |
| Yang et al. | A comparison of the automated blood bank system IH-500 and manual tube method for anti-blood group antibody titration: a quantitative approach | |
| CN102128927B (zh) | 一种艾滋病检测质控血清的配制方法 | |
| Caraguel et al. | A modelling approach to predict the variation of repeatability and reproducibility of a RT-PCR assay for infectious salmon anaemia virus across infection prevalences and infection stages | |
| CN116837142B (zh) | 一种断裂点簇集区艾贝尔逊白血病病毒校准品及其制备方法和应用 | |
| CN107312863A (zh) | 利用数字PCR鉴定B淋巴细胞的Ig基因分布的方法 | |
| Thomas et al. | Unlocking the reversal potential of solid supported membrane electrophysiology to determine transport stoichiometry | |
| Wang et al. | Integrating Patient-Based Real-Time Quality Control (PBRTQC) in a New Field: Inter-Comparison between Biochemical Instrumentations with LDL-C |
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 |