SUMMARYWorker honey bees undergo a socially regulated, highly stable lipid loss as part of their behavioral maturation. We used largescale transcriptomic and proteomic experiments, physiological experiments and RNA interference to explore the mechanistic basis for this lipid loss. Lipid loss was associated with thousands of gene expression changes in abdominal fat bodies. Many of these genes were also regulated in young bees by nutrition during an initial period of lipid gain. Surprisingly, in older bees, which is when maximum lipid loss occurs, diet played less of a role in regulating fat body gene expression for components of evolutionarily conserved nutrition-related endocrine systems involving insulin and juvenile hormone signaling. By contrast, fat body gene expression in older bees was regulated more strongly by evolutionarily novel regulatory factors, queen mandibular pheromone (a honey bee-specific social signal) and vitellogenin (a conserved yolk protein that has evolved novel, maturationrelated functions in the bee), independent of nutrition. These results demonstrate that conserved molecular pathways can be manipulated to achieve stable lipid loss through evolutionarily novel regulatory processes. Supplementary material available online at

Depression symptoms following immune response to a challenge have been reported after the recovery from sickness. A RNA-Seq study of the dysregulation of the microglia transcriptome in a model of inflammation-associated depressive behavior was undertaken. The transcriptome of microglia from mice at day 7 after Bacille Calmette Guérin (BCG) challenge was compared to that from unchallenged Control mice and to the transcriptome from peripheral macrophages from the same mice. Among the 562 and 3,851 genes differentially expressed between BCG-challenged and Control mice in microglia and macrophages respectively, 353 genes overlapped between these cells types. Among the most differentially expressed genes in the microglia, serum amyloid A3 (Saa3) and cell adhesion molecule 3 (Cadm3) were over-expressed and coiled-coil domain containing 162 (Ccdc162) and titin-cap (Tcap) were under-expressed in BCG-challenged relative to Control. Many of the differentially expressed genes between BCG-challenged and Control mice were associated with neurological disorders encompassing depression symptoms. Across cell types, S100 calcium binding protein A9 (S100A9), interleukin 1 beta (Il1b) and kynurenine 3-monooxygenase (Kmo) were differentially expressed between challenged and control mice. Immune response, chemotaxis, and chemokine activity were among the functional categories enriched by the differentially expressed genes. Functional categories enriched among the 9,117 genes differentially expressed between cell types included leukocyte regulation and activation, chemokine and cytokine activities, MAP kinase activity, and apoptosis. More than 200 genes exhibited alternative splicing events between cell types including WNK lysine deficient protein kinase 1 (Wnk1) and microtubule-actin crosslinking factor 1(Macf1). Network visualization revealed the capability of microglia to exhibit transcriptome dysregulation in response to immune challenge still after resolution of sickness symptoms, albeit lower than that observed in macrophages. The persistent transcriptome dysregulation in the microglia shared patterns with neurological disorders indicating that the associated persistent depressive symptoms share a common transcriptome basis.

Behavioral indicators in the murine Bacille Calmette Guérin (BCG) model of inflammation have been studied individually; however, the variability of the behaviors across BCG levels and the mouse-to-mouse variation within BCG-treatment group are only partially understood. The objectives of this study were: 1) to gain a comprehensive understanding of sickness and depression-like behaviors in a BCG model of inflammation using multivariate approaches, and 2) to explore behavioral differences between BCG-treatment groups and among mice within group. Adult mice were challenged with either 0mg (saline), 5mg or 10mg of BCG (BCG-treatment groups: BCG0, BCG5, or BCG10, respectively) at Day 0 of the experiment. Sickness indicators included body weight changes between Day 0 and Day 2 and between Day 2 and Day 5, and horizontal locomotor activity and vertical activity (rearing) measured at Day 6. Depression-like indicators included duration of immobility in the forced swim test and in the tail suspension test at Day 6 and sucrose consumption in the sucrose preference test at Day 7. The simultaneous consideration of complementary sickness and depression-like indicators enabled a more precise characterization of behavioral changes associated with BCG-treatment and of mouse-to-mouse variation, relative to the analysis of indicators individually. Univariate and multivariate analyses confirmed differences between BCG-treatment groups in weight change early on the trial. Significant differences between BCG-treatment groups in depression-like behaviors were still measurable after Day 5. The potential for multivariate models to account for the correlation between behavioral indicators and to augment the analytical precision relative to univariate models was demonstrated both for sickness and for depression-like indicators. Unsupervised learning approaches revealed the complementary information provided by the sickness and depression-like indicators considered. Supervised learning approaches using cross-validation confirmed subtle differences between BCG-treatment groups and among mice within group identified by the consideration of sickness and depression-like indicators. These findings support the recommendation for multivariate and multidimensional analyses of sickness and depression-like indicators to augment the systemic understanding of the behavioral changes associated with infection.

Microglia in the brain and macrophages in peripheral organs are cell types responsible for immune response to challenges. Indoleamine 2,3-dioxygenase 1 (IDO1) is an immunomodulatory enzyme of the tryptophan pathway that is expressed in the brain. The higher activity of IDO1 in response to immune challenge has been implicated in behavioral disorders. The impact of IDO1 depletion on the microglia transcriptome has not been studied. An investigation of the transcript networks in the brain microglia from IDO1-knockout (IDO1-KO) mice was undertaken, relative to peripheral macrophages and to wild-type (WT) mice under unchallenged conditions. Over 105 transcript isoforms were differentially expressed between WT and IDO1-KO within cell type. Within microglia, Saa3 and Irg1 were over-expressed in IDO1-KO relative to WT. Within macrophages, Csf3 and Sele were over-expressed in IDO1-KO relative to WT. Among the genes differentially expressed between strains, enriched biological processes included ion homeostasis and ensheathment of neurons within microglia, and cytokine and chemokine expression within macrophages. Over 11,110 transcript isoforms were differentially expressed between microglia and macrophages and of these, over 10,800 transcripts overlapped between strains. Enriched biological processes among the genes over- and under-expressed in microglia relative to macrophages included cell adhesion and apoptosis, respectively. Detected only in microglia or macrophages were 421 and 43 transcript isoforms, respectively. Alternative splicing between cell types based on differential transcript isoform abundance was detected in 210 genes including Phf11d, H2afy, and Abr. Across strains, networks depicted a predominance of genes under-expressed in microglia relative to macrophages that may be a precursor for the different response of both cell types to challenges. The detected transcriptome differences enhance the understanding of the role of IDO1 in the microglia transcriptome under unchallenged conditions.

Comprehensive and simultaneous analysis of all genes in a biological sample is a capability of RNA-Seq technology. Analysis of the entire transcriptome benefits from summarization of genes at the functional level. As a cellular response of interest not previously explored with RNA-Seq, peritoneal macrophages from mice under two conditions (control and immunologically challenged) were analyzed for gene expression differences. Quantification of individual transcripts modeled RNA-Seq read distribution and uncertainty (using a Beta Negative Binomial distribution), then tested for differential transcript expression (False Discovery Rate-adjusted p-value < 0.05). Enrichment of functional categories utilized the list of differentially expressed genes. A total of 2079 differentially expressed transcripts representing 1884 genes were detected. Enrichment of 92 categories from Gene Ontology Biological Processes and Molecular Functions, and KEGG pathways were grouped into 6 clusters. Clusters included defense and inflammatory response (Enrichment Score = 11.24) and ribosomal activity (Enrichment Score = 17.89). Our work provides a context to the fine detail of individual gene expression differences in murine peritoneal macrophages during immunological challenge with high throughput RNA-Seq.