SnRNA-seq analyse van frontale gliale transcriptoomprofielen en de annotatie van celtypen
In totaal werden 220.095 kernen en 32.077 genen in de frontale cortex van 17 mannelijke AD en 17 vrouwelijke AD verkregen (Figuur 1A). UMAP-plot visualiseerde de totale frontale transcriptomen met één kern die verschillende soorten kernen vertoonden na analyse van dimensiereductie (Figuur 1B). Het totale aantal geannoteerde kernen dat per geslacht was vastgelegd, werd getoond, wat de som maakte van 58.902 astrocyten, 14.265 microglia, 77.466 oligodendrocyten, 3.520 endotheel, 25.252 exciterende neuronen, 31.268 remmende neuronen en 9.422 oligodendrocyt-voorlopercellen (Figuur 1C). Gemiddelde expressies van de bekende celtypemarkers voor elke glia werden geprojecteerd op de UMAP-plots om de celpopulaties te identificeren (Figuur 1D).
De geslachtsspecifieke DEG's in astrocyten
58.902 astrocytaire kernen werden geanalyseerd (Figuur 2A), waarvan 27.504 van mannelijke AD (46,69%) en 31.398 van vrouwelijke AD (53,31%). Geslachtsspecifieke DEG's onthulden opregulatie van 138 genen, waaronder DST, CACNA2D3 en AC016831.7, neerwaartse regulatie van 105 genen, waaronder CNTN5, RORA, RASSF8 en CADM2, en onveranderde 4995 genen (Figuur 2B). Verder toonden GO- en KEGG-analyses aan dat die DEG's zich voornamelijk concentreerden in de routes van neuronen, synapsen en hormonen, met neuronale routes waaronder regulatie van de ontwikkeling van neuronenprojecten, neuronwervelkolom, enz., synaptische routes waaronder synaptische organisatie en glutamaterge/cholinerge synaps, en hormoonroutes waaronder schildklierhormoonsynthese en insulinesecretie (Figuur 2C,D). Dertig genen met een topfrequentie werden verkregen (Figuur 2E), waarbij PLCB1 de eerste was. PPI-netwerken werden geconstrueerd om de relatie tussen deze genen te onderzoeken (Figuur 2F) en identificeerden DLG2, CAMK2D, CALM2 en PRKACB als kerngenen. UMAP-plots vertoonden verschillen van geselecteerde DEG's: KCND2, CAMK2D, MT3, LINC00278 en XIST waren hoger bij vrouwelijke AD en lager bij mannelijke AD, terwijl NLGN4Y, DLG2, PRKACB, CALM2, UTY en TTTY14 tegenovergesteld waren (Figuur 2G). MT3, CALM2, DLG2, KCND2, PRKACB, CAMK2D en NLGN4Y werden uiteindelijk bepaald als geslachtsspecifieke DEG's van astrocyten, zoals weergegeven in tabel 1.
De geslachtsspecifieke DEG's bij microglia
14.265 microgliale kernen (Figuur 3A) werden geanalyseerd, waarvan 5.327 (37,34%) van mannelijke AD en 8.938 (62,66%) van vrouwelijke AD. Geslachtsspecifieke DEG's onthulden opregulatie van 224 genen, waaronder KCNIP4 en LRRTM4, en neerwaartse regulatie van 930 genen, waaronder APOE, MT-CO3 en FTL, en onveranderde 13,111 genen (Figuur 3B). Verder toonden GO- en KEGG-analyses aan dat die DEG's zich voornamelijk concentreerden in de routes van neuronen, fagosoom, hormonen en andere, met neuronale routes waaronder neuron-naar-neuronsynaps, fagocytische routes waaronder fagosoom en regulatie van fagocytose, hormoonroutes waaronder oestrogeensignaleringsroute, oxytocine-signaleringsroute, enz., en andere, waaronder regulatie van ontstekingsreactie, leren of geheugen, amyloïde-bèta-klaring, enz. (Figuur 3C,D). Dertig genen met een topfrequentie werden verkregen, waarbij TLR2 en TREM2 op de tweede plaats eindigden (Figuur 3E). PPI-netwerken werden geconstrueerd om de relatie tussen deze genen te onderzoeken (Figuur 3F) en identificeerden ACTB, APP en FYN als kerngenen. UMAP-plots vertoonden verschillen van geselecteerde DEG's: APP, FOS, XIST en CTSD waren hoger bij vrouwelijke AD en lager bij mannelijke AD, terwijl NLGN4Y, TREM2, LINC0028, APOE, UTY en TTTY14 tegengesteld waren (Figuur 3G). TREM2, FOS, APOE, APP en NLGN4Y werden uiteindelijk bepaald als geslachtsspecifieke DEG's van microglia, zoals weergegeven in Tabel 2.
De geslachtsspecifieke DEG's in oligodendrocyten
77.466 oligodendrocytische kernen werden geanalyseerd (Figuur 4A), waarvan 42.469 van mannelijke AD (54,82%) en 34.997 van vrouwelijke AD (45,18%). Geslachtsspecifieke DEG's toonden opregulatie van 384 genen, waaronder PCDH9, MT-CO1, NEAT1 en NPAS3, neerwaartse regulatie van 188 genen, waaronder FRMD4A, PLP1 en LSAMP, en de overige 76.894 genen onveranderd (Figuur 4B). Verder bleek uit GO- en KEGG-analyses dat die DEG's zich voornamelijk concentreerden in de routes van neuronen, synapsen en hormonen, met neuronale routes waaronder neuronwervelkolom, synaptische routes waaronder neuron-naar-neuronsynaps en glutamaterge/dopaminerge synaps, en hormoonroutes waaronder neurotrofine-signaleringsroute, aldosteronsynthese en -secretie, calciumsignaleringsroute, enz. (Figuur 4C, D). Dertig genen met een topfrequentie werden verkregen (Figuur 4E), waarbij GRIN2A en PSEN1 op de tweede plaats stonden. PPI-netwerken werden geconstrueerd om de relatie tussen deze genen te onderzoeken (Figuur 4F) en identificeerden GRIN2A en GRIA2 als kerngenen. UMAP-plots vertoonden verschillen in geselecteerde DEG's: GRIN2A, ITPR2, GNAS en XIST waren hoger in AD bij vrouwen en lager bij AD bij mannen, terwijl NLGN4Y, UTY en TTTY14 tegengesteld waren (Figuur 4G). GRIN2A, ITPR2, GNAS en NLGN4Y werden uiteindelijk bepaald als geslachtsspecifieke DEG's van oligodendrocyten, zoals weergegeven in Tabel 3.
Interactie tussen sleutelgenen en de top 30 genen van gliacellen, en NLGN4Y als het gemeenschappelijke gedeelde gen
De Venn-diagrammen en de PPI-netwerken gaven een overzicht van de nauwe interactie tussen sleutelgenen (Figuur 5A,B) en de top 30 genen (Figuur 5C, D) van elke gliacel. De resultaten toonden aan dat ACYB, APP, JUN, PRKACB en DLG2 de kern van het PPI-netwerk vormden en dat NLGN4Y een gemeenschappelijk gedeeld gen was voor geslachtsspecifieke DEG's in alle gliacellen.
Constructie van het Gene-TF-miRNA-netwerk
Het NLGN4Y-TF-miRNA-netwerk bevatte 13 knooppunten en 12 randen (Figuur 6A). NLGN4Y werd gereguleerd door 1 TF, namelijk CTCF, en 11 miRNA's, waaronder has-miR-185, has-miR-137 en has-miR-9.
Weergave van doelgeneesmiddel en TCM's van NLGN4Y met netwerk
In totaal werden 1 NLGN4Y-doelgeneesmiddel en 64 indirecte doel-TCM's opgehaald in Coremine Medical. TCM's met statistische significantie in de resultaten werden gemarkeerd met blauw. Geneesmiddel Antitrombine III en vijf TCM's met het netwerk, namelijk Heikunbu, Wulingzhi, Xiazhicao, Shuizhi en Mahuang, werden gevisualiseerd, die als statistisch significant werden beschouwd (Figuur 6B).
Effecten van doel-TCM's en overeenkomstige actieve ingrediënten op AD
Voor Kunbu werden 10 ingrediënten in TCMIP en 48 ingrediënten in TCMSP opgehaald. Door de PubMed-database te doorzoeken, werden 5 ingrediënten gevonden die verband houden met AD: fucosterol, saringosterol, thiamine, stearidonzuur en phlorofucofuroeckol-A, waarvan de orale biologische beschikbaarheid (OB) van de eerste twee ≥30% was en vergelijkbaar met het medicijngeslacht (DL) ≥ 0,18 was. Details worden weergegeven in tabel 4. Wat Mahuang betreft, werden 28 ingrediënten in TCMIP en 363 ingrediënten in TCMSP opgehaald. Met dezelfde methode als Kunbu werden 25 ingrediënten met betrekking tot AD opgehaald. Onder hen staan 6 ingrediënten met OB ≥30% en DL ≥0,18: quercetine, eriodictyol, naringenine, taxifoline, stigmasterol en luteoline bovenaan tabel 5.
Hurb-gen-ziektenetwerk en verrijkingsanalyse
Het Hurb-gen-ziektenetwerk wordt weergegeven in figuur 6C. Het gedeelde doelgen van Kunbu en Mahuang was ACHE, waarvan de geassocieerde ziekten AD en cognitieve stoornissen waren. De doelgenen van Kunbu alleen waren ALKBH3 en ELOVL4, ziekten die verband hielden met respectievelijk veroudering en hersenatrofie.
De resultaten van GO-verrijkingen (BP, MF en CC) voor Kunbu (Figuur 7A, Figuur 8A, Figuur 9A) en Mahuang (Figuur 7B, Figuur 8B, Figuur 9B) worden weergegeven met staafdiagrammen. Reactoomroutes voor Kunbu (Figuur 10A) en Mahuang (Figuur 10B) worden ook getoond. De verrijkte termen gemarkeerd met pijlen waren gerelateerd aan de "hormoon-synaps-neuron" -as, zoals steroïde hormoonreceptoractiviteit, neuronprojectie en steroïde hormoon-gemedieerde signaalroute in Kunbu en chemische synaptische transmissie, oestrogeenafhankelijke genexpressie en zenuwstelselproces in Mahuang.

Figuur 1: Verkrijgen van genexpressiegegevens, single-nuclei RNA-seq-profilering en karakterisering van celtypen. (A) Monsters verkregen uit GEO DataSets voor analysevoorbereiding. (B) 2-dimensionale UMAP-grafiek van somkernen (N = 116.101 voor mannen; N = 103.994 voor vrouwen). (C) Verhoudingen voor elk type cel, uitgesplitst naar geslacht. (D) Gemiddelde expressie van 5 gevestigde celtypemarkers geprojecteerd op de UMAP-plot. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: Astrocyten zijn heterogeen en hebben geslachtsspecifieke transcriptomische veranderingen bij de ziekte van Alzheimer. (A) UMAP-plot van astrocytaire kernen (N = 59.010). (B) Geslachtsspecifieke DEG's. (C,D) Cirkeldiagrammen illustreerden de significante functioneel verrijkte termen van astrocytaire DEG's verkregen uit GO- en KEGG-databases (GO: C, KEGG: D). (E) De top 30 van meest voorkomende astrocyt-geslacht DEG's in GO- en KEGG-routes. (F) PPI-netwerk van de top 30 van meest voorkomende astrocyt-geslacht DEG's in GO- en KEGG-routes. (G) Gemiddelde expressie van opmerkelijke geslachtsspecifieke DEG's geprojecteerd op de UMAP-plots. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Microglia zijn heterogeen en hebben geslachtsspecifieke transcriptomische veranderingen bij de ziekte van Alzheimer. (A) UMAP-diagram van microgliakernen (N = 14.265). (B) Geslachtsspecifieke DEG's. (C,D) Cirkeldiagrammen illustreerden de significante functioneel verrijkte termen van microgliale DEG's verkregen uit GO- en KEGG-databases (GO: C, KEGG: D). (E) De top 30 van meest voorkomende microglia-geslacht DEG's in GO en KEGGpathways. (F) PPI-netwerk van de top 30 van meest voorkomende microglia-geslacht DEG's in GO- en KEGG-routes. (G) Gemiddelde expressie van opmerkelijke geslachtsspecifieke DEG's geprojecteerd op de UMAP-plots. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: Oligodendrocyten zijn heterogeen en vertonen geslachtsspecifieke transcriptomische veranderingen bij de ziekte van Alzheimer. (A) UMAP-diagram van oligodendrocytische kernen (N = 77.466). (B) Geslachtsspecifieke DEG's. (C,D) Cirkeldiagrammen illustreerden de significante functioneel verrijkte termen van oligodendrocytische DEG's verkregen uit GO- en KEGG-databases (GO: C, KEGG: D). (E) De top 30 van meest voorkomende oligodendrocyten-geslacht DEG's in GO- en KEGG-routes. (F) PPI-netwerk van de top 30 van meest voorkomende oligodendrocyt-geslacht DEG's in GO- en KEGG-routes. (G) Gemiddelde expressie van opmerkelijke geslachtsspecifieke DEG's geprojecteerd op de UMAP-plots. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5: Overlappingen van DEG's en PPI-netwerken van DEG's met topfrequentie. (A,B) Venndiagram ter illustratie van de sleutelgenen voor elke glia en het bijbehorende PPI-netwerk. (C,D) Venndiagram ter illustratie van de top 30 frequente geslachtsspecifieke DEG's voor elke glia verrijkt in GO- en KEGG-routes en het bijbehorende PPI-netwerk. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: Netwerkdiagram met NLGN4Y en de TCM's als kern. (A) Gen-TF-miRNA coregulerend netwerk. (B) Gen-geneesmiddel-TCM netwerk voor NLGN4Y. (C) TCM-gen-ziekte netwerk van Kunbu en Mahuang. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7: GO (BP) Verrijkingsanalyse voor Kunbu en Mahuang. (A) Top 20 verrijkte biologische processen van Kunbu. (B) Top 20 verrijkte biologische processen van Mahuang. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 8: GO (MF) Verrijkingsanalyse voor Kunbu en Mahuang. (A) Top 20 verrijkte moleculaire functies van Kunbu. (B) Top 20 verrijkte moleculaire functies van Mahuang. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 9: GO (CC) Verrijkingsanalyse voor Kunbu en Mahuang. (A) Top 20 verrijkte cellulaire componenten van Kunbu. (B) Top 20 verrijkte cellulaire componenten van Mahuang. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 10: Reactoom routes voor Kunbu en Mahuang. (A) Reactoom routes van Kunbu. (B) Reactoom-routes van Mahuang. Klik hier om een grotere versie van deze figuur te bekijken.
Tabel 1: Sleutelgenen van de geslachtsspecifieke DEG's in astrocyten. Klik hier om deze tabel te downloaden.
Tabel 2: Sleutelgenen van de geslachtsspecifieke DEG's in microglia. Klik hier om deze tabel te downloaden.
Tabel 3: Sleutelgenen van de geslachtsspecifieke DEG's in oligodendrocyten. Klik hier om deze tabel te downloaden.
Tabel 4: Kunbu' actieve ingrediënten op AD. Klik hier om deze tabel te downloaden.
Tabel 5: Actieve ingrediënten van Mahuang op AD. Klik hier om deze tabel te downloaden.
Aanvullende afbeelding 1: Screenshot voor het gebruik van de R-software. Klik hier om dit bestand te downloaden.
Aanvullende tabel 1: Voorbeeldinformatie. Klik hier om dit bestand te downloaden.
Aanvullende tabel 2: Markergenen. Klik hier om dit bestand te downloaden.
Aanvullend bestand 1: De R-CODE. Klik hier om dit bestand te downloaden.