Beringian origins and cryptic speciation events in the fly agaric (Amanita muscaria)

Amanita muscaria sensu lato has a wide geographic distribution, occurring in Europe, Asia, Africa, Australia, New Zealand, and North, Central and South America. Previous phylogenetic work by others indicates three geographic clades (i.e. 'Eurasian', 'Eurasian-alpine' and 'North American' groups) within A. muscaria. However, the historical dispersal patterns of A. muscaria remained unclear. In our project, we collected specimens from arctic, boreal and humid temperate regions in Alaska, and generated DNA sequence data from the protein-coding beta-tubulin gene and the internal transcribed spacer (ITS) and large subunit (LSU) regions of the ribosomal DNA repeat. Homologous sequences from additional A. muscaria isolates were downloaded from GenBank. We conducted phylogenetic and nested clade analyses (NCA) to reveal the phylogeographic history of the species complex. Although phylogenetic analyses confirmed the existence of the three above-mentioned clades, representatives of all three groups were found to occur sympatrically in Alaska, suggesting that they represent cryptic phylogenetic species with partially overlapping geographic distributions rather than being allopatric populations. All phylogenetic species share at least two morphological varieties with other species, suggesting ancestral polymorphism in pileus and wart colour pre-dating their speciations. The ancestral population of A. muscaria likely evolved in the Siberian-Beringian region and underwent fragmentation as inferred from NCA and the coalescent analyses. The data suggest that these populations later evolved into species, expanded their range in North America and Eurasia. In addition to range expansions, populations of all three species remained in Beringia and adapted to the cooling climate.     

WHY DOES A METHOD THAT FAILS CONTINUE TO BE USED?

As a critical framework for addressing a diversity of evolutionary and ecological questions, any method that provides accurate and detailed phylogeographic inference would be embraced. What is difficult to understand is the continued use of a method that not only fails, but also has never been shown to work--nested clade analysis is applied widely even though the conditions under which the method will provide reliable results have not yet been demonstrated. This contradiction between performance and popularity is even more perplexing given the recent methodological and computational advances for making historical inferences, which include estimating population genetic parameters and testing different biogeographic scenarios. Here I briefly review the history of criticisms and rebuttals that focus specifically on the high rate of incorrect phylogeographic inference of nested-clade analysis, with the goal of understanding what drives its unfettered popularity. In this case, the appeal of what nested-clade analysis claims to do--not what the method actually achieves--appears to explain its paradoxical status as a favorite method that fails. What a method promises, as opposed to how it performs, must be considered separately when evaluating whether the method represents a valuable tool for historical inference.     

Statistical hypothesis testing in intraspecific phylogeography: nested clade phylogeographical analysis vs. approximate Bayesian computation

Nested clade phylogeographical analysis (NCPA) and approximate Bayesian computation (ABC) have been used to test phylogeographical hypotheses. Multilocus NCPA tests null hypotheses, whereas ABC discriminates among a finite set of alternatives. The interpretive criteria of NCPA are explicit and allow complex models to be built from simple components. The interpretive criteria of ABC are ad hoc and require the specification of a complete phylogeographical model. The conclusions from ABC are often influenced by implicit assumptions arising from the many parameters needed to specify a complex model. These complex models confound many assumptions so that biological interpretations are difficult. Sampling error is accounted for in NCPA, but ABC ignores important sources of sampling error that creates pseudo-statistical power. NCPA generates the full sampling distribution of its statistics, but ABC only yields local probabilities, which in turn make it impossible to distinguish between a good fitting model, a non-informative model, and an over-determined model. Both NCPA and ABC use approximations, but convergences of the approximations used in NCPA are well defined whereas those in ABC are not. NCPA can analyse a large number of locations, but ABC cannot. Finally, the dimensionality of tested hypothesis is known in NCPA, but not for ABC. As a consequence, the 'probabilities' generated by ABC are not true probabilities and are statistically non-interpretable. Accordingly, ABC should not be used for hypothesis testing, but simulation approaches are valuable when used in conjunction with NCPA or other methods that do not rely on highly parameterized models.     

Phylogeographical patterns in chloroplast DNA variation within the Acacia acuminata (Leguminosae: Mimosoideae) complex in Western Australia

Abstract The Acacia acuminata complex includes three taxa, A. acuminata ssp. acuminata, A. acuminata ssp. burkittii and A. oldfieldii, along with several informal variants of A. acuminata. It is widespread throughout southern Australia with the centre of diversity in south‐west Western Australia. Phylogeographical patterns in the complex were investigated using a nested clade analysis of cpDNA RFLPs from 25 populations in Western Australia. Except for A. oldfieldii that was clearly identified as a distinct entity, haplotypes were not restricted to sub‐specific taxa or variants within A. acuminata. There was significant association between phylogenetic position of many haplotypes and their geographical distribution. The fine‐scale phylogeographical patterns were complex but at deeper levels in the phylogeny there was evidence of divergence between two lineages. The pattern of shared haplotypes between lineages suggests retention of ancestral polymorphism as a result of incomplete lineage sorting. The divergence of these lineages is consistent with fragmentation caused by climatic instability during the Pleistocene.     

Impact of Plio-Pleistocene arid cycling on the population history of a southwestern Australian frog

Southwestern Australia is regarded as a global biodiversity hotspot. The region contains a high number of endemic species, ranging from Gondwanan relicts to much more recently evolved plant and animal species. Myobatrachid frogs are diverse in southwestern Australia, and while we know they have speciated in situ in the southwest, we know little about the temporal and geographical patterning of speciation events. Crinia georgiana is an ideal subject to test hypotheses concerning the effect of climatic history on southwestern Australian anurans, as it is an old lineage with a broad distribution covering the entire region. We compiled an extensive phylogeographical data set based on 1085 bp of the mitochondrial gene ND2 for 68 individuals from 18 sites across the species' range. Two major genetic clades were identified which were largely confined to the high rainfall and southeast coastal biogeographical zones, respectively. The clades appear to have diverged around the Plio-Pleistocene border (1.26-1.72 million years ago), concordant with increasing intensity and frequency of arid climate cycles. Subsequent phylogeographical structure appears to have developed primarily during the Pleistocene climatic fluctuations that also have been integral in generating species diversity in the endemic southwestern Australian flora. Phylogeographical analyses identified several dispersal routes, possible refugial areas within the range of the species and also regions of secondary contact. Dispersal routes identified may now be closed to the species because of habitat destruction and salinity problems in inland regions, posing concerns about the evolutionary potential of the species in light of predicted climate change.     

DetSel 1.0: a computer program to detect markers responding to selection

Estimating population parameters from polymorphism frequency data requires neutral genetic markers. Any departure from neutrality may invalidate the inferences drawn from such analyses. We recently discussed the possibility of identifying markers that show deviation from neutral expectations in pairwise comparisons of diverging populations. We are now releasing a user-friendly software package that implements all the necessary steps to identify the signature of selection among molecular markers in a set of polymorphism data. This software can be downloaded free of charge at http://www.univ-montp2.fr/~genetix/detsel/detsel.html.     

Phylogeography of subterranean and surface populations of water lice Asellus aquaticus (Crustacea: Isopoda)

The water louse Asellus aquaticus is a widespread, euryoecious species, mostly uniform throughout its range. However, six subspecies are known from the Dinaric karst in the northwestern Balkans. They include some specialized subterranean populations. The pattern of genetic variation among subterranean and surface populations in this hydrographically highly fragmented karst region was investigated using a 653 bp fragment of the mitochondrial gene (COI). Sequencing of 168 individuals from 25 localities revealed 72 haplotypes. amova and methods of phylogenetic reconstruction all uncovered hydrographic structuring of genetic variation of the populations. Nested clade analysis pointed out several fragmentation events, along with some range expansions within hydrographical systems. By superimposing the subterranean mode of life on the phylogeographical pattern, three independent cave colonizations could be inferred within a distance of < 100 km. Caves were invaded after the ancestral surface populations became isolated through vicariant fragmentation. A possible scenario of hydrographic history of the region was constructed combining the molecular data with palaeogeographical information.