텔로미어 길이와 갑상선 독성증의 연관성: 두 표본 멘델 무작위화 연구에서 얻은 통찰

요약

이 프로토콜은 공개 유전 요약 데이터를 사용하여 텔로미어 길이가 갑상선독성증 위험에 인과적으로 영향을 미치는지 평가하는 두 표본 멘델 무작위화 파이프라인을 설명합니다. 기기 선택, 조화, 1차 추정, 민감도 분석, 재현 가능한 R 코드를 포함하며, 투명한 보고를 지원하기 위한 수치 출력을 제공합니다.

초록

갑상선독성증은 갑상선 호르몬이 과도하게 분비되는 내분비 질환이지만, 전신 노화 생물학과의 원인학적 연관성은 아직 완전히 정의되지 않았습니다. 텔로미어 길이(TL)는 세포 노화와 유전체 안정성을 반영하며, 여러 복잡한 질병과 관련이 있는 것으로 알려져 있습니다. 우리는 유전적으로 예측된 TL이 갑상선독성증 위험에 미치는 인과관계를 평가하기 위해 두 표본 멘델 무작위화(MR) 연구를 수행하였습니다. TL을 위한 유전 도구는 유럽 혈통에 대한 대규모 유전체 전체 연관 연구(GWAS)에서 도출되었습니다(n > 470,000명). 갑상선독성증 요약 통계는 최신 FinnGen 발표(≈4,000건, 대조군 >210,000건)에서 얻었습니다. 1차 역분산 가중 분석에서는 유전적으로 대체된 TL이 길수록 갑상선독성증 위험이 낮아지며, 효과의 방향과 정도는 보완 추정량(MR-Egger, 가중 중앙값/최대우도, MR-PRESSO, MR-RAPS)에서 일관됨을 보여주었습니다. 감도 분석 결과 방향성 다재성의 증거는 없었으며, 코크란의 Q를 사용해 이질성 평가를 수행하였다. Steiger 방향성 검사는 TL에서 갑상선독성증으로의 인과 흐름을 뒷받침하였다.

우리가 아는 한, 이 연구는 최신 GWAS 자료를 활용해 전체 갑상선독성증 위험과 TL 간의 인과관계를 평가한 최초의 MR 분석 중 하나로, 갑상선기능항진증 관련 표현형에 초점을 맞춘 기존 증거를 확장한 것입니다. 이러한 발견은 TL이 지표하는 세포 노화 과정이 갑상선독성증 감수성에 기여할 수 있음을 시사하며, 갑상선 기능 장애의 텔로미어 생물학에 관한 향후 종단적 및 기전적 연구를 촉진합니다.

서론

갑상선독성증은 과도한 순환 갑상선 호르몬이 특징인 흔한 내분비 질환입니다. 대표적인 증상으로는 체중 감소, 빈맥, 과민성, 불안, 떨림이 있으며, 심한 경우 부정맥, 심부전,^{골 손실 소}인이 있습니다. 임상 특징과 치료 전략은 비교적 확립되어 있지만, 개인 간 감수성과 진행을 형성하는 생물학적 메커니즘은 아직 완전히 정의되지 않았습니다. 원인은 자가면역 과기능(예: 그레이브스병)부터 독성 결절, 갑상선염, 외인성 호르몬 노출, 약물 반응, 종양^{관련 원인에} 이르기까지 다양합니다. 이들 기관에서 면역 활성화, 산화 스트레스, 미토콘드리아 생합성 변화가 호르몬 과잉과 조직 손상 및 전신 합병증을 연결하는 그럴듯한 기전적 주제로 반복됩니다.

텔로미어 길이(TL)는 염색체 말단을 덮는 연쇄 반복 DNA-단백질 구조로, 세포 노화, 복제 역사, 유전체 안정성을 통합합니다. TL은 세포 분열과 함께 짧아지며, 산화 스트레스와 만성 염증에 의해 가속화되는데, 이 과정들은 갑상선 기능 장애에서도 두드러지게 나타납니다 ^3,4. 짧은 TL은 심장대사 장애, 특정 암, 신경퇴행^{성 질환과} 연관되어 있으며, 최근 연구된 데이터는 비만과 제2형 당뇨병과 같은 호르몬 및 면역-대사 축이 복제 스트레스와 교차하는 내분비 표현형과의 연관성을 시사합니다⁶. 갑상선 호르몬이 기저 대사율, 미토콘드리아 기능, 증식 신호를 조절하는 점을 고려할 때, TL 역학은 전신 과대사와 면역 조절 장애가 텔로미어 생물학에 영향을 미치거나 영향을 받을 수 있는 갑상선 독성 상태에서 특히 중요할 수 있습니다 ^7,8.

그러나 TL이 갑상선독성증에 인과적 역할을 하는지는 여전히 불확실하다. 기존의 관찰 연구는 교란(예: 생활습관, 동반 질환)과 역인과관계(TL에 영향을 미치는 질병)에 취약합니다. 이러한 한계를 해결하기 위해 우리는 TL과 강하게 연관된 생식세포 유전 변이체를 혼란 없는 도구로 활용하여^{질병 위험에} 미치는 영향을 검증하는 2표본 멘델 무작위화(MR)를 사용한다. 유전체 전체 연관 연구는 이제 충분히 강력한 TL 도구와 크고 표현형이 잘 정의된 결과 데이터셋을 제공하여 엄격한 인과 평가를 가능하게 합니다¹⁰. 유전학 역학과 텔로미어 생물학을 통합함으로써, 본 연구는 유전적으로 외형제가 갑상선독성증 위험에 영향을 미치는지 평가하고, 기전적 경로를 명확히 하며, 위험 계층화 및 향후 개입 개발을 위한 후보 생체지표로서 TL을 탐구합니다.

우리가 아는 한, 이 분석은 현대의 대규모 GWAS 자료를 활용해 TL과 관련된 갑상선독성증의 전반적인 위험을 조사한 최초의 MR 연구 중 하나로, 주로 갑상선 기능 항진증 관련 표현형이나 자가면역 갑상선 질환에 초점을 맞춘 연구를 확장한 것입니다^.

프로토콜

본 연구는 두 표본 멘델 무작위화(MR) 설계를 적용하여 유전적으로 반응하는 텔로미어 길이(TL)가 갑상선독성증 위험에 인과적으로 영향을 미치는지 평가합니다. 식별 해제된 공개 GWAS 요약 통계만을 사용했으며, 개인 수준의 데이터는 접근하지 않았습니다. 칭다오 시립병원의 기관심사위원회는 이 분석이 전적으로 공개된 요약 데이터에 의존하기 때문에 추가 검토에서 제외된다고 결정했습니다. 모든 기여하는 전유전체 연관 연구는 원래 프로토콜의 일환으로 사전 동의와 윤리적 승인을 받았습니다. 분석은 TL과 관련된 생식계 변이를 교란 없는 도구로 사용하고, 다변성, 이질성, 효과 방향을 조사하는 민감도 절차를 구현하여 교란 및 인과관계를 최소화하도록 구조화됩니다.

데이터 출처
TL의 노출 데이터셋은 IEU OpenGWAS에서 식별자 ieu-b-4879로 획득되었으며, 약 472,174명의 유럽 혈통 참가자를 포함하고 있습니다. 갑상선독성증의 결과 데이터셋은 FinnGen 컨소시엄에서 가져왔으며, 2021년 엔드포인트인 finn-b-thyROTOXICOSIS를 사용하여 4,142건의 사례와 213,693명의 대조군을 포함했습니다. 이 자료들은 SNP 수준 도구를 구성하고, 대립유전자를 조화시키며, 확립된 MR 추정량과 인과적 효과를 추정하는 데 필요한 효과적인 추정치와 표준오차를 제공합니다.

소프트웨어 및 컴퓨팅 환경
모든 분석은 R(버전 4.3.1)에서 수행되었습니다. 주요 분석 패키지는 TwoSampleMR(버전 0.5.7)으로, gwas 자원에 프로그래밍 접근을 위한 ieugwasr, 이상치 검출 및 왜곡 테스트를 위한 MRPRESSO, 약한 기기와 특이한 다면성에서 강건한 추정을 위한 mr.raps, 방사형 시각화를 위한 RadialMR, 그리고 tidyverse와 data.table 같은 범용 패키지가 보완되었습니다. 세션 정보와 패키지 버전은 엄격한 재현성을 보장하기 위해 파일에 기록됩니다.

멘델 무작위화의 핵심 가정
MR 프레임워크는 유전적 기구가 TL과 강하게 연관되어 있으며, TL과 갑상선독성증 간의 관계를 혼동하는 요인과 독립적이며, 대안 경로가 아닌 TL을 통해서만 갑상선독성증에 영향을 미친다는 전제를 제시합니다. 분석 계획은 각 변이체 강도를 정량화하고, 절개 기반 방법과 이상치 스크리닝을 이용한 불균형 수평 편면성 검사, 변이 간 이질성 평가, 그리고 TL에서 갑상선 독성증까지의 효과 방향을 공식 방향^{성 검사로} 확인함으로써 이러한 가정을 실질화합니다.

기기 선택 및 품질 관리
유전 기기는 TL GWAS에서 유전체 유의치 p< 5×^10-8로 선정되었습니다. 독립성을 보장하기 위해, 1000 게놈 프로젝트의 유럽 참조 데이터를 사용하여 10,000 kb 창 내에서 r² = 0.001에서 뭀집을 형성하여 연관 불균형을 해결했습니다; 변이체가 상관관계가 있는 경우, TL과 연관된 p값이 더 작은 변이가 유지되었습니다. 기기 강도는 각 변이별로 통계 F_i = \beta_{E,i}^{2}/\mathrm{SE}(\beta_{E,i})^{2}를 사용하여 요약되었으며, 기존의 F > 10 기준^13,14에 따라 평가되었습니다. 적절한 경우, 크기 N 샘플에서 여러 변종에 걸친 종합 기기 강도를 설명하기 위해 샌더슨-윈드마이어 근사가 인용되었으며, K 기기는 노출 분산 R^{2}: F_{\mathrm{SW}} = \{[R^{2}/(1-R^{2})]\,(N-K-1)\}/K ¹⁵를 설명하였다. 대립유전자 빈도가 0.5에 가까운 모호한 회문 변이는 대립유전자 빈도 정보가 명확한 정렬을 허용하지 않는 한 제외하였다.

조화화
각 기구별로 FinnGen에서 갑상선독성증과의 연관성을 추출하여 노출과 결과 전반에 걸쳐 동일한 효과 대립유전자를 나타내도록 정렬하였습니다. 조화는 대립유전자 불일치를 제거하고, 가닥 문제를 교정하며, 해결되지 않은 모호성을 가진 회문 SNP를 배제하여 유효한 Wald 비율 구성에 적합한 데이터셋을 생성했습니다.

1차 인과 추정
1차 분석에서는 역분산 가중(IVW) 추정량을 사용하여 SNP 특이적 Wald 비율을 통합된 인과 효과로 메타분석하였다¹⁶. 추정치는 TL 기준 1표준편차 증가당 오즈비로 보고되며, 고정 및 무작위 효과 모델 하에서 95% 신뢰구간을 도출합니다. 모델 선택은 이균성 진단에 의해 결정되었으며, 두 명세 모두 견고한 해석을 용이하게 하기 위해 제시되었습니다.

민감도 분석 및 강건성 검사
견고성은 MR-Egger 회귀를 사용해 인터셉트 항을 사용하여 방향성 불균형 다변성을, 균질 하에서 효율을 높이고 측정 오차를 완화하기 위한 최대우도 추정, 전역 이상치 검정과 왜곡 및 이상치 보정 효과를 추정하기 위한 MR-PRESSO, 약한 기기와 특이한 다면성에 대한 탄력성 추정을 위한 MR-RAPS를 사용해 평가했다¹⁷. 변이 간 이질성은 IVW^{프레임워크 18} 하에서 코크란의 Q 통계량을 사용하여 정량화하였다. 효과 방향은 Steiger 검정으로 조사되었으며, 이 검정은 노출과 결과 중 설명된 분산 비율을 비교하여 데이터가 TL이 갑상선독성증을 유발하는 것과 그 반대 중 어느 쪽인지 더 적합한지 판단합니다. MR-에거 요각은 방향성 다면형성에 대한 공식 시험으로 사용되었습니다. 전체 연관성이 단일 변이에 의해 주도되지 않았는지 확인하기 위해 생략 분석을 검토하였습니다.

역멘델 무작위화
역방향 인과관계를 탐색하기 위해 분석 파이프라인을 반복하여 갑상선독성증을 노출, TL을 결과로 삼았습니다. 동일한 기기 선택 기준, 조화 절차, IVW 1차 추정기, 민감도 분석이 적용되어 방향성에 관한 결론이 동일한 인과 틀 내에서 도출되었습니다.

다중 테스트, 전원 및 보고
주요 가설 검정은 TL이 갑상선독성증에 미치는 영향을 위한 IVW 추정량에 관한 것입니다. 민감도 추정기와 진단 검사는 지지 증거로 해석되며; p 값은 명확성을 위해 과학적 표기법으로 보고되며, 결론은 고립된 유의성 임계값보다는 방법 간 일관성을 강조합니다. 기기 강도 요약과 설명된 분산 비율은 표준 비중심성 공식에서 근사력을 제공하며, 검정력은 표본 크기, 기기 강도, 그리고 진정한 효과 크기에 따라 달라진다는 점을 인정합니다.

계산 재현성과 정확한 명령어
재현성은 기기 선택, 결과 추출, 조화, 1차 및 민감도 분석, 역방향 MR, 진단 출력, 분석 준비 파일 내보내기 등 R 명령어 전체 순서를 제공함으로써 보장됩니다. 스크립트는 기기 목록, 조화된 데이터셋, MR 추정 및 진단 요약에 대응하는 안정적이고 사람이 읽을 수 있는 CSV 파일을 작성합니다.

# R 4.3.1; TwoSampleMR 0.5.7
# 선택적 설치:
# install.packages(c("TwoSampleMR", "ieugwasr", "MRPRESSO", "mr.raps",
"RadialMR", "tidyverse", "data.table"))

라이브러리(TwoSampleMR)
도서관(Ieugwasr)
도서관(MRPRESSO)
도서관(Mr.Raps)
도서관(RadialMR)
도서관(타이디버스)
library(data.table)

# 노출: 텔로미어 길이 (IEU OpenGWAS)
exposure_id <- "Ieu-B-4879"

# 결과: 보고된 분석에 사용된 FinnGen 2021 갑상선 독성증 종말점
outcome_id <- "핀-b-갑상선 독성증"

# 유전체 전역 임계값과 엄격한 LD 응집을 가진 기기 선택
EXP <- extract_instruments(결과 = exposure_id, p1 = 5e-8, 클럼프 = TRUE, r2 = 0.001, kb = 10000)
익스피크F_stat <- (익스피크로스$베타.노출^2) / (익스피크시에이시오오스피어^2)
fwrite(경험치, "S1_instruments_TL.csv")

# 결과 추출과 조화화
아웃 <- extract_outcome_data(SNPS = exp$SNP, 결과 = outcome_id)
dat <- harmonise_data(경험치, 아웃, 행동 = 2)
fwrite(dat, "S2_harmonised_TL_vs_thyrotoxicosis.csv")

# 주요 MR 및 민감도 추정기
Res <- Mr(DAT, method_list = C("mr_ivw", "mr_ivw_fe", "mr_egger_regression",
"mr_weighted_median", "mr_raps", "mr_maxlik"))

het <- mr_heterogeneity(dat) # 코크란의 Q
플레오 <- mr_pleiotropy_test(dat) # 에거 가로채기
슈타이거 <- directionality_test(dat) # 슈타이거 방향성

# MR-PRESSO 전 세계 및 이상치 보정 추정치
mrpresso <- mr_presso(BetaOutresult = "beta.outcome", BetaExposure = "beta.exposure",
SdOutresult = "se.outcome", SdExposure = "se.노출",
OUTLIERtest = 참, DISTORTIONtest = 참,
data = dat, NbDistribution = 1000, SignifThreshold = 0.05) #

# 역방향 MR: 갑상선독성증(노출) -> TL (결과)
rev_exp <- extract_instruments(결과 = outcome_id, p1 = 5e-8, 클럼프 = 참, r² = 0.001, kb = 10000)
rev_out <- extract_outcome_data(SNPS = rev_exp$SNP, 결과 = exposure_id)
rev_dat <- harmonise_data(rev_exp, rev_out, 행동 = 2)

rev_res <- mr(rev_dat, method_list = c("mr_ivw","mr_ivw_fe","mr_egger_regression",

"mr_weighted_median", "mr_raps", "mr_maxlik"))
rev_het <- mr_heterogeneity(rev_dat)
rev_pleio <- mr_pleiotropy_test(rev_dat)
rev_steiger <- directionality_test(rev_dat)

# 아카이빙 및 피규/테이블 생성용 내보내기
write.csv(bind_rows(res), "S3_mr_results_primary.csv", row.names = FALSE)
write.csv(het, "S3_mr_heterogeneity.csv", row.names = FALSE)
write.csv(플레오, "S3_mr_pleiotropy_egger.csv", 행 이름 = 거짓)
write.csv(steiger, "S3_mr_steiger.csv", row.names = FALSE)
write.csv(bind_rows(rev_res), "S3_reverse_mr_results.csv", row.names = FALSE)
write.csv(rev_het, "S3_reverse_heterogeneity.csv", row.names = FALSE)
write.csv(rev_pleio, "S3_reverse_pleiotropy_egger.csv", row.names = FALSE)
write.csv(rev_steiger, "S3_reverse_steiger.csv", row.names = FALSE)

결과

기기 특성 및 가정 검사
텔로미어 길이 GWAS에서 129개의 유전체 전역 유의미하고 LD 독립적인 기기(p < 5×^10-8; r² < 0.001; 10,000 kb 창)가 유지되었습니다. per SNP F-통계량은 대다수 변이에서 기존 임계값을 초과하여 약한 기기 편향 위험이 제한적임을 나타냈다. 조화는 모호한 회문 변이체를 제거하고 대립유전자를 정렬하여 노출 효과와 결과 효과가 동일한 효과 대립유전자를 참조하도록 하여 Wald 비율 구성 및 통합 추정에 적합한 데이터셋을 만들었습니다.

주요 인과 추정치
모든 추정치는 방향성이 일관되어, 유전적으로 프록시 텔로미어 길이가 길수록 갑상선독성증 위험이 낮음을 시사하였다. 1차 IVW 분석은 OR = 0.689, 95% 신뢰구간(CI) 0.564-0.842, p < 1.0×^10-3을 보였다...

토론

이 멘델 무작위화(MR) 연구는 유전적으로 예측된 텔로미어가 길수록 갑상선독성증 위험을 감소시킨다는 증거를 제공합니다¹⁹. IVW, MR-Egger, 최대 우도, MR-PRESSO, MR-RAPS 전반에 걸친 방향성 일치 추정치는 강력한 보호 연관성을 나타내며, 방향성 검사는 역방향²⁰이 아닌 TL → 갑상선독성증 경로를 지지한다. 생물학적 타당성은 텔로미어가 염색체 완전성을 보존하고 세포 수명을 조절하는 역할과 일치합니다. 진행성 텔로미어 단축은 노화와^{세포 기능 장애}를 가속화하는 반면, 갑상선 중독증은 대사 교체 증가와 증식 신호 전달 증가로 특징지어집니다²². 이 높은 수요 상태에서는 더 긴 텔로미어가 산화 스트레스를 완충하고 DNA 손상을 제한하여 상피 및 난포 세포의 무결성을 유지하고 병리학적...

재료

이 논문에 사용된 재료 목록

이름	회사	카탈로그 번호	댓글
1000 게놈 프로젝트	연관 비평형(LD) 응집에 대한 참고 데이터.	1000 게놈 프로젝트	-
코크런' s Q 통계량	Statistical method used for heterogeneity testing in MR.	해당 없음	-
FinnGen 데이터셋	갑상선독성증에 대한 GWAS 요약 통계(사례 및 대조군).	핀젠 컨소시엄	2021년 발표, 핀-b-상선 독성증
IEU OpenGWAS 데이터셋	텔로미어 길이에 대한 GWAS 요약 통계.	IEU OpenGWAS	데이터셋 ID: ieu-b-4879
Ieugwasr 패키지	GWAS 데이터에 프로그래밍 접근을 위해 사용되는 R 패키지입니다.	크랜	버전 0.1.5
미스터랩스 패키지	약한 도구에서 멘델 무작위화의 강점 추정을 위한 R 패키지.	크랜	버전 0.1.2
MR-에거 회귀	방향성 다면형성을 검사하기 위한 감도 분석 방법.	해당 없음	-
MR-PRESSO 이상치 검사	MR 분석에서 이상치를 식별하고 교정하기 위해 사용되는 통계적 검정.	해당 없음	-
MRPRESSO 패키지	MR 분석에서 이상치 검출 및 왜곡 테스트에 사용되는 R 패키지.	크랜	버전 1.0
R 소프트웨어	모든 분석에 사용되는 통계 컴퓨팅 및 그래픽 환경.	R 통계 컴퓨팅 재단	버전 4.3.1
RadialMR 패키지	MR 결과의 방사형 시각화를 위한 R 패키지.	크랜	버전 1.1.2
TwoSampleMR 패키지	멘델 무작위화 분석에 사용되는 R 패키지입니다.	크랜	버전 0.5.7