Drosophila 25 species phylogeny

We used the peptide sequences for 25 Drosophila species along with an outgroup species, <i>Musca domestica</i>, that were available on NCBI (ftp://ftp.ncbi.nlm.nih.gov/genomes/) and FlyBase (ftp://flybase.net/12_species_analysis/genomes/) as of April 21, 2017. We employed the clustering algorithm MCL to predict orthologous groups of peptides and used the 1,634 ortho-groups that were predicted as being single-copy in all species for phylogeny reconstruction. These groups were aligned with PASTA and the alignments were concatenated and given to RAxML with the PROTGAMMAJTTF model to estimate the species topology and relative branch lengths. Other methods based on individual gene trees (average consensus and ASTRAL) resulted in the exact same topology. We then used r8s to estimate branch lengths in terms of millions of years with three calibration points: 49-59 mya for the common ancestor of <i>D. pseudoobscura</i> and <i>D. melanogaster</i>, 64-74 mya for the common ancestor of <i>D. grimshawi</i> and <i>D. melanogaster</i>, and 156 mya for the common ancestor of <i>D. melanogaster</i> and <i>Musca domestica</i>.<br><br>File summary:<br>1. drosophila-25spec-concordance-factor-tree.png<br>An image of the phylogeny with branches scaled by relative number of substitutions. Nodes are labeled with concordance factors. Concordance factors are calculated for each node in the species phylogeny as: (# of gene trees with that node / total # of gene trees).<br><br>2. drosophila-25spec-node-age-tree.png<br>An image of the phylogeny with branches scaled by millions of years. Nodes are labeled with their ages in millions of years.<br><br>3. drosophila-25spec-time-tree.tre<br>The Newick string of the phylogeny with branches scaled by millions of years. Nodes are labeled with concordance factors.<br><br>4. drosophila-25spec-tree.tre<br>The Newick string of the phylogeny with branches scaled by relative number of substitutions. Nodes are labeled with concordance factors.<br>