Genetic signatures of lineage fusion closely resemble population decline

Abstract

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F _S and Z _nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Keywords: Bayesian skyline plot; demographic history; fission‐fusion; neutrality test; phylogeography; population merging.