Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice.

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models.

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation.

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

The Contribution of Genetic and Environmental Effects to Bergmann's Rule and Allen's Rule in House Mice.

AbstractDistinguishing between genetic, environmental, and genotype × environment effects is central to understanding geographic variation in phenotypic clines. Two of the best-documented phenotypic clines are Bergmann's rule and Allen's rule, which describe larger body sizes and shortened extremities in colder climates, respectively. Although numerous studies have found inter- and intraspecific evidence for both ecogeographic patterns, we still have a poor understanding of the extent to which these patterns are driven by genetics, environment, or both. Here, we measured the genetic and environmental contributions to Bergmann's rule and Allen's rule across introduced populations of house mice (Mus musculus domesticus) in the Americas. First, we documented clines for body mass, tail length, and ear length in natural populations and found that these conform to both Bergmann's rule and Allen's rule. We then raised descendants of wild-caught mice in the lab and showed that these differences persisted in a common environment and are heritable, indicating that they have a genetic basis. Finally, using a full-sib design, we reared mice under warm and cold conditions. We found very little plasticity associated with body size, suggesting that Bergmann's rule has been shaped by strong directional selection in house mice. However, extremities showed considerable plasticity, as both tails and ears grew shorter in cold environments. These results indicate that adaptive phenotypic plasticity as well as genetic changes underlie major patterns of clinal variation in house mice and likely facilitated their rapid expansion into new environments across the Americas.

Are Lethal Alleles Too Abundant in Humans?

Across species, many individuals carry one or more recessive lethal alleles, posing an evolutionary conundrum for their persistence. Using a population genomic approach, Amorim et al. studied the abundance of lethal disease-causing mutations in humans and found that, while appearing more common than expected, most may nonetheless persist at frequencies predicted by mutation-selection balance.

Gene regulation underlies environmental adaptation in house mice.

Changes in cis-regulatory regions are thought to play a major role in the genetic basis of adaptation. However, few studies have linked cis-regulatory variation with adaptation in natural populations. Here, using a combination of exome and RNA-seq data, we performed expression quantitative trait locus (eQTL) mapping and allele-specific expression analyses to study the genetic architecture of regulatory variation in wild house mice (Mus musculus domesticus) using individuals from five populations collected along a latitudinal cline in eastern North America. Mice in this transect showed clinal patterns of variation in several traits, including body mass. Mice were larger in more northern latitudes, in accordance with Bergmann's rule. We identified 17 genes where cis-eQTLs were clinal outliers and for which expression level was correlated with latitude. Among these clinal outliers, we identified two genes (Adam17 and Bcat2) with cis-eQTLs that were associated with adaptive body mass variation and for which expression is correlated with body mass both within and between populations. Finally, we performed a weighted gene co-expression network analysis (WGCNA) to identify expression modules associated with measures of body size variation in these mice. These findings demonstrate the power of combining gene expression data with scans for selection to identify genes involved in adaptive phenotypic evolution, and also provide strong evidence for cis-regulatory elements as essential loci of environmental adaptation in natural populations.

Nature's fat-burning machine: brown adipose tissue in a hibernating mammal.

Brown adipose tissue (BAT) is a unique thermogenic tissue in mammals that rapidly produces heat via nonshivering thermogenesis. Small mammalian hibernators have evolved the greatest capacity for BAT because they use it to rewarm from hypothermic torpor numerous times throughout the hibernation season. Although hibernator BAT physiology has been investigated for decades, recent efforts have been directed toward understanding the molecular underpinnings of BAT regulation and function using a variety of methods, from mitochondrial functional assays to 'omics' approaches. As a result, the inner-workings of hibernator BAT are now being illuminated. In this Review, we discuss recent research progress that has identified players and pathways involved in brown adipocyte differentiation and maturation, as well as those involved in metabolic regulation. The unique phenotype of hibernation, and its reliance on BAT to generate heat to arouse mammals from torpor, has uncovered new molecular mechanisms and potential strategies for biomedical applications.

Enhanced oxidative capacity of ground squirrel brain mitochondria during hibernation.

During hibernation, thirteen-lined ground squirrels (Ictidomys tridecemlineatus) regularly cycle between bouts of torpor and interbout arousal (IBA). Most of the brain is electrically quiescent during torpor but regains activity quickly upon arousal to IBA, resulting in extreme oscillations in energy demand during hibernation. We predicted increased functional capacity of brain mitochondria during hibernation compared with spring to accommodate the variable energy demands of hibernation. To address this hypothesis, we examined mitochondrial bioenergetics in the ground squirrel brain across three time points: spring (SP), torpor (TOR), and IBA. Respiration rates of isolated brain mitochondria through complex I of the electron transport chain were more than twofold higher in TOR and IBA than in SP (P < 0.05). We also found a 10% increase in membrane potential between hibernation and spring (P < 0.05), and that proton leak was lower in TOR and IBA than in SP. Finally, there was a 30% increase in calcium loading in SP brain mitochondria compared with TOR and IBA (P < 0.01). To analyze brain mitochondrial abundance between spring and hibernation, we measured the ratio of copy number in a mitochondrial gene (ND1) vs. a nuclear gene (B2M) in frozen cerebral cortex samples. No significant differences were observed in DNA copies between SP and IBA. These data show that brain mitochondrial bioenergetics are not static across the year and suggest that brain mitochondria function more effectively during the hibernation season, allowing for rapid production of energy to meet demand when extreme physiological changes are occurring.

Seasonal changes in brown adipose tissue mitochondria in a mammalian hibernator: from gene expression to function.

Brown adipose tissue (BAT) is a thermogenic organ that is vital for hibernation in mammals. Throughout the hibernation season, BAT mitochondrial uncoupling protein 1 (UCP1) enables rapid rewarming from hypothermic torpor to periodic interbout arousals (IBAs), as energy is dissipated as heat. However, BAT's unique ability to rewarm the body via nonshivering thermogenesis is not necessary outside the hibernation season, suggesting a potential seasonal change in the regulation of BAT function. Here, we examined the BAT mitochondrial proteome and mitochondrial bioenergetics in the thirteen-lined ground squirrel (Ictidomys tridecemlineatus) across four time points: spring, fall, torpor, and IBA. Relative mitochondrial content of BAT was estimated by measuring BAT pad mass, UCP1 protein content, and mitochondrial DNA (mtDNA) copy number. BAT mtDNA content was significantly lower in spring compared with torpor and IBA (P < 0.05). UCP1 mRNA and protein levels were highest during torpor and IBA. Respiration rates of isolated BAT mitochondria were interrogated at each complex of the electron transport chain. Respiration at complex II was significantly higher in torpor and IBA compared with spring (P < 0.05), suggesting an enhancement in mitochondrial respiratory capacity during hibernation. Additionally, proteomic iTRAQ labeling identified 778 BAT mitochondrial proteins. Proteins required for mitochondrial lipid translocation and ß-oxidation were upregulated during torpor and IBA and downregulated in spring. These data imply that BAT bioenergetics and mitochondrial content are not static across the year, despite the year-round presence of UCP1.

Melatonin receptor signaling contributes to neuroprotection upon arousal from torpor in thirteen-lined ground squirrels.

The brain of mammalian hibernators is naturally protected. Hibernating ground squirrels undergo rapid and extreme changes in body temperature and brain perfusion as they cycle between lengthy torpor bouts and brief periods of euthermia called interbout arousals (IBAs). Arousal from torpor to IBA occurs rapidly, but there is no evidence of brain injury accompanying this extreme physiological transition. Production of the hormone melatonin accompanies arousal, suggesting that it plays a protective role at this time. Here, we investigated mechanisms of melatonin receptor-mediated protection in the brain of the hibernating ground squirrel. We administered the competitive melatonin receptor antagonist luzindole (30 mg/kg ip) to ground squirrels at the predicted end of a torpor bout, triggering an arousal. We found that luzindole-treated animals exhibited caspase-3 activity two times higher than vehicle-treated animals in the hypothalamus at midarousal (P = 0.01), suggesting that melatonin receptor signaling is important for protection in this brain region. We also found a 30% decline in succinate-fueled mitochondrial respiration in luzindole-treated animals compared with vehicle-treated animals (P = 0.019), suggesting that melatonin receptor signaling is important for optimal mitochondrial function during arousal from torpor. The mitochondrial effects of luzindole treatment were seen only during the hibernation season, indicating that this effect is specifically important for arousal from torpor. These data provide evidence for the protective role of melatonin receptor signaling during the extreme physiological transition that occurs when a hibernating mammal arouses from torpor and provide further evidence for regional and seasonal changes in the hibernator brain.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models.

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

The genomic basis of high-elevation adaptation in wild house mice (Mus musculus domesticus) from South America.

Understanding the genetic basis of environmental adaptation in natural populations is a central goal in evolutionary biology. The conditions at high elevation, particularly the low oxygen available in the ambient air, impose a significant and chronic environmental challenge to metabolically active animals with lowland ancestry. To understand the process of adaptation to these novel conditions and to assess the repeatability of evolution over short timescales, we examined the signature of selection from complete exome sequences of house mice (Mus musculus domesticus) sampled across two elevational transects in the Andes of South America. Using phylogenetic analysis, we show that house mice colonized high elevations independently in Ecuador and Bolivia. Overall, we found distinct responses to selection in each transect and largely nonoverlapping sets of candidate genes, consistent with the complex nature of traits that underlie adaptation to low oxygen availability (hypoxia) in other species. Nonetheless, we also identified a small subset of the genome that appears to be under parallel selection at the gene and SNP levels. In particular, three genes (Col22a1, Fgf14, and srGAP1) bore strong signatures of selection in both transects. Finally, we observed several patterns that were common to both transects, including an excess of derived alleles at high elevation, and a number of hypoxia-associated genes exhibiting a threshold effect, with a large allele frequency change only at the highest elevations. This threshold effect suggests that selection pressures may increase disproportionately at high elevations in mammals, consistent with observations of some high-elevation diseases in humans.