High sugar diet alters immune function and the gut microbiome in juvenile green iguanas (Iguana iguana).

The present work aimed to study whether a high sugar diet can alter immune responses and the gut microbiome in green iguanas. Thirty-six iguanas were split into four treatment groups using a 2x2 design. Iguanas either received a sugar supplemented diet or a control diet, and either received a lipopolysaccharide (LPS) injection or a phosphate buffer solution (PBS) injection. Iguanas were given their respective diet treatment through the entire study (∼3 months) and received a primary immune challenge one month and two months into the experiment. Blood samples and cloacal swabs were taken at various points in the experiment and used to measure changes in the immune system (bacterial killing ability, lysis and agglutination scores, LPS specific IgY concentrations), and alterations in the gut microbiome. We found that sugar diet reduces bacterial killing ability following an LPS challenge, and sugar and the immune challenge temporarily alters gut microbiome composition while reducing alpha diversity. While sugar did not directly reduce lysis and agglutination following the immune challenge, the change in these scores over a 24-hour period following an immune challenge was more drastic (it decreased) relative to the control diet group. Moreover, sugar increased constitutive agglutination outside of the immune challenges (i.e., pre-challenge levels). In this study, we provide evidence that a high sugar diet affects the immune system of green iguanas (in a disruptive manner) and alters the gut microbiome.

The expansion of agriculture has shaped the recent evolutionary history of a specialized squash pollinator.

The expansion of agriculture is responsible for the mass conversion of biologically diverse natural environments into managed agroecosystems dominated by a handful of genetically homogeneous crop species. Agricultural ecosystems typically have very different abiotic and ecological conditions from those they replaced and create potential niches for those species that are able to exploit the abundant resources offered by crop plants. While there are well-studied examples of crop pests that have adapted into novel agricultural niches, the impact of agricultural intensification on the evolution of crop mutualists such as pollinators is poorly understood. We combined genealogical inference from genomic data with archaeological records to demonstrate that the Holocene demographic history of a wild specialist pollinator of Cucurbita (pumpkins, squashes, and gourds) has been profoundly impacted by the history of agricultural expansion in North America. Populations of the squash bee Eucera pruinosa experienced rapid growth in areas where agriculture intensified within the past 1,000 y, suggesting that the cultivation of Cucurbita in North America has increased the amount of floral resources available to these bees. In addition, we found that roughly 20% of this bee species' genome shows signatures of recent selective sweeps. These signatures are overwhelmingly concentrated in populations from eastern North America where squash bees were historically able to colonize novel environments due to human cultivation of Cucurbita pepo and now exclusively inhabit agricultural niches. These results suggest that the widespread cultivation of crops can prompt adaptation in wild pollinators through the distinct ecological conditions imposed by agricultural environments.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees.

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata).

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.

Age-related mushroom body expansion in male sweat bees and bumble bees.

A well-documented phenomenon among social insects is that brain changes occur prior to or at the onset of certain experiences, potentially serving to prime the brain for specific tasks. This insight comes almost exclusively from studies considering developmental maturation in females. As a result, it is unclear whether age-related brain plasticity is consistent across sexes, and to what extent developmental patterns differ. Using confocal microscopy and volumetric analyses, we investigated age-related brain changes coinciding with sexual maturation in the males of the facultatively eusocial sweat bee, Megalopta genalis, and the obligately eusocial bumble bee, Bombus impatiens. We compared volumetric measurements between newly eclosed and reproductively mature males kept isolated in the lab. We found expansion of the mushroom bodies-brain regions associated with learning and memory-with maturation, which were consistent across both species. This age-related plasticity may, therefore, play a functionally-relevant role in preparing male bees for mating, and suggests that developmentally-driven neural restructuring can occur in males, even in species where it is absent in females.

Composition and acquisition of the microbiome in solitary, ground-nesting alkali bees.

Increasing evidence suggests the microbiome plays an important role in bee ecology and health. However, the relationship between bees and their bacterial symbionts has only been explored in a handful of species. We characterized the microbiome across the life cycle of solitary, ground-nesting alkali bees (Nomia melanderi). We find that feeding status is a major determinant of microbiome composition. The microbiome of feeding larvae was similar to that of pollen provisions, but the microbiome of post-feeding larvae (pre-pupae) was similar to that of the brood cell walls and newly-emerged females. Feeding larvae and pollen provisions had the lowest beta diversity, suggesting the composition of larval diet is highly uniform. Comparisons between lab-reared, newly-emerged, and nesting adult females suggest that the hindgut bacterial community is largely shaped by the external environment. However, we also identified taxa that are likely acquired in the nest or which increase or decrease in relative abundance with age. Although Lactobacillus micheneri was highly prevalent in pollen provisions, it was only detected in one lab-reared female, suggesting it is primarily acquired from environmental sources. These results provide the foundation for future research on metagenomic function and development of probiotics for these native pollinators.

Experience, but not age, is associated with volumetric mushroom body expansion in solitary alkali bees.

In social insects, changes in behavior are often accompanied by structural changes in the brain. This neuroplasticity may come with experience (experience-dependent) or age (experience-expectant). Yet, the evolutionary relationship between neuroplasticity and sociality is unclear, because we know little about neuroplasticity in the solitary relatives of social species. We used confocal microscopy to measure brain changes in response to age and experience in a solitary halictid bee (Nomia melanderi). First, we compared the volume of individual brain regions among newly emerged females, laboratory females deprived of reproductive and foraging experience, and free-flying, nesting females. Experience, but not age, led to significant expansion of the mushroom bodies - higher-order processing centers associated with learning and memory. Next, we investigated how social experience influences neuroplasticity by comparing the brains of females kept in the laboratory either alone or paired with another female. Paired females had significantly larger olfactory regions of the mushroom bodies. Together, these experimental results indicate that experience-dependent neuroplasticity is common to both solitary and social taxa, whereas experience-expectant neuroplasticity may be an adaptation to life in a social colony. Further, neuroplasticity in response to social chemical signals may have facilitated the evolution of sociality.

Brain microRNAs among social and solitary bees.

Evolutionary transitions to a social lifestyle in insects are associated with lineage-specific changes in gene expression, but the key nodes that drive these regulatory changes are unknown. We examined the relationship between social organization and lineage-specific microRNAs (miRNAs). Genome scans across 12 bee species showed that miRNA copy-number is mostly conserved and not associated with sociality. However, deep sequencing of small RNAs in six bee species revealed a substantial proportion (20-35%) of detected miRNAs had lineage-specific expression in the brain, 24-72% of which did not have homologues in other species. Lineage-specific miRNAs disproportionately target lineage-specific genes, and have lower expression levels than shared miRNAs. The predicted targets of lineage-specific miRNAs are not enriched for genes with caste-biased expression or genes under positive selection in social species. Together, these results suggest that novel miRNAs may coevolve with novel genes, and thus contribute to lineage-specific patterns of evolution in bees, but do not appear to have significant influence on social evolution. Our analyses also support the hypothesis that many new miRNAs are purged by selection due to deleterious effects on mRNA targets, and suggest genome structure is not as influential in regulating bee miRNA evolution as has been shown for mammalian miRNAs.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

Combining Experimental Evolution and Genomics to Understand How Seed Beetles Adapt to a Marginal Host Plant.

Genes that affect adaptive traits have been identified, but our knowledge of the genetic basis of adaptation in a more general sense (across multiple traits) remains limited. We combined population-genomic analyses of evolve-and-resequence experiments, genome-wide association mapping of performance traits, and analyses of gene expression to fill this knowledge gap and shed light on the genomics of adaptation to a marginal host (lentil) by the seed beetle Callosobruchus maculatus. Using population-genomic approaches, we detected modest parallelism in allele frequency change across replicate lines during adaptation to lentil. Mapping populations derived from each lentil-adapted line revealed a polygenic basis for two host-specific performance traits (weight and development time), which had low to modest heritabilities. We found less evidence of parallelism in genotype-phenotype associations across these lines than in allele frequency changes during the experiments. Differential gene expression caused by differences in recent evolutionary history exceeded that caused by immediate rearing host. Together, the three genomic datasets suggest that genes affecting traits other than weight and development time are likely to be the main causes of parallel evolution and that detoxification genes (especially cytochrome P450s and beta-glucosidase) could be especially important for colonization of lentil by C. maculatus.

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.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee.

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Coevolution of Genome Architecture and Social Behavior.

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

A split sex ratio in solitary and social nests of a facultatively social bee.

A classic prediction of kin selection theory is that a mixed population of social and solitary nests of haplodiploid insects should exhibit a split sex ratio among offspring: female biased in social nests, male biased in solitary nests. Here, we provide the first evidence of a solitary-social split sex ratio, using the sweat bee Megalopta genalis (Halictidae). Data from 2502 offspring collected from naturally occurring nests across 6 years spanning the range of the M. genalis reproductive season show that despite significant yearly and seasonal variation, the offspring sex ratio of social nests is consistently more female biased than in solitary nests. This suggests that split sex ratios may facilitate the evolutionary origins of cooperation based on reproductive altruism via kin selection.

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.

Phylogeographic and diversification patterns of the white-nosed coati (Nasua narica): Evidence for south-to-north colonization of North America.

White-nosed coatis (Nasua narica) are widely distributed throughout North, Central, and South America, but the patterns of temporal and spatial diversification that have contributed to this distribution are unknown. In addition, the biogeographic history of procyonid species in the Americas remains contentious. Using sequences from three mitochondrial loci (Cytochrome b, NAHD5 and 16S rRNA; 2201 bp) and genotypes from 11 microsatellite loci, we analyzed genetic diversity to determine phylogeographic patterns, genetic structure, divergence times, and gene flow among Nasua narica populations throughout the majority of the species' range. We also estimated the ancestral geographic range of N. narica and other procyonid species. We found a high degree of genetic structure and divergence among populations that conform to five evolutionarily significant units. The most southerly distributed population (Panama) branched off much earlier (∼3.8 million years ago) than the northern populations (<1.2 million years ago). Estimated gene flow among populations was low and mostly northwards and westwards. The phylogeographic patterns within N. narica are associated with geographic barriers and habitat shifts likely caused by Pliocene-Pleistocene climate oscillations. Significantly, our findings suggest the dispersal of N. narica was south-to-north beginning in the Pliocene, not in the opposite direction during the Pleistocene as suggested by the fossil record, and that the most recent common ancestor for coati species was most likely distributed in South or Central America six million years ago. Our study implies the possibility that the diversification of Nasua species, and other extant procyonid lineages, may have occurred in South America.

Investigating the viral ecology of global bee communities with high-throughput metagenomics.

Bee viral ecology is a fascinating emerging area of research: viruses exert a range of effects on their hosts, exacerbate impacts of other environmental stressors, and, importantly, are readily shared across multiple bee species in a community. However, our understanding of bee viral communities is limited, as it is primarily derived from studies of North American and European Apis mellifera populations. Here, we examined viruses in populations of A. mellifera and 11 other bee species from 9 countries, across 4 continents and Oceania. We developed a novel pipeline to rapidly and inexpensively screen for bee viruses. This pipeline includes purification of encapsulated RNA/DNA viruses, sequence-independent amplification, high throughput sequencing, integrated assembly of contigs, and filtering to identify contigs specifically corresponding to viral sequences. We identified sequences for (+)ssRNA, (-)ssRNA, dsRNA, and ssDNA viruses. Overall, we found 127 contigs corresponding to novel viruses (i.e. previously not observed in bees), with 27 represented by >0.1% of the reads in a given sample, and 7 contained an RdRp or replicase sequence which could be used for robust phylogenetic analysis. This study provides a sequence-independent pipeline for viral metagenomics analysis, and greatly expands our understanding of the diversity of viruses found in bee communities.

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.

Juvenile hormone, but not nutrition or social cues, affects reproductive maturation in solitary alkali bees (Nomia melanderi).

Eusocial insect colonies are defined by extreme variation in reproductive activity among castes, but the ancestral conditions from which this variation arose are unknown. Investigating the factors that contribute to variation in reproductive physiology among solitary insects that are closely related to social species can help to fill this gap. We experimentally tested the role of nutrition, juvenile hormone (JH) and social cues on reproductive maturation in solitary alkali bees (Halictidae: Nomia melanderi). We found that alkali bee females emerge from overwintering with small Dufour's glands and small ovaries, containing oocytes in the early stages of development. Oocyte maturation occurs rapidly, and is staggered between the two ovaries. Lab-reared females reached reproductive maturity without access to mates or nesting opportunities, and many had resorbed oocytes. Initial activation of these reproductive structures does not depend on pollen consumption, though dietary protein or lipids may be necessary for long-term reproductive activity. JH is likely to be a limiting factor in alkali bee reproductive activation, as females treated with JH were more likely to develop mature oocytes and Dufour's glands. Unlike for related social bees, the effects of JH were not suppressed by the presence of older, reproductive females. These results provide valuable insight into the factors that influence reproductive activity in an important native pollinator, and those that may have been particularly influential in the evolution of reproductive castes.

Carrion fly-derived DNA metabarcoding is an effective tool for mammal surveys: Evidence from a known tropical mammal community.

Metabarcoding of vertebrate DNA derived from carrion flies has been proposed as a promising tool for biodiversity monitoring. To evaluate its efficacy, we conducted metabarcoding surveys of carrion flies on Barro Colorado Island (BCI), Panama, which has a well-known mammal community, and compared our results against diurnal transect counts and camera trapping. We collected 1,084 flies in 29 sampling days, conducted metabarcoding with mammal-specific (16S) and vertebrate-specific (12S) primers, and sequenced amplicons on Illumina MiSeq. For taxonomic assignment, we compared blast with the new program protax, and we found that protax improved species identifications. We detected 20 mammal, four bird, and one lizard species from carrion fly metabarcoding, all but one of which are known from BCI. Fly metabarcoding detected more mammal species than concurrent transect counts (29 sampling days, 13 species) and concurrent camera trapping (84 sampling days, 17 species), and detected 67% of the number of mammal species documented by 8 years of transect counts and camera trapping combined, although fly metabarcoding missed several abundant species. This study demonstrates that carrion fly metabarcoding is a powerful tool for mammal biodiversity surveys and has the potential to detect a broader range of species than more commonly used methods.