Evolution of hindlimb bone dimensions and muscle masses in house mice selectively bred for high voluntary wheel-running behavior

A meta-analysis of whole-body and heart mass effect sizes from a long-term artificial selection experiment for high voluntary exercise

Selection experiments play an increasingly important role in comparative and evolutionary physiology. However, selection experiments can be limited by relatively low statistical power, in part because replicate line is the experimental unit for analyses of direct or correlated responses (rather than number of individuals measured). One way to increase the ability to detect correlated responses is through a meta-analysis of studies for a given trait across multiple generations. To demonstrate this, we applied meta-analytic techniques to two traits (body mass and heart ventricle mass, with body mass as a covariate) from a long-term artificial selection experiment for high voluntary wheel-running behavior. In this experiment, all four replicate High Runner (HR) lines reached apparent selection limits around generations 17-27, running approximately 2.5- to 3-fold more revolutions per day than the four non-selected Control (C) lines. Although both traits would also be expected to change in HR lines (relative heart size expected to increase, expected direction for body mass is less clear), their statistical significance has varied, despite repeated measurements. We compiled information from 33 unique studies and calculated a measure of effect size (Pearson's R). Our results indicate that, despite a lack of statistical significance in most generations, HR mice have evolved larger hearts and smaller bodies relative to controls. Moreover, plateaus in effect sizes for both traits coincide with the generational range during which the selection limit for wheel-running behavior was reached. Finally, since reaching the selection limit, absolute effect sizes for body mass and heart ventricle mass have become smaller (i.e. closer to 0).

Large changes in detected selection signatures after a selection limit in mice bred for voluntary wheel-running behavior

In various organisms, sequencing of selectively bred lines at apparent selection limits has demonstrated that genetic variation can remain at many loci, implying that evolution at the genetic level may continue even if the population mean phenotype remains constant. We compared selection signatures at generations 22 and 61 of the "High Runner" mouse experiment, which includes 4 replicate lines bred for voluntary wheel-running behavior (HR) and 4 non-selected control (C) lines. Previously, we reported multiple regions of differentiation between the HR and C lines, based on whole-genome sequence data for 10 mice from each line at generation 61, which was >31 generations after selection limits had been reached in all HR lines. Here, we analyzed pooled sequencing data from ~20 mice for each of the 8 lines at generation 22, around when HR lines were reaching limits. Differentiation analyses of allele frequencies at ~4.4 million SNP loci used the regularized T-test and detected 258 differentiated regions with FDR = 0.01. Comparable analyses involving pooling generation 61 individual mouse genotypes into allele frequencies by line produced only 11 such regions, with almost no overlap among the largest and most statistically significant peaks between the two generations. These results implicate a sort of "genetic churn" that continues at loci relevant for running. Simulations indicate that loss of statistical power due to random genetic drift and sampling error are insufficient to explain the differences in selection signatures. The 13 differentiated regions at generation 22 with strict culling measures include 79 genes related to a wide variety of functions. Gene ontology identified pathways related to olfaction and vomeronasal pathways as being overrepresented, consistent with generation 61 analyses, despite those specific regions differing between generations. Genes Dspp and Rbm24 are also identified as potentially explaining known bone and skeletal muscle differences, respectively, between the linetypes.

Does Behavior Evolve First? Correlated Responses to Selection for Voluntary Wheel-Running Behavior in House Mice

AbstractHow traits at multiple levels of biological organization evolve in a correlated fashion in response to directional selection is poorly understood, but two popular models are the very general "behavior evolves first" (BEF) hypothesis and the more specific "morphology-performance-behavior-fitness" (MPBF) paradigm. Both acknowledge that selection often acts relatively directly on behavior and that when behavior evolves, other traits will as well but most with some lag. However, this proposition is exceedingly difficult to test in nature. Therefore, we studied correlated responses in the high-runner (HR) mouse selection experiment, in which four replicate lines have been bred for voluntary wheel-running behavior and compared with four nonselected control (C) lines. We analyzed a wide range of traits measured at generations 20-24 (with a focus on new data from generation 22), coinciding with the point at which all HR lines were reaching selection limits (plateaus). Significance levels (226 P values) were compared across trait types by ANOVA, and we used the positive false discovery rate to control for multiple comparisons. This meta-analysis showed that, surprisingly, the measures of performance (including maximal oxygen consumption during forced exercise) showed no evidence of having diverged between the HR and C lines, nor did any of the life history traits (e.g., litter size), whereas body mass had responded (decreased) at least as strongly as wheel running. Overall, results suggest that the HR lines of mice had evolved primarily by changes in motivation rather than performance ability at the time they were reaching selection limits. In addition, neither the BEF model nor the MPBF model of hierarchical evolution provides a particularly good fit to the HR mouse selection experiment.

Effects of Development on Bone Mineral Density and Mechanical Properties in the Aquatic Frog, Xenopus Laevis, and a Terrestrial Frog, Lithobates Catesbianus

The limb bones of vertebrates have a critical role in supporting the weight of the body and transmitting forces that power locomotion. The loads that limb bones experience can vary in association with a range of factors, including locomotor environment or developmental stage. Limbed vertebrates that are habitually found in environments with low locomotor loads (e.g., water) might be predicted to also exhibit limb bones with less elevated mechanical properties, such as yield stiffness and yield stress. Frogs provide a distinctive case, in which these ideas can be tested as they experience changes in both locomotor style and habitat as they develop. However, while many frog taxa shift from aquatic to terrestrial habitats as they metamorphose, some lineages, such as pipids, maintain an aquatic lifestyle even after metamorphosis, providing a comparative framework for the effects of habitat shifts on developing limbs in vertebrates. This study compares the material composition and mechanical properties of the femur between frog species that are aquatic specialists (Xenopus laevis) vs generalists that spend considerable time both on land and in water (Lithobates catesbeianus) as they transition from metamorphic tadpoles to fully grown adults. MicroCT scanning was used to determine changes in bone density related to developmental stage and hindlimb use during swimming. Microindentation was then used to collect hardness values from the cortical bone of each femur, which was used to evaluate bone material properties. We found that aquatic frogs had less overall bone mineral density (BMD) than terrestrial frogs and that BMD was more elevated in the cortical region of the diaphysis than trabeculae and distal and proximal epiphyses. Despite its less elevated BMD, bone mechanical properties were not significantly different in aquatic specialist X. laevis than in more terrestrial L. catesbeianus. Our results suggest that the limb bones of aquatic frogs may experience compensatory effects through development to offset their lower BMD. Furthermore, changes in bone density and material properties across development may help to explain some of the differences in locomotor performance found between aquatic and terrestrial metamorphic frogs, providing insight into how environmental factors might correlate with bone ossification.

Effects of selective breeding for voluntary exercise, chronic exercise, and their interaction on muscle attachment site morphology in house mice

Skeletal muscles attach to bone at their origins and insertions, and the interface where tendon meets bone is termed the attachment site or enthesis. Mechanical stresses at the muscle/tendon-bone interface are proportional to the surface area of the bony attachment sites, such that a larger attachment site will distribute loads over a wider area. Muscles that are frequently active and/or are of larger size should cause attachment sites to hypertrophy (training effect); however, experimental studies of animals subjected to exercise have provided mixed results. To enhance our ability to detect training effects (a type of phenotypic plasticity), we studied a mouse model in which 4 replicate lines of High Runner (HR) mice have been selectively bred for 57 generations. Selection is based on the average number of wheel revolutions on days 5 & 6 of a 6-day period of wheel access as young adults (6-8 weeks old). Four additional lines are bred without regard to running and serve as non-selected controls (C). On average, mice from HR lines voluntarily run ~3 times more than C mice on a daily basis. For this study, we housed 50 females (half HR, half C) with wheels (Active group) and 50 (half HR, half C) without wheels (Sedentary group) for 12 weeks starting at weaning (~3 weeks old). We tested for evolved differences in muscle attachment site surface area between HR and C mice, plastic changes resulting from chronic exercise, and their interaction. We used a precise, highly repeatable method for quantifying the three-dimensional (3D) surface area of four muscle attachment sites: the humerus deltoid tuberosity (the insertion point for the spinodeltoideus, superficial pectoralis, and acromiodeltoideus), the femoral third trochanter (the insertion point for the quadratus femoris), the femoral lesser trochanter (the insertion point for the iliacus muscle), and the femoral greater trochanter (insertion point for the middle gluteal muscles). In univariate analyses, with body mass as a covariate, mice in the Active group had significantly larger humerus deltoid tuberosities than Sedentary mice, with no significant difference between HR and C mice and no interaction between exercise treatment and linetype. These differences between Active and Sedentary mice were also apparent in the multivariate analyses. Surface areas of the femoral third trochanter, femoral lesser trochanter, and femoral greater trochanter were unaffected by either chronic wheel access or selective breeding. Our results, which used robust measurement protocols and relatively large sample sizes, demonstrate that muscle attachment site morphology can be (but is not always) affected by chronic exercise experienced during ontogeny. However, contrary to previous results for other aspects of long bone morphology, we did not find evidence for evolutionary coadaptation of muscle attachments with voluntary exercise behavior in the HR mice.

Rapid and longer-term effects of selective breeding for voluntary exercise behavior on skeletal morphology in house mice

Selection experiments can elucidate the varying course of adaptive changes across generations. We examined the appendicular skeleton of house mice from four replicate High Runner (HR) lines bred for physical activity on wheels and four non-selected Control (C) lines. HR mice reached apparent selection limits between generations 17 and 27, running ~3-fold more than C. Studies at generations 11, 16, and 21 found that HR mice had evolved thicker hindlimb bones, heavier feet, and larger articular surface areas of the knee and hip joint. Based on biomechanical theory, any or all of these evolved differences may be beneficial for endurance running. Here, we studied mice from generation 68, plus a limited sample from generation 58, to test whether the skeleton continued to evolve after selection limits were reached. Contrary to our expectations, we found few differences between HR and C mice for these later generations, and some of the differences in bone dimensions identified in earlier generations were no longer statistically significant. We hypothesize that the loss of apparently coadapted lower-level traits reflects (1) deterioration related to a gradual increase in inbreeding and/or (2) additional adaptive changes that replace the functional benefits of some skeletal changes.

Early-life effects of juvenile Western diet and exercise on adult gut microbiome composition in mice

Alterations to the gut microbiome caused by changes in diet, consumption of antibiotics, etc., can affect host function. Moreover, perturbation of the microbiome during critical developmental periods potentially has long-lasting impacts on hosts. Using four selectively bred high runner and four non-selected control lines of mice, we examined the effects of early-life diet and exercise manipulations on the adult microbiome by sequencing the hypervariable internal transcribed spacer region of the bacterial gut community. Mice from high runner lines run ∼3-fold more on wheels than do controls, and have several other phenotypic differences (e.g. higher food consumption and body temperature) that could alter the microbiome, either acutely or in terms of coevolution. Males from generation 76 were given wheels and/or a Western diet from weaning until sexual maturity at 6 weeks of age, then housed individually without wheels on standard diet until 14 weeks of age, when fecal samples were taken. Juvenile Western diet reduced bacterial richness and diversity after the 8-week washout period (equivalent to ∼6 human years). We also found interactive effects of genetic line type, juvenile diet and/or juvenile exercise on microbiome composition and diversity. Microbial community structure clustered significantly in relation to both line type and diet. Western diet also reduced the relative abundance of Muribaculum intestinale These results constitute one of the first reports of juvenile diet having long-lasting effects on the adult microbiome after a substantial washout period. Moreover, we found interactive effects of diet with early-life exercise exposure, and a dependence of these effects on genetic background.

Genetic Basis of Aerobically Supported Voluntary Exercise: Results from a Selection Experiment with House Mice

The biological basis of exercise behavior is increasingly relevant for maintaining healthy lifestyles. Various quantitative genetic studies and selection experiments have conclusively demonstrated substantial heritability for exercise behavior in both humans and laboratory rodents. In the "High Runner" selection experiment, four replicate lines of Mus domesticus were bred for high voluntary wheel running (HR), along with four nonselected control (C) lines. After 61 generations, the genomes of 79 mice (9-10 from each line) were fully sequenced and single nucleotide polymorphisms (SNPs) were identified. We used nested ANOVA with MIVQUE estimation and other approaches to compare allele frequencies between the HR and C lines for both SNPs and haplotypes. Approximately 61 genomic regions, across all somatic chromosomes, showed evidence of differentiation; 12 of these regions were differentiated by all methods of analysis. Gene function was inferred largely using Panther gene ontology terms and KO phenotypes associated with genes of interest. Some of the differentiated genes are known to be associated with behavior/motivational systems and/or athletic ability, including Sorl1, Dach1, and Cdh10 Sorl1 is a sorting protein associated with cholinergic neuron morphology, vascular wound healing, and metabolism. Dach1 is associated with limb bud development and neural differentiation. Cdh10 is a calcium ion binding protein associated with phrenic neurons. Overall, these results indicate that selective breeding for high voluntary exercise has resulted in changes in allele frequencies for multiple genes associated with both motivation and ability for endurance exercise, providing candidate genes that may explain phenotypic changes observed in previous studies.

Trade-offs, Pleiotropy, and Shared Molecular Pathways: A Unified View of Constraints on Adaptation

The concept of trade-offs permeates our thinking about adaptive evolution because they are exhibited at every level of biological organization, from molecular and cellular processes to organismal and ecological functions. Trade-offs inevitably arise because different traits do not occur in isolation, but instead are imbedded within complex, integrated systems that make up whole organisms. The genetic and mechanistic underpinning of trade-offs can be found in the pleiotropic nodes that occur in the biological pathways shared between traits. Yet, often trade-offs are only understood as statistical correlations, limiting the ability to evaluate the interplay between how selection and constraint interact during adaptive evolution. Here, we first review the classic paradigms in which physiologists and evolutionary biologists have studied trade-offs and highlight the ways in which network and molecular pathway approaches unify these paradigms. We discuss how these approaches allow researchers to evaluate why trade-offs arise and how selection can act to overcome trait correlations and evolutionary constraints. We argue that understanding how the conserved molecular pathways are shared between different traits and functions provides a conceptual framework for evolutionary biologists, physiologists, and molecular biologists to meaningfully work together toward the goal of understanding why correlations and trade-offs occur between traits. We briefly highlight the melanocortin system and the hormonal control of osmoregulation as two case studies where an understanding of shared molecular pathways reveals why trade-offs occur between seemingly unrelated traits. While we recognize that applying such approaches poses challenges and limitations particularly in the context of natural populations, we advocate for the view that focusing on the biological pathways responsible for trade-offs provides a unified conceptual context accessible to a broad range of integrative biologists.

Patterns of fluctuating asymmetry in the limbs of freshwater turtles: Are more functionally important limbs more symmetrical?

Understanding how selective forces influence patterns of symmetry remains an active area of research in evolutionary biology. One hypothesis, which has received relatively little attention, suggests that the functional importance of morphological characters may influence patterns of symmetry. Specifically, it posits that for structures that display bilateral symmetry, those with greater functional importance should display lower levels of asymmetry. The aim of this study was to examine the patterns of fluctuating asymmetry (FA) present in the limb bones of freshwater turtles in the family Emydidae. Aquatic emydid turtles of the subfamily Deirochelyinae employ a hindlimb-dominant swimming style, suggesting that hindlimbs should display lower levels of FA. Consistent with the morpho-functional hypothesis of symmetry, we found a strong, clade-wise pattern of humeral-biased FA in aquatic Deirochelyinae. In contrast, some emydids of the subfamily Emydinae possess more terrestrial tendencies. As terrestrial locomotion places more equal importance on fore- and hindlimbs, we predicted that such behaviors may minimize differences in FA. No clade-wise pattern was detected in the subfamily Emydinae. We also detected a phylogenetic signal in FA within the femur and discovered that FA has evolved at vastly different rates between the fore- and hindlimbs.

Patterns of fluctuating asymmetry in the limbs of anurans

It has been hypothesized that fluctuating asymmetry (FA) may provide an indication of the functional importance of structures within an organism, with structures that more strongly impact fitness being more symmetric. Based on this idea, we predicted that for tetrapods in which the forelimbs and hindlimbs play an unequal role in locomotion, the less functionally important limb set should display higher levels of FA. We conducted a multispecies test of this hypothesis in anurans (frogs and toads), whose saltatory locomotor mode is powered by the hindlimbs. We also tested whether FA in the forelimbs, which play a more important role during landing, differed between families that differ in the degree of forelimb use in locomotion (Bufonidae vs. Ranidae). We calculated FA from the lengths of humeri and femora measured from disarticulated skeletal specimens of four anuran taxa (Bufonidae: Anaxyrus americanus, Rhinella marina; Ranidae: Lithobates catesbeianus, Lithobates clamitans). Our findings were consistent with the hypothesis that natural selection for increased locomotor performance may influence patterns of FA seen in vertebrate limbs, with all species displaying lower mean FA in the hindlimbs. More subtle functional roles between the forelimbs of bufonids and ranids, however, did not elicit different levels of FA.

Effects of selective breeding for high voluntary wheel-running behavior on femoral nutrient canal size and abundance in house mice

Bone modeling and remodeling are aerobic processes that entail relatively high oxygen demands. Long bones receive oxygenated blood from nutrient arteries, epiphyseal-metaphyseal arteries, and periosteal arteries, with the nutrient artery supplying the bulk of total blood volume in mammals (~ 50-70%). Estimates of blood flow into these bones can be made from the dimensions of the nutrient canal, through which nutrient arteries pass. Unfortunately, measuring these canal dimensions non-invasively (i.e. without physical sectioning) is difficult, and thus researchers have relied on more readily visible skeletal proxies. Specifically, the size of the nutrient artery has been estimated from dimensions (e.g. minimum diameters) of the periosteal (external) opening of the nutrient canal. This approach has also been utilized by some comparative morphologists and paleontologists, as the opening of a nutrient canal is present long after the vascular soft tissue has degenerated. The literature on nutrient arteries and canals is sparse, with most studies consisting of anatomical descriptions from surgical proceedings, and only a few investigating the links between nutrient canal morphology and physiology or behavior. The primary objective of this study was to evaluate femur nutrient canal morphology in mice with known physiological and behavioral differences; specifically, mice from an artificial selection experiment for high voluntary wheel-running behavior. Mice from four replicate high runner (HR) lines are known to differ from four non-selected control (C) lines in both locomotor and metabolic activity, with HR mice having increased voluntary wheel-running behavior and maximal aerobic capacity (VO₂ max) during forced treadmill exercise. Femora from adult mice (average age 7.5 months) of the 11th generation of this selection experiment were μCT-scanned and three-dimensional virtual reconstructions of nutrient canals were measured for minimum cross-sectional area as a skeletal proxy of blood flow. Gross observations revealed that nutrient canals varied far more in number and shape than prior descriptions would indicate, regardless of sex or genetic background (i.e. HR vs. C lines). Canals adopted non-linear shapes and paths as they traversed from the periosteal to endosteal borders through the cortex, occasionally even branching within the cortical bone. Additionally, mice from both HR and C lines averaged more than four nutrient canals per femur, in contrast to the one to two nutrient canals described for femora from rats, pigs, and humans in prior literature. Mice from HR lines had significantly larger total nutrient canal area than C lines, which was the result not of an increase in the number of nutrient canals, but rather an increase in their average cross-section size. This study demonstrates that mice with an evolutionary history of increased locomotor activity and maximal aerobic metabolic rate have a concomitant increase in the size of their femoral nutrient canals. Although the primary determinant of nutrient canal size is currently not well understood, the present results bolster use of nutrient canal size as a skeletal indicator of aerobically supported levels of physical activity in comparative studies.