A molecular phylogeny of the Old World stingless bees (Hymenoptera: Apidae: Meliponini) and the non-monophyly of the large genus Trigona

Abstract.  We examined the inter- and infrageneric relationships of Old World Meliponini with a near-complete sampling of supra-specific taxa. DNA sequences for the taxa were collected from four genes (mitochondrial 16S rRNA, nuclear long-wavelength rhodopsin copy 1 (opsin), elongation factor-1α copy F2 and arginine kinase). Additional sampling of New World taxa indicated that Trigona sensu lato is not monophyletic: Trigona from the Indo-Malayan/Australasian Regions forms a large clade distantly related to the Neotropical Trigona . A separate clade comprises the Afrotropical meliponines, and includes the 'minute' species found in the Afrotropical, Indo-Malayan and Australasian Regions. The Neotropical genus Melipona , by contrast with previous investigations, is not the sister lineage to the remaining stingless bees, but falls within the strongly supported Neotropical clade. These results constitute the framework for a revised classification and ongoing biological investigations of Meliponini. A single taxonomic change, Heterotrigona bakeri stat.n. , is proposed on the basis of sequence divergence.     

Task-partitioned nectar transfer in stingless bees: work organisation in a phylogenetic context

Abstract 1. The eusocial corbiculate bee tribes comprise the Apini (honey bees), Bombini (bumble bees), and Meliponini (stingless bees). Honey bee foragers ( Apis ) transfer nectar to receiver bees within the nest. This is an example of task partitioning, in which a task is split into sub-tasks connected by material transfer. Nectar transfer does not occur in Bombini. Although it is reported in some species of Meliponini, it has not been subject to detailed study.
2. Nectar transfer was investigated in five genera of Meliponini from Yucatan, Mexico ( Melipona , Trigona , Scaptotrigona , Nannotrigona , and Plebeia ). Nectar transfer occurred in all species and for > 99% of foragers. Multiple transfer, in which a forager unloads nectar to more than one receiver, occurred but at a lower level than in Apis . In M. beecheii , multiple transfer was associated strongly with putative recruitment dances.
3. The data provide some support for the hypothesis that task partitioning is favoured by large colony size, in that the Meliponini never have small colonies because colonies are swarm founded. This ensures that colonies are always large enough to prevent delays in finding a transfer partner imposing high costs. Further tests of this hypothesis are suggested.
4. Viewed in a phylogenetic context, the most parsimonious interpretation is that nectar transfer evolved once in the clade (Apini + Meliponini).     

Phylogenetics and biogeography of the parasitic genus Thesium L. (Santalaceae), with an emphasis on the Cape of South Africa

Thesium is a large genus of parasitic shrubs belonging to tribe Thesieae of Santalaceae. It has a principally Old World distribution, with the greatest diversity being found in southern Africa. Little is known about the relationships within Thesium or its relationships with its closest relatives. In this article, we present a first estimate of species‐level phylogenetic relationships in Thesium based on internal transcribed spacer (ITS) and trnL–trnF sequence data, and use this to explore the biogeographical history of the group. One hundred and four samples representing 72 Thesium spp. were included in a phylogenetic analysis. Plastid and combined data resolve Thesium as paraphyletic relative to Thesidium and Austroamericium with high posterior probability and bootstrap support. ITS sequence data place Thesidium as sister to a large Thesium clade, but with weak support. Ancestral range reconstruction and dating analysis suggest a southern African origin for the group, with a crown age of 39.1 ± 11.9 Mya, followed by dispersal into Europe and South America. A large clade of Cape species split in the Miocene from a clade comprising tropical species (25.5 ± 7.3 Mya) with the diversification of extant species beginning at 16.7 ± 6.3 Mya. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 162 , 435–452.     

Phylogeny,diversification patterns and historical biogeography of euglossine orchid bees (Hymenoptera: Apidae)

The orchid bees constitute a clade of prominent insect pollinators distributed throughout the Neotropical region. Males of all species collect fragrances from natural sources, including flowers, decaying vegetation and fungi, and store them in specialized leg pockets to later expose during courtship display. In addition, orchid bees provide pollination services to a diverse array of Neotropical angiosperms when foraging for food and nesting materials. However, despite their ecological importance, little is known about the evolutionary history of orchid bees. Here, we present a comprehensive molecular phylogenetic analysis based on ～4.0 kb of DNA from four loci [cytochrome oxidase (CO1), elongation factor 1‐α (EF1‐α), arginine kinase (ArgK) and RNA polymerase II (Pol‐II)] across the entire tribe Euglossini, including all five genera, eight subgenera and 126 of the approximately 200 known species. We investigated lineage diversification using fossil‐calibrated molecular clocks and the evolution of morphological traits using disparity‐through‐time plots. In addition, we inferred past biogeographical events by implementing model‐based likelihood methods. Our dataset supports a new view on generic relationships and indicates that the cleptoparasitic genus Exaerete is sister to the remaining orchid bee genera. Our divergence time estimates indicate that extant orchid bee lineages shared a most recent common ancestor at 27–42 Mya. In addition, our analysis of morphology shows that tongue length and body size experienced rapid disparity bursts that coincide with the origin of diverse genera (Euglossa and Eufriesea). Finally, our analysis of historical biogeography indicates that early diversification episodes shared a history on both sides of Mesoamerica, where orchid bees dispersed across the Caribbean, and through a Panamanian connection, thus reinforcing the hypothesis that recent geological events (e.g. the formation of the isthmus of Panama) contributed to the diversification of the rich Neotropical biota. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 552–572.     

Phylogeography of the whelk genus Cominella (Gastropoda: Buccinidae) suggests long‐distance counter‐current dispersal of a direct developer

The gastropod genus Cominella Gray, 1850 consists of approximately 20 species that inhabit a wide range of marine environments in New Zealand and Australia, including its external territory, the geographically isolated Norfolk Island. This distribution is puzzling, however, with apparently closely‐related species occurring either side of the Tasman Sea, even though all species are considered to have limited dispersal abilities. To determine how Cominella attained its current distribution, we derived a dated molecular phylogeny, which revealed a clade comprising all the Australian and Norfolk Island species nested within four clades of solely New Zealand species. This Australian clade diverged well after the vicariant separation of New Zealand from Australia, and implies two long‐distance dispersal events: a counter‐current movement across the Tasman Sea from New Zealand to Australia, occurring at the origination of the clade, followed by the colonization of Norfolk Island. The biology of Cominella suggests that the most likely method of long‐distance dispersal is rafting as egg capsules. Our robust phylogeny also means that the current Cominella classification requires revision. We propose that our clades be recognized as subgenera: Cominella (s.s.), Cominista, Josepha, Cominula, and Eucominia, with each subgenus comprising only of New Zealand or Australian species. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 315–332.     

Phylogenetics and taxonomy of the New World leafy spurges,Euphorbia section Tithymalus (Euphorbiaceae)

The 480 species of leafy spurges, Euphorbia subgenus Esula, represent the main temperate radiation in the large genus Euphorbia. This group is distributed primarily in temperate Eurasia, but with smaller, disjunct centres of diversity in the mountains of the Old World tropics, in temperate southern Africa and in the New World. The majority of New World diversity (32 species) occurs in a single section, section Tithymalus. We analysed sequences of the nrITS and plastid ndhF, trnH‐psbA, trnS‐trnG and trnD‐trnT regions to reconstruct the phylogeny of section Tithymalus and to examine the origins and diversification of the species native to the New World. Our results indicate that the New World species of section Tithymalus form a clade that is sister to the widespread, weedy E. peplus. The New World species fall into two primary groups: a ‘northern annual clade’ from eastern North America and a diverse clade of both annual and perennial species that is divided into three subgroups. Within the second group, there is a small ‘southern annual clade’ from Texas and northern Mexico, a perennial ‘Brachycera clade’ from the western United States and northern Mexico, and a perennial ‘Esuliformis clade’ from montane areas of Mexico, Guatemala, Honduras and the Caribbean island of Hispaniola. Ancestral state reconstructions indicate that the annual habit probably evolved in the ancestor of E. peplus and the New World clade, with a subsequent reversal to the perennial habit. In conjunction with this phylogenetic framework, the New World species of section Tithymalus are comprehensively reviewed. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 175 , 191–228.     

A comprehensive phylogeny of the bumble bees (Bombus)

Bumble bees ( Bombus Latreille) occupy a wide diversity of habitats, from alpine meadows to lowland tropical forest, yet they appear to be similar in morphology throughout their range, suggesting that behavioural adaptations play a more important role in colonizing diverse habitats. Notwithstanding their structural homogeneity, bumble bees exhibit striking inter- and intraspecific variation in colour pattern, purportedly the outcome of mimetic evolution. A robust phylogeny of Bombus would provide the framework for elucidating the history of their wide biogeographical distribution and the evolution of behavioural and morphological adaptations, including colour pattern. However, morphological studies of bumble bees have discovered too few phylogenetically informative characters to reconstruct a robust phylogeny. Using DNA sequence data, we report the first nearly complete species phylogeny of bumble bees, including most of the 250 known species from the 38 currently recognized subgenera. Bayesian analysis of nuclear (opsin, EF-1α, arginine kinase, PEPCK) and mitochondrial (16S) sequences results in a highly resolved and strongly supported phylogeny from base to tips, with clear-cut support for monophyly of most of the conventional morphology-based subgenera. Most subgenera fall into two distinct clades ( short-faced and long-faced ) associated broadly with differences in head morphology. Within the short-faced clade is a diverse New World clade, which includes nearly one-quarter of the currently recognized subgenera, many of which are restricted to higher elevations of Central and South America. The comprehensive phylogeny provides a firm foundation for reclassification and for evaluating character evolution in the bumble bees.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 161–188.     

Erratum to: Phylogeny and biogeography of the genus Stevia (Asteraceae: Eupatorieae): an example of diversification in the Asteraceae in the new world

The genus Stevia comprises approximately 200 species, which are distributed in North and South America, and are representative of the species diversity of the Asteraceae in the New World. We reconstructed the phylogenetic relationships using sequences of ITS and cpDNA and estimated the divergence times of the major clade of this genus. Our results suggested that Stevia originated in Mexico 7.0–7.3 million years ago (Mya). Two large clades, one with shrub species and another with herb species, were separated at about 6.6 Mya. The phylogenetic reconstruction suggested that an ancestor of Stevia was a small shrub in temperate pine–oak forests and the evolutionary change from a shrub state to a herb state occurred only once. A Brazilian clade was nested in a Mexican herb clade, and its origin was estimated to be 5.2 Mya, suggesting that the migration from North America to South America occurred after the formation of the Isthmus of Panama. The species diversity in Mexico appears to reflect the habitat diversity within the temperate pine–oak forest zone. The presence of many conspecific diploid–polyploid clades in the phylogenetic tree reflects the high frequency of polyploidization among the perennial Stevia species.

     

Phylogenetic analysis of the corbiculate Apinae based on morphology of the sting apparatus (Hymenoptera: Apidae)

This study aimed to test the various competing hypotheses regarding the relationships among the four tribes of corbiculate apine bees (Euglossini “orchid bees”, Bombini “bumble bees”, Meliponini “stingless bees”, and Apini “honey bees”) with a completely new set of previously unstudied morphological characters derived from the sting apparatus. The result was one most parsimonious tree of 49 steps, CI = 89, RI = 93 that is perfectly congruent with most studies based on morphological and combined morphological/molecular data, i.e., Euglossini + (Bombini + (Meliponini + Apini)), supporting a well accepted scenario of social evolution for these bees. This data matrix was then combined with other published matrices for this group in order to perform simultaneous analyses. The problem of how to best combine the multiple matrices that did not use the same exemplars was investigated. © The Willi Hennig Society 2007.     

Contrasting global‐scale evolutionary radiations: phylogeny,diversification, and morphological evolution in the major clades of iguanian lizards

Parallel evolutionary radiations in adjacent locations have been documented in many systems, but typically at limited geographical scales. Here, we compare patterns of evolutionary radiation at the global scale in iguanian lizards, the dominant clade of lizards. We generated a new time‐calibrated phylogeny including 153 iguanian species (based on mitochondrial and nuclear data) and obtained data on morphology and microhabitats. We then compared patterns of species diversification, morphological disparity, and ecomorphological relationships in the predominantly Old World and New World clades (Acrodonta and Pleurodonta, respectively), focusing on the early portions of these radiations. Acrodonts show relatively constant rates of species diversification and disparity over time. In contrast, pleurodonts show an early burst of species diversification and less‐than‐expected morphological disparity early in their history, and slowing diversification and increasing disparity more recently. Analyses including all species (with MEDUSA) suggest accelerated diversification rates in certain clades within both Acrodonta and Pleurodonta, which strongly influences present‐day diversity patterns. We also find substantial differences in ecomorphological relationships between these clades. Our results demonstrate that sister clades in different global regions can undergo very different patterns of evolutionary radiation over similar time frames. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●● , ●●–●●.     

Phylogenetics and morphological evolution of coral‐dwelling barnacles (Balanomorpha: Pyrgomatidae)

Pyrgomatid barnacles are a family of balanomorphs uniquely adapted to symbiosis on corals. The evolution of the coral‐dwelling barnacles is explored using a multi‐gene phylogeny (COI, 16S, 12S, 18S, and H3) and phenotypic trait‐mapping. We found that the hydrocoral associate Wanella should be excluded, while some archaeobalanids in the genus Armatobalanus should be included in the Pyrgomatidae. Three well supported clades were recovered: clade I is the largest group and is exclusively Indo‐West Pacific, clade II contains two plesiomorphic Indo‐West Pacific genera, while clade III is comprised of East and West Atlantic taxa. Some genera did not form reciprocally monophyletic groups, while the genus Trevathana was found to be paraphyletic and to include members of three other apomorphic genera/tribes. The highly unusual coral‐parasitic hoekiines appear to be of recent origin and rapidly evolving from Trevathana sensu lato. Pyrgomatids include six‐, four‐, and one‐plated forms, and exhibit convergent evolutionary tendencies towards skeletal reduction and fusion, loss of cirral armature, and increased host specificity. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 162–179.     

‘After Africa’: the evolutionary history and systematics of the genus Charaxes Ochsenheimer (Lepidoptera: Nymphalidae) in the Indo‐Pacific region

The predominantly Afrotropical genus Charaxes is represented by 31 known species outside of Africa (excluding subgenus Polyura Billberg). We explored the biogeographic history of the genus using every known non‐African species, with several African species as outgroup taxa. A phylogenetic hypothesis is proposed, based on molecular characters of the mitochondrial genes cytochrome oxidase subunit I (COI) and NADH dehydrogenase 5 (ND5), and the nuclear wingless gene. Phylogenetic analyses based on maximum parsimony and Bayesian inference of the combined dataset implies that the Indo‐Pacific Charaxes form a monophyletic assemblage, with the exception of Charaxes solon Fabricius. Eight major lineages are recognized in the Indo‐Pacific, here designated the solon (+African), elwesi, harmodius, amycus, mars, eurialus, latona, nitebis, and bernardus clades. Species group relationships are concordant with morphology and, based on the phylogeny, we present the first systematic appraisal and classification of all non‐African species. A biogeographical analysis reveals that, after the genus originated in Africa, the evolutionary history of Charaxes in the Indo‐Pacific, in particular Wallacea, may be correlated with the inferred geological and climatic history of the region. We propose that Wallacea was the area of origin of all Charaxes (excluding C. solon) occurring to the east of Wallace's [1863] Line. The earliest Indo‐Pacific lineages appear to have diverged subsequent to the initial fragmentation of a palaeo‐continent approximately 13 million years ago. Further diversification in Indo‐Pacific Charaxes appears primarily related to climatic changes during the Pliocene and possibly as recently as the Pleistocene. Although both dispersal and vicariance have played important roles in the evolution of the genus within the region, the latter has been particularly responsible for diversification of Charaxes in Wallacea. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 457–481.     

Biogeography of the grasses (Poaceae): a phylogenetic approach to reveal evolutionary history in geographical space and geological time

The grasses (Poaceae) are the fifth most diverse family of angiosperms, including 800 genera and more than 10 000 species. Few phylogenetic studies have tried to investigate palaeo‐biogeographical and palaeo‐ecological scenarios that may have led to present‐day distribution and diversity of grasses at the family level. We produced a dated phylogenetic tree based on combined plastid DNA sequences and a comprehensive sample of Poaceae. Furthermore, we produced an additional tree using a supermatrix of morphological and molecular data that included all 800 grass genera so that ancestral biogeography and ecological habitats could be inferred. We used a likelihood‐based method, which allows the estimation of ancestral polymorphism in both biogeographical and ecological analyses for large data sets. The origin of Poaceae was retrieved as African and shade adapted. The crown node of the BEP + PACCMAD clade was dated at 57 Mya, in the early Eocene. Grasses dispersed to all continents by approximately 60 million years after their Gondwanan origin in the late Cretaceous. PACCMAD taxa adapted to open habitats as early as the late Eocene, a date consistent with recent phytolith fossil data for North America. C₄ photosynthesis first originated in Africa, at least for Chloridoideae in the Eocene at c. 30 Mya. The BEP clade members adapted to open habitats later than PACCMAD members; this was inferred to occur in Eurasia in the Oligocene. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 162 , 543–557.     

Phylogenetic relationships of the pygmy rice rats of the genus Oligoryzomys Bangs, 1900 (Rodentia: Sigmodontinae)

Sequences from two mitochondrial genes (cytochrome b and NADH1) were used to produce a molecular phylogeny for 12 named and two undescribed species of the genus Oligoryzomys. All analyses placed Oligoryzomys microtis as the most basal taxon, a finding consistent with previous studies that suggested the west‐central Amazon as a centre of origin for the tribe Oryzomyini to which Oligoryzomys belongs. Biogeographically, this suggests that Oligoryzomys had a South American origin, and later advanced northwards, entering Central America and Mexico more recently. Different analyses have provided consistent support for several additional clades that did not necessarily agree with the species groups hypothesized by previous studies. A molecular clock derived for these data suggests an origin for the genus of 6.67 Mya, with most speciation within the genus occurring between 3.7 and 1.5 Mya. © 2010 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 160 , 551–566.     

Uncovering an overlooked radiation: molecular phylogeny and biogeography of Madagascar's endemic river snails (Caenogastropoda: Pachychilidae: Madagasikara gen. nov.)

By analysing key morphological characters (with emphasis on shell, radula and stomach anatomy) and a partial fragment of the mitochondrial 16S rRNA gene (alignment length 860 bp), we examined patterns of diversity and differentiation of a previously overlooked radiation of Madagascan pachychilid freshwater gastropods. These analyses resulted in the discovery of three new species in addition to the two species that were already recognized. The complex nomenclatural and taxonomic implications are discussed and the finding of a viviparous reproductive mode in at least one among otherwise oviparous species is reported. Using a mitochondrial phylogeny that includes all currently accepted pachychilid genera and a strict molecular clock approach, we address the historical biogeography of the Madagascan species with respect to vicariant versus dispersalist biogeographical models. Using two alternative calibrations that were previously suggested for other gastropods, the molecular clock tree suggested that the origin of the Pachychilidae dates back to no more than 50 Mya, whereas the origin of the Madagascan lineage is estimated to date to a period between 15.6–31.5 Mya. These estimates are approximately concurrent with the dating of colonization events in a number of other Madagascan animal taxa. The pachychilid radiation on Madagascar appears not to be older than 3–5 Mya. Thus, although the global patterns of pachychilid distribution have earlier been interpreted to suggest a Gondwanan origin of the family, the present study does not support this postulate. Neither the topology of the molecular phylogeny, nor the timing of events as suggested from a molecular clock were found to be congruent with a vicariance scenario within the framework of Gondwanan fragmentation during the Mesozoic but, instead, imply overseas dispersal during the Cenozoic. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 867–894.     

The evolution of Metriorhynchoidea (mesoeucrocodylia,thalattosuchia): an integrated approach using geometric morphometrics,analysis of disparity,and biomechanics

Metriorhynchoid crocodylians represent the pinnacle of marine specialization within Archosauria. Not only were they a major component of the Middle Jurassic–Early Cretaceous marine ecosystems, but they provide further examples that extinct crocodilians did not all resemble their modern extant relatives. Here, we use a varied toolkit of techniques, including phylogenetic reconstruction, geometric morphometrics, diversity counts, discrete character disparity analysis, and biomechanical finite‐element analysis (FEA), to examine the macroevolutionary history of this clade. All analyses demonstrate that this clade became more divergent, in terms of biodiversity, form, and function, up until the Jurassic–Cretaceous boundary, after which there is no evidence for recovery or further radiations. A clear evolutionary trend towards hypercarnivory in Dakosaurus is supported by phylogenetic character optimization, morphometrics, and FEA, which also support specialized piscivory within Rhacheosaurus and Cricosaurus. Within Metriorhynchoidea, there is a consistent trend towards increasing marine specialization, with the hypermarine Cricosaurus exhibiting numerous convergences with other Mesozoic marine reptiles (e.g. loss of the deltopectoral crest and retracted external nares). In addition, biomechanics, morphometrics, and character‐disparity analyses consistently distinguish the two newly erected metriorhynchid subfamilies. This study illustrates that together with phylogeny, quantitative assessment of diversity, form, and function help elucidate the macroevolutionary pattern of fossil clades. © 2010 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 158 , 801–859.     

Molecular phylogenetics and biogeography of the eastern Asian–eastern North American disjunct Mitchella and its close relative Damnacanthus (Rubiaceae,Mitchelleae)

Mitchella is a small genus of the Rubiaceae with only two species. It is the only herbaceous semishrub of the family showing a disjunct distribution in eastern Asia and eastern North America, extending to Central America. Its phylogeny and biogeographical diversification remain poorly understood. In this study, we conducted phylogenetic and biogeographical analyses for Mitchella and its close relative Damnacanthus based on sequences of the nuclear internal transcribed spacer (ITS) and four plastid markers (rbcL, atpB‐rbcL, rps16 and trnL‐F). Mitchella is monophyletic, consisting of an eastern Asian M. undulata clade and a New World M. repens clade. Our results also support Michella as the closest relative to the eastern Asian Damnacanthus. The divergence time between the two intercontinental disjunct Mitchella species was dated to 7.73 Mya, with a 95% highest posterior density (HPD) of 3.14?12.53 Mya, using the Bayesian relaxed clock estimation. Ancestral area reconstructions suggest that the genus originated in eastern Asia. The semishrub Mitchella appears to have arisen from its woody ancestor in eastern Asia and then migrated to North America via the Bering land bridge in the late Miocene. © 2013 The Linnean Society of London     

Phylogenetic relationships within the corbiculate Apinae (Hymenoptera) and the evolution of eusociality

We provide a comparison of 520 base pairs of the mitochondrial cytochrome b gene from exemplars of the Meliponini, the Apini, the Bombini and the Euglossini to determine the phylogenetic relationships within the corbiculate Apinae. Our results strongly suggest (91–97% according to bootstrap resampling) that the Meliponini and the Bombini are sister groups. This finding agrees with those of other molecular studies, but is discordant with previous hypotheses based on morphology and the combination of molecular and morphological data. If the Bombini and Meliponini are sister groups and reversal of advanced eusociality is unlikely, then advanced eusociality arose twice within this clade. However, if reversion of eusociality occurred, then it is not possible to discriminate between single or dual origins of advanced eusociality within this group.     

Behavioral Phylogeny of Corbiculate Apidae (Hymenoptera; Apinae), with Special Reference to Social Behavior

The phylogenetic relationships among the four tribes of corbiculate bees (Euglossini, Bombini, Meliponini, and Apini) are controversial. There is substantial incongruence between morphological and molecular data, and the single origin of eusociality is questionable. The use of behavioral characters by previous workers has been restricted to some typological definitions, such as solitary and eusocial. Here, I expand the term "social" to 42 characters and present a tree based only on behavioral characters. The reconstructed relationships were similar to those observed in morphological and "total evidence" analyses, i.e., Euglossini + (Bombini + (Meliponini + Apini)), all of which support a single origin of eusociality.