Uczestnicy
Uczestników zrekrutowano z Oddziału Neurologii Szpitala Miejskiego w Ningde oraz z Uniwersytetu Normalnego w Ningde w okresie od listopada 2022 do listopada 2023 roku. Zgodnie z kryteriami włączenia i wykluczenia, do badania włączono 35 pacjentów (16 z AD, 19 z MCI) oraz 25 zdrowych osób starszych stanowiących grupę kontrolną (NC). Od wszystkich uczestników oraz od rodzin osób z zaburzeniami poznawczymi uzyskano świadomą zgodę. Badanie zostało przeprowadzone zgodnie z zasadami Deklaracji Helsińskiej. Zgodę wyraziła Komisja Etyki Szpitala Miejskiego w Ningde (numer zgody: 20200901).
Kryteria włączenia obejmowały wiek od 50 do 85 lat, gotowość do poddania się badaniu MRI i wyrażenie świadomej zgody oraz możliwość ukończenia obrazowania z pomocą rodziny lub w ramach wielu sesji, jeśli było to konieczne. Kryteria wykluczenia obejmowały obecność guzów mózgu; ciężką chorobę sercowo-naczyniową; historię chorób psychiatrycznych; przebyty uraz czaszkowo-mózgowy; główne zaburzenia neurologiczne, takie jak infekcje ośrodkowego układu nerwowego, choroba Parkinsona, choroba Huntingtona, epilepsja lub stwardnienie rozsiane; ciężkie choroby ogólnoustrojowe; nienachalne implanty metaliczne; ciężkie upośledzenia sensoryczne uniemożliwiające ocenę poznawczą; artefakty obrazowania oraz leworęczność.
Diagnozy choroby Alzheimera (AD), łagodnych zaburzeń poznawczych (MCI) oraz funkcji poznawczych w normie (NC) zostały postawione przez dwóch neurologów zgodnie z kryteriami National Institute on Aging–Alzheimer’s Association (NIA-AA). Prawdopodobna AD wymagała postępującego pogorszenia funkcji poznawczych i funkcjonalnych; MCI wymagała obiektywnego upośledzenia funkcji poznawczych przy zachowaniu codziennego funkcjonowania; NC wymagała braku nieprawidłowości poznawczych lub neurologicznych. W przypadku wszystkich uczestników osiągnięto konsensus diagnostyczny. Potwierdzenie biomarkerami nie było dostępne w sposób jednolity, co stanowi ograniczenie badania. Wszystkie narzędzia i materiały użyte w tym badaniu są wymienione w Tabeli materiałów.
Dane kliniczne
Zebrano dane demograficzne i kliniczne, w tym wiek, płeć, wykształcenie, historię medyczną (nadciśnienie tętnicze, cukrzyca, hiperlipidemia) oraz wyniki badania fizykalnego. Funkcje poznawcze oceniono za pomocą mini-badania stanu psychicznego (MMSE) oraz Montrealnej skali oceny funkcji poznawczych (MoCA).
Nabywanie obrazów
Obrazowanie wykonano za pomocą skanera MRI 3.0T z zastosowaniem dedykowanej cewki głowy i szyi. Protokół obrazowania obejmował sekwencje T2-weighted fluid-attenuated inversion recovery (T2-FLAIR), obrazowanie dyfuzyjne (DWI), obrazowanie T1-weighted w płaszczyźnie strzałkowej oraz MRI z kontrastem fazowym (obejmujące 2D cine phase-contrast oraz sekwencje 4D flow). Wszyscy uczestnicy zostali ułożeni w pozycji na plecach, głową do przodu, z głową unieruchomioną za pomocą piankowych podkładek i z zastosowaniem ochrony słuchu w celu zminimalizowania artefaktów ruchowych.
Dwuwymiarowe obrazowanie MR w trybie cine z kontrastem fazowym
Dwuwymiarową sekwencję cine z kontrastem fazowym (Fast PC Cine Slice) wykonano prostopadle do środkowej części wodociągu mózgu, zgodnie z oznaczeniem na obrazach lokalizacyjnych T1-zależnych w płaszczyźnie strzałkowej. Parametry obrazowania były następujące: pole widzenia, 24 × 24 cm; grubość warstwy, 5.0 mm; czas repetycji/czas echa, 20.2/6.1 ms; macierz, 288 × 256; rozmiar piksela, 0.83 × 0.94 mm; szerokość pasma, 31.25 kHz; liczba uśrednień, 1.0; kąt odchylenia, 20°; prędkość kodowania, 12 cm/s; synchronizacja obwodowa; oraz 20–30 faz na cykl serca. Prędkość kodowania dobrano tak, aby odpowiednio zarejestrować przepływ płynu mózgowo-rdzeniowego w wodociągu bez występowania artefaktów aliasingu. Czas akwizycji wynosił od 1 do 3 min, w zależności od tętna uczestnika.
Czterowymiarowy rezonans magnetyczny przepływu (4D flow MRI)
W celu objęcia obszaru szyjnego tętnic szyjnych wykonano sekwencję 4D flow (Sag 4D Flow). Parametry obrazowania były następujące: pole widzenia (FOV), 24 × 24 cm; grubość warstwy, 1.0 mm; odstęp między warstwami, 2.0 mm; 40–50 warstw; czas repetycji/czas echa (TR/TE), 6.6/2.8 ms; macierz, 180 × 180; rozmiar piksela, 1.33 × 1.33 mm; szerokość pasma, 62.50 kHz; liczba ekscytacji, 1.0; kąt odchylenia, 15°; prędkość kodowania, 50 cm/s; oraz 20–30 faz na cykl serca. Prędkość kodowania 50 cm/s wybrano, aby odpowiednio zarejestrować przepływ w tętnicach szyjnych, minimalizując jednocześnie zjawisko aliasingu. Czas akwizycji wynosił od 1 do 3 minut, w zależności od częstości akcji serca. Wszystkie obrazy przeanalizowano przy użyciu specjalistycznego oprogramowania do postprocessingu (Tabela materiałów).
Przetwarzanie i analiza obrazów
Analiza parametrów związanych z płynem mózgowo-rdzeniowym
Dwuwymiarowe obrazy amplitudy i fazy cine PC analizowano przy użyciu specjalistycznego oprogramowania do postprocessingu (Rysunek 1A–D). Dwóch radiologów (każdy z trzyletnim doświadczeniem w neuroobrazowaniu) niezależnie wyznaczył obszary zainteresowania (ROI) wzdłuż zewnętrznej krawędzi akweduktu mózgu w miejscu jego maksymalnego przekroju poprzecznego, zidentyfikowanego na lokalizatorze T1-zależnym w płaszczyźnie strzałkowej i potwierdzonym na obrazie amplitudy. Każdy radiolog wyznaczył ROI dwukrotnie, z odstępem co najmniej dwóch tygodni między sesjami. Korekcję fazy tła przeprowadzono przy użyciu wbudowanego w oprogramowanie algorytmu korekcji tkanek stacjonarnych, który odejmuje średnią fazę regionu tkanki statycznej (umieszczonego w miąższu śródmózgowia sąsiadującego z akweduktem) od obrazów fazy kodowanej prędkością, aby zminimalizować efekty prądów wirowych i składnika Maxwella. Dla każdego uczestnika wygenerowano następujące parametry: maksymalną prędkość skurczową (peak systolic velocity, PSV), maksymalną prędkość rozkurczową (peak diastolic velocity, PDV), objętość wyrzutową netto oraz krzywe przepływu. Zapisano powierzchnię przekroju poprzecznego akweduktu (Rysunek 2). Niezawodność międzyobserwacyjną i wewnątrzobserwacyjną oceniono za pomocą dwukierunkowych współczynników korelacji wewnątrzklasowej (ICC) z efektami losowymi dla zgodności absolutnej; w przypadku wszystkich parametrów wartości te przekroczyły 0,90.

Rycina 1: Akwizycja i przetwarzanie obrazów przepływu płynu mózgowo-rdzeniowego.(A) Strzałkowy obraz lokalizacyjny T1-zależny pokazujący ustawienie 2D cine w sekwencji kontrastu fazowego, ustawionego prostopadle do środka wodociągu mózgu (linia przerywana). (B) Obraz amplitudy sekwencji cine 2D pokazujący strukturę anatomiczną wodociągu mózgu. (C) Obraz fazy sekwencji cine 2D pokazujący intensywność sygnału przepływu CSF. Skróty: CA = wodociąg mózgu; CSF = płyn mózgowo-rdzeniowy. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Reprezentatywna krzywa prędkości przepływu płynu mózgowo-rdzeniowego. Krzywa zależności prędkości przepływu od czasu, przedstawiająca dwufazowy, quasi-sinusoidalny wzorzec przepływu CSF w ciągu jednego cyklu serca. Krzywa ilustruje zmiany prędkości podczas skurczu i rozkurczu. Skróty: CSF = płyn mózgowo-rdzeniowy. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Analiza parametrów związanych z przepływem w tętnicy szyjnej wewnętrznej
Dane przepływu czterowymiarowego zrekonstruowano i przeanalizowano przy użyciu modułu 4D Flow w specjalistycznym oprogramowaniu do postprocessingu. Korekcję przesunięcia fazowego zastosowano, wykorzystując dopasowanie wielomianu pierwszego stopnia do obszarów tkanek statycznych. Lewą tętnicę szyjną wewnętrzną (LICA) wyodrębniono za pomocą półautomatycznego narzędzia do segmentacji opartego na progu z ręcznym doprecyzowaniem, co pozwoliło na wygenerowanie angiograficznej objętości 3D. Automatycznie obliczono linię środkową naczynia, a następnie w segmencie C1 (około 2–3 cm powyżej rozwidlenia tętnicy szyjnej) wyznaczono trzy płaszczyzny przekroju poprzecznego ortogonalne do linii środkowej, zachowując 5-mm odstęp między sąsiednimi płaszczyznami. Do analizy wybrano lewą tętnicę szyjną wewnętrzną ze względu na jej spójny przebieg anatomiczny i charakterystykę przepływu; prawą tętnicą szyjną wewnętrzną nie zajmowano się, aby zminimalizować czas skanowania i obciążenie uczestników. Choć w tym badaniu eksploracyjnym, ze względów praktycznych, przyjęto podejście jednostronne, w przyszłych badaniach zasadna będzie ocena obustronna. Dla każdej płaszczyzny obrysowano światło naczynia za pomocą półautomatycznego algorytmu wykrywania krawędzi z ręczną korektą w każdej ramce czasowej cyklu serca (około 20 ramek na cykl) (Rysunek 3A–D).
Z konturów światła naczynia rozdzielonych w czasie obliczono następujące parametry: maksymalną prędkość przepływu, objętość przepływu w jednym cyklu serca oraz średnie naprężenie ścinające ścianek (WSS). WSS wyznaczono za pomocą następującego wzoru:
WSS = μ × (∂v/∂r)|_wall (1)
gdzie µ = 0.004 Pa·s jest przyjętą lepkością krwi, a ∂v/∂r to gradient prędkości prostopadły do ściany naczynia, uzyskany poprzez dopasowanie profilu parabolicznego do danych o prędkości w pobliżu ściany w każdym położeniu obwodowym i uśrednienie wzdłuż obwodu światła naczynia. Wygenerowano mapy prędkości przepływu, wektorów oraz WSS, a także wykreślono krzywe zależności przepływu i WSS od czasu. Wartości z trzech płaszczyzn uśredniono, aby uzyskać jedną reprezentatywną wartość dla każdego uczestnika. Cały proces przetwarzania obrazów został przeprowadzony niezależnie przez tych samych dwóch radiologów, a do analizy statystycznej wykorzystano średnią z ich pomiarów (Rysunek 4).

Rysunek 3: Rekonstrukcja trójwymiarowa i analiza hemodynamiczna lewej tętnicy szyjnej wewnętrznej. (A) Trójwymiarowa rekonstrukcja z renderowaniem objętościowym lewej tętnicy szyjnej po segmentacji naczyń. (B) Położenie trzech płaszczyzn analizy przekroju poprzecznego (Flow 1–3) wzdłuż segmentu C1. (C) Mapy prędkości przepływu. (D) Diagramy wektorowe przedstawiające kierunek przepływu. (E) Mapy naprężenia ścinającego ścianę. Skróty: LICA = lewa tętnica szyjna wewnętrzna; WSS = naprężenie ścinające ścianę. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 4: Krzywe postprocessingu lewej tętnicy szyjnej wewnętrznej. Krzywe zależności przepływu od czasu (panel lewy) oraz krzywe zależności naprężenia ścinającego ścianę od czasu (panel prawy) w obrębie jednego cyklu serca. Różne kolory reprezentują trzy płaszczyzny próbkowania (czerwony, żółty, zielony). Skróty: LICA = lewa tętnica szyjna wewnętrzna; WSS = naprężenie ścinające ścianę. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
Analiza statystyczna
Dane analizowano przy użyciu standardowego oprogramowania statystycznego (Tabela materiałów). Zmienne o rozkładzie normalnym, oceniane testem Shapiro–Wilka, przedstawiono jako średnia ± SD i porównano za pomocą jednoczynnikowej analizy wariancji (ANOVA). Zmienne o rozkładzie innym niż normalny przedstawiono jako mediana (IQR) i porównano testem Kruskala–Wallisa. Zmienne kategoryczne porównano testem chi-kwadrat. Do porównań międzygrupowych wykorzystano analizę kowariancji (ANCOVA) z wiekiem i wykształceniem jako kowariantami, stosując poprawkę Bonferroniego dla porównań wielokrotnych. Korelację Spearmana oraz wielokrotną regresję liniową wykorzystano do oceny zależności między obrazowaniem a miarami poznawczymi, z uwzględnieniem czynników ryzyka naczyniowego (nadciśnienie tętnicze, cukrzyca, hiperlipidemia). Przeprowadzono analizę charakterystyki operacyjnej odbiornika (ROC) w celu oceny wydajności diagnostycznej każdego parametru obrazowania. Obliczono pole pod krzywą (AUC), czułość oraz swoistość. Analizy ROC wykonano bez korekty o kowarianty, aby ocenić wewnętrzną zdolność dyskryminacyjną każdego parametru. Jednakże, ze względu na zaobserwowane różnice grupowe w wieku i wykształceniu, uzupełniające analizy z uwzględnieniem tych kowariantów przyniosły spójne wyniki. Za statystycznie istotną uznano wartość p < 0.05 (test obustronny). Pełna lista materiałów i sprzętu użytych w tym badaniu znajduje się w Tabeli materiałów.