Molecular evidence for a founder effect in invasive house finch (Carpodacus mexicanus) populations experiencing an emergent disease epidemic

The impact of founder events on levels of genetic variation in natural populations remains a topic of significant interest. Well-documented introductions provide a valuable opportunity to examine how founder events influence genetic diversity in invasive species. House finches (Carpodacus mexicanus) are passerine birds native to western North America, with the large eastern North American population derived from a small number of captive individuals released in the 1940s. Previous comparisons using amplified fragment length polymorphism (AFLP) markers found equivalent levels of diversity in eastern and western populations, suggesting that any genetic effects of the founder event were ameliorated by the rapid growth of the newly established population. We used an alternative marker system, 10 highly polymorphic microsatellites, to compare levels of genetic diversity between four native and five introduced house finch populations. In contrast to the AFLP comparisons, we found significantly lower allelic richness and heterozygosity in introduced populations across all loci. Three out of five introduced populations showed significant reductions in the ratio of the number of alleles to the allele size range, a within-population characteristic of recent bottlenecks. Finally, native and introduced populations showed significant pairwise differences in allele frequencies in every case, with stronger isolation by distance within the introduced than native range. Overall, our results provide compelling molecular evidence for a founder effect during the introduction of eastern house finches that reduced diversity levels at polymorphic microsatellite loci and may have contributed to the emergence of the Mycoplasma epidemic which recently swept the eastern range of this species.     

Allozyme, mitochondrial-DNA, and morphometrie variability indicate cryptic species of anchovy (Engraulis encrasicolus)

Previous surveys of population structure in the Atlantic-Mediterranean anchovy Engraulis encrasicolus L. have reported heterogeneity in morphology, allozyme frequencies, and mitochondrial DNA haplotype frequencies at a regional scale. In particular, two stocks of anchovy have been detected in the Adriatic Sea. In this paper, the available data is reviewed with the aim to relate genetic variation to geography at the widest possible geographical scale, for investigating the evolutionary mechanisms underlying stock structure in anchovy. Correspondence analysis of allozyme frequencies (24 samples, three polymorphic loci) compiled from the literature indicates three distinct entities in the Mediterranean Sea. Open-sea or oceanic anchovy populations are genetically different from inshore-water populations within a region (Nei's ^ G _ST = 0.035–0.067), while broadscale geographical variation is weak for each of these two habitat-specific forms (^ G _ST = 0.005–0.006). Mitochondrial-DNA haplotype frequencies support the distinction between an inshore form and an oceanic form (^ G _ST = 0.067–0.107), with virtually no genetic differences among oceanic populations across the Gulf of Biscay, the western Mediterranean and the Ionian Sea (^ G _ST = −0.001). If natural selection on marker loci is unimportant, these results indicate the occurrence of two parapatric, genetically distinct, habitat-specific forms that are widely distributed throughout the Mediterranean Sea. Persistent allele and haplotype-frequency differences between these forms indicate reproductive isolation and the presence of an E. encrasicolus species complex in the Mediterranean. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society 75 : 261–269.     

Insights from 180 years of mitochondrial variability in the endangered Mediterranean monk seal (Monachus monachus)

Mediterranean monk seals (MMS) are among the most endangered marine mammals on Earth. We screened mitochondrial variability (control region [CR1] and mitogenomes) of the species through a 180‐yr timeframe and extended by 20% (n = 205) the number of samples from a previous investigation, including historical specimens from 1833 to 1975. Although we detected two new, rare CR1 haplotypes, genetic diversity remained extremely low. Fully resolved haplotype median network and rarefaction analysis both suggested low probability for further unscreened haplotypes. There was no clear phylogeographic structure across the 12 marine subdivisions covered by the species’ range. Haplotypes previously considered diagnostic of the extant North Atlantic and eastern Mediterranean populations had their distributions extended into the western Mediterranean and the North Atlantic, respectively, by both historical and recent samples. Our study suggests that MMS have been genetically depauperate since at least the mid‐19th century, and that the massive 1997 die‐off in Western Sahara (North Atlantic) could have caused local haplotype extinctions. Our results support the hypothesis of past metapopulation dynamics across the species range, where the current segregation into geographically distant and genetically depauperate breeding populations (i.e., North Atlantic and eastern Mediterranean Sea) derives from the combined effects of historical extinctions, genetic drift on small breeding groups, and persistently low levels of genetic diversity.     

Genetic diversity and gene flow between the wild olive (oleaster, Olea europaea L.) and the olive: several Plio-Pleistocene refuge zones in the Mediterranean basin suggested by simple sequence repeats analysis

Aim The oleaster is believed to have originated in the eastern Mediterranean, implying that those in the western Mediterranean basin could be feral. Several studies with different molecular markers (isozymes, random amplified polymorphic DNA, amplified fragment length polymorphism) have shown a cline between the eastern and the western populations, which supports this hypothesis. To reconstruct the post‐glacial colonization history and establish a relationship between olive and oleaster populations in the Mediterranean basin, analyses were carried out on the genetic variation of chloroplast DNA (chlorotype) and at 12 unlinked simple sequence repeat (SSR) loci, sampling a total of 20 oleaster groves. Location This is the first known large‐scale molecular study of SSR loci based on samples of both oleasters and cultivars from the entire Mediterranean basin. Methods Samples were taken from 166 oleasters in 20 groves of modern populations, and 40 cultivars to represent molecular diversity in the cultivated olive. The Bayesian method and admixture analysis were used to construct the ancestral populations (RPOP; reconstructed panmictic oleaster populations) and to estimate the proportion of each RPOP in each tree. If one tree can be assigned to two or more RPOPs, it can be regarded as a product of hybridization between trees from different populations (i.e. admix origin). Results On this first examination of the SSR genetic diversity in the olive and oleaster, it was found to be structured in seven RPOPs in both eastern and western populations. Based on different population genetic methods, it was shown that: (1) oleasters are equally present in the eastern and the western Mediterranean, (2) are native, and (3) are not derived from cultivars. Chlorotypes (one and three in the eastern and western Mediterranean, respectively) revealed fruit displacement for the oleasters. Main conclusions Oleaster genetic diversity is divided into seven regions that could overlay glacial refuges. The gradient, or cline, of genetic diversity revealed by chloroplast and SSR molecular markers was explained by oleaster recolonization of the Mediterranean basin from refuges after the last glacial event, located in both eastern and western regions. It is likely that gene flow has occurred in oleasters mediated by cultivars spread by human migration or through trade. Animals may have helped spread oleasters locally, but humans have probably transported olives but not oleaster fruits over long distances. We found that cultivars may have originated in several RPOPs, and thus, some may have a more complex origin than expected initially.     

A microsatellite-based estimation of clonal diversity and population subdivision in Zostera marina, a marine flowering plant

We examined the genetic population structure in eelgrass (Zostera marina L.), the dominant seagrass species of the northern hemisphere, over spatial scales from 12 km to 10 000 km using the polymorphism of DNA microsatellites. Twelve populations were genotyped for six loci representing a total of 67 alleles. Populations sampled included the North Sea (four), the Baltic Sea (three), the western Atlantic (two), the eastern Atlantic (one), the Mediterranean Sea (one) and the eastern Pacific (one). Microsatellites revealed substantial genetic variation in a plant group with low allozyme diversity. Average expected heterozygosities per population (monoclonal populations excluded) ranged from 0.32 to 0.61 (mean = 0. 48) and allele numbers varied between 3.3 and 6.7 (mean = 4.7). Using the expected frequency of multilocus genotypes within populations, we distinguished ramets from genetic individuals (i.e. equivalent to clones). Differences in clonal diversity among populations varied widely and ranged from maximal diversity (i.e. all ramets with different genotype) to near or total monoclonality (two populations). All multiple sampled ramets were excluded from further analysis of genetic differentiation within and between populations. All but one population were in Hardy-Weinberg equilibrium, indicating that Zostera marina is predominantly outcrossing. From a regression of the pairwise population differentiation with distance, we obtained an effective population size Ne of 2440-5000. The overall genetic differentiation among eelgrass populations, assessed as rho (a standardized estimate of Slatkin's RST) was 0.384 (95% CI 0.34-0.44, P < 0.001). Genetic differentiation was weak among three North Sea populations situated 12-42 km distant from one another, suggesting that tidal currents result in an efficient exchange of propagules. In the Baltic and in Nova Scotia, a small but statistically significant fraction of the genetic variance was distributed between populations (rho = 0.029-0. 053) at scales of 15-35 km. Pairwise genetic differentiation between European populations were correlated with distance between populations up to a distance of 4500 km (linear differentiation-by-distance model, R2 = 0.67). In contrast, both Nova Scotian populations were genetically much closer to North Sea and Baltic populations than expected from their geographical distance (pairwise rho = 0.03-0.08, P < 0.01). A biogeographical cluster of Canadian with Baltic/North Sea populations was also supported using a neighbour-joining tree based on Cavalli-Sforza's chord distance. Relatedness between populations may be very different from predictions based on geographical vicinity.     

A Genome-Wide Survey of Genetic Variation in Gorillas Using Reduced Representation Sequencing

All non-human great apes are endangered in the wild, and it is therefore important to gain an understanding of their demography and genetic diversity. Whole genome assembly projects have provided an invaluable foundation for understanding genetics in all four genera, but to date genetic studies of multiple individuals within great ape species have largely been confined to mitochondrial DNA and a small number of other loci. Here, we present a genome-wide survey of genetic variation in gorillas using a reduced representation sequencing approach, focusing on the two lowland subspecies. We identify 3,006,670 polymorphic sites in 14 individuals: 12 western lowland gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla beringei graueri). We find that the two species are genetically distinct, based on levels of heterozygosity and patterns of allele sharing. Focusing on the western lowland population, we observe evidence for population substructure, and a deficit of rare genetic variants suggesting a recent episode of population contraction. In western lowland gorillas, there is an elevation of variation towards telomeres and centromeres on the chromosomal scale. On a finer scale, we find substantial variation in genetic diversity, including a marked reduction close to the major histocompatibility locus, perhaps indicative of recent strong selection there. These findings suggest that despite their maintaining an overall level of genetic diversity equal to or greater than that of humans, population decline, perhaps associated with disease, has been a significant factor in recent and long-term pressures on wild gorilla populations.     

Genetic variation and evolutionary history of holly oak: a circum-Mediterranean species-complex [Quercus coccifera L./Q. calliprinos (Webb) Holmboe, Fagaceae]

Holly oak is the only evergreen oak possessing a circum-Mediterranean range; it has two predominant morphological forms, calliprinos and coccifera, described in the eastern and western Mediterranean Basin respectively. The concordance of allozyme and morphotype variation was analysed in the whole holly oak range, and the most plausible historical events responsible for the current geographic pattern of genetic variation were investigated. Individuals from 24 populations were scored for allozyme variation at seven polymorphic loci. Multi-locus genotypes were analysed by using a correspondence analysis (CA) and a Bayesian clustering approach. The relative positions of the populations were obtained from multi-dimensional scaling coupled with UPGMA treatment. A continuous genotype distribution was observed in the CA, and two groups were identified using the Bayesian approach. With a 0.95 threshold, 66 and 69% of the individuals showing the calliprinos and the coccifera morphotypes respectively were assigned to the corresponding genetic groups, which differed by private alleles. As compared to coccifera, the genetic group calliprinos was characterized by higher allelic richness and a strong geographical genetic structure. In agreement with fossil records, the most parsimonious explanation for lack of geographical structure in coccifera is a substantial regression of holly oak in the western Mediterranean Basin during the glaciations and a fast westward post-glacial expansion of coccifera populations, probably from Greece. Two population groups were obtained from the scaling/UPGMA treatment. One included all of the calliprinos populations and a Greek coccifera population, suggesting that the two morphotypes are closely related genetically and constitute two components of the same species.     

Allozyme variation in cork oak (Quercus suber L.): the role of phylogeography and genetic introgression by other Mediterranean oak species and human activities

 Genetic variation in the cork oak (Quercus suber L.) was investigated using 11 loci from seven enzyme systems in 40 populations sampled over the entire distribution of this species in the western Mediterranean Basin. Mean heterozygosity values over the polymorphic loci (Ho=0.283), the percentage of polymorphic populations (M=0.76), and the total genetic diversity (Ht=0.31) from which 11% was accounted for among-population variation, are among the highest recorded in oak species. In contrast to previous results in holm oak (Q. ilex L.), another evergreen species in the same area, cork oak possessed a smaller allele pool and a lower average number of alleles per locus and per population (A=2.0). More particularly, very few low-frequency alleles were observed in cork oak except for eight populations in which allozyme polymorphism at locus Pgi 1, diagnostic between both species, indicates that these low-frequency alleles are introgressed from holm oak. On the basis of the genetic distance estimated from allozyme frequencies, 32 of the 40 cork oak populations studied were classified into two very distinct sets which also corresponded to distinct geographic areas. One set gathered together the 18 populations from the Iberian peninsula and two adjacent areas in France, i.e. the centre of origin of cork oak, according to paleobotanical data. This set was characterized by a larger allele pool, a higher within-population genetic diversity and a lower differentiation between populations than was observed in the other set, which comprised the populations from North Africa, Sicily, Sardinia, Corsica, continental Italy and the region of Provence (southeastern France). In these more southern and eastern disjunct areas, cork oak migration from Iberia may have occurred at different periods since the end of the Tertiary. The possible effect of human activity on cork oak genetic structure, i.e. the selection of good-quality cork, acorn over-use for animal food, and even human nutrition, is discussed. Received: 3 March 1998 / Accepted: 19 March 1998     

Microsatellite markers reveal shallow genetic differentiation between cohorts of the common sea urchin Paracentrotus lividus (Lamarck) in northwest Mediterranean

Temporal variability was studied in the common sea urchin Paracentrotus lividus through the analysis of the genetic composition of three yearly cohorts sampled over two consecutive springs in a locality in northwestern Mediterranean. Individuals were aged using growth ring patterns observed in tests and samples were genotyped for five microsatellite loci. No reduction of genetic diversity was observed relative to a sample of the adult population from the same location or within cohorts across years. F _ST and amova results indicated that the differentiation between cohorts is rather shallow and not significant, as most variability is found within cohorts and within individuals. This mild differentiation translated into estimates of effective population size of 90–100 individuals. When the observed excess of homozygotes was taken into account, the estimate of the average number of breeders increased to c . 300 individuals. Given our restricted sampling area and the known small-scale heterogeneity in recruitment in this species, our results suggest that at stretches of a few kilometres of shoreline, large numbers of progenitors are likely to contribute to the larval pool at each reproduction event. Intercohort variation in our samples is six times smaller than spatial variation between adults of four localities in the western Mediterranean. Our results indicate that, notwithstanding the stochastic events that take place during the long planktonic phase and during the settlement and recruitment processes, reproductive success in this species is high enough to produce cohorts genetically diverse and with little differentiation between them. Further research is needed before the link between genetic structure and underlying physical and biological processes can be well established.     

Demographic Genetics of Brown Trout (Salmo Trutta) and Estimation of Effective Population Size from Temporal Change of Allele Frequencies

We studied temporal allele frequency shifts over 15 years and estimated the genetically effective size of four natural populations of brown trout (Salmo trutta L.) on the basis of the variation at 14 polymorphic allozyme loci. The allele frequency differences between consecutive cohorts were significant in all four populations. There were no indications of natural selection, and we conclude that random genetic drift is the most likely cause of temporal allele frequency shifts at the loci examined. Effective population sizes were estimated from observed allele frequency shifts among cohorts, taking into consideration the demographic characteristics of each population. The estimated effective sizes of the four populations range from 52 to 480 individuals, and we conclude that the effective size of natural brown trout populations may differ considerably among lakes that are similar in size and other apparent characteristics. In spite of their different effective sizes all four populations have similar levels of genetic variation (average heterozygosity) indicating that excessive loss of genetic variability has been retarded, most likely because of gene flow among neighboring populations.     

Development of microsatellite markers for the critically endangered Limonium dufourii (Girard) Kuntze (Plumbaginaceae)

Limonium dufourii is an endemic plant from the eastern Mediterranean coast of Spain with a triploid chromosome number and apomictic reproduction. We have isolated and characterized 13 polymorphic microsatellite loci from an enriched library in order to investigate its population genetic structure. Simple sequence repeat (SSR) loci were screened in 120 individuals from the six extant populations of this species. They show an average of 5.76 alleles per locus, ranging from 2 to 18, with seven loci exhibiting heterozygosities larger than 0.60. Three loci present one single allele in each individual, whereas one locus presents three alleles in every individual analysed.     

Microsatellite analysis of population structure in Canadian polar bears

Attempts to study the genetic population structure of large mammals are often hampered by the low levels of genetic variation observed in these species. Polar bears have particularly low levels of genetic variation with the result that their genetic population structure has been intractable. We describe the use of eight hypervariable microsatellite loci to study the genetic relationships between four Canadian polar bear populations: the northern Beaufort Sea, southern Beaufort Sea, western Hudson Bay, and Davis Strait - Labrador Sea. These markers detected considerable genetic variation, with average heterozygosity near 60% within each population. Interpopulation differences in allele frequency distribution were significant between all pairs of populations, including two adjacent populations in the Beaufort Sea. Measures of genetic distance reflect the geographic distribution of populations, but also suggest patterns of gene flow which are not obvious from geography and may reflect movement patterns of these animals. Distribution of variation is sufficiently different between the Beaufort Sea populations and the two more eastern ones that the region of origin for a given sample can be predicted based on its expected genotype frequency using an assignment test. These data indicate that gene flow between local populations is restricted despite the long-distance seasonal movements undertaken by polar bears.     

Differential population structuring of two closely related fish species, the mackerel (Scomber scombrus) and the chub mackerel (Scomber japonicus), in the Mediterranean Sea

Population genetic structures of the mackerel (Scomber scombrus) and chub mackerel (Scomber japonicus) were studied in the Mediterranean Sea. Fragments of 272 bp (S. scomber) and 387 bp (S. japonicus) of the 5'-end of the mitochondrial control region were sequenced from spawning individuals collected off the coasts of Greece, Italy, Spain, and Portugal. High levels of mitochondrial control region haplotypic diversity (> 0.98) were found for both Scomber species. Nucleotide diversity was higher in the mackerel (0.022) than in the chub mackerel (0.017). Global F(ST) values were also higher and significant in the mackerel (0.024, P < 0.0001) as opposed to the chub mackerel (0.003, P > 0.05). Molecular variance analyses showed differential genetic structuring for these two closely related species. There is extensive gene flow between Mediterranean Sea and Atlantic Ocean populations of chub mackerel, which are organized into a larger panmictic unit. In contrast, Mediterranean Sea populations of mackerel show some degree of genetic differentiation and are structured along an east-west axis. The analysed eastern Mediterranean Sea mackerel populations (Greece, Italy) are clearly separated from that of the western Mediterranean Sea (Barcelona), which forms a panmictic unit with eastern Atlantic Ocean populations. The genetic structures of both species showed asymmetric migration patterns and indicated population expansion.     

Population genetic structure of the endangered crayfish Austropotamobius pallipes in France based on microsatellite variation: biogeographical inferences and conservation implications

1. One important goal in conservation biology is to characterise evolutionary lineages within endangered species before management decisions are taken. Here, we assess population differentiation in the freshwater crayfish Austropotamobius pallipes, an endangered species endemic to western Europe and provide valuable information for the conservation of French populations. 2. Analysis of five microsatellite loci in 44 populations revealed very different within population levels of genetic diversity (0.000 < H₀ < 0.564). Two groups, corresponding to northern and southern French populations, showed a high degree of genetic differentiation in both allele frequencies and allele sizes. Comparison of these results with previous studies of A. pallipes strongly suggests that the divergence between northern and southern populations could have occurred during the last glaciation period of the Pleistocene from one Atlantic and one Mediterranean refuge. 3. Evidence for genetic admixture between these two lineages was revealed by correspondence analyses in southern populations, probably as the result of artificial translocations. 4. French populations appeared significantly differentiated among the different river drainages and were highly structured within rivers. The impact of population size, population bottlenecks and founder events on the population genetic differentiation are discussed. 5. Based on these results, we propose the designation of two evolutionarily significant units for A. pallipes in France. Our data also support the maintenance of separate demographic management strategies for crayfish inhabiting different river systems. However, genetic analyses will have to be combined with demographic and ecological data for sustainable conservation programmes.     

Genetic variation at nine short tandem repeat loci among islanders of the eastern Adriatic coast of Croatia

We have analyzed the extent of genetic variation at nine autosomal short tandem repeat loci (D3S1358, VWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, D7S820) among six populations from Croatia: five distributed in the islands of the eastern Adriatic coast and one from the mainland. The purpose is to investigate the usefulness of these loci in detecting regional genetic differentiation in the studied populations. Significant heterogeneity among the island and mainland populations is revealed in the distributions of allele frequencies; however, the absolute magnitude of the coefficient of gene differentiation is small but significant. The summary measures of genetic variation, namely, heterozygosity, number of alleles, and allele size variance, do not indicate reduced genetic variation in the island populations compared to the mainland population. In contrast to the two measures of genetic variation, allele size variance and within-locus heterozygosity, the imbalance index (beta) indicates evidence of recent expansion of population sizes in all islands and in the mainland. High mutation rates of the studied loci together with local drift effects are likely explanations for interisland genetic variation and the observed lack of reduced genetic diversity among the island populations.     

Around the Mediterranean: an extreme example of loop migration in a long‐distance migratory passerine

An important issue in migration research is how small‐bodied passerines pass over vast geographical barriers; in European–African avian migration, these are represented by the Mediterranean Sea and the Sahara Desert. Eastern (passing eastern Mediterranean), central (passing Apennine Peninsula) and western (via western Mediterranean) major migration flyways are distinguished for European migratory birds. The autumn and spring migration routes may differ (loop migration) and there could be a certain level of individual flexibility in how individuals navigate themselves during a single migration cycle. We used light‐level loggers to map migration routes of barn swallows Hirundo rustica breeding in the centre of a wide putative contact zone between the northeastern and southernwestern European populations that differ in migration flyways utilised and wintering grounds. Our data documented high variation in migration patterns and wintering sites of tracked birds (n = 19 individuals) from a single breeding colony, with evidence for loop migration in all but one of the tracked swallows. In general, two migratory strategies were distinguished. In the first, birds wintering in a belt stretching from southcentral to southern Africa that used an eastern route for both the spring and autumn migration, then shifted their spring migration eastwards (anti‐clockwise loops, n = 12). In the second, birds used an eastern or central route to their wintering grounds in central Africa, shifting the spring migration route westward (clockwise loops, n = 7). In addition, we observed an extremely wide clockwise loop migration encompassing the entire Mediterranean, with one individual utilising both the eastern (autumn) and western (spring) migratory flyway during a single annual migration cycle. Further investigation is needed to ascertain whether clockwise migratory loops encircling the entire Mediterranean also occur other small long‐distance passerine species.     

Conservation genetics of the short-beaked common dolphin (<Emphasis Type="Italic">Delphinus delphis</Emphasis>) in the Mediterranean Sea and in the eastern North Atlantic Ocean

Mediterranean Sea common dolphins have recently been listed as ‘endangered’ in the IUCN Red list, due to their reported decline since the middle of the 20th century. However, little is know about the number or distribution of populations in this region. We analysed 118 samples from the Black Sea, Mediterranean Sea and eastern North Atlantic at nine microsatellite nuclear loci and for 428 bps of the mtDNA control region. We found small but significant population differentiation across the basin between the eastern and the western Mediterranean populations at both nuclear and mtDNA markers (microsatellite F _ST = 0.052, mtDNA F _ST = 0.107, P values ≤ 0.001). This matched the differential distribution and habitat use patterns exhibited by this species in the eastern and the western parts of the Mediterranean Sea. The assignment test of a small number of samples from the central Mediterranean could not exclude further population structure in the central area of the basin. No significant genetic differentiation at either marker was observed among the eastern north Atlantic populations, though the Alboran population (inhabiting the Mediterranean waters immediately adjacent the Atlantic ocean) showed significant mtDNA genetic differentiation compared to the Atlantic populations. Directional estimates of gene flow suggested movement of females out of the Mediterranean, which may be relevant to the population decline. Phylogenetic analysis suggested that the observed population structure evolved recently.     

Genetic diversity and phylogeography of the European roe deer: the refuge area theory revisited

The extant taxa of central and northern Europe are commonly believed to derive from Pleistocene ancestors, who moved to the north from three separate glacial refugia: the Iberian and Italian peninsulae, as well as the southern Balkans. The issue of postglacial dispersal patterns was addressed through the investigation of population structure and phylogeography of the European roe deer, Capreolus capreolus . The genetic diversity in 376 individuals representing 14 allegedly native populations across their European range was assessed, using ten autosomal microsatellite loci and restriction fragment length polymorphisms of the mitochondrial D-loop and NADH dehydrogenase 1 gene segments. Our results suggest the existence of three major genetic lineages of roe deer in Europe. One comprises populations in the south-western limit of the species' distribution (i.e. Iberia), where an internal substructure splits a northern from a southern sublineage. A second lineage includes populations of southern and eastern Europe, as well as a separate sublineage sampled in central-southern Italy, where the existence of the subspecies Capreolus c. italicus was supported. In central-northern Europe, a third lineage is present, which appeared genetically rather homogeneous, although admixed, and equally divergent from both the eastern and western lineages. Current patterns of intraspecific genetic variation suggest that postglacial recolonization routes of this cervid to northern Europe could be due to range expansion from one or more refugia in central-eastern Europe, rather than proceeding from the Mediterranean areas.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 85–100.