$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Różnicowanie
Mysie komórki nabłonka dróg oddechowych pomyślnie różnicowały się po hodowli na granicy faz powietrze-ciecz z pożywką różnicującą przez 28 dni. Obecność komórek rzęskowych i kubkowych wykazano za pomocą testu immunofluorescencyjnego markera rzęsek acetylowanej α-tubuliny (zielony; Rysunek 3A) i marker komórki kubkowej Mucin5AC, odpowiednio6 (czerwony; Rysunek 3B).
Oznaczanie stężenia CSE dla różnicowania komórek
Stymulować izolowane komórki nabłonkowe różnymi stężeniami CSE przez 24 godziny. Wyniki wykazały, że żywotność komórek nabłonkowych spadała, gdy stężenie CSE było wyższe niż 6% (p <0,05, Figura 4A). W celu utrzymania aktywności komórek w warunkach długotrwałej stymulacji, pożywki hodowlane zawierające 2%, 4% i 6% CSE stosowano w celu stymulacji różnicowania komórek nabłonka przez 28 dni. Co więcej, wyniki wykazały, że wraz ze wzrostem stężenia CSE nie doszło do znaczącej śmierci komórek nabłonkowych (Figura 4B), ale nastąpiła stopniowa zmiana liczby międzykomórkowej odporności transbłonowej (TEER) (Figura 4C) i złuszczonych komórek (Rysunek 4D). Ponadto zaobserwowano ogólne zmniejszenie liczby komórek zróżnicowanych i komórek rzęskowych, gdy hodowle mysich komórek nabłonka dróg oddechowych były przewlekle narażone na CSE (Figura 4E-H). W związku z tym wybrano pożywkę różnicującą zawierającą 6% CSE, aby rozregulować funkcję bariery nabłonkowej i zmniejszyć liczbę komórek rzęskowych bez spadku żywotności komórek podczas różnicowania komórek nabłonkowych przez 28 dni.
Analiza multiomiczna
Pobrano dwie grupy komórek stymulowanych pożywką różnicową zawierającą 0% lub 6% CSE przez 28 dni. Całkowite mRNA wyekstrahowano do analizy transkryptomicznej. Po przygotowaniu biblioteki, sekwencjonowaniu, kontroli jakości, mapowaniu odczytów do genomu referencyjnego i kwantyfikacji ekspresji genów, przeanalizowano różnicową ekspresję genów między obiema grupami. Pakiet clusterProfiler R został wykorzystany do przeprowadzenia analizy wzbogacenia ontologii genów (GO) dla genów, które wykazywały zróżnicowaną ekspresję, z poprawkami dokonanymi w celu uwzględnienia rozbieżności w długości genów (Figura 5A)". Encyklopedia Genów i Genomów z Kioto (KEGG) służy jako obszerna baza danych, która ułatwia zrozumienie złożonych funkcji i ról systemów biologicznych na różnych poziomach, w tym komórkowym, organizmowym i ekologicznym, poprzez integrację danych na poziomie molekularnym pochodzących z obszernych zestawów danych, takich jak te uzyskane przez sekwencjonowanie genomu i wysokoprzepustowe metodologie eksperymentalne (Rysunek 5B). W przypadku proteomiki białka ogółem wyekstrahowano z dwóch grup komórek po 28 dniach od hodowli z pożywkami różnicowymi zawierającymi 0% lub 6% CSE. Po badaniu jakości białka, usuwaniu soli, hydrolizie enzymatycznej białek, znakowaniu peptydów Tandem Mass Tags (TMT), oddzieleniu frakcji, analizie LC-MS, identyfikacji białek i ilościowym określeniu ilości białek, można zdefiniować białka o zróżnicowanej ekspresji między dwiema grupami. Analizę GO przeprowadzono przy użyciu programu interproscan w stosunku do nieredundantnej bazy danych białek (Figura 6A), a KEGG wykorzystano do analizy szlaku (Figura 6B). Do celów metabolomicznych supernatanty hodowlane pobrano z dwóch grup komórek nabłonka dróg oddechowych po 28 dniach od hodowli z pożywką różnicującą zawierającą 0% lub 6% CSE. Po analizie HPLC-MS/MS, identyfikacji i ilościowym określeniu metabolitów można określić metabolity różnicujące między tymi dwiema grupami. Baza danych KEGG (Rysunek 7A), Baza danych metabolomu człowieka (Figura 7B) oraz Baza danych Lipidmaps (Rysunek 7C) zostały wykorzystane do opisu wszystkich metabolitów. Baza danych KEGG została wykorzystana do analizy funkcji tych zróżnicowanych metabolitów i szlaków metabolicznych (Rysunek 7D). W przypadku epigenomiki ATAC-seq wykonano zgodnie z wcześniejszymi raportami9. Jądra ekstrahowano z dwóch grup komórek po 28 dniach od hodowli z pożywką różnicową zawierającą 0% lub 6% CSE, a osad jądra zawieszono ponownie w mieszaninie reakcyjnej transpozazy Tn5. Reakcję transpozycji inkubowano w temperaturze 37 °C przez 30 minut. Po dodaniu adaptera, przygotowaniu biblioteki, ocenie jakości, grupowaniu i sekwencjonowaniu można zdefiniować epigenomikę różnicową między dwiema grupami. Analiza wzbogacenia GO genów związanych z różnicowo pikami bezpośrednio odzwierciedla rozkład liczby genów związanych z różnicami związanymi z pikami w pozycjach GO wzbogaconych w procesy biologiczne, składniki komórkowe i funkcję molekularną (Rysunek 8A). Baza danych KEGG została wykorzystana do analizy funkcji tych zróżnicowanych genów i szlaków (Rysunek 8B).
Geny kodujące białka o ścisłym połączeniu
Dane transkryptomiczne wykazały, że po długotrwałej stymulacji CSE nastąpiła zmniejszona ekspresja wielu genów kodujących białka połączeń ścisłych, w tym Cgn, Cldn2, Cldn3, Cldn8, Cldn10, Cldn20, Cldn23, Jam2 i Ocln (Figura 9A), w których zmniejszona ekspresja Cldn3 i Ocln została dodatkowo potwierdzona przez Western-Blot i immunohistochemię (IHC) (Ryciny 9B-E). Dane te sugerują, że narażenie na CSE znacznie zmniejsza ekspresję białek połączeń ścisłych komórek nabłonka dróg oddechowych.
Wyrażenie cytokin i chemokin
Oprócz rozpadu ścisłych połączeń, transkryptomika wykazała również, że po długotrwałej stymulacji CSE, ekspresja niektórych cytokin i chemokin, w tym Cxcl5, Csf3, Il1a, Il34, Ccl20 i Il33, wzrosła, ale ekspresja Ccl5 zmniejszyła się (Figura 10A). Test Luminex potwierdził ponadto, że w supernatancie mysich komórek nabłonka dróg oddechowych wystąpiły zwiększone stężenia CXCL-5, CSF-3 i IL-1α (Rycina 10B-D), ale zmniejszone CCL-5 (Figura 10E). Wyniki te sugerują, że ekspozycja na CSE selektywnie wpływa na ekspresję cytokin i chemokin.

Rysunek 1: Wzrost pierwotnych komórek nabłonka mysich dróg oddechowych. (A) Nowo wyizolowane komórki wydają się okrągłe i przezroczyste. (B) Do 2 dnia po izolacji, rozszerzone komórki uformowały się w małe wyspy komórkowe. (C) Do 4 dnia, tworzenie się większych wysp komórkowych. (D) Do 6 dnia ponad 90% fuzji z morfologią kostki brukowej. (4-krotne powiększenie; podziałka liniowa: 250 μm). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Test immunofluorescencyjny markerów komórek nabłonkowych. Ekspresja cytokeratyny w pierwotnych mysich komórkach nabłonka dróg oddechowych, A549 (kontrola pozytywna) i Raw264.7 (kontrola ujemna). Jądra barwione DAPI (powiększenie 63x; podziałka skali: 40 μm). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Test immunofluorescencyjny dla markerów komórek rzęskowych i kubkowych. Stosowanymi markerami były acetylowana (A) α-tubulina i (B) mucyna 5AC do potwierdzania odpowiednio komórek rzęsek i komórek kubkowych (powiększenie 100; podziałka skali: 10 μm w panelu A i 50 μm w panelu B). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Wpływ stężeń CSE na komórki nabłonkowe. (A) Żywotność komórek po 24 godzinach przy różnych stężeniach CSE. (B) Żywotność komórek po 28 dniach. (C) Pomiar TEER. (D) Liczba złuszczonych komórek. (E-H) Immunofluorescencja dla komórek rzęskowych (zielonych) i kubkowych (fioletowych) przy stężeniach 0%, 2%, 4% i 6% CSE (powiększenie 20x; podziałka skali: 50 μm). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Analizy wzbogacenia GO i KEGG pod kątem ekspresji genów. (A) Znaczenie terminu GO we wzbogaconym procesie biologicznym, komponencie komórkowym i funkcji molekularnej. (B) Odsetek genów o zróżnicowanej ekspresji w szlakach KEGG. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Analiza ekspresji białek. (A) Wykres słupkowy wzbogacający GO pokazujący różnicową ekspresję białka. (B) Analiza szlaku KEGG w oparciu o białka o zróżnicowanej ekspresji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Analiza szlaku metabolitów KEGG. (A) Adnotacja wszystkich metabolitów za pomocą bazy danych KEGG. (B) Opisywanie wszystkich metabolitów przy użyciu bazy danych metabolomu ludzkiego. (C) Opisywanie wszystkich metabolitów za pomocą bazy danych Lipidmaps. (D) Dystrybucja metabolitów o zróżnicowanej ekspresji w różnych szlakach KEGG. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 8: Analiza ekspresji genów związanych z pikami różnicowymi. (A) Wykres słupkowy wzbogacania GO pokazujący różnicową ekspresję białek. (B) Analiza szlaku KEGG w oparciu o białka o zróżnicowanej ekspresji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 9: Ekspresja białka ścisłego połączenia w komórkach nabłonka mysich dróg oddechowych ze stymulacją lub bez stymulacji CSE. (A) Mapa cieplna ekspresji białka połączenia ścisłego w transkryptomie. (B-E) Poziomy ekspresji (B,C) Cldn3 i (D,E) Ocln w mysich komórkach nabłonka dróg oddechowych wykryto za pomocą Western Blot i immunohistochemii dla Cldn3 i Ocln (powiększenie 80; paski skali: 50 μm). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 10: Analiza ekspresji cytokin. (A) Mapa cieplna transkryptomu ekspresji cytokin. (B-E) Test Luminex dla poziomów CXCL-5, CSF-3, IL-1α i CCL-5 w komórkach z/bez stymulacji CSE. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.