Artykuł metodologiczny

Praktyczny przewodnik po filogenetyce dla laików

DOI:

10.3791/50975

5 lutego 2014

W tym artykule

Podsumowanie

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Tutaj opisujemy krok po kroku proces generowania wiarygodnych filogenez z zestawów danych sekwencji nukleotydów lub aminokwasów. Niniejszy przewodnik ma na celu służyć naukowcom lub studentom, którzy nie mają doświadczenia w analizie filogenetycznej.

Streszczenie

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Wielu badaczy, z niezwykle różnych ognisk, stosuje filogenetykę do swoich pytań badawczych. Jednak wielu badaczy jest nowicjuszami w tym temacie, więc wiąże się to z nieodłącznymi problemami. W tym miejscu opracowujemy praktyczne wprowadzenie do filogenetyki dla laików. Przedstawiamy krok po kroku proces generowania wiarygodnych filogenez na podstawie zestawów danych sekwencji genów. Zaczynamy od podręcznika użytkownika narzędzi do wyszukiwania podobieństw za pośrednictwem interfejsów online, a także lokalnych plików wykonywalnych. Następnie badamy programy do generowania wielu dopasowań sekwencji, a następnie protokoły korzystania z oprogramowania do określania najlepiej dopasowanych modeli ewolucji. Następnie nakreślamy protokoły rekonstrukcji powiązań filogenetycznych za pomocą maksymalnego prawdopodobieństwa i kryteriów bayesowskich, a na koniec opisujemy narzędzia do wizualizacji drzew filogenetycznych. Chociaż w żadnym wypadku nie jest to wyczerpujący opis podejść filogenetycznych, dostarcza czytelnikowi praktycznych informacji wyjściowych na temat kluczowych aplikacji komputerowych powszechnie stosowanych przez filogenetyków. Wizja tego artykułu polega na tym, że może on służyć jako praktyczne narzędzie szkoleniowe dla naukowców rozpoczynających badania filogenetyczne, a także służyć jako zasób edukacyjny, który można włączyć do klasy lub laboratorium dydaktycznego.

Wprowadzenie

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Aby zrozumieć, jak wyewoluowały dwa (lub więcej) gatunków, najpierw konieczne jest uzyskanie danych sekwencyjnych lub morfologicznych z każdej próbki; dane te reprezentują wielkości, które możemy wykorzystać do zmierzenia ich relacji w przestrzeni ewolucyjnej. Podobnie jak w przypadku pomiaru odległości liniowej, posiadanie większej ilości dostępnych danych (np. mil, cali, mikronów) będzie równoznaczne z dokładniejszym pomiarem. Ergo, dokładność, z jaką badacz może wydedukować odległość ewolucyjną, zależy w dużym stopniu od ilości danych informacyjnych dostępnych do pomiaru relacji. Co więcej, ponieważ różne próbki ewoluują w różnym tempie i za pomocą różnych mechanizmów, metoda, której używamy do pomiaru relacji między dwoma taksonami, również bezpośrednio wpływa na dokładność pomiarów ewolucyjnych. Dlatego, ponieważ relacje ewolucyjne nie są bezpośrednio obserwowane, ale zamiast tego są ekstrapolowane na podstawie danych sekwencyjnych lub morfologicznych, problem wnioskowania o związkach ewolucyjnych staje się problemem statystycznym. Filogenetyka to gałąź biologii zajmująca się stosowaniem modeli statystycznych do wzorców ewolucji w celu optymalnej rekonstrukcji historii ewolucji między taksonami. Ta rekonstrukcja między taksonami jest określana jako filogeneza taksonów.

Aby wypełnić lukę w wiedzy specjalistycznej między biologami molekularnymi a biologami ewolucyjnymi, opisujemy tutaj krok po kroku proces wnioskowania filogenezy z zestawu sekwencji. Po pierwsze, szczegółowo opisujemy kroki związane z przeszukiwaniem bazy danych przy użyciu algorytmu Basic Local Alignment Search Tool (BLAST1) za pośrednictwem interfejsu internetowego, a także przy użyciu lokalnych plików wykonywalnych; jest to często pierwszy krok do uzyskania listy sekwencji podobnych do niezidentyfikowanego zapytania, chociaż niektórzy badacze mogą być również zainteresowani gromadzeniem danych dla pojedynczej grupy za pośrednictwem interfejsów internetowych, takich jak Phylota (http://www.phylota.net/). BLAST to algorytm służący do porównywania danych sekwencji pierwszorzędowych aminokwasów lub nukleotydów z bazą danych sekwencji w celu wyszukania "trafień" przypominających sekwencję zapytania. Program BLAST został opracowany przez Stephena Altschula i wsp. w Narodowym Instytucie Zdrowia (NIH)1. Serwer BLAST składa się z wielu różnych programów, a oto lista niektórych z najpopularniejszych programów BLAST:

i) NUKLEOTYD-nukleotyd BLAST (blastn): Ten program wymaga wprowadzenia sekwencji DNA i zwraca najbardziej podobne sekwencje DNA z bazy danych DNA, które użytkownik określa (np. dla konkretnego organizmu).

ii) BLAST białko-białko (blastp): Tutaj użytkownik wprowadza sekwencję białka, a program zwraca najbardziej podobne sekwencje białek z bazy danych białek, które określa użytkownik.

iii) Iteracyjny BLAST SPECYFICZNY DLA POZYCJI (PSI-BLAST) (blastpgp): Dane wejściowe użytkownika to sekwencja białek, która zwraca zestaw blisko powiązanych białek, a z tego zestawu danych generowany jest konserwatywny profil. Następnie generowane jest nowe zapytanie przy użyciu tylko tych konserwatywnych "motywów", które jest używane do przeszukiwania bazy danych białek, a to zwraca większą grupę białek, z których ekstrahuje się nowy zestaw konserwatywnych "motywów", a następnie wykorzystuje się do przeszukiwania bazy danych białek, aż do przestrojenia jeszcze większego zestawu białek i wygenerowania kolejnego profilu i powtórzenia procesu. Uwzględniając powiązane białka do zapytania w każdym kroku, program ten pozwala użytkownikowi zidentyfikować sekwencje, które są bardziej rozbieżne.

iv) Białko translacyjne nukleotydu w 6 ramkach (blastx): Tutaj użytkownik dostarcza dane wejściowe sekwencji nukleotydów, które są przekształcane w sześcioramkowe produkty translacji koncepcyjnej (tj. obie nici) w bazie danych sekwencji białek.

v) Tłumaczenie 6-ramkowe nukleotydów-6-ramkowe translacje nukleotydów (tblastx): Ten program pobiera dane wejściowe sekwencji nukleotydów DNA i tłumaczy je na wszystkie sześcioramkowe produkty translacji koncepcyjnej, które porównuje z sześcioramkowymi translacjami bazy danych sekwencji nukleotydów.

vi) Translacja 6-ramek białko-nukleotyd (tblastn): Ten program wykorzystuje dane wejściowe sekwencji białka do porównania ze wszystkimi sześcioma ramkami odczytu bazy danych sekwencji nukleotydów.

Następnie opisujemy powszechnie używane programy do generowania Multiple Sequence Alignment (MSA) z zestawu danych sekwencji, a następnie dołączamy przewodnik użytkownika do programów, które określają najlepiej dopasowane modele ewolucji dla zestawu danych sekwencji. Rekonstrukcja filogenetyczna jest problemem statystycznym i z tego powodu metody filogenetyczne muszą uwzględniać ramy statystyczne. Ta struktura statystyczna staje się modelem ewolucyjnym, który uwzględnia zmiany sekwencji w zbiorze danych. Ten model ewolucyjny składa się z zestawu założeń dotyczących procesu podstawiania nukleotydów lub aminokwasów, a najlepiej dopasowany model dla konkretnego zestawu danych można wybrać za pomocą testów statystycznych. Dopasowanie do danych różnych modeli można porównać za pomocą testów ilorazu prawdopodobieństwa (LRT) lub kryteriów informacyjnych w celu wybrania najlepiej dopasowanego modelu z zestawu możliwych. Dwa powszechne kryteria informacyjne to kryterium informacyjne Akaickiego (AIC)2 i bayesowskie kryterium informacyjne (BIC)3. Po wygenerowaniu optymalnego wyrównania istnieje wiele różnych metod tworzenia filogenezy na podstawie wyrównanych danych. Istnieje wiele metod wnioskowania o związkach ewolucyjnych; Ogólnie rzecz biorąc, można je podzielić na dwie kategorie: metody oparte na odległości i metody oparte na sekwencjach. Metody oparte na odległości obliczają odległości parami z sekwencji, a następnie używają tych odległości do uzyskania drzewa. Metody oparte na sekwencjach wykorzystują wyrównanie sekwencji bezpośrednio i zwykle przeszukują przestrzeń drzewa przy użyciu kryterium optymalności. Przedstawiamy dwie oparte na sekwencjach metody rekonstrukcji relacji filogenetycznych: są to PhyML4, która implementuje ramy maksymalnego prawdopodobieństwa, oraz MrBayes5, która wykorzystuje bayesowskie wnioskowanie Monte Carlo łańcucha Markowa. Metody prawdopodobieństwa i bayesowskie zapewniają ramy statystyczne dla rekonstrukcji filogenetycznej. Dostarczając użytkownikowi informacji o powszechnie używanych narzędziach do budowania drzew, zapoznajemy czytelnika z niezbędnymi danymi wymaganymi do wnioskowania o pokrewieństwie filogenetycznym.

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Protokół

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Podstawowe narzędzie do wyszukiwania wyrównania lokalnego (BLAST): Interfejs online

  1. Kliknij ten link, aby odwiedzić serwer internetowy BLAST1 w Narodowym Centrum Informacji Biotechnologicznej (NCBI). - http://blast.ncbi.nlm.nih.gov/Blast.cgi (Rysunek 1).
  2. Wprowadź sekwencję tekstową sformatowaną w formacie FASTA (patrz na przykład Rysunek 2) w polu zapytania.
  3. Kliknij odpowiedni program BLAST i odpowiednią bazę danych lub poszczególne gatunki do wykorzystania w wyszukiwaniu, a następnie kliknij "BLAST".
    nuta: Sekwencja sformatowana w formacie FASTA rozpoczyna się od wiersza opisu oznaczonego znakiem ">". Opis musi następować bezpośrednio po znaku ">", sekwencja (tj. nukleotydy lub aminokwasy) jest zgodna z opisem w następnym wierszu. Dane wyjściowe wyszukiwania BLAST są wyświetlane jako HTML, zwykły tekst, XML lub tabele trafień (Text lub csv) z wartością domyślną ustawioną na HTML (Rysunek 3).

2. Podstawowe narzędzie do wyszukiwania wyrównania lokalnego (BLAST): Lokalne pliki wykonywalne

  1. Pobierz najnowsze pliki wykonywalne BLAST z wiersza poleceń BLAST z tego linku:
    ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/ -
  2. em>Dla użytkowników komputerów PC: kliknij dwukrotnie najnowszy plik blast win32.exe i zaakceptuj umowę licencyjną, a następnie kliknij install.
    Uwaga: Domyślnym katalogiem instalacyjnym jest C:\ncbi-blast-2.2.27+.
  3. Skonfiguruj zmienną środowiskową PC w następujący sposób:
    1. Kliknij przycisk "start" komputera, a następnie kliknij prawym przyciskiem myszy "komputer",
    2. Kliknij "Właściwości" i w wyskakującym okienku kliknij zakładkę "zaawansowane"
    3. Kliknij przycisk "Zmienne środowiskowe" i w nowym wyskakującym okienku kliknij przycisk "nowy" w sekcji "Zmienne użytkownika dla użytkownika"
    4. W wyskakującym okienku dodaj nazwę zmiennej "Path" i wartość zmiennej "C:\ncbi-blast-2.2.27+\bin.
      Uwaga: katalog bin zawiera plik wykonywalny (np. blastp itp.).
  4. em>Dla użytkowników komputerów Mac: Otwórz aplikację Terminal (aby to zrobić, po prostu otwórz "Finder" i wyszukaj "Terminal", a spowoduje to wyświetlenie ikony "terminala"). W oknie terminala wpisz:
    >ftp ftp.ncbi.nih.gov
    >Uwaga: może również wpisać adres URL użyty powyżej w przykładzie na PC
  5. Aby uzyskać dostęp do witryny ftp NCBI, wpisz "anonimowy" w polu Nazwa i hasło, a następnie wpisz:
    >cd blast/pliki wykonywalne/NAJNOWSZE
  6. Wymień pliki wykonywalne, wpisując:
    >ls
  7. Pobierz najnowszą wersję, wpisując następujące polecenie (lub inną najnowszą wersję):
    >Skorzystaj ncbi-blast-2.2.7-macosx.tar.gz
  8. Wyjdź z witryny serwera ftp NCBI, wpisując "exit".
  9. Rozpakuj pobrane pliki, wpisując:
    >tar -xzf ncbi-blast-2.2.7-macosx.tar.gz
  10. Dodaj lokalizację plików binarnych dla pliku wykonywalnego blast do swojej ścieżki, aby powłoka mogła przeszukiwać ten katalog podczas szukania poleceń, wpisując:
    >PATH=$PATH:new_folder_location
  11. Sprawdź, czy lokalizacja została dodana do Twojej ścieżki, wpisując:
    > $PATH echa
  12. Pobierz wstępnie sformatowane bazy danych BLAST (które są aktualizowane codziennie), klikając tutaj:
    ftp://ftp.ncbi.nlm.nih.gov/blast/db/
  13. Umieść bazę danych w folderze "db".
  14. em>Na komputerze: otwórz wiersz polecenia MS-DOS (aby to zrobić, kliknij "start" i wpisz "cmd" w pasku wyszukiwania) i zmień katalog na folder ncbi-blast, wpisując:
    C:\Użytkownicy>cd .. \ [przenosi o jeden folder w górę]
    C:\>cd ncbi-blast-2.2.27+
    Spowoduje to zmianę katalogu na:
    C:\ncbi-blast-2.2.27+>
  15. Utwórz bazę danych za pomocą następującego polecenia "makedb":
    >makedb –in db/briggsae.fasta –dbtype prot –out db/briggsae
    >nuta: W poniższym przykładzie (ryc. 4) baza danych nosi nazwę "briggsae" i składa się z jednej grupy powiązań z organizmu Caenorhabditis briggsae.
  16. Utwórz sekwencję białka zapytania o nazwie "test", wstawiając sekwencję tekstową białka sformatowaną w formacie FASTA do folderu "db".
  17. Przeszukaj bazę danych za pomocą wyszukiwania blastp, wpisując następujące polecenie:
    >blastp –query db/test.txt –db db/briggsae –out text.txt
  18. em>Na komputerze Mac: pobierz bazę danych do lokalnego wyszukiwania w Blast, wchodząc na stronę NCBI ftp zgodnie z powyższymi instrukcjami (krok 2.4), a następnie wpisz:
    >lcd .. /Baz danych/
  19. Pobierz interesujący Cię genom lub sekwencję, wpisując:
    >Pobierz NC_[Akces #].fna
    >Uwaga: ".fna" odnosi się do sekwencji nukleotydowej sformatowanej w FADA, a ".faa" odnosi się do sekwencji aminokwasowych sformatowanych w FADA.
  20. Wpisz "quit", aby wyjść ze strony ftp.
  21. Utwórz bazę danych, wpisując:
    >makeblastdb -in db/mouse.faa -out mouse -dbtype prot
  22. Wstaw sekwencję zapytań sformatowaną w formacie FAST do folderu "bin" i przeszukaj bazę danych za pomocą następującego polecenia:
    > blastp -query "twoje zapytanie.fasta" -db "twoja baza danych" -out results.txt

3. Generowanie wielu wyrównań sekwencji

  1. Kliknij poniższe łącza, aby uzyskać dostęp do powszechnie używanych programów do wyrównywania wielu sekwencji (MSA):
    ClustalW6 http://www.clustal.org/
    Kalign7 http://msa.sbc.su.se/cgi-bin/msa.cgi
    MAFFT8,9 http://mafft.cbrc.jp/alignment/software/
    MIĘSIEŃ10 http://www.drive5.com/muscle/
    T-Coffee11 http://www.tcoffee.org/Projects/tcoffee/
    PROBCONS12 http://toolkit.tuebingen.mpg.de/probcons
  2. Kliknij na ten link - http://tcoffee.crg.cat/apps/tcoffee/do:regular - i wprowadź dane sekwencji w formacie FASTA w polu zapytania
    nuta: Próbka wyjściowa z T-Coffee jest widoczna na rysunku 5, podobne pozostałości są oznaczone kolorami.
  3. Pobierz Clustal MSA jako wersję wiersza poleceń (ClustalW) lub wersję graficzną (ClustalX) klikając w ten link: http://www.clustal.org/clustal2/ - następnie kliknij na odpowiedni plik wykonywalny (np. win, Linux, Mac OS X).
  4. Prześlij dane jako tekst sekwencji w formacie FASTA i wyrównaj (Rysunek 6).

4. Określanie najlepiej dopasowanych modeli ewolucji

  1. Kliknij tutaj, aby pobrać program ProtTest13:
    http://darwin.uvigo.es/our-software/
  2. Po pobraniu ProtTest kliknij dwukrotnie plik ProtTest.jar
  3. Po uruchomieniu ProtTest kliknij "wybierz plik" i załaduj dane sekwencji (Rysunek 7).
  4. Następnie kliknij "start", a program rozpocznie się (Rysunek 8).
    nuta: Po zakończeniu przebiegu (Rysunek 8) program wskaże najlepszy model na podstawie kryteriów, np. "Najlepszy model według AIC: WAG+I+G"

5. Wnioskowanie filogenezy opartej na sekwencjach na podstawie maksymalnego prawdopodobieństwa lub wnioskowania bayesowskiego

  1. Pobrano PhyML4 stąd:
    https://code.google.com/p/phyml/
  2. Uruchom plik wykonywalny, klikając dwukrotnie odpowiednią aplikację (np. phyml Windows, phyml Linux itp.), a pojawi się okno interfejsu (Rysunek 9).
  3. Załaduj sekwencję wejściową jako sekwencję sformatowaną przez PHYLIP, wpisując:
    >"nazwa pliku".phy
    >Uwaga: Aby konwertować między formatami sekwencji, użyj programu internetowego "Readseq" dostępnego pod adresem - http://iubio.bio.indiana.edu/cgi-bin/readseq.cgi.
  4. Uruchom program, wpisując "Y".
  5. Pobierz MrBayes5 stąd:
    http://mrbayes.sourceforge.net/download.php
  6. Aby uruchomić program, kliknij plik wykonywalny i wczytaj dane sekwencji sformatowane w formacie NEXUS do programu, wpisując:
    >wykonaj "nazwa pliku".nex
  7. Ustal model ewolucyjny.
  8. Wybierz liczbę generacji do uruchomienia, wpisując:
    >mcmcp ngen = 1000000 [ustawia liczbę pokoleń na 1000000]
    >wypalenie studzienki =10000 [ustawia wypalenie na 10000]
  9. Zapisz długości gałęzi w pliku wyników, wpisując:
    >MCMCP Savebrlens = Tak
  10. Uruchom analizę, wpisując:
    >mcmc
  11. powiedział:
  12. Podsumuj drzewa za pomocą polecenia "sumt".

6. Wizualizacja filogenezy

  1. Zobacz listę programów do przeglądania drzew tutaj:
    http://www.treedyn.org/overview/editors.html
  2. Pobierz program TreeView14 tutaj:
    http://taxonomy.zoology.gla.ac.uk/rod/treeview.html

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Wyniki

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Znalezienie podobieństw do zapytania pozwala badaczom przypisać potencjalną tożsamość nowym sekwencjom, a także wnioskować o związkach między sekwencjami. Typem wejściowym pliku dla BLAST1 jest sekwencja tekstu w formacie FASTA lub numer dostępu GenBank. Sekwencja w formacie FASTA rozpoczyna się od wiersza opisu oznaczonego znakiem ">" (Rysunek 2). Opis musi następować bezpośrednio po znaku ">", sekwencja (tj. nukleotydy lub aminokwasy) jest zgodna z opisem w następnym wierszu. ...

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Dyskusja

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Mamy nadzieję, że ten artykuł posłuży jako punkt wyjścia do prowadzenia badaczy lub studentów, którzy są nowicjuszami w filogenetyce. Projekty sekwencjonowania genomu stały się tańsze w ciągu ostatnich kilku lat, w wyniku czego rośnie zapotrzebowanie użytkowników na tę technologię, a obecnie produkcja dużych zbiorów danych sekwencji jest powszechna w małych laboratoriach. Te zestawy danych często dostarczają naukowcom zestawów genów, które wymagają ram filogenetycznych, aby zacząć rozumieć ich funkcję. Ponadto, ponieważ ...

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Oświadczenia

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Nie mamy nic do ujawnienia.

Podziękowania

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Dziękujemy członkom laboratorium O'Halloran za komentarze do manuskryptu. Dziękujemy Wydziałowi Nauk Biologicznych Uniwersytetu George'a Washingtona i Columbian College of Arts and Sciences za dofinansowanie dla D. O'Hallorana.

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Materiały

Lista materiałów użytych w tym artykule
NazwaFirmaNumer katalogowyKomentarze
Strona internetowa BLAST http://blast.ncbi.nlm.nih.gov/Blast.cgi
BLAST pliki wykonywalne ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/
Wstępnie sformatowane bazy danychBLASTftp://ftp.ncbi.nlm.nih.gov/blast/db/
Clustalhttp://www.clustal.org/
Kalignhttp://msa.sbc.su.se/cgi-bin/msa.cgi
MAFFThttp://mafft.cbrc.jp/alignment/software/
MUSCLEhttp://www.drive5.com/muscle/
T-Coffeehttp://www.tcoffee.org/Projects/tcoffee/
PROBCONSa href="http://toolkit.tuebingen.mpg.de/probcons">http://toolkit.tuebingen.mpg.de/probcons 
Se-Al http://tree.bio.ed.ac.uk/software/seal/
BSEdit http://www.bsedit.org/
JalViewhttp://www.jalview.org/
SeaViewhttp://pbil.univ-lyon1.fr/software/seaview.html
ProtTest& nbsp;https://code.google.com/p/prottest3/
Java Runtime http://www.java.com/en/download/chrome.jsp
Readseqhttp://iubio.bio.indiana.edu/cgi-bin/readseq.cgi
jModelTesthttps://code.google.com/p/jmodeltest2/
PhyMLa href="https://code.google.com/p/phyml/">https://code.google.com/p/phyml/
MrBayeshttp://mrbayes.sourceforge.net/download.php
TreeViewhttp://taxonomy.zoology.gla.ac.uk/rod/treeview.html
TreeDynhttp://www.treedyn.org/
<<

Bibliografia

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. Altschul, S. F., Carroll, R. J., Lipman, D. J. Weights for data related by a tree. J. Mol. Biol. 207 (4), 647-653 (1989).
  2. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automat. Contr. 19 (6), 706-723 (1974).
  3. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 6 (2), 461-464 (1978).
  4. Guindon, S., Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52 (5), 696-704 (2003).
  5. Huelsenbeck, J. P., Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 17 (8), 754-755 (2001).
  6. Thompson, J. D., Higgins, D. G., Gibson, T. J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22 (22), 4673-4680 (1994).
  7. Lassmann, T., Sonnhammer, E. L. Kalign--an accurate and fast multiple sequence alignment algorithm. BMC Bioinformatics. 6, 298 (2005).
  8. Katoh, K., Kuma, K., Toh, H., Miyata, T. MAFFT version 5: Improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 33 (2), 511-518 (2005).
  9. Katoh, K., Misawa, K., Kuma, K., Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 30 (14), 3059-3066 (2002).
  10. Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32 (5), 1792-1797 (2004).
  11. Notredame, C., Higgins, D. G., Heringa, J. T-coffee: A novel method for fast and accurate multiple sequence alignment. J. Mol. Biol. 302 (1), 205-217 (2000).
  12. Do, C. B., Mahabhashyam, M. S., Brudno, M., Batzoglou, S. ProbCons: Probabilistic consistency-based multiple sequence alignment. Genome Res. 15 (2), 330-340 (2005).
  13. Darriba, D., Taboada, G. L., Doallo, R., Posada, D. ProtTest 3: Fast selection of best-fit models of protein evolution. Bioinformatics. 27 (8), 1164-1165 (2011).
  14. Page, R. D. TreeView: An application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12 (4), 357-358 (1996).
  15. Darriba, D., Taboada, G. L., Doallo, R., Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods. 9 (8), 772 (2012).
  16. Chevenet, F., Brun, C., Banuls, A. L., Jacq, B., Christen, R. TreeDyn: Towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics. 7, 439 (2006).

Dostęp ograniczony. Zaloguj się lub rozpocznij wersję próbną, aby wyświetlić tę treść.

Przedruki i uprawnienia

Poproś o pozwolenie na ponowne wykorzystanie tekstu lub ilustracji tego artykułu JoVE

Poproś o pozwolenie

Tagi

PhylogeneticsSequence AlignmentBlast SearchMaximum LikelihoodBayesian InferenceMultiple Sequence AlignmentTree VisualizationModel SelectionEvolutionary AnalysisGene Sequences

Powiązane artykuły