Abstract:Abstract 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.]

Modelling gene content across a phylogeny to determine when genes become associated

Modeling gene content across a phylogeny to determine when genes become associated

The Inference of Gene Trees with Species Trees

Determining the evolutionary history of gene families

Counting and sampling gene family evolutionary histories in the duplication-loss and duplication-loss-transfer models

A Probabilistic Model for Gene Content Evolution with Duplication, Loss, and Horizontal Transfer

A tale of too many trees: a conundrum for phylogenetic regression

Computing the probability of gene trees concordant with the species tree in the multispecies coalescent

PhyloAcc-GT: A Bayesian method for inferring patterns of substitution rate shifts on targeted lineages accounting for gene tree discordance

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Genes with Minimal Phylogenetic Information Are Problematic for Coalescent Analyses when Gene Tree Estimation is Biased

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Anomaly zones for uniformly sampled gene trees under the gene duplication and loss model

Probabilities of unranked and ranked anomaly zones under birth-death models

The Theory of Gene Family Histories

Bayesian Inference of Species Trees from Multilocus Data

Linking phenotypic and genotypic variation: a relaxed phylogenetic approach using the probabilistic programming language Stan

Bayesian hidden Markov tree models for clustering genes with shared evolutionary history

Maximum-frequency gene tree: a simplified genome-scale approach to overcoming incongruence in molecular phylogenies

Properties of Consensus Methods for Inferring Species Trees from Gene Trees

Complex Models of Sequence Evolution Improve Fit, but not Gene Tree Discordance, for Tetrapod Mitogenomes