Alle Verfahren mit menschlichem Gewebe entsprachen den institutionellen Richtlinien und der Erklärung von Helsinki und wurden vom Institutional Review Board der Fujian Medical University genehmigt (Genehmigungsnr. 2021KYB089). Vor der Gewebebeschaffung wurde von allen Teilnehmern eine schriftliche informierte Zustimmung eingeholt.
Genexpression und Überlebensanalyse
RNA-Sequenzierungsdaten und entsprechende klinische Informationen wurden aus mehreren öffentlichen Datenbanken bezogen. 1) TCGA-Kohorte: RNA-seq (FPKM)-Daten für 175 Glioblastom-Multiforme- (GBM)- und 534 Low-Grade-Gliom-(LGG)-Proben wurden aus dem Cancer Genome Atlas (https://portal.gdc.cancer.gov/) heruntergeladen; 2) Normale Kontrollgruppen: Expressionsprofile von 211 normalen Hirngeweben und 662 Gliomgeweben wurden aus der UCSC Xena-Datenbank heruntergeladen (https://xenabrowser.net/datapages/); 3) Externe Validierung: Daten aus CGGA693 und CGGA325 Kohorten wurden aus dem Chinese Glioma Genome Atlas (http://www.cgga.org.cn) gewonnen; 4) GEO-Datensatz: Der GSE43378-Datensatz, der Expressions- und klinische Daten für 50 Gliomproben enthält, wurde aus dem Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) heruntergeladen. Alle Rohzähldaten wurden in Transkripte pro Million (TPM) umgewandelt und log2 transformiert. Für bereits normalisierte Datensätze wurden Expressionsmatrizen untersucht, um vergleichbare Verteilungen sicherzustellen. Gene mit TPM-Werten < 1 von mehr als 80 % der Proben wurden ausgeschlossen. Fehlende klinische Informationen (Alter, IDH-Status, 1p/19q-Codeletion, MGMT-Methylierung) wurden mittels vollständiger Fall-Filterung entfernt. Batch-Effekte in Datensätzen wurden mit dem ComBat-Algorithmus angepasst, der im R-Paket sva implementiert wurde. Die Ausdruckswerte wurden durch die Z-Score-Transformation innerhalb jedes Datensatzes standardisiert. Überlebensanalysen wurden mit den R-Paketen survival und survminer durchgeführt. Die Patienten wurden entsprechend dem medianen Expressionsniveau von IRAIN in Hoch- und Niedrigausdrucksgruppen eingeteilt. Kaplan-Meier-Überlebenskurven wurden generiert, und die statistische Signifikanz wurde durch den Log-Rank-Test bewertet. Hazard-Verhältnisse (HRs) und 95%-Konfidenzintervalle (KI) wurden mittels Cox-Proportionalhazard-Regressionsmodellen geschätzt.
Definition von Immun- und Stoffwechselgensets
Immunbezogene Gene (IRGs, n = 2.483) wurden aus der ImmPort-Datenbank (https://www.immport.org/shared/) gewonnen, und metabolisch verwandte Gene (MRGs, n = 948) wurden aus der Molecular Signatures Database (MSigDB, https://www.gsea-msigdb.org/) gewonnen. Der kombinierte Satz dieser Gene wurde als immunmetabolisch verwandte Gene (IMRGs) definiert. Diese Genlisten dienten als Referenzen für nachfolgende differentielle Expressions- und Netzwerkanalysen.
Analyse des differentiellen Expressionsnetzwerks und des gewichteten Gen-Koexpressionsnetzwerks
Differenziell exprimierte Gene (DEGs) zwischen normalem Gehirn- und Gliomgewebe wurden mit dem R-Paket-Limma identifiziert. Die Expressionsdaten wurden mit einem linearen Modell ausgestattet, gefolgt von empirischer Bayes-Moderation. Gene mit |log₂ Faltenänderung| > 1,5 und eine False Discovery Rate (FDR) < 0,05 wurden als signifikant unterschiedlich ausgedrückt angesehen. Die Analyse des gewichteten Gen-Koexpressionsnetzwerks (WGCNA) wurde mit dem R-Paket WGCNA durchgeführt. Ausreißerstichproben wurden durch hierarchische Clusterung ausgeschlossen. Die Soft-Thresholding-Potenz wurde auf β = 8 gesetzt, um einen skalenfreien Topologie-Fit-Index (R2 ≥ 0,85) zu erreichen, während eine ausreichende mittlere Konnektivität erhalten bleibt. Topologische Überlappungsmatrizen (TOM) wurden konstruiert, und Gene wurden mithilfe des dynamischen Baumschnitt-Algorithmus in Module mit einer Mindestgröße von 50 Modulen gruppiert. Die Eigengene des Moduls korrelierten mit klinischen Merkmalen, und das am stärksten mit Gliom assoziierte Modul (Pearsons r > 0,7, P < 1×10-10) wurde für die Hubgen-Identifikation ausgewählt.
Maschinell-Learning-basierte prognostische Modellkonstruktion
Ein umfassendes Leave-One-out Cross-Validation (LOOCV)-Framework, das zehn maschinelle Lernalgorithmen integriert, wurde angewandt, um prognostische Modelle zu erstellen und zu evaluieren. Insgesamt wurden 101 kombinatorische Workflows unter Verwendung der TCGA-Kohorte als Trainingsdatensatz implementiert. Prognosassoziierte immunmetabolisch bezogene Gene (IMRGs) wurden erstmals durch eine univariate Cox-Regression (P < 0,05) identifiziert. Das optimale Modell wurde ermittelt, indem der mittlere Harrell-Konkordanzindex (C-Index) über drei Validierungsdatensätze (CGGA693, CGGA325 und GSE43378) maximiert wurde. Das daraus resultierende RSF-Enet-Modell (α = 0,3) zeigte die höchste prädiktive Leistung und bewahrte eine robuste Generalisierbarkeit über unabhängige Kohortenhinweg 22.
TME und Immuninfiltration
Um die immunogenomische Landschaft umfassend zu charakterisieren, haben wir einen mehrstufigen analytischen Ansatz angewandt. Zunächst wurden die Immun- und Strominfiltrationsniveaus mit dem ESTIMATE-Algorithmus23 quantifiziert. Die differenzielle Expression wichtiger Immunkontrollpunktmoleküle, darunter PDCD1, CTLA4 und LAG3, wurde anschließend mittels Limma-basierter Analyse bewertet, und die Korrelationen zwischen Kontrollpunktgenen wurden mithilfe von Korrelationsmatrizen visualisiert. Somatische Mutationsprofile aus 903 Gliomproben der TCGA-Kohorte wurden verwendet, um Tumormutationslast (TMB), Mikrosatelliteninstabilität (MSI) sowie Tumor-Immunfunktionsstörungen und -ausschlusswerte (TIDE) zu berechnen, um potenzielle Reaktionen auf Immuntherapie vorherzusagen. Die Patienten wurden anschließend in vier prognostische Gruppen nach kombiniertem TMB-Status (hoch/niedrig) und Risikowerten (hoch/niedrig) eingeteilt, und die Überlebensergebnisse wurden mittels der Kaplan-Meier-Analyse verglichen.
Funktionale Anreicherungsanalyse
Gene Ontology (GO) und Kyoto Encyclopedia of Genes and Genomes (KEGG) Signalweganreicherungsanalysen wurden mit dem R-Paket clusterProfiler durchgeführt. Anreicherungsergebnisse mit bereinigten P-Werten < 0,05 wurden als statistisch signifikant angesehen. Überrepräsentierte biologische Prozesse, zelluläre Komponenten und molekulare Funktionen wurden mithilfe von Punktdiagrammen und Balkendiagrammen visualisiert. Protein-Protein-Interaktionsnetzwerke (PPI) wurden mit der STRING-Datenbank (≥ 0.4) konstruiert und in Cytoscape visualisiert. Funktionale Module innerhalb des PPI-Netzwerks wurden mithilfe des MCODE-Algorithmus identifiziert. Gen-Gen-Interaktions- und Koexpressionsnetzwerke wurden mit GeneMANIA (https://string-db.org; Konfidenzwert ≥ 0,4) weiter analysiert und in Cytoscape visualisiert. Funktionale Module innerhalb des PPI-Netzwerks wurden mithilfe des MCODE-Algorithmus identifiziert. Gen-Gen-Interaktions- und Koexpressionsnetzwerke wurden mit GeneMANIA (https://genemania.org) weiter analysiert, das Informationen über physische und genetische Interaktionen, gemeinsame Signalwege und Koexpressionsmuster integriert, um potenzielle funktionelle Assoziationen abzuleiten.
Klinische Proben
Frisches Gliomgewebe (n = 6) und gepaarte benachbarte nicht-tumoröse Hirngewebe (n = 6; mindestens 3 cm vom Tumorrand entfernt und histologisch als tumorfrei bestätigt) wurden von Patienten entnommen, die sich einer primären Gliomresektion im Zhangzhou Affiliated Hospital der Fujian Medical University unterzogen. Keiner der Patienten hatte vor der Operation eine Chemo- oder Strahlentherapie erhalten. Alle pathologischen Diagnosen wurden unabhängig von zwei Neuropathologen gemäß der Klassifikation von Tumoren im zentralen Nervensystem der Weltgesundheitsorganisation (WHO) 2021 bestätigt. Unmittelbar nach der chirurgischen Exzision wurden Gewebeproben mit eiskaltem phosphatgepuffertem Kochsalzlösungslösungszustand (PBS) abgespült, um Restblut zu entfernen, in flüssigem Stickstoff (-196 °C) gelassen und bei -80 °C bis zur RNA-Extraktion gelagert.
Zelllinien und Zellkultur
Die humanen Glioblastomzelllinien SHG44, U251, A172 und T98G sowie normale humane Gliazellen (HEB) wurden aus authentifizierten Repositorien gewonnen und vor der Anwendung als frei von Mycoplasma-Kontamination bestätigt. Die Zellen wurden in Dulbeccos Modified Eagle's Medium (DMEM, hoher Blutzucker) mit 10 % fetalem Rinderserum (FBS), 2 mM L-Glutamin und 1 % Penicillin-Streptomycin bei 37 °C in einem befeuchteten Inkubator mit 5 % CO₂ erhalten. Die Zellen wurden alle 4-5 Tage nach 80-90 % Konfluenz durchführt. Um IRAIN-überexpressive und Kontrollzelllinien zu etablieren, wurden Zellen mit lentiviralen Vektoren transduziert, die das vollständige IRAIN-Transkript oder einen leeren Vektor als Kontrollvektor trugen. Stabile Klone wurden 14 Tage lang mit Puromycin (2 μg/mL) ausgewählt. Die Überexpressionseffizienz wurde durch quantitative Reverse-Transkriptions-PCR (qRT-PCR) vor den nachgelagerten Tests bestätigt.
3- (4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazoliumbromid (MTT) Zellproliferationstest
Die Zellen wurden in 96-Well-Platten mit einer Dichte von 1 × 104 Zellen pro Bohrloch in 100 μL vollständigem Kulturmedium ausgesaat. Nach 24, 48 und 72 Stunden nach der Aussaat wurden 20 μL MTT-Lösung (5 mg/ml in phosphatgepufferter Kochsalzlösung) zu jedem Brunnen hinzugefügt und 4 Stunden bei 37 °C inkubiert. Anschließend wurde das Supernatant entfernt, und 150 μL Dimethylsulfoxid (DMSO) wurden hinzugefügt, um die Formazankristalle zu lösen. Der Teller wurde 10 Minuten lang vorsichtig bewegt, um eine vollständige Lösungslösung sicherzustellen. Die Absorption wurde bei 490 nm mit einem Mikroplattenspektrophotometer gemessen. Hintergrundmessungen aus leeren Quellen wurden abgezogen. Die Zelllebensfähigkeit wurde relativ zur 24-Stunden- oder Kontrollgruppe berechnet (gesetzt als 1,0). Alle Experimente wurden mit sechs technischen Replikaten und drei unabhängigen biologischen Replikaten durchgeführt. Die Daten werden als Mittelwert ± Standardabweichung (SD) angegeben, und die statistische Signifikanz wurde mit einem zweiseitigen t-Test bestimmt.
Durchflusszytometrie für Apoptose (Annexin V - FITC/PI-Färbung)
Die Zellen wurden bei 60–70 % Konfluenz ausgesät und unter den angegebenen Bedingungen 24 Stunden lang behandelt. Schwimmende und haftende Zellen wurden mit EDTA-freiem Trypsin gesammelt, kombiniert und zweimal mit eiskalten PBS gewaschen. Zellpellets wurden im Annexin-V-Bindungspuffer (10 mM HEPES pH 7,4, 140 mM NaCl, 2,5 mM CaCl2) bei 1 × 106 Zellen/ml wieder suspendiert. Für jede Probe wurden 100 μL Suspension mit 5 μL Annexin V-FITC und 5 μL Propidiumiodid (PI; 50 μg/mL Stock) im Dunkeln für 15 Minuten bei Raumtemperatur inkubiert. Nach der Zugabe von 400 μL Bindungspuffer wurden die Proben auf Eis gelegt und innerhalb von 1 Stunde auf einem Durchflusszytometer analysiert (488 nm Anregung; 530/30 nm für FITC und >585 nm für PI). Geeignete Einzelfärbungs- und Fluoreszenz-minus-eins-Kontrollen wurden zur Kompensation einbezogen. Mindestens 10.000 Ereignisse pro Stichprobe wurden aufgezeichnet. Die Daten wurden nach Quadrantengating analysiert: lebende (Annexin V⁻/PI⁻), frühe apoptotische (Annexin V⁺/PI⁻), spätapoptotische (Annexin V⁺/PI⁺) und nekrotische (Annexin V⁻/PI⁺) Populationen. Prozentsätze früher + später apoptotischer Zellen wurden berichtet (Mittelwert ± SD, n = 3).
Quantitative Echtzeit-PCR (qRT-PCR)
Total RNA wurde gemäß dem Protokoll des Herstellers mit einem sauren Phenol-Guanidinium-Reagenz isoliert. Die RNA-Reinheit wurde durch Spektrophotometrie (A₂₆₀/A₂₈₀ = 1,8-2,1) bestätigt, und die Integrität wurde durch Gelelektrophorese bestätigt (RNA-Integritätsnummer ≥ 7). Ein Mikrogramm Gesamt-RNA wurde mit DNase I behandelt und in einer 20-μL-Reaktion mit zufälligen Hexamern und Oligo(dT)-Primern rückwärts transkribiert. Die Reaktion erfolgte bei 25 °C für 10 Minuten, 50 °C für 30 Minuten und 85 °C für 5 Minuten. Die quantitative PCR wurde in einem 10-μL-System durchgeführt, das 5 μL 2× SYBR Green Master Mix, je 0,3 μM Primer und 1 μL cDNA (≈ 20 ng RNA-Äquivalent) enthielt. Die thermischen Zyklusbedingungen betrugen 95 °C für 5 Minuten, gefolgt von 40 Zyklen mit 95 °C für 15 Sekunden und 60 °C für 30 Sekunden, gefolgt von einer Schmelzkurvenanalyse von 65 °C bis 95 °C in 0,3 °C-Schritten. Alle Reaktionen wurden dreifach durchgeführt, zusammen mit No-Template- und Minus-RT-Steuerungen. CT-Werte > 35 oder technisch replikierte SD > 0,5 wurden ausgeschlossen. Die relative Expression wurde mit der 2⁻ΔΔCt-Methode berechnet, wobei GAPDH die interne Kontrolle war. Mittlere ± SD-Werte von drei unabhängigen biologischen Replikaten wurden berichtet, und Gruppenunterschiede wurden mit einem zweiseitigen t-Test analysiert.
Western-Blot-Analyse
Die Zellen wurden auf Eis im RIPA-Puffer (50 mM Tris-HCl, pH 7,4, 150 mM NaCl, 1 % NP-40, 0,5 % Natriumdeoxycholat, 0,1 % SDS) mit Protease- und Phosphataseinhibitoren aufgewertet. Lysate wurden 30 Minuten auf Eis mit intermittierendem Vortex inkubiert und durch Zentrifugation bei 12.000 × g für 15 Minuten bei 4 °C beseitigt. Die Proteinkonzentrationen wurden durch einen BCA-Test gemessen, angepasst auf 1-2 μg/μL und 1:3 mit dem 4× Laemmli-Puffer (der letzte 1×-Puffer mit 100 mM DTT) gemischt. Die Proben wurden bei 95 °C für 5 Minuten denaturiert. Gleiche Mengen an Protein (50 μg) wurden durch 12 % SDS-PAGE bei 100 V für 90 Minuten aufgelöst und für 90 Minuten auf PVDF-Membranen bei 250 mA elektrotransferiert. Die Membranen wurden mit 5 % fettfreier Milch in TBST (0,1 % Tween-20) für 1 Stunde bei Raumtemperatur (oder 5 % BSA für Phosphoproteine) blockiert und über Nacht bei 4 °C mit primären Antikörpern gegen IGF1, IGF1R, JAK2, p-JAK2 (Y1007/1008), STAT3, p-STAT3 (Y705), BIRC5 und β-Aktin (typische Verdünnung 1:1000, β-Aktin 1:5000) inkubiert. Nach drei 10-minütigen Waschen in TBST wurden die Membranen mit HRP-konjugierten Sekundärantikörpern (1:5000) bei Raumtemperatur für 1 Stunde inkubiert, erneut gewaschen und mit chemilumineszentem Substrat entwickelt. Die Bandintensitäten wurden mit ImageJ quantifiziert, auf β-Aktin oder Gesamtprotein normalisiert und als Mittelwert ± SD aus drei unabhängigen Experimenten ausgedrückt.
Immunzytochemie
Zellen, die auf sterilen Glasdeckfolien gezüchtet wurden, wurden zweimal mit PBS ausgespült und bei Zimmertemperatur 15 Minuten lang in 4 % Paraformaldehyd fixiert. Nach drei PBS-Waschungen wurden die Zellen mit 0,2 % Triton X-100 für 10 Minuten permeabilisiert, 1 Stunde lang mit 5 % Rinderserumalbumin (BSA) blockiert und über Nacht bei 4 °C mit dem primären Anti-CD31-Antikörper (1:200 Verdünnung in 1 % BSA) inkubiert. Nach drei PBS-Waschverfahren wurden die Zellen mit Alexa-fluor-konjugiertem Sekundärantikörper (1:500-Verdünnung) für 1 Stunde im Dunkeln inkubiert, mit DAPI (1 μg/mL, 5 Min) gekonnt gefärbt und in Antifade-Medium montiert. Die Bilder wurden mit einem Fluoreszenzmikroskop unter identischen Belichtungs- und Verstärkungseinstellungen aufgenommen. Der Prozentsatz der CD31-positiven Fläche wurde in fünf zufällig ausgewählten, nicht überlappenden Feldern pro Stichprobe mit der ImageJ-Software quantifiziert. Dieser Test wurde in Zellmodellen statt in Gewebeschnitten durchgeführt.
Statistische Analyse
Statistische Analysen wurden unter Verwendung von R-Version 4.3.0 sowie den zugehörigen Paketen durchgeführt. Um kategorische Variablen zu vergleichen, wurde der Chi-Quadrat-Test verwendet, während kontinuierliche Variablen entweder mit dem Wilcoxon-Rangsummentest oder dem T-Test bewertet wurden. Die Bewertung kontinuierlicher Variablen erfolgte über den Pearsonschen Korrelationskoeffizienten. Überlebensanalysen wurden mit dem Überlebenspaket durchgeführt, das Cox-Proportionalgefahrmodellierungen und die Generierung von Kaplan-Meier-Kurven umfasste, wobei optimale Stratifikationsschwellenwerte durch das Survminer-Paket und die Formel Riskscore =
festgelegt wurden. Das CompareC-Paket wurde verwendet, um die C-Indizes verschiedener Variablen zu bewerten. Die Empfänger-Betriebscharakteristiekkurve (ROC), die darauf abzielt, binäre Kategorienvariablen vorherzusagen, wurde mit dem pROC-Paket generiert. Zusätzlich wurde die zeitabhängige Fläche unter der ROC-Kurve (AUC) für Überlebensmetriken mit dem timeROC-Paket analysiert. Alle statistischen Tests wurden mit einem zweiseitigen Ansatz durchgeführt. Ein Signifikanzniveau von P < 0,05 wurde als statistisch signifikant angesehen.