Gegevensverzameling
De genexpressiematrix en de bijbehorende klinische gegevens voor de datasets GSE142025 en GSE96804, welke gerelateerd zijn aan DN, werden opgehaald uit de Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/)15. De trainingsset (GSE142025) bevatte niertestmonsters van 27 patiënten met DN en negen controles, die werden gesequenced met het GPL20301 platform. De validatieset (GSE96804) bestond uit sequencinggegevens van 41 DN-patiënten en 20 controles, verwerkt met het GPL17586 platform. Beide datasets richten zich op nierweefsel; de GSE96804 dataset onderzoekt specifiek de glomerulus, de primaire filtratie-eenheid van de nier (Figuur 1). De GSE142025 dataset (trainingsset) bestaat uit monsters van volledig nierweefsel en biedt een uitgebreid overzicht van het transcriptomische landschap van DN. De GSE96804 dataset (validatieset) richt zich daarentegen specifiek op glomerulair weefsel, wat de primaire locatie is van glomerulaire filtratieschade. Aangezien deze twee datasets niet direct werden samengevoegd vanwege verschillen in platform en weefsel, is er geen batch-effectcorrectie toegepast. In plaats daarvan werd cross-dataset validatie onafhankelijk uitgevoerd. Genen met consistente directionele veranderingen en een area under the curve (AUC) groter dan 0,7 in beide datasets werden geselecteerd als robuuste kandidaten, wat de generaliseerbaarheid over verschillende niervakken ondersteunt.
In totaal werden 1.136 mitochondrië-gerelateerde genen (MRG's) geëxtraheerd uit de MitoCarta3.0-database (https://www.broadinstitute.org/mitocarta). De actieve ingrediënten van LLF werden voorspeld met behulp van de Traditional Chinese Medicine Systems Pharmacology (TCMSP)-database (http://sm.nwsuaf.edu.cn/lsp/tcmsp.php), op basis van een drempelwaarde voor orale biobeschikbaarheid (OB) van ≥30% en een drempelwaarde voor drug-likeness (DL) van ≥0,18. Vervolgens werden potentiële doelwitgenen voor de actieve bestanddelen voorspeld met behulp van de Swiss Target Prediction-database (http://www.swisstargetprediction.ch/).
Analyse van differentiële expressie
De analyse van differentiële expressie van GSE142025 (DN vs. controle) werd uitgevoerd met het limma-pakket (v3.54.1), waarbij de significantiecriteria werden gesteld op P.adj < 0.05 en |log2FoldChange| > 0.516. Volcano-plots en heatmaps werden respectievelijk gevisualiseerd met de pakketten ggplot2 (v 3.3.6) en ComplexHeatmap (v 2.14.0)17,18. DEGs, MRGs en potentiële targetgenen van de actieve ingrediënten werden geintersecteerd, en de overlappende genen werden gedefinieerd als kandidaatgenen. Het netwerk dat de actieve ingrediënten en kandidaatgenen verbindt, werd geconstrueerd met behulp van Cytoscape-software (v 3.9.0)19.
Functionele verrijkingsanalyse en constructie van een eiwit-eiwitinteractie (PPI) netwerk
Gene Ontology (GO) en Kyoto Encyclopedia of Genes and Genomes (KEGG) verrijkingsanalyses voor kandidaatgenen werden uitgevoerd met het clusterProfiler-pakket (v 4.6.2) om hun biologische functies en geassocieerde signaleringspaden te onderzoeken (P.adjust < 0,05). De kandidaatgenen werden vervolgens ingevoerd in de STRING-database (https://cn.string-db.org/) om PPI-relaties op te halen (betrouwbaarheidsniveau ≥ 0,4), gevolgd door de constructie van een PPI-netwerk met behulp van Cytoscape (v 3.9.0)20.
Machine learning
Vier machine-learning-algoritmen, waaronder random forest (RF), k-nearest neighbor (KNN), partial least squares (PLS) en support vector machine met een radial basis kernel (SVM), werden geïmplementeerd met het caret-pakket (v6.0-93) op basis van de GSE142025-dataset. De kandidaatgenen die in de voorafgaande analyse waren geïdentificeerd, werden gebruikt als voorspellende variabelen, en de ziektestatus (DN of controle) werd gebruikt als uitkomst. Voor het KNN-model werd 10-voudige kruisvalidatie uitgevoerd met de trainControl-functie, waarbij tuneLength = 10. Het RF-model werd gefit met 20 bomen (ntree = 20); de PLS- en SVM-modellen werden gefit met de instellingen die waren geïmplementeerd in de oorspronkelijke caret-workflow.
Modelresidudistributies werden geëvalueerd met het DALEX-pakket (v2.4.3). Vervolgens werd een permutatiegebaseerde variabelebeoordeling uitgevoerd met de DALEX variable_importance-functie, waarbij de root mean square error (RMSE) als verliesfunctie werd gespecificeerd. De resulterende dropout_loss representeert de RMSE die is verkregen na permutatie van de overeenkomstige variabele21,22. In de oorspronkelijke analytische workflow werden variabelen met een dropout loss < 0,281 over alle modellen behouden, en genen die in alle vier de modellen voorkwamen, werden gedefinieerd als consensus-featuregenen voor daaropvolgende validatie.
De vier machine-learningmodellen werden primair gebruikt voor feature-prioritering in plaats van voor het construeren van een definitieve klinische classifier. Bijgevolg werd de diagnostische discriminatie vervolgens op individueel genniveau geëvalueerd met behulp van ROC-analyse in zowel de discovery- als de validatiedatasets.
Identificatie van kandidaat-biomarkers
In GSE142025 en GSE96804 werden expressieverschillen in kenmerkgenen tussen DN- en controlemonsters geëvalueerd met de Wilcoxon-test. Alleen genen met een significante differentiële expressie (P < 0,05) en concordante expressietrends over de twee datasets werden geselecteerd voor analyse van de receiver operating characteristic (ROC)-curve. Het pROC-pakket (v1.18.0) werd gebruikt om ROC-curves te genereren en de area under the curve (AUC) te berekenen, waarbij genen met een AUC > 0,7 in beide datasets werden geclassificeerd als kandidaat-biomarkers23.
Gene set enrichment analysis (GSEA)
De biologische functies en signaleringspaden die geassocieerd zijn met de kandidaat-biomarkers werden verder onderzocht met behulp van GSEA op de GSE142025-dataset. Eerst werd een Spearman-correlatieanalyse van de kandidaat-biomarkers met alle overige genen uitgevoerd met het psych-pakket (v2.2.9)24. De correlatiecoëfficiënten werden berekend en gerangschikt (van hoog naar laag). De referentie-genenset was c2.cp.kegg.v2023.1.Hs.symbols.gmt uit de Molecular Signatures Database (MSigDB, https://www.gsea-msigdb.org/gsea/msigdb/). Vervolgens werd GSEA uitgevoerd om de verrijking van de gerangschikte genen in de achtergrond-genenset te evalueren met behulp van het clusterProfiler-pakket (v4.6.2). Correctie voor meervoudig testen werd toegepast via de FDR-methode, en aangepaste P-waarden (aangeduid als P.adjust) werden als significant beschouwd indien deze < 0,05 waren.
m6A-modificatieanalyse
Om RNA-methyleringsmodificaties van kandidaat-biomarkers te onderzoeken, is de SRAMP-database (http://www.cuilab.cn/sramp/) gebruikt om m6A-modificatieplaatsen op kandidaat-biomarkers te voorspellen, waarbij de focus lag op posities met een hoge betrouwbaarheid binnen hun secundaire structuren. Vervolgens is de ENCORI-database (https://starbase.sysu.edu.cn/) gebruikt om m6A-gemodificeerde eiwitten te identificeren die interageren met kandidaat-biomarkers, waarbij de parameter |HepG2 (shRNA)| > 1 is gebruikt om sleuteleiwitten te screenen. Daarna is de RPISeq-database (http://pridb.gdcb.iastate.edu/RPISeq/) gebruikt om de waarschijnlijkheid van interacties tussen sleuteleiwitten en kandidaat-biomarkers te voorspellen. De RNA-sequenties van beide zijn in platte tekst geüpload om RF- en SVM-classifier-voorspellingsscores te genereren. Een interactie werd als significant beschouwd wanneer de score boven de 0,5 uitkwam25. De SRAMP-analyse werd uitgevoerd met de 'High'-drempelwaarde voor m6A-plaatsvoorspelling, gebruikmakend van de 'Transcript'-modus en standaardparameters. Voor de ENCORI-analyse is de 'miRNA-mRNA'-interactiefunctie gebruikt met de parameter 'HepG2 (shRNA)' > 1. Voor de RPISeq-analyse is de RF-classifier met standaardparameters gebruikt; scores > 0,5 wezen op een positieve interactie. Dit zijn computationele voorspellingen en geen experimenteel bewijs van m6A-modificatie of RNA-eiwitinteracties in nierweefsel. Het HepG2-shRNA-criterium is afgeleid van de vooraf berekende datasets van ENCORI en weerspiegelt mogelijk geen nierspecifieke regulatie.
Analyse van immueninfiltratie
Het CIBERSORT-algoritme werd toegepast om de proporties van 22 immuunceltypen te schatten in zowel controle- als DN-monsters uit GSE142025, met visualisatie via een heatmap gegenereerd met het ggplot2-pakket (v3.3.6)26. CIBERSORT werd uitgevoerd met gebruik van de LM22-signatuurmatrix, met 1.000 permutaties en uitgeschakelde kwantielnormalisatie (zoals aanbevolen voor microarray-data). Monsters met een CIBERSORT p < 0,05 werden behouden voor verdere analyse. CIBERSORT schat immuuncelfracties op basis van bulkexpressie van nierweefsel, wat compartimentspecifieke infiltratie (bijv. glomerulair versus tubulo-interstitieel) niet kan oplossen en infiltrerende leukocyten niet kan onderscheiden van residente immuuncellen. De gerapporteerde correlaties bevinden zich dus op weefselniveau en dienen te worden gevalideerd met spatiale methoden. De Spearman-correlatieanalyse tussen immuuncellen met een differentieel abundantieniveau en kandidaat-biomarkers werd uitgevoerd met het psych-pakket.
Constructie van netwerken en moleculaire docking
Micro-RNA's (miRNA's) die interageren met kandidaat-biomarkers werden voorspeld met behulp van de miRNet-database (https://www.mirnet.ca). Vervolgens werden long noncoding RNA's (lncRNA's) die gericht zijn op de geïdentificeerde miRNA's voorspeld via de TarBase- (http://www.diana.pcbi.upenn.edu/tarbase) en starbase-databases (http://starbase.sysu.edu.cn/). lncRNA's die in beide databases voorkwamen, werden geselecteerd voor de netwerkconstructie. Vervolgens werd een lncRNA-miRNA-mRNA-regulatienetwerk gebouwd met behulp van Cytoscape-software. Potentiële actieve ingrediënten die gericht zijn op kandidaat-biomarkers werden gekozen om een netwerk van actieve ingrediënten en biomarkers te construeren. Daarnaast werden actieve ingrediënten, kandidaat-biomarkers en pathways die in GSEA waren geïdentificeerd, in Cytoscape geïntegreerd om een netwerk van actieve ingrediënten-biomarkers-pathways te creëren.
Er werd een moleculaire docking-analyse uitgevoerd om de bindingsaffiniteit tussen potentiële actieve bestanddelen en kandidaat-biomarkers te beoordelen. De 3D-structuren van de biomarker-eiwitten werden verkregen in PDB-bestandsformaat uit de Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB, https://www.rscb.org/pdb). De 2D-structuren van de potentiële actieve bestanddelen werden in SDF-formaat opgehaald uit de PubChem-database (http://pubchem.ncbi.nlm.nih.gov). De moleculaire docking werd uitgevoerd met behulp van het CB-Dock-platform (http://clab.labshare.cn/cb-dock/php/blinddock.php). Een bindingsenergie van minder dan -5 kcal/mol duidde op een sterke bindingsaffiniteit27.
Bereiding en authenticatie van Ligustri Lucidi Fructus
Hier verwijst Ligustri Lucidi Fructus (LLF) naar de gedroogde rijpe vrucht van Ligustrum lucidum W. T. Aiton (Oleaceae). Het botanisch materiaal werd geauthenticeerd volgens de Chinese Farmacopee, en een bewijsexemplaar met nummer 20240506,20240911,20241103 werd gedeponeerd bij de Shanxi University of Traditional Chinese Medicine.
Voor de bereiding van het decoct werden 200 g gekwalificeerde LLF-plakjes gedurende 30 min bij kamertemperatuur geweekt in 1.000 mL gedestilleerd water. Het mengsel werd krachtig gekookt en vervolgens 60 min zachtjes laten sudderen. Het filtraat werd verzameld en het resterende kruidenmateriaal werd herhaaldelijk gedecocteerd met opnieuw 1.000 mL gedestilleerd water gedurende 60 min. De twee filtraten werden gecombineerd, gefilterd, gecentrifugeerd en onder vacuüm geconcentreerd om een uiteindelijke voorraadconcentratie van 1 g ruw drugmiddel/mL te verkrijgen (totaal volume van 100 mL). Het bereide decoct werd bewaard bij 4 °C voor kortstondig gebruik of bij −20 °C voor langdurige bewaring. De kwaliteit van LLF en het decoct werd strikt geïdentificeerd en geverifieerd volgens de normen van de Chinese Farmacopee om de betrouwbaarheid en reproduceerbaarheid van de experimenten te waarborgen.
Voor kwalitatieve identificatie werd dunne-laagchromatografie uitgevoerd. Kort samengevat werd een geschikt volume van het bereide decoct gecentrifugeerd en werd het supernatant geëxtraheerd met methanol. Na filtratie werden de monsteroplossing en de referentiestandaardoplossing van specnuezhenide op dezelfde silica gel G-plaat gespot. Na ontwikkeling, droging en inspectie onder ultraviolet licht vertoonde de vlek van de monsteroplossing een overeenstemmende kleur en positie met de referentieverbinding, waarmee het bestaan van de karakteristieke bioactieve component van LLF werd bevestigd.
Voor kwantitatieve kwaliteitscontrole werd detectie via hogedrukvloeistofchromatografie uitgevoerd. De analyse werd uitgevoerd met een C18-kolom en methanol-water als mobiele fase. De detectiegolflengte werd ingesteld op 224 nm. Het gehalte aan specnuezhenide in het decoct werd bepaald op basis van de standaardcurve. De resultaten toonden een stabiele en uniforme chemische samenstelling van het bereide decoct aan, waardoor een consistente drugskwaliteit tijdens het dierexperiment werd gewaarborgd.
Dierproeven
Twaalf mannelijke db/db-muizen van SPF-kwaliteit (8–9 weken oud) en zes leeftijd-gematched db/m-muizen werden gehouden in de SPF-dierfaciliteit van de Shanxi University of Traditional Chinese Medicine. Voorafgaand aan de experimenten werden de dieren 7 dagen geacclimatiseerd onder een cyclus van 12 u licht/12 u donker met ad libitum toegang tot voedsel en water. De studie werd goedgekeurd door de Ethiekcommissie van de Shanxi University of Traditional Chinese Medicine (Goedkeuringsnr. 2022DW167). Dieren met een gewichtsverlies van >20%, een moribunde toestand of onvermogen om toegang te krijgen tot voedsel of water werden humaan geëuthanaseerd vóór het geplande eindpunt van de studie. Aan het einde van de studie werden alle overgebleven muizen geëuthanaseerd via intraperitoneale injectie van pentobarbitalnatrium, gevolgd door cervicale dislocatie in overeenstemming met de institutionele protocollen.
Na de acclimatisatieperiode werd de vestiging van een DN-model in db/db-muizen bevestigd door een bloedglucosegehalte in de staartader van ≥ 16 mmol/L en microalbuminurie, aangetoond door een positieve urinestrip voor microalbumine. Na de succesvolle vestiging van het DN-model werden de db/db-muizen willekeurig verdeeld over twee groepen (n = 6 per groep): de DN-modelgroep (DN) en de LLF-behandelingsgroep (Treatment). Bovendien werden db/m-muizen (n = 6) gebruikt als controlegroep (Control). De dosering werd gekozen op basis van eerdere farmacodynamische studies van LLF bij diabetische ratten en werd omgerekend naar de menselijke equivalente dosis met behulp van normalisatie van het lichaamsoppervlak28. De Control- en DN-groepen kregen gedestilleerd water, terwijl de behandelingsgroep gedurende 8 weken 3,5 g/kg LLF ontving. Na 8 weken van toediening werden alle muizen geëuthanaseerd om serum, urine en nierweefsel te verzamelen voor vervolganalyses.
Bloed- en urine-indicatoren
Serumglucosewaarden werden geanalyseerd met een volledig automatische bloedbiochemie-analyser. De urinaire microalbumineconcentratie werd gemeten volgens de instructies van de kit (Aanvullend bestand 1).
Pathologische observatie van muisnierweefsels
Nierweefsels werden verwerkt voor histopathologisch onderzoek. Na fixatie in 4% paraformaldehyde werden de weefsels gewassen, gedehydrateerd, ingebed in paraffine en in coupes gesneden. Vervolgens werd hematoxyleen en eosine (HE) kleuring uitgevoerd, waarna de gekleurde coupes onder een optische microscoop werden onderzocht om pathologische veranderingen te evalueren.
Kwantitatieve reverse transcriptie polymerasekettingreactie (RT-qPCR)
RT-qPCR werd gebruikt om de expressie van de kandidaat-biomarkers in muisnierweefsel te bepalen. Totaal RNA werd geëxtraheerd volgens de instructies van de fabrikant, waarna de RNA-concentratie en -kwaliteit werden beoordeeld (Tabel 1). cDNA werd gesynthetiseerd uit het geëxtraheerde RNA met behulp van de cDNA Synthesis Kit. Amplificatie werd uitgevoerd met de primerparen vermeld in Tabel 1, waarbij GAPDH als referentiegen werd gebruikt. Relatieve expressieniveaus werden berekend met de 2−ΔΔCt methode13,26.
Statistische analyse
Alle statistische analyses zijn uitgevoerd met R-software (versie 4.2.2) in combinatie met de software die vereist was voor de overeenkomstige experimentele procedures. Tenzij anders aangegeven, waren alle statistische tests tweezijdig en werden verschillen als statistisch significant beschouwd bij P < 0,05. De transcriptomische analyse van differentiële expressie werd uitgevoerd met het limma-pakket. Genen met een aangepaste P < 0,05 en een absolute log2fold change groter dan 0,5 werden gedefinieerd als differentieel tot expressie gebracht.
Voor vergelijkingen van de expressie van kandidaatgenen tussen onafhankelijke DN- en controlesamples werden non-parametrische Wilcoxon rank-sum tests gebruikt, zoals aangegeven in de oorspronkelijke analytische workflow. Correlaties tussen kandidaatbiomarkers en immuuncelfracties werden beoordeeld met behulp van de rankcorrelatie van Spearman.
Experimentele gegevens worden weergegeven als het gemiddelde ± SD. Vergelijkingen tussen de drie onafhankelijke diergroepen werden uitgevoerd met behulp van een eenweg-variantieanalyse wanneer aan de aannames voor parametrische analyse werd voldaan. Voor post-hocvergelijkingen werd de test voor het kleinst significante verschil gebruikt wanneer de varianties homogeen waren, terwijl de T3-test van Dunnett werd gebruikt wanneer de varianties ongelijk waren. De baseline- en week-8-farmacodynamische metingen werden afzonderlijk geanalyseerd en gepresenteerd; er is geen conclusie getrokken met betrekking tot een interactie tussen groep en tijd. De R-scripts en brondata die zijn gebruikt voor de bio-informatica- en machine-learninganalyses zijn beschikbaar in Aanvullend bestand 1.