Wszystkie procedury na zwierzętach były przeprowadzane zgodnie z wytycznymi instytucjonalnymi i zatwierdzone przed rozpoczęciem badań przez Komitet Etyczny Uniwersytetu Tradycyjnej Medycyny Chińskiej w Chengdu (nr Zatwierdzenia 2022-126). Do analiz bioinformatycznych wykorzystywano wyłącznie publicznie dostępne zbiory danych GEO i GWAS; dlatego nie wymagano dodatkowej zgody świadomej ani zatwierdzenia etyki dla komponentu obliczeniowego tego badania. Wszystkie procedury dotyczące niebezpiecznych chemikaliów, próbek biologicznych i odpadów zwierzęcych były wykonywane zgodnie z instytucjonalnymi przepisami bezpieczeństwa laboratoriów oraz procedurami utylizacji odpadów niebezpiecznych.
Pozyskiwanie i wstępne przetwarzanie zbiorów danych GEO
Zbiory danych transkryptomicznych u ludzi związane z wrzodziejącym zapaleniem jelita grubego (UC) zostały pobrane z bazy GEO przy użyciu słowa kluczowego "wrzodziejące zapalenie okrężnicy". Po przesiewie zestawów danych profilujących ekspresję z jasno zdefiniowanymi grupami kontrolnymi UC i zdrowymi grupami kontrolnymi oraz odpowiednią wielkością próby, wybranoGSE87466 20 iGSE92415 21, oba wygenerowane na platformie GPL13158, jako kohortę treningową. Te zbiory danych zawierały odpowiednio 87 próbek UC i 21 zdrowych kontroli, a także 53 próbki UC i 21 zdrowych kontroli. Do walidacji zewnętrznej GSE8747322 (106 próbek UC i 21 zdrowych kontroli) oraz GSE1687923 (48 próbek UC i 12 zdrowych kontroli); platformy GPL570). Pliki surowe lub macierzyste szeregowe były pobierane i importowane do R do wstępnego przetwarzania. Identyfikatory sond były konwertowane na symbole genów zgodnie z odpowiadającymi plikami adnotacji platformy. Gdy wiele sond przypisywano temu samemu symbolowi genu, średnia wartość ekspresji była zachowywana. Próbki bez jasnych informacji grupowych zostały wykluczone przed analizą dalszą. Oba zbiory danych treningowych zostały połączone przez wspólne symbole genów. Wartości wyrażeń były przekształcane log2, gdy było to potrzebne, normalizowane za pomocą standardowej procedury normalizacji międzytablicowej oraz korygowane pod kątem efektów wsadowych metodą korekty partii po zdefiniowaniu pochodzenia zbioru danych jako zmiennej wsadowej. Analiza głównych komponentów oraz inspekcja boxplot były stosowane przed i po normalizacji/korekcie wsadowej jako punkty kontrolne jakości. Skuteczne wstępne przetwarzanie zostało zdefiniowane przez lepsze nakładanie się rozkładów próbek w różnych zbiorach danych oraz tłumienie klastrowania opartego na zbiorach danych. Łącznie 902 geny związane z PANoptosis zostały zebrane z opublikowanej literatury i wykorzystane jako referencyjny zestaw genów do analizy integracyjnej.
Identyfikacja genów ekspresji różnicowo i budowa sieci WGCNA
Scalona macierz treningowa była analizowana w R przy użyciu workflow różniczkowego ekspresji. Geny z|log 2-krotną zmianą| ≥ 0,585 oraz skorygowane P < 0,05 zostały zdefiniowane jako geny różnicowo ekspresyjne (DEGs). Wykresy wulkanów i mapy ciepła zostały wygenerowane jako wyniki pośrednie, aby potwierdzić, że kryteria filtracji dają biologicznie interpretowalne różnice ekspresji między grupami UC a kontrolnymi. W analizie ważonej sieci współekspresji genów (WGCNA) geny były klasyfikowane według wariancji w różnych próbkach, a 25% najczęstszych genów o najbardziej zmiennych zachowano jako dane wejściowe. Najpierw zbadano drzewo klasteryzujące próbki w celu identyfikacji potencjalnych próbek odstających od rzeczywistości; Nie zachowano oczywistych odstających punktów odchylających do budowy sieci, chyba że usunięcie było uzasadnione kryteriami kontroli jakości. Wybrano moc o miękkim progu na podstawie indeksu dopasowania topologii bezskalowej, a najmniejsza moc osiągająca sieć w miarę bezskalowaną użyto do skonstruowania macierzy sąsiedztwa. Macierz sąsiedztwa została następnie przekształcona w topologiczną macierz nakładania się, a hierarchiczne klasteryzowanie przeprowadzono przy użyciu nierówności opartej na TOM. Moduły zostały zidentyfikowane za pomocą dynamicznego cięcia drzew o minimalnym rozmiarze modułu wynoszącym 100 genów. Blisko spokrewnione moduły były łączone, gdy ich korelacja genów własnych przekraczała ustalone kryterium scalania. Obliczano wartości genów własnych, istotności modułów oraz istotności genu, aby zidentyfikować moduły najsilniej związane z fenotypem UC. Udany krok WGCNA został zdefiniowany przez stabilne rozdzielenie modułów, biologicznie prawdopodobne rozmiary modułów oraz jeden lub więcej modułów wykazujących wyraźną korelację ze stanem choroby.
Identyfikacja kandydatów na cele w GQD
Kandydatskie związki i cele odpowiadające czterem składnikom ziołowym wywaru Gegen Qinlian (GQD) – Pueraria lobata, Scutellaria baicalensis, Coptis chinensis i Glycyrrhiza uralensis – zostały zebrane z TCMSP24 i BATMAN-TCM25. Progi przesiewowe zostały ustalone następująco: doustna biodostępność ≥ 30%, a podobność do leku ≥ 0,18 dla rekordów pochodzących z jednej bazy danych; wynik zaufania ≥ 0,84, a skorygowany P < 0,05 dla rekordów pochodzących z innej bazy danych. Dodatkowe związki uzupełniano z bazy danych ziołnictwa, gdy odpowiednie zapisy nie zostały zarejestrowane podczas początkowych poszukiwań. Kanoniczne łańcuchy związków zachowanych SMILES zostały przesłane do internetowej platformy oceny ADME. Związki były zatrzymywane, gdy wykazywały wysoką absorpcję w przewodzie pokarmowym i spełniały co najmniej dwa zestawy zasad podobności do leków u Lipinskiego, Ghose, Vebera, Egana i Muegge. Potencjalne cele były następnie przewidywane za pomocą platformy cel-predykcyjne z prawdopodobieństwem > 0,1 i standaryzowane do oficjalnych symboli genów przy użyciu bazy danych adnotacji białek. Zduplikowane związki i zduplikowane cele zostały usunięte po scaleniu baz danych. Oczekiwanym wynikiem tego etapu był nieredundantny zbiór danych złożonych celów GQD odpowiedni do analizy nakładania się.
Identyfikacja nakładających się genów i wykonanie analizy funkcjonalnej
Potencjalne cele terapeutyczne zostały zdefiniowane jako przecięcie celów związanych z GQD, genów związanych z UC zidentyfikowanych na podstawie różnicowej ekspresji i WGCNA oraz zestawu genów związanych z PANoptozą. Relacje nakładania się zobrazowano za pomocą diagramu Venna. Następnie przecinające się geny zostały przekazane do bazy danych interakcji białko-białko (PPI), a pary interakcji z oceną ufności ≥ 0,4 zostały zachowane. Powstała sieć PPI została wyeksportowana i zwizualizowana w platformie analizy sieci. Udany krok PPI definiował obecność połączonej sieci zawierającej większość przecinających się genów, a nie tylko izolowane węzły. Analiza wzbogacenia funkcjonalnego, w tym analiza ontologii genów (GO) oraz analizy Encyklopedii Genów i Genomów z Kioto (KEGG), została przeprowadzona w języku R, z q < 0,05 jako progiem istotności. Terminy lub ścieżki niespełniające ten próg zostały wykluczone. Podsumowano najwyżej wzbogacone procesy biologiczne i szlaki sygnalizacyjne, aby wywnioskować potencjalny mechanizm, przez który GQD może regulować UC. Wykresy kreskowe i bąbelkowe zostały generowane jako wyniki pośrednie, aby zweryfikować, czy wyniki wzbogacenia są zdominowane przez szlaki związane z uszkodzeniem nabłonka lub stanu zapalnego, odporności lub nabłonkowego, co można się spodziewać w kontekście biologicznym.
Analiza wzorców ekspresji, lokalizacji chromosomowej i korelacji
Poziomy ekspresji wspólnych genów porównano między grupami UC a kontrolnymi przy użyciu standardowego workflow wizualizacji różnicowej ekspresji w R. Opracowano wykresy pudełkowe i mapy cieplne, a P < 0,05 uznano za statystycznie istotne. Zachowano tylko geny o interpretowalnym kierunku ekspresji i spójnym zachowaniu w całej kohorcie treningowej do późniejszej oceny mechanicznej. Lokalizacja chromosomów wybranych genów rdzeniowych została mapowana po sformatowaniu plików adnotacji genów do wymaganej struktury wejściowej, a pozycje genomowe zwizualizowano za pomocą workflow wykresów chromosomów kołowych. Następnie obliczono macierze korelacji między genami rdzeniowymi za pomocą analizy korelacji parowej. Oczekiwanym rezultatem tego kroku była identyfikacja par genów wykazujących skoordynowane wzorce ekspresji, które mogłyby wspierać wspólną regulację lub interakcję funkcjonalną.
Konstrukcja i ocena modeli uczenia maszynowego
Rozwój modeli diagnostycznych realizowano przy użyciu połączonej kohorty treningowej (GSE87466 + GSE92415), natomiast GSE16879 i GSE87473 były wykorzystywane jako niezależne zewnętrzne zbiory danych walidacyjnej. Dwanaście algorytmów uczenia maszynowego – Lasso, Ridge, Stepglm, XGBoost, Random Forest, Enet, plsRglm, GBM, NaiveBayes, LDA, glmBoost i SVM – zostało użytych do wygenerowania 113 kombinacji modeli26. W kohorcie treningowej zastosowano dziesięciokrotną walidację krzyżową. W każdym etapie treningu modelu, wyboru cech oraz estymacji wydajności przeprowadzano przy użyciu tej samej logiki partycjonowania treningowego/testowego, aby uniknąć wycieku informacji. Wydajność modelu oceniano głównie na podstawie powierzchni pod krzywą charakterystyki pracy odbiornika (AUC). Model optymalny został zdefiniowany jako model z najwyższym średnim AUC w zbiorach treningowych i walidacyjnych, a nie sam z najwyższym AUC treningowym. Modele wykazujące bardzo wysokie AUC treningowe, ale słabą wydajność walidacji zewnętrznej nie zostały zachowane. Dla każdego zachowanego modelu wygenerowano krzywe ROC, a z ostatecznych wybranych biomarkerów skonstruowano nomogram. Etap modelowania definiował stabilna wydajność walidacji krzyżowej oraz zachowanie zdolności dyskryminacyjnej w zewnętrznych zbiorach danych.
Wydajność GSVA i ekspresja genów rdzeniowych profilu
Analiza zmienności zestawu genów (GSVA) została przeprowadzona w celu porównania wzbogacenia szlaków między grupami o wysokiej i niskiej ekspresji, zdefiniowanych przez medianę wartości ekspresji każdego genu rdzeniowego. Skorygowany P < 0,05 został użyty jako próg dla istotnych różnic w szlakach. Szlaki niespełniające tego progu zostały wykluczone z interpretacji. Opracowano mapy ciepła i wykresy korelacji, aby zilustrować aktywność szlaków i wzorce współekspresji. Pośredni sukces tego etapu wskazywały profile wzbogacenia zgodne z zaburzeniami zapalnymi, immunologicznymi lub metabolicznymi istotnymi dla UC.
Analiza krajobrazu odpornościowego i konstrukcja sieci ceRNA
Skład komórek immunologicznych oszacowano za pomocą algorytmu dekonwolucyjnego z 1000 permutacji. Próbki spełniające kryterium istotności wewnętrznej dekonwolucji zostały wykluczone z późniejszej analizy korelacji immunologicznej. Aby uzupełnić tę analizę, zastosowano również analizę wzbogacenia zestawów genów na jednej próbie (ssGSEA) do oceny wzorców infiltracji immunologicznej. Następnie obliczono i zobrazowano korelacje między frakcjami komórek odpornościowych a ekspresją genów rdzeniowych. Do budowy sieci ceRNA zidentyfikowano miRNA atakujące rdzenie mRNA z publicznie dostępnych baz danych interakcji, a kandydaty lncRNA oddziałujące z tymi miRNA zostały następnie przesiewane. Zachowano tylko interakcje wspierane przez rekordy predykcji lub adnotacji bazy danych. Sieć regulacyjna mRNA-miRNA-lncRNA została zwizualizowana na platformie analizy sieciowej. Oczekiwanym efektem tego etapu była ustrukturyzowana sieć ceRNA z biologicznie wiarygodnymi regulatorami upstream powiązanymi z genami rdzeniowymi.
Wydajność analizy pojedynczej komórki
Zestaw danych pojedynczych komórek GSE214695 został zaimportowany do R i przekształcony w obiekty pojedynczych komórek przy użyciu standardowego workflow analizy pojedynczej komórki. Kontrola jakości była przeprowadzana poprzez filtrowanie komórek według standardowych wskaźników, w tym komórek o niskich cechach, potencjalnych pustych kropelek oraz komórek o słabej jakości transkryptu. Komórki przechodzące kontrolę jakości zostały znormalizowane, a do analizy głównych komponentów wybrano 2000 genów o wysokiej zmienności. Główne komponenty przyczyniające się do istotnej wariancji biologicznej zostały zachowane dla klasteryzacji bez nadzoru, a UMAP wykorzystano do wizualizacji dwuwymiarowej. Geny markerów klastrowych zidentyfikowano za pomocą analizy różnicowej ekspresji między klastrami. Adnotacje typu komórek były następnie wykonywane przez klasyfikację opartą na odniesieniach, a adnotacje były weryfikowane względem kanonicznego wyrażenia markerów, gdy to możliwe. Skuteczny krok pojedynczej komórki został zdefiniowany przez wyraźne rozdzielenie głównych populacji komórek na wykresie UMAP, interpretowalne geny markerowe dla każdego klastra27 oraz spójne adnotacje głównych typów komórek immunologicznych i stromalnych istotnych dla UC.
Wykonanie analizy randomizacji mendelowskiej
Przeprowadzono randomizację mendelowskią (MR) na dwóch próbach w celu oceny możliwych związków przyczynowych między ekspresją genów rdzeniowych a podatnością na UC. Zmienne instrumentalne zostały wyodrębnione z ilościowych zbiorów danych o cechach ilościowych, a jako zbiór danych wyników użyto podsumowań badań asocjacyjnych UC w całym genomie. SNP były testowane zgodnie ze standardowymi kryteriami kontroli jakości MR, w tym pod kątem znaczenia cechy ekspozycji oraz usuwaniem niejednoznacznych lub duplikatów wariantów, gdy wymagała tego harmonizacja. Metoda ważona odwrotną wariancją była używana jako główny estymator przyczynowy. Test Q Cochrana został zastosowany do oceny heterogeniczności. Analizy MR-Egger i MR-PRESSO zostały wykorzystane jako analizy czułości do oceny i, gdy było to konieczne, korekty poziomej plejotropii. Pomyślny krok MR został zdefiniowany przez prawidłową harmonizację instrumentów, brak dużej heterogeniczności lub plejotropii28 oraz spójne kierunkowo szacunki pomiędzy komplementarnymi metodami MR.
Wydajność dokowania molekularnego związków kandydatów
Trójwymiarowe struktury białek docelowych pozyskano z Protein Data Bank, a struktury ligandów z publicznej bazy danych małych cząsteczek. Przed dokowaniem białka i ligandy były wstępnie przetwarzane poprzez usuwanie cząsteczek wody, gdy było to stosowne, dodawanie wodorów, definiowanie typów atomów oraz konwersję formatów plików wymaganych do dokowania. Następnie dokowanie przeprowadzono za pomocą procesu dokowania molekularnego, a dla każdej pary ligand-cel obliczono energie wiązania. Powinowactwo wiązań poniżej -5 kcal/mol interpretowano jako korzystne wiązanie, natomiast wartości poniżej -7 kcal/mol jako stosunkowo silne wiązanie. Mapy termiczne dokowania zostały narysowane, aby porównać globalne wzorce wiązania między celami. Najwyżej oceniane konformacje zostały zwizualizowane w trzech wymiarach w celu zbadania wiązań wodorowych, interakcji hydrofobowych oraz bliskości pozostałości. Pomyślne zakończenie tego etapu definiowano na podstawie stabilnych plików wyjściowych dokowania, wiarygodnych pozycji ligandów w obszarze wiązania docelowego oraz energii wiązania zgodne z priorytetyzacją kandydatów.
Wydajność walidacji eksperymentalnej in vivo
Przygotowanie GQD
GQD przygotowywano z czterech tradycyjnych chińskich ziół leczniczych w następujących proporcjach surowego leku: 24 g Puerariae Lobatae Radix (Gegen), 9 g Scutellariae Radix (Huangqin), 9 g Coptidis Rhizoma (Huanglian) oraz 6 g Glycyrrhizae Radix et Rhizoma (Gancao). Szczegóły dotyczące materiałów botanicznych i ich źródeł znajdują się w Tabeli Materiałów. Zioła mielono na drobny proszek, dokładnie mieszamy i moczymy w 390 mL wody destylowanej przez 1 godzinę. Mieszanka była wyklejana dwukrotnie. Do pierwszego wywaru zioła gotowano na wysokim ogniu, a następnie gotowano na niskim ogniu przez 30 minut. Po przefiltrowaniu przez gazę resztki były ponownie gotowane z dziesięciokrotną objętością wody. Oba filtraty zostały połączone i skoncentrowane w kąpieli wodnej do ostatecznego stężenia surowego leku 5 g/mL. Następnie ekstrakt był filtrowany, chłodzony i przechowywany w temperaturze 4 °C do czasu użycia.
Ustanowienie modelu myszy UC indukowanego przez DSS
W badaniu wykorzystano osiemnaście samców myszy BALB/c (8 tygodni), a wszystkie procedury zwierzęce zostały zatwierdzone przez Komitet Etyki Uniwersytetu Tradycyjnej Medycyny Chińskiej w Chengdu (nr zatwierdzenia 2022-126). Szczegóły dotyczące źródła zwierzęcego znajdują się w Tabeli Materiałów. Myszy trzymano w określonych warunkach wolnych od patogenów w temperaturze 20–22 °C i wilgotności 55% w cyklu światła/ciemności 12 godzin. Po tygodniowym okresie aklimatyzacji zwierzęta zostały losowo podzielone na trzy grupy (n = 6 na grupę): NC, UC i GQD. Grupa NC otrzymywała przez cały eksperyment sterylną wodę pitną, natomiast grupy UC i GQD otrzymywały 1,5% (w/v) dektran-siarczanu sodu (MW 36-50 kDa) w wodzie pitnej przez 5 kolejnych dni, a następnie przez 2 dni regularną wodę. Ten cykl powtórzono trzy razy, aby wykryć przewlekłe zapalenie jelita grubego. W okresie interwencji grupy NC i UC były podawane dwa razy dziennie 200 μL karboksymetylocelulozy sodu o stężeniu 0,5%, natomiast grupa GQD otrzymywała GQD doustnie dwa razy dziennie.
Kliniczna recepta GQD obejmowała 24 g Gegen, 9 g Huangqin, 9 g Huanglian oraz 6 g Gancao, co odpowiadało całkowitej surowej dawce leku 48 g/dzień dla dorosłych. Zakładając standardową masę ciała dorosłego wynoszącą 60 kg, równoważna dawka dla dorosłych wynosiła 0,8 g/kg/dzień. Na podstawie konwencjonalnego współczynnika konwersji człowieka na mysz wynoszący 1:12, odpowiadająca dawka ekwiwalentna myszy wynosiła 9,6 g/kg/dzień. Zgodnie z ustalonym stosunkiem dawek niskich, średnich i wysokich wynoszących 1:2:4, obliczone dawki wynosiły odpowiednio 9,6, 19,2 i 38,4 g/kg/dzień. Biorąc pod uwagę nasilenie i utrzymujący się stan zapalny jelita grubego w modelu powtarzającego się zapalenia jelita grubego wywołanego DSS, wybrano program wysokich dawek, aby zapewnić odpowiednią interwencję farmakologiczną podczas walidacji in vivo. Aby zapewnić odpowiednią ekspozycję farmakologiczną podczas aktywnej fazy zapalnej, całkowita dawka dzienna wynosząca 38,4 g/kg/dzień została podzielona na dwie równe dawki. Przez cały eksperyment codziennie monitorowano masę ciała, konsystencję stolca oraz krwawienie z odbytu. Na koniec protokołu myszy zostały uśpione, a tkanki jelita grubego pobrane do kolejnych analiz.
Ocena wskaźnika aktywności masy ciała i choroby
Masa ciała była rejestrowana codziennie i wyrażana jako procent początkowej masy ciała w dniu 0. Ciężkość choroby oceniano za pomocą wskaźnika aktywności choroby (DAI), który obejmował trzy parametry: utratę masy ciała, spójność stolca oraz krwawienie z odbytu. Każdy parametr był oceniany niezależnie w skali od 0 do 4 w zależności od nasilenia objawów klinicznych, przy czym wyższe wyniki wskazywały na bardziej ciężką aktywność chorobową. Ostateczny DAI został obliczony jako średnia trzech indywidualnych wyników, co dało łączny wynik w zakresie od 0 do 4. Krótko mówiąc, wynik 0 oznaczał brak nieprawidłowości, podczas gdy wyniki 1–4 odzwierciedlały stopniowo większą utratę masy ciała, luźniejszą konsystencję stolca oraz silniejsze krwawienia z odbytu. Szczegółowe kryteria punktacji przedstawiono w Tabeli 1.
Wydajność barwienia H&E i obliczanie indeksu aktywności histologicznej
Do oceny histologicznej tkanki jelita grubego były utrwalane w 4% paradehydu, odwodnione przez rozwarstwiany etanol, osadzone w parafinie i przecięte przy grubości 4 μm. Fragmenty zostały odparafinizowane, ponownie nawodnione, barwione hematoksyliną, kontrbarwowane eozyną, odwodnione, oczyszczone i zamontowane. Obrazy histologiczne obserwowano za pomocą mikroskopu jasnego pola i skanowano za pomocą systemu szkieletowego skanowania. Szczegółowe informacje o sprzęcie znajdują się w Tabeli Materiałów. Uszkodzenie histopatologiczne oceniano na podstawie trzech parametrów, w tym wyczerpania komórek pucharowych, uszkodzeń architektonicznych krypt oraz infiltracji komórek zapalnych. Każdy parametr był punktowany według stopnia uszkodzenia tkanki, a suma trzech parametrów była definiowana jako indeks aktywności histologicznej (HAI), w zakresie od 0 do 10, jak podsumowano w Tabeli 2. Do ilościowej analizy komórek pucharowych przeanalizowano losowo pięć niepokrywających się pól na sekcję przy powiększeniu 400×. Liczba komórek pucharowych na pole była liczona ręcznie, w sposób ślepy, a średnia liczba na pole była obliczana dla każdego reprezentatywnego odcinka okrężnika.
Wydajność barwienia immunohistochemicznego
Po deparafinizacji i nawodnieniu przekroje tkanek zostały poddane pobieraniu antygenów z użyciem 0,1% trypsyny w temperaturze 37 °C przez 30 minut, a następnie inkubacji z 3% nadtlenkiem wodoru w celu zablokowania endogennej aktywności nadoksydazy. Następnie sekcje były inkubowane przez noc w temperaturze 4 °C z pierwotnymi przeciwciałami królików przeciwko STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP) oraz HSPA5 (1:1000, Cell Signaling Technology, #3177), rozcieńczonych w 1% BSA. Szczegóły dotyczące przeciwciał i ich źródeł podano w Tabeli Materiałów. Po umyciu PBS sekcje były inkubowane odpowiadającym przeciwciałem wtórnym (1:2000) w temperaturze pokojowej przez 1,5 godziny. Immunoreaktywność zobrazowano za pomocą chromogenicznego odczynnika wykrywającego, a jądra były kontrbarwowane hematoksyliną. Analiza ilościowa była przeprowadzana za pomocą oprogramowania do analizy obrazów wymienionego w Tabeli Materiałów.
Wydajność analizy qRT-PCR
Całkowite RNA zostało wyekstrahowane za pomocą odczynnika do ekstrakcji RNA opartego na fenolach i odwrócone transkrypcji za pomocą zestawu do syntezy cDNA pierwszej nici. Transkrypcja odwrotna była wykonywana w temperaturze 25 °C przez 10 minut, 55 °C przez 15 minut oraz 85 °C przez 5 minut. następnie wykonano qRT-PCR przy użyciu systemu PCR w czasie rzeczywistym oraz mieszanki qPCR SYBR Green. Warunki amplifikacji były następujące: początkowa denaturacja w 95 °C przez 60 s, następnie 40 cykli przy 95 °C przez 15 s, 60 °C przez 15 s oraz 72 °C przez 45 s. Sekwencje primerów zostały zaprojektowane za pomocą oprogramowania do projektowania primerów i są wymienione w Tabeli 3. Jako wewnętrzny punkt odniesienia użyto GAPDH, a względne poziomy ekspresji mRNA obliczono metodą 2^-ΔΔCt . Każda próbka była analizowana w trzech egzemplarzach, a każdy eksperyment powtarzano co najmniej trzy razy. Szczegóły dotyczące odczynników, instrumentów i oprogramowania znajdują się w Tabeli Materiałów.
Wydajność analizy western blot
Tkanki jelita grubego były homogenizowane w buforze lizycznej zawierającej inhibitor proteazy na lodzie i dokładnie lizowane przez 2 godziny. Stężenia białek zostały określone za pomocą zestawu do analizy białek. Równe ilości białka były mieszane z buforem obciążającym, denaturowane przed elektroforezą, rozdzielane przez 10% SDS-PAGE i przenoszone na membrany PVDF pod napięciem 100 V przez 1 godzinę. Po blokowaniu mlekiem beztłuszczowym 5% przez 1 godzinę w temperaturze pokojowej, błony były inkubowane przez noc w temperaturze 4 °C z głównymi przeciwciałami króliczym przeciwko STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP), HSPA5 (1:1000, Cell Signaling Technology, #3177) oraz GAPDH w odpowiednim rozcieńczeniu roboczym. Po myciu TBST (20 mM Tris, 150 mM NaCl, 0,1% Tween-20; pH 7,4) błony były inkubowane z odpowiadającym przeciwciałem wtórnym (1:5000) w temperaturze pokojowej przez 1 godzinę. Pasma białkowe były wizualizowane za pomocą odczynnika wykrywającego chemiluminescencji. Intensywności pasmowe analizowano za pomocą oprogramowania do analizy obrazów, a poziomy ekspresji białek docelowych znormalizowano do GAPDH. Oczekiwane masy cząsteczkowe wykrytych białek wynosiły około 88 kDa dla STAT3, 21 kDa dla TIMP1, 28 kDa dla SPHK2, 70 kDa dla HSPA5 oraz 36 kDa dla GAPDH. Szczegóły dotyczące przeciwciał, membran, odczynników i oprogramowania znajdują się w Tabeli Materiałów. Konkretne nazwy produktów, producenci, numery katalogowe, modele instrumentów oraz szczegóły oprogramowania są podane w Tabeli Materiałów, natomiast ogólne opisy są używane w głównym tekście, aby poprawić neutralność i ogólną użyteczność.
Wykonanie analizy statystycznej
Dane przedstawiane są jako średnia ± odchylenie standardowe. Testy statystyczne przeprowadzano za pomocą oprogramowania statystycznego. Jednokierunkowa ANOVA, a następnie test post hoc Tukeya została użyta do porównań wielogrupowych. Wartość P ≤ 0,05 została uznana za statystycznie istotną.