The draft genome of a socially polymorphic halictid bee, Lasioglossum albipes

Identification of chemosensory genes in the stingless bee Tetragonisca fiebrigi

Reception of chemical information from the environment is crucial for insects' survival and reproduction. The chemosensory reception mainly occurs by the antennae and mouth parts of the insect, when the stimulus contacts the chemoreceptors located within the sensilla. Chemosensory receptor genes have been well-studied in some social hymenopterans such as ants, honeybees, and wasps. However, although stingless bees are the most representative group of eusocial bees, little is known about their odorant, gustatory, and ionotropic receptor genes. Here, we analyze the transcriptome of the proboscis and antennae of the stingless bee Tetragonisca fiebrigi. We identified and annotated 9 gustatory and 15 ionotropic receptors. Regarding the odorant receptors, we identified 204, and we were able to annotate 161 of them. In addition, we compared the chemosensory receptor genes of T. fiebrigi with those annotated for other species of Hymenoptera. We found that T. fiebrigi showed the largest number of odorant receptors compared with other bees. Genetic expansions were identified in the subfamilies 9-exon, which was also expanded in ants and paper wasps; in G02A, including receptors potentially mediating social behavior; and in GUnC, which has been related to pollen and nectar scent detection. Our study provides the first report of chemosensory receptor genes in T. fiebrigi and represents a resource for future molecular and physiological research in this and other stingless bee species.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

The evolutionary history of bees in time and space

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

The gut microbiome of solitary bees is mainly affected by pathogen assemblage and partially by land use

Pollinators, including solitary bees, are drastically declining worldwide. Among the factors contributing to this decline, bee pathogens and different land uses are of relevance. The link between the gut microbiome composition and host health has been recently studied for social pollinators (e.g. honeybees), whereas the information related to solitary bees is sparse. This work aimed at the characterization of the gut microbiome of the solitary bees Xylocopa augusti, Eucera fervens and Lasioglossum and attempted to correlate the gut microbial composition with the presence and load of different pathogens and land uses. Solitary bees were sampled in different sites (i.e. a farm, a natural reserve, and an urban plant nursery) showing different land uses. DNA was extracted from the gut, 16S rRNA gene amplified and sequenced. Eight pathogens, known for spillover from managed bees to wild ones, were quantified with qPCR. The results showed that the core microbiome profile of the three solitary bees significantly varied in the different species. Pseudomonas was found as the major core taxa in all solitary bees analyzed, whereas Lactobacillus, Spiroplasma and Sodalis were the second most abundant taxa in X. augusti, E. fervens and Lasioglossum, respectively. The main pathogens detected with qPCR were Nosema ceranae, Nosema bombi and Crithidia bombi, although differently abundant in the different bee species and sampling sites. Most microbial taxa did not show any correlation with the land use, apart from Snodgrassella and Nocardioides, showing higher abundances on less anthropized sites. Conversely, the pathogens species and load strongly affected the gut microbial composition, with Bifidobacterium, Apibacter, Serratia, Snodgrassella and Sodalis abundance that positively or negatively correlated with the detected pathogens load. Therefore, pathogens presence and load appear to be the main factor shaping the gut microbiome of solitary bees in Argentina.

Testing the systematic status of Homalictus and Rostrohalictus with weakened cross-vein groups within Halictini (Hymenoptera: Halictidae) using low-coverage whole-genome sequencing

The halictid genus Lasioglossum, as one of the most species-rich bee groups with persistently contentious subgeneric boundaries, is one of the most challenging bee groups from a systematic standpoint. An enduring question is the relationship of Lasioglossum and Homalictus, whether all halictine bees with weakened distal wing venation comprise one or multiple genera. Here, we analyzed the phylogenetic relationships among the subgroups within Lasioglossum s.l. based on thousands of single-copy orthologs and ultraconserved elements, which were extracted from 23 newly sequenced low-coverage whole genomes alongside a published genome (22 ingroups plus 2 outgroups). Both marker sets provided consistent results across maximum likelihood and coalescent-based species tree approaches. The phylogenetic and topology test results show that the Lasioglossum and Hemihalictus series are reciprocally monophyletic and Homalictus and Rostrohalictus are valid subgenera of Lasioglossum. Consequently, we lower Homalictus to subgenus status within Lasioglossum again, and we also raise Rostrohalictus to subgenus status from its prior synonymy with subgenus Hemihalictus. Lasioglossum przewalskyi is also transferred to the subgenus Hemihalictus. Ultimately, we redefine Lasioglossum to include all halictine bees with weakened distal wing venation.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.

Wild Bee Nutritional Ecology: Integrative Strategies to Assess Foraging Preferences and Nutritional Requirements

Bees depend on flowering plants for their nutrition, and reduced availability of floral resources is a major driver of declines in both managed and wild bee populations. Understanding the nutritional needs of different bee species, and how these needs are met by the varying nutritional resources provided by different flowering plant taxa, can greatly inform land management recommendations to support bee populations and their associated ecosystem services. However, most bee nutrition research has focused on the three most commonly managed and commercially reared bee taxa—honey bees, bumble bees, and mason bees—with fewer studies focused on wild bees and other managed species, such as leafcutting bees, stingless bees, and alkali bees. Thus, we have limited information about the nutritional requirements and foraging preferences of the vast majority of bee species. Here, we discuss the approaches traditionally used to understand bee nutritional ecology: identification of floral visitors of selected focal plant species, evaluation of the foraging preferences of adults in selected focal bee species, evaluation of the nutritional requirements of focal bee species (larvae or adults) in controlled settings, and examine how these methods may be adapted to study a wider range of bee species. We also highlight emerging technologies that have the potential to greatly facilitate studies of the nutritional ecology of wild bee species, as well as evaluate bee nutritional ecology at significantly larger spatio-temporal scales than were previously feasible. While the focus of this review is on bee species, many of these techniques can be applied to other pollinator taxa as well.

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.

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Gene Tree Estimation Error with Ultraconserved Elements: An Empirical Study on Pseudapis Bees

Summarizing individual gene trees to species phylogenies using two-step coalescent methods is now a standard strategy in the field of phylogenomics. However, practical implementations of summary methods suffer from gene tree estimation error, which is caused by various biological and analytical factors. Greatly understudied is the choice of gene tree inference method and downstream effects on species tree estimation for empirical data sets. To better understand the impact of this method choice on gene and species tree accuracy, we compare gene trees estimated through four widely used programs under different model-selection criteria: PhyloBayes, MrBayes, IQ-Tree, and RAxML. We study their performance in the phylogenomic framework of $>$800 ultraconserved elements from the bee subfamily Nomiinae (Halictidae). Our taxon sampling focuses on the genus Pseudapis, a distinct lineage with diverse morphological features, but contentious morphology-based taxonomic classifications and no molecular phylogenetic guidance. We approximate topological accuracy of gene trees by assessing their ability to recover two uncontroversial, monophyletic groups, and compare branch lengths of individual trees using the stemminess metric (the relative length of internal branches). We further examine different strategies of removing uninformative loci and the collapsing of weakly supported nodes into polytomies. We then summarize gene trees with ASTRAL and compare resulting species phylogenies, including comparisons to concatenation-based estimates. Gene trees obtained with the reversible jump model search in MrBayes were most concordant on average and all Bayesian methods yielded gene trees with better stemminess values. The only gene tree estimation approach whose ASTRAL summary trees consistently produced the most likely correct topology, however, was IQ-Tree with automated model designation (ModelFinder program). We discuss these findings and provide practical advice on gene tree estimation for summary methods. Lastly, we establish the first phylogeny-informed classification for Pseudapis s. l. and map the distribution of distinct morphological features of the group. [ASTRAL; Bees; concordance; gene tree estimation error; IQ-Tree; MrBayes, Nomiinae; PhyloBayes; RAxML; phylogenomics; stemminess].

Local Variation in Recombination Rates of the Honey Bee (Apis mellifera) Genome among Samples from Six Disparate Populations

Meiotic recombination is an essential component of eukaryotic sexual reproduction but its frequency varies within and between genomes. Although it is well-established that honey bees have a high recombination rate with about 20 cM/Mbp, the proximate and ultimate causes of this exceptional rate are poorly understood. Here, we describe six linkage maps of the Western Honey Bee Apis mellifera that were produced with consistent methodology from samples from distinct parts of the species' near global distribution. We compared the genome-wide rates and distribution of meiotic crossovers among the six maps and found considerable differences. Overall similarity of local recombination rates among our samples was unrelated to geographic or phylogenetic distance of the populations that our samples were derived from. However, the limited sampling constrains the interpretation of our results because it is unclear how representative these samples are. In contrast to previous studies, we found only in two datasets a significant relation between local recombination rate and GC content. Focusing on regions of particularly increased or decreased recombination in specific maps, we identified several enriched gene ontologies in these regions and speculate about their local adaptive relevance. These data are contributing to an increasing comparative effort to gain an understanding of the intra-specific variability of recombination rates and their evolutionary role in honey bees and other social insects.

Metabolomics-based biomarker discovery for bee health monitoring: A proof of concept study concerning nutritional stress in Bombus terrestris

Bee pollinators are exposed to multiple natural and anthropogenic stressors. Understanding the effects of a single stressor in the complex environmental context of antagonistic/synergistic interactions is critical to pollinator monitoring and may serve as early warning system before a pollination crisis. This study aimed to methodically improve the diagnosis of bee stressors using a simultaneous untargeted and targeted metabolomics-based approach. Analysis of 84 Bombus terrestris hemolymph samples found 8 metabolites retained as potential biomarkers that showed excellent discrimination for nutritional stress. In parallel, 8 significantly altered metabolites, as revealed by targeted profiling, were also assigned as candidate biomarkers. Furthermore, machine learning algorithms were applied to the above-described two biomarker sets, whereby the untargeted eight components showed the best classification performance with sensitivity and specificity up to 99% and 100%, respectively. Based on pathway and biochemistry analysis, we propose that gluconeogenesis contributed significantly to blood sugar stability in bumblebees maintained on a low carbohydrate diet. Taken together, this study demonstrates that metabolomics-based biomarker discovery holds promising potential for improving bee health monitoring and to identify stressor related to energy intake and other environmental stressors.

Population Genomics Insights into Adaptive Evolution and Ecological Differentiation in Streptomycetes

Deciphering the genomic variation that represents microevolutionary processes toward species divergence is key to understanding microbial speciation, which has long been under debate. Streptomycetes are filamentous bacteria that are ubiquitous in nature and the richest source of antibiotics; however, their speciation processes remain unknown. To tackle this issue, we performed a comprehensive population genomics analysis on Streptomyces albidoflavus residing in different habitats and with a worldwide distribution and identified and characterized the foundational changes within the species. We detected three well-defined phylogenomic clades, of which clades I and III mainly contained free-living (soil/marine) and insect-associated strains, respectively, and clade II had a mixed origin. By performing genome-wide association studies (GWAS), we identified a number of genetic variants associated with free-living or entomic (denoting or relating to insects) habitats in both the accessory and core genomes. These variants contributed collectively to the population structure and had annotated or confirmed functions that likely facilitate differential adaptation of the species. In addition, we detected higher levels of homologous recombination within each clade and in the free-living group than within the whole species and in the entomic group. A subset of the insect-associated strains (clade III) showed a relatively independent evolutionary trajectory with more symbiosis-favorable genes but little genetic interchange with the other lineages. Our results demonstrate that ecological adaptation promotes genetic differentiation in S. albidoflavus, suggesting a model of ecological speciation with gene flow in streptomycetes.IMPORTANCE Species are the fundamental units of ecology and evolution, and speciation leads to the astounding diversity of life on Earth. Studying speciation is thus of great significance to understand, protect, and exploit biodiversity, but it is a challenge in the microbial world. In this study, using population genomics, we placed Streptomyces albidoflavus strains in a spectrum of speciation and showed that the genetic differences between phylogenomic clusters evolved mainly by environmental selection and gene-specific sweeps. These findings highlight the role of ecology in structuring recombining bacterial species, making a step toward a deeper understanding of microbial speciation. Our results also raise concerns of an underrated microbial diversity at the intraspecies level, which can be utilized for mining of ecologically relevant natural products.

Draft Genome Assembly and Population Genetics of an Agricultural Pollinator, the Solitary Alkali Bee (Halictidae: Nomia melanderi)

Alkali bees (Nomia melanderi) are solitary relatives of the halictine bees, which have become an important model for the evolution of social behavior, but for which few solitary comparisons exist. These ground-nesting bees defend their developing offspring against pathogens and predators, and thus exhibit some of the key traits that preceded insect sociality. Alkali bees are also efficient native pollinators of alfalfa seed, which is a crop of major economic value in the United States. We sequenced, assembled, and annotated a high-quality draft genome of 299.6 Mbp for this species. Repetitive content makes up more than one-third of this genome, and previously uncharacterized transposable elements are the most abundant type of repetitive DNA. We predicted 10,847 protein coding genes, and identify 479 of these undergoing positive directional selection with the use of population genetic analysis based on low-coverage whole genome sequencing of 19 individuals. We found evidence of recent population bottlenecks, but no significant evidence of population structure. We also identify 45 genes enriched for protein translation and folding, transcriptional regulation, and triglyceride metabolism evolving slower in alkali bees compared to other halictid bees. These resources will be useful for future studies of bee comparative genomics and pollinator health research.

The red bayberry genome and genetic basis of sex determination

Morella rubra, red bayberry, is an economically important fruit tree in south China. Here, we assembled the first high-quality genome for both a female and a male individual of red bayberry. The genome size was 313-Mb, and 90% sequences were assembled into eight pseudo chromosome molecules, with 32 493 predicted genes. By whole-genome comparison between the female and male and association analysis with sequences of bulked and individual DNA samples from female and male, a 59-Kb region determining female was identified and located on distal end of pseudochromosome 8, which contains abundant transposable element and seven putative genes, four of them are related to sex floral development. This 59-Kb female-specific region was likely to be derived from duplication and rearrangement of paralogous genes and retained non-recombinant in the female-specific region. Sex-specific molecular markers developed from candidate genes co-segregated with sex in a genetically diverse female and male germplasm. We propose sex determination follow the ZW model of female heterogamety. The genome sequence of red bayberry provides a valuable resource for plant sex chromosome evolution and also provides important insights for molecular biology, genetics and modern breeding in Myricaceae family.

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.

The genetic basis of a social polymorphism in halictid bees

The emergence of eusociality represents a major evolutionary transition from solitary to group reproduction. The most commonly studied eusocial species, honey bees and ants, represent the behavioral extremes of social evolution but lack close relatives that are non-social. Unlike these species, the halictid bee Lasioglossum albipes produces both solitary and eusocial nests and this intraspecific variation has a genetic basis. Here, we identify genetic variants associated with this polymorphism, including one located in the intron of syntaxin 1a (syx1a), a gene that mediates synaptic vesicle release. We show that this variant can alter gene expression in a pattern consistent with differences between social and solitary bees. Surprisingly, syx1a and several other genes associated with sociality in L. albipes have also been implicated in autism spectrum disorder in humans. Thus, genes underlying behavioral variation in L. albipes may also shape social behaviors across a wide range of taxa, including humans.

The halictid bee Lasioglossum albipes has both solitary and eusocial individuals, making it a model for social evolution. Here, Kocher et al. identify a genetic variation associated with this social polymorphism, including a variant that can regulate the expression of an autism-associated gene, syntaxin 1a.

Back to the roots: the importance of using simple insect societies to understand the molecular basis of complex social life

The evolutionary trajectories toward insect eusociality come in two broad forms. In species like wasps, bees, and ants, the first helpers remained at the nest primarily to help with brood care. In species like aphids and termites, on the other hand, nest defense was initially the primary ecological driving force. To better understand the molecular basis of these two alternative evolutionary trajectories, it is therefore important to study the mechanistic basis of brood care and nest defense behavior. So far, most studies have compared morphologically distinct castes in advanced eusocial species of ants, bees, wasps, and termites. However, the interpretation of such comparisons is limited by multiple confounding factors and the fact that castes are typically fixed and cannot be manipulated at the adult stage. In this review, we argue that conducting molecular studies of brood care and nest defense in simpler, more flexible insect societies may complement studies of advanced eusocial insects and provide avenues toward more functional analyses. We review the available literature and propose candidate study systems for future molecular investigations of brood care and nest defense in social insects.

Positive selection on sociobiological traits in invasive fire ants

The fire ant Solenopsis invicta and its close relatives are highly invasive. Enhanced social cooperation may facilitate invasiveness in these and other invasive ant species. We investigated whether invasiveness in Solenopsis fire ants was accompanied by positive selection on sociobiological traits by applying a phylogenomics approach to infer ancient selection, and a population genomics approach to infer recent and ongoing selection in both native and introduced S. invicta populations. A combination of whole-genome sequencing of 40 haploid males and reduced-representation genomic sequencing of 112 diploid workers identified 1,758,116 and 169,682 polymorphic markers, respectively. The resulting high-resolution maps of genomic polymorphism provide high inference power to test for positive selection. Our analyses provide evidence of positive selection on putative ion channel genes, which are implicated in neurological functions, and on vitellogenin, which is a key regulator of development and caste determination. Furthermore, molecular functions implicated in pheromonal signalling have experienced recent positive selection. Genes with signatures of positive selection were significantly more often those overexpressed in workers compared with queens and males, suggesting that worker traits are under stronger selection than queen and male traits. These results provide insights into selection pressures and ongoing adaptation in an invasive social insect and support the hypothesis that sociobiological traits are under more positive selection than nonsocial traits in such invasive species.

Variation in DNA Methylation Is Not Consistently Reflected by Sociality in Hymenoptera

Changes in gene regulation that underlie phenotypic evolution can be encoded directly in the DNA sequence or mediated by chromatin modifications such as DNA methylation. It has been hypothesized that the evolution of eusocial division of labor is associated with enhanced gene regulatory potential, which may include expansions in DNA methylation in the genomes of Hymenoptera (bees, ants, wasps, and sawflies). Recently, this hypothesis garnered support from analyses of a commonly used metric to estimate DNA methylation in silico, CpG content. Here, we test this hypothesis using direct, nucleotide-level measures of DNA methylation across nine species of Hymenoptera. In doing so, we generated new DNA methylomes for three species of interest, including one solitary and one facultatively eusocial halictid bee and a sawfly. We demonstrate that the strength of correlation between CpG content and DNA methylation varies widely among hymenopteran taxa, highlighting shortcomings in the utility of CpG content as a proxy for DNA methylation in comparative studies of taxa with sparse DNA methylomes. We observed strikingly high levels of DNA methylation in the sawfly relative to other investigated hymenopterans. Analyses of molecular evolution suggest the relatively distinct sawfly DNA methylome may be associated with positive selection on functional DNMT3 domains. Sawflies are an outgroup to all ants, bees, and wasps, and no sawfly species are eusocial. We find no evidence that either global expansions or variation within individual ortholog groups in DNA methylation are consistently associated with the evolution of social behavior.

Notos - a galaxy tool to analyze CpN observed expected ratios for inferring DNA methylation types

BACKGROUND

DNA methylation patterns store epigenetic information in the vast majority of eukaryotic species. The relatively high costs and technical challenges associated with the detection of DNA methylation however have created a bias in the number of methylation studies towards model organisms. Consequently, it remains challenging to infer kingdom-wide general rules about the functions and evolutionary conservation of DNA methylation. Methylated cytosine is often found in specific CpN dinucleotides, and the frequency distributions of, for instance, CpG observed/expected (CpG o/e) ratios have been used to infer DNA methylation types based on higher mutability of methylated CpG.

RESULTS

Predominantly model-based approaches essentially founded on mixtures of Gaussian distributions are currently used to investigate questions related to the number and position of modes of CpG o/e ratios. These approaches require the selection of an appropriate criterion for determining the best model and will fail if empirical distributions are complex or even merely moderately skewed. We use a kernel density estimation (KDE) based technique for robust and precise characterization of complex CpN o/e distributions without a priori assumptions about the underlying distributions.

CONCLUSIONS

We show that KDE delivers robust descriptions of CpN o/e distributions. For straightforward processing, we have developed a Galaxy tool, called Notos and available at the ToolShed, that calculates these ratios of input FASTA files and fits a density to their empirical distribution. Based on the estimated density the number and shape of modes of the distribution is determined, providing a rational for the prediction of the number and the types of different methylation classes. Notos is written in R and Perl.

Limited social plasticity in the socially polymorphic sweat bee Lasioglossum calceatum

Abstract

Eusociality is characterised by a reproductive division of labour, where some individuals forgo direct reproduction to instead help raise kin. Socially polymorphic sweat bees are ideal models for addressing the mechanisms underlying the transition from solitary living to eusociality, because different individuals in the same species can express either eusocial or solitary behaviour. A key question is whether alternative social phenotypes represent environmentally induced plasticity or predominantly genetic differentiation between populations. In this paper, we focus on the sweat bee Lasioglossum calceatum, in which northern or high-altitude populations are solitary, whereas more southern or low-altitude populations are typically eusocial. To test whether social phenotype responds to local environmental cues, we transplanted adult females from a solitary, northern population, to a southern site where native bees are typically eusocial. Nearly all native nests were eusocial, with foundresses producing small first brood (B1) females that became workers. In contrast, nine out of ten nests initiated by transplanted bees were solitary, producing female offspring that were the same size as the foundress and entered directly into hibernation. Only one of these ten nests became eusocial. Social phenotype was unlikely to be related to temperature experienced by nest foundresses when provisioning B1 offspring, or by B1 emergence time, both previously implicated in social plasticity seen in two other socially polymorphic sweat bees. Our results suggest that social polymorphism in L. calceatum predominantly reflects genetic differentiation between populations, and that plasticity is in the process of being lost by bees in northern populations.

Significance statement

Phenotypic plasticity is thought to play a key role in the early stages of the transition from solitary to eusocial behaviour, but may then be lost if environmental conditions become less variable. Socially polymorphic sweat bees exhibit either solitary or eusocial behaviour in different geographic populations, depending on the length of the nesting season. We tested for plasticity in the socially polymorphic sweat bee Lasioglossum calceatum by transplanting nest foundresses from a northern, non-eusocial population to a southern, eusocial population. Plasticity would be detected if transplanted bees exhibited eusocial behaviour. We found that while native bees were eusocial, 90% of transplanted bees and their offspring did not exhibit traits associated with eusociality. Environmental variables such as time of offspring emergence or temperatures experienced by foundresses during provisioning could not explain these differences. Our results suggest that the ability of transplanted bees to express eusociality is being lost, and that social polymorphism predominantly reflects genetic differences between populations.

Electronic supplementary material

The online version of this article (10.1007/s00265-018-2475-9) contains supplementary material, which is available to authorized users.

Genomes of the Hymenoptera

Hymenoptera is the second-most sequenced arthropod order, with 52 publically archived genomes (71 with ants, reviewed elsewhere), however these genomes do not capture the breadth of this very diverse order (Figure 1, Table 1). These sequenced genomes represent only 15 of the 97 extant families. Although at least 55 other genomes are in progress in an additional 11 families (see Table 2), stinging wasps represent 35 (67%) of the available and 42 (76%) of the in progress genomes. A more comprehensive catalog of hymenopteran genomes is needed for research into the evolutionary processes underlying the expansive diversity in terms of ecology, behavior, and physiological traits within this group. Additional sequencing is needed to generate an assembly for even 0.05% of the estimated 1 million hymenopteran species, and we recommend premier level assemblies for at least 0.1% of the >150,000 named species dispersed across the order. Given the haplodiploid sex determination in Hymenoptera, haploid male sequencing will help minimize genome assembly issues to enable higher quality genome assemblies.

Building a new research framework for social evolution: intralocus caste antagonism

The breeding and non‐breeding ‘castes’ of eusocial insects provide a striking example of role‐specific selection, where each caste maximises fitness through different morphological, behavioural and physiological trait values. Typically, queens are long‐lived egg‐layers, while workers are short‐lived, largely sterile foragers. Remarkably, the two castes are nevertheless produced by the same genome. The existence of inter‐caste genetic correlations is a neglected consequence of this shared genome, potentially hindering the evolution of caste dimorphism: alleles that increase the productivity of queens may decrease the productivity of workers and vice versa, such that each caste is prevented from reaching optimal trait values. A likely consequence of this ‘intralocus caste antagonism’ should be the maintenance of genetic variation for fitness and maladaptation within castes (termed ‘caste load’), analogous to the result of intralocus sexual antagonism. The aim of this review is to create a research framework for understanding caste antagonism, drawing in part upon conceptual similarities with sexual antagonism. By reviewing both the social insect and sexual antagonism literature, we highlight the current empirical evidence for caste antagonism, discuss social systems of interest, how antagonism might be resolved, and challenges for future research. We also introduce the idea that sexual and caste antagonism could interact, creating a three‐way antagonism over gene expression. This includes unpacking the implications of haplodiploidy for the outcome of this complex interaction.

Hymenoptera Genome Database: Using HymenopteraMine to Enhance Genomic Studies of Hymenopteran Insects

The Hymenoptera Genome Database (HGD; http://hymenopteragenome.org ) is a genome informatics resource for insects of the order Hymenoptera, which includes bees, ants and wasps. HGD provides genome browsers with manual annotation tools (JBrowse/Apollo), BLAST, bulk data download, and a data mining warehouse (HymenopteraMine). This chapter focuses on the use of HymenopteraMine to create annotation data sets that can be exported for use in downstream analyses. HymenopteraMine leverages the InterMine platform to combine genome assemblies and official gene sets with data from OrthoDB, RefSeq, FlyBase, Gene Ontology, UniProt, InterPro, KEGG, Reactome, dbSNP, PubMed, and BioGrid, as well as precomputed gene expression information based on publicly available RNAseq. Built-in template queries provide starting points for data exploration, while the QueryBuilder tool supports construction of complex custom queries. The List Analysis and Genomic Regions search tools execute queries based on uploaded lists of identifiers and genome coordinates, respectively. HymenopteraMine facilitates cross-species data mining based on orthology and supports meta-analyses by tracking identifiers across gene sets and genome assemblies.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

Nutritional Regulation of Phenotypic Plasticity in a Solitary Bee (Hymenoptera: Megachilidae)

Phenotypic plasticity involves adaptive responses to predictable environmental fluctuations and may promote evolutionary change. We studied the regulation of phenotypic plasticity in an important agricultural pollinator, the solitary alfalfa leafcutting bee (Megachile rotundata F.). Specifically, we investigated how larval nutrition affects M. rotundata diapause plasticity and how diapause plasticity affects adult female reproductive behavior. Field surveys and laboratory manipulations of aspects of larval diet demonstrated nutritional regulation of M. rotundata diapause plasticity. Manipulation of larval diet quality through the addition of royal jelly, the caste-determining substance of the honey bee Apis mellifera L., increased the probability of diapause in M. rotundata. We also found that larval nutrition and diapause status affected M. rotundata adult female reproductive behavior. Nutritional effects on larval diapause that also impact adult fitness have intriguing implications for the evolution of developmental plasticity in bees. In particular, as the solitary lifestyle of M. rotundata is considered to be the ancestral condition in bees, nutritionally regulated plasticity may have been an ancestral condition in all bees that facilitated the evolution of other forms of phenotypic plasticity, such as the castes of social bees.

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.

Solitary bees reduce investment in communication compared with their social relatives

Social animals must communicate to define group membership and coordinate social organization. For social insects, communication is predominantly mediated through chemical signals, and as social complexity increases, so does the requirement for a greater diversity of signals. This relationship is particularly true for advanced eusocial insects, including ants, bees, and wasps, whose chemical communication systems have been well-characterized. However, we know surprisingly little about how these communication systems evolve during the transition between solitary and group living. Here, we demonstrate that the sensory systems associated with signal perception are evolutionarily labile. In particular, we show that differences in signal production and perception are tightly associated with changes in social behavior in halictid bees. Our results suggest that social species require a greater investment in communication than their solitary counterparts and that species that have reverted from eusociality to solitary living have repeatedly reduced investment in these potentially costly sensory perception systems.

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.

Evolution of DNA Methylation across Insects

DNA methylation contributes to gene and transcriptional regulation in eukaryotes, and therefore has been hypothesized to facilitate the evolution of plastic traits such as sociality in insects. However, DNA methylation is sparsely studied in insects. Therefore, we documented patterns of DNA methylation across a wide diversity of insects. We predicted that underlying enzymatic machinery is concordant with patterns of DNA methylation. Finally, given the suggestion that DNA methylation facilitated social evolution in Hymenoptera, we tested the hypothesis that the DNA methylation system will be associated with presence/absence of sociality among other insect orders. We found DNA methylation to be widespread, detected in all orders examined except Diptera (flies). Whole genome bisulfite sequencing showed that orders differed in levels of DNA methylation. Hymenopteran (ants, bees, wasps and sawflies) had some of the lowest levels, including several potential losses. Blattodea (cockroaches and termites) show all possible patterns, including a potential loss of DNA methylation in a eusocial species whereas solitary species had the highest levels. Species with DNA methylation do not always possess the typical enzymatic machinery. We identified a gene duplication event in the maintenance DNA methyltransferase 1 (DNMT1) that is shared by some Hymenoptera, and paralogs have experienced divergent, nonneutral evolution. This diversity and nonneutral evolution of underlying machinery suggests alternative DNA methylation pathways may exist. Phylogenetically corrected comparisons revealed no evidence that supports evolutionary association between sociality and DNA methylation. Future functional studies will be required to advance our understanding of DNA methylation in insects.

Transcriptomics and neuroanatomy of the clonal raider ant implicate an expanded clade of odorant receptors in chemical communication

A major aim of sociogenomic research is to uncover common principles in the molecular evolution of sociality. This endeavor has been hampered by the small number of specific genes currently known to function in social behavior. Here we provide several lines of evidence suggesting that ants have evolved a large and novel clade of odorant receptor (OR) genes to perceive hydrocarbon-based pheromones, arguably the most important signals in ant communication. This genomic expansion is also mirrored in the ant brain via a corresponding expansion of a specific cluster of glomeruli in the antennal lobe. We show that in the clonal raider ant, hydrocarbon-sensitive basiconic sensilla are found only on the ventral surface of the female antennal club. Correspondingly, nearly all genes in a clade of 180 ORs within the 9-exon subfamily of ORs are expressed exclusively in females and are highly enriched in expression in the ventral half of the antennal club. Furthermore, we found that across species and sexes, the number of 9-exon ORs expressed in antennae is tightly correlated with the number of glomeruli in the antennal lobe region innervated by odorant receptor neurons from basiconic sensilla. Evolutionary analyses show that this clade underwent a striking gene expansion in the ancestors of all ants and slower but continued expansion in extant ant lineages. This evidence suggests that ants have evolved a large clade of genes to support pheromone perception and that gene duplications have played an important role in the molecular evolution of ant communication.

Identification of Complete Repertoire of Apis florea Odorant Receptors Reveals Complex Orthologous Relationships with Apis mellifera

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea. Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs. RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication.

Why do social insect queens live so long? Approaches to unravel the sociality-aging puzzle

Social insects are characterized by an apparent reshaping of the fecundity/longevity trade-off with sociality. Currently, we have only sketchy information about the potential underlying causes and mechanisms of aging and senescence which in addition are restricted to few model insect organisms (mainly the fruit fly Drosophila melanogaster and the honey bee Apis mellifera). How can we gain a more thorough understanding how sociality shapes senescence and the fecundity/longevity trade-off? By reviewing available literature, I propose a comparative approach that offers the opportunity to gain fundamental insights into uncovering the basis for this life history trade-off and its reshaping with sociality.

Facultative social insects can provide insights into the reversal of the longevity/fecundity trade-off across the eusocial insects

In eusocial insects, reversal of the fecundity/longevity trade-off and extreme differences in life histories between castes of the same species garner scientific and public interest. Facultative social species at the threshold of sociality, in which individuals are socially plastic, provide an excellent opportunity to understand the causes and mechanisms underlying this reversal in life history trade-off associated with eusociality. We briefly present the ultimate factors favoring sociality and the association between fecundity and longevity in facultative eusocial insects, including kin selection and disposable soma, as well as proximate mechanisms observed in such species, such as differences in hormone titers and functions. Potential genetic underpinnings of lifespan and fecundity differences between castes are discussed and future research directions are proposed.

ArthropodaCyc: a CycADS powered collection of BioCyc databases to analyse and compare metabolism of arthropods

Arthropods interact with humans at different levels with highly beneficial roles (e.g. as pollinators), as well as with a negative impact for example as vectors of human or animal diseases, or as agricultural pests. Several arthropod genomes are available at present and many others will be sequenced in the near future in the context of the i5K initiative, offering opportunities for reconstructing, modelling and comparing their metabolic networks. In-depth analysis of these genomic data through metabolism reconstruction is expected to contribute to a better understanding of the biology of arthropods, thereby allowing the development of new strategies to control harmful species. In this context, we present here ArthropodaCyc, a dedicated BioCyc collection of databases using the Cyc annotation database system (CycADS), allowing researchers to perform reliable metabolism comparisons of fully sequenced arthropods genomes. Since the annotation quality is a key factor when performing such global genome comparisons, all proteins from the genomes included in the ArthropodaCyc database were re-annotated using several annotation tools and orthology information. All functional/domain annotation results and their sources were integrated in the databases for user access. Currently, ArthropodaCyc offers a centralized repository of metabolic pathways, protein sequence domains, Gene Ontology annotations as well as evolutionary information for 28 arthropod species. Such database collection allows metabolism analysis both with integrated tools and through extraction of data in formats suitable for systems biology studies.

Database URL:http://arthropodacyc.cycadsys.org/

The Genome and Methylome of a Subsocial Small Carpenter Bee, Ceratina calcarata

Understanding the evolution of animal societies, considered to be a major transition in evolution, is a key topic in evolutionary biology. Recently, new gateways for understanding social evolution have opened up due to advances in genomics, allowing for unprecedented opportunities in studying social behavior on a molecular level. In particular, highly eusocial insect species (caste-containing societies with nonreproductives that care for siblings) have taken center stage in studies of the molecular evolution of sociality. Despite advances in genomic studies of both solitary and eusocial insects, we still lack genomic resources for early insect societies. To study the genetic basis of social traits requires comparison of genomes from a diversity of organisms ranging from solitary to complex social forms. Here we present the genome of a subsocial bee, Ceratina calcarata This study begins to address the types of genomic changes associated with the earliest origins of simple sociality using the small carpenter bee. Genes associated with lipid transport and DNA recombination have undergone positive selection in C. calcarata relative to other bee lineages. Furthermore, we provide the first methylome of a noneusocial bee. Ceratina calcarata contains the complete enzymatic toolkit for DNA methylation. As in the honey bee and many other holometabolous insects, DNA methylation is targeted to exons. The addition of this genome allows for new lines of research into the genetic and epigenetic precursors to complex social behaviors.

Genome, transcriptome and methylome sequencing of a primitively eusocial wasp reveal a greatly reduced DNA methylation system in a social insect

Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect Polistes dominula, accompanied by analysis of caste-related transcriptome and methylome sequence data for adult queens and workers. Polistes dominula possesses a fairly typical hymenopteran genome, but shows very low genomewide GC content and some evidence of reduced genome size. We found numerous caste-related differences in gene expression, with evidence that both conserved and novel genes are related to caste differences. Most strikingly, these -omics data reveal a major reduction in one of the major epigenetic mechanisms that has been previously suggested to be important for caste differences in social insects: DNA methylation. Along with a conspicuous loss of a key gene associated with environmentally responsive DNA methylation (the de novo DNA methyltransferase Dnmt3), these wasps have greatly reduced genomewide methylation to almost zero. In addition to providing a valuable resource for comparative analysis of social insect evolution, our integrative -omics data for this important behavioural and evolutionary model system call into question the general importance of DNA methylation in caste differences and evolution in social insects.

Erratum to: The draft genome of a socially polymorphic halictid bee, Lasioglossum albipes

During the type-setting of the final version of the article [1] some of the additional files were swapped, and several were completely replaced. The correct files are republished in this Erratum.

Hymenoptera Genome Database: integrating genome annotations in HymenopteraMine

We report an update of the Hymenoptera Genome Database (HGD) (http://HymenopteraGenome.org), a model organism database for insect species of the order Hymenoptera (ants, bees and wasps). HGD maintains genomic data for 9 bee species, 10 ant species and 1 wasp, including the versions of genome and annotation data sets published by the genome sequencing consortiums and those provided by NCBI. A new data-mining warehouse, HymenopteraMine, based on the InterMine data warehousing system, integrates the genome data with data from external sources and facilitates cross-species analyses based on orthology. New genome browsers and annotation tools based on JBrowse/WebApollo provide easy genome navigation, and viewing of high throughput sequence data sets and can be used for collaborative genome annotation. All of the genomes and annotation data sets are combined into a single BLAST server that allows users to select and combine sequence data sets to search.

The Genome and Methylome of a Beetle with Complex Social Behavior, Nicrophorus vespilloides (Coleoptera: Silphidae)

Testing for conserved and novel mechanisms underlying phenotypic evolution requires a diversity of genomes available for comparison spanning multiple independent lineages. For example, complex social behavior in insects has been investigated primarily with eusocial lineages, nearly all of which are Hymenoptera. If conserved genomic influences on sociality do exist, we need data from a wider range of taxa that also vary in their levels of sociality. Here, we present the assembled and annotated genome of the subsocial beetle Nicrophorus vespilloides, a species long used to investigate evolutionary questions of complex social behavior. We used this genome to address two questions. First, do aspects of life history, such as using a carcass to breed, predict overlap in gene models more strongly than phylogeny? We found that the overlap in gene models was similar between N. vespilloides and all other insect groups regardless of life history. Second, like other insects with highly developed social behavior but unlike other beetles, does N. vespilloides have DNA methylation? We found strong evidence for an active DNA methylation system. The distribution of methylation was similar to other insects with exons having the most methylated CpGs. Methylation status appears highly conserved; 85% of the methylated genes in N. vespilloides are also methylated in the hymentopteran Nasonia vitripennis. The addition of this genome adds a coleopteran resource to answer questions about the evolution and mechanistic basis of sociality and to address questions about the potential role of methylation in social behavior.

Developmental Transcriptome for a Facultatively Eusocial Bee, Megalopta genalis

Transcriptomes provide excellent foundational resources for mechanistic and evolutionary analyses of complex traits. We present a developmental transcriptome for the facultatively eusocial bee Megalopta genalis, which represents a potential transition point in the evolution of eusociality. A de novo transcriptome assembly of Megalopta genalis was generated using paired-end Illumina sequencing and the Trinity assembler. Males and females of all life stages were aligned to this transcriptome for analysis of gene expression profiles throughout development. Gene Ontology analysis indicates that stage-specific genes are involved in ion transport, cell-cell signaling, and metabolism. A number of distinct biological processes are upregulated in each life stage, and transitions between life stages involve shifts in dominant functional processes, including shifts from transcriptional regulation in embryos to metabolism in larvae, and increased lipid metabolism in adults. We expect that this transcriptome will provide a useful resource for future analyses to better understand the molecular basis of the evolution of eusociality and, more generally, phenotypic plasticity.

Chemoreceptor Evolution in Hymenoptera and Its Implications for the Evolution of Eusociality

Eusocial insects, mostly Hymenoptera, have evolved unique colonial lifestyles that rely on the perception of social context mainly through pheromones, and chemoreceptors are hypothesized to have played important adaptive roles in the evolution of sociality. However, because chemoreceptor repertoires have been characterized in few social insects and their solitary relatives, a comprehensive examination of this hypothesis has not been possible. Here, we annotate ∼3,000 odorant and gustatory receptors in recently sequenced Hymenoptera genomes and systematically compare >4,000 chemoreceptors from 13 hymenopterans, representing one solitary lineage (wasps) and three independently evolved eusocial lineages (ants and two bees). We observe a strong general tendency for chemoreceptors to expand in Hymenoptera, whereas the specifics of gene gains/losses are highly diverse between lineages. We also find more frequent positive selection on chemoreceptors in a facultative eusocial bee and in the common ancestor of ants compared with solitary wasps. Our results suggest that the frequent expansions of chemoreceptors have facilitated the transition to eusociality. Divergent expression patterns of odorant receptors between honeybee and ants further indicate differential roles of chemoreceptors in parallel trajectories of social evolution.

Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants

BACKGROUND

The use of transcriptomic and genomic datasets for phylogenetic reconstruction has become increasingly common as researchers attempt to resolve recalcitrant nodes with increasing amounts of data. The large size and complexity of these datasets introduce significant phylogenetic noise and conflict into subsequent analyses. The sources of conflict may include hybridization, incomplete lineage sorting, or horizontal gene transfer, and may vary across the phylogeny. For phylogenetic analysis, this noise and conflict has been accommodated in one of several ways: by binning gene regions into subsets to isolate consistent phylogenetic signal; by using gene-tree methods for reconstruction, where conflict is presumed to be explained by incomplete lineage sorting (ILS); or through concatenation, where noise is presumed to be the dominant source of conflict. The results provided herein emphasize that analysis of individual homologous gene regions can greatly improve our understanding of the underlying conflict within these datasets.

RESULTS

Here we examined two published transcriptomic datasets, the angiosperm group Caryophyllales and the aculeate Hymenoptera, for the presence of conflict, concordance, and gene duplications in individual homologs across the phylogeny. We found significant conflict throughout the phylogeny in both datasets and in particular along the backbone. While some nodes in each phylogeny showed patterns of conflict similar to what might be expected with ILS alone, the backbone nodes also exhibited low levels of phylogenetic signal. In addition, certain nodes, especially in the Caryophyllales, had highly elevated levels of strongly supported conflict that cannot be explained by ILS alone.

CONCLUSION

This study demonstrates that phylogenetic signal is highly variable in phylogenomic data sampled across related species and poses challenges when conducting species tree analyses on large genomic and transcriptomic datasets. Further insight into the conflict and processes underlying these complex datasets is necessary to improve and develop adequate models for sequence analysis and downstream applications. To aid this effort, we developed the open source software phyparts ( https://bitbucket.org/blackrim/phyparts ), which calculates unique, conflicting, and concordant bipartitions, maps gene duplications, and outputs summary statistics such as internode certainy (ICA) scores and node-specific counts of gene duplications.