Recent advances in optogenetics and gene therapy have led to promising new treatment strategies for blindness caused by retinal photoreceptor loss. Preclinical studies often rely on the retinal degeneration 1 (rd1 or Pde6b(rd1)) retinitis pigmentosa (RP) mouse model. The rd1 founder mutation is present in more than 100 actively used mouse lines. Since secondary genetic traits are well-known to modify the phenotypic progression of photoreceptor degeneration in animal models and human patients with RP, negligence of the genetic background in the rd1 mouse model is unwarranted. Moreover, the success of various potential therapies, including optogenetic gene therapy and prosthetic implants, depends on the progress of retinal degeneration, which might differ between rd1 mice. To examine the prospect of phenotypic expressivity in the rd1 mouse model, we compared the progress of retinal degeneration in two common rd1 lines, C3H/HeOu and FVB/N.

Vesico-ureteric reflux is the most common congenital anomaly of the urinary tract, characterized by a defective uretero-vesical junction with retrograde urine flow from the bladder toward the kidneys. Because there is strong evidence for a genetic basis for some cases of vesico-ureteric reflux, we screened 11 inbred mouse strains for reflux and kidney size and identified one strain, C3H/HeJ, that has a 100 percent incidence of vesico-ureteric reflux with otherwise normal kidneys at birth. These mice are predisposed to reflux as a result of a defective uretero-vesical junction characterized by a short intravesical ureter. This defect results from a delay in urinary tract development initially manifested by a ureteric bud arising from a more caudal location along the mesonephric duct. In contrast, C57BL/6J mice (resistant to reflux at birth) have long intravesical ureters, normally positioned ureteric buds, and no delay in urinary tract development. Genome-wide and additional fine mapping of backcross mice, derived from C3H/HeJ and C57BL/6J crosses, identified a significant reflux susceptibility locus, Vurm1, on chromosome 12 (peak logarithm of the odds=7.39). The C3H/HeJ mouse is a model of vesico-ureteric reflux without renal malformation, and further characterization of this model will allow for the identification of a pathway important for urinary tract development, a finding that will serve as a model for the human disorder.

Insertions of endogenous retroviruses cause a significant fraction of mutations in inbred mice but not all strains are equally susceptible. Notably, most new Intracisternal A particle (IAP) ERV mutagenic insertions have occurred in C3H mice. We show here that strain-specific insertional polymorphic IAPs accumulate faster in C3H/HeJ mice, relative to other sequenced strains, and that IAP transcript levels are higher in C3H/HeJ embryonic stem (ES) cells compared to other ES cells. To investigate the mechanism for high IAP activity in C3H mice, we identified 61 IAP copies in C3H/HeJ ES cells enriched with H3K4me3 (a mark of active promoters) and, among those tested, all are unmethylated in C3H/HeJ ES cells. Notably, 13 of the 61 are specific to C3H/HeJ and are members of the non-autonomous 1Δ1 IAP subfamily that is responsible for nearly all new insertions in C3H. One copy is full length with intact open reading frames and hence potentially capable of providing proteins in trans to other 1Δ1 elements. This potential "master copy" is present in other strains, including 129, but its 5' long terminal repeat (LTR) is methylated in 129 ES cells. Thus, the unusual IAP activity in C3H may be due to reduced epigenetic repression coupled with the presence of a master copy.

Variation at regulatory elements, identified through hypersensitivity to digestion by DNase I, is believed to contribute to variation in complex traits, but the extent and consequences of this variation are poorly characterized. Analysis of terminally differentiated erythroblasts in eight inbred strains of mice identified reproducible variation at approximately 6% of DNase I hypersensitive sites (DHS). Only 30% of such variable DHS contain a sequence variant predictive of site variation. Nevertheless, sequence variants within variable DHS are more likely to be associated with complex traits than those in non-variant DHS, and variants associated with complex traits preferentially occur in variable DHS. Changes at a small proportion (less than 10%) of variable DHS are associated with changes in nearby transcriptional activity. Our results show that whilst DNA sequence variation is not the major determinant of variation in open chromatin, where such variants exist they are likely to be causal for complex traits.

The genetic factors underlying female infertility in humans are only partially understood. Here, we performed a genome-wide association study of female infertility in 25 inbred mouse strains by using publicly available SNP data. As a result, a total of four SNPs were identified after chromosome-wise multiple test correction. The first SNP rs29972765 is located in a gene desert on chromosome 18, about 72 kb upstream of Skor2 (SKI family transcriptional corepressor 2). The second SNP rs30415957 resides in the intron of Plce1 (phospholipase C epsilon 1). The remaining two SNPs (rs30768258 and rs31216810) are close to each other on chromosome 19, in the vicinity of Sorbs1 (sorbin and SH3 domain containing 1). Using quantitative RT-PCR, we found that Sorbs1 is highly expressed in the mouse uterus during embryo implantation. Knockdown of Sorbs1 by siRNA attenuates the induction of differentiation marker gene Prl8a2 (decidual prolactin-related protein) in an in vitro model of decidualization using mouse endometrial stromal cells, suggesting that Sorbs1 may be a potential candidate gene for female infertility in mice. Our results may represent an opportunity to further understand female infertility in humans.

Genetic make-up strongly influences the skeleton's susceptibility to the loss of weight bearing with some inbred mouse strains experiencing great amounts of bone loss while others lose bone at much smaller rates. At young adulthood, female inbred C3H/HeJ (C3H) mice are largely resistant to catabolic pressure induced by unloading. Here, we tested whether the depressed responsivity to unloading is inherent to the C3H genetic make-up or whether a younger age facilitates a robust skeletal response to unloading. Nine-week-old, skeletally immature, female C3H mice were subjected to 3wk of hindlimb unloading (HLU, n = 12) or served as normal baseline controls (BC, n = 10) or age-matched controls (AC, n = 12). In all mice, cortical and trabecular architecture of the femur, as well as levels of bone formation and resorption, were assessed with μCT, histomorphometry, and histology. Changes in bone marrow progenitor cell populations were determined with flow cytometry. Following 21d of unloading, HLU mice had 52% less trabecular bone in the distal femur than normal age-matched controls. Reflecting a loss of trabecular tissue compared to baseline controls, trabecular bone formation rates (BFR/BS) in HLU mice were 40% lower than in age-matched controls. Surfaces undergoing osteoclastic resorption were not significantly different between groups. In the mid-diaphysis, HLU inhibited cortical bone growth leading to 14% less bone area compared to age-matched controls. Compared to AC, BFR/BS of HLU mice were 53% lower at the endo-cortical surface and 49% lower at the periosteal surface of the mid-diaphysis. The enriched osteoprogenitor cell population (OPC) comprised 2% of the bone marrow stem cells in HLU mice, significantly different from 3% OPC in the AC group. These data show that bone tissue in actively growing C3H mice is lost rapidly, or fails to grow, during the removal of functional weight bearing-in contrast to the insignificant response previously demonstrated in female young adult C3H mice. Thus, the attributed low sensitivity of the C3H mouse strain to the loss of mechanical signals is not apparent at a young age and this trait therefore does not reflect a genetic regulation throughout the life span of this strain. These results highlight the significance of age in modulating the contribution of genetics in orchestrating bone's response to unloading and that the skeletal unresponsiveness of young adult C3H mice to the loss of weight bearing is not genetically hard-wired.

The interfrontal bone (IF) is a minor skeletal trait residing between the frontal bones. IF is considered a quasi-continuous trait. Genetic and environmental factors appear to play roles in its development. The mechanism(s) underlying IF bone development are poorly understood. We sought to survey inbred strains of mice for the prevalence of IF and to perform QTL mapping studies. Archived mouse skulls from a mouse phenome project (MPP) were available for this study. 27 inbred strains were investigated with 6-20 mice examined for each strain. Skulls were viewed dorsally and the IF measured using a zoom stereomicroscope equipped with a calibrated reticle. A two generation cross between C3H/HeJ and C57BL/6J mice was performed to generate a panel of 468 F2 mice. F2 mice were phenotyped for presence or absence of IF bone and among mice with the IF bone maximum widths and lengths were measured. F2 mice were genotyped for 573 SNP markers informative between the two strains and subjected to linkage map construction and interval QTL mapping. Results: Strain dependent differences in the prevalence of IF bones were observed. Overall, 77.8% or 21/27, of the inbred strains examined had IF bones. Six strains (C3H/HeJ, MOLF/EiJ, NZW/LacJ, SPRET/EiJ, SWR/J, and WSB/EiJ) lack IF bones. Among the strains with IF bones, the prevalence ranged from 100% for C57BL/6J, C57/LJ, CBA/J, and NZB/B1NJ and down to 5% for strains such as CAST/Ei. QTL mapping for IF bone length and widths identifies for each trait one strong QTL detected on chromosome 14 along with several other significant QTLs on chromosomes 3, 4, 7, and 11. Strain dependent differences in IF will facilitate investigation of genetic factors contributing to IF development. IF bone formation may be a model to understand intrasutural bone formation.

To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis.

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Numerous candidate genes have been suggested in the recent literature with proposed roles in regulation of voluntary physical activity, with little evidence of these genes' functional roles. This study compared the haplotype structure and expression profile in skeletal muscle and brain of inherently high- (C57L/J) and low- (C3H/HeJ) active mice. Expression of nine candidate genes [Actn2, Actn3, Casq1, Drd2, Lepr, Mc4r, Mstn, Papss2, and Glut4 (a.k.a. Slc2a4)] was evaluated via RT-qPCR. SNPs were observed in regions of Actn2, Casq1, Drd2, Lepr, and Papss2; however, no SNPs were located in coding sequences or associated with any known regulatory sequences. In mice exposed to a running wheel, Casq1 (P = 0.0003) and Mstn (P = 0.002) transcript levels in the soleus were higher in the low-active mice. However, when these genes were evaluated in naïve animals, differential expression was not observed, demonstrating a training effect. Among naïve mice, no genes in either tissue exhibited differential expression between strains. Considering that no obvious SNP mechanisms were determined or differential expression was observed, our results indicate that genomic structural variation or gene expression data alone is not adequate to establish any of these genes' candidacy or causality in relation to regulation of physical activity.

For over a century, inbred mice have been used in many areas of genetics research to gain insight into the genetic variation underlying traits of interest. The generalizability of any genetic research study in inbred mice is dependent upon all individual mice being genetically identical, which in turn is dependent on the breeding designs of companies that supply inbred mice to researchers. Here, we compare whole-genome sequences from individuals of four commonly used inbred strains that were procured from either the colony nucleus or from a production colony (which can be as many as ten generations removed from the nucleus) of a large commercial breeder, in order to investigate the extent and nature of genetic variation within and between individuals. We found that individuals within strains are not isogenic, and there are differences in the levels of genetic variation that are explained by differences in the genetic distance from the colony nucleus. In addition, we employ a novel approach to mutation rate estimation based on the observed genetic variation and the expected site frequency spectrum at equilibrium, given a fully inbred breeding design. We find that it provides a reasonable per nucleotide mutation rate estimate when mice come from the colony nucleus (~7.9 × 10-9 in C3H/HeN), but substantially inflated estimates when mice come from production colonies.

Vomeronasal receptors (VRs), expressed in sensory neurons of the vomeronasal organ, are thought to bind pheromones and mediate innate behaviours. The mouse reference genome has over 360 functional VRs arranged in highly homologous clusters, but the vast majority are of unknown function. Differences in these receptors within and between closely related species of mice are likely to underpin a range of behavioural responses. To investigate these differences, we interrogated the VR gene repertoire from 17 inbred strains of mice using massively parallel sequencing.

Absence epilepsy, characterized by spike-wave discharges (SWD) in the electroencephalogram, arises from aberrations within the circuitry of the cerebral cortex and thalamus that regulates awareness. The inbred mouse strain C3H/HeJ is prone to absence seizures, with a major susceptibility locus, spkw1, accounting for most of the phenotype. Here we find that spkw1 is associated with a hypomorphic retroviral-like insertion mutation in the Gria4 gene, encoding one of the four amino-3-hydroxy-5-methyl-4isoxazolepropionic acid (AMPA) receptor subunits in the brain. Consistent with this, Gria4 knockout mice also have frequent SWD and do not complement spkw1. In contrast, null mutants for the related gene Gria3 do not have SWD, and Gria3 loss actually lowers SWD of spkw1 homozygotes. Gria3 and Gria4 encode the predominant AMPA receptor subunits in the reticular thalamus, which is thought to play a central role in seizure genesis by inhibiting thalamic relay cells and promoting rebound burst firing responses. In Gria4 mutants, synaptic excitation of inhibitory reticular thalamic neurons is enhanced, with increased duration of synaptic responses-consistent with what might be expected from reduction of the kinetically faster subunit of AMPA receptors encoded by Gria4. These results demonstrate for the first time an essential role for Gria4 in the brain, and suggest that abnormal AMPA receptor-dependent synaptic activity can be involved in the network hypersynchrony that underlies absence seizures.

Despite established health benefits of regular exercise, the majority of Americans do not meet the recommended levels of physical activity. While it is known that voluntary activity levels are largely heritable, the genetic mechanisms that regulate activity are not well understood. MicroRNAs (miRNAs) are small non-coding RNAs that inhibit transcription by binding to a target gene, inhibiting protein production. The purpose of this study was to investigate differential miRNA expression between inherently high- (C57L/J) and low- (C3H/HeJ) active inbred mice in soleus, extensor digitorum longus (EDL), and nucleus accumbens tissues. Expression was initially determined by miRNA microarray analysis, and selected miRNAs were validated by qRT-PCR. Expression of 13 miRNAs varied between strains in the nucleus accumbens, 20 in soleus, and eight in EDL, by microarray analysis. Two miRNAs were validated by qRT-PCR in the nucleus accumbens; miR-466 was downregulated (~4 fold; P < 0.0004), and miR-342-5p was upregulated (~115 fold; P < 0.0001) in high-active mice. MiR-466 was downregulated (~5 fold; P < 0.0001) in the soleus of high-active mice as well. Interestingly, miR-466 is one of several miRNA families with sequence located in intron 10 of Sfmbt2; miRNAs at this locus are thought to drive imprinting of this gene. "Pathways in cancer" and "TGFβ signaling" were the most significant pathways of putative target genes in both the soleus and nucleus accumbens. Our results are the first to consider differential miRNA expression between high- and low-active mice, and suggest that miRNAs may play a role in regulation of physical activity.

Insulin resistance plays a central role in type 2 diabetes and obesity, which develop as a consequence of genetic and environmental factors. Dietary changes including high fat diet (HFD) feeding promotes insulin resistance in rodent models which present useful systems for studying interactions between genetic background and environmental influences contributing to disease susceptibility and progression. We applied a combination of classical physiological, biochemical and hormonal studies and plasma (1)H NMR spectroscopy-based metabonomics to characterize the phenotypic and metabotypic consequences of HFD (40%) feeding in inbred mouse strains (C57BL/6, 129S6, BALB/c, DBA/2, C3H) frequently used in genetic studies. We showed the wide range of phenotypic and metabonomic adaptations to HFD across the five strains and the increased nutrigenomic predisposition of 129S6 and C57BL/6 to insulin resistance and obesity relative to the other strains. In contrast mice of the BALB/c and DBA/2 strains showed relative resistance to HFD-induced glucose intolerance and obesity. Hierarchical metabonomic clustering derived from (1)H NMR spectral data of the strains provided a phylometabonomic classification of strain-specific metabolic features and differential responses to HFD which closely match SNP-based phylogenetic relationships between strains. Our results support the concept of genomic clustering of functionally related genes and provide important information for defining biological markers predicting spontaneous susceptibility to insulin resistance and pathological adaptations to fat feeding.

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

The expanding set of genomics tools available for inbred mouse strains has renewed interest in phenotyping larger sets of strains. The present study aims to explore phenotypic variability among six commonly-used inbred mouse strains to both the rewarding and locomotor stimulating effects of cocaine in a place conditioning task, including several strains or substrains that have not yet been characterized for some or all of these behaviors.

Many strains of mice are utilized in mouse models of cerebrovascular diseases. Variations in vascular anatomy between these strains has been documented and may influence the phenotype in stroke models. To address inter-strain variations in the circle of Willis anatomy, the diameters of internal carotid, posterior communicating, anterior cerebral, and middle cerebral arteries in 144 mice from 32 inbred strains were measured. Arterial diameters were analyzed as a function of animal weight, age, and strain. Variations in the structure of the circle of Willis across strains were observed and noted. While right-sided anterior cerebral arteries were significantly greater in diameter than their left-sided counterparts across most strains, variations in arterial diameter are strain specific. Adult mouse weight was not found to be associated with arterial diameter across strains, suggesting that cerebral artery size is associated with strain independently of weight. This study demonstrates strain dependent variations in the murine circle of Willis, which should be taken into consideration when studying mouse models of cerebrovascular diseases.

Inbred strains of mice vary in their frequency of liver tumors initiated by a mutation in the Hras1 (H-ras) proto-oncogene. We sequenced 4.5 kb of the Hras1 gene on distal Chr 7 in a diverse set of 12 commonly used laboratory inbred strains of mice and detected no sequence variation to account for strain-specific differences in Hras1 mutation prevalence. Furthermore, the Hras1 sequence is essentially monoallelic for an ancestral gene derived from the M. m. domesticus species. To determine if the monoallelism and associated low rate of polymorphism are unique to Hras1 or representative of the general chromosomal locale, we extended the sequence analysis to 12 genes in the final 8 Mb of distal Chr 7. A region of at least 2.5 Mb that encompasses several genes, including Hras1 and the H19/Igf2 loci, demonstrates virtually no sequence variation. The 12 inbred strains share one dominant haplotype derived from the M. m. domesticus allele. Chromosomal regions flanking the monoallelic segment exhibit a significantly higher rate of variation and multiple haplotypes, a majority of which are attributed to M. m. domesticus or M. m. musculus ancestry.

To investigate the role of complement C3 (C3) convertase on the strain difference for C3 protein expression in three inbred mice strains, we compared the levels of C2, C3 and C4 mRNA, as well as C3 protein and C3 convertase activity in the serum and liver tissue of FVB/N, C3H/HeN and C57BL/6N mice. The level of mRNA, inactive form (InACF) and active form (ACF) for C3 showed a regular pattern, which they were higher in the FVB/N and C57BL/6N mice than C3H/HeN mice. However, the level of C3b fragments (C3bα and β) derived from C3 protein were constantly maintained in the liver of FVB/N, C3H/HeN and C57BL/6N mice in spite of the strain difference on the transcriptional and translation level of C3. Especially, a reverse pattern of the level of mRNA, InACF and ACF for C3 was observed on the activity level of C3 convertase activity. The highest level of C3 convertase activity was measured in C3H/HeN mice, followed by C57BL/6N and FVB/N mice. In case of C3 convertase components, the level of C2 mRNA was higher in C3H/HeN mice than FVB/N and C57BL/6 N mice, while levels of C4 mRNA were higher in FVB/N and C57BL/6N mice than C3H/HeN mice. The current study results provide the first scientific evidence that C3 convertase may play complementary role to overcome the strain difference on the C3 protein expression in FVB/N, C3H/HeN and C57BL/6N mice.