Obrazy wentylacji
Jeśli procedury przygotowania zwierząt i wentylacji zostaną prawidłowo wdrożone, trójwymiarowe obrazowanie radialne może skutecznie uchwycić wzorce wentylacyjne, gdy akwizycja danych odbywa się podczas wdechowego lub wydechowego etapu cyklu (Rysunek 7). Chociaż obrazy te są gromadzone przez wiele oddechów, opisana tutaj metoda jest podobna do metody obrazowania pojedynczego oddechu stosowanej u ludzi. Wynika to z faktu, że wszystkie linie przestrzeni k są zbierane w konkretnym czasie cyklu oddechowego (np. wszystkie w końcowej fazie wdechu), co w sposób kumulatywny tworzy obraz reprezentujący pojedynczy wdech. Na obrazach tych główne drogi oddechowe (tchawica i oskrzela) są wyraźnie widoczne, nie występują żadne oczywiste artefakty obrazu, a miąższ płuc wydaje się w pełni rozwinięty z wysokim SNR (~15 lub wyższym, w zależności od polaryzacji początkowej). Choć może wydawać się to nieintuicyjne, obrazy uzyskane podczas wdechu (Rysunek 7A, SNR = 15) wykazują niższy SNR niż obrazy uzyskane podczas wydechu (Rysunek 7B, SNR = 27). Jest to związane z niewielkim ciśnieniem napędowym i powolnym przepływem stosowanym 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 koncentrację 129Xe (a tym samym sygnału) w większych drogach oddechowych. Efekt ten jest wizualnie potęgowany przez skalowanie przestrzeni k podczas rekonstrukcji obrazu, ponieważ wysoki sygnał z dróg oddechowych dominuje w zakresie dynamicznym danych. Niemniej jednak analiza powinna być prosta, jeśli miąższ zostanie poddany segmentacji w taki sposób, aby pominąć duże drogi oddechowe.

Rysunek 7: Reprezentatywne trójwymiarowe obrazy wentylacji radialnej u dorosłej myszy C57BL/6J, uzyskane przy użyciu opublikowanych wcześniej metod39. Choć oba obrazy charakteryzują się stosunkowo wysokim stosunkiem sygnału do szumu (SNR), obrazy przy (A) inspiracji mają niższy stosunek sygnału do szumu (SNR) niż obrazy przy (B) wygaśnięcia głównie dlatego, że 129Gaz Xe nie miał jeszcze czasu, aby przedostać się do płuc. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Ważne jest odróżnienie patologicznej niedrożności dróg oddechowych od innych przyczyn niskiego stosunku sygnału do szumu (SNR). Na Ryc. 7A nie przewiduje się niedrożności wentylacji (tj. niskiego VDP), ponieważ zwierzę było zdrową myszą typu dzikiego. W przypadku niedrożności wentylacji defekty zazwyczaj objawiałyby się jako obszary o niskim lub żadnym sygnale 129Xe. Niepatologiczne przyczyny niskiego sygnału obejmują zazwyczaj:
(i) Niska polaryzacja 129Xe: Powoduje ona niski SNR w miąższu płuc i drogach oddechowych bez towarzyszących artefaktów. Niską polaryzację można zazwyczaj przypisać nieprawidłowemu przeprowadzeniu procedur polaryzacji, depolaryzacji podczas transportu lub nadmiernej zawartości tlenu O2 w przewodach doprowadzających HP 129Xe (np. z powodu nieprawidłowo skalibrowanego respiratora).
(ii) Zbyt głębokie wprowadzenie kaniuli intubacyjnej: Jeśli kaniula zostanie wprowadzona zbyt głęboko do tchawicy, może ona przesunąć się w stronę lewego lub prawego oskrzela. W przypadku silnego zaklinowania, objętość oddechowa całego płuca może zostać ograniczona do zaledwie kilku płatów i dróg oddechowych, podczas gdy reszta płuca nie będzie otrzymywać powietrza. Jeśli tak się stanie, nadmiernie rozdęte płaty mogą ulec uszkodzeniu, co może wymagać eutanazji zwierzęcia. Alternatywnie, końcówka kaniuli może wejść w bezpośredni kontakt z tkanką. W takim przypadku 129Xe będzie w głównej mierze ograniczony do wnętrza kaniuli, co spowoduje uzyskanie silnego sygnału z widocznie prostej rurki i minimalnego sygnału w obrębie miąższu. W takiej sytuacji należy niezwłocznie skorygować położenie kaniuli i monitorować zwierzę pod kątem objawów stresu (np. nieprawidłowy, brak lub asynchroniczny ruch klatki piersiowej w stosunku do respiratora).
(iii) Zbyt płytkie wprowadzenie kaniuli intubacyjnej: Gdy bolus gazu jest podawany zbyt płytko w obrębie gardła, opór przepływu tego bolusa do płuc może być większy niż opór jego przepływu wokół zewnętrznej części kaniuli i wydostawania się przez usta. Może to skutkować niewystarczającym napowietrzeniem płuc. Stosunek sygnału do szumu (SNR) w tchawicy może być wysoki, podczas gdy SNR w miąższu płucnym będzie niski. Ze względu na przepływ gazu w górę tchawicy podczas akwizycji danych mogą wystąpić artefakty ruchowe. Charakter artefaktu będzie zależeć od rodzaju sekwencji (GRE, radialna itp.).
(iv) Nieprawidłowy moment akwizycji: Jeśli impulsy RF i akwizycja danych odbywały się w okresach wdechowego lub wydechowego przepływu gazów, zazwyczaj wystąpią silne artefakty ruchowe — szczególnie przy zastosowaniu sekwencji kartezjańskich — a SNR miąższu będzie niski.
(v) Oddychanie spontaniczne: Jeśli zwierzę jest niewystarczająco uspokojone podczas akwizycji danych, ruchy oddechowe spontaniczne będą nakładać się asynchronicznie na wzorzec wentylacji mechanicznej. Będzie to skutkować zmienną podażą 129Xe z oddechu na oddech, umiarkowanymi lub silnymi artefaktami ruchowymi oraz niskim stosunkiem sygnału do szumu (SNR) w miąższu. Kluczowe jest uważne monitorowanie dawki i czasu podania środków sedacyjnych.
(vi) Śmierć zwierzęcia: W rzadkich przypadkach myszy mogą zginąć (np. z powodu nieprawidłowego umiejscowienia kaniuli, niewystarczającej podaży O2, przedawkowania środków uspokajających lub ciężkiej patologii eksperymentalnej). Jeśli tak się stanie, tkanka płuca szybko sztywnieje, a miąższ płuc będzie słabo wypełniał się powietrzem. Górny odcinek tchawicy może pozostać jasny, ale sygnał z miąższu płuca będzie bardzo niski.
Spektroskopia wymiany gazowej
U ludzi parametry pochodne ze spektroskopii (np. intensywność sygnału, szerokość widma, przesunięcie chemiczne i faza) dostarczają informacji o dynamice krążeniowo-oddechowej i stanie chorobowym52,53. W przeciwieństwie do ludzi86, psów10 i szczurów20,101, widma całych płuc myszy otrzymywane, gdy namagnesowanie 129Xe jest ograniczone do miąższu płuca, wykazują pojedynczy, szeroki rezonans fazy rozpuszczonej, odpowiadający 129Xe rozpuszczonemu w tkance błony kapilarnej, osoczu krwi i erytrocytach54 (Rycina 8A, 197 ppm). Gdy 129Xe jest transportowany przez przepływ krwi do tkanek bardziej dystalnych, generuje odrębne częstotliwości rezonansowe z różnych tkanek, w tym z tkanki tłuszczowej54, mięśni, istoty szarej i białej102,103. Jednak ze względu na stosunkowo długie czasy przejścia (kilka sekund) i umiarkowaną rozpuszczalność104 (rozpuszczalność Ostwalda = 0,1 - 0,3), sygnały te są zazwyczaj słabe i obserwuje się je jedynie podczas dedykowanych, starannie przeprowadzonych eksperymentów. Wyjątkiem jest 129Xe rozpuszczony w tkankach tłuszczowych (Rycina 8A, 189 ppm). Z tych tkanek sygnał 129Xe może pojawić się w ciągu kilku sekund (np. podczas pierwszych kilku oddechów podczas przygotowania zwierzęcia). Ponadto intensywność sygnału może być wysoka, ponieważ rozpuszczalność Ostwalda 129Xe w tkankach tłuszczowych jest niemal 20-krotnie większa niż w osoczu krwi (tj. odpowiednio 1,7 vs 0,094 w 37 ˚C104). Przewiduje się również, że czasy relaksacji podłużnej będą stosunkowo długie (T1 in vivo wynosi ~20 s w lipidach105 w porównaniu do 3-8 s we krwi106). W rezultacie sygnał 129Xe z tkanki tłuszczowej może być stosunkowo wysoki w porównaniu z innymi tkankami.
Ze względu na odmienną dynamikę sygnału, eksperymenty z spektroskopii 129Xe mogą być dostosowane do badania tkanki tłuszczowej (np. w celu pomiaru termometrii tkanki tłuszczowej brunatnej107) lub wchłaniania gazu w płucach poprzez dobór odpowiednich parametrów akwizycji. Gdy widma 129Xe z płuc są rejestrowane przy ustawieniach opisanych w Tabeli 4 (tzn. TR = 50 ms i kąt flip 90˚ wycentrowany na piku fazy rozpuszczonej), obserwowane są piki zarówno fazy gazowej, jak i fazy rozpuszczonej (Rysunek 8A,B). Warto zauważyć, że objętość gazu w płucach myszy jest 2-3x większa od objętości tkanki, a rozpuszczalność 129Xe w tej tkance wynosi zaledwie ~10%104,108,109,10,1. Pomimo tego amplituda sygnału rozpuszczonego przekracza amplitudę fazy gazowej około 2-krotnie, ponieważ wzbudzenie fazy gazowej odbywa się daleko poza rezonansem (~15 kHz przy 7 T; Rysunek 8B, SNR odpowiednio 102 vs. 43). Zatem 20- do 30-krotnie większa pula magnetyzacji w fazie gazowej może być traktowana jako stabilny (tzn. niewyczerpany) rezerwuar dla sygnału z fazy rozpuszczonej podczas statycznych części cyklu oddechowego.
Przy nowym cyklu namagnesowania co ok. 60 ms i TR wynoszącym 50 ms, intensywność sygnału z fazy gazowej oraz fazy rozpuszczonej 129Xe oscyluje. W przypadku sygnału gazowego wynika to ze zmiany całkowitej objętości gazu (tj. wysokiej podczas wdechu i niskiej podczas wydechu). Sytuacja jest bardziej złożona dla sygnału rozpuszczonego, ponieważ namagnesowanie gazu źródłowego, T2* oraz objętość krwi w naczyniach włosowatych mogą fluktuować w zależności od stanu nadmuchania. W rezultacie stosunek amplitudy sygnału rozpuszczonego do gazowego również oscyluje w odniesieniu do cyklu oddechowego (Rysunek 8C, średni stosunek rozpuszczony:gaz = 1,9 ± 0,39).

Rycina 8: Reprezentatywna spektroskopia dynamiczna w ramach dorosła mysz szczepu C57BL/6J. (A) Widmo pojedynczego impulsu wykazuje sygnał pochodzący z 129Xe w fazie gazowej, rozpuszczony w tkance tłuszczowej (189 ppm) oraz rozpuszczony w błonie kapilarnej i tkankach krwi (197 ppm).B) Spektroskopia dynamiczna wykazująca zmienną w czasie intensywność sygnału zarówno dla formy rozpuszczonej, jak i gazowej 129Xe (czas powtórzenia = 50 ms). Aby ułatwić wizualizację, wyświetlono tylko pięć widm. Linie ukośne na osi czasu oznaczają przerwy w skali. (CStosunek amplitudy sygnału fazy rozpuszczonej do fazy gazowej w odstępie 2 s (odpowiadający czerwonej przerwie w skali na panelu B) wykazuje wyraźną zmienność w zależności od cyklu oddechowego (0,6 s na oddech). Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Poniżej znajduje się lista czynników, które mogą wpłynąć negatywnie na wyniki spektroskopii wymiany gazowej 129Xe. W każdym z tych przypadków prawdopodobnym skutkiem jest niski stosunek sygnału do szumu (SNR), co może sprawić, że dopasowanie widma i wynikowe parametry (np. stosunek amplitudy sygnału rozpuszczonego do gazowego, przesunięcie chemiczne, szerokość widma oraz faza) będą niedokładne.
(i) Błędna częstotliwość robocza: Ze względu na znacznie większą objętość gazu i niską rozpuszczalność w tkankach, impulsy RF wycentrowane na częstotliwości fazy gazowej nie wywołają sygnału z fazy rozpuszczonej. Impulsy RF powinny być wycentrowane na częstotliwości fazy rozpuszczonej lub w jej pobliżu.
(ii) Zbyt szerokie pasmo: Pasmo wzbudzenia powinno być wystarczająco szerokie, aby wzbudzić część 129Xe w fazie gazowej w celu pomiaru stosunku amplitudy sygnału fazy rozpuszczonej do gazowej. Jeśli jednak pasmo jest zbyt szerokie, impuls wzbudzenia zniszczy magnetyzację w fazie gazowej, zanim ta rozpuści się w tkance. Przejawia się to jako bardzo wysoki sygnał gazowy oraz niski lub brak sygnału z fazy rozpuszczonej, pomimo że impuls został wycentrowany na częstotliwości fazy rozpuszczonej.
(iii) Nieoptymalizowany TR: Aby zbadać dynamikę sygnału w trakcie jednego oddechu, twierdzenie Nyquista-Shannona o próbkowaniu wymaga, aby sygnał był próbkowany z częstotliwością co najmniej dwukrotnie większą niż jego częstotliwość (tutaj 10 oddechów na minutę próbkowano co 50 ms, czyli ~12x na oddech). W zależności od warunków eksperymentalnych możliwe jest dalsze badanie dynamiki sygnału kardiopulmonarnego przy szybszej częstotliwości próbkowania. Wyniki tych eksperymentów zależałyby od wielu czynników, w tym od szerokości pasma wzbudzenia, częstości oddychania i tętna. Należy również wziąć pod uwagę czas wymiany ksenonu (tzn. czas rozpuszczania się ksenonu w tkance; ~56 ms54) oraz czas przejścia płucnego (tzn. czas przepływu krwi z prawej do lewej komory; ~830 ms przy tętnie ~436 uderzeń/min112). Zbyt częste próbkowanie doprowadzi do przedwczesnego zniszczenia namagnesowania i niskiego poziomu sygnału, ale jednocześnie zmniejszy wpływ przepływu krwi. Rzadsze próbkowanie pozwoli na dalsze rozpuszczanie namagnesowania i uzyskanie silniejszego sygnału, jednak stanie się ono bardziej zależne od przepływu krwi i mniej zlokalizowane w obszarach wymiany gazowej w pęcherzykach płucnych.
Obrazy dyfuzyjne
Obrazy ważone dyfuzją, stosowane do ilościowego określania wielkości przestrzeni powietrznych pęcherzyków płucnych, zostały rygorystycznie zwalidowane w modelach małych zwierząt za pomocą tradycyjnych metod histomorfometrii98,13,14. W niniejszej pracy ograniczyliśmy dyskusję do obrazowania 2D z selekcją warstw, aby zminimalizować czas skanowania i zużycie gazu, zachowując zgodność z wcześniej zwalidowanymi protokołami98. Do obliczenia ADC wystarczające są co najmniej 2 obrazy o różnych wartościach b, natomiast do zaawansowanych obliczeń morfometrycznych (takich jak średnia odległość między ściankami pęcherzyków, tzw. mean linear intercept) niezbędne są 4 lub więcej obrazy o różnych wartościach b. Pozyskanie wielu obrazów o różnych wartościach b oferuje znaczącą zaletę w postaci zwiększenia dokładności obliczeń. W przeciwieństwie do obrazowania ludzi podczas pojedynczego wstrzymania oddechu, obrazowanie przedkliniczne obejmujące wiele oddechów zapewnia stale odnawiany sygnał. Dzięki temu możliwe jest pozyskanie wielu obrazów o wysokim SNR dla różnych wartości b, ograniczonych jedynie objętością dostępnego gazu HP oraz czasem sedacji zwierzęcia. W tym przypadku pozyskano 7 obrazów o różnych wartościach b, każdy z SNR >15, wykorzystując 400 mL gazu 129Xe w czasie skanowania wynoszącym 18 min (Rysunek 9). Sygnał w każdym obrazie zanika wraz ze wzrostem natężenia gradientu dyfuzji (wartości b). Ten zanik sygnału jest proporcjonalny do ADC 129Xe w płucach. Duże drogi oddechowe można odróżnić od tkanki pęcherzykowej, porównując obrazy b0 i b6. W obszarach, które anatomicznie prawdopodobnie zawierają drogi oddechowe, piksele będą bardzo jasne na obrazie b0, ale ciemne na obrazie b6 (czerwone strzałki, Rysunek 9). Piksele te powinny zostać wykluczone z analizy miąższu płuc.

Rysunek 9Reprezentatywne obrazy obrazowania rozproszonego (DWI) u dorosłej myszy szczepu C57BL/6J. Duże drogi oddechowe (czerwone strzałki), przez które 129Xe dyfunduje swobodnie podczas rozcieńczenia tlenem, wykazując zwiększone tłumienie sygnału w stosunku do miąższu płuc. Skutkuje to jasnymi drogami oddechowymi w b0 obrazowanie i drogi oddechowe o niskim sygnale w b6 obraz. Myszy i szczury mają podobnie niskie wartości ADC w stosunku do ludzi, jednak wartości te będą zależeć od czynników biologicznych i eksperymentalnych, takich jak wiek, objętość tidalna, ciśnienie i parametry dyfuzyjne. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
Współczynnik dyfuzji pozornej (ADC) tej zdrowej, dorosłej myszy typu dzikiego jest zgodny z oczekiwanym wynikiem na podstawie literatury, wieku zwierzęcia, stosunku sygnału do szumu (SNR), parametrów akwizycji i innych czynników. Czynniki, które mogą negatywnie wpływać na precyzję i dokładność wyników ADC, obejmują:
(i) Niski stosunek sygnału do szumu (SNR): Ważne jest, że niski stosunek sygnału do szumu (SNR) w obrazach powoduje przesunięcie obliczeń ADC w stronę wartości niższych.99,115Z tego powodu wykluczamy voksele o SVNR0 < 2,5-krotność szumu obrazu97.
(ii) Nieprawidłowa segmentacja: Niewłaściwa segmentacja (np. stosowanie progowania sygnału bez wizualnej kontroli jakości) może zniekształcić wyniki. Jeśli obrazy zawierają artefakty (np. spowodowane ruchem), mogą one przekroczyć próg i zostać błędnie zakwalifikowane jako miąższ. Ponadto szczególnie ważne jest wykluczenie pikseli w pobliżu dużych dróg oddechowych, ponieważ efekty objętości cząstkowej mogą przesunąć wartości ADC w stronę wartości wyższych.
(iii) Źle skalibrowany respirator: Ogólnie rzecz biorąc, rozcieńczenie 129Xe lżejszymi gazami (np. O2 lub N2) może przesunąć ADC w stronę wyższych wartości, dlatego w celu porównań ilościowych niezbędna jest wentylacja ze znanymi mieszankami gazowymi. W przypadku wentylacji mieszaniną 129Xe/O2, zwiększenie stężenia tlenu przed podaniem zwiększy relaksację T1 i zmniejszy namagnesowanie 129Xe w przewodach doprowadzających gaz, co obniży SNR.
(iv) Derekrutacja płuc: U zwierząt zanieczulonych i wentylowanych mechanicznie z czasem może dojść do zapadnięcia się pęcherzyków płucnych (miażdrycy). Gdy zjawisko to występuje podczas podczas wentylacji o stałej objętości ta sama objętość oddechowa jest nieuchronnie redystrybuowana do pęcherzyków płucnych, które pozostają otwarte. W miarę rozwoju umiarkowanej lub ciężkiej atelektazy pozostała część płuca może ulec nadmiernemu rozdęciu, co prowadzi do uszkodzeń i zwiększenia wartości ADC w tych obszarach. Regularne manewry rekrutacyjne oraz zastosowanie PEEP zapobiegają zapadaniu się pęcherzyków płucnych, przeciwdziałają zmniejszeniu podatności płuc i minimalizują ryzyko uszkodzenia płuc, które mogłoby zniekształcić wyniki ADC.
(v) Słabe dopasowanie danych: Należy wziąć pod uwagę metodę szacowania ADC. Dopasowanie log-liniowe jest najbardziej wydajne obliczeniowo, ale wprowadza błąd systematyczny w obrazach o niskim SNR. Dlatego najlepiej stosować je w obrazach o wysokim SNR (>20). Metoda liniowa ważona zmniejsza ten błąd i jest odpowiednia dla obrazów z SNR > 15. Estymacja bayesowska może dostarczać wartości ADC z niską niepewnością nawet przy SNR <10, ale jest kosztowna obliczeniowo9.
(vi) Nieoptymalny czas dyfuzji: Optymalny czas dyfuzji zależy od maksymalnej siły gradientu skanera MRI oraz promienia dróg oddechowych w obrębie acinusów u myszy (promień ~10 µm49), który jest około trzy razy mniejszy niż u ludzi. Aby zapewnić optymalne warunki, czas dyfuzji ustawia się tak, aby osiągnąć długość dyfuzji większą niż średni promień pęcherzyka płucnego, ale mniejszą niż średnia długość przewodów pęcherzykowych. Optymalny czas dyfuzji dla maksymalnej siły gradientu w danej lokalizacji można wyznaczyć stosując metody opisane przez Sukstanskiego i Yablonskiego49. Zastosowanie nieoptymalnego czasu dyfuzji może prowadzić do nieliniowych spadków ADC 129Xe, co komplikuje 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 różnych wartościach b. Jeśli poświęci się CNR, stosując wyłącznie niskie wartości b (tj. spadek sygnału dyfuzyjnego, jak w przypadku pozyskania tylko pierwszych 3 obrazów o różnych wartościach b na Rysunku 9), oszacowanie ADC będzie niedokładne ze względu na słabe dopasowanie. Jeśli ważenie dyfuzyjne jest zbyt silne, wynikowy SNR będzie niewystarczający, a wartość ADC zostanie przesunięta w stronę niższych wartości.