Évaluation basée sur l’élastographie de la réponse au traitement après une ablation par radiofréquence des nodules thyroïdiens

Résumé

Un protocole standardisé de suivi par ondes de cisaillement et élastographie de déformation a été appliqué avant et après l’ablation radiofréquence des nodules thyroïdiens afin de quantifier les variations longitudinales de rigidité. En utilisant la modélisation à effets mixtes et la régression multivariée, les changements précoces de l’élastographie post-ablation étaient indépendamment associés à une réponse de réduction de volume sur 12 mois.

Résumé

L’ablation par radiofréquence (RFA) est un traitement mini-invasif établi pour les nodules thyroïdiens bénins, mais l’évaluation de la réponse repose principalement sur une réduction retardée du volume, limitant la prise de décision post-procédurale précoce. Dans cette cohorte rétrospective, les participants présentant un seul nodule thyroïdien index ont subi une échographie standardisée en mode B et une élastographie au départ, immédiatement après l’ablation, et à 1, 3, 6 et 12 mois. La rigidité et le rapport de déformation de l’élastographie par ondes de cisaillement (SWE) ont été obtenus à l’aide de régions d’intérêt prédéfinies et de critères de contrôle qualité. La réponse à douze mois a été définie par un ratio de réduction de volume (VRR) ≥ 50 %. Les changements longitudinaux ont été évalués à l’aide de modèles linéaires à effets mixtes, et les modifications de l’élastographie immédiatement post-RFA (ligne de départ immédiate négative) ont été liées à la VRR à 12 mois en utilisant corrélation et régression multivariable ajustée au volume nodulique de base, à la composition des nodules, à l’énergie délivrée et à l’âge. Parmi les 187 participants inclus, 146 étaient des intervenants et 41 n’étaient pas intervenants. L’EE a augmenté immédiatement après l’RFA chez les répondants (31,5 ± 5,1 kPa à 41,7 ± 5,6 kPa) et les non-répondants (35,1 ± 5,3 kPa à 43,4 ± 5,8 kPa), suivie d’un déclin progressif avec une ESE à 12 mois plus faible chez les répondants que chez les non-répondants (21,6 ± 4,7 contre 31,2 ± 5,2 kPa). Le VRR s’est séparé tôt, avec un VRR sur 12 mois de 73,6 ± 6,9 % chez les répondants contre 38,5 ± 8,7 % chez les non-répondants. Dans les modèles à effets mixtes, l’augmentation immédiate de l’EAT était significative (estimation +10,22 kPa, P < 0,001), et une divergence tardive des groupes était évidente (non-répondant × estimation sur 12 mois +9,41 kPa, P < 0,001). Le changement précoce de l’élastographie était corrélé au VRR à 12 mois (rapport de déformation Δ : r = 0,445, P < 0,001 ; ΔSWE : r = 0,302, P < 0,001) et restait indépendamment associé au VRR dans les modèles ajustés (rapport de déformation Δ β = 14,80 % par +1,0, P < 0,001 ; ΔSWE β = 1,80 % par +5 kPa, P = 0,011). Les changements d’élastographie immédiatement post-RFA capturent les premiers effets tissulaires et fournissent des informations pronostiques indépendantes pour la réponse volumétrique à 12 mois.

Introduction

Les nodules thyroïdiens sont devenus une « maladie de détection » définissante dans la pratique endocrinienne moderne, non pas parce que leur biologie a brusquement changé, mais parce que l’imagerie a rendu la thyroïde extrêmement visible ^1,2,3. La prévalence mesurée dépend fortement de la manière dont les nodules sont recherchés : les synthèses épidémiologiques classiques estiment 2 % à 6 % par palpation, environ 19 % à 35 % par échographie, et jusqu’à environ 8 % à 65 % dans les séries d’autopsies, illustrant comment la sensibilité à la méthode modifie l’épidémiologie^{apparente 1}. Dans des contextes pauvres en iode et avec des transducteurs à haute fréquence, les enquêtes ultrasonores peuvent rapporter des taux de détection encore plus élevés ; par exemple, une étude basée sur la population utilisant la technologie 13 MHz a identifié des nodules chez 68 % des adultes, mettant en avant à la fois l’influence de l’équipement et la susceptibilité sous-jacente de la thyroïde aux changements nodulaires tout au long de la^{vie 2}. La Chine offre une vision particulièrement vivante de l’échelle car l’échographie est fréquemment intégrée aux examens de santé. Dans une cohorte nationale d’examen de santé comprenant 6 985 956 participants, la prévalence globale des nodules thyroïdiens était de 36,9 %, avec une prévalence standardisée par âge et sexe de 38,0 %, et une prédominance féminine^{marquée de 3}. Même lorsque la plupart des lésions détectées sont cliniquement indolentes, les conséquences globales sont substantielles, incluant des imageries répétées, des cascades procédurales, une vigilance accrue des symptômes et une tension persistante entre réassurance et intervention. Cette tension est particulièrement aiguë pour les nodules bénignes qui sont symptomatiques, cosmétiquement visibles ou manifestement dilatants, où « bénin » ne signifie pas « insignifiant », et où les résultats centrés sur le patient dépendent souvent de la capacité d’un traitement à atténuer les effets locaux tout en préservant la fonction et en évitant la morbidité chirurgicale.

Au cours de la dernière décennie, l’ablation thermique guidée par échographie est passée d’une adoption de niche à une alternative de plus en plus établie à la chirurgie pour certains nodules thyroïdiens, en particulier les nodules solides bénins ou principalement solides et, dans des contextes soigneusement choisis, les lésions^{récurrentes 4,5}. L’ablation par radiofréquence (RFA) est devenue la technique phare car elle associe la visualisation échographique en temps réel à une livraison d’énergie contrôlable et un profil de complication favorable aux mains^{expérimentées 4,5}. Au-delà de l’attrait évident d’éviter une cicatrice cervicale et d’une anesthésie générale dans de nombreux cas, l’ARF s’inscrit dans une philosophie d’épargne d’organes selon laquelle la lésion est remodelée tout en préservant la thyroïde. Les rapports contemporains des résultats ont convergé vers le ratio de réduction de volume (VRR), généralement évalué à 6 à 12 mois, avec le succès du traitement dans de nombreuses séries et analyses groupées définies selon des seuils tels que une réduction volumique d’au moins 50 % à 1 an⁶. Pourtant, la réduction de volume est intrinsèquement un substitut tardif. L’effet biologique de l’ablation commence immédiatement avec la dénaturation des protéines et la nécrose coagulative, puis évolue par inflammation, résorption et fibrose sur des semaines à⁶ mois. En conséquence, la première question cliniquement décisive, à savoir fournir des informations supplémentaires sur la réponse tissulaire précoce après l’ablation, ne peut souvent pas être répondue avec confiance en se basant uniquement^{sur la taille 4,6}. La vascularisation Doppler et l’échotexture en niveaux de gris ajoutent un contexte utile, mais elles restent des mesures indirectes de la transformation tissulaire, et peuvent être confondues dans la première fenêtre post-procédurale par une hyperémie réactive, un œdème ou un changement nécrotique hétérogène. Un signal précoce et objectif d’un traitement adéquat pourrait rationaliser l’intensité du suivi, aider à anticiper la probabilité d’un retard de la réponse volumique, et aider à identifier les patients présentant des schémas de réponse longitudinale moins favorables avant que des mois d’attente attentive ne se consolident en symptômes persistants ou en^{repousse 2,3}.

L’élastographie est captivante dans ce contexte car elle mesure une propriété physique que l’ablation est conçue pour changer. L’élastographie par ultrasons, englobant les techniques basées sur la déformation et les ondes de cisaillement, caractérise la raideur tissulaire de manière semi-quantitative par des motifs de déformation ou quantitativement par la vitesse de l’onde de cisaillement et dérive du module^{élastique 7}. En imagerie thyroïdienne, l’élastographie a été largement évaluée comme complément de la stratification du risque de cancer. Cependant, des tests rigoureux ont tempéré l’enthousiasme initial concernant sa valeur diagnostique incrémentale par rapport à l’échographie conventionnelle dans de larges populations cliniques, rappelant au domaine que la physique élégante ne garantit pas un gain clinique dans tous les^{contextes 8}. Le problème post-ablation est différent. Ici, la lésion cible est connue et traitée. La tâche n’est pas une discrimination histologique mais une évaluation de la réponse, et la transition biologique centrale est d’un tissu viable et déformable vers un composite nécrotique et fibrosique moins docile. Les premières investigations cliniques à travers les modalités thermiques soutiennent ce postulat mécaniste. Chez les patients subissant une ablation micro-ondes, le score basé sur l’élastographie a montré une altération structurelle mesurable après le traitement et a été proposé comme outil de suivi^{réalisable 9}. Dans les cohortes RFA, l’élastographie a été explorée comme méthode pour quantifier les variations de rigidité dans les nodules traités et les tissus de référence adjacents, offrant une fenêtre potentielle sur le remodelage induit par l’ablation qui pourrait précéder le retrait^{macroscopique 10}. Le travail technique sur la courbe d’apprentissage a également intégré l’élastographie de la déformation dans la caractérisation échographie autour de la procédure, reflétant une reconnaissance croissante, bien que toujours incohérente, que les schémas de rigidité peuvent fournir des informations au-delà de^{la taille 11}. Pris ensemble, ces études suggèrent que l’élastographie pourrait évoluer d’un accessoire diagnostique à un biomarqueur de réponse. Cependant, l’absence de standardisation des protocoles et la variabilité des méthodes d’acquisition restent des obstacles importants à une mise en œuvre clinique plus^{large 10,11}. Les rapports existants varient en termes de modalité d’élastographie, de paramètres d’acquisition, de stratégie de région d’intérêt, de sélection des tissus de référence et de calendrier d’évaluation, ce qui rend difficile la traduction des variations de rigidité en règles cliniquement robustes pouvant guider la prise de décision précoce entre opérateurs et plateformes.

Cette étude a été conçue pour combler ce manque translationnel en évaluant une approche basée sur l’élastographie pour évaluer la réponse au traitement après l’AF des nodules thyroïdiens, avec une attention explicite à la mise en œuvre pratique lors du suivi de routine. Le postulat central de la présente étude est que la réponse après l’RFA est fondamentalement un processus biomécanique qui commence par une lésion thermique et progresse par des changements structurels et de rigidité longitudinaux après l’ablation, et que les métriques de rigidité doivent donc évoluer de manière cohérente et interprétable sur le plan clinique. Il a été émis l’hypothèse que les indices dérivés de l’élastographie démontreraient des déplacements systématiques post-ablation cohérents avec l’effet du traitement, suivraient la trajectoire de réponse sur le suivi, et fourniraient des informations complémentaires aux critères conventionnels tels que le VRR et la vascularité Doppler, en particulier dans la période initiale où l’évaluation basée sur la taille est la moins^{informative 4,6,7}. La nouveauté de la présente étude réside dans la reformulation de l’élastographie, passant d’une caractéristique supplémentaire de la discrimination des cancers vers une lecture objective et quantifiable de la transformation tissulaire post-ablation, et dans l’alignement des mesures d’élastographie avec les questions cliniques qui guident la prise en charge après l’RFA, y compris l’évaluation longitudinale de la réponse au traitement et la caractérisation du changement de raideur post-ablation lors du suivi. Si l’élastographie peut opérationnaliser de manière fiable ces questions, elle pourrait raccourcir la boucle de rétroaction entre l’intervention et le résultat, permettant potentiellement une caractérisation plus précoce des schémas de réponse longitudinale après l’RFA tout en préservant les avantages mini-invasifs qui ont fait de l’RFA thyroïdien une option de plus en plus importante pour les patients présentant des nodules bénins mais lourds ^4,6.

Protocole

Déclaration d’éthique
Toutes les procédures impliquant des participants humains ont été menées conformément aux directives institutionnelles de l’hôpital de Xingtai et en conformité avec les principes de la Déclaration d’Helsinki et de ses amendements ultérieurs. L’étude a été examinée et approuvée par le Comité d’éthique de la recherche humaine de l’hôpital de Xingtai ; comme tous les participants ont fourni un consentement clinique écrit standard pour les examens RFA et échographiques dans le cadre des soins de routine, le comité a approuvé l’étude avec une dérogation à un consentement écrit supplémentaire à la recherche et a permis le consentement verbal pour la participation à l’étude lorsque nécessaire. La confidentialité des participants a été protégée par une gestion sécurisée des données à accès restreint, avec la désidentification de toutes les données de l’étude avant l’analyse.

Conception de l’étude et schéma général
Cette étude a été conçue comme une cohorte clinique rétrospective à centre unique pour évaluer la réponse au traitement basée sur l’élastographie après une RFA guidée par échographie des nodules thyroïdiens. Les patients éligibles consécutifs ayant subi une RFA dans notre établissement ont été dépistés et inscrits jusqu’à ce qu’un objectif prédéfini de 187 participants soit atteint. Pour garantir l’indépendance analytique et éviter le regroupement interne au patient, un seul nodule indice par participant a été inclus, défini a priori comme le nodule sélectionné pour le traitement lors de la séance d’ablation indice. Les patients atteints de maladies thyroïdiennes diffuses pouvant affecter considérablement les mesures de raideur thyroïdienne de fond, y compris la thyroïdite diffuse de Hashimoto avec atteinte parenchymale hétérogène marquée, la maladie de Basedow non contrôlée, la thyroïdite diffuse ou la goïde diffuse, ont été exclus. Une positivité isolée des autoanticorps thyroïdiens sans anomalie structurelle thyroïdienne diffuse à l’échographie n’a pas été considérée comme un critère d’exclusion.

L’étude suivait un flux de travail longitudinal standardisé. L’évaluation clinique de base et l’imagerie ont été réalisées dans les deux semaines précédant l’ablation et comprenaient une échographie conventionnelle, une évaluation Doppler et l’acquisition d’élastographie selon un protocole prédéfini. La RFA a été réalisée sous guidage par échographie en temps réel selon les normes institutionnelles. L’imagerie immédiatement post-ablation a été réalisée dans les 30 minutes suivant l’administration d’énergie pour documenter le succès technique et les changements aigus des tissus. Les visites de suivi programmées ont été effectuées 1, 3, 6 et 12 mois après l’ablation, avec une échographie et une élastographie répétées à chaque moment en utilisant les mêmes réglages d’acquisition et stratégies de mesure afin de minimiser la variabilité technique.

Calcul de la taille de l’échantillon et hypothèses statistiques
Le critère principal pour l’estimation de la taille de l’échantillon était le changement intra-nodule de l’élastographie quantitative à un mois après l’ablation par rapport à la ligne de base, exprimé par la différence de rigidité moyenne (par exemple, le module élastique en kPa) mesurée à l’aide d’une stratégie standardisée de la région d’intérêt. Comme l’analyse primaire compare les mesures appariées au sein du même nodule, la taille de l’échantillon a été estimée à l’aide d’un cadre de différence moyenne appariée. Soit Δ la variation moyenne attendue de la rigidité de la base à 1 mois, et soit σ_d l’écart-type des différences par paires. Pour un test bilatéral au niveau de signification α avec une puissance de 1−β, la taille d’échantillon requise est :

En utilisant α = 0,05 (à double face) et la puissance 1−β = 0,80, Z1−α/2 = 1,96 et Z1−β = 0,84. Sur la base de mesures pilotes et de la variabilité post-ablation publiée en élastographie thyroïdienne, cette étude a supposé β_d = 35 kPa et un changement conservateur cliniquement significatif Δ = 8 kPa. Cela donne une taille d’échantillon minimale requise d’au moins 150 participants pour l’analyse par paires primaire de la rigidité dérivée de l’élastographie avant et après RFA.

Évaluation clinique de référence avant l’ablation
Lors de l’inscription, un historique clinique standardisé a été obtenu, en mettant l’accent sur les symptômes liés aux nodules et les facteurs liés à la sécurité des procédures. L’âge, le sexe, les mesures anthropométriques et les comorbidités pertinentes, y compris les maladies cardiovasculaires et pulmonaires, le diabète sucré, la maladie rénale chronique et une chirurgie cervique ou une irradiation antérieure, ont été enregistrés. Les médicaments actuels ont été documentés, avec une attention particulière aux antiplaquettaires et aux anticoagulants, et tout ajustement péri-procédural prévu selon la politique institutionnelle a été enregistré. L’indication pour l’ARF a été confirmée, précisant si le principal facteur était des symptômes compressifs, une inquiétude esthétique, une croissance documentée lors de l’imagerie en série, ou une combinaison de ces facteurs.

La charge symptomatique de base a été établie à l’aide d’un instrument prédéfini. Des symptômes compressifs locaux tels que la dysphagie, la pression de la gorge, la sensation corporelle étrangère, la dyspnée et les troubles vocaux ont été évalués, et leur gravité a été documentée à l’aide d’une échelle d’évaluation numérique ou d’un score symptomatique adopté institutionnellement. Lorsque l’impact esthétique était un indicateur, l’apparence esthétique était évaluée à l’aide d’une échelle esthétique standardisée, et l’insatisfaction rapportée par les patients était enregistrée. Les antécédents thyroïdiens pertinents, y compris la durée de la conscience des nodules, les résultats antérieurs d’aspiration par aiguille fine, la stratégie de confirmation de la cytologie bénigne antérieure (échantillonnage simple versus répété selon les besoins), ainsi que toute thérapie médicale antérieure pour les maladies thyroïdiennes, ont été recueillies.

Des tests de laboratoire de base ont été réalisés selon la pratique clinique de routine, incluant l’hormone thyroïdienne stimulante et des paramètres supplémentaires de fonction thyroïdienne lorsque cela était nécessaire. Les rapports de cytologie et de pathologie confirmant le statut bénin du nodule indice, ainsi que toute documentation sur la stratification du risque par échographie si disponible, ont été examinés et enregistrés. Une évaluation de sécurité de base, incluant les signes vitaux et une évaluation ciblée des voies respiratoires lorsque le nodule jouxait la trachée ou causait des symptômes suggérant une compromission des voies respiratoires, a été établie. Le calendrier de suivi prévu a été documenté, et des instructions standardisées avant l’imagerie, y compris l’évitement de la déglutition lors de l’acquisition de l’élastographie, ont été renforcées afin d’améliorer la cohérence des mesures entre les visites.

Acquisition conventionnelle par échographie et mesure du volume nodulique
L’échographie en niveaux de gris a été réalisée avec le patient en position allongée et le cou légèrement étendu. Un préréglage d’imagerie standardisé était utilisé, et les réglages clés étaient maintenus constants sur toutes les visites, notamment la plage de fréquences, la profondeur, la mise au point, le gain et la plage dynamique. Le taux d’images, la plage d’échelle d’élastographie, la persistance et les réglages de lissage sont restés constants lors des examens de base et de suivi. Les examens par ultrasons ont été réalisés à l’aide d’un transducteur à réseau linéaire haute fréquence (5–14 MHz) équipé de fonctions d’onde de cisaillement et d’élastographie par déformation. Tous les examens d’échographie et d’élastographie ont été réalisés à l’aide de la même plateforme d’échographie et de la même version logicielle tout au long de l’étude afin de minimiser la variabilité intersystémique des mesures de raideur. Des vues transversales et longitudinales qui délimitaient clairement les marges des nodules et les repères anatomiques adjacents ont été acquises. Les diamètres orthogonaux maximaux du nodule indice — antéropostérieur (A), transverse (B) et longitudinal (C) — ont été mesurés. Des étriers électroniques étaient placés sur le bord extérieur du nodule sur les images figées, évitant l’inclusion de parenchyme ou d’auréole périnodulaire environnante. Toutes les mesures étaient enregistrées en millimètres et, lorsque cela était possible, chaque dimension était répétée et la valeur moyenne était utilisée. Le volume des nodules a été calculé à l’aide de la formule ellipsoïde :

La même approche de mesure a été répétée à chaque visite de suivi. Le VRR était calculé comme suit :

Acquisition de l’élastographie et contrôle qualité
L’élastographie a été réalisée immédiatement après l’échographie conventionnelle, avec le patient allongé et le cou étendu. Le patient a été prié d’éviter d’avaler et de suspendre brièvement la respiration lors de l’acquisition afin de minimiser les artefacts de mouvement. Le même mode et réglages d’élastographie étaient utilisés à chaque visite.

Pour l’élastographie par ondes de cisaillement, la boîte d’échantillonnage a été placée pour inclure la portion solide du nodule indice et, lorsque cela était possible, le tissu thyroïdien normal adjacent à une profondeur similaire. Les régions circulaires d’intérêt étaient positionnées préférentiellement dans la composante centrale solide du nodule tout en évitant la dégénérescence kystique, la macrocalcification, les régions de halo périphériques et les changements visibles de la cavité post-ablation. Après le placement par région d’intérêt, l’acquisition de l’image n’a été lancée que lorsque la carte colorimétrique de l’élastographie est restée visuellement stable pendant au moins 3 à 5 secondes. L’évaluation de la raideur périnodulaire n’a pas été systématiquement intégrée au protocole d’acquisition prédéfini car l’objectif principal de l’étude portait sur l’évolution longitudinale de la raideur au sein du nodule traité lui-même. Au moins trois mesures techniquement adéquates ont été obtenues, et la valeur médiane de rigidité a été enregistrée comme métrique principale ; L’intervalle interquartile était documenté lorsqu’il était disponible.

Pour l’élastographie par déformation, une compression douce et constante était appliquée, guidée par l’indicateur de qualité du système. La compression était appliquée manuellement à une fréquence lente et rythmique tout en surveillant l’indicateur de qualité en temps réel fourni par le système d’échographie. Une pression excessive de la sonde, provoquant une déformation visible de la capsule thyroïdienne ou une cartographie colorimétrique instable de l’élastographie, a été évitée. Une région de référence d’intérêt a été définie dans un tissu thyroïdien normal adjacent à une profondeur similaire, et le rapport de déformation a été calculé. Chaque séquence d’acquisition d’élastographie a été maintenue pendant environ 5 à 10 secondes afin d’assurer la stabilisation de l’image et un échantillonnage reproductible. Au moins trois mesures acceptables ont été obtenues et la valeur médiane a été enregistrée.

Les mesures n’étaient acceptées que lorsque la carte d’élastographie était stable, que le placement par région d’intérêt suivait les règles prédéfinies, et qu’il n’y avait pas d’artefact de mouvement brut ou d’artefact de compression. Les mesures avec un remplissage de couleur incomplet, des vides évidents de signal dans la région d’intérêt, ou des fluctuations marquées d’un référentiel à l’autre ont été rejetées et reprises. Des examens par échographie et élastographie ont été réalisés dans une salle d’examen à température contrôlée afin de minimiser les influences environnementales sur les mesures de raideur tissulaire. Des images d’élastographie représentatives avec les images en mode B correspondantes étaient stockées à chaque moment pour la traçabilité et la qualité de contrôle.

Procédure RFA
L’RFA a été réalisée sous échographie continue par des opérateurs ayant suivi une formation institutionnelle en techniques d’ablation thyroïdienne, ayant plus de 5 ans d’expérience en interventions guidées par échographie thyroïdienne, et ayant chacun réalisé plus de 200 interventions RFA thyroïdiennes avant le début de l’étude. Le patient était positionné sur le dos, le cou étendu. La région cervicale a été préparée et drapée selon la technique stérile standard. Une anesthésie locale était administrée selon la pratique courante ; La sédation consciente était utilisée lorsqu’elle était cliniquement indiquée, avec un suivi continu des signes vitaux.

L’électrode à radiofréquence a été introduite dans le nodule indice sous guidage par échographie en temps réel, utilisant de préférence une approche transisthmique pour améliorer la stabilité de l’électrode et minimiser les déplacements involontaires. Lorsque le nodule indice était adjacent à la trachée, à la gaine carotide, à l’œsophage ou à la région du nerf laryngé récurrent, l’hydrodissection était réalisée par injection lente de liquide stérile sous guidage par échographie afin de maintenir une séparation visible entre la zone d’ablation et les structures critiques adjacentes tout au long de la délivrance d’énergie.

L’énergie radiofréquence était délivrée à l’aide d’une technique de tir mobile, partant de la partie la plus profonde du nodule et progressant systématiquement vers les régions superficielles pour obtenir des zones d’ablation qui se chevauchent. Chaque unité d’ablation a été traitée jusqu’à ce qu’un changement hyperéchoïque transitoire occupe la zone ciblée avant que la pointe de l’électrode ne soit repositionnée sur la partie adjacente non traitée. La puissance de sortie a été ajustée en réponse aux changements échogéniques en temps réel et à la tolérance du patient, et l’ablation a été poursuivie jusqu’à ce que le nodule ciblé montre le changement hyperéchoïque transitoire attendu et que la vascularisation intranodulaire ait été nettement réduite lors de l’évaluation Doppler. Une ablation adéquate a également été soutenue par la visualisation d’une zone hyperéchoïque continue couvrant l’unité de traitement ciblée sans signaux vasculaires intranodulaires résiduels évidents.

Des paramètres procéduraux clés, incluant le temps d’ablation, la plage de puissance, l’estimation de l’énergie livrée, l’utilisation de l’hydrodissection et tout événement intraprocédural, ont été enregistrés. Après l’ablation, le patient a été observé selon le protocole institutionnel, et des instructions post-procédurales standardisées ont été fournies. L’ablation était temporairement interrompue ou interrompue lorsque les patients présentaient une douleur importante, un changement de voix, une toux incontrôlée, un gonflement progressif ou d’autres constatations suscitant des inquiétudes en cas de lésion thermique ou de complications procédurales.

Imagerie et documentation immédiatement post-ablation
Dans les 30 minutes suivant la fin de la RFA, un examen échographique post-procédural immédiat a été effectué en utilisant les mêmes préréglages d’imagerie que la ligne de départ. Des images en niveaux de gris ont été acquises pour documenter l’étendue de la zone traitée et identifier les caractéristiques typiques post-ablation, y compris l’hyperéchogénicité transitoire et les changements d’architecture interne. La vascularisation intranodulaire et périnodule a été évaluée à l’aide de Doppler coloré ou de puissance afin de confirmer la réduction ou l’absence du flux sanguin dans la zone ablée.

L’acquisition de l’élastographie a été répétée après la résolution du mouvement transitoire lié à la procédure et la capacité du patient à maintenir une suspension respiratoire stable, suivant le même protocole d’acquisition qu’au départ. Les valeurs de rigidité étaient enregistrées, et des cartes d’élastographie représentatives ainsi que des images en mode B correspondantes étaient stockées. Toute limitation technique rencontrée lors de l’acquisition, comme l’inconfort du patient ou les artefacts de mouvement, a été documentée. Ces images immédiatement post-ablation servaient de référence procédurale et n’étaient pas utilisées comme point d’arrêt principal, mais fournissaient un contexte pour les modifications longitudinales ultérieures.

Calendrier de suivi et points de terminaison standardisés
Les visites de suivi étaient programmées à 1 mois, 3 mois, 6 mois et 12 mois après l’ablation. À chaque visite, une évaluation standardisée a été réalisée, incluant l’évaluation des symptômes, la classification esthétique lorsque applicable, l’échographie conventionnelle et l’élastographie en utilisant les mêmes paramètres d’acquisition et la même stratégie de mesure qu’au départ. Le principal critère de l’élastographie était le changement de rigidité du nodule indice par rapport à la ligne de base à des moments de suivi prédéfinis, avec un accent particulier sur un changement précoce à un mois. Le principal point d’épreuve d’imagerie conventionnel était le VRR à 12 mois. Les critères secondaires comprenaient la VRR lors des visites de suivi antérieures, les trajectoires longitudinales de raideur, et la présence de tissu vascularisé résiduel ou récurrent dans la zone traitée. Les événements indésirables et toute intervention supplémentaire lors du suivi ont été consignés de manière systématique.

Gestion des données et aveuglement
Chaque participant s’est vu attribuer un identifiant d’étude unique afin d’assurer l’anonymisation. Les données cliniques et les fichiers d’imagerie étaient stockés dans une base de données institutionnelle sécurisée, l’accès étant réservé au personnel autorisé de l’étude. L’élastographie et les mesures volumétriques sont réalisées par des lecteurs formés, aveugles aux scores cliniques des symptômes et des résultats antérieurs à l’élastographie, lorsque cela est possible. Pour garantir la cohérence, tenez un manuel de mesure écrit spécifiant les règles de placement par région d’intérêt, les critères de qualité d’image et les procédures de calcul. Dans les cas d’écart important entre lecteurs dépassant un seuil prédéfini, les mesures sont examinées conjointement et résolues par consensus. Toutes les modifications apportées à l’ensemble de données sont consignées afin de préserver une trace d’audit. Toutes les acquisitions d’échographie et d’élastographie ont été réalisées par des opérateurs ayant une expérience préalable en échographie thyroïdienne et en évaluation par imagerie RFA, qui ont suivi une formation spécifique au protocole avant le début de l’étude. Les images stockées étaient examinées hors ligne en utilisant les mêmes règles prédéfinies de placement par région d’intérêt et les mêmes critères de contrôle qualité appliqués lors de l’acquisition.

Analyse statistique
Les variables continues ont été évaluées pour les caractéristiques de distribution à l’aide d’une inspection visuelle des histogrammes et du test de Shapiro–Wilk. Les variables normalement distribuées sont présentées en moyenne ± écart-type, tandis que les variables non normalement distribuées sont résumées en médiane avec une plage interquartile. Les variables catégorielles sont rapportées sous forme de comptages et de pourcentages. L’objectif analytique principal était d’évaluer les modifications intra-nodules de la rigidité dérivée de l’élastographie après l’RFA. Pour la comparaison primaire par paires entre les points de référence et post-ablation, les différences de valeurs de rigidité ont été analysées à l’aide de tests t appariés lorsque la distribution des différences appariées approchait la normalité ; sinon, le test du grade signé Wilcoxon était appliqué. Les tailles d’effet sont rapportées comme des variations moyennes ou médianes avec des intervalles de confiance correspondants à 95 % lorsque cela est approprié. Pour caractériser les variations longitudinales des mesures d’élastographie sur toutes les visites de suivi (référence, immédiatement après l’ablation, 1 mois, 3 mois, 6 mois et 12 mois), des modèles linéaires à effets mixtes ont été construits avec des interceptions aléatoires au niveau du participant pour tenir compte de la corrélation interne-sujet. Le temps a été modélisé comme un effet fixe catégorique plutôt que comme une variable continue, car les changements biomécaniques post-ablation n’étaient pas censés suivre une trajectoire linéaire à travers les visites de suivi. En particulier, le plan de l’étude anticipait une augmentation aiguë immédiate de la rigidité post-ablation, suivie d’un remodelage progressif et d’un déclin tardif, rendant la modélisation catégorique plus appropriée pour capturer des schémas de réponse temporelle potentiellement non linéaires sans imposer d’hypothèses de linéarité. Les hypothèses du modèle ont été évaluées par inspection des parcelles résiduelles. Lorsque cela est approprié, les valeurs de rigidité ont été transformées en logarthmique pour améliorer l’ajustement du modèle. Les associations entre les variations de raideur dérivées de l’élastographie et les métriques de réponse conventionnelles au traitement ont été explorées dans des analyses secondaires. La corrélation entre le changement précoce de rigidité et la VRR lors du suivi ultérieur a été évaluée à l’aide des coefficients de corrélation de Pearson ou de Spearman, selon la distribution des données. Des modèles de régression multivariée ont été utilisés pour examiner si les changements précoces de l’élastographie prédisaient indépendamment la réponse volumétrique ultérieure après ajustement du volume nodulique de base, de la composition des nodules (solide versus principalement solide), de l’énergie délivrée et de l’âge. Les termes d’interaction entre variables de sous-groupes et temps ont été testés dans le cadre de modélisation à effets mixtes. Pour les analyses longitudinales, les modèles à effets mixtes prennent intrinsèquement en compte les observations manquantes sous l’hypothèse qu’elles manquent au hasard. Tous les tests statistiques étaient bidirectionnels, et une valeur P < 0,05 a été considérée comme indiquant une signification statistique.

Résultats

Conception de l’étude et schéma général
Les participants exclus ne remplissaient pas les critères d’éligibilité prédéfinis, disposaient de données d’imagerie de base incomplètes ou manquaient d’un suivi longitudinal standardisé. Parmi les 187 participants inclus, une classification complète des résultats sur 12 mois était disponible pour tous les cas inclus dans l’analyse primaire, bien que des observations manquantes intermittentes isolées lors des suivis individuels aient été prises en compte dans le cadre de modélisation à effets mixtes (Figure 1). En utilisant la définition prédéfinie de la réponse de 12 mois basée sur le VRR, 146 nodules ont été classés comme répondants (VRR ≥ 50 %) et 41 comme non-répondants (VRR < 50 %) (Figure 1). Cette conception a permis des comparaisons appariées intra-nodule (la référence étant le contrôle interne) et une évaluation longitudinale standardisée sur des points temporels de suivi prédéfinis.

Caractéristiques de base du participant et du nodule
Les caractéristiques cliniques de base sont résumées dans le tableau 1. Les non-répondants étaient plus âgés que les répondants (49,63 ± 8,01 contre 45,94 ± 7,88 ans, P = 0,009) et présentaient un taux de base de symptômes et une charge esthétique légèrement plus élevés (score de symptômes de pression 3,43 ± 0,57 contre 3,19 ± 0,63, P = 0,029 ; score de préoccupation esthétique 3,01 ± 0,47 contre 2,79 ± 0,49, P = 0,011). La fonction thyroïdienne était largement similaire entre les groupes, bien que la TSH ait été légèrement plus élevée chez les non-répondants (1,86 ± 0,31 mUI/L contre 1,74 ± 0,33 mIU/L, P = 0,039) (Tableau 1). Ces différences de base indiquent un profil clinique moins favorable chez les non-répondants avant l’intervention. L’échographie de base a montré des nodules significativement plus importants chez les non-répondeurs, notamment un diamètre maximal plus important (33,87 ± 4,03 mm contre 29,31 ± 3,94 mm, P < 0,001) et un volume de base (12,04 ± 2,41 mL contre 8,33 ± 1,90 mL, P < 0,001) (Tableau 1). La composition des nodules différait également, avec un déplacement vers des catégories de composantes kystiques plus élevées chez les non-répondants (P < 0,001) (Tableau 1). Ces caractéristiques échographies de base fournissent un contexte anatomique pour les schémas de réponse volumétriques et biomécaniques ultérieurs. Aucun participant n’a montré d’anomalies parenchymeuses thyroïdiennes diffuses suffisamment graves pour répondre aux critères d’exclusion prédéfinis pour la maladie thyroïdienne diffuse.

Les mesures d’élastographie ont été obtenues dans des conditions standardisées avec des règles de région d’intérêt et des critères d’acceptation prédéfinis afin de minimiser les artefacts de mouvement et de compression. L’utilisation de mesures de base dans les nodules internes comme références internes peut aider à réduire l’influence de la variabilité interindividuelle et soutenir l’interprétation des variations longitudinales de rigidité. Les caractéristiques procédurales sont indiquées dans le Tableau 1. Conformément à une charge nodulique de base plus élevée, les non-répondants ont eu une durée de procédure plus longue (21,13 ± 3,01 min contre 18,42 ± 2,79 min, P < 0,001) et une énergie délivrée plus élevée (13,18 ± 2,63 kJ contre 11,04 ± 1,98 kJ, P < 0,001) (Tableau 1). Ces différences procédurales étaient attendues, compte tenu des différences de taille et de composition de base, et ont été prises en compte dans des analyses ajustées de prédiction de réponse.

Trajectoires d’élastographie longitudinale après RFA
Immédiatement après l’RFA, les deux modalités d’élastographie ont montré une augmentation marquée de la raideur aiguë (Tableau 2 ; Figure 2), cohérente avec l’effet biomécanique immédiat supposé de la coagulation thermique. L’élastographie par ondes de cisaillement (SWE) est passée de 31,5 ± 5,1 à 41,7 ± 5,6 kPa chez les répondants et de 35,1 ± 5,3 à 43,4 ± 5,8 kPa chez les non-répondants (Tableau 2 ; Figure 2). Le rapport de déformation a augmenté de manière concordante (répondants 2,41 ± 0,37 à 3,14 ± 0,41 ; non-répondants 2,77 ± 0,40 à 3,35 ± 0,42) (Tableau 2 ; Figure 2). En revanche, le volume immédiatement post-ablation est resté proche de la référence dans les deux groupes (Tableau 2), soutenant l’hypothèse de l’étude selon laquelle l’élastographie fournit un indicateur d’imagerie précoce et précoce du changement tissulaire lié au traitement par rapport à la volumetrie à court intervalle.

Résultats volumétriques longitudinaux et VRR
Au cours du suivi, l’élastographie a montré un remodelage progressif après la montée immédiatement après l’ablation (Tableau 2 ; Figure 2). Chez les répondants, l’EES est tombée en dessous du seuil initial de 3 mois et a continué de diminuer pendant 12 mois (31,5 ± 5,1 kPa au départ à 21,6 ± 4,7 kPa à 12 mois), tandis que les non-répondants ont montré une réduction atténuée et sont restés nettement plus rigides à 12 mois (31,2 ± 5,2 kPa) (Tableau 2 ; Figure 2). Le ratio de déformation a suivi le même schéma, diminuant à 1,87 ± 0,35 à 12 mois chez les répondants contre 2,55 ± 0,41 chez les non-répondants (Tableau 2 ; Figure 2).

Les résultats volumétriques ont évolué en parallèle avec les changements biomécaniques (Tableau 2 ; Figure 3). Le volume des nodules diminuait régulièrement dans les deux groupes, avec une séparation persistante entre les points temporels ; à 12 mois, les répondants atteignaient 2,20 ± 0,45 mL contre 7,35 ± 1,80 mL chez les non-répondants (Tableau 2 ; Figure 3). Le VRR s’est également séparé tôt et est resté distinct : à 1 mois, le VRR était de 31,5 ± 6,3 % chez les répondants contre 16,1 ± 3,8 % chez les non-répondants, et à 12 mois, le VRR était de 73,6 ± 6,9 % contre 38,5 ± 8,7 %, respectivement (Tableau 2 ; Figure 3). Ensemble, ces trajectoires alignées d’élastographie et de VRR soutiennent l’hypothèse que la remodelisation biomécanique post-RFA suit et différencie la réponse volumétrique durable.

Changement précoce de l’élastographie comme prédicteur indépendant de la réponse volumétrique à 12 mois
Les modèles à effets mixtes corroboraient les trajectoires descriptives et quantifiaient les différences de groupe au fil du temps (Tableau 3). Pour le SWE, les effets temporels ont indiqué une augmentation immédiate significative par rapport à la ligne de base (estimation +10,22 kPa, P < 0,001), suivie de valeurs progressivement plus faibles lors de visites ultérieures (par exemple, 12 mois contre estimation de base -9,02 kPa, P < 0,001). Il est important de noter que les interactions temps par groupe ont montré une divergence lors d’un suivi ultérieur, les non-répondants présentant une SE plus élevée par rapport aux répondants à partir de 3 mois (non-répondant × estimation à 3 mois +4,20 kPa, P < 0,001 ; non-répondant × estimation à 12 mois +9,41 kPa, P < 0,001) (Tableau 3). Les modèles de rapport de déformation ont montré le même schéma, avec une divergence tardive significative (estimation des non-répondants × 12 mois +0,68, P < 0,001) (Tableau 3).

Dans les analyses de corrélation, les changements précoces de l’élastographie post-ablation (immédiat moins la ligne de départ) étaient positivement associés à une réponse volumétrique à 12 mois (Figure 4). Le rapport de déformation Δ a montré une corrélation modérée avec la VRR à 12 mois (Pearson r = 0,445, P < 0,001), tandis que ΔSWE a montré une corrélation plus faible mais statistiquement significative ( Pearson r = 0,302, P < 0,001). Ces résultats indiquent que des changements aigus plus importants de raideur post-ablation étaient associés à une réduction accrue du volume ultérieur, soutenant l’hypothèse que l’élastographie capture des effets biomécaniques précoces qui sont informatifs pour une réponse à plus long terme au traitement.

Dans la régression multivariable en ajustant selon le volume nodulique de base, la composition, l’énergie délivrée et l’âge, les changements précoces de l’élastographie post-ablation prédisaient indépendamment un VRR à 12 mois (Tableau 4). Le changement précoce du rapport de déformation (immédiat moins la ligne de base) était fortement associé à un VRR à 12 mois plus élevé (β = 14,80 % par +1,0, P < 0,001), et le changement précoce de l’EAT a également contribué indépendamment (β = 1,80 % par +5 kPa, P = 0,011). Un volume de référence plus élevé était associé à un VRR plus faible sur 12 mois (β = −1,05 % par +1 mL, P < 0,001), tandis que la composition solide et l’énergie délivrée plus élevée étaient associées à un VRR plus élevé ( P < 0,05) (Tableau 4).

Lors du suivi, aucun nodule traité n’a été reclassé comme pathologie maligne ou alternative sur la base d’imagerie répétée, de cytologie ou d’évaluation clinique. Bien que certains nodules non réactifs aient montré un volume résiduel persistant ou une réduction atténuée du volume, une repousse manifeste nécessitant une ablation ou une chirurgie répétée n’a pas été observée durant la période de suivi de 12 mois.

DISPONIBILITÉ DES DONNÉES :
Les données brutes sous-jacentes aux analyses présentées dans cette étude ont été préparées et téléchargées sous forme de Tableau Supplémentaire 1 afin d’améliorer la transparence et la reproductibilité. L’ensemble de données comprend les variables démographiques, cliniques, échographies, élastographiques, procédurales et de suivi longitudinales au niveau des participants non identifiées, utilisées dans les analyses.

Figure 1 : Conception de l’étude et flux de travail analytique. Les participants dépistés pendant la période d’étude ont été évalués pour leur éligibilité et inclus après des exclusions prédéfinies. Un nodule thyroïdien indice par participant a été analysé. Le statut de réponse à 12 mois était défini par le rapport de réduction de volume (VRR) : répondant (VRR ≥ 50 %) et non-répondant (VRR < 50 %). Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 2 : Rigidité de l’élastographie longitudinale à ondes de cisaillement après ablation par radiofréquence. La raideur moyenne de l’élastographie par onde de cisaillement (SWE) est indiquée au départ, immédiatement après l’ablation, et à 1, 3, 6 et 12 mois pour la cohorte globale, et stratifiée selon le niveau de réponse à 12 mois. Les barres d’erreur indiquent l’écart-type. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 3 : Variations longitudinales du volume des nodules et du rapport de réduction du volume. (A) Volume moyen des nodules (mL) et (B) rapport moyen de réduction du volume (VRR, %) au départ, immédiatement après l’ablation, et à 1, 3, 6 et 12 mois, stratifié selon le statut de réponse à 12 mois. Les barres d’erreur indiquent l’écart-type. Les crochets indiquent des comparaisons entre groupes à chaque moment. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 4 : Relation entre les changements précoces de l’élastographie et la réponse volumétrique à 12 mois. Les diagrammes de points montrent des associations entre les variations précoces de l’élastographie post-ablation (immédiat moins la ligne de base) et la VRR à 12 mois : (A) rapport de déformation et (B) raideur ΔWE. Les positions des points de dispersion ont été conservées de l’analyse globale initiale pour faciliter la comparaison visuelle entre les groupes. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Section	Variable	Répondant (n = 146)	Non-répondant (n = 41)	Statistique de test (t/χ2)	Valeur p
Démographie et état clinique de base
	Âge, années	45,94 ± 7,88	49,63 ± 8,01	-2.640	0.009
	Sexe	Femmes : 119 (81,5 %) ; Hommes : 27 ans (18,5 %)	Femmes : 29 ans (70,7 %) ; Hommes : 12 (29,3 %)	2.252	0.133
	IMC, kg/m 2	23,74 ± 1,96	24,74 ± 1,94	-1.914	0.057
	Fumer (actuel/jamais jamais)	N° : 126 (86,3 %) ; Oui : 20 (13,7 %)	No : 32 (78,0 %) ; Oui : 9 (22,0 %)	1.664	0.197
	Hypertension	N° : 119 (81,5 %) ; Oui : 27 ans (18,5 %)	No : 30 (73,2 %) ; Oui : 11 (26,8 %)	1.374	0.241
	Diabète sucré	N° : 136 (93,2 %) ; Oui : 10 (6,8 %)	N° : 37 (90,2 %) ; Oui : 4 (9,8 %)	0.391	0.532
	Maladie coronarienne	N° : 140 (95,9 %) ; Oui : 6 (4,1 %)	N° : 38 (92,7 %) ; Oui : 3 (7,3 %)	0.719	0.397
	Maladie rénale chronique	N° : 143 (97,9 %) ; Oui : 3 (2,1 %)	No : 40 (97,6 %) ; Oui : 1 (2,4 %)	0.023	0.881
	BPCO/asthme	N° : 140 (95,9 %) ; Oui : 6 (4,1 %)	N° : 39 (95,1 %) ; Oui : 2 (4,9 %)	0.046	0.830
	Score des symptômes de pression, 0–10	3.19 ± 0.63	3,43 ± 0,57	-2.199	0.029
	Score de préoccupation cosmétique, 1–4	2,79 ± 0,49	3.01 ± 0.47	-2.562	0.011
Analyses thyroïdiennes et auto-immunité
	TSH, mIU/L	1,74 ± 0,33	1,86 ± 0,31	-2.084	0.039
	T4 libre, pmol/L	15.06 ± 2.03	14,95 ± 1,92	0.310	0.757
	TPOAb positif	N° : 129 (88,4 %) ; Oui : 17 ans (11,6 %)	N° : 35 (85,4 %) ; Oui : 6 (14,6 %)	0.265	0.606
	TgAb positif	N° : 132 (90,4 %) ; Oui : 14 (9,6 %)	N° : 36 (87,8 %) ; Oui : 5 (12,2 %)	0.238	0.626
	Utilisation de la lévothyroxine	N° : 134 (91,8 %) ; Oui : 12 (8,2 %)	N° : 36 (87,8 %) ; Oui : 5 (12,2 %)	0.612	0.434
	Utilisation des antiplaquettes/anticoagulants	N° : 139 (95,2 %) ; Oui : 7 (4,8 %)	N° : 38 (92,7 %) ; Oui : 3 (7,3 %)	0.402	0.526
Caractéristiques d’échographie de base
	Emplacement du nodule	Gauche : 66 (45,2 %) ; Isthme : 6 (4,1 %) ; À droite : 74 (50,7 %)	Gauche : 19 (46,3 %) ; Isthme : 2 (4,9 %) ; À droite : 20 (48,8 %)	0.076	0.963
	Diamètre maximal, mm	29.31 ± 3.94	33,87 ± 4,03	-6.516	<0,001
	Volume de nodule de base, mL	8,33 ± 1,90	12.04 ± 2.41	-10.385	<0,001
	Composition (solide)	N° : 44 (30,1 %) ; Oui : 102 (69,9 %)	N° : 17 (41,5 %) ; Oui : 24 ans (58,5 %)	1.868	0.172
	Composante kystique	<10 % : 89 (61,0 %) ; 10–25%: 47 (32.2%); >25 % : 10 (6,8 %)	<10 % : 10 (24,4 %) ; 10–25%: 20 (48.8%); >25 % : 11 (26,8 %)	21.923	<0,001
	Calcification	Grossier : 10 (6,8 %) ; Micro : 35 (24,0 %) ; Aucun : 101 (69,2 %)	Grossières : 3 (7,3 %) ; Micro : 8 (19,5 %) ; Aucun : 30 (73,2 %)	0.360	0.835
	Grade de vascularisation de base	0: 27 (18.5%); 1: 47 (32.2%); 2: 53 (36.3%); 3: 19 (13.0%)	0: 5 (12.2%); 1: 13 (31.7%); 2: 15 (36.6%); 3: 8 (19.5%)	1.681	0.641
Élastographie de base
	Rigidité de base de l’équipement logiciel, kPa	31,50 ± 5,11	35.08 ± 5.34	-3.925	<0,001
	Rapport de déformation de base	2,41 ± 0,37	2,77 ± 0,40	-5.407	<0,001
Paramètres procéduraux et sécurité
	Temps de la procédure, minimum	18.42 ± 2.79	21.13 ± 3.01	-5.401	<0,001
	Puissance d’ablation, W	35.02 ± 4.78	37.05 ± 4.63	-2.419	0.017
	Énergie délivrée, kJ	11.04 ± 1.98	13.18 ± 2.63	-5.665	<0,001
	Hydrodissection utilisée	No : 100 (68,5 %) ; Oui : 46 (31,5 %)	N° : 25 (61,0 %) ; Oui : 16 (39,0 %)	0.816	0.366
	Sédation utilisée	N° : 115 (78,8 %) ; Oui : 31 (21,2 %)	No : 30 (73,2 %) ; Oui : 11 (26,8 %)	0.576	0.448
	Douleur intra-procédurale, AV 0–10	2,60 ± 0,59	2,78 ± 0,61	-1.713	0.088
	Saignement/hématome mineur	N° : 139 (95,2 %) ; Oui : 7 (4,8 %)	N° : 37 (90,2 %) ; Oui : 4 (9,8 %)	1.423	0.233
	Changement de voix transitoire	N° : 143 (97,9 %) ; Oui : 3 (2,1 %)	N° : 41 (100,0 %) ; Oui : 0 (0,0 %)	0.856	0.355

Tableau 1 : Caractéristiques démographiques de base, cliniques, de laboratoire, d’échographie, d’élastographie et de procédures stratifiées par réponse volumétrique à 12 mois après ablation radiofréquence (VRR ≥50 % vs <50 %). Les données sont indiquées en moyenne ± DE (ou médiane [IQR]) et n ( %). Les comparaisons entre groupes utilisaient le test t de Student ou le test de somme de rang de Wilcoxon pour les variables continues et χ² ou le test exact de Fisher pour les variables catégorielles, selon le cas. Un P < à deux faces de 0,05 indiquait une signification statistique. Abréviations : IMC, indice de masse corporelle ; TSH, hormone thyroïdienne ; FT4, thyroxine libre ; TPOAb, anticorps de peroxydase thyroïdienne ; TgAb, anticorps thyroglobuliniques ; SWE, élastographie par ondes de cisaillement ; ROI, région d’intérêt ; VAS, échelle analogique visuelle ; RFA, ablation par radiofréquence ; VRR, ratio de réduction de volume.

Visite	Groupe	n	Volume (mL), moyenne ± SD	VRR ( %), moyenne ± SD	SWE (kPa), moyenne ± SD	Rapport de déformation, moyenne ± SD
Référence	Répondeur	146	8,33 ± 1,90	—	31,5 ± 5,1	2,41 ± 0,37
	Non-répondant	41	12.04 ± 2.41	—	35.1 ± 5.3	2,77 ± 0,40
Immédiat	Répondeur	146	8,75 ± 1,99	—	41,7 ± 5,6	3.14 ± 0.41
	Non-répondant	41	12,66 ± 2,54	—	43,4 ± 5,8	3,35 ± 0,42
1 mois	Répondeur	146	5,69 ± 1,35	31,5 ± 6,3	33,6 ± 5,0	2,51 ± 0,39
	Non-répondant	41	10.08 ± 2.36	16.1 ± 3.8	36,1 ± 5,3	2,82 ± 0,41
3 mois	Répondeur	146	3,72 ± 0,85	55.3 ± 7.1	26,4 ± 4,8	2,05 ± 0,36
	Non-répondant	41	8,70 ± 2,10	28,1 ± 6,8	32,6 ± 5,1	2,62 ± 0,39
6 mois	Répondeur	146	2,99 ± 0,70	64.1 ± 7.0	23,8 ± 4,7	1,91 ± 0,35
	Non-répondant	41	7,93 ± 1,95	34.2 ± 8.1	31,6 ± 5,3	2,56 ± 0,40
12 mois	Répondeur	146	2,20 ± 0,45	73,6 ± 6,9	21,6 ± 4,7	1,87 ± 0,35
	Non-répondant	41	7,35 ± 1,80	38,5 ± 8,7	31.2 ± 5.2	2,55 ± 0,41

Tableau 2 : Variations longitudinales du volume nodulaire, du ratio de réduction de volume et des métriques d’élastographie de la base à 12 mois suivant une ablation par radiofréquence, stratifiées par réponse au traitement. Les valeurs sont résumées à chaque visite ; n désigne les observations disponibles. Le VRR a été calculé par rapport au volume de référence. L’inférence longitudinale a été réalisée à l’aide de modèles à effets mixtes. Abréviations : VRR, ratio de réduction de volume ; SWE, élastographie par ondes de cisaillement ; RFA, ablation par radiofréquence.

Conséquences	Terme à effet fixe	Estimation	SE	t	Valeur p
Rigidité SWE (kPa)
	(Interception) Référence, répondeur	31.51	0.55	57.273	<0,001
	Groupe : Non-répondant (vs répondant)	3.58	0.92	3.891	<0,001
	Heure : Immédiate (vs Baseline)	10.22	0.62	16.452	<0,001
	Durée : 1 mois (vs Baseline)	2.17	0.65	3.502	<0,001
	Durée : 3 mois (par rapport à la base)	-5.14	0.61	-8.361	<0,001
	Durée : 6 mois (contre la base)	-7.73	0.64	-12.419	<0,001
	Durée : 12 mois (par rapport à la base)	-9.02	0.63	-14.286	<0,001
	Temps × groupe : Non-répondant × immédiat	-0.80	0.98	-0.816	0.415
	Période × groupe : Non-répondant × 1 mois	0.45	0.95	0.474	0.636
	Temps × groupe : Non-répondant × 3 mois	4.20	0.96	4.375	<0,001
	Groupe × Temps : Non-répondant × 6 mois	6.10	0.98	6.224	<0,001
	Groupe × Temps : Non-répondant × 12 mois	9.41	1.02	9.216	<0,001
Rapport de déformation
	(Interception) Référence, répondeur	2.41	0.04	60.252	<0,001
	Groupe : Non-répondant (vs répondant)	0.36	0.07	5.143	<0,001
	Heure : Immédiate (vs Baseline)	0.73	0.05	14.601	<0,001
	Durée : 1 mois (vs Baseline)	0.10	0.07	2.012	0.046
	Durée : 3 mois (par rapport à la base)	-0.36	0.08	-7.211	<0,001
	Durée : 6 mois (contre la base)	-0.51	0.06	-10.032	<0,001
	Durée : 12 mois (par rapport à la base)	-0.54	0.05	-10.802	<0,001
	Temps × groupe : Non-répondant × immédiat	-0.08	0.08	-1.000	0.318
	Période × groupe : Non-répondant × 1 mois	0.05	0.09	0.625	0.532
	Temps × groupe : Non-répondant × 3 mois	0.32	0.04	4.031	<0,001
	Groupe × Temps : Non-répondant × 6 mois	0.55	0.08	6.875	<0,001
	Groupe × Temps : Non-répondant × 12 mois	0.68	0.07	8.521	<0,001

Tableau 3 : Résultats du modèle linéaire à effets mixtes pour les trajectoires d’élastographie longitudinale sur les visites de suivi, incluant les effets fixes pour l’interaction temps et temps par réponse avec les interceptions aléatoires au niveau du participant. Les modèles linéaires à effets mixtes incluaient des interceptions aléatoires au niveau du participant et des effets en temps catégoriel ; Les termes de groupe de réponse et d’interaction temps par groupe ont été utilisés pour comparer les trajectoires. L’ajustement du modèle a été évalué à l’aide de diagnostics résiduels ; la transformation logarithmique était appliquée lorsque nécessaire. Les estimations sont rapportées avec des IC de 95 %. Abréviations : SWE, élastographie à ondes de cisaillement ; IC, intervalle de confiance ; VRR, ratio de réduction de volume.

Conséquences	Mandat	Beta	SE	t	df	Valeur p	IC à 95 %
VRR à 12 mois ( %)
	(Interception)	82.00	4.20	19.524	180	<0,001	73,71 à 90,29
	Δ Rapport de déformation (immédiat - Base), par 1,0	14.80	2.40	6.167	180	<0,001	10h06 à 19h54
	Δ SWE (immédiat - Référence), par 5 kPa	1.80	0.70	2.571	180	0.011	0,42 à 3,18
	Volume de nodule de base, par 1 mL	-1.05	0.28	-3.752	180	<0,001	-1,60 à -0,50
	Composition solide (Oui contre Non)	3.20	1.40	2.286	180	0.023	0,44 à 5,96
	Énergie livrée, par 1 kJ	0.55	0.22	2.512	180	0.013	0,12 à 0,98
	Âge, par 10 ans	-1.10	0.70	-1.571	180	0.118	-2,48 à 0,28

Tableau 4 : Analyses de régression multivariées évaluant si les changements précoces de l’élastographie post-ablation prédisent indépendamment la réponse volumétrique à 12 mois après ajustement pour les nodules de base et les covariables procédurales. La régression linéaire multivariée a évalué les associations entre le changement précoce de l’élastographie (immédiat moins la ligne de base) et la VRR à 12 mois après ajustement des covariables. Des coefficients (β) avec SE, t, df, P bilatéral et IC à 95 % sont rapportés ; df représente des degrés de liberté résiduels. Abréviations : VRR, ratio de réduction de volume ; SWE, élastographie par ondes de cisaillement ; SE, erreur standard ; IC, intervalle de confiance ; DF, degrés de liberté.

Fichier supplémentaire 1 : Données brutes de cette étude.Veuillez cliquer ici pour télécharger ce fichier.

Discussion

Cette étude a évalué une approche basée sur l’élastographie pour évaluer la réponse au traitement après RFA des nodules thyroïdiens en utilisant l’élastographie comme biomarqueur d’imagerie quantitative associé à la modification de la raideur post-ablation. La proposition centrale est que la rigidité, mesurée par l’élastographie logicielle et la déformation, peut servir d’indicateur précoce de l’effet d’ablation et de marqueur longitudinal du remodelage post-ablation, complétant les extrémités conventionnelles basées sur le volume. Les résultats soutiennent ce concept à deux niveaux. Premièrement, la raideur et le rapport de déformation de l’éticulation SWE ont augmenté immédiatement après l’ablation, un schéma cohérent avec une coagulation thermique aiguë et une réaction tissulaire précoce. Deuxièmement, la raideur a ensuite diminué au fil du suivi, les répondants montrant une réduction tardive plus importante que les non-répondeurs, ce qui reflète la séparation observée dans les résultats volumétriques et la VRR. De plus, les changements précoces post-ablation dans l’élastographie ont montré des associations statistiquement significatives avec la VRR à 12 mois et sont restés associés indépendamment après ajustement pour les nodules de base et les covariables procédurales. Ensemble, ces résultats soutiennent l’hypothèse selon laquelle l’élastographie capture les variations longitudinales de rigidité associées à la réponse au traitement par RFA, et que les premiers déplacements de raideur portent des informations pronostiques qui ne sont pas entièrement expliquées uniquement par la taille de base, la composition ou la délivrance^d’énergie 12,13,14,15,16.

Les trajectoires observées sont cohérentes avec les changements tissulaires post-ablation décrits précédemment et offrent un récit cliniquement utile pour l’imagerie^{post-ablation} 12,13,14. L’augmentation immédiate de la raideur reflète probablement une nécrose coagulative dans la zone d’ablation, ainsi que l’œdème et un changement inflammatoire précoce, qui augmentent la résistance à la déformation et modifient la propagation^{des ondes 15,16}. Au fil du temps, la baisse de la raideur est cohérente avec la résorption progressive du tissu nécrotique, l’organisation de la zone d’ablation et le remodelage qui aboutissent finalement à un rétrécissement et à une amélioration des^symptômes 15,17,18,19. Le fait que les intervenants aient montré une réduction tardive plus importante de la raideur suggère qu’une ablation réussie n’est pas simplement un événement aigu mais un processus de remodelage soutenu, et que le traitement incomplet ou le tissu résiduel viable peuvent s’exprimer par une trajectoire d’adoucissement^{émoussée 20}. Cette perspective aide à concilier une limitation courante de l’évaluation précoce du volume. Immédiatement après l’RFA, le volume peut rester proche de la ligne de base ou fluctuer transitoirement en raison d’œdème, d’hémorragie ou de changement réactif, ce qui fait de la volumetrie précoce un indicateur imparfait de la suffisance de la zone d’ablation. L’élastographie, en revanche, répond à la principale conséquence physique de l’ablation à un moment où les métriques basées sur la taille sont les moins^{discriminantes 21}. La divergence dans les trajectoires tardives de rigidité, quantifiée dans les modèles à effets mixtes, renforce le concept que la raideur capture des différences significatives dans l’évolution tissulaire qui correspondent à une réduction durable du volume.

Du point de vue des avancées scientifiques, ce travail contribue en faisant passer l’élastographie d’un rôle principalement diagnostique à un paradigme de surveillance thérapeutique. Dans les soins bénins des nodules thyroïdiens, le domaine adopte de plus en plus l’ablation guidée par imagerie, mais les stratégies de suivi dépendent encore fortement de résultats volumétriques retardés tels que le VRR entre 6 et 12 mois. Un biomarqueur fournissant une lecture précoce objective de l’effet du traitement pourrait fournir des informations d’imagerie supplémentaires durant la période post-procédurale précoce et soutenir des recherches supplémentaires sur les stratégies de suivi adaptées au^{risque 22,23}. Des approches similaires basées sur l’élastographie ont également été explorées dans d’autres contextes d’ablation^{thermique 24}. La cartographie de la rigidité est attrayante car elle n’est pas invasive, peut être intégrée dans les flux de travail échographiques de routine et peut être répétée sans agents de contraste ni^{rayonnement 25}. D’un point de vue de flux de travail, l’élastographie peut être intégrée dans les examens échographiques post-ablation de routine sans prolongation substantielle du temps d’examen lorsque des protocoles d’acquisition standardisés sont utilisés. Cependant, une mise en œuvre plus large nécessitera une attention particulière à la variabilité multiplateforme, à l’expérience de l’opérateur et à l’harmonisation des paramètres d’acquisition entre différents systèmes d’échographie. Pour réduire cette source de variabilité dans la présente étude, toutes les mesures d’élastographie ont été prises à l’aide de la même plateforme d’échographie et du même environnement logiciel tout au long du suivi. L’utilisabilité clinique peut dépendre moins des valeurs de rigidité absolues que de schémas de variation longitudinale cohérents obtenus dans des conditions d’acquisition stables. En conséquence, la formation standardisée des opérateurs, les règles de contrôle qualité prédéfinies et l’harmonisation des protocoles restent essentielles avant que l’élastographie puisse être intégrée plus largement dans les voies de surveillance post-ablation de routine. Les résultats présents suggèrent qu’un protocole standardisé d’acquisition de l’élastographie pourrait permettre d’obtenir des schémas de rigidité longitudinale plus cohérents, adaptés aux applications de recherche clinique. Il est important de noter que l’approche actuelle considère les valeurs de base des nodules internes comme des contrôles internes, mettant ainsi l’accent sur des indicateurs de changement qui peuvent être plus robustes chez les individus que la rigidité absolue seule. Un autre point à considérer est que le protocole actuel d’élastographie se concentrait principalement sur les mesures de raideur au sein de la composante solide du nodule traité, tandis que le tissu périnodulaire n’était pas évalué de manière systématique. Cela représente une limitation potentiellement importante car la zone de traitement périphérique peut contenir des informations biologiquement significatives liées à la réponse inflammatoire, à la diffusion thermique, à l’évolution de la fibrose, à la présence de tissu résiduel viable et au remodelage précoce des tissus après ablation. La littérature émergente en imagerie met de plus en plus l’accent sur la pertinence diagnostique et biologique des caractéristiques périphériques et microenvironnementales des lésions dans l’évaluation traitement-réponse. En conséquence, de futures études intégrant la cartographie combinée de la rigidité intranodulaire et périnodule pourraient fournir une caractérisation biomécanique plus complète de l’évolution tissulaire post-ablation et améliorer encore l’évaluation précoce de la réponse après l’RFA thyroïdienne.

Dans la mise en œuvre pratique, plusieurs pièges techniques doivent être pris en compte lors de l’acquisition de l’élastographie et de l’évaluation post-ablation. Une compression excessive du transducteur, un mouvement d’englutition, une suspension respiratoire instable, une dégénérescence kystique hétérogène et un remplissage de couleur incomplet peuvent tous introduire une instabilité de mesure et réduire la reproductibilité. L’évaluation immédiatement post-ablation peut également être influencée par un œdème, une hémorragie transitoire ou une modification inflammatoire réactive. Maintenir un positionnement stable des sondes, minimiser la compression externe, utiliser des règles de placement prédéfinies par région d’intérêt, et réacquérir des mesures lorsque des fluctuations marquées d’image à image ou de vides de signal sont des stratégies de dépannage importantes pour améliorer la fiabilité de l’acquisition.

Plusieurs variables procédurales semblent particulièrement importantes pour l’évaluation reproductible de l’élastographie longitudinale. Cela inclut une pression constante du transducteur, un positionnement du cou stable, un timing standardisé d’acquisition, l’évitement des régions fortement calcifiées ou kystiques, et le maintien d’une profondeur comparable de la région d’intérêt lors des examens en série. La reproductibilité peut également dépendre du maintien de préréglages d’élastographie identiques et de la minimisation de la variation entre opérateurs grâce à une formation spécifique au protocole et des critères de contrôle qualité prédéfinis.

Plusieurs applications de recherche spécifiques découlent de ces résultats. Premièrement, l’élastographie peut aider à stratifier les patients selon des trajectoires de réponse ultérieures présentant une réponse biomécanique précoce sous-optimale, qui pourraient bénéficier d’une surveillance plus rapprochée, d’une imagerie complémentaire plus précoce ou d’une ablation supplémentaire. En conséquence, les seuils de rigidité dérivés d’un protocole ou d’une plateforme unique doivent être interprétés avec prudence jusqu’à ce qu’une standardisation multicentrique plus large soit disponible. Deuxièmement, les trajectoires de rigidité peuvent servir de points intermédiaires dans les essais cliniques de techniques d’ablation, de stratégie d’électrodes ou de méthodes procédurales complémentaires, réduisant potentiellement le temps nécessaire pour détecter des différences significatives entre approches. Troisièmement, l’élastographie peut aider à interpréter des cas ambigus où la réduction du volume est retardée mais les symptômes changent, ou où une réduction du volume se produit mais une raideur résiduelle suggère un remodelage incomplet. Quatrièmement, intégrer l’élastographie dans les modèles prédictifs pourrait améliorer la pronostication individualisée en combinant les indicateurs de changement précoce avec des déterminants établis tels que le volume de base, la composition et l’énergie délivrée. Les résultats multivariés de cette étude soutiennent la faisabilité d’une telle intégration, puisque les changements précoces du rapport de déformation et les changements SWE ont apporté des informations indépendantes au-delà des covariables conventionnelles. Ces résultats suggèrent la faisabilité potentielle d’intégrer les métriques d’élastographie dans les futurs modèles prédictifs ou de stratification par risque.

Il existe également d’autres approches pour étudier cette hypothèse, et elles méritent d’être considérées. L’imagerie basée sur la perfusion, incluant l’échographie avec contraste et l’évaluation de la vascularité, peut identifier le tissu résiduel viable et prédire la réponse à long terme en estimant la fraction effective^{traitée 26}. Des métriques morphologiques telles que le rapport d’ablation initiale et le volume de la zone d’ablation précoce ont également été utilisées pour prévoir le risque de rétrécissement et de retraitement à^{long terme 27,28}. L’imagerie transversale peut caractériser davantage les zones d’ablation et les complications dans certains cas. Les résultats rapportés par les patients, le soulagement des symptômes et l’amélioration esthétique peuvent fournir des mesures plus directes du bénéfice clinique et doivent être intégrés aux biomarqueurs d’imagerie. Dans ce contexte, l’élastographie doit être considérée comme complémentaire plutôt que comme exclusive. Une orientation future forte est la comparaison directe de l’élastographie avec les méthodes basées sur la perfusion et les substituts morphologiques au sein d’une même cohorte, en utilisant des calendriers de suivi harmonisés et des définitions standardisées des résultats. Bien que la présente étude ait démontré des associations significatives entre le changement précoce de l’élastographie et la VRR ultérieure, la cohorte actuelle n’a pas été spécifiquement conçue ni propulsée pour l’élaboration de valeurs de seuil cliniquement exploitables avec des estimations validées de sensibilité et de spécificité. Les études prospectives multicentriques avec des analyses de seuil prédéfinies basées sur les résultats seront donc importantes pour la future traduction clinique. Enfin, le domaine bénéficierait d’études prospectives dans lesquelles les décisions de gestion précoces guidées par l’élastographie sont testées selon les stratégies de suivi standard, avec des résultats incluant le VRR, l’amélioration des symptômes et de l’esthétique, ainsi que les taux de retraitement.

Les limites de cette étude doivent être prises en compte lors de l’interprétation des résultats et de la planification des travaux ultérieurs. Premièrement, l’élastographie est intrinsèquement sensible aux conditions d’acquisition, notamment la pression du transducteur, la profondeur, la région d’intérêt et la reconstruction spécifique à chaque dispositif, et ces facteurs peuvent introduire une variabilité de la mesure malgré la protocolisation. Pour réduire les risques de confusion causés par des anomalies diffuses de raideur thyroïdienne, les patients présentant une maladie thyroïdienne diffuse cliniquement significative et une atteinte parenchymatose hétérogène diffuse ont été exclus lors de l’inscription. L’analyse formelle de la reproductibilité inter-observateur n’a pas été réalisée dans la présente étude, bien que des protocoles d’acquisition standardisés et des procédures de formation des opérateurs aient été appliqués pour améliorer la cohérence des mesures. Deuxièmement, le calendrier de suivi capture des étapes clés mais peut ne pas caractériser pleinement l’évolution très précoce au-delà de l’intervalle immédiatement post-ablation, durant lequel l’œdème et les changements réactifs peuvent influencer la raideur et compliquer la séparation de l’effet aigu du traitement de la physiologie post-procédurale transitoire. Troisièmement, la classification de la réponse basée sur un seuil de VRR, bien que largement utilisée, n’englobe pas toutes les dimensions du succès clinique, car l’amélioration des symptômes, les résultats esthétiques et la repousse à long terme ne suivent pas toujours parfaitement le VRR. De plus, la durée actuelle du suivi a été limitée à 12 mois, tandis que les nodules thyroïdiens traités peuvent continuer à se remodeler, à réduire ou à présenter une recroissance retardée sur des périodes d’observation plus longues. Bien que la repousse ouverte nécessitant une ablation ou une chirurgie répétée n’ait pas été observée durant l’intervalle de suivi actuel, les futures études avec suivi longitudinal prolongé devraient intégrer des schémas de repousse, des taux de retraitement et des exigences d’ablation répétée comme mesures d’issue cliniquement pertinentes supplémentaires. Quatrièmement, bien que les modèles multivariés s’ajustent aux principales covariables de référence et procédurales, la confusion résiduelle reste possible, y compris des nuances techniques non mesurées telles que la complétude des marges, la dynamique du positionnement des électrodes et l’expérience de l’opérateur, qui peuvent affecter à la fois le changement précoce de rigidité et la réponse à long terme. De plus, comme toutes les procédures ont été réalisées par des opérateurs expérimentés de RFA thyroïdien dans un seul centre, la généralisation externe à des réglages de faible volume nécessite une évaluation supplémentaire. Cinquièmement, l’étude n’a pas inclus de comparaisons directes avec d’autres outils de réponse précoce, tels que l’imagerie basée sur la perfusion ou les métriques de viabilité formelle de la zone d’ablation, au sein du même protocole, de sorte que la valeur incrémentale de l’élastographie par rapport à ces méthodes ne peut pas être pleinement établie ici. Comme l’étude était basée sur l’observation et l’imagerie, les associations identifiées ne doivent pas être interprétées comme une preuve d’une causalité mécaniste directe entre les modifications de l’élastographie et la biologie tissulaire sous-jacente. De futures études pourraient également évaluer individuellement la performance prédictive indépendante de l’élastographie de base, les mesures de raideur immédiatement post-ablation et l’élastographie de suivi d’un mois, ainsi que la variation relative en pourcentage entre les valeurs de base et immédiatement post-RFA. De telles analyses peuvent aider à clarifier si les mesures de rigidité absolue ou les métriques dynamiques de variation de rigidité offrent une plus grande utilité pronostique pour l’évaluation à long terme de la réponse volumétrique après une RFA thyroïdienne.

L’évaluation basée sur l’élastographie offre une approche réalisable et standardisée par protocole pour l’évaluation longitudinale après une ablation par radiofréquence par nodule thyroïdien. Les premiers changements de rigidité post-ablation étaient indépendamment associés à une réponse volumétrique à 12 mois et peuvent fournir des informations cliniquement utiles avant qu’une réduction substantielle du volume ne devienne apparente. Ces résultats soutiennent le rôle potentiel des protocoles standardisés d’élastographie comme biomarqueurs complémentaires d’imagerie pour l’évaluation de la réponse post-ablation et le suivi longitudinal après l’RFA thyroïdienne.

Matériaux

Liste des matériaux utilisés dans cet article

Nom	Entreprise	Numéro de catalogue	Commentaires
3M Tegaderm Film Dressing	Tegaderm	https://www.3m.com/3M/en_US/p/c/medical/bandages-dressings/film/i/health-care/medical/	Pansement en film adhésif; Couverture du site post-intervention.
Seringue à usage unique	Becton Dickinson	BD 309628	Hydrodissection et administration d'anesthésique local
Électrodes de mise à la terre pour RFA	3M Healthcare	Modèle 9547	Électrodes de retour à usage unique compatibles avec RF.
Transducteur à matrice linéaire haute fréquence (5–14 MHz)	SuperSonic Imagine	Sonde linéaire SL15-4	Échographie thyroïdienne et acquisition par élastographie
Transducteur à matrice linéaire ultrasonique à haute fréquence (7–15 MHz)	Mindray	L11-3U	Sonde linéaire optimisée pour l'imagerie thyroïdienne et l'élastographie.
IntelliVue	Philips Healthcare	MX450	Système de surveillance du patient; ECG, SpO2, moniteur de tension artérielle non invasif pendant la RFA.
IntelliVue	Philips Healthcare	MP50	Moniteur de signes vitaux; Surveillance continue péri-procédurale
Injection de chlorhydrate de lidocaïne	Hospira Inc.	NDC 0409-4276-17	Anesthésie locale avant l'ablation
Package lme4	CRAN	RRID:SCR_015654	Modélisation à effets mixtes linéaire en R.
Package lmerTest	CRAN	RRID:SCR_015656	Tests de modèles mixtes et valeurs p dans R.
Gants sans poudre Microflex	Ansell Healthcare		Barrière aseptique pour l'équipe de procédure.
Gel pour échographie non stérile	EcoGel		Médium de couplage pour l'imagerie de base et de suivi.
Système d'archivage et de communication d'images (PACS)	GE Healthcare	Centricity PACS	Système de stockage et de récupération d'images de qualité professionnelle.
Povidone-iode 10%	Betadine		Antiseptique pré-procédural.
Station de travail Precision	Dell	T5820	Pour la révision d'images et l'analyse de données.
Aiguille PrecisionGlide 22G–27G	BD (Becton Dickinson)	https://www.bd.com/en-us/products-and-solutions/products	Ensemble d'aiguilles pour l'anesthésie locale/hydrodissection.
Sac de sérum physiologique pré-mélangé 500 mL (0,9% NaCl)	Baxter Healthcare		Fluide pour l'hydrodissection et le rinçage.
Logiciel de calcul statistique R	R Foundation	RRID:SCR_001905	Environnement d'analyse statistique.
Électrode RFA percutanée refroidie	STARmed Co., Ltd.	Électrode Star RF 18-10s10F	Électrode RFA active pour l'ablation thyroïdienne (TaeWoong Medical USA)
Environnement de développement intégré RStudio	Posit (anciennement RStudio)	RRID:SCR_000432	IDE pour les analyses statistiques R.
Module logiciel d'élastographie par ondes de cisaillement	Mindray	Module SWE (série RESONA)	Cartographie quantitative de la rigidité en kPa dans le système US.
Kit de drapeau stérile	Medline Industries	DYNJP2500	Préparation stérile pour la procédure RFA thyroïdienne
Gaze et tampons stériles	Covidien	Kit d'éponges stériles	Utilisé pour la préparation de la peau et l'hémostase.
Sérum physiologique stérile	Baxter Healthcare	2B1324X	Hydrodissection pendant la RFA thyroïdienne
Draps de procédure stériles (Pack de draps bleus)	Medline Industries		Drap standard pour la configuration RFA percutanée stérile.
Gel pour échographie stérile	Aquasonic 100	Gel conducteur 100 mL	Médium de couplage stérile pour l'orientation par US.
Housse de sonde à ultrasons stérile	Parker Laboratories	Modèle 1001	Housse stérile pour l'imagerie intra-procédurale.
Module logiciel d'élastographie par déformation	Mindray	Module SE (série RESONA)	Calcul du rapport de déformation avec indicateur de qualité intégré.
Seringues (5/10/20 mL)	BD (Becton Dickinson)	Seringues Monoject	Pour l'administration d'anesthésique et d'hydrodissection.
Robinet à trois voies	BD (Becton Dickinson)	Modèle 305129	Robinet et tuyauterie pour une gestion fluide contrôlée.
Électrode RF thyroïdienne	STARmed Co., Ltd.	https://en.starmed4u.com/	Électrode RF refroidie en interne; Ablation par radiofréquence des nodules thyroïdiens
Gel de couplage pour échographie	Parker Laboratories	Aquasonic 100	Acquisition par ultrasons et élastographie
Système d'imagerie par échographie avec capacité d'élastographie	SuperSonic Imagine	Aixplorer	Utilisé pour l'échographie en mode B, l'imagerie Doppler, la SWE et l'élastographie par déformation
Générateur RF VIVA Combo	STARmed Co., Ltd.		Générateur RF pour l'ablation thyroïdienne percutanée (TaeWoong Medical USA)
Générateur RF VIVA	STARmed Co., Ltd.		Ablation par radiofréquence thyroïdienne guidée par échographie
Xylocaine 1% 10 mL	AstraZeneca		Anesthésique local (lidocaïne) pour l'analgésie procédurale.