Comparative phylogeographic summary statistics for testing simultaneous vicariance

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

EVIDENCE FOR CLIMATE‐DRIVEN DIVERSIFICATION? A CAUTION FOR INTERPRETING ABC INFERENCES OF SIMULTANEOUS HISTORICAL EVENTS

Approximate Bayesian computation (ABC) is rapidly gaining popularity in population genetics. One example, msBayes , infers the distribution of divergence times among pairs of taxa, allowing phylogeographers to test hypotheses about historical causes of diversification in co‐distributed groups of organisms. Using msBayes , we infer the distribution of divergence times among 22 pairs of populations of vertebrates distributed across the Philippine Archipelago. Our objective was to test whether sea‐level oscillations during the Pleistocene caused diversification across the islands. To guide interpretation of our results, we perform a suite of simulation‐based power analyses. Our empirical results strongly support a recent simultaneous divergence event for all 22 taxon pairs, consistent with the prediction of the Pleistocene‐driven diversification hypothesis. However, our empirical estimates are sensitive to changes in prior distributions, and our simulations reveal low power of the method to detect random variation in divergence times and bias toward supporting clustered divergences. Our results demonstrate that analyses exploring power and prior sensitivity should accompany ABC model selection inferences. The problems we identify are potentially mitigable with uniform priors over divergence models (rather than classes of models) and more flexible prior distributions on demographic and divergence‐time parameters.     

Idiosyncratic responses to climate‐driven forest fragmentation and marine incursions in reed frogs from Central Africa and the Gulf of Guinea Islands

Organismal traits interact with environmental variation to mediate how species respond to shared landscapes. Thus, differences in traits related to dispersal ability or physiological tolerance may result in phylogeographic discordance among co‐distributed taxa, even when they are responding to common barriers. We quantified climatic suitability and stability, and phylogeographic divergence within three reed frog species complexes across the Guineo‐Congolian forests and Gulf of Guinea archipelago of Central Africa to investigate how they responded to a shared climatic and geological history. Our species‐specific estimates of climatic suitability through time are consistent with temporal and spatial heterogeneity in diversification among the species complexes, indicating that differences in ecological breadth may partly explain these idiosyncratic patterns. Likewise, we demonstrated that fluctuating sea levels periodically exposed a land bridge connecting Bioko Island with the mainland Guineo‐Congolian forest and that habitats across the exposed land bridge likely enabled dispersal in some species, but not in others. We did not find evidence that rivers are biogeographic barriers across any of the species complexes. Despite marked differences in the geographic extent of stable climates and temporal estimates of divergence among the species complexes, we recovered a shared pattern of intermittent climatic suitability with recent population connectivity and demographic expansion across the Congo Basin. This pattern supports the hypothesis that genetic exchange across the Congo Basin during humid periods, followed by vicariance during arid periods, has shaped regional diversity. Finally, we identified many distinct lineages among our focal taxa, some of which may reflect incipient or unrecognized species.     

Phylogeographic Patterns and Evolution of the Mitochondrial DNA Control Region in Two Neotropical Cats (Mammalia, Felidae)

The ocelot (Leopardus pardalis) and margay (L. wiedii) are sister-species of Neotropical cats which evolved from a lineage that migrated into South America during the formation of the Panamanian land bridge 3–5 million years ago. Patterns of population genetic divergence of each species were studied by phylogenetic analyses of mitochondrial DNA (mtDNA) control region sequences in individuals sampled across the distribution of these taxa. Abundant genetic diversity and remarkably concordant phylogeographic partitions for both species were observed, identifying parallel geographic regions which likely reflect historical faunal barriers. Inferred aspects of phylogeography, population genetic structure, and demographic history were used to formulate conservation recommendations for these species. In addition, observed patterns of sequence variation provided insight into the molecular evolution of the mtDNA control region in closely related felids. Received: 26 January 1998 / Accepted: 14 May 1998     

PHYLOGEOGRAPHY OF THE JUMPING SPIDER HABRONATTUS PUGILLIS (ARANEAE: SALTICIDAE): RECENT VICARIANCE OF SKY ISLAND POPULATIONS?

In island systems with diverging populations, the history of island formation and genealogical estimates of divergence dates can be mutually informative. In the "sky islands" of southeastern Arizona, climate-induced contraction of woodlands appears to have fragmented populations of woodland-dwelling species onto disjunct mountain ranges. Montane populations of the jumping spider, Habronattus pugillis, display striking amounts of phenotypic divergence among ranges. Paleoclimatic estimates date woodland fragmentation at approximately 10,000 years ago, suggesting that phenotypic divergence has been extraordinarily rapid in these spiders. This phylogeographic study of populations of H. pugillis attempts to clarify the species' history of isolation and divergence and to address the suitability of available paleoclimatic data for dating divergences among populations of the region's woodland-dwelling organisms. Mitochondrial sequence data of spiders from 13 mountain ranges was used to reconstruct genealogical relationships. Gene trees show that small mountain ranges tend to have populations whose sequences form monophyletic groups, whereas larger ranges do not. Paraphyly among genes from larger ranges could result from either recent migration or incomplete lineage sorting. I use phylogenetic and geographic information to test these alternatives, and conclude that incomplete lineage sorting best explains the observed paraphyly. Gene trees are concordant with some of the predictions of vegetation history generated by examination of topography. Dates estimated for divergence of populations vary from 30,000 years to more than 2 million years ago, suggesting multiple vicariance events that are older than would be inferred from paleoclimatic studies. These findings illustrate that use of any single paleontological dataset to calibrate molecular clocks can potentially greatly underestimate actual divergence times.     

MULTILOCUS COALESCENCE ANALYSES SUPPORT A mtDNA-BASED PHYLOGEOGRAPHIC HISTORY FOR A WIDESPREAD PALEARCTIC PASSERINE BIRD, SITTA EUROPAEA

Our understanding of species phylogeography in much of the Palearctic is incomplete. In addition, many existing studies based solely on mitochondrial DNA (mtDNA) can provide a biased view of phylogeographic history because of the effects of lineage sorting, natural selection, or hybridization. We analyzed 13 introns to assess a mtDNA study of the Eurasian nuthatch (Sitta europaea) that suggested a seemingly contemporaneous origin of distinct taxa in the Caucasus, Europe, and Asia. Neutrality tests showed no evidence of selection on either the mtDNA or nuclear sequences. Most nuclear gene trees, except for Z-linked ones, did not recover the three lineages, which we attribute to recent splitting. Analyses of the 13 introns combined revealed the same three groups as did the mtDNA and suggested that nuthatches experienced a trichotomous (or two indistinguishable) split(s) 1-2 million years ago (Mya) and have remained isolated with trifling if not zero gene flow since then, and the Asian group increased in population size. This result demonstrates the usefulness of mtDNA in discovering phylogeographic patterns. The use of multiple nuclear loci facilitated detection of an introgressed individual and improved estimates of process parameters such as divergence time and population expansion. We recommend that phylogeographic studies should be based on both mtDNA and nuclear genes.     

Phylogeography and population history of the crab-eating fox (Cerdocyon thous)

The crab-eating fox is a medium-sized Neotropical canid with generalist habits and a broad distribution in South America. We have investigated its genetic diversity, population structure and demographic history across most of its geographic range by analysing 512 base pairs (bp) of the mitochondrial DNA (mtDNA) control region, 615 bp of the mtDNA cytochrome b gene and 1573 total nucleotides from three different nuclear fragments. MtDNA data revealed a strong phylogeographic partition between northeastern Brazil and other portions of the species' distribution, with complete separation between southern and northern components of the Atlantic Forest. We estimated that the two groups diverged from each other c. 400,000-600,000 years ago, and have had contrasting population histories. A recent demographic expansion was inferred for the southern group, while northern populations seem to have had a longer history of large population size. Nuclear sequence data did not support this north-south pattern of subdivision, likely due at least in part to secondary male-mediated historical gene flow, inferred from multilocus coalescent-based analyses. We have compared the inferred phylogeographic patterns to those observed for other Neotropical vertebrates, and report evidence for a major north-south demographic discontinuity that seems to have marked the history of the Atlantic Forest biota.     

Evolutionary history of the greater white-toothed shrew (Crocidura russula) inferred from analysis of mtDNA, Y, and X chromosome markers

We investigate the evolutionary history of the greater white-toothed shrew across its distribution in northern Africa and mainland Europe using sex-specific (mtDNA and Y chromosome) and biparental (X chromosome) markers. All three loci confirm a large divergence between eastern (Tunisia and Sardinia) and western (Morocco and mainland Europe) lineages, and application of a molecular clock to mtDNA divergence estimates indicates a more ancient separation (2.25 M yr ago) than described by some previous studies, supporting claims for taxonomic revision. Moroccan ancestry for the mainland European population is inconclusive from phylogenetic trees, but is supported by greater nucleotide diversity and a more ancient population expansion in Morocco than in Europe. Signatures of rapid population expansion in mtDNA, combined with low X and Y chromosome diversity, suggest a single colonization of mainland Europe by a small number of Moroccan shrews >38 K yr ago. This study illustrates that multilocus genetic analyses can facilitate the interpretation of species' evolutionary history but that phylogeographic inference using X and Y chromosomes is restricted by low levels of observed polymorphism.     

Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America

We constructed phylogenetic hypotheses for Mesoamerican Rhamdia, the only genus of primary freshwater fish represented by sympatric species across Central America. Phylogenetic relationships were inferred from analysis of 1990 base pairs (bp) of mitochondrial DNA (mtDNA), represented by the complete nucleotide sequences of the cytochrome b (cyt b) and the ATP synthase 8 and 6 (ATPase 8/6) genes. We sequenced 120 individuals from 53 drainages to provide a comprehensive geographic picture of Central American Rhamdia systematics and phylogeography. Phylogeographic analysis distinguished multiple Rhamdia mtDNA lineages, and the geographic congruence across evolutionarily independent Rhamdia clades indicated that vicariance has played a strong role in the Mesoamerican diversification of this genus. Phylogenetic analyses of species-level relationships provide strong support for the monophyly of a trans-Andean clade of three evolutionarily equivalent Rhamdia taxa: R. guatemalensis, R. laticauda, and R. cinerascens. Application of fish-based mitochondrial DNA clocks ticking at 1.3-1.5% sequence divergence per million years (Ma), suggests that the split between cis- and trans-Andean Rhamdia extends back about 8 Ma, and the three distinct trans-Andean Rhamdia clades split about 6 Ma ago. Thus the mtDNA divergence observed between cis- and trans-Andean Rhamdia species is too low to support an ancient colonization of Central America in the Late Cretaceous or Paleocene as had been hypothesized in one colonization model for Mesoamerican fishes. Rather the mtDNA data indicate that Rhamdia most likely colonized Central America in the late Miocene or Pliocene, promoting a strong role for the Isthmus of Panamá in the Mesoamerican expansion of this genus. Basal polytomies suggest that both the R. laticauda and R. guatemalensis clades spread rapidly across the Central American landscape, but differences in the average mtDNA genetic distances among clades comprising the two species, indicate that the R. laticauda spread and diversified across Mesoamerica about 1 million years before R. guatemalensis.     

Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach

Background  

Marine allopatric speciation is an enigma because pelagic larval dispersal can potentially connect disjunct populations thereby preventing reproductive and morphological divergence. Here we present a new hierarchical approximate Bayesian computation model (HABC) that tests two hypotheses of marine allopatric speciation: 1.) "soft vicariance", where a speciation involves fragmentation of a large widespread ancestral species range that was previously connected by long distance gene flow; and 2.) peripatric colonization, where speciations in peripheral archipelagos emerge from sweepstakes colonizations from central source regions. The HABC approach analyzes all the phylogeographic datasets at once in order to make across taxon-pair inferences about biogeographic processes while explicitly allowing for uncertainty in the demographic differences within each taxon-pair. Our method uses comparative phylogeographic data that consists of single locus mtDNA sequences from multiple co-distributed taxa containing pairs of central and peripheral populations. We use the method on two comparative phylogeographic data sets consisting of cowrie gastropod endemics co-distributed in the Hawaiian (11 taxon-pairs) and Marquesan archipelagos (7 taxon-pairs).     

Evidence of multiple vicariance in a marine suture‐zone in the Gulf of Mexico

Aim Our aim was to assess whether single or multiple vicariance event(s) have occurred in a marine suture‐zone in the northern Gulf of Mexico. We estimated the divergence time of two genetically distinct groups of lane snapper (Lutjanus synagris), distributed on either side of this zone, and compared it to the timing of other vicariance events described previously in the area. Location The northern Gulf of Mexico from the Texas coastline to the Atlantic coast of South Florida, USA. The marine suture‐zone is proximal to Mobile Bay in Alabama. Methods A Bayesian Markov chain Monte Carlo (MCMC) approach was used to estimate mutation‐corrected divergence and genetic migration between eastern and western groups of lane snapper, based on genotypic differences at 13 nuclear‐encoded microsatellites obtained previously from 77 and 171 individual samples, respectively. A second estimate of divergence time for eastern and western groups of lane snapper was generated, based on 590 base pairs of DNA sequence from the mitochondrially encoded NADH dehydrogenase subunit 4 (ND4) obtained previously from 45 and 93 individual samples, respectively. Results Both classes of genetic markers (microsatellites and mtDNA) indicated a fairly recent divergence. Confidence intervals for microsatellite‐based estimates suggested that divergence began less than 21 thousand years ago (ka), while confidence intervals for mitochondrial DNA‐based estimates suggested divergence began less than 130 ka. These estimates were not consistent with vicariance events in the literature, which are hypothesized to have occurred well before 135 ka. Main conclusions These results indicate that the marine suture‐zone in the northern Gulf of Mexico features multiple vicariance events and may be characterized by a complex geological/environmental history. We suggest that processes during or after the Wisconsin glaciation (c. 110–18 ka) may have created a previously unrecognized barrier for lane snapper and possibly other species as well.     

Molecular evidence for a recent demographic expansion in the puma (Puma concolor) (Mammalia,Felidae)

The puma is an iconic predator that ranges throughout the Americas, occupying diverse habitats. Previous phylogeographic analyses have revealed that it exhibits moderate levels of genetic structure across its range, with few of the classically recognized subspecies being supported as distinct demographic units. Moreover, most of the species’ molecular diversity was found to be in South America. To further investigate the phylogeographic structure and demographic history of pumas we analyzed mtDNA sequences from 186 individuals sampled throughout their range, with emphasis on South America. Our objectives were to refine the phylogeographic assessment within South America and to investigate the demographic history of pumas using a coalescent approach. Our results extend previous phylogeographic findings, reassessing the delimitation of historical population units in South America and demonstrating that this species experienced a considerable demographic expansion in the Holocene, ca. 8,000 years ago. Our analyses indicate that this expansion occurred in South America, prior to the hypothesized re-colonization of North America, which was therefore inferred to be even more recent. The estimated demographic history supports the interpretation that pumas suffered a severe demographic decline in the Late Pleistocene throughout their distribution, followed by population expansion and re-colonization of the range, initiating from South America.     

Living with the genetic signature of Miocene induced change: evidence from the phylogeographic structure of the endemic angulate tortoise Chersina angulata

The phylogeographic structure of the monotypic endemic southern African angulate tortoise Chersina angulata was investigated throughout its distribution with the use of partial sequences from three mtDNA loci (COI, cyt b and ND4). Phylogeographic and phylogenetic structuring obtained for the three mtDNA markers were highly congruent and suggested the presence of two genetically distinct, reciprocally monophyletic evolutionary lineages. Group one contained two subclades with haplotypes from the north-western Cape and south-western Cape, respectively, while haplotypes from the southern Cape comprised group two. The two major clades were separated by nine and eight mutational steps for COI and ND4, respectively. Of the three mtDNA gene regions examined, the ND4 partial sequence contained the most phylogenetic signal. Haplotype diversity was generally low and we recovered 34 haplotypes for the 125 animals sequenced for the ND4 subunit. Nested clade analyses performed on the variable ND4 partial sequences suggested the presence of two major refugial areas for this species. The demographic history of the taxon was characterised by range expansion and prolonged historical fragmentation. Divergence time estimates suggest that the temporal and spatial distribution of the taxon was sculpted by changes in temperature and rainfall patterns since the late Miocene. Corroborative evidence from other reptiles is also suggestive of a late Miocene divergence, indicating that this was a major epoch for cladogenesis in southern Africa. Apart from the genetic differences between the two major clades, we also note morphometric and behavioural differences, alluding to the presence of two putative taxa nested within C. angulata.     

Mechanisms of global diversification in the brown booby (Sula leucogaster) revealed by uniting statistical phylogeographic and multilocus phylogenetic methods

Recent theoretical and empirical research suggests that statistical models based on coalescent theory can improve both phylogeographic and phylogenetic inference. An approach that involves elements of both statistical phylogeography (e.g. Isolation with Migration analyses) and multilocus phylogenetic inference (e.g. *beast) may be particularly useful when applied to populations with relatively old divergence times. Here, we use such an approach in the globally distributed brown booby (Sula leucogaster). We sampled 215 individuals from all major breeding areas and genotyped them at eight microsatellite and three nuclear intron loci. We found that brown booby populations were highly differentiated and that colonies can be grouped into four major genetic populations (Caribbean Sea, Central Atlantic Ocean, Indo-Central Pacific and Eastern Pacific). These populations apparently diverged in the absence of gene flow and, with one exception, currently exchange few to no migrants. The Eastern Pacific population diverged from all other populations approximately one million years ago [90% highest posterior density: 330,000-2,000,000 years ago] and exhibits a distinct male plumage, relative to other populations. However, recent gene flow from the Indo-Central Pacific into the Eastern Pacific appears to have occurred, suggesting that approximately one million years of genetic isolation and divergence in male plumage colour are not sufficient to prevent interbreeding. Gene flow following secondary contact of the Indo-Central Pacific and Eastern Pacific populations was not detected in previous mitochondrial DNA (mtDNA) studies, and the contrast between the mtDNA results and our current results highlights the advantage of a multilocus phylogeographic approach.     

msBayes: Pipeline for testing comparative phylogeographic histories using hierarchical approximate Bayesian computation

Background  

Although testing for simultaneous divergence (vicariance) across different population-pairs that span the same barrier to gene flow is of central importance to evolutionary biology, researchers often equate the gene tree and population/species tree thereby ignoring stochastic coalescent variance in their conclusions of temporal incongruence. In contrast to other available phylogeographic software packages, msBayes is the only one that analyses data from multiple species/population pairs under a hierarchical model.     

Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

The Western Lyresnake (Trimorphodon biscutatus) is a widespread, polytypic taxon inhabiting arid regions from the warm deserts of the southwestern United States southward along the Pacific versant of Mexico to the tropical deciduous forests of Mesoamerica. This broadly distributed species provides a unique opportunity to evaluate a priori biogeographical hypotheses spanning two major distinct biogeographical realms (the Nearctic and Neotropical) that are usually treated separately in phylogeographical analyses. I investigated the phylogeography of T. biscutatus using maximum likelihood and Bayesian phylogenetic analysis of mitochondrial DNA (mtDNA) from across this species' range. Phylogenetic analyses recovered five well-supported clades whose boundaries are concordant with existing geographical barriers, a pattern consistent with a model of vicariant allopatric divergence. Assuming a vicariance model, divergence times between mitochondrial lineages were estimated using Bayesian relaxed molecular clock methods calibrated using geological information from putative vicariant events. Divergence time point estimates were bounded by broad confidence intervals, and thus these highly conservative estimates should be considered tentative hypotheses at best. Comparison of mtDNA lineages and taxa traditionally recognized as subspecies based on morphology suggest this taxon is comprised of multiple independent lineages at various stages of divergence, ranging from putative secondary contact and hybridization to sympatry of 'subspecies'.     

Diversification across the Palaearctic desert belt throughout the Pleistocene: phylogeographic history of the Houbara–Macqueen's bustard complex (Otididae: Chlamydotis) as revealed by mitochondrial DNA

Studies on the influence of Pleistocene climatic fluctuations and associated habitat changes on arid‐adapted bird species living in the Holarctic region are comparatively rare. In contrast to temperate species, the populations of arid‐adapted avian species might be characterized by low genetic differentiation because periods of population isolation were associated with the short interglacial periods, while population expansion events might have occurred during the longer glacial periods when steppe‐like vegetation might have been prevalent. In this study, we tested this hypothesis in a widespread arid‐adapted taxon of the Palaearctic desert belt, the Houbara–Macqueen's bustard complex. The later includes the Houbara bustard Chlamydotis undulata, comprising the North African subspecies Chlamydotis u. undulata and Chlamydotis u. fuertaventurae from the Canary Islands, and the Asian Macqueen's bustard Chlamydotis macqueenii. A long fragment (1042 bp) of the Cyt‐b gene was investigated in 39 representatives of the two species to assess phylogenetic and phylogeographic patterns, and demographic history and to compute divergence time estimates using a Bayesian relaxed molecular clock approach based on different coalescent priors. While the two species are genetically distinct, we found little intraspecific genetic differentiation. The divergence time of the two species falls within a period of extreme aridity at around 0.9 million years ago, which most likely resulted in an east–west vicariance along the Arabo‐Saharan deserts. Differentiation within Houbara and Macqueen's bustard occurred later during the Middle to Upper Pleistocene, and as we have predicted, periods of range expansion were associated to the last glacial period at least in the Macqueen's bustard.     

Geographical structuring of mitochondrial DNA in Chalcides sexlineatus within the island of Gran Canaria

The skink Chalcides sexlineatus shows clear north-south ecology-correlated variation in morphology within the island of Gran Canaria. Detailed study of this variation previously supported the hypothesis of in situ differentiation arising through climate-mediated variation in selection pressures. The more recent discovery of substantial within-island mitochondrial DNA (mtDNA) sequence divergence suggested an alternative historical vicariance scenario. This paper provides a detailed analysis of geographical structuring in mtDNA, based on a fragment amplified from the 12S rRNA region in 96 individuals. A very clear phylogeographic pattern was detected: three ''deep'' lineages were evident corresponding to northern, south-eastern and south-western parts of Gran Canaria with generally low levels of mtDNA introgression. The pattern of among-site differentiation was highly concordant with the pattern of morphological variation. It also provided a close fit to a simple microevolutionary model based on population vicariance during volcanic eruptions known to have occurred during the past 2.8 million years in the north of the island. The minimum number of historical migration events required to explain the mtDNA tree was calculated and a novel randomization test used to show that there were a lower number of putative migrations across the morphology transition zone, relative to migrations within northern and southern zones. Differential adaptation to the two major habitat types and selection against hybrids may explain why the morphology transition zone coincides with the main ecotone and is relatively narrow, respectively.     

Pleistocene effects on North American songbird evolution

Recent studies have used comparisons of mitochondrial DNA (mtDNA) sequence divergence among populations and species to test existing hypotheses about avian evolution during the Pleistocene epoch. In 1998, Avise and Walker concluded that the Pleistocene was an important time for avian evolution, including the initiation of phylogeographic separations and the completion of speciation events that began in the Pliocene. The study implied that these conclusions conflicted with the study, in 1997, by Klicka and Zink, which concluded that most species pairs previously thought to have originated in the past 250 000 years were much older. The two studies are complementary in the sense that Avise and Walker dealt primarily with phylogeographic (intraspecific) separations. Furthermore, Klicka and Zink concentrated on the inception of divergences whereas the Avise and Walker focused on the timing of the completion of speciation. To accomplish this, Avise and Walker analysed ''phylogroups'', geographically coherent subsets of biological species in which mtDNA haplotypes exhibit reciprocal monophyly. The study used the average interphylogroup mtDNA distance (0.027), calibrated at 2% per million years, to conclude that speciation required on average one million years to complete. Hence, speciation events begun in the Late Pliocene would have been completed in the mid- to late Pleistocene. Although we appreciate the extended nature of the speciation process and Avise and Walker''s insightful attempt to estimate its duration, we conclude that their value was an overestimate by a factor of two. In particular we question whether phylogroups can be used in the novel evolutionary role that Avise and Walker envisioned, because of the vagaries of taxonomic practices and lack of consensus regarding species concepts. To extend their analysis of intraspecific, phylogeographic separations, we compiled previously analysed and newly available data for divergence times for North American songbird (order Passeriformes) phylogroups. More than 80% were initiated at least one million years ago, which is inconsistent with the late Pleistocene origins model previously rejected by Klicka and Zink mentioned above. Although some divergence events can be traced to the late Pleistocene, the significance of the distribution must be judged with reference to a null model. Whether the Pleistocene was a profound time for avian phylogeographic differentiation is at present unknown.     

Nuclear and mitochondrial phylogeography of the Atlantic forest endemic Xiphorhynchus fuscus (Aves: Dendrocolaptidae): biogeography and systematics implications

We studied the intraspecific evolutionary history of the South American Atlantic forest endemic Xiphorhynchusfuscus (Aves: Dendrocolaptidae) to address questions such as: Was the diversification of this bird's populations associated to areas of avian endemism? Which models of speciation (i.e., refuges, river as barriers or geotectonism) explain the diversification within X. fuscus? Does the genetic data support subspecies as independent evolutionary units (species)? We used mitochondrial (n=34) and nuclear (n=68) DNA sequences of X. fuscus to study temporal and spatial relationships within and between populations. We described four main monophyletic lineages that diverged during the Pleistocene. The subspecies taxonomy did not match all the evolutionary lineages; subspecies atlanticus was the only one that represented a monophyletic and isolated lineage. The distribution of these lineages coincided with some areas of endemism for passerines, suggesting that those areas could be regions of biotic differentiation. The ancestor of X. fuscus diverged approximately 3 million years ago from Amazonian taxa and the phylogeographic pattern suggested that X. fuscus radiated from northeastern Brazil. Neither the riverine nor the geotectonic vicariance models are supported as the primary cause for diversification of geographic lineages, but rainforest contractions and expansions (ecological vicariance) can explain most of the spatial divergence observed in this species. Finally, analyses of gene flow and divergence time estimates suggest that the endangered subspecies atlanticus (from northeastern Brazil) can be considered a full species under the general lineage species concept.