REMARQUE: Dans ce protocole, l’utilisation de JUMPn est illustrée par l’utilisation d’un ensemble de données publié de profilage du protéome entier au cours de la différenciation des cellules B quantifié par le réactif d’étiquette isobare TMT27.
1. Configuration du logiciel JUMPn
REMARQUE: Deux options sont fournies pour configurer le logiciel JUMPn: (i) installation sur un ordinateur local pour un usage personnel; et (ii) le déploiement de JUMPn sur un serveur Shiny distant pour plusieurs utilisateurs. Pour l’installation locale, un ordinateur personnel avec accès à Internet et ≥4 Go de RAM est suffisant pour exécuter l’analyse JUMPn pour un jeu de données avec une petite taille d’échantillon (n < 30); une RAM plus importante (p. ex., 16 Go) est nécessaire pour l’analyse de grandes cohortes (p. ex., n = 200 échantillons).
- Installez le logiciel sur un ordinateur local. Après l’installation, autorisez le navigateur Web à lancer JUMPn et laissez l’analyse s’exécuter sur l’ordinateur local.
- Installez anaconda42 ou miniconda43 en suivant les instructions en ligne.
- Téléchargez le code sourceJUMPn 41. Double-cliquez pour décompresser le fichier téléchargé JUMPn_v_1.0.0.zip; un nouveau dossier nommé JUMPn_v_1.0.0 sera créé.
- Ouvrez le terminal de ligne de commande. Sous Windows, utilisez l’invite Anaconda. Sous MacOS, utilisez l’application Terminal intégrée.
- Créer l’environnement JUMPn Conda : obtenez le chemin absolu du dossier JUMPn_v_1.0.0 (par exemple, /path/to/JUMPn_v_1.0.0). Pour créer et activer un environnement Conda vide, tapez les commandes suivantes sur le terminal
conda create -p /path/to/JUMPn_v_1.0.0/JUMPn -y
conda activer /path/to/JUMPn_v_1.0.0/JUMPn
- Installer les dépendances JUMPn : Installez R (sur le terminal, tapez conda install -c conda-forge r=4.0.0 -y), remplacez le répertoire actuel par le dossier JUMPn_v_1.0.0 (sur le terminal, tapez cd path/to/JUMPn_v_1.0.0) et installez les packages de dépendance (sur le terminal, tapez Rscript bootstrap. R)
- Lancez JUMPn sur le navigateur Web : Remplacez le répertoire actuel par le dossier d’exécution (sur le terminal, tapez cd execution) et lancez JUMPn (sur le terminal, tapez R -e « shiny::runApp() »)
- Une fois ce qui précède exécuté, l’écran du terminal apparaîtra Écoute sur http://127.0.0.1: XXXX (ici XXXX indique 4 nombres aléatoires). Copiez et collez http://127.0.0.1:XXXX sur le navigateur Web, sur lequel la page d’accueil JUMPn apparaîtra (Figure 2).
- Déploiement sur Shiny Server. Des exemples de Shiny Server incluent le serveur commercial shinyapps.io ou tout serveur Shiny pris en charge par un établissement.
- Téléchargez et installez RStudio en suivant les instructions44.
- Obtenez l’autorisation de déploiement pour le serveur Shiny. Pour le serveur shinyapps.io, configurez le compte d’utilisateur en suivant l’instruction45. Pour le serveur Shiny institutionnel, contactez l’administrateur du serveur pour demander des autorisations.
- Téléchargez le code sourceJUMPn 41 sur la machine locale ; l’installation n’est pas nécessaire. Ouvrez l’un ou l’autre serveur. R ou ui. R dans RStudio et cliquez sur le menu déroulant Publier sur le serveur en haut à droite de l’IDE RStudio.
- Dans le panneau Publier dans le compte , tapez l’adresse du serveur. Appuyez sur le bouton Publier . Le déploiement réussi est validé lors de la redirection automatique de RStudio vers le serveur RShiny où l’application a été déployée.
2. Exécution de la démo à l’aide d’un exemple de jeu de données
REMARQUE: JUMPn propose une démonstration à l’aide de l’ensemble de données protéomique des cellules B publié. L’exécution de démonstration illustre un flux de travail rationalisé qui prend la matrice de quantification des protéines exprimées différentiellement comme entrée et effectue le regroupement de co-expression, l’enrichissement des voies et l’analyse séquentielle du réseau PPI.
- Sur la page d’accueil de JUMPn (Figure 2), cliquez sur le bouton Commencer l’analyse pour lancer l’analyse JUMPn.
- Dans le coin inférieur gauche de la page Commencer l’analyse (Figure 3), cliquez sur le bouton Télécharger les données protéomiques des cellules B de démonstration ; une boîte de dialogue apparaîtra pour notifier la réussite du téléchargement des données.
- Dans le coin inférieur droit de la page, cliquez sur le bouton Soumettre l’analyse JUMPn pour lancer l’exécution de la démonstration en utilisant les paramètres par défaut; une barre de progression apparaîtra indiquant le déroulement de l’analyse. Attendez que la barre de progression soit remplie (3 min prévus).
- Une fois l’exécution de la démo terminée, une boîte de dialogue apparaîtra avec le message de réussite de l’exécution et le chemin absolu vers le dossier de résultats. Cliquez sur Continuer vers les résultats pour continuer.
- La page Web guidera d’abord l’utilisateur vers les résultats du cluster de co-expression par WGCNA. Cliquez sur Afficher les résultats dans la fenêtre de dialogue pour continuer.
- Trouvez les modèles de co-expression des protéines à gauche de la page de résultats 1 : Page de sortie WGCNA . Cliquez sur la liste déroulante Sélectionner le format d’expression pour naviguer entre deux formats de figure :
- Sélectionnez Tendances pour afficher le graphique des tendances, chaque ligne représentant l’abondance de protéines individuelles dans les échantillons. La couleur de chaque ligne représente la proximité du modèle d’expression avec le consensus du cluster de co-expression (c’est-à-dire « eigengene » tel que défini par l’algorithme WGCNA).
- Sélectionnez Boxplot pour afficher des modèles de co-expression au format boxplot pour chaque exemple.
- Consultez la carte thermique d’enrichissement des voies/ontologies à droite de la page de sortie WGCNA. Les voies les plus enrichies pour chaque cluster sont affichées ensemble dans une carte thermique, l’intensité des couleurs reflétant la valeur p ajustée de Benjamini-Hochberg.
- Faites défiler la page Web vers le bas pour afficher le modèle d’expression des protéines individuelles.
- Utilisez la liste déroulante Sélectionnez le cluster de co-expression pour afficher les protéines de chaque cluster (la valeur par défaut est cluster 1). Sélectionnez une protéine spécifique dans le tableau, sur laquelle le graphique à barres sous le tableau sera automatiquement mis à jour pour refléter son abondance de protéines.
- Recherchez des noms de protéines spécifiques à l’aide de la zone de recherche sur le côté droit du tableau pour une protéine spécifique.
- Pour afficher les résultats des IPP, cliquez sur la page résultats 2 : Sortie PPI en haut.
- Cliquez sur Sélectionner le cluster de co-expression pour afficher les résultats d’un cluster de co-expression spécifique (la valeur par défaut est le cluster 1). Les affichages de tous les panneaux de figure sur cette page seront mis à jour pour le cluster nouvellement sélectionné.
- Affichez les réseaux PPI pour le cluster de coexpression sélectionné dans le panneau de figure de gauche :
- Cliquez sur la liste déroulante Sélectionner par groupe pour mettre en surbrillance les modules PPI individuels du réseau. Cliquez sur la liste déroulante Sélectionner un format de disposition réseau pour modifier la disposition du réseau (la valeur par défaut est par Fruchterman Reingold).
- Utilisez la souris et le trackpad pour effectuer les étapes 2.11.3 à 2.11.5.
- Effectuez un zoom avant ou arrière sur le réseau PPI selon les besoins. Les noms de gènes de chaque nœud du réseau seront affichés lorsqu’ils seront suffisamment zoomés.
- Lorsque vous effectuez un zoom avant, sélectionnez et cliquez sur une certaine protéine pour mettre en évidence cette protéine et ses voisins du réseau.
- Faites glisser un certain nœud (protéine) dans le réseau pour changer sa position dans la disposition; ainsi, la disposition du réseau peut être réorganisée par l’utilisateur.
- Dans le panneau de droite de la page de résultats de l’IPP, affichez les informations au niveau du cluster de coexpression qui facilitent l’interprétation des résultats de l’IPP :
- Affichez le modèle de coexpression du cluster sélectionné sous forme de diagramme de boîte par défaut.
- Cliquez sur la liste déroulante Sélectionner le format d’expression pour plus d’informations ou affiche comme mentionné aux étapes 2.12.3-2.12.5.
- Sélectionnez Tendances pour afficher le tracé des tendances pour le modèle de co-expression.
- Sélectionnez Path Barplot pour afficher les chemins enrichis de manière significative pour le cluster de co-expression.
- Sélectionnez Tracé circulaire de trajectoire pour afficher les voies considérablement enrichies pour le groupe de coexpression au format de tracé circulaire.
- Faites défiler la page de résultats 2 : Sortie PPI vers le bas pour afficher les résultats au niveau du module PPI individuel. Cliquez sur la liste déroulante Sélectionner le module pour sélectionner un module PPI spécifique à afficher (Cluster1: Module 1 est affiché par défaut).
- Affichez le module PPI sur le panneau de gauche. Pour manipuler l’affichage réseau, suivez les étapes 2.11.2 à 2.11.5.
- Affichez les résultats de l’enrichissement des voies/ontologies sur le panneau de droite. Cliquez sur la liste déroulante Sélectionner le style d’annotation de chemin pour plus d’informations et affiche :
- Sélectionnez Barplot pour afficher les voies enrichies de manière significative pour le module PPI sélectionné.
- Sélectionnez Tracé circulaire pour afficher les tracés significativement enrichis pour le module PPI sélectionné au format d’un tracé circulaire.
- Sélectionnez Heatmap pour afficher les voies significativement enrichies et les noms de gènes associés à partir du module PPI sélectionné.
- Sélectionnez Tableau pour afficher les résultats détaillés de l’enrichissement des voies, y compris le nom des voies/termes d’ontologie, les noms des gènes et la valeur P par le test exact de Fisher.
- Affichez le tableau de publication sous forme de feuille de calcul : suivez le chemin absolu (imprimé en haut des deux pages de résultats) et recherchez le tableau de feuille de calcul de publication nommé ComprehensiveSummaryTables.xlsx.
3. Préparation du fichier d’entrée et téléchargement sur JUMPn
NOTE: JUMPn prend comme entrée la matrice de quantification des protéines exprimées différentiellement (méthode supervisée) ou des protéines les plus variables (méthode non supervisée). Si l’objectif du projet est de comprendre les protéines modifiées dans plusieurs conditions (p. ex., différents groupes de maladies ou analyse chronologique du processus biologique), la méthode supervisée d’analyse de l’DE est préférable; sinon, une approche non supervisée de sélection des protéines les plus variables peut être utilisée à des fins exploratoires.
- Générez le tableau de quantification des protéines, avec chaque protéine sous forme de lignes et chaque échantillon sous forme de colonnes. Atteignez cet objectif grâce à la suite logicielle de protéomique moderne basée sur la spectrométrie de masse (par exemple, JUMP suite 13,14,39, Proteome Discoverer, Maxquant15,46).
- Définissez la variable protéome.
- Utilisez les résultats d’analyse statistique fournis par la suite logicielle de protéomique pour définir les protéines exprimées différentiellement (DE) (par exemple, avec une valeur p ajustée < 0,05).
- Alternativement, les utilisateurs peuvent suivre l’exemple du code R47 pour définir soit DE, soit la plupart des protéines variables.
- Formatez le fichier d’entrée à l’aide de la variable protéome définie.
REMARQUE : Le format de fichier d’entrée requis (Figure 4) comprend une ligne d’en-tête ; les colonnes comprennent l’accession des protéines (ou toute identification unique), gn (symboles officiels des gènes), la description des protéines (ou toute information fournie par l’utilisateur), suivie de la quantification des protéines d’échantillons individuels.
- Suivez l’ordre des colonnes spécifié à l’étape 3.1, mais les noms de colonne de l’en-tête sont flexibles pour l’utilisateur.
- Pour le protéome quantifié TMT (ou similaire), utilisez l’intensité du rapporteur TMT résumée comme valeurs de quantification d’entrée. Pour les données sans étiquette, utilisez des comptes spectraux normalisés (p. ex., NSAF48) ou une méthode fondée sur l’intensité (p. ex., intensité LFQ ou intensité des protéines iBAQ déclarée par Maxquant46).
- Les valeurs manquantes sont autorisées pour l’analyse JUMPn. Assurez-vous de les étiqueter comme NA dans la matrice de quantification. Cependant, il est recommandé de n’utiliser que des protéines avec quantification dans plus de 50% des échantillons.
- Enregistrez le fichier d’entrée résultant au format .txt, .xlsx ou .csv (les trois sont pris en charge par JUMPn).
- Télécharger le fichier d’entrée :
- Cliquez sur le bouton Navigateur et sélectionnez le fichier de saisie (Figure 3, panneau de gauche) ; le format de fichier (xlsx, csv et txt sont pris en charge) sera automatiquement détecté.
- Si le fichier d’entrée contient des valeurs de quantification de type intensité (par exemple, celles générées par JUMP suite39) ou de type ratio (par exemple, à partir de Proteome Discoverer), sélectionnez Oui pour l’option Exécuter Log2-Transformation des données ; sinon, les données ont peut-être déjà été transformées par journal, alors sélectionnez Non pour cette option.
4. Analyse du clustering de co-expression
NOTE: Notre groupe 25,26,27 et d’autres 28,29,31 ont prouvé WGCNA 49 une méthode efficace pour l’analyse groupée de co-expression de la protéomique quantitative. JUMPn suit une procédure en 3 étapes pour l’analyse WGCNA25,50: (i) définition initiale des groupes gènes/protéines de co-expression par coupe dynamique d’arbre51 basée sur la matrice de chevauchement topologique (TOM; déterminée par des similitudes de quantification entre gènes/protéines); ii) la fusion de groupes similaires pour réduire la redondance (sur la base du dendrogramme des similitudes degènes propres); et (iii) l’attribution finale de gènes/protéines à chaque grappe qui dépassent le seuil minimal de corrélation de Pearson.
- Configurez les paramètres WGCNA (Figure 3, panneau central). Les trois paramètres suivants contrôlent respectivement les trois étapes :
- Définissez la taille minimale du cluster sur 30. Ce paramètre définit le nombre minimal de protéines requises pour chaque groupe de co-expression dans l’étape initiale (i) de la coupe dynamique dynamique hybride basée sur TOM. Plus la valeur est grande, plus le nombre de clusters renvoyés par l’algorithme est petit.
- Définissez la distance minimale du cluster sur 0,2. L’augmentation de cette valeur (p. ex., de 0,2 à 0,3) peut entraîner une plus grande fusion de grappes au cours de l’étape ii), ce qui entraîne une diminution du nombre de grappes.
- Définissez le kME minimum sur 0,7. Les protéines seront affectées au groupe le plus corrélé défini à l’étape ii), mais seules les protéines dont la corrélation de Pearson dépasse ce seuil seront conservées. Les protéines qui échouent à cette étape ne seront affectées à aucun groupe (cluster « NA » pour les protéines défaillantes dans le rapport final).
- Lancez l’analyse. Il existe deux façons de soumettre l’analyse de clustering de co-expression :
- Cliquez sur le bouton Soumettre l’analyse JUMPn dans le coin inférieur droit pour lancer automatiquement l’analyse complète de WGCNA, suivie d’une analyse du réseau PPI.
- Vous pouvez également choisir d’exécuter l’étape WGCNA uniquement (en particulier à des fins de réglage des paramètres ; voir les étapes 4.2.3 à 4.2.4) :
- Cliquez sur le bouton Paramètres avancés au bas de la page Commencer l’analyse ; une nouvelle fenêtre de paramètre apparaîtra. Dans le widget du bas, Sélectionnez le mode d’analyse, sélectionnez WGCNA uniquement, puis cliquez sur Ignorer pour continuer.
- Sur la page Commencer l’analyse , cliquez sur le bouton Soumettre l’analyse JUMPn .
- Dans les deux cas ci-dessus, une barre de progression apparaîtra lors de la soumission de l’analyse.
REMARQUE : Une fois l’analyse terminée (généralement < 1 min pour l’analyse WGCNA uniquement et <3 min pour l’analyse complète), une boîte de dialogue s’affiche avec un message d’exécution réussie et le chemin d’accès absolu au dossier de résultats.
- Examinez les résultats du WGCNA comme illustré aux étapes 2.4 à 2.8 (Figure 5). Notez que le chemin absolu vers le fichier co_exp_clusters_3colums.txt est mis en surbrillance en haut de la page de résultats : Sortie WGCNA pour enregistrer l’appartenance au cluster de chaque protéine et l’utiliser comme entrée pour l’analyse PPI uniquement .
- Dépannage. Les trois cas courants suivants sont discutés. Une fois les paramètres mis à jour comme indiqué ci-dessous, suivez les étapes 4.2.2-4.2.4 pour générer de nouveaux résultats WGCNA.
- Si un modèle de co-expression important est attendu des données mais manqué par l’algorithme, suivez les étapes 4.4.2-4.4.4
- Un groupe manquant est particulièrement probable pour les petits groupes de coexpression, c’est-à-dire seulement un nombre limité (p. ex., <30) de protéines présentant ce modèle. Avant la réanalyse, réexaminez le fichier d’entrée de la matrice de quantification des protéines et localisez plusieurs protéines témoins positives qui adhèrent à cet important modèle de co-expression.
- Pour sauver les petits clusters, diminuez la taille minimale du cluster (par exemple, 10 ; la taille du cluster inférieure à 10 peut ne pas être robuste et n’est donc pas recommandée) et diminuez la distance minimale du cluster (par exemple, 0,1 ; ici, le paramètre 0 est également autorisé, ce qui signifie que la fusion automatique des clusters sera ignorée).
- Après avoir exécuté l’étape de clustering de co-expression avec les paramètres mis à jour, vérifiez d’abord si le cluster est sauvé des diagrammes de modèle de co-expression, puis vérifiez les contrôles positifs en recherchant leurs accessions de protéines dans Detailed Protein Quantification (assurez-vous de sélectionner le cluster de co-expression approprié dans le widget déroulant de gauche avant la recherche).
REMARQUE: Plusieurs itérations de réglage et de réexécution des paramètres peuvent être nécessaires pour le sauvetage.
- S’il y a trop de protéines qui ne peuvent être attribuées à aucun groupe, suivez les étapes 4.4.6 à 4.4.7.
REMARQUE: Habituellement, un petit pourcentage (généralement <10%) de protéines peut ne pas être attribué à un groupe car il peut s’agir de protéines aberrantes qui n’ont suivi aucun des modèles d’expression courants de l’ensemble de données. Cependant, si ce pourcentage est significatif (p. ex., >30 %), cela suggère qu’il existe d’autres modèles de coexpression qui ne peuvent être ignorés.
- Diminuez les paramètres Taille minimale du cluster et Distance minimale du cluster pour atténuer cette situation en détectant les « nouveaux » clusters de coexpression.
- De plus, diminuez le paramètre de corrélation minimale de Pearson (kME) pour réduire ces protéines du « cluster NA ».
REMARQUE: Le réglage de ce paramètre ne générera pas de nouveaux clusters, mais augmentera plutôt la taille des clusters « existants » en acceptant plus de protéines précédemment échouées avec le seuil inférieur; cependant, cela augmentera également l’hétérogénéité de chaque groupe, car plus de protéines bruyantes sont maintenant autorisées.
- Deux groupes ont une très petite différence de modèles; les fusionner en un seul cluster en suivant les étapes 4.4.9 à 4.4.11.
- Augmentez le paramètre Distance minimale du cluster pour résoudre le problème.
- Cependant, dans certaines situations, l’algorithme peut ne jamais renvoyer le modèle souhaité; dans un tel instant, ajustez ou modifiez manuellement l’appartenance au cluster dans le fichier co_exp_clusters_3colums.txt (fichier de l’étape 4.3) à fusionner.
- Prenez le fichier post-édité comme entrée pour l’analyse du réseau PPI en aval. En cas d’édition manuelle, justifiez les critères d’affectation de cluster et enregistrez la procédure d’édition manuelle.
5. Analyse du réseau d’interaction protéine-protéine
REMARQUE : En superposant des clusters de co-expression sur le réseau PPI, chaque cluster de co-expression est stratifié en modules PPI plus petits. L’analyse est effectuée pour chaque groupe de co-expression et comprend deux étapes : dans la première étape, JUMPn superpose des protéines du cluster de co-expression au réseau IPP et trouve tous les composants connectés (c’est-à-dire plusieurs groupes de nœuds/protéines connectés; à titre d’exemple, voir la figure 6A); ensuite, des communautés ou des modules (de nœuds densément connectés) seront détectés pour chaque composant connecté de manière itérative à l’aide de la méthode52 de la matrice de chevauchement topologique (TOM).
- Configurez les paramètres pour l’analyse du réseau PPI (Figure 3, panneau de droite).
- Définissez la taille minimale du module PPI sur 2. Ce paramètre définit la taille minimale des composants déconnectés à partir de la première étape de l’analyse. Tout composant plus petit que le paramètre spécifié sera supprimé des résultats finaux.
- Définissez la taille maximale du module PPI sur 40. Les composants volumineux et déconnectés qui dépassent ce seuil feront l’objet d’une analyse TOM de deuxième étape. L’analyse de la deuxième étape divisera davantage chaque grand composant en modules plus petits: chaque module contient probablement des protéines plus densément connectées que le composant d’origine dans son ensemble.
- Lancez l’analyse. Il existe deux façons de soumettre l’analyse du réseau PPI :
- Appuyez sur le bouton Soumettre l’analyse JUMPn pour effectuer automatiquement l’analyse PPI après l’analyse WGCNA par défaut.
- Vous pouvez également télécharger des résultats de cluster de co-expression personnalisés et effectuer une analyse PPI uniquement en suivant les étapes 5.2.3 à 5.2.5.
- Préparez le fichier d’entrée en suivant le format du fichier co_exp_clusters_3colums.txt (voir la sous-section 4.4).
- Cliquez sur le bouton Paramètres avancés au bas de la page Commencer l’analyse ; une nouvelle fenêtre de paramètre apparaîtra. Dans la session supérieure Upload Co-Expression Cluster Result for 'PPI Only' Analysis, cliquez sur Browser pour télécharger le fichier d’entrée préparé à l’étape 5.2.3.
- Dans le widget du bas, Sélectionnez le mode d’analyse, sélectionnez PPI uniquement, puis cliquez sur Ignorer pour continuer. Sur la page Commencer l’analyse , cliquez sur le bouton Soumettre l’analyse JUMPn .
- Une fois l’analyse terminée (généralement <3 min), examinez les résultats de l’IPP comme illustré aux étapes 2.10-2.15 (Figure 6).
- Étape avancée en option) Ajustez la modularisation PPI en réglant les paramètres:
- Augmentez le paramètre Taille maximale du module pour autoriser plus de protéines incluses dans les résultats de l’IPP. Téléchargez un réseau PPI personnalisé pour couvrir les interactions non documentées, en suivant les étapes 5.4.2 à 5.4.3.
- Cliquez sur le bouton Paramètres avancés au bas de la page Commencer l’analyse ; une nouvelle fenêtre de paramètre apparaîtra. Préparez le fichier PPI personnalisé, qui contient trois colonnes au format , C onnection et ; ici sont présentés par les noms de gènes officiels de chaque protéine.
- Dans Télécharger une base de données PPI, cliquez sur le bouton Parcourir pour télécharger le fichier PPI personnalisé.
6. Analyse de l’enrichissement des voies
REMARQUE : Les structures hiérarchiques dérivées de JUMPn des clusters de co-expression et des modules PPI à l’intérieur sont automatiquement annotées avec des voies surreprésentées à l’aide du test exact de Fisher. Les bases de données de voies/topologies utilisées incluent Gene Ontology (GO), KEGG, Hallmark et Reactome. Les utilisateurs peuvent utiliser des options avancées pour télécharger des bases de données personnalisées pour l’analyse (par exemple, dans le cas de l’analyse de données provenant d’espèces non humaines).
- Par défaut, l’analyse d’enrichissement des voies est lancée automatiquement avec le clustering de co-expression et l’analyse du réseau PPI.
- Voir les résultats de l’enrichissement de la voie :
- Suivez les étapes 2.7, 2.12 et 2.15 pour visualiser différents formats sur les pages de résultats. Affichez les résultats détaillés dans le tableau de publication de feuille de calcul dans le fichier .xlsx ComplèteSummaryTables (étape 2.16).
- (Étape avancée facultative) Téléchargez une base de données personnalisée pour l’analyse de l’enrichissement des voies :
- Préparez le fichier d’arrière-plan des gènes, qui contient généralement les noms officiels des gènes de tous les gènes d’une espèce.
- Préparez le fichier de bibliothèque d’ontologie en suivant les étapes 6.3.3 à 6.3.4.
- Téléchargez les fichiers de la bibliothèque d’ontologie à partir de sites Web publics, notamment EnrichR53 et MSigDB54. Par exemple, téléchargez l’ontologie de la drosophile à partir du site Web EnrichR55.
- Modifiez le fichier téléchargé pour le format requis avec deux colonnes: le nom de la voie comme première colonne, puis les symboles officiels du gène (séparés par « / ») comme deuxième colonne. Le format de fichier détaillé est décrit dans la page d’aide du logiciel JUMPn R shiny.
REMARQUE : Recherchez des exemples de fichiers d’arrière-plan génétique et de bibliothèque d’ontologie (en utilisant Drosophila comme instance) dans le site GitHubJUMPn 56.
- Cliquez sur le bouton Paramètres avancés au bas de la page Commencer l’analyse; une nouvelle fenêtre de paramètre apparaîtra.
- Trouvez Télécharger un fichier d’arrière-plan pour l’élément Analyse d’enrichissement de chemin et cliquez sur Navigateur pour télécharger le fichier d’arrière-plan préparé à l’étape 6.3.1. Ensuite, dans la session, sélectionnez l’arrière-plan à utiliser pour l’analyse de l’enrichissement des voies, cliquez sur Arrière-plan fourni par l’utilisateur.
- Trouvez Télécharger un fichier de bibliothèque d’ontologie pour l’élément Analyse d’enrichissement de cheminement et cliquez sur Navigateur pour télécharger le fichier de bibliothèque d’ontologie préparé aux étapes 6.3.2-6.3.4. Ensuite, dans la session, Sélectionnez Bases de données pour l’analyse d’enrichissement des voies, cliquez sur Base de données fournie par l’utilisateur dans .xlsx Format.
- Cliquez sur le bouton Soumettre l’analyse JUMPn dans le coin inférieur droit pour lancer l’analyse à l’aide de la base de données personnalisée.
7. Analyse de l’ensemble de données avec un échantillon de grande taille
REMARQUE: JUMPn prend en charge l’analyse d’un ensemble de données avec une grande taille d’échantillon (jusqu’à 200 échantillons testés). Pour faciliter la visualisation d’un échantillon de grande taille, un fichier supplémentaire (appelé « méta fichier ») qui spécifie le groupe d’échantillons est nécessaire pour faciliter l’affichage des résultats du clustering de co-expression.
- Préparez et téléchargez le méta-fichier.
- Préparez le métafique qui spécifie les informations de groupe (par exemple, les groupes témoins et de maladies) pour chaque échantillon en suivant les étapes 7.1.2 à 7.1.3.
- Assurez-vous que le métafichier contient au moins deux colonnes : la colonne 1 doit contenir les noms d’échantillons identiques aux noms de colonne et à l’ordre du fichier de matrice de quantification des protéines (tel que préparé à l’étape 3.3) ; La colonne 2 sera utilisée pour l’affectation de groupe pour n’importe quel nombre de fonctionnalités définies par l’utilisateur. Le nombre de colonnes est flexible.
- Assurez-vous que la première ligne du méta-fichier contient les noms de colonne de chaque colonne ; à partir de la deuxième rangée, des échantillons individuels de groupes ou d’autres caractéristiques (p. ex., sexe, âge, traitement, etc.) doivent être énumérés.
- Téléchargez le méta-fichier en cliquant sur le bouton Paramètres avancés au bas de la page Commencer l’analyse ; une nouvelle fenêtre de paramètre apparaîtra. Passez à l’étape 7.1.5
- Trouvez Télécharger un élément Meta File et cliquez sur Navigateur pour télécharger le fichier d’arrière-plan. Si le format inattendu ou les noms d’échantillons non correspondants sont détectés par JUMPn, un message d’erreur apparaîtra pour une mise en forme ultérieure du méta-fichier (étapes 7.1.1 à 7.1.3).
- Ajustez les paramètres pour l’analyse de clustering de co-expression : définissez la corrélation minimale de Pearson sur 0,2. Ce paramètre doit être assoupli en raison de la taille plus grande de l’échantillon.
- Cliquez sur le bouton Soumettre l’analyse JUMPn dans le coin inférieur droit pour soumettre l’analyse.
- Afficher les résultats de l’analyse : toutes les sorties de données sont les mêmes, à l’exception de l’affichage des modèles de cluster de co-expression.
- Dans la page Résultats Page 1 : Sortie WGCNA , visualisez les clusters de co-expression sous forme de diagrammes en boîte avec des échantillons stratifiés par les groupes d’échantillons ou les entités définis par l’utilisateur. Chaque point du diagramme représente le gène propre (c’est-à-dire le modèle de consensus du cluster) calculé par l’algorithme WGCNA.
- Si l’utilisateur a fourni plusieurs fonctionnalités (par exemple, âge, sexe, traitement, etc.) pour regrouper les échantillons, cliquez sur la liste déroulante Sélectionner le format d’expression pour sélectionner une autre fonctionnalité pour regrouper les échantillons.