Eurasian house mouse (Mus musculus L.) differentiation at microsatellite loci identifies the Iranian plateau as a phylogeographic hotspot

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Natural variation in the zinc-finger-encoding exon of Prdm9 affects hybrid sterility phenotypes in mice

PRDM9-mediated reproductive isolation was first described in the progeny of Mus musculus musculus (MUS) PWD/Ph and Mus musculus domesticus (DOM) C57BL/6J inbred strains. These male F1 hybrids fail to complete chromosome synapsis and arrest meiosis at prophase I, due to incompatibilities between the Prdm9 gene and hybrid sterility locus Hstx2. We identified 14 alleles of Prdm9 in exon 12, encoding the DNA-binding domain of the PRDM9 protein in outcrossed wild mouse populations from Europe, Asia, and the Middle East, 8 of which are novel. The same allele was found in all mice bearing introgressed t-haplotypes encompassing Prdm9. We asked whether 7 novel Prdm9 alleles in MUS populations and the t-haplotype allele in 1 MUS and 3 DOM populations induce Prdm9-mediated reproductive isolation. The results show that only combinations of the dom2 allele of DOM origin and the MUS msc1 allele ensure complete infertility of intersubspecific hybrids in outcrossed wild populations and inbred mouse strains examined so far. The results further indicate that MUS mice may share the erasure of PRDM9msc1 binding motifs in populations with different Prdm9 alleles, which implies that erased PRDM9 binding motifs may be uncoupled from their corresponding Prdm9 alleles at the population level. Our data corroborate the model of Prdm9-mediated hybrid sterility beyond inbred strains of mice and suggest that sterility alleles of Prdm9 may be rare.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Environmental implications of Ptolemaic Period rodents and shrews from the Sacred Falcon Necropolis at Quesna, Egypt (Mammalia: Muridae and Soricidae)

BACKGROUND

Assemblages of mummified and preserved animals in necropoleis of Ptolemaic Period Egypt (ca. 332-30 BC) document some aspects of the ceremonial and religious practices of the ancient Egyptians, but study of these animal remains can also provide insight into the local environments in which the animals and humans lived.

RESULTS

Excavations of the Sacred Falcon Necropolis at Quesna in the Nile Delta have yielded many thousands of animal remains, mostly of raptors, but also of a lesser number of small, wild mammals. Among the latter, we identified four species of murid rodents (Rodentia: Muridae) and five species of shrews (Eulipotyphla: Soricidae). The soricids are of particular interest because they represent a more diverse assemblage of species than occurs in the delta today. They include one species, Crocidura gueldenstaedtii (Pallas, 1811), that no longer occurs in the delta and another, C. fulvastra (Sundevall, 1843), that is now extirpated from Egypt.

CONCLUSIONS

The coexistence of this diverse small mammal community suggests that a greater availability and variety of mesic habitats were present during the Ptolemaic Period than occur there now. The local mammal faunas recovered at Quesna and other well-studied ancient Egyptian sites together provide evidence of a richer, more complex regional environment along the Nile Valley. They also provide important insight regarding the biogeography of the individual species comprising the faunas and about the extent of faunal turnover since the Ptolemaic Period.

Taxonomic assessment of two wild house mouse subspecies using whole-genome sequencing

The house mouse species complex (Mus musculus) is comprised of three primary subspecies. A large number of secondary subspecies have also been suggested on the basis of divergent morphology and molecular variation at limited numbers of markers. While the phylogenetic relationships among the primary M. musculus subspecies are well-defined, relationships among secondary subspecies and between secondary and primary subspecies remain less clear. Here, we integrate de novo genome sequencing of museum-stored specimens of house mice from one secondary subspecies (M. m. bactrianus) and publicly available genome sequences of house mice previously characterized as M. m. helgolandicus, with whole genome sequences from diverse representatives of the three primary house mouse subspecies. We show that mice assigned to the secondary M. m. bactrianus and M. m. helgolandicus subspecies are not genetically differentiated from M. m. castaneus and M. m. domesticus, respectively. Overall, our work suggests that the M. m. bactrianus and M. m. helgolandicus subspecies are not well-justified taxonomic entities, emphasizing the importance of leveraging whole-genome sequence data to inform subspecies designations. Additionally, our investigation provides tailored experimental procedures for generating whole genome sequences from air-dried mouse skins, along with key genomic resources to inform future genomic studies of wild mouse diversity.

Population structure and inbreeding in wild house mice (Mus musculus) at different geographic scales

House mice (Mus musculus) have spread globally as a result of their commensal relationship with humans. In the form of laboratory strains, both inbred and outbred, they are also among the most widely used model organisms in biomedical research. Although the general outlines of house mouse dispersal and population structure are well known, details have been obscured by either limited sample size or small numbers of markers. Here we examine ancestry, population structure, and inbreeding using SNP microarray genotypes in a cohort of 814 wild mice spanning five continents and all major subspecies of Mus, with a focus on M. m. domesticus. We find that the major axis of genetic variation in M. m. domesticus is a south-to-north gradient within Europe and the Mediterranean. The dominant ancestry component in North America, Australia, New Zealand, and various small offshore islands are of northern European origin. Next we show that inbreeding is surprisingly pervasive and highly variable, even between nearby populations. By inspecting the length distribution of homozygous segments in individual genomes, we find that inbreeding in commensal populations is mostly due to consanguinity. Our results offer new insight into the natural history of an important model organism for medicine and evolutionary biology.

Phylogeography of the White-Toothed Shrews Crocidura suaveolens and Crocidura sibirica: Searching for the Geographical Homeland

Abstract

We studied the polymorphism of the cytb gene in two forms of the Lesser White-toothed Shrew species complex: Crocidura suaveolens s. stricto and C. sibirica. The haplotypes of C. sibirica are found to be very similar to those of Crocidura suaveolens. They do not belong to a distinct haplogroup. The molecular diversity of the populations in the Asian part of the range is higher than in Eastern Europe. For the combined sample from Asia and Europe together, we revealed a significant signal of population expansion. Analysis of the expansion time showed that the Asian territory was colonized earlier (before the last glacial maximum) than the Eastern Europe (at the very end of the Late Pleistocene and in the early Holocene). The results of the ancestral area reconstruction are consistent with the hypothesis of a Middle Asian origin of the C. suaveolens/C. sibirica group, recent colonization of Inner Asia and later penetration into Eastern Europe.

Insights into Mus musculus population structure across Eurasia revealed by whole-genome analysis

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000–500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000–226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.

Bidirectional introgression between Mus musculus domesticus and Mus spretus

Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles—including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.

Tracing the origin and evolutionary fate of recent gene retrocopies in natural populations of the house mouse

Although the contribution of retrogenes to the evolution of genes and genomes has long been recognized, the evolutionary patterns of very recently derived retrocopies that are still polymorphic within natural populations have not been much studied so far. We use here a set of 2,025 such retrocopies in nine house mouse populations from three subspecies (Mus musculus domesticus, M. m. musculus, and M. m. castaneus) to trace their origin and evolutionary fate. We find that ancient house-keeping genes are significantly more likely to generate retrocopies than younger genes and that the propensity to generate a retrocopy depends on its level of expression in the germline. Although most retrocopies are detrimental and quickly purged, we focus here on the subset that appears to be neutral or even adaptive. We show that retrocopies from X-chromosomal parental genes have a higher likelihood to reach elevated frequencies in the populations, confirming the notion of adaptive effects for “out-of-X” retrogenes. Also, retrocopies in intergenic regions are more likely to reach higher population frequencies than those in introns of genes, implying a more detrimental effect when they land within transcribed regions. For a small subset of retrocopies, we find signatures of positive selection, indicating they were involved in a recent adaptation process. We show that the population-specific distribution pattern of retrocopies is phylogenetically informative and can be used to infer population history with a better resolution than with SNP markers.

Selection shapes the landscape of functional variation in wild house mice

BACKGROUND

Through human-aided dispersal over the last ~ 10,000 years, house mice (Mus musculus) have recently colonized diverse habitats across the globe, promoting the emergence of new traits that confer adaptive advantages in distinct environments. Despite their status as the premier mammalian model system, the impact of this demographic and selective history on the global patterning of disease-relevant trait variation in wild mouse populations is poorly understood.

RESULTS

Here, we leveraged 154 whole-genome sequences from diverse wild house mouse populations to survey the geographic organization of functional variation and systematically identify signals of positive selection. We show that a significant proportion of wild mouse variation is private to single populations, including numerous predicted functional alleles. In addition, we report strong signals of positive selection at many genes associated with both complex and Mendelian diseases in humans. Notably, we detect a significant excess of selection signals at disease-associated genes relative to null expectations, pointing to the important role of adaptation in shaping the landscape of functional variation in wild mouse populations. We also uncover strong signals of selection at multiple genes involved in starch digestion, including Mgam and Amy1. We speculate that the successful emergence of the human-mouse commensalism may have been facilitated, in part, by dietary adaptations at these loci. Finally, our work uncovers multiple cryptic structural variants that manifest as putative signals of positive selection, highlighting an important and under-appreciated source of false-positive signals in genome-wide selection scans.

CONCLUSIONS

Overall, our findings highlight the role of adaptation in shaping wild mouse genetic variation at human disease-associated genes. Our work also highlights the biomedical relevance of wild mouse genetic diversity and underscores the potential for targeted sampling of mice from specific populations as a strategy for developing effective new mouse models of both rare and common human diseases.

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.

Non-linear genetic diversity and notable population differentiation caused by low gene flow of bermudagrass [Cynodon dactylon (L.) Pers.] along longitude gradients

Background

Environmental variation related to ecological habitat is the main driver of plant adaptive divergence. Longitude plays an important role in the formation of plant population structure, indicating that environmental differentiation can significantly shape population structure.

Methods

Genetic diversity and population genetic structure were estimated using 105 expressed sequence tag-derived simple sequence repeat (EST-SSR) loci. A total of 249 C. dactylon (L.) Pers. (common bermudagrass) individuals were sampled from 13 geographic sites along the longitude (105°57′34″–119°27′06″E).

Results

There was no obvious linear trend of intra-population genetic diversity along longitude and the intra-population genetic diversity was not related to climate in this study. Low gene flow (Nm = 0.7701) meant a rich genetic differentiation among populations of C. dactylon along longitude gradients. Significantly positive Mantel correlation (r = 0.438, P = 0.001) was found between genetic distance and geographical interval while no significant partial Mantel correlation after controlling the effect of mean annual precipitation, which indicated geographic distance correlated with mean annual precipitation affect genetic distance. The genetic diversity of C. dactylon with higher ploidy level was higher than that with lower ploidy level and groups of individuals with higher ploidy level were separated further away by genetic distance from the lower ploidy levels. Understanding the different genetic bases of local adaptation comparatively between latitude and longitude is one of the core findings in the adaptive evolution of plants.

Forest Saccharomyces paradoxus are robust to seasonal biotic and abiotic changes

Microorganisms are famous for adapting quickly to new environments. However, most evidence for rapid microbial adaptation comes from laboratory experiments or domesticated environments, and it is unclear how rates of adaptation scale from human-influenced environments to the great diversity of wild microorganisms. We examined potential monthly-scale selective pressures in the model forest yeast Saccharomyces paradoxus. Contrary to expectations of seasonal adaptation, the S. paradoxus population was stable over four seasons in the face of abiotic and biotic environmental changes. While the S. paradoxus population was diverse, including 41 unique genotypes among 192 sampled isolates, there was no correlation between S. paradoxus genotypes and seasonal environments. Consistent with observations from other S. paradoxus populations, the forest population was highly clonal and inbred. This lack of recombination, paired with population stability, implies that selection is not acting on the forest S. paradoxus population on a seasonal timescale. Saccharomyces paradoxus may instead have evolved generalism or phenotypic plasticity with regard to seasonal environmental changes long ago. Similarly, while the forest population included diversity among phenotypes related to intraspecific interference competition, there was no evidence for active coevolution among these phenotypes. At least ten percent of the forest S. paradoxus individuals produced "killer toxins," which kill sensitive Saccharomyces cells, but the presence of a toxin-producing isolate did not predict resistance to the toxin among nearby isolates. How forest yeasts acclimate to changing environments remains an open question, and future studies should investigate the physiological responses that allow microbial cells to cope with environmental fluctuations in their native habitats.

The mutational load in natural populations is significantly affected by high primary rates of retroposition

The phenomenon of retroposition (the reintegration of reverse-transcribed RNA into the genome) has been well studied in comparisons between species and has been identified as a source of evolutionary innovation. However, less attention has been paid to possible negative effects of retroposition. To trace the evolutionary dynamics of these negative effects, our study uses a unique genomic dataset of house mouse populations. It reveals that the initial retroposition rate is very high and that most of these newly transposed retrocopies have a deleterious impact, apparently through modifying the expression of their parental genes. In humans, this effect is expected to cause disease alleles, and we propose that genetic screening should include the search for newly transposed retrocopies.

Gene retroposition is known to contribute to patterns of gene evolution and adaptations. However, possible negative effects of gene retroposition remain largely unexplored since most previous studies have focused on between-species comparisons where negatively selected copies are mostly not observed, as they are quickly lost from populations. Here, we show for natural house mouse populations that the primary rate of retroposition is orders of magnitude higher than the long-term rate. Comparisons with single-nucleotide polymorphism distribution patterns in the same populations show that most retroposition events are deleterious. Transcriptomic profiling analysis shows that new retroposed copies become easily subject to transcription and have an influence on the expression levels of their parental genes, especially when transcribed in the antisense direction. Our results imply that the impact of retroposition on the mutational load has been highly underestimated in natural populations. This has additional implications for strategies of disease allele detection in humans.

Ancient mitochondrial DNA connects house mice in the British Isles to trade across Europe over three millennia

Background

The earliest records in Britain for the western European house mouse (Mus musculus domesticus) date from the Late Bronze Age. The arrival of this commensal species in Britain is thought to be related to human transport and trade with continental Europe. In order to study this arrival, we collected a total of 16 ancient mouse mandibulae from four early British archaeological sites, ranging from the Late Bronze Age to the Roman period.

Results

From these, we obtained ancient mitochondrial DNA (mtDNA) house mouse sequences from eight house mice from two of the sites dating from the Late Bronze to Middle Iron Age. We also obtained five ancient mtDNA wood mouse (Apodemus spp.) sequences from all four sites. The ancient house mouse sequences found in this study were from haplogroups E (N = 6) and D (N = 2). Modern British house mouse mtDNA sequences are primarily characterised by haplogroups E and F and, much less commonly, haplogroup D.

Conclusions

The presence of haplogroups D and E in our samples and the dating of the archaeological sites provide evidence of an early house mouse colonisation that may relate to Late Bronze Age/Iron Age trade and/or human expansion. Our results confirm the hypothesis, based on zooarchaeological evidence and modern mtDNA predictions, that house mice, with haplogroups D and E, were established in Britain by the Iron Age and, in the case of haplogroup E, possibly as early as the Late Bronze Age.

Tracking the Near Eastern origins and European dispersal of the western house mouse

The house mouse (Mus musculus) represents the extreme of globalization of invasive mammals. However, the timing and basis of its origin and early phases of dispersal remain poorly documented. To track its synanthropisation and subsequent invasive spread during the develoment of complex human societies, we analyzed 829 Mus specimens from 43 archaeological contexts in Southwestern Asia and Southeastern Europe, between 40,000 and 3,000 cal. BP, combining geometric morphometrics numerical taxonomy, ancient mitochondrial DNA and direct radiocarbon dating. We found that large late hunter-gatherer sedentary settlements in the Levant, c. 14,500 cal. BP, promoted the commensal behaviour of the house mouse, which probably led the commensal pathway to cat domestication. House mouse invasive spread was then fostered through the emergence of agriculture throughout the Near East 12,000 years ago. Stowaway transport of house mice to Cyprus can be inferred as early as 10,800 years ago. However, the house mouse invasion of Europe did not happen until the development of proto urbanism and exchange networks - 6,500 years ago in Eastern Europe and 4000 years ago in Southern Europe - which in turn may have driven the first human mediated dispersal of cats in Europe.

The amylase gene cluster in house mice (Mus musculus) was subject to repeated introgression including the rescue of a pseudogene

Background

Amylase gene clusters have been implicated in adaptive copy number changes in response to the amount of starch in the diet of humans and mammals. However, this interpretation has been questioned for humans and for mammals there is a paucity of information from natural populations.

Results

Using optical mapping and genome read information, we show here that the amylase cluster in natural house mouse populations is indeed copy-number variable for Amy2b paralogous gene copies (called Amy2a1 - Amy2a5), but a direct connection to starch diet is not evident. However, we find that the amylase cluster was subject to introgression of haplotypes between Mus musculus sub-species. A very recent introgression can be traced in the Western European populations and this leads also to the rescue of an Amy2b pseudogene. Some populations and inbred lines derived from the Western house mouse (Mus musculus domesticus) harbor a copy of the pancreatic amylase (Amy2b) with a stop codon in the first exon, making it non-functional. But populations in France harbor a haplotype introgressed from the Eastern house mouse (M. m. musculus) with an intact reading frame. Detailed analysis of phylogenetic patterns along the amylase cluster suggest an additional history of previous introgressions.

Conclusions

Our results show that the amylase gene cluster is a hotspot of introgression in the mouse genome, making it an evolutionary active region beyond the previously observed copy number changes.

Population Genomics of the House Mouse and the Brown Rat

Mice (Mus musculus) and rats (Rattus norvegicus) have long served as model systems for biomedical research. However, they are also excellent models for studying the evolution of populations, subspecies, and species. Within the past million years, they have spread in various waves across large parts of the globe, with the most recent spread in the wake of human civilization. They have developed into commensal species, but have also been able to colonize extreme environments on islands free of human civilization. Given that ample genomic and genetic resources are available for these species, they have thus also become ideal mammalian systems for evolutionary studies on adaptation and speciation, particularly in the combination with the rapid developments in population genomics. The chapter provides an overview of the systems and their history, as well as of available resources.

Quantifying the efficiency and biases of forest Saccharomyces sampling strategies

Saccharomyces yeasts are emerging as model organisms for ecology and evolution, and researchers need environmental Saccharomyces isolates to test ecological and evolutionary hypotheses. However, methods for isolating Saccharomyces from nature have not been standardized, and isolation methods may influence the genotypes and phenotypes of studied strains. We compared the effectiveness and potential biases of an established enrichment culturing method against a newly developed direct plating method for isolating forest floor Saccharomyces spp. In a European forest, enrichment culturing was both less successful at isolating Saccharomyces paradoxus per sample collected and less labour intensive per isolated S. paradoxus colony than direct isolation. The two methods sampled similar S. paradoxus diversity: The number of unique genotypes sampled (i.e., genotypic diversity) per S. paradoxus isolate and average growth rates of S. paradoxus isolates did not differ between the two methods, and growth rate variances (i.e., phenotypic diversity) only differed in one of three tested environments. However, enrichment culturing did detect rare Saccharomyces cerevisiae in the forest habitat and also found two S. paradoxus isolates with outlier phenotypes. Our results validate the historically common method of using enrichment culturing to isolate representative collections of environmental Saccharomyces. We recommend that researchers choose a Saccharomyces sampling method based on resources available for sampling and isolate screening. Researchers interested in discovering new Saccharomyces phenotypes or rare Saccharomyces species from natural environments may also have more success using enrichment culturing. We include step-by-step sampling protocols in the supplemental materials.

Dealing with the adaptive immune system during de novo evolution of genes from intergenic sequences

Background

The adaptive immune system of vertebrates has an extraordinary potential to sense and neutralize foreign antigens entering the body. De novo evolution of genes implies that the genome itself expresses novel antigens from intergenic sequences which could cause a problem with this immune system. Peptides from these novel proteins could be presented by the major histocompatibility complex (MHC) receptors to the cell surface and would be recognized as foreign. The respective cells would then be attacked and destroyed, or would cause inflammatory responses. Hence, de novo expressed peptides have to be introduced to the immune system as being self-peptides to avoid such autoimmune reactions. The regulation of the distinction between self and non-self starts during embryonic development, but continues late into adulthood. It is mostly mediated by specialized cells in the thymus, but can also be conveyed in peripheral tissues, such as the lymph nodes and the spleen. The self-antigens need to be exposed to the reactive T-cells, which requires the expression of the genes in the respective tissues. Since the initial activation of a promotor for new intergenic transcription of a de novo gene could occur in any tissue, we should expect that the evolutionary establishment of a de novo gene in animals with an adaptive immune system should also involve expression in at least one of the tissues that confer self-recognition.

Results

We have studied this question by analyzing the transcriptomes of multiple tissues from young mice in three closely related natural populations of the house mouse (M. m. domesticus). We find that new intergenic transcription occurs indeed mostly in only a single tissue. When a second tissue becomes involved, thymus and spleen are significantly overrepresented.

Conclusions

We conclude that the inclusion of de novo transcripts in the processes for the induction of self-tolerance is indeed an important step in the evolution of functional de novo genes in vertebrates.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1232-z) contains supplementary material, which is available to authorized users.

Meta-populational demes constitute a reservoir for large MHC allele diversity in wild house mice (Mus musculus)

Background

The MHC class I and II loci mediate the adaptive immune response and belong to the most polymorphic loci in vertebrate genomes. In fact, the number of different alleles in a given species is often so large that it remains a challenge to provide an evolutionary model that can fully account for this.

Results

We provide here a general survey of MHC allele numbers in house mouse populations and two sub-species (M. m. domesticus and M. m. musculus) for H2 class I D and K, as well as class II A and E loci. Between 50 and 90% of the detected different sequences constitute new alleles, confirming that the discovery of new alleles is indeed far from complete. House mice live in separate demes with small effective population sizes, factors that were proposed to reduce, rather than enhance the possibility for the maintenance of many different alleles. To specifically investigate the occurrence of alleles within demes, we focused on the class II H2-Aa and H2-Eb exon 2 alleles in nine demes of M. m. domesticus from two different geographic regions. We find on the one hand a group of alleles that occur in different sampling regions and three quarters of these are also found in both sub-species. On the other hand, the larger group of different alleles (56%) occurs only in one of the regions and most of these (89%) only in single demes. We show that most of these region-specific alleles have apparently arisen through recombination and/or partial gene conversion from already existing alleles.

Conclusions

Demes can act as sources of alleles that outnumber the set of alleles that are shared across the species range. These findings support the reservoir model proposed for human MHC diversity, which states that large pools of rare MHC allele variants are continuously generated by neutral mutational mechanisms. Given that these can become important in the defense against newly emerging pathogens, the reservoir model complements the selection based models for MHC diversity and explains why the exceptional diversity exists.

Electronic supplementary material

The online version of this article (10.1186/s12983-018-0266-9) contains supplementary material, which is available to authorized users.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Mitochondrial Genetic Diversity and Phylogeography of Mus musculus castaneus in Northern Punjab, Pakistan

Regions of Iran, Afghanistan, Pakistan and northwestern India have been proposed as the place of origin of Mus musculus castaneus. But despite the fact that Pakistan encompasses an important part of its range, M. m. castaneus populations in Pakistan have not been the subject of intensive genetic and biogeographic studies, except for a very small number of samples included in past studies. We studied genetic variation in M. m. castaneus (CAS) from northern Punjab Province, Pakistan, by using cytochrome b (Cytb) analysis in a sample of 98 individuals. Median-joining network revealed four well differentiated CAS sub-lineages coexisting within a small geographical region; these had previously been thought to have largely non-overlapping geographic distributions. Moreover, haplotypes from Pakistan occupied a central position in the network and all identified global haplotypes were also present in Pakistan. All identified CAS sub-lineages proved to be highly diverse on the basis of haplotype and nucleotide diversity indices. Tajima's D test and Fu's Fs tests of neutrality suggest recent population expansions in all sub-lineages. Expansion times were estimated as 21,760-134,930, 10,800-64,400 and 4950-30,665 ybp using substitution rates of 2.5%, 5% and 11%, respectively. Our results support the hypothesis that northern Punjab Province in Pakistan is the most likely source area for M. m. castaneus, and that the CAS sub-lineages in this region have undergone rapid population expansion events at different time periods, which appear to have benefitted from human-mediated transport, although one of them clearly predates the establishment of human settlements in this region.

Whole exome sequencing of wild-derived inbred strains of mice improves power to link phenotype and genotype

The house mouse is a powerful model to dissect the genetic basis of phenotypic variation, and serves as a model to study human diseases. Despite a wealth of discoveries, most classical laboratory strains have captured only a small fraction of genetic variation known to segregate in their wild progenitors, and existing strains are often related to each other in complex ways. Inbred strains of mice independently derived from natural populations have the potential to increase power in genetic studies with the addition of novel genetic variation. Here, we perform exome-enrichment and high-throughput sequencing (~8× coverage) of 26 wild-derived strains known in the mouse research community as the "Montpellier strains." We identified 1.46 million SNPs in our dataset, approximately 19% of which have not been detected from other inbred strains. This novel genetic variation is expected to contribute to phenotypic variation, as they include 18,496 nonsynonymous variants and 262 early stop codons. Simulations demonstrate that the higher density of genetic variation in the Montpellier strains provides increased power for quantitative genetic studies. Inasmuch as the power to connect genotype to phenotype depends on genetic variation, it is important to incorporate these additional genetic strains into future research programs.

First molecular evidence of Mus musculus bactrianus in Nepal inferred from the mitochondrial DNA cytochrome B gene sequences

To identify the house mice collected in Pokhara and Lumbini of Nepal at the subspecies level, morphological and molecular analyses were carried out. Morphologically, two populations collected in Pokhara and Lumbini were distinguished by fur colour, but there was no significant difference in external measurements (p > .05). The phylogenetic analysis results revealed that the haplotypes sequences of mitochondrial DNA (mtDNA) Cytochrome B (CytB) gene distinguished into two distinct clades on a phylogenetic tree representing two subspecies, Mus musculus bactrianus and M. m. castaneus in Pokhara and Lumbini, respectively. In Nepal, the subspecies M. m. bactrianus was not reported before this study. These findings concluded that at least two subspecies, M. m. bactrianus and M. m. castaneus currently exist in Nepal. We estimated that these two subspecies could have introduced together with human migration, while further study is required to understand their evolutionary history and current distribution.

Phylogeny and adaptation shape the teeth of insular mice

By accompanying human travels since prehistorical times, the house mouse dispersed widely throughout the world, and colonized many islands. The origin of the travellers determined the phylogenetic source of the insular mice, which encountered diverse ecological and environmental conditions on the various islands. Insular mice are thus an exceptional model to disentangle the relative role of phylogeny, ecology and climate in evolution. Molar shape is known to vary according to phylogeny and to respond to adaptation. Using for the first time a three-dimensional geometric morphometric approach, compared with a classical two-dimensional quantification, the relative effects of size variation, phylogeny, climate and ecology were investigated on molar shape diversity across a variety of islands. Phylogeny emerged as the factor of prime importance in shaping the molar. Changes in competition level, mostly driven by the presence or absence of the wood mouse on the different islands, appeared as the second most important effect. Climate and size differences accounted for slight shape variation. This evidences a balanced role of random differentiation related to history of colonization, and of adaptation possibly related to resource exploitation.

Genomic resources for wild populations of the house mouse, Mus musculus and its close relative Mus spretus

Wild populations of the house mouse (Mus musculus) represent the raw genetic material for the classical inbred strains in biomedical research and are a major model system for evolutionary biology. We provide whole genome sequencing data of individuals representing natural populations of M. m. domesticus (24 individuals from 3 populations), M. m. helgolandicus (3 individuals), M. m. musculus (22 individuals from 3 populations) and M. spretus (8 individuals from one population). We use a single pipeline to map and call variants for these individuals and also include 10 additional individuals of M. m. castaneus for which genomic data are publically available. In addition, RNAseq data were obtained from 10 tissues of up to eight adult individuals from each of the three M. m. domesticus populations for which genomic data were collected. Data and analyses are presented via tracks viewable in the UCSC or IGV genome browsers. We also provide information on available outbred stocks and instructions on how to keep them in the laboratory.

Fast turnover of genome transcription across evolutionary time exposes entire non-coding DNA to de novo gene emergence

Deep sequencing analyses have shown that a large fraction of genomes is transcribed, but the significance of this transcription is much debated. Here, we characterize the phylogenetic turnover of poly-adenylated transcripts in a comprehensive sampling of taxa of the mouse (genus Mus), spanning a phylogenetic distance of 10 Myr. Using deep RNA sequencing we find that at a given sequencing depth transcriptome coverage becomes saturated within a taxon, but keeps extending when compared between taxa, even at this very shallow phylogenetic level. Our data show a high turnover of transcriptional states between taxa and that no major transcript-free islands exist across evolutionary time. This suggests that the entire genome can be transcribed into poly-adenylated RNA when viewed at an evolutionary time scale. We conclude that any part of the non-coding genome can potentially become subject to evolutionary functionalization via de novo gene evolution within relatively short evolutionary time spans.

DOI:http://dx.doi.org/10.7554/eLife.09977.001

Traditionally, the genome – the sum total of DNA within a cell – was thought to be divided into genes and ‘non-coding’ regions. Genes are copied, or “transcribed”, into molecules called RNA that perform essential tasks in the cell. The roles of the non-coding regions were often less clear, although it has since become apparent that some are also transcribed and generate low levels of RNA molecules. However, many debate how significant this transcription is to living organisms.

Neme and Tautz have now used a technique called deep RNA sequencing to study the RNA molecules produced in several different species and types of mice whose last common ancestor lived 10 million years ago. Different species produced RNA molecules from different portions – both genes and non-coding regions – of their genomes. Comparing these RNA sequences suggests that changes to the regions that are transcribed occur relatively quickly for a large portion of the genome. Furthermore, there have been no significant areas of the common ancestor’s genome that have not been transcribed at some point in at least one of its descendent species.

This therefore suggests that over a relatively short evolutionary period, any part of the genome can acquire the ability to be transcribed and potentially form a new gene. The next challenge is to find out how often these transcribed non-coding parts of the genome show important biochemical activities, and how they find their way into becoming new genes.

DOI:http://dx.doi.org/10.7554/eLife.09977.002

Selective sweeps versus introgression - population genetic dynamics of the murine leukemia virus receptor Xpr1 in wild populations of the house mouse (Mus musculus)

Background

The interaction between viruses and their receptors in the host can be expected to lead to an evolutionary arms race resulting in cycles of rapid adaptations. We focus here on the receptor gene Xpr1 (xenotropic and polytropic retrovirus receptor 1) for murine leukemia viruses (MLVs). In a previous screen for selective sweeps in mouse populations we discovered that a population from Germany was almost monomorphic for Xpr1 haplotypes, while a population from France was polymorphic.

Results

Here we analyze Xpr1 sequences and haplotypes from a broad sample of wild mouse populations of two subspecies, M. m. domesticus and M. m. musculus, to trace the origins of this distinctive polymorphism pattern. We show that the high polymorphism in the population in France is caused by a relatively recent invasion of a haplotype from a population in Iran, rather than a selective sweep in Germany. The invading haplotype codes for a novel receptor variant, which has itself undergone a recent selective sweep in the Iranian population.

Conclusions

Our data support a scenario in which Xpr1 is frequently subject to positive selection, possibly as a response to resistance development against recurrently emerging infectious viruses. During such an infection cycle, receptor variants that may convey viral resistance can be captured from another population and quickly introgress into populations actively dealing with the infectious virus.

Electronic supplementary material

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

Divergence patterns of genic copy number variation in natural populations of the house mouse (Mus musculus domesticus) reveal three conserved genes with major population-specific expansions

Copy number variation represents a major source of genetic divergence, yet the evolutionary dynamics of genic copy number variation in natural populations during differentiation and adaptation remain unclear. We applied a read depth approach to genome resequencing data to detect copy number variants (CNVs) ≥1 kb in wild-caught mice belonging to four populations of Mus musculus domesticus. We complemented the bioinformatics analyses with experimental validation using droplet digital PCR. The specific focus of our analysis is CNVs that include complete genes, as these CNVs could be expected to contribute most directly to evolutionary divergence. In total, 1863 transcription units appear to be completely encompassed within CNVs in at least one individual when compared to the reference assembly. Further, 179 of these CNVs show population-specific copy number differences, and 325 are subject to complete deletion in multiple individuals. Among the most copy-number variable genes are three highly conserved genes that encode the splicing factor CWC22, the spindle protein SFI1, and the Holliday junction recognition protein HJURP. These genes exhibit population-specific expansion patterns that suggest involvement in local adaptations. We found that genes that overlap with large segmental duplications are generally more copy-number variable. These genes encode proteins that are relevant for environmental and behavioral interactions, such as vomeronasal and olfactory receptors, as well as major urinary proteins and several proteins of unknown function. The overall analysis shows that genic CNVs contribute more to population differentiation in mice than in humans and may promote and speed up population divergence.

Insights into mammalian biology from the wild house mouse Mus musculus

The house mouse, Mus musculus, was established in the early 1900s as one of the first genetic model organisms owing to its short generation time, comparatively large litters, ease of husbandry, and visible phenotypic variants. For these reasons and because they are mammals, house mice are well suited to serve as models for human phenotypes and disease. House mice in the wild consist of at least three distinct subspecies and harbor extensive genetic and phenotypic variation both within and between these subspecies. Wild mice have been used to study a wide range of biological processes, including immunity, cancer, male sterility, adaptive evolution, and non-Mendelian inheritance. Despite the extensive variation that exists among wild mice, classical laboratory strains are derived from a limited set of founders and thus contain only a small subset of this variation. Continued efforts to study wild house mice and to create new inbred strains from wild populations have the potential to strengthen house mice as a model system.

DOI:http://dx.doi.org/10.7554/eLife.05959.001