Phylogeographic study of the California vole, Microtus californicus

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co-evolution of microbes and their hosts as well as current host-specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial-nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome-wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β-diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among-population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).

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.

Isotope values of California vole (Microtus californicus) hair relate to historical drought and land use patterns in California, USA

Increased drought frequency and intensity and agricultural intensification have been key stressors to ecological systems over the past century. Biological proxies (e.g., pollen, tree rings) have been used to track this environmental change; however, linking these changes to the ecology of organisms remains challenging. Here, we link historical drought records to conditions of high water-stress in grassland habitats through the stable isotope analysis of California vole museum specimens (Microtus californicus). Using museum collections spanning 118-years (1891-2009), isotope values of dated hair tissues were associated with statewide drought metrics on the Palmer Drought Severity Index. We observed a positive correlation between δ¹⁵N and δ¹⁸O values and drought severity. The range in δ¹⁵N values (~ 18‰) is greater than what would be expected as a result of dietary shifts across the landscape (~ 3‰), and is likely attributed to the combined effects of physiological responses of M. californicus and isotopic shifts in plant resources with increased water-stress. Geospatial patterns in δ³⁴S values of hair tissues reflect higher baseline isotope values in coastal habitats. However, comparably high δ³⁴S values in the southern-most inland localities suggest sulfur fertilization of croplands and subsequent transfer to surrounding grassland habitats in ³⁴S enriched forms. A broad δ¹³C range (- 28.7 to - 14.3‰) further suggests the consumption of C₃ and C₄ plant-based dietary proteins. As shown here, stable isotope analysis of museum collections can provide a climate and land use record based on the physiological performance and ecology of a study species in a region affected intensely by anthropogenic activities.

Mito-nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi-permeable. Using mitochondrial cytochrome b sequences and genome-wide loci developed via ddRAD-seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5-54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

Comprehensive cross-amplification of microsatellite multiplex sets across the rodent genus Microtus

We developed four multiplex panels comprising 19 microsatellite loci and tested their amplification in 21 rodent species important for agricultural and conservation management (Microtus, Arvicola, Chionomys). On average, 17.6 loci amplified per species. Number of alleles ranged from 1 to 19 per locus. We report an additional locus polymorphic in 15 vole species.