Dataset Information

ASTRAL: genome-scale coalescent-based species tree estimation.

ABSTRACT:

Motivation

Species trees provide insight into basic biology, including the mechanisms of evolution and how it modifies biomolecular function and structure, biodiversity and co-evolution between genes and species. Yet, gene trees often differ from species trees, creating challenges to species tree estimation. One of the most frequent causes for conflicting topologies between gene trees and species trees is incomplete lineage sorting (ILS), which is modelled by the multi-species coalescent. While many methods have been developed to estimate species trees from multiple genes, some which have statistical guarantees under the multi-species coalescent model, existing methods are too computationally intensive for use with genome-scale analyses or have been shown to have poor accuracy under some realistic conditions.

Results

We present ASTRAL, a fast method for estimating species trees from multiple genes. ASTRAL is statistically consistent, can run on datasets with thousands of genes and has outstanding accuracy-improving on MP-EST and the population tree from BUCKy, two statistically consistent leading coalescent-based methods. ASTRAL is often more accurate than concatenation using maximum likelihood, except when ILS levels are low or there are too few gene trees.

Availability and implementation

ASTRAL is available in open source form at https://github.com/smirarab/ASTRAL/. Datasets studied in this article are available at http://www.cs.utexas.edu/users/phylo/datasets/astral.

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Mirarab S

PROVIDER: S-EPMC4147915 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Taxonomic identifications in some groups of lichen-forming fungi have been challenge largely due to the scarcity of taxonomically relevant features and limitations of morphological and chemical characters traditionally used to distinguish closely related taxa. Delineating species boundaries in closely related species or species complexes often requires a range of multisource data sets and comprehensive analytical methods. Here we aim to examine species boundaries in a group of saxicolous lichen forming fungi, the Aspiciliella intermutans complex (Megasporaceae), widespread mainly in the Mediterranean. We gathered DNA sequences of the nuclear ribosomal internal transcribed spacer (nuITS), the nuclear large subunit (nuLSU), the mitochondrial small subunit (mtSSU) ribosomal DNA, and the DNA replication licensing factor MCM7 from 80 samples mostly from Iran, Caucasia, Greece and eastern Europe. We used a combination of phylogenetic strategies and a variety of empirical, sequence-based species delimitation approaches to infer species boundaries in this group. The latter included: the automatic barcode gap discovery (ABGD), the multispecies coalescent approach *BEAST and Bayesian Phylogenetics and Phylogeography (BPP) program. Different species delimitation scenarios were compared using Bayes factors species delimitation analysis. Furthermore, morphological, chemical, ecological and geographical features of the sampled specimens were examined. Our study uncovered cryptic species diversity in A. intermutans and showed that morphology-based taxonomy may be unreliable, underestimating species diversity in this group of lichens. We identified a total of six species-level lineages in the A. intermutans complex using inferences from multiple empirical operational criteria. We found little corroboration between morphological and ecological features with our proposed candidate species, while secondary metabolite data do not corroborate tree topology. The present study on the A. intermutans species-complex indicates that the genus Aspiciliella, as currently circumscribed, is more diverse in Eurasia than previously expected.

Project description:Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs-paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.].