Phylogeography of the western jumping mouse (Zapus princeps) detects deep and persistent allopatry with expansion

Phylogeography within the Peromyscus maniculatus species group: Understanding past distribution of genetic diversity and areas of refugia in western North America

The effects of anthropogenic climate change on biodiversity have been recognized on every continent, ocean, and across different taxonomic groups. Here, we study the range dynamics and demography of a cosmopolitan species: the deer mouse, Peromyscus maniculatus. We generated a multilocus SNP dataset using the ddRADseq protocol for 218 individuals across the geographic range within three western North American lineages of this species group. We evaluated population structure using several methods and explored the correlation between geographic and genetic distances. We modeled the demographic history using a site frequency spectrum approach and used a machine learning algorithm to infer current and past (Last Glacial Maximum; LGM) environmental suitability. Lastly, we explored the origin of population expansion for the identified lineages. The genome-wide SNP dataset was able to identify-three regionally distinct groups- 1) P. m. gambelii (southern California); 2) P. keeni (Pacific Northwest); 3) P. m. sonoriensis (a broad population spanning the Pacific Northwest through central California and across the Rocky Mountains into the Great Plains). Demographic analysis indicated the splits between the three populations occurred within the last 500 thousand years, with one very recent (late Holocene) split. Ecological niche models for each of these lineages predicted suitable environment present throughout their known ranges for current conditions, and a severe reduction of northern habitat in the past. The deer mouse has responded to past climate changes by expanding its range during interglacial periods and contracting its range during glacial periods leading to strong population differentiation. But lower magnitude climate change or other processes within the Holocene interglacial period led to population differentiation as well, which is likely still ongoing today given the substantial anthropogenic climate change and other landscape transformations caused by humans during the Anthropocene. By understanding the historical processes that led to the contemporary geographic distribution of biodiversity, we can determine the relative importance of different factors that shape biodiversity, now and into the future.

A precarious future for distinctive peripheral populations of meadow voles (Microtus pennsylvanicus)

Conservation efforts rely on robust taxonomic assessments that should be based on critical assessment of interspecific boundaries, infraspecific variation, and potentially distinctive peripheral populations. The meadow vole (Microtus pennsylvanicus) is widely distributed across North America, including 28 morphologically defined subspecies and numerous isolated populations. Because some subspecies are of high conservation concern, we examined geographic variation across the range of the species to test existing infraspecific taxonomy in terms of local and regional diversification. We sequenced mitochondrial DNA (mtDNA) from 20 subspecies of M. pennsylvanicus and contextualized infraspecific variation through comparison of pairwise genetic distances derived from an extended data set of 63 species of Microtus. We found strong support for at least three divergent clades within M. pennsylvanicus, with observed intraspecific clade divergence exceeding that between several pairwise comparisons of sister species within Microtus. Six nuclear genes were then sequenced to test the validity of mtDNA structure and to further evaluate the possibility of cryptic, species-level diversity using Bayes factor species delimitation (BFD) analyses. BFD consistently and decisively supported multiple species based on the multilocus approach. We propose that taxonomic revision of the meadow vole is required, with the eastern clade now identified as M. pennsylvanicus (Ord 1815), the western clade as M. drummondii (Audubon and Bachman 1853), and the coastal Florida clade as M. dukecampbelli (Woods, Post, and Kilpatrick 1982). We suggest that such an arrangement would more closely reflect evolutionary history and provide critical context for further examination of distinctive southern peripheral populations that harbor novel evolutionary legacies and adaptive potential.

Molecular phylogeographic analyses and species delimitations reveal that Leopoldamys edwardsi (Rodentia: Muridae) is a species complex

Leopoldamys edwardsi is a species with wide distribution ranges in southern China but is not discussed in studies on geographic variation and species differentiation. We used 2 mitochondrial (Cytb, CO1) and 3 nuclear (GHR, IRBP and RAG1) genes to clarify species phylogeography and geographical differentiation. Maximum likelihood (ML) and Bayesian phylogenetic inference (BI) trees consistently indicated that L. edwardsi is a species complex containing 3 main lineages with high Kimura-2-parameter (K2P) divergences (i.e. lineages L_N , L_S and L_HN ) found in the northern and southern China and Hainan Island, respectively. The 3 species delimitation methods, automated barcoding gap discovery, Bayesian poisson tree process analysis and Bayesian phylogenetics and phylogeography, consistently supported the existence of cryptic species. Divergence times among the main lineages were inferred to be during the Pleistocene, with L_HN /L_S split at 1.33 Ma and L_N /(L_HN +L_S ) at 2.61 Ma; the diversifications of L. edwardsi complex might be caused by the rapid uplifts of Tibetan Plateau, paleoclimate change and complex topography. The divergence between L_HN and L_S was probably related to the separation of Hainan Island from the mainland via the formation of the Qiongzhou Strait. Lineages L_N and (L_S +L_HN ) likely diverged due to the Wuyi-Nanling mountain range forming a dispersal barrier. Our results suggested that L. edwardsi complex contains at least 3 distinct species: L_HN represents L. hainanensis, endemic to Hainan Island and previously considered as a subspecies L. e. hainanensis; L_S represents a cryptic species distributed throughout the southern Chinese continent; and L_N represents the nominotypical species L. edwardsi.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

Integrative species delimitation of the widespread North American jumping mice (Zapodinae)

Delimiting species can be challenging, but is a key step for the critical examination of evolutionary history and for prioritizing conservation efforts. Because systematic relationships are often determined iteratively using tests based on taxonomy, such methods can fail to detect cryptic variation and result in biased conclusions. Conversely, discovery-based approaches provide a powerful way to define operational taxonomic units and test species boundaries. We compare both approaches (taxonomy-based delimitation - TBD and discovery-based delimitation - DBD) within North American jumping mice (Zapodinae) using broad sampling, multilocus analyses, and ecological tests. This group diversified through the dynamic glacial-interglacial periods of the Quaternary and phylogeographic tests reveal 28 lineages that correspond poorly with current taxonomy (4 species, 32 nominal subspecies). However, neither the 4-species or 28-lineage hypotheses are optimal for species-level classification. Rather, information theoretic approaches (Bayes Factors) indicate a 15-species hypothesis is best for characterizing genetic variation in this group, with subsequent iterative pairwise ecological tests failing to confirm four species pairs. Taken together, evolutionary and ecological tests capture divergence among 11 putative species that, if upheld by additional tests, will lead to taxonomic revision and reevaluation of conservation plans.

Temporal and spatial mosaics: deep host association and shallow geographic drivers shape genetic structure in a widespread pinworm, Rauschtineria eutamii

Climate and host demographic cycling often shape both parasite genetic diversity and host distributions, processes that transcend a history of strict host-parasite association. We explored host associations and histories based on an evaluation of mitochondrial and nuclear sequences to reveal the underlying history and genetic structure of a pinworm, Rauschtineria eutamii, infecting 10 species of western North American chipmunks (Rodentia:Tamias, subgenus Neotamias). Rauschtineria eutamii contains divergent lineages influenced by the diversity of hosts and variation across the complex topography of western North America. We recovered six reciprocally monophyletic R. eutamii mitochondrial clades, largely supported by nuclear gene trees, exhibiting divergence levels comparable to intraspecific variation reported for other nematodes. Phylogenetic relationships among pinworm clades suggest that R. eutamii colonized an ancestral lineage of western chipmunks and lineages persisted during historical isolation in diverging Neotamias species or species groups. Pinworm diversification, however, is incongruent and asynchronous relative to host diversification. Secondarily, patterns of shallow divergence were shaped by geography through events of episodic colonization reflecting an interaction of taxon pulses and ecological fitting among assemblages in recurrent sympatry. Pinworms occasionally infect geographically proximal host species; however, host switching may be unstable or ephemeral, as there is no signal of host switching in the deeper history of R. eutamii.

Variation in phenology of hibernation and reproduction in the endangered New Mexico meadow jumping mouse (Zapus hudsonius luteus)

Hibernation is a key life history feature that can impact many other crucial aspects of a species’ biology, such as its survival and reproduction. I examined the timing of hibernation and reproduction in the federally endangered New Mexico meadow jumping mouse (Zapus hudsonius luteus), which occurs across a broad range of latitudes and elevations in the American Southwest. Data from museum specimens and field studies supported predictions for later emergence and shorter active intervals in montane populations relative to lower elevation valley populations. A low-elevation population located at Bosque del Apache National Wildlife Refuge (BANWR) in the Rio Grande valley was most similar to other subspecies of Z. hudsonius: the first emergence date was in mid-May and there was an active interval of 162 days. In montane populations of Z. h. luteus, the date of first emergence was delayed until mid-June and the active interval was reduced to ca 124–135 days, similar to some populations of the western jumping mouse (Z. princeps). Last date of immergence into hibernation occurred at about the same time in all populations (mid to late October). In montane populations pregnant females are known from July to late August and evidence suggests that they have a single litter per year. At BANWR two peaks in reproduction were expected based on similarity of active season to Z. h. preblei. However, only one peak was clearly evident, possibly due to later first reproduction and possible torpor during late summer. At BANWR pregnant females are known from June and July. Due to the short activity season and geographic variation in phenology of key life history events of Z. h. luteus, recommendations are made for the appropriate timing for surveys for this endangered species.

Complex phylogeography and historical hybridization between sister taxa of freshwater sculpin (Cottus)

Species ranges that span different geographic landscapes frequently contain cryptic species- or population-level structure. Identifying these possible diversification factors can often be accomplished under a comparative phylogeographic framework. However, comparisons suffer if previous studies are limited to a particular group or habitat type. In California, a complex landscape has led to several phylogeographic breaks, primarily in terrestrial species. However, two sister taxa of freshwater fish, riffle sculpin (Cottus gulosus) and Pit sculpin (Cottus pitensis), display ranges based on morphological identifications that do not coincide with these breaks. Using a comprehensive sampling and nuclear, mitochondrial and microsatellite markers, we hypothesized that proposed species ranges are erroneous based on potential hybridization/gene flow between species. Results identified a phylogeographic signature consistent with this hypothesis, with breaks at the Coast Range Mountains and Sacramento/San Joaquin River confluence. Coastal locations of C. gulosus represent a unique lineage, and 'true' C. gulosus were limited to the San Joaquin basin, both regions under strong anthropogenic influence and potential conservation targets. C. pitensis limits extended historically throughout the Sacramento/Pit River basin but currently are restricted to the Pit River. Interestingly, locations in the Sacramento River contained low levels of ancestral hybridization and gene flow from C. gulosus but now appear to be a distinct population. The remaining population structure was strongly correlated with Sierra Nevada presence (high) or absence (low). This study stresses the importance of testing phylogeographic breaks across multiple taxa/habitats before conservation decisions are made, but also the potential impact of different geographic landscapes on evolutionary diversification.