İnsan katılımcıları içeren tüm prosedürler Chongqing Chenjiaqiao Hastanesi Etik Komitesi tarafından incelenip onaylandı (onay No. 20240131). Örnek toplamadan önce tüm katılımcılardan yazılı bilgilendirilmiş onay alındı. Tüm dışkı örnekleri ve katılımcı bilgileri, Materyaller Tablosu'nda listelenen kimliği yok edilmiş çalışma kodları kullanılarak ele alındı.
Katılımcı taraması ve kayıt
Haziran 2023 ile Aralık 2023 arasında Chongqing Chenjiaqiao Hastanesi'ni ziyaret eden yeni teşhisli T2DM hastaları uygunluk açısından tarandı. Katılımcılar, 20 ile 65 yaş arasındaysa, yeni T2DM tanısı almış, taze dışkı örneği vermeyi kabul eden ve gerekli tıbbi geçmiş değerlendirmesi ile klinik muayeneleri tamamlamışlarsa DM grubuna dahil ediliyordu. Aynı dönemde işe alınan sağlıklı gönüllüler, diyabetik olmayan eşleşen (NM) grup olarak dahil edildi ve yaş, cinsiyet ve genel demografik özellikler açısından DM grubuna mümkün olduğunca yakın eşleştirildi. Katılımcılar, hematolojik bozukluklar, merkezi sinir sistemi hastalıkları, aktif romatizmal hastalık, otoimmün hastalık, akut veya kronik gastrointestinal enfeksiyon, kronik ishal, kabızlık, aktif veya iyileşmekçi gastrointestinal ülser, inflamatuar bağırsak hastalığı, irritabl bağırsak sendromu, bağırsak tüberkülozu, gastrointestinal tümörler, diğer kanserler, şiddetli kalp yetmezliği, şiddetli karaciğer yetmezliği, şiddetli böbrek yetmezliği ve diğer şiddetli metabolik durumlara sahiplerse eleme alındı. hastalıklar, yetersiz beslenme, bağışıklık eksikliği, doğuştan metabolik bozukluklar, psikiyatrik hastalıklar, sedatif-hipnotik kullanım veya uyuşturucu kötüye kullanımı. Katılımcılar, dışkı toplamadan önce 1 ay içinde antibiyotik, asit baskılayıcı ajanlar, gastrointestinal motilite ilaçları, probiyotikler, glukokortikoidler veya immünosupresantlar almışlarsa yarış dışı bırakıldı. Ayrıca, dışkı toplanmadan önce 1 ay içinde ishal yaşamış, gastrointestinal cerrahi veya gastrointestinal endoskopi geçiren veya yaşam ortamında veya beslenme alışkanlıklarında ani değişiklikler bildiren kişiler dışlanmıştır. Örnek toplamadan önce her uygun katılımcıya benzersiz bir çalışma kodu atandı. Uygun kohorttan, örnek kullanılabilirliği ve dizileme kalite gereksinimlerine göre 12 yeni teşhis konan T2DM ve 12 NM katılımcısı, örnek kullanılabilirliği ve dizileme kalite gereksinimlerine göre 16S rDNA dizileme ve mikrobiyom analizleri için seçildi.
Dışkı örnekleri toplama ve depolama
Her katılımcıya, kimliğinin silinmesini sağlamak için yalnızca atanan çalışma kodu ile etiketlenmiş steril dışkı toplama kapı sağlandı. Katılımcılara, idrar kontaminasyonunu en aza indirmek için dışkıdan önce mesanelerini boşaltmaları ve tuvalet suyu, idrar, dezenfektan veya diğer potansiyel kirleticilerle temas etmeden taze dışkıyı doğrudan steril kaba toplamaları talimatı verildi. Toplandıktan hemen sonra, yaklaşık 1–2 g dışkı, steril tek kullanımlık örnekleme kaşığı kullanılarak steril kriyojenik tüpe transfer edildi. Tüp sıkıca kapatıldı ve hemen kuru buz üzerine veya önceden soğutulmuş transfer kabına kondu. Tüm örnekler toplandıktan sonra 2 saat içinde laboratuvara taşındı. Laboratuvara varıldığında, çalışma kodu doğrulandı ve her tüp sızıntı veya görünür kontaminasyon açısından incelendi. Tüm örnekler için toplama ve depolama süreleri kaydedildi. Dışkı örnekleri genomik DNA çıkarımı yapılana kadar −80 °C'de depolandı. Numuneler, depolama ve taşıma sırasında tamamen donmuş kalıyorsa ve boru sızıntısı veya dış kirlenme belirtisi bulunmadığında aşağı akış analizi için kabul edilebilir sayılırdı.
Genomik DNA çıkarımı
Dışkı örnekleri −80 °C sıcaklıkta saklanırken işlemden önce buz üzerine kondu. Her örnek, tekrar eden donma-erime döngüleriyle ilişkili bozulmayı en aza indirmek için homojenizasyon mümkün olana kadar çözüldü. Yaklaşık 200 mg dışkı, steril bir mikrosantrifüj tüpüne aktarıldı ve dışkı DNA çıkarma kitinde sağlanan liz tamponu, Malzemeler Tablosu'nda listelenen üretici talimatlarına göre eklendi. Örnekler tam homojenizasyon için 30 saniye boyunca şiddetle girdablandırıldı ve ardından hücre lizisini kolaylaştırmak için oda sıcaklığında 5–10 dakika kuluçka altına alındı. Lizisin ardından, örnekler 12.000 × g. sıcaklıkta 4 °C'de 10 dakika boyunca santrifüjlendi. Süpernatant, pelleti rahatsız etmeden dikkatlice yeni steril bir mikrosantrifüj tüpüne aktarıldı. Genomik DNA, üreticinin protokolüne göre kit ile sağlanan elusyon tamponu veya nükleaz içermeyen su kullanılarak elüasyonlandı. Çıkarılan DNA, kısa süreli kullanım için −20 °C'de veya uzun süreli koruma için −80 °C'de depolanmıştır. Aşağı akış analizi için uygun olarak değerlendirilen DNA örnekleri net ve görünür partikül madde içermeyen görünüyordu.
DNA kalite değerlendirmesi
DNA kalitesi, agaroz jel elektroforezi ve floresans tabanlı niceleme kullanılarak değerlendirildi. %1 agararoz jeli elektroforez tamponunda hazırlandı ve çıkarılan DNA örnekleri, yükleme tamponu ile karıştırıldı ve jel kuyularına yüklendi. Elektroforez işlemi yapıldı, ta ki DNA bantları yeterince ayrılana kadar, ardından jel ultraviyole veya mavi ışık altında incelendi. DNA bütünlüğü, belirgin bozulma veya aşırı bulaşma olmadan sağlam bir genomik DNA bandı gözlemlendiğinde kabul edilebilir sayıldı. DNA konsantrasyonu daha sonra floresans tabanlı bir DNA nicelendirme sistemi kullanılarak nicel edildi. Numuneler, sonraki PCR amplifikasyonu için gereken konsantrasyona kadar seyreltildi. Yüksek kaliteli DNA örnekleri genellikle agaroz jel elektroforezisinde belirgin bir bant gösterir ve sonraki amplifikasyon işlemleri için yeterince konsantre olurdu.
16S rDNA V3–V4 bölgesinin PCR amplifikasyonu
Bakteriyel 16S rDNA'nın V3–V4 hiperdeğişken bölgesi, barkodlu primerlerle amplifikasyon yapıldı. Astarlar şu şekildeydi: ön primer 5′-ACTCCTACGGGAGGCAG-3′ ve ters primer 5′-GGACTACHVGGGTWTCTAAT-3′. PCR reaksiyonları, 4 μL PCR tamponu, 2 μL nükleotid karışımı (2,5 mmol/L), 0,8 μL ileri ve geri primerler (5 μmol/L), 0,4 μL yüksek isabetli DNA polimeraz, 10 ng şablon DNA ve nukleazsız su içeren 20 μL nihai hacimde buz üzerinde hazırlandı. Reaksiyon karışımları pipetle nazikçe karıştırılır ve kısa bir süreliğine santrifüjlenerek sıvının tüpün dibinde toplanmasına izin verilmiştir. PCR amplifikasyonu aşağıdaki döngü koşullarıyla gerçekleştirildi: başlangıç denatürasyonu 95 °C'de 5 dakika yapıldı, ardından 95 °C'de 30 saniye boyunca 25 döngü denatürasyon, 30 saniye boyunca 55 °C'de tavlama ve 30 saniye boyunca 72 °C'de uzatma ile son uzatma adımı 72 °C'de 10 dakika sürdü. Amplifikasyon yanlılığını en aza indirmek için her örnek üçlü reaksiyonlarda amplifikasyon yapıldı ve aynı örnekten elde edilen PCR ürünleri daha sonra tek bir tüpte toplandı. Başarılı amplifikasyon, agaroz jel elektroforezide beklenen boyutta net bir amplikon bandının varlığıyla doğrulandı.
PCR ürün doğrulaması ve arıtma
PCR ürünleri %2 agaroz jel elektroforezi kullanılarak doğrulandı. Her örnekten alınan birleşik PCR ürünleri jel üzerine yüklenir ve hedef amplikon bantları net şekilde ayrılana kadar elektrofore edilirdi. Beklenen amplikon boyutuna uygun bantlar, temiz steril bir bıçakla çıkarıldı ve Malzemeler Tablosu'nda listelenen üretici talimatlarına göre jel çıkarma setiyle saflaştırıldı. Arıtılmış PCR ürünleri daha sonra elusyon tamponu veya nükleaz içermeyen suda elüasyon yapıldı. Her saflaştırılmış PCR ürününün konsantrasyonu, floresans tabanlı DNA nicelendirme sistemi kullanılarak nicelendirildi ve konsantrasyon değerleri aşağı akış kütüphane hazırlığı için kaydedildi. Başarılı amplifikasyon ve saflaştırma, agaroz jel elektroforezide beklenen amplikon boyutunda tek, net bir bant ile doğrulandı.
Kütüphane hazırlığı ve havuzlama
Dizileme kütüphaneleri, üreticinin Materyal Tablosu'nda listelenen talimatlarına göre bir DNA kütüphane hazırlık kiti kullanılarak hazırlanmıştır. Dizileme adaptörü dizileri, standart kütüphane hazırlama iş akışına uygun olarak saflaştırılmış amplikonlara bağlandı. Adaptör ile ligasyonlu ürünler, seçilen kütüphane hazırlama protokolüne bağlı olarak jel ekstraksiyonu veya boncuk tabanlı arıtma yöntemiyle daha sonra saflaştırıldı. Arıtılmış kütüphane ürünleri, uygun parça dağılımını doğrulamak için %2 agararoz jel elektroforezi kullanılarak değerlendirildi. Kütüphane konsantrasyonları, floresans tabanlı DNA nicelik sistemi kullanılarak nicelendirildi. DNA konsantrasyonu ve parça boyutuna dayanarak, her kütüphane aynı azı dişleri konsantrasyonuna normalize edildi ve örnekler arasında dengeli dizileme derinliği sağlandı. Bireysel kütüphanelerin eşit azı miktarları 1:1 oranında toplanarak nihai dizileme kütüphane havuzu oluşturuldu. Dizilemeden önce, havuzlu kütüphane, örnek yüklemeden önce sıralama platformu protokolüne göre denatüre edilmiştir. Dizileme için uygun kabul edilen kütüphaneler, beklenen boyut aralığında net bir parça dağılımı ve aşağı akış dizileme analizi için yeterli konsantrasyon sergilemiştir.
Çift uçlu dizileme
Son havuzlanmış kütüphane, üreticinin standart işletim prosedürlerine uygun olarak eşleştirilmiş uçlu bir dizileme platformuna yüklendi. Eşleştirilmiş uçlu dizileme, her örnekte mikrobiyal çeşitliliğin kapsamlı karakterizasyonu için yeterli dizileme derinliği sağlamak amacıyla ya 2 × 250 bp ya da 2 × 300 bp koşu konfigürasyonları kullanılarak gerçekleştirildi. Dizileme kalitesi, platform tarafından oluşturulan kalite metrikleri kullanılarak çalışma boyunca izlendi. Dizileme tamamlandıktan sonra, her örnek için ham eşleştirilmiş uçlu FASTQ dosyaları aşağı akış biyoinformatik analizi için dışa aktarıldı. Dizileme veri setleri, örnek başına yeterli okuma derinliği, uygun Q30 kalite puanları ve başarılı barkod ataması gösterdiklerinde kabul edilebilir kabul edildi.
Ham dizi işleme
Ham eşleştirilmiş uçlu FASTQ dosyaları, aşağı akış işlemleri için mikrobiyom biyoinformatik analiz boru hattına aktarıldı. Diziler, her numuneye atanan barkodlara göre demultiplex edildi. Daha sonra, düşük kaliteli puanlar, belirsiz tabanlar veya dizileme artefaktları içeren okumalar kaldırılmak için kalite filtreleme gerçekleştirildi. Okumaların sonunda bulunan düşük kaliteli bazlar, genel dizi güvenilirliğini artırmak için kalite puan dağılımına göre kesildi. Kalite kontrolünden sonra, eşleştirilmiş uçlu okumalar örtüşen bölgeler arasında birleştirilerek tam uzunlukta diziler yeniden oluşturuldu. Sonraki mikrobiyom analizleri için yalnızca yüksek kaliteli birleşik diziler tutuldu.
Amplikon dizisi varyantı (ASV) üretimi ve taksonomik açıklama
Diziyle denoizasyon, dizileme hatalarını düzeltmek ve dizis doğruluğunu artırmak için DADA2 veya Deblur dahil olmak üzere doğrulanmış bir denozizasyon algoritması kullanılarak gerçekleştirildi. Kimerik diziler denoizasyon sürecinde tespit edilip çıkarıldı. Her örnek için temsilci ASV dizileri ve ilgili ASV bolluk tabloları oluşturuldu. ASV dizilerinin taksonomik ataması Naive Bayes sınıflandırıcısı kullanılarak yapıldı ve taksonomik notasyon SILVA 138 referans veritabanına göre yapıldı. Taksonomik bolluk tabloları, sonraki analiz için filüm, aile ve cins seviyelerinde ihraç edildi. Son ASV veri seti, her örnek için temsil dizileri, bolluk bilgileri ve ilgili taksonomik açıklamaları içeriyordu.
Alfa ve beta çeşitlilik analizi
Chao, ACE, Shannon, Simpson ve Coverage endeksleri dahil olmak üzere alfa çeşitlilik endeksleri, örnek içindeki mikrobiyal çeşitlilik ve zenginliği değerlendirmek için hesaplandı. Alfa çeşitlilik indeksleri ortalama ± standart sapma olarak ifade edildi. Veri normalliği Shapiro–Wilk testi kullanılarak değerlendirildi. Normal dağılımlı değişkenler DM ve NM grupları arasında Student's t-testleri kullanılarak karşılaştırılırken, normal dağılım olmayan değişkenler Wilcoxon sıralama toplamı testleriyle analiz edildi. Beta çeşitlilik mesafeleri, örnekler arasında mikrobiyal topluluk bileşimindeki farklılıkları değerlendirmek için uygun bir mesafe matrisi kullanılarak hesaplandı. Gruplar arasındaki mikrobiyal topluluk yapısındaki farkları görselleştirmek için ana koordinat analizi (PCoA) gerçekleştirildi. Ayrıca, gruplar arası farkların grup içi varyasyonu aşıp aşmadığını belirlemek için benzerliklerin analizi (ANOSIM) yapıldı ve ilgili ANOSIM R ve p değerleri raporlandı. Alfa çeşitliliği analizi, örnek içi mikrobiyal çeşitliliği yansıtırken, beta çeşitlilik analizi gruplar arasındaki mikrobiyal topluluk bileşimindeki farklılıkları değerlendirdi.
Diferansiyel taksonomik analiz
Mikrobiyal topluluk bileşimi filüm, aile ve cins seviyelerinde özetlenmiştir. Baskın taksonların bireysel örnekler ve çalışma grupları arasında göreceli bolluğunu görselleştirmek için üst üste koymuş çubuk grafikleri oluşturuldu. DM ve NM grupları arasında paylaşılan ve gruba özgü ASV'leri karşılaştırmak için Pan/Core eğrileri ve Venn diyagramları da oluşturuldu. Gruplar arasındaki takson bolluğundaki farklar Wilcoxon sıralama toplamı testleriyle değerlendirildi. Daha sonra LEfSe analizi yapılarak gruplar arası ayrımcılığa en çok katkıda bulunan mikrobiyal taksonlar belirlenmiştir. Diferansiyel taksonomik analiz, DM veya NM grubunda zenginleşmiş mikrobiyal taksonların tanımlanmasını mümkün kıldı.
Fonksiyonel tahmin
PICRUSt2 veya eşdeğer doğrulanmış fonksiyonel tahmin hattı, 16S rDNA dizileme profillerinden mikrobiyal fonksiyonel yollar çıkarmak için kullanıldı. ASV bolluğu verileri, seçilen analitik boru hattının gereksinimlerine göre normalleştirildi ve tahmin edilen fonksiyonlar KEGG veya MetaCyc yol notasyonlarına eşlendi. Tahmin edilen yol bollukları DM ve NM grupları arasında karşılaştırıldı. Önemli ölçüde farklı yollar, ısı haritaları veya diğer uygun grafiksel yaklaşımlarla görselleştirildi. Tahmin edilen fonksiyonlar, doğrudan ölçülen metabolit bolluğu yerine çıkarılan mikrobiyal metabolik potansiyel olarak yorumlandı. Fonksiyonel tahmin analizi, gruplar arasında farklılık gösteren aday mikrobiyal yolların tanımlanmasını mümkün kılmıştır; bunlara karbonhidrat ve aminokislota metabolizmasıyla ilişkili yollar da dahildir.
Mendel rastgeleleştirme analizi
Genom çapında ilişki çalışması (GWAS) özeti bağırsak mikrobiyotası için kamuya açık veri setlerindenalınmıştır 10. Mikrobiyal taksonlar için genom genelinde anlamlı ilişki verileri, NHGRI-EBI GWAS Kataloğu'ndan (https://www.ebi.ac.uk/gwas/) GCST90032172'den GCST90032644'e kadar değişen erişim sayılarıyla alındı. Ek metagenomik veriler, FINRISK 2002 kohortundan Avrupa Genom-Fenom Arşivi (Araştırma Kimliği: EGAS00001005020) aracılığıyla erişildi. Tip 2 diyabet hastalığı için GWAS özet istatistikleri NHGRI-EBI GWAS Kataloğundan da alınmıştır (erişim numarası: EBI-A-GCST006867). Mikrobiyal taksonlarla ilişkili genetik varyantlar, önceden tanımlanmış istatistiksel eşiklere göre arstüman değişkenleri olarak seçildi ve bağlantı dengesizliğindeki varyantlar hariç tutuldu. Maruz kalma ve sonuç veri setleri, tutarlı allel yönelimi sağlamak için uyumlandırıldı. Mendel randomizasyon analizi daha sonra birincil analitik yaklaşım olarak ters varyans ağırlıklı (IVW) yöntemiyle gerçekleştirilmiştir. Duyarlılık analizleri MR-Egger ve ağırlıklı medyan yöntemleri kullanılarak gerçekleştirildi. Alet gücü F-istatistikleriyle değerlendirilirken, enstrümantal değişkenler arasındaki heterojenlik Cochran'ın Q testi kullanılarak değerlendirildi. Yatay pleiotropi MR-Egger kesici testi kullanılarak incelendi. Bonferroni düzeltmesi birden fazla karşılaştırma için uygulanmıştır. Mendel randomizasyon bulguları, nedenselliğe dair kesin kanıt değil, genetik olarak tahmin edilen ilişkiler olarak yorumlandı.
Veri çıkışı ve uç nokta
Son analitik çıktılar arasında ASV bolluk tabloları, taksonomik bileşim grafikleri, alfa çeşitlilik metrikleri, beta çeşitlilik analizleri, farklı taksonomik sonuçlar, öngörülen fonksiyonel yol profilleri ve Mendel rastgeleleştirme tahminleri yer aldı. Dizileme veri setlerindeki tüm örnek tanımlayıcılar, ilgili kimliği yok edilmiş çalışma kodlarıyla tutarlılık sağlamak için doğrulandı. Ham dizileme dosyaları, işlenen ASV tabloları, istatistiksel çıktılar ve şekil kaynak dosyaları aşağı akış analizi ve veri yönetimi için arşivlendi. Yüksek kaliteli dizileme verileri, taksonomik profiller, çeşitlilik metrikleri, öngörülen fonksiyonel yol analizleri ve Mendel randomizasyon sonuçları başarıyla üretilip doğrulandıktan sonra protokol tamamlanmış sayıldı.