Neotropical bee microbiomes point to a fragmented social core and strong species-level effects

BACKGROUND

Individuals that band together create new ecological opportunities for microorganisms. In vertical transmission, theory predicts a conserved microbiota within lineages, especially social bees. Bees exhibit solitary to social behavior among and/or within species, while life cycles can be annual or perennial. Bee nests may be used over generations or only once, and foraging ecology varies widely. To assess which traits are associated with bee microbiomes, we analyzed microbial diversity within solitary and social bees of Apidae, Colletidae, and Halictidae, three bee families in Panama's tropical forests. Our analysis considered the microbiome of adult gut contents replicated through time, localities, and seasons (wet and dry) and included bee morphology and comparison to abdominal (dissected) microbiota. Diversity and distribution of tropical bee microbes (TBM) within the corbiculate bee clade were emphasized.

RESULTS

We found the eusocial corbiculate bees tended to possess a more conserved gut microbiome, attributable to vertical transmission, but microbial composition varied among closely related species. Euglossine bees (or orchid bees), corbiculates with mainly solitary behavior, had more variable gut microbiomes. Their shorter-tongued and highly seasonal species displayed greater diversity, attributable to flower-visiting habits. Surprisingly, many stingless bees, the oldest corbiculate clade, lacked bacterial genera thought to predate eusociality, while several facultatively social, and solitary bee species possessed those bacterial taxa. Indeed, nearly all bee species displayed a range of affinities for single or multiple variants of the "socially associated" bacterial taxa, which unexpectedly demonstrated high sequence variation.

CONCLUSIONS

Taken together, these results call into question whether specific bacterial associates facilitate eusocial behavior, or are subsequently adopted, or indicate frequent horizontal transmission between perennial eusocial colonies and other social, facultatively social, and solitary bees. Video Abstract.

Stingless Bee (Apidae: Apinae: Meliponini) Ecology

Stingless bees form perennial colonies of honey-making insects. The >600 species of stingless bees, mainly Neotropical, live throughout tropical latitudes. Foragers influence floral biology, plant reproduction, microbe dispersal, and diverse ecosystem functions. As tropical forest residents since the upper Cretaceous, they have had a long evolutionary history without competition from honey bees. Most stingless bees are smaller than any Apis species and recruit nest mates to resources, while their defense strategies exclude stinging behavior but incorporate biting. Stingless bees have diversified ecologically; excel in nesting site selection and mutualisms with plants, arthropods, and microbes; and display opportunism, including co-opting plant defenses. As their biology becomes better known, applications to human endeavors are imposing selective pressures from exploitation and approaches to conservation that entail colony extraction from wildlands. Although some meliponines can adjust to new conditions, their populations shall require tropical diversity for survival and reproduction.

Size-dependent Scaling of Stingless Bee Flight Metabolism Reveals an Energetic Benefit to Small Body Size

Understanding the effect of body size on flight costs is critical for development of models of aerodynamics and animal energetics. Prior scaling studies that have shown that flight costs scale hypometrically have focused primarily on larger (> 100 mg) insects and birds, but most flying species are smaller. We studied the flight physiology of thirteen stingless bee species over a large range of body sizes (1-115 mg). Metabolic rate during hovering scaled hypermetrically (scaling slope = 2.11). Larger bees had warm thoraxes while small bees were nearly ecothermic; however, even controlling for body temperature variation, flight metabolic rate scaled hypermetrically across this clade. Despite having a lower mass-specific metabolic rate during flight, smaller bees could carry the same proportional load. Wingbeat frequency did not vary with body size, in contrast to most studies that find wingbeat frequency increases as body size decreases. Smaller stingless bees have greater relative forewing surface area which may help them reduce the energy requirements needed to fly. Further, we hypothesize that the relatively larger heads of smaller species may change their body pitch in flight. Synthesizing across all flying insects, we demonstrate that the scaling of flight metabolic rate changes from hypermetric to hypometric at approximately 58 mg body mass with hypermetic scaling below (slope = 1.2) and hypometric scaling (slope = 0.67) above 58 mg in body mass. The reduced cost of flight likely provides selective advantages for the evolution of small body size in insects. The biphasic scaling of flight metabolic rates and wingbeat frequencies in insects supports the hypothesis that the scaling of metabolic rate is closely related to the power requirements of locomotion and cycle frequencies.

The importance of insects on land and in water: a tropical view

Tropical insects are astonishingly diverse and abundant yet receive only marginal scientific attention. In natural tropical settings, insects are involved in regulating and supporting ecosystem services including seed dispersal, pollination, organic matter decomposition, nutrient cycling, herbivory, food webs and water quality, which in turn help fulfill UN Sustainable Development Goals (SDGs). Current and future global changes that affect insect diversity and distribution could disrupt key ecosystem services and impose important threats on ecosystems and human well-being. A significant increase in our knowledge of tropical insect roles in ecosystem processes is thus vital to ensure sustainable development on a rapidly changing planet.

Changing of the guard: mixed specialization and flexibility in nest defense (Tetragonisca angustula)

Task allocation is a central challenge of collective behavior in a variety of group-living species, and this is particularly the case for the allocation of social insect workers for group defense. In social insects, both benefits and considerable costs are associated with the production of specialized soldiers. We asked whether colonies mitigate costs of production of specialized soldiers by simultaneously employing behavioral flexibility in nonspecialist workers that can augment defense capabilities at short time scales. We studied colonies of the stingless bee Tetragonisca angustula, a species that has 2 discrete nest-guarding tasks typically performed by majors: hovering guarding and standing guarding. Majors showed age polyethism across nest-guarding tasks, first hovering and then changing to the task of standing guarding after 1 week. Colonies were also able to reassign minors to guarding tasks when majors were experimentally removed. Replacement guards persisted in nest defense tasks until colonies produced enough majors to return to their initial state. Tetragonisca angustula colonies thus employed a coordinated set of specialization strategies in nest defense: morphologically specialized soldiers, age polyethism among soldiers within specific guarding tasks, and rapid flexible reallocation of nonspecialists to guarding during soldier loss. This mixed strategy achieves the benefits of a highly specialized defensive force while maintaining the potential for rapid reinforcement when soldiers are lost or colonies face unexpectedly intense attack.

Colony-level non-associative plasticity of alarm responses in the stingless honey bee, Tetragonisca angustula

In ants, bees, and other social Hymenoptera alarm pheromones are widely employed to coordinate colony nest defense. In that context, alarm pheromones elicit innate species-specific defensive behaviors. Therefore, in terms of classical conditioning, an alarm pheromone could act as an unconditioned stimulus (US). Here we test this hypothesis by establishing whether repeated exposure to alarm pheromone in different testing contexts modifies the alarm response. We evaluate colony level alarm responses in the stingless bee, Tetragonisca angustula, which has a morphologically distinct guard caste. First, we describe the overall topology of defense behaviors in the presence of an alarm pheromone. Second, we show that repeated, regular exposure to synthetic alarm pheromone reduces different components of the alarm response, and memory of that exposure decays over time. This observed decrease followed by recovery occurs over different time frames and is consistent with behavioral habituation. We further tested whether the alarm pheromone can act as a US to classically condition guards to modify their defense behaviors in the presence of a novel (conditioned) stimulus (CS). We found no consistent changes in the response to the CS. Our study demonstrates the possibility that colony-level alarm responses can be adaptively modified by experience in response to changing environmental threats. Further studies are now needed to reveal the extent of these habituation-like responses in regard to other pheromones, the potential mechanisms that underlie this phenomenon, and the range of adaptive contexts in which they function at the colony level.

Setting the pace of life: membrane composition of flight muscle varies with metabolic rate of hovering orchid bees

Patterns of metabolic rate variation have been documented extensively in animals, but their functional basis remains elusive. The membrane pacemaker hypothesis proposes that the relative abundance of polyunsaturated fatty acids in membrane phospholipids sets the metabolic rate of organisms. Using species of tropical orchid bees spanning a 16-fold range in body size, we show that the flight muscles of smaller bees have more linoleate (%18 : 3) and stearate (%18 : 0), but less oleate (%18 : 1). More importantly, flight metabolic rate (FlightMR) varies with the relative abundance of 18 : 3 according to the predictions of the membrane pacemaker hypothesis. Although this relationship was found across large differences in metabolic rate, a direct association could not be detected when taking phylogeny and body mass into account. Higher FlightMR, however, was related to lower %16 : 0, independent of phylogeny and body mass. Therefore, this study shows that flight muscle membrane composition plays a significant role in explaining diversity in FlightMR, but that body mass and phylogeny are other factors contributing to their variation. Multiple factors are at play to modulate metabolic capacity, and changing membrane composition can have gradual and stepwise effects to achieve a new range of metabolic rates. Orchid bees illustrate the correlated evolution between membrane composition and metabolic rate, supporting the functional link proposed in the membrane pacemaker hypothesis.

Impact of Striped-Squirrel Nectar-Robbing Behaviour on Gender Fitness in Alpinia roxburghii Sweet (Zingiberaceae)

Nectar-robbing has the potential to strongly affect male and female reproductive fitness of plants. One example of nectar theft is that shown by striped-squirrels (Tamiops swinhoei) on a number of ginger species, including Alpinia roxburghii and A. kwangsiensis (Zingiberaceae). In this study, we used a fluorescent dye as a pollen analogue, and measured fruit and seed output, to test the effect of squirrel nectar-robbing on A. roxburghii reproductive fitness. Pollen transfer between robbed and unrobbed flowers was assessed by comparing 60 randomly established plots containing robbed and unrobbed flowers. The frequency of squirrel robbing visits and broken styles were recorded from a number of flowers for five consecutive days. Two bee species (Bombus eximius and Apis cerana), were the primary pollinators, and their visitation frequency was recorded for six consecutive days. The results showed that fluorescent powder from unrobbed flowers was dispersed further, and to a greater number of flowers than that placed on robbed flowers. Additionally, robbing flowers caused significant damage to reproductive organs, resulting in lower fruit and seed sets in robbed than in unrobbed flowers and influencing both male and female fitness. The frequency of the primary pollinator visits (B. eximius) was significantly higher for unrobbed plants than for robbed plants. The present study clearly shows the negative impact of squirrel robbing on A. roxburghii male reproductive fitness and neutral impact on female reproductive fitness.

Arthropod Distribution in a Tropical Rainforest: Tackling a Four Dimensional Puzzle

Quantifying the spatio-temporal distribution of arthropods in tropical rainforests represents a first step towards scrutinizing the global distribution of biodiversity on Earth. To date most studies have focused on narrow taxonomic groups or lack a design that allows partitioning of the components of diversity. Here, we consider an exceptionally large dataset (113,952 individuals representing 5,858 species), obtained from the San Lorenzo forest in Panama, where the phylogenetic breadth of arthropod taxa was surveyed using 14 protocols targeting the soil, litter, understory, lower and upper canopy habitats, replicated across seasons in 2003 and 2004. This dataset is used to explore the relative influence of horizontal, vertical and seasonal drivers of arthropod distribution in this forest. We considered arthropod abundance, observed and estimated species richness, additive decomposition of species richness, multiplicative partitioning of species diversity, variation in species composition, species turnover and guild structure as components of diversity. At the scale of our study (2 km of distance, 40 m in height and 400 days), the effects related to the vertical and seasonal dimensions were most important. Most adult arthropods were collected from the soil/litter or the upper canopy and species richness was highest in the canopy. We compared the distribution of arthropods and trees within our study system. Effects related to the seasonal dimension were stronger for arthropods than for trees. We conclude that: (1) models of beta diversity developed for tropical trees are unlikely to be applicable to tropical arthropods; (2) it is imperative that estimates of global biodiversity derived from mass collecting of arthropods in tropical rainforests embrace the strong vertical and seasonal partitioning observed here; and (3) given the high species turnover observed between seasons, global climate change may have severe consequences for rainforest arthropods.

Detecting insect pollinator declines on regional and global scales

Recently there has been considerable concern about declines in bee communities in agricultural and natural habitats. The value of pollination to agriculture, provided primarily by bees, is >$200 billion/year worldwide, and in natural ecosystems it is thought to be even greater. However, no monitoring program exists to accurately detect declines in abundance of insect pollinators; thus, it is difficult to quantify the status of bee communities or estimate the extent of declines. We used data from 11 multiyear studies of bee communities to devise a program to monitor pollinators at regional, national, or international scales. In these studies, 7 different methods for sampling bees were used and bees were sampled on 3 different continents. We estimated that a monitoring program with 200-250 sampling locations each sampled twice over 5 years would provide sufficient power to detect small (2-5%) annual declines in the number of species and in total abundance and would cost U.S.$2,000,000. To detect declines as small as 1% annually over the same period would require >300 sampling locations. Given the role of pollinators in food security and ecosystem function, we recommend establishment of integrated regional and international monitoring programs to detect changes in pollinator communities.

Color dependent learning in restrained Africanized honey bees

Associative color learning has been demonstrated to be very poor using restrained European honey bees unless the antennae are amputated. Consequently, our understanding of proximate mechanisms in visual information processing is handicapped. Here we test learning performance of Africanized honey bees under restrained conditions with visual and olfactory stimulation using the proboscis extension response (PER) protocol. Restrained individuals were trained to learn an association between a color stimulus and a sugar-water reward. We evaluated performance for "absolute" learning (learned association between a stimulus and a reward) and "discriminant" learning (discrimination between two stimuli). Restrained Africanized honey bees (AHBs) readily learned the association of color stimulus for both blue and green LED stimuli in absolute and discriminatory learning tasks within 7 presentations, but not with violet as the rewarded color. Additionally, 24-hour memory improved considerably during the discrimination task, compared to absolute association (15%-55%). We found that antennal amputation was unnecessary and reduced performance in AHBs. Thus color learning can now be studied using the PER protocol with intact AHBs. This finding opens the way toward investigating visual and multimodal learning with application of neural techniques commonly used in restrained honey bees.

Arthropod diversity in a tropical forest

Most eukaryotic organisms are arthropods. Yet, their diversity in rich terrestrial ecosystems is still unknown. Here we produce tangible estimates of the total species richness of arthropods in a tropical rainforest. Using a comprehensive range of structured protocols, we sampled the phylogenetic breadth of arthropod taxa from the soil to the forest canopy in the San Lorenzo forest, Panama. We collected 6144 arthropod species from 0.48 hectare and extrapolated total species richness to larger areas on the basis of competing models. The whole 6000-hectare forest reserve most likely sustains 25,000 arthropod species. Notably, just 1 hectare of rainforest yields >60% of the arthropod biodiversity held in the wider landscape. Models based on plant diversity fitted the accumulated species richness of both herbivore and nonherbivore taxa exceptionally well. This lends credence to global estimates of arthropod biodiversity developed from plant models.

Large-range movements of neotropical orchid bees observed via radio telemetry

Neotropical orchid bees (Euglossini) are often cited as classic examples of trapline-foragers with potentially extensive foraging ranges. If long-distance movements are habitual, rare plants in widely scattered locations may benefit from euglossine pollination services. Here we report the first successful use of micro radio telemetry to track the movement of an insect pollinator in a complex and forested environment. Our results indicate that individual male orchid bees (Exaerete frontalis) habitually use large rainforest areas (at least 42-115 ha) on a daily basis. Aerial telemetry located individuals up to 5 km away from their core areas, and bees were often stationary, for variable periods, between flights to successive localities. These data suggest a higher degree of site fidelity than what may be expected in a free living male bee, and has implications for our understanding of biological activity patterns and the evolution of forest pollinators.

Allometric scaling of flight energetics in Panamanian orchid bees: a comparative phylogenetic approach

The relationship between body size and flight energetics was studied in the clade of tropical orchid bees, in order to investigate energy metabolism and evolution. Body mass, which varied from 47 to 1065 mg, was found to strongly affect hovering flight mass-specific metabolic rates, which ranged from 114 ml CO(2) h(-1) g(-1) in small species to 37 ml CO(2) h(-1) g(-1) in large species. Similar variation of wingbeat frequency in hovering flight occurred among small to large species, and ranged from 250 to 86 Hz. The direct relationship between such traits was studied by the comparative method of phylogenetically independent contrasts (PIC), using a new molecular phylogeny generated from the cytochrome b gene partial sequences. We found wingbeat frequency variation is satisfactorily explained by variation in wing loading, after corrections for body mass and phylogeny. The correlated evolution of mass-specific metabolic rate, wingbeat frequency and wing loading was also revealed after correcting for phylogeny and body mass. Further, the effect of body size on flight energetics can be understood in terms of a relationship between wing form and kinematics, which directly influence and explain the scaling of metabolic rate in this group of bees.

Energy metabolism in orchid bee flight muscles: carbohydrate fuels all

The widely accepted idea that bees fuel flight through the oxidation of carbohydrate is based on studies of only a few species. We tested this hypothesis as part of our research program to investigate the size-dependence of flight energetics in Panamanian orchid bees. We succeeded in measuring rates of O(2) consumption and CO(2) production in vivo during hovering flight, as well as maximal activities (V(max) values) in vitro of key enzymes in flight muscle energy metabolism in nine species belonging to four genera. Respiratory quotients (ratios of rates of CO(2) production to O(2) consumption) in all nine species are close to 1.0. This indicates that carbohydrate is the main fuel used for flight. Trehalase, glycogen phosphorylase and hexokinase activities are sufficient to account for the glycolytic flux rates estimated from rates of CO(2) production. High activities of other glycolytic enzymes, as well as high activities of mitochondrial oxidative enzymes, are consistent with the estimated rates of carbohydrate-fueled oxidative metabolism. In contrast, hydroxyacylCoA dehydrogenase, an enzyme involved in fatty acid oxidation, was not detectable in any species. Thoracic homogenates displayed ADP-stimulated oxidition of pyruvate + proline, but did not oxidize palmitoyl l-carnitine + proline as substrates. A metabolic map, based on data reported herein and information from the literature, is presented. The evidence available supports the hypothesis that carbohydrate serves as the main fuel for flight in bees.

Allometric scaling of flight energetics in orchid bees: evolution of flux capacities and flux rates

The evolution of metabolic pathways involved in energy production was studied in the flight muscles of 28 species of orchid bees. Previous work revealed that wingbeat frequencies and mass-specific metabolic rates decline in parallel by threefold as body mass increases interspecifically over a 20-fold range. We investigated the correlated evolution of metabolic rates during hovering flight and the flux capacities, i.e. Vmaxvalues, of flight muscle enzymes involved in substrate catabolism, the Krebs cycle and the electron transport chain. Vmax at the hexokinase (HK) step scales allometrically with an exponent almost identical to those obtained for wingbeat frequency and mass-specific metabolic rate. Analysis of this relationship using phylogenetically independent contrasts supports the hypothesis of correlated evolution between HK activity and mass-specific metabolic rate. Although other enzymes scale allometrically with respect to body mass, e.g. trehalase, glycogen phosphorylase and citrate synthase, no other enzyme activities were correlated with metabolic rate after controlling for phylogenetic relatedness. Pathway flux rates were used with enzyme Vmax values to estimate fractional velocities(fraction of Vmax at which enzymes operate) for various reactions to gain insights into enzyme function and how this varies with body mass. Fractional velocity is highly conserved across species at the HK step,but varied at all other steps examined. These results are discussed in the context of the regulation and evolution of pathways of energy metabolism.

Long-distance gene flow and cross-Andean dispersal of lowland rainforest bees (Apidae: Euglossini) revealed by comparative mitochondrial DNA phylogeography

Euglossine bees (Apidae; Euglossini) exclusively pollinate hundreds of orchid species and comprise up to 25% of bee species richness in neotropical rainforests. As one of the first studies of comparative phylogeography in a neotropical insect group, we performed a mitochondrial DNA (mtDNA)-based analysis of 14 euglossine species represented by populations sampled across the Andes and/or across the Amazon basin. The mtDNA divergences within species were consistently low; across the 12 monophyletic species the mean intraspecific divergence among haplotypes was 0.9% (range of means, 0-1.9%). The cytochrome oxidase 1 (CO1) divergence among populations separated by the Andes (N = 11 species) averaged 1.1% (range 0.0-2.0%). The mtDNA CO1 data set displayed homogeneous rates of nucleotide substitution, permitting us to infer dispersal across the cordillera long after the final Andean uplift based on arthropod molecular clocks of 1.2-1.5% divergence per million years. Gene flow across the 3000-km breadth of the Amazon basin was inferred from identical cross-Amazon haplotypes found in five species. Although mtDNA haplotypes for 12 of the 14 euglossine species were monophyletic, a reticulate CO1 phylogeny was recovered in Euglossa cognata and E. mixta, suggesting large ancestral populations and recent speciation. Reference to closely related outgroups suggested recent speciation for the majority of species. Phylogeographical structure across a broad spatial scale is weaker in euglossine bees than in any neotropical group previously examined, and may derive from a combination of Quaternary speciation, population expansion and/or long-distance gene flow.

Use of diploid male frequency data as an indicator of pollinator decline

Pollination deficits in agricultural and natural systems are suggestive of large reductions in pollinator populations. However, actual declines are difficult to demonstrate using census data. Here, we show census data to be misleading because many abundant pollinators exhibit high levels of production of sterile diploid males usually found only in small inbred hymenopteran populations; Euglossa imperialis exhibits high levels of diploid male production induced by low effective population sizes (Ne approximately 15), despite being the most abundant orchid bee in lowland tropical forests in Panama. We caution that although some pollinators appear abundant on the basis of census data, their long-term persistence may be highly tenuous based on genetic evidence. We propose the use of diploid male frequency data as a metric for assessing the sustainability of bee populations.

Odor compound detection in male euglossine bees

Male euglossine bees collect fragrances from various sources, which they store and use for as yet unknown purposes. They are attracted, often specifically, to single odor compounds and blends thereof. We used gas chromatography with electroantennographic detection (GC-EAD) and electroantennography (EAG) to investigate the response to 8 odor compounds by males of two euglossine species, Euglossa cybelia Moure and Eulaema polychroma (Mocsàry). In E. cybelia, we recorded EAD reactions in response to 1,8-cineole, methyl benzoate, benzyl actetate, methyl salicylate, eugenol, and methyl cinnamate. E. polychroma responded to the same compounds in EAG experiments, while (1s)(-)alpha-pinene and beta-pinene failed to trigger EAD or EAG responses in the bees. Blends of two compounds triggered larger responses than single compounds in EAG experiments with E. polychroma, however, when alpha-pinene was added, reactions decreased. In the light of existing data on the bees' behavior towards these odor compounds, our work indicates that both peripheral and central nervous processes influence the attraction of euglossine bees to odors.

The value of bees to the coffee harvest

The self-pollinating African shrub Coffea arabica, a pillar of tropical agriculture, was considered to gain nothing from insect pollinators. But I show here that naturalized, non-native honeybees can augment pollination and boost crop yields by over 50%. These findings, together with world coffee-harvest statistics and results from field studies of organically shade-grown coffee, indicate that coffee plants would benefit from being grown in habitats that are suitable for sustaining valuable pollinators.

A Framework for Melittology

The Bees of the World . Charles D. Michener. Johns Hopkins University Press, Baltimore, MD, 2000. 929 pp. $135, £93. ISBN 0-8018-6133-0.

As the capstone of his research over more than a half century, Michener offers a comprehensive taxonomic and phylogenetic framework for bee biology.

Obligate Necrophagy in a Social Bee

The social bee Trigona hypogea uses carrion instead of pollen as a protein source. Nests lack stored pollen, pollen is absent in larval provisions, and corbiculae for pollen transport on worker hind legs are reduced. Glandular secretions of 20 percent protein content appear to replace stored pollen. Toothed mandibles, pheromonal recruitment to resources, and aggressive foraging behavior facilitated evolution of necrophagy in this tropical forest bee lineage.