The Multiple Origins of Roe Deer Populations in Western Iberia and Their Relevance for Conservation

Simple Summary

The Iberian Peninsula is considered a reservoir of genetic diversity and the source for the recolonization of continental Europe by several mammal species, after the last glaciation period. Here, we intended to characterize the genetic patterns and origins of one of those species, the roe deer, through the analysis of different genetic markers, as there is a major knowledge gap about the species’ populations at the western edge of roe deer distribution in the Iberian Peninsula. We found that western Iberia is a diversity hotspot for roe deer, with shared gene pools with other European and Iberian regions, but also with unique genetic elements, particularly the case of the relict population of Peneda Gerês National Park. Due to the distinct genetic diversity that was observed in western Iberian populations, we highlight the importance of these populations as sources of resilience against global changes. Our results provide useful information for the management and conservation of this species in the Iberian Peninsula. We advise transboundary management between Portugal and Spain as a rule, as well as careful evaluation of reintroduction actions, that should take in account the genetic data, in order to maintain the genetic heritage of roe deer in Europe.

Abstract

The roe deer (Capreolus capreolus) is native and widespread in Europe and its phylogeography has been clarified in the last decades. Southern peninsulas are considered as reservoirs of genetic diversity and the source for the recolonization of Europe after the last glacial maximum. Even though roe deer populations have been genetically characterized, there is a major knowledge gap about the populations at the western edge of its distribution. To fill this caveat, and based on mitochondrial and nuclear DNA data, we aim to: (i) characterize the genetic diversity and structure of roe deer in western Iberia; (ii) clarify the origins and phylogeographical affinities of these populations, namely the relict population from Peneda Gerês National Park (PNPG, Portugal) and the likely allochthonous populations from central and south (CS) Portugal; (iii) discuss the implications of our findings for the management and conservation of the roe deer. Three major genetic clusters were inferred based on nuclear genotypes and were structured in a similar way as the three major mtDNA clades present in Iberia. Patterns inferred with nuclear markers confirmed PNPG as a relict population. Roe deer from CS Portugal share haplotypes with Central Europe rather than with other western Iberian populations, confirming its mainly allochthonous origin. Our results highlight western Iberia as a diversity hotspot for roe deer. We highlight the role of intraspecific genetic diversity as a source of resilience against ongoing global changes; the need for transboundary management and the importance of genetic data to inform management and conservation. When considered, repopulation or translocation measures should follow the IUCN Law of Reintroductions and meticulously conducted in order to preserve the genetic heritage of the species.

Evolutionary history and phylogeographic relationships of shrews from Sorex araneus group

Shrews of the Sorex genus are an evolutionarily successful group that includes more than 77 species widely distributed in Eurasia and North America. The genus is one of the rare cases where karyotypic changes reflect well the evolutionary relationships among its species. The taxa showing the greatest variation in karyotype are usually classified into the Sorex araneus group. Its evolution was associated with chromosomal rearrangements, which could have promoted fast diversification of this group into many chromosomal races and species. These processes were additionally complicated by introgressions of mitochondrial DNA, which made the evolutionary history of this group quite complex and difficult to infer. To tackle the problem, we performed multi-method phylogenetic analyses based on mitochondrial cytochrome b that is considered a good molecular marker available for many representatives of Sorex. The results were compared with phylogenies based on chromosomal rearrangement data and put into temporal and spatial context using molecular dating and historical biogeography methods. We complemented the study with the estimation of diversification rates within the S. araneus group as well as comparing the results with paleontological records and climatic oscillations within the last 4 million years. Based on the gathered data, we proposed a hypothetical scenario for the evolution and geographic dispersion of species belonging to the S. araneus group. The shrews began to diversify about 2.7 million years ago in Eurasia and then migrated at least twice to North America. The evolution of shrews was driven by Pleistocene glacial and interglacial cycles, which increased their speciation rate and the emergence of new lineages. The migrations of populations were accompanied by introgressions of mitochondrial DNA into native shrews and occurred at least twice.

Altitudinal Barrier to the Spread of an Invasive Species: Could the Pyrenean Chain Slow the Natural Spread of the Pinewood Nematode?

Mountain ranges may delimit the distribution of native species as well as constitute potential barriers to the spread of invasive species. The invasive pinewood nematode, Bursaphelenchus xylophilus, is a severe forest pest inducing pine wilt disease. It is vectored in Europe by a native long-horned beetle, Monochamus galloprovincialis. This study explored the potential of the Pyrenean chain to slow or prevent the natural spread of nematode-infested beetles from the Iberian Peninsula, where the nematode is established and is expanding its range, towards France and the rest of Europe. An analysis of the genetic structure and migration patterns of the beetle populations throughout the Pyrenean mountain range was combined with a spread model simulating the potential movements of nematode-infested beetles across it. The central part of the Pyrenees, which corresponds to the highest elevation zone, was shown to prevent gene flow between the French and Spanish populations of M. galloprovincialis on each side of the mountains. Conversely, strong admixture was detected between populations located on both sides of low elevation hills, and especially at the east and west extremities of the mountain range. Simulations of the spread of nematode-infested beetles under various thresholds of beetle survival and pine wilt disease expression gave results consistent with the variation in genetic make-up, suggesting that western and eastern hillsides may represent corridors favoring natural spread of the nematode from the Iberian Peninsula to France. Simulations also showed that temperature rise due to climate change may significantly reduce the extent of the barrier formed by highest elevations. Our results support the hypothesis that the Pyrenean chain represents a partial barrier to the natural spread of nematode-infested beetles. These results, which have to be considered together with potential human-assisted long-distance spread of the nematode, highlight priority zones for future pest monitoring and management programs. More generally, such an integrated approach could be used to assess the role of mountain chains in the potential spread of other invasive pests.

The telomeric sync model of speciation: species-wide telomere erosion triggers cycles of transposon-mediated genomic rearrangements, which underlie the saltatory appearance of nonadaptive characters

Charles Darwin knew that the fossil record is not overwhelmingly supportive of genetic and phenotypic gradualism; therefore, he developed the core of his theory on the basis of breeding experiments. Here, I present evidence for the existence of a cell biological mechanism that strongly points to the almost forgotten European concept of saltatory evolution of nonadaptive characters, which is in perfect agreement with the gaps in the fossil record. The standard model of chromosomal evolution has always been handicapped by a paradox, namely, how speciation can occur by spontaneous chromosomal rearrangements that are known to decrease the fertility of heterozygotes in a population. However, the hallmark of almost all closely related species is a differing chromosome complement and therefore chromosomal rearrangements seem to be crucial for speciation. Telomeres, the caps of eukaryotic chromosomes, erode in somatic tissues during life, but have been thought to remain stable in the germline of a species. Recently, a large human study spanning three healthy generations clearly found a cumulative telomere effect, which is indicative of transgenerational telomere erosion in the human species. The telomeric sync model of speciation presented here is based on telomere erosion between generations, which leads to identical fusions of chromosomes and triggers a transposon-mediated genomic repatterning in the germline of many individuals of a species. The phenotypic outcome of the telomere-triggered transposon activity is the saltatory appearance of nonadaptive characters simultaneously in many individuals. Transgenerational telomere erosion is therefore the material basis of aging at the species level.

Chromosomal rearrangements do not seem to affect the gene flow in hybrid zones between karyotypic races of the common shrew (Sorex araneus)

Chromosomal rearrangements are proposed to promote genetic differentiation between chromosomally differentiated taxa and therefore promote speciation. Due to their remarkable karyotypic polymorphism, the shrews of the Sorex araneus group were used to investigate the impact of chromosomal rearrangements on gene flow. Five intraspecific chromosomal hybrid zones characterized by different levels of karyotypic complexity were studied using 16 microsatellites markers. We observed low levels of genetic differentiation even in the hybrid zones with the highest karyotypic complexity. No evidence of restricted gene flow between differently rearranged chromosomes was observed. Contrary to what was observed at the interspecific level, the effect of chromosomal rearrangements on gene flow was undetectable within the S. araneus species.

Mitochondrial and chromosomal insights into karyotypic evolution of the pygmy mouse, Mus minutoides, in South Africa

The African pygmy mouse, Mus minutoides, displays extensive Robertsonian (Rb) diversity. The two extremes of the karyotypic range are found in South Africa, with populations carrying 2n = 34 and 2n = 18. In order to reconstruct the scenario of chromosomal evolution of M. minutoides and test the performance of Rb fusions in resolving fine-scale phylogenetic relationships, we first describe new karyotypes, and then perform phylogenetic analyses by two independent methods, using respectively mitochondrial cytochrome b sequences and chromosomal rearrangements as markers. The molecular and chromosomal phylogenies were in perfect congruence, providing strong confidence both for the tree topology and the chronology of chromosomal rearrangements. The analysis supports a division of South African specimens into two clades showing opposite trends of chromosomal evolution, one containing all specimens with 34 chromosomes (karyotypic stasis) and the other grouping all mice with 18 chromosomes that have further diversified by the fixation of different Rb fusions (extensive karyotypic reshuffling). The results confirm that Rb fusions are by far the predominant rearrangement in M. minutoides but strongly suggest that recurrent whole-arm reciprocal translocations have also shaped this genome.

Chromosomal rearrangements and gene flow over time in an inter-specific hybrid zone of the Sorex araneus group

Most hybrid zones have existed for hundreds or thousands of years but have generally been observed for only a short time period. Studies extending over periods long enough to track evolutionary changes in the zones or assess the ultimate outcome of hybridization are scarce. Here, we describe the evolution over time of the level of genetic isolation between two karyotypically different species of shrews (Sorex araneus and Sorex antinorii) at a hybrid zone located in the Swiss Alps. We first evaluated hybrid zone movement by contrasting patterns of gene flow and changes in cline parameters (centre and width) using 24 microsatellite loci, between two periods separated by 10 years apart. Additionally, we tested the role of chromosomal rearrangements on gene flow by analysing microsatellite loci located on both rearranged and common chromosomes to both species. We did not detect any movement of the hybrid zone during the period analysed, suggesting that the zone is a typical tension zone. However, the gene flow was significantly lower among the rearranged than the common chromosomes for the second period, whereas the difference was only marginally significant for the first period. This further supports the role of chromosomal rearrangements on gene flow between these taxa.

Using a Bayesian method to assign individuals to karyotypic taxa in shrew hybrid zones

Individuals sampled in hybrid zones are usually analysed according to their sampling locality, morphology, behaviour or karyotype. But the increasing availability of genetic information more and more favours its use for individual sorting purposes and numerous assignment methods based on the genetic composition of individuals have been developed. The shrews of the Sorex araneus group offer good opportunities to test the genetic assignment on individuals identified by their karyotype. Here we explored the potential and efficiency of a Bayesian assignment method combined or not with a reference dataset to study admixture and individual assignment in the difficult context of two hybrid zones between karyotypic species of the Sorex araneus group. As a whole, we assigned more than 80% of the individuals to their respective karyotypic categories (i.e. 'pure' species or hybrids). This assignment level is comparable to what was obtained for the same species away from hybrid zones. Additionally, we showed that the assignment result for several individuals was strongly affected by the inclusion or not of a reference dataset. This highlights the importance of such comparisons when analysing hybrid zones. Finally, differences between the admixture levels detected in both hybrid zones support the hypothesis of an impact of chromosomal rearrangements on gene flow.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples. Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene. For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula. In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period. The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples.

Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene.

For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula.

In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period.

The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.

Lack of mitochondrial DNA divergence between chromosome races of the common shrew, Sorex araneus, in Sweden. Implications for interpreting chromosomal evolution and colonization history

The common shrew, Sorex araneus, has one of the most variable karyotypes among mammals, displaying numerous chromosome races throughout its distribution. The six chromosome races present in Sweden can be categorized in two different karyotypic groups, the west and north European karyotypic groups (western and northern). Three races belonging to the western group are considered to have arisen through whole arm reciprocal translocations (WARTs). Race formation through this process requires a bottleneck event. In the present study we sequenced a part of the mitochondrial DNA (mtDNA) genome to investigate molecular differences between the chromosome races in Sweden. We found no mtDNA differentiation between the mainland chromosome races or the karyotypic groups. Genetic variation is as large between populations within a race as between populations among the races or karyotypic groups, suggesting that the karyotypic groups might have originated in a common glacial refugium. The noticeable exception is the Oland race, which shows higher mtDNA diversity compared to the other Swedish races, indicating a divergent origin difficult to explain. Mitochondrial DNA variation in Sweden suggests that most haplotypes arose in situ and that the populations has undergone a rapid size expansion. Altogether, the mtDNA data are in agreement with the WART hypothesis, which still holds as the most plausible variant of karyotype evolution for three of the chromosome races of the common shrew in Sweden.

Phylogeographic structure and postglacial evolutionary history of water voles (Arvicola terrestris) in the United Kingdom

The phylogeographic pattern of mitochondrial DNA variation in water voles (Arvicola terrestris) from 57 localities across the United Kingdom and representative samples from Spain, France, Switzerland and Finland was determined from sequence variation in the central portion of the control region. Twenty-seven different haplotypes were resolved which formed two distinct phylogenetic clades. This major division separated haplotypes found in Scotland from those found in England and Wales. Nested clade analysis of haplotypes indicated that such a division was a consequence of allopatric fragmentation. The haplotypes found in Switzerland, France and Spain clustered with Scottish haplotypes, whereas the haplotype from Finland clustered with the English/Welsh haplotypes. These patterns indicate that contemporary Scottish populations are derived from an Iberian glacial refugium, whereas English and Welsh populations are derived from an eastern European refugium. As such, the postglacial recolonization of the United Kingdom must have involved two colonization events, either in different localities with no subsequent contact, or as two waves separated over time, with the second wave of colonizers displacing the first. An analysis of molecular variance (AMOVA) identified significant population genetic divergence within both the major clades, indicative of restricted gene flow and regional population isolation. The implications of both phylogeographical and population genetic structure are discussed in context with the conservation of water voles in Britain.

Mitochondrial phylogeography of the Woodmouse (Apodemus sylvaticus) in the Western Palearctic region

We sequenced 965 base pairs of the mitochondrial DNA cytochrome b from 102 woodmice (Apodemus sylvaticus) collected from 40 European localities. The aims of the study were to answer the following questions. (i) Did the Mediterranean peninsulas play a role as refuge for woodmice? (ii) Is genetic variability of A. sylvaticus higher in the Mediterranean region compared with northern Europe? (iii) Are the patterns of the postglacial colonization of Europe by woodmice similar to those presently recognized for other European species? The results provide a clear picture of the impact of the Quaternary glaciations on the genetic and geographical structure of the woodmouse. Our analyses indicate a higher genetic variability of woodmice in the Mediterranean peninsulas compared to northern Europe, suggesting a role of the former as refuge regions for this small mammal. An original pattern of postglacial colonization is proposed where the Iberian and southern France refuge populations colonized almost all European regions. The Sicilian population appears to be very differentiated and highly variable. This emphasizes the importance of this island as a 'hot spot' for the intraspecific genetic diversity of the woodmouse. Finally, woodmice in North Africa originated from southwestern Europe, most probably as a result of a recent anthropogenic introduction.

Wind-dispersed pollen mediates postglacial gene flow among refugia

A long-term genetic legacy of refugial isolation has been postulated and was demonstrated for maternal refugial lineages for numerous plant and animal species. The lineages were assumed to have remained separated from each other for several glacial periods. The conifer Abies alba Miller, silver fir, is an excellent model to test whether pollen-mediated gene flow may eliminate the genetic imprints of Pleistocene refugial isolation. Two DNA markers with contrasting modes of inheritance were applied to 100 populations covering the entire range of silver fir in Europe. The markers exhibited each two highly conserved alleles based on an insertion/deletion of 80 bp in the fourth intron of the mitochondrial nad5 gene and on a synonymous substitution in the chloroplast psbC gene. The geographical distribution of the maternally inherited mitochondrial variation supported the existence of at least two refugia with two recolonizing maternal lineages remaining largely separated throughout the range. The cline of the nad5 allele frequencies was much steeper than the one of the two psbC alleles. The psbC cline was as wide as the whole range of the species. Our results provide striking evidence that even a species with very long generation times and heavy pollen grains was able to establish a highly efficient pollen-mediated gene flow between refugia. Therefore we postulate that an exchange of genetic information between refugia by range-wide paternal introgression is possible in wind-pollinated plant species.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus)

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

The evolutionary history of the two karyotypic groups of the common shrew, Sorex araneus, in Poland

Genetic variability within and among two karyotypic groups and five chromosome races of the common shrew (Sorex araneus) in Poland was assayed by sequencing a 1023 bp part of the cytochrome b gene (mtDNA) from 28 individuals. Thirty-four variable positions defined 21 distinct haplotypes with a maximum sequence divergence of 0.88%. No significant differentiation in the cytochrome b gene between Western and Eastern Karyotypic groups was found. Haplotype diversity estimates within the races and groups sampled were high (h = 0.800-0.928), while nucleotide diversity estimates were low (pi = 0.0034-0.0053). The distribution of pairwise nucleotide differences fits well with expectations of a "sudden expansion" model. High haplotype diversity was accompanied by relatively high expected heterozygosity (H(E)) values in nuclear genes (calculated over 47 enzyme loci: H(E) = 0.031 - 0.049), giving no evidence for a recent bottleneck after the process of post-Pleistocene recolonization of Poland by the shrews. Thus, for S. araneus chromosome races in Poland, the data on the cytochrome b gene variability support the hypothesis assuming the Robertsonian fusions having spread into an ancestral acrocentric distribution.

Chromosome races with Pliocene origins: evidence from mtDNA

There are eight distinct chromosomal races of the New Zealand weta Hemideina thoracica. We used mtDNA sequence data to test the hypothesis that these races originated on islands during the early Pliocene (7--4 million years ago). Nine major mitochondrial lineages were identified from 65 cytochrome oxidase I sequences. Phylogenetic analysis of these lineages suggests that they arose at approximately the same time. The geographical distribution of some lineages coincides with areas that were islands during the Pliocene. Overall, hierarchical AMOVA analysis shows that chromosomal races and Pliocene islands describe only 28% and 24%, respectively, of the total current mtDNA variation. However, removing one widespread (A) and one putatively introgressed (F) lineage increases these estimates to 65% and 80%, respectively. Intraspecific sequence divergence was very high, reaching a maximum of 9.5% (uncorrected distance) and GC content was high compared to other insect mtDNA sequences. Average corrected distance among mtDNA lineages supports the Pliocene origins of this level of genetic diversity. In the southern part of the species range there is reduced mtDNA variation, probably related to local extinction of H. thoracica populations from recent volcanic activity and subsequent re-colonization from a leading edge. In contrast, in this southern part there are five chromosome races, suggesting that chromosome races here may be younger than those in the north.

On the mtDNA restriction patterns variation of the Iberian wood mouse (Apodemus sylvaticus). Comparison with other west Mediterranean populations

In previous studies, the presence of three main lineages of wood mice mtDNA was described in western Europe: a first one distributed from the Pyrenees to Scandinavia. a Thyrrenian one occurring in peninsular Italy, Elba, Corsica and Sardinia and a third one restricted to Sicily and Marettimo. Do the Iberian wood mice belong to one of these lineages? In order to answer this question, animals were trapped all over the Iberian peninsula (11 sites) as well as in three of the Balearic islands. Comparisons with specimens from the above mentioned lineages were made. From 158 animals trapped in 30 sites, 78 different mtDNA restriction patterns were obtained and compared using the NEI and LI index ((Nei M and Li WH, (1979). Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. USA 76: 5269-5273) of nucleotide divergence. A neighbour-joining tree, with a specimen of Apodemus flavicollis as outgroup, was then computed from the similarity matrix. Apodemus flavicollis is well separated from all the A. sylvaticus and all the Iberian restriction patterns are clustered with those of continental France, showing a great similarity level between the Iberian and north-west European animals. This group is well separated from the Sicilian and from the Tyrrhenian ones. These results suggest that the Pyrenees are not a biogeographic barrier for the wood mouse and that the postglacial recolonisation of western Europe by that species has its origin in populations which, during the latest Ice age, were living in refuges situated in southern France or in the Iberian peninsula. Because the Balearic restriction patterns are clustered in a lineage joining Iberian pattern at a low divergence level, we suggest that their origin is continental Spain or southern France.

Mediterranean Europe as an area of endemism for small mammals rather than a source for northwards postglacial colonization

There is a general perception that central and northern Europe were colonized by range expansion from Mediterranean refugia at the end of the last glaciation. Data from various species support this scenario, but we question its universality. Our mitochondrial DNA studies on three widespread species of small mammal suggest that colonization may have occurred from glacial refugia in central Europe-western Asia. The haplotypes on the Mediterranean peninsulae are distinctive from those found elsewhere. Rather than contributing to the postglacial colonization of Europe, Mediterranean populations of widespread small mammals may represent long-term isolates undergoing allopatric speciation. This could explain the high endemism of small mammals associated with the Mediterranean peninsulae.

Comparative phylogeography and postglacial colonization routes in Europe

The Quaternary cold periods in Europe are thought to have heavily influenced the amount and distribution of intraspecific genetic variation in both animals and plants. The phylogeographies of 10 taxa, including mammals (Ursus arctos, Sorex spp., Crocidura suaveolens, Arvicola spp.), amphibians (Triturus spp.), arthropods (Chorthippus parallelus), and plants (Abies alba, Picea abies, Fagus sylvatica, Quercus spp.), were analysed to elucidate general trends across Europe. Only a small degree of congruence was found amongst the phylogeographies of the 10 taxa, but the likely postglacial colonization routes exhibit some similarities. A Brooks parsimony analysis produced an unrooted area phylogram, showing that: (i) the northern regions were colonized generally from the Iberic and Balkanic refugia; and (ii) the Italian lineages were often isolated due to the presence of the Alpine barrier. The comparison of colonization routes highlighted four main suture-zones where lineages from the different refugia meet. Some of the intraspecific genetic distances among lineages indicated a prequaternary divergence that cannot be connected to any particular cold period, but are probably related mainly to the date of arrival of each taxon in the European continent. As a consequence, molecular genetics so far appears to be of limited use in dating Quaternary events.

Mitochondrial DNA polymorphism, phylogeography, and conservation genetics of the brown bear Ursus arctos in Europe

Some small European populations of the brown bear (Ursus arctos) are threatened by the risk of extinction in the near future. The reinforcement of these populations with bears from other regions might provide a solution to their future survival. However, before any population transfer, the different conservation units must be identified. The phylogeographic approach has been advocated for this purpose. The different European populations were assayed for mitochondrial (mt) DNA polymorphism. A remarkable degree of concordance was found between the geographic distribution and the mtDNA haplotypes. Two clearly distinct lineages differing by more than 7% in mtDNA control region sequences were found and, furthermore, the western lineage appears to be organized into two clades which correspond to two different ancestral refugia. The potential conservation units can be deduced from these results, and a management policy can consequently be inferred. This study clearly demonstrates the relevance of the molecular phylogeographic approach to the identification of conservation units.