ペル・ポリフルオロアルキル物質(PFAS)は、フライパンの焦げ付き防止加工や防水・防汚加工など幅広い工業製品・日用品に使われてきた化学物質群です。化学的・熱的に極めて安定しているため環境中でほとんど分解されず、「フォーエバー・ケミカル(永遠の化学物質)」とも呼ばれています。飲料水や食品、大気中の粉塵などを介して体内に取り込まれ、肝臓や腎臓、血液のほか尿・羊水・卵胞液・胎盤組織からも検出されており、胎盤を通過し母乳にも移行することが報告されています。神経毒性、免疫毒性、発がん性、内分泌かく乱作用、生殖毒性など複数の臓器にまたがる健康影響が指摘されており、代表的な物質であるペルフルオロオクタン酸(PFOA)とペルフルオロオクタンスルホン酸(PFOS)は、ストックホルム条約の残留性有機汚染物質として国際的に規制が進んでいます。
卵巣予備能低下(DOR)は、卵子の数や質が年齢相応の水準を下回った状態を指し、女性不妊の主要な原因の一つです。診断には、抗ミュラー管ホルモン(AMH)が1.1 ng/mL未満、胞状卵胞数(AFC)が7個未満、基礎卵胞刺激ホルモン(FSH)が10 IU/L以上、といった基準が用いられます。原因は遺伝的素因や自己免疫疾患、感染症、喫煙、薬剤使用など多因子にわたるとされ、近年は環境汚染物質への曝露も一因として注目されています。PFASの曝露が多囊胞性卵巣症候群や子宮内膜症、早発卵巣不全、DORのリスク上昇と関連することは複数の疫学研究で示されてきましたが、PFASが顆粒膜細胞のどの分子を狂わせてDORにつながるのかという経路は、十分に解明されていませんでした。
この問いに、ネットワーク毒性学・機械学習・分子ドッキング・分子動力学(MD)シミュレーションという計算科学的手法と、KGN細胞(ヒト卵巣顆粒膜細胞に由来)を用いた培養実験を組み合わせて挑んだ研究が、Scientific Reports誌(2026年、16巻:26191、doi:10.1038/s41598-026-56614-2)に発表されました。代表的な4種のPFAS(PFOS、PFOA、PFHxS、PFNA)を対象にした網羅的な標的探索と、そのうち最も曝露量の多いPFOAを用いた実証実験から、PFASが卵巣予備能を削る道筋が描き出されています。
方法──計算科学と培養細胞実験を組み合わせた二段構え
研究は大きく二段階で構成されています。第一段階では、PFOS、PFOA、PFHxS、PFNAの4化合物についてChEMBL、STITCH、SwissTargetPrediction、PharmMapperから毒性標的を収集し、GeneCardsとOMIMから集めたDOR関連標的932個と突き合わせて候補を絞り込みました。得られた候補でPPI(タンパク質間相互作用)ネットワークを構築し、GEOデータベース(GSE232306)のDOR患者6例と健常対照6例の発現データとも突き合わせたうえで、ランダムフォレスト・LASSO回帰・XGBoostという3種の機械学習を段階的に適用し、核心となる標的遺伝子を選び出しています。第二段階では、標的タンパク質とPFASとの結合を分子ドッキングと100ナノ秒間のMDシミュレーションで検証したうえで、KGN細胞(ヒト卵巣顆粒膜細胞由来)にPFOAを曝露させ、増殖・アポトーシス・酸化ストレス・ホルモン分泌への影響をCCK-8、EdU、フローサイトメトリー、RT-qPCR、ウエスタンブロット、ELISAなどで測定しています。
35個の候補から、たった3つの遺伝子へ
4種のPFASに共通する毒性標的は435個に絞り込まれ、DOR関連標的932個と照合した結果、両者に重なる候補標的は35個でした。この35個のPPIネットワークは33ノード・140エッジで構成され、平均次数(一つの分子がつながる相手の数)は8でした。次数の高い上位10標的には、アルブミン(ALB)やエストロゲン受容体1(ESR1)、PPARγ、ERBB2などが含まれ、GO解析では482件の生物学的プロセス、25件の細胞構成要素、81件の分子機能が、KEGG経路解析では72件の経路が有意に抽出されました。関与が目立ったのは、PI3K-Akt、MAPK、Rasなど細胞の増殖・分化に関わる経路、Th17細胞分化などの免疫・炎症関連経路、プロラクチンやエストロゲンなどのホルモンシグナル伝達経路でした。
続いて、DOR患者6例と健常対照6例の顆粒膜細胞の発現データを解析したところ、5,836個の発現変動遺伝子(上昇2,593個、低下3,243個)が同定され、先の35標的と重なる遺伝子は19個に絞られました。この19個にランダムフォレスト(11個)、LASSO回帰(4個)、XGBoostという3段階の機械学習を順に適用した結果、最終的にグルタチオンS-転移酵素P1(GSTP1)、ペルオキシレドキシン5(PRDX5)、TGF-β受容体1型(TGFBR1)という3つの核心遺伝子が残りました。DOR患者の顆粒膜細胞では、健常対照と比べてGSTP1とPRDX5の発現が有意に低下し(P<0.01)、TGFBR1の発現は有意に上昇していました(P<0.01)。これら3遺伝子は、発見データセットでの診断精度(AUC。1に近いほど精度が高い指標)が0.9を超え、2つの独立した外部データセットでも、GSE315266でAUC 0.75(GSTP1)・0.78(PRDX5)・0.75(TGFBR1)、GSE296744でAUC 0.639(GSTP1)・0.667(PRDX5)・0.819(TGFBR1)と一定以上の精度を保っていました。特にTGFBR1はAUC 0.819と高く、コホートを超えて再現性のある指標であることが示唆されています。
結合エネルギー最大-9.8 kcal/mol――PFASは標的に強くくっついていた
3つの核心遺伝子が実際にPFASと結合しうるかを確かめるため、GSTP1、PRDX5、TGFBR1それぞれのタンパク質立体構造と4種のPFASとの間で分子ドッキングが行われました。結合エネルギーは低いほど(マイナス側に大きいほど)安定な結合を示す指標で、一般に-5 kcal/mol未満であれば安定した結合とみなされます。今回の解析では、3標的と4種のPFASの組み合わせすべてで結合エネルギーが-9.8~-6.5 kcal/mol(平均-7.9 kcal/mol)となり、いずれも安定した結合を示しました。なかでもTGFBR1は4種すべてに対して-8.0 kcal/molを下回る、最も強い結合を示しています。さらに代表的なPFASであるPFOAを用いた100ナノ秒間の分子動力学(MD)シミュレーションでも、GSTP1・PRDX5・TGFBR1の各PFOA複合体は10~65ナノ秒ほどで構造が安定し、水素結合もGSTP1で0~3本、PRDX5で0~5本、TGFBR1で0~6本と安定的に形成されました。あくまで計算科学的な予測ですが、PFOAが3つの標的タンパク質に持続的に結合しうることを支持しています。
培養顆粒膜細胞で確かめられたこと
予測された機序が実際の細胞でも起きているかを調べるため、KGN細胞にPFOAを曝露する実験が行われました。0.2~200 µMのPFOAを24~48時間作用させたところ、濃度と曝露時間に依存して生存率と増殖能(EdU取り込み)が低下し、200 µMでは顕著な増殖抑制が確認されました。フローサイトメトリーでは2~20 µMの曝露でアポトーシス(プログラムされた細胞死)を起こした細胞の割合が有意に増加し(P<0.001)、これを促進するBAXタンパク質の発現上昇と、抑制するBcl-2の発現低下が確認されました。
ホルモン分泌への影響としては、48時間の曝露で高濃度のPFOAが、卵巣予備能の指標であるAMHと、インヒビンBを構成するINHBA・INHBBのmRNA発現、およびAMHタンパク質の分泌量を有意に低下させました(いずれもP<0.05)。あわせて、先の計算科学的解析で見出された核心遺伝子についても、48時間の高濃度曝露でGSTP1とPRDX5の発現が有意に低下し、TGFBR1の発現が有意に上昇しており(いずれもP<0.05)、予測と実験結果が一致しました。
100 nMという、臨床検体に近い濃度でも起きていたこと
これらの結果は比較的高い濃度での曝露に基づくものですが、この研究では10~1,000 nMという、より環境曝露に近い濃度での72時間曝露も検討されています。この濃度域では生存率は92~94%程度に保たれていたものの、統計学的には有意な低下が認められました。さらに、100 nM以上で活性酸素種(ROS。酸化ストレスの原因となる分子)が有意に増加し、AMHの分泌量も100 nMおよび1,000 nMで有意に低下していました(10 nMでは変化なし)。100 nMは実際に一部の卵胞液や血清で検出される濃度に近く、細胞が死に至るよりも前の段階で、卵巣予備能に関わるホルモン分泌機能がすでに損なわれている可能性を示す結果です。
臨床的意味と日本への含意
日本では、卵巣予備能の評価にAMHやAFC、FSHが日常的に用いられており、体外受精を検討する女性にとって卵巣予備能の見立ては治療方針を左右する重要な情報です。今回の研究が示したGSTP1・PRDX5・TGFBR1は、まだ臨床応用が確立した検査項目ではありませんが、環境曝露と卵巣予備能低下を分子レベルで結ぶ手がかりといえます。卵胞液中のPFOA濃度は中国のコホートで中央値5.56~5.85 ng/mL、オーストラリア・欧州・米国では平均1.80~2.40 ng/mL程度と、地域差も報告されています。日本人女性の詳細なデータは今回の対象外ですが、診察室で不妊の原因を説明する医師にとっては、遺伝的要因や年齢要因に加え、環境曝露という視点が一つの補助線になりうる知見です。今回示された関連はあくまで確率としてのリスクであり、PFASに曝露したからといって必ずDORになるわけではありません。
研究の強みと限界
この研究の強みは、計算科学的な標的探索にとどまらず、機械学習による絞り込みと2つの外部データセットによる検証、さらに培養細胞を用いた実証実験までを一貫して行った点にあります。既存研究の多くは、代謝物との相関解析や一経路の検討にとどまるか、計算科学的な予測にとどまり細胞実験による裏づけを欠くものでした。今回の研究は、標的の同定から細胞レベルでの機能確認までを結びつけている点で、これらを補っています。
一方で限界も少なくありません。発見に用いたデータセットの症例数はDOR群・対照群とも6例ずつと小規模であり、同定された遺伝子は確定した診断バイオマーカーというより機序解明のための候補と位置づけるべきです。分子ドッキングとMDシミュレーションはあくまで計算科学的な予測であり、実際の生化学的な結合を証明するものではありません。培養細胞実験で用いられたのはPFOA一種類のみで、他のPFASにそのまま当てはめることはできません。KGN細胞は平面(2次元)培養であり、卵巣組織が持つ立体構造や線維化のような組織レベルの変化は再現できておらず、TGFBR1の上昇もコラーゲン沈着などの直接的な線維化指標は確認されていません。さらに、遺伝子の働きを人為的に操作してDORが実際に起きるかを確かめる機能的な実験は行われておらず、原因であることを直接証明するには至っていません。
おわりに
今回の研究は、ネットワーク毒性学・機械学習・分子ドッキング・MDシミュレーションという計算科学的アプローチと、培養細胞を用いた実証実験を組み合わせ、PFAS、とりわけPFOAが卵巣予備能を低下させうる分子的な道筋の一端を示しました。GSTP1とPRDX5の低下による酸化ストレスへの脆弱性の増大、TGFBR1の上昇による線維化シグナルの活性化、BAX/Bcl-2バランスの変化によるアポトーシスの促進、さらにAMHとインヒビンBの分泌低下という複数の変化が、環境的に現実的な低濃度でも生じうることが示されています。今後は、より大規模な臨床データセットでの検証や機能的な実験、in vivoモデルでの組織学的な確認を通じて、臨床応用につながるかどうかが問われることになります。
参考文献
1. Li, H., Hu, X., Guo, Z. & Sun, Z.『Integrated computational and in vitro analysis reveals the toxic mechanisms of PFAS in female diminished ovarian reserve』Scientific Reports 2026;16(1):26191. doi:10.1038/s41598-026-56614-2