Wilke:Molecular Evolution

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
(Sequence Alignment)
(Sequence Alignment)
Line 13: Line 13:
<pre>
<pre>
mafft --auto infile > outfile
mafft --auto infile > outfile
 +
</pre>
 +
 +
 +
== Phylogenetic Inference ==
 +
Several options exist for creating a phylogeny. One can either infer a Maximum Likelihood or Bayesian tree. Accuracy is comparable between these two methods, but maximum likelihood inference may be faster. We recommend using the software RAxML for maximum likelihood phylogenetic inference. In the event that you have a large number of sequences (>500), we recommend using FastTree, which employs maximum likelihood techniques but has significant speed improvements.
 +
 +
'''RAxML Usage'''
 +
* Requires input alignment in phylip format
 +
* A model must always be specified using the -m option. While many options are available, for nucleotide data we recommend "GTRGAMMA" and for amino acid data "PROTJTTGAMMA" or "PROTWAGGAMMA." The following examples will all use GTRGAMMA, but note that other models are possible.
 +
 +
* Basic inference, yielding a single tree inference.
 +
<pre>
 +
raxmlHPC -m GTRGAMMA -s <infile> -n <outfile>
</pre>
</pre>

Revision as of 15:06, 10 July 2013

The Basics

Evolutionary rate (dN/dS) analyses involve a basic pipeline: sequence alignment, phylogenetic inference, and finally evolutionary rate inference. When dealing with protein coding sequences, always align using amino acid data in order to preserve codons. Then, back-translate into nucleotide data, as is required for the final step in the pipeline. Phylogenies may be made either with amino acid or nucleotide data, although an amino acid tree may be slightly more accurate.

Additionally note that a minimum of 10 sequences are recommended to achieve well-supported results in an evolutionary rates analysis.


Sequence Alignment

Commonly aligners include mafft, muscle, and prank. While prank is probably the most accurate, it is also very time consuming. We recommend using mafft for sequence alignments.

To align sequences in mafft, we additionally recommend using the "--auto" option. This will allow mafft to select the optimal alignment algorithm to use on your data. Mafft accepts a variety of file formats, including fasta and phylip (sequential and/or interleaved).

The infile should contain unaligned amino acid sequences, and aligned sequences will be sent to the outfile name provided.

mafft --auto infile > outfile


Phylogenetic Inference

Several options exist for creating a phylogeny. One can either infer a Maximum Likelihood or Bayesian tree. Accuracy is comparable between these two methods, but maximum likelihood inference may be faster. We recommend using the software RAxML for maximum likelihood phylogenetic inference. In the event that you have a large number of sequences (>500), we recommend using FastTree, which employs maximum likelihood techniques but has significant speed improvements.

RAxML Usage

  • Requires input alignment in phylip format
  • A model must always be specified using the -m option. While many options are available, for nucleotide data we recommend "GTRGAMMA" and for amino acid data "PROTJTTGAMMA" or "PROTWAGGAMMA." The following examples will all use GTRGAMMA, but note that other models are possible.
  • Basic inference, yielding a single tree inference.
raxmlHPC -m GTRGAMMA -s <infile> -n <outfile>
Personal tools