$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Sekwencje DNA lub RNA zawierają informacje genetyczne dla podstawowych fenotypów, które mogą być analizowane metodami fizjologicznymi i biochemicznymi lub obserwowane za pomocą dowodów morfologicznych i kopalnych. W pewnym sensie informacja genetyczna jest bardziej wiarygodna niż ocena fenotypów zewnętrznych, ponieważ ta pierwsza jest podstawą dla drugiej. W badaniach ewolucyjnych dowody kopalne są bardzo bezpośrednie i przekonujące. Jednak wiele organizmów, takich jak mikroorganizmy, ma niewielkie szanse na uformowanie skamieniałości podczas długich epok geologicznych. Dlatego informacje molekularne, takie jak sekwencje nukleotydowe i sekwencje aminokwasów z pokrewnych istniejących organizmów, są cenne dla badania relacji ewolucyjnych1. W niniejszym badaniu zapewniono proste wprowadzenie podstawowej wiedzy filogenetycznej i łatwy do nauczenia się protokół dla nowicjuszy, którzy muszą samodzielnie skonstruować drzewo filogenetyczne
.
Zarówno sekwencje DNA (nukleotydowe), jak i białkowe (aminokwasowe) mogą być używane do wnioskowania o filogenetycznych związkach między homologicznymi genami, organellami, a nawet organizmami2. Bardziej prawdopodobne jest, że sekwencje DNA ulegną zmianom zachodzącym podczas ewolucji. W przeciwieństwie do tego, sekwencje aminokwasów są znacznie bardziej stabilne, biorąc pod uwagę, że synonimiczne mutacje w sekwencjach nukleotydowych nie powodują mutacji w sekwencjach aminokwasów. W rezultacie, sekwencje DNA są przydatne do porównywania homologicznych genów z blisko spokrewnionych organizmów, podczas gdy sekwencje aminokwasów są odpowiednie dla homologicznych genów z odległych organizmów3.
Analiza filogenetyczna rozpoczyna się od dopasowania sekwencji aminokwasów lub nukleotydów4 pobranych z bazy danych sekwencjonowania genomu z adnotacjami5 wymienionych w formacie FASTA, tj. przypuszczalnych lub wyrażonych sekwencji białek, sekwencji RNA lub sekwencji DNA. Warto zauważyć, że bardzo ważne jest zebranie wysokiej jakości sekwencji do analizy, a do analizy powiązań filogenetycznych można wykorzystać tylko sekwencje homologiczne. Do wyrównywania sekwencji można wykorzystać wiele różnych platform, takich jak Clustal W, Clustal X, Muscle, T-coffee, MAFFT. Najczęściej używanym jest Clustal Omega6,7 (http://www.ebi.ac.uk/Tools/msa/clustalo/), z którego można korzystać online lub można go pobrać bezpłatnie. Narzędzie do wyrównywania ma wiele parametrów, które użytkownik może dostosować przed rozpoczęciem wyrównania, ale parametry domyślne działają dobrze w większości przypadków. Po zakończeniu procesu wyrównane sekwencje powinny zostać zapisane w odpowiednim formacie dla następnego kroku. Powinny być one następnie edytowane lub przycinane za pomocą oprogramowania do edycji, takiego jak BioEdit, ponieważ konstrukcja drzewa filogenetycznego przez MEGA wymaga, aby sekwencje były równej długości (w tym zarówno skróty aminokwasów, jak i łączniki). W wyrównanej sekwencji dowolna pozycja bez aminokwasu lub nukleotydu jest reprezentowana przez łącznik "-"). Ogólnie rzecz biorąc, wszystkie wystające aminokwasy lub nukleotydy na obu końcach wyrównania powinny zostać usunięte. Ponadto kolumny zawierające źle wyrównane sekwencje w wyrównaniu mogą zostać usunięte, ponieważ przekazują niewiele cennych informacji, a czasami mogą dawać mylące lub fałszywe informacje3. Kolumny zawierające jeden lub więcej łączników można usunąć w tym momencie lub na późniejszym etapie budowy drzewa. Alternatywnie można je wykorzystać do obliczeń filogenetycznych. Po zakończeniu wyrównywania i przycinania sekwencji wyrównane sekwencje należy zapisać w formacie FASTA lub żądanym formacie do późniejszego wykorzystania.
Wiele platform oprogramowania udostępnia funkcje budowy drzew przy użyciu różnych metod lub algorytmów. Ogólnie rzecz biorąc, metody te można sklasyfikować jako metody macierzy odległości lub metody danych dyskretnych. Metody macierzy odległości są proste i szybkie do obliczenia, podczas gdy metody danych dyskretnych są skomplikowane i czasochłonne. W przypadku bardzo blisko spokrewnionych taksonów o wysokim stopniu współdzielenia identyczności sekwencji aminokwasów lub nukleotydów stosuje się metodę macierzy odległości (Neighbor Joining: NJ; Metoda grup par nieważonych ze średnią arytmetyczną: UPGMA) jest odpowiednia; w przypadku daleko spokrewnionych taksonów – metodę danych dyskretnych (maksymalne prawdopodobieństwo: ML; Maksymalna oszczędność: MP; Wnioskowanie bayesowskie) jest optymalne3,8. W tym badaniu metody uczenia maszynowego w MEGA (6.0.6) i wnioskowaniu bayesowskim (MrBayes 3.2) zostały zastosowane do konstruowania drzew filogenetycznych9. Idealnie, gdy stosuje się odpowiedni model i parametry, wyniki uzyskane za pomocą różnych metod mogą być spójne, a tym samym bardziej wiarygodne i przekonujące.
Dla drzewa filogenetycznego ML zbudowanego przy użyciu MEGA10, plik sekwencji wyrównanej w formacie FASTA musi zostać przesłany do programu. Pierwszym krokiem jest więc wybór optymalnego modelu podstawiania dla wgrywanych danych. Wszystkie dostępne modele substytucji są porównywane na podstawie przesłanych sekwencji, a ich końcowe wyniki zostaną pokazane w tabeli wyników. Wybierz model z najmniejszym wynikiem Bayesowskiego kryterium informacyjnego (BIC) (wymienionym jako pierwszy w tabeli), ustaw parametry uczenia maszynowego zgodnie z zalecanym modelem i rozpocznij obliczenia. Czas obliczeń waha się od kilku minut do kilku dni, w zależności od złożoności załadowanych danych (długość sekwencji i liczba taksonów) oraz wydajności komputera, na którym uruchamiane są programy. Po zakończeniu obliczeń drzewo filogenetyczne zostanie wyświetlone w nowym oknie. Zapisz plik jako "FileName.mat". Po ustawieniu parametrów określających wygląd drzewa, zapisz jeszcze raz. Korzystając z tej metody, MEGA może generować figury drzewa filogenetycznego klasy publikacji.
Do budowy drzewa za pomocą MrBayes11, pierwszym krokiem jest przekształcenie wyrównanej sekwencji, która zwykle jest wyświetlana w formacie FASTA, w format nexus (.nex jako typ pliku). Przekształcanie plików FASTA do formatu nexus można przetwarzać w MEGA. Następnie wyrównaną sekwencję w formacie nexus można przesłać do MrBayes. Po pomyślnym przekazaniu pliku określ szczegółowe parametry obliczeń drzewa. Parametry te obejmują szczegóły, takie jak model podstawiania aminokwasów, współczynniki zmienności, numer łańcucha dla sprzężenia Monte Carlo łańcucha Markowa (MCMC), liczba ngen, średnie odchylenie standardowe częstości podziału i tak dalej. Po określeniu tych parametrów rozpocznij obliczenia. Na koniec na ekranie zostaną wyświetlone dwie figury drzewa w kodzie ASC II, jedna pokazująca wiarygodność kladu, a druga pokazująca długość gałęzi.
Wynik drzewa zostanie automatycznie zapisany jako "FileName.nex.con". Ten plik drzewa może być otwierany i edytowany przez FigTree, a rysunek wyświetlany w FigTree może być dalej modyfikowany, aby był bardziej odpowiedni do publikacji.
W tym badaniu, jako przykład, analizowano 228 białek SWEET, w tym 35 SWEETs od jednokomórkowych eukariontów i 57 SemiSWEETs od prokariontów. Zarówno SWEETs, jak i SemiSWEET zostały scharakteryzowane jako transportery glukozy, fruktozy lub sacharozy przez błony12,13. Analiza filogenetyczna sugeruje, że dwie domeny MtN3/śliny zawierające SWEETs mogą pochodzić z ewolucyjnej fuzji bakterii SemiSWEET i archeon14.