What Is Peromyscus? Evidence from nuclear and mitochondrial DNA sequences suggests the need for a new classification

Three distinct forms of Pneumocystis coexist in individuals of two species of deer mice (genus Peromyscus)

As emerging zoonoses represent a significant public health threat, understanding how pathogens' host ranges evolve is critical to protect human and wildlife health. Closely related hosts infected with host-specific pathogens provide valuable opportunities for clear inferences of host range evolution, as they allow for the examination of early diversification patterns in their resident pathogens. Pneumocystis, an obligate lung symbiont that is believed to be ubiquitous in mammals, exemplifies such a model. To explore the early stages of divergence in Pneumocystis, we collected geographically dispersed samples from two sister species of deer mice: Peromyscus leucopus (white-footed mice) and Peromyscus gossypinus (cotton mice). We sequenced two nuclear and two mitochondrial loci of Pneumocystis sampled from the lungs of these mice. These sequences revealed three distinct Pneumocystis taxa, two of which were found to cross-infect both host species and were often found coexisting within the same individual. Genetic diversity and phylogenetic analysis suggest that the three Pneumocystis taxa represent separate species. Further analysis of the mitochondrial large subunit rRNA gene from the most common taxon of these three revealed that host geographic origins influenced Pneumocystis genetic structure more than host species identity. Nevertheless, the results also suggest an overall interconnectedness of the symbiont metapopulation.

Unveiling hidden diversity: Phylogenomics of neotomine rodents and taxonomic implications for the genus Peromyscus

Neotomine rodents (Cricetidae, Neotominae) represent one of the most commonly encountered and diverse group of rodents in North America, yet phylogenetic relationships within this group remain uncertain. This subfamily is known for its rapid evolution, adding more complexity to our efforts to unravel their evolutionary history. The main debate revolves around the recognition of the genus Peromyscus as monophyletic or paraphyletic due to its relationship with other genera such as Habromys, Megadontomys, Podomys, Neotomodon, and Osgoodomys. Here, we aim to resolve phylogenetic relationships within Neotominae, to further explore their evolutionary history and taxonomic boundaries. We used target capture and high-throughput sequencing of complete mitogenomes and thousands of genome-wide ultraconserved elements loci (UCEs). Our comprehensive analyses encompassed 53 species of Neotominae spanning 12 previously described genera, along with one yet-undescribed genus. We also investigated 12 out of the 13 species groups within Peromyscus. Our analyses, including Maximum Likelihood and Bayesian Inference with both mitogemomes and UCEs, as well as the coalescent species-tree-based approach with UCEs, consistently recovered concordant and well-resolved phylogenies with high levels of nodal support. We identified seven main clades within Neotominae that could potentially be recognized at the generic level, mostly to categorize the genus Peromyscus as a monophyletic group, including one species group within "Peromyscus". Furthermore, our divergence dating estimates place the crown age of Neotominae to be around the late Miocene at ca. 7.9 - 10.7 mya. While generic level diversification continued through the Pliocene, species level diversification predominantly occurred during the late Pliocene, extending through the Pleistocene and Holocene. These epochs have been recognized as periods with significant changes in flora and fauna, driving ecological transformations on a global scale. We hypothesized that climatic and vegetation shifts during the Neogene and Quaternary, coupled with geological events, topographical features, and the presence of biogeographical corridors played a pivotal role in the speciation and diversification of Neotominae. Recognizing the importance of generating genomic-scale data coupled with a broad taxonomic sampling, our study, for the first time, offers resolution of the relationships among the main lineages of Neotominae. We expect that the phylogeny presented here will serve as a foundational resource for future systematic and evolutionary studies. This includes facilitating a proper comprehensive taxonomic revision of the group and the formal description and naming of new genera.

A Systematic Review of the Distribution and Prevalence of Viruses Detected in the Peromyscus maniculatus Species Complex (Rodentia: Cricetidae)

The North American Deermouse, Peromyscus maniculatus, is one of the most widespread and abundant mammals on the continent. It is of public health interest as a known host of several viruses that are transmissible to humans and can cause illness, including the acute respiratory disease Hantavirus pulmonary syndrome (HPS). However, recent taxonomic studies indicate that P. maniculatus is a complex of multiple species, raising questions about how to identify and interpret three decades of hantavirus monitoring data. We conducted a systematic review investigating the prevalence and spatial distribution of viral taxa detected in wild populations allocated to P. maniculatus. From the 49 relevant studies published from 2000 to 2022, we extracted and analyzed spatial occurrence data to calculate weighted populational prevalences for hantaviruses. We found that detection efforts have been concentrated in the Western United States and Mexico with a focus on the spread of Sin Nombre virus (Orthohantavirus sinnombreense), the primary causative agent of HPS. There are significant gaps in the existing literature both geographically and in regard to the types of viruses being sampled. These results are significantly impacted by a recent taxonomic split of P. maniculatus into four species, and we were able to update 94% of hantavirus observations to reflect this change. Investigating the uncertain, and likely multiple, phylogenetic histories of these viral hosts should be a key emphasis of future modeling efforts.

Hearing in Two Closely Related Peromyscus Species (Peromyscus maniculatus and P. leucopus)

The genus Peromyscus has been extensively used as a model for ecological, behavioral, and evolutionary investigations. We used auditory brainstem responses (ABRs), craniofacial morphology, and pinna measurements to compare characteristics that impact hearing in two wild-caught species, P. leucopus P. maniculatus. We observed significant statistical differences in craniofacial and pinna attributes between species with P. leucopus overall exhibiting larger features than P. maniculatus. ABR recordings indicated that both species showed similar best frequency thresholds between 8-24 kHz. We found significant effects of intensity on amplitude ratio of wave I and IV for P. maniculatus, but not P. leucopus and effects of wave number on slope of the latency-intensity function with higher wave IV and shorter wave I slope of latency intensity function in P. leucopus. Finally, the data showed significant differences in latency shift of the DN1 component of the BIC in relation to ITD between species, while no significant differences were observed across relative DN1 amplitude. This study supports the used of P. leucopus and P. maniculatus as future model species for auditory research.

Mating system variation and gene expression in the male reproductive tract of Peromyscus mice

Genes involved in reproduction often evolve rapidly at the sequence level due to postcopulatory sexual selection (PCSS) driven by male-male competition and male-female sexual conflict, but the impact of PCSS on gene expression has been under-explored. Further, though multiple tissues contribute to male reproductive success, most studies have focused on the testes. To explore the influence of mating system variation on reproductive tract gene expression in natural populations, we captured adult males from monogamous Peromyscus californicus and polygynandrous P. boylii and P. maniculatus. We generated RNAseq libraries, quantified gene expression in the testis, seminal vesicle, epididymis, and liver, and identified 3627 mating system-associated differentially expressed genes (MS-DEGs), where expression shifted in the same direction in P. maniculatus and P. boylii relative to P. californicus. Gene expression variation was most strongly associated with mating behaviour in the seminal vesicles, where 89% of differentially expressed genes were MS-DEGs, including the key seminal fluid proteins Svs2 and Pate4. We also used published rodent genomes to test for positive and relaxed selection on Peromyscus-expressed genes. Though we did not observe more overlap than expected by chance between MS-DEGs and positively selected genes, 203 MS-DEGs showed evidence of positive selection. Fourteen reproductive genes were under tree-wide positive selection but convergent relaxed selection in P. californicus and Microtus ochrogaster, a distantly related monogamous species. Changes in transcript abundance and gene sequence evolution in association with mating behaviour suggest that male mice may respond to sexual selection intensity by altering aspects of sperm motility, sperm-egg binding and copulatory plug formation.

Novel genomic resources contribute to the systematics of threatened arboreal deer mice of the genus Habromys Hooper & Musser, 1964 (Cricetidae, Neotominae) within a neotomine-peromyscine phylogeny

The Crested-tailed deer mouse, Habromyslophurus, is one of seven arboreal species within the genus Habromys. Species of this genus are monotypic, relatively rare, and occur in low densities. Their geographical distribution is highly fragmented due to being restricted to montane cloud forest in Mesoamerica and they are of conservation concern. All Habromys species are endemic to Mexico, except H.lophurus, which is also distributed in Guatemala and El Salvador. In this study, we obtained and characterized the first mitogenome and several thousand nuclear ultraconserved elements (UCEs) of H.lophurus to determine its phylogenetic position within neotomine-peromyscine mice. Its mitogenome sequence (16,509 bp) is only the second complete mitogenome obtained for this poorly known genus. We also obtained the first nuclear genomic data for H.lophurus, including 3,654 UCE loci, as well as a partial mitogenome of H.simulatus (6,349 bp), and 2,186 UCE for the outgroup Holochilussciureus. Phylogenetic analyses that included our newly generated genomic data coupled with previously published data from other neotomine-peromyscine mice confirm the placement of H.lophurus, H.simulatus, and H.ixtlani within a highly supported clade. The Habromys clade was nested within a clade that also contains members of the genus Peromyscus and provides further support for the hypothesis of the paraphyly of Peromyscus. These genomic resources will contribute to future phylogenomic studies that aim to further elucidate the evolutionary history of this rare and critically endangered genus of rodents.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900 bp full-length element of mys-related sequences (mysRS), an 8000 bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800 bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

Evolution in reproductive tempo and investment across the Peromyscus radiation

Mammals display diverse reproductive strategies, however, the ultimate and proximate mechanisms that underlie this diversity and its composite traits remain poorly understood from both evolutionary and physiological perspectives. The Peromyscus genus of rodents, which is found throughout the north and central Americas, has diversified along life history gradients, varying both within and among species in reproductive strategies. This variation provides a useful model for studying reproductive diversity. Here, we combine a literature review with new analyses of captive colony breeding records from six Peromyscus species to assess our current understanding of how plasticity and local adaptation contribute to diversity in two classes of reproductive traits: phenology and litter investment. There is substantial evidence that many traits underlying phenology and litter investment have diverged among populations in ways that are likely to be locally adaptive, though plasticity in these traits remains common. However, these conclusions are largely based on data collected from the two most widespread Peromyscus species: P. maniculatus and P. leucopus. The majority of Peromyscus species diversity remains understudied regarding reproductive phenology and litter traits. We conclude by discussing key challenges and considerations relevant to using Peromyscus as a mammalian model for reproductive trait diversity and evolution moving forward.

Morphology and genetics of grasshopper mice revisited in a paleontological framework: reinstatement of Onychomyini (Rodentia, Cricetidae)

Grasshopper mice of the genus Onychomys, represented by three living species in North America, have a long and controversial taxonomic history. Usually allocated to either the cricetine or neotomine cricetids, they also have been considered to represent a distinct tribe. Since the discovery and description of the extinct grasshopper mouse relative Acrolophomys rhodopetros from the late Miocene of the upper Dove Spring Formation of California, dated at 9.3–8.8 Ma, it has become apparent that the grasshopper mouse clade has a long, distinct evolutionary history. Using a combination of morphological (including paleontological material) and molecular data, we reassessed the phylogenetic position of grasshopper mice. A morphological phylogenetic analysis was done on fossil and modern specimens of all recognized neotomine tribes, including craniodental, phallic, and soft tissue characters. A DNA-based matrix was constructed including 72 species representing all known living genera of Neotominae and 13 outgroup taxa belonging mostly to cricetid subfamilies. DNA sampling covered the mitochondrial protein-coding gene cytochrome-b (Cytb), and seven nuclear loci. The morphological analysis yielded a single most parsimonious tree of 42 steps, placing Ochrotomys (Ochrotomyini), Baiomys (Baiomyini), Reithrodontomys (Reithrodontomyini), and an Onychomys–Acrolophomys clade as successive sister clades to a Peromyscus clade, respectively. The molecular phylogenetic analyses recovered seven major clades: (1) a clade including Habromys, Megadontomys, Neotomodon, Osgoodomys, Podomys, and a paraphyletic Peromyscus clade, sister to (2) a second clade containing extant Onychomys species, (3) a Reithrodontomys clade, (4) an Isthmomys clade, (5) a clade including Baiomys and Scotinomys, (6) an Ochrotomys clade, and (7) a well-­supported clade containing Hodomys, Neotoma, and Xenomys. A Bayesian combined morphological and molecular analysis recovered the same major phylogenetic associations as the molecular analyses. The sum of molecular markers and morphological traits expressed by Acrolophomys and Onychomys leads to a phylogenetic position supporting their recognition as a distinct tribe.

Museomics and the holotype of a critically endangered cricetid rodent provide key evidence of an undescribed genus

Historical DNA obtained from voucher specimens housed in natural history museums worldwide have allowed the study of elusive, rare or even extinct species that in many cases are solely represented by museum holdings. This has resulted in the increase of taxonomic representation of many taxa, has led to the discovery of new species, and has yielded stunning novel insights into the evolutionary history of cryptic or even undescribed species. Peromyscus mekisturus, is a critically endangered cricetid rodent endemic to Mexico and is only known from two museum specimens collected in 1898 and 1947. Intensive field work efforts to attempt to determine if viable populations still exist have failed, suggesting that this species is extinct or is nearing extinction. In addition, a recent study using mitogenomes demonstrated that P. mekisturus forms a well-supported clade outside the genus Peromyscus and hypothesized that this taxon is the sister group of the genus Reithrodontomys. Here, we used target enrichment and high-throughput sequencing of several thousand nuclear ultraconserved elements and mitogenomes to reconstruct dated phylogenies to test the previous phylogenetic hypothesis. We analyzed the holotype and the only other known specimen of P. mekisturus and museum samples from other peromyscine rodents to test the phylogenetic position of the species. Our results confirm that the only two specimens known to science of P. mekisturus belong to the same species and support the hypothesis that this species belongs to an undescribed genus of cricetid rodents that is sister to the genus Reithrodontomys. We dated the origin of P. mekisturus together with other speciation events in peromyscines during the late Pliocene – early Pleistocene and related these events with the Pleistocene climatic cycles. In light of our results, we recommend a taxonomic re-evaluation of this enigmatic species to properly recognize its taxonomic status as a new genus. We also acknowledge the relevance of generating genomic data from type specimens and highlight the need and importance of continuing to build the scientific heritage of the collections to study and better understand past, present, and future biodiversity.

Two new species of Peromyscus (Cricetidae: Neotominae) from the Transverse Volcanic Belt of Mexico

Specimens of the Peromyscus boylii species group distributed in the western and northeastern montane regions of Michoacán, México, historically have been assigned to P. levipes. Previous studies indicated that these specimens possessed mitochondrial DNA haplotypes that were distinct from both P. levipes and P. kilpatricki, a recently named species in the P. boylii species group from northeastern Michoacán and western Morelos. Herein karyotypic, DNA sequence, and morphological data were analyzed from those populations to evaluate their taxonomic affinity. Karyotypic data indicated that individuals from western Michoacán (Dos Aguas and Aguililla) and from a newly discovered population in northeastern Michoacán (Zinapécuaro) were chromosomally similar to P. carletoni (FN = 68) but distinct from other taxa assigned to the P. boylii species group. Analyses of cranial characteristics indicated that, relative to other species in the P. boylii species group, two morphologically distinct groups were present that corresponded to the Dos Aguas/Aguililla and Zinapécuaro populations, respectively. The latter population, although represented by a small sample size (n = 5 specimens), appeared to exhibit some trenchant morphological distinctions compared with other cryptic species in the P. boylii group. Phylogenetic analyses (parsimony, Bayesian, and likelihood) of DNA sequences obtained from the mitochondrial cytochrome-b gene indicated that although the individuals from Dos Aguas/Aguililla and Zinapécuaro formed a sister group relationship, they formed monophyletic clades that differed genetically (2.54%)-a level approaching that seen between other sister species of Peromyscus. Further, the Dos Aguas/Aguililla and Zinapécuaro clade was more closely aligned with a clade containing representatives of P. carletoni and P. levipes instead of with those from closer geographic proximities (P. kilpatricki) located in eastern Michoacán. Together, these results indicated that these two populations seemingly represent two undescribed species in the P. boylii species group for which we propose the names Peromyscus greenbaumi for populations in western Michoacán (circa Dos Aguas and Aguililla) and Peromyscus ensinki for populations in northeastern Michoacán (circa Zinapécuaro).

Mitochondrial DNA and other lines of evidence clarify species diversity in the Peromyscus truei species group (Cricetidae: Neotominae)

Deer mice (genus Peromyscus) are among the commonest small mammals in the Nearctic zoogeographic region. Nevertheless, systematic relationships are only partially settled and numerous taxonomic questions await resolution. For instance, researchers have found that some members of the Peromyscus truei species group contain high levels of genetic divergence that could indicate the presence of cryptic species. We analyzed the systematics and phylogenetic relationships of the P. truei group using new and previously published mitochondrial cytochrome b sequences. Our analyses verify several earlier conclusions, but we also detected new clades that deserve recognition. Considering their mitochondrial distinctiveness, allopatric ranges, and previously reported molecular, biochemical, chromosomal, morphological, and ecological differences, we elevate three previously described taxa to species. We support the recognition of two subgroupings. The first comprises P. gratus, P. truei, and possibly P. cf. martirensis and P. cf. zapotecae. The second contains to P. amplus, P. attwateri, P. collinus, P. difficilis, P. felipensis, P. laceianus, P. nasutus, P. ochraventer, and P. pectoralis. Placement of P. bullatus will likely remain unknown until genetic data are available. Further research could improve our understanding of the evolutionary history of Peromyscus, but in some cases taxonomic issues must be resolved first.

Molecular systematics of the Reithrodontomys tenuirostris group (Rodentia: Cricetidae) highlighting the Reithrodontomys microdon species complex

The Reithrodontomys tenuirostris species group is considered "the most specialized" within the genus Reithrodontomys from morphological and ecological perspectives. Previous studies based on molecular data recommended changes in the taxonomy of the group. In particular, R. microdon has been the most taxonomically questioned, with the suggestion that it constitutes a complex of cryptic species. We analyzed the phylogenetic relationships of the R. tenuirostris species group using DNA sequences from the mitochondrial Cytochrome b gene and Intron 7 of the nuclear beta fibrinogen gene. In addition, divergence times were estimated, and possible new taxa delimited with three widely used species delimitation methods. Finally, possible connectivity routes based on shared haplotypes were tested among the R. microdon populations. All species were recovered as monophyletic with the exception of R. microdon, whose individuals were grouped into four different haplogroups, one of which included specimens of R. bakeri. Diversification within the R. tenuirostris species group began about 3 Ma, in the Pleistocene. The bGMYC and STACEY delimitation methods were congruent with each other, delimiting at the species-level each haplogroup within R. microdon, while the mPTP suggested a greater number of species. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a higher number of species in the R. tenuirostris group, because we show that there are four species within what is currently recognized as R. microdon.

Mitochondrial capture and subsequent genetic divergence generates a novel haplogroup: evidence from ancient and ongoing hybridization in mule and white-tailed deer

Odocoileus virginianus (white-tailed deer) and O. hemionus (mule deer) are sympatric across much of North America. Molecular evidence suggests that up to 24% of individuals in some populations are a product of hybrid ancestry. Several studies have alluded to ancient and recent introgression between Odocoileus spp.; however, no divergence dates were proposed. Herein, phylogenetic analyses of DNA sequences obtained from the mitochondrial Cytochrome b gene in 690 individuals identified three clades corresponding to black-tailed deer, white-tailed deer, or a unique combination of both white-tailed deer and mule deer. White-tailed deer and mule deer diverged from a common ancestor of approximately 3.13 mya followed by an ancient hybridization event of approximately 1.32 mya, in which the white-tailed deer mitochondrial genome was “captured” by mule deer. This hybridization event produced a novel haplogroup for white-tailed deer and mule deer located west of the Appalachian Mountains and east of the Cascade Range, south to Veracruz, Mexico, and north to the Yukon Territory, Canada. The ancestral mule deer-like mitochondrial genome appears to be restricted to black-tailed deer distributed along the western portion of the Cascade and Sierra Nevada Ranges of the United States and Canada, whereas the ancestral white-tailed deer-like mitochondrial genome is restricted to the eastern United States and portions of Latin America and Caribbean regions. The “captured mitochondrial genome” has continued on an independent evolutionary trajectory and represents a unique and broadly distributed haplogroup that is 7.25% and 2.84% different from the ancestral mule deer and ancestral white-tailed deer haplogroups, respectively.

Variation and Selection in the Putative Sperm-Binding Region of ZP3 in Muroid Rodents: A Comparison between Cricetids and Murines

In mammals, the zona pellucida glycoprotein 3 (ZP3) is considered a primary sperm receptor of the oocyte and is hypothesized to be involved in reproductive isolation. We investigated patterns of diversity and selection in the putative sperm-binding region (pSBR) of mouse ZP3 across Cricetidae and Murinae, two hyperdiverse taxonomic groups within muroid rodents. In murines, the pSBR is fairly conserved, in particular the serine-rich stretch containing the glycosylation sites proposed as essential for sperm binding. In contrast, cricetid amino acid sequences of the pSBR were much more variable and the serine-rich motif, typical of murines, was generally substantially modified. Overall, our results suggest a general lack of species specificity of the pSBR across the two muroid families. We document statistical evidence of positive selection acting on exons 6 and 7 of ZP3 and identified several amino acid sites that are likely targets of selection, with most positively selected sites falling within or adjacent to the pSBR.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

Taxonomic and phylogenetic beta diversity of cricetid rodents in Oaxaca, southern Mexico

Understanding the ecological and historical causes and processes that shape biodiversity distribution patterns remains a challenging and fundamental task in biogeography, ecology, and evolution. To address this issue, taxonomic and phylogenetic β diversity can help us to assess the importance of ecological and historical factors that structure these biotic patterns. To make inferences about the processes underlying current spatial patterns in communities of Cricetidae across the state of Oaxaca, Mexico, their taxonomic and phylogenetic β diversity were assessed jointly. Our aims were: 1) to examine taxonomic and phylogenetic β diversity and their turnover and nestedness components among physiographic subprovinces; 2) to test for statistical significance of observed phylogenetic β diversity against the expected values of a null model; and 3) to evaluate if these metrics were correlated with geographical distance. We obtained the species composition for 12 subprovinces based on distribution models for 49 cricetid species present in Oaxaca, then carried out a maximum likelihood analysis to estimate their phylogenetic relationships. Our results show that the taxonomic and phylogenetic dissimilarities mainly were explained by the turnover component of species and lineages. In almost all pairwise comparisons, the null model approach revealed random patterns for phylogenetic β diversity values and its components. Mantel correlation models showed that the values of total taxonomic and phylogenetic diversity and their components are correlated with the geographical distances between subprovinces. Our results suggest that both taxonomic and phylogenetic β diversity are explained by the interplay between biogeographical history from southern Mexico, and the recent speciation processes in cricetid rodents. Given that speciation processes are allopatric for most cricetid taxa, the high values of spatial turnover can be explained by the small ranges of species, coupled with current abiotic conditions that act as filters, promoting specialization of species on particular conditions. Our results show the importance of the phylogenetic approach to unravel the multidimensional spatial patterns of biodiversity.

Genetic and morphological variation of Oxymycterus (Rodentia: Sigmodontinae) in the Brazilian Atlantic Forest

We present a new assessment of the genetic and morphological variation within Oxymycterus quaestor Thomas, 1903, which currently includes the junior synonyms O. judex Thomas, 1909 and O. misionalis Sanborn, 1931. We integrate distinct lines of evidence, including variation of mitochondrial (Cytochrome b [Cytb]) and nuclear (intron 7 of beta fibrinogen gene [Fgb]) sequences, and the assessment of skull quantitative traits based on geometric morphometrics, throughout the Atlantic Forest of Southeastern-Southern Brazil, Argentina, and Paraguay. Phylogenetic relationships based on Cytb indicate that O. quaestor is structured in four well-supported clades (lineages A–D), one of them (lineage C) including topotypes of a previously associated nominal form (O. judex). However, these Cytb lineages exhibit lower levels of differentiation based on the Fgb locus, and are not recovered in the genealogies of this nuclear marker, representing a case of mitonuclear discordance. The Cytb lineages also broadly overlapped in the morphospace both in skull shape and size, which sustain the current wider concept of O. quaestor as one single young species (0.947 Myr) that is recently expanding, and ultimately branching out, in the Atlantic Forest.

Gross stomach morphology in akodontine rodents (Cricetidae: Sigmodontinae: Akodontini): a reappraisal of its significance in a phylogenetic context

Akodontini, the second largest tribe within sigmodontine rodents, encompasses several stomach morphologies. This is striking because most sigmodontine groups of comparable taxonomic rank are very conservative in this respect. Based on extensive sampling of newly dissected specimens (213 stomachs representing 36 species), as well as published examples, covering almost all akodontine living genera (15 of 16), we undertook a reappraisal of the gross morphology of this organ. We then mapped this information, together with gallbladder occurrence, in a refined multilocus molecular phylogeny of the tribe. We surveyed three different configurations of stomachs in akodontines, according to the degree of development and location of the glandular epithelium; in addition, two minor variations of one of these types were described. Of the five major clades that integrate Akodontini, four are characterized by a single stomach morphology, while one clade exhibits two morphologies. Mapping stomach type on the phylogeny recovered two configurations for the most recent ancestor of Akodontini. A revised survey of gallbladder evidence also revealed overlooked congruencies. The observed stomach diversity and its arrangement in the phylogeny, along with additional morphological characters and the genetic diversity among the main clades, supports the necessity of changes in the current classification of the tribe. Recognition of subtribes or partitioning of Akodontini into several additional tribes of equal rank could be suitable options.

A new mouse of the Peromyscus maniculatus species complex (Cricetidae) from the highlands of central Mexico

The Peromyscus maniculatus species complex is one of the most widespread group of small mammals in North America. However, the taxonomy and phylogenetic relationships among its constituent taxa remain unclear. As part of a revision of Peromyscus specimens from the highlands of the Trans-Mexican Volcanic Belt in central Mexico, we identified five individuals collected in 1968 that differed externally from other Peromyscus specimens, although morphologically similar to P. labecula and P. melanotis, both latter in the P. maniculatus species complex. Based on cranial measurements and mitochondrial DNA sequences, we aimed to more accurately determine the phylogenetic relationships and the taxonomic status of these individuals. Molecular phylogenetic analyses showed that the specimens formed a monophyletic clade sister to the P. maniculatus species complex. Pairwise genetic distances between those specimens and other species within the P. maniculatus species complex were greater than 7.91%. In addition, morphological analyses clearly distinguished the test specimens from P. melanotis and P. labecula. Based on the results of our molecular and morphological analyses, we conclude that these specimens represent an undescribed species of the P. maniculatus species complex, which we describe herein.

Differential Expression in Testis and Liver Transcriptomes from Four Species of Peromyscus (Rodentia: Cricetidae)

The genus Peromyscus represents a rapidly diverged clade of Cricetid rodents that contains multiple cryptic species and has a propensity for morphologic conservation across its members. The unresolved relationships in previously proposed phylogenies reflect a suspected rapid adaptive radiation. To identify functional groups of genes that may be important in reproductive isolation in a reoccurring fashion across the Peromyscus phylogeny, liver and testis transcriptomes from four species (P. attwateri, P. boylii, P. leucopus, and P. maniculatus) were generated and differential expression (DE) tests were conducted. Taxa were selected to represent members diverged from a common ancestor: P. attwateri + P. boylii (clade A), and P. leucopus + P. maniculatus (clade B). Comparison of clades (A vs. B) suggested that 252 transcripts had significant DE in the liver data set, whereas significant DE was identified for 657 transcripts in the testis data set. Further, 45 genes had DE isoforms in the 657 testis transcripts and most of these functioned in major reproductive roles such as acrosome assembly, spermatogenesis, and cell cycle processes (meiosis). DE transcripts in the liver mapped to more broad gene ontology terms (metabolic processes, catabolic processes, response to chemical, and regulatory processes), and DE transcripts in the testis mapped to gene ontology terms associated with reproductive processes, such as meiosis, sperm motility, acrosome assembly, and sperm–egg fusion. These results suggest that a suite of genes that conduct similar functions in the testes may be responsible for the adaptive radiation events and potential reoccurring speciation of Peromyscus in terms of reproduction through varying expression levels.

The Melanin-Concentrating Hormone (MCH) System: A Tale of Two Peptides

The melanin-concentrating hormone (MCH) system is a robust integrator of exogenous and endogenous information, modulating arousal and energy balance in mammals. Its predominant function in teleosts, however, is to concentrate melanin in the scales, contributing to the adaptive color change observed in several teleost species. These contrasting functions resulted from a gene duplication that occurred after the teleost divergence, which resulted in the generation of two MCH-coding genes in this clade, which acquired distinctive sequences, distribution, and functions, examined in detail here. We also describe the distribution of MCH immunoreactivity and gene expression in a large number of species, in an attempt to identify its core elements. While initially originated as a periventricular peptide, with an intimate relationship with the third ventricle, multiple events of lateral migration occurred during evolution, making the ventrolateral and dorsolateral hypothalamus the predominant sites of MCH in teleosts and mammals, respectively. Substantial differences between species can be identified, likely reflecting differences in habitat and behavior. This observation aligns well with the idea that MCH is a major integrator of internal and external information, ensuring an appropriate response to ensure the organism’s homeostasis. New studies on the MCH system in species that have not yet been investigated will help us understand more precisely how these habitat changes are connected to the hypothalamic neurochemical circuits, paving the way to new intervention strategies that may be used with pharmacological purposes.

A centromere satellite concomitant with extensive karyotypic diversity across the Peromyscus genus defies predictions of molecular drive

A common feature of eukaryotic centromeres is the presence of large tracts of tandemly arranged repeats, known as satellite DNA. However, these centromeric repeats appear to experience rapid evolution under forces such as molecular drive and centromere drive, seemingly without consequence to the integrity of the centromere. Moreover, blocks of heterochromatin within the karyotype, including the centromere, are hotspots for chromosome rearrangements that may drive speciation events by contributing to reproductive isolation. However, the relationship between the evolution of heterochromatic sequences and the karyotypic dynamics of these regions remains largely unknown. Here, we show that a single conserved satellite DNA sequence in the order Rodentia of the genus Peromyscus localizes to recurrent sites of chromosome rearrangements and heterochromatic amplifications. Peromyscine species display several unique features of chromosome evolution compared to other Rodentia, including stable maintenance of a strict chromosome number of 48 among all known species in the absence of any detectable interchromosomal rearrangements. Rather, the diverse karyotypes of Peromyscine species are due to intrachromosomal variation in blocks of repeated DNA content. Despite wide variation in the copy number and location of repeat blocks among different species, we find that a single satellite monomer maintains a conserved sequence and homogenized tandem repeat structure, defying predictions of molecular drive. The conservation of this satellite monomer results in common, abundant, and large blocks of chromatin that are homologous among chromosomes within one species and among diverged species. Thus, such a conserved repeat may have facilitated the retention of polymorphic chromosome variants within individuals and intrachromosomal rearrangements between species-both factors that have previously been hypothesized to contribute towards the extremely wide range of ecological adaptations that this genus exhibits.

Evolution of the rodent Trim5 cluster is marked by divergent paralogous expansions and independent acquisitions of TrimCyp fusions

Evolution of cellular innate immune genes in response to viral threats represents a rich area of study for understanding complex events that shape mammalian genomes. One of these genes, TRIM5, is a retroviral restriction factor that mediates a post-entry block to infection. Previous studies on the genomic cluster that contains TRIM5 identified different patterns of gene amplification and the independent birth of CypA gene fusions in various primate species. However, the evolution of Trim5 in the largest order of mammals, Rodentia, remains poorly characterized. Here, we present an expansive phylogenetic and genomic analysis of the Trim5 cluster in rodents. Our findings reveal substantial evolutionary changes including gene amplifications, rearrangements, loss and fusion. We describe the first independent evolution of TrimCyp fusion genes in rodents. We show that the TrimCyp gene found in some Peromyscus species was acquired about 2 million years ago. When ectopically expressed, the P. maniculatus TRIMCyp shows anti-retroviral activity that is reversed by cyclosporine, but it does not activate Nf-κB or AP-1 promoters, unlike the primate TRIMCyps. These results describe a complex pattern of differential gene amplification in the Trim5 cluster of rodents and identify the first functional TrimCyp fusion gene outside of primates and tree shrews.

Melanin-concentrating hormone peptidergic system: Comparative morphology between muroid species

Melanin-concentrating hormone (MCH) is a conserved neuropeptide, predominantly located in the diencephalon of vertebrates, and associated with a wide range of functions. While functional studies have focused on the use of the traditional mouse laboratory model, critical gaps exist in our understanding of the morphology of the MCH system in this species. Even less is known about the nontraditional animal model Neotomodon alstoni (Mexican volcano mouse). A comparative morphological study among these rodents may, therefore, contribute to a better understanding of the evolution of the MCH peptidergic system. To this end, we employed diverse immunohistochemical protocols to identify key aspects of the MCH system, including its spatial relationship to another neurochemical population of the tuberal hypothalamus, the orexins. Three-dimensional (3D) reconstructions were also employed to convey a better sense of spatial distribution to these neurons. Our results show that the distribution of MCH neurons in all rodents studied follows a basic plan, but individual characteristics are found for each species, such as the preeminence of a periventricular group only in the rat, the lack of posterior groups in the mouse, and the extensive presence of MCH neurons in the anterior hypothalamic area of Neotomodon. Taken together, these data suggest a strong anatomical substrate for previously described functions of the MCH system, and that particular neurochemical and morphological features may have been determinant to species-specific phenotypes in rodent evolution.

Spatial incongruence in the species richness and functional diversity of cricetid rodents

Biodiversity is multidimensional and different mechanisms can influence different dimensions. The spatial distribution of these dimensions can help in conservation decisions through the location of complementary areas with high diversity. We analyzed congruence in spatial patterns of species richness and functional diversity of cricetid rodents in the state of Oaxaca, southern Mexico, at different scales, and environmental variables related. Potential distribution models were produced for 49 species of cricetids in Maxent and superimposed to obtain potential communities in cells of 25, 50,100, 200 and 400 km2. We estimated species richness (SR) and functional diversity (SES.FD) eliminating the species richness effect through null models. The patterns and spatial congruence of species richness and functional diversity are described. The relationships between the environmental variables (elevation, temperature, precipitation, net primary productivity and potential evapotranspiration) and the SR and SES.FD were explored using Generalized Linear Models (GLMs) and Generalized Additive Models (GAMs). The highest species richness was found in mountainous ecosystems while the highest functional diversity was in tropical forests, revealing a spatial incongruence among these components of biodiversity (r = -0.14, p = 0.42; Pearson correlation). The locations of the cells of low congruence varied according to spatial resolution. In univariate models, elevation was the variable that best explained species richness (R2 = 0.77). No single variable explained the functional diversity; however, the models that included multiple environmental variables partially explained both the high and low functional diversity. The different patterns suggest that different historic, ecological and environmental processes could be responsible for the community structure of cricetid rodents in Oaxaca. These results indicate that one great challenge to be met to achieve more effective planning for biological conservation is to integrate knowledge regarding the spatial distribution of different dimensions of biodiversity.

Maladaptive phenotypic plasticity in cardiac muscle growth is suppressed in high-altitude deer mice

How often phenotypic plasticity acts to promote or inhibit adaptive evolution is an ongoing debate among biologists. Recent work suggests that adaptive phenotypic plasticity promotes evolutionary divergence, though several studies have also suggested that maladaptive plasticity can potentiate adaptation. The role of phenotypic plasticity, adaptive, or maladaptive, in evolutionary divergence remains controversial. We examined the role of plasticity in evolutionary divergence between two species of Peromyscus mice that differ in native elevations. We used cardiac mass as a model phenotype, since ancestral hypoxia-induced responses of the heart may be both adaptive and maladaptive at high-altitude. While left ventricle growth should enhance oxygen delivery to tissues, hypertrophy of the right ventricle can lead to heart failure and death. We compared left- and right-ventricle plasticity in response to hypoxia between captive-bred P. leucopus (representing the ancestral lowland condition) and P. maniculatus from high-altitude. We found that maladaptive ancestral plasticity in right ventricle hypertrophy is reduced in high-altitude deer mice. Analysis of the heart transcriptome suggests that changes in expression of inflammatory signaling genes, particularly interferon regulatory factors, contribute to the suppression of right ventricle hypertrophy. We found weak evidence that adaptive plasticity of left ventricle mass contributes to evolution. Our results suggest that selection to suppress ancestral maladaptive plasticity plays a role in adaptation.

Molecular systematics and phylogeography of the endemic Osgood's deermouse Osgoodomys banderanus (Rodentia: Cricetidae) in the lowlands of western Mexico

Osgoodomys banderanus is a recognized and endemic rodent species of western Mexico, an area known for its high biodiversity and number of endemisms. Phylogeographical relationships within this taxon were analyzed based on mitochondrial (ND3, tRNA-Arginine, ND4L and partial ND4) and nuclear (GHR) nucleotide sequences. We obtained a total of 112 samples from 22 localities, covering the complete distribution of the species. Phylogenetic analyses using Maximum Likelihood and Bayesian inference confirmed that Osgoodomys is a monophyletic group. In addition, phylogenetic and phylogeographic analyses detected three major clades, which do not coincide with the recognized subspecies of O. banderanus. The genetic lineages detected are the western clade (Nayarit, Jalisco and northern Colima), the central clade (Colima, Michoacán, and northern Guerrero) and the eastern clade (central and southern Guerrero). Genetic distances among clades (5-9%) and nucleotide substitutions (30-88) among haplogroups were high, especially in the southern group. Mountain ranges such as the Transmexican Volcanic Belt and the Sierra Madre del Sur, as well as the Balsas River act as geographical barriers for Osgoodomys. Our results suggest the presence of three independent species, which need to be characterized morphologically to confirm our hypothesis.

Sexual imprinting and speciation between two Peromyscus species

Sexual isolation, a reproductive barrier, can prevent interbreeding between diverging populations or species. Sexual isolation can have a clear genetic basis; however, it may also result from learned mate preferences that form via sexual imprinting. Here, we demonstrate that two sympatric species of mice-the white-footed mouse (Peromyscus leucopus) and its sister species, the cotton mouse (P. gossypinus)-hybridize only rarely in the wild despite co-occurrence in the same habitat and lack of any measurable intrinsic postzygotic barriers in laboratory crosses. We present evidence that strong conspecific mating preferences in each species result in significant sexual isolation. We find that these preferences are learned in at least one species: P. gossypinus sexually imprints on its parents, but in P. leucopus, additional factors influence mating preferences. Our study demonstrates that sexual imprinting contributes to reproductive isolation that reduces hybridization between otherwise interfertile species, supporting the role for learning in mammalian speciation.

Microbial communities exhibit host species distinguishability and phylosymbiosis along the length of the gastrointestinal tract

Host-associated microbial communities consist of stable and transient members that can assemble through purely stochastic processes associated with the environment or by interactions with the host. Phylosymbiosis predicts that if host-microbiota interactions impact assembly patterns, then one conceivable outcome is concordance between host evolutionary histories (phylogeny) and the ecological similarities in microbial community structures (microbiota dendrogram). This assembly pattern has been demonstrated in several clades of animal hosts in laboratory and natural populations, but in vertebrates, it has only been investigated using samples from faeces or the distal colon. Here, we collected the contents of five gut regions from seven rodent species and inventoried the bacterial communities by sequencing the 16S rRNA gene. We investigated how community structures varied across gut regions and whether the pattern of phylosymbiosis was present along the length of the gut. Gut communities varied by host species and gut region, with Oscillospira and Ruminococcus being more abundant in the stomach and hindgut regions. Gut microbial communities were highly distinguishable by host species across all gut regions, with the strength of the discrimination increasing along the length of the gut. Last, the pattern of phylosymbiosis was found in all five gut regions, as well as faeces. Aspects of the gut environment, such as oxygen levels, production of antimicrobials or other factors, may shift microbial communities across gut regions. However, regardless of these differences, host species maintain distinguishable, phylosymbiotic assemblages of microbes that may have functional impacts for the host.

Pleistocene climatic oscillations in Neotropical open areas: Refuge isolation in the rodent Oxymycterus nasutus endemic to grasslands

Pleistocene climatic oscillations favoured the expansion of grassland ecosystems and open vegetation landscapes throughout the Neotropics, and influenced the evolutionary history of species adapted to such environments. In this study, we sampled populations of the rodent Oxymycterus nasutus endemic to open areas in the Pampas and Atlantic Forest biomes to assess the tempo and mode of population divergence using an integrative approach, including coalescence theory, ecological niche models, and morphometry. Our results indicated that these O. nasutus populations exhibited high levels of genetic structure. Six major mtDNA clades were found, structuring these biomes into distinct groups. Estimates of their divergence times was indicated to be 0.571 myr. The high degree of genetic structure is reflected in the analyses of geometric morphometric; skull differences between lineages in the two ecoregions were detected. During the last glacial maximum, there was a strong increase in suitable abiotic conditions for O. nasutus. Distinct molecular markers revealed a population expansion over time, with a possible demographic retraction during the post-glacial period. Considering that all clades coalesce with the last interglacial maximum, our results indicated that reduction in suitable conditions during this period may have resulted in a possible vicariance associated with refuge isolation.

Muroid rodent phylogenetics: 900-species tree reveals increasing diversification rates

We combined new sequence data for more than 300 muroid rodent species with our previously published sequences for up to five nuclear and one mitochondrial genes to generate the most widely and densely sampled hypothesis of evolutionary relationships across Muroidea. An exhaustive screening procedure for publically available sequences was implemented to avoid the propagation of taxonomic errors that are common to supermatrix studies. The combined data set of carefully screened sequences derived from all available sequences on GenBank with our new data resulted in a robust maximum likelihood phylogeny for 900 of the approximately 1,620 muroids. Several regions that were equivocally resolved in previous studies are now more decisively resolved, and we estimated a chronogram using 28 fossil calibrations for the most integrated age and topological estimates to date. The results were used to update muroid classification and highlight questions needing additional data. We also compared the results of multigene supermatrix studies like this one with the principal published supertrees and concluded that the latter are unreliable for any comparative study in muroids. In addition, we explored diversification patterns as an explanation for why muroid rodents represent one of the most species-rich groups of mammals by detecting evidence for increasing net diversification rates through time across the muroid tree. We suggest the observation of increasing rates may be due to a combination of parallel increases in rate across clades and high average extinction rates. Five increased diversification-rate-shifts were inferred, suggesting that multiple, but perhaps not independent, events have led to the remarkable species diversity in the superfamily. Our results provide a phylogenetic framework for comparative studies that is not highly dependent upon the signal from any one gene.

The ultimate and proximate mechanisms driving the evolution of long tails in forest deer mice

Understanding both the role of selection in driving phenotypic change and its underlying genetic basis remain major challenges in evolutionary biology. Here, we use modern tools to revisit a classic system of local adaptation in the North American deer mouse, Peromyscus maniculatus, which occupies two main habitat types: prairie and forest. Using historical collections, we find that forest‐dwelling mice have longer tails than those from nonforested habitat, even when we account for individual and population relatedness. Using genome‐wide SNP data, we show that mice from forested habitats in the eastern and western parts of their range form separate clades, suggesting that increased tail length evolved independently. We find that forest mice in the east and west have both more and longer caudal vertebrae, but not trunk vertebrae, than nearby prairie forms. By intercrossing prairie and forest mice, we show that the number and length of caudal vertebrae are not correlated in this recombinant population, indicating that variation in these traits is controlled by separate genetic loci. Together, these results demonstrate convergent evolution of the long‐tailed forest phenotype through two distinct genetic mechanisms, affecting number and length of vertebrae, and suggest that these morphological changes—either independently or together—are adaptive.

Phylosymbiosis: Relationships and Functional Effects of Microbial Communities across Host Evolutionary History

Phylosymbiosis was recently proposed to describe the eco-evolutionary pattern, whereby the ecological relatedness of host-associated microbial communities parallels the phylogeny of related host species. Here, we test the prevalence of phylosymbiosis and its functional significance under highly controlled conditions by characterizing the microbiota of 24 animal species from four different groups (Peromyscus deer mice, Drosophila flies, mosquitoes, and Nasonia wasps), and we reevaluate the phylosymbiotic relationships of seven species of wild hominids. We demonstrate three key findings. First, intraspecific microbiota variation is consistently less than interspecific microbiota variation, and microbiota-based models predict host species origin with high accuracy across the dataset. Interestingly, the age of host clade divergence positively associates with the degree of microbial community distinguishability between species within the host clades, spanning recent host speciation events (~1 million y ago) to more distantly related host genera (~108 million y ago). Second, topological congruence analyses of each group's complete phylogeny and microbiota dendrogram reveal significant degrees of phylosymbiosis, irrespective of host clade age or taxonomy. Third, consistent with selection on host–microbiota interactions driving phylosymbiosis, there are survival and performance reductions when interspecific microbiota transplants are conducted between closely related and divergent host species pairs. Overall, these findings indicate that the composition and functional effects of an animal's microbial community can be closely allied with host evolution, even across wide-ranging timescales and diverse animal systems reared under controlled conditions.

Studies on the assembly and function of host-microbiota symbioses are inherently complicated by the diverse effects of diet, age, sex, host genetics, and endosymbionts. Central to unraveling one effect from the other is an experimental framework that reduces confounders. Using common rearing conditions across four animal groups (deer mice, flies, mosquitoes, and wasps) that span recent host speciation events to more distantly related host genera, this study tests whether microbial community assembly is generally random with respect to host relatedness or "phylosymbiotic," in which the phylogeny of the host group is congruent with ecological relationships of their microbial communities. Across all four animal groups and one external dataset of great apes, we apply several statistics for analyzing congruencies and demonstrate phylosymbiosis to varying degrees in each group. Moreover, consistent with selection on host–microbiota interactions driving phylosymbiosis, transplanting interspecific microbial communities in mice significantly decreased their ability to digest food. Similarly, wasps that received transplants of microbial communities from different wasp species had lower survival than those given their own microbiota. Overall, this experimental and statistical framework shows how microbial community assembly and functionality across related species can be linked to animal evolution, health, and survival.