Skip to main content

Sequence divergence of Mus spretus and Mus musculus across a skin cancer susceptibility locus



Mus spretus diverged from Mus musculus over one million years ago. These mice are genetically and phenotypically divergent. Despite the value of utilizing M. musculus and M. spretus for quantitative trait locus (QTL) mapping, relatively little genomic information on M. spretus exists, and most of the available sequence and polymorphic data is for one strain of M. spretus, Spret/Ei. In previous work, we mapped fifteen loci for skin cancer susceptibility using four different M. spretus by M. musculus F1 backcrosses. One locus, skin tumor susceptibility 5 (Skts5) on chromosome 12, shows strong linkage in one cross.


To identify potential candidate genes for Skts5, we sequenced 65 named and unnamed genes and coding elements mapping to the peak linkage area in outbred spretus, Spret/EiJ, FVB/NJ, and NIH/Ola. We identified polymorphisms in 62 of 65 genes including 122 amino acid substitutions. To look for polymorphisms consistent with the linkage data, we sequenced exons with amino acid polymorphisms in two additional M. spretus strains and one additional M. musculus strain generating 40.1 kb of sequence data. Eight candidate variants were identified that fit with the linkage data. To determine the degree of variation across M. spretus, we conducted phylogenetic analyses. The relatedness of the M. spretus strains at this locus is consistent with the proximity of region of ascertainment of the ancestral mice.


Our analyses suggest that, if Skts5 on chromosome 12 is representative of other regions in the genome, then published genomic data for Spret/EiJ are likely to be of high utility for genomic studies in other M. spretus strains.


Mus spretus mice are derived from wild mice collected in Southern France, Spain and Northern Africa [1, 2]. M. spretus diverged from Mus musculus one to three million years ago[3, 4]. Recent data suggest that these mice are reported to have an average of one sequence variant every 50 bp [5]. Strains of M. spretus are frequently utilized in combination with M. musculus strains in quantitative trait loci (QTL) studies due to the high degree of sequence and phenotypic diversity between the strains[6, 7]. Indeed, M. spretus mice have been valuable for the identification of loci contributing to differences in immune response and inflammation [810]. M. spretus strains are also cancer resistant in comparison to several M. musculus strains and have been utilized for mapping of cancer susceptibility loci [1116].

Unlike some of the other commonly used strains of mice for QTL mapping, none of the M. spretus strains have been fully sequenced. Few reports of comparisons of sequence diversity within M. spretus strains exist, and most have been generally limited to one gene or small genomic region [1719]. The majority of the sequence, microsatellite, and single nucleotide polymorphism (SNP) data available for M. spretus mice are from Spret/Ei; for researchers using non-Spret/Ei strains of M. spretus, such as SEG/PAS mice, there is limited strain-specific information available in public databases.

In previous studies, we observed a large difference between M. musculus and M. spretus strains in susceptibility to chemically-induced skin cancer. M. spretus mice chemically treated with 7, 12-dimethylbenz [a]anthracene (DMBA) and 12-O-tetradecanoylphorbol-13-acetate (TPA) develop no skin papillomas after 20 weeks of treatment but M. musculus strains, such as NIH/Ola and FVB/N, develop an average of 20 or more papillomas [11]., To identify loci for skin cancer susceptibility we created four F1 backcrosses between M. spretus (Spret/EiJ, STF/PAS, outbred spretus) and M. musculus (FVB/NT, NIH/Ola) mice. Fifteen regions of linkage were identified [1113, 20, 21]. One locus, skin tumor susceptibility 5 (Skts5), a region on chromosome 12 with peak linkage at 17 cM showed strong evidence of linkage in one F1 backcross NIH/Ola × F1(NIH/Ola × outbred spretus).

As one means of identifying potential skin cancer susceptibility candidate genes at Skts5, we sequenced all predicted genes in four of the strains used for linkage NIH/Ola, FVB/NJ, outbred spretus, and Spret/EiJ. Exons that contained amino acid alterations between the strains were sequenced in two additional strains of M. Spretus (SEG/PAS and STF/PAS) as well as FVB/NT (FVB/N mice, Taconics). Here, we describe the results from over 140 kb of sequence data from Skts5, 40.1 kb of which was sequenced in seven strains, and a phylogenetic comparison of the relatedness of four M. spretus strains.


Confirmation of Skts5 as a skin tumor susceptibility locus

Skin tumor susceptibility 5 (Skts5) was originally mapped as a skin tumor susceptibility locus to mouse chromosome 12 in a genome wide linkage analysis for skin cancer in an F1 backcross between NIH/Ola and outbred spretus mice [11, 12]. We performed additional crosses between skin tumor resistant mice (Spret/EiJ, STF/PAS) and skin tumor susceptible mice (NIH/Ola, FVB/NT) to further test this locus for linkage [13, 20, 21]. We observed no evidence for linkage at Skts5 in the interspecific backcross NIH/Ola × [STF/PAS × NIH/Ola] or in the Spret/EiJ by FVB/NT F1 interspecific backcross. The linkage analysis for the Spret/EiJ by NIH/Ola F1 interspecific backcross (approximately 80 mice) was inconclusive for linkage at Skts5. To refine the linkage region in the NIH/Ola × outbred spretus backcross, we genotyped additional markers on chromosome 12. Linkage analyses showed a peak LOD score of 3.2 at approximately 17 cM (35 Mb; Figure 1).

Figure 1

Linkage data at Skts5. Eighteen markers mapping to chromosome 12 were genotyped in 351 NIH/Ola × (NIH/Ola × spretus outbred) backcross mice. Analysis using WinQTLcart shows two peak linkages. Skts15 maps to proximal chromosome 12 and Skts5 maps to chromosome 12 with a peak linkage at ABI marker E_29.924 (35 Mb, 17 cM). The one negative lod confidence interval is marked by a dotted line.

Outbred spretus mice are not bred to homogeneity and segregate different polymorphisms and haplotypes across the genome. However, there was no evidence from the original linkage data generated from nine different outbred spretus fathers using 18 polymorphic markers for heterogeneity at Skts5. To further rule out the possibility of heterogeneity within these mice at Skts5, we genotyped seven additional microsatellite markers mapping to the peak linkage region of Skts5 in DNA from six non-littermate outbred spretus male mice. We observed no evidence of genetic heterogeneity (data not shown). From these results we felt reasonably confident that there was no genetic heterogeneity within our outbred spretus mouse colony at this locus.

Identification of variants between M. musculus and M. spretus

In order to identify polymorphic coding variants that fit the linkage data for Skts5, we first identified all known coding and non-coding genes at the peak region for linkage on chromosome 12 (between 31.7 and 47.9 Mb, Ensembl build 50) using genetic maps from the UCSC Genome Browser, Ensembl and Entrez Gene databases. Maps were updated and throughout our project to include newly mapped genes. At the conclusion of this study, 65 coding elements including ten hypothetical genes/coding elements, nine snoRNAs and one microRNA were identified as mapping within Skts5 (see Additional file 1). We designed PCR primers flanking all exons for the genes or coding elements mapping to Skts5. When multiple splice forms of a gene were predicted, we sequenced exons representing all splice variants.

To carry out a screen for coding differences between the resistant strain of mice (outbred spretus) and the susceptible strain (NIH/Ola) that could be candidate variants for Skts5, we initially sequenced DNA from NIH/Ola, FVB/NJ, Spret/EiJ, and one outbred spretus animal. Exons containing amino acid variants were then sequenced in FVB/NT, SEG/PAS and STF/PAS. Variants shared between outbred spretus and STF/PAS were eliminated as candidates since STF/PAS does not show evidence for linkage at Skts5 in the NIH/Ola backcross. PCR product was obtained from both M. spretus and M. musculus DNA for all exons which rules out the possibility of large deletions in one strain of mice at this locus. Our analysis did not rule out the possibility of causal variants in regulatory regions or duplication or inversion copy number polymorphisms at this locus.

Analysis of sequence alignments of forward and reverse reads of 140,448 bp from genes and coding elements mapping to Skts5 led to the identification of 1,123 polymorphisms between M. musculus and M. spretus strains (Tables 1, 2 and Additional file 2). Of these, 122 resulted in amino acid substitutions, 121 occurred between M. musculus and M. spretus strains and one occurred between NIH/Ola and FVB/NJ. There were six sequence differences identified between FVB/NJ and NIH/Ola, but no differences between FVB/NT and NIH/Ola.

Table 1 Amino acid variants within M. spretus strains
Table 2 Potentially functionally significant amino acid variants

Characterization of nucleotide substitutions at Skts5

To further characterize the nucleotide variants, we analyzed polymorphisms on a gene by gene basis. The number of polymorphisms identified per gene ranged from no changes to 65. The rate of substitution per nucleotide ranged from 3/3593 (0.08%) in Nova1 to 8/133 (6%) in SNORNA17 (Ensembl ENSMUSG00000077294). Because synonomous changes are less likely to affect protein function, these were not evaluated further in the other strains of mice used in the linkage crosses [2224]. A catalogue of the types and locations of polymorphisms (Table 3: see Additional file 2) at Skts5 include the following changes: nonsense (1), missense (122), small insertion/deletions (14), 5'UTR (89), 3'UTR (404) and potential splice-site alterations (23). To determine if the type of substitutions identified at Skts5 could shed light on the evolutionary forces between M. spretus and M. musculus, we looked at the number and percentages of nucleotide changes and the type of mutations (transitions versus transversions) in this data set (see Additional file 2). The transition to transversion ratio was approximately 2.0 (Figure 2). Given that C to T (G to A) transitions are a common mutation type due to the CpG methylation deamination process; these data are consistent with other large sequence comparisons[25, 26].

Table 3 Characteristics of nucleotide substitutions at Skts5
Figure 2

Type and number of nucleotide substitutions. The number of each type of single nucleotide polymorphism identified at Skts5 between NIH/Ola and Spret/EiJ is listed. The transition to transversion rate is roughly 2-fold.

Identification of variants within M. spretus

We rationalized that in order for a bi-allelic polymorphism to be considered a candidate for Skts5 that the cancer resistant strain (outbred spretus) would have one allele and the susceptible strain (NIH/Ola) would have the other. In addition, closely related strain STF/PAS that does not show evidence for linkage at Skts5 should not share alleles with outbred spretus. Of the 121 amino acid substitutions identified between M. musculus and M. spretus, three variants were present in Spret/EiJ but not the other M. spretus strains (Table 1). To rule out a sequencing or PCR artifact for those variants, new PCR product was generated and sequenced. All variants were confirmed as being unique to Spret/EiJ. Of the other 118 amino acid variants which were polymorphic between NIH/Ola and M. spretus, 99 were shared in all M. spretus strains analyzed and 19 were polymorphic between different M. spretus strains (Table 1). Eight variants fit our criteria for candidacy based on the linkage data (Table 4).

Table 4 Candidate variants for Skts5

Conservation of amino acid alterations

Just as non-synonomous changes have a greater likelihood of affecting protein function than synonomous changes, non-synonomous changes that occur in a evolutionarily conserved residue are more likely to be functionally significant [27, 28]. To assess amino acid polymorphisms for potential functional significance, proteins were compared using Ensembl ortholog prediction alignments between M. musculus and Homo sapiens, Pan troglyodytes, Gallus gallus, Rattus norvegicus, Canis familiaris, Bos Taurus, Xenopus tropicalis, Takifugu rubripes, and Monodelphis domestica when protein homologies existed If only one amino acid in the orthologous protein position was present in at least eight of the nine non-musuclus species examined, we considered that position to be strongly conserved. Twenty-four amino acids showed strong conservation by this criterion (see Additional file 3).

In addition to sequence alignments, we assessed amino acid polymorphisms for potential functionality through a sequence homology-based tool Sorting Intolerant From Tolerant (SIFT)[28]. Using SIFT, the majority of amino acid polymorphisms were scored as tolerant substitutions. Surprisingly, a few amino acids that appear to be highly conserved across species, such as Pik3cg Ala405Thr, were scored as well-tolerated substitutions using SIFT. Eight variants were identified that were scored as intolerant substitutions by SIFT (scores 0.05 or less). These variants include Ahr Phe348Leu, Ankmy2 Ile369Thr, Sostdc1 Cys15Ser, Ifrd1 Leu449Phe, EG629820 Gly47Ser, Gly70Ser and Asn72Tyr and F730043M19Rik Asn133Ile. Of these, we identified three (Ahr Phe348Leu, Sostdc1 Cys15Ser and Ifrd1 Leu449Phe) in our initial ortholog comparison studies (see Additional file 3). Seven additional variants in other genes had SIFT scores of less than 0.15 and may have mild functional effects (Table 2).

Of the 22 amino acid alterations that were polymorphic within M. spretus strains, four showed amino acid substitutions that were considered to be intolerant using SIFT (Table 5). These changes all occurred in genes for which there is limited homology data for other species which may lead to over prediction of intolerance by SIFT. Five other substitutions resulted in SIFT scores of 0.15 or lower including 1110049B09Rik Pro691Ala and Lrrn3 Ala69Gln. Polymorphisms in highly conserved amino acids may result in functional differences in these proteins and could account for phenotypic differences between M. spretus strains.

Table 5 Evolutionary conservation of variants seen between M. Spretus strains

Polymorphisms at splice sites

Variants that fall at highly conserved splice acceptor or splice donor sites may also affect protein function due to alterations in splicing. Twenty-four variants mapped within 3 bp of a splice site. To test whether these variants affected splicing, we isolated tail RNA from NIH/Ola, outbred spretus and Spret/EiJ and performed RT-PCR using primers from flanking exons. Identically sized bands of roughly the same intensity were observed in all of the strains for nineteen of the primer sets tested indicating that these variants are unlikely to affect splicing (data not shown). Five variants were not tested for effects on splice site alterations because the variations mapped to the first or last exon of the gene.

Comparison of sequence divergence within M. Spretus strains

Because this was, to our knowledge, the most sequence data generated from multiple M. spretus strains, we were interested in using this information to determine the degree of relatedness between M. spretus strains. We therefore conducted maximum parsimony analyses on exon sequence for Spret/EiJ, SEG/PAS, STF/PAS and outbred spretus to FVB/NJ, and NIH/Ola utilizing Phylip 3.66 [30]. We did not include FVB/NT in the analysis as we did not identify any sequence differences between FVB/NT and NIH/Ola. We included 40,122 bp from coding and flanking exonic regions across Skts5 on mouse chromosome 12 which represented sequence from 90 exons from 35 genes with the identified amino acid alterations and the non-coding polymorphisms within them. DNA variants were all weighted equally. A maximum parsimony tree was generated from the data which was identical to a consensus tree from 100 bootstrap reiterations (Figure 3). As expected, the greatest phylogenetic distance occurred between the M. musculus and M. spretus groups. The fewest number of polymorphisms occurred between the M. musculus strains FVB/NJ and NIH/Ola (FVB/NT) mice. These data suggest that the inbred Swiss-derived strains are more closely related than the wild-derived M. spretus strains. Within M. spretus strains, Spret/EiJ and outbred spretus were the most closely related, and STF/PAS mice showed the most divergence from the other M. spretus strains. Nonetheless, of the roughly 13,300 amino acids assessed in all four M. spretus strains only 0.17% (22 amino acids) were polymorphic between M. spretus strains. If Skts5 is representative of the genome as a whole, then most single nucleotide polymorphism data available for Spret/EiJ from coding regions will have utility for other M. spretus strains.

Figure 3

Phylogenetic tree from sequence data at Skts5. Phylogenetic tree generated by parsimony analysis of four M. Spretus and two Swiss-derived M. musculus strains. DNAPars program from Phylip 3.66 was used to generate a parsimony tree from 40.1 kb of sequence data from mouse chromosome 12 derived from six of the seven mouse strains. Because no sequence differences existed between FVB/NT and NIH/Ola these mice were grouped on the tree. The four M. spretus strains cluster together. The STF/PAS mice show the greatest divergence within M. spretus.

Discussion and Conclusion

To our knowledge, this is the largest sequence-based comparison of multiple strains of M. spretus. Based on sequence comparison of over 40 kb across a 15-Mb region on chromosome 12, M. spretus strains are more closely related to each other than to the Swiss-inbred strains of M. musculus. Outbred spretus and Spret/EiJ have the highest sequence similarity. This is not unexpected as both of these M. spretus mice strains were derived from wild mice isolated from the same region of Spain. STF/PAS mice are derived from mice isolated in Tunisia. This strain showed the highest divergence from the other M. spretus strains and was the M. spretus strain most closely related to the Swiss-derived M. musculus mice. These data fit with the geographical ascertainment of the lines of M. spretus as the other M. spretus strains all come from Spain.

In this dataset we observed an average of one polymorphism every 130 bp between M. musculus and M. spretus. This is less than half as frequent as previous reports of sequence variance between M. musculus and M. spretus [5]. However, this data set was almost entirely restricted to coding regions with only minimal amounts of flanking intronic sequence. Genes are under different mutation constraints than non-functional intragenic regions and therefore show fewer variants [31]. Within the 65 genes analyzed in this study, 31 had no amino acid polymorphisms, 19 genes had one to three amino acid variants and 15 genes had greater than four amino acid variants. Despite the number of non-synonomous changes, the majority of amino acid variants identified in this study were not predicted to show functional consequences using SIFT. Genes and amino acids that showed lower SIFT tolerance scores are being considered as candidates for skin cancer susceptibility.

We conducted this analysis to identify candidate genes for Skts5. From our criteria based on the linkage, a number polymorphisms were identified that were consistent with the linkage (Table 4, Additional file 2). These include eight amino acid changes that were found in outbred spretus that were not seen in STF/PAS including three that have SIFT scores below 0.10. These variants will be considered as candidates for the linkage. However, the nucleotide driving linkage at Skts5 may also be in a non-sequenced regulatory or intragenic region and may affect gene expression or stability. There is also the small possibility that undetected heterogeneity at Skts5 exists between the outbred spretus mice and we have not sequenced the genome containing the resistance allele. Future studies will address these possibilities.

The five genes that showed the highest number of amino acid polymorphisms (Cog5, Slc26a3, Ahr, Zfp277, 1110049B09Rik) do not share any known function and do not map next to each other at Skts5. Ahr (aryl hydrocarbon receptor) mediates the toxicity of dioxin compounds, particularly polycylic aromatic hydrocarbons [32]. It has a large number of reported amino acid variants in both mouse and humans. Some variants within Ahr have been associated with differences in chemical toxicity in mouse and man [3335] and in modest effects on cancer risk in humans [36, 37]. One of the amino acids identified in our study shows a very low SIFT score and may therefore affect Ahr function, however this variant was observed in all M. spretus strains and is therefore not considered to be a candidate for Skts5. Two Ahr variants, Asn533Ser and Asn544Asp, fit with the linkage data but are not predicted to be functionally significant by SIFT. Cog5 is part of the conserved oligomeric Golgi (cog) complex and has been postulated to mediate Golgi trafficking. Two amino acid polymorphisms have been reported in human but none in mouse. Little is known about Zfp27. Slc26a3 is a transmembrane glycoprotein that transports sulfate in the lower gut. Slc26a3 has two amino acids reported in human and none in mouse. The function of 1110049B09Rik is not known, but the purported human homologue, FLJ23834, has a large number of reported amino acid alterations as well (seven in mouse; six in humans). It is possible that the higher rate of amino acid substitution in these genes reflects some underlying selection pressures between M. musculus and M. spretus.

A number of potentially functional interesting variants were identified that do not fit with the linkage data. The only species in which the Sostdc1 amino acid variant Cys15Ser was observed was in M. spretus. It is not known if substitution of a cysteine at position 15 with serine in the M. spretus isoform results in altered disulfide bond formation, disruption of other secondary structure or other functional consequence. Another interesting polymorphism is in 4921508M14Rik. This gene contains a stop codon at position 109 in all of the M. spretus strains theoretically resulting in a truncated protein missing the 37 terminal amino acids. The function of this gene is unknown. EG629820, which contains several amino acid variants thought to be functional by SIFT, is a predicted gene that has shown expression differences of unknown biological consequence in a number of murine studies [38].

There was a relatively low degree of sequence variation between the M. spretus strains; however we identified twenty-two changes that resulted in intra-spretus amino acid polymorphisms. Thus, based on this dataset, nearly twenty percent of the amino acid alterations identified between M. spretus and M. musculus are polymorphic within M. spretus. In contrast, sequencing results from the two Swiss-derived M. musculus strains identified only one amino acid that was polymorphic between the FVB/NJ strain and NIH/Ola. Sequence comparison of three closely related M. musculus s trains: NIH/Ola, FVB/NT and FVB/NJ revealed only six of 140,448 basepairs (0.004%) that were polymorphic between NIH/Ola and FVB/NJ. There were no sequence differences identified between FVB/NT and NIH/Ola. Because sequence variants exist between FVB/NT and FVB/NJ, researchers comparing results across studies should be aware of which strain of FVB/N is used. Our results confirm previous studies in which variations in coding regions were more common in wild-mice than in laboratory-derived strains [7], perhaps reflecting differences in breeding and selection pressure.

Our data comparing sequence diversity between M. spretus strains may not be representative of genome as a whole since this study was limited to one locus on chromosome 12 and included only genes, predicted coding elements, and small amounts of flanking intronic sequences. In this study, we conducted our sequence comparisons by analyzing only exons in which we initially identified an amino acid variation between outbred spretus, Spret/EiJ, FVB/N, and NIH/Ola. There may be additional amino acid variations at Skts5 in which Spret/EiJ, outbred spretus and NIH/Ola share an allele that is different from the other strains of M. spretus. We also cannot rule out the possibility that Skts5 may have unique genomic properties compared to other loci. Previous studies have reported different phylogenetic trees when comparing sequence samples from across the genome (global comparison) with gene-specific phylogenetic trees (local comparison)[39]. Skts5 may represent a region with unique local characteristics between M. spretus strains. For instance, Skts5 may be under differential evolutionary constraints against sequence divergence in M. spretus strains. Nonetheless, our data suggest that genomic data and single nucleotide markers in the database for Spret/EiJ, particularly intragenic SNPs, are more likely to be shared with other M. spretus strains than with M. musculus strains, increasing the utility of the publicly available SNP maps for studies utilizing any M. spretus mouse.


Animal material

The ancestral and laboratory origins of each of the mice in the study are as follows: SEG/PAS (Grenada, Spain, inbred line maintained by the Institute Pasteur), STF/PAS (Fondouk el Djedid, Tunisia, inbred line maintained by the Institute Pasteur), Spret/EiJ (Puerto Real, Ibiza, Spain, inbred line maintained by the Jackson Laboratories), outbred spretus (Spain, outbred line maintained by Stephen Brown, PhD of the Medical Research Council, Harwell, England), FVB/NT (Switzerland, inbred line maintained by Taconics, US), FVB/NJ (Switzerland, inbred line maintained by the Jackson laboratories) and NIH/Ola (Switzerland, inbred line maintained by Harlan Olac, UK). DNA was isolated from tails or spleens of outbred spretus, NIH/Ola and FVB/NT mice using standard methods [40]. Spret/EiJ and FVB/NJ DNA was purchased from the Jackson Laboratories. STF/PAS and SEG/PAS DNAs were a gift from Xavier Montagutelli, DVM, PhD of the Institut Pasteur.

Linkage analysis and genotyping of markers at Skts5

To refine linkage at Skts5 in the outbred spretus × NIH/Ola backcross ten additional single nucleotide polymorphisms from chromosome 12 were genotyped in mice at Applied Biosystems. Genotypes of the new markers were combined with previously genotyped microsatellite markers D12Mit182, D12Mit242, D12Mit60, D12Mit153, D12Mit154, D12Mit2, D12Mit68, and D12Mit5. Quantitative trait locus analysis using WinQTLCart [41] was conducted on the genotypes of the 18 markers on 351 intercross mice using skin papilloma number as the trait. To rule out genetic heterogeneity within outbred spretus mice at Skts5 (15–19 cM), we genotyped seven additional markers (D12Mit243, D12Mit187, D12Mit146, D12Mit186, D12Mit222, D12Mit61 and D12Mit62) in six different outbred spretus fathers, Spret/EiJ, NIH/Ola, and NIH/Ola × outbred spretus F1 and STF/PAS.

PCR primer design, sequencing and analysis

Intron/exon boundaries of all genes and coding regions mapping to Skts5 were identified using the Ensembl database, builds 35–48. When multiple splice forms of a gene were present, exons from all splice forms were included in the analyses. We designed PCR primers using Integrated DNA Technology's SciTools PrimerQuest web-based program [42]. Primer sequences and PCR conditions are listed in Additional file 4. PCR products were treated with Exo/SAP-IT to remove single stranded DNA (USB). Automated sequencing of PCR products was conducted on an ABI 3700 by standard methods. Primers used for PCR were also used for the sequencing. Forward and reverse sequences were analyzed and compared using DNAstar 3.0. The sequence traces were inspected visually whenever a nucleotide substitution was indicated. When a polymorphism was identified between FVB/NJ and Spret/EiJ, PCR products from additional strains NIH/Ola, FVB/NT, outbred spretus, SEG/PAS and STF/PAS were sequenced.

Evolutionary comparison

Sequence alignments between M. musculus and other species were obtained from the Ensembl database when homologies were present. Tolerance scores were generated using the SIFT program and entering in the NCBI protein G numbers for each M. musculus gene at Skts5 [29]. SIFT sorts intolerant amino acid substitutions from tolerant substitutions by protein based alignments from other species. A tolerance score of 0 to 1.0 is generated for each amino acid variant. Low tolerance indices (less than 0.05) are suggested of non-tolerated (deleterious or functional) amino acid substitutions. Higher tolerance indices indicate that an amino acid substitution is less likely to have functional consequences.

Splicing Studies

RNA was isolated from tails of NIH/Ola, Spret/EiJ and outbred spretus mice using standard Trizol methods per manufacturer's recommendations (Invitrogen). One microgram of RNA was reversed transcribed using the Iscript cDNA synthesis kit (Biorad) according to manufacturer's directions. PCR primers were designed in exons flanking the potentially aberrantly spliced exon. Primer sequences and PCR conditions are listed in Additional file 4. Genomic DNA from the mice was used as a negative control. PCR products were run on 1.5% agarose gels for comparison of product sizes.

Phylogenetic comparison

Maximum DNA parsimony analyses were conducted using sequence from exons that were sequenced in all strains of mice using the DNAPars program included in the Phylogeny Interference Package (Phylip) version 3.66 [30]. All nucleotide changes were weighted equally. To assess support for the most parsimonious tree, we also performed a bootstrap analysis and compared the resulting consensus tree to the most parsimonious tree. This resulted in an identical tree. We used the program Drawgram in Phylip 3.66 to construct a phylogenetic tree from this data without using branch length considerations.



Phylogeny Interference Package


Quantitative trait loci

Skts5 :

Skin tumor susceptibility 5


Sorting tolerant from intolerant.


  1. 1.

    Bonhomme F, Guenet JL: The laboratory mouse and its wild relatives. Genetic variants and strains of the laboratory mouse. Edited by: Lyon MF, Rastan S, Brown SDM. 1996, Oxford, England: Oxford University Press, 2: 1577-1596. 3

    Google Scholar 

  2. 2.

    Hochepied T, Schoonjans L, Staelens J, Kreemers V, Danloy S, Puimege L, Collen D, Van Roy F, Libert C: Breaking the species barrier: derivation of germline-competent embryonic stem sells from Mus spretus × C57BL/6 Hybrids. Stem Cells. 2004, 22 (4): 441-447. 10.1634/stemcells.22-4-441.

    PubMed  Article  Google Scholar 

  3. 3.

    She JX, Bonhomme F, Boursot P, Thaler L, Catzeflis FM: Molecular phylogenics in the genus Mus: comparative analysis of electrophoretic scnDNA hybridization and mtDNA RFLP data. Biol J Linn Soc. 1990, 41 (1): 83-103. 10.1111/j.1095-8312.1990.tb00823.x.

    Article  Google Scholar 

  4. 4.

    Suzuki H, Shimada T, Terashima M, Tsuchiya K, Aplin K: Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences. Mol Phylogenet Evol. 2004, 33 (3): 626-646. 10.1016/j.ympev.2004.08.003.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    Zhang J, Wheeler DA, Yakub I, Wei S, Sood R, Rowe W, Liu PP, Gibbs RA, Buetow KH: SNPdetector: A Software Tool for Sensitive and Accurate SNP Detection. PLoS Comput Biol. 2005, 1 (5): e53-10.1371/journal.pcbi.0010053.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    de Gouyon B, Melanitou E, Richard MF, Requarth M, Hahn IH, Guenet JL, Demenais F, Julier C, Lathrop GM, Boitard C: Genetic analysis of diabetes and insulitis in an interspecific cross of the nonobese diabetic mouse with Mus spretus. Proc Natl Acad Sci USA. 1993, 90 (5): 1877-1881. 10.1073/pnas.90.5.1877.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  7. 7.

    Guenet JL, Bonhomme F: Wild mice: an ever-increasing contribution to a population mammalian model. Trends Genet. 2003, 19 (1): 24-31. 10.1016/S0168-9525(02)00007-0.

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Staelens J, Wielockx B, Puimege L, Van Roy F, Guenet JL, Libert C: Hyporesponsiveness of SPRET/Ei mice to lethal shock induced by tumor necrosis factor and implications for a TNF-based antitumor therapy. Proc Natl Acad Sci USA. 2002, 99 (14): 9340-9345. 10.1073/pnas.142293699.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  9. 9.

    Staelens J, Puimege L, Mahieu T, Pynaert G, Hochepied T, Vandenabeele A, Grooten J, Kontoyiannis D, Van Roy F, Kollias G: Response of TNF-hyporesponsive SPRET/Ei mice in models of inflammatory disorders. Mamm Genome. 2004, 15 (7): 537-543. 10.1007/s00335-004-3002-z.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Stephan K, Smirnova I, Jacque B, Poltorak A: Genetic analysis of the innate immune responses in wild-derived inbred strains of mice. Eur J Immunol. 2007, 37 (1): 212-223. 10.1002/eji.200636156.

    PubMed  CAS  Article  Google Scholar 

  11. 11.

    Nagase H, Bryson S, Cordell H, Kemp CJ, Fee F, Balmain A: Distinct genetic loci control development of benign and malignant skin tumors in mice. Nat Genet. 1995, 10 (4): 424-429. 10.1038/ng0895-424.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Nagase H, Mao JH, Balmain A: A subset of skin tumor modifier loci determines survival time of tumor-bearing mice. Proc Natl Acad Sci USA. 1999, 96 (26): 15032-15037. 10.1073/pnas.96.26.15032.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  13. 13.

    Ewart-Toland A, Briassouli P, de Koning JP, Mao JH, Yuan J, Chan F, MacCarthy-Morrogh L, Ponder BA, Nagase H, Burn J: Identification of Stk6/STK15 as a candidate modifier of cancer risk in mouse and man. Nat Genet. 2003, 34 (4): 403-412. 10.1038/ng1220.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Benelli R, Peissel B, Manenti G, Gariboldi M, Vanzetto C, Albini A, Dragani TA: Allele-specific patterns of the mouse parathyroid hormone-related protein: influences on cell adhesion and migration. Oncogene. 2003, 22 (49): 7711-7715. 10.1038/sj.onc.1207088.

    PubMed  CAS  Article  Google Scholar 

  15. 15.

    Santos J, Montagutelli X, Acevedo A, Lopez P, Vaquero C, Fernandez M, Szatanik M, Salido E, Guenet JL, Fernandez-Piqueras J: A new locus for resistance to gamma-radiation-induced thymic lymphoma identified using inter-specific consomic and inter-specific recombinant congenic strains of mice. Oncogene. 2005, 21 (43): 6680-6683. 10.1038/sj.onc.1205846.

    Article  Google Scholar 

  16. 16.

    To MD, Perez-Losada J, Mao JH, Hsu J, Jacks T, Balmain A: A functional switch from lung cancer resistance to susceptibility at the Pas1 locus in Kras2LA2 mice. Nat Genet. 2006, 38 (8): 926-930. 10.1038/ng1836.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  17. 17.

    Greene-Till R, Zhao Y, Hardies SC: Gene flow of unique sequences between Mus musculus domesticus and Mus Spretus. Mamm Genome. 2000, 11 (3): 225-230. 10.1007/s003350010041.

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Prager EM, Tichy H, Sage RD: Mitochondrial DNA sequence variation in the eastern house mouse, Mus Musculus: comparison with other house mice and report of a 75-bp tandem repeat. Genetics. 1996, 143 (1): 427-446.

    PubMed  CAS  PubMed Central  Google Scholar 

  19. 19.

    Villa-Morales M, Santos J, Fernandez-Piqueras J: Functional fas (cd95/apo-1) promoter polymorphisms in inbred mouse strains exhibiting different susceptibility to gamma-radiation induced thymic lymphoma. Oncogene. 2006, 25 (14): 2022-2029. 10.1038/sj.onc.1209234.

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Nagase H, Mao JH, de Koning JP, Minami T, Balmain A: Epistatic interactions between skin tumor modifier loci in interspecific (Spretus/musculus) backcross mice. Cancer Res. 2001, 61 (4): 1305-1308.

    PubMed  CAS  Google Scholar 

  21. 21.

    Mao JH, Saunier EF, de Koning JP, McKinnon MM, Higgins MN, Nicklas K, Yang HT, Balmain A, Akhurst RJ: Genetic variants of Tgfb1 act as context-dependent modifiers of mouse skin tumor susceptibility. Proc Natl Acad Sci USA. 2006, 103 (21): 8125-8130. 10.1073/pnas.0602581103.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  22. 22.

    Chasman D, Adams RM: Predicting the functional consequences of non-synonymous single nucleotide polymorphisms: structure-based assessment of amino acid variation. J Mol Bio. 2001, 307 (2): 683-706. 10.1006/jmbi.2001.4510.

    CAS  Article  Google Scholar 

  23. 23.

    Chamary JV, Parmley JL, Hurst LD: Hearing silence: non-neutral evolution at synonymous sites in mammals. Nat Rev Genet. 2006, 7 (2): 98-108. 10.1038/nrg1770.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Sunyaev S, Ramensky V, Koch I, Lathe W, Kondrashov AS, Bork P: Prediction of deleterious human alleles. Hum Mol Genet. 2001, 10 (6): 591-597. 10.1093/hmg/10.6.591.

    PubMed  CAS  Article  Google Scholar 

  25. 25.

    Sjoblom T, Jones S, Wood LD, Parson DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N: The consensus coding sequences of human breast and colorectal cancers. Science. 2006, 314 (5797): 268-274. 10.1126/science.1133427.

    PubMed  Article  Google Scholar 

  26. 26.

    Zhang Z, Gerstein M: Patterns of nucleotide substitution, insertion and deletion in the human genome inferred from pseudogenes. Nucl Acids Res. 2003, 31 (18): 5338-5348. 10.1093/nar/gkg745.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  27. 27.

    Savas S, Taylor IW, Wrana JL, Ozcelik H: Functional non-synonymous single nucleotide polymorphisms from the TGF-{beta} protein interaction network. Phyisiol Genomics. 2007, 29 (2): 109-117. 10.1152/physiolgenomics.00226.2006.

    CAS  Article  Google Scholar 

  28. 28.

    Zhu Y, Spitz MR, Amos CI, Lin J, Schabath MB, Wu X: An Evolutionary Perspective on Single-Nucleotide Polymorphism Screening in Molecular Cancer Epidemiology. Cancer Res. 2004, 64 (6): 2251-2257. 10.1158/0008-5472.CAN-03-2800.

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Ng PC, Henikoff SS: Accounting for human polymorphisms predicted to affect protein function. Genome Res. 2002, 12 (13): 436-46. 10.1101/gr.212802.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  30. 30.

    PHYLLIP. []

  31. 31.

    Liu GE, Matukumalli LK, Sonstegard TS, Shade LL, Van Tassel CP: Genomic divergences among cattle, dog and human estimated from large-scale alignments of genomic sequences. BMC Genomics. 2006, 7: 140-153. 10.1186/1471-2164-7-140.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Shimuzu Y, Nakatsuru Y, Ichinose M, Takahashi Y, Kume H, Mimura J, Fujii-Kuriyama Y, Ishikawa T: Benzo[a]pyrene carcinogenicity is lost in mice lacking the aryl hydrocarbon receptor. Proc Natl Acad Sci USA. 2000, 97 (2): 779-782. 10.1073/pnas.97.2.779.

    Article  Google Scholar 

  33. 33.

    Connor KT, Aylward LL: Human response to dioxin: aryl hydrocarbon receptor (AhR) molecular structure, function, and dose-response data for enzyme induction indicate an impaired human AhR. J Toxicol Environ Health B Crit Rev. 2006, 9 (2): 147-171. 10.1080/15287390500196487.

    PubMed  CAS  Article  Google Scholar 

  34. 34.

    Okey AB, Franc MA, Moffat ID, Tijet N, Boutros PC, Korkalainen M, Tuomisto J, Pohjanvirta R: Toxicological implications of polymorphisms in receptors for xenobiotic chemicals: The case of the aryl hydrocarbon receptor. Toxicol Appl Pharmacol. 2005, 207 (2 Suppl): 43-51. 10.1016/j.taap.2004.12.028.

    PubMed  Article  Google Scholar 

  35. 35.

    Abnet CC, Fagundes RB, Strickland PT, Kamangar F, Roth MJ, Taylor PR, Dawsey SM: The influence of genetic polymorphisms in Ahr, CYP1A1, CYP1A2, CYP1B1, GST M1, GST T1 and UGT1A1 on urine 1-hydroxypyrene glucuronide concentrations in healthy subjects from Rio Grande do Sul, Brazil. Carcinogenesis. 2007, 28 (1): 112-117. 10.1093/carcin/bgl131.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Kim JH, Kim H, Lee KY, Kang JW, Lee KH, Park SY, Yoon HI, Jheon SH, Sung SW, Hong YC: Aryl hydrocarbon receptor gene polymorphisms affect lung cancer risk. Lung Cancer. 2007, 56 (1): 9-15. 10.1016/j.lungcan.2006.11.010.

    PubMed  Article  Google Scholar 

  37. 37.

    Long JR, Egan KM, Dunning L, Shu XO, Cai Q, Cai H, Dai Q, Holtzman J, Gao YT, Zheng W: Population-based case-control study of AhR (aryl hydrocarbon receptor) and CYP1A2 polymorphisms and breast cancer risk. Pharmacogenet Genomics. 2006, 16 (4): 237-243. 10.1097/01.fpc.0000189803.34339.ed.

    PubMed  CAS  Article  Google Scholar 

  38. 38.

    Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, Suzuki H, Carninci P, Hayashizaki Y, Wells C, Frith M, Ravasi T, Pang KC, Hallinan J, Mattick J, Hume DA, Lipovich L, Batalov S, Engström PG, Mizuno Y, Faghihi MA, Sandelin A, Chalk AM, Mottagui-Tabar S, Liang Z, Lenhard B, Wahlestedt C, RIKEN Genome Exploration Research Group; Genome Science Group (Genome Network Project Core Group); FANTOM Consortium: Antisense transcription in the mammalian transcriptome. Science. 2005, 309 (5740): 1564-1566. 10.1126/science.1112009.

    PubMed  Article  Google Scholar 

  39. 39.

    Park YG, Lukes L, Yang H, Debies MT, Samant RS, Welch DR, Lee M, Hunter KW: Comparative sequence analysis in eight inbred strains of the metastasis modifier QTL candidate gene Brms1. Mamm Genome. 2002, 13 (6): 289-292. 10.1007/s00335-001-2151-6.

    PubMed  CAS  Article  Google Scholar 

  40. 40.

    Laird PW, Zijderveld A, Linders K, Rudnicki MA, Jaenisch R, Burns A: Simplified mammalian DNA isolation procedure. Nucl Acids Res. 1991, 19 (15): 4293-10.1093/nar/19.15.4293.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  41. 41.

    Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer 2.5. 2007, Department of Statistics, North Carolina State University, Raleigh, NC, []

    Google Scholar 

  42. 42.

    IDT SciTools Primer Quest. []

Download references


We thank Dr. Laura Kubatko for her thoughtful advice on phylogenetic analyses. We thank the Ohio State University Comprehensive Cancer Center Nucleic Acids Shared Resource for sequencing support (supported by a NIH/NCI cancer center support grant). This work was supported in part by the Concern Foundation and the American Cancer Society RSG-07-083-01-MGO (AET).

Author information



Corresponding author

Correspondence to Amanda Ewart Toland.

Additional information

Authors' contributions

KLM participated in sequence analysis, phylogenetic analysis and organization and interpretation of data. JLF, HYC, NG and AMD contributed to sequence analysis and organization of data. AB and JHM contributed to linkage data. AET conceived the study, participated in data analysis and coordination and helped to draft the manuscript. All authors participating in editing of the final manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table 1: Genes mapping to Skts5. This table shows a list of all genes mapping to Skts5 including their Ensembl reference number. (XLS 18 KB)


Additional file 2: Table 2: All sequence variants identified at Skts5. This table shows all the sequence variants identified in the study, the strain of mouse in which it was found, the amino acid change if applicable, and the rs number if previously entered in dbSNP. (XLS 220 KB)

Additional file 3: Table 3: Genes showing strong conservation across species. This table shows a list of amino acid changes that appear to be strongly conserved across species. In red are the amino acids that match M. musculus and in blue are the amino acids that match M. Spretus. The corresponding amino acids in nine species as well as the mouse strains sequenced in the study are shown. (XLS 23 KB)

Additional file 4: Table 4: Primers for each exon for sequencing and splice site analysis. This table lists all forward and reverse primers used for sequencing in this study. (XLS 128 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Mahler, K.L., Fleming, J.L., Dworkin, A.M. et al. Sequence divergence of Mus spretus and Mus musculus across a skin cancer susceptibility locus. BMC Genomics 9, 626 (2008).

Download citation


  • Quantitative Trait Locus
  • Linkage Data
  • Amino Acid Variant
  • Peak Linkage
  • Sort Intolerant From Tolerant