Colonization and/or mitochondrial selective sweeps across the North Atlantic intertidal assemblage revealed by multi‐taxa approximate Bayesian computation

Intertidal and subtidal communities of the western and eastern coasts of the North Atlantic Ocean were greatly affected by Pleistocene glaciations, with some taxa persisting on both coasts, and others recolonizing after being extirpated on one coast during the Last Glacial Maximum. In the original spirit of comparative phylogeography, we conducted a comparative analysis using mtDNA sequence data and a hierarchical approximate Bayesian computation (ABC) approach for testing these two scenarios across 12 intertidal and subtidal coastal invertebrates spanning the North Atlantic to determine the temporal dynamics of species membership of these two ephemeral communities. Conditioning on a low gene‐flow model, our results suggested that a colonization or mitochondrial selective sweep history was predominant across all taxa, with only the bivalve mollusc Mytilus edulis showing a history of trans‐Atlantic persistence. Conditioning on a high gene‐flow model weakened the support for this assemblage‐level demographic history. The predominance of a colonization‐type history also highlights concerns about analyses based on single‐locus data where genetic hitchhiking may be incorrectly inferred as colonization. In conclusion, driving factors in shifting species range distributions and membership of ephemeral coastal communities could be species‐specific environmental tolerances, species interactions, and/or stochastic demographic extinction. Through a re‐examination of a long‐standing question of North Atlantic phylogeography, we highlight the flexibility and statistical honesty of using a model‐based ABC approach.     

Comparative phylogeographic summary statistics for testing simultaneous vicariance

Testing for simultaneous vicariance across comparative phylogeographic data sets is a notoriously difficult problem hindered by mutational variance, the coalescent variance, and variability across pairs of sister taxa in parameters that affect genetic divergence. We simulate vicariance to characterize the behaviour of several commonly used summary statistics across a range of divergence times, and to characterize this behaviour in comparative phylogeographic datasets having multiple taxon-pairs. We found Tajima's D to be relatively uncorrelated with other summary statistics across divergence times, and using simple hypothesis testing of simultaneous vicariance given variable population sizes, we counter-intuitively found that the variance across taxon pairs in Nei and Li's net nucleotide divergence (pi(net)), a common measure of population divergence, is often inferior to using the variance in Tajima's D across taxon pairs as a test statistic to distinguish ancient simultaneous vicariance from variable vicariance histories. The opposite and more intuitive pattern is found for testing more recent simultaneous vicariance, and overall we found that depending on the timing of vicariance, one of these two test statistics can achieve high statistical power for rejecting simultaneous vicariance, given a reasonable number of intron loci (> 5 loci, 400 bp) and a range of conditions. These results suggest that components of these two composite summary statistics should be used in future simulation-based methods which can simultaneously use a pool of summary statistics to test comparative the phylogeographic hypotheses we consider here.     

Coalescent-based hypothesis testing supports multiple Pleistocene refugia in the Pacific Northwest for the Pacific giant salamander (Dicamptodon tenebrosus)

Phylogeographic patterns of many taxa are explained by Pleistocene glaciation. The temperate rainforests within the Pacific Northwest of North America provide an excellent example of this phenomenon, and competing phylogenetic hypotheses exist regarding the number of Pleistocene refugia influencing genetic variation of endemic organisms. One such endemic is the Pacific giant salamander, Dicamptodon tenebrosus. In this study, we estimate this species' phylogeny and use a coalescent modeling approach to test five hypotheses concerning the number, location and divergence times of purported Pleistocene refugia. Single refugium hypotheses include: a northern refugium in the Columbia River Valley and a southern refugium in the Klamath-Siskiyou Mountains. Dual refugia hypotheses include these same refugia but separated at varying times: last glacial maximum (20,000 years ago), mid-Pleistocene (800,000 years ago) and early Pleistocene (1.7 million years ago). Phylogenetic analyses and inferences from nested clade analysis reveal distinct northern and southern lineages expanding from the Columbia River Valley and the Klamath-Siskiyou Mountains, respectively. Results of coalescent simulations reject both single refugium hypotheses and the hypothesis of dual refugia with a separation date in the late Pleistocene but not hypotheses predicting dual refugia with separation in early or mid-Pleistocene. Estimates of time since divergence between northern and southern lineages also indicate separation since early to mid-Pleistocene. Tests for expanding populations using mismatch distributions and 'g' distributions reveal demographic growth in the northern and southern lineages. The combination of these results provides strong evidence that this species was restricted into, and subsequently expanded from, at least two Pleistocene refugia in the Pacific Northwest.     

Ginkgo: spatially-explicit simulator of complex phylogeographic histories

We present Ginkgo, a software package for agent-based, forward-time simulations of genealogies of multiple unlinked loci from diploid populations. Ginkgo simulates the evolution of one or more species on a spatially explicit landscape of cells. The user of the software can specify the geographical and environmental characteristics of the landscape, and these properties can change according to a prespecified schedule. The geographical elements modelled include the arrangement of cells and movement rates between particular cells. Each species has a function that can calculate a fitness score for any combination of an individual organism's phenotype and environmental characteristics. The user can control the number of fitness factors (the dimensionality of the cell-specific fitness factors and the individuals phenotypic vectors) and the weighting of each of these dimensions in the fitness calculation. Cell-specific fitness trait optima can be specified across the landscape to mimic differences in habitat. In addition to their differing fitness functions, species can differ in terms of their vagility and fecundity. Genealogies and occurrence data can be produced at any time during the simulation in NEXUS and ESRI Ascii Grid formats, respectively.     

The anatomy of a ‘suture zone’ in Amazonian butterflies: a coalescent‐based test for vicariant geographic divergence and speciation

Attempts by biogeographers to understand biotic diversification in the Amazon have often employed contemporary species distribution patterns to support particular theories, such as Pleistocene rainforest refugia, rather than to test among alternative hypotheses. Suture zones, narrow regions where multiple contact zones and hybrid zones between taxa cluster, have been seen as evidence for past expansion of whole biotas that have undergone allopatric divergence in vicariant refuges. We used coalescent analysis of mutilocus sequence data to examine population split times in 22 pairs of geminate taxa in ithomiine and heliconiine butterflies. We test a hypothesis of simultaneous divergence across a suture zone in NE Peru. Our results reveal a scattered time course of diversification in this suture zone, rather than a tight cluster of split times. Additionally, we find rapid diversification within some lineages such as Melinaea contrasting with older divergence within lineages such as the Oleriina (Hyposcada and Oleria). These results strongly reject simple vicariance as a cause of the suture zone. At the same time, observed lineage effects are incompatible with a series of geographically coincident vicariant events which should affect all lineages similarly. Our results suggest that Pleistocene climatic forcing cannot readily explain this Peruvian suture zone. Lineage‐specific biological traits, such as characteristic distances of gene flow or varying rates of parapatric divergence, may be of greater importance.     

叶绿体DNA标记在谱系地理学中的应用研究进展

谱系地理学研究旨在探究历史上发生的影响目前遗传谱系系统发育和空间分布格局关系的生态与进化过程。叶绿体DNA具有单亲遗传、低突变率、单倍体等特征, 其分子标记不同程度地保留着植物长期进化的历史遗传痕迹, 有助于深度解析谱系地理变异的形成机制。本文探讨了上述特征是怎样影响分子标记的选择、扩大或缩小群体遗传结构分化、延长或缩短空间基因溯祖时间、促进或阻碍种间基因渐渗及谱系分选(复系、并系和单系形成)进程, 重点阐述了这些影响过程的理论基础, 并结合实际例子阐述谱系地理研究进展。由于位点间在突变率、选择强度及它们与漂变互作等方面存在异质性, 今后一个研究重点就是基于叶绿体全基因组序列分析谱系地理变化格局, 包括分析DNA位点间的基因渐渗或基因流动程度差异分布及沿着叶绿体DNA序列上谱系分选差异分布。     

A species assemblage approach to comparative phylogeography of birds in southern Australia

We present a novel approach to investigating the divergence history of biomes and their component species using single-locus data prior to investing in multilocus data. We use coalescent-based hierarchical approximate Bayesian computation (HABC) methods (MsBayes) to estimate the number and timing of discrete divergences across a putative barrier and to assign species to their appropriate period of co-divergence. We then apply a coalescent-based full Bayesian model of divergence (IMa) to suites of species shown to have simultaneously diverged. The full Bayesian model results in reduced credibility intervals around divergence times and allows other parameters associated with divergence to be summarized across species assemblages. We apply this approach to 10 bird species that are wholly or patchily discontinuous in semi-arid habitats between Australia's southwest (SW) and southeast (SE) mesic zones. There was substantial support for up to three discrete periods of divergence. HABC indicates that two species wholly restricted to more mesic habitats diverged earliest, between 594,382 and 3,417,699 years ago, three species from semi-arid habitats diverged between 0 and 1,508,049 years ago, and four diverged more recently, between 0 and 396,843 years ago. Eight species were assigned to three periods of co-divergence with confidence. For full Bayesian analyses, we accounted for uncertainty in the two remaining species by analyzing all possible suites of species. Estimates of divergence times from full Bayesian divergence models ranged between 429,105 and 2,006,355; 67,172 and 663,837; and 24,607 and 171,085 for the earliest, middle, and most recent periods of co-divergence, respectively. This single-locus approach uses the power of multitaxa coalescent analyses as an efficient means of generating a foundation for further, targeted research using multilocus and genomic tools applied to an understudied biome.     

Beringian origins and cryptic speciation events in the fly agaric (Amanita muscaria)

Amanita muscaria sensu lato has a wide geographic distribution, occurring in Europe, Asia, Africa, Australia, New Zealand, and North, Central and South America. Previous phylogenetic work by others indicates three geographic clades (i.e. 'Eurasian', 'Eurasian-alpine' and 'North American' groups) within A. muscaria. However, the historical dispersal patterns of A. muscaria remained unclear. In our project, we collected specimens from arctic, boreal and humid temperate regions in Alaska, and generated DNA sequence data from the protein-coding beta-tubulin gene and the internal transcribed spacer (ITS) and large subunit (LSU) regions of the ribosomal DNA repeat. Homologous sequences from additional A. muscaria isolates were downloaded from GenBank. We conducted phylogenetic and nested clade analyses (NCA) to reveal the phylogeographic history of the species complex. Although phylogenetic analyses confirmed the existence of the three above-mentioned clades, representatives of all three groups were found to occur sympatrically in Alaska, suggesting that they represent cryptic phylogenetic species with partially overlapping geographic distributions rather than being allopatric populations. All phylogenetic species share at least two morphological varieties with other species, suggesting ancestral polymorphism in pileus and wart colour pre-dating their speciations. The ancestral population of A. muscaria likely evolved in the Siberian-Beringian region and underwent fragmentation as inferred from NCA and the coalescent analyses. The data suggest that these populations later evolved into species, expanded their range in North America and Eurasia. In addition to range expansions, populations of all three species remained in Beringia and adapted to the cooling climate.