Using whole-genome sequences of the LG/J and SM/J inbred mouse strains to prioritize quantitative trait genes and nucleotides
BMC Genomics volume 16, Article number: 415 (2015)
The laboratory mouse is the most commonly used model for studying variation in complex traits relevant to human disease. Here we present the whole-genome sequences of two inbred strains, LG/J and SM/J, which are frequently used to study variation in complex traits as diverse as aging, bone-growth, adiposity, maternal behavior, and methamphetamine sensitivity.
We identified small nucleotide variants (SNVs) and structural variants (SVs) in the LG/J and SM/J strains relative to the reference genome and discovered novel variants in these two strains by comparing their sequences to other mouse genomes. We find that 39% of the LG/J and SM/J genomes are identical-by-descent (IBD). We characterized amino-acid changing mutations using three algorithms: LRT, PolyPhen-2 and SIFT. We also identified polymorphisms between LG/J and SM/J that fall in regulatory regions and highly informative transcription factor binding sites (TFBS). We intersected these functional predictions with quantitative trait loci (QTL) mapped in advanced intercrosses of these two strains. We find that QTL are both over-represented in non-IBD regions and highly enriched for variants predicted to have a functional impact. Variants in QTL associated with metabolic (231 QTL identified in an F16 generation) and developmental (41 QTL identified in an F34 generation) traits were interrogated and we highlight candidate quantitative trait genes (QTG) and nucleotides (QTN) in a QTL on chr13 associated with variation in basal glucose levels and in a QTL on chr6 associated with variation in tibia length.
We show how integrating genomic sequence with QTL reduces the QTL search space and helps researchers prioritize candidate genes and nucleotides for experimental follow-up. Additionally, given the LG/J and SM/J phylogenetic context among inbred strains, these data contribute important information to the genomic landscape of the laboratory mouse.
The LG/J (large) and SM/J (small) strains of inbred mice were independently derived from selection experiments for large and small body size at 60 days, respectively . LG/J was created from a pool of albino mice obtained from a commercial breeder over a nine-month period and selected for increased body size for over 50 generations [2,3].
SM/J was created from a pool of mice derived from four crosses of seven inbred strains: dilute brown agouti (dba), silver chocolate (sv ba), black and tan (a t), pink-eyed, short-eared dilute brown agouti (ps e dba), albino (c), cinnamon spotted (bs), and agouti (a) . It is unclear whether any of these strains are related to current laboratory inbred strains. The SM/J strain is kept heterozygous at the agouti locus by mating black animals (a/a) with their white-bellied agouti (A w /a) siblings. Attempts to fix the SM/J strain or its Recombinant Inbred Strain offspring for the SM/J allele at agouti have resulted in strain failure .
Subsequent to selection, the LG/J and SM/J strains are fully inbred and have been maintained at the Jackson Laboratory since the 1950s by brother-sister mating. They are at the extremes of the adult body size distribution among the common laboratory inbred strains and have been profitably studied for genetic variation in adult body size and growth [6-18]. Early genetic studies of these two strains found the differences in body size are caused by many genes of individually small effects . This result has been confirmed in many quantitative trait locus (QTL) mapping studies. The strains have remained phenotypically stable over time, except that the SM/J strain is now about 6g heavier than it was in Chai’s studies .
LG/J and SM/J differ in many complex traits in addition to size and growth. The parental strains and their crosses differ in skeletal morphology, including the size and shape of the skull [20-22], mandible [23-30], tooth morphology [31,32], long-bone lengths [33-37], and a variety of other skeletal elements as well as bone biomechanical and structural properties [38,39].
They also differ for a variety of metabolic traits including obesity, diabetes, and serum cholesterol, triglycerides, and free fatty acids levels [40-51]. The SM/J strain responds more strongly than LG/J to a high-fat diet for these metabolic traits. The two strains also differ for maternal genetic effects on offspring growth and offspring adult metabolic traits [52-56] and their cross has been useful in mapping parent-of-origin genetic effects on metabolic traits [48-51,53,57].
In addition to these diverse metabolic and skeletal phenotypes, the LG/J strain shows the rare ability to regenerate tissues after injury. The LG/J strain regenerates ear pinna tissues after a 2mm hole-punch [58-61] while the SM/J strain does not. The LG/J strain also can regenerate damaged articular cartilage  and is protected from post-traumatic osteoarthritis .
Others have mapped a variety of behavioral phenotypes using a LG x SM cross, including prepulse inhibition  and methamphetamine sensitivity and locomotor activity . Additional studies have investigated blood cell parameters  and skeletal muscle weight and fiber types [67,68].
Here we describe the whole-genome sequences of the LG/J and SM/J inbred mouse strains, adding two more sequences to the collection of whole-genomes of laboratory mice [69,70]. We integrate these sequences with quantitative trait loci (QTL), and illustrate how top-down (phenotype to genomic sequence) and bottom-up (genomic sequence to phenotype) approaches can be used to identify candidate quantitative trait genes (QTG) and prioritize positional candidate quantitative trait nucleotides (QTN) for further mechanistic studies. Results from QTL studies using crosses of LG/J and SM/J will be more interpretable given these whole-genome sequences.
Results and discussion
LG/J and SM/J whole-genome sequences
Sequence was generated from DNA isolated from 1 LG/J female and 1 SM/J female as described below. Variants (comprised of small nucleotide variants (SNVs) including single nucleotide polymorphisms (SNPs) and small insertion-deletions (Indels) and structural variants (SVs) including deletions, insertions and inversions) were discovered based on whether they were the same as, or different from, the C57BL/6J reference sequence (mm10, NCBI build 38). Greater than 90% of reads for each strain could be uniquely mapped to the reference genome and, based on a genome size of 2.5G, our coverage is approximately 35X and 30X for LG/J and SM/J, respectively. Overviews of annotated variants identified for each strain individually, as well as those that are polymorphic between LG/J and SM/J, are described in Tables 1 and 2 and in Additional files 1 and 2. SNVs identified for the LG/J strain are available in Additional files 3 and 4 and for the SM/J strain in Additional files 5 and 6. Structural variant positions and classifications for each strain are available in Additional file 7. Polymorphic variants and their genomic context are illustrated in Figure 1. Quality was assessed by comparing our SNP calls with SNPs in dbSNP that were called from unpublished low-coverage sequence generated independently via two separate library preparations for LG/J (≈18X) and SM/J (≈14X) (Table 3).
We compared SNPs identified in LG/J and SM/J to those identified in other mouse strains sequenced at high coverage and identified novel variants (Table 1). Pairwise distances were calculated for LG/J and SM/J from each of these strains and a phylogenetic tree was constructed using all 20 mouse sequences (Additional files 8 and 9). The LG/J and SM/J ancestries are expected given what we know about their origins and previous studies of inbred mouse phylogenies . It has been suggested that the SM/J strain may be related to the current DBA strains, DBA/1J and DBA/2J, as one of the seven strains contributing to the base population that was selected for smaller body weight was a dba strain named after its dilute brown agouti coat color . Comparison of SNPs between LG/J or SM/J and the sequenced DBA/2J shows no specific relationship between SM/J and DBA/2J. However, the percent shared SNPs between inbred mouse strains can be biased by the use of C57BL/6J as the reference genome since variants may be missed due to alignment quality favoring bases that match the reference.
Deleterious amino acid mutations in LG/J and SM/J were predicted using three independent methods: LRT , PolyPhen-2  and SIFT  as implemented in VEP . The LRT compares the probability that a codon has evolved under a conserved model to the probability that a codon has evolved under a neutral model. The conserved model allows a codon to have evolved under negative selection and the neutral model assumes the rates of nonsynonymous and synonymous mutations are not significantly different. PolyPhen-2 uses both sequence- and structure-based predictive features to characterize the functional importance of a mutation, making use of non-redundant protein databases. SIFT is based on the principles of protein evolution and has been applied to a variety of organisms ranging from bacteria to humans. The algorithm is based on sequence homology and uses a median conservation score to measure protein conservation. The numbers of deleterious amino acid predictions for each strain are listed in Table 4 and predicted scores are provided in Additional file 10. The intersections of these predictions for SNPs that are polymorphic between LG/J and SM/J are illustrated in Figure 2. It is worth noting that the three different methods give largely non-overlapping results, which is consistent with a similar analysis using human genomes . This result implies that our ability to reliably annotate whole genome sequences is still somewhat lacking.
In addition to SIFT predictions, VEP identifies variants falling in noncoding and potentially regulatory regions that may impact a gene’s expression. An overview of potentially functional regulatory variants is listed in Table 4 and positions and scores are provided in Additional files 11 and 12.
LG/J, SM/J identical by descent regions
We identified genomic regions that are highly conserved between LG/J and SM/J using a hidden Markov model (Figure 1 and Additional file 13). These regions are defined as stretches of sequenced DNA ≥ 50kbp in length containing little to no polymorphism between the strains. These regions are likely to be identical by descent (IBD), and are non-randomly distributed throughout the genome, being more clustered than expected by chance (p < 2.2e−16, Wald-Wolfowitz test). We classified 39% of the LG/J and SM/J sequenced genomes as IBD.
LG x SM QTL and whole-genome sequences
We integrated all of these data – IBD regions, polymorphisms and functional predictions – with 272 published QTL for metabolic traits (obesity, diabetes and serum-lipids) and for bone-growth traits (femur, humerus, radius and tibia length) mapped in very advanced intercrosses of LGxSM [37,48-51]. LOD scores range from 2.4 to 11.89 and are significant for their respective study. Additional file 14 describes these QTL, including their genomic coordinates, and their integration with various sequence characterizations. Figure 1 illustrates these QTL in relation to the sequence parameters we have generated. We find that 27% of all QTL interrogated in this study are covered by IBD bases, whereas 39% would be expected if they were randomly distributed. We find that both trait-specific QTL peaks (the point of the QTL with the lowest probability of a chance association with a specific trait) and QTL-specific peaks (accounting for pleiotropy, where more than one trait maps to a locus) are over-represented in non-IBD regions (χ 2 = 42.6, df = 1, p = 6.6x10−11 for trait-specific peaks and χ 2 = 16.8, df = 1, p = 4.03x10−05 for QTL-specific peaks). Further, we find that QTL regions are more likely to harbor variants that are predicted to have a functional impact by comparing the number of amino-acid changing variants predicted to be damaging by at least one algorithm and the number of variants predicted to fall in high-impact positions of transcription factor binding sites to empirical distributions of numbers generated from 1000 sets of randomly chosen non-QTL, non-IBD regions of similar size (p < 2.2e−16 for both amino acid changing and TFBS variants, Additional file 15).
Integrating these data provides an opportunity to examine variation between LG/J and SM/J from both bottom-up (genomic sequence to phenotype) and top-down (phenotype to genomic sequence) perspectives with the goal of identifying plausible quantitative trait nucleotides (QTN) for testing mechanistic hypotheses. To illustrate a bottom-up approach, we focus on a SNP identified in the SM/J strain (chr13:104,041,257 A→C) falling in a cMyb TFBS (JASPAR ID MA01001) with a predicted high-impact (Figure 3). This position overlaps regulatory elements including H3K36 and H3K4 histone marks and a DNase1 hypersensitive site (ENSMUSR00000276453). The SNP falls in an intronic region of the oligopeptidase neurolysin, Nln (NM_029447). Nln knockout mice have been shown to be more insulin sensitive and glucose tolerant and to have increased gluconeogenesis relative to littermate controls . In mammals, the liver is the main site of gluconeogenesis. A microarray analysis of hepatic gene expression shows Nln to be highly significantly differentially expressed between LG/J and SM/J, with the SM/J strain’s expression barely detectable . This gene falls within the support intervals of a QTL associated with basal serum glucose levels in an F16 generation of an advanced intercross (AI) between the LG/J and SM/J strains . Thus Nln is an attractive candidate quantitative trait gene (QTG) and this SNP is an attractive quantitative trait nucleotide (QTN) for further mechanistic studies of Nln involvement in glucose metabolism.
To illustrate a top-down approach, we focus on a QTL mapped in an F34 generation of the LGxSM AI that is associated with variation in tibia length at chr6: 20,650,821-23,746,386 (Figure 4). There are 12 protein coding genes and 10 RNA genes falling within the QTL support interval . Two of these genes, Wnt16 (NM_053116) and Ptprz1 (NM_001081306), are involved in limb development and/or bone formation [79,80]. Variants in a third gene, Cped1 (NM_001081351), have recently been associated with variation in bone mineral density in a GWAS . 23% of this QTL is covered by bases falling in IBD regions, but these three genes are located in non-IBD portions of the QTL. Three nonsynonymous SNPs fall in Cped1 and are predicted to have functional consequences by at least one of the algorithms used: LRT, PolyPhen or SIFT. A fourth nonsynonymous SNP is predicted to have functional consequences by all three algorithms making it a highly attractive QTN. The SNP falls in exon 14 of Ptprz1 (ENSMUSE00000619494, position 23016230 C→A; P1676H). This exon is highly conserved. Sanger sequencing confirmed the SM/J variant and the amino acid changing polymorphism between the strains. Ptprz1 encodes the protein R-PTP-Z, which is thought to modulate osteoblast metabolism through dephosphorylating Src, which plays a key role in osteoblast activities such as adhesion and differentiation [82,83]. The amino acid change occurs in the SM/J strain, and the only other strain of sequenced mouse carrying this variant is NZO/HlLtJ. NZO/HlLtJ is a common laboratory strain used as a model of metabolic syndrome because of its extreme obesity and hyperglycemic phenotype . It bears no special relationship with SM/J, despite sharing this particular variant and a similar metabolic phenotype on a high-fat diet.
Here we describe the whole-genome sequences for the LG/J and SM/J inbred mouse strains. LG/J and SM/J are frequently compared in QTL mapping studies because of their great phenotypic diversity, and because this diversity is normally distributed in intercrosses of these two strains. This makes LG x SM an ideal model system for studying the genetic architecture of normal variation in complex traits, because this most closely mimics that found in human populations, where most variation is the result of many interacting genes of individually small effects. Comparison with previous imputation methods using SNPs called from low-coverage sequence for LG/J and SM/J indicated that there were many unreported variants in these strains [51,85]. Our sequence provides higher resolution, which allows us to capture more of the genomic sequence variation found in these strains with greater certainty.
Because we are interested in identifying candidate variants for functional follow-up, our first order of business was to identify regions of the LG/J and SM/J genomes that do not vary between the strains. Identification of IBD genomic regions is a powerful way to narrow existing QTL support intervals, as the regions that do not vary between the strains are unlikely to cause variation in the mapped trait. Further, identification of IBD and non-IBD regions between two strains is a powerful tool for testing mechanistic hypotheses and for planning focused candidate gene studies . We classified 39% of the LG/J and SM/J sequenced genomes as IBD, and we are able to use this classification to narrow QTL support intervals by orders of magnitude when eliminating stretches of DNA that contain little to no variation between the parental strains within QTL support intervals. Eliminating these regions does not explicitly rule them out as harboring potentially causal variants, however focusing on the non-IBD sequence within QTL support intervals – especially on the variants predicted to have functional consequences – allows one to prioritize the so-called ‘low-hanging fruit’. Other commonly used strain pairs, such as C57BL/6J and DBA/2J, have been sequenced but have not been subject to such close evaluation of IBD genomic regions within QTL.
We integrated the rich source of QTL results generated by very advanced intercrosses of LG x SM for metabolic traits (231 QTL associated with variation in obesity, diabetes and serum-lipids) mapped in the F16 and for bone-growth traits (41 QTL associated with variation in femur, humerus, radius and tibia) mapped in the F34 generations with the LG/J and SM/J genomic sequences. The median QTL interval in the F16 data is ≈ 3.5mB and for the F34 it is ≈ 2.0mB (Additional file 14). We find that ≈ 29% of F16 and ≈ 19% of F34 QTL intervals are covered by IBD regions. Subtracting these IBD regions from QTL reduces the genomic search space for both generations, resulting in a comparable ≈ 1.2mB median number of per QTL bases for both F16 and F34 generations. We find that QTL from both generations have proportionally equivalent numbers of SNVs, bases covered by SVs, and bases predicted to have functional consequences (Additional file 14). Further, we find that QTL regions as a class are more likely to harbor variants predicted to have functional consequences relative to randomly selected non-QTL and non-IBD genomic regions of similar size (Additional file 15).
The traits associated with the QTL interrogated here belong to different classes of quantitative phenotype, namely metabolic traits mapped in the F16 generation and developmental traits mapped in the F34 generation. Having QTL for the same trait replicate in multiple generations of a mapping study increases the probability that the QTL is causal. However, given the high costs (both financial and time-wise) of breeding and maintaining large populations of mice for many generations, incorporating whole-genome sequence parameters as we have illustrated here can be used to extract compelling information from QTL mapped in earlier generations of intercrosses, even when the support intervals span many mega-bases of sequence. Thus whole-genome sequence data should be made part of the toolkit used to inform the design, execution and follow-up of QTL mapping studies.
We have highlighted variants falling in a QTL on chromosome 13 associated with basal glucose levels and a QTL on chromosome 6 associated with tibia length (Figures 3 and 4) to illustrate how all of these data – whole genome sequence variants, IBD regions, QTL, and functional predictions – can be integrated with each other and with other public datasets to identify and prioritize QTG and QTN. Identifying such variants within QTL in mice can facilitate translational research for correlated traits in human studies, and potentially uncover genetic underpinnings of disease phenotypes [87,88]. Our initial description and analysis of the LG/J and SM/J whole genomes offers new data that can be used to address fundamental questions about the molecular nature of quantitative phenotypic variation in an important model system.
All animal care and handling procedures conformed to IACUC guidelines.
DNA isolation and library construction
DNA was isolated from the livers of one adult female LG/J and one adult female SM/J mouse using the Qiagen DNAeasy Blood and Tissue kit (Qiagen, West Sussex, UK). Genomic DNA was sonicated to an average size of 175 bp. The fragments were blunt ended, had addition of “A” base to 3’ end, and had Illumina’s sequencing adapters ligated to the ends. The ligated fragments underwent amplification incorporating a unique indexing sequence tag. The resulting libraries were sequenced on 6 lanes using the Illumina HiSeq-2500 as paired end reads extending 101 bases from both ends of the fragments.
Alignment and SNV detection
Reads were aligned to the GRC m38 (mm10) mouse reference genome using NovoAlign-2.08.02 and a BAM file was produced for each strain (http://www.novocraft.com).
The gene annotation model used was Ensembl Mus musculus GRC 38.72. Variants were called from each strain’s BAM file using two variant discovery tools: SAMtools pileup and FreeBayes Bayesian genetic variant detector [89,90]. Results were merged and only variants with a minimum read depth of 3 and a quality score of at least 20 for SNPs and at least 50 for indels were included in the final set of SNVs for analysis.
Quality assessment and Sanger sequencing candidate variant
SNPs identified for LG/J and SM/J were compared to SNPs called from independently generated, low-coverage sequence available on dbSNP  using custom python scripts. DNA was isolated from the livers of 4 female LG/J and 4 female SM/J animals to validate the candidate amino acid changing mutation P1676H in exon 14 of Ptprz1. The following primers were designed: Forward 5’-GCT CCA TGG CCA CTA TCT TTA CTC-3’ and reverse 5’-CAA TTC ATG CCT CAA GGT GAC TGC-3’. Sanger sequencing was performed by Genewiz (South Plainfield, NJ).
Comparison with variants identified in other sequenced mouse strains
VCF files for SNPs discovered in 18 mouse strains’ whole-genomes were downloaded (http://www.sanger.ac.uk/resources/mouse/genomes) and compared to LG/J and SM/J SNPs using custom python scripts.
Structural variant prediction
Structural variants (SV) were identified using a local implementation of the SVmerge pipeline . SVmerge simplifies SV discovery and integration using multiple SV discovery algorithms, producing a single consensus SV callset. The SVmerge pipeline was run with 4 calling algorithms: Breakdancer 1.0, pindel 0.2.4q, SECluster and cnD [92-95]. After merging the results of these 4 callers into a single consensus callset, the Velvet genome assembler was used to attempt breakpoint assembly for all SVs in the callset .
Deleterious mutation prediction
Deleterious amino-acid changing mutations were predicted using a likelihood ratio test (LRT) as described in Chun and Fay (2009). The algorithm was modified to compare the genomes of LG/J and SM/J each to the reference C57BL/6 J. P-values were calculated by comparing twice the log-likelihood ratio of the two models to a χ2 distribution with df = 1. Deleterious amino acid changing mutations were predicted using an LRT cutoff of p < 0.001 while controlling for false positives and false negatives as previously described and modified for mouse genomic sequence .
PolyPhen-2 was downloaded and run locally to predict deleterious mutations in the LG/J and SM/J genomes. The algorithm supports analysis of mouse proteins using prebuilt human models if the option ‘-n mouse’ is specified when preparing the local copy of UniProtKB (www.uniprot.org). Deleterious amino acid changing mutations were predicted for both strains using a recommended false discovery rate (FDR) cutoff of 20% .
VEP and SIFT
VEP was downloaded and run locally to predict deleterious mutations in the LG/J and SM/J genomes using alignments built using the TrEMBL 39.8 protein database. Deleterious amino acid changing mutations were predicted for both strains using a score of < 0.05 and a median conservation cutoff of 3.25. Regulatory annotations are based on the Ensembl regulatory build, which integrates data from ENCODE and several other large scale projects. Transcription factor binding site (TFBS) variant scoring is provided for regulatory regions that have ChIP-seq data to support binding predictions. This is done using the MOODS software , which assigns significance scores by matching polymorphisms against motifs in the JASPAR database . The output was filtered for TFBS sites that are classified as ‘Highly Informative’ by the software.
IBD region identification
To identify regions of shared ancestral background, we clustered segments of observed polymorphisms using a two state hidden Markov model. For each state, we modeled two types of observations: 1) the number of polymorphisms in a 50kbp window; and 2) the observation of a SV in a 50kbp window. The count of polymorphisms is expressed as a Poisson variable, while the occurrence of a SV is a Binomial variable. Parameters for this model were estimated using the EM algorithm implemented in the depmixS4 package in the R programming language . A Wald-Wolfowitz test of randomness was performed using the adehabitatLT package in R .
The QTL interrogated in this study are from previously published studies that mapped to the mm9 (NCBI build 37) mouse reference genome. QTL peaks and support intervals were converted to mm10 (GRC m38) using the Batch Coordinate Conversion (liftOver) tool . Sequence data and results generated here were integrated with other published and publically available data as indicated using custom python and R scripts.
We have made these data available to the community through multiple data portals: LG/J (SAMN03075510) and SM/J (SAMN03075514) raw reads have been submitted to the NCBI Sequence Read Archive and BAM files have been submitted to the Wellcome Trust mouse genomes portal (http://www.sanger.ac.uk/resources/mouse/genomes). Other results from this study are available as Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, and 15.
Small nucleotide variants
Transcription factor binding sites
Quantitative trait loci
- candidate QTG:
Quantitative trait gene
Quantitative trait nucleotide
Single nucleotide polymorphisms
Likelihood ratio test
False discovery rate
Beck JA, Lloyd S, Hafezparast M, Lennon-Pierce M, Eppig JT, Festing MF, et al. Genealogies of mouse inbred strains. Nat Genet. 2000;24(1):23–5.
Goodale HD. A study of the inheritance of body weight in the albino mouse by selection. J Hered. 1938;29:101–12.
Wilson SP, Goodale HD, Kyle WH, Godfrey EF. Long term selection for body weight in mice. The Journal of Heredity. 1971;62:228–34.
MacArthur JW. Genetics of body size and related characters. I.Selecting small and large races of the laboratory mouse. Am Nat. 1944;78(775):142–57.
Hrbek T, de Brito RA, Wang B, Pletscher LS, Cheverud JM. Genetic characterization of a new set of recombinant inbred lines (LGXSM) formed from the inter-cross of SM/J and LG/J inbred mouse strains. Mamm Genome. 2006;17(5):417–29.
Chai CK. Analysis of quantitative inheritance of body size in mice. I: hybridization and maternal influence. Genetics. 1956;41:157–64.
Chai CK. Analysis of quantitative inheritance of body size in mice. III. Dominance. Genetics. 1957;42(5):601–7.
Chai CK. Analysis of quantitative inheritance of body size in mice. II: gene action and segregation. Genetics. 1956;41:167–78.
Chai CK. Analysis of quantitative inheritance of body size in mice. IV. An attempt to isolate polygenes. Genet Res. 1961;2:25–32.
Chai CK. Analysis of quantitative inheritance of body size in mice. V. Effects of small numbers of polygenes on similar genetic backgrounds. Genet Res. 1968;11(3):239–46.
Cheverud JM, Routman EJ, Duarte FA, van Swinderen B, Cothran K, Perel C. Quantitative trait loci for murine growth. Genetics. 1996;142(4):1305–19.
Vaughn TT, Pletscher LS, Peripato A, King-Ellison K, Adams E, Erikson C, et al. Mapping quantitative trait loci for murine growth: a closer look at genetic architecture. Genet Res. 1999;74(3):313–22.
Wu R, Wang Z, Zhao W, Cheverud JM. A mechanistic model for genetic machinery of ontogenetic growth. Genetics. 2004;168(4):2383–94.
Zhao W, Ma C, Cheverud JM, Wu R. A unifying statistical model for QTL mapping of genotype x sex interaction for developmental trajectories. Physiol Genomics. 2004;19(2):218–27.
Lin M, Ma CX, Zhao W, Cheverud JM, Wu R. Mechanistic mapping of ontogenetic growth based on biological principles. Growth Dev Aging. 2005;69(1):31–7.
Zhao W, Chen YQ, Casella G, Cheverud JM, Wu R. A non-stationary model for functional mapping of complex traits. Bioinformatics. 2005;21(10):2469–77.
Long F, Chen YQ, Cheverud JM, Wu R. Genetic mapping of allometric scaling laws. Genet Res. 2006;87(3):207–16.
Hager R, Cheverud JM, Wolf JB. Relative contribution of additive, dominance, and imprinting effects to phenotypic variation in body size and growth between divergent selection lines of mice. Evolution. 2009;63(5):1118–28.
Kramer M, Vaughn TC, Pletscher LS, King-Ellison K, Adams E, Erikson C, et al. Genetic variation in body weight growth and composition in the intercross of large (LG/J) and small (SM/J) inbred strains of mice. Genet Mol Biol. 1998;21:211–8.
Leamy LJ, Routman E, Cheverud J. Quantitative trait loci for early and late developing skull characters in mice: a test of the genetic independence model of morphological integration. Am Nat. 1999;153:201–14.
Wolf JB, Leamy LJ, Routman EJ, Cheverud JM. Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice. Genetics. 2005;171(2):683–94.
Roseman C, Kenney-Hunt J, Cheverud J. Phenotypic integration without modularity: testing hypotheses about the distribution of pleiotropic quantitative trait loci in a continuous space. Evol Biol. 2009;36:282–91.
Cheverud J, Routman E, Irschick D. Pleiotropic effects of individual gene loci on mandibular morphology. Evolution. 1997;51:2004–14.
Cheverud JM, Ehrich TH, Vaughn TT, Koreishi SF, Linsey RB, Pletscher LS. Pleiotropic effects on mandibular morphology II: differential epistasis and genetic variation in morphological integration. J Exp Zool B Mol Dev Evol. 2004;302(5):424–35.
Klingenberg CP, Leamy LJ, Routman EJ, Cheverud JM. Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics. Genetics. 2001;157(2):785–802.
Ehrich TH, Vaughn TT, Koreishi SF, Linsey RB, Pletscher LS, Cheverud JM. Pleiotropic effects on mandibular morphology I. Developmental morphological integration and differential dominance. J Exp Zool B Mol Dev Evol. 2003;296(1):58–79.
Leamy LJ, Pomp D, Eisen EJ, Cheverud JM. Quantitative trait loci for directional but not fluctuating asymmetry of mandible characters in mice. Genet Res. 2000;76(1):27–40.
Leamy LJ, Pomp D, Eisen EJ, Cheverud JM. Pleiotropy of quantitative trait loci for organ weights and limb bone lengths in mice. Physiol Genomics. 2002;10(1):21–9.
Leamy LJ, Klingenberg CP, Sherratt E, Wolf JB, Cheverud JM. A search for quantitative trait loci exhibiting imprinting effects on mouse mandible size and shape. Heredity (Edinb). 2008;101(6):518–26.
Willmore KE, Roseman CC, Rogers J, Cheverud JM, Richtsmeier JT. Comparison of mandibular phenotypic and genetic integration between baboon and mouse. Evol Biol. 2009;36(1):19–36.
Workman MS, Leamy LJ, Routman EJ, Cheverud JM. Analysis of quantitative trait locus effects on the size and shape of mandibular molars in mice. Genetics. 2002;160(4):1573–86.
Leamy LJ, Workman MS, Routman EJ, Cheverud JM. An epistatic genetic basis for fluctuating asymmetry of tooth size and shape in mice. Heredity (Edinb). 2005;94(3):316–25.
Kenney-Hunt JP, Vaughn TT, Pletscher LS, Peripato A, Routman E, Cothran K, et al. Quantitative trait loci for body size components in mice. Mamm Genome. 2006;17(6):526–37.
Norgard EA, Roseman CC, Fawcett GL, Pavlicev M, Morgan CD, Pletscher LS, et al. Identification of quantitative trait loci affecting murine long bone length in a two-generation intercross of LG/J and SM/J mice. J Bone Miner Res. 2008;23(6):887–95.
Pavlicev M, Kenney-Hunt JP, Norgard EA, Roseman CC, Wolf JB, Cheverud JM. Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight. Evolution. 2008;62(1):199–213.
Norgard EA, Jarvis JP, Roseman CC, Maxwell TJ, Kenney-Hunt JP, Samocha KE, et al. Replication of long-bone length QTL in the F9–F10 LG, SM advanced intercross. Mamm Genome. 2009;20(4):224–35.
Norgard EA, Lawson HA, Pletscher LS, Wang B, Brooks VR, Wolf JB, et al. Genetic factors and diet affect long-bone length in the F34 LG, SM advanced intercross. Mamm Genome. 2011;22(3–4):178–96.
Reich MS, Jarvis JP, Silva MJ, Cheverud JM. Genetic relationships between obesity and osteoporosis in LGXSM recombinant inbred mice. Genet Res (Camb). 2008;90(5):433–44.
Carson EA, Kenney-Hunt JP, Pavlicev M, Bouckaert KA, Chinn AJ, Silva MJ, et al. Weak genetic relationship between trabecular bone morphology and obesity in mice. Bone. 2012;51(1):46–53.
Cheverud JM, Vaughn TT, Pletscher LS, Peripato AC, Adams ES, Erikson CF, et al. Genetic architecture of adiposity in the cross of LG/J and SM/J inbred mice. Mamm Genome. 2001;12(1):3–12.
Cheverud JM, Ehrich TH, Hrbek T, Kenney JP, Pletscher LS, Semenkovich CF. Quantitative trait loci for obesity- and diabetes-related traits and their dietary responses to high-fat feeding in LGXSM recombinant inbred mouse strains. Diabetes. 2004;53(12):3328–36.
Cheverud JM, Ehrich TH, Kenney JP, Pletscher LS, Semenkovich CF. Genetic evidence for discordance between obesity- and diabetes-related traits in the LGXSM recombinant inbred mouse strains. Diabetes. 2004;53(10):2700–8.
Fawcett GL, Roseman CC, Jarvis JP, Wang B, Wolf JB, Cheverud JM. Genetic architecture of adiposity and organ weight using combined generation QTL analysis. Obesity (Silver Spring). 2008;16(8):1861–8.
Fawcett GL, Jarvis JP, Roseman CC, Wang B, Wolf JB, Cheverud JM. Fine-mapping of obesity-related quantitative trait loci in an F9/10 advanced intercross line. Obesity (Silver Spring). 2010;18(7):1383–92.
Parker CC, Cheng R, Sokoloff G, Lim JE, Skol AD, Abney M, et al. Fine-mapping alleles for body weight in LG/J x SM/J F and F(34) advanced intercross lines. Mamm Genome. 2011;22(9–10):563–71.
Ehrich TH, Kenney JP, Vaughn TT, Pletscher LS, Cheverud JM. Diet, obesity, and hyperglycemia in LG/J and SM/J mice. Obes Res. 2003;11(11):1400–10.
Ehrich TH, Kenney-Hunt JP, Pletscher LS, Cheverud JM. Genetic variation and correlation of dietary response in an advanced intercross mouse line produced from two divergent growth lines. Genet Res. 2005;85(3):211–22.
Cheverud JM, Lawson HA, Fawcett GL, Wang B, Pletscher LS ARF, Maxwell TJ, et al. Diet-dependent genetic and genomic imprinting effects on obesity in mice. Obesity (Silver Spring). 2011;19(1):160–70.
Lawson HA, Zelle KM, Fawcett GL, Wang B, Pletscher LS, Maxwell TJ, et al. Genetic, epigenetic, and gene-by-diet interaction effects underlie variation in serum lipids in a LG/JxSM/J murine model. J Lipid Res. 2010;51(10):2976–84.
Lawson HA, Lee A, Fawcett GL, Wang B, Pletscher LS, Maxwell TJ, et al. The importance of context to the genetic architecture of diabetes-related traits is revealed in a genome-wide scan of a LG/J x SM/J murine model. Mamm Genome. 2011;22(3–4):197–208.
Lawson HA, Cady JE, Partridge C, Wolf JB, Semenkovich CF, Cheverud JM. Genetic effects at pleiotropic loci are context-dependent with consequences for the maintenance of genetic variation in populations. PLoS Genet. 2011;7(9), e1002256.
Jarvis JP, Kenney-Hunt J, Ehrich TH, Pletscher LS, Semenkovich CF, Cheverud JM. Maternal genotype affects adult offspring lipid, obesity, and diabetes phenotypes in LGXSM recombinant inbred strains. J Lipid Res. 2005;46(8):1692–702.
Hager R, Cheverud JM, Wolf JB. Maternal effects as the cause of parent-of-origin effects that mimic genomic imprinting. Genetics. 2008;178(3):1755–62.
Wolf JB, Vaughn TT, Pletscher LS, Cheverud JM. Contribution of maternal effect QTL to genetic architecture of early growth in mice. Heredity (Edinb). 2002;89(4):300–10.
Wolf JB, Leamy LJ, Roseman CC, Cheverud JM. Disentangling prenatal and postnatal maternal genetic effects reveals persistent prenatal effects on offspring growth in mice. Genetics. 2011;189(3):1069–82.
Wolf J, Cheverud JM. Detecting maternal-effect Loci by statistical cross-fostering. Genetics. 2012;191(1):261–77.
Wolf JB, Cheverud JM, Roseman C, Hager R. Genome-wide analysis reveals a complex pattern of genomic imprinting in mice. PLoS Genet. 2008;4(6), e1000091.
Cheverud JM, Lawson HA, Funk R, Zhou J, Blankenhorn EP, Heber-Katz E. Healing quantitative trait loci in a combined cross analysis using related mouse strain crosses. Heredity (Edinb). 2012;108(4):441–6.
Cheverud JM, Lawson HA, Bouckaert K, Kossenkov AV, Showe LC, Cort L, et al. Fine-mapping quantitative trait loci affecting murine external ear tissue regeneration in the LG/J by SM/J advanced intercross line. Heredity (Edinb). 2014;112:508–18.
Blankenhorn EP, Bryan G, Kossenkov AV, Clark LD, Zhang XM, Chang C, et al. Genetic loci that regulate healing and regeneration in LG/J and SM/J mice. Mamm Genome. 2009;20(11–12):720–33.
Bedelbaeva K, Snyder A, Gourevitch D, Clark L, Zhang XM, Leferovich J, et al. Lack of p21 expression links cell cycle control and appendage regeneration in mice. Proc Natl Acad Sci U S A. 2010;107(13):5845–50.
Rai MF, Hashimoto S, Johnson EE, Janiszak KL, Fitzgerald J, Heber-Katz E, et al. Heritability of articular cartilage regeneration and its association with ear-wound healing. Arthritis Rheum. 2012;64:2300–10.
Hashimoto S, Rai MF, Janiszak KL, Cheverud JM, Sandell LJ. Cartilage and bone changes during development of post-traumatic osteoarthritis in selected LGXSM recombinant inbred mice. Osteoarthritis Cartilage. 2012;20(6):562–71.
Samocha KE, Lim JE, Cheng R, Sokoloff G, Palmer AA. Fine mapping of QTL for prepulse inhibition in LG/J and SM/J mice using F(2) and advanced intercross lines. Genes Brain Behav. 2010;9(7):759–67.
Cheng R, Lim JE, Samocha KE, Sokoloff G, Abney M, Skol AD, et al. Genome-wide association studies and the problem of relatedness among advanced intercross lines and other highly recombinant populations. Genetics. 2010;185(3):1033–44.
Bartnikas TB, Parker CC, Cheng R, Campagna DR, Lim JE, Palmer AA, et al. QTLs for murine red blood cell parameters in LG/J and SM/J F(2) and advanced intercross lines. Mamm Genome. 2012;23:356–66.
Carroll AM, Palmer AA, Lionikas A. QTL analysis of type I and type IIA fibers in soleus muscle in a cross between LG/J and SM/J mouse strains. Front Genet. 2011;2:99.
Lionikas A, Cheng R, Lim JE, Palmer AA, Blizard DA. Fine-mapping of muscle weight QTL in LG/J and SM/J intercrosses. Physiol Genomics. 2010;42A(1):33–8.
Yalcin B, Wong K, Agam A, Goodson M, Keane TM, Gan X, et al. Sequence-based characterization of structural variation in the mouse genome. Nature. 2011;477(7364):326–9.
Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, et al. Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011;477(7364):289–94.
Petkov PM, Graber JH, Churchill GA, DiPetrillo K, King BL, Paigen K. Evidence of a large-scale functional organization of mammalian chromosomes. PLoS Genet. 2005;1(3), e33.
Eppig JT, Blake JA, Bult CJ, Kadin JA, Richardson JE, The_Mouse_Genome_Database_Group. The mouse genome database (MGD): comprehensive resource for genetics and genomics of the laboratory mouse. Nucleic Acids Res. 2012;40(1):881–6.
Chun S, Fay JC. Identification of deleterious mutations within three human genomes. Genome Res. 2009;19(9):1553–61.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81.
McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, Cunningham F. Deriving the consequences of genomic variants with the ensembl API and SNP effect predictor. Bioinformatics. 2010;26(16):2069–70.
Cavalcanti DM, Castro LM, Rosa Neto JC, Seelaender M, Neves RX, Oliveira V, et al. Neurolysin knockout mice generation and initial phenotype characterization. J Biol Chem. 2014;289(22):15426–40.
Partridge CG, Fawcett GL, Wang B, Semenkovich CF, Cheverud JM. The effect of dietary fat intake on hepatic gene expression in LG/J AND SM/J mice. BMC Genomics. 2014;15:99.
Schinke T, Gebauer M, Schilling AF, Lamprianou S, Priemel M, Mueldner C, et al. The protein tyrosine phosphatase rptpzeta is expressed in differentiated osteoblasts and affects bone formation in mice. Bone. 2008;42(3):524–34.
Witte F, Dokas J, Neuendorf F, Mundlos S, Stricker S. Comprehensive expression analysis of all Wnt genes and their major secreted antagonists during mouse limb development and cartilage differentiation. Gene Expr Patterns. 2009;9(4):215–23.
Kemp JP, Medina-Gomez C, Estrada K, St Pourcain B, Heppe DH, Warrington NM, et al. Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment. PLoS Genet. 2014;10(6), e1004423.
Toledo SR, Oliveira ID, Okamoto OK, Zago MA, de Seixas Alves MT, Filho RJ, et al. Bone deposition, bone resorption, and osteosarcoma. J Orthop Res. 2010;28(9):1142–8.
Zambuzzi WF, Milani R, Teti A. Expanding the role of Src and protein-tyrosine phosphatases balance in modulating osteoblast metabolism: lessons from mice. Biochimie. 2010;92(4):327–32.
Bray GA, York DA. Genetically transmitted obesity in rodents. Physiol Rev. 1971;51(3):598–646.
Wang JR, de Villena FP, Lawson HA, Cheverud JM, Churchill GA, McMillan L. Imputation of single-nucleotide polymorphisms in inbred mice using local phylogeny. Genetics. 2012;190(2):449–58.
DiPetrillo K, Wang X, Stylianou IM, Paigen B. Bioinformatics toolbox for narrowing rodent quantitative trait loci. Trends Genet. 2005;21(12):683–92.
Kraja AT, Lawson HA, Arnett DK, Borecki IB, Broeckel U, de Las Fuentes L, et al. Obesity-insulin targeted genes in the 3p26-25 region in human studies and LG/J and SM/J mice. Metabolism. 2012;61:1129–41.
Lawson HA. Animal models of metabolic syndrome. In: Conn M, editor. Animal models of human disease. Cambridge: Elsevier; 2013. p. 243–64. vol. 1.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Garrison E, Marth G: Haplotype-based variant detection from short-read sequencing. In.: arXiv preprint arXiv:1207.3907[q-bio.GN]; 2012.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11.
Wong K, Keane TM, Stalker J, Adams DJ. Enhanced structural variant and breakpoint detection using SVMerge by integration of multiple detection methods and local assembly. Genome Biol. 2010;11(12):R128.
Chen K, Wallis JW, McLellan MD, Larson DE, Kalicki JM, Pohl CS, et al. BreakDancer: an algorithm for high-resolution mapping of genomic structural variation. Nat Methods. 2009;6(9):677–81.
Ye K, Schulz MH, Long Q, Apweiler R, Ning Z. Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics. 2009;25(21):2865–71.
Simpson JT, McIntyre RE, Adams DJ, Durbin R. Copy number variant detection in inbred strains from short read sequence data. Bioinformatics. 2010;26(4):565–7.
Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de bruijn graphs. Genome Res. 2008;18(5):821–9.
Korhonen J, Martinmaki P, Pizzi C, Rastas P, Ukkonen E. MOODS: fast search for position weight matrix matches in DNA sequences. Bioinformatics. 2009;25(23):3181–2.
Bryne JC, Valen E, Tang MH, Marstrand T, Winther O, da Piedade I, et al. JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update. Nucleic Acids Res. 2008;36(Database issue):D102–106.
Visser: Depmix: An R-package for fitting mixture models on mixed multivariate data with Markov dependencies. In: R-package manual and introduction into Dependent Mixture models. 2007. [dssm.unipa.it/CRAN/web/packages/depmix/citation.html]
Calenge C. The package adehabitat for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9.
Karolchik D, Barber GP, Casper J, Clawson H, Cline MS, Diekhans M, et al. The UCSC genome browser database: 2014 update. Nucleic Acids Res. 2014;42(Database issue):D764–770.
This work was supported by NIH-NIDDK K01 DK095003 and P30 DK056341 to HAL and in part by NIMH R01MH10180 to DFC and by NIH-NCRR UL1 RR024992 to the genome technology access center at Washington university in St Louis. The content of this paper is solely the responsibility of the authors and does not necessarily represent the official views of the national institutes of health.
The authors declare that they have no competing interests.
IN, DC, SC, JF, and HL generated and analyzed data. JC and HL wrote the manuscript. All authors edited and approved the final manuscript.
LG/J, SM/J and LG/J x SM/J polymorphic SNPs overview.
LG/J, SM/J and LG/J x SM/J polymorphic indels overview.
SNPs in LG/J discovered by integrating output from SamTools pileup and FreeBayes variant discovery tools.
Indels in LG/J discovered by integrating output from SamTools pileup and FreeBayes variant discovery tools.
SNPs in SM/J discovered by integrating output from SamTools pileup and FreeBayes variant discovery tools.
Indels in SM/J discovered by integrating output from SamTools pileup and FreeBayes variant discovery tools.
Structural variants found in LG/J and SM/J.
Distances between the LG/J and SM/J strains and other sequenced strains.
Phylogenetic tree of 20 sequenced mouse strains.
Amino acid changing mutations and corresponding deleterious predicted scores using SIFT/polyphen/LRT for both LG/J and SM/J strains.
VEP output of variants in predicted regulatory features for both LG/J and SM/J.
VEP output of variants in predicted motifs for both LG/J and SM/J.
Predicted IBD regions in LG/J and SM/J sequenced genomes.
QTL metrics and sequence integration.
Distribution of the number of variants predicted to have functional consequences in QTL intervals and in 1000 sets of randomly selected non-QTL regions of similar size.
About this article
Cite this article
Nikolskiy, I., Conrad, D.F., Chun, S. et al. Using whole-genome sequences of the LG/J and SM/J inbred mouse strains to prioritize quantitative trait genes and nucleotides. BMC Genomics 16, 415 (2015). https://doi.org/10.1186/s12864-015-1592-3