$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Trening AI w rozpoznawaniu rzęsek
Pomiar i ocena długości strukturalnej i składu rzęsek może być żmudnym, czasochłonnym i podatnym na błędy procesem. Tutaj używamy Ai do dokładnej segmentacji rzęsek z dużej puli obrazów i analizujemy ich długość i intensywność za pomocą narzędzia analitycznego (Rysunek 1). Wszystkie podejścia do sztucznej inteligencji wymagają etapów szkoleniowych w celu ich wdrożenia. Ustaliliśmy proces treningowy rozpoznawania rzęsek, który został przeprowadzony poprzez ręczne nakładanie masek binarnych na struktury rzęskowe. Informacje te są następnie wykorzystywane do trenowania sztucznej inteligencji na podstawie charakterystyki pikseli w zastosowanych plikach binarnych. Ogólnie rzecz biorąc, szkolenie obejmuje oprogramowanie przechodzące przez kilka iteracji, około 1000, i jest uważane za optymalne, jeśli strata lub wskaźnik błędów uczenia jest mniejszy niż 1%. Jednak liczba iteracji i błędów w procesie trenowania może się różnić w zależności od przykładowych obrazów używanych do trenowania. Na przykład, po naszych sesjach treningowych z wykorzystaniem obrazów rzęsek neuronalnych in vitro, poziom błędu wynosił 1,378% w porównaniu z 3,36% w przypadku obrazów przekrojów mózgu in vivo (rysunek uzupełniający 1). Po zakończeniu trenowania sztuczna inteligencja może być następnie wykorzystana do segmentacji rzęsek z obrazów eksperymentalnych w ciągu kilku sekund, a powstałe w ten sposób maski binarne są wykorzystywane do pomiaru parametrów strukturalnych. Eliminuje to konieczność segmentacji obiektów przy użyciu tradycyjnej metody progowania intensywności, co może być trudne w przypadku obrazów o wysokim poziomie szumów tła lub gdy obiekty znajdują się w bliskiej odległości od siebie. Sztuczna inteligencja zmniejsza również ryzyko błędów i stronniczości, stosując ten sam algorytm do wszystkich obrazów, niezależnie od użytkownika.
Pomiar długości rzęsek za pomocą GA3
Długość rzęsek jest ściśle regulowana i wiąże się z funkcjonalnym wpływem na sygnalizację rzęskową16,19. W tym przypadku zmierzyliśmy długości rzęsek za pomocą potoku analizy w oprogramowaniu NIS Elements o nazwie General Analysis 3 lub GA3. GA3 pomaga w łączeniu wielu narzędzi w jednym przepływie pracy w celu tworzenia niestandardowych procedur dla każdego eksperymentu. Zaczęliśmy od zmierzenia długości rzęsek w linii komórkowej. Rzęski na komórkach wewnętrznego kanału zbiorczego rdzenia myszy (IMCD-3) znakowano immunologicznie acetylowaną tubuliną i obrazowano za pomocą mikroskopu konfokalnego. Zmierzyliśmy długości rzęsek za pomocą GA3 po segmentacji za pomocą segment.ai (rysunek uzupełniający 3A). Podczas gdy acetylowana α-tubulina znajduje się preferencyjnie w rzęsce pierwotnej, znajduje się również w innych regionach bogatych w mikrotubule, takich jak cytoszkielet, a także mostek cytokinekinetyczny. Wytrenowana sztuczna inteligencja prawidłowo zidentyfikowała rzęski na obrazie, ale nie inne nierzęskowe, acetylowane struktury dodatnie tubuliny. Rzęski na komórkach IMCD wahały się od 0,5 μm do 4,5 μm o średniej długości 1,8 ± 0,04 μm (Ryc. 2A). Następnie przetestowaliśmy zdolność sztucznej inteligencji do pomiaru długości rzęsek w pierwotnych kulturach neuronalnych. Hodowaliśmy neurony z podwzgórza i hipokampa nowonarodzonych myszy przez 10 dni i znakowaliśmy je immunologicznie markerem rzęsek cyklazy adenylanowej III (ACIII)21,41. Analizując hodowle neuronów, stwierdziliśmy, że przydatne jest zastosowanie filtra przed analizą statystyczną długości. Ze względu na niższy stosunek sygnału do szumu zidentyfikowano kilka obiektów mniejszych niż 1 μm, które nie były rzęskami. Dlatego przefiltrowaliśmy dane, aby wyeliminować wszelkie obiekty o długości mniejszej niż 1 μm, aby upewnić się, że analizowane są tylko rzęski. W hodowanych neuronach podwzgórza długość rzęsek wahała się od 2 μm do 7 μm przy średniej długości 3,8 ± 0,19 μm (Ryc. 2B). Co ciekawe, hodowane rzęski neuronalne hipokampa były dłuższe i miały średnią długość 6,73 ±0,15 μm (Rysunek 2C). Doniesiono, że różne jądra neuronalne w podwzgórzu wykazują różne długości rzęsek i że te rzęski zmieniają swoją długość w odpowiedzi na zmiany fizjologiczne w sposób specyficzny dla jądra19,23. Dlatego oznaczyliśmy również skrawki mózgu podwzgórza od dorosłych samców myszy C57BL / 6J z ACIII i zobrazowaliśmy jądro łukowate (ARC) i jądro przykomorowe (PVN). Używając GA3 do pomiaru długości rzęsek, zaobserwowaliśmy, że rzęski podwzgórza in vivo wydawały się dłuższe niż rzęski in vitro. W szczególności rzęski podwzgórza in vivo wahają się od 1 μm do około 15 μm (Ryc. 3). Nie stwierdzono istotnych różnic między długościami rzęsek w PVN (5,54 ± 0,0,42 μm) a tymi w ARC (6,16 ± 0,27 μm) (Rysunek 3C)23. Podobnie, rzęski w obszarze rogu-amonisa (CA1) hipokampa mają węższy zakres długości od 1 μm do 10 μm ze średnią długością 5,28 ± 0,33 μm (Rysunek 3). Zgodnie z wcześniej opublikowanymi badaniami, nasza analiza z wykorzystaniem narzędzi Ai i GA3 wykazała, że rzęski z różnych regionów mózgu wykazują różnorodność długości19,23. Co więcej, korzystając z tego podejścia AI, jesteśmy w stanie szybko ocenić dużą liczbę rzęsek.
Pomiar składu rzęsek za pomocą GA3
Rzęska pierwotna jest węzłem sygnalizacyjnym dla wielu szlaków, które wykorzystują różne rodzaje białek do wykonywania unikalnych funkcji, takich jak białka motoryczne, białka transportu wewnątrzwiciowego i GPCR, aby wymienić tylko kilka3,24,42,43. Utrzymanie odpowiedniego poziomu tych białek w rzęsce jest ważne dla prawidłowego funkcjonowania i często wydaje się zależne od kontekstu komórkowego. Znakowanie fluorescencyjne tych białek nie tylko umożliwiło nam ich wizualizację, ale także ilościowe określenie ich intensywności jako miary ilości znakowanego białka w stosunkowo małym przedziale20. Dlatego staraliśmy się określić intensywność rzęskowego GPCR, receptora hormonu koncentracji melaniny 1 (MCHR1), in vivo zarówno w ARC, jak i PVN podwzgórza dorosłych samców myszy24,44. Korzystając z Ai i GA3, zmierzyliśmy długości rzęsek dodatnich MCHR1 wraz z intensywnością, aby upewnić się, że liczone obiekty były rzęskami (rysunek uzupełniający 3A). Wyeliminowaliśmy obiekty po analizie, które miały długość mniejszą niż 2 μm i przeanalizowaliśmy intensywności pozostałych masek binarnych. Co ciekawe, odkryliśmy, że intensywność rzęskowego MCHR1 w PVN jest znacznie wyższa niż w ARC, co wskazuje na silniejszą obecność rzęskowego MCHR1 w PVN (Figura 4). Konieczne są dalsze badania w celu określenia znaczenia rzęskowego MCHR1 w tych obwodach neuronalnych. Zmierzyliśmy również intensywność rzęskowego MCHR1 w pierwotnych hodowanych neuronach podwzgórza i hipokampa. Rzęski z obu kultur wykazują szeroki rozkład intensywności MCHR1, co sugeruje obecność heterogenicznych populacji neuronów (ryc. uzupełniająca 2). Tak więc stosowanie zaawansowanych narzędzi analitycznych, takich jak Ai i GA3, pozwala na ocenę niejednorodności rzęsek w tej samej tkance lub między wieloma tkankami. Interesujące będzie sprawdzenie, czy inne neuronalne GPCR wykazują podobne różnice w swojej lokalizacji w neuronach tej samej tkanki i czy zmienia się to w odpowiedzi na zmiany fizjologiczne.
Kolokalizacja
Podczas gdy pomiar intensywności fluorescencji w całym polu obrazu może dać wrażenie białka, nie dostarcza informacji, takich jak rozmieszczenie przestrzenne lub bliskość innych pobliskich białek i struktur komórkowych. Tutaj zmierzyliśmy nakładanie się MCHR1 z ACIII jako markerem rzęsek, wykreślając intensywność MCHR1 w stosunku do intensywności ACIII dla każdej maski binarnej (Rysunek 5). Wykres pokazuje, że większość rzęsek jest dodatnia zarówno dla ACIII, jak i MCHR1, chociaż niektóre rzęski wykazują silniejszą ekspresję jednego kanału w stosunku do drugiego. Ponadto istnieją rzęski, które wykazują obecność ACIII lub MCHR1, co jest oczywiste w punktach, które leżą bezpośrednio odpowiednio na osi x i y. Aby określić ilościowo to nakładanie się, zmierzyliśmy współczynnik nakładania się Mandera i porównaliśmy zakres ekspresji MCHR1 w rzęskach neuronalnych ARC i PVN40. Co ciekawe, nasza analiza wykazała, że nastąpił znaczny wzrost współczynników PVN (0,6382 ± 0,0151) niż tych w ARC (0,5430 ± 0,0181) (Rysunek 5C). Jest to zgodne z naszymi poprzednimi danymi, w których zaobserwowaliśmy wyższą intensywność MCHR1 w PVN w porównaniu z ARC (Rysunek 4). Dane te sugerują, że podobnie jak długość rzęsek, wzorzec ekspresji MCHR1 w przedziale rzęskowym różni się w różnych obszarach mózgu. Korzystając z tego samego toku analizy, możliwe będzie określenie, czy inne rzęskowe GPCR, takie jak receptor neuropeptydu Y typu 2 (NPY2R) i receptor somatostatyny typu 3 (SSTR3) wykazują podobną różnorodność.
Pomiar profilu intensywności wzdłuż rzęski
Po zidentyfikowaniu rzęsek za pomocą segment.ai, receptura GA3 może zostać zmodyfikowana, aby połączyć analizę rzęsek z identyfikacją innych struktur interesujących na obrazie. Na przykład znakowanie za pomocą podstawowych markerów ciała jest przydatne do identyfikacji polaryzacji rzęsek. Aby przeprowadzić tę analizę, zobrazowaliśmy podwzgórzowe skrawki mózgu myszy P0, które wyrażają ARL13B-mCherry i Centrin2-GFP oraz zobrazowaliśmy ARC i PVN34. W tym przypadku rzęski zostały zidentyfikowane za pomocą Ai, tak jak wcześniej, ale teraz zmodyfikowana receptura GA3 obejmuje identyfikację Centrin2-GFP, białka dośrodkowego znajdującego się u podstawy rzęsek (rysunek uzupełniający 3B). Znakując Centrin2-GFP, podstawę rzęsek można odróżnić od końcówek rzęsek dodatnich ARL13B-mCherry (Rysunek 6A). Następnie, zamiast mierzyć intensywność w obrębie całych rzęsek, jesteśmy w stanie zmierzyć zmiany intensywności ARL13B na całej długości rzęsek (Rysunek 6B). Możemy również porównać różnice w intensywności ARL13B między końcami bliższymi a dalszymi końcami rzęsek. Aby to zrobić, podzieliliśmy długość rzęski na pojemniki o wielkości 1 mikrona, zaczynając od podstawy i wyznaczyliśmy pierwszy pojemnik mikronowy jako koniec proksymalny i ostatni pojemnik mikronowy jako koniec dystalny. Nasza analiza wykazała, że zarówno w ARC, jak i PVN, zarówno w ARC, jak i PVN, jest znacznie więcej ARL13B bliżej podstawy niż na czubku rzęski, co jest zgodne z wcześniej opublikowanymi badaniami na ludzkich chondrocytach45 (Figura 6C). W tego typu analizie, zamiast stosować filtr długości w celu wykluczenia z analizy małych obiektów nierzęskowych, analizowane są tylko rzęski związane z znakowaniem Centrin2-GFP. Może to być korzystne w sytuacjach, gdy mutacje genetyczne powodują, że rzęski są bardzo krótkie lub jeśli zaangażowane są zmiany w subdomenach rzęsek, takich jak strefa przejściowa lub wierzchołek. Identyfikacja rzęsek za pomocą analizy Ai i GA3 jest wysoce elastyczna i może być dostosowana do różnych złożonych pytań badawczych.

Rysunek 1. Przepływ pracy do pomiaru długości i intensywności rzęsek za pomocą Ai.(A) Aby wytrenować Ai, pliki binarne są rysowane wokół obiektów zainteresowania (rzęsek) na surowych obrazach treningowych. Korzystając z narysowanych plików binarnych, Segment Ai jest szkolony w rozpoznawaniu kształtu i intensywności pikseli rzęsek. (B) Następnie wytrenowany Segment Ai jest stosowany do surowych obrazów eksperymentalnych. Rysuje pliki binarne na obiektach, które rozpoznaje jako rzęski. Te pliki binarne można udoskonalić, aby upewnić się, że analizowane są wszystkie i tylko rzęski. (C) Program GA3 jest skonstruowany do analizy intensywności i długości obiektów rozpoznawanych przez sztuczną inteligencję. (D) Rekordy są importowane do tabeli w oprogramowaniu. Tabelę tę można następnie wyeksportować do dalszej analizy. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2. Pomiary długości rzęsek in vitro. Reprezentatywne obrazy rzęsek w (A) komórkach IMCD (zielony, acetylowana tubulina) (B) pierwotnych kulturach podwzgórza (zielony, ACIII) i (C) kulturach hipokampa (zielony, ACIII). Do rozpoznania rzęsek użyto wytrenowanej sztucznej inteligencji, jak pokazano na masce binarnej (magenta), a następnie GA3 użyto do pomiaru długości rzęsek. Rozkład długości rzęsek jest przedstawiony na wykresie jako procent rzęsek w pojemnikach 0,5 lub 1,0 mikrona. * wskazuje na mostek cytokinetyczny prawidłowo nierozpoznany przez Ai. n=225 rzęsek w komórkach IMCD z 3 powtórzeń, 54 rzęski w podwzgórzu i 139 rzęsek w hodowlach hipokampa od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3. Pomiary długości rzęsek in vivo. (A) Reprezentatywne obrazy rzęsek (zielony, ACIII) w ARC, PVN i CA1 sekcji mózgu dorosłej myszy. (B) Do rozpoznania rzęsek użyto wytrenowanej sztucznej inteligencji w elementach NIS, jak pokazano na masce binarnej (magenta), a następnie użyto GA3 do pomiaru długości rzęsek. (C) Rozkład długości rzęsek jest przedstawiony na wykresie jako procent rzęsek w pojemnikach jednomikronowych. n = 68 rzęsek w ARC, 36 w PVN i 29 w CA1 od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4. Pomiary intensywności barwienia rzęsek rzęsek neuronalnych w podwzgórzu wspomagane sztuczną. (A) Reprezentatywne obrazy rzęsek (MCHR1, czerwony) w ARC i PVN sekcji mózgu dorosłej myszy. Do rozpoznania rzęsek wykorzystano wytrenowaną sztuczną inteligencję w elementach NIS, jak pokazano w masce binarnej (cyjan), a następnie GA3 użyto do pomiaru intensywności barwienia MCHR1 w rzęskach. (B) Intensywności MCHR1 są przedstawione na wykresie jako średnia ± S.E.M. Każda kropka reprezentuje rzęskę. * p < 0,05, test t-studenta. (C) Rozkład intensywności MCHR1 jest przedstawiony na wykresie jako procent rzęsek w przedziałach 0,2 x 107 dowolnych jednostek (A). n = 53 rzęski w ARC, 78 w PVN od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5. Analiza kolokalizacji rzęsek wspomagana sztuczną inteligencją. (A, B) Reprezentatywne obrazy rzęsek odpowiednio w ARC i PVN. Rzęski są oznaczone ACIII (na zielono) i MCHR1 (na czerwono). Do rozpoznania rzęsek wykorzystano wytrenowaną sztuczną inteligencję w elementach NIS, jak pokazano na masce binarnej (magenta dla rzęsek znakowanych ACIII, cyjan dla rzęsek oznaczonych MCHR1). GA3 był używany do rozpoznawania rzęsek, które zawierały zarówno ACIII, jak i MCHR1. (C) Wartości współczynnika nakładania się Manders (MOC) są przedstawione na wykresie jako średnia ± S.E.M. Każda kropka reprezentuje rzęskę. * p < 0,05, test t-studenta. (D) Wykres punktowy intensywności MCHR1 w funkcji intensywności ACIII w ARC i PVN. Każda kropka reprezentuje rzęskę. n = 72 rzęski w ARC, 47 w PVN od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6. Analiza rzęsek i ciała podstawnego. (A) Reprezentatywne obrazy rzęsek (czerwony, ARL13B-mCherry) i markera podstawy ciała (zielony, Centrin2-GFP) w ARC i PVN myszy P0. Wytrenowana sztuczna inteligencja została użyta do rozpoznania rzęsek, jak pokazano w masce binarnej (cyjan). Maska binarna dla ciała podstawowego (magenta) została narysowana przez progowanie w przepisie GA3. (B) Reprezentatywna intensywność skanowania linii rzęski. (C) Natężenie ARL13B na bliższych i dalszych końcach Ai zidentyfikowało rzęski przedstawione na wykresie jako średnia ± S.E.M. Końce proksymalne i dystalne definiuje się jako obszar o długości 1 μm i długości ostatniej 1 μm odpowiednio od podstawy rzęski. Każda kropka reprezentuje rzęskę. * p < 0,05. n = 6 rzęsek w ARC od 2 zwierząt i 21 rzęsek w PVN od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Rysunek uzupełniający 1. Wykresy strat związanych z trenowaniem AI. (A, B) Wykresy przedstawiające utratę segment.ai treningowej na rzęskach neuronalnych odpowiednio in vitro i in vivo. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający 2. Pomiary intensywności barwienia rzęsek neuronalnych in vitro wspomagane sztuczną intesencją. (A, B) Reprezentatywne obrazy rzęsek (MCHR1, czerwony) odpowiednio w pierwotnych kulturach podwzgórza i hipokampa. Do rozpoznania rzęsek wykorzystano wytrenowaną sztuczną inteligencję w elementach NIS, jak pokazano w masce binarnej (cyjan), a następnie GA3 użyto do pomiaru intensywności barwienia MCHR1 w rzęskach. Rozkład intensywności MCHR1 jest przedstawiony na wykresie jako procent rzęsek w 1000 przedziałach AU dla posiewów podwzgórza i 2000 AU dla kultur hipokampa. n = 30 rzęsek w podwzgórzu i 106 rzęsek w kulturach hipokampa od 3 zwierząt. Podziałka 10 μm. Kliknij tutaj, aby pobrać ten plik.
Rysunek uzupełniający 3. Analiza ogólna 3 przepisy na analizę rzęsek. (A) Prosta Analiza Ogólna (GA3) receptura do pomiaru długości rzęsek, intensywności i współczynnika Mandera. (B) Złożona receptura GA3 do pomiaru intensywności wzdłuż długości rzęski przy użyciu markera dla ciała podstawnego. Kliknij tutaj, aby pobrać ten plik.