DN tedavisinde LLF için potansiyel mitokondri ile ilişkili aday biyobelirteçleri sistematik olarak araştırmak amacıyla dört aşamalı bir analitik iş akışı tasarladık (Şekil 1). I. Aşama'da, GSE142025 veri setinden (eğitim seti, tüm böbrek, n=36) ve GSE96804 veri setinden (doğrulama seti, glomerulus, n = 61) alınan transkriptomik verileri, MitoCarta 3.0 veri tabanındaki 1.136 mitokondri ile ilişkili gen ve TCMSP veri tabanındaki 9 aktif bileşenin 517 öngörülen hedefi ile entegre ettik. Bu üç gen setinin kesişimi sonucunda 9 aday gen elde edildi. II. Aşama'da, RMSE < 0,281 eşik değeri kullanılarak özellik genlerini önceliklendirmek için dört makine öğrenimi modeli (RF, KNN, PLS ve SVM) uygulandı. ROC analizi (her iki veri setinde AUC > 0,7) kullanılarak yapılan veri setleri arası doğrulama ile dört aday biyobelirteç belirlendi: CAT, FABP1, MAOB ve MAOA. III. Aşama'da, zenginleşmiş KEGG yolaklarını belirlemek için GSEA, CIBERSORT kullanarak immün infiltrasyon analizi ve m6A modifikasyon öngörüsü gerçekleştirdik; ardından lncRNA-miRNA-mRNA, aktif bileşen-biyobelirteç ve aktif bileşen-biyobelirteç-yolak ağları oluşturup moleküler kenetleme yaptık. IV. Aşama'da, LLF'nin farmakodinamik etkileri ve dört aday biyobelirtecin mRNA ekspresyonundaki değişiklikler, DN'nin bir db/db fare modelinde değerlendirildi.
DN tedavisi için LLF aday genlerinin taranması
GSE142025 veri setinde, DN ve kontrol grupları arasında 1.904'ü artmış ve 1.906'sı azalmış olmak üzere toplam 3.810 DEG belirlenmiştir (Şekil 2A,B). TCMSP veri tabanı kullanılarak LLF'den on üç aktif bileşen öngörülmüştür; bunlar beta-sitosterol, kaempferol, taxifolin, Lucidumoside D, Lucidumoside D_qt, (20S)-24-ene-3,20-diol-3-acetate, eriodictyol, syringaresinol diglucoside_qt, Lucidusculine, Olitoriside, Olitoriside_qt, luteolin ve quercetin'dir (Tablo 2). Dört aktif bileşen—Lucidumoside D_qt, (20S)-24-ene-3,20-diol-3-acetate, syringaresinol diglucoside_qt ve Olitoriside_qt—herhangi bir potansiyel hedef gen öngörmemiş, geri kalan dokuz bileşen ise 517 potansiyel hedef gen öngörmüştür. 3.810 DEG, 1.136 MRG ve 517 potansiyel hedef genin çakıştırılmasıyla dokuz aday gen belirlenmiştir: GPX1, BAX, CASP8, MAOA, MAOB, CAT, AKR1B10, ALDH2 ve FABP1 (Şekil 2C). Ardından bir aktif bileşen-aday gen ağı oluşturulmuştur (Şekil 2D). Bu dokuz aday gen; toksik maddelere yanıt, organik hidroksi bileşik katabolik süreci ve hücresel detoksifikasyon dahil olmak üzere 341 GO teriminde zenginleşmiştir (Şekil 2E). Ayrıca, triptofan metabolizması, nörodejeneratif yolaklar ve histidin metabolizması gibi 52 KEGG yolağı ile ilişkilendirilmişlerdir (Şekil 2F).
LLF'de DN tedavisi için aday biyobelirteçlerin taranması
PPI ağı, MAOA, ALDH2, MAOB ve AKR1B10'in etkileşime girdiği yedi düğüm ve sekiz kenar ortaya çıkarmıştır (Şekil 3A). Dört makine öğrenimi modeli genelinde RMSE değerleri 0,281'den küçük olan genler, özellik genleri olarak tanımlanmıştır: CAT, MAOB, MAOA, BAX ve FABP1 (Şekil 3B-E). Ekspresyon analizi, CAT, FABP1, MAOB ve MAOA'nın DN ve kontrol grupları arasında anlamlı derecede farklı olduğunu ve hem GSE142025 hem de GSE96804 veri setlerinde tutarlı olduğunu göstermiştir (Şekil 3F,G). Ayrıca, ROC eğrisi analizindeki AUC değerleri her iki veri setinde de 0,7'yi aşmış olup, bu durum söz konusu genlerin DN örneklerini kontrol örneklerinden etkili bir şekilde ayırt edebileceğini ve LLF'de DN tedavisi için aday biyobelirteçler olarak hizmet edebileceğini göstermektedir (Şekil 4A-H).
İnflamatuvar ve immün ilişkili yolaklarda aday biyobelirteçlerin anlamlı zenginleşmesi
GSEA, kemokin sinyal yolu ve sitokin-sitokin reseptör etkileşimlerinde belirgin şekilde zenginleşmiş dört aday biyobelirteç tanımlamıştır (Şekil 5A-D). Bunlar arasında peroksidaz sinyal yolu; CAT, MAOA ve MAOB ile anlamlı bir ilişki göstermiştir.
Aday biyobelirteçlerin bağışıklık hücreleri ile korelasyonu
DN ve kontrol örnekleri arasında dokuz bağışıklık hücresi tipinin—naif B hücreleri, M0 Makrofaj, M1 Makrofaj, M2 Makrofaj, aktive mast hücresi, aktive NK hücresi, dinlenme halindeki hafıza CD4+ T hücresi, naif CD4+ T hücresi ve CD8+ T hücresi—ekspresyonunda belirgin farklılıklar gözlemlenmiştir (P < 0.05) (Şekil 6A,B). Naif B hücreleri ile aktive NK hücreleri arasında anlamlı pozitif bir korelasyon (cor = 0.6) bulunurken, naif B hücreleri ile aktive mast hücreleri arasında anlamlı negatif bir korelasyon (cor = -0.69) saptanmıştır (Şekil 6C). Tüm aday biyobelirteçler, CD8+ T hücreleri ve aktive mast hücreleri ile güçlü negatif korelasyonlar; aktive NK hücreleri ve naif B hücreleri ile ise pozitif korelasyonlar sergilemiştir (Şekil 6D).
Temel modifiye m6A proteinlerinin aday biyobelirteçlerle etkileşimi
The m6Bir RNA metilasyon modifikasyonu, RNA sentezini ve metabolizmasını derinlemesine etkiler ve çeşitli hastalıkların patogenezinde rol oynar29. m'nin konumları6Aday biyobelirteçlerdeki modifikasyon bölgeleri ve bunların ikincil yapılardaki yüksek güvenilirlikli konumları şurada gösterilmiştir Şekil 7A-HYapılan daha detaylı analizler, temel m6A-modifiye proteinlerin etkileşime girdiği Katalaz AQR ve RBM2'yi içerirken, FABP1 hem SF3A3 hem de AQR ile etkileşime girmiştir. MAOA IGF2BP3 ve IGF2BP2 ile etkileşime girdiği saptanmış olup, MAOB TIA1 ile (Tablo 3).
DN tedavisinde LLF için taksifolin, beta-sitosterol ve eriodiktiyolün uygun in silico bağlanma öngörüleri
miRNet'te, CAT'nin 24 miRNA ile etkileşime gireceği öngörülürken, FABP1 beş miRNA ile ilişkilendirilmiştir. Ayrıca, MAOB ve MAOA sırasıyla 29 ve 26 miRNA ile bağlantılı bulunmuştur. Bunlar arasında, hem TarBase hem de Starbase veri tabanlarında 23 lncRNA tanımlanmıştır. Daha sonra, dört aday biyobelirteç, 74 miRNA ve 23 lncRNA'yı içeren bir lncRNA-miRNA-mRNA düzenleyici ağı oluşturulmuştur (Şekil 8A). Aday biyobelirteçleri hedefleyen potansiyel aktif bileşenler arasında luteolin, beta-sitosterol, eriodiktiyol, kaempferol, kuersetin ve taksifolin yer almıştır (Şekil 8B). Ayrıca; aktif bileşenler, aday biyobelirteçler ve GSEA'da tanımlanan ilk beş yolak temel alınarak bir aktif bileşen-biyobelirteç-yolak ağı kurulmuştur (Şekil 8C). Örneğin, taksifolin peroksizom yolağındaki CAT'yi hedeflemiştir. CAT ile taksifolin (-8.8 kcal/mol), FABP1 ile beta-sitosterol (-8.1 kcal/mol) ve MAOB ile eriodiktiyol (-9.8 kcal/mol) arasındaki bağlanma enerjilerinin tamamı -5 kcal/mol'in altında çıkmış olup, bu durum söz konusu aday biyobelirteçler ile ilgili aktif bileşenleri arasında güçlü afiniteler olduğunu göstermektedir27. Taksifolin, beta-sitosterol ve eriodiktiyol, DN tedavisinde LLF içerisinde uygun in silico bağlanma öngörüleri olan potansiyel aktif bileşenler olarak tanımlanmıştır (Şekil 8D-F). Ancak bunlar, gözlemlenen in vivo etkilerin doğrulanmış biyoaktif ara ürünleri değil, veri tabanı tarafından öngörülen bileşenler olarak sunulmuştur.
DN fare modelindeki aday biyobelirteçleri doğrulayın
DN farelerinin tedavisinde LLF'nin farmakodinamik değerlendirmesi
Uygulama süresi boyunca, farelerde kan glukozu ve idrar mikroalbumin seviyeleri izlenmiştir (Şekil 9A-D). Kontrol grubuyla karşılaştırıldığında, DN model grubunda kan glukozu ve idrar mikroalbumin seviyeleri anlamlı düzeyde artmıştır (P < 0.01); DN model grubuyla karşılaştırıldığında, tedavi grubundaki farelerin kan glukozu 4 haftalık uygulamadan sonra anlamlı derecede azalmış (P < 0.01) ve tedavi grubundaki farelerin idrar mikroalbumin seviyeleri 8 haftalık uygulamadan sonra anlamlı derecede azalmıştır (P < 0.05). Sonuçlar, LLF'nin DN tedavisinde faydalı olabileceğini göstermektedir.
DN farelerinin tedavisinde LLF'nin patolojik değerlendirmesi
HE boyamasının ardından, kontrol grubunun böbrek dokusunda belirgin glomerüler yapılar görülmüştür. Buna karşılık, DN model grubunda, normal grupla karşılaştırıldığında glomerüler nükleer piknoz ve hiperkromazinin yanı sıra glomerüller çevresinde inflamatuar hücre infiltrasyonu gözlenmiştir. LLF tedavisi, db/db farelerinin böbreklerindeki patolojik hasarı iyileştirmiştir (Şekil 9E).
DN farelerinde aday biyobelirteçlerin ekspresyonunun RT-PCR analizi
DN fare modelinin başarıyla oluşturulması ve LLF tedavisiyle semptomlarda belirgin iyileşme gözlemlenmesinin ardından, aday biyobelirteçlerdeki değişiklikleri analiz etmek için RT-qPCR kullanılmıştır. Kontrol grubuyla karşılaştırıldığında, DN grubunda CAT ve MAOA ekspresyonu önemli ölçüde azalmıştır (P < 0.05 veya P < 0.01). Buna karşılık, tedavi grubu DN grubuna göre anlamlı derecede daha yüksek CAT ve MAOA ekspresyonu göstermiştir (P < 0.05). Ancak, gruplar arasında MAOB ve FABP1 ekspresyonunda istatistiksel olarak anlamlı bir fark gözlemlenmemiştir (Şekil 9F-I).
Veri kullanılabilirliği
Bu çalışmada analiz edilen gen ekspresyon veri setlerine, Gene Expression Omnibus (GEO) üzerinden GSE142025 ve GSE96804 erişim numaralarıyla kamuya açık olarak ulaşılabilir. Biyoenformatik analizler için kullanılan R scriptleri, deneysel kaynak verilerle (kan glukozu, üriner mikroalbümin ve RT-qPCR verileri) birlikte Ek Dosya 1'de sunulmuştur. Bu çalışmada kullanılan diğer tüm veri tabanları, yazılımlar ve web kaynakları Materyal Tablosu'nda listelenmiştir.

Şekil 1: Çalışmanın iş akışı. Aday genleri belirlemek için transkriptomik veri setleri, mitokondri ile ilişkili genler ve Ligustri Lucidi Fructus'un öngörülen hedefleri entegre edilmiştir. Ardından, özellik genlerini önceliklendirmek için dört makine öğrenimi algoritması kullanılmış, bunu veri seti arası doğrulama, fonksiyonel karakterizasyon ve db/db farelerde deneysel doğrulama izlemiştir. Kısaltmalar: DN = diyabetik nefropati; DEGs = diferansiyel olarak ifade edilen genler; MRGs = mitokondri ile ilişkili genler; LLF = Ligustri Lucidi Fructus; RF = random forest; KNN = k-en yakın komşu; PLS = kısmi en küçük kareler; SVM = destek vektör makinesi; RMSE = kare ortalama karekök hata; GSEA = gen seti zenginleştirme analizi; RT-qPCR = ters transkripsiyon kantitatif polimeraz zincir reaksiyonu. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 2: DN'nin LLF tedavisi için aday genlerin taranması ve fonksiyonel karakterizasyonu. (A) GSE142025'teki DN ve kontrol örnekleri arasında farklı şekilde eksprese edilen genleri gösteren yanardağ grafiği. (B) |log2FC| değerine göre sıralanmış, en yüksek eksprese olan (upregulated) ilk 10 ve en düşük eksprese olan (downregulated) ilk 10 genin ısı haritası. (C) DEG'lerin, MRG'lerin ve öngörülen LLF hedef genlerinin kesişimini gösteren Venn diyagramı. (D) Aktif bileşen–aday gen ağı. (E) Aday genlerin Gen Ontolojisi zenginleştirme analizi. Çubuk yüksekliği zenginleştirme anlamlılığını temsil eder ve z-skoru fonksiyonel regülasyonun öngörülen yönünü belirtir. (F) Aday genlerin Kyoto Gen ve Genom Ansiklopedisi yolak zenginleştirme analizi. Kısaltmalar: DN = diyabetik nefropati; LLF = Ligustri Lucidi Fructus; DEGs = farklı şekilde eksprese edilen genler; MRGs = mitokondri ile ilişkili genler; GO = Gen Ontolojisi; KEGG = Kyoto Gen ve Genom Ansiklopedisi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 3: Makine öğrenimine dayalı aday biyobelirteçlerin tanımlanması. (A) Aday genler tarafından kodlanan proteinlerin protein-protein etkileşim ağı. (B) RF, KNN, PLS ve SVM modelleri için artıkların ters kümülatif dağılımı. (C) Dört modelin artık dağılımlarını gösteren kutu grafikleri; kırmızı nokta ortalama karekök hatasını belirtir. (D) Dört makine öğrenimi modeli genelinde aday genlerin RMSE tabanlı önemi. (E) Dört modelin tamamında RMSE < 0.281 kriterini karşılayan özellik genlerinin kesişimi. (F,G) Sırasıyla GSE142025 ve GSE96804'teki seçili özellik genlerinin ekspresyonu. Kısaltmalar: RF = rastgele orman; KNN = k-en yakın komşu; PLS = kısmi en küçük kareler; SVM = destek vektör makinesi; RMSE = ortalama karekök hata. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 4: Dört aday biyobelirtecin alıcı işletim karakteristik eğrileri. (A-D) GSE142025 eğitim veri seti ve (E-H) GSE96804 doğrulama veri setindeki CAT, FABP1, MAOB ve MAOA için ROC eğrileri. AUC, alıcı işletim karakteristik eğrisinin altında kalan alanı temsil eder. Kısaltmalar: ROC = alıcı işletim karakteristiği; AUC = eğri altında kalan alan. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 5: Aday biyobelirteçlerin gen seti zenginleştirme analizi. GSE142025 veri setinde (A) CAT, (B) FABP1, (C) MAOA ve (D) MAOB ile ilişkili, anlamlı derecede zenginleşmiş KEGG yolaklarını gösteren GSEA. Kısaltmalar: GSEA = gen seti zenginleştirme analizi; KEGG = Kyoto Gen ve Genom Ansiklopedisi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 6: DN'de immün hücre infiltrasyonu ve bunun aday biyobelirteçlerle ilişkisi. (A) DN ve kontrol örneklerinde CIBERSORT ile tahmin edilen 22 immün hücre tipinin göreceli oranları. (B) DN ve kontrol grupları arasında anlamlı derecede farklı olan immün hücre fraksiyonlarının karşılaştırılması. (C) Diferansiyel bolluğa sahip immün hücre tipleri arasındaki korelasyon matrisi. (D) CAT, FABP1, MAOA ve MAOB ekspresyonları ile diferansiyel bolluğa sahip immün hücre tipleri arasındaki Spearman korelasyonları. Kısaltmalar: DN = diyabetik nefropati. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 7: Aday biyobelirteç transkriptlerinin öngörülen m6A modifikasyon bölgeleri ve RNA ikincil yapıları. (A) CAT, (B) FABP1, (C) MAOA ve (D) MAOB'deki öngörülen m6A modifikasyon bölgeleri. (E) CAT, (F) FABP1, (G) MAOA ve (H) MAOB'nin yüksek güvenilirlikli m6A ile ilişkili bölgelerini gösteren öngörülen RNA ikincil yapıları. Sarı ile vurgulanan bölgeler, m6A modifikasyon bölgelerini içeren öngörülen sekans bölgelerini göstermektedir. Kısaltma: m6A = N6-metiladenozin. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 8: LLF'nin potansiyel aktif bileşenlerinin düzenleyici ağları ve moleküler kenetlenmeleri. (A) Aday biyobelirteçleri içeren öngörülen lncRNA–miRNA–mRNA düzenleyici ağı. (B) Potansiyel LLF aktif bileşenleri ve aday biyobelirteçlerin ağı. (C) GSEA sonuçlarına dayalı aktif bileşen–biyobelirteç–yolak ağı. (D-F) (D) taksifolin ile CAT, (E) beta-sitosterol ile FABP1 ve (F) eriodiktiyol ile MAOB'nin öngörülen moleküler kenetlenme konformasyonları. Kısaltmalar: LLF = Ligustri Lucidi Fructus; lncRNA = uzun kodlamayan RNA; miRNA = mikroRNA; GSEA = gen seti zenginleştirme analizi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 9: db/db farelerde LLF tedavisinin biyokimyasal göstergeler, renal histopatoloji ve aday biyobelirteç ekspresyonu üzerindeki etkileri. (A,B) Sırasıyla başlangıçtaki ve 8. haftadaki kan glukoz seviyeleri. (C,D) Sırasıyla başlangıçtaki ve 8. haftadaki idrar mikroalbümin seviyeleri. (E) Kontrol, DN ve Tedavi gruplarından temsilci hematoksilin ve eosin boyalı böbrek kesitleri (büyütme, ×40; ölçek çubuğu = 25 µm). (F-I) RT-qPCR ile ölçülen sırasıyla Cat, Maoa, Maob ve Fabp1'in bağıl renal mRNA ekspresyon seviyeleri. Kontrol grubuna kıyasla #P < 0.05, ##P < 0.01 ve ###P < 0.01; DN grubuna kıyasla *P < 0.05, **P < 0.01 ve ***P < 0.01. Kısaltmalar: LLF = Ligustri Lucidi Fructus; DN = diyabetik nefropati; RT-qPCR = ters transkripsiyon kantitatif polimeraz zincir reaksiyonu. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
| primer | diziler |
| Katalaz F | TCACTGACGAGATGGCACAC |
| KAT R | ATCGAACGGCAATAGGGGTC |
| FABP1 F | CAATAGGTCTGCCCGAGGAC |
| FABP1 R | GTCATGGTCTCCAGTTCGCA |
| MAOB F | GCACTGAAACAGCCTCACAC |
| MAOB R | TCGTGCAGGGACATCCAAAG |
| MAOA F | ACTTACCCATTCCGTGGTGC |
| MAOA R | ACCACAGGGCAGATACCTCA |
| M-GAPDH F | CCTTCCGTGTTCCTACCCC |
| M-GAPDH R | GCCCAAGATGCCCTTCAGT |
Tablo 1: Fare böbrek dokularının RT-qPCR analizi için kullanılan primer dizileri. Kısaltmalar: F = ileri primer; R = geri primer; RT-qPCR = ters transkripsiyon kantitatif polimeraz zincir reaksiyonu.
| Molekül Kimliği | Molekül adı | OB (%) | DL | Hedef sayı |
| MOL00358 | beta-sitosterol | 36.91 | 0.75 | 100 |
| MOL0042 | kaempferol | 41.88 | 0.24 | 103 |
| MOL04576 | taksifolin | 57.84 | 0.27 | 92 |
| MOL05146 | Lucidumoside D | 48.87 | 0.71 | 104 |
| MOL05147 | Lucidumoside D_qt | 54.41 | 0.47 | 0 |
| MOL05169 | (20S)-24-en-3,20-diol-3-asetat | 40.23 | 0.82 | 0 |
| MOL05190 | eriodictyol | 71.79 | 0.24 | 101 |
| MOL05195 | syringaresinol diglukozid_qt | 83.12 | 0.8 | 0 |
| MOL05209 | Lucidusculine | 30.11 | 0.75 | 105 |
| MOL0521 | Olitoriside | 65.45 | 0.23 | 100 |
| MOL05212 | Olitorisid_qt | 103.23 | 0.78 | 0 |
| MOL06 | luteolin | 36.16 | 0.25 | 102 |
| MOL0098 | kuersetin | 46.43 | 0.28 | 103 |
Tablo 2: TCMSP veri tabanı kullanılarak tanımlanan Ligustri Lucidi Fructus'nun on üç aktif bileşeni. Kısaltmalar: OB = oral biyoyararlanım; DL = ilaç benzerliği.
| mRNA | Protein | RF | SVM |
| CAT | AQR | 0.7 | 0.98 |
| CAT | RBM2 | 0.8 | 0.97 |
| FABP1 | AQR | 0.65 | 0.94 |
| FABP1 | SF3A3 | 0.7 | 0.8 |
| MAOA | IGF2BP2 | 0.75 | 0.97 |
| MAOA | IGF2BP3 | 0.75 | 0.97 |
| MAOB | TIA1 | 0.85 | 0.89 |
Tablo 3: Dört mitokondriyal biyobelirteç mRNA'sı ile m6A ile ilişkili RNA-bağlayıcı proteinler arasındaki öngörülen etkileşimler. CAT, FABP1, MAOA ve MAOB insan biyobelirteç mRNA'larını; AQR, RBM2, SF3A3, IGF2BP2, IGF2BP3 ve TIA1 ise RNA-bağlayıcı proteinleri ifade eder. RF ve SVM skorları > 0,5 tahmin edilen RNA–protein etkileşimlerini belirtir. Kısaltmalar: RF = rastgele orman; SVM = destek vektör makinesi.
Ek Dosya 1. Biyoenformatik betikleri ve deneysel kaynak veriler. Bu arşiv; veri işleme, diferansiyel ekspresyon analizi, fonksiyonel zenginleştirme, makine öğrenimi, alıcı işletim karakteristiği (ROC) analizi, gen seti zenginleştirme analizi, Spearman korelasyon analizi ve CIBERSORT immün hücre infiltrasyon analizi için kullanılan R betiklerinin yanı sıra kan glukozu, idrar mikroalbumin ve RT-qPCR deneylerine ait kaynak verileri içermektedir. Bu dosyayı indirmek için lütfen buraya tıklayınız.