Mito‐nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW–NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW–NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito‐nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within‐species discord. Male‐mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.     

Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.)

Many species have mitochondrial DNA lineages that are phylogenetically intermixed with other species, but studies have rarely tested the cause of such paraphyly. In this study, we tested two hypotheses that could explain mitochondrial paraphyly of Holarctic gadwalls (Anas strepera) with respect to Asian falcated ducks (A. falcata). First, hybridization could have resulted in falcated duck mitochondrial DNA (mtDNA) introgressing into the gadwall gene pool. Second, gadwalls and falcated ducks could have diverged so recently that mtDNA lineages have not sorted to reciprocal monophyly. We used coalescent analyses of three independent loci to distinguish between these two hypotheses. Two lines of evidence support introgression. First, analyses of the three loci combined show that some introgression is necessary to explain current genetic diversity in gadwalls. Second, we generated alternative predictions regarding time since divergence estimated from mtDNA: falcated ducks and gadwalls would have diverged between 65,000 and 700,000 years before present (ybp) under the introgression hypothesis and between 11,000 and 76,000 ybp under the incomplete lineage sorting hypothesis. The two independent nuclear introns indicated that these species diverged between 210,000 and 5,200,000 ybp, which did not overlap the predicted time for incomplete lineage sorting. These analyses also suggested that ancient introgression ( approximately 14,000 ybp) has resulted in the widespread distribution and high frequency of falcated-like mtDNA (5.5% of haplotypes) in North America. This is the first study to use a rigorous quantitative framework to reject incomplete lineage sorting as the cause of mitochondrial paraphyly.     

How and when did Old World ratsnakes disperse into the New World?

To examine Holarctic snake dispersal, we inferred a phylogenetic tree from four mtDNA genes and one scnDNA gene for most species of the Old World (OW) and New World (NW) colubrid group known as ratsnakes. Ancestral area distributions are estimated for various clades using divergence-vicariance analysis and maximum likelihood on trees produced using Bayesian inference. Dates of divergence for the same clades are estimated using penalized likelihood with statistically crosschecked calibration references obtained from the Miocene fossil record. With ancestral areas and associated dates estimated, various hypotheses concerning the age and environment associated with the origin of ratsnakes and the dispersal of NW taxa from OW ancestors were tested. Results suggest that the ratsnakes originated in tropical Asia in the late Eocene and subsequently dispersed to the Western and Eastern Palearctic by the early Oligocene. These analyses also suggest that the monophyletic NW ratsnakes (the Lampropeltini) diverged from OW ratsnakes and dispersed through Beringia in the late Oligocene/early Miocene when this land bridge was mostly composed of deciduous and coniferous forests.     

Range‐wide multilocus phylogeography of the red fox reveals ancient continental divergence,minimal genomic exchange and distinct demographic histories

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization time frame suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y‐chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long‐term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.     

Genetic similarity among Eurasian subspecies of boreal owls Aegolius funereus

Boreal owls Aegolius funereus (referred to as Tengmalm's owls in Europe) breed in boreal forests throughout the Holarctic region and in high-elevation subalpine forests further south. They are currently classified as seven subspecies; six found throughout Eurasia, and one in North America. The geographic distribution of boreal owls in North America and Eurasia is similar, as are their patterns of dispersal and irruption. Because a recent genetic study of boreal owls in North America found very little genetic differentiation among widely disparate locations, we expected that boreal owls in Eurasia similarly would have very little genetic differentiation. Using seven microsatellite markers, we analyzed genetic samples from 275 boreal owls in North America, 36 in Norway, and five in eastern Russia. We found no detectable genetic differentiation between Norwegian and Russian owls, but notable differentiation between North American and Eurasian owls. Low intra-continental genetic differentiation likely results from high rates of long-distance dispersal among subpopulations of boreal owls. In light of these results, we recommend further genetic sampling of boreal owls throughout Eurasia in order to determine whether six separate subspecies here are warranted.     

Genetic signatures of intermediate divergence: population history of Old and New World Holarctic ravens (Corvus corax)

Many studies of phylogeography, speciation, and species limits restrict their focus to a narrow issue: gene tree monophyly. However, reciprocal monophyly does not provide an ideal touchstone criterion of any aspect of evolutionary divergence. There is a continuum of divergence stages as isolated populations go from initial allele frequency differences to well-differentiated species. Studying intermediate stages of divergence will increase our understanding of geographical speciation, species limits, and conservation priorities. We develop a conceptual framework and terminology for thinking about the stages of 'intermediate polyphyly'. The Holarctic clade of common ravens (Corvus corax), found throughout much of Eurasia and North America, provides a case study of these stages of intermediate divergence. We used coalescent, phylogenetic, and population genetic methods to investigate the history and current status of this Old World-New World distribution using 107 mitochondrial control region sequences. Phylogenetically, New World and Old World samples are intermixed. However, most samples are grouped into small subclades that are restricted to either the New World or the Old World, and only one haplotype is shared between the hemispheres. Analysis of moleculalr variance (amova) results reflect this low haplotype sharing between hemispheres (Phi(ST) = 0.13, P < 0.01). Isolation with Migration (im) coalescent results suggest a sustained period of divergence between the hemispheres and low levels of maternal gene flow. Although there has not been sufficient time to evolve reciprocal monophyly and some gene flow may occur, New World and Old World ravens are genetically quite distinct. We use this example to demonstrate these early stages of divergence as populations go from sharing only internal haplotypes, to sharing no haplotypes, to having population specific subclades. Studies of phylogeography, speciation and systematics will benefit from increased attention to these stages of intermediate polyphyly.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Parasites, fossils and geologic history: Historical biogeography of the South American freshwater croakers, Plagioscion spp. (Teleostei, Sciaenidae)

A hypothesis on the historical biogeography of the freshwater croakers, Plagioscion spp. (Teleostei: Sciaenidae), is developed based on data from the phylogeny of their monogenoidean parasites (Platyhelminthes), the geology of South America and the fossil record. Analyses suggest that the common ancestor of Plagioscion spp. colonized freshwater concomitantly with the common ancestor of their parasites, Euryhaliotrema spp. Colonization probably occurred via a marine transgression through western Venezuela that developed about 20 million years ago (Mya). This transgression with its postulated highly variable temporal and spatial salinity conditions most likely represented the facilitating event associated with freshwater colonization. A less likely alternative, that colonization occurred via the Mar del Plata in southern South America, is not supported by the geographical distributions and putative phylogeny of extant species of Plagioscion .     

New Guinea highland origin of a widespread arthropod supertramp

The biologically and geologically extremely diverse archipelagos of Wallacea, Australasia and Oceania have long stimulated ecologists and evolutionary biologists. Yet, few molecular phylogenetic analyses of the terrestrial fauna have been carried out to understand the evolutionary patterns. We use dense taxon and character sampling of more than 7000 bp DNA sequence data for a group of diving beetles ranging from the Holarctic throughout Asia to as far east as French Polynesia. We here show that an ecologically diverse, common and widespread (Portugal to New Zealand) arthropod supertramp species originated in the highlands of New Guinea, ca 6.0–2.7 Myr ago. The approximately 25 closely related species are narrow endemics in Australasia/Oceania. The ancestor of this clade colonized that region from Eurasia ca 9–7 Myr ago. Our finding contradicts the widely held view of local endemism as an evolutionary dead end, as we find multiple peripatric speciation events within the Pleistocene and complex colonization patterns between the Oriental and Australian zoogeographic regions, including the recolonization of Eurasia, jumping across Wallace''s line and colonization of continental Australia out of New Guinea. Our study strongly highlights the importance of dispersal over water gaps in shaping biogeographic patterns.     

Invasion history of North American Canada thistle,Cirsium arvense

Aim Canada thistle (Cirsium arvense– Cardueae, Asteraceae) is one of the worst invasive plants world‐wide. Native to Eurasia, its unintentional introduction into North America now threatens the native flora and is responsible for enormous agricultural losses. The goals of this study are to: (1) reconstruct the evolutionary history of C. arvense and estimate how often it may have colonized North America, (2) compare the genetic diversity between European and North American populations to detect signs of demographic bottlenecks and/or patterns of population admixture, and (3) conduct bioclimatic comparisons to infer eventual niche shifts following this species’ introduction into North America. Location Europe and North America. Methods A total of 1522 individuals from 58 populations were investigated with six microsatellite markers. Estimates of heterozygosity (H_E) and allelic richness (R_S) were quantified for each population, and population structure was inferred via analyses of molecular variance (AMOVAs), principal components analyses (PCAs), Mantel tests and Bayesian clustering analyses. Climatic niche spaces were based on 19 bioclimatic variables extracted from approximately 32,000 locations covering the entire range, and compared using PCA and hierarchical cluster analysis. Results Although there is evidence of multiple introductions from divergent European lineages, North American populations of C. arvense exhibited significantly lower levels of genetic diversity than their putative ancestors. Bioclimatic comparisons pointed to a high degree of niche conservatism during invasion, but indicated that genotypes from the former USSR and Central European mountain chains were probably best adapted to invade North America upon entry into the continent. Main conclusions Genetic and historical data suggest that C. arvense first entered North America from Western Europe with the first European settlers, and was later introduced from Eastern Europe into the prairie states during the agricultural boom. The species went through a significant bottleneck following its introduction into the New World, but the level of genetic diversity remained high owing to admixture between genetically differentiated lineages and to a highly efficient outcrossing breeding system.     

Realized vs apparent reduction in enemies of the European starling

Release from parasites, pathogens or predators (i.e. enemies) is a widely cited ‘rule of thumb’ to explain the proliferation of nonindigenous species in their introduced regions (i.e. the ‘enemy release hypothesis’, or ERH). Indeed, profound effects of some parasites and predators on host populations are well documented. However, some support for the ERH comes from studies that find a reduction in the species richness of enemies in the introduced range, relative to the native range, of particular hosts. For example, data on helminth parasites of the European starling in both its native Eurasia and in North America support a reduction of parasites in the latter. However, North American ‘founder’ starlings were likely not chosen randomly from across Eurasia. This could result in an overestimation of enemy release since enemies affect their hosts on a population level. We control for the effects of subsampling colonists and find, contrary to previous reports, no evidence that introduced populations of starlings experienced a reduction in the species richness of helminth parasites after colonization of North America. These results highlight the importance of choosing appropriate contrast groups in biogeographical analyses of biological invasions to minimize the confounding effects of ‘propagule biases’.     

Refugial origins of reindeer (Rangifer tarandus L.) inferred from mitochondrial DNA sequences

The glacial-interglacial cycles of the upper Pleistocene have had a major impact on the recent evolutionary history of Arctic species. To assess the effects of these large-scale climatic fluctuations to a large, migratory Arctic mammal, we assessed the phylogeography of reindeer (Rangifer tarandus) as inferred from mitochondrial DNA (mtDNA) sequence variation in the control region. Phylogenetic relationships among haplotypes seem to reflect historical patterns of fragmentation and colonization rather than clear-cut relationships among extant populations and subspecies. Three major haplogroups were detected, presumably representing three separate populations during the last glacial. The most influential one has contributed to the gene pool of all extant subspecies and seems to represent a large and continuous glacial population extending from Beringia and far into Eurasia. A smaller, more localized refugium was most likely isolated in connection with ice expansion in western Eurasia. A third glacial refugium was presumably located south of the ice sheet in North America, possibly comprising several separate refugial populations. Significant demographic population expansion was detected for the two haplogroups representing the western Eurasian and Beringian glacial populations. The former apparently expanded when the ice cap retreated by the end of the last glacial. The large continuous one, in contrast, seems to have expanded by the end of the last interglacial, indicating that the warm interglacial climate accompanied by marine transgression and forest expansion significantly confined population size on the continental mainland. Our data demonstrate that the current subspecies designation does not reflect the mtDNA phylogeography of the species, which in turn may indicate that morphological differences among subspecies have evolved as adaptive responses to postglacial environmental change.     

Phylogeographical lineages of Arctic grayling (Thymallus arcticus) in North America: divergence, origins and affinities with Eurasian Thymallus

The number and location of Arctic glacial refugia utilized by taxa during the Pleistocene are continuing uncertainties in Holarctic phylogeography. Arctic grayling (Thymallus arcticus) are widely distributed in freshwaters from the eastern side of Hudson Bay (Canada) west to central Asia. We studied mitochondrial DNA (mtDNA) and microsatellite DNA variation in North American T. arcticus to test for genetic signatures of survival in, and postglacial dispersal from, multiple glacial refugia, and to assess their evolutionary affinities with Eurasian Thymallus. In samples from 32 localities, we resolved 12 mtDNA haplotypes belonging to three assemblages that differed from each other in sequence by between 0.75 and 2.13%: a 'South Beringia' lineage found from western Alaska to northern British Columbia, Canada; a 'North Beringia' lineage found on the north slope of Alaska, the lower Mackenzie River, and to eastern Saskatchewan; and a 'Nahanni' lineage confined to the Nahanni River area of the upper Mackenzie River drainage. Sequence analysis of a portion of the control region indicated monophyly of all North American T. arcticus and their probable origin from eastern Siberian T. arcticus at least 3 Mya. Arctic grayling sampled from 25 localities displayed low allelic diversity and expected heterozygosity (H(E)) across five microsatellite loci (means of 2.1 alleles and 0.27 H(E), respectively) and there were declines in these measures of genetic diversity with distance eastward from the lower Yukon River Valley. Assemblages defined by mtDNA divergences were less apparent at microsatellite loci, but again the Nahanni lineage was the most distinctive. Analysis of molecular variance indicated that between 24% (microsatellite DNA) and 81% (mtDNA) of the variance was attributable to differences among South Beringia, North Beringia and Nahanni lineages. Our data suggest that extant North American Arctic grayling are more diverse phylogeographically than previously suspected and that they consist of at least three major lineages that originated in distinct Pleistocene glacial refugia. T. arcticus probably originated and dispersed from Eurasia to North America in the late to mid-Pliocene, but our data also suggest more recent (mid-late Pleistocene) interactions between lineages across Beringia.     

Contemporary richness of holarctic trees and the historical pattern of glacial retreat

The length of time land has been available for colonization by plants and other organisms could provide a partial explanation of the contemporary richness gradients of trees. According to this hypothesis, increasing times of land availability entail higher chances of recolonization, which eventually have positive effects on tree richness. To test this, we generated a dataset of the Holarctic trees and evaluated the influence of cell age, a measure of the time since an area became free of ice, on the observed tree richness gradients. We found that cell age is associated with richness in both Europe and North America, after controlling for contemporary climate patterns, suggesting that the historical pattern of glacial retreat in response to post-Pleistocene global warming has left a signal still detectable after at least 14 000 yr. The results were consistent using a range of modelling approaches or whether Europe and North America were analyzed separately or in concert. We conclude that, although secondary to contemporary climate, the post-glacial recolonization hypothesis is broadly supported at temperate latitudes.     

Mitochondrial population genomics supports a single pre-Clovis origin with a coastal route for the peopling of the Americas

It is well accepted that the Americas were the last continents reached by modern humans, most likely through Beringia. However, the precise time and mode of the colonization of the New World remain hotly disputed issues. Native American populations exhibit almost exclusively five mitochondrial DNA (mtDNA) haplogroups (A-D and X). Haplogroups A-D are also frequent in Asia, suggesting a northeastern Asian origin of these lineages. However, the differential pattern of distribution and frequency of haplogroup X led some to suggest that it may represent an independent migration to the Americas. Here we show, by using 86 complete mitochondrial genomes, that all Native American haplogroups, including haplogroup X, were part of a single founding population, thereby refuting multiple-migration models. A detailed demographic history of the mtDNA sequences estimated with a Bayesian coalescent method indicates a complex model for the peopling of the Americas, in which the initial differentiation from Asian populations ended with a moderate bottleneck in Beringia during the last glacial maximum (LGM), around approximately 23,000 to approximately 19,000 years ago. Toward the end of the LGM, a strong population expansion started approximately 18,000 and finished approximately 15,000 years ago. These results support a pre-Clovis occupation of the New World, suggesting a rapid settlement of the continent along a Pacific coastal route.     

Barn swallows before barns: population histories and intercontinental colonization

The barn swallow (Hirundo rustica) is one of most widely distributed swallows, owing in part to its recent switch from natural nest sites to human structures. We conducted phylogenetic analysis of mitochondrial (mt) and nuclear DNA to explore the recent evolutionary history of this species. Strongly supported mtDNA clades corresponded to Europe, Asia and North America plus the Baikal region of Asia. Analysis of sequence data from a sex-linked nuclear gene was unable to recover the phylogenetic splits in the mtDNA tree, confirming that the main clades evolved recently. The phylogenetic pattern suggests that the ancestral area of the barn swallow was the holarctic; most divergence events are consistent with vicariance. Most unexpectedly, analyses show that barn swallows from North America colonized the Baikal region in the recent past (one fixed substitution). This dispersal direction is opposite of that for most nearctic-palearctic taxon exchanges. Although this invasion was envisioned to coincide with the appearance of new types of human dwelling in the Baikal region, calibration of molecular divergence suggests an older dispersal event. A recent history of gene flow within the main palearctic clades is consistent with range and population expansion owing to new nesting opportunities provided by human settlements. Contrary to expectation, populations in North America appear historically larger and more stable than those in the palearctic. The Baikal population apparently has not increased greatly since colonization.     

Microsatellites reveal origin and genetic diversity of Eurasian invasions by one of the world’s most notorious marine invader,Mnemiopsis leidyi (Ctenophora)

Marine invasions are taking place at an increasing rate. When occurring in blooms, zooplanktivorous comb jellies of the genus Mnemiopsis are able to cause pelagic regime shifts in coastal areas and may cause the collapse of commercially important fish populations. Using microsatellites, developed for the first time in the phylum Ctenophora, we show that Mnemiopsis leidyi has colonized Eurasia from two source regions. Our preliminary data set included four sites within the putative source region (US East Coast and Gulf of Mexico) and 10 invaded locations in Eurasian waters. Bayesian clustering and phylogeographic approaches revealed the origin of earlier invasions of the Black and Caspian Sea in the 1980s/1990s within or close to the Gulf of Mexico, while the 2006 invasion of the North and Baltic Seas can be directly traced to New England (pairwise F_ST = 0). We found no evidence for mixing among both gene pools in the invaded areas. While the genetic diversity (allelic richness) remained similar in the Baltic Sea compared to the source region New England, it was reduced in the North Sea, supporting the view of an initial invasion of Northern Europe to a Baltic Sea port. In Black and Caspian Sea samples, we found a gradual decline in allelic richness compared to the Gulf of Mexico region, supporting a stepping‐stone model of colonization with two sequential genetic founder events. Our data also suggest that current practices of ballast water treatment are insufficient to prevent repeated invasions of gelatinous zooplankton.     

Genetic structure and the North American postglacial expansion of the barnacle, Semibalanus balanoides

Population genetic characteristics are shaped by the life-history traits of organisms and the geologic history of their habitat. This study provides a neutral framework for understanding the population dynamics and opportunities for selection in Semibalanus balanoides, a species that figures prominently in ecological and evolutionary studies in the Atlantic intertidal. We used mitochondrial DNA (mtDNA) control region (N = 131) and microsatellite markers (～40 individuals/site/locus) to survey populations of the broadly dispersing acorn barnacle from 8 sites spanning 800 km of North American coast and 1 site in Europe. Patterns of mtDNA sequence evolution were consistent with larger population sizes in Europe and population expansion at the conclusion of the last ice age, approximately 20?000 years ago, in North America. A significant portion of mitochondrial diversity was partitioned between the continents (?(ST) = 0.281), but there was only weak structure observed from mtDNA within North America. Microsatellites showed significant structuring between the continents (F(ST) = 0.021) as well as within North America (F(ST) = 0.013). Isolation by distance in North America was largely driven by a split between populations south of Cape Cod and all others (P < 10(-4)). The glacial events responsible for generating allelic diversity at mtDNA and microsatellites may also be responsible for generating selectable variation at metabolic enzymes in S. balanoides.     

Historical biogeography, divergence times, and diversification patterns of bumble bees (Hymenoptera: Apidae: Bombus)

Bumble bees (Bombus) are a cold-adapted, largely alpine group that can elucidate patterns of Holarctic historical biogeography, particularly in comparison to the alpine plants with which they likely coevolved. A recently published molecular phylogeny of bumble bees provides uniquely comprehensive species sampling for exploring historical patterns of distribution and diversification. Using this phylogeny and detailed data on extant distributions, I reconstruct the historical distribution of bumble bees in a temporal framework, estimating divergence times using fossil data and molecular rates derived from the literature. The nearly comprehensive phylogeny allows assessment of the tempo of diversification within the bumble bees using lineage-through-time plots and diversification statistics, which have been performed with special consideration to confidence intervals. These analyses reveal movements of Bombus concordant with geographic and climatic events of the late Cenozoic. The initial diversification of extant bumble bee lineages was estimated at around 25 to 40 Ma, near the Eocene-Oligocene boundary 34 Ma, a period of dramatic global cooling. Dispersal-vicariance analysis (DIVA) predicted an Old World Bombus ancestor, with early diversification events largely restricted to the eastern Old World. The numerous intercontinental dispersal events occurred mostly in the direction of Old World to New World and North America to South America. Early movements from the Palearctic into the Nearctic most likely took place after 20 Ma and may have coincided with a period of Miocene cooling that gave rise to taiga habitat across Beringia. Subsequent dispersal between these regions is estimated to have occurred among boreal and tundra-adapted species mostly in the last 5 million years. Radiations are estimated in both Nearctic and Neotropical regions at approximately 6 to 8 Ma and after 3.5 Ma, concordant with the opening of land corridors between the continents.