Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies

Molecular methods as applied to the biogeography of single species (phylogeography) or multiple codistributed species (comparative phylogeography) have been productively and extensively used to elucidate common historical features in the diversification of the Earth's biota. However, only recently have methods for estimating population divergence times or their confidence limits while taking into account the critical effects of genetic polymorphism in ancestral species become available, and earlier methods for doing so are underutilized. We review models that address the crucial distinction between the gene divergence, the parameter that is typically recovered in molecular phylogeographic studies, and the population divergence, which is in most cases the parameter of interest and will almost always postdate the gene divergence. Assuming that population sizes of ancestral species are distributed similarly to those of extant species, we show that phylogeographic studies in vertebrates suggest that divergence of alleles in ancestral species can comprise from less than 10% to over 50% of the total divergence between sister species, suggesting that the problem of ancestral polymorphism in dating population divergence can be substantial. The variance in the number of substitutions (among loci for a given species or among species for a given gene) resulting from the stochastic nature of DNA change is generally smaller than the variance due to substitutions along allelic lines whose coalescence times vary due to genetic drift in the ancestral population. Whereas the former variance can be reduced by further DNA sequencing at a single locus, the latter cannot. Contrary to phylogeographic intuition, dating population divergence times when allelic lines have achieved reciprocal monophyly is in some ways more challenging than when allelic lines have not achieved monophyly, because in the former case critical data on ancestral population size provided by residual ancestral polymorphism is lost. In the former case differences in coalescence time between species pairs can in principle be explained entirely by differences in ancestral population size without resorting to explanations involving differences in divergence time. Furthermore, the confidence limits on population divergence times are severely underestimated when those for number of substitutions per site in the DNA sequences examined are used as a proxy. This uncertainty highlights the importance of multilocus data in estimating population divergence times; multilocus data can in principle distinguish differences in coalescence time (T) resulting from differences in population divergence time and differences in T due to differences in ancestral population sizes and will reduce the confidence limits on the estimates. We analyze the contribution of ancestral population size (theta) to T and the effect of uncertainty in theta on estimates of population divergence (tau) for single loci under reciprocal monophyly using a simple Bayesian extension of Takahata and Satta's and Yang's recent coalescent methods. The confidence limits on tau decrease when the range over which ancestral population size theta is assumed to be distributed decreases and when tau increases; they generally exclude zero when tau/(4Ne) > 1. We also apply a maximum-likelihood method to several single and multilocus data sets. With multilocus data, the criterion for excluding tau = 0 is roughly that l tau/(4Ne) > 1, where l is the number of loci. Our analyses corroborate recent suggestions that increasing the number of loci is critical to decreasing the uncertainty in estimates of population divergence time.     

Lizards as model organisms for linking phylogeographic and speciation studies

Lizards have been model organisms for ecological and evolutionary studies from individual to community levels at multiple spatial and temporal scales. Here we highlight lizards as models for phylogeographic studies, review the published population genetics/phylogeography literature to summarize general patterns and trends and describe some studies that have contributed to conceptual advances. Our review includes 426 references and 452 case studies: this literature reflects a general trend of exponential growth associated with the theoretical and empirical expansions of the discipline. We describe recent lizard studies that have contributed to advances in understanding of several aspects of phylogeography, emphasize some linkages between phylogeography and speciation and suggest ways to expand phylogeographic studies to test alternative pattern‐based modes of speciation. Allopatric speciation patterns can be tested by phylogeographic approaches if these are designed to discriminate among four alternatives based on the role of selection in driving divergence between populations, including: (i) passive divergence by genetic drift; (ii) adaptive divergence by natural selection (niche conservatism or ecological speciation); and (iii) socially‐mediated speciation. Here we propose an expanded approach to compare patterns of variation in phylogeographic data sets that, when coupled with morphological and environmental data, can be used to to discriminate among these alternative speciation patterns. [Correction made after online publication (28/07/2010): (minor deletion in the last line of the abstract)].     

Environmental DNA phylogeography: Successful reconstruction of phylogeographic patterns of multiple fish species from cups of water

Phylogeography is an integrative field of science linking micro- and macro-evolutionary processes, contributing to the inference of vicariance, dispersal, speciation, and other population-level processes. Phylogeographic surveys usually require considerable effort and time to obtain numerous samples from many geographical sites covering the distribution range of target species; this associated high cost limits their application. Recently, environmental DNA (eDNA) analysis has been useful not only for detecting species but also for assessing genetic diversity; hence, there has been growing interest in its application to phylogeography. As the first step of eDNA-based phylogeography, we examined (1) data screening procedures suitable for phylogeography and (2) whether the results obtained from eDNA analysis accurately reflect known phylogeographic patterns. For these purposes, we performed quantitative eDNA metabarcoding using group-specific primer sets in five freshwater fish species belonging to two taxonomic groups from a total of 94 water samples collected from western Japan. As a result, three-step data screening based on the DNA copy number of each haplotype detected successfully eliminated suspected false positive haplotypes. Furthermore, eDNA analysis could almost perfectly reconstruct the phylogenetic and phylogeographic patterns obtained for all target species with the conventional method. Despite existing limitations and future challenges, eDNA-based phylogeography can significantly reduce survey time and effort and is applicable for simultaneous analysis of multiple species in single water samples. eDNA-based phylogeography has the potential to revolutionize phylogeography.     

Testing for shared biogeographic history in the lower Central American freshwater fish assemblage using comparative phylogeography: concerted,independent, or multiple evolutionary responses?

A central goal of comparative phylogeography is determining whether codistributed species experienced (1) concerted evolutionary responses to past geological and climatic events, indicated by congruent spatial and temporal patterns (“concerted‐response hypothesis”); (2) independent responses, indicated by spatial incongruence (“independent‐response hypothesis”); or (3) multiple responses (“multiple‐response hypothesis”), indicated by spatial congruence but temporal incongruence (“pseudocongruence”) or spatial and temporal incongruence (“pseudoincongruence”). We tested these competing hypotheses using DNA sequence data from three livebearing fish species codistributed in the Nicaraguan depression of Central America (Alfaro cultratus, Poecilia gillii, and Xenophallus umbratilis) that we predicted might display congruent responses due to co‐occurrence in identical freshwater drainages. Spatial analyses recovered different subdivisions of genetic structure for each species, despite shared finer‐scale breaks in northwestern Costa Rica (also supported by phylogenetic results). Isolation‐with‐migration models estimated incongruent timelines of among‐region divergences, with A. cultratus and Xenophallus populations diverging over Miocene–mid‐Pleistocene while P. gillii populations diverged over mid‐late Pleistocene. Approximate Bayesian computation also lent substantial support to multiple discrete divergences over a model of simultaneous divergence across shared spatial breaks (e.g., Bayes factor [B₁₀] = 4.303 for Ψ [no. of divergences] > 1 vs. Ψ = 1). Thus, the data support phylogeographic pseudoincongruence consistent with the multiple‐response hypothesis. Model comparisons also indicated incongruence in historical demography, for example, support for intraspecific late Pleistocene population growth was unique to P. gillii, despite evidence for finer‐scale population expansions in the other taxa. Empirical tests for phylogeographic congruence indicate that multiple evolutionary responses to historical events have shaped the population structure of freshwater species codistributed within the complex landscapes in/around the Nicaraguan depression. Recent community assembly through different routes (i.e., different past distributions or colonization routes), and intrinsic ecological differences among species, has likely contributed to the unique phylogeographical patterns displayed by these Neotropical fishes.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

INTRASPECIFIC PHYLOGEOGRAPHY ACROSS THE POINT CONCEPTION BIOGEOGRAPHIC BOUNDARY

Recent studies of intraspecific phylogeography have suggested that the geographic location of genetic discontinuities, or phylogeographic breaks, may frequently coincide with biogeographic boundaries. The concordance is hypothesized to reflect similarity in the processes governing species boundaries and intraspecific lineage boundaries. This concordance has not, however, been widely tested. In the case of the Point Conception biogeographic boundary between the Oregonian and Californian marine biotas, only the supralittoral copepod Tigriopus californicus has been found to have a coincident phylogeographic break. Here I show that the apparent relationship between this break and Point Conception was, in fact, an artifact of insufficient geographic sampling. Mitochondrial DNA analyses of T. californicus populations between Morro Bay and San Diego reveal at least five equally deep phylogeographic breaks in the region (where only one biogeographic boundary is recognized). Limited nuclear DNA sequence data and allozyme data also support the occurrence of multiple genetic discontinuities along this geographic range. Lack of one-to-one correspondence between intraspecific phylogeography and biogeographic boundaries indicates that the processes affecting the genetic differentiation of populations of T. californicus differ from those responsible for determining species distributional limits at the Point Conception biogeographic boundary. A review of genetic data from other species also fails to provide evidence for concordance of biogeography and intraspecific phylogeography across Point Conception. I suggest that the concordance of phylogeography with biogeography will only be pronounced where the biogeographic boundary separates biotas that are phylogenetically related. The numerous cases of interspecific hybrid zones in the region of Cape Canaveral, for example, indicate that many sister-species pairs occur across this biogeographic boundary. Such hybrid zones are not common at Point Conception, and there appears to be no cases of intraspecific phylogeographic breaks associated with this well-recognized biogeographic boundary.     

Pleistocene divergence of Dinaric Drusus endemics (Trichoptera, Limnephilidae) in multiple microrefugia within the Balkan Peninsula

The Balkan Peninsula is one of three major European refugial areas. It has high biodiversity and endemism, but data on the age and origin of its fauna, especially endemics, are limited. Mitochondrial sequence data (COI and 16S genes) were used to study the population structure and phylogeography of the caddisfly Drusus croaticus and the phylogeny and divergence of seven other Drusus species, mostly range-restricted endemics of the Dinaric region of the Balkan Peninsula. The divergence of D. croaticus populations in Croatia and allopatric Drusus species in Bosnia dated to the Pleistocene, showing the importance of this time period for the origin and diversification of Balkan endemic taxa. The divergence of more distantly related species dated to the Late Miocene/Early Pliocene. Population genetic and phylogeographic analysis of 115 individuals from 11 populations of D. croaticus revealed a high level of genetic differentiation and absence of gene flow between populations separated by more than 10 km. The existence of allopatrically fragmented lineages in D. croaticus and the endemic Bosnian species is most likely the result of long-term isolation in multiple microrefugia, probably due to the specific habitat requirements and life-history traits of Drusinae coupled with the topographic complexity and historical changes in geomorphology of the region. Overall, these findings shed light on the processes generating the high genetic complexity of this refugial region that parallels the 'refugia within refugia' pattern widely reported from the Iberian refugium.     

Is phylogeography helpful for invasive species risk assessment? The case study of the bark beetle genus Dendroctonus

Despite evidence that conspecific lineages may display different climatic tolerances, most invasion risk assessment tools are calibrated without considering phylogeographic information. This study aims to investigate the existence of intraspecific niche divergence within a group of insect pests and to explore how the inclusion of phylogeographic information into species distribution models may alter the estimation of the potential distribution of a species. We studied North American bark beetles belonging to the genus Dendroctonus, a group of pests of conifers that are listed as quarantine species in numerous countries. Most Dendroctonus species exhibit strong genetic divergence that appears to be geographically structured and shaped by historical events and biotic factors. We modeled all lineage distributions within five species, using MaxEnt and Boosted Regression Trees, and compared the results with the models fitted at the species scale. Multivariate analysis and niche similarity and equivalency tests were additionally performed to investigate the existence and magnitude of climatic niche divergence between conspecific lineages. We also tested the ability of lineage‐based models to predict the region invaded by D. valens in China. Conspecific lineages showed a climatic niche more similar than expected by chance, but displayed different climatic envelopes in their native range and, consequently, different estimates of potential distributions. We also observed that classical models calibrated using the entire range of the species could potentially under‐ or overestimate the potential range of the species when compared to a global prediction built by aggregating lineage‐based projections. This study showed that the invasive phylogeographic lineage of D. valens has invaded regions characterized by climatic conditions highly similar to those encountered in its native range suggesting that preadaptations to environment might have played a role in this invasion. This study highlights how our perception of the invasion risk of pests may be altered when integrating phylogeographic information.     

Phylogeographic evidence for two mesic refugia in a biodiversity hotspot

Phylogeographic studies of flora in species-rich south-western Australia point to complex evolutionary histories, reflecting patterns of persistence and resilience to climatic changes during the Pleistocene. We asked whether coastal areas of the mid-west and south, as well as granite outcrops and inland ranges, have acted as major refugia within this region during Pleistocene climatic fluctuations by analysing phylogeographic patterns in the shrub Calothamnus quadrifidus R.Br. (Myrtaceae). We determined variation in chloroplast DNA data for 41 populations across the geographic range. Relationships and major clades were resolved using parsimony and Bayesian analyses. We tested for demographic and spatial expansion of the major clades and estimated clade divergence dates using an uncorrelated, lognormal relaxed clock based on two conservative chloroplast mutation rates. Two distinct phylogeographic clades were identified showing divergence during the Pleistocene, consistent with other phylogeographic studies of south-west Australian flora, emphasising the impact of climatic oscillations in driving divergence in this landscape. The southern clade was more diverse, having higher haplotype diversity and greater genetic structure, while the northern clade showed evidence of fluctuation in population size. Regions of high haplotype diversity with adjacent areas of low diversity observed in each clade indicated the locations of two coastal refugia: one on the south coast and another along the mid-west coast. This is the first evidence for major Pleistocene refugia using chloroplast genetic data in a common, widespread species from this region.     

Integrating coalescent and ecological niche modeling in comparative phylogeography

Understanding the factors that contribute to the formation of population genetic structure is a central goal of phylogeographic research, but achieving this goal can be complicated by the stochastic variance inherent to genetic processes. Statistical approaches to testing phylogeographic hypotheses accommodate this stochasticity by evaluating competing models of putative historical population structure, often by simulating null distributions of the expected variance. The effectiveness of these tests depends on the biological realism of the models. Information from the fossil record can aid in reconstructing the historical distributions of some taxa. However, for the majority of taxa, which lack sufficient fossils, paleodistributional modeling can provide valuable spatial-geographic data concerning ancestral distributions. Paleodistributional models are generated by projecting ecological niche models, which predict the current distribution of each species, onto a model of past climatic conditions. Here, we generate paleodistributional models describing the suitable habitat during the last glacial maximum for lineages from the mesic forests of the Pacific Northwest of North America, and use these models to generate alternative phylogeographic hypotheses. Coalescent simulations are then used to test these hypotheses to improve our understanding of the historical events that promoted the formation of population genetic structure in this ecosystem. Results from Pacific Northwest mesic forest organisms demonstrate the utility of these combined approaches. Paleodistribution models and population genetic structure are congruent across three amphibian lineages, suggesting that they have responded in a concerted manner to environmental change. Two other species, a willow and a water vole, despite being currently codistributed and having similar population genetic structure, were predicted by the paleodistributional model to have had markedly different distributions during the last glacial maximum. This suggests that congruent phylogeographic patterns can arise from incongruent ancestral distributions. Paleodistributional models introduce a much-needed spatial-geographic perspective to statistical phylogeography. In conjunction with coalescent models of population genetic structure, they have the potential to improve our understanding of the factors that promote population divergence and ultimately produce regional patterns of biodiversity.     

Application of non-coding DNA regions in intraspecific analyses

In this review we discuss the use of non-coding DNA at the intraspecific level in plants. Both nuclear and organelle non-coding regions are widely used in interspecific phylogenetic approaches. However, they are also valuable in analyses on the intraspecific level. Besides taxonomy, that is, defining subspecies or varieties, large fields for the application of non-coding DNA are population genetic and phylogeographic studies. Population genetics tries to explain the genetic patterns within species mostly by the amount of extant gene flow among populations, while phylogeography explicitly tries to reconstruct historic events. Depending on the study different molecular markers can be used, varying between very fast evolving microsatellites or some more slowly changing regions like intergenic spacers and introns. Here, we focus mainly on the use of non-coding regions in phylogeographic analyses. Mostly used in this context are regions of the genomes of the chloroplasts and mitochondria. In phylogeography, the correct estimation of allele or haplotype relationships is particularly important. As tree-based methods are mostly insufficient to depict relationships within species, network approaches are better suitable to infer gene or locus genealogies. Problematic for phylogeographic studies are alleles shared among multiple species, which could result from either hybridization or incomplete lineage sorting. Especially the latter can severely influence the interpretation of the phylogeographic patterns. Therefore, it seems necessary for us to also include close relatives of the species under study in phylogeographic analyses. Not only the sample design but also the analysis methods are currently changing, as some new methods such as statistical phylogeography were emerging recently and widely used methods like nested clade analysis might not be reliable in every case. During the last few years, a multitude of studies were published, which mainly analyzed phylogeographic patterns in European and North American plants. Phylogeographic studies in other regions of the earth are still comparably rare, although questions like the influence of the ice age on the vegetation in the tropics or southern hemisphere are still open and phylogeography provides an excellent remedy to answer them.     

Phylogeographic patterns in widespread corvid birds

Intraspecific genetic diversity and phylogeography of Corvus corone was investigated using the mitochondrial (mt) control region as a molecular marker. A split into two distinct mt lineages was observed. One represents individuals from a wide geographic range spanning from England to the Russian Far East (Kamchatka), while the other one was found in the Primorye and Khabarovsk regions (southern parts of Russian Far East) as well as Japan. For comparison, we investigated several widespread Palearctic corvid taxa with respect to their phylogeographic patterns. A deep split into two lineages was revealed in five cases: Besides C. corone, within Corvus frugilegus, Pica pica, and between the species pairs Corvus monedula-Corvus dauuricus and Cyanopica cyanus-Cyanopica cooki. Although these taxa display a variety of distribution patterns, from disjunct, para/allopatric to continuous, the genetic pattern and level of divergence between clades is very similar. This implies that the differentiation started in about the same time range. In contrast, no differentiation into highly divergent lineages was detected in Corvus corax, Perisoreus infaustus, and Nucifraga caryocatactes. We try to explain the two phylogeographic patterns in corvid birds with ecological factors accompanying the changing climatic conditions during the Pleistocene. The deep genetic splits within several widely distributed Palearctic corvids are discussed with respect to taxonomic questions.     

Fine-scale phylogeographic congruence despite demographic incongruence in two low-mobility saproxylic springtails

Evolutionary trajectories of codistributed taxa with comparable ecological preferences and dispersal abilities may be similarly impacted by historical landscape-level processes. Species' responses to changes in a shared biogeographic landscape may be purely concerted, completely independent, or classified as falling within an intermediate part of the continuum bounded by these two extremes. With sufficient molecular data, temporal contrasts of congruence among taxa with respect to these responses can be made. Such contrasts provide insights into the relative influence of ancient versus more recent climatic (and other) impacts on genetic structuring. Using phylogenetic, allele frequency, and genotypic data from two low-mobility, rotting-log-adapted (saproxylic) springtail species (Collembola) from an isolated 100-km-long section of the Great Dividing Range in southeastern Australia, we tested the concerted-response hypothesis over three timescales. Tests of phylogeographic, demographic, and contemporary population-genetic congruence were performed using an integrative approach that draws on both direct (pattern-based) and indirect (scenario-based) analyses. Our data revealed a general pattern of broad-scale similarities in species' responses to the interaction between Pleistocene climatic cycles and landscape setting, overlaid with some species-specific differences on local geographic and more recent temporal scales. This general pattern of phylogeographic congruence was accompanied by evidence for contemporaneous demographic incongruence indicating that, even at relatively small spatial scales, biogeographic context can exert an overarching influence on genetic structuring.     

Significant Asia‐Europe divergence in the middle spotted woodpecker (Aves,Picidae)

Population divergence could be strongly affected by species’ ecology and might not be a direct response to climate‐driven environmental change. We tested this in the middle spotted woodpecker (Dendrocoptes medius), a non‐migratory, low‐dispersal habitat specialist associated with old deciduous forests of the Western Palearctic. We present the first phylogeographic study of this species integrating genetic data (three mitochondrial loci, one autosomal and one Z‐linked intron) with species distribution modelling. Based on this species’ ecology, we predicted that the middle spotted woodpecker could have colonized its current range from multiple Last Glacial Maximum (LGM) refugia and that strongly structured populations could be expected. Indeed, we discovered a strong genetic divergence between Asian and European populations, with a split estimated at around one million of years ago. This was surprising given only slight intraspecific variation in plumage and morphology. Although there was no significant phylogeographic structure within the Asian and European groups, we cannot exclude the possibility of multiple refugia within either group during the LGM. This has to be further investigated with more extensive geographic sampling and larger number of variable independently evolving markers. Future studies should also investigate potential differences in vocalizations and ecology between the two groups. Lineages showing similar level of genetic differentiation including woodpeckers are often treated as species‐level taxa. Comparison of our results with the phylogeographic history of other woodpeckers, suggests that sympatric species with similar life‐histories might have idiosyncratic phylogeographic patterns probably resulting from different ecological requirements or historic stochasticity.     

Phylogeographic breaks and the mechanisms of their formation in the Sino-Japanese floristic region

Due to combinations of diverse geography and climate, and complex geo-climate histories and sea level fluctuations, the Sino-Japanese floristic region has extremely high species diversity. Phylogeography is an effective method to identify the factors triggering the formation and differentiation of species diversity. Previous studies showed that phylogeographic breaks, the genetic discontinuity between different gene genealogies, were ubiquitously present. From the west to the east, seven general phylogeographic breaks occur, including the Mekong- Salween Divide, the Tanaka-Kaiyong Line, the Sichuan Basin, ca. 105° E, the boundary between the Second and Third ladders, the North China, and the East China Sea and Korea Strait. These phylogeographic breaks are mainly attributable to both historical and ecological factors, which are generally due to a combined effect of the isolation by distance (IBD) and the isolation by environment (IBE). Geological events and climate changes are the historical factors, mainly including the uplift of Qinghai-Xizang Plateau, the formation and intensification of the Asian monsoon and the Asian interior aridification, the redevelopment of the arid belt, and the Quaternary climate oscillations and sea level fluctuations. Adaptive divergence, namely the divergence induced by different selective pressures under different environments, is responsible for the ecological factors. Adaptive divergence could obstacle gene flow among populations, resulting in the formation of phylogeographic break. However, an identical phylogeographic break is not shared by all the plants because of their various intrinsic biological characteristics, among which the difference in dispersal ability is most important. Finally, we envisaged the future development of phylogeographic break studies based on accurate divergence time estimation, relative contribution of IBD and IBE, and also the utilization of comparative phylogeography.     

Contrasting insights provided by single and multispecies data in a regional comparative phylogeographic study

Many single‐species freshwater phylogeographic studies have been carried out in south‐east Queensland; however comparative phylogeography requires multiple lines of evidence to infer deep, significant relationships between landscape and biota. The present study aimed to test conclusions resulting from single taxon studies in a multispecies comparative framework: (1) how influential are river basins in the genetic structure of freshwater species; (2) are there biogeographic frontiers between groups of basins; and (3) could deep intraspecific lineages be explained by a single event? New and existing data from 33 freshwater species (23 fishes and 10 crustaceans) were combined, and both standard single‐species analyses (haplotype networks, genetic distances, Φ_ST) and multispecies methods (hierarchical ABC) were carried out for 1814 sequences from eight basins. More than half of the species displayed a high phylogeographic structure and contained at least two distinct lineages. Almost all of the lineage divergences displayed an element of north/south geographic breaks, with the most influential boundary being between the Mary and Brisbane rivers. Of the 11 basin‐pair multispecies coalescent analyses, four implied a single divergence as being most likely. A regional analysis of deep lineages within 16 taxon‐pairs resulted in a strongly supported inference of a single divergence, probably dating to the Pleistocene. Basin boundaries are a key determinant of phylogeographic patterns for most of these freshwater species, although the specific biogeographic relationship between basins often varies depending on the species. There are a number of influential biogeographic frontiers, with the Brisbane‐Mary being the most important. The finding that a single event may be responsible for multiple deep lineages across the region implies that a highly influential climate change event may have been detected. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 554–569.     

Differences in Quaternary co-divergence reveals community-wide diversification in the mountains of southwest China varied among species

The mountains of southwest China (MSWC) is a biodiversity hotspot with highly complex and unusual terrain. However, with the majority of studies focusing on the biogeographic consequences of massive mountain building, the Quaternary legacy of biodiversity for the MSWC has long been overlooked. Here, we took a statistical comparative phylogeography approach to examine factors that shaped community-wide diversification. With data from 30 vertebrate species, the results reveal spatially concordant genetic structure, and temporally clustered co-divergence events associated with river barriers during severe glacial cycles. This indicates the importance of riverine barriers in the phylogeographic history of the MSWC vertebrate community. We conclude that the repeated glacial cycles are associated with co-divergences that are themselves structured by the heterogeneity of the montane landscape of the MSWC. This orderly process of diversification has profound implications for conservation by highlighting the relative independence of different geographical areas in which some, but not all species in communities have responded similarly to climate change and calls for further comparative phylogeographic investigations to reveal the connection between biological traits and divergence pulses in this biodiversity hotspot.     

SuperFine: fast and accurate supertree estimation

Many research groups are estimating trees containing anywhere from a few thousands to hundreds of thousands of species, toward the eventual goal of the estimation of a Tree of Life, containing perhaps as many as several million leaves. These phylogenetic estimations present enormous computational challenges, and current computational methods are likely to fail to run even on data sets in the low end of this range. One approach to estimate a large species tree is to use phylogenetic estimation methods (such as maximum likelihood) on a supermatrix produced by concatenating multiple sequence alignments for a collection of markers; however, the most accurate of these phylogenetic estimation methods are extremely computationally intensive for data sets with more than a few thousand sequences. Supertree methods, which assemble phylogenetic trees from a collection of trees on subsets of the taxa, are important tools for phylogeny estimation where phylogenetic analyses based upon maximum likelihood (ML) are infeasible. In this paper, we introduce SuperFine, a meta-method that utilizes a novel two-step procedure in order to improve the accuracy and scalability of supertree methods. Our study, using both simulated and empirical data, shows that SuperFine-boosted supertree methods produce more accurate trees than standard supertree methods, and run quickly on very large data sets with thousands of sequences. Furthermore, SuperFine-boosted matrix representation with parsimony (MRP, the most well-known supertree method) approaches the accuracy of ML methods on supermatrix data sets under realistic conditions.     

Complex phylogeographic patterns indicate Central American origin of two widespread Mesoamerican <Emphasis Type="Italic">Quercus</Emphasis> (Fagaceae) species

The northern Neotropical region is characterized by a heterogeneous geological and climatic history. Recent studies have shown contrasting patterns regarding the role of geographic elements as barriers that could have determined phylogeographic structure in various species. Recently, the phylogeography and biogeography of Quercus species have been studied intensively, and the patterns observed so far suggest contrasting evolutionary histories for Neotropical species in comparison with their Holarctic relatives. The goal of this study was to describe the phylogeographic structure of two Neotropical oak species (Quercus insignis and Quercus sapotifolia) in the context of the geological and palaeoclimatic history of the northern Neotropics. Populations through the distribution range of both species were collected and characterized using nine chloroplast DNA microsatellite loci. Both oak species showed high levels of genetic diversity and strong phylogeographic structure. The distribution of genetic variation in Q. insignis suggested an influence of two major barriers, the Isthmus of Tehuantepec and the Nicaraguan Depression, while Q. sapotifolia exhibited a genetic structure defined by the heterogeneity of the Chortis highlands. The haplotype networks of both species indicated complex histories, suggesting that colonization from the Sierra Madre de Chiapas to central Mexico and from the north of the Nicaraguan Depression to the Costa Rican mountains may have occurred during different stages, and apparently more than one time. In conclusion, the phylogeographic structure of Neotropical oak species seems to be defined by a combination of geological and climatic events.