Comparative metagenomics reveals expanded insights into intra- and interspecific variation among wild bee microbiomes.

The holobiont approach proposes that species are most fully understood within the context of their associated microbiomes, and that both host and microbial community are locked in a mutual circuit of co-evolutionary selection. Bees are an ideal group for this approach, as they comprise a critical group of pollinators that contribute to both ecological and agricultural health worldwide. Metagenomic analyses offer comprehensive insights into an organism's microbiome, diet, and viral load, but remain largely unapplied to wild bees. Here, we present metagenomic data from three species of carpenter bees sampled from around the globe, representative of the first ever carpenter bee core microbiome. Machine learning, co-occurrence, and network analyses reveal that wild bee metagenomes are unique to host species. Further, we find that microbiomes are likely strongly affected by features of their local environment, and feature evidence of plant pathogens previously known only in honey bees. Performing the most comprehensive comparative analysis of bee microbiomes to date we discover that microbiome diversity is inversely proportional to host species social complexity. Our study helps to establish some of the first wild bee hologenomic data while offering powerful empirical insights into the biology and health of vital pollinators.

Social divergence: molecular pathways underlying castes and longevity in a facultatively eusocial small carpenter bee.

Unravelling the evolutionary origins of eusocial life is a longstanding endeavour in the field of evolutionary-developmental biology. Descended from solitary ancestors, eusocial insects such as honeybees have evolved ontogenetic division of labour in which short-lived workers perform age-associated tasks, while a long-lived queen produces brood. It is hypothesized that (i) eusocial caste systems evolved through the co-option of deeply conserved genes and (ii) longevity may be tied to oxidative damage mitigation capacity. To date, however, these hypotheses have been examined primarily among only obligately eusocial corbiculate bees. We present brain transcriptomic data from a Japanese small carpenter bee, Ceratina japonica (Apidae: Xylocopinae), which demonstrates both solitary and eusocial nesting in sympatry and lives 2 or more years in the wild. Our dataset captures gene expression patterns underlying first- and second-year solitary females, queens and workers, providing an unprecedented opportunity to explore the molecular mechanisms underlying caste-antecedent phenotypes in a long-lived and facultatively eusocial bee. We find that C. japonica's queens and workers are underpinned by divergent gene regulatory pathways, involving many differentially expressed genes well-conserved among other primitively eusocial bee lineages. We also find support for oxidative damage reduction as a proximate mechanism of longevity in C. japonica.

Sociality sculpts similar patterns of molecular evolution in two independently evolved lineages of eusocial bees.

While it is well known that the genome can affect social behavior, recent models posit that social lifestyles can, in turn, influence genome evolution. Here, we perform the most phylogenetically comprehensive comparative analysis of 16 bee genomes to date: incorporating two published and four new carpenter bee genomes (Apidae: Xylocopinae) for a first-ever genomic comparison with a monophyletic clade containing solitary through advanced eusocial taxa. We find that eusocial lineages have undergone more gene family expansions, feature more signatures of positive selection, and have higher counts of taxonomically restricted genes than solitary and weakly social lineages. Transcriptomic data reveal that caste-affiliated genes are deeply-conserved; gene regulatory and functional elements are more closely tied to social phenotype than phylogenetic lineage; and regulatory complexity increases steadily with social complexity. Overall, our study provides robust empirical evidence that social evolution can act as a major and surprisingly consistent driver of macroevolutionary genomic change.

Social modularity: conserved genes and regulatory elements underlie caste-antecedent behavioural states in an incipiently social bee.

The evolutionary origins of advanced eusociality, one of the most complex forms of phenotypic plasticity in nature, have long been a focus within the field of sociobiology. Although eusocial insects are known to have evolved from solitary ancestors, sociogenomic research among incipiently social taxa has only recently provided empirical evidence supporting theories that modular regulation and deeply conserved genes may play important roles in both the evolutionary emergence and elaboration of insect sociality. There remains, however, a paucity of data to further test the biological reality of these and other evolutionary theories among taxa in the earliest stages of social evolution. Here, we present brain transcriptomic data from the incipiently social small carpenter bee, Ceratina calcarata, which captures patterns of cis-regulation and gene expression associated with female maturation, and underlying two well-defined behavioural states, foraging and guarding, concurrently demonstrated by mothers and daughters during early autumn. We find that an incipiently social nest environment may dramatically affect gene expression. We further reveal foraging and guarding behaviours to be putatively caste-antecedent states in C. calcarata, and offer strong empirical support for the operation of modular regulation, involving deeply conserved and differentially expressed genes in the expression of early social forms.

Development of Multiple Polymorphic Microsatellite Markers for Ceratina calcarata (Hymenoptera: Apidae) Using Genome-Wide Analysis.

The small carpenter bee, Ceratina calcarata (Robertson), is a widespread native pollinator across eastern North America. The behavioral ecology and nesting biology of C. calcarata has been relatively well-studied and the species is emerging as a model organism for both native pollinator and social evolution research. C. calcarata is subsocial: reproductively mature females provide extended maternal care to their brood. As such, studies of C. calcarata may also reveal patterns of relatedness and demography unique to primitively social Hymenoptera. Here, we present 21 microsatellite loci, isolated from the recently completed C. calcarata genome. Screening in 39 individuals across their distribution revealed that no loci were in linkage disequilibrium, nor did any deviate significantly from Hardy-Weinberg following sequential Bonferroni correction. Allele count ranged from 2 to 14, and observed and expected heterozygosities ranged from 0.08 to 0.82 (mean 0.47) and 0.26 to 0.88 (mean 0.56), respectively. These markers will enable studies of population-wide genetic structuring across C. calcarata's distribution. Such tools will also allow for exploration of between and within-colony relatedness in this subsocial native pollinator.