Name: 代谢通路信息 GWAS 分析的路径关联研究工具 (视频) | JoVE Journal
Uploaded: 2020-07-01T14:00:00.000+00:00
Duration: 5 min 1 s
Description: 了解如何使用PAST工具整合代谢通路信息与GWAS结果，识别传统分析遗漏的性状相关基因。本方法通过R代码和Shiny界面实现数据加载、连锁不平衡分析及显著通路可视化，为遗传改良提供新靶点。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

1. 设置
<ol>
	<li>如果尚未安装 R，请先安装 R。 
	注意：PAST 使用 R 编写，因此要求用户安装 R。在撰写本文时，直接从 Bioconductor 安装 PAST 需要 R4.0。对于 R3.6，可从 Bioconductor 安装较早版本的 PAST；而对于 R3.5 用户，可从 Github 安装 PAST。R 的安装说明可从以下链接下载：https://www.r-project.org/。</li>
	<li>安装最新版本的 RStudio Desktop 或更新现有 RStudio（可选）。 
	注意：RStudio 是一个有助于使用 R 语言的开发环境，建议安装，特别是对于选择通过命令行而非 Shiny 图形用户界面运行 PAST 的用户。RStudio 及其安装说明可在以下链接获取：<a href="https://rstudio.com/products/rstudio/">https://rstudio.com/products/rstudio/</a>。</li>
	<li>按照 Bioconductor 上的说明，从 Bioconductor 安装 PAST11。 
	注意：通过 Bioconductor 安装将自动处理 PAST 所需依赖包的安装。此外，PAST 也可从 Github12 安装，但从 Github 安装不会自动安装依赖包。</li>
	<li>安装 PAST Shiny（可选）。从 Github 代码仓库的 Releases 页面下载文件 "app.R"：https://github.com/IGBB/PAST/releases/，并记住下载文件的存储位置。 
	注意：用户可通过在 R 中直接调用 PAST 的方法来使用该工具，但对于不熟悉 R 的用户，可运行 PAST Shiny 应用程序，其提供引导式用户界面。PAST Shiny 是一个 R 脚本，位于 PAST Github 仓库的 shiny_app 分支中。PAST Shiny 在首次运行时会尝试自动安装其依赖包。</li>
	<li>通过以下三种方式之一启动应用程序，开始分析。
	<ol>
		<li>PAST Shiny 与 RStudio 配合使用
		<ol>
			<li>在 RStudio 中，于 app.R 所在的文件夹中创建一个新项目。点击文件 | 新建项目，并选择该文件夹。</li>
			<li>创建新项目后，打开之前下载的 app.R 文件。RStudio 会识别出 app.R 是一个 Shiny 应用，并在显示的源代码上方工具栏中生成一个 运行应用 按钮。点击 运行应用。RStudio 随即会启动一个窗口，显示 PAST Shiny 应用程序。</li>
		</ol>
		</li>
		<li>PAST Shiny 与 R 控制台配合使用
		<ol>
			<li>启动 R，并运行以下代码以启动 PAST Shiny 应用程序：shiny::runApp('path/to/folder/with/shiny/app.R')。请将引号内的文本替换为 app.R 文件的下载路径，并保留引号。</li>
		</ol>
		</li>
		<li>不使用 R Shiny 的 PAST 使用方式
		<ol>
			<li>在 R 控制台中运行 library(PAST) 以加载 PAST。</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
2. 自定义 Shiny 分析（可选）
<ol>
	<li>将分析标题从“New Analysis”更改为更能反映所执行分析类型的名称，以便于跟踪多个分析（见图1）。</li>
</ol>
<img alt="遗传数据分析的PAST软件界面；输入选项、分析设置、基因通路、实验设置。" class="xfigimg" src="/files/ftp_upload/61268/61268fig1.jpg" data-altupdatedat="1762425657000" data-alt="图1"> 
图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>使用 PAST Shiny 加载 GWAS 数据。
	<ol>
		<li>使用关联文件和效应文件选择框分别选择一个关联文件和一个效应文件。在文件选择框下方的关联列名称和效应列名称输入框中，修改列名称以匹配数据中的实际列名。</li>
	</ol>
	</li>
</ol>
<img alt="用于通路绘制和基因分配实验的遗传分析软件界面。" class="xfigimg" src="/files/ftp_upload/61268/61268fig2.jpg" data-altupdatedat="1762425657000" data-alt="图 2"> 
图 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", "path/to/effects_file.tsv", association_columns = c("Trait", "Marker", "Locus", "Site", "p", "marker_R2"), effects_columns = c("Trait", "Marker", "Locus", "Site", "Effect")</li>
	</ol>
	</li>
	<li>注意：请将路径更改为 GWAS 文件的实际存储位置。association_columns 和 effects_columns 所提供的值为默认值。如果列名与默认值不一致，则需明确指定列名；否则可省略这些参数。</li>
</ol>
4. 加载连锁不平衡（LD）数据
注意：请验证连锁不平衡（LD）数据是否为制表符分隔，且包含以下类型的数据：位点（Locus）、位置1（Position1）、位点1（Site1）、位置2（Position2）、位点2（Site2）、位置1与位置2之间的碱基对距离，以及 R2 值。
<ol>
	<li>使用 PAST Shiny 加载 LD 数据。
	<ol>
		<li>选择包含 LD 数据的文件。如有必要，在文件选择框下方的 LD 列名称 输入框中修改列名称，使其与 LD 数据中的列名称一致。</li>
	</ol>
	</li>
</ol>
<img alt="基因数据分析界面、示意图、通路发现设置、基因分配、连锁研究。" class="xfigimg" src="/files/ftp_upload/61268/61268fig3.jpg" data-altupdatedat="1762425657000" data-alt="图3"> 
图3。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig3large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<ol start="2">
	<li>在 R 控制台中使用 PAST 加载 LD 数据。
	<ol>
		<li>修改并运行以下代码以加载 LD 数据： 
		LD = load_LD("path/to/LD.tsv", LD_columns = c("Locus1", "Position1", "Site1", "Position2", "Site2", "Dist_bp", "R.2") 
		注意：请将路径更改为 LD 文件的实际位置。LD_columns 提供的值为默认值。如果列名与这些默认值不匹配，请指定正确的列名；否则可省略此参数。</li>
	</ol>
	</li>
</ol>
5. 将SNP分配给基因
注意：下载或以其他方式获取 GFF 格式的注释文件。这些注释通常可在特定生物体的在线数据库中找到。需谨慎对待质量较低的注释，因为注释数据的质量将直接影响通路分析的结果质量。请确认这些注释文件的第一列（染色体）与关联分析、效应值及 LD 数据中的位点/染色体命名格式一致。例如，如果 GWAS 和 LD 数据文件中将第一条染色体记为“1”，则注释文件不应将其记为“chr1”。
<ol>
	<li>使用 PAST Shiny 将 SNP 分配至基因。 
	注意：有关确定合适 R2 截断值的更多信息，请参见 Tang 等人6 的文章中名为“通路分析的 SNP 至基因算法”部分。
	<ol>
		<li>选择包含 GFF 注释的文件。根据所研究物种的特点，考虑最合适的窗口大小和 R2 截断值，若默认值不适用于上传的数据，则进行相应修改。 
		注意：PAST 中的默认值主要适用于玉米。PAST Shiny 分析开始时设置的核数（步骤 2.2）将在此步骤中使用。</li>
	</ol>
	</li>
</ol>
<img alt="GWAS、通路发现、基因注释的生物信息学分析界面；输入与结果。" class="xfigimg" src="/files/ftp_upload/61268/61268fig4.jpg" data-altupdatedat="1762425657000" data-alt="图4"> 
图4. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig4large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<ol start="2">
	<li>在 R 控制台中使用 PAST 将 SNP 分配给基因。
	<ol>
		<li>修改并运行以下代码，以将 SNP 分配给基因： 
		genes = assign_SNPs_to_genes(gwas_data, LD, "path/to/annotations.gff", c("gene"), 1000, 0.8, 2) 
		注意：在此示例代码中，提供了若干默认参数：1000 表示在 SNP 周围搜索基因的窗口大小；0.8 为 R2 的阈值；2 为用于并行处理的计算核心数。注释文件的路径也应更改为注释文件的实际存储位置。</li>
	</ol>
	</li>
</ol>
6. 发现显著性通路
注意：请验证通路文件是否包含以下以制表符分隔的数据，每条通路中的每个基因对应一行：通路 ID——例如“PWY-6475-1”的标识符；通路描述——对通路功能的详细说明，例如“反式番茄红素生物合成”；基因——通路中的一个基因，其名称应与注释中提供的名称一致。通路信息通常可在特定生物的在线数据库中找到，例如 MaizeGDB。第二个用户自定义选项是模式。“增加”指表型测量值越高越有利的情况，例如产量；而“减少”指测量值越低越有利的性状，例如虫害评分。通路的显著性通过先前描述的方法进行检验4,6,14。
<ol>
	<li>使用 PAST Shiny 发现显著通路。
	<ol>
		<li>选择包含通路数据的文件，并确保在分析选项中已选择相应的模式。如有必要，可更改保留在分析中的通路所必须包含的基因数量，以及用于构建零分布以检验效应显著性的置换次数。</li>
	</ol>
	</li>
</ol>
<img alt="基因组数据分析软件界面，显示全基因组关联分析、连锁不平衡和基因分配的选项。" class="xfigimg" src="/files/ftp_upload/61268/61268fig5.jpg" data-altupdatedat="1762425657000" data-alt="图 5"> 
图 5．<a href="https://www.jove.com/files/ftp_upload/61268/61268fig5large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
注意：此步骤使用在 PAST Shiny 分析开始时设置的核心数量和模式（步骤 2.2）。默认的基因数量目前设定为 5 个基因，因此已知基因数少于该值的通路将被剔除。用户可将此值降低至 4 或 3 以包含较短的通路，但这样做会增加假阳性结果的风险。提高该值可增强分析的效力，但会导致更多通路被排除在分析之外。改变所用置换次数的数值会相应地提高或降低检验的效力。
<ol start="2">
	<li>在 R 控制台中使用 PAST 发现显著通路。
	<ol>
		<li>修改并运行以下代码以发现显著通路： 
		rugplots_data &lt;- find_pathway_significance(基因, "path/to/pathways.tsv", 5, "increasing", 1000, 2) 
		​注意：在此示例代码中，提供了若干建议的默认参数。其中，5 表示通路中必须包含的最少基因数，以便在分析中保留该通路；"increasing" 表示所测性状呈递增趋势（建议用户无论性状如何，均分别运行递增和递减方向；但两种方向的数据解释将有所不同）；1000 表示为确定零分布而进行抽样的次数；2 表示用于并行处理的计算核心数。请将路径修改为通路文件的实际存储位置。</li>
	</ol>
	</li>
</ol>
7. 查看 Rug图
<ol>
	<li>使用 PAST Shiny 查看 Rugplots。
	<ol>
		<li>上传并设置所有输入后，点击 开始分析。此时将显示一个进度条，指示分析流程中最后完成的步骤。分析完成后，PAST Shiny 将切换至“结果”标签页。左侧列（标注为“通路”）将显示结果表格，右侧列（标注为“图谱”）将显示 Rugplots。</li>
		<li>使用滑块调节筛选参数。当筛选水平达到满意时，点击左下角的 下载结果 按钮，将所有图像和表格分别下载为一个以分析标题命名的 ZIP 文件。该 ZIP 文件包含筛选后的表格、未筛选的表格，以及筛选后表格中每个通路对应的单张图像。</li>
	</ol>
	</li>
</ol>
<img alt="SNP 分配至基因的进度条显示 75%，遗传数据分析界面，生物信息学。" class="xfigimg" src="/files/ftp_upload/61268/61268fig6.jpg" data-altupdatedat="1762425657000" data-alt="图 6"> 
图 6。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig6large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<img alt="通路显著性筛选界面与基因合成图谱；数据分析结果。" class="xfigimg" src="/files/ftp_upload/61268/61268fig7.jpg" data-altupdatedat="1762425657000" data-alt="图7"> 
图7。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig7large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<ol start="2">
	<li>在 R 控制台中使用 PAST 查看 Rugplots
	<ol>
		<li>修改并运行以下代码以保存结果： 
		plot_pathways(rugplots_data, "pvalue", 0.02, "increasing", "output_folder") 
		注意：在此示例代码中提供了若干建议的默认参数。pvalue 提供可用于在用户选择显著性阈值后过滤不显著通路的数据；0.02 是用于过滤的默认值；increasing 表示所测量性状的增加趋势（建议用户无论性状如何，均分别运行 increasing 和 decreasing 两种情况；然而，两者的数据解释将有所不同）；output_folder 是图像和表格将被写入的文件夹（该文件夹必须在运行函数前已存在）。过滤后的结果表格、未过滤的结果以及每个过滤后通路的单独图像将被写入此文件夹。</li>
	</ol>
	</li>
</ol>

<ol><li>Rafalski, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Association+genetics+in+crop+improvement%5BTitle%5D)+AND+%22Current+Opinion+in+Plant+Biology%22%5BJournal%5D)">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/pubmed?term=((Association+Mapping+for+Enhancing+Maize+%28Zea+mays+L.%29+Genetic+Improvement%5BTitle%5D)+AND+%22Crop+Science%22%5BJournal%5D)">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/pubmed?term=((Genome-wide+Association+Studies+in+Maize%3A+Praise+and+Stargaze%5BTitle%5D)+AND+%22Molecular+Plant%22%5BJournal%5D)">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/pubmed?term=((Pathway-Based+Approaches+for+Analysis+of+Genomewide+Association+Studies%5BTitle%5D)+AND+%22The+American+Journal+of+Human+Genetics%22%5BJournal%5D)">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/pubmed?term=((SNP-based+pathway+enrichment+analysis+for+genome-wide+association+studies%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">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/pubmed?term=((Using+genome-wide+associations+to+identify+metabolic+pathways+involved+in+maize+aflatoxin+accumulation+resistance%5BTitle%5D)+AND+%22BMC+Genomics%22%5BJournal%5D)">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/pubmed?term=((Genome-Wide+Association+Mapping+of+Aspergillus+flavus+and+Aflatoxin+Accumulation+Resistance+in+Maize%5BTitle%5D)+AND+%22Crop+Science%22%5BJournal%5D)">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/pubmed?term=((Genome-Wide+Association+and+Metabolic+Pathway+Analysis+of+Corn+Earworm+Resistance+in+Maize%5BTitle%5D)+AND+%22The+Plant+Genome%22%5BJournal%5D)">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/pubmed?term=((Leveraging+GWAS+data+to+identify+metabolic+pathways+and+networks+involved+in+maize+lipid+biosynthesis%5BTitle%5D)+AND+%22The+Plant+Journal%22%5BJournal%5D)">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/pubmed?term=((PAST%3A+The+Pathway+Association+Studies+Tool+to+Infer+Biological+Meaning+from+GWAS+Datasets%5BTitle%5D)+AND+%22Plants%22%5BJournal%5D)">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://scholar.google.com/scholar?hl=en&safe=off&q=author%3aAdam%2C+T+%22PAST%3A+Pathway+Association+Study+Tool+%28PAST%29.+Bioconductor+version%3A+Release+%283.10%29%22">PAST: Pathway Association Study Tool (PAST). Bioconductor version: Release (3.10)</a>. , (2020).</li>
<li>Thrash, A., DeOrnellis, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aThrash%2C+A+%22IGBB%2FPAST%22">IGBB/PAST</a>. , at https://github.com/IGBB/PAST (2019).</li>
<li>Bradbury, P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TASSEL%3A+software+for+association+mapping+of+complex+traits+in+diverse+samples%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">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/pubmed?term=((Gene+set+enrichment+analysis%3A+A+knowledge-based+approach+for+interpreting+genome-wide+expression+profiles%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+U.S.A%22%5BJournal%5D)">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 时，易用性是首要目标，这体现在选择开发 Shiny 应用程序以及采用 R 和 Bioconductor 来发布该应用程序。用户无需学习如何编译程序即可使用 PAST。
与大多数分析软件一样，PAST 的分析结果质量取决于输入数据的质量；如果输入数据存在错误或格式不正确，PAST 将无法运行或产生无意义的结果。确保全基因组关联研究（GWAS）数据、连锁不平衡（LD）数据、注释信息和通路文件的格式正确，是获得准确输出结果的关键。PAST 仅能分析双等位基因标记，并且每组输入数据只能分析一个性状。此外，由低质量基因分型或错误、不精确的表型鉴定所产生的 GWAS 数据，也不太可能产生清晰或可重复的结果。PAST 可辅助对 GWAS 结果进行生物学解释，但如果环境变异、实验误差或群体结构未得到适当校正，则难以澄清混乱的数据集。
用户可以选择更改分析的某些参数，既可以在 Shiny 应用程序中进行，也可以通过将这些参数传递给 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。因此，尽管可能保留许多与性状相关的SNP，但每个SNP的效应可能非常微弱。过多缺失的SNP调用会导致对性状生物学意义及遗传结构的解释不完整且令人困惑。代谢通路分析可通过关注按生物学功能分组的基因的综合效应，帮助解决其中部分问题4,5,6。
已有若干研究采用了本文所述方法的先前版本。先前版本曾被用于研究黄曲霉毒素积累7、玉米穗虫抗性8以及油脂生物合成9。尽管这些分析取得了成功，但分析过程复杂、耗时且繁琐，原因是分析工具由 R、Perl 和 Bash 多种语言编写，且流程未实现自动化。由于修改该方法以适应不同分析需要专业知识，目前开发出了一种新方法，可供其他研究人员共享使用。
通路关联分析工具（Pathway Association Study Tool, PAST）10 旨在克服先前方法的不足，其设计减少了对编程语言知识的需求，并可在更短时间内完成分析。尽管该方法已在玉米中进行了测试，但PAST并不依赖于特定物种的假设。PAST可通过R控制台运行，也可作为Shiny应用程序使用，且预计不久后将在MaizeGDB上推出在线版本。

<table><tbody><tr><td>计算机</td><td>NA</td><td>NA</td><td>任何具有 8GB 内存的计算机均应足够</td></tr><tr><td>R</td><td>R Project</td><td>NA</td><td>需要 R 4.0 或更高版本才能从 Bioconductor 3.11 安装</td></tr></tbody></table>

一种用于代谢通路信息全基因组关联分析的通路关联分析工具

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

最近，一种基于代谢通路分析解读全基因组关联研究（GWAS）数据的已有方法得到了新的实现并已发布。该通路关联分析工具（PAST）的开发旨在解决原有工具在用户友好性方面的不足以及分析速度缓慢的问题。这一新型用户友好型工具已发布于Bioconductor和Github平台。在测试中，PAST将原本需要二十四小时或更长时间的分析任务缩短至一小时以内完成。本文中，我们介绍了使用Shiny应用程序或R控制台运行PAST的操作方案。

如果运行 PAST 软件工具后未生成结果，请检查所有输入文件是否格式正确。使用 PAST 软件包中提供的示例数据（基于玉米籽粒颜色的全基因组关联分析）成功运行的结果如图8所示。该表格及生成的图像可通过“下载结果”按钮进行下载。下载后的图像示例如图210所示。不正确的设置可能导致结果在生物学上无意义，但判断结果是否正确应由研究人员自行完成，需仔细核对所选设置的有效性，并综合考虑与目标性状相关的所有已知证据。
图910展示了利用包含288个自交系的玉米群体进行全基因组关联分析（GWAS）所得结果的通路分析生成的rugplot。该示例较为简单，表型仅为“白色”或“黄色”，因为已知负责产生亮黄色类胡萝卜素色素的代谢通路，并且该通路应能解释大部分表型变异。因此，我们预期会观察到反式番茄红素生物合成通路（即类胡萝卜素合成通路）与籽粒颜色显著相关，而结果确实如此。图中顶部列出了通路的编号和名称。图的横轴表示分析中包含的所有基因，按其对性状影响效应从大到小从左至右排列。然而，仅标出了属于反式番茄红素生物合成通路的基因（在图顶部以短横线表示，其位置对应于该基因在所有分析基因中的效应排序）。该通路共包含7个基因。纵轴表示累积富集分数（ES）。每个基因的ES值按照效应大小顺序累加至总分中，并根据分析的基因总数进行调整。因此，随着沿横轴向右移动，分数不断变化，当纳入效应较大的基因时，总分趋于上升；但在某一节点后，新增基因带来的效应增量小于因增加基因数而进行的调整量，导致总分开始下降。累积ES曲线的峰值由一条虚线垂直标出；该值即为整个通路的ES，程序据此判断该通路是否被选中并以rugplot形式展示。&#160;
<img alt="图8" class="xfigimg" src="/files/ftp_upload/61268/61268fig8.jpg" /> 
图8：PAST Shiny 完成运行结果。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig8large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>
<img alt="图9" class="xfigimg" src="/files/ftp_upload/61268/61268fig9.jpg" /> 
图9：PAST 完整运行（或从 Shiny 下载）生成的通路图。 本图引自 Thrash 等10。 <a href="https://www.jove.com/files/ftp_upload/61268/61268fig9large.jpg" target="_blank">请点击此处查看此图的放大版本。</a>

通过使用 Shiny 应用程序或 R 控制台运行通路关联分析工具（PAST），研究人员可以通过研究相关的代谢通路，更深入地理解全基因组关联研究（GWAS）结果的生物学意义。

最近，一种基于代谢通路分析全基因组关联研究（GWAS）数据的已有方法得到了新的应用……

全基因组关联分析是发现与性状相关基因的强大工具，但会遗漏其中许多基因。PAST 能够找到更多的这类基因，从而更完整地揭示整体图景。PAST 以全基因组关联分析的输出结果为基础，加入细胞用于合成代谢物及形成生物体表型的通路信息。仅通过全基因组关联研究（GWAS）通常会遗漏这些通路中的基因。在确认GWAS数据为制表符分隔格式后，使用“关联文件”选择框选取关联文件，并使用“效应文件”选择框选取效应文件。要在R控制台中使用PAST加载GWAS数据，请按照提示修改并运行相应代码。使用 PAST Shiny 加载连锁不平衡数据时，首先确认连锁不平衡数据为制表符分隔，且包含适当类型的数据，然后选择包含连锁不平衡数据的文件。若要在 R 控制台中使用 PAST 加载连锁不平衡数据，请按说明修改并运行代码。使用 PAST Shiny 将 SNP 分配至基因时，需定位通用特征格式（general feature format）的注释文件，选择包含通用特征格式注释的文件，并根据所研究物种的适宜窗口大小和 R 平方截断值，按需调整设置。在 R 控制台中使用 PAST 将 SNP 分配给基因时，按指示修改并运行代码。在确认通路文件包含以制表符分隔格式的适当数据后，选择包含通路数据的文件，并确认分析选项中已选中相应的模式。如有必要，更改通路中必须包含的基因数量，以确定保留用于分析的通路，并更改用于创建零分布的置换次数，以检验效应的显著性。要在 R 控制台中使用 PAST 发现重要通路，请按照指示修改并运行代码。要在 PAST Shiny 中查看地毯图，在上传并设置所有输入后，点击“开始分析”。将出现一个进度条，指示分析的最后完成步骤。分析完成后，PAST Shiny 将切换到结果选项卡，并显示结果表格和地毯图。使用滑块来调整过滤参数。当过滤水平达到满意时，点击下载结果按钮，即可将所有图像和表格单独下载为一个以分析标题命名的压缩文件（zip 格式）。要在 R 控制台中使用 PAST 查看 rug 图，请修改并运行代码以保存结果。此处展示的是使用 PAST 软件包中提供的示例数据（基于籽粒颜色的全基因组关联研究）成功运行的结果。同时展示了点击“下载结果”按钮后可下载的图像示例。该峰线图基于利用288个近交系组成的迷宫面板进行的全基因组关联分析（GWAS）结果生成，这些近交系均具有籽粒颜色表型数据。正如预期，产生类胡萝卜素的反式番茄红素生物合成通路与籽粒颜色显著相关。在基因效应排序中，该通路中的基因以及反式番茄红素生物合成通路本身也相对于分析中的所有其他基因进行了标注。整个通路的富集评分曲线最高点由虚线标示。程序利用该点来判断该通路是否被选中并以地毯图形式展示。PAST 可应用于所有数据，但仅对高质量数据产生良好结果。序列必须包含较少的错误，且表型数据应在严格的统计设计下进行收集。通过PAST结果，可以鉴定出控制关键位点和重要通路的基因，并最终通过选择育种、基因编辑和药物治疗等手段对这些基因进行操控，以改良目标性状。

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

Research

科研

JoVE Journal

Genetics

遗传学

代谢通路信息 GWAS 分析的路径关联研究工具

结合基因组学、共表达基因分析以及通过代谢途径鉴定靶向化合物，可实现基因功能注释。...

通过整合基因组学、代谢组学与生物信息学注释植物基因功能

一种用于生成高质量机器可读生物通路的网络工具

我们提供了一种标准化的转录组数据基因集富集分析方案，用于确定适用于转化研究的理想小鼠模型。<br...

使用基因集富集分析确定转化研究中合适动物模型的方案

在本实验方案中，我们介绍了一种优化的工作流程，该流程结合了高效且快速的大量样本制备方法。此外，我们还提供了逐步操作指南，以减少分析变异，实现对代谢组全基因组关联研究的高通量评估。...

大规模多组学全基因组关联研究（Mo-GWAS）：样本制备与标准化操作指南

我们介绍了 CorrelationCalculator 和 Filigree 两种用于代谢组学数据驱动网络构建与分析的工具。CorrelationCalculator 支持基于表达数据构建单一的代谢物相互作用网络，而 Filigree...

CorrelationCalculator 和 Filigree：面向代谢组学数据的基于数据驱动的网络分析工具

DeepOmicsAE：利用蛋白质组学、代谢组学和临床数据的深度学习分析揭示阿尔茨海默病中的信号模块

胰岛β细胞死亡发生在1型糖尿病发病之前，检测这一过程可能有助于早期进行治疗干预。本文详细描述了如何测定人血清中胰岛素基因（INS）不同甲基化DNA序列，作为β细胞死亡的生物标志物。...

人血清样本中差异甲基化INS DNA 物种的检测作为胰岛β细胞死亡的生物标志物

在微观生物的基因组测序过程中，污染问题仍然十分严重。本文介绍一种方法，仅需50...

单个缓步动物标本的超低起始量基因组测序文库制备

本文介绍了一种可准确量化线粒体DNA（mtDNA）甲基化水平的实验方案。该方案描述了使用BamHI酶对DNA进行消化，并结合生物信息学分析流程，以避免因mtDNA二级结构导致的mtDNA甲基化水平高估问题。...

线粒体DNA甲基化的精确定量检测方法

鉴定基因与调控元件之间的物理相互作用具有挑战性，但染色体构象捕获技术为此提供了便利。该4C-seq方案的改进通过减少PCR模板的过度扩增，减轻了PCR偏好性，同时通过增加限制性内切酶消化步骤，提高了测序读段的可比对性。...

利用下一代测序技术对染色体构象捕获环化文库（4C-seq）进行基因调控序列的高通量鉴定

了解如何使用PAST工具整合代谢通路信息与GWAS结果，识别传统分析遗漏的性状相关基因。本方法通过R代码和Shiny界面实现数据加载、连锁不平衡分析及显著通路可视化，为遗传改良提供新靶点。

一种用于代谢通路信息全基因组关联分析的通路关联分析工具

探索更多视频