本プロトコルでは、個別のリスク層別化のために、歯周組織パラメータと従来のリスク因子を統合し、脳小血管病の負荷量を推定する探索的なノモグラムの開発および内部妥当性の検証について記述します。
本プロトコルでは、個別のリスク層別化のために、歯周組織パラメータと従来のリスク因子を統合し、脳小血管病の負荷量を推定する探索的なノモグラムの開発および内部妥当性の検証について記述します。
脳小血管病(CSVD)は脳卒中および認知機能低下の主要な原因ですが、CSVDの全体的な負荷に対する歯周病の寄与は依然として不明です。本研究では、高いCSVD負荷を推定するために、歯周パラメータを組み込んだ探索的な臨床ノモグラムの開発および内部妥当性の検証を目的としました。計234名を対象に、磁気共鳴画像法(MRI)を用いてCSVDの総負荷(スコア0–4)を評価しました。また、歯周炎の重症度と残存歯数を記録しました。least absolute shrinkage and selection operator(LASSO)回帰を用いて選択された最適な予測因子を多変量ロジスティック回帰に投入してノモグラムを構築し、その識別能、較正能、および臨床的有用性を評価しました。制限付き三次スプライン解析により、歯の喪失と高いCSVD負荷との間に線形の用量反応関係が認められました(非線形性のP = 0.332)。多変量解析の結果、高齢および高血圧が独立した予後因子として特定されましたが、重症歯周炎(P = 0.580)および重度の歯の喪失(P = 0.112)との関連は減弱し、調整後は高いCSVD負荷と独立して関連していませんでした。ノモグラムは中程度の識別能(曲線下面積 = 0.675)を示し、ブートストラップ法で検証された良好な較正能(平均絶対誤差 = 0.038)を示しました。決定曲線分析では、選択したリスク閾値全体で潜在的な臨床的有用性が示唆されましたが、外部妥当性の検証が行われていないため、解釈は探索的なものに留まります。重症歯周炎および歯の喪失は、単変量解析では高いCSVD負荷との関連を示しましたが、多変量調整後の独立した予測因子ではありませんでした。この探索的なノモグラムは、個別化されたリスク層別化のための予備的な視覚化フレームワークを提供しますが、臨床導入の前にはさらなる外部妥当性の検証が必要です。
脳小血管病(CSVD)は、脳の慢性的かつ潜行的な微小血管障害であり、高齢者における脳卒中、認知機能低下、および運動機能障害の主要な血管性原因であると広く認識されています1。従来、白質高信号(WMH)、ラクナ梗塞、脳微小出血(CMBs)、および血管周囲腔拡大(EPVS)を含むCSVDの個々の磁気共鳴画像(MRI)マーカーは、独立して評価されてきました。しかし、これらの神経画像表現型は頻繁に共存しており、特に内皮機能不全や微小血管漏出といった共通の病理学的メカニズムによって引き起こされます2。そのため、近年のコンセンサスステートメントでは、CSVD総負荷量スコアリングシステムの利用が推奨されています3。この統合的な指標は、個々の画像マーカーよりも脳微小血管損傷の累積的な影響をより包括的に反映し、全般的な神経学的悪化の優れた指標となります4。
CSVDの病因はまだ完全には解明されていませんが、慢性の低グレードな全身性炎症が、脳血管内皮傷害および血液脳関門破綻の主要な要因であるとの認識が高まっています5。このような背景から、歯周組織に影響を及ぼす、バイオフィルム誘発性のジスバイオーシスを伴う有病率の高い炎症性疾患である歯周炎が、口腔-脳軸の枠組みにおいて、神経血管疾患の修正可能なリスク要因となる可能性が示されています6。重度の歯周炎は、進行性の歯槽骨吸収と最終的な歯の喪失をもたらすだけでなく、歯周病原体および炎症性サイトカインの持続的な貯蔵庫としても機能します7。これらの炎症性メディエーターは全身循環に入り、それによって血管リモデリングや動脈硬化の形成に寄与する可能性があります8。
近年の疫学研究により、重度歯周炎と虚血性脳卒中のリスク増加との関連が示されています9。それにもかかわらず、大脳小血管病(CSVD)の累積的な負担との関係に関するエビデンスは、依然として限定的かつ断片的です。さらに、歯周病の最終的な臨床的結果である喪歯は、認知機能の低下や心血管死と関連していることが報告されています10。しかし、喪歯と脳血管損傷との用量反応関係については、十分に解明されていません。具体的には、残存歯数が臨界閾値に達した後にCSVDのリスクが増加するのか、あるいは連続的な線形用量反応パターンに従うのかは不明です。
歯周組織の破壊と脳微小血管損傷との関連性は生物学的に妥当であると考えられるが、この関係の根底にある全身性の免疫炎症メカニズムについてはさらなる特性解析が必要である。単離した白血球サブセット数や高感度C反応性蛋白を含む従来の炎症バイオマーカーでは、臨床現場における宿主の免疫炎症反応の複雑さを十分に捉えられない可能性がある。末梢血の血小板数、好中球数、およびリンパ球数から算出される全身性免疫炎症指数(SII)は、近年、包括的な炎症バイオマーカーとして注目されている11。SIIは、全身性の炎症活性と免疫状態のバランスを反映する。SII値の上昇は、急性脳卒中の重症度や不良な機能的転帰を予測することで、さまざまな脳血管障害において予後予測能を示すことが証明されており、また、全般的なCSVD負荷、認知機能障害、および歯周炎の段階的診断とも関連している12,13。歯周炎は慢性の全身性炎症を誘発するため、歯周パラメータとともにSIIを評価することで、潜在的な脳微小血管損傷のリスクがある個人の特定を改善できる可能性がある。
口腔-脳軸への関心が高まっているものの、現在の文献にはいくつかの方法論的な限界がある。ほとんどの研究は、個々のCSVD画像マーカーまたは単一の口腔健康パラメータに焦点を当てており、歯周炎の重症度、その最終的な臨床的結果である歯の喪失、および全身的な免疫炎症反応を統合した評価は行われていない。さらに重要なことに、臨床現場では実用的なリスク評価ツールが必要とされているが、口腔および全身の炎症プロファイルの組み合わせに基づいてCSVDの高リスク患者を特定するためのノモグラムのような視覚的ツールは依然として限られている。既存のCSVDリスク評価手法は主に個々のバイオマーカーや従来の人口統計学的要因に依存しており、全身性炎症と口腔健康との間の多因子的な相互作用を十分に捉えられていない可能性がある。さらに、信頼性の高い予測モデルの開発には、強固な変数選択法が必要である。ステップワイズ回帰などの従来の手法は、高次元の臨床データセットに適用した場合に不安定になる可能性があり、相関のある予測因子の影響を受けることがある。LASSO(Least Absolute Shrinkage and Selection Operator)回帰アプローチは、係数へのペナルティを適用することで過剰適合を抑制し、モデル開発において最も情報量の多い予測因子を特定することで、これらの課題を解決する14。
これらの知識の空白に着目し、本レトロスペクティブ横断研究は2つの主要な目的をもって設計されました。第一に、制限付き3次スプライン(RCS)解析を用いて、残存歯数と高度なCSVD負荷との間の用量反応関係を特徴付け、潜在的な非線形閾値効果および線形傾向を評価することを目指しました。第二に、LASSOベースの変数選択を用いて探索的な臨床ノモグラムを構築し、従来のリスク因子に加えて歯周パラメータによるリスク再分類の増分的な寄与を評価することを目的としました。私たちは、歯周状態、歯の喪失、およびSIIを多次元モデルに統合することで、累積的なCSVD負荷を推定するための探索的な可視化フレームワークを提供し、個別化されたリスク層別化を改善できるという仮説を立てました。臨床現場において、提案されたノモグラムは、早期の脳血管評価や多職種による歯科評価の恩恵を受ける可能性のある患者を神経内科医が特定するための、迅速かつ非侵襲的な補助ツールとして意図されています。しかし、本ノモグラムは探索的な性質を持ち、識別能が控えめであること、および外部妥当性の検証が行われていないことから、磁気共鳴画像法(MRI)に基づく診断の代わりとして検討されるべきではありません。加えて、活動性の急性感染症、重度の自己免疫疾患、または最近の外傷がある患者では、これらの状態が全身性炎症バイオマーカーに大幅に影響を与え、リスク推定に影響を及ぼす可能性があるため、その適用は制限される場合があります。
この後方視的横断的関連研究では、202年1月から2024年12月の間に、めまい、頭痛、脳卒中スクリーニング、または認知機能評価のために常州第二人民病院の神経内科に入院した連続的な患者を登録しました。本研究が後方視的かつ非介入的であり、完全に匿名化されたデータを使用しているため、インフォームドコンセントは免除されました。本研究は常州第二人民病院の臨床医療技術倫理委員会(IRB承認ID:[2023] YLJSA069)によって承認され、施設ガイドラインに従って実施されました。
研究デザインおよび患者集団
当初、850人の患者を評価しました。除外基準は次のように厳格に定義されました:(1) 頭部MRIシーケンスが不完全である(特に、微小出血の正確な評価を妨げる感受性強調画像[SWI]の欠如)、(2) 詳細な歯周ポケット測定記録または全顎の歯数データがない、(3) 血液学的データが不完全でSIIの算出ができない、(4) 大規模脳卒中、脳腫瘍、外傷性脳損傷、または中枢神経系感染症の既往がある、(5) 全身性炎症マーカーへの混絡影響を最小限にするため、過去1ヶ月以内に活動性急性感染症、重症自己免疫疾患、または悪性腫瘍がある。データの整合性を確保するため、主要な臨床放射線学的パラメータのいずれかが欠落している個体を除外する完全ケース分析手法を用いました。スクリーニング後、適格な234人の参加者を分析に含めました。129件の高負荷イベントを含むこのサンプルは、安定した多変量モデリングのためのイベント数/変数(EPV)>10という経験則を満たしていました。
歯周状態および炎症の評価
校正済みの2名の歯周病専門医が、標準化された歯周プローブを用いてすべての臨床口腔評価を実施した。
研究開始前に、両検者は標準化トレーニングを受け、検者間信頼性は極めて良好であった(Cohen’s κ = 0.82)。プロービングポケット深さ(PPD)および臨床的付着レベル(CAL)を、1歯あたり6箇所で測定した。Tonettiら15によって記載された2018年の分類フレームワークに基づき、歯周病を軽度/なし、中等度、または重度の3段階に分類した。具体的に、参加者は以下のようにカテゴリー分けした:(1) 軽度/なしの歯周炎(歯周組織の健康およびステージI歯周炎を含む)。これは、最大付着喪失部位における歯間臨床的付着レベル(CAL)が2 mm以下かつプロービングポケット深さ(PPD)が4 mm以下で、歯周炎に関連する歯の喪失がないと定義した。(2) 中等度歯周炎(ステージII歯周炎)。これは、歯間CALが3–4 mm、最大PPDが5 mm以下で、歯周炎による喪失歯が4本以下であると定義した。(3) 重度歯周炎(ステージIII/IV歯周炎)。これは、歯間CALが5 mm以上、PPDが6 mm以上、および/または歯周組織破壊により4本以上の歯を喪失していることと定義した15。
また、全参加者について、残存歯数についても記録した。機能的な歯を最低20本維持することは、基本的な口腔機能を保持し、口腔の健康な老化を実現するための臨床的な指標として広く認められているためである。16,17,18未処置の歯数(raw tooth count)は、二値変数(バイナリ変数)に変換された:重度歯牙喪失(<重度歯欠損(20本未満)対非重度歯欠損(20本以上)。
全身性炎症に関しては、入院時に採取した空腹時血液検体からSIIを算出した。静脈血をEDTAチューブに採取し、病院の標準的な臨床検査プロトコルに従って処理した後、自動血球計数装置を用いて好中球数、リンパ球数、および血小板数を測定した。指標は以下のように算出した:

SII値の分布が著しく偏っていることを考慮し、統計解析に組み込む前に自然対数変換(Log_SII)を適用した。
MRI撮像と総CSVD負荷
すべての被験者を対象に、標準ヘッドコイルを備えた3.0 T MRIスキャナーを用いて、標準化された頭部MRI検査を実施した。撮像プロトコルには、スライス厚5 mmで取得したT1強調画像、T2強調画像、液体減衰反転回復(FLAIR)画像、および磁気感受性強調画像(SWI)のシーケンスを含めた。標準的な撮像パラメータは以下の通りに設定した:T1強調画像(繰り返し時間 [TR] = 2,0 ms、エコー時間 [TE] = 9 ms、視野 [FOV] = 230 × 230 mm2、マトリクスサイズ = 256 × 256)、T2強調画像(TR = 4,50 ms, TE = 85 ms, FOV = 230 × 230 mm2, マトリクスサイズ = 256 × 256)、液体減衰反転回復(FLAIR)画像(TR = 8,50 ms, TE = 120 ms、反転時間 [TI] = 2,40 ms, FOV = 230 × 230 mm2, マトリクスサイズ = 256 × 256)、および磁気感受性強調画像(SWI)(TR = 28 ms, TE = 20 ms, フリップ角 = 15°, FOV = 230 × 230 mm2, マトリクスサイズ = 256 × 256)。すべてのシーケンスは、スライス厚5 mm、スライス間ギャップ1.0 mmで取得した。使用した個別のスキャナーについては、「材料一覧」に記載している。CSVDの撮像マーカーは、参加者の臨床データを伏せられた2名の神経放射線科医によって独立して評価された。CSVD総負荷スコアの検者間信頼性は高く(Cohen’s κ = 0.85)、不一致が生じた場合は、3人目の上級医師との合意によって解決した。
確立された国際的なコンセンサス1,3に基づき、4つの神経画像所見を評価することで、累積的CSVD負荷を0から4のスケールで定量化した。以下のMRI所見ごとに1点ずつ割り当てた:(1) 1つ以上のラクナ梗塞、(2) 1つ以上の脳微小出血(CMBs)、(3) 基底核における中等度から重度の拡大血管周囲腔(EPVS)(グレード ≥2)、(4) 重度の白質高信号(WMH)(深部白質でFazekasスコア ≥2、または脳室周囲領域で3と定義)。診断定義は、神経画像における血管変化の報告基準(STRIVE-1)国際コンセンサス1に厳密に従った。白質高信号(WMH)はFazekasスケール19を用いてグレード分けし、基底核の拡大血管周囲腔(EPVS)はPotterらによって記述され、妥当性が確認された4段階の視覚的評価スケール20を用いてグレード分けした。リスク評価モデルの開発にあたり、CSVD負荷の合計スコアが ≥2 の患者をCSVD高負荷群として分類した。この閾値は、認知機能の低下および死亡率の加速と一貫して関連しているためである21,2,23。
共変量データの収集
患者の人口統計学的特性および臨床歴は、施設内の電子健康記録から抽出した。収集した変数には、年齢、性別、体格指数(BMI)、喫煙状況、およびアルコール摂取量が含まれる。さらに、脂質異常症、糖尿病、高血圧、冠動脈疾患(CAD)、心筋梗塞(MI)の既往、および虚血性脳卒中の既往を含む、心代謝性および血管系の併存疾患を各参加者について記録した。
すべての併存疾患は、確立された臨床ガイドラインに従い、入院診療録に記載された医師の診断、薬剤歴、および入院時のルーチンな臨床検査および画像評価の確認を通じて確定された。高血圧は、収縮期血圧140 mmHg以上、拡張期血圧90 mmHg以上、または現在血圧降下薬による治療を受けていることと定義した。糖尿病は、空腹時血糖値7.0 mmol/L以上、糖化ヘモグロビン(HbA1c)6.5%以上、または血糖降下薬の使用と定義した。脂質異常症は、空腹時血清脂質の異常または現在脂質低下療法を受けていることによって定義した。CAD、心筋梗塞(MI)の既往、および虚血性脳卒中の既往は、電子カルテシステムで利用可能な記載済みの臨床歴および過去の神経血管または心臓画像診断報告書を通じて確認された。
統計解析およびモデル開発
統計解析を行う前に、連続変数の正規性をShapiro–Wilk検定を用いて評価した。カテゴリー変数の記述統計は件数(パーセンテージ)としてまとめ、適宜、ピアソンのカイ二乗検定またはフィッシャーの正確確率検定を用いて比較した。正規分布を示す連続変数は平均値 ± 標準偏差で表し、Studentのt検定を用いて比較した一方、非正規分布の変数は中央値(四分位範囲)で表し、Mann-Whitney U検定を用いて比較した。残存歯数と高度CSVD負荷のリスクとの間に潜在的な非線形用量反応パターンがあるかを調べるため、年齢、高血圧、およびLog_SIIで調整したRCSモデルを構築し、歯数分布の第5、35、65、95パーセンタイルに4つのノットを配置した。非線形性の正式な検定は分散分析(ANOVA)を用いて行い、閾値効果と連続的な線形トレンドのどちらが観察された関連性を最も適切に説明するかを判定した。
相関のある予測変数が存在する場合の変数選択を行うため、すべてのベースライン変数をLASSO回帰モデルに投入しました。最適なチューニングパラメータ(λ)は10分割交差検証を用いて決定し、最小二項偏差基準(λmin)に従って選択しました。係数が非ゼロである変数を、その後の解析のために保持しました。選択された変数はその後、多変数ロジスティック回帰モデルに投入され、高度なCSVD負荷との関連におけるオッズ比(OR)および95%信頼区間(CI)を推定しました。最終的な多変数モデルを用いて、個別のリスク推定のための探索的な臨床ノモグラムを構築しました。
歯周病パラメータの増分価値を評価するため、ベースラインモデル(年齢、高血圧、およびLog_SIIを含む)を、歯周病のグレードと歯数を組み込んだ拡張モデルと比較しました。モデルの性能は、複数の相補的な指標を用いて評価しました。識別能は、受信者動作特性曲線の下面積(AUC)を用いて評価し、モデル間の差はDeLong検定を用いて比較しました24。最大Youden指数によって決定された最適カットオフ値において、感度、特異度、陽性的中率(PPV)、および陰性的中率(NPV)も算出しました。拡張モデルによってもたらされるリスク再分類の増分を定量化するため、連続的な正味再分類改善率(NRI)および統合識別改善度(IDI)を算出しました25。較正能は、rmsパッケージを用いて実施した1,00回のブートストラップ再標本抽出から生成された較正プロットを用いて評価し、予測リスクと観察リスクの間の平均絶対誤差を算出しました。潜在的な過学習を考慮し、モデル性能の不偏な評価を行うため、ブートストラップ検証には楽観度補正を含め、バイアス補正後のモデル性能と見かけ上のモデル性能の両方を報告しました。決定曲線分析(DCA)を行い、一連の閾値確率における正味の利益を推定することで、潜在的な臨床的有用性を評価しました。すべての統計解析および視覚化は、統計ソフトウェアおよびrms、glmnet、pROC、PredictABEL、dcurvesパッケージを用いて行いました。統計的有意性は、両側 P 値 <0.05 と定義しました。
研究デザインおよび患者集団
重大な混絡疾患があること、または臨床放射線学的データが欠損していることを理由に616名を除外した結果、最終的な解析には神経学的評価のために入院した適格な参加者234名が含まれた (Figure 1)。参加者は、合計CSVDスコアに基づいて、高負荷群(n = 129)と低負荷群(n = 105)に層別化された。両群のベースラインの人口統計学的および臨床的特性をTable 1に示す。CSVD負荷が高い参加者は、CSVD負荷が低い参加者よりも高齢であり(67.71 ± 7.41 years vs. 64.91 ± 7.74 years, P = 0.05)、高血圧の有病率が高かった(63.6% vs. 48.6%, P = 0.024)。既往の虚血性脳卒中の有病率も高負荷群で高かったが、その差は統計的有意性に達しなかった(21.7% vs. 12.4%, P = 0.076)。口腔健康パラメータに関しては、高負荷群で重度歯周炎の有病率が高く(52.7% vs. 16.2%, P < 0.01)、それに伴い平均プロービングポケット深度(PPD; 4.8 mm vs. 2.54 mm, P < 0.01)および臨床的アタッチメントレベル(CAL; 4.54 mm vs. 2.12 mm, P < 0.01)が大きかった。また、高負荷群の参加者は保持天然歯数が少なかった(18.39 ± 8.15 vs. 24.36 ± 6.21, P < 0.01)。さらに、Log_SII値は低負荷群よりも高負荷群で高かった(6.48 ± 0.52 vs. 6.27 ± 0.45, P = 0.01)。喫煙状況、体格指数(BMI)、および性の分布については、群間に有意な差は認められなかった。

図1研究参加者選定のフローチャート。 後方視的横断研究における参加者選定を示すフローチャート。202年1月から2024年12月の間に、めまい、頭痛、脳卒中スクリーニング、または認知機能評価のために常州第二人民病院神経内科に入院した連続的な患者をスクリーニングした。除外基準を順次適用した後、参加者を脳小血管病(CSVD)負荷の低負荷群と高負荷群に層別化した。 この図の拡大表示を見るには、ここをクリックしてください。
| 変数 | 低CSVD負荷 (n = 105) | 高CSVD負荷 (n = 129) | P値 |
| 年齢 (years)、平均 ± SD | 64.91 ± 7.74 | 67.71 ± 7.41 | 0.05 |
| 男性、n (%) | 49 (46.7) | 69 (53.5) | 0.358 |
| ボディマス指数 (kg/m²)、平均 ± SD | 24.24 ± 2.92 | 23.62 ± 3.12 | 0.122 |
| 高血圧、n (%) | 51 (48.6) | 82 (63.6) | 0.024 |
| 糖尿病、n (%) | 27 (25.7) | 30 (23.3) | 0.76 |
| 脂質異常症、n (%) | 5 (52.4) | 59 (45.7) | 0.358 |
| 冠動脈疾患、n (%) | 20 (19.0) | 26 (20.2) | 0.87 |
| 心筋梗塞、n (%) | 4 (3.8) | 4 (3.1) | >0.99* |
| 脳虚血性卒中の既往、n (%) | 13 (12.4) | 28 (21.7) | 0.076 |
| 歯周炎グレード、n (%) | <0.01 | ||
| 健康/軽度 | 45 (42.9) | 21 (16.3) | |
| 中等度 | 43 (41.0) | 40 (31.0) | |
| 重度 | 17 (16.2) | 68 (52.7) | |
| 残存歯数、平均 ± SD | 24.36 ± 6.21 | 18.39 ± 8.15 | <0.01 |
| 平均プロービングポケット深さ (PPD, mm)、中央値 (IQR) | 2.54 (2.12–4.13) | 4.88 (3.65–6.21) | <0.01 |
| 平均臨床的付着レベル (CAL, mm)、中央値 (IQR) | 2.12 (0.89–3.76) | 4.54 (3.12–6.45) | <0.01 |
| Log_SII、平均 ± SD | 6.27 ± 0.45 | 6.48 ± 0.52 | 0.01 |
| 高感度C反応性タンパク (hs-CRP, mg/L)、中央値 (IQR) | 1.45 (0.76–2.54) | 2.89 (1.56–4.32) | <0.01 |
| ホモシステイン (Hcy, µmol/L)、中央値 (IQR) | 12.34 (9.87–15.65) | 15.76 (12.45–19.82) | <0.01 |
表1:脳小血管病(CSVD)負荷に基づく研究集団のベースライン特性。連続変数は、適宜、平均値 ± 標準偏差(SD)または中央値(四分位範囲 [IQR])で示す。カテゴリー変数は、数(パーセンテージ)で示す。* 細胞数が少ないため、Fisherの直接検定を用いて算出したP値。特に断りのない限り、その他のすべてのカテゴリー比較はPearsonのカイ二乗検定を用いて行った。略語:BMI、身体質量指数;CAL、臨床的アタッチメントレベル;CAD、冠動脈疾患;CSVD、脳小血管病;DM、糖尿病;Hcy、ホモシステイン;hs-CRP、高感度C反応性タンパク質;HTN、高血圧;IQR、四分位範囲;Log_SII、全身性免疫炎症指数の自然対数;MI、心筋梗塞;PPD、プロービングポケット深さ;SD、標準偏差;SII、全身性免疫炎症指数。
歯周状態および炎症の評価
歯周炎のグレード別における総CSVD負荷スコアの分布を図2に示す。CSVDスコアが高い参加者の割合は、歯周炎の重症度が増すにつれて増加し、コクラン・アーミテージ検定により有意な傾向が示された(P < 0.01)。残存歯数と高CSVD負荷の可能性との関連について、さらにRCS分析を用いて評価した(図3)。年齢、高血圧、およびLog_SIIで調整した後、正式な検定により、非線形の閾値効果ではなく線形の用量反応パターンが支持された(非線形性のP = 0.32)。残存歯数の減少に伴い、高CSVD負荷の推定確率は漸進的に増加した。ただし、残存歯数が最も少ない群における推定値については、信頼区間が広いため、慎重に解釈する必要がある。

図2歯周病のグレード別の総脳小血管病(CSVD)負荷の分布。 全CSVD負荷スコア(0~4)における、健康・軽度、中等度、および重度歯周炎の割合の分布を示す積層棒グラフ。Cochran-Armitageトレンド検定により、CSVD負荷の増大に伴い歯周炎の重症度が増加する有意な傾向が示された(P トレンドについて <0.001). こちらの図の拡大版を表示するには、ここをクリックしてください。

図3残存歯数と高度脳小血管病(CSVD)負荷との関連に関する制限付き三次スプライン分析。 年齢、高血圧、および全身性免疫炎症指数(Systemic Immune-Inflammation Index)の自然対数(Log_SII)で調整後における、残存歯数と高度な脳小血管疾患(CSVD)負荷の対数オッズとの関連を示す制限付き立方スプライン(RCS)モデル。赤の実線は推定された対数オッズを、網掛け領域は95%信頼区間(CI)を表す。非線形性の検定は統計的に有意ではなかった(P (非線形性のp値 = 0.332であり)、線形相関が支持された。歯数が最も少ない群の推定値については、信頼区間が広いため、解釈に注意が必要である。 この図の拡大版を表示するには、ここをクリックしてください。
統計解析およびモデル開発
ベースラインパラメータの間で変数選択を行うため、10分割交差検証を用いた最小絶対収縮選択演算子(LASSO)回帰を使用した(補足図1)。この手順では、年齢、高血圧、歯周炎のグレード、および残存歯数をコア特徴量として保持した。それに続く多変量ロジスティック回帰モデル(表2)、年齢(オッズ比 [OR] = 1.06、95% 信頼区間 [CI]: 1.02–1.10、 P = 0.08) および高血圧(OR = 2.02, 95% CI: 1.16–3.55, P = 0.014) は、高いCSVD負荷と独立して関連していた。重度の歯周炎(OR = 1.73, 95% CI: 0.25–12.19, P = 0.580) および残存歯数(OR = 0.94, 95% CI: 0.87–1.01, P (p = 0.12)は統計的有意性に達しておらず、多変数調整後にそれらの関連性が減弱したことが示された。これらの歯周病パラメータは、探索的なプロファイリング項目として最終モデルに保持された。重症歯周炎と高いCSVD負荷との関連に関するサブグループ解析の結果を、以下に提示する。 補足図2関連性の方向性は、事前に規定したサブグループ間で概ね一貫していた。具体的には、オッズ比(OR)は、女性で1.8(95% CI: 0.19–18.57)、男性で1.76(95% CI: 0.12–25.14)であった。また、高血圧を有する参加者では2.4(95% CI: 0.21–28.53)、高血圧のない参加者では1.48(95% CI: 0.1–19.34)であり、65歳以上の参加者では2.12(95% CI: 0.23–19.82)、それ未満の参加者では1.45(95% CI: 0.10–21.05)であった。 <65歳。性別に関して統計的に有意な相互作用は認められなかった(P相互作用 = 0.916)、高血圧の状態(P相互作用 = 0.627)、または年齢(P相互作用 = 0.684)であり、これらの探索的サブグループ解析において、効果修飾の根拠は認められなかった。
| 変数 | オッズ比 (OR) | 95%信頼区間 (CI) | P 値 |
| 年齢 | 1.06 | 1.02–1.10 | 0.008 |
| 高血圧 | 2.02 | 1.16–3.55 | 0.014 |
| Log_SII | 0.35 | 0.07–1.81 | 0.214 |
| 保持歯数 | 0.94 | 0.87–1.01 | 0.112 |
| 歯周炎のグレード | |||
| 健常/軽症(対照) | 1 | — | — |
| 中程度 | 1.41 | 0.51–3.90 | 0.51 |
| 重度 | 1.73 | 0.25–12.19 | 0.58 |
表2:高度な脳小血管疾患(CSVD)負荷に関連する因子の多変数ロジスティック回帰分析。多変数ロジスティック回帰モデルに含まれる変数は、LASSO(least absolute shrinkage and selection operator)回帰を用いて選択された。健康/軽度歯周炎を対照群とした。残存歯数のオッズ比(OR)は、残存歯が1本増えるごとに高度なCSVD負荷となるオッズがどのように変化するかを示している。略語:CI, 信頼区間; CSVD, 脳小血管疾患; LASSO, least absolute shrinkage and selection operator; Log_SII, 全身性免疫炎症指数(SII)の自然対数; OR, オッズ比。
個別化されたCSVDリスク評価のための探索的視覚化フレームワークとして、最終的な多変数モデルから臨床ノモグラムを構築した(Figure 4)。このノモグラムは、年齢、高血圧、Log_SII、歯周炎グレード、および残存歯数が、高CSVD負荷の推定確率にどのように寄与するかを示している。歯周組織パラメータの潜在的な増分価値を評価するため、年齢、高血圧、Log_SIIで構成されるベースラインモデルと、さらに歯周炎グレードと残存歯数を組み込んだ拡張モデルを比較した(Table 3)。受信者動作特性分析の結果、AUCは拡張モデルで0.68(95% CI: 0.621–0.75)、ベースラインモデルで0.656(95% CI: 0.587–0.726)であり、DeLong検定によればその差は統計的に有意ではなかった(P = 0.263)(Figure 5A)。拡張モデルでは、連続的純再分類改善度(NRI)が0.454(P < 0.01)、IDIが0.052(P < 0.01)となった。校正(キャリブレーション)は、1,00回のブートストラップ再サンプリングを用いて評価した。校正プロットでは、推定イベント確率と観察イベント確率の一致が示され、平均絶対誤差は0.038であった(Figure 5B)。DCAは、選択した閾値確率において拡張モデルに潜在的なネットベネフィットがあることを示した(Figure 6)。しかし、AUCの改善が統計的に有意ではなく、外部妥当性の検証が行われていないため、この明らかなベネフィットは慎重に解釈されるべきである。
| 指標 | ベースラインモデル | 拡張モデル | P 値 / 改善 |
| AUC (95% CI) | 0.656 (0.587–0.726) | 0.688 (0.621–0.755) | 0.263 (DeLong法) |
| 感度 | — | 0.512 | — |
| 特異性 | — | 0.8 | — |
| 陽性的中率 (PPV) | — | 0.759 | — |
| 陰性予測値 (NPV) | — | 0.571 | — |
| 連続的正味再分類改善度 (Continuous NRI) | — | 0.454 | <0.001 |
| 統合識別改善量 (IDI) | — | 0.052 | <0.001 |
表3:ベースラインモデルと拡張モデルの予測性能の比較。 ベースラインモデルには、年齢、高血圧、およびLog_SIIを含めた。拡張モデルには、さらに歯周炎のグレードと残存歯数を組み込んだ。モデルの識別能は、受信者動作特性曲線下面積(AUC)を用いて評価し、AUC間の差はDeLong検定を用いて比較した。感度、特異度、陽性的中率(PPV)、および陰性的中率(NPV)は、最大Youden指数によって決定された最適なカットオフ値を用いて算出した。拡張モデルの予測能の向上を評価するために、連続ネット再分類改善度(NRI)および統合識別改善度(IDI)を用いた。略語:AUC、受信者動作特性曲線下面積;CI、信頼区間;IDI、統合識別改善度;Log_SII、全身性免疫炎症指標の自然対数;NPV、陰性的中率;NRI、ネット再分類改善度;PPV、陽性的中率。

図 4高度な脳小血管病(CSVD)負荷の確率を推定するためのノモグラム。 最終的な多変数ロジスティック回帰モデルから作成されたノモグラム。個々の脳小血管病(CSVD)負荷が高い確率を推定するには、各予測因子の患者値を特定し、「Points(点数)」軸で対応するスコアを割り当て、個々のスコアを合計して合計点数を算出します。その合計点数を確率スケールに投影することで、CSVD負荷が高いという予測リスクを推定します。 この図の拡大版を表示するには、こちらをクリックしてください。

図5予測モデルの判別能および校正。 (A) ベースラインモデルと拡張モデルを比較した受信機動作特性(ROC)曲線。モデルの判別能はROC曲線下面積(AUC)を用いて評価し、モデル間の差異はDeLongテストを用いて評価した。(B) 1,0回のブートストラップ再標本抽出を用いて作成した拡張モデルのキャリブレーションプロット。外見上のキャリブレーション曲線、バイアス補正済みキャリブレーション曲線、および理想的なキャリブレーション曲線を、信頼限界(CL)と共に示した。キャリブレーション性能は、予測確率と観察確率の間の平均絶対誤差を用いて評価した。 この図の拡大版を表示するには、ここをクリックしてください。

図6予測モデルの決定曲線分析。 決定曲線分析 (DCA閾値確率の範囲におけるベースラインモデルと拡張モデルの正味の便益(net benefit)の比較。 「全員治療(Treat All)」および「誰も治療しない(Treat None)」戦略をリファレンス曲線として示している。 この図の拡大版を表示するには、ここをクリックしてください。
全体として、結果は歯周病の重症度、歯の喪失、および高いCSVD負荷との間に有意な未調整の関連があること、ならびに、残存歯数の減少とCSVDリスクの増加との間に線形的な関連があることを示した。それにもかかわらず、多変数調整後において、重度の歯周炎および残存歯数は、高いCSVD負荷と独立した関連を示さなかった。これらの知見は、予備的なリスク可視化モデルにおける歯周パラメータの探索的な利用を支持するが、それらを独立した予測因子として解釈することや、ノモグラムを臨床的に活用可能なスクリーニングツールとして利用することは支持しない。
データの可用性:
本研究で報告された結果の根拠となる、匿名化された参加者レベルのデータセットは、以下として提供されています。 補足表 1 また、本原稿とともに提出されています。
付録図1. Least Absolute Shrinkage and Selection Operator (LASSO) 回帰を用いた変数選択。 二項偏差に基づく10分割交差検証を用いて、チューニングパラメータ (λ) を最適化した。垂直の点線は、最小基準および1標準誤差 (1-SE) 基準に対応する最適な λ 値を示す。こちらのリンクからファイルをダウンロードしてください。
補足図2. 重度歯周炎と脳小血管疾患(CSVD)の高負荷との関連に関するサブグループ解析。事前定義されたサブグループにおける重度歯周炎と高CSVD負荷との関連を示すオッズ比(ORs)および95%信頼区間(CIs)のフォレストプロット。垂直の点線はORが1.0であることを示す。こちらをクリックしてファイルをダウンロードしてください。
補足表1. 本研究で行われた解析を裏付ける、匿名化された参加者レベルのデータセット。このデータセットには、統計解析に使用された、匿名化された人口統計学的特性、血管リスク要因、歯周検査変数、臨床検査測定値、磁気共鳴画像(MRI)所見、脳小血管病(CSVD)負荷スコア、および評価者間の一致度変数が含まれています。変数の定義は、「プロトコール」セクションで述べられているものと一致しています。こちらをクリックしてファイルをダウンロードしてください。
本研究では、歯周組織の破壊、残存歯数、およびCSVDの累積負荷の関連性を検討し、最終的に探索的な臨床ノモグラムを開発した。単変量解析では、重度歯周炎、歯の喪失、および高いCSVD負荷の間に相関が認められたが、多変量ロジスティック回帰分析の結果、年齢や高血圧などの従来の心血管リスク因子で調整後、重度歯周炎(P = 0.580)および残存歯数(P = 0.112)は、高いCSVD負荷と独立した関連を示さなかった。したがって、本コホートにおいて、これらの歯周パラメータは高いCSVD負荷の独立した予測因子として解釈されるべきではない。むしろ、観察された関連性は多変量調整後に減弱した。口腔健康と脳血管疾患の関連を検討した先行研究の多くは、ラクナ梗塞や白質高信号域の負荷などの個別の神経画像マーカーに焦点を当てていた26。対照的に、本研究では、現代の神経画像診断の推奨事項および脳小血管疾患に関する現在の概念的枠組みに従い、総CSVD負荷スコアを主要アウトカムとして採用し、びまん性脳微小血管損傷のより包括的な評価を行った1,3,27。
RCS解析による注目すべき知見は、統計的に有意な非線形関連が認められなかったこと(非線形性のP = 0.332)であり、これは、残存歯数の減少と高負荷のCSVD確率の上昇との間に、明確な閾値効果ではなく連続的な線形関係があることを支持している。層別化に使用した残存歯数20本という閾値は、スプライン解析から導出されたものではなく、機能的歯列および成功した口腔老化に関する確立された臨床基準に基づいて選択された16,17,18。極めて重要な点として、二値層別化に用いた残存歯数20本という閾値は、データ解析前に事前指定されていた。このカットオフ値は、現在のデータセットから導出されたものではなく、基本的な口腔機能を維持するために一般的に必要とされる最小歯数を表すために、機能的歯列と成功した口腔老化に関する確立された臨床的コンセンサス(例:WHOの「8020」イニシアチブ)に基づいて選択された。歯の喪失と慢性歯周病は、口腔機能の低下および全身健康への関連があるため、依然として重要な公衆衛生上の懸念事項である28。また、新たなエビデンスは、口腔-脳軸における潜在的な双方向性の関係を示唆しており、脳小血管病自体が歯周組織の悪化を加速させることに関連している可能性がある29。予測因子の選択を最適化するため、座標降下最適化を用いたLASSO回帰モデルを実装した30。従来のステップワイズ回帰の手法と比較して、LASSOは正則化フレームワークを提供し、相関のある臨床変数を同時に検討する場合にモデルの安定性を向上させ、過学習を抑制することができる31。LASSO選択後に保持された変数は、その後、ノモグラム作成のための多変数ロジスティック回帰モデルに投入された。Log_SIIは自動的なLASSO変数選択手順では保持されなかったが、脳血管疾患および歯周炎における全身性炎症のマーカーとしての確立された臨床的関連性があるため、ベースラインおよび拡張予測モデルの両方に事前指定された臨床共変量として意図的に組み込まれた。このアプローチにより、歯周病変数の増分的な寄与の評価を可能にしつつ、全身性炎症状態の調整を確実にした。LASSOは因果推論の手法ではなく予測的な変数選択手法であるため、個々の変数の組み入れまたは除外を、多重共線性や因果的重要性の存在または不在の根拠として解釈すべきではない。
生物学的に、重度歯周炎は局所的なディスバイオーシスを伴う感染症として定義され、全身性の炎症反応に寄与する可能性がある32。歯周治療は、全身性の炎症マーカーを減少させ、末梢血管内皮機能を改善させることが示されており、口腔健康と血管機能との間の関連性を支持している33。Porphyromonas gingivalisを含む歯周病原菌およびそのリポ多糖成分は、咀嚼や口腔衛生処置の際に全身循環に入ることがあり、細菌成分が神経組織内で検出されている34。これらの炎症刺激に対する慢性的な全身曝露は、血管内皮機能不全および血液脳関門の障害に寄与する一つのメカニズムとして提案されており、これらのプロセスは脳小血管疾患の病態生理に関与している35。得られた予測モデルの臨床的有用性を評価するため、さまざまな閾値確率における潜在的な正味の利益を評価する目的でDCAを実施した36。DCAはデフォルトの戦略に対する潜在的な優位性を示唆したが、これは理論的なモデル性能を表しており、現実世界での臨床的利益や患者のアウトカムの改善を確立するものではない。透明性のある報告を促進し、独立した評価と再現を容易にするため、この探索的なワークフローの開発と報告は、個人の予後または診断のための多変数予測モデルの透明な報告(TRIPOD)ステートメントに従った37。慢性の低グレードな口腔内炎症が、脳の神経変性に関与する全身的な代謝経路と相互作用するという仮説が立てられているため、この報告フレームワークは適切である38。モデル評価には、確立された検証尺度に加え、予測性能の追加指標を組み込んだ39。より広範には、このアプローチは局所的な粘膜炎症が全身に及ぼす潜在的な影響を反映している40。それにもかかわらず、本モデルの識別能が中程度であること(AUC = 0.688)、DeLong検定によるベースラインモデルと比較したAUCの統計的に有意な改善が見られないこと(P = 0.263)、および口腔健康変数を取り巻く信頼区間が広いことから、提案されたノモグラムは、日常的な神経科診療における臨床的に実行可能なツールではなく、あくまで予備的な探索的視覚化フレームワークとして厳格に解釈されるべきである。
提案されたワークフローの再現性と正常な実施には、いくつかの方法論的なステップが重要であり、それらが報告されたモデル性能に影響を与える可能性があります。第一に、歯周組織検査は高度に標準化されている必要があり、校正を受けた検査者が、2018年EFP/AAP分類フレームワーク15に従って、1歯あたり6部位で全顎的なプロービングポケット深さと臨床的付着レベルの測定を行う必要があります。第二に、神経画像診断の再現性は、FLAIR(流体減衰反転回復)およびSWI(感受性強調画像)シーケンスを含む標準化された3.0 T磁気共鳴画像法(MRI)の撮像プロトコルと、全脳小血管病(CSVD)負荷量に関する検証済みの合意に基づくスコアリング基準に依存します1。第三に、全身性免疫炎症指数(SII)の自然対数(Log_SII)を信頼性高く算出するため、空腹時の静脈血検体を用い、一貫した検体処理手順を用いて採血および臨床検査処理を標準化する必要があります。さらに、手入力によるデータエントリーミスを減らすための電子歯科チャートの導入や、欠損データの管理に関する標準化された手順などの対策が、実施精度を高める可能性があります。
これらの方法論的な検討事項はあるものの、いくつかの限界を認める必要がある。単一施設でのレトロスペクティブな横断的研究デザインであるため、因果関係に関する結論を出すことはできない。さらに、高CSVD負荷のイベント129件を含めたことで、安定したモデル構築のための変数あたりイベント数の基準は満たされたが、サンプルサイズが比較的小さかったため、多変量解析の統計学的パワーが制限された。また、本モデルでは、歯磨きやフロスの頻度などの口腔衛生習慣や、過去の歯周治療歴など、歯周状態に影響を及ぼし得るいくつかの潜在的な交絡因子が考慮されていない。加えて、進行したCSVD患者や脳卒中後の後遺症がある患者は、認知機能、運動機能、および手指の巧緻性が低下している可能性があり、それが口腔衛生習慣を損ない41、定期的な歯科ケアの利用を減少させている42可能性があるため、逆の因果関係を排除できない。したがって、進行した歯周炎は、神経血管障害の原因ではなく、一部にはその結果である可能性がある。最後に、レトロスペクティブなデザインであったため、16SリボソームRNAシーケンシングやメタゲノムプロファイリングのための歯下プラークサンプルの収集ができず、口腔マイクロバイオームと血中炎症性メディエーターまたはアミロイド関連バイオマーカーとの潜在的な関係の調査が制限された43。独立した外部検証コホートが存在しないことは依然として大きな限界であり、提案したモデルの他の集団への汎用性はまだ確立されていない。今後は、正式なサンプルサイズ設計、外部検証、および追加の生物学的マーカーの組み込みを行った、多施設共同プロスペクティブ研究が必要である。
結論として、重度の歯周炎および残存歯数が20本未満であることは、高いCSVD負荷量と有意な単変量相関を示しましたが、多変量調整後にはこれらの相関は減弱し、本コホートにおいて高いCSVD負荷量と独立して相関してはいませんでした。歯周パラメータと従来のリスク要因を統合した提案済みの臨床ノモグラムは、臨床的に検証されたスクリーニングツールではなく、個別化されたリスク可視化のための予備的な非侵襲的探索的枠組みを提供するものです。日常的な臨床応用の検討には、独立したコホートでの外部妥当性検証が必要です。口腔健康とCSVDの間の因果関係を明らかにするには縦断的研究が必要ですが、歯周健康の維持と機能的な歯列の保存は、加齢に伴う脳血管疾患に関する今後の多角的共同研究における重要な領域であり続けます。
利益相反:
著者らは、本研究に関して利益相反がないことを宣言します。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| 自動血球計数装置 | Sysmex | XN-1000 | 全身性免疫炎症指数(SII)の算出に必要な末梢血細胞数を測定するために使用。 |
| EDTA採血管 (K2 EDTA) | Becton, Dickinson and Company | 367841 | 空腹時静脈血採血に使用。 |
| 電子カルテ(EMR)システム | Winning Health Technology Group | Winning EMR v6.0 | 患者の属性および臨床情報を抽出するために使用。 |
| 磁気共鳴画像法(MRI)スキャナー (3.0 T) | Siemens Healthineers | MAGNETOM Prisma 3.0T | 標準的な頭部MRI検査に使用。 |
| MRI頭部コイル | Siemens Healthineers | 64-channel Head/Neck Coil | MRI撮像に使用される標準的な頭部コイル。 |
| 歯周プローブ (UNC-15) | Hu-Friedy | PCPUNC15 | 歯周ポケット深さ(PPD)および臨床的アタッチメントレベル(CAL)の測定に使用。 |
| Rパッケージ: dcurves | CRAN | Version 0.4.0 | 決定曲線分析(DCA)に使用。 |
| Rパッケージ: glmnet | CRAN | Version 4.1-8 | LASSO(least absolute shrinkage and selection operator)回帰および変数選択に使用。 |
| Rパッケージ: pROC | CRAN | Version 1.18.5 | 受信者動作特性(ROC)曲線分析およびDeLongテストに使用。 |
| Rパッケージ: PredictABEL | CRAN | Version 1.2-4 | 連続ネット再分類改善(NRI)および統合識別改善(IDI)の算出に使用。 |
| Rパッケージ: rms | CRAN | Version 6.7-1 | 制限付き立方スプライン(RCS)分析、ノモグラム作成、およびブートストラップ校正に使用。 |
| 統計計算用Rソフトウェア | R Foundation for Statistical Computing | Version 4.3.1 | すべての分析に使用した統計計算環境。 |