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

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

No abstract

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns, and worker discrimination behaviors, all evidently regulated by an inversion-based supergene that spans more than 13Mb of a "social chromosome," contains over 400 protein-coding genes, and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore effects of social chromosome genotype and natal colony social form on gene expression in virgin queens sampled as they embarked on nuptial flights, using RNA-sequencing of two important tissues. We observe relatively minor effects of natal social form, that is, of the social/developmental environment, on gene expression profiles, but substantial effects of genotype, including i) supergene-associated gene upregulation, ii) allele-specific expression, and iii) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following evolutionary divergence of the inversion-mediated Sb haplotype from its homolog that largely retained the ancestral gene order.The distinctive social chromosome-linked gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic underpinnings of the alternate social syndromes.

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion-based supergene that spans more than 13 Mb of a "social chromosome," contains over 400 protein-coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA-sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene-associated gene upregulation, (b) allele-specific expression and (c) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion-mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene-associated gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.