Çalışma, Helsinki Bildirgesi'ne uygun olarak yürütülmüştür. Protokol, Anhui Halk Sağlığı Klinik Merkezi Etik Komitesi tarafından 19 Eylül 2025'te onaylandı (onay kimliği: PJ-YX2025-062). Kan alımı öncesi tüm katılımcılardan yazılı bilgilendirilmiş onay alındı. Yerel kohort, akut miyokard enfarktüsü () olan sekiz hasta ve sekiz sağlıklı kontrol grubundan oluşuyordu. Protokolde kullanılan araştırma araçları Materyaller Tablosu'nda listelenmiştir.
1. Veri kaynakları ve işleme
Miyokard enfarktüsüyle ilgili toplu RNA dizileme veri setleri Gen İfadesi Omnibus'tan (GEO) alındı. GSE59867 ve GSE48060 veri setleri toplu transkriptomik analizler için kullanılmış, GSE269269 tek hücreli RNA dizileme veri seti ise hücresel düzeyde analizler için kullanılmıştır (Tablo 1). İnsan Protein Atlas'ından sonraki gen seti analizleri için 255 kalmodulinle ilişkili gen seti alındı.
| Veri seti | Örnek türü | Örnek (kontroller) | Örnek (hastalar) | Dizileme platformu |
| GSE59867 | Toplu RNA-dizisi | 46 | 111 | GPL6244 |
| GSE48060 | Toplu RNA-dizisi | 21 | 31 | GPL570 |
| GSE269269 | scRNA-seq (periferik kan) | | 10 | GPL24676 |
Tablo 1: Çalışmada kullanılan veri setlerinin özellikleri. Tablo, toplu RNA ve tek hücreli RNA dizileme veri setleri için veri seti erişim sayılarını, örnek türlerini, kontrol ve hasta örneklerinin sayılarını ve dizileme platformlarını listeler. RNA dizileri, RNA dizileme; scRNA-dizi, tek hücreli RNA dizilemesi.
Toplu transkriptomik veri setlerindeki örnek aralarındaki değişkenlik, limma paketindeki 3.60.6 sürümündeki normalizeBetweenArrays fonksiyonu kullanılarak düzeltildi. Daha sonra limma kullanılarak diferansiyel gen ekspresyonu analizi yapıldı. Farklı şekilde ifade edilen genler (DEG), P < 0.05 ve |log₂ katlama değişimi| > 0.5. Ortaya çıkan DEG'ler, volkan grafikleri ve ısı haritaları kullanılarak görselleştirildi ve önemli ölçüde yükseltilmiş, önemli ölçüde aşağı düzenlenmiş veya önemli ölçüde değişmemiş olarak sınıflandırıldı.
2. Tek örnek gen seti zenginleştirme analizi ve ağırlıklı gen ortak ifade ağı analizi
Tek örneklem gen seti zenginleştirme analizi (ssGSEA) 255 kalmodülinle ilişkili gen kullanılarak gerçekleştirildi. GSVA paketi, her ve kontrol örneği için kalmodülinle ilişkili gen skoru hesaplamak için kullanıldı; bu skor Calmodulin_score olarak belirlendi. ve kontrol grupları arasındaki Calmodulin_score farkları, Wilcoxon sıralama toplamı testi kullanılarak değerlendirildi.
Ağırlıklı gen ortak ekspresyon ağı analizi (WGCNA), hastalarından alınan toplu transkriptomik veriler kullanılarak gerçekleştirildi. Milyon haritalanmış okuma başına kilobaz başına ortalama parça değeri ≤0.5 olan genler hariç tutuldu. Örneklemeler istisnaları belirlemek ve çıkarmak için kümelendi.
R² > 0.8 ölçeksiz topoloji uyumu sağlayan yumuşak eşik güç seçildi. Daha sonra bir topolojik örtüşme matrisi oluşturuldu. Gen modülleri, dinamik ağaç kesimi algoritması kullanılarak minimum modül boyutu 200 ile tanımlandı. Oldukça benzer özgenlere sahip modüller, modül birleştirme eşiği 0.25'e karşılık gelen >0.75 korelasyon eşiği kullanılarak birleştirildi.
Modül özgenleri ile klinik özellikler arasındaki ilişkiler, Calmodulin_score dahil, Pearson korelasyon analiziyle değerlendirildi. Ortaya çıkan modül-özellik ilişkileri, korelasyon katsayıları ve karşılık gelen P değerleriyle açıklanan bir ısı haritasında gösterildi. Her gen için modül üyeliği ve gen anlamlılığı hesaplandı. Modül üyeliğinin gen anlamlılığına karşı saçılım grafikleri, yüksek intramodüler bağlantıya ve özellik alakalılığına sahip genleri belirlemek için oluşturuldu.
3. ile ilişkili kalmodülin ilişkili genlerin tanımlanması
ile ilişkili kalmodülin ile ilişkili genler, DEG'lerin WGCNA modüllerinden gelen genlerle kesişmesiyle tanımlandı ve Calmodulin_score ile anlamlı korelasyon sağlandı. Örtüşen genler, sonraki analizler için korundu.
Fonksiyonel zenginleştirme analizi clusterProfiler paketi kullanılarak gerçekleştirildi. Gen Ontolojisi ve Kyoto Genler ve Genomlar Ansiklopedisi açıklamaları, örtüşen genlerle ilişkili biyolojik süreçleri, moleküler fonksiyonları, hücresel bileşenleri ve sinyal yollarını tanımlamak için kullanıldı.
4. Makine öğrenimi kullanarak anahtar genlerin tanımlanması ve doğrulanması
Tek değişkenli lojistik regresyon analizi, örtüşen genler kullanılarak gerçekleştirildi. Daha sonra üç makine öğrenimi algoritması, aşağıdaki R paketleri ve parametreleri kullanılarak bağımsız olarak uygulandı: rastgele orman, XGBoost ve destek vektör makinesi.
Her algoritma, için öngörücü değere sahip genleri önceliklendirmek için kullanıldı. Aday anahtar genler, üç algoritmanın tamamı tarafından tanımlanan genler olarak tanımlandı. Hem GSE59867 eğitim veri setinde hem de GSE48060 dış doğrulama veri setinde anlamlı ve yönsel tutarlı diferansiyel ifade gösteren genler, nihai anahtar genler olarak korundu.
5. Tanı modellerinin inşası ve değerlendirilmesi
Tanımlanan anahtar genler ve lrm fonksiyonu kullanılarak lojistik bir regresyon modeli oluşturuldu. Regplot fonksiyonu kullanılarak her özellik geninin tahmin edilen olasılığına katkısını göstermek için bir nomogram oluşturuldu.
Model ayırtıcılığı, alıcının çalışma karakteristik analiziyle pROC paketi kullanılarak değerlendirildi. Alıcı çalışma karakteristik eğrisinin altındaki alan, modelin'yi kontrol örneklerinden ayırt etme yeteneğini değerlendirmek için hesaplandı.
Tahmin edilen olasılıkları gözlemlenen sonuçlarla karşılaştırmak için kalibrasyon eğrileri oluşturuldu. Karar eğrisi analizi, modelin net klinik faydasını çeşitli eşik olasılıklar aralığında tahmin etmek için gerçekleştirildi.
6. Gen seti zenginleştirme analizi ve rakip endojen RNA ağı yapısı
Gen seti zenginleştirme analizi, her anahtar gen için ayrı ayrı genler, gen-korelasyon matrisi ve clusterProfiler paketi kullanılarak gerçekleştirildi. Kyoto Genler ve Genomlar Ansiklopedisi yol zenginleştirme sonuçları, mutlak normalize zenginleştirme puanına göre sıralandı. Her gen için en yüksek beş yol gösterildi.
Fonksiyonel ilişkiler ve gen-etkileşim ağları GeneMANIA kullanılarak incelendi. Ana genlerin potansiyel mikroRNA düzenleyicileri miRanda, miRTarBase, TargetScan ve miRDB kullanılarak tahmin edilmiştir. Aday mikroRNA-mRNA etkileşimleri, dört veritabanındaki tahminlerin kesişmesiyle tanımlandı.
SpongeScan'den uzun kodlamayan RNA-mikroRNA etkileşimleri elde edildi. Uzun kodlamayan RNA-mikroRNA ve mikroRNA-mRNA ilişkileri, rakip bir endogen RNA düzenleyici ağ oluşturmak için entegre edildi. Ağ, ggalluvial paket kullanılarak bir Sankey diyagramı olarak görselleştirildi.
7. İlaç tahmini ve moleküler kenetlenme
İlaç-gen etkileşimleri, İlaç-Gen Etkileşim Veritabanı kullanılarak tahmin edilmiştir. Ortaya çıkan etkileşim ağı, ağ analiz yazılımı kullanılarak görselleştirildi.
CCL4 için UniProt protein tanımlayıcısı P13236 olarak alındı. İlgili üç boyutlu protein yapısı, kenetlenmek üzere seçilen 1HUM (insan MIP-1β, X-ışını kırınım yapısı) ile Protein Veri Bankası formatında (PDB) elde edildi. Biyolojik olarak ilgili monomeri temsil eden A zinciri kenetlenmek üzere seçildi. Protein hazırlığı, CB-Dock2'deki Prepare Protein modülü kullanılarak gerçekleştirildi; bu modül su moleküllerinin çıkarılmasını, polar hidrojenlerin eklenmesini ve Gasteiger yüklerinin atanmasını içerir. Aday bileşiklerin (klodronik asit ve epoetin alfa) üç boyutlu kimyasal yapıları, PubChem veritabanından Yapı-Veri Dosyası (SDF) formatında alındı. Kenetlenme simülasyonları, kör kenetlenme için AutoDock Vina algoritmasını kullanan CB-Dock2 çevrimiçi platformu kullanılarak gerçekleştirildi. Kenetlenme alanı, potansiyel bağlanma ceplerinin tarafsız tanımlanmasını sağlamak için tüm protein yüzeyini kapsayacak şekilde ayarlandı. Bağlanma afinitesi, kcal/mol cinsinden tahmin edilen serbest bağlanma enerjisi (ΔG) olarak hesaplandı. Son kenetlenme pozları ve protein-ligand etkileşimleri (örneğin, hidrojen bağları, hidrofobik kontaklar) PyMOL ve CB-Dock2'nin yerleşik etkileşim görüntüleyicisi kullanılarak görselleştirildi.
8. Tek hücreli RNA dizileme veri ön işlemesi
Kalite kontrol, aşağı akış tek hücreli RNA dizileme analizinden önce gerçekleştirildi. Hücreler, tespit edilen gen sayısı 200 ile 10.000 arasında olduğunda, toplam benzersiz moleküler tanımlayıcı sayısı ≥1.000 ve mitokondriyal transkript oranı %≤20 olduğunda tutuldu.
200'den az gen ifade eden hücreler ve üçten az hücrede tespit edilen genler dışlandı. Bu filtreler, düşük kaliteli hücrelerin ve teknik gürültünün dahil edilmesini azaltmak için uygulandı. Gen ifade değerleri, Seurat paketindeki NormalizeData fonksiyonu kullanılarak normalleştirildi. FindVariableFeatures fonksiyonu kullanılarak yüksek değişkenlikli genler tanımlandı. Yüksek değişkenlikli genlerin ifade değerleri ScaleData fonksiyonu kullanılarak merkezlenip standartlaştırıldı.
Deneysel veya dizileme varyasyonuyla ilişkili parti etkileri, Harmony entegrasyonçerçevesi 17'den RunHarmony fonksiyonu kullanılarak düzeltildi.
9. Tek hücre boyutluluk azaltma, kümeleme ve açıklama
Ana bileşen analizi ilk olarak tek hücreli RNA dizileme veri setinin boyutluluğunu azaltmak için uygulanmıştır. Uniform manifold yaklaşımı ve projeksiyonu ile t-dağılımlı stokastik komşu gömülü, hücresel heterojenliği görselleştirmek için daha sonra kullanılmıştır.
Transkripsiyonel olarak benzer hücreler, Seurat'taki FindNeighbors ve FindClusters fonksiyonları kullanılarak gruplandırıldı. Her küme için farklı şekilde ifade edilen belirteç genleri, FindAllMarkers fonksiyonu kullanılarak her küme kalan tüm kümelerle karşılaştırıldığında tanımlandı.
Hücre tipleri, yayımlanmış literatürden ve kurulmuş hücre-belirteç veri tabanlarından alınan kanonik belirteç genleri kullanılarak atanmıştır. Ana genlerin mekansal dağılımı ve ifade seviyeleri, FeaturePlotfonksiyonu 18 kullanılarak görselleştirildi.
10. Nicel polimeraz zincir reaksiyon analizi
Anhui Halk Sağlığı Klinik Merkezi'nde 8 hastası ve 8 sağlıklı kontrolden periferik kan örnekleri alındı. grubu, Miyokard İnfarktüsü'nün Dördüncü Evrensel Tanımı'na göre teşhis konmuş hastaları içeriyordu; semptomları miyokard iskemisi ile uyumlu ve 99. perdelik üst referans sınırının üzerinde yükselen kardiyak troponin I seviyeleri bulunuyordu. Kontrol grubu, kardiyovasküler hastalık öyküsü olmayan, normal elektrokardiyogram ve rutin kan testlerinde, karaciğer fonksiyonunda veya böbrek fonksiyonunda anormallik olmayan yaş ve cinsiyet uyumlu sağlıklı bireylerden oluşuyordu. hastaları için, hastaneye yatıştan sonraki 24 saat içinde 3 mL etilen diamintetraasetik asit antikoagüllenmiş kan toplandı. Aynı hacim, ilgili çalışma döneminde sağlıklı kontrol kontrollerinden de toplanmıştır.
Kan RNA izolasyon kiti ile birlikte verilen protokole göre periferik kandan toplam RNA izole edildi. RNA konsantrasyonu ve saflığı NanoDrop spektrofotometresi kullanılarak değerlendirildi ve RNA bütünlüğü agaroz jel elektroforezi ile doğrulandı. Sonraki analizler için yalnızca A260/A280 oranı 1.8 ile 2.1 arasında olan örnekler kullanıldı. Toplamda 500 ng RNA, birinci iplik tamamlayıcı DNA sentezi reaktifi kullanılarak tamamlayıcı DNA'ya ters-transkribe edildi. Ortaya çıkan tamamlayıcı DNA, son konsantrasyon seviyesi 150 ng/mL'ye kadar seyreltilmiştir. Nicel polimeraz zincir reaksiyon amplifikasyonu, pasif referans boyası olmadan SYBR Yeşil bazlı ana karışım kullanılarak toplam 10 μL reaksiyon hacminde gerçekleştirilmiştir. Tüm qPCR reaksiyonları teknik kopyalar halinde iki kez gerçekleştirildi ve sonraki hesaplamalar ortalama CT değerlerine dayandı.
Amplifikasyon, gerçek zamanlı polimeraz zincir reaksiyon cihazı kullanılarak gerçekleştirildi. Döngü koşulları, ilk denatürasyondan oluşuyordu (95 °C'de 5 dakika), ardından 40 döngü denaturasyon (95 °C) sıcaklıkta 10 saniye, 60 °C'de 30 saniye boyunca tavlama ve 30 saniye boyunca 72 °C'de uzatma yapıldı. Erik eğrisi analizi amplifikasyondan sonra yapılmıştır.
Gen ekspresyonu seviyeleri β-aktine normalleştirildi. Göreli ifade 2−ΔΔCt yöntemiyle hesaplandı.
11. İstatistiksel analiz
İstatistiksel analizler R'de yapıldı. Ağ görselleştirmeleri ağ analiz yazılımı kullanılarak oluşturuldu. İki grup arasındaki farklar, aksi belirtilmedikçe Wilcoxon testi kullanılarak değerlendirildi. Normal dağılıma sahip sürekli değişkenler Öğrencinin t-testi kullanılarak karşılaştırıldı. Normal dağılım olmayan sürekli değişkenler, Mann-Whitney U testi ile karşılaştırıldı; bu test Wilcoxon sıra-toplamı testi olarak da bilinir. Tüm istatistiksel testler iki kuyrukluydu. <0,05 P değeri istatistiksel olarakanlamlı olarak 19 olarak kabul edildi.