Higher differentiation among subspecies of the house mouse (Mus musculus) in genomic regions with low recombination

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species are enhanced in these regions. Here, we study the patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with collinear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation‐with‐gene‐flow model to estimate divergence times and effective population sizes (Ne) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time approximately 350 000 years ago. Overall, Ne for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest Ne and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest Ne and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either were from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. FST was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.

[1]  M. Nachman,et al.  Recombination and Speciation: Loci Near Centromeres Are More Differentiated Than Loci Near Telomeres Between Subspecies of the European Rabbit (Oryctolagus cuniculus) , 2009, Genetics.

[2]  G. A. Watterson.  On the number of segregating sites in genetical models without recombination. , 1975, Theoretical population biology.

[3]  F. Bonhomme,et al.  Molecular phylogenies in the genus Mus: Comparative analysis of electrophoretic, scnDNA hybridization, and mtDNA RFLP data , 1990 .

[4]  Keith A. Gardner,et al.  Hybrid zones and the genetic architecture of a barrier to gene flow between two sunflower species. , 1999, Genetics.

[5]  M. Nei,et al.  Mathematical model for studying genetic variation in terms of restriction endonucleases. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[6]  R. Nielsen,et al.  Distinguishing migration from isolation: a Markov chain Monte Carlo approach. , 2001, Genetics.

[7]  M. Nachman,et al.  Asymmetry and Polymorphism of Hybrid Male Sterility During the Early Stages of Speciation in House Mice , 2007, Evolution; international journal of organic evolution.

[8]  N. Saitou,et al.  Genetic variation versus recombination rate in a structured population of mice. , 2003, Molecular biology and evolution.

[9]  M. Nachman,et al.  Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula , 2008, Molecular ecology.

[10]  M. Noor,et al.  Chromosomal inversions and the reproductive isolation of species , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Annie Orth,et al.  Hybridation naturelle entre deux sous-espèces de souris domestique, Mus musculus domesticus et Mus musculus castaneus, près du lac Casitas (Californie) , 1998 .

[12]  B. Charlesworth,et al.  The effect of deleterious mutations on neutral molecular variation. , 1993, Genetics.

[13]  Peter Donnelly,et al.  A comparison of bayesian methods for haplotype reconstruction from population genotype data. , 2003, American journal of human genetics.

[14]  Joel Dudley,et al.  MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences , 2008, Briefings Bioinform..

[15]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[16]  M. Nachman,et al.  THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE , 2010, Evolution; international journal of organic evolution.

[17]  B. Payseur,et al.  Extensive recombination rate variation in the house mouse species complex inferred from genetic linkage maps. , 2011, Genome research.

[18]  M. Nei.  Analysis of gene diversity in subdivided populations. , 1973, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Hyuna Yang,et al.  On the subspecific origin of the laboratory mouse , 2007, Nature Genetics.

[20]  R. Faria,et al.  Chromosomal speciation revisited: rearranging theory with pieces of evidence. , 2010, Trends in ecology & evolution.

[21]  August E. Woerner,et al.  Recombination-filtered genomic datasets by information maximization , 2007, Bioinform..

[22]  R. Nielsen,et al.  Multilocus Methods for Estimating Population Sizes, Migration Rates and Divergence Time, With Applications to the Divergence of Drosophila pseudoobscura and D. persimilis , 2004, Genetics.

[23]  S WRIGHT,et al.  Genetical Structure of Populations , 1950, British medical journal.

[24]  J. Mattick.  Genome research , 1990, Nature.

[25]  C. A. Machado,et al.  Inferring the history of speciation from multilocus DNA sequence data: the case of Drosophila pseudoobscura and close relatives. , 2002, Molecular biology and evolution.

[26]  F. Bonhomme,et al.  Inferences of selection and migration in the Danish house mouse hybrid zone , 2005 .

[27]  M. Nachman,et al.  Inferring the history of speciation in house mice from autosomal, X‐linked, Y‐linked and mitochondrial genes , 2008, Molecular ecology.

[28]  Martin S. Taylor,et al.  A High-Resolution Single Nucleotide Polymorphism Genetic Map of the Mouse Genome , 2006, PLoS biology.

[29]  Jody Hey,et al.  Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics , 2007, Proceedings of the National Academy of Sciences.

[30]  Eleazar Eskin,et al.  A sequence-based variation map of 8.27 million SNPs in inbred mouse strains , 2007, Nature.

[31]  M. Slatkin,et al.  Estimation of levels of gene flow from DNA sequence data. , 1992, Genetics.

[32]  Gary A. Churchill,et al.  A New Standard Genetic Map for the Laboratory Mouse , 2009, Genetics.

[33]  Colin N. Dewey,et al.  Fine-Scale Phylogenetic Discordance across the House Mouse Genome , 2009, PLoS genetics.

[34]  P. Tucker,et al.  ABRUPT CLINE FOR SEX CHROMOSOMES IN A HYBRID ZONE BETWEEN TWO SPECIES OF MICE , 1992, Evolution; international journal of organic evolution.

[35]  M. Nachman,et al.  Contrasting Patterns of Introgression at X-Linked Loci Across the Hybrid Zone Between Subspecies of the European Rabbit (Oryctolagus cuniculus) , 2006, Genetics.

[36]  J. Forejt,et al.  Genetic analysis of X-linked hybrid sterility in the house mouse , 2004, Mammalian Genome.

[37]  M. Noor,et al.  Islands of speciation or mirages in the desert? Examining the role of restricted recombination in maintaining species , 2009, Heredity.

[38]  J. Forejt.  Hybrid sterility in the mouse , 1996 .

[39]  M. Whitlock,et al.  Indirect measures of gene flow and migration: FST≠1/(4Nm+1) , 1999, Heredity.

[40]  J. M. Smith,et al.  The hitch-hiking effect of a favourable gene. , 1974, Genetical research.

[41]  D. Nickerson,et al.  PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. , 1997, Nucleic acids research.

[42]  N. Barton,et al.  ACCUMULATING POSTZYGOTIC ISOLATION GENES IN PARAPATRY: A NEW TWIST ON CHROMOSOMAL SPECIATION , 2003, Evolution; international journal of organic evolution.

[43]  Richard R. Hudson,et al.  Generating samples under a Wright-Fisher neutral model of genetic variation , 2002, Bioinform..

[44]  D. Reich,et al.  Human Population Differentiation Is Strongly Correlated with Local Recombination Rate , 2010, PLoS genetics.

[45]  C. Aquadro,et al.  African and North American populations of Drosophila melanogaster are very different at the DNA level , 1993, Nature.

[46]  M. Przeworski,et al.  A new approach to estimate parameters of speciation models with application to apes. , 2007, Genome research.

[47]  M. Nei.  Molecular Evolutionary Genetics , 1987 .

[48]  L H. Rieseberg,et al.  Chromosomal rearrangements and speciation. , 2001, Trends in ecology & evolution.

[49]  T. A. Hall,et al.  BIOEDIT: A USER-FRIENDLY BIOLOGICAL SEQUENCE ALIGNMENT EDITOR AND ANALYSIS PROGRAM FOR WINDOWS 95/98/ NT , 1999 .

[50]  K. Aplin,et al.  Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences. , 2004, Molecular phylogenetics and evolution.

[51]  Cestmir Vlcek,et al.  A Mouse Speciation Gene Encodes a Meiotic Histone H3 Methyltransferase , 2009, Science.

[52]  R. Harrison,et al.  Genealogical Discordance and Patterns of Introgression and Selection Across a Cricket Hybrid Zone , 2009, Evolution; international journal of organic evolution.

[53]  F. Tajima.  Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. , 1989, Genetics.

[54]  J. Wakeley,et al.  Estimating ancestral population parameters. , 1997, Genetics.

[55]  Jared L. Strasburg,et al.  How robust are "isolation with migration" analyses to violations of the im model? A simulation study. , 2010, Molecular biology and evolution.

[56]  J. Piálek,et al.  GENETIC ANALYSIS OF AUTOSOMAL AND X-LINKED MARKERS ACROSS A MOUSE HYBRID ZONE , 2007, Evolution; international journal of organic evolution.

[57]  P. Green,et al.  Consed: a graphical tool for sequence finishing. , 1998, Genome research.

[58]  N. Saitou,et al.  The neighbor-joining method: a new method for reconstructing phylogenetic trees. , 1987, Molecular biology and evolution.

[59]  F. Bronson,et al.  The Reproductive Ecology of the House Mouse , 1979, The Quarterly Review of Biology.

[60]  P. Boursot,et al.  Postzygotic isolation between the two European subspecies of the house mouse: estimates from fertility patterns in wild and laboratory-bred hybrids , 2005 .

[61]  Maria A Sans-Fuentes,et al.  Genome-wide patterns of gene flow across a house mouse hybrid zone. , 2007, Genome research.

[62]  F. Bonhomme,et al.  THE EVOLUTION OF HOUSE MICE , 1993 .

[63]  M. Noor,et al.  Islands of speciation or mirages in the desert? Examining the role of restricted recombination in maintaining species , 2010, Heredity.

[64]  P. Donnelly,et al.  A new statistical method for haplotype reconstruction from population data. , 2001, American journal of human genetics.

[65]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[66]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[67]  W. Li,et al.  Statistical tests of neutrality of mutations. , 1993, Genetics.

[68]  B. Charlesworth.  Measures of divergence between populations and the effect of forces that reduce variability. , 1998, Molecular biology and evolution.