Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

Systems genetics in diversity outbred mice inform BMD GWAS and identify determinants of bone strength.

Genome-wide association studies (GWASs) for osteoporotic traits have identified over 1000 associations; however, their impact has been limited by the difficulties of causal gene identification and a strict focus on bone mineral density (BMD). Here, we use Diversity Outbred (DO) mice to directly address these limitations by performing a systems genetics analysis of 55 complex skeletal phenotypes. We apply a network approach to cortical bone RNA-seq data to discover 66 genes likely to be causal for human BMD GWAS associations, including the genes SERTAD4 and GLT8D2. We also perform GWAS in the DO for a wide-range of bone traits and identify Qsox1 as a gene influencing cortical bone accrual and bone strength. In this work, we advance our understanding of the genetics of osteoporosis and highlight the ability of the mouse to inform human genetics.

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

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

Genetic Architecture Modulates Diet-Induced Hepatic mRNA and miRNA Expression Profiles in Diversity Outbred Mice.

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Lastly, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes.

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants 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 and facial 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.

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.

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.

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.

Improving Metabolic Health Through Precision Dietetics in Mice.

The incidence of diet-induced metabolic disease has soared over the last half-century, despite national efforts to improve health through universal dietary recommendations. Studies comparing dietary patterns of populations with health outcomes have historically provided the basis for healthy diet recommendations. However, evidence that population-level diet responses are reliable indicators of responses across individuals is lacking. This study investigated how genetic differences influence health responses to several popular diets in mice, which are similar to humans in genetic composition and the propensity to develop metabolic disease, but enable precise genetic and environmental control. We designed four human-comparable mouse diets that are representative of those eaten by historical human populations. Across four genetically distinct inbred mouse strains, we compared the American diet's impact on metabolic health to three alternative diets (Mediterranean, Japanese, and Maasai/ketogenic). Furthermore, we investigated metabolomic and epigenetic alterations associated with diet response. Health effects of the diets were highly dependent on genetic background, demonstrating that individualized diet strategies improve health outcomes in mice. If similar genetic-dependent diet responses exist in humans, then a personalized, or "precision dietetics," approach to dietary recommendations may yield better health outcomes than the traditional one-size-fits-all approach.

Biological/Genetic Regulation of Physical Activity Level: Consensus from GenBioPAC.

PURPOSE: Physical activity unquestionably maintains and improves health; however, physical activity levels globally are low and not rising despite all the resources devoted to this goal. Attention in both the research literature and the public policy domain has focused on social-behavioral factors; however, a growing body of literature suggests that biological determinants play a significant role in regulating physical activity levels. For instance, physical activity level, measured in various manners, has a genetic component in both humans and nonhuman animal models. This consensus article, developed as a result of an American College of Sports Medicine-sponsored round table, provides a brief review of the theoretical concepts and existing literature that supports a significant role of genetic and other biological factors in the regulation of physical activity. CONCLUSIONS: Future research on physical activity regulation should incorporate genetics and other biological determinants of physical activity instead of a sole reliance on social and other environmental determinants.

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.

Aerobic exercise prevents rarefaction of pial collaterals and increased stroke severity that occur with aging.

Variation in extent of the brain's collateral circulation is an important determinant of variation in the severity of stroke and efficacy of revascularization therapies. However, the number and diameter of pial collateral "arterioles" decrease with aging in associated with reduced eNOS and increased oxidative stress. We tested whether exercise reduces this aging-induced rarefaction. Twelve-month-old mice were randomized to sedentary or voluntary wheel-running. At 26 months' age, permanent MCA occlusion was followed 72 h later by determination of infarct volume and vascular casting after maximal dilation. The decline in collateral number and diameter and 2.4-fold increase in infarct volume evident in 26-versus 3-month-old sedentary mice were prevented by exercise-training. In contrast, number and diameter of the posterior communicating collateral "arteries" were unaffected by aging or exercise. Interestingly, diameter of the primary intracranial arteries increased with aging. Mechanistically, genetic overexpression of eNOS inhibited age-induced collateral rarefaction, and exercise increased eNOS and SOD2 and decreased the inflammatory marker NFkB assessed in hindlimb arteries. In conclusion, exercise prevented age-induced rarefaction of pial collaterals and reduced infarct volume. Aging also promoted outward remodeling of intracranial arteries. These effects were associated with increased eNOS and reduced markers of inflammation and aging in the vascular wall.

Prevention of tumorigenesis in mice by exercise is dependent on strain background and timing relative to carcinogen exposure.

Among cancer diagnoses, colorectal cancer (CRC) is prevalent, with a lifetime risk of developing CRC being approximately 5%. Population variation surrounding the mean risk of developing CRCs has been associated with both inter-individual differences in genomic architecture and environmental exposures. Decreased risk of CRC has been associated with physical activity, but protective responses are variable. Here, we utilized a series of experiments to examine the effects of genetic background (strain), voluntary exercise (wheel running), and their interaction on azoxymethane (AOM)-induced intestinal tumor number and size in mice. Additionally, we investigated how the timing of exercise relative to AOM exposure, and amount of exercise, affected tumor number and size. Our results indicated that voluntary exercise significantly reduced tumor number in a strain dependent manner. Additionally, among strains where exercise reduced tumor number (A/J, CC0001/Unc) the timing of voluntary exercise relative to AOM exposure was crucial. Voluntary exercise prior to or during AOM treatment resulted in a significant reduction in tumor number, but exercise following AOM exposure had no effect. The results indicate that voluntary exercise should be used as a preventative measure to reduce risk for environmentally induced CRC with the realization that the extent of protection may depend on genetic background.

Long-term exercise in mice has sex-dependent benefits on body composition and metabolism during aging.

Aging is associated with declining exercise and unhealthy changes in body composition. Exercise ameliorates certain adverse age-related physiological changes and protects against many chronic diseases. Despite these benefits, willingness to exercise and physiological responses to exercise vary widely, and long-term exercise and its benefits are difficult and costly to measure in humans. Furthermore, physiological effects of aging in humans are confounded with changes in lifestyle and environment. We used C57BL/6J mice to examine long-term patterns of exercise during aging and its physiological effects in a well-controlled environment. One-year-old male (n = 30) and female (n = 30) mice were divided into equal size cohorts and aged for an additional year. One cohort was given access to voluntary running wheels while another was denied exercise other than home cage movement. Body mass, composition, and metabolic traits were measured before, throughout, and after 1 year of treatment. Long-term exercise significantly prevented gains in body mass and body fat, while preventing loss of lean mass. We observed sex-dependent differences in body mass and composition trajectories during aging. Wheel running (distance, speed, duration) was greater in females than males and declined with age. We conclude that long-term exercise may serve as a preventive measure against age-related weight gain and body composition changes, and that mouse inbred strains can be used to characterize effects of long-term exercise and factors (e.g. sex, age) modulating these effects. These findings will facilitate studies on relationships between exercise and health in aging populations, including genetic predisposition and genotype-by-environment interactions.

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.

The 'Omics' of Voluntary Exercise: Systems Approaches to a Complex Phenotype.

Predisposition to engage in exercise is highly variable and simultaneously influenced by the environment, complex genomics, and their interactions. Given the importance of exercise to health, understanding the underlying influences of variability is crucial. Here, we discuss murine systems approaches, focusing on 'omics', relevant to revealing the architecture of voluntary activity.

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.

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.