1. セットアップ
<ol><li>R がインストールされていない場合は、R をインストールします。 注: PAST は R で記述されているため、ユーザーには R がインストールされている必要があります。この書き込みの時点で、生物伝導体から直接PASTをインストールするにはR4.0が必要です。 PAST の古いバージョンは、R3.6 のバイオコンダクタからインストールでき、また、R3.5 を使用するユーザーは Github からインストールできます。R のインストール手順は、次のリンクからダウンロードできます https://www.r-project.org/。</li>
	<li>RStudio デスクトップの最新バージョンをインストールするか、RStudio を更新します (省略可能)。 注: RStudio は、R 言語を使用する場合に役立つ環境です。特に、シンニー GUI アプリケーションを使用するのではなく、コマンド ラインで PAST を実行することを選択したユーザーには、インストールすることをお勧めします。RStudio とそのインストール手順は、次のリンク https://rstudio.com/products/rstudio/<a href="https://rstudio.com/products/rstudio/">にあります。</a></li>
	<li>生体伝導体11 から、生体伝導体の指示に従ってPASTをインストールします。 注: 生体伝導体によるインストールは、PAST の依存関係のインストールを処理する必要があります。さらに、PAST は Github12からインストールできますが、Github からインストールすると依存関係は自動的にインストールされません。</li>
	<li>過去のシャイニー(オプション)をインストールします。ファイル"アプリ"をダウンロードします。Github リポジトリの [リリース] ページの R" https://github.com/IGBB/PAST/releases/、ダウンロードしたファイルの場所を覚えておいてください。 注: PAST は R でメソッドを直接呼び出すことによって使用できますが、R にあまり慣れていないユーザーは、ガイド付きユーザー インターフェイスを提供する PAST シャイニー アプリケーションを実行できます。PAST シャイニーは、過去の Github リポジトリのshiny_appブランチで利用できる R スクリプトです。PAST シャイニーは、最初の実行時にその依存関係をインストールしようとします。</li>
	<li>次に説明する 3 つの方法のいずれかで、アプリケーションを開始して分析を開始します。<ol>
<li>RStudioと過去の光沢<ol>
<li>RStudio を使用して、アプリのフォルダーに新しいプロジェクトを作成します。R が配置されます。[ ファイル] |新しいプロジェクト をクリックし、そのフォルダを選択します。</li>
			<li>新しいプロジェクトが作成されたら、アプリを開きます。R ファイルは以前にダウンロードされました。RStudio はそのアプリを認識します。R は、シャイニー アプリで、表示されたソース コードの上のバーに [アプリの実行 ] ボタンを作成します。[ アプリケーションの実行] をクリックします。その後、RStudio は PAST のシャイニー アプリケーションを表示するウィンドウを起動します。</li>
		</ol>
</li>
		<li>Rコンソールで過去の光沢<ol>
<li>Rを起動し、PASTシャイニーアプリケーションを起動するために次のコードを実行します: 光沢::runApp('パス/フォルダ/ウィズ/シャイニー/アプリ。R'。引用符で囲まれたテキストを、アプリのフォルダーに置き換えます。Rがダウンロードされ、引用符を保持します。</li>
		</ol>
</li>
		<li>Rシャイニーなしの過去<ol>
<li>ライブラリ (PAST)を Rコンソールで実行して、PASTをロードします。</li>
		</ol>
</li>
	</ol>
</li>
</ol>2. シャイニー分析のカスタマイズ (オプション)
<ol><li>解析のタイトルを「新規分析」から、複数の分析を追跡するのに役立つ、実行されている分析のタイプをよりよく反映したものに変更します ( 図 1を参照)。</li>
</ol><img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/61268/61268fig1.jpg"> 図1. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig1large.jpg" target="_blank">ここをクリックして、この図の大きなバージョンを表示してください。</a>
<ol start="2"><li>コアの数とモードを変更します。コアの数を 1 からコンピューターの合計数の間の任意の数に設定しますが、PAST にリソースを追加すると、マシン上の他の操作が遅くなる可能性があることに注意してください。セクション 6 の説明に基づいてモードを設定します。</li>
</ol>3. GWAS データの読み込み
注: GWAS データがタブ区切りであることを確認します。関連ファイルに、特性、マーカー名、軌跡または染色体、染色体上の位置、p値、およびマーカーの R2 値が含まれていることを確認します。エフェクトファイルに、特性、マーカー名、軌跡または染色体、染色体上の位置、および効果の列が含まれていることを確認します。ユーザーはデータを読み込むときに列の名前を指定できるため、これらの列の順序は重要ではありません。追加の列は無視されます。TASSEL13 は、これらのファイルを生成するために使用できます。
<ol><li>過去の光沢のあるGWASデータをロードします。<ol>
<li>[関連付けファイル] および [エフェクト ファイル] 選択ボックスを使用して、 関連付けファイル と エフェクト ファイル を選択します。ファイル選択ボックスの下にある [関連列名 ]および[ エフェクト列名] 入力ボックスの列名を変更して、データ内の列名を反映します。</li>
	</ol>
</li>
</ol><img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/61268/61268fig2.jpg"> 図2. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig2large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
<ol start="2"><li>R コンソールで PAST を使用して GWAS データをロードします。<ol>
<li>次のコードを変更して実行します。 gwas_data =load_GWAS_data("path/to/association_file.tsv"、"パス effects_file/association_file.tsv"、association_columns=c(「形質」、「マーカー」、「ローカス」、「サイト」、「サイト」、「p」、「marker_R2」)、effects_columns=c(「形質」、「マーカー」、「ローカ」、「サイト」、「効果」)
</li>
	</ol>
</li>
	<li>注: GWAS ファイルの実際の場所へのパスを変更します。association_columnsとeffects_columnsに指定された値は、既定値です。名前がデフォルト値と一致しない場合は、列名を指定します。それ以外の場合は省略できます。</li>
</ol>4. リンケージの不一定衡 (LD) データの読み込み
注: リンケージの不衡 (LD) データがタブ区切りで、次の種類のデータが含まれていることを確認します: 軌跡、位置 1、サイト 1、位置 2、サイト 2、位置 1 と位置 2 の間のベース ペアの距離、および R2 値。
<ol><li>過去のシャイニーで LD データをロードします。<ol>
<li>LD データを含むファイルを選択します。必要に応じて 、LD データの列名 と一致するように、ファイル選択ボックスの下にある [LD 列名] 入力ボックスの列名を変更します。</li>
	</ol>
</li>
</ol><img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/61268/61268fig3.jpg"> 図 3. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig3large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
<ol start="2"><li>LD データを PAST を使用して R コンソールにロードします。<ol>
<li>LD データを読み込むには、次のコードを変更して実行します。 LD = load_LD("パス/先/LD.tsv"、LD_columns=c(「ローカス1」、「位置1」、「サイト1」、「位置2」、「サイト2」、「Dist_bp」、「R.2」) メモ: LD ファイルの実際の場所へのパスを変更します。LD_columnsに指定された値は、デフォルト値です。名前がこれらのデフォルト値と一致しない場合は、列の正しい名前を指定します。それ以外の場合は省略できます。</li>
	</ol>
</li>
</ol>5. 遺伝子へのスナップの割り当て
注: GFF 形式の注釈をダウンロードまたは検索します。これらの注釈は、多くの場合、特定の生物のためのオンラインデータベースで見つけることができます。注釈データの品質が経路解析の品質に影響を与えるため、低品質の注釈については注意が必要です。これらの注釈の最初の列(染色体)が、関連付け、効果、LD データの軌跡/染色体の形式と一致することを確認します。たとえば、GWAS と LD データ ファイルが最初の染色体 "1" を呼び出す場合、注釈は最初の染色体 "chr1" を呼び出す必要があります。
<ol><li>PASTシャイニーを使用して遺伝子にスナップを割り当てます。 注: 適切な R2 カットオフの決定に関する詳細については、Tang et al.6の「経路解析のための SNP to 遺伝子アルゴリズム」のセクションで見つけることができます。<ol>
<li>GFF 注釈を含むファイルを選択します。考慮する種に最適なウィンドウ サイズと R2 カットオフを検討し、デフォルトがアップロードされたデータに適合しない場合は変更します。 注: PAST のデフォルト値は、主にトウモロコシに適した値を反映します。このステップでは、PAST シャイニー解析の最初に設定されたコアの数 (ステップ 2.2) を使用します。</li>
	</ol>
</li>
</ol><img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/61268/61268fig4.jpg"> 図 4. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig4large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
<ol start="2"><li>R コンソールで PAST を使用して、SNPs を遺伝子に割り当てます。<ol>
<li>次のコードを変更して実行し、SNPs を遺伝子に割り当てます。 遺伝子 = assign_SNPs_to_genes(gwas_data、LD、「パス/パス/注釈.gff」、c(「遺伝子」)、1000、0.8、2) 注: このサンプル コードでは、いくつかの既定の提案が提供されています: 1000 は、遺伝子を検索する SNP の周囲のウィンドウのサイズです。0.8 は R2のカットオフ値です。2 は、並列処理に使用されるコアの数です。注釈へのパスも、注釈ファイルの実際の場所に変更する必要があります。</li>
	</ol>
</li>
</ol>6. 重要な経路を発見する
注: 経路ファイルに、各経路の各遺伝子に対して 1 行のデータがタブ区切り形式で含まれていることを確認します。経路記述 - 経路が「トランスリコピン生合成」などの何をするかのより長い説明;gene - アノテーションで提供される名前と一致する経路の遺伝子。経路情報は、MaizeGDB などの特定の生物のオンライン データベースで見つかる可能性があります。2 番目のユーザー指定オプションはモードです。「増加」とは、歩留まりなど測定された形質の増加値が望ましい場合に反射する表現型を指し、「減少」とは、昆虫の損傷評価などの測定値の減少が有益である形質を指す。経路の重要性は、前に説明した方法4,6,14を使用してテストされます。
<ol><li>PASTシャイニーで重要な経路を発見してください。<ol>
<li>パス データを含むファイルを選択し、解析オプションでモードが選択されていることを確認します。必要に応じて、分析のためにそれを保持するために経路に必要な遺伝子の数と、効果の有意性をテストするためにヌル分布を作成するために使用される順列の数を変更します。</li>
	</ol>
</li>
</ol><img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/61268/61268fig5.jpg"> 図 5. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig5large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
注: このステップでは、PAST シャイニー解析の開始時に設定されたコア数とモードセットを使用します。遺伝子のデフォルト数は現在5遺伝子に設定されているので、既知の遺伝子が少ない経路は除去されます。ユーザーはこの値を 4 または 3 に下げて短い経路を含めることができますが、これを行うと誤検出の結果が生じる危険性があります。この値を大きくすると、解析の力が増えますが、解析からより多くの経路が除去されます。使用される順列の数を変更すると、テストの能力が増加し、減少します。
<ol start="2"><li>RコンソールでPASTを使用して重要な経路を発見してください。<ol>
<li>重要な経路を検出するには、次のコードを変更して実行します。 rugplots_data <-find_pathway_significance(遺伝子、「パス/パス/パス/パス.tsv」、5、「増加」、1000、2) 注: このサンプル コードでは、推奨される既定値がいくつか用意されています。5は、分析中の経路を維持するために経路内になければならない遺伝子の最小数であり、増加は測定された形質の増加量を指し(ユーザーは形質に関係なく増加と減少の両方を実行することが推奨され、データ解釈は2つの場合は異なるが、1000はヌル分布を決定する効果をサンプリングする回数である。 2 は並列処理に使用されるコアの数です。パスをパスパスファイルの実際の場所に変更します。</li>
	</ol>
</li>
</ol>7. ラグプロットを表示する
<ol><li>過去のシャイニーでラグプロットを表示します。<ol>
<li>すべての入力がアップロードされ、設定されたら、[ 分析の開始] をクリックします。進行状況バーが表示され、解析のどのステップが最後に完了したかを示します。解析が完了すると、PAST シャイニーが [結果] タブに切り替わります。結果の表が左側の列(「経路」とラベル付け)に表示され、ラグプロットが右側の列(「プロット」とラベル付け)に表示されます。</li>
		<li>スライダーを使用して、フィルター処理パラメーターを制御します。フィルタリングレベルが十分な場合は、左下の [結果のダウンロード ]ボタンをクリックして、すべての画像とテーブルを、解析タイトルで指定されたZIPファイルに個別にダウンロードします。この ZIP ファイルには、フィルター処理されたテーブル、フィルタ処理されていないテーブル、およびフィルター選択されたテーブルのパスごとに 1 つのイメージが含まれています。</li>
	</ol>
</li>
</ol><img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/61268/61268fig6.jpg"> 図 6. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig6large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/61268/61268fig7.jpg"> 図 7. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig7large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>
<ol start="2"><li>R コンソールで PAST を使用してラグプロットを表示する<ol>
<li>次のコードを変更して実行し、結果を保存します。 plot_pathways(rugplots_data、"pvalue"、0.02、「増加」、「output_folder」) 注: このサンプル コードでは、推奨される既定値がいくつか用意されています。pvalue は、有意性しきい値がユーザーによって選択された後に、重要でない経路をフィルタリングするために使用できるデータを提供します。0.02 はフィルター処理で使用される既定値で、増加とは、測定された特性の増加量を指します (ユーザーは、特性に関係なく、増加と減少の両方を実行することをお勧めします。output_folderは、イメージとテーブルが書き込まれるフォルダです (このフォルダは、関数を実行する前に存在している必要があります)。フィルター処理された結果、フィルター処理されていない結果、およびフィルター処理された結果内のすべての経路の個々のイメージの表がこのフォルダーに書き込まれます。</li>
	</ol>
</li>
</ol>

<ol>
	<li>Rafalski, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Association+genetics+in+crop+improvement.">Association genetics in crop improvement.</a> Current Opinion in Plant Biology. 13 (2), 174-180 (2010).</li><li>Yan, J., Warburton, M., Crouch, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Association+Mapping+for+Enhancing+Maize+(Zea+mays+L.)+Genetic+Improvement.">Association Mapping for Enhancing Maize (Zea mays L.) Genetic Improvement.</a> Crop Science. 51 (2), 433-449 (2011).</li><li>Xiao, Y., Liu, H., Wu, L., Warburton, M., Yan, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome-wide+Association+Studies+in+Maize:+Praise+and+Stargaze.">Genome-wide Association Studies in Maize: Praise and Stargaze.</a> Molecular Plant. 10 (3), 359-374 (2017).</li><li>Wang, K., Li, M., Bucan, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pathway-Based+Approaches+for+Analysis+of+Genomewide+Association+Studies.">Pathway-Based Approaches for Analysis of Genomewide Association Studies.</a> The American Journal of Human Genetics. 81 (6), 1278-1283 (2007).</li><li>Weng, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SNP-based+pathway+enrichment+analysis+for+genome-wide+association+studies.">SNP-based pathway enrichment analysis for genome-wide association studies.</a> BMC Bioinformatics. 12 (1), 99 (2011).</li><li>Tang, J., Perkins, A., Williams, W., Warburton, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Using+genome-wide+associations+to+identify+metabolic+pathways+involved+in+maize+aflatoxin+accumulation+resistance.">Using genome-wide associations to identify metabolic pathways involved in maize aflatoxin accumulation resistance.</a> BMC Genomics. 16 (1), 673 (2015).</li><li>Warburton, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome-Wide+Association+Mapping+of+Aspergillus+flavus+and+Aflatoxin+Accumulation+Resistance+in+Maize.">Genome-Wide Association Mapping of Aspergillus flavus and Aflatoxin Accumulation Resistance in Maize.</a> Crop Science. 55 (5), 1857-1867 (2015).</li><li>Warburton, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Genome-Wide+Association+and+Metabolic+Pathway+Analysis+of+Corn+Earworm+Resistance+in+Maize.">Genome-Wide Association and Metabolic Pathway Analysis of Corn Earworm Resistance in Maize.</a> The Plant Genome. 11 (1), 170069 (2018).</li><li>Li, H., Thrash, A., Tang, J., He, L., Yan, J., Warburton, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Leveraging+GWAS+data+to+identify+metabolic+pathways+and+networks+involved+in+maize+lipid+biosynthesis.">Leveraging GWAS data to identify metabolic pathways and networks involved in maize lipid biosynthesis.</a> The Plant Journal. 98 (5), 853-863 (2019).</li><li>Thrash, A., Tang, J., DeOrnellis, M., Peterson, D., Warburton, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=PAST:+The+Pathway+Association+Studies+Tool+to+Infer+Biological+Meaning+from+GWAS+Datasets.">PAST: The Pathway Association Studies Tool to Infer Biological Meaning from GWAS Datasets.</a> Plants. 9 (1), 58 (2020).</li><li>Adam, T., Mason, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> PAST: Pathway Association Study Tool (PAST). Bioconductor version: Release (3.10). , (2020).</li><li>Thrash, A., DeOrnellis, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> IGBB/PAST. , (2019).</li><li>Bradbury, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=TASSEL:+software+for+association+mapping+of+complex+traits+in+diverse+samples.">TASSEL: software for association mapping of complex traits in diverse samples.</a> Bioinformatics. 23 (19), 2633-2635 (2007).</li><li>Subramanian, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Gene+set+enrichment+analysis:+A+knowledge-based+approach+for+interpreting+genome-wide+expression+profiles.">Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles.</a> Proceedings of the National Academy of Sciences U.S.A. 102, 15545-15550 (2005).</li></ol>

著者らは開示するものは何もない。

PASTの主な目標は、GWASデータの代謝経路分析を、特に非ヒトおよび非動物生物に対してより多くの聴衆に持ち込むことである。PAST の代替方法は、多くの場合、人間や動物に焦点を当てたコマンドラインプログラムです。ユーザーフレンドリさは、光沢のあるアプリケーションを開発することを選択することと、Rとバイオコンダクターを使用してアプリケーションをリリースするという点の両方で、PASTの開発における主要な目標でした。ユーザーは PAST を使用するためにプログラムをコンパイルする方法を学ぶ必要はありません。
ほとんどのタイプの解析ソフトウェアと同様に、PAST の結果は入力データと同じ結果にすぎません。入力データにエラーがある場合や、形式が誤っている場合、PAST は実行に失敗するか、または情報が得られない結果を生成します。GWAS データ、LD データ、アノテーション、およびパス・ファイルが正しくフォーマットされていることを確認することは、PAST から正しい出力を受け取るうえで重要です。PAST は、バイアラリック マーカーのみを分析し、入力データのセットごとに 1 つの特性のみを実行できます。さらに、遺伝子型入力が不十分、誤った、または不正確な型定義によって生成されたGWASデータは、明確または再現可能な結果を生み出す可能性は低い。PASTはGWASの結果の生物学的解釈に役立ちますが、環境変動、実験的エラー、または人口構造が適切に考慮されなかった場合、混沌としたデータセットを明確にすることはほとんどありません。
ユーザーは、シャイニーアプリケーションと、それらのパラメータをRコンソールのPASTの関数に渡すことによって、分析の一部のパラメータを変更することができます。これらのパラメータは PAST によって報告される結果を変更する可能性があり、ユーザーはデフォルトからこれらの結果を変更する際に注意する必要があります。LDはユーザによって測定されるので、通常はGWASでも使用されていたのと同じマーカーデータセットを使用するため、LD測定値は母集団に固有である。すべての研究、特にトウモロコシ以外の種(特に自己受粉、多発性、または非常に異種種)については、デフォルトの変化が保証される可能性があります。

An erratum was issued for: A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information. One of the affiliations was updated.
The second affiliation was updated from:
USDA-ARS Corn Host Plant Resistance Research Unit, Mississippi State University
to:
Corn Host Plant Resistance Research Unit, USDA-ARS

An erratum was issued for:&#160;<a href="https://www.jove.com/t/61268">A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information</a>. One of the affiliations was updated.

Erratum: A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information

ゲノムワイド関連研究(GWAS)は、複雑な形質とそれらに関連するゲノム領域を1、2、3に関連するゲノム領域を研究する一般的な方法である。このタイプの研究では、数十万個の一塩基多型(SNP)マーカーが形質との関連についてテストされ、関連の有意性が評価される。誤検出率 (FDR) のしきい値 (またはその他の種類の有意性しきい値) を満たすマーカー特性アソシエーションは、調査のために保持されますが、真の関連付けは除外される場合があります。複雑なポリジェニック形質の場合、各遺伝子の効果は小さく(したがって除外される)可能性があり、一部の対立遺伝子は研究3に存在しないかもしれない特定の条件でのみ発現される。したがって、多くのSNPsは形質に関連付けられているように保持されるかもしれませんが、それぞれが非常に小さな効果を有する可能性があります。あまりにも多くのSNP呼び出しが欠落し、形質の生物学的意味と遺伝的アーキテクチャの解釈が不完全で混乱する可能性があります。代謝経路解析は、生物学的機能4、5、6に従ってグループ化された遺伝子の複合効果に焦点を当てることで、これらの問題のいくつかに対処するのに役立ちます。
この記事で説明したメソッドの以前の実装を使用して、いくつかの調査が完了しました。アフラトキシン蓄積7、コーンミミズ耐性8、及び油生合成9 はいずれも、以前の実施で検討した。これらの解析は成功しましたが、分析ツールは R、Perl、Bash の組み合わせで作成され、パイプラインは自動化されていないため、分析プロセスは複雑で時間がかかり、面倒でした。解析ごとにこの方法を変更するために必要な専門的な知識のため、他の研究者と共有できる新しい方法が開発されました。
経路関連スタディツール(PAST)10 は、プログラミング言語の知識を減らし、より短い期間で分析を実行することで、以前の方法の欠点に対処するように設計されました。この方法はトウモロコシでテストされましたが、PASTは種固有の仮定をしません。PASTは、シャイニーアプリとしてRコンソールを介して実行することができ、オンラインバージョンはまもなくMaizeGDBで利用可能になると予想されます。

<table border="0" cellpadding="0" cellspacing="0" style="width:719px;" width="717">
	<tbody>
		<tr height="21">
			<td height="21" style="height:21px;width:223px;">Computer</td>
			<td style="width:74px;">NA</td>
			<td style="width:123px;">NA</td>
			<td style="width:300px;">Any computer with 8GB RAM should be sufficient</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:223px;">R</td>
			<td style="width:74px;">R Project</td>
			<td style="width:123px;">NA</td>
			<td>R 3.6 or greater is required to install from Bioconductor</td>
		</tr>
	</tbody>
</table>

a pathway association study tool for gwas analyses of metabolic pathway information

近年、代謝経路解析を用いたゲノムワイド関連解析(GWAS)データの解釈に関する新たな実装が開発され、発表された。経路協会スタディツール(PAST)は、ユーザーの使いやすさと遅い分析の懸念に対処するために開発されました。この新しいユーザーフレンドリーなツールは、バイオコンダクターとGithubでリリースされました。テストでは、PASTは以前は24時間以上必要な分析を1時間以内に実行しました。この記事では、書き込みアプリケーションまたは R コンソールを使用して PAST を実行するためのプロトコルを紹介します。

PAST ソフトウェアツールの実行後に結果が生成されない場合は、すべての入力ファイルが正しくフォーマットされていることを確認してください。図 8に、グレイン カラーのトウモロコシ GWAS に基づく PAST パッケージのサンプル データを使用した正常な実行を示します。この表と結果の画像は、[結果のダウンロード] ボタンを使用してダウンロードできます。ダウンロードした画像の例を 図 210に示します。誤った設定は生物学的に意味をなさない結果につながる可能性がありますが、誤りを判断することは研究者次第であり、選択した設定の妥当性を再確認し、関心のある特性に関するすべての既知の証拠を考慮する必要があります。
図910 は、粒色に対してフェノタイプされた288本の近交系ラインのトウモロコシパネルで作成されたGWAS結果の経路解析から生成されたラグプロットを示す。この単純な例は、表現型が「白」または「黄色」であったが、明るい黄色のカロテノイド顔料を作成する経路が知られており、表現型の大部分を担うべきであるため、使用された。したがって、トランスリコピン生合成経路(カロテノイドを産生する)が穀物色と有意に関連することを期待した。パス ID と名前は、グラフの上部に表示されます。グラフの横軸は、分析に含まれていたすべての遺伝子を、形質に対する最大の効果から最小の順に左から右に並べ替えます。しかし、トランスリコピン生合成経路の遺伝子のみがマークされています(グラフの上部には、ハッチマークとして、分析中の他のすべての遺伝子と比較してその効果の遺伝子ランクに現れます)。この経路には7つの遺伝子があります。ランニングエンリッチメントスコア(ES)は、垂直軸に沿ってプロットされます。各遺伝子のESは効果の順序で走行合計に加え、合計は分析された遺伝子の数に合わせて調整されます。したがって、スコアは水平軸に沿って右に移動するにつれて変化し、より大きな効果遺伝子が含まれるにつれて増加する傾向がありますが、ある時点で、効果の増加は別の遺伝子を追加した場合の調整よりも小さく、スコア全体が減少し始めます。実行中の ES 線の頂点は、垂直点線で示されます。これは経路全体のESであり、経路が選択され、ラグプロットとして提示されているかどうかを判断するためにプログラムによって使用されます。
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/61268/61268fig8.jpg"> 図8:PASTシャイニーの実行が完了しました。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig8large.jpg" target="_blank"></a>
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/61268/61268fig9.jpg"> 図9:PASTの実行が完了した(またはシャイニーからダウンロードされた)からの経路画像。 この図は、スラッシュら10から引用されています。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig9large.jpg" target="_blank">この図の大きなバージョンを表示するには、ここをクリックしてください。</a>

Watch this Scientific Journal Video about A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information at JoVE.com

A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information

シャイニーアプリケーションまたはRコンソールを通じて経路協会研究ツール(PAST)を実行することにより、研究者は関連する代謝経路を調査することによって、ゲノム全体の関連研究(GWAS)の結果の生物学的意味をより深く理解することができます。

代謝経路情報のGWAS解析のための経路関連研究ツール

Recently, a new implementation of a previously described method for interpreting genome-wide association study (GWAS) data using metabolic pathway ...

a-pathway-association-study-tool-for-gwas-analyses-metabolic-pathway

Research

JoVE Journal

Genetics

3.1K ビュー.  Mississippi State University. シャイニーアプリケーションまたはRコンソールを通じて経路協会研究ツール(PAST)を実行することにより、研究者は関連する代謝経路を調査することによって、ゲノム全体の関連研究(GWAS)の結果の生物学的意味をより深く理解することができます。 

ビデオ: A Pathway Association Study Tool for GWAS Analyses of Metabolic Pathway Information

Ultralow Input Genome Sequencing Library Preparation from a Single Tardigrade Specimen

Tardigrades are microscopic animals that enter an ametabolic state called anhydrobiosis when facing desiccation and can return to their original state when water is supplied. The genomic sequencing of microscopic animals such as tardigrades risks bacterial contamination that sometimes leads to erroneous interpretations, for example, regarding the extent of horizontal gene transfer in these animals. Here, we provide an ultralow input method to sequence the genome of the tardigrade, Hypsibius dujardini, from a single specimen. By employing rigorous washing and contaminant exclusion along with an efficient extraction of the 50 ~ 200 pg genomic DNA from a single individual, we constructed a library sequenced with a DNA sequencing instrument. These libraries were highly reproducible and unbiased, and an informatics analysis of the sequenced reads with other H. dujardini genomes showed a minimal amount of contamination. This method can be applied to unculturable tardigrades that could not be sequenced using previous methods.

Contamination during the genomic sequencing of microscopic organisms remains a large problem. Here, we show a method to sequence the genome of a tardigrade from a single specimen with as little as 50 pg of genomic DNA without whole genome amplification to minimize the risk of contamination.

超低入力ヒトゲノム単一 Tardigrade の標本からのライブラリの準備

Tardigrades are microscopic animals that enter an ametabolic state called anhydrobiosis when facing desiccation and can return to their original state ...

遺伝学

Methodology for Accurate Detection of Mitochondrial DNA Methylation

Quantification of DNA methylation can be achieved using bisulfite sequencing, which takes advantage of the property of sodium bisulfite to convert unmethylated cytosine into uracil, in a single-stranded DNA context. Bisulfite sequencing can be targeted (using PCR) or performed on the whole genome and provides absolute quantification of cytosine methylation at the single base-resolution. Given the distinct nature of nuclear- and mitochondrial DNA, notably in the secondary structure, adaptions of bisulfite sequencing methods for investigating cytosine methylation in mtDNA should be made. Secondary and tertiary structure of mtDNA can indeed lead to bisulfite sequencing artifacts leading to false-positives due to incomplete denaturation poor access of bisulfite to single-stranded DNA. Here, we describe a protocol using an enzymatic digestion of DNA with BamHI coupled with bioinformatic analysis pipeline to allow accurate quantification of cytosine methylation levels in mtDNA. In addition, we provide guidelines for designing the bisulfite sequencing primers specific to mtDNA, in order to avoid targeting undesirable NUclear MiTochondrial segments (NUMTs) inserted into the nuclear genome.

Here, we present a protocol to allow accurate quantification of mitochondrial DNA (mtDNA) methylation. In this protocol, we describe an enzymatic digestion of DNA with BamHI coupled with a bioinformatic analysis pipeline which can be used to avoid overestimation of mtDNA methylation levels caused by the secondary structure of mtDNA.

ミトコンドリア DNA のメチル化の正確な検出のための方法論

Quantification of DNA methylation can be achieved using bisulfite sequencing, which takes advantage of the property of sodium bisulfite to convert ...

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture 4C-seq

The identification of physical interactions between genes and regulatory elements is challenging but has been facilitated by chromosome conformation capture methods. This modification to the 4C-seq protocol mitigates PCR bias by minimizing over-amplification of PCR templates and maximizes the mappability of reads by incorporating an addition restriction enzyme digest step.

円の染色体構造の次世代シーケンシングによる遺伝子制御配列の高スループット識別キャプチャ (4 C seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and ...

Chromogenic In Situ Hybridization as a Tool for HPV-Related Head and Neck Cancer Diagnosis

Human papillomavirus (HPV) infection is a major risk factor for a subtype of oropharyngeal squamous cell carcinoma (OPSCC), which tends to be associated with a better outcome than alcohol- and tobacco-related OPSCC. Chromogenic in situ hybridization (CISH) of HPV viral RNA could allow the semiquantitative evaluation of viral transcripts of the oncogenic proteins E6 and E7 and an in situ visualization with a good spatial resolution. This technique allows the diagnosis of an active infection with the visualization of HPV transcription in the tumoral HPV-infected cells. An advantage of this technique is the avoidance of contamination from nonneoplastic HPV-infected cells adjacent to the tumor. Overall, its good diagnosis performances have it considered to be the gold standard for active HPV infection identification. Since E6 and E7 viral protein interaction with cell proteins pRb and p53 is mandatory for cell transformation, HPV RNA CISH is functionally relevant and acutely reflects active oncogenic HPV infection. This technique is clinically relevant as well since &#34;low&#34; or &#34;high&#34; HPV transcription levels helped the identification of two prognosis groups among HPV-related p16-positive head and neck cancer patients. Here we present the protocol for manual HPV RNA CISH performed on formalin-fixed paraffin-embedded (FFPE) slides with a kit obtained from the manufacturer. Instead of chromogenic revelation, RNA in situ hybridization may also be performed with fluorescent revelation (RNA FISH). It may also be combined with conventional immunostaining.

Human papillomavirus (HPV) RNA chromogenic in situ hybridization is considered to be one of the gold standards for active human papillomavirus infection detection within tumors. It allows the visualization of HPV E6-E7 mRNA expression with localization and semiquantitative evaluation of its signal.

発色性の発色性ハイブリダイゼーションは、HPV関連の頭頸部癌診断のためのツールを提供する

Human papillomavirus (HPV) infection is a major risk factor for a subtype of oropharyngeal squamous cell carcinoma (OPSCC), which tends to be associated ...

Navigating MARRVEL, a Web-Based Tool that Integrates Human Genomics and Model Organism Genetics Information

Through whole-exome/genome sequencing, human geneticists identify rare variants that segregate with disease phenotypes. To assess if a specific variant is pathogenic, one must query many databases to determine whether the gene of interest is linked to a genetic disease, whether the specific variant has been reported before, and what functional data is available in model organism databases that may provide clues about the gene&#8217;s function in human. MARRVEL (Model organism Aggregated Resources for Rare Variant ExpLoration) is a one-stop data collection tool for human genes and variants and their orthologous genes in seven model organisms including in mouse, rat, zebrafish, fruit fly, nematode worm, fission yeast, and budding yeast. In this Protocol, we provide an overview of what MARRVEL can be used for and discuss how different datasets can be used to assess whether a variant of unknown significance (VUS) in a known disease-causing gene or a variant in a gene of uncertain significance (GUS) may be pathogenic. This protocol will guide a user through searching multiple human databases simultaneously starting with a human gene with or without a variant of interest. We also discuss how to utilize data from OMIM, ExAC/gnomAD, ClinVar, Geno2MP, DGV and DECHIPHER. Moreover, we illustrate how to interpret a list of ortholog candidate genes, expression patterns, and GO terms in model organisms associated with each human gene. Furthermore, we discuss the value protein structural domain annotations provided and explain how to use the multiple species protein alignment feature to assess whether a variant of interest affects an evolutionarily conserved domain or amino acid. Finally, we will discuss three different use-cases of this website. MARRVEL is an easily accessible open access website designed for both clinical and basic researchers and serves as a starting point to design experiments for functional studies.

Here, we present a protocol to access and analyze many human and model organism databases efficiently. This protocol demonstrates the use of MARRVEL to analyze candidate disease-causing variants identified from next-generation sequencing efforts.

人間のゲノムとモデル生物遺伝学情報を統合するWebベースのツール、MARRVELのナビゲート

Through whole-exome/genome sequencing, human geneticists identify rare variants that segregate with disease phenotypes. To assess if a specific variant is ...

Mass Spectrometry-Based Proteomics Analyses Using the OpenProt Database to Unveil Novel Proteins Translated from Non-Canonical Open Reading Frames

Genome annotation is central to today's proteomic research as it draws the outlines of the proteomic landscape. Traditional models of open reading frame (ORF) annotation impose two arbitrary criteria: a minimum length of 100 codons and a single ORF per transcript. However, a growing number of studies report expression of proteins from allegedly non-coding regions, challenging the accuracy of current genome annotations. These novel proteins were found encoded either within non-coding RNAs, 5' or 3' untranslated regions (UTRs) of mRNAs, or overlapping a known coding sequence (CDS) in an alternative ORF. OpenProt is the first database that enforces a polycistronic model for eukaryotic genomes, allowing annotation of multiple ORFs per transcript. OpenProt is freely accessible and offers custom downloads of protein sequences across 10 species. Using OpenProt database for proteomic experiments enables novel proteins discovery and highlights the polycistronic nature of eukaryotic genes. The size of OpenProt database (all predicted proteins) is substantial and need be taken in account for the analysis. However, with appropriate false discovery rate (FDR) settings or the use of a restricted OpenProt database, users will gain a more realistic view of the proteomic landscape. Overall, OpenProt is a freely available tool that will foster proteomic discoveries.

OpenProt is a freely accessible database that enforces a polycistronic model of eukaryotic genomes. Here, we present a protocol for the use of OpenProt databases when interrogating mass spectrometry datasets. Using OpenProt database for analysis of proteomic experiments allows for discovery of novel and previously undetectable proteins.

非正規の開いたリーディング ・ フレームから翻訳された蛋白質を発表する OpenProt データベースを利用した質量分析を用いたプロテオミクス解析

Genome annotation is central to today's proteomic research as it draws the outlines of the proteomic landscape. Traditional models of open reading frame ...

A Method to Study de novo Formation of Chromatin Domains

The organization and structure of chromatin domains are unique to individual cell lineages. Their misregulation might lead to a loss in cellular identity and/or disease. Despite tremendous efforts, our understanding of the formation and propagation of chromatin domains is still limited. Chromatin domains have been studied under steady-state conditions, which are not conducive to following the initial events during their establishment. Here, we present a method to inducibly reconstruct chromatin domains and follow their re-formation as a function of time. Although, first applied to the case of PRC2-mediated repressive chromatin domain formation, it could be easily adapted to other chromatin domains. The modification of and/or the combination of this method with genomics and imaging technologies will provide invaluable tools to study the establishment of chromatin domains in great detail. We believe that this method will revolutionize our understanding of how chromatin domains form and interact with each other.

A Method to Study de novo Formation of Chromatin Domains

This method is designed to follow formation of PRC2-mediated chromatin domains in cell lines, and the method can be adapted to many other systems.

クロマチンドメインのノボ形成を研究する方法

The organization and structure of chromatin domains are unique to individual cell lineages. Their misregulation might lead to a loss in cellular identity ...

3D Multicolor DNA FISH Tool to Study Nuclear Architecture in Human Primary Cells

A major question in cell biology is genomic organization within the nuclear space and how chromatin architecture can influence processes such as gene expression, cell identity and differentiation. Many approaches developed to study the 3D architecture of the genome can be divided into two complementary categories: chromosome conformation capture based technologies (C-technologies) and imaging. While the former is based on capturing the chromosome conformation and proximal DNA interactions in a population of fixed cells, the latter, based on DNA fluorescence in situ hybridization (FISH) on 3D-preserved nuclei, allows contemporary visualization of multiple loci at a single cell level (multicolor), examining their interactions and distribution within the nucleus (3D multicolor DNA FISH). The technique of 3D multicolor DNA FISH has a limitation of visualizing only a few predetermined loci, not permitting a comprehensive analysis of the nuclear architecture. However, given the robustness of its results, 3D multicolor DNA FISH in combination with 3D-microscopy and image reconstruction is a possible method to validate C-technology based results and to unambiguously study the position and organization of specific loci at a single cell level. Here, we propose a step by step method of 3D multicolor DNA FISH suitable for a wide range of human primary cells and discuss all the practical actions, crucial steps, notions of 3D imaging and data analysis needed to obtained a successful and informative 3D multicolor DNA FISH within different biological contexts.

3D multicolor DNA FISH represents a tool to visualize multiple genomic loci within 3D preserved nuclei, unambiguously defining their reciprocal interactions and localization within the nuclear space at a single cell level. Here, a step by step protocol is described for a wide spectrum of human primary cells.

ヒト初等細胞の核アーキテクチャを研究する3D多色DNA FISHツール

A major question in cell biology is genomic organization within the nuclear space and how chromatin architecture can influence processes such as gene ...

An Approach to Study Shape-Dependent Transcriptomics at a Single Cell Level

Different types of cardiac hypertrophy have been associated with an increased volume of cardiac myocytes (CMs), along with changes in CM morphology. While the effects of cell volume on gene expression are well known, the effects of cell shape are not well understood. This paper describes a method that has been designed to systematically analyze the effects of CM morphology on gene expression. It details the development of a novel single-cell trapping strategy that is then followed by single-cell mRNA sequencing. A micropatterned chip has also been designed, which contains 3000 rectangular-shaped fibronectin micropatterns. This makes it possible to grow CMs in distinct length:width aspect ratios (AR), corresponding to different types of heart failure (HF). The paper also describes a protocol that has been designed to pick up single cells from their pattern, using a semi-automated micro-pipetting cell picker, and individually inject them into a separate lysis buffer. This has made it possible to profile the transcriptomes of single CMs with defined geometrical morphotypes and characterize them according to a range of normal or pathological conditions: hypertrophic cardiomyopathy (HCM) or afterload/concentric versus dilated cardiomyopathy (DCM) or preload/eccentric. In summary, this paper presents methods for growing CMs with different shapes, which represent different pathologies, and sorting these adherent CMs based on their morphology at a single-cell level. The proposed platform provides a novel approach to high throughput and drug screening for different types of HF.

This paper presents methods for growing cardiac myocytes with different shapes, which represent different pathologies, and sorting these adherent cardiac myocytes based on their morphology at a single cell level. The proposed platform provides a novel approach to high throughput and drug screening for different types of heart failure.

単一細胞レベルでの形状依存性転写論の研究方法

Different types of cardiac hypertrophy have been associated with an increased volume of cardiac myocytes (CMs), along with changes in CM morphology. While ...

Sample Preparation to Bioinformatics Analysis of DNA Methylation: Association Strategy for Obesity and Related Trait Studies

Obesity is directly connected to lifestyle and has been associated with DNA methylation changes that may cause alterations in the adipogenesis and lipid storage processes contributing to the development of the disease. We demonstrate a complete protocol from selection to epigenetic data analysis of patients with and without obesity. All steps from the protocol were tested and validated in a pilot study. 32 women participated in the study, in which 15 individuals were classified with obesity according to Body Mass Index (BMI) (45.1 &#177; 5.4 kg/m2); and 17 individuals were classified without obesity according to BMI (22.6 &#177; 1.8 kg/m2). In the group with obesity, 564 CpG sites related to fat mass were identified by linear regression analysis. The CpG sites were in the promoter regions. The differential analysis found 470 CpGs hypomethylated and 94 hypermethylated sites in individuals with obesity. The most hypomethylated enriched pathwayswere in the RUNX, WNT signaling, and response to hypoxia. The hypermethylated pathways were related to insulin secretion, glucagon signaling, and Ca2+. We conclude that the protocol effectively identified DNA methylation patterns and trait-related DNA methylation. These patterns could be associated with altered gene expression, affecting adipogenesis and lipid storage. Our results confirmed that an obesogenic lifestyle could promote epigenetic changes in human DNA.

The present study describes the workflow to manage DNA methylation data obtained by microarray technologies. The protocol demonstrates steps from sample preparation to data analysis. All procedures are described in detail, and the video shows the significant steps.

DNAメチル化のバイオインフォマティクス分析へのサンプル調製:肥満および関連形質研究の関連戦略

Obesity is directly connected to lifestyle and has been associated with DNA methylation changes that may cause alterations in the adipogenesis and lipid ...