A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains

Background

The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms.

Results

We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems.

Conclusions

Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.     

Morphometry and strain distribution of the C57BL/6 mouse aorta

The goal of the present study was to obtain a systematic set of data along the length of the mouse aorta to study variations of morphometry (diameter, wall thickness, and curvature), strain, and stress of the mouse aorta. Five mice were imaged with a 13-MHz ultrasound probe to determine the in vivo diameter along the aorta. A cast was made of these aortas to validate the ultrasonic diameter measurements. The root mean squared and systematic errors for these measurements were 12.6% and 6.4% of the mean ultrasound diameter, respectively. The longitudinal variations of geometry, stress, and strain from the aortic valve to the common iliac bifurcation were documented. Our results show that the residual circumferential strain leads to a uniformity of transmural strain of the aorta in the loaded state along the entire length of the aorta. Furthermore, we validated the incompressibility condition along the length of the aorta. These data of normal mice will serve as a reference state for the study of disease in future knockout models.     

Lipidomic analysis and electron transport chain activities in C57BL/6J mouse brain mitochondria

The objective of this study was to characterize the lipidome and electron transport chain activities in purified non-synaptic (NS) and synaptic (Syn) mitochondria from C57BL/6J mouse cerebral cortex. Contamination from subcellular membranes, especially myelin, has hindered past attempts to accurately characterize the lipid composition of brain mitochondria. An improved Ficoll and sucrose discontinuous gradient method was employed that yielded highly enriched mitochondrial populations free of myelin contamination. The activities of Complexes I, II, III, and II/III were lower in Syn than in NS mitochondria, while Complexes I/III and IV activities were similar in both populations. Shotgun lipidomics showed that levels of cardiolipin (Ptd₂Gro) were lower, whereas levels of ceramide and phosphatidylserine were higher in Syn than in NS mitochondria. Coenzyme Q₉ and Q₁₀ was also lower in Syn than in NS mitochondria. Gangliosides, phosphatidic acid, sulfatides, and cerebrosides were undetectable in brain mitochondria. The distribution of Ptd₂Gro molecular species was similar in both populations and formed a unique pattern, consisting of seven major molecular species groups, when arranged according to mass to charge ratios. Remodeling involving choline and ethanolamine phosphoglycerides could explain Ptd₂Gro heterogeneity. NS and Syn mitochondrial lipidomic heterogeneity could influence energy metabolism, which may contribute to metabolic compartmentation of the brain.     

Genetic polymorphisms among C57BL/6 mouse inbred strains

Mice from the inbred C57BL/6 strain have been commonly used for the generation and analysis of transgenic and knockout animal models. However, several C57BL/6 substrains exist, and these are genetically and phenotypically different. In addition, each of these substrains can be purchased from different animal providers and, in some cases, they have maintained their breeding stocks separated for a long time, allowing genetic differences to accumulate due to individual variability and genetic drift. With the aim of describing the differences in the genotype of several C57BL/6 substrains, we applied the Illumina^® Mouse Medium Density Linkage Mapping panel, with 1,449 single nucleotide polymorphisms (SNPs), to individuals from ten C57BL/6-related strains: C57BL/6JArc, C57BL/6J from The Jackson Lab, C57BL/6J from Crl, C57BL6/JRccHsd, C57BL/6JOlaHsd, C57BL/6JBomTac, B6(Cg)-Tyr ^c?2j /J, C57BL/6NCrl, C57BL/6NHsd and C57BL/6NTac. Twelve SNPs were found informative to discriminate among the mouse strains considered. Mice derived from the original C57BL/6J: C57BL/6JArc, C57BL/6J from The Jackson Lab and C57BL/6J from Crl, were indistinguishable. Similarly, all C57BL/6N substrains displayed the same genotype, whereas the additional substrains showed intermediate cases with substrain-specific polymorphisms. These results will be instrumental for the correct genetic monitoring and appropriate mouse colony handling of different transgenic and knockout mice produced in distinct C57BL/6 inbred substrains.     

Embryonic stem cells derived from C57BL/6J and C57BL/6N mice

Mouse embryonic stem (ES) cells with the C57BL/6 genetic background allow the generation of knockout mice without the need to backcross to C57BL/6. However, C57BL/6 ES cells whose pluripotency after homologous recombination has been confirmed are not yet available from public cell banks. To facilitate the use of ES cells derived from C57BL/6 sublines in both biologic and medical research, we demonstrated that the use of knockout serum replacement as a medium supplement and 8-cell blastomeres as recipient embryos allowed establishment of ES cells and production of germline chimeric mice, respectively. Under effective conditions, a large number of ES cell lines were established from C57BL/6J and C57BL/6N blastocysts. The majority of ES cells in many cell lines obtained from both strains showed a normal chromosome number. Germline chimeric mice were generated from C57BL/6J and C57BL/6N ES cells. Finally, the ES cell line B6J-S1UTR, derived from C57BL/6J, was used for successful production of gene knockout mice. C57BL/6J ES (B6J-S1UTR and B6J-23UTR) and C57BL/6N ES (B6N-22UTR) cells are available from the cell bank of the BioResource Center at RIKEN Tsukuba Institute (http://www.brc.riken.jp/lab/cell/english/).     

Nuclear proteome profile of C57BL/6J mouse liver

The liver proteome can serve as a reference to better understand both disease mechanisms and possible therapeutics,since the liver is an important organ in the body that performs a large number of tasks.Here we identify the organelle proteome of C57BL/6J mouse liver nuclei as a promising strategy to enrich low abundance proteins,in the sense that analysis of whole liver cells is rather complex for current techniques and may not be suitable for proteins with low abundance.Evaluation of nucleus integrity and purity was performed to demonstrate the effectiveness of the optimized isolation procedure.The extracted nuclear proteins were identified by 2-DE MS analyses,and a total of 748 proteins were identified.Bioinformatic analyses were performed to demonstrate the physicochemical properties,cellular locations and functions of the proteins.     

Pavlovian biconditional discrimination learning in the C57BL/6J mouse

Four experiments using mice examined acquisition of Pavlovian biconditional discriminations in which two stimulus compounds were paired with food (AX+ and BY+) and two were not (AY- and BX-). Temporally asynchronous compounds were generated by using contextual stimuli (Experiment 1) and 15-s discrete visual cues (Experiments 2A, 2B and 3) to disambiguate when embedded noise or tone stimuli would be paired with food. When food pellets followed both reinforced compounds, successful acquisition was obtained in Experiment 1 but not in Experiments 2A and 2B even though the order of trials was modeled after that used in Experiment 1. However, when differential outcomes followed the reinforced compounds in Experiment 3, acquisition was obtained with discrete cue stimulus compounds. The implications of these results for modulatory models of conditional discrimination learning in animals are discussed.     

Genetic differentiation within the inbred C57BL/6J mouse strain

The C57BL/6J (B6) inbred mouse strain is commonly used in biomedical researches. However, some unexpected inconsistency was reported compared with previous studies, and in most cases, it can be attributed to environmental, epigenetic or stochastic differences. The goal of this study was to investigate the genetic stability of the B6 strain maintained in different breeders. B6 mice purchased from five Chinese commercial breeders were examined, and mitochondrial D-loop sequence and 18 microsatellite loci were genotyped. There is no difference in the D-loop sequences, but variations exist in the nucleic microsatellite markers. Combining the data from MGI_4.01, a significant divergence is observed among those mice. The present study indicates that different B6 mice share the common maternal lineage and are still inbred in each breeder, but subline divergence occurs.     

Purification of the Ah receptor from C57BL/6J mouse liver

The photoaffinity ligand for the Ah receptor, [125I]-2-azido-3-iodo-7,8-dibromodibenzo-p-dioxin, previously has been shown to selectively label two peptides in the cytosol fraction of C57BL/6J mouse liver: a 95-kDa peptide, the ligand binding moiety of the Ah receptor, and a 70-kDa proteolytic fragment formed from the larger peptide (Poland, A., Glover, E., Ebetino, F. H., and Kende, A.S. (1986) J. Biol. Chem. 261, 6352-6365). These two peptides were partially purified to an approximately 20,000-fold enrichment with a 15-20% yield by the following scheme: 1) photoaffinity labeling of the 35-55% ammonium sulfate fraction of liver cytosol; 2) chromatography on polyethyleneimine-Sepharose coupled at low charge density and heparin/Mn2+ precipitation of the dilute column eluate; 3) DEAE-Sepharose chromatography to remove heparin; 4) chromatography on heparin-Sepharose; 5) preparative sodium dodecyl sulfate-polyacrylamide slab gel electrophoresis followed by electroelution of the protein and ion pair extraction to remove sodium dodecyl sulfate; and 6) high performance liquid chromatography on a reverse-phase C-4 column. Following initial chromatography on polyethyleneimine Sepharose, it was found that substantial subsequent purification could only be achieved under denaturing conditions.     

Sequence and characterization of the Ig heavy chain constant and partial variable region of the mouse strain 129S1

Although the entire mouse genome has been sequenced, there remain challenges concerning the elucidation of particular complex and polymorphic genomic loci. In the murine Igh locus, different haplotypes exist in different inbred mouse strains. For example, the Igh(b) haplotype sequence of the Mouse Genome Project strain C57BL/6 differs considerably from the Igh(a) haplotype of BALB/c, which has been widely used in the analyses of Ab responses. We have sequenced and annotated the 3' half of the Igh(a) locus of 129S1/SvImJ, covering the C(H) region and approximately half of the V(H) region. This sequence comprises 128 V(H) genes, of which 49 are judged to be functional. The comparison of the Igh(a) sequence with the homologous Igh(b) region from C57BL/6 revealed two major expansions in the germline repertoire of Igh(a). In addition, we found smaller haplotype-specific differences like the duplication of five V(H) genes in the Igh(a) locus. We generated a V(H) allele table by comparing the individual V(H) genes of both haplotypes. Surprisingly, the number and position of D(H) genes in the 129S1 strain differs not only from the sequence of C57BL/6 but also from the map published for BALB/c. Taken together, the contiguous genomic sequence of the 3' part of the Igh(a) locus allows a detailed view of the recent evolution of this highly dynamic locus in the mouse.     

Differences in anxiety-related behavior between apolipoprotein E-deficient C57BL/6 and wild type C57BL/6 mice

The influence of ApoE gene deletion on the anxiety state has not been previously investigated. The elevated plus maze was used in this study to determine differences in anxiety-related behavior between apoE-deficient and wild type C57BL/6 mice. The apoE-deficient mice demonstrated less anxiety on the elevated plus maze by spending more time in the open arms of the elevated plus maze compared to wild type mice (p<0.001). Additionally, female apoE-deficient mice visited the open arm of the maze more often than their apoE-deficient male counterpart (p<0.05). The anxiety state and/or sex are possible variables to be considered when designing physiological and/or behavioral studies involving mice that are apoE-deficient.     

Characterization of the serotoninergic system in the C57BL/6 mouse skin.

We showed expression of the tryptophan hydroxylase gene and of tryptophan hydroxylase protein immunoreactivity in mouse skin and skin cells. Extracts from skin and melanocyte samples acetylated serotonin to N-acetylserotonin and tryptamine to N-acetyltryptamine. A different enzyme from arylalkylamine N-acetyltransferase mediated this reaction, as this gene was defective in the C57BL6 mouse, coding predominantly for a protein without enzymatic activity. Serotonin (but not tryptamine) acetylation varied according to hair cycle phase and anatomic location. Serotonin was also metabolized to 5-hydroxytryptophol and 5-hydroxyindole acetic acid, probably through stepwise transformation catalyzed by monoamine oxidase, aldehyde dehydrogenase and aldehyde reductase. Activity of the melatonin-forming enzyme hydroxyindole-O-methyltransferase was notably below detectable levels in all samples of mouse corporal skin, although it was detectable at low levels in the ears and in Cloudman melanoma (derived from the DBA/2 J mouse strain). In conclusion, mouse skin has the molecular and biochemical apparatus necessary to produce and metabolize serotonin and N-acetylserotonin, and its activity is determined by topography, physiological status of the skin, cell type and mouse strain.     

A model of sleep-disordered breathing in the C57BL/6J mouse.

To investigate the pathophysiological sequelae of sleep-disordered breathing (SDB), we have developed a mouse model in which hypoxia was induced during periods of sleep and was removed in response to arousal or wakefulness. An on-line sleep-wake detection system, based on the frequency and amplitude of electroencephalograph and electromyograph recordings, served to trigger intermittent hypoxia during periods of sleep. In adult male C57BL/6J mice (n = 5), the sleep-wake detection system accurately assessed wakefulness (97.2 +/- 1.1%), non-rapid eye movement (NREM) sleep (96.0 +/- 0.9%) and rapid eye movement (REM) sleep (85.6 +/- 5.0%). After 5 consecutive days of SDB, 554 +/- 29 (SE) hypoxic events were recorded over a 24-h period at a rate of 63.6 +/- 2.6 events/h of sleep and with a duration of 28.2 +/- 0.7 s. The mean nadir of fraction of inspired O(2) (FI(O(2))) on day 5 was 13.2 +/- 0.1%, and 137.1 +/- 13.2 of the events had a nadir FI(O(2)) <10% O(2). Arterial blood gases confirmed that hypoxia of this magnitude lead to a significant degree of hypoxemia. Furthermore, 5 days of SDB were associated with decreases in both NREM and REM sleep during the light phase compared with the 24-h postintervention period. We conclude that our murine model of SDB mimics the rate and magnitude of sleep-induced hypoxia, sleep fragmentation, and reduction in total sleep time found in patients with moderate to severe SDB in the clinical setting.