The evolution of thermogenesis in mammals.

Comparative genomics elucidates the steps enabling heat production in fat tissue.

Dynamic RNA Regulation in the Brain Underlies Physiological Plasticity in a Hibernating Mammal.

Hibernation is a physiological and behavioral phenotype that minimizes energy expenditure. Hibernators cycle between profound depression and rapid hyperactivation of multiple physiological processes, challenging our concept of mammalian homeostasis. How the hibernator orchestrates and survives these extremes while maintaining cell to organismal viability is unknown. Here, we enhance the genome integrity and annotation of a model hibernator, the 13-lined ground squirrel. Our new assembly brings this genome to near chromosome-level contiguity and adds thousands of previously unannotated genes. These new genomic resources were used to identify 6,505 hibernation-related, differentially-expressed and processed transcripts using RNA-seq data from three brain regions in animals whose physiological status was precisely defined using body temperature telemetry. A software tool, squirrelBox, was developed to foster further data analyses and visualization. SquirrelBox includes a comprehensive toolset for rapid visualization of gene level and cluster group dynamics, sequence scanning of k-mer and domains, and interactive exploration of gene lists. Using these new tools and data, we deconvolute seasonal from temperature-dependent effects on the brain transcriptome during hibernation for the first time, highlighting the importance of carefully timed samples for studies of differential gene expression in hibernation. The identified genes include a regulatory network of RNA binding proteins that are dynamic in hibernation along with the composition of the RNA pool. In addition to passive effects of temperature, we provide evidence for regulated transcription and RNA turnover during hibernation. Significant alternative splicing, largely temperature dependent, also occurs during hibernation. These findings form a crucial first step and provide a roadmap for future work toward defining novel mechanisms of tissue protection and metabolic depression that may 1 day be applied toward improving human health.

Genetic variation drives seasonal onset of hibernation in the 13-lined ground squirrel.

Hibernation in sciurid rodents is a dynamic phenotype timed by a circannual clock. When housed in an animal facility, 13-lined ground squirrels exhibit variation in seasonal onset of hibernation, which is not explained by environmental or biological factors. We hypothesized that genetic factors instead drive variation in timing. After increasing genome contiguity, here, we employ a genotype-by-sequencing approach to characterize genetic variation in 153 ground squirrels. Combined with datalogger records (n = 72), we estimate high heritability (61-100%) for hibernation onset. Applying a genome-wide scan with 46,996 variants, we identify 2 loci significantly (p < 7.14 × 10-6), and 12 loci suggestively (p < 2.13 × 10-4), associated with onset. At the most significant locus, whole-genome resequencing reveals a putative causal variant in the promoter of FAM204A. Expression quantitative trait loci (eQTL) analyses further reveal gene associations for 8/14 loci. Our results highlight the power of applying genetic mapping to hibernation and present new insight into genetics driving its onset.

Comparative tissue transcriptomics highlights dynamic differences among tissues but conserved metabolic transcript prioritization in preparation for arousal from torpor.

During the hibernation season, 13-lined ground squirrels spend days to weeks in torpor with body temperatures near freezing then spontaneously rewarm. The molecular drivers of the drastic physiological changes that orchestrate and permit torpor are not well understood. Although transcription effectively ceases at the low body temperatures of torpor, previous work has demonstrated that some transcripts are protected from bulk degradation in brown adipose tissue (BAT), consistent with the importance of their protein products for metabolic heat generation during arousal from torpor. We examined the transcriptome of skeletal muscle, heart, and liver to determine the patterns of differentially expressed genes in these tissues, and whether, like BAT, a subset of these were relatively increased during torpor. EDGE-tags were quantified from five distinct physiological states representing the seasonal and torpor-arousal cycles of 13-lined ground squirrels. Supervised clustering on relative transcript abundances with Random Forest separated the two states bracketing prolonged torpor, entrance into and aroused from torpor, in all three tissues. Independent analyses identified 3347, 6784, and 2433 differentially expressed transcripts among all sampling points in heart, skeletal muscle, and liver, respectively. There were few differentially expressed genes in common across all three tissues; these were enriched in mitochondrial and apoptotic pathway components. Divisive clustering of these data revealed unique cohorts of transcripts that increased across the torpor bout in each tissue with patterns reflecting various combinations of cycling within and between seasons as well as between torpor and arousal. Transcripts that increased across the torpor bout were likewise tissue specific. These data shed new light on the biochemical pathways that alter in concert with hibernation phenotype and provide a rich resource for further hypothesis-based studies.

Proteomics approaches shed new light on hibernation physiology.

The broad phylogenetic distribution and rapid phenotypic transitions of mammalian hibernators imply that hibernation is accomplished by differential expression of common genes. Traditional candidate gene approaches have thus far explained little of the molecular mechanisms underlying hibernation, likely due to (1) incomplete and imprecise sampling of a complex phenotype, and (2) the forming of hypotheses about which genes might be important based on studies of model organisms incapable of such dynamic physiology. Unbiased screening approaches, such as proteomics, offer an alternative means to discover the cellular underpinnings that permit successful hibernation and may reveal previously overlooked, important pathways. Here, we review the findings that have emerged from proteomics studies of hibernation. One striking feature is the stability of the proteome, especially across the extreme physiological shifts of torpor-arousal cycles during hibernation. This has led to subsequent investigations of the role of post-translational protein modifications in altering protein activity without energetically wasteful removal and rebuilding of protein pools. Another unexpected finding is the paucity of universal proteomic adjustments across organ systems in response to the extreme metabolic fluctuations despite the universality of their physiological challenges; rather each organ appears to respond in a unique, tissue-specific manner. Additional research is needed to extend and synthesize these results before it will be possible to address the whole body physiology of hibernation.

Enhanced stability and polyadenylation of select mRNAs support rapid thermogenesis in the brown fat of a hibernator.

During hibernation, animals cycle between torpor and arousal. These cycles involve dramatic but poorly understood mechanisms of dynamic physiological regulation at the level of gene expression. Each cycle, Brown Adipose Tissue (BAT) drives periodic arousal from torpor by generating essential heat. We applied digital transcriptome analysis to precisely timed samples to identify molecular pathways that underlie the intense activity cycles of hibernator BAT. A cohort of transcripts increased during torpor, paradoxical because transcription effectively ceases at these low temperatures. We show that this increase occurs not by elevated transcription but rather by enhanced stabilization associated with maintenance and/or extension of long poly(A) tails. Mathematical modeling further supports a temperature-sensitive mechanism to protect a subset of transcripts from ongoing bulk degradation instead of increased transcription. This subset was enriched in a C-rich motif and genes required for BAT activation, suggesting a model and mechanism to prioritize translation of key proteins for thermogenesis.

Metabolic changes associated with the long winter fast dominate the liver proteome in 13-lined ground squirrels.

Small-bodied hibernators partition the year between active homeothermy and hibernating heterothermy accompanied by fasting. To define molecular events underlying hibernation that are both dependent and independent of fasting, we analyzed the liver proteome among two active and four hibernation states in 13-lined ground squirrels. We also examined fall animals transitioning between fed homeothermy and fasting heterothermy. Significantly enriched pathways differing between activity and hibernation were biased toward metabolic enzymes, concordant with the fuel shifts accompanying fasting physiology. Although metabolic reprogramming to support fasting dominated these data, arousing (rewarming) animals had the most distinct proteome among the hibernation states. Instead of a dominant metabolic enzyme signature, torpor-arousal cycles featured differences in plasma proteins and intracellular membrane traffic and its regulation. Phosphorylated NSFL1C, a membrane regulator, exhibited this torpor-arousal cycle pattern; its role in autophagosome formation may promote utilization of local substrates upon metabolic reactivation in arousal. Fall animals transitioning to hibernation lagged in their proteomic adjustment, indicating that the liver is more responsive than preparatory to the metabolic reprogramming of hibernation. Specifically, torpor use had little impact on the fall liver proteome, consistent with a dominant role of nutritional status. In contrast to our prediction of reprogramming the transition between activity and hibernation by gene expression and then within-hibernation transitions by posttranslational modification (PTM), we found extremely limited evidence of reversible PTMs within torpor-arousal cycles. Rather, acetylation contributed to seasonal differences, being highest in winter (specifically in torpor), consistent with fasting physiology and decreased abundance of the mitochondrial deacetylase, SIRT3.

Daily Rhythms of PERIOD protein in the eyestalk of the American lobster, Homarus americanus.

The daily rhythm of PERIOD protein (PER) expression is an integral component of the circadian clock, which is found among a broad range of animal species including fruit flies, marine mollusks and even humans. The use of antibodies directed against PER has provided a helpful tool in the discovery of PER homologues and the labeling of putative pacemaker cells, especially in animals for which an annotated genome is not readily available. In this study, DrosophilaPER antibodies were used to probe for PER in the American lobster, Homarus americanus. This species exhibits robust endogenous circadian rhythms but the circadian clock has yet to be located or characterized. PER was detected in the eyestalks of the lobster but not in the brain. Furthermore, a significant effect of the LD cycle on daily PER abundance was identified, and PER was significantly more abundant at mid dark than in early light or mid light hours. Our results suggest that PER is a part of the molecular machinery of the circadian clock located in the eyestalk of the lobster.

Multistate proteomics analysis reveals novel strategies used by a hibernator to precondition the heart and conserve ATP for winter heterothermy.

The hibernator's heart functions continuously and avoids damage across the wide temperature range of winter heterothermy. To define the molecular basis of this phenotype, we quantified proteomic changes in the 13-lined ground squirrel heart among eight distinct physiological states encompassing the hibernator's year. Unsupervised clustering revealed a prominent seasonal separation between the summer homeotherms and winter heterotherms, whereas within-season state separation was limited. Further, animals torpid in the fall were intermediate to summer and winter, consistent with the transitional nature of this phase. A seasonal analysis revealed that the relative abundances of protein spots were mainly winter-increased. The winter-elevated proteins were involved in fatty acid catabolism and protein folding, whereas the winter-depleted proteins included those that degrade branched-chain amino acids. To identify further state-dependent changes, protein spots were re-evaluated with respect to specific physiological state, confirming the predominance of seasonal differences. Additionally, chaperone and heat shock proteins increased in winter, including HSPA4, HSPB6, and HSP90AB1, which have known roles in protecting against ischemia-reperfusion injury and apoptosis. The most significant and greatest fold change observed was a disappearance of phospho-cofilin 2 at low body temperature, likely a strategy to preserve ATP. The robust summer-to-winter seasonal proteomic shift implies that a winter-protected state is orchestrated before prolonged torpor ensues. Additionally, the general preservation of the proteome during winter hibernation and an increase of stress response proteins, together with dephosphorylation of cofilin 2, highlight the importance of ATP-conserving mechanisms for winter cardioprotection.

Replication and further characterization of a Type 1 diabetes-associated locus at the telomeric end of the major histocompatibility complex.

BACKGROUND: We recently reported an association between Type 1 diabetes and the telomeric major histocompatibility complex (MHC) single nucleotide polymorphism (SNP) rs1233478. As further families have been analyzed in the Type 1 Diabetes Genetics Consortium (T1DGC), we tested replication of the association and, with more data, analyzed haplotypic associations. METHODS: An additional 2717 case and 1315 control chromosomes have been analyzed from the T1DGC, with human leukocyte antigen (HLA) typing and data for 2837 SNPs across the MHC region. RESULTS: We confirmed the association of rs1233478 (new data only: P=2.2E-5, OR=1.4). We also found two additional SNPs nearby that were significantly associated with Type 1 diabetes (new data only rs3131020: P=8.3E-9, OR=0.65; rs1592410: P=2.2E-8, OR=1.5). For studies of Type 1 diabetes in the MHC region, it is critical to account for linkage disequilibrium with the HLA genes. Logistic regression analysis of these new data indicated that the effects of rs3131020 and rs1592410 on Type 1 diabetes risk are independent of HLA alleles (rs3131020: P=2.3E-3, OR=0.73; rs1592410: P=2.1E-3, OR=1.4). Haplotypes of 12 SNPs (including the three highly significant SNPs) stratify diabetes risk (high risk, protective, and neutral), with high-risk haplotypes limited to approximately 20,000 bp in length. The 20,000-bp region is telomeric of the UBD gene and contains LOC729653, a hypothetical gene. CONCLUSIONS: We believe that polymorphisms of the telomeric MHC locus LOC729653 may confer risk for Type 1 diabetes.