Het eenvoudigste 3D PREFUL-beeldvormingsprotocol bestaat uit twee sequenties. Eén anatomische localizer + één meting tijdens vrije ademhaling met behulp van een golden-angle 3D stack-of-stars acquisitie. Het stack-of-stars traject is weergegeven in Figuur 1Deze traject combineert radiale bemonstering (in het vlak) met lineaire cartesiaanse bemonstering (in de z-richting). Aanvankelijk worden monsters in de k-z de partitierichtingen worden verkregen, gevolgd door de rotatie met de gouden hoek om de k-ruimte te samplenHet spijt me, maar u heeft geen brontekst opgegeven om te vertalen. Voer a.u.b. de Engelse tekst in die u naar het Nederlands vertaald wilt hebben.-ky vlak. Het acquisitieschema met de gouden hoek bevordert self-gating met een uniforme dekking van de k-ruimte, wat beeldreconstructie met sterk ondergesamplede gegevens vergemakkelijkt. Het gehele onderzoek duurt niet langer dan 10 minuten. Verdere details van de 3D PREFUL-acquisitie zijn opgenomen in de Tabel met materialen.

Figuur 1: Trajectorie van een geroteerde stapel sterren volgens de gulden snede. Let op de radiale bemonstering toegepast langs de in-plane dimensie (kx-ky) en de Cartesiaanse bemonstering toegepast langs de slice-richting (kz). Dit schema resulteert in een cilindrische dekking van de verworven k-ruimte. In dit voorbeeld worden drie partities afgebeeld, elk met 32 radiale projecties per partitie. Klik hier om een grotere versie van deze figuur te bekijken.
Een lijst met typische MRI-parameters die worden gebruikt voor 3D PREFUL-beeldvorming op een 1.5 T scanner wordt gepresenteerd in Tabel 1B. Voor de 3T scanner worden de sequentieparameters vermeld in Tabel 2 voorgesteld.
Nadat de gegevensverwerving en gegevensoverdracht zijn voltooid, begint het deel van de beeldreconstructie. Deze stap wordt uitgevoerd door een MATLAB-code die deel uitmaakt van het intellectueel eigendom (IP) van een bedrijf (BioVisioneers GmbH) en niet openbaar gedeeld kan worden. Het script is volledig geautomatiseerd zodra het pad naar de ruwe gegevens is opgegeven. Een schematisch overzicht van de reconstructieprocedure is afgebeeld in Figuur 2.

Figuur 2: Schematisch overzicht van de 3D Phase Resolved Functional Lung (PREFUL) MRI-methode. Eerst worden de gegevens verzameld met een 8 minuten durende MR-acquisitie bij vrije ademhaling met een stack-of-stars traject. Nadat de gegevens van de MR-scanner naar een rekeneenheid zijn overgebracht, worden 3D-beelden met een lage resolutie en een temporele resolutie van ongeveer 100 ms gereconstrueerd. Het longparenchym van elk beeld wordt gesegmenteerd en het longvolume wordt berekend. De informatie over het longvolume wordt vervolgens gebruikt als gatingsignaal. Op basis van de amplitude en fase van het gatingsignaal worden de radiale projecties gesorteerd in respiratoire bins die één ademhalingscyclus beslaan. Het binning-gedeelte wordt gevolgd door de dynamische reconstructie van beelden met volledige resolutie, die vervolgens worden geregistreerd naar het niveau van de maximale inspiratie. Na verschillende postprocessing-stappen wordt de regionale ventilatiecyclus (RVent) berekend en wordt de beeldanalyse, inclusief parameterextractie, uitgevoerd. De RVent-cyclus wordt beoordeeld door de flow-volumelus (FVL) voor elke voxel te berekenen. De FVL-CM ventilatiekaarten worden geëxtraheerd door de FVL van elke voxel te beoordelen ten opzichte van een gezonde referentie-FVL, waarbij de gelijkenis wordt geëvalueerd met behulp van een cross-correlatiemetriek. Voor zowel RVent als FVL-CM worden de globale waarden van het totale ventilatiedefectpercentage (VDP) gekwantificeerd. Bovendien wordt de dynamiek van de RVent-cyclus geanalyseerd via time-to-peak analyse, wat resulteert in een ventilatie time-to-peak (VTTP) parametekaart. Daarnaast wordt de afwijking van de verwachte piekventilatie bij 50% van de RVent-cyclus gekwantificeerd in de VTTPDev kaart. Klik hier om een grotere versie van deze figuur te bekijken.
Na de procedure voor beeldreconstructie worden alle beelden ruimtelijk uitgelijnd naar één vaste toestand. De registratie zelf wordt op de CPU uitgevoerd met Advanced Normalization Tools (ANTs27) of het registratiepakket van Forsberg28. Hoewel het ANTs-pakket de gouden standaard blijft voor beeldregistratietaken bij MRI- en CT-beeldvorming, maakt de registratie van Forsberg een 10 keer snellere registratieprocedure mogelijk van maximaal 9 minuten met vergelijkbare resultaten19. Het registratiepakket kan worden gekozen afhankelijk van de prioriteiten van de gebruiker en de beschikbare computereenheden. Nadat de registratie is voltooid, worden de geregistreerde morfologische beelden opnieuw geïnverteerd, zodat de long donker verschijnt in het grijswaardenbeeld.
In Figuur 3 wordt de volgende stap van de pijplijn weergegeven, die bestaat uit de segmentatie van het longparenchym. Eerst wordt het longparenchym gesegmenteerd uit het beeld aan het einde van de inspiratie met behulp van een convolutioneel neuraal netwerk met nnUnet-architectuur29. De segmentatie van het longparenchym wordt gevolgd door vaatherkenning30, waarbij de vaten worden uitgesloten van het uiteindelijke segmentatiemasker.

Figuur 3: Voorbeeldresultaten van deep learning-gebaseerde segmentatie voor een 32-jarige mannelijke proefpersoon. De bovenste rij toont morfologische beelden van acht representatieve coronale coupes. In de tweede en derde rij kunnen respectievelijk het bijbehorende masker van het longparenchym en het uiteindelijke masker met uitsluiting van vaten worden waargenomen. Klik hier om een grotere versie van deze figuur te bekijken.
Een voorbeeld van een onjuiste segmentatie van het longparenchym is afgebeeld in Figuur 4. Het is essentieel om de op deep learning gebaseerde segmentaties visueel te inspecteren; indien deze onbevredigend worden geacht, moeten handmatige correcties worden overwogen om de nauwkeurigheid van het uiteindelijke longparenchymmasker te verbeteren.

Figuur 4: Voorbeeld van onjuiste deep learning-gebaseerde segmentatie bij een mannelijke proefpersoon van 57 jaar. De bovenste rij toont morfologische beelden van acht representatieve coronale coupes. In de tweede en derde rij zijn respectievelijk het bijbehorende masker van het longparenchym en het definitieve masker met uitsluiting van vaten te zien. Zoals duidelijk is, worden verschillende fibrotische regio's onterecht herkend als vaten of niet-pulmonale structuren. Deze onjuistheden zijn handmatig gecorrigeerd, zoals gedemonstreerd in de vierde rij. Klik hier om een grotere versie van deze figuur te bekijken.
Na verschillende filtratiestappen die in het protocol worden beschreven, worden surrogaten voor ventilatie berekend. 3D PREFUL MRI genereert kwantitatieve kaarten voor statische regionale ventilatie (RVent), dynamische flow-volumecorrelatiemetrieken (FVL-CM) en twee parameters gebaseerd op ventilatie tijd-tot-piek (VTTP) analyse, zoals geïllustreerd in Figuur 5. Deze figuur toont acht coronale coupes van een gezonde man van 32 jaar oud. Let op de verwachte homogene distributie van alle ventilatieparameters.

Figuur 5: 3D PREFUL MRI van een gezonde vrijwilliger (man van 32 jaar oud). Representatieve morfologische (bovenste rij) en 3D PREFUL MRI ventilatieparameterkaarten (tweede tot vijfde rij) voor een gezonde vrijwilliger (man van 32 jaar oud). De statische regionale ventilatie wordt weergegeven door de regionale ventilatie (RVent), terwijl de ventilatiedynamiek wordt beoordeeld met de flow-volumeloop correlatiemetriek (FVL-CM); de ventilatietijd-tot-piek (VTTP) en de afwijking van VTTP (VTTPDev) zijn ventilatieparameters die de ventilatiedynamiek beoordelen. Zoals verwacht worden homogene ventilatiewaarden waargenomen voor alle ventilatieparameterkaarten. Klik hier om een grotere versie van deze figuur te bekijken.
Om de ventilatiekaarten te vereenvoudigen, worden ventilatiedefect-kaarten (VD-kaarten) afgeleid voor RVent en FVL-CM, wat een snellere interpretatie van de resultaten mogelijk maakt. De voorbeeld-VD-kaarten worden gepresenteerd in Figuur 6. Voor zowel VDRVent als VDFVL-CM bleken de VDP-waarden respectievelijk 3,6% en 3,0% te zijn, wat binnen het gezonde normale bereik valt. Idealiter, en afhankelijk van de leeftijd van de gezonde vrijwilligers, mag de VDP-waarde 10% niet overschrijden. De bovengenoemde parameters (RVent, FVL-CM en hun VD-kaarten) zijn in verschillende studies gevalideerd11,12,35,36 en zijn gevoelig voor de detectie van ziekten alsovel als voor de detectie van therapie-geïnduceerde effecten20,34,37,38,39.

Figuur 6: Representatieve RVent- en FVL-CM-kaarten, inclusief de bijbehorende ventilatiedefectkaarten. De RVent- (bovenste rij) en FVL-CM-kaarten (tweede rij), inclusief de ventilatiedefectkaarten (derde en vierde rij) afgeleid via 3D PREFUL MRI voor een gezonde vrijwilliger (man van 32 jaar), worden hier weergegeven. In de VD-kaarten zijn gezonde regio's groen en zijn gebieden met een ventilatiedefect in rood gemarkeerd. Klik hier om een grotere versie van deze figuur te bekijken.
Verschillende pulmonale MRI-studies zijn gerapporteerd bij magnetische veldsterktes van zowel 1,5 T als 3 T. Hoewel 3 T een theoretisch voordeel biedt vanwege een hogere signaal-ruisverhouding (SNR), kan dit voordeel tenietgedaan worden door meer uitgesproken magnetische susceptibiliteitseffecten bij 3 T. De exacte invloed van magnetische veldsterktes op de beeldkwaliteit van 3D PREFUL-ventilatiekaarten is momenteel onbekend. Hier presenteren we haalbaarheidsresultaten (Figuur 7) verkregen voor een gezonde vrijwilliger (man van 35 jaar) met gebruik van een 3 T MR-scanner. Let op het meer heterogene uiterlijk van de ventilatieparameters bij 3 T in vergelijking met 1,5 T (Figuur 5).

Figuur 7: 3D PREFUL MRI van een gezonde vrijwilliger (man van 35 jaar). Representatieve morfologische (bovenste rij) en 3D PREFUL MRI ventilatieparameterkaarten (tweede tot vijfde rij) voor een gezonde vrijwilliger (man van 35 jaar) bij 3T. De statische regionale ventilatie wordt weergegeven door de regionale ventilatie (RVent), terwijl de ventilatiedynamiek wordt beoordeeld met behulp van de flow-volume loop correlatiemetriek (FVL-CM), ventilatietijd-tot-piek (VTTP), en afwijking van VTTP (VTTPDev), welke ventilatieparameters zijn voor de beoordeling van de ventilatiedynamiek. Klik hier om een grotere versie van deze figuur te bekijken.
Een vergelijking/validatie van 3D PREFUL, afgeleid van directere metingen, is momenteel niet gepubliceerd. Er zijn verschillende studies die een positieve correlatie rapporteren tussen VDP-waarden en de ruimtelijke overlap van ventilatiedefecten tussen de 2D PREFUL-techniek en 129Xe11,12,14 en 19F MRI35. Een recente vergelijking tussen ventilatie afgeleid van 3D PREFUL en directe ventilatiemetingen met 19F MRI bij patiënten met chronische obstructieve longziekte (COPD), astma en gezonde vrijwilligers toonde een matige tot sterke correlatie op globaal niveau40. Er wordt een voorbeeldvergelijking gepresenteerd (Figuur 8) voor een patiënte met de diagnose chronische obstructieve longziekte (54-jarige vrouw, FEV1 = 42% van de voorspelde waarde), die werd onderzocht met behulp van 3D PREFUL en 19F wash-in MRI.

Figuur 8: 3D PREFUL MRI van een 54-jarige vrouwelijke COPD-patiënt. Getoond worden morfologische beelden (bovenste rij), ventilatieparameters en de bijbehorende ventilatiedefectkaarten van een 54-jarige vrouwelijke COPD-patiënt (FEV1 = 42 % pred., FVC = 102% pred.) verkregen via 3D PREFUL MRI (2e en 3e rij) en 19F MRI (4e en 5e rij), samen met een vergelijking van de ventilatiedefectkaarten van beide methoden (onderste rij). In totaal, over alle coupes heen, was de Sørensen-Dice-coëfficiënt 25,5% in defecte gebieden en 80,6% in gezonde regio's, wat resulteerde in een totale ruimtelijke overlap van 69,2%. Opvallend is de visuele correlatie tussen overeenkomstige ventilatiegebieden in gezonde en defecte regio's (wegegeven in donkergroen). Klik hier om een grotere versie van deze figuur te bekijken.
3D PREFUL MRI kan worden gebruikt om regionale reacties op therapie te meten20. Als voorbeeld toont Figuur 9 een vergelijking van drie coupes van de long van een subject met cystische fibrose (43-jarige vrouw) bij baseline (links, FEV1 = 94% van de voorspelde waarde) en na CFTR-modulatortherapie (rechts, FEV1 = 112% van de voorspelde waarde). Omdat beide metingen ruimtelijk overeenkomen, maakt 3D PREFUL regionale analyse van de behandelingsrespons mogelijk, zoals geïllustreerd onderaan Figuur 9. Let op de verhoogde waarden van de flow-volume-loop correlatiemetriek alsook de afname van ventilatiedefecten na therapie.

Figuur 9: 3D PREFUL-metingen van een 43-jarige vrouwelijke patiënt met cystische fibrose (CF). Voorbeeld van ventilatiemarkerkaarten van baseline (links) en na behandeling (rechts) 3D PREFUL-metingen van een 43-jarige vrouwelijke CF-patiënt. VDPFVL-CM nam af van 18,0% (baseline) naar 3,6% (na behandeling). FEV1 % pred. baseline: 94%, FEV1 % pred. na behandeling: 112%. LCI baseline: 10,48, LCI na behandeling: 9,39. De bijbehorende responskaarten op de behandeling zijn onderaan afgebeeld. Let op de groene gebieden die de herstelde ventilatie na de therapie aangeven. Klik hier om een grotere versie van deze figuur te bekijken.
Tabel 1: Lijst van typische MRI-parameters gebruikt voor de localizer en 3D PREFUL-acquisitie op een 1.5T scanner. (A) Localizer en (B) 3D PREFUL-acquisitie. Klik hier om deze tabel te downloaden.
Tabel 2: Voorgestelde MRI-parameters voor 3D PREFUL-acquisitie op een 3T-scanner. Klik hier om deze tabel te downloaden.
Tabel 3: Voorbeeld van een statistisch rapport van de 3D PREFUL-parameters voor een gezonde vrijwilliger (man van 35 jaar) Klik hier om deze tabel te downloaden.