ゲノムワイド往復半角化分析によるサッカロマイセス酵母の種間の熱耐性の遺伝学的マッピング

サマリー

シーケンシング(RH-seq)を介した相互異形性は、種間の形質の違いの遺伝的基礎をマッピングする強力な新しい方法である。ヘミジゴットのプールはトランスポゾン変異によって生成され、その適合性は、高いシーケンシングを使用して競争力のある成長を通じて追跡されます。得られたデータの分析は、形質の基礎となる遺伝子を特定します。

要約

現代遺伝学の中心的な目標は、野生の生物が表現型で異なる方法と理由を理解することです。現在までに、この分野は、種の個体間の交配から組換え子孫間のDNA配列変異体と表現型との関係をトレースするリンケージおよびアソシエーションマッピング法の強度に大きく進歩した。これらのアプローチは、強力ですが、生殖的に分離された種間の形質の違いに適していません。ここでは、互換性のない種に容易に適用できる自然形質変動のゲノム全体の解剖のための新しい方法について説明する。私たちの戦略、RH-seqは、相互ヘミジゴット試験のゲノム全体の実装です。我々はそれを利用して、その姉妹種S.パラドックスに対する酵母サッカロマイセスセレビシエの顕著な高温増殖の原因となる遺伝子を同定した。RH-seqはトランスポゾン変異を利用して相互ヘミジゴットのプールを作成し、ハイスループットシーケンシングを介して高温競争を通じて追跡されます。ここで説明した RH-seq ワークフローは、新芽酵母クレードの古代の複雑な形質を解剖する厳格で公平な方法を提供し、遺伝的マッピングのためのゲノムカバレッジを確保するためにリソース集約的な深いシーケンシングが必要であることを警告します。シーケンシングコストが低下するにつれて、このアプローチは真核生物全体での将来の使用に大きな期待を持っています。

概要

フィールドの夜明け以来、それは野生の個体間の変動の機械的基礎を理解するために遺伝学の主要な目標となっています。関心のある形質の根底にある遺伝子座をマッピングすると、新しい遺伝子は診断や薬剤の標的として即座に使用され、進化の原理を明らかにすることができます。この目的に向けた業界標準は、リンケージまたはアソシエーション1を介して集団全体の遺伝子型と表現型の関係をテストすることです。これらのアプローチは強力であり、彼らは1つの主要な制限を持っています - 彼らは相互不妊の個人間の交差から組換え子孫の大きなパネルに依存しています。そもそも子孫を形成することができない種の研究では役に立たない。したがって、この分野は、生殖的に単離された種2の間の形質の違いを公平に解剖する能力がほとんど持っていなかった。

本研究では、種間の形質変動の遺伝的基盤のゲノムスケール調査に関する新しい方法RH-seq³の技術的基盤を報告する。このアプローチは、相互ヘミジゴット試験4、5の非常に平行なバージョン^であり、これは最初に2つの遺伝的に異なる背景の間の対立性の差異の表現効果を評価する方法として考案された特定の軌跡(図1A)。このスキームでは、2つの異なる個体が最初にハイブリッドを形成するために交配され、そのゲノムの半分はそれぞれの親から来ている。この背景には、複数の株が生成され、それぞれが各親の軌跡の対立遺伝子の中断または削除されたコピーを含む。これらの株は、それらは、関心のある軌跡を除いてゲノムのいたるところに浸留し、ハプロイドとみなされ、それぞれが1つの親の対立遺伝子しか欠け、残りの対立遺伝子が由来するので、多角質である。他の親。これらの相互ヘミジゴット株のフェノタイプを比較することにより、遺伝子座の変異体が相互の間の唯一の遺伝的差異であるため、操作された遺伝子座におけるDNA配列変異体が関心のある形質に寄与するかどうかを結論付けることができる。ヘミジゴット株。このようにして、種間の遺伝的差異を、よく制御された実験セットアップでそれらの間の見型的な違いに結び付けが可能である。現在までに、この試験の応用は候補遺伝子フレームワーク、すなわち候補遺伝子座の自然変動が形質に影響を与える可能性があるという仮説が既に手に入っている場合である。

以下では、酵母をモデルシステムとして用いて、ゲノムスケールの相互半角膜スクリーンのプロトコルをレイアウトする。我々の方法は、種間で生存可能な無菌F1ハイブリッドを生成し、トランスポゾン変異形成を受けることによって、ヘミジゴット変異体のゲノム補体を作成する。我々は、ヘミジゴットをプールし、シーケンシングベースのアッセイでそれらの表現型を測定し、所定の遺伝子の2つの親の対立遺伝子を有するプールのクローン間の頻度の違いをテストする。その結果、種間の変異体が関心のある形質に影響を与える遺伝子座のカタログが得られます。我々は、2つの新芽酵母種、サッカロマイセスセレビシエとS.パラドックスの間の熱耐性の違いの遺伝的基礎を解明するためにRH-seqワークフローを実装し、6年前に約500万年前に発散した。

プロトコル

1. 形質転換のためのピギーバク含有プラスミドの調製

1. 単一のコロニーにストリークは、LB +カルベニシリン寒天プレート上にプラスミドpJR487を収容する大腸菌株を収容する。37°Cで1泊、または単一のコロニーが現れるまでインキュベートします。
   注:プラスミドpJR487がクローン化された方法の説明は、前の作品3で見つけることができます。
2. 2LガラスフラスコにpJR487を含有する大腸菌の単一コロニーを有するLB+カルベニシリンの1Lを100μg/mLで接種する。飽和するまで200 rpmで振りながら37°Cで一晩成長します(OD₆₀₀ ≥ 1.0)。
3. 製造元の公開プロトコルで指示されているように、大規模なプラスミド準備キットを使用して培養物からプラスミドDNAを精製します(詳細は材料の表を参照)。溶出バッファーの5 mLを37°Cに温めた5mLで10分間インキュベーションした後にDNAを溶出させる。
4. 分光光度計を用いてプラスミドDNAの量と質を測定します(詳細は材料表を参照)。
5. ステップ1.2~1.4を繰り返し、A260で少なくとも11mgのプラスミドDNAの合計が単離されるまで:A₂₈₀比が少なくとも1.8個単離されるまでである。これは、効率に応じて、いくつかの準備を取る場合があります。
6. すべてのプラスミドの準備物を1つのチューブに混ぜ、溶出バッファーまたは水で20 mLまでの総容積を持って来る。分光光度計で最終的な量と品質を再度測定します。プラスミドの濃度は、この最終20 mL体積において少なくとも538 ng/μLでなければなりません。濃度が538 ng/μLより高い場合は、溶出バッファーまたは水でプラスミドを538 ng/μLに希釈します。プラスミドは、使用するまで数週間まで4°Cで保存することができる。

2. 標的化されていないゲノム全体の相互ヘミジゴットのプールを作成する

1. 形質転換のためのハイブリッド酵母細胞の調製
   1. -80 °Cの冷凍庫ストック株からYPD寒天プレート上の単一のコロニーにJR507をストリークアウト。26 °Cで2日間、またはコロニーが現れるまでインキュベートします。
      注:JR507は、S.セレビシアDBVPG1373およびS.パラドキスZ1(四面体解剖顕微鏡を用いて)3のハプロイド胞子の単細胞交配を介して作られたハイブリッド株である。
   2. JR507の単一のコロニーで250 mLのガラスフラスコで液体YPDの100 mLを接種し、28°C、200 rpmで24時間、または静止相に達するまで振ります。
   3. 翌日、一晩培養の600nm(OD600)で光学密度を測定する。新鮮な液体YPDで一晩培養の一部を0.2のOD 600と500 mLのOD₆₀₀に新しい1Lガラスフラスコに戻すことによって、新しい文化を作成します。
      注: 一晩培養に OD₆₀₀が 5.0 で、C が光学密度で V がボリュームである場合の逆希釈の計算例:
      
      したがって、飽和一晩培養の20mLを液体YPDの480 mLに添加し、0.2のOD₆₀₀で合計500mLの培養を行う。
   4. ステップ 2.1.3 をさらに 3 回繰り返し、4 つの 1 L ガラスフラスコで 0.2 の OD₆₀₀で合計 4 つの 500 mL 培養を行い、4 つの新しい培養物すべてに同じ一晩培養を使用します。28°Cで6時間(2~3世代)振り、200rpmでインキュベートします。
   5. 500 mL 培養物のうち 2 つを組み合わせて、1 L 培養を作成します。残りの 2 つの 500 mL 培養を組み合わせて、別の 1 L 培養を作成します。この時点で、2 つの 1 L カルチャがあります。これらの 1 つの L 培養のそれぞれは、次の手順で pJR487 を使用して変換されます。
2. pJR487をハイブリッド酵母細胞に変換
   1. 1L培養物を20本のプラスチック円錐管に20個の50mLのアリコットに分割し、合計40本のチューブを使用します。20本のチューブを脇に置き、一度に20本のチューブで次の手順を実行します。
   2. 酵母細胞をペレットに1,000 x gで3分間20本のチューブのそれぞれを遠心分離する。上清を捨てます。
   3. 渦によって無菌H2 Oの25 mLで各ペレットを再懸濁する。1,000 x gで3分間遠心分離機。上清を捨てます。
   4. 渦によって1x TE、0.1 M LiOAcバッファーの5 mLで各ペレットを再中断します。1,000 x gで3分間遠心分離機。上清を捨てます。
   5. 手順 2.2.4 を繰り返します。細胞が遠心分離している間、39.52%ポリエチレングリコール、0.12 M LiOAcおよび1.2xトリス-EDTAバッファー(12 mMトリス-HClおよび1.2 mM EDTA)の溶液の少なくとも120 mLを調製する。氷の上に保管してください。
   6. 形質転換のためにプラスミドDNAを調用するために、まず100°Cでサケ精子DNAの4 mLを5分間沸騰させ、すぐに氷の上で5分間冷却する。次いで、pJR487の20mL(セクション1で得られた)を538ng/μLの濃度で混合し、24mLの総体積に対して冷却されたサケ精子DNAの4mLと混合する。使用するまで氷の上に保管してください。
   7. 各細胞ペレットの上にサケ精子DNAと混合プラスミドDNAの600 μLを追加します。まだ再中断しないでください。
   8. ステップ2.2.5で作られたPEG-LiOAc-TE溶液の3 mLを各ペレットに追加します。上下にピペッティングし、渦を巻くことによってペレットを再中断します。
   9. 室温で10分間各チューブをインキュベートします。
   10. 39°Cに設定された水浴で26分間、各チューブにヒートショックを与えます。
       注:数分ごとに、各チューブを反転して、細胞がチューブの底部に沈着するのを防ぎます。
   11. 1,000 x gで3分間各チューブを遠心分離します。上清を廃棄し、渦によってYPDの10 mLで各ペレットを再懸濁する。新しいガラスフラスコにすべての20チューブを組み合わせます。細胞の総体積は〜200 mLでなければなりません。
   12. 66.6 mLのセルを新しい1Lガラスフラスコに移し、液体YPDで500mLの体積まで持ち上げます。形質転換細胞の200 mL全体を使用するには、さらに2回繰り返します。新しい 500 mL 培養ごとに OD₆₀₀を測定します (OD₆₀₀ ~0.35-4 を期待します)。
   13. 3つのフラスコを28°Cで2時間振り、200rpmで回復(<1世代)。
   14. 3つのフラスコのそれぞれに300mg/mL G418の0.5 mLを加え、最終的な濃度300 μg/mL G418に加え、28°C、200rpmで振り返ります。
       注: この手順の前に、変換されたハイブリッド セルは変換から回復しています。G418の添加に伴い、プラスミドpJR487の存在が選択される。形質転換中にプラスミドを取り上げなかった細胞は死に始める。
   15. 残りの 20 本の円錐形のチューブで、手順 2.2.2 ~ 2.2.14 を繰り返します。この時点で 6 つの 1 L ガラスフラスコがあり、それぞれに G418 が追加された 500 mL のセルがあります。
   16. 28°Cで細胞の6つのフラスコをすべてインキュベートし、200rpmで振り、約2日間、または各フラスコで約2.3_のOD600に達するまで。6 つのフラスコをすべて組み合わせて、単一のカルチャを作成します。
       注: このカルチャ内のすべてのセルは下流のステップでは使用されませんが、このような大容量を使用する目的は、できるだけ多くの固有の変換イベントを作成し、それらをすべてプールすることによって、単一の変換全体のバイアスを正規化することです。一緒に。
   17. 2.2.16 で作成されたカルチャを使用して、YPD + G418 (300 μg/mL) の 500 mL で 2 つの新しい 1 L フラスコを 0.2 の OD₆₀₀に接種します。廃棄できる残りカルチャがあります。
   18. 28°Cで1Lフラスコを一晩でインキュベートし、200rpmで振り、それぞれが〜2.2(〜3.5世代)のOD₆₀₀に達するまで。両方のカルチャを単一のカルチャに結合し、結合されたカルチャの OD₆₀₀を再度測定します。
       注:この時点で、培養物は、プラスミドpJR487を収容する細胞からほぼ完全に構成されるべきである。細胞の集団の一部では、PiggyBacトランスポゾンはプラスミドからプラスミドを発現するトランスポザーゼによってプラスミドからゲノムに転置されるであろう。しかし、トランスポザーゼの発現を継続すると、選択の過程で転置が起こし、遺伝子型と表現型の関係が不明瞭になります。次のいくつかのステップの目的は、プラスミドの存在に対してカウンターセレクションを実行し、トランスポザーゼの発現がこれ以上ないことを確実にすることである。得られたプールは、トランスポゾンをゲノムに統合する有無にかかわらず細胞の混合物であるが、その後のマッピングステップ中にトランスポゾンを含む細胞のみが検出される。トランスポザーゼが発現する変換の時間は、プラスミド符号化が失われる前に、変異後の所定のクローンが複数のトランスポジソン挿入を持つ可能性を支配し得る。一度に1つの遺伝子の分析で「二次的」変異として現れるこれらの頻度は、変異発生後に定義されたコロニー数を配列し、そのDNAを組み合わせ、独立した挿入数を配列することによって推定することができる。プール内の位置。
   19. 1,000 x gで3分間、この培養の遠心分離機25 mL。25 mL にあるセルの合計 OD₆₀₀単位の数を計算します (下の計算例を参照)。上清を廃棄し、渦によって1.85 OD₆₀₀/mLの細胞懸濁液を作成するのに十分なH2 Oで再懸濁する。
       注: 結合培養の OD₆₀₀が 2.2 であった場合の水中の細胞の再懸濁の計算例:
       
       そこで、細胞培養の25mLを紡糸し、上清を廃棄した後、細胞と水の総体積を最大29.7mLまで持ち込むのに十分なH2Oを細胞に加える(細胞ペレットも体積を持つため、H2Oの29.7mL未満を添加する)。
   20. ガラスビーズを使用して、5-FOAを持つ12の大きな正方形の完全な人工寒天プレートのそれぞれに水中に再懸濁した細胞のプレート1 mL。各プレートを28°Cで1~2日間、またはプレート上に芝生が形成されるまでインキュベートします。
   21. 小さな無菌スクイージを使用して、6枚のプレートのそれぞれから細胞をこすり落とし、35 mLの無菌水を含むチューブに入れます。他の6枚のプレートで、細胞と水の合計2本のチューブを繰り返します。すべてのセルサスペンションを1つのチューブに組み合わせます。水をブランクとして使用して、この懸濁液のOD₆₀₀を測定します。OD₆₀₀/mLの細胞濃度を水で44.4 OD₆₀₀単位/mLにします。私たちの経験では、転置効率(URA-であるKAN+細胞の割合)は平均50%です。
   22. 保存するセルの-80 °C冷凍庫ストックの数を決定します。各アリコートは、将来的に単一の実験に使用することができます。
       注: プールの生成にかかる時間を考えると、偶発的な誤用や反復実験の実行時に複数のバイアルを保管してください。20~30株は妥当な数です。
   23. 各冷凍庫のストックは10%DMSOの1 mLの細胞の40 OD₆₀₀単位を含んでいる。DMSOの100 μLに900 μLの細胞を加えます。作成されたフリーザーストックの合計数に対して繰り返します。将来の使用のために-80 °Cでそれぞれを保管してください。

3. プール形式での相互ヘミジゴットの選択

1. -80°C冷凍庫から、室温でセクション2からプールされた相互ヘミジゴットの単一のアリコートを解凍する。
   注:アリコートが解凍したら室温で長時間座らせないようにし、すぐに使用してください。
2. 1 mL アリコート全体を使用して、250 mL のガラスフラスコに 150 mL の液体 YPD を接種します。この培養物_のOD600を測定し、次いで28°Cでインキュベートし、200rpmで振り、〜7時間、または培養が2〜3個の集団倍増を経るまで。この時点で、培養は選択中の培養物を接種するために使用される準備ができています。
   注: 計算例: 元のフラスコの OD₆₀₀が 0.25 を測定する場合は、少なくとも 1.0 の OD₆₀₀に達するまでカルチャをインキュベートします。「タイムゼロ」(T-0)でサンプルポイントが望まれる場合、選択前にヘミジゴテ集団を調べるための方法として、ペレット当たりの培養量を1,000xgで3分間遠心分離し、上清を廃棄し、 -80 °C.
3. 成長したヘミジゴットプールを使用して、高温(39°C)と許容温度(28°C)の両方で、適切なレプリケートスキームで選択するための培養物を接種します。少なくとも、各温度で3つの生物学的複製選択培養物を、合計6つの選択培養物に設定する。
   1. 液体YPDが付いている2 Lガラスフラスコの合計500 mLの各選択培養を作成し、0.02のOD₆₀₀に接種する。各選択培養物を28°Cまたは39°Cのいずれかで100rpmで振り、6~7個の集団倍増が起こった(~1.28-2.56のOD₆₀₀に相当)。すべての選択カルチャの最終的な OD₆₀₀を可能な限り密接に一致するようにしてください。
      注:28 °Cの選択培養物は、39°Cの選択培養物よりも速く成長します。その結果、39°Cの選択培養物は、インキュベーターでより長い期間を過ごします。インキュベーターで費やされた合計時間数に関係なく、準備が整った各フラスコで次の手順を進めます。私たちの経験では、28°Cまたは39°Cの培養物は、それぞれ〜12または〜18時間かかり、約2.0のODに達しました。長い選択は、小さなフィットネス効果を増幅する利点を持つことができますが、また、任意の1つの遺伝子/対立遺伝子のトランスポゾン変異体間のフィットネスの最終的な分布にノイズを導入するde novoの背景変異が発生することを可能にします。そのため、RH-seq実験では選択時間を制限することが重要です。
4. 各選択培養から細胞ペレットを収穫する。ライブラリの準備とシーケンシングのための技術的な複製として、各選択培養物からこのボリュームの少なくとも4つのペレットから、1,000 x gで7 OD₆₀₀単位の細胞と遠心分離機を得るために必要な体積を計算します(セクション4および5を参照)。、以下)。上清を捨て、-80°Cで保管してください。
   注: 選択フラスコの最終 OD₆₀₀が 2.0 の場合の例:
   

4. Tn-seqライブラリ構築とイルミナシーケンシングにより、トランスポゾン変異型ヘミジゴットの豊富さを判定

1. シークエンスされるセクション3から各細胞ペレットを氷上で融解する。
2. メーカーの指示に従って酵母gDNA精製キットを使用して、各細胞ペレットから総ゲノムDNA(gDNA)を分離します。65°Cに温めた溶出バッファーの50μLでDNAを再中断します。
3. フッ素計を用いて各ペレットからのgDNA量を定量します。次の手順を使用して Tn-seq の次世代シーケンシング (NGS) ライブラリーを作成するために各セル ペレットに必要な gDNA の最小総量は 1 μg です。
   注: ライブラリの作成には 1 μg 未満の gDNA を使用できますが、ライブラリの最終的な量と品質は損なわれます。
4. Tn-seq ライブラリ⁷を作成するための確立されたプロトコルに従ってください。このプロトコルに固有の次の関連情報に注意してください。
   1. gDNAせん断、修理およびアダプタライゲーションを終了した後、PCRを介してトランスポゾンを含むgDNAを増幅します。この PCR では、PiggyBac トランスポゾンアダプターと NGS アダプターに固有の次のフォワード プライマーとリバース プライマーをそれぞれ使用します。
      フォワード (N – ランダムヌクレオチド)
      5' ATGATACCGACCACCGAGATCTACCTCTTCTACACGACG
      CTTCCGATCTNNNNNNNナグナガットカアットカアガガットGCGTCAATCTCAATCTCAATCTCAATAGCGTCAAT3'
      リバース(Nsのストレッチは、多重化に使用される一意の6bpインデックスを表します。インデックスの詳細については、以下を参照してください)
      5' CAAGCAGAAGAGGGGAGAGAGAGGAGGAGATNNNNNNGTGACTGGAGTCAG
      ACGTGTGTGCTCTTCCGCT 3'
   2. サイズ選択ビーズを含む付属のクリーンアップ手順を使用して、最終的なライブラリ内のクローンフラグメントの比率を最小限に抑えます。
      注: 選択培養に関するこれまでの最小反復要件に従った結果、シーケンス用に 24 個の gDNA サンプルが存在します。シーケンスの現在のコストを考えると、各サンプルが単独で実行されることはほとんどありません。同じレーン上のサンプルを組み合わせるには、複数のリバースプライマーを作成し、それぞれに固有の 6 ベース ペア インデックスを作成します。異なるインデックスを持つサンプルは、同じシーケンス レーンに結合し、その後計算で分離できます。
5. 8 レーンにわたる NGS テクノロジを使用して、各ライブラリからシーケンス シングルエンド 150 bp 読み取りします。
   注:必要なシーケンシング読み取りの量は、前のステップで調製したライブラリの品質に大きく依存します(すなわち、相互ヘミジゴット由来のDNAを表すトランスポゾンDNAを実際に含むライブラリ内のDNAの割合)。これに寄与する主な要因は2つあります。第 1 に、統合トランスポゾンのないセルはプール作成中に対して逆選択されないため、各培養はトランスポゾンの有無にかかわらずセルの組み合わせになります。第二に、トランスポゾン含有逆半血のゲノム内であっても、ゲノムのほとんどは配列を含むトランスポゾンではなく、このgDNAはやむを得ずライブラリー製剤の一部となる。トランスポゾン含有DNAの最終的なPCR増幅の目的は、これら2つのバックグラウンドgDNA源に対するトランスポゾン含有DNAの比率を増加させる。この増幅がより効率的であるほど、読み取りの割合が高いほど、ダウンストリーム解析に使用できるようになります。ライブラリの品質が低いほど、読み取りの割合が増えるほどトランスポゾン DNA は含まれないため、より多くのシーケンシングを行う必要があります。上記の制約を考えると、8レーンのシーケンシングは、相互ヘミジゴットの豊富さを合理的な程度に追跡することができる。シーケンシングを増やすと、より深い分析が可能になります。

5. トランスポゾン挿入とRH-seq解析の位置のマッピング

注: 次のデータ分析は、カスタム Python スクリプト (https://github.com/weiss19/rh-seq でオンラインで見つかった) で行われましたが、他のスクリプト言語を使用してやり取りできます。以下に、プロセスの主要な手順の概要を示します。個々のレプリケート読み取りファイルを組み合わせることに注意しない限り、次の手順を実行します。

1. アダプターのシーケンスを読み取りから除外し、インデックスに従って各反復の読み取り値を分離します。
2. トランスポゾンゲノム接点を含む読み取りを検索します。これを実現するには、トランスポゾン CAGACTATCTTTTTagGGTTAA の最後の 20 個の基本ペアを各読み取り内で検索します。このシーケンスを含まないすべての読み取りは破棄します。
   注:私たちの経験では、トランスポゾンの終わりにマッピングする読み取りの割合は83-95%です。
3. トランスポゾンの3'終端の下流のシーケンスのみを含むトランスポゾン含有読み取りをトリミングします。この配列を酵母ゲノムにマッピングすることにより、読み取りごとにトランスポジット挿入のゲノムコンテキストを決定します(下図5.4)。
4. BLAT または同等のマッピング ツールを使用して、トランスポゾンの下流のシーケンスを S. cerevisiae DBVPG1373 x S.パラドカスZ1 ハイブリッド ゲノムにマッピングします (スクリプト名: map_and_pool_BLAT.py)。
   1. トランスポゾンの3'終端の下流で使用可能なシーケンスの50未満の基本ペアがある読み取りを破棄します。短いシーケンスは、一意にマップするのが困難です。
   2. BLAT を使用する場合は、ID = 95、タイル サイズ = 12 というパラメーターを使用します。
   3. S. cerevisiae S288c および S. パラドックス CBS432 の参照ゲノムの最新バージョンを連結してマッピングに使用する基本的なハイブリッドゲノムを作成します。
      注:ハイブリッドゲノム全体の個々の遺伝子のゲノム境界を説明する基本的な注釈ファイルは、上記のGithubリポジトリで見つけることができます(ファイル名:YS2+CBS432+plasmid_clean)。ハイブリッドゲノム内の単一の場所にマップする読み取りのみを使用します(つまり、S. cerevisiae または S. パラドックスのいずれかに固有です)。ゲノム全体の挿入イベントの均一な頻度が期待されます。ゲノム間の挿入位置の分布は、他の場所で報告^{されています 3}.
5. 一意のトランスポゾン挿入位置にマッピングされる読み取りの合計数を集計し、1 つのトランスポジット挿入変異クローンのセルから発生したと推測します。1 つのライブラリーからのすべてのこのような値の合計は、そのライブラリーのマップされた読み取りの合計数と呼ばれます。
6. 互いに 3 つの基本ペア内に複数の挿入マッピングがある場合は、これらすべてを 1 つの挿入ポイントに結合し、すべての読み取り値を最も高い読み取り数を持つ単一の場所に割り当てります。この値n_insertは、gDNAが配列された細胞ペレット内の挿入クローンの豊富さを表す。この時点で、n_{インサート}のリストがあり、それぞれが一意のマッピングトランスポジット挿入の豊富さ、シーケンスされたセルペレットごとに1つのリストがあります。
   注:PiggyBacトランスポゾンは、ゲノム内のTTAA配列、4塩基配列に挿入します。したがって、互いの 3 つの基本ペア内の挿入マッピングは、同じ TTAA サイトから発生している必要があると推測します。
7. シーケンスされた各ライブラリからの合計読み取りの数は若干異なるので、比較する場合は、すべてのファイル間で n_挿入値を正規化します。この操作は、個々のライブラリ (n_{ペレット)}からマップされた読み取りの合計数を表し、すべてのライブラリ 、ペレット>のすべての n_ペレットの平均を取ります。個々のライブラリのデータ内の各n_{インサート}にペレット>/n_ペレットの比率を掛け、_{インサート}を計算し、所定のトランスポゾン挿入クローンの正規化された豊富さを計算します。
   
   または、DESeq2⁸ (スクリプト名: total_reads_and_normalize.py) などの利用可能なツールを使用して、ライブラリサイズを推定することもできます。
8. すべてのライブラリにマップされたすべての挿入のセットを集計します。一部のライブラリで見つかったが、他のライブラリでは見つからない挿入の場合は、下流の計算に_insert = 1 を設定します。
9. 読み取りをフィルタリングして、アニテーション ファイルに従って遺伝子内に含まれる挿入を見つけます (スクリプト名: remove_NC_and_plasmid_inserts.py)。
10. 一意の挿入ごとに、各選択の技術的な反復 (各カルチャは 28 °C または 39 °C)、インサート>_{テクニカル}(スクリプト名: combine_tech_reps_V2.py) を計算します。
11. 一意の挿入ごとに、各温度ですべての insert_>技術的な平均を取_{ることによって、各}温度の生物学的反復全体の平均量を計算します。同時に、挿入ごとに変動係数を計算し、CV_{挿入、 インサート}>_{テクニカル}(スクリプト名: combine_bio_reps.py)全体を計算します。
    注:この時点で、各温度、28°Cおよび39°Cについて、固有のトランスポジット挿入、それらの平均豊富さと各生物学的反復間の変動係数のリストがあります。私たちの実験のためのこれらのデータは、他の場所^{で報告されています 3}.
12. 28 °C または 39 °C、挿入>_合計> 1.1、および CV_挿入、_合計≤ 1.5 (スクリプト名: filter_insert.py) を持つもののすべての挿入のリストをフィルター処理します。
13. 一意の挿入ごとに、ログ₂(<挿入>_{合計,28 °C} / インサーt>_{合計,39 °C)}を計算します。この値は、特定のトランスポゾン挿入変異クローンの「熱耐性」を表します(スクリプト名: fitness_ratios.py)。
14. 遺伝子と対立遺伝子(S.セレビシエまたはS.パラドックス)によってユニークな挿入物をすべてソートし、各対立遺伝子の挿入数を集計する。各対立遺伝子に少なくとも5個の挿入を持つ遺伝子のみが分析できるように遺伝子をフィルタリングします(スクリプト名:organize_and_filter_genes.py)。
    注:各対立遺伝子に対する複数のユニークな挿入は、その相互ヘミジゴットの熱耐性のより正確な尺度を可能にします。対立遺伝子あたりに必要な挿入数を減らすことは可能ですが、この測定の精度を損ない、より多くの遺伝子をテストできるようにすることで、複数の検査負担を増大させます。さらに、対立遺伝子あたりの挿入数が少なすぎる遺伝子をフィルタリングすると、非常に異なる表現型を付与する二次部位変異を抱える個々のヘミジゴットクローンの試験結果への影響を軽減するのに役立ちます。
15. 上記のフィルタリング後にデータセットに残っている各遺伝子について、S.セレビシアの全挿入物の耐熱性(対数2比)を、マン・ホイットニーU検定を用いたS.パラドックス対立遺伝子のサーチ耐性(log2比)と比較する。あるいは、回帰モデルを実装し、DESeq2⁸ (スクリプト名: mann_whitney_u.py) から調整することもできます。
16. Benjamini-Hochberg法を用いた複数のテストのp値を補正する。
17. 有意なp値を持つ遺伝子(例えば、≤ 0.01)は、2つの種間の熱耐性の違いに重要な遺伝子の候補である。

結果

我々はS.セレビシエとS.パラドックスを交配して無菌ハイブリッドを形成し、トランスポゾン変異形成を行った。各変異型クローンは、1つの遺伝子の1つの対立遺伝子が破壊されるヒプロイドハイブリッドであるヘミジゴットであった(図1A、図2)。我々は39°Cでの成長によって互いにヘミツゴテスを競い合い、対照として別々の実験で、28°C(図1B)で、各培養からDNAを分離した。各ヘミジゴットの適合性を報告するために、DNAが断片化され、アダプターにライゲーションされたプロトコルを使用して、バルクシーケンシングを介して豊富さを定量し、続いてトランスポジット挿入位置の増幅を行いました(図1C)。この増幅のプライマーがプロトコルで提供されるものと異なり、効率が悪い場合、バックグラウンド読み取りはシーケンスデータで優勢になり、使用可能な読み取りが少なく...

ディスカッション

以前の統計的遺伝的方法に比べてRH-seqの利点は、数倍です。リンケージおよびアソシエーション分析とは対照的に、RH-seqは単一遺伝子マッピング分解能を提供します。したがって、特定の種の個体間の形質変動の研究や、特異的な違いについても重要な有用性を持つ可能性が高い。また、ゲノム全体の相互異端異端分析における以前の試みは、遺伝子欠失変異体のコレクションを使用し、そのうちのいくつかは偽陽性結果9、10につながる可能性のある二次変異を保有する。RH-seq戦略は、個々の変異クローンの背景が最終的な結果にわずかに寄与するように、各遺伝子に多くのヘミジゴット変異体を順番に生成し、表現することによって、この問題を脇に立たします。原理的には、RH-seqは遺伝子座をコード化しない研究も行っていますが、現在の研究では遺伝子のみに焦点を当てています。

RH-seq、いくつかの生物学的およびいくつかの技術的には、成功した開業医がアプローチの有用性を最大化し、最良の結果への道を加速する...

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
1-2 プラスミド ギガプレップキット	ザイモリサーチ	D4204	必要なキットの数は、各キットのプレップの効率によって異なります。このキットには、5つの個別のプラスミド調製カラムが付属しています。1カラムあたり1 Lを超えると、調製フィルターが目詰まりし、DNAの収量が低く、品質が低下する可能性があるため、各調製カラムに1 Lの飽和大腸菌培養液を流します。
10X Tris-EDTA(TE)バッファー(100 mM Tris-HClおよび10 mM EDTA)	Any	N / A	フィルターは、0.22 μ使用前にmフィルター。
1M LiOAc	任意の	N/A	フィルターは 0.22 μ使用前にmフィルター。
300 mg/mL ジェネチシン (G418)	Gibco	11811023
52% ポリエチレングリコール (PEG) 3350	シグマ	1546547	水に溶解し、0.22 μ;使用前にmフィルター。1X trafo mix: 228 uL 52% PEG, 36 uL 1M LiOAc, 36 uL 10X TE buffer
オートクレーブLB液体ブロス	BD Difco	244620	任意のブランドの粉末とミリクウォーターを使用してLB液体ブロスを作ります。使用前にオートクレーブで
カルベニシリン水(100 mg / mL)	任意の	N / A	フィルターで0.22 μ使用前にmフィルター。
5-フルオロロチン酸(5-FOA)[ウラシルまたは酵母窒素ベース(YNB)を含まない0.2%ドロップアウトアミノ酸ミックス]、0.005%ウラシル、2%D-グルコース、アミノ酸なし0.67%YNB、0.075%5-FOA]	5-FOA:Zymo Research、ドロップアウトミックス:US Biological、Uracil:Sigma、D-グルコース:Sigm)、YNB:Difco	5-FOA:F9001-5、ドロップアウトミックス: D9535、ウラシル:U0750、D-グルコース:G8270、YNB:DF0919
pJR487(CEN-/ARS+piggyBac含有プラスミド)	を
ハイブリッド酵母 JR507 (S. cerevisiae DBVPG1373 x S. paradoxus Z1, URA-/URA-)	該当なし		ム研究室からのお願い
SE-150	イルミナHiseq 2500
可変温度設定の大型振とうインキュベーター	任意の	N/A
LB + カルベニシリン寒天プレート(100 μg/mL)	寒天:BD Difco	寒天:214010	LB寒天プレートを通常通り作り、カルベニシリンを100 μ乾燥前のg / mL。
ナノドロップ分光光度計	Thermo Scientific	ND-2000
量子ビット蛍光計	Thermo Scientific	Q33240
サケ精子DNA	Invitrogen	15632011
39°Cの水浴;C	任意の	N/A
酵母菌 gDNA 調製キット	Zymo Research	D6005
酵母ペプトンデキストロース (YPD) 液体培地	BD Difco	ペプトン: 211677, 酵母 抽出物: 212750	オートクレーブ後にフィルター滅菌した D-グルコースを 2% に加えます。
YPD + G418寒天プレート(300 μg/mL)	寒天:BD Difco	寒天:214010	YPD寒天プレートを通常通り作成し、G418を300 μ乾燥前のg / mL。
YPD寒天プレート	寒天:BD Difco	寒天:214010