Regulatory Evolution of TEM β-lactamases through Large-Scale Genomic Analysis Reveals Promoter Diversity and its Impact on β-lactam Resistance in Escherichia coli

技術要約：大規模ゲノム解析を通じたTEM β-ラクタマーゼの調節進化

問題提起
大腸菌（Escherichia coli）における薬剤耐性（AMR）は、主にβ-ラクタマーゼ遺伝子、特にbla^TEMファミリーの獲得と発現によって駆動されている。bla^TEMのコーディング配列の多様化は十分に文書化されているが、上流の調節変異（具体的にはプロモーターの多様性）がどの程度耐性表現型に寄ari寄与しているかは依然として不明である。先行研究では、様々な転写強度を持つ標準的なプロモーター型（P3、Pa/Pb、P4、P5）が特定されているが、大規模なゲノムデータセットにおけるそれらの普及率、進化関係、および耐性への因果的影響を体系的に調査した研究は不足していた。さらに、非標準的なプロモーターモチーフの役割と、ハウスキーピング・シグマ因子σ⁷⁰との構造的適合性についても、さらなる解明が必要である。

手法
本研究では、大規模ゲノム解析、因果推論、および構造生物学を組み合わせた統合的なフレームワークを採用した：

1. ゲノムデータの収集と処理：

   • データセット： BV-BRC（Bacterial and Viral Bioinformatics Resource Center）から60,087個の大腸菌ゲノムを取得し、そのうち明示的にアノテーションされたプラスミド配列を含む2,548個の分離株を対象とした。
   • 遺伝子の特定： ABRicateを用いてCARDデータベースに対してbla^TEM遺伝子をスクリーニングした（同一性およびカバレッジの閾値は90%）。これにより、染色体に関連する20,317個およびプラスミドに関連する965個の分離株から、それぞれ21,367個および1,242個のbla^TEM配列が得られた。
   • 上流解析： カスタムPythonスクリプトを用いて上流の調節領域（最大500 bp）を抽出し、-35および-10プロモーターモチーフを特定した。1-kbの隣接領域内の挿入配列（IS）はISEScanを用いて特性評価を行った。
   • 表現型のキュレーション： 17種類のβ-ラクタム系抗菌薬に関する耐性表現型をメタデータからキュレーションし、中間的な耐性を耐性と再分類した。

2. 因果推論フレームワーク：

   • 目的： 基底となるbla^TEM変異プロファイルを一定に保ったまま、プロモーター遷移が耐性確率に与える平均処置効果（ATE）を推定すること。
   • モデル： 反事実的回帰アプローチを用いた。プロモーター型を処置変数、抗菌薬耐性を二値の出力、特定のbla^TEM変異の有無を調整共変量とした。
   • 検証： 多項プロペンシティスコアを推定し、ポジティビティ（重複）仮定を評価した。ATEはロジスティック回帰を通じて計算され、安定性（一貫性≥80%）を確保するための1,000回のブートストラップ反復と、統計的有意性（p ≤ 0.05）を決定するための500回の置換テストが行われた。

3. 構造および分子動力学（MD）解析：

   • 選択： 他のβ-ラクタマーゼによる交絡を最小限に抑えるため、単一のbla^TEMを持つ分離株から特定された新規プロモーターモチーフを、その出現頻度に基づいて優先的に選択した。
   • シミュレーション： 選択されたモチーフを、大腸菌RNAポリメラーゼσ⁷⁰オープン複合体（PDB: 6CA0）の文脈内でモデリングした。Desmondを用い、OPLS4力場を使用して100 nsのMDシミュレーションを実施した。
   • 相互作用解析： タンパク質-DNA相互作用（水素結合、塩橋、π-カチオン、π–πスタッキング）を、軌跡の占有率10%以上を閾値として定量化した。

主な結果

• ゲノム分布と多様性：

  • bla^TEM-1およびbla^TEM-150は、染色体およびプラスミドの両方のコンテキストにおいて最も一般的な変異であった。
  • シャノン・エントロピー解析により、上流の調節領域（−250から−50 bp）は、高度に保存されたコーディング領域と比較して、有意に高い配列多様性を示すことが明らかになった。プラスミドに関連する上流配列は、染色体に関連するもの（〜0.55 bits）よりも高い多様性（〜0 ways 0.75 bits）を示した。
  • P3プロモーターは両方のコンテキストにおいて支配的な標準型であったが、その次にPa/Pbが続いた。しかし、多数の新規な-10/-35の組み合わせが検出され、これらは主に染色体配列に見られた。

• 表現型とプロモーターの関連性：

  • bla^TEM遺伝子は耐性株と感受性株の両方で見出されたことは、遺伝子の存在自体が耐性を決定づけるわけではないことを示している。
  • 耐性株は、感受性株と比較して、より強力な標準的プロモーター（Pa/Pb、P4、P5）や新規モチーフを含む、より大きな調節多様性を示した。
  • bla^TEM-1のような祖先変異は、コーディング配列の変化なしに、異なるプロモーター型間を移行する「マルチプロモーター利用」を示した。

• 因果推論の結果：

  • 反事実的解析により、標準的なP3プロモーターから他の標準的（Pa/Pb、P4、P5）および新規プロモーター型への遷移は、複数の抗菌薬に対して正のATE（0.16から0.46の範囲）をもたらすことが示された。
  • これは、bla^TEM変異とは独立して、P3からこれらの代替プロモーターへの切り替えが、予測される耐性確率を最大46パーセントポイント上昇させることを示している。

• 新規モチーフの構造的検証：

  • 12個の耐性・単一bla^TEM分離株から、5つの新規プロモーターモチーフが特定された。
  • MDシミュレーションにより、選択された新規モチーフ（特にM1、M2、M4）がσ⁷⁰ドメインとの安定した相互作用を維持していることが確認された。
  • 特筆すべきは、新規-35モチーフである**M2 (TGCAAC)**が、25個の保持された相互作用（13個のH結合、12個の塩橋）を示し、標準的なモチーフと同等またはそれを上回る性能を持ち、高発現のP5様プロモーターと重要な相互作用残基を共有していたことである。これは、機能的な調節要素としての構造的適合性を支持している。

• 挿入配列（IS）のコンテキスト：

  • IS要素は、染色体に関連するbla^TEM遺伝子の約14.5%（20,317個中2,948個）およびプラスミドに関連するものの26.1%（965個中252個）で特定された。
  • 特定されたIS要素のうち、染色体に随伴するものの30%が「完全（complete）」に分類されたのに対し、プラスミドに関連する要素の73%が完全であった。
  • IS要素（例：IS91、IS1182、IS6）はbla^TEMとしばしば関連していたが、ほとんどの上流重複は遺伝子開始点から400 bp以上離れた位置にあり、近位のプロモーターモチーフを直接破壊することはないことが示唆されたが、構造には影響を与える可能性がある。

意義と主張
本論文は、本研究がbla^TEMの調節進化を解明するために、ゲノミクス、因果推論、および構造解析を大規模に統合した初めての試みであることを主張している。主な意義は以下の点にある：

1. 調節の多様化が主要な駆動力である： プロモーターの変異は、コーディング配列の変異とは異なる、β-ラクタム耐性を高めるための主要な進化経路である。
2. 新規モチーフは機能的である： 本研究は、標準的なプロモーターのみが高レベルの発現を駆動するという見解に疑問を投げかけ、機能的な耐性レギュレーターとして機能する可能性が高い新規プロモーターモチーフ（特にM2）を特定し、構造的に検証した。
3. 手法の有用性： 反事実的な因果フレームワークは、遺伝的背景からプロモーター型の効果を分離することに成功し、耐性表現体を単なる遺伝子の保有ではなく、調節の変化に帰属させるための堅牢な手法を提供した。
4. コンテキストの違い： 本研究は、染色体とプラスミドのコンテキスト間における異なる進化圧を浮き彫りにした。プラスミドは、その移動性とコピー数効果により、より高い上流の多様性を示している。

著者らは、耐性の表現型は遺伝子の存在そのものだけでなく、耐性遺伝子、変異、および調節要素の間の複雑な相互作用によって形成されることを示唆しており、これらの調節メカニズムを理解することがAMR進化の完全な理解には不可欠であると結論付けている。