Unique spatially and temporary-regulated/sex-specific expression of a long ncRNA, Nb-1, suggesting its pleiotropic functions associated with honey bee lifecycle.

Honey bees are social insects, and each colony member has unique morphological and physiological traits associated with their social tasks. Previously, we identified a long non-coding RNA from honey bees, termed Nb-1, whose expression in the brain decreases associated with the age-polyethism of workers and is detected in some neurosecretory cells and octopaminergic neurons, suggesting its role in the regulation of worker labor transition. Herein, we investigated its spatially and temporary-regulated/sex-specific expression. Nb-1 was expressed as an abundant maternal RNA during oogenesis and embryogenesis in both sexes. In addition, Nb-1 was expressed preferentially in the proliferating neuroblasts of the mushroom bodies (a higher-order center of the insect brain) in the pupal brains, suggesting its role in embryogenesis and mushroom body development. On the contrary, Nb-1 was expressed in a drone-specific manner in the pupal and adult retina, suggesting its role in the drone visual development and/or sense. Subcellular localization of Nb-1 in the brain during development differed depending on the cell type. Considering that Nb-1 is conserved only in Apidae, our findings suggest that Nb-1 potentially has pleiotropic functions in the expression of multiple developmental, behavioral, and physiological traits, which are closely associated with the honey bee lifecycle.

regeneration factors expressed on myeloid expression in macrophage-like cells is required for tail regeneration in Xenopus laevis tadpoles.

Xenopus laevis tadpoles can regenerate whole tails after amputation. We have previously reported that interleukin 11 (il11) is required for tail regeneration. In this study, we have screened for genes that support tail regeneration under Il11 signaling in a certain cell type and have identified the previously uncharacterized genes Xetrov90002578m.L and Xetrov90002579m.S [referred to hereafter as regeneration factors expressed on myeloid.L (rfem.L) and rfem.S]. Knockdown (KD) of rfem.L and rfem.S causes defects of tail regeneration, indicating that rfem.L and/or rfem.S are required for tail regeneration. Single-cell RNA sequencing (scRNA-seq) revealed that rfem.L and rfem.S are expressed in a subset of leukocytes with a macrophage-like gene expression profile. KD of colony-stimulating factor 1 (csf1), which is essential for macrophage differentiation and survival, reduced rfem.L and rfem.S expression levels and the number of rfem.L- and rfem.S-expressing cells in the regeneration bud. Furthermore, forced expression of rfem.L under control of the mpeg1 promoter, which drives rfem.L in macrophage-like cells, rescues rfem.L and rfem.S KD-induced tail regeneration defects. Our findings suggest that rfem.L or rfem.S expression in macrophage-like cells is required for tail regeneration.

Identification of ecdysone receptor target genes in the worker honey bee brains during foraging behavior.

Ecdysone signaling plays central roles in morphogenesis and female ovarian development in holometabolous insects. In the European honey bee (Apis mellifera L.), however, ecdysone receptor (EcR) is expressed in the brains of adult workers, which have already undergone metamorphosis and are sterile with shrunken ovaries, during foraging behavior. Aiming at unveiling the significance of EcR signaling in the worker brain, we performed chromatin-immunoprecipitation sequencing of EcR to search for its target genes using the brains of nurse bees and foragers. The majority of the EcR targets were common between the nurse bee and forager brains and some of them were known ecdysone signaling-related genes. RNA-sequencing analysis revealed that some EcR target genes were upregulated in forager brains during foraging behavior and some were implicated in the repression of metabolic processes. Single-cell RNA-sequencing analysis revealed that EcR and its target genes were expressed mostly in neurons and partly in glial cells in the optic lobes of the forager brain. These findings suggest that in addition to its role during development, EcR transcriptionally represses metabolic processes during foraging behavior in the adult worker honey bee brain.

Evolutionary dynamics of mushroom body Kenyon cell types in hymenopteran brains from multifunctional type to functionally specialized types.

Evolutionary dynamics of diversification of brain neuronal cell types that have underlain behavioral evolution remain largely unknown. Here, we compared transcriptomes and functions of Kenyon cell (KC) types that compose the mushroom bodies between the honey bee and sawfly, a primitive hymenopteran insect whose KCs likely have the ancestral properties. Transcriptome analyses show that the sawfly KC type shares some of the gene expression profile with each honey bee KC type, although unique gene expression profiles have also been acquired in each honey bee KC type. In addition, functional analysis of two sawfly genes suggested that the functions in learning and memory of the ancestral KC type were heterogeneously inherited among the KC types in the honey bee. Our findings strongly suggest that the functional evolution of KCs in Hymenoptera involved two previously hypothesized processes for evolution of cell function: functional segregation and divergence.

Mblk-1/E93, an ecdysone related-transcription factor, targets synaptic plasticity-related genes in the honey bee mushroom bodies.

Among hymenopteran insects, aculeate species such as bees, ants, and wasps have enlarged and morphologically elaborate mushroom bodies (MBs), a higher-order brain center in the insect, implying their relationship with the advanced behavioral traits of aculeate species. The molecular bases leading to the acquisition of complicated MB functions, however, remains unclear. We previously reported the constitutive and MB-preferential expression of an ecdysone-signaling related transcription factor, Mblk-1/E93, in the honey bee brain. Here, we searched for target genes of Mblk-1 in the worker honey bee MBs using chromatin immunoprecipitation sequence analyses and found that Mblk-1 targets several genes involved in synaptic plasticity, learning, and memory abilities. We also demonstrated that Mblk-1 expression is self-regulated via Mblk-1-binding sites, which are located upstream of Mblk-1. Furthermore, we showed that the number of the Mblk-1-binding motif located upstream of Mblk-1 homologs increased associated with evolution of hymenopteran insects. Our findings suggest that Mblk-1, which has been focused on as a developmental gene transiently induced by ecdysone, has acquired a novel expression pattern to play a role in synaptic plasticity in honey bee MBs, raising a possibility that molecular evolution of Mblk-1 may have partly contributed to the elaboration of MB function in insects.

Effective enrichment of stem cells in regenerating Xenopus laevis tadpole tails using the side population method.

Xenopus laevis tadpoles have a strong regenerative ability and can regenerate their whole tails after tail amputation. Lineage-restricted tissue stem cells are thought to provide sources for the regenerating tissues by producing undifferentiated progenitor cells in response to tail amputation. However, elucidating the behavioral dynamics of tissue stem cells during tail regeneration is difficult because of their rarity, and there are few established methods of isolating these cells in amphibians. Here, to detect and analyze rare tissue stem cells, we attempted to enrich tissue stem cells from tail regeneration buds. High Hoechst dye efflux capacity is thought to be a common characteristic of several types of mammalian tissue stem cells; these stem cells, designated as the "side population (SP)," may be enriched by flow cytometry (SP method). To evaluate the effectiveness of stem cell enrichment using the SP method in regenerating X. laevis tadpole tails, we performed single-cell RNA sequencing (scRNA-seq) of SP cells from regeneration buds and analyzed the frequency of satellite cells, which are muscle stem/progenitor cells expressing pax7. The pax7-expressing cells were enriched in the SP compared with whole normal tails and regeneration buds. Furthermore, hes1-expressing cells, which are assumed to be neural stem/progenitor cells, were also enriched in the SP. Our findings suggest that the SP method is efficient for successfully enriching tissue stem cells in regenerating X. laevis tadpole tails, indicating that the combination of the SP method and scRNA-seq is useful for studying tissue stem cells that contribute to tail regeneration.

Cellular responses in the FGF10-mediated improvement of hindlimb regenerative capacity in Xenopus laevis revealed by single-cell transcriptomics.

Xenopus laevis tadpoles possess regenerative capacity in their hindlimb buds at early developmental stages (stages ~52-54); they can regenerate complete hindlimbs with digits after limb bud amputation. However, they gradually lose their regenerative capacity as metamorphosis proceeds. Tadpoles in late developmental stages regenerate fewer digits (stage ~56), or only form cartilaginous spike without digits or joints (stage ~58 or later) after amputation. Previous studies have shown that administration of fibroblast growth factor 10 (FGF10) in late-stage (stage 56) tadpole hindlimb buds after amputation can improve their regenerative capacity, which means that the cells responding to FGF10 signaling play an important role in limb bud regeneration. In this study, we performed single-cell RNA sequencing (scRNA-seq) of hindlimb buds that were amputated and administered FGF10 by implanting FGF10-soaked beads at a late stage (stage 56), and explored cell clusters exhibiting a differential gene expression pattern compared with that in controls treated with phosphate-buffered saline. The scRNA-seq data showed expansion of fgf8-expressing cells in the cluster of the apical epidermal cap of FGF10-treated hindlimb buds, which was reported previously, indicating that the administration of FGF10 was successful. On analysis, in addition to the epidermal cluster, a subset of myeloid cells and a newly identified cluster of steap4-expressing cells showed remarkable differences in their gene expression profiles between the FGF10- or phosphate-buffered saline-treatment conditions, suggesting a possible role of these clusters in improving the regenerative capacity of hindlimbs via FGF10 administration.

Xenopus laevis il11ra.L is an experimentally proven interleukin-11 receptor component that is required for tadpole tail regeneration.

Xenopus laevis tadpoles possess high regenerative ability and can regenerate functional tails after amputation. An early event in regeneration is the induction of undifferentiated cells that form the regenerated tail. We previously reported that interleukin-11 (il11) is upregulated immediately after tail amputation to induce undifferentiated cells of different cell lineages, indicating a key role of il11 in initiating tail regeneration. As Il11 is a secretory factor, Il11 receptor-expressing cells are thought to mediate its function. X. laevis has a gene annotated as interleukin 11 receptor subunit alpha on chromosome 1L (il11ra.L), a putative subunit of the Il11 receptor complex, but its function has not been investigated. Here, we show that nuclear localization of phosphorylated Stat3 induced by Il11 is abolished in il11ra.L knocked-out culture cells, strongly suggesting that il11ra.L encodes an Il11 receptor component. Moreover, knockdown of il11ra.L impaired tadpole tail regeneration, suggesting its indispensable role in tail regeneration. We also provide a model showing that Il11 functions via il11ra.L-expressing cells in a non-cell autonomous manner. These results highlight the importance of il11ra.L-expressing cells in tail regeneration.

Low-temperature incubation improves both knock-in and knock-down efficiencies by the CRISPR/Cas9 system in Xenopus laevis as revealed by quantitative analysis.

The recent development of the CRISPR/Cas9-mediated gene editing technique has provided various gene knock-down and knock-in methods for Xenopus laevis. Gene-edited F0 individuals created by these methods, however, are mosaics with both mutated/knocked-in and unedited wild-type cells, and therefore precise determination and higher efficiency of knock-down and knock-in methods are desirable, especially for analyses of F0 individuals. To clarify the ratio of cells that are gene-edited by CRISPR/Cas9 methods to the whole cells in F0 individuals, we subjected Inference of CRISPR Edits analysis for knock-down experiments and flow cytometry for knock-in experiments to the F0 individuals. With these quantitative methods, we showed that low-temperature incubation of X. laevis embryos after microinjection improved the mutation rate in the individuals. Moreover, we applied low-temperature incubation when using a knock-in method with long single-strand DNA and found improved knock-in efficiency. Our results provide a simple and useful way to evaluate and improve the efficiency of gene editing in X. laevis.

Publisher Correction: Developmental stage-specific distribution and phosphorylation of Mblk-1, a transcription factor involved in ecdysteroid-signaling in the honey bee brain.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Neural activity mapping of bumble bee (Bombus ignitus) brains during foraging flight using immediate early genes.

Honey bees and bumble bees belong to the same family (Apidae) and their workers exhibit a division of labor, but the style of division of labor differs between species. The molecular and neural bases of the species-specific social behaviors of Apidae workers have not been analyzed. Here, we focused on two immediate early genes, hormone receptor 38 (HR38) and early growth response gene-1 (Egr1), and late-upregulated ecdysone receptor (EcR), all of which are upregulated by foraging flight and expressed preferentially in the small-type Kenyon cells of the mushroom bodies (MBs) in the honey bee brain. Gene expression analyses in Bombus ignitus revealed that HR38 and Egr1, but not EcR, exhibited an immediate early response during awakening from CO2 anesthesia. Both premature mRNA for HR38 and mature mRNA for Egr1 were induced during foraging flight, and mRNAs for HR38 and Egr1 were sparsely detected inside the whole MB calyces. In contrast, EcR expression was higher in forager brains than in nurse bees and was expressed preferentially in the small-type Kenyon cells inside the MBs. Our findings suggest that Kenyon cells are active during foraging flight and that the function of late-upregulated EcR in the brain is conserved among these Apidae species.

Developmental stage-specific distribution and phosphorylation of Mblk-1, a transcription factor involved in ecdysteroid-signaling in the honey bee brain.

In the honey bee, the mushroom bodies (MBs), a higher-order center in insect brain, comprise interneurons termed Kenyon cells (KCs). We previously reported that Mblk-1, which encodes a transcription factor involved in ecdysteroid-signaling, is expressed preferentially in the large-type KCs (lKCs) in the pupal and adult worker brain and that phosphorylation by the Ras/MAPK pathway enhances the transcriptional activity of Mblk-1 in vitro. In the present study, we performed immunoblotting and immunofluorescence studies using affinity-purified anti-Mblk-1 and anti-phosphorylated Mblk-1 antibodies to analyze the distribution and phosphorylation of Mblk-1 in the brains of pupal and adult workers. Mblk-1 was preferentially expressed in the lKCs in both pupal and adult worker brains. In contrast, some Mblk-1 was phosphorylated almost exclusively in the pupal stages, and phosphorylated Mblk-1 was preferentially expressed in the MB neuroblasts and lKCs in pupal brains. Immunofluorescence studies revealed that both Mblk-1 and phosphorylated Mblk-1 are located in both the cytoplasm and nuclei of the lKC somata in the pupal and adult worker brains. These findings suggest that Mblk-1 plays a role in the lKCs in both pupal and adult stages and that phosphorylated Mblk-1 has pupal stage-specific functions in the MB neuroblasts and lKCs in the honey bee brain.

Genetics in the Honey Bee: Achievements and Prospects toward the Functional Analysis of Molecular and Neural Mechanisms Underlying Social Behaviors.

The European honey bee is a model organism for studying social behaviors. Comprehensive analyses focusing on the differential expression profiles of genes between the brains of nurse bees and foragers, or in the mushroom bodies-the brain structure related to learning and memory, and multimodal sensory integration-has identified candidate genes related to honey bee behaviors. Despite accumulating knowledge on the expression profiles of genes related to honey bee behaviors, it remains unclear whether these genes actually regulate social behaviors in the honey bee, in part because of the scarcity of genetic manipulation methods available for application to the honey bee. In this review, we describe the genetic methods applied to studies of the honey bee, ranging from classical forward genetics to recently developed gene modification methods using transposon and CRISPR/Cas9. We then discuss future functional analyses using these genetic methods targeting genes identified by the preceding research. Because no particular genes or neurons unique to social insects have been found yet, further exploration of candidate genes/neurons correlated with sociality through comprehensive analyses of mushroom bodies in the aculeate species can provide intriguing targets for functional analyses, as well as insight into the molecular and neural bases underlying social behaviors.

Kenyon Cell Subtypes/Populations in the Honeybee Mushroom Bodies: Possible Function Based on Their Gene Expression Profiles, Differentiation, Possible Evolution, and Application of Genome Editing.

Mushroom bodies (MBs), a higher-order center in the honeybee brain, comprise some subtypes/populations of interneurons termed as Kenyon cells (KCs), which are distinguished by their cell body size and location in the MBs, as well as their gene expression profiles. Although the role of MBs in learning ability has been studied extensively in the honeybee, the roles of each KC subtype and their evolution in hymenopteran insects remain mostly unknown. This mini-review describes recent progress in the analysis of gene/protein expression profiles and possible functions of KC subtypes/populations in the honeybee. Especially, the discovery of novel KC subtypes/populations, the "middle-type KCs" and "KC population expressing FoxP," necessitated a redefinition of the KC subtype/population. Analysis of the effects of inhibiting gene function in a KC subtype-preferential manner revealed the function of the gene product as well as of the KC subtype where it is expressed. Genes expressed in a KC subtype/population-preferential manner can be used to trace the differentiation of KC subtypes during the honeybee ontogeny and the possible evolution of KC subtypes in hymenopteran insects. Current findings suggest that the three KC subtypes are unique characteristics to the aculeate hymenopteran insects. Finally, prospects regarding future application of genome editing for the study of KC subtype functions in the honeybee are described. Genes expressed in a KC subtype-preferential manner can be good candidate target genes for genome editing, because they are likely related to highly advanced brain functions and some of them are dispensable for normal development and sexual maturation in honeybees.

Detection of Phospholipase C Activity in the Brain Homogenate from the Honeybee.

The honeybee is a model organism for evaluating complex behaviors and higher brain function, such as learning, memory, and division of labor. The mushroom body (MB) is a higher brain center proposed to be the neural substrate of complex honeybee behaviors. Although previous studies identified genes and proteins that are differentially expressed in the MBs and other brain regions, the activities of the proteins in each region are not yet fully understood. To reveal the functions of these proteins in the brain, pharmacologic analysis is a feasible approach, but it is first necessary to confirm that pharmacologic manipulations indeed alter the protein activity in these brain regions. We previously identified a higher expression of genes encoding phospholipase C (PLC) in the MBs than in other brain regions, and pharmacologically assessed the involvement of PLC in honeybee behavior. In that study, we biochemically tested two pharmacologic agents and confirmed that they decreased PLC activity in the MBs and other brain regions. Here, we present a detailed description of how to detect PLC activity in honeybee brain homogenate. In this assay system, homogenates derived from different brain regions are reacted with a synthetic fluorogenic substrate, and fluorescence resulting from PLC activity is quantified and compared between brain regions. We also describe our evaluation of the inhibitory effects of certain drugs on PLC activity using the same system. Although this system is likely affected by other endogenous fluorescence compounds and/or the absorbance of the assay components and tissues, the measurement of PLC activity using this system is safer and easier than that using the traditional assay, which requires radiolabeled substrates. The simple procedure and manipulations allow us to examine PLC activity in the brains and other tissues of honeybees involved in different social tasks.

mKast is dispensable for normal development and sexual maturation of the male European honeybee.

The European honeybee (Apis mellifera L.) exhibits various social behaviors. The molecular and neural mechanisms underlying these behaviors have long been explored, but causal relations between genes or neurons and behaviors remain to be elucidated because effective gene manipulation methods in the honeybee have not been available until recently. We recently established a basic technology to produce mutant honeybee drones using CRISPR/Cas9. Here we produced mutant drones using CRISPR/Cas9 targeting mKast, which is preferentially expressed in a certain subtype of class I Kenyon cells that comprise the mushroom bodies in the honeybee brain. By immunoblot analysis, we showed that mKast protein expression was completely lost in the mutant drone heads. In addition, during the production process of homozygous mutant workers, we demonstrated that heterozygous mutant workers could be produced by artificial insemination of wild-type queens with the sperm of mutant drones, indicating that mKast mutant drones were sexually mature. These results demonstrate that mKast is dispensable for normal development and sexual maturation in drone honeybees, and allow us to proceed with the production of homozygous mutant workers for the analysis of a particular gene by gene knockout in the future.

Shotgun proteomics deciphered age/division of labor-related functional specification of three honeybee (Apis mellifera L.) exocrine glands.

The honeybee (Apis mellifera L.) uses various chemical signals produced by the worker exocrine glands to maintain the functioning of its colony. The roles of worker postcerebral glands (PcGs), thoracic glands (TGs), and mandibular glands (MGs) and the functional changes they undergo according to the division of labor from nursing to foraging are not as well studied. To comprehensively characterize the molecular roles of these glands in workers and their changes according to the division of labor of workers, we analyzed the proteomes of PcGs, TGs, and MGs from nurse bees and foragers using shotgun proteomics technology. We identified approximately 2000 proteins from each of the nurse bee or forager glands and highlighted the features of these glands at the molecular level by semiquantitative enrichment analyses of frequently detected, gland-selective, and labor-selective proteins. First, we found the high potential to produce lipids in PcGs and MGs, suggesting their relation to pheromone production. Second, we also found the proton pumps abundant in TGs and propose some transporters possibly related to the saliva production. Finally, our data unveiled candidate enzymes involved in labor-dependent acid production in MGs.

Pharmacologic inhibition of phospholipase C in the brain attenuates early memory formation in the honeybee (Apis mellifera L.).

Although the molecular mechanisms involved in learning and memory in insects have been studied intensively, the intracellular signaling mechanisms involved in early memory formation are not fully understood. We previously demonstrated that phospholipase C epsilon (PLCe), whose product is involved in calcium signaling, is almost selectively expressed in the mushroom bodies, a brain structure important for learning and memory in the honeybee. Here, we pharmacologically examined the role of phospholipase C (PLC) in learning and memory in the honeybee. First, we identified four genes for PLC subtypes in the honeybee genome database. Quantitative reverse transcription-polymerase chain reaction revealed that, among these four genes, three, including PLCe, were expressed higher in the brain than in sensory organs in worker honeybees, suggesting their main roles in the brain. Edelfosine and neomycin, pan-PLC inhibitors, significantly decreased PLC activities in homogenates of the brain tissues. These drugs injected into the head of foragers significantly attenuated memory acquisition in comparison with the control groups, whereas memory retention was not affected. These findings suggest that PLC in the brain is involved in early memory formation in the honeybee. To our knowledge, this is the first report of a role for PLC in learning and memory in an insect.

Increased complexity of mushroom body Kenyon cell subtypes in the brain is associated with behavioral evolution in hymenopteran insects.

In insect brains, the mushroom bodies (MBs) are a higher-order center for sensory integration and memory. Honeybee (Apis mellifera L.) MBs comprise four Kenyon cell (KC) subtypes: class I large-, middle-, and small-type, and class II KCs, which are distinguished by the size and location of somata, and gene expression profiles. Although these subtypes have only been reported in the honeybee, the time of their acquisition during evolution remains unknown. Here we performed in situ hybridization of tachykinin-related peptide, which is differentially expressed among KC subtypes in the honeybee MBs, in four hymenopteran species to analyze whether the complexity of KC subtypes is associated with their behavioral traits. Three class I KC subtypes were detected in the MBs of the eusocial hornet Vespa mandarinia and the nidificating scoliid wasp Campsomeris prismatica, like in A. mellifera, whereas only two class I KC subtypes were detected in the parasitic wasp Ascogaster reticulata. In contrast, we were unable to detect class I KC subtype in the primitive and phytophagous sawfly Arge similis. Our findings suggest that the number of class I KC subtypes increased at least twice - first with the evolution of the parasitic lifestyle and then with the evolution of nidification.

Correction: Correction: Gene expression and immunohistochemical analyses of mKast suggest its late pupal and adult-specific functions in the honeybee brain.

[This corrects the article DOI: 10.1371/journal.pone.0183522.].