Analyse SnRNA-seq des profils de transcriptome glial frontal et annotation des types cellulaires
Au total, 220 095 noyaux et 32 077 gènes dans le cortex frontal de 17 MA mâles et 17 DA femelles ont été obtenus (Figure 1A). Le graphique UMAP a visualisé l’ensemble des transcriptomes frontaux des noyaux uniques affichant des types distincts de noyaux après analyse de réduction de dimension (Figure 1B). Le nombre total de noyaux annotés capturés par sexe a été indiqué, ce qui fait la somme de 58 902 astrocytes, 14 265 microglies, 77 466 oligodendrocytes, 3 520 endothéliaux, 25 252 neurones excitateurs, 31 268 neurones inhibiteurs et 9 422 cellules progénitrices oligodendrocytes (Figure 1C). Les expressions moyennes des marqueurs de type cellulaire connus pour chaque glie ont été projetées sur les graphiques UMAP afin d’identifier les populations cellulaires (Figure 1D).
Les DEG spécifiques au sexe dans les astrocytes
58 902 noyaux astrocytaires ont été analysés (Figure 2A), dont 27 504 chez les mâles (46,69%) et 31 398 chez les femelles (53,31%). Les DEG spécifiques au sexe ont révélé une régulation à la hausse de 138 gènes, y compris l’heure DST, CACNA2D3 et AC016831,7, une régulation à la baisse de 105 gènes, y compris CNTN5, RORA, RASSF8 et CADM2, et 4995 gènes inchangés (Figure 2B). De plus, les analyses GO et KEGG ont révélé que ces DEG se concentraient principalement dans les voies des neurones, des synapses et des hormones, avec des voies neuronales comprenant la régulation du développement du projet neuronal, de la colonne vertébrale neuronale, etc., des voies synaptiques, y compris l’organisation synaptique et la synapse glutamatergique/cholinergique, et des voies hormonales comprenant la synthèse des hormones thyroïdiennes et la sécrétion d’insuline (Figure 2C, D). Trente gènes avec la fréquence maximale ont été obtenus (Figure 2E), PLCB1 se classant en premier. Des réseaux d’IPP ont été construits pour explorer la relation entre ces gènes (figure 2F) et ont identifié DLG2, CAMK2D, CALM2 et PRKACB comme gènes de base. Les graphiques UMAP ont montré des différences entre les DEG sélectionnés : KCND2, CAMK2D, MT3, LINC00278 et XIST étaient plus élevés chez les femelles et plus faibles chez les mâles, tandis que NLGN4Y, DLG2, PRKACB, CALM2, UTY et TTTY14 étaient opposés (Figure 2G). MT3, CALM2, DLG2, KCND2, PRKACB, CAMK2D et NLGN4Y ont finalement été déterminés comme des DEG spécifiques au sexe des astrocytes, comme le montre le tableau 1.
Les DEG spécifiques au sexe dans la microglie
14 265 noyaux microgliaux (figure 3A) ont été analysés, dont 5 327 (37,34 %) chez les mâles et 8 938 (62,66 %) chez les femelles. Les DEG spécifiques au sexe ont révélé une régulation à la hausse de 224 gènes, y compris KCNIP4 et LRRTM4, et une régulation à la baisse de 930 gènes, dont APOE, MT-CO3 et FTL, et 13 111 gènes inchangés (figure 3B). De plus, les analyses GO et KEGG ont identifié que ces DEG se concentraient principalement dans les voies des neurones, des phagosomes, des hormones et autres, avec des voies neuronales comprenant la synapse de neurone à neurone, les voies phagocytaires y compris le phagosome et la régulation de la phagocytose, les voies hormonales y compris la voie de signalisation des œstrogènes, la voie de signalisation de l’ocytocine, etc., et d’autres y compris la régulation de la réponse inflammatoire, de l’apprentissage ou de la mémoire, la clairance de la bêta-amyloïde, etc. (Figure 3C,D). Trente gènes avec la fréquence maximale ont été obtenus, TLR2 et TREM2 étant à égalité pour la deuxième place (Figure 3E). Des réseaux d’IPP ont été construits pour explorer la relation entre ces gènes (figure 3F) et ont identifié ACTB, APP et FYN comme gènes de base. Les graphiques UMAP ont montré des différences entre les DEG sélectionnés : APP, FOS, XIST et CTSD étaient plus élevés chez les femmes et plus faibles chez les mâles, tandis que NLGN4Y, TREM2, LINC0028, APOE, UTY et TTTY14 étaient opposés (figure 3G). TREM2, FOS, APOE, APP et NLGN4Y ont finalement été déterminés comme des DEG de la microglie spécifiques au sexe, comme le montre le tableau 2.
Les DEG spécifiques au sexe dans les oligodendrocytes
77 466 noyaux oligodendrocytaires ont été analysés (figure 4A), dont 42 469 chez les mâles (54,82 %) et 34 997 chez les femelles (45,18 %). Les DEG spécifiques au sexe ont révélé une régulation à la hausse de 384 gènes, y compris PCDH9, MT-CO1, NEAT1 et NPAS3, une régulation négative de 188 gènes, y compris FRMD4A, PLP1 et LSAMP, et les 76 894 gènes restants inchangés (figure 4B). De plus, les analyses GO et KEGG ont révélé que ces DEG se concentraient principalement dans les voies des neurones, des synapses et des hormones, avec des voies neuronales comprenant la colonne vertébrale des neurones, des voies synaptiques comprenant la synapse de neurone à neurone et la synapse glutamatergique/dopaminergique, et les voies hormonales, y compris la voie de signalisation de la neurotrophine, la synthèse et la sécrétion d’aldostérone, la voie de signalisation du calcium, etc. Trente gènes avec la fréquence maximale ont été obtenus (Figure 4E), GRIN2A et PSEN1 étant à égalité en deuxième position. Des réseaux d’IPP ont été construits pour explorer la relation entre ces gènes (figure 4F) et ont identifié GRIN2A et GRIA2 comme gènes de base. Les graphiques UMAP ont montré des différences dans certaines DEG : GRIN2A, ITPR2, GNAS et XIST étaient plus élevés chez les femmes et plus faibles chez les mâles, tandis que NLGN4Y, UTY et TTTY14 étaient opposés (figure 4G). GRIN2A, ITPR2, GNAS et NLGN4Y ont finalement été déterminés en tant que DEG spécifiques au sexe des oligodendrocytes, comme le montre le tableau 3.
Interaction entre les gènes clés et les 30 principaux gènes des cellules gliales, et NLGN4Y en tant que gène commun partagé
Les diagrammes de Venn et les réseaux PPI ont fourni une vue d’ensemble de l’interaction étroite entre les gènes clés (figures 5A, B) et les 30 premiers gènes (figures 5C, D) de chaque cellule gliale. Les résultats ont montré que ACYB, APP, JUN, PRKACB et DLG2 étaient au cœur du réseau IPP, et que NLGN4Y était un gène commun partagé pour les DEG spécifiques au sexe dans toutes les cellules gliales.
Construction d’un réseau gène-TF-miARN
Le réseau de miARN NLGN4Y-TF contenait 13 nœuds et 12 arêtes (Figure 6A). NLGN4Y a été régulé par 1 TF, à savoir CTCF, et 11 miARN, dont has-miR-185, has-miR-137 et has-miR-9.
Affichage du médicament cible et des MTC de NLGN4Y avec réseau
Au total, 1 médicament cible NLGN4Y et 64 MTC cibles indirectes ont été récupérés dans Coremine Medical. Les MTC ayant une signification statistique dans les résultats ont été marqués en bleu. Le médicament Antithrombine III et cinq MTC du réseau, à savoir Heikunbu, Wulingzhi, Xiazhicao, Shuizhi et Mahuang, ont été visualisés, ce qui a été considéré comme statistiquement significatif (figure 6B).
Effets des MTC cibles et des principes actifs correspondants sur la DA
Pour Kunbu, 10 ingrédients du TCMIP et 48 ingrédients du TCMSP ont été récupérés. En effectuant une recherche dans la base de données PubMed, 5 ingrédients liés à la MA ont été trouvés : le fucostérol, le saringostérol, la thiamine, l’acide stéaridonique et le phlorofucofuroeckol-A, dont la biodisponibilité orale (OB) des deux premiers était de ≥30 % et similaire au sexe du médicament (DL) était de ≥0,18. Le tableau 4 donne plus de détails. En ce qui concerne Mahuang, 28 ingrédients du TCMIP et 363 ingrédients du TCMSP ont été récupérés. Avec la même méthode que Kunbu, 25 ingrédients liés à la DA ont été récupérés. Parmi eux, 6 ingrédients avec OB ≥30% et DL ≥0,18 : quercétine, ériodictyol, naringénine, taxifoline, stigmastérol et lutéoline, sont énumérés en haut du tableau 5.
Analyse du réseau et de l’enrichissement des gènes de Hurb
Le réseau Hurb-gène-maladie est illustré à la figure 6C. Le gène cible partagé par Kunbu et Mahuang était ACHE, dont les maladies associées étaient la MA et les déficits cognitifs. Les gènes cibles de Kunbu seul étaient ALKBH3 et ELOVL4, des maladies liées au vieillissement et à l’atrophie cérébrale, respectivement.
Les résultats des enrichissements en GO (BP, MF et CC) pour Kunbu (figure 7A, figure 8A, figure 9A) et Mahuang (figure 7B, figure 8B, figure 9B) sont affichés à l’aide de diagrammes à barres. Les trajectoires du réactome pour Kunbu (figure 10A) et Mahuang (figure 10B) sont également illustrées. Les termes enrichis marqués de flèches étaient liés à l’axe « hormone-synapse-neurone », tels que l’activité des récepteurs des hormones stéroïdiennes, la projection des neurones et la voie de signalisation médiée par les hormones stéroïdes à Kunbu et la transmission synaptique chimique, l’expression génique dépendante des œstrogènes et le processus du système nerveux à Mahuang.

Figure 1 : Acquisition de données d’expression génique, profilage du séquençage de l’ARN à noyaux uniques et caractérisation du type cellulaire. (A) Échantillons obtenus à partir de GEO DataSets pour la préparation de l’analyse. (B) Graphique UMAP en 2 dimensions de la somme des noyaux (N = 116 101 pour les mâles ; N = 103 994 pour les femmes). (C) Proportions pour chaque type de cellule réparties par sexe. (D) Expression moyenne de 5 marqueurs de type cellulaire bien établis projetés sur le graphique UMAP. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 2 : Les astrocytes sont hétérogènes et présentent des modifications transcriptomiques spécifiques au sexe dans la maladie d’Alzheimer. (A) Graphique UMAP des noyaux astrocytaires (N = 59 010). (B) DEG spécifiques au sexe. (C, D) Les diagrammes circulaires ont illustré les termes significativement enrichis fonctionnellement des DEG astrocytaires obtenus à partir des bases de données GO et KEGG (GO : C, KEGG : D). (E) Les 30 DEG astrocytes-sexe les plus fréquents dans les voies GO et KEGG. (F) Réseau d’IPP des 30 DEG astrocytaires-sexe les plus fréquents dans les voies GO et KEGG. (G) Expression moyenne de DEG remarquables spécifiques au sexe projetées sur les graphiques UMAP. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Les microglies sont hétérogènes et présentent des modifications transcriptomiques spécifiques au sexe dans la maladie d’Alzheimer. (A) Graphique UMAP des noyaux microgliaux (N = 14 265). (B) DEG spécifiques au sexe. (C, D) Des diagrammes circulaires ont illustré les termes significativement enrichis fonctionnellement des DEG microgliaux obtenus à partir des bases de données GO et KEGG (GO : C, KEGG : D). (E) Les 30 DEG microglia-sexe les plus fréquents inGO et KEGGpathways. (F) Réseau d’IPP des 30 DEG microglia-sexe les plus fréquents dans les voies GO et KEGG. (G) Expression moyenne de DEG remarquables spécifiques au sexe projetées sur les graphiques UMAP. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Les oligodendrocytes sont hétérogènes et présentent des modifications transcriptomiques spécifiques au sexe dans la maladie d’Alzheimer. (A) Graphique UMAP des noyaux oligodendrocytaires (N = 77 466). (B) DEG spécifiques au sexe. (C, D) Des diagrammes circulaires ont illustré les termes significativement enrichis fonctionnellement des DEG oligodendrocytaires obtenus à partir des bases de données GO et KEGG (GO : C, KEGG : D). (E) Les 30 DEG les plus fréquentes des oligodendrocytes et du sexe dans les voies GO et KEGG. (F) Réseau IPP des 30 DEG oligodendrocytaires-sexe les plus fréquents dans les voies GO et KEGG. (G) Expression moyenne de DEG remarquables spécifiques au sexe projetées sur les graphiques UMAP. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Chevauchements des DEG et des réseaux d’IPP des DEG de fréquence supérieure. (A,B) Diagramme de Venn illustrant les gènes clés de chaque glie et le réseau d’IPP correspondant. (C, D) Diagramme de Venn illustrant les 30 principaux DEG spécifiques au sexe fréquents pour chaque glie enrichie dans les voies GO et KEGG et le réseau IPP correspondant. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : Schéma de réseau avec NLGN4Y et ses MTC ciblant comme noyau. (A) Réseau de corégulation gène-TF-miARN. (B) Réseau gène-médicament-MTC pour NLGN4Y. (C) Réseau MTC-gène-maladie de Kunbu et Mahuang. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 7 : Analyse de l’enrichissement de GO (BP) pour Kunbu et Mahuang. (A) Les 20 principaux processus biologiques enrichis de Kunbu. (B) Les 20 principaux processus biologiques enrichis de Mahuang. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 8 : Analyse d’enrichissement GO (MF) pour Kunbu et Mahuang. (A) Les 20 principales fonctions moléculaires enrichies de Kunbu. (B) Les 20 principales fonctions moléculaires enrichies de Mahuang. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 9 : Analyse de l’enrichissement GO (CC) pour Kunbu et Mahuang. (A) Les 20 principaux composants cellulaires enrichis de Kunbu. (B) Les 20 principaux composants cellulaires enrichis de Mahuang. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 10 : Trajectoires du réacatome de Kunbu et de Mahuang. (A) Trajectoires du réacatome de Kunbu. (B) Voies du réacatome de Mahuang. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Tableau 1 : Gènes clés des DEG spécifiques au sexe dans les astrocytes. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 2 : Gènes clés des DEG spécifiques au sexe dans la microglie. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 3 : Gènes clés des DEG spécifiques au sexe dans les oligodendrocytes. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 4 : Principes actifs de Kunbu' sur la maladie d’Alzheimer. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 5 : Principes actifs de Mahuang' sur la maladie d’Alzheimer. Veuillez cliquer ici pour télécharger ce tableau.
Figure supplémentaire 1 : Capture d’écran de l’utilisation du logiciel R. Veuillez cliquer ici pour télécharger ce fichier.
Tableau supplémentaire 1 : Renseignements sur l’échantillon. Veuillez cliquer ici pour télécharger ce fichier.
Tableau supplémentaire 2 : Gènes marqueurs. Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 1 : Le code R. Veuillez cliquer ici pour télécharger ce fichier.