Opisana powyżej metoda ilościowej RT-PCR została wdrożona w celu wykrywania i oznaczania ilościowego wirusa herpes koni typu 2 (EHV-2) w płynach z dróg oddechowych. Rysunek 1 przedstawia schematyczny wykres przebiegu prac nad opracowaniem i walidacją metody ilościowej RT-PCR zgodnie z normą AFNOR NF U47-600. Specyficzność starterów i sond została zweryfikowana podczas etapowego opracowywania procedury PCR. W tym systemie amplifikowano wyłącznie szczepy EHV-2. Następnie należało scharakteryzować wydajność metody qRT-PCR.
W pierwszej kolejności, w celu oszacowania LODPCR, wykonano 6 dziesięciokrotnych rozcieńczeń seryjnych, aby wyznaczyć strefę wygasania (Ryc. 2). W tym przykładzie, w celu oszacowania LODPCR, przygotowano 6 dziesięciokrotnych rozcieńczeń seryjnych w zakresie od 10-5 do 10-10 (od 26 000 do 0,26 kopii/2,5 µl próbki). Strefa wygasania znajduje się pomiędzy rozcieńczeniami 10-9 a 10-10 (między 2,6 a 0,26 kopii/2,5 µl próbki). Aby w tym przypadku wyznaczyć wartość LODPCR , w obrębie tej strefy wygasania wykonano 6 dwukrotnych rozcieńczeń seryjnych plazmidu w zakresie od 5,2 do 0,16 kopii/2,5 µl próbki. Wartość LODPCR wyniosła 2,6 kopii/2,5 µl próbki.
Aby określić zakres liniowości oraz LOQPCR, wartość LODPCR została wykorzystana jako punkt wyjścia dla serii 6 rozcieńczeń dziesięciokrotnych, w zakresie od 2,6 (LODPCR) do 260 000 kopii/2,5 µl próbki. Rycina 3 przedstawia regresję liniową dla qRT-PCR EHV2 z jednego badania. Parametry regresji liniowej (Rycina 4) zostały zweryfikowane w czterech powtórzeniach przy użyciu obliczeń opisanych w Tabeli 3. Obliczenia przeprowadzono w celu zdefiniowania zakresu liniowości zgodnie z kryterium wartości bezwzględnej Biasi ≤0,25 log10, niezależnie od poziomu i obciążenia plazmidowego. W tym przypadku zakres liniowości mieścił się między 2,6 a 260 000 kopii/2,5 µl próbki. LOQPCR to najniższe stężenie w zakresie liniowości (t.j. w tym przypadku 2,6 kopii/2,5 µl próbki). ULIN określono na poziomie 0,12 log10 w zakresie 2,6-260 000 kopii/2,5 µl DNA.
Po opracowaniu (Rycina 1, kolor niebieski) i charakterystyce qRT-PCR (Rycina 1, kolor żółty), norma AFNOR NF U47-600 zaleca charakterystykę całej metody analitycznej, od ekstrakcji DNA po qRT-PCR (Rycina 1, kolor pomarańczowy). Czułość i swoistość diagnostyczną obliczono zgodnie z opisem w Tabeli 4. Wydajność ilościowa całej metody analitycznej qRT-PCR została oceniona i zwalidowana za pomocą profilu dokładności (Rycina 5).
Niniejszy protokół, wykorzystujący najnowocześniejszą technologię molekularną, pozwolił nam wykryć i określić ilość kopii genomu wirusa EHV-2 w 172 próbkach wymazów z nozdrzy pobranych od koni z zaburzeniami oddechowymi i/lub klinicznym podejrzeniem infekcji. Częstość występowania EHV-2 w próbkach terenowych (biologicznych) w tej populacji wyniosła 50% (86/172). Analizy ilościowe wykazały, że ładunek genomu wirusa EHV-2 był znacznie wyższy u młodych koni, a rozkład ładunków genomu wirusa malał wraz z wiekiem (Rycina 6). W niniejszym badaniu najwyższy ładunek genomu wirusa EHV-2 (1.9 x 1011 kopii/ml) stwierdzono u źrebiąt (Rycina 6).

Rysunek 1: Schemat blokowy opracowania (kolor niebieski), charakterystyki ilościowej RT-PCR (kolor żółty) oraz charakterystyki całej metody analitycznej od ekstrakcji DNA do qRT-PCR (kolor pomarańczowy) zgodnie z normą AFNOR NF U47-600-2. Schemat blokowy podsumowuje poszczególne etapy opracowania, charakterystyki ilościowej RT-PCR oraz charakterystyki całej metody analitycznej od ekstrakcji DNA do qRT-PCR. Dla każdego etapu schemat określa liczbę wymaganych serii pomiarowych, rozcieńczenia do wykonania oraz liczbę wymaganych analityków. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 2: Wyznaczenie obszaru zaniku na podstawie reprezentatywnych wyników krzywych PCR w czasie rzeczywistym uzyskanych dla 6 dziesięciokrotnych rozcieńczeń seryjnych plazmidu. Aby oszacować obszar zaniku, przygotowano 6 dziesięciokrotnych rozcieńczeń seryjnych w zakresie od 10-5 (26 000 kopii/2,5 µl próbki) do 10-10 (0,26 kopii/2,5 µl próbki). Obszar zaniku znajduje się pomiędzy rozcieńczeniami 10-9 (2,6 kopii/2,5 µl próbki) a 10-10 (0,26 kopii/2,5 µl próbki). W tym przypadku w obrębie tego obszaru zaniku wykonano 6 dwukrotnych rozcieńczeń seryjnych plazmidu, aby wyznaczyć LOD 95% PCR, pomiędzy 5,2 a 0,16 kopii/2,5 µl próbki. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Regresja liniowa dla qRT-PCR EHV2. Liniowość badania ilościowego to zdolność do generowania wyników proporcjonalnych do stężenia celu w określonym zakresie. Można ją modelować za pomocą regresji liniowej (y = ax + b) pomiędzy odpowiedzią aparaturową (próg cyklu lub Ct) a logarytmem ilości celu (liczba kopii celu/2.5 µl próbki). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Wyniki regresji liniowej qPCR dla EHV-2. Średni błąd systematyczny reprezentuje średnią różnicę między zmierzoną ilością plazmidu (
) a teoretyczną ilością plazmidu (x'i) dla każdego poziomu plazmidu. Pionowe słupki reprezentują niepewność liniowości (ULINi) wyznaczoną według wzoru

gdzie SD'i to odchylenie standardowe zmierzonej ilości plazmidu. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Profile dokładności opracowane na podstawie wyników walidacji metody qRT-PCR dla EHV-2. Zielona linia (kółka) reprezentuje prawdziwość danych (błąd systematyczny, czyli obciążenie). Granice akceptowalności zostały zdefiniowane przez laboratorium na poziomie ± 0.75 Log10 (linie przerywane). Dolne i górne granice dokładności zostały określone dla każdego poziomu zawartości plazmidu jako średnie obciążenie ± dwukrotność odchylenia standardowego danych niezawodności (linie czerwone). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 6: Ilościowe oznaczenie ładunków genomu wirusa EHV-2 w zależności od wieku. Przedstawiono rozkład ładunku genomu wirusa EHV-2 wykrytego w próbkach wymazów z nosa dla różnych grup wiekowych. Linie poziome reprezentują wartości mediany wraz z odchyleniem standardowym (m = miesiące). * Różnica istotna statystycznie w analizie ANOVA z testem post-hoc Newmana-Keulsa (p <0.05). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Gen docelowy | Sekwencje starterów, sondy i plazmidu (5'-3') | Pozycja nukleotydu | Wielkość produktu (nukleotydy) | Warunki cyklu termicznego | Bibliografia |
EHV2 gB
(HQ247755.1) | Forward: GTGGCCAGCGGGGTGTTC | 2113-2130 | 78 | 95 °C 5 min | | 11 |
| Odwrócona sekwencja: AAGTTYTT AGGGAAACCCCC | 2189-2170 | 95 °C 15 s | 45 cykli |
| Sonda: FAM-CCCTCTTTGGGAGCATAGTCTCGGGG-MGB | 2132-2157 | 60 °C 1 min |
Plazmid:
ACCTGGGCACCATAGGCAAGGTGGTGGTCA
ATGTGGCCAGCGGGGTGTTCTCCCTCTTTG
GGAGCATAGTCTCGGGGGTGATAAGCTTTTT
CAAAAATCCCTTTGGGGGCATGCTGCTCATA
GTCCTCATCATAGCCGGGGTAGTGGTGGTG
TACCTGTTTATGACCAGGTCCAGGAGCATAT
ACTCTGCCCCCATTAGAATGCTCTACCCCGG
GGTGGAGAGGGCGGCCCAGGAGCCGGGCG
CGCACCCGGTGTCAGAAGACCAAATCAGGA
ACATCCTGATGGGAATGCACCAATTTCAG | 2081-2381 | | |
Tabela 1: Sekwencje starterów, sond oraz pozytywnych syntetycznych kontroli DNA użytych w tym protokole. Sekwencja plazmidu (pozytywnej syntetycznej kontroli DNA) odpowiada pozycjom nukleotydowym 2081-2381 sekwencji EHV2gB (HQ247755.1). Projekt starterów i sond użytych w tym protokole został opracowany przy użyciu specjalistycznego oprogramowania.
| PATOGENY | Referencja (pochodzenie) | Liczba szczepów | WYNIKI |
| EHV-2 |
| EHV-2 | VR701 (ATCC) | 20 | Dodatnie |
| 20 próbek (kolekcja FDL) |
| EHV-5 | KD05 (GERC) | 20 | Ujemne |
| 20 próbek (kolekcja FDL) |
| EHV-3 | VR352 (ATCC) | 2 | Ujemne |
| T934 WSV (GERC) |
| EHV-1 | Szczep Kentucky Ky A (ATCC) | 3 | Ujemne |
| 2 próbki (kolekcja FDL) |
| EHV-4 | VR2230 (ATCC) | 1 | Ujemne |
| Herpesvirus osłów AHV5 | Kolekcja FDL | 1 | Ujemne |
| Wirus grypy koni | A/equine/Jouars/4/2006 (H3N8) | 1 | Ujemne |
| (Numer dostępu JX091752) |
| Wirus arteritis koni | VR796 (ATCC) | 2 | Ujemne |
| Rhodococcus equi | Kolekcja FDL | 1 | Ujemne |
| Streptococcus equi subsp. Zooepidemicus | Kolekcja FDL | 1 | Ujemne |
| Streptococcus equi subsp. equi | Kolekcja FDL | 1 | Ujemne |
| Coxiella burnetii | ADI-142-100 (Adiagene) | 1 | Ujemne |
| Chlamydophila abortus | ADI-211-50 (Adiagene) | 1 | Ujemne |
| Klebsiella pneumoniae | Kolekcja FDL | 1 | Ujemne |
Tabela 2: Specyficzność analityczna qRT-PCR dla EHV-2.

Tabela 3: Obliczanie błędu systematycznego (bias) i niepewności liniowości (adaptowano z NF U47-600-212). Dla każdej próby parametry regresji liniowej (y = ax+b) są walidowane przy użyciu tabeli, gdzie y to uzyskany próg cyklu (cycle threshold); a to uzyskany nachylenie; x to poziom plazmidu, a b to punkt przecięcia z osią rzędnych. i oznacza poziom plazmidu (i przyjmuje wartości od 1 do k poziomów); k to liczba użytych poziomów plazmidu (np. k = 6 w tej tabeli); j to próba (j przyjmuje wartości od 1 do I prób); I to liczba prób, mieszcząca się w przedziale od 3 do 6 (np. I = 4 w tej tabeli). xi to szacowana ilość plazmidu dla każdego poziomu plazmidu i. x'i to teoretyczna ilość plazmidu uzyskana za pomocą równania x'i = log10(xi) dla każdego poziomu plazmidu i. Podczas każdej próby j próg cyklu uzyskany dla każdego poziomu plazmidu i jest obliczany za pomocą regresji liniowej yi,j = ajxi,j + bj.
to zmierzona ilość plazmidu podczas próby j. Biasi to różnica zaobserwowana między zmierzoną ilością plazmidu a teoretyczną ilością plazmidu dla każdej próby i każdego poziomu plazmidu.
to wartość średnia z
dla każdego poziomu plazmidu i; SD'i to odchylenie standardowe zmierzonej ilości
dla każdego poziomu plazmidu i; Mean bias to średnia z Biasi; ULINi to niepewność liniowości określona dla każdego poziomu plazmidu i, obliczona z SD'i oraz mean bias. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ilustracji.
| | Rzeczywisty stan próbki |
| | Dodatni | Ujemny |
| Wyniki uzyskane przy zastosowaniu pełnej metody | Dodatni | RP (prawdziwie dodatni) | FP (wynik fałszywie dodatni) |
| Negatywny | wynik fałszywie ujemny (FN) | RN (prawdziwie ujemny) |
| Suma | RP+FN | FP+RN |
| | Se = RP/(RP+FN) | Sp = RN/(RN+FP) |
Tabela 4: Obliczenie czułości diagnostycznej (Se) i swoistości (Sp) całej metody. Do obliczenia 95% przedziału ufności dla czułości i swoistości całej metody, zgodnie z opisem w NF U47-600-2, wykorzystano tabelę Schwartza.