|
|
|
|
Published online before print
February 6, 2007, 10.1101/gr.5941007 Genome Res. 17:293-298, 2007 ©2007 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/07 $5.00
Letter mtDNA phylogeny and evolution of laboratory mouse strains1 Instituto de Patologia e Imunologia Molecular da Universidade do Porto (IPATIMUP), 4200-465 Porto, Portugal; 2 Faculdade de Ciências da Universidade do Porto, 4099-002 Porto, Portugal; 3 Medical Faculty, University of Porto, 4200-319 Porto, Portugal; 4 The Jackson Laboratory, Bar Harbor, Maine 04609, USA; 5 Department of Statistics, University of Glasgow, G12 8QQ Glasgow, United Kingdom
Inbred mouse strains have been maintained for more than 100 years, and they are thought to be a mixture of four different mouse subspecies. Although genealogies have been established, female inbred mouse phylogenies remain unexplored. By a phylogenetic analysis of newly generated complete mitochondrial DNA sequence data in 16 strains, we show here that all common inbred strains descend from the same Mus musculus domesticus female wild ancestor, and suggest that they present a different mitochondrial evolutionary process than their wild relatives with a faster accumulation of replacement substitutions. Our data complement forthcoming results on resequencing of a group of priority strains, and they follow recent efforts of the Mouse Phenome Project to collect and make publicly available information on various strains.
Long before the beginning of mouse genetics, humans in Eastern Asia were already breeding mice in order to obtain animals with different coat colors. Modern mouse genetics, however, did not start until the early 20th century with William Castles studies on inheritance. Most of his mice derived from collections of mice fanciers, and they were the ancestors of many inbred strains that are still used today (Rader 2004
Mouse strains are known to have mixtures of various ancestral genomes from different Mus musculus (house mouse) subspecies (for review, see Yoshiki and Moriwaki 2006
More than 450 mouse inbred strains have been established since the firstDBA/2 (dilute, brown, non-agouti)was developed by Castles student Clarence Cook Little in what would become The Jackson Laboratory (Beck et al. 2000
A mouse strain is defined as inbred if the animals have been crossed brother x sister for
Beck et al. (2000)
Mammalian mitochondrial DNA (mtDNA) is a circular double-stranded molecule that encodes 13 genes of the respiratory chain. Defects in these molecules have been associated with a variety of disorders that may affect different tissues in different ways (for review, see Wallace 1999
The analysis of complete mtDNA sequences of inbred mice is also useful for the establishment of mutation/substitution rates. Given that, in principle, these restricted animal populations have a reasonably well documented history, with no inclusion of foreign DNA, variation that occurs can be generated only by mutation. Moreover, study of mice mitochondrial phylogenies may become helpful in clarifying differences that have been reported between mutation rates estimated from pedigrees and substitution rates calculated from phylogenies (Howell et al. 2003
Two main issues were addressed while performing this work: (1) validation of published genealogies and clarification of the matrilineal origin of common inbred strains and (2) evolutionary analysis of inbred strains. It was developed by focusing on 16 selected strains that are part of the Mouse Phenome Project (Bogue 2003
mtDNA diversity We have sequenced the complete mtDNA of 32 mice belonging to 16 inbred strains (two of each strain): 12 common inbred (129S1/SvImJ, A/J, AKR/J, BALB/cByJ, BTBR T+tf/J, C3H/HeJ, C57BL/6J, DBA/2J, FVB/NJ, KK/HlJ, NOD/LtJ, NZW/LacJ) and four wild-derived inbred strains, each representing a different subspecies (CAST/EiJ, MOLF/EiJ, PWD/PhJ, and WSB/EiJ). Out of the 16 strains sequenced we obtained 12 different haplotypes. For all inbred strains, including the wild-derived, sequencing of the second animal of the same strain always confirmed the first sequence (haplotype) obtained.
We combined our data with 13 complete sequences available online, giving a total of 20 different haplotypes consisting of (1) 11 common inbred strains plus (2) two more divergent (NZB and MilP) common inbred haplotypes, (3) two cell lines, and (4) five wild-derived inbred strains haplotypes. The mtDNAs of all 11 common inbred strains were very similar. Only 15 substitutions were detected in the 16,299 mtDNA nucleotides, 14 inside genes and one in a ribosomal RNA (Supplemental Table 1). No substitutions were observed in the control region. Insertions in transfer RNA genes were also observed in sites previously reported as highly polymorphic (Bayona-Bafaluy et al. 2003
For strains C3H and C57BL/6, we obtained different haplotypes from the published ones. C3H/HeJ differs from C3H/He (Kiebish and Seyfried 2005
The four strains derived from wild mice that we sequenced are highly divergent from each other, and three of them also diverge from the reference C57BL/6 sequence: CAST (M. m. castaneus) differs from the reference in 377 positions, PWD (M. m. musculus) and MOLF differ (M. m. molossinus) in 386 positions. The MOLO strain (M. m. molossinus) mtDNA that was already reported (Akimoto et al. 2005
Phylogenetic analyses
We estimated coalescence times for these branches assuming a divergence time between Mus and Rattus of 12 million yr (paleontological data from Jacobs et al. 1989
Based on the same assumption, we obtained an overall substitution rate for the mouse coding mtDNA of 3.7 x 108 substitutions per site per yr. This value is rather similar to the 1.26 x 108 substitutions per site per yr obtained by Mishmar et al. (2003)
We produced a second tree (not shown) with complete mtDNA data from all common inbred strains, cell lines, and wild-derived mice, but excluding the rat. This analysis confirmed that all common inbred mice group with the M. m. domesticus representative. Based on the coalescence time for the main subspecies of Mus calculated in the above analysis, we obtained similar substitution rates to those derived from the previous group of strains (Table 1). Moreover, we could estimate a substitution rate for the control region of 5.6 x 108 per site per yr, which is about three to four times higher than the substitution rate of the first codon positions and eight to 11 times higher than that of the second codon positions. Third codon positions, however, remain the ones with a highest substitution rate, which is twice the value obtained for the control region.
We estimated the divergence time for these common inbred mice and obtained older times than expected. The divergence of the common inbred strains (excluding the most divergent strains NZB and MilP) was dated at 3000 ± 1400 yr ago (Table 2), although it is documented that the first inbred strain was established only
The mtDNA sequences of the common inbred mice show a star-like network (Fig. 2). The domesticus-derived WSB is the furthest away from the center. Except for the SAM strains, all the strains that belong to groups A (Swiss mice), B (Castles mice), and E (C57-related strains) are in the central groups or less than one polymorphism apart. Four inbred strains are two or more polymorphisms apart from the central groups: KK belongs to group C (strains derived from colonies from China and Japan), cell-line LA9, and the SAM strains. Although these belong to group B, they were crossed with mice from an unknown strain in 1968 (Beck et al. 2000
The maximum likelihood analysis of the wild-derived strains did not permit the establishment of the first split among the three major subspecies, and so we present a tree with a tricotomy at the base (Fig. 1): (1) M. m. castaneus, (2) M. m. musculus, and (3) M. m. domesticus. Not surprisingly, this ambiguity had already been unresolved in previous analyses with mitochondrial genes Cyt b and 12S (Lundrigan et al. 2002
Our results from complete mtDNA sequences agree with others based on particular genes (Lundrigan et al. 2002
A single female origin for lab mice strains Genealogies do not suggest any relationship between the different groups (except for C57-related mice and Castles mice, which both descend from Abbie Lathrops stocks). Our results show, however, that all the common inbred strains appear clustered in the network, irrespectively of the group to which they belong. This finding suggests a common origin for all analyzed groups of common inbred strains. Furthermore, the fact that, in the phylogenetic tree, common inbred strains appear clustered with the M. m. domesticus-derived strain, and far from representatives of other subspecies, shows that the common mitochondrial ancestor of the inbred strains used today belonged to the M. m. domesticus wild subspecies.
Taking together our results and the history of inbred mice, it is likely that the same stocks of M. m. domesticus fancy mice were at the origin of all three groups; the maintenance of separate colonies in different regions would have resulted in the few differences that are now observed. The most divergent inbred strains were NZB and MilP as already reported by Bayona-Bafaluy et al. (2003)
Although few polymorphisms were observed between the different common inbred strains, a couple of comments are worth making: (1) A polymorphism has occurred in the split between two substrains of C3H. C3H/He and C3H/HeJ can therefore be distinguished by the nucleotide at position 8889. (2) SAMP strains, although established only 30 yr ago, appear to have accumulated more substitutions than any other strain. This could be related to their senescence-accelerated phenotype, since mtDNA mutations have been proposed to be implicated in the aging process (Chinnery et al. 2002
Different mtDNA evolution in wild and inbred strains
Although data on RFLPs have placed the divergence inside the Mus musculus species 12 million yr ago (Yonekawa et al. 1981 Taken together, these estimates show a faster mtDNA evolution in inbred than in wild mice with a higher accumulation of replacement substitutions. Three hypotheses may account for these discrepancies: (1) the calibration point that we used (12 million yr for the Mus/Rattus divergence) is too distant relative to the young splits that we are trying to date; (2) mtDNA evolution in inbred strains, with such controlled conditions and restricted crosses, occurs faster than in wild environments; or (3) there are differences in generation times of inbred and wild strains.
The issue of the calibration points has been under debate in recent publications, and it has been suggested that substitution rates extrapolated across different time scales result in invalid date estimates (Ho and Larson 2006
Therefore, a higher substitution rate in inbred mice may still stand, despite this rough time estimate. Curiously, mice cell lines included in our analyses have accumulated a similar number of substitutions to that observed for mouse inbred strains, although it had been reported that relaxation of negative selection resulted in a high substitution rate in HeLa cell lines, when compared with phylogenetically derived divergence of mtDNAs (Herrnstadt et al. 2002
A higher substitution rate could also be apparent if generation times of wild mice are longer than in inbred strains. It is possible that mice breeders cross mice before they would mate in the wild. For C57BL/6 mice it is documented (Genetic Background Resource Manual 2006, The Jackson Laboratory, http://jaxmice.jax.org/literature/manuals/mouse_genetics_resource_manual.pdf) that
Final remarks
Samples and sequencing Livers from 32 mice belonging to 16 inbred strains (129S1/SvImJ, A/J, AKR/J, BALB/cByJ, BTBR T+tf/J, C3H/HeJ, C57BL/6J, DBA/2J, FVB/NJ, KK/HlJ, NOD/LtJ, NZW/LacJ, CAST/EiJ, MOLF/EiJ, PWD/PhJ, and WSB/EiJ) were provided as part of a collaboration with the Mouse Phenome Project (The Jackson Laboratory, Bar Harbor, ME). All mice were male and of about the same age. We sequenced DNA from two animals of each strain, belonging to proximate generations. Potential differences that could have been detected between the two mice would result in one of the following conclusions: (1) there was a sequencing error and/or contamination, (2) the strain was misidentified, or (3) the strain presented mtDNA diversity.
We extracted DNA from livers following a standard phenol-chloroform protocol (digested with proteinase K in buffer TE containing 0.5% SDS, purified with phenol-chlorophorm-isoamyl alcohol, and precipitated with ethanol). We amplified 34 overlapping fragments of
Construction of phylogenetic trees
We then performed maximum likelihood (PAML 3.15; Yang 1997
Due to the difficulty in aligning the D-loop of rat and mouse sequences, the first analysis did not include the control region. After the tree was established, and because of the difference in evolutionary scale between the divergence of wild and inbred strains, we redid the analysis including only (1) the rat; (2) the wild-derived inbred strains; (3) two common inbred strains reported to be more divergent than the others, NZB and MilP; and (4) the C57BL/6 strain, as a representative of all other common inbred strains. Assuming from paleontological data (Jacobs et al. 1989 We also constructed a maximum parsimony phylogenetic network including data from the complete mtDNAs from all common inbred strains and the M. m. domesticus representative (WSB). This summarizes all sequence data and illustrates the extreme lack of variation among the different common inbred strains.
Accession numbers
This work was partially supported by Fundação para a Ciência e a Tecnologia through a research grant to A.G. (SFRH/BD/16518/2004) and by "Programa Operacional Ciência, e Inovação 2010" (POCI 2010), VI Programa-Quadro (20022006).
6 Corresponding author.
E-mail aalmeida{at}ipatimup.pt; fax +351-22-5570799. [Supplemental material is available online at www.genome.org. The sequence data from this study have been submitted to GenBank under accession nos. EF108330EF108345.] Article published online before print. Article and publication date are online at http://www.genome.org/cgi/doi/10.1101/gr.5941007
Akimoto, M., Niikura, M., Ichikawa, M., Yonekawa, H., Nakada, K., Honma, Y., and Hayashi, J. 2005. Nuclear DNA but not mtDNA controls tumor phenotypes in mouse cells. Biochem. Biophys. Res. Commun. 327: 10281035.[CrossRef][Medline] Bandelt, H.-J., Forster, P., Sykes, B., and Richards, M. 1995. Mitochondrial portraits of human populations using median networks. Genetics 141: 743753.[Abstract] Battersby, B.J. and Shoubridge, E.A. 2001. Selection of a mtDNA sequence variant in hepatocytes of heteroplasmic mice is not due to differences in respiratory chain function or efficiency of replication. Hum. Mol. Genet. 10: 24692479. Bayona-Bafaluy, M.P., Acín-Pérez, R., Mullikin, J.C., Park, J.S., Moreno-Loshuertos, R., Hu, P., Pérez-Martos, A., Fernández-Silva, P., Bai, Y., and Enríquez, J.A. 2003. Revisiting the mouse mitochondrial DNA sequence. Nucleic Acids Res. 31: 53495355. Beck, J.A., Lloyd, S., Hafezparast, M., Lennon-Pierce, M., Eppig, J.T., Festing, M.F.W., and Fisher, E.M.C. 2000. Genealogies of mouse inbred strains. Nat. Genet. 24: 2325.[CrossRef][Medline] Bishop, C.E., Boursot, P., Baron, B., Bonhomme, F., and Hatat, D. 1985. Most classical Mus musculus domesticus laboratory mouse strains carry a Mus musculus musculus Y chromosome. Nature 315: 7072.[CrossRef][Medline] Bogue, M.. 2003. Mouse Phenome Project: Understanding human biology through mouse genetics and genomics. J. Appl. Physiol. 95: 13351337. Chenna, R., Sugawara, H., Koike, T., Lopez, R., Gibson, T.J., Higgins, D.G., and Thompson, J.D. 2003. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 31: 34973500. Chinnery, P.F., Samuels, D.C., Elson, J., and Turnbull, D.M. 2002. Accumulation of mitochondrial DNA mutations in ageing, cancer, and mitochondrial disease: Is there a common mechanism? Lancet 360: 13231325.[CrossRef][Medline] Chomyn, A. and Attardi, G. 2003. MtDNA mutations in aging and apoptosis. Biochem. Biophys. Res. Commun. 304: 519529.[CrossRef][Medline] Eppig, J.T., Bult, C.J., Kadin, J.A., Richardson, J.E., and Blake, J.A. the members of the Mouse Genome Database Group 2005. The Mouse Genome Database (MGD): From genes to miceA community resource for mouse biology. Nucleic Acids Res. 33: D471D475. Goios, A., Pereira, L., Bogue, M., Macaulay, V., and Amorim, A. 2006. Mitochondrial DNA sequence of 16 inbred strains Mouse Phenome Database, Accession MPD: 202. The Jackson Laboratory, Bar Harbor, ME, USA. Herrnstadt, C., Preston, G., Andrews, R., Chinnery, P., Lightowlers, R.N., Turnbull, D.M., Kubacka, I., and Howell, N. 2002. A high frequency of mtDNA polymorphisms in HeLa cell sublines. Mutat. Res. 501: 1928.[Medline] Ho, S.Y.W. and Larson, G. 2006. Molecular clocks: When times are a-changin. Trends Genet. 22: 7983.[CrossRef][Medline] Ho, S.Y.W., Phillips, M.J., Cooper, A., and Drummond, A.J. 2005. Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol. Biol. Evol. 22: 15611568. Howell, N., Smejkal, C.B., Mackey, D.A., Chinnery, P.F., Turnbull, D.M., and Herrnstadt, C. 2003. The pedigree rate of sequence divergence in the human mitochondrial genome: There is a difference between phylogenetic and pedigree rates. Am. J. Hum. Genet. 72: 659670.[CrossRef][Medline] Jacobs, L.L., Winkler, D.A., and Murry, P.A. 1989. Modern mammal origins: Evolutionary grades in the Early Cretaceous of North America. Proc. Natl. Acad. Sci. 86: 49924995.[CrossRef][Medline] Jaeger, J.-J., Tong, H., and Buffetaut, E. 1986. The age of MusRattus divergence: Paleontological data compared with the molecular clock. Compt. Rend. Acad. Sci. Paris Ser. 2302: 917922. Kiebish, M.A. and Seyfried, T.N. 2005. Absence of pathogenic mitochondrial DNA mutations in mouse brain tumors. BMC Cancer 5: 102109.[CrossRef][Medline] Loveland, B., Wang, C.R., Yonekawa, H., Hermel, E., and Lindahl, K.F. 1990. Maternally transmitted histocompatibility antigen of mice: A hydrophobic peptide of a mitochondrially encoded protein. Cell 60: 971980.[CrossRef][Medline] Lundrigan, B.L., Jansa, S.A., and Tucker, P.K. 2002. Phylogenetic relationships in the genus Mus, based on paternally, maternally, and biparentally inherited characters. Syst. Biol. 51: 410431.[CrossRef][Medline] McKenzie, M., Trounce, I.A., Cassar, C.A., and Pinkert, C.A. 2004. Production of homoplasmic xenomitochondrial mice. Proc. Natl. Acad. Sci. 101: 16851690. Michaux, J.R., Chevret, P., Filippucci, M.-G., and Macholand, M. 2002. Phylogeny of the genus Apodemus with a special emphasis on the subgenus Sylvaemus using the nuclear IRBP gene and two mitochondrial markers: Cytochrome b and 12S rRNA. Mol. Phylogenet. Evol. 23: 123136.[CrossRef][Medline] Mishmar, D., Ruiz-Pesini, E., Golik, P., Macaulay, V., Clark, A.G., Hosseini, S., Brandon, M., Easley, K., Chen, E., and Brown, M.D., et al. 2003. Natural selection shaped regional mtDNA variation in humans. Proc. Natl. Acad. Sci. 100: 171176. Mizutani, J., Chiba, T., Tanaka, M., Higuchi, K., and Mori, M. 2001. Unique mutations in mitochondrial DNA of senescence-accelerated mouse (SAM) strains. J. Hered. 92: 352355. Moreno-Loshuertos, R., Acín-Pérez, R., Fernández-Silva, P., Movilla, N., Pérez-Martos, A., de Cordoba, S.R., Gallardo, M.E., and Enríquez, J.A. 2006. Differences in reactive oxygen species production explain the phenotypes associated with common mouse mitochondrial DNA variants. Nat. Genet. 38: 12611268.[CrossRef][Medline] Pearson, H.. 2004. Mouse sequencing plan aims to boost models. Nature 432: 5.[Medline] Rader, K.. 2004. Making mice: Standardizing animals for American biomedical research, 19001955. Princeton University Press, Princeton, NJ. 1st ed. Sakai, T., Kikkawa, Y., Miura, I., Inoue, T., Moriwaki, K., Shiroishi, T., Satta, Y., Takahata, N., and Yonekawa, H. 2005. Origins of mouse inbred strains deduced from whole-genome scanning by polymorphic microsatellite loci. Mamm. Genome 16: 1119.[CrossRef][Medline] Sligh, J.E., Levy, S.E., Waymire, K.G., Allard, P., Dillehay, D.L., Nusinowitz, S., Heckenlively, J.R., MacGregor, G.R., and Wallace, D.C. 2000. Maternal germ-line transmission of mutant mtDNAs from embryonic stem cell-derived chimeric mice. Proc. Natl. Acad. Sci. 97: 1446114466. Suzuki, H., Shimadaa, T., Terashimaa, M., Tsuchiyab, K., and Aplin, K. 2004. Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences. Mol. Phylogenet. Evol. 33: 626646.[CrossRef][Medline] Takeda, T., Hosokawa, M., and Higuchi, K. 1997. Senescence-accelerated mouse (SAM): A novel murine model of senescence. Exp. Gerontol. 32: 105109.[CrossRef][Medline] Wade, C.M., Kulbokas, E.J., Kirby, A.W., Zody, M.C., Mullikin, J.C., Lander, E.S., Lindblad-Toh, K., and Daly, M.J. 2002. The mosaic structure of variation in the laboratory mouse genome. Nature 420: 574578.[CrossRef][Medline] Wallace, D.C.. 1999. Mitochondrial diseases in man and mouse. Science 283: 14821488. Yang, Z.. 1997. PAML: A program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13: 555556. Yonekawa, H., Moriwaki, K., Gotoh, O., Hayashi, J.-I., Watanabe, J., Miyashita, N., Petras, M.L., and Tagashira, Y. 1981. Evolutionary relationships among five subspecies of Mus musculus based on restriction enzyme cleavage patterns of mitochondrial DNA. Genetics 98: 801816. Yonekawa, H., Moriwaki, K., Gotoh, O., Miyashita, N., Migita, S., Bonhomme, F., Hjorth, J.P., Petras, M.L., and Tagashira, Y. 1982. Origins of laboratory mice deduced from restriction patterns of mitochondrial DNA. Differentiation 22: 222226.[CrossRef][Medline] Yoshiki, A. and Moriwaki, K. 2006. Mouse phenome research: Implications of genetic background. ILAR J. 47: 94102.[Medline]
Received September 12, 2006; accepted in revised format December 8, 2006. This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||