テロメア長と甲状腺毒性症の関連:2サンプルメンデルランダム化研究からの洞察

サマリー

このプロトコルは、公開された遺伝的概要データを用いて、テロメア長が甲状腺毒性リスクに因果関係を与えるかどうかを評価する2サンプルのメンデルランダム化パイプラインを記述しています。機器選択、調和、一次推定、感度解析、再現可能なRコードを含み、透明性のある報告を支援する図形出力も含めてカバーしています。

要約

甲状腺中毒症は甲状腺ホルモンの過剰を特徴とする内分泌障害ですが、その全身的老化生物学との病因的関連は完全には解明されていません。テロメア長(TL)は細胞の老化とゲノムの安定性を反映し、複数の複雑な疾患に関与していることが示唆されています。遺伝的に予測されたTLが甲状腺中毒症リスクに与える因果関係を評価するため、2サンプルのメンデルランダム化(MR)試験を実施しました。TLの遺伝的手法は、ヨーロッパ系の祖先に関する大規模なゲノムワイド関連研究(GWAS)から導き出されました(n> 470,000人)。甲状腺毒性の要約統計は最新のFinnGenリリース(≈4,000件、対照群>210,000件)から取得されました。一次逆分散加重解析では、遺伝的に代理されたTLが長いほど甲状腺中毒リスクが低く、その効果の方向と大きさは補完推定量(MR-Egger、加重中央値/最大尤度、MR-PRESSO、MR-RAPS)で一貫していました。感度解析では方向性多面性の証拠は見られず、異質性を評価するためにコクランQを用いました。シュタイガー方向性検定はTLから甲状腺中毒への因果的流れを支持しました。

私たちの知る限り、本研究は現代のGWAS資源を用いて甲状腺毒性リスクとTLの因果関係を評価する初期のMR解析の一つであり、甲状腺機能亢進症関連表現型に焦点を当てた既存の証拠を拡張しています。これらの発見は、TLがインデックス化した細胞老化過程が甲状腺中毒の感受性に寄与する可能性を示唆し、甲状腺機能障害におけるテロメア生物学に関する今後の縦断的および機構的研究の動機付けとなります。

概要

甲状腺毒性症は、過剰な循環中の甲状腺ホルモンを特徴とする一般的な内分泌疾患です。典型的な症状には、体重減少、頻脈、イライラ、不安、振戦があり、重度の場合は不整脈、心不全、^骨量減少の素因があります。臨床的特徴や治療戦略は比較的確立されていますが、個体間の感受性や進行を形作る生物学的メカニズムは完全には定義されていません。病因は自己免疫性の高機能症(例:バセドウ病)から、毒性結節、甲状腺炎、外因性ホルモン曝露、薬物反応、腫瘍^{関連原因まで}多岐にわたります。これらの存在において、免疫活性化、酸化ストレス、ミトコンドリア生合成の変化が、ホルモン過剰と組織損傷や全身合併症を結びつけるもっともらしいメカニズム的テーマとして繰り返し現れます。

テロメア長(TL)は染色体末端を覆う二重のDNAタンパク質構造で、細胞の老化、複製履歴、ゲノム安定性を統合します。TLは細胞分裂とともに短縮され、酸化ストレスや慢性炎症によって加速されます。これらは^{甲状腺機能障害}でも顕著なプロセスです^3,4。短縮TLは心代謝障害、特定のがん、神経変^性疾患と関連しており、新たなデータは肥満や2型糖尿病などの内分泌表現型との関連を示唆しています。これらの表現型では、ホルモンおよび免疫代謝軸が複製ストレスと交差する^6。甲状腺ホルモンは基礎代謝率、ミトコンドリア機能、増殖シグナルを調節するため、TLダイナミクスは甲状腺毒性状態で特に重要であり、全身性高代謝や免疫調節障害がテロメア生物学に影響を与え、また影響を受ける可能性があります^7,8。

しかし、TLが甲状腺中毒症に因果関係を果たすかどうかは依然として不確かです。従来の観察研究は、混絡(例:生活習慣、併存症)や逆因果(TLに影響を与える疾患)に脆弱です。これらの制限に対処するため、2標本メンデルランダム化(MR)を用い、TLに強く関連した生殖細胞系遺伝変異を交絡なしの手法として利用し、^{疾患リスクへの}影響を検証します。ゲノムワイド関連研究は、厳密な因果評価を可能にするために、十分に強力なTLツールと大規模で表現型がよく分かれたアウトカムデータセットを提供しています。¹⁰。遺伝的疫学とテロメア生物学を統合することで、本研究は遺伝的に代理されたTLが甲状腺中毒症のリスクに影響を与えるかどうかを評価し、メカニズム経路を明確にし、TLをリスク階層化および将来の介入開発の候補バイオマーカーとして探求します。

私たちの知る限り、この分析は、当時の大規模なGWAS資源を用いてTLに関連する甲状腺中毒症の全体的なリスクを検証した初期のMR研究の一つであり、主に甲状腺機能亢進症関連表現型や自己免疫^{性甲状腺疾患}に焦点を当てた先行研究を拡張しています。

プロトコル

本研究は、遺伝的プロキシテロメア長(TL)が甲状腺中毒症リスクに因果的に影響を与えるかどうかを評価するために、2標本メンデルランダム化(MR)デザインを用いています。非識別化された公開されているGWASの要約統計のみを使用し、個人レベルのデータはアクセスされませんでした。青島市立病院の機関審査委員会は、この分析が公開された要約データのみに依拠しているため、さらなる審査の対象外と判断しました。すべての寄与ゲノムワイド関連研究は、元のプロトコルの一環としてインフォームドコンセントおよび倫理的承認を取得しました。分析は、TLに関連する生殖細胞系変異を交絡のないツールとして用い、多数性、異質性、効果の方向を問う感度手法を実装することで交絡を最小限に抑え、因果関係の逆転を図る構成が施されています。

データソース
TLの曝露データセットはIEU OpenWASから識別子ieu-b-4879で取得され、約472,174人のヨーロッパ系参加者を含みました。甲状腺中毒のアウトカムデータセットはFinnGenコンソーシアムから取得され、2021年のエンドポイントであるfinn-b-thyROTOXICOSISを用いて、4,142件の症例と213,693人の対照群を対象としました。これらのリソースは、SNPレベルの計測器を構築し、対立遺伝子を調和させ、確立されたMR推定量と因果効果を推定するために必要な効果的な推定と標準誤差を提供します。

ソフトウェアおよびコンピューティング環境
すべての解析はR(バージョン4.3.1)で実施されました。主な解析パッケージはTwoSampleMR(バージョン0.5.7)で、GIWASリソースへのプログラム的アクセス用に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に照らして判断されました。適切な場合、Sanderson–Windmeijer近似が参照され、サイズNのサンプルにおける複数のバリアントにわたる計測器強度を記述し、K計測器が曝露分散R^{2}: F_{\mathrm{SW}} = \{[R^{2}/(1-R^{2})]\,(N-K-1)\}/K ¹⁵を説明しました。対立遺伝子頻度が0.5付近の曖昧な回文変異は、対立遺伝子頻度情報が明確なアラインメントを認めない限り除外した。

調和
各器具について、FinnGenから甲状腺毒性症との関連を抽出し、曝露とアウトカムで効果量が同じ効果対立遺伝子を表すようにアラインメントしました。調和により対立遺伝子の不一致が除去され、鎖の問題が修正され、未解決の曖昧さを持つ回文SNPが除外され、有効なWald比構築に適したデータセットが生成されました。

一次因果推定
一次解析では逆分散加重(IVW)推定量を用いて、SNP特異的なWald比を統合因果効果にメタ解析しました¹⁶。推定値は、目標時間線(TL)で標準偏差1回増加するオッズ比として報告され、固定効果モデルおよびランダム効果モデルで95%信頼区間が導出されます。モデルの選択は異質性診断に基づいており、両方の仕様は堅牢な解釈を促進するために提示されています。

感度解析とロバストネスチェック
ロバスト性は、不均衡な方向性多数性の検定のためのMR-Egger回帰(インターセプト項付き)、均質性下での効率向上と測定誤差の軽減のための最尤推定、グローバル外れ値検定と歪みおよび外れ値補正効果の推定のためのMR-PRESSO、弱い機器や特異な多面性に対する耐性のある推定を提供するMR-RAPSを用いて評価しました¹⁷.変異間の異質性は、IVWフレームワーク¹⁸の下でコクランのQ統計量を用いて定量化されました。効果の方向は、曝露とアウトカムにおける分散の説明割合を比較し、TLによる甲状腺毒性の原因と比較するデータを判断するSteiger検定で検証されました。MR-エッガー切断は方向性多面力の正式な試験として用いられました。リーブワンアウト分析を用いて、全体の関連が単一の変異に駆動されていないことを確認しました。

逆メンデルランダム化
逆因果の可能性を探るために、曝露を甲状腺中毒、結果をTLとした分析パイプラインを繰り返した。同じ機器選択基準、調和手順、IVW一次推定量、感度分析が適用され、方向性に関する結論は同一の因果的枠組み内で導き出されました。

複数のテスト、電源、報告
主要な仮説検定は、甲状腺毒性に対するTLの影響に関するIVW推定量に関するものです。感度推定器や診断検査は支持的な証拠として解釈されます。p値は明確さのために科学的記譜法で報告され、結論は孤立した有意性閾値よりも手法間の一貫性を強調しています。機器の強度の要約と説明された分散の割合は、標準的な非中心性定式化における近似パワーを判断し、パワーはサンプルサイズ、機器の強度、真の効果の大きさに依存することを認めています。

計算再現性と正確なコマンド
再現性は、機器選択、結果抽出、調和、一次および感度解析、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)
ライブラリー(イエウグヴァースル)
図書館(MRPRESSO)
ライブラリー(Mr.Raps)
図書館(RadialMR)
ライブラリー(Tidyverse)
library(data.table)

# 露光:テロメア長(IEU OpenGWAS)
exposure_id <- 「IU-B-4879」

# アウトカム:報告解析に使用されたFinnGen 2021の甲状腺中毒症エンドポイント
outcome_id <- 「フィン・B・甲状腺中毒症」

# ゲノムワイド閾値と厳格なLD凝集による器具選択
経験値<- extract_instruments(結果=exposure_id、p1=5e-8、塊=真、r²=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と感度推定器
解像度 <- MR(DAT, method_list = c("mr_ivw","mr_ivw_fe","mr_egger_regression",
「mr_weighted_median」「mr_raps」「mr_maxlik」))

ヘット<- mr_heterogeneity(dat) # コクランのQ
プレオ <- mr_pleiotropy_test(DAT) # エッガー迎撃
シュタイガー <- directionality_test(dat) # シュタイガー方向性

# MR-PRESSO グローバルおよび外れ値補正推定
mrpresso <- mr_presso(BetaOutcome = "beta.outcome", BetaExposure = "beta.exposure",
SdOutcome = 「se.outcome」、SdExposure = 「se.exposure」、
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, action = 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からは、ゲノム全体で有意でLDに依存しない129個の機器(p < 5×^10-8;r² < 0.001;10,000 kbウィンドウ)が保持されました。Per-SNP F統計量は大多数のバリアントで従来の閾値を超えており、弱いインストゥルメントバイアスのリスクが限定的であることを示しました。調和により曖昧な回文変異が除去され、曝露効果とアウトカム効果が同じ効果対立遺伝子を参照するようにアライントされ、Wald比率構築およびプール推定に適したデータセットが得られました。

一次因果推定
すべての推定値は方向性が一貫しており、遺伝的に代理されるテロメア長が長いほど甲状腺中毒症のリスクが低いことが示されました。一次IVW分析ではOR = 0.689、95%信頼区間0.564-0.842、p<1.0×^10-3となりました。最尤推定量はOR = 0.689、95%信頼区間は0.576-0.8...

ディスカッション

このメンデルランダム化(MR)研究は、遺伝的に予測されるテロメアがより長いほど甲状腺毒性症のリスクを減らすという証拠を提供しています¹⁹。IVW、MR-Egger、最大尤度、MR-PRESSO、MR-RAPSにおける方向性一致の推定値は強固な保護関連を示し、方向性チェックは逆²⁰ではなくTL→甲状腺毒性経路を支持しています。生物学的な妥当性は、染色体の完全性を保ち細胞寿命を調節するテロメアの役割と一致しています。進行性のテロメア短縮は老化と細胞機能障害を加速させる²¹に対し、甲状腺毒性症は代謝ターンオーバーの加速と増殖シグナル伝達を特徴とする²²。この高需要状態では、より長いテロメアが酸化ストレスを緩衝しDNA損傷を抑えることで、上皮細胞や濾胞の完全性を維持し、病理^学的再構築を遅らせることができます。長いテロメアと心血管疾患や一部のが...

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
1000ゲノムプロジェクト	連鎖不平衡(LD)凝集の参考データ。	1000ゲノムプロジェクト	-
コクランs Q統計量	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