Rycina 2. Wpływ dwóch różnych filtrów przestrzennych HYPR na wygładzenie danych aktywności w czasie dla pojedynczego woksela prążkowia. Górny rząd: obrazy emisji PET z 11C-rakloprydy z 3-minutowej ramki wycentrowanej na 46,5 minucie (bez filtrowania, filtrowane jądrem 3 x 3 x 3 woksele, filtrowane jądrem 5 x 5 x 5 woksele). Środkowy rząd: obrazy emisji PET z 11C-rakloprydy z 3-minutowej ramki wycentrowanej na 61,5 minucie (bez filtrowania, filtrowane jądrem 3 x 3 x 3 woksele, filtrowane jądrem 5 x 5 x 5 woksele). Dolny rząd: odpowiadające im krzywe aktywności w czasie dla tego samego pojedynczego woksela w lewym grzbietowym jądrze ogonowym. Należy zauważyć, że widoczny spadek wychwytu 11C-rakloprydy (wynikający z uwalniania dopaminy) w momencie palenia zostaje zachowany, mimo że szum jest zmniejszony przy większym rozmiarze filtra.
Rysunek 3. Wybór reprezentatywnych funkcji odpowiedzi dopaminergicznej, które zostały obliczone wcześniej w celu dopasowania modelu lp-ntPET do danych czasowo-aktywnościowych PET w każdym wokselu zgodnie z Normandin et al. 7. W przypadku naszego paradygmatu palenia, palenie papierosa rozpoczyna się 45 minut po rozpoczęciu wstrzyknięcia znacznika. Nawet jeśli odpowiedzi dopaminergiczne w prążkowiu kodują oczekiwanie na palenie – np. z powodu kontaktu z papierosem lub innych wskazówek zapowiadających palenie – uznaliśmy, że funkcje odpowiedzi można bezpiecznie ograniczyć do krzywych, które odchylają się od linii bazowej nie wcześniej niż 5 minut przed paleniem (a). Podobnie ograniczono krzywe do czasów rozpoczęcia odpowiedzi nie późniejszych niż 15 minut po rozpoczęciu palenia. Krzywe z czasem rozpoczęcia odpowiedzi na poziomie 40 minut reprezentują możliwe odpowiedzi dopaminergiczne wynikające z oczekiwania. (b) Reprezentatywne funkcje odpowiedzi, wszystkie odchylające się od linii bazowej w 45 minucie; czyli w czasie rozpoczęcia palenia. Wygenerowano 500 różnych prawdopodobnych funkcji odpowiedzi. Dla celów ilustracyjnych wykresy w (a) i (b) pokazują jedynie próbkę kształtów krzywych i czasów rozpoczęcia odpowiedzi.
Rycina 4. (a) Równanie operacyjne dla modelu lp-ntPET. Model jest liniowy w parametrach (R1, k2, k2a, γ), co pozwala na szybkie obliczanie szacunków parametrów dla każdego woksela w obrębie maski prążkowia. (b) Obrazy parametryczne (R1, k2, k2a, γ) dla pojedynczego przekroju coronalnego mózgu u jednego badanego. Chociaż tylko γ jest parametrem kodującym wielkość odpowiedzi dopaminergicznej, jednoczesna estymacja wszystkich 4 parametrów znacznika jest niezbędna do dopasowania modelu do danych czasowo-aktywnościowych w każdym wokselu.
Rysunek 5. Dopasowanie konwencjonalnego modelu (MRTM) oraz nowego modelu (lp-ntPET) do danych aktywności w czasie z woksela w lewym jądrze ogoniastym. Dopasowanie MRTM zaznaczono kolorem niebieskim. Dopasowanie lp-ntPET zaznaczono kolorem czerwonym.
Rysunek 6. (a) przedstawia ważoną sumę kwadratów reszt (WSSR) z dopasowań MRTM, a (b) z dopasowań lp-ntPET do danych w każdym wokselu prążkowia. Dwa obrazy WSSR wygenerowane z tych samych danych są porównywane w celu stworzenia mapy stosunku F w każdym wokselu (tj. mapy F), pokazanej na (c). (d) Mapa F zostaje poddana progowaniu przy p < 0,05 w celu utworzenia binarnej mapy istotności (patrz krok 2.10 protokołu). Dla n ramek czasowych i 4 parametrów modelu lp-ntPET, próg statystyki F odpowiadający poziomowi prawdopodobieństwa p < 0,05 (dla danych z 90 minut podzielonych na 3-minutowe ramki próg wynosi 4,23) (e) Mapa istotności jest filtrowana za pomocą filtra morfologicznego („otwieranie”), aby wyeliminować niewielkie skupiska wokseli, które najprawdopodobniej reprezentują szum. Końcowa Maska Istotności zachowuje tylko te woksele w prążkowiu, których TAC są (statystycznie) lepiej dopasowane do modelu lp-ntPET niż do konwencjonalnego modelu MRTM, a zatem uznaje się, że zawierają odpowiedź dopaminergiczną na palenie. Ten próg nie uwzględnia korekty na wielokrotne porównania. Zamiast tego, aby zapobiec wynikom fałszywie dodatnim, tworzymy Końcowe Maski Istotności również dla warunku kontrolnego (patrz Rysunek 7 oraz kroki protokołu 1.8 – 1.10).
Rysunek 7. (a) przedstawia jeden przekrój wieńcowy Końcowej Maski Istotności dla warunku palenia u pojedynczego badanego. Rysunek (b) przedstawia Końcową Maskę Istotności dla odpowiadającego badanego i przekroju w warunku bazowym. Obecność skupisk zachowanych wokseli w masce palenia w przeciwieństwie do niemal całkowitego braku skupisk w masce kontrolnej wspiera tezę, że filmy dopaminowe (patrz poniżej) nie są po prostu zdarzeniami przypadkowymi ani zdarzeniami związanymi z szumem w danych. (Uwaga: wstrzyknięta aktywność – a tym samym stosunek sygnału do szumu – w warunkach bazowych i palenia była porównywalna.)
Rysunek 8. Film obrazujący poziom dopaminy w pojedynczym przekroju mózgu w orientacji wieńcowej przedstawia klatka po klatce poziom dopaminy w stosunku do podstawowego (spoczynkowego) poziomu dopaminy. (a) przedstawia film dla stanu bazowego, a (b) przedstawia film dla stanu palenia. Poziomy dopaminy są zakodowane kolorami. W szczególności kolory – przedstawione na pasku kolorów z odpowiadającymi im wartościami liczbowymi – reprezentują zmianę poziomu dopaminy powyżej poziomu bazowego jako procent wartości bazowej. Ponownie, poziomy dopaminy są pokazane tylko dla wokseli w końcowej masce istotności (Final Significance Mask), które przekraczają poziom istotności p < 0,05.
Rycina 9. Wieloprzekrojowy i wielowariantowy film obrazujący poziom dopaminy u tego samego osobnika co na Rycynie 8, z jednoczesnym wyświetleniem wszystkich przekrojów prążkowia brzusznego dla warunków bazowych i palenia.

Rycina 1. Schemat blokowy procedur eksperymentalnych i analizy obrazów (a-c). Kliknij tutaj, aby wyświetlić powiększoną rycinę.

Rycina 2. Wpływ filtrów HYPR o różnych rozmiarach jądra na obrazy (góra i środek) oraz na krzywe aktywności czasowej (dół) w pojedynczym wokselu. Kliknij tutaj, aby wyświetlić powiększoną rycinę.

Rysunek 3. Przykłady funkcji odpowiedzi dopaminy, które rozpoczynają się w (a) 40 min lub (b) 45 min po wstrzyknięciu znacznika.

Rycina 4. Obrazy parametryczne wygenerowane poprzez dopasowanie równania operacyjnego lp-ntPET (a) do danych PET. (b) Obrazy odpowiadające 4 parametrom modelu, R1, k2, k2a, γ, zostały ocenione dla prążkowia i przedstawiono je na nałożonych odpowiadających im przekrojach MR. Kliknij tutaj, aby wyświetlić większą rycinę.

Rysunek 5. Dopasowania modeli MRTM (niebieski) i lp-ntPET (czerwony) do danych czasowo-aktywnościowych z pojedynczego woksela.

Rysunek 6. Obrazy parametryczne WSSR dla (a) MRTM oraz (b) lp-ntPET. Odpowiednie mapy WSSR są porównywane w celu stworzenia mapy F (c), która z kolei jest poddawana progowaniu w celu uzyskania maski binarnej (d), a następnie filtrowaniu w celu wygenerowania końcowej maski istotności (Final Significance Mask). Kliknij tutaj, aby wyświetlić powiększony rysunek.

Rysunek 7. Porównanie końcowych masek istotności dla warunku palenia (a) oraz kontrolnego (b) u tego samego badanego.
Rysunek 8. Film przedstawiający poziom dopaminy w jednym przekroju dla jednego badanego w warunkach kontrolnych („Odpoczynek”) i palenia. Kliknij tutaj, aby wyświetlić Rysunek 8.
Rycina 9. Wieloprzekrojowy film dopaminowy dla pojedynczego obiektu w warunkach (na górze) palenia i (na dole) kontrolnych („Odpoczynek”). Kliknij tutaj, aby wyświetlić Rycina 9.