Comparative Phylogeography Highlights the Double-Edged Sword of Climate Change Faced by Arctic- and Alpine-Adapted Mammals

Recent studies suggest that alpine and arctic organisms may have distinctly different phylogeographic histories from temperate or tropical taxa, with recent range contraction into interglacial refugia as opposed to post-glacial expansion out of refugia. We use a combination of phylogeographic inference, demographic reconstructions, and hierarchical Approximate Bayesian Computation to test for phylodemographic concordance among five species of alpine-adapted small mammals in eastern Beringia. These species (Collared Pikas, Hoary Marmots, Brown Lemmings, Arctic Ground Squirrels, and Singing Voles) vary in specificity to alpine and boreal-tundra habitat but share commonalities (e.g., cold tolerance and nunatak survival) that might result in concordant responses to Pleistocene glaciations. All five species contain a similar phylogeographic disjunction separating eastern and Beringian lineages, which we show to be the result of simultaneous divergence. Genetic diversity is similar within each haplogroup for each species, and there is no support for a post-Pleistocene population expansion in eastern lineages relative to those from Beringia. Bayesian skyline plots for four of the five species do not support Pleistocene population contraction. Brown Lemmings show evidence of late Quaternary demographic expansion without subsequent population decline. The Wrangell-St. Elias region of eastern Alaska appears to be an important zone of recent secondary contact for nearctic alpine mammals. Despite differences in natural history and ecology, similar phylogeographic histories are supported for all species, suggesting that these, and likely other, alpine- and arctic-adapted taxa are already experiencing population and/or range declines that are likely to synergistically accelerate in the face of rapid climate change. Climate change may therefore be acting as a double-edged sword that erodes genetic diversity within populations but promotes divergence and the generation of biodiversity.     

The ‘Expansion–Contraction’ model of Pleistocene biogeography: rocky shores suffer a sea change?

Approximately 20 000 years ago the last glacial maximum (LGM) radically altered the distributions of many Northern Hemisphere terrestrial organisms. Fewer studies describing the biogeographic responses of marine species to the LGM have been conducted, but existing genetic data from coastal marine species indicate that fewer taxa show clear signatures of post-LGM recolonization. We have assembled a mitochondrial DNA (mtDNA) data set for 14 co-distributed northeastern Pacific rocky-shore species from four phyla by combining new sequences from ten species with previously published sequences from eight species. Nuclear sequences from four species were retrieved from GenBank, plus we gathered new elongation factor 1-α sequences from the barnacle Balanus glandula . Results from demographic analyses of mtDNA for five (36%) species ( Evasterias troschelii, Pisaster ochraceus, Littorina sitkana, L. scutulata, Xiphister mucosus ) were consistent with large population expansions occurring near the LGM, a pattern expected if these species recently recolonized the region. However, seven (50%) species ( Mytilus trossulus, M. californianus, B. glandula, S. cariosus, Patiria miniata, Katharina tunicata , X. atropurpureus ) exhibited histories consistent with long-term stability in effective population size, a pattern indicative of regional persistence during the LGM. Two species of Nucella with significant mtDNA genetic structure showed spatially variable demographic histories. Multilocus analyses for five species were largely consistent with mtDNA: the majority of multilocus interpopulation divergence times significantly exceeded the LGM. Our results indicate that the LGM did not extirpate the majority of species in the northeastern Pacific; instead, regional persistence during the LGM appears a common biogeographic history for rocky-shore organisms in this region.     

Utility of DNA viruses for studying human host history: case study of JC virus

Microbial pathogens, and viruses in particular, can serve as important complements to traditional genetic markers when investigating the population histories of their human host. The range of mutation rates for DNA viruses suggests that DNA viruses can be useful markers of both recent and ancient events in their host histories. Here, we assess the utility of a well known DNA virus, JC virus (JCV), for investigating human history and demography. Using complete coding viral genomes, we confirm the phylogeographic structure of JCV in populations worldwide and provide coalescent estimates of its evolutionary rate under two alternative models of its history. Using these rate estimates, we compare Bayesian skyline plots of population size changes for JCV to those of its human host as estimated with coding mitochondrial genomes of the latter. These comparisons, when combined with other evidence including a log Bayes Factor model test, show that JCV is evolving rapidly and is therefore tracking the recent history of its human host. These results support the hypothesis that post-World War II societal changes are most likely responsible for the recent demographic patterns observed among different regional JCV populations. In sum, fast evolving DNA viruses, such as JCV, can complement RNA viruses to provide novel insights about the recent history and demography of their human host.     

Do complex population histories drive higher estimates of substitution rate in phylogenetic reconstructions?

Our curiosity about biodiversity compels us to reconstruct the evolutionary past of species. Molecular evolutionary theory now allows parameterization of mathematically sophisticated and detailed models of DNA evolution, which have resulted in a wealth of phylogenetic histories. But reconstructing how species and population histories have played out is critically dependent on the assumptions we make, such as the clock-like accumulation of genetic differences over time and the rate of accumulation of such differences. An important stumbling block in the reconstruction of evolutionary history has been the discordance in estimates of substitution rate between phylogenetic and pedigree-based studies. Ancient genetic data recovered directly from the past are intermediate in time scale between phylogenetics-based and pedigree-based calibrations of substitution rate. Recent analyses of such ancient genetic data suggest that substitution rates are closer to the higher, pedigree-based estimates. In this issue, Navascués & Emerson (2009) model genetic data from contemporary and ancient populations that deviate from a simple demographic history (including changes in population size and structure) using serial coalescent simulations. Furthermore, they show that when these data are used for calibration, we are likely to arrive at upwardly biased estimates of mutation rate.     

Reconstructing the demographic history of orang‐utans using Approximate Bayesian Computation

Investigating how different evolutionary forces have shaped patterns of DNA variation within and among species requires detailed knowledge of their demographic history. Orang‐utans, whose distribution is currently restricted to the South‐East Asian islands of Borneo (Pongo pygmaeus) and Sumatra (Pongo abelii), have likely experienced a complex demographic history, influenced by recurrent changes in climate and sea levels, volcanic activities and anthropogenic pressures. Using the most extensive sample set of wild orang‐utans to date, we employed an Approximate Bayesian Computation (ABC) approach to test the fit of 12 different demographic scenarios to the observed patterns of variation in autosomal, X‐chromosomal, mitochondrial and Y‐chromosomal markers. In the best‐fitting model, Sumatran orang‐utans exhibit a deep split of populations north and south of Lake Toba, probably caused by multiple eruptions of the Toba volcano. In addition, we found signals for a strong decline in all Sumatran populations ~24 ka, probably associated with hunting by human colonizers. In contrast, Bornean orang‐utans experienced a severe bottleneck ~135 ka, followed by a population expansion and substructuring starting ~82 ka, which we link to an expansion from a glacial refugium. We showed that orang‐utans went through drastic changes in population size and connectedness, caused by recurrent contraction and expansion of rainforest habitat during Pleistocene glaciations and probably hunting by early humans. Our findings emphasize the fact that important aspects of the evolutionary past of species with complex demographic histories might remain obscured when applying overly simplified models.     

Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle‐Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X‐linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb‐tRNA‐control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species’ demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre‐dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero‐Maghreb region and an eastern location (Anatolia/Middle East)–that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species’ movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero‐Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies.     

Storage and absorption in the digestive system of carnivorous and herbivorous prickleback fishes (Teleostei: Stichaeidae): ontogenetic, dietary, and phylogenetic effects

The effects of ontogeny, diet, and phylogeny on glycogen storage levels and esterase and alkaline phosphatase activities in four related prickleback fishes were determined in situ using quantitative histochemistry. Of these species, Cebidichthys violaceus and Xiphister mucosus shift from carnivory to herbivory at approximately 45 mm standard length (SL), whereas Xiphister atropurpureus and Anoplarchus purpurescens remain carnivores. Comparisons between small (30-40 mm SL) and larger (60-75 mm SL) wild-caught juveniles showed that glycogen storage levels and alkaline phosphatase activity were unchanged with ontogeny. Comparisons between the larger wild-caught juveniles and juveniles of the same size that had been raised on a high-protein animal diet revealed that glycogen storage level and alkaline phosphatase activity increased in all species in response to this diet. Esterase activity also increased in response to the high-protein animal diet in all four species but increased with ontogeny only in C. violaceus, X. mucosus, and X. atropurpureus, in the xiphisterine clade, and not in A. purpurescens, in the adjacent alectriine clade. Xiphister mucosus and X. atropurpureus showed indistinguishable responses in esterase activity to ontogeny and diet despite their divergent natural diets. Overall, glycogen storage level and alkaline phosphatase activity responded primarily to diet, whereas esterase activity was also influenced by ontogeny and phylogeny and differed between intestinal regions among the species.     

Inferring the phylogeny of disjunct populations of the azure-winged magpie Cyanopica cyanus from mitochondrial control region sequences

The azure-winged magpie (AWM), Cyanopica cyanus, is found in Asia and Iberia. This remarkable disjunct distribution has been variously explained by either the sixteenth-century introduction of birds into Iberia from the Far East, or by the loss of individuals from the central part of their range as a result of Pleistocene glaciations. We have used the mitochondrial control region to undertake a molecular phylogenetic analysis of the AWM, with sequences examined from individuals collected from across the current distribution range and incorporating representatives of all currently defined subspecies. The Western birds are genetically distinct from their Asian congeners and their divergence is basal in the phylogenetic tree. This indicates that the AWM is native to Iberia and not the result of a recent introduction from Asia. In Asia, two major mitochondrial DNA lineages were identified. These correspond to an Inland Asia group and a Pacific Seaboard group, and are separated topographically by the Da Hingan Ling mountains and the Yellow Sea. Molecular clock estimates suggest that these divergences are associated with Pleistocene glaciations. Furthermore, our data do not support the current classification of the AWM into 10 subspecies, as defined based on morphology and geographical distribution.     

ACCOUNTING FOR RATE VARIATION AMONG LINEAGES IN COMPARATIVE DEMOGRAPHIC ANALYSES

Genetic analyses of contemporary populations can be used to estimate the demographic histories of species within an ecological community. Comparison of these demographic histories can shed light on community responses to past climatic events. However, species experience different rates of molecular evolution, and this presents a major obstacle to comparative demographic analyses. We address this problem by using a Bayesian relaxed‐clock method to estimate the relative evolutionary rates of 22 small mammal taxa distributed across northwestern North America. We found that estimates of the relative molecular substitution rate for each taxon were consistent across the range of sampling schemes that we compared. Using three different reference rates, we rescaled the relative rates so that they could be used to estimate absolute evolutionary timescales. Accounting for rate variation among taxa led to temporal shifts in our skyline‐plot estimates of demographic history, highlighting both uniform and idiosyncratic evolutionary responses to directional climate trends for distinct ecological subsets of the small mammal community. Our approach can be used in evolutionary analyses of populations from multiple species, including comparative demographic studies.     

Nuclear vs. plastid data: complex Pleistocene history of a circumpolar key species

To fully understand the contemporary genetic structure of plants, both nuclear and plastid markers are needed. Three chloroplast DNA (cpDNA) lineages, which probably diverged before the major Pleistocene glaciations, have been identified in the circumpolar/circumboreal Vaccinium uliginosum. Here we investigate its nuclear DNA variation using nuclear ribosomal internal transcribed spacer (ITS) sequences, DNA ploidy level measurements and amplified fragment length polymorphisms (AFLPs). We also extend the cpDNA dataset. Two ITS lineages, corresponding to diploids and tetraploids, respectively, were identified. However, both main sequence types apparently occurred in most individual plants but showed ploidy-biased homogenization and possibly reflect paralogy predating the origin of V. uliginosum. The ploidy levels were largely consistent with the cpDNA lineages, suggesting that the initial cpDNA divergence followed early polyploidizations. Five main AFLP groups were identified, consistent with recent glacial refugia in Beringia, western Siberia, the southern European mountains and areas south/east of the Scandinavian and Laurentide ice sheets. Except from the southern European mountains, there has been extensive expansion from all refugia, resulting in several contact zones. Surprisingly, the presumably older ploidy and cpDNA patterns were partly inconsistent with the main AFLP groups and more consistent with AFLP subgroups. A likely major driver causing the inconsistencies is recent nuclear gene flow via unreduced pollen from diploids to tetraploids. This may prevent cytoplasmic introgression and result in overlayed patterns formed by processes dominating at different time scales. The data also suggest more recent polyploidizations, as well as several chloroplast capture events, further complicating this scenario. This study highlights the importance of combining different marker systems to unravel intraspecific histories.     

Human populations show reduced DNA sequence variation at the factor IX locus

Levels and patterns of human DNA sequence variation vary widely among loci. However, some of this variation may be due to the different populations used in different studies. So far, few studies of diverse human populations have compared different genetic loci for the same samples of populations and individuals. Here, we present new polymorphism data from intron 4 of the Factor IX gene (FIX) sequenced in diverse Old World populations. An explicit comparison is made with another X-linked gene, PDHA1, for which the sampling of individuals was very similar. Despite having a similar amount of divergence from chimpanzees, as do other nuclear genes, FIX has comparatively much less DNA sequence variation among humans. Nucleotide diversity at FIX is the lowest among the existing non-Y chromosome nuclear gene datasets and is less than 10% of the diversity found at PDHA1. Estimates of effective population size based on FIX are 8,558, about half of the value obtained for PDHA1, and the time to the most recent common ancestry among human FIX gene copies (282,000 years) is one of the most recent estimates reported for human genes. Analyses presented here suggest a history for the FIX region that includes recent positive directional selection, or background, selection. The general conclusion emerging is that very large variations can exist between the histories of similar genomic regions, even when sampling differences are minimized.     

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

ABSTRACT: BACKGROUND: The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs). RESULTS: To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are reciprocally monophyletic in the mtgenome phylogeny but a H. pileatus individual falls into the H. lar Y chromosome clade. CONCLUSIONS: Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.     

Phylogeographic history of grey wolves in Europe

Background  

While it is generally accepted that patterns of intra-specific genetic differentiation are substantially affected by glacial history, population genetic processes occurring during Pleistocene glaciations are still poorly understood. In this study, we address the question of the genetic consequences of Pleistocene glaciations for European grey wolves. Combining our data with data from published studies, we analysed phylogenetic relationships and geographic distribution of mitochondrial DNA haplotypes for 947 contemporary European wolves. We also compared the contemporary wolf sequences with published sequences of 24 ancient European wolves.     

Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag

Human artificial chromosomes have been used to model requirements for human chromosome segregation and to explore the nature of sequences competent for centromere function. Normal human centromeres require specialized chromatin that consists of alpha satellite DNA complexed with epigenetically modified histones and centromere-specific proteins. While several types of alpha satellite DNA have been used to assemble de novo centromeres in artificial chromosome assays, the extent to which they fully recapitulate normal centromere function has not been explored. Here, we have used two kinds of alpha satellite DNA, DXZ1 (from the X chromosome) and D17Z1 (from chromosome 17), to generate human artificial chromosomes. Although artificial chromosomes are mitotically stable over many months in culture, when we examined their segregation in individual cell divisions using an anaphase assay, artificial chromosomes exhibited more segregation errors than natural human chromosomes (P < 0.001). Naturally occurring, but abnormal small ring chromosomes derived from chromosome 17 and the X chromosome also missegregate more than normal chromosomes, implicating overall chromosome size and/or structure in the fidelity of chromosome segregation. As different artificial chromosomes missegregate over a fivefold range, the data suggest that variable centromeric DNA content and/or epigenetic assembly can influence the mitotic behavior of artificial chromosomes.     

Two gene families clustered in a small region of the Drosophila genome

Three Drosophila genes that are clustered within 8 X 10(3) bases of DNA at the chromosomal region 44D have been identified and mapped, and the gene cluster entirely sequenced. The three genes are 55 to 60% homologous in DNA sequence. One gene contains an intron in its 5'-proximal protein coding sequence while the other two have none at this position; similarly, another gene has an intron in its 3'-proximal protein coding sequence which is not found in the other genes. All three genes are abundantly expressed together in Drosophila first, second, and early third instar larval stages and in adults, but they are not abundantly expressed in either embryonic, late third instar larval, or pupal stages. This gene family lies 11 X 10(3) bases away from another cluster containing four Drosophila larval cuticle protein genes plus a pseudogene. The cuticle genes are all abundantly expressed throughout third instar larval development. Thus, at least seven protein-coding genes and one pseudogene lie within 27 X 10(3) bases of DNA. Moreover, two small gene families can lie adjacent on a chromosome and exhibit different patterns of developmental regulation, even though individual genes within each clustered family are co-ordinately expressed.     

Responses of plants to changes in Qinghai–Tibetan Plateau and glaciations: Evidence from phylogeography of a Sibiraea (Rosaceae) complex

The genetic structure and phylogeographical history of the alpine shrubs Sibiraea angustata (Rosaceae) and Sibiraea laevigata from the Qinghai–Tibetan Plateau (QTP) were investigated to identify alpine plant responses to changes in the QTP and glaciations. Fifty-five populations were analyzed using four chloroplast DNA (cpDNA) regions and (nuclear ribosomal internal transcribed spacer) nrITS sequence data. In all, 21 cpDNA haplotypes and 13 nrITS sequence types were detected. Analyses of the genetic diversity and phylogenetic relationships detected two rarely reported glacial refugia. One was the Yushu–Nangqian area, and the other consisted of the area from the Songpan Plateau to the southeastern margin of the QTP. Sibiraea species populations experienced divergent evolution and founder effects when they recolonized the QTP platform and adjacent high-altitude regions following glaciations. The divergence times of the main lineages and haplotypes were in the range of 1.60–2.58 Ma. The population size of Sibiraea species in the QTP decreased approximately 23-fold during the last 0.12 Ma, indicating that Sibiraea species were significantly affected by environmental changes in the QTP. Therefore, the rapid uplift of the QTP and subsequent glaciations likely played an important role in driving genetic divergence and population size changes of Sibiraea species in the QTP.     

Functionally distinct nucleic acid binding sites for a group I intron encoded RNA maturase/DNA homing endonuclease

A large number of group I introns encode a family of homologous proteins that either promote intron splicing (maturases) or are site-specific DNA endonucleases that function in intron mobility (a process called "homing"). Genetic studies have shown that some of these proteins have both activities, yet how a single protein carries out both functions remains obscure. The similarity between respective DNA-binding sites and the RNA structure near the 5' and 3' splice sites has fueled speculation that such proteins may use analogous interactions to perform both functions. The Aspergillus nidulans mitochondrial COB group I intron encodes a bi-functional protein, I-AniI, that has both RNA maturase and site-specific DNA endonuclease activities in vitro. Here, we show that I-AniI shows distinctive features of the endonuclease family to which it belongs, including highly specific, tight binding and sequential DNA strand cleavage. Competition experiments demonstrate that I-AniI binds the COB intron RNA even in saturating concentrations of its DNA target site substrate, suggesting that the protein has a separate binding site for RNA. In addition, we provide evidence that two different DNA-binding site mutants of I-AniI have little effect on the protein's RNA maturation activity. Since RNA splicing is likely a secondary adaptation of the protein, these observations support a model in which homing endonucleases may have developed maturase function by utilizing a hitherto "non-functional" protein surface.     

Phylogeography of sable (Martes zibellina L. 1758) in the southeast portion of its range based on mitochondrial DNA variation: highlighting the evolutionary history of the sable

To reveal phylogeographic features of sable (Martes zibellina) in the southeast part of its range, we analyzed variability of the mitochondrial DNA (mtDNA) cytochrome b gene, tRNA (Pro), tRNA (Thr) and control region (D-loop) sequences from 78 specimens in populations of the Russian Far East, northeast China, and Mongolia. Our results revealed the presence of 49 different haplotypes split into two major phylogenetic groups—clades A and B, the latter separated into two clades, B1 and B2. Comparative analysis of D-loop haplotypes in populations originating from the southeast (Russian Far East, China and Mongolia) and the west (northern Urals) portions of sable range indicated that all three mtDNA clades were present in different regional groups. However, highest diversity of clade B1 in northeast China and its nearly complete absence from the Urals suggest that the southeast sable range, being a refuge during Pleistocene glacial periods, can be considered the center of genetic diversification and possibly origin of this species. All divergence estimates fall within the Pleistocene suggesting that Quarternary glaciations played an important role in phylogeographic differentiation of sable.     

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

ABSTRACT: BACKGROUND: Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescence-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies. RESULTS: We found that divergence estimates, from both concatenated gene tree and coalescence-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages. CONCLUSIONS: Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.     

Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.