Molecular signatures of plastic phenotypes in two eusocial insect species with simple societies

Amino acid synthesis loss in parasitoid wasps and other hymenopterans

Insects utilize diverse food resources which can affect the evolution of their genomic repertoire, including leading to gene losses in different nutrient pathways. Here, we investigate gene loss in amino acid synthesis pathways, with special attention to hymenopterans and parasitoid wasps. Using comparative genomics, we find that synthesis capability for tryptophan, phenylalanine, tyrosine, and histidine was lost in holometabolous insects prior to hymenopteran divergence, while valine, leucine, and isoleucine were lost in the common ancestor of Hymenoptera. Subsequently, multiple loss events of lysine synthesis occurred independently in the Parasitoida and Aculeata. Experiments in the parasitoid Cotesia chilonis confirm that it has lost the ability to synthesize eight amino acids. Our findings provide insights into amino acid synthesis evolution, and specifically can be used to inform the design of parasitoid artificial diets for pest control.

High-Quality Assemblies for Three Invasive Social Wasps from the Vespula Genus

Social wasps of the genus Vespula have spread to nearly all landmasses worldwide and have become significant pests in their introduced ranges, affecting economies and biodiversity. Comprehensive genome assemblies and annotations for these species are required to develop the next generation of control strategies and monitor existing chemical control. We sequenced and annotated the genomes of the common wasp (Vespula vulgaris), German wasp (Vespula germanica), and the western yellowjacket (Vespula pensylvanica). Our chromosome-level Vespula assemblies each contain 176–179 Mb of total sequence assembled into 25 scaffolds, with 10–200 unanchored scaffolds, and 16,566–18,948 genes. We annotated gene sets relevant to the applied management of invasive wasp populations, including genes associated with spermatogenesis and development, pesticide resistance, olfactory receptors, immunity and venom. These genomes provide evidence for active DNA methylation in Vespidae and tandem duplications of venom genes. Our genomic resources will contribute to the development of next-generation control strategies, and monitoring potential resistance to chemical control.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

DNA methylation suppresses chitin degradation and promotes the wing development by inhibiting Bmara-mediated chitinase expression in the silkworm, Bombyx mori

BACKGROUND

DNA methylation, as an essential epigenetic modification found in mammals and plants, has been implicated to play an important role in insect reproduction. However, the functional role and the regulatory mechanism of DNA methylation during insect organ or tissue development are far from being clear.

RESULTS

Here, we found that DNA methylation inhibitor (5-aza-dC) treatment in newly molted pupae decreased the chitin content of pupal wing discs and adult wings and resulted in wing deformity of Bombyx mori. Transcriptome analysis revealed that the up-regulation of chitinase 10 (BmCHT10) gene might be related to the decrease of chitin content induced by 5-aza-dC treatment. Further, the luciferase activity assays demonstrated that DNA methylation suppressed the promoter activity of BmCHT10 by down-regulating the transcription factor, homeobox protein araucan (Bmara). Electrophoretic mobility shift assay, DNA pull-down and chromatin immunoprecipitation demonstrated that Bmara directly bound to the BmCHT10 promoter. Therefore, DNA methylation is involved in keeping the structural integrity of the silkworm wings from unwanted chitin degradation, as a consequence, it promotes the wing development of B. mori.

CONCLUSIONS

This study reveals that DNA methylation plays an important role in the wing development of B. mori. Our results support that the indirect transcriptional repression of a chitin degradation-related gene BmCHT10 by DNA methylation is necessary to keep the proper wing development in B. mori.

The gene expression network regulating queen brain remodeling after insemination and its parallel use in ants with reproductive workers

Caste differentiation happens early in development to produce gynes as future colony germlines and workers as present colony soma. However, gynes need insemination to become functional queens, a transition that initiates reproductive role differentiation relative to unmated gynes. Here, we analyze the anatomy and transcriptomes of brains during this differentiation process within the reproductive caste of Monomorium pharaonis Insemination terminated brain growth, whereas unmated control gynes continued to increase brain volume. Transcriptomes revealed a specific gene regulatory network (GRN) mediating both brain anatomy changes and behavioral modifications. This reproductive role differentiation GRN hardly overlapped with the gyne-worker caste differentiation GRN, but appears to be also used by distantly related ants where workers became germline individuals after the queen caste was entirely or partially lost. The genes corazonin and neuroparsin A in the anterior neurosecretory cells were overexpressed in individuals with reduced or nonreproductive roles across all four ant species investigated.

Bumblebee Workers Show Differences in Allele-Specific DNA Methylation and Allele-Specific Expression

Allele-specific expression is when one allele of a gene shows higher levels of expression compared with the other allele, in a diploid organism. Recent work has identified allele-specific expression in a number of Hymenopteran species. However, the molecular mechanism which drives this allelic expression bias remains unknown. In mammals, DNA methylation is often associated with genes which show allele-specific expression. DNA methylation systems have been described in species of Hymenoptera, providing a candidate mechanism. Using previously generated RNA-Seq and whole-genome bisulfite sequencing from reproductive and sterile bumblebee (Bombus terrestris) workers, we have identified genome-wide allele-specific expression and allele-specific DNA methylation. The majority of genes displaying allele-specific expression are common between reproductive and sterile workers and the proportion of allele-specific expression bias generally varies between genetically distinct colonies. We have also identified genome-wide allele-specific DNA methylation patterns in both reproductive and sterile workers, with reproductive workers showing significantly more genes with allele-specific methylation. Finally, there is no significant overlap between genes showing allele-specific expression and allele-specific methylation. These results indicate that cis-acting DNA methylation does not directly drive genome-wide allele-specific expression in this species.

Social reprogramming in ants induces longevity-associated glia remodeling

In social insects, workers and queens arise from the same genome but display profound differences in behavior and longevity. In Harpegnathos saltator ants, adult workers can transition to a queen-like state called gamergate, which results in reprogramming of social behavior and life-span extension. Using single-cell RNA sequencing, we compared the distribution of neuronal and glial populations before and after the social transition. We found that the conversion of workers into gamergates resulted in the expansion of neuroprotective ensheathing glia. Brain injury assays revealed that activation of the damage response gene Mmp1 was weaker in old workers, where the relative frequency of ensheathing glia also declined. On the other hand, long-lived gamergates retained a larger fraction of ensheathing glia and the ability to mount a strong Mmp1 response to brain injury into old age. We also observed molecular and cellular changes suggestive of age-associated decline in ensheathing glia in Drosophila.

Brain microRNAs among social and solitary bees

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

DNA methyltransferase 3 participates in behavioral phase change in the migratory locust

DNA methylation plays important roles in the behavioral plasticity of animals. The migratory locust, Locusta migratoria, displays striking density-dependent phenotypic plasticity that can reversely transit between solitarious and gregarious phases. However, the role and the mechanism through which DNA methylation is involved in locust phase change remain unknown. Here, we investigated the expression levels of three DNA methyltransferase genes and their roles in the regulation of locust phase changes. All three Dnmt genes, namely, Dnmt1, Dnmt2 and Dnmt3 showed high expression levels in the brains of gregarious locusts. By contrast, only Dnmt3 transcript rapidly responded to population density changes, decreasing during the isolation of gregarious locusts and steadily increasing upon the crowding of solitarious locusts. Dnmt3 knockdown significantly reduced the phase-related locomotor activity, rather than the attraction index, in gregarious and crowded solitarious locusts. Transcriptome analysis showed that Dnmt3 knockdown upregulated the genes related to metabolism and transporting activity and downregulated those associated with oxidative stress response. The expression level of the phase-core transcriptional factor, hormone receptor HR3, was significantly suppressed in the brain after Dnmt3 knockdown. Moreover, there was significant overlap in the differentially expressed genes between Dnmt3 RNAi and HR3 RNAi data sets, suggesting HR3 may act as key transcriptional factor mediating Dnmt3-controlled gene expression profiles in locust brains. These findings suggest that Dnmt3 transcription is involved in locust behavioral transition, implying the possible roles of DNA methylation in phase-related phenotypic plasticity in locusts.

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

Developmental processes are an important source of phenotypic variation, but the extent to which this variation contributes to evolutionary change is unknown. We used integrative genomic analyses to explore the relationship between developmental and social plasticity in a bee species that can adopt either a social or solitary lifestyle. We find genes regulating this social flexibility also regulate development, and positive selection on these genes is influenced by their function during development. This suggests that developmental plasticity may influence the evolution of sociality. Our additional finding of genetic variants linked to differences in social behavior sheds light on how phenotypic variation derived from development may become encoded into the genome, and thus contribute to evolutionary change.

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

The Roles of DNA Methyltransferases 1 (DNMT1) in Regulating Sexual Dimorphism in the Cotton Mealybug, Phenacoccus solenopsis

The cotton mealybug, Phenacoccus solenopsis, is an invasive pest that can cause massive damage to many host plants of agricultural importance. P. solenopsis is highly polyphagous, and shows extreme sexual dimorphism between males and females. The functions of DNA methyltransferase (DNMT) enzymes in the cotton mealybug have not been well studied. Here, we carried out an investigation of DNMTs in cotton mealybug to study their roles in sexual dimorphism. We found that the cotton mealybug has two copies of PsDnmt1, but Dnmt3 is absent. We then amplified the full-length cDNAs of PsDnmt1A (2,225 bp) and PsDnmt1B (2,862 bp) using rapid amplification cDNA ends (RACE). Quantitative reverse transcriptase PCR shows that both PsDnmt1A and PsDnmt1B are highly expressed in adult males, while the expression of PsDnmt1B is 30-fold higher in gravid females than in virgin females. We knocked down PsDnmt1A and PsDnmt1B with small interfering RNAs (siRNAs), and both genes were successfully down-regulated after 24 h or 72 h in adult females and pupa (t-test, p < 0.05). Down-regulating the expression of these two DNMT genes led to offspring lethality and abnormal body color in adult females. Furthermore, the silencing of PsDnmt1B induced abnormal wing development in emerged adult males. Our results provide evidence that PsDnmt1 plays a crucial role in regulating sexual dimorphism in the cotton mealybug.

Evolutionary dynamics of recent selection on cognitive abilities

Cognitive abilities can vary dramatically among species. The relative importance of social and ecological challenges in shaping cognitive evolution has been the subject of a long-running and recently renewed debate, but little work has sought to understand the selective dynamics underlying the evolution of cognitive abilities. Here, we investigate recent selection related to cognition in the paper wasp Polistes fuscatus-a wasp that has uniquely evolved visual individual recognition abilities. We generate high quality de novo genome assemblies and population genomic resources for multiple species of paper wasps and use a population genomic framework to interrogate the probable mode and tempo of cognitive evolution. Recent, strong, hard selective sweeps in P. fuscatus contain loci annotated with functions in long-term memory formation, mushroom body development, and visual processing, traits which have recently evolved in association with individual recognition. The homologous pathways are not under selection in closely related wasps that lack individual recognition. Indeed, the prevalence of candidate cognition loci within the strongest selective sweeps suggests that the evolution of cognitive abilities has been among the strongest selection pressures in P. fuscatus' recent evolutionary history. Detailed analyses of selective sweeps containing candidate cognition loci reveal multiple cases of hard selective sweeps within the last few thousand years on de novo mutations, mainly in noncoding regions. These data provide unprecedented insight into some of the processes by which cognition evolves.

DNA Methylation Is Correlated with Gene Expression during Diapause Termination of Early Embryonic Development in the Silkworm (Bombyx mori)

DNA modification is a naturally occurring DNA modification in prokaryotic and eukaryotic organisms and is involved in several biological processes. Although genome-wide methylation has been studied in many insects, the understanding of global and genomic DNA methylation during insect early embryonic development, is lacking especially for insect diapause. In this study, we analyzed the relationship between DNA methylomes and transcriptomes in diapause-destined eggs compared to diapause-terminated eggs in the silkworm, Bombyx mori (B. mori). The results revealed that methylation was sparse in this species, as previously reported. Moreover, methylation levels in diapause-terminated eggs (HCl-treated) were 0.05% higher than in non-treated eggs, mainly due to the contribution of CG methylation sites. Methylation tends to occur in the coding sequences and promoter regions, especially at transcription initiation sites and short interspersed elements. Additionally, 364 methylome- and transcriptome-associated genes were identified, which showed significant differences in methylation and expression levels in diapause-destined eggs when compared with diapause-terminated eggs, and 74% of methylome and transcriptome associated genes showed both hypermethylation and elevated expression. Most importantly, Kyoto Encyclopaedia of Genes and Genomes (KEGG) analyses showed that methylation may be positively associated with Bombyx mori embryonic development, by regulating cell differentiation, metabolism, apoptosis pathways and phosphorylation. Through analyzing the G2/M phase-specific E3 ubiquitin-protein ligase (G2E3), we speculate that methylation may affect embryo diapause by regulating the cell cycle in Bombyx mori. These findings will help unravel potential linkages between DNA methylation and gene expression during early insect embryonic development and insect diapause.

A partial genome assembly of the miniature parasitoid wasp, Megaphragma amalphitanum

Body size reduction, also known as miniaturization, is an important evolutionary process that affects a number of physiological and phenotypic traits and helps animals conquer new ecological niches. However, this process is poorly understood at the molecular level. Here, we report genomic and transcriptomic features of arguably the smallest known insect-the parasitoid wasp, Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae). In contrast to expectations, we find that the genome and transcriptome sizes of this parasitoid wasp are comparable to other members of the Chalcidoidea superfamily. Moreover, compared to other chalcid wasps the gene content of M. amalphitanum is remarkably conserved. Intriguingly, we observed significant changes in M. amalphitanum transposable element dynamics over time, in which an initial burst was followed by suppression of activity, possibly due to a recent reinforcement of the genome defense machinery. Overall, while the M. amalphitanum genomic data reveal certain features that may be linked to the unusual biological properties of this organism, miniaturization is not associated with a large decrease in genome complexity.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

Methylation and gene expression differences between reproductive and sterile bumblebee workers

Phenotypic plasticity is the production of multiple phenotypes from a single genome and is notably observed in social insects. Multiple epigenetic mechanisms have been associated with social insect plasticity, with DNA methylation being explored to the greatest extent. DNA methylation is thought to play a role in caste determination in Apis mellifera, and other social insects, but there is limited knowledge on its role in other bee species. In this study, we analyzed whole genome bisulfite sequencing and RNA-seq data sets from head tissue of reproductive and sterile castes of the eusocial bumblebee Bombus terrestris. We found that genome-wide methylation in B. terrestris is similar to other holometabolous insects and does not differ between reproductive castes. We did, however, find differentially methylated genes between castes, which are enriched for multiple biological processes including reproduction. However, we found no relationship between differential methylation and differential gene expression or differential exon usage between castes. Our results also indicate high intercolony variation in methylation. These findings suggest that methylation is associated with caste differences but may serve an alternate function, other than direct caste determination in this species. This study provides the first insights into the nature of a bumblebee caste-specific methylome as well as its interaction with gene expression and caste-specific alternative splicing, providing greater understanding of the role of methylation in phenotypic plasticity within social bee species. Future experimental work is needed to determine the function of methylation and other epigenetic mechanisms in insects.

Sex-specific changes in the aphid DNA methylation landscape

Aphids present an ideal system to study epigenetics as they can produce diverse, but genetically identical, morphs in response to environmental stimuli. Here, using whole genome bisulphite sequencing and transcriptome sequencing of the green peach aphid (Myzus persicae), we present the first detailed analysis of cytosine methylation in an aphid and investigate differences in the methylation and transcriptional landscapes of male and asexual female morphs. We found that methylation primarily occurs in a CG dinucleotide (CpG) context and that exons are highly enriched for methylated CpGs, particularly at the 3' end of genes. Methylation is positively associated with gene expression, and methylated genes are more stably expressed than unmethylated genes. Male and asexual female morphs have distinct methylation profiles. Strikingly, these profiles are divergent between the sex chromosome and the autosomes; autosomal genes are hypomethylated in males compared to asexual females, whereas genes belonging to the sex chromosome, which is haploid in males, are hypermethylated. Overall, we found correlated changes in methylation and gene expression between males and asexual females, and this correlation was particularly strong for genes located on the sex chromosome. Our results suggest that differential methylation of sex-biased genes plays a role in aphid sexual differentiation.

Genome-wide effects of social status on DNA methylation in the brain of a cichlid fish, Astatotilapia burtoni

BACKGROUND

Successful social behavior requires real-time integration of information about the environment, internal physiology, and past experience. The molecular substrates of this integration are poorly understood, but likely modulate neural plasticity and gene regulation. In the cichlid fish species Astatotilapia burtoni, male social status can shift rapidly depending on the environment, causing fast behavioral modifications and a cascade of changes in gene transcription, the brain, and the reproductive system. These changes can be permanent but are also reversible, implying the involvement of a robust but flexible mechanism that regulates plasticity based on internal and external conditions. One candidate mechanism is DNA methylation, which has been linked to social behavior in many species, including A. burtoni. But, the extent of its effects after A. burtoni social change were previously unknown.

RESULTS

We performed the first genome-wide search for DNA methylation patterns associated with social status in the brains of male A. burtoni, identifying hundreds of Differentially Methylated genomic Regions (DMRs) in dominant versus non-dominant fish. Most DMRs were inside genes supporting neural development, synapse function, and other processes relevant to neural plasticity, and DMRs could affect gene expression in multiple ways. DMR genes were more likely to be transcription factors, have a duplicate elsewhere in the genome, have an anti-sense lncRNA, and have more splice variants than other genes. Dozens of genes had multiple DMRs that were often seemingly positioned to regulate specific splice variants.

CONCLUSIONS

Our results revealed genome-wide effects of A. burtoni social status on DNA methylation in the brain and strongly suggest a role for methylation in modulating plasticity across multiple biological levels. They also suggest many novel hypotheses to address in mechanistic follow-up studies, and will be a rich resource for identifying the relationships between behavioral, neural, and transcriptional plasticity in the context of social status.

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.

Patterns of reproductive differentiation and reproductive plasticity in the major evolutionary transition to superorganismality

Major evolutionary transitions in individuality are characterised by the formation of new levels of biological complexity from the cooperation of previously independent lower-level units. The evolution of superorganismality in insects is one such major transition, and is characterised by an extreme division of reproductive labour between ancestrally autonomous units, in the form of queen and worker castes. Here, we discuss the nature of plasticity in the emergence of castes across the major transition to superorganismality in insects. We identify key changes in plasticity which act at different levels of selection: a loss of reproductivity plasticity at the individual level is matched by a gain in plasticity at the colony level. Taking multi-level selection into consideration has important implications for formulating testable hypotheses regarding the nature of plasticity in a major transition from a lower to a higher level of biological complexity.

Bottom-Up Proteomic Analysis of Polypeptide Venom Components of the Giant Ant Dinoponera Quadriceps

Ant species have specialized venom systems developed to sting and inoculate a biological cocktail of organic compounds, including peptide and polypeptide toxins, for the purpose of predation and defense. The genus Dinoponera comprises predatory giant ants that inoculate venom capable of causing long-lasting local pain, involuntary shaking, lymphadenopathy, and cardiac arrhythmias, among other symptoms. To deepen our knowledge about venom composition with regard to protein toxins and their roles in the chemical-ecological relationship and human health, we performed a bottom-up proteomics analysis of the crude venom of the giant ant D. quadriceps, popularly known as the "false" tocandiras. For this purpose, we used two different analytical approaches: (i) gel-based proteomics approach, wherein the crude venom was resolved by denaturing sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and all protein bands were excised for analysis; (ii) solution-based proteomics approach, wherein the crude venom protein components were directly fragmented into tryptic peptides in solution for analysis. The proteomic data that resulted from these two methodologies were compared against a previously annotated transcriptomic database of D. quadriceps, and subsequently, a homology search was performed for all identified transcript products. The gel-based proteomics approach unequivocally identified nine toxins of high molecular mass in the venom, as for example, enzymes [hyaluronidase, phospholipase A1, dipeptidyl peptidase and glucose dehydrogenase/flavin adenine dinucleotide (FAD) quinone] and diverse venom allergens (homologous of the red fire ant Selenopsis invicta) and venom-related proteins (major royal jelly-like). Moreover, the solution-based proteomics revealed and confirmed the presence of several hydrolases, oxidoreductases, proteases, Kunitz-like polypeptides, and the less abundant inhibitor cysteine knot (ICK)-like (knottin) neurotoxins and insect defensin. Our results showed that the major components of the D. quadriceps venom are toxins that are highly likely to damage cell membranes and tissue, to cause neurotoxicity, and to induce allergic reactions, thus, expanding the knowledge about D. quadriceps venom composition and its potential biological effects on prey and victims.

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).

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

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

The effects of the neonicotinoid imidacloprid on gene expression and DNA methylation in the buff-tailed bumblebee Bombus terrestris

Neonicotinoids are effective insecticides used on many important arable and horticultural crops. They are nicotinic acetylcholine receptor agonists which disrupt the function of insect neurons and cause paralysis and death. In addition to direct mortality, there are numerous sublethal effects of low doses of neonicotinoids on bees. We hypothesize that some of these large array of effects could be a consequence of epigenetic changes in bees induced by neonicotinoids. We compared whole methylome (BS-seq) and RNA-seq libraries of the brains of buff-tailed bumblebee Bombus terrestris workers exposed to field-realistic doses of the neonicotinoid imidacloprid to libraries from control workers. We found numerous genes which show differential expression between neonicotinoid-treated bees and control bees, but no differentially methylated cytosines in any context. We found CpG methylation to be focused mainly in exons and associated with highly expressed genes. We discuss the implications of our results for future legislation.

Convergent eusocial evolution is based on a shared reproductive groundplan plus lineage-specific plastic genes

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.

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.

Changes in gene DNA methylation and expression networks accompany caste specialization and age-related physiological changes in a social insect

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, 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.

Extending r/K selection with a maternal risk-management model that classifies animal species into divergent natural selection categories

Reproduction is a defining process of biological systems. Every generation, across all species, breeding females repopulate ecosystems with offspring. r/K selection was the first theory to classify animal species by linking the rates with which breeding females repopulated ecosystems, to the stability of ecosystems. Here, I introduce a species classification scheme that extends the reach of r-K selection and CSR selection by linking breeder investments in offspring quantity, quality, and diversity to specific natural selection pressures. The species classification scheme is predicated on the assumption that high rates of predation favor breeders that invest more in offspring quantity than quality; and that spatiotemporal scarcity favors breeders that investment more in offspring quality than quantity. I present equations that convert the species classification scheme into a maternal risk-management model. Thereafter, using the equations, I classify eighty-seven animal species into the model's natural selection categories. Species of reptiles, fish, and marine invertebrates clustered in the predation selection category. Species of birds and mammals clustered in the scarcity selection category. Several species of apex predators clustered in the weak selection category. Several species of social insects and social mammals clustered in the convergent selection category. In summary, by acknowledging breeding females as the individuals upon which natural selection acts to repopulate ecosystems with offspring, the proposed maternal risk-management model offers a testable, theoretical framework for the field of ecology.

The role of juvenile hormone in regulating reproductive physiology and dominance in Dinoponera quadriceps ants

Unequal reproductive output among members of the same sex (reproductive skew) is a common phenomenon in a wide range of communally breeding animals. In such species, reproductive dominance is often acquired during antagonistic interactions between group members that establish a reproductive hierarchy in which only a few individuals reproduce. Rank-specific syndromes of behavioural and physiological traits characterize such hierarchies, but how antagonistic behavioural interactions translate into stable rank-specific syndromes remains poorly understood. The pleiotropic nature of hormones makes them prime candidates for generating such syndromes as they physiologically integrate environmental (social) information, and often affect reproduction and behaviour simultaneously. Juvenile hormone (JH) is one of several hormones that occupy such a central regulatory role in insects and has been suggested to regulate reproductive hierarchies in a wide range of social insects including ants. Here we use experimental manipulation to investigate the effect of JH levels on reproductive physiology and social dominance in high-ranked workers of the eusocial ant Dinoponera quadriceps, a species that has secondarily reverted to queenless, simple societies. We show that JH regulated reproductive physiology, with ants in which JH levels were experimentally elevated having more regressed ovaries. In contrast, we found no evidence of JH levels affecting dominance in social interactions. This could indicate that JH and ovary development are decoupled from dominance in this species, however only high-ranked workers were investigated. The results therefore confirm that the regulatory role of JH in reproductive physiology in this ant species is in keeping with its highly eusocial ancestors rather than its secondary reversion to simple societies, but more investigation is needed to disentangle the relationships between hormones, behaviour and hierarchies.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.

Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

Changes of gene expression but not cytosine methylation are associated with male parental care reflecting behavioural state, social context and individual flexibility

Behaviour is often a front line response to changing environments. Recent studies show behavioural changes are associated with changes of gene expression; however, these studies have primarily focused on discrete behavioural states. We build on these studies by addressing additional contexts that produce qualitatively similar behavioural changes. We measured levels of gene expression and cytosine methylation, which is hypothesized to regulate the transcriptional architecture of behavioural transitions, within the brain during male parental care of the burying beetle Nicrophorus vespilloides in a factorial design. Male parenting is a suitably plastic behaviour because although male N. vespilloides typically do not provide direct care (i.e. feed offspring) when females are present, levels of feeding by a male equivalent to the female can be induced by removing the female. We examined three different factors: behavioural state (caring versus non-caring), social context (with or without a female mate) and individual flexibility (if a male switched to direct care after his mate was removed). The greatest number of differentially expressed genes were associated with behavioural state, followed by social context and individual flexibility. Cytosine methylation was not associated with changes of gene expression in any of the factors. Our results suggest a hierarchical association between gene expression and the different factors, but that this process is not controlled by cytosine methylation. Our results further suggest that the extent a behaviour is transient plays an underappreciated role in determining its underpinning molecular mechanisms.

Dnmt1 is essential for egg production and embryo viability in the large milkweed bug, Oncopeltus fasciatus

BACKGROUND

The function of cytosine (DNA) methylation in insects remains inconclusive due to a lack of mutant and/or genetic studies.

RESULTS

Here, we provide evidence for the functional role of the maintenance DNA methyltransferase 1 (Dnmt1) in an insect using experimental manipulation. Through RNA interference (RNAi), we successfully posttranscriptionally knocked down Dnmt1 in ovarian tissue of the hemipteran Oncopeltus fasciatus (the large milkweed bug). Individuals depleted for Dnmt1, and subsequently DNA methylation, failed to reproduce. Eggs were inviable and declined in number, and nuclei structure of follicular epithelium was aberrant. Erasure of DNA methylation from gene or transposon element bodies did not reveal a direct causal link to steady-state mRNA levels in somatic cells. These results reveal an important function of Dnmt1 seemingly not contingent on directly controlling gene expression.

CONCLUSIONS

This study provides direct experimental evidence for a functional role of Dnmt1 in egg production and embryo viability and uncovers a trivial role, if any, for DNA methylation in control of gene expression in O. fasciatus.

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.

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

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

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.

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

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

Social wasps as models to study the major evolutionary transition to superorganismality

The major evolutionary transition to superorganismality has taken place several times in the insects. Although there has been much consideration of the ultimate evolutionary explanations for superorganismality, we know relatively little about what proximate mechanisms constrain or promote this major transition. Here, we propose that Vespid wasps represent an understudied, but potentially very useful, model system for studying the mechanisms underpinning superorganismality. We highlight how there is an abundance of behavioural data for many wasp species, confirming their utility in studies of social evolution; however, there is a sparsity of genomic data from which we can test proximate and ultimate hypotheses on this major evolutionary transition.

Do social insects support Haig's kin theory for the evolution of genomic imprinting?

Although numerous imprinted genes have been described in several lineages, the phenomenon of genomic imprinting presents a peculiar evolutionary problem. Several hypotheses have been proposed to explain gene imprinting, the most supported being Haig's kinship theory. This theory explains the observed pattern of imprinting and the resulting phenotypes as a competition for resources between related individuals, but despite its relevance it has not been independently tested. Haig's theory predicts that gene imprinting should be present in eusocial insects in many social scenarios. These lineages are therefore ideal for testing both the theory's predictions and the mechanism of gene imprinting. Here we review the behavioral evidence of genomic imprinting in eusocial insects, the evidence of a mechanism for genomic imprinting and finally we evaluate recent results showing parent of origin allele specific expression in honeybees in the light of Haig's theory.

Social regulation of insulin signaling and the evolution of eusociality in ants

Queens and workers of eusocial Hymenoptera are considered homologous to the reproductive and brood care phases of an ancestral subsocial life cycle. However, the molecular mechanisms underlying the evolution of reproductive division of labor remain obscure. Using a brain transcriptomics screen, we identified a single gene, insulin-like peptide 2 (ilp2), which is always up-regulated in ant reproductives, likely because they are better nourished than their nonreproductive nestmates. In clonal raider ants (Ooceraea biroi), larval signals inhibit adult reproduction by suppressing ilp2, thus producing a colony reproductive cycle reminiscent of ancestral subsociality. However, increasing ILP2 peptide levels overrides larval suppression, thereby breaking the colony cycle and inducing a stable division of labor. These findings suggest a simple model for the origin of ant eusociality via nutritionally determined reproductive asymmetries potentially amplified by larval signals.

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.

Tensor decomposition-based and principal-component-analysis-based unsupervised feature extraction applied to the gene expression and methylation profiles in the brains of social insects with multiple castes

Background

Even though coexistence of multiple phenotypes sharing the same genomic background is interesting, it remains incompletely understood. Epigenomic profiles may represent key factors, with unknown contributions to the development of multiple phenotypes, and social-insect castes are a good model for elucidation of the underlying mechanisms. Nonetheless, previous studies have failed to identify genes associated with aberrant gene expression and methylation profiles because of the lack of suitable methodology that can address this problem properly.

Methods

A recently proposed principal component analysis (PCA)-based and tensor decomposition (TD)-based unsupervised feature extraction (FE) can solve this problem because these two approaches can deal with gene expression and methylation profiles even when a small number of samples is available.

Results

PCA-based and TD-based unsupervised FE methods were applied to the analysis of gene expression and methylation profiles in the brains of two social insects, Polistes canadensis and Dinoponera quadriceps. Genes associated with differential expression and methylation between castes were identified, and analysis of enrichment of Gene Ontology terms confirmed reliability of the obtained sets of genes from the biological standpoint.

Conclusions

Biologically relevant genes, shown to be associated with significant differential gene expression and methylation between castes, were identified here for the first time. The identification of these genes may help understand the mechanisms underlying epigenetic control of development of multiple phenotypes under the same genomic conditions.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2068-7) contains supplementary material, which is available to authorized users.

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

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

Comparison of whole-genome bisulfite sequencing library preparation strategies identifies sources of biases affecting DNA methylation data

Background

Whole-genome bisulfite sequencing (WGBS) is becoming an increasingly accessible technique, used widely for both fundamental and disease-oriented research. Library preparation methods benefit from a variety of available kits, polymerases and bisulfite conversion protocols. Although some steps in the procedure, such as PCR amplification, are known to introduce biases, a systematic evaluation of biases in WGBS strategies is missing.

Results

We perform a comparative analysis of several commonly used pre- and post-bisulfite WGBS library preparation protocols for their performance and quality of sequencing outputs. Our results show that bisulfite conversion per se is the main trigger of pronounced sequencing biases, and PCR amplification builds on these underlying artefacts. The majority of standard library preparation methods yield a significantly biased sequence output and overestimate global methylation. Importantly, both absolute and relative methylation levels at specific genomic regions vary substantially between methods, with clear implications for DNA methylation studies.

Conclusions

We show that amplification-free library preparation is the least biased approach for WGBS. In protocols with amplification, the choice of bisulfite conversion protocol or polymerase can significantly minimize artefacts. To aid with the quality assessment of existing WGBS datasets, we have integrated a bias diagnostic tool in the Bismark package and offer several approaches for consideration during the preparation and analysis of WGBS datasets.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1408-2) contains supplementary material, which is available to authorized users.

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.

A molecular concept of caste in insect societies

The term 'caste' is used to describe the division of reproductive labour that defines eusocial insect societies. The definition of 'caste' has been debated over the last 50 years, specifically with respect to the simplest insect societies; this raises the question of whether a simple categorisation of social behaviour by reproductive state alone is helpful. Gene-level analyses of behaviours of individuals in hymenopteran social insect societies now provide a new empirical base-line for defining caste and understanding the evolution and maintenance of a reproductive division of labour. We review this literature to identify a set of potential molecular signatures that, combined with behavioural, morphological and physiological data, help define caste more precisely; these signatures vary with the type of society, and are likely to be influenced by ecology, life-history, and stage in the colony cycle. We conclude that genomic approaches provide us with additional ways to help quantify and categorise caste, and behaviour in general.