A recent origin of marmots

Ultraconserved elements improve resolution of marmot phylogeny and offer insights into biogeographic history

Marmots (Marmota spp.) comprise a lineage of large-bodied ground squirrels that diversified rapidly in the Pleistocene, when the planet quickly transitioned to a drier, colder, and highly seasonal climate-particularly at high latitudes. Fossil evidence indicates the genus spread from North America, across Beringia, and into the European Alps over the course of only a few million years, beginning in the late Pliocene. Marmots are highly adapted to survive long and severely cold winters, and this likely favored their expansion and diversification over this time period. Previous phylogenetic studies have identified two major subgenera of marmots, but the timing of important speciation events and some species relationships have been difficult to resolve. Here we use ultraconserved elements and mitogenomes, with samples from all 15 extant species, to more precisely retrace how and when marmots came to inhabit a vast Holarctic range. Our results indicate marmots arose in North America in the mid Miocene (∼16.3 Mya) and dispersed across the Bering Land Bridge in the late Pliocene (∼3-4 Mya); in addition, our fossil-calibrated timeline is suggestive of the rise and spread of open grasslands as being particularly important to marmot diversification. The woodchuck (M. monax) and the Alaska marmot (M. broweri) are found to be more closely related to the Eurasian species than to the other North American species. Paraphyly is evident in the bobak marmot (M. bobak) and the hoary marmot (M. caligata), and in the case of the latter the data are highly suggestive of a second, cryptic species in the Cascade Mountains of Washington.

Gene expression shifts in yellow-bellied marmots prior to natal dispersal

The causes and consequences of vertebrate natal dispersal have been studied extensively, yet little is known about the molecular mechanisms involved. We used RNA-seq to quantify transcriptomic gene expression in blood of wild yellow-bellied marmots (Marmota flaviventer) prior to dispersing from or remaining philopatric to their natal colony. We tested 3 predictions. First, we hypothesized dispersers and residents will differentially express genes and gene networks since dispersal is physiologically demanding. Second, we expected differentially expressed genes to be involved in metabolism, circadian processes, and immune function. Finally, in dispersing individuals, we predicted differentially expressed genes would change as a function of sampling date relative to dispersal date. We detected 150 differentially expressed genes, including genes that have critical roles in lipid metabolism and antigen defense. Gene network analysis revealed a module of 126 coexpressed genes associated with dispersal that was enriched for extracellular immune function. Of the dispersal-associated genes, 22 altered expression as a function of days until dispersal, suggesting that dispersal-associated genes do not initiate transcription on the same time scale. Our results provide novel insights into the fundamental molecular changes required for dispersal and suggest evolutionary conservation of functional pathways during this behavioral process.

Hypoxic and Cold Adaptation Insights from the Himalayan Marmot Genome

The Himalayan marmot (Marmota himalayana) is a hibernating mammal that inhabits the high-elevation regions of the Himalayan mountains. Here we present a draft genome of the Himalayan marmot, with a total assembly length of 2.47 Gb. Phylogenetic analyses showed that the Himalayan marmot diverged from the Mongolian marmot approximately 1.98 million years ago. Transcriptional changes during hibernation included genes responsible for fatty acid metabolism in liver and genes involved in complement and coagulation cascades and stem cell pluripotency pathways in brain. Two selective sweep genes, Slc25a14 and ψAamp, showed apparent genotyping differences between low- and high-altitude populations. As a processed pseudogene, ψAamp may be biologically active to influence the stability of Aamp through competitive microRNA binding. These findings shed light on the molecular and genetic basis underlying adaptation to extreme environments in the Himalayan marmot.

Insights into the evolution of pathogenicity of Escherichia coli from genomic analysis of intestinal E. coli of Marmota himalayana in Qinghai-Tibet plateau of China

Escherichia coli is both of a widespread harmless gut commensal and a versatile pathogen of humans. Domestic animals are a well-known reservoir for pathogenic E. coli. However, studies of E. coli populations from wild animals that have been separated from human activities had been very limited. Here we obtained 580 isolates from intestinal contents of 116 wild Marmot Marmota himalayana from Qinghai–Tibet plateau, China, with five isolates per animal. We selected 125 (hereinafter referred to as strains) from the 580 isolates for genome sequencing, based on unique pulse field gel electrophoresis patterns and at least one isolate per animal. Whole genome sequence analysis revealed that all 125 strains carried at least one and the majority (79.2%) carried multiple virulence genes based on the analysis of 22 selected virulence genes. In particular, the majority of the strains carried virulence genes from different pathovars as potential 'hybrid pathogens'. The alleles of eight virulence genes from the Marmot E. coli were found to have diverged earlier than all known alleles from human and other animal E. coli. Phylogenetic analysis of the 125 Marmot E. coli genomes and 355 genomes selected from 1622 human and other E. coli strains identified two new phylogroups, G and H, both of which diverged earlier than the other phylogroups. Eight of the 12 well-known pathogenic E. coli lineages were found to share a most recent common ancestor with one or more Marmot E. coli strains. Our results suggested that the intestinal E. coli of the Marmots contained a diverse virulence gene pool and is potentially pathogenic to humans. These findings provided a new understanding of the evolutionary origin of pathogenic E. coli.

Complex history of isolation and gene flow in hoary, Olympic, and endangered Vancouver Island marmots

Climate change resulting in a reduction of alpine habitat is believed to pose a considerable risk to alpine-dependent species, including many marmots. Hoary marmots (Marmota caligata) range throughout much of the mountainous Pacific Northwest (PNW) and Rocky Mountains while the closely related Olympic and Vancouver Island marmots (M. olympus and M. vancouverensis, respectively) are restricted to small isolated regions of the PNW. The endemic Vancouver Island marmot is currently classified as Critically Endangered and the Olympic marmot has recently experienced dramatic population declines. Previous phylogenetic studies of PNW marmot species have had limited power as they focused on resolving interspecific relationships, implicitly assumed an absence of gene flow among currently recognized species, included relatively few individuals, and relied heavily or entirely on mitochondrial DNA. We sequenced 2 mitochondrial and 4 nuclear markers from 167 hoary, 4 Vancouver Island, and 5 Olympic marmots in order to investigate phylogenetic relationships and historic gene flow among these species. We recovered 2 monophyletic (and predominantly allopatric) mitochondrial clades of hoary marmots that are not sister groups. Instead, Vancouver Island marmots formed a monophyletic mitochondrial sister clade to 1 of the hoary marmot clades. Nuclear loci did not recover the 2 mitochondrial clades of hoary marmots and suggest that Vancouver Island marmots may have experienced mitochondrial introgression from coastal mainland hoary marmots. Additionally, our nuclear results suggest possible gene flow between hoary and Olympic marmots despite different chromosomal formulas. Rather than resolving what has previously been considered a straightforward 3-taxon phylogenetic question, our findings suggest a complicated history of rapid divergence of the 3 species followed by intermittent and possibly ongoing gene flow between hoary marmots and both Olympic and Vancouver Island marmots. These results therefore have significant implications for the conservation of the latter 2 species, both of which are conservation concerns.

Does sociality drive the evolution of communicative complexity? A comparative test with ground-dwelling sciurid alarm calls

While sociality has been hypothesized to drive the evolution of communicative complexity, the relationship remains to be formally tested. We derive a continuous measure of social complexity from demographic data and use this variable to explain variation in alarm repertoire size in ground-dwelling sciurid rodents (marmots, Marmota spp.; prairie dogs, Cynomys spp.; and ground squirrels, Spermophilus spp.). About 40% of the variation in alarm call repertoire size was explained by social complexity in the raw data set. To determine the degree to which this relationship may have been influenced by historical relationships between species, we used five different phylogenetic hypotheses to calculate phylogenetically independent contrasts. Less variation was significantly explained in contrast-based analyses, but a general positive relationship remained. Social complexity explained more variation in alarm call repertoire size in marmots, while sociality explained no variation in repertoire size in prairie dogs and no variation in phylogenetically based analyses of squirrels. In most cases, substantial variation remained unexplained by social complexity. We acknowledge that factors other than social complexity, per se, may contribute to the evolution of alarm call repertoire size in sciurid rodents, and we discuss alternative hypotheses. Our measure of social complexity could be used by other researchers to test explicit evolutionary hypotheses that involve social complexity.

Phylogeographical structure and regional history of the dusky-footed woodrat, Neotoma fuscipes

The dusky-footed woodrat, Neotoma fuscipes, is a medium-sized rodent that inhabits low elevation woodland habitats along the Pacific coast of North America from Oregon, throughout California and into Baja California. Analyses of mitochondrial sequence variation throughout the distribution reveal substantial phylogeographical structure within N. fuscipes. The major mitochondrial lineages are largely concordant with previously identified morphological subdivisions within the taxon. The geographical distribution of distinct clades suggests that a combination of topographic barriers and the expansion and contraction of suitable habitat during the past 2 million years, especially along particular mountain ranges, have played a major role in the diversification of N. fuscipes. Furthermore, relatively low levels of genetic variation across the northern half of the distribution suggest that dusky-footed woodrats may have only recently expanded into this region.

Molecular phylogeny of the chipmunks inferred from Mitochondrial cytochrome b and cytochrome oxidase II gene sequences

There are currently 25 recognized species of the chipmunk genus Tamias. In this study we sequenced the complete mitochondrial cytochrome b (cyt b) gene of 23 Tamias species. We analyzed the cyt b sequence and then analyzed a combined data set of cyt b along with a previous data set of cytochrome oxidase subunit II (COII) sequence. Maximum-likelihood was used to further test the fit of models of evolution to the cyt b data. Other sciurid cyt b sequence was added to examine the evolution of Tamias in the context of other sciurids. Relationships among Tamias species are discussed, particularly the possibility of a current sorting event among taxa of the southwestern United States and the extreme divergences among the three subgenera (Neotamias, Eutamias, and Tamias).

Molecular phylogeny of European muroid rodents based on complete cytochrome b sequences

Phylogenetic relationships among 18 species of mainly European muroid rodents that belong to three subfamilies were estimated using complete sequences of the mitochondrial cytochrome b gene. The inferred monophyly of the subfamilies Murinae (mice and rats) and Arvicolinae (voles, lemmings, and muskrats) is in agreement with previous studies. Within the Murinae, the morphology-based division of the genus Apodemus into three subgenera is supported by these DNA sequence data. The relationships among the different genera of the Murinae were generally poorly resolved, and the relationships of Micromys and Acomys to the other murine genera remained unresolved. Within the subfamily Arvicolinae, the relations of the genera Arvicola, Clethrionomys, and Microtus remained tentative with our data. However, within the Microtus group, there is a good molecular support for the phylogenetic relationships. These findings suggest that the origin of the different murine and arvicoline lineages was rapid, indicating an adaptive radiation with fast speciation.

A test of the acoustic adaptation hypothesis in four species of marmots

Acoustic signals must be transmitted from a signaller to a receiver during which time they become modified. The acoustic adaptation hypothesis suggests that selection should shape the structure of long-distance signals to maximize transmission through different habitats. A specific prediction of the acoustic adaptation hypothesis is that long-distance signals of animals in their native habitat are expected to change less during transmission than non-native signals within that habitat. This prediction was tested using the alarm calls of four species of marmots that live in acoustically different habitats and produce species-specific, long-distance alarm vocalizations: yellow-bellied marmot, Marmota flaviventris; Olympic marmot, M. olympus; hoary marmot, M. caligata; and woodchuck, M. monax. By doing so, we evaluated the relative importance the acoustic environment plays on selecting for divergent marmot alarm calls. Representative alarm calls of the four species were broadcast and rerecorded in each species' habitat at four distances from a source. Rerecorded, and therefore degraded alarm calls, were compared to undegraded calls using spectrogram correlation. If each species' alarm call was transmitted with less overall degradation in its own environment, a significant interaction between species' habitat and species' call type would be expected. Transmission fidelity at each of four distances was treated as a multivariate response and differences among habitat and call type were tested in a two-way MANOVA. Although significant overall differences in the transmission properties of the habitats were found, and significant overall differences in the transmission properties of the call types were found, there was no significant interaction between habitat and call type. Thus, the evidence did not support the acoustic adaptation hypothesis for these marmot species. Factors other than maximizing long-distance transmission through the environment may be important in the evolution of species-specific marmot alarm calls. (c) 1998 The Association for the Study of Animal Behaviour.

Mammalian mitochondrial DNA evolution: a comparison of the cytochrome b and cytochrome c oxidase II genes

The evolution of two mitochondrial genes, cytochrome b and cytochrome c oxidase subunit II, was examined in several eutherian mammal orders, with special emphasis on the orders Artiodactyla and Rodentia. When analyzed using both maximum parsimony, with either equal or unequal character weighting, and neighbor joining, neither gene performed with a high degree of consistency in terms of the phylogenetic hypotheses supported. The phylogenetic inconsistencies observed for both these genes may be the result of several factors including differences in the rate of nucleotide substitution among particular lineages (especially between orders), base composition bias, transition/transversion bias, differences in codon usage, and different constraints and levels of homoplasy associated with first, second, and third codon positions. We discuss the implications of these findings for the molecular systematics of mammals, especially as they relate to recent hypotheses concerning the polyphyly of the order Rodentia, relationships among the Artiodactyla, and various interordinal relationships.