Hibernation in brown bears is an annual process involving multiple physiologically distinct seasons—hibernation, active, and hyperphagia. While recent studies have characterized broad patterns of differential gene regulation and isoform usage between hibernation and active seasons, patterns of gene and isoform expression during hyperphagia remain relatively poorly understood. The hyperphagia stage occurs between active and hibernation seasons and involves the accumulation of large fat reserves in preparation for hibernation. Here, we use time-series analyses of gene expression and isoform usage to interrogate transcriptomic regulation associated with all three seasons. We identify a large number of genes with significant differential isoform usage (DIU) across seasons and show that these patterns of isoform usage are largely tissue-specific. We also show that DIU and differential gene-level expression responses are generally non-overlapping, with only a small subset of multi-isoform genes showing evidence of both gene-level expression changes and changes in isoform usage across seasons. Additionally, we investigate nuanced regulation of candidate genes involved in the insulin signaling pathway and find evidence of hyperphagia-specific gene expression and isoform regulation that may enhance fat accumulation during hyperphagia. Our findings highlight the value of using temporal analyses of both gene- and isoform-level gene expression when interrogating complex physiological phenotypes and provide new insight into the mechanisms underlying seasonal changes in bear physiology.
This content will become publicly available on August 23, 2024
Feeding during hibernation shifts gene expression toward active season levels in brown bears (Ursus arctos)
Hibernation in bears involves a suite of metabolical and physiological changes, including the onset of insulin resistance, that are driven in part by sweeping changes in gene expression in multiple tissues. Feeding bears glucose during hibernation partially restores active season physiological phenotypes, including partial resensitization to insulin, but the molecular mechanisms underlying this transition remain poorly understood. Here, we analyze tissue-level gene expression in adipose, liver, and muscle to identify genes that respond to midhibernation glucose feeding and thus potentially drive postfeeding metabolical and physiological shifts. We show that midhibernation feeding stimulates differential expression in all analyzed tissues of hibernating bears and that a subset of these genes responds specifically by shifting expression toward levels typical of the active season. Inferences of upstream regulatory molecules potentially driving these postfeeding responses implicate peroxisome proliferator-activated receptor gamma (PPARG) and other known regulators of insulin sensitivity, providing new insight into high-level regulatory mechanisms involved in shifting metabolic phenotypes between hibernation and active states.
more »
« less
- Award ID(s):
- 2138649
- NSF-PAR ID:
- 10469444
- Publisher / Repository:
- American Physiological Society
- Date Published:
- Journal Name:
- Physiological genomics
- Volume:
- 55
- Issue:
- 9
- ISSN:
- 1531-2267
- Subject(s) / Keyword(s):
- ["differential isoform usage","gene expression","hibernation","insulin resistance","PPAR signaling"]
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Hibernation is a highly seasonal physiological adaptation that allows brown bears (Ursus arctos) to survive extended periods of low food availability. Similarly, daily or circadian rhythms conserve energy by coordinating body processes to optimally match the environmental light/dark cycle. Brown bears express circadian rhythms in vivo and their cells do invitro throughout the year, suggesting that these rhythms may play important roles during periods of negative energy balance. Here, we use time-series analysis of RNA sequencing data and timed measurements of ATP production in adipose-derived fibroblasts from active and hibernation seasons under two temperature conditions to confirm that rhythmicity was present. Culture temperature matching that of hibernation body temperature (34°C) resulted in a delay of daily peak ATP production in comparison with active season body temperatures (37°C). The timing of peaks of mitochondrial gene transcription was altered as were the amplitudes of transcripts coding for enzymes of the electron transport chain. Additionally, we observed changes in mean expression and timing of key metabolic genes such as SIRT1 and AMPK which are linked to the circadian system and energy balance. The amplitudes of several circadian gene transcripts were also reduced. These results reveal a link between energy conservation and a functioning circadian system in hibernationmore » « less
-
Objectives Complex physiological adaptations often involve the coordination of molecular responses across multiple tissues. Establishing transcriptomic resources for non-traditional model organisms with phenotypes of interest can provide a foundation for understanding the genomic basis of these phenotypes, and the degree to which these resemble, or contrast, those of traditional model organisms. Here, we present a one-of-a-kind gene expression dataset generated from multiple tissues of two hibernating brown bears (Ursus arctos). Data description This dataset is comprised of 26 samples collected from 13 tissues of two hibernating brown bears. These samples were collected opportunistically and are typically not possible to attain, resulting in a highly unique and valuable gene expression dataset. In combination with previously published datasets, this new transcriptomic resource will facilitate detailed investigation of hibernation physiology in bears, and the potential to translate aspects of this biology to treat human disease.more » « less
-
more » « less
Abstract Objectives Complex physiological adaptations often involve the coordination of molecular responses across multiple tissues. Establishing transcriptomic resources for non-traditional model organisms with phenotypes of interest can provide a foundation for understanding the genomic basis of these phenotypes, and the degree to which these resemble, or contrast, those of traditional model organisms. Here, we present a one-of-a-kind gene expression dataset generated from multiple tissues of two hibernating brown bears (
Ursus arctos ).Data description This dataset is comprised of 26 samples collected from 13 tissues of two hibernating brown bears. These samples were collected opportunistically and are typically not possible to attain, resulting in a highly unique and valuable gene expression dataset. In combination with previously published datasets, this new transcriptomic resource will facilitate detailed investigation of hibernation physiology in bears, and the potential to translate aspects of this biology to treat human disease.
-
INTRODUCTION Diverse phenotypes, including large brains relative to body size, group living, and vocal learning ability, have evolved multiple times throughout mammalian history. These shared phenotypes may have arisen repeatedly by means of common mechanisms discernible through genome comparisons. RATIONALE Protein-coding sequence differences have failed to fully explain the evolution of multiple mammalian phenotypes. This suggests that these phenotypes have evolved at least in part through changes in gene expression, meaning that their differences across species may be caused by differences in genome sequence at enhancer regions that control gene expression in specific tissues and cell types. Yet the enhancers involved in phenotype evolution are largely unknown. Sequence conservation–based approaches for identifying such enhancers are limited because enhancer activity can be conserved even when the individual nucleotides within the sequence are poorly conserved. This is due to an overwhelming number of cases where nucleotides turn over at a high rate, but a similar combination of transcription factor binding sites and other sequence features can be maintained across millions of years of evolution, allowing the function of the enhancer to be conserved in a particular cell type or tissue. Experimentally measuring the function of orthologous enhancers across dozens of species is currently infeasible, but new machine learning methods make it possible to make reliable sequence-based predictions of enhancer function across species in specific tissues and cell types. RESULTS To overcome the limits of studying individual nucleotides, we developed the Tissue-Aware Conservation Inference Toolkit (TACIT). Rather than measuring the extent to which individual nucleotides are conserved across a region, TACIT uses machine learning to test whether the function of a given part of the genome is likely to be conserved. More specifically, convolutional neural networks learn the tissue- or cell type–specific regulatory code connecting genome sequence to enhancer activity using candidate enhancers identified from only a few species. This approach allows us to accurately associate differences between species in tissue or cell type–specific enhancer activity with genome sequence differences at enhancer orthologs. We then connect these predictions of enhancer function to phenotypes across hundreds of mammals in a way that accounts for species’ phylogenetic relatedness. We applied TACIT to identify candidate enhancers from motor cortex and parvalbumin neuron open chromatin data that are associated with brain size relative to body size, solitary living, and vocal learning across 222 mammals. Our results include the identification of multiple candidate enhancers associated with brain size relative to body size, several of which are located in linear or three-dimensional proximity to genes whose protein-coding mutations have been implicated in microcephaly or macrocephaly in humans. We also identified candidate enhancers associated with the evolution of solitary living near a gene implicated in separation anxiety and other enhancers associated with the evolution of vocal learning ability. We obtained distinct results for bulk motor cortex and parvalbumin neurons, demonstrating the value in applying TACIT to both bulk tissue and specific minority cell type populations. To facilitate future analyses of our results and applications of TACIT, we released predicted enhancer activity of >400,000 candidate enhancers in each of 222 mammals and their associations with the phenotypes we investigated. CONCLUSION TACIT leverages predicted enhancer activity conservation rather than nucleotide-level conservation to connect genetic sequence differences between species to phenotypes across large numbers of mammals. TACIT can be applied to any phenotype with enhancer activity data available from at least a few species in a relevant tissue or cell type and a whole-genome alignment available across dozens of species with substantial phenotypic variation. Although we developed TACIT for transcriptional enhancers, it could also be applied to genomic regions involved in other components of gene regulation, such as promoters and splicing enhancers and silencers. As the number of sequenced genomes grows, machine learning approaches such as TACIT have the potential to help make sense of how conservation of, or changes in, subtle genome patterns can help explain phenotype evolution. Tissue-Aware Conservation Inference Toolkit (TACIT) associates genetic differences between species with phenotypes. TACIT works by generating open chromatin data from a few species in a tissue related to a phenotype, using the sequences underlying open and closed chromatin regions to train a machine learning model for predicting tissue-specific open chromatin and associating open chromatin predictions across dozens of mammals with the phenotype. [Species silhouettes are from PhyloPic]more » « less
