Although mice are the most widely used mammalian model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of ...crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay in comparison to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping by sequencing (GBS) to obtain genotypes at 92,734 SNPs. We also measured gene expression using RNA sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA sequencing constitutes a powerful approach to GWAS in mice.
Two bottlenecks impeding the genetic analysis of complex traits in rodents are access to mapping populations able to deliver gene-level mapping resolution and the need for population-specific ...genotyping arrays and haplotype reference panels. Here we combine low-coverage (0.15×) sequencing with a new method to impute the ancestral haplotype space in 1,887 commercially available outbred mice. We mapped 156 unique quantitative trait loci for 92 phenotypes at a 5% false discovery rate. Gene-level mapping resolution was achieved at about one-fifth of the loci, implicating Unc13c and Pgc1a at loci for the quality of sleep, Adarb2 for home cage activity, Rtkn2 for intensity of reaction to startle, Bmp2 for wound healing, Il15 and Id2 for several T cell measures and Prkca for bone mineral content. These findings have implications for diverse areas of mammalian biology and demonstrate how genome-wide association studies can be extended via low-coverage sequencing to species with highly recombinant outbred populations.
We examine some of the genetic features of neuroticism (N) taking as an animal model the Maudsley Reactive (MR) and Maudsley Nonreactive (MNR) rat strains which were selectively bred, respectively, ...for high and low open-field defecation (OFD) starting in the late 1950s. To draw analogies with human genetic studies, we explore the genetic correlation of N with irritable bowel syndrome (IBS). We review progress with the rat model and developments in the field of human complex trait genetics, including genetic association studies that relate to current understanding of the genetics of N. The widespread differences in the tone of the peripheral sympathetic nervous system that have been found between the Maudsley strains, particularly those observed in the colon, may underly the differences in OFD (MNR, higher sympathetic tone and zero defecation). In humans, a large genome-wide association study (GWAS) reported six genes contributing to IBS, four of which were implicated in mood and anxiety disorders or were expressed in the brain, with three of the four also expressed in the nerve fibers and ganglia of the gut. Heritability of N is estimated at around 50% in twin and family studies, and GWASs identified hundreds of loci, enabling estimation of genome-wide correlations (r
) with other traits. Significantly, the estimate for r
between risk of IBS, anxiety, N, and depression was >0.5 and suggested genetic pleiotropy without evidence for causal mechanisms. Findings on the adrenergic pharmacology of the colon, coupled with new understanding of the role of the locus ceruleus in modifying afferent information from this organ, generate hypotheses that challenge traditional cause/effect notions about the relationship of the central nervous system to peripheral events in response to stress, suggest specific targets for gene action in the Maudsley model and emphasize the value of reciprocal evaluation of genetic architecture underlying N in rodents and humans.
We compared the rate of acquisition and strength of retention of conditioned context aversion (CCA) with conditioned taste aversion (CTA) using pigmented, genetically heterogeneous mice (derived from ...Large and Small strains). Extending previous findings, in Experiment
1
, mice accustomed to drinking from large glass bottles in the colony room learned to avoid graduated tubes after a single conditioning trial when drinking from these novel tubes was paired with injections of LiCl. The results also showed that CCA could be developed even when there was a 30-minute delay between conditioned stimulus and unconditioned stimulus. Retention of the aversion lasted for 4 weeks in both Immediate and Delay groups. Studies of conditioned saccharin aversion were conducted in Experiment
2
. CTA acquisition was very similar to that observed in CCA and duration of aversion retention was similar in the CCA and CTA Delay groups, although at least 2 weeks longer in the Immediate group. Thus, CCA acquisition and retention characteristics are closer to those seen for CTA than has previously been reported. In Experiment
3
, we examined whether albino mice (which are known to have weaker visual abilities compared to pigmented mice) would develop CCA comparable to those of pigmented mice. The development of conditioned aversion and its duration of retention was similar in albinos and pigmented mice. Nonspecific aversion emerged as an important contributor to strength of aversion during retention trials in both CCA and CTA paradigms with pigmented (but not albino) mice and deserves additional scrutiny in this field of inquiry.
It is well known that pairing of large contextual changes with illness can cause conditioned context aversion in laboratory rats. The aim of present study was to develop a paradigm to study this ...phenomenon in laboratory mice, a species widely employed in neurobehavioral studies. Genetically heterogeneous mice, drinking from plastic bottles in the colony room, learned to avoid glass bottles after a single conditioning trial when drinking from these was paired with injections of lithium chloride. The aversion was independent of any difference in the taste of water in plastic vs. glass bottles. When the variation in the visual stimulus was less distinct, development of a strong aversion required two conditioning trials and was not retained as well. The results also showed that conditioned context aversion, just like conditioned taste aversion, could also be developed across a 30–minute CS–UCS delay. The fact that taste was not a factor in distinguishing drinking from glass and plastic water bottles raises the possibility that, contextual stimuli, not taste, may have been the CS when rats (in Garcia’s original experiments) avoided drinking from plastic bottles that had been paired with radiation. The development of contextual aversion conditioning protocols for mice will enable the molecular resources available for this species to be exploited. Furthermore, representation of the CS by discrete rather than the multimodal CSs typically used in most studies on contextual conditioning offers more focus when considering its neuroanatomical basis.
Studies of inbred strains of rats and mice have suggested a positive association between strain variations in sweet taste and ethanol intake. However, strain associations by themselves are ...insufficient to support a functional link between taste and ethanol intake. We used conditioned taste aversion (CTA) to explore the sweet and bitter taste of ethanol and ability to detect sucrose, quinine and ethanol in C57BL/6J (B6) and DBA/2J (D2) mouse strains that are frequently used in alcohol research. The present study showed that C57BL/6J mice generalized taste aversions from sucrose and quinine solutions to 10% ethanol and, reciprocally, aversions to 10% ethanol generalized to each of these solutions presented separately. Only conditioned aversions to quinine generalized to ethanol in the DBA/2J strain but an aversion conditioned to ethanol did not generalize reciprocally to quinine. Thus, considering these two gustatory qualities, 10% ethanol tastes both sweet and bitter to B6 mice but only bitter to D2. Both strains were able to generalize taste aversions across different concentrations of the same compound. B6 were able to detect lower concentrations of quinine than D2 but both strains were able to detect sucrose and (in contrast to previous findings) ethanol at similar concentrations. The strain-dependent gustatory profiles for ethanol may make an important contribution to the understanding of the undoubtedly complex mechanisms influencing high ethanol preference of B6 and pronounced ethanol avoidance of D2 mice.
We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence ...intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.
13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).
Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.
The Maudsley Reactive and Maudsley Non-Reactive strains have been selectively bred for differences in open-field defecation (OFD), a putative index of stress. We investigated whether variations in ...the hypothalamic–pituitary–adrenal (HPA) axis are correlated with strain differences in OFD in the Maudsley model. Exposure to the open-field test did not result in increases in ACTH in male rats of either strain and there were no strain differences in the large increases in ACTH and corticosteroid that occurred in response to intermittent footshock. Parallel studies of prolactin showed that Maudsley Reactive rats had greater response to the open-field and to footshock than Maudsley Non-Reactive rats. The lack of correlation between strain differences in OFD and reactivity of the HPA axis is consistent with the idea that HPA response to stress and OFD reflect the output of different neural systems and that individual differences in emotionality, as indexed by OFD do not influence other measures of stress-reactivity in a simple manner, if at all. The reactivity of the prolactin system to the open-field test and lack of response of ACTH to the same situation is consistent with the idea that the prolactin system is sensitive to lower levels of stress than the HPA axis, a finding at variance with the presumed extreme sensitivity of the latter system. Earlier comparisons of the HPA axis in these strains implicate local factors such as neuropeptide-Y peptide in the adrenal in attenuating the response of the adrenal cortex to ACTH and hints at the complexity of regulation of the HPA axis.
Berlin high (BEH) and Berlin low (BEL) strains selected for divergent growth differ threefold in body weight. We aimed at examining muscle mass, which is a major contributor to body weight, by ...exploring morphological characteristics of the soleus muscle (fiber number and cross sectional area; CSA), by analyzing the transcriptome of the gastrocnemius and by initiating quantitative trait locus (QTL) mapping. BEH muscles were four to eight times larger than those of BEL. In substrain BEH+/+, mutant myostatin was replaced with a wild-type allele; however, BEH+/+muscles still were two to four times larger compared with BEL. BEH soleus muscle fibers were two times more numerous (P < 0.0001) and CSA was two times larger (P < 0.0001) compared with BEL. In addition, soleus femoral attachment anomaly (SFAA) was observed in all BEL mice. One significant (Chr 1) and four suggestive (Chr 3, 4, 6, and 9) muscle weight QTLs were mapped in a 21-day-old F2 intercross (n = 296) between BEH and BEL strains. The frequency of SFAA incidence in the F2 and in the backcross to BEL strain (BCL) suggested the presence of more than one causative gene. Two suggestive SFAA QTLs were mapped in BCL; however, their peak markers were not associated with the phenotype in F2. RNA-Seq analysis revealed 2,148 differentially expressed (P < 0.1) genes and 45,673 single nucleotide polymorphisms and >2,000 indels between BEH+/+ and BEL males. In conclusion, contrasting muscle traits and genomic and gene expression differences between BEH and BEL strains provide a promising model for the search for genes involved in muscle growth and musculoskeletal morphogenesis.
Correlations among bone strength, muscle mass, and physical activity suggest that these traits may be modulated by each other and/or by common genetic and/or environmental mechanisms. This study used ...structural equation modeling (SEM) to explore the extent to which select genetic loci manifest their pleiotropic effects through functional adaptations commonly referred to as Wolff's law. Quantitative trait locus (QTL) analysis was used to identify regions of chromosomes that simultaneously influenced skeletal mechanics, muscle mass, and/or activity‐related behaviors in young and aged B6×D2 second‐generation (F2) mice of both sexes. SEM was used to further study relationships among select QTLs, bone mechanics, muscle mass, and measures of activity. The SEM approach provided the means to numerically decouple the musculoskeletal effects of mechanical loading from the effects of other physiological processes involved in locomotion and physical activity. It was found that muscle mass was a better predictor of bone mechanics in young females, whereas mechanical loading was a better predictor of bone mechanics in older females. An activity‐induced loading factor positively predicted the mechanical behavior of hindlimb bones in older males; contrarily, load‐free locomotion (i.e., the remaining effects after removing the effects of loading) negatively predicted bone performance. QTLs on chromosomes 4, 7, and 9 seem to exert some of their influence on bone through actions consistent with Wolff's Law. Further exploration of these and other mechanisms through which genes function will aid in development of individualized interventions able to exploit the numerous complex pathways contributing to skeletal health.