Geographical overlap of two mitochondrial genomes in Spanish honeybees (Apis mellifera iberica)

Restriction enzyme cleavage maps of mitochondrial DNA from the Spanish honeybee, Apis mellifera iberica (Hymenoptera: Apidae), were compared with those from the European subspecies A. m. mellifera, A. m. ligustica, and A. m. carnica, and the African subspecies A. m. intermissa and A. m. scutellata. The mitochondrial DNA (mtDNA) of the two African subspecies can be distinguished by restriction fragment polymorphisms revealed by Hinf I digests. Two distinct mtDNA types were found among Spanish honeybees: a west European mellifera-like type, which predominates in the north of Spain, and an African intermissa-like type, which predominates in the south. Spain appears to be a region of contact and hybridization between the two subspecies A. m. intermissa and A. m. mellifera, which respectively represent African and west European honeybee lineages. This natural boundary between European and African honeybee populations in the Old World may provide a model for predicting the eventual outcome of the colonization of North America by introduced African honeybees.     

Varying degrees of Apis mellifera ligustica introgression in protected populations of the black honeybee, Apis mellifera mellifera, in northwest Europe

The natural distribution of honeybee subspecies in Europe has been significantly affected by human activities during the last century. Non-native subspecies of honeybees have been introduced and propagated, so that native black honeybee (Apis mellifera mellifera) populations lost their identity by gene-flow or went extinct. After previous studies investigated the remaining gene-pools of native honeybees in France and Spain, we here assess the genetic composition of eight northwest European populations of the black honeybee, using both mitochondrial (restriction fragment length polymorphisms of the intergenic transfer RNAleu-COII region) and nuclear (11 microsatellite loci) markers. Both data sets show that A. m. mellifera populations still exist in Norway, Sweden, Denmark, England, Scotland and Ireland, but that they are threatened by gene flow from commercial honeybees. Both Bayesian admixture analysis of the microsatellite data and DraI-RFLP (restriction fragment length polymorphism) analysis of the intergenic region indicated that gene-flow had hardly occurred in some populations, whereas almost 10% introgression was observed in other populations. The most introgressed population was found on the Danish Island of Laeso, which is the last remaining native Danish population of A. m. mellifera and the only one of the eight investigated populations that is protected by law. We discuss how individual admixture analysis can be used to monitor the restoration of honeybee populations that suffer from unwanted hybridization with non-native subspecies.     

The impact of apiculture on the genetic structure of wild honeybee populations (<Emphasis Type="Italic">Apis mellifera</Emphasis>) in Sudan

Apiculture often relies on the importation of non-native honeybees (Apis mellifera) and large distance migratory beekeeping. These activities can cause biodiversity conflicts with the conservation of wild endemic honeybee subspecies. We studied the impact of large scale honeybee imports on managed and wild honeybee populations in Sudan, a centre of biodiversity of A. mellifera, using as set of linked microsatellite DNA loci and mitochondrial DNA markers. The mitochondrial DNA analyses showed that all wild honey bees exclusively belonged to African haplotypes. However, we revealed non-native haplotypes in managed colonies on apiaries reflecting unambiguous evidence of imports from European stock. Moreover, we found significantly higher linkage disequilibria for microsatellite markers in wild populations in regions with apiculture compared to wild populations which had no contact to beekeeping. Introgression of imported honeybees was only measurable at the population level in close vicinity to apicultural activities but not in wild populations which represent the vast majority of honeybees in Sudan.     

中国意大利蜜蜂微卫星遗传多态性

为了分析我国不同生态条件下意大利蜜蜂Apis mellifera ligustica（简称意蜂）的基因特征,对品种间基因交流或基因渗入现象进行监测,完成对我国主要地区意蜂多样性水平的评价。本研究利用10对微卫星（SSR）引物分析了来自我国黑龙江等14个省市自治区31个采集点的意蜂样品的遗传多态性,并与来自德国的喀尼阿兰蜂A. m. carnica（简称喀蜂）和来自美国的意蜂A. m. ligustica进行比对。通过NTSYS-pc2.1和popGene32软件对SSR数据进行多态性分析。结果表明：10对SSR引物中8对表现出较高的遗传稳定性,可作为意蜂遗传多态性标记。遗传一致度在0.61以上时,浙江、四川、甘肃和云南地区的样品出现较严重分化,其他地区的意蜂样品则保持了美国原种意蜂的典型基因特征。东北和新疆地区的样品与德国喀蜂血统亲缘关系较近（遗传距离在0.18149以内）,分析主要原因可能是该地区的意蜂样品杂有黑蜂A. m. mellifera血统。据此提出应在保护各地意蜂原种特性及种群多态性的基础上,科学有效地进行高产品种的繁育。     

Evolutionary history of the honey bee Apis mellifera inferred from mitochondrial DNA analysis

Variability of mitochondrial DNA (mtDNA) of the honey bee Apis mellifera L. has been investigated by restriction and sequence analyses on a sample of 68 colonies from ten different subspecies. The 19 mtDNA types detected are clustered in three major phylogenetic lineages. These clades correspond well to three groups of populations with distinct geographical distributions: branch A for African subspecies (intermissa, monticola, scutellata, andansonii and capensis), branch C for North Mediterranean subspecies (caucasica, carnica and ligustica) and branch M for the West European populations (mellifera subspecies). These results partially confirm previous hypotheses based on morphometrical and allozymic studies, the main difference concerning North African populations, now assigned to branch A instead of branch M. The pattern of spatial structuring suggests the Middle East as the centre of dispersion of the species, in accordance with the geographic areas of the other species of the same genus. Based on a conservative 2% divergence rate per Myr, the separation of the three branches has been dated at about 1 Myr BP.     

Microsatellite Variation in Honey Bee (Apis Mellifera L.) Populations: Hierarchical Genetic Structure and Test of the Infinite Allele and Stepwise Mutation Models

Samples from nine populations belonging to three African (intermissa, scutellata and capensis) and four European (mellifera, ligustica, carnica and cecropia) Apis mellifera subspecies were scored for seven microsatellite loci. A large amount of genetic variation (between seven and 30 alleles per locus) was detected. Average heterozygosity and average number of alleles were significantly higher in African than in European subspecies, in agreement with larger effective population sizes in Africa. Microsatellite analyses confirmed that A. mellifera evolved in three distinct and deeply differentiated lineages previously detected by morphological and mitochondrial DNA studies. Dendrogram analysis of workers from a given population indicated that super-sisters cluster together when using a sufficient number of microsatellite data whereas half-sisters do not. An index of classification was derived to summarize the clustering of different taxonomic levels in large phylogenetic trees based on individual genotypes. Finally, individual population X loci data were used to test the adequacy of the two alternative mutation models, the infinite allele model (IAM) and the stepwise mutation models. The better fit overall of the IAM probably results from the majority of the microsatellites used including repeats of two or three different length motifs (compound microsatellites).     

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M‐lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C‐lineage). Effective conservation actions require reliable molecular tools to identify pure‐bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q‐values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.     

Morphological and microsatellite differentiation in Melospiza melodia (Aves) at a microgeographic scale

Geographical variation in microsatellite allele frequencies and morphology were compared for five subspecies of Melospiza melodia (song sparrow; M. m. samuelis, M. m. maxillaris, M. m. pusillula, M. m. gouldii, and M. m. heermanni) in 14 populations in the San Francisco Bay region to (a) assess divergence based on these estimates and (b) test the hypothesis that drift is responsible for morphological and genetic divergence. Morphological differentiation between subspecies was high despite low differentiation at microsatellite loci, indicating high gene flow and large effective population sizes. Low concordance of morphological and genetic estimates of divergence suggests that selection or phenotypic plasticity in morphology has caused morphological differentiation among the subspecies.     

Phylogeography and subspecies taxonomy of dunlins (Calidris alpina) in western Palearctic analysed by DNA microsatellites and amplified fragment length polymorphism markers

In a circumpolar wader, the dunlin ( Calidris alpina ), there are 11 named subspecies, but only five mitochondrial DNA (mtDNA) lineages have been found. In the present study, we investigated the genetic structure of dunlins in western Palearctic (from East Greenland to Taimyr peninsula) using DNA microsatellites and amplified fragment length polymorphism (AFLP) markers that may detect more recent differentiation than mtDNA. In this region, we consider four described subspecies; alpina , schinzii , arctica and centralis , together comprising two mtDNA lineages. We analyse seven polymorphic microsatellite loci and 91 AFLP markers in 287 and 152 unrelated individuals, respectively, originating from 17 populations. Neither microsatellites nor AFLPs reveal distinct groups that correspond to currently recognized subspecies. There is a clear pattern of isolation by distance in microsatellites. Our results do not contradict the former mtDNA results that there are two phylogenetic lineages (approximately corresponding to schinzii and centralis ) that have met and formed a cline ( alpina ). We find no difference between schinzii and arctica (East Greenland). We conclude that, given the lack of distinct groups and the gradual changes in microsatellite allele frequencies, these markers provide little genetic support for the dunlin subspecies taxonomy in the western Palearctic.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 713–726.     

Mitochondrial Sequence Characterisation of Australian Commercial and Feral Honeybee Strains, Apis mellifera L. (Hymenoptera: Apidae), in the Context of the Species Worldwide

The history of honeybee ( Apis mellifera ) importations and management in Australia is largely anecdotal and therefore a survey and characterisation of this agriculturally important insect is of value. We give information on the genetic composition of 42 feral and commercial strains by sequencing sections of the ATPase 6, cytochrome oxidase III, cytochrome b and ND2 mitochondrial genes to determine the relationship of the strains to each other. Our phylogenetic analysis shows novel associations between A. m. mellifera, A. m. scutellata, A. m. ligustica and A. m. caucasica.     

Gene flow across secondary contact zones of the Emys orbicularis complex in the Western Mediterranean and evidence for extinction and re‐introduction of pond turtles on Corsica and Sardinia (Testudines: Emydidae)

European pond turtles represent a phylogeographically deeply structured complex of distinct taxa. Here, we use mitochondrial DNA sequences (cytochrome b gene) and eight polymorphic microsatellite loci to investigate genetic differentiation and gene flow of Sicilian, Corsican and Sardinian pond turtles and of subspecies involved in two secondary contact zones in the Pyrenean region and Southern Italy. Mitochondrial and microsatellite differentiation is largely concordant in populations from the core regions of the distribution ranges of the studied taxa. Both marker systems provide no evidence for gene flow between Sicilian pond turtles (Emys trinacris) and Southern Italian subspecies of E. orbicularis. By contrast, in the contact zones limited gene flow occurs between distinct subspecies of E. orbicularis. Although the Southern Italian contact zone is significantly older than the Pyrenean contact zone of Holocene age, patterns of asymmetric introgression are similar. Introgressive hybridization leads to the exchange of mitochondria, but microsatellite data indicate only a few individuals with mixed ancestry. This suggests that incipient isolating mechanisms maintain largely discrete nuclear genomic gene pools. Furthermore, this implies that Southern Italy acted as a hotspot rather than as a melting pot of genetic diversity during the last glacial. Pond turtles from Corsica and Sardinia are not differentiated from continental populations of the subspecies E. o. galloitalica, neither in the mitochondrial nor in the quickly evolving microsatellite markers. As the fossil record argues for a continuous presence of pond turtles on both islands since the Middle Pleistocene, this suggests that the native island populations became extinct and the extant turtles were later introduced by prehistoric settlers. The lack of genetic differentiation of pond turtles from Corsica and Sardinia supports the view that the subspecies described from these islands are not valid.     

Discrimination of Iranian honeybee populations (Apis mellifera meda) from commercial subspecies of Apis mellifera L. using morphometric and genetic methods

The morphological characters of honeybees have an important role for discriminating honeybee subspecies. In the present research, Iranian populations of honeybee (Apis mellifera) were collected from 19 areas in Iran. The samples were collected from stationary beekeeping sites. Moreover, pictures of honeybee forewings held in the Bee Data Bank in Oberursel were compared with Iranian honeybee populations. 19 morphological characters were measured for each forewing of worker honeybee to evaluate differentiation of Iranian honeybee populations from the commercial honeybee subspecies A. m. mellifera, A. m. carnica, A. m. caucasica and A. m. ligustica. Additionally, part of the tRNA^leu gene, an intergenic region and part of COII was used to confirm differentiation of the commercial subspecies from Iranian honeybee populations. Results of the cluster analyses showed that 19 morphological characters of forewings differentiated Iranian populations from the commercial subspecies. Moreover, the phylogenetic tree of part of the tRNA^leu gene, an intergenic region and part of COII differentiated the commercial subspecies from Iranian honeybee populations. Results of the discriminant function analyses (DFA) indicated that the references samples of A. m. meda overlapped with Iranian populations.     

Microsatellite analysis of non-migratory colonies of Apis mellifera iberica from south-eastern Spain

Forty‐five unmanaged honeybee colonies from the south‐east of the Iberian Peninsula (Apis mellifera iberica) were selected for analysing their genetic structure using eight microsatellite loci. These colonies were not subjected to management for queen replacement, rearing or migratory movements and previous studies showed that they had mitochondrial DNA (mtDNA) of African origin. Six of the microsatellite loci show intermediate levels of polymorphism with a total number of alleles detected per locus ranging from 4 to 10. Microsatellite data relate these Iberian populations to the African A. m. intermissa, although the presence of some alleles and the observed heterozygosity are characteristic of the European A. m. mellifera, thus corroborating the postulated hybrid origin of A. m. iberica. The results suggest that no recent introgression from Africa has happened and that the populations of A. m. iberica are differentiated in many demes.     

Genetic variation across the historical range of the wild turkey (Meleagris gallopavo)

Genetic differences within and among naturally occurring populations of wild turkeys (Meleagris gallopavo) were characterized across five subspecies' historical ranges using amplified fragment length polymorphism (AFLP) analysis, microsatellite loci and mitochondrial control region sequencing. Current subspecific designations based on morphological traits were generally supported by these analyses, with the exception of the eastern (M. g. silvestris) and Florida (M. g. osceola) subspecies, which consistently formed a single unit. The Gould's subspecies was both the most genetically divergent and the least genetically diverse of the subspecies. These genetic patterns were consistent with current and historical patterns of habitat continuity. Merriam's populations showed a positive association between genetic and geographical distance, Rio Grande populations showed a weaker association and the eastern populations showed none, suggesting differing demographic forces at work in these subspecies. We recommend managing turkeys to maintain subspecies integrity, while recognizing the importance of maintaining regional population structure that may reflect important adaptive variation.     

Hybrid origin of Japanese mice "Mus musculus molossinus": evidence from restriction analysis of mitochondrial DNA

The Japanese mouse, Mus musculus molossinus, has long been considered an independent subspecies of the house mouse. A survey of restriction- site haplotypes of mitochondrial DNA (mtDNA) showed that Japanese mice have two main maternal lineages. The most common haplotype is closely related to the mtDNA of the European subspecies M. m. musculus. The other common haplotype and two minor ones are closely related to each other and to the mtDNA of an Asiatic subspecies, M. m. castaneus. Two other rare variants are probably the result of recent contamination by European M. m. domesticus. The musculus type of mtDNA is found in the southern two-thirds of Japan, whereas the common castaneus type is found in the northern third and the minor variants are found sporadically throughout Japan. The castaneus mtDNA lineage had a few minor variants, whereas the musculus lineage was completely monomorphic. By contrast, the native population of M. m. castaneus and the Chinese and Korean musculus populations were highly polymorphic. These results suggest that M. m. molossinus is a hybrid between ancestral colonies, possibly very small, of M. m. musculus and M. m. castaneus, rather than an independent subspecies.      

Molecular phylogeography of European Sciurus vulgaris: refuge within refugia?

The red squirrel ( Sciurus vulgaris ) is a well-known forest animal distributed all over Europe. Still, we are far from having a firm knowledge of the species' phylogeography. This study investigates the genetic differentiation of S. vulgaris across the species' Eurasian range, using sequence data from the mitochondrial DNA gene (D-loop, 252 base pairs, cytochrome b , 359 base pairs), and eight variable autosomal microsatellite loci genotyped for 236 individuals. The results reveal the presence of two main mitochondrial phylogroups. The first clade comprises the individuals from the region of Calabria in southern Italy, belonging to the subspecies S. v . meridionalis , while the second clade contains the remainder of the studied individuals. Bayesian analysis of microsatellite genotypes resulted in three main clusterings corresponding to the three S. vulgaris subspecies: infuscatus, meridionalis and fuscoater . Geographical distribution of mtDNA haplotypes and mismatch analysis suggest a common refugium for the red squirrel across most of its present range from which expansion happened rather rapidly. The genotype mixing of italicus with northern populations could be a residual of postglacial expansion. The lack of mixing between the Calabrian lineage and the rest of European red squirrel haplotypes can be seen as evidence for distinct histories throughout the Pleistocene. Calabrian mtDNA probably diverged in an ice age contraction and remained isolated from the neighbouring squirrel populations until very recent times.     

New approach to the mitotype classification in black honeybee <Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> and Iberian honeybee <Emphasis Type="Italic">Apis mellifera iberiensis</Emphasis>

The black honeybee Apis mellifera mellifera L. is today the only subspecies of honeybee which is suitable for commercial breeding in the climatic conditions of Northern Europe with long cold winters. The main problem of the black honeybee in Russia and European countries is the preservation of the indigenous gene pool purity, which is lost as a result of hybridization with subspecies, A. m. caucasica, A. m. carnica, A. m. carpatica, and A. m. armeniaca, introduced from southern regions. Genetic identification of the subspecies will reduce the extent of hybridization and provide the gene pool conservation of the black honeybee. Modern classification of the honeybee mitotypes is mainly based on the combined use of the DraI restriction endonuclease recognition site polymorphism and sequence polymorphism of the mtDNA COI–COII region. We performed a comparative analysis of the mtDNA COI–COII region sequence polymorphism in the honeybees of the evolutionary lineage M from Ural and West European populations of black honeybee A. m. mellifera and Spanish bee A. m. iberiensis. A new approach to the classification of the honeybee M mitotypes was suggested. Using this approach and on the basis of the seven most informative SNPs of the mtDNA COI–COII region, eight honeybee mitotype groups were identified. In addition, it is suggested that this approach will simplify the previously proposed complicated mitotype classification and will make it possible to assess the level of the mitotype diversity and to identify the mitotypes that are the most valuable for the honeybee breeding and rearing.     

Phylogeography of the threatened crayfish (genus Austropotamobius) in Italy: implications for its taxonomy and conservation

A nucleotide sequence analysis of a portion of the mitochondrial large ribosomal subunit was performed to define the phylogeography of the threatened crayfish Austropotamobius (Decapoda; Astacidae) in Italy. We collected 61 specimens from 31 localities across the Italian peninsula. For the phylogenetic inference, we combined the 61 Austropotamobius spp sequences obtained from this study with 18 sequences deposited in GenBank and corresponding to Italian, French, Irish, Swiss, and Slovenian locations. Among the analysed sequences, 34 distinct haplotypes were detected. Our results confirmed the presence of both A. pallipes and A. italicus in the Italian peninsula and the existence within the latter species of a strong intraspecific genetic variation, due to the occurrence of four subspecies with a well-defined geographic distribution. From a conservation viewpoint, Italy, with its high haplotype variability, may be considered a 'hot spot' for the genetic diversity of the European native crayfish Austropotamobius. We suggest that re-introduction programs should be conducted with extreme caution in Italy, since not only the two Austropotamobius species but also the four A. italicus subspecies are genetically and taxonomically separate units and require independent conservation plans.     

Geographical patterns of mitochondrial DNA variation in Apis mellifera iberiensis (Hymenoptera: Apidae)

An extensive survey of mitochondrial haplotypes in honeybee colonies from the Iberian Peninsula has corroborated previous hypotheses about the existence of a joint clinal variation of African (A) and west European (M) evolutionary lineages. It has been found that the Iberian Peninsula is the European region with the highest haplotype diversity (12 haplotypes detected of the M lineage and 10 of the A lineage). The frequency of A haplotypes decreases in a SW-NE trend, while that of M haplotypes increases. These results are discussed in relation to hypotheses about the African origin of Apis mellifera and an early colonization of west Europe during intermediate Pleistocene glaciation events, followed by a regional differentiation. The extant pattern of haplotype frequency and distribution seems to be influenced at a regional scale by adaptation to local climatic conditions and the mobile beekeeping that has become a large-scale practice during the last decades. Other previous anthropogenic influences (Greek, Roman and Arab colonizations) are thought to be of minor importance in present day populations.     

Fine-scale phylogeographical analysis of Mediterranean Anacamptis palustris (Orchidaceae) populations based on chloroplast minisatellite and microsatellite variation

The phylogeographical history of the rare marsh orchid Anacamptis palustris (Orchidaceae) was reconstructed using highly polymorphic chloroplast minisatellite and microsatellite loci. Allelic variation at chloroplast microsatellite loci was due to length variation in poly(A/T) repeats and was informative on a regional scale, but was not sufficient to unravel relationships among populations on a local geographical scale. The minisatellite locus, however, was found to be highly variable. Nine distinct repeat types were found and variation in repeat number occurred in five repeat types. The distribution of chloroplast haplotypes, combining microsatellite and minisatellite repeat type variation, provided a clear phylogeographical picture on a large geographical scale, whereas length variation in one highly polymorphic minisatellite repeat type provided fine-scale phylogeographical information. Mediterranean populations could be divided into four main lineages, a western European lineage, a northern and central Italian lineage, a well-isolated southern Italian (Apulian) lineage, and an eastern European lineage. Variation at the most variable minisatellite repeat type N revealed 19 alleles and allowed the study of seed-mediated gene flow and an estimation of the ratio of pollen to seed flow among neighbouring populations.