Methylation and worker reproduction in the bumble-bee (Bombus terrestris)

A role for DNA methylation in bumblebee morphogenesis hints at female-specific developmental erasure

Epigenetic mechanisms, such as DNA methylation, are crucial factors in animal development. In some mammals, almost all DNA methylation is erased during embryo development and re-established in a sex- and cell-specific manner. This erasure and re-establishment is thought to primarily be a vertebrate-specific trait. Insects are particularly interesting in terms of development as many species often undergo remarkable morphological changes en route to maturity, that is, morphogenesis. However, little is known about the role of epigenetic mechanisms in this process across species. We have used whole-genome bisulfite sequencing to track genome-wide DNA methylation changes through the development of an economically and environmentally important pollinator species, the bumblebee Bombus terrestris (Hymenoptera:Apidae Linnaeus). We find overall levels of DNA methylation vary throughout development, and we find developmentally relevant differentially methylated genes throughout. Intriguingly, we have identified a depletion of DNA methylation in ovaries/eggs and an enrichment of highly methylated genes in sperm. We suggest this could represent a sex-specific DNA methylation erasure event. To our knowledge, this is the first suggestion of possible developmental DNA methylation erasure in an insect species. This study lays the required groundwork for functional experimental work to determine if there is a causal nature to the DNA methylation differences identified. Additionally, the application of single-cell methylation sequencing to this system will enable more accurate identification of if or when DNA methylation is erased during development.

Tracking Existing Factors Directly Affecting the Reproduction of Bumblebees: Current Knowledge

Bumblebees are primary social insects and a vital class of pollinating insects. Their distinctive reproductive mode is characterized by the independent initiation and construction of the nest by the queen and the subsequent production of sufficient workers, males, and gynes following colony development. After successful mating, the queen transitions to the first phase of its annual life cycle. The reproductive processes are directly influenced by environmental factors, including floral resources and pesticides. Moreover, the reproductive level is regulated by biological factors, particularly the role of workers, who participate in egg laying and pass on their genetic material to the next generation of queens. Successful reproduction can only be achieved by maintaining colony development under natural or artificial breeding conditions. Consequently, understanding the known factors that influence bumblebee reproduction is essential for developing conservation strategies for wild bumblebees and for successfully breeding diverse bumblebee species. Breeding various bumblebee species is crucial for in-depth research into known factors and for further exploration of other potential factors, which will also help to meet the demand for pollination in agricultural facilities globally.

DNA methylation is associated with codon degeneracy in a species of bumblebee

Social insects display extreme phenotypic differences between sexes and castes even though the underlying genome can be almost identical. Epigenetic processes have been proposed as a possible mechanism for mediating these phenotypic differences. Using whole genome bisulfite sequencing of queens, males, and reproductive female workers we have characterised the sex- and caste-specific methylome of the bumblebee Bombus terrestris. We have identified a potential role for DNA methylation in histone modification processes which may influence sex and caste phenotypic differences. We also find differentially methylated genes generally show low levels of DNA methylation which may suggest a separate function for lowly methylated genes in mediating transcriptional plasticity, unlike highly methylated genes which are usually involved in housekeeping functions. We also examined the relationship between the underlying genome and the methylome using whole genome re-sequencing of the same queens and males. We find DNA methylation is enriched at zero-fold degenerate sites. We suggest DNA methylation may be acting as a targeted mutagen at these sites, providing substrate for selection via non-synonymous changes in the underlying genome. However, we did not see any relationship between DNA methylation and rates of positive selection in our samples. In order to fully assess a possible role for DNA methylation in adaptive processes a specifically designed study using natural population data is needed.

Long-term sperm storage in eusocial Hymenoptera

In internally fertilizing species, sperm transfer is not always immediately followed by egg fertilization, and female sperm storage (FSS) may occur. FSS is a phenomenon in which females store sperm in a specialized organ for periods lasting from a few hours to several years, depending on the species. Eusocial hymenopterans (ants, social bees, and social wasps) hold the record for FSS duration. In these species, mating takes place during a single nuptial flight that occurs early in adult life for both sexes; they never mate again. Males die quickly after copulation but survive posthumously as sperm stored in their mates' spermathecae. Reproductive females, also known as queens, have a much longer life expectancy, up to 20 years in some species. Here, we review what is currently known about the molecular adaptations underlying the remarkable FSS capacities in eusocial hymenopterans. Because sperm quality is crucial to the reproductive success of both sexes, we also discuss the mechanisms involved in sperm storage and preservation in the male seminal vesicles prior to ejaculation. Finally, we propose future research directions that should broaden our understanding of this unique biological phenomenon.

DNA methylation and expression of the egfr gene are associated with worker size in monomorphic ants

The reproductive division of labour is a hallmark of eusocial Hymenoptera. Females are either reproductive queens or non-reproductive workers. In ants, workers often display further task specialisation that is associated with variation in size and/or morphology. Because female polyphenism is typically under environmental control, it is thought epigenetic mechanisms (such as DNA methylation) play a central role since they mediate gene-by-environment interactions. Methylation of the growth-promoting gene epidermal growth factor receptor (egfr) was indeed shown to control worker size in a highly polymorphic ant. However, it remains unknown if egfr methylation could also regulate worker size in monomorphic species. By combining experimental pharmacology and molecular biology, we show that worker size is associated with egfr methylation in two monomorphic ants. Furthermore, we functionally demonstrate that EGFR signalling affects worker size. These results indicate that worker size regulation by egfr methylation has been mechanistically conserved in ants but remains unexploited in monomorphic species.

Complex regulatory role of DNA methylation in caste- and age-specific expression of a termite

The reproductive castes of eusocial insects are often characterized by extreme lifespans and reproductive output, indicating an absence of the fecundity/longevity trade-off. The role of DNA methylation in the regulation of caste- and age-specific gene expression in eusocial insects is controversial. While some studies find a clear link to caste formation in honeybees and ants, others find no correlation when replication is increased across independent colonies. Although recent studies have identified transcription patterns involved in the maintenance of high reproduction throughout the long lives of queens, the role of DNA methylation in the regulation of these genes is unknown. We carried out a comparative analysis of DNA methylation in the regulation of caste-specific transcription and its importance for the regulation of fertility and longevity in queens of the higher termite Macrotermes natalensis. We found evidence for significant, well-regulated changes in DNA methylation in mature compared to young queens, especially in several genes related to ageing and fecundity in mature queens. We also found a strong link between methylation and caste-specific alternative splicing. This study reveals a complex regulatory role of fat body DNA methylation both in the division of labour in termites, and during the reproductive maturation of queens.

Molecular Identification of Two DNA Methyltransferase Genes and Their Functional Characterization in the Anti-Bacterial Immunity of Antheraea pernyi

Under different physiological conditions, such as microbial infection, epigenetic mechanisms regulate genes at the transcription level in living organisms. DNA methylation is a type of epigenetic mechanism in which DNA methyltransferases modify the expression of target genes. Here, we identified a full-length sequence of DNMT-1 and DNMT-2 from the Chinese oak silkworm, A. pernyi, which was highly similar to the homologous sequences of Bombyx mori. ApDNMT-1 and ApDNMT-2 have unique domain architectures of insect DNMTs, highlighting their conserved functions in A. pernyi. ApDNMT-1 and ApDNMT-2 were found to be widely expressed in various tissues, with the highest levels of expression in hemocytes, the ovary, testis, and fat bodies. To understand the biological role of these genes in microbial resistance, we challenged the fifth instar larvae of A. pernyi by administrating Gram-positive and Gram-negative bacteria and fungi. The results revealed that transcript levels of ApDNMT-1 and ApDNMT-2 were increased compared to the control group. The inhibition of these genes by a DNMTs inhibitor [5-azacytidine (5-AZA)] significantly reduced bacterial replication and larvae mortality. In addition, 5-AZA treatment modified the expression patterns of antimicrobial peptides (AMPs) in the A. pernyi larvae. Our results suggest that ApDNMT-1 and ApDNMT-2 seem to have a crucial role in innate immunity, mediating antimicrobial peptide responses against bacterial infection in A. pernyi.

DNA Methylation and Detoxification in the Earthworm Lumbricus terrestris Exposed to Cadmium and the DNA Demethylation Agent 5-aza-2'-deoxycytidine

Earthworms are well-established model organisms for testing the effects of heavy metal pollution. How DNA methylation affects cadmium (Cd) detoxification processes such as the expression of metallothionein 2 (MT2), however, is largely unknown. We therefore exposed Lumbricus terrestris to 200 mg concentrations of Cd and 5-aza-2'-deoxycytidine (Aza), a demethylating agent, and sampled tissue and coelomocytes, cells of the innate immune system, for 48 h. MT2 transcription significantly increased in the Cd- and Cd-Aza-treated groups. In tissue samples, a significant decrease in MT2 in the Aza-treated group was detected, showing that Aza treatment inhibits basal MT2 gene activity but has no effect on Cd-induced MT2 levels. Although Cd repressed the gene expression of DNA-(cytosine-5)-methyltransferase-1 (DNMT1), which is responsible for maintaining DNA methylation, DNMT activity was unchanged, meaning that methylation maintenance was not affected in coelomocytes. The treatment did not influence DNMT3, which mediates de novo methylation, TET gene expression, which orchestrates demethylation, and global levels of hydroxymethylcytosine (5hmC), a product of the demethylation process. Taken together, this study indicates that Aza inhibits basal gene activity, in contrast to Cd-induced MT2 gene expression, but does not affect global DNA methylation. We therefore conclude that Cd detoxification based on the induction of MT2 does not relate to DNA methylation changes.

Phenotypic Plasticity: What Has DNA Methylation Got to Do with It?

How does one genome give rise to multiple, often markedly different, phenotypes in response to an environmental cue? This phenomenon, known as phenotypic plasticity, is common amongst plants and animals, but arguably the most striking examples are seen in insects. Well-known insect examples include seasonal morphs of butterfly wing patterns, sexual and asexual reproduction in aphids, and queen and worker castes of eusocial insects. Ultimately, we need to understand how phenotypic plasticity works at a mechanistic level; how do environmental signals alter gene expression, and how are changes in gene expression translated into novel morphology, physiology and behaviour? Understanding how plasticity works is of major interest in evolutionary-developmental biology and may have implications for understanding how insects respond to global change. It has been proposed that epigenetic mechanisms, specifically DNA methylation, are the key link between environmental cues and changes in gene expression. Here, we review the available evidence on the function of DNA methylation of insects, the possible role(s) for DNA methylation in phenotypic plasticity and also highlight key outstanding questions in this field as well as new experimental approaches to address these questions.

Reproductive Potential Impacts Body Maintenance Parameters and Global DNA Methylation in Honeybee Workers (Apis mellifera L.)

Simple Summary

The queens and sterile workers arise from genetically identical eggs but as imagoes, they differ in their life span, DNA methylation, and their functions. In the absence of the queen, the larvae develop into rebels, i.e., workers with increased reproductive potential. We assumed that since rebels are similar to the queen in many anatomical and behavioral features, they live longer and have lower levels of global DNA methylation, even when infected, e.g., by Nosema spp. Rebels always lived longer in comparison in normal workers and unexpectedly extended longevity of normal workers when they were together, similarly as the presence of a queen did. Rebels became infected more easily but tolerated the infection better. They also had lower level of global DNA methylation than normal workers. These features expand possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.

Abstract

The widely accepted hypothesis in life history evolution about the trade-off between fecundity and longevity is not confirmed by long-living and highly fecund queens in eusocial insects. The fact that the queens and facultatively sterile workers usually arise from genetically identical eggs but differ in DNA methylation makes them a good model for studies on senescence, eusocial evolution, and epigenetics. Therefore, honeybees seem to be especially useful here because of long living rebel-workers (RW) with high reproductive potential recently described. Longevity, ovariole number, nosema tolerance, and global DNA methylation have been assayed in normal workers (NW) versus RW in hives and cages. RW always lived longer than NW and unexpectedly extended longevity of NW when they were together, similarly as the presence of a queen did. RW lived longer despite the fact that they had higher Nosema spore load; surprisingly they became infected more easily but tolerated the infection better. Global DNA methylation increased with age, being lower in RW than in NW. Therefore, RW are queen-like considering global DNA methylation and the link between fecundity, longevity, and body maintenance. Presented features of RW expands possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.

Worker-Born Males Are Smaller but Have Similar Reproduction Ability to Queen-Born Males in Bumblebees

Queen-worker conflict over the reproduction of males exists in the majority of haplodiplioidy hymenpteran species such as bees, wasps, and ants, whose workers lose mating ability but can produce haploid males in colony. Bumblebee is one of the representatives of primitively eusocial insects with plastic division labor and belongs to monandrous and facultative low polyandry species that have reproductive totipotent workers, which are capable of competing with mother queen to produce haploid males in the queenright colony compared to higher eusocial species, e.g., honeybees. So, bumblebees should be a better material to study worker reproduction, but the reproductive characteristics of worker-born males (WMs) remain unclear. Here, we choose the best-studied bumblebee Bombus terrestris to evaluate the morphological characteristics and reproductive ability of WMs from the queenless micro-colonies. The sexually matured WMs showed smaller in forewing length and weight, relatively less sperm counts but equally high sperm viability in comparison with the queen-born males (QMs) of the queenright colony. Despite with smaller size, the WMs are able to successfully mate with the virgin queens in competition with the QMs under laboratory conditions, which is quite different from the honeybees reported. In addition, there was no difference in the colony development, including the traits such as egg-laying rate, colony establishment rate, and populations of offspring, between the WM- and the QM-mated queens. Our study highlights the equivalent reproductive ability of worker-born males compared to that of queens, which might exhibit a positive application or special use of bumblebee rearing, especially for species whose males are not enough for copulation. Further, our finding contributes new evidence to the kin selection theory and suggests worker reproduction might relate to the evolution of sociality in bees.

Independent variations in genome-wide expression, alternative splicing, and DNA methylation in brain tissues among castes of the buff-tailed bumblebee, Bombus terrestris

Caste differentiation in social hymenopterans is an intriguing example of phenotypic plasticity. However, the co-ordination among gene regulatory factors to mediate caste differentiation remains inconclusive. In this study, we determined the role of gene regulation and related epigenetic processes in pre-imaginal caste differentiation in the primitively eusocial bumblebee Bombus terrestris. By combining RNA-Seq data from Illumina and PacBio and accurately quantifying methylation at whole-genomic base pair resolution, we found that queens, workers, and drones mainly differentiate in gene expression but not in alternative splicing and DNA methylation. Gynes are the most distinct with the lowest global level of whole-genomic methylation and with the largest number of caste-specific transcripts and alternative splicing events. By contrast, workers exhibit few uniquely expressed or alternatively spliced genes. Moreover, several genes involved in hormone and neurotransmitter metabolism are related to caste differentiation, whereas several neuropeptides are linked with sex differentiation. Despite little genome-wide association among differential gene expression, splicing, and differential DNA methylation, the overlapped gene ontology (GO) terms point to nutrition-related activity. Therefore, variations in gene regulation correlate with the behavioral differences among castes and highlight the specialization of toolkit genes in bumblebee gynes at the beginning of the adult stage.

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

10-hydroxy-2E-decenoic acid (10HDA) does not promote caste differentiation in Melipona scutellaris stingless bees

In bees from genus Melipona, differential feeding is not enough to fully explain female polyphenism. In these bees, there is a hypothesis that in addition to the environmental component (food), a genetic component is also involved in caste differentiation. This mechanism has not yet been fully elucidated and may involve epigenetic and metabolic regulation. Here, we verified that the genes encoding histone deacetylases HDAC1 and HDAC4 and histone acetyltransferase KAT2A were expressed at all stages of Melipona scutellaris, with fluctuations between developmental stages and castes. In larvae, the HDAC genes showed the same profile of Juvenile Hormone titers-previous reported-whereas the HAT gene exhibited the opposite profile. We also investigated the larvae and larval food metabolomes, but we did not identify the putative queen-fate inducing compounds, geraniol and 10-hydroxy-2E-decenoic acid (10HDA). Finally, we demonstrated that the histone deacetylase inhibitor 10HDA-the major lipid component of royal jelly and hence a putative regulator of honeybee caste differentiation-was unable to promote differentiation in queens in Melipona scutellaris. Our results suggest that epigenetic and hormonal regulations may act synergistically to drive caste differentiation in Melipona and that 10HDA is not a caste-differentiation factor in Melipona scutellaris.

The effect of DNA methylation on bumblebee colony development

BACKGROUND

Although around 1% of cytosines in bees' genomes are known to be methylated, less is known about methylation's effect on bee behavior and fitness. Chemically altered DNA methylation levels have shown clear changes in the dominance and reproductive behavior of workers in queen-less colonies, but the global effect of DNA methylation on caste determination and colony development remains unclear, mainly because of difficulties in controlling for genetic differences among experimental subjects in the parental line. Here, we investigated the effect of the methylation altering agent decitabine on the developmental rate of full bumblebee colonies. Whole genome bisulfite sequencing was used to assess differences in methylation status.

RESULTS

Our results showed fewer methylated loci in the control group. A total of 22 CpG loci were identified as significantly differentially methylated between treated and control workers with a change in methylation levels of 10% or more. Loci that were methylated differentially between groups participated in pathways including neuron function, oocyte regulation and metabolic processes. Treated colonies tended to develop faster, and therefore more workers were found at a given developmental stage. However, male production followed the opposite trend and it tended to be higher in control colonies.

CONCLUSION

Overall, our results indicate that altered methylation patterns resulted in an improved cooperation between workers, while there were no signs of abnormal worker dominance or caste determination.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

Parent of origin gene expression in the bumblebee, Bombus terrestris, supports Haig's kinship theory for the evolution of genomic imprinting

Genomic imprinting is the differential expression alleles in diploid individuals, with the expression being dependent on the sex of the parent from which it was inherited. Haig's kinship theory hypothesizes that genomic imprinting is due to an evolutionary conflict of interest between alleles from the mother and father. In social insects, it has been suggested that genomic imprinting should be widespread. One recent study identified parent-of-origin expression in honey bees and found evidence supporting the kinship theory. However, little is known about genomic imprinting in insects and multiple theoretical predictions must be tested to avoid single-study confirmation bias. We, therefore, tested for parent-of-origin expression in a primitively eusocial bee. We found equal numbers of maternally and paternally biased expressed genes. The most highly biased genes were maternally expressed, offering support for the kinship theory. We also found low conservation of potentially imprinted genes with the honey bee, suggesting rapid evolution of genomic imprinting in Hymenoptera.

Multi-omics analysis reveals that natural hibernation is crucial for oocyte maturation in the female Chinese alligator

BACKGROUND

Hibernation in an appropriate environment not only is important for the survival of hibernators in winter, but also is crucial for breeding in the following season for many hibernating species. However, the genetic and epigenetic mechanism underlying this process remain unclear. In the current study, we performed an integrative multi-omics analysis of gonads collected from Chinese alligators that overwintered in wild cave and artificial warmroom to explore transcriptomic and epigenomic alternations in these organs.

RESULTS

The data revealed that in the breeding season, female alligators were more strongly affected in terms of gene expression than males by non-hibernation because of overwintering in a warm room, especially for genes related to oocyte maturation, and this effect commenced in winter with the downregulation of STAR, which is the rate limiting factor of steroid biosynthesis. Further, miRNAs were found to play essential roles in this negative effect of overwintering in the warm room on hibernation. The upregulated miRNAs likely were responsible for the suppression of oocyte maturation in the breeding season. Finally, DNA methylome changes, especially hypomethylation, were found to play an important role in the alterations in ovarian function-related gene expression induced by non-hibernation.

CONCLUSIONS

Our study revealed the crucial role of hibernation quality for oocyte maturation in the Chinese alligator and the underlying genetic and epigenetic mechanisms, and highlights the importance of habitat, and especially, the overwintering site, in the conservation of not only the Chinese alligator, but also other endangered hibernators.

A novel method for using RNA-seq data to identify imprinted genes in social Hymenoptera with multiply mated queens

Genomic imprinting results in parent-of-origin-dependent gene expression biased towards either the maternally or paternally derived allele at the imprinted locus. The kinship theory of genomic imprinting argues that this unusual expression pattern can be a manifestation of intra-genomic conflict between the maternally and paternally derived halves of the genome that arises because they are not equally related to the genomes of social partners. The theory thus predicts that imprinting may evolve wherever there are close interactions among asymmetrically related kin. The social Hymenoptera with permanent caste differentiation are suitable candidates for testing the kinship theory because haplodiploid sex determination creates strong relatedness asymmetries and nursing workers interact closely with kin. However, progress in the search for imprinted genes in the social Hymenoptera has been slow, in part because tests for imprinting rely on reciprocal crosses that are impossible in most species. Here, we develop a method to systematically search for imprinting in haplodiploid social insects without crosses, using instead samples of pooled individuals collected from natural colonies. We tested this protocol using data available for the leaf-cutting ant Acromyrmex echinatior, providing the first genome-wide search for imprinting in any ant. Although we identified several genes as potentially imprinted, none of the four genes tested could be verified as imprinted using digital droplet PCR, highlighting the need for higher quality genomic assemblies that accurately map duplicated genes.

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.

DNA methylation patterns respond to thermal stress in the viviparous cockroach Diploptera punctata

It is increasingly recognized that epigenetic mechanisms play a key role in acclimatization and adaptation to thermal stress in invertebrates. DNA methylation and its response to temperature variation has been poorly studied in insects. Here, we investigated DNA methylation and hydroxymethylation patterns in the viviparous cockroach Diploptera punctata at a global and gene specific level in response to variation in temperature. We specifically studied methylation percentage in the heat shock protein 70 (Hsp70), whose function is linked to thermal plasticity and resistance. We found high levels of DNA methylation in several tissues but only low levels of DNA hydroxymethylation in the brain. Hsp70 methylation patterns showed significant differences in response to temperature. We further found that global DNA methylation variation was considerably lower at 28°C compared to higher or lower temperatures, which may be indicative of the optimal temperature for this species. Our results demonstrate that DNA methylation could provide a mechanism for insects to dynamically respond to changing temperature conditions in their environment.

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.

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.

Induction of DNA methyltransferase genes in Helicoverpa armigera following injection of pathogenic bacteria modulates expression of antimicrobial peptides and affects bacterial proliferation

Following pathogen attack in a host, widespread changes are induced in the host's gene expression, in particular those involved in the immune system, growth and survival. Epigenetic mechanisms have been suggested to be involved in the regulation of these changes through a number of mechanisms. DNA methylation is one of the important epigenetic processes that is carried out by DNA (cytosine-5) methyltransferase (DNMT) and alters expression of target genes. Here, we identified two putative sequences of DNMT (i.e. DNMT1 and DNMT2) from the transcriptome dataset of Helicoverpa armigera that showed high similarity to the homologous sequences in Bombyx mori. Domain architectures of DNMT1 and DNMT2 exhibit the unique pattern of DNMTs that highlights conserved function of these genes in different insects. To see if these genes play any role in bacterial infection, we challenged the fifth instar larvae of H. armigera by injecting Bacillus thuringiensis and Serratia marcescens cells into the hemolymph. Transcript levels of the DNMTs were analyzed by RT-qPCR. The results showed that the expression levels of DNMT1 and DNMT2 increased in the bacteria-injected larvae. Injection of the heat-killed bacteria also induced the expression of the DNMTs, but lower than that of the live bacteria. To determine whether these genes function during bacterial infection, we injected the inhibitor of DNMTs, 5-azacytidine (5-AZA), into the larvae and 24 h later, the bacterial cells were also injected into the larvae. Bacterial replication and larval mortality were analyzed in the treated and control insects. We found that 5-AZA reduced bacterial replication and also mortality of the bacterial-injected larvae regardless of the pathogenic bacterial species. Interestingly, the expression levels of antimicrobial peptides (AMPs) were also modulated following 5-AZA treatment. In conclusion, we showed that upregulation of the DNMTs in H. armigera following bacterial infections modulates AMPs and thereby affects the insect-bacteria interactions.

Validating the Demethylating Effects of 5-aza-2'-deoxycytidine in Insects Requires a Whole-Genome Approach

We previously demonstrated that treatment with the demethylating agent 5-aza-2'-deoxycytidine (5-aza-dC) alters the offspring sex ratios produced by females of the parasitoid wasp Nasonia vitripennis. Females allocate offspring sex ratio in line with local mate competition theory, producing more or less female-biased sex ratios as the number of other females laying eggs on a patch varies, thereby reducing competition among their sons for mates. Interestingly, treatment with 5-aza-dC did not ablate the facultative sex allocation response. Instead, sex ratios became less female biased, a shift in the direction of the optimum sex ratio for paternally inherited alleles according to genomic conflict theory. This was the first (albeit indirect) experimental evidence for genomic conflict over sex allocation. In their comment, Ellers and colleagues assayed the effects of 5-aza-dC on DNA methylation in 10 Nasonia genes, finding no evidence of demethylation in these 10 genes, from which they conclude that 5-aza-dC has no demethylating capability in N. vitripennis. Quantifying the efficacy of 5-aza-dC in terms of demethylation is indeed crucial to in-depth interpretation of studies using 5-aza-dC to link phenotypes to epigenetic regulation. Here we outline the mode of action of 5-aza-dC and demonstrate that determining the efficacy of 5-aza-dC in insect systems requires a whole-genome approach.

The Importance of Validating the Demethylating Effect of 5-aza-2'-deoxycytidine in Model Species (A Comment on Cook et al., "DNA Methylation and Sex Allocation in the Parasitoid Wasp Nasonia vitripennis")

The use of DNA demethylating agents has been popular in epigenetic studies. Recently, Cook and colleagues, in a 2015 American Naturalist article, claimed an effect of 5-aza-2'-deoxycytidine (5-aza-dC) on the sex ratio of a parasitoid wasp without verifying its effect on DNA methylation. We repeated the 5-aza-dC feeding treatment to test its effectiveness. We used bisulfite amplicon sequencing of 10 genes that either were heavily methylated, previously showed a response to 5-aza-dC, or were suggested to regulate fatty acid synthesis epigenetically, and we demonstrate that wasps fed 5-aza-dC did not show reduced DNA methylation at these loci. Therefore, the conclusion that demethylation shifts sex ratios upward needs reconsideration.

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.

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.

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.

Queen Longevity and Fecundity Affect Conflict with Workers over Resource Inheritance in a Social Insect

Resource inheritance is a major source of conflict in animal societies. However, the assumptions and predictions of models of conflict over resource inheritance have not been systematically tested within a single system. We developed an inclusive fitness model for annual eusocial Hymenoptera that predicts a zone of conflict in which future reproductive workers are selected to enforce nest inheritance before the queen is selected to cede the nest. We experimentally tested key elements of this model in the bumblebee Bombus terrestris. In colonies from which queens were sequentially removed, queen tenure was significantly negatively associated with worker male production, confirming that workers gain direct fitness by usurping the queen. In unmanipulated colonies, queen fecundity decreased significantly over the latter part of the colony cycle, confirming that workers' indirect fitness from maintaining queens declines over time. Finally, in an experiment simulating loss of queen fecundity by removal of queens' eggs, worker-to-queen aggression increased significantly and aggressive workers were significantly more likely to become egg layers, consistent with workers monitoring queen fecundity to assess the net benefit of future reproduction. Overall, by upholding key assumptions and predictions of the model, our results provide novel empirical support for kin-selected conflict over resource inheritance.

Divergent methylation pattern in adult stage between two forms of Tetranychus urticae (Acari: Tetranychidae)

The two-spotted spider mite, Tetranychus urticae Koch has two forms: green form and red form. Understanding the molecular basis of how these two forms established without divergent genetic background is an intriguing area. As a well-known epigenetic process, DNA methylation has particularly important roles in gene regulation and developmental variation across diverse organisms that do not alter genetic background. Here, to investigate whether DNA methylation could be associated with different phenotypic consequences in the two forms of T. urticae, we surveyed the genome-wide cytosine methylation status and expression level of DNA methyltransferase 3 (Tudnmt3) throughout their entire life cycle. Methylation-sensitive amplification polymorphism (MSAP) analyses of 585 loci revealed variable methylation patterns in the different developmental stages. In particular, principal coordinates analysis (PCoA) indicates a significant epigenetic differentiation between female adults of the two forms. The gene expression of Tudnmt3 was detected in all examined developmental stages, which was significantly different in the adult stage of the two forms. Together, our results reveal the epigenetic distance between the two forms of T. urticae, suggesting that DNA methylation might be implicated in different developmental demands, and contribute to different phenotypes in the adult stage of these two forms.

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.

Identification and caste-dependent expression patterns of DNA methylation associated genes in Bombus terrestris

DNA methylation has been proposed to play critical roles in caste fate and behavioral plasticity in bumblebees, however, there is little information on its regulatory mechanisms. Here, we identified six important genes mediating the modification of DNA methylation and determined their expression patterns in the bumblebee Bombus terrestris. There is a complete functional DNA methylation system, including four DNA methyltransferases (DNMT1a, DNMT1b, DNMT2, and DNMT3), a DNA demethylase (Ten-eleven translocation), and a methyl-CpG-binding domain protein in B. terrestris. Most of these genes were highly expressed in fat bodies and gonads but lowly expressed in antennae and brains of bumblebee adults. Besides, these genes exhibited caste-specific expression patterns in bumblebees, with higher transcription levels in queens than workers and drones. Whereas their expression levels showed no remarkable difference in queenright and queenless workers. These results suggested potential roles of DNA methylation-related genes in caste differentiation in bumblebees.

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

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

Allele specific expression and methylation in the bumblebee, Bombus terrestris

The social hymenoptera are emerging as models for epigenetics. DNA methylation, the addition of a methyl group, is a common epigenetic marker. In mammals and flowering plants methylation affects allele specific expression. There is contradictory evidence for the role of methylation on allele specific expression in social insects. The aim of this paper is to investigate allele specific expression and monoallelic methylation in the bumblebee, Bombus terrestris. We found nineteen genes that were both monoallelically methylated and monoallelically expressed in a single bee. Fourteen of these genes express the hypermethylated allele, while the other five express the hypomethylated allele. We also searched for allele specific expression in twenty-nine published RNA-seq libraries. We found 555 loci with allele-specific expression. We discuss our results with reference to the functional role of methylation in gene expression in insects and in the as yet unquantified role of genetic cis effects in insect allele specific methylation and expression.

Insights into the biochemical defence and methylation of the solitary bee Osmia rufa L: A foundation for examining eusociality development

We examined age-related biochemical and histological changes in the fat bodies and hemolymph of Osmia rufa males and females. We analysed solitary bees during diapause, in October and in April; as well as the flying insects following diapause, in May and June. The trophocyte sizes, as well as the numbers of lipid droplets were the greatest at the beginning of diapause. Subsequently, they decreased along with age. Triglyceride and glucose concentrations systematically decreased in fat body cells but increased in the hemolymph from October to June. Concentrations/activities of (enzymatic and non-enzymatic) antioxidant and proteolytic systems, as well as phenoloxidase, aspartate aminotransferase, alanine aminotransferase and alkaline phosphatase levels were constant during the diapause, usually lower in the males than the females. Prior to the diapause/overwintering, the concentrations/activities of all the compounds were higher in the fat bodies than in the hemolymph. Later in the spring and in the summer, they increased in the hemolymph and on the body surfaces, while decreasing in the fat bodies. The global DNA methylation levels increased with age. Higher levels were always observed in the males than in the females. The study will promote better understanding of bee evolution and will be useful for the protection and management of solitary bees, with benefits to the environment and agriculture.

Population genomics reveals a candidate gene involved in bumble bee pigmentation

Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.

Sexual selection modulates genetic conflicts and patterns of genomic imprinting

Recent years have seen a surge of interest in linking the theories of kin selection and sexual selection. In particular, there is a growing appreciation that kin selection, arising through demographic factors such as sex-biased dispersal, may modulate sexual conflicts, including in the context of male-female arms races characterized by coevolutionary cycles. However, evolutionary conflicts of interest need not only occur between individuals, but may also occur within individuals, and sex-specific demography is known to foment such intragenomic conflict in relation to social behavior. Whether and how this logic holds in the context of sexual conflict-and, in particular, in relation to coevolutionary cycles-remains obscure. We develop a kin-selection model to investigate the interests of different genes involved in sexual and intragenomic conflict, and we show that consideration of these conflicting interests yields novel predictions concerning parent-of-origin specific patterns of gene expression and the detrimental effects of different classes of mutation and epimutation at loci underpinning sexually selected phenotypes.

Queen pheromones modulate DNA methyltransferase activity in bee and ant workers

DNA methylation is emerging as an important regulator of polyphenism in the social insects. Research has concentrated on differences in methylation between queens and workers, though we hypothesized that methylation is involved in mediating other flexible phenotypes, including pheromone-dependent changes in worker behaviour and physiology. Here, we find that exposure to queen pheromone affects the expression of two DNA methyltransferase genes in Apis mellifera honeybees and in two species of Lasius ants, but not in Bombus terrestris bumblebees. These results suggest that queen pheromones influence the worker methylome, pointing to a novel proximate mechanism for these key social signals.

Genome methylation patterns across castes and generations in a parasitoid wasp

Environmental influences shape phenotypes within and across generations, often through DNA methylations that modify gene expression. Methylations were proposed to mediate caste and task allocation in some eusocial insects, but how an insect's environment affects DNA methylation in its offspring is yet unknown. We characterized parental effects on methylation profiles in the polyembryonic parasitoid wasp Copidosoma koehleri, as well as methylation patterns associated with its simple caste system. We used methylation‐sensitive amplified fragment length polymorphism (MS‐AFLP) to compare methylation patterns, among (1) reproductive and soldier larvae; and (2) offspring (larvae, pupae, and adults) of wasps that were reared at either high or low larval density and mated in the four possible combinations. Methylation frequencies were similar across castes, but the profiles of methylated fragments differed significantly. Parental rearing density did not affect methylation frequencies in the offspring at any developmental stage. Principal coordinate analysis indicated no significant differences in methylation profiles among the four crossbreeding groups and the three developmental stages. Nevertheless, a clustering analysis, performed on a subset of the fragments, revealed similar methylation patterns in larvae, pupae, and adults in two of the four parental crosses. Nine fragments were methylated at two cytosine sites in all larvae, and five others were methylated at two sites in all adults. Thus, DNA methylations correlate with within‐generation phenotypic plasticity due to caste. However, their association with developmental stage and with transgenerational epigenetic effects is not clearly supported.

Chapter 6. Comparative genomic approaches to investigate molecular traits specific to social insects

Ageing is a feature of nearly all known organisms and, by its connection to survival, appears to trade off with fecundity. However, in some organisms such as in queens of social insects, this negative relation appears reversed and individuals live long and reproduce much. Since new experimental techniques, transcriptomes and genomes of many social insects have recently become available, a comparison of these data in a phylogenetic framework becomes feasible. This allows the study of general trends, species specific oddities and evolutionary dynamics of the molecular properties and changes which underlie ageing, fecundity and the reversal of this negative association. In the framework of social insect evolution, we review the most important recent insights, computational methods, their applications and data resources which are available.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

Exploring complex pheromone biosynthetic processes in the bumblebee male labial gland by RNA sequencing

Male marking pheromones (MPs) are used by the majority of bumblebee species (Hymenoptera: Apidae), including a commercially important greenhouse pollinator, the buff-tailed bumblebee (Bombus terrestris), to attract conspecific females. MP biosynthetic processes in the cephalic part of the bumblebee male labial gland (LG) are of extraordinary complexity, involving enzymes of fatty acid and isoprenoid biosynthesis, which jointly produce more than 50 compounds. We employed a differential transcriptomic approach to identify candidate genes involved in MP biosynthesis by sequencing Bombus terrestris LG and fat body (FB) transcriptomes. We identified 12 454 abundantly expressed gene products (reads per kilobase of exon model per million mapped reads value > 1) that had significant hits in the GenBank nonredundant database. Of these, 876 were upregulated in the LG (> 4-fold difference). We identified more than 140 candidate genes potentially involved in MP biosynthesis, including esterases, fatty acid reductases, lipases, enzymes involved in limited fatty acid chain shortening, neuropeptide receptors and enzymes involved in biosynthesis of triacylglycerols, isoprenoids and fatty acids. For selected candidates, we confirmed their abundant expression in LG using quantitative real-time reverse transcription-PCR (qRT-PCR). Our study shows that the Bombus terrestris LG transcriptome reflects both fatty acid and isoprenoid MP biosynthetic processes and identifies rational gene targets for future studies to disentangle the molecular basis of MP biosynthesis. Additionally, LG and FB transcriptomes enrich the available transcriptomic resources for Bombus terrestris.

DNA Methylation and Potential for Epigenetic Regulation in Pygospio elegans

Transitions in developmental mode are common evolutionarily, but how and why they occur is not understood. Developmental mode describes larval phenotypes, including morphology, ecology and behavior of larvae, which typically are generalized across different species. The polychaete worm Pygospio elegans is one of few species polymorphic in developmental mode, with multiple larval phenotypes, providing a possibility to examine the potential mechanisms allowing transitions in developmental mode. We investigated the presence of DNA methylation in P. elegans, and, since maternal provisioning is a key factor determining eventual larval phenotype, we compared patterns of DNA methylation in females during oogenesis in this species. We demonstrate that intragenic CpG site DNA methylation and many relevant genes necessary for DNA methylation occur in P. elegans. Methylation-sensitive AFLP analysis showed that gravid females with offspring differing in larval developmental mode have significantly different methylation profiles and that the females with benthic larvae and non-reproductive females from the same location also differ in their epigenetic profiles. Analysis of CpG sites in transcriptome data supported our findings of DNA methylation in this species and showed that CpG observed/expected ratios differ among females gravid with embryos destined to different developmental modes. The differences in CpG site DNA methylation patterns seen among the samples suggest a potential for epigenetic regulation of gene expression (through DNA methylation) in this species.

Intragenomic Conflict over Soldier Allocation in Polyembryonic Parasitoid Wasps

Understanding the selection pressures that have driven the evolution of sterile insect castes has been the focus of decades of intense scientific debate. An amenable empirical test bed for theory on this topic is provided by the sterile-soldier caste of polyembryonic parasitoid wasps. The function of these soldiers has been a source of controversy, with two basic hypotheses emerging: the "brood-benefit" hypothesis that they provide an overall benefit for their siblings and the "sex-ratio-conflict" hypothesis that the soldiers mediate a conflict between brothers and sisters by killing their opposite-sex siblings. Here, we investigate the divergent sex-ratio optima of a female embryo's maternal-origin and paternal-origin genes, to determine the potential for, and direction of, intragenomic conflict over soldiering. We then derive contrasting empirically testable predictions concerning the patterns of genomic imprinting that are expected to arise out of this intragenomic conflict, for the brood-benefit versus the sex-ratio-conflict hypothesis of soldier function.

Gene expression differences in relation to age and social environment in queen and worker bumble bees

Eusocial insects provide special insights into the genetic pathways influencing aging because of their long-lived queens and flexible aging schedules. Using qRT-PCR in the primitively eusocial bumble bee Bombus terrestris (Linnaeus), we investigated expression levels of four candidate genes associated with taxonomically widespread age-related pathways (coenzyme Q biosynthesis protein 7, COQ7; DNA methyltransferase 3, Dnmt3; foraging, for; and vitellogenin, vg). In Experiment 1, we tested how expression changes with queen relative age and productivity. We found a significant age-related increase in COQ7 expression in queen ovary. In brain, all four genes showed higher expression with increasing female (queen plus worker) production, with this relationship strengthening as queen age increased, suggesting a link with the positive association of fecundity and longevity found in eusocial insect queens. In Experiment 2, we tested effects of relative age and social environment (worker removal) in foundress queens and effects of age and reproductive status in workers. In this experiment, workerless queens showed significantly higher for expression in brain, as predicted if downregulation of for is associated with the cessation of foraging by foundress queens following worker emergence. Workers showed a significant age-related increase in Dnmt3 expression in fat body, suggesting a novel association between aging and methylation in B. terrestris. Ovary activation was associated with significantly higher vg expression in fat body and, in younger workers, in brain, consistent with vitellogenin's ancestral role in regulating egg production. Overall, our findings reveal a mixture of novel and conserved features in age-related genetic pathways under primitive eusociality.

Genomic sources of phenotypic novelty in the evolution of eusociality in insects

Genomic resources are now available for closely related species that vary in social behavior, providing insight on the genomics of social evolution. Changes in the architecture of gene regulatory networks likely influence the evolutionary trajectory of social traits. Evolutionarily novel genes are likely important in the evolution of social diversity among insects, but it is unclear whether new genes played a driving role in the advent or elaboration of eusociality or if they were instead a result of other genomic features of eusociality. The worker phenotype appears to be the center of genetic novelty, but the mechanisms for this remain unresolved. Future studies are needed to understand how genetic novelty arises, becomes incorporated into existing gene regulatory networks, and the effects this has on the evolution of social traits in closely related social and solitary species.