The location of theMus musculus/M. domesticus hybrid zone in the Balkans: clues from morphology

Disrupted Gene Networks in Subfertile Hybrid House Mice

The Dobzhansky–Muller (DM) model provides a widely accepted mechanism for the evolution of reproductive isolation: incompatible substitutions disrupt interactions between genes. To date, few candidate incompatibility genes have been identified, leaving the genes driving speciation mostly uncharacterized. The importance of interactions in the DM model suggests that gene coexpression networks provide a powerful framework to understand disrupted pathways associated with postzygotic isolation. Here, we perform weighted gene coexpression network analysis to infer gene interactions in hybrids of two recently diverged European house mouse subspecies, Mus mus domesticus and M. m. musculus, which commonly show hybrid male sterility or subfertility. We use genome-wide testis expression data from 467 hybrid mice from two mapping populations: F₂s from a laboratory cross between wild-derived pure subspecies strains and offspring of natural hybrids captured in the Central Europe hybrid zone. This large data set enabled us to build a robust consensus network using hybrid males with fertile phenotypes. We identify several expression modules, or groups of coexpressed genes, that are disrupted in subfertile hybrids, including modules functionally enriched for spermatogenesis, cilium and sperm flagellum organization, chromosome organization, and DNA repair, and including genes expressed in spermatogonia, spermatocytes, and spermatids. Our network-based approach enabled us to hone in on specific hub genes likely to be influencing module-wide gene expression and hence potentially driving large-effect DM incompatibilities. A disproportionate number of hub genes lie within sterility loci identified previously in the hybrid zone mapping population and represent promising candidate barrier genes and targets for future functional analysis.

Sperm quality, aggressiveness and generation turnover may facilitate unidirectional Y chromosome introgression across the European house mouse hybrid zone

The widespread and locally massive introgression of Y chromosomes of the eastern house mouse (Mus musculus musculus) into the range of the western subspecies (M. m. domesticus) in Central Europe calls for an explanation of its underlying mechanisms. Given the paternal inheritance pattern, obvious candidates for traits mediating the introgression are characters associated with sperm quantity and quality. We can also expect traits such as size, aggression or the length of generation cycles to facilitate the spread. We have created two consomic strains carrying the non-recombining region of the Y chromosome of the opposite subspecies, allowing us to study introgression in both directions, something impossible in nature due to the unidirectionality of introgression. We analyzed several traits potentially related to male fitness. Transmission of the domesticus Y onto the musculus background had negative effects on all studied traits. Likewise, domesticus males possessing the musculus Y had, on average, smaller body and testes and lower sperm count than the parental strain. However, the same consomic males tended to produce less- dissociated sperm heads, to win more dyadic encounters, and to have shorter generation cycles than pure domesticus males. These data suggest that the domesticus Y is disadvantageous on the musculus background, while introgression in the opposite direction can confer a recognizable, though not always significant, selective advantage. Our results are thus congruent with the unidirectional musculus → domesticus Y chromosome introgression in Central Europe. In addition to some previous studies, they show this to be a multifaceted phenomenon demanding a multidisciplinary approach.

Intensity of infection with intracellular Eimeria spp. and pinworms is reduced in hybrid mice compared to parental subspecies

Genetic diversity in animal immune systems is usually beneficial. In hybrid recombinants, this is less clear, as the immune system could also be impacted by genetic conflicts. In the European house mouse hybrid zone, the long-standing impression that hybrid mice are more highly parasitized and less fit than parentals persists despite the findings of recent studies. Working across a novel transect, we assessed infections by intracellular protozoans, Eimeria spp., and infections by extracellular macroparasites, pinworms. For Eimeria, we found lower intensities in hybrid hosts than in parental mice but no evidence of lowered probability of infection or increased mortality in the centre of the hybrid zone. This means ecological factors are very unlikely to be responsible for the reduced load of infected hybrids. Focusing on parasite intensity (load in infected hosts), we also corroborated reduced pinworm loads reported for hybrid mice in previous studies. We conclude that intensity of diverse parasites, including the previously unstudied Eimeria, is reduced in hybrid mice compared to parental subspecies. We suggest caution in extrapolating this to differences in hybrid host fitness in the absence of, for example, evidence for a link between parasitemia and health.

Small terrestrial mammals of Albania: distribution and diversity (Mammalia, Eulipotyphla, Rodentia)

In this paper new records are reported for 23 species of small terrestrial mammals (STM) of Albania collected during the field work campaigns organised in the framework of the project “Strengthening capacity in National Nature Protection – preparation for Natura 2000 network” (NaturAL) in Albania during the summer and autumn of 2016 and 2017 Data on small mammals were primarily collected through Sherman live-trapping campaigns in six high priority protected areas of Albania: Korab-Koritnik, Bredhi i Hotovës, Tomorri, Llogara-Karaburun, Divjakë-Karavasta, Liqeni i Shkodrës (Skadar lake), Lëpushë-Vermosh. Other data were obtained by analysis of owl pellets or by direct observation of individuals (dead or alive) in the field. For 21 species Erinaceusroumanicus, Neomysanomalus, Crocidurasuaveolens, Crociduraleucodon, Suncusetruscus, Talpastankovici/caeca, Myocastorcoypus, Sciurusvulgaris, Glisglis, Dryomysnitedula, Muscardinusavellanarius, Microtuslevis/arvalis, Microtussubterraneus, Microtusthomasi, Microtusfelteni, Myodesglareolus, Apodemussylvaticus, Apodemusflavicollis, Apodemusepimelas, Musmusculus, and Musmacedonicus additional records are provide and their distributions reviewed, while the presence of two new species of shrews (Sorexaraneus and Sorexminutus) for Albania is reported for the first time. A comprehensive review of the published and unpublished distribution records of STM species of the country is made, together with an updated checklist and distribution maps of the species.

Coevolution of Cryptosporidium tyzzeri and the house mouse (Mus musculus)

Two house mouse subspecies occur in Europe, eastern and northern Mus musculus musculus (Mmm) and western and southern Mus musculus domesticus (Mmd). A secondary hybrid zone occurs where their ranges meet, running from Scandinavia to the Black Sea. In this paper, we tested a hypothesis that the apicomplexan protozoan species Cryptosporidium tyzzeri has coevolved with the house mouse. More specifically, we assessed to what extent the evolution of this parasite mirrors divergence of the two subspecies. In order to test this hypothesis, we analysed sequence variation at five genes (ssrRNA, Cryptosporidium oocyst wall protein (COWP), thrombospondin-related adhesive protein of Cryptosporidium 1 (TRAP-C1), actin and gp60) in C. tyzzeri isolates from Mmd and Mmm sampled along a transect across the hybrid zone from the Czech Republic to Germany. Mmd samples were supplemented with mice from New Zealand. We found two distinct isolates of C. tyzzeri, each occurring exclusively in one of the mouse subspecies (C. tyzzeri-Mmm and C. tyzzeri-Mmd). In addition to genetic differentiation, oocysts of the C. tyzzeri-Mmd subtype (mean: 4.24×3.69μm) were significantly smaller than oocysts of C. tyzzeri-Mmm (mean: 4.49×3.90 μm). Mmm and Mmd were susceptible to experimental infection with both C. tyzzeri subtypes; however, the subtypes were not infective for the rodent species Meriones unguiculatus, Mastomys coucha, Apodemus flavicollis or Cavia porcellus. Overall, our results support the hypothesis that C. tyzzeri is coevolving with Mmm and Mmd.

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

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

Reinforcement selection acting on the European house mouse hybrid zone

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent-specific mate-recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate-recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.

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

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

A candidate subspecies discrimination system involving a vomeronasal receptor gene with different alleles fixed in M. m. domesticus and M. m. musculus

Assortative mating, a potentially efficient prezygotic reproductive barrier, may prevent loss of genetic potential by avoiding the production of unfit hybrids (i.e., because of hybrid infertility or hybrid breakdown) that occur at regions of secondary contact between incipient species. In the case of the mouse hybrid zone, where two subspecies of Mus musculus (M. m. domesticus and M. m. musculus) meet and exchange genes to a limited extent, assortative mating requires a means of subspecies recognition. We based the work reported here on the hypothesis that, if there is a pheromone sufficiently diverged between M. m. domesticus and M. m. musculus to mediate subspecies recognition, then that process must also require a specific receptor(s), also sufficiently diverged between the subspecies, to receive the signal and elicit an assortative mating response. We studied the mouse V1R genes, which encode a large family of receptors in the vomeronasal organ (VNO), by screening Perlegen SNP data and identified one, Vmn1r67, with 24 fixed SNP differences most of which (15/24) are nonsynonymous nucleotide substitutions between M. m. domesticus and M. m. musculus. We observed substantial linkage disequilibrium (LD) between Vmn1r67 and Abpa27, a mouse salivary androgen-binding protein gene that encodes a proteinaceous pheromone (ABP) capable of mediating assortative mating, perhaps in conjunction with its bound small lipophilic ligand. The LD we observed is likely a case of association rather than residual physical linkage from a very recent selective sweep, because an intervening gene, Vmn1r71, shows significant intra(sub)specific polymorphism but no inter(sub)specific divergence in its nucleotide sequence. We discuss alternative explanations of these observations, for example that Abpa27 and Vmn1r67 are coevolving as signal and receptor to reinforce subspecies hybridization barriers or that the unusually divergent Vmn1r67 allele was not a product of fast positive selection, but was derived from an introgressed allele, possibly from Mus spretus.

Signalling components of the house mouse mate recognition system

Subspecies-specific mate recognition may represent significant barrier to gene flow between diverged genomes potentially leading to speciation. In the house mouse, assortative mating involves the coevolution of several signals and receptors. We compared signalling ability of bedding material, faeces, urine, saliva, salivary androgen binding proteins (ABP) and combinations of urine with saliva and urine with ABP in mate choice in two wild-derived inbred strains (one of Mus musculus musculus and one of Mus musculus domesticus origin). We observed high levels of variation in assortative preferences between the two strains and sexes. The strongest preferences were observed in M. m. musculus-derived individuals in tests where urine was present either alone or as part of a composite signal target. M. m. domesticus-derived mice displayed strain-specific preferences for faeces. Saliva was the least preferred stimulus in both strains and sexes. No effect of two-compound cues was detected. We conclude that there is divergence across both the stimulus and preference parts of the recognition system for both house mouse strains. Of the tested stimuli, those that have the capacity to carry a signal for extended periods under natural conditions (such as urine and faeces) seem to be the most important substances in strain-specific recognition.

Genome-wide patterns of gene flow across a house mouse hybrid zone

Hybrid zones between closely related species or subspecies provide useful settings for studying the genetic architecture of speciation. Using markers distributed throughout the mouse genome, we use a hybrid zone between two recently diverged species of house mice (Mus musculus and Mus domesticus) as a natural mapping experiment to identify genomic regions that may be involved in reproductive isolation. Using cline analysis we document a nearly 50-fold variation in level of introgression among markers. Some markers have extremely narrow cline widths; these genomic regions may contribute to reproductive isolation. Biological processes associated with these narrow clines include physiological and immune responses to the environment as well as physiological and behavioral aspects of reproduction. Other autosomal markers exhibit asymmetrically broad clines, usually with high frequencies of M. domesticus alleles on the M. musculus side of the hybrid zone. These markers identify genome regions likely housing genes with alleles that are spreading from one species to the other. Biological processes associated with these wide clines include cell signaling, olfaction, and pheromone response. These processes play important roles in survival and reproduction, and associated genes are likely targets of selection. Patterns of linkage disequilibrium in the center of the hybrid zone suggest that isolation may be caused by multiple epistatic interactions between sets of genes. These data highlight the complex genetic architecture underlying speciation even at early stages of divergence and point to some of the biological processes that may govern this architecture.

Genetic analysis of autosomal and X-linked markers across a mouse hybrid zone

In this paper, we present results of the first comprehensive study of the introgression of both autosomal and sex-chromosome markers across the central European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and M. m. domesticus. More than 1800 individuals sampled from 105 sites were analyzed with a set of allozyme loci (hopefully representing neutral or nearly neutral markers) and X-linked loci (which are assumed to be under selection). The zone center is best modeled as a single straight line independent of fine-scale local geographic or climatic conditions, being maintained by a balance between dispersal and selection against hybrids. The width (w) of the multilocus autosomal cline was estimated as 9.6 km whereas the estimate for the compound X-chromosome cline was about 4.6 km only. As the former estimate is comparable to that of the Danish portion of the zone (assumed to be much younger than the central European one), zone width does not appear to be related to its age. The strength (B) of the central barrier was estimated as about 20 km; with dispersal (sigma) of about 1 km/gen(1/2), this means effective selection (s*) is approximately 0.06-0.09 for autosomal loci and about 0.25 for X-linked loci. The number of loci under selection was estimated as N= 56-99 for autosomes and about 380 for X-linked loci. Finally, we highlight some potential pitfalls in hybrid zone analyses and in comparisons of different transects. We suggest that conclusions about parts of the mouse genome involved in reproductive isolation and speciation should be drawn with caution and that analytical approaches always providing some estimates should not be used without due care regarding the support or confidence of such estimates, especially if conclusions are based on the difference between these estimates. Finally, we recommend that analysis in two-dimensional space, dense sampling, and rigorous treatment of data, including inspection of likelihood profiles, are essential for hybrid zone studies.