De belangrijkste bevindingen van de voorgestelde studie worden hier gepresenteerd, waaronder verschillende analyses die zijn uitgevoerd om de moleculaire mechanismen die ten grondslag liggen aan NAFLD en MI op te helderen.
Identificatie van DEG's
In de GSE89632-dataset werden 76 opgereguleerde en 20 neerwaarts gereguleerde genen geïdentificeerd als NAFLD-DEG's (Figuur 2B,D), terwijl de GSE66360-dataset 118 opgereguleerde en 8 neerwaarts gereguleerde genen onthulde als MI-DEG's (Figuur 2C,E). Vervolgens werden co-DEG's afgeleid uit de overlap tussen NAFLD-DEG's en MI-DEG's, waarbij 13 opwaarts gereguleerde DEG's en 1 neerwaarts gereguleerde DEG onder de co-DEG's werden geïdentificeerd (Figuur 2A).
Het identificeren van DEG's die uniek zijn voor NAFLD en MI, evenals die welke beide aandoeningen gemeen hebben, biedt waardevolle inzichten in de moleculaire processen die ten grondslag liggen aan deze ziekten. Inzicht in de verschillende en gedeelde genexpressiepatronen die verband houden met beide aandoeningen helpt de pathofysiologie van NAFLD en MI te verduidelijken en kan potentiële moleculaire doelen voor therapeutische interventie identificeren. Verder onderzoek naar de co-DEG's zou kritische regulatoren of routes kunnen onthullen die essentieel zijn voor de interactie tussen deze twee aandoeningen, waardoor de ontwikkeling van gerichte therapieën wordt vergemakkelijkt die gericht zijn op het verminderen van de negatieve gezondheidsresultaten die verband houden met NAFLD en MI.
Verrijking analyse
GO-analyse toonde aan dat NAFLD DEG's voornamelijk werden verrijkt in de volgende biologische processen (BP): vetceldifferentiatie, respons op lipopolysaccharide en meercellige organismeprocessen (p < 0,05, figuur 3A). Voor cellulaire componenten (CC) waren NAFLD DEG's voornamelijk geassocieerd met het fosfatidylinositol 3-kinasecomplex, secretoir granulelumen en cytoplasmatisch blaasjeslumen (p < 0,05, figuur 3A). In termen van moleculaire functie (MF) omvatten de belangrijkste activiteiten fosfatidylinositol-3-kinaseregulatoractiviteit, DNA-bindende transcriptieactivatoractiviteit (RNA-polymerasespecifiek) en algemene DNA-bindende transcriptieactivatoractiviteit (p < 0,05, Figuur 3A).
Voor MI-DEG's toonde de BP-analyse verrijking in leukocytenchemotaxis, activering van myeloïde leukocyten en celchemotaxis (p < 0,05, figuur 3B). De CC-analyse benadrukte verrijking in tertiaire korrels, ficoline-1-rijke korrels en secretoir korrellumen (p < 0,05, Figuur 3B). MF-analyse onthulde verrijking in patroonherkenningsreceptoractiviteit, RAGE-receptorbinding en immuunreceptoractiviteit (p < 0,05, Figuur 3B).
KEGG-verrijkingsanalyse gaf aan dat NAFLD-DEG's significant verrijkt waren in de volgende routes: IL-17, AGE-RAGE, TNF, osteoclastdifferentiatie, malaria, reumatoïde artritis, ziekte van Chagas en JAK-STAT-signaleringsroutes (p < 0,05, figuur 3C). MI-DEG's waren verrijkt in IL-17, TNF, lipidenmetabolisme en atherosclerose, toll-like receptorsignalering, C-type lectinereceptorsignalering, legionellose, reumatoïde artritis en NF-KB-signaleringsroutes (p < 0,05, Figuur 3D). Disease Ontology (DO) verrijkingsanalyse toonde aan dat de top drie DEG's voor zowel NAFLD als MI significant geassocieerd waren met atherosclerose, arteriosclerotische hart- en vaatziekten en arteriosclerose (p < 0,05, figuur 3E-F).
Analyse van PPI-netwerkdiagrammen
In het PPI-netwerkdiagram voor NAFLD werden 92 knooppunten en 250 randen geïdentificeerd, wat een weerspiegeling is van een complex netwerk van interacties tussen eiwitten die verband houden met deze aandoening. Evenzo onthulde het PPI-netwerkdiagram voor MI 117 knooppunten en 587 randen, wat wijst op een ingewikkelder netwerk dat kenmerkend is voor het moleculaire landschap van een hartinfarct. Bovendien identificeerden de PPI-netwerkdiagrammen voor co-DEG's 14 knooppunten en 21 randen, wat de nadruk legde op mogelijke sleutelinteracties tussen eiwitten die worden gedeeld tussen NAFLD en MI (Figuur 4A-F).
PPI-netwerken zijn waardevol voor het construeren en analyseren van therapeutische doelen en bieden belangrijke inzichten in de moleculaire verbindingen die ten grondslag liggen aan NAFLD en MI. De unieke netwerktopologieën die in NAFLD en MI worden waargenomen, benadrukken de verschillende pathofysiologische mechanismen die bij deze ziekten betrokken zijn. Bovendien suggereert de ontdekking van gedeelde verbindingen tussen co-DEG's potentiële gedeelde routes of regulerende netwerken die de interactie tussen MI en NAFLD bemiddelen. Verder onderzoek naar deze verbanden zou nieuwe moleculaire doelen voor interventie en therapieontwikkeling aan het licht kunnen brengen, die mogelijk nieuwe behandelingsopties voor deze complexe ziekten kunnen bieden.
Screening van kandidaat-hub DEG's
Het SVM-RFE-algoritme17 werd gebruikt om te screenen op potentiële cruciale DEG's, wat resulteerde in de identificatie van 6 genen in GSE89632 en 4 genen in GSE66360 (Figuur 5A,B). Evenzo werd het LASSO-algoritme19 gebruikt om kandidaat-cruciale DEG's te identificeren, wat 4 genen opleverde in GSE89632 en 5 genen in GSE66360 (Figuur 5C,D). Een Venn-diagram toonde met name de kruising van de hub-DEG (THBS1) geïdentificeerd door SVM-RFE-, LASSO- en co-PPI-netwerkanalyses, wat het belang ervan in zowel NAFLD- als MI-pathogenese onderstreepte (Figuur 5F).
Deze geavanceerde machine learning-algoritmen boden een robuust raamwerk voor het identificeren van belangrijke DEG's die verband houden met MI en NAFLD en boden waardevolle inzichten in hun potentieel als therapeutische of diagnostische doelen. De identificatie van THBS1 als een hub-gen door middel van meerdere analytische technieken benadrukt het belang ervan in de moleculaire routes die aan beide ziekten ten grondslag liggen. Verder onderzoek naar de functionele betekenis van THBS1 en de interacties ervan binnen het co-PPI-netwerk zou nieuwe mechanismen kunnen ophelderen die NAFLD en MI met elkaar verbinden, wat mogelijk de weg vrijmaakt voor gerichte therapieën en strategieën voor precisiegeneeskunde.
Diagnostische waarde van één hub DEG's
ROC-analyse werd uitgevoerd om de diagnostische efficiëntie van de geïdentificeerde hub-DEG's in zowel NAFLD- als MI-datasets te beoordelen. Zoals weergegeven in figuur 6A,B, vertoonde THBS1 een hoog onderscheidend vermogen, met een AUC van 0,981 (95% BI: 0,949-1,000) bij NAFLD-patiënten en 0,900 (95% BI: 0,834-0,956) bij MI-patiënten. De sterke prestaties van THBS1 als diagnostische biomarker benadrukken het potentiële nut ervan bij het differentiëren tussen NAFLD, MI en gezonde controles.
Deze resultaten suggereren dat THBS1 zou kunnen dienen als een waardevolle marker voor de ziektetoestand, en clinici een nuttig hulpmiddel biedt voor risicobeoordeling en vroege identificatie. Om de integratie van THBS1 in de klinische praktijk te vergemakkelijken en de patiëntenzorg en -behandeling te verbeteren, zijn verdere validatiestudies nodig om de diagnostische nauwkeurigheid ervan te bevestigen in diverse patiëntenpopulaties en klinische omgevingen.
Analyse van immuuninfiltratie
ssGSEA werd uitgevoerd om 29 immuungerelateerde genen tussen NAFLD/MI en de controlegroep (CON) te analyseren. Vergeleken met de CON-groep was de infiltratie van immuuncellen significant verhoogd in de NAFLD- en MI-groepen, met name voor CCR-, MHC-klasse I, neutrofielen, para-inflammatie en Tfh-cellen (p < 0,05, figuur 7A,B). Omgekeerd nam de infiltratie van immuuncellen aanzienlijk af in de NAFLD- en MI-groepen ten opzichte van de CON-groep, vooral voor CD8+ T-cellen, cytolytische activiteit en tumor-infiltrerende lymfocyten (TIL) (p < 0,05, Figuur 7C,D).
Deze resultaten suggereren ontregeling van de infiltratie van immuuncellen in zowel NAFLD als MI, wat het complexe immunologische landschap benadrukt dat met deze aandoeningen gepaard gaat. De waargenomen veranderingen in de samenstelling van immuuncellen bieden waardevolle inzichten in de immunologische processen die ten grondslag liggen aan de ontwikkeling en pathofysiologie van deze ziekten. Een beter begrip van deze immuunkenmerken zou nieuwe therapeutische doelen voor interventie kunnen onthullen en de ontwikkeling van geïndividualiseerde behandelplannen kunnen ondersteunen op basis van de immunologische profielen van patiënten met NAFLD en MI.

Figuur 1: De ontwerpvolgorde van het medische beslissingsondersteunende systeem voor het identificeren van THBS1. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: Venn-diagrammen, vulkaandiagrammen en heatmaps. (A) Venn-diagram voor GSE89632 en GSE66360. (B) Een vulkaandiagram van differentieel tot expressie gebrachte genen (DEG's) in GSE89632. (C) Een vulkaandiagram van DEG's in GSE66360. (D) Een heatmap-diagram van DEG's in GSE89632. (E) Een heatmap-diagram van DEG's in GSE66360. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Verrijkingsanalyse van differentieel tot expressie gebrachte genen (DEG's). (A) De top vijf NALFD-DEG's verrijkte genontologie (GO) in biologisch proces (BP), cellulaire component (CC) en moleculaire functie (MF). (B) De top vijf MI-DEG's verrijkte GO in BP, CC en MF. (C) De top acht NALFD-DEG's verrijkte Kyoto-encyclopedie van genen en genomen (KEGG) signaalroutes. (D) De top acht MI-DEG's verrijkten KEGG-signaleringsroutes. (E) De top vijf NALFD-DEG's verrijkte ziekteontologie (DO). (F) De top vijf MI-DEG's verrijkte DO. Klik hier om een grotere versie van dit cijfer te bekijken.

Figuur 4: Eiwit-eiwit interactie (PPI) netwerkdiagrammen. (A) Een PPI-netwerkdiagram van niet-alcoholische leververvetting (NAFLD) verkregen van het online platform String. (B) Een PPI-netwerkdiagram van NAFLD verkregen van Cytoscape. (C) Een PPI-netwerkdiagram van een myocardinfarct (MI) verkregen van het online platform String. (D) Een PPI-netwerkdiagram van MI verkregen Cytoscape. (E) Een PPI-netwerkdiagram van eiwitinteracties. (F) Een PPI-netwerkdiagram van IL1B en FOS. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5: Screening van DEG's in de kandidatenhub. (A) Zes cruciale DEG's verkregen door SVM-RFE in NAFLD. (B) Vier cruciale DEG's verkregen door SVM-RFE bij een myocardinfarct (MI). (C) Vier cruciale DEG's verkregen door LASSO in NAFLD. (D) Vijf cruciale DEG's verkregen door LASSO in MI. (E) Slechts één DEG verkregen uit het bovenstaande resulteert in het Venn-diagram. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: Diagnostische waarde van één hub differentieel tot expressie gebrachte genen (DEG's). (A) ROC-curven in NAFLD. (B) ROC-curven in MI. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7: Verdeling van infiltrerende immuuncellen en correlatie tussen THBS1 en infiltrerende immuuncellen bij niet-alcoholische leververvetting (NAFLD) en myocardinfarct (MI). (A) Verdeling van infiltrerende immuuncellen bij NAFLD. (B) Distributie van infiltrerende immuuncellen in MI. (C) Correlatie tussen THBS1 en infiltrerende immuuncellen in NAFLD. (D) Correlatie tussen THBS1 en infiltrerende immuuncellen in MI. Klik hier om een grotere versie van deze figuur te bekijken.
| Afkorting | Nomenclatuur |
| AUC | Gebied onder kromme |
| BP | Biologisch proces |
| CC | Cellulaire component |
| BEDRIEGEN | Besturingselementen |
| DEG | Differentieel tot expressie gebrachte genen |
| DOEN | Ontologie van de ziekte |
| GAAN | Gen ontologie |
| KEGG | Kyoto Encyclopedie van Genen en Genomen |
| LASSO | Operator Minste absolute krimp en selectie |
| MF | Moleculaire functie |
| MI | Myocardinfarct |
| NAFLD | Niet-alcoholische leververvetting |
| PPI | Eiwit-eiwit interactie |
| ROC | Ontvanger Operator Karakteristieken |
| ssGSEA | verrijkingsanalyse van genenset met één steekproef |
| SVM-RFE | Ondersteuning voor eliminatie van vectormachine-recursieve functies |
| THBS1 | Trombospondine 1 |
Tabel 1: Simulatieparameters die bij de experimenten zijn gebruikt.