Search

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real time and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after it was brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate fixation of traits beneficial for exploitation of the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

Background The honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity. Results For both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4 °C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms. Conclusions Ethanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

The study of island community assembly has been fertile ground for developing and testing theoretical ideas in ecology and evolution. The ecoevolutionary trajectory of lineages after colonization has been a particular interest, as this is a key component of understanding community assembly. In this system, existing ideas, such as the taxon cycle, posit that lineages pass through a regular sequence of ecoevolutionary changes after colonization, with lineages shifting toward reduced dispersal ability, increased ecological specialization, and declines in abundance. However, these predictions have historically been difficult to test. Here, we integrate phylogenomics, population genomics, and X‐ray microtomography/3D morphometrics, to test hypotheses for whether the ecomorphological diversity of trap‐jaw ants (Strumigenys) in the Fijian archipelago is assembled primarily through colonization or postcolonization radiation, and whether species show ecological shifts toward niche specialization, toward upland habitats, and decline in abundance after colonization. We infer that most Fijian endemic Strumigenys evolved in situ from a single colonization and have diversified to fill a large fraction of global morphospace occupied by the genus. Within this adaptive radiation, lineages trend to different degrees toward high elevation, reduced dispersal ability, and demographic decline, and we find no evidence of repeated colonization that displaces the initial radiation. Overall these results are only partially consistent with taxon cycle and associated ideas, while highlighting the potential role of priority effects in assembling island communities.

While ectoparasitic Varroa mites cause minimal damage to their co-evolved ancestral host, the eastern honey bee (Apis cerana), they devastate their novel host, the western honey bee (Apis mellifera). Over several decades, the host switch caused worldwide population collapses, threatening global food security. Varroa management strategies have focused on breeding bees for tolerance. But, can Varroa overcome these counter-adaptations in a classic coevolutionary arms race? Despite increasing evidence for Varroa genetic diversity and evolvability, this eventuality has largely been neglected. We therefore suggest a more holistic paradigm for studying this host-parasite interaction, one in which ‘Varroa-tolerant’ bee traits should be viewed as a shared phenotype resulting from Varroa and honey bee interaction.

Microbiomes provide a range of benefits to their hosts which can lead to the coevolution of a joint ecological niche. However, holometabolous insects, some of the most successful organisms on Earth, occupy different niches throughout development, with larvae and adults being physiologically and morphologically highly distinct. Furthermore, transition between the stages usually involves the loss of the gut microbiome since the gut is remodeled during pupation. Most eusocial organisms appear to have evolved a workaround to this problem by sharing their communal microbiome across generations. However, whether this vertical microbiome transmission can overcome perturbations of the larval microbiome remains untested. Honey bees have a relatively simple, conserved, coevolved adult microbiome which is socially transmitted and affects many aspects of their biology. In contrast, larval microbiomes are more variable, with less clear roles. Here, we manipulated the gut microbiome of in vitro-reared larvae, and after pupation of the larvae, we inoculated the emerged bees with adult microbiome to test whether adult and larval microbiome stages may be coupled (e.g., through immune priming). Larval treatments differed in bacterial composition and abundance, depending on diet, which also drove larval gene expression. Nonetheless, adults converged on the typical core taxa and showed limited gene expression variation. This work demonstrates that honey bee adult and larval stages are effectively microbiologically decoupled, and the core adult microbiome is remarkably stable to early developmental perturbations. Combined with the transmission of the microbiome in early adulthood, this allows the formation of long-term host-microbiome associations. IMPORTANCE This work investigated host-microbiome interactions during a crucial developmental stage—the transition from larvae to adults, which is a challenge to both, the insect host and its microbiome. Using the honey bee as a tractable model system, we showed that microbiome transfer after emergence overrides any variation in the larvae, indicating that larval and adult microbiome stages are effectively decoupled. Together with the reliable vertical transfer in the eusocial system, this decoupling ensures that the adults are colonized with a consistent and derived microbiome after eclosion. Taken all together, our data provide additional support that the evolution of sociality, at least in the honey bee system tested here, is linked with host-microbiome relationships.

In this study, we report the first complete mitochondrial genome sequence of the Aquatic Coralsnake Micrurus surinamensis. The mitochondrial genome lengthis 17,375 bp, comprising 13 protein-coding genes, 2 rRNA (12S and 16S) and 22 tRNA, as well as two typical control regions. Phylogenetic analysis based upon 13 protein-coding genes showed clusters based on terrestrial and marine species.

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities., Maeva Techer et al. report genome assemblies of two honeybee parasitic mites, Varroa destructor and V. jacobsoni. They find that 4.8% of orthologous genes show evidence of positive selection in at least one species, though the genes under selection are distinct between species indicating divergent evolution.

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of host-symbiont interactions between the northernmost fungus-farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome-wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant-fungus genome-genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population-genetic structure was concordant between the ants and one cultivar type (M-fungi, concordant clines) but discordant for the other cultivar type (T-fungi). Discordance in population-genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between-nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant-fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.

Traditional static publication formats make visualization, exploration, and sharing of massive phylogenetic trees difficult. A phylogenetic study often involves hundreds of taxa, and the resulting tree has to be split across multiple journal pages, or be shrunk onto one, which jeopardizes legibility. Furthermore, additional data layers, such as species-specific information or time calibrations are often displayed in separate figures, making the entire picture difficult for readers to grasp. Web-based technologies, such as the Data Driven Document (D3) JavaScript library, were created to overcome such challenges by allowing interactive displays of complex data sets. The new phylogeny.IO web server (https://phylogeny.io) overcomes this issue by allowing users to easily import, annotate, and share interactive phylogenetic trees. It allows a range of static (e.g. such as shapes and colors) and dynamic (e.g. pop-up text and images) annotations. Annotated trees can be saved on the server for subsequent modification or they may be shared as IFrame HTML objects, easily embeddable in any web page. The principal goal of phylogeny.IO is not to produce publication-ready figures, but rather to provide a simple and intuitive annotation interface that allows easy and rapid sharing of figures in blogs, lecture notes, press releases, etc.

BackgroundThe ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybeeApis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.ResultsEach of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor >98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds ofAvaIandAluIrepeats associated with these features.ConclusionsThe improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Evolutionary innovations underlie the rise of diversity and complexity—the 2 long-term trends in the history of life. How does natural selection redesign multiple interacting parts to achieve a new emergent function? We investigated the evolution of a biomechanical innovation, the latch-spring mechanism of trap-jaw ants, to address 2 outstanding evolutionary problems: how form and function change in a system during the evolution of new complex traits, and whether such innovations and the diversity they beget are repeatable in time and space. Using a new phylogenetic reconstruction of 470 species, and X-ray microtomography and high-speed videography of representative taxa, we found the trap-jaw mechanism evolved independently 7 to 10 times in a single ant genus (Strumigenys), resulting in the repeated evolution of diverse forms on different continents. The trap mechanism facilitates a 6 to 7 order of magnitude greater mandible acceleration relative to simpler ancestors, currently the fastest recorded acceleration of a resettable animal movement. We found that most morphological diversification occurred after evolution of latch-spring mechanisms, which evolved via minor realignments of mouthpart structures. This finding, whereby incremental changes in form lead to a change of function, followed by large morphological reorganization around the new function, provides a model for understanding the evolution of complex biomechanical traits, as well as insights into why such innovations often happen repeatedly., The mousetrap-like mandibles of trap-jaw ants are a functional innovation capable of record-breaking acceleration. This study uses phylogenomics, X-ray microtomography and high-speed videography to reveal that they evolved multiple times around the world through a pathway of incremental steps.

BackgroundThe honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity.ResultsFor both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4°C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms.ConclusionEthanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

Leafcutter ants propagate co-evolving fungi for food. The nearly 50 species of leafcutter ants (Atta, Acromyrmex) range from Argentina to the United States, with the greatest species diversity in southern South America. We elucidate the biogeography of fungi cultivated by leafcutter ants using DNA sequence and microsatellite-marker analyses of 474 cultivars collected across the leafcutter range. Fungal cultivars belong to two clades (Clade-A and Clade-B). The dominant and widespread Clade-A cultivars form three genotype clusters, with their relative prevalence corresponding to southern South America, northern South America, Central and North America. Admixture between Clade-A populations supports genetic exchange within a single species, Leucocoprinus gongylophorus. Some leafcutter species that cut grass as fungicultural substrate are specialized to cultivate Clade-B fungi, whereas leafcutters preferring dicot plants appear specialized on Clade-A fungi. Cultivar sharing between sympatric leafcutter species occurs frequently such that cultivars of Atta are not distinct from those of Acromyrmex. Leafcutters specialized on Clade-B fungi occur only in South America. Diversity of Clade-A fungi is greatest in South America, but minimal in Central and North America. Maximum cultivar diversity in South America is predicted by the Kusnezov-Fowler hypothesis that leafcutter ants originated in subtropical South America and only dicot-specialized leafcutter ants migrated out of South America, but the cultivar diversity becomes also compatible with a recently proposed hypothesis of a Central American origin by postulating that leafcutter ants acquired novel cultivars many times from other nonleafcutter fungus-growing ants during their migrations from Central America across South America. We evaluate these biogeographic hypotheses in the light of estimated dates for the origins of leafcutter ants and their cultivars.

The Lord Howe Island stick insect, Dryococelus australis, was once common on the island but was driven to extinction after the arrival of ship rats in the early 20th century [1, 2]. It was thought to be extinct for decades, until a tiny population of similar-looking stick insects was discovered 20 km away, on the islet of Ball's Pyramid, in 2001 [2]. Individuals from this population are currently being reared in Australia and elsewhere in the world, with the eventual goal of recolonizing Lord Howe Island [3]. Recent surveys of the wild population on Ball's Pyramid suggest that it is among the world's rarest species. However, there are significant morphological differences between Ball's Pyramid and museum specimens, and there has never been a genetic confirmation of the rediscovered population's species identity. Because Dryococelus is monotypic, there are also no known extant relatives for comparison. Using shotgun genomic data from the Ball's Pyramid population, we assembled a draft genome and the complete mitochondrial genome. We found that the genome is massive, over 4 Gb in size, and is most likely hexaploid. We re-sequenced mitochondrial genomes from historic museum specimens collected on Lord Howe Island before the extinction event. Sequence divergence between the two populations is less than 1% and is within the range of intraspecific differences between the museum specimens, suggesting that they are conspecific and that D. australis has successfully evaded extinction so far. This work highlights the importance of museum collections for taxonomic validation in the context of ongoing conservation efforts.

Venoms are among the most biologically active secretions known, and are commonly believed to evolve under extreme positive selection. Many venom gene families, however, have undergone duplication, and are often deployed in doses vastly exceeding the LD50 for most prey species, which should reduce the strength of positive selection. Here, we contrast these selective regimes using snake venoms, which consist of rapidly evolving protein formulations. Though decades of extensive studies have found that snake venom proteins are subject to strong positive selection, the greater action of drift has been hypothesized, but never tested. Using a combination of de novo genome sequencing, population genomics, transcriptomics, and proteomics, we compare the two modes of evolution in the pitviper, Protobothrops mucrosquamatus. By partitioning selective constraints and adaptive evolution in a McDonald–Kreitman-type framework, we find support for both hypotheses: venom proteins indeed experience both stronger positive selection, and lower selective constraint than other genes in the genome. Furthermore, the strength of selection may be modulated by expression level, with more abundant proteins experiencing weaker selective constraint, leading to the accumulation of more deleterious mutations. These findings show that snake venoms evolve by a combination of adaptive and neutral mechanisms, both of which explain their extraordinarily high rates of molecular evolution. In addition to positive selection, which optimizes efficacy of the venom in the short term, relaxed selective constraints for deleterious mutations can lead to more rapid turnover of individual proteins, and potentially to exploration of a larger venom phenotypic space.

Development of polymorphic phenotypes from similar genomes requires gene expression differences. However, little is known about how morph-specific gene expression patterns vary on a broad phylogenetic scale. We hypothesize that evolution of morph-specific gene expression, and consequently morph-specific phenotypic evolution, may be constrained by gene essentiality and the amount of pleiotropic constraints. Here, we use comparative transcriptomics of queen and worker morphs, that is, castes, from 15 ant species to understand the constraints of morph-biased gene expression. In particular, we investigate how measures of evolutionary constraints at the sequence level (expression level, connectivity, and number of gene ontology [GO] terms) correlate with morph-biased expression. Our results show that genes indeed vary in their potential to become morph-biased. The existence of genes that are constrained in becoming caste-biased potentially limits the evolutionary decoupling of the caste phenotypes, that is, it might result in "caste load" occasioning from antagonistic fitness variation, similarly to sexually antagonistic fitness variation between males and females. On the other hand, we suggest that genes under low constraints are released from antagonistic variation and thus more likely to be co-opted for morph specific use. Overall, our results suggest that the factors that affect sequence evolutionary rates and evolution of plastic expression may largely overlap.

Kin selection is thought to drive the evolution of cooperation and conflict, but the specific genes and genome-wide patterns shaped by kin selection are unknown. We identified thousands of genes associated with the sterile ant worker caste, the archetype of an altruistic phenotype shaped by kin selection, and then used population and comparative genomic approaches to study patterns of molecular evolution at these genes. Consistent with population genetic theoretical predictions, worker-upregulated genes showed relaxed adaptive evolution compared to genes upregulated in reproductive castes. Worker-upregulated genes included more taxonomically-restricted genes, indicating that the worker caste has recruited more novel genes, yet these genes also showed relaxed selection. Our study identifies a putative genomic signature of kin selection and helps to integrate emerging sociogenomic data with longstanding social evolution theory.

Museums and herbaria represent vast repositories of biological material. Until recently, working with these collections has been difficult, due to the poor condition of historical DNA. However, recent advances in next-generation sequencing technology, and subsequent development of techniques for preparing and sequencing historical DNA, have recently made working with collection specimens an attractive option. Here we describe the unique technical challenges of working with collection specimens, and innovative molecular methods developed to tackle them. We also highlight possible applications of collection specimens, for taxonomy, ecology and evolution. The application of next-generation sequencing methods to museum and herbaria collections is still in its infancy. However, by giving researchers access to billions of specimens across time and space, it holds considerable promise for generating future discoveries across many fields.

Island biodiversity has long fascinated biologists as it typically presents tractable systems for unpicking the eco-evolutionary processes driving community assembly. In general, two recurring themes are of central theoretical interest. First, immigration, diversification, and extinction typically depend on island geographical properties (e.g., area, isolation, and age). Second, predictable ecological and evolutionary trajectories readily occur after colonization, such as the evolution of adaptive trait syndromes, trends toward specialization, adaptive radiation, and eventual ecological decline. Hypotheses such as the taxon cycle draw on several of these themes to posit particular constraints on colonization and subsequent eco-evolutionary dynamics. However, it has been challenging to examine these integrated dynamics with traditional methods. Here, we combine phylogenomics, population genomics and phenomics, to unravel community assembly dynamics among Pheidole (Hymenoptera, Formicidae) ants in the isolated Fijian archipelago. We uphold basic island biogeographic predictions that isolated islands accumulate diversity primarily through in situ evolution rather than dispersal, and population genomic support for taxon cycle predictions that endemic species have decreased dispersal ability and demography relative to regionally widespread taxa. However, rather than trending toward island syndromes, ecomorphological diversification in Fiji was intense, filling much of the genus-level global morphospace. Furthermore, while most endemic species exhibit demographic decline and reduced dispersal, we show that the archipelago is not an evolutionary dead-end. Rather, several endemic species show signatures of population and range expansion, including a successful colonization to the Cook islands. These results shed light on the processes shaping island biotas and refine our understanding of island biogeographic theory.

Eusociality has convergently evolved multiple times, but the genomic basis of caste-based division of labor and degree to which independent origins of eusociality have utilized common genes remain largely unknown. Here we characterize caste-specific transcriptomic profiles across development and adult body segments from pharaoh ants (Monomorium pharaonis) and honey bees (Apis mellifera), representing two independent origins of eusociality. We identify a substantial shared core of genes upregulated in the abdomens of queen ants and honey bees that also tends to be upregulated in mated female flies, suggesting that these genes are part of a conserved insect reproductive groundplan. Outside of this shared groundplan, few genes are differentially expressed in common. Instead, the majority of the thousands of caste-associated genes are plastically expressed, rapidly evolving, and relatively evolutionarily young. These results emphasize that the recruitment of both highly conserved and lineage-specific genes underlie the convergent evolution of novel traits such as eusociality., Eusocial caste systems have evolved independently multiple times. Here, Warner et al. investigate the amount of shared vs. lineage-specific genes involved in the evolution of caste in pharaoh ants and honey bees by comparing transcriptomes across tissues, developmental stages, and castes.

Development is often strongly regulated by interactions among close relatives, but the underlying molecular mechanisms are largely unknown. In eusocial insects, interactions between caregiving worker nurses and larvae regulate larval development and resultant adult phenotypes. Here, we begin to characterize the social interactome regulating ant larval development by collecting and sequencing the transcriptomes of interacting nurses and larvae across time. We find that the majority of nurse and larval transcriptomes exhibit parallel expression dynamics across larval development. We leverage this widespread nurse-larva gene co-expression to infer putative social gene regulatory networks acting between nurses and larvae. Genes with the strongest inferred social effects tend to be peripheral elements of within-tissue regulatory networks and are often known to encode secreted proteins. This includes interesting candidates such as the nurse-expressed giant-lens, which may influence larval epidermal growth factor signaling, a pathway known to influence various aspects of insect development. Finally, we find that genes with the strongest signatures of social regulation tend to experience relaxed selective constraint and are evolutionarily young. Overall, our study provides a first glimpse into the molecular and evolutionary features of the social mechanisms that regulate all aspects of social life., Author summary Social interactions are fundamental to all forms of life, from single-celled bacteria to complex plants and animals. Despite their obvious importance, little is known about the molecular causes and consequences of social interactions. In this paper, we study the molecular basis of nurse-larva social interactions that regulate larval development in the pharaoh ant Monomorium pharaonis. We infer the effects of social interactions on gene expression from samples of nurses and larvae collected in the act of interaction across a developmental time series. Gene expression appears to be closely tied to these interactions, such that we can identify genes expressed in nurses with putative regulatory effects on larval gene expression. Genes which we infer to have strong social regulatory effects tend to have weak regulatory effects within individuals, and highly social genes tend to experience relatively weaker natural selection in comparison to fewer social genes. This study represents a novel approach and foundation upon which future studies at the intersection of genetics, behavior, and evolution can build.

Queen pheromones are chemical signals that mediate reproductive division of labor in eusocial animals. Remarkably, queen pheromones are composed of identical or chemically similar compounds in some ants, wasps and bees, even though these taxa diverged >150MYA and evolved queens and workers independently. Here, we measure the transcriptomic consequences of experimental exposure to queen pheromones in workers from two ant and two bee species (genera: Lasius, Apis, Bombus), and test whether they are similar across species. Queen pheromone exposure affected transcription and splicing at many loci. Many genes responded consistently in multiple species, and the set of pheromone-sensitive genes was enriched for functions relating to lipid biosynthesis and transport, olfaction, production of cuticle, oogenesis, and histone (de)acetylation. Pheromone-sensitive genes tend to be evolutionarily ancient, positively selected, peripheral in the gene coexpression network, hypomethylated, and caste-specific in their expression. Our results reveal how queen pheromones achieve their effects, and suggest that ants and bees use similar genetic modules to achieve reproductive division of labor., Queen pheromones are used by eusocial insects to regulate all aspects of colony life. Here, Holman et al. compare the effects of queen pheromone on gene expression and splicing in four eusocial insect species, giving insight into the mechanism and evolution of division of reproductive labour.

Social insects provide systems for studying epigenetic regulation of phenotypes, particularly with respect to differentiation of reproductive and worker castes, which typically arise from a common genetic background. The role of gene expression in caste specialization has been extensively studied, but the role of DNA methylation remains controversial. Here, we perform well-replicated, integrated analyses of DNA methylation and gene expression in brains of an ant (Formica exsecta) with distinct female castes using traditional approaches (tests of differential methylation) combined with a novel approach (analysis of co-expression and co-methylation networks). We found differences in expression and methylation profiles between workers and queens at different life stages, as well as some overlap between DNA methylation and expression at the functional level. Large portions of the transcriptome and methylome are organized into‘modules’of genes, some significantly associated with phenotypic traits of castes and developmental stages. Several gene co-expression modules are preserved in co-methylation networks, consistent with possible regulation of caste-specific gene expression by DNA methylation. Surprisingly, brain co-expression modules were highly preserved when compared with a previous study that examined whole-body co-expression patterns in 16 ant species (Morandin et al. 2016), suggesting that these modules are evolutionarily conserved and for specific functions in various tissues. Altogether, these results suggest that DNA methylation participates in regulation of caste specialization and age-related physiological changes in social insects.

From human crowds to cells in tissue, the detection and efficient tracking of multiple objects in dense configurations is an important and unsolved problem. In the past, limitations of image analysis have restricted studies of dense groups to tracking a single or subset of marked individuals, or to coarse-grained group-level dynamics, all of which yield incomplete information. Here, we combine convolutional neural networks (CNNs) with the model environment of a honeybee hive to automatically recognize all individuals in a dense group from raw image data. We create new, adapted individual labeling and use the segmentation architecture U-Net with a loss function dependent on both object identity and orientation. We additionally exploit temporal regularities of the video recording in a recurrent manner and achieve near human-level performance while reducing the network size by 94% compared to the original U-Net architecture. Given our novel application of CNNs, we generate extensive problem-specific image data in which labeled examples are produced through a custom interface with Amazon Mechanical Turk. This dataset contains over 375,000 labeled bee instances across 720 video frames at 2 FPS, representing an extensive resource for the development and testing of tracking methods. We correctly detect 96% of individuals with a location error of ~7% of a typical body dimension, and orientation error of 12 degrees, approximating the variability of human raters. Our results provide an important step towards efficient image-based dense object tracking by allowing for the accurate determination of object location and orientation across time-series image data efficiently within one network architecture., 15 pages, including supplementary figures. 1 supplemental movie available as an ancillary file

Early in the process of adaptive radiation, allopatric disruption of gene flow followed by ecological specialization is key for speciation; but, do adaptive radiations occur on small islands without internal geographical barriers? Island populations sometimes harbour polymorphism in ecological specializations, but its significance remains unclear. On one hand, morphs may correspond to ‘cryptic’ species. Alternatively, they could result from population, developmental or behavioural plasticity. The spider Wendilgarda galapagensis (Araneae, Theridiosomatidae) is endemic to the small Isla del Coco and unique in spinning three different web types, each corresponding to a different microhabitat. We tested whether this variation is associated with ‘cryptic’ species or intraspecific behavioural plasticity. Despite analysing 36 803 loci across 142 individuals, we found no relationship between web type and population structure, which was only weakly geographically differentiated. The same pattern holds when looking within a sampling site or considering only F st outliers. In line with genetic data, translocation experiments showed that web architecture is plastic within an individual. However, not all transitions between web types are equally probable, indicating the existence of individual preferences. Our data supports the idea that diversification on small islands might occur mainly at the behavioural level producing an intraspecific niche partition without speciation.

The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

A key requirement for social cooperation is the mitigation and/or social regulation of aggression towards other group members. Populations of the harvester ant Pogonomyrmex californicus show the alternate social phenotypes of queens founding nests alone (haplometrosis) or in groups of unrelated yet cooperative individuals (pleometrosis). Pleometrotic queens display an associated reduction in aggression. To understand the proximate drivers behind this variation, we placed foundresses of the two populations into social environments with queens from the same or the alternate population, and measured their behaviour and head gene expression profiles. A proportion of queens from both populations behaved aggressively, but haplometrotic queens were significantly more likely to perform aggressive acts, and conflict escalated more frequently in pairs of haplometrotic queens. Whole-head RNA sequencing revealed variation in gene expression patterns, with the two populations showing moderate differentiation in overall transcriptional profile, suggesting that genetic differences underlie the two founding strategies. The largest detected difference, however, was associated with aggression, regardless of queen founding type. Several modules of coregulated genes, involved in metabolism, immune system and neuronal function, were found to be upregulated in highly aggressive queens. Conversely, nonaggressive queens exhibited a striking pattern of upregulation in chemosensory genes. Our results highlight that the social phenotypes of cooperative vs. solitary nest founding tap into a set of gene regulatory networks that seem to govern aggression level. We also present a number of highly connected hub genes associated with aggression, providing opportunity to further study the genetic underpinnings of social conflict and tolerance.

Corals belong to the most basal class of the Phylum Cnidaria, which is considered the sister group of bilaterian animals, and thus have become an emerging model to study the evolution of developmental mechanisms. Although cell renewal, differentiation, and maintenance of pluripotency are cellular events shared by multicellular animals, the cellular basis of these fundamental biological processes are still poorly understood. To understand how changes in gene expression regulate morphogenetic transitions at the base of the eumetazoa, we performed quantitative RNA-seq analysis during Acropora digitifera’s development. We collected embryonic, larval, and adult samples to characterize stage-specific transcription profiles, as well as broad expression patterns. Transcription profiles reconstructed development revealing two main expression clusters. The first cluster grouped blastula and gastrula and the second grouped subsequent developmental time points. Consistently, we observed clear differences in gene expression between early and late developmental transitions, with higher numbers of differentially expressed genes and fold changes around gastrulation. Furthermore, we identified three coexpression clusters that represented discrete gene expression patterns. During early transitions, transcriptional networks seemed to regulate cellular fate and morphogenesis of the larval body. In late transitions, these networks seemed to play important roles preparing planulae for switch in lifestyle and regulation of adult processes. Although developmental progression in A. digitifera is regulated to some extent by differential coexpression of well-defined gene networks, stage-specific transcription profiles appear to be independent entities. While negative regulation of transcription is predominant in early development, cell differentiation was upregulated in larval and adult stages.

Stony corals (Scleractinia) form the basis for some of the most diverse ecosytems on Earth, but we have much to learn about their evolutionary history and systematic relationships. In order to improve our understanding of species in corals we here investigated phylogenetic relationships between morphologically defined species and genetic lineages in the genus Galaxea (Euphyllidae) using a combined phylogenomic and phylogeographic approach. Previous studies revealed the nominal species G. fascicularis included three genetically well-differentiated lineages (L, S & L+) in the western Pacific, but their distribution and relationship to other species in the genus was unknown. Based on genomic (RAD-seq) and mitochondrial sequence data (non-coding region between cytb and ND2) we investigated whether the morphological taxa represent genetically coherent entities and what is the phylogenetic relationship and spatial distribution of the three lineages of G. fascicularis throughout the observed species range. Using the RAD-seq data, we find that the genus Galaxea is monophyletic and contains three distinct clades: an Indo-Pacific, a Pacific, and a small clade restricted to the Chagos Archipelago. The three lineages of G. fascicularis were associated with different RAD-seq clades, with the ‘L’ lineage showing some morphological distinction from the other two lineages (larger more asymmetrical polyps). In addition to these, three more genetic lineages in G. fascicularis may be distinguished – a Chagossian, an Ogasawaran, and one from the Indian-Red Sea. Among nominal taxa for which we have multiple samples, G. horrescens was the only monophyletic species. The mitochondrial non-coding region is highly conserved apart of the length polymorphism used to define L, S & L+ lineages and lacks the power to distinguish morphological and genetic groups resolved with genomic RAD-sequencing. The polyphyletic nature of most species warrants a careful examination of the accepted taxonomy of this group with voucher collections and their comparison to type specimens to resolve species boundaries. Further insight to the speciation process in corals will require international cooperation for the sharing of specimens to facilitate scientific discovery.

Summary A basic expectation of evolution by natural selection is that species morphologies will adapt to their ecological niche. In social organisms, this may include selective pressure from the social environment. Many non-ant parasites of ant colonies are known to mimic the morphology of their host species, often in striking fashion [ 1 , 2 ], indicating there is selection on parasite morphology to match the host (Batesian and/or Wasmannian mimicry [ 3 ]). However, ants that parasitize other ant societies are usually closely related to their hosts (Emery’s rule) [ 4 , 5 , 6 , 7 , 8 ] and expected to be similar due to common ancestry, making any kind of mimicry difficult to detect [ 9 ]. Here, we investigate the diversification of the hyperdiverse ant genus Pheidole in Madagascar, including the evolution of 13 putative social parasite species within a broader radiation of over 100 ant species on the island. We find that the parasitic species are monophyletic and that their associated hosts are spread across the Malagasy Pheidole radiation. This provides an opportunity to test for selection on morphological similarity and divergence between parasites and hosts. Using X-ray microtomography and both linear measurements and three-dimensional (3D) geometric morphometrics, we show that ant social parasite worker morphologies feature a mix of “host-matching” and “parasitic” traits, where the former converge on the host phenotype and the latter diverge from typical Pheidole phenotypes to match a common parasitic syndrome. This finding highlights the role of social context in shaping the evolution of phenotypes and raises questions about the role of morphological sensing in nestmate recognition.

Microbes ubiquitously inhabit animals and plants, often affecting their host’s phenotype. As a result, even in a constant genetic background, the host’s phenotype may evolve through indirect selection on the microbiome. ‘Microbiome engineering’ offers a promising novel approach for attaining desired host traits but has been attempted only a few times. Building on the known role of the microbiome on development in fruit flies, we attempted to evolve earlier-eclosing flies by selecting on microbes in the growth media. We carried out parallel evolution experiments in no- and high-sugar diets by transferring media associated with fast-developing fly lines over the course of four selection cycles. In each cycle, we used sterile eggs from the same inbred population, and assayed mean fly eclosion times. Ultimately, flies eclosed seven to twelve hours earlier, depending on the diet, but microbiome engineering had no effect relative to a random-selection control treatment. 16S rRNA gene sequencing showed that the microbiome did evolve, particularly in the no sugar diet, with an increase in Shannon diversity over time. Thus, while microbiome evolution did affect host eclosion times, these effects were incidental. Instead, any experimentally enforced selection effects were swamped by uncontrolled microbial evolution, likely resulting in its adaptation to the media. These results imply that selection on host phenotypes must be strong enough to overcome other selection pressures simultaneously operating on the microbiome. The independent evolutionary trajectories of the host and the microbiome may limit the extent to which indirect selection on the microbiome can ultimately affect host phenotype. Random-selection lines accounting for independent microbial evolution are essential for experimental microbiome engineering studies.

Key innovations provide ecological opportunity by enabling access to new resources, colonization of new environments, and are associated with adaptive radiation. The most well-known pattern associated with adaptive radiation is an early burst of phenotypic diversification. Venoms facilitate prey capture and are widely believed to be key innovations leading to adaptive radiation. However, few studies have estimated their evolutionary rate dynamics. Here, we test for patterns of adaptive evolution in venom gene expression data from 52 venomous snake species. By identifying shifts in tempo and mode of evolution along with models of phenotypic evolution, we show that snake venom exhibits the macroevolutionary dynamics expected of key innovations. Namely, all toxin families undergo shifts in their rates of evolution, likely in response to changes in adaptive optima. Furthermore, we show that rapid-pulsed evolution modelled as a Lévy process better fits snake venom evolution than conventional early burst or Ornstein–Uhlenbeck models. While our results support the idea of snake venom being a key innovation, the innovation of venom chemistry lacks clear mechanisms that would lead to reproductive isolation and thus adaptive radiation. Therefore, the extent to which venom directly influences the diversification process is still a matter of contention.

To explore landscape genomics at the range limit of an obligate mutualism, we used genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of hostsymbiont interactions between the northernmost fungus-farming leafcutter antAtta texanaand its two main types of cultivated fungus. At local scales, genome-wide differentiation between ants associated with either of the two fungal types is greater than the differentiation associated with the abiotic factors temperature and precipitation, suggesting that specific ant-fungus genome-genome combinations may have been favored by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. In contrast, genetic structure was patchy in the cultivated fungi, with no consistent reduction of fungal genetic diversity at the range margins. This discordance in population-genetic structure between ant hosts and fungal symbionts is surprising because the ant farmers co-disperse with their vertically-transmitted fungal symbionts, but apparently the fungi disperse occasionally also through between-nest horizontal transfer or other unknown dispersal mechanisms. The discordance in populationgenetic structure indicates that genetic drift and gene flow differ in magnitude between each partner in this leafcutter mutualism. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influences adaptive genetic differentiation of the host during range evolution in this obligate mutualism.

We dekaryotized the multinucleate fungus Leucocoprinus gongylophorus, a symbiotic fungus cultivated vegetatively by leafcutter ants as their food. To track genetic changes resulting from dekaryotization (elimination of some nuclei from the multinuclear population), we developed two multiplex microsatellite fingerprinting panels (15 loci total), then characterized the allele profiles of 129 accessions generated by dekaryotization treatment. Genotype profiles of the 129 accessions confirmed allele loss expected by dekaryotization of the multinucleate fungus. We found no evidence for haploid and single-nucleus strains among the 129 accessions. Microscopy of fluorescently stained dekaryotized accessions revealed great variation in nuclei number between cells of the same vegetative mycelium, with cells containing typically between 3 and 15 nuclei/cell (average = 9.4 nuclei/cell; mode = 8). We distinguish four mycelial morphotypes among the dekaryotized accessions; some of these morphotypes had lost the full competence to produce gongylidia (nutritive hyphal-tip swellings consumed by leafcutter ants as food). In mycelial growth confrontations between different gongylidia-incompetent accessions, allele profiles suggest exchange of nuclei between dekaryotized accessions, restoring full gongylidia competence in some of these strains. The restoration of gongylidia competence after genetic exchange between dekaryotized strains suggests the hypothesis that complementary nuclei interact, or nuclear and cytoplasmic factors interact, to promote or enable gongylidia competence.

A central goal of biology is to uncover the genetic basis for the origin of new phenotypes. A particularly effective approach is to examine the genomic architecture of species that have secondarily lost a phenotype with respect to their close relatives. In the eusocial Hymenoptera, queens and workers have divergent phenotypes that may be produced via either expression of alternative sets of caste-specific genes and pathways or differences in expression patterns of a shared set of multifunctional genes. To distinguish between these two hypotheses, we investigated how secondary loss of the worker phenotype in workerless ant social parasites impacted genome evolution across two independent origins of social parasitism in the ant generaPogonomyrmex and Vollenhovia. We sequenced the genomes of three social parasites and their most-closely related eusocial host species and compared gene losses in social parasites with gene expression differences between host queens and workers. Virtually all annotated genes were expressed to some degree in both castes of the host, with most shifting in queen-worker bias across developmental stages. As a result, despite 4 1M y of divergence from the last common ancestor that had workers, the social parasite ss howed strikingly little evidence of gene loss, damaging mutations, or shifts in selection regime resulting from loss of the worker caste. This suggests that regulatory changes within a multifunctional genome, rather than sequence differences, have played a predominant role in the evolution of social parasitism, and perhaps also in the many gains and losses of phenotypes in the social insects.

We present Bioconda (https://bioconda.github.io), a distribution of bioinformatics software for the lightweight, multiplatform and language-agnostic package manager Conda. Currently, Bioconda offers a collection of over 3000 software packages, which is continuously maintained, updated, and extended by a growing global community of more than 200 contributors. Bioconda improves analysis reproducibility by allowing users to define isolated environments with defined software versions, all of which are easily installed and managed without administrative privileges.

Social insects are emerging models to study how gene regulation affects behavior because their colonies comprise individuals with the same genomes but greatly different behavioral repertoires. To investigate the molecular mechanisms that activate distinct behaviors in different castes, we exploited a natural behavioral plasticity in Harpegnathos saltator, where adult workers can transition to a reproductive, queen-like state called gamergate. Analysis of brain transcriptomes during the transition revealed that corazonin, a neuropeptide homologous to the vertebrate gonadotropin-releasing hormone, was downregulated as workers became gamergates. Corazonin was also preferentially expressed in workers and/or foragers from other social insect species. Injection of corazonin in transitioning Harpegnathos individuals suppressed expression of vitellogenin in the brain and stimulated worker-like hunting behaviors, while inhibiting gamergate behaviors such as dueling and egg deposition. We propose that corazonin is a central regulator of caste identity and behavior in social insects.

While decades of research have focused on snake venom proteins, far less attention has been paid to small organic venom constituents. Using mostly pooled samples, we surveyed 31 venoms (six elapid, six viperid, and 19 crotalid) for spermine, spermidine, putrescine, and cadaverine. Most venoms contained all four polyamines, although some in essentially trace quantities. Spermine is a potentially significant component of many viperid and crotalid venoms (≤0.16% by mass, or 7.9 µmol/g); however, it is almost completely absent from elapid venoms assayed. All elapid venoms contained larger molar quantities of putrescine and cadaverine than spermine, but still at levels that are likely to be biologically insignificant. As with venom purines, polyamines impact numerous physiological targets in ways that are consistent with the objectives of prey envenomation, prey immobilization via hypotension and paralysis. Most venoms probably do not contain sufficient quantities of polyamines to induce systemic effects in prey; however, local effects seem probable. A review of the pharmacological literature suggests that spermine could contribute to prey hypotension and paralysis by interacting with N-methyl-d-aspartate (NMDA) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors, nicotinic and muscarinic acetylcholine receptors, γ-Aminobutyric acid (GABA) receptors, blood platelets, ryanodine receptors, and Ca2+-ATPase. It also blocks many types of cation-permeable channels by interacting with negatively charged amino acid residues in the channel mouths. The site of envenomation probably determines which physiological targets assume the greatest importance; however, venom-induced liberation of endogenous, intracellular stores of polyamines could potentially have systemic implications and may contribute significantly to envenomation sequelae.

Queen-worker differentiation in eusocial organisms may have originated from decoupling of maternal care and reproductive behaviours. Recent advances in sequencing techniques have begun to elucidate the molecular basis of queen-worker differentiation. However, current knowledge of the molecular basis of caste differentiation is limited, especially to species with morphological castes. It seems likely that at the dawn of eusociality morphologically undifferentiated, monomorphic females underwent physiological differentiation that yielded egg-laying and caretaking castes. The molecular basis of such physiological differentiation may provide evolutionary insight into the emergent state of eusociality. In this study, we identify genes associated with monomorphic caste differentiation, specifically focusing on the onset of queen-worker differentiation, using a monomorphic queenless ant, Diacamma sp., that secondarily lost morphological castes. Using individuals experimentally manipulated to become sterile or reproductive, we identified 1546 caste-biased transcripts in brain and 10 in gaster. Because caste differentiation occurs in Diacamma soon after eclosion via behavioural dominance, identified transcripts are interpreted as molecular agents responding immediately to dominance rank formation. Among identified genes, expression levels of genes involved in nutrition processing and storage, such as insulin signalling genes and hexamerins, were strongly altered soon after dominance rank formation. We conclude that the rapid modification of nutrition-related genes in response to social rank may be the fundamental mechanism underlying caste differentiation in Diacamma. Together with functional evidence from the literature, we show that a specific set of genes frequently plays a role in reproductive differentiation across systems with and without morphological castes.

Biological invasions are generally thought to occur after human aided migration to a new range. However, human activities prior to migration may also play a role. We studied here the evolutionary genetics of introduced populations of the invasive ant Wasmannia auropunctata at a worldwide scale. Using microsatellite markers, we reconstructed the main routes of introduction of the species. We found three main routes of introduction, each of them strongly associated to human history and trading routes. We also demonstrate the overwhelming occurrence of male and female clonality in introduced populations of W. auropunctata, and suggest that this particular reproduction system is under selection in human-modified habitats. Together with previous researches focused on native populations, our results suggest that invasive clonal populations may have evolved within human modified habitats in the native range, and spread further from there. The evolutionarily most parsimonious scenario for the emergence of invasive populations of the little fire ant might thus be a two-step process. The W. auropunctata case illustrates the central role of humans in biological change, not only due to changes in migration patterns, but also in selective pressures over species.

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2–6 toxin classes that account for 91–99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A2 (PLA2s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA2s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1–2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (

The sophisticated colony organization of eusocial insects is primarily maintained through the utilization of pheromones. The regulation of these complex social interactions requires intricate chemoreception systems. The recent publication of the genome of Zootermopsis nevadensis opened a new avenue to study molecular basis of termite caste systems. Although there has been a growing interest in the termite chemoreception system that regulates their sophisticated caste system, the relationship between division of labor and expression of chemoreceptor genes remains to be explored. Using high-throughput mRNA sequencing (RNA-seq), we found several chemoreceptors that are differentially expressed among castes and between sexes in a subterranean termite Reticulitermes speratus. In total, 53 chemoreception-related genes were annotated, including 22 odorant receptors, 7 gustatory receptors, 12 ionotropic receptors, 9 odorant-binding proteins, and 3 chemosensory proteins. Most of the chemoreception-related genes had caste-related and sex-related expression patterns; in particular, some chemoreception genes showed king-biased or queen-biased expression patterns. Moreover, more than half of the genes showed significant age-dependent differences in their expression in female and/or male reproductives. These results reveal a strong relationship between the evolution of the division of labor and the regulation of chemoreceptor gene expression, thereby demonstrating the chemical communication and underlining chemoreception mechanism in social insects., 役割によって異なるシロアリの化学受容体の発現 : 高度な分業体制を支えるコミュニケーション機構解明へ期待. 京都大学プレスリリース. 2016-02-17.

Surveys of leucocoprinaceous fungi (Lepiotaceae, Agaricales, Basidiomycota) in the rain-forests of Panama and Brazil revealed several free-living counterparts of fungi cultivated by primitive attine ants (the lower Attini, Formicidae, Hymenoptera), adding to two such collections identified in a survey by Mueller et al (1998). The accumulated evidence supports the hypothesis that perhaps all fungi of lower attine ants have close free-living relatives. Free-living counterparts of ant-cultivated fungi are collected most readily during the early rainy season; in particular these are free-living mushrooms of fungal counterparts that are cultivated as yeasts in gardens of ants in the Cyphomyrmex rimosus group. Free-living and symbiotic fungi of these yeast-cultivating ant species might represent a promising study system to compare the biology of sympatric, conspecific fungi existing outside versus inside the attine symbiosis.

Background Reproductive division of labor in eusocial insects is a striking example of a shared genetic background giving rise to alternative phenotypes, namely queen and worker castes. Queen and worker phenotypes play major roles in the evolution of eusocial insects. Their behavior, morphology and physiology underpin many ecologically relevant colony-level traits, which evolved in parallel in multiple species. Results Using queen and worker transcriptomic data from 16 ant species we tested the hypothesis that conserved sets of genes are involved in ant reproductive division of labor. We further hypothesized that such sets of genes should also be involved in the parallel evolution of other key traits. We applied weighted gene co-expression network analysis, which clusters co-expressed genes into modules, whose expression levels can be summarized by their ‘eigengenes’. Eigengenes of most modules were correlated with phenotypic differentiation between queens and workers. Furthermore, eigengenes of some modules were correlated with repeated evolution of key phenotypes such as complete worker sterility, the number of queens per colony, and even invasiveness. Finally, connectivity and expression levels of genes within the co-expressed network were strongly associated with the strength of selection. Although caste-associated sets of genes evolve faster than non-caste-associated, we found no evidence for queen- or worker-associated co-expressed genes evolving faster than one another. Conclusions These results identify conserved functionally important genomic units that likely serve as building blocks of phenotypic innovation, and allow the remarkable breadth of parallel evolution seen in ants, and possibly other eusocial insects as well.

Background Understanding how dissimilar phenotypes are produced from a similar set of genes is a crucial topic in evolutionary biology. In social insects, both workers and queens develop from diploid eggs, and the extent of phenotypic dimorphism between castes varies between species. Gene expression differences underlying the two female castes, queens and workers, are crucial for social insect evolution. However, it is unknown how much the caste biased proportion of the genome varies in quantity or quality and how this relates to the level of caste polymorphism. This project investigates the differences in gene expression between female castes in Wasmannia auropunctata using deep transcriptional sequencing to understand the role of coding, and potentially non-coding sequences on caste differentiation.