Un pipeline standardisé pour l’examen de la morphométrie de la matière grise cérébelleuse humaine à l’aide de l’imagerie par résonance magnétique structurelle

Résumé

Un pipeline normalisé est présenté pour l’examen de la morphométrie de la matière grise du cervelet. Le pipeline combine des approches de pointe à haute résolution pour une parcellation optimisée et automatisée du cervelet et l’enregistrement du cervelet à base de voxel pour la quantification volumétrique.

Résumé

De multiples axes de recherche fournissent des preuves convaincantes du rôle du cervelet dans un large éventail de fonctions cognitives et affectives, allant bien au-delà de son association historique avec le contrôle moteur. Les études de neuroimagerie structurale et fonctionnelle ont affiné la compréhension de la neuroanatomie fonctionnelle du cervelet au-delà de ses divisions anatomiques, soulignant la nécessité d’examiner les sous-unités cérébelleuses individuelles dans la variabilité saine et les maladies neurologiques. Cet article présente un pipeline standardisé pour l’examen de la morphométrie de la matière grise du cervelet qui combine des approches de pointe à haute résolution pour une parcellation optimisée et automatisée du cervelet (Automatic Cerebellum Anatomical Parcellation using U-Net Locally Constrained Optimization; ACAPULCO) et l’enregistrement du cervelet à base de voxel (modèle infra-tentorial spatialement non biaisé; SUIT) pour la quantification volumétrique.

Le pipeline a une large applicabilité à une gamme de maladies neurologiques et est entièrement automatisé, avec une intervention manuelle uniquement nécessaire pour le contrôle de la qualité des résultats. Le pipeline est disponible gratuitement, accompagné d’une documentation d’accompagnement substantielle, et peut être exécuté sur les systèmes d’exploitation Mac, Windows et Linux. Le pipeline est appliqué dans une cohorte de personnes atteintes d’ataxie de Friedreich (FRDA), et des résultats représentatifs, ainsi que des recommandations sur les analyses statistiques inférentielles au niveau du groupe, sont fournis. Ce pipeline pourrait faciliter la fiabilité et la reproductibilité dans l’ensemble du domaine, fournissant ainsi une approche méthodologique puissante pour caractériser et suivre les changements structurels cérébelleux dans les maladies neurologiques.

Introduction

Le cervelet est une partie du cerveau historiquement associée au contrôle moteur ^1,2,3 et on pense qu’il n’est pleinement impliqué que dans un petit ensemble de maladies rares, telles que les ataxies héréditaires⁴. Cependant, les lignes de recherche convergentes issues d’études de traçage anatomique chez des primates non humains, ainsi que d’études sur les lésions humaines et la neuroimagerie, fournissent des preuves convaincantes du rôle du cervelet dans un large éventail de fonctions cognitives ^5,6,7, affectives ^8,9,10,11 et d’autres fonctions non ^{motrices 7,12} (voir 6  pour examen). En outre, les anomalies du cervelet sont de plus en plus impliquées dans un large éventail de troubles neurologiques et psychiatriques, y compris la maladie de Parkinson¹³, la maladie d’Alzheimer^14,15, l’épilepsie^16,17, la schizophrénie¹⁸ et les troubles du spectre autistique¹⁹ . Par conséquent, il est devenu essentiel d’incorporer le cervelet dans les modèles fonctionnels et structurels des maladies du cerveau humain et de la variabilité comportementale normative.

Anatomiquement, le cervelet peut être divisé le long de son axe supérieur à inférieur en trois lobes: antérieur, postérieur et floculonodulaire. Les lobes sont subdivisés en 10 lobules désignés par des chiffres romains I-X^20,21 (Figure 1). Le cervelet peut également être regroupé en zones médianes (vermis) et latérales (hémisphère), qui reçoivent respectivement des entrées de la moelle épinière et du cortex cérébral. Le lobe antérieur, composé de lobules I-V, a traditionnellement été associé à des processus moteurs et présente des connexions réciproques avec les cortex moteurs cérébraux²². Le lobe postérieur, composé de lobules VI-IX, est principalement associé aux processus non moteurs¹¹ et présente des connexions réciproques avec le cortex préfrontal, le cortex pariétal postérieur et le cortex cérébral temporal supérieur ^8,23. Enfin, le lobe floculonodulaire, comprenant le lobule X, a des connexions réciproques avec les noyaux vestibulaires qui régissent les mouvements oculaires et l’équilibre du corps pendant la posture et la démarche²¹.

Un nombre croissant de travaux récents utilisant la neuroimagerie fonctionnelle a affiné la compréhension de la neuroanatomie fonctionnelle du cervelet au-delà de ses divisions anatomiques. Par exemple, des techniques d’imagerie par résonance magnétique fonctionnelle (IRMf) à l’état de repos ont été utilisées pour cartographier le modèle d’interactions fonctionnelles entre le cervelet et le cerveau²⁴. De plus, en utilisant une approche de parcellation basée sur les tâches, King et ses collègues⁷ ont démontré que le cervelet présente un modèle riche et complexe de spécialisation fonctionnelle dans toute son étendue, mis en évidence par des limites fonctionnelles distinctes associées à une variété de tâches motrices, affectives, sociales et cognitives. Collectivement, ces études soulignent l’importance d’examiner les sous-unités cérébelleuses individuelles afin de développer des caractérisations biologiques complètes de l’implication du cervelet dans la variabilité saine et les maladies neurologiques caractérisées par des altérations de la structure et/ou de la fonction cérébelleuses.

Le présent travail se concentre sur les méthodes de quantification des changements locaux dans le volume cérébelleux à l’aide de l’IRM structurelle chez l’homme. En général, il existe deux approches fondamentales pour quantifier le volume cérébral régional à l’aide de données IRM : la segmentation basée sur les caractéristiques et l’enregistrement basé sur le voxel. Les approches de segmentation basées sur les caractéristiques utilisent des repères anatomiques et des atlas standardisés pour identifier automatiquement les frontières entre les sous-régions. Les progiciels grand public pour la segmentation incluent FreeSurfer²⁵, BrainSuite²⁶ et FSL-FIRST²⁷. Cependant, ces emballages ne fournissent que des parcellations grossières du cervelet (par exemple, étiqueter toute la matière grise et toute la substance blanche dans chaque hémisphère), négligeant ainsi les lobules cérébelleux individuels. Ces approches sont également sujettes à une mauvaise segmentation, en particulier à la surinclusion du système vasculaire environnant.

De nouveaux algorithmes d’apprentissage automatique et d’étiquetage multi-atlas ont été développés, qui fournissent une parcellation plus précise et plus fine du cervelet, y compris la classification automatique de l’algorithme de lobules cérébelleux à l’aide de l’évolution multilimite implicite (ACCLAIM ^28,29), la boîte à outils d’analyse cérébelleuse (CATK³⁰), les modèles générés automatiquement multiples (MAGeT³¹), la segmentation automatique rapide du cervelet humain et de ses lobules (RASCAL³² ), segmentation graphique³³ et segmentation CEREbellum (CERES³⁴). Dans un article récent comparant les approches de parcellation du cervelet entièrement automatisées de pointe, CERES2 s’est avéré surpasser les autres approches par rapport à la segmentation manuelle de référence des lobules cérébelleux³⁵. Plus récemment, Han et ses collègues^{de 36} ont développé un algorithme d’apprentissage profond appelé ACAPULCO (Automatic Cerebellum Anatomical Parcellation using U-Net with locally constrained optimization), qui fonctionne sur un pied d’égalité avec CERES2, a une large applicabilité aux cervelets sains et atrophiés, est disponible en format de conteneur Docker et Singularity open-source pour une mise en œuvre « prête à l’emploi » et est plus rapide que d’autres approches. ACAPULCO répartit automatiquement le cervelet en 28 régions anatomiques.

Contrairement à la segmentation basée sur les fonctionnalités, les approches d’enregistrement basées sur le voxel fonctionnent en mappant précisément une IRM à une image modèle. Pour réaliser ce mappage, les voxels de l’image originale doivent être déformés en taille et en forme. L’ampleur de cette distorsion fournit effectivement une mesure du volume à chaque voxel par rapport au modèle d’étalon-or. Cette forme d’évaluation volumétrique est connue sous le nom de « morphométrie à base de voxel »³⁷. Les approches d’enregistrement basées sur le voxel du cerveau entier, telles que FSL-FLIRT³⁸/FNIRT³⁹, SPM unified segmentation⁴⁰ et CAT12⁴¹, sont couramment utilisées pour la morphométrie à base de voxel. Cependant, ces approches ne tiennent pas bien compte du cervelet, ce qui entraîne une fiabilité et une validité médiocres dans les régions infratentorielles (cervelet, tronc cérébral⁴²). Pour tenir compte de ces limitations, l’algorithme SUIT (Spatially Unbiased Infra-tentorial Template) a été développé pour optimiser l’enregistrement du cervelet et améliorer la précision de la morphométrie à base de voxel^42,43.

La segmentation basée sur les caractéristiques et les approches d’enregistrement basées sur le voxel pour l’estimation du volume cérébelleux régional ont des forces et des faiblesses fondamentales. Les approches de segmentation sont beaucoup plus précises pour quantifier le volume des zones anatomiquement définies (p. ex., lobules³⁵). Cependant, les limites entre les modules fonctionnels distincts du cervelet ne correspondent pas à ses folias anatomiques et à ses fissures (équivalents aux gyri et sulci du cerveau⁷). Comme les approches basées sur l’enregistrement ne sont pas limitées par des repères anatomiques, une inférence spatiale plus fine et une cartographie structure-fonction à haute dimension du cervelet sont^{possibles 44}. Prises ensemble, les approches de segmentation et d’enregistrement sont complémentaires et peuvent être utilisées pour répondre à différentes questions de recherche.

Ici, un nouveau pipeline standardisé est présenté, qui intègre ces approches existantes et validées pour fournir une parcellation optimisée et automatisée (ACAPULCO) et un enregistrement du cervelet à base de voxel (SUIT) pour la quantification volumétrique (Figure 2). Le pipeline s’appuie sur les approches établies pour inclure des protocoles de contrôle de la qualité, en utilisant la visualisation qualitative et la détection quantitative des valeurs aberrantes, ainsi qu’une méthode rapide pour obtenir une estimation du volume intracrânien (ICV) à l’aide de Freesurfer. Le pipeline est entièrement automatisé, avec une intervention manuelle uniquement nécessaire pour vérifier les sorties de contrôle qualité, et peut être exécuté sur les systèmes d’exploitation Mac, Windows et Linux. Le pipeline est disponible gratuitement sans restriction quant à son utilisation à des fins non commerciales et peut être consulté à partir de la page Web ENIGMA Consortium Imaging Protocols (sous « ENIGMA Cerebellum Volumetrics Pipeline »), après avoir rempli un bref formulaire d’inscription⁴⁵.

Tous les logiciels requis sont répertoriés dans la table des matériaux, et des didacticiels détaillés, y compris une démonstration en direct, sont disponibles lors du téléchargement du pipeline, en plus du protocole décrit ci-dessous. Enfin, des résultats représentatifs sont fournis, issus de la mise en œuvre du pipeline dans une cohorte de personnes atteintes d’ataxie de Friedreich (FRDA) et de témoins sains appariés selon l’âge et le sexe, ainsi que des recommandations pour des analyses inférentielles statistiques au niveau du groupe.

Protocole

REMARQUE : Les données utilisées dans cette étude faisaient partie d’un projet approuvé par le Comité d’éthique de la recherche humaine de l’Université Monash (projet 7810). Les participants ont donné leur consentement éclairé par écrit. Alors que le pipeline peut être exécuté sur les systèmes d’exploitation Mac, Windows ou Linux, ACAPULCO, SUIT et les pipelines QC ont été explicitement testés sur les systèmes d’exploitation Linux (Ubuntu) et Mac (Catalina, Big Sur v11.0.1).

1. Module 1 : ACAPULCO (parcellation anatomique)

1. Collecte de données
   1. Recueillez des images IRM pondérées en T1 3D de l’ensemble du cerveau à une résolution de 1 mm³ ou moins. Les dimensions isotropes du voxel (généralement 1 mm x 1 mm x 1 mm) et un scanner 3-Telsa (ou supérieur) sont recommandés. Consultez un spécialiste en imagerie de son centre de radiographie pour mettre en place et acquérir des données qui répondent à ces spécifications.
      REMARQUE: Les images pondérées en T2 sont parfois utiles pour les analyses volumétriques; toutefois, le pipeline présenté ici repose uniquement sur des données pondérées T1, et certains des outils utilisés sont exclusifs à ce type de données. En tant que telles, les images pondérées en T2 ne peuvent pas être utilisées.
   2. Entreprendre une évaluation de la qualité visuelle des images afin d’exclure les malformations cérébelleuses grossières (p. ex., les grandes lésions) ou les artefacts de mouvement importants qui empêchent l’identification des principaux repères cérébelleux (p. ex., les principales fissures anatomiques). N’excluez pas automatiquement la cervelet atrophié, même s’il est substantiel.
   3. Pour les études de groupe, envisagez également des évaluations quantitatives de la qualité à l’aide d’outils normalisés et disponibles gratuitement, tels que MRIQC⁴⁶, afin d’identifier davantage les données problématiques.
   4. Convertissez toutes les données au format NIFTI-GZ à l’aide d’un outil tel que dcm2niix⁴⁷.
2. Organisation des données recommandée
   1. Procurez-vous tous les logiciels nécessaires comme indiqué dans le tableau des matériaux. Assurez-vous que Docker⁴⁸ ou Singularity⁴⁹, Matlab⁵⁰ et SPM12⁵¹ sont installés avant d’exécuter le pipeline.
      REMARQUE : Des didacticiels écrits et vidéo détaillés décrivant le pipeline sont également disponibles (voir le Tableau des matériaux).
   2. Une fois tous les logiciels nécessaires installés, créez des dossiers dans le répertoire de travail et étiquetez-les « acapulco », « suit » et « freesurfer ». Pour ce faire, utilisez la commande mkdir de la ligne de commande.
   3. Dans le répertoire 'acapulco', créez un dossier de sortie . Dans le dossier de sortie , créez un répertoire pour chaque sujet de l’étude contenant l’image pondérée T1 au format NIFTI-GZ.
      REMARQUE: Il est recommandé de conserver une copie des données originales ailleurs.
3. Parcellation cérébelleuse anatomique à l’aide d’ACAPULCO
   1. Accédez à la table des matériaux et téléchargez les scripts et conteneurs pertinents requis pour exécuter ACAPULCO (sous les fichiers de pipeline acapulco). Dans le répertoire 'acapulco', placez le conteneur (i) ACAPULCO Docker OR Singularity ('acapulco_0.2.1.tar.gz' ou '.sif', respectivement), (ii) le contenu de l’archive QC_scripts (3 fichiers: 'QC_Master.R', 'QC_Plots.Rmd' et 'QC_Image_Merge.Rmd'), et (iii) 'R.sif' (singularité) OU 'calculate_icv.tar' (docker).
   2. Ouvrez un terminal et, à partir de la ligne de commande, exécutez le conteneur ACAPULCO sur une seule image (remplacez <> dans ce qui suit). Attendez environ 5 minutes pour que le traitement se termine.
      1. À l’aide de Docker, tapez la commande :
         docker load --input acapulco_0.2.1.tar.gz
         docker run -v $PWD:$PWD -w $PWD -t --user $(id -u):$(id -g) --rm acapulco:latest -i output/<>/<>.nii.gz -o output/<>
      2. À l’aide de Singularité, tapez la commande :
         singularity run --cleanenv -B $PWD:$PWD acapulco-0.2.1.sif -i output/<>/<>.nii.gz -o output/<>
   3. Boucle à travers tous les sujets / scans de la cohorte. Consultez la table des matériaux pour un lien vers le site Web ENIGMA Imaging Protocols pour télécharger le pipeline (sous ENIGMA Cerebellum Volumetrics Pipeline) et le manuel du didacticiel contenant des exemples de création d’une boucle for pour le traitement en série de plusieurs sujets.
   4. Après le traitement, recherchez les fichiers suivants générés dans les dossiers spécifiques à l’objet :
      1. Identifiez « _n4_mni_seg_post_inverse.nii.gz »: masque de cervelet parcellé dans l’original (espace sujet).
      2. Identifier « _n4_mni_seg_post_volumes.csv »: volumes (en mm³) pour chacune des 28 sous-unités générées par l’acapulco;
      3. Identifiez les images représentatives (dans le répertoire 'pics'): sagittales, axiales et coronales.
4. Détection statistique des valeurs aberrantes et contrôle de la qualité (CQ)
   1. Depuis le terminal et dans le répertoire 'acapulco', assurez-vous que le contenu de QC_scripts se trouve dans le répertoire 'acapulco'. Pour exécuter les scripts QC :
      1. À l’aide de Docker, tapez la commande :
         calculate_icv.tar de chargement docker
         docker run -v $PWD:$PWD -w $PWD --rm -it luhancheng/calculate_icv:latest Rscript
         sortie QC_Master.R/
      2. À l’aide de Singularité, tapez la commande :
         singularité exec -B $PWD:$PWD R.sif Rscript /path/to/QC_Master.R /path/to/acapulco/output
5. Examen des images QC générées par ACAPULCO
   REMARQUE: Il existe un processus en 3 étapes pour vérifier la qualité des images parcellées ACAPULCO.
   1. Ouvrez le 'QC_Images.html' dans un navigateur Web et faites défiler rapidement (~ 10 s par sujet) les images pour identifier les échecs évidents ou les problèmes systématiques. Notez les ID d’objet des images parcellées défaillantes ou suspectes pour le suivi.
      REMARQUE: Voir la figure 3 pour un guide sur la neuroanatomie des lobules cérébelleux et la figure 4, la figure 5 et la figure 6 dans la section des résultats représentatifs ci-dessous pour des exemples de « bonnes » parcellations, de « mauvaises parcellations subtiles » et de « défaillance globale ».
   2. Ouvrez le 'Plots_for_Outliers.html' pour cocher les cases des tracés pour les valeurs statistiques quantitatives aberrantes. Recherchez les valeurs aberrantes (2,698 s.d au-dessus ou en dessous de la moyenne) au-dessus ou au-dessous des moustaches des diagrammes de boîte. Passez le curseur de la souris sur les points de données pour afficher l’ID d’objet. Identifiez les valeurs aberrantes signalées par un « 1 » dans la colonne correspondante du fichier « Valeurs aberrantes .csv » et notez le nombre total de segments identifiés comme valeurs aberrantes pour chaque sujet dans la dernière colonne de « Valeurs aberrantes .csv ».
   3. Inspectez manuellement chaque image ayant une ou plusieurs valeurs aberrantes. CRITIQUE : À l’aide d’une visionneuse d’images NIFTI standard (par exemple, FSLEyes ou MRICron), superposez le masque ACAPULCO sur l’image T1w d’origine pour vérifier la qualité de la parcellation tranche par tranche.
      1. Pour générer des superpositions pour le QC détaillé à partir de la ligne de commande à l’aide de FSLEyes, i) remplacez le répertoire par le répertoire 'acapulco', ii) spécifiez l’objet à afficher (remplacez ) :
         subj=
      2. Copiez/collez le code suivant sur le terminal (sans modifier manuellement {subj} car il a été défini par la ligne précédente :
         t1_image=output/${subj}/${subj}.nii.gz
         acapulco_image=output/${subj}/${subj}_n4_mni_seg_post_inverse.nii.gz
         fsleyes ${t1_image} ${acapulco_image} --overlayType label --lut random_big --outline --outlineWidth 3 ${acapulco_image} --overlayType volume --alpha 50 --cmap aléatoire
         REMARQUE: Il faudra déterminer s’il faut inclure le segment anormal ou non, c’est-à-dire s’il y a une erreur de parcellation ou s’il s’agit simplement d’une variabilité normale de l’anatomie de l’individu? Chaque région parcellée est considérée individuellement, de sorte que quelques régions peuvent être exclues pour une image, tandis que le reste peut être conservé si elle est correcte.
      3. Une ou plusieurs régions parcellées doivent-elles être exclues de l’ensemble de données final ?
         Si oui (valeur aberrante confirmée), exclure cette (ces) parcellation(s) de l’analyse en remplaçant l’estimation du volume par NA dans la cellule correspondante du fichier « Cerebel_vols.csv » pour ce sujet.
      4. Les erreurs de parcellation entraînent-elles l’exclusion d’une partie du cervelet du masque?
         Si oui (par exemple, si des lobules cérébelleux particuliers sont absents du masque ou semblent « coupés »), excluez immédiatement le sujet des analyses ultérieures (c’est-à-dire ne pas exécuter le module SUIT sur ces sujets).

2. Module 2 : Morphométrie à base de voxel optimisée pour le cervelet SUIT

1. Analyses morphométriques à base de voxel à l’aide de SUIT
   CRITIQUE: Ce pipeline nécessite que le module ACAPULCO ait déjà été exécuté, car il repose sur la génération d’un masque cérébelleux spécifique au sujet pour optimiser l’enregistrement et la normalisation du cervelet dans le modèle SUIT. Si le masque spécifique au sujet généré par ACAPULCO n’inclut pas l’ensemble du cervelet, cela justifie l’exclusion du module SUIT. Pour obtenir des instructions sur l’exécution de SUIT autonome, reportez-vous à^{la section 52}.
   1. Procurez-vous tous les logiciels nécessaires répertoriés dans le tableau des matériaux. Assurez-vous que le dossier SPM12 et tous les sous-dossiers se trouvent dans le chemin d’accès MATLAB. Assurez-vous enigma_suit scripts sont enregistrés dans le répertoire 'spm12/toolbox' et ajoutés au chemin MATLAB. Pour vérifier le chemin d’accès MATLAB, tapez pathtool dans la fenêtre de commande MATLAB, puis cliquez sur Ajouter avec des sous-dossiers pour ajouter les dossiers appropriés.
   2. Exécutez le pipeline SUIT pour un ou plusieurs sujets. Attendez ~15-20 min (si vous utilisez l’interface utilisateur graphique [GUI]) et ~5-7 min si vous exécutez à partir du terminal (bash/shell) pour que le traitement se termine.
      1. Pour utiliser l’interface graphique (les sujets seront exécutés en série), dans la fenêtre de commande MATLAB, tapez la commande :
         suit_enigma_all
      2. Dans la première fenêtre contextuelle, sélectionnez les dossiers d’objet dans le répertoire 'acapulco/output' à inclure dans l’analyse. Cliquez sur les dossiers individuels sur le côté droit de la fenêtre ou cliquez avec le bouton droit de la souris et sélectionnez tout. Appuyez sur Terminé. Dans la deuxième fenêtre contextuelle, sélectionnez le répertoire SUIT, où les analyses seront écrites.
      3. OU Appelez la fonction à partir de la ligne de commande MATLAB pour un seul sujet, tapez la commande :
         suit_enigma_all('/path/to/acapulco/output/subjdir','/path/to/suitoutputdir')
      4. OU Appelez la fonction depuis la fenêtre du terminal, en dehors de MATLAB, pour un seul sujet en tapant la commande :
         matlab -nodisplay -nosplash -r « suit_enigma_all('/path/to/acapulco/output/subjdir','/path/to/suitoutputdir'), exit »
   3. Consultez la table des matériaux pour un lien vers le site Web ENIGMA Imaging Protocols pour télécharger le pipeline (sous ENIGMA Cerebellum Volumetrics Pipeline) et le manuel du didacticiel contenant des exemples de création d’une boucle for pour le traitement en série de plusieurs sujets.
   4. Recherchez les points suivants concernant le script.
      1. Assurez-vous que le script copie l’image T1 corrigée des biais N4, alignée sur MNI (corps rigide) et le masque de cervelet ACAPULCO dans le répertoire de sortie.
      2. Assurez-vous que le script segmente la matière grise et blanche du cervelet.
      3. Assurez-vous que le script corrige les erreurs de surinclusion dans la parcellation à l’aide du masque ACAPULCO.
      4. Assurez-vous que le script DARTEL normalise et recoupe les données dans l’espace SUIT avec modulation jacobienne afin que la valeur de chaque voxel soit proportionnelle à son volume d’origine.
      5. Vérifiez dans le dossier de chaque sujet les résultats finaux suivants : 'wd_seg1.nii' (matière grise) et 'wd_seg2.nii' (substance blanche).
2. Détection statistique des valeurs aberrantes et contrôle de la qualité
   1. Inspectez visuellement les images normalisées et modulées (wd*) pour détecter les défaillances majeures. Dans MATLAB, tapez la commande :
      spm_display_4D
   2. Sélectionnez manuellement les images 'wd*seg1' dans les sous-dossiers de la combinaison ou accédez au répertoire 'suit'; insérez '^wd.*seg1' dans la zone Filtre (pas de guillemets) et appuyez sur le bouton Rec . Appuyez sur Terminé.
   3. Faites défiler les images pour vous assurer qu’elles sont toutes bien alignées. Voir la figure 7 pour des images correctement normalisées provenant de témoins sains (A, B) et d’une personne atteinte d’un cervelet fortement atrophique (D).
      REMARQUE: À ce stade, l’anatomie entre les sujets est très similaire (car ils ont été enregistrés dans le même modèle), et les différences de volume sont plutôt codées par des intensités de voxel différentes. Les défaillances majeures seront évidentes, p. ex. images vierges, grandes zones de tissu manquant, gradients d’intensité inhabituels (c.-à-d. voxels brillants tout en haut, voxels sombres tout en bas). Ces images doivent être exclues des étapes suivantes.
   4. Vérifiez la covariance spatiale pour les valeurs aberrantes. Dans MATLAB, tapez la commande :
      check_spatial_cov
      1. Sélectionnez les images 'wd*seg1' comme à l’étape précédente. Lorsque vous y êtes invité, sélectionnez les options suivantes : Mise à l’échelle des accessoires : Oui ; Variable à covary out: Non; Tranche (mm): - 48 , Espace: 1.
      2. Regardez le boxplot affichant la covariance spatiale moyenne de chaque image par rapport à toutes les autres de l’échantillon. Identifiez les points de données >2s.d. sous la moyenne dans la fenêtre de commande MATLAB. Pour ce faire, inspectez l’image « _n4_mni.nii.gz » dans le dossier SUIT à la recherche d’artefacts (mouvement, anomalies anatomiques), de problèmes de qualité d’image ou d’erreurs de prétraitement.
      3. Si la qualité de l’image et le prétraitement sont acceptables et que l’inspection visuelle des images modulées à l’étape précédente n’indique pas de problème de segmentation et de normalisation, conservez ces données dans l’échantillon. Sinon, excluez ces données.

3. MODULE 3 (facultatif) : Estimation du volume intracrânien (ICV) à l’aide de FreeSurfer

REMARQUE : Ce module utilisera le pipeline FreeSurfer pour calculer ICV. Il n’est pas nécessaire de le réexécuter s’il existe des sorties Freesurfer pour la cohorte (n’importe quelle version).

1. Configuration de FreeSurfer
   1. Assurez-vous que FreeSurfer est téléchargé et installé⁵³. Accédez à la table des matériaux et téléchargez les scripts pertinents pour exécuter ce module (sous fichiers de pipeline ICV). Lorsque vous travaillez avec FreeSurfer, définissez les variables suivantes :
      export FREESURFER_HOME=répertoire
      source $FREESURFER_HOME/SetUpFreeSurfer.sh
   2. Remplacez comme suit :
      export SUBJECTS_DIR=
2. Exécution de Freesurfer autorecon1
   1. Pour un seul sujet, à partir du répertoire 'freesurfer' (temps de traitement ~20 min), tapez la commande:
      cd
   2. Consultez le manuel du didacticiel pour obtenir des exemples de création d’une boucle for pour le traitement en série de plusieurs sujets.
3. Calcul de l’ICV
   1. Organisation des données
      1. Dans le répertoire 'freesurfer', placez le conteneur (i) Docker OR Singularity utilisé dans le module 1 ('calculate_icv.tar' ou 'R.sif', respectivement) et (ii) le script xfm2det (voir la table des matériaux). Ensuite, faites un clone git pour cloner le script ICV requis :
         https://github.com/Characterisation-Virtual-Laboratory/calculate_icv de clone git
   2. Extraction ICV en cours d’exécution (temps de traitement ~ 5 min)
      1. Dans le répertoire 'freesurfer', avec conteneur de singularité ('R.sif'), tapez :
         singularité exec --cleanenv -B $PWD:$PWD R.sif calculate_icv/calculate_icv.py --freesurfer_dir=/path/to/freesurfer --acapulco_dir=/path/to/acapulco/QC/Cerebelvolsfile --output_csv_name=Cerebel_vols.csv calculate_icv
      2. Dans le répertoire 'freesurfer', avec conteneur docker, tapez :
         docker run -v $PWD:$PWD -w $PWD -rm -it luhancheng/calculate_icv:latest
         calculate_icv/calculate_icv.py --freesurfer_dir=/chemin/vers/Freesurfer --
         acapulco_dir=/chemin/vers/acapulco/QC/Cerebelvolsfile --output_csv_name=Cerebel_vols.csv
         calculate_icv
      3. Exécution d’un script sans conteneur - consultez la table des matériaux pour plus de logiciels et de dépendances requis. Dans le répertoire 'freesurfer', tapez :
         ./calculate_icv/ calculate_icv.py ---freesurfer_dir=/chemin/vers/freesurfer --
         acapulco_dir=/chemin/vers/acapulco/QC/Cerebelvolsfile --
         output_csv_name=Cerebel_vols.csv calculate_icv
         REMARQUE: Cela calculera le ICV pour chaque sujet et ajoutera une colonne avec ICV à la fin du fichier 'Cerebel_vols.csv'.

Résultats

Parcellation du cervelet (ACAPULCO)

Contrôle de la qualité des masques parcellés au cervelet:
Les exemples suivants illustrent les extrants parcellés d’ACAPULCO et guident la prise de décision concernant a) la qualité du masque parcellé au niveau individuel et b) l’inclusion ou l’exclusion ultérieure d’un ou de plusieurs lobules particuliers des analyses statistiques. En fin de compte, la décision d’inclure ou d’exclure un sujet est sub...

Discussion

Le cervelet est essentiel à un large éventail de fonctions motrices humaines³,^{cognitives 58},^{affectives 10} et^{linguistiques 7,59} et est impliqué dans de nombreuses maladies neurologiques et psychiatriques. La disponibilité d’une approche normalisée et facilement réalisable pour la quantification des volumes cérébelleux régionaux contribuera à une cartographie de plus en plus détail...

matériels

Name	Company	Catalog Number	Comments
ACAPULCO pipeline files	0.2.1	http://enigma.ini.usc.edu/protocols/imaging-protocols/	Please make sure to use acapulco version 0.2.1
Docker for Mac		https://docs.docker.com/desktop/mac/install/	macOS must be version 10.14 or newer Docker requires sudo priviledges Docker imposes a memory (RAM) constraint on Mac OS. To increase the RAM, open Docker Desktop, go to Preferences and click on resources. Increase the Memory to the maximum
Docker for Windows		https://docs.docker.com/docker-for-windows/install/
ENIGMA SUIT scripts		http://enigma.ini.usc.edu/protocols/imaging-protocols/
FreeSurfer	7	https://surfer.nmr.mgh.harvard.edu/fswiki/DownloadAndInstall	Following variables need to be set everytime you work with Freesurfer: export FREESURFER_HOME=freesurfer _installation_directory source $FREESURFER_HOME/SetUpFreeSurfer.sh
export SUBJECTS_DIR=path/enigma/Freesurfer
FSL (for FSLeyes). Optional	6	https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FslInstallation
ICV pipeline files		http://enigma.ini.usc.edu/protocols/imaging-protocols/	ICV pipeline can be run in two ways: 1) with docker/singularity. You will not require additionl software; 2) without docker/singularity- this involves running the ICV script (calculate_icv.py) manually. You will require the following additional software: Python version =3.5 Python module pandas Python module fire Python module tabulate Python module Colorama
		https://github.com/Characterisation-Virtual-Laboratory/calculate_icv
MATLAB*	2019 or newer	https://au.mathworks.com/	An academic license is required
Singularity	3.7 or newer	https://www.sylabs.io/docs/	Prefered for high performance computing (HPC) clusters
SPM	12	http://www.fil.ion.ucl.ac.uk/spm/software/spm12/	Make sure spm12 and all subfolders are in your MATLAB path
SUIT Toolbox	3.4	http://www.diedrichsenlab.org/imaging/suit_download.htm	Make sure you place SUIT toolbox in spm12/toolbox directory
Troubleshooting manual and segmentation output examples		http://enigma.ini.usc.edu/protocols/imaging-protocols/
Tutorial manual and video		http://enigma.ini.usc.edu/protocols/imaging-protocols/	Manual and accompanying live demonstration provide detailed step-by-step instructions on how to run the pipeline from start to finish.
*Not freely available; an academic license is required