KAV generuje wizualne reprezentacje struktury sieci
Kontrast pomiędzy ECFC a tłem macierzy na obrazach w kontraście fazowym umożliwia programowi KAV identyfikację struktur specyficznych dla komórek. Sieci ECFC zidentyfikowane przez KAV są przedstawione graficznie w formie szkieletu oraz maski, aby zilustrować struktury wykorzystywane przez oprogramowanie do kwantyfikacji (Rysunek 2A). Co istotne, jakościowa ocena reprezentacji szkieletu i maski pozwala na szybką weryfikację czułości i dokładności KAV, co może być pomocne przy ustalaniu optymalnych ustawień progowych oraz interpretacji wyników analizy. Gdy jakość obrazów w kontraście fazowym jest wysoka i osiągnięto wystarczający kontrast, KAV precyzyjnie identyfikuje sieci ECFC, o czym świadczy podobieństwo pomiędzy obrazem w kontraście fazowym a wygenerowanym przez KAV szkieletem i maską (Rysunek 2A oraz Wideo 1).
Alternatywnie, jeśli obrazy w kontraście fazowym nie charakteryzują się wysokim kontrastem i/lub występują artefakty obrazowania, takie jak prążkowanie, dokładność detekcji sieci ulega zmniejszeniu, a wyniki stają się niejednoznaczne (Rycina 2B). Ponadto powstanie pęcherzyków powietrza w pożywce komórkowej może również utrudniać dokładność detekcji (Rycina 2). Problemy z jakością obrazu można jednak często rozwiązać poprzez wybór innych metod progowania dostępnych w oprogramowaniu KAV. Na przykład obraz w kontraście fazowym na Rycynie 2B został przeanalizowany przy użyciu identycznych ustawień przetwarzania obrazu, z wyjątkiem metody progowania. Z wersji szkieletowej (Skeleton) i maski (Mask) wynika, że progowanie metodą Otsu doprowadziło do dokładniejszej detekcji sieci ECFC widocznych na obrazie w kontraście fazowym (Rycina 2B). Zatem jakość obrazu ma kluczowe znaczenie dla uzyskania dokładnych i znaczących wyników w tym teście. Niemniej jednak różne metody progowania zawarte w interfejsie użytkownika KAV umożliwiają dostosowanie analizy obrazu w zależności od jakości obrazów wejściowych.
Mikroskopia poklatkowa identyfikuje różnice jakościowe i ilościowe w waskulogenezie ECFC po płodowej ekspozycji na cukrzycę ciążową
Niedawno podejście analityczne KAV zostało zastosowane do oceny funkcji płodowych ECFC po ekspozycji na macierną cukrzycę typu 2 (T2DM) in utero16. Dzięki zastosowaniu KAV zidentyfikowano zmienioną kinetykę waskulogenezy w płodowych ECFC eksponowanych na T2DM. Jednakże, poza ekspozycją na T2DM, która występuje przez cały okres ciąży, cukrzyca ciążowa (GDM), czyli rozwój nietolerancji glukozy zazwyczaj w trzecim trymestrze ciąży, również upośledza funkcję ECFC13. W związku z tym zastosowano KAV, aby ustalić, czy ECFC eksponowane na GDM również wykazują zmienioną kinetykę tworzenia sieci ECFC. Fazę 1 (0-5 h) i Fazę 2 (5+ h) oceniono za pomocą mikroskopii poklatkowej w połączeniu z analizą KAV (Rycina 3). Reprezentatywne obrazy w kontraście fazowym, pozyskane na początku akwizycji obrazu (0.50 h) oraz w trakcie przebiegu czasowego (5.00 i 10.00 h), przedstawiają tworzenie sieci ECFC w czterech próbkach badanych w jednym dniu eksperymentalnym (Rycina 3 oraz Wideo 1). Pomimo równoważnego zasiewu komórek, co zaobserwowano na obrazach w kontraście fazowym w 0.50 h, różnice jakościowe w strukturze sieci są widoczne w punktach czasowych 5.00 i 10.00 h. ECFC z niepowikłanej ciąży (UC) tworzą złożoną i gęstą sieć 5 godzin po wysianiu, co jest zgodne z naszymi wcześniej opublikowanymi danymi16. Przeciwnie, ECFC z próbki GDM 1 (GDM1) tworzą bardzo niewiele struktur sieciowych, które nie są ze sobą połączone. Wzór ten nie odzwierciedla jednak wszystkich próbek ECFC uzyskanych z ciąż z GDM, co wskazuje na heterogeniczność między próbkami. Próbki GDM2 i GDM3 wykazują większe tworzenie sieci w porównaniu z GDM1, choć wzorce łączności wydają się zmienione w stosunku do próbki UC. Co ważne, KAV mierzy kilka komponentów strukturalnych sieci ECFC, aby zidentyfikować zarówno oczywiste, jak i subtelne fenotypy.
Oprócz generowania obrazów szkieletowych i masek struktury sieci, KAV kwantyfikuje dziesięć parametrów struktury sieci, w tym liczbę poszczególnych struktur sieci, węzłów, węzłów trójdzielnych, węzłów czterodzielnych, odgałęzień, całkowitą długość odgałęzień, średnią długość odgałęzień, stosunek odgałęzień do węzłów, całkowitą liczbę sieci zamkniętych oraz powierzchnię sieci16. KAV kompiluje dane dla każdej próbki w jednej tabeli wartości dla wszystkich obrazów z przebiegu czasowego. Tabela 1 zawiera reprezentatywne surowe dane z jednego badania obrazowania dla jednej próbki ECFC. Parametry struktury sieci mierzone przez KAV są uporządkowane w kolumnach, a dane z kolejnych obrazów pozyskanych w czasie są uporządkowane w wierszach. Na przykład w naszych badaniach obraz 1 został uzyskany 30 min po wysianiu, a kolejne obrazy gromadzono co 15 min. Surowe wartości wygenerowane przez KAV mogą być następnie przedstawione na wykresach lub poddane dalszej analizie przy użyciu bardziej złożonych metod statystycznych16.
Przedstawiono pięć wykresów obrazujących średnie wartości struktury sieci z trzech oddzielnych eksperymentów dla pojedynczej próbki kontrolnej (UC) oraz trzech próbek GDM (Rysunek 4). Próbka UC w tych eksperymentach zachowywała się podobnie do wcześniej analizowanych próbek UC16. Wcześniej zidentyfikowano, że waskulogeneza ECFC in vitro jest dwufazowa i składa się z Fazy 1 (0 - 5 h) oraz Fazy 2 (5 - 10 h)16. Wzorzec ten jest spójny w obecnych badaniach, co potwierdza wykres przedstawiający dane dotyczące zamkniętych sieci (Rysunek 4). Próbka UC utworzyła większą liczbę zamkniętych sieci w porównaniu do wszystkich trzech próbek GDM. Co ciekawe, trzy z czterech próbek wykazują podobny czas osiągnięcia maksymalnej liczby zamkniętych sieci, co następuje między 2,5 a 3 godzinami. Jednak w próbce GDM2 osiągnięcie maksymalnej liczby zamkniętych sieci nastąpiło wolniej, tj. 5 godzin po wysiewie. Ponadto, mimo że próbka GDM2 utworzyła mniej sieci maksymalnych w porównaniu do próbki UC, utworzone sieci były utrzymywane w czasie podobnie jak w próbce UC. Przeciwnie, w przypadku GDM1 i GDM3, które utworzyły mniej sieci w porównaniu do próbki UC, odnotowano również ogólny spadek liczby sieci w czasie. Podsumowując, z wykresu przedstawiającego liczbę zamkniętych sieci wynika, że wszystkie próbki ECFC wykazały dwufazowy wzorzec tworzenia sieci, jednak tempo ich tworzenia oraz osiągnięta maksymalna liczba sieci różnią się w zależności od próbki.
Pole sieci reprezentuje średnią powierzchnię wewnątrz zamkniętych sieci utworzonych przez ECFCs. W związku z tym, im większa liczba zamkniętych sieci, tym mniejsze są obszary sieci. Wzorzec ten jest odzwierciedlony na wykresach pola sieci, gdzie próbka UC, która utworzyła większą liczbę zamkniętych sieci, wykazywała z czasem mniejsze obszary sieci w porównaniu do trzech próbek GDM (Rycina 4). Próbka GDM2, w której maksymalna liczba zamkniętych sieci została osiągnięta wolniej, wykazywała początkowo duże pole sieci, jednak z czasem obszar ten się ustabilizował, co było najbardziej zbliżone do próbki UC w 15. godzinie. Z biegiem czasu średnie pole sieci we wszystkich próbkach wzrosło z powodu destabilizacji sieci. Jednak niektóre próbki, takie jak GDM1 i GDM3, wykazały szybszy wzrost pola sieci, co prawdopodobnie wskazuje na mniejszą stabilność w porównaniu do pozostałych próbek, w których wzrost był bardziej stopniowy.
ECFCs eksponowane na GDM wykazują obniżoną stabilność sieci
Stosunek odgałęzień do węzłów jest nowym fenotypem obliczanym przez KAV i zidentyfikowanym w naszych poprzednich badaniach, który wskazuje na łączność sieci16. W niniejszym badaniu próbka UC wykazała niski stosunek, co reprezentowało wysoki poziom łączności sieci utrzymywany w czasie (Rysunek 4). Przeciwnie, trzy próbki GDM miały wyższy stosunek odgałęzień do węzłów, szczególnie w fazie 2 tworzenia sieci. Obserwacja ta świadczy o zmniejszonej łączności i destabilizacji struktur sieciowych.
Ogólnie GDM1 utworzyło mniej węzłów i odgałęzień w porównaniu z pozostałymi próbkami, a redukcja ta utrzymała się w trakcie trwania eksperymentu. GDM2 i GDM3 utworzyły i utrzymały większą liczbę odgałęzień w porównaniu z próbką UC, szczególnie między 5 a 15 h. Liczba węzłów wykrytych w sieciach GDM2 i GDM3 była bardziej zbliżona do liczby węzłów w próbce UC, zwłaszcza między 10 a 15 h. Większa liczba odgałęzień przy podobnej liczbie węzłów może wyjaśniać zwiększony stosunek liczby odgałęzień do węzłów, widoczny w późniejszych punktach czasowych w próbkach GDM. Co ważne, jednoczesne zmiany liczby odgałęzień i węzłów mogą być trudne do zinterpretowania na oddzielnych wykresach. Jednakże nowy stosunek liczby odgałęzień do węzłów oferuje sposób oceny tego, jak zmiany w liczbie zarówno odgałęzień, jak i węzłów – w tym niewielkie zmiany trudne do wykrycia na pojedynczych wykresach – wpływają na zmienioną łączność sieci.

Rysunek 1: Schemat przedstawiający 10 parametrów określanych ilościowo za pomocą kinetycznej analizy waskulogenezy (KAV). Metoda KAV określa ilościowo dziesięć różnych parametrów struktury sieci. Wszystkie parametry są oznaczone kolorami i ponumerowane na schemacie (1-10). Parametry obejmują liczbę odgałęzień (zielony, 1), sieci zamkniętych (niebieski, 2) i węzłów (czerwony, 3), średnią powierzchnię sieci (pomarańczowy, 4), liczbę struktur sieciowych (czarny, 5), węzły trójdzielne (żółty, 6) i węzły czterodzielne (fioletowy, 7), a także całkowitą (9) i średnią (10) długość odgałęzień oraz stosunek liczby odgałęzień do liczby węzłów (10). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 2: Analiza kinetyczna waskulogenezy (KAV) ilościowo określa strukturę sieci. (A) Przedstawienia szkieletowe i maskowe struktur sieci zidentyfikowanych przez KAV na obrazach w kontraście fazowym stanowią wizualizację struktur sieci wykrytych i zmierzonych za pomocą KAV. Pasek skali = 500 µm. (B) Reprezentatywny obraz w kontraście fazowym sieci ECFC 10 h po wysianiu, charakteryzujący się niskim kontrastem i widocznymi liniami siatki powstałymi w wyniku łączenia pojedynczych obrazów. W celu poprawy dokładności pomiarów, w przypadku niskiego kontrastu obrazów fazowych lub wystąpienia artefaktów siatki, wtyczka KAV umożliwia wybór różnych metod progowania, takich jak metoda średniej (Mean) lub metoda Otsu. Przedstawiono szkielety i maski obrazu w kontraście fazowym dla obu metod progowania (Mean oraz Otsu). Obrazy w kontraście fazowym wykonano przy użyciu obiektywu 10X. Pasek skali = 500 µm. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 3: Wewnątrzmaciczna ekspozycja na GDM zmienia tworzenie sieci ECFC. Obrazy tworzenia sieci ECFC rejestrowano w odstępach 15 min przez 15 h za pomocą mikroskopii kontrastowej. Reprezentatywne obrazy kontrastowe w odstępach 5 h, zaczynając od momentu wysiewu (0,5 h), przedstawiono dla UC oraz trzech próbek GDM. Obrazy kontrastowe wykonano przy użyciu obiektywu 10X. Pasek skali = 500 µm. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Rysunek 4: ECFC narażone na wewnątrzmaciczny GDM wykazują zmienioną kinetykę tworzenia sieci. ECFC pobrano z ciąży niepowikłanej (UC, kolor czarny) oraz z trzech ciąż powikłanych cukrzycą ciążową (GDM 1-3, kolor czerwony). Obrazy w kontraście fazowym rejestrowano w odstępach 15 min przez 15 h. Oprogramowanie Kinetic Analysis of Vasculogenesis (KAV) służyło do ilościowego oznaczania zamkniętych sieci, obszarów sieci, odgałęzień, węzłów oraz stosunku odgałęzień do węzłów. Przedstawione dane obrazują średnią ± błąd standardowy średniej (SEM) z trzech oddzielnych eksperymentów dla każdej próbki. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Wideo 1: Wewnątrzmaciczna ekspozycja na GDM zmienia kinetykę tworzenia sieci ECFC. Obrazy tworzenia sieci ECFC rejestrowano co 15 min przez 15 h za pomocą mikroskopii kontrastowej. Przedstawiono obrazy kontrastowe dla grupy UC oraz trzech próbek GDM w czasie 15 h, zaczynając od 0,5 h. Pasek skali reprezentuje 500 µm. Kliknij tutaj, aby obejrzeć to wideo. (Kliknij prawym przyciskiem myszy, aby pobrać.)
| Obraz | Całkowite sieci rozgałęzień | Węzły | Triplety | Czwórki | Gałęzie | Całkowita długość rozgałęzień* | Średnia długość gałęzi* | Stosunek liczby odgałęzień do liczby węzłów | Sieci zamknięte | Obszar sieciowy** |
| 1 | 382 | 290 | 278 | 11 | 943 | 47210 | 124 | 3.25 | 9 | 480214 |
| 2 | 284 | 345 | 337 | 8 | 940 | 50699 | 179 | 2.72 | 16 | 264469 |
| 3 | 205 | 376 | 366 | 9 | 910 | 55728 | 272 | 2.42 | 27 | 150621 |
| 4 | 162 | 422 | 407 | 15 | 947 | 59692 | 368 | 2.24 | 40 | 98713 |
| 5 | 132 | 454 | 441 | 12 | 967 | 61923 | 469 | 2.13 | 53 | 72951 |
| 6 | 122 | 435 | 419 | 14 | 907 | 62587 | 513 | 2.09 | 68 | 56948 |
| 7 | 88 | 429 | 411 | 17 | 852 | 63983 | 727 | 1.99 | 74 | 51429 |
| 8 | 68 | 451 | 437 | 13 | 874 | 63960 | 941 | 1.94 | 74 | 51780 |
| 9 | 93 | 437 | 417 | 18 | 905 | 60901 | 655 | 2.07 | 49 | 82919 |
| 10 | 126 | 511 | 487 | 24 | 1075 | 61981 | 492 | 2.1 | 37 | 116857 |
| 11 | 49 | 397 | 384 | 12 | 737 | 61823 | 1262 | 1.86 | 79 | 48805 |
| 12 | 81 | 376 | 364 | 10 | 751 | 56817 | 701 | 2 | 63 | 64431 |
| 13 | 98 | 509 | 482 | 26 | 1028 | 60799 | 620 | 2.02 | 44 | 98056 |
| 14 | 93 | 449 | 416 | 31 | 909 | 58282 | 627 | 2.02 | 45 | 94081 |
| *zaokrąglono do najbliższego mikrona, **zaokrąglono do najbliższego mikrona^2 |
Tabela 1: Reprezentatywne surowe dane z jednego badania obrazowania dla jednej próbki ECFC.