$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Obrazy wentylacji
Jeśli procedury przygotowania i wentylacji zwierząt są prawidłowo wdrożone, obrazowanie radialne 3D może z powodzeniem uchwycić wzorce wentylacji, gdy akwizycja danych odbywa się zarówno przy wdechu, jak i wydechu (Rysunek 7). Chociaż obrazy te są zbierane przez wiele oddechów, opisana tutaj metoda jest podobna do metody obrazowania pojedynczego oddechu stosowanej u ludzi. Dzieje się tak dlatego, że wszystkie linie k-przestrzeni są zbierane w określonym momencie cyklu oddechu (np. wszystkie przy końcowym wdechu), kumulatywnie tworząc obraz reprezentatywny dla pojedynczego wdechu. Na tych obrazach główne drogi oddechowe (tchawica i oskrzela) są wyraźnie widoczne, nie ma żadnych oczywistych artefaktów obrazu, a miąższ płuc wydaje się w pełni napompowany z wysokim SNR (~15 lub wyższym, w zależności od początkowej polaryzacji). Chociaż jest to sprzeczne z intuicją, obrazy uzyskane podczas wdechu (Rysunek 7A, SNR = 15) wykazują niższy SNR niż te uzyskane w momencie wygaśnięcia (Rysunek 7B, SNR = 27). Jest to związane ze skromnym ciśnieniem napędowym i powolnym przepływem używanym podczas wentylacji tej myszy. W przypadku obrazów wdechowych, gaz HP nie został jeszcze w pełni wypłukany do miąższu płuc, co spowodowało, że 129Xe (a tym samym sygnał) był skoncentrowany w większych drogach oddechowych. Efekt ten jest wizualnie wzmocniony przez skalowanie k-space podczas rekonstrukcji obrazu, ponieważ wysoki sygnał dróg oddechowych dominuje w zakresie dynamicznym danych. Niemniej jednak analiza powinna być prosta, jeśli miąższ jest podzielony na segmenty w taki sposób, że duże drogi oddechowe są pominięte.

Rysunek 7: Reprezentatywne obrazy wentylacji radialnej 3D u dorosłej myszy C57BL/6J uzyskane przy użyciu wcześniej opublikowanych metod39. Podczas gdy oba obrazy mają stosunkowo wysoki SNR, obrazy przy wdechu (A) mają niższy SNR niż obrazy przy wydechu (B), głównie dlatego, że gaz 129Xe nie zdążył jeszcze przedostać się do płuc. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Ważne jest, aby odróżnić patologicznie utrudnioną wentylację od innych przyczyn niskiego SNR. W Rysunek 7A, nie przewiduje się żadnych przeszkód wentylacyjnych (tj. niskiego VDP), ponieważ zwierzę było zdrową, dziką myszą. Jeśli wentylacja jest zakłócona, defekty zwykle pojawiają się jako obszary z niskim lub zerowym sygnałem 129Xe. Niepatologiczne przyczyny niskiego sygnału często obejmują:
(i) Niska polaryzacja 129Xe: Powoduje to niski SNR w miąższu płuc i drogach oddechowych bez towarzyszących artefaktów. Niska polaryzacja może być zwykle przypisana nieprawidłowemu przestrzeganiu procedur polaryzacji, depolaryzacji podczas transportu lub nadmiernej zawartości gazuO2 w przewodzie doprowadzającym HP 129Xe (np. z nieprawidłowo skalibrowanego wentylatora).
(ii) Kaniula intubacyjna włożona zbyt głęboko: Jeśli zostanie włożona zbyt głęboko do tchawicy, kaniula może odchylić się w dół lewego lub prawego oskrzela. Jeśli jest wystarczająco ciasno osadzony, objętość oddechowa całego płuca może być ograniczona tylko do kilku płatów i dróg oddechowych, podczas gdy reszta płuca nie otrzymuje żadnego. Jeśli tak się stanie, nadmiernie napompowane płaty mogą zostać uszkodzone i może być konieczne poddanie zwierzęcia eutanazji. Alternatywnie, końcówka kaniuli może mieć bezpośredni kontakt z tkanką. W tym przypadku 129Xe będzie przede wszystkim ograniczony do wnętrza kaniuli, co spowoduje jasny sygnał z widocznie prostej rurki i minimalny sygnał w miąższu. W takim przypadku należy niezwłocznie skorygować pozycję kaniuli, a zwierzę powinno być monitorowane pod kątem oznak niepokoju (np. nieprawidłowych, nieobecnych lub asynchronicznych ruchów klatki piersiowej w porównaniu z respiratorem).
(iii) Kaniula intubacyjna umieszczona zbyt płytko: Gdy bolus gazu jest dostarczany zbyt płytko do jamy ustnej gardła, opór dla tego bolusa wpływającego do płuc może być większy niż opór dla tego bolusa przepływającego po zewnętrznej stronie kaniuli i na zewnątrz ust. Może to spowodować słabe napompowanie płuc. Stosunek sygnału do ssawicy może być wysoki, podczas gdy SNR w miąższu będzie niski. Artefakty ruchu mogą być obecne z powodu gazu przepływającego w górę tchawicy podczas pozyskiwania danych. Charakter artefaktu będzie się różnił w zależności od typu sekwencji (GRE, radialna itp.).
(iv) Nieprawidłowy czas akwizycji: Jeśli pulsowanie RF i akwizycja danych wystąpiły w okresach wdechowego lub wydechowego przepływu gazu, zwykle występują poważne artefakty ruchu - szczególnie podczas korzystania z sekwencji kartezjańskich - a miąższowy SNR będzie niski.
(v) Spontaniczne oddychanie: Jeśli zwierzę nie jest wystarczająco uspokojone podczas zbierania danych, spontaniczny ruch oddechowy zostanie nałożony asynchronicznie na wzorzec wentylacji mechanicznej. Spowoduje to zmienny przepływ 129Xe z oddechu do oddechu, umiarkowane do ciężkich artefakty ruchowe i słaby SNR miąższu. Uważne śledzenie dawki i czasu stosowania środków uspokajających ma kluczowe znaczenie.
(vi) Śmierć zwierzęcia: W rzadkich przypadkach myszy mogą umrzeć (np. z powodu nieprawidłowego ułożenia kaniuli, niewystarczającego dostarczeniaO2, przedawkowania środków uspokajających lub ciężkiej patologii eksperymentalnej). Jeśli tak się stanie, tkanka płucna szybko sztywnieje, a miąższ płuc słabo się napompuje. Górna tchawica może pozostać jasna, ale miąższ płuc będzie miał bardzo słaby sygnał.
Spektroskopia wymiany gazowej
U ludzi, parametry pochodzące ze spektroskopii (np. intensywność sygnału, szerokość spektralna, przesunięcie chemiczne i faza) dają wgląd w dynamikę krążeniowo-oddechową i stan choroby52,53. W przeciwieństwie do ludzi86, dogs100 i rats20,101, widma całych płuc od myszy uzyskane, gdy namagnesowanie 129Xe jest ograniczone do miąższu płuca, wykazują pojedynczy, szeroki, rozpuszczony rezonans fazowy odpowiadający 129Xe rozpuszczonemu w tkance błony naczyń włosowatych, osoczu krwi i krwinkach czerwonych54 (Rysunek 8A, 197 ppm). Po przetransportowaniu przez przepływ krwi do bardziej dystalnych tkanek, 129Xe wytwarza wyraźne częstotliwości rezonansowe z różnych tkanek, w tym tkanki tłuszczowej54, mięśnie, istota szara i istota biała102,103. Jednak ze względu na stosunkowo długie czasy tranzytu (kilka sekund) i skromną rozpuszczalność104 (rozpuszczalność Ostwalda = 0,1 - 0,3), sygnały te są zazwyczaj słabe i obserwowane tylko podczas dedykowanych, ostrożnych eksperymentów. Wyjątkiem jest 129Xe rozpuszczony w tkankach tłuszczowych (Rysunek 8A, 189 ppm). Z tych tkanek może pojawić się sygnał 129Xe w ciągu kilku sekund (np. w ciągu pierwszych kilku oddechów podczas konfiguracji zwierzęcia). Co więcej, intensywność sygnału może być wysoka, ponieważ rozpuszczalność Ostwalda 129Xe jest prawie 20-krotnie większa w tkance tłuszczowej niż w osoczu krwi (tj. odpowiednio 1,7 vs 0,094 w temperaturze 37 °C104). Oczekuje się również, że podłużne czasy relaksacji będą stosunkowo długie (in vivo T1 wynosi ~20 s w lipidach105 vs. 3-8 s w blood106). W rezultacie sygnał 129Xe z tkanki tłuszczowej może być stosunkowo wysoki w porównaniu z innymi tkankami.
Ze względu na różną dynamikę sygnału, eksperymenty spektroskopii 129Xe mogą być dostosowane do badania tkanki tłuszczowej (np. do celu termometrii brunatnej tkanki tłuszczowej107) lub poboru gazu w płucach poprzez wybór odpowiednich parametrów akwizycji. Gdy widma 129Xe z płuc są rejestrowane przy ustawieniach opisanych w Tabeli 4 (tj. TR = 50 ms i kąt obrotu 90° wyśrodkowany na piku fazy rozpuszczonej), obserwuje się zarówno piki fazy gazowej, jak i fazy rozpuszczonej (Rysunek 8A,B). Warto zauważyć, że objętość gazu w płucu myszy jest 2-3 razy większa niż objętość tkanki, a rozpuszczalność 129Xe w tej tkance wynosi tylko ~10%104,108,109,110,111. Mimo to amplituda rozpuszczonego sygnału jest około 2-krotnie większa niż w fazie gazowej, ponieważ wzbudzenie fazy gazowej jest dalekie od rezonansu (~15 kHz przy 7 T; Rysunek 8B, SNR = odpowiednio 102 vs. 43). W związku z tym 20-30-krotnie większą pulę namagnesowania w fazie gazowej można przybliżyć jako stabilny (tj. niewyczerpany) zbiornik dla sygnału fazy rozpuszczania podczas statycznych części cyklu oddechowego.
Z nowym oddechem namagnesowania co ~600 ms i TR 50 ms, intensywność sygnału zarówno z fazy gazowej, jak i rozpuszczonej 129Xe oscyluje. W przypadku sygnału gazowego można to przypisać zmianie całkowitej objętości gazu (tj. wysoka przy wdechu i niska przy wydechu). Sytuacja jest bardziej złożona w przypadku rozpuszczonego sygnału, ponieważ namagnesowanie gazu źródłowego,T2* i objętość krwi włośniczkowej mogą się wahać wraz ze stanem nadmuchiwania. W rezultacie, stosunek amplitudy sygnału rozpuszczonego do gazu oscyluje również w odniesieniu do cyklu oddechowego (Rysunek 8C, średni stosunek rozpuszczonego do gazu = 1,9 ± 0,39).

Rysunek 8: Reprezentatywna spektroskopia dynamiczna u dorosłej myszy C57BL/6J. (A) Pojedyncze widmo impulsów pokazuje sygnał pochodzący od 129Xe w fazie gazowej, rozpuszczony w tkance tłuszczowej (189 ppm) i rozpuszczony zarówno w błonie włosowatej, jak i tkankach krwi (197 ppm). (B) Spektroskopia dynamiczna pokazująca zmienne w czasie natężenie sygnału zarówno rozpuszczonego, jak i gazowego 129Xe (czas powtarzania = 50 ms). Aby ułatwić wizualizację, wyświetlanych jest tylko pięć widm. Ukośne linie na osi czasu reprezentują podziały skali. (C) Stosunek amplitudy sygnału rozpuszczonego do gazowego w przedziale 2 s (odpowiadający czerwonej podziałce skali w panelu B wykazuje wyraźną zmienność w cyklu oddechowym (0,6 s na oddech). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Poniżej znajduje się lista czynników, które mogą zagrozić wynikom spektroskopii wymiany gazowej 129Xe. W każdym przypadku prawdopodobnym wynikiem jest niski SNR, co może spowodować, że dopasowanie spektralne i wynikowe wskaźniki (np. stosunek amplitudy sygnału rozpuszczonego do gazu, przesunięcie chemiczne, szerokość widma i faza) będą niedokładne.
(i) Niewłaściwa częstotliwość pracy: Ze względu na znacznie większą objętość gazu i niską rozpuszczalność w tkankach, impulsy RF wyśrodkowane na częstotliwości fazy gazowej nie wytwarzają sygnału fazy rozpuszczonej. Impulsy RF powinny być wyśrodkowane na lub w pobliżu częstotliwości fazy rozpuszczonej.
(ii) Zbyt szerokie pasmo: Pasmo wzbudzenia powinno być na tyle szerokie, aby wzbudzić część fazy gazowej 129Xe, aby zmierzyć stosunek amplitudy sygnału rozpuszczonego do gazowego. Ale jeśli szerokość pasma jest zbyt szeroka, impuls wzbudzenia zniszczy namagnesowanie w fazie gazowej, zanim rozpuści się w tkance. Może to być postrzegane jako bardzo wysoki sygnał gazowy i niski lub zerowy sygnał fazy rozpuszczonej, mimo że impuls jest wyśrodkowany na częstotliwości fazy rozpuszczonej.
(iii) Niezoptymalizowany TR: Aby zbadać dynamikę sygnału podczas oddechu, twierdzenie Nyquista-Shannona o próbkowaniu wymaga, aby sygnał był próbkowany z co najmniej dwukrotnie większą częstotliwością (tutaj próbkowano 100 oddechów na minutę co 50 ms, czyli ~12x na oddech). W zależności od warunków eksperymentalnych możliwe jest dalsze badanie dynamiki sygnału krążeniowo-oddechowego przy szybszej częstotliwości próbkowania. Wyniki tych eksperymentów zależałyby od wielu czynników, w tym szerokości pasma wzbudzenia, częstości oddechów i tętna. Należy również wziąć pod uwagę czas wymiany ksenonu (tj. czas rozpuszczenia się ksenonu w tkance; ~56 ms54) i czas przejścia przez płuca (tj. czas przepływu krwi z prawej do lewej komory; ~830 ms dla tętna ~436 uderzeń/min112). Zbyt częste pobieranie próbek spowoduje wczesne zniszczenie namagnesowania i słaby sygnał, ale zmniejszy również wpływ przepływu krwi. Rzadsze pobieranie próbek pozwoli na dalsze rozpuszczenie namagnesowania i wytworzenie większej ilości sygnału, ale stanie się ono bardziej zależne od przepływu krwi i mniej zlokalizowane w obszarach wymiany gazowej w pęcherzykach płucnych.
Obrazy ważone dyfuzją
Obrazy ważone dyfuzją, używane do ilościowego określania wielkości pęcherzyków płucnych i przestrzeni powietrznej, zostały rygorystycznie zweryfikowane na małych modelach zwierzęcych za pomocą tradycyjnych metod histomorfometrii98,113,114. W tym miejscu ograniczyliśmy dyskusję do obrazowania 2D, selektywnego obrazowania warstwowego, aby zminimalizować czas skanowania i zużycie gazu, jednocześnie dopasowując wcześniej zatwierdzone protokoły98. Do obliczenia ADC wystarczą co najmniej 2 obrazy o wartości b, a do zaawansowanych obliczeń morfometrycznych (takich jak średnia wartość liniowa) niezbędne są 4 lub więcej obrazów o wartości b. Akwizycja wielu obrazów o wartości b oferuje znaczącą zaletę poprzez zwiększenie dokładności obliczeń. W przeciwieństwie do obrazowania ludzi podczas pojedynczego wstrzymania oddechu, obrazowanie przedkliniczne na wielu oddechach zapewnia stale uzupełniany sygnał. Z tego powodu możliwe jest uzyskanie wielu obrazów o wysokiej wartości SNR b, ograniczonych jedynie ilością dostępnego gazu HP i czasem sedacji zwierząt. W tym przypadku uzyskaliśmy 7 obrazów o wartości b, każdy z SNR >15, przy użyciu 400 ml gazu 129Xe w ciągu 18 minut skanowania (Rysunek 9). Sygnał na każdym obrazie zanika wraz ze wzrostem siły gradientu dyfuzji (wartość b). Ten zanik sygnału jest proporcjonalny do ADC 129Xe w płucach. Duże drogi oddechowe można odróżnić od tkanki groniastej, porównując obrazy b0 i b6. W regionach, w których anatomicznie istnieje prawdopodobieństwo, że zawierają drogi oddechowe, piksele będą bardzo jasne na obrazie b0, ale będą ciemne na obrazie b6 (czerwone strzałki, Rysunek 9). Piksele te powinny być wyłączone z analizy miąższu płuc.

Rysunek 9: Reprezentatywne obrazy ważone dyfuzją w dorosłej myszy C57BL/6J. Duże drogi oddechowe (czerwone strzałki), przez które 129Xe swobodnie dyfunduje po rozcieńczeniu tlenem, wykazują zwiększone tłumienie sygnału w stosunku do miąższu płuc. Powoduje to jasne drogi oddechowe na obrazie b0 i drogi oddechowe o słabym sygnale na obrazie b6. Myszy i szczury mają podobnie małe wartości ADC w porównaniu z ludźmi, ale wartości będą zależeć od czynników biologicznych i eksperymentalnych, takich jak wiek, objętość oddechowa, ciśnienie i parametry dyfuzji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
ADC tej zdrowej, dorosłej myszy typu dzikiego podąża za oczekiwanym wynikiem na podstawie literatury, jej wieku, SNR, parametrów akwizycji i innych. Czynniki, które mogą mieć szkodliwy wpływ na precyzję i dokładność wyników ADC, obejmują:
(i) Low SNR: Co ważne, niski SNR na obrazach skłania obliczenia ADC do niższych wartości99,115. Z tego powodu wykluczamy woksele o SVNR0 < 2,5-krotności szumu obrazu97.
(ii) Niewłaściwa segmentacja: Słaba segmentacja (np. użycie progowania sygnału bez wizualnej kontroli jakości) może zniekształcić wyniki. Jeśli obrazy zawierają artefakty (np. te spowodowane ruchem), mogą one przekroczyć próg i zostać błędnie sklasyfikowane jako miąższ. Co więcej, szczególnie ważne jest wykluczenie pikseli w pobliżu dużych dróg oddechowych, ponieważ częściowe efekty głośności mogą skłaniać ich ADC do wyższych wartości.
(iii) Źle skalibrowany wentylator: Ogólnie rzecz biorąc, rozcieńczenie 129Xe lżejszymi gazami (np. O2 lub N2) może przesunąć ADC w kierunku wyższych wartości, więc wentylacja za pomocą znanych mieszanin gazów jest konieczna do porównań ilościowych. W przypadku wentylacji za pomocą mieszaniny 129Xe/O2 wzrost stężenia tlenu przed podaniem zwiększy relaksację T1 i zmniejszy namagnesowanie 129Xe w rurze doprowadzającej gaz, zmniejszając w ten sposób SNR.
(iv) Derekrutacja płuc: U zwierząt znieczulonych i wentylowanych mechanicznie, z czasem może dojść do zapadnięcia się pęcherzyków płucnych (niedodmy). Kiedy dzieje się to podczas wentylacji o stałej objętości, ta sama objętość oddechowa jest z konieczności redystrybuowana do pęcherzyków płucnych, które pozostają rekrutowane. W miarę rozwoju umiarkowanej do ciężkiej niedodmy reszta płuc może stać się nadmiernie napompowana, powodując uraz i zwiększony ADC w tych regionach. Regularne manewry rekrutacyjne i PEEP zapobiegają zapadaniu się pęcherzyków płucnych, pozwalają uniknąć zmniejszonej podatności płuc i minimalizują ryzyko uszkodzenia płuc, które może zniekształcić wyniki ADC.
(v) Słabe dopasowanie danych: Metoda szacowania ADC musi być brana pod uwagę. Dopasowanie logarytmiczno-liniowe jest najbardziej wydajne obliczeniowo, ale wprowadza odchylenie w obrazach o niskim SNR. Dlatego najlepiej nadaje się do obrazów o wysokim stosunku do ceny (SNR) (>20). Metoda ważonej linii zmniejsza to odchylenie i jest odpowiednia dla obrazów o SNR > 15. Estymacja bayesowska może dać oszacowania ADC z niską niepewnością nawet przy SNR <10, ale jest kosztowna obliczeniowo99.
(vi) Nieoptymalny czas dyfuzji: Na optymalny czas dyfuzji wpływa maksymalna siła gradientu skanera MRI oraz promień ziarniastych dróg oddechowych u myszy (promień ~100 μm49), który jest w przybliżeniu trzy razy mniejszy niż u ludzi. Aby zapewnić optymalne warunki, czas dyfuzji ustala się tak, aby osiągnąć długość dyfuzji większą niż średni promień wyrostka zębodołowego, ale mniejszą niż średnia długość przewodów pęcherzykowych. Optymalny czas dyfuzji dla maksymalnej siły gradientu specyficznego dla danego miejsca można określić za pomocą metod opisanych przez Sukstanskii i Yablonskiy49. Użycie nieoptymalnego czasu dyfuzji może prowadzić do nieliniowych spadków 129Xe ADC, komplikując interpretację wyników.
(vii) Suboptymalne wartości b: Istnieje kompromis między maksymalizacją SNR a maksymalizacją stosunku kontrastu do szumu (CNR) generowanego na obrazach o wartości b. Jeśli CNR jest poświęcany przez użycie tylko niskich wartości b (tj. dyfuzyjny zanik sygnału, tak jakby uzyskiwano tylko pierwsze 3 obrazy o wartości b w Rysunek 9), to oszacowanie ADC będzie niedokładne z powodu złego dopasowania. Jeśli waga dyfuzyjna jest zbyt silna, wynikowy SNR będzie niewystarczający, a ADC będzie skłaniał się ku niższym wartościom.