Phylogeographic breaks without geographic barriers to gene flow

Abstract.— The spatial distribution of genetic markers can be useful both in estimating patterns of gene flow and in reconstructing biogeographic history, particularly when gene genealogies can be estimated. Genealogies based on nonrecombining genetic units such as mitochondrial and chloroplast DNA often consist of geographically separated clades that come into contact in narrow regions. Such phylogeographic breaks are usually assumed to be the result of long-term barriers to gene flow. Here I show that deep phylogeographic breaks can form within a continuously distributed species even when there are no barriers to gene flow. The likelihood of observing phylogeographic breaks increases as the average individual dispersal distance and population size decrease. Those molecular markers that are most likely to show evidence of real geographic barriers are also most likely to show phylogeographic breaks that formed without any barrier to gene flow. These results might provide an explanation as to why some species, such as the greenish warblers ( Phylloscopus trochiloides ), have phylogeographic breaks in mitochondrial or chloroplast DNA that do not coincide with sudden changes in other traits.     

Phylogeographic structure, gene flow and species status in blue grouse (Dendragapus obscurus)

We investigated the genetic population structure and species status of a relatively sedentary bird that is a permanent resident of western North American forests, the blue grouse (Dendragapus obscurus). Phylogenetic analysis of complete mitochondrial control region DNA sequences resulted in the identification of three basal clades of haplotypes that were largely congruent with well-known biogeographical regions. These clades corresponded to the parapatric sooty (D. o. fuliginosus) and dusky (D. o. obscurus) subspecies groups of blue grouse plus a previously unrecognized division between northern and southern dusky grouse populations; the latter does not correspond closely to any currently recognized subspecies boundary. Approximately 66% of the total genetic variation was distributed among these three regions. Maximum likelihood estimates of gene flow between the regions were low or asymmetric; gene flow has been insufficient to prevent genetic divergence between dusky and sooty grouse. Estimates of gene flow among populations within sooty grouse were large except across the Columbia River valley. Among populations of dusky grouse, estimates of gene flow were heterogeneous and asymmetrical, reflecting large-scale fragmentation of the distribution due to landscape features and associated vegetation. Genetic, morphological and behavioural evidence suggest that sooty and dusky grouse are species-level taxa; the specific status of a third clade remains ambiguous.     

Revealing the demographic histories of species using DNA sequences

Various methodological approaches using molecular sequence data have been developed and applied across several fields, including phylogeography, conservation biology, virology and human evolution. The aim of these approaches is to obtain predictive estimates of population history from DNA sequence data that can then be used for hypothesis testing with empirical data. This recent work provides opportunities to evaluate hypotheses of constant population size through time, of population growth or decline, of the rate of growth or decline, and of migration and growth in subdivided populations. At the core of many of these approaches is the extraction of information from the structure of phylogenetic trees to infer the demographic history of a population, and underlying nearly all methods is coalescent theory. With the increasing availability of DNA sequence data, it is important to review the different ways in which information can be extracted from DNA sequence data to estimate demographic parameters.     

Phylogeographic patterns and high levels of chloroplast DNA diversity in four Packera (asteraceae) species in southwestern Alberta

Chloroplast DNA (cpDNA) haplotype variation is compared among alpine and prairie/montane species of Packera from a region in southwestern Alberta that straddles the boundary of Pleistocene glaciation. The phylogeny of the 15 haplotypes identified reveals the presence of two groups: one generally found in coastal and northern species and the other from species in drier habitats. The presence of both groups in all four species and most populations from southwestern Alberta is evidence of past hybridization involving species or lineages that may no longer be present in the region. With the exception of the alpine P. subnuda (phiST = 1.0), interpopulational subdivision of haplotype variation is low (phiST < 0.350), suggesting that interpopulational gene flow is high. However, based on haplotype distribution patterns, we propose that Pleistocene hybridization and incomplete lineage sorting have resulted in reduced subdivision of interpopulational variation so that gene flow may not be as high as indicated. Drift has been more important in the alpine species populations, especially P. subnuda.     

Gene trees and species trees – mutual influences and interdependences of population genetics and systematics

Both population genetics and systematics are core disciplines of evolutionary biology. While systematics deals with genealogical relationships among taxa, population genetics has mainly been based on allele frequencies and the distribution of genetic variants whose genealogical relations could for a long time, due mainly to methodological constraints, not be inferred. The advent of mitochondrial DNA analyses and modern sequencing techniques in the 1970s revolutionized evolutionary genetic studies and gave rise to molecular phylogenetics. In the wake of this new development systematic approaches and principles were incorporated into intraspecific studies at the population level, e.g. the concept of monophyly which is used to delineate evolutionarily significant units in conservation biology. A new discipline combining phylogenetic analyses of genetic lineages with their geographic distribution ('phylogeography') was introduced as an explicit synthesis of population genetics and systematics. On the other hand, it has increasingly become obvious that discordances between gene trees and species trees not only result from spurious data sets or methodological flaws in phylogenetic analyses, but that they often reflect real population genetic processes such as lineage sorting or hybridization. These processes have to be taken into account when evaluating the reliability of gene trees to avoid wrong phylogenetic conclusions. The present review focuses on the phenomenon of non-phylogenetic sorting of ancestral polymorphisms, its probability and its consequences for molecular systematics.     

Barriers to gene flow and ring species formation

Ring species are groups of organisms that dispersed along a ring‐shaped region in such a way that the two ends of the population that meet after many generations are reproductively isolated. They provide a rare opportunity to understand the role of spatial structuring in speciation. Here, we simulate the evolution of ring species assuming that individuals become sexually isolated if the genetic distance between them is above a certain threshold. The model incorporates two forms of dispersal limitation: exogenous geographic barriers that limit the population range and endogenous barriers that result in genetic structuring within the population range. As expected, species' properties that reduce gene flow within the population range facilitate the evolution of reproductive isolation and ring species formation. However, if populations are confined to narrow ranges by geographic barriers, ring species formation increases when local mating is less spatially restricted. Ring species are most likely to form if a population expands while confined to a quasi‐unidimensional range but preserving high mobility in the direction of the range expansion. These conditions are unlikely to be met or persist in real populations and may explain why ring species are rare.     

North American origin and recent European establishments of the amphi-Atlantic peat moss Sphagnum angermanicum

Genetic and morphological similarity between populations separated by large distances may be caused by frequent long-distance dispersal or retained ancestral polymorphism. The frequent lack of differentiation between disjunct conspecific moss populations on different continents has traditionally been explained by the latter model, and has been cited as evidence that many or most moss species are extremely ancient and slowly diverging. We have studied intercontinental differentiation in the amphi-Atlantic peat moss Sphagnum angermanicum using 23 microsatellite markers. Two major genetic clusters are found, both of which occur throughout the distributional range. Patterns of genetic structuring and overall migration patterns suggest that the species probably originated in North America, and seems to have been established twice in Northern Europe during the past 40,000 years. We conclude that similarity between S. angermanicum populations on different continents is not the result of ancient vicariance and subsequent stasis. Rather, the observed pattern can be explained by multiple long-distance dispersal over limited evolutionary time. The genetic similarity can also partly be explained by incomplete lineage sorting, but this appears to be caused by the short time since separation. Our study adds to a growing body of evidence suggesting that Sphagnum, constituting a significant part of northern hemisphere biodiversity, may be more evolutionary dynamic than previously assumed.     

Comparative phylogeography in a genus of coral reef fishes: biogeographic and genetic concordance in the Caribbean

Geographic barriers that limit the movement of individuals between populations may create or maintain phylogenetically discrete lineages. Such barriers are often inferred from geographic surveys of a single mitochondrial marker to identify phylogenetic splits. Mitochondrial DNA, however, has an effective population size one-fourth that of nuclear DNA, which can facilitate the rapid evolution of monophyletic mtDNA lineages in the absence of geographic barriers. The identification of geographic barriers will thus be more robust if barriers are proposed a priori, and tested with multiple independent genetic markers in multiple species. Here, we tested two proposed marine biogeographic breaks located at the Mona Passage in the Caribbean Sea and at the southern end of Exuma Sound in the Bahamas. We sequenced mitochondrial cytochrome b (400 bp) and nuclear rag1 (573 bp) for nine species and colour forms (183 individuals total) within the teleost genus Elacatinus (Gobiidae) that span the proposed breaks. Our results showed that Mona Passage separated mtcyb and rag1 lineages, with no genetic exchange between populations separated by just 23 km. However, the Central Bahamas barrier was only weakly supported by our data. Importantly, neither barrier coincided with deep genetic splits. This suggests that these two barriers did not initially isolate regional populations, but instead disrupt ongoing gene flow between regions. Our inferred relationships further suggested a division of the Caribbean region into northwestern and southeastern regions, a pattern reflected by some freshwater and terrestrial vertebrates. Our results, coupled with genetic and demographic data from other reef fishes and corals, provide robust support for the Mona Passage as a long-term biogeographic barrier for Caribbean animals.     

Isolation with asymmetric gene flow during the nonsynchronous divergence of dry forest birds

Dry forest bird communities in South America are often fragmented by intervening mountains and rainforests, generating high local endemism. The historical assembly of dry forest communities likely results from dynamic processes linked to numerous population histories among codistributed species. Nevertheless, species may diversify in the same way through time if landscape and environmental features, or species ecologies, similarly structure populations. Here we tested whether six co‐distributed taxon pairs that occur in the dry forests of the Tumbes and Marañón Valley of northwestern South America show concordant patterns and modes of diversification. We employed a genome reduction technique, double‐digest restriction site‐associated DNA sequencing, and obtained 4407–7186 genomewide SNPs. We estimated demographic history in each taxon pair and inferred that all pairs had the same best‐fit demographic model: isolation with asymmetric gene flow from the Tumbes into the Marañón Valley, suggesting a common diversification mode. Overall, we also observed congruence in effective population size (N_e) patterns where ancestral N_e were 2.9–11.0× larger than present‐day Marañón Valley populations and 0.3–2.0× larger than Tumbesian populations. Present‐day Marañón Valley N_e was smaller than Tumbes. In contrast, we found simultaneous population isolation due to a single event to be unlikely as taxon pairs diverged over an extended period of time (0.1–2.9 Ma) with multiple nonoverlapping divergence periods. Our results show that even when populations of codistributed species asynchronously diverge, the mode of their differentiation can remain conserved over millions of years. Divergence by allopatric isolation due to barrier formation does not explain the mode of differentiation between these two bird assemblages; rather, migration of individuals occurred before and after geographic isolation.     

Speciational history of Australian grass finches (Poephila) inferred from thirty gene trees

Multilocus genealogical approaches are still uncommon in phylogeography and historical demography, fields which have been dominated by microsatellite markers and mitochondrial DNA, particularly for vertebrates. Using 30 newly developed anonymous nuclear loci, we estimated population divergence times and ancestral population sizes of three closely related species of Australian grass finches (Poephila) distributed across two barriers in northern Australia. We verified that substitution rates were generally constant both among lineages and among loci, and that intralocus recombination was uncommon in our dataset, thereby satisfying two assumptions of our multilocus analysis. The reconstructed gene trees exhibited all three possible tree topologies and displayed considerable variation in coalescent times, yet this information provided the raw data for maximum likelihood and Bayesian estimation of population divergence times and ancestral population sizes. Estimates of these parameters were in close agreement with each other regardless of statistical approach and our Bayesian estimates were robust to prior assumptions. Our results suggest that black-throated finches (Poephila cincta) diverged from long-tailed finches (P. acuticauda and P. hecki) across the Carpentarian Barrier in northeastern Australia around 0.6 million years ago (mya), and that P. acuticauda diverged from P. hecki across the Kimberley Plateau-Arnhem Land Barrier in northwestern Australia approximately 0.3 mya. Bayesian 95% credibility intervals around these estimates strongly support Pleistocene timing for both speciation events, despite the fact that many gene divergences across the Carpentarian region clearly predated the Pleistocene. Estimates of ancestral effective population sizes for the basal ancestor and long-tailed finch ancestor were large (about 521,000 and about 384,000, respectively). Although the errors around the population size parameter estimates are considerable, they are the first for birds taking into account multiple sources of variance.     

Phylogeographic population structure of great reed warblers: an analysis of mtDNA control region sequences

A portion of the mitochondrial control region (494 bp) was sequenced in 106 great reed warblers sampled in six breeding populations in Europe and one wintering population in Africa. In total, 33 different haplotypes were found. There was little evidence of divergence between populations in northern and western Europe whereas the sample from Greece differed significantly from the other European breeding populations. The lowest haplotypc diversity was found near the distribution range limit in Sweden and in The Netherlands suggesting recent effects of bottlenecks/founder events in these areas. A neighbour-joining analysis of the different haplotypes placed the haplotypes into two distinctive clades, A and B. The divergence of the two clades was on average 1.29%. Accounting for the within cladc variation suggested a divergence time between these lines approximately 70 000 years BP. The frequency of the two clades changed longitudinally across Europe with the A haplotypc in the west and the B haplotypc in the east. All birds from Kenya carried the B haplotypc suggesting an origin of these birds east of Latvia/Greece. The long-term female effective population size was estimated to be 20 000 individuals, which is approximately 2% of current population size.     

Phylogeographic concordance factors quantify phylogeographic congruence among co‐distributed species in the Sarracenia alata pitcher plant system

Comparative phylogeographic investigations have identified congruent phylogeographic breaks in co‐distributed species in nearly every region of the world. The qualitative assessments of phylogeographic patterns traditionally used to identify such breaks, however, are limited because they rely on identifying monophyletic groups across species and do not account for coalescent stochasticity. Only long‐standing phylogeographic breaks are likely to be obvious; many species could have had a concerted response to more recent landscape events, yet possess subtle signs of phylogeographic congruence because ancestral polymorphism has not completely sorted. Here, we introduce Phylogeographic Concordance Factors (PCFs), a novel method for quantifying phylogeographic congruence across species. We apply this method to the Sarracenia alata pitcher plant system, a carnivorous plant with a diverse array of commensal organisms. We explore whether a group of ecologically associated arthropods have co‐diversified with the host pitcher plant, and identify if there is a positive correlation between ecological interaction and PCFs. Results demonstrate that multiple arthropods share congruent phylogeographic breaks with S. alata, and provide evidence that the level of ecological association can be used to predict the degree of similarity in the phylogeographic pattern. This study outlines an approach for quantifying phylogeographic congruence, a central concept in biogeographic research.     

Phylogeographic histories of representative herpetofauna of the southwestern U.S.: mitochondrial DNA variation in the desert iguana (Dipsosaurus dorsalis) and the chuckwalla (Sauromalus obesus)

To determine whether genetic variation in representative reptiles of the southwestern U.S. may have been similarly molded by the geologic history of the lower Colorado River, we examined restriction site polymorphisms in the mitochondrial DNA (mtDNA) of desert iguanas (Dipsosaurus dorsalis) and chuckwallas (Sauromalus obesus). Observed phylogeographic structure in these lizards was compared to that reported for the desert tortoise (Xerobates agassizi), whose mtDNA phylogeny demonstrates a striking genetic break at the Colorado River. Both the desert iguana and chuckwalla exhibit extensive mtDNA polymorphism, with respective genotypic diversities G = 0.963 and 0.983, close to the maximum possible value of 1.0. Individual mtDNA clones, as well as clonal assemblages defined by specific levels of genetic divergence, showed pronounced geographic localization. Nonetheless, for each species the distributions of certain clones and most major clonal groupings encompass both sides of the Colorado River valley, and hence are clearly incongruent with the phylogeographic pattern of the desert tortoise. Overall, available molecular evidence provides no indication that the intraspecific phylogenies of the southwestern U.S. herpetofauna have been concordantly shaped by a singular vicariant factor of overriding significance.     

Bayesian estimation of concordance among gene trees

Multigene sequence data have great potential for elucidating important and interesting evolutionary processes, but statistical methods for extracting information from such data remain limited. Although various biological processes may cause different genes to have different genealogical histories (and hence different tree topologies), we also may expect that the number of distinct topologies among a set of genes is relatively small compared with the number of possible topologies. Therefore evidence about the tree topology for one gene should influence our inferences of the tree topology on a different gene, but to what extent? In this paper, we present a new approach for modeling and estimating concordance among a set of gene trees given aligned molecular sequence data. Our approach introduces a one-parameter probability distribution to describe the prior distribution of concordance among gene trees. We describe a novel 2-stage Markov chain Monte Carlo (MCMC) method that first obtains independent Bayesian posterior probability distributions for individual genes using standard methods. These posterior distributions are then used as input for a second MCMC procedure that estimates a posterior distribution of gene-to-tree maps (GTMs). The posterior distribution of GTMs can then be summarized to provide revised posterior probability distributions for each gene (taking account of concordance) and to allow estimation of the proportion of the sampled genes for which any given clade is true (the sample-wide concordance factor). Further, under the assumption that the sampled genes are drawn randomly from a genome of known size, we show how one can obtain an estimate, with credibility intervals, on the proportion of the entire genome for which a clade is true (the genome-wide concordance factor). We demonstrate the method on a set of 106 genes from 8 yeast species.     

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.     

Phylogeographic patterns of two tiger beetle species at both sides of the strait of Gibraltar (Coleoptera: Cicindelini)

The Strait of Gibraltar is one of the major barriers to gene flow between land masses of Europe and Africa at the western end of the Mediterranean. Since the opening of the Strait at the end of the Miocene 5.33 million years ago (Mya) it has exerted a strong influence on the dispersal and colonization of whole biotas, particularly with regard to the retreat of organisms during glaciation peaks and the northward colonization events during the warm periods of Pleistocene. The aim of this study is to elucidate the influence of the Strait of Gibraltar in the gene flow among populations of two tiger beetle species collected in Morocco, Cicindela campestris (L. 1758) and Lophyra flexuosa (Fabricius 1787), with regard to both new and published data from populations of southern Iberia. The phylogeographic analysis showed that Moroccan haplotypes of L. flexuosa belonged to a single coastal mitochondrial clade and that populations at both sides of` the Strait of Gibraltar were genetically well connected. The haplotype network of C. campestris showed that Moroccan populations made up a robust cluster clearly differentiated from those of Iberian and other European populations. These differences are thought to reflect the distinct evolutionary history (dispersal capacity, ecological strategies) of both species, as L. flexuosa shows an almost continuous distribution on the coasts located at both sides of the West Mediterranean, whereas C. campestris shows a patchy distribution and prefers montane habitats in the Western Palaearctic.     

Calibrating a molecular clock from phylogeographic data: moments and likelihood estimators

Abstract We present moments and likelihood methods that estimate a DNA substitution rate from a group of closely related sister species pairs separated at an assumed time, and we test these methods with simulations. The methods also estimate ancestral population size and can test whether there is a significant difference among the ancestral population sizes of the sister species pairs. Estimates presented in the literature often ignore the ancestral coalescent prior to speciation and therefore should be biased upward. The simulations show that both methods yield accurate estimates given sample sizes of five or more species pairs and that better likelihood estimates are obtained if there is no significant difference among ancestral population sizes. The model presented here indicates that the larger than expected variation found in multitaxa datasets can be explained by variation in the ancestral coalescence and the Poisson mutation process. In this context, observed variation can often be accounted for by variation in ancestral population sizes rather than invoking variation in other parameters, such as divergence time or mutation rate. The methods are applied to data from two groups of species pairs (sea urchins and Alpheus snapping shrimp) that are thought to have separated by the rise of Panama three million years ago.     

Phylogeographic studies in plants: problems and prospects

Genetic structuring of plant populations is strongly influenced by both common ancestry and current patterns of interpopulation genetic exchange. The interaction of these two forces is particularly confounding and hence interesting in plants. This complexity of plant genetic structures is due in part to a diversity of reproductive ecologies affecting genetic exchange and the fact that reproductive barriers are often weak between otherwise morphologically well-defined species. Phylogeographic methods provide a means of examining the history of genetic exchange among populations, with the potential to distinguish biogeographic patterns of genetic variation caused by gene flow from those caused by common ancestry. With regard to plants, phylogeography will be most useful when applied broadly across the entire spectrum of potential genetic exchange. Although current phylogeographic studies of plants show promise, widespread application of this approach has been hindered by a lack of appropriate molecular variation; this problem is discussed and possible solutions considered.     

The effect of gene flow on the coalescent time in the human-chimpanzee ancestral population

The coalescent process in the human-chimpanzee ancestral population is investigated using a model, which incorporates a certain time period of gene flow during the speciation process. a is a parameter to represent the degree and time of gene flow, and the model is identical to the null model with an instantaneous species split when a=infinity. A maximum likelihood (ML) method is developed to estimate a, and its power and reliability is investigated by coalescent simulations. The ML method is applied to nucleotide divergence data between human and chimpanzee. It is found that the null model with an instantaneous species split explains the data best, and no strong evidence for gene flow is detected. The result is discussed in the view of the mode of speciation. Another ML method is developed to estimate the male-female ratio (alpha) of mutation rate, in which the coalescent process in the ancestral population is taken into account.