Differentiation at mitochondrial and nuclear loci between the blesbok (Damaliscus pygargus phillipsi) and bontebok (D. p. pygargus): implications for conservation strategy

Damaliscus pygargus is an endemic species in South Africa belonging to the contemporary antelope tribe Alcelaphine. The species has been subdivided into two subspecies based upon phenotypic differences and historical geographic isolation. We found D. p. pygargus (bontebok) and D. p. phillipsi (blesbok) to be significantly differentiated based on microsatellite and mtDNA control region markers. Bontebok genetic diversity was depauperate with just one mtDNA haplotype, not found in blesbok populations, and reduced heterozygosity and polymorphism at microsatellite DNA loci compared to the blesbok. Erosion of molecular genetic variation in bontebok may have been due to anthropogenic impact causing isolation and reduced population size. Our data based on 34 bontebok and 42 blesbok indicate that the classification of alpha taxonomy should be reconsidered in this genus (in support of earlier similar suggestions), and that careful management should seek to avoid hybridization and retain the remaining diversity in the bontebok.     

Characterization of major histocompatibility complex DRB diversity in the endemic South African antelope Damaliscus pygargus: a comparison in two subspecies with different demographic histories

Major histocompatibility complex (MHC) class II locus DRB was investigated by single-strand conformation polymorphism analysis (SSCP) and sequence analysis in the endemic South African antelope, Damaliscus pygargus, of which there are two subspecies. Greater polymorphism was found in the blesbok (D. p. phillipsi) subspecies (n = 44; 22 alleles) than in the bontebok (D. p. pygargus) subspecies (n = 45; 6 alleles). Erosion of allelic diversity in bontebok was most likely the result of two severe bottleneck events caused by hunting pressure and parasitic infection. A majority of the polymorphism observed was found within the peptide binding region (PBR) where dN/dS ratios were higher than for the non-PBR region. This, and the apparent trans-species relationship among alleles in a bovid phylogeny, suggest the evolution of diversity by heterosis or frequency-dependent selection.     

Quantitative evaluation of hybridization and the impact on biodiversity conservation

Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.     

Chromosome homologies between tsessebe (Damaliscus lunatus) and Chinese muntjac (Muntiacus reevesi) facilitate tracing the evolutionary history of Damaliscus (Bovidae, Antilopinae, Alcelaphini)

Genome-wide homologies between the tsessebe (Damaliscus lunatus, 2n = 36) and Chinese muntjac (Muntiacus reevesi, 2n = 46) have been established by cross-species painting with Chinese muntjac chromosome paints. Twenty-two autosomal painting probes detected 35 orthologous segments in the tsessebe. Hybridization results confirmed that: (i) D. lunatus carries the (9;14) reciprocal translocation that has been proposed to be a derived chromosomal landmark shared by all species of the Antilopinae; (ii) the karyotype of D. lunatus can be derived almost exclusively from the bovid ancestral karyotype through 12 Robertsonian translocations involving 24 ancestral acrocentric autosomes; (iii) in addition to the Rb fusions, pericentric heterochromatic amplification has shaped the morphology of several of the D. lunatus chromosomes. Integrated analysis of these and published cytogenetic data on pecorans has allowed us to accurately discern the karyotype history of Damaliscus (D. lunatus; D. pygargus, 2n = 38; D. hunteri, 2n = 44). The phylogenomic relationships of 3 species reflected by specific chromosomal rearrangements were consistent with published phylogenies based on morphology, suggesting that chromosomal rearrangements have played an important role in speciation within the Alcelaphini, and that karyotype characters are valuable phylogenetic markers in this group.     

Ten species specific microsatellite loci for Lahontan cutthroat trout,Oncorhynchus clarki henshawi

Ten polymorphic microsatellite loci (containing tri and tetra‐nucleotide repeats) were developed for the Lahontan cutthroat trout (Oncorhynchus clarki henshawi), a subspecies of cutthroat trout listed as threatened under the United States Endangered Species Act. Polymorphism was assessed for 445 individuals from 12 populations representing eight watersheds spread throughout the three Distinct Population Segments defined for this subspecies. All loci were polymorphic (X? = 19, range 7–26 alleles). All loci were in Hardy–Weinberg equilibrium (HWE) except for one locus (OCH 9) in a single population (P < 0.00014 after Bonferroni correction for multiple tests).     

Genetic diversity and population divergence in fragmented habitats: Conservation of Idaho ground squirrels

The Idaho ground squirrel, which consists of a northern (Spermophilus brunneus brunneus) and a southern subspecies (S. b. endemicus), has suffered from habitat loss and fragmentation, resulting in a reduction in both numbers and geographic range of the species. The northern Idaho ground squirrel (NIDGS) is listed as a threatened subspecies under the Endangered Species Act, and the southern Idaho ground squirrel (SIDGS) is a candidate. Because Idaho ground squirrel populations are small and often isolated, they are susceptible to inbreeding and loss of genetic diversity through drift. This research evaluates levels of genetic diversity and patterns of population divergence in both subspecies of Idaho ground squirrels. We hypothesized that NIDGS would exhibit lower genetic diversity and greater population divergence due to a longer period of population isolation relative to most SIDGS populations. Genetic diversity and divergence were quantified using 8 microsatellite loci. Contrary to expectations, SIDGS populations exhibited consistently lower levels of microsatellite diversity. Additionally, NIDGS exhibited only modest divergence among populations, while divergence levels among SIDGS populations were highly varied. Preliminary evaluations of mitochondrial DNA diversity and structure revealed lower diversity in NIDGS and some differences in gene flow that warrant further study. Based on our results, we suggest different management strategies for the two subspecies. Habitat restoration appears to be the most desirable conservation strategy for NIDGS populations. In contrast, low genetic diversity observed in SIDGS may warrant supplementation of isolated populations through translocations or captive breeding to mitigate further loss of genetic variability.     

Abundance,genetic diversity and conservation of Louisiana black bears (<Emphasis Type="Italic">Ursus americanus luteolus</Emphasis>) as detected through noninvasive sampling

Although the Louisiana black bear (Ursus americanus luteolus) is currently listed as threatened under the Endangered Species Act, there have been no attempts to estimate range-wide abundance. This subspecies was thought to occupy a near contiguous range across southern Mississippi, Louisiana and east Texas but is now restricted to three isolated areas in Louisiana. In 1964, Louisiana initiated a restocking program in which black bears from Minnesota were introduced into two of these areas. It is not clear how the additions affected population structure or if substantial breeding occurred between native and introduced bears. Using baited sites to snare hair samples, and microsatellite DNA analysis to distinguish individuals, we estimated abundance of two geographically isolated bear populations in south central Louisiana: Inland and Coastal. Additionally, we examined genetic variation both within and between the two populations. Mark recapture analysis of the distribution of individual captures during two primary sampling periods resulted in population estimates of 77 ± 9 for Coastal and 41 ± 6 for Inland. Genetic analysis revealed significant population differentiation (F _ST = 0.206) between the two populations. The apparently smaller Inland population exhibited more diversity than the Coastal, which suggests that the genetic structure of the Inland population has been influenced by the reintroduction. Both of these populations are isolated and face considerable demographic and genetic threats, thus conservation measures to protect both are warranted. However, the Coastal population is more representative of Louisiana black bears prior to reintroduction and special consideration should be given to insure its integrity.     

Quantitative analysis of genetic parameters in populations of European (Capreolus capreolus L.) and Siberian (Capreolus pygargus Pall.) roe deer with RAPD markers

A technique for obtaining unbiased estimates of genetic parameters (allelic frequencies of RAPD loci, heterozygosity (H), Wright's F statistic, and Nei's genetic distances) in populations of the European (Capreolus capreolus L.) and Siberian (Capreolus pygargus Pall.) roe deer is presented. The technique employs jackknifing and multiple comparative analysis based on a modified Holmes's procedure for Bonferroni's test. It was demonstrated that samples from local groups of roe deer in the Trans-Ural region did not differ significantly in allelic frequencies (0.8, 0.81, and 0.78; P > 0.447) or Nei's genetic distances (0.0056, 0.0273, and 0.0218; P = 0.26), but they could be differentiated based on Wright's F statistic (0.0346, 0.0519, and 0.0450; P = 10(-9)). The parameters of intrapopulation heterozygosity (from 0.18 to 0.042) formed a gradient from the east to the west. Calibration estimates of molecular evolution rate in the family Cervidae obtained based on published data and Jukes-Cantor genetic distances estimated in this study demonstrated that the Siberian roe deer has split into two subspecies, C. pygargus pygargus Pall. and C. pygargus tianschanicus Satunin in the interval between 229 and 462.3 thousand years ago. The species formation of the Siberian and European roe deer was dated between 1.375 and 2.75 Myr ago. Based on the results obtained we recommend the approaches used in the study for analysis of population genetic structure and phylogenetic relationships between populations, subspecies, species, and higher taxa.     

Assessing introgressive hybridization between blue wildebeest (<Emphasis Type="Italic">Connochaetes taurinus</Emphasis>) and black wildebeest (<Emphasis Type="Italic">Connochaetes gnou</Emphasis>) from South Africa

Introgressive hybridization poses a threat to the genetic integrity of black wildebeest (Connochaetes gnou) and blue wildebeest (Connochaetes taurinus) populations in South Africa. Black wildebeest is endemic to South Africa and was driven to near extinction in the early 1900s due to habitat destruction, hunting pressure and disease outbreaks. Blue wildebeest on the other hand are widely distributed in southern and east Africa. In South Africa the natural distribution ranges of both species overlap, however, extensive translocation of black wildebeest outside of its normal distribution range in South Africa have led to potential hybridization between the two species. The molecular identification of pure and admixed populations is necessary to design viable and sustainable conservation strategies, since phenotypic evidence of hybridization is inconclusive after successive generations of backcrossing. The aim of this study was to assess levels of hybridization in wildebeest using both species-specific and cross-species microsatellite markers. Black wildebeest (157) and blue wildebeest (122) from provincial and national parks and private localities were included as reference material, with 180 putative hybrid animals also screened. A molecular marker panel consisting of 13 cross-species and 11 species-specific microsatellite markers was developed. We used a Bayesian clustering model to confirm the uniqueness of blue- and black wildebeest reference groups, assign individuals to each of the two clusters, and determine levels of admixture. Results indicated a clear partition between black wildebeest and blue wildebeest (the average proportions of membership to black wildebeest and blue wildebeest clusters were QI?=?0.994 and QI?=?0.955 respectively). From the putative hybrid samples, only five hybrid individuals were confirmed. However, high levels of linkage disequilibrium were observed in the putative hybrid populations which may indicate historical hybridization. Measures of genetic diversity in the black wildebeest populations were found to be lower than that of the blue wildebeest. The observed lower level of genetic diversity was expected due to the demographic history of the specie. This study will make a significant contribution to inform a national conservation strategy to conserve the genetic integrity of both species.     

Reduced SNP Panels for Genetic Identification and Introgression Analysis in the Dark Honey Bee (Apis mellifera mellifera)

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham’s F_ST, an F_ST-based outlier test, Delta, informativeness (I_n), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.     

Evolution of body size in the genus Damaliscus: a comparison with hartebeest Alcelaphus spp.

In species with low levels of sexual size dimorphism, it may be relatively easy to detect the role of natural selection in the evolution of body size. Habitat primary production (HPP) appears to be a key factor in the divergence of size in the hartebeest clade ( Alcelaphus spp.), such that subspecies in less productive savannahs are smaller than those in richer ones. Here I test whether a similar pattern exists within the genus Damaliscus (topi and their allies). Basal skull length was used as a surrogate of body size and measured in the seven allopatric subspecies of Damaliscus . Means for each subspecies and sex were regressed against climatic factors as surrogates of HPP. Variation in skull length across Damaliscus taxa was less than in hartebeest. Two clusters were present in both sexes and corresponded to the distinction between the species, Damaliscus dorcas and Damaliscus lunatus . This may reflect differences in productivity between edaphic grasslands, occupied by all D. lunatus , and dry grasslands, occupied by D. dorcas . Mean annual rainfall was the best predictor of body size in males and showed a non-significant positive tendency in females. After accounting for phylogenetic effects, these correlations were both non-significant. Edaphic grasslands might be less dependent on precipitation for primary production because the impeded drainage of their soil prolongs water availability after the end of the rains. Furthermore, they are probably more consistent in productivity across African regions than secondary grasslands and savannah woodlands, which rely on rainfall for grass growth. These properties of edaphic grasslands may explain why size in Damaliscus appears to be less sensitive to variation in rainfall and less variable across subspecies than in Alcelaphus .     

Development of microsatellite markers in Gonystylus bancanus (Ramin) useful for tracing and tracking of wood of this protected species

Ten polymorphic microsatellite markers have been developed for Gonystylus bancanus (Ramin), a protected tree species of peat swamp forests in Malaysia and Indonesia. Eight markers were also shown to be polymorphic in other Gonystylus species. The markers will enable assessing the amount of genetic variation within and among populations and the degree of population differentiation, such that donor populations can be selected for reforestation projects. They may be used for tracing and tracking of wood in the production chain, so that legal trade in this Convention on International Trade in Endangered Species of Wild Fauna and Flora-protected timber species, derived from specifically described origins, can be distinguished from illegally logged timber.     

Intraspecific comparison of population structure, genetic diversity, and dispersal among three subspecies of Townsend’s big-eared bats, Corynorhinus townsendii townsendii, C. t. pallescens, and the endangered C. t. virginianus

Townsend’s big-eared bat, Corynorhinus townsendii, is distributed broadly across western North America and in two isolated, endangered populations in central and eastern United States. There are five subspecies of C. townsendii; C. t. pallescens, C. t. australis, C. t. townsendii, C. t. ingens, and C. t. virginianus with varying degrees of concern over the conservation status of each. The aim of this study was to use mitochondrial and microsatellite DNA data to examine genetic diversity, population differentiation, and dispersal of three C. townsendii subspecies. C. t. virginianus is found in isolated populations in the eastern United States and was listed as endangered under the Endangered Species Act in 1979. Concern also exists about declining populations of two western subspecies, C. t. pallescens and C. t. townsendii. Using a comparative approach, estimates of the genetic diversity within populations of the endangered subspecies, C. t. virginianus, were found to be significantly lower than within populations of the two western subspecies. Further, both classes of molecular markers revealed significant differentiation among regional populations of C. t. virginianus with most genetic diversity distributed among populations. Genetic diversity was not significantly different between C. t. townsendii and C. t. pallescens. Some populations of C. t. townsendii are not genetically differentiated from populations of C. t. pallescens in areas of sympatry. For the western subspecies gene flow appears to occur primarily through male dispersal. Finally, geographic regions representing significantly differentiated and genetically unique populations of C. townsendii virginianus are recognized as distinct evolutionary significant units.     

Population genetics and conservation of the threatened southeastern beach mouse (Peromyscus polionotus niveiventris): subspecies and evolutionary units

We investigated genetic diversity within the southeastern beach mouse (SEBM-Peromyscus polionotus niveiventris) and also tested the hypothesis that the subspecies recognition of P.p. niveiventris, based on size and color differences, is congruent with this taxon representing a discrete evolutionary lineage. We used ten polymorphic microsatellite loci and mitochondrial cytochrome-b gene DNA sequences to investigate genetic diversity and population structure within the SEBM, and to determine the level of divergence between the SEBM and the nearest known inland subspecies of the oldfield mouse (Peromyscus polionotus rhoadsi). Moderate genetic distances were observed between the SEBM and the inland oldfield mouse based on microsatellite data, with F _ST values ranging from 0.11 to 0.22 between these taxa. Additionally, mitochondrial DNA haplotypes of the SEBM formed a distinct monophyletic group relative to haplotypes sampled from P. p. rhoadsi. Based on previous estimates of rates of mitochondrial DNA evolution in rodents, we inferred that Pleistocene sea-level fluctuations are likely responsible for the historical isolation of the SEBM lineage from mainland P. polionotus. Our data demonstrate the genetic distinctiveness of the SEBM, justifying the current subspecies designation for the SEBM and its continued protection under the United States Endangered Species Act. We classify the Cape Canaveral and Smyrna Dunes Park populations of SEBM as a single evolutionary significant unit. The two known extant allopatric populations of the SEBM showed some differentiation in microsatellite frequencies and were moderately reciprocally distinguishable based on assignment to distinct genetic clusters by a Bayesian admixture procedure. These results justify the classification of these two extant SEBM populations as distinct management units that should be independent targets of management and conservation attention.     

Variability of multilocus DNA markers in populations of the Siberian (Capreolus pygargus Pall.) and European (C. capreolus L.) roe deer

Two types of nuclear DNA markers, M13 minisatellites and RAPD, were used to examine intraspecific and interspecific variation in closely related roe deer species, Capreolus capreolus L. and C. pygargus Pall. The roe deer populations studied were highly polymorphic for minisatellite DNA markers (S = 0.12-0.36). Heterozygosities of the RAPD loci were 0.185 (Russian Far East), 0.145, 0.131, 0.088 (Cis-Ural), and 0.06 (France). They correlated with karyotypic variation of B chromosomes (r = 0.975, P < 0.02; Spearman's correlation coefficient r = 1, P = 0.1 x 10(-5)), which indicated a contribution of microchromosomes to genetic variation of the species. The genetic distance D between the closely related species C. capreolus and C. pygargus was 50 times greater than the distance between populations within a species. The estimates of heterozygosity and genetic distance between local populations of Cis-Ural and the Far East suggest their specific spatial organization within this geographical range and reveal features of their historical development.     

Subspecific affinity of black bears in the White River National Wildlife Refuge.

The black bear population of the White River National Wildlife Refuge (NWR) is adjacent to populations of black bear in Louisiana (Urusus americanus luteolus) which are listed as threatened under the U.S. Endangered Species Act. Wildlife management plans can pose restrictions on bear harvests and timber extraction; therefore the management plan for the White River NWR is sensitive to subspecific classification of its bear population. The objective of this study was to analyze genetic variation in the White River NWR and seven adjacent populations of black bears to assess the subspecific affinity of the White River NWR population. Here we report the variation at seven microsatellite DNA loci among eight black bear populations. The patterns of genetic variation gave strong support for distinguishing a southern group of black bears comprised of the White River, Arkansas; Tensas River, Louisiana; Upper Atchafalaya, Louisiana; Lower Atchafalaya, Louisiana; and Alabama/Mississippi populations. Phylogenetic analysis of individual variation suggested that historical black bear introductions into Arkansas and Louisiana affected gene pools of certain southern receiving populations, but did not significantly change interpopulation relatedness. Phylogenetic inferences at both the population and individual levels support the hypothesis that the White River NWR population of black bears belongs to the U. a. luteolus subspecies.     

Comprehensive genetic analyses reveal evolutionary distinction of a mouse (Zapus hudsonius preblei) proposed for delisting from the US Endangered Species Act

Zapus hudsonius preblei, listed as threatened under the US Endangered Species Act (ESA), is one of 12 recognized subspecies of meadow jumping mice found in North America. Recent morphometric and phylogenetic comparisons among Z. h. preblei and neighbouring conspecifics questioned the taxonomic status of selected subspecies, resulting in a proposal to delist the Z. h. preblei from the ESA. We present additional analyses of the phylogeographic structure within Z. hudsonius that calls into question previously published data (and conclusions) and confirms the original taxonomic designations. A survey of 21 microsatellite DNA loci and 1380 base pairs from two mitochondrial DNA (mtDNA) regions (control region and cytochrome b) revealed that each Z. hudsonius subspecies is genetically distinct. These data do not support the null hypothesis of a homogeneous gene pool among the five subspecies found within the southwestern portion of the species' range. The magnitude of the observed differentiation was considerable and supported by significant findings for nearly every statistical comparison made, regardless of the genome or the taxa under consideration. Structuring of nuclear multilocus genotypes and subspecies-specific mtDNA haplotypes corresponded directly with the disjunct distributions of the subspecies investigated. Given the level of correspondence between the observed genetic population structure and previously proposed taxonomic classification of subspecies (based on the geographic separation and surveys of morphological variation), we conclude that the nominal subspecies surveyed in this study do not warrant synonymy, as has been proposed for Z. h. preblei, Z. h. campestris, and Z. h. intermedius.     

Gene flow in admixed populations and implications for the conservation of the Western honeybee, Apis mellifera

Anthropogenic activity, especially modern apiculture, has considerable impact on the natural distribution of the Western honeybee, Apis mellifera, leading to the spread, replacement and fragmentation of many subspecies. This creates demand for the conservation of some subspecies, in particular, Apis mellifera mellifera, which once was widely distributed in Western Europe and nowadays is endangered through habitat loss and fragmentation. Moreover, A. m. mellifera may be further endangered by hybridisation in populations that now occur in artificial sympatry with other subspecies. Here, we quantify and compare individual hybridisation between sympatric and allopatric honeybee populations of A. m. mellifera and A. m. carnica using microsatellite markers and a Bayesian model-based approach. We had a special focus on pure breeding populations, which are a major tool in honeybee conservation. Our results demonstrate that subspecies are still highly differentiated, but gene flow is not prevented by the current management strategies, creating urgent demand for an improved conservation management of A. m. mellifera. However, the occurrence of a high number of pure individuals might suggest that some sort of hybrid barrier acts against the complete admixture of the two subspecies.     

Characterization of ten microsatellite marker loci in the Komodo Monitor (<Emphasis Type="Italic">Varanus komodoensis</Emphasis>)

The Komodo monitor (Varanus komodoensis) is a classic example of a species that has been fragmented into small isolated populations over a limited range. This species, classified as endangered under Appendix I of the Convention on the International Trade of Endangered Species. We describe 10 novel species-specific microsatellite loci characterized in representatives from three of the endemic island populations. One locus, 12HDZ820 appears to be fixed in one population at an allele size not found in the other two. This microsatellite suite should be helpful in augmenting the marker selection currently available for Komodo Monitor population studies.     

Hybrid origins of honeybees from italy (Apis mellifera ligustica) and sicily (A. m. sicula)

The genetic variability of honeybee populations Apis mellifera ligustica, in continental Italy, and of A. m. sicula, in Sicily, was investigated using nuclear (microsatellite) and mitochondrial markers. Six populations (236 individual bees) and 17 populations (664 colonies) were, respectively, analysed using eight microsatellite loci and DraI restriction fragment length polymorphism (RFLP) of the cytochrome oxidase I (COI)-cytochrome oxidase II (COII) region. Microsatellite loci globally confirmed the southeastern European heritage of both subspecies (evolutionary branch C). However, A. m. ligustica mitochondrial DNA (mtDNA) appeared to be a composite of the two European (M and C) lineages over most of the Italian peninsula, and only mitotypes from the African (A) lineage were found in A. m. sicula samples. This demonstrates a hybrid origin for both subspecies. For A. m. ligustica, the most widely exported subspecies, this hybrid origin has long been obscured by the fact that in the main area of queen production (from which most of the previous ligustica bee samples originated) the M mitochondrial lineage is absent, whereas it is present almost everywhere else in Italy. This presents a new view of the evolutionary history of European honeybees. For instance, the Iberian peninsula was considered as the unique refuge for the M branch during the quaternary ice periods. Our results show that the Apennine peninsula played a similar role. The differential distribution of nuclear and mitochondrial markers observed in Italy seems to be a general feature of introgressed honeybee populations. Presumably, it stems from the social nature of the species in which both genome compartments are differentially affected by the two (individual and colonial) reproduction levels.