Impacts of genetic correlation on the independent evolution of body mass and skeletal size in mammals

Shared Selection and Genetic Architecture Drive Strikingly Repeatable Evolution in Long-Term Experimental Hybrid Populations

The degree to which evolution repeats itself has implications regarding the major forces driving evolution and the potential for evolutionary biology to be a predictive (vs. solely historical) science. To understand the factors that control evolutionary repeatability, we experimentally evolved four replicate hybrid populations of sunflowers at natural sites for up to 14 years and tracked ancestry across the genome. We found that there was very strong negative selection against introgressed ancestry in several chromosomes, but positive selection for introgressed ancestry in one chromosome. Further, the strength of selection was influenced by recombination rate. High recombination regions had lower selection against introgressed ancestry due to more frequent recombination away from incompatible backgrounds. Strikingly, evolution was highly parallel across replicates, with shared selection driving 88% of variance in introgressed allele frequency change. Parallel evolution was driven by both high levels of sustained linkage in introgressed alleles and strong selection on large-effect quantitative trait loci. This work highlights the repeatability of evolution through hybridization and confirms the central roles that natural selection, genomic architecture, and recombination play in the process.

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).

Impact of hindlimb length variation on jumping dynamics in the Longshanks mouse

Distantly related mammals (e.g. jerboa, tarsiers, kangaroos) have convergently evolved elongated hindlimbs relative to body size. Limb elongation is hypothesized to make these species more effective jumpers by increasing their kinetic energy output (through greater forces or acceleration distances), thereby increasing take-off velocity and jump distance. This hypothesis, however, has rarely been tested at the population level, where natural selection operates. We examined the relationship between limb length, muscular traits and dynamics using Longshanks mice, which were selectively bred over 22 generations for longer tibiae. Longshanks mice have approximately 15% longer tibiae and 10% longer femora compared with random-bred Control mice from the same genetic background. We collected in vivo measures of locomotor kinematics and force production, in combination with behavioral data and muscle morphology, to examine how changes in bone and muscle structure observed in Longshanks mice affect their hindlimb dynamics during jumping and clambering. Longshanks mice achieved higher mean and maximum lunge-jump heights than Control mice. When jumping to a standardized height (14 cm), Longshanks mice had lower maximum ground reaction forces, prolonged contact times and greater impulses, without significant differences in average force, power or whole-body velocity. While Longshanks mice have longer plantarflexor muscle bodies and tendons than Control mice, there were no consistent differences in muscular cross-sectional area or overall muscle volume; improved lunge-jumping performance in Longshanks mice is not accomplished by simply possessing larger muscles. Independent of other morphological or behavioral changes, our results point to the benefit of longer hindlimbs for performing dynamic locomotion.

Histone deacetylases regulate organ-specific growth in a horned beetle

BACKGROUND

Nutrient availability is among the most widespread means by which environmental variability affects developmental outcomes. Because almost all cells within an individual organism share the same genome, structure-specific growth responses must result from changes in gene regulation. Earlier work suggested that histone deacetylases (HDACs) may serve as epigenetic regulators linking nutritional conditions to trait-specific development. Here we expand on this work by assessing the function of diverse HDACs in the structure-specific growth of both sex-shared and sex-specific traits including evolutionarily novel structures in the horned dung beetle Onthophagus taurus.

RESULTS

We identified five HDAC members whose downregulation yielded highly variable mortality depending on which HDAC member was targeted. We then show that HDAC1, 3, and 4 operate in both a gene- and trait-specific manner in the regulation of nutrition-responsiveness of appendage size and shape. Specifically, HDAC 1, 3, or 4 knockdown diminished wing size similarly while leg development was differentially affected by RNAi targeting HDAC3 and HDAC4. In addition, depletion of HDAC3 transcript resulted in a more rounded shape of genitalia at the pupal stage and decreased the length of adult aedeagus across all body sizes. Most importantly, we find that HDAC3 and HDAC4 pattern the morphology and regulate the scaling of evolutionarily novel head and thoracic horns as a function of nutritional variation.

CONCLUSION

Collectively, our results suggest that both functional overlap and division of labor among HDAC members contribute to morphological diversification of both conventional and recently evolved appendages. More generally, our work raises the possibility that HDAC-mediated scaling relationships and their evolution may underpin morphological diversification within and across insect species broadly.

The genetic architecture and evolution of the human skeletal form

The human skeletal form underlies bipedalism, but the genetic basis of skeletal proportions (SPs) is not well characterized. We applied deep-learning models to 31,221 x-rays from the UK Biobank to extract a comprehensive set of SPs, which were associated with 145 independent loci genome-wide. Structural equation modeling suggested that limb proportions exhibited strong genetic sharing but were independent of width and torso proportions. Polygenic score analysis identified specific associations between osteoarthritis and hip and knee SPs. In contrast to other traits, SP loci were enriched in human accelerated regions and in regulatory elements of genes that are differentially expressed between humans and great apes. Combined, our work identifies specific genetic variants that affect the skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.

Adapting the pantograph limb: Differential robustness of fore- and hindlimb kinematics against genetically induced perturbation in the neural control networks and its evolutionary implications

The evolutionary transformation of limb morphology to the four-segmented pantograph of therians is among the milestones of mammalian evolution. But, it is still unknown if changes of the mechanical limb function were accompanied by corresponding changes in development and sensorimotor control. The impressive locomotor performance of mammals leaves no doubt about the high integration of pattern formation, neural control and mechanics. But, deviations from normal intra- and interlimb coordination (spatial and temporal) become evident in the presence of perturbations. We induced a perturbation in the development of the neural circuits of the spinal cord of mice (Mus musculus) using a deletion of the Wilms tumor suppressor gene Wt1 in a subpopulation of dI6 interneurons. These interneurons are assumed to participate in the intermuscular coordination within the limb and in left-right-coordination between the limbs. We describe the locomotor kinematics in mice with conditional Wt1 knockout and compare them to mice without Wt1 deletion. Unlike knockout neonates, knockout adult mice do not display severe deviations from normal (=control group) interlimb coordination, but the coordinated protraction and retraction of the limbs is altered. The forelimbs are more affected by deviations from the control than the hindlimbs. This observation appears to reflect a different degree of integration and resistance against the induced perturbation between the limbs. Interestingly, the observed effects are similar to locomotor deficits reported to arise when sensory feedback from proprioceptors or cutaneous receptors is impaired. A putative participation of Wt1 positive dI6 interneurons in sensorimotor integration is therefore considered.

The genetic architecture of the human skeletal form

The human skeletal form underlies our ability to walk on two legs, but unlike standing height, the genetic basis of limb lengths and skeletal proportions is less well understood. Here we applied a deep learning model to 31,221 whole body dual-energy X-ray absorptiometry (DXA) images from the UK Biobank (UKB) to extract 23 different image-derived phenotypes (IDPs) that include all long bone lengths as well as hip and shoulder width, which we analyzed while controlling for height. All skeletal proportions are highly heritable (∼40-50%), and genome-wide association studies (GWAS) of these traits identified 179 independent loci, of which 102 loci were not associated with height. These loci are enriched in genes regulating skeletal development as well as associated with rare human skeletal diseases and abnormal mouse skeletal phenotypes. Genetic correlation and genomic structural equation modeling indicated that limb proportions exhibited strong genetic sharing but were genetically independent of width and torso proportions. Phenotypic and polygenic risk score analyses identified specific associations between osteoarthritis (OA) of the hip and knee, the leading causes of adult disability in the United States, and skeletal proportions of the corresponding regions. We also found genomic evidence of evolutionary change in arm-to-leg and hip-width proportions in humans consistent with striking anatomical changes in these skeletal proportions in the hominin fossil record. In contrast to cardiovascular, auto-immune, metabolic, and other categories of traits, loci associated with these skeletal proportions are significantly enriched in human accelerated regions (HARs), and regulatory elements of genes differentially expressed through development between humans and the great apes. Taken together, our work validates the use of deep learning models on DXA images to identify novel and specific genetic variants affecting the human skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.

Trade-offs in muscle physiology in selectively bred high runner mice

A trade-off between locomotor speed and endurance occurs in various taxa, and is thought to be underpinned by a muscle-level trade-off. Among four replicate high runner (HR) lines of mice, selectively bred for voluntary wheel-running behavior, a negative correlation between average running speed and time spent running has evolved. We hypothesize that this trade-off is due to changes in muscle physiology. We studied the HR lines at generation 90, at which time one line (L3) is fixed for the mini-muscle phenotype, another is polymorphic (L6) and the others (L7, L8) lack mini-muscle individuals. We used in situ preparations to quantify the contractile properties of the triceps surae muscle complex. Maximal shortening velocity varied significantly, being lowest in mini-muscle mice (L3 mini=25.2 mm s-1, L6 mini=25.5 mm s-1), highest in normal-muscle mice L6 and L8 (40.4 and 50.3 mm s-1, respectively) and intermediate in normal-muscle L7 mice (37.2 mm s-1). Endurance, measured both as the slope of the decline in force and the proportion of initial force that could be sustained, also varied significantly. The slope was shallowest in mini-muscle mice (L3 mini=-0.00348, L6 mini=-0.00238), steepest in lines L6 and L8 (-0.01676 and -0.01853), and intermediate in L7 (-0.01145). Normalized sustained force was highest in mini-muscle mice (L3 mini=0.98, L6 mini=0.92) and lowest in L8 (0.36). There were significant, negative correlations between velocity and endurance metrics, indicating a muscle-level trade-off. However, this muscle-level trade-off does not seem to underpin the organismal-level speed and endurance trade-off previously reported as the ordering of the lines is reversed: the lines that run the fastest for the least time have the lowest muscle complex velocity and highest endurance.

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.

A Na+/K+ ATPase Pump Regulates Chondrocyte Differentiation and Bone Length Variation in Mice

The genetic and developmental mechanisms involved in limb formation are relatively well documented, but how these mechanisms are modulated by changes in chondrocyte physiology to produce differences in limb bone length remains unclear. Here, we used high throughput RNA sequencing (RNAseq) to probe the developmental genetic basis of variation in limb bone length in Longshanks, a mouse model of experimental evolution. We find that increased tibia length in Longshanks is associated with altered expression of a few key endochondral ossification genes such as Npr3, Dlk1, Sox9, and Sfrp1, as well reduced expression of Fxyd2, a facultative subunit of the cell membrane-bound Na⁺/K⁺ ATPase pump (NKA). Next, using murine tibia and cell cultures, we show a dynamic role for NKA in chondrocyte differentiation and in bone length regulation. Specifically, we show that pharmacological inhibition of NKA disrupts chondrocyte differentiation, by upregulating expression of mesenchymal stem cell markers (Prrx1, Serpina3n), downregulation of chondrogenesis marker Sox9, and altered expression of extracellular matrix genes (e.g., collagens) associated with proliferative and hypertrophic chondrocytes. Together, Longshanks and in vitro data suggest a broader developmental and evolutionary role of NKA in regulating limb length diversity.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Selection for increased tibia length in mice alters skull shape through parallel changes in developmental mechanisms

Bones in the vertebrate cranial base and limb skeleton grow by endochondral ossification, under the control of growth plates. Mechanisms of endochondral ossification are conserved across growth plates, which increases covariation in size and shape among bones, and in turn may lead to correlated changes in skeletal traits not under direct selection. We used micro-CT and geometric morphometrics to characterize shape changes in the cranium of the Longshanks mouse, which was selectively bred for longer tibiae. We show that Longshanks skulls became longer, flatter, and narrower in a stepwise process. Moreover, we show that these morphological changes likely resulted from developmental changes in the growth plates of the Longshanks cranial base, mirroring changes observed in its tibia. Thus, indirect and non-adaptive morphological changes can occur due to developmental overlap among distant skeletal elements, with important implications for interpreting the evolutionary history of vertebrate skeletal form.

Multiple haplotype reconstruction from allele frequency data

Because haplotype information is of widespread interest in biomedical applications, effort has been put into their reconstruction. Here, we propose an efficient method, called haploSep, that is able to accurately infer major haplotypes and their frequencies just from multiple samples of allele frequency data. Even the accuracy of experimentally obtained allele frequencies can be improved by re-estimating them from our reconstructed haplotypes. From a methodological point of view, we model our problem as a multivariate regression problem where both the design matrix and the coefficient matrix are unknown. Compared to other methods, haploSep is very fast, with linear computational complexity in the haplotype length. We illustrate our method on simulated and real data focusing on experimental evolution and microbial data.

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.

Estimating the genome-wide contribution of selection to temporal allele frequency change

Rapid phenotypic adaptation is often observed in natural populations and selection experiments. However, detecting the genome-wide impact of this selection is difficult since adaptation often proceeds from standing variation and selection on polygenic traits, both of which may leave faint genomic signals indistinguishable from a noisy background of genetic drift. One promising signal comes from the genome-wide covariance between allele frequency changes observable from temporal genomic data (e.g., evolve-and-resequence studies). These temporal covariances reflect how heritable fitness variation in the population leads changes in allele frequencies at one time point to be predictive of the changes at later time points, as alleles are indirectly selected due to remaining associations with selected alleles. Since genetic drift does not lead to temporal covariance, we can use these covariances to estimate what fraction of the variation in allele frequency change through time is driven by linked selection. Here, we reanalyze three selection experiments to quantify the effects of linked selection over short timescales using covariance among time points and across replicates. We estimate that at least 17 to 37% of allele frequency change is driven by selection in these experiments. Against this background of positive genome-wide temporal covariances, we also identify signals of negative temporal covariance corresponding to reversals in the direction of selection for a reasonable proportion of loci over the time course of a selection experiment. Overall, we find that in the three studies we analyzed, linked selection has a large impact on short-term allele frequency dynamics that is readily distinguishable from genetic drift.

Can behaviour impede evolution? Persistence of singing effort after morphological song loss in crickets

Evolutionary loss of sexual signals is widespread. Examining the consequences for behaviours associated with such signals can provide insight into factors promoting or inhibiting trait loss. We tested whether a behavioural component of a sexual trait, male calling effort, has been evolutionary reduced in silent populations of Hawaiian field crickets (Teleogryllus oceanicus). Cricket song requires energetically costly wing movements, but 'flatwing' males have feminized wings that preclude song and protect against a lethal, eavesdropping parasitoid. Flatwing males express wing movement patterns associated with singing but, in contrast with normal-wing males, sustained periods of wing movement cannot confer sexual selection benefits and should be subject to strong negative selection. We developed an automated technique to quantify how long males spend expressing wing movements associated with song. We compared calling effort among populations of Hawaiian crickets with differing proportions of silent males and between male morphs. Contrary to expectation, silent populations invested as much in calling effort as non-silent populations. Additionally, flatwing and normal-wing males from the same population did not differ in calling effort. The lack of evolved behavioural adjustment following morphological change in silent Hawaiian crickets illustrates how behaviour might sometimes impede, rather than facilitate, evolution.

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

Endochondral ossification and the evolution of limb proportions

Mammals have remarkably diverse limb proportions hypothesized to have evolved adaptively in the context of locomotion and other behaviors. Mechanistically, evolutionary diversity in limb proportions is the result of differential limb bone growth. Longitudinal limb bone growth is driven by the process of endochondral ossification, under the control of the growth plates. In growth plates, chondrocytes undergo a tightly orchestrated life cycle of proliferation, matrix production, hypertrophy, and cell death/transdifferentiation. This life cycle is highly conserved, both among the long bones of an individual, and among homologous bones of distantly related taxa, leading to a finite number of complementary cell mechanisms that can generate heritable phenotype variation in limb bone size and shape. The most important of these mechanisms are chondrocyte population size in chondrogenesis and in individual growth plates, proliferation rates, and hypertrophic chondrocyte size. Comparative evidence in mammals and birds suggests the existence of developmental biases that favor evolutionary changes in some of these cellular mechanisms over others in driving limb allometry. Specifically, chondrocyte population size may evolve more readily in response to selection than hypertrophic chondrocyte size, and extreme hypertrophy may be a rarer evolutionary phenomenon associated with highly specialized modes of locomotion in mammals (e.g., powered flight, ricochetal bipedal hopping). Physical and physiological constraints at multiple levels of biological organization may also have influenced the cell developmental mechanisms that have evolved to produce the highly diverse limb proportions in extant mammals. This article is categorized under: Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Comparative Development and Evolution > Regulation of Organ Diversity Comparative Development and Evolution > Organ System Comparisons Between Species.

Testing for correlation between traits under directional evolution

Being confounding factors, directional trends are likely to make two quantitative traits appear as spuriously correlated. By determining the probability distributions of independent contrasts when traits evolve following Brownian motions with linear trends, we show that the standard independent contrasts can not be used to test for correlation in this situation. We propose a multiple regression approach which corrects the bias caused by directional evolution. We show that our approach is equivalent to performing a Phylogenetic Generalized Least Squares (PGLS) analysis with tip times as covariables by providing a new and more general proof of the equivalence between PGLS and independent contrasts methods. Our approach is assessed and compared with three previous correlation tests on data simulated in various situations and overall outperforms all the other methods. The approach is next illustrated on a real dataset to test for correlation between hominin cranial capacity and body mass.

Morphology and development of a novel murine skeletal dysplasia

Background

Limb bones develop and grow by endochondral ossification, which is regulated by specific cell and molecular pathways. Changes in one or more of these pathways can have severe effects on normal skeletal development, leading to skeletal dysplasias. Many skeletal dysplasias are known to result from mis-expression of major genes involved in skeletal development, but the etiology of many skeletal dysplasias remains unknown. We investigated the morphology and development of a mouse line with an uncharacterized mutation exhibiting a skeletal dysplasia-like phenotype (Nabo).

Methods

We used µCT scanning and histology to comprehensively characterize the phenotype and its development, and to determine the developmental stage when this phenotype first appears.

Results

Nabo mice have shorter limb elements compared to wildtype mice, while clavicles and dermal bones of the skull are not affected. Nabo embryos at embryonic stage E14 show shorter limb cartilage condensations. The tibial growth plate in Nabo mice is wider than in wildtype, particularly in the proliferative zone, however proliferative chondrocytes show less activity than wildtype mice. Cell proliferation assays and immunohistochemistry against the chondrogenic marker Sox9 suggest relatively lower, spatially-restricted, chondrocyte proliferation activity in Nabo. Bone volume and trabecular thickness in Nabo tibiae are also decreased compared to wildtype.

Discussion

Our data suggest that the Nabo mutation affects endochondral ossification only, with the strongest effects manifesting in more proximal limb structures. The phenotype appears before embryonic stage E14, suggesting that outgrowth and patterning processes may be affected. Nabo mice present a combination of skeletal dysplasia-like characteristics not present in any known skeletal dysplasia. Further genomic and molecular analysis will help to identify the genetic basis and precise developmental pathways involved in this unique skeletal dysplasia.

An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Selection for longer limbs in mice increases bone stiffness and brittleness, but does not alter bending strength

The ability of a bone to withstand loads depends on its structural and material properties. These tend to differ among species with different modes of locomotion, reflecting their unique loading patterns. The evolution of derived limb morphologies, such as the long limbs associated with jumping, may compromise overall bone strength. We evaluated bone mechanical properties in the Longshanks mouse, which was selectively bred for increased tibia length relative to body mass. We combined analyses of 3D shape and cross-sectional geometry of the tibia, with mechanical testing and bone composition assays, to compare bone strength, elastic properties and mineral composition in Longshanks mice and randomly bred controls. Our data show that, despite being more slender, cortical geometry and predicted bending strength of the Longshanks tibia were similar to controls. In whole bone bending tests, measures of bone bending strength were similar across groups; however, Longshanks tibiae were significantly more rigid, more brittle, and required less than half the energy to fracture. Tissue-level elastic properties were also altered in Longshanks mice, but the bones did not differ from the control in water content, ash content or density. These results indicate that while Longshanks bones are as strong as control tibiae, selection for increased tibia length has altered its elastic properties, possibly through changes in organic bony matrix composition. We conclude that selection for certain limb morphologies, and/or selection for rapid skeletal growth, can lead to tissue-level changes that can increase the risk of skeletal fracture, which in turn may favor the correlated evolution of compensatory mechanisms to mitigate increased fracture risk, such as delayed skeletal maturity.

The genetic architecture of teosinte catalyzed and constrained maize domestication

Crop domestication is a well-established system for understanding evolution. We interrogated the genetic architecture of maize domestication from a quantitative genetics perspective. We analyzed domestication-related traits in a maize landrace and a population of its ancestor, teosinte. We observed strong divergence in the underlying genetic architecture including change in the genetic correlations among traits. Despite striking divergence, selection intensities were low for all traits, indicating that selection under domestication can be weaker than natural selection. Analyses suggest total grain weight per plant was not improved and that genetic correlations placed considerable constraint on selection. We hope our results will motivate crop evolutionists to perform similar work in other crops.

The process of evolution under domestication has been studied using phylogenetics, population genetics–genomics, quantitative trait locus (QTL) mapping, gene expression assays, and archaeology. Here, we apply an evolutionary quantitative genetic approach to understand the constraints imposed by the genetic architecture of trait variation in teosinte, the wild ancestor of maize, and the consequences of domestication on genetic architecture. Using modern teosinte and maize landrace populations as proxies for the ancestor and domesticate, respectively, we estimated heritabilities, additive and dominance genetic variances, genetic-by-environment variances, genetic correlations, and genetic covariances for 18 domestication-related traits using realized genomic relationships estimated from genome-wide markers. We found a reduction in heritabilities across most traits, and the reduction is stronger in reproductive traits (size and numbers of grains and ears) than vegetative traits. We observed larger depletion in additive genetic variance than dominance genetic variance. Selection intensities during domestication were weak for all traits, with reproductive traits showing the highest values. For 17 of 18 traits, neutral divergence is rejected, suggesting they were targets of selection during domestication. Yield (total grain weight) per plant is the sole trait that selection does not appear to have improved in maize relative to teosinte. From a multivariate evolution perspective, we identified a strong, nonneutral divergence between teosinte and maize landrace genetic variance–covariance matrices (G-matrices). While the structure of G-matrix in teosinte posed considerable genetic constraint on early domestication, the maize landrace G-matrix indicates that the degree of constraint is more unfavorable for further evolution along the same trajectory.

Education, Smoking, and Cohort Change: Forwarding a Multidimensional Theory of the Environmental Moderation of Genetic Effects

We introduce a genetic correlation by environment interaction model [(rG)xE] which allows for social environmental moderation of the genetic relationship between two traits. To empirically demonstrate the significance of the (rG)xE perspective, we use genome wide information from respondents of the Health and Retirement Study (HRS; n = 8,181; birth years 1920-1959) and the National Longitudinal Study of Adolescent to Adult Health (Add Health; n = 4,347; birth years 1974-1983) to examine whether the genetic correlation (rG) between education and smoking has increased over historical time. Genetic correlation estimates (rG_HRS = -0.357; rG_{Add Health} = -0.729) support this hypothesis. Using polygenic scores for educational attainment, we show that this is not due to latent indicators of intellectual capacity, and we highlight the importance of education itself as an explanation of the increasing genetic correlation. Analyses based on contextual variation the milieus of the Add Health respondents corroborate key elements of the birth cohort analyses. We argue that the increasing overlap with respect to genes associated with educational attainment and smoking is a fundamentally social process involving complex process of selection based on observable behaviors that may be linked to genotype.

Artificial selection sheds light on developmental mechanisms of limb elongation

Species diversity in limb lengths and proportions is thought to have evolved adaptively in the context of locomotor and habitat specialization, but the heritable cellular processes that drove this evolution within species are poorly understood. In this study, we take a novel "micro-evo-devo" approach, using artificial selection on relative limb length to amplify phenotypic variation in a population of mice, known as Longshanks, to examine the cellular mechanisms of postnatal limb development that contribute to intraspecific limb length variation. Cross-sectional growth data indicate that differences in bone length between Longshanks and random-bred controls are not due to prolonged growth, but to accelerated growth rates. Histomorphometric and cell proliferation assays on proximal tibial growth plates show that Longshanks' increased limb bone length is associated with an increased number of proliferative chondrocytes. In contrast, we find no differences in other growth plate cellular features known to underlie interspecific differences in limb bone size and shape, such as the rates of chondrocyte proliferation or the size and number of hypertrophic cells in the growth plate. These data suggest that small differences among individuals in the number of proliferating chondrocytes are a potentially important determinant of selectable intraspecific variation in individual limb bone lengths, independent of body size.

Cortical and trabecular morphology is altered in the limb bones of mice artificially selected for faster skeletal growth

Bone strength is influenced by mineral density and macro- and microstructure. Research into factors that contribute to bone morphology and strength has focused on genetic, environmental and morphological factors (e.g., body mass index), but little is known regarding the impact of rates of skeletal elongation on adult skeletal morphology and strength. Using micro-CT, we examined the impact of rates of skeletal elongation on bone cortical and trabecular morphology, and on rates of estrogen-dependent bone loss in the tibia in CD-1 mice, and in mice with accelerated skeletal growth (Longshanks). Groups of adult mice (n = 7/group) were subjected to ovariectomy or sham surgeries, scanned for 6 weeks, and indices of bone morphology were collected. Results show that Longshanks mice had significantly less trabecular bone at skeletal maturity, characterized by fewer, thinner trabeculae, and furthermore lost trabecular bone more slowly in response to ovariectomy. Artificial selection for rapid skeletal growth relative to somatic growth thus had a significant impact on trabecular bone morphology in Longshanks. Our data do not unequivocally demonstrate a causal relationship between rapid bone growth and reduced trabecular bone quality, but suggest that rapid linear bone growth may influence the risk of cancellous bone fragility.

Gait changes in a line of mice artificially selected for longer limbs

In legged terrestrial locomotion, the duration of stance phase, i.e., when limbs are in contact with the substrate, is positively correlated with limb length, and negatively correlated with the metabolic cost of transport. These relationships are well documented at the interspecific level, across a broad range of body sizes and travel speeds. However, such relationships are harder to evaluate within species (i.e., where natural selection operates), largely for practical reasons, including low population variance in limb length, and the presence of confounding factors such as body mass, or training. Here, we compared spatiotemporal kinematics of gait in Longshanks, a long-legged mouse line created through artificial selection, and in random-bred, mass-matched Control mice raised under identical conditions. We used a gait treadmill to test the hypothesis that Longshanks have longer stance phases and stride lengths, and decreased stride frequencies in both fore- and hind limbs, compared with Controls. Our results indicate that gait differs significantly between the two groups. Specifically, and as hypothesized, stance duration and stride length are 8–10% greater in Longshanks, while stride frequency is 8% lower than in Controls. However, there was no difference in the touch-down timing and sequence of the paws between the two lines. Taken together, these data suggest that, for a given speed, Longshanks mice take significantly fewer, longer steps to cover the same distance or running time compared to Controls, with important implications for other measures of variation among individuals in whole-organism performance, such as the metabolic cost of transport.

How does epistasis influence the response to selection?

Much of quantitative genetics is based on the 'infinitesimal model', under which selection has a negligible effect on the genetic variance. This is typically justified by assuming a very large number of loci with additive effects. However, it applies even when genes interact, provided that the number of loci is large enough that selection on each of them is weak relative to random drift. In the long term, directional selection will change allele frequencies, but even then, the effects of epistasis on the ultimate change in trait mean due to selection may be modest. Stabilising selection can maintain many traits close to their optima, even when the underlying alleles are weakly selected. However, the number of traits that can be optimised is apparently limited to ~4Ne by the 'drift load', and this is hard to reconcile with the apparent complexity of many organisms. Just as for the mutation load, this limit can be evaded by a particular form of negative epistasis. A more robust limit is set by the variance in reproductive success. This suggests that selection accumulates information most efficiently in the infinitesimal regime, when selection on individual alleles is weak, and comparable with random drift. A review of evidence on selection strength suggests that although most variance in fitness may be because of alleles with large Nes, substantial amounts of adaptation may be because of alleles in the infinitesimal regime, in which epistasis has modest effects.

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

Integrating Evolutionary Game Theory into Mechanistic Genotype-Phenotype Mapping

Natural selection has shaped the evolution of organisms toward optimizing their structural and functional design. However, how this universal principle can enhance genotype-phenotype mapping of quantitative traits has remained unexplored. Here we show that the integration of this principle and functional mapping through evolutionary game theory gains new insight into the genetic architecture of complex traits. By viewing phenotype formation as an evolutionary system, we formulate mathematical equations to model the ecological mechanisms that drive the interaction and coordination of its constituent components toward population dynamics and stability. Functional mapping provides a procedure for estimating the genetic parameters that specify the dynamic relationship of competition and cooperation and predicting how genes mediate the evolution of this relationship during trait formation.

Changes in shape and cross-sectional geometry in the tibia of mice selectively bred for increases in relative bone length

Limb bone size and shape in terrestrial mammals scales predictably with body mass. Weight-bearing limb bones in these species have geometries that enable them to withstand deformations due to loading, both within and between species. Departures from the expected scaling of bone size and shape to body mass occur in mammals that have become specialized for different types of locomotion. For example, mammals adapted for frequent running and jumping behaviors have hind limb bones that are long in relation to body mass, but with narrower cross-sections than predicted for their length. The Longshanks mouse was recently established, a selectively bred line of mice with ~12-13% longer tibiae relative to body mass. This increased limb length resembles superficially the derived limb proportions of rodents adapted for hopping and jumping. Here, 3D geometric morphometrics and analyses of bone cross-sectional geometry were combined to determine whether selection for increased relative tibia length in Longshanks mice has altered the scaling relationship of size and shape, and/or bone robusticity, relative to the tibiae of random-bred control mice from the same genetic background. The results suggest that the Longshanks tibia is not a geometrically scaled version of the control tibiae. Instead, the Longshanks tibia has become narrower in cross-section in relation to its increased length, leading to a decrease in overall bending strength when compared with control tibiae. These changes in bone shape and robusticity resemble the derived morphology of mammals adapted for running and jumping, with important implications for the material properties and strength of bone in these mammals.