MP-MRI u dzieci z mukowiscydozą: nowa analiza rozdziela defekty od obszarów niskiej funkcji płuca

Funkcjonalny rezonans magnetyczny płuc wykonywany podczas swobodnego oddychania może jednocześnie dostarczać map regionalnej wentylacji i perfuzji bez podawania środka kontrastowego oraz bez stosowania wziewnych gazów znacznikowych. Badacze związani z ośrodkami w Bernie, Bazylei i Lucernie zaadaptowali metodę linear binning do obrazowania MP-MRI u dzieci. W analizie obejmującej zdrowe dzieci oraz dzieci z mukowiscydozą metoda pozwoliła nie tylko różnicować obie grupy, ale przede wszystkim oddzielać obszary odpowiadające najgłębszym defektom funkcji od regionów o mniej nasilonym zmniejszeniu wentylacji lub perfuzji. Praca została opublikowana w Magnetic Resonance in Medicine.

MP-MRI umożliwia funkcjonalną ocenę płuc podczas swobodnego oddychania

Autorami pracy są Marion Curdy, Carmen Streibel, Oliver Bieri, Elisabeth Kieninger, Florian Wyler, Lukas Ebner, Grzegorz Bauman i Philipp Latzin. Badacze byli związani m.in. z Pediatric Respiratory Medicine, Department of Pediatrics, Inselspital, Bern University Hospital, University of Bern; Graduate School for Health Sciences, University of Bern; Division of Radiological Physics, Department of Radiology, University of Basel Hospital; Department of Biomedical Engineering, University of Basel; Department of Radiology and Nuclear Medicine, Lucerne Cantonal Hospital, University of Lucerne oraz Department of Diagnostic, Interventional and Paediatric Radiology, Inselspital, Bern University Hospital, University of Bern.

Publikacja ukazała się online 13 sierpnia 2026 roku w czasopiśmie Magnetic Resonance in Medicine. Celem autorów było dostosowanie i ocena metody linear binning, wcześniej opracowanej dla rezonansu magnetycznego z hiperpolaryzowanym ksenonem-129, do funkcjonalnego protonowego MRI płuc wykorzystującego matrix pencil decomposition.

Funkcjonalne techniki MRI takie jak SENCEFUL, PREFUL czy MP-MRI mają szczególne znaczenie w pediatrii, ponieważ mogą być wykonywane podczas swobodnego oddychania i nie wymagają podawania środka kontrastowego ani wziewnych gazów znacznikowych. Dodatkową zaletą tych technik jest możliwość uzyskiwania map ważonych zarówno wentylacją, jak i perfuzją.

MP-MRI jest pochodną metod wykorzystujących dekompozycję Fouriera, przy czym matrix pencil decomposition ma być bardziej odporna na artefakty wynikające z krótkich szeregów czasowych. Wcześniejsze badania cytowane przez autorów wykazały powtarzalność MP-MRI u zdrowych dzieci oraz dzieci z mukowiscydozą, a parametry wentylacji i perfuzji uzyskiwane tą metodą korelowały z wynikami badań czynności płuc w wybranych populacjach pediatrycznych.

Klasyczny próg 70% mediany ma istotne ograniczenia

Do standardowych parametrów funkcjonalnego MRI płuc należą odsetek defektu wentylacji oraz odsetek defektu perfuzji. Ponieważ MRI nie zapewnia bezwzględnej skali sygnału, jedną z metod klasyfikowania obszarów nieprawidłowych jest zastosowanie progu odpowiadającego 70% mediany wartości sygnału.

W tym podejściu próg obliczany jest oddzielnie dla każdej warstwy obrazu, a woksele znajdujące się poniżej 70% mediany zostają zakwalifikowane jako należące do defektu.

Autorzy wskazują jednak kilka ograniczeń. Po pierwsze, metoda jest matematycznie ograniczona: ponieważ mediana dzieli rozkład na dwie równe części, odsetek defektu wyznaczony w ten sposób nie może osiągnąć 50% ani więcej. Po drugie, próg zależy od rozkładu wartości w konkretnej warstwie obrazu. Przy względnie jednorodnym upośledzeniu funkcji płuca cały rozkład może przesunąć się w kierunku niższych wartości, a wraz z nim przesuwa się próg.

Po trzecie, klasyfikacja jest zasadniczo binarna: dany woksel należy albo do obszaru defektu, albo poza niego. Takie rozstrzygnięcie nie oddaje w pełni heterogeniczności regionalnej fizjologii płuca.

Badanie było retrospektywną analizą danych z wcześniejszych projektów

Badanie miało charakter eksploracyjnej, retrospektywnej analizy danych zgromadzonych podczas wcześniejszych badań prowadzonych w laboratorium autorów. Projekty uzyskały zgodę Ethics Committee of Bern, Switzerland, a uczestnicy lub ich opiekunowie prawni wyrazili pisemną świadomą zgodę.

Analizę podzielono na dwie główne części.

Najpierw metodę linear binning dostosowano do MP-MRI. Do stworzenia referencyjnych przedziałów wykorzystano 154 badania zdrowych dzieci, odpowiadające około 80% losowo wybranych danych z Bern-Basel Infant Lung Development cohort. Na podstawie tego zbioru określono cztery kategorie: defekt, niska funkcja, funkcja prawidłowa i wysoka.

Następnie zdolność metody do różnicowania dzieci zdrowych i dzieci z mukowiscydozą oceniono w zbiorze obejmującym pozostałe 41 zdrowych dzieci z tej kohorty oraz 30 dzieci z mukowiscydozą pochodzących z wcześniej opublikowanego badania.

W drugiej części wykorzystano wcześniej opublikowane dane 24 dzieci z mukowiscydozą leczonych elexacaftorem, tezacaftorem i ivacaftorem. Pierwsze MRI wykonano średnio 8,75 ± 4,69 miesiąca przed rozpoczęciem ETI, a kolejne średnio 4,70 ± 2,59 miesiąca po rozpoczęciu terapii.

Jak wykonywano MP-MRI?

Badania wykonywano przy użyciu klinicznego skanera całego ciała 1,5 T MAGNETOM Aera. Obrazowanie przeprowadzano w pozycji leżącej podczas swobodnego oddychania, wykorzystując wielowarstwową, dwuwymiarową sekwencję czasowo rozdzielczą ultra-fast balanced steady-state free precession. Do objęcia całej klatki piersiowej uzyskiwano około 8–14 przekrojów czołowych.

Obrazy przetwarzano algorytmem MP, uzyskując mapy miąższu płuc ważone wentylacją i perfuzją. Ostateczne maski płuc były oceniane wzrokowo przez doświadczonego badacza i w razie potrzeby ręcznie korygowane.

Sygnał struktur naczyniowych identyfikowano na mapach perfuzyjnych. Woksele o sygnale przekraczającym 90. percentyl klasyfikowano jako należące do naczyń, a następnie sygnał naczyniowy usuwano z map perfuzyjnych oraz wentylacyjnych przed dalszą analizą.

Linear binning wykorzystuje cztery kategorie funkcjonalne

W przeciwieństwie do progowania zależnego od mediany konkretnego obrazu linear binning opiera końcową klasyfikację na progach uzyskanych z dużej grupy zdrowych dzieci.

Nie oznacza to jednak, że cały proces jest całkowicie niezależny od badanego obrazu. Każda mapa wentylacyjna i perfuzyjna jest najpierw indywidualnie normalizowana. Autorzy przeprowadzili grid search różnych sposobów normalizacji, dążąc zarówno do uzyskania odpowiedniego kształtu histogramu referencyjnego, jak i dobrej zdolności różnicowania dzieci zdrowych i dzieci z mukowiscydozą.

Ostatecznie dla map perfuzyjnych wybrano punkt odcięcia odpowiadający medianie powiększonej o czterokrotność MAD, natomiast dla map wentylacyjnych – medianie powiększonej o pięciokrotność MAD. Wartości przekraczające odpowiedni punkt normalizacyjny ustawiano na 1.

Po normalizacji wartości wszystkich obrazów należących do zdrowego zbioru referencyjnego połączono we wspólne histogramy. Kategorie definiowano względem maksimum rozkładu i jego odchylenia standardowego:

• funkcja prawidłowa – wartości w zakresie maksimum ± 1 SD,
• funkcja wysoka – wartości powyżej maksimum + 1 SD,
• funkcja niska – wartości pomiędzy maksimum − 2 SD i maksimum − 1 SD,
• defekt – wartości poniżej maksimum − 2 SD.

Po wyznaczeniu referencyjnych granic każde kolejne badanie poddawano własnej normalizacji, a następnie poszczególne woksele przypisywano według stałych granic wyprowadzonych ze zdrowego zbioru referencyjnego.

Jak analizowano wyniki statystycznie?

Różnice pomiędzy zdrowymi dziećmi i dziećmi z mukowiscydozą oceniano jednostronnym testem Welcha, przy założeniu wyższych wartości parametrów nieprawidłowej funkcji w mukowiscydozie. Po korekcji Bonferroniego przyjęto poziom alfa 0,013. Przeprowadzona wcześniej analiza mocy wskazywała na minimalną liczebność 34 zdrowych dzieci i 24 dzieci z mukowiscydozą; rzeczywiste grupy liczyły odpowiednio 41 i 30 osób.

Zdolność różnicowania grup oceniano za pomocą krzywych ROC i pola pod krzywą – AUC.

W analizie danych przed i po rozpoczęciu ETI zastosowano jednostronny sparowany test t z poziomem alfa 0,01 po korekcji Bonferroniego. Ze względu na trudność wcześniejszego oszacowania wielkości efektu autorzy wykonali analizę czułości post hoc. Dla 24 osób, mocy 80% i alfa 0,05 minimalna wykrywalna wielkość efektu wyniosła 0,52.

Dzieci z mukowiscydozą miały więcej obszarów defektu i niskiej funkcji

Dzieci w zbiorze referencyjnym miały średnio 14,3 ± 5,3 roku, a 52,6% stanowiły dziewczęta. W grupie wykorzystywanej do oceny różnicowania zdrowe dzieci miały średnio 13,8 ± 5,5 roku, natomiast dzieci z mukowiscydozą 12,5 ± 3,2 roku. Dziewczęta stanowiły odpowiednio 46% i 66,7% tych grup.

W 24-osobowej kohorcie ETI średni wiek podczas badania poprzedzającego rozpoczęcie terapii wynosił 13,0 ± 2,27 roku, a podczas drugiego MRI 14,1 ± 2,30 roku. Dziewczęta stanowiły 41,2% grupy.

Histogramy zdrowych dzieci były zbliżone do rozkładów referencyjnych. U dzieci z mukowiscydozą obserwowano przesunięcie rozkładów sygnału w kierunku niższych wartości zarówno w mapach perfuzji, jak i wentylacji.

Perfuzja

Średni QDP_B, czyli odsetek defektu perfuzji określony metodą linear binning, wynosił 0,20% u zdrowych dzieci i 1,82% u dzieci z mukowiscydozą.

Odsetek obszarów o niskiej perfuzji LQP wynosił odpowiednio 9,37% i 18,82%. Dla obu porównań p<0,0001.

AUC wyniosło 0,88 dla QDP_B i 0,85 dla LQP. Dla klasycznego QDP_M, opartego na progu 70% mediany, AUC wynosiło 0,89.

Wentylacja

Średni VDP_B wynosił 1,18% u zdrowych dzieci i 7,05% u dzieci z mukowiscydozą. Średni LVP wynosił odpowiednio 11,23% oraz 21,17%. Również w obu tych porównaniach p<0,0001.

AUC wyniosło 0,90 dla VDP_B i 0,87 dla LVP. Dla klasycznego VDP_M wartość AUC wynosiła 0,89.

Linear binning nie wykazał więc wyraźnej przewagi nad klasycznym progiem 70% mediany w samym rozdzielaniu grup zdrowych dzieci i dzieci z mukowiscydozą. Potencjalna wartość metody polega przede wszystkim na większej granularności opisu regionalnych zaburzeń.

Podobny VDP, ale inna struktura zaburzeń

Autorzy przedstawili dwa celowo wybrane przykłady dzieci z kohorty przed rozpoczęciem ETI, u których klasyczna metoda dawała względnie podobne wartości VDP, natomiast linear binning ujawniał różną strukturę nieprawidłowości.

U 12-letniej dziewczynki z LCI 13,0 VDP obliczony metodą 70% mediany wynosił 37,6%. Linear binning wykazał VDP_B równy 30,4% oraz LVP 28,6%. Łącznie kategorie defektu i niskiej wentylacji obejmowały 59% wokseli miąższu płucnego.

U 14-letniego chłopca z LCI 10,3 klasyczny VDP wynosił 32,6%. VDP_B wynosił natomiast 19,9%, a LVP 31,3%. Łącznie kategorie te obejmowały 48,2% analizowanych wokseli.

Przypadki te ilustrują podstawową ideę metody: podobny całkowity VDP wyznaczony na podstawie mediany nie musi oznaczać podobnego rozkładu nasilenia regionalnych zaburzeń.

Zmiany obserwowane po rozpoczęciu ETI

W analizie 24 dzieci z mukowiscydozą porównano parametry MP-MRI przed rozpoczęciem ETI i podczas leczenia. Po rozpoczęciu terapii średnie histogramy zarówno wentylacji, jak i perfuzji przesuwały się w kierunku wartości charakterystycznych dla zdrowej referencji.

Zmiany wentylacji

W zakresie wentylacji obie metody wskazywały na poprawę analizowanych parametrów.

Średni VDP_B zmniejszył się z 9,17% przed rozpoczęciem ETI do 3,49% w drugim badaniu. Średnia różnica wynosiła 5,68 punktu procentowego, p<0,0001.

LVP zmniejszył się z 21,75% do 15,90%. Średnia różnica wyniosła 5,85 punktu procentowego, p=0,0019.

Klasyczny VDP_M zmniejszył się z 24,21% do 18,88%, co odpowiadało średniej różnicy 5,34 punktu procentowego, p<0,0001.

Perfuzja ujawniła różnicę między metodami

W przypadku perfuzji klasyczne progowanie wskazywało na zmniejszenie QDP_M z 21,64% do 19,00%. Średnia różnica wyniosła 2,64 punktu procentowego, p=0,0029.

Linear binning pozwolił jednak rozdzielić tę zmianę na dwa komponenty.

LQP, czyli udział obszarów o niskiej perfuzji, zmniejszył się z 17,14% do 13,03%. Średnia różnica wyniosła 4,11 punktu procentowego, p=0,002.

QDP_B, czyli udział obszarów należących do kategorii najgłębszego defektu perfuzji, zmniejszył się natomiast jedynie z 1,95% do 1,53%. Średnia różnica wynosiła 0,43 punktu procentowego i nie osiągnęła istotności statystycznej – p=0,21.

Jest to jeden z najważniejszych rezultatów publikacji. Nie oznacza on, że klasyczny sposób analizy jest błędny. Obie metody stosują inne reguły klasyfikacji wokseli i dlatego opisują nieprawidłowość w odmienny sposób. Linear binning pozwala jednak rozdzielić część zmian ukrytą wcześniej we wspólnej kategorii „defektu”.

Przykład pojedynczej pacjentki

Autorzy pokazali także przypadek dziewczynki ocenianej przed i po rozpoczęciu terapii. Podczas pierwszego badania miała 11 lat i LCI 7,1, natomiast podczas kolejnego 12 lat i LCI 6,7.

Na mapach perfuzji klasyczny QDP_M był bardzo podobny przed i po rozpoczęciu leczenia: 24,4% i 24,8%. Linear binning ujawnił jednak zmniejszenie LQP z 23,53% do 15,82%, przy jednoczesnym wzroście QDP_B z 1,98% do 4,46%.

W mapach wentylacji klasyczny VDP_M zmniejszył się z 31,1% do 17,2%. VDP_B zmniejszył się z 8,57% do 0,74%, a LVP z 36,77% do 9,35%.

Jest to jednak opis pojedynczego przypadku i nie należy na jego podstawie wnioskować o typowej odpowiedzi wszystkich dzieci leczonych ETI.

Co rzeczywiście można wywnioskować z badania?

Największą potencjalną zaletą linear binning jest możliwość bardziej szczegółowego opisania regionalnej heterogeniczności płuca. Zamiast jednej kategorii „defektu” uzyskuje się przynajmniej dwie klasy odpowiadające różnemu nasileniu obniżenia sygnału funkcjonalnego.

Może mieć to znaczenie szczególnie w badaniach longitudinalnych. Obszar płuca może bowiem zmieniać kategorię z defektu do niskiej funkcji albo z niskiej funkcji do zakresu prawidłowego, a klasyczny binarny parametr nie zawsze pokaże charakter takiej zmiany.

Jednocześnie publikacja nie wykazuje jeszcze klinicznej przewagi linear binning nad dotychczas stosowanymi metodami. Nie ustalono klinicznie istotnych wartości granicznych, nie przeprowadzono zewnętrznej walidacji oraz nie wykazano, że zastosowanie nowej analizy poprawia decyzje kliniczne lub wyniki leczenia.

Potrzebna jest zewnętrzna walidacja

Autorzy podkreślają, że brak standaryzacji protokołów akwizycji i analizy funkcjonalnego MRI pomiędzy ośrodkami może wpływać na uzyskiwane wartości. Nie wiadomo więc, czy opracowane progi referencyjne można zastosować bez zmian do obrazów pochodzących z innych skanerów, protokołów i laboratoriów.

Drugą nierozwiązaną kwestią jest wiek. Wartości referencyjne opracowano na podstawie populacji pediatrycznej i autorzy zastrzegają, że mogą one mieć zastosowanie wyłącznie w tym zakresie wieku. Nie wiadomo również, czy nawet w obrębie pediatrii potrzebne będą odrębne normy zależne od wieku.

Kolejnym etapem powinno być porównanie zaadaptowanego linear binning z wynikami pochodzącymi z innych ośrodków wykonujących MP-MRI, innymi wariantami linear binning stosowanymi w funkcjonalnym MRI płuc, metodami takimi jak adaptive k-means oraz klinicznymi badaniami czynności płuc.

Badanie zostało wsparte przez Cystische Fibrose Schweiz oraz Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung, granty 10002556 i 10004139. Philipp Latzin ujawnił granty lub umowy z Vertex i OM Pharma, honoraria m.in. od Vertex, Vifor i OM Pharma oraz udział w gremiach doradczych kilku firm; autorzy zaznaczyli, że relacje te nie były związane z prezentowanym badaniem. Elisabeth Kieninger zgłosiła honoraria od Sanofi-Aventis i Vertex poza zakresem przedstawionej pracy.

Źródło: Magnetic Resonance in Medicine, Adapted Linear Binning Method to Assess Pulmonary Ventilation and Perfusion in Children Using Lung Matrix Pencil Decomposition (MP)-MRI.
DOI: https://doi.org/10.1002/mrm.70544