Germline genomic and methylomic dynamics following three generations of early-life metabolic challenges

技術要約：3世代にわたる早期生命における代謝的挑戦に伴う生殖細胞系列のゲノムおよびメチロームのダイナミクス

問題提起
環境および食事因子、特に過栄養のような早期生命における代謝的挑戦は、健康と発達に多世代的な影響を及ぼすことが知られている。これまでの研究により、代謝異常が（miRNAやDNAメチル化などの）エピジェネティックなメカニズムを通じて父系生殖細胞系列を通じて伝達されることは確立されているが、多世代的（連続的な）曝露による長期的なゲノムへの影響については、依然として十分に理解されていない。具体的には、3世代にわたる持続的な早期生命の代謝ストレスが、生殖細胞ゲノム（SNPおよびCNV）にどのように影響を与えるのか、また、それらのゲノム変化がエピゲノムの変化（DNAメチル化）によって引き起こされるのか、あるいは相関しているのかは不明である。既存研究の多くは、単発の開発期曝露（世代間伝達）に焦点を当てており、環境性肥満誘発物質（オベソゲン）がより現実的なシナリオを示す「多世代的」な連続曝露については扱っていない。

方法論
本研究では、リターサイズ（産仔数）の減少によって早期生命の代謝的挑戦をシミュレートするために、アウトブレッド（外来交雑）系統であるICR (CD1) マウスモデルを用いた。

• 実験デザイン: 3世代（F0–F2）を追跡した。過栄養（ON）群は、母体あたりの産仔数を4匹に減少させたリターで構成され、対照（CT）群は8匹とした。このモデルは、新生児期の肥満およびその後の代謝症候群を誘発するために確立されている。各世代のオスは、父系血統の影響を分離するために、外部のナイーブな雌と交配させた。
• 表現型解析: 生後4ヶ月時点での生理学的および代謝的パラメータ（体重、グルコース、インスリン、肝臓および脂肪組織重量）を評価した。
• シーケンシング戦略: 精子DNAを抽出し、ゲノム変異の解析にはGenotyping-by-Sequencing (GBS) を、メチローム・プロファイリングにはGBSとMethyl-Immunoprecipitation (GBS-MeDIP) を併用して並行解析を行った。
• ゲノム解析:
  • SNP: GATKのベストプラクティスを用いて同定した。主成分分析（PCA）および階層的クラスタリング（Identity by State）を用い、遺伝的類似性と処理の影響を評価した。
  • CNV: GATKの GermlineCNVCaller を用いて生殖細胞系列のコピー数多型（CNV）をコールした。偽陽性を排除するため、ランダムなサブサンプリング制御を含む厳格なフィルタリングプロセスを適用した。
  • アノテーション: 置換配列（TEs）および反復配列（REs）に対するCNVのアノテーションには、濃縮を決定するための置換テストを用いた。
• エピゲノム解析:
  • DMR: 世代内（CT vs. ON）および世代間（F1 vs. F2）の比較において、Mann-Whitney検定を用いて差次的メチル化領域（DMR）を特定した。
  • パスウェイ解析: メチル化されたウィンドウに関連する遺伝子について、KEGGおよびPathBIXを用いたパスウェイ濃縮解析を行った。
• 統合: 父系のメチル化レベルと、次世代における新規SNPまたはCNVの出現との関係を調査するために、Spearmanの順位相関を算出した。

主な結果

1. ゲノムへの影響 (SNP): SNPデータのPCAにより、無関係なONファミリーがCTファミリー（血縁関係に基づきクラスター化）とは異なる形で共にクラスター化していることが明らかになった。これは、代謝的挑戦が共通のゲノムシグネチャーを誘導したことを示唆している。処理による遺伝的分散の説明率は約4.6%であった。さらに、ON群はCT群と比較して、次世代における新規SNPの出現が「減少」しており、これは遺伝的多様性の制約を示唆している。
2. ゲノム不安定性 (CNV): ON群では、特定のCNVイベントの出現が見られた（F1からF2まで維持された7つの推定イベント）。これらのイベントは、長鎖インタースパーシング核要素（LINEs）および長末端反復配列（LTRs）に有意に濃縮されていた。欠失は世代を超えて安定していたが、重複は拡大する傾向があった。
3. エピゲノムへの影響: メチロームのグローバルなPCAでは、ON群とCT群の間に明確な分離は見られず、広範なグローバルなメチル化シフトは存在しないことが示された。しかし、パスウェイ濃縮解析の結果、ON群はCT群と比較して、多糖類代謝、リポリシス（脂肪分解）、P53、mTOR、およびヘッジホッグシグナルリングに関連するパスウェイにおいて、濃縮シグナルが弱い、あるいは欠如していることが判明した。
4. 領域的メチル化変化: DMR解析により、F2世代における特異的な変化が特定された。これには、LINEsおよびSINEsにおける低メチル化、ならびに Fhod3 および Etl4 遺伝子における高メチル化が含まれる。世代間比較においても、Fhod3 およびガンマ・サテライトリピートにおけるDMRが強調された。
5. ゲノム-エピゲノムのダイナミクス:
   • SNP: CT群では、父系のメチル化レベルが高いほど、次世代におけるSNPの出現率が高かった（メチル化されたシトシンが変異しやすいという既知の性質と一致する）。対照的に、ON群では逆の相関が観察された。これは、代謝ストレス下では、酸化損傷に対する保護メカニズムにより、メチル化されたCpGが変異に対してより抵抗性を持つ可能性を示唆している。
   • CNV: CNVにおけるLINEs/LTRsの濃縮と、これらの要素における差次的メチル化は、反復配列（RE）のメチル化異常がCNV形成を駆動している可能性を示唆している。

意義および主張
著者らは、本研究が、多世代的な代謝的挑戦が、生殖細胞におけるゲノムおよびエピゲノムの両方の変化を誘発し得るという証拠を提供すると主張している。主な主張は以下の通りである：

• 制約 vs 不安定性: 代謝的挑戦は、SNPレベルでは遺伝的多様性を制約する一方で、構造レベル（CNV）ではゲノムの不安定性を誘発するように見える。
• メカニズム: 観察されたゲノム不安定性（CNV）は、広範なグローバルDNAメチル化の変化ではなく、転移因子（特にLINEsおよびLTRs）の調節不全によって媒介されている可能性が高い。
• 現実性: 3世代にわたる連続的な曝露を追跡することで、本研究は、単発曝露の世代間伝達モデルと比較して、環境性肥満誘発物質が進化のポテンシャルやゲノム構造にどのように影響を与えるかについて、より現実的な描写を提供している。
• 新規性: 本研究は、多世代曝露シナリオにおいて、精子のDNAメチル化の変化がゲノムにどのように影響するかを扱った初の制御された研究であり、エピジェネティックな摂動がどのように構造的なゲノム変異へとつながるかという複雑な相互作用を浮き彫りにしている。

結論として、早期生命における代謝的挑戦は、ゲノム構造に対して長期的な影響を及ぼし、それは転移因子の活動によって媒介される可能性があり、ゲノムとエピゲノムのダイナミクスを現実的な曝露シナリオの下で共に研究することの必要性を強調している。