The effect of host admixture on wild house mouse gut microbiota is weak when accounting for spatial autocorrelation.

The question of how interactions between the gut microbiome and vertebrate hosts contribute to host adaptation and speciation is one of the major problems in current evolutionary research. Using bacteriome and mycobiome metabarcoding, we examined how these two components of the gut microbiota vary with the degree of host admixture in secondary contact between two house mouse subspecies (Mus musculus musculus and M. m. domesticus). We used a large data set collected at two replicates of the hybrid zone and model-based statistical analyses to ensure the robustness of our results. Assuming that the microbiota of wild hosts suffers from spatial autocorrelation, we directly compared the results of statistical models that were spatially naive with those that accounted for spatial autocorrelation. We showed that neglecting spatial autocorrelation can strongly affect the results and lead to misleading conclusions. The spatial analyses showed little difference between subspecies, both in microbiome composition and in individual bacterial lineages. Similarly, the degree of admixture had minimal effects on the gut bacteriome and mycobiome and was caused by changes in a few microbial lineages that correspond to the common symbionts of free-living house mice. In contrast to previous studies, these data do not support the hypothesis that the microbiota plays an important role in host reproductive isolation in this particular model system.

New Perspective on the Geographic Distribution and Evolution of Lymphocytic Choriomeningitis Virus, Central Europe.

Lymphocytic choriomeningitis virus (LCMV) is an Old World mammarenavirus found worldwide because of its association with the house mouse. When LCMV spills over to immunocompetent humans, the virus can cause aseptic meningitis; in immunocompromised persons, systemic infection and death can occur. Central Europe is a strategic location for the study of LCMV evolutionary history and host specificity because of the presence of a hybrid zone (genetic barrier) between 2 house mouse subspecies, Mus musculus musculus and M. musculus domesticus. We report LCMV prevalence in natural mouse populations from a Czech Republic-Germany transect and genomic characterization of 2 new LCMV variants from the Czech Republic. We demonstrate that the main division in the LCMV phylogenetic tree corresponds to mouse host subspecies and, when the virus is found in human hosts, the mouse subspecies found at the spillover location. Therefore, LCMV strains infecting humans can be predicted by the genetic structure of house mice.

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.

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.

The strength of gut microbiota transfer along social networks and genealogical lineages in the house mouse.

The gut microbiota of vertebrates is acquired from the environment and other individuals, including parents and unrelated conspecifics. In the laboratory mouse, a key animal model, inter-individual interactions are severely limited and its gut microbiota is abnormal. Surprisingly, our understanding of how inter-individual transmission impacts house mouse gut microbiota is solely derived from laboratory experiments. We investigated the effects of inter-individual transmission on gut microbiota in two subspecies of house mice (Mus musculus musculus and M. m. domesticus) raised in a semi-natural environment without social or mating restrictions. We assessed the correlation between microbiota composition (16S rRNA profiles), social contact intensity (microtransponder-based social networks), and mouse relatedness (microsatellite-based pedigrees). Inter-individual transmission had a greater impact on the lower gut (colon and cecum) than on the small intestine (ileum). In the lower gut, relatedness and social contact independently influenced microbiota similarity. Despite female-biased parental care, both parents exerted a similar influence on their offspring's microbiota, diminishing with the offspring's age in adulthood. Inter-individual transmission was more pronounced in M. m. domesticus, a subspecies, with a social and reproductive network divided into more closed modules. This suggests that the transmission magnitude depends on the social and genetic structure of the studied population.

Aberrant microbiomes are associated with increased antibiotic resistance gene load in hybrid mice.

Antibiotic resistance is a priority public health problem resulting from eco-evolutionary dynamics within microbial communities and their interaction at a mammalian host interface or geographical scale. The links between mammalian host genetics, bacterial gut community, and antimicrobial resistance gene (ARG) content must be better understood in natural populations inhabiting heterogeneous environments. Hybridization, the interbreeding of genetically divergent populations, influences different components of the gut microbial communities. However, its impact on bacterial traits such as antibiotic resistance is unknown. Here, we present that hybridization might shape bacterial communities and ARG occurrence. We used amplicon sequencing to study the gut microbiome and to predict ARG composition in natural populations of house mice (Mus musculus). We compared gastrointestinal bacterial and ARG diversity, composition, and abundance across a gradient of pure and hybrid genotypes in the European House Mouse Hybrid Zone. We observed an increased overall predicted richness of ARG in hybrid mice. We found bacteria-ARG interactions by their co-abundance and detected phenotypes of extreme abundances in hybrid mice at the level of specific bacterial taxa and ARGs, mainly multidrug resistance genes. Our work suggests that mammalian host genetic variation impacts the gut microbiome and chromosomal ARGs. However, it raises further questions on how the mammalian host genetics impact ARGs via microbiome dynamics or environmental covariates.

The Expansion of House Mouse Major Urinary Protein Genes Likely Did Not Facilitate Commensalism with Humans.

Mouse wild-derived strains (WDSs) combine the advantages of classical laboratory stocks and wild animals, and thus appear to be promising tools for diverse biomedical and evolutionary studies. We employed 18 WDSs representing three non-synanthropic species (Mus spretus, Mus spicilegus, and M. macedonicus) and three house mouse subspecies (Mus musculus musculus, M. m. domesticus, M. m. castaneus), which are all important human commensals to explore whether the number of major urinary protein (MUP) genes and their final protein levels in urine are correlated with the level of commensalism. Contrary to expectations, the MUP copy number (CN) and protein excretion in the strains derived from M. m. castaneus, which is supposed to be the strongest commensal, were not significantly different from the non-commensal species. Regardless of an overall tendency for higher MUP amounts in taxa with a higher CN, there was no significant correlation at the strain level. Our study thus suggests that expansion of the Mup cluster, which appeared before the house mouse diversification, is unlikely to facilitate commensalism with humans in three house mouse subspecies. Finally, we found considerable variation among con(sub)specific WDSs, warning against generalisations of results based on a few strains.

A gene copy number arms race in action: X,Y-chromosome transmission distortion across a species barrier.

A remarkable gene copy number (CN) arms race system has recently been described in laboratory mice, where Slx;Slxl1 and Sly genes compete over transmission by altering the fertilization success of X and Y chromosome-bearing sperm, respectively. Here, we focus on this system in nature, where natural selection can counter CN/gene product escalation. Our model is house mouse subspecies hybridizing in Europe. In some regions, Y chromosomes of the Eastern subspecies have introgressed onto Western genomic backgrounds, accompanied by sex ratio distortion in favor of males, consistent with the inbred lines suggested mechanism: Overabundance of SLY protein expressed by invading Y chromosomes. We take Slx as representative of the X side of this arms race and measure Slx|Sly CN and expression across an "Invasion" transect where Ys introgress and a "Control" transect with negligible introgression. Since we found similar Slx|Sly ratios in both transects, SLY overabundance is unlikely to explain the introgression. However, Slx CN is relatively low in the introgression area, suggesting that Slx is less able to combat Sly effects here. Furthermore, deterministic changes in Slx;Sly expression proportions versus CN proportions suggest standing variation for trans regulation of Slx|Sly is being co-opted in nature where their arms race reduces population fitness.

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse.

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane's rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.

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.

House mouse subspecies do differ in their social structure.

It is widely acknowledged that population structure can have a substantial impact on evolutionary trajectories. In social animals, this structure is strongly influenced by relationships among the population members, so studies of differences in social structure between diverging populations or nascent species are of prime interest. Ideal models for such a study are two house mouse subspecies, Mus musculus musculus and M. m. domesticus, meeting in Europe along a secondary contact zone. Though the latter subspecies has usually been supposed to form tighter and more isolated social units than the former, the evidence is still inconclusive. Here, we carried out a series of radiofrequency identification experiments in semi-natural enclosures to gather large longitudinal data sets on individual mouse movements. The data were summarized in the form of uni- and multi-layer social networks. Within them, we could delimit and describe the social units ("modules"). While the number of estimated units was similar in both subspecies, domesticus revealed a more "modular" structure. This subspecies also showed more intramodular social interactions, higher spatial module separation, higher intramodular persistence of parent-offspring contacts, and lower multiple paternity, suggesting more effective control of dominant males over reproduction. We also demonstrate that long-lasting modules can be identified with basic reproductive units or demes. We thus provide the first robust evidence that the two subspecies differ in their social structure and dynamics of the structure formation.

Divergent gut microbiota in two closely related house mouse subspecies under common garden conditions.

The gastrointestinal microbiota (GM) is considered an important component of the vertebrate holobiont. GM-host interactions influence the fitness of holobionts and are, therefore, an integral part of evolution. The house mouse is a prominent model for GM-host interactions, and evidence suggests a role for GM in mouse speciation. However, previous studies based on short 16S rRNA GM profiles of wild house mouse subspecies failed to detect GM divergence, which is a prerequisite for the inclusion of GM in Dobzhansky-Muller incompatibilities. Here, we used standard 16S rRNA GM profiling in two mouse subspecies, Mus musculus musculus and M. m. domesticus, including the intestinal mucosa and content of three gut sections (ileum, caecum, and colon). We reduced environmental variability by sampling GM in the offspring of wild mice bred under seminatural conditions. Although the breeding conditions allowed a contact between the subspecies, we found a clear differentiation of GM between them, in all three gut sections. Differentiation was mainly driven by several Helicobacters and two H. ganmani variants showed a signal of codivergence with their hosts. Helicobacters represent promising candidates for studying GM-host coadaptations and the fitness effects of their interactions.

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.

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.

Ontogeny of social hierarchy in two European house mouse subspecies and difference in the social rank of dispersing males.

In social species such as house mouse, being dominant is vital. Determination of dominance may start early in life and vary during ontogeny. We asked whether pre-pubertal and adolescent behaviour predicts the rank a male mouse finally obtains. Moreover, we asked how dominant vs. subordinate adults differ in exploration and propensity to emigrate. We studied fraternal pairs as the simple social units, from weaning to full-grown adulthood. By utilizing two mouse subspecies known to differ in many behavioural traits, we take into account any potential subspecific idiosyncrasies. We did not find any significant effect of future social status on any behavioural type displayed before adulthood, but the subspecies themselves differ in behaviours prevailing in particular ontogeny phases. While musculus males start as more pro-social, they later became significantly more passive. Conversely, domesticus are slightly less passive at the beginning but significantly more proactive close to adulthood and rapidly establishing hierarchy through overt conflicts. We found no difference in exploration between ranks, however, domesticus males were significantly more active in an unknown area than musculus. Most importantly, while dominant domesticus males seem to be more prone to emigration, in musculus it was the subordinate males who left base significantly more often. This is consistent with extended contests of musculus males over dominance found in this study as well as with differences in endocrinological changes we have reported previously.

Loss of Sprouty Produces a Ciliopathic Skeletal Phenotype in Mice Through Upregulation of Hedgehog Signaling.

The Sprouty family is a highly conserved group of intracellular modulators of receptor tyrosine kinase (RTK)-signaling pathways, which have been recently linked to primary cilia. Disruptions in the structure and function of primary cilia cause inherited disorders called ciliopathies. We aimed to evaluate Sprouty2 and Sprouty4 gene-dependent alterations of ciliary structure and to focus on the determination of its association with Hedgehog signaling defects in chondrocytes. Analysis of the transgenic mice phenotype with Sprouty2 and Sprouty4 deficiency revealed several defects, including improper endochondral bone formation and digit patterning, or craniofacial and dental abnormalities. Moreover, reduced bone thickness and trabecular bone mass, skull deformities, or chondroma-like lesions were revealed. All these pathologies might be attributed to ciliopathies. Elongation of the ciliary axonemes in embryonic and postnatal growth plate chondrocytes was observed in Sprouty2-/- and Sprouty2+/- /Sprouty4-/- mutants compared with corresponding littermate controls. Also, cilia-dependent Hedgehog signaling was upregulated in Sprouty2/4 mutant animals. Ptch1 and Ihh expression were upregulated in the autopodium and the proximal tibia of Sprouty2-/- /Sprouty4-/- mutants. Increased levels of the GLI3 repressor (GLI3R) form were detected in Sprouty2/4 mutant primary fibroblast embryonic cell cultures and tissues. These findings demonstrate that mouse lines deficient in Sprouty proteins manifest phenotypic features resembling ciliopathic phenotypes in multiple aspects and may serve as valuable models to study the association between overactivation of RTK and dysfunction of primary cilia during skeletogenesis. © 2021 American Society for Bone and Mineral Research (ASBMR).

Polarized Sonic Hedgehog Protein Localization and a Shift in the Expression of Region-Specific Molecules Is Associated With the Secondary Palate Development in the Veiled Chameleon.

Secondary palate development is characterized by the formation of two palatal shelves on the maxillary prominences, which fuse in the midline in mammalian embryos. However, in reptilian species, such as turtles, crocodilians, and lizards, the palatal shelves of the secondary palate develop to a variable extent and morphology. While in most Squamates, the palate is widely open, crocodilians develop a fully closed secondary palate. Here, we analyzed developmental processes that underlie secondary palate formation in chameleons, where large palatal shelves extend horizontally toward the midline. The growth of the palatal shelves continued during post-hatching stages and closure of the secondary palate can be observed in several adult animals. The massive proliferation of a multilayered oral epithelium and mesenchymal cells in the dorsal part of the palatal shelves underlined the initiation of their horizontal outgrowth, and was decreased later in development. The polarized cellular localization of primary cilia and Sonic hedgehog protein was associated with horizontal growth of the palatal shelves. Moreover, the development of large palatal shelves, supported by the pterygoid and palatine bones, was coupled with the shift in Meox2, Msx1, and Pax9 gene expression along the rostro-caudal axis. In conclusion, our results revealed distinctive developmental processes that contribute to the expansion and closure of the secondary palate in chameleons and highlighted divergences in palate formation across amniote species.

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.

Genetic conflict outweighs heterogametic incompatibility in the mouse hybrid zone?

BACKGROUND: The Mus musculus musculus/M. m. domesticus contact zone in Europe is characterised by sharp frequency discontinuities for sex chromosome markers at the centre of wider clines in allozyme frequencies. RESULTS: We identify a triangular area (approximately 330 km2) where the musculus Y chromosome introgresses across this front for up to 22 km into domesticus territory. Introgression of the Y chromosome is accompanied by a perturbation of the census sex ratio: the sex ratio is significantly female biased in musculus localities and domesticus localities lacking Y chromosome introgression. In contrast, where the musculus Y is detected in domesticus localities, the sex ratio is close to parity, and significantly different from both classes of female biased localities. The geographic position of an abrupt cline in an X chromosome marker, and autosomal clines centred on the same position, seem unaffected by the musculus Y introgression. CONCLUSION: We conclude that sex ratio distortion is playing a role in the geographic separation of speciation genes in this section of the mouse hybrid zone. We suggest that clines for genes involved in sex-ratio distortion have escaped from the centre of the mouse hybrid zone, causing a decay in the barrier to gene flow between the two house mouse taxa.

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.