W niniejszym rozdziale opisujemy praktyczne zastosowanie procedury kwantyfikacji obrazów i genotypowania embrionów, opublikowanej w innym opracowaniu3. Schemat tej metody przedstawiono na Ryc. 1. Aby zilustrować sposób zastosowania tej metody, przeprowadzono hybrydyzację in situ (ISH) dla dnmt3bb.1 w embrionach o wieku 3 hpf z krzyżówki wewnątrzrasowej runx1W84X/+ 22 (Ryc. 2). Wykonano obrazowanie 130 embrionów, stosując te same warunki oświetlenia, które szczegółowo opisano w protokole, i nadano im unikalne numery. Po obrazowaniu każdy embrion przeniesiono do probówki PCR w celu genotypowania. Następnie przeprowadzono analizę obrazów, aby przypisać wartość intensywności pikseli do każdego zdjęcia. Genotyp przypisano do odpowiadającego mu obrazu, a wartości intensywności pikseli pogrupowano zgodnie z genotypem w celu przeprowadzenia analizy statystycznej. Wykryto spadek ekspresji dnmt3bb.1 u mutantów runx1W84X/W84X (Ryc. 2A,B)3, co jest zgodne z wcześniejszymi obserwacjami5. Co ciekawe, embriony heterozygnotyczne runx1W84X/+ nie wykazały istotnych różnic w ekspresji dnmt3bb.1 (Ryc. 2A,B) w porównaniu z rodzeństwem o genotypie dzikim, co sugeruje, że jedna kopia Runx1 jest wystarczająca do utrzymania ekspresji dnmt3bb.1 na odpowiednim poziomie.
Wiele mutantów rybek danio nie wykazuje fenotypu embrionalnego, który można wykryć za pomocą innych technologii utraty funkcji, takich jak oligonukleotydy morfolino (MO). Rozbieżność ta może wynikać z szeregu przyczyn, w tym efektów poza docelowych (off-target), kompensacji białek matczynych, obecności allelu hypomorficznego23 lub niedawno odkrytego zjawiska kompensacji genetycznej24,25,26,27. W tym przykładzie badaliśmy, czy ekspresja runx1 była zmniejszona lub utracona u mutantów lmo4auob10, ponieważ wcześniej opublikowane dane z użyciem MO przeciwko lmo4a sugerowały, że poziom runx1 jest obniżony u morfantów lmo4a28. Tutaj analiza nie wykazała istotnych różnic w ekspresji runx1 pomiędzy typem dzikim a homozygotycznymi mutantami lmo4auob103 (Rysunek 3A,B). Dalsza analiza za pomocą qPCR pojedynczych embrionów wykazała niewielki, ale istotny spadek ekspresji runx1 u mutantów lmo4auob10 (Rysunek 3C). Zatem możliwe jest, że kwantyfikacja obrazowa może nie być w stanie wykryć niewielkich różnic w poziomach ekspresji. Alternatywnie, brak różnic między genotypami, który stwierdziliśmy, jest rzeczywisty, a eksperymenty qPCR wykrywają zmiany w ekspresji runx1 w innych tkankach, takich jak telencefalon, gdzie ekspresowane są zarówno lmo4a, jak i runx1. Badacze powinni zawsze weryfikować swoje wyniki niezależną metodą, taką jak qPCR, najlepiej poprzedzając ją np. wzbogaceniem frakcji interesującej tkanki za pomocą cytometrii przepływowej.
W rzadkich przypadkach, gdy ISH wykazuje wysokie tło (Rycyna 3D), wartość intensywności pikseli w tym obszarze jest tak wysoka, że odjęcie jej od wartości sygnału daje liczbę ujemną; w takich sytuacjach te embriony zostają wykluczone z analizy. W naszym doświadczeniu zdarzyło się to w około 0,4% embrionów badanych sondą runx13, jednak odsetek ten może różnić się w zależności od eksperymentów, sond lub partii odczynników. Choć może to stanowić ograniczenie metody, niska częstotliwość występowania wysokiego tła jest bardzo mało prawdopodobna, aby wpłynęła na ogólne wyniki.
Aby przetestować wpływ wyboru różnych obszarów do korekcji tła, zmierzyliśmy najpierw intensywność pikseli sygnału ISH dla runx1 w embrionach w 28 hpf, wykorzystując różne regiony do korekcji tła (Rycyna 4). Wybrano cztery różne regiony: dwa w obszarze tułowia (R1 i R2), jeden w obszarze żółtka (niebarwiony, ale prawdopodobnie akumulujący barwienie tła) oraz mniejszy obszar przed ROI (R4, Rycyna 4B). Pomiar intensywności pikseli w tych regionach wykazał relatywnie stabilną różnicę w intensywności między ROI a dowolnym obszarem tła (Rycyna 4C). Jednak region R3 zawsze wykazywał bardzo wysokie wartości (powyżej tych w ROI). Po inwersji i konwersji do formatu 8-bitowego obszar żółtka wydaje się bardzo jasny, a zatem nie nadaje się do użytku w korekcji tła. R2 znajdował się bliżej ROI, ale zawierał częściowo sygnał ISH, a jego użycie do korekcji obniżyło średnią intensywność pikseli w porównaniu z regionem R1 (znajdującym się dalej grzbietowo, z dala od sygnału ISH) lub R4. Zatem zarówno R1, jak i R4 są odpowiednimi obszarami, które można wykorzystać do korekcji tła (mimo że powierzchnia R4 jest mniejsza niż R1). Następnie chcieliśmy porównać, w jaki sposób użycie R1 lub R4 wpływa na wyniki przy porównywaniu ekspresji runx1. W tym celu przeprowadziliśmy krzyżowanie heterozygotów dll4+/-29 i przeanalizowaliśmy ekspresję runx1 w losowo wybranych embrionach typu dzikiego oraz dll4-/- (Rycyna 4E). Chociaż użycie R1 lub R4 do korekcji tła wpłynęło na wartości indywidualne, średnie intensywności pikseli w obrębie tego samego genotypu nie różniły się istotnie (Rycyna 4E). Co więcej, porównanie ekspresji runx1 nadal daje podobne średnie wartości intensywności między genotypami przy użyciu obszarów R1 lub R4 jako korekcji tła (odpowiednio µR1=16.3 i µR4=18.2). Podsumowując, stwierdziliśmy, że choć wybór obszaru tła jest ważny, głównym kryterium jest to, aby nie obejmował on obszarów żółtka (podatnych na akumulację barwienia tła) oraz aby nie zawierał żadnego (specyficznego) barwienia, które mogłoby zafałszować wartości intensywności pikseli tła.

Rysunek 1: Schemat procedury równoległego oznaczania ilościowego obrazów i genotypowania. Zarodki pobrane z inbreedingu ryb heterozygotycznych pod względem allelu mutacyjnego poddaje się hybrydyzacji in situ (ISH) w celu wykrycia mierzonego genu zgodnie ze standardowym protokołem ISH. Po obrazowaniu izoluje się DNA genomowe zgodnie z protokołem HotSHOT, dodając bufor lizujący bezpośrednio do zarodka w probówce PCR o pojemności 0,2 ml, a następnie inkubując przez 30 min w temperaturze 95 °CWyizolowane DNA jest wykorzystywane do genotypowania embrionów za pomocą metody PCR, PCR i polimorfizmu długości fragmentów restrykcyjnych (RFLP), testów KASP lub każdej innej odpowiedniej metody. Równolegle obrazy każdego embrionu są odwracane i konwertowane do 8-bitowej skali szarości. Ręcznie wybierane i mierzone są obszary zainteresowania (ROI) o identycznym kształcie i rozmiarze, zawierające sygnał ISH (żółty) oraz tło (niebieskie). Pomiary przypisane do odpowiednich genotypów są poddawane analizie statystycznej. Rysunek zaadaptowany z Dobrzycki et al.3 Kliknij tutaj, aby wyświetlić powiększoną wersję tej figury.

Rysunek 2: Ilościowa analiza obrazów w runx1 mutanty wykazują obniżone poziomy dnmt3bb.1 ekspresja metodą ISH. (A) Przykładowe obrazy ISH w embrionach typu dzikiego w 3 hpf (niebieski), runx1+/W84X (zielony) i runx1W84X/W84X zarodki (pomarańczowe), wykazujące dnmt3bb.1 ekspresja w aorcie grzbietowej. (B) Wartości intensywności pikseli dnmt3bb.1 mRNA w runx1W84X/W84Xliczba zarodków (n=36) jest istotnie mniejsza w porównaniu do typów dzikich (n=32) oraz heterozygot (n=62) (ANOVA, p < 0,01). Współczynnikiy zmienności wynoszą odpowiednio 24%, 2% i 21% dla grup typu dzikiego, heterozygnot i mutantów. Niebieskie, zielone i pomarańczowe punkty danych odpowiadają przykładowym obrazom z panelu A. Słupki reprezentują średnią ± odchylenie standardowe (s.d.). ***p<0,01 (test post-hoc Gamesa-Howella). Rysunek zaadaptowano z Dobrzycki i wsp.3 Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Pomiar runx1 poziomy ekspresji metodą ISH w lmo4auob100mutanty. (A) Reprezentatywne obrazy ISH dla runx1 w 28 hpf typu dzikiego (niebieski), heterozygoty (zielony) oraz lmo4auob10/uob10 (pomarańczowe) embriony, wykazujące ekspresję w aorcie grzbietowej. (B) Ilościowe oznaczenie runx1 sygnał mRNA, wykryty za pomocą ISH, z embrionów typu dzikiego (n=15) i heterozygot w 28 hpf lmo4a+/- (het) (n=34) i lmo4auob10/uob10 zarodki mutanty (n=18) z jednego cyklu nie wykazują istotnych różnic w runx1 intensywność pikseli pomiędzy różnymi genotypami (ANOVA,> p 0,6). Niebieskie, zielone i pomarańczowe punkty danych odpowiadają przykładowym obrazom z panelu A. Słupki reprezentują średnią ± odchylenie standardowe (s.d.C) Wykresy pudełkowe przedstawiające znormalizowane runx1 poziomy mRNA (2-ΔCt) u pojedynczych osobników typu dzikiego (kolor niebieski; n=12) oraz lmo4auob100/uob10 zarodki (mut, pomarańczowe; n=12), w których pomiarom qRT-PCR towarzyszył spadek poziomu runx1 u mutantów w porównaniu z typem dzikim.p < 0.05 (t test). (D) Przykład eksperymentu ISH na embrionie w 28 hpf (barwionym na runx1, żółte groty strzałek) wykazujące wysokie tło. Rysunek zaadaptowano z pracy Dobrzycki i wsp.3 Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4: Wpływ korekty natężenia tła na wyniki pomiarów. (A) Reprezentatywny obraz runx1 Barwienie ISH w embrionie typu dzikiego w 28. godzinie po zapłodnieniu (hpf).B) Ten sam obraz po inwersji i konwersji do formatu 8-bitowego. Obszar zainteresowania (ROI) zaznaczono kolorem zielonym, a cztery różne obszary wykorzystane do korekcji tła (R1-R4) zaznaczono kolorem żółtym. (C) Pomiary surowej intensywności pikseli we wszystkich obszarach przedstawionych na panelu B. Zwróć uwagę, że intensywność w R3 (żółtko) jest konsekwentnie wyższa niż rzeczywisty sygnał ISH w ROI (n=1). (D) Runx1 poziomy ekspresji w ROI z wykorzystaniem obszarów tła R1, R2 i R4. Powierzchnie dla ROI, R1, R2 i R3 wyniosły ok. 28 50 pikseli; R4 ok. 850 pikseli. Należy zauważyć, że tło R3 nie zostało użyte do tego porównania, ponieważ korekcja tła (ROI-R3) konsekwentnie dawała wartości ujemne. (E) Runx1 poziomy ekspresji w typie dzikim i dll4-/- mutanty, wykorzystując R1 lub R4 do korekcji tła (n=10 dla każdej próbki). Analiza statystyczna w panelach D i E została przeprowadzona przy użyciu nieparametrycznego testu Kruskala-Wallisa, przy założeniu, że wartości intensywności pikseli nie mają rozkładu normalnego. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.