Różnicowanie
Mysie komórki nabłonka dróg oddechowych uległy pomyślnemu różnicowaniu po hodowli na granicy faz powietrze-ciecz w medium do różnicowania przez 28 dni. Obecność komórek rzęskowych i kubkowych wykazano odpowiednio za pomocą testu immunofluorescencyjnego z użyciem markera rzęsek acetylowanej α-Tubuliny (kolor zielony; Rysunek 3A) oraz markera komórek kubkowych Mucin5AC6 (kolor czerwony; Rysunek 3B).
Wyznaczenie stężenia CSE dla różnicowania komórek
Stymulowano wyizolowane komórki nabłonkowe różnymi stężeniami CSE przez 24 h. Wyniki wykazały, że żywotność komórek nabłonkowych spadła, gdy stężenie CSE było wyższe niż 6% (p <0,05, Rysunek 4A). Aby zachować aktywność komórek podczas długotrwałej stymulacji, do stymulowania różnicowania komórek nabłonkowych przez 28 dni wykorzystano pożywki hodowlane zawierające 2%, 4% i 6% CSE. Ponadto wyniki wykazały, że wraz ze wzrostem stężenia CSE nie odnotowano istotnej śmierci komórek nabłonkowych (Rysunek 4B), ale nastąpiła stopniowa zmiana wartości oporu przezbłonowego (TEER) (Rysunek 4C) oraz liczby złuszczonych komórek (Rysunek 4D). Dodatkowo odnotowano ogólną redukcję liczby komórek zróżnicowanych i komórek rzęsatych w przypadku przewlekłej ekspozycji kultur komórek nabłonkowych dróg oddechowych myszy na CSE (Rysunek 4E-H). W związku z tym, do dysregulacji funkcji bariery nabłonkowej i zmniejszenia liczby komórek rzęsatych bez obniżenia żywotności komórek podczas 28-dniowego różnicowania komórek nabłonkowych, wybrano pożywkę do różnicowania zawierającą 6% CSE.
Analiza multiomiczna
Zabrano dwie grupy komórek stymulowanych różnymi mediami zawierającymi 0% lub 6% CSE przez 28 dni. Do analizy transkryptomicznej wyekstrahowano całkowite mRNA. Po przygotowaniu biblioteki, sekwencjonowaniu, kontroli jakości, mapowaniu odczytów do genomu referencyjnego i kwantyfikacji ekspresji genów przeanalizowano różnicową ekspresję genów pomiędzy dwiema grupami. Do przeprowadzenia analizy wzbogacenia Gene Ontology (GO) dla genów wykazujących różnicową ekspresję wykorzystano pakiet R clusterProfiler, wprowadzając korekty w celu uwzględnienia rozbieżności w długości genów (Ryc. 5A). Kyoto Encyclopedia of Genes and Genomes (KEGG) służy jako kompleksowa baza danych ułatwiająca zrozumienie funkcjonalności i ról złożonych systemów biologicznych na różnych poziomach, w tym komórkowym, organizmalnym i ekologicznym, poprzez integrację danych na poziomie molekularnym pochodzących z rozległych zbiorów danych, takich jak te wytworzone w drodze sekwencjonowania genomicznego i wysokoprzepustowych metod eksperymentalnych (Ryc. 5B). W przypadku proteomiki, całkowite białka wyekstrahowano z dwóch grup komórek w 28. dniu po hodowli w różnych mediach zawierających 0% lub 6% CSE. Po teście jakości białek, usunięciu soli, enzymatycznej hydrolizie białek, znakowaniu peptydów za pomocą Tandem Mass Tags (TMT), rozdzieleniu frakcji, analizie LC-MS, identyfikacji białek i kwantyfikacji białek, zdefiniowano białka o różnicowej ekspresji pomiędzy dwiema grupami. Analizę GO przeprowadzono przy użyciu programu interproscan w odniesieniu do nieredundantnej bazy danych białek (Ryc. 6A), a KEGG wykorzystano do analizy szlaków (Ryc. 6B). W przypadku metabolomiki, nadnadrobia hodowlane zebrano z dwóch grup nabłonkowych komórek dróg oddechowych w 28. dniu po hodowli w medium różnicującym zawierającym 0% lub 6% CSE. Po analizie HPLC-MS/MS, identyfikacji i kwantyfikacji metabolitów, zdefiniowano metabolity różnicowe pomiędzy dwiema grupami. Do adnotacji wszystkich metabolitów wykorzystano bazę danych KEGG (Ryc. 7A), Human Metabolome Database (Ryc. 7B) oraz bazę danych Lipidmaps (Ryc. 7C). Bazę danych KEGG wykorzystano do analizy funkcji tych różnicowych metabolitów oraz szlaków metabolicznych (Ryc. 7D). W przypadku epigenomiki, przeprowadzono ATAC-seq zgodnie z wcześniejszymi doniesieniami9. Jądra komórkowe wyekstrahowano z dwóch grup komórek w 28. dniu po hodowli w medium różnicującym zawierającym 0% lub 6% CSE, a osad z jąder resuspenderowano w mieszaninie reakcyjnej transpozazy Tn5. Reakcję transpozycji inkubowano w 37 °C przez 30 min. Po dodaniu adaptera, przygotowaniu biblioteki, ocenie jakości, klastrowaniu i sekwencjonowaniu, zdefiniowano różnicową epigenomikę pomiędzy dwiema grupami. Analiza wzbogacenia GO genów powiązanych z różnicowymi pikami bezpośrednio odzwierciedla rozkład liczby genów powiązanych z różnicowymi pikami w elementach GO wzbogaconych w procesy biologiczne, komponenty komórkowe i funkcje molekularne (Ryc. 8A). Bazy danych KEGG użyto do analizy funkcji tych różnicowych genów i szlaków (Ryc. 8B).
Geny kodujące białka połączeń ścisłych
Dane transkrypcyjne wykazały, że po długotrwałej stymulacji CSE nastąpił spadek ekspresji licznych genów kodujących białka połączeń ścisłych, w tym Cgn, Cldn2, Cldn3, Cldn8, Cldn10, Cldn20, Cldn23, Jam2 i Ocln (Rysunek 9A), przy czym obniżona ekspresja Cldn3 i Ocln została dodatkowo potwierdzona za pomocą Western-Blot oraz immunohistochemii (IHC) (Rysunki 9B-E). Dane te sugerują, że ekspozycja na CSE znacznie zmniejsza ekspresję białek połączeń ścisłych w komórkach nabłonkowych dróg oddechowych.
Ekspresja cytokin i chemokin
Oprócz rozpadu połączeń ścisłych, analiza transkryptomiczna wykazała również, że po długotrwałej stymulacji CSE ekspresja niektórych cytokin i chemokin, w tym Cxcl5, Csf3, Il1a, Il34, Ccl20 oraz Il33, wzrosła, natomiast ekspresja Ccl5 uległa zmniejszeniu (Rysunek 10A). Test Luminex dodatkowo potwierdził zwiększone stężenia CXCL-5, CSF-3 i IL-1α (Rysunek 10B-D) oraz zmniejszone stężenie CCL-5 (Rysunek 10E) w supernatancie mysich komórek nabłonka dróg oddechowych. Wyniki te sugerują, że ekspozycja na CSE selektywnie wpływa na ekspresję cytokin i chemokin.

Rysunek 1: Wzrost pierwotnych komórek nabłonka dróg oddechowych myszy. (A) Nowo wyizolowane komórki mają okrągły kształt i są przezroczyste. (B) W 2. dniu po izolacji namnożone komórki tworzą małe wyspy komórkowe. (C) W 4. dniu obserwuje się tworzenie większych wysp komórkowych. (D) W 6. dniu następuje fuzja powyżej 90% komórek z morfologią brukową (4x powiększenie; paski skali: 250 µm). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 2: Analiza immunofluorescencyjna markerów komórek nabłonkowych. Ekspresja cytokeratyny w pierwotnych komórkach nabłonkowych dróg oddechowych myszy, A549 (kontrola pozytywna) oraz Raw264.7 (kontrola negatywna). Jądra komórkowe zabarwione DAPI (powiększenie 63x; paski skali: 40 µm). Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 3: Analiza immunofluorescencyjna markerów komórek rzęskowych i kubkowych. Wykorzystane markery to acetylowana (A) α-Tubulina oraz (B) Mucyna 5AC w celu potwierdzenia obecności odpowiednio komórek rzęskowych i komórek kubkowych (powiększenie 100; paski skali: 10 µm na panelu A i 50 µm na panelu B). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Wpływ stężeń CSE na komórki nabłonkowe. (A) Żywotność komórek po 24 h 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ęsatych (zielone) i komórek kubkowych (fioletowe) przy stężeniach CSE 0%, 2%, 4% i 6% (powiększenie 20x; paski skali: 50 µm). Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Analizy wzbogacenia GO i KEGG dla ekspresji genów. (A) Istotność terminów GO we wzbogaconych procesach biologicznych, komponentach komórkowych i funkcjach molekularnych. (B) Proporcja genów różnicowo wyrażonych w szlakach KEGG. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 6: Analiza ekspresji białek. (A) Wykres słupkowy wzbogacenia GO przedstawiający różnicową ekspresję białek. (B) Analiza ścieżek KEGG na podstawie różnicowo eksponowanych białek. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 7: Analiza szlaków metabolicznych KEGG dla metabolitów. (A) Adnotowanie wszystkich metabolitów przy użyciu bazy danych KEGG. (B) Adnotowanie wszystkich metabolitów przy użyciu Human Metabolome Database. (C) Adnotowanie wszystkich metabolitów przy użyciu bazy danych Lipidmaps. (D) Rozkład różnicowo wyrażonych metabolitów w różnych szlakach KEGG. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rycina 8: Analiza ekspresji genów związanych z różnicami w peakach. (A) Wykres słupkowy wzbogacenia GO przedstawiający różnicową ekspresję białek. (B) Analiza szlaków KEGG oparta na różnicowo eksponowanych białkach. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 9: Ekspresja białek połączeń ścisłych w komórkach nabłonka dróg oddechowych myszy z lub bez stymulacją CSE. (A) Mapa cieplna ekspresji białek połączeń ścisłych w transkryptomie. (B-E) Poziomy ekspresji (B,C) Cldn3 oraz (D,E) Ocln w komórkach nabłonka dróg oddechowych myszy zostały określone metodą Western Blot oraz immunohistochemii dla Cldn3 i Ocln (powiększenie 80x; paski skali: 50 µm). Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rysunek 10: Analiza ekspresji cytokin. (A) Mapa ciepła transkryptomu ekspresji cytokin. (B-E) Analiza Luminex dla poziomów CXCL-5, CSF-3, IL-1α oraz CCL-5 w komórkach ze stymulacją CSE i bez niej. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.