Veri toplama
DN ile ilişkili olan GSE142025 ve GSE96804 veri setleri için gen ekspresyon matrisi ve ilgili klinik veriler, Gene Expression Omnibus (GEO) veri tabanından (https://www.ncbi.nlm.nih.gov/geo/) alındı15. Eğitim seti (GSE142025), GPL20301 platformu kullanılarak sekanslanan 27 DN'li hasta ve dokuz kontrolden alınan böbrek dokusu örneklerini içermekteydi. Doğrulama seti (GSE96804), GPL17586 platformu kullanılarak işlenen 41 DN'li hasta ve 20 kontrole ait sekanslama verilerinden oluşmaktaydı. Her iki veri seti de böbrek dokusuna odaklanmaktadır; GSE96804 veri seti özellikle böbreğin birincil filtrasyon birimi olan glomerulü incelemektedir (Şekil 1). GSE142025 veri seti (eğitim seti), tüm böbrek dokusu örneklerinden oluşmakta ve DN'nin transkriptomik manzarasına dair kapsamlı bir genel bakış sunmaktadır. Diğer yandan GSE96804 veri seti (doğrulama seti), özellikle glomerüler filtrasyon hasarının birincil bölgesi olan glomerüler dokuya odaklanmaktadır. Platform ve doku farklılıkları nedeniyle bu iki veri seti doğrudan birleştirilmediğinden, grup etkisi (batch effect) düzeltmesi uygulanmamıştır. Bunun yerine, veri setleri arası doğrulama bağımsız olarak gerçekleştirilmiştir. Her iki veri setinde tutarlı yönlü değişimler gösteren ve eğri altında kalan alan (AUC) değeri 0,7'den büyük olan genler, farklı böbrek kompartmanları arasındaki genelleştirilebilirliği destekleyen güçlü adaylar olarak seçilmiştir.
MitoCarta3.0 veri tabanından (https://www.broadinstitute.org/mitocarta) toplam 1.136 mitokondri ile ilişkili gen (MRG) çıkarılmıştır. LLF'nin aktif bileşenleri, ≥30% oral biyoyararlanım (OB) eşiği ve ≥0.18 ilaç benzerliği (DL) eşiği temel alınarak Geleneksel Çin Tıbbı Sistem Farmakolojisi (TCMSP) veri tabanı (http://sm.nwsuaf.edu.cn/lsp/tcmsp.php) kullanılarak öngörülmüştür. Ardından, aktif bileşenler için potansiyel hedef genler Swiss Target Prediction veri tabanı (http://www.swisstargetprediction.ch/) kullanılarak öngörülmüştür.
Diferansiyel ekspresyon analizi
GSE142025'in (DN ile kontrol grubu) diferansiyel ekspresyon analizi, anlamlılık kriterleri P.adj < 0.05 ve |log2FoldChange| > 0.5 olarak belirlenerek limma paketi (v3.54.1) kullanılarak gerçekleştirilmiştir16. Volcano grafikleri ve ısı haritaları, sırasıyla ggplot2 (v 3.3.6) ve ComplexHeatmap (v 2.14.0) paketleri kullanılarak görselleştirilmiştir17,18. DEG'ler, MRG'ler ve aktif bileşenlerin potansiyel hedef genleri kesiştirilmiş ve örtüşen genler aday genler olarak tanımlanmıştır. Aktif bileşenleri ve aday genleri birbirine bağlayan ağ, Cytoscape yazılımı (v 3.9.0) kullanılarak oluşturulmuştur19.
Fonksiyonel zenginleştirme analizi ve protein-protein etkileşim (PPI) ağı oluşturma
Aday genlerin biyolojik fonksiyonlarını ve ilişkili sinyal yollarını araştırmak amacıyla, clusterProfiler paketi (v 4.6.2) kullanılarak Gen Ontolojisi (GO) ve Kyoto Gen ve Genom Ansiklopedisi (KEGG) zenginleştirme analizleri gerçekleştirilmiştir (P.adjust < 0.05). Ardından, PPI ilişkilerini belirlemek için aday genler STRING veri tabanına (https://cn.string-db.org/) girilmiş (güven düzeyi ≥ 0.4) ve Cytoscape (v 3.9.0) kullanılarak bir PPI ağı oluşturulmuştur20.
Makine öğrenimi
GSE142025 veri setine dayalı olarak; rastgele orman (RF), k-en yakın komşu (KNN), kısmi en küçük kareler (PLS) ve radyal tabanlı çekirdekli destek vektör makinesi (SVM) dahil olmak üzere dört makine öğrenimi algoritması, caret paketi (v6.0-93) kullanılarak uygulanmıştır. Önceki analizde tanımlanan aday genler tahmin değişkenleri olarak, hastalık durumu (DN veya kontrol) ise sonuç değişkeni olarak kullanılmıştır. KNN modeli için, tuneLength = 10 olacak şekilde trainControl fonksiyonu kullanılarak 10 katlı çapraz doğrulama uygulanmıştır. RF modeli 20 ağaçla (ntree = 20) kurulmuştur; PLS ve SVM modelleri ise orijinal caret iş akışında uygulanan ayarlar kullanılarak kurulmuştur.
Model artık dağılımları DALEX paketi (v2.4.3) kullanılarak değerlendirildi. Ardından, kayıp fonksiyonu olarak ortalama karekök hata (RMSE) belirlenerek DALEX variable_importance fonksiyonu ile permütasyon tabanlı değişken değerlendirmesi gerçekleştirildi. Elde edilen dropout_loss, ilgili değişkenin permütasyonundan sonra elde edilen RMSE değerini temsil eder21,2. Orijinal analitik iş akışında, tüm modellerde dropout kaybı < 0.281 olan değişkenler tutuldu ve dört modelin tamamında ortak olan genler, sonraki doğrulama aşamaları için konsensüs özellik genleri olarak tanımlandı.
Dört makine öğrenimi modeli, nihai bir klinik sınıflandırıcı oluşturmaktan ziyade öncelikle özellik önceliklendirmesi için kullanılmıştır. Buna bağlı olarak, tanısal ayrım daha sonra hem keşif hem de doğrulama veri setlerinde ROC analizi kullanılarak bireysel gen düzeyinde değerlendirilmiştir.
Aday biyobelirteçlerin belirlenmesi
GSE142025 ve GSE96804'te, DN ve kontrol örnekleri arasındaki özellik genlerinin ekspresyon farkları Wilcoxon testi kullanılarak değerlendirilmiştir. Sadece anlamlı diferansiyel ekspresyona (P < 0.05) sahip olan ve her iki veri setinde uyumlu ekspresyon trendleri gösteren genler, alıcı işletim karakteristiği (ROC) eğrisi analizi için seçilmiştir. ROC eğrilerinin oluşturulması ve eğri altındaki alanın (AUC) hesaplanması için pROC paketi (v1.18.0) kullanılmış olup, her iki veri setinde AUC > 0.7 değerini gösteren genler aday biyobelirteçler olarak sınıflandırılmıştır23.
Gen seti zenginleştirme analizi (GSEA)
Aday biyobelirteçlerle ilişkili biyolojik fonksiyonlar ve sinyal yolları, GSE142025 veri seti üzerinde GSEA kullanılarak daha ayrıntılı şekilde incelenmiştir. İlk olarak, aday biyobelirteçlerin diğer tüm genlerle olan Spearman korelasyon analizi psych paketi (v2.2.9) kullanılarak gerçekleştirilmiştir.24Korelasyon katsayıları hesaplanmış ve (yüksekten düşüğe doğru) sıralanmıştır. Referans gen seti olarak Molecular Signatures Database (MSigDB, https://www.gsea-msigdb.org/gsea/msigdb/) üzerinden c2.cp.kegg.v2023.1.Hs.symbols.gmt kullanılmıştır. Ardından, clusterProfiler paketinin (v4.6.2) kullanımıyla, sıralanan genlerin arka plan gen setindeki zenginleşmesini değerlendirmek amacıyla GSEA gerçekleştirilmiştir. Çoklu test düzeltmesi uygulanmıştır. aracılığıyla FDR yöntemi ve düzeltilmiş P-değerler (şöyle belirtilmiştir: P.adjust) değerleri, şu durumda anlamlı kabul edildi: < 0.05.
m6A modifikasyon analizi
Aday biyobelirteçlerin RNA metilasyon modifikasyonlarını incelemek için, m'yi öngörmek amacıyla SRAMP veri tabanı (http://www.cuilab.cn/sramp/) kullanılmıştır.6Aday biyobelirteçlerin ikincil yapıları içindeki yüksek güvenilirlikli pozisyonlara odaklanılarak modifikasyon bölgeleri belirlenmiştir. Ardından, |HepG2 (shRNA)| kullanılarak aday biyobelirteçlerle etkileşime giren m6A-modifiye proteinlerin tanımlanması için ENCORI veri tabanından (https://starbase.sysu.edu.cn/) yararlanılmıştır. > Kilit proteinleri taramak için 1 parametre kullanıldı. Ardından, kilit proteinler ile aday biyobelirteçler arasındaki etkileşim olasılığını tahmin etmek için RPISeq veri tabanı (http://pridb.gdcb.iastate.edu/RPISeq/) kullanıldı. RF ve SVM sınıflandırıcı tahmin skorlarını oluşturmak için her ikisinin RNA dizileri düz metin formatında yüklendi. Skor 0,5'i geçtiğinde etkileşim anlamlı kabul edildi.25SRAMP analizi, m6A bölgesi tahmini için ‘High’ eşiği, ‘Transcript’ modu ve varsayılan parametreler kullanılarak gerçekleştirilmiştir. ENCORI analizi, ‘HepG2 (shRNA)’ parametresi ile ‘miRNA-mRNA’ etkileşim fonksiyonu kullanılarak yapılmıştır. > 1. RPISeq analizi, varsayılan parametrelere sahip RF sınıflandırıcı kullanılarak gerçekleştirilmiştir; skorlar > 0,5 pozitif etkileşime işaret etmektedir. Bunlar hesaplamalı öngörülerdir; böbrek dokusundaki m6A modifikasyonunun veya RNA-protein etkileşimlerinin deneysel kanıtları değildir. HepG2 shRNA kriteri, ENCORI'nin önceden hesaplanmış veri setlerinden türetilmiştir ve böbreğe özgü düzenlemeyi yansıtmayabilir.
İmmün infiltrasyon analizi
GSE142025'teki hem kontrol hem de DN örneklerinde 2 immün hücre tipinin oranlarını tahmin etmek için CIBERSORT algoritması uygulanmış ve sonuçlar ggplot2 paketi (v3.3.6) kullanılarak oluşturulan bir ısı haritası ile görselleştirilmiştir26. CIBERSORT, LM2 imza matrisi kullanılarak, 1.0 permutasyon ile ve (mikrodiziler için önerildiği üzere) kuantil normalizasyonu devre dışı bırakılarak çalıştırılmıştır. CIBERSORT p < 0.05 olan örnekler ileri analizler için saklanmıştır. CIBERSORT, immün hücre fraksiyonlarını dökme (bulk) böbrek dokusu ekspresyonundan tahmin eder; bu yöntem kompartman spesifik infiltrasyonu (örneğin, glomerüler ile tubulointerstisyel arası) çözemez ve infiltre olan lökositleri resident immün hücrelerden ayıramaz. Bu nedenle, bildirilen korelasyonlar doku düzeyindedir ve uzamsal yöntemlerle doğrulanmalıdır. Diferansiyel bolluğa sahip immün hücreler ile aday biyobelirteçler arasındaki Spearman korelasyon analizi psych paketi kullanılarak gerçekleştirilmiştir.
Ağların oluşturulması ve moleküler kenetlenme
Aday biyobelirteçlerle etkileşime giren mikroRNA'lar (miRNA'lar), miRNet veri tabanı (https://www.mirnet.ca) kullanılarak öngörülmüştür. Ardından, belirlenen miRNA'ları hedefleyen uzun kodlamayan RNA'lar (lncRNA'lar), TarBase (http://www.diana.pcbi.upenn.edu/tarbase) ve starbase (http://starbase.sysu.edu.cn/) veri tabanları aracılığıyla öngörülmüştür. Ağ oluşturma için her iki veri tabanında ortak olan lncRNA'lar seçilmiştir. Daha sonra Cytoscape yazılımı kullanılarak bir lncRNA-miRNA-mRNA düzenleyici ağı kurulmuştur. Aday biyobelirteçleri hedefleyen potansiyel aktif bileşenler, bir aktif bileşen-biyobelirteç ağı oluşturmak üzere seçilmiştir. Ek olarak, aktif bileşenler, aday biyobelirteçler ve GSEA'da belirlenen yolaklar, bir aktif bileşen-biyobelirteç-yolak ağı oluşturmak amacıyla Cytoscape'e dahil edilmiştir.
Potansiyel aktif bileşenler ile aday biyobelirteçler arasındaki bağlanma afinitesini değerlendirmek için moleküler docking analizi gerçekleştirilmiştir. Biyobelirteç proteinlerinin 3D yapıları, PDB dosya formatında Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB, https://www.rscb.org/pdb) üzerinden temin edilmiştir. Potansiyel aktif bileşenlerin 2D yapıları, PubChem veri tabanından (http://pubchem.ncbi.nlm.nih.gov) SDF formatında alınmıştır. Moleküler docking, CB-Dock platformu (http://clab.labshare.cn/cb-dock/php/blinddock.php) kullanılarak yürütülmüştür. -5 kcal/mol'den düşük bir bağlanma enerjisi, güçlü bir bağlanma afinitesine işaret etmiştir27.
Ligustri Lucidi Fructus hazırlama ve kimliklendirme
Burada, Ligustri Lucidi Fructus (LLF), Ligustrum lucidum W. T. Aiton (Oleaceae) bitkisinin kurutulmuş olgun meyvelerini ifade etmektedir. Botanik materyal, Çin Farmakopesine göre kimliklendirilmiş olup, 20240506, 20240911, 20241103 numaralı voucher örnekleri Shanxi Geleneksel Çin Tıbbı Üniversitesi'ne teslim edilmiştir.
Dekoksiyon hazırlığı için, nitelikli LLF dilimlerinin 20 g'ı, oda sıcaklığında 30 dk boyunca 1.0 mL distile suda bekletildi. Karışım şiddetli bir şekilde kaynatıldı ve ardından 60 dk boyunca kısık ateşte kaynatılmaya devam edildi. Filtrat toplandı ve geride kalan bitki materyalleri, başka bir 1.0 mL distile su ile 60 dk boyunca tekrar dekoksiyona tabi tutuldu. İki filtrat birleştirildi, filtrelendi, santrifüj edildi ve indirgenmiş basınç altında konsantre edilerek 1 g ham ilaç/mL nihai stok konsantrasyonu (toplam hacim 10 mL) elde edildi. Hazırlanan dekoksiyon, kısa süreli kullanım için 4 °C'de veya uzun süreli saklama için −20 °C'de muhafaza edildi. Deneysel güvenilirliği ve tekrarlanabilirliği sağlamak amacıyla, LLF'nin ve dekoksiyonunun kalitesi Çin Farmakopesi standartlarına göre titizlikle tanımlandı ve doğrulandı.
Kalitatif tanımlama için ince tabaka kromatografisi uygulanmıştır. Kısaca, hazırlanan dekoksiyonun uygun bir hacmi santrifüjlenmiş ve süpernatan metanol ile ekstrakte edilmiştir. Filtrasyonun ardından, numune çözeltisi ve specnuezhenide referans standart çözeltisi aynı silika jel G plakasına damlatılmıştır. Geliştirme, kurutma ve ultraviyole ışık altında inceleme işlemlerinden sonra, numune çözeltisinin lekesi referans bileşik ile tutarlı renk ve konum göstermiş, böylece LLF'nin karakteristik biyoaktif bileşeninin varlığı doğrulanmıştır.
Kantitatif kalite kontrolü için yüksek performanslı sıvı kromatografisi deteksiyonu gerçekleştirilmiştir. Analiz, mobil faz olarak metanol–su kullanılan bir C18 kolonu ile yapılmıştır. Deteksiyon dalga boyu 24 nm olarak ayarlanmıştır. Dekoksiyondaki specnuezhenide içeriği, standart eğriye dayanarak belirlenmiştir. Sonuçlar, hazırlanan dekoksiyonun kimyasal bileşiminin stabil ve homojen olduğunu göstermiş ve hayvan müdahale deneyi boyunca tutarlı ilaç kalitesi sağlamıştır.
Hayvan deneyleri
On iki adet SPF sınıfı erkek db/db fare (8-9 haftalık) ve yaşla uyumlu altı adet db/m fare, Shanxi Geleneksel Çin Tıbbı Üniversitesi'nin SPF hayvan tesisinde tutulmuştur. Deneylerden önce hayvanlar, 12 saat ışık/12 saat karanlık döngüsü altında 7 gün boyunca adapte edilmiştir. ad libitum gıda ve suya erişim. Çalışma, Shanxi Geleneksel Çin Tıbbı Üniversitesi Etik Kurulu tarafından onaylanmıştır (Onay No. 2022DW167). Belirtiler gösteren hayvanlar >Vücut ağırlığının %20'sini kaybeden, moribund duruma gelen veya yiyecek ve suya erişemeyen fareler, planlanan çalışma bitiş noktasından önce insani yöntemlerle ötanize edilmiştir. Çalışmanın sonunda, geri kalan tüm fareler, kurum protokollerine uygun olarak intraperitonal pentobarbital sodyum enjeksiyonu ve ardından servikal dislokasyon ile ötanize edilmiştir.
Alıştırma dönemini takiben, db/db farelerde DN modelinin oluşturulduğu, kuyruk veni kan glikoz seviyesinin ≥ 16 mmol/L olması ve pozitif idrar mikroalbumin test şeridi ile belirtilen mikroalbüminüri ile doğrulanmıştır. DN modelinin başarıyla oluşturulmasının ardından, db/db fareler rastgele iki gruba (grup başına n = 6) ayrılmıştır: DN model grubu (DN) ve LLF tedavi grubu (Treatment). Ayrıca, kontrol grubu (Control) olarak db/m fareler (n = 6) kullanılmıştır. Dozaj, diyabetik sıçanlarda LLF ile yapılan önceki farmakodinamik çalışmalara dayanarak seçilmiş ve vücut yüzey alanı normalizasyonu kullanılarak insan eşdeğer dozuna dönüştürülmüştür28. Kontrol ve DN gruplarına distile su verilirken, tedavi grubu 8 hafta boyunca 3.5 g/kg LLF almıştır. 8 haftalık uygulamadan sonra, sonraki testler için serum, idrar ve böbrek dokusu toplamak amacıyla tüm fareler ötanize edilmiştir.
Kan ve idrar göstergeleri
Serum glukoz seviyeleri tam otomatik bir Kan Biyokimya Analizörü kullanılarak analiz edildi. İdrardaki mikroalbümin konsantrasyonu kit talimatlarına göre ölçüldü (Ek Dosya 1).
Fare böbrek dokularının patolojik gözlemi
Böbrek dokuları histopatolojik inceleme için işlendi. %4 paraformaldehit ile fikse edildikten sonra dokular yıkandı,yıkandı,dehidre edildi,parafine gömüldü ve kesitleri alındı. Ardından hematoksilen ve eozin (HE) boyaması yapıldı ve patolojik değişiklikleri değerlendirmek için boyalı kesitler optik mikroskop altında incelendi.
Ters transkripsiyon kantitatif polimeraz zincir reaksiyonu (RT-qPCR)
Fare böbrek dokusundaki aday biyobelirteçlerin ekspresyonunu belirlemek için RT-qPCR kullanılmıştır. Toplam RNA, üretici talimatlarına göre ekstrakte edilmiş, ardından RNA konsantrasyonu ve kalitesi değerlendirilmiştir (Tablo 1). Ekstrakt RNA'dan, cDNA Sentez Kiti kullanılarak cDNA sentezlenmiştir. Amplifikasyon, şurada listelenen primer çiftleri ile gerçekleştirilmiştir: Tablo 1, kullanarak GAPDH referans gen olarak kullanılmıştır. Göreceli ekspresyon seviyeleri 2 ile hesaplanmıştır−ΔΔCt yöntem13,26.
İstatistiksel analiz
Tüm istatistiksel analizler, ilgili deneysel prosedürler için gerekli yazılımlarla birlikte R yazılımı (versiyon 4.2.2) kullanılarak gerçekleştirilmiştir. Aksi belirtilmedikçe, tüm istatistiksel testler çift yönlü yapılmış ve farklar P < 0.05 değerinde istatistiksel olarak anlamlı kabul edilmiştir. Transkriptomik diferansiyel ekspresyon analizi limma paketi kullanılarak gerçekleştirilmiştir. Ayarlanmış P < 0.05 ve 0.5'ten büyük mutlak log2 kat değişim değerine sahip genler, diferansiyel eksprese olmuş olarak tanımlanmıştır.
Bağımsız DN ve kontrol örnekleri arasındaki aday gen ekspresyonlarının karşılaştırılması için, orijinal analitik iş akışında belirtildiği yerlerde parametrik olmayan Wilcoxon sıra toplamı testleri kullanılmıştır. Aday biyobelirteçler ile bağışıklık hücresi fraksiyonları arasındaki korelasyonlar, Spearman sıra korelasyonu kullanılarak değerlendirilmiştir.
Deneysel veriler ortalama ± SD olarak sunulmuştur. Üç bağımsız hayvan grubu arasındaki karşılaştırmalar, parametrik analiz varsayımları karşılandığında tek yönlü varyans analizi kullanılarak gerçekleştirilmiştir. Post hoc karşılaştırmalar için, varyanslar homojen olduğunda en düşük anlamlı fark testi, varyanslar eşit olmadığında ise Dunnett's T3 testi kullanılmıştır. Bazal ve 8. hafta farmakodinamik ölçümleri ayrı analiz edilmiş ve sunulmuştur; grup-zaman etkileşimine ilişkin herhangi bir çıkarım yapılmamıştır. Biyoenformatik ve makine öğrenimi analizleri için kullanılan R betikleri ve kaynak veriler Ek Dosya 1'de sunulmuştur.