Genome assembly and variant analysis of two Saccharomyces cerevisiae strains isolated from stingless bee pollen

Products from stingless bees are rich reservoirs of microbial diversity, including yeasts with fermentative potential. Previously, two Saccharomyces cerevisiae strains, JP14 and IP9, were isolated from Jataí (Tetragonisca angustula) and Iraí (Nannotrigona testaceicornis) bees, respectively, aiming at mead production. Both strains presented great osmotic and sulfite tolerance, and ethanol production, although they have a high free amino nitrogen demand. Herein, their genomes were sequenced, assembled, and annotated, and the variants were compared to the S. cerevisiae S288c reference strain. The final assembly of IP9 and JP14 presented high N50 and BUSCO scores, and more than 6430 protein-coding genes. Additionally, nQuire predicted the ploidy of IP9 as diploid, but the results were not enough to determine the ploidy of JP14. The mitochondrial genomes of IP9 and JP14 presented the same gene content as S288c but the genes were rearranged and fragmented in different patterns. Meanwhile, the genes with mutations of high impact (e.g., indels, gain of stop codon) for both yeasts were enriched for transmembrane transport, electron transfer, oxidoreductase, heme binding, fructose, mannose, and glucose transport, activities related to the respiratory chain and sugar metabolism. The IP9 strain presented copy number gains in genes related to sugar transport and cell morphogenesis; in JP14, genes were enriched for disaccharide metabolism and transport, response to reactive oxygen species, and polyamine transport. On the other hand, IP9 presented copy number losses related to disaccharide, thiamine, and aldehyde metabolism, while JP14 presented depletions related to disaccharide, oligosaccharide, asparagine, and aspartate metabolism. Notably, both strains presented a killer toxin gene, annotated from the assembling of unmapped reads, representing a potential mechanism for the control of other microorganisms population in the environment. Therefore, the annotated genomes of JP14 and IP9 presented a high selective pressure for sugar and nitrogen metabolism and stress response, consistent with their isolation source and fermentative properties.

Utilization of systematic error-assessment software to improve phylogenetic accuracy

Unlike classical systems based on the use of morphological data, modern phylogenetic analyses use genetic information to construct phylogenetic trees. Ongoing research in the field of phylogenetics is evaluating the accuracy of phylogenetic estimation results and the reliability of phylogenetic trees to explain evolutionary relationships. Recently, the probability of stochastic errors in large-scale phylogenetic datasets has decreased, while the probability of systematic errors has increased. Therefore, before constructing a phylogenetic tree, it is necessary to assess the causes of systematic bias to improve the accuracy of phylogenetic estimates. We performed analyses of three datasets (Terebelliformia, Daphniid, and Glires clades) using bioinformatics software to assess systematic error and improve phylogenetic tree accuracy. Then, we proposed a combination of systematic biases capable of discerning the most suitable gene markers within a series of taxa and generating conflicting phylogenetic topologies. Our findings will help improve the reliability of phylogenetic software to estimate phylogenies more accurately by exploiting systematic bias.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

Further evidences of an emerging stingless bee-yeast symbiosis

Symbiotic interactions between microorganisms and social insects have been described as crucial for the maintenance of these multitrophic systems, as observed for the stingless bee Scaptotrigona depilis and the yeast Zygosaccharomyces sp. SDBC30G1. The larvae of S. depilis ingest fungal filaments of Zygosaccharomyces sp. SDBC30G1 to obtain ergosterol, which is the precursor for the biosynthesis of ecdysteroids that modulate insect metamorphosis. In this work, we find a similar insect-microbe interaction in other species of stingless bees. We analyzed brood cell samples from 19 species of stingless bees collected in Brazil. The osmophilic yeast Zygosaccharomyces spp. was isolated from eight bee species, namely Scaptotrigona bipunctata, S. postica, S. tubiba, Tetragona clavipes, Melipona quadrifasciata, M. fasciculata, M. bicolor, and Partamona helleri. These yeasts form pseudohyphae and also accumulate ergosterol in lipid droplets, similar to the pattern observed for S. depilis. The phylogenetic analyses including various Zygosaccharomyces revealed that strains isolated from the brood cells formed a branch separated from the previously described Zygosaccharomyces species, suggesting that they are new species of this genus and reinforcing the symbiotic interaction with the host insects.

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.

Positive selection in cytochrome P450 genes is associated with gonad phenotype and mating strategy in social bees

The honey bee, Apis mellifera differs from all other social bees in its gonad phenotype and mating strategy. Honey bee queens and drones have tremendously enlarged gonads, and virgin queens mate with several males. In contrast, in all the other bees, the male and female gonads are small, and the females mate with only one or very few males, thus, suggesting an evolutionary and developmental link between gonad phenotype and mating strategy. RNA-seq comparisons of A. mellifera larval gonads revealed 870 genes as differentially expressed in queens versus workers and drones. Based on Gene Ontology enrichment we selected 45 genes for comparing the expression levels of their orthologs in the larval gonads of the bumble bee Bombus terrestris and the stingless bee, Melipona quadrifasciata, which revealed 24 genes as differentially represented. An evolutionary analysis of their orthologs in 13 solitary and social bee genomes revealed four genes with evidence of positive selection. Two of these encode cytochrome P450 proteins, and their gene trees indicated a lineage-specific evolution in the genus Apis, indicating that cytochrome P450 genes may be involved in the evolutionary association of polyandry and the exaggerated gonad phenotype in social bees.

Deep Divergence and Genomic Diversification of Gut Symbionts of Neotropical Stingless Bees

Social bees harbor conserved gut microbiotas that may have been acquired in a common ancestor of social bees and subsequently codiversified with their hosts. However, most of this knowledge is based on studies on the gut microbiotas of honey bees and bumblebees. Much less is known about the gut microbiotas of the third and most diverse group of social bees, the stingless bees. Specifically, the absence of genomic data from their microbiotas presents an important knowledge gap in understanding the evolution and functional diversity of the social bee microbiota. Here, we combined community profiling with culturing and genome sequencing of gut bacteria from six neotropical stingless bee species from Brazil. Phylogenomic analyses show that most stingless bee gut isolates form deep-branching sister clades of core members of the honey bee and bumblebee gut microbiota with conserved functional capabilities, confirming the common ancestry and ecology of their microbiota. However, our bacterial phylogenies were not congruent with those of the host, indicating that the evolution of the social bee gut microbiota was not driven by strict codiversification but included host switches and independent symbiont gain and losses. Finally, as reported for the honey bee and bumblebee microbiotas, we found substantial genomic divergence among strains of stingless bee gut bacteria, suggesting adaptation to different host species and glycan niches. Our study offers first insights into the genomic diversity of the stingless bee microbiota and highlights the need for broader samplings to understand the evolution of the social bee gut microbiota. IMPORTANCE Stingless bees are the most diverse group of the corbiculate bees and represent important pollinator species throughout the tropics and subtropics. They harbor specialized microbial communities in their gut that are related to those found in honey bees and bumblebees and that are likely important for bee health. Few bacteria have been cultured from the gut of stingless bees, which has prevented characterization of their genomic diversity and functional potential. Here, we established cultures of major members of the gut microbiotas of six stingless bee species and sequenced their genomes. We found that most stingless bee isolates belong to novel bacterial species distantly related to those found in honey bees and bumblebees and encoding similar functional capabilities. Our study offers a new perspective on the evolution of the social bee gut microbiota and presents a basis for characterizing the symbiotic relationships between gut bacteria and stingless bees.

Experimental disruption of social structure reveals totipotency in the orchid bee, Euglossa dilemma

Eusociality has evolved multiple times across the insect phylogeny. Social insects with greater levels of social complexity tend to exhibit specialized castes with low levels of individual phenotypic plasticity. In contrast, species with simple social groups may consist of totipotent individuals that transition among behavioral and reproductive states. However, recent work has shown that in simple social groups, there can still be constraint on individual plasticity, caused by differences in maternal nourishment or social interaction. It is not well understood how these constraints arise, ultimately leading to the evolution of nonreproductive workers. Some species of orchid bees form social groups of a dominant and-one to two subordinate helpers where all individuals are reproductive. Females can also disperse to start their own nest as a solitary foundress, which includes a nonreproductive phase characterized by ovary inactivation, not typically expressed by subordinates. Little is known about individual flexibility across these trajectories. Here, using the orchid bee Euglossa dilemma, we assess the plasticity of subordinate helpers, finding that they are capable of the same behavioral, physiological, transcriptomic, and chemical changes seen in foundresses. Our results suggest that the lack of nonreproductive workers in E. dilemma is not due to a lack of subordinate plasticity.

Convergent Loss of Prothoracicotropic Hormone, A Canonical Regulator of Development, in Social Bee Evolution

The evolution of insect sociality has repeatedly involved changes in developmental events and their timing. Here, we propose the hypothesis that loss of a canonical regulator of moulting and metamorphosis, prothoracicotropic hormone (PTTH), and its receptor, Torso, is associated with the evolution of sociality in bees. Specifically, we posit that the increasing importance of social influences on early developmental timing in social bees has led to their decreased reliance on PTTH, which connects developmental timing with abiotic cues in solitary insects. At present, the evidence to support this hypothesis includes the absence of genes encoding PTTH and Torso from all fully-sequenced social bee genomes and its presence in all available genomes of solitary bees. Based on the bee phylogeny, the most parsimonious reconstruction of evolutionary events is that this hormone and its receptor have been lost multiple times, across independently social bee lineages. These gene losses shed light on possible molecular and cellular mechanisms that are associated with the evolution of social behavior in bees. We outline the available evidence for our hypothesis, and then contextualize it in light of what is known about developmental cues in social and solitary bees, and the multiple precedences of major developmental changes in social insects.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

The Emerging Phylogenetic Perspective on the Evolution of Actinopterygian Fishes

The emergence of a new phylogeny of ray-finned fishes at the turn of the twenty-first century marked a paradigm shift in understanding the evolutionary history of half of living vertebrates. We review how the new ray-finned fish phylogeny radically departs from classical expectations based on morphology. We focus on evolutionary relationships that span the backbone of ray-finned fish phylogeny, from the earliest divergences among teleosts and nonteleosts to the resolution of major lineages of Percomorpha. Throughout, we feature advances gained by the new phylogeny toward a broader understanding of ray-finned fish evolutionary history and the implications for topics that span from the genetics of human health to reconsidering the concept of living fossils. Additionally, we discuss conceptual challenges that involve reconciling taxonomic classification with phylogenetic relationships and propose an alternate higher-level classification for Percomorpha. Our review highlights remaining areas of phylogenetic uncertainty and opportunities for comparative investigations empowered by this new phylogenetic perspective on ray-finned fishes.

Exploring systematic biases, rooting methods and morphological evidence to unravel the evolutionary history of the genus Ficus (Moraceae)

Despite many attempts in the Sanger sequencing era, the phylogeny of fig trees remains unresolved, which limits our ability to analyze the evolution of key traits that may have contributed to their evolutionary and ecological success. We used restriction-site-associated DNA sequencing (c. 420 kb) and 102 morphological characters to elucidate the relationships between 70 species of Ficus. To increase phylogenetic information for higher-level relationships, we targeted conserved regions and assembled paired reads into long loci to enable the retrieval of homologous loci in outgroup genomes. We compared morphological and molecular results to highlight discrepancies and reveal possible inference bias. For the first time, we recovered a monophyletic subgenus Urostigma (stranglers) and a clade with all gynodioecious Ficus. However, we show, with a new approach based on iterative principal component analysis, that it is not (and will probably never be) possible to homogenize evolutionary rates and GC content for all taxa before phylogenetic inference. Four competing positions for the root of the molecular tree are possible. The placement of section Pharmacosycea as sister to other fig trees is not supported by morphological data and considered a result of a long-branch attraction artefact to the outgroups. Regarding morphological features and indirect evidence from the pollinator tree of life, the topology that divides Ficus into monoecious versus gynodioecious species appears most plausible. It seems most likely that the ancestor of fig trees was a freestanding tree and active pollination is inferred as the ancestral state, contrary to previous hypotheses. However, ambiguity remains on the ancestral breeding system. Despite morphological plasticity, we advocate restoring a central role to morphology in our understanding of the evolution of Ficus, as it can help detect systematic errors that appear more pronounced with larger molecular datasets.

Segmentation of the subcuticular fat body in Apis mellifera females with different reproductive potentials

Evolution has created different castes of females in eusocial haplodiploids. The difference between them lies in their functions and vulnerability but above all in their reproductive potentials. Honeybee queens are highly fertile. On the other hand, the workers are facultatively sterile. However, rebel workers, i.e. workers that develop in a queenless colony, reproduce more often than normal workers. As a result, the fat body of these bees, which apart from acting as the energy reserve, is also the site of numerous metabolic processes, had to specialize in different functions perfected over millions of years of eusocial evolution. Assuming that the variety of functions manifests itself in the pleomorphic structure of the fat body cells, we predicted that also different parts of the fat body, e.g. from different segments of the abdomen, contain different sets of cells. Such differences could be expected between queens, rebels and normal workers, i.e. females with dramatically different reproductive potentials. We confirmed all these expectations. Although all bees had the same types of cells, their proportion and segmental character corresponded with the caste reproductive potential and physiological characteristics shaped in the evolutionary process. The females with an increased reproductive potential were characterized by the presence of oenocytes in the third tergite and high concentrations of compounds responsible for energy reserves, like glucose, glycogen and triglycerides. Queens had very large trophocytes, especially in the third tergite. Only in workers did we observe intercellular spaces in all the segments of the fat body, as well as high protein concentrations-especially in the sternite. As expected, the rebels combined many features of the queens and normal workers, what with other findings can help understand the ways that led to the origin of different castes in females of eusocial Hymenoptera.

Stingless bees and microbial interactions

Stingless bees (Meliponini) are a monophyletic group of eusocial insects inhabiting tropical and subtropical regions. These insects represent the most abundant and diversified group of corbiculate bees. Meliponini mostly rely on fermentation by symbiont microbes to preserve honey and transform pollen in stored food. Bee nests harbor diverse microbiota that includes bacteria, yeasts, filamentous fungi, and viruses. These microorganisms may interact with the bees through symbiotic relationships, or they may act as food for the insects, or produce biomolecules that aid in the biotransformation of bee products, such as honey and bee bread. Certain microbial species can also produce antimicrobial compounds that inhibit opportunistic bee pathogens.

Investigating Morphological Complexes Using Informational Dissonance and Bayes Factors: A Case Study in Corbiculate Bees

It is widely recognized that different regions of a genome often have different evolutionary histories and that ignoring this variation when estimating phylogenies can be misleading. However, the extent to which this is also true for morphological data is still largely unknown. Discordance among morphological traits might plausibly arise due to either variable convergent selection pressures or else phenomena such as hemiplasy. Here, we investigate patterns of discordance among 282 morphological characters, which we scored for 50 bee species particularly targeting corbiculate bees, a group that includes the well-known eusocial honeybees and bumblebees. As a starting point for selecting the most meaningful partitions in the data, we grouped characters as morphological modules, highly integrated trait complexes that as a result of developmental constraints or coordinated selection we expect to share an evolutionary history and trajectory. In order to assess conflict and coherence across and within these morphological modules, we used recently developed approaches for computing Bayesian phylogenetic information allied with model comparisons using Bayes factors. We found that despite considerable conflict among morphological complexes, accounting for among-character and among-partition rate variation with individual gamma distributions, rate multipliers, and linked branch lengths can lead to coherent phylogenetic inference using morphological data. We suggest that evaluating information content and dissonance among partitions is a useful step in estimating phylogenies from morphological data, just as it is with molecular data. Furthermore, we argue that adopting emerging approaches for investigating dissonance in genomic datasets may provide new insights into the integration and evolution of anatomical complexes. [Apidae; entropy; morphological modules; phenotypic integration; phylogenetic information.].

Phylogenomics of Parasitic and Nonparasitic Lice (Insecta: Psocodea): Combining Sequence Data and Exploring Compositional Bias Solutions in Next Generation Data Sets

The insect order Psocodea is a diverse lineage comprising both parasitic (Phthiraptera) and nonparasitic members (Psocoptera). The extreme age and ecological diversity of the group may be associated with major genomic changes, such as base compositional biases expected to affect phylogenetic inference. Divergent morphology between parasitic and nonparasitic members has also obscured the origins of parasitism within the order. We conducted a phylogenomic analysis on the order Psocodea utilizing both transcriptome and genome sequencing to obtain a data set of 2370 orthologous genes. All phylogenomic analyses, including both concatenated and coalescent methods suggest a single origin of parasitism within the order Psocodea, resolving conflicting results from previous studies. This phylogeny allows us to propose a stable ordinal level classification scheme that retains significant taxonomic names present in historical scientific literature and reflects the evolution of the group as a whole. A dating analysis, with internal nodes calibrated by fossil evidence, suggests an origin of parasitism that predates the K-Pg boundary. Nucleotide compositional biases are detected in third and first codon positions and result in the anomalous placement of the Amphientometae as sister to Psocomorpha when all nucleotide sites are analyzed. Likelihood-mapping and quartet sampling methods demonstrate that base compositional biases can also have an effect on quartet-based methods.[Illumina; Phthiraptera; Psocoptera; quartet sampling; recoding methods.].

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

STELAR: a statistically consistent coalescent-based species tree estimation method by maximizing triplet consistency

Background

Species tree estimation is frequently based on phylogenomic approaches that use multiple genes from throughout the genome. However, estimating a species tree from a collection of gene trees can be complicated due to the presence of gene tree incongruence resulting from incomplete lineage sorting (ILS), which is modelled by the multi-species coalescent process. Maximum likelihood and Bayesian MCMC methods can potentially result in accurate trees, but they do not scale well to large datasets.

Results

We present STELAR (Species Tree Estimation by maximizing tripLet AgReement), a new fast and highly accurate statistically consistent coalescent-based method for estimating species trees from a collection of gene trees. We formalized the constrained triplet consensus (CTC) problem and showed that the solution to the CTC problem is a statistically consistent estimate of the species tree under the multi-species coalescent (MSC) model. STELAR is an efficient dynamic programming based solution to the CTC problem which is highly accurate and scalable. We evaluated the accuracy of STELAR in comparison with SuperTriplets, which is an alternate fast and highly accurate triplet-based supertree method, and with MP-EST and ASTRAL – two of the most popular and accurate coalescent-based methods. Experimental results suggest that STELAR matches the accuracy of ASTRAL and improves on MP-EST and SuperTriplets.

Conclusions

Theoretical and empirical results (on both simulated and real biological datasets) suggest that STELAR is a valuable technique for species tree estimation from gene tree distributions.

Mating status and the evolution of eusociality: Oogenesis is independent of mating status in the solitary bee Osmia bicornis

The fundamental trait underlying eusociality is the reproductive division of labour. In honeybees (Apis mellifera), queens lay eggs while workers forage, defend and care for brood. The division of labour is maintained by pheromones including queen mandibular pheromone (QMP) produced by the queen. QMP constrains reproduction in adult honeybee workers, but in the absence of their queen workers can activate their ovaries and, although they cannot mate, they lay haploid male eggs. The reproductive ground plan hypothesis suggests that reproductive constraint may have evolved by co-opting mechanisms of reproductive control in solitary ancestors. In many insects mating is required to activate or accelerate oogenesis. Here, we use the solitary bee Osmia bicornis (Megachilidae) to test whether reproductive constraint evolved from ancestral control of reproduction by mating status. We present a structural study of the O. bicornis ovary, and compare key stages of oogenesis with honeybee workers. Importantly, we show that mating does not affect any aspect of the reproductive physiology of O. bicornis. We therefore conclude that mechanisms governing reproductive constraint in honeybees were unlikely to have been co-opted from mechanisms pertaining to mating status.

Model Choice, Missing Data, and Taxon Sampling Impact Phylogenomic Inference of Deep Basidiomycota Relationships

Resolving deep divergences in the tree of life is challenging even for analyses of genome-scale phylogenetic data sets. Relationships between Basidiomycota subphyla, the rusts and allies (Pucciniomycotina), smuts and allies (Ustilaginomycotina), and mushroom-forming fungi and allies (Agaricomycotina) were found particularly recalcitrant both to traditional multigene and genome-scale phylogenetics. Here, we address basal Basidiomycota relationships using concatenated and gene tree-based analyses of various phylogenomic data sets to examine the contribution of several potential sources of bias. We evaluate the contribution of biological causes (hard polytomy, incomplete lineage sorting) versus unmodeled evolutionary processes and factors that exacerbate their effects (e.g., fast-evolving sites and long-branch taxa) to inferences of basal Basidiomycota relationships. Bayesian Markov Chain Monte Carlo and likelihood mapping analyses reject the hard polytomy with confidence. In concatenated analyses, fast-evolving sites and oversimplified models of amino acid substitution favored the grouping of smuts with mushroom-forming fungi, often leading to maximal bootstrap support in both concatenation and coalescent analyses. On the contrary, the most conserved data subsets grouped rusts and allies with mushroom-forming fungi, although this relationship proved labile, sensitive to model choice, to different data subsets and to missing data. Excluding putative long-branch taxa, genes with high proportions of missing data and/or with strong signal failed to reveal a consistent trend toward one or the other topology, suggesting that additional sources of conflict are at play. While concatenated analyses yielded strong but conflicting support, individual gene trees mostly provided poor support for any resolution of rusts, smuts, and mushroom-forming fungi, suggesting that the true Basidiomycota tree might be in a part of tree space that is difficult to access using both concatenation and gene tree-based approaches. Inference-based assessments of absolute model fit strongly reject best-fit models for the vast majority of genes, indicating a poor fit of even the most commonly used models. While this is consistent with previous assessments of site-homogenous models of amino acid evolution, this does not appear to be the sole source of confounding signal. Our analyses suggest that topologies uniting smuts with mushroom-forming fungi can arise as a result of inappropriate modeling of amino acid sites that might be prone to systematic bias. We speculate that improved models of sequence evolution could shed more light on basal splits in the Basidiomycota, which, for now, remain unresolved despite the use of whole genome data.

Loss of developmental diapause as prerequisite for social evolution in bees

Diapause is a physiological arrest of development ahead of adverse environmental conditions and is a critical phase of the life cycle of many insects. In bees, diapause has been reported in species from all seven taxonomic families. However, they exhibit a variety of diapause strategies. These different strategies are of particular interest since shifts in the phase of the insect life cycle in which diapause occurs have been hypothesized to promote the evolution of sociality. Here we provide a comprehensive evaluation of this hypothesis with phylogenetic analysis and ancestral state reconstruction (ASR) of the ecological and evolutionary factors associated with diapause phase. We find that social lifestyle, latitude and voltinism are significant predictors of the life stage in which diapause occurs. ASR revealed that the most recent common ancestor of all bees likely exhibited developmental diapause and shifts to adult, reproductive, or no diapause have occurred in the ancestors of lineages in which social behaviour has evolved. These results provide fresh insight regarding the role of diapause as a prerequisite for the evolution of sociality in bees.

Social trait definitions influence evolutionary inferences: a phylogenetic approach to improving social terminology for bees

The comparative method relies not only on a good understanding of the phylogenetic relationships among taxa, but also on consistent terminology for describing phenotypes. Clear and consistent terminology allows similar phenotypes to be described and phylogenetically analyzed in different organisms, whereas inconsistent terminology is a major impediment to comparisons, even for taxonomically restricted groups such as bees. Here, I propose that the usefulness of social terminology can be judged by its value in phylogenetic trait-mapping aimed at uncovering evolutionary transitions between solitary and social behavior. I propose a four-step approach to evaluate and update social terminology, in which definitions are first updated based on recent behavioral studies (step 1), mapped onto a phylogeny (step 2), evaluated for their utility in the trait-mapping exercise (step 3), and then, if necessary, revised (step 4). To demonstrate the approach, I define four terms important for understanding social evolution in bees (solitary, social, eusocial, and hypersocial) and map them onto a very recent phylogeny of Apidae. This not only illustrates an objective method for evaluating social terminology, but also provides novel inferences about social evolution in Apidae, including support for a parasocial origin of eusociality and at least two Major Evolutionary Transitions to hypersociality.

Sociality emerges from solitary behaviours and reproductive plasticity in the orchid bee Euglossa dilemma

The evolution of eusociality and sterile worker castes represents a major transition in the history of life. Despite this, little is known about the mechanisms involved in the initial transition from solitary to social behaviour. It has been hypothesized that plasticity from ancestral solitary life cycles was coopted to create queen and worker castes in insect societies. Here, we tested this hypothesis by examining gene expression involved in the transition from solitary to social behaviour in the orchid bee Euglossa dilemma. To this end, we conducted observations that allowed us to classify bees into four distinct categories of solitary and social behaviour. Then, by sequencing brain and ovary transcriptomes from these behavioural phases, we identified gene expression changes overlapping with socially associated genes across multiple eusocial lineages. We find that genes involved in solitary E. dilemma ovarian plasticity overlap extensively with genes showing differential expression between fertile and sterile workers-or between queens and workers in other eusocial bees. We also find evidence that sociality in E. dilemma reflects gene expression patterns involved in solitary foraging and non-foraging nest care behaviours. Our results provide strong support for the hypothesis that eusociality emerges from plasticity found across solitary life cycles.

A comparative study of the pharyngeal plate of Apoidea (Hymenoptera: Aculeata), with implications for the understanding of phylogenetic relationships of bees

The pharyngeal plate is a morphological complex with extensive anatomical variation among bees and, therefore, potential as a source of phylogenetic information. The pharyngeal plate of bees is divided into four morphologically distinct regions: sitophore, hypopharyngeal lobe, pharyngeal rods, and median oral plate. In this work we illustrate and document in detail for the first time the pharyngeal plate of 43 bee species, providing descriptions of the morphological variation and contrasting these findings with representatives of apoid wasps (Crabronidae and Sphecidae). We evaluate and discuss the potential of this structure as a rich source of morphological information in the context of bee phylogeny and any research potentially impacted by comparative morphological data. The shape of the hypopharyngeal lobe is highly variable among suprageneric taxa of bees and can be readily employed to characterise taxa at various levels. We argue that the global patterns in the variation of the pharyngeal plate can provide information for phylogenetic inference within bees and constructed and coded 10 characters that encompass the most noticeable morphological differences discussed herein.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Evidence for Stabilizing Selection Driving Mutational Turnover of Short Motifs in the Eukaryotic Complementary Sex Determiner (Csd) Protein

Short linear motifs (SLiMs) can play pivotal functional roles in proteins, such as targeting proteins to specific subcellular localizations, modulating the efficiency of translation and tagging proteins for degradation. Until recently we had little knowledge about SLiM evolution. Only a few amino acids in these motifs are functionally important, making them likely to evolve ex nihilo and suggesting that they can play key roles in protein evolution. Several reports now suggest that these motifs can appear and disappear while their function in the protein is preserved, a process sometimes referred to as “turnover”. However, there has been a lack of specific experiments to determine whether independently evolved motifs do indeed have the same function, which would conclusively determine whether the process of turnover actually occurs. In this study, we experimentally detected evidence for such a mutational turnover process for nuclear localization signals (NLS) during the post-duplication divergence of the Complementary sex determiner (Csd) and Feminizer (Fem) proteins in the honeybee (Apis mellifera) lineage. Experiments on the nuclear transport activity of protein segments and those of the most recent common ancestor (MRCA) sequences revealed that three new NLS motifs evolved in the Csd protein during the post-duplication divergence while other NLS motifs were lost that existed before duplication. A screen for essential and newly evolved amino acids revealed that new motifs in the Csd protein evolved by one or two missense mutations coding for lysine. Amino acids that were predating the duplication were also essential in the acquisition of the C1 motif suggesting that the ex nihilo origin was constrained by preexisting amino acids in the physical proximity. Our data support a model in which stabilizing selection maintains the constancy of nuclear transport function but allowed mutational turnover of the encoding NLS motifs.

Genetic accommodation and the role of ancestral plasticity in the evolution of insect eusociality

For over a century, biologists have proposed a role for phenotypic plasticity in evolution, providing an avenue for adaptation in addition to 'mutation-first' models of evolutionary change. According to the various versions of this idea, the ability of organisms to respond adaptively to their environment through phenotypic plasticity may lead to novel phenotypes that can be screened by natural selection. If these initially environmentally induced phenotypes increase fitness, then genetic accommodation can lead to allele frequency change, influencing the expression of those phenotypes. Despite the long history of 'plasticity-first' models, the importance of genetic accommodation in shaping evolutionary change has remained controversial - it is neither fully embraced nor completely discarded by most evolutionary biologists. We suggest that the lack of acceptance of genetic accommodation in some cases is related to a lack of information on its molecular mechanisms. However, recent reports of epigenetic transgenerational inheritance now provide a plausible mechanism through which genetic accommodation may act, and we review this research here. We also discuss current evidence supporting a role for genetic accommodation in the evolution of eusociality in social insects, which have long been models for studying the influence of the environment on phenotypic variation, and may be particularly good models for testing hypotheses related to genetic accommodation. Finally, we introduce 'eusocial engineering', a method by which novel social phenotypes are first induced by environmental modification and then studied mechanistically to understand how environmentally induced plasticity may lead to heritable changes in social behavior. We believe the time is right to incorporate genetic accommodation into models of the evolution of complex traits, armed with new molecular tools and a better understanding of non-genetic heritable elements.

Towards reconstructing the ancestral brain gene-network regulating caste differentiation in ants

Specialized queens and life-time unmated workers evolved once in the common ancestor of all ants, but whether caste development across ants continues to be at least partly regulated by a single core set of genes remains obscure. We analysed brain transcriptomes from five ant species (three subfamilies) and reconstructed the origins of genes with caste-biased expression. Ancient genes predating the Neoptera were more likely to regulate gyne (virgin queen) phenotypes, while caste differentiation roles of younger, ant-lineage-specific genes varied. Transcriptome profiling showed that the ancestral network for caste-specific gene-regulation has been maintained, but that signatures of common ancestry are obscured by later modifications. Adjusting for such differences, we identified a core gene-set that: 1. consistently displayed similar directions and degrees of caste-differentiated expression, and 2. have mostly not been reported as being involved in caste differentiation. These core regulatory genes exist in the genomes of ant species that secondarily lost the queen caste, but expression differences for reproductive and sterile workers are minor and similar to social paper wasps that lack differentiated castes. Many caste-biased ant genes have caste-differentiated expression in honeybees, but directions of caste bias were uncorrelated, as expected when permanent castes evolved independently in both lineages.

Combining transcriptomes and ultraconserved elements to illuminate the phylogeny of Apidae

Two increasingly popular approaches to reconstruct the Tree of Life involve whole transcriptome sequencing and the target capture of ultraconserved elements (UCEs). Both methods can be used to generate large, multigene datasets for analysis of phylogenetic relationships in non-model organisms. While targeted exon sequencing across divergent lineages is now a standard method, it is still not clear if UCE data can be readily combined with published transcriptomes. In this study, we evaluate the combination of UCEs and transcriptomes in a single analysis using genome-, transcriptome-, and UCE data for 79 bees in the largest and most biologically diverse bee family, Apidae. Using existing tools, we first developed a workflow to assemble phylogenomic data from different sources and produced two large nucleotide matrices of combined data. We then reconstructed the phylogeny of the Apidae using concatenation- and coalescent-based methods, and critically evaluated the resulting phylogenies in the context of previously published genetic, genomic, and morphological data sets. Our estimated phylogenetic trees are robustly supported and largely congruent with previous molecular hypotheses, from deep nodes to shallow species-level phylogenies. Moreover, the combined approach allows us to resolve controversial nodes of the apid Tree of Life, by clarifying the relationships among the genera of orchid bees (Euglossini) and the monophyly of the Centridini. Additionally, we present novel phylogenetic evidence supporting the monophyly of the diverse clade of cleptoparasitic Apidae and the placement of two enigmatic, oil-collecting genera (Ctenoplectra and Tetrapedia). Lastly, we propose a revised classification of the family Apidae that reflects our improved understanding of apid higher-level relationships.

Phylogenomic Evidence Overturns Current Conceptions of Social Evolution in Wasps (Vespidae)

The hypothesis that eusociality originated once in Vespidae has shaped interpretation of social evolution for decades and has driven the supposition that preimaginal morphophysiological differences between castes were absent at the outset of eusociality. Many researchers also consider casteless nest-sharing an antecedent to eusociality. Together, these ideas endorse a stepwise progression of social evolution in wasps (solitary → casteless nest-sharing → eusociality with rudimentary behavioral castes → eusociality with preimaginal caste-biasing (PCB) → morphologically differentiated castes). Here, we infer the phylogeny of Vespidae using sequence data generated via anchored hybrid enrichment from 378 loci across 136 vespid species and perform ancestral state reconstructions to test whether rudimentary and monomorphic castes characterized the initial stages of eusocial evolution. Our results reject the single origin of eusociality hypothesis, contest the supposition that eusociality emerged from a casteless nest-sharing ancestor, and suggest that eusociality in Polistinae + Vespinae began with castes having morphological differences. An abrupt appearance of castes with ontogenetically established morphophysiological differences conflicts with the current conception of stepwise social evolution and suggests that the climb up the ladder of sociality does not occur through sequential mutation. Phenotypic plasticity and standing genetic variation could explain how cooperative brood care evolved in concert with nest-sharing and how morphologically dissimilar castes arose without a rudimentary intermediate. Furthermore, PCB at the outset of eusociality implicates a subsocial route to eusociality in Polistinae + Vespinae, emphasizing the role of mother-daughter interactions and subfertility (i.e. the cost component of kin selection) in the origin of workers.

Histochemistry, immunohistochemistry and cytochemistry of the anterior midgut region of the stingless bee Melipona quadrifasciata and honey bee Apis mellifera (Hymenoptera: Apidae)

The anterior midgut region of stingless bees is anatomically differentiated with tall and narrow cells, whereas in other social and solitary bees this anatomical gut region is lacking. The objective of the present study was to describe the histochemistry, immunohistochemistry and cytochemistry of the anterior midgut region of the stingless bee Melipona quadrifasciata in comparison with the honey bee Apis mellifera. The anterior midgut region of both species was evaluated for identification of the enzymes β- galactosidase, glucose-6-phosphatase, acid phosphatase, and alkaline phosphatase, the membrane transporter aquaporin, the hormone FMRF-amide, and lysosomes. Histology of the anterior midgut region showed that this region in M. quadrifasciata workers did not present external folds of the wall, whereas the following midgut wall presented many. In A. mellifera, folds in the midgut wall occur starting from the fore- midgut transition region. Despite these morphological differences, the tests evaluated were similar in both species. β-galactosidase was not found in the anterior midgut cells. Glucose-6-phosphatase and acid phosphatase occurred in the apical region of the gut epithelium. Alkaline phosphatase occurred in vesicles in apical cytoplasm and in the basal plasma membrane infoldings of the epithelial cells. Aquaporin was found in the basal region of the midgut epithelium and in the associated visceral muscles. FMRF-amide was found only in nerve endings in the anterior midgut region. All cells in the anterior midgut region were rich in lysosomes. These results suggest that in both bee species, although they have anatomically different anterior midgut regions, these regions present high metabolic activity and function in cellular homeostasis, lipid absorption and are under neurohormone control.

Insects with similar social complexity show convergent patterns of adaptive molecular evolution

Eusociality has independently evolved multiple times in the hymenoptera, but the patterns of adaptive molecular evolution underlying the evolution and elaboration of eusociality remain uncertain. Here, we performed a population genomics study of primitively eusocial Polistes (paper wasps), and compared their patterns of molecular evolution to two social bees; Bombus (bumblebees), and Apis (honey bees). This species triad allowed us to study molecular evolution across a gradient of social complexity (Polistes < Bombus < Apis) and compare species pairs that have similar (i.e. Polistes and Bombus) or different (i.e. Polistes and Apis) life histories, while controlling for phylogenetic distance. We found that regulatory genes have high levels of positive selection in Polistes; consistent with the prediction that adaptive changes in gene regulation are important during early stages of social evolution. Polistes and Bombus exhibit greater similarity in patterns of adaptive evolution including greater overlap of genes experiencing positive selection, and greater positive selection on queen-biased genes. Our findings suggest that either adaptive evolution of a few key genes underlie the evolution of simpler forms of eusociality, or that the initial stages of social evolution lead to selection on a few key traits orchestrated by orthologous genes and networks.

The Evolutionary Dynamics of the Odorant Receptor Gene Family in Corbiculate Bees

Insects rely on chemical information to locate food, choose mates, and detect potential predators. It has been hypothesized that adaptive changes in the olfactory system facilitated the diversification of numerous insect lineages. For instance, evolutionary changes of Odorant Receptor (OR) genes often occur in parallel with modifications in life history strategies. Corbiculate bees display a diverse array of behaviors that are controlled through olfaction, including varying degrees of social organization, and manifold associations with floral resources. Here we investigated the molecular mechanisms driving the evolution of the OR gene family in corbiculate bees in comparison to other chemosensory gene families. Our results indicate that the genomic organization of the OR gene family has remained highly conserved for ∼80 Myr, despite exhibiting major changes in repertoire size among bee lineages. Moreover, the evolution of OR genes appears to be driven mostly by lineage-specific gene duplications in few genomic regions that harbor large numbers of OR genes. A selection analysis revealed that OR genes evolve under positive selection, with the strongest signals detected in recently duplicated copies. Our results indicate that chromosomal translocations had a minimal impact on OR evolution, and instead local molecular mechanisms appear to be main drivers of OR repertoire size. Our results provide empirical support to the longstanding hypothesis that positive selection shaped the diversification of the OR gene family. Together, our results shed new light on the molecular mechanisms underlying the evolution of olfaction in insects.

Phylogenomics of palearctic Formica species suggests a single origin of temporary parasitism and gives insights to the evolutionary pathway toward slave-making behaviour

Background

The ants of the Formica genus are classical model species in evolutionary biology. In particular, Darwin used Formica as model species to better understand the evolution of slave-making, a parasitic behaviour where workers of another species are stolen to exploit their workforce. In his book “On the Origin of Species” (1859), Darwin first hypothesized that slave-making behaviour in Formica evolved in incremental steps from a free-living ancestor.

Methods

The absence of a well-resolved phylogenetic tree of the genus prevent an assessment of whether relationships among Formica subgenera are compatible with this scenario. In this study, we resolve the relationships among the 4 palearctic Formica subgenera (Formica str. s., Coptoformica, Raptiformica and Serviformica) using a phylogenomic dataset of 945 genes for 16 species.

Results

We provide a reference tree resolving the relationships among the main Formica subgenera with high bootstrap supports.

Discussion

The branching order of our tree suggests that the free-living lifestyle is ancestral in the Formica genus and that parasitic colony founding could have evolved a single time, probably acting as a pre-adaptation to slave-making behaviour.

Conclusion

This phylogenetic tree provides a solid backbone for future evolutionary studies in the Formica genus and slave-making behaviour.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1159-4) contains supplementary material, which is available to authorized users.

The Nuclear and Mitochondrial Genomes of the Facultatively Eusocial Orchid Bee Euglossa dilemma

Bees provide indispensable pollination services to both agricultural crops and wild plant populations, and several species of bees have become important models for the study of learning and memory, plant–insect interactions, and social behavior. Orchid bees (Apidae: Euglossini) are especially important to the fields of pollination ecology, evolution, and species conservation. Here we report the nuclear and mitochondrial genome sequences of the orchid bee Euglossa dilemma Bembé & Eltz. E. dilemma was selected because it is widely distributed, highly abundant, and it was recently naturalized in the southeastern United States. We provide a high-quality assembly of the 3.3 Gb genome, and an official gene set of 15,904 gene annotations. We find high conservation of gene synteny with the honey bee throughout 80 MY of divergence time. This genomic resource represents the first draft genome of the orchid bee genus Euglossa, and the first draft orchid bee mitochondrial genome, thus representing a valuable resource to the research community.

Queens and Workers Contribute Differently to Adaptive Evolution in Bumble Bees and Honey Bees

Eusociality represents a major transition in evolution and is typified by cooperative brood care and reproductive division of labor between generations. In bees, this division of labor allows queens and workers to phenotypically specialize. Worker traits associated with helping are thought to be crucial to the fitness of a eusocial lineage, and recent studies of honey bees (genus Apis) have found that adaptively evolving genes often have worker-biased expression patterns. It is unclear however if worker-biased genes are disproportionately acted on by strong positive selection in all eusocial insects. We undertook a comparative population genomics study of bumble bees (Bombus) and honey bees to quantify natural selection on queen- and worker-biased genes across two levels of social complexity. Despite sharing a common eusocial ancestor, genes, and gene groups with the highest levels of positive selection were often unique within each genus, indicating that life history and the environment, but not sociality per se, drives patterns of adaptive molecular evolution. We uncovered differences in the contribution of queen- and worker-biased genes to adaptive evolution in bumble bees versus honey bees. Unlike honey bees, where worker-biased genes are enriched for signs of adaptive evolution, genes experiencing positive selection in bumble bees were predominately expressed by reproductive foundresses during the initial solitary-founding stage of colonies. Our study suggests that solitary founding is a major selective pressure and that the loss of queen totipotency may cause a change in the architecture of selective pressures upon the social insect genome.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as "GC biased gene conversion". As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

Dynamic microbiome evolution in social bees

The highly social (eusocial) corbiculate bees, comprising the honey bees, bumble bees, and stingless bees, are ubiquitous insect pollinators that fulfill critical roles in ecosystem services and human agriculture. Here, we conduct wide sampling across the phylogeny of these corbiculate bees and reveal a dynamic evolutionary history behind their microbiota, marked by multiple gains and losses of gut associates, the presence of generalist as well as host-specific strains, and patterns of diversification driven, in part, by host ecology (for example, colony size). Across four continents, we found that different host species have distinct gut communities, largely independent of geography or sympatry. Nonetheless, their microbiota has a shared heritage: The emergence of the eusocial corbiculate bees from solitary ancestors appears to coincide with the acquisition of five core gut bacterial lineages, supporting the hypothesis that host sociality facilitates the development and maintenance of specialized microbiomes.

Analytical Biases Associated with GC-Content in Molecular Evolution

Molecular evolution is being revolutionized by high-throughput sequencing allowing an increased amount of genome-wide data available for multiple species. While base composition summarized by GC-content is one of the first metrics measured in genomes, its genomic distribution is a frequently neglected feature in downstream analyses based on DNA sequence comparisons. Here, we show how base composition heterogeneity among loci and taxa can bias common molecular evolution analyses such as phylogenetic tree reconstruction, detection of natural selection and estimation of codon usage. We then discuss the biological, technical and methodological causes of these GC-associated biases and suggest approaches to overcome them.

Foraging in an unsteady world: bumblebee flight performance in field-realistic turbulence

Natural environments are characterized by variable wind that can pose significant challenges for flying animals and robots. However, our understanding of the flow conditions that animals experience outdoors and how these impact flight performance remains limited. Here, we combine laboratory and field experiments to characterize wind conditions encountered by foraging bumblebees in outdoor environments and test the effects of these conditions on flight. We used radio-frequency tags to track foraging activity of uniquely identified bumblebee (Bombus impatiens) workers, while simultaneously recording local wind flows. Despite being subjected to a wide range of speeds and turbulence intensities, we find that bees do not avoid foraging in windy conditions. We then examined the impacts of turbulence on bumblebee flight in a wind tunnel. Rolling instabilities increased in turbulence, but only at higher wind speeds. Bees displayed higher mean wingbeat frequency and stroke amplitude in these conditions, as well as increased asymmetry in stroke amplitude-suggesting that bees employ an array of active responses to enable flight in turbulence, which may increase the energetic cost of flight. Our results provide the first direct evidence that moderate, environmentally relevant turbulence affects insect flight performance, and suggest that flying insects use diverse mechanisms to cope with these instabilities.

A Genome-Scale Investigation of How Sequence, Function, and Tree-Based Gene Properties Influence Phylogenetic Inference

Molecular phylogenetic inference is inherently dependent on choices in both methodology and data. Many insightful studies have shown how choices in methodology, such as the model of sequence evolution or optimality criterion used, can strongly influence inference. In contrast, much less is known about the impact of choices in the properties of the data, typically genes, on phylogenetic inference. We investigated the relationships between 52 gene properties (24 sequence-based, 19 function-based, and 9 tree-based) with each other and with three measures of phylogenetic signal in two assembled data sets of 2,832 yeast and 2,002 mammalian genes. We found that most gene properties, such as evolutionary rate (measured through the percent average of pairwise identity across taxa) and total tree length, were highly correlated with each other. Similarly, several gene properties, such as gene alignment length, Guanine-Cytosine content, and the proportion of tree distance on internal branches divided by relative composition variability (treeness/RCV), were strongly correlated with phylogenetic signal. Analysis of partial correlations between gene properties and phylogenetic signal in which gene evolutionary rate and alignment length were simultaneously controlled, showed similar patterns of correlations, albeit weaker in strength. Examination of the relative importance of each gene property on phylogenetic signal identified gene alignment length, alongside with number of parsimony-informative sites and variable sites, as the most important predictors. Interestingly, the subsets of gene properties that optimally predicted phylogenetic signal differed considerably across our three phylogenetic measures and two data sets; however, gene alignment length and RCV were consistently included as predictors of all three phylogenetic measures in both yeasts and mammals. These results suggest that a handful of sequence-based gene properties are reliable predictors of phylogenetic signal and could be useful in guiding the choice of phylogenetic markers.