Search

Summary Objective To investigate the transcriptomic differences in chondrocytes obtained from LG/J (large, healer) and SM/J (small, non-healer) murine strains in an attempt to discern the molecular pathways implicated in cartilage regeneration and susceptibility to osteoarthritis (OA). Design We performed RNA-sequencing on chondrocytes derived from LG/J (n = 16) and SM/J (n = 16) mice. We validated the expression of candidate genes and compared single nucleotide polymorphisms (SNPs) between the two mouse strains. We also examined gene expression of positional candidates for ear pinna regeneration and long bone length quantitative trait loci (QTLs) that display differences in cartilaginous expression. Results We observed a distinct genetic heterogeneity between cells derived from LG/J and SM/J mouse strains. We found that gene ontologies representing cell development, cartilage condensation, and regulation of cell differentiation were enriched in LG/J chondrocytes. In contrast, gene ontologies enriched in the SM/J chondrocytes were mainly related to inflammation and degeneration. Moreover, SNP analysis revealed that multiple validated genes vary in sequence between LG/J and SM/J in coding and highly conserved noncoding regions. Finally, we showed that most QTLs have 20–30% of their positional candidates displaying differential expression between the two mouse strains. Conclusions While the enrichment of pathways related to cell differentiation, cartilage development and cartilage condensation infers superior healing potential of LG/J strain, the enrichment of pathways related to cytokine production, immune cell activation and inflammation entails greater susceptibility of SM/J strain to OA. These data provide novel insights into chondrocyte transcriptome and aid in identification of the quantitative trait genes and molecular differences underlying the phenotypic differences associated with individual QTLs.

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Parent-of-origin effects are unexpectedly common in complex traits, including metabolic and neurological diseases. Parent-of-origin effects can be modified by the environment, but the architecture of these gene-by-environmental effects on phenotypes remains to be unraveled. Previously, quantitative trait loci (QTL) showing context-specific parent-of-origin effects on metabolic traits were mapped in the F16generation of an advanced intercross between LG/J and SM/J inbred mice. However, these QTL were not enriched for known imprinted genes, suggesting another mechanism is needed to explain these parent-of-origin effects phenomena. We propose that non-imprinted genes can generate complex parent-of-origin effects on metabolic traits through interactions with imprinted genes. Here, we employ data from mouse populations at different levels of intercrossing (F0, F1, F2, F16) of the LG/J and SM/J inbred mouse lines to test this hypothesis. Using multiple populations and incorporating genetic, genomic, and physiological data, we leverage orthogonal evidence to identify networks of genes through which parent-of-origin effects propagate. We identify a network comprised of 3 imprinted and 6 non-imprinted genes that show parent-of-origin effects. This epistatic network forms a nutritional responsive pathway and the genes comprising it jointly serve cellular functions associated with growth. We focus on 2 genes,NnatandF2r, whose interaction associates with serum glucose levels across generations in high fat-fed females. Single-cell RNAseq reveals thatNnatandF2rare negatively correlated in pre-adipocytes along an adipogenic trajectory, a result that is consistent with our observations in bulk white adipose tissue.

Objective Recombinant inbred mouse strains generated from an LG/J and SM/J intercross offer a unique resource to study complex genetic traits such as osteoarthritis (OA). We undertook this study to determine the susceptibility of 14 strains to various phenotypes characteristic of posttraumatic OA. We hypothesized that phenotypic variability is associated with genetic variability. Methods Ten-week-old male mice underwent surgical destabilization of the medial meniscus (DMM) to induce posttraumatic OA. Mice were killed 8 weeks after surgery, and knee joints were processed for histology to score cartilage degeneration and synovitis. Micro-computed tomography was used to analyze trabecular bone parameters including subchondral bone plate thickness and synovial ectopic calcifications. Gene expression in the knees was assessed using a QuantiGene Plex assay. Results Broad-sense heritability ranged from 0.18 to 0.58, which suggested that the responses to surgery were moderately heritable. The LGXSM-33, LGXSM-5, LGXSM-46, and SM/J strains were highly susceptible to OA, while the LGXSM-131b, LGXSM-163, LGXSM-35, LGXSM-128a, LGXSM-6, and LG/J strains were relatively OA resistant. This study was the first to accomplish measurement of genetic correlations of phenotypes that are characteristic of posttraumatic OA. Cartilage degeneration was significantly positively associated with synovitis (r = 0.83-0.92), and subchondral bone plate thickness was negatively correlated with ectopic calcifications (r = -0.59). Moreover, we showed that 40 of the 78 genes tested were significantly correlated with various OA phenotypes. However, unlike the OA phenotypes, there was no evidence for genetic variation in differences in gene expression levels between DMM-operated and sham-operated knees. Conclusion For these mouse strains, various characteristics of posttraumatic OA varied with genetic composition, which demonstrated a genetic basis for susceptibility to posttraumatic OA. The heritability of posttraumatic OA was established. Phenotypes exhibited various degrees of correlations; cartilage degeneration was positively correlated with synovitis, but not with the formation of ectopic calcifications. Further investigation of the genome regions that contain genes implicated in OA, as well as further investigation of gene expression data, will be useful for studying mechanisms of OA and identifying therapeutic targets.

Obesity is generally protective against osteoporosis and bone fracture. However, recent studies indicate that the influence of obesity on the skeleton is complex and can be detrimental. We evaluated the effects of a high-fat, obesogenic diet on the femur and radius of 1100 mice (males and females) from the Large-by-Small advanced intercross line (F34 generation). At age 5 months, bone morphology was assessed by microCT and mechanical properties by three-point bending. Mice raised on a high-fat diet had modestly greater cortical area, bending stiffness, and strength. Size-independent material properties were unaffected by a high-fat diet, indicating that diet influenced bone quantity but not quality. Bone size and mechanical properties were strongly correlated with body mass. However, the increases in many bone traits per unit increase in body mass were less in high-fat diet mice than low-fat diet mice. Thus, although mice raised on a high-fat diet have, on average, bigger and stronger bones than low-fat-fed mice, a high-fat diet diminished the positive relationship between body mass and bone size and whole-bone strength. The findings support the concept that there are diminishing benefits to skeletal health with increasing obesity. © 2019 American Society for Bone and Mineral Research.

Intervertebral disc degeneration (IDD) causes back pain and sciatica, affecting quality of life and resulting in high economic/social burden. The etiology of IDD is not well understood. Along with aging and environmental factors, genetic factors also influence the onset, progression and severity of IDD. Genetic studies of risk factors for IDD using human cohorts are limited by small sample size and low statistical power. Animal models amenable to genetic and functional studies of IDD provide desirable alternatives. Despite differences in size and cellular content as compared to human intervertebral discs (IVDs), the mouse is a powerful model for genetics and assessment of cellular changes relevant to human biology. Here, we provide evidence for early onset disc degeneration in SM/J relative to LG/J mice with poor and good tissue healing capacity respectively. In the first few months of life, LG/J mice maintain a relatively constant pool of notochordal-like cells in the nucleus pulposus (NP) of the IVD. In contrast, chondrogenic events are observed in SM/J mice beginning as early as one-week-old, with progressive fibrotic changes. Further, the extracellular matrix changes in the NP are consistent with IVD degeneration. Leveraging on the genomic data of two parental and two recombinant inbred lines, we assessed the genetic contribution to the NP changes and identified processes linked to the regulation of ion transport systems. Significantly, "transport" system is also in the top three gene ontology (GO) terms from a comparative proteomic analysis of the mouse NP. These findings support the potential of the SM/J, LG/J and their recombinant inbred lines for future genetic and biological analysis in mice and validation of candidate genes and biological relevance in human cohort studies. The proteomic data has been deposited to the ProteomeXchange Consortium via the PRIDE [1] partner repository with the dataset identifier PXD008784.

Background Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression. Results Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape. Conclusions Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

The biology of how faces are built and come to differ from one another is complex. Discovering the genes that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL’s peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.Author SummaryThe mammalian face is a complex structure serving many functions. We studied the genetic basis for facial skeletal diversity in a large sample of mice from an experimental population designed for the study of complex traits. We quantified the contribution of genetic variation to variation in three-dimensional facial shape across more than 55,000 genetic markers spread throughout the mouse genome. We found 37 genetic regions which are very likely to contribute to differences in facial shape. We then conducted a more detailed analysis of the genetic regions associated with the most variable aspects of facial shape. For these regions, a disproportionately large number of genes are known to be important to growth and to skeletal and facial development. The magnitude of these genetic contributions to differences in facial shape are consistently small. Our results therefore support the notion that facial skeletal diversity is influenced by many genes of small effect, but that some of these small effects may be related to genes that are fundamental to skeletal and facial development.

Summary Objective The etiology of osteoarthritis (OA) is unknown, however, there appears to be a significant contribution from genetics. We have identified recombinant inbred strains of mice derived from LG/J (large) and SM/J (small) strains that vary significantly in their ability to repair articular cartilage and susceptibility to post-traumatic OA due to their genetic composition. Here, we report cartilage repair phenotypes in the same strains of mice in which OA susceptibility was analyzed previously, and determine the genetic correlations between phenotypes. Design We used 12 recombinant inbred strains, including the parental strains, to test three phenotypes: ear-wound healing (n = 263), knee articular cartilage repair (n = 131), and post-traumatic OA (n = 53) induced by the surgical destabilization of the medial meniscus (DMM). Genetic correlations between various traits were calculated as Pearson's correlation coefficients of strain means. Results We found a significant positive correlation between ear-wound healing and articular cartilage regeneration (r = 0.71; P = 0.005). We observed a strong inverse correlation between articular cartilage regeneration and susceptibility to OA based on maximum (r = −0.54; P = 0.036) and summed Osteoarthritis Research Society International (OARSI) scores (r = −0.56; P = 0.028). Synovitis was not significantly correlated with articular cartilage regeneration but was significantly positively correlated with maximum (r = 0.63; P = 0.014) and summed (r = 0.70; P = 0.005) OARSI scores. Ectopic calcification was significantly positively correlated with articular cartilage regeneration (r = 0.59; P = 0.021). Conclusions Using recombinant inbred strains, our study allows, for the first time, the measurement of genetic correlations of regeneration phenotypes with degeneration phenotypes, characteristic of OA (cartilage degeneration, synovitis). We demonstrate that OA is positively correlated with synovitis and inversely correlated with the ability to repair cartilage. These results suggest an addition to the risk paradigm for OA from a focus on degeneration to regeneration.

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high‐resolution quantitative trait locus (QTL) mapping in the F34 generation of an advanced intercross experiment (LG,SM‐G34) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome‐ or chromosome‐wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation‐seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Parent-of-origin effects (POE) are unexpectedly common in complex traits, including metabolic and neurological diseases. POE can also be modified by the environment, but the architecture of these gene-by-environmental effects on phenotypes remains to be unraveled. Previously, quantitative trait loci (QTL) showing context-specific POE on metabolic traits were mapped in the F16 generation of an advanced intercross between LG/J and SM/J inbred mice. However, these QTL were not enriched for known imprinted genes, suggesting another mechanism is needed to explain these POE phenomena. Here, we use a simple yet powerful F1 reciprocal cross model to test the hypothesis that non-imprinted genes can generate complex POE on metabolic traits through genetic interactions with imprinted genes. Male and female mice from a F1 reciprocal cross of LG/J and SM/J strains were fed either high or low fat diets. We generated expression profiles from three metabolically-relevant tissues: hypothalamus, white adipose, and liver. We identified two classes of parent-of-origin expression biases: genes showing parent-of-origin-dependent allele-specific expression and biallelic genes that are differentially expressed by reciprocal cross. POE patterns of both gene classes are highly tissue-and context-specific, sometimes occurring only in one sex and/or diet cohort in a particular tissue. We then constructed tissue-specific interaction networks among genes from these two classes of POE. A key subset of gene pairs show significant epistasis in the F16 LG/J x SM/J advanced intercross data in cases where the biallelic gene fell within a previously-identified metabolic POE QTL interval. We highlight one such interaction in adipose, between Nnat and Mogat1, which associates with POE on multiple adiposity traits. Both genes localize to the endoplasmic reticulum of adipocytes and play a role in adipogenesis. Additionally, expression of both genes is significantly correlated in human visceral adipose tissue. The genes and networks we present here represent a set of actionable interacting candidates that can be probed to further identify the machinery driving POE on complex traits.

BackgroundIn newborns, severe congenital heart defects are rarer than mild ones. The reason why is unknown, but presumably related to a liability threshold that rises with the severity of a defect. Because the same genetic mutation can cause different defects, other variables may contribute to pushing an individual past a defect-specific liability threshold. We consider here how variables in the genetic architecture of a heart defect depend upon its fitness cost, as defined by the likelihood of survival to reproductive age in natural history studies.MethodsWe phenotyped ~10,000Nkx2-5+/-newborn mice, a model of human congenital heart disease, from two inbred strain crosses. Genome-wide association analyses detected loci that modify the risk of an atrial septal defect, membranous or muscular ventricular septal defect, or atrioventricular septal defect. The number of loci, heritability and quantitative effects on risk of pairwise (G×GNkx) and higher-order (G×G×GNkx) epistasis between the loci andNkx2-5mutation were examined as a function of the fitness cost of a defect.ResultsNkx2-5+/-mice have pleiotropic heart defects; about 70% have normal hearts. The model recapitulates the epidemiological relationship between the severity and incidence of a heart defect. Neither the number of modifier loci nor heritability depends upon the severity of a defect, but G×GNkxand G×G×GNkxeffects on risk do. Interestingly, G×G×GNkxeffects are three times more likely to suppress risk when the genotypes at the first two loci are homozygous and from the same, rather than opposite strains in a cross. Syn- and anti-homozygous genotypes at G×G×GNkxinteractions can have an especially large impact on the risk of an atrioventricular septal defect.ConclusionsGiven a modestly penetrant mutation, epistasis contributes more to the risk of severe than mild congenital heart defect. Conversely, genetic compatibility between interacting genes, as indicated by the protective effects of syn-homozygosity at G×G×GNkxinteractions, plays a newfound role in the robustness of cardiac development. The experimental model offers practical insights into the nature of genetic risk in congenital heart disease. The results more fundamentally address a longstanding question regarding how mutational robustness could arise from natural selection.

Evolutionary studies have long emphasized that the genetic architecture of traits holds important microevolutionary consequences. Yet, studies comparing the genetic architecture of traits across species are rare, and discussions of the evolution of genetic systems are made on theoretical arguments rather than on empirical evidence. Here, we compared the genetic architecture of cranial traits in two different mammalian model organisms: the gray short-tailed opossum, Monodelphis domestica, and the laboratory mouse, Mus musculus. We show that both organisms share a highly polygenic genetic architecture for craniofacial traits, with many loci of small effect. However, these two model species differ significantly in the overall degree of pleiotropy, N, of the genotype-to-phenotype map, with opossums presenting a higher average N. They also diverge in their degree of genetic modularity, with opossums presenting less modular patterns of genetic association among traits. We argue that such differences highlight the context dependency of gene effects, with developmental systems shaping the variational properties of genetic systems. Finally, we also demonstrate based on the opossum data that current measurements for the relationship between the mutational effect size and N need to be re-evaluated in relation to the importance of the cost of pleiotropy for mammals.

Modularity has emerged as a central concept for evolutionary biology, thereby providing the field with a theory of organismal structure and variation. This theory has reframed long-standing questions and serves as a unified conceptual framework for genetics, developmental biology, and multivariate evolution. Research programs in systems biology and quantitative genetics are bridging the gap between these fields. Although this synthesis is ongoing, some major themes have emerged, and empirical evidence for modularity has become abundant. In this review, we look at modularity from a historical perspective, highlighting its meaning at different levels of biological organization and the different methods that can be used to detect it. We then explore the relationship between quantitative genetic approaches to modularity and developmental genetic studies. We conclude by investigating the dynamic relationship between modularity and the adaptive landscape and how this relationship potentially shapes evolution and can help bridge the gap between micro- and macroevolution.

Despite decades of progress, congenital heart disease remains a major cause of mortality and suffering in children and young adults. Prevention would be ideal, but formidable biological and technical hurdles face any intervention that seeks to target the main causes, genetic mutations in the embryo. Other factors, however, significantly modify the total risk in individuals who carry mutations. Investigation of these factors could lead to an alternative approach to prevention. To define the risk modifiers, our group has taken an "experimental epidemiologic" approach via inbred mouse strain crosses. The original intent was to map genes that modify an individual's risk of heart defects caused by an Nkx2-5 mutation. During the analysis of >2000 Nkx2-5(+/-) offspring from one cross we serendipitously discovered a maternal-age associated risk, which also exists in humans. Reciprocal ovarian transplants between young and old mothers indicate that the incidence of heart defects correlates with the age of the mother and not the oocyte, which implicates a maternal pathway as the basis of the risk. The quantitative risk varies between strain backgrounds, so maternal genetic polymorphisms determine the activity of a factor or factors in the pathway. Most strikingly, voluntary exercise by the mother mitigates the risk. Therefore, congenital heart disease can in principle be prevented by targeting a maternal pathway even if the embryo carries a causative mutation. Further mechanistic insight is necessary to develop an intervention that could be implemented on a broad scale, but the physiology of maternal-fetal interactions, aging, and exercise are notoriously complex and undefined. This suggests that an unbiased genetic approach would most efficiently lead to the relevant pathway. A genetic foundation would lay the groundwork for human studies and clinical trials.

The rapid divergence of male genitalia is a preeminent evolutionary pattern. This rapid divergence is especially striking in the baculum, a bone that occurs in the penis of many mammalian species. Closely related species often display diverse baculum morphology where no other morphological differences can be discerned. While this fundamental pattern of evolution has been appreciated at the level of gross morphology, nearly nothing is known about the genetic basis of size and shape divergence. Quantifying the genetic basis of baculum size and shape variation has been difficult because these structures generally lack obvious landmarks, so comparing them in three dimensions is not straightforward. Here, we develop a novel morphometric approach to quantify size and shape variation from three-dimensional micro-CT scans taken from 369 bacula, representing 75 distinct strains of the BXD family of mice. We identify two quantitative trait loci (QTL) that explain ∼50% of the variance in baculum size, and a third QTL that explains more than 20% of the variance in shape. Together, our study demonstrates that baculum morphology may diverge relatively easily, with mutations at a few loci of large effect that independently modulate size and shape. Based on a combination of bioinformatic investigations and new data on RNA expression, we prioritized these QTL to 16 candidate genes, which have hypothesized roles in bone morphogenesis and may enable future genetic manipulation of baculum morphology.

With the advent of modern imaging and measurement technology, complex phenotypes are increasingly represented by large numbers of measurements, which may not bear biological meaning one by one. For such multivariate phenotypes, studying the pairwise associations between all measurements and all alleles is highly inefficient and prevents insight into the genetic pattern underlying the observed phenotypes. We present a new method for identifying patterns of allelic variation (genetic latent variables) that are maximally associated—in terms of effect size—with patterns of phenotypic variation (phenotypic latent variables). This multivariate genotype–phenotype mapping (MGP) separates phenotypic features under strong genetic control from less genetically determined features and thus permits an analysis of the multivariate structure of genotype–phenotype association, including its dimensionality and the clustering of genetic and phenotypic variables within this association. Different variants of MGP maximize different measures of genotype–phenotype association: genetic effect, genetic variance, or heritability. In an application to a mouse sample, scored for 353 SNPs and 11 phenotypic traits, the first dimension of genetic and phenotypic latent variables accounted for >70% of genetic variation present in all 11 measurements; 43% of variation in this phenotypic pattern was explained by the corresponding genetic latent variable. The first three dimensions together sufficed to account for almost 90% of genetic variation in the measurements and for all the interpretable genotype–phenotype association. Each dimension can be tested as a whole against the hypothesis of no association, thereby reducing the number of statistical tests from 7766 to 3—the maximal number of meaningful independent tests. Important alleles can be selected based on their effect size (additive or nonadditive effect on the phenotypic latent variable). This low dimensionality of the genotype–phenotype map has important consequences for gene identification and may shed light on the evolvability of organisms.

Bilateral symmetry is a fundamental property of the vertebrate central nervous system. Local deviations from symmetry provide various types of information about the development, evolution, and function of elements within the CNS, especially the cerebral hemispheres. Here, we quantify the pattern and extent of asymmetry in cortical folding within the cerebrum of Papio baboons and assess the evolutionary and developmental implications of the findings. Analyses of directional asymmetry show a population-level trend in length measurements indicating that baboons are genetically predisposed to be asymmetrical, with the right side longer than the left in the anterior cerebrum while the left side is longer than the right posteriorly. We also find a corresponding bias to display a right frontal petalia (overgrowth of the anterior pole of the cerebral cortex on the right side). By quantifying fluctuating asymmetry, we assess canalization of brain features and the susceptibility of the baboon brain to developmental perturbations. We find that features are differentially canalized depending on their ontogenetic timing. We further deduce that development of the two hemispheres is to some degree independent. This independence has important implications for the evolution of cerebral hemispheres and their separate specialization. Asymmetry is a major feature of primate brains and is characteristic of both brain structure and function.

Analysis of wing variation within and between fly species reveals an unexpectedly slow evolutionary rate. Variations due to mutation and interspecific differences are similar, perhaps as a result of complex genetic interactions. See Letter p.447 Genetic mutation is the engine room of evolution. Without genetic mutation, genes would not vary, and natural selection would be powerless to act. But if genetic mutation provides the raw material for adaption, the extent to which mutation further constrains evolution over geological time is a matter of debate. Mutation rates are typically very small and are thought to be swamped by the amount of genetic divergence seen over evolutionary timescales. David Houle and colleagues present results that could be a fly in the ointment. Several thousand flies in fact. They examined the wings of more than 50,000 Drosophila and its close relatives to reveal an unexpected and strong positive relationship between mutation, genetic variation and evolutionary rate over the past 40 million years. This challenges several models of the role of mutation in evolution.

Background While the genetics of obesity has been well defined, the epigenetics of obesity is poorly understood. Here, we used a genome-wide approach to identify genes with differences in both DNA methylation and expression associated with a high-fat diet in mice. Results We weaned genetically identical Small (SM/J) mice onto a high-fat or low-fat diet and measured their weights weekly, tested their glucose and insulin tolerance, assessed serum biomarkers, and weighed their organs at necropsy. We measured liver gene expression with RNA-seq (using 21 total libraries, each pooled with 2 mice of the same sex and diet) and DNA methylation with MRE-seq and MeDIP-seq (using 8 total libraries, each pooled with 4 mice of the same sex and diet). There were 4356 genes with expression differences associated with diet, with 184 genes exhibiting a sex-by-diet interaction. Dietary fat dysregulated several pathways, including those involved in cytokine-cytokine receptor interaction, chemokine signaling, and oxidative phosphorylation. Over 7000 genes had differentially methylated regions associated with diet, which occurred in regulatory regions more often than expected by chance. Only 5–10% of differentially methylated regions occurred in differentially expressed genes, however this was more often than expected by chance (p = 2.2 × 10− 8). Conclusions Discovering the gene expression and methylation changes associated with a high-fat diet can help to identify new targets for epigenetic therapies and inform about the physiological changes in obesity. Here, we identified numerous genes with altered expression and methylation that are promising candidates for further study. Electronic supplementary material The online version of this article (10.1186/s12864-018-5327-0) contains supplementary material, which is available to authorized users.

Quantitative genetical theory suggests that characters existing in developmentally or functionally integrated groups are expected to be genetically correlated because they share a common inheritance. The genetic independence model for the cause of this genetic integration predicts that pleiotropic effects of single genes are mostly restricted to the characters in these integrated groups. We tested this model by estimating the additive and dominance effects of quantitative trait loci (QTLs) affecting early‐ (cranial vault) and late‐developing (face) skeletal characters in F2 house mice originally derived from a cross of the Large and Small inbred strains. Interval mapping procedures were used that resulted in the identification of 26 QTLs on 17 of the 19 autosomes that significantly affected these characters. Additive, but not dominance, genotypic effects of many of these QTLs predominantly affected either the cranial vault or face characters, which supports the genetic independence model. Only t...

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this question is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will depend on conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Evolutionary constraint due to pleiotropy refers to a situation in which mutations in genes shared among traits generate trait covariance; therefore, traits that are not directly exposed to selective challenge show a correlated response. When such a correlated response is deleterious, it may constrain the trait from evolving. Here, we argue that the idea of absolute constraints draws from the perception that gene effects are inherent to alleles and thus invariant across genetic and environmental backgrounds. However, evidence from studies involving genetic effects on multiple traits, observed across different genetic backgrounds and environments, supports the notion that genes' effects on traits change. Consequently, pleiotropy also varies across backgrounds. We argue for a stronger emphasis on interaction effects when describing a trait's genetic basis and its evolutionary potential. By discussing different cases of trait individuation, we demonstrate how this approach can lead to new insights.

This study reports on genetic susceptibility to ectopic calcification in the LG/J and SM/J advanced intercross mice. Using 347 mice in 98 full-sibships, destabilization of medial meniscus (DMM) was performed to induce joint injury. We found that joint destabilization instigated ectopic calcifications as detected and quantified by micro-CT. We performed quantitative trait locus (QTL) analysis to map ectopic calcification phenotypes to discrete genomic locations. To validate the functional significance of the selected QTL candidate genes, we compared mRNA expression between parental LG/J and SM/J inbred strains. Overall, we detected 20 QTLs affecting synovial and meniscus calcification phenotypes with 11 QTLs linked to synovial calcification. Functional and bioinformatic analyses of single nucleotide polymorphism (SNP) identified functional classifications relevant to angiogenesis (Myo1e, Kif26b, Nprl3, Stab2, Fam105b), bone metabolism/calcification (Tle3, Tgfb2, Lipc, Nfe2l1, Ank, Fam105b), arthritis (Stab2, Tbx21, Map4k4, Hoxb9, Larp6, Col1a2, Adam10, Timp3, Nfe2l1, Trpm3), and ankylosing-spondylitis (Ank, Pon1, Il1r2, Tbkbp1) indicating that ectopic calcification involves multiple mechanisms. Furthermore, the expression of 11 out of 78 candidate genes was significantly different between LG/J and SM/J. Correlation analysis showed that Aff3, Fam81a, Syn3, and Ank were correlated with synovial calcification. Taken together, our findings of multiple genetic loci suggest the involvement of multiple genes contributing to ectopic calcification.

Folding of the primate brain cortex allows for improved neural processing power by increasing cortical surface area for the allocation of neurons. The arrangement of folds (sulci) and ridges (gyri) across the cerebral cortex is thought to reflect the underlying neural network. Gyrification, an adaptive trait with a unique evolutionary history, is affected by genetic factors different from those affecting brain volume. Using a large pedigreed population of ∼1000 Papio baboons, we address critical questions about the genetic architecture of primate brain folding, the interplay between genetics, brain anatomy, development, patterns of cortical–cortical connectivity, and gyrification’s potential for future evolution. Through Mantel testing and cluster analyses, we find that the baboon cortex is quite evolvable, with high integration between the genotype and phenotype. We further find significantly similar partitioning of variation between cortical development, anatomy, and connectivity, supporting the predictions of tension-based models for sulcal development. We identify a significant, moderate degree of genetic control over variation in sulcal length, with gyrus-shape features being more susceptible to environmental effects. Finally, through QTL mapping, we identify novel chromosomal regions affecting variation in brain folding. The most significant QTL contain compelling candidate genes, including gene clusters associated with Williams and Down syndromes. The QTL distribution suggests a complex genetic architecture for gyrification with both polygeny and pleiotropy. Our results provide a solid preliminary characterization of the genetic basis of primate brain folding, a unique and biomedically relevant phenotype with significant implications in primate brain evolution.

We tested the hypothesis that the rate of marsupial cranial evolution is dependent on the distribution of genetic variation in multivariate space. To do so, we carried out a genetic analysis of cranial morphological variation in laboratory strains of Monodelphis domestica and used estimates of genetic covariation to analyze the morphological diversification of the Monodelphis brevicaudata species group. We found that within-species genetic variation is concentrated in only a few axes of the morphospace and that this strong genetic covariation influenced the rate of morphological diversification of the brevicaudata group, with between-species divergence occurring fastest when occurring along the genetic line of least resistance. Accounting for the geometric distribution of genetic variation also increased our ability to detect the selective regimen underlying species diversification, with several instances of selection only being detected when genetic covariances were taken into account. Therefore, this work directly links patterns of genetic covariation among traits to macroevolutionary patterns of morphological divergence. Our findings also suggest that the limited distribution of Monodelphis species in morphospace is the result of a complex interplay between the limited dimensionality of available genetic variation and strong stabilizing selection along two major axes of genetic variation.

Robustness to perturbation is a fundamental feature of complex organisms. Mutations are the raw material for evolution, yet robustness to their effects is required for species survival. The mechanisms that produce robustness are poorly understood. Nonlinearities are a ubiquitous feature of development that may link variation in development to phenotypic robustness. Here, we manipulate the gene dosage of a signaling molecule, Fgf8, a critical regulator of vertebrate development. We demonstrate that variation in Fgf8 expression has a nonlinear relationship to phenotypic variation, predicting levels of robustness among genotypes. Differences in robustness are not due to gene expression variance or dysregulation, but emerge from the nonlinearity of the genotype–phenotype curve. In this instance, embedded features of development explain robustness differences. How such features vary in natural populations and relate to genetic variation are key questions for unraveling the origin and evolvability of this feature of organismal development., Developmental processes often involve nonlinearities, but the consequences for translating genotype to phenotype are not well characterized. Here, Green et al. vary Fgf8 signaling across allelic series of mice and show that phenotypic robustness in craniofacial shape is explained by a nonlinear effect of Fgf8 expression.

Objectives: Determining the genetic architecture of quantitative traits and genetic correlations among them is important for understanding morphological evolution patterns. We address two questions regarding papionin evolution: (1) what effect do body and cranial size, age, and sex have on phenotypic (VP) and additive genetic (VA) variation in baboon crania, and (2) how might additive genetic correlations between craniofacial traits and body mass affect morphological evolution? Materials and Methods: We use a large captive pedigreed baboon sample to estimate quantitative genetic parameters for craniofacial dimensions (EIDs). Our models include nested combinations of the covariates listed above. We also simulate the correlated response of a given EID due to selection on body mass alone. Results: Covariates account for 1.2-91% of craniofacial VP. EID VA decreases across models as more covariates are included. The median genetic correlation estimate between each EID and body mass is 0.33. Analysis of the multivariate response to selection reveals that observed patterns of craniofacial variation in extant baboons cannot be attributed solely to correlated response to selection on body mass, particularly in males. Discussion: Because a relatively large proportion of EID VA is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual VP patterns. This may conflate direct selection effects on craniofacial variation with those resulting from a correlated response to body mass selection. This shared genetic variation may partially explain how selection for increased body mass in two different papionin lineages produced remarkably similar craniofacial phenotypes.

Relationship loci (rQTL) exist when the correlation between multiple traits varies by genotype. rQTL often occur due to gene-by-gene (G × G) or gene-by-environmental interactions, making them a powerful tool for detecting G × G. Here we present an empirical analysis of apolipoprotein E (APOE) with respect to lipid traits and incident CHD leading to the discovery of loci that interact with APOE to affect these traits. We found that the relationship between total cholesterol (TC) and triglycerides (ln TG) varies by APOE isoform genotype in African-American (AA) and European-American (EA) populations. The e2 allele is associated with strong correlation between ln TG and TC while the e4 allele leads to little or no correlation. This led to a priori hypotheses that APOE genotypes affect the relationship of TC and/or ln TG with incident CHD. We found that APOE*TC was significant (P = 0.016) for AA but not EA while APOE*ln TG was significant for EA (P = 0.027) but not AA. In both cases, e2e2 and e2e3 had strong relationships between TC and ln TG with CHD while e2e4 and e4e4 results in little or no relationship between TC and ln TG with CHD. Using ARIC GWAS data, scans for loci that significantly interact with APOE produced four loci for African Americans (one CHD, one TC, and two HDL). These interactions contribute to the rQTL pattern. rQTL are a powerful tool to identify loci that modify the relationship between risk factors and disease and substantially increase statistical power for detecting G × G.

Ontogenetic and static allometries describe how a character changes in size when the size of the organism changes during ontogeny and among individuals measured at the same developmental stage, respectively. Understanding the relationship between these two types of allometry is crucial to understanding the evolution of allometry and, more generally, the evolution of shape. However, the effects of ontogenetic allometry on static allometry remain largely unexplored. Here, we first show analytically how individual variation in ontogenetic allometry and body size affect static allometry. Using two longitudinal data sets on ontogenetic and static allometry, we then estimate variances and covariances for the different parameters of the ontogenetic allometry defined in our model and assess their relative contribution to the static allometric slope. The mean ontogenetic allometry is the main parameter that determines the static allometric slope, while the covariance between the ontogenetic allometric slope and body size generates most of the discrepancies between ontogenetic and static allometry. These results suggest that the apparent evolutionary stasis of the static allometric slope is not generated by internal (developmental) constraints but more likely results from external constraints imposed by selection. © 2013 by The University of Chicago. All rights reserved.

Nonalcoholic fatty liver disease, a condition in which excess fat accumulates in the liver, is strongly associated with the metabolic syndrome, including obesity and other related conditions. This disease has the potential to progress from steatosis to steatohepatitis, fibrosis, and cirrhosis. The recent increase in the prevalence of the metabolic syndrome is largely driven by changes in diet and activity levels. Individual variation in the response to this obesogenic environment, however, is attributable in part to genetic variation between individuals, but very few mammalian genetic loci have been identified with effects on fat accumulation in the liver. To study the genetic basis for variation in liver fat content in response to dietary fat, liver fat proportion was determined using quantitative magnetic resonance imaging in 478 mice from 16 LG/J X SM/J recombinant inbred strains fed either a high-fat (42% kcal from fat) or low-fat (15% kcal from fat) diet. An analysis of variance confirmed that there is a genetic basis for variation in liver fat content within the population with significant effects of sex and diet. Three quantitative trail loci that contribute to liver fat content also were mapped.

We present full-genome genotype imputations for 100 classical laboratory mouse strains, using a novel method. Using genotypes at 549,683 SNP loci obtained with the Mouse Diversity Array, we partitioned the genome of 100 mouse strains into 40,647 intervals that exhibit no evidence of historical recombination. For each of these intervals we inferred a local phylogenetic tree. We combined these data with 12 million loci with sequence variations recently discovered by whole-genome sequencing in a common subset of 12 classical laboratory strains. For each phylogenetic tree we identified strains sharing a leaf node with one or more of the sequenced strains. We then imputed high- and medium-confidence genotypes for each of 88 nonsequenced genomes. Among inbred strains, we imputed 92% of SNPs genome-wide, with 71% in high-confidence regions. Our method produced 977 million new genotypes with an estimated per-SNP error rate of 0.083% in high-confidence regions and 0.37% genome-wide. Our analysis identified which of the 88 nonsequenced strains would be the most informative for improving full-genome imputation, as well as which additional strain sequences will reveal more new genetic variants. Imputed sequences and quality scores can be downloaded and visualized online.

Inbred mouse strains MRL and LG share the ability to fully heal ear hole punches with the full range of appropriate tissues without scarring. They also share a common ancestry, MRL being formed from a multi-strain cross with two final backcrosses to LG before being inbred by brother-sister mating. Many gene-mapping studies for healing ability have been performed using these two strains, resulting in the location of about 20 quantitative trait loci (QTLs). Here, we combine two of these crosses (N = 638), MRL/lpr × C57BL/6NTac and LG/J × SM/J, in a single combined cross analysis to increase the mapping power, decrease QTL support intervals, separate multiple QTLs and establish allelic states at individual QTL. The combined cross analysis located 11 QTLs, 6 affecting only one cross (5 LG × SM and 1 MRL × B6) and 5 affecting both crosses, approximately the number of common QTLs expected given strain SNP similarity. Amongst the five QTLs mapped in both crosses, three had significantly different genetic effects, additive in one cross and over or underdominant in the other. It is possible that allelic states at these three loci are different in SM and B6 because they lead to differences in dominance interactions with the LG and MRL alleles. QTL support intervals are 40% smaller in the combined cross analysis than in either of the single crosses. Combined cross analysis was successful in enhancing the interpretation of earlier QTL results for these strains.

It is well known that standard population genetic theory predicts decreased additive genetic variance (V(a) ) following a population bottleneck and that theoretical models including interallelic and intergenic interactions indicate such loss may be avoided. However, few empirical data from multicellular model systems are available, especially regarding variance/covariance (V/CV) relationships. Here, we compare the V/CV structure of seventeen traits related to body size and composition between control (60 mating pairs/generation) and bottlenecked (2 mating pairs/generation; average F = 0.39) strains of mice. Although results for individual traits vary considerably, multivariate analysis indicates that V(a) in the bottlenecked populations is greater than expected. Traits with patterns and amounts of epistasis predictive of enhanced V(a) also show the largest deviations from additive expectations. Finally, the correlation structure of weekly weights is not significantly different between control and experimental lines but correlations between necropsy traits do differ, especially those involving the heart, kidney and tail length.

Reports from studies of twins, disease aggregation in families, animal models for periodontal disease, and various genetic-analysis studies have determined that genetics plays a role in the susceptibility to periodontal disease. The purpose of this pilot study is to evaluate the effect of genetics on periodontal disease by evaluating the heritability of alveolar bone loss in a captive baboon population.A collection of baboon skulls from a pedigreed colony (for which scientists and veterinarians maintain complete genealogic and veterinary records) was obtained from the Southwest National Primate Research Center, San Antonio, Texas and used in this pilot study. Measurements of alveolar bone loss were performed on 390 dry baboon skulls. A periodontal probe was used to measure alveolar bone loss. Maximum likelihood methods (designed to handle complex genealogies) were used to determine the heritability of alveolar bone loss. This software used known pedigrees in the captive baboon sample and tested the relationship between pairwise kinship and alveolar bone loss data to determine the heritability of alveolar bone loss from periodontal disease.Genetic data were available for 347 of the 390 specimens. Using age and sex as covariates, genetic analysis indicated a heritability of 35% (standard error = 20%; P = 0.01). Although gender was not a significant factor in periodontal disease (P = 0.96), age was highly significantly associated with periodontal disease (P0.0001).In this pilot study, analysis of alveolar bone loss measurements from captive baboons indicates that bone loss increases with age and that a portion of periodontal disease risk may be caused by genetic variance. These findings provide evidence that periodontal disease is heritable in captive baboons and indicate that a larger, more-detailed study is warranted.

Genome-wide mapping analyses are now commonplace in many species and several networks of interacting loci have been reported. However, relatively few details regarding epistatic interactions and their contribution to complex trait variation in multicellular organisms are available and the identification of positional candidate loci for epistatic QTL (epiQTL) is hampered, especially in mammals, by the limited genetic resolution inherent in most study designs. Here we further investigate the genetic architecture of reproductive fatpad weight in mice using the F10 generation of the LG,SM advanced intercross (AI) line. We apply multiple mapping techniques including a single-locus model, locus-specific composite interval mapping (CIM), and tests for multiple QTL per chromosome to the 12 chromosomes known to harbor single-locus QTL (slQTL) affecting obesity in this cross. We also perform a genome-wide scan for pairwise epistasis. Using this combination of approaches we detect 199 peaks spread over all 19 autosomes, which potentially contribute to trait variation including all eight original F2 loci (Adip1-8), novel slQTL peaks on chromosomes 7 and 9, and several novel epistatic loci. Extensive epistasis is confirmed involving both slQTL confidence intervals (C.I.) as well as regions that show no significant additive or dominance effects. These results provide important new insights into mapping complex genetic architectures and the role of epistasis in complex trait variation.

Variations in diabetic phenotypes are caused by complex interactions of genetic effects, environmental factors, and the interplay between the two. We tease apart these complex interactions by examining genome-wide genetic and epigenetic effects on diabetes-related traits among different sex, diet, and sex-by-diet cohorts in a Mus musculus model. We conducted a genome-wide scan for quantitative trait loci that affect serum glucose and insulin levels and response to glucose stress in an F(16) Advanced Intercross Line of the LG/J and SM/J intercross (Wustl:LG,SM-G16). Half of each sibship was fed a high-fat diet and half was fed a relatively low-fat diet. Context-dependent genetic (additive and dominance) and epigenetic (parent-of-origin imprinting) effects were characterized by partitioning animals into sex, diet, and sex-by-diet cohorts. We found that different cohorts often have unique genetic effects at the same loci, and that genetic signals can be masked or erroneously assigned to specific cohorts if they are not considered individually. Our data demonstrate that the effects of genes on complex trait variation are highly context-dependent and that the same genomic sequence can affect traits differently depending on an individual's sex and/or dietary environment. Our results have important implications for studies of complex traits in humans.

Previous studies on the LG,SM advanced intercross line have identified approximately 40 quantitative trait loci (QTL) for long -bone (humerus, ulna, femur, and tibia) lengths. In this study, long-bone-length QTL were fine-mapped in the F(34) generation (n = 1424) of the LG,SM advanced intercross. Environmental effects were assessed by dividing the population by sex between high-fat and low-fat diets, producing eight sex/diet cohorts. We identified 145 individual bone-length QTL comprising 45 pleiotropic QTL; 69 replicated QTL from previous studies, 35 were new traits significant at previously identified loci, and 41 were novel QTL. Many QTL affected only a subset of the population based on sex and/or diet. Eight of ten known skeletal growth genes were upregulated in 3-week-old LG/J male proximal tibial growth plates relative to SM/J. The sequences of parental strains LG/J and SM/J indicated the presence of over half a million polymorphisms in the confidence intervals of these 45 QTL. We examined 526 polymorphisms and found that 97 represented radical changes to amino acid composition while 40 were predicted to be deleterious to protein function. Additional experimentation is required to understand how changes in gene regulation or protein function can alter the genetic architecture and interact with the environment to produce phenotypic variation.

Directional epistasis describes a situation in which epistasis consistently increases or decreases the effect of allele substitutions, thereby affecting the amount of additive genetic variance available for selection in a given direction. This study applies a recent parameterization of directionality of epistasis to empirical data. Data stems from a QTL mapping study on an intercross between inbred mouse (Mus musculus) strains LG/J and SM/J, originally selected for large and small body mass, respectively. Results show a negative average directionality of epistasis for body-composition traits, predicting a reduction in additive allelic effects and in the response to selection for increased size. Focusing on average modification of additive effect of single loci, we find a more complex picture, whereby the effects of some loci are enhanced consistently across backgrounds, while effects of other loci are decreased, potentially contributing to either enhancement or reduction of allelic effects when selection acts at single loci. We demonstrate and discuss how the interpretation of the overall measurement of directionality depends on the complexity of the genotype–phenotype map. The measure of directionality changes with the power of scale in a predictable way; however, its expected effect with respect to the modification of additive genetic effects remains constant.

Animals capable of regenerating multiple tissue types, organs, and appendages after injury are common yet sporadic and include some sponge, hydra, planarian, and salamander (i.e., newt and axolotl) species, but notably such regenerative capacity is rare in mammals. The adult MRL mouse strain is a rare exception to the rule that mammals do not regenerate appendage tissue. Certain commonalities, such as blastema formation and basement membrane breakdown at the wound site, suggest that MRL mice may share other features with classical regenerators. As reported here, MRL fibroblast-like cells have a distinct cell-cycle (G2/M accumulation) phenotype and a heightened basal and wound site DNA damage/repair response that is also common to classical regenerators and mammalian embryonic stem cells. Additionally, a neutral and alkaline comet assay displayed a persistent level of intrinsic DNA damage in cells derived from the MRL mouse. Similar to mouse ES cells, the p53-target p21 was not expressed in MRL ear fibroblasts. Because the p53/p21 axis plays a central role in the DNA damage response and cell cycle control, we directly tested the hypothesis that p21 down-regulation could functionally induce a regenerative response in an appendage of an otherwise nonregenerating mouse strain. Using the ear hole closure phenotype, a genetically mapped and reliable quantitative indicator of regeneration in the MRL mouse, we show that the unrelated Cdkn1a tmi/Tyj /J p21 −/− mouse (unlike the B6129SF2/J WT control) closes ear holes similar to MRL mice, providing a firm link between cell cycle checkpoint control and tissue regeneration.

Animal models have enriched understanding of the physiological basis of metabolic disorders and advanced identification of genetic risk factors underlying the metabolic syndrome (MetS). Murine models are especially appropriate for this type of research, and are an excellent resource not only for identifying candidate genomic regions, but also for illuminating the possible molecular mechanisms or pathways affected in individual components of MetS. In this review, we briefly discuss findings from mouse models of metabolic disorders, particularly in light of issues raised by the recent flood of human genome-wide association studies (GWAS) results. We describe how mouse models are revealing that genotype interacts with environment in important ways, indicating that the underlying genetics of MetS is highly context dependant. Further we show that epistasis, imprinting and maternal effects each contribute to the genetic architecture underlying variation in metabolic traits, and mouse models provide an opportunity to dissect these aspects of the genetic architecture that are difficult if not impossible to ascertain in humans. Finally we discuss how knowledge gained from mouse models can be used in conjunction with comparative genomic methods and bioinformatic resources to inform human MetS research.

We investigated maternal obesity in inbred SM/J mice by assigning females to a high-fat diet or a low-fat diet at weaning, mating them to low-fat-fed males, cross-fostering the offspring to low-fat-fed SM/J nurses at birth, and weaning the offspring onto a high-fat or low-fat diet. A maternal high-fat diet exacerbated obesity in the high-fat-fed daughters, causing them to weigh more, have more fat, and have higher serum levels of leptin as adults, accompanied by dozens of gene expression changes and thousands of DNA methylation changes in their livers and hearts. Maternal diet particularly affected genes involved in RNA processing, immune response, and mitochondria. Between one-quarter and one-third of differentially expressed genes contained a differentially methylated region associated with maternal diet. An offspring high-fat diet reduced overall variation in DNA methylation, increased body weight and organ weights, increased long bone lengths and weights, decreased insulin sensitivity, and changed the expression of 3,908 genes in the liver. Although the offspring were more affected by their own diet, their maternal diet had epigenetic effects lasting through adulthood, and in the daughters these effects were accompanied by phenotypic changes relevant to obesity and diabetes.

MRL mice display unusual healing properties. When MRL ear pinnae are hole punched, the holes close completely without scarring, with regrowth of cartilage and reappearance of both hair follicles and sebaceous glands. Studies using (MRL/lpr x C57BL/6)F(2) and backcross mice first showed that this phenomenon was genetically determined and that multiple loci contributed to this quantitative trait. The lpr mutation itself, however, was not one of them. In the present study we examined the genetic basis of healing in the Large (LG/J) mouse strain, a parent of the MRL mouse and a strain that shows the same healing phenotype. LG/J mice were crossed with Small (SM/J) mice and the F(2) population was scored for healing and their genotypes determined at more than 200 polymorphic markers. As we previously observed for MRL and (MRL x B6)F(2) mice, the wound-healing phenotype was sexually dimorphic, with female mice healing more quickly and more completely than male mice. We found quantitative trait loci (QTLs) on chromosomes (Chrs) 9, 10, 11, and 15. The heal QTLs on Chrs 11 and 15 were linked to differential healing primarily in male animals, whereas QTLs on Chrs 9 and 10 were not sexually dimorphic. A comparison of loci identified in previous crosses with those in the present report using LG/J x SM/J showed that loci on Chrs 9, 11, and 15 colocalized with those seen in previous MRL crosses, whereas the locus on Chr 10 was not seen before and is contributed by SM/J.

Genomic imprinting, where the effects of alleles depend on their parent-of-origin, can be an important component of the genetic architecture of complex traits. Although there has been a rapidly increasing number of studies of genetic architecture that have examined imprinting effects, none have examined whether imprinting effects depend on genetic background. Such effects are critical for the evolution of genomic imprinting because they allow the imprinting state of a locus to evolve as a function of genetic background. Here we develop a two-locus model of epistasis that includes epistatic interactions involving imprinting effects and apply this model to scan the mouse genome for loci that modulate the imprinting effects of quantitative trait loci (QTL). The inclusion of imprinting leads to nine orthogonal forms of epistasis, five of which do not appear in the usual two-locus decomposition of epistasis. Each form represents a change in the imprinting status of one locus across different classes of genotypes at the other locus. Our genome scan identified two different locus pairs that show complex patterns of epistasis, where the imprinting effect at one locus changes across genetic backgrounds at the other locus. Thus, our model provides a framework for the detection of genetic background-dependent imprinting effects that should provide insights into the background dependence and evolution of genomic imprinting. Our application of the model to a genome scan supports this assertion by identifying pairs of loci that show reciprocal changes in their imprinting status as the background provided by the other locus changes.

Autoresuscitation (AR) is a highly conserved response among mammals, which allows survival from transient extreme hypoxia. During hypoxia, bradycardia, and hypoxic gasping develop after a brief period of hyperactivity. Normally, AR occurs if oxygen is restored during the gasping period where an initial heart rate increase is rapidly followed resumption or eupneic breathing. Humans and other mammals can survive multiple immediately repeated AR. A defective AR capacity has been implicated in Sudden Infant Death Syndrome. We had reported earlier that inbred strains of mice such as BALB/cJ could survive a characteristic number of immediately repeated AR trials, but that SWR/J mice failed to AR from a single hypoxic episode. We now report that strains closely related to SWR/J, FVB/N and SJL/J exhibit partial resuscitation defects relative to BALB/cJ or other mouse strains, establishing a genetic basis for variation in AR failure. The AR trial phenotype of BALB/cJ x SWR/J intercross F(1) and F(2) mice was consistent with BALB/cJ dominance and a discrete number of loci. Genome-wide mapping conducted with 60 intercross F(2) animals linked two loci to the number of AR trials survived, including one sex-specific locus with male expression, consistent with the observed 50% male bias for Sudden Infant Death Syndrome in humans. A locus carried on SWR/J chromosome 10 seems to be particularly important in AR failure and was confirmed in a partial consomic line. These results establish a genetic basis for AR failure phenotype in mice, with relevance to Sudden Infant Death Syndrome.

Theories of phenotypic integration have relied heavily on the concept of modularity in order to model the ways in which traits in an organism correlate and covary. Recent investigations suggest that, while some functional and developmental processes may be morphologically and ontogenetically localized, and thus modular in a developmental sense, there is a great deal of overlap among these influences on patterns of integration in the adult form. This can result in blurry boundaries between hypothesized modules constructed to test hypotheses about phenotypic integration. This investigation tests hypotheses about the contribution of pleiotropic quantitative trait loci (QTL) to phenotypic integration in the mouse mandible without using a priori categorical hypotheses about which traits constitute a module. We ask two main questions: (1) Are the effects of pleiotropic QTL localized to highly correlated traits or more spread out among traits than one might expect by chance? (2) Does the pattern of trait influence when all pleiotropic QTL are considered together deviate from what we might expect if QTL affect traits without regard for the correlations among traits? We find that a large proportion of pleiotropic QTL affect traits that are more highly correlated than we expect by chance with the remainder having effects that are distributed as if by chance. Furthermore, the overall distribution of the effects of pleiotropic QTL differs significantly from the null distribution of no association between pleiotropic effects on traits and correlations among traits. The main modular hypothesis used by earlier studies often does not predict the distribution of sets of traits sharing a common QTL. These results suggest that there is a clear tendency for pleiotropic effects of QTL to be localized but that the localization may be best thought of as occurring in a continuous space rather being clustered in discrete modules.