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.     

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.     

Adrenal function and the effect of a high-fat diet on C57BL/6J and C57BL/6J-ob/ob mice.

In C57BL/6J mice and the ob/+ and ob/ob mutants total plasma corticosterone levels were found to be statistically different. In C57BL/6J mice the level was 1.9 +/- 0.2 mug/100 ml plasma, in ob/+ mice 8.6 +/- 1.6 mug/100 ml and in ob/ob mice 13.7 +/- 1.5 mug/100 ml. The percentage of protein-bound corticosterone as well as the free endogenous corticosterone levels were also different. Feeding a high-fat diet to young C57BL/6J and C57BL/6J-ob/ob mice for a period of 4 weeks had no effect upon blood glucose, plasma insulin and plasma corticosterone levels. The significantly higher increase in body weight of the high-fat diet groups of both lines of mice was mainly due to fat cell hypertrophy.     

Hippocampal lipids linked to spatial memory in the C57BL/6J mouse

Although the role of individual brain lipids for learning and memory has been reported, no systematic approach associating brain lipids with spatial memory has been carried out. It was therefore the aim of the study to determine brain lipids in hippocampus of mice forming and yoked controls that did not form spatial memory using the probe trial as the endpoint. 10 animals were trained in the Morris water maze (MWM) and 10 mice were serving as yoked controls i.e. no platform was used during the whole experiment. Hippocampal tissue lipids were extracted and data were acquired with Fourier transformation ion cyclotron resonance mass spectrometry (LTQ-FT) coupled to HPLC. Glycerophosphatidylethanolamines (18:0/22:6, 18:0/20:4 and 18:1/18:1), plasmalogens (16:0-10/22:6 and 18:0-10/22:6) and ceramides (18:0) showed higher levels in the trained group, while glycerolysophosphatidylcholines (16:0, 18:1, 18:0, 20:4), sphingomyelins (16:0, 24:1), ether linked glycerophosphatidylcholines (16:0-10/18:0), glycerophosphatidylcholines (16:0/18:1, 16:0/18:0, 18:0/18:1, 38:7, 18:1/20:1, 20:4/20:4, 22:1/18:1, 22:0/18:1, 20:4/22:6, 22:6/22:6), glucosylceramide (24:1) and plasmalogen (18:0-10/20:1) revealed lower levels in the trained group. Decreased levels of certain species of lysophosphatidylcholine, sphingomyelin, plasmenylphosphatidylcholine, phosphatidylcholine, glycosylceramide and plasmalogen at the probe trial for spatial memory may indicate catabolism in terms of consumption during this process. Increased hippocampal levels of long chain highly unsaturated phosphatidylethanolamines, plasmalogens and ceramides may reflect increased synthesis or decreased degradation at the endpoint of memory testing, probably representing interactions in the brain lipid pathways. The study shows pathways involved in spatial memory, may propose the use of individual brain lipids as probable cognitive enhancers and forms the basis for further studies on the role of brain lipids per se.     

Genetic variation in C57BL/6 ES cell lines and genetic instability in the Bruce4 C57BL/6 ES cell line

Genetically modified mouse strains derived from embryonic stem (ES) cells are powerful tools for gene function analysis. ES cells from the C57BL/6 mouse strain are not widely used to generate mouse models despite the advantage of a defined genetic background. We assessed genetic variation in six such ES cell lines with 275 SSLP markers. Compared to C57BL/6, Bruce4 differed at 34 SSLP markers and had significant heterozygosity on three chromosomes. BL/6#3 and Dale1 ES cell lines differed at only 3 SSLP makers. The C2 and WB6d ES cell lines differed at 6 SSLP markers. It is important to compare the efficiency of producing mouse models with available C57BL/6 ES cells relative to standard 129 mouse strain ES cells. We assessed genetic stability (the tendency of cells to become aneuploid) in 110 gene-targeted ES cell clones from the most widely used C57BL/6 ES cell line, Bruce4, and 710 targeted 129 ES cell clones. Bruce4 clones were more likely to be aneuploid and unsuitable for ES cell-mouse chimera production. Despite their tendency to aneuploidy and consequent inefficiency, use of Bruce4 ES cells can be valuable for models requiring behavioral studies and other mouse models that benefit from a defined C57BL/6 background. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.