Phylogeography of the barbel (Barbus barbus) assessed by mitochondrial DNA variation

Using the phylogeographic framework, we assessed the DNA sequence variation at the mitochondrial cytochrome b gene across the distribution range of the barbel Barbus barbus, a widely distributed European cyprinid. Reciprocal monophyly of non-Mediterranean European and Balkan/Anatolian populations is taken as evidence for a long-term barrier to gene flow, and interpreted as a consequence of survival of the species in two separate refugia during several later glacial cycles. Lack of profound genealogical divergence across Europe from western France to the northwestern Black Sea basin is consistent with recent colonization of this area from a single glacial refuge, which was probably located in the Danube River basin. This may have occurred in two steps: into the Western European river basins during the last interglacial, and throughout the Central European river basins after the last glacial. The populations from the Balkans and Anatolia apparently did not contribute mitochondrial DNA to the post-Pleistocene colonization of non-Mediterranean Europe. Lack of detectable variation within the Balkans/Anatolia is attributed mainly to recent expansion throughout these regions, facilitated by the freshwater conditions and seashore regression in the Black Sea during the Late Pleistocene and Early Holocene.     

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

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

Separating historic events from recurrent processes in cryptic species: phylogeography of mud snails (Hydrobia spp.)

The present study combines methods that were designed to infer intraspecific relationships (e.g. nested-clade analysis (NCA), mismatch distributions and maximum likelihood gene flow analysis) to analyse historic events and recurrent processes in the cryptic mud snail species Hydrobia acuta and H. glyca. Specifically, we test the proposed allopatry of cryptic species and whether the peculiar range-subdivision of the putative subspecies H. a. acuta and H. a. neglecta is a result of long-distance dispersal or continuous range expansion. The NCA indicates a past fragmentation of the two H. acuta subspecies as well as past fragmentations within H. glyca. Gene-flow analyses show extensive gene flow in an E-W direction (towards the Atlantic) in the Mediterranean H. a. acuta, generally low gene flow in a W-E direction in the Atlantic H. a. neglecta and complex gene-flow pattern in a N-S but also in a S-N direction (against the Gulf Stream) in H. glyca. Based on these data and supportive ecological and oceanographical data, we hypothesize that the separation of the two H. acuta subspecies was not caused by long-distance dispersal but by a range shift and/or range expansion of the closely related competitor H. glyca as a result of an interglacial warming with a subsequent range shift in H. acuta. Moreover, our data do not show evidence for a long-term, stable sympatry of Hydrobia species, supporting the concept of allopatric relationships within cryptic radiations. NCA and gene-flow analyses indicate that the only sympatric population found in our study is the result of a recent dispersal event from the nearby Mediterranean. It is assumed that allopatric relationships in ephemeral Hydrobia populations constitute an evolutionary advantage relative to competition, recruitment and re-establishment of habitats. Mechanisms that could be of relevance for maintaining allopatry are discussed.     

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.     

Phylogeography and genetic structure of northern populations of the yellow warbler (Dendroica petechia)

Phylogeographic patterns of intraspecific variation can provide insights into the population-level processes responsible for speciation and yield information useful for conservation purposes. To examine phylogeography and population structure in a migratory passerine bird at both continental and regional geographical scales, we analysed 344 bp of mitochondrial DNA (mtDNA) control region sequence from 155 yellow warblers (Dendroica petechia) collected from seven locations across Canada and from Alaska. There is a major subdivision between eastern (Manitoba to Newfoundland) and western (Alaska and British Columbia) populations which appears to have developed during the recent Pleistocene. Some localities within these two regions also differ significantly in their genetic composition, suggesting further subdivision on a regional geographical scale. Eastern and western birds form distinct phylogeographic entities and the clustering of all western haplotypes with two eastern haplotypes suggests that the western haplotypes may be derived from an eastern lineage. Analyses based on coalescent models support this explanation for the origin of western haplotypes. These results are consistent with important features of Mengel's model of warbler diversification. From a conservation perspective they also suggest that individual populations of migrant birds may form demographically isolated management units on a smaller scale than previously appreciated.     

A genealogical approach to quantifying lineage divergence

We introduce a statistic, the genealogical sorting index (gsi), for quantifying the degree of exclusive ancestry of labeled groups on a rooted genealogy and demonstrate its application. The statistic is simple, intuitive, and easily calculated. It has a normalized range to facilitate comparisons among different groups, trees, or studies and it provides information on individual groups rather than a composite measure for all groups. It naturally handles polytomies and accommodates measures of uncertainty in phylogenetic relationships. We use coalescent simulations to explore the behavior of the gsi across a range of divergence times, with the mean value increasing to 1, the maximum value when exclusivity within a group reached monophyly. Simulations also demonstrate that the power to reject the null hypothesis of mixed genealogical ancestry increased markedly as sample size increased, and that the gsi provides a statistically more powerful measure of divergence than FST. Applications to data from published studies demonstrated that the gsi provides a useful way to detect significant exclusivity even when groups are not monophyletic. Although we describe this statistic in the context of divergence, it is more broadly applicable to quantify and assess the significance of clustering of observations in labeled groups on any tree.     

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.     

Phylogeography and speciation in the <Emphasis Type="Italic">Pseudocrenilabrus philander</Emphasis> species complex in Zambian Rivers

Haplochromine cichlids form the most species-rich lineage of cichlid fishes that both colonized almost all river systems in Africa and radiated to species flocks in several East African lakes. The enormous diversity of lakes is contrasted by a relatively poor albeit biogeographically clearly structured species diversity in rivers. The present study analyzed the genetic structure and phylogeographic history of species and populations of the genus Pseudocrenilabrus in Zambian rivers that span two major African drainage systems, the Congo- and the Zambezi-system. The mtDNA phylogeny identifies four major lineages, three of which occur in the Congo-system and one in the Zambezi system. Two of the Congo-clades (Lake Mweru and Lunzua River) comprise distinct albeit yet undescribed species, while the fish of the third Congo-drainage clade (Chambeshi River and Bangweulu swamps), together with the fish of the Zambezi clade (Zambezi and Kafue River) are assigned to Pseudocrenilabrus philander. Concerning the intraspecific genetic diversity observed in the sampled rivers, most populations are highly uniform in comparison to lacustrine haplochromines, suggesting severe founder effects and/or bottlenecking during their history. Two bursts of diversification are reflected in the structure of the linearized tree. The first locates at about 3.9% mean sequence divergence and points to an almost simultaneous colonization of the sampled river systems. Subsequent regional diversification (with about 1% mean sequence divergence) occurred contemporaneously within the Kafue River and the Zambezi River. The clear-cut genetic biogeographic structure points to the dominance of geographic speciation in this lineage of riverine cichlid fishes, contrasting the importance of in situ diversification observed in lake cichlids.     

Circum Black Sea phylogeography of Barbus freshwater fishes: divergence in the Pontic glacial refugium

We used DNA sequence variation at the mitochondrial cytochrome b gene (1141 bp) to assess the phylogeography of Barbus fishes in the Black Sea region. Our aim was to test whether the recent ( approximately 22 000-7500 years ago) freshwater phase of the Black Sea was a conduit for gene flow among freshwater fishes that today are found in streams entering the saltwater Black Sea. Deep phylogeographical breaks suggestive of allopatric divergence were observed between four regional groups of populations. Coalescent simulation used to distinguish between this and an alternative scenario that the phylogeographical structure was due to random lineage sorting showed that the contemporary populations were unlikely (P < 0.001) to have been founded by a single ancestral population. Divergences between the lineages (0.86-2.54%) were dated to the Middle to Late Pleistocene using distances and a molecular clock corrected for superimposed substitutions. Taken together, this evidence suggests that multiple refugial populations survived over several later glaciations in the vicinity of the Black Sea. This Pontic refugium served as the primary source for the postglacial expansion throughout Europe as far as the Atlantic basin. However, only one of the phylogeographical lineages contributed to this dispersion, whereas the others remain restricted to the Black Sea region and followed independent evolutionary trajectories.     

Estimating the number of founder lineages from haplotypes of closely linked SNPs

We consider an isolated population founded by a small number of individuals randomly chosen from a source population of known genetic composition at a known time in the past. We develop a Monte-Carlo maximum-likelihood method for estimating the number of founding individuals from the haplotype frequencies at several SNP (single nucleotide polymorphism) loci in a sample. We assume the isolated population was founded recently enough that that mutation can be ignored and that haplotype frequencies in the source population have not changed. We apply the method to simulated data and show that it is unbiased. With a reasonable number of individuals sampled, it is possible to estimate the number of founders within a factor of 2. We show that the performance of the method is not degraded substantially if the frequencies of the rare haplotypes in the source are not known precisely and if there is some recombination. We illustrate the use of our method by applying it to a previously published data set from a recently founded population of wolves (Canis lupus) in Scandinavia.     

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.     

Evolutionary history of bulldog bats (genus Noctilio): recent diversification and the role of the Caribbean in Neotropical biogeography

We present phylogenetic relationships and phylogeographic patterns of the two species of bulldog bats, genus Noctilio. Using a comprehensive sampling of 118 individuals throughout the species distribution, we investigated the distribution of molecular variation in one nuclear and two mitochondrial markers. Phylogenetic trees do not recover Noctilio albiventris as a monophyletic group and point to three similar‐age intraspecific genetic lineages, suggesting cryptic diversity in this taxon. These lineages correspond to the subspecies previously proposed, and are strongly associated with major river basins in South America. Analyses also suggest a very recent origin for the fishing bat Noctilio leporinus, which probably originates from N. albiventris, with a population expansion corresponding to its invasion in South America. Based on our analysis, the speciation event was dated in the Pleistocene epoch and seems to be associated with the variation of the sea level in the Caribbean islands. The present work indicates how phylogeographic studies support the identification of independent evolutionary lineages, driving new systematic/taxonomic investigations, while at the same time shed light on the role of the Caribbean in shaping Neotropical bat fauna diversity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

PHYLOGEOGRAPHY OF SPOTTED OWL (STRIX OCCIDENTALIS) POPULATIONS BASED ON MITOCHONDRIAL DNA SEQUENCES: GENE FLOW,GENETIC STRUCTURE,AND A NOVEL BIOGEOGRAPHIC PATTERN

Mitochondrial DNA control region sequences of spotted owls (Strix occidentalis) allowed us to investigate gene flow, genetic structure, and biogeographic relationships among these forest-dwelling birds of western North America Estimates of gene flow based on genetic partitioning and the phylogeography of haplotypes indicate substantial dispersal within three long-recognized subspecies. However, patterns of individual phyletic relationships indicate a historical absence of gene flow among the subspecies, which are essentially monophyletic. The pattern of haplotype coalescence enabled us to identify the approximate timing and direction of a recent episode of gene flow from the Sierra Nevada to the northern coastal ranges. The three subspecies comprise phylogenetic species, and the northern spotted owl (S. o. caurina) is sister to a clade of California (S. o. occidentalis) plus Mexican spotted owls (S o lucida); this represents a novel biogeographic pattern within birds. The California spotted owl had substantially lower nucleotide diversity than the other two subspecies; this result is inconsistent with present patterns of population density A causal explanation requires postulating a severe bottleneck or a selective sweep, either of which was confined to only one geographic region.     

Mitochondrial genealogy of Western Australian barramundi: applications of inbreeding coefficients and coalescent analysis for separating temporal population processes

In Australian populations of barramundi Lates calcarifer , phylogenetic reconstruction of mtDNA sequences provided evidence of significant historical levels of gene flow, despite the substantial structuring of contemporary populations. The geographical pattern of mtDNA sequences among the populations was not congruent with previous evidence of a major disjunction between western and eastern populations of barramundi in Australia.     

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.     

Cryptic variation in mole voles Ellobius (Arvicolinae,Rodentia) of Mongolia

The mole vole subgenus Ellobius is currently considered to include three species: Ellobius talpinus (distributed from SE Europe and Turkmenistan through Kazakhstan to SW Siberia), Ellobius alaicus (S Tianshan, Pamir-Alay) and Ellobius tancrei (East and West Central Asia, from the Amu-Darya to Mongolia and N China). A study focusing on the genetic variation in Ellobius from Mongolia was conducted using one mitochondrial and three nuclear markers. Two divergent allopatric lineages endemic to East Central Asia were revealed. The first lineage occurs from Dzungaria eastwards to central Mongolia and represents E. tancrei sensu stricto. The second lineage is found in East Gobi only and corresponds to a taxon described as Ellobius orientalis, which has been traditionally treated as a subspecies of E. tancrei. However, molecular and chromosome data indicate that orientalis is related not to E. tancrei but to E. talpinus, which is separated from the former by a distribution gap of ~2,000 km. The taxonomic status of the East Gobi mole vole is ambiguous, and its genetic distance from E. talpinus s. str. falls into the range characteristic for closely related vole species or semi-species. According to molecular estimates, the two taxa have been isolated since the late Middle Pleistocene. A similar divergence is observed between the East and West Central Asian lineages of E. tancrei. E. alaicus is placed as sister to the latter rendering E. tancrei sensu lato paraphyletic. The revealed phylogeographic pattern implies that East Central Asia was colonized by mole voles through multiple eastward dispersal events.     

Genetic variation in gorillas

This review summarizes what is currently known concerning genetic variation in gorillas, on both inter- and intraspecific levels. Compared to the human species, gorillas, along with the other great apes, possess greater genetic variation as a consequence of a demographic history of rather constant population size. Data and hence conclusions from analysis of mitochondrial DNA (mtDNA), the usual means of describing intraspecific patterns of genetic diversity, are limited at this time. An important task for future studies is to determine the degree of confidence with which gorilla mtDNA can be analyzed, in view of the risk that one will inadvertently analyze artifactual rather than genuine sequences. The limited information available from sequences of nuclear genomic segments does not distinguish western from eastern gorillas, and, in comparison with results from the two chimpanzee species, suggests a relatively recent common ancestry for all gorillas. In the near future, the greatest insights are likely to come from studies aimed at genetic characterization of all individual members of social groups. Such studies, addressing topics such as behavior of individuals with kin and non-kin, and the actual success of male reproductive strategies, will provide a link between behavioral and genetic studies of gorillas.     

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.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

Statistical phylogeography

While studies of phylogeography and speciation in the past have largely focused on the documentation or detection of significant patterns of population genetic structure, the emerging field of statistical phylogeography aims to infer the history and processes underlying that structure, and to provide objective, rather than ad hoc explanations. Methods for parameter estimation are now commonly used to make inferences about demographic past. Although these approaches are well developed statistically, they typically pay little attention to geographical history. In contrast, methods that seek to reconstruct phylogeographic history are able to consider many alternative geographical scenarios, but are primarily nonstatistical, making inferences about particular biological processes without explicit reference to stochastically derived expectations. We advocate the merging of these two traditions so that statistical phylogeographic methods can provide an accurate representation of the past, consider a diverse array of processes, and yet yield a statistical estimate of that history. We discuss various conceptual issues associated with statistical phylogeographic inferences, considering especially the stochasticity of population genetic processes and assessing the confidence of phylogeographic conclusions. To this end, we present some empirical examples that utilize a statistical phylogeographic approach, and then by contrasting results from a coalescent-based approach to those from Templeton's nested cladistic analysis (NCA), we illustrate the importance of assessing error. Because NCA does not assess error in its inferences about historical processes or contemporary gene flow, we performed a small-scale study using simulated data to examine how our conclusions might be affected by such unconsidered errors. NCA did not identify the processes used to simulate the data, confusing among deterministic processes and the stochastic sorting of gene lineages. There is as yet insufficient justification of NCA's ability to accurately infer or distinguish among alternative processes. We close with a discussion of some unresolved problems of current statistical phylogeographic methods to propose areas in need of future development.