Asymmetrical reproductive character displacement in the house mouse

Our study addressed reproductive character displacement between two subspecies of the house mouse, Mus musculus musculus and Mus musculus domesticus that hybridize in Europe along a zone where selection against hybridization is known to occur. Based on a multi-population approach, we investigated spatial patterns of divergence of mate preference in the two taxa. Mate preference was significantly higher in the contact zone than in allopatry in both subspecies, suggesting that reproductive character displacement occurs. Moreover, patterns of preference were stronger in M. m. musculus than in M. m. domesticus, indicating an asymmetrical divergence between the two. In the context of selection against hybridization, our results may provide empirical support for the hypothesis of reinforcement in a parapatric hybrid zone. We discuss factors that could explain the asymmetrical pattern of divergence and the possible impact of a unimodal structure on the maintenance of premating divergence between the two subspecies.     

Intrinsic reproductive isolating mechanisms in the maintenance of a hybrid zone between ecologically divergent subspecies

Understanding factors involved in maintaining stable hybrid zones is important for predicting the ultimate fate of the interacting taxa, but the relative importance of mechanisms such as ecological selection and intrinsic reproductive isolation remains unclear. Most studies of reproductive isolation in hybrid zones have focused either on zones with strongly bimodal patterns in genotype or phenotype frequencies, with relatively strong isolation, or unimodal zones with relatively weak isolation, whereas less is known about more intermediate classes of hybrid zone. Here, we utilize a hybrid zone of this intermediate type occurring between northern and southern subspecies of Atlantic killifish, Fundulus heteroclitus, to identify isolating mechanisms playing a role in maintaining this type of zone. The two subspecies differ in environmental tolerance, and we found some evidence of microhabitat preference between subspecies within a small tidal creek at the centre of the hybrid zone. There was also an association between sex, mitochondrial genotype and habitat within this creek. Fertilization success did not differ between consubspecific and heterosubspecific crosses, but hatching success was significantly lower for crosses involving southern males and northern females, and crosses between southern females and northern males had altered developmental rates. Southern females and northern males showed patterns consistent with positive assortative mating. Together, these results indicate a role for a combination of factors including assortative mating and/or early hybrid inviability in the maintenance of this hybrid zone and suggest that hybrid zones with intermediate levels of reproductive isolation are likely to be maintained by multiple interacting isolating mechanisms.     

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.     

Development of unique house mouse resources suitable for evolutionary studies of speciation

Two house mouse subspecies, Mus musculus domesticus and Mus musculus musculus, form a hybrid zone in Europe and represent a suitable model for inferring the genes contributing to isolation barriers between parental taxa. Despite long-term intensive studies of this hybrid zone, we still know relatively little about the causes and mechanisms maintaining the 2 taxa as separate subspecies; therefore, to gain insight into this process, we developed 8 wild-derived inbred house mouse strains. In order to produce strains as pure domesticus or musculus genomes as possible, the individuals used to establish the breeding colony for the 3 domesticus and 2 of the musculus strains were captured in the Czech Republic from wild populations at extreme western and eastern edges of the subspecific contact zone, respectively. The remaining 3 musculus strains were bred from mice captured about 250 km east of the hybrid zone. Genetic analysis based on 361 microsatellite loci showed that 82% of these markers are diagnostic for either the musculus or the domesticus strains. In order to demonstrate the potential utility of this genetic differentiation in such strains, phenotypic variation was scored for 2 strains from opposite edges of the hybrid zone and significant differences in morphology, reproductive performance, in vitro immune responses, mate choice based on urinary signals, and aggressiveness were found. In addition, the 3 strains derived from musculus populations far from the hybrid zone display significant differences in polymorphism in hybrid male sterility when crossed with the laboratory strains C57BL/6 or C57BL/10, which have a predominantly domesticus genome. Although further studies will be necessary to demonstrate intersubspecific differences, all analyses presented here indicate that these newly developed house mouse strains represent a powerful tool for elucidating the genetic basis of isolation barriers in hybrid zones and for studying speciation in general.     

The role of salivary androgen-binding protein in reproductive isolation between two subspecies of house mouse: Mus musculus musculus and Mus musculus domesticus

Previous behavioural studies using inbred lines have suggested that the gene ( Abpa ) for the alpha subunit of salivary androgen-binding protein (ABP) plays a role in prezygotic isolation between house mouse Mus musculus subspecies. We tested this hypothesis in animals from wild allopatric (121 individuals from four samples) and parapatric (320 animals from 15 samples) populations sampled on the Czech–Bavarian transect across the hybrid zone between M. m. domesticus and M. m. musculus . The study did not reveal a consistent statistically significant trend of homosubspecific preferences in individual allopatric and parapatric populations. Nonetheless, the whole pattern of preference was skewed toward homosubspecific preference mostly on the M. m. musculus side of the hybrid zone. The pattern of homosubspecific preferences was stronger for the time spent sniffing than it was for the first choice of the signal (the ratio of homosubspecific vs. heterosubspecific preferences for both sexes was 6 : 2 in allopatric and 21 : 9 in parapatric populations, while the same rates were 4 : 4 and 16 : 14 for the first choice). To the extent that Y-maze tests reflect preference under wild conditions, we suggest that this slight preference may not in itself be sufficient to impede gene flow between the two subspecies and thus act as a reproductive barrier. ABP most probably participates in a complex system of subspecies-specific recognition in the hybrid zone, but the picture is far too complex at this time to allow a conclusive evaluation of the importance of this role.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 349–361.     

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.     

GENOMIC AND MORPHOLOGICAL ANALYSIS OF A SEMIPERMEABLE AVIAN HYBRID ZONE SUGGESTS ASYMMETRICAL INTROGRESSION OF A SEXUAL SIGNAL

Hybrid zones are geographic regions where differentiated taxa meet and potentially exchange genes. Increasingly, genomic analyses have demonstrated that many hybrid zones are semipermeable boundaries across which introgression is highly variable. In some cases, certain alleles penetrate across the hybrid zone in only one direction, recombining into the alternate genome. We investigated this phenomenon using genomic (genotyping‐by‐sequencing) and morphological (plumage reflectance spectrophotometry) analyses of the hybrid zone between two subspecies of the red‐backed fairy‐wren (Malurus melanocephalus) that differ conspicuously in a sexual signal, male back plumage color. Geographic cline analyses revealed a highly variable pattern of differential introgression, with many narrow coincident clines combined with several significantly wider clines, suggesting that the hybrid zone is a semipermeable tension zone. The plumage cline was shifted significantly into the genomic background of the orange subspecies, consistent with sexual selection driving asymmetrical introgression of red plumage alleles across the hybrid zone. This interpretation is supported by previous experimental work demonstrating an extra‐pair mating advantage for red males, but the role of genetic dominance in driving this pattern remains unclear. This study highlights the potential for sexual selection to erode taxonomic boundaries and promote gene flow, particularly at an intermediate stage of divergence.     

Environmental factors, microgeographic patterns of endemism and hybrid zones in Monobella grassei (Insecta: Collembola: Neanuridae)

A complex pattern of narrow parapatric distributions with hybrid zones is described among subspecies of the Collembola species Monobella grassei in the eastern Pyrenees. Though hybrid zones are usually narrow, they vary from simple contact zones without recorded hybrids to relatively wide hybrid swarms, suggesting that large differences exist locally in the level of genetic divergence between subspecies. The course and structure of contact zones are analysed on morphological grounds, in relation to environmental (rainfall, vegetation, substrate) and topographical (altitude, hydrography) parameters. The results show that, unlike most cases documented so far, the courses of the contact zones in M. grassei are not connected to major ecological boundaries nor to ecological ecotones, on either the regional or the habitat scale. It is suggested that this is due to the great ecological similarity between the different parapatric subspecies of M. grassei , and their broad ecological requirements. There is no conclusive evidence that the distribution observed is linked to recent palaeoclimatic events such as glaciation. Endogenous selection could therefore be prevalent in this model to explain parapatry and the hybrid zone pattern.     

Hybrids between Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) do not show endogenous postzygotic isolation

Studies of prezygotic and postzygotic isolation in Drosophila have shown in general that species in sympatry tend to evolve prezygotic barriers earlier than do species in allopatry. However, postzygotic barriers tend to evolve at the same evolutionary rate in both sympatric and allopatric species. In contrast to these observations, the grasshoppers Chorthippus parallelus parallelus and C. p. erythropus show complete hybrid male sterility but only limited prezygotic isolation after an estimated 0.5 millions years of divergence. Like their congeners, C. brunneus and C. jacobsi form a hybrid zone where their ranges meet in northern Spain. However, the hybrid zone is mosaic and bimodal and, in contrast to the high levels of postzygotic isolation between C. parallelus subspecies, these two species showed no significant reduction in hybrid fitness in F₁ or backcross generations relative to the parental generations. The level of prezygotic isolation in laboratory tests was comparable to that between C. parallelus subspecies. These results suggest that endogenous postzygotic isolation does not play an important role in the reproductive isolation between C. brunneus and C. jacobsi , or in determining the structure of the hybrid zone. Exogenous postzygotic isolation may be present and should be tested in future studies.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 195–203.     

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography-dependent asymmetrical dispersal or selection favouring one subspecies over the other.     

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.     

Testing for selection on the androgen-binding protein in the Danish mouse hybrid zone

The three peripheral subspecies of the house mouse Mus musculus have fixed specific variants of the androgen-binding protein (ABP) that have been proposed to be part of a recognition mechanism that could participate in sexual isolation between the subspecies. We tested for selection on Abpa by characterizing the pattern of Abpa introgression across a transect of the hybrid zone between M. m. musculus and M. m. domesticus in Jutland. On the musculus side, the cline for Abpa resembled that of a nearly neutral allozyme more than that of strongly selected X and Y chromosome markers. However, the high central step which displaces the tail of introgression of Abpa to higher frequencies was best accounted for by linkage to a locus under strong selection against hybrids. Still, we cannot exclude that this pattern results from weak selection against Abpa in the tail of introgression, which would be compatible with an assortative choice mechanism. On the domesticus side there was little introgression close to the hybrid zone, presumably due to a geographical barrier to migration. However, substantial frequencies of musculus alleles occurred further away, suggesting mixed colonization patterns as well as flow across the hybrid zone.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 447–459.     

Strikingly variable divergence times inferred across an Amazonian butterfly 'suture zone'

'Suture zones' are areas where hybrid and contact zones of multiple taxa are clustered. Such zones have been regarded as strong evidence for allopatric divergence by proponents of the Pleistocene forest refugia theory, a vicariance hypothesis frequently used to explain diversification in the Amazon basin. A central prediction of the refugia and other vicariance theories is that the taxa should have a common history so that divergence times should be coincident among taxa. A suture zone for Ithomiinae butterflies near Tarapoto, NE Peru, was therefore studied to examine divergence times of taxa in contact across the zone. We sequenced 1619bp of the mitochondrial COI/COII region in 172 individuals of 31 species from across the suture zone. Inferred divergence times differed remarkably, with divergence between some pairs of widespread species (each of which may have two or more subspecies interacting in the zone, as in the genus Melinaea) being considerably less than that between hybridizing subspecies in other genera (for instance in Oleria). Our data therefore strongly refute a simple hypothesis of simultaneous vicariance and suggest that ongoing parapatric or other modes of differentiation in continuous forest may be important in driving diversification in Amazonia.     

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.     

Geographic variation in the structure of oak hybrid zones provides insights into the dynamics of speciation

Studying geographic variation in the rate of hybridization between closely related species could provide a useful window on the evolution of reproductive isolation. Reinforcement theory predicts greater prezygotic isolation in areas of prolonged contact between recently diverged species than in areas of recent contact, which implies that old contact zones would be dominated by parental phenotypes with few hybrids (bimodal hybrid zones), whereas recent contact zones would be characterized by hybrid swarms (unimodal hybrid zones). Here, we investigate how the hybrid zones of two closely related Chinese oaks, Quercus mongolica and Q. liaotungensis, are structured geographically using both nuclear and chloroplast markers. We found that populations of Q. liaotungensis located around the Changbai Mountains in Northeast China, an inferred glacial refugium, were introgressed by genes from Q. mongolica, suggesting historical contact between the two species in this region. However, these introgressed populations form sharp bimodal hybrid zones with Q. mongolica. In contrast, populations of Q. liaotungensis located in North China, which show no sign of ancient introgression with Q. mongolica, form unimodal hybrid zones with Q. mongolica. These results are consistent with the hypothesis that selection against hybrids has had sufficient time to reinforce the reproductive barriers between Q. liaotungensis and Q. mongolica in Northeast China but not in North China.     

Behavioural response to song and genetic divergence in two subspecies of white‐crowned sparrows (Zonotrichia leucophrys)

Divergence in sexual signals may drive reproductive isolation between lineages, but behavioural barriers can weaken in contact zones. Here, we investigate the role of song as a behavioural and genetic barrier in a contact zone between two subspecies of white‐crowned sparrows (Zonotrichia leucophrys). We employed a reduced genomic data set to assess population structure and infer the history underlying divergence, gene flow and hybridization. We also measured divergence in song and tested behavioural responses to song using playback experiments within and outside the contact zone. We found that the subspecies form distinct genetic clusters, and demographic inference supported a model of secondary contact. Song phenotype, particularly length of the first note (a whistle), was a significant predictor of genetic subspecies identity and genetic distance along the hybrid zone, suggesting a close link between song and genetic divergence in this system. Individuals from both parental and admixed localities responded significantly more strongly to their own song than to the other subspecies song, supporting song as a behavioural barrier. Putative parental and admixed individuals were not significantly different in their strength of discrimination between own and other songs; however, individuals from admixed localities tended to discriminate less strongly, and this difference in discrimination strength was explained by song dissimilarity as well as genetic distance. Therefore, we find that song acts as a reproductive isolating mechanism that is potentially weakening in a contact zone between the subspecies. Our findings also support the hypothesis that intraspecific song variation can reduce gene flow between populations.     

Using the Mus musculus hybrid zone to assess covariation and genetic architecture of limb bone lengths

Two subspecies of the house mouse, Mus musculus domesticus and Mus musculus musculus, meet in a narrow contact zone across Europe. Mice in the hybrid zone are highly admixed, representing the full range of mixed ancestry from the two subspecies. Given the distinct morphologies of these subspecies, these natural hybrids can be used for genomewide association mapping at sufficiently high resolution to directly infer candidate genes. We focus here on limb bone length differences, which is of special interest for understanding the evolution of developmentally correlated traits. We used 172 first‐generation descendants of wild‐caught mice from the hybrid zone to measure the length of stylopod (humerus/femur), zeugopod (ulna/tibia) and autopod (metacarpal/metatarsal) elements in skeletal CT scans. We find phenotypic covariation between limb elements in the hybrids similar to patterns previously described in Mus musculus domesticus inbred strains, suggesting that the hybrid genotypes do not influence the covariation pattern in a major way. Mapping was performed using 143,592 SNPs and identified several genomic regions associated with length differences in each bone. Bone length was found to be highly polygenic. None of the candidate regions include the canonical genes known to control embryonic limb development. Instead, we are able to identify candidate genes with known roles in osteoblast differentiation and bone structure determination, as well as recently evolved genes of, as yet, unknown function.     

Reduced male fertility is common but highly variable in form and severity in a natural house mouse hybrid zone

Barriers to gene flow between naturally hybridizing taxa reveal the initial stages of speciation. Reduced hybrid fertility is a common feature of reproductive barriers separating recently diverged species. In house mice (Mus musculus), hybrid male sterility has been studied extensively using experimental crosses between subspecies. Here, we present the first detailed picture of hybrid male fertility in the European M. m. domesticus-M. m. musculus hybrid zone. Complete sterility appears rare or absent in natural hybrids but a large proportion of males (~30%) have sperm count or relative testis weight below the range in pure subspecies, and likely suffer reduced fertility. Comparison of a suite of traits related to fertility among subfertile males indicates reduced hybrid fertility in the contact zone is highly variable among individuals and ancestry groups in the type, number, and severity of spermatogenesis defects present. Taken together, these results suggest multiple underlying genetic incompatibilities are segregating in the hybrid zone, which likely contribute to reproductive isolation between subspecies.     

Phenotypic effects of the Y chromosome are variable and structured in hybrids among house mouse recombinant lines

Hybrid zones between divergent populations sieve genomes into blocks that introgress across the zone, and blocks that do not, depending on selection between interacting genes. Consistent with Haldane's rule, the Y chromosome has been considered counterselected and hence not to introgress across the European house mouse hybrid zone. However, recent studies detected massive invasion of M. m. musculus Y chromosomes into M. m. domesticus territory. To understand mechanisms facilitating Y spread, we created 31 recombinant lines from eight wild‐derived strains representing four localities within the two mouse subspecies. These lines were reciprocally crossed and resulting F1 hybrid males scored for five phenotypic traits associated with male fitness. Molecular analyses of 51 Y‐linked SNPs attributed ~50% of genetic variation to differences between the subspecies and 8% to differentiation within both taxa. A striking proportion, 21% (frequencies of sperm head abnormalities) and 42% (frequencies of sperm tail dissociations), of phenotypic variation was explained by geographic Y chromosome variants. Our crossing design allowed this explanatory power to be examined across a hierarchical scale from subspecific to local intrastrain effects. We found that divergence and variation were expressed diversely in different phenotypic traits and varied across the whole hierarchical scale. This finding adds another dimension of complexity to studies of Y introgression not only across the house mouse hybrid zone but potentially also in other contact zones.