A single introduction of wild rabbits triggered the biological invasion of Australia.

Biological invasions are a major cause of environmental and economic disruption. While ecological factors are key determinants of their success, the role of genetics has been more challenging to demonstrate. The colonization of Australia by the European rabbit is one of the most iconic and devastating biological invasions in recorded history. Here, we show that despite numerous introductions over a 70-y period, this invasion was triggered by a single release of a few animals that spread thousands of kilometers across the continent. We found genetic support for historical accounts that these were English rabbits imported in 1859 by a settler named Thomas Austin and traced the origin of the invasive population back to his birthplace in England. We also find evidence of additional introductions that established local populations but have not spread geographically. Combining genomic and historical data we show that, contrary to the earlier introductions, which consisted mostly of domestic animals, the invasive rabbits had wild ancestry. In New Zealand and Tasmania, rabbits also became a pest several decades after being introduced. We argue that the common denominator of these invasions was the arrival of a new genotype that was better adapted to the natural environment. These findings demonstrate how the genetic composition of invasive individuals can determine the success of an introduction and provide a mechanism by which multiple introductions can be required for a biological invasion.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil.

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

A loss-of-function mutation in RORB disrupts saltatorial locomotion in rabbits.

Saltatorial locomotion is a type of hopping gait that in mammals can be found in rabbits, hares, kangaroos, and some species of rodents. The molecular mechanisms that control and fine-tune the formation of this type of gait are unknown. Here, we take advantage of one strain of domesticated rabbits, the sauteur d'Alfort, that exhibits an abnormal locomotion behavior defined by the loss of the typical jumping that characterizes wild-type rabbits. Strikingly, individuals from this strain frequently adopt a bipedal gait using their front legs. Using a combination of experimental crosses and whole genome sequencing, we show that a single locus containing the RAR related orphan receptor B gene (RORB) explains the atypical gait of these rabbits. We found that a splice-site mutation in an evolutionary conserved site of RORB results in several aberrant transcript isoforms incorporating intronic sequence. This mutation leads to a drastic reduction of RORB-positive neurons in the spinal cord, as well as defects in differentiation of populations of spinal cord interneurons. Our results show that RORB function is required for the performance of saltatorial locomotion in rabbits.

African climate and geomorphology drive evolution and ghost introgression in sable antelope.

The evolutionary history of African ungulates has been explained largely in the light of Pleistocene climatic oscillations and the way these influenced the distribution of vegetation types, leading to range expansions and/or isolation in refugia. In contrast, comparatively fewer studies have addressed the continent's environmental heterogeneity and the role played by its geomorphological barriers. In this study, we performed a range-wide analysis of complete mitogenomes of sable antelope (Hippotragus niger) to explore how these different factors may have contributed as drivers of evolution in southcentral Africa. Our results supported two sympatric and deeply divergent mitochondrial lineages in west Tanzanian sables, which can be explained as the result of introgressive hybridization of a mitochondrial ghost lineage from an archaic, as-yet-undefined, congener. Phylogeographical subdivisions into three main lineages suggest that sable diversification may not have been driven solely by climatic events affecting populations differently across a continental scale. Often in interplay with climate, geomorphological features have also clearly shaped the species' patterns of vicariance, where the East Africa Rift System and the Eastern Arc Mountains acted as geological barriers. Subsequent splits among southern populations may be linked to rearrangements in the Zambezi system, possibly framing the most recent time when the river attained its current drainage profile. This work underlines how the use of comprehensive mitogenomic data sets on a model species with a wide geographical distribution can contribute to a much-enhanced understanding of environmental, geomorphological and evolutionary patterns in Africa throughout the Quaternary.

Changes in brain architecture are consistent with altered fear processing in domestic rabbits.

The most characteristic feature of domestic animals is their change in behavior associated with selection for tameness. Here we show, using high-resolution brain magnetic resonance imaging in wild and domestic rabbits, that domestication reduced amygdala volume and enlarged medial prefrontal cortex volume, supporting that areas driving fear have lost volume while areas modulating negative affect have gained volume during domestication. In contrast to the localized gray matter alterations, white matter anisotropy was reduced in the corona radiata, corpus callosum, and the subcortical white matter. This suggests a compromised white matter structural integrity in projection and association fibers affecting both afferent and efferent neural flow, consistent with reduced neural processing. We propose that compared with their wild ancestors, domestic rabbits are less fearful and have an attenuated flight response because of these changes in brain architecture.

Unravelling population processes over the Late Pleistocene driving contemporary genetic divergence in Palearctic buzzards.

Population range expansions and contractions as a response to climate and habitat change throughout the Quaternary are known to have contributed to complex phylogenetic and population genetic events. Speciation patterns and processes in Palearctic buzzards (genus Buteo) are a long-standing example of morphological and genetic data incongruence, attributed to panmixia, habitat range shifts, contact zones, and climate change. Here we assess the systematics, phylogeography and population genetic structure of three nominal species of Palearctic buzzards, Buteo buteo (including B. b. vulpinus), B. rufinus (including B. r. cirtensis) and B. hemilasius. Phylogenetic analyses inferred from mitochondrial data recover B. hemilasius as sister to the sister clades B. r. rufinus and B. buteo complex (B. b. buteo, B. b. vulpinus, but also including B. r. cirtensis). In contrast, we find an unresolved genetic delimitation inferred from four nuclear loci, suggesting an ancestral genetic pool for all species. Time-trees suggest population contractions and expansions throughout the Pleistocene, which likely reflect habitat change and contrasting ecological niche requirements between species. Microsatellite-based extended Bayesian skyline plots reveal relatively constant population sizes for B. hemilasius, B. r. rufinus, and B. b. vulpinus, in contrast to a dramatic population expansion in B. r. cirtensis within the last 3 kya. Overall, our study illustrates how complex population processes over the Late Pleistocene have shaped the patterns of genetic divergence in Palearctic buzzards, due to the joint effects of shared ancestral polymorphisms, population expansions and contractions, with hybridization at contact zones leading to admixture and introgression.

A genomic map of clinal variation across the European rabbit hybrid zone.

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

The genomic architecture of population divergence between subspecies of the European rabbit.

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

Understanding the mechanisms of antitropical divergence in the seabird White-faced Storm-petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach.

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White-faced Storm-petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial ΦST = 0.89 (P < 0.01) and global nuclear ΦST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long-term antitropical divergence in isolation. Coalescent-based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000-1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000-600 000), both probably associated with the profound palaeo-oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long-distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.

Molecular evidence for cryptic candidate species in Iberian Pelodytes (Anura, Pelodytidae).

Species delineation is a central topic in evolutionary biology, with current efforts focused on developing efficient analytical tools to extract the most information from molecular data and provide objective and repeatable results. In this paper we use a multilocus dataset (mtDNA and two nuclear markers) in a geographically comprehensive population sample across Iberia and Western Europe to delineate candidate species in a morphologically cryptic species group, Parsley frogs (genus Pelodytes). Pelodytes is the sole extant representative of an ancient, historically widely distributed anuran clade that currently includes three species: P. caucasicus in the Caucasus; P. punctatus in Western Europe, from Portugal to North-Western Italy; and P. ibericus in Southern Iberia. Phylogenetic analyses recovered four major well-supported haplotype clades in Western Europe, corresponding to well demarcated geographical subdivisions and exhibiting contrasting demographic histories. Splitting times date back to the Plio-Pleistocene and are very close in time. Species-tree analyses recovered one of these species lineages, corresponding to P. ibericus (lineage B), as the sister taxon to the other three major species lineages, distributed respectively in: western Iberian Peninsula, along the Atlantic coast and part of central Portugal (lineage A); Central and Eastern Spain (lineage C); and North-eastern Spain, France and North-western Italy (lineage D). The latter is in turn subdivided into two sub-clades, one in SE France and NW Italy and the other one from NE Spain to NW France, suggesting the existence of a Mediterranean-Atlantic corridor along the Garonne river. An information theory-based validation approach implemented in SpedeSTEM supports an arrangement of four candidate species, suggesting the need for a taxonomic revision of Western European Pelodytes.

Candidate genes underlying heritable differences in reproductive seasonality between wild and domestic rabbits.

Reproductive seasonality is a trait that often differs between domestic animals and their wild ancestors, with domestic animals showing prolonged or even continuous breeding seasons. However, the genetic basis underlying this trait is still poorly understood for most species, and because environmental factors and resource availability are known to play an important role in determining breeding seasons, it is also not clear in most cases to what extent this phenotypic shift is determined by the more lenient captive conditions or by genetic factors. Here, using animals resulting from an initial cross between wild and domestic rabbits followed by two consecutive backcrosses (BC1 and BC2) to wild rabbits, we evaluated the yearly distribution of births for the different generations. Similar to domestic rabbits, F1 animals could be bred all year round but BC1 and BC2 animals showed a progressive and significant reduction in the span of the breeding season, providing experimental evidence that reduced seasonal breeding in domestic rabbits has a clear genetic component and is not a simple by-product of rearing conditions. We then took advantage of a recently published genome-wide scan of selection in the domesticated lineage and searched for candidate genes potentially associated with this phenotypic shift. Candidate genes located within regions targeted by selection include well-known examples of genes controlling clock functions (CRY1 and NR3C1) and reproduction (PRLR).

A comparison of brain gene expression levels in domesticated and wild animals.

Domestication has led to similar changes in morphology and behavior in several animal species, raising the question whether similarities between different domestication events also exist at the molecular level. We used mRNA sequencing to analyze genome-wide gene expression patterns in brain frontal cortex in three pairs of domesticated and wild species (dogs and wolves, pigs and wild boars, and domesticated and wild rabbits). We compared the expression differences with those between domesticated guinea pigs and a distant wild relative (Cavia aperea) as well as between two lines of rats selected for tameness or aggression towards humans. There were few gene expression differences between domesticated and wild dogs, pigs, and rabbits (30-75 genes (less than 1%) of expressed genes were differentially expressed), while guinea pigs and C. aperea differed more strongly. Almost no overlap was found between the genes with differential expression in the different domestication events. In addition, joint analyses of all domesticated and wild samples provided only suggestive evidence for the existence of a small group of genes that changed their expression in a similar fashion in different domesticated species. The most extreme of these shared expression changes include up-regulation in domesticates of SOX6 and PROM1, two modulators of brain development. There was almost no overlap between gene expression in domesticated animals and the tame and aggressive rats. However, two of the genes with the strongest expression differences between the rats (DLL3 and DHDH) were located in a genomic region associated with tameness and aggression, suggesting a role in influencing tameness. In summary, the majority of brain gene expression changes in domesticated animals are specific to the given domestication event, suggesting that the causative variants of behavioral domestication traits may likewise be different.

Multilocus assessment of phylogenetic relationships in Alytes (Anura, Alytidae).

With the advent of large multilocus datasets, molecular systematics is experiencing very rapid progress, but important challenges remain regarding data analysis and interpretation. Midwife toads (genus Alytes) exemplify two of the most widespread problems for accurate phylogenetic reconstruction: discerning the causes of discordance between gene trees, and resolving short internodes produced during rapid, successive splitting events. The three species in subgenus Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), the sister group to A. obstetricans, have disjunct and highly restricted geographical ranges, which are thought to result from old vicariant events affecting their common ancestor, but their phylogenetic relationships are still unresolved. In this study we re-address the phylogeny of Alytes with a special focus on the relationships in Baleaphryne with a multilocus dataset including >9000 base pairs of mitochondrial DNA and four nuclear markers (3142bp) in all recognized taxa, including all subspecies of A. obstetricans. Both concatenation and species tree analyses suggest that A. muletensis, endemic to the Balearic island of Mallorca, is the sister taxon to a clade comprising the southeastern Iberian endemic A. dickhilleni and the North African A. maurus. This scenario is consistent with palaeogeological evidence associated with the fragmentation of the Betic-Rifean Massif, followed by the opening of the Strait of Gibraltar. On the other hand, analyses of intraspecific variation in A. obstetricans are inconclusive regarding relationships between major clades and conflict with current subspecific taxonomy.

Selection against domestication alleles in introduced rabbit populations.

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome.

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit.

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X-chromosome, Y-chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95-4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co-localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X- and Y-chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.

Gene-culture coevolution between cattle milk protein genes and human lactase genes.

Milk from domestic cows has been a valuable food source for over 8,000 years, especially in lactose-tolerant human societies that exploit dairy breeds. We studied geographic patterns of variation in genes encoding the six most important milk proteins in 70 native European cattle breeds. We found substantial geographic coincidence between high diversity in cattle milk genes, locations of the European Neolithic cattle farming sites (>5,000 years ago) and present-day lactose tolerance in Europeans. This suggests a gene-culture coevolution between cattle and humans.

Assessment of genetic diversity among native Algerian rabbit populations using microsatellite markers.

Having higher adaptability against abiotic stress, which is characterized in rural areas in developing countries, local farm animal genetic resources (FAGRs) are increasingly precarious for random and unsystematic crossing with exotic breeds. In this study, 85 microsatellite loci were utilized to assess genetic diversity among native Algerian rabbits (NARs) sampled from an area of 753 km (from north to south) and 919 km (from east to west). Those distances covered 25 significant geographical points in seven rural areas (El Taref, Mostaganem, Sidi Bel Abbès, M'Sila, Dar Chioukh, Faidh El Botma, and Laghouat). A total of 558 alleles were observed in this study. The highest genetic diversity was registered in the southern direction among NAR populations. The mean number of alleles per locus (MNa) and the inbreeding coefficient (FIS) were highest in Laghouat (4.482 and 0.232), while they were lowest in El Taref (4.000 and 0.149). In the current study, the number of private alleles (Pa) ranged from 9 to 23. In addition, the average of observed heterozygosity (0.427) was lower than the expected value (0.524) due to high levels of inbreeding. The discriminant analysis of principal components (DAPC), the neighbor-joining tree (NJ), and the analysis of STRUCTURE software confirmed the classification of populations according to geographical zones into four main groups (east, west, south, and middle). The results of the current study are useful for breeding improvement and conservation plan research in relation to local animal genetic resources in Algeria.

The genetic structure of domestic rabbits.

Understanding the genetic structure of domestic species provides a window into the process of domestication and motivates the design of studies aimed at making links between genotype and phenotype. Rabbits exhibit exceptional phenotypic diversity, are of great commercial value, and serve as important animal models in biomedical research. Here, we provide the first comprehensive survey of nucleotide polymorphism and linkage disequilibrium (LD) within and among rabbit breeds. We resequenced 16 genomic regions in population samples of both wild and domestic rabbits and additional 35 fragments in 150 rabbits representing six commonly used breeds. Patterns of genetic variation suggest a single origin of domestication in wild populations from France, supporting historical records that place rabbit domestication in French monasteries. Levels of nucleotide diversity both within and among breeds were ~0.2%, but only 60% of the diversity present in wild populations from France was captured by domestic rabbits. Despite the recent origin of most breeds, levels of population differentiation were high (F(ST) = 17.9%), but the majority of polymorphisms were shared and thus transferable among breeds. Coalescent simulations suggest that domestication began with a small founding population of less than 1,200 individuals. Taking into account the complex demographic history of domestication with two successive bottlenecks, two loci showed deviations that were consistent with artificial selection, including GPC4, which is known to be associated with growth rates in humans. Levels of diversity were not significantly different between autosomal and X-linked loci, providing no evidence for differential contributions of males and females to the domesticated gene pool. The structure of LD differed substantially within and among breeds. Within breeds, LD extends over large genomic distances. Markers separated by 400 kb typically showed r(2) higher than 0.2, and some LD extended up to 3,200 kb. Much less LD was found among breeds. This advantageous LD structure holds great promise for reducing the interval of association in future mapping studies.

Postglacial colonization of Europe by the barbastelle bat: agreement between molecular data and past predictive modelling.

The barbastelle (Barbastella barbastellus) is a rare forest bat with a wide distribution in Europe. Here, we combine results from the analysis of two mtDNA fragments with species distribution modelling to determine glacial refugia and postglacial colonization routes. We also investigated whether niche conservatism occurs in this species. Glacial refugia were identified in the three southern European peninsulas: Iberia, Italy and the Balkans. These latter two refugia played a major role in the postglacial colonization process, with their populations expanding to England and central Europe, respectively. Palaeo-distribution models predicted that suitable climatic conditions existed in the inferred refugia during the last glacial maximum (LGM). Nevertheless, the overlap between the current and the LGM distributions was almost inexistent in Italy and in the Balkans, meaning that B. barbastellus populations were forced to shift range between glacial and interglacial periods, a process that probably caused some local extinctions. In contrast, Iberian populations showed a 'refugia within refugium' pattern, with two unconnected areas containing stable populations (populations that subsisted during both glacial and interglacial phases). Moreover, the match between LGM models and the refugial areas determined by molecular analysis supported the hypothesis of niche conservatism in B. barbastellus. We argue that geographic patterns of genetic structuring, altogether with the modelling results, indicate the existence of four management units for conservation: Morocco, Iberia, Italy and UK, and Balkans and central Europe. In addition, all countries sampled possessed unique gene pools, thus stressing the need for the conservation of local populations.