Characterization of a centromeric marker on mouse Chromosome 11 and its introgression in a domesticus/musculus hybrid zone

Transmission distortion and genetic incompatibilities between alleles in a multigenerational mouse advanced intercross line

While direct additive and dominance effects on complex traits have been mapped repeatedly, additional genetic factors contributing to the heterogeneity of complex traits have been scarcely investigated. To assess genetic background effects, we investigated transmission ratio distortions (TRDs) of alleles from parent to offspring using an advanced intercross line (AIL) of an initial cross between the mouse inbred strains C57BL/6NCrl (B6N) and BFMI860-12 [Berlin Fat Mouse Inbred (BFMI)]. A total of 341 males of generation 28 and their respective 61 parents and 66 grandparents were genotyped using Mega Mouse Universal Genotyping Arrays. TRDs were investigated using allele transmission asymmetry tests, and pathway overrepresentation analysis was performed. Sequencing data were used to test for overrepresentation of nonsynonymous SNPs (nsSNPs) in TRD regions. Genetic incompatibilities were tested using the Bateson-Dobzhansky-Muller two-locus model. A total of 62 TRD regions were detected, many in close proximity to the telocentric centromere. TRD regions contained 44.5% more nsSNPs than randomly selected regions (182 vs 125.9 ± 17.0, P < 1 × 10-4). Testing for genetic incompatibilities between TRD regions identified 29 genome-wide significant incompatibilities between TRD regions [P(BF) < 0.05]. Pathway overrepresentation analysis of genes in TRD regions showed that DNA methylation, epigenetic regulation of RNA, and meiotic/meiosis regulation pathways were affected independent of the parental origin of the TRD. Paternal BFMI TRD regions showed overrepresentation in the small interfering RNA biogenesis and in the metabolism of lipids and lipoproteins. Maternal B6N TRD regions harbored genes involved in meiotic recombination, cell death, and apoptosis pathways. The analysis of genes in TRD regions suggests the potential distortion of protein-protein interactions influencing obesity and diabetic retinopathy as a result of disadvantageous combinations of allelic variants in Aass, Pgx6, and Nme8. Using an AIL significantly improves the resolution at which we can investigate TRD. Our analysis implicates distortion of protein-protein interactions as well as meiotic drive as the underlying mechanisms leading to the observed TRD in our AIL. Furthermore, genes with large amounts of nsSNPs located in TRD regions are more likely to be involved in pathways that are related to the phenotypic differences between the parental strains. Genes in these TRD regions provide new targets for investigating genetic adaptation, protein-protein interactions, and determinants of complex traits such as obesity.

Evidence of functional Cd94 polymorphism in a free-living house mouse population

The CD94 receptor, expressed on natural killer (NK) and CD8+ T cells, is known as a relatively non-polymorphic receptor with orthologues in humans, other primates, cattle, and rodents. In the house mouse (Mus musculus), a single allele is highly conserved among laboratory strains, and reports of allelic variation in lab- or wild-living mice are lacking, except for deficiency in one lab strain (DBA/2J). The non-classical MHC-I molecule Qa-1b is the ligand for mouse CD94/NKG2A, presenting alternative non-americ fragment of leader peptides (Qa-1 determinant modifier (Qdm)) from classical MHC-I molecules. Here, we report a novel allele identified in free-living house mice captured in Norway, living among individuals carrying the canonical Cd94 allele. The novel Cd94^LocA allele encodes 12 amino acid substitutions in the extracellular lectin-like domain. Flow cytometric analysis of primary NK cells and transfected cells indicates that the substitutions prevent binding of CD94 mAb and Qa-1b/Qdm tetramers. Our data further indicate correlation of Cd94 polymorphism with the two major subspecies of house mice in Europe. Together, these findings suggest that the Cd94^LocA/NKG2A heterodimeric receptor is widely expressed among M. musculus subspecies musculus, with ligand-binding properties different from mice of subspecies domesticus, such as the C57BL/6 strain.

Molecular and phenotypic distinction of the very recently evolved insular subspecies Mus musculus helgolandicus ZIMMERMANN, 1953

Background

Populations and subspecies of the house mouse Mus musculus were able to invade new regions worldwide in the wake of human expansion. Here we investigate the origin and colonization history of the house mouse inhabiting the small island of Heligoland on the German Bight - Mus musculus helgolandicus. It was first described by Zimmermann in 1953, based on morphological descriptions which were considered to be a mosaic between the subspecies M. m. domesticus and M. m. musculus. Since mice on islands are excellent evolutionary model systems, we have focused here on a molecular characterization and an extended phenotype analysis.

Results

The molecular data show that the mice from Heligoland are derived from M. m. domesticus based on mitochondrial D-loop sequences as well as on four nuclear diagnostic markers, including one each from the sex-chromosomes. STRUCTURE analysis based on 21 microsatellite markers assigns Heligoland mice to a distinct population and D-loop network analysis suggests that they are derived from a single colonization event. In spite of mice from the mainland arriving by ships, they are apparently genetically refractory against further immigration. Mutation frequencies in complete mitochondrial genome sequences date the colonization age to approximately 400 years ago. Complete genome sequences from three animals revealed a genomic admixture with M. m. musculus genomic regions with at least 6.5 % of the genome affected. Geometric morphometric analysis of mandible shapes including skull samples from two time points during the last century suggest specific adaptations to a more carnivorous diet.

Conclusions

The molecular and morphological analyses confirm that M. m. helgolandicus consists of a distinct evolutionary lineage with specific adaptations. It shows a remarkable resilience against genetic mixture with mainland populations of M. m. domesticus despite major disturbances in the past century and a high ship traffic. The genomic admixture with M. m. musculus genetic material may have contributed to the genomic distinction of the Heligoland mice. In spite of its young age, M. m. helgolandicus may thus be considered as a true subspecies of Mus, whose evolution was triggered through fast divergence on a small island.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0439-5) contains supplementary material, which is available to authorized users.

Norwegian house mice (Mus musculus musculus/domesticus): distributions, routes of colonization and patterns of hybridization

We investigated the distributions and routes of colonization of two commensal subspecies of house mouse in Norway: Mus musculus domesticus and M. m. musculus. Five nuclear markers (Abpa, D11 cenB2, Btk, SMCY and Zfy2) and a morphological feature (tail length) were used to differentiate the two subspecies and assess their distributions, and mitochondrial (mt) D-loop sequences helped to elucidate their colonization history. M. m. domesticus is the more widespread of the two subspecies, occupying the western and southern coast of Norway, while M. m. musculus is found along Norway's southeastern coast and east from there to Sweden. Two sections of the hybrid zone between the two subspecies were localized in Norway. However, hybrid forms also occur well away from that hybrid zone, the most prevalent of which are mice with a M. m. musculus-type Y chromosome and an otherwise M. m. domesticus genome. MtDNA D-loop sequences of the mice revealed a complex phylogeography within M. m. domesticus, reflecting passive human transport to Norway, probably during the Viking period. M. m. musculus may have colonized earlier. If so, that leaves open the possibility that M. m. domesticus replaced M. m. musculus from much of Norway, with the widely distributed hybrids a relict of this process. Overall, the effects of hybridization are evident in house mice throughout Norway.

The diverse origins of New Zealand house mice

Molecular markers and morphological characters can help infer the colonization history of organisms. A combination of mitochondrial (mt) D-loop DNA sequences, nuclear DNA data, external measurements and skull characteristics shows that house mice (Mus musculus) in New Zealand and its outlying islands are descended from very diverse sources. The predominant genome is Mus musculus domesticus (from western Europe), but Mus musculus musculus (from central Europe) and Mus musculus castaneus (from southern Asia) are also represented genetically. These subspecies have hybridized to produce combinations of musculus and domesticus nuclear DNA coupled with domesticus mtDNA, and castaneus or musculus mtDNA with domesticus nuclear DNA. The majority of the mice with domesticus mtDNA that we sampled had D-loop sequences identical to two haplotypes common in Britain. This is consistent with long-term British-New Zealand cultural linkages. The origins of the castaneus mtDNA sequences widespread in New Zealand are less easy to identify.

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.