Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

Leveraging technological innovations to investigate evolutionary transitions to eusociality

The study of the major transition to eusociality presents several challenges to researchers, largely resulting from the importance of complex behavioral phenotypes and the shift from individual to group level selection. These challenges are being met with corresponding technological improvements. Advances in resource development for non-model taxa, behavioral tracking, nucleic acid sequencing, and reverse genetics are facilitating studies of hypotheses that were previously intractable. These innovations are resulting in the development of new model systems tailored to the exploration of specific behavioral phenotypes and the querying of underlying molecular mechanisms that drive eusocial behaviors. Here, we present a brief overview of how methodological innovations are advancing our understanding of the evolution of eusociality.

Rate variation in the evolution of non-coding DNA associated with social evolution in bees

The evolutionary origins of eusociality represent increases in complexity from individual to caste-based, group reproduction. These behavioural transitions have been hypothesized to go hand in hand with an increased ability to regulate when and where genes are expressed. Bees have convergently evolved eusociality up to five times, providing a framework to test this hypothesis. To examine potential links between putative gene regulatory elements and social evolution, we compare alignable, non-coding sequences in 11 diverse bee species, encompassing three independent origins of reproductive division of labour and two elaborations of eusocial complexity. We find that rates of evolution in a number of non-coding sequences correlate with key social transitions in bees. Interestingly, while we find little evidence for convergent rate changes associated with independent origins of social behaviour, a number of molecular pathways exhibit convergent rate changes in conjunction with subsequent elaborations of social organization. We also present evidence that many novel non-coding regions may have been recruited alongside the origin of sociality in corbiculate bees; these loci could represent gene regulatory elements associated with division of labour within this group. Thus, our findings are consistent with the hypothesis that gene regulatory innovations are associated with the evolution of eusociality and illustrate how a thorough examination of both coding and non-coding sequence can provide a more complete understanding of the molecular mechanisms underlying behavioural evolution. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Core transcriptional signatures of phase change in the migratory locust

Phenotypic plasticity plays fundamental roles in successful adaptation of animals in response to environmental variations. Here, to reveal the transcriptome reprogramming in locust phase change, a typical phenotypic plasticity, we conducted a comprehensive analysis of multiple phase-related transcriptomic datasets of the migratory locust. We defined PhaseCore genes according to their contribution to phase differentiation by the adjustment for confounding principal components analysis algorithm (AC-PCA). Compared with other genes, PhaseCore genes predicted phase status with over 87.5% accuracy and displayed more unique gene attributes including the faster evolution rate, higher CpG content and higher specific expression level. Then, we identified 20 transcription factors (TFs) named PhaseCoreTF genes that are associated with the regulation of PhaseCore genes. Finally, we experimentally verified the regulatory roles of three representative TFs (Hr4, Hr46, and grh) in phase change by RNAi. Our findings revealed that core transcriptional signatures are involved in the global regulation of locust phase changes, suggesting a potential common mechanism underlying phenotypic plasticity in insects. The expression and network data are accessible in an online resource called LocustMine (http://www.locustmine.org:8080/locustmine).

A Maternal Effect on Queen Production in Honeybees

Influences from the mother on offspring phenotype, known as maternal effects, are an important cause of adaptive phenotypic plasticity [1, 2]. Eusocial insects show dramatic phenotypic plasticity with morphologically distinct reproductive (queen) and worker castes [3, 4]. The dominant paradigm for honeybees (Apis mellifera) is that castes are environmentally rather than genetically determined, with the environment and diet of young larvae causing caste differentiation [5-9]. A role for maternal effects has not been considered, but here we show that egg size also influences queen development. Queens laid significantly bigger eggs in the larger queen cells than in the worker cells. Eggs laid in queen cells (QE), laid in worker cells (WE), and 2-day old larvae from worker cells (2L) were transferred to artificial queen cells to be reared as queens in a standardized environment. Newly emerged adult queens from QE were heavier than those from the other two groups and had more ovarioles, indicating a consequence of egg size for adult queen morphology. Gene expression analyses identified several significantly differentially expressed genes between newly emerged queens from QE and those from the other groups. These included a disproportionate number of genes involved in hormonal signaling, body development, and immune pathways, which are key traits differing between queens and workers. That egg size influences emerging queen morphology and physiology and that queens lay larger eggs in queen cells demonstrate both a maternal effect on the expression of the queen phenotype and a more active role for the queen in gyne production than has been realized previously.

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.

The first draft genomes of the ant Formica exsecta, and its Wolbachia endosymbiont reveal extensive gene transfer from endosymbiont to host

Background

Adapting to changes in the environment is the foundation of species survival, and is usually thought to be a gradual process. However, transposable elements (TEs), epigenetic modifications, and/or genetic material acquired from other organisms by means of horizontal gene transfer (HGTs), can also lead to novel adaptive traits. Social insects form dense societies, which attract and maintain extra- and intracellular accessory inhabitants, which may facilitate gene transfer between species. The wood ant Formica exsecta (Formicidae; Hymenoptera), is a common ant species throughout the Palearctic region. The species is a well-established model for studies of ecological characteristics and evolutionary conflict.

Results

In this study, we sequenced and assembled draft genomes for F. exsecta and its endosymbiont Wolbachia. The F. exsecta draft genome is 277.7 Mb long; we identify 13,767 protein coding genes, for which we provide gene ontology and protein domain annotations. This is also the first report of a Wolbachia genome from ants, and provides insights into the phylogenetic position of this endosymbiont. We also identified multiple horizontal gene transfer events (HGTs) from Wolbachia to F. exsecta. Some of these HGTs have also occurred in parallel in multiple other insect genomes, highlighting the extent of HGTs in eukaryotes.

Conclusion

We present the first draft genome of ant F. exsecta, and its endosymbiont Wolbachia (wFex), and show considerable rates of gene transfer from the symbiont to the host. We expect that especially the F. exsecta genome will be valuable resource in further exploration of the molecular basis of the evolution of social organization.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5665-6) contains supplementary material, which is available to authorized users.

Caste-specific microRNA expression in termites: insights into soldier differentiation

Eusocial insects have polyphenic caste systems in which each caste exhibits characteristic morphology and behaviour. In insects, caste systems arose independently in different lineages, such as Isoptera and Hymenoptera. Although partial molecular mechanisms for the development of eusociality in termites have been clarified by the functional analysis of genes and hormones, the contribution of microRNAs (miRNAs) to caste differentiation is unknown. To understand the role of miRNAs in termite caste polyphenism, we performed small RNA sequencing in a subterranean termite (Reticulitermes speratus) and identified the miRNAs that were specifically expressed in the soldier and worker castes. Of the 550 miRNAs annotated in the R. speratus genome, 74 were conserved in insects and 174 were conserved in other termite species. We found that eight miRNAs (mir-1, mir-125, mir-133, mir-2765, mir-87a and three termite-specific miRNAs) are differentially expressed (DE) in soldiers and workers of R. speratus. This differential expression was experimentally verified for five miRNAs by real-time quantitative PCR. Further, four of the eight DE miRNAs in soldier and worker termite castes were also differentially expressed in hymenopteran castes. The finding that Isoptera and Hymenoptera shared several DE miRNAs amongst castes suggests that these miRNAs evolved independently in these phylogenetically distinct lineages.

Evolutionary Potential of Cis-Regulatory Mutations to Cause Rapid Changes in Transcription Factor Binding

Transcriptional regulation is crucial for all biological processes and well investigated at the molecular level for a wide range of organisms. However, it is quite unclear how innovations, such as the activity of a novel regulatory element, evolve. In the case of transcription factor (TF) binding, both a novel TF and a novel-binding site would need to evolve concertedly. Since promiscuous functions have recently been identified as important intermediate steps in creating novel specific functions in many areas such as enzyme evolution and protein-protein interactions, we ask here how promiscuous binding of TFs to TF-binding sites (TFBSs) affects the robustness and evolvability of this tightly regulated system. Specifically, we investigate the binding behavior of several hundred TFs from different species at unprecedented breadth. Our results illustrate multiple aspects of TF-binding interactions, ranging from correlations between the strength of the interaction bond and specificity, to preferences regarding TFBS nucleotide composition in relation to both domains and binding specificity. We identified a subset of high A/T binding motifs. Motifs in this subset had many functionally neutral one-error mutants, and were bound by multiple different binding domains. Our results indicate that, especially for some TF-TFBS associations, low binding specificity confers high degrees of evolvability, that is that few mutations facilitate rapid changes in transcriptional regulation, in particular for large and old TF families. In this study we identify binding motifs exhibiting behavior indicating high evolutionary potential for innovations in transcriptional regulation.

Universality of the DNA methylation codes in Eucaryotes

Genetics and epigenetics are tightly linked heritable information classes. Question arises if epigenetics provides just a set of environment dependent instructions, or whether it is integral part of an inheritance system. We argued that in the latter case the epigenetic code should share the universality quality of the genetic code. We focused on DNA methylation. Since availability of DNA methylation data is biased towards model organisms we developed a method that uses kernel density estimations of CpG observed/expected ratios to infer DNA methylation types in any genome. We show here that our method allows for robust prediction of mosaic and full gene body methylation with a PPV of 1 and 0.87, respectively. We used this prediction to complement experimental data, and applied hierarchical clustering to identify methylation types in ~150 eucaryotic species covering different body plans, reproduction types and living conditions. Our analysis indicates that there are only four gene body methylation types. These types do not follow phylogeny (i.e. phylogenetically distant clades can have identical methylation types) but they are consistent within clades. We conclude that the gene body DNA methylation codes have universality similar to the universality of the genetic code and should consequently be considered as part of the inheritance system.

Genomic Signals of Adaptation towards Mutualism and Sociality in Two Ambrosia Beetle Complexes

Mutualistic symbiosis and eusociality have developed through gradual evolutionary processes at different times in specific lineages. Like some species of termites and ants, ambrosia beetles have independently evolved a mutualistic nutritional symbiosis with fungi, which has been associated with the evolution of complex social behaviors in some members of this group. We sequenced the transcriptomes of two ambrosia complexes (Euwallacea sp. near fornicatus⁻Fusarium euwallaceae and Xyleborus glabratus⁻Raffaelea lauricola) to find evolutionary signatures associated with mutualism and behavior evolution. We identified signatures of positive selection in genes related to nutrient homeostasis; regulation of gene expression; development and function of the nervous system, which may be involved in diet specialization; behavioral changes; and social evolution in this lineage. Finally, we found convergent changes in evolutionary rates of proteins across lineages with phylogenetically independent origins of sociality and mutualism, suggesting a constrained evolution of conserved genes in social species, and an evolutionary rate acceleration related to changes in selective pressures in mutualistic lineages.

Differential Gene Expression in Red Imported Fire Ant (Solenopsis invicta) (Hymenoptera: Formicidae) Larval and Pupal Stages

Solenopsis invicta Buren is an invasive ant species that has been introduced to multiple continents. One such area, the southern United States, has a history of multiple control projects using chemical pesticides over varying ranges, often resulting in non-target effects across trophic levels. With the advent of next generation sequencing and RNAi technology, novel investigations and new control methods are possible. A robust genome-guided transcriptome assembly was used to investigate gene expression differences between S. invicta larvae and pupae. These life stages differ in many physiological processes; of special importance is the vital role of S. invicta larvae as the colonies' "communal gut". Differentially expressed transcripts were identified related to many important physiological processes, including digestion, development, cell regulation and hormone signaling. This dataset provides essential developmental knowledge that reveals the dramatic changes in gene expression associated with social insect life stage roles, and can be leveraged using RNAi to develop effective control methods.

Social environment affects the transcriptomic response to bacteria in ant queens

Social insects have evolved enormous capacities to collectively build nests and defend their colonies against both predators and pathogens. The latter is achieved by a combination of individual immune responses and sophisticated collective behavioral and organizational disease defenses, that is, social immunity. We investigated how the presence or absence of these social defense lines affects individual-level immunity in ant queens after bacterial infection. To this end, we injected queens of the ant Linepithema humile with a mix of gram+ and gram- bacteria or a control solution, reared them either with workers or alone and analyzed their gene expression patterns at 2, 4, 8, and 12 hr post-injection, using RNA-seq. This allowed us to test for the effect of bacterial infection, social context, as well as the interaction between the two over the course of infection and raising of an immune response. We found that social isolation per se affected queen gene expression for metabolism genes, but not for immune genes. When infected, queens reared with and without workers up-regulated similar numbers of innate immune genes revealing activation of Toll and Imd signaling pathways and melanization. Interestingly, however, they mostly regulated different genes along the pathways and showed a different pattern of overall gene up-regulation or down-regulation. Hence, we can conclude that the absence of workers does not compromise the onset of an individual immune response by the queens, but that the social environment impacts the route of the individual innate immune responses.

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

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

Signatures of DNA Methylation across Insects Suggest Reduced DNA Methylation Levels in Holometabola

It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.

Conserved Genes Underlie Phenotypic Plasticity in an Incipiently Social Bee

Despite a strong history of theoretical work on the mechanisms of social evolution, relatively little is known of the molecular genetic changes that accompany transitions from solitary to eusocial forms. Here, we provide the first genome of an incipiently social bee that shows both solitary and social colony organization in sympatry, the Australian carpenter bee Ceratina australensis. Through comparative analysis, we provide support for the role of conserved genes and cis-regulation of gene expression in the phenotypic plasticity observed in nest-sharing, a rudimentary form of sociality. Additionally, we find that these conserved genes are associated with caste differences in advanced eusocial species, suggesting these types of mechanisms could pave the molecular pathway from solitary to eusocial living. Genes associated with social nesting in this species show signatures of being deeply conserved, in contrast to previous studies in other bees showing novel and faster-evolving genes are associated with derived sociality. Our data provide support for the idea that the earliest social transitions are driven by changes in gene regulation of deeply conserved genes.

Clonal raider ant brain transcriptomics identifies candidate molecular mechanisms for reproductive division of labor

BACKGROUND

Division of labor between reproductive queens and workers that perform brood care is a hallmark of insect societies. However, studies of the molecular basis of this fundamental dichotomy are limited by the fact that the caste of an individual cannot typically be experimentally manipulated at the adult stage. Here we take advantage of the unique biology of the clonal raider ant, Ooceraea biroi, to study brain gene expression dynamics during experimentally induced transitions between reproductive and brood care behavior.

RESULTS

Introducing larvae that inhibit reproduction and induce brood care behavior causes much faster changes in adult gene expression than removing larvae. In addition, the general patterns of gene expression differ depending on whether ants transition from reproduction to brood care or vice versa, indicating that gene expression changes between phases are cyclic rather than pendular. Finally, we identify genes that could play upstream roles in regulating reproduction and behavior because they show large and early expression changes in one or both transitions.

CONCLUSIONS

Our analyses reveal that the nature and timing of gene expression changes differ substantially depending on the direction of the transition, and identify a suite of promising candidate molecular regulators of reproductive division of labor that can now be characterized further in both social and solitary animal models. This study contributes to understanding the molecular regulation of reproduction and behavior, as well as the organization and evolution of insect societies.

The evolution of sociality in termites from cockroaches: A taxonomic and phylogenetic perspective

Despite multiple studies and advances, sociality still puzzles evolutionary biologists in numerous ways, which might be partly addressed with the advent of sociogenomics. In insects, the majority of sociogenomic studies deal with Hymenoptera, one of the two groups that evolved eusociality with termites. But, to fully grasp the evolution of sociality, studies must obviously not restrict to eusocial lineages. Multiple kinds of social system transitions have been recorded and they all bring complementary insights. For instance, cockroaches, the closest relatives to termites, display a wide range of social interactions and evolved convergently subsocial behaviors (i.e., brood care). In this context, we emphasize the need for natural history, taxonomic, and phylogenetic studies. Natural history studies provide the foundations on which building hypotheses, whereas taxonomy provides the taxa to sample to test these hypotheses, and phylogenetics brings the historical framework necessary to test evolutionary scenarios of sociality evolution.

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

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

Levels of DNA cytosine methylation in the Drosophila genome

Insects provide an accessible system to study the contribution of DNA methylation to complex epigenetic phenotypes created to regulate gene expression, chromatin states, imprinting and dosage compensation. The members of genus Drosophila have been used as a model system to study aspects of biology like development, behaviour and genetics. Despite the popularity of Drosophila melanogaster as a genetic and epigenetic model organism, DNA methylation studies are limited due to low levels of genomic 5-methylcytosine. Our study employs a sensitive liquid chromatography-mass spectrometry (LCMS) based method to quantify the levels of 5-methylcytosine from the genomic DNA in different members of the genus Drosophila. Our results reveal that, despite being phylogenetically related, there is a marked variation in the levels of 5-methylcytosine between the genomes of the members of genus Drosophila. Also, there is a change in the genomic levels of 5-methylcytosine through each life cycle stage of holometabolous development in D. melanogaster.

Oxytocin-like signaling in ants influences metabolic gene expression and locomotor activity

Ants are emerging model systems to study cellular signaling because distinct castes possess different physiologic phenotypes within the same colony. Here we studied the functionality of inotocin signaling, an insect ortholog of mammalian oxytocin (OT), which was recently discovered in ants. In Lasius ants, we determined that specialization within the colony, seasonal factors, and physiologic conditions down-regulated the expression of the OT-like signaling system. Given this natural variation, we interrogated its function using RNAi knockdowns. Next-generation RNA sequencing of OT-like precursor knock-down ants highlighted its role in the regulation of genes involved in metabolism. Knock-down ants exhibited higher walking activity and increased self-grooming in the brood chamber. We propose that OT-like signaling in ants is important for regulating metabolic processes and locomotion.-Liutkevičiūtė, Z., Gil-Mansilla, E., Eder, T., Casillas-Pérez, B., Di Giglio, M. G., Muratspahić, E., Grebien, F., Rattei, T., Muttenthaler, M., Cremer, S., Gruber, C. W. Oxytocin-like signaling in ants influences metabolic gene expression and locomotor activity.

FoxP in bees: A comparative study on the developmental and adult expression pattern in three bee species considering isoforms and circuitry

Mutations in the transcription factors FOXP1, FOXP2, and FOXP4 affect human cognition, including language. The FoxP gene locus is evolutionarily ancient and highly conserved in its DNA-binding domain. In Drosophila melanogaster FoxP has been implicated in courtship behavior, decision making, and specific types of motor-learning. Because honeybees (Apis mellifera, Am) excel at navigation and symbolic dance communication, they are a particularly suitable insect species to investigate a potential link between neural FoxP expression and cognition. We characterized two AmFoxP isoforms and mapped their expression in the brain during development and in adult foragers. Using a custom-made antiserum and in situ hybridization, we describe 11 AmFoxP expressing neuron populations. FoxP was expressed in equivalent patterns in two other representatives of Apidae; a closely related dwarf bee and a bumblebee species. Neural tracing revealed that the largest FoxP expressing neuron cluster in honeybees projects into a posterior tract that connects the optic lobe to the posterior lateral protocerebrum, predicting a function in visual processing. Our data provide an entry point for future experiments assessing the function of FoxP in eusocial Hymenoptera.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Hemimetabolous genomes reveal molecular basis of termite eusociality

Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.

Developmental DNA methyltransferase expression in the fire ant Solenopsis invicta

DNA methylation is accomplished in animals by 2 classes of enzymes known as DNA methyltransferases, DNMT3 and DNMT1, which perform de novo methylation and maintenance methylation, respectively. Several studies of hymenopteran eusocial insects suggest that DNA methylation is capable of influencing developmental plasticity. However, fundamental questions remain about the patterning of DNA methylation during the course of insect development. In this study, we performed quantitative real-time PCR (qPCR) on transcripts from the single-copy orthologs of DNMT1 and DNMT3 in the red imported fire ant, Solenopsis invicta. In particular, we assessed the expression of S. invicta Dnmt1 and Dnmt3 mRNA during 7 stages of worker development, among behaviorally distinct adults, and among male and female gonads. Dnmt3 was most highly expressed during embryonic development, whereas Dnmt1 was similarly expressed throughout the course of development. Moreover, Dnmt1 and Dnmt3 were highly expressed in testes and ovaries. Neither Dnmt was significantly differentially expressed among heads of behaviorally distinct adult castes. Our results support the hypothesis that extensive patterning of DNA methylation occurs during gametogenesis and embryogenesis in the insect order Hymenoptera.

Comparative analyses of co-evolving host-parasite associations reveal unique gene expression patterns underlying slavemaker raiding and host defensive phenotypes

The transition to parasitism is a drastic shift in lifestyle, involving rapid changes in gene structure, function, and expression. After the establishment of antagonistic relationships, parasites and hosts co-evolve through reciprocal adaptations, often resulting in evolutionary arms-races. Repeated evolution of social parasitism and slavery among Temnothorax ants allows us to examine those gene expression patterns that characterize slavemaker raiding and reciprocal host defensive phenotypes. Previous behavioural studies have established that raiding strategies between Temnothorax slavemakers diverge, while host defense portfolios shift similarly under parasite pressure. We are the first to confirm this at the molecular level, revealing that slavemaking species exhibit a wider variety of genes with species-specific patterns of expression within their raiding phenotypes, whereas expression similarity is commonly found during the non-raiding phenotype. Host species response to slavemaker aggression, however, is indicated by strong changes in the expression of a relatively few number genes. Additionally, the expression of individual genes such as Acyl-CoA-Delta(11) desaturase and Trypsin-7 is strongly associated with the raiding phenotype of all three slavemaking species. Here, we provide novel insight into the gene expression patterns associated with raiding and nest defense behavior in Temnothorax ants, suggesting lineage-specific evolutionary patterns among both slavemakers and hosts.

The genomic and transcriptomic analyses of serine proteases and their homologs in an endoparasitoid, Pteromalus puparum

In insects, serine proteases (SPs) and serine protease homologs (SPHs) constitute a large family of proteins involved in multiple physiological processes such as digestion, development, and immunity. Here we identified 145 SPs and 38 SPHs in the genome of an endoparasitoid, Pteromalus puparum. Gene duplication and tandem repeats were observed in this large SPs/SPHs family. We then analyzed the expression profiles of SP/SPH genes in response to different microbial infections (Gram-positive bacterium Micrococcus luteus, Gram-negative bacterium Escherichia coli, and entomopathogenic fungus Beauveria bassiana), as well as in different developmental stages and tissues. Some SPs/SPHs also displayed distinct expression patterns in venom gland, suggesting their specific physiological functions as venom proteins. Our finding lays groundwork for further research of SPs and SPHs expressed in the venom glands.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

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

Deconstructing Superorganisms and Societies to Address Big Questions in Biology

Social insect societies are long-standing models for understanding social behaviour and evolution. Unlike other advanced biological societies (such as the multicellular body), the component parts of social insect societies can be easily deconstructed and manipulated. Recent methodological and theoretical innovations have exploited this trait to address an expanded range of biological questions. We illustrate the broadening range of biological insight coming from social insect biology with four examples. These new frontiers promote open-minded, interdisciplinary exploration of one of the richest and most complex of biological phenomena: sociality.

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

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

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0136-2) contains supplementary material, which is available to authorized users.

Sociomics: Using Omic Approaches to Understand Social Evolution

All of life is social, from genes cooperating to form organisms, to animals cooperating to form societies. Omic approaches offer exceptional opportunities to solve major outstanding problems in the study of how sociality evolves. First, omics can be used to clarify the extent and form of sociality in natural populations. This is especially useful in species where it is difficult to study social traits in natural populations, such as bacteria and other microbes. Second, omics can be used to examine the consequences of sociality for genome evolution and gene expression. This is especially useful in cases where there is clear variation in the level of sociality, such as the social insects. Major tasks for the future are to apply these approaches to a wider range of non-model organisms, and to move from exploratory analyses to the testing of evolutionary theory.

Genomic exaptation enables Lasius niger adaptation to urban environments

Background

The world is rapidly urbanizing, and only a subset of species are able to succeed in stressful city environments. Efficient genome-enabled stress response appears to be a likely prerequisite for urban adaptation. Despite the important role ants play in the ecosytem, only the genomes of ~13 have been sequenced so far. Here, we present the draft genome assembly of the black garden ant Lasius niger – the most successful urban inhabitant of all ants – and we compare it with the genomes of other ant species, including the closely related Camponotus floridanus.

Results

Sequences from 272 M Illumina reads were assembled into 41,406 contigs with total length of 245 MB, and N50 of 16,382 bp, similar to other ant genome assemblies enabling comparative genomic analysis. Remarkably, the predicted proteome of L. niger is significantly enriched relative to other ant genomes in terms of abundance of domains involved in nucleic acid binding, DNA repair, and nucleotidyl transferase activity, reflecting transposable element proliferation and a likely genomic response. With respect to environmental stress, we note a proliferation of various detoxification genes, including glutatione-S-transferases and those in the cytochrome P450 families. Notably, the CYP9 family is highly expanded with 19 complete and 21 nearly complete members - over twice as many compared to other ants. This family exhibits the signatures of strong directional selection, with eleven positively selected positions in ligand-binding pockets of enzymes. Gene family contraction was detected for several components of the olfactory system, accompanied by instances of both directional selection and relaxation.

Conclusions

Our results suggest that the success of L. niger in urbanized areas may be the result of fortuitous coincidence of several factors, including the expansion of the CYP9 cytochrome family due to coevolution with parasitic fungi, the diversification of DNA repair systems as an answer to proliferation of retroelements, and the reduction of olfactory system and behavioral preadaptations from non-territorial subdominant life strategies found in natural environments. Diversification of cytochromes and DNA repair systems along with reduced odorant communication are the basis of L. niger pollutant resistance and polyphagy, while non-territorial and mobilization strategies allows more efficient exploitation of large but patchy food sources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0867-x) contains supplementary material, which is available to authorized users.

Molecular Evolution of Insect Sociality: An Eco-Evo-Devo Perspective

The evolution of eusociality is a perennial issue in evolutionary biology, and genomic advances have fueled steadily growing interest in the genetic changes underlying social evolution. Along with a recent flurry of research on comparative and evolutionary genomics in different eusocial insect groups (bees, ants, wasps, and termites), several mechanistic explanations have emerged to describe the molecular evolution of eusociality from solitary behavior. These include solitary physiological ground plans, genetic toolkits of deeply conserved genes, evolutionary changes in protein-coding genes, cis regulation, and the structure of gene networks, epigenetics, and novel genes. Despite this proliferation of ideas, there has been little synthesis, even though these ideas are not mutually exclusive and may in fact be complementary. We review available data on molecular evolution of insect sociality and highlight key biotic and abiotic factors influencing social insect genomes. We then suggest both phylogenetic and ecological evolutionary developmental biology (eco-evo-devo) perspectives for a more synthetic view of molecular evolution in insect societies.

Gene Expression Dynamics in Major Endocrine Regulatory Pathways along the Transition from Solitary to Social Life in a Bumblebee, Bombus terrestris

Understanding the social evolution leading to insect eusociality requires, among other, a detailed insight into endocrine regulatory mechanisms that have been co-opted from solitary ancestors to play new roles in the complex life histories of eusocial species. Bumblebees represent well-suited models of a relatively primitive social organization standing on the mid-way to highly advanced eusociality and their queens undergo both, a solitary and a social phase, separated by winter diapause. In the present paper, we characterize the gene expression levels of major endocrine regulatory pathways across tissues, sexes, and life-stages of the buff-tailed bumblebee, Bombus terrestris, with special emphasis on critical stages of the queen's transition from solitary to social life. We focused on fundamental genes of three pathways: (1) Forkhead box protein O and insulin/insulin-like signaling, (2) Juvenile hormone (JH) signaling, and (3) Adipokinetic hormone signaling. Virgin queens were distinguished by higher expression of forkhead box protein O and downregulated insulin-like peptides and JH signaling, indicated by low expression of methyl farnesoate epoxidase (MFE) and transcription factor Krüppel homolog 1 (Kr-h1). Diapausing queens showed the expected downregulation of JH signaling in terms of low MFE and vitellogenin (Vg) expressions, but an unexpectedly high expression of Kr-h1. By contrast, reproducing queens revealed an upregulation of MFE and Vg together with insulin signaling. Surprisingly, the insulin growth factor 1 (IGF-1) turned out to be a queen-specific hormone. Workers exhibited an expression pattern of MFE and Vg similar to that of reproducing queens. Males were characterized by high Kr-h1 expression and low Vg level. The tissue comparison unveiled an unexpected resemblance between the fat body and hypopharyngeal glands across all investigated genes, sexes, and life stages.

Hoarding and horizontal transfer led to an expanded gene and intron repertoire in the plastid genome of the diatom, Toxarium undulatum (Bacillariophyta)

Although the plastid genomes of diatoms maintain a conserved architecture and core gene set, considerable variation about this core theme exists and can be traced to several different processes. Gene duplication, pseudogenization, and loss, as well as intracellular transfer of genes to the nuclear genome, have all contributed to variation in gene content among diatom species. In addition, some noncoding sequences have highly restricted phylogenetic distributions that suggest a recent foreign origin. We sequenced the plastid genome of the marine diatom, Toxarium undulatum, and found that the genome contains three genes (chlB, chlL, and chlN) involved in light-independent chlorophyll a biosynthesis that were not previously known from diatoms. Phylogenetic and syntenic data suggest that these genes were differentially retained in this one lineage as they were repeatedly lost from most other diatoms. Unique among diatoms and other heterokont algae sequenced so far, the genome also contains a large group II intron within an otherwise intact psaA gene. Although the intron is most similar to one in the plastid-encoded psaA gene of some green algae, high sequence divergence between the diatom and green algal introns rules out recent shared ancestry. We conclude that the psaA intron was likely introduced into the plastid genome of T. undulatum, or some earlier ancestor, by horizontal transfer from an unknown donor. This genome further highlights the myriad processes driving variation in gene and intron content in the plastid genomes of diatoms, one of the world's foremost primary producers.