microRNA-146a-5p association with the cardiometabolic disease risk factor TMAO.

Trimethylamine-N-oxide (TMAO), a microbial choline metabolism byproduct that is processed in the liver and excreted into circulation, is associated with increased atherosclerotic lesion formation and cardiovascular disease risk. Genetic regulators of TMAO levels are largely unknown. In the present study, we used 288 mice from a genetically heterogeneous mouse population [Diversity Outbred (DO)] to determine hepatic microRNA associations with TMAO in the context of an atherogenic diet. We also validated findings in two additional animal models of atherosclerosis: liver-specific insulin receptor knockout mice fed a chow diet (LIRKO) and African green monkeys fed high-fat/high-cholesterol diet. Small RNA-sequencing analysis in DO mice, LIRKO mice, and African green monkeys identified only one hepatic microRNA (miR-146a-5p) that is aberrantly expressed across all three models. Moreover, miR-146a-5p levels are associated with circulating TMAO after atherogenic diet in each of these models. We also performed high-resolution genetic mapping and identified a novel quantitative trait locus on Chromosome 12 for TMAO levels. This interval includes two genes, Numb and Dlst, which are inversely correlated with both miR-146a and TMAO and are predicted targets of miR-146a. Both of these genes have been validated as direct targets of miR-146a, though in other cellular contexts. This is the first report to our knowledge of a link between miR-146 and TMAO. Our findings suggest that miR-146-5p, as well as one or more genes at the Chromosome 12 QTL (possibly Numb or Dlst), is strongly linked to TMAO levels and likely involved in the control of atherosclerosis.

Correction to: Quantitative trait mapping in Diversity Outbred mice identifies two genomic regions associated with heart size.

The original article has been published with an incorrect text in Materials and Methods section. The corrected text should read as.

Quantitative trait mapping in Diversity Outbred mice identifies two genomic regions associated with heart size.

Heart size is an important factor in cardiac health and disease. In particular, increased heart weight is predictive of adverse cardiovascular outcomes in multiple large community-based studies. We use two cohorts of Diversity Outbred (DO) mice to investigate the role of genetics, sex, age, and diet on heart size. DO mice (n = 289) of both sexes from generation 10 were fed a standard chow diet, and analyzed at 12-15 weeks of age. Another cohort of female DO mice (n = 258) from generation 11 were fed either a high-fat, cholesterol-containing (HFC) diet or a low-fat, high-protein diet, and analyzed at 24-25 weeks. We did not observe an effect of diet on body or heart weight in generation 11 mice, although we previously reported an effect on other cardiovascular risk factors, including cholesterol, triglycerides, and insulin. We do observe a significant genetic effect on heart weight in this population. We identified two quantitative trait loci for heart weight, one (Hwtf1) at a genome-wide significance level of p ≤ 0.05 on MMU15 and one (Hwtf2) at a genome-wide suggestive level of p ≤ 0.1 on MMU10, that together explain 13.3% of the phenotypic variance. Hwtf1 contained collagen type XXII alpha 1 chain (Col22a1), and the NZO/HlLtJ and WSB/EiJ haplotypes were associated with larger hearts. This is consistent with heart tissue Col22a1 expression in DO founders and SNP patterns within Hwtf1 for Col22a1. Col22a1 has been previously associated with cardiac fibrosis in mice, suggesting that Col22a1 may be involved in pathological cardiac hypertrophy.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

Systems genetics identifies a co-regulated module of liver microRNAs associated with plasma LDL cholesterol in murine diet-induced dyslipidemia.

Chronically altered levels of circulating lipids, termed dyslipidemia, is a significant risk factor for a number of metabolic and cardiovascular morbidities. MicroRNAs (miRNAs) have emerged as important regulators of lipid balance, have been implicated in dyslipidemia, and have been proposed as candidate therapeutic targets in lipid-related disorders including atherosclerosis. A major limitation of most murine studies of miRNAs in lipid metabolic disorders is that they have been performed in just one (or very few) inbred strains, such as C57BL/6. Moreover, although individual miRNAs have been associated with lipid phenotypes, it is well understood that miRNAs likely work together in functional modules. To address these limitations, we implemented a systems genetics strategy using the Diversity Outbred (DO) mouse population. Specifically, we performed gene and miRNA expression profiling in the livers from ~300 genetically distinct DO mice after 18 wk on either a high-fat/high-cholesterol diet or a high-protein diet. Large-scale correlative analysis of these data with a wide range of cardio-metabolic end points revealed a co-regulated module of miRNAs significantly associated with circulating low-density lipoprotein cholesterol (LDL-C) levels. The hubs of this module were identified as miR-199a, miR-181b, miR-27a, miR-21_-_1, and miR-24. In sum, we demonstrate that a high-fat/high-cholesterol diet robustly rewires the miRNA regulatory network, and we identify a small group of co-regulated miRNAs that may exert coordinated effects to control circulating LDL-C.

R2d2 Drives Selfish Sweeps in the House Mouse.

A selective sweep is the result of strong positive selection driving newly occurring or standing genetic variants to fixation, and can dramatically alter the pattern and distribution of allelic diversity in a population. Population-level sequencing data have enabled discoveries of selective sweeps associated with genes involved in recent adaptations in many species. In contrast, much debate but little evidence addresses whether "selfish" genes are capable of fixation-thereby leaving signatures identical to classical selective sweeps-despite being neutral or deleterious to organismal fitness. We previously described R2d2, a large copy-number variant that causes nonrandom segregation of mouse Chromosome 2 in females due to meiotic drive. Here we show population-genetic data consistent with a selfish sweep driven by alleles of R2d2 with high copy number (R2d2(HC)) in natural populations. We replicate this finding in multiple closed breeding populations from six outbred backgrounds segregating for R2d2 alleles. We find that R2d2(HC) rapidly increases in frequency, and in most cases becomes fixed in significantly fewer generations than can be explained by genetic drift. R2d2(HC) is also associated with significantly reduced litter sizes in heterozygous mothers, making it a true selfish allele. Our data provide direct evidence of populations actively undergoing selfish sweeps, and demonstrate that meiotic drive can rapidly alter the genomic landscape in favor of mutations with neutral or even negative effects on overall Darwinian fitness. Further study will reveal the incidence of selfish sweeps, and will elucidate the relative contributions of selfish genes, adaptation and genetic drift to evolution.

Maternal exercise before and during pregnancy does not impact offspring exercise or body composition in mice.

BACKGROUND: The genome, the environment, and their interactions simultaneously regulate complex traits such as body composition and voluntary exercise levels. One such environmental influence is the maternal milieu (i.e., in utero environment or maternal care). Variability in the maternal environment may directly impact the mother, and simultaneously has the potential to influence the physiology and/or behavior of offspring in utero, post birth, and into adulthood. Here, we utilized a murine model to examine the effects of the maternal environment in regard to voluntary exercise (absence of wheel running, wheel running prior to gestation, and wheel running prior to and throughout gestation) on offspring weight and body composition (% fat tissue and % lean tissue) throughout development (~3 to ~9 weeks of age). Additionally, we examined the effects of ~6 weeks of maternal exercise (prior to and during gestation) on offspring exercise levels at ~9 weeks of age. RESULTS: We observed no substantial effects of maternal exercise on subsequent male or female offspring body composition throughout development, or on the propensity of offspring to engage in voluntary wheel running. At the level of the individual, correlational analyses revealed some statistically significant relationships between maternal and offspring exercise levels, likely reflecting previously known heritability estimates for such traits. CONCLUSIONS: The current results conflict with previous findings in human and mouse models demonstrating that maternal exercise has the potential to alter offspring phenotypes. We discuss our negative findings in the context of the timing of the maternal exercise and the level of biological organization of the examined phenotypes within the offspring.

Functional characterization of cytochrome P450-derived epoxyeicosatrienoic acids in adipogenesis and obesity.

Adipogenesis plays a critical role in the initiation and progression of obesity. Although cytochrome P450 (CYP)-derived epoxyeicosatrienoic acids (EETs) have emerged as a potential therapeutic target for cardiometabolic disease, the functional contribution of EETs to adipogenesis and the pathogenesis of obesity remain poorly understood. Our studies demonstrated that induction of adipogenesis in differentiated 3T3-L1 cells (in vitro) and obesity-associated adipose expansion in high-fat diet (HFD)-fed mice (in vivo) significantly dysregulate the CYP epoxygenase pathway and evoke a marked suppression of adipose-derived EET levels. Subsequent in vitro experiments demonstrated that exogenous EET analog administration elicits potent anti-adipogenic effects via inhibition of the early phase of adipogenesis. Furthermore, EET analog administration to mice significantly mitigated HFD-induced weight gain, adipose tissue expansion, pro-adipogenic gene expression, and glucose intolerance. Collectively, these findings suggest that suppression of EET bioavailability in adipose tissue is a key pathological consequence of obesity, and strategies that promote the protective effects of EETs in adipose tissue offer enormous therapeutic potential for obesity and its downstream pathological consequences.

Quantitative genomics of voluntary exercise in mice: transcriptional analysis and mapping of expression QTL in muscle.

Motivation and ability both underlie voluntary exercise, each with a potentially unique genetic architecture. Muscle structure and function are one of many morphological and physiological systems acting to simultaneously determine exercise ability. We generated a large (n = 815) advanced intercross line of mice (G4) derived from a line selectively bred for increased wheel running (high runner) and the C57BL/6J inbred strain. We previously mapped quantitative trait loci (QTL) contributing to voluntary exercise, body composition, and changes in body composition as a result of exercise. Using brain tissue in a subset of the G4 (n = 244), we have also previously reported expression QTL (eQTL) colocalizing with the QTL for the higher-level phenotypes. Here, we examined the transcriptional landscape of hind limb muscle tissue via global mRNA expression profiles. Correlations revealed an ∼1,168% increase in significant relationships between muscle transcript expression levels and the same exercise and body composition phenotypes examined previously in the brain. The exercise trait most often significantly correlated with gene expression in the brain was running duration while in the muscle it was maximum running speed. This difference may indicate that time spent engaging in exercise behavior may be more influenced by central (neurobiological) mechanisms, while intensity of exercise may be largely controlled by peripheral mechanisms. Additionally, we used subsets of cis-acting eQTL, colocalizing with QTL, to identify candidate genes based on both positional and functional evidence. We discuss three plausible candidate genes (Insig2, Prcp, Sparc) and their potential regulatory role.

Genetic analysis of complex traits in the emerging Collaborative Cross.

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors.

In vertebrates, including humans, individuals harbor gut microbial communities whose species composition and relative proportions of dominant microbial groups are tremendously varied. Although external and stochastic factors clearly contribute to the individuality of the microbiota, the fundamental principles dictating how environmental factors and host genetic factors combine to shape this complex ecosystem are largely unknown and require systematic study. Here we examined factors that affect microbiota composition in a large (n = 645) mouse advanced intercross line originating from a cross between C57BL/6J and an ICR-derived outbred line (HR). Quantitative pyrosequencing of the microbiota defined a core measurable microbiota (CMM) of 64 conserved taxonomic groups that varied quantitatively across most animals in the population. Although some of this variation can be explained by litter and cohort effects, individual host genotype had a measurable contribution. Testing of the CMM abundances for cosegregation with 530 fully informative SNP markers identified 18 host quantitative trait loci (QTL) that show significant or suggestive genome-wide linkage with relative abundances of specific microbial taxa. These QTL affect microbiota composition in three ways; some loci control individual microbial species, some control groups of related taxa, and some have putative pleiotropic effects on groups of distantly related organisms. These data provide clear evidence for the importance of host genetic control in shaping individual microbiome diversity in mammals, a key step toward understanding the factors that govern the assemblages of gut microbiota associated with complex diseases.

Lack of both bradykinin B1 and B2 receptors enhances nephropathy, neuropathy, and bone mineral loss in Akita diabetic mice.

An insertion polymorphism of the angiotensin-I converting enzyme gene (ACE) is common in humans and the higher expressing allele is associated with an increased risk of diabetic complications. The ACE polymorphism does not significantly affect blood pressure or angiotensin II levels, suggesting that the kallikrein-kinin system partly mediates the effects of the polymorphism. We have therefore explored the influence of lack of both bradykinin receptors (B1R and B2R) on diabetic nephropathy, neuropathy, and osteopathy in male mice heterozygous for the Akita diabetogenic mutation in the insulin 2 gene (Ins2). We find that all of the detrimental phenotypes observed in Akita diabetes are enhanced by lack of both B1R and B2R, including urinary albumin excretion, glomerulosclerosis, glomerular basement membrane thickening, mitochondrial DNA deletions, reduction of nerve conduction velocities and of heat sensation, and bone mineral loss. Absence of the bradykinin receptors also enhances the diabetes-associated increases in plasma thiobarbituric acid-reactive substances, mitochondrial DNA deletions, and renal expression of fibrogenic genes, including transforming growth factor beta1, connective tissue growth factor, and endothelin-1. Thus, lack of B1R and B2R exacerbates diabetic complications. The enhanced renal injury in diabetic mice caused by lack of B1R and B2R may be mediated by a combination of increases in oxidative stress, mitochondrial DNA damage and over expression of fibrogenic genes.

A suite of engineered mice for interrogating psychedelic drug actions.

Psychedelic drugs like lysergic acid diethylamide (LSD) and psilocybin have emerged as potentially transformative therapeutics for many neuropsychiatric diseases, including depression, anxiety, post-traumatic stress disorder, migraine, and cluster headaches. LSD and psilocybin exert their psychedelic effects via activation of the 5-hydroxytryptamine 2A receptor (HTR2A). Here we provide a suite of engineered mice useful for clarifying the role of HTR2A and HTR2A-expressing neurons in psychedelic drug actions. We first generated Htr2a-EGFP-CT-IRES-CreERT2 mice (CT:C-terminus) to independently identify both HTR2A-EGFP-CT receptors and HTR2A-containing cells thereby providing a detailed anatomical map of HTR2A and identifying cell types that express HTR2A. We also generated a humanized Htr2a mouse line and an additional constitutive Htr2A-Cre mouse line. Psychedelics induced a variety of known behavioral changes in our mice validating their utility for behavioral studies. Finally, electrophysiology studies revealed that extracellular 5-HT elicited a HTR2A-mediated robust increase in firing of genetically-identified pyramidal neurons--consistent with a plasma membrane localization and mode of action. These mouse lines represent invaluable tools for elucidating the molecular, cellular, pharmacological, physiological, behavioral, and other actions of psychedelic drugs in vivo.

Increased energy expenditure, dietary fat wasting, and resistance to diet-induced obesity in mice lacking renin.

An overactive renin-angiotensin system is associated with obesity and the metabolic syndrome. However, the mechanisms behind it are unclear. Cleaving angiotensinogen to angiotensin I by renin is a rate-limiting step of angiotensin II production, but renin is suggested to have angiotensin-independent effects. We generated mice lacking renin (Ren1c) using embryonic stem cells from C57BL/6 mice, a strain prone to diet-induced obesity. Ren1c(-/-) mice are lean, insulin sensitive, and resistant to diet-induced obesity without changes in food intake and physical activity. The lean phenotype is likely due to a higher metabolic rate and gastrointestinal loss of dietary fat. Most of the metabolic changes in Ren1c(-/-) mice were reversed by angiotensin II administration. These results support a role for angiotensin II in the pathogenesis of diet-induced obesity and insulin resistance.

Exercise and diet affect quantitative trait loci for body weight and composition traits in an advanced intercross population of mice.

Driven by the recent obesity epidemic, interest in understanding the complex genetic and environmental basis of body weight and composition is great. We investigated this by searching for quantitative trait loci (QTLs) affecting a number of weight and adiposity traits in a G(10) advanced intercross population produced from crosses of mice in inbred strain C57BL/6J with those in a strain selected for high voluntary wheel running. The mice in this population were fed either a high-fat or a control diet throughout the study and also measured for four exercise traits prior to death, allowing us to test for pre- and postexercise QTLs as well as QTL-by-diet and QTL-by-exercise interactions. Our genome scan uncovered a number of QTLs, of which 40% replicated QTLs previously found for similar traits in an earlier (G(4)) generation. For those replicated QTLs, the confidence intervals were reduced from an average of 19 Mb in the G(4) to 8 Mb in the G(10). Four QTLs on chromosomes 3, 8, 13, and 18 were especially prominent in affecting the percentage of fat in the mice. About of all QTLs showed interactions with diet, exercise, or both, their genotypic effects on the traits showing a variety of patterns depending on the diet or level of exercise. It was concluded that the indirect effects of these QTLs provide an underlying genetic basis for the considerable variability in weight or fat loss typically found among individuals on the same diet and/or exercise regimen.

Exercise, weight loss, and changes in body composition in mice: phenotypic relationships and genetic architecture.

The regulation of body weight and composition is complex, simultaneously affected by genetic architecture, the environment, and their interactions. We sought to analyze the complex phenotypic relationships between voluntary exercise, food consumption, and changes in body weight and composition and simultaneously localize quantitative trait loci (QTL) controlling these traits. A large (n = 815) murine advanced intercross line (G(4)) was created from a reciprocal cross between a high-running line and the inbred strain C57BL/6J. Body weight and composition (% fat, % lean) were measured at 4, 6, and 8 wk of age. After measurements at 8 wk of age, mice were given access to running wheels, during which food consumption was quantified and after which body weight and composition were assessed to evaluate exercise-induced changes. Phenotypic correlations indicated that the relationship between exercise and overall change in weight and adiposity depended on body composition before the initiation of exercise. Interval mapping revealed QTL for body weight, % fat, and % lean at 4, 6, and 8 wk of age. Furthermore, QTL were observed for food consumption and changes in weight, % fat, and % lean in response to short-term exercise. Here we provide some clarity for the relationship between weight loss, reduction in adiposity, food consumption, and exercise. Simultaneously, we reinforce the genetic basis for body weight and composition with some independent loci controlling growth at different ages. Finally, we present unique QTL providing insight regarding variation in weight loss and reduction in adiposity in response to exercise.

Parent-of-origin effects on voluntary exercise levels and body composition in mice.

Despite the health-related benefits of exercise, many people do not engage in enough activity to realize the rewards, and little is known regarding the genetic or environmental components that account for this individual variation. We created and phenotyped a large G(4) advanced intercross line originating from reciprocal crosses between mice with genetic propensity for increased voluntary exercise (HR line) and the inbred strain C57BL/6J. G(4) females (compared to males) ran significantly more when provided access to a running wheel and were smaller with a greater percentage of body fat pre- and postwheel access. Change in body composition resulting from a 6-day exposure to wheels varied between the sexes with females generally regulating energy balance more precisely in the presence of exercise. We observed parent-of-origin effects on most voluntary wheel running and body composition traits, which accounted for 3-13% of the total phenotypic variance pooled across sexes. G(4) individuals descended from progenitor (F(0)) crosses of HRfemale symbol and C57BL/6Jmale symbol ran greater distances, spent more time running, ran at higher maximum speeds/day, and had lower percent body fat and higher percent lean mass than mice descended from reciprocal progenitor crosses (C57BL/6Jfemale symbol x HRmale symbol). For some traits, significant interactions between parent of origin and sex were observed. We discuss these results in the context of sex dependent activity and weight loss patterns, the contribution of parent-of-origin effects to predisposition for voluntary exercise, and the genetic (i.e., X-linked or mtDNA variations), epigenetic (i.e., genomic imprinting), and environmental (i.e., in utero environment or maternal care) phenomena potentially modulating these effects.

Genetic architecture of voluntary exercise in an advanced intercross line of mice.

Exercise is essential for health, yet the amount, duration, and intensity that individuals engage in are strikingly variable, even under prescription. Our focus was to identify the locations and effects of quantitative trait loci (QTL) controlling genetic predisposition for exercise-related traits, utilizing a large advanced intercross line (AIL) of mice. This AIL (G(4)) population originated from a reciprocal cross between mice with genetic propensity for increased voluntary exercise [high-runner (HR) line, selectively bred for increased wheel running] and the inbred strain C57BL/6J. After adjusting for family structure, we detected 32 significant and 13 suggestive QTL representing both daily running traits (distance, duration, average speed, and maximum speed) and the mean of these traits on days 5 and 6 (the selection criteria for HR) of a 6-day test conducted at 8 wk of age, with many co-localizing to similar genomic regions. Additionally, seven significant and five suggestive QTL were observed for the slope and intercept of a linear regression across all 6 days of running, some representing a combination of the daily traits. We also observed two significant and two suggestive QTL for body mass before exercise. These results, from a well-defined animal model, reinforce a genetic basis for the predisposition to engage in voluntary exercise, dissect this predisposition into daily segments across a continuous time period, and present unique QTL that may provide insight into the initiation, continuation, and temporal pattern of voluntary activity in mammals.

CC002/Unc females are mouse models of exercise-induced paradoxical fat response.

Exercise results in beneficial health outcomes and protects against a variety of chronic diseases. However, U.S. exercise guidelines recommend identical exercise programs for everyone, despite individual variation in responses to these programs, including paradoxical fat gain. Experimental models of exercise-induced paradoxical outcomes may enable the dissection of underlying physiological mechanisms as well as the evaluation of potential interventions. Whereas several studies have identified individual mice exhibiting paradoxical fat gain following exercise, no systematic effort has been conducted to identify and characterize models of paradoxical response. Strains from the Collaborative Cross (CC) genetic reference population were used due to its high levels of genetic variation, its reproducible nature, and the observation that the CC is a rich source of novel disease models, to assess the impact genetic background has on exercise responses. We identified the strain CC002/Unc as an exercise-induced paradoxical fat response model in a controlled voluntary exercise study across multiple ages in female mice. We also found sex and genetic differences were consistent with this pattern in a study of forced exercise programs. These results provide a novel model for studies to determine the mechanisms behind paradoxical metabolic responses to exercise, and enable development of more rational personalized exercise recommendations based on factors such as age, sex, and genetic background.

Developmental constraint through negative pleiotropy in the zygomatic arch.

BACKGROUND: Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology. RESULTS: We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2. CONCLUSIONS: A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.