Abstract
Contrasting almost all other mammalian wintering strategies, Eurasian common shrews, Sorex araneus, endure winter by shrinking their brain, skull, and most organs, only to then regrow to breeding size the following spring. How such tiny mammals achieve this unique brain size plasticity while maintaining activity through the winter remains unknown. To discover potential adaptations underlying this trait, we analyzed seasonal differential expression in the shrew hypothalamus, a brain region that both regulates metabolic homeostasis and drastically changes size and compared hypothalamus expression across species. We discovered seasonal variation in suites of genes involved in energy homeostasis and apoptosis, shrew-specific upregulation of genes involved in the development of the hypothalamic blood brain barrier and calcium signaling, as well as overlapping seasonal and comparative gene expression divergence in genes implicated in the development and progression of human neurological and metabolic disorders, including CCDC22, FAM57B, and GPR3. With high metabolic rates and facing harsh winter conditions, Sorex araneus have evolved both adaptive and plastic mechanisms to sense and regulate its energy budget. Many of these expression changes mirrored those identified in human neurological and metabolic disease, highlighting the interactions between metabolic homeostasis, brain size plasticity, and longevity.
Introduction
Typical mammalian brain development consists of unidirectional postnatal growth until reaching an adulthood maximum (1, 2), but the Eurasian common shrew, Sorex araneus, seasonally changes the size of its body, most organs, skull, and —especially— its brain (3–6). In the most acute case of Dehnel’s phenomenon (DP), or seasonal size plasticity, the common shrew first grows to an initial summer maximum as a juvenile, then reduces size in autumn losing up to 18% body mass, ∼20% skull volume, and 26% brain mass, reaching its winter size minimum (3–6). Then, in the spring prior to mating, these shrews either partially or fully regrow these organs, reaching a second maximum. DP is hypothesized to be a plastic adaptation to decrease energy demand without slowing metabolism, improving shrew fitness during the low temperatures and limited food supply of winter (3, 4, 6–12). Corroborating this adaptive plasticity hypothesis, DP anticipates winter conditions and is both modulated by environmental temperature (3) and geographically variable (4). Although we have begun to reveal the molecular mechanisms of DP (12), the comparative context of proposed adaptations, and commonality of processes across brain regions remain unknown. Elucidating both the genetic basis and evolutionary mechanisms underlying DP will illuminate the adaptations involved in this unusual case of mammalian brain plasticity.
As a hypothesized wintering strategy, energy metabolism is key to the evolution of DP in S. araneus (11). Tiny mammals have limited options to survive winter conditions (13), and DP is exceedingly rare as a wintering strategy. While many mammals seasonally migrate away from low-temperature, low- productivity environments, some endure these conditions by reducing their energy requirements through hibernation; a phenotype that has independently evolved many times (14–16). The evolution and regulation of hibernation relies upon a suite of metabolic genes (17, 18), including those associated with thermogenesis, circadian rhythms, and feeding behavior. Compared to most species, S. araneus has astronomically large energetic demands, with one the highest basal metabolic rates per unit of body mass identified in any mammal (19, 20). To meet these energy requirements, Sorex araneus forage and feed constantly. Together with brief lifespans (∼1 year) that nonetheless require surviving a single winter (21, 22), these extreme metabolic demands constrain survival options for S. araneus. Thus, by decreasing resources devoted to energetically expensive tissue such as the brain (23), DP bypasses these constraints and expands the typical niche of a small mammal during winter while allowing this fast-living, predatory shrew to remain active year round.
Critical metabolic shifts have been implicated in the regulation of seasonal size plasticity at a molecular level (10, 12, 24). Previous stage-by-stage analyses of DP characterized gene expression and network topology in three regions of the body that undergo seasonal size change: both the cortex and hippocampus in the brain, and the metabolically active liver (12). While those analyses found decreased cholesterol efflux in the cortex and hippocampus during brain shrinkage, there were also profound metabolic changes across the body. Liver expression revealed reversible regulation of key metabolic transcription factors (FOXOs, PPARs, RXRs), indicating a shift away from lipid- toward glucose metabolism during autumn shrinkage. Metabolomic data supported this pattern, and a potential internal circadian clock regulating metabolic changes was also identified. Together, those results indicated a brain-liver crosstalk is a pivotal avenue of communication to coordinate metabolic modifications underlying seasonal size change across organs. Although disruption of the brain-liver axis has been implicated in neurodegenerative diseases in other mammals (25–28), how a highly functional brain-liver crosstalk which coordinates the shrew’s natural seasonal size change operates and evolved remains unknown.
We hypothesize DP, which involves metabolic shifts and inter-organ communication (12), has evolved through the combination of both physiological plasticity and genetic adaptation. This idea is twofold; phenotypic plasticity is not only adaptive in itself as the evolved spectrum of a trait, but can also allow populations of a species to persist in varying environments, leading to ancillary evolutionary adaptation when given enough time and strength of selection (29). Thus, while temporal analysis of shrew expression can elucidate putatively adaptive plastic regulation, it will miss adaptive canalization of gene expression that has contributed to the evolution of DP. Therefore, plastic regulatory mechanisms and proposed adaptations can be further tested via cross-species comparisons of expression divergence. These cross-species comparisons treat expression as a trait and have been used to characterize diverse vertebrate adaptations including thermal tolerance and vision loss in fish (30–32), the effects of whole genome duplication on expression (33–37), sex specific alternative splicing in birds (38), and patterns of expression evolution in mammalian organs (39–42). By incorporating these comparative approaches with seasonally varying gene expression, we can further elucidate the role of plasticity and adaption in both the regulation and evolution of DP.
The hypothalamus is an intriguing candidate brain region for comparative analysis, as it is both the brain region with the most intense seasonal size change in shrews, and the center of bodily homeostatic maintenance across mammals. In shrews, the hypothalamus undergoes drastic seasonal size change, with a 31.6% volume reduction in autumn through winter, followed by a 47.8% increase in spring (43). As in a previous study (12), we analyze hypothalamic gene expression across size change and compare these molecular processes of brain size change to those identified in the hippocampus and cortex, such as decreased cholesterol efflux or parallels to neurological disease. However, across mammals, the hypothalamus plays a pivotal role in the maintenance of an energy budgets, with functions influencing: 1) energy intake and feeding behavior, 2) energy expenditure and metabolic rate, and 3) energy deficits and storage. While such functions are important for all mammals, these are critical for wintering mammals to maintain stable internal body temperatures as external temperatures decrease. For example, the hypothalamus can activate the sympathetic nervous system as temperatures decrease, stimulating thermogenesis in winter. Thus, shrews may deploy typical mammalian hypothalamic plasticity through DP. Alternatively, and considering how rare DP is, shrews may have evolved divergent approaches, akin to naked mole rat traits, where this species lives under consistent hypoxic conditions resulting in adaptations to decreased oxygen and blood flow across the brain (44). Comparative approaches, then, may be able identify novel adaptations that contribute to the evolution of DP in the shrew hypothalamus.
As the brain region able to quickly respond to change, the hypothalamus has remarkable plastic capabilities associated with metabolism (45), especially as environmental conditions vary between seasons and can be unpredictable. While the central nervous system (CNS) acts somewhat independently of bodily metabolism through the protective effects of the blood brain barrier (BBB), the CNS is still sensitive to bodily metabolic change. We hypothesize hypothalamic plasticity is central to the evolution and regulation of DP through both environmental sensing mechanisms and signaling responses to these stimuli. Specifically, we hypothesize adaptations in the shrew: 1) hypothalamic BBB associated with dynamic metabolic fluctuations, 2) sensing of metabolic state and seasonal size through hormonal signals, mediated by BBB-crossing of insulin, ghrelin, and leptin, and 3) responses to metabolic fluctuations involving ion-dependent signaling. With limited energy inertia, shrews must continuously sense their peripheral metabolism, and a combination of receptors and selective BBB permeability in the hypothalamus allows certain molecules to relay information from the peripherals to the brain. Among the best characterized of these are the peripheral hormones insulin, ghrelin, and leptin, the latter known to excite anorexigenic and inhibit orexigenic neurons in the arcuate nucleus of the hypothalamus. Finally, according to Ramon y Cajal’s neuronal doctrine (46), differences in brain functional responses may stem from altered synaptic firing. Thus, inter- and intracellular ion concentrations play a large role in the communication between neuronal networks by propagating signals in response to perceived environmental stimuli.
To test these hypotheses, we analyzed both the seasonal and phylogenetic variation in shrew hypothalamic expression to detect signals of adaptive plasticity. First, we aimed to identify differential expression of genes across five stages of DP that might promote regulatory responses to seasonal variation and can be functionally validated with cell line perturbations. Second, this analysis was paired with a comparative transcriptomics approach using hypothalamic expression data from 15 additional mammal species. These analyses infer adaptation by testing for branch-specific expression shifts using Ornstein-Uhlenbeck models. The objective of these evolutionary analyses was to quantify lineage-specific hypothalamic expression changes in S. araneus, with evolutionary expression divergence consistent with selection. Finally, by comparing individual genes and associated pathways from these two analyses, we can determine potential adaptive plasticity, indicating mechanisms that both potentially regulate DP and were selected for higher expression in the evolution of DP. Our expression results implicate several key processes in DP, including seasonal plasticity in feeding behavior, adaptive modulation of both hypothalamic blood brain barrier and downstream signaling, and plastic apoptosis responses.
Results
RNA sequencing, mapping, and quantification
Between and within species, measurements of RNA quality and alignment, quantification, and normalization procedures were suitable for analyzing differential evolutionary and temporal expression associated with DP (Supplemental Data). First, a single “hypothalamus” from our autumn data was removed, as the expression profile resembled that of a cortex sample. Thus, our novel shrew hypothalamus sequences (n=23) consisted of five summer juveniles, three autumn juveniles, five winter juveniles, five spring adults, and five summer adults. RNA extracted across these seasons had a mean RNA Integrity Number (RIN) of 6.6 (range 5-7.8) and a mapping rate mean of 50.8% (range 45.1-57.4%). Although four samples had lower than initially expected (<6) RIN values, mapping rate (an indicator of sequencing quality) was only slightly impacted and resembled previously published cortex (mean 53.7%; range 48.1-61.5%) and hippocampus (mean 54.0%; range 48.1-61.6%) data. For evolutionary analyses, hypothalamic RNA sequencing data sets were identified for 19 species comprising 5 mammal orders.
Four species were removed from this experiment during filtration. Oryctolagus cuniculus and Pan paniscus were prepubescent, Phodopus sungorus was removed due to a low-quality genome assembly at the time of analysis, and Papio anubis was removed from this experiment because of its extremely low mapping percentage (mean 16.75%; range 9.5%-23.7%), which had a large effect on the count distribution. Remaining species had an average mapping rate of 56.4% (range 33.1-77.9%). Although RINs were not available for every species in the dataset, roughly similar mapping percentages suggest adequate RNA quality. Furthermore, normalization by library size should reduce potential biases from differences in mapping rate. We identified 6,496 single-copy orthologs found in all species with OrthoFinder and filtered normalized expression (TPMs) by these genes (oTPMs). Mean expression of orthologs across all species was 31.44 oTPM (range 20.66-44.91), with the shrew lineage having an average 39.42oTPM. The distributions of oTPM by species was visualized in Supplemental Figure 1, illustrating similar frequency distributions.
Temporal Expression
We identified several known neural signaling pathways correlated with patterns of seasonal hypothalamus expression that may be associated with changes during DP. We began by hierarchically clustering expression to identify patterns of gene expression through time. Most genes did not exhibit much variation through time, indicating constitutive gene expression in the hypothalamus, as only 786 of 19,296 genes passed our filters (>0.5-fold change between any two stages and 2 samples > 10 normalized reads). We identified 12 distinct clusters of gene expression patterns by bootstrapping (n=20) a gap statistic within cluster distances. Of these 12 clusters, five clusters consisting of 392 genes resembled a large divergence between summer juveniles and the remaining individuals (Supplemental Figure 2; Clusters 2, 3, 8, 11, 12). These genes likely represent a large developmental shift between recently postnatal shrews and the remainder of individuals, rather than a shift associated with seasonality or Dehnel’s phenomenon. We then functionally characterized the remaining 394 genes (Figure 1A), which exhibited seasonal shifts, with a GO pathway enrichment using the DAVID Gene Functional Classification Tool. Although no enrichment pathway was significant after a Bonferroni correction, 14 pathways were enriched prior to correction (p < 0.05); with the five pathways with the lowest p-values including fluid shear stress and atherosclerosis, relaxin signaling, neuroactive ligand-receptor interaction, HIF-1 signaling, and PI3K-Akt signaling (Figure 1B). Many of these pathways have been implicated in various physiological processes, which suggests that variation identified in the hypothalamus are likely the result of autonomic processes in homeostatic maintenance.
We also discovered hundreds of differentially expressed genes between autumn (Stage 2) and spring (Stage 4) shrews that may mediate phenotypic divergence in both size (shrinkage vs. regrowth) and metabolic changes (liver expression shifts) associated with DP. By comparing the RNA expression of autumn and spring shrews, we found 333 differentially expressed genes, with 194 upregulated and 139 downregulated in spring (Figure 2A). Of these genes, 57 upregulated and 34 downregulated were differentially expressed to a higher degree (> absolute fold change 3). The hypothalamus had a similar number of differentially expressed genes during DP compared to previous experiments elucidating expression differences in the cortex (DEG=540) and hippocampus (DEG=266). We ran a GO pathway enrichment to functionally characterize the 333 DEGs in the hypothalamus and found only five pathways were enriched (p<0.05; Modified Fischer’s Exact Test) with significant genes prior to multiple comparison correction. These pathways include four upregulated pathways; Fanconi anemia (9.3 fold- enrichment, p<0.01), GABAergic synapse (5.6 fold-enrichment, p<0.05), spliceosome (4.3 fold- enrichment, p<0.05), and nicotine addiction (9.42 fold-enrichment, p<0.05), with one downregulated pathway; apoptosis (5.6 fold-enrichment, p<0.05) (Figure 2B). The apoptosis pathway consists of significantly downregulated genes (from DESeq2 Wald test), BCL2L1 (−0.5 Log-Fold Change [LFC] padj<0.05), NGF (-2.3 LFC, padj<0.05), NFKB1A (−2.3 LFC, padj<0.05), FOS (−1.7 LFC, padj<0.05), and CTSK (−1.2 LFC, padj<0.05), while the GABAergic pathways consists of upregulated GABRE (5.0 LFC, padj<0.01), GABRQ (4.4 LFC, padj<0.05), CACNA1D (0.6 LFC, padj<0.05), and GNB3 (3.5 LFC padj<0.05) (Figure 2C). Differential expression for all genes can be found in Supplemental Data. Overall, the processes identified suggest regulation of cell death and synaptic plasticity in the hypothalamus may be associated with phenotypic changes during DP.
Cell Viability Analysis with BCL2L1 Overexpression
We examined the potential functional effects of apoptosis regulating gene, BCL2L1, which was differentially expressed between seasons, by propagating an in vitro model of domesticated ferret brain cells (MPF-CRL-1656, ATCC). We chose this cell line because there is no established cell line of Sorex araneus. Mustela putorius furo also undergoes a Dehnel’s-like phenotype (47) and is more closely related to shrews than mice are. MPF cells were transfected with either the anti-apoptotic BCL2L1 RNA or a scrambled version of BCL2L1 (sBCL2L1) RNA. 24 hours after transfection, cells were exposed to heat stress followed by crystal violet staining. A significant difference in cell viability was observed between the heat scrambled compared to the control scrambled groups (df=3.56, t3.56=8.78, p<0.01), however, there was no significant difference between cells transfected with BCL2L1 compared to those transfected with sBCL2L1 when exposed to heat (df=3.37, t3.37=-0.99, p=0.39).
Evolutionary Divergence in Expression
Analyses quantifying branch-specific shifts in expression identified hundreds of genes that were upregulated in the shrew hypothalamus compared to other mammals, suggesting adaptation for increased expression in these genes (Figure 3). EVE models found 222 genes significantly rejected (padj < 0.05) a single expression optimum for all species in favor of a second expression optimum for the shrew lineage. While nonsignificant genes had a mean 0.93-fold change, genes experiencing an expression branch-shift in shrews had a mean 9.66 oTPM fold change compared to the other species. To validate genes with expression branch-shifts as shrew-specific, we ran a dropout test for the shrew lineage, testing for phylogeny-wide evidence of selection for each gene. Of the 6,496 genes tested, 81 genes had significantly (padj < 0.05) lower β (within population variance to between species variance), suggesting selection elsewhere in the phylogeny. None of these genes overlapped with those identified to be differentially expressed in the shrew, further validating the specificity of expression shifts.
We identified several changes associated with processes that may underlie DP. Pathway enrichment analysis identified two significantly enriched pathways: calcium signaling (2.8 fold-enrichment, p<0.05) and autophagy (9.1 fold-enrichment, p<0.05) (Figure 3). Upon manual inspection of the list, we also identified several other potential adaptive processes, including blood-brain-barrier formation and function, feeding behavior and leptin responses, metabolism, neuroprotection, and GABAergic neuron development (Figure 3, Table 1).
Many of these putative adaptations mediate responses to environmental cues centered around energy demands. Five of the genes experiencing branch-shift changes in expression were also differentially expressed between autumn and spring individuals, suggesting not only an adaptive shift in the shrew lineage associated with the evolution of Dehnel’s phenomenon, but also a direct molecular mechanism for brain changes. The five genes that overlap between both analyses include (Figure 4): CCDC22, which plays an important role in endosomal recycling of membrane proteins (48); FAM57B, which mediates synaptic cell membrane architecture and function (49), GPR3, an orphan G protein-coupled receptor linked to both Alzheimer’s disease (50) and obesity (51); LMX1A, a transcription factor essential for dopaminergic neuron development (52); and PAQR4, which appears to regulate growth and apoptosis through sphingolipid synthesis in human cancers (53).
Discussion
By characterizing both seasonal and between-species differential expression of the hypothalamus, we generated and probed a unique data set, discovering expression shifts associated with extreme brain size change in Sorex araneus. While the focus on seasonal hypothalamus expression provides insights into both metabolic regulation and processes underlying brain size change, evolutionary shifts in expression can reveal the underpinnings of mammalian brain degeneration and regeneration in a natural system. Our analyses identified a suite of genes related to energy homeostasis that were both seasonally plastic and adaptively upregulated in the shrew, reinforcing the role of metabolism in both the evolution and regulation of DP. We also found shrew-specific upregulation of genes associated with the development of the hypothalamic blood brain barrier, which we hypothesize improves the metabolic sensing capabilities of the hypothalamus, further highlighting the importance of brain-liver crosstalk in promoting seasonal size change (12). Other molecular mechanisms, including adaptive calcium signaling and plastic apoptotic responses, were also implicated in DP. Many of these results, including genes that were both seasonally plastic and evolutionarily upregulated, resemble expression changes found in human neurological and metabolic disease, and thus may prove important for therapeutic treatment of those diseases.
Plastic Adaptations of Metabolism, Feeding Behavior, and Leptin Responses
Seasonal expression variation in the shrew hypothalamus was largely associated with the maintenance of homeostasis, especially processes related to energy balance, corroborating the metabolic hypothesis of DP evolution in the shrew. Although DP has been proposed to reduce shrew energy requirements by limiting energy devoted to the maintenance of larger tissues (4, 7, 20), the extent of system-wide metabolic shifts during shrinkage and regrowth, particularly in the liver, has only recently been revealed (12). We found three pathways enriched for genes with seasonal variation in the hypothalamus which also have metabolic associations: relaxin signaling, phosphoinositide-3 kinase (PI3K) signaling, and neuron ligand receptor interactions. Relaxin, a neuropeptide, plays a pivotal role in many physiological functions in model organisms, including food and water intake (54–57); food restriction upregulates relaxin-3 (RLN3) in vivo (58), while administering RLN3 protein into the rat hypothalamus vastly increases food intake (59–62).
Relaxin has also been found to activate the PI3K pathway (63, 64), which mediates feeding behavior and glucose homeostasis (65–69), as well as angiogenesis via vascular endothelial growth factor (VEGF) production (64, 70, 71). PI3Ks activity influences feeding behavior and weight gain in response to key metabolic hormones, including leptin and insulin (67–69, 72). Lastly, the neuroactive ligand-receptor interaction pathway, a subclass of KEGG Environmental Information Processing pathway, is associated with responses in model organisms to shifts in diet (73), stress (74), and temperature (75, 76), all of which the shrew experiences during seasonal changes. Although exploration of gene function in S. araneus has only begun, seasonal expression in the hypothalamus indicates maintenance of metabolic homeostasis, likely reflecting an extreme physiology that pairs disproportionately high metabolic rates (20) with lifespans shorter than expected by body size (∼1 year) (21, 22). Facing scarcity in the winter, yet maintaining costly activity in their environment, shrews must efficiently regulate feeding behavior, food intake, and energy homeostasis, and we find signatures of such regulation in their hypothalami.
We also found evidence for adaptive expression relating to maintenance of energy balance, specifically in feeding behavior and leptin responses, further emphasizing the importance of selection to meet high energetic demands in wintering environments in the evolution of DP. Key genes associated with the hypothalamic leptin response, such as NPY, AgRP, or POMC, were not universally orthologous and so could not be analyzed, but there was evolutionary upregulation of many genes involved in downstream leptin responses (FOXO4, GRB2, IRX5, AGT) (Figure 3, Table 1). Unlike their upstream effectors, these genes may have subtle effects on continuous feeding behavior. In the hypothalamus, FOXO1 negatively regulates POMC and promotes AgRP expression (77), which can contribute to leptin resistance resulting in a long term appetite suppression (78–80). However, much less is known about the hypothalamic functions of the FOXO1 paralog, FOXO4. In mouse adipocytes, leptin- dependent expression of FOXO4 can rapidly clear blood glucose levels (81), which can lead to energy storage and decreased satiety. We hypothesize that, in shrews, upregulation of FOXO4 could either increase blood glucose clearance or leptin resistance, reducing long term satiation and promoting continuous feeding. Further supporting this hypothesis, we also found shrew-specific upregulation of other genes associated with leptin sensitivity, such as GRB2, which in hamsters regulates seasonal leptin sensitivity and body weight (82), and IRX5, whose mice knockouts exhibit decreased weight and food intake due to an enhanced hypothalamic leptin response (83, 84). Functionally, evolutionary upregulation of these genes in the shrew hypothalamus may result in evolved leptin insensitivity, instead of the development of leptin resistance. Reduced leptin sensitivity in shrews could reduce overall satiation, promoting foraging even when the short-term energy budget is balanced thereby improving anticipatory storage before winter.
Evolution of BBB Permeability and Energy Sensing
Comparative analyses also indicate shrews may have evolved adaptive environmental sensing. We found upregulation of several proangiogenic genes associated with the formation and function of the blood brain barrier, including VEGFA and WNT7A. In normal neural development, neural progenitor cells and surrounding astro- and pericytes express these two genes defining an expression gradient that guides blood vessels to form the blood brain barrier (85–87).
While this suggests VEGFA and WNT7A upregulation in the shrew may increase vascularization of the brain, continued expression of VEGFA beyond development increases BBB permeability, leading to a breakdown of its protective effects (88, 89). BBB permeability is a common symptom of many disease- like states including neuroinflammation (90) and neurodegeneration (90–92). The hypothalamus BBB consists of tanycytes and highly fenestrated capillaries with less tight junctions (93), replacing the less permeable BBB of other brain regions with a specialized barrier that improves nutrient and energy sensing through dynamic passage of molecules from peripheral circulation (94–96). Experiments on mice indicate reduced blood glucose levels from fasting promote capillary fenestration in the hypothalamic BBB via upregulation of VEGFA, eliciting feeding responses through increased hypothalamic exposure to glucose levels (94). In shrews, adaptive upregulation of VEGFA and WNT7A may constitutively increase capillary fenestration, an adaptation of the BBB to improve hypothalamus metabolic sensing and signaling in response to high energy demands. Notably, adaptive expression of VEGFA has also been identified in the seasonal vascular plasticity of other mammal species, including reproductive sheep (97, 98), bison (99), and squirrels (100). Reproduction is an extremely metabolically demanding process, thus, constitutive shrew specific upregulation of VEGFA to meet unusually high metabolic demands may parallel mechanisms used in seasonal reproduction.
Calcium Signaling and Apoptosis
Calcium signaling is a ubiquitous form of neural communication (101, 102) that occurs through the regulation of intra and extracellular calcium concentrations, transmitting signals that can release neurotransmitters (103, 104) and promote gene expression (105–108). Thus, it has been implicated in numerous functions including dendrite growth (109, 110), synaptogenesis (111–113), and initiation and maintenance of long-term potentiation (114, 115). In shrews, adaptive calcium signaling through the upregulation of several calcium-responsive genes (VEGFA, CASQ1, HRH2, TACR2, PPIF, SPHK2, GRIN1) may improve environmental information relay in the hypothalamus, both internally and throughout the central nervous system.
We found that HRH2 is upregulated in the hypothalamus of shrews, while maintaining strong sequence homology with other mammals. Although adaptive calcium signaling could have numerous physiological effects associated with shrew fitness or DP, many neuroadaptations (116, 117) have pinpointed the effects of calcium signaling by identifying the local functions of calcium responsive molecules and genes (118, 119). One such neuroadaptation involves the histamine receptor HRH2, which depresses neuronal action potential when bound to histamine by blocking calcium-dependent potassium channels (120, 121). In the hypothalamus, HRH2 and other histamine receptors regulate circadian rhythms and feeding behavior (122). Alternatively, agonism of HRH1 in rodents, an excitatory receptor in contrast to HRH2 (123), reduces food consumption by ∼50% (124), while a genetic knockout of HRH1 both increases feeding and disrupts circadian feeding behavior (125). Data on HRH2 are sparse in comparison to those on HRH1 (126), but agonism of HRH2 did not significantly affect rodent feeding behavior (124) and mice HRH2 knockouts reduce nocturnal mice activity (126). We propose the ∼12-fold upregulation of HRH2 receptors in shrews, despite being evolutionary divergent from rodents, also disrupts light entrainment, with the potential to increase food intake.
Calcium signaling also regulates opposing processes, including both cell survival and programmed cell death (118), and thus evolutionary upregulation of genes involved in these processes may be associated with seasonal brain shrinkage and regrowth. Autophagy is the degradation of cellular parts to provide energy and materials to maintain cellular homeostasis under stressful conditions. We found an evolutionary upregulation of autophagy genes in the shrew hypothalamus, including MLST8, an integral component of the mTOR protein complex 1 (mTORC1), which regulates sugar and lipid metabolism(127), but also has recently been found to modulate autophagy in response to nutrient deprivation (128). However, sensing calcium levels is also vital for the regulation of apoptosis and efflux of damaged or dying cells. Two evolutionarily upregulated genes PPIF and SPHK2, can promote apoptosis associated with mitochondrial calcium influx. PPIF is an important component of the mitochondrial permeability transition, which responds to cellular calcium overload and oxidative stress by promoting cell death (129–132). Similar to PPIF, overexpression of SPHK2 in stressed cells induces cell cycle arrest and promotes apoptosis through a series of signaling cascades that increase cytosolic calcium and promote the release of pro-apoptotic cytochrome C from the mitochondria (133, 134). While both SPHK2 and PPIF are overexpressed and likely contribute to neuronal death in both human patients with Huntington’s or Alzheimer’s disease and rodent models of these diseases, knockout or inhibition of these genes protects against neurodegeneration (135–138). In the shrew hypothalamus, canalized upregulation of these genes likely influences cell survival and death, but precisely how requires further investigation.
Validating the critical role of apoptosis regulation in Dehnel’s phenomenon, several apoptosis-regulating genes are differentially expressed between seasons. We found apoptosis pathway genes, including NGF, NFKBIA, and BCL2L1, are all upregulated during shrinkage and downregulated in regrowth, despite previous analyses that identified no changes in brain cell numbers in most brain regions through DP (139). Regulation of apoptosis may represent a convergent mechanism across wintering strategies, as natural experiments found seasonal reduction of apoptosis, including upregulation of BCL2L1, occurs in various organs of hibernators (140, 141), including the brain (142). Specific to DP, this result indicates tight regulation of cell survival during brain size changes and periods of metabolic and environmental stress: a balance between brain shrinkage and maintenance of cell numbers.
Preserving cell numbers despite size change requires both pro- and anti-apoptotic regulation. Known functions of genes cycling during DP bolster this view: NGF is an upstream activator of neuronal cell death through interactions with the p75 neurotrophin “death” receptor (143), but NGF can promote cell survival by activating the NFKB pathway that blocks apoptosis further downstream (144–146). For example, in primary rat hippocampal neurons, NGF promotes expression of anti-apoptotic BCL2L1 through an NFKB-dependent pathway (147) whose upregulation reduces apoptosis and promotes cell survival in rodents (148, 149). To test if autumn upregulation of these genes acts to reduce apoptosis, we overexpressed BCL2L1 in a related mammalian neural cell line, but we did not find a rescue against apoptosis from heat. This result defied our expectations, with several possible explanations.
Methodologically, heat-induced apoptosis may not accurately reflect the challenges to cell survival present in DP. Alternatively, the anti-apoptotic function of BCL2L1 may be specific to the unique neural environment of shrews, with apoptosis being propagated through alternative pathways or species- specific gene interactions. DP-related apoptosis regulation as an adaptation may be convergent across wintering mammals, but the regulatory mechanisms of this process through shrew size change require further investigation.
Disease-associated Adaptive Plasticity
A subset of genes, CCDC22, FAM57B, GPR3, LMX1A, and PAQR4, are both seasonally plastic and evolutionarily upregulated in the shrew hypothalamus (Figure 4). The combination of seasonal and evolutionary upregulation suggests these genes play integral roles in DP, likely through adaptive plasticity. Many of these genes are directly related to the development and progression of neurological disorders in humans. While some of these genes may modulate DP, others may be a byproduct of brain size change, or act as neuroprotection against the negative consequences of brain size and structure changes.
CCDC22
CCDC22 regulation in shrews may provide protection against harm incurred during seasonal size changes. The coiled-coil domain containing 22 protein (CCDC22) plays a key role in endosomal recycling of proteins and is a novel candidate gene for several neurological disorders. CCDC22 is a subunit of the CCC (CCDC22, CCDC93, COMMD) complex that aids in trafficking and recycling of endosomal membrane proteins (48, 150, 151). Improper recycling and resulting protein aggregation readily occurs in neurons, as they do not proliferate (152, 153), and has also been implicated in neurodegenerative diseases, including cytotoxicity in amyotrophic lateral sclerosis (154), amyloid-B aggregation in Alzheimer’s disease (155–157), and MUL1 degradation in Parkinson’s disease (158). Missense mutations and subsequent downregulation of CCDC22, are potential mechanisms underlying X- linked intellectual disability (159, 160) and Ritscher-Schinzel syndrome (160–162) in humans. Both disorders are characterized by macrocephaly, malformations of the brain, craniofacial abnormalities, and intellectual disabilities. In pubescent spring shrews, we found CCDC22 has a ∼7-fold upregulation compared to other mammals. This evolutionary upregulation may protect against neurodegeneration, as brain size and structure are altered during seasonal cycling. We also found a significant upregulation during spring compared to autumn shrinkage. Neuroprotective functions may be especially important in spring, as brain mass returns increasing the need for endosomal recycling of proteins and may explain seasonal plasticity in CCDC22 expression in tandem with adaptive upregulation.
FAM57B
Combined evolutionary upregulation and seasonal plasticity was also found in the expression of family with sequence similarity 57 member b gene, FAM57B, which has been associated with 16p11.2 deletion syndrome (163). This deletion leads to copy number variations and haploinsufficiency of 25 genes, including FAM57B (164), resulting in autism spectrum disorders (165). Disease models of this deletion have a variety of symptoms that parallel both shrew ecology and DP. In rodents, heterozygous 16p11.2 deletion leads to hyperactivity, craniofacial defects, reduced brain size, and altered brain morphology and shape (163, 166, 167), with duplications showing the opposite effects on weight, activity, and craniofacial morphology (166, 167). FAM57B haploinsufficiency in the more evolutionarily divergent zebrafish creates similarly substantively modified phenotypes, with heterozygous and homozygous knockouts of FAM57B resulting in abnormal axon development (168) and increased head and body size (164, 167).
Brain and body phenotypes related to FAM57B expression may be associated with the gene’s role in lipid metabolism and synaptic composition (169, 170). Across animals, the dosage and compensation effects of FAM57B expression have conserved impacts on metabolism, behavior, brain size, and body size. In the shrews, FAM57B has a higher baseline expression than in model organisms. We hypothesize upregulation of FAM57B in shrews contributes to seasonal plasticity of the skull, brain and body.
GPR3
As with CCDC22 and FAM57B, G-protein receptor 3, GPR3, influences both metabolism and neural development, but unlike them, its effects are age dependent. In mouse models of apoptosis, GPR3 has been associated with neurite outgrowth and maturation by promoting neuronal survival pathways (171). Despite these benefits, overexpression of GPR3 also been implicated in Alzheimer’s disease and contributes to amyloid beta generation (50, 172). The metabolic effects of GPR3 also indicate age- dependent benefits and outcomes. GPR3 knockout mice exhibit late-onset obesity associated with reduced adipose thermogenesis (51). If the benefits of GPR3 are age dependent, downregulation of this gene during brain regrowth in the shrew may be neuroprotective. Although evolutionary upregulation may increase neuronal cell survival during size change, plastic regulation is also important to reduce its negative consequences as shrews age. We propose GPR3 expression likely influences the adaptive and plastic processes important in DP, including apoptosis/cell survival, calcium signaling, and metabolism.
Conclusion
We discovered evidence for adaptive evolution of gene expression involved in the regulation of metabolism, cell survival, and size plasticity associated with DP. Not only does the hypothalamus shrink to reduce energetic load during winter, but it must also control metabolic homeostasis during these changes. Our results reveal numerous conserved and novel gene expression mechanisms that likely underlie the central functions of the hypothalamus during DP. While previous work had hinted at the metabolic nature of DP, by pairing seasonal differential expression with comparisons to other species, we identified adaptations likely involved in the evolution of this unique brain size plasticity. Metabolism, body size, and longevity are intrinsically linked life history factors. In the shrew, unusually high metabolic rates force winter activity, leading to intense metabolic sensing and regulation and changes in cell survival pathways, and –through both the evolution of DP and these sensing and regulatory adaptations– shortened lifespans. Together, these gene expression adaptations both likely underlie drastic seasonal size change and mitigate its detrimental effects.
Materials and Methods
Shrew Sample Collection
We used tissues from shrews collected for experiments analyzing gene expression and metabolic shifts in brain regions through DP (12). Briefly, S. araneus were collected from a single German population (47.9684N, 8.9761 E) across five different stages of DP from June 2020-June 2021 (n=25, nper_stage = 4-5). Shrews were trapped with insulated wooden traps that contained mealworms, which were checked every two hours. This protocol was used to minimize trap-related stress from heat/cold shock or lack of food and thus reduce stress-related variation on gene expression. Shrews were then euthanized via vascular perfusion of PAXgene Tissue Fixative. Brain regions were dissected into individual regions in cold fixative and then incubated in PAXgene Tissue Stabilizer for 2-24 hours before long-term storage.
Samples were preserved in stabilizer at −180°C in liquid nitrogen until RNA extraction.
RNA Extraction, Library Preparation, and Sequencing
We extracted the S. araneus hypothalamus and olfactory bulb using the same methods described for the cortex and hippocampus in previous work (12). These extractions reduce degradation due to heat and improve RIN from standard Qiagen Micro RNAeasy protocols by disrupting tissue on dry ice using glass, instead of plastic, mortar and pestles. A full list of changes to the Qiagen Micro RNAeasy protocol can be found in our previous manuscript (12). RNA was sent to Azenta Life Sciences for quality control (nanodrop and RNA ScreenTape), library preparation, and sequencing. Hypothalamus libraries were prepared with poly-A selection and sequenced for approximately 15-25 million reads per sample in 150bp PE reads.
RNA sequences from other species were collected from the NIH National Center for Biotechnology Information’s (NCBI) Sequence Read Archive data sets. First, we searched the Sequencing Read Archives (SRA) for the keyword hypothalamus, filtered by both RNA and mammal species. As models described below rely on using the variance both between and within species, species required three or more biological samples to be included. This removed any species that had only one hypothalamus sample sequenced. Species were excluded from the study if a genome assembly was not readily available, as this information was needed for an unbiased alignment of reads. Our final criterion was to only include species for which post-pubescent individuals were available to reduce noise from age effects and the onset of puberty. If the onset of puberty was not specifically stated for the dataset, AnAge was used to determine if samples were pubescent. For the remaining species, if multiple hypothalamus RNA-seq data sets were available, we used less stringent rules to determine which data to select. First, to reduce domestication or captivity effects, we chose wild or non-domesticated individuals over laboratory-raised animals. Second, to reduce the variance from different extraction protocols and sequencing methods, we used samples from the same experiments or larger data sets, available for multiple species or brain regions. Using samples from datasets with multiple regions or species also suggests intimate knowledge of neural anatomy, increasing confidence in dissections. Third, we used a maximum of eight samples per species and attempted to maintain a 1:1 sex ratio when possible. We did not filter by sequencing method (e.g., used 50bp SE as well as 150 PE reads) or sequencing depth, as filtering poor reads and normalizing for library depth and content should account for both these factors.
RNA Quantification, Normalization, and Orthology
Adapter sequences were trimmed from the raw reads using the default parameters of fastp (173), which is able to autodetect adapters regardless of different library preparations across species. This program also corrected and pruned low quality bases and reads using a sliding window approach. Processed reads were aligned and quantified with Kallisto 0.46.2 (174), which probabilistically estimates gene counts through pseudo-alignment to each species-specific genome assembly. Samples were removed from the data set if mapping rate was below 30%, indicating either poor sequencing or low assembly quality. Additionally, novel shrew hypothalamus sequences were compared against previously published shrew region data to verify tissue type. Gene counts for all species were then normalized for total library size into Transcripts Per Million (TPM). Orthologous genes between species were inferred with OrthoFinder (175) using default parameters, retaining only single-copy orthologs identified across all species. Visualizations of the frequency distribution of the orthologous (oTPMs; Supplemental Figure 1) were used to identify outlier species.
Branch-shift Changes in Expression
We modeled gene expression as a trait using phylogenetic comparative methods, as these account for evolutionary divergence among species (176). We ran Expression Variance and Evolution models (EVE) to test for significant changes of expression in the shrew lineage (39). EVE models both parameterize and estimate the ratio (β) of population (within species) to evolutionary (between species) variance, such that high β ratios indicate expression plasticity, and low β ratios indicate differential expression between species, necessary for inferring selection. We ran two EVE models to test for divergence in expression in the shrew lineage. First, we tested for branch-specific expression level shifts by contrasting the likelihood of two Ornstein-Uhlenbeck models for the data, one with a single expression optimum (null; stabilizing selection), and another with a second optimum on the shrew branch (selection). A likelihood ratio test between the null and selective hypotheses for each gene (χ21) was conducted using these likelihoods. A Bonferroni correction was used to account for multiple hypothesis tests to identify candidate genes under selection (padj<.05). For the second model, we ran a dropout test to further validate the specificity of the expression change in the shrew lineage. After removing the shrew expression data, we identified genes with high expression divergence across the phylogeny (significantly low β ratios), which would indicate relaxed or diversifying selection. This gene list was compared to that from first model to prune genes not specifically associated with divergence in the shrew. Both models described above ran on the Bayesian molecular-clock mammalian phylogeny of (177) pruned to match our species sample. Finally, we used the DAVID functional annotation tool (178) to determine enriched KEGG Gene Ontology (GO) pathways from the candidate gene set.
Temporal Clustering
To analyze the temporal variation of shrew expression in the hypothalamus, we temporally clustered our data using the package TCseq (179). We began by further normalizing our data across DP stages using the median of ratios from DESeq2 (180), which normalized for library size and content. Genes with consistent expression across time with little variation between stages would not be associated with observed phenotypic changes. These genes were filtered from this data set by only selecting genes with an absolute fold change of 0.5 between any two stages. We also removed genes with low expression that would appear to have high fold changes despite lacking enough transcript expression to influence phenotype, retaining only reads with two samples > 10 normalized reads. After filtering, counts for each stage were converted to mean z-scores and then clustered into groups of similar gene expression profiles using fuzzy cmeans clustering. The number of resulting gene clusters was calculated a priori, by bootstrapping (n=20) a gap statistic that minimizes within-cluster distances. Genes with highest membership in clusters associated with variation across all stages (seasonal effect; Clusters 1, 4, 5, 6, 7, 9, 10), as compared to those that just show divergent expression in Stage 1 (developmental effect; Clusters 2, 3, 8, 11, 12), were analyzed for KEGG GO enrichment using DAVID functional annotation (178).
Differential Gene Expression
We tested for differential gene expression between autumn (Stage 2) and spring (Stage 4) individuals, as these seasons differ in hypothalamic size phenotype (shrinkage vs. regrowth) (43), and were previously identified divergence in liver expression related to metabolism (12). To test for significant differential expression, we fit a negative binomial generalized linear model to the normalized (median of ratios) gene counts using DESeq2 (180). We then tested for significant differences in gene expression between autumn and spring using a Wald test, followed by multiple testing correction on resulting p-values with the Benjamini and Hochberg procedure (181). Significantly differentially expressed genes were binned by fold change (>1.58 log-fold change), to quantify differentially expressed genes of high effect. Significant differentially expressed genes were also used to identify KEGG GO enrichment (178), and compared against the candidate list of genes with branch shift changes in expression to identify genes both plastic across seasons and consistent with expression adaptation in the shrew.
Cell culture
Mustelo putoris furo (MPF) brain cells (ATCC, CRL-1656) were cultured in 10mL of Basal Medium Eagle (BME) supplemented with 14% sheep serum. Cells were cultured in 60.8-cm2 treated tissue culture dishes at 37°C in 5% CO2 atmosphere and 95% humidity. With induced cell death, cells were seeded into treated flat bottom six-well plates with 2mL/well of complete culture media.
Cell death induction
Three experiments were conducted to induce cell death in MPF cells (cold induced, peroxide induced, and heat induced), with only heat causing significant decrease in cell viability. Cold induced: Two 60.8-cm2 tissue culture dishes of MPF cells were seeded into two six-well plates for 24 hours before reaching 80% confluency. In one experiment, “short cold” apoptosis was performed by immersing the six-well plate into a 0°C ice-water bath for 2 hours, followed by 3 hours of rewarming to 37°C in an incubator. Cells were then treated with crystal violet staining. In a similar experiment, “long cold” apoptosis was performed by immersing the six-well plate into a 0°C water bath for 4 hours, followed by 24 hours of rewarming to 37°C before crystal violet staining. Peroxide-induced cell death: Two 60.8-cm2 tissue culture dishes of MPF cells were seeded into two six-well plates for 24 hours before reaching 80% confluency. Hydrogen peroxide (H2O2) induced cell death was performed by treating the plate with 125µM H2O2 followed by 3 hours of incubation at 37°C. Cells were then treated with crystal violet staining. Heat-induced cell death: Two 60.8-cm2 tissue culture dishes of MPF cells were seeded into two six-well plates for 24 hours before reaching 80% confluency. Heat cell death was attempted by immersing the six-well plate into a 45°C water bath for 2 hours, followed by re-cooling to 37°C prior to crystal violet staining.
In vitro transcription and transfection of cells
The AmpliScribe T7 High Yield Transcription kit (Lucigen) was used for RNA in vitro transcription with a 2-hour incubation at 37°C using gblock sequences available in supplemental. Sequences were designed with a T7 promoter inserted at the beginning to facilitate transcription (Integrated DNA Technologies).
RNA transcripts were purified using the Qiagen miRNA kit and concentration determined using a Nanodrop2000 (Thermo Fisher Scientific). Two 60.8-cm2 tissue culture dishes of MPF cells were seeded into two six-well plates for 24 hours before reaching 80% confluency. 300ng of BCL2L1 RNA or a scrambled form of BCL2L1 RNA were transfected with Lipofectamine 3000 into each well for 24 hours according to manufacturer’s instructions (Thermo Fisher Scientific). In short, two tubes were used to make a master mix for each RNA sample. Tube 1 contained 125µL/well of Opti-mem (Thermo Fisher Scientific) with 7.5 µL/well of Lipo3000. Tube 2 contained 125µL/well of Opti-mem, 5µL/well of P3000 and 300ng/well of RNA. Tube 2 was added to tube 1 and incubated 15 minutes at room temperature. The solution was then gently mixed and dispersed evenly to cells.
Crystal violet cell viability assays
Crystal violet solution was made using 500mg crystal violet powder in 100mL of 50% methanol. Cells were cultured to 80% confluency in two six-well plates. Complete media was aspirated in wells and washed three times in 1mL of 1X DPBS before 500µL of crystal violet solution was added to each well. Plates were wrapped in foil and placed on a shaker for 30 minutes. After time elapsed, crystal violet was removed, and cells were washed with tap water until free color was no longer visible. Plates were left at room temperature for 10-15 minutes until dry. 500µL of 100% methanol was added to each well and plates were put on a shaker for one hour at room temperature. 100µL of solution was taken from each well in triplicate and added to a 96-well plate. 100µL of 100% methanol was added in triplicate to a row for normalization of background absorbance. Wells were read at 570nm using a spectrophotometer (BioTek). For viability, data are reported as a percentage of control or optical density for absorbance. For statistical analysis, differences were judged to be statistically significant when p<0.05 by a Welch’s two sample t-test.
Data Availability
Supplementary tables, data, results, and code deposited and found on Github https://github.com/wrthomas315/Sorex_Hypothalamus_Transcriptomics. Raw sequencing data located in the NCBI Sequencing Read Archive (BioProject PRJNA941271). Supplemental Figures. S1–S2: DOI: https://doi.org/10.6084/m9.figshare.26049739.v1.
Additional information
Grants
Human Frontiers of Science Program, award: RGP0013/2019 (to Dina K. N. Dechmann, John Nieland, Liliana M. Dávalos).
WRT was supported in part by a Stony Brook University Presidential Innovation and Excellence award to LMD.
Disclosures
Authors declare no competing interests.
Author Contributions
Conceptualization: DD RGH WRT Methodology: ETO LMD RGH TAR WRT Software: WRT
Validation: APC ETO
Formal analysis: APC ETO TAR WRT Investigation: CB ETO TAR WRT Resources: CB DD ETO RGH TAR Data Curation: WRT
Writing – Original Draft: ETO TAR WRT
Writing – Review and Editing: APC CB DD DE ETO JDN LMD TAR Visualization: ETO TAR WRT
Supervision: LMD RGH
Project administration: DD LMD Funding acquisition: APC DD LMD JDN
References
- 1.Maturation of the adolescent brainNeuropsychiatr. Dis. Treat 9:449–461
- 2.Allometry and size in ontogeny and phylogenyBiol. Rev. Camb. Philos. Soc. 41:587–640
- 3.Seasonal reversible size changes in the braincase and mass of common shrews are flexibly modified by environmental conditionsSci. Rep 9:1–10
- 4.Seasonal and age change in shrews as an adaptive processSymp. zool. Soc. Lond 26:189–207
- 5.Water contents and seasonal changes of the brain-weight in shrewsActa Theriol. (Warsz 10:353–367
- 6.Seasonal and age changes in the weight of internal organs of shrewsActa Theriol. (Warsz 10:369–438
- 7.Stable carbon isotopes in breath reveal fast metabolic incorporation rates and seasonally variable but rapid fat turnover in the common shrew (Sorex araneus)J. Exp. Biol 220:2834–2841
- 8.Food resources and foraging habits of the common shrew, Sorex araneus: Does winter food shortage explain Dehnel’s phenomenon?Oikos 121:1593–1602
- 9.Winter reduction in body mass in a very small, nonhibernating mammal: Consequences for heat loss and metabolic ratesPhysiol. Biochem. Zool 86:9–18
- 10.Wintering strategy of voles and shrews in FindlandWinter Ecol. Small Mamm 10:139–148
- 11.Metabolic rate in common shrews is unaffected by temperature, leading to lower energetic costs through seasonal size reductionR. Soc. Open Sci 7https://doi.org/10.1098/rsos.191989
- 12.Molecular mechanisms of seasonal brain shrinkage and regrowth in Sorex araneusbioRxiv :1–16
- 13.Decimated little brown bats show potential for adaptive changeSci. Rep 10:1–10
- 14.Ontogeny and phylogeny of endothermy and torpor in mammals and birdsComp. Biochem. Physiol. - A Mol. Integr. Physiol 150:176–180
- 15.Parallel Accelerated Evolution in Distant Hibernators Reveals Candidate Cis Elements and Genetic Circuits Regulating Mammalian ObesityCell Rep 29:2608–2620
- 16.Hibernation is associated with increased survival and the evolution of slow life histories among mammalsProc. R. Soc. B Biol. Sci 278:3355–3363
- 17.Gene expression profiling in the hibernating primate, Cheirogaleus mediusGenome Biol. Evol 8:2413–2426
- 18.Comparative genomics of mammalian hibernators using gene networksIntegr. Comp. Biol 54:452–462
- 19.Comparative analyses of basal rate of metabolism in mammals : data selection does matterBiol. Rev 93:404–438
- 20.Evolution of Energetic Strategies in ShrewsEvolution of Shrews :309–346
- 21.Shrews, Chromosomes and SpeciationCambridge: Cambridge University Press
- 22.Ecology and mode-of-life explain lifespan variation in birds and mammalsProc. R. Soc. B 281https://doi.org/10.1098/rspb.2014.0298
- 23.Sugar for the brain: The role of glucose in physiological and pathological brain functionTrends Neurosci 36:587–597
- 24.Seasonality of Leptin Levels in the BAT of the Common Shrew (Sorex araneus)Z. Naturforsch 55:455–460
- 25.Neurodegeneration in juvenile Iberian pigs with diet-induced nonalcoholic fatty liver diseaseAm. J. Physiol. - Endocrinol. Metab 319:E592–E606
- 26.Impaired hepatic function and central dopaminergic denervation in a rodent model of Parkinson’s disease: A self-perpetuating crosstalk?Biochim. Biophys. Acta - Mol. Basis Dis 1822:176–184
- 27.Targeting Alzheimer’s Disease: The Critical Crosstalk between the Liver and BrainNutrients 14https://doi.org/10.3390/nu14204298
- 28.Human physiomimetic model integrating microphysiological systems of the gut, liver, and brain for studies of neurodegenerative diseasesSci. Adv 7https://doi.org/10.1126/SCIADV.ABD1707
- 29.Genotype-environment interaction and the evolution of phenotypic plasticityEvolution 39:505–522
- 30.Comparative ecological transcriptomics and the contribution of gene expression to the evolutionary potential of a threatened fishMol. Ecol 26:6841–6856
- 31.Species-specific molecular responses of wild coral reef fishes during a marine heatwaveSci. Adv 6:1–12
- 32.The evolution of gene expression underlying vision loss in cave animalsMol. Biol. Evol 35:2005–2014
- 33.Comparative regulomics supports pervasive selection on gene dosage following whole genome duplicationGenome Biol 22:1–18
- 34.The lasting after-effects of an ancient polyploidy on the genomes of teleostsPLoS One 15:1–26
- 35.The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisonsNat. Genet 48:427–437
- 36.The Atlantic salmon genome provides insights into rediploidizationNature 533:200–205
- 37.Subfunctionalization versus neofunctionalization after whole-genome duplicationNat. Genet 50:908–909
- 38.Sex-Specific Selection Drives the Evolution of Alternative Splicing in BirdsMol. Biol. Evol 38:519–530
- 39.Phylogenetic ANOVA: The expression variance and evolution model for quantitative trait evolutionSyst. Biol 64:695–708
- 40.A quantitative model for characterizing the evolutionary history of mammalian gene expressionbioRxiv :53–63
- 41.Convergen evolution of primate testis transcriptomes reflects mating strategybioRxiv :1–25
- 42.An atlas of the protein-coding genes in the human, pig, and mouse brainScience 367:1–16
- 43.Profound seasonal changes in brain size and architecture in the common shrewBrain Struct. Funct 223:2823–2840
- 44.Adaptations to a hypoxic lifestyle in naked mole-ratsJ. Exp. Biol 225https://doi.org/10.1242/jeb.196725
- 45.Hypothalamic control of energy balance: Insights into the role of synaptic plasticityTrends Neurosci 36:65–73
- 46.Sur la fine structure du lobe optique des oiseaux et sur l’origine réelle des nerfs optiquesJ. Int. Anat. Physiol 7:1–30
- 47.Postnatal development of the brain of the ferret Mustela putorius F Furo (Mustelidae Mammalia)Int. J. Mamm. Biol 44:127–131
- 48.Endosomal PI(3)P regulation by the COMMD/CCDC22/CCDC93 (CCC) complex controls membrane protein recyclingNat. Commun 10https://doi.org/10.1038/s41467-019-12221-6
- 49.FAM57B is a modulator of ceramide synthesis that regulates sphingolipid homeostasis and synaptic composition in the developing brainbioRxiv :1–23
- 50.The orphan G protein-coupled receptor 3 modulates amyloid-beta peptide generation in neuronsScience 323:946–951
- 51.Mice lacking GPR3 receptors display late-onset obese phenotype due to impaired thermogenic function in brown adipose tissueSci. Rep 5:1–5
- 52.Lmx1a Encodes a Rostral Set of Mesodiencephalic Dopaminergic Neurons Marked by the Wnt/B-Catenin Signaling Activator R-spondin 2PLoS One 8:1–12
- 53.Golgi-localized PAQR4 mediates antiapoptotic ceramidase activity in breast cancerCancer Res 80:2163–2174
- 54.Relaxin-3 and its role in neuroendocrine functionAnnals of the New York Academy of Sciences 1160:250–255
- 55.Relaxin-3 systems in the brain — The first 10 yearsJ. Chem. Neuroanat. jo 42:262–275
- 56.Centrally administered relaxin-3 induces Fos expression in the osmosensitive areas in rat brain and facilitates water intakePeptides 31:1124–1130
- 57.Central relaxin-3 administration causes hyperphagia in male wistar ratsEndocrinology 146:3295–3300
- 58.Sex differences in the effects of chronic stress and food restriction on body weight gain and brain expression of CRF and relaxin-3 in rats. GenesBrain Behav 12:370–387
- 59.Differential effects of relaxin-3 and a selective relaxin-3 receptor agonist on food and water intake and hypothalamic neuronal activity in ratsBehav. Brain Res 336:135–144
- 60.Differential effects of central administration of relaxin-3 on food intake and hypothalamic neuropeptides in male and female rats. GenesBrain Behav 14:550–563
- 61.Sex-specific effects of relaxin-3 on food intake and brain expression of corticotropin-releasing factor in ratsEndocrinology 156:523–533
- 62.Hypothalamic mapping of orexigenic action and Fos-like immunoreactivity following relaxin-3 administration in male Wistar ratsAm. J. Physiol. - Endocrinol. Metab 292:913–919
- 63.Relaxin induces matrix-metalloproteinases-9 and -13 via RXFP1: Induction of MMP-9 involves the PI3K, ERK, Akt and PKC-ζ pathwaysMol. Cell. Endocrinol 363:46–61
- 64.Relaxin stimulates multiple signaling pathways: Activation of cAMP, PI3K, and PKCζ in THP-1 cellsAnn. N. Y. Acad. Sci. 1041:272–279
- 65.Behavior Apolipoprotein E reduces food intake via PI3K / Akt signaling pathway in the hypothalamusPhysiol. Behav 105:124–128
- 66.Biochemical and Biophysical Research Communications Apolipoprotein A-IV exerts its anorectic action through a PI3K / Akt signaling pathway in the hypothalamusBiochem. Biophys. Res. Commun 494:152–157
- 67.Intracellular signalling: Key enzyme in leptin-induced anorexiaNature 413:794–795
- 68.Insulin Action in AgRP-Expressing Neurons Is Required for Suppression of Hepatic Glucose ProductionCell Metab 5:438–449
- 69.Phosphatidyl inositol 3-kinase signaling in hypothalamic proopiomelanocortin neurons contributes to the regulation of glucose homeostasisEndocrinology 150:4874–4882
- 70.PI3K/AKT/mTOR Pathway in AngiogenesisFront. Mol. Neurosci 4:1–8
- 71.Epidermal growth factor receptor transcriptionally up-regulates vascular endothelial growth factor expression in human glioblastoma cells via a pathway involving phosphatidylinositol 3’-kinase and distinct from that induced by hypoxiaCancer Res 60:5879–5886
- 72.Critical nodes in signalling pathways: insights into insulin actionNat. Rev. Mol. Cell Biol 7:85–96
- 73.Comparative transcriptome analysis of hypothalamus-regulated feed intake induced by exogenous visfatin in chicksBMC Genomics 19:1–17
- 74.RNA-seq and differential expression analysis of the duck transcriptome : The effect of short-term cage-rearingbioRxiv https://doi.org/10.1101/2021.05.13.444049
- 75.Identification of the acclimation genes in transcriptomic responses to heat stress of White Pekin duckCell Stress Chaperones 22:787–797
- 76.Identification of differentially expressed genes in hypothalamus of chicken during cold stressMiol Biol Rep 41:2243–2248
- 77.FoxO1 negatively regulates leptin- induced POMC transcription through its direct interaction with STAT3Biochem. J 466:291–298
- 78.The role of leptin in the control of insulin-glucose axisFront. Neurosci 7:1–12
- 79.Induction of hypothalamic Sirt1 leads to cessation of feeding via agouti- related peptideEndocrinology 151:2556–2566
- 80.The LXXLL motif of murine forkhead transcription factor FoxO1 mediates Sirt1-dependent transcriptional activityJ. Clin. Invest 116:2473–2483
- 81.Overexpression of the transcription factor foxo4 is associated with rapid glucose clearanceMol. Cell. Endocrinol 307:217–223
- 82.Seasonal leptin resistance is associated with impaired signalling via JAK2-STAT3 but not ERK, possibly mediated by reduced hypothalamic GRB2 proteinJ. Comp. Physiol. B Biochem. Syst. Environ. Physiol 182:553–567
- 83.Deficiency of Irx5 protects mice from obesity and associated metabolic abnormalitiesInt. J. Obes 46:2029–2039
- 84.Irx3 and Irx5 in Ins2-Cre + cells regulate hypothalamic postnatal neurogenesis and leptin responseNat. Metab 3:701–713
- 85.Impaired brain angiogenesis and neuronal apoptosis induced by conditional homozygous inactivation of vascular endothelial growth factorThromb. Haemost 91:595–605
- 86.Abnormal blood vessel development and lethality in embryos lacking a single vascular endothelial growth factor alleleNature 380:435–439
- 87.Wnt/β-catenin signaling is required for CNS, but not non-CNS, angiogenesisProc. Natl. Acad. Sci. U. S. A 106:641–646
- 88.Astrocyte-derived VEGF-A drives blood-brain barrier disruption in CNS inflammatory diseaseJ. Clin. Invest 122:2454–2468
- 89.Dexamethasone coordinately regulates angiopoietin-1 and VEGF: A mechanism of glucocorticoid-induced stabilization of blood-brain barrierBiochem. Biophys. Res. Commun 372:243–248
- 90.VEGF-mediated disruption of endothelial CLN-5 promotes blood-brain barrier breakdownProc. Natl. Acad. Sci. U. S. A 106:1977–1982
- 91.Astrocytic VEGFA: An essential mediator in blood–brain-barrier disruption in Parkinson’s diseaseGlia 70:337–353
- 92.Development, maintenance and disruption of the blood-brain barrierNat. Med 19:1584–1596
- 93.Circumventricular organs: A novel site of neural stem cells in the adult brainMol. Cell. Neurosci 41:337–347
- 94.Tanycytic VEGF-A boosts blood-hypothalamus barrier plasticity and access of metabolic signals to the arcuate nucleus in response to fastingCell Metab 17:607–617
- 95.Development and function of the blood-brain barrier in the context of metabolic controlFront. Neurosci 11:1–12
- 96.Hypothalamic sensing of fatty acidsNat. Neurosci 8:579–584
- 97.Seasonal vascular plasticity in the mediobasal hypothalamus of the adult eweHistochem. Cell Biol 157:581–593
- 98.Mechanisms regulating angiogenesis underlie seasonal control of pituitary functionProc. Natl. Acad. Sci. U. S. A 114:E2514–E2523
- 99.Local regulators of seasonal reproduction processes in uterus masculinus of an adult male european bison (bison bonasus, linnaeus 1758)J. Physiol. Pharmacol 69:747–753
- 100.Seasonal expressions of VEGF and its receptors VEGFR1 and VEGFR2 in the prostate of the wild ground squirrels (Spermophilus dauricus)Eur. J. Histochem 65
- 101.Whole-brain functional imaging at cellular resolution using light-sheet microscopyNat. Methods 10:413–420
- 102.Whole-Brain Calcium Imaging Reveals an Intrinsic Functional Network in DrosophilaCurr. Biol 27:2389–2396
- 103.Calcium regulation of the neuronal growth coneTrends Neurosci 11:315–321
- 104.Neurotransmitters in the regulation of neuronal cytoarchitectureBrain Res. Rev 472:179–212
- 105.Signaling from synapse to nucleus: Postsynaptic CREB phosphorylation during multiple forms of hippocampal synaptic plasticityNeuron 16:89–101
- 106.Ca2+-dependent regulation in neuronal gene expressionCurr. Opin. Neurobiol 7:419–429
- 107.Neuronal calciumNeuron. 21:13–26
- 108.Nuclear calcium signaling controls expression of a large gene pool: Identification of a gene program for acquired neuroprotection induced by synaptic activityPLoS Genet 5https://doi.org/10.1371/journal.pgen.1000604
- 109.Regulation of dendritic development by calcium signalingCell Calcium 37:411–416
- 110.Regulation of dendritic growth and plasticity by local and global calcium dynamicsCell Calcium 37:403–409
- 111.Synaptic components necessary for retrograde signaling triggered by calcium/calmodulin-dependent protein kinase II during synaptogenesisNeuroscience 145:1007–1015
- 112.Spatial changes in calcium signaling during the establishment of neuronal polarity and synaptogenesisJ. Cell Biol 126:1527–1536
- 113.Calcium dynamics at developing synapses: Mechanisms and functionsEur. J. Neurosci 32:218–223
- 114.Postsynaptic Calcium Is Sufficient for Potentiation of Hippocampal Synaptic TransmissionScience 242:81–84
- 115.Role of Ca2+-stimulated adenylyl cyclases in LTP and memory formationInt. J. Dev. Neurosci 19:387–394
- 116.Nuclear calcium signalling in the regulation of brain functionNat. Rev. Neurosci 14:593–608
- 117.Transcription-dependent neuronal plasticityEur. J. Biochem. 267:5280–5283
- 118.Local Calcium Signaling in NeuronsNeuron 40:331–346
- 119.Spatial features of calcium-regulated gene expressionBioEssays 19:657–660
- 120.Histamine H2-receptors on single central neuronesNature 255:634–635
- 121.Histamine and noradrenaline decrease calcium-activated potassium conductance in hippocampal pyramidal cellsNature 302:432–434
- 122.Histamine in the nervous systemPhysiol. Rev 88:1183–1241
- 123.Inhibitory and excitatory effects of histamine on suprachiasmatic neurons in rat hypothalamic slice preparationNeurosci. Lett 41:109–113
- 124.Effects of intracerebroventricularly infused histamine and selective H1, H2 and H3 agonists on food and water intake and urine flow in Wistar ratsBrain Res. 793:279–288
- 125.Involvement of hypothalamic histamine H1 receptor in the regulation of feeding rhythm and obesityDiabetes 53:2250–2260
- 126.Modulation of behavior by the histaminergic system: Lessons from H1R-and H2R-deficient miceNeurosci. Biobehav. Rev 42:252–266
- 127.Ablation in Mice of the mTORC Components raptor, rictor, or mLST8 Reveals that mTORC2 Is Required for Signaling to Akt-FOXO and PKCα DevCell 11:859–871
- 128.mTOR Signaling in Growth, Metabolism, and DiseaseCell 168:960–976
- 129.Loss of cyclophilin D reveals a critical role for mitochondrial permeability transition in cell deathNature 434:984–984
- 130.Cyclophilin D- dependentmitochondrial permeabilitytransition regulates some necroticbut not apoptotic cell deathNature 434:648–652
- 131.Role of cyclophilin D-dependent mitochondrial permeability transition in glutamate-induced calcium deregulation and excitotoxic neuronal deathExp. Neurol 218:171–182
- 132.A pore way to dieNature 430:984–984
- 133.SphK1 and SphK2, sphingosine kinase isoenzymes with opposing functions in sphingolipid metabolismJ. Biol. Chem 280:37118–37129
- 134.Sphingosine kinase type 2 is a putative BH3-only protein that induces apoptosisJ. Biol. Chem 278:40330–40336
- 135.Inhibiting sphingosine kinase 2 mitigates mutant Huntingtin-induced neurodegeneration in neuron models of Huntington diseaseHum. Mol. Genet 26:1305–1317
- 136.BACE1 activity is modulated by cell-associated sphingosine-1-phosphateJ. Neurosci 31:6850–6857
- 137.Cyclophilin D deficiency attenuates mitochondrial and neuronal perturbation and ameliorates learning and memory in Alzheimer’s diseaseNat. Med 14:1097–1105
- 138.Cyclophilin D deficiency improves mitochondrial function and learning/memory in aging Alzheimer disease mouse modelNeurobiol. Aging 32:398–406
- 139.ThesisMax Planck Institute for Ornithology
- 140.Modulation of apoptotic pathways in intestinal mucosa during hibernationAm. J. Physiol. - Regul. Integr. Comp. Physiol 289:586–595
- 141.Response of the JAK-STAT pathway to mammalian hibernation in 13-lined ground squirrel striated muscleMol. Cell. Biochem 414:115–127
- 142.Anti-apoptotic signaling as a cytoprotective mechanism in mammalian hibernationPeerJ :1–21
- 143.Neuron-type-specific signaling by the p75NTR death receptor is regulated by differential proteolytic cleavageJ. Cell Sci 128:1507–1517
- 144.A prosurvival function for the p75 receptor death domain mediated via the caspase recruitment domain receptor- interacting protein 2J. Neurosci 21:5854–5863
- 145.Selective Activation of NF- B by Nerve Growth Factor Through the Neurotrophin Receptor p75Science 272:542–545
- 146.Constitutive nuclear factor- κB activity is required for central neuron survivalJ. Neurosci 22:8466–8475
- 147.Activation of nuclear factor κB and bcl-x survival gene expression by nerve growth factor requires tyrosine phosphorylation of IκBαJ. Cell Biol 152:753–763
- 148.Massive Cell Death of Immature Hematopoietic Cells and Neurons in Bcl-x-Deficient MiceScience 267:1506–1510
- 149.Bcl-XL can inhibit apoptosis in cells that have undergone Fas- induced protease activationProc. Natl. Acad. Sci. U. S. A 94:3759–3764
- 150.COMMD1 is linked to the WASH complex and regulates endosomal trafficking of the copper transporter ATP7AMol. Biol. Cell 26:91–103
- 151.CCC- and WASH-mediated endosomal sorting of LDLR is required for normal clearance of circulating LDLNat. Commun 7:1–11
- 152.The roles of intracellular protein-degradation pathways in neurodegenerationNature 443:780–786
- 153.Explosion in the complexity of membrane protein recyclingAm. J. Physiol. - Cell Physiol 320:C483–C494
- 154.Retromer stabilization results in neuroprotection in a model of Amyotrophic Lateral SclerosisNat. Commun 11:1–17
- 155.Model-guided microarray implicates the retromer complex in Alzheimer’s diseaseAnn. Neurol 58:909–919
- 156.The neuronal sortilin-related receptor SORL1 is genetically associated with Alzheimer diseaseNat. Genet 39:168–177
- 157.Diabetes-associated SorCS1 regulates Alzheimer’s amyloid-β metabolism: Evidence for involvement of SorL1 and the retromer complexJ. Neurosci 30:13110–13115
- 158.VPS35 Deficiency or Mutation Causes Dopaminergic Neuronal Loss by Impairing Mitochondrial Fusion and FunctionCell Rep 12:1631–1643
- 159.CCDC22: A novel candidate gene for syndromic X-linked intellectual disabilityMol. Psychiatry 17:4–7
- 160.Missense variant in CCDC22 causes X-linked recessive intellectual disability with features of Ritscher-Schinzel/3C syndromeEur. J. Hum. Genet 23:633–638
- 161.Expanding the pre- and postnatal phenotype of WASHC5 and CCDC22 -related Ritscher-Schinzel syndromesEur. J. Med. Genet 65:1–7
- 162.Expansion of the CCDC22 associated ritscher-schinzel/3C syndrome and review of the literature: Should the minimal diagnostic criteria be revised?Eur. J. Med. Genet 64https://doi.org/10.1016/j.ejmg.2021.104246
- 163.Behavioral abnormalities and circuit defects in the basal ganglia of a mouse model of 16p11.2 deletion syndromeCell Rep 7:1077–1092
- 164.The 16p11.2 homologs fam57ba and doc2a generate certain brain and body phenotypesHum. Mol. Genet 26:3699–3712
- 165.Association between Microdeletion and Microduplication at 16p11.2 and AutismN. Engl. J. Med 358:667–675
- 166.Reciprocal Effects on Neurocognitive and Metabolic Phenotypes in Mouse Models of 16p11.2 Deletion and Duplication SyndromesPLoS Genet 12:1–35
- 167.Oligogenic Effects of 16p11.2 Copy-Number Variation on Craniofacial DevelopmentCell Rep 28:3320–3328
- 168.Zebrafish homologs of genes within 16p11.2, a genomic region associated with brain disorders, are active during brain development, and include two deletion dosage sensor genesDis. Model. Mech 5:834–851
- 169.Fam57b (Family with sequence similarity 57, member B), a novel peroxisome proliferator-activated receptor y target gene that regulates adipogenesis through ceramide synthesisJ. Biol. Chem 288:4522–4537
- 170.16pdel lipid changes in iPSC-derived neurons and function of FAM57B in lipid metabolism and synaptogenesisiScience 25
- 171.Developmental expression of GPR3 in rodent cerebellar granule neurons is associated with cell survival and protects neurons from various apoptotic stimuliNeurobiol. Dis 68:215–227
- 172.Loss of GPR3 reduces the amyloid plaque burden and improves memory in Alzheimer’s disease mouse modelsSci. Transl. Med 7https://doi.org/10.1126/scitranslmed.aab3492
- 173.Fastp: An ultra-fast all-in-one FASTQ preprocessorBioinformatics 34:i884–i890
- 174.Near-optimal probabilistic RNA-seq quantificationNat. Biotechnol 34:525–527
- 175.OrthoFinder: Phylogenetic orthology inference for comparative genomicsGenome Biol 20:1–14
- 176.Pairwise comparisons across species are problematic when analyzing functional genomic dataProc. Natl. Acad. Sci. U. S. A 115:E409–E417
- 177.A species-level timeline of mammal evolution integrating phylogenomic dataNature 602:263–267
- 178.Systematic and integrative analysis of large gene lists using DAVID bioinformatics resourcesNat. Protoc 4:44–57
- 179.TCseq: Time course sequencing data analysisR package
- 180.Moderated estimation of fold change and dispersion for RNA- seq data with DESeq2Genome Biol 15:1–21
- 181.Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple TestingJ. R. Stat. Soc 57:289–300
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Copyright
© 2024, Thomas et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 263
- downloads
- 0
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.