Hexamerin-based regulation of juvenile hormone-dependent gene expression underlies phenotypic plasticity in a social insect

Transcriptomic correlates of nutritional manipulation in a facultatively social bee

Subsocial behaviour in insects consists of extended parental care and may set the stage for the evolution of cooperation through manipulation of offspring. Manipulation of brood nutrition may produce differences in developmental or adult gene regulation, but it also produces smaller offspring which may be coerced into cooperation. The eastern small carpenter bee Ceratina calcarata frequently produces a smaller under-provisioned dwarf eldest daughter (DED). These DEDs are the only offspring to forage and feed siblings. To test whether nutritional manipulation of DEDs alters gene expression, inducing cooperative sibling care, we conducted a transcriptomic study, using whole heads, to assess differences in brain gene expression among naturally provisioned regular daughters and DEDs, experimentally under-provisioned regular daughters, and experimentally supplemented DEDs, prior to social interaction. Differences in gene expression were minimal among groups but were dramatic as a function of body size as a continuous variable, suggesting that differences in gene expression are more associated with absolute differences in body size, not discrete castes or order of eclosion. Enrichment for GO terms related to hormonal regulation in small bees points to hormonal regulation of transcription factors in behavioural differences that emerge in DEDs. Subordinate behaviours thus likely involve experience and social environment, though other developmental mechanisms, such as parental care, and later adult social interactions after eclosion, may act on differences in body size and gene expression to produce the distinct behaviour of DEDs.

The proteomic content of Varroa destructor gut varies according to the developmental stage of its host

The nutritional physiology of parasites is often overlooked although it is at the basis of host-parasite interactions. In the case of Varroa destructor, one of the major pests of the Western honey bee Apis mellifera, the nature of molecules and tissues ingested by the parasite is still not completely understood. Here, the V. destructor feeding biology was explored through artificial feeding, dissection of the mite's gut and proteomic analyses. More specifically, the proteome of guts extracted from starved mites and honey bee-fed mites was compared to highlight both the parasite proteins likely involved in food processing and the honey bee proteins actually ingested by the mite. We could identify 25 V. destructor candidate proteins likely involved in the parasite digestion. As the host developmental stages infested by the mite are diverse, we also focused on the identity and on the origin of honey bee proteins ingested by the mite when it feeds on larvae, pupae or adults. We highlighted profiles of consumed honey bee proteins and their variations throughout the V. destructor life cycle. These variations matched the ones observed in the honey bee hemolymph, showing that this tissue is an important part of the mite's diet. Based on the variations of abundance of the most consumed honey bee proteins and on their functions, the potential implication of these key candidate nutrients in V. destructor reproduction is also discussed.

Body part-specific development in termite caste differentiation: crosstalk between hormonal actions and developmental toolkit genes

In social insects, interactions among colony members trigger caste differentiation with morphological modifications. During caste differentiation in termites, body parts and caste-specific morphologies are modified during postembryonic development under endocrine controls such as juvenile hormone (JH) and ecdysone. In addition to endocrine factors, developmental toolkit genes such as Hox- and appendage-patterning genes also contribute to the caste-specific body part modifications. These toolkits are thought to provide spatial information for specific morphogenesis. During social evolution, the complex crosstalks between physiological and developmental mechanisms should be established, leading to the sophisticated caste systems. This article reviews recent studies on these mechanisms underlying the termite caste differentiation and addresses implications for the evolution of caste systems in termites.

Nutrition- and hormone-controlled developmental plasticity in Blattodea

Blattodea, which includes cockroaches and termites, possesses high developmental plasticity that is mainly controlled by nutritional conditions and insect hormones. Insulin/insulin-like growth factor signaling (IIS), target of rapamycin complex 1 (TORC1), and adenosine monophosphate-activated protein complex are the three primary nutrition-responsive signals. Juvenile hormone (JH) and 20-hydroxyecdysone (20E) constitute the two most vital insect hormones that might interact with each other through the Met, Kr-h1, E93 (MEKRE93) pathway. Nutritional and hormonal signals interconnect to create a complex regulatory network. Here we summarize recent progress in our understanding of how nutritional and hormonal signals coordinately control the developmental plasticity of metamorphosis, reproduction, and appendage regeneration in cockroaches as well as caste differentiation in termites. We also highlight several perspectives that should be further emphasized in the studies of developmental plasticity in Blattodea. This review provides a general landscape in the field of nutrition- and hormone-controlled developmental plasticity in insects.

Brain Gene Expression of Foraging Behavior and Social Environment in Ceratina calcarata

Rudimentary social systems have the potential to both advance our understanding of how complex sociality may have evolved and our understanding of how changes in social environment may influence gene expression and cooperation. Recently, studies of primitively social Hymenoptera have greatly expanded empirical evidence for the role of social environment in shaping behavior and gene expression. Here, we compare brain gene expression profiles of foragers across social contexts in the small carpenter bee, Ceratina calcarata. We conducted experimental manipulations of field colonies to examine gene expression profiles among social contexts including foraging mothers, regular daughters, and worker-like dwarf eldest daughters in the presence and absence of mother. Our analysis found significant differences in gene expression associated with female age, reproductive status, and social environment, including circadian clock gene dyw, hexamerin, and genes involved in the regulation of juvenile hormone and chemical communication. We also found that candidate genes differentially expressed in our study were also associated with division of labor, including foraging, in other primitively and advanced eusocial insects. Our results offer evidence for the role of the regulation of key developmental hormones and circadian rhythms in producing cooperative behavior in rudimentary insect societies.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

Fat enough for the winter? Does nutritional status affect diapause?

Many insects enter a dormant state termed diapause in anticipation of seasonal inhospitable conditions. Insects drastically reduce their feeding during diapause. Their reduced nutrient intake is paired with substantial nutrient costs: maintaining basal metabolism during diapause, repairing tissues damaged by adverse conditions, and resuming development after diapause. Many investigators have asked "Does nutrition affect diapause?" In this review, we survey the studies that have attempted to address this question. We propose the term nutritional status, a holistic view of nutrition that explicitly includes the perception, intake, and storage of the great breadth of nutrients. We examine the studies that have sought to test if nutrition affects diapause, trying to identify specific facets of nutritional status that affect diapause phenotypes. Curiously, low quality host plants during the diapause induction phase generally induce diapause, but food deprivation during the same phase generally averts diapause. Using the geometric framework of nutrition to identify specific dietary components that affect diapause may reconcile these contrasting findings. This framework can establish nutritionally permissive space, distinguishing nutrient changes that affect diapause from changes that induce other dormancies. Refeeding is another important experimental technique that distinguishes between diapause and quiescence, a non-diapause dormancy. We also find insufficient evidence for the hypothesis that nutrient stores regulate diapause length and suggest manipulations to investigate the role of nutrient stores in diapause termination. Finally, we propose mechanisms that could interface nutritional status with the diapause program, focusing on combined action of the nutritional axis between the gut, fat body, and brain.

Myofilaments promote wing expansion and maintain genitalia morphology in the American cockroach, Periplaneta americana

Insects are the most widely distributed and successful animals on the planet. A large number of insects are capable of flight with functional wings. Wing expansion is an important process for insects to achieve functional wings after eclosion and healthy genital morphology is crucial for adult reproduction. Myofilaments are functional units that constitute sarcomeres and trigger muscle contraction. Here, we identified four myofilament proteins, including Myosin, Paramyosin, Tropomyosin and Troponin T, from the wing pads of nymphs in the American cockroach, Periplaneta americana. RNAi-mediated knockdown of Myosin, Paramyosin, Tropomyosin and Troponin T in the early stage of final instar nymphs caused a severely curly wing phenotype in the imaginal moult, especially in the Paramyosin and Troponin T knockdown groups, indicating that these myofilament proteins are involved in controlling wing expansion behaviours during the nymph-adult transition. In addition, the knockdown resulted in abnormal external genitalia, caused ovulation failure, and affected male accessory gland development. Interestingly, the expression of myofilament genes was induced by methoprene, a juvenile hormone (JH) analogue, and decreased by the depletion of the JH receptor gene Met. Altogether, we have determined that myofilament genes play an important role in promoting wing expansion and maintaining adult genitalia morphology, and their expression is induced by JH signalling. Our data reveal a novel mechanism by which wing expansion is regulated by myofilaments and the functions of myofilaments are involved in maintaining genitalia morphology.

Protein storage and reproduction increase in grasshoppers on a diet matched to the amino acids of egg yolk protein

The diets of animals are essential to support development, and protein is key. Accumulation of stored nutrients can support developmental events such as molting and initiation of reproduction. Agricultural studies have addressed how dietary protein quality affects growth, but few studies have addressed the effects of dietary protein quality on developmental transitions. Studies on how dietary quality may affect protein storage and development are possible in arthropods, which store proteins in the hemolymph. We hypothesized that diets with a composition of amino acids that matches the precursor of egg yolk protein (vitellogenin, Vg) will be high quality and support both egg production and accumulation of storage proteins. Grasshoppers were fed one of two isonitrogenous solutions of amino acids daily: Vg-balanced (matched to Vg) or Unbalanced (same total moles of amino acids, but not matched to egg yolk). We measured reproduction and storage protein levels in serial hemolymph samples from individuals. The Vg-balanced group had greater reproduction and greater cumulative levels of storage proteins than did the Unbalanced group. This occurred even though amino acids fed to the Vg-balanced group were not a better match to storage protein than were the amino acids fed to the Unbalanced group. Further, oviposition timing was best explained by a combination of diet, age at the maximum level of storage protein hexamerin-270 and accumulation of hexamerin-90. Our study tightens the link between storage proteins and commitment to reproduction, and shows that dietary protein quality is vital for protein storage and reproduction.

Drone honey bees are disproportionately sensitive to abiotic stressors despite expressing high levels of stress response proteins

Drone honey bees (Apis mellifera) are the obligate sexual partners of queens, and the availability of healthy, high-quality drones directly affects a queen's fertility and productivity. Yet, our understanding of how stressors affect adult drone fertility, survival, and physiology is presently limited. Here, we investigated sex biases in susceptibility to abiotic stressors (cold stress, topical imidacloprid exposure, and topical exposure to a realistic cocktail of pesticides). We found that drones (haploid males) were more sensitive to cold and imidacloprid exposure than workers (sterile, diploid females), but the cocktail was not toxic at the concentrations tested. We corroborated this lack of cocktail toxicity with in-hive exposures via pollen feeding. We then used quantitative proteomics to investigate protein expression profiles in the hemolymph of topically exposed workers and drones, and found that 34 proteins were differentially expressed in exposed drones relative to controls, but none were differentially expressed in exposed workers. Contrary to our hypothesis, we show that drones express surprisingly high baseline levels of putative stress response proteins relative to workers. This suggests that drones' stress tolerance systems are fundamentally rewired relative to workers, and susceptibility to stress depends on more than simply gene dose or allelic diversity.

miRNA-34 and miRNA-210 target hexamerin genes enhancing their differential expression during early brain development of honeybee (Apis mellifera) castes

During the honeybee larval stage, queens develop larger brains than workers, with morphological differentiation appearing at the fourth larval phase (L4), just after a boost in nutritional difference both prospective females experience. The molecular promoters of this caste-specific brain development are already ongoing in previous larval phases. Transcriptomic analyses revealed a set of differentially expressed genes in the L3 brains of queens and workers, which represents the early molecular response to differential feeding females receive during larval development. Three genes of this set, hex70b, hex70c and hex110, are more highly transcribed in the brain of workers than in queens. The microRNAs miR-34, miR-210 and miR-317 are in higher levels in the queens' brain at the same phase of larval development. Here, we tested the hypothesis that the brain of workers expresses higher levels of hexamerins than that of queens during key phases of larval development and that this differential hexamerin genes expression is further enhanced by the repressing activity of miR-34, miR-210 and miR-317. Our transcriptional analyses showed that hex70b, hex70c and hex110 genes are differentially expressed in the brain of L3 and L4 larval phases of honeybee queens and workers. In silico reconstructed miRNA-mRNA interaction networks were validated using luciferase assays, which showed miR-34 and miR-210 negatively regulate hex70b and hex110 genes by directly and redundantly binding their 3'UTR (untranslated region) sequences. Taken together, our results suggest that miR-34 and miR-210 act together promoting differential brain development in honeybee castes by downregulating the expression of the putative antineurogenic hexamerin genes hex70b and hex110.

Regulation of host phenotypic plasticity by gut symbiont communities in the eastern subterranean termite (Reticulitermes flavipes)

Termites are eusocial insects that host a range of prokaryotic and eukaryotic gut symbionts and can differentiate into a range of caste phenotypes. Soldier caste differentiation from termite workers follows two successive molts (worker-presoldier-soldier) that are driven at the endocrine level by juvenile hormone (JH). Although physiological and eusocial mechanisms tied to JH signaling have been studied, the role of gut symbionts in the caste differentiation process is poorly understood. Here, we used the JH analog methoprene in combination with the antibiotic kanamycin to manipulate caste differentiation and gut bacterial loads in Reticulitermes flavipes termites via four bioassay treatments: kanamycin, methoprene, kanamycin+methoprene, and an untreated (negative) control. Bioassay results demonstrated a significantly higher number of presoldiers in the methoprene treatment, highest mortality in kanamycin+methoprene treatment, and significantly reduced protist numbers in all treatments except the untreated control. Bacterial 16S rRNA gene sequencing provided alpha and beta diversity results that mirrored bioassay findings. From ANCOM analysis, we found that several bacterial genera were differentially abundant among treatments. Finally, follow-up experiments showed that if methoprene and kanamycin or untreated termites are placed together, zero or rescued presoldier initiation, respectively, occurs. These findings reveal that endogenous JH selects for symbiont compositions required to successfully complete presoldier differentiation. However, if the gut is voided before the influx of JH, it cannot select for the necessary symbionts that are crucial for molting. Based on these results, we are able to provide a novel example of linkages between gut microbial communities and host phenotypic plasticity.

Characterization and Developmental Expression Patterns of Four Hexamerin Genes in the Bumble Bee, Bombus terrestris (Hymenoptera: Apidae)

Hexamerins are members of the hemocyanin superfamily and play essential roles in providing amino acids and energy for the nonfeeding stages of insects. In this study, we cloned and analyzed the expression patterns of four hexamerin genes (hex 70a, hex 70b, hex 70c, and hex 110) at different worker development stages and queen diapause statuses in the bumble bee, Bombus terrestris. The results of this study showed that hex 110 has the longest open reading frame (ORF; 3,297 bp) compared to the ORFs of hex 70a (2,034 bp), hex 70b (2,067 bp), and hex 70c (2,055 bp). The putative translation product of Hex 70a, Hex 70b, Hex70c, and Hex 110 has 677, 688, 684, and 1,098aa with predicted molecular mass of 81.13, 79.69, 81.58, and 119 kDa. In the development stages of workers, the expression levels of hex 70a, hex 70b, and hex 70c increased gradually from the larval stage and exhibited high expression levels at the pink eyed and brown eyed pupae stage, whereas hex 110 exhibited the highest expression level at the larval period. Four hexamerin genes were highly expressed at the prediapause status of queen (P < 0.05), and compared to the eclosion queen, the lowest upregulation was 3.7-fold, and the highest upregulation was 1,742-fold. The expression levels of hex 70b, hex 70c, and hex 110 at diapause were significantly higher than those at postdiapause (P < 0.05). In conclusion, hexamerins may play important roles in queen diapause and metamorphosis of larval and pupal stages.

Molecular underpinnings of the early brain developmental response to differential feeding in the honey bee Apis mellifera

Brain differential morphogenesis in females is one of the major phenotypic manifestations of caste development in honey bees. Brain diphenism appears at the fourth larval phase as a result of the differential feeding regime developing females are submitted during early phases of larval development. Here, we used a forward genetics approach to test the early brain molecular response to differential feeding leading to the brain diphenism observed at later developmental phases. Using RNA sequencing analysis, we identified 53 differentially expressed genes (DEGs) between the brains of queens and workers at the third larval phase. Since miRNAs have been suggested to play a role in caste differentiation after horizontal and vertical transmission, we tested their potential participation in regulating the DEGs. The miRNA-mRNA interaction network, including the DEGs and the royal- and worker-jelly enriched miRNA populations, revealed a subset of miRNAs potentially involved in regulating the expression of DEGs. The interaction of miR-34, miR-210, and miR-317 with Takeout, Neurotrophin-1, Forked, and Masquerade genes was experimentally confirmed using a luciferase reporter system. Taken together, our results reconstruct the regulatory network that governs the development of the early brain diphenism in honey bees.

Comparative analysis of phenotypic plasticity sheds light on the evolution and molecular underpinnings of locust phase polyphenism

Locusts exhibit one of nature's most spectacular examples of complex phenotypic plasticity, in which changes in density cause solitary and cryptic individuals to transform into gregarious and conspicuous locusts forming large migrating swarms. We investigated how these coordinated alternative phenotypes might have evolved by studying the Central American locust and three closely related non-swarming grasshoppers in a comparative framework. By experimentally isolating and crowding during nymphal development, we induced density-dependent phenotypic plasticity and quantified the resulting behavioural, morphological, and molecular reaction norms. All four species exhibited clear plasticity, but the individual reaction norms varied among species and showed different magnitudes. Transcriptomic responses were species-specific, but density-responsive genes were functionally similar across species. There were modules of co-expressed genes that were highly correlated with plastic reaction norms, revealing a potential molecular basis of density-dependent phenotypic plasticity. These findings collectively highlight the importance of studying multiple reaction norms from a comparative perspective.

Transcript Assembly and Quantification by RNA-Seq Reveals Significant Differences in Gene Expression and Genetic Variants in Mosquitoes of the Culex pipiens (Diptera: Culicidae) Complex

The taxonomy of Culex pipiens complex of mosquitoes is still debated, but in North America it is generally regarded to include Culex pipiens pipiens, Culex pipiens molestus, and Culex quinquefasciatus (or Culex pipiens quinquefasciatus). Although these mosquitoes have very similar morphometry, they each have unique life strategies specifically adapted to their ecological niche. Differences include the capability for overwintering diapause, bloodmeal preference, mating behaviors, and reliance on blood meals to produce eggs. Here, we used RNA-seq transcriptome analysis to investigate the differential gene expression and nucleotide polymorphisms that may link to the divergent traits specifically between Cx. pipiens pipiens and Cx. pipiens molestus.

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.

Characterization of Dendrolimus houi Lajonquiere (Lepidoptera: Lasiocampidae) Transcriptome across All Life Stages

Dendrolimus houi Lajonquiere is a phytophagous caterpillar infesting many economically important coniferous tree species in China, causing serious economic and ecological environment losses. Based on previous research, it has one generation per year in South China and East China in contrast to two generations per year in Yunnan province in southwestern China. The species is potentially resilient to climatic extremes in these regions with the eggs and 1st instar larvae surviving in the winter (5 °C), older instar larvae and pupae surviving high temperatures in the summer (35 °C), suggesting some temperature stress tolerance during different developmental stages. However, little is known in this species at the genetic and genomic level. In this study, we used high throughput sequencing to obtain transcriptome data from different developmental stages (eggs, 1st-3rd instar larvae, 4th-5th instar larvae, 6th-7th instar larvae, pupae, male and female adults), which were collected from Fujian province. In total, we obtained approximately 90 Gb of data, from which 33,720 unigenes were assembled and 17,797 unigenes were annotated. We furtherly analyzed the differentially expressed genes (DGEs) across all stages, the largest number between the eggs and 1st instar larvae stage and gene expression varied significantly in different developmental stages. Furthermore, 4138 SSR genes and 114,977 SNP loci were screened from transcriptome data. This paper will be a foundation for further study towards improved integrated pest management strategies for this species.

Transcriptomic and Functional Analyses of Phenotypic Plasticity in a Higher Termite, Macrotermes barneyi Light

Eusocial termites have a complex caste system, which leads to the division of labor. Previous studies offered some insight into the caste differentiation in lower termites; however, few studies were focusing on the molecular mechanisms of higher termites with sophisticated societies. Comparative transcriptomic analyses of five immature castes of a higher termite, Macrotermes barneyi Light, suggest that phenotypic plasticity is modulated by an array of transcriptional changes, including differentially expressed genes (e.g., caste-biased genes Vtg and TnC), co-expression networks (e.g., genes associated with nymph reproduction), and alternative splicing (e.g., events related to muscle development in presoldiers). Transcriptional (RT-PCR and RT-qPCR) and functional (in vivo RNAi) validation studies reveal multiple molecular mechanisms contributing to the phenotypic plasticity in eusocial termites. Molecular mechanisms governing the phenotypic plasticity in M. barneyi could be a rule rather than an exception in the evolution of sociality.

Developmental expression and evolution of hexamerin and haemocyanin from Folsomia candida (Collembola)

Haemocyanins constitute a group of copper-containing respiratory proteins, and hexamerins were derived from hexapod haemocyanin but lost the ability to transport oxygen and serve as storage proteins. Although hexamerins have been reported in most insect species, none of them has been identified in Collembola, one of the most primitive hexapod lineages, thereby preventing us from exploring relevant evolutionary scenarios regarding the origin and evolution of hexamerins in hexapods. Here we report on collembolan hexamerins for the first time, and investigated the temporal expression profiles of hexamerin and haemocyanin in the collembolan Folsomia candida. Haemocyanin was expressed over the entire life cycle, with higher expression at the embryonic stage than at other stages, whereas hexamerin expression was restricted to embryos, unlike insect hexamerins, which are generally expressed from larval to adult stages. A phylogenetic analysis and molecular clock estimation suggested that all investigated hexapod hexamerins have a single and ancient origin (~423 Ma), coincident with the rise of atmospheric oxygen levels in the Silurian-Devonian period, indicating a physiological link between molecular evolution and Palaeozoic oxygen changes.

Transcriptome changes reveal the genetic mechanisms of the reproductive plasticity of workers in lower termites

Background

The reproductive plasticity of termite workers provides colonies with tremendous flexibility to respond to environmental changes, which is the basis for evolutionary and ecological success. Although it is known that all colony members share the same genetic background and that differences in castes are caused by differences in gene expression, the pattern of the specific expression of genes involved in the differentiation of workers into reproductives remains unclear. In this study, the isolated workers of Reticulitermes labralis developed into reproductives, and then comparative transcriptomes were used for the first time to reveal the molecular mechanisms underlying the reproductive plasticity of workers.

Results

We identified 38,070 differentially expressed genes and found a pattern of gene expression involved in the differentiation of the workers into reproductives. 12, 543 genes were specifically upregulated in the isolated workers. Twenty-five signal transduction pathways classified into environmental information processing were related to the differentiation of workers into reproductives. Ras functions as a signalling switch regulates the reproductive plasticity of workers. The catalase gene which is related to longevity was up-regulated in reproductives.

Conclusion

We demonstrate that workers leaving the natal colony can induce the expression of stage-specific genes in the workers, which leads to the differentiation of workers into reproductives and suggests that the signal transduction along the Ras-MAPK pathway crucially controls the reproductive plasticity of the workers. This study also provides an important model for revealing the molecular mechanism of longevity changes.

Electronic supplementary material

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

Caste-based differential transcriptional expression of hexamerins in response to a juvenile hormone analog in the red imported fire ant (Solenopsis invicta)

The reproductive ground plan hypothesis proposes that gene networks regulating foraging behavior and reproductive female physiology in social insects emerged from ancestral gene and endocrine factor networks. Expression of storage proteins such as vitellogenins and hexamerins is an example of this co-option. Hexamerins, through their role modulating juvenile hormone availability, are involved in caste determination in termites. The genome of the fire ant (Solenopsis invicta) encodes four hexamerin genes, hexamerin-like (LOC105192919, hereafter called hexamerin 1), hexamerin (LOC105204474, hereafter called hexamerin 2), arylphorin subunit alpha-like, and arylphorin subunit beta. In this study, a phylogenetic analysis of the S. invicta hexamerins determined that each predicted protein clustered with one of the orthologous Apis mellifera hexamerins. Gene expression analyses by RT-qPCR revealed differential expression of the hexamerins between queens and workers, and between specific task-allocated workers (nurses and foragers). Queens and nurses had significantly higher expression of all genes when compared to foragers. Hexamerin 1 was expressed at higher levels in queens, while hexamerin 2 and arylphorin subunit beta were expressed at significantly higher levels in nurses. Arylphorin subunit alpha-like showed no significant difference in expression between virgin queens and nurses. Additionally, we analyzed the relationship between the expression of hexamerin genes and S-hydroprene, a juvenile hormone analog. Significant changes in hexamerin expression were recorded in nurses, virgin queens, and foragers 12 h after application of the analog. Hexamerin 1 and arylphorin subunit alpha-like expression were significantly lower after analog application in virgin queens. In foragers, hexamerin 2 and arylphorin subunit beta were significantly lower after analog application, while in nurses expression of all genes were significantly lower after analog application. Our results suggest that in S. invicta hexamerin genes could be associated with reproductive division of labor and task-allocation of workers.

Juvenile hormone and sesquiterpenoids in arthropods: Biosynthesis, signaling, and role of MicroRNA

Arthropod molting and reproduction are precisely controlled by the levels of sesquiterpenoids, a class of C15 hormones derived from three isoprene units. The two major functional arthropod sesquiterpenoids are juvenile hormone (JH) and methyl farnesoate (MF). In hemimetabolous insects (such as the aphids, bugs, and cockroaches) and holometabolous insects (such as beetles, bees, butterflies, and flies), dramatic decrease in the titers of JH and/or MF promote metamorphosis from larvae to adults either directly or through an intermediate pupal stage, respectively. JH is absent in crustaceans (lobster, shrimp, crab) and other arthropods (chelicerates such as ticks, mites, spiders, scorpions and myriapods such as millipede and centipedes). In some crustaceans, molting and reproduction is dependent on changing levels of MF. The regulation of sesquiterpenoid production is thus crucial in the life cycle of arthropods. Dynamic and complex mechanisms have evolved to regulate sesquiterpenoid production. Noncoding RNAs such as the microRNAs are primary regulators. This article provides an overview of microRNAs that are known to regulate sesquiterpenoid production in arthropods.

Chronic Nosema ceranae infection inflicts comprehensive and persistent immunosuppression and accelerated lipid loss in host Apis mellifera honey bees

Nosema ceranae is an intracellular microsporidian parasite of the Asian honey bee Apis cerana and the European honey bee Apis mellifera. Until relatively recently, A. mellifera honey bees were naïve to N. ceranae infection. Symptoms of nosemosis, or Nosema disease, in the infected hosts include immunosuppression, damage to gut epithelium, nutrient and energetic stress, precocious foraging and reduced longevity of infected bees. Links remain unclear between immunosuppression, the symptoms of nutrient and energetic stress, and precocious foraging behavior of hosts. To clarify physiological connections, we inoculated newly emerged A. mellifera adult workers with N. ceranae spores, and over 21 days post inoculation (21 days pi), gauged infection intensity and quantified expression of genes representing two innate immune pathways, Toll and Imd. Additionally, we measured each host's whole-body protein, lipids, carbohydrates and quantified respirometric and activity levels. Results show sustained suppression of genes of both humorally regulated immune response pathways after 6 days pi. At 7 days pi, elevated protein levels of infected bees may reflect synthesis of antimicrobial peptides from an initial immune response, but the lack of protein gain compared with uninfected bees at 14 days pi may represent low de novo protein synthesis. Carbohydrate data do not indicate that hosts experience severe metabolic stress related to this nutrient. At 14 days pi infected honey bees show high respirometric and activity levels, and corresponding lipid loss, suggesting lipids may be used as fuel for increased metabolic demands resulting from infection. Accelerated lipid loss during nurse honey bee behavioral development can have cascading effects on downstream physiology that may lead to precocious foraging, which is a major factor driving colony collapse.

Multiple functions of CREB-binding protein during postembryonic development: identification of target genes

Background

Juvenile hormones (JH) and ecdysteroids control postembryonic development in insects. They serve as valuable targets for pest management. Hence, understanding the molecular mechanisms of their action is of crucial importance. CREB-binding protein (CBP) is a universal transcriptional co-regulator. It controls the expression of several genes including those from hormone signaling pathways through co-activation of many transcription factors. However, the role of CBP during postembryonic development in insects is not well understood. Therefore, we have studied the role of CBP in postembryonic development in Tribolium, a model coleopteran insect.

Results

CBP is ubiquitously expressed in the red flour beetle, Tribolium castaneum. RNA interference (RNAi) mediated knockdown of CBP resulted in a decrease in JH induction of Kr-h1 gene expression in Tribolium larvae and led to a block in their development. Moreover, the injection of CBP double-stranded RNA (dsRNA) showed lethal phenotypes within 8 days of injection. RNA-seq and subsequent differential gene expression analysis identified CBP target genes in Tribolium. Knockdown of CBP caused a decrease in the expression of 1306 genes coding for transcription factors and other proteins associated with growth and development. Depletion of CBP impaired the expression of several JH response genes (e.g., Kr-h1, Hairy, early trypsin) and ecdysone response genes (EcR, E74, E75, and broad complex). Further, GO enrichment analyses of the downregulated genes showed enrichment in different functions including developmental processes, pigmentation, anatomical structure development, regulation of biological and cellular processes, etc.

Conclusion

These data suggest diverse but crucial roles for CBP during postembryonic development in the coleopteran model insect, Tribolium. It can serve as a target for RNAi mediated pest management of this stored product pest.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4373-3) contains supplementary material, which is available to authorized users.

Juvenile hormone and hemimetabolan eusociality: a comparison of cockroaches with termites

Termites are social Dictyoptera that evolved eusociality independently from social Hymenoptera. They are characterized by unique developmental plasticity that is the basis of caste differentiation and social organization. As developmental plasticity is a result of endocrine regulation, in order to understand the evolution of termite sociality it is helpful to compare the endocrine underpinning of development between termites and cockroaches. Nijhout and Wheeler (1982) proposed that varying JH titers determine caste differentiation in termites. Based on current results, we extend this model by adding the importance of social interactions. High JH titers in the presence of soldiers lead to regressive development (decrease in body size, apparent regression in development), while an absence of soldiers induces (pre-)soldier differentiation. On the opposite side, low JH titers in colonies headed by reproductives result in progressive molts toward adults, while an absence of reproductives induces development of replacement reproductives. In cockroaches, transcription factors involved in JH signaling, including the adult specifier E93 (the co-called MEKRE93 pathway) regulate the morphogenetic transition between the nymph and the adult. In termites, we speculate that castes might be determined by social effects playing a modulatory action of JH in the MEKRE93 pathway.

The Genetic Architecture of Ovariole Number in Drosophila melanogaster: Genes with Major, Quantitative, and Pleiotropic Effects

Ovariole number has a direct role in the number of eggs produced by an insect, suggesting that it is a key morphological fitness trait. Many studies have documented the variability of ovariole number and its relationship to other fitness and life-history traits in natural populations of Drosophila However, the genes contributing to this variability are largely unknown. Here, we conducted a genome-wide association study of ovariole number in a natural population of flies. Using mutations and RNAi-mediated knockdown, we confirmed the effects of 24 candidate genes on ovariole number, including a novel gene, anneboleyn (formerly CG32000), that impacts both ovariole morphology and numbers of offspring produced. We also identified pleiotropic genes between ovariole number traits and sleep and activity behavior. While few polymorphisms overlapped between sleep parameters and ovariole number, 39 candidate genes were nevertheless in common. We verified the effects of seven genes on both ovariole number and sleep: bin3, blot, CG42389, kirre, slim, VAChT, and zfh1 Linkage disequilibrium among the polymorphisms in these common genes was low, suggesting that these polymorphisms may evolve independently.

RNA interference unveils the importance of Pseudotrichonympha grassii cellobiohydrolase, a protozoan exoglucanase, in termite cellulose degradation

Based on prior work, a cellulase from glycosyl hydrolase family 7 (GHF7) was identified and found to be expressed at a high level in Coptotermes formosanus. To determine the function of GHF7 family members in vivo, we used RNA interference (RNAi) to functionally analyse the exoglucanase gene Pseudotrichonympha grassii cellobiohydrolase gene (PgCBH), which was highly expressed in Pseudotrichonympha grassii, a flagellate found in the hindgut of C. formosanus. In this study, the expression level of PgCBH was down-regulated by RNAi, causing the death of P. grassii, but no effect was observed for other flagellates found in C. formosanus. RNAi also resulted in significantly reduced exoglucanase activity, and no effect was observed for endoglucanase and β-glucosidase activities. This result demonstrated that the PgCBH gene plays a role in the protist lignocellulolytic process and is also important for host survival. PgCBH can be used as a target gene and has potential as a bioinsecticide for use against termites.

How ants send signals in saliva

Adult ants use saliva to transfer juvenile hormone and other chemical signals to their larvae.

Differential gene expression underlying ovarian phenotype determination in honey bee, Apis mellifera L., caste development

Adult honey bee queens and workers drastically differ in ovary size. This adult ovary phenotype difference becomes established during the final larval instar, when massive programmed cell death leads to the degeneration of 95-99% of the ovariole anlagen in workers. The higher juvenile hormone (JH) levels in queen larvae protect the ovaries against such degeneration. To gain insights into the molecular architecture underlying this divergence critical for adult caste fate and worker sterility, we performed a microarray analysis on fourth and early fifth instar queen and worker ovaries. For the fourth instar we found nine differentially expressed genes (DEGs) with log₂FC > 1.0, but this number increased to 56 in early fifth-instar ovaries. We selected 15 DEGs for quantitative PCR (RT-qPCR) analysis. Nine differed significantly by the variables caste and/or development. Interestingly, genes with enzyme functions were higher expressed in workers, while those related to transcription and signaling had higher transcript levels in queens. For the RT-qPCR confirmed genes we analyzed their response to JH. This revealed a significant up-regulation for two genes, a short chain dehydrogenase reductase (sdr) and a heat shock protein 90 (hsp90). Five other genes, including hsp60 and hexamerin 70b (hex70b), were significantly down-regulated by JH. The sdr gene had previously come up as differentially expressed in other transcriptome analyses on honey bee larvae and heat shock proteins are frequently involved in insect hormone responses, this making them interesting candidates for further functional assays.

Oral transfer of chemical cues, growth proteins and hormones in social insects

Social insects frequently engage in oral fluid exchange - trophallaxis - between adults, and between adults and larvae. Although trophallaxis is widely considered a food-sharing mechanism, we hypothesized that endogenous components of this fluid might underlie a novel means of chemical communication between colony members. Through protein and small-molecule mass spectrometry and RNA sequencing, we found that trophallactic fluid in the ant Camponotus floridanus contains a set of specific digestion- and non-digestion related proteins, as well as hydrocarbons, microRNAs, and a key developmental regulator, juvenile hormone. When C. floridanus workers' food was supplemented with this hormone, the larvae they reared via trophallaxis were twice as likely to complete metamorphosis and became larger workers. Comparison of trophallactic fluid proteins across social insect species revealed that many are regulators of growth, development and behavioral maturation. These results suggest that trophallaxis plays previously unsuspected roles in communication and enables communal control of colony phenotypes.

Transcriptional Analysis of The Adaptive Digestive System of The Migratory Locust in Response to Plant Defensive Protease Inhibitors

Herbivorous insects evolved adaptive mechanisms to compensate for the presence of plant defensive protease inhibitors (PI) in their food. The underlying regulatory mechanisms of these compensatory responses remain largely elusive. In the current study, we investigated the initiation of this adaptive response in the migratory locust, Locusta migratoria, via microarray analysis of gut tissues. Four hours after dietary uptake of PIs, 114 and 150 transcripts were respectively found up- or downregulated. The results suggest a quick trade-off between compensating for potential loss of digestive activity on the one hand, and stress tolerance, defense, and structural integrity of the gut on the other hand. We additionally addressed the role of a group of related upregulated hexamerin-like proteins in the PI-induced response. Simultaneous knockdown of corresponding transcripts by means of RNA interference resulted in a reduced capacity of the locust nymphs to cope with the effects of PI. Moreover, since insect hexamerins have been shown to bind Juvenile Hormone (JH), we also investigated the effect of JH on the proteolytic digestion in L. migratoria. Our results indicate that JH has a stimulatory effect on the expression of three homologous chymotrypsin genes, while knocking down the JH receptor (methoprene tolerant) led to opposite effects.

Physiological profiles associated with ceasing growth of unfertilized eggs produced by unmated queens in the subterranean termite Reticulitermes chinensis

In Reticulitermes chinensis, a close relative of R.speratus with asexual queen succession, unfertilized eggs can be produced but do not hatch as larvae. To explain this phenomenon, we analyzed the physiological differences between unfertilized eggs/unmated queens and fertilized eggs/mated queens. Fertilized eggs had significantly lower quantities of five amino acids (Cys, Met, Ile, Leu and Tyr), Ca, protein and cholesterol during development. The higher levels of four trace elements (Na, K, Zn and Fe) in fertilized eggs and their lower levels in mated queens indicated that mated queens might transfer these trace elements to fertilized eggs to aid development. The higher levels of Mn, triglycerides and serotonin in mated queens and higher levels of Mn and glucose in fertilized eggs suggested that these substances are very important for normal ovarian and embryonic growth. The different expression of three reproductive genes (vtg 1, rab 11 and JHE 1) suggested that they might be involved in the regulation of ovarian and embryonic growth. Overall, changes in these physiological indices may substantially affect ovarian and embryonic growth and inhibit development of unfertilized eggs in R. chinensis.

Summary: Changes in amino acids, trace elements, nutrient content, serotonin and reproductive genes may substantially affect ovarian and embryonic growth and inhibit development of unfertilized eggs in Reticulitermes chinensis.

Comparative Transcriptomics Reveals Key Gene Expression Differences between Diapausing and Non-Diapausing Adults of Culex pipiens

Diapause is a critical eco-physiological adaptation for winter survival in the West Nile Virus vector, Culex pipiens, but little is known about the molecular mechanisms that distinguish diapause from non-diapause in this important mosquito species. We used Illumina RNA-seq to simultaneously identify and quantify relative transcript levels in diapausing and non-diapausing adult females. Among 65,623,095 read pairs, we identified 41 genes with significantly different transcript abundances between these two groups. Transcriptome divergences between these two phenotypes include genes related to juvenile hormone synthesis, anaerobic metabolism, innate immunity and cold tolerance.

Methyl farnesoate synthesis is necessary for the environmental sex determination in the water flea Daphnia pulex

Sex-determination systems can be divided into two groups: genotypic sex determination (GSD) and environmental sex determination (ESD). ESD is an adaptive life-history strategy that allows control of sex in response to environmental cues in order to optimize fitness. However, the molecular basis of ESD remains largely unknown. The micro crustacean Daphnia pulex exhibits ESD in response to various external stimuli. Although methyl farnesoate (MF: putative juvenile hormone, JH, in daphnids) has been reported to induce male production in daphnids, the role of MF as a sex-determining factor remains elusive due to the lack of a suitable model system for its study. Here, we establish such a system for ESD studies in D. pulex. The WTN6 strain switches from producing females to producing males in response to the shortened day condition, while the MFP strain only produces females, irrespective of day-length. To clarify whether MF has a novel physiological role as a sex-determining factor in D. pulex, we demonstrate that a MF/JH biosynthesis inhibitor suppressed male production in WTN6 strain reared under the male-inducible condition, shortened day-length. Moreover, we show that juvenile hormone acid O-methyltransferase (JHAMT), a critical enzyme of MF/JH biosynthesis, displays MF-generating activity by catalyzing farnesoic acid. Expression of the JHAMT gene increased significantly just before the MF-sensitive period for male production in the WTN6 strain, but not in the MFP strain, when maintained under male-inducible conditions. These results suggest that MF synthesis regulated by JHAMT is necessary for male offspring production in D. pulex. Our findings provide novel insights into the genetic underpinnings of ESD and they begin to shed light on the physiological function of MF as a male-fate determiner in D. pulex.

Proteomic analysis reveals that COP9 signalosome complex subunit 7A (CSN7A) is essential for the phase transition of migratory locust

The migratory locust displays a reversible, density-dependent transition between the two phases of gregaria and solitaria. This phenomenon is a typical kind of behavior plasticity. Here, we report that COP9 signalosome complex subunit 7A (CSN7A) is involved in the regulation of locust phase transition. Firstly, 90 proteins were identified to express differentially between the two phases by quantitative proteomic analysis. Gregaria revealed higher levels in proteins related to structure formation, melanism and energy metabolism, whereas solitaria had more abundant proteins related to digestion, absorption and chemical sensing. Subsequently, ten proteins including CSN7A were found to reveal differential mRNA expression profiles between the two phases. The CSN7A had higher mRNA level in the gregaria as compared with the solitaria, and the mRNA amount in the gregaria decreased remarkably during the 32 h-isolation. However, the mRNA level in the solitaria kept constant during the crowding rearing. Finally and importantly, RNA interference of CSN7A in gregaria resulted in obvious phase transition towards solitaria within 24 h. It suggests that CSN7A plays an essential role in the transition of gregaria towards solitaria in the migratory locust. To our knowledge, it's the first time to report the role of CSN in behavior plasticity of animals.

Expression and evolution of hexamerins from the tobacco hornworm, Manduca sexta, and other Lepidoptera

Hexamerins are large hemolymph-proteins that accumulate during the late larval stages of insects. Hexamerins have emerged from hemocyanin, but have lost the ability to bind oxygen. Hexamerins are mainly considered as storage proteins for non-feeding stages, but may also have other functions, e.g. in cuticle formation, transport and immune response. The genome of the hornworm Manduca sexta harbors six hexamerin genes. Two of them code for arylphorins (Msex2.01690, Msex2.15504) and two genes correspond to a methionine-rich hexamerin (Msex2.10735) and a moderately methionine-rich hexamerin (Msex2.01694), respectively. Two other genes do not correspond to any known hexamerin and distantly resemble the arylphorins (Msex2.01691, Msex2.01693). Five of the six hexamerin genes are clustered within ∼45 kb on scaffold 00023, which shows conserved synteny in various lepidopteran genomes. The methionine-rich hexamerin gene is located at a distinct site. M. sexta and other Lepidoptera have lost the riboflavin-binding hexamerin. With the exception of Msex2.01691, which displays low mRNA levels throughout the life cycle, all hexamerins are most highly expressed during pre-wandering phase of the 5th larval instar of M. sexta, supporting their role as storage proteins. Notably, Msex2.01691 is most highly expressed in the brain, suggesting a divergent function. Phylogenetic analyses showed that hexamerin evolution basically follows insect systematics. Lepidoptera display an unparalleled diversity of hexamerins, which exceeds that of other hexapod orders. In contrast to previous analyses, the lepidopteran hexamerins were found monophyletic. Five distinct types of hexamerins have been identified in this order, which differ in terms of amino acid composition and evolutionary history: i. the arylphorins, which are rich in aromatic amino acids (∼20% phenylalanine and tyrosine), ii. the distantly related arylphorin-like hexamerins, iii. the methionine-rich hexamerins, iv. the moderately methionine rich hexamerins, and v. the riboflavin-binding hexamerins.

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.

Termites as targets and models for biotechnology

Termites have many unique evolutionary adaptations associated with their eusocial lifestyles. Recent omics research has created a wealth of new information in numerous areas of termite biology (e.g., caste polyphenism, lignocellulose digestion, and microbial symbiosis) with wide-ranging applications in diverse biotechnological niches. Termite biotechnology falls into two categories: (a) termite-targeted biotechnology for pest management purposes, and (b) termite-modeled biotechnology for use in various industrial applications. The first category includes several candidate termiticidal modes of action such as RNA interference, digestive inhibition, pathogen enhancement, antimicrobials, endocrine disruption, and primer pheromone mimicry. In the second category, termite digestomes are deep resources for host and symbiont lignocellulases and other enzymes with applications in a variety of biomass, industrial, and processing applications. Moving forward, one of the most important approaches for accelerating advances in both termite-targeted and termite-modeled biotechnology will be to consider host and symbiont together as a single functional unit.

Phenotypic plasticity of HSP70s gene expression during diapause: signs of evolutionary responses to cold stress among Soybean Pod Borer populations (Leguminivora glycinivorella) in Northeast of China

The soybean pod borer (Leguminivora glycinivorella Matsumura) successfully survives the winter because of its high expression of 70-kDa heat shock proteins (HSP70s) during its overwintering diapause. The amount of HSP70s is different under different environmental stresses. In this study, inducible heat shock protein 70 and its constitutive heat shock cognate 70 were cloned by RT-PCR and RACE. These genes were named Lg-hsp70 and Lg-hsc70, respectively. Gene transcription and protein expression after cold stress treatment (5°C to −5°C) were analyzed by western blotting and by qRT-PCR for four populations that were sampled in the northeast region of China, including Shenyang, Gongzhuling, Harbin and Heihe, when the soybean pod borer was in diapause. As the cold shock temperature decreased, the levels of Lg-HSP70s were significantly up-regulated. The amount of cold-induced Lg-HSP70s was highest in the southernmost population (Shenyang, 41°50′N) and lowest in the northernmost population (Heihe, 50°22′N). These results support the hypothesis that the soybean pod borer in the northeast region of China displays phenotypic plasticity, and the accumulation of Lg-HSP70s is a strategy for overcoming environmental stress. These results also suggest that the induction of HSP70 synthesis, which is a complex physiological adaptation, can evolve quickly and inherit stability.

Endocrine regulation of predator-induced phenotypic plasticity

Elucidating the developmental and genetic control of phenotypic plasticity remains a central agenda in evolutionary ecology. Here, we investigate the physiological regulation of phenotypic plasticity induced by another organism, specifically predator-induced phenotypic plasticity in the model ecological and evolutionary organism Daphnia pulex. Our research centres on using molecular tools to test among alternative mechanisms of developmental control tied to hormone titres, receptors and their timing in the life cycle. First, we synthesize detail about predator-induced defenses and the physiological regulation of arthropod somatic growth and morphology, leading to a clear prediction that morphological defences are regulated by juvenile hormone and life-history plasticity by ecdysone and juvenile hormone. We then show how a small network of genes can differentiate phenotype expression between the two primary developmental control pathways in arthropods: juvenoid and ecdysteroid hormone signalling. Then, by applying an experimental gradient of predation risk, we show dose-dependent gene expression linking predator-induced plasticity to the juvenoid hormone pathway. Our data support three conclusions: (1) the juvenoid signalling pathway regulates predator-induced phenotypic plasticity; (2) the hormone titre (ligand), rather than receptor, regulates predator-induced developmental plasticity; (3) evolution has favoured the harnessing of a major, highly conserved endocrine pathway in arthropod development to regulate the response to cues about changing environments (risk) from another organism (predator).

Olive fly transcriptomics analysis implicates energy metabolism genes in spinosad resistance

Background

The olive fly, Bactrocera oleae, is the most devastating pest of cultivated olives. Its control has been traditionally based on insecticides, mainly organophosphates and pyrethroids. In recent years, the naturalyte spinosad is used against the olive fly. As with other insecticides, spinosad is subject to selection pressures that have led to resistance development. Mutations in the α6 subunit of the nicotinic acetylcholine receptor (nAChR) have been implicated in spinosad resistance in several species (e.g., Drosophila melanogaster) but excluded in others (e.g., Musca domestica). Yet, additional mechanisms involving enhanced metabolism of detoxification enzymes (such as P450 monooxygenases or mixed function oxidases) have also been reported. In order to clarify the spinosad resistance mechanisms in the olive fly, we searched for mutations in the α6-subunit of the nAChR and for up-regulated genes in the entire transcriptome of spinosad resistant olive flies.

Results

The olive fly α6-subunit of the nAChR was cloned from the laboratory sensitive strain and a spinosad selected resistant line. The differences reflected silent nucleotide substitutions or conserved amino acid changes. Additionally, whole transcriptome analysis was performed in the two strains in order to reveal any underlying resistance mechanisms. Comparison of over 13,000 genes showed that in spinosad resistant flies nine genes were significantly over-expressed, whereas ~40 were under-expressed. Further functional analyses of the nine over-expressed and eleven under-expressed loci were performed. Four of these loci (Yolk protein 2, ATP Synthase FO subunit 6, Low affinity cationic amino acid transporter 2 and Serine protease 6) showed consistently higher expression both in the spinosad resistant strain and in wild flies from a resistant California population. On the other side, two storage protein genes (HexL1 and Lsp1) and two heat-shock protein genes (Hsp70 and Hsp23) were unfailingly under-expressed in resistant flies.

Conclusion

The observed nucleotide differences in the nAChR-α6 subunit between the sensitive and spinosad resistant olive fly strains did not advocate for the involvement of receptor mutations in spinosad resistance. Instead, the transcriptome comparison between the two strains indicated that several immune system loci as well as elevated energy requirements of the resistant flies might be necessary to lever the detoxification process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-714) contains supplementary material, which is available to authorized users.

Transcriptome analysis reveals signature of adaptation to landscape fragmentation

We characterize allelic and gene expression variation between populations of the Glanville fritillary butterfly (Melitaea cinxia) from two fragmented and two continuous landscapes in northern Europe. The populations exhibit significant differences in their life history traits, e.g. butterflies from fragmented landscapes have higher flight metabolic rate and dispersal rate in the field, and higher larval growth rate, than butterflies from continuous landscapes. In fragmented landscapes, local populations are small and have a high risk of local extinction, and hence the long-term persistence at the landscape level is based on frequent re-colonization of vacant habitat patches, which is predicted to select for increased dispersal rate. Using RNA-seq data and a common garden experiment, we found that a large number of genes (1,841) were differentially expressed between the landscape types. Hexamerin genes, the expression of which has previously been shown to have high heritability and which correlate strongly with larval development time in the Glanville fritillary, had higher expression in fragmented than continuous landscapes. Genes that were more highly expressed in butterflies from newly-established than old local populations within a fragmented landscape were also more highly expressed, at the landscape level, in fragmented than continuous landscapes. This result suggests that recurrent extinctions and re-colonizations in fragmented landscapes select a for specific expression profile. Genes that were significantly up-regulated following an experimental flight treatment had higher basal expression in fragmented landscapes, indicating that these butterflies are genetically primed for frequent flight. Active flight causes oxidative stress, but butterflies from fragmented landscapes were more tolerant of hypoxia. We conclude that differences in gene expression between the landscape types reflect genomic adaptations to landscape fragmentation.

Evolutionary implications of dipluran hexamerins

Hexamerin, as a member of the highly conserved arthropod hemocyanin superfamily, has been shown to be a good marker for the phylogenetic study of insects. However, few studies have been conducted on hexamerins in basal hexapods. The first Diplura hexamerin CspHex1 was reported only recently (Pick and Burmester, 2009). Remarkably, CspHex1 was suggested to have evolved from hexapod hemocyanin subunit type 2, which is very different from all insect hexamerins originated from hexapod hemocyanin subunit type 1. Does this finding suggest double or even multiple origins of hexamerins in Hexapoda? To find more evidence on the evolution of dipluran hexamerins, eight putative hexamerin gene sequences were obtained from three dipluran species, as were three hemocyanin genes from two collembolan species. Unexpectedly, after adding the new sequences into the phylogenetic analyses, all dipluran hexamerins including CspHex1 grouped together and as sister to the insect hexamerins, with high likelihood and Bayesian support. Our analysis supports a single origin of the hexamerins in Hexapoda, and suggests the close relationship between Diplura and Insecta. In addition, our study indicates that a relatively comprehensive taxa sampling is essential to solve some problems in phylogenetic reconstruction.

Endocrine control of cuticular hydrocarbon profiles during worker-to-soldier differentiation in the termite Reticulitermes flavipes

The social organization of termites, unlike that of other social insects, is characterized by a highly plastic caste system. With the exception of the alates, all other individuals in a colony remain at an immature stage of development. Workers in particular remain developmentally flexible; they can switch castes to become soldiers or neotenics. Juvenile hormone (JH) is known to play a key role in turning workers into soldiers. In this study, we analyzed differences in cuticular hydrocarbon (CHC) profiles among castes, paying particular attention to the transition of workers to soldiers, in the subterranean termite species Reticulitermes flavipes. CHCs have a fundamental function in social insects as they serve as cues in inter- and intraspecific recognition. We showed that (1) the CHC profiles of the different castes (workers, soldiers, nymphs and neotenics) are different and (2) when workers were experimentally exposed to a JH analog and thus induced to become soldiers, their CHC profiles were modified before and after the worker-presoldier molt and before and after the presoldier-soldier molt.

Integrating phenotypic plasticity within an Ecological Genomics framework: recent insights from the genomics, evolution, ecology, and fitness of plasticity

E.B. Ford's 1964 book Ecological Genetics was a call for biologists to engage in multidisciplinary work in order to elucidate the link between genotype, phenotype, and fitness for ecologically relevant traits. In this review, we argue that the integration of an ecological genomics framework in studies of phenotypic plasticity is a promising approach to elucidate the causal links between genes and the environment, particularly during colonization of novel environments, environmental change, and speciation. This review highlights some of the questions and hypotheses generated from a mechanistic, evolutionary, and ecological perspective, in order to direct the continued and future use of genomic tools in the study of phenotypic plasticity.

Biased gene expression in early honeybee larval development

Background

Female larvae of the honeybee (Apis mellifera) develop into either queens or workers depending on nutrition. This nutritional stimulus triggers different developmental trajectories, resulting in adults that differ from each other in physiology, behaviour and life span.

Results

To understand how these trajectories are established we have generated a comprehensive atlas of gene expression throughout larval development. We found substantial differences in gene expression between worker and queen-destined larvae at 6 hours after hatching. Some of these early changes in gene expression are maintained throughout larval development, indicating that caste-specific developmental trajectories are established much earlier than previously thought. Within our gene expression data we identified processes that potentially underlie caste differentiation. Queen-destined larvae have higher expression of genes involved in transcription, translation and protein folding early in development with a later switch to genes involved in energy generation. Using RNA interference, we were able to demonstrate that one of these genes, hexamerin 70b, has a role in caste differentiation. Both queen and worker developmental trajectories are associated with the expression of genes that have alternative splice variants, although only a single variant of a gene tends to be differentially expressed in a given caste.

Conclusions

Our data, based on the biases in gene expression early in development together with published data, supports the idea that caste development in the honeybee consists of two phases; an initial biased phase of development, where larvae can still switch to the other caste by differential feeding, followed by commitment to a particular developmental trajectory.