Alle procedures waarbij menselijke deelnemers betrokken waren, werden beoordeeld en goedgekeurd door de Ethische Commissie van het Chongqing Chenjiaqiao Ziekenhuis (goedkeuringsnr. 20240131). Schriftelijke geïnformeerde toestemming werd van alle deelnemers verkregen vóór de monsterafname. Alle ontlastingsmonsters en deelnemersinformatie werden behandeld met behulp van gedeïdentificeerde studiecodes, zoals vermeld in de Materiaaltabel.
Selectie en inschrijving van deelnemers
Nieuw gediagnosticeerde patiënten met T2DM die tussen juni 2023 en december 2023 het Chongqing Chenjiaqiao Ziekenhuis bezochten, werden gescreend op geschiktheid. Deelnemers werden ingeschreven in de DM-groep als ze tussen de 20 en 65 jaar oud waren, een nieuwe diagnose T2DM hadden, akkoord gingen met het afnemen van een verse ontlastingsmonster en de vereiste medische anamnesbeoordeling en klinische onderzoeken afrondden. Gezonde vrijwilligers die in dezelfde periode werden gerekruteerd, werden opgenomen als de niet-diabetisch gematchte (NM) groep en zo nauw mogelijk gekoppeld aan de DM-groep wat betreft leeftijd, geslacht en algemene demografische kenmerken. Deelnemers werden uitgesloten als zij hematologische aandoeningen, aandoeningen van het centrale zenuwstelsel, actieve reumatische ziekte, auto-immuunziekte, acute of chronische gastro-intestinale infectie, chronische diarree, constipatie, actieve of genezende gastro-intestinale ulcer, inflammatoire darmziekte, prikkelbare darmsyndroom, darmtuberculose, gastro-intestinale tumoren, andere maligniteiten, ernstige hartinsufficiëntie, ernstige leverinsufficiëntie, ernstige nierinsufficiëntie, andere ernstige metabole aandoeningen ziekten, ondervoeding, immunodeficiëntie, aangeboren stofwisselingsstoornissen, psychiatrische aandoeningen, gebruik van kalmeringsmiddelen-hypnose of drugsgebruik. Deelnemers werden ook uitgesloten als zij binnen 1 maand voor de ontlasting antibiotica, zuuronderdrukkende middelen, gastro-intestinale motiliteitsmedicijnen, probiotica, glucocorticoïden of immunosuppressiva hadden ontvangen. Daarnaast werden personen die diarree hadden gehad, een gastro-intestinale operatie of gastro-intestinale endoscopie hadden ondergaan, of binnen een maand voor de ontlasting plotselinge veranderingen in leefomgeving of eetgewoonten meldden, uitgesloten. Een unieke studiecode werd aan elke in aanmerking komende deelnemer toegekend vóór het verzamelen van monsters. Uit de in aanmerking komende cohort werden 12 nieuw gediagnosticeerde T2DM- en 12 NM-deelnemers geselecteerd voor downstream 16S rDNA-sequencing en microbiome-analyses op basis van de beschikbaarheid van monsters en de kwaliteit van de sequencing.
Ontlastingsmonsterverzameling en opslag
Elke deelnemer kreeg een container voor het verzamelen van steriele ontlasting voorzien van alleen het label van de toegewezen studiecode om deïdentificatie te waarborgen. Deelnemers kregen de instructie om hun blaas te legen vóór het ontlasten om urinevervuiling te minimaliseren en verse ontlasting direct in de steriele container te verzamelen zonder contact met toiletwater, urine, desinfectiemiddelen of andere mogelijke verontreinigingen. Direct na het verzamelen werd ongeveer 1–2 g ontlasting met een steriele wegwerplepel in een steriele cryogene buis overgebracht. De buis werd stevig afgesloten en direct op droogijs of in een voorgekoelde transfercontainer geplaatst. Alle monsters werden binnen 2 uur na de verzameling naar het laboratorium vervoerd. Bij aankomst in het laboratorium werd de studiecode geverifieerd en werd elke buis geïnspecteerd op lekkage of zichtbare besmetting. Verzamel- en opslagtijden werden voor alle monsters geregistreerd. Ontlastingsmonsters werden vervolgens opgeslagen bij −80 °C tot genomische DNA-extractie. Monsters werden als acceptabel beschouwd voor downstream analyse alleen als ze volledig bevroren bleven tijdens opslag en transport en geen bewijs van buislekkage of externe besmetting vertoonden.
Extractie van genomisch DNA
Ontlastingsmonsters werden uit een −80 °C opslag gehaald en op ijs gelegd vóór de verwerking. Elk monster werd alleen ontdooid totdat homogenisatie mogelijk was om de afbraak die samenhangt met herhaalde vries-dooicycli te minimaliseren. Ongeveer 200 mg ontlasting werd overgebracht in een steriele microcentrifugebuis, en de lysisbuffer die in de ontlasting DNA-extractiekit werd geleverd, toegevoegd volgens de instructies van de fabrikant in de Materiaaltabel. Monsters werden 30 seconden lang krachtig gevortexed om volledige homogenisatie te bereiken en vervolgens 5–10 minuten bij kamertemperatuur geïncubeerd om cellyse te bevorderen. Na de lyse werden monsters gecentrifugeerd op 12.000 × g gedurende 10 minuten bij 4 °C. Het supernatant werd voorzichtig overgebracht in een nieuwe steriele microcentrifugebuis zonder de pellet te verstoren. Genomisch DNA werd vervolgens geëlueerd met behulp van de eluatiebuffer die bij de kit of nucleasevrij water werd geleverd, volgens het protocol van de fabrikant. Het geëxtraheerde DNA werd opgeslagen bij −20 °C voor kortdurend gebruik of bij −80 °C voor langdurige conservering. DNA-monsters die geschikt werden geacht voor downstream analyse leken helder en vrij van zichtbaar deeltjesmateriaal.
DNA-kwaliteitsbeoordeling
De DNA-kwaliteit werd beoordeeld met behulp van agarosegel-elektroforese en fluorescentie-gebaseerde kwantificatie. Een 1% agarosegel werd bereid in een elektroforesebuffer, en de geëxtraheerde DNA-monsters, gemengd met de laadbuffer, werden in de gelputten geladen. Elektroforese werd uitgevoerd totdat de DNA-banden voldoende gescheiden waren, waarna de gel werd onderzocht onder ultraviolet of blauw licht. DNA-integriteit werd als acceptabel beschouwd wanneer een intacte genomische DNA-band werd waargenomen zonder duidelijke degradatie of overmatige uitsmeering. De DNA-concentratie werd vervolgens gekwantificeerd met behulp van een fluorescentiegebaseerd DNA-kwantificatiesysteem. Monsters werden vervolgens verdund tot de concentratie die nodig is voor downstream PCR-versterking. Hoogwaardige DNA-monsters vertoonden doorgaans een duidelijk band op agarosegelelektroforese en waren voldoende geconcentreerd voor latere amplificatieprocedures.
PCR-amplificatie van het 16S rDNA V3–V4-gebied
Het V3–V4 hypervariabele gebied van bacteriële 16S rDNA werd versterkt met barcodeprimers. De primersequenties waren als volgt: voorste primer 5′-ACTCCTACGGGAGGCAG-3′ en achterste primer 5′-GGACTACHVGGGTWTCTAAT-3′. PCR-reacties werden op ijs voorbereid in een eindvolume van 20 μL met 4 μL PCR-buffer, 2 μL nucleotidmengsel (2,5 mmol/L), 0,8 μL elk van voor- en omgekeerde primers (5 μmol/L), 0,4 μL hoog-fideliteit DNA-polymerase, 10 ng sjabloon-DNA en nucleasevrij water per volume. Reactiemengsels werden voorzichtig gemengd door pipetten en kort gecentrifugeerd om de vloeistof op de bodem van de buis op te vangen. PCR-versterking werd uitgevoerd met de volgende cyclische condities: initiële denaturatie bij 95 °C gedurende 5 minuten, gevolgd door 25 cycli denaturatie bij 95 °C gedurende 30 seconden, annealing bij 55 °C gedurende 30 s, en extensie bij 72 °C gedurende 30 s, met een laatste extensiestap bij 72 °C voor 10 minuten. Om amplificatiebias te minimaliseren, werd elk monster in driedubbele reacties versterkt, en werden de resulterende PCR-producten uit hetzelfde monster vervolgens samengevoegd in één buis. Succesvolle versterking werd bevestigd door de aanwezigheid van een duidelijke ampliconband van de verwachte grootte bij agarosegelelektroforese.
PCR-productverificatie en zuivering
PCR-producten werden geverifieerd met 2% agarosegelelektroforese. De gecombineerde PCR-producten van elk monster werden op de gel geladen en elektrofore geplaatst totdat de doel-ampliconbanden duidelijk gescheiden waren. Banden die overeenkwamen met de verwachte amplicongrootte werden verwijderd met een schoon steriel mes en gezuiverd met een gelextractiekit volgens de instructies van de fabrikant in de Materiaaltabel. Gezuiverde PCR-producten werden vervolgens geëlueerd in eluatiebuffer of nucleasevrij water. De concentratie van elk gezuiverd PCR-product werd gekwantificeerd met behulp van een fluorescentiegebaseerd DNA-kwantificatiesysteem, en de concentratiewaarden werden geregistreerd voor de downstream bibliotheekvoorbereiding. Succesvolle versterking en zuivering werden bevestigd door een enkele, heldere band bij de verwachte amplicongrootte bij agarosegel-elektroforese.
Bibliotheekvoorbereiding en pooling
Sequencing-bibliotheken werden voorbereid met behulp van een DNA-bibliotheekvoorbereidingskit volgens de instructies van de fabrikant in de Materiaallijst. Sequencing adapter-sequenties werden geligeerd aan de gezuiverde amplicons volgens de standaard workflow voor bibliotheekvoorbereiding. Adapter-geligeerde producten werden vervolgens gezuiverd met behulp van gelextractie of een op kralen gebaseerde zuiveringsmethode, afhankelijk van het geselecteerde bibliotheekvoorbereidingsprotocol. Gezuiverde bibliotheekproducten werden geëvalueerd met 2% agarosegelelektroforese om de juiste fragmentverdeling te bevestigen. Bibliotheekconcentraties werden gekwantificeerd met behulp van een fluorescentie-gebaseerd DNA-kwantificatiesysteem. Op basis van DNA-concentratie en fragmentgrootte werd elke bibliotheek genormaliseerd naar dezelfde molaire concentratie om een evenwichtige sequencingdiepte over de monsters te waarborgen. Gelijke molaire hoeveelheden individuele bibliotheken werden vervolgens samengevoegd in een verhouding van 1:1 om de uiteindelijke sequencing library pool te genereren. Voorafgaand aan het sequencen werd de gepoolde bibliotheek gedenatureerd volgens het sequencingplatformprotocol voordat de sample werd geladen. Bibliotheken die geschikt werden geacht voor sequencing vertoonden een duidelijke fragmentverdeling binnen het verwachte groottebereik en voldoende concentratie voor downstream sequencinganalyse.
Paired-end sequencing
De uiteindelijke gepoolde bibliotheek werd geladen op een paired-end sequencingplatform volgens de standaard operationele procedures van de fabrikant. Paired-end sequencing werd uitgevoerd met een runconfiguratie van 2 × 250 bp of 2 × 300 bp om voldoende sequencingdiepte te garanderen voor een uitgebreide karakterisering van microbiële diversiteit in elk monster. De sequencingkwaliteit werd gedurende de hele run gemonitord met behulp van door het platform gegenereerde kwaliteitsmetrics. Na voltooiing van de sequencing werden ruwe FASTQ-bestanden met gekoppelde eind voor elk monster geëxporteerd voor downstream bio-informatica-analyse. Sequencingdatasets werden als acceptabel beschouwd wanneer ze voldoende leesdiepte per monster aantoonden, geschikte Q30-kwaliteitsscores en succesvolle barcode-toewijzing.
Ruwe sequentieverwerking
Ruwe FASTQ-bestanden met gepaarde eind werden geïmporteerd in een microbiome-bioinformatica-analysepijplijn voor downstream-verwerking. Sequenties werden gedemultiplexeerd op basis van de barcodes die aan elk monster waren toegewezen. Kwaliteitsfiltering werd vervolgens uitgevoerd om reads te verwijderen met scores van lage kwaliteit, ambigue bases of sequencing-artefacten. Laagwaardige bases aan het einde van reads werden getrimd op basis van kwaliteitsscoreverdeling om de algehele betrouwbaarheid van de sequentie te verbeteren. Na kwaliteitscontrole werden gepaarde uiteindelijke reads samengevoegd over overlappende regio's om volledige reeksen te reconstrueren. Voor latere microbioomanalyses werden alleen hoogwaardige samengevoegde sequenties behouden.
Ampliconsequentievariant (ASV) generatie en taxonomische annotatie
Sequentie-denoising werd uitgevoerd met een gevalideerd denoising-algoritme, waaronder DADA2 of Deblur, om sequencingfouten te corrigeren en de sequentienauwkeurigheid te verbeteren. Chimerische sequenties werden geïdentificeerd en verwijderd tijdens het denoisatieproces. Vervolgens werden voor elk monster representatieve ASV-sequenties en bijbehorende ASV-abundantietabellen gegenereerd. De taxonomische toewijzing van ASV-sequenties werd uitgevoerd met een Naive Bayes-classifier, en taxonomische annotatie werd uitgevoerd op basis van de SILVA 138 referentiedatabase. Taxonomische abundantietabellen werden geëxporteerd op het niveau van het fylum, familie en geslacht voor downstream analyse. De uiteindelijke ASV-dataset bevatte representatieve sequenties, abundantie-informatie en bijbehorende taxonomische annotaties voor elke steekproef.
Alfa- en bètadiversiteitsanalyse
Alfa-diversiteitsindexen, waaronder Chao, ACE, Shannon, Simpson en Coverage-indices, werden berekend om de microbiële diversiteit en rijkdom binnen de steekproef te evalueren. Alfa-diversiteitsindices werden uitgedrukt als gemiddelde ± standaarddeviatie. De normaliteit van de gegevens werd beoordeeld met behulp van de Shapiro–Wilk-test. Normaal verdeelde variabelen werden vergeleken tussen de DM- en NM-groepen met behulp van Student's t-tests, terwijl niet-normaal verdeelde variabelen werden geanalyseerd met Wilcoxon rangsomtests. Bètadiversiteitsafstanden werden vervolgens berekend met behulp van een geschikte afstandsmatrix om verschillen in de samenstelling van microbiële gemeenschappen tussen monsters te evalueren. Principal coördinatenanalyse (PCoA) werd uitgevoerd om verschillen in de structuur van microbiële gemeenschappen tussen groepen te visualiseren. Daarnaast werd een analyse van overeenkomsten (ANOSIM) uitgevoerd om te bepalen of de intergroepsverschillen de intragroepsvariatie overschreden, en werden de bijbehorende ANOSIM R- en p-waarden gerapporteerd. Alpha-diversiteitsanalyse weerspiegelde de microbiële diversiteit binnen de steekproef, terwijl de beta-diversiteitsanalyse verschillen in de samenstelling van microbiële gemeenschappen tussen groepen evalueerde.
Differentiële taxonomische analyse
De samenstelling van microbiële gemeenschappen werd samengevat op het niveau van het stam, de familie en het geslacht. Gestapelde bargrafieken werden gegenereerd om de relatieve overvloed van dominante taxa over individuele steekproeven en studiegroepen te visualiseren. Pan/Core-curves en Venn-diagrammen werden verder geconstrueerd om gedeelde en groepsspecifieke ASV's tussen de DM- en NM-groepen te vergelijken. Verschillen in taxonabundantie tussen groepen werden geëvalueerd met behulp van Wilcoxon rangsomtests. Vervolgens werd LEfSe-analyse uitgevoerd om microbiële taxa te identificeren die het sterkst bijdroegen aan intergroepsdiscriminatie. Differentiële taxonomische analyse maakte het mogelijk microbiële taxa te identificeren die verrijkt zijn in de DM- of NM-groep.
Functionele voorspelling
PICRUSt2 of een equivalente gevalideerde functionele voorspellingspijplijn werd gebruikt om microbiële functionele routes af te leiden uit 16S rDNA-sequencingprofielen. ASV-abundantiegegevens werden genormaliseerd volgens de eisen van de geselecteerde analytische pijplijn, en voorspelde functies werden gekoppeld aan KEGG- of MetaCyc-route-annotaties. Voorspelde route-abundanties werden vervolgens vergeleken tussen de DM- en NM-groepen. Significant verschillende paden werden gevisualiseerd met behulp van heatmaps of andere geschikte grafische benaderingen. Voorspelde functies werden geïnterpreteerd als het afgeleide microbiële metabole potentieel in plaats van als direct gemeten metabolietabundantie. Functionele voorspellingsanalyse maakte het mogelijk om kandidaat-microbiële routes te identificeren die tussen groepen verschilden, waaronder routes die geassocieerd zijn met koolhydraat- en aminozuurmetabolisme.
Mendeliaanse randomisatieanalyse
Genome-wide association study (GWAS) samenvattende statistieken voor darmmicrobiota zijn verkregen uit publiek beschikbare datasets10. Genoombrede significante associatiegegevens voor microbiële taxa werden gehaald uit de NHGRI-EBI GWAS-catalogus (https://www.ebi.ac.uk/gwas/) met toegangsnummers variërend van GCST90032172 tot GCST90032644. Extra metagenomische gegevens werden verkregen van de FINRISK 2002-cohort via het European Genome-Phenome Archive (Research ID: EGAS00001005020). GWAS samenvattende statistieken voor type 2 diabetes mellitus zijn ook verkregen uit de NHGRI-EBI GWAS Catalogus (toegangsnummer: EBI-A-GCST006867). Genetische varianten die geassocieerd zijn met microbiële taxa werden geselecteerd als instrumentele variabelen volgens vooraf gedefinieerde statistische drempels, en varianten in koppelingsdisequilibrium werden uitgesloten. Blootstellings- en uitkomstdatasets werden geharmoniseerd om een consistente alleloriëntatie te waarborgen. Mendeliaanse randomisatieanalyse werd vervolgens uitgevoerd met behulp van de inverse variantiegewogen (IVW) methode als primaire analytische benadering. Gevoeligheidsanalyses werden uitgevoerd met behulp van MR-Egger en gewogen mediaanmethoden. De instrumentsterkte werd geëvalueerd met behulp van F-statistieken, terwijl heterogeniteit tussen instrumentele variabelen werd beoordeeld met behulp van de Q-test van Cochran. Horizontale pleiotropie werd onderzocht met behulp van de MR-Egger interceptietest. Bonferroni-correctie werd toegepast om rekening te houden met meerdere vergelijkingen. Mendeliaanse randomisatiebevindingen werden geïnterpreteerd als genetisch voorspelde associaties in plaats van als definitief bewijs van causaliteit.
Data-uitvoer en eindpunt
De uiteindelijke analytische output omvatten ASV-abundantietabellen, taxonomische samenstellingsgrafieken, alfa-diversiteitsmetriek, bètadiversiteitsanalyses, differentiële taxonomische resultaten, voorspelde functionele padprofielen en Mendeliaanse randomisatieschattingen. Alle steekproefidentificaties in de sequencingdatasets werden geverifieerd om consistentie te waarborgen met de bijbehorende gedeïdentificeerde studiecodes. Ruwe sequencingbestanden, verwerkte ASV-tabellen, statistische outputs en bronbestanden van de figuur werden gearchiveerd voor downstream analyse en databeheer. Het protocol werd als compleet beschouwd zodra hoogwaardige sequencinggegevens, taxonomische profielen, diversiteitsmetrics, voorspelde functionele padanalyses en Mendeliaanse randomisatieresultaten succesvol waren gegenereerd en geverifieerd.