The influence of Pleistocene glacial refugia on tawny owl genetic diversity and phylogeography in western Europe

The glacial refugia hypothesis indicates that during the height of the Pleistocene glaciations the temperate species that are today widespread in western Europe must have survived in small and climatically favourable areas located in the southern peninsulas of Iberia, Italy and Balkans. One such species is the tawny owl, a relatively sedentary, nonmigratory bird presently distributed throughout Europe. It is a tree-nesting species closely associated with deciduous and mixed coniferous woodlands. In this study I used control region mtDNA sequences from 187 individuals distributed among 14 populations to determine whether current genetic patterns in tawny owl populations were consistent with postglacial expansion from peninsular refugia. European, North African and Asian tawny owls were found to represent three distinct lineages, where North Africa is the sister clade to all European owls. Within Europe, I found three well-supported clades that correspond to each of the three allopatric refugia. Expansion patterns indicate that owls from the Balkan refugium repopulated most of northern Europe, while expansion out of Iberia and Italy had only regional effects leading to admixture in France. Estimates of population divergence times between refugia populations are roughly similar, but one order of magnitude smaller between Greece and northern Europe. Based on a wide range of mutation rates and generation times, divergence between refugia appears to date to the Pleistocene.     

A long‐standing Pleistocene refugium in southern Africa and a mosaic of refugia in East Africa: insights from mtDNA and the common eland antelope

Aim Previous genetic studies of African savanna ungulates have indicated Pleistocene refugial areas in East and southern Africa, and recent palynological, palaeovegetation and fossil studies have suggested the presence of a long‐standing refugium in the south and a mosaic of refugia in the east. Phylogeographic analysis of the common eland antelope, Taurotragus oryx (Bovidae), was used to assess these hypotheses and the existence of genetic signatures of Pleistocene climate change. Location The sub‐Saharan savanna biome of East and southern Africa. Methods Mitochondrial DNA control‐region fragments (414 bp) from 122 individuals of common eland were analysed to elucidate the phylogeography, genetic diversity, spatial population structuring, historical migration and demographic history of the species. The phylogeographic split among major genetic lineages was dated using Bayesian coalescent‐based methods and a calibrated fossil root of 1.6 Ma for the split between the common eland and the giant eland, Taurotragus derbianus. Results Two major phylogeographic lineages comprising East and southern African localities, respectively, were separated by a net nucleotide distance of 4.7%. A third intermediate lineage comprised only three haplotypes, from Zimbabwe in southern Africa. The estimated mutation rate of 0.097 Myr^?1 revealed a more recent common ancestor for the eastern lineage (0.21 Ma; 0.07–0.37) than for the southern lineage (0.35 Ma; 0.10–0.62). Compared with the latter, the eastern lineage showed pronounced geographic structuring, lower overall nucleotide diversity, higher population differentiation, and isolation‐by‐distance among populations. Main conclusions The data support the hypothesis of Pleistocene refugia occurring in East and southern Africa. In agreement with palynological, palaeovegetation and fossil studies, our data strongly support the presence of a longer‐standing population in the south and a mosaic of Pleistocene refugia in the east, verifying the efficacy of genetic tools in addressing such questions. The more recent origin of the common eland inhabiting East Africa could result from colonization following extinction from the region. Only two other dated African ungulate phylogenies have been published, applying different methods, and the similarity of dates obtained from the three distinct approaches indicates a significant event c. 200 ka, which left a strong genetic signature across a range of ungulate taxa.     

Putative glacial refugia of Cedrus atlantica deduced from Quaternary pollen records and modern genetic diversity

Aim To investigate the impact of past environmental changes on Cedrus atlantica and its current genetic diversity, and to predict its future distribution. Location Morocco, Algeria and Tunisia. Methods Eleven fossil pollen records from these three countries were used to locate putative glacial refugia and to reconstruct past climate changes. A mechanistic vegetation distribution model was used to simulate the distribution of C. atlantica in the year 2100. In addition, a genetic survey was carried out on modern Moroccan C. atlantica. Results Pollen records indicate that Cedrus was present during the last glacial period, probably in scattered refugia, in Tunisia, Algeria and Morocco. In the Tunisian and Algerian sites, cedar expanded during the late glacial and the early Holocene, then disappeared after c. 8000 yr bp . Reconstructed mean annual precipitation and January temperature show that the last glacial period in Morocco was cooler by 10–15°C and drier by c. 300–400 mm year^?1 than the climate today. Modern chloroplast microsatellites of 15 C. atlantica populations in Morocco confirm the presence of multiple refugia and indicate that cedar recolonized the Moroccan mountains fairly recently. Model simulation indicates that by the year 2100 the potential distribution of C. atlantica will be much restricted with a potential survival area located in the High Atlas. Main conclusions Environmental changes in northern Africa since the last glacial period have had an impact on the geographical distribution of C. atlantica and on its modern genetic diversity. It is possible that by the end of this century C. atlantica may be unable to survive in its present‐day locations. To preserve the species, we suggest that seedlings from modern C. atlantica populations located in the High Atlas mountains, where a high genetic diversity is found, be transplanted into the western High Atlas.     

Pliocene and Pleistocene events shaping the genetic diversity within the central corridor of the Brazilian Atlantic Forest

Dinoponera lucida (Formicidae; Ponerinae) is an extinction‐threatened species of ant which is endemic in the central corridor of the Atlantic Forest. We used mitochondrial sequences of the Cox1, Cox2 and Cytb genes in order to infer some aspects of the evolutionary history and phylogeography of this ant. High genetic divergence and population structure were observed for the whole species. The current pattern of D. lucida diversity seems to be shaped during different geological times: middle Pliocene, early Pleistocene and mainly late Pleistocene, when the reduction of populations generated a structure pattern of the genetic variation of this species. Our data show that this structure results from the maintenance of populations of D. lucida within very small putative refuges to the south of the central Bahia refugium. We thus argue that, for some Atlantic forest endemic species, especially those resistant to very small fragments of forest, such as D. lucida, the small putative refuges were as important as, or even more important than, larger and stable refuges for the creation and maintenance of diversity, adding another piece to the puzzle of the mechanisms underlying local endemism. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 949–960.     

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa

Evidence from numerous Pan‐African savannah mammals indicates that open‐habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open‐habitat formations throughout sub‐Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.     

Evaluating mechanisms of diversification in a Guineo‐Congolian tropical forest frog using demographic model selection

The accumulation of biodiversity in tropical forests can occur through multiple allopatric and parapatric models of diversification, including forest refugia, riverine barriers and ecological gradients. Considerable debate surrounds the major diversification process, particularly in the West African Lower Guinea forests, which contain a complex geographic arrangement of topographic features and historical refugia. We used genomic data to investigate alternative mechanisms of diversification in the Gaboon forest frog, Scotobleps gabonicus, by first identifying population structure and then performing demographic model selection and spatially explicit analyses. We found that a majority of population divergences are best explained by allopatric models consistent with the forest refugia hypothesis and involve divergence in isolation with subsequent expansion and gene flow. These population divergences occurred simultaneously and conform to predictions based on climatically stable regions inferred through ecological niche modelling. Although forest refugia played a prominent role in the intraspecific diversification of S. gabonicus, we also find evidence for potential interactions between landscape features and historical refugia, including major rivers and elevational barriers such as the Cameroonian Volcanic Line. We outline the advantages of using genomewide variation in a model‐testing framework to distinguish between alternative allopatric hypotheses, and the pitfalls of limited geographic and molecular sampling. Although phylogeographic patterns are often species‐specific and related to life‐history traits, additional comparative studies incorporating genomic data are necessary for separating shared historical processes from idiosyncratic responses to environmental, climatic and geological influences on diversification.     

Mitochondrial differentiation in a polymorphic land snail: evidence for Pleistocene survival within the boundaries of permafrost

The genetic differentiation of populations having colonized formerly unsuitable habitats after the Pleistocene glaciations depends to a great extent on the speed of expansion. Slow dispersers maintain their refugial diversity whereas fast dispersal leads to a reduction of diversity in the newly colonized areas. During the Pleistocene, almost the entire current range of the land snail Arianta arbustorum has repeatedly been covered with ice or been subjected to permafrost. Owing to the low potential for dispersal of land snails, slow (re)colonization of the wide range from southern refugia can be excluded. Alternatively, fast, passive dispersal from southern refugia or survival in and expansion from multiple refugia within the area subjected to permafrost may account for the current distribution. To distinguish between these scenarios we reconstructed a phylogeography based on the sequences of a fragment of the cytochrome oxidase I from 133 individuals collected at 45 localities and analysed the molecular variance. Seventy-five haplotypes were found that diverged on average at 7.52% of positions. This high degree of diversity suggests that A. arbustorum is an old species in which the population structure, isolation and the hermaphroditic nature have reduced the probability of lineage extinction. The genetic structure was highly significant with the highest variance partition found among regions. Geographic distance and mitochondrial differentiation were not congruent. Lineages had overlapping ranges. The clear genetic differentiation and the patchy pattern of haplotype distribution suggest that colonization of formerly unsuitable habitats was mainly achieved from multiple populations from within the permafrost area.     

Mitochondrial gene diversity in the common vole Microtus arvalis shaped by historical divergence and local adaptations

The phylogeography of the common vole (Microtus arvalis) was examined by analysing mitochondrial DNA (mtDNA) sequence variation in 1044 base pairs (bp) of the cytochrome b (cytb) gene and in 322 bp of the control region (ctr) among 106 individuals from 58 locations. The geographical distribution of four previously recognized cytb evolutionary lineages in Europe was refined and a new lineage was found in southern Germany. All lineages were distributed allopatrically, except in one sample that was probably located in a contact zone. The occurrence of several lineages in the Alps is in keeping with their recent recolonization from distinct sources. The translation of 84 cytb DNA sequences produced 33 distinct proteins with relationships that differed from those of the DNA haplotypes, suggesting that the mtDNA lineages did not diverge in response to selection. In comparison with M. agrestis, a neutrality test detected no overall evidence for selection in the cytb gene, but a closer examination of a structural model showed that evolutionarily conserved and functionally important positions were often affected. A new phylogeographical test of random accumulation of nonsynonymous mutations generated significant results in three lineages. We therefore conclude that the molecular diversity of cytb in M. arvalis is overall the result of the demographic history of the populations, but that there have been several episodes of local adaptation to peculiar environments.     

Multiple genetic divergences and population expansions of a Mediterranean sandfly,Phlebotomus ariasi,in Europe during the Pleistocene glacial cycles

Phlebotomus ariasi is one of the two sandflies transmitting the causative agent of zoonotic leishmaniasis, Leishmania infantum, in France and Iberia, and provides a rare case study of the postglacial re-colonization of France by a Mediterranean species. Four DNA sequences were analysed—mitochondrial cytochrome b (cyt b), nuclear elongation factor-1α (EF-1α) and two anonymous nuclear loci—for 14–15 French populations and single populations from northeast Spain, northwest Spain, Portugal and Morocco. The presence of cryptic sibling species was not revealed by phylogenetic analyses and testing for reproductive isolation between sympatric populations defined by the two most divergent cyt b haplogroups. No locus was shown to be under positive directional or balancing selection and, therefore, molecular variation was explained demographically. Each nuclear locus showed shallow isolation by distance from Portugal to the French Pyrenees, but for both cyt b and EF-1α there was then a step change to the upland Massif Central, where leading-edge populations showed low diversity at all loci. Multiple genetic divergences and population expansions were detected by analyses of cyt b and dated to the Pleistocene. Endemicity of one cyt b sub-lineage suggested the presence of a refuge north of the Pyrenees during the last glacial period. Monopolization of the Massif Central by genetically differentiated populations of P. ariasi might possibly hinder the northwards spread of leishmaniasis.     

Phylogeographic structure and regional history of Lemniscomys striatus (Rodentia: Muridae) in tropical Africa

Aim This study aims to elucidate the phylogeography of the murid rodent Lemniscomys striatus and to evaluate the relative roles of ecological change, habitat patchiness, rivers and geological barriers in structuring patterns of diversity. Location Sub‐Saharan Africa. Methods The extent of phylogeographic patterns and molecular genetic diversity (cytochrome b gene) were addressed in a survey of 128 individuals of L. striatus from 42 localities. Using maximum parsimony, maximum likelihood, Bayesian, network and genetic structure analyses, we inferred intraspecific relationships and tested hypotheses for historical patterns of gene flow within L. striatus. Results Our results identified four major geographical clades within L. striatus: a West African clade, a Benin‐Nigeria clade, a Central African clade, and an East African clade. Several subclades were identified within these four major clades. Restricted gene flow with isolation by distance was recorded, which is congruent with the low dispersal ability of such a small murid rodent. No clear signal of population expansion was detected within clades or subclades. Main conclusions The western rift system and the Volta and Niger rivers may have acted as long‐term extrinsic barriers to gene flow, resulting in the emergence of the four main clades of L. striatus. The observed pattern of mitochondrial variation observed within each clade probably results from late Pleistocene climatic and vegetation changes: during adverse conditions (forest expansion), L. striatus may have survived only in refugia, and then experienced range expansion under favourable conditions (savanna expansion).     

The ancient tropical rainforest tree Symphonia globulifera L. f. (Clusiaceae) was not restricted to postulated Pleistocene refugia in Atlantic Equatorial Africa

Understanding the history of forests and their species'' demographic responses to past disturbances is important for predicting impacts of future environmental changes. Tropical rainforests of the Guineo-Congolian region in Central Africa are believed to have survived the Pleistocene glacial periods in a few major refugia, essentially centred on mountainous regions close to the Atlantic Ocean. We tested this hypothesis by investigating the phylogeographic structure of a widespread, ancient rainforest tree species, Symphonia globulifera L. f. (Clusiaceae), using plastid DNA sequences (chloroplast DNA [cpDNA], psbA-trnH intergenic spacer) and nuclear microsatellites (simple sequence repeats, SSRs). SSRs identified four gene pools located in Benin, West Cameroon, South Cameroon and Gabon, and São Tomé. This structure was also apparent at cpDNA. Approximate Bayesian Computation detected recent bottlenecks approximately dated to the last glacial maximum in Benin, West Cameroon and São Tomé, and an older bottleneck in South Cameroon and Gabon, suggesting a genetic effect of Pleistocene cycles of forest contraction. CpDNA haplotype distribution indicated wide-ranging long-term persistence of S. globulifera both inside and outside of postulated forest refugia. Pollen flow was four times greater than that of seed in South Cameroon and Gabon, which probably enabled rapid population recovery after bottlenecks. Furthermore, our study suggested ecotypic differentiation—coastal or swamp vs terra firme—in S. globulifera. Comparison with other tree phylogeographic studies in Central Africa highlighted the relevance of species-specific responses to environmental change in forest trees.     

Population genetic structure of Grauer's Swamp Warbler Bradypterus graueri,an Albertine Rift endemic

The endangered warbler Bradypterus graueri is endemic to the Albertine Rift, where it is restricted to montane swamps above 1900 m across the region. We studied genetic structure among six populations sampled across the species' distribution in northern Burundi, Rwanda, Uganda and the eastern Democratic Republic of Congo. A total of 2117 base pairs of mitochondrial data were sequenced. Phylogenetic analyses and network reconstruction of B. graueri haplotypes recovered three clades with a defined geographical pattern: clade 1, Virunga Volcanoes and Kigezi Highlands; clade 2, Rugege Highlands; and clade 3, Kahuzi‐Biega Highlands; clades 2 and 3 are sisters to each other. Both landscape dynamics and historical climate are likely to have played a role in the diversification of this species. The divergence between clade 1 and clades 2 and 3 (168.5 ka, 95% HPD 108.5, 244.4) coincides with a prolonged period of aridity in tropical Africa between 130 and 270 ka. Similarly, the divergence between clades 2 and 3 (99.4 ka, 95% HPD 55.4, 153.8) corresponds with a period of aridity just prior to 94 ka. Populations sampled from the eastern arm of the central Albertine Rift (Kigezi and Rugege Highlands) show a coincident increase in effective population size after the Last Glacial Maximum at c. 15 ka, whereas those sampled from Kahuzi‐Biega on the western arm of the rift do not. Despite the perceived higher vagility of bird species relative to other vertebrates, the degree of phylogeographical structure among populations of B. graueri is similar to that reported for small mammals (Hylomyscus vulcanorum, Lophuromys woosnami, Sylvisorex vulcanorum) and a frog Hyperolius castaneus sampled across the central Albertine Rift. Collectively our results suggest that climate dynamics associated with late Pleistocene cycles had a significant influence on driving the population genetic structure and associated levels of genetic diversity in B. graueri and other small terrestrial vertebrates. Our results have implications for the conservation of B. graueri and other endemics to the Albertine Rift, particularly in the context of other phylogegeographical studies centred on this biodiversity hotspot.     

Phylogeography of sugar kelp: Northern ice‐age refugia in the Gulf of Alaska

Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (Φ_ST = 0.834, p < .0001). Microsatellites showed moderate levels of heterozygosity within populations (mean H _E = 0.391). Kelps in the same organellar lineage tended to cluster together, regardless of geographic origins, as indicated in a principal coordinate analysis (PCoA) of microsatellite genotypes. The PCoA also showed evidence of nuclear hybridizations between co‐occurring organellar lineages. Individual admixture plots with population clusters of K = 2, 6, and 9 showed increasing complexity with considerable historical admixture between some clusters. A time‐calibrated phylogeny placed divergences between rbcL‐COI lineages at 1.4 million years at most. The time frames of mutation in the rbcL‐COI lineages and microsatellite population clusters differed among locations. The existence of ancient lineages in the Gulf of Alaska, moderate levels of genetic diversity, and the absence of departures from neutrality are consistent with northern refugia during multiple Croll‐Milankovitch climate cycles in the Pleistocene Epoch.