İnsan dokularıyla ilgili tüm prosedürler, kurumsal yönergeler ve Helsinki Bildirgesi'ne uygun olup, Fujian Tıp Üniversitesi Kurumsal İnceleme Kurulu tarafından onaylandı (Onay No. 2021KYB089). Tüm katılımcılardan doku alımından önce yazılı bilgilendirilmiş onay alındı.
Gen ifadesi ve hayatta kalma analizi
RNA dizileme verileri ve ilgili klinik bilgiler birden fazla kamuya açık veri tabanından elde edilmiştir. 1) TCGA kohortu: 175 glioblastoma multiforme (GBM) ve 534 düşük dereceli glioma (LGG) örneği için RNA-seq (FPKM) verileri The Cancer Genome Atlas (https://portal.gdc.cancer.gov/)'dan indirildi; 2) Normal kontroller: UCSC Xena veritabanından 211 normal beyin dokusu ve 662 glioma dokusunun ifade profilleri indirildi (https://xenabrowser.net/datapages/); 3) Dış doğrulama: CGGA693 ve CGGA325 kohortlarından veriler Çin Glioma Genom Atlası'ndan (http://www.cgga.org.cn) alındı; 4) GEO veri seti: 50 glioma örneği için ifade ve klinik verileri içeren GSE43378 veri seti, Gen İfadesi Omnibus'tan (https://www.ncbi.nlm.nih.gov/geo/) indirildi. Tüm ham sayım verileri milyonda transkriptlere (TPM) dönüştürüldü ve log2-dönüşümlü oldu. Zaten normalleştirilmiş veri setleri için, karşılaştırılabilir dağılımlar sağlamak amacıyla ifade matrisleri incelendi. TPM değeri %80'den fazla %80'den fazlasında 1 < genler hariç tutuldu. Eksik klinik bilgiler (yaş, IDH durumu, 1p/19q kodelesiyonu, MGMT metilasyonu) tam vaka filtreleme ile kaldırıldı. Veri setleri arasındaki toplu etkiler, R paketi sva'da uygulanan ComBat algoritması kullanılarak ayarlandı. İfade değerleri, her veri setinde z-skor dönüşümü ile standartlaştırıldı. Sağ kalma analizleri, R paketleri olan hayatta kalma ve survminer kullanılarak gerçekleştirildi. Hastalar, IRAIN'in medyan ekspresyon seviyesine göre yüksek ve düşük ekspresyonlu gruplara ayrıldı. Kaplan-Meier hayatta kalma eğrileri oluşturuldu ve istatistiksel anlamlılık log-rank testi ile değerlendirildi. Tehlike oranları (HR'ler) ve %95 güven aralıkları (CI'ler) Cox orantılı tehlike regresyon modelleri kullanılarak tahmin edilmiştir.
Bağışıklık ve metabolik gen kümelerinin tanımı
Bağışıklık ilişkili genler (IRG, n = 2.483) ImmPort veritabanından (https://www.immport.org/shared/) alındı, metabolik ilgili genler (MRGs, n = 948) ise Moleküler İmza Veritabanı'ndan (MSigDB, https://www.gsea-msigdb.org/) alındı. Bu genlerin birleşimi immünometabolik ilişkili genler (IMRG) olarak tanımlandı. Bu gen listeleri, sonraki diferansiyel ifade ve ağ analizleri için referans olarak hizmet etti.
Diferansiyel ifade ve ağırlıklı gen ortak ekspresyon ağı analizi
Normal beyin ve glioma dokuları arasında farklı ekspresyoner genler (DEG) R paketi limma kullanılarak tanımlandı. İfade verileri, doğrusal bir model ve ardından ampirik Bayes moderasyonu ile donatıldı. |log₂ katlama değişimi olan genler| > 1.5 ve yanlış keşf oranı (FDR) < 0.05 anlamlı şekilde farklı olarak ifade edilmişti. Ağırlıklı gen ortak ekspresyon ağı analizi (WGCNA), R paketi WGCNA kullanılarak gerçekleştirildi. Dış örnekler hiyerarşik kümeleme yoluyla dışlandı. Yumuşak eşik gücü, ölçeksiz bir topoloji uyum indeksi (R2 ≥ 0.85) elde etmek için β = 8 olarak ayarlandı ve yeterli ortalama bağlantıyı korudu. Topolojik örtüşme matrisleri (TOM) oluşturuldu ve genler, dinamik ağaç kesim algoritması kullanılarak minimum 50 boyutlu modüllere gruplandı. Modül özgenleri klinik özelliklerle ilişkilendirildi ve glioma ile en güçlü ilişkili modül (Pearson r > 0.7, P < 1×10-10) hub gen tanımlaması için seçildi.
Makine öğrenimi tabanlı prognostik model inşası
Prognostik modeller oluşturmak ve değerlendirmek için on makine öğrenimi algoritmasını entegre eden kapsamlı bir baskı (LOOCV) çerçevesi uygulandı. Toplamda, TCGA kohortu eğitim veri seti olarak kullanılarak 101 kombinatoryal iş akışı uygulandı. Prognozla ilişkili immünometabolik genler (IMRG'ler) ilk olarak tek değişkenli Cox regresyonuyla (P < 0.05) tanımlanmıştır. Optimal model, üç doğrulama veri setinde (CGGA693, CGGA325 ve GSE43378) ortalama Harrell'in uyum indeksinin (C-indeksi) maksimize edilmesiyle belirlendi. Ortaya çıkan RSF-Enet modeli (α = 0.3) en yüksek tahmin performansını gösterdi ve bağımsızkohortlar 22 arasında sağlam genellenebilirliği korudu.
TME ve bağışıklık sızması
İmmünogenomik dünyayı kapsamlı bir şekilde tanımlamak için çok katmanlı bir analitik yaklaşım kullandık. İlk olarak, bağışıklık ve stromal infiltrasyon seviyeleri ESTIMATEalgoritması 23 kullanılarak nicelendirildi. PDCD1, CTLA4 ve LAG3 dahil olmak üzere anahtar bağışıklık kontrol noktası moleküllerinin diferansiyel ekspresyonu limma tabanlı analizle değerlendirildi ve kontrol noktası genleri arasındaki korelasyonlar korelasyon matrisleri kullanılarak görselleştirildi. TCGA kohortunda 903 glioma örneğinden alınan somatik mutasyon profilleri, tümör mutasyon yükü (TMB), mikrosatellit instabilitesi (MSI) ve tümör bağışıklık disfonksiyonu ve dışlanma (TIDE) skorlarını hesaplamak için kullanıldı ve immünoterapiye potansiyel yanıtlar tahmin edildi. Hastalar daha sonra birleşik TMB durumu (yüksek/düşük) ve risk puanları (yüksek/düşük) olarak dört prognostik gruba ayrıldı ve sağlıktaki sonuçlar Kaplan-Meier analiziyle karşılaştırıldı.
Fonksiyonel zenginleştirme analizi
Gen Ontolojisi (GO) ve Kyoto Genler ve Genomlar Ansiklopedisi (KEGG) yol zenginleştirme analizleri R paket kümesi Profiler kullanılarak gerçekleştirildi. Ayarlanmış P değerleri < 0,05 ile zenginleştirme sonuçları istatistiksel olarak anlamlı olarak kabul edildi. Aşırı temsil edilen biyolojik süreçler, hücresel bileşenler ve moleküler fonksiyonlar, nokta grafikleri ve çubuk grafikleri kullanılarak görselleştirildi. Protein-protein etkileşimi (PPI) ağları STRING veritabanı (≥ 0.4) kullanılarak oluşturulmuş ve Cytoscape'te görselleştirilmiştir. PPI ağı içindeki fonksiyonel modüller MCODE algoritması kullanılarak tanımlandı. Gen-gen etkileşimi ve ortak ifade ağları GeneMANIA ile daha fazla analiz edildi (https://string-db.org; güven puanı ≥ 0.4) ve Cytoscape'de görselleştirildi. PPI ağı içindeki fonksiyonel modüller MCODE algoritması kullanılarak tanımlandı. GeneMANIA (https://genemania.org) ile gen-gen etkileşimi ve ortak ifade ağları, fiziksel ve genetik etkileşimler, paylaşılan yollar ve ortak ifade kalıpları hakkında bilgileri entegre ederek potansiyel fonksiyonel ilişkiler çıkaran GeneMANIA () kullanılarak daha fazla analiz edilmiştir.
Klinik örnekler
Fujian Tıp Üniversitesi'nin Zhangzhou Bağlı Hastanesi'nde primer glioma rezeksiyonu geçiren hastalardan taze glioma dokuları (n = 6) ve çift bitişik tümör olmayan beyin dokuları (n = 6; tümör kenarından en az 3 cm uzakta bulunan ve histolojik olarak tümörden arınmış olduğu doğrulanmış) toplandı. Hastaların hiçbiri ameliyattan önce kemoterapi veya radyoterapi almamıştı. Tüm patolojik tanılar, 2021 Dünya Sağlık Örgütü (WHO) merkezi sinir sistemi tümörleri sınıflandırmasına göre iki nöropatolog tarafından bağımsız olarak doğrulandı. Cerrahi eksizyondan hemen sonra, doku örnekleri buz gibi fosfatla tamponlanmış tuzlu su (PBS) ile durulandırılarak kalan kan alındı, sıvı azotta (-196 °C) donduruldu ve -80 °C'de saklanarak RNA çıkarımı yapıldı.
Hücre hatları ve hücre kültürü
İnsan glioblastoma hücre hatları SHG44, U251, A172 ve T98G ile normal insan glial hücreleri (HEB), doğrulamalı depolardan alınmış ve kullanımdan önce mikoplazma kontaminasyonundan arındırılmış olduğu doğrulanmıştır. Hücreler, %10 fetal sığır serumu (FBS), 2 mM L-glutamin ve %1 penisillin-streptomisin ile takviye edilmiş Dulbecco'nun Modified Eagle's Medium (DMEM, yüksek glikoz) içinde nemlendirilmiş bir kuluçka içinde, %5 CO₂ ile nemlendirilmiş bir kuluçka makinesinde tutuldu. Hücreler %80-90 birleşime ulaştıklarında her 4-5 günde bir geçiş yapıldı. IRAAN aşırı ekspresyona ve kontrol hücre hatlarını oluşturmak için, hücreler tam uzunlukta IRAIN transkripti taşıyan lentiviral vektörlerle veya kontrol olarak boş bir vektörle transdüksiyon yapıldı. Stabil klonlar 14 gün boyunca puromisin (2 μg/mL) kullanılarak seçildi. Aşırı ifade verimliliği, sonraki testlerden önce nicel ters transkripsiyon PCR (qRT-PCR) ile doğrulandı.
3- (4,5-dimetiltiazol-2-il)-2,5-difeniltetrazolyum bromid (MTT) hücre proliferasyon testi
Hücreler, 100 μL tam kültür ortamında 1 ×10-4 hücre yoğunluğunda 96 kuyu tabakalarına tohumlanmıştır. Tohumlamadan 24, 48 ve 72 saat sonra, her kuyuya 20 μL MTT çözeltisi (fosfatla tamponlanmış tuzlu salinde 5 mg/mL) eklendi ve 37 °C'de 4 saat kuluçka edildi. Süpernatant çıkarıldı ve formazan kristallerini çözmek için 150 μL dimetil sülfoksit (DMSO) eklendi. Plaka tam çözünürme sağlamak için 10 dakika boyunca hafifçe çalkalandı. Absorbans 490 nm olarak mikroplaka spektrofotometresi kullanılarak ölçüldü. Boş kuyulardan alınan arka plan okumaları çıkarıldı. Hücre canlılığı 24 saat veya kontrol grubuna göre hesaplandı (1.0 olarak ayarlandı). Tüm deneyler altı teknik kopya ve üç bağımsız biyolojik replika ile gerçekleştirildi. Veriler ortalama ± standart sapma (SD) olarak ifade edilir ve istatistiksel anlamlılık iki kuyruklu t-testi kullanılarak belirlenmiştir.
Apoptoz için akış sitometrisi (Ek V - FITC/PI boyama)
Hücreler %60-70 birleşimde tohumlanarak belirtilen koşullarda 24 saat boyunca tedavi edildi. Yüzen ve yapışkan hücreler, EDTA-sız tripsin kullanılarak toplandı, birleştirildi ve buz gibi PBS ile iki kez yıkandı. Hücre pelletleri, 1 × 106 hücre/mL aralığında Annexin V bağlanma tamponunda (10 mM HEPES pH 7.4, 140 mM NaCl, 2.5 mM CaCl2) yeniden süspansiyona getirildi. Her örnek için, 100 μL süspansiyon, 5 μL Annexin V-FITC ve 5 μL propidium iyodid (PI; 50 μg/mL stok) ile oda sıcaklığında 15 dakika boyunca karanlıkta kuluçka altına alındı. 400 μL bağlanma tamponu eklendikten sonra, örnekler buzda tutuldu ve 1 saat içinde akış sitometresinde analiz edildi (488 nm uyarılma; FITC için 530/30 nm ve PI için >585 nm). Uygun tek boyalı ve floresans-eksi bir kontroller telaffuz için dahil edildi. Her örneklem için en az 10.000 olay kaydedildi. Veriler dörtgen kapıştırmalara göre analiz edildi: canlı (Annexin V⁻/PI⁻), erken apoptotik (Annexin V⁺/PI⁻), geç apoptotik (Annexin V⁺/PI⁺) ve nekrotik (Annexin V⁻/PI⁺) popülasyonlar. Erken + geç apoptotik hücrelerin yüzdeleri rapor edildi (ortalama ± SD, n = 3).
Nicel gerçek zamanlı PCR (qRT-PCR)
Toplam RNA, üreticinin protokolüne göre asit fenol-guanidinyum reaktifi kullanılarak izole edildi. RNA saflığı spektrofotometri ile doğrulandı (A₂₆₀/A₂₈₀ = 1.8-2.1), bütünlük ise jel elektroforezi ile doğrulandı (RNA bütünlük numarası ≥ 7). Toplam RNA'nın bir mikrogramı DNase I ile işlendi ve rastgele heksamerler ve oligo(dT) primerleri kullanılarak 20 μL reaksiyonunda ters transkribe edildi. Reaksiyon 25 °C'de 10 dakika, 50 °C'de 30 dakika ve 85 °C'de 5 dakika gerçekleştirildi. Nicel PCR, 5 μL 2× SYBR Green Master Mix, 0,3 μM primer ve 1 μL cDNA (≈ 20 ng RNA eşdeğeri) içeren 10 μL sistemde gerçekleştirildi. Termal döngü koşulları 5 dakika 95 °C, ardından 15 saniye için 95 °C ve 30 saniye 60 °C 40 döngü, ardından 0,3 °C aralıklarla 65 °C'den 95 °C'ye erime eğrisi analizi yapıldı. Tüm reaksiyonlar üçlü olarak, şablonsuz ve eksi-RT kontrolleri ile birlikte gerçekleştirildi. CT değerleri > 35 veya teknik kopya SD > 0.5 hariç tutuldu. Göreli ifade, GAPDH ile iç kontrol olarak 2⁻ΔΔCt yöntemiyle hesaplandı. Üç bağımsız biyolojik kopyadan ortalama ± SD değerleri raporlandı ve grup farklılıkları iki kuyruklu bir t-testi kullanılarak analiz edildi.
Western Blot analizi
Hücreler, proteaz ve fosfataz inhibitörleriyle takviye edilmiş RIPA tamponunda (50 mM Tris-HCl, pH 7,4, 150 mM NaCl, %1 NP-40, %0,5 sodyum deoksikolat, %0,1 SDS) üzerinde lizize edildi. Lizatlar aralıklı girdaplarla 30 dakika boyunca buz üzerinde kuluçka yapıldı ve 12.000 × g sıcaklıkta 4 °C'de 15 dakika boyunca santrifüjle temizlendi. Protein konsantrasyonları BCA testiyle ölçüldü, 1-2 μg/μL olarak ayarlandı ve 1:3 oranında 4× Laemmli tamponu (son 1× tamponu 100 mM DTT içerir) karıştırıldı. Örnekler 95 °C'de 5 dakika boyunca denatüre edildi. Eşit miktarda protein (50 μg) %100 V'da %12 SDS-PAGE ile 90 dakika çözüldü ve 250 mA'da PVDF membranlarına 90 dakika elektro olarak aktarıldı. Zarlar, oda sıcaklığında 1 saat boyunca TBST'de %5 yağsız süt (%0,1 Tween-20) ile bloklandı (fosfoproteinler için %5 BSA) ve IGF1, IGF1R, JAK2, p-JAK2 (Y1007/1008), STAT3, p-STAT3 (Y705), BIRC5 ve β-aktin (tipik seyreltme 1:1000, β-aktin 1:5000) karşısında primer antikorlarla gece boyunca 4 °C'de inkübe edildi. TBST'de üç 10 dakikalık yıkamadan sonra, zarlar HRP ile konjuge edilmiş ikincil antikor (1:5000) ile oda sıcaklığında 1 saat kuluçka altına alındı, tekrar yıkandı ve kemilyünesan substrat kullanılarak geliştirildi. Bant yoğunlukları ImageJ ile nicelendirildi, β-aktin veya toplam protein olarak normalize edildi ve üç bağımsız deneyden ortalama SD ± olarak ifade edildi.
İmmünositokimya
Steril cam kapak örtülerinde yetiştirilen hücreler PBS ile iki kez duruldu ve oda sıcaklığında 15 dakika boyunca %4 paraformaldehit ile sabitlendi. Üç PBS yıkamasından sonra, hücreler %0,2 Triton X-100 ile 10 dakika permeabilize edildi, %5 sığır serum albümini (BSA) ile 1 saat bloklandı ve gece boyunca 4 °C'de primer anti-CD31 antikoru ile inkubasyon edildi (%1 BSA'da 1:200 seyreltme). Üç PBS yıkamasının ardından, hücreler Alexa Fluor ile konjuge edilmiş ikincil antikor (1:500 seyreltme) ile 1 saat karanlıkta inkübe edildi, DAPI (1 μg/mL, 5 dakika) ile karşı boyandı ve antiifade ortamına monte edildi. Görüntüler, aynı pozlama ve kazanç ayarlarında floresan mikroskobu kullanılarak çekildi. CD31-pozitif alan yüzdesi, ImageJ yazılımı kullanılarak her örneklemde rastgele seçilmiş beş örtüşmeyen alanda nicelendirildi. Bu test, doku kesitleri yerine hücre modellerinde gerçekleştirildi.
İstatistiksel analiz
İstatistiksel analizler, R sürüm 4.3.0 ve ilgili paketleri kullanılarak gerçekleştirildi. Kategorik değişkenleri karşılaştırmak için ki-kare testi kullanıldı, sürekli değişkenler ise Wilcoxon sıra-toplamı testi veya T testi kullanılarak değerlendirildi. Sürekli değişkenlerin değerlendirilmesi Pearson korelasyon katsayısı ile gerçekleştirilmiştir. Hayatta kalma analizleri, Cox orantılı tehlike modellemesi ve Kaplan-Meier eğrilerinin oluşturulmasını içeren hayatta kalma paketi kullanılarak gerçekleştirildi; optimal tabakalaşma eşikleri survminer paketi ve Riskscore =
formülü ile belirlendi. CompareC paketi, çeşitli değişkenlerin C-indekslerini değerlendirmek için kullanıldı. İkili kategorik değişkenleri tahmin etmeyi amaçlayan alıcı çalışma karakteristik eğrisi (ROC), pROC paketi kullanılarak oluşturuldu. Ayrıca, hayatta kalma ölçütleri için ROC eğrisinin (AUC) altındaki zamana bağlı alan timeROC paketi kullanılarak analiz edildi. Tüm istatistiksel testler iki taraflı bir yaklaşımla gerçekleştirildi. P < 0.05 anlamlılık seviyesi istatistiksel olarak anlamlı olarak kabul edildi.