ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2011, том 21, № 2, c. 27-33
= ОБЗОРЫ, ИССЛЕДОВАНИЯ, ПРИБОРЫ =
УДК 543.426; 543.9
© Ю. В. Белов, А. И. Петров, В. В. Лавров, В. Е. Курочкин
ИЗУЧЕНИЕ ВЛИЯНИЯ ШУМОВ ДЕТЕКТОРА НА ПОГРЕШНОСТИ КОЛИЧЕСТВЕННЫХ АНАЛИЗОВ НУКЛЕИНОВЫХ КИСЛОТ
НА ПРИБОРАХ ПЦР-РВ
В статье выполнен сравнительный анализ эффективности способов фильтрации и аппроксимации сигналов ПЦР-РВ при количественных анализах нуклеиновых кислот. Показано, что разные способы аппроксимации могут быть дополнены предварительной фильтрацией сигналов ПЦР. Даны рекомендации по выбору уровня порога и участков для нормировки сигналов. Полученные результаты могут быть использованы для выбора объема реактивов в пробирках при количественных анализах ПЦР на основе критериев допустимого уровня шума и допустимой погрешности измерения.
Кл. сл.: ПЦР в реальном времени, фильтрация сигналов, пороговый цикл, количественные анализы
ВВЕДЕНИЕ
При разработке анализаторов нуклеиновых кислот основной задачей является уменьшение объема реактивов в пробирках для количественных анализов методом полимеразной цепной реакции в реальном времени (ПЦР-РВ). При меньшем объеме реактивов уменьшается стоимость каждого анализа и повышается быстродействие амплифи-катора за счет уменьшения тепловой инерции. Нижний допустимый предел объема реактивов в пробирках при количественных анализах методом ПЦР-РВ ограничен чувствительностью флуори-метрического детектора. В работе [1] был выполнен практический анализ погрешностей количественных измерений методом ПЦР-РВ при различных способах расчета порогового цикла. Выяснены вклады различных составляющих погрешностей при количественных измерениях. Показано, что из-за разброса свойств реактивов возникают трудности при выделении вкладов в погрешности измерений, вносимых шумами флуориметрическо-го детектора. В настоящей статье выполнено моделирование сигналов ПЦР-РВ при различных уровнях шума. Модельные сигналы использованы для сравнения эффективности способов фильтрации и аппроксимации сигналов при количественных анализах нуклеиновых кислот.
ЭФФЕКТИВНОСТЬ ФИЛЬТРАЦИИ СИГНАЛОВ ПЦР-РВ
В основе метода ПЦР-РВ лежит наблюдение в реальном времени сигналов флуоресценции в ходе реакции. Концентрация исходных специфических фрагментов ДНК увеличивается приблизи-
2п
, где п — текущее количество циклов. Процесс завершения ПЦР обусловлен расходом реагентов при накоплении специфических продуктов амплификации. Поэтому в процессе увеличения количества циклов можно условно выделить 3 участка. На первом участке сигнал флуоресценции специфических фрагментов ДНК находится на уровне, ниже шума. Второй участок можно назвать информационным: на нем наблюдается рост сигнала флуоресценции. На третьем участке увеличения сигнала флуоресценции не происходит.
Исходный модельный сигнал, приведенный на рис. 1, был получен из реальных сигналов ПЦР-РВ анализатора нуклеиновых кислот АНК-32, который серийно выпускается ИАП РАН [2]. Для этого был вычислен усредненный нормированный сигнал ПЦР от 31 пробирки с одинаковыми реактивами; первому участку в течение 14 температурных циклов усредненного сигнала присвоены нулевые значения, а третьему участку от 30 до 50 циклов были присвоены значения 1000 относительных единиц (о.е.).
Если обозначить значение модельного сигнала флуоресценции 7п (п — номер цикла), то в качестве полезного аналитического сигнала ПЦР можно выбрать приращение интенсивности флуоресценции L\ или L5 за 1 или 5 циклов в середине графика, где наблюдается максимальная крутизна.
Максимальное приращение исходного модельного сигнала ПЦР за один цикл (цикл № 24) равно
Ll = 724 - ^23 = 598 - 416 = 182 о.е.
Приращение за пять циклов (циклы № 22-26) равно
L5 = 726 - 721 = 864 - 150 = 714 о.е.
^ 1200
<и о
1| 1000
(и
^ 800
и
400 200 0
0 20 40 60
Номер цикла
Рис. 1. Исходный модельный сигнал ПЦР-РВ
В большинстве случаев в результате фильтрации приращение сигнала флуоресценции L1f или L5f уменьшается. В качестве исходных модельных шумовых сигналов использовались 10 последовательностей 50 случайных чисел (рис. 2). Интенсивности исходных шумовых сигналов определялись путем вычисления стандартного отклонения 5". При максимальной амплитуде однополярного шума, равной одной относительной единице,
стандартное отклонение 5 = 0.31 о.е., а после фильтрации величина 5т < 5.
С целью определения улучшения отношения сигнала к шуму К (эффективность фильтрации) предложены следующие предварительные критерии:
- относительное уменьшение полезного сигнала ПЦР L0 = Lf / L;
- относительное изменение интенсивности шума 50 = 5т / 5.
Применение относительных величин позволяет исключить влияние масштабов исходных величин. В результате эффективность фильтрации К (о.е.) можно выразить следующим образом:
К = L / 5.
Желательно, чтобы после фильтрации величина L0 незначительно отличалась от 1, а величина 50 была минимальной.
С целью сравнения эффективности были использованы следующие фильтры [3-6]:
- скользящее среднее значение (прямоугольное окно) по 3 циклам (далее — С3);
- последовательное применение этого фильтра 2 и 3 раза (2С3 и 3С3);
- скользящее среднее значение по 5 циклам (С5);
- фильтр Savitzky-Golay (Савицкого—Голея, СГ);
- спектральный фильтр (СП);
- медиана по 3 и 5 точкам (М3 и М5).
1.2
0 10 20 30 40 50 60
Номер цикла
Рис. 2. Один из исходных модельных шумовых сигналов. По горизонтальной оси — номер цикла п, по вертикальной оси — величина шума в относительных единицах. Для наглядности дискретные значения шумового сигнала условно соединены линиями
Результаты фильтрации сигнала 7п для цикла с номером п могут быть выражены в следующем виде:
7С3, п = (7п-1 + 7п +7п+1)/3;
72С3, п = (7п—2+ 27п-1 + 37п + 27п+1 + 7п+2 )/9;
73С3, п = (7п_3 + 37п_2 + 67п_1 + 77п + 67п+1 + +37п+2 + 7п+3)/27;
7С5, п = (7п—2 + 7п-1 + 7п + 7п+1 + 7п+2)/5;
7сг, п = (-27п_3 + 37п_2 + 67п_1 + 77п + 67п+1 + + 37п+2 -2 7п+3)/21.
Фильтр СГ называют также фильтром полиномиального цифрового сглаживания на основе метода наименьших квадратов (реализован в программе МайаЬ).
Спектральный фильтр был реализован следующим способом. После фурье-преобразования действительная и мнимая части спектра умножались на спадающую экспоненциальную функцию ЕХР(-0.1(/- 1)), где / — порядковый номер спек-
тральной точки; например, при / = 2 значение ЕХР(-0.1) = 0.9048 соответствует ослаблению величины сигнала второй точки по сравнению с величиной сигнала первой точки спектра. В качестве фильтрованного сигнала ПЦР использовалась реальная часть обратного фурье-преобразования.
Медиана выделяет число, которое является серединой упорядоченного множества значений нечетного числа членов.
Результаты применения фильтрации к модельным сигналам ПЦР и шума приведены в табл. 1.
Видно, что медиана совершенно не изменяет приращение интенсивности флуоресценции (¿0 = = 1). Это свойство медианы означает, что шумы фильтруются только на части сигнала вне этого участка. Лучшие результаты (значения К больше 2) получены при использовании 3-кратной фильтрации прямоугольным окном по 3 точкам. Можно предположить, что величина К1 будет иметь существенное значение при последующей аппроксимации графика по 2 точкам, а К5 — при аппроксимации графика по многим точкам.
Табл. 1. Сравнение эффективности фильтров К (о.е.)
Фильтр Lol 50 К1 К5
С3 0.92 0.96 0.52 1.77 1.85
2С3 0.85 0.92 0.43 1.98 2.14
3С3 0.80 0.89 0.38 2.11 2.34
С5 0.79 0.88 0.40 1.98 2.20
СГ 0.93 0.99 0.57 1.63 1.74
СП 0.8 0.86 0.44 1.82 1.95
М3 1 1 0.71 1.41 1.41
М5 1 1 0.62 1.61 1.61
Рис. 3. Графики модельных сигналов ПЦР с разными уровнями шума
ПОГРЕШНОСТИ РАСЧЕТА ПОРОГОВОГО ЦИКЛА БЕЗ ФИЛЬТРАЦИИ
Для оценки влияния шума на погрешности измерения пороговых циклов исходный модельный сигнал, содержащий одно числовое значение при каждом температурном цикле, суммировался со случайными числами. При амплитуде исходного модельного сигнала 1000 относительных единиц (о.е.) выбирались следующие максимальные величины однополярных случайных чисел: 25, 50 и 100 о.е. В результате получались 3 блока суммарных модельных сигналов с разными величинами случайных чисел. Каждый из 10 сигналов в одном блоке отличается от другого последовательностями случайных чисел. Три блока суммарных модельных сигналов и исходный модельный сигнал вводились в программу анализатора АНК-32 в качестве имитации сигналов от пробирок.
В программе анализатора была выполнена нормировка: привязка к нулю среднего значения сигналов от 1 до 14 циклов и привязка к 1000 среднего значения сигналов от 30 до 50 циклов. При нормировке на начальном участке шумовые сигналы становятся симметричными относительно нуля. При этом максимальные амплитуды шума
для трех наборов сигналов соответствуют 12.5, 25 и 50 о.е. На рис. 3 приведены одновременно графики всех суммарных модельных сигналов и линия на уровне 500 о.е.
Анализаторы АНК позволяют рассчитать пороговые циклы методом порога, используя разные способы аппроксимации: способ кусочно-линейной аппроксимации по двум точкам (в дальнейшем — АДТ), способ аппроксимации по всем точкам графика (АВТ), а также способ аппроксимации на базе S-образной функции (АСФ) [7]. Для сравнения этих методов выполнена оценка погрешности измерения путем расчета значений стандартного отклонения 5 пороговых циклов 3 блоков, содержащих по 10 суммарных модельных сигналов.
Результаты расчета погрешностей пороговых циклов приведены в табл. 2. Погрешности пороговых циклов в табл. 2 и следующей таблице вычислены как отношение к одному циклу; например, 5 = 0.029 соответствует погрешности 2.9 % от одного цикла.
Видно, что величины стандартного отклонения при различных способах расчета пороговых циклов имеют близкие значения и приблизительно линейно зависят от уровня шума.
Табл. 2. Значения стандартного отклонения 5 при различных способах расчета пороговых циклов в зависимости от уровня шума при пороге 500 о.е.
Амплитуда шума (о.е.) Способы расчета пороговых циклов
5 при АДТ 5 при АВТ 5 при АСФ
12.5 0.029 0.022 0.021
25 0.056 0.047 0.047
50 0.117 0.095 0.099
Табл. 3. Величины стандартного отклонения 5 и отклонения среднего значения D при различных способах фильтрации и способах расчета пороговых циклов
Фильтры Критерии АДТ АВТ АСФ
Нет 5 0.123 0.094 0.080
Б 0.036 0.017 0.019
3С3 5 0.092 0.083 0.079
Б 0.002 -0.014 0.012
СГ 5 0.102 0.093 0.084
Б 0.001 0.015 0.008
S
0.7 т
0.6
0.5
0,4 -
0,3
0.2
0.1
—Р?щ 1 -»-Ряд 2 —Ü— Ряд 3 -Ö— Рад 4
Ii Iii i i i i
0 1 00 200 300 400 500 600 700 800 900 1 000 Уровень порога (o.e.)
Рис. 4. Зависимости значений стандартного отклонения от уровня порога.
Ряд 1 — АДТ без фильтра; ряд 2 — АДТ с фильтром 3С3; ряд 3 — АВТ с фильтром 3С3; ряд 4 — АСФ с фильтром 3С3
ВЛИЯНИЕ ФИЛЬТРАЦИИ НА ПОГРЕШНОСТИ РАСЧЕТА ПОРОГОВЫХ ЦИКЛОВ
Для достижения минимальной погрешности пороговых циклов необходимо оптимизировать выбор трех уровней на графиках сигналов ПЦР: уровни привязки к нулю и к максимальному значению, а также уровень порога.
Для определения уровней привязки к нулю и к максимальному значению можно рекомендовать участки, на которых даже после фильтрации не заметно изменение сигнала флуоресценции. Необходимо также исключать несколько точек в начале и конце графика, которые не могут быть правильно отфильтрованы (некоторые фильтры усредняют значения по 7 температурным циклам). Уровни привязки вычисляются как средняя величина сигналов на выбранном участке, т. е. добавляется еще одно звено фильтрации шумов. Рекомендации для выбора оптимального уровня порога (500 о.е.) приведены ниже.
С целью сравнения эффективности фильтров при разных методах расчета были использованы 10 модельных сигналов при уровне шума 50 о.е.
При расчете значений 5 пороговых циклов для привязки к нулю выбран участок от 5 до 10 циклов, для привязки к 1000 — участок от 35 до 40 циклов. Уровень порога выбран 500 о.е. Результаты расчета приведены в табл. 3. Видно, что наибольшая величина 5 наблюдается при способе АДТ без фильтрации, а наименьшая величина — при АСФ с фильтром 3С3. Фильтр 3С3 уменьшает величину 5 при АДТ на 25 %, а при АСФ — только на 1 %.
В этой таблице приведен еще один критерий для сравнения влияния фильтров — отклонение среднего значения D вычисленных пороговых циклов с использованием фильтров от порогового цикла исходного сигнала. Можно отметить, что величина D значительно меньше 5 и поэтому не вносит заметного вклада в погрешности измерения.
В табл. 3 приведены значения 5 при пороге 500 о.е., однако при изменении уровня порога величины стандартного отклонения изменяются (рис. 4).
На графике ряд 1 (АДТ без фильтра) величины 5 сильно зависят от уровня порога. При использовании фильтрации сигналов ПЦР зависимости ве-
личин 5 от уровня порога для разных способов аппроксимации практически совпадают, при этом минимальные величины 5 наблюдаются при пороге 500 о.е. Этот уровень порога можно считать оптимальным.
ВЫВОДЫ
С помощью предложенных модельных сигналов установлено, что при увеличении уровня шума происходит приблизительно линейное увеличение стандартного отклонения пороговых циклов.
На основе предложенных критериев выполнено сравнение эффективности нескольких типов фильтров сигналов ПЦР-РВ. Лучшие результаты получены при использовании 3-кратной фильтрации с помощью скользящего среднего значения по 3 пороговым циклам.
Исследовано влияние фильтрации на погрешности расчета пороговых циклов. Даны рекомендации по выбору участков для нормировки сигналов. Показано, что разные способы аппроксимации могут быть дополнены предварительной фильтрацией сигналов ПЦР, при этом способ аппроксимации по двум точкам по величине стандартного отклонения приближается к двум другим более сложным способам аппроксимации.
Исследована зависимость погрешности расчета пороговых циклов от уровня порога. В результате при расчетах может быть рекомендован уровень порога 500 о.е., при котором получаются минимальные погрешности измерения.
Полученные результаты позволяют пользователю выбрать объем реактивов в пробирках при количественных анализах ПЦР на основе критериев: допустимого уровня шума и допустимой погрешности измерения. Например, в соответствии с табл. 1 при максимальной амплитуде шума 12.5 о.е. и сигнале 1000 о.е. стандартное отклонение 5 имеет приблизительно величину 0.02 или 2 % порогового цикла. Примерно такие результаты получаются при реальных анализах ПЦР с объемом реактивов 23 мкл и объемом пробы 2 мкл.
СПИСОК ЛИТЕРАТУРЫ
1. Белов Ю.В., Петров А.И., Лавров В.В., Куроч-кин В.Е. Особенности количественных измерений содержания нуклеиновых кислот методом полиме-разной цепной реакции в реальном времени // Научное приборостроение. 2011. Т. 21, № 1. С. 44-49.
2. Алексеев Я.И., Белов Ю.В., Варламов Д.А. и др. Приборы для диагностики биологических объектов на основе метода полимеразной цепной реакции в реальном времени (ПЦР-РВ) // Научное приборостроение. 2006. Т. 16, № 3. С. 132-136.
3. Родионова О.Е., Померанцев А.Л. Хемометрика в аналитической химии. Институт химической физики им. Н.Н. Семенова РАН.
URL: (http://www. chemometrics.ru/materials/articles/ chemometrics_review.pdf).
4. Savitzky A., Golay M.J.E. // Anal. Chem. 1964. V. 36. P. 1627.
5. Крылов А.С., Втюрин А.Н., Герасимова Ю.В. Обработка данных инфракрасной Фурье-спектроскопии. Методическое пособие. Препринт № 832 Ф. Красноярск, Институт физики СО РАН, 2005.
URL: (http://www.kirensky.ru/papers/preprint/832.pdf).
6. Фурье-спектроскопия ядерного магнитного резонанса. Лабораторная работа СПбГУ.
URL: (http://lab2.phys.spbu.ru/pdf_to/NMR.pdf).
7. 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.
URL: ( http://www.biomedcentral.com/1472-6750/8/47).
Институт аналитического приборостроения РАН, г. Санкт-Петербург
Контакты: Белов Юрий Васильевич, [email protected]
Материал поступил в редакцию 15.02.2011.
INFLUENCE OF DETECTOR NOISE ON ERRORS OF QUANTITATIVE ANALYSIS OF NUCLEIC ACIDS USING REAL-TIME PCR
Yu. V. Belov, A. I. Petrov, V. V. Lavrov, V. E. Kurochkin
Institute for Analytical Instrumentation RAS, Saint-Petersburg
Comparative analysis of efficiency for filtering methods and approximation of real-time PCR signals in studying of quantitative analysis of nucleic acids are discussed. It's shown that different approximation methods
could be supplemented with preliminary PCR signals filtration. Recommendations for choice of threshold level and sections for signals normalization are given. Results received by us will be used for determination of volume of reagents for PCR quantitative analysis based on criteria of noise permissible level and permissible errors of measuring.
Keywords: real-time PCR, signal filtering, threshold cycle, quantitative analyis