$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Identyfikacja kompletnego zestawu elementów kodujących w genomie była głównym celem od momentu rozpoczęcia Projektu Poznania Ludzkiego Genomu i pozostaje głównym celem w kierunku zrozumienia systemów biologicznych i etiologii chorób genetycznych1,2,3,4. Postępy w technikach NGS doprowadziły do wytworzenia całych sekwencji genomu dla dużej liczby organizmów, w tym kręgowców, bezkręgowców, drożdży i roślin5. Ponadto wysokoprzepustowe metody sekwencjonowania transkrypcji jeszcze bardziej ujawniły złożoność transkryptomu komórkowego i zidentyfikowały tysiące nowych cząsteczek RNA o funkcjach zarówno kodowania białek, jak i niekodowania6,7. Dekodowanie tej ogromnej ilości informacji o sekwencji jest procesem ciągłym, a wyzwania związane z kompleksowymi wysiłkami w zakresie adnotacji genów nadal istnieją8.
Niedawny rozwój metod profilowania translacyjnego, w tym profilowania rybosomów9,10 i sekwencjonowanie polirybosomów11, dostarczyły dowodów wskazujących, że setki niekanonicznych zdarzeń translacji mapują na obecnie nieoznaczone sORFy w całym genomie, z potencjałem do generowania małych białek zwanych mikroproteinami lub mikropeptydami12,13,14,15,16,17. Mikroproteiny pojawiły się jako nowa klasa wszechstronnych białek, wcześniej pomijanych przez standardowe metody adnotacji genów ze względu na ich mały rozmiar (<100 aminokwasów) i brak klasycznych cech genów kodujących białka8,12,18,19,20. Mikroproteiny zostały opisane praktycznie we wszystkich organizmach, w tym drożdżach21,22, flies17,23,24, and mammals25,26,27,28i wykazano, że odgrywają kluczową rolę w różnych procesach, w tym w rozwoju, metabolizmie i sygnalizacji stresu19,20,29,30,31,32,33,34. W związku z tym konieczne jest kontynuowanie eksploracji genomu w poszukiwaniu dodatkowych członków tej długo pomijanej klasy funkcjonalnych małych białek.
Pomimo powszechnego uznania biologicznego znaczenia mikrobiałek, ta klasa genów pozostaje zdecydowanie niedostatecznie reprezentowana w adnotacjach genomu, a ich dokładna identyfikacja nadal stanowi ciągłe wyzwanie, które hamuje postęp w tej dziedzinie. Ostatnio opracowano różne narzędzia obliczeniowe i metody eksperymentalne w celu przezwyciężenia trudności związanych z identyfikacją sekwencji kodujących mikrobiałka (omówione obszernie w kilku obszernych przeglądach8,35,36,37). Wiele ostatnich badań nad identyfikacją mikrobiałek38,39,40,41,42,43,44,45,46,47 w dużym stopniu polegały na użyciu jednego z takich algorytm o nazwie PhyloCSF48,49, potężne podejście genomiki porównawczej, które można wykorzystać do odróżnienia konserwatywnych regionów genomu kodujących białka od tych, które nie kodują.
PhyloCSF porównuje częstość substytucji kodonów (CSF) za pomocą wielogatunkowych dopasowań nukleotydów i modeli filogenetycznych do wykrywania ewolucyjnych sygnatur genów kodujących białka. To empiryczne podejście oparte na modelu opiera się na założeniu, że białka są przede wszystkim konserwatywne na poziomie aminokwasów, a nie sekwencji nukleotydów. Dlatego substytucje synonimiczne kodonów, które kodują ten sam aminokwas, lub substytucje kodonów do aminokwasów o konserwatywnych właściwościach (tj. ładunek, hydrofobowość, polarność) są oceniane pozytywnie, podczas gdy substytucje niesynonimiczne, w tym podstawienia missense i nonsensowne, są oceniane negatywnie. PhyloCSF jest szkolony na danych z całego genomu i okazał się skuteczny w ocenianiu krótkich fragmentów sekwencji kodującej (CDS) w izolacji od pełnej sekwencji, co jest niezbędne podczas analizy mikrobiałek lub pojedynczych eksonów standardowych genów kodujących białka48,49.
Warto zauważyć, że niedawna integracja węzłów śledzenia PhyloCSF w przeglądarce genomów Uniwersytetu Kalifornijskiego w Santa Cruz (UCSC)49,50,51 umożliwia badaczom ze wszystkich środowisk łatwy dostęp do przyjaznego dla użytkownika interfejsu w celu wyszukiwania informacji o potencjał kodowania białek w regionach genomu. Protokół opisany poniżej zawiera szczegółowe instrukcje dotyczące ładowania koncentratorów śledzenia PhyloCSF w przeglądarce genomu UCSC, a następnie przesłuchiwania interesujących regionów genomowych w celu zbadania regionów kodujących białka o wysokim poziomie ufności (lub ich braku). Dodatkowo, w przypadku zaobserwowania dodatniego wyniku PhyloCSF, wyznaczane są kroki w celu dalszej analizy potencjału kodowania mikrobiałek i wydajnego generowania wielogatunkowych dopasowań zidentyfikowanych sekwencji aminokwasów w celu zilustrowania zachowania sekwencji międzygatunkowych. Wreszcie, w dyskusji przedstawiono kilka dodatkowych publicznie dostępnych zasobów i narzędzi w celu zbadania zidentyfikowanych cech mikrobiałek, w tym przewidywanych struktur domenowych i wglądu w przypuszczalną funkcję mikroprotein.