Conference materials UDC 57.087.1
DOI: https://doi.org/10.18721/JPM.153.269
Processing of qPCR signals obtained on microfluidic chips in the measurement sequence disorder event
D. Yu. Klimenko \ N. A. Esikova 2, A. L. Bulyanitsa 2, R. V. Davydov 1 3, D. A. Belov 2e
1 Peter the Great St. Petersburg Polytechnic University, St. Petersburg , Russia; 2 Institute of Analytical Instrumentation of the Russian Academy of Sciences, St. Petersburg , Russia; 3 All-Russian Research Institute of Phytopathology, Moscow Region, Russia H [email protected]
Abstract. The formation and movement of gas bubbles in microfluidic chips leads to the measurement sequence disorder in the form of a signal distortion. Two new methods for automatic processing of distorted qPCR signals were developed and implemented. Methods are based on the qPCR signals approximation by a sigmoid function and make it possible to successfully perform a microfluidic chips qPCR analysis in the event of bubbles, which is demonstrated on experimental and simulation curves.
Keywords: microfluidics, microfluidic chip, qPCR, bubbles, approximation
Funding: The work was carried out within the framework of the state task of the Ministry of Science and Higher Education of the Russian Federation 075-00761-22-00.
Citation: Klimenko D. Yu., Esikova N. A., Bulyanitsa A. L., Davydov R. V., Belov D. A., Processing of qPCR signals obtained on microfluidic chips in the measurement sequence disorder event, St. Petersburg State Polytechnical University Journal. Physics and Mathematics. 15 (3.2) (2022) 375-380. DOI: https://doi.org/10.18721/JPM.153.269
This is an open access article under the CC BY-NC 4.0 license (https://creativecommons. org/licenses/by-nc/4.0/)
Материалы конференции УДК 57.087.1
DOI: https://doi.org/10.18721/JPM.153.269
Обработка полученных на микрофлюидных чипах сигналов ПЦР-РВ при разладке последовательности измерений
Д. Ю. Клименко 1И, Н. А. Есикова 2, А. Л. Буляница 2, Р. В. Давыдов 13, Д. А. Белов 2 1 Санкт-Петербургский политехнический университет Петра Великого, Санкт-Петербург, Россия;
2 Институт аналитического приборостроения Российской академии наук, Санкт-Петербург, Россия;
3 Всероссийский научно-исследовательский институт фитопатологии, Московская область, Россия
н [email protected] Аннотация. Одним из основных недостатков микрофлюидных чипов является образование пузырей, приводящее к возникновению разладки в последовательности измерений в форме искажения сигналов ПЦР-РВ. Разработаны и реализованы две новые методики автоматической обработки искаженных сигналов, основанные на аппроксимации сигмоидальной функцией. Выполнена апробация методов на экспериментальных и имитационных сигналах. Данные методики позволяют получать удовлетворяющие значения погрешностей при определении порогового цикла и успешно выполнять количественный анализ ПЦР-РВ в микрофлюидных чипах при возникновении пузырей.
Ключевые слова: микрофлюидика, микрофлюидный чип, ПЦР-РВ, пузыри, аппроксимация
© Klimenko D. Yu., Esikova N. A., Bulyanitsa A. L., Davydov R. V., Belov D. A., 2022. Published by Peter the Great St. Petersburg Polytechnic University.
Финансирование: Работа выполнена в рамках Госзадания Минобрнауки РФ 07500761-22-00.
Ссылка при цитировании: Клименко Д. Ю., Есикова Н. А., Буляница А. Л., Давыдов Р. В., Белов Д. А. Обработка полученных на микрофлюидных чипах сигналов ПЦР-РВ при разладке последовательности измерений // Научно-технические ведомости СПбГПУ. Физико-математические науки. Т. 15. № 3.2. С. 375-380. DOI: https://doi.org/10.18721/ JPM.153.269
Статья открытого доступа, распространяемая по лицензии CC BY-NC 4.0 (https:// creativecommons.org/licenses/by-nc/4.0/)
Introduction
Recently, a large number of methods [1-4] and instruments [5-7] for diagnosing condensed matter have been developed in the world. The main disadvantages of such instruments are their large size, high consumption of reagents and samples, etc. [8, 9]. The use of microfluidic chips for the study of liquids and solutions in microsystems allows solving many of these problems. These chips provide reproducible control of the laminar flows of nano- and picoliter volumes of liquid in microsized channels, high sensitivity and low power, reagents and samples consumption [10-12]. So they are widely used in chemical engineering, pharmaceuticals, biotechnology, and medicine [13].
Quantitative polymerase chain reaction (qPCR) is a popular method of genetic analysis, which makes it possible to quantify the content of the target nucleotide sequence in a sample by determining quantification cycle values C [14]. Microfluidic qPCR chips are being actively developed and improved in the world [15-20]. A significant disadvantage of performing PCR in a chip chamber compared to test tubes is the formation of gas bubbles during thermal cycling. It has no noticeable effect on the reaction result, but prevents correct optical detection. Several methods have been used to inhibit the bubble generation: (i) the design of PCR chamber, (ii) the surface treatment, (iii) the sealing pressurization of the PCR chamber (iv) degasification of the PCR sample and (v) the addition of high boiling-point biocompatible reagents to the PCR sample [13]. Mathematical processing of qPCR signals makes it possible to evaluate the analysis result, despite the gas bubbles formation influence on signal. The PCR efficiency evaluation depending on the materials and design of the chips, the reagents used, and the reaction conditions is complex due to the occurrence of discord in the measurement sequence due to the formation and movement of bubbles. This results in a fluorescence intensity decrease over several cycles. Therefore, the development of methods for processing of qPCR signals distorted due to the formation of bubbles is an extremely urgent task.
Materials and Methods
Microchips for qPCR analysis were made of polycarbonate Novattro (SafPlast, Russia) and polypropylene PP 4445S (PJSC Nizhnekamskneftekhim, Russia) by thermal pressing in a MM-100 hydraulic press (MTDI, Korea) on a stainless steel master mold made by laser microprocessing. Microchip topology represents 3 chambers with supply channels [21]. Polymer film P-500 (LLC PKF Modern technologies, Russia) was used for sealing microchips.
The obtained microchips were used for qPCR of soybean DNA isolated manually using a set of M-sorb-OOM (LLC Sintol, Russia). The experiments were carried out for the initial concentration of isolated nucleic acids, for dilutions of 1, 2 and 3 orders, and for a positive control in three repeats for each variant. R6G and Cy5 dyes were used to detect soybean DNA samples and confirm a positive control respectively.
The measurements were carried out on a specially designed model that provides a thermal cycling mode for qPCR, dyes excitation at wavelengths of 530, 570 and 685 nm and signal registration at 580, 630 and 660 nm. The microchip chambers were filled with a reaction mixture with the addition of nucleic acids; the inputs/outputs were sealed with PCR film. The microchip was placed on the heating element with the film facing down. Overall, 114 experimental qPCR plots were obtained, about 70 % had distortion caused by bubbles.
© Клименко Д. Ю., Есикова Н. А., Буляница А. Л., Давыдов Р. В., Белов Д. А., 2022. Издатель: Санкт-Петербургский политехнический университет Петра Великого.
Two signal processing methods, based on the qPCR plots approximation by a known sigmoid function [22] (formula 1) were developed and implemented in the MATLAB software to identify the correct C values:
F„ =■
F
Cq -C
+ fo,
(1)
1 + e
where Fc is the fluorescence signal at cycle C, RFU; Fm is the maximal reaction fluorescence, RFU; C is the fractional cycle at which reaction fluorescence reaches half of F ; k is related to
" q J m~
the slope of the curve, f is the fluorescence background, RFU.
Each cycle of the program, when implementing the developed methods, the graph is approximated. In method No. 1, the signal value with the largest absolute deviation from the approximating function (1) AF is assigned the value corresponding to the approximating dependence point. In method No. 2, the value corresponding to the deviation AF is excluded from the experimental set. The program end is determined automatically.
The use of the program on experimental data does not allow estimating the errors in determining the threshold cycle values since there is no a priori information about the true Cq values. The experimental curves errors study was carried out, which revealed the following:
- the mathematical expectation of errors close to zero;
- the distribution of errors is more consistent with the normal law: information discrepancy of the Kullback—Leibler histogram of the real distribution of interference over 6 intervals with a uniform law of 0.0710 nits, with a normal one — 0.0422 nits, a coefficient of kurtosis — 3.07;
- the error is close to multiplicative, its standard deviation is possibly proportional to the value of the average signal.
Based on the results of the error study, Set A was created from 10000 simulation curves (Fig. 1,a) with parameters Fm = 1000; Cq = 30; k = 1, f = 600, the dependence law of errors normally distributed standard deviations (RMS) on the signal value o = 0.0024Fc + 1.1368, uniform distributions of the duration of signal distortion from 3 to 7 cycles, and the intensity of the signal drop when a bubble occurs (0.22-0.24)-Fc. The beginning of the bubble is also random, its distribution is even.
We separately created five sets of 10000 curves: with signal distortion caused by bubbles that occur in the ranges: Set B from 1 to 24 cycles, Set C from 25 to 35 cycle and Set D from 36 to 50 cycle; as well as with signal distortions of duration: Set E 3 cycles and Set F 7 cycles.
Results and Discussion
As a result of applying the method No. 1 (Fig. 1,b) on Set A errors 76.41 % of values C lie in the range ±0.15 cycle from the true value, as a result of applying method No. 2, 77.39 %. The mean Cq scores for the methods were 30.36 ± 0.84 and 30.52 ± 1.30, respectively. The estimate bias is caused by large positive errors, the proportion of which is about 20 %. The acceptable value of 0.15 cycles corresponds to a quantification error of less than 10 % and was chosen as the deviation in the determination of Cq, which is acceptable for most tasks of qPCR analysis. According to this criterion, method No. 2 has advantages.
a)
b)
Cycles number
Cycles number
Fig. 1. Graphs for qPCR: model graph (a); model graph with eliminated distortion according to Method No. 1 (b)
k
Fig. 2. Results of Method No. 1 application in Set C: threshold cycle position estimates histogram (a); possible probabilistic model for the distribution of scores (b), which is the sum of two Gaussian distributions
In case of distortion in the initial (Set B) or final (Set D) section of the kinetic curve, the average errors in estimating the position of the inflection point are at an acceptable level: 0.013 for the initial section and 0.022 for the final section when choosing any of the two methods. Permissible error in determining C obtained in 100 % of cases.
The main errors occur when the signal is distorted in the central section of the qPCR graph. 11% of the curves fall within the allowable range of ±0.15 cycles (Fig. 2).
The mathematical expectation (M) of the threshold cycle position in method No. 1 was 31.18 cycles. As follows from the histogram, the distribution is close to bell-shaped. According to the Novickij-Zograf rule, at least 90 % of the sample must fall within the interval M ± m^o standard deviation at m = (1.65 ± 0.05) [23], in our case in the interval {29.81; 32.55} with m = 1.65.
Set analysis shows that 9329 out of 10000 results fell into this interval (i.e., 93.3 %), which correlates well with the Novickij-Zograf rule. Depending on the method for constructing an approximating dependence and estimating the inflection point of the kinetic curve in the presence of bubbles in the section 25—35 cycles, the estimation error are given in Table 1. The average error according to method No. 1 is 1.18 cycles, according to method No. 2 0.99 cycles.
With method No. 2, the proportion of small (up to 0.3 cycles) errors is much larger (29.5 versus 18.5 %). However, the proportion of super-large estimation errors is almost the same: 20 % of the largest errors with method No. 1 - 1.89 cycles or more, with method No. 2 - 1.80 cycles or more. Approximately every 10th quantitative qPCR analysis, when distortion is on the central part, we will mistakenly exclude a reasonable and correct result.
Another stage of the simulation involved introducing a discord corresponding to a 3 or 7 distortion duration (Sets E and F).
Table 1
Comparison of the errors distribution in determining C by two methods when processing the Set C
Method No. 1 Method No. 2
Level, cycle Probability, % Level, cycle Probability, %
0.30 18.5 0.30 29.5
0.50 25.2 0.50 38.0
1.00 42.6 0.85 50.0
1.12 50.0 1.00 56.8
1.89 80.0 1.80 80.0
According to method No. 1 with the replacement of values with a distortion duration of 3 cycles, the average error in threshold cycle determining was 0.135; with an alternative method of exclusion, it was 0.065 cycles. That is, for a "small" bubble, it is more reasonable to exclude the corresponding 3 measurements from the signal.
When simulating a bubble with distortion duration of 7 cycles, the average errors in estimating the position of the threshold cycle are practically the same with an accuracy of several thousandths and are approximately equal to 0.254 cycles.
Conclusion
Two automated methods have been developed that allow to obtain a satisfying RMS value when determining the threshold cycle and successfully perform a quantitative analysis of qPCR in microfluidic chips when bubbles occur. Approbation on a set of 10000 created simulation curves showed that the acceptable error in determining Cq values, method No. 1 provides in 76.41 % of cases, method No. 2 - in 77.39 %. The main errors occur when the distortion is in the central section of the curve, which is characterized by the most intense fluorescence intensity change. When using the developed methods, 11 % of the processed curves fall within the allowable range of Cq value errors. Methods allow identifying such cases and recommending the operator to discard the values as false.
As a result, developed methods approbation on sets of simulation curves with different signal distortion lengths revealed that for short distortion durations, it is more reasonable to exclude the corresponding measurements from the signal, for large ones, the choice of the method is not fundamental.
REFERENCES
1. Zheng Q., Chu W., Zhao C., Zhang L., Guo H., Wang Y., Jiang X., Zhao J., Ab initio nonadiabatic molecular dynamics investigations on the excited carriers in condensed matter systems. Wiley Interdisciplinary Reviews: Computational Molecular Science. 9(8) (2019) e1411.
2. Makeev S. S., Grevtzeva A. S., Glinushkin A. P., Matorin D. N., Possibilities of using spectral analysis in method of nuclear magnetic spectroscopy for condensed media investigation, Journal of Physics: Conference Series. 1695(1) (2020) 012112.
3. Kuzmin M. S., Rogov S. A., On the use of a multi-raster input of one-dimensional signals in two-dimensional optical correlators, Computer Optics. 43(3) (2019) 391—396.
4. Mazing М. S., Zaitceva A. Yu., Kislyakov Yu. Ya., Kondakov N. S., Avdushenko S. A., Davydov V. V., Monitoring of Oxygen Supply of Human Tissues Using A Noninvasive Optical System Based on A MultiChannel Integrated Spectrum Analyzer, International Journal of Pharmaceutical Research. 12(2) (2020) 1974-1978.
5. Grebenikova N. M., Smirnov K. J., Rud V. Yu., Artemiev V. V., Features of monitoring the state of the liquid medium by refractometer, Journal of Physics: Conference Series. 1135(1) (2018) 012055.
6. Myazin N. S., Yushkova V. V., Rud V. Y., On the possibility of recording absorption spectra in weak magnetic fields by the method of nuclear magnetic resonance, Journal of Physics: Conference Series. 1038(1) (2018) 012088.
7. Grebenikova N. M., Smirnov K. J., Features of optical signals processing for monitoring the state of the flowing liquid medium with a refractometer, Journal of Physics: Conference Series. 1368(2) (2019) 022057.
8. Davydov V.V., Dudkin V.I., Karseev A.Y., A Compact Nuclear Magnetic Relaxometer for the Express Monitoring of the State of Liquid and Viscous Media, Measurement Techniques. 57(8) (2014) 912-918.
9. Karseev A., Vologdin V., Davydov V., Features of nuclear magnetic resonance signals registration in weak magnetic fields for express - Control of biological solutions and liquid medium by nuclear magnetic spectroscopy method, Journal of Physics: Conference Series. 643(1) (2015) 012108.
10. Evstrapov A.A., Microfluidic chips for biological research and research, Russian Chemical Journal. 55 (2) (2011) 99-110.
11. Bruus H., Theoretical microfluidics. Lecture notes third edition. MIC Department of Micro and Nanotechnology Technical University of Denmark, 2006.
12. Fiorini G. S., Chiu D. T., Disposable microfluidic devices: fabrication, function and application, BioTechniques. 38(3) (2005) 429-446.
13. Zhang C., Xing D., Miniaturized PCR chips for nucleic acid amplification and analysis: latest advances and future trends, Nucleic Acids Research. 35(13) (2007) 4223—4237.
14. Nakayama T., Kurosawa Y., Furui S., Kerman K., Kobayashi M., Rao S. R., Tamiya E., Circumventing air bubbles in microfluidic systems and quantitative continuous-flow PCR applications. Analytical and Bioanalytical Chemistry, 386(5) (2006) 1327-1333.
15. Fedorov A. A., Berdnikov A. S., Kurochkin V. E., The polymerase chain reaction model analyzed by the homotopy perturbation method, Journal of Mathematical Chemistry. 57(4) (2019) 971-985.
16. Natyrov A. N., Vlasova N. A., Matvienko I. V., Kurochkin V. E., Alexeev J. I., Synthesis of Unsymmetrical Polymethine Cyanine Fluorescent Dyes for Nucleic Acid Analysis by Real-Time PCR, Russian Journal of Bioorganic Chemistry. 44(5) (2018) 562-571.
17. Reznik V. S., Kruglov V. A., Petrov A. I., Glinuchkin A. P., Rud V. Y., Development of a measuring device for the study of thermal processes during the polymerase chain reaction, Journal of Physics: Conference Series. 1410(1) (2019) 012078.
18. Kiselev I. G., Belov D. A., Physical processes simulation in a precision device for liquid samples thermal cycling, Journal of Physics: Conference Series. 2131(2) (2021) 022061.
19. Matvienko I. V., Bayramov V. M., Parygina N. A., Kurochkin V. E., Alekseev Y. I., Synthesis of Dihydroquinoline-Based Derivatives of Fluorescent Rhodamine Dyes for Nucleic Acid Analysis by a Real-Time Polymerase Chain Reaction, Russian Journal of Bioorganic Chemistry. 46(3) (2020) 349-359.
20. Zubik A. N., Rudnitskaya G. E., Bulyanitsa A. L., Lukashenko T. A., Evstrapov A.A., Microfluidic chips for real-time PCR, Journal of Physics: Conference Series. 2086 (2021) 012124.
21. Yesikova N.A., Germash N.N., Evstrapov A.A., Operational production of microchips for PCR analysis from polymer materials in laboratory conditions, Scientific instrumentation. 30 (4) (2020) 21-6.
22. Novickij P.V., Zograf I.A., Ocenka pogreshnostej rezul'tatov izmerenij [Assessment of errors of observed data], Nauka, Leningrad, 1991.
23. Rutledge R.G., Stewart D. A., Kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR. BMC Biotechnol. (2008) 8:47.
THE AUTHORS
KLIMENKO Daria
[email protected] ORCID: 0000-0002-1216-2192
[email protected] ORCID: 0000-0003-1958-4221
DAVYDOV Roman
ESIKOVA Nadezhda
[email protected] ORCID: 0000-0002-0145-6451
BELOV Dmitrii
ORCID: 0000-0003-3219-0446
BULYANITSA Anton
[email protected] ORCID: 0000-0002-9235-8549
Received 14.08.2022. Approved after reviewing 16.08.2022. Accepted 16.08.2022.
© Peter the Great St. Petersburg Polytechnic University, 2022