Alle procedures met betrekking tot menselijk weefsel voldeden aan de institutionele richtlijnen en de Verklaring van Helsinki en werden goedgekeurd door de Institutionele Beoordelingscommissie van de Fujian Medische Universiteit (Goedkeuringsnummer 2021KYB089). Schriftelijke geïnformeerde toestemming werd van alle deelnemers verkregen voorafgaand aan de inkoop van weefsel.
Genexpressie en overlevingsanalyse
RNA-sequencinggegevens en bijbehorende klinische informatie werden verkregen uit meerdere openbare databases. 1) TCGA-cohort: RNA-seq (FPKM) gegevens voor 175 glioblastoma multiforme (GBM) en 534 laaggradige gliommonsters (LGG) zijn gedownload van The Cancer Genome Atlas (https://portal.gdc.cancer.gov/); 2) Normale controles: Expressieprofielen van 211 normale hersenweefsels en 662 gliomaweefsels werden gedownload uit de UCSC Xena-database (https://xenabrowser.net/datapages/); 3) Externe validatie: Gegevens van CGGA693 en CGGA325 cohorten zijn verkregen uit de Chinese Glioma Genome Atlas (http://www.cgga.org.cn); 4) GEO-dataset: De GSE43378-dataset, met expressie- en klinische gegevens voor 50 gliomamonsters, werd gedownload van de Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/). Alle ruwe telgegevens werden omgezet naar transcripts per miljoen (TPM) en log2 getransformeerd. Voor datasets die al genormaliseerd waren, werden expressiematrices onderzocht om vergelijkbare verdelingen te waarborgen. Genen met TPM-waarden < 1 op meer dan 80% van de steekproeven werden uitgesloten. Ontbrekende klinische informatie (leeftijd, IDH-status, 1p/19q-codeletie, MGMT-methylering) werd verwijderd met volledige casusfiltering. Batcheffecten tussen datasets werden aangepast met behulp van het ComBat-algoritme dat is geïmplementeerd in het R-pakket sva. Expressiewaarden werden gestandaardiseerd door z-score transformatie binnen elke dataset. Overlevingsanalyses werden uitgevoerd met behulp van de R-pakketten survival en survminer. Patiënten werden gedichotomiseerd in groepen met hoge en lage expressie op basis van het mediane expressieniveau van IRAIN. Kaplan-Meier overlevingscurves werden gegenereerd en de statistische significantie werd geëvalueerd met de log-rank test. Hazardratio's (HR's) en 95% betrouwbaarheidsintervallen (BI's) werden geschat met behulp van Cox-regressiemodellen met proportionele gevaren.
Definitie van immuun- en metabole gensets
Immuungerelateerde genen (IRG's, n = 2.483) werden verkregen uit de ImmPort-database (https://www.immport.org/shared/), en metabolisch gerelateerde genen (MRG's, n = 948) werden verkregen uit de Molecular Signatures Database (MSigDB, https://www.gsea-msigdb.org/). De gecombineerde set van deze genen werd gedefinieerd als immunometabolisch-gerelateerde genen (IMRG's). Deze genlijsten dienden als referenties voor latere differentiële expressie- en netwerkanalyses.
Analyse van differentiële expressie en gewogen genco-expressienetwerkanalyse
Differentieel tot expressie gebrachte genen (DEGs) tussen normaal hersen- en gliomaweefsel werden geïdentificeerd met behulp van het R-pakket limma. De expressiegegevens werden voorzien van een lineair model gevolgd door empirische Bayes-moderatie. Genen met |log₂ vouwverandering| > 1,5 en een valse ontdekkingsrate (FDR) < 0,05 werden als significant verschillend geexpresseerd beschouwd. Gewogen genco-expressienetwerkanalyse (WGCNA) werd uitgevoerd met het R-pakket WGCNA. Uitschieters werden uitgesloten door hiërarchische clustering. De soft-thresholding power werd ingesteld op β = 8 om een schaalvrije topologie-fitindex (R2 ≥ 0,85) te bereiken, terwijl de gemiddelde connectiviteit behouden bleef. Topologische overlappingsmatrices (TOM) werden geconstrueerd en genen werden gegroepeerd in modules met een minimale grootte van 50 met behulp van het dynamische boom-cut-algoritme. Module-eigengenen waren gecorreleerd met klinische kenmerken, en de module die het sterkst geassocieerd is met glioom (Pearson's r > 0,7, P < 1×10-10) werd geselecteerd voor hubgenidentificatie.
Op basis van prognosetische modelconstructie gebaseerd op machine learning
Een uitgebreid leave-one-out cross-validation (LOOCV) kader dat tien machine-learning algoritmen integreert, werd toegepast om prognostische modellen te construeren en te evalueren. In totaal werden 101 combinatorische workflows geïmplementeerd met de TCGA-cohort als trainingsdataset. Prognose-geassocieerde immunometabolisch-gerelateerde genen (IMRG's) werden voor het eerst geïdentificeerd door univariate Cox-regressie (P < 0,05). Het optimale model werd bepaald door het gemiddelde Harrell's concordantie-index (C-index) te maximaliseren over drie validatiedatasets (CGGA693, CGGA325 en GSE43378). Het resulterende RSF-Enet model (α = 0,3) toonde de hoogste voorspellende prestaties en behield robuuste generaliseerbaarheid over onafhankelijke cohorten22.
TME en immuuninfiltratie
Om het immunogenomische landschap volledig te karakteriseren, hebben we een meerlaagse analytische aanpak toegepast. Ten eerste werden immuun- en stromale infiltratieniveaus gekwantificeerd met behulp van het ESTIMATE-algoritme23. Differentiële expressie van belangrijke immuuncontrolepuntmoleculen, waaronder PDCD1, CTLA4 en LAG3, werd vervolgens beoordeeld via limma-gebaseerde analyse, en de correlaties tussen checkpointgenen werden gevisualiseerd met behulp van correlatiematrices. Somatische mutatieprofielen uit 903 glioma-monsters in de TCGA-cohort werden gebruikt om tumor mutatiebelasting (TMB), microsatellietinstabiliteit (MSI) en tumor immuundisfunctie en -uitsluiting (TIDE) scores te berekenen om mogelijke responsen op immunotherapie te voorspellen. Patiënten werden vervolgens gestratificeerd in vier prognostische groepen op basis van gecombineerde TMB-status (hoog/laag) en risicoscores (hoog/laag), en de overlevingsuitkomsten werden vergeleken met behulp van de Kaplan-Meier-analyse.
Functionele verrijkingsanalyse
Gene Ontology (GO) en Kyoto Encyclopedia of Genes and Genomes (KEGG) routeverrijkingsanalyses werden uitgevoerd met behulp van het R-pakket clusterProfiler. Verrijkingsresultaten met aangepaste P-waarden < 0,05 werden als statistisch significant beschouwd. Oververtegenwoordigde biologische processen, cellulaire componenten en moleculaire functies werden gevisualiseerd met behulp van dotplots en staafdiagrammen. Eiwit-eiwitinteractie (PPI) netwerken werden geconstrueerd met behulp van de STRING-database (≥ 0.4) en gevisualiseerd in Cytoscape. Functionele modules binnen het PPI-netwerk werden geïdentificeerd met behulp van het MCODE-algoritme. Gen-gen-interactie en co-expressienetwerken werden verder geanalyseerd met behulp van GeneMANIA (https://string-db.org; betrouwbaarheidsscore ≥ 0,4) en gevisualiseerd in Cytoscape. Functionele modules binnen het PPI-netwerk werden geïdentificeerd met behulp van het MCODE-algoritme. Gen-gen-interactie en co-expressienetwerken werden verder geanalyseerd met behulp van GeneMANIA (https://genemania.org), dat informatie integreert over fysieke en genetische interacties, gedeelde routes en co-expressiepatronen om potentiële functionele associaties af te leiden.
Klinische monsters
Verse gliomaweefsels (n = 6) en gekoppelde aangrenzende niet-tumoreuze hersenweefsels (n = 6; ten minste 3 cm van de tumorrand en histologisch bevestigd als tumorvrij) werden verzameld van patiënten die een primaire glioma-resectie ondergingen in het Zhangzhou Affiliated Hospital van de Fujian Medical University. Geen van de patiënten had vóór de operatie chemotherapie of radiotherapie ondergaan. Alle pathologische diagnoses zijn onafhankelijk geverifieerd door twee neuropathologen volgens de classificatie van tumoren van het centrale zenuwstelsel door de Wereldgezondheidsorganisatie (WHO) uit 2021. Direct na de chirurgische excisie werden weefselmonsters afgespoeld met ijskoude fosfaatgebufferde zoutoplossing (PBS) om restbloed te verwijderen, ingevroren in vloeibare stikstof (-196 °C) en bewaard bij -80 °C tot RNA-extractie.
Cellijnen en celkweek
Menselijke glioblastomcellijnen SHG44, U251, A172 en T98G, evenals normale menselijke gliacellen (HEB), werden verkregen uit geauthenticeerde bewaarplaatsen en bevestigd vrij van mycoplasma-besmetting vóór gebruik. Cellen werden behouden in Dulbecco's Modified Eagle's Medium (DMEM, hoge glucose) aangevuld met 10% foetaal rundereserum (FBS), 2 mM L-glutamine en 1% penicilline-streptomycine, bij 37 °C in een bevochtigde incubator met 5% CO₂. Cellen werden elke 4-5 dagen doorgevoerd bij het bereiken van 80-90% confluenz. Om IRAIN-overexpressieve en controlecellijnen vast te stellen, werden cellen getransduceerd met lentivirale vectoren die het volledige IRAIN-transcript droegen of een lege vector als controle. Stabiele klonen werden geselecteerd met puromycine (2 μg/mL) gedurende 14 dagen. Overexpressie-efficiëntie werd bevestigd door kwantitatieve reverse transcription PCR (qRT-PCR) voorafgaand aan downstream assays.
3- (4,5-dimethylthiazol-2-yl)-2,5-difenyltetrazoliumbromide (MTT) celproliferatietest
Cellen werden gezaaid in 96-put platen met een dichtheid van 1 × 104 cellen per put in 100 μL volledig kweekmedium. Na 24, 48 en 72 uur na het zaaien werd 20 μL MTT-oplossing (5 mg/mL in fosfaatgebufferde zoutoplossing) aan elke put toegevoegd en 4 uur geïncubeerd bij 37 °C. Het supernatant werd vervolgens verwijderd en 150 μL dimethylsulfoxide (DMSO) toegevoegd om de formazankristallen op te lossen. Het bord werd voorzichtig 10 minuten geroerd om volledige oplosbaarheid te garanderen. De absorptie werd gemeten op 490 nm met een microplaatspectrofotometer. Achtergrondmetingen van lege putten werden afgetrokken. De cellevensvatbaarheid werd berekend ten opzichte van de 24-uurs- of controlegroep (ingesteld als 1,0). Alle experimenten werden uitgevoerd met zes technische replicaten en drie onafhankelijke biologische replicaten. De gegevens worden uitgedrukt als gemiddelde ± standaarddeviatie (SD), en de statistische significantie werd bepaald met behulp van een tweezijdige t-test.
Flowcytometrie voor apoptose (Annexine V - FITC/PI-kleuring)
Cellen werden gezaaid bij 60-70% samenvloeiing en 24 uur behandeld onder de aangegeven omstandigheden. Zwevende en hechtende cellen werden verzameld met EDTA-vrije trypsine, gecombineerd en twee keer gewassen met ijskoude PBS. Celpellets werden opnieuw opgehangen in Annexin V bindingsbuffer (10 mM HEPES pH 7,4, 140 mM NaCl, 2,5 mM CaCl2) bij 1 × 106 cellen /mL. Voor elk monster werd 100 μL suspensie geïncubeerd met 5 μL Annexin V-FITC en 5 μL propidiumjodide (PI; 50 μg/mL voorraad) in het donker gedurende 15 minuten bij kamertemperatuur. Na toevoeging van 400 μL bindingsbuffer werden monsters op ijs gelegd en binnen 1 uur geanalyseerd op een flowcytometer (488 nm excitatie; 530/30 nm voor FITC en >585 nm voor PI). Geschikte enkelkleurige en fluorescentie-minus-één controles werden opgenomen ter compensatie. Er werden per steekproef minstens 10.000 gebeurtenissen geregistreerd. De gegevens werden geanalyseerd per kwadrantgating: levende (Annexin V⁻/PI⁻), vroege apoptotische (Annexin V⁺/PI⁻), late apoptotische (Annexin V⁺/PI⁺) en necrotische (Annexin V⁻/PI⁺) populaties. Percentages vroege + late apoptotische cellen werden gerapporteerd (gemiddelde ± SD, n = 3).
Kwantitatieve real-time PCR (qRT-PCR)
Total RNA werd geïsoleerd met behulp van een zuur fenol-guanidiniumreagens volgens het protocol van de fabrikant. De zuiverheid van RNA werd geverifieerd door spectrofotometrie (A₂₆₀/A₂₈₀ = 1,8-2,1), en de integriteit werd bevestigd door gelelektroforese (RNA-integriteitsnummer ≥ 7). Eén microgram totaal RNA werd behandeld met DNase I en omgekeerd getranscribeerd in een 20 μL-reactie met behulp van willekeurige hexamers en oligo(dT)-primers. De reactie werd uitgevoerd bij 25 °C gedurende 10 minuten, 50 °C gedurende 30 minuten en 85 °C gedurende 5 minuten. Kwantitatieve PCR werd uitgevoerd in een 10 μL-systeem met 5 μL 2× SYBR Green Master Mix, elk 0,3 μM primer, en 1 μL cDNA (≈ 20 ng RNA-equivalent). De thermische cycluscondities waren 95 °C gedurende 5 minuten, gevolgd door 40 cycli van 95 °C gedurende 15 seconden en 60 °C gedurende 30 seconden, gevolgd door een smeltcurveanalyse van 65 °C tot 95 °C in stappen van 0,3 °C. Alle reacties werden driemaal uitgevoerd, samen met no-template en minus-RT controles. CT-waarden > 35 of technisch replicate SD > 0,5 werden uitgesloten. De relatieve expressie werd berekend met de 2⁻ΔΔCt-methode, met GAPDH als interne controle. Gemiddelde ± SD-waarden van drie onafhankelijke biologische replicaten werden gerapporteerd, en groepsverschillen werden geanalyseerd met een tweezijdige t-test.
Western blot-analyse
Cellen werden gelyseerd op ijs in een RIPA-buffer (50 mM Tris-HCl, pH 7,4, 150 mM NaCl, 1% NP-40, 0,5% natriumdeoxycholaat, 0,1% SDS) aangevuld met protease- en fosfataseremmers. Lysaten werden 30 minuten op ijs geïncubeerd met intermitterende vortexen en verwijderd door centrifugatie bij 12.000 × g gedurende 15 minuten bij 4 °C. De eiwitconcentraties werden gemeten met een BCA-assay, aangepast naar 1-2 μg/μL, en gemengd 1:3 met 4× Laemmli-buffer (de laatste 1× buffer met 100 mM DTT). De monsters werden 5 minuten gedenatureerd bij 95 °C. Gelijke hoeveelheden eiwit (50 μg) werden opgelost met 12% SDS-PAGE bij 100 V gedurende 90 minuten en elektro-overgebracht naar PVDF-membranen bij 250 mA gedurende 90 minuten. De membranen werden geblokkeerd met 5% magere melk in TBST (0,1% Tween-20) gedurende 1 uur bij kamertemperatuur (of 5% BSA voor fosfoproteïnen) en 's nachts geïncubeerd bij 4 °C met primaire antilichamen tegen IGF1, IGF1R, JAK2, p-JAK2 (Y1007/1008), STAT3, p-STAT3 (Y705), BIRC5 en β-actine (typische verdunning 1:1000, β-actine 1:5000). Na drie wasbeurten van 10 minuten in TBST werden membranen geïncubeerd met HRP-geconjugeerde secundaire antistoffen (1:5000) gedurende 1 uur bij kamertemperatuur, opnieuw gewassen en ontwikkeld met chemiluminescente substraat. Bandintensiteiten werden gekwantificeerd met ImageJ, genormaliseerd tot β-actine of totaal eiwit, en uitgedrukt als gemiddelde ± SD uit drie onafhankelijke experimenten.
Immunocytochemie
Cellen die op steriele glazen dekfolies werden gekweekt, werden twee keer afgespoeld met PBS en 15 minuten bij kamertemperatuur vastgezet in 4% paraformaldehyde. Na drie PBS-wasbeurten werden de cellen permeabiliseerd met 0,2% Triton X-100 gedurende 10 minuten, geblokkeerd met 5% rundereserumalbumine (BSA) gedurende 1 uur, en 's nachts geïncubeerd bij 4 °C met het primaire anti-CD31-antilichaam (1:200 verdunning in 1% BSA). Na drie PBS-washes werden cellen geïncubeerd met Alexa-fluor-geconjugeerde secundaire antilichamen (1:500 verdunning) gedurende 1 uur in het donker, tegengekleurd met DAPI (1 μg/mL, 5 min) en gemonteerd in antifade medium. Beelden werden vastgelegd met een fluorescentiemicroscoop onder identieke belichtings- en versterkingsinstellingen. Het percentage CD31-positieve gebied werd gekwantificeerd in vijf willekeurig geselecteerde niet-overlappende velden per steekproef met behulp van ImageJ-software. Deze assay werd uitgevoerd in celmodellen in plaats van in weefseldoorsneden.
Statistische analyse
Statistische analyses werden uitgevoerd met behulp van R-versie 4.3.0 samen met de bijbehorende pakketten. Om categorische variabelen te vergelijken werd de chi-kwadraattest gebruikt, terwijl continue variabelen werden beoordeeld met behulp van ofwel de Wilcoxon-rangsomtest of de T-test. De evaluatie van continue variabelen werd uitgevoerd via Pearsons correlatiecoëfficiënt. Overlevingsanalyses werden uitgevoerd met behulp van het survivalpakket, dat Cox-proportionele gevarenmodellering en het genereren van Kaplan-Meier-curves omvatte, met optimale stratificatiedrempels vastgesteld door het survminer-pakket en de formule Riskscore =
. Het CompareC-pakket werd gebruikt om de C-indices van verschillende variabelen te beoordelen. De ontvanger-operating characteristic curve (ROC), bedoeld om binaire categorische variabelen te voorspellen, werd gegenereerd met het pROC-pakket. Daarnaast werd het tijdsafhankelijke gebied onder de ROC-curve (AUC) voor survivalmetrieken geanalyseerd met behulp van het timeROC-pakket. Alle statistische tests werden uitgevoerd met een tweezijdige aanpak. Een significantieniveau van P < 0,05 werd als statistisch significant beschouwd.