Francis Morais Franco Nunes

Genomic signatures of evolutionary transitions from solitary to group living

For bees, many roads lead to social harmony Eusociality, where workers sacrifice their reproductive rights to support the colony, has evolved repeatedly and represents the most evolved form of social evolution in insects. Kapheim et al. looked across the genomes of 10 bee species with varying degrees of sociality to determine the underlying genomic contributions. No one genomic path led to eusociality, but similarities across genomes were seen in features such as increases in gene regulation and methylation. It also seems that selection pressures relaxed after the emergence of complex sociality. Science, this issue p. 1139 Social evolution in bees has followed diverse genomic paths but shares genomic patterns. The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Caste development and reproduction: a genome-wide analysis of hallmarks of insect eusociality

The honey bee queen and worker castes are a model system for developmental plasticity. We used established expressed sequence tag information for a Gene Ontology based annotation of genes that are differentially expressed during caste development. Metabolic regulation emerged as a major theme, with a caste‐specific difference in the expression of oxidoreductases vs. hydrolases. Motif searches in upstream regions revealed group‐specific motifs, providing an entry point to cis‐regulatory network studies on caste genes. For genes putatively involved in reproduction, meiosis‐associated factors came out as highly conserved, whereas some determinants of embryonic axes either do not have clear orthologs (bag of marbles, gurken, torso), or appear to be lacking (trunk) in the bee genome. Our results are the outcome of a first genome‐based initiative to provide an annotated framework for trends in gene regulation during female caste differentiation (representing developmental plasticity) and reproduction.

The gene vitellogenin affects microRNA regulation in honey bee (Apis mellifera) fat body and brain

SUMMARY In honey bees, vitellogenin (Vg) is hypothesized to be a major factor affecting hormone signaling, food-related behavior, immunity, stress resistance and lifespan. MicroRNAs, which play important roles in post-transcriptional gene regulation, likewise affect many biological processes. The actions of microRNAs and Vg are known to intersect in the context of reproduction; however, the role of these associations on social behavior is unknown. The phenotypic effects of Vg knockdown are best established and studied in the forager stage of workers. Thus, we exploited the well-established RNA interference (RNAi) protocol for Vg knockdown to investigate its downstream effects on microRNA population in honey bee foragers brain and fat body tissue. To identify microRNAs that are differentially expressed between tissues in control and knockdown foragers, we used μParaflo microfluidic oligonucleotide microRNA microarrays. Our results showed that 76 and 74 microRNAs were expressed in the brain of control and knockdown foragers whereas 66 and 69 microRNAs were expressed in the fat body of control and knockdown foragers, respectively. Target prediction identified potential seed matches for a differentially expressed subset of microRNAs affected by Vg knockdown. These candidate genes are involved in a broad range of biological processes including insulin signaling, juvenile hormone (JH) and ecdysteroid signaling previously shown to affect foraging behavior. Thus, here we demonstrate a causal link between the Vg knockdown forager phenotype and variation in the abundance of microRNAs in different tissues, with possible consequences for the regulation of foraging behavior.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera).

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A+ transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

Exoskeleton formation in Apis mellifera: Cuticular hydrocarbons profiles and expression of desaturase and elongase genes during pupal and adult development

Cuticular hydrocarbons (CHCs) are abundant in the superficial cuticular layer (envelope) of insects where they play roles as structural, anti-desiccation and semiochemical compounds. Many studies have investigated the CHC composition in the adult insects. However, studies on the profiles of these compounds during cuticle formation and differentiation are scarce and restrict to specific stages of a few insect species. We characterized the CHCs developmental profiles in the honeybee workers during an entire molting cycle (from pupal-to-adult ecdyses) and in mature adults (forager bees). Gas chromatography/mass spectrometry (GC/MS) analysis revealed remarkable differences in the relative quantities of CHCs, thus discriminating pupae, developing and newly-ecdysed adults, and foragers from each other. In parallel, the honeybee genome database was searched for predicted gene models using known amino acid sequences of insect enzymes catalyzing lipid desaturation (desaturases) or elongation (elongases) as queries in BLASTP analysis. The expression levels of six desaturase genes and ten elongase genes potentially involved in CHC biosynthesis were determined by reverse transcription and real time polymerase chain reaction (RT-qPCR) in the developing integument (cuticle and subjacent epidermis). Aiming to predict roles for these genes in CHC biosynthesis, the developmental profiles of CHCs and desaturase/elongase transcript levels were evaluated using Spearman correlation coefficient. This analysis pointed to differential roles for these gene products in the biosynthesis of certain CHC classes. Based on the assumption that homologous proteins may share a similar function, phylogenetic trees were reconstructed as an additional strategy to predict functions and evolutionary relationships of the honeybee desaturases and elongases. Together, these approaches highlighted the molecular complexity underlying the formation of the lesser known layer of the cuticular exoskeleton, the envelope.

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.).

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep‐sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR‐1, miR‐31a, miR‐13b, miR‐125, let‐7 RNA, miR‐100, miR‐276, miR‐12, miR‐263a, miR‐306, miR‐317, miR‐92a and miR‐9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera - A positive loop linking JH, ftz-f1, and vitellogenin

Ftz‐f1 is an orphan member of the nuclear hormone receptor superfamily. A 20‐hydroxyecdysone pulse allows ftz‐f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone‐esterase during late pharate‐adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz‐f1 is known for mediating intracellular JH signalling, we hypothesized that ftz‐f1 could mediate JH action during the pharate‐adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz‐f1 has caste‐specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi‐mediated knock down of ftz‐f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz‐f1 expression. Finally, a double‐repressor hypothesis‐inspired vg gene knock‐down experiment suggests the existence of a positive molecular loop between JH, ftz‐f1 and vg.