Glyoxalase 1 (Glo1) has been implicated in anxiety-like behavior in mice and in multiple psychiatric diseases in humans. We used mouse Affymetrix exon arrays to detect copy number variants (CNV) among inbred mouse strains and thereby identified a approximately 475 kb tandem duplication on chromosome 17 that includes Glo1 (30,174,390-30,651,226 Mb; mouse genome build 36). We developed a PCR-based strategy and used it to detect this duplication in 23 of 71 inbred strains tested, and in various outbred and wild-caught mice. Presence of the duplication is associated with a cis-acting expression QTL for Glo1 (LOD>30) in BXD recombinant inbred strains. However, evidence for an eQTL for Glo1 was not obtained when we analyzed single SNPs or 3-SNP haplotypes in a panel of 27 inbred strains. We conclude that association analysis in the inbred strain panel failed to detect an eQTL because the duplication was present on multiple highly divergent haplotypes. Furthermore, we suggest that non-allelic homologous recombination has led to multiple reversions to the non-duplicated state among inbred strains. We show associations between multiple duplication-containing haplotypes, Glo1 expression and anxiety-like behavior in both inbred strain panels and outbred CD-1 mice. Our findings provide a molecular basis for differential expression of Glo1 and further implicate Glo1 in anxiety-like behavior. More broadly, these results identify problems with commonly employed tests for association in inbred strains when CNVs are present. Finally, these data provide an example of biologically significant phenotypic variability in model organisms that can be attributed to CNVs.

Patterns of genetic differentiation among taxa at early stages of divergence provide an opportunity to make inferences about the history of speciation. Here, we conduct a survey of DNA-sequence polymorphism and divergence at loci on the autosomes, X chromosome, Y chromosome and mitochondrial DNA in samples of Mus domesticus, M. musculus and M. castaneus. We analyzed our data under a divergence with gene flow model and estimate that the effective population size of M. castaneus is 200,000-400,000, of M. domesticus is 100,000-200,000 and of M. musculus is 60,000-120,000. These data also suggest that these species started to diverge approximately 500,000 years ago. Consistent with this recent divergence, we observed considerable variation in the genealogical patterns among loci. For some loci, all alleles within each species formed a monophyletic group, while at other loci, species were intermingled on the phylogeny of alleles. This intermingling probably reflects both incomplete lineage sorting and gene flow after divergence. Likelihood ratio tests rejected a strict allopatric model with no gene flow in comparisons between each pair of species. Gene flow was asymmetric: no gene flow was detected into M. domesticus, while significant gene flow was detected into both M. castaneus and M. musculus. Finally, most of the gene flow occurred at autosomal loci, resulting in a significantly higher ratio of fixed differences to polymorphisms at the X and Y chromosomes relative to autosomes in some comparisons, or just the X chromosome in others, emphasizing the important role of the sex chromosomes in general and the X chromosome in particular in speciation.

We have studied different subspecies of the house mouse and their reciprocal F(1) hybrids to estimate the within-locus mode of inheritance for subspecies differences in gene expression in three tissues (brain, liver, and testis) of male mice. This study investigates the mode of inheritance in crosses at a larger taxonomic distance than has been previously systematically investigated. We found the vast majority of transcripts to be additively expressed with only a few transcripts showing dominance or overdominance in expression, except for one direction of one cross, which showed large mis-expression in the testis. We suggest that, as time passes, more genes come to influence expression, and if there is no directional dominance, additivity becomes increasingly more likely, up to a threshold beyond which there is F(1) hybrid breakdown. Some previous studies on different organisms have found a large degree of dominance, commonly at shorter taxonomic differences. We surveyed these findings and show that the most consistent association exists between the amount of additivity detected in a study and the expression analysis method (in particular microarray platform), suggesting that at least some of the differences among studies might be methodological. Most studies agree with ours in that within-locus additivity seems to be general mode of inheritance for transcript expression. Differentially expressed transcripts identified in our screen among subspecies of house mice are candidate genes that could be involved in reproductive isolation between these subspecies.

Contrasting patterns of X-linked vs. autosomal diversity may be indicative of the mode of selection operating in natural populations. A number of observations have shown reduced X-linked (or Z-linked) diversity relative to autosomal diversity in various organisms, suggesting a large impact of genetic hitchhiking. However, the relative contribution of other forces such as population bottlenecks, variation in reproductive success of the two sexes, and differential introgression remains unclear. Here, we survey 13 loci, 6 X-linked and 7 autosomal, in natural populations of the house mouse (Mus musculus) subspecies complex. We studied seven populations of three different subspecies, the eastern house mouse M. musculus castaneus, the central house mouse M. m. musculus, and the western house mouse M. m. domesticus, including putatively ancestral and derived populations for each. All populations display lower diversity on the X chromosomes relative to autosomes, and this effect is most pronounced in derived populations. To assess the role of demography, we fit the demographic parameters that gave the highest likelihood of the data using coalescent simulations. We find that the reduction in X-linked diversity is too large to be explained by a simple demographic model in at least two of four derived populations. These observations are also not likely to be explained by differences in reproductive success between males and females. They are consistent with a greater impact of positive selection on the X chromosome, and this is supported by the observation of an elevated K(A) and elevated K(A)/K(S) ratios on the rodent X chromosome. A second contribution may be that the X chromosome less readily introgresses across subspecies boundaries.

Regulatory changes in genes involved in reproduction are thought to be prime targets for divergence during speciation, since they are expected to play an important role in sexual selection and sexual conflict. We used microarray analysis of RNA from different wild populations of house mouse subspecies (including Mus m. musculus, Mus m. domesticus, and Mus m. castaneus) and from the sister species Mus spretus to test this assumption. A comparison of expression divergence in brain, liver/kidney, and testis shows a major difference in the evolutionary dynamics of testis-related genes. While the comparison between species confirms an excess in divergence in testis genes, we find that all comparisons between subspecies yield only a very small number of genes with significantly different expression levels in the testis. These results suggest that the early phase of the speciation process may not be driven by regulatory changes in genes that are potential targets of sexual selection, and that the divergence in these genes is only established during a later phase of the speciation process.

Changes in gene expression are known to occur between closely related species, but it is not yet clear how many of these are due to random fixation of allelic variants or due to adaptive events. In a microarray survey between subspecies of the Mus musculus complex, we identified the mitogen-activated protein-kinase-kinase MKK7 as a candidate for change in gene expression. Quantitative PCR experiments with multiple individuals from each subspecies confirmed a specific and significant up-regulation in the testis of M. m. domesticus. Northern blot analysis shows that this is due to a new transcript that is not found in other tissues, nor in M. m. musculus. A cis-trans test via allele specific expression analysis of the MKK7 gene in F1 hybrids between domesticus and musculus shows that the expression change is mainly caused by a mutation located in cis. Nucleotide diversity was found to be significantly reduced in a window of at least 20 kb around the MKK7 locus in domesticus, indicative of a selective sweep. Because the MKK7 gene is involved in modulating a kinase signalling cascade in a stress response pathway, it seems a plausible target for adaptive differences between subspecies, although the functional role of the new testis-specific transcripts will need to be further studied.

Understanding the genes that contribute to reproductive isolation is essential to understanding speciation, but isolating such genes has proven very difficult. In this study I apply a multilocus test statistic to >10,000 SNP markers assayed in wild-derived inbred strains of house mice to identify genomic regions of elevated differentiation between two subspecies of house mice, Mus musculus musculus and M. m. domesticus. Differentiation was high through approximately 90% of the X chromosome. In addition, eight regions of high differentiation were identified on the autosomes, totaling 7.5% of the autosomal genome. Regions of high differentiation were confirmed by direct sequencing of samples collected from the wild. Some regions of elevated differentiation have an overrepresentation of genes with host-pathogen interactions and olfaction. The most strongly differentiated region on the X has previously been shown to fail to introgress across a hybrid zone between the two subspecies. This survey indicates autosomal regions that should also be examined for differential introgression across the hybrid zone, as containing potential genes causing hybrid unfitness.