갑상선 결절의 고주파 소작술 후 치료 반응에 대한 탄성촬영 기반 평가

Summary

갑상선 결절 고주파 소작술 전후에 표준화된 전단파 및 변형 탄성학 추적 프로토콜을 적용하여 강성 변화를 종단적으로 정량화하였습니다. 혼합 효과 모델링과 다변수 회귀분석을 사용하여, 초기 절제 후 탄성학 변화는 12개월 부피 감소 반응과 독립적으로 연관되어 있었습니다.

Abstract

고주파 절제술(RFA)은 양성 갑상선 결절에 대한 확립된 최소 침습 치료법이지만, 반응 평가는 주로 지연된 부피 감소에 의존하여 초기 시술 후 의사결정이 제한됩니다. 이 회고적 코호트에서 단일 지표 갑상선 결절을 가진 참가자들은 기저선, 절제 직후, 그리고 1, 3, 6, 12개월에 표준화된 B모드 초음파와 탄성촬영술을 받았습니다. 전단파 탄성학(SWE)의 강성과 변형률은 미리 지정된 관심 영역과 품질 관리 기준을 사용하여 획득되었습니다. 12개월 반응은 50% ≥ 부피 감소 비율(VRR)으로 정의되었습니다. 종단 변화는 선형 혼합 효과 모델을 사용하여 평가했으며, RFA 직후 탄성학 변화(즉각적으로 기저선 기준)는 상관관계 및 다변수 회귀분석을 기준 결절 부피, 결절 조성, 전달 에너지, 연령을 조정하여 12개월 VRR과 연관시켰습니다. 187명의 참가자 중 146명이 응답자였고 41명이 비응답자였습니다. 반응자 31.5 ± 5.1 kPa에서 41.7 ± 5.6 kPa, 비응답자(35.1 ± 5.3 kPa에서 43.4 ± 5.8 kPa)에서 SWE가 즉시 증가했으며, 응답자 중 12개월 SWE가 비응답자보다 낮아졌습니다(21.6 ± 4.7 대 31.2 ± 5.2 kPa). VRR은 조기 분리되었으며, 응답자 12개월 VRR은 73.6± 6.9%였으며, 비응답자는 38.5± 8.7%였습니다. 혼합 효과 모델에서는 즉각적인 SWE 증가가 유의미했으며(추정치 +10.22 kPa, P < 0.001), 후기 그룹 분기가 나타났다(비응답×자 12개월 추정치 +9.41 kPa, P < 0.001). 초기 탄성학 변화는 12개월 VRR과 상관관계가 있었으며(Δstrain ratio: r = 0.445, P < 0.001; ΔSWE: r = 0.302, P < 0.001), 조정 모델에서는 VRR과 독립적으로 연관되어 있었습니다(Δstrain ratio β = 14.80% / +1.0, P < 0.001; ΔSWE β = 1.80% / +5 kPa, P = 0.011). RFA 직후 탄성검사 변화는 초기 조직 효과를 포착하고 12개월 체적 반응에 대한 독립적인 예후 정보를 제공합니다.

Introduction

갑상선 결절은 현대 내분비계 진료에서 '검출 질환'으로 정의되는 것이 되었는데, 이는 그 생물학적 특성이 갑자기 변했기 때문이 아니라, 영상 촬영 덕분에 갑상선이 매우 명확하게 보였기 때문입니다 ^1,2,3. 측정된 유병률은 결절을 찾는 방법에 크게 좌우되는데, 고전적인 역학 종합에서는 촉진 시 2%–6%, 초음파 약 19%–35%, 부검 시리즈에서 약 8%–65%까지 추정되며, 이는 방법 민감도가 겉보기 역학을 어떻게 재구성하는지를 보여줍니다^. 요오드 결핍 환경과 고주파 변환기가 있을 때는 초음파 조사가 더 높은 검출률을 보고할 수 있습니다; 예를 들어, 13MHz 기술을 사용한 인구 기반 조사에서는 성인의 68%에서 결절이 발견되었으며, 이는 장비의 영향과 갑상선이 생애 전반에 걸친 결절 변화에 대한 근본적인 취약성을 강조했다². 중국은 초음파가 건강 검진에 자주 포함되기 때문에 규모에 대한 생생한 시각을 제공합니다. 6,985,956명의 참가자가 참여한 국가 건강검진 코호트에서 갑상선 결절의 전체 유병률은 36.9%였으며, 연령 및 성별 표준화 유병률은 38.0%, 여성 우세가 뚜렷했습니다^. 대부분의 병변이 임상적으로 무기력하더라도, 반복적인 영상 촬영, 절차적 연쇄 반응, 증상 경계 강화, 안심과 개입 사이의 지속적인 긴장 등 전체적인 결과는 상당합니다. 이 긴장은 증상이 있거나 미용적으로 눈에 띄거나 명백히 커지는 양성 결절에서 가장 두드러지며, 여기서 '양성'이 '중요하지 않다'는 의미가 아니며, 환자 중심의 결과가 치료가 국소 효과를 완화하면서 기능을 유지하고 수술적 이환을 피할 수 있는지에 달려 있습니다.

지난 10년간 초음파 유도 열절제술은 특수적 채택에서 점점 확립된 수술 대안으로 자리 잡았으며, 특히 양성 고형 또는 주로 고형 결절, 그리고 신중하게 선택된 맥락에서는 재발성 병변에 대해 수술의 대안으로 자리 잡았습니다 ^4,5. 고주파 절제술(RFA)은 실시간 초음파 시각화와 제어 가능한 에너지 전달, 그리고 경험 많은 손에 유리한 합병증 프로필을 결합하여 대표적인 기법이 되었습니다 ^4,5. 많은 경우 경추 흉터와 전신마취를 피하는 명백한 매력 외에도, RFA는 갑상선 예비력을 보존하면서 병변을 재형성하는 장기 보존 철학과 일치합니다. 현대의 결과 보고는 일반적으로 6개월에서 12개월에 평가되는 부피 감소 비율(VRR)에 수렴하며, 치료 성공률은 여러 시리즈에서 이루어졌으며, 1년⁶개월에 최소 50% 용피 감소와 같은 임계값을 사용하여 통합 분석이 정의되었습니다. 하지만 볼륨 감소는 본질적으로 늦은 대체물입니다. 소작술의 생물학적 효과는 단백질 변성과 응고성 괴사로 즉시 시작되며, 이후 염증, 흡수, 섬유화를 거쳐 몇 주에서 몇 개월에 걸쳐 진행됩니다^. 그 결과, 초기 절제 후 조직 반응에 관한 추가 정보를 제공하는 임상적으로 가장 초기 결정적인 질문은 크기 ^4,6만으로는 확신을 가지고 답하기 어려운 경우가 많습니다. 도플러 혈관과 회색조 반향 질감은 유용한 맥락을 제공하지만, 여전히 조직 변형의 간접 측정에 불과하며, 수술 후 초기 시기에는 반응성 괴혈, 부종, 이질성 괴사 변화로 인해 혼동될 수 있습니다. 적절한 치료의 조기 객관적 신호는 추적 치료 강도를 합리화하고, 체피 반응 지연 가능성을 예측하며, 수개월간의 관찰 대기가 지속 증상이나 재발배로 확고화되기 전에 덜 우호적인 종단 반응 패턴을 가진 환자를 식별하는 데 도움을 줄 수 있다 ^2,3.

탄성학은 이 환경에서 매력적인데, 이는 절제술이 바꾸려는 물리적 특성을 측정하기 때문입니다. 초음파 탄성술은 변형 기반 및 전단파 기법을 포함하며, 변형 패턴을 통해 조직 강성을 반정량적으로, 전단파 속도와 탄성 계수⁷을 통해 정량적으로 특성화합니다. 갑상선 영상에서 탄성촬영은 악성 위험 계층화를 위한 보조 요법으로 광범위하게 평가되어 왔습니다. 그러나 엄격한 검사는 광범위한 임상 집단에서 기존 초음파보다 점진적 진단 가치를 높여준다는 초기 열정을 누그러뜨렸으며, 우아한 물리학이 모든 맥락에서 임상적 이점을 보장하지 않는다는 점을 상기시켰습니다⁸. 절제 후 문제는 다릅니다. 여기서는 표적 병변이 알려져 있고 치료됩니다. 이 과제는 조직학적 구별이 아니라 반응 평가이며, 핵심 생물학적 전환은 생존 가능하고 변형 가능한 조직에서 덜 순응적인 괴사성 및 섬유화 복합체로의 전환입니다. 열 모달리티 전반에 걸친 초기 임상 연구들은 이 기전적 전제를 뒷받침합니다. 마이크로파 절제술을 받은 환자에서 엘라스그래피 기반 점수 측정은 치료 후 측정 가능한 구조적 변화를 보여주었으며, 실현 가능한^{추적 관찰 도구}로 제안되었습니다. RFA 코호트에서는 탄성학이 치료된 결절과 인접한 기준 조직의 강성 변화를 정량화하는 방법으로 탐구되어 있으며, 이는 거시적 수축 이전에 발생할 수 있는 소작 유도 재형성을 잠재적으로 엿볼 수 있는 창을^{제공합니다.} 기술적 학습 곡선 연구에서는 초음파 특성화에 변형 탄성술이 포함되어, 강성 패턴이^{크기 11}을 넘어 정보를 더할 수 있다는 인식이 점점 커지고 있지만 여전히 일관되지는 않다. 이 연구들을 종합해 보면, 탄성학이 진단 보조 도구에서 반응 바이오마커로 발전할 수 있음을 시사합니다. 그러나 프로토콜 표준화 부족과 획득 방법의 변동성은 여전히 광범위한 임상 적용에 중요한 장애물로 남아 있습니다^10,11. 기존 보고서들은 탄성학 모달리티, 획득 매개변수, 관심 영역 전략, 참조 조직 선정, 평가 시기 등 다양한 면에서 차이가 있어, 강성 변화를 임상적으로 견고한 규칙으로 전환하여 운영자와 플랫폼 간 초기 의사결정을 안내하는 데 어려움을 줍니다.

이 연구는 갑상선 결절의 RFA 후 치료 반응을 평가하기 위한 탄성학 기반 접근법을 평가하고, 일상 추적 관찰에서의 실용적 적용에 명시적으로 주의를 기울여 이러한 중개 격차를 해소하기 위해 설계되었습니다. 본 연구의 핵심 전제는 RFA 후 반응이 근본적으로 열손상에서 시작하여 절제 후 종방향 구조 및 강성 변화를 거치는 생체역학적 과정이며, 따라서 강성 지표는 방향성 있고 임상적으로 해석 가능한 방식으로 변화해야 한다는 것입니다. 탄성소파 유도 지수가 치료 효과와 일치하는 체계적인 소작 후 이동을 보여주고, 추적 관찰 기간 동안 반응 경로를 추적하며, 특히 크기 기반 평가가 가장 정보 제공이 적은 초기 단계에서 VRR 및 도플러 혈관성과 같은 기존 평가 지표에 보완적인 정보를 제공할 것이라는 가설이 세워졌다 ^4,6,7. 본 연구의 새로움은 악성 분별을 위한 추가 특징에서 엘라스그래피를 객관적이고 정량적으로 측정 가능한 사후 조직 변형 판독으로 재구성하고, RFA 이후 관리에 영향을 미치는 임상 질문들, 즉 종단 치료 반응 평가 및 추적 관찰 중 절제 후 강성 변화의 특성화와 일치시키는 데 있습니다. 탄성촬영이 이러한 질문들을 신뢰성 있게 운영화할 수 있다면, 중재와 결과 간의 피드백 루프를 단축하여 RFA 이후 종단 반응 패턴의 조기 특성화를 가능하게 할 수 있으며, 양성이지만 부담이 많은 결절 환자들에게 갑상선 RFA가 점점 더 중요한 선택이 된 최소 침습적 이점을 유지할 수 있습니다 ^4,6.

Protocol

윤리 선언
인간 참여자가 포함된 모든 절차는 싱타이 병원의 기관 지침에 따라 수행되었으며, 헬싱키 선언 및 이후 개정 원칙을 준수하였습니다. 이 연구는 싱타이 병원 인간 연구 윤리 위원회의 검토 및 승인을 받았으며; 모든 참가자가 일상 진료의 일환으로 RFA 및 초음파 검사에 대해 표준 서면 동의서를 제공했기 때문에, 위원회는 추가 서면 연구 동의를 포기하고 필요 시 구두 동의를 허용하며 연구를 승인했습니다. 참가자 기밀성은 안전하고 접근 제한 데이터 관리를 통해 보호되었으며, 분석 전 모든 연구 데이터의 식별 해제가 이루어졌습니다.

연구 설계 및 전체 도식도
본 연구는 갑상선 결절의 초음파 유도 RFA 후 치료 반응에 대한 탄성촬영 기반 평가를 평가하기 위한 후향적 단일 기관 임상 코호트로 설계되었습니다. 우리 기관에서 연속적으로 RFA를 받는 적격 환자들을 선별하여 사전 정해진 목표인 187명까지 등록하였습니다. 분석적 독립성을 보장하고 환자 내 군집화를 피하기 위해, 참가자당 한 개의 지시 결절만 포함하였으며, 사전에 지시 절제술 시 치료 대상이 선택된 결절로 정의되었습니다. 배경 갑상선 경직성 측정에 상당한 영향을 미칠 수 있는 확산성 갑상선 질환, 예를 들어 뚜렷한 이질성 실질 침범이 있는 하시모토 갑상선염, 조절 불가능한 그레이브스병, 확산성 갑상선염, 또는 확산성 갑상선종 등이 포함된다. 초음파에서 확산성 구조적 갑상선 이상이 없는 고립된 갑상선 자가항체 양성은 배제 기준으로 간주되지 않았습니다.

연구는 표준화된 종단적 워크플로우를 따랐습니다. 기초 임상 평가 및 영상 검사는 절제술 2주 이내에 수행되었으며, 기존 초음파, 도플러 평가, 사전 정해진 프로토콜에 따른 탄성촬영 획득이 포함되었습니다. RFA는 기관 기준에 따라 실시간 초음파 지침 하에 진행되었습니다. 에너지 전달 후 30분 이내에 즉시 영상 촬영을 통해 기술적 성공과 급성 조직 변화를 기록하였습니다. 정기적 추적 방문은 절제술 후 1, 3개월, 6개월, 12개월에 실시되었으며, 각 시점에서 동일한 획득 설정과 측정 전략을 사용하여 반복 초음파와 탄성촬영을 수행하여 기술적 변동성을 최소화하였습니다.

표본 크기 계산 및 통계적 가정
표본 크기 추정의 주요 평가 평가 대상은 절제 후 1개월 시점의 정량적 탄성검사 변화로, 표준화된 관심 영역 전략을 사용해 평균 강성(예: 탄성 계법 kPa 단위)의 차이로 표현되었습니다. 1차 분석은 동일한 결절 내 쌍 측정값을 비교하기 때문에, 표본 크기는 쌍 평균 차이 프레임워크를 사용하여 추정되었습니다. Δ 를 기준선에서 1개월까지의 강성 평균 변화를 나타내고, σ_d 를 쌍 차이의 표준편차로 정의한다. 유의수준 α 에서 멱 1−β의 양측 검정에 대해 필요한 표본 크기는 다음과 같습니다:

α = 0.05(양면)와 거듭제곱 1−β = 0.80, Z1−α/2 = 1.96, Z1−β = 0.84를 사용한다. 기초 측정과 갑상선 탄성조영술의 소작 후 변동성 발표를 바탕으로, 본 연구는 _{β d} = 35 kPa, 보존적 임상적으로 의미 있는 변화 Δ = 8 kPa를 가정하였습니다. 이로 인해 RFA 전후 탄성학 강성에 대한 1차 쌍 분석에서 최소 150명의 참가자가 필요합니다.

절제 전 기초 임상 평가
등록 시점에는 결절 관련 증상과 시술 안전성에 중요한 요인에 중점을 둔 표준화된 임상 병력이 작성되었습니다. 나이, 성별, 인체 측정 수치, 그리고 심혈관 및 폐 질환, 당뇨병, 만성 신장 질환, 이전 경부 수술 또는 방사선 치료 등 관련 동반 질환이 기록되었습니다. 현재 복용 중인 약물을 문서화했으며, 특히 항혈소판제와 항응고제에 주의를 기울였으며, 기관 정책에 따른 계획된 수술 후 조정도 기록하였습니다. RFA의 적응증이 확인되었으며, 주요 원인이 압박 증상, 미용 문제, 연쇄 영상에서 기록된 성장 또는 이 요인들의 조합인지 명시되었습니다.

기초 증상 부담은 미리 지정된 기기를 사용하여 확립되었습니다. 연하곤란, 목 압박감, 이물질 감각, 호흡곤란, 음성 관련 불만과 같은 국소 압박 증상을 평가했으며, 그 심각도는 수치 평가 척도 또는 기관에서 채택한 증상 점수를 사용해 문서화되었습니다. 미용 영향이 지표일 경우, 미용 외관은 표준화된 미용 척도로 평가되었고, 환자 보고된 불만족도 기록되었습니다. 결절 인식 기간, 이전의 미세주사 흡인 결과, 이전 양성 세포학 확인 전략(필요에 따라 단일 샘플링 vs 반복 샘플링), 갑상선 질환에 대한 이전 약물 치료 등 관련 갑상선 병력을 기록하였습니다.

기초 실험실 검사는 일상적인 임상 관행에 따라 이루어졌으며, 갑상선 자극 호르몬 및 필요한 경우 추가 갑상선 기능 지표를 포함하였습니다. 지시 결절의 양성 상태를 확인하는 세포학 및 병리학 보고서와 가능한 경우 초음파 위험 층화 문서가 검토 및 기록되었습니다. 결절이 기관에 부딪히거나 기도 손상을 시사하는 증상을 일으켰을 때 생체 징후와 집중 기도 평가를 포함한 기초 안전성 평가가 확립되었습니다. 계획된 추적 관찰 일정이 문서화되었고, 탄성촬영 중 삼킴 방지를 포함한 표준화된 영상 전 지침이 강화되어 방문 간 측정 일관성을 개선하였습니다.

기존 초음파 획득 및 결절 부피 측정
환자를 누운 채 목을 부드럽게 펴는 상태에서 회색 초음파를 시행했습니다. 표준화된 영상 프리셋이 사용되었고, 주파수 범위, 깊이, 초점, 게인, 다이내믹 레인지 등 주요 설정은 방문마다 일관되게 유지되었습니다. 프레임 속도, 탄성학 척도 범위, 지속도, 평탄화 설정은 기저 및 추적 검사 전반에서 일정하게 유지되었습니다. 초음파 검사는 전단파 및 변형 탄성술 기능이 탑재된 고주파 선형 배열 변환기(5–14 MHz)를 사용하여 수행되었습니다. 모든 초음파 및 탄성학 검사는 연구 기간 내내 동일한 초음파 플랫폼과 소프트웨어 버전을 사용하여 강성 측정의 시스템 간 변동성을 최소화하였습니다. 결절 가장자리와 인접한 해부학적 랜드마크를 명확히 구분한 횡방향 및 종방향 사진이 획득되었습니다. 집정 결절의 최대 직교 직경—전후(A), 횡단(B), 종(C)—을 측정하였다. 전자 캘리퍼스는 결절의 외곽 가장자리에 부착되어 주변 실질이나 주위 후광을 포함하지 않도록 했습니다. 모든 측정값은 밀리미터 단위로 기록되었으며, 가능할 경우 각 치수를 반복하고 평균값을 사용했습니다. 결절 부피는 타원체 공식을 사용하여 계산되었습니다:

동일한 측정 방식을 각 추적 방문에서 반복하였다. VRR은 다음과 같이 계산되었습니다:

탄성촬영 획득 및 품질 관리
엘라스그래피는 환자를 누운 상태로 목을 뻗은 상태로 기존 초음파 검사 직후에 시행되었습니다. 환자는 삼키지 말고 획득 시 잠시 호흡을 멈추어 움직임 아티팩트를 최소화하도록 요청받았습니다. 모든 방문에서 동일한 탄성학 모드와 설정이 사용되었습니다.

전단파 탄성술에서는 샘플링 박스가 인덱스 결절의 고체 부분과 가능할 경우 유사한 깊이의 인접한 정상 갑상선 조직을 포함하도록 배치되었습니다. 관심 있는 원형 영역은 낭포성 퇴행, 대석회화, 말초 후광 영역, 명백한 소작 후 기동 변화를 피하면서 결절의 중심 고체 구성 내에 우선적으로 배치되었습니다. 관심 영역 배치 후, 엘라스토그래피 색상 지도가 최소 3–5초 동안 시각적으로 안정적일 때만 이미지 획득이 시작되었습니다. 주위 강성 평가는 치료 결절 자체의 종방향 강성 변화에 초점이 맞춰져 있었기 때문에 사전 정의된 획득 프로토콜에 체계적으로 통합되지 않았습니다. 최소 세 가지 기술적으로 적절한 측정이 이루어졌으며, 중앙 강성 값이 주요 지표로 기록되었습니다; 사분위 범위는 가능할 때 문서화되었습니다.

변형 탄성술에서는 시스템의 품질 지표에 따라 부드럽고 일관된 압축이 적용되었습니다. 압축은 느리고 리드미컬한 주파수로 수동으로 가해졌으며, 초음파 시스템이 제공하는 실시간 품질 지표를 모니터링했습니다. 갑상선 캡슐의 눈에 띄는 변형이나 불안정한 탄성 색 매핑을 유발하는 과도한 프로브 압력은 피하였습니다. 인접한 정상 갑상선 조직에서 유사한 깊이의 기준 관심 영역을 정의하고 긴장 비율을 계산하였다. 각 탄성학 획득 시퀀스는 이미지 안정화와 재현 가능한 샘플링을 보장하기 위해 약 5–10초간 유지되었습니다. 최소 세 가지 허용 가능한 측정값이 확보되었고, 중앙값이 기록되었습니다.

측정은 탄성 지도가 안정적이고, 관심 영역 위치가 미리 정해진 규칙을 따르며, 대운동이나 압축 아티팩트가 없을 때만 허용되었습니다. 불완전한 색상 채우기, 관심 영역 내 명백한 신호 빈틈, 또는 프레임 간 뚜렷한 변동이 있는 측정은 거부되고 재획득되었습니다. 초음파 및 탄성촬영 검사는 조직 강성 측정에 대한 환경 영향을 최소화하기 위해 온도 조절된 진료실에서 수행되었습니다. 각 시점마다 해당 B-모드 이미지와 대표적인 탄성 이미지가 추적성과 품질 검토를 위해 저장되었습니다.

RFA 절차
RFA는 갑상선 소작술 기법에 대한 기관 교육을 마쳤고, 갑상선 초음파 유도 중재 분야에서 5년 이상의 경험을 가진 운영자들이 연속 초음파 지도 하에 수행했으며, 연구 시작 전 각각 200건 이상의 갑상선 RFA 절차를 수행한 경험이 있었습니다. 환자는 목을 펴고 누운 자세로 있었다. 경추 부위는 표준 멸균 기법으로 준비되고 덮였습니다. 국소 마취는 일상적인 절차에 따라 시행되었으며; 임상적으로 필요 시 의식진정을 사용하고, 활력징후를 지속적으로 모니터링하였다.

고주파 전극은 실시간 초음파 유도 하에 인덱스 결절에 삽입되었으며, 전극 안정성을 높이고 의도치 않은 변위를 최소화하기 위해 트랜지스믹 방식을 선호했습니다. 검지 결절이 기관, 경동맥 피, 식도 또는 반발 후두신경 부위에 인접해 있을 때, 에너지 전달 전반에 걸쳐 소작부와 인접 중계 구조물 사이의 가시적 분리를 유지하기 위해 초음파 유도 하에 멸균액을 천천히 주입하여 수액 절리술을 시행했습니다.

고주파 에너지는 결절의 가장 깊은 부분에서 시작해 체계적으로 표면 영역으로 이동하여 겹치는 소작 부위를 달성하는 무빙 샷 기법을 사용해 전달되었습니다. 각 소작 유닛은 일시적인 고반 초음파 변화가 목표 부위를 차지할 때까지 처리한 후, 전극 끝을 인접한 미처리 부위로 재배치했습니다. 출력 출력은 실시간 반향 변화와 환자 내성에 따라 조정되었고, 도플러 평가에서 목표 결절이 일시적이고 과도한 고반향 변화를 보일 때까지 절제술을 계속했습니다. 또한 명확한 잔류 결절 내 혈관 신호 없이 표적 치료 단위를 덮는 연속적인 고에코 영역의 시각화로 충분한 절제술이 지원되었습니다.

소작 시간, 전력 범위, 추정 전달 에너지, 수압 절부 사용, 시술 중 사건 등 주요 절차 매개변수가 기록되었습니다. 소작술이 완료된 후, 환자는 기관 프로토콜에 따라 관찰되었고, 표준화된 사후 지침이 제공되었습니다. 환자들이 심각한 통증, 목소리 변화, 통제 불가능한 기침, 진행성 부종 또는 열손상이나 시술 합병증에 우려를 불러일으키는 기타 소견이 나타날 경우 절제술은 일시적으로 중단되거나 중단되었습니다.

소작 직후 영상 및 문서 기록
RFA 완료 후 30분 이내에 기저선과 동일한 영상 프리셋을 사용하여 즉시 시술 후 초음파 검사를 실시하였습니다. 회색 이미지를 촬영하여 치료 부위의 범위를 기록하고, 일시적 초에코 발생성 및 내부 구조 변화 등 전형적인 소작 후 특징을 식별하였습니다. 색감 또는 배율 도플러를 사용해 소화부(ablated) 영역 내 혈류 감소 또는 부재를 확인하기 위해 외주혈관 및 주막혈관을 평가하였다.

일시적인 절차 관련 움직임이 해결되고 환자가 안정적인 호흡 정지기를 유지할 수 있게 된 후에도 기저라인에서 사용된 동일한 획득 프로토콜에 따라 엘라스그래피 획득이 반복되었습니다. 강성 값이 기록되었고, 대표적인 탄성학 지도와 해당 B-모드 이미지가 저장되었습니다. 획득 과정에서 발생한 환자의 불편함이나 운동 허위 등 기술적 한계는 문서화되었습니다. 이 소작 직후 이미지는 절차적 참고 자료로 사용되었으며 주요 종말점으로 사용되지는 않았으나 이후 종단적 변화에 대한 맥락을 제공했습니다.

추적 관찰 일정 및 표준화된 평가 지표
추적 방문은 절제 후 1개월, 3개월, 6개월, 12개월에 예약되었습니다. 각 방문마다 증상 평가, 미용 등급 평가, 기존 초음파, 기저 시와 동일한 획득 매개변수와 측정 전략을 사용한 탄성촬영 등 표준화된 평가가 수행되었습니다. 주요 탄성학 평가 지점은 사전 정해진 추적 시점에서 기준선 대비 지상 결절의 강성 변화였으며, 특히 1개월경의 조기 변화에 중점을 두었습니다. 주요 기존 영상 평가 평가변수는 12개월의 VRR이었습니다. 부차평가 평가 대상에는 이전 추적 방문에서의 VRR, 종방향 강성 궤적, 치료 부위 내 잔류 또는 재발성 혈관 조직의 존재가 포함되었습니다. 추적 관찰 중 부작용과 추가 개입은 체계적으로 기록하였다.

데이터 관리 및 블라인딩
각 참가자에게는 익명성을 보장하기 위해 고유 연구 식별자가 할당되었습니다. 임상 데이터와 영상 파일은 허가된 연구진만 접근할 수 있는 안전한 기관 데이터베이스에 저장되었습니다. 탄성학 및 체적 측정은 임상 증상 점수와 이전 탄성학 결과에 눈가림이 없는 훈련된 독자에 의해 수행됩니다. 일관성을 위해 관심 지역 배치 규칙, 이미지 품질 기준, 계산 절차를 명시한 서면 측정 매뉴얼을 유지하세요. 사전 정해진 기준을 초과하는 상당한 독자 간 불일치가 발생한 경우, 측정값을 공동 검토하여 합의하여 해결합니다. 데이터셋에 대한 모든 변경 사항은 감사 기록을 위해 기록됩니다. 모든 초음파 및 탄성촬영 획득은 갑상선 초음파 및 RFA 영상 평가 경험이 있는 운영자가 연구 시작 전에 프로토콜별 교육을 받은 상태에서 수행하였습니다. 저장된 이미지는 채취 시 적용된 동일한 사전 정의된 관심 영역 배치 규칙과 품질 관리 기준을 사용하여 오프라인에서 검토되었습니다.

통계 분석
연속 변수들은 히스토그램의 시각적 검사와 샤피로-윌크 검정을 통해 분포 특성을 평가하였다. 정규 분포 변수는 평균 ± 표준편차로 제시되며, 비정규 분포 변수는 중간값으로 요약되어 4분위 범위가 있습니다. 범주별 변수는 개수와 백분율로 보고됩니다. 주요 분석 목적은 RFA 후 결절 내 탄성학 유도 강성 변화를 평가하는 것이었습니다. 기본 시점과 절제 후 시점 간의 1차 쌍 비교에서는, 쌍 간 차이 분포가 정규성을 근사할 때 쌍 t검정을 사용해 강성 값 차이를 분석하였다; 그렇지 않으면 윌콕슨 부호 랭크 테스트가 적용되었다. 효과량은 평균 또는 중앙값 변화와 적절한 경우 95% 신뢰구간으로 보고됩니다. 모든 추적 방문(기저선, 절제 직후 방문, 1개월, 3개월, 6개월, 12개월)에서의 탄성학 측정 종단적 변화를 특성화하기 위해, 참가자 수준에서 무작위 절개를 적용한 선형 혼합 효과 모델을 구축하여 피험자 내 상관관계를 반영하였다. 시간은 연속 변수가 아닌 범주적 고정 효과로 모델링되었는데, 이는 사후 생체역학적 변화가 추적 방문 간 선형 경로를 따르지 않을 것으로 예상되었기 때문입니다. 특히, 연구 설계는 절제 직후 강성 급성 증가와 점진적 재형성 및 후기 감소를 예상하여, 선형성 가정을 강요하지 않고 비선형 시간적 반응 패턴을 포착하는 범주적 모델링이 더 적합하다고 판단했습니다. 잔류 플롯 검사를 통해 모델 가정을 평가하였다. 적절한 경우, 강성 값은 모델 적합도를 높이기 위해 로그 변환되었습니다. 2차 분석에서 탄성학에서 유도된 강성 변화와 기존 치료 반응 지표 간의 연관성을 탐구하였다. 초기 강성 변화와 이후 추적 관찰 시 VRR 간의 상관관계는 데이터 분포에 따라 Pearson 또는 Spearman 상관계수를 사용하여 평가하였다. 다변수 회귀 모델을 사용하여 초기 탄성학 변화가 기초 결절 부피, 결절 조성(고체 대 주로 고체 상태), 전달 에너지, 연령에 대한 조정 후 이후 부피 반응을 독립적으로 예측하는지 조사하였다. 하위 그룹 변수 간의 상호작용 항과 시간은 혼합 효과 모델링 프레임워크 내에서 테스트되었습니다. 종단 분석의 경우, 혼합 효과 모델은 무작위로 누락된 관측을 본질적으로 수용합니다. 모든 통계 검정은 양측 검사였으며, P 값 0.05< 통계적 유의성을 나타내는 것으로 간주되었습니다.

Results

연구 설계 및 전체 도식도
제외된 참가자들은 사전에 정해진 적격 기준을 충족하지 못했거나, 불완전한 기초 영상 데이터가 있었거나, 표준화된 종단 추적 관찰이 부족했습니다. 포함된 187명 참가자 중 1차 분석에 포함된 모든 사례에 대해 완전한 12개월 결과 분류가 제공되었으나, 개별 추적 방문에서 간헐적으로 누락된 관찰은 혼합 효과 모델링 프레임워크 내에 수용되었습니다(그림 1). VRR을 기반으로 한 사전 정해진 12개월 반응 정의를 사용하여, 146개의 결절이 반응자(VRR ≥ 50%), 41개가 비반응자(VRR < 50%)로 분류되었습니다(그림 1). 이 설계는 결절 내 쌍 비교(내부 대조군으로 기준선)와 사전에 정해진 추적 관찰 시점에 대한 표준화된 종단 평가를 가능하게 했습니다.

기초 참가자 및 결절 특성
기초 임상 특성은 표 1에 요약되어 있습니다. 비응답자는 응답자보다 나이가 많았으며(7.88세± 49.63 ± 8.01 대 45.94, P = 0.009), 기초 증상과 미용 부담이 약간 더 높았다(압력 증상 점수 3.43 ± 0.57 대 3.19 ± 0.63, P = 0.029; 미용 문제 점수 3.01 ± 0.47 대 2.79 ± 0.49, P = 0.011). 갑상선 기능은 그룹 간 대체로 유사했으나, 비반응자들 사이에서 TSH가 다소 높았다(1.86 ± 0.31 mIU/L 대 1.74 ± 0.33 mIU/L, P = 0.039)(표 1). 이러한 기저선 차이는 개입 전 비응답자들 중 임상 프로필이 더 불리했음을 나타냅니다. 기저 초음파에서는 비반응 환자에서 유의하게 더 큰 결절이 나타났으며, 최대 직경(33.87 ± 4.03 mm 대 29.31 ± 3.94 mm, P < 0.001)과 기저 부피(12.04 ± 2.41 mL 대 8.33 ± 1.90 mL, P < 0.001)가 포함된다(표 1). 결절 구성도 달랐는데, 비반응자들은 낭포성 구성 범주가 더 높은 쪽으로 이동하는 것으로 나타났습니다(P < 0.001)(표 1). 이러한 기본 초음파 특성은 이후 부피 및 생체역학적 반응 패턴에 대한 해부학적 맥락을 제공합니다. 어떤 참가자도 미상성 갑상선 질환의 사전 정해진 제외 기준을 충족할 만큼 심각한 확산성 갑상선 실질의 이상을 보였습니다.

탄성학 측정은 운동 및 압축 인공물을 최소화하기 위해 사전에 정해진 관심 영역 규칙과 수용 기준이 적용된 표준화된 조건 하에서 이루어졌습니다. 결절 내 기준선 측정을 내부 기준으로 사용하면 개인 간 변동성의 영향을 줄이고 종방향 강성 변화 해석을 지원할 수 있습니다. 절차적 특성은 표 1에 나와 있습니다. 더 큰 기초 결절 부담과 일치하게, 비반응자는 시술 지속 시간이 더 길었고(21.13 ± 3.01분 대 18.42 ± 2.79분, P < 0.001), 전달 에너지도 더 높았다(13.18 ± 2.63 kJ 대 11.04 ± 1.98 kJ, P < 0.001)(표 1). 이러한 절차적 차이는 기준선 크기와 조성 차이를 고려할 때 예상된 것이었으며, 반응 예측의 조정 분석에서 고려되었습니다.

RFA 이후 종방향 탄성학 궤적
RFA 직후, 두 가지 탄성검사 모두 급성 강성 증가를 나타냈습니다(표 2; 그림 2)는 열 응고의 즉각적인 생체역학적 효과라고 가정하는 것과 일치합니다. 전단파 탄성검사(SWE)는 응답자의 31.5 ± 5.1에서 41.7 ± 5.6 kPa로, 비응답자의 경우 35.1 ± 5.3에서 43.4 ± 5.8 kPa로 증가했습니다(표 2; 그림 2). 긴장 비율도 일치하여 증가했다(응답자 2.41 ± 0.37에서 3.14 ± 0.41; 비응답자 2.77 ± 0.40에서 3.35 ± 0.42) (표 2; 그림 2). 반면, 두 그룹 모두 즉시 절제술 후 부피는 기준선과 거의 일치하여(표 2), 탄성학이 단기간의 부피보다 치료 관련 조직 변화를 더 일찍, 조기에 영상으로 판단할 수 있다는 연구 가설을 뒷받침합니다.

종단 체적 결과 및 VRR
추적 관찰 중, 탄성학은 소작술 직후 상승 후 점진적인 재형성을 보여주었습니다(표 2; 그림 2). 응답자들의 경우 SWE는 기준선 3개월 이하로 떨어졌고 12개월 동안 계속 감소했습니다(기준선 기준치 31.5 ± 5.1 kPa에서 12개월 시 21.6 ± 4.7 kPa), 비응답자들은 감소가 약해졌고 12개월 기준 상당히 더 강하게 유지되었습니다(31.2 ± 5.2 kPa)(표 2; 그림 2). 긴장 비율도 같은 패턴을 보였으며, 12개월 때 반응자는 1.87± 0.35로 감소한 반면, 비응답자는 2.55± 0.41이었습니다(표 2; 그림 2).

체적 결과는 생체역학적 변화와 병행하여 진화했습니다(표 2; 그림 3). 두 그룹 모두 결절 부피가 꾸준히 감소했으며, 시간대에 걸쳐 지속적인 분리가 지속되었습니다; 12개월에 반응자는 2.20 ± 0.45 mL에 도달했으며, 비반응자는 7.35 ± 1.80 mL에 도달했습니다(표 2; 그림 3). VRR도 조기 분리되었고 구별되었습니다: 1개월 시점에 응답자 중 VRR은 31.5± 6.3%였으며, 비응답자는 16.1± 3.8%였고, 12개월 시점의 VRR은 각각 73.6± 6.9%± 38.5 8.7%였습니다(표 2; 그림 3). 이러한 정렬된 탄성학과 VRR 궤적은 RFA 이후 생체역학적 재모델링이 지속 가능한 체적 반응을 추적하고 분화한다는 가설을 뒷받침합니다.

12개월 체적 반응의 독립 예측 변수로서의 조기 탄성학 변화
혼합 효과 모델은 기술적 경로를 확인하고 시간에 따른 집단 차이를 정량화했습니다(표 3). SWE의 경우, 시간 효과는 기준선에 비해 즉각적으로 유의미한 증가를 나타냈습니다(추정치 +10.22 kPa, P < 0.001), 이후 방문에서는 점차 낮은 수치가 나타났습니다(예: 12개월 대비 기준 추정치 -9.02 kPa, P < 0.001). 중요하게도, 그룹별 시간별 상호작용에서는 후속 추적 관찰에서 차이가 나타났으며, 3개월 이후 비응답자는 응답자에 비해 더 높은 SWE를 보였다(3개월 × 비응답자 추정치 +4.20 kPa, P < 0.001; 12개월 × 비응답자 추정치 +9.41 kPa, P < 0.001)(표 3). 변형률 모델 역시 동일한 패턴을 보였으며, 유의미한 후기 발산(12개월 추정치 +0.68, P < 0.001× 비응답자)이 나타났다(표 3).

상관분석 결과, 초기 절제 후 탄성학 변화(즉각적 마이너스 베이스라인)는 12개월 체적 반응과 양의 상관관계가 있었습니다(그림 4). δ스트레인 비율은 12개월 VRR과 중간 정도의 상관관계를 보였으며(Pearson r = 0.445, P < 0.001), 반면 ΔSWE는 약하지만 통계적으로 유의미한 상관관계를 보였습니다(Pearson r = 0.302, P < 0.001). 이 결과는 더 큰 급성 절제 후 강성 변화가 이후 부피 감소와 연관되어 있어, 탄성학이 장기 치료 반응에 유용한 초기 생체역학적 효과를 포착한다는 가설을 뒷받침합니다.

기초 결절 부피, 조성, 전달 에너지, 연령을 보정한 다변수 회귀에서, 초기 절제 후 탄성학 변화는 독립적으로 12개월 VRR을 예측했습니다(표 4). 초기 변형률 변화(즉각적으로 기초선 손실)는 12개월 VRR 증가와 강하게 연관되어 있었으며(β = +1.0 14.80%, P < 0.001), 초기 SWE 변화도 독립적으로 기여했습니다(β = +5 kPa당 1.80%, P = 0.011). 기저 부피가 클수록 12개월 VRR이 낮아졌고(β = −1.05%/+1 mL, P < 0.001), 고체 조성과 더 높은 전달 에너지는 더 높은 VRR과 관련이 있었으며( 두 < 모두 0.05) (표 4).

추적 관찰 중, 반복 영상, 세포학 또는 임상 평가를 통해 악성 또는 대체 병리로 재분류된 치료된 결절은 없었습니다. 일부 반응하지 않는 결절은 지속적인 잔류 부피나 감소된 부피 감소를 보였으나, 12개월 추적 기간 동안 반복 소작술이나 수술이 필요한 명백한 재성장은 관찰되지 않았습니다.

데이터 가용성:
본 연구에서 제시된 분석의 원시 데이터는 투명성과 재현성을 높이기 위해 보충 표 1 로 준비 및 업로드되었습니다. 데이터셋에는 분석에 사용된 식별 해제된 참가자 수준 인구통계학, 임상, 초음파, 탄성학, 절차적 및 종단적 추적 관찰 변수가 포함되어 있습니다.

그림 1: 연구 설계 및 분석 워크플로우. 연구 기간 동안 선별된 참가자들은 적격성 평가를 받고 사전 정해진 제외 기준 후 포함되었습니다. 참가자당 1개의 갑상선 지수 결절을 분석하였다. 12개월 시점의 반응 상태는 볼륨 감소 비율(VRR)으로 정의되었으며, 반응자(VRR ≥ 50%)와 비반응자(VRR < 50%)로 정의되었습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

그림 2: 고주파 소작 후 종방향 전단파 탄성도 강성. 평균 전단파 탄성도(SWE) 강성도(kPa)는 전체 코호트에서 기저선, 절제 직후, 1, 3, 6, 12개월 시점에 나타나며, 12개월 반응 상태별로 층층화됩니다. 오차 바는 표준편차를 나타냅니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

그림 3: 결절 부피 및 부피 감소 비율의 종단적 변화. (A) 기저라인, 절제 직후, 1, 3, 6, 12개월 시점의 평균 부피 감소 비율(VRR, %)을 12개월 반응 상태별로 층별로 나누어 나누었습니다. 오차 바는 표준편차를 나타냅니다. 괄호는 각 시점에서의 그룹 간 비교를 나타냅니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

그림 4: 초기 탄성학 변화와 12개월 체적 반응 간의 관계. 산점도는 초기 절제 후 탄성학 변화(즉각적 기저선 감소)와 12개월 VRR 간의 연관성을 보여줍니다: (A) Δstrain ratio와 (B) ΔSWE 강성도. 그룹 간 시각적 비교를 용이하게 하기 위해 산란점 위치는 원래 전체 분석에서 보존되었습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

섹션	가변	응답자 (n = 146)	비응답자 (n = 41)	검정 통계량 (t/χ2)	P 값
인구 통계 및 기초 임상 상태
	나이, 연도	45.94 ± 7.88	49.63 ± 8.01	-2.640	0.009
	성	여성: 119명 (81.5%); 남성: 27명 (18.5%)	여성: 29명 (70.7%); 남성: 12명 (29.3%)	2.252	0.133
	BMI, kg/m2	23.74 ± 1.96	24.74 ± 1.94	-1.914	0.057
	흡연 (현재/과거)	No: 126명 (86.3%); 예: 20명 (13.7%)	번호: 32명 (78.0%); 찬성: 9명 (22.0%)	1.664	0.197
	고혈압	번호: 119명 (81.5%); 찬성: 27명 (18.5%)	번호: 30명 (73.2%); 찬성: 11석 (26.8%)	1.374	0.241
	당뇨병	No: 136명 (93.2%); 예: 10명 (6.8%)	No: 37명(90.2%); 찬성: 4명 (9.8%)	0.391	0.532
	관상동맥 질환	번호: 140명 (95.9%); 찬성: 6석 (4.1%)	No: 38명(92.7%); 예: 3명 (7.3%)	0.719	0.397
	만성 신장 질환	No: 143명 (97.9%); 예: 3석 (2.1%)	No: 40명(97.6%); 예: 1명 (2.4%)	0.023	0.881
	COPD/천식	번호: 140명 (95.9%); 찬성: 6석 (4.1%)	No: 39명(95.1%); 찬성: 2명 (4.9%)	0.046	0.830
	압박 증상 점수, 0–10	3.19 ± 0.63	3.43 ± 0.57	-2.199	0.029
	화장품 문제 점수, 1–4	2.79 ± 0.49	3.01 ± 0.47	-2.562	0.011
갑상선 검사와 자가면역
	TSH, mIU/L	1.74 ± 0.33	1.86 ± 0.31	-2.084	0.039
	자유 T4, pmol/L	15.06 ± 2.03	14.95 ± 1.92	0.310	0.757
	TPOAb 양성	번호: 129명 (88.4%); 찬성: 17명 (11.6%)	번호: 35명 (85.4%); 찬성: 6명 (14.6%)	0.265	0.606
	TgAb 양성	번호: 132명 (90.4%); 찬성: 14명 (9.6%)	번호: 36명 (87.8%); 찬성: 5명 (12.2%)	0.238	0.626
	레보티록신 사용	번호: 134명 (91.8%); 찬성: 12명 (8.2%)	번호: 36명 (87.8%); 찬성: 5명 (12.2%)	0.612	0.434
	항혈소판/항응고제 사용	번호: 139명 (95.2%); 찬성: 7석 (4.8%)	No: 38명(92.7%); 예: 3명 (7.3%)	0.402	0.526
기본 초음파 특징
	결절 위치	좌측: 66명 (45.2%); 지협: 6개 (4.1%); 우측: 74명 (50.7%)	왼쪽: 19명 (46.3%); 지협: 2개 (4.9%); 우측: 20명 (48.8%)	0.076	0.963
	최대 직경, mm	29.31 ± 3.94	33.87 ± 4.03	-6.516	<0.001
	기준 결절 부피, mL	8.33 ± 1.90	12.04 ± 2.41	-10.385	<0.001
	조성 (고체)	번호: 44명 (30.1%); 찬성: 102명 (69.9%)	No: 17명 (41.5%); 찬성: 24명 (58.5%)	1.868	0.172
	낭성 구성 요소	<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
	석회화	굵기: 10명 (6.8%); 마이크로: 35 (24.0%); 없음: 101명 (69.2%)	거칠게: 3마리(7.3%); 마이크로: 8명 (19.5%); 없음: 30명 (73.2%)	0.360	0.835
	기저혈관 등급	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
기저 탄성학
	기본 SWE 강성, kPa	31.50 ± 5.11	35.08 ± 5.34	-3.925	<0.001
	기준 변형률	2.41 ± 0.37	2.77 ± 0.40	-5.407	<0.001
절차적 매개변수 및 안전성
	시술 시간 1, 분	18.42 ± 2.79	21.13 ± 3.01	-5.401	<0.001
	소작력 W	35.02 ± 4.78	37.05 ± 4.63	-2.419	0.017
	전달 에너지, kJ	11.04 ± 1.98	13.18 ± 2.63	-5.665	<0.001
	수압해부 사용	No: 100명 (68.5%); 찬성: 46명 (31.5%)	번호: 25명 (61.0%); 찬성: 16명 (39.0%)	0.816	0.366
	진정제 사용	번호: 115명 (78.8%); 찬성: 31명 (21.2%)	번호: 30명 (73.2%); 찬성: 11석 (26.8%)	0.576	0.448
	시술 중 통증, VAS 0–10	2.60 ± 0.59	2.78 ± 0.61	-1.713	0.088
	경미한 출혈/혈종	번호: 139명 (95.2%); 찬성: 7석 (4.8%)	No: 37명(90.2%); 찬성: 4명 (9.8%)	1.423	0.233
	일시적인 음성 변화	No: 143명 (97.9%); 예: 3석 (2.1%)	번호: 41명 (100.0%); 예: 0 (0.0%)	0.856	0.355

표 1: 12개월 부피 반응에 따른 기초 인구통계학, 임상, 실험실, 초음파, 탄성학 및 절차적 특성(VRR ≥50% 대 <50%)에 따라 계층화됨. 데이터는 평균 ± SD(또는 중앙값 [IQR])와 n(%)으로 표시됩니다. 그룹 간 비교에서는 연속 변수에 대해 학생 t 검정 또는 윌콕슨 순위합 검정을, 범주 변수에 대해서는 χ² 또는 피셔 정확 검정을 사용했습니다. 양측 P < 0.05는 통계적 유의성을 나타냅니다. 약어: BMI, 체질량지수; 갑상선 자극 호르몬인 TSH; FT4, 자유 갑상선; TPOAb, 갑상선 과산화효소 항체; TgAb, 갑상선글로불린 항체; SWE, 전단파 탄성학; ROI, 관심 지역; VAS, 시각적 아날로그 스케일; RFA, 고주파 소작; VRR, 볼륨 감소비.

방문	그룹	n	부피(mL), 평균 ± SD	VRR (%), 평균 ± SD	SWE (kPa), 평균 ± SD	변형률, 평균 ± SD
기본 기준	대응자	146	8.33 ± 1.90	—	31.5 ± 5.1	2.41 ± 0.37
	응답하지 않는 사람	41	12.04 ± 2.41	—	35.1 ± 5.3	2.77 ± 0.40
즉각	대응자	146	8.75 ± 1.99	—	41.7 ± 5.6	3.14 ± 0.41
	응답하지 않는 사람	41	12.66 ± 2.54	—	43.4 ± 5.8	3.35 ± 0.42
1개월	대응자	146	5.69 ± 1.35	31.5 ± 6.3	33.6 ± 5.0	2.51 ± 0.39
	응답하지 않는 사람	41	10.08 ± 2.36	16.1 ± 3.8	36.1 ± 5.3	2.82 ± 0.41
3개월	대응자	146	3.72 ± 0.85	55.3 ± 7.1	26.4 ± 4.8	2.05 ± 0.36
	응답하지 않는 사람	41	8.70 ± 2.10	28.1 ± 6.8	32.6 ± 5.1	2.62 ± 0.39
6개월	대응자	146	2.99 ± 0.70	64.1 ± 7.0	23.8 ± 4.7	1.91 ± 0.35
	응답하지 않는 사람	41	7.93 ± 1.95	34.2 ± 8.1	31.6 ± 5.3	2.56 ± 0.40
12개월	대응자	146	2.20 ± 0.45	73.6 ± 6.9	21.6 ± 4.7	1.87 ± 0.35
	응답하지 않는 사람	41	7.35 ± 1.80	38.5 ± 8.7	31.2 ± 5.2	2.55 ± 0.41

표 2: 고주파 소작술 후 12개월까지의 결절 부피, 부피 감소 비율 및 탄성학 지표의 종단 변화를 치료 반응에 따라 층별로 나누어 나타났습니다. 각 방문마다 가치 사항이 요약되며; n은 사용 가능한 관측을 나타냅니다. VRR은 기준 부피를 기준으로 계산되었습니다. 종단 추론은 혼합 효과 모델을 사용하여 수행되었습니다. 약어: VRR, 부피 감소 비율; SWE, 전단파 탄성학; RFA, 고주파 소작술.

결과	고정 효과 항	추정치	SE	t	P 값
SWE 강성 (kPa)
	(가로채기) 기본 상태, 응답자	31.51	0.55	57.273	<0.001
	그룹: 비응답자 (vs 응답자)	3.58	0.92	3.891	<0.001
	시간: 즉각적(기초 기준 대비)	10.22	0.62	16.452	<0.001
	시간: 1개월 (기준선 대비)	2.17	0.65	3.502	<0.001
	시간: 3개월 (기준 대비)	-5.14	0.61	-8.361	<0.001
	시간: 6개월 (기준 대비)	-7.73	0.64	-12.419	<0.001
	기간: 12개월 (기준선 기준 대)	-9.02	0.63	-14.286	<0.001
	그룹 × 시간: 즉시 × 응답하지 않음	-0.80	0.98	-0.816	0.415
	그룹 × 시간: 응답 없음 × 1개월	0.45	0.95	0.474	0.636
	그룹 × 시간: 3개월 × 응답하지 않는 경우	4.20	0.96	4.375	<0.001
	그룹 × 시간: 6개월 × 응답하지 않는 경우	6.10	0.98	6.224	<0.001
	그룹 × 시간: 12개월 × 응답하지 않는 경우	9.41	1.02	9.216	<0.001
변형률
	(가로채기) 기본 상태, 응답자	2.41	0.04	60.252	<0.001
	그룹: 비응답자 (vs 응답자)	0.36	0.07	5.143	<0.001
	시간: 즉각적(기초 기준 대비)	0.73	0.05	14.601	<0.001
	시간: 1개월 (기준선 대비)	0.10	0.07	2.012	0.046
	시간: 3개월 (기준 대비)	-0.36	0.08	-7.211	<0.001
	시간: 6개월 (기준 대비)	-0.51	0.06	-10.032	<0.001
	기간: 12개월 (기준선 기준 대)	-0.54	0.05	-10.802	<0.001
	그룹 × 시간: 즉시 × 응답하지 않음	-0.08	0.08	-1.000	0.318
	그룹 × 시간: 응답 없음 × 1개월	0.05	0.09	0.625	0.532
	그룹 × 시간: 3개월 × 응답하지 않는 경우	0.32	0.04	4.031	<0.001
	그룹 × 시간: 6개월 × 응답하지 않는 경우	0.55	0.08	6.875	<0.001
	그룹 × 시간: 12개월 × 응답하지 않는 경우	0.68	0.07	8.521	<0.001

표 3: 추적 방문 간 종방향 탄성학 궤적에 대한 선형 혼합 효과 모델 결과, 참가자 수준 무작위 인터셉트와 시간 및 반응별 상호작용에 대한 고정 효과 포함. 선형 혼합 효과 모델에는 참가자 수준의 무작위 절편과 범주별 시간 효과가 포함되었으며; 반응 그룹과 시간별 상호작용 항을 사용하여 궤적을 비교하였다. 모델 적합도는 잔여 진단을 사용하여 평가되었으며; 필요할 때 로그 변환이 적용되었습니다. 추정치는 95% 신뢰감시구간(CI)으로 보고됩니다. 약어: SWE, 전단파 탄성학; 신뢰구간(CI), 신뢰구간; VRR, 볼륨 감소비.

결과	용어	베타	SE	t	DF	P 값	신뢰감시구간(95%)
12개월 VRR (%)
	(가로채기)	82.00	4.20	19.524	180	<0.001	73.71에서 90.29
	Δ 변형률 (즉각형 - 기준선), 1.0 기준	14.80	2.40	6.167	180	<0.001	10.06부터 19.54까지
	Δ SWE(즉각 - 기준선), 5 kPa당	1.80	0.70	2.571	180	0.011	0.42에서 3.18
	기준 결절 부피, 1 mL당	-1.05	0.28	-3.752	180	<0.001	-1.60에서 -0.50
	탄탄한 작문 (예 vs 노)	3.20	1.40	2.286	180	0.023	0.44에서 5.96
	1 kJ당 전달된 에너지	0.55	0.22	2.512	180	0.013	0.12에서 0.98
	10년당 나이	-1.10	0.70	-1.571	180	0.118	-2.48에서 0.28

표 4: 초기 절제 후 탄성학 변화가 기저선 결절 및 절차적 공변량 조정 후 12개월 체적 반응을 독립적으로 예측하는지 평가하는 다변수 회귀분석. 다변수 선형 회귀는 조기 탄성도 변화(즉각 마이너스 베이스라인)와 공변량 조정 후 12개월 VRR 간의 연관성을 평가하였습니다. SE, t, df, 양측 P, 95% CI를 가진 계수(β)가 보고됩니다; DF는 잔여 자유도를 나타냅니다. 약어: VRR, 부피 감소 비율; SWE, 전단파 탄성학; SE, 표준오차; 신뢰구간(CI), 신뢰구간; DF, 자유도.

보충 파일 1: 이 연구의 원시 데이터.이 파일을 다운로드하려면 여기를 클릭해 주세요.

Discussion

본 연구는 갑상선 결절의 RFA 후 치료 반응을 평가하기 위한 탄성학 기반 접근법을 평가하였으며, 절제 후 강성 변화와 관련된 정량적 영상 바이오마커로 엘라스그래피를 활용하였습니다. 핵심 제안은 SWE와 변형 탄성술로 측정되는 강성이 소작 효과의 조기 지표이자 소작 후 재형성의 종단적 지표로서 기존 부피 기반 평가단과를 보완할 수 있다는 것입니다. 연구 결과는 두 가지 수준에서 이 개념을 지지합니다. 첫째, 소작술 직후 SWE 강성과 변형률 모두 증가했으며, 이는 급성 열응고와 초기 조직 반응과 일치하는 패턴이다. 둘째, 추적 관찰 과정에서 강직성이 감소했으며, 반응자가 비응답자보다 늦게 더 큰 감소를 보였으며, 이는 부피 결과와 VRR에서 관찰된 분리와 유사합니다. 또한, 초기 절삭술 후 탄성소매 변화는 12개월 VRR과 통계적으로 유의한 연관성을 보였으며, 기저선 결절 및 절차적 공변량 조정 후에도 독립적으로 연관성을 유지하였습니다. 이 결과들은 탄성학이 RFA 치료 반응과 관련된 종방향 강성 변화를 포착하며, 초기 강성 변화가 기준선 크기, 조성, 에너지 전달만으로는 완전히 설명되지 않는 예후 정보를 담고 있다는 가설을 뒷받침합니다.^{12, 13, 14, 15, 16}.

관찰된 궤적은 이전에 기술된 소작술 후 조직 변화와 일치하며, 12,13,14 소작 후 영상에 대한 임상적으로 유용한 서사를 제공합니다. 즉각적인 강성 상승은 소작부위 내 응고성 괴사와 부종, 초기 염증 변화를 반영하는 것으로 보이며, 이들은 변형 저항성을 증가시키고 파동 전파를 변화시킵니다^15,16. 시간이 지남에 따라 뻣뻣함 감소는 괴사 조직의 점진적 흡수, 절제부의 조직화, 그리고 궁극적으로 수축과 증상 개선을 가져오는 재형성과 일치한다 15,17,18,19. 반응자들이 늦게 경직이 더 크게 감소한 사실은 성공적인 절제술이 단순한 급성 사건이 아니라 지속적인 재형성 과정임을 시사하며, 불완전한 치료나 잔류 생존 조직이 둔해진 부화 궤적으로 표현될 수 있음을 시사한다²⁰. 이 관점은 초기 부피 평가의 공통된 한계를 조화시키는 데 도움을 줍니다. RFA 직후에는 부종, 출혈, 반응성 변화로 인해 부피가 기준선 근처에 머무르거나 일시적으로 변동할 수 있어, 초기 부피 측정은 절제 구역이 적절한지 여부를 완전히 판단하지 못합니다. 반면 탄성학은 크기 기반 지표가 가장 구^{별력이 낮은} 시기에 소작술의 주요 물리적 결과에 반응합니다. 혼합 효과 모델에서 정량화된 후기 강성 궤적의 발산은 강성이 지속 가능한 부피 감소와 일치하는 조직 진화의 의미 있는 차이를 포착한다는 개념을 더욱 뒷받침합니다.

과학적 발전의 관점에서 본 연구는 탄성학을 주로 진단적 역할에서 치료적 모니터링 패러다임으로 전환시키는 데 기여합니다. 양성 갑상선 결절 치료에서는 영상 유도 절제술이 점점 더 많이 채택되고 있지만, 추적 치료 전략은 여전히 6개월에서 12개월의 VRR과 같은 지연된 부피 결과에 크게 의존하고 있습니다. 치료 효과를 객관적으로 조기 판독하는 바이오마커는 시술 후 초기 기간 동안 추가 영상 정보를 제공하고 위험 적응 추적 전략에 대한 추가 연구를 지원할 수 있습니다^22,23. 유사한 탄성학 기반 접근법은 다른 열 소작 환경에서도 탐구되었습니다²⁴. 강성 매핑은 비침습적이고, 일상적인 초음파 워크플로우에 통합할 수 있으며, 조영제^{나 방사}선 없이 반복할 수 있어 매력적입니다. 워크플로우 관점에서 탄성학술은 표준화된 획득 프로토콜을 사용할 경우 검사 시간을 크게 늘리지 않고도 일상적인 소작술 후 초음파 검사에 통합할 수 있습니다. 그러나 더 넓은 구현에는 플랫폼 간 변동성, 운영자 경험, 그리고 서로 다른 초음파 시스템 간 획득 매개변수의 조화에 주의를 기울여야 합니다. 본 연구에서 이러한 변동성 원천을 줄이기 위해, 모든 탄성촬영 측정은 추적 관찰 기간 동안 동일한 초음파 플랫폼과 소프트웨어 환경을 사용하여 수행되었습니다. 임상적 사용성은 절대 강성 값보다는 안정적인 획득 조건에서 얻는 일관된 종방향 변화 패턴에 더 의존할 수 있습니다. 따라서 탄성학이 일상적인 소작 후 감시 경로에 보다 널리 통합되기 위해서는 표준화된 운영자 교육, 사전 정의된 품질 관리 규칙, 프로토콜 조화가 필수적입니다. 현재 결과는 표준화된 탄성학 획득 프로토콜이 임상 연구에 적합한 보다 일관된 종방향 강성 패턴을 생성하는 데 도움이 될 수 있음을 시사합니다. 중요한 점은, 본 접근법은 결절 내 기준 값을 내부 통제로 간주하여, 절대적 강직성만보다 개인별로 더 견고할 수 있는 변화 지표를 강조한다는 것이다. 추가로 고려할 점은 현재의 탄성학 프로토콜이 주로 처리된 결절의 고체 구성 요소 내 강성 측정에 초점을 맞추었으며, 주위 조직은 체계적으로 평가되지 않았다는 점입니다. 이는 말초 치료 부위가 염증 반응, 열 확산, 진화하는 섬유증, 잔류 생존 조직, 소작술 후 초기 조직 재형성과 관련된 생물학적으로 의미 있는 정보를 포함할 수 있기 때문에 잠재적으로 중요한 한계를 나타냅니다. 신흥 영상 문헌은 치료 반응 평가에서 병변 말초와 미세환경 특성의 진단 및 생물학적 중요성을 점점 더 강조하고 있습니다. 따라서 결절 내 및 주위 강성 매핑을 결합한 향후 연구는 절제 후 조직 진화에 대한 보다 포괄적인 생체역학적 특성화를 제공하며, 갑상선 RFA 후 초기 반응 평가를 더욱 향상시킬 수 있습니다.

실무 구현 시 엘라스트그래피 획득 및 소작 후 평가 시 여러 기술적 함정을 고려해야 합니다. 과도한 트랜스듀서 압축, 삼킴 동작, 불안정한 호흡 정지, 이질성 낭포성 퇴행, 불완전한 색채 채움은 측정 불안정성을 유발하고 재현성을 저하시킬 수 있습니다. 소작술 후 즉각적인 평가는 부종, 일시적 출혈, 반응성 염증 변화에 의해 영향을 받을 수도 있습니다. 프로브 위치 유지, 외부 압축 최소화, 미리 정의된 관심 영역 배치 규칙 사용, 프레임 간 변동이나 신호 공백이 표시될 때 측정값을 재획득하는 것은 획득 신뢰성을 높이기 위한 중요한 문제 해결 전략입니다.

재현 가능한 종단 탄성학 평가에 특히 중요한 절차적 변수들이 있습니다. 여기에는 일관된 트랜스듀서 압력, 안정적인 목 위치, 표준화된 획득 타이밍, 심하게 석회화된 또는 낭포성 부위 회피, 그리고 연속 검사 간 관심 영역 깊이 유지가 포함됩니다. 재현성은 또한 동일한 엘라스그래피 프리셋을 유지하고, 프로토콜별 교육과 사전 정의된 품질 관리 기준을 통해 운영자 간 변동을 최소화하는 것에 의존할 수 있습니다.

이 발견에서 여러 구체적인 연구 응용 분야가 도출됩니다. 첫째, 탄성학은 초기 생체역학적 반응이 최적이 아닌 환자를 후속 반응 경로에 따라 계층화하는 데 도움을 줄 수 있으며, 이들은 더 면밀한 감시, 조기 보조 영상 촬영, 또는 추가 절제술 고려가 도움이 될 수 있습니다. 따라서 단일 프로토콜이나 플랫폼에서 도출된 강성 임계값은 더 넓은 다중센터 표준화가 이루어질 때까지 신중하게 해석되어야 합니다. 둘째, 강성 궤적은 절제 기법, 전극 전략 또는 보조 절차 방법의 임상시험에서 중간 평가 지표로 작용할 수 있어, 치료 간 의미 있는 차이를 감지하는 데 소요되는 시간을 단축할 수 있습니다. 셋째, 탄성학은 부피 감소가 지연되었지만 증상이 변하거나, 부피 감소는 일어나지만 잔여 강성이 불완전한 재형성을 시사하는 모호한 사례를 해석하는 데 도움을 줄 수 있습니다. 넷째, 엘라스그래피를 예측 모델에 통합하면 초기 변화 지표와 기초 부피, 조성, 전달 에너지 등 확립된 결정 요인을 결합하여 개별화된 예측을 향상시킬 수 있습니다. 본 연구의 다변수 결과는 초기 변형률 변화와 SWE 변화가 기존 공변량을 넘어 독립적인 정보를 제공했기 때문에 이러한 통합의 실현 가능성을 뒷받침합니다. 이러한 결과는 향후 예측 또는 위험 층화 모델에 탄성학 지표를 통합하는 잠재적 가능성을 시사합니다.

이 가설을 연구하는 다른 방법들도 있으며, 이들도 고려할 가치가 있다. 조영제 강화 초음파 및 혈관 평가를 포함한 관류 기반 영상은 잔류 생존 조직을 식별할 수 있으며, 효과적인 치료 분획을 추정하여 장기적인 반응을 예측할 수 있습니다²⁶. 초기 소작 비율과 초기 소작 부피와 같은 형태학적 지표도 장기적인 수축 및 재치료 위험 예측에 사용되었습니다^27,28. 단면 영상은 선택된 경우에 대해 소작 부위와 합병증을 추가로 특성화할 수 있습니다. 환자 보고한 결과, 증상 완화, 미용 개선은 임상적 이점의 보다 직접적인 지표를 제공할 수 있으므로, 영상 바이오마커와 함께 포함되어야 합니다. 이 맥락에서 탄성학은 배타적이기보다는 보완적인 것으로 봐야 합니다. 강력한 미래 방향은 동일한 코호트 내에서 관류 기반 방법 및 형태학적 대체체를 엘라스토그래피와 직접 비교하는 것이며, 이는 조화로운 추적 관찰 일정과 표준화된 결과 정의를 사용하는 것입니다. 본 연구는 초기 탄성학 변화와 이후 VRR 사이에 유의미한 연관성을 보여주었으나, 본 코호트는 검증된 민감도 및 특이도 추정치를 포함한 임상적으로 실행 가능한 컷오프 값을 도출하도록 특별히 설계되거나 검증된 기준이 아니었습니다. 따라서 사전 정의된 결과 기반 임계값 분석을 가진 전향적 다기관 연구는 향후 임상 번역에 중요할 것입니다. 마지막으로, 조기 탄성학 중심 관리 결정을 표준 추적 관찰 전략과 비교하여 VRR, 증상 및 미용 개선, 재치료율 등의 결과를 반영하는 전향적 연구가 이 분야에 도움이 될 것입니다.

이 연구의 한계는 연구 결과를 해석하고 추가 연구를 계획할 때 반드시 고려해야 합니다. 첫째, 엘라스그래피는 변환기 압력, 깊이, 관심 영역 위치, 장치별 재구성 등 획득 조건에 본질적으로 민감하며, 이러한 요소들은 프로토콜화에도 불구하고 측정 변동성을 초래할 수 있습니다. 확산성 갑상선 뻣뻣함 이상으로 인한 잠재적 혼란을 줄이기 위해, 임상적으로 유의한 확산성 갑상선 질환과 확산성 이질성 실질성 침범이 있는 환자는 등록 시 제외하였습니다. 본 연구에서는 공식적인 상호관찰자 재현성 분석이 수행되지 않았으나, 표준화된 획득 프로토콜과 조작자 교육 절차가 적용되어 측정 일관성을 개선하였습니다. 둘째, 추적 관찰 일정은 주요 이정표를 포착하지만, 부종과 반응성 변화가 경직에 영향을 미치고 급성 치료 효과와 일시적인 시술 후 생리학적 분리를 복잡하게 만들 수 있는 아주 초기 단계의 변화를 완전히 특징짓지는 못할 수 있습니다. 셋째, VRR 임계값을 기반으로 한 반응 분류는 널리 사용되고 있지만, 증상 개선, 미용 결과, 장기 재성장이 VRR을 완벽하게 추적하지 못할 수 있어 임상 성공의 모든 차원을 포함하지는 못합니다. 또한, 현재의 추적 관찰 기간은 12개월로 제한되어 있었으나, 치료된 갑상선 결절은 더 긴 관찰 기간에 걸쳐 재형성, 축소 또는 재생 지연이 계속될 수 있습니다. 현재 추적 관찰 기간 동안 반복 절제술이나 수술이 필요한 명백한 재성장은 관찰되지 않았으나, 장기 종단 추적 관찰이 가능한 향후 연구에서는 재성장 패턴, 재치료율, 반복 절제술 요구사항을 추가적인 임상적 결과 측정으로 포함해야 합니다. 넷째, 다변수 모델은 주요 기준선 및 절차적 공변량을 조정하지만, 여진 완전성, 전극 위치 역학, 조작자 경험과 같은 측정되지 않은 기술적 뉘앙스를 포함한 잔여 교란 요소는 초기 강성 변화와 장기 반응 모두에 영향을 미칠 수 있습니다. 또한, 모든 절차가 경험 많은 갑상선 RFA 시술사가 단일 센터에서 수행했기 때문에, 저부피 환경으로의 외부 일반화는 추가 평가가 필요합니다. 다섯째, 연구에서는 관류 기반 영상이나 공식 절제 구역 생존 지표와 같은 대체 초기 반응 도구와의 직접적인 비교를 포함하지 않아, 엘라스그래피가 이들 방법에 비해 추가적인 가치를 완전히 확립할 수 없습니다. 연구가 관찰 및 영상 기반이었기 때문에, 확인된 연관성은 탄성학 변화와 근본적인 조직 생물학 간의 직접적인 기전적 인과관계 증거로 해석되어서는 안 됩니다. 향후 연구에서는 기저선 엘라스그래피, 즉각적인 절제 후 강성 측정, 1개월 추적 엘라스그래피의 독립적인 예측 성능과 기저선과 RFA 직후 값 간의 상대적 비율 변화를 추가로 평가할 수 있습니다. 이러한 분석은 절대 강성 측정이나 동적 강성 변화 지표가 갑상선 RFA 후 장기 체적 반응 평가에 더 큰 예후 유용성을 제공하는지 명확히 하는 데 도움이 될 수 있습니다.

탄성촬영 기반 평가는 갑상선 결절 고주파 소작술 후 종단 평가를 위한 실행 가능하고 프로토콜 표준화된 접근법을 제공합니다. 초기 절제 후 강성 변화는 12개월 부피 반응과 독립적으로 연관되어 있으며, 상당한 부피 감소가 나타나기 전에 임상적으로 유용한 정보를 제공할 수 있습니다. 이러한 발견은 표준화된 탄성촬영 프로토콜이 갑상선 RFA 후 절제 반응 평가 및 종단적 추적 관찰을 위한 보완적인 영상 바이오마커로서 역할을 할 가능성을 뒷받침합니다.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
3M Tegaderm Film Dressing	Tegaderm	https://www.3m.com/3M/en_US/p/c/medical/bandages-dressings/film/i/health-care/medical/	접착식 드레싱; 수술 후 부위 덮개.
일회용 주사기	Becton Dickinson	BD 309628	하이드로디섹션 및 국소 마취제 투여
RFA용 접지 패드	3M Healthcare	모델 9547	RF와 호환되는 일회용 반환 전극.
고주파 선형 어레이 초음파 변환기 (5–14 MHz)	SuperSonic Imagine	SL15-4 선형 프로브	갑상선 초음파 및 탄성도 획득
고주파 선형 어레이 초음파 변환기 (7–15 MHz)	Mindray	L11-3U	갑상선 이미징 및 탄성도에 최적화된 선형 프로브.
IntelliVue	Philips Healthcare	MX450	환자 모니터링 시스템; RFA 중 ECG, SpO2, 비침습성 혈압 모니터
IntelliVue	Philips Healthcare	MP50	생체 신호 모니터; 시술 중 지속적 모니터링
리도카인 하이드로클로라이드 주사	Hospira Inc.	NDC 0409-4276-17	절제 전 국소 마취
lme4 패키지	CRAN	RRID:SCR_015654	R에서 선형 혼합 효과 모델링.
lmerTest 패키지	CRAN	RRID:SCR_015656	R에서 혼합 모델 검정 및 p-값.
Microflex 분말 무료 장갑	Ansell Healthcare		절차 팀을 위한 무균 장벽.
비 멸균 초음파 젤	EcoGel		기준 및 후속 이미징을 위한 커플링 매체.
Picture Archiving and Communication System (PACS)	GE Healthcare	Centricity PACS	엔터프라이즈급 이미지 저장 및 검색 시스템.
포비돈-요오드 10%	Betadine		시술 전 살균제.
Precision Workstation	Dell	T5820	이미지 검토 및 데이터 분석용.
PrecisionGlide 22G–27G	BD (Becton Dickinson)	https://www.bd.com/en-us/products-and-solutions/products	국소 마취/하이드로디섹션 용 바늘 세트.
사전 혼합된 식염수 백 500 mL (0.9% NaCl)	Baxter Healthcare		하이드로디섹션 및 세척용 유체.
R 통계 컴퓨팅 소프트웨어	R Foundation	RRID:SCR_001905	통계 분석 환경.
RFA 경피 냉각 전극	STARmed Co., Ltd.	Star RF Electrode 18-10s10F	갑상선 절제용 활성 RFA 전극 (TaeWoong Medical USA)
RStudio 통합 개발 환경	Posit (구 RStudio)	RRID:SCR_000432	R 통계 분석을 위한 IDE.
전단파 탄성도 소프트웨어 모듈	Mindray	SWE 모듈 (RESONA 시리즈)	US 시스템 내 kPa의 정량적 강성 매핑.
멸균 드레이프 키트	Medline Industries	DYNJP2500	갑상선 RFA 시술을 위한 멸균 준비
멸균 거즈 및 면봉	Covidien	멸균 스펀지 키트	피부 준비 및 지혈에 사용.
멸균 생리 식염수	Baxter Healthcare	2B1324X	갑상선 RFA 중 하이드로디섹션
멸균 절차 드레이프 (블루 드레이프 팩)	Medline Industries		멸균 경피 RFA 설정을 위한 표준 드레이프.
멸균 초음파 젤	Aquasonic 100	전도성 젤 100 mL	US 안내용 멸균 커플링 매체.
멸균 초음파 프로브 커버	Parker Laboratories	모델 1001	절차 중 초음파 획득용 멸균 커버.
변형 탄성도 소프트웨어 모듈	Mindray	SE 모듈 (RESONA 시리즈)	통합 품질 지표를 갖춘 변형 비 계산.
주사기 (5/10/20 mL)	BD (Becton Dickinson)	Monoject 주사기	마취제 및 하이드로디섹션 투여용.
삼방향 스톱콕	BD (Becton Dickinson)	모델 305129	제어된 유체 관리를 위한 스톱콕 및 튜브.
갑상선 RF 전극	STARmed Co., Ltd.	https://en.starmed4u.com/	내부 냉각 RF 전극