SnRNA-seq-Analyse von frontalen glialen Transkriptomprofilen und die Annotation von Zelltypen
Insgesamt wurden 220.095 Kerne und 32.077 Gene im frontalen Kortex von 17 männlichen und 17 weiblichen AD gewonnen (Abbildung 1A). Das UMAP-Diagramm visualisierte die Gesamtzahl der frontalen Einzelkerntranskriptome, die nach der Dimensionsreduktionsanalyse unterschiedliche Kerntypen zeigten (Abbildung 1B). Es wurde eine Gesamtanzahl der annotierten Zellkerne gezeigt, die nach Geschlecht erfasst wurden, was die Summe von 58.902 Astrozyten, 14.265 Mikroglia, 77.466 Oligodendrozyten, 3.520 endothelialen, 25.252 exzitatorischen Neuronen, 31.268 inhibitorischen Neuronen und 9.422 Oligodendrozyten-Vorläuferzellen ergab (Abbildung 1C). Die durchschnittliche Expression der bekannten Zelltypmarker für jede Gliazelle wurde auf die UMAP-Diagramme projiziert, um die Zellpopulationen zu identifizieren (Abbildung 1D).
Die geschlechtsspezifischen DEGs in Astrozyten
Es wurden 58.902 astrozytäre Kerne analysiert (Abbildung 2A), davon 27.504 von männlichen AD (46,69%) und 31.398 von weiblichen AD (53,31%). Geschlechtsspezifische DEGs zeigten eine Hochregulierung von 138 Genen, einschließlich DST, CACNA2D3 und AC016831,7, eine Herunterregulierung von 105 Genen, einschließlich CNTN5, RORA, RASSF8 und CADM2, und unveränderte 4995 Gene (Abbildung 2B). Darüber hinaus zeigten GO- und KEGG-Analysen, dass sich diese DEGs hauptsächlich auf die Signalwege von Neuronen, Synapsen und Hormonen konzentrierten, mit neuronalen Signalwegen einschließlich der Regulation der Neuronenprojektentwicklung, der Neuronenwirbelsäule usw., synaptischen Signalwegen einschließlich der synaptischen Organisation und der glutamatergen/cholinergen Synapse und Hormonwegen einschließlich der Schilddrüsenhormonsynthese und der Insulinsekretion (Abbildung 2C, D). Es wurden dreißig Gene mit der höchsten Häufigkeit erhalten (Abbildung 2E), wobei PLCB1 an erster Stelle steht. PPI-Netzwerke wurden konstruiert, um die Beziehung zwischen diesen Genen zu untersuchen (Abbildung 2F) und DLG2, CAMK2D, CALM2 und PRKACB als Kerngene identifiziert. UMAP-Diagramme zeigten Unterschiede zwischen ausgewählten DEGs: KCND2, CAMK2D, MT3, LINC00278 und XIST waren bei weiblichen AD höher und bei männlichen AD niedriger, während NLGN4Y, DLG2, PRKACB, CALM2, UTY und TTTY14 entgegengesetzt waren (Abbildung 2G). MT3, CALM2, DLG2, KCND2, PRKACB, CAMK2D und NLGN4Y wurden schließlich als geschlechtsspezifische DEGs von Astrozyten bestimmt, wie in Tabelle 1 gezeigt.
Die geschlechtsspezifischen DEGs in Mikroglia
14.265 Mikrogliakerne (Abbildung 3A) wurden analysiert, davon 5.327 (37,34%) von männlichen AD und 8.938 (62,66%) von weiblichen AD. Geschlechtsspezifische DEGs zeigten eine Hochregulation von 224 Genen, einschließlich KCNIP4 und LRRTM4, und eine Herunterregulierung von 930 Genen, einschließlich APOE, MT-CO3 und FTL, und unveränderte 13.111 Gene (Abbildung 3B). Darüber hinaus zeigten GO- und KEGG-Analysen, dass sich diese DEGs hauptsächlich auf die Signalwege von Neuronen, Phagosomen, Hormonen und anderen konzentrierten, mit neuronalen Signalwegen einschließlich Neuron-zu-Neuron-Synapse, phagozytären Signalwegen einschließlich Phagosomen und Regulation der Phagozytose, Hormonwegen einschließlich Östrogensignalweg, Oxytocin-Signalweg usw. und anderen einschließlich der Regulation von Entzündungsreaktion, Lernen oder Gedächtnis, Amyloid-Beta-Clearance, usw. (Abbildung 3C,D). Es wurden dreißig Gene mit der höchsten Häufigkeit erhalten, wobei TLR2 und TREM2 den zweiten Platz belegten (Abbildung 3E). PPI-Netzwerke wurden konstruiert, um die Beziehung zwischen diesen Genen zu untersuchen (Abbildung 3F) und identifizierten ACTB, APP und FYN als Kerngene. UMAP-Diagramme zeigten Unterschiede zwischen ausgewählten DEGs: APP, FOS, XIST und CTSD waren bei weiblichen AD höher und bei männlichen AD niedriger, während NLGN4Y, TREM2, LINC0028, APOE, UTY und TTTY14 entgegengesetzt waren (Abbildung 3G). TREM2, FOS, APOE, APP und NLGN4Y wurden schließlich als geschlechtsspezifische DEGs von Mikroglia bestimmt, wie in Tabelle 2 gezeigt.
Die geschlechtsspezifischen DEGs in Oligodendrozyten
Es wurden 77.466 oligodendrozytäre Kerne analysiert (Abbildung 4A), davon 42.469 von männlichen AD (54,82%) und 34.997 von weiblichen AD (45,18%). Geschlechtsspezifische DEGs zeigten eine Hochregulierung von 384 Genen, darunter PCDH9, MT-CO1, NEAT1 und NPAS3, eine Herunterregulierung von 188 Genen, einschließlich FRMD4A, PLP1 und LSAMP, und die restlichen 76.894 Gene unverändert (Abbildung 4B). Darüber hinaus zeigten GO- und KEGG-Analysen, dass sich diese DEGs hauptsächlich auf die Signalwege von Neuronen, Synapsen und Hormonen konzentrierten, mit neuronalen Signalwegen einschließlich Neuronenwirbelsäule, synaptischen Signalwegen einschließlich Neuron-zu-Neuron-Synapse und glutamaterger/dopaminerger Synapse und Hormonwegen einschließlich Neurotrophin-Signalweg, Aldosteronsynthese und -sekretion, Kalzium-Signalweg usw. (Abbildung 4C, D). Es wurden dreißig Gene mit der höchsten Häufigkeit erhalten (Abbildung 4E), wobei GRIN2A und PSEN1 den zweiten Platz belegten. PPI-Netzwerke wurden konstruiert, um die Beziehung zwischen diesen Genen zu untersuchen (Abbildung 4F) und GRIN2A und GRIA2 als Kerngene identifiziert. UMAP-Diagramme zeigten Unterschiede in ausgewählten DEGs: GRIN2A, ITPR2, GNAS und XIST waren bei weiblichen AD höher und bei männlichen AD niedriger, während NLGN4Y, UTY und TTTY14 entgegengesetzt waren (Abbildung 4G). GRIN2A, ITPR2, GNAS und NLGN4Y wurden schließlich als geschlechtsspezifische DEGs von Oligodendrozyten bestimmt, wie in Tabelle 3 gezeigt.
Interaktion zwischen Schlüsselgenen und den 30 wichtigsten Genen von Gliazellen und NLGN4Y als gemeinsames Gen
Die Venn-Diagramme und die PPI-Netzwerke gaben einen Überblick über die enge Interaktion zwischen Schlüsselgenen (Abbildung 5A, B) und den 30 wichtigsten Genen (Abbildung 5C, D) jeder Gliazelle. Die Ergebnisse zeigten, dass ACYB, APP, JUN, PRKACB und DLG2 den Kern des PPI-Netzwerks bildeten und NLGN4Y ein gemeinsames Gen für geschlechtsspezifische DEGs in allen Gliazellen war.
Aufbau eines Gen-TF-miRNA-Netzwerks
Das NLGN4Y-TF-miRNA-Netzwerk enthielt 13 Knoten und 12 Kanten (Abbildung 6A). NLGN4Y wurde durch 1 TF, nämlich CTCF, und 11 miRNAs, einschließlich has-miR-185, has-miR-137 und has-miR-9, reguliert.
Anzeige des Zielmedikaments und der TCMs von NLGN4Y mit Netzwerk
Insgesamt wurden 1 NLGN4Y-Zielmedikament und 64 indirekte Ziel-TCMs in Coremine Medical gewonnen. TCMs mit statistischer Signifikanz in den Ergebnissen wurden blau markiert. Das Medikament Antithrombin III und fünf TCMs mit dem Netzwerk, nämlich Heikunbu, Wulingzhi, Xiazhicao, Shuizhi und Mahuang, wurden visualisiert, die als statistisch signifikant angesehen wurden (Abbildung 6B).
Auswirkungen von Ziel-TCMs und entsprechenden Wirkstoffen auf AD
Für Kunbu wurden 10 Inhaltsstoffe in TCMIP und 48 Inhaltsstoffe in TCMSP gewonnen. Durch die Suche in der PubMed-Datenbank wurden 5 Inhaltsstoffe im Zusammenhang mit AD gefunden: Fucosterol, Saringosterol, Thiamin, Stearidonsäure und Phlorofucofuroeckol-A, von denen die orale Bioverfügbarkeit (OB) der ersten beiden ≥30% und ähnlich wie das Medikament Sex (DL) ≥0,18 betrug. Die Details sind in Tabelle 4 dargestellt. Was Mahuang betrifft, so wurden 28 Inhaltsstoffe in TCMIP und 363 Inhaltsstoffe in TCMSP gefunden. Mit der gleichen Methode wie bei Kunbu wurden 25 Inhaltsstoffe gewonnen, die mit Alzheimer in Verbindung stehen. Unter ihnen sind 6 Inhaltsstoffe mit OB ≥30% und DL ≥0,18: Quercetin, Eriodictyol, Naringenin, Taxifolin, Stigmasterol und Luteolin, oben in Tabelle 5 aufgeführt.
Hurb-Gen-Krankheitsnetzwerk und Anreicherungsanalyse
Das Hurb-Gen-Krankheitsnetzwerk ist in Abbildung 6C dargestellt. Das gemeinsame Zielgen von Kunbu und Mahuang war ACHE, dessen assoziierte Krankheiten Alzheimer und kognitive Defizite waren. Die Zielgene von Kunbu allein waren ALKBH3 und ELOVL4, Krankheiten, die mit dem Altern bzw. der Hirnatrophie in Verbindung stehen.
Die Ergebnisse der GO-Anreicherungen (BP, MF und CC) für Kunbu (Abbildung 7A, Abbildung 8A, Abbildung 9A) und Mahuang (Abbildung 7B, Abbildung 8B, Abbildung 9B) werden in Balkendiagrammen angezeigt. Die Reaktomwege für Kunbu (Abbildung 10A) und Mahuang (Abbildung 10B) sind ebenfalls dargestellt. Die angereicherten Begriffe, die mit Pfeilen markiert waren, standen in Beziehung zur "Hormon-Synapse-Neuron"-Achse, wie z.B. der Steroidhormonrezeptoraktivität, der Neuronenprojektion und dem Steroidhormon-vermittelten Signalweg in Kunbu und der chemischen synaptischen Übertragung, der östrogenabhängigen Genexpression und dem Prozess des Nervensystems in Mahuang.

Abbildung 1: Erfassung von Genexpressionsdaten, Erstellung von RNA-Seq-Profilen mit Einzelkernen und Charakterisierung des Zelltyps. (A) Proben, die aus GEO-DataSets für die Analysevorbereitung gewonnen wurden. (B) 2-dimensionales UMAP-Diagramm der Summenkerne (N = 116.101 für männlich; N = 103.994 für Frauen). (C) Anteile für jeden Zelltyp, aufgeteilt nach Geschlecht. (D) Durchschnittliche Expression von 5 gut etablierten Zelltypmarkern, die auf das UMAP-Diagramm projiziert wurden. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 2: Astrozyten sind heterogen und weisen geschlechtsspezifische transkriptomische Veränderungen bei der Alzheimer-Krankheit auf. (A) UMAP-Diagramm der astrozytären Kerne (N = 59.010). (B) Geschlechtsspezifische DEGs. (C,D) Die Kreisdiagramme veranschaulichten die signifikanten funktionell angereicherten Terme astrozytärer DEGs, die aus GO- und KEGG-Datenbanken (GO: C, KEGG: D) gewonnen wurden. (E) Die Top 30 der häufigsten Astrozyten-Geschlechts-DEGs in GO- und KEGG-Signalwegen. (F) PPI-Netzwerk der Top 30 der häufigsten Astrozyten-Geschlechts-DEGs in GO- und KEGG-Signalwegen. (G) Durchschnittliche Expression bemerkenswerter geschlechtsspezifischer DEGs, projiziert auf die UMAP-Diagramme. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 3: Mikroglia sind heterogen und weisen geschlechtsspezifische transkriptomische Veränderungen bei der Alzheimer-Krankheit auf. (A) UMAP-Diagramm von Mikrogliakernen (N = 14.265). (B) Geschlechtsspezifische DEGs. (C,D) Kreisdiagramme veranschaulichten die signifikanten funktionell angereicherten Terme von Mikroglia-DEGs, die aus GO- und KEGG-Datenbanken (GO: C, KEGG: D) gewonnen wurden. (E) Die Top 30 der häufigsten Mikroglia-Geschlechts-DEGs inGO- und KEGG-Signalwegen. (F) PPI-Netzwerk der Top 30 der häufigsten Mikroglia-Geschlechts-DEGs in GO- und KEGG-Signalwegen. (G) Durchschnittliche Expression bemerkenswerter geschlechtsspezifischer DEGs, projiziert auf die UMAP-Diagramme. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 4: Oligodendrozyten sind heterogen und weisen geschlechtsspezifische transkriptomische Veränderungen bei der Alzheimer-Krankheit auf. (A) UMAP-Diagramm der oligodendrozytären Kerne (N = 77.466). (B) Geschlechtsspezifische DEGs. (C,D) Die Kreisdiagramme veranschaulichten die signifikanten funktionell angereicherten Terme oligodendrozytärer DEGs, die aus GO- und KEGG-Datenbanken (GO: C, KEGG: D) gewonnen wurden. (E) Die 30 häufigsten oligodendrozyten-geschlechtlichen DEGs in GO- und KEGG-Signalwegen. (F) PPI-Netzwerk der 30 häufigsten Oligodendrozyten-Geschlechts-DEGs in GO- und KEGG-Signalwegen. (G) Durchschnittliche Expression bemerkenswerter geschlechtsspezifischer DEGs, projiziert auf die UMAP-Diagramme. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 5: Überlappungen von DEGs und PPI-Netzwerken von hochfrequenten DEGs. (A,B) Venn-Diagramm, das die Schlüsselgene für jede Glia und das entsprechende PPI-Netzwerk veranschaulicht. (C,D) Venn-Diagramm, das die 30 häufigsten geschlechtsspezifischen DEGs für jede Glia, die in GO- und KEGG-Signalwegen angereichert ist, und das entsprechende PPI-Netzwerk veranschaulicht. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 6: Netzwerkdiagramm mit NLGN4Y und seinen Targeting-TCMs als Kern. (A) Gen-TF-miRNA-koregulatorisches Netzwerk. (B) Gen-Wirkstoff-TCM-Netzwerk für NLGN4Y. (C) TCM-Gen-Krankheitsnetzwerk von Kunbu und Mahuang. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 7: GO (BP) Anreicherungsanalyse für Kunbu und Mahuang. (A) Top 20 angereicherte biologische Prozesse von Kunbu. (B) Top 20 angereicherte biologische Prozesse von Mahuang. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 8: GO (MF) Anreicherungsanalyse für Kunbu und Mahuang. (A) Top 20 angereicherte molekulare Funktionen von Kunbu. (B) Top 20 angereicherte molekulare Funktionen von Mahuang. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 9: GO (CC) Anreicherungsanalyse für Kunbu und Mahuang. (A) Top 20 angereicherte zelluläre Komponenten von Kunbu. (B) Top 20 angereicherte zelluläre Bestandteile von Mahuang. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 10: Reaktomwege für Kunbu und Mahuang. (A) Reaktomwege von Kunbu. (B) Reaktomwege von Mahuang. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Tabelle 1: Schlüsselgene der geschlechtsspezifischen DEGs in Astrozyten. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 2: Schlüsselgene der geschlechtsspezifischen DEGs in Mikroglia. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 3: Schlüsselgene der geschlechtsspezifischen DEGs in Oligodendrozyten. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 4: Kunbu'-Wirkstoffe bei AD. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Tabelle 5: Mahuang'-Wirkstoffe bei AD. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.
Ergänzende Abbildung 1: Screenshot für die Nutzung der R-Software. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 1: Informationen zum Beispiel. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 2: Marker-Gene. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 1: Der R-CODE. Bitte klicken Sie hier, um diese Datei herunterzuladen.