Karyotypes of the Anopheles 2La chromosomal inversion carry functionally distinct transcriptional enhancers

技術要約：Anopheles 2La 染色体逆転が機能的に異なる転写エンハンサーを保持していることについて

問題提起
染色体逆転は、共進化したアレル（対立遺伝子）のブロックを捕捉し、組換えを抑制して、明確に異なる表現型に関連する安定したハプロタイプを形成するための、重要な適応メカニズムである。Anopheles gambiae 種複合体において、2La 傍中心逆転（ゲノムの約10%）は、マラリア伝播効率、殺虫剤耐性、および環境適応と強く関連している。しかし、これらの表現型の違いを駆動する分子メカニズムは、依然として大部分が解明されていない。この理解における大きな障壁は、ゲノム研究における「アサーテインメント・バイアス（採録バイアス）」である。標準的な解析は単一のリファレンスゲノム（一方の逆転カリオタイプに対してホモ接合体）に依存しており、そのため、もう一方のカリオタイプに特有、あるいは高度に分岐した調節要素が不可視となってしまう。その結果、逆転カリオタイプに特異的な新規調節要素の進化や、祖先的な要素の喪失をゲノムワイドなスケールで特徴付けることは困難であった。

方法論
著者らは、Anopheles coluzzii における逆転特異的な転写エンハンサーを同定および検証するために、2つの独立した補完的な計算および実験的アプローチを採用した：

1. リファレンス・アセンブリ比較 (STARR-seq マッピング):

   • 既報の STARR-seq データ（30 個の 2La/2La 型雌から得られたもの）を MOPTI リファレンスゲノム（2La/2La）に再マッピングし、元の PEST リファレンスゲノム（2L+a/2L+a）へのマッピングと比較した。
   • MOPTI アセンブリに特有のピークを、2La カリオタイプに特異的な候補エンハンサーとして特定した。
   • アセンブリのアーティファクトを軽減するため、候補をユニークな配列（トランスポゾンや単純反復配列を除く）でフィルタリングし、野生のホモ接合体蚊を用いて PCR およびサンガーシーケンシングによって検証した。

2. 集団再シーケンシング (Ag1000G データ):

   • 2L 染色体を解析し、2La/2La と 2L+a/2L+a のホモ接合体間で差別的に固定されている一塩基多型（SNP）を特定した（定義：一方のホモ接合体で頻度が >95%、他方で <5%）。
   • これにより、18,072 個の「逆転特異的 SNP」が得られた。
   • これらの SNP を 788 個の STARR-seq エンハンサーのカタログに重ね合わせ、逆転特異的な変異の負担が高い（定義：≥6 個の逆転特異的 SNP を持つ）エンハンサーを特定した。

3. 機能検証:

   • 9 つの候補エンハンサー（アセンブリによるアプローチから 1 つ、集団アプローチから 8 つ）を Firefly ルシフェラーゼレポーターベクターにクローニングした。
   • アレルは、細胞株（2La 用の Sua4.0、2L+a 用の Ag55）および野生の蚊から増幅した。
   • 両方の細胞株において、相対的なエンハンサー活性を測定するためにルシフェラーゼアッセイを実施した。
   • 部位特異的変異導入 (SDM): 最も有望な候補について、特定の逆転特異的 SNP および挿入/欠失（indel）を異なる遺伝的背景間でスワップし、それらの変異が活性の違いを説明するために必要、あるいは十分であるかどうかを決定した。
   • マルチオミクス統合: ATAC-seq データを用いてクロマチンアクセシビリティを評価し、Micro-C データ（近接ライゲーション法）を用いて、3D クロマチン相互作用を介した候補標的遺伝子を特定した。

主な結果

• 候補の同定: アセンブリによるアプローチでは 28 個の 2La 特異的候補エンハンサーが同定され、ユニークな配列およびリプリケートによるフィルタリングを経て、3 つが検証用に選択された。集団アプローチでは、逆転特異的 SNP の負担が高い 8 つのエンハンサーが同定された。
• 機能的分離: テストされた 9 つの候補のうち、2 つが明確かつ堅牢な逆転特異的な機能差を示した：
  • ISEN2L_21.9MB: 2La アレルは 2L+a アレルよりも有意に高い活性を示した。変異導入の結果、13 個の逆転特異的 SNP は活性の差を完全に再現するために「必要」であるが、「十分」ではないことが明らかになった。また、8 bp の逆転特異的 indel も同定されたが、単独では活性に有意な影響を与えなかった。
  • ISEN2L_25.8MB: このエンハンサーには、大きな（32–40 bp）逆転特異的 indel が含まれていた。部位特異的変異導入により、この indel がカリオタイプ間のエンハンサー活性の差を説明するために「必要かつ十分」であることが示された。この indel には、予測される転写因子（例：combgap）の結合サイトが含まれている。
• ゲノム分布: 逆転特異的 SNP は遺伝子間領域に濃縮されており、イントロン領域では減少している。これらは 2La 逆転全体にわたって U 字型の分布を示し、ブレイクポイント付近で密度が高く、中央部で減少している（おそらく二重組換えイベントによる）。
• 標的遺伝子: Micro-C 解析により、ISEN2L_25.8MB の標的遺伝子として AGAP006058（ヘパラン硫酸 2-O-スルfオトランスフェラーゼ）などの候補が示唆されたが、発現レベルは概して低かった。
• 保存性: 逆転特異的な変異（SNP および indel）は、Anopheles gambiae 種複合体全体で保存されており、A. arabiensis や A. merus といった種においても 2La カリオタイプと共分離している。

意義および主張
本論文は、染色体逆転が、個別の転写エンハンサーの進化を通じて表現型の違いを駆動し得るという、最初の機能的な証拠を提供すると主張している。本研究は以下のことを示している：

1. 調節的分離: 2 つのカリオタイプは、それぞれのカリオタイプ内に固定された特定の塩基変異（SNP および indel）によって駆動される、機能的に異なる活性を持つエンハンサーを保持している。
2. 表現型変異のメカニズム: これらのエンハンサー活性の機能的な違いは、集団としてのゲノム特性（例：休息行動やマラリア感受性などの既知の大きな表現型の違い）に寄与する、カリオタイプ特異的な遺伝子調節プロファイルをもたらすと考えられる。
3. 手法上の進歩: リファレンス・アセンブリ比較と集団再シーケンシングを組み合わせることで、単一リファレンスゲノム研究に固有のアサーテインメント・バイアスを克服し、そうでなければ「不可視」であったままとなる調節変異を明らかにした。
4. 媒介制御への示唆: これらの具体的な調節メカニズムを理解することは、媒介生物の集団を単一のモノリスとして扱うのではなく、特定のセグメント（例：「スーパースプレッダー」表現型に関連する 2L+a グループ）を標的とした媒介制御戦略を策定する上で役立つ可能性がある。

著者らは、特定の機能的エンハンサーを同定したものの、複雑な表現型（マラリア感受性など）に対するこれらの特定の要素の正確な寄与や完全な調節ネットワークの解明には、さらなる調査が必要であることを認め、控えめなトーンを維持している。また、現在の技術では、その移動性や反復的な性質に起因するトランスポゾン（TE）関連エンハンサーを完全に解明することには限界があることも指摘している。