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.

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.

Reproductive workers show queenlike gene expression in an intermediately eusocial insect, the buff-tailed bumble bee Bombus terrestris.

Bumble bees represent a taxon with an intermediate level of eusociality within Hymenoptera. The clear division of reproduction between a single founding queen and the largely sterile workers is characteristic for highly eusocial species, whereas the morphological similarity between the bumble bee queen and the workers is typical for more primitively eusocial hymenopterans. Also, unlike other highly eusocial hymenopterans, division of labour among worker subcastes is plastic and not predetermined by morphology or age. We conducted a differential expression analysis based on RNA-seq data from 11 combinations of developmental stage and caste to investigate how a single genome can produce the distinct castes of queens, workers and males in the buff-tailed bumble bee Bombus terrestris. Based on expression patterns, we found males to be the most distinct of all adult castes (2411 transcripts differentially expressed compared to nonreproductive workers). However, only relatively few transcripts were differentially expressed between males and workers during development (larvae: 71 and pupae: 162). This indicates the need for more distinct expression patterns to control behaviour and physiology in adults compared to those required to create different morphologies. Among female castes, reproductive workers and their nonreproductive sisters displayed differential expression in over ten times more transcripts compared to the differential expression found between reproductive workers and their mother queen. This suggests a strong shift towards a more queenlike behaviour and physiology when a worker becomes fertile. This contrasts with eusocial species where reproductive workers are more similar to nonreproductive workers than the queen.

Differential gene expression and alternative splicing in insect immune specificity.

BACKGROUND: Ecological studies routinely show genotype-genotype interactions between insects and their parasites. The mechanisms behind these interactions are not clearly understood. Using the bumblebee Bombus terrestris/trypanosome Crithidia bombi model system (two bumblebee colonies by two Crithidia strains), we have carried out a transcriptome-wide analysis of gene expression and alternative splicing in bees during C. bombi infection. We have performed four analyses, 1) comparing gene expression in infected and non-infected bees 24 hours after infection by Crithidia bombi, 2) comparing expression at 24 and 48 hours after C. bombi infection, 3) determining the differential gene expression associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours after infection and 4) determining the alternative splicing associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours post infection. RESULTS: We found a large number of genes differentially regulated related to numerous canonical immune pathways. These genes include receptors, signaling pathways and effectors. We discovered a possible interaction between the peritrophic membrane and the insect immune system in defense against Crithidia. Most interestingly, we found differential expression and alternative splicing of immunoglobulin related genes (Dscam and Twitchin) are associated with the genotype-genotype interactions of the given bumblebee colony and Crithidia strain. CONCLUSIONS: In this paper we have shown that the expression and alternative splicing of immune genes is associated with specific interactions between different host and parasite genotypes in this bumblebee/trypanosome model.

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

Insects are at the dawn of an epigenetics era. Numerous social insect species have been found to possess a functioning methylation system, previously not thought to exist in insects. Methylation, an epigenetic tag, may be vital for the sociality and division of labour for which social insects are renowned. In the bumble-bee Bombus terrestris, we found methylation differences between the genomes of queenless reproductive workers and queenless non-reproductive workers. In a follow up experiment, queenless workers whose genomes had experimentally altered methylation were more aggressive and more likely to develop ovaries compared with control queenless workers. This shows methylation is important in this highly plastic reproductive division of labour. Methylation is an epigenetic tag for genomic imprinting (GI). It is intriguing that the main theory to explain the evolution of GI predicts that GI should be important in this worker reproduction behaviour.

Lack of paternal silencing and ecotype-specific expression in head and body lice hybrids.

Paternal genome elimination (PGE) is a non-Mendelian inheritance system, described in numerous arthropod species, in which males develop from fertilized eggs, but their paternally inherited chromosomes are eliminated before or during spermatogenesis. Therefore, PGE males only transmit their maternally inherited set of chromosomes to their offspring. In addition to the elimination of paternal chromosomes, diverse PGE species have also repeatedly evolved the transcriptional silencing of the paternal genome, making males effectively haploid. However, it is unclear if this paternal chromosome silencing is mechanistically linked to the chromosome elimination or has evolved at a later stage, and if so, what drives the haploidization of males under PGE. In order to understand these questions, here we study the human louse, Pediculus humanus, which represents an ideal model system, as it appears to be the only instance of PGE where males eliminate, but not silence their paternal chromosomes, although the latter remains to be shown conclusively. In this study, we analyzed parent-of-origin allele-specific expression patterns in male offspring of crosses between head and body lice ecotypes. We show that hybrid adult males of P. humanus display biparental gene expression, which constitutes the first case of a species with PGE in which genetic activity of paternal chromosomes in the soma is not affected by embryonic silencing or (partial or complete) elimination. We did however also identify a small number of maternally biased genes (potentially imprinted genes), which may be involved in the elimination of paternal chromosomes during spermatogenesis. Finally, we have identified genes that show ecotype-specific expression bias. Given the low genetic diversity between ecotypes, this is suggestive for a role of epigenetic processes in ecotype differences.

Sexual morph specialisation in a trioecious nematode balances opposing selective forces.

The coexistence of different mating strategies, whereby a species can reproduce both by selfing and outcrossing, is an evolutionary enigma. Theory predicts two predominant stable mating states: outcrossing with strong inbreeding depression or selfing with weak inbreeding depression. As these two mating strategies are subject to opposing selective forces, mixed breeding systems are thought to be a rare transitory state yet can persist even after multiple speciation events. We hypothesise that if each mating strategy plays a distinctive role during some part of the species life history, opposing selective pressures could be balanced, permitting the stable co-existence of selfing and outcrossing sexual morphs. In this scenario, we would expect each morph to be specialised in their respective roles. Here we show, using behavioural, physiological and gene expression studies, that the selfing (hermaphrodite) and outcrossing (female) sexual morphs of the trioecious nematode Auanema freiburgensis have distinct adaptations optimised for their different roles during the life cycle. A. freiburgensis hermaphrodites are known to be produced under stressful conditions and are specialised for dispersal to new habitat patches. Here we show that they exhibit metabolic and intestinal changes enabling them to meet the cost of dispersal and reproduction. In contrast, A. freiburgensis females are produced in favourable conditions and facilitate rapid population growth. We found that females compensate for the lack of reproductive assurance by reallocating resources from intestinal development to mate-finding behaviour. The specialisation of each mating system for its role in the life cycle could balance opposing selective forces allowing the stable maintenance of both mating systems in A. freiburgensis.

Evidence of capacitation in the parasitoid wasp, Nasonia vitripennis, and its potential role in sex allocation.

The allocation of resources to the production of one sex or another has been observed in a large variety of animals. Its theoretical basis allows accurate predictions of offspring sex ratios in many species, but the mechanisms by which sex allocation is controlled are poorly understood. Using previously published data, we investigated whether alternative splicing, combined with differential gene expression, was involved with sex allocation in the parasitoid wasp, Nasonia vitripennis. We found that sex allocation is not controlled by alternative splicing but changes in gene and transcript-specific expression, which were identified to be involved with oviposition, were shown to be similar to those involved in sperm motility and capacitation. Genes involved in cholesterol efflux, a key component of capacitation, along with calcium transport, neurotransmission, trypsin, and MAPKinase activity were regulated in ovipositing wasps. The results show evidence for regulation of sperm motility and of capacitation in an insect which, in the context of the physiology of the N. vitripennis spermatheca, could be important for sex allocation.

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.

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.

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.

In silico Identification of a Molecular Circadian System With Novel Features in the Crustacean Model Organism Parhyale hawaiensis.

The amphipod Parhyale hawaiensis is a model organism of growing importance in the fields of evolutionary development and regeneration. A small, hardy marine crustacean that breeds year-round with a short generation time, it has simple lab culture requirements and an extensive molecular toolkit including the ability to generate targeted genetic mutant lines. Here we identify canonical core and regulatory clock genes using genomic and transcriptomic resources as a first step in establishing this species as a model in the field of chronobiology. The molecular clock of P. hawaiensis lacks orthologs of the canonical circadian genes cryptochrome 1 and timeless, in common with the mammalian system but in contrast to many arthropods including Drosophila melanogaster. Furthermore the predicted CLOCK peptide is atypical and CRY2 shows an extended 5' region of unknown function. These results appear to be shared by two other amphipod species.

Deleterious Mutation Accumulation in Arabidopsis thaliana Pollen Genes: A Role for a Recent Relaxation of Selection.

In many studies, sex-related genes have been found to evolve rapidly. We therefore expect plant pollen genes to evolve faster than sporophytic genes. In addition, pollen genes are expressed as haploids which can itself facilitate rapid evolution because recessive advantageous and deleterious alleles are not masked by dominant alleles. However, this mechanism is less straightforward to apply in the model plant species Arabidopsis thaliana. For 1 Myr, A. thaliana has been self-compatible, a life history switch that has caused: a reduction in pollen competition, increased homozygosity, and a dilution of masking in diploid expressed, sporophytic genes. In this study, we have investigated the relative strength of selection on pollen genes compared with sporophytic genes in A. thaliana. We present two major findings: 1) before becoming self-compatible, positive selection was stronger on pollen genes than sporophytic genes for A. thaliana and 2) current polymorphism data indicate that selection is weaker on pollen genes compared with sporophytic genes. This weaker selection on pollen genes can in part be explained by their higher tissue specificity, which in outbreeding plants can be outweighed by the effects of haploid expression and pollen competition. These results indicate that since A. thaliana has become self-compatible, selection on pollen genes has become more relaxed. This has led to higher polymorphism levels and a higher build-up of deleterious mutations in pollen genes compared with sporophytic genes.

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.

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.

An immune response in the bumblebee, Bombus terrestris leads to increased food consumption.

BACKGROUND: The concept of a costly immune system that must be traded off against other important physiological systems is fundamental to the burgeoning field of ecological immunity. Bumblebees have become one of the central models in this field. Although previous work has demonstrated costs of immunity in numerous life history traits, estimates of the more direct costs of bumblebee immunity have yet to be made. RESULTS: Here we show a 7.5% increase in energy consumption in response to non-pathogenic immune stimulation. CONCLUSION: This increase in energy consumption along with other results suggests that immunity is one of the most important physiological systems, with other systems being sacrificed for its continuing efficiency. This increased consumption and maintained activity contrasts with the sickness-induced anorexia and reduced activity found in vertebrates.

Acute and chronic gregarisation are associated with distinct DNA methylation fingerprints in desert locusts.

Desert locusts (Schistocerca gregaria) show a dramatic form of socially induced phenotypic plasticity known as phase polyphenism. In the absence of conspecifics, locusts occur in a shy and cryptic solitarious phase. Crowding with conspecifics drives a behavioural transformation towards gregariousness that occurs within hours and is followed by changes in physiology, colouration and morphology, resulting in the full gregarious phase syndrome. We analysed methylation-sensitive amplified fragment length polymorphisms (MS-AFLP) to compare the effect of acute and chronic crowding on DNA methylation in the central nervous system. We find that crowd-reared and solitary-reared locusts show markedly different neural MS-AFLP fingerprints. However, crowding for a day resulted in neural MS-AFLP fingerprints that were clearly distinct from both crowd-reared and uncrowded solitary-reared locusts. Our results indicate that changes in DNA methylation associated with behavioural gregarisation proceed through intermediate states that are not simply partial realisations of the endpoint states.

Allele specific expression in worker reproduction genes in the bumblebee Bombus terrestris.

Methylation has previously been associated with allele specific expression in ants. Recently, we found methylation is important in worker reproduction in the bumblebee Bombus terrestris. Here we searched for allele specific expression in twelve genes associated with worker reproduction in bees. We found allele specific expression in Ecdysone 20 monooxygenase and IMP-L2-like. Although we were unable to confirm a genetic or epigenetic cause for this allele specific expression, the expression patterns of the two genes match those predicted for imprinted genes.

Large-scale structure of brown rat (Rattus norvegicus) populations in England: effects on rodenticide resistance.

The brown rat (Rattus norvegicus) is a relatively recent (<300 years) addition to the British fauna, but by association with negative impacts on public health, animal health and agriculture, it is regarded as one of the most important vertebrate pest species. Anticoagulant rodenticides were introduced for brown rat control in the 1950s and are widely used for rat control in the UK, but long-standing resistance has been linked to control failures in some regions. One thus far ignored aspect of resistance biology is the population structure of the brown rat. This paper investigates the role population structure has on the development of anticoagulant resistance. Using mitochondrial and microsatellite DNA, we examined 186 individuals (from 15 counties in England and one location in Wales near the Wales-England border) to investigate the population structure of rural brown rat populations. We also examined individual rats for variations of the VKORC1 gene previously associated with resistance to anticoagulant rodenticides. We show that the populations were structured to some degree, but that this was only apparent in the microsatellite data and not the mtDNA data. We discuss various reasons why this is the case. We show that the population as a whole appears not to be at equilibrium. The relative lack of diversity in the mtDNA sequences examined can be explained by founder effects and a subsequent spatial expansion of a species introduced to the UK relatively recently. We found there was a geographical distribution of resistance mutations, and relatively low rate of gene flow between populations, which has implications for the development and management of anticoagulant resistance.