Lifestyles of Gypsy-family transposons shape their regulatory mechanisms

問題提起
トランスポゾン（TE）は真核生物ゲノムの大部分を占めており、ゲノムの安定性に脅威をもたらす。Drosophila melanogaster（キイロショウジョウバエ）におけるTEと宿主ゲノムの共進化、特にPIWI結合RNA（piRNA）によるサイレンシング経路についてはよく研究されているが、これらの要素の種を越えた進化パターンに関する理解は不足している。具体的には、ドロソフィラ属（drosophilid）の広範な系統におけるGypsyファミリーのLTRレトロトランスポゾン（これらはDrosophilaにおいて最も活性の高いグループである）の進化動態は依然として不明である。主要な生物学的問いは、これらの要素の「ライフスタイル」に関するものである。すなわち、どのようにして生殖細胞限定的な活性から、体細胞から生殖細胞への感染へと移行するのかという点であり、このプロセスはENV（Envelope）やsORF2といったウイルス由来のタンパク質によって媒介されると仮定されている。

手法
著者らは、313種のドロソフィラ属を代表する、公開されている450個のロングリードゲノムアセンブリのリソースを活用した。品質管理フィルタリングの後、約5,000万年の進化にわたる10科を網羅する249個の高品質ゲノムを精査した。

研究では多段階の計算パイプラインを用いた：

1. TEアノテーション: Extensive de novo TE Annotator (EDTA) および RepeatMasker を使用し、グローバルなTEランドスケープを評価した。
2. Gypsy特異的な精査: 特定のGypsyファミリーLTRレトロトランスポゾンを同定するために、精査されたD. melanogasterのTEライブラリに基づく保存されたポリメラーゼ（POL）オープンリーディングフレーム（ORF）を用いた隠れマルコフモデル（HMM）を構築した。これらのHMMを用いて、249個のゲノムにおける相同なPOL配列を同定した。
3. コンセンサス構築: 相補的ベストヒットを用いてオーソログを割り当てた。フルレングスのコンセンサス配列は、ヒットの拡張およびMAFFTやCIAlignなどのツールを用いた多重配列アライメント（MSA）によって生成された。
4. 検証: LTRの同一性（最近の挿入における5'および3' LTRは同一であるはずである）を確認すること、および小細胞RNAシーケンシング（sRNA-seq）データをマッピングしてpiRNA経路の標的化（アンチセンスマッピング、5' Uバイアス、および23-29 ntの長さ）を確認することによって、パイプラインの精度を検証した。
5. 進化解析: 非同義置換と同義置換の比率（dN/dS）および遷移・転換比（Ts/Tv）を算出し、選択圧を評価した。また、水平トランスポゾン転移（HTT）を示す系統的不一致を検出するために、POL配列の系統樹を構築した。
6. 制御解析: D. melanogasterのモチーフに対するFIMOを用いたスクリーニングにより、5' LTR領域の転写因子結合部位（TFBS）を探索した。TFBSの存在とENV ORFの状態（フルレングスか、あるいは欠失・断片化しているか）との相関を調べるために、混合効果ロジスティック回帰を用いた。

主な結果

• ゲノムランドスケープ: レトロエレメントは、DNAトランスポゾンよりも有意に大きな割合でドロソフィラ属のゲノムを占めており、GypsyファミリーはLTRレトロトランスポゾンの中で最も高いゲノムカバレッジを示した。D. melanogasterに類似したGypsy TEは広く分布しているが、その存在量は系統（clade）によって大きく異なる。
• ORFの進化動態:
  • sORF2: このタンパク質は、Gypsy/mdg1ファミリーに関連しており、高度に保存されており、ドロソフィラ属のほぼすべての（>80%）コンセンサス配列に存在するため、この系統における普遍的な構成要素であることが示唆される。
  • ENV: 対照的に、Gypsy/gypsyファミリーに関連するENVタンパク質は、系統全体を通じて繰り返される消失と断片化のパターンを示している。多くのGypsy/gypsy要素は機能的なENV ORFを失っており、これは体細胞から生殖細胞への感染能から、生殖細胞限定的な活性へのシフトを示唆している。
  • 選択: すべての主要なORF（POL, GAG, sORF2, ENV）は、浄化選択（dN/dS < 1）の証拠を示しているが、ENVは他のORFよりも高いdN/dS比を示しており、これはその頻繁な不活性化と一致している。
• 水平転移: TEのPOL配列の系統樹は、ホストの種系統樹としばしば矛盾しており、特にD. melanogaster系統とZaprionus属の間での広範な水平トランスポゾン転移（HTT）イベントの強力な証拠を提供している。
• 転写制御: 研究では、LTR領域における特定のTFBSモチーフがENVの有無と相関していることが特定された。
  • 体細胞発現: フルレングスのENVの存在は、卵巣の体細胞で発現する転写因子であるtraffic jam (tj) および CHES-1-like の結合部位と有意に関連している。
  • 生殖細胞発現: ENVの欠如は、生殖細胞で発現する Cf2, GATAd, CrebB, および BEAF-32 の結合部位と有意に関連している。
  • これらの関連性は、Gypsy TEが組織特異的な転写活性化を確実にするために、ホストの発達経路（特に卵巣の発達）をハイジャックしていることを示唆している。

意義と主張
本論文は、ドロソフィラ属におけるGypsyファミリーLTRレトロトランスポゾンのライフスタイルに関する「前例のない洞察」を提供すると主張している。249種から3,438個のコンセンサス配列を精査することで、著者らは以前には到達不可能であった規模でのTE進化の研究のための枠組みを確立した。

本研究は、これらの要素の進化史が、垂直継承と水平転移の複雑な相互作用によって形作られていることを強調している。決定的なことは、Gyencyのサブファミリーにおける進化戦略の二分法を明らかにしている点である。すなわち、Gypsy/mdg1ファミリーは、体細胞から生殖細胞への転移のための安定し保存されたメカニズム（sORF2）を維持している一方で、Gypsy/gypsyファミリーは、ENVタンパク質の獲得とその後の消失という動的な歴史を示している。この消失は転写制御の変化を伴っているようであり、要素は体細胞特異的な調節モチーフ（tj, CHES-1-like）を失い、生殖細胞特異的なモチーフを保持または獲得している。

著者らは、本データセットと知見が、TEがホストのサイレンシング防御を回避し、異なる細胞ニッチ間での増殖を確実にするために、どのようにその転写制御メカニズムを適応させるかという、TEとホストゲノムの共進化に関する将来の研究を促進すると結論づけている。本研究は、HTTと不完全な系統分離（incomplete lineage sorting）などの他の現象を区別することが困難であることを認めつつ、決定的な解決策ではなく、基礎的なリソースとして提示されている。