Testing nested phylogenetic and phylogeographic hypotheses in the Plethodon vandykei species group

Mesic forests in the North American Pacific Northwest occur in two disjunct areas: along the coastal and Cascade ranges of Oregon, Washington, and British Columbia as well as the Northern Rocky Mountains of Idaho, Montana, and British Columbia. Over 150 species or species complexes have disjunct populations in each area, and a priori hypotheses based on phytogeography and geology potentially explain the disjunction via either dispersal or vicariance. Here, we test these hypotheses in the disjunct salamander complex Plethodon vandykei and P. idahoensisby collecting genetic data (669 bp of Cyt b) from 262 individuals. Maximum likelihood analysis indicated reciprocal monophyly of these species, supporting the ancient vicariance hypothesis, whereas parametric bootstrap and Bayesian hypothesis testing allow rejection of the dispersal hypothesis. The coalescent estimate of the time since population divergence (estimated using MDIV) is 3.75 x 106 years, and the 95%credibility interval of this value overlaps with the geological estimate of vicariance, but not the hypothesized dispersal. These results are congruent with the pattern seen in other mesic forest amphibian lineages and suggest disjunction in amphibians may be a concerted response to a geological/climatological event. WithinP. idahoensis, we tested the corollary hypothesis of an inland Pleistocene refugium in the Clearwater drainage with nested clade analysis and coalescent estimates of population growth rate (g). Both analyses support post-Pleistocene expansion from the Clearwater refugium. We corroborated this result by calculating Tajima's Dand mismatch distribution within each drainage, showing strong evidence for recent population expansion within most drainages. This work demonstrates the utility of statistical phylogeography and contributes two novel analytical tools: tests of stationarity with respect to topology in the Bayesian estimation, and the use of coalescent simulations to test the significance of the population growth-rate parameter.     

Testing hypotheses of speciation timing in Dicamptodon copei and Dicamptodon aterrimus (Caudata: Dicamptodontidae)

Giant salamanders of the genus Dicamptodon are members of the mesic forest ecosystem that occurs in the Pacific Northwest of North America. We estimate the phylogeny of the genus to test several hypotheses concerning speciation and the origin of current species distributions. Specifically, we test competing a priori hypotheses of dispersal and vicariance to explain the disjunct inland distribution of the Idaho giant salamander (D. aterrimus) and to test the hypothesis of Pleistocene speciation of Cope's giant salamander (D. copei) using Bayesian hypothesis testing. We determined that available outgroups were too divergent to root the phylogeny effectively, and we calculated Bayesian posterior probabilities for each of the 15 possible root placements for this four-taxon group. This analysis placed the root on the branch leading to D. aterrimus, indicating that current distribution and speciation of D. aterrimus fit the ancient vicariance hypothesis and are attributable to the orogeny of the Cascade Mountains rather than recent inland dispersal. Furthermore, test results indicate that D. copei is distantly related to other coastal lineages and likely originated much earlier than the Pleistocene. These results suggest that speciation within the genus is attributable to ancient geologic events, while more recent Pleistocene glaciation has shaped genetic variation and distributions within the extant species.     

Phylogeographic structure and cryptic speciation in the trans-Antarctic moss Pyrrhobryum mnioides

Abstract Many bryophyte species have distributions that span multiple continents. The hypotheses historically advanced to explain such distributions rely on either long-distance spore dispersal or slow rates of morphological evolution following ancient continental vicariance events. We use phylogenetic analyses of DNA sequence variation at three chloroplast loci ( atpB-rbcL spacer, rps4 gene, and trnL intron and 3'spacer) to examine these two hypotheses in the trans-Antarctic moss Pyrrhobryum mnioides. We find: (1) reciprocal monophyly of Australasian and South American populations, indicating a lack of intercontinental dispersal; (2) shared haplotypes between Australia and New Zealand, suggesting recent or ongoing migration across the Tasman Sea; and (3) reciprocal monophyly among Patagonian and neotropical populations, suggesting no recent migration along the Andes. These results corroborate experimental work suggesting that spore features may be critical determinants of species range. We use the mid-Miocene development of the Atacama Desert, 14 million years ago, to calibrate a molecular clock for the tree. The age of the trans-Antarctic disjunction is estimated to be 80 million years ago, consistent with Gondwanan vicariance, making it among the most ancient documented cases of cryptic speciation. These data are in accord with niche conservatism, but whether the morphological stasis is a product of stabilizing selection or phylogenetic constraint is unknown.     

Origin and evolution of the northern hemisphere disjunction in the moss genus Homalothecium (Brachytheciaceae)

Competing hypotheses that rely either on a stepping-stone dispersal via the North Atlantic or the Bering land bridges, or more recent transoceanic dispersal, have been proposed to explain the disjunct distribution of Mediterranean flora in southern Europe and western North America. These hypotheses were tested with molecular dating using a phylogeny of the moss genus Homalothecium based on ITS, atpB-rbcL, and rpl16 sequence data. The monophyly of two main lineages in Western Palearctic (Europe, central Asia and north Africa) and North America is consistent with the ancient vicariance hypothesis. The monophyly of Madeiran H. sericeum accessions supports the recognition of the Macaronesian endemic H. mandonii. A range of absolute rates of molecular evolution documented in land plants was used as probabilistic calibration prior by a Bayesian inference implementing a relaxed-clock model to derive ages for the nodes of interest. Our age estimates for the divergence of the American and Western Palearctic Homalothecium clade (5.7 Ma, IC 3.52-8.26) and the origin of H. mandonii (2.52 Myr IC 0.86-8.25) are not compatible with the ancient vicariance hypothesis. Age estimates suggests that species distributions result from rare instances of dispersal and subsequent sympatric diversification. The calibrated phylogeny indicates that Homalothecium has undergone a fast radiation during the last 4 Myr, which is consistent with the low levels of morphological divergence among sibling species.     

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.     

The evolution of Queensland spiny mountain crayfish of the genus Euastacus. I. Testing vicariance and dispersal with interspecific mitochondrial DNA

The upland mesic rainforests of eastern Australia have been described as a "mesothermal archipelago" where a chain of cool mountain "islands" arise from a warm "sea" of tropical and subtropical lowlands. An endemic freshwater crayfish belonging to the genus Euastacus is found on each of these mountain "islands." The Euastacus are particularly suitable for the study of evolution because each mountain harbors a unique species, there are many taxa present providing replication within the group and, most importantly, their distribution is linear, extending along a south-north axis. This group could have evolved by "simultaneous vicariance" where there was one vicariant separation event of a widespread ancestor, or by "south to north stepping stone dispersal" where there were long distance dispersal events from neighboring mountain islands, starting in the south and proceeding north in a dispersal-colonization wave. We used pairwise genetic distances between nearest geographic neighbors as a novel way to test the two hypotheses. If diversification was due to "south to north stepping stone dispersal," then pairwise genetic distances between nearest geographic neighbors should decrease progressively the farther north the taxon pairs are found, reflecting the decreasing periods of isolation. In this case there should be a negative correlation between the south to north rank order of nearest neighbors and pairwise genetic distances. A Spearman's correlation on 16S mtDNA pairwise genetic distances and geographic rank order was not significant, indicating there was no support for the south to north stepping stone dispersal hypothesis. If simultaneous vicariance was responsible for diversification then all nearest geographic neighbor taxon pairs should have similar genetic distances and, therefore, the variance in nearest neighbor distances should be zero, or close to it. To test if the observed variance was tending towards zero we developed a randomization test where nearest neighbor taxon pairs were assigned random genetic distances and the variances calculated. The observed variance lay in the < 0.05 range of the simulated variances, providing support for the simultaneous vicariance hypothesis. The data also suggest there was simultaneous vicariance of at least two ancestral Queensland lineages. The timing of this vicariant event was probably in the Pliocene, which is consistent with the divergence times reported for other Australian mesic rainforest restricted taxa.     

Decline of a biome: evolution,contraction, fragmentation,extinction and invasion of the Australian mesic zone biota

Aim The mesic biome, encompassing both rain forest and open sclerophyllous forests, is central to understanding the evolution of Australia’s terrestrial biota and has long been considered the ancestral biome of the continent. Our aims are to review and refine key hypotheses derived from palaeoclimatic data and the fossil record that are critical to understanding the evolution of the Australian mesic biota. We examine predictions arising from these hypotheses using available molecular phylogenetic and phylogeographical data. In doing so, we increase understanding of the mesic biota and highlight data deficiencies and fruitful areas for future research. Location The mesic biome of Australia, along the eastern coast of Australia, and in the south‐east and south‐west, including its rain forest and sclerophyllous, often eucalypt‐dominated, habitats. Methods We derived five hypotheses based on palaeoclimatic and fossil data regarding the evolution of the Australian mesic biota, particularly as it relates to the mesic biome. We evaluated predictions formulated from these hypotheses using suitable molecular phylogenies of terrestrial plants and animals and freshwater invertebrates. Results There was support for the ancestral position of mesic habitat in most clades, with support for rain forest habitat ancestry in some groups, while evidence of ancestry in mesic sclerophyllous habitats was also demonstrated for some plants and herpetofauna. Contraction of mesic habitats has led to extinction of numerous lineages in many clades and this is particularly evident in the rain forest component. Species richness was generally higher in sclerophyllous clades than in rain forest clades, probably due to higher rates of net speciation in the former and extinction in the latter. Although extinction has been prominent in rain forest communities, tropical rain forests appear to have experienced extensive immigration from northern neighbours. Pleistocene climatic oscillations have left genetic signatures at multiple levels of divergence and with complex geographical structuring, even in areas with low topographical relief and few obvious geographical barriers. Main conclusions Our review confirms long‐held views of the ancestral position of the Australian mesic biome but also reveals new insights into the complexity of the processes of contraction, fragmentation, extinction and invasion during the evolution of this biome.     

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

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

The biogeographic history of the relictual Gondwanan lineage of Australian burrowing crayfish

Biogeographic investigations of Gondwanan mesic Australian fauna are scarce. The burrowing clade of Australian freshwater crayfish represent an ideal group to provide biogeographic inferences, due to their extensive distribution across the continent and their presumed ancient origin. This study tested the competing hypotheses of a ‘early’ versus ‘late’ origin of this clade, coinciding with the early or late fragmentation of Gondwana, respectively. The biogeographic history of this group was investigated through: (a) examination of the phylogenetic relationships between the seven extant taxon groups; (b) reconstruction of four species trees, each using a different calibration method; and (c) reconstruction of ancestral ranges and correlation of estimated dispersal and vicariance events with historical geological data to propose plausible mechanisms responsible for driving diversification. The phylogenetic relationships between the taxon groups were generally well supported (although some uncertainty exists for the oldest genera), and all calibration methods produced concordant results. The hypothesis that the clade arose during the early fragmentation of Gondwana in southern Australia is supported. Divergence between the extant taxa likely resulted from a combination of both short- and long-distance dispersal events (often followed by later vicariance), coincident with phases of sea level oscillation and changing climate continuing into the Eocene.

     

Intraspecific phylogeography of Percina evides (Percidae: Etheostomatinae): an additional test of the Central Highlands pre-Pleistocene vicariance hypothesis

North America exhibits the most diverse freshwater fish fauna among temperate regions of the world. Species diversity is concentrated in the Central Highlands, drained by the Mississippi, Gulf Slope and Atlantic Slope river systems. Previous investigations of Central Highlands biogeography have led to conflicting hypotheses involving dispersal and vicariance to explain the diversity and distribution of the freshwater fish fauna. In this investigation predictions of the Central Highlands pre-Pleistocene vicariance hypothesis are tested with a phylogeographic analysis of the percid species Percina evides, which is widely distributed in several disjunct areas of the Central Highlands. Phylogenetic analysis of complete gene sequences of mitochondrially encoded cytochrome b recover three phylogroups, with very low levels of sequence polymorphism within groups. The two western phylogroups are monophyletic with respect to the eastern phylogroup. The recovery of two monophyletic lineages with an eastern and western distribution in the disjunct highland areas is a pattern expected from vicariance, but is not predicted by the Central Highlands pre-Pleistocene vicariance hypothesis. The recovery of very limited mitochondrial DNA polymorphism and lack of phylogeographic structuring across the entire range of the eastern clade, very shallow polymorphism between the disjunct Missouri River and upper Mississippi River populations, and lack of sequence polymorphism in the upper Mississippi River populations, support a hypothesis of dispersal during or following the Pleistocene. The present distribution of P. evides is best explained by both vicariant and dispersal events.     

Timing the eastern Asian-eastern North American floristic disjunction: molecular clock corroborates paleontological estimates

Sequence data of the chloroplast gene rbcL were used to estimate the time of the well-known eastern Asian-eastern North American floristic disjunction. Sequence divergence of rbcL was examined for 22 species of 11 genera (Campsis, Caulophyllum, Cornus, Decumaria, Liriodendron, Menispermum, Mitchella, Pachysandra, Penthorum, Podophyllum, and Phryma) representing a diverse array of flowering plants occurring disjunctly in eastern Asia and eastern North America. Divergence times of putative disjunct species pairs were estimated from synonymous substitutions, using rbcL molecular clocks calibrated for Cornus. Relative rate tests were performed to assess rate constancy of rbcL evolution among lineages. Corrections of estimates of divergence times for each species pair were made based on rate differences of rbcL between Cornus and other species pairs. Results of these analyses indicate that the time of divergence of species pairs examined ranges from 12.56 +/- 4.30 million years to recent (<0.31 million years), with most within the last 10 million years (in the late Miocene and Pliocene). These results suggest that the isolation of most morphologically similar disjunct species in eastern Asia and eastern North America occurred during the global climatic cooling period that took place throughout the late Tertiary and Quaternary. This estimate is closely correlated with paleontological evidence and in agreement with the hypothesis that considers the eastern Asian-eastern North American floristic disjunction to be the result of the range restriction of a once more or less continuously distributed mixed mesophytic forest of the Northern Hemisphere that occurred during the late Tertiary and Quaternary. This implies that in most taxa the disjunction may have resulted from vicariance events. However, long-distance dispersal may explain the disjunct distribution of taxa with low divergence, such as Menispermum.     

The potential of paleoecology for functional forest restoration planning: lessons from Late Holocene Italian pollen records

We describe forest landscape transformations during the last two millennia in the Italian peninsula by analyzing local (Rieti basin – Lago Lungo) and regional (RF93-30 Adriatic Sea) sediment cores. We identify a dynamic forest ecosystem through paleoecologic reconstruction and consider potential interventions for effective restoration of the most ancient, least disturbed forest ecosystem. The most degraded ecosystems in consequence of human activities were hygrophilous (wet) and mesic forests. In the Rieti Basin, degraded forest ecosystems on mountain slopes are undergoing some degree of forest succession and have less need of restoration. However, management plans for biodiversity, ecosystem services and resources conservation are needed to achieve more sustainable development. In Rieti, the paleoecological investigation revealed a dramatic decrease of deciduous wet and mesic tree taxa through time due to human landscape transformation. The starting point for restoration of a Mediterranean forest ecosystem that preserves natural biodiversity and associated ecosystem services requires recreating some portion of the floodplain wetland ecological niche. Once floodplain forest ecological niche has been recreated, the original ecosystem composed of Alnus, Fraxinus excelsior, Tilia spp., Carpinus betulus and Acer spp., all species which today are rare, should be planted on the basis of microsite characteristics and tree autoecology.     

A multi-compartmented glacial refugium in the northern Rocky Mountains: Evidence from the phylogeography of <Emphasis Type="Italic">Cardamine constancei</Emphasis> (Brassicaceae)

The age and origin of the mesic coniferous forest ecosystem of the Pacific Northwest of North America have long been the subject of debate by biogeographers. Cardamine constancei, an endemic of the Rocky Mountain segment of this ecosystem, was subjected to phylogeographic analysis to test explicit hypotheses on the age of the ecosystem. We have predicted genetic homogeneity among river drainages if C. constancei and other associated species migrated into the region after glaciation, in contrast to the genetic differentiation that may have accrued if the species and its ecosystem have long survived in the relatively warm river canyons south of glaciation. We detected 19 haplotypes with divergence up to 1.5%, and they comprise 4 well-differentiated cpDNA clades. These clades are allopatric except for two haplotypes from the lower Clearwater clade that appear to have dispersed north into partial sympatry with the clade endemic to St.␣Joe River. The divergence and distribution of these clades is consistent with the existence of a complex glacial refugium with at least four compartments. The surprisingly high cpDNA diversity within this species suggests that conservation of mesic coniferous forest ecosystems in the region warrant a conservation plan that accounts for the historically imposed spatial structure of genetic diversity. We are currently testing our phylogeographic hypotheses by the comparative analyses of a suite of plants, animals and fungi.     

Trans-oceanic and endemic origins of the small minnow mayflies (Ephemeroptera, Baetidae) of Madagascar

We investigated the relative importance of dispersal and vicariance in forming the Madagascar insect fauna, sequencing approximately 2300bp from three rRNA gene regions to investigate the phylogeny of Afrotropical small minnow mayflies (Ephemeroptera: Baetidae). Six lineages contained trans-oceanic sister taxa, and variation in genetic divergence between sister taxa revealed relationships that range from very recent dispersal to ancient vicariance. Dispersal was most recent and frequent in species that spend the larval stage in standing water, adding to evidence that these evolutionarily unstable habitats may select for ecological traits that increase dispersal in insects. Ancestral state likelihood analysis suggested at least one Afrotropical lineage had its origin in Madagascar, demonstrating that unidirectional dispersal from a continental source may be too simplistic. We conclude that the Malagasy mayfly fauna should be considered in a biogeographical context that extends beyond Madagascar itself, encompassing trans-oceanic dispersal within multiple lineages.     

Molecular systematics of the trans-Pacific alpine genus Oreomyrrhis (Apiaceae): phylogenetic affinities and biogeographic implications

The alpine ecosystem is the only terrestrial biogeographic unit that is distributed globally. Studying phylogenetics of the plant species in this widespread ecosystem can provide insights into the historical biogeographic processes that have shaped the global biodiversity. The trans-Pacific disjunct alpine genus Oreomyrrhis (Apiaceae) was investigated using nrDNA ITS sequences to test the taxonomic and biogeographic hypotheses. Phylogenetic analyses using maximum parsimony, maximum likelihood, and Bayesian inference revealed that species of Oreomyrrhis form a weakly supported monophyletic clade that is nested within Chaerophyllum sect. Chaerophyllum (subtribe Scandicinae, tribe Scandiceae). The optimal solutions of dispersal-vicariance analysis indicate that the ancestor of Chaerophyllum sect. Chaerophyllum (including Oreomyrrhis) was distributed in Eurasia and subsequently dispersed to North America and southern Pacific Rim. Based on dating using ITS sequence variation, these dispersal events were most likely recent, probably during late Tertiary to Quaternary. The structure of the ITS haplotype network suggests that a rapid range expansion via long-distance dispersal had been crucial in generating the trans-Pacific disjunction of Oreomyrrhis. Furthermore, evolution toward smaller mericarp size and a transition from outcrossing to selfing during Oreomyrrhis's evolution might have increased the chances for long-distance dispersal, facilitating its range expansion and occupation on alpine environments.     

A peninsula as an island: multiple forms of evidence for overwater colonization of Baja California by the gartersnake Thamnophis validus

The recent shift toward dispersal rather than vicariant explanations of disjunct distributions has been driven by the use of molecular data to estimate divergence dates between lineages. However, other kinds of evidence can also be critical in evaluating such biogeographic hypotheses. In the present study, we used a multifaceted approach employing diverse analyses of mitochondrial DNA sequences to assess explanations for the disjunct distribution of the gartersnake Thamnophis validus. The occurrence of this species in the Cape Region of the Baja California peninsula, isolated from mainland populations that occur along the west coast of Mexico, might be explained by: (1) separation of the peninsula from mainland Mexico through rifting 4–8 Mya (tectonic vicariance); (2) fragmentation of the range of this semi‐aquatic species because of post‐Pleistocene aridification (vicariance by aridification); (3) natural overwater dispersal across the Gulf of California; or (4) human introduction. Divergence dating indicates that peninsular and mainland T. validus separated from each other within the last 0.5 Myr, thus rejecting tectonic vicariance. In addition, the estimated closest mainland relatives of peninsular snakes are found farther north than expected under this hypothesis. Three findings argue against vicariance by aridification: (1) peninsular snakes and their closest mainland relatives are more genetically similar than predicted; (2) the location of closest mainland relatives is farther south than predicted; and (3) the species is absent from areas where one might expect to find relict populations. Taken together, refutation of the vicariance hypotheses and the fact that the estimated closest mainland relatives are found almost directly across the Gulf from the Cape Region supports some form of overwater colonization. Various additional arguments suggest that natural dispersal is more likely than human introduction. The present study emphasizes the need for multiple kinds of evidence, beyond divergence dates, to discriminate among hypotheses and to provide independent sources of corroboration or refutation in historical biogeography. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 409–424.     

Paleogene radiation of a plant pathogenic mushroom

Background

The global movement and speciation of fungal plant pathogens is important, especially because of the economic losses they cause and the ease with which they are able to spread across large areas. Understanding the biogeography and origin of these plant pathogens can provide insights regarding their dispersal and current day distribution. We tested the hypothesis of a Gondwanan origin of the plant pathogenic mushroom genus Armillaria and the currently accepted premise that vicariance accounts for the extant distribution of the species.

Methods

The phylogeny of a selection of Armillaria species was reconstructed based on Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI). A timeline was then placed on the divergence of lineages using a Bayesian relaxed molecular clock approach.

Results

Phylogenetic analyses of sequenced data for three combined nuclear regions provided strong support for three major geographically defined clades: Holarctic, South American-Australasian and African. Molecular dating placed the initial radiation of the genus at 54 million years ago within the Early Paleogene, postdating the tectonic break-up of Gondwana.

Conclusions

The distribution of extant Armillaria species is the result of ancient long-distance dispersal rather than vicariance due to continental drift. As these finding are contrary to most prior vicariance hypotheses for fungi, our results highlight the important role of long-distance dispersal in the radiation of fungal pathogens from the Southern Hemisphere.     

Biogeography and diversification of colletid bees (Hymenoptera: Colletidae): emerging patterns from the southern end of the world

Aim The evolutionary history of bees is presumed to extend back in time to the Early Cretaceous. Among all major clades of bees, Colletidae has been a prime example of an ancient group whose Gondwanan origin probably precedes the complete break‐up of Africa, Antarctica, Australia and South America, because modern lineages of this family occur primarily in southern continents. In this paper, we aim to study the temporal and spatial diversification of colletid bees to better understand the processes that have resulted in the present southern disjunctions. Location Southern continents. Methods We assembled a dataset comprising four nuclear genes of a broad sample of Colletidae. We used Bayesian inference analyses to estimate the phylogenetic tree topology and divergence times. Biogeographical relationships were investigated using event‐based analytical methods: a Bayesian approach to dispersal–vicariance analysis, a likelihood‐based dispersal–extinction–cladogenesis model and a Bayesian model. We also used lineage through time analyses to explore the tempo of radiations of Colletidae and their context in the biogeographical history of these bees. Results Initial diversification of Colletidae took place at the Late Cretaceous (≥ 70 Ma). Several (6–14) lineage exchanges between Australia and South America via Antarctica during the Late Cretaceous and Eocene epochs could explain the disjunctions observed between colletid lineages today. All biogeographical methods consistently indicated that there were multiple lineage exchanges between South America and Australia, and these approaches were valuable in exploring the degree of uncertainty inherent in the ancestral reconstructions. Biogeographical and dating results preclude an explanation of Scrapterinae in Africa as a result of vicariance, so one dispersal event is assumed to explain the disjunction in relation to Euryglossinae. The net diversification rate was found to be highest in the recent history of colletid evolution. Main conclusions The biogeography and macroevolutionary history of colletid bees can be explained by a combination of Cenozoic vicariance and palaeoclimatic changes during the Neogene. The austral connection and posterior break‐up of South America, Antarctica and Australia resulted in a pattern of disjunct sister lineages. Increased biome aridification coupled with floristic diversification in the southern continents during the Neogene may have contributed to the high rates of cladogenesis in these bees in the last 25–30 million years.     

Glacial vicariance in the Pacific Northwest: evidence from a lodgepole pine mitochondrial DNA minisatellite for multiple genetically distinct and widely separated refugia

The Canadian side of the Pacific Northwest was almost entirely covered by ice during the last glacial maximum, which has induced vicariance and genetic population structure for several plant and animal taxa. Lodgepole pine ( Pinus contorta Dougl. ex. Loud.) has a wide latitudinal and longitudinal distribution in the Pacific Northwest. Our main objective was to identify relictual signatures of glacial vicariance in the population structure of the species and search for evidence of distinct glacial refugia in the Pacific Northwest. A maternally inherited mitochondrial DNA minisatellite-like marker was used to decipher haplotype diversity in 91 populations of lodgepole pine located across the natural range. Overall population differentiation was sizeable ( G _ST = 0.365 and R _ST = 0.568). Four relatively homogeneous groups of populations, possibly representative of as many genetically distinct glacial populations, were identified for the two main subspecies, ssp. latifolia and ssp. contorta . For ssp. contorta , one glacial lineage is suggested to have been located at high latitudes and possibly off the coast of mainland British Columbia (BC), while the other is considered to have been located south of the ice sheet along the Pacific coast. For ssp. latifolia , two genetically distinct glacial populations probably occurred south of the ice sheet: in the area bounded by the Cascades and Rocky Mountains ranges, and on the eastern side of the Rockies. A possible fifth refugium located in the Yukon may have also been present for ssp. latifolia . Zones of contact between these ancestral lineages were also apparent in interior and northern BC. These results indicate the role of the Queen Charlotte Islands and the Alexander Archipelago as a refugial zone for some Pacific Northwest species and the vicariant role played by the Cascades and the American Rocky Mountains during glaciation.     

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.