IDENTIFICATION OF CANDIDATE GENES UNDERLYING A NOVEL QTL FOR ROOT LENGTH AND DIAMETER IN BARLEY IN RESPONSE TO SALINITY

技術要約：塩ストレスに対するオオムギの根の長さおよび直径に関する新規QTLの基礎となる候補遺伝子の同定

問題提起
土壌の塩分は、世界の作物生産性を制限する重要な非生物的ストレス要因であり、気候変動によってこの問題はさらに悪化すると予想されている。塩耐性は浸透圧およびイオンの複雑なメカニズムを伴うが、根系構造（特に、塩条件下で根の伸長と養分吸収を維持する能力）の遺伝的制御については、依然として理解が進んでいない。塩分耐性に関する多数の量的形質遺伝子座（QTL）がオオムギ（Hordeum vulgare L.）において特定されているものの、根の成長反応の基礎となる因果遺伝子や分子メカニズムは、大部分が不明なままである。本研究は、塩耐性の向上を目指した育種を促進するために、塩ストレス下での根の形態形成（初生根の長さや直径など）を制御する特定の遺伝子座および候補遺伝子を同定する必要性に応えるものである。

方法論
本研究では、対照的な塩応答を示す2つのオオムギ系統（根の伸長を維持する耐塩性のClipperと、感受性のSahara 3771）を交配して得られた倍数体半数体（DH）マッピング集団を利用した。

1. 表現型解析： 144のDH系統と親系統を、制御環境下においてコントロール条件および塩ストレス（100 mM NaCl）条件下で栽培した。測定した形質は、初生根の長さ（SRL）、平均根径（RD）、バイオマス、およびイオン含有量（Na+, K+）である。空間的な変動を考慮するために線形混合モデルを用いてデータを分析し、QTL解析のための最良線形不偏推定値（BLUEs）を生成した。
2. QTLマッピング： 初期のQTL解析は、レガシー・リンケージマップを用いて行われた。有意なQTLが染色体2H上に検出された。この領域を精査するため、RNAシーケンシングによって発見された16個の一塩基多型（SNP）に対して、競合アレル特異的PCR（KASP）アッセイを設計した。これらのSNPを用いて染色体2Hの濃縮リンケージマップを作成し、より高解像度なQTLマッピングを可能にした。
3. バルク分離RNAシーケンシング（BSR-Seq）： 根の表現型に基づき、20のDH系統を選択した（耐塩性系統［長い根］10系統、感受性系統［短い根］10系統）。これらを、コントロールおよび塩処理の両条件下で4つのバルク（耐性用2つ、感受性用2つ）にプールした。根端をシーケンシングし、表現型の違いに関連する差発現遺伝子（DEG）および配列多型を特定した。
4. 候補遺伝子の同定： 本研究では、QTLマッピングの結果とBSR-Seqデータを統合した。染色体2Hの短腕上の2.5-Mbの候補領域（500,000 ～ 3,000,000 bp）を定義した。以下の基準に基づいて、遺伝子の優先順位を決定した：
   • 耐性バルクと感受性バルク間での差発現。
   • 塩分に対する転写応答。
   • 親系統間における高影響または中程度のインパクトを持つ配列変異（SNP/indel）の存在。
   • 機能注釈。

主な結果

• 表現型の変異： DH集団は、根の形質において連続的な変異を示した。耐性系統（Clipper由来）は、感受性系統（Sahara 3771由来）よりも塩ストレス下で根の長さをより良く維持した。
• QTLの局在化： SNP濃縮マップを用いたQTL解析により、初生根の長さと平均根径の両方に有意に影響を与える、染色体2Hの短腕遠位領域（約2.1 cM）における主要な遺伝子座が特定された。このQTLは、塩条件下でより明確になった。Sahara 3771のアレルは、より短く太い根と関連しており、Clipperのアレルは、より長く細い根と関連していた。
• BSR-Seqの知見：
  • 塩応答遺伝子： 耐性バルクと感受性バルクの両方において、ストレス応答経路（グルタチオン代謝、MAPKシグナル伝達、フェニルプロパノイド生合成）の濃縮が見られた。
  • 表現型依存的遺伝子： 耐性バルクと感受性バルク間で差発現した遺伝子は、DNA結合転写因子およびRNAバイオ合成の調節に富んでいた。
  • 候補領域： 2.5-Mbの区間内に、76の高信頼性遺伝子モデルが特定された。そのうち12個がバルク間で差発現しており、31個の遺伝子が高/中程度のインパクトを持つ変異を含んでいた。
• 優先された候補遺伝子： 6つの遺伝子が強力な候補として優先された。そのうち2つは、オーキシン恒常性の役割において特に注目されるものである：
  1. HORVU.MOREX.r3.2HG0096800： オーキシン生合成の鍵となる酵素であるYUCCA5に相同な、フラビン含有モノオキシゲナーゼをコードしている。これはミスセンス変異を含み、Saharaアレルと比較してClipperアレルにおいてダウンレギュレートされている（log2FC = -1.27）。
  2. HORVU.MOREX.r3.2HG0096550： 極性オーキシン輸送に関与するABCBトランスポーターをコードしている。これもまた、Saharaアレルと比較してClipperアレルにおいてダウンレギュレートされており（log2FC = -4.33）、ミスセンス多型を含んでいる。
     その他の候補には、セルロース合成、グリコシル転移酵素活性、チトクロームP450機能、およびペタントリックペプチドリピート（PPR）タンパク質に関与する遺伝子が含まれる。

意義と主張
本論文は、コントロール条件および塩条件下（塩ストレス下）の両方において、オオムギの根の形態形成の主要な決定因子である新規のゲノム領域を染色体2H上に特定したと主張している。QTLマッピングとBSR-Seqを統合することで、研究チームは探索範囲を2.5-Mbの区間に絞り込み、特にオーキシン生合成（YUCCA様）および輸送（ABCB）に関与する特定の候補遺伝子を浮き彫りにした。

著者らは、ClipperとSahara 3771の間で見られる対照的な根の表現型（長く細い vs 短く太い）は、これらの候補遺伝子によって制御される、特定の塩誘導的な転写応答ではなく、オーキシン恒常性の構成的な差異によって駆動されている可能性があると示唆している。具体的には、本研究は、YUCCA様遺伝子およびABCBトランスポーター遺伝子の両方が、Sahara 3771アレル（短く太い根に関連）を持つ系統においてより高い発現を示し、Clipper系統（長い細い根に関連）ではダウンレギュレートされていることを指摘している。本研究は、これらの遺伝子の因果的役割を検証するための将来的な機能検証の基礎を提供し、塩分を含む土壌における改良された根系を持つ作物の育種に向けた標的を提示している。本研究は、塩耐性は複雑であるが、根の形態形成とその基礎となる遺伝的メカニズムに焦点を当てることは、作物改良のための実行可能な道筋であることを強調している。