The Balkans and the colonization of Europe: the post‐glacial range expansion of the wild boar,Sus scrofa

Aim We focus on the biogeographical role of the Balkan Peninsula as a glacial refugium and source of northward post‐glacial dispersal for many European taxa. Specifically, we analysed the genetic structure and variation of wild boar (Sus scrofa) samples primarily from Greece, a region that has repeatedly served as a glacial refugium within the Balkan Peninsula. Location Continental Greece, the Aegean island of Samos and Bulgaria. Methods We analysed wild boar samples from 18 localities. Samples from common domestic breeds were also examined to take into account interactions between wild and domesticated animals. Phylogenetic analyses were carried out on a 637‐bp fragment of the mitochondrial DNA control region in 200 wild boar and 27 domestic pigs. The sequences were also compared with 791 Eurasian wild boar and domestic pig D‐loop sequences obtained from GenBank. Results Ninety‐four haplotypes were identified in the European wild boar data set, of which 68 were found in the Balkan samples and assigned to two previously described clades: the E1 European and Near Eastern clades. All of the continental samples clustered in the E1 clade and the samples from Samos fell into the Near Eastern clade, consistent with the island’s proximity to Asia Minor. Intriguingly, 62 novel haplotypes were identified and are found exclusively in the Balkans. Only six haplotypes were shared between wild boar and domestic pigs. Main conclusions Our data reveal numerous novel and geographically restricted haplotypes in wild boar populations, suggesting the presence of separate refugia in the Balkans. Our analyses support the hypothesis of a post‐glacial wild boar expansion consistent with the leading edge model, north and west from modern day Greece, and suggest little maternal introgression of Near Eastern and domestic haplotypes into wild Balkan populations.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

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

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

Molecular markers provide evidence for a broad‐fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

• Calcareous grasslands belong to the most species‐rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post‐glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands.
• The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post‐glacial migration routes to Central Europe.
• Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post‐glacial recolonisation.
• Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post‐glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.

     

Mitochondrial population structure and post‐glacial dispersal of longnose sucker Catostomus catostomus in Labrador,Canada: evidence for multiple refugial origins and limited ongoing gene flow

Two hundred and eighty‐seven longnose sucker Catostomus catostomus were collected from 14 lakes in Labrador, 52 from three lakes in Ontario, 43 from two lakes in British Columbia and 32 from a lake in Yukon; a total of 414 in all. The resulting 34 haplotypes (20 in Labrador) contained moderate haplotypic diversity (h = 0·657) and relatively low nucleotide diversity (π = 3·730 × 10^?3. Mean ?_ST (0·453, P < 0·05) over all populations revealed distinct genetic structuring among C. catostomus populations across Canada, based on province, which was validated by the analysis and spatial analysis of molecular variance (c. 80% variation between provinces). These results probably reflect the historical imprint of recolonization from different refugia and possibly indicate limited ongoing gene flow within provinces. A haplotype network revealed one major and two minor clades within Labrador that were assigned to the Atlantic, Beringian and Mississippian refugia, respectively, with tests of neutrality and mismatch distribution indicative of a recent population expansion in Labrador, dated between c. 3500 and 8300 years ago.     

Ancient DNA identifies post‐glacial recolonisation,not recent bottlenecks,as the primary driver of contemporary mtDNA phylogeography and diversity in Scandinavian brown bears

Aim

Brown bear populations in Scandinavia show a strong mitochondrial DNA (mtDNA) phylogeographic structure and low diversity relative to other parts of Europe. Identifying the timing and origins of this mtDNA structure is important for conservation programs aimed at restoring populations to a natural state. Therefore, it is essential to identify whether contemporary genetic structure is linked to post‐glacial recolonisation from divergent source populations or an artefact of demographic impacts during recent population bottlenecks. We employed ancient DNA techniques to investigate the timing and potential causes of these patterns.

Location

Scandinavia and Europe.

Methods

Ancient mtDNA sequences from 20 post‐glacial Scandinavian bears were used to investigate phylogeographic structure and genetic diversity over the last 6000 years. MtDNA from 19 Holocene Norwegian bears was compared with 499 sequences from proximate extant populations in Sweden, Finland, Estonia and western Russia. A single mtDNA sequence from a Holocene Denmark sample was compared with 149 ancient and modern bears from Western Europe.

Results

All nineteen Holocene Norwegian samples are identical to or closely related to the most common mtDNA haplotype found in northern Europe today. MtDNA diversity was low and not significantly different from extant populations in northern Europe. In Denmark, we identified a single mtDNA haplotype that is previously unrecorded from Scandinavia.

Main conclusions

The current discrete phylogeographic structure and lack of mtDNA diversity in Scandinavia is attributed to serial founder effects during post‐glacial recolonisation from divergent source populations rather than an artefact of recent anthropogenic impacts. In contrast to previous interpretations, the recolonisation of southern Scandinavia may not have been limited to bears from a single glacial refugium. Results highlight the importance of conserving the long‐term evolutionary separation between northern and southern populations and identify southern Scandinavia as an important reservoir of mtDNA diversity that is under threat in other parts of Europe.

     

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow‐spotted mountain newt,Neurergus derjugini (Caudata: Salamandridae) in the mid‐Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.     

Palaeodistribution modelling and genetic evidence highlight differential post‐glacial range shifts of a rain forest conifer distributed across a latitudinal gradient

Aim We examine the range expansion/contraction dynamics during the last glacial cycle of the late‐successional tropical rain forest conifer Podocarpus elatus using a combination of modelling and molecular marker analyses. Specifically, we test whether distributional changes predicted by environmental niche modelling are in agreement with (1) the glacial maximum contractions inferred from the southern fossil record, and (2) population genetic‐based estimates of range disjunctions and demographic dynamics. In addition, we test whether northern and southern ranges are likely to have experienced similar expansion/contraction dynamics. Location Eastern Australian tropical and subtropical rain forests. Methods Environmental niche modelling was completed for three time periods during the last glacial cycle and was interpreted in light of the known palynology. We collected 109 samples from 32 populations across the entire range of P. elatus. Six microsatellite loci and Bayesian coalescence analysis were used to infer population expansion/contraction dynamics, and five sequenced loci (one plastid and four nuclear) were used to quantify genetic structure/diversity. Results Environmental niche modelling suggested that the northern and southern ranges of P. elatus experienced different expansion/contraction dynamics. In the northern range, the habitat suitable for P. elatus persisted in a small refugial area during the Last Glacial Maximum (LGM, 21 ka) and then expanded during the post‐glacial period. Conversely, in the south suitable habitat was widespread during the LGM but subsequently contracted. These differential dynamics were supported by Bayesian analyses of the population genetic data (northern dispersal) and are consistent with the greater genetic diversity in the south compared with the north. A contact zone between the two genetically divergent groups (corresponding to the Macleay Overlap Zone) was supported by environmental niche modelling and molecular analyses. Main conclusions The climatic fluctuations of the Quaternary have differentially impacted the northern and southern ranges of a broadly distributed rain forest tree in Australia. Recurrent contraction/expansion cycles contributed to the genetic distinction between northern and southern distributions of P. elatus. By combining molecular and environmental niche modelling evidence, this unique study undermines the general assumption that broadly distributed species respond in a uniform way to climate change.     

Genetic diversity and population demography of narrow‐ridged finless porpoises from South Korea on the basis of mitochondrial DNA variation: Implications for its conservation in East Asia

The narrow‐ridged finless porpoise, Neophocaena asiaeorientalis, is a small cetacean inhabiting the coastal waters of Taiwan, East China, Korea, and Japan, and is threatened by anthropogenic pressures. We analyzed the mitochondrial control region of N. asiaeorientalis from Korean waters to better understand the genetic diversity, population structure, and demographic history through comparison with previously reported sequences from Japan and China. The network and phylogenetic analyses illustrated that few common haplotypes were centrally located, and that most of the unique haplotypes were on the edges of the branches. Moreover, the estimates of fixation indices and gene flow revealed significant genetic differentiation among the regions in East Asia. The mismatch distribution suggested a population expansion in the late Pleistocene. We propose that the current genetic variation in N. asiaeorientalis in East Asia was shaped by historical demographic expansion during the late Pleistocene. The N. asiaeorientalis of Korean waters has the highest diversity observed with the sampling available at this time. However, the abundance of N. asiaeorientalis in Korean waters has rapidly declined and the species needs protection from bycatch.