Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Gene expression plasticity and desert adaptation in house mice

Understanding how organisms adapt to new environments is a key problem in evolution, yet it remains unclear whether phenotypic plasticity generally facilitates or hinders this process. Here we studied evolved and plastic responses to water-stress in lab-born descendants of wild house mice (Mus musculus domesticus) collected from desert and non-desert environments and measured gene expression and organismal phenotypes under control and water-stressed conditions. After many generations in the lab, desert mice consumed significantly less water than mice from other localities, indicating that this difference has a genetic basis. Under water-stress, desert mice maintained more weight than non-desert mice, and exhibited differences in blood chemistry related to osmoregulatory function. Gene expression in the kidney revealed evolved differences between mice from different environments as well as plastic responses between hydrated and dehydrated mice. Desert mice showed reduced expression plasticity under water-stress compared to non-desert mice. Importantly, non-desert mice under water-stress generally showed shifts toward desert-like expression, consistent with adaptive plasticity. Finally, we identify several co-expression modules linked to phenotypes of interest. These findings provide evidence for local adaptation after a recent invasion and suggest that adaptive plasticity may have facilitated colonization of the desert environment.

Haldane's rule in the placenta: Sex-biased misregulation of the Kcnq1 imprinting cluster in hybrid mice

Hybrid phenotypes that contribute to postzygotic reproductive isolation often exhibit pronounced asymmetry, both between reciprocal crosses and between the sexes in accordance with Haldane's rule. Inviability in mammalian hybrids is associated with parent-of-origin placental growth abnormalities for which misregulation of imprinted gene (IGs) is the leading candidate mechanism. However, direct evidence for the involvement of IGs in hybrid growth dysplasia is limited. We used transcriptome and reduced representation bisulfite sequencing to conduct the first genome-scale assessment of the contribution of IGs to parent-of-origin placental growth dysplasia in the cross between the house mouse (Mus musculus domesticus) and the Algerian mouse (Mus spretus). IGs with transgressive expression and methylation were concentrated in the Kcnq1 cluster, which contains causal genes for prenatal growth abnormalities in mice and humans. Hypermethylation of the cluster's imprinting control region, and consequent misexpression of the genes Phlda2 and Ascl2, is a strong candidate mechanism for transgressive placental undergrowth. Transgressive placental and gene regulatory phenotypes, including expression and methylation in the Kcnq1 cluster, were more extreme in hybrid males. Although consistent with Haldane's rule, male-biased defects are unexpected in rodent placenta because the X-chromosome is effectively hemizygous in both sexes. In search of an explanation, we found evidence of leaky imprinted (paternal) X-chromosome inactivation in hybrid female placenta, an epigenetic disturbance that may buffer females from the effects of X-linked incompatibilities to which males are fully exposed. Sex differences in chromatin structure on the X and sex-biased maternal effects are nonmutually exclusive alternative explanations for adherence to Haldane's rule in hybrid placenta. The results of this study contribute to understanding the genetic basis of hybrid inviability in mammals, and the role of IGs in speciation.

Zoonotic Gastrointestinal Helminths in Rodent Communities in Southern Guatemala

Rodents are reservoirs and hosts of several pathogens around the world, including zoonotic parasite species. This study aimed to determine the occurrence of zoonotic gastrointestinal helminths in rodents captured inside households in a rural community from southern Guatemala. Sixty-nine rodents were captured in 33% (49/148) of the surveyed households, including Rattus rattus, Rattus norvegicus, Mus musculus, and Sigmodon hispidus. Thirty-six percent (25/69) of these rodents (3 Rattus and 22 Mus musculus), from 45% (22/49) of the households, were parasitized with at least 1 gastrointestinal helminth species. Helminths from 6 species were identified: Hymenolepis diminuta, Moniliformis moniliformis, Heterakis spumosa, Nippostrongylus sp., Strongyloides sp., and Syphacia sp. Two zoonotic species were found in Rattus, H. diminuta in R. rattus (1/6), and M. moniliformis in R. norvegicus (1/1). Coinfection with other non-zoonotic helminth parasites, such as He. spumosa and Strongyloides sp., also was observed in the Rattus genus. Mus musculus had only non-zoonotic helminths: He. spumosa, Nippostrongylus sp., and Syphacia sp. being the most common, and He. spumosa (96%) followed by Nippostrongylus sp. (48%), with a higher presence in males than females, with a similar proportion in adult and young individuals. This is the first report of zoonotic and non-zoonotic helminths parasites in rodents from Guatemala.

Distinct evolutionary trajectories of V1R clades across mouse species

BACKGROUND

Many animals rely heavily on olfaction to navigate their environment. Among rodents, olfaction is crucial for a wide range of social behaviors. The vomeronasal olfactory system in particular plays an important role in mediating social communication, including the detection of pheromones and recognition signals. In this study we examine patterns of vomeronasal type-1 receptor (V1R) evolution in the house mouse and related species within the genus Mus. We report the extent of gene repertoire turnover and conservation among species and clades, as well as the prevalence of positive selection on gene sequences across the V1R tree. By exploring the evolution of these receptors, we provide insight into the functional roles of receptor subtypes as well as the dynamics of gene family evolution.

RESULTS

We generated transcriptomes from the vomeronasal organs of 5 Mus species, and produced high quality V1R repertoires for each species. We find that V1R clades in the house mouse and relatives exhibit distinct evolutionary trajectories. We identify putative species-specific gene expansions, including a large clade D expansion in the house mouse. While gene gains are abundant, we detect very few gene losses. We describe a novel V1R clade and highlight candidate receptors for future study. We find evidence for distinct evolutionary processes across different clades, from largescale turnover to highly conserved repertoires. Patterns of positive selection are similarly variable, as some clades exhibit abundant positive selection while others display high gene sequence conservation. Based on clade-level evolutionary patterns, we identify receptor families that are strong candidates for detecting social signals and predator cues. Our results reveal clades with receptors detecting female reproductive status are among the most conserved across species, suggesting an important role in V1R chemosensation.

CONCLUSION

Analysis of clade-level evolution is critical for understanding species' chemosensory adaptations. This study provides clear evidence that V1R clades are characterized by distinct evolutionary trajectories. As receptor evolution is shaped by ligand identity, these results provide a framework for examining the functional roles of receptors.

Improved mitochondrial coupling as a response to high mass-specific metabolic rate in extremely small mammals

Mass-specific metabolic rate negatively co-varies with body mass from the whole-animal to the mitochondrial levels. Mitochondria are the mainly consumers of oxygen inspired by mammals to generate ATP or compensate energetic losses dissipated as the form of heat (proton leak) during oxidative phosphorylation. Consequently, ATP synthesis and proton leak thus compete for the same electrochemical gradient. Because proton leak co-varies negatively with body mass, it is unknown if extremely small mammals further decouple their mitochondria to maintain their body temperature or if they implement metabolic innovations to ensure cellular homeostasis. The present study investigates the impact of body mass variation on cellular and mitochondrial functioning in small mammals, comparing the two extremely small African pygmy mice (Mus mattheyi, approx. 5 g and Mus minutoides, approx. 7 g) with the larger house mouse (Mus musculus, approx. 22 g). Oxygen consumption rates were measured from the animal to the mitochondrial levels. We also measured mitochondrial ATP synthesis in order to appreciate the mitochondrial efficiency (ATP/O). At the whole-animal scale, mass- and surface-specific metabolic rates co-varied negatively with body mass, whereas this was not necessarily the case at cellular and mitochondrial levels. M. mattheyi had generally the lowest cellular and mitochondrial fluxes, depending on the tissue considered (liver or skeletal muscle), as well as having higher efficient muscle mitochondria than the other two species. M. mattheyi presents metabolic innovations to ensure its homeostasis, by generating more ATP per oxygen consumed.

Demography of the salt marsh harvest mouse (Reithrodontomys raviventris halicoetes) and associated rodents in tidal and managed wetlands

Suisun Marsh (Solano County, California) is the largest contiguous marsh remaining on the West Coast of the United States, and makes up approximately 10% of the wetlands remaining in the San Francisco Estuary. Suisun Marsh has been safeguarded from development through the operation of over 100 privately owned waterfowl hunting clubs, which manage for diked waterfowl habitat. However, this management—and the subsequent loss of tidal influence—has been considered harmful for some species, including the endangered salt marsh harvest mouse (SMHM; Reithrodontomys raviventris). To determine the value of tidal wetlands relative to those managed for waterfowl, we performed periodic surveys for rodents in managed and tidal wetlands over 5 years, and used capture-mark-recapture analyses to estimate demographic parameters and abundance for the three most common rodents—the northern SMHM (R. r. halicoetes), the western harvest mouse (a sympatric native species; R. megalotis, WHM), and the house mouse (a sympatric invasive species; Mus musculus). Wetland type had no effect on detection, temporary emigration, or survival for any of these species. However, fecundity and population growth for all three species were affected by an interaction of season and wetland type, although none of these parameters was consistently superior in either habitat type. Estimated abundance of SMHM and Mus was similar in both wetland types, whereas WHM were more abundant in managed wetlands. Salt marsh harvest mice also showed no affinity for any microhabitat characteristics associated with tidal wetlands. Managed wetlands in Suisun Marsh support SMHM and Mus equally, and abundances of WHM were greater than in tidal wetlands, suggesting managed wetlands may be superior in terms of supporting native rodents. As climate change and sea level rise are predicted to threaten coastal marshes, these results suggest the recovery strategy for SMHM could incorporate managed wetlands.

Fossils know it best: Using a new set of fossil calibrations to improve the temporal phylogenetic framework of murid rodents (Rodentia: Muridae)

Murid rodents (Rodentia: Muridae) represent the most diverse and abundant mammalian family. In this study, we provide a refined set of fossil calibrations which is used to reconstruct a dated phylogeny of the family using a multilocus dataset (six nuclear and nine mitochondrial gene fragments) encompassing 161 species representing 82 murid genera from four extant subfamilies (Deomyinae, Gerbillinae, Lophiomyinae and Murinae). In comparison with previous studies on murid or muroid rodents, our work stands out for the implementation of nine robust fossil constraints within the Muridae thanks to a thorough review of the fossil record. Before being assigned to specific nodes of the phylogeny, all potential fossil constraints were carefully assessed; they were also subjected to several cross-validation analyses. The resulting phylogeny is consistent with previous phylogenetic studies on murids, and recovers the monophyly of all sampled murid subfamilies and tribes. Based on nine controlled fossil calibrations, our inferred temporal timeframe indicates that the murid family likely originated in the course of the Early Miocene, 22.0-17.0 million years ago (Ma), and that most major lineages (i.e. tribes) started diversifying ca. 10 Ma. Historical biogeography analyses support the tropical origin for the family, with an initial internal split (vicariance event) between Afrotropical and Oriental (Indomalaya and Philippines) lineages. During the course of their diversification, the biogeographic pattern of murids is marked by several dispersal events toward the Australasian and the Palearctic regions. The Afrotropical region was also secondarily colonized at least three times from the Indomalaya, indicating that the latter region has acted as a major centre of diversification for the family.

Environmental Factors and Zoonotic Pathogen Ecology in Urban Exploiter Species

Knowledge of pathogen ecology, including the impacts of environmental factors on pathogen and host dynamics, is essential for determining the risk that zoonotic pathogens pose to people. This review synthesizes the scientific literature on environmental factors that influence the ecology and epidemiology of zoonotic microparasites (bacteria, viruses and protozoa) in globally invasive urban exploiter wildlife species (i.e., rock doves [Columba livia domestica], European starlings [Sturnus vulgaris], house sparrows [Passer domesticus], Norway rats [Rattus norvegicus], black rats [R. rattus] and house mice [Mus musculus]). Pathogen ecology, including prevalence and pathogen characteristics, is influenced by geographical location, habitat, season and weather. The prevalence of zoonotic pathogens in mice and rats varies markedly over short geographical distances, but tends to be highest in ports, disadvantaged (e.g., low income) and residential areas. Future research should use epidemiological approaches, including random sampling and robust statistical analyses, to evaluate a range of biotic and abiotic environmental factors at spatial scales suitable for host home range sizes. Moving beyond descriptive studies to uncover the causal factors contributing to uneven pathogen distribution among wildlife hosts in urban environments may lead to targeted surveillance and intervention strategies. Application of this knowledge to urban maintenance and planning may reduce the potential impacts of urban wildlife-associated zoonotic diseases on people.

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.

EVOLUTION OF β-GLUCURONIDASE REGULATION IN THE GENUS MUS

Despite the central role suggested for regulatory mutations in many evolutionary scenarios, there is relatively little information available about the type and extent of regulatory differences between species, or to what extent differences between species are independent of variation within species. To address this issue we have studied the regulatory system of β-glucuronidase, a gene implicated in a murine androgen-inducible pheromone-signalling system. We examined the changes in β-glucuronidase hormonal regulation which have occurred during the radiation of a group of 12 closely related species of mice by assaying β-glucuronidase activity in six different tissues after treatment with estrogen, and with androgen alone and in combination with either estrogen or growth hormone. We also examined in some detail the extent of variation in regulatory responses within species. We found extensive variation in regulatory phenotypes both within and among the species surveyed, suggesting that many of the species examined are currently polymorphic for various regulatory factors that affect inducibility of β-glucuronidase. The variation we observed reflects changes in the ability of the β-glucuronidase gene to respond to hormonal influences, rather than changes in aspects of the hormonal signalling system exterior to the gene. The marked differences among species in the renal and uterine responses to hormonal induction of β-glucuronidase are not easily related to the phylogeny of the genus Mus. If hormonal induction of the gene for β-glucuronidase is subject to natural selection, it appears to be subject to widely fluctuating selective forces. We review evidence that the apparently disorderly evolution of the hormonal responsiveness of β-glucuronidase does not appear to be a unique property of this regulatory system. In contrast to the evolution of many protein sequences, which are tightly correlated with phylogeny and proceed at a relatively constant rate, some, perhaps many, regulatory phenotypes are in rapid evolutionary flux, providing an extensive range of phenotypes upon which selection can act.

The laboratory mouse and wild immunology

The laboratory mouse, Mus musculus domesticus, has been the workhorse of the very successful laboratory study of mammalian immunology. These studies--discovering how the mammalian immune system can work--have allowed the development of the field of wild immunology that is seeking to understand how the immune responses of wild animals contributes to animals' fitness. Remarkably, there have hardly been any studies of the immunology of wild M. musculus domesticus (or of rats, another common laboratory model), but the general finding is that these wild animals are more immunologically responsive, compared with their laboratory domesticated comparators. This difference probably reflects the comparatively greater previous exposure to antigens of these wild-caught animals. There are now excellent prospects for laboratory mouse immunology to make major advances in the field of wild immunology.

Lassa serology in natural populations of rodents and horizontal transmission

Lassa virus causes hemorrhagic fever in West Africa. Previously, we demonstrated by PCR screening that only the multimammate mouse, Mastomys natalensis, hosts Lassa virus in Guinea. In the present study, we used the same specimen collection from 17 villages in Coastal, Upper, and Forest Guinea to investigate the Lassa virus serology in the rodent population. The aim was to determine the dynamics of antibody development in M. natalensis and to detect potential spillover infections in other rodent species. Immunoglobulin G (IgG) antibody screening was performed using the indirect immunofluorescence assay with the Guinean Lassa virus strain Bantou 289 as antigen. The overall seroprevalence was 8% (129/1551) with the following rodents testing positive: 109 M. natalensis, seven Mastomys erythroleucus, four Lemniscomys striatus, four Praomys daltoni, three Mus minutoides, and two Praomys rostratus. Nearly all of them (122/129) originated from Bantou, Tanganya, and Gbetaya, where Lassa virus is highly endemic in M. natalensis. The antibody seroprevalence in M. natalensis from this high-endemic area (27%; 108/396) depended on the village, habitat, host age, and host abundance. A main positive factor was age; the maximum seroprevalence reached 50% in older animals. Our data fit with a model implicating that most M. natalensis rodents become horizontally infected, clear the virus within a period significantly shorter than their life span, and develop antibodies. In addition, the detection of antibodies in other species trapped in the habitats of M. natalensis suggests spillover infections.

Nasal aerodynamics protects brain and lung from inhaled dust in subterranean diggers, Ellobius talpinus

Inhalation of air-dispersed sub-micrometre and nano-sized particles presents a risk factor for animal and human health. Here, we show that nasal aerodynamics plays a pivotal role in the protection of the subterranean mole vole Ellobius talpinus from an increased exposure to nano-aerosols. Quantitative simulation of particle flow has shown that their deposition on the total surface of the nasal cavity is higher in the mole vole than in a terrestrial rodent Mus musculus (mouse), but lower on the olfactory epithelium. In agreement with simulation results, we found a reduced accumulation of manganese in olfactory bulbs of mole voles in comparison with mice after the inhalation of nano-sized MnCl2 aerosols. We ruled out the possibility that this reduction is owing to a lower transportation from epithelium to brain in the mole vole as intranasal instillations of MnCl2 solution and hydrated nanoparticles of manganese oxide MnO · (H2O)x revealed similar uptake rates for both species. Together, we conclude that nasal geometry contributes to the protection of brain and lung from accumulation of air-dispersed particles in mole voles.

Multiple paternity does not depend on male genetic diversity

Polyandry is common in many species and it has been suggested that females engage in multiple mating to increase the genetic diversity of their offspring (genetic diversity hypothesis). Multiple paternity occurs in 30% of litters in wild populations of house mice, Mus musculus musculus, and multiple-sired litters are genetically more diverse than single-sired ones. Here, we aimed to test whether female house mice produce multiple-sired litters when they have the opportunity to produce genetically diverse litters. We assessed the rates of multiple paternity when females could choose to mate with two males that were genetically dissimilar to each other (i.e. nonsiblings and MHC dissimilar) compared with when females could choose to mate with two males that were genetically similar to each other (i.e. siblings and shared MHC alleles). Multiple mating may depend upon a female's own condition, and, therefore, we also tested whether inbred (from full-sibling matings) females were more likely to produce multiple-sired progeny than outbred controls. Overall we found that 29% of litters had multiple sires, but we found no evidence that females were more likely to produce multiple-sired litters when they had the opportunity to mate with genetically dissimilar males compared with controls, regardless of whether females were inbred or outbred. Thus, our findings do not support the idea that female mice increase multiple paternity when they have the opportunity to increase the genetic diversity of their offspring, as expected from the genetic diversity hypothesis.

The complex ontogenetic trajectory of mandibular shape in a laboratory mouse

The mouse mandible is a popular model system that continues to be the focus of studies in evo-devo and other fields. Yet, little attention has been given to the role of postnatal growth in producing the adult form. Using cleared and stained specimens, we describe the timing of tooth and jaw development and changes in jaw size and shape from postnatal day 1 (p1) through weaning to adulthood. We found that tooth development is relatively advanced at birth, and that the functional adult dentition is in place by p15 (just before the start of weaning). Shape analysis showed that the trajectory of mandible shape changes direction at least twice between birth and adulthood, at p7 and p15. At each stage there are changes in shape to all tooth- and muscle-bearing regions and, at each change of direction, all of these regions change their pattern of growth. The timing of the changes in direction in Mus suggests there are signals that redirect growth patterns independently of changes in function and loading associated with weaning and jaw muscle growth. A better understanding of these signals and how they produce a functionally integrated mandible may help explain the mechanisms guiding evolutionary trends and patterns of plasticity and may also provide valuable clues to therapeutic manipulation of growth to alleviate the consequences of trauma or disease.