Phylogeny, new generic-level classification, and historical biogeography of the Eucera complex (Hymenoptera: Apidae).

The longhorn bee tribe Eucerini (Hymenoptera: Apidae) is a diverse, widely distributed group of solitary bees that includes important pollinators of both wild and agricultural plants. About half of the species in the tribe are currently assigned to the genus Eucera and to a few other related genera. In this large genus complex, comprising ca. 390 species, the boundaries between genera remain ambiguous due to morphological intergradation among taxa. Using ca. 6700 aligned nucleotide sites from six gene fragments, 120 morphological characters, and more than 100 taxa, we present the first comprehensive molecular, morphological, and combined phylogenetic analyses of the 'Eucera complex'. The revised generic classification that we propose is congruent with our phylogeny and maximizes both generic stability and ease of identification. Under this new classification most generic names are synonymized under an expanded genus Eucera. Thus, Tetralonia, Peponapis, Xenoglossa, Cemolobus, and Syntrichalonia are reduced to subgeneric rank within Eucera, and Synhalonia is retained as a subgenus of Eucera. Xenoglossodes is reestablished as a valid subgenus of Eucera while Tetraloniella is synonymized with Tetralonia and Cubitalia with Eucera. In contrast, we suggest that the venusta-group of species, currently placed in the subgenus Synhalonia, should be recognized as a new genus. Our results demonstrate the need to evaluate convergent loss or gain of important diagnostic traits to minimize the use of potentially homoplasious characters when establishing classifications. Lastly, we show that the Eucera complex originated in the Nearctic region in the late Oligocene, and dispersed twice into the Old World. The first dispersal event likely occurred 24.2-16.6 mya at a base of a clade of summer-active bees restricted to warm region of the Old World, and the second 13.9-12.3 mya at the base of a clade of spring-active bees found in cooler regions of the Holarctic. Our results further highlight the role of Beringia as a climate-regulated corridor for bees.

A fossil bee from Early Cretaceous Burmese amber.

The bee fossil record is fragmentary, making it difficult to accurately estimate the antiquity of bee-mediated pollination. Here, we describe a bee fossil [Melittosphex burmensis (new species), Melittosphecidae (new family)] from Early Cretaceous Burmese amber (approximately 100 million years before the present). The fossil provides insights into the morphology of the earliest bees and provides a new minimum date for the antiquity of bees and bee-mediated pollination.

Phylogenetic utility of the major opsin in bees (Hymenoptera: Apoidea): a reassessment.

Major opsin (LW Rh) DNA sequence has been reported to provide useful data for resolving phylogenetic relationships among tribes of corbiculate bees based on analyses of 502 bp of coding sequence. However, the corbiculate tribes are believed to be of Cretaceous age, and strong support for insect clades of this age from small data sets of nucleotide sequence data has rarely been demonstrated. To more critically assess opsin's phylogenetic utility we generated an expanded LW Rh data set by sequencing the same gene fragment from 52 additional bee species from 24 tribes and all six extant bee families. Analyses of this data set failed to provide substantial support for monophyly of corbiculate bees, for relationships among corbiculate tribes, or for most other well-established higher-level relationships among long-tongued bees. However, monophyly of nearly all genera and tribes is strongly supported, indicating that LW Rh provides useful phylogenetic signal at lower taxonomic levels. When our expanded LW Rh data set is combined with a morphological and behavioral data set for corbiculate bees, the results unambiguously support the traditional phylogeny of the corbiculate bee tribes: (Euglossini + (Bombini + (Meliponini + Apini))). This implies a single origin of advanced eusocial behavior among bees rather than dual origins, as proposed by several recent studies.

Australian Lasioglossum + Homalictus form a monophyletic group: resolving the "Australian enigma".

The bee genus Lasioglossum includes > 1,000 species of bees distributed on all continents except Antarctica. Lasioglossum is a major component of the bee fauna in the Holarctic, Ethiopian, and Asian regions and is an important group for investigating the evolution of social behavior in bees. Given its cosmopolitan distribution, the historical biogeography of the genus is of considerable interest. We reconstructed phylogenetic relationships among the subgenera and species within Lasioglossum s.s., using DNA sequence data from a slowly evolving nuclear gene, elongation factor-1 alpha. The entire data set includes > 1,604 aligned nucleotide sites (including three exons plus two introns) for 89 species (17 outgroups plus 72 ingroups). Parsimony and maximum likelihood analyses provide strong evidence that the primarily Indoaustralian subgenera (Homalictus, Chilalictus, Parasphecodes) form a monophyletic group. Bootstrap support for the Australian clade ranged from 73% to 77%, depending on the method of analysis. Monophyly of the Australian Lasioglossum suggests that a single colonization event (by way of Southeast Asia and New Guinea) gave rise to a lineage of > 350 native Indoaustralian bees. We discuss the implications of Australian monophyly for resolving the "Australian enigma"--the similarity in social behavior among the Australian halictine bees relative to that of Holarctic groups.

Phylogeny of the bee genus Halictus (Hymenoptera: halictidae) based on parsimony and likelihood analyses of nuclear EF-1alpha sequence data.

We investigated higher-level phylogenetic relationships within the genus Halictus based on parsimony and maximum likelihood (ML) analysis of elongation factor-1alpha DNA sequence data. Our data set includes 41 OTUs representing 35 species of halictine bees from a diverse sample of outgroup genera and from the three widely recognized subgenera of Halictus (Halictus s.s., Seladonia, and Vestitohalictus). We analyzed 1513 total aligned nucleotide sites spanning three exons and two introns. Equal-weights parsimony analysis of the overall data set yielded 144 equally parsimonious trees. Major conclusions supported in this analysis (and in all subsequent analyses) included the following: (1) Thrincohalictus is the sister group to Halictus s.l., (2) Halictus s.l. is monophyletic, (3) Vestitohalictus renders Seladonia paraphyletic but together Seladonia + Vestitohalictus is monophyletic, (4) Michener's Groups 1 and 3 are monophyletic, and (5) Michener's Group 1 renders Group 2 paraphyletic. In order to resolve basal relationships within Halictus we applied various weighting schemes under parsimony (successive approximations character weighting and implied weights) and employed ML under 17 models of sequence evolution. Weighted parsimony yielded conflicting results but, in general, supported the hypothesis that Seladonia + Vestitohalictus is sister to Michener's Group 3 and renders Halictus s.s. paraphyletic. ML analyses using the GTR model with site-specific rates supported an alternative hypothesis: Seladonia + Vestitohalictus is sister to Halictus s.s. We mapped social behavior onto trees obtained under ML and parsimony in order to reconstruct the likely historical pattern of social evolution. Our results are unambiguous: the ancestral state for the genus Halictus is eusociality. Reversal to solitary behavior has occurred at least four times among the species included in our analysis.

Elongation factor-1 alpha occurs as two copies in bees: implications for phylogenetic analysis of EF-1 alpha sequences in insects.

We report the complete sequence of a paralogous copy of elongation factor-1 alpha (EF-1 alpha) in the honeybee, Apis mellifera (Hymenoptera: Apidae). This copy differs from a previously described copy in the positions of five introns and in 25% of the nucleotide sites in the coding regions. The existence of two paralogous copies of EF-1 alpha in Drosophila and Apis suggests that two copies of EF-1 alpha may be widespread in the holometabolous insect orders. To distinguish between a single, ancient gene duplication and parallel, independent fly and bee gene duplications, we performed a phylogenetic analysis of hexapod EF-1 alpha sequences. Unweighted parsimony analysis of nucleotide sequences suggests an ancient gene duplication event, whereas weighted parsimony analysis of nucleotides and unweighted parsimony analysis of amino acids suggests the contrary: that EF-1 alpha underwent parallel gene duplications in the Diptera and the Hymenoptera. The hypothesis of parallel gene duplication is supported both by congruence among nucleotide and amino acid data sets and by topology-dependent permutation tail probability (T-PTP) tests. The resulting tree topologies are also congruent with current views on the relationships among the holometabolous orders included in this study (Diptera, Hymenoptera, and Lepidoptera). More sequences, from diverse orders of holometabolous insects, will be needed to more accurately assess the historical patterns of gene duplication in EF-1 alpha.

Secondarily solitary: the evolutionary loss of social behavior.

Studies of social behavior frequently assume that evolution proceeds from a solitary state to a social one, and social to social lineages give rise to line are also social, excluding parasitic taxa. Recent phylogenetic studies of some bees contradict this assumption, and more examples are known or hypothesized in other animals. Social behaviour can be lost to give rise to species that are secondarily solitary. Studies of the conditions to the suppression or loss of social behavior can help to illuminate those factors that lead to its origins and maintenance.

DNA fingerprinting data and the problem of non-independence among pairwise comparisons.

Multilocus DNA fingerprinting is commonly used to assess genetic similarity within and between geographically disjunct populations. Typically, the proportion of DNA fingerprinting bands shared between two individuals (SXY) is calculated for all possible pairwise comparisons and the resulting data analysed parametrically to test differences in mean band-sharing among groups. The degree to which covariation among interdependent SXY values (S(ab)-Sbc) biases the analyses is often unknown. Here, we assess the extent of covariation in four DNA fingerprinting studies and evaluate the effectiveness of two corrective procedures, a permutation test and a subsampling routine using only independent pairwise comparisons drawn without replacement from the overall data. Covariation among interdependent SXY values was significantly greater than zero in every data set examined, including those from a bee, a rodent, and two passerine birds. Permutation tests did not correct for interdependence and yielded significance values nearly identical to those derived from uncorrected parametric procedures. In contrast, the subsampling procedure yielded corrected estimates of the standard error that were two to four times larger than those derived parametrically. As a result, comparisons that were significant using parametric tests were either non-significant or only marginally significant with the subsampling routine. We conclude that interdependence among SXY values poses a substantial obstacle to hypothesis testing that must be addressed in future studies.

Mechanisms of spontaneous mutagenesis: an analysis of the spectrum of spontaneous mutation in the Escherichia coli lacI gene.

We have obtained via DNA sequence analysis a spectrum of 174 spontaneous mutations occurring in the lac I gene of Escherichia coli. The spectrum comprised base substitution, frameshift, deletion, duplication and insertion mutations, of which the relative contributions to spontaneous mutation could be estimated. Two thirds of all lacI mutations occurred in the frameshift hotspot site. An analysis of the local DNA sequence suggested that the intensity of this hotspot may depend on structural features of the DNA that extend beyond those permitted by the repeated tetramer at this site. Deletions comprised the largest non-hotspot class (37%). They could be divided into two subclasses, depending on whether they included the lac operator sequence; the latter was found to be a preferred site for deletion endpoints. Most of the deletions internal to the lacI gene were associated with the presence of directly or invertedly repeated sequences capable of accounting for their endpoints. Base substitutions comprised 34% of the non-hotspot events. Unlike the base substitution spectrum obtained via nonsense mutations, G . C----A . T transitions do not predominate. A new base substitution hotspot was discovered at position +6 in the lac operator; its intensity may reflect specific features of the operator DNA. IS1 insertion mutations contributed 12% of the non-hotspot mutations and occurred dispersed throughout the gene in both orientations. Since the lacI gene is not A + T-rich, the contribution of IS1 insertion to spontaneous mutation in general might be underestimated. Single-base frameshift mutations were found only infrequently. In general, they did not occur in runs of a common base. Instead, their occurrence seemed based on the "perfection" of direct or inverted repeats in the local DNA sequence. Three (tandem) duplication events were recovered. No repeated sequences were found that might have determined their endpoints.

Rapid repeated cloning of mutant lac repressor genes.

We have developed a procedure to efficiently recover lac repressor mutations (lacI-) from F'lac onto a single-stranded M13 phage vector. The recovery is based on homologous recombination between F'lac and an M13lac vector. This vector, mRS81, carries the entire Escherichia coli lacI gene as well as the adjacent alpha-complementation region of the lacZ gene, inserted in the AvaI site of the M13 ori region. It also carries a single point mutation in lacZ- alpha which abolishes its alpha-complementing ability. Recovery of lacI- genes from F is based on the conversion of this lacI+Z- alpha phage to lacI-Z+ alpha by recombination with F'lacI-Z+. This double exchange restores its alpha-complementing ability in the absence of any inducer of the lac operon. Detection requires a lacI- alpha-complementation host, which was also constructed in this study. The procedure was developed to obtain rapid nucleotide sequence information on large collections of lacI mutants for the purpose of studying mutational mechanisms and specificities.