Transcriptome comparison between honey bee queen- and worker-destined larvae

The Juvenile-Hormone-Responsive Factor AmKr-h1 Regulates Caste Differentiation in Honey Bees

Honey bees are typical model organisms for the study of caste differentiation, and the juvenile hormone (JH) is a crucial link in the regulatory network of caste differentiation in honey bees. To investigate the mechanism of JH-mediated caste differentiation, we analyzed the effect of the JH response gene AmKr-h1 on this process. We observed that AmKr-h1 expression levels were significantly higher in queen larvae than in worker larvae at the 48 h, 84 h, and 120 h larval stages, and were regulated by JH. Inhibiting AmKr-h1 expression in honey bee larvae using RNAi could lead to the development of larvae toward workers. We also analyzed the transcriptome changes in honey bee larvae after AmKr-h1 RNAi and identified 191 differentially expressed genes (DEGs) and 682 differentially expressed alternative splicing events (DEASEs); of these, many were related to honey bee caste differentiation. Our results indicate that AmKr-h1 regulates caste differentiation in honey bees by acting as a JH-responsive gene.

A key gene for the climatic adaptation of Apis cerana populations in China according to selective sweep analysis

BACKGROUND

Apis cerana is widely distributed in China and, prior to the introduction of western honeybees, was the only bee species kept in China. During the long-term natural evolutionary process, many unique phenotypic variations have occurred among A. cerana populations in different geographical regions under varied climates. Understanding the molecular genetic basis and the effects of climate change on the adaptive evolution of A. cerana can promote A. cerana conservation in face of climate change and allow for the effective utilization of its genetic resources.

RESULT

To investigate the genetic basis of phenotypic variations and the impact of climate change on adaptive evolution, A. cerana workers from 100 colonies located at similar geographical latitudes or longitudes were analyzed. Our results revealed an important relationship between climate types and the genetic variation of A. cerana in China, and a greater influence of latitude compared with longitude was observed. Upon selection and morphometry analyses combination for populations under different climate types, we identified a key gene RAPTOR, which was deeply involved in developmental processes and influenced the body size.

CONCLUSION

The selection of RAPTOR at the genomic level during adaptive evolution could allow A. cerana to actively regulate its metabolism, thereby fine-tuning body sizes in response to harsh conditions caused by climate change, such as food shortages and extreme temperatures, which may partially elucidate the size differences of A. cerana populations. This study provides crucial support for the molecular genetic basis of the expansion and evolution of naturally distributed honeybee populations.

Detection and pH-Thermal Characterization of Proteinases Exclusive of Honeybee Worker-Fate Larvae (Apis mellifera L.)

The occurrence of the honeybee caste polyphenism arises when a change in diet is transduced into cellular metabolic responses, resulting in a developmental shift mediated by gene expression. The aim of this investigation was to detect and describe the expression profile of water-soluble proteases during the ontogenesis of honeybee worker-fate larvae. The extraction of insect homogenates was followed by the electrophoretic separation of the protein extract in polyacrylamide gels under semi-denaturing condition, precast with gelatin, pollen, or royal jelly protein extracts. The worker-fate honeybee larva showed a proteolytic pattern that varied with aging, and a protease with the highest activity at 72 h after hatching was named PS4. PS4 has a molecular weight of 45 kDa, it remained active until cell sealing, and its enzymatic properties suggest a serine-proteinase nature. To define the process that originates a queen-fate larvae, royal jelly and pollen were analysed, but PS4 was not detected in either of them. The effect of food on the PS4 was investigated by mixing crude extracts of queen and worker-fate larvae with pollen and royal jelly, respectively. Only royal jelly inhibited PS4 in worker-fate larvae. Taken together, our data suggest that PS4 could be involved in caste differentiation.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

OUP accepted manuscript

Many species have separate haploid and diploid phases. Theory predicts that each phase should experience the effects of evolutionary forces (like selection) differently. In the haploid phase, all fitness-affecting alleles are exposed to selection, whereas in the diploid phase, those same alleles can be masked by homologous alleles. This predicts that selection acting on genes expressed in haploids should be more effective than diploid-biased genes. Unfortunately, in arrhenotokous species, this prediction can be confounded with the effects of sex-specific expression, as haploids are usually reproductive males. Theory posits that, when accounting for ploidal- and sex-specific expression, selection should be equally efficient on haploid- and diploid-biased genes relative to constitutive genes. Here, we used a multiomic approach in honey bees to quantify the evolutionary rates of haploid-biased genes and test the relative effects of sexual- and haploid-expression on molecular evolution. We found that 16% of the honey bee’s protein-coding genome is highly expressed in haploid tissue. When accounting for ploidy and sex, haploid- and diploid-biased genes evolve at a lower rate than expected, indicating that they experience strong negative selection. However, the rate of molecular evolution of haploid-biased genes was higher than diploid-based genes. Genes associated with sperm storage are a clear exception to this trend with evidence of strong positive selection. Our results provide an important empirical test of theory outlining how selection acts on genes expressed in arrhenotokous species. We propose the haploid life history stage affects genome-wide patterns of diversity and divergence because of both sexual and haploid selection.

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.

Amsp3 may act upstream of Amdnmt3 in female caste differentiation in the honeybee (Apis mellifera)

In honey bees, the process of producing two female castes, including queens and workers, is nutritionally controlled by differential feeding royal jelly to newly emerged larvae. Although they have almost identical genetic blueprints, these castes show striking differences in their morphologies, longevities and reproductive capabilities. DNA methyltransferase 3 (Amdnmt3) gene is involved in the regulatory network for honeybee caste differentiation. Due to the role of two zinc fingers containing transcription factors, SP1 and SP3 in controlling mammalian Dnmts, this study aimed to determine a similar interaction of SPs with Amdnmt3 in the honeybee. We confirmed that the promoter region of Amdnmt3 contained multiple predicted SP1/SP3 binding sites and then investigated the role of AmSP3 in queen-worker differentiation network. We observed that the expression level of Amsp3 was significantly higher in worker larvae than that in queen larvae at 48 h, 84 h and 120 h. Knockdown of Amsp3 expression by RNAi in worker larvae significantly reduced the expression level of Amdnmt3 and caused morphological changes in adult bees towards a queen-like phenotype. However, the expression levels of Amsp3 and Amdnmt3 were repressed by juvenile hormone (JH). Our results suggest that AmSP3 is an important part of the queen-worker differentiation network and supports the role of Amdnmt3 in determining the phenotypic outcome of developing larvae.

Female developmental environment delays development of male honeybee (Apis mellifera)

Background

Nutrition and cell size play an important role in the determination of caste differentiation in queen and worker of honeybees (Apis mellifera), whereas the haploid genome dominates the differentiation of drones. However, the effects of female developmental environment on the development of males remain unclear. In this study, young drone larvae were transferred into worker cells (WCs) or remained in drone cells (DCs) to rear drones. The drone larvae were also grafted into queen cells (QCs) for 48 h and then transplanted into drone cells until emerging. Morphological indexes and reproductive organs of these three types of newly emerged drones were measured. Newly emerged drones and third instar drone larvae from WCs, DCs and QCs were sequenced by RNA sequencing (RNA-Seq).

Results

The amount of food remaining in cells of the QC and WC groups was significantly different to that in the DC group at the early larval stage. Morphological results showed that newly emerged DC drones had bigger body sizes and more well-developed reproductive tissues than WC and QC drones, whereas the reproductive tissues of QC drones were larger than those of WC drones. Additionally, whole body gene expression results showed a clear difference among three groups. At larval stage there were 889, 1761 and 1927 significantly differentially expressed genes (DEGs) in WC/DC, QC/DC and WC/QC comparisons, respectively. The number of DEGs decreased in adult drones of these three comparisons [678 (WC/DC), 338 (QC/DC) and 518 (WC/QC)]. A high number of DEGs were involved in sex differentiation, growth, olfaction, vision, mammalian target of rapamycin (mTOR), Wnt signaling pathways, and other processes.

Conclusions

This study demonstrated that the developmental environment of honeybee females can delay male development, which may serve as a model for understanding the regulation of sex differentiation and male development in social insects by environmental factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08014-1.

Age and Behavior-Dependent Differential miRNAs Expression in the Hypopharyngeal Glands of Honeybees (Apis mellifera L.)

This study aims to investigate the expression differences of miRNAs in the hypopharyngeal glands (HPGs) of honeybees at three developmental stages and to explore their regulation functions in the HPGs development. Small RNA sequencing was employed to analyze the miRNA profiles of HPGs in newly-emerged bees (NEB), nurse bees (NB), and forager bees (FB). Results showed that a total of 153 known miRNAs were found in the three stages, and ame-miR-276-3p, ame-miR-375-3p, ame-miR-14-3p, ame-miR-275-3p, and ame-miR-3477-5p were the top five most abundant ones. Furthermore, the expression of 11 miRNAs, 17 miRNAs, and 18 miRNAs were significantly different in NB vs. FB comparison, NB vs. NEB comparison, and in FB vs. NEB comparison, respectively, of which ame-miR-184-3p and ame-miR-252a-5p were downregulated in NB compared with that in both the FB and NEB, while ame-miR-11-3p, ame-miR-281-3p, and ame-miR-31a-5p had lower expression levels in FB compared with that in both the NB and NEB. Bioinformatic analysis showed that the potential target genes of the differentially expressed miRNAs (DEMs) were mainly enriched in several key signaling pathways, including mTOR signaling pathway, MAPK signaling pathway-fly, FoxO signaling pathway, Hippo signaling pathway-fly. Overall, our study characterized the miRNA profiles in the HPGs of honeybees at three different developmental stages and provided a basis for further study of the roles of miRNAs in HPGs development.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

The trophocytes and oenocytes of worker and queen honey bees (Apis mellifera) exhibit distinct age-associated transcriptome profiles

Despite the identical genomic context, trophocytes and oenocytes in worker bees exhibit aging-related phenotypes, in contrast to the longevity phenotypes in queen bees. To explore this phenomenon at the molecular level, we evaluated the age-associated transcriptomes of trophocytes and oenocytes in worker bees and queen bees using high-throughput RNA-sequencing technology (RNA-seq). The results showed that (i) while gene expression profiles were different between worker and queen bees, they remained similar between young and old counterparts; (ii) worker bees express a high proportion of low-abundance genes, whereas queen bee transcriptomes display a high proportion of moderate-expression genes; (iii) genes were upregulated to a greater extent in queen bees vs. worker bees; and (iv) distinct aging-related and longevity-related candidate genes were found in worker and queen bees. These results provide new insights into the cellular aging and longevity of trophocytes and oenocytes in honey bees. Identification of aging-associated biomarker genes also constitutes a basis for translational research of aging in higher organisms.

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.

High-quality chromosome-level genome assembly and full-length transcriptome analysis of the pharaoh ant Monomorium pharaonis

Background

Ants with complex societies have fascinated scientists for centuries. Comparative genomic and transcriptomic analyses across ant species and castes have revealed important insights into the molecular mechanisms underlying ant caste differentiation. However, most current ant genomes and transcriptomes are highly fragmented and incomplete, which hinders our understanding of the molecular basis for complex ant societies.

Findings

By hybridizing Illumina, Pacific Biosciences, and Hi-C sequencing technologies, we de novo assembled a chromosome-level genome for Monomorium pharaonis, with a scaffold N50 of 27.2 Mb. Our new assembly provides better resolution for the discovery of genome rearrangement events at the chromosome level. Analysis of full-length isoform sequencing (ISO-seq) suggested that ∼15 Gb of ISO-seq data were sufficient to cover most expressed genes, but the number of transcript isoforms steadily increased with sequencing data coverage. Our high-depth ISO-seq data greatly improved the quality of gene annotation and enabled the accurate detection of alternative splicing isoforms in different castes of M. pharaonis. Comparative transcriptome analysis across castes based on the ISO-seq data revealed an unprecedented number of transcript isoforms, including many caste-specific isoforms. We also identified a number of conserved long non-coding RNAs that evolved specifically in ant lineages and several that were conserved across insect lineages.

Conclusions

We produced a high-quality chromosome-level genome for M. pharaonis, which significantly improved previous short-read assemblies. Together with full-length transcriptomes for all castes, we generated a highly accurate annotation for this ant species. These long-read sequencing results provide a useful resource for future functional studies on the genetic mechanisms underlying the evolution of social behaviors and organization in ants.

Transgenerational accumulation of methylome changes discovered in commercially reared honey bee (Apis mellifera) queens

Whether a female honey bee (Apis mellifera) develops into a worker or a queen depends on her nutrition during development, which changes the epigenome to alter the developmental trajectory. Beekeepers typically exploit this developmental plasticity to produce queen bee by transplanting worker larvae into queen cells to be reared as queens, thus redirecting a worker developmental pathway to a queen developmental pathway. We studied the consequences of this manipulation for the queen phenotype and methylome over four generations. Queens reared from worker larvae consistently had fewer ovarioles than queens reared from eggs. Over four generations the methylomes of lines of queens reared from eggs and worker larvae diverged, accumulating increasing differences in exons of genes related to caste differentiation, growth and immunity. We discuss the consequences of these cryptic changes to the honey bee epigenome for the health and viability of honey bee stocks.

First Draft Genome Assembly of the Malaysian Stingless Bee, Heterotrigona itama (Apidae, Meliponinae)

The Malaysian stingless bee industry is hugely dependent on wild colonies. Nevertheless, the availability of new queens to establish new colonies is insufficient to meet the growing demand for hives in the industry. Heterotrigona itama is primarily utilized for honey production in the region and the major source of stingless bee colonies comes from the wild. To propagate new colonies domestically, a fundamental understanding of the biology of queen development, especially from the genomics aspect, is necessary. The whole genome was sequenced using a paired-end 150 strategy on the Illumina HiSeq X platform. The shotgun sequencing generated approximately 89 million raw pair-end reads with a total output of 13.37 Gb and a GC content of 37.31%. The genome size of the species was estimated to be approximately 272 Mb. Phylogenetic analysis showed H. itama are much more closely related to the bumble bee (Bombus spp.) than they are to the modern honey bee (Apis spp.). The genome data provided here are expected to contribute to a better understanding of the genetic aspect of queen differentiation as well as of important molecular pathways which are crucial for stingless bee biology, management and conservation.

Genome-Wide Differential DNA Methylation in Reproductive, Morphological, and Visual System Differences Between Queen Bee and Worker Bee (Apis mellifera)

There are many differences in external morphology and internal physiology between the Apis mellifera queen bee and worker bee, some of which are relevant to beekeeping production. These include reproductive traits, body size, royal jelly secreting properties, and visual system development, among others. The identification of candidate genes that control the differentiation of these traits is critical for selective honeybee breeding programs. In this study, we compared the genomic methylation of queen bee and worker bee larvae at 3, 4, and 5 days of age by whole-genome bisulfite sequencing, and found that the basic characteristics of genomic methylation in queen and worker larvae were the same. There were approximately 49 million cytosines in the Apis larvae genome, of which about 90,000 were methylated. Methylated CpG sites accounted for 99% of the methylated cytosines, and methylation mainly occurred in exons. However, methylation levels of queen and worker larvae showed different trends with age: the methylation level of queen larvae varied with age in an inverted parabola, while the corresponding trend for worker larvae with resembled an exponential curve with a platform. The methylation level of queen larvae was higher than that of worker larvae at 3 days of age, lower than that of worker larvae at 4 days of age, and similar to that of worker larvae at 5 days old. The top 10 differentially methylated genes (DMGs) and 13 caste-specific methylated genes were listed, and correlations with caste determination were speculated. We additionally screened 38 DMGs between queen larvae and worker larvae involved in specific organ differentiation as well as reproduction, morphology, and vision differentiation during caste determination. These genes are potential molecular markers for selective breeding of A. mellifera to improve fecundity, royal jelly production, body size, and foraging, and represent candidate genes for investigating specialized functional segregation during the process of caste differentiation.

Leveraging technological innovations to investigate evolutionary transitions to eusociality

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

A comparison of honeybee (Apis mellifera) queen, worker and drone larvae by RNA-Seq

Honeybees (Apis mellifera) have haplodiploid sex determination: males develop from unfertilized eggs and females develop from fertilized ones. The differences in larval food also determine the development of females. Here we compared the total somatic gene expression profiles of 2-day and 4-day-old drone, queen and worker larvae by RNA-Seq. The results from a co-expression network analysis on all expressed genes showed that 2-day-old drone and worker larvae were closer in gene expression profiles than 2-day-old queen larvae. This indicated that for young larvae (2-day-old) environmental factors such as larval diet have a greater effect on gene expression profiles than ploidy or sex determination. Drones had the most distinct gene expression profiles at the 4-day larval stage, suggesting that haploidy, or sex dramatically affects the gene expression of honeybee larvae. Drone larvae showed fewer differences in gene expression profiles at the 2-day and 4-day time points than the worker and queen larval comparisons (598 against 1190 and 1181), suggesting a different pattern of gene expression regulation during the larval development of haploid males compared to diploid females. This study indicates that early in development the queen caste has the most distinct gene expression profile, perhaps reflecting the very rapid growth and morphological specialization of this caste compared to workers and drones. Later in development the haploid male drones have the most distinct gene expression profile, perhaps reflecting the influence of ploidy or sex determination on gene expression.

Aedes aegypti HPX8C modulates immune responses against viral infection

Mosquitoes act as vectors of numerous pathogens that cause human diseases. Dengue virus (DENV) transmitted by mosquito, Aedes aegypti, is responsible for dengue fever epidemics worldwide with a serious impact on human health. Currently, disease control mainly relies on vector targeted intervention strategies. Therefore, it is imperative to understand the molecular mechanisms underlying the innate immune response of mosquitoes against pathogens. In the present study, the expression profiles of immunity-related genes in the midgut responding to DENV infection by feeding were analyzed by transcriptome and quantitative real-time PCR. The level of Antimicrobial peptides (AMPs) increased seven days post-infection (d.p.i.), which could be induced by the Toll immune pathway. The expression of reactive oxygen species (ROS) genes, including antioxidant genes, such as HPX7, HPX8A, HPX8B, HPX8C were induced at one d.p.i. and peaked again at ten d.p.i. in the midgut. Interestingly, down-regulation of the antioxidant gene HPX8C by RNA interference led to reduction in the virus titer in the mosquito, probably due to the elevated levels of ROS. Application of a ROS inhibitor and scavenger molecules further established the role of oxygen free radicals in the modulation of the immune response to DENV infection. Overall, our comparative transcriptome analyses provide valuable information about the regulation of immunity related genes in the transmission vector in response to DENV infection. It further allows us to identify novel molecular mechanisms underlying the host-virus interaction, which might aid in the development of novel strategies to control mosquito-borne diseases.

HPX8C is a heme-containing peroxidase, which can move reactive oxygen species (ROS) damage to the organism by reducing H₂O₂ to H₂O. Previously, the peroxidase gene has been shown to modulate midgut immunity and regulate anti-malarial response in mosquitoes. In this study, the classical immune signaling pathways, Toll and IMD genes might be late responses against the viruses. HPX8C was demonstrated here to play a role in antiviral immunity against DENV infection in Ae. Aegypti mosquitoes. HPX8C expression was induced by DENV infection and continued to increase with an elevated virus titer. In HPX8C-depleted mosquitoes, the ROS level was found to be increased with a corresponding decrease in the DENV and ZIKV virus titer. Therefore, it was speculated that HPX8C mediated immune responses against the DENV in the mosquito in the late stage of viral infection, which could be controlled by Toll pathway.

Morphological characterization and staging of bumble bee pupae

Bumble bees (Hymenoptera: Apidae, Bombus) are important pollinators and models for studying mechanisms underlying developmental plasticity, such as factors influencing size, immunity, and social behaviors. Research on such processes, as well as expanding use of gene-manipulation and gene expression technologies, requires a detailed understanding of how these bees develop. Developmental research often uses time-staging of pupae, however dramatic size differences in these bees can generate variation in developmental timing. To study developmental mechanisms in bumble bees, appropriate staging of developing bees using morphology is necessary. In this study, we describe morphological changes across development in several bumble bee species and use this to establish morphology-based staging criteria, establishing 20 distinct illustrated stages. These criteria, defined largely by eye and cuticle pigmentation patterns, are generalizable across members of the subgenus Pyrobombus, and can be used as a framework for study of other bumble bee subgenera. We examine the effects of temperature, caste, size, and species on pupal development, revealing that pupal duration shifts with each of these factors, confirming the importance of staging pupae based on morphology rather than age and the need for standardizing sampling.

Phenotypically distinct female castes in honey bees are defined by alternative chromatin states during larval development

The capacity of the honey bee to produce three phenotypically distinct organisms (two female castes; queens and sterile workers, and haploid male drones) from one genotype represents one of the most remarkable examples of developmental plasticity in any phylum. The queen-worker morphological and reproductive divide is environmentally controlled during post-embryonic development by differential feeding. Previous studies implicated metabolic flux acting via epigenetic regulation, in particular DNA methylation and microRNAs, in establishing distinct patterns of gene expression underlying caste-specific developmental trajectories. We produce the first genome-wide maps of chromatin structure in the honey bee at a key larval stage in which developmental canalization into queen or worker is virtually irreversible. We find extensive genome-wide differences in H3K4me3, H3K27ac, and H3K36me3, many of which correlate with caste-specific transcription. Furthermore, we identify H3K27ac as a key chromatin modification, with caste-specific regions of intronic H3K27ac directing the worker caste. These regions may harbor the first examples of caste-specific enhancer elements in the honey bee. Our results demonstrate a key role for chromatin modifications in the establishment and maintenance of caste-specific transcriptional programs in the honey bee. We show that at 96 h of larval growth, the queen-specific chromatin pattern is already established, whereas the worker determination is not, thus providing experimental support for the perceived timing of this critical point in developmental heterochrony in two types of honey bee females. In a broader context, our study provides novel data on environmentally regulated organismal plasticity and the molecular foundation of the evolutionary origins of eusociality.

Extensive Differential Splicing Underlies Phenotypically Plastic Aphid Morphs

Phenotypic plasticity results in a diversity of phenotypes from a single genotype in response to environmental cues. To understand the molecular basis of phenotypic plasticity, studies have focused on differential gene expression levels between environmentally determined phenotypes. The extent of alternative splicing differences among environmentally determined phenotypes has largely been understudied. Here, we study alternative splicing differences among plastically produced morphs of the pea aphid using RNA-sequence data. Pea aphids express two separate polyphenisms (plasticity with discrete phenotypes): a wing polyphenism consisting of winged and wingless females and a reproduction polyphenism consisting of asexual and sexual females. We find that pea aphids alternatively splice 34% of their genes, a high percentage for invertebrates. We also find that there is extensive use of differential spliced events between genetically identical, polyphenic females. These differentially spliced events are enriched for exon skipping and mutually exclusive exon events that maintain the open reading frame, suggesting that polyphenic morphs use alternative splicing to produce phenotype-biased proteins. Many genes that are differentially spliced between polyphenic morphs have putative functions associated with their respective phenotypes. We find that the majority of differentially spliced genes is not differentially expressed genes. Our results provide a rich candidate gene list for future functional studies that would not have been previously considered based solely on gene expression studies, such as ensconsin in the reproductive polyphenism, and CAKI in the wing polyphenism. Overall, this study suggests an important role for alternative splicing in the expression of environmentally determined phenotypes.

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

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

The Role of Brood in Eusocial Hymenoptera

Study of social traits in offspring traditionally reflects on interactions in simple family groups, with famous examples including parent-offspring conflict and sibling rivalry in birds and mammals. In contrast, studies of complex social groups such as the societies of ants, bees, and wasps focus mainly on adults and, in particular, on traits and interests of queens and workers. The social role of developing individuals in complex societies remains poorly understood. We attempt to fill this gap by illustrating that development in social Hymenoptera constitutes a crucial life stage with important consequences for the individual as well as the colony. We begin by describing the complex social regulatory network that modulates development in Hymenoptera societies. By highlighting the inclusive fitness interests of developing individuals, we show that they may differ from those of other colony members. We then demonstrate that offspring have evolved specialized traits that allow them to play a functional, cooperative role within colonies and give them the potential power to act toward increasing their inclusive fitness. We conclude by providing testable predictions for investigating the role of brood in colony interactions and giving a general outlook on what can be learned from studying offspring traits in hymenopteran societies.

A molecular concept of caste in insect societies

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

Transcriptomics of an extended phenotype: parasite manipulation of wasp social behaviour shifts expression of caste-related genes

Parasites can manipulate host behaviour to increase their own transmission and fitness, but the genomic mechanisms by which parasites manipulate hosts are not well understood. We investigated the relationship between the social paper wasp, Polistes dominula, and its parasite, Xenos vesparum (Insecta: Strepsiptera), to understand the effects of an obligate endoparasitoid on its host's brain transcriptome. Previous research suggests that X. vesparum shifts aspects of host social caste-related behaviour and physiology in ways that benefit the parasitoid. We hypothesized that X. vesparum-infested (stylopized) females would show a shift in caste-related brain gene expression. Specifically, we predicted that stylopized females, who would normally be workers, would show gene expression patterns resembling pre-overwintering queens (gynes), reflecting gyne-like changes in behaviour. We used RNA-sequencing data to characterize patterns of brain gene expression in stylopized females and compared these with those of unstylopized workers and gynes. In support of our hypothesis, we found that stylopized females, despite sharing numerous physiological and life-history characteristics with members of the worker caste, show gyne-shifted brain expression patterns. These data suggest that the parasitoid affects its host by exploiting phenotypic plasticity related to social caste, thus shifting naturally occurring social behaviour in a way that is beneficial to the parasitoid.

Genetics and Evolution of Social Behavior in Insects

The study of insect social behavior has offered tremendous insight into the molecular mechanisms mediating behavioral and phenotypic plasticity. Genomic applications to the study of eusocial insect species, in particular, have led to several hypotheses for the processes underlying the molecular evolution of behavior. Advances in understanding the genetic control of social organization have also been made, suggesting an important role for supergenes in the evolution of divergent behavioral phenotypes. Intensive study of social phenotypes across species has revealed that behavior and caste are controlled by an interaction between genetic and environmentally mediated effects and, further, that gene expression and regulation mediate plastic responses to environmental signals. However, several key methodological flaws that are hindering progress in the study of insect social behavior remain. After reviewing the current state of knowledge, we outline ongoing challenges in experimental design that remain to be overcome in order to advance the field.

Queen Quality and the Impact of Honey Bee Diseases on Queen Health: Potential for Interactions between Two Major Threats to Colony Health

Western honey bees, Apis mellifera, live in highly eusocial colonies that are each typically headed by a single queen. The queen is the sole reproductive female in a healthy colony, and because long-term colony survival depends on her ability to produce a large number of offspring, queen health is essential for colony success. Honey bees have recently been experiencing considerable declines in colony health. Among a number of biotic and abiotic factors known to impact colony health, disease and queen failure are repeatedly reported as important factors underlying colony losses. Surprisingly, there are relatively few studies on the relationship and interaction between honey bee diseases and queen quality. It is critical to understand the negative impacts of pests and pathogens on queen health, how queen problems might enable disease, and how both factors influence colony health. Here, we review the current literature on queen reproductive potential and the impacts of honey bee parasites and pathogens on queens. We conclude by highlighting gaps in our knowledge on the combination of disease and queen failure to provide a perspective and prioritize further research to mitigate disease, improve queen quality, and ensure colony health.

MicroRNAs Associated with Caste Determination and Differentiation in a Primitively Eusocial Insect

In eusocial Hymenoptera (ants, bees and wasps), queen and worker adult castes typically arise via environmental influences. A fundamental challenge is to understand how a single genome can thereby produce alternative phenotypes. A powerful approach is to compare the molecular basis of caste determination and differentiation along the evolutionary trajectory between primitively and advanced eusocial species, which have, respectively, relatively undifferentiated and strongly differentiated adult castes. In the advanced eusocial honeybee, Apis mellifera, studies suggest that microRNAs (miRNAs) play an important role in the molecular basis of caste determination and differentiation. To investigate how miRNAs affect caste in eusocial evolution, we used deep sequencing and Northern blots to isolate caste-associated miRNAs in the primitively eusocial bumblebee Bombus terrestris. We found that the miRNAs Bte-miR-6001-5p and -3p are more highly expressed in queen- than in worker-destined late-instar larvae. These are the first caste-associated miRNAs from outside advanced eusocial Hymenoptera, so providing evidence for caste-associated miRNAs occurring relatively early in eusocial evolution. Moreover, we found little evidence that miRNAs previously shown to be associated with caste in A. mellifera were differentially expressed across caste pathways in B. terrestris, suggesting that, in eusocial evolution, the caste-associated role of individual miRNAs is not conserved.

Linking Genes and Brain Development of Honeybee Workers: A Whole-Transcriptome Approach

Honeybees live in complex societies whose capabilities far exceed those of the sum of their single members. This social synergism is achieved mainly by the worker bees, which form a female caste. The worker bees display diverse collaborative behaviors and engage in different behavioral tasks, which are controlled by the central nervous system (CNS). The development of the worker brain is determined by the female sex and the worker caste determination signal. Here, we report on genes that are controlled by sex or by caste during differentiation of the worker's pupal brain. We sequenced and compared transcriptomes from the pupal brains of honeybee workers, queens and drones. We detected 333 genes that are differently expressed and 519 genes that are differentially spliced between the sexes, and 1760 genes that are differentially expressed and 692 genes that are differentially spliced between castes. We further found that 403 genes are differentially regulated by both the sex and caste signals, providing evidence of the integration of both signals through differential gene regulation. In this gene set, we found that the molecular processes of restructuring the cell shape and cell-to-cell signaling are overrepresented. Our approach identified candidate genes that may be involved in brain differentiation that ensures the various social worker behaviors.

Genome, transcriptome and methylome sequencing of a primitively eusocial wasp reveal a greatly reduced DNA methylation system in a social insect

Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect Polistes dominula, accompanied by analysis of caste-related transcriptome and methylome sequence data for adult queens and workers. Polistes dominula possesses a fairly typical hymenopteran genome, but shows very low genomewide GC content and some evidence of reduced genome size. We found numerous caste-related differences in gene expression, with evidence that both conserved and novel genes are related to caste differences. Most strikingly, these -omics data reveal a major reduction in one of the major epigenetic mechanisms that has been previously suggested to be important for caste differences in social insects: DNA methylation. Along with a conspicuous loss of a key gene associated with environmentally responsive DNA methylation (the de novo DNA methyltransferase Dnmt3), these wasps have greatly reduced genomewide methylation to almost zero. In addition to providing a valuable resource for comparative analysis of social insect evolution, our integrative -omics data for this important behavioural and evolutionary model system call into question the general importance of DNA methylation in caste differences and evolution in social insects.

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.

MicroRNAs in Honey Bee Caste Determination

The cellular mechanisms employed by some organisms to produce contrasting morphological and reproductive phenotypes from the same genome remains one of the key unresolved issues in biology. Honeybees (Apis mellifera) use differential feeding and a haplodiploid sex determination system to generate three distinct organismal outcomes from the same genome. Here we investigate the honeybee female and male caste-specific microRNA and transcriptomic molecular signatures during a critical time of larval development. Both previously undetected and novel miRNAs have been discovered, expanding the inventory of these genomic regulators in invertebrates. We show significant differences in the microRNA and transcriptional profiles of diploid females relative to haploid drone males as well as between reproductively distinct females (queens and workers). Queens and drones show gene enrichment in physio-metabolic pathways, whereas workers show enrichment in processes associated with neuronal development, cell signalling and caste biased structural differences. Interestingly, predicted miRNA targets are primarily associated with non-physio-metabolic genes, especially neuronal targets, suggesting a mechanistic disjunction from DNA methylation that regulates physio-metabolic processes. Accordingly, miRNA targets are under-represented in methylated genes. Our data show how a common set of genetic elements are differentially harnessed by an organism, which may provide the remarkable level of developmental flexibility required.

Developmental Transcriptome for a Facultatively Eusocial Bee, Megalopta genalis

Transcriptomes provide excellent foundational resources for mechanistic and evolutionary analyses of complex traits. We present a developmental transcriptome for the facultatively eusocial bee Megalopta genalis, which represents a potential transition point in the evolution of eusociality. A de novo transcriptome assembly of Megalopta genalis was generated using paired-end Illumina sequencing and the Trinity assembler. Males and females of all life stages were aligned to this transcriptome for analysis of gene expression profiles throughout development. Gene Ontology analysis indicates that stage-specific genes are involved in ion transport, cell-cell signaling, and metabolism. A number of distinct biological processes are upregulated in each life stage, and transitions between life stages involve shifts in dominant functional processes, including shifts from transcriptional regulation in embryos to metabolism in larvae, and increased lipid metabolism in adults. We expect that this transcriptome will provide a useful resource for future analyses to better understand the molecular basis of the evolution of eusociality and, more generally, phenotypic plasticity.

How Many Genes Are Expressed in a Transcriptome? Estimation and Results for RNA-Seq

RNA-seq experiments estimate the number of genes expressed in a transcriptome as well as their relative frequencies. However, an undetermined number of genes can remain undetected due to their low expression relative to the sample size (sequence depth). Estimation of the true number of genes expressed in a transcriptome is essential in order to determine which genes are exclusively expressed in specific tissues or under particular conditions. A reliable estimate of the true number of expressed genes is also required to accurately measure transcriptome changes and to predict the sequencing depth needed to increase the proportion of detected genes. This problem is analogous to ecological sampling problems such as estimating the number of species at a given site. Here we present a non-parametric estimator for the number of undetected genes as well as for the extra sample size needed to detect a given proportion of the undetected genes. Our estimators are superior to ones already published by having smaller standard errors and biases. We applied our method to a set of 32 publicly available RNA-seq experiments, including the evaluation of 311 individually sequenced libraries. We found that in the majority of the cases more than one thousand genes are undetected, and that on average approximately 6% of the expressed genes per accession remain undetected. This figure increases to approximately 10% if individual sequencing libraries are analyzed. Our method is also applicable to metagenomic experiments. Using our method, the number of undetected genes as well as the sample size needed to detect them can be calculated, leading to more accurate and complete gene expression studies.

A Critical Role for CLSP2 in the Modulation of Antifungal Immune Response in Mosquitoes

Entomopathogenic fungi represent a promising class of bio-insecticides for mosquito control. Thus, detailed knowledge of the molecular mechanisms governing anti-fungal immune response in mosquitoes is essential. In this study, we show that CLSP2 is a modulator of immune responses during anti-fungal infection in the mosquito Aedes aegypti. With a fungal infection, the expression of the CLSP2 gene is elevated. CLSP2 is cleaved upon challenge with Beauveria bassiana conidia, and the liberated CLSP2 CTL-type domain binds to fungal cell components and B. bassiana conidia. Furthermore, CLPS2 RNA interference silencing significantly increases the resistance to the fungal challenge. RNA-sequencing transcriptome analysis showed that the majority of immune genes were highly upregulated in the CLSP2-depleted mosquitoes infected with the fungus. The up-regulated immune gene cohorts belong to melanization and Toll pathways, but not to the IMD or JAK-STAT. A thioester-containing protein (TEP22), a member of α2-macroglobulin family, has been implicated in the CLSP2-modulated mosquito antifungal defense. Our study has contributed to a greater understanding of immune-modulating mechanisms in mosquitoes.

Climbing the social ladder: the molecular evolution of sociality

Genomic tools are allowing us to dissect the roles of genes and genetic architecture in social evolution, and eusocial insects are excellent models. Numerous hypotheses for molecular evolution of eusociality have been proposed, ranging from regulatory shifts in 'old' genes to rapid evolution of 'new' genes. A broad model to explain this major transition in evolution has been lacking. We provide a synthetic framework centered on the idea that different evolutionary processes dominate during different transitional stages, beginning with changes in gene regulation and culminating in novel genes later on. By considering multiple mechanisms as we 'climb the social ladder', we can test whether the transitions from solitary to simple sociality to complex sociality represent incremental changes or genetic revolutions.

Sphingolipids, Transcription Factors, and Conserved Toolkit Genes: Developmental Plasticity in the Ant Cardiocondyla obscurior

Developmental plasticity allows for the remarkable morphological specialization of individuals into castes in eusocial species of Hymenoptera. Developmental trajectories that lead to alternative caste fates are typically determined by specific environmental stimuli that induce larvae to express and maintain distinct gene expression patterns. Although most eusocial species express two castes, queens and workers, the ant Cardiocondyla obscurior expresses diphenic females and males; this provides a unique system with four discrete phenotypes to study the genomic basis of developmental plasticity in ants. We sequenced and analyzed the transcriptomes of 28 individual C. obscurior larvae of known developmental trajectory, providing the first in-depth analysis of gene expression in eusocial insect larvae. Clustering and transcription factor binding site analyses revealed that different transcription factors and functionally distinct sets of genes are recruited during larval development to induce the four alternative trajectories. In particular, we found complex patterns of gene regulation pertaining to sphingolipid metabolism, a conserved molecular pathway involved in development, obesity, and aging.

Nourishment level affects caste-related gene expression in Polistes wasps

Background

Social insects exhibit striking phenotypic plasticity in the form of distinct reproductive (queen) and non-reproductive (worker) castes, which are typically driven by differences in the environment during early development. Nutritional environment and nourishment during development has been shown to be broadly associated with caste determination across social insect taxa such as bees, wasps, and termites. In primitively social insects such as Polistes paper wasps, caste remains flexible throughout adulthood, but there is evidence that nourishment inequalities can bias caste development with workers receiving limited nourishment compared to queens. Dominance and vibrational signaling are behaviors that have also been linked to caste differences in paper wasps, suggesting that a combination of nourishment and social factors may drive caste determination. To better understand the molecular basis of nutritional effects on caste determination, we used RNA-sequencing to investigate the gene expression changes in response to proteinaceous nourishment deprivation in Polistes metricus larvae.

Results

We identified 285 nourishment-responsive transcripts, many of which are related to lipid metabolism and oxidation-reduction activity. Via comparisons to previously identified caste-related genes, we found that nourishment restriction only partially biased wasp gene expression patterns toward worker caste-like traits, which supports the notion that nourishment, in conjunction with social environment, is a determinant of developmental caste bias. In addition, we conducted cross-species comparisons of nourishment-responsive genes, and uncovered largely lineage-specific gene expression changes, suggesting few shared nourishment-responsive genes across taxa.

Conclusion

Overall, the results from this study highlight the complex and multifactorial nature of environmental effects on the gene expression patterns underlying plastic phenotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1410-y) contains supplementary material, which is available to authorized users.

Transcriptomic study of the red palm weevil Rhynchophorus ferrugineus embryogenesis

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-β, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

Genetic architecture of a hormonal response to gene knockdown in honey bees

Variation in endocrine signaling is proposed to underlie the evolution and regulation of social life histories, but the genetic architecture of endocrine signaling is still poorly understood. An excellent example of a hormonally influenced set of social traits is found in the honey bee (Apis mellifera): a dynamic and mutually suppressive relationship between juvenile hormone (JH) and the yolk precursor protein vitellogenin (Vg) regulates behavioral maturation and foraging of workers. Several other traits cosegregate with these behavioral phenotypes, comprising the pollen hoarding syndrome (PHS) one of the best-described animal behavioral syndromes. Genotype differences in responsiveness of JH to Vg are a potential mechanistic basis for the PHS. Here, we reduced Vg expression via RNA interference in progeny from a backcross between 2 selected lines of honey bees that differ in JH responsiveness to Vg reduction and measured JH response and ovary size, which represents another key aspect of the PHS. Genetic mapping based on restriction site-associated DNA tag sequencing identified suggestive quantitative trait loci (QTL) for ovary size and JH responsiveness. We confirmed genetic effects on both traits near many QTL that had been identified previously for their effect on various PHS traits. Thus, our results support a role for endocrine control of complex traits at a genetic level. Furthermore, this first example of a genetic map of a hormonal response to gene knockdown in a social insect helps to refine the genetic understanding of complex behaviors and the physiology that may underlie behavioral control in general.

Transcriptome comparison between inactivated and activated ovaries of the honey bee Apis mellifera L

Ovarian activity not only influences fertility, but is also involved with the regulation of division of labour between reproductive and behavioural castes of female honey bees. In order to identify candidate genes associated with ovarian activity, we compared the gene expression patterns between inactivated and activated ovaries of queens and workers by means of high-throughput RNA-sequencing technology. A total of 1615 differentially expressed genes (DEGs) was detected between ovaries of virgin and mated queens, and more than 5300 DEGs were detected between inactivated and activated worker ovaries. Intersection analysis of DEGs amongst five libraries revealed that a similar set of genes (824) participated in the ovary activation of both queens and workers. A large number of these DEGs were predominantly related to cellular, cell and cell part, binding, biological regulation and metabolic processes. In addition, over 1000 DEGs were linked to more than 230 components of Kyoto Encyclopedia of Genes and Genomes pathways, including 25 signalling pathways. The reliability of the RNA-sequencing results was confirmed by means of quantitative real-time PCR. Our results provide new insights into the molecular mechanisms involved in ovary activation and reproductive division of labour.

Proteomics analysis reveals protein expression differences for hypopharyngeal gland activity in the honeybee, Apis mellifera carnica Pollmann

BACKGROUND

Most of the proteins contained in royal jelly (RJ) are secreted from the hypopharyngeal glands (HG) of young bees. Although generic protein composition of RJ has been investigated, little is known about how age-dependent changes on HG secretion affect RJ composition and their biological consequences. In this study, we identified differentially expressed proteins (DEPs) during HG development by using the isobaric tag for relative and absolute quantification (iTRAQ) labeling technique. This proteomic method increases the potential for new protein discovery by improving the identification of low quantity proteins.

RESULTS

A total of 1282 proteins were identified from five age groups of worker bees, 284 of which were differentially expressed. 43 (15.1%) of the DEPs were identified for the first time. Comparison of samples at day 6, 9, 12, and 16 of development relative to day 3 led to the unambiguous identification of 112, 117, 127, and 127 DEPs, respectively. The majority of these DEPs were up-regulated in the older worker groups, indicating a substantial change in the pattern of proteins expressed after 3 days. DEPs were identified among all the age groups, suggesting that changes in protein expression during HG ontogeny are concomitant with different states of worker development. A total of 649 proteins were mapped to canonical signaling pathways found in the Kyoto Encyclopedia of Genes and Genomes (KEGG), which were preferentially associated with metabolism and biosynthesis of secondary metabolites. More than 10 key high-abundance proteins were involved in signaling pathways related to ribosome function and protein processing in the endoplasmic reticulum. The results were validated by qPCR.

CONCLUSION

Our approach demonstrates that HG experienced important changes in protein expression during its ontogenic development, which supports the secretion of proteins involved in diverse functions in adult workers beyond its traditional role in royal jelly production.

De novo transcriptome assembly from fat body and flight muscles transcripts to identify morph-specific gene expression profiles in Gryllus firmus

Wing polymorphism is a powerful model for examining many aspects of adaptation. The wing dimorphic cricket species, Gryllus firmus, consists of a long-winged morph with functional flight muscles that is capable of flight, and two flightless morphs. One (obligately) flightless morph emerges as an adult with vestigial wings and vestigial flight muscles. The other (plastic) flightless morph emerges with fully-developed wings but later in adulthood histolyzes its flight muscles. Importantly both flightless morphs have substantially increased reproductive output relative to the flight-capable morph. Much is known about the physiological and biochemical differences between the morphs with respect to adaptations for flight versus reproduction. In contrast, little is known about the molecular genetic basis of these morph-specific adaptations. To address this issue, we assembled a de novo transcriptome of G. firmus using 141.5 million Illumina reads generated from flight muscles and fat body, two organs that play key roles in flight and reproduction. We used the resulting 34,411 transcripts as a reference transcriptome for differential gene expression analyses. A comparison of gene expression profiles from functional flight muscles in the flight-capable morph versus histolyzed flight muscles in the plastic flight incapable morph identified a suite of genes involved in respiration that were highly expressed in pink (functional) flight muscles and genes involved in proteolysis highly expressed in the white (histolyzed) flight muscles. A comparison of fat body transcripts from the obligately flightless versus the flight-capable morphs revealed differential expression of genes involved in triglyceride biosynthesis, lipid transport, immune function and reproduction. These data provide a valuable resource for future molecular genetics research in this and related species and provide insight on the role of gene expression in morph-specific adaptations for flight versus reproduction.