L’imagerie par résonance magnétique (IRM) en gaz hyperpolarisé est apparue pour la première fois comme une nouvelle modalité d’imagerie pulmonaire fonctionnelle pour évaluer la distribution de la ventilation pulmonaire il y a près de trois décennies1. Depuis lors, des études de recherche utilisant l’IRM en gaz hyperpolarisé ont révélé de nombreuses informations sur la nature de la fonction pulmonaire chez les patients atteints de maladies pulmonaires chroniques telles que l’asthme, la bronchopneumopathie chronique obstructive (BPCO) et la fibrose kystique 2,3,4,5,6. Le gaz 3He hyperpolarisé et le gaz 129Xe ont été utilisés historiquement ; cependant, le 129Xe est maintenant le principal agent inhalé en raison de la disponibilité limitée du gaz 3He. 129Xe se diffuse également librement à travers la membrane alvéolaire et est absorbé par les globules rouges dans les capillaires pulmonaires ; dans cette « phase dissoute », le 129Xe résonne à des fréquences uniques qui permettent de mesurer les échanges gazeux régionaux en un seul balayage en apnée 4,7,8. Pour la quantification, des images IRM anatomiques 1H appariées au volume sont généralement acquises en même temps pour être co-enregistrées avec 129Xe afin de délimiter les limites de la cavité thoracique. L’IRM conventionnelle de 1H, cependant, ne fournit pas d’informations supplémentaires sur la structure pulmonaire. L’impulsion pour la traduction clinique de l’IRM hyperpolarisée 129Xe s’est accrue ces dernières années avec l’approbation du NHS britannique en 2015 et l’approbation de la FDA américaine fin 2022 5,9, mais la caractérisation structurelle avancée est encore largement absente de l’arsenal de l’IRM pulmonaire.
La tomodensitométrie thoracique (TDM) reste le pilier de l’évaluation de l’imagerie clinique des poumons, fournissant des images tridimensionnelles à haute résolution de la structure pulmonaire à l’aide de protocoles d’imagerie conventionnels. Des approches quantitatives ont permis de mesurer rapidement et reproductiblement l’intégrité du parenchyme, comme l’emphysème et les anomalies pulmonaires interstitielles, la morphologie des grandes voies respiratoires et la vascularisation pulmonaire, ainsi que la caractérisation anatomique régionale par identification et segmentation des lobes pulmonaires10,11. Dans le domaine de la recherche, la TDM quantitative a été largement utilisée pour mieux comprendre les altérations structurelles et leurs relations avec les résultats des patients atteints d’asthme et de MPOC dans le cadre de vastes études observationnelles telles que le Programme de recherche sur l’asthme sévère (SARP)12, l’épidémiologie génétique de la MPOC (COPDGene)13, l’étude sur les sous-populations et les résultats intermédiaires dans l’étude sur la MPOC (SPIROMICS)14, l’évaluation longitudinale de la MPOC pour identifier les paramètres de substitution prédictifs (ECLIPSE)15et la Cohorte canadienne de pneumopathie obstructive (CanCOLD)16. D’autres méthodes de TDM telles que l’imagerie expiratoire17,18 ou les modèles informatiques19 peuvent dériver des informations fonctionnelles, mais ces méthodes sont indirectes et la TDM conventionnelle ne fournit pas grand-chose pour la caractérisation fonctionnelle des poumons.
Ensemble, la TDM et l’IRM 129Xe fournissent des informations complémentaires sur la structure et la fonction pulmonaires qui peuvent être exploitées pour l’analyse régionale à l’aide de l’enregistrement d’images. Les lobes pulmonaires identifiés par TDM ont permis de caractériser lobaire les schémas de ventilation par IRM dans l’asthme20,21,22, la BPCO23,24, la bronchectasie25 et le cancer du poumon26,27. Les anomalies de la ventilation par IRM dans l’asthme ont également été directement appariées spatialement à de grandes voies respiratoires anormalement remodelées 28,29,30,31 et à un piégeage de l’air indiquant un dysfonctionnement des petites voies respiratoires 20,32 mesurées par tomodensitométrie, et pour sonder les réponses de traitement régionales après une thermoplastie bronchique du poumon entier 33. Dans la BPCO, les anomalies de la ventilation par IRM ont été associées à un dysfonctionnement des petites voies respiratoires dans les maladies plus bénignes et à l’emphysème dans les maladies plus graves 34,35,36. Au-delà de l’imagerie de ventilation dans la pneumopathie obstructive, des relations spatiales hétérogènes entre les anomalies pulmonaires interstitielles CT et les modèles d’échange gazeux 129Xe MRI ont également été mises en évidence dans la fibrose pulmonaire idiopathique37. De telles études ont permis de mieux comprendre la structure et la fonction pulmonaires régionales dans une gamme de maladies pulmonaires qui peuvent être utilisées pour éclairer de futures interventions guidées par l’image.
L’enregistrement direct de la tomodensitométrie anatomique et de l’IRM fonctionnelle en gaz hyperpolarisé est toutefois difficile en raison du contraste d’imagerie fondamentalement différent entre les deux méthodes, de l’absence de signal de gaz hyperpolarisé dans les régions présentant des anomalies de la ventilation et des volumes pulmonaires potentiellement différents. La figure 1 montre quatre exemples d’IRM et de TDM 129Xe et anatomiques 1 Happariées chez un volontaire en bonne santé (figure 1A) et trois participants atteints de bronchopneumopathie chronique obstructive (BPCO ; Figure 1B-D), mettant en évidence des modèles de ventilation 129Xe hétérogènes et des limites pulmonaires manquantes variables dans les cas de BPCO. La clé pour surmonter ces défis a été d’utiliser l’IRM anatomique 1H acquise en même temps que l’IRM en gaz hyperpolarisé comme étape intermédiaire pour enregistrer l’IRM à gaz hyperpolarisé à la TDM indirectement34,38. Les premiers travaux ont utilisé la comparaison visuelle côte à côte et la segmentation manuelle des structures de la tomodensitométrie, telles que les lobes pulmonaires, sur l’espace IRM20. Les progrès des ressources de calcul et des outils de traitement d’images open source ont permis l’enregistrement tridimensionnel de la TDM et de l’IRM des gaz hyperpolarisés, par exemple, à l’aide du descripteur de voisinage indépendant de la modalité (MIND)23,30,34,39,40,41 ou de l’enregistrement de la boîte à outils de normalisation avancée (ANTs) 21,22,27 ,31,32,37,38,42,43, qui ont tous deux été les plus performants dans un défi d’enregistrement d’images pulmonaires44. Une nouvelle méthode a couplé les deux enregistrements plutôt que de les traiter indépendamment45, qui a été mise en œuvre dans un pipeline complet d’analyse d’images pulmonaires conçu pour le phénotypage des maladies pulmonaires46. Dans l’ensemble, la précision de l’enregistrement de l’IRM en TDM en gaz hyperpolarisé a été améliorée en utilisant l’étape intermédiaire 1H38 et en utilisant des approches déformables par rapport aux approches affines uniquement38,45.
L’objectif ici est de s’appuyer sur la littérature existante et de fournir un protocole pour le recalage d’images 129Xe MR à CT à l’aide de plates-formes open source 47,48,49. Le protocole est mis en œuvre à l’aide d’ANTsPy et, conformément aux travaux précédents38, enregistre un masque pulmonaire à étiquetteunique à partir de l’IRM 1 H au masque pulmonaire à étiquette unique de la TDM ; la transformation résultante est ensuite appliquée à l’image 129Xe pour la mapper à l’espace d’image CT. Le protocole décrit est destiné à être approprié pour la recherche ou les milieux cliniques, le cas échéant, et l’IRM hyperpolarisée 129Xe est disponible.
Pour mettre en contexte, l’acquisition et l’analyse d’images pour les exemples fournis ici ont été effectuées comme suit. La tomodensitométrie thoracique a été acquise à pleine inspiration (capacité pulmonaire totale, CCM) selon un protocole de recherche à faible doseétabli 50 avec des paramètres : collimation 64 x 0,625, kilotension de crête de 120, courant de tube 100 mA, temps de révolution de 0,5 s, pas de spirale de 1,0, épaisseur de coupe de 1,25 mm, espacement de tranche de 0,80 mm, noyau de reconstruction standard, champ de vision d’affichage limité aux étendues les plus latérales des poumons (pour maximiser la résolution spatiale). La segmentation et l’analyse de la tomodensitométrie ont été effectuées à l’aide d’un logiciel commercial (voir la table des matériaux).
129L’IRM Xe et l’IRM 1H appariée au volume ont été réalisées conformément aux directives publiées9. Pour plus de détails et le protocole d’acquisition de l’IRM, les lecteurs sont dirigés vers un autre article de cette collection51. La segmentation et l’enregistrement de l’IRM ont été effectués à l’aide d’un pipeline personnalisé semi-automatisé utilisant le clustering k-means pour la segmentation 129Xe, la croissance de la région ensemencée pour la segmentation 1H et l’enregistrement affine basé sur des points de repère pour mapper l’image 1H à l’image 129Xe52. L’enregistrement affine est généralement suffisant pour l’enregistrement1 H-129Xe MR pour tenir compte de la plupart des gonflements pulmonaires ou des différences de position du patient entre les acquisitions ; L’enregistrement déformable n’est généralement pas nécessaire. L’étape d’enregistrement 1 H-129Xe peut être éliminée avec l’IRM 129Xe et 1H acquises simultanément dans la même apnée53,54.