Name: Un outil d’étude d’association de voies pour les analyses GWAS de l’information sur les voies métaboliques (Vidéo) | JoVE Journal
Uploaded: 2020-07-01T14:00:00.000+00:00
Duration: 5 min 1 s
Description: Découvrez l'analyse d'association de voies métaboliques avec l'outil PAST pour compléter les résultats GWAS. Cette méthode bioinformatique identifie les gènes et voies significatifs liés aux traits d'intérêt.

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. Configuration
<ol><li>Installez R, s’il n’est pas déjà installé. REMARQUE: PAST est écrit en R et, par conséquent, nécessite que ses utilisateurs aient installé R. Au moment d’écrire ces lignes, l’installation de PAST directement à partir de Bioconductor nécessite R4.0.  Les anciennes versions de PAST peuvent être installées à partir de Bioconductor pour R3.6, et PAST peut être installé à partir de Github pour les utilisateurs de R3.5. Les instructions d’installation de R peuvent être téléchargées à partir du lien suivant : https://www.r-project.org/.</li>
	<li>Installez la dernière version de RStudio Desktop ou mettez à jour RStudio (facultatif). REMARQUE: RStudio est un environnement utile pour travailler avec le langage R. Son installation est recommandée, en particulier pour ceux qui choisissent d’exécuter PAST dans la ligne de commande plutôt que via l’application Shiny GUI. RStudio et ses instructions d’installation peuvent être trouvés sur le lien suivant: <a href="https://rstudio.com/products/rstudio/">https://rstudio.com/products/rstudio/</a>.</li>
	<li>Installez PAST à partir de Bioconductor11 en suivant les instructions sur Bioconductor. REMARQUE: L’installation via Bioconductor doit gérer l’installation des dépendances de PAST. En outre, PAST peut être installé à partir de Github12, mais l’installation à partir de Github n’installera pas automatiquement les dépendances.</li>
	<li>Installez PAST Shiny (facultatif). Téléchargez le fichier « app. R » à partir de la page Versions du référentiel Github : https://github.com/IGBB/PAST/releases/ et rappelez-vous où se trouve le fichier téléchargé. REMARQUE : PAST peut être utilisé en appelant ses méthodes directement avec R, mais les utilisateurs qui sont moins familiers avec R peuvent exécuter l’application PAST Shiny, qui fournit une interface utilisateur guidée. PAST Shiny est un script R disponible dans la branche shiny_app du référentiel GITHUB PAST. PAST Shiny tentera d’installer ses dépendances lors de la première exécution.</li>
	<li>Commencez l’analyse en démarrant l’application de l’une des trois manières décrites ci-dessous.<ol>
<li>PAST Shiny avec RStudio<ol>
<li>À l’aide de RStudio, créez un nouveau projet dans le dossier où se trouve l’application. R est situé. Cliquez sur | fichier Nouveau projet et sélectionnez ce dossier.</li>
			<li>Une fois qu’un nouveau projet a été créé, ouvrez l’application. R téléchargé précédemment. RStudio reconnaît cette application. R est une application Shiny et crée un bouton Exécuter l’application sur la barre au-dessus du code source affiché. Cliquez sur Exécuter l’application. RStudio lancera alors une fenêtre qui affiche l’application PAST Shiny.</li>
		</ol>
</li>
		<li>PAST Shiny avec console R<ol>
<li>Lancez R et exécutez le code suivant pour démarrer l’application PAST Shiny : shiny::runApp('path/to/folder/with/shiny/app. R'. Remplacez le texte entre guillemets par le dossier dans lequel l’application. R a été téléchargé et conservez les citations.</li>
		</ol>
</li>
		<li>PASSÉ sans R Shiny<ol>
<li>Exécutez library(PAST) dans une console R pour charger PAST.</li>
		</ol>
</li>
	</ol>
</li>
</ol>2. Personnaliser l’analyse Shiny (facultatif)
<ol><li>Remplacez le titre de l’analyse de « Nouvelle analyse » par quelque chose qui reflète mieux le type d’analyse en cours d’exécution, ce qui permet de suivre plusieurs analyses (voir la figure 1).</li>
</ol><img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/61268/61268fig1.jpg"> Figure 1. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig1large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<ol start="2"><li>Modifiez le nombre de cœurs et le mode. Définissez le nombre de cœurs sur n’importe quel nombre compris entre 1 et le nombre total sur la machine, mais sachez que consacrer plus de ressources à PAST peut ralentir d’autres opérations sur la machine. Définissez le mode en fonction de la description de la section 6.</li>
</ol>3. Charger les données GWAS
Remarque : Vérifiez que les données GWAS sont délimitées par des tabulations. Assurez-vous que le fichier d’association contient les colonnes suivantes : trait, nom du marqueur, locus ou chromosome, position sur le chromosome, valeur p et valeur R2 pour le marqueur. Assurez-vous que le fichier d’effets contient les colonnes suivantes : trait, nom du marqueur, locus ou chromosome, position sur le chromosome et effet. L’ordre de ces colonnes n’est pas important, car l’utilisateur peut spécifier les noms des colonnes lors du chargement des données. Toutes les colonnes supplémentaires sont ignorées. TASSEL13 peut être utilisé pour produire ces fichiers.
<ol><li>Chargez les données GWAS avec PAST Shiny.<ol>
<li>Sélectionnez un fichier d’association et un fichier d’effets à l’aide des zones de sélection Fichier d’association et Fichier d’effets. Modifiez les noms des colonnes dans les zones de saisie Nom de la colonne d’association et Nom des colonnes d’effets sous les zones de sélection de fichier pour refléter les noms de colonne dans les données. </li>
	</ol>
</li>
</ol><img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/61268/61268fig2.jpg"> Figure 2. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig2large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<ol start="2"><li>Chargez les données GWAS avec PAST dans la console R.<ol>
<li>Modifiez et exécutez le code suivant : 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>REMARQUE : Modifiez les chemins d’accès à l’emplacement réel des fichiers GWAS. Les valeurs fournies pour association_columns et effects_columns sont les valeurs par défaut. Si les noms ne correspondent pas aux valeurs par défaut, spécifiez les noms des colonnes. Sinon, ceux-ci peuvent être omis.</li>
</ol>4. Données sur le déséquilibre de liaison de charge (LD)
Remarque : Vérifiez que les données de déséquilibre de liaison (LD) sont délimitées par des tabulations et contiennent les types de données suivants : Locus, Position1, Site1, Position2, Site2, Distance dans les paires de bases entre Position1 et Position2 et valeur R2.
<ol><li>Chargez les données LD avec PAST Shiny.<ol>
<li>Sélectionnez le fichier contenant les données LD. Modifiez les noms des colonnes dans les zones de saisie Noms des colonnes LD sous la zone de sélection de fichier pour qu’elles correspondent aux noms des colonnes dans les données LD si nécessaire.</li>
	</ol>
</li>
</ol><img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/61268/61268fig3.jpg"> Figure 3. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig3large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<ol start="2"><li>Chargez LD Data avec PAST dans la console R.<ol>
<li>Modifiez et exécutez le code suivant pour charger les données LD : LD = load_LD(« path/to/LD.tsv », LD_columns = c(« Locus1 », « Position1 », « Site1 », « Position2 », « Site2 », « Dist_bp », « R.2 ») Remarque : Modifiez le chemin d’accès à l’emplacement réel du fichier LD. Les valeurs fournies pour LD_columns sont les valeurs par défaut. Si les noms ne correspondent pas à ces valeurs par défaut, spécifiez les noms corrects des colonnes ; sinon, ceux-ci peuvent être omis.</li>
	</ol>
</li>
</ol>5. Assigner des SLP aux gènes
REMARQUE: Téléchargez ou recherchez des annotations au format GFF. Ces annotations peuvent souvent être trouvées dans des bases de données en ligne pour des organismes spécifiques. Soyez prudent avec les annotations de faible qualité, car la qualité des données d’annotations affectera la qualité de l’analyse des voies. Vérifiez que la première colonne de ces annotations (le chromosome) correspond au format du locus/chromosome dans les données d’association, d’effets et de LD. Par exemple, les annotations ne doivent pas appeler le premier chromosome « chr1 » si les fichiers de données GWAS et LD appellent le premier chromosome « 1 ».
<ol><li>Attribuez des SLP aux gènes avec PAST Shiny. REMARQUE: Plus d’informations sur la détermination d’un seuil R2 approprié peuvent être trouvées dans Tang et al.6, dans la section intitulée « SNP to gene algorithm for the pathway analysis ».<ol>
<li>Sélectionnez le fichier contenant les annotations GFF. Considérez quelle taille de fenêtre et quelle coupure R2 conviennent le mieux à l’espèce considérée et modifiez-la si les valeurs par défaut ne conviennent pas aux données téléchargées. REMARQUE: Les valeurs par défaut dans PAST reflètent principalement les valeurs appropriées pour le maïs. Le nombre de cœurs défini au début de l’analyse PAST Shiny (étape 2.2) est utilisé dans cette étape.</li>
	</ol>
</li>
</ol><img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/61268/61268fig4.jpg"> Figure 4. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig4large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<ol start="2"><li>Attribuez des SLP aux gènes avec PAST dans la console R.<ol>
<li>Modifiez et exécutez le code suivant pour affecter des SMP aux gènes : gènes = assign_SNPs_to_genes(gwas_data, LD, « path/to/annotations.gff », c(« gene »), 1000, 0.8, 2) REMARQUE: Dans cet exemple de code, plusieurs suggestions par défaut sont fournies: 1000 est la taille de la fenêtre autour du SNP pour rechercher des gènes; 0,8 est la valeur seuil pour R2; 2 est le nombre de cœurs utilisés pour le traitement parallèle. Le chemin d’accès aux annotations doit également être remplacé par l’emplacement réel du fichier d’annotations.</li>
	</ol>
</li>
</ol>6. Découvrez des voies importantes
REMARQUE: Vérifiez que le fichier de voies contient les données suivantes dans un format délimité par des tabulations, avec une ligne pour chaque gène dans chaque voie: ID de voie - un identificateur tel que « PWY-6475-1 »; description de la voie - une description plus longue de ce que font les voies, comme la « biosynthèse trans-lycopène »; gène - un gène dans la voie, qui doit correspondre aux noms fournis dans les annotations. Les informations sur les voies d’accès peuvent probablement être trouvées dans des bases de données en ligne pour des organismes spécifiques, tels que MaizeGDB. La deuxième option spécifiée par l’utilisateur est le mode. « Croissant » fait référence aux phénotypes qui reflètent quand une valeur croissante du trait mesuré est souhaitable, comme le rendement, tandis que « décroissant » fait référence à un trait où une diminution des valeurs mesurées est bénéfique, comme les cotes de dommages causés par les insectes. L’importance des voies est testée à l’aide des méthodesdécritesprécédemment 4,6,14.
<ol><li>Découvrez des voies importantes avec PAST Shiny.<ol>
<li>Sélectionnez le fichier contenant les données des chemins d’accès et assurez-vous que le mode est sélectionné dans les options d’analyse. Si nécessaire, modifiez le nombre de gènes qui doivent se présenter dans une voie pour le conserver pour l’analyse et le nombre de permutations utilisées pour créer la distribution nulle afin de tester la signification de l’effet.</li>
	</ol>
</li>
</ol><img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/61268/61268fig5.jpg"> Graphique 5. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig5large.jpg" target="_blank">Veuillez cliquer ici pour voir une version agrandie de cette figure.</a>
REMARQUE : Le nombre de cœurs et le mode défini au début de l’analyse PAST Shiny (étape 2.2) sont utilisés dans cette étape. Le nombre par défaut de gènes est actuellement fixé à 5 gènes, de sorte que les voies avec moins de gènes connus seront supprimées. L’utilisateur peut abaisser cette valeur à 4 ou 3, pour inclure des voies plus courtes, mais cela risque de fausses résultats positifs. L’augmentation de cette valeur peut augmenter la puissance de l’analyse, mais supprimera davantage de voies de l’analyse. Changer le nombre de permutations utilisées augmente et diminue la puissance du test.
<ol start="2"><li>Découvrez des voies importantes avec PAST dans la console R.<ol>
<li>Modifiez et exécutez le code suivant pour découvrir des chemins d’accès significatifs : rugplots_data &lt;- find_pathway_significance(gènes, « path/to/pathways.tsv », 5, « increasing », 1000, 2) Remarque : Dans cet exemple de code, plusieurs valeurs par défaut suggérées sont fournies. 5 est le nombre minimum de gènes qui doivent être dans une voie afin de maintenir la voie dans l’analyse, l’augmentation fait référence à une quantité croissante du trait mesuré (il est recommandé que l’utilisateur exécute à la fois en augmentation et en diminution, quel que soit le trait; l’interprétation des données différera pour les deux, cependant), 1000 est le nombre de fois pour échantillonner les effets pour déterminer la distribution nulle, et 2 est le nombre de cœurs utilisés pour le traitement parallèle. Modifiez le chemin d’accès à l’emplacement réel du fichier de chemins d’accès.</li>
	</ol>
</li>
</ol>7. Voir Rugplots
<ol><li>Voir Rugplots avec PAST Shiny.<ol>
<li>Une fois toutes les entrées téléchargées et définies, cliquez sur Commencer l’analyse. Une barre de progression apparaîtra et indiquera quelle étape de l’analyse a été terminée pour la dernière fois. Une fois l’analyse terminée, PAST Shiny passe à l’onglet Résultats. Un tableau des résultats sera affiché dans la colonne de gauche (intitulée « chemins ») et les Rugplots seront affichés dans la colonne de droite (intitulée « tracés »).</li>
		<li>Utilisez le curseur pour contrôler les paramètres de filtrage. Lorsque le niveau de filtrage est satisfaisant, cliquez sur le bouton Télécharger les résultats en bas à gauche pour télécharger toutes les images et tous les tableaux individuellement dans un fichier ZIP nommé avec le titre de l’analyse. Ce fichier ZIP contient la table filtrée, la table non filtrée et une image par chemin d’accès dans la table filtrée.</li>
	</ol>
</li>
</ol><img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/61268/61268fig6.jpg"> Figure 6. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig6large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/61268/61268fig7.jpg"> Figure 7. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig7large.jpg" target="_blank">Veuillez cliquer ici pour l’agrandir.</a>
<ol start="2"><li>Afficher rugplots avec PAST dans la console R<ol>
<li>Modifiez et exécutez le code suivant pour enregistrer les résultats : plot_pathways(rugplots_data, « pvalue », 0,02, « augmentation », « output_folder ») Remarque : Dans cet exemple de code, plusieurs valeurs par défaut suggérées sont fournies. pvalue fournit les données qui peuvent être utilisées pour filtrer les voies insignifiantes après qu’un seuil de signification a été choisi par l’utilisateur; 0,02 est la valeur par défaut utilisée dans le filtrage, et l’augmentation fait référence à une quantité croissante du trait mesuré (il est recommandé que l’utilisateur exécute à la fois en augmentation et en diminution, quel que soit le trait; l’interprétation des données différera pour les deux, cependant); output_folder est le dossier dans lequel les images et les tables seront écrites (ce dossier doit exister avant d’exécuter la fonction). Un tableau des résultats filtrés, des résultats non filtrés et des images individuelles pour chaque chemin d’accès dans les résultats filtrés sont écrits dans ce dossier.</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>

Les auteurs n’ont rien à divulguer.

L’un des principaux objectifs de PAST est d’apporter des analyses des voies métaboliques des données GWAS à un public plus large, en particulier pour les organismes non humains et non animaux. Les méthodes alternatives à PAST sont souvent des programmes en ligne de commande qui se concentrent sur les humains ou les animaux. La convivialité était l’un des principaux objectifs du développement de PAST, à la fois en choisissant de développer une application Shiny et en choisissant d’utiliser R et Bioconductor pour libérer l’application. Les utilisateurs n’ont pas besoin d’apprendre à compiler des programmes pour utiliser PAST.
Comme avec la plupart des types de logiciels d’analyse, les résultats de PAST ne sont aussi bons que les données d’entrée; si les données d’entrée présentent des erreurs ou sont mal formatées, PAST ne s’exécutera pas ou ne produira pas de résultats non informatifs. Il est essentiel de s’assurer que les données GWAS, les données LD, les annotations et les fichiers de chemins d’accès sont correctement formatés est essentiel pour recevoir une sortie correcte de PAST. PAST n’analyse que les marqueurs bi-alléliques et ne peut exécuter qu’un seul trait pour chaque ensemble de données d’entrée. En outre, les données GWAS produites par un mauvais génotypage ou un phénotypage incorrect ou imprécis ne sont pas non plus susceptibles de produire des résultats clairs ou reproductibles. PAST peut aider à l’interprétation biologique des résultats du GWAS, mais il est peu probable qu’il clarifie les ensembles de données chaotiques si la variation environnementale, l’erreur expérimentale ou la structure de la population n’ont pas été correctement prises en compte.
Les utilisateurs peuvent choisir de modifier certains paramètres de l’analyse, à la fois dans l’application Shiny et en transmettant ces paramètres aux fonctions de PAST dans la console R. Ces paramètres peuvent modifier les résultats rapportés par PAST, et les utilisateurs doivent faire attention lorsqu’ils les modifient à partir des valeurs par défaut. Étant donné que la LD est mesurée par les utilisateurs, généralement à l’aide du même ensemble de données de marqueurs qui a également été utilisé dans le GWAS, les mesures de LD sont spécifiques à la population. Pour toutes les études, en particulier pour les espèces autres que le maïs (en particulier les espèces autopolinisatrices, polyploïdes ou très hétérogènes), des changements dans les valeurs par défaut peuvent être justifiés.

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

Les études d’association à l’échelle du génome (GWAS) sont une méthode populaire d’étude des traits complexes et des régions génomiques qui leur sont associées1,2,3. Dans ce type d’étude, des centaines de milliers de marqueurs de polymorphisme nucléotidique unique (SNP) sont testés pour leur association avec le trait, et la signification des associations est évaluée. Les associations marqueurs-traits qui atteignent le seuil du taux de fausse découverte (FDR) (ou un autre type de seuil de signification) sont conservées pour l’étude, mais les associations vraies peuvent être filtrées. Pour les caractères polygéniques complexes, l’effet de chaque gène peut être faible (et donc filtré), et certains allèles ne sont exprimés que dans des conditions spécifiques qui pourraient ne pas être présentes dans l’étude3. Ainsi, bien que de nombreux SLP puissent être conservés comme associés au trait, chacun peut avoir un très faible effet. Trop d’appels SNP seront manquants, et une interprétation de la signification biologique et de l’architecture génétique du trait peut être incomplète et déroutante. L’analyse des voies métaboliques peut aider à résoudre certains de ces problèmes en se concentrant sur les effets combinés des gènes regroupés en fonction de leur fonction biologique4,5,6.
Plusieurs études ont été réalisées à l’aide d’une implémentation antérieure de la méthode décrite dans cet article. L’accumulation d’aflatoxines7,la résistance au ver de l’oreille du maïs8et la biosynthèse de l’huile9 ont toutes été étudiées avec la mise en œuvre précédente. Bien que ces analyses aient été couronnées de succès, le processus d’analyse était compliqué, long et fastidieux, car les outils d’analyse étaient écrits dans une combinaison de R, Perl et Bash, et le pipeline n’était pas automatisé. En raison des connaissances spécialisées nécessaires pour modifier cette méthode pour chaque analyse, une nouvelle méthode a maintenant été développée qui peut être partagée avec d’autres chercheurs.
Le Pathway Association Study Tool (PAST)10 a été conçu pour remédier aux lacunes de la méthode précédente en exigeant moins de connaissances en langages de programmation et en exécutant des analyses dans un délai plus court. Bien que la méthode ait été testée avec du maïs, PAST ne fait aucune hypothèse spécifique à l’espèce. PAST peut être exécuté via la console R, en tant qu’application Shiny, et une version en ligne devrait bientôt être disponible sur MaizeGDB.

<table><tbody><tr><td>Ordinateur</td><td>NA</td><td>NA</td><td>Tout ordinateur doté de 8 Go de RAM doit suffire</td></tr><tr><td></td><td></td><td></td><td>à l’installation de Bioconductor 3.11</td></tr></tbody></table>

un outil d’étude d’association de voies pour les analyses gwas de l’information sur les voies métaboliques

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

Récemment, une nouvelle mise en œuvre d’une méthode décrite précédemment pour interpréter les données d’étude d’association à l’échelle du génome (GWAS) à l’aide de l’analyse des voies métaboliques a été développée et publiée. L’outil d’étude pathway association (PAST) a été développé pour répondre aux préoccupations concernant la convivialité et les analyses lentes. Ce nouvel outil convivial a été publié sur Bioconductor et Github. Lors des tests, PAST a exécuté des analyses en moins d’une heure qui nécessitaient auparavant vingt-quatre heures ou plus. Dans cet article, nous présentons le protocole permettant d’utiliser l’application Shiny ou la console R pour exécuter PAST.

Si les résultats ne sont pas produits après l’exécution de l’outil logiciel PAST, vérifiez que tous les fichiers d’entrée sont correctement formatés. Une exécution réussie à l’aide des exemples de données du package PAST, qui sont basées sur un GWAS de maïs de couleur de grain, est illustrée à la figure 8. Ce tableau et l’image résultante peuvent être téléchargés à l’aide du bouton Télécharger les résultats. Un exemple de l’image téléchargée est illustré à la Figure 210. Des paramètres incorrects peuvent conduire à des résultats qui n’ont pas de sens biologique, mais la détermination de l’incorrect doit être du choix du chercheur, qui doit vérifier la validité des paramètres choisis et prendre en compte toutes les preuves connues concernant le trait d’intérêt.
La figure 910 montre le rugplot produit à partir de l’analyse des voies des résultats GWAS créés avec un panneau de maïs de 288 lignées consanguines qui avaient été phénotypées pour la couleur des grains. Cet exemple simpliste, où les phénotypes étaient soit « blancs » soit « jaunes », a été utilisé parce que la voie responsable de la création des pigments caroténoïdes jaune vif est connue et devrait être responsable de la plupart des phénotypes. Ainsi, nous nous attendions à ce que la voie de biosynthèse trans-lycopène (qui produit des caroténoïdes) soit significativement associée à la couleur du grain, ce qui est le cas. L’ID et le nom du chemin d’accès sont répertoriés en haut du graphique. L’axe horizontal du graphique classe tous les gènes qui ont été inclus dans l’analyse, disposés de gauche à droite dans l’ordre du plus grand effet sur le trait au plus petit. Cependant, seuls les gènes de la voie de biosynthèse du trans-lycopène sont marqués (en haut du graphique, sous forme de marques d’éclosion, apparaissant dans le rang génétique de leur effet par rapport à tous les autres gènes de l’analyse). Il y a 7 gènes dans cette voie. Le score d’enrichissement en cours d’exécution (ES) est tracé le long de l’axe vertical. L’ES pour chaque gène est ajouté au total courant par ordre d’effet et le total est ajusté au nombre de gènes analysés. Ainsi, le score change à mesure que l’on se déplace le long de l’axe horizontal et tend à augmenter à mesure que les gènes à effet plus importants sont inclus, mais à un moment donné, l’augmentation de l’effet est plus petite que l’ajustement pour avoir ajouté un autre gène, et le score entier commence à diminuer. Le sommet de la ligne ES en cours d’exécution est marqué d’une ligne verticale pointillée; il s’agit de l’ES pour l’ensemble du parcours et est utilisé par le programme pour déterminer si le chemin est choisi et présenté comme un rugplot.
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/61268/61268fig8.jpg"> Figure 8: Exécution terminée de PAST Shiny. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig8large.jpg" target="_blank">Veuillez cliquer ici pour afficher une version agrandie de cette figure.</a>
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/61268/61268fig9.jpg"> Figure 9: Image du chemin d’accès à partir de l’exécution terminée de PAST (ou téléchargée à partir de Shiny). Ce chiffre a été cité dans Thrash et al.10. <a href="https://www.jove.com/files/ftp_upload/61268/61268fig9large.jpg" target="_blank">Veuillez cliquer ici pour voir une version agrandie de cette figure.</a>

En exécutant l’outil d’étude pathway association (PAST), soit via l’application Shiny, soit via la console R, les chercheurs peuvent mieux comprendre la signification biologique des résultats de leur étude d’association à l’échelle du génome (GWAS) en étudiant les voies métaboliques impliquées.

Un outil d’étude d’association de voies pour les analyses GWAS de l’information sur les voies métaboliques

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

L’analyse d’association à l’échelle du génome est un outil puissant pour découvrir des gènes associés à un caractère, mais en manque beaucoup. PAST trouve plus de ces gènes et complète mieux le tableau. PAST prend les résultats des GWAS et ajoute des informations sur les voies que les cellules utilisent pour construire des métabolites et pour créer des phénotypes dans un organisme.Les gènes de ces voies sont souvent manqués par les GWAS seuls. Après avoir vérifié que les données GWAS sont délimitées par des tabulations, utilisez la zone de fichier d’association pour sélectionner un fichier d’association et la zone de sélection de fichier d’effets pour sélectionner un fichier d’effets. Pour charger des données GWAS avec PAST dans la console R, modifiez et exécutez le code comme indiqué.Pour charger des données de déséquilibre de liaison avec PAST Shiny, après avoir vérifié que les données de déséquilibre de liaison sont délimitées par des tabulations et contiennent les types de données appropriés, sélectionnez le fichier contenant les données de déséquilibre de liaison. Pour charger des données de déséquilibre de liaison avec PAST dans la console R, modifiez et exécutez le code comme indiqué. Pour attribuer des SNP à des gènes avec PAST Shiny, en localisant les annotations dans le format de caractéristique générale, sélectionnez le fichier contenant les annotations de format de caractéristique générale et modifiez les paramètres en fonction de la taille de la fenêtre et de la coupure R carré qui conviennent le mieux à l’espèce considérée comme nécessaire.Pour attribuer des SNP à des gènes avec PAST dans la console R, modifiez et exécutez le code comme indiqué. Après avoir vérifié que le fichier de parcours contient les données appropriées au format délimité par des tabulations, sélectionnez le fichier contenant les données de parcours et vérifiez que le mode est sélectionné dans les options d’analyse. Si nécessaire, modifiez le nombre de gènes qui doivent se trouver dans une voie pour conserver la voie d’analyse et modifiez le nombre de permutations utilisées pour créer la distribution nulle afin de tester la signification de l’effet.Pour découvrir les chemins significatifs avec PAST dans la console R, modifiez et exécutez le code comme indiqué. Pour afficher les parcelles de tapis avec PAST Shiny, une fois que toutes les entrées ont été téléchargées et définies, cliquez sur commencer l’analyse. Une barre de progression s’affiche pour indiquer la dernière étape de l’analyse.Une fois l’analyse terminée, PAST Shiny passera à l’onglet résultats et un tableau des résultats et les tracés du tapis s’afficheront. Utilisez le curseur pour contrôler les paramètres de filtrage. Lorsque le niveau de filtrage est satisfaisant, cliquez sur le bouton Télécharger les résultats pour télécharger individuellement toutes les images et tables dans un fichier zip nommé avec le titre de l’analyse.Pour afficher les tracés de tapis avec PAST dans la console R, modifiez et exécutez le code pour enregistrer les résultats. Ici, une exécution réussie utilisant les données d’exemple du package PAST, qui sont basées sur un labyrinthe GWAS de couleur de grain, est montrée. Et ici, un exemple d’image qui peut être téléchargée en cliquant sur le bouton de téléchargement des résultats est montré.Ce tracé de tapis a été produit à partir de l’analyse de la voie des résultats GWAS créés avec le panel de labyrinthe de 288 lignées consanguines qui avaient été phénotypiques pour la couleur du grain. Comme prévu, la voie de biosynthèse du trans-lycopène, qui produisait des caroténoïdes, était significativement associée à la couleur du grain. Les gènes et la voie de biosynthèse du trans-lycopène sont également marqués dans le rang des gènes de leur effet par rapport à tous les autres gènes de l’analyse.Le sommet de la ligne de score d’enrichissement courant pour l’ensemble du parcours est indiqué par la ligne pointillée. Il est utilisé par le programme pour déterminer si le chemin est choisi et présenté sous la forme d’un graphique de tapis. PAST sur toutes les données, mais ne produit de bons résultats que pour des données de haute qualité.Les séquences doivent contenir peu d’erreurs et les données phénotypiques doivent être recueillies dans le cadre d’un plan statistique rigoureux. Grâce aux résultats de l’étude PAST, les gènes contrôlant des points importants et des voies significatives peuvent être identifiés et éventuellement manipulés par sélection, édition de gènes et traitement médicamenteux pour améliorer un trait d’intérêt.

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

Research

Recherche

JoVE Journal

Genetics

Génétique

Combinaison de la génomique, l'analyse du gène co-expression et l'identification de composés cibles par métabolisme donner annotation des gènes fonctionnels.

Annotation de la fonction des gènes des plantes par l&#39;intermédiaire de génomique, la métabolomique et combinés informatique

PathWhiz est un outil de dessin de voie en ligne complète pour générer des voies biochimiques et biologiques. Il utilise des bases de données accessibles au public et facilement extensibles palettes constituées de composants pré-dessinés voie. Ce protocole décrit comment construire facilement de nouvelles voies, reproduire et modifier les voies existantes, et de propager les voies précédemment tirées à différents organismes.

Un outil Web pour la génération de haute qualité lisible par machine voies biologiques

Nous fournissons un protocole normalisé pour l’utilisation du gène enrichissement set analyse des données transcriptomiques pour identifier un modèle souris idéal pour la recherche translationnelle.Ce protocole peut être utilisé avec les puces à ADN et des données de séquençage de l’ARN et encore peut être étendu à d’autres données omics si des données sont disponibles.

Un protocole pour l’utilisation de gènes série d’analyses d’enrichissement pour identifier le modèle Animal approprié de recherche translationnelle

Dans ce protocole, nous présentons un flux de travail optimisé, qui combine une préparation efficace et rapide des échantillons de nombreux échantillons. En outre, nous fournissons un guide étape par étape pour réduire les variations analytiques pour l’évaluation à haut débit des études GWAS métaboliques.

Études d’association multi-omiques à grande échelle à l’échelle du génome (Mo-GWAS): Lignes directrices pour la préparation et la normalisation des échantillons

Nous présentons CorrelationCalculator et Filigree, deux outils pour la construction de réseaux basés sur les données et l’analyse des données métabolomiques.  CorrelationCalculator prend en charge la création d’un réseau d’interaction unique de métabolites basé sur des données d’expression, tandis que Filigree permet de construire un réseau différentiel, suivi d’une analyse de clustering et d’enrichissement du réseau.

CorrelationCalculator et Filigree : outils pour l’analyse de réseau basée sur les données métabolomiques

DeepOmicsAE est un flux de travail centré sur l’application d’une méthode d’apprentissage profond (c’est-à-dire un auto-encodeur) pour réduire la dimensionnalité des données multi-omiques, fournissant une base pour les modèles prédictifs et les modules de signalisation représentant plusieurs couches de données omiques.

DeepOmicsAE : Représentation de modules de signalisation dans la maladie d’Alzheimer avec l’analyse de l’apprentissage profond de la protéomique, de la métabolomique et des données cliniques

β la mort cellulaire précède le développement du diabète de type 1, et la détection de ce processus peut permettre une intervention thérapeutique précoce. Nous fournissons ici une description détaillée de la façon de mesurer les espèces d’ADN INS méthylées de manière différentielle dans le sérum humain en tant que biomarqueur de la mort cellulaire β.

Mesure de Différentiellement méthylé<em&gt; INS</em&gt; ADN Espèce dans des échantillons de sérum humain en tant que mort cellulaire Biomarker de l&#39;Islet β

Contamination pendant le séquençage génomique des organismes microscopiques demeure un grand problème. Ici, nous montrons une méthode pour séquencer le génome d’un tardigrade un seul spécimen avec aussi peu que 50 pg d’ADN génomique sans amplification du génome entier pour minimiser le risque de contamination.

Ultra basse génome d’entrée séquençage bibliothèque préparation à partir d’un seul spécimen Tardigrade

Nous présentons ici un protocole pour permettre une quantification précise de la méthylation de l’ADN (ADN mitochondrial) mitochondriale. Dans ce protocole, nous décrivons une digestion enzymatique de l’ADN avec BamHI couplé avec un pipeline d’analyse bio-informatique qui peut être utilisé pour éviter une surestimation des niveaux de méthylation de l’ADNmt causée par la structure secondaire de l’ADNmt.

Méthodologie pour une détection précise de la méthylation de l’ADN Mitochondrial

L’identification des interactions physiques entre les gènes et les éléments de régulation est difficile, mais a été facilitée par les méthodes de capture de conformation du chromosome. Cette modification du protocole 4C-seq atténue la partialité de la PCR en réduisant au minimum l’amplification excessive des modèles PCR et maximise la mappability se lit en y incorporant une étape digest d’enzyme de restriction addition.

High-throughput Identification des séquences régulatrices du gène à l’aide de prochaine génération séquençage du Chromosome circulaire Conformation Capture (4C-seq)

Découvrez l'analyse d'association de voies métaboliques avec l'outil PAST pour compléter les résultats GWAS. Cette méthode bioinformatique identifie les gènes et voies significatifs liés aux traits d'intérêt.

Un outil d’étude d’association de voies pour les analyses GWAS de l’information sur les voies métaboliques

Explorer plus de vidéos