Genomic resolution of cryptic species diversity in chipmunks

Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species diversity within an ecologically well-known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence-capture data for 21,551 autosomal and 493 X-linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow-pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population-level sequence data revealed that the so-called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow-pine chipmunk but does show evidence for recurrent hybridization with this and other species.

Diversification, Introgression, and Rampant Cytonuclear Discordance in Rocky Mountains Chipmunks (Sciuridae: Tamias)

Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon-capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (T. canipes, T. cinereicollis, T. dorsalis, T. quadrivittatus, T. rufus, and T. umbrinus) that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group.

Disentangling lousy relationships: Comparative phylogenomics of two sucking louse lineages parasitizing chipmunks

The evolution of obligate parasites is often interpreted in light of their hosts' evolutionary history. An expanded approach is to examine the histories of multiple lineages of parasites that inhabit similar environments on a particular host lineage. Western North American chipmunks (genus Tamias) have a broad distribution, a history of divergence with gene flow, and host two species of sucking lice (Anoplura), Hoplopleura arboricola and Neohaematopinus pacificus. From total genomic sequencing, we obtained sequences of over 1100 loci sampled across the genomes of these lice to compare their evolutionary histories and examine the roles of host association in structuring louse relationships. Within each louse species, clades are largely associated with closely related chipmunk host species. Exceptions to this pattern appear to have a biogeographic component, but differ between the two louse species. Phylogenetic relationships among these major louse clades, in both species, are not congruent with chipmunk relationships. In the context of host associations, each louse lineage has a different evolutionary history, supporting the hypothesis that host-parasite assemblages vary both across the landscape and with the taxa under investigation. In addition, the louse Hoplopleura erratica (parasitizing the eastern Tamias striatus) is embedded within H. arboricola, rendering it paraphyletic. This phylogenetic result, together with comparable divergences within H. arboricola, indicate a need for taxonomic revision. Both host divergence and biogeographic components shape parasite diversification as demonstrated by the distinctive diversification patterns of these two independently evolving lineages that parasitize the same hosts.

No evidence for phylosymbiosis in western chipmunk species

Phylosymbiosis refers to a congruent pattern between the similarity of microbiomes of different species and the branching pattern of the host phylogeny. Phylosymbiosis has been detected in a variety of vertebrate and invertebrate hosts, but has only been assessed in geographically isolated populations. We tested for phylosymbiosis in eight (sub)species of western chipmunks with overlapping ranges and ecological niches; we used a nuclear (Acrosin) and a mitochondrial (CYTB) phylogenetic marker because there are many instances of mitochondrial introgression in chipmunks. We predicted that similarity among microbiomes increases with: (1) increasing host mitochondrial relatedness, (2) increasing host nuclear genome relatedness and (3) decreasing geographic distance among hosts. We did not find statistical evidence supporting phylosymbiosis in western chipmunks. Furthermore, in contrast to studies of other mammalian microbiomes, similarity of chipmunk microbiomes is not predominantly determined by host species. Sampling site explained most variation in microbiome composition, indicating an important role of local environment in shaping microbiomes. Fecal microbiomes of chipmunks were dominated by Bacteroidetes (72.2%), followed by Firmicutes (24.5%), which is one of the highest abundances of Bacteroidetes detected in wild mammals. Future work will need to elucidate the effects of habitat, ecology and host genomics on chipmunk microbiomes.

Directional selection effects on patterns of phenotypic (co)variation in wild populations

Phenotypic (co)variation is a prerequisite for evolutionary change, and understanding how (co)variation evolves is of crucial importance to the biological sciences. Theoretical models predict that under directional selection, phenotypic (co)variation should evolve in step with the underlying adaptive landscape, increasing the degree of correlation among co-selected traits as well as the amount of genetic variance in the direction of selection. Whether either of these outcomes occurs in natural populations is an open question and thus an important gap in evolutionary theory. Here, we documented changes in the phenotypic (co)variation structure in two separate natural populations in each of two chipmunk species (Tamias alpinus and T. speciosus) undergoing directional selection. In populations where selection was strongest (those of T. alpinus), we observed changes, at least for one population, in phenotypic (co)variation that matched theoretical expectations, namely an increase of both phenotypic integration and (co)variance in the direction of selection and a re-alignment of the major axis of variation with the selection gradient.

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.

Morphological and dietary responses of chipmunks to a century of climate change

Predicting how individual taxa will respond to climatic change is challenging, in part because the impacts of environmental conditions can vary markedly, even among closely related species. Studies of chipmunks (Tamias spp.) in Yosemite National Park provide an important opportunity to explore the reasons for this variation in response. While the alpine chipmunk (T. alpinus) has undergone a significant elevational range contraction over the past century, the congeneric and partially sympatric lodgepole chipmunk (T. speciosus) has not experienced an elevational range shift during this period. As a first step toward identifying the factors underlying this difference in response, we examined evidence for dietary changes and changes in cranial morphology in these species over the past century. Stable isotope analyses of fur samples from modern and historical museum specimens of these species collected at the same localities indicated that signatures of dietary change were more pronounced in T. alpinus, although diet breadth did not differ consistently between the study species. Morphometric analyses of crania from these specimens revealed significant changes in cranial shape for T. alpinus, with less pronounced changes in shape for T. speciosus; evidence of selection on skull morphology was detected for T. alpinus, but not for T. speciosus. These results are consistent with growing evidence that T. alpinus is generally more responsive to environmental change than T. speciosus, but emphasize the complex and often geographically variable nature of such responses. Accordingly, future studies that make use of the taxonomically and spatially integrative approach employed here may prove particularly informative regarding relationships between environmental conditions, range changes, and patterns of phenotypic variation.

Unlocking the vault: next-generation museum population genomics

Natural history museum collections provide unique resources for understanding how species respond to environmental change, including the abrupt, anthropogenic climate change of the past century. Ideally, researchers would conduct genome-scale screening of museum specimens to explore the evolutionary consequences of environmental changes, but to date such analyses have been severely limited by the numerous challenges of working with the highly degraded DNA typical of historic samples. Here, we circumvent these challenges by using custom, multiplexed, exon capture to enrich and sequence ~11,000 exons (~4 Mb) from early 20th-century museum skins. We used this approach to test for changes in genomic diversity accompanying a climate-related range retraction in the alpine chipmunks (Tamias alpinus) in the high Sierra Nevada area of California, USA. We developed robust bioinformatic pipelines that rigorously detect and filter out base misincorporations in DNA derived from skins, most of which likely resulted from postmortem damage. Furthermore, to accommodate genotyping uncertainties associated with low-medium coverage data, we applied a recently developed probabilistic method to call single-nucleotide polymorphisms and estimate allele frequencies and the joint site frequency spectrum. Our results show increased genetic subdivision following range retraction, but no change in overall genetic diversity at either nonsynonymous or synonymous sites. This case study showcases the advantages of integrating emerging genomic and statistical tools in museum collection-based population genomic applications. Such technical advances greatly enhance the value of museum collections, even where a pre-existing reference is lacking and points to a broad range of potential applications in evolutionary and conservation biology.