Isolation and characterization of microsatellite markers in two species of the neotropical Epicharis (Hymenoptera, Apidae, Centridini) genus and cross-amplification in related species

Genetics can provide invaluable tools for management and conservation of bee populations, which are declining worldwide. Among these tools, microsatellite are very useful molecular markers for population analyses. The aim of this study was to isolate and characterize microsatellites for Epicharis (Anepicharis) dejeanii and Epicharis (Epicharis) nigrita, two Neotropical species of solitary bees, both exhibiting the habit of nesting in aggregations. Microsatellite loci were identified from two enriched genomic libraries. The characterization and analysis of loci were carried out using 35 females of E. dejeanii and 34 of E. nigrita. In total, we report the development of 12 microsatellite loci for E. dejeanii and 13 for E. nigrita. For E. dejeanii, all loci were polymorphic, the number of alleles per locus ranged from 2 to 12, averaging 8.7 and, observed and expected heterozygosity were 0.485 (range 0.229-0.857) and 0.633 (range 0.288-0.843), respectively. For E. nigrita, only nine out of 13 loci amplified were polymorphic, the number of alleles per locus ranged from 2 to 12, averaging 5.5. For this species, the observed and expected heterozygosity were 0.440 (range 0.118-0.676) and 0.545 (range 0.167-0.814), respectively. Cross-amplification of primers was successful in other Centridini species. The two sets of loci described for E. dejeanii and E. nigrita species are polymorphic and informative and show promising applicability for both population genetic approaches and relatedness on these and other Centridini species.

Landscape genomics to the rescue of a tropical bee threatened by habitat loss and climate change

Habitat degradation and climate change are currently threatening wild pollinators, compromising their ability to provide pollination services to wild and cultivated plants. Landscape genomics offers powerful tools to assess the influence of landscape modifications on genetic diversity and functional connectivity, and to identify adaptations to local environmental conditions that could facilitate future bee survival. Here, we assessed range-wide patterns of genetic structure, genetic diversity, gene flow, and local adaptation in the stingless bee Melipona subnitida, a tropical pollinator of key biological and economic importance inhabiting one of the driest and hottest regions of South America. Our results reveal four genetic clusters across the species' full distribution range. All populations were found to be under a mutation-drift equilibrium, and genetic diversity was not influenced by the amount of reminiscent natural habitats. However, genetic relatedness was spatially autocorrelated and isolation by landscape resistance explained range-wide relatedness patterns better than isolation by geographic distance, contradicting earlier findings for stingless bees. Specifically, gene flow was enhanced by increased thermal stability, higher forest cover, lower elevations, and less corrugated terrains. Finally, we detected genomic signatures of adaptation to temperature, precipitation, and forest cover, spatially distributed in latitudinal and altitudinal patterns. Taken together, our findings shed important light on the life history of M. subnitida and highlight the role of regions with large thermal fluctuations, deforested areas, and mountain ranges as dispersal barriers. Conservation actions such as restricting long-distance colony transportation, preserving local adaptations, and improving the connectivity between highlands and lowlands are likely to assure future pollination services.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

Extreme genome diversity in the hyper-prevalent parasitic eukaryote Blastocystis

Blastocystis is the most prevalent eukaryotic microbe colonizing the human gut, infecting approximately 1 billion individuals worldwide. Although Blastocystis has been linked to intestinal disorders, its pathogenicity remains controversial because most carriers are asymptomatic. Here, the genome sequence of Blastocystis subtype (ST) 1 is presented and compared to previously published sequences for ST4 and ST7. Despite a conserved core of genes, there is unexpected diversity between these STs in terms of their genome sizes, guanine-cytosine (GC) content, intron numbers, and gene content. ST1 has 6,544 protein-coding genes, which is several hundred more than reported for ST4 and ST7. The percentage of proteins unique to each ST ranges from 6.2% to 20.5%, greatly exceeding the differences observed within parasite genera. Orthologous proteins also display extreme divergence in amino acid sequence identity between STs (i.e., 59%-61% median identity), on par with observations of the most distantly related species pairs of parasite genera. The STs also display substantial variation in gene family distributions and sizes, especially for protein kinase and protease gene families, which could reflect differences in virulence. It remains to be seen to what extent these inter-ST differences persist at the intra-ST level. A full 26% of genes in ST1 have stop codons that are created on the mRNA level by a novel polyadenylation mechanism found only in Blastocystis. Reconstructions of pathways and organellar systems revealed that ST1 has a relatively complete membrane-trafficking system and a near-complete meiotic toolkit, possibly indicating a sexual cycle. Unlike some intestinal protistan parasites, Blastocystis ST1 has near-complete de novo pyrimidine, purine, and thiamine biosynthesis pathways and is unique amongst studied stramenopiles in being able to metabolize α-glucans rather than β-glucans. It lacks all genes encoding heme-containing cytochrome P450 proteins. Predictions of the mitochondrion-related organelle (MRO) proteome reveal an expanded repertoire of functions, including lipid, cofactor, and vitamin biosynthesis, as well as proteins that may be involved in regulating mitochondrial morphology and MRO/endoplasmic reticulum (ER) interactions. In sharp contrast, genes for peroxisome-associated functions are absent, suggesting Blastocystis STs lack this organelle. Overall, this study provides an important window into the biology of Blastocystis, showcasing significant differences between STs that can guide future experimental investigations into differences in their virulence and clarifying the roles of these organisms in gut health and disease.

Overexpression of GSK3-like Kinase 5 (OsGSK5) in rice (Oryza sativa) enhances salinity tolerance in part via preferential carbon allocation to root starch

Rice (Oryza sativa L.) is very sensitive to soil salinity. To identify endogenous mechanisms that may help rice to better survive salt stress, we studied a rice GSK3-like isoform (OsGSK5), an orthologue of a Medicago GSK3 previously shown to enhance salinity tolerance in Arabidopsis by altering carbohydrate metabolism. We wanted to determine whether OsGSK5 functions similarly in rice. OsGSK5 was cloned and sequence, expression, evolutionary and functional analyses were conducted. OsGSK5 was expressed highest in rice seedling roots and was both salt and sugar starvation inducible in this tissue. A short-term salt-shock (150mM) activated OsGSK5, whereas moderate (50mM) salinity over the same period repressed the transcript. OsGSK5 response to salinity was due to an ionic effect since it was unaffected by polyethylene glycol. We engineered a rice line with 3.5-fold higher OsGSK5 transcript, which better tolerated cultivation on saline soils (EC=8 and 10dSm-2). This line produced more panicles and leaves, and a higher shoot biomass under high salt stress than the control genotypes. Whole-plant 14C-tracing and correlative analysis of OsGSK5 transcript with eco-physiological assessments pointed to the accelerated allocation of carbon to the root and its deposition as starch, as part of the tolerance mechanism.

Beekeeping practices and geographic distance, not land use, drive gene flow across tropical bees

Across the globe, wild bees are threatened by ongoing natural habitat loss, risking the maintenance of plant biodiversity and agricultural production. Despite the ecological and economic importance of wild bees and the fact that several species are now managed for pollination services worldwide, little is known about how land use and beekeeping practices jointly influence gene flow. Using stingless bees as a model system, containing wild and managed species that are presumed to be particularly susceptible to habitat degradation, here we examine the main drivers of tropical bee gene flow. We employ a novel landscape genetic approach to analyse data from 135 populations of 17 stingless bee species distributed across diverse tropical biomes within the Americas. Our work has important methodological implications, as we illustrate how a maximum-likelihood approach can be applied in a meta-analysis framework to account for multiple factors, and weight estimates by sample size. In contrast to previously held beliefs, gene flow was not related to body size or deforestation, and isolation by geographic distance (IBD) was significantly affected by management, with managed species exhibiting a weaker IBD than wild ones. Our study thus reveals the critical importance of beekeeping practices in shaping the patterns of genetic differentiation across bee species. Additionally, our results show that many stingless bee species maintain high gene flow across heterogeneous landscapes. We suggest that future efforts to preserve wild tropical bees should focus on regulating beekeeping practices to maintain natural gene flow and enhancing pollinator-friendly habitats, prioritizing species showing a limited dispersal ability.

From dusk till dawn: the Arabidopsis thaliana sugar starving responsive network

Plant growth and development are tightly controlled by photosynthetic carbon availability. The understanding of mechanisms governing carbon partitioning in plants will be a valuable tool in order to satisfy the rising global demand for food and biofuel. The goal of this study was to determine if sugar starvation responses were transcriptionally coordinated in Arabidopsis thaliana. A set of sugar-starvation responsive (SSR) genes was selected to perform a co-expression network analysis. Posteriorly, a guided-gene approach was used to identify the SSR-network from public data and to discover candidate regulators of this network. In order to validate the SSR network, a global transcriptome analysis was realized on three A. thaliana starch-deficient mutants. The starch-deficient phenotype in leaves induces sugar starvation syndrome at the end of the night due to the absence of photosynthesis. Promoter sequences of genes belonging to the SSR-network were analyzed in silico reveling over-represented motifs implicated in light, abscisic acid, and sugar responses. A small cluster of protein encoding genes belonging to different metabolic pathways, including three regulatory proteins, a protein kinase, a transcription factor, and a blue light receptor, were identified as the cornerstones of the SSR co-expression network. In summary, a large transcriptionally coordinated SSR network was identified and was validated with transcriptional data from three starch-deficient mutant lines. Candidate master regulators of this network were point out.

Cytochrome c oxidase I primers for corbiculate bees: DNA barcode and mini-barcode

Bees (Apidae), of which there are more than 19 900 species, are extremely important for ecosystem services and economic purposes, so taxon identity is a major concern. The goal of this study was to optimize the DNA barcode technique based on the Cytochrome c oxidase (COI) mitochondrial gene region. This approach has previously been shown to be useful in resolving taxonomic inconsistencies and for species identification when morphological data are poor. Specifically, we designed and tested new primers and standardized PCR conditions to amplify the barcode region for bees, focusing on the corbiculate Apids. In addition, primers were designed to amplify small COI amplicons and tested with pinned specimens. Short barcode sequences were easily obtained for some Bombus century-old museum specimens and shown to be useful as mini-barcodes. The new primers and PCR conditions established in this study proved to be successful for the amplification of the barcode region for all species tested, regardless of the conditions of tissue preservation. We saw no evidence of Wolbachia or numts amplification by these primers, and so we suggest that these new primers are of broad value for corbiculate bee identification through DNA barcode.

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2011-31 January 2012

This article documents the addition of 473 microsatellite marker loci and 71 pairs of single-nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Barteria fistulosa, Bombus morio, Galaxias platei, Hematodinium perezi, Macrocentrus cingulum Brischke (a.k.a. M. abdominalis Fab., M. grandii Goidanich or M. gifuensis Ashmead), Micropogonias furnieri, Nerita melanotragus, Nilaparvata lugens Stål, Sciaenops ocellatus, Scomber scombrus, Spodoptera frugiperda and Turdus lherminieri. These loci were cross-tested on the following species: Barteria dewevrei, Barteria nigritana, Barteria solida, Cynoscion acoupa, Cynoscion jamaicensis, Cynoscion leiarchus, Cynoscion nebulosus, Cynoscion striatus, Cynoscion virescens, Macrodon ancylodon, Menticirrhus americanus, Nilaparvata muiri and Umbrina canosai. This article also documents the addition of 116 sequencing primer pairs for Dicentrarchus labrax.

Mitochondrial DNA diversity of orchid bee Euglossa fimbriata (Hymenoptera: Apidae) populations assessed by PCR-RFLP

Euglossa fimbriata is a euglossine species widely distributed in Brazil and occurring primarily in Atlantic Forest remnants. In this study, the genetic mitochondrial structure of E. fimbriata from six Atlantic Forest fragments was studied by RFLP analysis of three PCR-amplified mtDNA gene segments (16S, COI-COII, and cyt b). Ten composite haplotypes were identified, six of which were exclusive and represented singleton mitotypes. Low haplotype diversity (0.085-0.289) and nucleotide diversity (0.000-0.002) were detected within samples. AMOVA partitioned 91.13% of the overall genetic variation within samples and 8.87% (phi(st) = 0.089; P < 0.05) among samples. Pairwise comparisons indicated high levels of differentiation among some pairs of samples (phi(st) = 0.161-0.218; P < 0.05). These high levels indicate that these populations of E. fimbriata, despite their highly fragmented landscape, apparently have not suffered loss of genetic variation, suggesting that this particular population is not currently endangered.

Morphometric and genetic changes in a population of Apis mellifera after 34 years of Africanization

Though the replacement of European bees by Africanized honey bees in tropical America has attracted considerable attention, little is known about the temporal changes in morphological and genetic characteristics in these bee populations. We examined the changes in the morphometric and genetic profiles of an Africanized honey bee population collected near where the original African swarms escaped, after 34 years of Africanization. Workers from colonies sampled in 1968 and in 2002 were morphometrically analyzed using relative warps analysis and an Automatic Bee Identification System (ABIS). All the colonies had their mitochondrial DNA identified. The subspecies that mixed to form the Africanized honey bees were used as a comparison for the morphometric analysis. The two morphometric approaches showed great similarity of Africanized bees with the African subspecies, Apis mellifera scutellata, corroborating with other markers. We also found the population of 1968 to have the pattern of wing venation to be more similar to A. m. scutellata than the current population. The mitochondrial DNA of European origin, which was very common in the 1968 population, was not found in the current population, indicating selective pressure replacing the European with the African genome in this tropical region. Both morphometric methodologies were very effective in discriminating the A. mellifera groups; the non-linear analysis of ABIS was the most successful in identifying the bees, with more than 94% correct classifications.

Mitochondrial DNA polymorphism among populations of Melipona quadrifasciata quadrifasciata Lepeletier (Apidae: Meliponini) from southern Brazil

The geographical distribution of the Brazilian endemic stingless bee Melipona quadrifasciata quadrifasciata Lepeletier ranges from Rio Grande do Sul to Minas Gerais states. The objective of the present study was to verify mtDNA polymorphisms among samples of M. q. quadrifasciata collected in southern Brazil. Twenty nine colonies from three localities (Blumenau and Mafra/SC and Prudentópolis/ PR) were sampled. Seven mtDNA regions were amplified and further digested with 15 restriction enzymes (PCR-RFLP). Five composite haplotypes were identified, with two unique to samples from Prudentópolis and the remaining three to samples from Mafra and/or Blumenau.

Cytochrome-b variation in Apis mellifera samples and its association with COI-COII patterns

Five Mbo I (Mbo-A, Mbo-M, Mbo-C(1), Mbo-C(2) and Mbo-C(3)) and Hinf I (Hinf-1 to Hinf-5) patterns were observed in Apis mellifera samples after restriction of a 485 bp fragment of the mitochondrial cytochrome-b (cyt-b) gene. Associating the cyt-b Restriction fragment length polymorphism (RFLP) pattern of each sample to its respective previously established COI-COII (Dra I sites) pattern, five restriction patterns (Mbo-C(1), Mbo-C(2), Mbo-C(3), Hinf-1 and Hinf-4) were observed in samples of maternal origin associated to the evolutionary branch C. No deletions or insertions were observed and the nucleotide substitution rate was estimated at 5.4%. Higher nucleotide diversity was observed among the branch C-haplotypes when compared with A and M lineages. Further studies are needed to confirm if the cyt-b + COI-COII haplotypes help to assign certain phylogeographic patterns to the branch C and to clarify phylogenetic relationships among A. mellifera subspecies.

Genetic structure of Africanized honeybee populations (Apis mellifera L.) from Brazil and Uruguay viewed through mitochondrial DNA COI-COII patterns

Mitochondrial genotypes of Africanized honeybees from Brazil and Uruguay were surveyed by DraI restriction of the COI-COII region. Eleven mitotypes were found, three of which had not previously been described (A28-A30). Out of 775 samples (725 from Brazil, 50 from Uruguay), 197 were A1 and 520 were A4. A1 frequency increases toward the north of Brazil, whereas A4 frequency increases toward the south, a pattern echoing the African distribution. The origin of the A4 and most of the A1 African patterns can be attributed to the introduction of Apis mellifera scutellata into Brazil in 1956. The A29 and A30 patterns have the P1 sequence observed in many Iberian Peninsula samples, which represent the traces of the introductions into Brazil and Uruguay by settlers.

Alelle number and heterozigosity for microsatellite loci in different stingless bee species (Hymenoptera: Apidae, Meliponini)

In the present study we compare genetic characteristics (allele diversity and observed heterozygosity) of microsatellite loci, from three stingless bee species (Plebeia remota Holmberg, Partamona mulata Moure In Camargo and Partamona helleri Friese), amplified by using heterospecific primers originally designed for Melipona bicolor Lepeletier and Scaptotrigona postica Latreille. We analyzed 360 individuals of P. remota from 72 nests, 58 individuals of R. mulata from 58 nests, and 47 individuals of P. helleri from 47 nests. The three species studied showed low level of polymorphism for the loci amplified with primers derived from M. bicolor. However, for the loci amplified with primers derived from S. postica, only P. remota presented low level of polymorphism.

EIMERIA SPECIES FROM THE EUROPEAN POND TURTLE, EMYS ORBICULARIS (REPTILIA: TESTUDINES), IN GALICIA (NW SPAIN), WITH DESCRIPTION OF TWO NEW SPECIES

Parasitological examination of feces from 44 Emys orbicularis from Galicia (NW Spain) revealed the presence of 2 new eimerian species, Eimeria gallaeciaensis sp. n. and E. emydis sp. n., as well as E. mitraria (Laveran and Mesnil, 1902) Doflein, 1909. Oocysts of E. gallaeciaensis n. sp. were found in 20 of 44 (45.4%) turtles and are subspherical to lightly ovoid-ellipsoid, 19.3 x 16.0 (17-22 x 15-18) microm, shape index 1.2 (1.1-1.3), with a smooth, single-layered wall. Micropyle and polar granule are absent, but an oocyst residuum is present. Sporocysts are ellipsoid, 9.7 x 5.1 (9-10 x 5-6) microm, shape index 1.9 (1.7-2.0), each with a sporocyst residuum and a conical Stieda body usually bearing 1-4 short and thin projections. Oocysts of E. emydis n. sp. were found in the feces of 5 of 44 (11.4%) turtles and are ovoid, rarely pear-shaped, 22.6 x 17.0 (20-25 x 15.5-18) microm, shape index 1.3 (1.2-1.5), with a smooth, single-layered wall with a slight thinning at the pointed end. Micropyle and polar granule are absent, and an oocyst residuum is present. Sporocysts are ellipsoid, 11.4 x 6.0 (9-13 x 5-7) microm, shape index 1.9 (1.6-2.2), each with sporocyst residuum and a prominent Stieda body bearing 3-5 club-shaped projections. In addition to the new species described, this is the first report of E. mitraria parasitizing E. orbicularis.

Phylogenetic relationships of honey bees (Hymenoptera:Apinae:Apini) inferred from nuclear and mitochondrial DNA sequence data

Two different genomic regions (ND2 mitochondrial gene and EF1-alpha intron) were PCR amplified, cloned and sequenced for the ten known honey bee species collected within their natural range distribution. DNA sequences were analyzed using parsimony, distance and maximum likelihood methods to investigate phylogenetic relationships within Apis. The phylogenetic analyses strongly supported the basic topology recoverable from morphometric analysis, grouping the honey bees into three major clusters: giant bees (A. dorsata, A. binghami, and A. laboriosa), dwarf bees (A. andreniformis and A. florea), and cavity-nesting bees (A. mellifera, A. cerana, A. koschevnikovi, A. nuluensis, and A. nigrocincta). However, the clade of Asian cavity-nesting bees included paraphyletic taxa. Exemplars of Apis cerana collected from divergent portions of its range were less related to each other than were sympatric A. cerana, A. nuluensis, and A. nigrocincta taxa. Nucleotide sequence divergence between allopatrically distributed western (A. mellifera) and eastern (A. cerana, A. koschevnikovi, A. nigrocincta, and A. nuluensis) cavity-nesting species, around 18% for the mitochondrial gene and 10-15% for the nuclear intron, suggested an earlier divergence for these groups than previously estimated from morphometric and behavioral studies. This latter finding neccessitates reevaluation of the hypothesized origin of extant European, African, and west Asian Apis mellifera. Sequence divergence between A. laboriosa and A. dorsata was consistent with behavioral data and supports the species status of A. laboriosa.

A cytoplasmically inherited mutant controlling early chloroplast development in barley seedlings

Cytoplasmic line 2 (CL2) has been previously reported as a cytoplasmically inherited chlorophyll-deficient mutant selected from a chloroplast-mutator genotype of barley. It was characterized by a localized effect on the upper part of the first-leaf blade. At emergence the CL2 seedlings-phenotype varied from a grainy light green to an albino color. They gradually greened during the following days, starting from the base of the blade and extending to cover most of its surface when it was fully grown. The present results, from both light microscopy and transmission electron microscopy (TEM), confirmed the previously described positional and time-dependent expression of the CL2 syndrome along the first-leaf blade. During the first days after emergence, light microscopy showed a normally developed chloroplast at the middle part of the CL2 first-leaf blade, meanwhile at the tip only small plastids were observed. TEM showed that the shapes and the internal structure of the small plastids were abnormal, presenting features of proplastids, amyloplasts and/or senescent gerontoplasts. Besides, they lack plastid ribosomes, contrasting with what was observed inside chloroplasts from normal tips, which presented abundant ribosomes. Phenotypic observations and spectrophotometric analysis of seedlings produced by mother plants that had been grown under different temperatures indicated that higher temperatures during seed formation were negatively associated with pigment content in CL2 seedlings. In contrast, higher temperatures during the growth of CL2 seedlings have been associated with increased pigment content. Aqueous solution with kanamycin and streptomycin, which are antibiotics known to interfere with plastid gene translation, were used for imbibition of wild-type and CL2 seeds. Antibiotic treatments differentially reduced the chlorophyll content in the upper part of the first-leaf blade in CL2, but not in wild-type seedlings. These results suggest that in the wild-type, plastid-gene proteins which are necessary for chloroplast development and chlorophyll synthesis in the upper part of the first-leaf blade are usually synthesized during embryogenesis. However, under certain circumstances, in CL2 seedlings, they would be synthesized after germination. In addition, a shortening of the sheath has been observed in association with pigment decrease suggesting the existence of plastid factors affecting the expression of some nuclear genes. We consider the CL2 mutant a unique experimental material useful to study biological phenomena and external factors regulating plastid, and nuclear gene expression during embryogenesis and early seedling development.

[Scintigraphic pattern in a case of Tarlov cyst]

We report the case of a 40-year-old man remitted to our department with a history of lower back pain and sciatica with no history of trauma. The laboratory analyses showed normal values whereas plain radiographs showed a sacrum rarefaction area. A 99mTc-MDP bone scintigraphy was performed to evaluate the lumbosacral area. Planar images did not show any abnormality. SPECT images revealed photopenic abnormality in the second sacral vertebral right hemibody, with no peripherally increased radiotracer accumulation. Subsequent MRI and CT myelography demonstrated the nature of the photopenic area as secondary to vertebral erosion by sacral perineurial cyst (Tarlov cyst).

Phylogeographical autocorrelation of phenotypic evolution in honey bees (Apis mellifera L.)

The analysis of phenotypic divergence among local populations within a species has been traditionally performed in a spatial context, although advances in genetic analysis using mtDNA have permitted a simultaneous evaluation of geographical and historical patterns of variation, so-called phylogeographical analysis. In this paper, we combine these two dimensions of variation (geographical space and phylogenetic history) to evaluate patterns of phenotypic evolution in honey bees (Apis mellifera L.). Data on 39 phenotypic traits, derived from 417 colonies grouped into 14 subspecies, were analysed using autocorrelation methods. Mantel tests indicated that the relationship between phenotypic divergence, estimated by Euclidean distances among subspecies' morphological centroids, was significant both when compared to geographical distance (r=0.371; P < 0.01) and to genetic distance (estimated as sequence divergence (%) in a mtDNA region encompassing part of the NADH dehydrogenase subunit 2 and isoleucine transfer RNA (r=0.329; P < 0.01)). For the analysis of each trait, the effects of the geographical co-ordinates (latitude and longitude of subspecies geographical range) and of the phylogenetic patterns (defined by eigenvectors of the genetic distance matrix) on phenotypic variation were simultaneously analysed using an extension of a recently developed model, called Phylogenetic Eigenvector Regression (PVR). In general terms, the partial regression slopes indicated that the variation in the characters traditionally associated with adaptive processes, such as body and wing size, were better explained by geographical position. However, characters usually thought to be neutral, such as wing venation angle, were more associated with phylogeny. This is expected because PVR can be interpreted as a partition model, in which adaptive variation tends to be independent of phylogeny (and, in this case, associated with geography). In addition, the first principal component derived from the expected values of the model for each trait, which can be interpreted as the phenotypic variation predicted by phylogeny, is more structured in a north-south cline than are the original data, supporting an adaptive interpretation. The phylogeographical autocorrelation analyses performed in this study show that different traits are more related to one of the two dimensions of variation (geography and phylogeny), and these patterns can furnish insights into the nature of phenotypic evolution in these organisms.

Molecular phylogenetics of honey bee subspecies (Apis mellifera L.) inferred from mitochondrial DNA sequence

A mitochondrial DNA region encompassing part of the NADH dehydrogenase subunit 2 and isoleucine transfer RNA genes was PCR amplified, cloned, and sequenced for 14 morphometrically identified Apis mellifera subspecies and the New World "Africanized" honeybee. Twenty different haplotypes were detected and phylogenetic analyses supported the existence of 3 or 4 major subspecies groups similar to those based on morphometric measurements. However, some discrepancies are reported concerning the subspecies composition of each group. Based on the sequence divergence of Drosophila (2% per Myr) we found that the four lineages may have diverged around 0.67 Myr. The variability found in this region enables us to infer phylogenetic relationships and test hypotheses concerning subspecies origin, dispersion, and biogeography.