この記事では、 すなわち 、テロワールの概念を洞察するために、ベリー品質形質に対する環境の影響をブドウ果実の転写産物および代謝物への非標的メタボロミクス、トランスクリプトミクスと多変量統計解析のアプリケーションについて説明します。
このコンテンツを表示するには、JoVEへの購読が必要です。 または、無料トライアルをお申し込みください。
方法論記事
* These authors contributed equally
この記事では、 すなわち 、テロワールの概念を洞察するために、ベリー品質形質に対する環境の影響をブドウ果実の転写産物および代謝物への非標的メタボロミクス、トランスクリプトミクスと多変量統計解析のアプリケーションについて説明します。
テロワールは、特定の生息地や経営慣行に従って、このようなブドウ( ブドウ )などの作物の特性に影響を与える環境要因の組み合わせを指します。この記事では、特定のテロワールの署名が多変量統計解析を用いて、ブドウ品種のコルビナのベリーメタボロームとトランスクリプトームで検出することができる方法を示しています。この方法は、まず、適切なサンプリング計画が必要です。このケーススタディでは、コルビナ品種の特定のクローンは遺伝的差異を最小限に抑えるように選択し、試料は、3つの異なる栽培期間中に三つの異なるマクロゾーンを表す7ブドウ畑から収集しました。非標的LC-MSメタボロミクスのアプローチはMZmineソフトウェアおよびフラグメンテーションツリー分析に基づいて、代謝物同定の戦略を使用して効率的なデータ処理を伴う、その高感度のために推奨されます。包括的トランスクリプトーム解析は、マイクロアレイを用いて達成することができますすべての〜99%をカバー含むプローブは異なるテロワールのコンテキスト内のすべての示差的に発現する遺伝子の同時分析を可能にする、ブドウの遺伝子を予測しました。最後に、投影法に基づいて、多変量データ解析は、メタボロミクスとトランスクリプトミクスデータを統合し、有益な相関を同定するために詳細に分析することを可能にする、強力なヴィンテージ特有の効果を克服するために使用することができます。
植物のゲノム、トランスクリプトーム、プロテオームおよびメタボロームに基づいて、大規模データ解析は、ブドウの植物とその環境との相互作用を反映したワインのテロワールの特性のような複雑なシステムの挙動に前例のない洞察を提供しています。ワインのテロワールが同一のブドウのクローンが異なるブドウ畑で栽培された場合でも明瞭な可能性があるので、クローンのゲノムが同じであるため、ゲノミクス解析はほとんど役に立ちません。その代わりに、遺伝子発現とワインの品質特性を決定ベリーの代謝特性との相関関係を調べることが必要です。マイクロアレイ上に固定化されたプローブに、そのようなハイブリダイゼーションのような普遍的な特性を利用することによって定量分析を容易にする全ての転写物の類似の化学的特性から、トランスクリプトームの利益のレベルでの遺伝子発現の解析。プロテオミクスaの対照的に、ユニバーサル分析方法NDメタボロミクスが原因個々のタンパク質および代謝産物の巨大な物理的および化学的多様性より挑戦です。メタボロミクスの場合、個々の代謝産物がサイズ、極性、豊かさとボラティリティに大幅に異なるため、この多様性はさらに極端なので、単一の抽出方法や分析方法は、全体的なアプローチを提供しています。
不揮発性代謝物のために適切な分析プラットフォームのうち、高速液体クロマトグラフィーに基づくものを質量分析に連結された(HPLC-MS)は、はるかに敏感例えば紫外線またはダイオードアレイ検出器を用いたHPLC(HPLC-UV、HPLC-DADなどの代替よりも)または核磁気共鳴(NMR)分光法が、HPLC-MSによる定量分析は、マトリックス効果とイオン抑制/強化1-3のような現象に影響され得ます。 (HPLC-Eをエレクトロスプレーイオン化源を用いてHPLC-MSによるコルビナブドウ果実の解析中にこのような効果の調査SI-MS)は、最小の保持時間を有する糖および他の分子を強くおそらくこの領域中の分子の多くを反映し、過少報告されたことを示し、他の分子の存在量は、マトリックス効果によって過大評価または影響を受けず、過小評価することができることしかし、マトリックス効果のためのデータの正規化は、全体の結果4,5への影響は限定的を持っているように見えました。本明細書に記載される方法は、熟成中ブドウ果実において高レベルで蓄積媒体極性代謝物の分析のために最適化され、そしてこれはかなりテロワールによって影響を受けます。彼らは一緒にワインの色、味と健康関連の特性を決定する、アントシアニン、フラボノール、フラバン-3-オール、プロシアニジン、他のフラボノイド、レスベラトロール、スチルベン、ヒドロキシ酸およびヒドロキシ安息香酸が含まれます。 HPLC-MSによる定量の信頼性が低いので、そのような糖や脂肪族有機酸などの他の代謝産物は、原因effec行列に無視されますトンとイオン抑制現象5。この方法により、選択された極性の範囲内に、アプローチは、可能な6限り多くの異なった代謝物を検出することを目的とするという点で非標的です。
ブドウの転写産物の数千人が同時に監視できるようにするトランスクリプトーム法は、完全なブドウのゲノム配列7,8の利用可能性によって促進されます。ハイスループットのcDNA塩基配列決定に基づいて、初期のトランスクリプトミクス手法をまとめ、RNAのSeq技術として記載された手順のコレクションに次世代シーケンサーの登場と進化してきました。このアプローチは、急速にトランスクリプトミクス研究のための最適な方法になってきています。しかし、転写産物の数千人は、ハイブリダイゼーションによって並列に定量化することを可能にするマイクロアレイに基づく文学の大きな体は、ブドウのために蓄積してきました。 RNA-seqが主流の技術になる前に実際に、多くの専用の市販のマイクロアレイプラットフォームをしていましたブドウトランスクリプトームは、非常に詳細に検査することを可能にする開発。プラットフォームの多種多様な中で、2つだけは、ゲノムワイドなトランスクリプトーム解析9を可能にしました。最も進化した配列は、このように各実験のコストを削減し、単一のデバイスに最大12の独立したサンプルのハイブリダイゼーションを可能にしました。 12サブアレイはそれぞれ29549ブドウの転写産物を表す135,000 60-merのプローブを含んでいました。この装置は、試験10-24多数のに使用されています。これら2つのプラットフォームは、現在廃止されましたが、新しいカスタムマイクロアレイは、最近設計されていて、それが追加新たに発見されたブドウの遺伝子25を表すプローブのさらに多くが含まれているとして、より多くの最近の発展を表しています。
トランスクリプトミクスとメタボロミクス解析によって生成された大販売データセットは、別のフォームとの間の相関を決定するために、多変量技術を含むデータ分析に適した統計的手法を必要とし、データの秒。最も広く使用されている多変量技術は、投影に基づくものであり、これらは、主成分分析(PCA)、又は、教師なしとすることができるような潜在構造判別分析(O2PLS-DA)26への双方向直交射影として、監督します。この記事で紹介したプロトコルは、サンプルのグループ間の差異を識別するために、探索的データ解析のためのPCAとO2PLS-DAを使用しています。
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
1.適切な材料を選択し、サンプリング計画を構築
| AM | BA | BM | CS | FA | MN | PM | |
| マクロゾーン | ソアーヴェ | ガルダ湖 | ヴァルポリチェッラ | ガルダ湖 | ヴァルポリチェッラ | ヴァルポリチェッラ | ソアーヴェ |
| 身長(m) | 250 | 120 | 450 | 100 | 130 | 250 | 130 |
| 台木 | 41B | S04 | K5BB | 420A | 420A | K5BB | 41B |
| 行方向 | EW | NS | EW | EW | EW | NS | NS |
| トレーニングシステム | オーバーヘッドシステム(パーゴラ) | オーバーヘッドシステム(パーゴラ) | 垂直シュートポジショニング(ギヨー) | オーバーヘッドシステム(パーゴラ) | Overhe広告システム(パーゴラ) | 垂直シュートポジショニング(ギヨー) | 垂直シュートポジショニング(ギヨー) |
| 土壌型 | シルト質粘土 | ローム | 粘土 | ローム | 埴壌土 | シルトローム | 埴壌土 |
| 植栽のレイアウト(メートル) | 3.20のx 1.00 | 4.50のx 0.80 | 4.00のx 1.25 | 3.50のx 1.20 | 3.50のx 0.75 | 2.80のx 1.00 | 1.80のx 0.80 |
| 総ライム% | 3.9 | 19.3 | 18.3 | 14.4 | 31 | 5.9 | 27.9 |
| アクティブ石灰% | 0.5 | 2.6 | 9.4 | 6.3 | 11.3 | 3.1 | 8.3 |
| 砂 % | 15 | 47 | 66 | 42 | 29 | 13 | 36 |
| ローム% | 43 | 36 | 21 | 37 | 39 | 67 | 36 |
| 粘土 % | 42 | 17 | 13 | 21 | 32 | 20 | 28 |
| 土壌pH | 8.3 | 7.9 | 7.8 | 8.2 | 8.2 | 7.8 | 7.9 |
| 有機物質(%) | 2.9 | 2.5 | 2.2 | 1.2 | 2.9 | 1.6 | 2.5 |
| 交換可能なリン(mgを/キログラム) | 26 | 73 | 73 | 68 | 48 | 47 | 64 |
| 交換性カリウム(ミリグラム/キログラム) | 190 | 376 | 620 | 230 | 168 | 154 | 126 |
| 交換可能なマグネシウム(Mg / kg)を | 272 | 468 | 848 | 623 | 294 | 293 | 183 |
| 交換性カルシウム(ミリグラム/キログラム) | 6500 | 5380 | 7358 | 6346 | 4652 | 10055 | 2878 |
| ベリー還元糖2006 | 211.25±1.20 | 176.20±0.42 | 187.40±0.00 | 203.70±1.13 | 212.55±0.64 | 195.20±0.00 | 211.65±0.64 |
| ベリー還元糖2007 | 190.00±1.27 | 165.25±0.49 | 153.00±0.42 | 203.60±0.71 | 210.90±0.71 | 192.25±0.64 | 188.70±1.84 |
| ベリー還元糖2008 | 191.35±0.64 | 178.90±0.57 | 170.05±0.49 | 205.15±1.48 | 188.70±0.57 | 169.35±0.49 | 108.05±1.06 |
| ベリーのpH 2006 | 3.01±0.01 | 2.96±0.01 | 2.84±0.00 | 2.9±0.00 | 2.98±0.00 | 3.02±0.00 | 3.06±0.01 |
| ベリーのpH 2007 | 2.97±0.00 | 3.00±0.00 | 2.74±0.00 | 3.07±0.01 | 2.98±0.00 | 2.87±0.01 | 3.09±0.00 |
| ベリーのpH 2008 | 2.83±0.00 | 3.04±0.01 | 2.71±0.00 | 2.98±0.01 | 2.98±0.00 | 2.82±0.00 | 3.11±0.00 |
| 収穫日 | 2006 | 2007 | 2008 | ||||
| Veraison | 8 8月 | 18 - 7月 | 12 - 8月 | ||||
| ミッド登熟 | 4 - 9月 | 8 8月 | 2 - 9月 | ||||
| 熟しました | 18 - 9月 | 29 - 8月 | 23 - 9月 |
表 1:各ブドウ畑の主な機能とサンプル収集日付。メートル=メートル、EW =、NS =南北ウエストを食べます。 
図1: サンプリング手順の模式図は、3つのワインの生産マクロゾーンはヴェローナ、ヴェネト地方、イタリアの街の周囲に配置されています。 3つの時点ではブドウ栽培に熟成の開始を表すveraison(V)であり、中期熟成(MR)と熟した果実(R)。 この図の拡大版をご覧になるにはこちらをクリックしてください。
2.ベリーパウダー抽出物を準備し、代謝物を分析し、データを処理
| 質量分析計のコンポーネント | 関数 | パラメーター |
| エレクトロスプレーイオン化ソース | 噴霧ガス | 50 psiで、350°C |
| 乾燥ガス | 10 L分-1 | |
| イオントラップと検出器 | スキャン | フルスキャンモード、13000メートル/秒あたりのz、範囲50-1,500 のm / z |
| ターゲット質量 | 400メートル/ zの | |
| 衝突ガス | ヘリウム | |
| 真空圧力 | 1.4×10 -5ミリバール | |
| 毛細管ソース | 4000 V | |
| エンドプレートオフセット | -500 V | |
| スキマー | -40 V | |
| キャップ出口 | -121 V | |
| 10月1日DC | -12 V | |
| 10月2日DC | -1.7 V | |
| レンズ1 | 5 V | |
| レンズ2 | 60 V | |
| 正イオン化モードのためのICC | 20.000 | |
| 負のイオン化モードのためのICC | 7,000 |
表2:質量スペクトルを取得するための主なパラメータセット。
| 操作 | 選択 | 関数 | パラメーター | 価値観 |
| ピーク検出 | 質量検出 | 重心 | 騒音レベル | 3500 |
| クロマトグラムビルダー | 最高のデータポイント | 分の期間 | 0.15 | |
| 分の高さ | 4,000 | |||
| M / Zトレランス | 0.3 | |||
| ピーク検出 | ピークデコンボリューション | ローカル最小検索 | クロマトグラフのしきい値 | 70 |
| RT範囲内検索の最小値(分) | 0時50分 | |||
| 最小の相対的な高さ | 15% | |||
| 最小の絶対高さ | 4,000 | |||
| ピークトップ/エッジの分率 | 2 | |||
| Durationの範囲(分) | 0-10 | |||
| 同位体 | 同位体ピークのハタ | - | M / Zトレランス | 1.2 |
| RT公差 | 0時50分 | |||
| 単調な形状 | いいえ | |||
| 最大充電 | 3 | |||
| 代表的同位体 | いいえ | |||
| アラインメント | アライナーに参加 | - | M / Zトレランス | 1.2 |
| m / zの重み | 10 | |||
| 保持時間の許容範囲 | 0時50分 | |||
| RTの重み | 5 | |||
| 同じ充電状態を必要とします | いいえ | |||
| 同じIDが必要です | いいえ | |||
| 同位体パターンの比較 | いいえ | |||
| ギャップを埋める | ピークファインダー | - | 強度の許容範囲 | 20% |
| M / Zトレランス | 0.9 | |||
| 保持時間の許容範囲 | 0時40分 | |||
| RT補正 | いいえ | |||
| フィルタリング | ピークフィルタを複製 | - | M / Zトレランス | 1.2 |
| RT公差 | 0時30分 | |||
| 同じ識別を必要とします | いいえ |
表3:負のLC-MSのブドウ果実のデータファイルを処理するために特定の値を持つMzmineワークフロー。
3.トランスクリプトーム解析のためのベリーパウダー抽出物を調製し、データを処理
| メトリック名 | 上限 | 下限 | 説明 |
| AnyColorPrcntFeatNonUnif | 1.00 | NA | いずれかのチャンネルでムラの外れ値を特徴としている機能の割合 |
| DetectionLimit | 2.00 | 0.10 | 平均プラス直線濃度範囲を下回るスパイクインの1標準偏差 |
| absGE1E1aSlope | 1.20 | 0.90 | E1aプローブの濃度対信号のためのフィットの傾きの絶対 |
| MedCVProcSignal | 8.00 | NA | 処理された信号の中央値%のCV |
| gNegCtrlAveBGSubSig | 5.00 | -10.00 | バックグラウンドの平均値は(BGSubSignalが特徴からBGUsedと呼ばれる値をsubstractingによって計算された信号を意味する)すべてのインライヤネガティブコントロールの信号を減算しました |
| gNegCtrlAveNetSig | 40.00 | NA | すべてインライヤ陰性対照の正味の信号の平均 |
| gNegCtrlSDevBGSubSig | 10.00 | NA | の標準偏差背景はすべてインライヤ陰性対照の信号を減算しました |
| gNonCntrlMedCVProcSignal | 8.00 | NA | 非対照プローブの処理された信号のための中央値%のCV |
| gSpatialDetrendRMSFilter | 15.00 | NA | 背景トレンド除去フィットの残留 |
表4:マイクロアレイハイブリダイゼーションの品質を確認するためにチェックされる主なパラメータ。
メタボロミクスとトランスクリプトミクスデータの詳細な統計分析を行う4.
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
この資料に記載されているケーススタディでは、比較的189サンプル間の定量化分子イオンプラスその同位体、付加物およびいくつかの断片を含む552の信号( のm / zの機能)、(7ブドウ畑×3熟成段階×3栽培期間xを含む最終データ行列が得られました3生物学的複製)。データポイントの合計数は、したがって、104328でした。断片化ツリー分析は、代謝産物に加えて付加体、同位体およびその断片に対応する282 のm / zの機能注釈をもたらしました。 PCAにより全データ行列の探索的分析は、熟成の段階に応じてサンプルがクラスタ化されたことが示された(veraison、半ば熟成と熟した)は、第1および第2主成分(T1-T2、 図2A)に沿って、かつ沿って成長する季節に応じて第3主成分(T1〜T3、 図2B)。
SRC = &...
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
この記事では、メタボロミクス、トランスクリプトミクスとブドウ果実のテロワールの概念を解釈するために使用される統計解析プロトコルについて説明します。 HPLC-ESI-MSによるメタボロミクス解析を同時に代謝産物の多くを検出するのに十分な感度であるが、相対的な定量は、マトリックス効果とイオン抑制/増強に影響されます。しかし、同様のアプローチが既にコルビナ果実の熟成及び収穫後の枯死を記述するために使用されており、マトリクス効果の補正が結果5への影響は限定的でした。さらに、同じ試料を用いて、非標的メタボロミクスのためのNMRおよびLC-MSの比較堅牢性を決定するために、最近の大規模なマルチインストゥルメント研究室間の研究では、さまざまな楽器や技術を6一貫した結果が得られたことを明らかにしました。本明細書中に記載のケーススタディでは、3マクロ-Zで7ブドウ園から3熟成の段階でサンプルの収集と分析を関与しました3生育期を超えるもの。ここで、テロワールの署名がPCAとO2PLS-DAの多変量統計を用いて、(3ブドウ畑で表さガルダ湖とソアーヴェ、2ブドウ畑によって表される各、およびヴァルポリチェッラ、)マク...
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
著者らは開示するものは何もない。
この作業は、EUが資金提供したCOST ACTION FA1106 "トマトとブドウの木の果肉質の品質を制御する開発メカニズムを決定するための統合システムアプローチ"内で調整されたネットワーク活動の恩恵を受けました。"この研究は、カリヴェローナ銀行財団が資金提供した「Completamento del Centro di Genomica Funzionale Vegetale」プロジェクトと、イタリア農林政策省が資金提供した「Valorizzazione dei Principali Vitigni Autoctoni Italiani e dei loro Terroir (Vigneto)」プロジェクトの支援を受けました。SDSは、イタリア大学研究省FIRBのRBFR13GHC5プロジェクト「環境相互作用による遺伝子型におけるブドウのエピゲノム可塑性」によって資金提供されました。
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| ミルグラインダー | IKA | IKA A11 ベーシック | |
| HPLC オートサンプラー | ベックマン・コールター | - | System Gold 508 オートサンプラー |
| HPLC システム | ベックマン・コールター | - | System Gold 127 ソルベントモジュール HPLC |
| C18 ガードカラム | グレース | - | Alltima HP C18 (7.5 mm x 2.1 mm; 5 μm) ガードコラム |
| C18 コラム | グレース | - | アルティマ HP C18 (150 mm x 2.1 mm; 3 μm) カラム |
| 質量分析計 | Bruker Daltonics | - | Bruker Esquire 6000;質量分析計にはESIソースが装備されており、分析装置はイオントラップでした。 |
| 抽出溶媒およびHPLCバッファー | Sigma | 34966 | メタノールLC-MSグレード |
| Sigma | 94318 | ギ酸LC-MSグレード | |
| Sigma | 34967 | アセトニトリルLC-MSグレード | |
| Sigma | 39253 | 水 LC-MS グレード | |
| Minisart RC 4 シリンジフィルター (0.2 μm) | Sartorius | 17764 | |
| データ収集(a)および処理のためのソフトウェア(b) | Bruker Daltonics-Bruker | Daltonics Esquire 5.2 Control(a);Esquire 3.2 Data Analysis and MzMine 2.2 software (b) | |
| Spectrum Plant Total RNA kit | Sigma-Aldrich | STRN250-1KT | For total RNA extractino from grape pericarps |
| Nanodrop 1000 | Thermo Scientific | 1000 | |
| BioAnalyzer 2100 | Agilent Technologies | G2939A | |
| RNA 6000 Nano Reagents | Agilent Technologies | 5067-1511 | |
| RNAチップ | Agilent Technologies | 5067-1511 | |
| Agilent Gene Expression Wash Buffer 1 | Agilent Technologies | 5188-5325 | |
| Agilent Gene Expression Wash Buffer 2 | Agilent Technologies | 5188-5326 | |
| LowInput QuickAmp Labeling kit One-Color | Agilent Technologies | 5190-2305 | |
| Kit RNA Spike In - One-Color | AgilentTechnologies | 5188-5282 | |
| 遺伝子発現ハイブリダイゼーションキット | Agilent Technologies | 5188-5242 | |
| RNeasy Mini Kit (50) | Qiagen | 74104 | for cRNA purification |
| Agilent SurePrint HD 4X44K 60-mer Microarray | Agilent Technologies | G2514F-048771 | |
| eArray | Agilent Technologies | - | https://earray.chem.agilent.com/earray/ |
| >Gasket スライド | Agilent Technologies | G2534-60012 | Agilent SurePrint Microarray 4-array Hybridization |
| Thermostatic bath | Julabo-Hybridization | ||
| Chamber | Agilent技術 | G2534-60001 | |
| マイクロアレイハイブリダイゼーションオーブン | アジレントテクノロジー | G2545A | |
| ハイブリダイゼーションオーブン ローテーターラック | アジレントテクノロジー | ズ G2530-60029 | |
| ローテーターラックコンバージョンロッド | アジレントテクノロジー | ズ G2530-60030 | |
| 染色キット | バイオオプティカ | 10-2000 | スライド染色皿とスライドラック |
| マグネチックスターラー装置 | アレックス ヒーティングマグネティックスターラー | F20540163 | |
| サーモスタットオーブン | Thermo Scientific | Heraeus - 6030 | |
| Agilent Microarray Scanner | Agilent Technologies | G2565CA | |
| スキャナーカルーセル、48ポジション | Agilent Technologies | G2505-60502 | |
| スライドホルダー | Agilent Technologies | G2505-60525 | |
| 機能抽出ソフトウェア v11.5 | AgilentTechnologies-Inside | the Agilent Microarray Scanner G2565CA | |
| SIMCA + V13 Software | Umetrics |
アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。