Skip to main content
  • Research article
  • Open access
  • Published:

Measurements of hybrid fertility and a test of mate preference for two house mouse races with massive chromosomal divergence

Abstract

Background

Western house mice Mus musculus domesticus are among the most important mammalian model species for chromosomal speciation. Hybrids between chromosomal races of M. m. domesticus suffer various degrees of fertility reduction between full fertility and complete sterility, depending on the complexity of the chromosomal differences between the races. This complexity presents itself in hybrids as meiotic configurations of chromosome chains and rings, with longer configurations having a stronger impact on fertility. While hybrids with short configurations have been intensively studied, less work has been done on hybrids with very long configurations. In this study, we investigated laboratory-reared wild mice from two chromosomally very different races in Switzerland found in close proximity. Hybrids between these races form a meiotic chain of fifteen chromosomes. We performed a detailed analysis of male and female hybrid fertility, including three generations of female backcrosses to one of the parental races. We also tested for possible divergence of mate preference in females.

Results

While all male F1 hybrids were sterile with sperm counts of zero, 48% of female F1 hybrids produced offspring. Their litter sizes ranged from one to three which is significantly lower than the litter size of parental race females. When hybrid females were backcrossed to a parental race, half of the offspring resembled the parental race in karyotype and fertility, while the other half resembled the F1 hybrids. In the preference test, females of both races indicated a lack of a preference for males of their own karyotype.

Conclusions

Although the fertility of the F1 hybrids was extremely low because of the complexity of the chromosomal differences between the races, reproductive isolation was not complete. As we did not find assortative female preferences, we expect that contact between these races would lead to the production of hybrids and that gene flow would occur eventually, as fertility can be restored fully after one backcross generation.

Background

Speciation is a central topic in evolutionary biology [1], with a need for detailed studies of incipient reproductive isolation of major genetic forms within species through an analysis of hybridisation either in the laboratory or in nature, i.e. where the races make contact at hybrid zones [2]. One characteristic that can define major genetic forms within species is the presence of chromosomal rearrangements; such forms may be termed chromosomal races [3]. Chromosomal rearrangements have the potential to be potent agents promoting reproductive isolation. Mechanical problems during meiosis in hybrids between chromosomal races may lead to subfertility or sterility (hybrid dysfunction) [4, 5]. These include incorrect pairing of parental chromosomes, malsegregation of chromosomes, and errors during crossing-over leading to deletion or duplication of chromosomal regions. In addition, chromosomal rearrangements can also suppress recombination and protect larger parts of the genome from introgression despite hybridisation, facilitating the accumulation of species-specific gene variants involved in Dobzhansky-Muller incompatibilities in these regions [6, 7]. Such incompatibilities may also be related to mate choice.

Chromosomal rearrangements between races and species of mammals include inversions, translocations and fusions/fissions. The western house mouse (Mus musculus domesticus) is an excellent model for the study of a very common rearrangement in mammalian karyotype evolution, the Robertsonian (Rb) fusion (= centric fusions): the joining of a pair of acrocentric chromosomes at their centromeres to form a metacentric [8]. While the ancestral karyotype of the western house mouse consists of 2n = 40 acrocentrics, Rb fusions have led to metacentric races—over 100 have been described—with reduced diploid numbers down to 2n = 22 [9, 10]. Hybrids between these races have been shown to suffer a decrease in fitness linked to the complexity of their chromosomal differences [10]. Similar observations have been made for hybrids between chromosomal races in other mammalian species [3, 11], although the house mouse is exceptional for the ease of laboratory studies and availability of genomic and other resources [12].

The individual metacentric races in house mice occur over very limited geographical areas, although groups of races may occur near to one another, such as in eastern Switzerland (Fig. 1), from where the metacentric races in the present study derive [9]. Two types of heterozygote can be formed on hybridisation of metacentric races in nature or in captivity: simple and complex [3]. Simple heterozygotes are formed when one race has a specific metacentric (e.g. 1.3, formed from chromosomes 1 and 3) and the other one has these chromosomes in the ancestral acrocentric form (1 and 3 in this case). Hybrids are heterozygous for this metacentric, which leads to the formation, during prophase I of meiosis, of a chain of three chromosomes—a trivalent (here 1–1.3–3)—instead of bivalents formed from homologues of either the 1.3 metacentrics or the acrocentrics 1 and 3. Complex heterozygotes are formed when hybridising races have metacentrics that share arms (monobrachial homology, e.g. 1.3 and 3.6). Hybrids in this case produce chains or rings involving four or more chromosomes during meiosis. These configurations are called multivalents (e.g. 1–1.3–3.6–6, a chain of four chromosomes).

Fig. 1
figure 1

Map of chromosomal variation in house mice in Switzerland and neighbouring areas. The trapping sites of the two races used in this study are marked with stars (pink for CHHN, red for CHBU: see text). The red star represents three closely located trapping sites (within 4 km). Information on the Swiss races was taken from Piálek et al. [9], Hübner [27] and Gropp et al. [64]

Simple heterozygotes with few trivalents (those with one to three have been well studied) show near-normal fertility, while those with a higher number of trivalents (up to nine) are more severely affected and show decreasing fertility with increasing number of trivalents [13, 14]. In contrast, complex heterozygotes usually suffer a more substantial decrease in fertility, with chain meiotic configurations being more detrimental than rings [15, 16]. Although the reduction in fertility is strongly linked to the complexity of the meiotic configurations, there also seems to be an effect of the genetic background [13, 17, 18]. Also, female heterozygotes are generally less severely affected than males [10, 15]. The most detailed studies of fertility of hybrids between chromosomal races, with substantial quantitative data, relate to simple heterozygotes and complex heterozygotes for short chain configurations. There are data for both males and females, but substantially more on males [13, 14, 17, 19,20,21,22,23,24].

Here we carry out detailed fertility studies on both male and female complex heterozygous F1 hybrids between two chromosomal races that produce a very long chain of fifteen chromosomes at meiosis. We test the degree of F1 hybrid sterility by performing sperm counts on males and backcrossing females. We repeated this procedure for the offspring of the backcrosses. As F1 hybrids are expected to produce only two types of balanced gametes (Fig. 2), F1 females should have offspring with two karyotypes: half should match the F1 karyotype, the other half that of the parental race used in the backcross.

Fig. 2
figure 2

Expected inheritance of a meiotic chain of fifteen chromosomes considering parental, F1 and backcross individuals. F1 female mice are expected to produce gametes containing either chain chromosomes from their CHBU parent (white), or chain chromosomes from their CHHN parent (black). The fusion of these oocytes with sperm from CHHN males (black) should lead to backcross offspring (BC) with the hybrid karyotype (2n = 23) or the CHHN karyotype (2n = 24). Chromosomes not involved in the chain need not segregate according to race of origin

In addition to the effect of chromosomal differences, geographic separation itself could lead to divergence of behaviour, which could act as a premating barrier if the races were to come into contact. Such premating divergence has been described between the western house mouse and its sister subspecies, the eastern house mouse Mus musculus musculus (e.g. [25, 26]). Therefore, we tested for female preference when given a choice between males of the two races.

Methods

Mice

Mice from two very different and previously little studied chromosomal races from Switzerland were used in this study (Fig. 1). Following the chromosome and race nomenclature in Piálek et al. [9], they are the Buchs race CHBU with 2n = 22 (1.18 2.5 3.6 4.12 7.15 8.16 9.14 10.17 11.13 [27]), and the Hünikon race CHHN with 2n = 24 (1.3 2.8 4.12 5.7 6.15 9.14 10.11 13.16 [27]). F1 hybrids between these races are expected to have the following configurations during meiosis I: a multivalent chain of fifteen chromosomes (18–18.1–1.3–3.6–6.15–15.7–7.5–5.2–2.8–8.16–16.13–13.11–11.10–10.17–17) and three autosomal bivalents formed from homologues of 4.12, 9.14 and 19, respectively, as well as the sex-bivalent (Fig. 3). All animals used in this study were laboratory-reared and are descended from wild mice (CHHN and CHBU: four and five generations in the laboratory, respectively). CHHN mice were trapped in Illnau, Switzerland (see König & Lindholm [28] for a description of the study population), while CHBU mice were trapped in farms and stables in the nearby locations of Grabs, Buchs and Haag, Switzerland, between 2010 and 2014 (trapping locations shown in Fig. 1). The approximate direct distance between Illnau and Buchs is 60 km. Previous sampling in this region found CHHN and CHBU races at sites as close as 30 km apart, but with no sampling in between [27].

Fig. 3
figure 3

Diakinesis preparation from testes of a male F1 hybrid between the races CHHN and CHBU. The chain of fifteen chromosomes is visible in the centre, with two sets of metacentric bivalents above and below (filled arrows), the chromosome 19 bivalent (dashed arrow) and the sex bivalent (arrowhead)

Mice were kept either singly in Makrolon Type II cages or in groups in Makrolon Type III cages (Indulab). Groups consisted of same-sex siblings or of a breeding pair. Cages contained bedding, cardboard housing and paper towels as nesting material. Food and water was available ad libitum. Animals in the fertility experiment were kept under a light:dark cycle of 14:10 with lights on at 05.30 CET. Animals in the mate preference test were kept under a reversed cycle with lights on at 17.30 CET.

Crosses

The first generation (P) consisted of 14 crosses between mice from CHBU and CHHN (five ♀CHBUx♂CHHN, nine ♀CHHNx♂CHBU). One to nine litters per breeding pair were obtained. We assessed the fertility of hybrid offspring. Where available, two male offspring per breeding pair were used for sperm counts. On average two female offspring (2.2 ± 0.7; mean ± SD) per breeding pair were backcrossed to unrelated CHHN males. This process was repeated for the following generations, with males sacrificed for sperm counts and females backcrossed to CHHN males (Fig. 4, Table 1). Females in backcrosses were kept with males for up to five months (12–159 d) in order to also assess litter size of females with very low fecundity. However, if breeding pairs showed high levels of aggression within those five months, they were separated prematurely but not excluded from the analysis. Thirty-two (74%) backcrossed females were separated prematurely (after 76 ± 37 d; details in Additional file 1). In addition, if females had two litters with litter sizes of at least five pups each, we considered them as fully fertile and terminated the cross. We performed sperm counts on two male offspring from such females but did not backcross any of their female offspring. Mice were checked for pups on a weekly basis.

Fig. 4
figure 4

Crossing scheme applied between mice that derive from the Buchs (CHBU) and Hünikon (CHHN) races. Male offspring were sacrificed and sperm counts were performed. Female offspring were backcrossed to a male of the CHHN race. After the first backcross, not all offspring were processed equally (stars). Fewer offspring were analysed if their mother was fully fertile (see text)

Table 1 Details of animal numbers in all experimental crosses (see Fig. 4)

Karyotype analysis

Mitotic chromosome spreads were prepared according to Ford [29] and stained with Giemsa. For each individual, at least twenty clear and complete metaphases were counted. According to our expectations outlined in the Introduction (Fig. 2), all F1 mice were expected to have 2n = 23 chromosomes. We verified this by counting chromosomes of seven randomly chosen F1. After backcrossing to a male of the parental CHHN race with 2n = 24, we expected offspring of either 2n = 23 (hybrid karyotype) or 2n = 24 (parental CHHN karyotype). Thus, we made chromosome counts for offspring of all female 2n = 23 hybrid mice (this included offspring from BC1, BC2 and BC3 crosses). Chromosomes from offspring from 2n = 24 hybrid females (this included offspring from BC2 and BC3 crosses) were not counted.

Offspring were classified into three categories: (1) F1 hybrids with 2n = 23, (2) BC.23 hybrids with 2n = 23, resulting from backcrosses to CHHN, and (3) BC.24 hybrids with 2n = 24 also resulting from backcrosses to CHHN.

Fertility of males

Male offspring were sacrificed at an age of nine weeks (mean ± SD, 65 ± 3 d). Twenty-five male F1 hybrids were analysed. As controls, 12 male CHBU and 17 male CHHN of similar age were used (62 ± 3 d, 64 ± 5 d, respectively). The right cauda epididymis was processed as described in Hauffe & Searle [22]. Spermatozoa were counted in both chambers of a Neubauer Improved haemocytometer. Sperm estimates were calculated by dividing the number of sperm by the volume of the haemocytometer squares counted (ten squares of 0.04 mm2 with 0.1 mm depth or, if no sperm was visible in that area, the whole 9 mm2) and then by multiplying with the total solution used to dilute the epididymis (2 ml). The results from both chambers were averaged for each male. The resulting value estimates the total number of sperm in one cauda epididymis. Other data collected from males were as follows: body mass after euthanasia, epididymis mass (mean of separate measurements of right and left complete epididymides), testis mass (mean of separate measurements of right and left testes), mass of seminal vesicle and coagulating gland (as these organs could not be reliably separated, they were weighed together; half the mass of combined left and right dissected as one entity) and preputial gland mass (half the mass of combined left and right dissected as one entity). Male offspring of backcrosses were processed in the same fashion (N = 21) and included male offspring from fully fertile BC2 and BC3 females (two per female). To increase sample size for these comparisons, we set up additional crosses in exactly the same way as the P crosses to provide 35 additional males (mean age 60 ± 3 d; 17 F1, 12 CHBU, 6 CHHN).

We applied linear mixed-effects models (LMM) with parent pair as a random effect in order to take into account that siblings were included in the analysis. Independent variables were category of male (CHBU, CHHN, F1, BC.23, BC.24) and body mass. We used the lme4, lsmeans and lmerTest packages in R 3.3.0 [30,31,32,33]. Denominator degrees of freedom were approximated using the Satterthwaite method. We tested the assumption of normality for all models in this study with a Normal Q-Q plot of the residuals. However, the assumption of normality was not well satisfied for sperm number. Thus we first applied a Fisher’s exact test to analyse the difference in prevalence of males with and without sperm depending on whether a meiotic chain of fifteen was expected or not (F1 and BC.23 vs CHBU, CHHN and BC.24). Then we applied the LMM described above on the male categories with sperm (CHBU, CHHN and BC.24) after performing a common logarithmic transformation on sperm number.

In order to demonstrate the presence of a meiotic chain configuration of fifteen chromosomes in F1 hybrids, we made diakinesis preparations from the testes of two male F1 hybrids according to Bulatova et al. [34].

Fertility of females

Female fertility was estimated by whether or not pups were born, and by the number of offspring in the female’s first litter (31 F1, 12 BC). For a subset of these females (19 F1, 12 BC), the number of implantation scars in the uterus was counted post mortem from fresh tissue [35, 36]. In order to attribute any decrease in fertility to the females rather than to the males, we performed sperm counts on the latter. As estimates for number of offspring in females of the parental races, we used data from the same breeding system of our laboratory that had furnished P mice and CHHN backcross males for this experiment. From all first litters born to females from within-race breeding pairs over the period 2013–2016, we randomly selected 12 litters from CHBUxCHBU crosses, and 12 litters from CHHNxCHHN crosses. We used a one-way ANOVA to compare litter sizes between experimental females and controls. Additionally, we performed a LMM on the data of the experimental females with litter size as the dependent variable, category of female and body mass as independent variables, and parent pair as a random effect to correct for the presence of full-siblings. The same R packages were used as for the analysis of male fertility. Plots were created using ggplot2 [37].

Female preference test

Female mate preference was tested in an apparatus of 75x55cm containing three chambers: two small ones for the males (each 25.5 × 19.5 cm) and a large one for the female (Fig. 5) [38,39,40]. The male chambers were separated by an opaque wall which extended into the female’s chamber. Male and female chambers were separated by a mesh (0.5 cm2 grids) that allowed vocal, visual and tactile interactions between the male and the female but prevented copulation. Males were unable to see or access each other. Fifty-eight mature but sexually inexperienced female mice (two to four months old) from the races CHBU (N = 30) and CHHN (N = 28) were tested for mate preference using pairs of CHBU and CHHN males, also sexually inexperienced. Males were allocated randomly to the left and right male chambers. Twenty-two of thirty-six (61%) male pairs were used twice, once with a CHBU female and once with a CHHN female. All experiments were performed between 11.00 and 15.00 CET in a room with reversed light cycle (light:dark cycle of 14:10, lights on at 17.30 CET; habituation period of two weeks) and filmed from above under red light in the dark without the presence of an observer. Males were first placed into their respective chambers during a habituation period of 15 min. Afterwards, a female was placed into the neutral zone (see Fig. 5) and allowed to roam freely within her chamber. We only used females at the pre-oestrus and oestrus stage of their cycle, verified by vaginal smears as described in Byers et al. [41]. After the experiment, the movements of female and male mice were scored blindly with respect to their race using the event recorder BORIS [42]. Measurements were started after the female had exited the neutral zone on both male sides once [43, 44], and lasted 30 min [39, 40]. We measured the time females spent with their nose within 1 cm of the mesh [39, 40]. In addition, we measured the time males spent with their nose within 1 cm of the mesh.

Fig. 5
figure 5

Diagram of the preference test apparatus. The thick dashed black line indicates the mesh grid. The dashed grey line demarks the neutral zone and does not represent a physical barrier

For the statistical analysis, the proportion of time P a female spent near the mesh of the CHHN male was calculated as (time at mesh of CHHN male)/(time at mesh of CHHN male + time at mesh of CHBU male). LMMs were performed following logit transformation of the response variable time proportion P. Transformations were performed using the boot package in R [45, 46]. Fixed effects were female race, difference in male body mass, difference in male age, and difference in time males spent near the mesh. All differences between males were calculated as (CHHN value) − (CHBU value). These differences were centred, and scaled to standardise variation. Random effects were ID of the male pair used and the chamber side the CHHN male was placed in (right or left, to correct for a potential preference of the female for a specific side).

Results

Karyotype analysis

The seven analysed F1 hybrids had a diploid chromosome number of 2n = 23 with 17 metacentrics and 6 acrocentrics. Offspring from F1 females had either 2n = 23 with 17 metacentrics and 6 acrocentrics (N = 4), or 2n = 24 with 16 metacentrics and 8 acrocentrics (N = 11). This was also the case for offspring from females with 2n = 23 from the BC1 and BC2 generation (offspring with 2n = 23: N = 4, offspring with 2n = 24: N = 3, Table 1).

Fertility of males

All F1 and BC.23 males had sperm counts of zero, which was significantly different from the males of the remaining categories (Fisher’s exact test, p < 10− 15; Table 2; Fig. 6a). For the males of these remaining categories, which all had sperm counts well above zero, the linear mixed effects model did not indicate a significant effect of male category but did indicate a significant effect of body mass (Tables 3 and 4). In the remaining models containing all five male categories, category of male was highly significant for testis mass and epididymis mass (Table 2; Fig. 6b). Category of male was marginally significant for the mass of the seminal vesicle with coagulating gland but not for the preputial gland mass (Table 3). Testis mass and epididymis mass were significantly lower for F1 and BC.23 males (Table 4). However, the mass of the seminal vesicle with coagulating gland and the preputial gland mass were not significantly lower for F1 or BC.23 males (Table 4). There was a significant positive relationship between body mass and testis mass, epididymis mass, mass of the seminal vesicle with coagulating gland, and preputial gland mass (Table 3). The relationship between body mass and the response variables was examined in detail using linear regression models, which indicated that for sperm number, testis mass and epididymis mass, the positive relationship was limited to sperm producing males (CHHN, CHBU, BC.24; Additional file 2).

Table 2 Mean and standard deviation for male fertility measurements for all groups
Fig. 6
figure 6

Boxplot representation of measurements taken from male mice for the comparison of fertility. a Number of sperm in the cauda epididymis and (b) the mass of epididymis, testis, seminal vesicle with coagulating gland, and preputial gland were compared between males of the CHBU race, the CHHN race, F1 hybrids between CHBU and CHHN, and offspring resulting from backcrosses of hybrids to CHHN (see Fig. 4). Offspring from backcrosses were grouped according to their diploid chromosome number (BC.23 for 2n = 23, BC.24 for 2n = 24). The boxplot boundaries indicate the minimum, the first quantile, the median, the third quantile and the maximum. Outliers are shown as circles. The sample sizes are listed below or above the corresponding boxplots. Note that in the male fertility analysis, the category BC.24 includes offspring from fully fertile BC2 and BC3 females (two per female)

Table 3 Results of ANOVA analyses of effect of male category and body mass on male fertility
Table 4 Pairwise comparisons of reproductive parameters for different categories of males

Fertility of females

Female hybrids were backcrossed to males of the CHHN chromosomal population over three generations (crosses BC1–BC3, Fig. 4, animal numbers in Table 1). They gave birth to their first litters between 3 and 20 weeks after pairing (mean ± SD, 8.0 ± 5.0 weeks). From the 31 F1 females which were all backcrossed to CHHN males, 15 (48%) of them produced offspring (Table 1). From the ten BC1 females and the two BC2 females backcrossed, all produced offspring. In the LMM for size of first litter, female category was significant (F4, 39.7 = 63.61, p < 10-15). F1 and BC.23 females both had significantly smaller first litters than females of the categories CHHN, CHBU and BC.24 but did not differ significantly from one another (CHHN: 6.3 ± 1.5, CHBU: 6.1 ± 1.4, F1: 1.1 ± 0.3, BC.23: 1.8 ± 1.0, BC.24: 7.5 ± 0.9; Table 5, Fig. 7).

Table 5 Pairwise comparisons of size of first litter for different categories of females using Dunnett’s method
Fig. 7
figure 7

Number of offspring in each litter for females that reproduced. For the experimental crosses, all litters born are shown. For the controls, only the first litter is shown. The boxplots are overlain with the individual data points

In the additional comparisons between experimental females, BC.24 females reproduced sooner than F1 and BC.23 females (mean ± SD time to first litter: BC.24: 3.5 ± 0.8 weeks, F1: 8.0 ± 5.0 weeks, BC.23: 8.5 ± 2.1 weeks; Wilcoxon-Mann-Whitney test, BC.24–F1: W = 99, p = 0.012; BC.24–BC.23: W = 32, p = 0.008). We also examined 31 females for the presence of scars in the uterus as an indication of successful embryo implantation. Generally, we expect fertile females to have the same number of scars as offspring if all implanted embryos survived. In contrast, an excess of scars would indicate that not all of the implanted embryos survived. Twenty-two (71%) of the females examined had more scars than offspring (Fig. 8). While the number of scars in BC.24 females was very similar to the number of total offspring, the number of scars in many F1 and BC.23 females greatly exceeded the number of offspring, with the number of scars ranging from zero to thirty-one (Fig. 8). Three females without offspring had multiple scars. Four females had one scar less than offspring number, suggesting that scar numbers may have been underestimated.

Fig. 8
figure 8

Total offspring number versus number of uterine scars for the experimental females (F1’s and backcrosses). Females are separated into three categories: F1 between the races CHBU and CHHN (N = 19), females from backcrosses with 2n = 23 (BC.23, N = 4), and females from backcrosses with 2n = 24 (BC.24, N = 8). The dashed line indicates where the number of scars equals the number of total offspring i.e. every embryo that was implanted survived until birth. Females above the dashed line had more scars than offspring, indicating the loss of implanted embryos. The comparison of cross duration for the different females shows that it was not a major explanatory factor for differences in scar and offspring number

The CHHN males used in the backcrosses had a mean ± SD number of sperm in the cauda epididymis of 23.1 ± 12.1 million, which is in the range of the sperm counts for CHHN control males in the male hybrid fertility analysis (Fig. 6).

Female preference test

Neither CHHN nor CHBU female mice spent significantly more than 50% of their time at the mesh with males of any race (CHHN females: t = − 0.35, p = 0.77; CHBU females: t = 1.65, p = 0.30). However, CHBU females spent a greater proportion of their time close to the CHHN male than did CHHN females (Table 6, Fig. 9). Difference of male age influenced time near the CHHN male: the older the CHHN male was relative to the CHBU male, the less time females spent near him. There was no effect of time males spent near the mesh, nor of differences in body mass (Table 6). Male age did not correlate with body mass (linear regression, F1, 114 = 0.30, p = 0.58). On average, females spent 34.5 ± 11.9% of the experimental time near the meshes.

Table 6 Linear mixed model of logit transformed relative time proportion spent by female near CHHN male
Fig. 9
figure 9

Relative time that females spent near the CHHN male when presented with males of both races. A value of 0.5 (dashed line) represents equal amount of time spent near either male. Positive values indicate more time spent near the CHHN male, negative values more time spent near the CHBU male. The grey circles represent the actual data points, while the bars represent the linear mixed model’s estimate of the means (thick horizontal bars) and the 95%-confidence intervals for females of either race

Discussion

Studying hybridisation between major genetic forms within species is one approach to better understand the process of diversification. In this study, we investigated the effect of chromosomal rearrangements, specifically Robertsonian fusions, on the fertility of hybrids between two such genetic forms—chromosomal races—of house mice. F1 hybrids (2n = 23) between the races CHBU (2n = 22) and CHHN (2n = 24), which show substantial chromosomal differences, form a multivalent chain of fifteen chromosomes during meiosis. We report in a detailed analysis that male F1 hybrids are sterile while female F1’s have low fertility. Backcrossing of female F1 to males of the CHHN race resulted in the restoration of the CHHN karyotype and full fertility in half of the offspring. However, we consider hybrid fertility overall to remain relatively low as the occurrence of fertile offspring would be strongly limited by the very low F1 fertility, and the probability of a sterile backcross offspring. In the preference test, females did not show a preference for males of either race.

Hybrid fertility

Spermatogenesis was disrupted in the male F1 hybrids, as they had sperm counts of zero with a lower mass of the testis and epididymis but normally sized seminal vesicle (measured with coagulating gland) and preputial gland. Low sperm counts linked to low testes mass have been previously described, e.g. in Flachs et al. [47]. This tight relationship between sperm count and mass of the testis and epididymis is expected for physiological reasons, as sperm is located in these organs and should thus influence their mass.

Unlike the males, F1 female hybrids were partially fertile. Half of the F1 females produced litters upon backcrossing. However, it took them longer than for parental race females until their first litter was born and the number of offspring was significantly decreased. Most hybrid females had litter sizes of only one offspring instead of the usual six to seven observed in mice of the parental races. F1 females also had an excess of uterine scars that indicated that they had lost multiple foetuses. The cause is most likely aneuploidy. Aneuploid embryos and foetuses have been previously found in pregnant mice heterozygous for metacentrics [48]. In addition, extremely high rates of nondisjunction, ranging up to 100%, have been described in female complex heterozygotes [15, 22]. Nondisjunction results in monosomic and trisomic zygotes after fertilisation. While monosomics die before implantation, trisomics usually survive past implantation and leave visible scars in the uterus [49, 50]. Thus, the scars in the F1 females could indicate the loss of trisomic foetuses.

The fertility results of these F1 hybrids concur with previous studies on complex heterozygotes with long chain configurations during meiosis [15, 16, 18, 51, 52]. In general, male-limited sterility with disrupted spermatogenesis has been observed. However, these studies are mostly qualitative, some without detailed information on sample sizes or the races that were crossed. In addition, the fertility of females was not always investigated. We aimed to contribute a detailed and complete study of complex heterozygote fertility. The most similar study to ours is that of Capanna et al. [15], in which fertility of F1 hybrids between the Poschiavo race (2n = 26) and the Cittaducale race (2n = 22) were studied. These complex heterozygotes possess a chain of seventeen and a chain of three chromosomes during meiosis. Preparation of spermatocytes from nine male hybrids and histological inspection of the testes of eight males indicated that the formation of mature sperm was absent. In contrast, 13 oocyte preparations of meiosis II could be collected from eight female hybrids. Of the ten female hybrids that were crossed to fertile males during four to five months, only two litters with one offspring each were obtained. Both of these offspring were sterile, although no information on karyotype was given. These results are similar to our findings, although female F1’s between the CHBU und CHHN races appear to have a higher fertility. This could be due to them having only one, slightly shorter, meiotic chain with fifteen chromosomes or to a more similar genetic background, as these races are geographically more closely located than the Poschiavo and Cittaducale races.

Thus in Capanna et al. [15] and the current study, gametogenesis in complex heterozygotes with very long chain configurations was affected differently in males and females. This accords with Haldane’s Rule [53], and matches findings in mice that mutations affecting meiotic chromosome alignment and segregation generally have stronger effects in sperm than in oocytes, with cessation of spermatogenesis and continued oogenesis [54]. Garagna et al. [19] studied oogenesis directly and found that in heterozygous metacentric mice, egg follicles that successfully pass the first stages of oogenesis also finish their development. In addition, they described a reduction in follicle number for complex heterozygote females with long meiotic chains, indicating a shorter reproductive lifespan. We thus expect CHBUxCHHN F1 females not only to suffer from reduced litters and an increased number of aneuploid offspring due to nondisjunction, but to also suffer a decrease in reproductive lifespan due to a smaller supply of follicles.

We expected that backcrossing female F1 (2n = 23) to CHHN males (2n = 24) would lead to two types of offspring (Fig. 2): those with a hybrid karyotype (2n = 23) and those with the CHHN karyotype (2n = 24). This pattern was corroborated in our findings. In addition, the fertility analyses showed that when the CHHN karyotype was restored in the offspring, their fertility was normal and equal to that of mice of the parental races. In contrast, offspring with the F1 hybrid karyotype displayed the fertility of the F1 hybrids with sterility in males and low fertility in females. Under the assumption of no recombination (a very simplified scenario), we would expect backcross offspring with 2n = 24 chromosomes to not only regain the parental karyotype, but also the genetic material from the chain chromosomes of that race, while offspring with 2n = 23 would be genetically heterozygous for all chain chromosomes (Fig. 2). However, there is strong support that recombination is obligate in mice and other mammals [55] and studies on simple heterozygotes between mouse chromosomal races have shown that a minimum of one recombination event is maintained [56,57,58]. While complex heterozygotes with long chains have not yet been studied, it would seem probable that recombination rates remain close to the seemingly obligate one recombination event per chromosome arm. Thus, under the assumption of recombination, offspring from backcrosses should be genetic hybrids, even if they have the karyotype of a parental race. As no intermediate fertility values were observed, karyotype and fertility restoration seem to be directly related, suggesting that the decrease in fertility in hybrids between CHBU and CHHN is entirely due to chromosomal heterozygosity instead of genic incompatibilities or other genic effects. Even closely related individuals such as full siblings showed opposite fertility patterns if their karyotypes differed. These results differ from findings of studies in which the fertility of hybrids between mouse chromosomal races varied and indicated a possible involvement of a genic effect [13, 17, 22]. Nevertheless, as chromosomal races in the house mouse all belong to the same subspecies, the accumulation of Dobzhansky-Muller incompatibilities seems unlikely.

Female preference

We did not find a preference for males of the same chromosomal race in our behavioural test, in which differences in the duration of proximity of the female to each male was used as a proxy for mate choice, a common experimental design [59], and in which differences in male behaviour, age and body mass were accounted for in the analysis. These results suggest that females would not discriminate between males of compatible versus incompatible karyotype in their mate choice. In a system of chromosomal races on Madeira, assortative preferences were also absent, although the races did show a divergence of their preference [60]. Furthermore, behavioural differences between neighbouring metacentric races have been previously described [61,62,63]. In our study, we found a trend for an effect of male behaviour and a significant effect of age on the females’ behaviour. Thus, females may have been responding to other cues, e.g. of male quality or interest. In addition, we cannot exclude that laboratory rearing had an effect on the expression of female preferences or male traits.

As there does not seem to be assortative mating preferences between the CHBU and CHHN races, we expect contact to lead to the formation of F1 hybrids and, as female F1’s are not sterile, to allow gene flow between the races, potentially opposing the evolution of assortative mating. As fertility is further restored in part of the backcross offspring, we expect gene flow to be higher than expected when only considering F1 fertility.

Conclusions

This study presents a detailed report on fertility of laboratory-reared hybrids between two chromosomal races of house mice with a high complexity of chromosomal differences. Despite the long meiotic multivalent of fifteen chromosomes in F1 hybrids and sterility of male F1’s, female F1’s remained fertile, although to a low degree. Backcrosses of these females to males of one of the parental races resulted for part of the offspring in restoration of the karyotype to that of the parental race and full fertility, indicating that the decrease in fertility of hybrids is likely entirely due to karyotype. The remaining offspring had the karyotype and low fertility of F1 hybrids, with males being sterile and females having low fertility. A preference test on females of the parental races indicated a lack of preference for chromosomally compatible males. Thus, these races are expected to hybridise in the event of contact. While the chromosomal differences have the potential to promote reproductive isolation between the races through the extremely low fertility of F1 hybrids, the marginal fertility of female F1’s combined with the restoration of fertility in part of the offspring from backcrosses could lead to more gene flow than expected based solely on the length of the chain configuration in hybrids.

References

  1. Coyne JA, Orr HA. Speciation. Sunderland, MA: Sinauer; 2004.

    Google Scholar 

  2. Barton NH, Hewitt GM. Analysis of hybrid zones. Annu Rev Ecol Syst. 1985;16:113–48.

    Article  Google Scholar 

  3. Searle JB. Chromosomal hybrid zones in eutherian mammals. In: Harrison RG, editor. Hybrid zones and the evolutionary process. New York: Oxford University Press; 1993. p. 309–53.

    Google Scholar 

  4. White MJD. Modes of speciation. San Francisco: Freeman; 1978.

  5. King M. Species evolution: the role of chromosome change. New York: Cambridge University Press; 1993.

    Google Scholar 

  6. Rieseberg LH. Chromosomal rearrangements and speciation. Trends Ecol Evol. 2001;16:351–8.

    Article  Google Scholar 

  7. Fishman L, Stathos A, Beardsley P, Williams C, Hill J. Chromosomal rearrangements and the genetics of reproductive barriers in Mimulus (monkey flowers). Evolution. 2013;67:2547–60.

    Article  Google Scholar 

  8. Qumsiyeh MB. Evolution of number and morphology of mammalian chromosomes. J Hered. 1994;85:455–65.

    Article  CAS  Google Scholar 

  9. Piálek J, Hauffe HC, Searle JB. Chromosomal variation in the house mouse. Biol J Linn Soc. 2005;84:535–63.

    Article  Google Scholar 

  10. Hauffe HC, Giménez MD, Searle JB. Chromosomal hybrid zones in the house mouse. In: Macholán M, Baird SJE, Munclinger P, Piálek J, editors. Evolution of the house mouse. Cambridge: Cambridge University Press; 2012. p. 407–530.

    Chapter  Google Scholar 

  11. Close RL, Bell JN, Dollin AE, Harding HR. Spermatogenesis and synaptonemal complexes of hybrid Petrogale (Marsupialia). J Hered. 1996;87:96–107.

    Article  CAS  Google Scholar 

  12. Phifer-Rixey M, Nachman MW. Insights into mammalian biology from the wild house mouse Mus musculus. elife. 2015;4:e05959.

    Article  Google Scholar 

  13. Chatti N, Britton-Davidian J, Catalan J, Auffray J-C, Saïd K. Reproductive trait divergence and hybrid fertility patterns between chromosomal races of the house mouse in Tunisia: analysis of wild and laboratory-bred males and females. Biol J Linn Soc. 2005;84:407–16.

    Article  Google Scholar 

  14. Castiglia R, Capanna E. Contact zone between chromosomal races of Mus musculus domesticus. 2. Fertility and segregation in laboratory-reared and wild mice heterozygous for multiple Robertsonian rearrangements. Heredity. 2000;85:147–56.

    Article  Google Scholar 

  15. Capanna E, Gropp A, Winking H, Noack G, Civitelli M-V. Robertsonian metacentrics in the mouse. Chromosoma. 1976;58:341–53.

    Article  CAS  Google Scholar 

  16. Malorni W, Capanna E, Cristaldi M, De Martino C. Changes of seminiferous epithelium in hybrids of mice carrying Robertsonian karyotype. Arch Androl. 1982;9:333–41.

    Article  CAS  Google Scholar 

  17. Saïd K, Sâad A, Auffray J-C, Britton-Davidian J. Fertility estimates in the Tunisian all-acrocentric and Robertsonian populations of the house mouse and their chromosomal hybrids. Heredity. 1993;71:532–8.

    Article  Google Scholar 

  18. Gropp A, Winking H. Robertsonian translocations: cytology, meiosis, segregation patterns and biological consequences of heterozygosity. Symp Zool Soc Lond. 1981;47:141–81.

    Google Scholar 

  19. Garagna S, Redi C, Zuccotti M, Britton-Davidian J, Winking H. Kinetics of oogenesis in mice heterozygous for Robertsonian translocations. Differentiation. 1990;42:167–71.

    Article  CAS  Google Scholar 

  20. Merico V, Pigozzi MI, Esposito A, Merani MS, Garagna S. Meiotic recombination and spermatogenic impairment in Mus musculus domesticus carrying multiple simple Robertsonian translocations. Cytogenet Genome Res. 2003;103:321–9.

    Article  CAS  Google Scholar 

  21. Nunes AC, Catalan J, Lopez J, Ramalhinho MG, Mathias ML, Britton Davidian J. Fertility assessment in hybrids between monobrachially homologous Rb races of the house mouse from the island of Madeira: implications for modes of chromosomal evolution. Heredity. 2011;106:348–56.

    Article  CAS  Google Scholar 

  22. Hauffe HC, Searle JB. Chromosomal heterozygosity and fertility in house mice (Mus musculus domesticus) from northern Italy. Genetics. 1998;150:1143–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Redi CA, Garagna S, Hilscher B, Winking H. The effects of some Robertsonian chromosome combinations on the seminiferous epithelium of the mouse. J Embryol Exp Morphol. 1985;85:1–19.

    CAS  PubMed  Google Scholar 

  24. Britton-Davidian J, Sonjaya H, Catalan J, Cattaneo-Berrebi G. Robertsonian heterozygosity in wild mice: fertility and transmission rates in Rb(16.17) translocation heterozygotes. Genetica. 1990;80:171–4.

    Article  CAS  Google Scholar 

  25. Smadja C, Catalan J, Ganem G. Strong premating divergence in a unimodal hybrid zone between two subspecies of the house mouse. J Evol Biol. 2004;17:165–76.

    Article  CAS  Google Scholar 

  26. Bímová B, Albrecht T, Macholán M, Piálek J. Signalling components of the house mouse mate recognition system. Behav Process. 2009;80:20–7.

    Article  Google Scholar 

  27. Hübner, RKP. Chromosomal and biochemical variation in wild mice from Switzerland: relevance for models of chromosomal evolution in European house mice. PhD thesis. Oxford: University of Oxford; 1992.

  28. König B, Lindholm AK. The complex social environment of female house mice (Mus domesticus). In: Macholán M, Baird SJE, Munclinger P, Pialek J, editors. Evolution of the house mouse. Cambridge: Cambridge University Press; 2012. p. 114–34.

    Chapter  Google Scholar 

  29. Ford CE. The use of chromosome markers. In: Micklem HS, Loutit JF, editors. Tissue grafting and radiation. New York: Academic Press; 1966. p. 197–206.

    Google Scholar 

  30. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016. https://www.R-project.org/

    Google Scholar 

  31. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

    Article  Google Scholar 

  32. Lenth RV. Least-squares means: the R package lsmeans. J Stat Softw. 2016;69:1–33.

    Article  Google Scholar 

  33. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest: tests in linear mixed effects models. R package version. 2015;2:–33 https://CRAN.R-project.org/package=lmerTest.

  34. Bulatova NS, Searle JB, Nadjafova RS, Pavlova SV, Bystrakova NV. Field protocols for the genomic era. Comp Cytogenet. 2009;3:57–62.

    Article  Google Scholar 

  35. Adamczewska-Andrzejewska KA. Estimations of young born on the basis of placental scars count in the laboratory mouse. Acta Theriol. 1969;14:263–71.

    Article  Google Scholar 

  36. Krackow S. Sex-specific embryonic mortality during concurrent pregnancy and lactation in house mice. J Exp Zool. 1990;256:106–12.

    Article  CAS  Google Scholar 

  37. Wickham H. ggplot2: Elegant graphics for data analysis. New York: Springer; 2009. http://ggplot2.org

  38. Ryan KK, Altmann J. Selection for male choice based primarily on mate compatibility in the oldfield mouse, Peromyscus polionotus rhoadsi. Behav Ecol Sociobiol. 2001;50:436–40.

    Article  Google Scholar 

  39. Firman R. Female social preference for males that have evolved via monogamy: evidence of a trade-off between pre- and post-copulatory sexually selected traits? Biol Lett. 2014;10:20140659.

    Article  Google Scholar 

  40. Lopes PC, König B. Choosing a healthy mate: sexually attractive traits as reliable indicators of current disease status in house mice. Anim Behav. 2016;111:119–26.

    Article  Google Scholar 

  41. Byers SL, Wiles MV, Dunn SL, Taft RA. Mouse estrous cycle identification tool and images. PLoS One. 2012;7:e35538.

    Article  CAS  Google Scholar 

  42. Friard O, Gamba M. BORIS: A free, versatile open-source event-logging software for video/audio coding and live observations. Methods Ecol Evol. 2016;7:1325–30.

    Article  Google Scholar 

  43. Drickamer LC, Gowaty PA, Holmes CM. Free female mate choice in house mice affects reproductive success and offspring viability and performance. Anim Behav. 2000;59:371–8.

    Article  CAS  Google Scholar 

  44. Manser A, König B, Lindholm AK. Female house mice avoid fertilization by t haplotype incompatible males in a mate choice experiment. J Evol Biol. 2015;28:54–64.

    Article  CAS  Google Scholar 

  45. Davison AC, Hinkley DV. Bootstrap methods and their applications. Cambridge: Cambridge University Press; 1997.

    Book  Google Scholar 

  46. Canty A, Ripley B. Boot: Bootstrap R (S-plus) functions. R package version. 2017;1:3–20.

    Google Scholar 

  47. Flachs P, Bhattacharyya T, Mihola O, Piálek J, Forejt J, Trachtulec Z. Prdm9 incompatibility controls oligospermia and delayed fertility but no selfish transmission in mouse intersubspecific hybrids. PLoS One. 2014;9:e95806.

    Article  Google Scholar 

  48. Harris MJ, Wallace ME, Evans EP. Aneuploidy in the embryonic progeny of females heterozygous for the Robertsonian chromosome (9.12) in genetically wild Peru-Coppock mice (Mus musculus). J Reprod Fertil. 1986;76:193–203.

    Article  CAS  Google Scholar 

  49. Epstein C. Mouse monosomies and trisomies as experimental systems for studying mammalian aneuploidy. Trends Genet. 1985;1:129–34.

    Article  Google Scholar 

  50. Magnuson T, Debrot S, Dimpfl J, Zweig A, Zamora T, Epstein C. The early lethality of autosomal monosomy in the mouse. J Exp Zool. 1985;236:353–60.

    Article  CAS  Google Scholar 

  51. Winking H, Dulíc B, Bulfield G. Robertsonian karyotype variation in the European house mouse, Mus musculus: survey of present knowledge and new observations. Z Säugetierkd. 1988;53:148–61.

    Google Scholar 

  52. Gropp A, Winking H, Redi C. Consequences of Robertsonian heterozygosity: segregational impairment of fertility versus male-limited sterility. In: Crosignani PG, Rubin BL, editors. Genetic control of gamete production and function. Orlando: Grune & Stratton; 1982. p. 115–34.

    Google Scholar 

  53. Haldane JBS. Sex ratio and unisexual sterility in hybrid animals. J Genet. 1922;12:101–9.

    Article  Google Scholar 

  54. Hunt P, Hassold T. Sex matters in meiosis. Science. 2002;296:2181–3.

    Article  CAS  Google Scholar 

  55. Baudat F, Imai Y, de Massy B. Meiotic recombination in mammals: localization and regulation. Nat Rev Genet. 2013;14:794.

    Article  CAS  Google Scholar 

  56. Castiglia R, Capanna E. Chiasma repatterning across a chromosomal hybrid zone between chromosomal races of Mus musculus domesticus. Genetica. 2002;114:35–40.

    Article  Google Scholar 

  57. Dumas D, Catalan J, Britton-Davidian J. Reduced recombination patterns in Robertsonian hybrids between chromosomal races of the house mouse: chiasma analyses. Heredity. 2015;114:56–64.

    Article  CAS  Google Scholar 

  58. Bidau CJ, Giménez MD, Palmer CL, Searle JB. The effects of Robertsonian fusions on chiasma frequency and distribution in the house mouse (Mus musculus domesticus) from a hybrid zone in northern Scotland. Heredity. 2001;87:305–13.

    Article  CAS  Google Scholar 

  59. Dougherty L, Shuker D. The effect of experimental design on the measurement of mate choice: a meta-analysis. Behav Ecol. 2015;26:311–9.

    Article  Google Scholar 

  60. Nunes AC, Mathias ML, Ganem G. Odor preference in house mice: influences of habitat heterogeneity and chromosomal incompatibility. Behav Ecol. 2009;20:1252–61.

    Article  Google Scholar 

  61. Capanna E, Corti M, Mainardi D, Parmigiani S, Brain PF. Karyotype and intermale aggression in wild house mice: ecology and speciation. Behav Genet. 1984;14:195–208.

    Article  CAS  Google Scholar 

  62. Ganem G, Searle JB. Behavioural discrimination among chromosomal races of the house mouse (Mus musculus domesticus). J Evol Biol. 1996;9:817–30.

    Article  Google Scholar 

  63. Carpineti M, Castiglia R. Analysis of behavioural discrimination mechanisms in a contact zone between two metacentric races of the house mouse, Mus musculus domesticus, in Central Italy. Rend Fis. 2004;15:31.

    Article  Google Scholar 

  64. Gropp A, Winking H, Zech L, Müller H. Robertsonian chromosomal variation and identification of metacentric chromosomes in feral mice. Chromosoma. 1972;39:265–88.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We warmly thank Barbara König for her support throughout and landowners who allowed us to trap mice at their properties and the many trapping volunteers.

Funding

Funding was provided by the University Research Priority Program “Evolution in Action” of the University of Zurich, the Claraz-Stiftung, the Promotor-Stiftung, and the University of Zurich.

Availability of data and materials

All data generated or analysed during this study are included in the Additional files 3, 4, 5, 6, 7.

Author information

Authors and Affiliations

Authors

Contributions

AL and SG designed the study. EW performed the preference test. EW and SG trapped the mice. SG performed the crosses, fertility analyses and statistics.

SG, JS and AL wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Anna K. Lindholm.

Ethics declarations

Ethics approval and consent to participate

Experiments were performed in accordance with the Swiss law and with the ethics approval of the Veterinary Office Kanton Zürich, Switzerland (licence no. 64/2014).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Duration of female backcrosses. (PDF 11 kb)

Additional file 2:

Linear regressions between male body mass and the variables sperm number, testis mass, epididymis mass, mass of seminal vesicle with coagulating gland and mass of preputial gland. (PDF 136 kb)

Additional file 3:

Data of male fertility. (XLSX 20 kb)

Additional file 4:

Data of female fertility. (XLSX 13 kb)

Additional file 5:

Data of litter sizes for all litters of hybrid and control females. (XLSX 14 kb)

Additional file 6:

Data of sperm counts of males used in female backcrosses. (XLSX 10 kb)

Additional file 7:

Data of female preference test. (XLSX 17 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Grize, S.A., Wilwert, E., Searle, J.B. et al. Measurements of hybrid fertility and a test of mate preference for two house mouse races with massive chromosomal divergence. BMC Evol Biol 19, 25 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-018-1322-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12862-018-1322-y

Keywords