Analiza snRNA-seq profili transkryptomicznych gleju czołowego i adnotacja typów komórek
Łącznie uzyskano 220 095 jąder i 32 077 genów z kory czołowej 17 samców i 17 samic z AD (Rycyna 1A). Wykres UMAP zwizualizował całkowite transkryptomy pojedynczych jąder czołowych, wykazując odrębne typy jąder po analizie redukcji wymiarowości (Rycina 1B). Przedstawiono całkowitą liczbę adnotowanych jąder z podziałem na płeć, co dało sumę obejmującą 58 902 astrocytów, 14 265 mikroglejów, 77 466 oligodendrocytów, 3 520 komórek śródbłonka, 25 252 neuronów pobudzających, 31 268 neuronów hamujących oraz 9 422 komórek progenitorowych oligodendrocytów (Rycina 1C). Średnia ekspresja znanych markerów typów komórek dla każdego rodzaju gleju została rzutowana na wykresy UMAP w celu identyfikacji populacji komórkowych (Rycina 1D).
Specyficzne dla płci geny różnicowo wyrażone (DEG) w astrocytach
Przeanalizowano 58 902 jądra astrocytów (Rysunek 2A), z czego 27 504 pochodziło od samców z AD (46,69%), a 31 398 od samic z AD (53,31%). Specyficzne dla płci DEG ujawniły nadekspresję 138 genów, w tym DST, CACNA2D3 oraz AC016831.7, niedoekspresję 105 genów, w tym CNTN5, RORA, RASSF8 oraz CADM2, a także 4995 genów bez zmian w ekspresji (Rysunek 2B). Ponadto analizy GO i KEGG wykazały, że DEG te koncentrują się głównie w szlakach dotyczących neuronów, synaps i hormonów; szlaki neuronalne obejmowały m.in. regulację rozwoju wypustek neuronowych i kolców dendrytycznych, szlaki synaptyczne obejmowały organizację synaptyczną oraz synapsy glutaminergiczne/cholinergiczne, a szlaki hormonalne obejmowały syntezę hormonów tarczycy i sekrecję insuliny (Rysunek 2C,D). Wyłoniono trzydzieści genów o najwyższej częstotliwości występowania (Rysunek 2E), wśród których na pierwszym miejscu znalazł się PLCB1. W celu zbadania relacji między tymi genami skonstruowano sieci PPI (Rysunek 2F), które pozwoliły zidentyfikować DLG2, CAMK2D, CALM2 oraz PRKACB jako geny kluczowe. Wykresy UMAP wykazały różnice w ekspresji wybranych DEG: KCND2, CAMK2D, MT3, LINC00278 oraz XIST wykazywały wyższy poziom u samic z AD i niższy u samców z AD, natomiast w przypadku NLGN4Y, DLG2, PRKACB, CALM2, UTY oraz TTTY14 zależności były odwrotne (Rysunek 2G). Ostatecznie, jak pokazano w Tabeli 1, jako specyficzne dla płci DEG astrocytów określono geny MT3, CALM2, DLG2, KCND2, PRKACB, CAMK2D oraz NLGN4Y.
Swoiste dla płci DEGs w mikrogleju
Przeanalizowano 14 265 jąder mikrogleju (Rycyna 3A), z czego 5 327 (37,34%) pochodziło od samców z AD, a 8 938 (62,66%) od samic z AD. Analiza swoistych dla płci DEGs wykazała nadekspresję 224 genów, w tym KCNIP4 i LRRTM4, niedoekspresję 930 genów, w tym APOE, MT-CO3 i FTL, oraz brak zmian w 13 111 genach (Rycyna 3B). Ponadto analizy GO i KEGG wykazały, że DEGs te koncentrowały się głównie w szlakach neuronów, fagosomów, hormonów i innych; szlaki neuronalne obejmowały synapsy neuron-neuron, szlaki fagocytarne obejmowały fagosom i regulację fagocytozy, szlaki hormonalne obejmowały szlak sygnalizacyjny estrogenów, szlak sygnalizacyjny oksytocyny itp., a pozostałe obejmowały regulację odpowiedzi zapalnej, uczenie się lub pamięć, usuwanie amyloidu-beta itp. (Rycyna 3C,D). Wyłoniono trzydzieści genów o najwyższej częstotliwości, przy czym TLR2 i TREM2 zajęły ex aequo drugie miejsce (Rycyna 3E). W celu zbadania relacji między tymi genami skonstruowano sieci PPI (Rycyna 3F), które pozwoliły zidentyfikować ACTB, APP i FYN jako geny kluczowe. Wykresy UMAP wykazały różnice w ekspresji wybranych DEGs: APP, FOS, XIST i CTSD miały wyższy poziom u samic z AD i niższy u samców z AD, natomiast w przypadku NLGN4Y, TREM2, LINC0028, APOE, UTY i TTTY14 obserwowano sytuację odwrotną (Rycyna 3G). Ostatecznie jako swoiste dla płci DEGs mikrogleju określono TREM2, FOS, APOE, APP i NLGN4Y, co przedstawiono w Tabeli 2.
Seksowo specyficzne DEGs w oligodendrocytach
Przeanalizowano 77 466 jąder oligodendrocytów (Rysunek 4A), z czego 42 469 pochodziło od mężczyzn z AD (54,82%), a 34 997 od kobiet z AD (45,18%). Analiza seksowo specyficznych DEGs ujawniła nadekspresję 384 genów, w tym PCDH9, MT-CO1, NEAT1 i NPAS3, niedoekspresję 188 genów, w tym FRMD4A, PLP1 i LSAMP, natomiast pozostałe 76 894 geny nie uległy zmianie (Rysunek 4B). Ponadto analizy GO i KEGG wykazały, że DEGs te koncentrowały się głównie w szlakach dotyczących neuronów, synaps i hormonów; szlaki neuronalne obejmowały kolce dendrytyczne, szlaki synaptyczne obejmowały synapsy neuron-neuron oraz synapsy glutaminergiczne/dopaminergiczne, a szlaki hormonalne obejmowały m.in. szlak sygnalizacyjny neurotrofin, syntezę i sekrecję aldosteronu oraz szlak sygnalizacji wapniowej (Rysunek 4C,D). Wyłoniono trzydzieści genów o najwyższej częstotliwości występowania (Rysunek 4E), przy czym GRIN2A i PSEN1 zajęły ex aequo drugie miejsce. Konstrukcje sieci PPI posłużyły do zbadania relacji między tymi genami (Rysunek 4F), co pozwoliło zidentyfikować GRIN2A i GRIA2 jako geny kluczowe. Wykresy UMAP wykazały różnice w wybranych DEGs: GRIN2A, ITPR2, GNAS i XIST charakteryzowały się wyższym poziomem u kobiet z AD i niższym u mężczyzn z AD, natomiast w przypadku NLGN4Y, UTY i TTTY14 obserwowano zjawisko przeciwne (Rysunek 4G). Ostatecznie jako seksowo specyficzne DEGs oligodendrocytów określono GRIN2A, ITPR2, GNAS oraz NLGN4Y, co przedstawiono w Tabeli 3.
Interakcje między kluczowymi genami a 30 najważniejszymi genami komórek glejowych oraz NLGN4Y jako wspólny gen
Diagramy Venna oraz sieci PPI zapewniły przegląd bliskich interakcji między kluczowymi genami (Rysunek 5A,B) a 30 najważniejszymi genami (Rysunek 5C, D) każdej komórki glejowej. Wyniki wykazały, że ACYB, APP, JUN, PRKACB oraz DLG2 znajdowały się w centrum sieci PPI, a NLGN4Y był wspólnym genem dla specyficznych dla płci DEG we wszystkich komórkach glejowych.
Konstrukcja sieci gen-TF-miRNA
Sieć NLGN4Y-TF-miRNA składała się z 13 węzłów i 12 krawędzi (Rycyna 6A). Ekspresja NLGN4Y była regulowana przez 1 TF, a mianowicie CTCF, oraz 11 miRNA, w tym hsa-miR-185, hsa-miR-137 i hsa-miR-9.
Wyświetlanie leku docelowego i TCM dla NLGN4Y wraz z siecią
W bazie Coremine Medical zidentyfikowano łącznie 1 lek docelowy dla NLGN4Y oraz 64 pośrednie docelowe TCM. TCM wykazujące istotność statystyczną w wynikach zostały zaznaczone na niebiesko. Zwizualizowano lek Antithrombin III oraz pięć TCM wchodzących w skład sieci, mianowicie Heikunbu, Wulingzhi, Xiazhicao, Shuizhi i Mahuang, które uznano za istotne statystycznie (Rysunek 6B).
Wpływ docelowych TCM i odpowiadających im składników aktywnych na AD
W przypadku Kunbu w bazie TCMIP wyszukano 10 składników, a w TCMSP 48 składników. Poprzez przeszukiwanie bazy danych PubMed zidentyfikowano 5 składników związanych z AD: fukosterol, saringosterol, tiaminę, kwas stearydonowy oraz phlorofucofuroeckol-A, z których dwa pierwsze charakteryzowały się biodostępnością doustną (OB) ≥30% i podobieństwem do leku (DL) ≥0,18. Szczegóły przedstawiono w Tabeli 4. W odniesieniu do Mahuang wyszukano 28 składników w TCMIP i 363 składniki w TCMSP. Stosując tę samą metodę co w przypadku Kunbu, zidentyfikowano 25 składników związanych z AD. Spośród nich 6 składników z OB ≥30% i DL ≥0,18: kwercetyna, eriodyktyol, naringenina, taksifolina, stigmasterol i luteolina, wymieniono w górnej części Tabeli 5.
Analiza wzbogacenia i sieć hub-gen-choroba
Sieć hub-gen-choroba została przedstawiona na Rysunku 6C. Wspólnym genem docelowym dla Kunbu i Mahuang był ACHE, z którym związane choroby to AD oraz deficyty poznawcze. Genami docelowymi wyłącznie dla Kunbu były ALKBH3 i ELOVL4, z którymi związane choroby to odpowiednio starzenie się i zanik mózgu.
Wyniki wzbogaceń GO (BP, MF i CC) dla Kunbu (Rysunek 7A, Rysunek 8A, Rysunek 9A) i Mahuang (Rysunek 7B, Rysunek 8B, Rysunek 9B) przedstawiono na wykresach słupkowych. Pokazano również ścieżki Reactome dla Kunbu (Rysunek 10A) i Mahuang (Rysunek 10B). Wzbogacone terminy zaznaczone strzałkami były związane z osią „hormon-synapsa-neuron”, takie jak aktywność receptora hormonów steroidowych, projekcja neuronów i szlak sygnalizacyjny zależny od hormonów steroidowych w przypadku Kunbu oraz chemiczna transmisja synaptyczna, estrogenowo-zależna ekspresja genów i procesy układu nerwowego w przypadku Mahuang.

Rycina 1: Pozyskiwanie danych o ekspresji genów, profilowanie RNA-seq pojedynczych jąder oraz charakterystyka typów komórek. (A) Próbki pobrane z GEO DataSets w celu przygotowania analizy. (B) Dwuwymiarowy wykres UMAP sumy jąder (N = 116,101 dla samców; N = 103,994 dla samic). (C) Proporcje poszczególnych typów komórek z podziałem na płeć. (D) Średnia ekspresja 5 dobrze ustanowionych markerów typów komórek rzutowana na wykres UMAP. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rysunek 2: Astrocyty są heterogeniczne i wykazują specyficzne dla płci zmiany transkrypcyjne w chorobie Alzheimera. (A) Wykres UMAP jąder astrocytów (N = 59 010). (B) Specyficzne dla płci DEGs. (C,D) Wykresy kołowe ilustrujące istotne funkcjonalnie wzbogacone terminy dla astrocytarnych DEGs uzyskane z baz danych GO i KEGG (GO: C, KEGG: D). (E) 30 najczęściej występujących astrocytarnych DEGs specyficznych dla płci w ścieżkach GO i KEGG. (F) Sieć PPI dla 30 najczęściej występujących astrocytarnych DEGs specyficznych dla płci w ścieżkach GO i KEGG. (G) Średnia ekspresja istotnych, specyficznych dla płci DEGs rzutowana na wykresy UMAP. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Mikroglej jest heterogenny i wykazuje specyficzne dla płci zmiany transkrypcyjne w chorobie Alzheimera. (A) Wykres UMAP jąder mikrogleju (N = 14 265). (B) Specyficzne dla płci DEGs. (C,D) Wykresy kołowe ilustrujące istotne funkcjonalnie wzbogacone terminy mikroglejowych DEGs uzyskane z baz danych GO i KEGG (GO: C, KEGG: D). (E) 30 najczęstszych mikroglejowych DEGs specyficznych dla płci w szlakach GO i KEGG. (F) Sieć PPI 30 najczęstszych mikroglejowych DEGs specyficznych dla płci w szlakach GO i KEGG. (G) Średnia ekspresja istotnych specyficznych dla płci DEGs rzutowana na wykresy UMAP. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 4: Oligodendrocyty są heterogeniczne i wykazują specyficzne dla płci zmiany transkrypcyjne w chorobie Alzheimera. (A) Wykres UMAP jąder oligodendrocytów (N = 77 466). (B) Specyficzne dla płci geny różnicowo wyrażone (DEG). (C,D) Wykresy kołowe ilustrujące istotne funkcjonalnie wzbogacone terminy dla DEG oligodendrocytów uzyskane z baz danych GO i KEGG (GO: C, KEGG: D). (E) 30 najczęstszych DEG oligodendrocytów specyficznych dla płci w szlakach GO i KEGG. (F) Sieć interakcji białko-białko (PPI) dla 30 najczęstszych DEG oligodendrocytów specyficznych dla płci w szlakach GO i KEGG. (G) Średnia ekspresja wybranych specyficznych dla płci DEG rzutowana na wykresy UMAP. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 5: Części wspólne DEG oraz sieci PPI dla najczęściej występujących DEG. (A,B) Diagram Venna ilustrujący kluczowe geny dla każdego typu gleju oraz odpowiadająca im sieć PPI. (C,D) Diagram Venna ilustrujący 30 najczęstszych płciowo specyficznych DEG dla każdego typu gleju, wzbogaconych w ścieżki GO i KEGG, oraz odpowiadająca im sieć PPI. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 6: Schemat sieci z NLGN4Y i celowanymi lekami TCM jako rdzeniem. (A) Sieć koregulacyjna gen-TF-miRNA. (B) Sieć gen-lek-TCM dla NLGN4Y. (C) Sieć TCM-gen-choroba dla Kunbu i Mahuang. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 7: Analiza wzbogacenia GO (BP) dla Kunbu i Mahuang. (A) 20 najlepiej wzbogaconych procesów biologicznych dla Kunbu. (B) 20 najlepiej wzbogaconych procesów biologicznych dla Mahuang. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 8: Analiza wzbogacenia GO (MF) dla Kunbu i Mahuang. (A) 20 najlepiej wzbogaconych funkcji molekularnych Kunbu. (B) 20 najlepiej wzbogaconych funkcji molekularnych Mahuang. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 9: Analiza wzbogacenia GO (CC) dla Kunbu i Mahuang. (A) 20 najlepiej wzbogaconych komponentów komórkowych Kunbu. (B) 20 najlepiej wzbogaconych komponentów komórkowych Mahuang. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 10: Szlaki Reactome dla Kunbu i Mahuang. (A) Szlaki Reactome dla Kunbu. (B) Szlaki Reactome dla Mahuang. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Tabela 1: Kluczowe geny specyficznych dla płci różnicowo wyrażonych genów (DEG) w astrocytach. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Kluczowe geny płciowo specyficznych DEGs w mikrogleju. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 3: Kluczowe geny wśród specyficznych dla płci DEGs w oligodendrocytach. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 4: Aktywne składniki Kunbu w chorobie AD. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 5: Aktywne składniki Mahuang w chorobie AD. Kliknij tutaj, aby pobrać tę tabelę.
Rycina uzupełniająca 1: Zrzut ekranu z użycia oprogramowania R. Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 1: Informacje o próbkach. Kliknij tutaj, aby pobrać ten plik.
Tabela uzupełniająca 2: Geny markerowe. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 1: KOD R. Kliknij tutaj, aby pobrać ten plik.