Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter‐subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.     

Seeking signatures of reinforcement at the genetic level: a hitchhiking mapping and candidate gene approach in the house mouse

Reinforcement is the process by which prezygotic isolation is strengthened as a response to selection against hybridization. Most empirical support for reinforcement comes from the observation of its possible phenotypic signature: an accentuated degree of prezygotic isolation in the hybrid zone as compared to allopatry. Here, we implemented a novel approach to this question by seeking for the signature of reinforcement at the genetic level. In the house mouse, selection against hybrids and enhanced olfactory‐based assortative mate preferences are observed in a hybrid zone between the two European subspecies Mus musculus musculus and M. m. domesticus, suggesting a possible recent reinforcement event. To test for the genetic signature of reinforcing selection and identify genes involved in sexual isolation, we adopted a hitchhiking mapping approach targeting genomic regions containing candidate genes for assortative mating in mice. We densely scanned these genomic regions in hybrid zone and allopatric samples using a large number of fast evolving microsatellite loci that allow the detection of recent selection events. We found a handful of loci showing the expected pattern of significant reduction in variability in populations close to the hybrid zone, showing assortative odour preference in mate choice experiments as compared to populations further away and displaying no such preference. These loci lie close to genes that we pinpoint as testable candidates for further investigation.     

Do changes in gene expression contribute to sexual isolation and reinforcement in the house mouse?

Expression divergence, rather than sequence divergence, has been shown to be important in speciation, particularly in the early stages of divergence of traits involved in reproductive isolation. In the two European subspecies of house mice, Mus musculus musculus and Mus musculus domesticus, earlier studies have demonstrated olfactory‐based assortative mate preference in populations close to their hybrid zone. It has been suggested that this behaviour evolved following the recent secondary contact between the two taxa (~3,000 years ago) in response to selection against hybridization. To test for a role of changes in gene expression in the observed behavioural shift, we conducted a RNA sequencing experiment on mouse vomeronasal organs. Key candidate genes for pheromone‐based subspecies recognition, the vomeronasal receptors, are expressed in these organs. Overall patterns of gene expression varied significantly between samples from the two subspecies, with a large number of differentially expressed genes between the two taxa. In contrast, only ~200 genes were found repeatedly differentially expressed between populations within M. m. musculus that did or did not display assortative mate preferences (close to or more distant from the hybrid zone, respectively), with an overrepresentation of genes belonging to vomeronasal receptor family 2. These receptors are known to play a key role in recognition of chemical cues that handle information about genetic identity. Interestingly, four of five of these differentially expressed receptors belong to the same phylogenetic cluster, suggesting specialization of a group of closely related receptors in the recognition of odorant signals that may allow subspecies recognition and assortative mating.     

Transgressive segregation in a behavioural trait? Explorative strategies in two house mouse subspecies and their hybrids

Hybrid zones between genetically diverged populations are widespread among animals and plants. Their dynamics usually depend on selection against admixture and dispersal of parental forms in the zone. Although indirect estimates of selection have been the target of many studies, dispersal has been neglected. In this study we carried out open field experiments to test whether males of two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus, differ in their propensity to disperse and in their character of exploration. We tested wild‐caught males and males of two wild‐derived inbred strains. In addition, we examined reciprocal F₁ crosses to test the prediction that these hybrids display intermediate behaviours. We revealed that M. m. musculus males were less hesitant to enter the experimental arena than were M. m. domesticus males, but once inside the arena their movements were more timid. F₁ males differed from both parental strains, with longer latencies to enter the arena, but explored the arena in a similar fashion as the M. m. domesticus males, thus displaying transgressive behavioural phenotypes. These results contribute to our knowledge of behavioural divergence between the mouse subspecies, and add a new facet to the study of speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●●, ●●–●●.     

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.     

The coevolution of male and female genitalia in a mammal: A quantitative genetic insight

Male genitalia are among the most phenotypically diverse morphological traits, and sexual selection is widely accepted as being responsible for their evolutionary divergence. Studies of house mice suggest that the shape of the baculum (penis bone) affects male reproductive fitness and experimentally imposed postmating sexual selection has been shown to drive divergence in baculum shape across generations. Much less is known of the morphology of female genitalia and its coevolution with male genitalia. In light of this, we used a paternal half-sibling design to explore patterns of additive genetic variation and covariation underlying baculum shape and female vaginal tract size in house mice (Mus musculus domesticus). We applied a landmark-based morphometrics approach to measure baculum size and shape in males and the length of the vaginal tract and width of the cervix in females. Our results reveal significant additive genetic variation in house mouse baculum morphology and cervix width, as well as evidence for genetic covariation between male and female genital measures. Our data thereby provide novel insight into the potential for the coevolutionary divergence of male and female genital traits in a mammal.     

Ancient onset of geographical divergence,interpopulation genetic exchange,and natural selection on the Mc1r coat‐colour gene in the house mouse (Mus musculus)

We examined genetic variation in house mice from India and Pakistan, a predominant part of the predicted homeland of this species and also the territory of the subspecies Mus musculus castaneus (CAS), using a nuclear marker for seven tandemly arranged genes (Fanca–Spire2–Tcf25–Mc1r–Def8–Afg3l1–Dbndd1) and compared them with those previously determined for mice from other parts of Eurasia. Construction of a network with the concatenate sequences yielded three distinct clusters representing the three major subspecies groups: CAS, Mus musculus domesticus (DOM) and Mus musculus musculus (MUS). STRUCTURE analysis provided evidence for further subdivision of CAS into two main haplogroups within the Indian subcontinent. Single‐gene networks revealed not only gene‐specific architecture for subgrouping in CAS, but also allelic exchange among subspecies. These results suggest the earlier onset of allopatric divergence in the predicted homeland (the Middle East and Indian subcontinent) and subsequent intermittent admixing via gene flow across the CAS haplogroups and among the three subspecies groups. A comparison of the levels of nucleotide diversity among the gene regions revealed a less divergent state in the chromosome region containing Mc1r and its adjacent genes, indicative of a selective sweep, suggesting the involvement of natural selection in the Mc1r allelic variation. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 778–794.     

DEVELOPMENTAL STABILITY,FITNESS, AND TRAIT SIZE IN LABORATORY HYBRIDS BETWEEN EUROPEAN SUBSPECIES OF THE HOUSE MOUSE

The effects of hybridization on developmental stability and size of tooth characters were investigated in intersubspecific crosses between random-bred wild strains of the house mouse (Mus musculus domesticus and M. m. musculus). Fluctuating asymmetry (FA) and trait size were compared within and between parental, F₁, backcross, and F₂ hybrid groups. The relationship between FA and reproductive fitness within the F₁ hybrids was also studied. The results indicated that both FA and character size levels differed significantly between the two subspecies. The F₁ hybrids and the recombined groups (backcrosses and F₂ hybrids) showed heterosis for both parameters. No significant differences in the FA of fertile and sterile F₁ hybrid individuals were found. Comparison of the FA levels obtained in this study with those found in wild populations from the hybrid zone in Denmark showed that the levels of FA were lower in laboratory-bred samples than in the wild populations. This study provides further evidence that, in hybrids, the developmental processes underlying most of the morphological traits we studied benefit from a heterotic effect, despite the genomic incompatibilities between the two European house mice revealed by previous genetical and parasitological studies.     

The uncharacterized gene 1700093K21Rik and flanking regions are correlated with reproductive isolation in the house mouse,Mus musculus

Reproductive barriers exist between the house mouse subspecies, Mus musculus musculus and M. m. domesticus, members of the Mus musculus species complex, primarily as a result of hybrid male infertility, and a hybrid zone exists where their ranges intersect in Europe. Using single nucleotide polymorphisms (SNPs) diagnostic for the two taxa, the extent of introgression across the genome was previously compared in these hybrid populations. Sixty-nine of 1316 autosomal SNPs exhibited reduced introgression in two hybrid zone transects suggesting maladaptive interactions among certain loci. One of these markers is within a region on chromosome 11 that, in other studies, has been associated with hybrid male sterility of these subspecies. We assessed sequence variation in a 20 Mb region on chromosome 11 flanking this marker, and observed its inclusion within a roughly 150 kb stretch of DNA showing elevated sequence differentiation between the two subspecies. Four genes are associated with this genomic subregion, with two entirely encompassed. One of the two genes, the uncharacterized 1700093K21Rik gene, displays distinguishing features consistent with a potential role in reproductive isolation between these subspecies. Along with its expression specifically within spermatogenic cells, we present various sequence analyses that demonstrate a high rate of molecular evolution of this gene, as well as identify a subspecies amino acid variant resulting in a structural difference. Taken together, the data suggest a role for this gene in reproductive isolation.     

Patterns of morphological evolution in the mandible of the house mouse Mus musculus (Rodentia: Muridae)

The worldwide distributed house mouse, Mus musculus, is subdivided into at least three lineages, Mus musculus musculus, Mus musculus domesticus, and Mus musculus castaneus. The subspecies occur parapatrically in a region considered to be the cradle of the species in Southern Asia (‘central region’), as well as in the rest of the world (‘peripheral region’). The morphological evolution of this species in a phylogeographical context is studied using a landmark‐based approach on mandible morphology of different populations of the three lineages. The morphological variation increases from central to peripheral regions at the population and subspecific levels, confirming a centrifugal sub‐speciation within this species. Furthermore, the outgroup comparison with sister species suggests that M. musculus musculus and populations of all subspecies inhabiting the Iranian plateau have retained a more ancestral mandible morphology, suggesting that this region may represent one of the relevant places of the origin of the species. Mus musculus castaneus, both from central and peripheral regions, is morphologically the most variable and divergent subspecies. Finally, the results obtained in the present study suggest that the independent evolution to commensalism in the three lineages is not accompanied by a convergence detectable on jaw morphology. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 635–647.     

Asymmetric size and shape variation in the Central European transect across the house mouse hybrid zone

We studied asymmetric variation of the mandible in the Central European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. Within introgression classes, defined by the share of diagnostic allozymes, we quantified the directional and fluctuating component of asymmetric variation, as well as skewness and kurtosis of individual asymmetry distributions. Furthermore, in the same manner we re‐analysed asymmetric variation of the ventral side of the skull. According to the quadratic polynomial model, the mandible shape‐fluctuating asymmetry, but not size‐fluctuating asymmetry, was significantly decreased in the centre of the hybrid zone (with a minimum predicted for a hybrid index of 0.41). On the contrary, the skull shape‐fluctuating asymmetry non‐monotonically increased towards the musculus side of the hybrid zone (with a peak predicted for a hybrid index of 0.86). Thus, the impact of hybridization on fluctuating asymmetry is trait‐specific in this portion of the house mouse hybrid zone. The only general feature of asymmetric variation we observed was the shift towards the platykurtosis of asymmetry distributions in the centre of the hybrid zone. Taken together, we suggest genetic variability for right–left asymmetries to be generally increased, but the developmental instability of mandible shape to be decreased, by hybridization. We hypothesize the decrease of developmental instability to be caused by overdominant effects on developmental dynamics rather than by increased heterozygosity. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 13–27.     

There is no heterotic effect upon developmental stability in the ventral side of the skull within the house mouse hybrid zone

Fluctuating asymmetry (FA) of a body is expected to be related to an organism's developmental instability. We studied patterns of FA in the ventral side of the skull along a transect across the central-European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. We found that FA was not significantly different between introgression classes relative to differences between localities within those classes and that the within-class differences were much larger than differences between individual-signed asymmetries within localities. However, if year was added to the same analysis as another factor, FA was not significantly different among localities within the same introgression class. When individual asymmetries were plotted against individual hybrid indices, hybrids appeared more asymmetric than individuals from outside of the zone. Thus contrary to previous studies, we did not find lower FA indicating heterotic effect in hybrids for the traits studied. It is suggested that the impact of hybridization on FA in the ventral side of the mouse skull is negligible or overwhelmed by other factors.     

The meek inherit the earth: less aggressive wild mice are more successful in challenging situations

Numerous studies have shown an association between aggressiveness and several other behavioural traits. For example, more aggressive animals were bold and active explorers tending to form persistent routines whereas less aggressive animals were shy, careful but more flexible. While the former are thought to be more successful under stable conditions the latter should have advantages in more dynamic situations. These differences can apply not only to individuals but also to populations, species or groups of species with important implications to species distributions and speciation rates. Here we utilized the Morris water task (MWT) to investigate how two subspecies, Mus musculus musculus and M. m. domesticus, known to differ in aggressiveness, cope with stressful situations. We found that less aggressive musculus males performed significantly better in solving the MWT than more aggressive domesticus males. This suggests that M. m. musculus is more flexible and could be more successful under stressful and/or dynamic situations typical of dispersal bouts. It seems plausible that this difference may have had an influence on the secondary contact between musculus and domesticus populations in the past and perhaps still can affect the dynamics of the European hybrid zone between the subspecies. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 310–319.     

Speciation and reduced hybrid female fertility in house mice

In mammals, intrinsic postzygotic isolation has been well studied in males but has been less studied in females, despite the fact that female gametogenesis and pregnancy provide arenas for hybrid sterility or inviability that are absent in males. Here, we asked whether inviability or sterility is observed in female hybrids of Mus musculus domesticus and M. m. musculus, taxa which hybridize in nature and for which male sterility has been well characterized. We looked for parent‐of‐origin growth phenotypes by measuring adult body weights in F1 hybrids. We evaluated hybrid female fertility by crossing F1 females to a tester male and comparing multiple reproductive parameters between intrasubspecific controls and intersubspecific hybrids. Hybrid females showed no evidence of parent‐of‐origin overgrowth or undergrowth, providing no evidence for reduced viability. However, hybrid females had smaller litter sizes, reduced embryo survival, fewer ovulations, and fewer small follicles relative to controls. Significant variation in reproductive parameters was seen among different hybrid genotypes, suggesting that hybrid incompatibilities are polymorphic within subspecies. Differences in reproductive phenotypes in reciprocal genotypes were observed and are consistent with cyto‐nuclear incompatibilities or incompatibilities involving genomic imprinting. These findings highlight the potential importance of reduced hybrid female fertility in the early stages of speciation.     

Ultrasonic Vocalizations of Male Mice Differ among Species and Females Show Assortative Preferences for Male Calls

Male house mice (Mus musculus) emit ultrasonic vocalizations (USVs) during courtship, which attract females, and we aimed to test whether females use these vocalizations for species or subspecies recognition of potential mates. We recorded courtship USVs of males from different Mus species, Mus musculus subspecies, and populations (F1 offspring of wild-caught Mus musculus musculus, Mus musculus domesticus (and F1 hybrid crosses), and Mus spicilegus), and we conducted playback experiments to measure female preferences for male USVs. Male vocalizations contained at least seven distinct syllable types, whose frequency of occurrence varied among species, subspecies, and populations. Detailed analyses of multiple common syllable types indicated that Mus musculus and Mus spicilegus could be discriminated based on spectral and temporal characteristics of their vocalizations, and populations of Mus musculus were also distinctive regardless of the classification model used. Females were able to discriminate USVs from different species, and showed assortative preferences for conspecific males. We found no evidence that females discriminate USVs of males from a different subspecies or separate populations of the same species, even though our spectral analyses identified acoustic features that differ between species, subspecies, and populations of the same species. Our results provide the first comparison of USVs between Mus species or between Mus musculus subspecies, and the first evidence that male USVs potentially facilitate species recognition.     

Detection of recombinant haplotypes in wild mice (Mus musculus) provides new insights into the origin of Japanese mice

Japanese house mice (Mus musculus molossinus) are thought to be a hybrid lineage derived from two prehistoric immigrants, the subspecies M. m. musculus of northern Eurasia and M. m. castaneus of South Asia. Mice of the western European subspecies M. m. domesticus have been detected in Japanese ports and airports only. We examined haplotype structuring of a 200 kb stretch on chromosome 8 for 59 mice from throughout Eurasia, determining short segments (≈ 370–600 bp) of eight nuclear genes (Fanca, Spire2, Tcf25, Mc1r, Tubb3, Def8, Afg3l1 and Dbndd1) which are intermittently arranged in this order. Where possible we identified the subspecies origin for individual gene alleles and then designated haplotypes for concatenated alleles. We recovered 11 haplotypes among 19 Japanese mice examined, identified either as ‘intact’ haplotypes derived from the subspecies musculus (57.9%), domesticus (7.9%), and castaneus (2.6%), or as ‘recombinant’ haplotypes (31.6%). We also detected recombinant haplotypes unique to Sakhalin. The complex nature of the recombinant haplotypes suggests ancient introduction of all three subspecies components into the peripheral part of Eurasia or complicated genomic admixture before the movement from source areas. ‘Intact’domesticus and castaneus haplotypes in other Japanese wild mice imply ongoing stowaway introductions. The method has general utility for assessing the history of genetic admixture and for disclosing ongoing genetic contamination.     

Divergence of odorant signals within and between the two European subspecies of the house mouse

The olfactory acuity of mice allows them to discriminate odorsof conspecifics differing by a few genes. This acuity is usedin habituation procedures where investigation of novel odorsby the mouse can be translated into relative difference or similaritybetween the stimuli. This study adapts these behavioral proceduresto address suprapopulation divergence among urinary odors inthe house mouse. Specifically, we investigate geographical patternsof odor divergence within and between 2 subspecies of the housemouse, Mus musculus musculus and Mus musculus domesticus, whichdiverged in allopatry and met secondarily in Europe where theyhybridize. Based on M. m. musculus perception, our study suggeststhat odors of the 2 subspecies differ in both allopatric andcontact zone populations and that divergence is more markedin the latter. Our earlier studies documented mate preferenceand signal divergence between the 2 subspecies. Hence, we considerthe role of the urinary odors as mating signals. We discusshow signal divergence between the 2 subspecies may relate toreproductive character displacement. This study validates theuse of habituation procedures to reconstruct geographical patternsof odorant signal divergence, providing a strong methodologicalplatform to address reproductive character displacement affectingcryptic mating signals in mammals.