Dopasowanie sekwencji

Search…
​← Spis treści​
Dopasowanie sekwencji
Kolejnym etapem budowania drzew filogenetycznych jest dopasowanie sekwencji. Jest to bardzo ważny proces, złe dopasowanie sekwencji skutkuje utworzeniem nieprawidłowych drzew.
Etapy i programy
Po zebraniu sekwencji w pliku FASTA przechodzimy do ich dopasowania. Dopasowanie można wykonywać przy pomocy programów uruchamianych z linii komend. Kolejnym etapem jest ręczne poprawianie dopasowania zwykle połączone z przycinaniem zebranych funkcji. Ten etap najwygodniej jest wykonywać używając programów ,,okienkowych''. Programy z interfejsem graficznym zwykle mają także funkcje dopasowania sekwencji, często zresztą poprzez uruchamianie programów działających w linii komend.
Długą listę takich programów można znaleźć np. na Wikipedii​
Programy dopasowujące wiele sekwencji uruchamiane z linii komend
Programów uruchamianych z linii komend jest wiele, zwykle każdy z nich używa osobnego algorytmu dopasowania lub jego modyfikacje. Należą do nich:
​Clustal W i Clustal Omega​
​MUSCLE​
​MAFFT​
​T-Coffee​
​PRANK​
Można je zainstalować na Debianie lub innych debianopodobnych dystybucjach komendą (usuń nazwy programów, których nie chcesz instalować):
1
sudo apt-get install clustalw clustalo muscle mafft t-coffee prank
Copied!
W internecie można także łatwo znaleźć strony, które umożliwiają procesu dopasowania online z wykorzystaniem tych (i innych) programów.
Programy z interfejsem graficznym do edycji dopasowań
W pracy z zestawami sekwencji praktyczne są programy z interfejsem graficznym, które posiadają szereg użytecznych funkcji ułatwiających ich porządkowanie i edycję. a także np. uruchamianie programów dopasowujących wstępnie sekwencje czy obliczanie i rysowanie drzew.
Praktyka: Proste dopasowanie sekwencji atp6
Pobierz z GenBank-u następujące sekwencje i zapisz je w pliku atp6-surowe.fasta:
1
AY007817 Daucus carota
2
AY847285 Brassica juncea
3
AY961627 Oryza sativa
4
EU882268 Sapria himalayana
5
FJ595983 Helianthus annuus
6
GU075810 Zea mays
7
HQ593780 Ajuga reptans
8
HQ593782 Mimulus guttatus
9
KC825300 Fragaria virginiana
10
KC879635 Magnolia stellata
11
KU180476 Centaurea scabiosa
12
KX524674 Lindenbergia siniaca
Copied!
Zmień domyślne nazwy ,,GenBankowe'' na nazwy wg. schematu: nr_rodzaj_gatunek, np: AY007817_Daucus_carota. Wykonaj poprawki w edytorze tekstu, lub używając odpowiedniej komendy, np:
1
  sed -i 's/>\([[:alnum:]]\+\)\.1 \([[:alnum:]]\+\) \([[:alnum:]]\+\).\+/>\1_\2_\3/' atp6-surowe.fasta
2
  
Copied!
Zapisz jako atp6-surowe.fasta.
Zanim przejdziemy do dopasowania sekwencji w linii komend, użyjemy programu AliView, który można pobrać ze strony domowej programu: http://ormbunkar.se/aliview/. AliView pozwala na wygodne przeglądanie i edycję zestawu sekwencji a także posiada szereg innych przydatnych funkcji do których z pewnością należy możliwość uruchamiania programów dopasowujących sekwencje. Warto zajrzeć na stronę programu http://ormbunkar.se/aliview/ gdzie przedstawione są jego możliwości.
Uruchom AliView i otwóż plik z sekwencjami.
AliView - interfejs
Aliview dopasowuje sekwencje używając zewnętrznych programów.
Można zmienić uruchamiany program wybierając z menu
AliView - zmiana programu dopasowującego sekwencje - menu
Pokazuje się okno:
AliView - zmiana programu dopasowującego sekwencje
Tu można wybrać program z domyślnie udostępnionych przez program, dodać inne i dostosować opcje ich uruchamiania. Instrukcje można znaleźć na stronie programu. Teraz zmień program na Mafft (ten wyżej).
Wybierz z paska narzędziowego Align->Realign everything
AliView - Dopasuj wszystko
Program zapyta czy na pewno dopasować sekwencje, zatwierdzamy. Pokazuje się okno z postępem pracy programu i w końcu widać wynik. Przesuń suwakiem na dole ekranu zestawienie sekwencji w prawo. Przyjrzyj się jak wyglądają dopasowane sekwencje.
Jak widać mają one różną długość, w końcu jednak znajduje się strefa dopasowania, która obejmuje odcinki wszystkich sekwencji, od 563 miejsca dopasowania. Przy okazji przyjrzyj się informacjom widocznym w pasku na dole okna. Są one dość jasne:
AliView - dopasowanie
Teraz kliknij na czwartą ikonę od lewej. Zaznaczone zostają kolorami kodony odpowiadające różnym aminokwasom.
AliView - ,,translacja''
Pojawiają się także nowe ikony pozwalające na zmianę m. in. sposobu wyświetlania a także ustawienia początkowego nukleotydu dla odczytu kodonów. Domyślnie ,,translacja'' rozpoczyna się od pierwszej pozycji zestawienia. Prawidłową pozycję startu odczytywania kodonu można odczytać np. z GenBanku. Sprawdź dane dla sekwencji KU180476. Sprawdź gdzie zaczyna się prawidłowe miejsce początku pierwszego kompletnego kodonu w części wspólnej dopasowanych sekwencji.
AliView - ,,translacja''
Teraz wróć do domyślnego widoku sekwencji i zaznacz wszystkie kolumny na lewo od miejsca, w którym rozpoczyna się wspólna część sekwencji:
AliView - interfejs
Następnie usuń je używając skrótu klawiszowego <Ctrl>+<Del> (warto go zapamiętać) lub odpowiedniej pozycji w pasku menu: Edit->Delete selected.
Podobnie przytnij prawą część dopasowania.
Teraz znów możesz sprawdzić ,,translację'' sekwencji, przy okazji zauważ, że (co nie jest zaskoczeniem) przy złym ustawieniu miejsca początku kodonu pojawia się wiele miejsc STOP oznaczonych szarym kolorem.
AliView - zły początek kodonu
Jeśli ustawimy podgląd za pomocą jednoliterowych oznaczeń aminokwasów, miejsca STOP będą wyświetlone jako litera x.
AliView - oznaczenia jednoliterowe
Teraz zastanów się, jak można znaleźć prawidłowe miejsce początku kodonów w działających genach, jeśli go nie znamy.
Wróć teraz do standardowego wyświetlania sekwencji i znajdź miejsce, w którym znajdują się indele w dwu sekwencjach. Zauważ, że pierwszym nukleotydem flankującym indel w sekwencji Oryza sativa jest G, który w zasadzie pasowałby równie dobrze po obu stronach indelu. Nieco inaczej jest w Ajuga reptans, tam nukleotyd C nie pasuje do żadnego nukleotydu w innych sekwencjach:
AliView - indele
W ich prawidłowym umieszczeniu może pomóc podświetlenie kodonów:
AliView - podświetlenie kodonów
Pamiętając, że indele w działających genach kodujących białka powinny obejmować trójki nukleotydów lub ich wielokrotność (inaczej zmienia się ramka odczytu), widać że G w pierwszej sekwencji zostało umieszczone prawidłowo, natomiast C przy drugim indelu należy przesunąć w lewo. Nukleotydy przesuwamy klikając na nie i przeciągając kursorem myszy, przy więcej niż jednym nukleotydzie należy wcześniej zaznaczyć odcinek, który chcemy przemieścić.
Teraz kodony są kompletne:
AliView - indele dopasowane
Zapisz plik jako atp6-dopasowane.fasta.
Przy pracy z AliView warto korzystać z licznych dostępnych opcji ułatwiających pracę. Przejrzyj te dostępne w menu, zazwyczaj łatwo domyślić się jakie mają zastosowanie.
Dopasowanie sekwencji z linii komend
Teraz przyszedł czas na dopasowanie sekwencji z użyciem linii komend. Najpierw użyjemy w tym celu programu muscle. Sprawdź jak go używać używając komendy:
1
muscle -h 
Copied!
lub
1
man muscle
Copied!
W najprostszej formie podajemy plik wejściowy i nazwę pliku wyjściowego:
1
muscle -in atp6-surowe.fasta -out atp6-wyrownanie-muscle.fasta
Copied!
Sprawdź plik wynikowy w Aliview.
Program mafft oferuje dużo więcej opcji niż muscle, uruchom:
1
man mafft
Copied!
Jak widać można w tym programie ustawiać np. wartości kar za otwieranie czy wydłużanie indeli, co może być przydatne przy sekwencjach, które mają ich wiele, zaraz się z takimi spotkamy, na razie uruchom mafft bez dodatkowych opcji:
1
mafft atp6-surowe.fasta > atp6-wyrownanie-mafft.fasta
Copied!
Po otwarciu pliku wynikowego w Aliview widać, mafft domyślnie zmienia wielkość liter w sekwencjach na małe, raczej tego nie chcemy, dlatego użyjemy opcji --preservecase
1
mafft --preservecase atp6-surowe.fasta > atp6-wyrownanie-mafft.fasta
Copied!
Teraz pobierz sekwencje o poniższych numerach i zapisz je w pliku trnL-trnF-GB.fasta:
1
JX073793
2
KU238878
3
GU818008
4
KX524675 
5
AY743457
6
AY575533
7
KU238873
8
KU238869
9
KU238874
10
KU238865
11
KU238864
12
KU238875
13
KU238866
14
KU238871
15
KU238876
16
KU238867
17
KU238877
18
KU238868
19
KU238872
20
KU238870
21
GU323356
Copied!
Sugeruję zrobić to za pomocą skryptu. W takim wypadku sprawdź czy pierwsza linia pliku FASTA nie jest pusta, jeśli tak to usuń ją.
Zmień nazwy sekwencji wg. schematu dla sekwencji atp6.
Teraz spróbuj dopasować sekwencje za pomocą programu muscle.
Następnie spróbuj mafft, także z opcją --auto oraz --op (sprawdź co oznaczają).
Porównaj wyniki i spróbuj poprawić wybrane dopasowanie ręcznie.
​← Spis treści​