Rebecca B. Harris

The Influence of Gene Flow on Species Tree Estimation: A Simulation Study

Gene flow among populations or species and incomplete lineage sorting (ILS) are two evolutionary processes responsible for generating gene tree discordance and therefore hindering species tree estimation. Numerous studies have evaluated the impacts of ILS on species tree inference, yet the ramifications of gene flow on species trees remain less studied. Here, we simulate and analyse multilocus sequence data generated with ILS and gene flow to quantify their impacts on species tree inference. We characterize species tree estimation errors under various models of gene flow, such as the isolation-migration model, the n-island model, and gene flow between non-sister species or involving ancestral species, and species boundaries crossed by a single gene copy (allelic introgression) or by a single migrant individual. These patterns of gene flow are explored on species trees of different sizes (4 vs. 10 species), at different time scales (shallow vs. deep), and with different migration rates. Species trees are estimated with the multispecies coalescent model using Bayesian methods (BEST and *BEAST) and with a summary statistic approach (MPEST) that facilitates phylogenomic-scale analysis. Even in cases where the topology of the species tree is estimated with high accuracy, we find that gene flow can result in overestimates of population sizes (species tree dilation) and underestimates of species divergence times (species tree compression). Signatures of migration events remain present in the distribution of coalescent times for gene trees, and with sufficient data it is possible to identify those loci that have crossed species boundaries. These results highlight the need for careful sampling design in phylogeographic and species delimitation studies as gene flow, introgression, or incorrect sample assignments can bias the estimation of the species tree topology and of parameter estimates such as population sizes and divergence times.

Comparative Species Divergence across Eight Triplets of Spiny Lizards (Sceloporus) Using Genomic Sequence Data

Species divergence is typically thought to occur in the absence of gene flow, but many empirical studies are discovering that gene flow may be more pervasive during species formation. Although many examples of divergence with gene flow have been identified, few clades have been investigated in a comparative manner, and fewer have been studied using genome-wide sequence data. We contrast species divergence genetic histories across eight triplets of North American Sceloporus lizards using a maximum likelihood implementation of the isolation–migration (IM) model. Gene flow at the time of species divergence is modeled indirectly as variation in species divergence time across the genome or explicitly using a migration rate parameter. Likelihood ratio tests (LRTs) are used to test the null model of no gene flow at speciation against these two alternative gene flow models. We also use the Akaike information criterion to rank the models. Hundreds of loci are needed for the LRTs to have statistical power, and we use genome sequencing of reduced representation libraries to obtain DNA sequence alignments at many loci (between 340 and 3,478; mean = 1,678) for each triplet. We find that current species distributions are a poor predictor of whether a species pair diverged with gene flow. Interrogating the genome using the triplet method expedites the comparative study of species divergence history and the estimation of genetic parameters associated with speciation.

apex: phylogenetics with multiple genes

Genetic sequences of multiple genes are becoming increasingly common for a wide range of organisms including viruses, bacteria and eukaryotes. While such data may sometimes be treated as a single locus, in practice, a number of biological and statistical phenomena can lead to phylogenetic incongruence. In such cases, different loci should, at least as a preliminary step, be examined and analysed separately. The r software has become a popular platform for phylogenetics, with several packages implementing distance‐based, parsimony and likelihood‐based phylogenetic reconstruction, and an even greater number of packages implementing phylogenetic comparative methods. Unfortunately, basic data structures and tools for analysing multiple genes have so far been lacking, thereby limiting potential for investigating phylogenetic incongruence. In this study, we introduce the new r package apex to fill this gap. apex implements new object classes, which extend existing standards for storing DNA and amino acid sequences, and provides a number of convenient tools for handling, visualizing and analysing these data. In this study, we introduce the main features of the package and illustrate its functionalities through the analysis of a simple data set.

Evidence for concerted movement of nuclear and mitochondrial clines in a lizard hybrid zone

Moving hybrid zones provide compelling examples of evolution in action, yet long‐term studies that test the assumptions of hybrid zone stability are rare. Using replicated transect samples collected over a 10‐year interval from 2002 to 2012, we find evidence for concerted movement of genetic clines in a plateau fence lizard hybrid zone (Sceloporus tristichus) in Arizona. Cline‐fitting analyses of SNP and mtDNA data both provide evidence that the hybrid zone shifted northward by approximately 2 km during the 10‐year interval. For each sampling period, the mtDNA cline centre is displaced from the SNP cline centre and maintaining an introgression distance of approximately 3 km. The northward expansion of juniper trees into the Little Colorado River Basin in the early 1900s provides a plausible mechanism for hybrid zone formation and movement, and a broadscale quantification of recent land cover change provides support for increased woody species encroachment at the southern end of the hybrid zone. However, population processes can also contribute to hybrid zone movement, and the current stability of the ecotone habitats in the centre of the hybrid zone suggests that movement could decelerate in the future.

Rape and the Prevalence of Hybrids in Broadly Sympatric Species: a Case Study using Albatrosses

Conspecific rape often increases male reproductive success. However, the haste and aggression of forced copulations suggests that males may sometimes rape heterospecific females, thus making rape a likely, but undocumented, source of hybrids between broadly sympatric species. We present evidence that heterospecific rape may be the source of hybrids between Black-footed and Laysan Albatrosses (Phoebastria nigripes, and P. immutabilis, respectively). Extensive field studies have shown that paired (but not unpaired) males of both of these albatross species use rape as a supplemental reproductive strategy. Between species differences in size, timing of laying, and aggressiveness suggest that Black-footed Albatrosses should be more successful than Laysan Albatrosses in heteropspecific rape attempts, and male Black-footed Albatrosses have been observed attempting to force copulations on female Laysan Albatrosses. Nuclear markers showed that the six hybrids we studied were F1s and mitochondrial markers showed that male Black-footed Albatrosses sired all six hybrids. Long-term gene exchange between these species has been from Black-footed Albatrosses into Laysan Albatrosses, suggesting that the siring asymmetry found in our hybrids has long persisted. If hybrids are sired in heterospecific rapes, they presumably would be raised and sexually imprinted on Laysan Albatrosses, and two unmated hybrids in a previous study courted only Laysan Albatrosses.

Discordance between genomic divergence and phenotypic variation in a rapidly evolving avian genus (Motacilla)

Generally, genotypes and phenotypes are expected to be spatially congruent; however, in widespread species complexes with few barriers to dispersal, multiple contact zones, and limited reproductive isolation, discordance between phenotypes and phylogeographic groups is more probable. Wagtails (Motacilla) are a genus of birds with striking plumage pattern variation across the Old World. Up to 13 subspecies are recognized within a single species, yet previous studies using mitochondrial DNA have supported polyphyletic phylogeographic groups that are inconsistent with subspecies plumage characteristics. In this study, we investigate the link between phenotypes and genotype by taking a phylogenetic approach. We use genome-wide SNPs, nuclear introns, and mitochondrial DNA to estimate population structure, isolation by distance, and species relationships. Together, our genetic sampling includes complete species-level sampling and comprehensive coverage of the three most phenotypically diverse Palearctic species. Our study provides strong evidence for species-level patterns of differentiation, however population-level differentiation is less pronounced. SNPs provide a robust estimate of species-level relationships, which are mostly corroborated by a combined analysis of mtDNA and nuclear introns (the first time-calibrated species tree for the genus). However, the mtDNA tree is strongly incongruent and is considered to misrepresent the species phylogeny. The extant wagtail lineages originated during the Pliocene and the Eurasian lineage underwent rapid diversification during the Pleistocene. Three of four widespread Eurasian species exhibit an east-west divide that contradicts both subspecies taxonomy and phenotypic variation. Indeed, SNPs fail to distinguish between phenotypically distinct subspecies within the M. alba and M. flava complexes, and instead support geographical regions, each of which is home to two or more different looking subspecies. This is a major step towards our understanding of wagtail phylogeny compared to previous analyses of fewer species and considerably less sequence data.