Il protocollo di imaging 3D PREFUL più semplice comprende due sequenze. Un localizzatore anatomico + una misurazione in libera respirazione utilizzando l'acquisizione 3D a stack-of-stars con angolo aureo. La traiettoria stack-of-stars è mostrata in Figura 1. Questa traiettoria combina un campionamento radiale (nel piano) con un campionamento cartesiano lineare (nella direzione z). Inizialmente vengono acquisiti i campioni nella direzione della partizione kz, seguiti dalla rotazione dell'angolo aureo per campionare il piano kx-ky. Lo schema di acquisizione con angolo aureo favorisce l'auto-gating con una copertura uniforme dello spazio-k, facilitando la ricostruzione delle immagini anche con dati fortemente sottocampionati. L'esame completo non supera i 10 minuti. Ulteriori dettagli sull'acquisizione 3D PREFUL sono riportati nella Tabella dei Materiali.

Figura 1: Traiettoria a pila di stelle ruotata secondo il rapporto aureo. Si noti il campionamento radiale applicato lungo la dimensione nel piano (kx-ky) e il campionamento cartesiano applicato lungo la direzione della fetta (kz). Questo schema determina una copertura cilindrica dello spazio-k acquisito. In questo esempio, sono rappresentate tre partizioni, ciascuna comprendente 32 proiezioni radiali per partizione. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Una lista dei parametri tipici di risonanza magnetica utilizzati per l'immagine 3D PREFUL su scanner a 1,5 T è presentata nella Tabella 1B. Per lo scanner a 3 T, si propongono i parametri di sequenza elencati nella Tabella 2.
Al termine dell'acquisizione e del trasferimento dei dati, inizia la fase di ricostruzione delle immagini. Questo passaggio viene eseguito mediante un codice MATLAB di proprietà intellettuale (IP) di un'azienda (BioVisioneers GmbH) e non può essere condiviso pubblicamente. Lo script è completamente automatizzato una volta fornito il percorso verso i dati grezzi. Una panoramica schematica della procedura di ricostruzione è mostrata in Figura 2.

Figura 2: Panoramica schematica del metodo di risonanza magnetica polmonare funzionale 3D in fase risolta (PREFUL). Inizialmente, i dati vengono acquisiti mediante un'acquisizione RM a libero respiro della durata di 8 minuti, utilizzando una traiettoria a stack-of-stars. Dopo il trasferimento dei dati dallo scanner RM a un'unità di calcolo, vengono ricostruite immagini 3D a bassa risoluzione con una risoluzione temporale di circa 100 ms. Viene segmentato il parenchima polmonare di ciascuna immagine e calcolato il volume polmonare. L'informazione sul volume polmonare viene successivamente utilizzata come segnale di gating. Sulla base dell'ampiezza e della fase del segnale di gating, le proiezioni radiali vengono suddivise in bin respiratori che coprono un intero ciclo respiratorio. A questa suddivisione segue la ricostruzione dinamica di immagini a risoluzione completa, che vengono quindi registrate al livello di fine inspirazione. Dopo diversi passaggi di post-elaborazione, viene calcolato il ciclo di ventilazione regionale (RVent) e viene eseguita l'analisi delle immagini, inclusa l'estrazione dei parametri. Il ciclo RVent viene valutato calcolando il loop flusso-volume (FVL) per ciascun voxel. Le mappe di ventilazione FVL-CM vengono estratte analizzando ciascun FVL del voxel rispetto a un FVL di riferimento sano, con la similarità valutata mediante metrica di cross-correlazione. Per entrambi RVent e FVL-CM, vengono quantificati i valori globali percentuali totali di difetto di ventilazione (VDP). Inoltre, la dinamica del ciclo RVent viene analizzata mediante l'analisi del tempo per raggiungere il picco, ottenendo una mappa parametrica del tempo di picco della ventilazione (VTTP). In aggiunta, viene quantificata la deviazione rispetto al picco di ventilazione atteso al 50% del ciclo RVent, attraverso la mappa VTTPDev. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Dopo la procedura di ricostruzione dell'immagine, tutte le immagini vengono allineate spazialmente a uno stato fisso singolo. Il processo di registrazione viene eseguito sulla CPU utilizzando Advanced Normalization Tools (ANTs27) o il pacchetto di registrazione Forsberg28. Mentre il pacchetto ANTs rimane lo standard di riferimento per i compiti di registrazione delle immagini in risonanza magnetica e tomografia computerizzata, il pacchetto Forsberg consente una procedura di registrazione fino a 10 volte più veloce, con tempi ridotti fino a 9 minuti, ottenendo risultati comparabili19. La scelta del pacchetto di registrazione può dipendere dalle priorità dell'utente e dalle unità di calcolo disponibili. Al termine della registrazione, le immagini morfologiche registrate vengono nuovamente reinvertite in modo che il polmone appaia scuro nell'immagine in scala di grigi.
In Figura 3 è rappresentato il passaggio successivo del flusso di lavoro, che prevede la segmentazione del parenchima polmonare. In primo luogo, il parenchima polmonare viene segmentato dall'immagine di fine inspirazione utilizzando una rete neurale convoluzionale con architettura nnUnet29. Alla segmentazione del parenchima polmonare segue un riconoscimento dei vasi30, che viene escluso dalla maschera di segmentazione finale.

Figura 3: Risultati esemplificativi della segmentazione basata sull'apprendimento profondo per un soggetto maschio di 32 anni. La riga superiore mostra immagini morfologiche di otto fette coronali rappresentative. Nella seconda e nella terza riga si possono osservare rispettivamente la maschera del parenchima polmonare corrispondente e la maschera finale con l'esclusione dei vasi. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Un esempio di segmentazione errata del parenchima polmonare è mostrato in Figura 4. È essenziale esaminare visivamente le segmentazioni basate sull'apprendimento profondo e, se ritenute insoddisfacenti, si dovrebbero considerare correzioni manuali per migliorare l'accuratezza della maschera finale del parenchima polmonare.

Figura 4: Risultati esemplificativi di una segmentazione basata su apprendimento profondo errata per un soggetto maschio di 57 anni. La riga superiore mostra immagini morfologiche di otto sezioni coronali rappresentative. Nella seconda e nella terza riga, si possono osservare rispettivamente la maschera del parenchima polmonare corrispondente e la maschera finale con l'esclusione dei vasi. Come evidente, diverse aree fibrotiche vengono erroneamente riconosciute come vasi o strutture non polmonari. Queste imprecisioni sono state corrette manualmente, come mostrato nella quarta riga. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Dopo diversi passaggi di filtraggio descritti nel protocollo, vengono calcolati i parametri surrogati della ventilazione. La risonanza magnetica 3D PREFUL genera mappe quantitative per la ventilazione regionale statica (RVent), le metriche dinamiche di correlazione flusso-volume (FVL-CM) e due parametri basati sull'analisi del tempo di picco della ventilazione (VTTP), come illustrato in Figura 5. Questa figura mostra otto sezioni coronali di un maschio sano di 32 anni. Si noti la distribuzione omogenea attesa di tutti i parametri di ventilazione.

Figura 5: Risonanza magnetica 3D PREFUL di un volontario sano (maschio di 32 anni). Mappe rappresentative dei parametri morfologici (prima riga) e della ventilazione mediante risonanza magnetica 3D PREFUL (dalla seconda alla quinta riga) per un volontario sano (maschio di 32 anni). La ventilazione regionale statica è rappresentata dalla ventilazione regionale (RVent), mentre la dinamica della ventilazione viene valutata mediante l'indice di correlazione del loop flusso-volume (FVL-CM), il tempo per raggiungere il picco di ventilazione (VTTP) e la deviazione del VTTP (VTTPDev), parametri che valutano la dinamica della ventilazione. Come previsto, si osservano valori di ventilazione omogenei in tutte le mappe dei parametri di ventilazione. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Per semplificare le mappe di ventilazione, vengono derivate mappe di difetto ventilatorio (VD) per RVent e FVL-CM, che consentono un'interpretazione più rapida dei risultati. Le mappe VD esemplificative sono presentate in Figura 6. Per entrambe VDRVent e VDFVL-CM, i valori di VDP sono risultati rispettivamente pari a 3,6% e 3,0%, rientrando nell'intervallo normale sano. Idealmente, e in base all'età dei volontari sani, il valore di VDP non dovrebbe superare il 10%. I parametri sopra menzionati (RVent, FVL-CM e le relative mappe VD) sono stati validati in diversi studi11,12,35,36 e risultano sensibili nel rilevamento delle malattie nonché degli effetti indotti dalla terapia20,34,37,38,39.

Figura 6: Mappe rappresentative di RVent e FVL-CM, incluse le relative mappe di difetto ventilatorio. Sono mostrate le mappe di RVent (prima riga) e FVL-CM (seconda riga), incluse le mappe di difetto ventilatorio (terza e quarta riga) ottenute mediante risonanza magnetica 3D PREFUL per un volontario sano (maschio di 32 anni). Nelle mappe di difetto ventilatorio (VD), le regioni sane sono in verde, mentre le aree con deficit ventilatorio sono evidenziate in rosso. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Diversi studi di risonanza magnetica polmonare sono stati riportati a intensità di campo magnetico di 1,5 T e 3 T. Sebbene esista un vantaggio teorico del 3 T dovuto al maggiore rapporto segnale-rumore (SNR), tale vantaggio potrebbe essere superato da effetti di suscettibilità magnetica più marcati a 3 T. L'influenza esatta delle diverse intensità di campo magnetico sulla qualità dell'immagine delle mappe di ventilazione 3D PREFUL è attualmente sconosciuta. Presentiamo qui i risultati di fattibilità (Figura 7) ottenuti in un volontario sano (maschio di 35 anni) utilizzando uno scanner RM a 3 T. Si noti l'aspetto più eterogeneo dei parametri di ventilazione a 3 T rispetto a 1,5 T (Figura 5).

Figura 7: Risonanza magnetica 3D PREFUL di un volontario sano (maschio di 35 anni). Mappe rappresentative di parametri morfologici (prima riga) e di ventilazione ottenuti con risonanza magnetica 3D PREFUL (dalla seconda alla quinta riga) per un volontario sano (maschio di 35 anni) a 3T. La ventilazione regionale statica è rappresentata dalla ventilazione regionale (RVent), mentre la dinamica ventilatoria viene valutata mediante l'indice di correlazione del loop flusso-volume (FVL-CM), il tempo di raggiungimento del picco di ventilazione (VTTP) e la deviazione del VTTP (VTTPDev), parametri che valutano la dinamica della ventilazione. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Al momento non sono stati pubblicati studi di confronto o validazione delle misurazioni tridimensionali di PREFUL con metodi più diretti. Diversi studi riportano una correlazione positiva tra i valori di VDP e la sovrapposizione spaziale dei difetti di ventilazione ottenuti con la tecnica bidimensionale di PREFUL e con la risonanza magnetica con 129Xe11,12,14 e con la risonanza magnetica con 19F35. Un recente confronto tra la ventilazione derivata da PREFUL tridimensionale e misurazioni dirette della ventilazione mediante risonanza magnetica con 19F, effettuato su pazienti affetti da malattia polmonare ostruttiva cronica (BPCO), asma e su volontari sani, ha evidenziato una correlazione da moderata a forte a livello globale40. Viene presentato un esempio di confronto (Figura 8) relativo a un paziente diagnosticato con malattia polmonare ostruttiva cronica (donna di 54 anni, FEV1 = 42% del valore predetto), sottoposto a esame mediante PREFUL tridimensionale e risonanza magnetica con wash-in di 19F.

Figura 8: Risonanza magnetica 3D PREFUL di una paziente donna di 54 anni affetta da BPCO. Vengono presentate immagini morfologiche (prima riga), parametri di ventilazione e rispettive mappe dei difetti di ventilazione di una paziente donna di 54 anni affetta da BPCO (FEV1 = 42 % pred., FVC = 102% pred.) ottenute mediante risonanza magnetica 3D PREFUL (seconda e terza riga) e risonanza magnetica con 19F (quarta e quinta riga), insieme a un confronto tra le mappe dei difetti di ventilazione ottenute con i due metodi (ultima riga). Nel complesso, considerando tutte le sezioni, il coefficiente di Sørensen-Dice è stato del 25,5% nelle aree difettose e dell'80,6% nelle regioni sane, con un sovrapposizione spaziale totale del 69,2%. È degno di nota il riscontro visivo di una correlazione tra le aree di ventilazione sovrapposte nelle regioni sane e quelle difettose (rappresentate in verde scuro). Cliccare qui per visualizzare una versione ingrandita di questa figura.
La risonanza magnetica 3D PREFUL può essere utilizzata per misurare le risposte regionali alla terapia20. Ad esempio, Figura 9 mostra un confronto tra tre sezioni del polmone di un soggetto affetto da fibrosi cistica (donna di 43 anni) allo stato basale (a sinistra, FEV1 = 94% del valore predetto) e dopo la terapia con modulatori del CFTR (a destra, FEV1 = 112% del valore predetto). Poiché entrambe le misurazioni sono spazialmente allineate, la tecnica 3D PREFUL consente un'analisi regionale della risposta al trattamento, come illustrato nella parte inferiore di Figura 9. Si noti l'aumento dei valori della metrica di correlazione del loop flusso-volume, nonché la riduzione dei difetti di ventilazione dopo la terapia.

Figura 9: Misure 3D PREFUL di una paziente femmina di 43 anni affetta da fibrosi cistica (FC). Mappe rappresentative dei marcatori di ventilazione pre-trattamento (a sinistra) e post-trattamento (a destra) di una paziente femmina di 43 anni affetta da FC. Il VDPFVL-CM è diminuito dal 18,0% (pre-trattamento) al 3,6% (post-trattamento). FEV1 % pred. pre-trattamento: 94%, FEV1 % pred. post-trattamento: 112%. LCI pre-trattamento: 10,48, LCI post-trattamento: 9,39. Le corrispondenti mappe di risposta al trattamento sono mostrate in basso. Si notino le aree verdi che indicano la risoluzione della ventilazione dopo la terapia. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Tabella 1: Elenco dei parametri tipici di risonanza magnetica utilizzati per l'acquisizione del localizzatore e della sequenza 3D PREFUL su uno scanner a 1,5 T. (A) Localizzatore e (B) acquisizione 3D PREFUL. Cliccare qui per scaricare questa tabella.
Tabella 2: Parametri proposti per l'acquisizione MRI 3D PREFUL su uno scanner 3T. Cliccare qui per scaricare questa tabella.
Tabella 3: Esempio di rapporto statistico dei parametri 3D PREFUL per un volontario sano (maschio di 35 anni) Cliccare qui per scaricare questa tabella.