Inferring the history of speciation in house mice from autosomal, X-linked, Y-linked and mitochondrial genes

A pronounced evolutionary shift of the pseudoautosomal region boundary in house mice

The pseudoautosomal region (PAR) is essential for the accurate pairing and segregation of the X and Y chromosomes during meiosis. Despite its functional significance, the PAR shows substantial evolutionary divergence in structure and sequence between mammalian species. An instructive example of PAR evolution is the house mouse Mus musculus domesticus (represented by the C57BL/6J strain), which has the smallest PAR among those that have been mapped. In C57BL/6J, the PAR boundary is located just ~700 kb from the distal end of the X chromosome, whereas the boundary is found at a more proximal position in Mus spretus, a species that diverged from house mice 2-4 million years ago. In this study we used a combination of genetic and physical mapping to document a pronounced shift in the PAR boundary in a second house mouse subspecies, Mus musculus castaneus (represented by the CAST/EiJ strain), ~430 kb proximal of the M. m. domesticus boundary. We demonstrate molecular evolutionary consequences of this shift, including a marked lineage-specific increase in sequence divergence within Mid1, a gene that resides entirely within the M. m. castaneus PAR but straddles the boundary in other subspecies. Our results extend observations of structural divergence in the PAR to closely related subspecies, pointing to major evolutionary changes in this functionally important genomic region over a short time period.

Genetics and evolution of hybrid male sterility in house mice

Comparative genetic mapping provides insights into the evolution of the reproductive barriers that separate closely related species. This approach has been used to document the accumulation of reproductive incompatibilities over time, but has only been applied to a few taxa. House mice offer a powerful system to reconstruct the evolution of reproductive isolation between multiple subspecies pairs. However, studies of the primary reproductive barrier in house mice-hybrid male sterility-have been restricted to a single subspecies pair: Mus musculus musculus and Mus musculus domesticus. To provide a more complete characterization of reproductive isolation in house mice, we conducted an F(2) intercross between wild-derived inbred strains from Mus musculus castaneus and M. m. domesticus. We identified autosomal and X-linked QTL associated with a range of hybrid male sterility phenotypes, including testis weight, sperm density, and sperm morphology. The pseudoautosomal region (PAR) was strongly associated with hybrid sterility phenotypes when heterozygous. We compared QTL found in this cross with QTL identified in a previous F(2) intercross between M. m. musculus and M. m. domesticus and found three shared autosomal QTL. Most QTL were not shared, demonstrating that the genetic basis of hybrid male sterility largely differs between these closely related subspecies pairs. These results lay the groundwork for identifying genes responsible for the early stages of speciation in house mice.

Genome-wide architecture of reproductive isolation in a naturally occurring hybrid zone between Mus musculus musculus and M. m. domesticus

Studies of a hybrid zone between two house mouse subspecies (Mus musculus musculus and M. m. domesticus) along with studies using laboratory crosses reveal a large role for the X chromosome and multiple autosomal regions in reproductive isolation as a consequence of disrupted epistasis in hybrids. One limitation of previous work has been that most of the identified genomic regions have been large. The goal here is to detect and characterize precise genomic regions underlying reproductive isolation. We surveyed 1401 markers evenly spaced across the genome in 679 mice collected from two different transects. Comparisons between transects provide a means for identifying common patterns that likely reflect intrinsic incompatibilities. We used a genomic cline approach to identify patterns that correspond to epistasis. From both transects, we identified contiguous regions on the X chromosome in which markers were inferred to be involved in epistatic interactions. We then searched for autosomal regions showing the same patterns and found they constitute about 5% of autosomal markers. We discovered substantial overlap between these candidate regions underlying reproductive isolation and QTL for hybrid sterility identified in laboratory crosses. Analysis of gene content in these regions suggests a key role for several mechanisms, including the regulation of transcription, sexual conflict and sexual selection operating at both the postmating prezygotic and postzygotic stages of reproductive isolation. Taken together, these results indicate that speciation in two recently diverged (c. 0.5 Ma) house mouse subspecies is complex, involving many genes dispersed throughout the genome and associated with distinct functions.

Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice

Recently diverged taxa may continue to exchange genes. A number of models of speciation with gene flow propose that the frequency of gene exchange will be lower in genomic regions of low recombination and that these regions will therefore be more differentiated. However, several population-genetic models that focus on selection at linked sites also predict greater differentiation in regions of low recombination simply as a result of faster sorting of ancestral alleles even in the absence of gene flow. Moreover, identifying the actual amount of gene flow from patterns of genetic variation is tricky, because both ancestral polymorphism and migration lead to shared variation between recently diverged taxa. New analytic methods have been developed to help distinguish ancestral polymorphism from migration. Along with a growing number of datasets of multi-locus DNA sequence variation, these methods have spawned a renewed interest in speciation models with gene flow. Here, we review both speciation and population-genetic models that make explicit predictions about how the rate of recombination influences patterns of genetic variation within and between species. We then compare those predictions with empirical data of DNA sequence variation in rabbits and mice. We find strong support for the prediction that genomic regions experiencing low levels of recombination are more differentiated. In most cases, reduced gene flow appears to contribute to the pattern, although disentangling the relative contribution of reduced gene flow and selection at linked sites remains a challenge. We suggest fruitful areas of research that might help distinguish between different models.

Where are the wormy mice? A reexamination of hybrid parasitism in the European house mouse hybrid zone

Wormy mice in a hybrid zone have been interpreted as evidence of low hybrid fitness, such that parasites contribute to species separation. However, because of its natural heterogeneity, observations of parasite load must be numerous with good field area coverage. We sampled 689 mice from 107 localities across the Bavaria-Bohemia region of the European house mouse hybrid zone and calculated their hybrid indices using 1401 diagnostic single nucleotide polymorphisms (SNPs). We tested whether hybrids have greater or lesser diversity and load of parasite helminths than additive expectations, performing load analyses on the four most common taxa. We found hybrids have significantly reduced diversity and load of each of the commonest helminths; rarer helminths further support reduced load. Although within-locality comparisons have little power, randomization tests show the repeated pattern is unlikely to be due to local parasite heterogeneity, and simulations show a patch of low parasite diversity is unlikely to fall by chance just so in the field area, such that it produces the observed effects. Our data therefore contradict the idea that helminths reduce hybrid fitness through increased load. We discuss a vicariant Red Queen model that implies immune genes tracking parasites will escape Dobzhansky-Muller incompatibilities, generating hybrid variants untargeted by parasites.

Adaptive evolution and effective population size in wild house mice

Estimates of the proportion of amino acid substitutions that have been fixed by selection (α) vary widely among taxa, ranging from zero in humans to over 50% in Drosophila. This wide range may reflect differences in the efficacy of selection due to differences in the effective population size (N(e)). However, most comparisons have been made among distantly related organisms that differ not only in N(e) but also in many other aspects of their biology. Here, we estimate α in three closely related lineages of house mice that have a similar ecology but differ widely in N(e): Mus musculus musculus (N(e) ∼ 25,000-120,000), M. m. domesticus (N(e) ∼ 58,000-200,000), and M. m. castaneus (N(e) ∼ 200,000-733,000). Mice were genotyped using a high-density single nucleotide polymorphism array, and the proportions of replacement and silent mutations within subspecies were compared with those fixed between each subspecies and an outgroup, Mus spretus. There was significant evidence of positive selection in M. m. castaneus, the lineage with the largest N(e), with α estimated to be approximately 40%. In contrast, estimates of α for M. m. domesticus (α = 13%) and for M. m. musculus (α = 12 %) were much smaller. Interestingly, the higher estimate of α for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection. These results support the hypothesis that differences in N(e) contribute to differences among species in the efficacy of selection.

Statistical guidelines for detecting past population shifts using ancient DNA

Populations carry a genetic signal of their demographic past, providing an opportunity for investigating the processes that shaped their evolution. Our ability to infer population histories can be enhanced by including ancient DNA data. Using serial-coalescent simulations and a range of both quantitative and temporal sampling schemes, we test the power of ancient mitochondrial sequences and nuclear single-nucleotide polymorphisms (SNPs) to detect past population bottlenecks. Within our simulated framework, mitochondrial sequences have only limited power to detect subtle bottlenecks and/or fast post-bottleneck recoveries. In contrast, nuclear SNPs can detect bottlenecks followed by rapid recovery, although bottlenecks involving reduction of less than half the population are generally detected with low power unless extensive genetic information from ancient individuals is available. Our results provide useful guidelines for scaling sampling schemes and for optimizing our ability to infer past population dynamics. In addition, our results suggest that many ancient DNA studies may face power issues in detecting moderate demographic collapses and/or highly dynamic demographic shifts when based solely on mitochondrial information.

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.

Chromosomal dynamics of nucleolar organizer regions (NORs) in the house mouse: micro-evolutionary insights

Variation in the number and chromosomal location of nucleolar organizer regions (NORs) was studied in the house mouse, Mus musculus (2n=40). From an origin in Western Asia, this species colonized the Middle East, Europe and Asia. This expansion was accompanied by diversification into five subspecies. NOR diversity was revealed by fluorescence in situ hybridization using 18S and 28S probes on specimens spanning Asia to Western Europe. The results showed that the house mouse genome possessed a large number of NOR-bearing autosomes and a surprisingly high rate of polymorphism for the presence/absence of rRNA genes on all these chromosomes. All NOR sites were adjacent to the centromere except for two that were telomeric. Subspecific differentiation established from the NOR frequency data was concordant with the overall pattern of radiation proposed from molecular studies, but highlighted several discrepancies that need to be further addressed. NOR diversity in M. musculus consisted of a large number of polymorphic NORs that were common to at least two subspecies, and a smaller number of NORs that were unique to one subspecies. The most parsimonious scenario argues in favor of a subspecific differentiation by lineage sorting of ancestral NOR polymorphisms; only the unique NORs would have appeared by inter-chromosomal transposition, except for the two telomeric ones that may have originated by hybridization with another species. Such a scenario provides an alternative view from the one prevailing in most systematic and phylogenetic analyses that NORs have a high transposition rate due to concerted evolution of rRNA genes.

Isolation and gene flow: inferring the speciation history of European house mice

Inferring the history of isolation and gene flow during species differentiation can inform us on the processes underlying their formation. Following their recent expansion in Europe, two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus) have formed a hybrid zone maintained by hybrid incompatibilities and possibly behavioural reinforcement, offering a good model of incipient speciation. We reconstruct the history of their divergence using an approximate Bayesian computation framework and sequence variation at 57 autosomal loci. We find support for a long isolation period preceding the advent of gene flow around 200,000 generations ago, much before the formation of the European hybrid zone a few thousand years ago. The duration of the allopatric episode appears long enough (74% of divergence time) to explain the accumulation of many post-zygotic incompatibilities expressed in the present hybrid zone. The ancient contact inferred could have played a role in mating behaviour divergence and laid the ground for further reinforcement. We suggest that both subspecies originally colonized the Middle East from the northern Indian subcontinent, domesticus settling on the shores of the Persian Gulf and musculus on those of the Caspian Sea. Range expansions during interglacials would have induced secondary contacts, presumably in Iran, where they must have also interacted with Mus musculus castaneus. Future studies should incorporate this possibility, and we point to Iran and its surroundings as a hot spot for house mouse diversity and speciation studies.

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.

A comparative assessment of mandible shape in a consomic strain panel of the house mouse (Mus musculus)--implications for epistasis and evolvability of quantitative traits

Background

Expectations of repeatedly finding associations between given genes and phenotypes have been borne out by studies of parallel evolution, especially for traits involving absence or presence of characters. However, it has rarely been asked whether the genetic basis of quantitative trait variation is conserved at the intra- or even at the interspecific level. This question is especially relevant for shape, where the high dimensionality of variation seems to require a highly complex genetic architecture involving many genes.

Results

We analyse here the genetic effects of chromosome substitution strains carrying M. m. musculus chromosomes in a largely M. m. domesticus background on mandible shape and compare them to the results of previously published QTL mapping data between M. m. domesticus strains. We find that the distribution of genetic effects and effect sizes across the genome is consistent between the studies, while the specific shape changes associated with the chromosomes are different. We find also that the sum of the effects from the different M. m. musculus chromosomes is very different from the shape of the strain from which they were derived, as well as all known wild type shapes.

Conclusions

Our results suggest that the relative chromosome-wide effect sizes are comparable between the long separated subspecies M. m. domesticus and M. m. musculus, hinting at a relative stability of genes involved in this complex trait. However, the absolute effect sizes and the effect directions may be allele-dependent, or are context dependent, i.e. epistatic interactions appear to play an important role in controlling shape.

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 (N(e) ) 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, N(e) for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest N(e) and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest N(e) 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. F(ST) 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.

Accelerated evolution of sex chromosomes in aphids, an x0 system

Sex chromosomes play a role in many important biological processes, including sex determination, genomic conflicts, imprinting, and speciation. In particular, they exhibit several unusual properties such as inheritance pattern, hemizygosity, and reduced recombination, which influence their response to evolutionary factors (e.g., drift, selection, and demography). Here, we examine the evolutionary forces driving X chromosome evolution in aphids, an XO system where females are homozygous (XX) and males are hemizygous (X0) at sex chromosomes. We show by simulations that the unusual mode of transmission of the X chromosome in aphids, coupled with cyclical parthenogenesis, results in similar effective population sizes and predicted levels of genetic diversity for X chromosomes and autosomes under neutral evolution. These results contrast with expectations from standard XX/XY or XX/X0 systems (where the effective population size of the X is three-fourths that of autosomes) and have deep consequences for aphid X chromosome evolution. We then localized 52 microsatellite markers on the X and 351 on autosomes. We genotyped 167 individuals with 356 of these loci and found similar levels of allelic richness on the X and on the autosomes, as predicted by our simulations. In contrast, we detected higher dN and dN/dS ratio for X-linked genes compared with autosomal genes, a pattern compatible with either positive or relaxed selection. Given that both types of chromosomes have similar effective population sizes and that the single copy of the X chromosome of male aphids exposes its recessive genes to selection, some degree of positive selection seems to best explain the higher rates of evolution of X-linked genes. Overall, this study highlights the particular relevance of aphids to study the evolutionary factors driving sex chromosomes and genome evolution.

Joint inference of population assignment and demographic history

A new approach to assigning individuals to populations using genetic data is described. Most existing methods work by maximizing Hardy-Weinberg and linkage equilibrium within populations, neither of which will apply for many demographic histories. By including a demographic model, within a likelihood framework based on coalescent theory, we can jointly study demographic history and population assignment. Genealogies and population assignments are sampled from a posterior distribution using a general isolation-with-migration model for multiple populations. A measure of partition distance between assignments facilitates not only the summary of a posterior sample of assignments, but also the estimation of the posterior density for the demographic history. It is shown that joint estimates of assignment and demographic history are possible, including estimation of population phylogeny for samples from three populations. The new method is compared to results of a widely used assignment method, using simulated and published empirical data sets.

Genetic dissection of a key reproductive barrier between nascent species of house mice

Reproductive isolation between species is often caused by deleterious interactions among loci in hybrids. Finding the genes involved in these incompatibilities provides insight into the mechanisms of speciation. With recently diverged subspecies, house mice provide a powerful system for understanding the genetics of reproductive isolation early in the speciation process. Although previous studies have yielded important clues about the genetics of hybrid male sterility in house mice, they have been restricted to F1 sterility or incompatibilities involving the X chromosome. To provide a more complete characterization of this key reproductive barrier, we conducted an F2 intercross between wild-derived inbred strains from two subspecies of house mice, Mus musculus musculus and Mus musculus domesticus. We identified a suite of autosomal and X-linked QTL that underlie measures of hybrid male sterility, including testis weight, sperm density, and sperm morphology. In many cases, the autosomal loci were unique to a specific sterility trait and exhibited an effect only when homozygous, underscoring the importance of examining reproductive barriers beyond the F1 generation. We also found novel two-locus incompatibilities between the M. m. musculus X chromosome and M. m. domesticus autosomal alleles. Our results reveal a complex genetic architecture for hybrid male sterility and suggest a prominent role for reproductive barriers in advanced generations in maintaining subspecies integrity in house mice.

Bonobos fall within the genomic variation of chimpanzees

To gain insight into the patterns of genetic variation and evolutionary relationships within and between bonobos and chimpanzees, we sequenced 150,000 base pairs of nuclear DNA divided among 15 autosomal regions as well as the complete mitochondrial genomes from 20 bonobos and 58 chimpanzees. Except for western chimpanzees, we found poor genetic separation of chimpanzees based on sample locality. In contrast, bonobos consistently cluster together but fall as a group within the variation of chimpanzees for many of the regions. Thus, while chimpanzees retain genomic variation that predates bonobo-chimpanzee speciation, extensive lineage sorting has occurred within bonobos such that much of their genome traces its ancestry back to a single common ancestor that postdates their origin as a group separate from chimpanzees.

Genetic analysis of genome-scale recombination rate evolution in house mice

The rate of meiotic recombination varies markedly between species and among individuals. Classical genetic experiments demonstrated a heritable component to population variation in recombination rate, and specific sequence variants that contribute to recombination rate differences between individuals have recently been identified. Despite these advances, the genetic basis of species divergence in recombination rate remains unexplored. Using a cytological assay that allows direct in situ imaging of recombination events in spermatocytes, we report a large (∼30%) difference in global recombination rate between males of two closely related house mouse subspecies (Mus musculus musculus and M. m. castaneus). To characterize the genetic basis of this recombination rate divergence, we generated an F2 panel of inter-subspecific hybrid males (n = 276) from an intercross between wild-derived inbred strains CAST/EiJ (M. m. castaneus) and PWD/PhJ (M. m. musculus). We uncover considerable heritable variation for recombination rate among males from this mapping population. Much of the F2 variance for recombination rate and a substantial portion of the difference in recombination rate between the parental strains is explained by eight moderate- to large-effect quantitative trait loci, including two transgressive loci on the X chromosome. In contrast to the rapid evolution observed in males, female CAST/EiJ and PWD/PhJ animals show minimal divergence in recombination rate (∼5%). The existence of loci on the X chromosome suggests a genetic mechanism to explain this male-biased evolution. Our results provide an initial map of the genetic changes underlying subspecies differences in genome-scale recombination rate and underscore the power of the house mouse system for understanding the evolution of this trait.

Homologous recombination is an indispensable feature of the mammalian meiotic program and an important mechanism for creating genetic diversity. Despite its central significance, recombination rates vary markedly between species and among individuals. Although recent studies have begun to unravel the genetic basis of recombination rate variation within populations, the genetic mechanisms of species divergence in recombination rate remain poorly characterized. In this study, we show that two closely related house mouse subspecies differ in their genomic recombination rates by ∼30%, providing an excellent model system for studying evolutionary divergence in this trait. Using quantitative genetic methods, we identify eight genomic regions that contribute to divergence in global recombination rate between these subspecies, including large effect loci and multiple loci on the X-chromosome. Our study uncovers novel genomic loci contributing to species divergence in global recombination rate and offers simple genetic explanations for rapid phenotypic divergence in this trait.

Subspecific origin and haplotype diversity in the laboratory mouse

Here we provide a genome-wide, high-resolution map of the phylogenetic origin of the genome of most extant laboratory mouse inbred strains. Our analysis is based on the genotypes of wild-caught mice from three subspecies of Mus musculus. We show that classical laboratory strains are derived from a few fancy mice with limited haplotype diversity. Their genomes are overwhelmingly Mus musculus domesticus in origin, and the remainder is mostly of Japanese origin. We generated genome-wide haplotype maps based on identity by descent from fancy mice and show that classical inbred strains have limited and non-randomly distributed genetic diversity. In contrast, wild-derived laboratory strains represent a broad sampling of diversity within M. musculus. Intersubspecific introgression is pervasive in these strains, and contamination by laboratory stocks has played a role in this process. The subspecific origin, haplotype diversity and identity by descent maps can be visualized using the Mouse Phylogeny Viewer (see URLs).

Speciation in the European rabbit (Oryctolagus cuniculus): islands of differentiation on the X chromosome and autosomes

Studies of gene flow between closely related taxa can provide insight into the genetic basis of speciation. To evaluate the importance of the X chromosome in reproductive isolation between subspecies of the European rabbit and to study the genomic scale over which islands of differentiation extend, we resequenced a total of 34 loci distributed along the X chromosome and chromosome 14. Previous studies based on few markers suggested that loci in centromeric regions were highly differentiated between rabbit subspecies, whereas loci in telomeric regions were less differentiated. Here, we confirmed this finding but also discovered remarkable variation in levels of differentiation among loci, with F(ST) values from nearly 0 to 1. Analyses using isolation-with-migration models suggest that this range appears to be largely explained by differential levels of gene flow among loci. The X chromosome was significantly more differentiated than the autosomes. On chromosome 14, differentiation decayed very rapidly at increasing distances from the centromere, but on the X chromosome distinct islands of differentiation encompassing several megabases were observed both at the centromeric region and along the chromosome arms. These findings support the idea that the X chromosome plays an important role in reproductive isolation between rabbit subspecies. These results also demonstrate the mosaic nature of the genome at species boundaries.

Systematic detection of polygenic cis-regulatory evolution

The idea that most morphological adaptations can be attributed to changes in the cis-regulation of gene expression levels has been gaining increasing acceptance, despite the fact that only a handful of such cases have so far been demonstrated. Moreover, because each of these cases involves only one gene, we lack any understanding of how natural selection may act on cis-regulation across entire pathways or networks. Here we apply a genome-wide test for selection on cis-regulation to two subspecies of the mouse Mus musculus. We find evidence for lineage-specific selection at over 100 genes involved in diverse processes such as growth, locomotion, and memory. These gene sets implicate candidate genes that are supported by both quantitative trait loci and a validated causality-testing framework, and they predict a number of phenotypic differences, which we confirm in all four cases tested. Our results suggest that gene expression adaptation is widespread and that these adaptations can be highly polygenic, involving cis-regulatory changes at numerous functionally related genes. These coordinated adaptations may contribute to divergence in a wide range of morphological, physiological, and behavioral phenotypes.

Darwin and the origin of interspecific genetic incompatibilities

Darwin's Origin of Species is often criticized for having little to say about speciation. The complaint focuses in particular on Darwin's supposed failure to explain the evolution of the sterility and inviability of interspecific hybrids. But in his chapter on hybridism, Darwin, working without genetics, got as close to the modern understanding of the evolution of hybrid sterility and inviability as might reasonably be expected. In particular, after surveying what was then known about interspecific crosses and the resulting hybrids, he established two facts that, while now taken for granted, were at the time radical. First, the sterility barriers between species are neither specially endowed by a creator nor directly favored by natural selection but rather evolve as incidental by-products of interspecific divergence. Second, the sterility of species hybrids results when their development is "disturbed by two organizations having been compounded into one." Bateson, Dobzhansky, and Muller later put Mendelian detail to Darwin's inference that the species-specific factors controlling development (i.e., genes) are sometimes incompatible. In this article, I highlight the major developments in our understanding of these interspecific genetic incompatibilities--from Darwin to Muller to modern theory--and review comparative, genetic, and molecular rules that characterize the evolution of hybrid sterility and inviability.

Differential expression of non-coding RNAs and continuous evolution of the X chromosome in testicular transcriptome of two mouse species

BACKGROUND

Tight regulation of testicular gene expression is a prerequisite for male reproductive success, while differentiation of gene activity in spermatogenesis is important during speciation. Thus, comparison of testicular transcriptomes between closely related species can reveal unique regulatory patterns and shed light on evolutionary constraints separating the species.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we compared testicular transcriptomes of two closely related mouse species, Mus musculus and Mus spretus, which diverged more than one million years ago. We analyzed testicular expression using tiling arrays overlapping Chromosomes 2, X, Y and mitochondrial genome. An excess of differentially regulated non-coding RNAs was found on Chromosome 2 including the intronic antisense RNAs, intergenic RNAs and premature forms of Piwi-interacting RNAs (piRNAs). Moreover, striking difference was found in the expression of X-linked G6pdx gene, the parental gene of the autosomal retrogene G6pd2.

CONCLUSIONS/SIGNIFICANCE

The prevalence of non-coding RNAs among differentially expressed transcripts indicates their role in species-specific regulation of spermatogenesis. The postmeiotic expression of G6pdx in Mus spretus points towards the continuous evolution of X-chromosome silencing and provides an example of expression change accompanying the out-of-the X-chromosomal retroposition.

Inference of selection and stochastic effects in the house mouse hybrid zone

We explored the transition of 13 X-linked markers across two separate portions of the house mouse hybrid zone, asking whether such a comparison can distinguish the effects of selection from random factors. A heuristic search in the likelihood landscape revealed more complex likelihood profiles for data sampled in two-dimensional (2D) space relative to data sampled along a linear transect. Randomized resampling of localities analyzed for individual loci showed that deletion of sites away from the zone center can decrease cline width estimates whereas deletion of sites close to the center can significantly increase the width estimates. Deleting localities for all loci resulted in wider clines if the number of samples from the center was limited. The results suggest that, given the great variation in width estimates resulting from inclusion/exclusion of sampling sites, the geographic sampling design is important in hybrid zone studies and that our inferences should take into account measures of uncertainty such as support intervals. The comparison of the two transects indicates cline widths are narrower for loci in the central part of the X chromosome, suggesting selection is stronger in this region and genetic incompatibilities may have at least partly common architecture in the house mouse hybrid zone.

Extraordinary sequence divergence at Tsga8, an X-linked gene involved in mouse spermiogenesis

The X chromosome plays an important role in both adaptive evolution and speciation. We used a molecular evolutionary screen of X-linked genes potentially involved in reproductive isolation in mice to identify putative targets of recurrent positive selection. We then sequenced five very rapidly evolving genes within and between several closely related species of mice in the genus Mus. All five genes were involved in male reproduction and four of the genes showed evidence of recurrent positive selection. The most remarkable evolutionary patterns were found at Testis-specific gene a8 (Tsga8), a spermatogenesis-specific gene expressed during postmeiotic chromatin condensation and nuclear transformation. Tsga8 was characterized by extremely high levels of insertion-deletion variation of an alanine-rich repetitive motif in natural populations of Mus domesticus and M. musculus, differing in length from the reference mouse genome by up to 89 amino acids (27% of the total protein length). This population-level variation was coupled with striking divergence in protein sequence and length between closely related mouse species. Although no clear orthologs had previously been described for Tsga8 in other mammalian species, we have identified a highly divergent hypothetical gene on the rat X chromosome that shares clear orthology with the 5' and 3' ends of Tsga8. Further inspection of this ortholog verified that it is expressed in rat testis and shares remarkable similarity with mouse Tsga8 across several general features of the protein sequence despite no conservation of nucleotide sequence across over 60% of the rat-coding domain. Overall, Tsga8 appears to be one of the most rapidly evolving genes to have been described in rodents. We discuss the potential evolutionary causes and functional implications of this extraordinary divergence and the possible contribution of Tsga8 and the other four genes we examined to reproductive isolation in mice.

The mouse "xenotropic" gammaretroviruses and their XPR1 receptor

The xenotropic/polytropic subgroup of mouse leukemia viruses (MLVs) all rely on the XPR1 receptor for entry, but these viruses vary in tropism, distribution among wild and laboratory mice, pathogenicity, strategies used for transmission, and sensitivity to host restriction factors. Most, but not all, isolates have typical xenotropic or polytropic host range, and these two MLV tropism types have now been detected in humans as viral sequences or as infectious virus, termed XMRV, or xenotropic murine leukemia virus-related virus. The mouse xenotropic MLVs (X-MLVs) were originally defined by their inability to infect cells of their natural mouse hosts. It is now clear, however, that X-MLVs actually have the broadest host range of the MLVs. Nearly all nonrodent mammals are susceptible to X-MLVs, and all species of wild mice and several common strains of laboratory mice are X-MLV susceptible. The polytropic MLVs, named for their apparent broad host range, show a more limited host range than the X-MLVs in that they fail to infect cells of many mouse species as well as many nonrodent mammals. The co-evolution of these viruses with their receptor and other host factors that affect their replication has produced a heterogeneous group of viruses capable of inducing various diseases, as well as endogenized viral genomes, some of which have been domesticated by their hosts to serve in antiviral defense.

Norwegian house mice (Mus musculus musculus/domesticus): distributions, routes of colonization and patterns of hybridization

We investigated the distributions and routes of colonization of two commensal subspecies of house mouse in Norway: Mus musculus domesticus and M. m. musculus. Five nuclear markers (Abpa, D11 cenB2, Btk, SMCY and Zfy2) and a morphological feature (tail length) were used to differentiate the two subspecies and assess their distributions, and mitochondrial (mt) D-loop sequences helped to elucidate their colonization history. M. m. domesticus is the more widespread of the two subspecies, occupying the western and southern coast of Norway, while M. m. musculus is found along Norway's southeastern coast and east from there to Sweden. Two sections of the hybrid zone between the two subspecies were localized in Norway. However, hybrid forms also occur well away from that hybrid zone, the most prevalent of which are mice with a M. m. musculus-type Y chromosome and an otherwise M. m. domesticus genome. MtDNA D-loop sequences of the mice revealed a complex phylogeography within M. m. domesticus, reflecting passive human transport to Norway, probably during the Viking period. M. m. musculus may have colonized earlier. If so, that leaves open the possibility that M. m. domesticus replaced M. m. musculus from much of Norway, with the widely distributed hybrids a relict of this process. Overall, the effects of hybridization are evident in house mice throughout Norway.

House mouse colonization patterns on the sub-Antarctic Kerguelen Archipelago suggest singular primary invasions and resilience against re-invasion

BACKGROUND

Starting from Western Europe, the house mouse (Mus musculus domesticus) has spread across the globe in historic times. However, most oceanic islands were colonized by mice only within the past 300 years. This makes them an excellent model for studying the evolutionary processes during early stages of new colonization. We have focused here on the Kerguelen Archipelago, located within the sub-Antarctic area and compare the patterns with samples from other Southern Ocean islands.

RESULTS

We have typed 18 autosomal and six Y-chromosomal microsatellite loci and obtained mitochondrial D-loop sequences for a total of 534 samples, mainly from the Kerguelen Archipelago, but also from the Falkland Islands, Marion Island, Amsterdam Island, Antipodes Island, Macquarie Island, Auckland Islands and one sample from South Georgia. We find that most of the mice on the Kerguelen Archipelago have the same mitochondrial haplotype and all share the same major Y-chromosomal haplotype. Two small islands (Cochons Island and Cimetière Island) within the archipelago show a different mitochondrial haplotype, are genetically distinct for autosomal loci, but share the major Y-chromosomal haplotype. In the mitochondrial D-loop sequences, we find several single step mutational derivatives of one of the major mitochondrial haplotypes, suggesting an unusually high mutation rate, or the occurrence of selective sweeps in mitochondria.

CONCLUSIONS

Although there was heavy ship traffic for over a hundred years to the Kerguelen Archipelago, it appears that the mice that have arrived first have colonized the main island (Grande Terre) and most of the associated small islands. The second invasion that we see in our data has occurred on islands that are detached from Grande Terre and were likely to have had no resident mice prior to their arrival. The genetic data suggest that the mice of both primary invasions originated from related source populations. Our data suggest that an area colonized by mice is refractory to further introgression, possibly due to fast adaptations of the resident mice to local conditions.

Extensive recombination rate variation in the house mouse species complex inferred from genetic linkage maps

The rate of recombination is a key genomic parameter that displays considerable variation among taxa. Species comparisons have demonstrated that the rate of evolution in recombination rate is strongly dependent on the physical scale of measurement. Individual recombination hotspots are poorly conserved among closely related taxa, whereas genomic-scale recombination rate variation bears a strong signature of phylogenetic history. In contrast, the mode and tempo of evolution in recombination rates measured on intermediate physical scales is poorly understood. Here, we conduct a detailed statistical comparison between two whole-genome F₂ genetic linkage maps constructed from experimental intercrosses between closely related house mouse subspecies (Mus musculus). Our two maps profile a common wild-derived inbred strain of M. m. domesticus crossed to distinct wild-derived inbred strains representative of two other house mouse subspecies, M. m. castaneus and M. m. musculus. We identify numerous orthologous genomic regions with significant map length differences between these two crosses. Because the genomes of these recently diverged house mice are highly collinear, observed differences in map length (centimorgans) are suggestive of variation in broadscale recombination rate (centimorgans per megabase) within M. musculus. Collectively, these divergent intervals span 19% of the house mouse genome, disproportionately aggregating on the X chromosome. In addition, we uncover strong statistical evidence for a large effect, sex-linked, site-specific modifier of recombination rate segregating within M. musculus. Our findings reveal considerable variation in the megabase-scale recombination landscape among recently diverged taxa and underscore the continued importance of genetic linkage maps in the post-genome era.

Widespread over-expression of the X chromosome in sterile F₁hybrid mice

The X chromosome often plays a central role in hybrid male sterility between species, but it is unclear if this reflects underlying regulatory incompatibilities. Here we combine phenotypic data with genome-wide expression data to directly associate aberrant expression patterns with hybrid male sterility between two species of mice. We used a reciprocal cross in which F₁ males are sterile in one direction and fertile in the other direction, allowing us to associate expression differences with sterility rather than with other hybrid phenotypes. We found evidence of extensive over-expression of the X chromosome during spermatogenesis in sterile but not in fertile F₁ hybrid males. Over-expression was most pronounced in genes that are normally expressed after meiosis, consistent with an X chromosome-wide disruption of expression during the later stages of spermatogenesis. This pattern was not a simple consequence of faster evolutionary divergence on the X chromosome, because X-linked expression was highly conserved between the two species. Thus, transcriptional regulation of the X chromosome during spermatogenesis appears particularly sensitive to evolutionary divergence between species. Overall, these data provide evidence for an underlying regulatory basis to reproductive isolation in house mice and underscore the importance of transcriptional regulation of the X chromosome to the evolution of hybrid male sterility.

The X chromosome plays an important role in the development of reproductive isolation between species, but the basis for this has remained unclear. One possible explanation is that sperm development is sensitive to disruption of X-linked gene regulation. In mice, evidence linking abnormal gene expression on the X chromosome with reproductive isolation has been lacking until now. Here we use experimental crosses within and between species of mice and genome-wide expression data to identify aberrant expression patterns associated with hybrid male sterility. We observed chromosome-wide over-expression of the X chromosome during spermatogenesis in sterile hybrid males and developmentally localized this breakdown to an apparent disruption of X-inactivation. Collectively, these results highlight the importance of gene regulation to the evolution of reproductive isolation and support the hypothesis that improper expression of the X chromosome during spermatogenesis is an important mechanism contributing to the rapid evolution of hybrid male sterility.

Genetic differentiation of the house mouse around the Mediterranean basin: matrilineal footprints of early and late colonization

The molecular signatures of the recent expansion of the western house mouse, Mus musculus domesticus, around the Mediterranean basin are investigated through the study of mitochondrial D-loop polymorphism on a 1313 individual dataset. When reducing the complexity of the matrilineal network to a series of haplogroups (HGs), our main results indicate that: (i) several HGs are recognized which seem to have almost simultaneously diverged from each other, confirming a recent expansion for the whole subspecies; (ii) some HGs are geographically delimited while others are widespread, indicative of multiple introductions or secondary exchanges; (iii) mice from the western and the eastern coasts of Africa harbour largely different sets of HGs; and (iv) HGs from the two shores of the Mediterranean are more similar in the west than in the east. This pattern is in keeping with the two-step westward expansion proposed by zooarchaeological data, an early one coincident with the Neolithic progression and limited to the eastern Mediterranean and a later one, particularly evident in the western Mediterranean, related to the generalization of maritime trade during the first millennium BC and onwards. The dispersal of mice along with humans, which continues until today, has for instance left complex footprints on the long ago colonized Cyprus or more simple ones on the much more recently populated Canary Islands.

Evolution of functional and sequence variants of the mammalian XPR1 receptor for mouse xenotropic gammaretroviruses and the human-derived retrovirus XMRV

Genetic conflicts between retroviruses and their receptors result in the evolution of novel host entry restrictions and novel virus envelopes, and such variants can influence trans-species transmission. We screened rodents and other mammals for sequence variation in the Xpr1 receptor for the mouse xenotropic or polytropic mouse leukemia viruses (X-MLVs or P-MLVs, respectively) of the gammaretrovirus family and for susceptibility to mouse-derived X/P-MLVs and to XMRV (xenotropic murine leukemia virus-related virus), an X-MLV-like virus isolated from humans with prostate cancer and chronic fatigue syndrome. We identified multiple distinct susceptibility phenotypes; these include the four known Xpr1 variants in Mus and a novel fifth Xpr1 gene found in Mus molossinus and Mus musculus. We describe the geographic and species distribution of the Mus Xpr1 variants but failed to find the X-MLV-restrictive laboratory mouse allele in any wild mouse. We used mutagenesis and phylogenetic analysis to evaluate the functional contributions made by constrained, variable, and deleted residues. Rodent Xpr1 is under positive selection, indicating a history of host-pathogen conflicts; several codons under selection have known roles in virus entry. All non-Mus mammals are susceptible to mouse X-MLVs, but some restrict other members of the X/P-MLV family, and the resistance of hamster and gerbil cells to XMRV indicates that XMRV has unique receptor requirements. We show that the hypervariable fourth extracellular XPR1 loop (ECL4) contains three evolutionarily constrained residues that do not contribute to receptor function, we identify two novel residues important for virus entry (I579 and T583), and we describe a unique pattern of ECL4 variation in the three virus-restrictive Xpr1 variants found in MLV-infected house mice; these mice carry different deletions in ECL4, suggesting either that these sites or loop size affects receptor function.

Maintenance of fungal pathogen species that are specialized to different hosts: allopatric divergence and introgression through secondary contact

Sympatry of species that lack complete prezygotic isolation is ideal for the study of how species can be maintained in the face of potential gene flow. This is particularly important in the context of emerging diseases on new hosts because pathogen adaptation is facilitated by reduced gene flow from ancestral populations. Here, we investigated divergence and gene flow between two closely related fungal species, Microbotryum lychnidis-dioicae and M. silenes-dioicae, causing anther-smut disease on the wide-spread plant species Silene latifolia and S. dioica, respectively. Using model-based clustering algorithms on microsatellite data from samples across Europe, we identified rare disease transmission between the host species and rare pathogen hybrids. Using a coalescent-based approach and an isolation-with-migration model, the age of divergence between the two fungal species was estimated at approximately 4.2 × 10(5) years. Levels of gene flow were low and concentrated in very recent times. In addition, gene flow appeared unidirectional from M. silenes-dioicae to M. lychnidis-dioicae. Altogether, our findings are consistent with a scenario of recurrent introgressive hybridization but at a very low level and through secondary contact following initial divergence in allopatry. Asymmetry in the direction of gene flow mirrors previous findings on introgression between the two host plants. Our study highlights the consequences of bringing closely related pathogens into contact, which is increasing through modern global changes and favors cross-species disease transmission, hybridization, and introgression by pathogens.

Evolutionary and functional properties of a two-locus beta-globin polymorphism in Indian house mice

Electrophoretic surveys of hemoglobin (Hb) polymorphism in house mice from South Asia and the Middle East have revealed that two alternative beta-globin haplotypes, Hbb(d) and Hbb(p), are often present at intermediate frequencies in geographically disparate populations. Both haplotypes harbor two functionally distinct beta-globin paralogs, HBB-T1 (which encodes the beta-chain subunits of the major Hb isoform) and HBB-T2 (which encodes the beta-chains of the minor Hb isoform). The Hbb(d) and Hbb(p) haplotypes share identical HBB-T1 alleles, but products of the alternative HBB-T2 alleles (d(minor) and p(minor)) are distinguished by two amino acid substitutions. To investigate the possible adaptive significance of the Hbb(d)/Hbb(p) polymorphism we conducted a population genetic analysis of the duplicated beta-globin genes of Indian house mice (Mus castaneus) in conjunction with experimental studies of Hb function in inbred strains of mice that carry the alternative Hbb(d) and Hbb(p) haplotypes. The main objectives of this study were (i) to characterize patterns of nucleotide polymorphism and linkage disequilibrium in the duplicated beta-globin genes of M. castaneus, (ii) to test the hypothesis that the Hbb(d) and Hbb(p) haplotypes are maintained as a balanced polymorphism, and (iii) to assess whether allelic differences in the alternative minor Hb isoforms (d(minor) and p(minor)) are associated with different O(2)-binding properties. A multilocus analysis of polymorphism and divergence revealed that levels of diversity at the HBB-T2 gene exceeded neutral expectations, and reconstructed haplotype networks for both beta-globin paralogs revealed extensive allele sharing with several other closely related species of Mus. However, despite this suggestive evidence for balancing selection, O(2)-equilibrium curves revealed no discernible functional differences between red cell lysates containing the d(minor) and p(minor) Hb isoforms. If the d(minor) and p(minor) alleles are maintained as a balanced polymorphism, our results indicate that the associated fitness variance is not directly related to respiratory functions of Hb.

Islands of speciation or mirages in the desert? Examining the role of restricted recombination in maintaining species

Over the past decade, many studies documented high genetic divergence between closely related species in genomic regions experiencing restricted recombination in hybrids, such as within chromosomal rearrangements or areas adjacent to centromeres. Such regions have been called 'islands of speciation' because of their presumed role in maintaining the integrity of species despite gene flow elsewhere in the genome. Here, we review alternative explanations for such patterns. Segregation of ancestral variation or artifacts of nucleotide diversity within species can readily lead to higher F(ST) in regions of restricted recombination than other parts of the genome, even in the complete absence of interspecies gene flow, and thereby cause investigators to erroneously conclude that islands of speciation exist. We conclude by discussing strengths and weaknesses of various means for testing the role of restricted recombination in maintaining species.

Evidence for pervasive adaptive protein evolution in wild mice

The relative contributions of neutral and adaptive substitutions to molecular evolution has been one of the most controversial issues in evolutionary biology for more than 40 years. The analysis of within-species nucleotide polymorphism and between-species divergence data supports a widespread role for adaptive protein evolution in certain taxa. For example, estimates of the proportion of adaptive amino acid substitutions (alpha) are 50% or more in enteric bacteria and Drosophila. In contrast, recent estimates of alpha for hominids have been at most 13%. Here, we estimate alpha for protein sequences of murid rodents based on nucleotide polymorphism data from multiple genes in a population of the house mouse subspecies Mus musculus castaneus, which inhabits the ancestral range of the Mus species complex and nucleotide divergence between M. m. castaneus and M. famulus or the rat. We estimate that 57% of amino acid substitutions in murids have been driven by positive selection. Hominids, therefore, are exceptional in having low apparent levels of adaptive protein evolution. The high frequency of adaptive amino acid substitutions in wild mice is consistent with their large effective population size, leading to effective natural selection at the molecular level. Effective natural selection also manifests itself as a paucity of effectively neutral nonsynonymous mutations in M. m. castaneus compared to humans.

Fine-scale phylogenetic discordance across the house mouse genome

Population genetic theory predicts discordance in the true phylogeny of different genomic regions when studying recently diverged species. Despite this expectation, genome-wide discordance in young species groups has rarely been statistically quantified. The house mouse subspecies group provides a model system for examining phylogenetic discordance. House mouse subspecies are recently derived, suggesting that even if there has been a simple tree-like population history, gene trees could disagree with the population history due to incomplete lineage sorting. Subspecies of house mice also hybridize in nature, raising the possibility that recent introgression might lead to additional phylogenetic discordance. Single-locus approaches have revealed support for conflicting topologies, resulting in a subspecies tree often summarized as a polytomy. To analyze phylogenetic histories on a genomic scale, we applied a recently developed method, Bayesian concordance analysis, to dense SNP data from three closely related subspecies of house mice: Mus musculus musculus, M. m. castaneus, and M. m. domesticus. We documented substantial variation in phylogenetic history across the genome. Although each of the three possible topologies was strongly supported by a large number of loci, there was statistical evidence for a primary phylogenetic history in which M. m. musculus and M. m. castaneus are sister subspecies. These results underscore the importance of measuring phylogenetic discordance in other recently diverged groups using methods such as Bayesian concordance analysis, which are designed for this purpose.

Molecular insights into the colonization and chromosomal diversification of Madeiran house mice

The colonization history of Madeiran house mice was investigated by analysing the complete mitochondrial (mt) D-loop sequences of 156 mice from the island of Madeira and mainland Portugal, extending on previous studies. The numbers of mtDNA haplotypes from Madeira and mainland Portugal were substantially increased (17 and 14 new haplotypes respectively), and phylogenetic analysis confirmed the previously reported link between the Madeiran archipelago and northern Europe. Sequence analysis revealed the presence of four mtDNA lineages in mainland Portugal, of which one was particularly common and widespread (termed the 'Portugal Main Clade'). There was no support for population bottlenecks during the formation of the six Robertsonian chromosome races on the island of Madeira, and D-loop sequence variation was not found to be structured according to karyotype. The colonization time of the Madeiran archipelago by Mus musculus domesticus was approached using two molecular dating methods (mismatch distribution and Bayesian skyline plot). Time estimates based on D-loop sequence variation at mainland sites (including previously published data from France and Turkey) were evaluated in the context of the zooarchaeological record of M. m. domesticus. A range of values for mutation rate (mu) and number of mouse generations per year was considered in these analyses because of the uncertainty surrounding these two parameters. The colonization of Portugal and Madeira by house mice is discussed in the context of the best-supported parameter values. In keeping with recent studies, our results suggest that mutation rate estimates based on interspecific divergence lead to gross overestimates concerning the timing of recent within-species events.

Rapid evolution of mouse Y centromere repeat DNA belies recent sequence stability

The Y centromere sequence of house mouse, Mus musculus, remains unknown despite our otherwise significant knowledge of the genome sequence of this important mammalian model organism. Here, we report the complete molecular characterization of the C57BL/6J chromosome Y centromere, which comprises a highly diverged minor satellite-like sequence (designated Ymin) with higher-order repeat (HOR) sequence organization previously undescribed at mouse centromeres. The Ymin array is approximately 90 kb in length and resides within a single BAC clone that provides sequence information spanning an endogenous animal centromere for the first time. By exploiting direct patrilineal inheritance of the Y chromosome, we demonstrate stability of the Y centromere DNA structure spanning at least 175 inbred generations to beyond the time of domestication of the East Asian M.m. molossinus "fancy" mouse through which the Y chromosome was first introduced into the classical inbred laboratory mouse strains. Despite this stability, at least three unequal genetic exchange events have altered Ymin HOR unit length and sequence structure since divergence of the ancestral Mus musculus subspecies around 900,000 yr ago, with major turnover of the HOR arrays driving rapid divergence of sequence and higher-order structure at the mouse Y centromere. A comparative sequence analysis between the human and chimpanzee centromeres indicates a similar rapid divergence of the primate Y centromere. Our data point to a unique DNA sequence and organizational architecture for the mouse Y centromere that has evolved independently of all other mouse centromeres.

Evolution of duplicated beta-globin genes and the structural basis of hemoglobin isoform differentiation in Mus

The functional diversification of multigene families may be strongly influenced by mechanisms of concerted evolution such as interparalog gene conversion. The beta-globin gene family of house mice (genus Mus) represents an especially promising system for evaluating the effects of gene conversion on the functional divergence of duplicated genes. Whereas the majority of mammalian species possess tandemly duplicated copies of the adult beta-globin gene that are identical in sequence, natural populations of house mice are often polymorphic for distinct two-locus haplotypes that differ in levels of functional divergence between duplicated beta-globin genes, HBB-T1 and HBB-T2. Here, we use a phylogenetic approach to unravel the complex evolutionary history of the HBB-T1 and HBB-T2 paralogs in a taxonomically diverse set of species in the genus Mus. The main objectives of this study were 1) to reconstruct the evolutionary history of the different HBB haplotypes of house mice, 2) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in promoting concerted evolution, 3) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in creating chimeric genes, and 4) to assess the structural basis of hemoglobin isoform differentiation in species that possess distinct HBB paralogs. Results of our phylogenetic survey revealed that the HBB-T1 and HBB-T2 genes in different species of Mus exhibit the full range of evolutionary outcomes with respect to levels of interparalog divergence. At one end of the spectrum, the two identical HBB paralogs on the Hbb(s) haplotype (shared by Mus domesticus, Mus musculus, and Mus spretus) represent a classic example of concerted evolution. At the other end of the spectrum, the two distinct HBB paralogs on the Hbb(d), Hbb(p), Hbb(w1), and Hbb(w2) haplotypes (shared by multiple species in the subgenus Mus) show no trace of gene conversion and are distinguished by a number of functionally important amino acid substitutions. Because the possession of distinct HBB paralogs expands the repertoire of functionally distinct hemoglobin isoforms that can be synthesized during fetal development and postnatal life, variation in the level of functional divergence between HBB-T1 and HBB-T2 may underlie important physiological variation within and among species.