Divergence du microbiote racinaire dans différents habitats sur la base de réseaux de corrélation pondérés

L’analyse du réseau a été appliquée pour évaluer l’association de diverses communautés microbiennes écologiques, telles que le sol, l’eau et la rhizosphère. Présenté ici est un protocole sur la façon d’utiliser l’algorithme WGCNA pour analyser différents réseaux de co-occurrence qui peuvent se produire dans les communautés microbiennes en raison de différents environnements écologiques.

Le microbiome racinaire joue un rôle important dans la croissance des plantes et l’adaptation environnementale. L’analyse des réseaux est un outil important pour étudier les communautés, qui peut explorer efficacement la relation d’interaction ou le modèle de cooccurrence de différentes espèces microbiennes dans différents environnements. Le but de ce manuscrit est de fournir des détails sur la façon d’utiliser l’algorithme de réseau de corrélation pondéré pour analyser différents réseaux de cooccurrence qui peuvent se produire dans des communautés microbiennes en raison de différents environnements écologiques. Toutes les analyses de l’expérience sont effectuées dans le package WGCNA. WGCNA est un package R pour l’analyse de réseau de corrélation pondérée. Les données expérimentales utilisées pour démontrer ces méthodes étaient des données de communauté microbienne de la base de données NCBI (National Center for Biotechnology Information) pour trois niches du système racinaire du riz (Oryza sativa). Nous avons utilisé l’algorithme du réseau de corrélation pondéré pour construire des réseaux de co-abondance de communauté microbienne dans chacune des trois niches. Ensuite, des réseaux de co-abondance différentielle entre le sol de l’endosphère, du rhizoplan et de la rhizosphère ont été identifiés. En outre, les genres de base en réseau ont été obtenus par le paquet « WGCNA », qui joue un rôle régulé important dans les fonctions de réseau. Ces méthodes permettent aux chercheurs d’analyser la réponse du réseau microbien aux perturbations environnementales et de vérifier différentes théories de la réponse écologique microbienne. Les résultats de ces méthodes montrent que les réseaux microbiens différentiels significatifs identifiés dans le sol de l’endosphère, du rhizoplan et de la rhizosphère du riz.

La recherche sur le microbiome a des implications importantes pour la compréhension et la manipulation des processus écosystémiques^1,2. Les populations microbiennes sont interconnectées par des réseaux écologiques en interaction, dont les caractéristiques peuvent affecter la réponse des micro-organismes aux changements environnementaux^3,4. De plus, les propriétés de ces réseaux affectent la stabilité des communautés microbiennes et sont étroitement associées à la fonction du sol⁵. L’analyse pondérée du réseau de corrélation génique a maintenant été largement appliquée à la recherche sur la relation entre les gènes et les communautés microbiennes⁶. Les études précédentes se sont principalement concentrées sur les associations entre les réseaux de différents gènes ou populations et le monde extérieur⁷. Cependant, les différences dans les réseaux de corrélation formés par les populations microbiennes dans différentes conditions environnementales ont été à peine étudiées. Le but de la recherche présentée dans cet article est de fournir des informations et des détails sur la mise en œuvre rapide de l’algorithme WGCNA pour construire un réseau de cooccurrence d’échantillons de microbiome collectés dans différentes conditions environnementales. Sur la base des résultats de l’analyse, nous avons évalué la composition et les différences de la population et discuté plus en détail de la relation entre les différentes populations microbiennes. Le flux de base suivant de l’algorithme de réseau de corrélation^{pondérée 8} a été appliqué. Tout d’abord, une matrice de similitude devait être construite en calculant le coefficient de corrélation de Pearson entre les profils d’expression des unités taxonomiques opérationnelles (OTU). Ensuite, les paramètres des fonctions d’adjacence (les fonctions de puissance ou d’adjacence sigmoïde) ont été adoptés avec un critère de topologie sans échelle, la matrice de similarité a été transformée en matrice d’adjacence, et chaque réseau de cooccurrence correspondait à une matrice d’adjacence. Nous avons utilisé le clustering hiérarchique de liaison moyenne couplé à la dissimilarité basée sur TOM pour regrouper les OTU avec des profils d’expression cohérents en modules. De plus, nous avons calculé la relation entre les statistiques conservatrices et les modules d’analyse des paramètres associés, identifiant enfin l’OTU du hub dans le module. Ces méthodes sont particulièrement adaptées à l’analyse des différences dans les structures de réseau entre diverses populations microbiennes dans des conditions environnementales divergentes. Dans ce manuscrit, nous avons décrit en détail la méthode de développement des réseaux de co-expression, l’analyse des dissemblances entre les modules, et avons fourni un bref aperçu des étapes de la procédure appliquée pour obtenir les espèces centrales dans différents réseaux de modules.

1. Téléchargement des données

1. Téléchargez les données de l’adhésion PRJNA386367 à la base de données NCBI. À partir des données de l’accession PRJNA386367, sélectionnez les données sur le microbiome de la rhizosphère, du rhizoplan et de l’endosphère des plants de riz cultivés pendant 14 semaines dans une rizière submergée à Arbuckle, en Californie, en 2014.
   REMARQUE: Les données sur le microbiome de la rhizosphère, du rhizoplan et de l’endosphère ont été présentées par le tableau des OTU dans l’accession PRJNA386367.

2. Détermination optimale de la valeur de puissance

Remarque : Le package WGCNA contient tous les paramètres fonctionnels suivants. WGCNA est un package R pour l’analyse de réseau de corrélation pondérée. Les lignes de commande clés se réfèrent au Supplément S1.

1. Dans l’environnement linguistique R, ouvrez le logiciel Rstudio et installez le package WGCNA.
2. Chargez les données et utilisez la fonction goodSamplesGenes pour vérifier l’exactitude des données. Exécutez les lignes de commande :
   « gsg = goodSamplesGenes(datExpr0, verbeux = 3)
   gsg$allOK »
   Cliquez sur Exécuter.
3. Vérifiez les valeurs aberrantes et stockez les échantillons qui répondent aux exigences. Lorsque le résultat de la vérification est TRUE, passez à l’étape suivante. Enregistrez le résultat.
4. Utilisez la fonction PickSoftThreshold pour calculer l’indice sans échelle R² des deux groupes de données sous des valeurs de puissance différentes. Exécutez la ligne de commande :
   « sft = pickSoftThreshold(datExpr0, powerVector = puissances, verbe = 5) »
   Cliquez sur Exécuter.
5. Visualisez les résultats (Figure 1). Exécutez la ligne de commande :
   « plot(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2],
   xlab="Soft Threshold (power) »,ylab="Scale Free Topology Model Fit,signé R^2 »,type="n »,
   main = paste(« indépendance ES_Scale »));
   text(sft$fitIndices[,1], -sign(sft$fitIndices[,3])*sft$fitIndices[,2],
   labels=powers,cex=cex1,col="rouge »);
   abline(h=0,9,col="rouge »)
   plot(sft$fitIndices[,1], sft$fitIndices[,5],
   xlab="Soft Threshold (power) »,ylab="Mean Connectivity », type="n »,
   main = paste(« connectivité ES_Mean »))
   text(sft$fitIndices[,1], sft$fitIndices[,5], labels=powers, cex=cex1,col="red ») »
   Cliquez sur Exécuter.
   REMARQUE: La prémisse de l’algorithme de réseau de corrélation pondérée est que la structure de réseau de co-expression établie est conforme aux normes du critère de topologie sans échelle, ce qui augmente sa robustesse. Un index sans échelle plus proche de 1 indique une structure de réseau plus proche du réseau sans échelle.
6. Sélectionnez la valeur de puissance lorsque l’indice sans échelle R² est supérieur à 0,9 et passez à l’étape suivante de l’analyse.
   REMARQUE : Lorsque l’index sans échelle est proche de 1, la structure du réseau est plus proche du réseau sans échelle. Lors de l’analyse de deux ou plusieurs réseaux, il est nécessaire de choisir de rendre chaque réseau proche de la valeur de puissance du réseau sans échelle pour satisfaire la comparabilité entre les réseaux co-exprimés.

3. Construction d’un réseau de co-expression et identification de module

REMARQUE: Sur la base de la valeur de puissance calculée ci-dessus, le réseau de cooccurrence est construit. Les lignes de commande clés se réfèrent au Supplément S2.

1. Utilisez la fonction d’adjacence dans le package WGCNA pour ajouter des paramètres signés pour la construction d’un réseau de cooccurrence symbolique. Exécutez la ligne de commande :
   « adjacency = adjacency(datExpr0, power = softPower) »
   Cliquez sur Exécuter.
2. Appliquer la fonction de similitude TOM pour développer un réseau topologique qui se chevauche et calculer le réseau de dissimilarité. Exécutez la ligne de commande :
   « TOM = TOMsimilarité (contiguïté);
   dissTOM = 1-TOM »
   Cliquez sur Exécuter.
   Remarque : Le paramètre signé a été ajouté pour définir le type de réseau de chevauchement de topologie.
3. Utilisez la fonction hclust pour sélectionner la méthode de clustering hiérarchique de liaison moyenne pour le clustering hiérarchique. Exécutez la ligne de commande :
   « geneTree = hclust(as.dist(dissTOM), méthode = « moyenne »); »
   Cliquez sur Exécuter.
4. Utilisez la fonction cutreeDynamic pour effectuer une découpe dynamique des branches et définissez le paramètre minClusterSize sur 30. Obtenez le résultat de la reconnaissance du module. Exécutez la ligne de commande :
   « dynamicMods = cutreeDynamic(dendro = geneTree, distM = dissTOM, deepSplit = 2, pamRespectsDendro = FALSE, minClusterSize = minModuleSize); »
   Cliquez sur Exécuter.
   REMARQUE : La taille minimale du module ne peut pas être inférieure à 30.
5. Calculez l’eigen de module de chaque module OTUs par la fonction moduleEigengenes. Exécutez la ligne de commande :
   « MEList = moduleEigengenes(datExpr0, colors = dynamicColors)
   ME = MEList$eigengenes »
   Cliquez sur Exécuter.
   REMARQUE : Le module eigen représentait le niveau d’expression OTU global dans le module. Il ne s’agissait pas d’un OTU spécifique, mais de la première composante principale de chaque cluster obtenue par décomposition de valeur de réseau singulière.
6. Exécuter la fonction de cluster en fonction du coefficient de corrélation du module eigen. Utilisez la fonction mergeCloseModules pour fusionner les modules dont la valeur est inférieure à 0,25. Exécutez la ligne de commande :
   « merge = mergeCloseModules(datExpr0, dynamicColors, cutHeight = MEDissThres, verbose = 3) »
   Cliquez sur Exécuter.
7. Enfin, utilisez la fonction plotDendroAndColors pour la visualisation afin d’obtenir le diagramme d’affichage de l’affectation de module de chaque réseau de co-expression (Figure 2). Utilisez la fonction table pour extraire l’attribution de module correspondant à chaque OTin de la table d’affectation de module. Exécutez la ligne de commande :
   « plotDendroAndColors(geneTree, mergedColors, « Merged dynamic »,dendroLabels = FALSE,
   hang = 0.03,addGuide = TRUE, guideHang = 0.05,
   main = « ES_Gene dendrogramme et couleurs de module ») »
   Cliquez sur Exécuter.
   REMARQUE : Dans le diagramme d’affectation de module du réseau de co-expression, différentes couleurs représentent différents modules et le gris représente les OTU qui ne peuvent être classés dans aucun module. Un plus grand nombre d’OTU dans le module gris indique que la qualité de prétraitement précoce de la matrice d’expression est médiocre.

4. Comparaison des modules

REMARQUE: Cette méthode peut être utilisée pour comparer les modules de réseau de deux communautés microbiennes écologiques. Dans cet article, comparez les différences de modules de réseau microbien entre l’endosphère et le rhizoplan, l’endosphère et la rhizosphère, la rhizosphère et le rhizoplan.

1. Test de conservation
   1. Chargez les paramètres et les résultats des deux ensembles de données enregistrés dans les étapes précédentes.
   2. Définissez le résultat de l’affectation du module réseau d’un groupe de données microbiennes comme groupe de référence, tandis que l’autre groupe est le groupe de test.
   3. Utilisez la fonction modulePreservation pour calculer les valeurs des paramètres statistiques de conservateur Z_summary et medianRank. Exécutez la ligne de commande :
      « system.time({mp=modulePreservation(multiExpr,
      multiColor,referenceNetworks=1,
      nPermutation=100, randomSeed=1,quickCor=0,verbose=3)}) »
      Cliquez sur Exécuter.
      REMARQUE: Ce résultat peut quantifier la conservateurité entre les modules. Z_summary>10 indique que deux modules sont hautement conservés, tandis que Z_summary<2 désigne des modules non conservés. medianRank exprime la préservation relative du module évalué par classement. Des valeurs medianRank plus élevées indiquent des modules non conservés. (Les lignes de commande clés se réfèrent au Supplément S3.)
   4. Utilisez la fonction de tracé pour visualiser les résultats (Figure 3). Obtenez les paramètres Z_summary et medianRank (Tableau 1).
      REMARQUE: Les modules de réseau qui satisfont à la fois à la valeur Z_summary inférieure à 2 et à la valeur médiane Rank en haut, est le module le plus non préservé dans les deux communautés microbiennes écologiques.
   5. Basé sur les résultats des deux paramètres statistiques susmentionnés pour identifier le module avec le module le plus hautement non conservé des deux réseaux.
2. Analyse de corrélation de l’appartenance au module
   1. Définir les résultats d’affectation de module des deux réseaux ont été définis comme référence et comme groupe de test, respectivement.
      REMARQUE: Les paramètres doivent être les mêmes que le test de préservation.
   2. Utilisez la fonction corPvalueStudent pour extraire la valeur kME (appartenance au module) de chaque OTU dans plusieurs modules candidats.
      Exécutez la ligne de commande :
      « Pvalue = as.data.frame(corPvalueStudent(as.matrix
      (ModuleAdhésion), Échantillons)) »
      Cliquez sur Exécuter.
      REMARQUE: kME signifie le degré d’appartenance au module. ME signifie module eigen, qui représente le niveau global d’expression OTU dans le module. kME est le coefficient de corrélation entre chaque OTU et le ME. Quantifier l’importance de l’OTU dans le réseau par la valeur kME de l’OTU. (Les lignes de commande clés se réfèrent au supplément S4.)
   3. Ensuite, utilisez la fonction verboseScatterplot pour calculer le coefficient de corrélation de la valeur kME des OTU correspondants dans les deux réseaux et dessinez le diagramme d’analyse de corrélation (Figure 4).
      Exécutez la ligne de commande :
      « verboseScatterplot(abs(TModuleMembership
      [TmoduleGenes, Tcolumn]),
      abs(NModuleMembership[NmoduleGenes, Ncolumn]),
      xlab = paste(« kME in », « ES »),
      ylab = paste(« kME in », « RP »),
      main = paste(« lightyellow »),
      cex.main = 1.7, cex.lab = 1.6, cex.axis = 1.6, col = modulecolor) »
      Cliquez sur Exécuter.
   4. Sélectionnez le module avec le plus petit coefficient de corrélation de la valeur kME de l’OTU des deux réseaux. Considérez ce module comme ayant la plus grande différence des deux réseaux.

5. Analyse du module de réseau différentiel microbien

1. Obtenir des données sur le phyla de la bactérie dominante grâce à l’analyse statistique de l’ensemble de séquences OTU du module avec la plus grande différence.
   REMARQUE: L’ensemble de séquences OTU du module présentant la plus grande différence est additionné par la taxonomie des phyla. La bactérie dominante phyla représentait plus de 10%.
2. Ensuite, utilisez la fonction exportNetworkToCytoscape pour obtenir le fichier contenant les informations de relation d’interaction de l’OTU dans le plus grand module différentiel.
   Exécutez la ligne de commande :
   « cyt = exportNetworkToCytoscape(modTOM,
   edgeFile = paste(« NEW-ES_CytoscapeInput-edges-« , modules , « .txt », sep=" »),
   nodeFile = paste(« NEW-ES_CytoscapeInput-nodes-« , modules, « .txt », sep=" »),
   pondéré = TRUE,threshold = 0.5, nodeNames = modProbes,
   altNodeNames = modGenes, nodeAttr = moduleColors[inModule]) »
   Cliquez sur Exécuter.
3. Importez le fichier dans Cytoscape. Définissez le seuil sur 0,5 et ajustez d’autres paramètres si nécessaire.
4. Construire un réseau de cooccurrence de micro-organismes différentiels (Figure 5).
5. Obtention des informations sur le genre central qui a le rôle régulateur le plus important dans le réseau.
   REMARQUE: Selon la valeur kME de OUT, le genre de base peut être défini.
6. Enfin, les fonctions du genre central ont été évaluées et son influence sur l’ensemble du réseau de différences a été analysée.

Les résultats représentatifs de cet article ont été téléchargés à partir des données de 2014 sur le microbiome radiculaire du riz Abaker de Californie dans la base de données NCBI (PRJNA386367)⁹. Les données comprennent des échantillons de microbiome de la rhizosphère, du rhizoplan et de l’endosphère provenant de plants de riz cultivés pendant 14 semaines dans une rizière submergée. Nous avons utilisé l’algorithme WGCNA pour sélectionner la valeur de puissance qui satisfaisait les trois réseaux proches du réseau sans échelle (Figure 1) et développé trois réseaux de co-expression (Figure 2). Dans le réseau de co-expression microbienne du sol de l’endosphère, du rhizoplan et de la rhizosphère, 23, 22 et 21 modules ont été identifiés, respectivement. Ces résultats indiquent que le nombre de réseaux d’interaction microbienne dans les trois niches était fondamentalement égal.

Nous avons ensuite comparé les différences dans les modules de réseau microbien dans le sol de l’endosphère, du rhizoplan et de la rhizosphère. Les résultats suivants des tests de préservation des modules des trois groupes de niches ont été obtenus. Trois modules extrêmement non conservés existaient entre le sol de la rhizosphère et le rhizoplan(Figure 3a,Tableau 1). De plus, neuf modules extrêmement non conservés étaient présents entre le sol et l’endosphère de la rhizosphère(Figure 3b,Tableau 2)et six modules extrêmement non conservés entre le rhizoplan et l’endosphère (Figure 3c, Tableau 3). En outre, des modules extrêmement non conservés parmi les trois niches ont été trouvés, indiquant la présence de grandes différences dans la composition des micro-organismes entre les trois niches. Les résultats de l’analyse de corrélation des appartenances aux modules obtenus des modules non conservateurs sont illustrés à la figure 4. À partir de la figure, un module significativement différent avec le moins de corrélation entre chacune des deux niches est visible parmi les trois niches, ce qui représente la différence la plus significative entre elles.

Dans le réseau de différence rhizosphère-rhizoplan(figure 5a),le phylum dominant était celui des protéobactéries (72,97 %). Dans le réseau de différence rhizosphère-endosphère (Figure 5b), les phyla dominants étaient les protéobactéries (66,36%), les actinobactéries (10,1%) et les bactéroïdes (10,9%). Dans le réseau de différence rhizoplan-endosphère (Figure 5c), les phyla dominants étaient les protéobactéries (41,41%), les Bacteroidetes (10,10%), les Firmicutes (12,12%) et les Verrucomicrobia.

Trois genres principaux (Figure 5a), dont Rhodobacter et Novosphingobium, six genres de noyau (Figure 5b), dont Blvii28 et Dechloromonas, et cinq genres de noyau (Figure 5c), y compris Cellvibrio et Geobacter, ont exercé d’importantes fonctions de régulation dans les trois réseaux de cooccurrence différentielle. Tous les genres de base, à l’exception de Dechloromonas, n’ont eu d’influence que sur un seul réseau, indiquant la disponibilité de différences considérables dans l’abondance relative des populations et des espèces microbiennes entre les trois niches de racines de riz, ce qui a eu un effet critique sur l’abondance et la diversité des communautés microbiennes racinaires existantes.

Le genre de base Azospirillum,présent dans le réseau de différence rhizosphère-endosphère du riz, a participé à la fixation de l’azote et favorisé la croissance des plantes^10. De plus, le genre Geobacter, qui s’est considérablement enrichi dans le réseau de différence rhizoplan-endosphère, peut être le principal facteur induisant la réduction des oxydes insolubles de Fe et de Mn dans de nombreux sols et sédiments¹¹. Ces micro-organismes interagissent avec une gamme de communautés microbiennes dans les niches racinaires et participent activement à la régulation des réseaux microbiens, ce qui pourrait être d’une importance cruciale pour la croissance et le développement des racines de riz.

Graphique 1. Évaluation de powerβ dans les jeux de données. a) Évaluation des β de puissance dans l’ensemble de données ES; b) Évaluation de la puissance β dans l’ensemble de données RS; c) Évaluation de la puissance β dans l’ensemble de données RP, distribution de l’indice sans échelle R2 (à gauche), distribution de la connectivité moyenne (à droite) le long de différents indices de puissance douce. La valeur de la meilleure puissance a été obtenue lorsque R2 avait tendance à saturation et n’était pas inférieure à 0,8. Les trois réseaux de co-expression devaient être réglés sur la même valeur de puissance pour assurer leur comparabilité. (RS : sol de la rhizosphère, RP : rhizoplan et ES : endosphère). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Graphique 2. Dendrogramme OTU obtenu par regroupement hiérarchique de liaison moyenne. a)dendrogramme OTU du réseau ES; b) Dendrogramme OTU du réseau RS; c) Dendrogramme OTU du réseau RP. La ligne de couleur sous le dendrogramme indique l’affectation de module déterminée par l’algorithme Dynamic Tree Cut. (RS : sol de la rhizosphère, RP : rhizoplan et ES : endosphère) Veuillez cliquer ici pour voir une version agrandie de cette figure.

Graphique 3. Résultats du test de conservation. a) Les résultats de l’analyse sont fondés sur l’allocation du module réseau de co-expression RS en tant que groupe de référence et les résultats de l’allocation du module réseau de co-expression RP en tant que groupe d’essai; b) Les résultats de l’analyse sont fondés sur l’allocation du module réseau de co-expression ES en tant que groupe de référence et les résultats de l’allocation du module réseau de co-expression RS en tant que groupe d’essai; c) Les résultats de l’analyse sont fondés sur l’allocation du module réseau de co-expression ES en tant que groupe de référence et les résultats de l’allocation du module réseau de co-expression RP en tant que groupe d’essai. Z_summary > 10 indique que deux modules sont très conservés, tandis que Z_summary < 2 désigne des modules non conservés. medianRank exprime la préservation relative du module évalué par classement. Des valeurs medianRank plus élevées indiquent des modules non conservés. (RS : sol de la rhizosphère, RP : rhizoplan et ES : endosphère) Veuillez cliquer ici pour voir une version agrandie de cette figure.

Graphique 4. Analyse de corrélation des appartenances aux modules. a) Corrélation modulaire de la valeur kME entre le réseau RS et RP; b) Corrélation modulaire de la valeur kME entre le réseau RS et RP; c) Corrélation modulaire de la valeur kME entre rp et réseau ES (RS : sol de la rhizosphère, RP : rhizoplan, ES : endosphère et KME : appartenance au module). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Graphique 5. Réseau de cooccurrence de population microbienne différentielle dans la racine de riz. a)Réseau de cooccurrence de populations microbiennes différentielles dans la RP-RS; b) Réseau de cooccurrence de populations microbiennes différentielles dans l’ES-RS; c) Réseau de cooccurrence de populations microbiennes différentielles dans l’ES-RP. L’analyse a été réalisée à l’aide du logiciel Cytoscape. Différentes couleurs représentent différentes portes dans la figure. (RS: Sol de la rhizosphère, RP: rhizoplan et ES: endosphère). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Tableau 1. Résultat de Zsummary et medianRank entre le sol de la rhizosphère et le rhizoplan.

Tableau 2. Résultat de Zsummary et medianRank entre le sol de la rhizosphère et l’endosphère.

Tableau 3. Résultat de Zsummary et medianRank entre le rhizoplane et l’endosphère.

Supplément S1: Veuillez cliquer ici pour télécharger ce fichier.

Supplément S2 : Veuillez cliquer ici pour télécharger ce fichier.

Supplément S3: Veuillez cliquer ici pour télécharger ce fichier.

Supplément S4: Veuillez cliquer ici pour télécharger ce fichier.

Les réseaux de corrélation sont de plus en plus utilisés dans les applications bioinformatiques. WGCNA est une méthode de biologie des systèmes pour l’analyse descriptive des relations entre divers éléments d’un système biologique¹². Le progiciel R a été utilisé dans des travaux antérieurs sur WGCNA^13,14,15. Le package comprend des fonctions pour la construction de réseaux, la détection de modules, le calcul des propriétés topologiques, la simulation de données, la visualisation et la capacité d’interfaçage avec un logiciel externe. WGCNA a été largement utilisé pour analyser les données d’expression génique du cancer du cerveau¹⁶, du cycle cellulaire de levure¹⁷, de la génétique de la souris^18,19, du tissu cérébral primate^20,21, du diabète²²et des plantes²³. Utiliser l’analyse pondérée du réseau de corrélation génique pour construire le réseau doit inclure au moins 8 échantillons. Dans cet article, nous nous sommes concentrés sur les réseaux de co-expression génique qui décrivent les interactions entre les populations microbiennes dans différents environnements. Nous avons obtenu des réseaux différentiels entre les populations microbiennes dans des environnements divergents et identifié les espèces clés dans chaque réseau. L’idée que les espèces clés sont importantes pour la communauté a été largement utilisée dans la recherche sur le réseau^{trophique 24}. Certaines espèces d’une communauté microbienne complexe peuvent être essentielles pour maintenir la stabilité et la fonctionnalité de la communauté, comme les Bacteroides dans la flore intestinale²⁵. L’analyse des communautés microbiennes peut être considérablement simplifiée en ciblant des espèces spécifiques d’importance potentielle.

Nos résultats représentatifs mettent en évidence les différences dans les communautés microbiennes, qui peuvent être identifiées à l’aide de la méthode décrite ci-dessus. Ici, les micro-organismes dans différentes niches du système racinaire du riz ont été soumis à WGCNA. La différence entre les trois créneaux a été identifiée à l’aide d’analyses conservatrices et d’analyses de l’appartenance aux modules. Nous avons déterminé les espèces clés dans les modules de différence et obtenu des informations sur les différences dans la composition des communautés microbiennes dans les trois niches. Pendant ce temps, le réseau de cooccurrence a révélé la présence d’une interaction significative entre les micro-organismes changeants dans les racines du riz. Nos résultats fournissent des preuves directes de l’importance et de la faisabilité de WGCNA dans l’évaluation des différences de communauté microbienne dans divers environnements.

Les auteurs n’ont rien à divulguer.

Le développement de ce manuscrit a été soutenu par des fonds de la Fondation nationale des sciences naturelles du projet de centre de recherche sur les sciences karstiques du gouvernement populaire de Chine-Guizhou (U1812401), du projet de recherche doctorale de l’Université normale du Guizhou (GZNUD [2017]1), du projet de soutien scientifique et technologique de la province du Guizhou (QKHZC [2021] YB459) et du projet scientifique et technologique du Guiyang ([2019]2-8).

Les auteurs tiennent à remercier Edwards J.A et al pour avoir fourni des données sur le microbiome du riz dans des bases de données publiques et le soutien de TopEdit (www.topeditsci.com) pour son assistance linguistique lors de la préparation de ce manuscrit.