Population genetic structure in migratory sandhill cranes and the role of Pleistocene glaciations

Previous studies of migratory sandhill cranes (Grus canadensis) have made significant progress explaining evolution of this group at the species scale, but have been unsuccessful in explaining the geographically partitioned variation in morphology seen on the population scale. The objectives of this study were to assess the population structure and gene flow patterns among migratory sandhill cranes using microsatellite DNA genotypes and mitochondrial DNA haplotypes of a large sample of individuals across three populations. In particular, we were interested in evaluating the roles of Pleistocene glaciation events and postglaciation gene flow in shaping the present-day population structure. Our results indicate substantial gene flow across regions of the Midcontinental population that are geographically adjacent, suggesting that gene flow for most of the region follows an isolation-by-distance model. Male-mediated gene flow and strong female philopatry may explain the differing patterns of nuclear and mitochondrial variation. Taken in context with precise geographical information on breeding locations, the morphologic and microsatellite DNA variation shows a gradation from the Arctic-nesting subspecies G. c. canadensis to the nonArctic subspecies G. c. tabida. Analogous to other Arctic-nesting birds, it is probable that the population structure seen in Midcontinental sandhill cranes reflects the result of postglacial secondary contact. Our data suggest that subspecies of migratory sandhills experience significant gene flow and therefore do not represent distinct and independent genetic entities.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Phylogeography of Aglais urticae (Lepidoptera) based on DNA sequences of the mitochondrial COI gene and control region

Mitochondrial DNA (mtDNA) sequences of the COI gene and the control region were used to examine the genetic population structure of Aglais urticae L. (Lepidoptera) over its entire geographic range, i.e., the Palaearctic. The phylogenetic relationships within and between A. urticae subspecies were determined and patterns of mtDNA divergence and ecological differentiation were compared. High gene flow together with a recent and sudden population expansion characterise the genetic population structure of this species. No geographically induced differentiation was observed, nor were subspecies identified as separate evolutionary units. The discrepancy between the genetic and ecological variation is most likely due to the slower rate of mtDNA evolution compared to ecological differentiation. The control region proved to be a less useful molecular marker for the population genetics and the phylogenetic reconstruction of closely related taxa in A. urticae than it has for other species. The extreme bias in adenine and thymine content (A+T=90.91%) probably renders this region highly susceptible to homoplasy, resulting in a less informative molecular marker.     

Phylogeography of the marbled crab Pachygrapsus marmoratus (Decapoda,Grapsidae) along part of the African Mediterranean coast reveals genetic homogeneity across the Siculo-Tunisian Strait versus heterogeneity across the Gibraltar Strait

We investigate the influence of previously postulated biogeographic barriers in the Mediterranean Sea on the population genetic structure of a highly dispersive and continuously distributed coastal species. In particular, we examine nuclear and mitochondrial genetic variation in the marbled crab, Pachygrapsus marmoratus, across part of the African Mediterranean coast in order to assess the influence of the Siculo-Tunisian Strait on its population genetic structure. Four polymorphic microsatellite loci were genotyped for 110 individuals, collected from eight locations covering parts of the Algerian, Tunisian and Libyan coasts. In addition, mtDNA corresponding to the Cox1 gene was sequenced for 80 samples. The corresponding results show contrasting patterns of genetic differentiation. While mtDNA results revealed a homogeneous haplotype composition in our study area, microsatellite data depicted genetic differentiation among populations, but not associated with any geographic barrier. This pattern, already recorded for this species from different geographic regions, may hint at the involvement of a complex series of abiotic and biotic factors in determining genetic structure. Demographic history reconstruction, inferred from mtDNA data, supports demographic and spatial expansion for the North African metapopulation dating back to the Mid-Pleistocene and following an historical bottleneck. Comparison of these African mitochondrial sequences with new sequences from a Turkish population and previously published sequences revealed a weak but significant separation of Atlantic and Mediterranean populations across the Gibraltar Strait, which was not recorded in previous studies of this grapsid species.     

Genetic structure in Iberian and Moroccan populations of the globally threatened great bustard Otis tarda: a microsatellite perspective

Patterns of genetic structure and gene flow among populations help us understand population dynamics and properly manage species of concern. Matrilineal mtDNA sequence data have been instrumental in revealing genetic structure at the intraspecific level, but bi‐parentally inherited markers are needed to confirm patterns at the genome level and to assess the potential role of sex‐biased dispersal on gene flow, particularly in species where males are known to be the main dispersing sex. Here we use microsatellite loci to examine patterns of genetic structure across the range of the great bustard in Iberia and Morocco, an area representing 70% of the world population of this globally threatened species. We used population differentiation statistics and Bayesian analysis of population structure to analyse data from 14 microsatellite loci. These data provide greater resolution than mtDNA sequences, and results reveal the existence of three main genetic units corresponding to Morocco, the northeastern part of Spain, and the rest of the Iberian Peninsula. Our results, together with previous mtDNA data, confirm the genetic differentiation of the northern Africa population and the importance of the Strait of Gibraltar as a barrier to gene flow for both males and females, rendering the Moroccan population a separate management unit of high conservation concern.     

Spatial patterns of genetic diversity across European subspecies of the mountain hare, Lepus timidus L

Fossil evidence shows that populations of species that currently inhabit arctic and boreal regions were not isolated in refugia during glacial periods, but instead maintained populations across large areas of central Europe. These species commonly display little reduction in genetic diversity in northern areas of their range, in contrast to many temperate species. The mountain hare currently inhabits both temperate and arctic-boreal regions. We used nuclear microsatellite and mtDNA sequence data to examine population structure and alternate phylogeographic hypotheses for the mountain hare, that is, temperate type (lower genetic diversity in northern areas) and arctic-boreal type (high northern genetic diversity). Both data sets revealed concordant patterns. Highest allelic richness, expected heterozygosity and mtDNA haplotype diversity were identified in the most northerly subspecies, indicating that this species more closely maps to phylogeographic patterns observed in arctic-boreal rather than temperate species. With regard to population structure, the Alpine and Fennoscandian subspecies were most genetically similar (F(ST) approximately 0.1). These subspecies also clustered together on the mtDNA tree and were assigned with highest likelihood to a common Bayesian cluster. This is consistent with fossil evidence for intermediate populations in the central European plain, persisting well into the postglacial period. In contrast, the geographically close Scottish and Irish populations occupied separate Bayesian clusters, distinct clades on the mtDNA maximum likelihood tree and were genetically divergent from each other (F(ST) > 0.4) indicating the influence of genetic drift, long isolation (possibly dating from the late glacial era) and/or separate postglacial colonisation routes.     

Genetic structure of muskrat (Ondatra zibethicus) and its concordance with taxonomy in North America

Extrinsic factors such as physical barriers play an important role in shaping population genetic structure. A reduction in gene flow leading to population structuring may ultimately lead to population divergence. These divergent populations are often considered subspecies. Because genetic differentiation may represent differences between subspecies, patterns of genetic structure should reflect subspecies groupings. In this study, we examine the contemporary population genetic structure of muskrat (n = 331) and assess the relevance of 4 geographically distinct subspecies designations across northern North America using 9 microsatellite loci. We predicted that patterns of gene flow and genetic structure would reflect the described subspecies. We found evidence of genetic differentiation between western and eastern regions, and muskrats from Newfoundland (NF) showed significantly lower genetic diversity than central regions. A strong isolation by distance pattern was also detected within the eastern cluster. Our results did not differentiate Ondatra zibethicus spatulus (northwest) from O. z. albus (central), but they suggest a distinction between O. z. obscurus (NF) and O. z. zibethicus (east). This study highlights the need for more phylogenetic studies in order to better understand intraspecific divergence and the genetic characterization of subspecies.     

Population differentiation within and among Asian elephant (Elephas maximus) populations in southern India

Southern India, one of the last strongholds of the endangered Asian elephant (Elephas maximus), harbours about one-fifth of the global population. We present here the first population genetic study of free-ranging Asian elephants, examining within- and among-population differentiation by analysing mitochondrial DNA (mtDNA) and nuclear microsatellite DNA differentiation across the Nilgiris-Eastern Ghats, Anamalai, and Periyar elephant reserves of southern India. Low mtDNA diversity and 'normal' microsatellite diversity were observed. Surprisingly, the Nilgiri population, which is the world's single largest Asian elephant population, had only one mtDNA haplotype and lower microsatellite diversity than the two other smaller populations examined. There was almost no mtDNA or microsatellite differentiation among localities within the Nilgiris, an area of about 15,000 km2. This suggests extensive gene flow in the past, which is compatible with the home ranges of several hundred square kilometres of elephants in southern India. Conversely, the Nilgiri population is genetically distinct at both mitochondrial and microsatellite markers from the two more southerly populations, Anamalai and Periyar, which in turn are not genetically differentiated from each other. The more southerly populations are separated from the Nilgiris by only a 40-km-wide stretch across a gap in the Western Ghats mountain range. These results variably indicate the importance of population bottlenecks, social organization, and biogeographic barriers in shaping the distribution of genetic variation among Asian elephant populations in southern India.     

Genetic structuring and recent demographic history of red pandas (Ailurus fulgens) inferred from microsatellite and mitochondrial DNA

Clarification of the genetic structure and population history of a species can shed light on the impacts of landscapes, historical climate change and contemporary human activities and thus enables evidence‐based conservation decisions for endangered organisms. The red panda (Ailurus fulgens) is an endangered species distributing at the edge of the Qinghai‐Tibetan Plateau and is currently subject to habitat loss, fragmentation and population decline, thus representing a good model to test the influences of the above‐mentioned factors on a plateau edge species. We combined nine microsatellite loci and 551 bp of mitochondrial control region (mtDNA CR) to explore the genetic structure and demographic history of this species. A total of 123 individuals were sampled from 23 locations across five populations. High levels of genetic variation were identified for both mtDNA and microsatellites. Phylogeographic analyses indicated little geographic structure, suggesting historically wide gene flow. However, microsatellite‐based Bayesian clustering clearly identified three groups (Qionglai‐Liangshan, Xiaoxiangling and Gaoligong‐Tibet). A significant isolation‐by‐distance pattern was detected only after removing Xiaoxiangling. For mtDNA data, there was no statistical support for a historical population expansion or contraction for the whole sample or any population except Xiaoxiangling where a signal of contraction was detected. However, Bayesian simulations of population history using microsatellite data did pinpoint population declines for Qionglai, Xiaoxiangling and Gaoligong, demonstrating significant influences of human activity on demography. The unique history of the Xiaoxiangling population plays a critical role in shaping the genetic structure of this species, and large‐scale habitat loss and fragmentation is hampering gene flow among populations. The implications of our findings for the biogeography of the Qinghai‐Tibetan Plateau, subspecies classification and conservation of red pandas are discussed.     

Phylogeography and genetic ancestry of tigers (Panthera tigris)

Eight traditional subspecies of tiger (Panthera tigris),of which three recently became extinct, are commonly recognized on the basis of geographic isolation and morphological characteristics. To investigate the species' evolutionary history and to establish objective methods for subspecies recognition, voucher specimens of blood, skin, hair, and/or skin biopsies from 134 tigers with verified geographic origins or heritage across the whole distribution range were examined for three molecular markers: (1) 4.0 kb of mitochondrial DNA (mtDNA) sequence; (2) allele variation in the nuclear major histocompatibility complex class II DRB gene; and (3) composite nuclear microsatellite genotypes based on 30 loci. Relatively low genetic variation with mtDNA,DRB,and microsatellite loci was found, but significant population subdivision was nonetheless apparent among five living subspecies. In addition, a distinct partition of the Indochinese subspecies P. t. corbetti in to northern Indochinese and Malayan Peninsula populations was discovered. Population genetic structure would suggest recognition of six taxonomic units or subspecies: (1) Amur tiger P. t. altaica; (2) northern Indochinese tiger P. t. corbetti; (3) South China tiger P. t. amoyensis; (4) Malayan tiger P. t. jacksoni, named for the tiger conservationist Peter Jackson; (5) Sumatran tiger P. t. sumatrae; and (6) Bengal tiger P. t. tigris. The proposed South China tiger lineage is tentative due to limited sampling. The age of the most recent common ancestor for tiger mtDNA was estimated to be 72,000-108,000 y, relatively younger than some other Panthera species. A combination of population expansions, reduced gene flow, and genetic drift following the last genetic diminution, and the recent anthropogenic range contraction, have led to the distinct genetic partitions. These results provide an explicit basis for subspecies recognition and will lead to the improved management and conservation of these recently isolated but distinct geographic populations of tigers.     

Phylogeography, intraspecific structure and sex-biased dispersal of Dall's porpoise, Phocoenoides dalli, revealed by mitochondrial and microsatellite DNA analyses

Mitochondrial DNA (mtDNA) control-region sequences and microsatellite loci length polymorphisms were used to estimate phylogeographical patterns (historical patterns underlying contemporary distribution), intraspecific population structure and gender-biased dispersal of Phocoenoides dalli dalli across its entire range. One-hundred and thirteen animals from several geographical strata were sequenced over 379 bp of mtDNA, resulting in 58 mtDNA haplotypes. Analysis using F(ST) values (based on haplotype frequencies) and phi(ST) values (based on frequencies and genetic distances between haplotypes) yielded statistically significant separation (bootstrap values P < 0.05) among most of the stocks currently used for management purposes. A minimum spanning network of haplotypes showed two very distinctive clusters, differentially occupied by western and eastern populations, with some common widespread haplotypes. This suggests some degree of phyletic radiation from west to east, superimposed on gene flow. Highly male-biased migration was detected for several population comparisons. Nuclear microsatellite DNA markers (119 individuals and six loci) provided additional support for population subdivision and gender-biased dispersal detected in the mtDNA sequences. Analysis using F(ST) values (based on allelic frequencies) yielded statistically significant separation between some, but not all, populations distinguished by mtDNA analysis. R(ST) values (based on frequencies of and genetic distance between alleles) showed no statistically significant subdivision. Again, highly male-biased dispersal was detected for all population comparisons, suggesting, together with morphological and reproductive data, the existence of sexual selection. Our molecular results argue for nine distinct dalli-type populations that should be treated as separate units for management purposes.     

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.     

Genetic structuring of Patagonian toothfish populations in the Southwest Atlantic Ocean: the effect of the Antarctic Polar Front and deep-water troughs as barriers to genetic exchange

The environmental and/or life history factors affecting genetic exchange in marine species with potential for high dispersal are of great interest, not only from an evolutionary standpoint but also with regard to effective management. Previous genetic studies have demonstrated substantial differentiation among populations of the Patagonian toothfish around the Southern Ocean, indicating breakdown of gene flow across large distances between inhabited shelf areas. The present study examined genetic structuring through analysis of microsatellite loci and restriction fragment length polymorphism (RFLP) of the mitochondrial ND2 gene and control region of the toothfish population in the SW Atlantic, allowing examination of the relative effects of the Antarctic Polar Front (APF), deep-water troughs and distance between sites. Mitochondrial DNA (mtDNA) data indicated a sharp genetic division between the Patagonian Shelf/North Scotia Ridge and the Shag Rocks/South Georgia samples, whereas microsatellite data showed much less distinct structuring and an intermediate position of the North Scotia Ridge samples. We suggest these data indicate that the APF, as a barrier to larval dispersal, is the major inhibitor of genetic exchange between toothfish populations, with deep-water troughs and distance between sites contributing to genetic differentiation by inhibiting migration of relatively sedentary adults. We also suggest that differences between mtDNA and nuclear DNA population patterns may reflect either genome population size effects or (putative) male-biased dispersal.     

Geographic Structure of Mitochondrial and Nuclear Gene Polymorphisms in Australian Green Turtle Populations and Male-Biased Gene Flow

The genetic structure of green turtle (Chelonia mydas) rookeries located around the Australian coast was assessed by (1) comparing the structure found within and among geographic regions, (2) comparing microsatellite loci vs. restriction fragment length polymorphism analyses of anonymous single copy nuclear DNA (ascnDNA) loci, and (3) comparing the structure found at nuclear DNA markers to that of previously analyzed mitochondrial (mtDNA) control region sequences. Significant genetic structure was observed over all regions at both sets of nuclear markers, though the microsatellite data provided greater resolution in identifying significant genetic differences in pairwise tests between regions. Inferences about population structure and migration rates from the microsatellite data varied depending on whether statistics were based on the stepwise mutation or infinite allele model, with the latter being more congruent with geography. Estimated rates of gene flow were generally higher than expected for nuclear DNA (nDNA) in comparison to mtDNA, and this difference was most pronounced in comparisons between the northern and southern Great Barrier Reef (GBR). The genetic data combined with results from physical tagging studies indicate that the lack of nuclear gene divergence through the GBR is likely due to the migration of sGBR turtles through the courtship area of the nGBR population, rather than male-biased dispersal. This example highlights the value of combining comparative studies of molecular variation with ecological data to infer population processes.     

A re-evaluation of management units based on gene flow of a rare waterbird in the Americas

The maintenance of gene flow in species that have experienced population contractions and are geographically fragmented is important to the maintenance of genetic variation and evolutionary potential; thus, gene flow is also important to conservation and management of these species. For example, the Reddish Egret (Egretta rufescens) has recovered after severe population reductions during the 19th and 20th centuries, but population numbers remain below historical levels. In this study, we characterized gene flow among management units of the Reddish Egret by using ten nuclear microsatellite markers and part of the mitochondrial (mtDNA) control region from 176 nestlings captured at eight localities in Mexico (Baja California, Chiapas, Tamaulipas, and Yucatan), the USA (Texas, Louisiana, and Florida), and the Bahamas. We found evidence of population structure and that males disperse more often and across longer distances compared with females, which is congruent with previous banding and telemetry data. The maternally inherited mtDNA and biparentally inherited microsatellite data supported slightly different MU models; however, when interpreted together, a four MU model that considered population structure and geographic proximity was most optimal. Namely, MU 1 (Baja California); MU 2 (Chiapas); MU 3 (Yucatan, Tamaulipas, Texas, and Louisiana); and MU 4 (Florida and the Bahamas). Regions outside our sampled localities (e.g., the Greater Antilles and South America) require additional sampling to fully understand gene flow and movement of individuals across the species’ entire range. However, the four MUs we have defined group nesting localities into genetically similar subpopulations, which can guide future management plans.     

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.     

PHYLOGEOGRAPHY OF SPOTTED OWL (STRIX OCCIDENTALIS) POPULATIONS BASED ON MITOCHONDRIAL DNA SEQUENCES: GENE FLOW,GENETIC STRUCTURE,AND A NOVEL BIOGEOGRAPHIC PATTERN

Mitochondrial DNA control region sequences of spotted owls (Strix occidentalis) allowed us to investigate gene flow, genetic structure, and biogeographic relationships among these forest-dwelling birds of western North America Estimates of gene flow based on genetic partitioning and the phylogeography of haplotypes indicate substantial dispersal within three long-recognized subspecies. However, patterns of individual phyletic relationships indicate a historical absence of gene flow among the subspecies, which are essentially monophyletic. The pattern of haplotype coalescence enabled us to identify the approximate timing and direction of a recent episode of gene flow from the Sierra Nevada to the northern coastal ranges. The three subspecies comprise phylogenetic species, and the northern spotted owl (S. o. caurina) is sister to a clade of California (S. o. occidentalis) plus Mexican spotted owls (S o lucida); this represents a novel biogeographic pattern within birds. The California spotted owl had substantially lower nucleotide diversity than the other two subspecies; this result is inconsistent with present patterns of population density A causal explanation requires postulating a severe bottleneck or a selective sweep, either of which was confined to only one geographic region.     

Flyway structure in the circumpolar greater white‐fronted goose

Dispersal and migratory behavior are influential factors in determining how genetic diversity is distributed across the landscape. In migratory species, genetic structure can be promoted via several mechanisms including fidelity to distinct migratory routes. Particularly within North America, waterfowl management units have been delineated according to distinct longitudinal migratory flyways supported by banding data and other direct evidence. The greater white‐fronted goose (Anser albifrons) is a migratory waterfowl species with a largely circumpolar distribution consisting of up to six subspecies roughly corresponding to phenotypic variation. We examined the rangewide population genetic structure of greater white‐fronted geese using mtDNA control region sequence data and microsatellite loci from 23 locales across North America and Eurasia. We found significant differentiation in mtDNA between sampling locales with flyway delineation explaining a significant portion of the observed genetic variation (~12%). This is concordant with band recovery data which shows little interflyway or intercontinental movements. However, microsatellite loci revealed little genetic structure suggesting a panmictic population across most of the Arctic. As with many high‐latitude species, Beringia appears to have played a role in the diversification of this species. A common Beringian origin of North America and Asian populations and a recent divergence could at least partly explain the general lack of structure at nuclear markers. Further, our results do not provide strong support for the various taxonomic proposals for this species except for supporting the distinctness of two isolated breeding populations within Cook Inlet, Alaska (A. a. elgasi) and Greenland (A. a. flavirostris), consistent with their subspecies status.     

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America

Tree-feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3-4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.     

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.