Genetic analysis of a hybrid zone between domesticus and musculus mice (Mus musculus complex): hemoglobin polymorphisms

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

SALIVARY ANDROGEN-BINDING PROTEIN (ABP) MEDIATES SEXUAL ISOLATION IN MUS MUSCULUS

We wanted to determine whether the microevolution of the mouse salivary androgen-binding protein (ABP) Alpha subunit gene (Abpa) could mediate sexual selection and thereby have a potential role in maintaining gene pool integrity where radiating mouse subspecies make secondary contact. This hypothesis is based upon previous work in this laboratory, which has shown that each subspecies apparently has its own allele and that these alleles have a 25-fold excess of nonsynonymous/synonymous base substitutions compared to an average protein under purifying selection. We provide direct evidence for ABP-assortative mate selection in a laboratory setting: Mus musculus domesticus and M. m. musculus female mice recognize and discriminate between the territories of male mice that essentially differ solely in their Abpa genotype and, when the males are present, the female prefers to mate with the one of her own ABP type. The observation that females could differentiate between the territories of the two males when those mice were absent suggests that the males marked their territories with ABP. In this study, we also detected ABP on the pelts of male mice and in their environment. It is likely that the animals apply the protein to their pelts by licking and that it is then deposited in their surroundings. We suggest that females of the two subspecies are able to discriminate between males of those subspecies on the basis of this protein molecule. Mouse salivary ABP might present a worthwhile system with which to study a prezygotic isolation mechanism in a mammal.

ABRUPT CLINE FOR SEX CHROMOSOMES IN A HYBRID ZONE BETWEEN TWO SPECIES OF MICE

We compared the patterns of movement of sex chromosomal and autosomal loci along a 160 km transect across a zone of hybridization between M. domesticus and M. musculus in southern Germany and western Austria using seven genetic markers. These included one Y-specific DNA sequence (YB10), two X-specific loci (DXWas68 and DXWas31), and four autosomal isozyme loci (Es-10, Es-1, Mpi-1, and Np-1). Random effects logistic regression analysis enabled us to examine the relationship between M. domesticus allele frequency and geographic distance from the western edge of the hybrid zone and allowed statistical evaluation of differences in cline midpoint and width among loci. More limited movement was observed for all three sex chromosomal markers across the zone compared with three of the four autosomal markers. If differential movement reflects fitness differences of specific alleles (or alleles at closely linked loci) on a hybrid background, then alleles that move to a limited extent across a hybrid zone may contribute to hybrid breakdown between two species. The limited flow of both X- and Y-specific alleles suggest that sex chromosomes have played an important role in Mus speciation.

Genome-wide mapping in a house mouse hybrid zone reveals hybrid sterility loci and Dobzhansky-Muller interactions

Mapping hybrid defects in contact zones between incipient species can identify genomic regions contributing to reproductive isolation and reveal genetic mechanisms of speciation. The house mouse features a rare combination of sophisticated genetic tools and natural hybrid zones between subspecies. Male hybrids often show reduced fertility, a common reproductive barrier between incipient species. Laboratory crosses have identified sterility loci, but each encompasses hundreds of genes. We map genetic determinants of testis weight and testis gene expression using offspring of mice captured in a hybrid zone between M. musculus musculus and M. m. domesticus. Many generations of admixture enables high-resolution mapping of loci contributing to these sterility-related phenotypes. We identify complex interactions among sterility loci, suggesting multiple, non-independent genetic incompatibilities contribute to barriers to gene flow in the hybrid zone.

Multilocus Characterization of a Woodrat (Genus Neotoma) Hybrid Zone

In order to investigate hybridization between 2 species of woodrats, Neotoma floridana and Neotoma micropus, 103 specimens were collected, in March of 1988, from a known area of sympatry, and compared with reference collections from areas of allopatry. Ten genetic markers, consisting of 7 microsatellite loci, 1 mitochondrial gene (cytochrome-b [Cytb]), and 2 nuclear introns (intron 2 of the vertebrate alcohol dehydrogenase gene [Adh1-I2] and intron 7 of the beta-fibrinogen gene [Fgb-I7]) were used to develop a composite genotype for each individual and for detection of hybridization. Six individuals were identified as pure parental N. micropus, 96 as hybrids, and 1 as pure parental N. floridana Hybrids were formed primarily through matings between complex genotypes, resulting in a high prevalence of individuals classified as backcrosses. The ratio of hybrid classes, population substructure, and presence of significant linkage disequilibrium within the zone of contact could not reject either the hybrid superiority or hybrid equilibrium model as responsible for maintenance of this hybrid zone. The collection date of this dataset (1988) provided not only a point in time assessment of the hybrid zone but also provided opportunities for future comparisons of temporal datasets with the purpose of examining hybrid zone characteristics over multiple generations.

Where are the wormy mice? A reexamination of hybrid parasitism in the European house mouse hybrid zone

Wormy mice in a hybrid zone have been interpreted as evidence of low hybrid fitness, such that parasites contribute to species separation. However, because of its natural heterogeneity, observations of parasite load must be numerous with good field area coverage. We sampled 689 mice from 107 localities across the Bavaria-Bohemia region of the European house mouse hybrid zone and calculated their hybrid indices using 1401 diagnostic single nucleotide polymorphisms (SNPs). We tested whether hybrids have greater or lesser diversity and load of parasite helminths than additive expectations, performing load analyses on the four most common taxa. We found hybrids have significantly reduced diversity and load of each of the commonest helminths; rarer helminths further support reduced load. Although within-locality comparisons have little power, randomization tests show the repeated pattern is unlikely to be due to local parasite heterogeneity, and simulations show a patch of low parasite diversity is unlikely to fall by chance just so in the field area, such that it produces the observed effects. Our data therefore contradict the idea that helminths reduce hybrid fitness through increased load. We discuss a vicariant Red Queen model that implies immune genes tracking parasites will escape Dobzhansky-Muller incompatibilities, generating hybrid variants untargeted by parasites.

Distinct lineages of Schistocephalus parasites in threespine and ninespine stickleback hosts revealed by DNA sequence analysis

Parasitic interactions are often part of complex networks of interspecific relationships that have evolved in biological communities. Despite many years of work on the evolution of parasitism, the likelihood that sister taxa of parasites can co-evolve with their hosts to specifically infect two related lineages, even when those hosts occur sympatrically, is still unclear. Furthermore, when these specific interactions occur, the molecular and physiological basis of this specificity is still largely unknown. The presence of these specific parasitic relationships can now be tested using molecular markers such as DNA sequence variation. Here we test for specific parasitic relationships in an emerging host-parasite model, the stickleback-Schistocephalus system. Threespine and ninespine stickleback fish are intermediate hosts for Schistocephalus cestode parasites that are phenotypically very similar and have nearly identical life cycles through plankton, stickleback, and avian hosts. We analyzed over 2000 base pairs of COX1 and NADH1 mitochondrial DNA sequences in 48 Schistocephalus individuals collected from threespine and ninespine stickleback hosts from disparate geographic regions distributed across the Northern Hemisphere. Our data strongly support the presence of two distinct clades of Schistocephalus, each of which exclusively infects either threespine or ninespine stickleback. These clades most likely represent different species that diverged soon after the speciation of their stickleback hosts. In addition, genetic structuring exists among Schistocephalus taken from threespine stickleback hosts from Alaska, Oregon and Wales, although it is much less than the divergence between hosts. Our findings emphasize that biological communities may be even more complex than they first appear, and beg the question of what are the ecological, physiological, and genetic factors that maintain the specificity of the Schistocephalus parasites and their stickleback hosts.

Evolutionary and functional properties of a two-locus beta-globin polymorphism in Indian house mice

Electrophoretic surveys of hemoglobin (Hb) polymorphism in house mice from South Asia and the Middle East have revealed that two alternative beta-globin haplotypes, Hbb(d) and Hbb(p), are often present at intermediate frequencies in geographically disparate populations. Both haplotypes harbor two functionally distinct beta-globin paralogs, HBB-T1 (which encodes the beta-chain subunits of the major Hb isoform) and HBB-T2 (which encodes the beta-chains of the minor Hb isoform). The Hbb(d) and Hbb(p) haplotypes share identical HBB-T1 alleles, but products of the alternative HBB-T2 alleles (d(minor) and p(minor)) are distinguished by two amino acid substitutions. To investigate the possible adaptive significance of the Hbb(d)/Hbb(p) polymorphism we conducted a population genetic analysis of the duplicated beta-globin genes of Indian house mice (Mus castaneus) in conjunction with experimental studies of Hb function in inbred strains of mice that carry the alternative Hbb(d) and Hbb(p) haplotypes. The main objectives of this study were (i) to characterize patterns of nucleotide polymorphism and linkage disequilibrium in the duplicated beta-globin genes of M. castaneus, (ii) to test the hypothesis that the Hbb(d) and Hbb(p) haplotypes are maintained as a balanced polymorphism, and (iii) to assess whether allelic differences in the alternative minor Hb isoforms (d(minor) and p(minor)) are associated with different O(2)-binding properties. A multilocus analysis of polymorphism and divergence revealed that levels of diversity at the HBB-T2 gene exceeded neutral expectations, and reconstructed haplotype networks for both beta-globin paralogs revealed extensive allele sharing with several other closely related species of Mus. However, despite this suggestive evidence for balancing selection, O(2)-equilibrium curves revealed no discernible functional differences between red cell lysates containing the d(minor) and p(minor) Hb isoforms. If the d(minor) and p(minor) alleles are maintained as a balanced polymorphism, our results indicate that the associated fitness variance is not directly related to respiratory functions of Hb.

Mosaic genealogy of the Mus musculus genome revealed by 21 nuclear genes from its three subspecies

Patterns of genetic variation provide insight into the evolutionary history of a species. Mouse (Mus musculus) is a good model for this purpose. Here we present the analysis of genealogies of the 21 nuclear loci and one mitochondrial DNA region in M. musculus based on our nucleotide sequences of nine inbred strains from three M. musculus subspecies (musculus, domesticus, and castaneus) and one M. spicilegus strain as an outgroup. The mitochondrial DNA gene genealogy of those strains confirmed the introgression pattern of one musculus strain. When all the nuclear DNA data were concatenated to produce a phylogenetic tree of nine strains, musculus and domesticus strains formed monophyletic clusters with each other, while the two castaneus strains were paraphyletic. When each DNA region was treated independently, the phylogenetic networks revealed an unnegligibly high level of subspecies admixture and the mosaic nature of their genome. Estimation of ancestral and derived population sizes and migration rates suggests the effects of ancestral polymorphism and gene flow on the pattern of genetic variation of the current subspecies. Gene genealogies of Fut4 and Dfy loci also suggested existence of the gene flow between M. musculus and M. spicilegus or other distant species.

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.

Complex signatures of selection and gene conversion in the duplicated globin genes of house mice

Results of electrophoretic surveys have suggested that hemoglobin polymorphism may be maintained by balancing selection in natural populations of house mice, Mus musculus. Here we report a survey of nucleotide variation in the adult globin genes of house mice from South America. We surveyed nucleotide polymorphism in two closely linked alpha-globin paralogs and two closely linked beta-globin paralogs to test whether patterns of variation are consistent with a model of long-term balancing selection. Surprisingly high levels of nucleotide polymorphism at the two beta-globin paralogs were attributable to the segregation of two highly divergent haplotypes, Hbbs (which carries two identical beta-globin paralogs) and Hbbd (which carries two functionally divergent beta-globin paralogs). Interparalog gene conversion on the Hbbs haplotype has produced a highly unusual situation in which the two paralogs are more similar to one another than either one is to its allelic counterpart on the Hbbd haplotype. Levels of nucleotide polymorphism and linkage disequilibrium at the two beta-globin paralogs suggest a complex history of diversity-enhancing selection that may be responsible for long-term maintenance of alternative protein alleles. The alternative two-locus beta-globin haplotypes are associated with pronounced differences in intraerythrocyte glutathione and nitric oxide metabolism, suggesting a possible mechanism for selection on hemoglobin function.

Region-specific rates of molecular evolution: a fourfold reduction in the rate of accumulation of "silent" mutations in transcribed versus nontranscribed regions of homologous DNA fragments derived from two closely related mouse species

We have sequenced homologous DNA fragments of 2.7 and 2.8 kbp derived from the closely related mouse species Mus musculus domesticus (M. domesticus) and Mus musculus musculus (M. musculus), respectively. These two species diverged approximately 1 million years ago. Each DNA fragment contains 1.35 kbp of the 3' end of the constitutively expressed 2.2-kbp aprt (adenine phosphoribosyltransferase) gene and a similarly sized nontranscribed region downstream of the aprt gene. The aprt gene region contains protein coding sequences (0.35 kbp), intronic sequences (0.75 kbp), and a 3' nontranslated sequence (0.25 kbp). Both the M. domesticus and M. musculus downstream regions share three partial copies of the B1 repetitive element with the M. musculus downstream region containing an additional complete copy of this element. A comparison of the 2.7- and 2.8-kbp DNA fragments revealed a total of 63 molecular alterations (i.e., mutations) that were approximately fourfold more abundant in the nontranscribed downstream region than in the transcribed aprt gene. Of the 11 mutations observed in the transcribed region, 7 were found in introns, 3 in the 3' untranslated sequence, and 1 was a synonymous change in an exon. A comparison of the human and M. domesticus aprt genes has previously revealed no homology in either the intronic or 3' nontranslated regions with the exception of a 26-bp sequence in intron 3 and sequences at the exon/intron boundaries necessary for correct mRNA splicing (Broderick et al., Proc. Natl. Acad. Sci. USA, 84:3349, 1987). Therefore, there does not appear to be selective pressure for sequences within these regions. We conclude that there is a lower rate of accumulation of "silent" mutations in the transcribed mouse aprt gene than in a contiguous nontranscribed downstream region. A possible molecular mechanism involving preferential DNA repair for the transcribed region is discussed.

Geographic origin of the Y Chromosomes in ?old? inbred strains of mice

Six distinct Y Chromosomes (Chr) were identified among 39 standard inbred strains of mice with five probes that identified Y Chr-specific restriction fragments on Southern blots. Three Y Chr types, distributed among 31 strains, were of Asian Mus musculus origin. The remaining three Y Chr types, distributed among eight strains, were of M. domesticus origin. The Asian source of the M. musculus Y Chr was confirmed by determining the DNA sequence of 221 bp from an open reading frame within the Sry (sex determining region Y) gene (Gubbay et al., Nature 346: 245-250, 1990) in three inbred strains (C57BL/6J, AKR/J, and SWR/J) and comparing the sequence to the homologous sequences derived from wild caught European and Asian M. musculus males. These data indicate that a minimum of six male mice contributed to the formation of the old inbred strains.

Sexual reproduction as an adaptation to resist parasites (a review)

Darwinian theory has yet to explain adequately the fact of sex. If males provide little or no aid to offspring, a high (up to 2-fold) extra average fitness has to emerge as a property of a sexual parentage if sex is to be stable. The advantage must presumably come from recombination but has been hard to identify. It may well lie in the necessity to recombine defenses to defeat numerous parasites. A model demonstrating this works best for contesting hosts whose defense polymorphisms are constrained to low mutation rates. A review of the literature shows that the predictions of parasite coevolution fit well with the known ecology of sex. Moreover, parasite coevolution is superior to previous models of the evolution of sex by supporting the stability of sex under the following challenging conditions: very low fecundity, realistic patterns of genotype fitness and changing environment, and frequent mutation to parthenogenesis, even while sex pays the full 2-fold cost.

Variations of a Y chromosome repeated sequence across subspecies of Mus musculus

The complex species Mus musculus is widespread in Eurasia and consists of four parapatric genetical entities (subspecies) that have recently radiated. Two of them (M. m. domesticus and M. m. musculus) are known to interact through a narrow zone of hybridisation across which autosomal and mitochondrial exchanges are very limited and Y chromosome exchange is absent. We extend here the study of this group by the genetical analysis of 22 Asian strains of various origins (China, Korea, Japan, Taiwan, Philippines and Indonesia). A survey of protein variation at ten polymorphic loci confirmed that these animals belong to either the subspecies M. m. musculus (northern type in Asia, ranging westwards to Eastern Europe) or to M. m. castaneus (southern Asian type) and revealed a certain degree of intergradation between the two taxa. Y chromosome variations were assessed in these strains using a Y specific DNA probe representing part of a small multigene family and also in four M. m. domesticus (the Western European house mouse) strains of various origins and one M. m. bactrianus (from Pakistan). Musculus and castaneus were identically monomorphic for one type of organisation of this Y repeated family, while domesticus and bactrianus were very similar to each other, showing slightly different types of organisation. Introgression of a bactrianus Y chromosome into the territory of castaneus was found in Indonesia. The present distribution of the Y types among the four subspecies is not phylogenetically concordant with the known distributions of autosomal and mitochondrial variants.(ABSTRACT TRUNCATED AT 250 WORDS)

Evolution of the immunoglobulin heavy chain variable region (Igh-V) locus in the genus Mus

The evolution of the mouse immunoglobulin heavy chain variable region (Igh-V) locus was investigated by the comprehensive analysis of variable region (Vh) gene family content and restriction fragment polymorphism in the genus Mus. The examination of natural Mus domesticus populations suggests an important role for recombination in the generation of the considerable restriction fragment polymorphism found at the Igh-V locus. Although the sizes of individual Vh gene families vary widely both within and between different Mus species, evolutionary trends of Vh gene family copy number are revealed by the analysis of homologues of mouse Vh gene families in Rattus and Peromyscus. Processes of duplication, deletion, and sequence divergence all contribute to the evolution of Vh gene copy number. Certain Vh gene families have expanded or contracted differently in the various muroid lineages examined. Collectively, these findings suggest that the evolution of individual Vh family size is not driven by strong selective pressure but is relatively neutral, and that gene flow, rather than selection, serves to maintain the high level of restriction fragment polymorphism seen in M. domesticus.

A unique alpha chain in hemoglobin of "Skive" Danish Mus musculus

The primary structures of the alpha chains in hemoglobins from three stocks of mice with the Hbaw2, Hbaw3, and Hbaw4 haplotypes were determined to establish whether the tentative alpha-chain assignments based on the results of isoelectric focusing patterns were correct. These Hba haplotypes were identified in laboratory descendants of feral mice captured in different parts of the world. Hemoglobin from "Centreville", Maryland, Mus musculus domesticus (Hbaw2) contains equal amounts of alpha chains 1 and 3. Hemoglobin from "Czech" Mus musculus musculus (Hbaw4) contains equal amounts of alpha chains 3 and 4. Amino acid analysis of the alpha-globins of "Skive" Danish Mus musculus musculus (Hbaw3) establishes that its hemoglobin is comprised of about one-third alpha chain 2 as expected plus a greater amount of a unique alpha chain that has not been described previously. This unique alpha chain has glycine at position 25, isoleucine at position 62, and serine at position 68; it is called chain 7. It may represent an intermediate in the evolution of genes that code for chain 2 (which has glycine, valine, and serine at positions 25, 62, and 68, respectively) and chain 4 (which has valine, isoleucine, and serine at positions 25, 62, and 62, respectively).

Diverse wild mouse origins of xenotropic, mink cell focus-forming, and two types of ecotropic proviral genes

We analyzed wild mouse DNAs for the number and type of proviral genes related to the env sequences of various murine leukemia viruses (MuLVs). Only Mus species closely related to laboratory mice carried these retroviral sequences, and the different subclasses of viral env genes tended to be restricted to specific taxonomic groups. Only Mus musculus molossinus carried proviral genes which cross-reacted with the inbred mouse ecotropic MuLV env gene. The ecotropic viral env sequence associated with the Fv-4 resistance gene was found in the Asian mice M. musculus molossinus and Mus musculus castaneus and in California mice from Lake Casitas (LC). Both M. musculus castaneus and LC mice carried many additional Fv-4 env-related proviruses, two of which are common to both mouse populations, which suggests that these mice share a recent common ancestry. Xenotropic and mink cell focus-forming (MCF) virus env sequences were more widely dispersed in wild mice than the ecotropic viral env genes, which suggests that nonecotropic MuLVs were integrated into the Mus germ line at an earlier date. Xenotropic MuLVs represented the major component of MuLV env-reactive genes in Asian and eastern European mice classified as M. musculus molossinus, M. musculus castaneus, and Mus musculus musculus, whereas Mus musculus domesticus from western Europe, the Mediterranean, and North America contained almost exclusively MCF virus env copies. M. musculus musculus mice from central Europe trapped near the M. musculus domesticus/M. musculus musculus hybrid zone carried multiple copies of both types of env genes. LC mice also carried both xenotropic and MCF viral env genes, which is consistent with the above conclusion that they represent natural hybrids of M. musculus domesticus and M. musculus castaneus.