skip to main content


Title: Postglacial colonization in the Great Lakes Region by the white-footed mouse ( Peromyscus leucopus ): conflicts between genomic and field data
Abstract

For regions that were covered by ice during the Pleistocene glaciations, species must have emigrated from unglaciated regions. However, it can be difficult to discern when and from what ancestral source populations such expansions took place, especially since warming climates introduce the possibility of very recent expansions. For example, in the Great Lakes region, pronounced climatic change includes past glaciations as well as recent, rapid warming. Here we evaluate different expansion hypotheses with a genomic study of the white-footed mouse (Peromyscus leucopus noveboracensis), which is one of the most common mammals throughout the Great Lakes region. Ecological surveys coupled with historical museum records suggest a recent range expansion of P. leucopus associated with the warming climate over the last decades. These detailed records have yet to be complemented by genomic data that provide the requisite resolution for detecting recent expansion, although some mitochondrial DNA (mtDNA) sequences have suggested possible hypotheses about the geography of expansion. With more than 7,000 loci generated using RADseq, we evaluate support for multiple hypotheses of a geographic expansion in the Upper Peninsula of Michigan (UP). Analysis of a single random single-nucleotide polymorphism per locus revealed a fine-scale population structure separating the Lower Peninsula (LP) population from all other populations in the UP. We also detected a genetic structure that reflects an evolutionary history of postglacial colonization from two different origins into the UP, one coming from the LP and one coming from the west. Instead of supporting a climate-driven range expansion, as suggested by field surveys, our results support more ancient postglacial colonization of the UP from two different ancestral sources. With these results, we offer new insights about P. leucopus geographic expansion history, as well as a more general phylogeographic framework for testing range shifts in the Great Lakes region.

 
more » « less
NSF-PAR ID:
10392853
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Journal of Mammalogy
Volume:
103
Issue:
2
ISSN:
0022-2372
Page Range / eLocation ID:
p. 243-254
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weedSilene vulgarisacross Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in theLGMon the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion forS. vulgaris. We know thatS. vulgarishas recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity ofS. vulgarisin Europe. Therefore, the history of range expansion ofS. vulgarislikely began with postglacial expansion after theLGM, followed by more recent global dispersal by humans.

     
    more » « less
  2. Abstract Aim

    We test the predictions of the Stockholm Paradigm, a synthesis of eco‐evolutionary theory explaining the nature of faunal assembly, host range and parasite diversification. Faunal diversification and assembly, manifested in patterns of host colonization, co‐adaptation and parasite speciation, are predicted to emerge as a consequence of alternating episodes of ecological disruption and stability. Specifically, for a diverse cestode genus (Arostrilepis), we evaluate the number and direction of Pleistocene dispersal events across Beringia, the number and relative timing of host colonization events and the relationship between host and parasite biogeographic histories and associations through time.

    Location

    Beringia and adjacent temperate to arctic biomes in North America and Eurasia.

    Taxon

    Arostrilepis(Cyclophyllidea: Hymenolepididae) and its rodent hosts.

    Methods

    Multi‐locus phylogenetic reconstruction and biogeographic ancestral range estimation.

    Results

    Arostrilepislineages crossed Beringia eastward into North America a minimum of four times and westward into Asia twice in association with temporally disjunct geographic expansions of three major tribes of cricetid rodents (Arvicolini, Myodini, Lemmini). Inferences of ancestral host associations support at least nine instances of host colonization involving shifts from one rodent tribe or family to another. Several previously unrecognized lineages ofArostrilepisare revealed.

    Main conclusions

    Consistent with expectations of the Stockholm Paradigm, episodes of intercontinental dispersal were both frequent in the history ofArostrilepisand preceded a majority of inferred host‐colonization events. Events of historical geographic expansion created numerous opportunities for development of novel host–parasite associations through ecological fitting, as parasites tracked historically conserved resources available across diverse host taxa. Beringia played a major role in shaping rodent/parasite assemblages by mediating dispersal between the northern continents during glacial episodes of the Pleistocene, rather than by serving as a zone of refugial isolation.

     
    more » « less
  3. Synopsis

    Populations at the leading front of a range expansion must rapidly adapt to novel conditions. Increased epigenetic diversity has been hypothesized to facilitate adaptation and population persistence via non-genetic phenotypic variation, especially if there is reduced genetic diversity when populations expand (i.e., epigenetic diversity compensates for low genetic diversity). In this study, we use the spatial distribution of genetic and epigenetic diversity to test this hypothesis in populations of the white-footed mouse (Peromyscus leucopus) sampled across a purported recent range expansion gradient. We found mixed support for the epigenetic compensation hypothesis and a lack of support for expectations for expansion populations of mice at the range edge, which likely reflects a complex history of expansion in white-footed mice in the Upper Peninsula of Michigan. Specifically, epigenetic diversity was not increased in the population at the purported edge of the range expansion in comparison to the other expansion populations. However, input from an additional ancestral source populations may have increased genetic diversity at this range edge population, counteracting the expected genetic consequences of expansion, as well as reducing the benefit of increased epigenetic diversity at the range edge. Future work will expand the focal populations to include expansion areas with a single founding lineage to test for the robustness of a general trend that supports the hypothesized compensation of reduced genetic diversity by epigenetic variation observed in the expansion population that was founded from a single historical source.

     
    more » « less
  4. The phenology of critical biological events in aquatic ecosystems are rapidly shifting due to climate change. Growing variability in phenological cues can increase the likelihood of trophic mismatches, causing recruitment failures in commercially, culturally, and recreationally important fisheries. We tested for changes in spawning phenology of regionally important walleye (Sander vitreus) populations in 194 Midwest US lakes in Minnesota, Michigan, and Wisconsin spanning 1939-2019 to investigate factors influencing walleye phenological responses to climate change and associated climate variability, including ice-off timing, lake physical characteristics, and population stocking history. Data from Wisconsin and Michigan lakes (185 and 5 out of 194 total lakes, respectively) were collected by the Wisconsin Department of Natural Resources (WDNR) and the Great Lakes Indian Fish and Wildlife Commission (GLIFWC) through standardized spring walleye mark-recapture surveys and spring tribal harvest season records. Standardized spring mark-recapture population estimates are performed shortly after ice-off, where following a marking event, a subsequent recapture sampling event is conducted using nighttime electrofishing (typically AC – WDNR, pulsed-DC – GLIFWC) of the entire shoreline including islands for small lakes and index stations for large lakes (Hansen et al. 2015) that is timed to coincide with peak walleye spawning activity (G. Hatzenbeler, WDNR, personal communication; M. Luehring, GLIFWC, personal communication; Beard et al. 1997). Data for four additional Minnesota lakes were collected by the Minnesota Department of Natural Resources (MNDNR) beginning in 1939 during annual collections of walleye eggs and broodstock (Schneider et al. 2010), where date of peak egg take was used to index peak spawning activity. For lakes where spawning location did not match the lake for which the ice-off data was collected, the spawning location either flowed into (Pike River) or was within 50 km of a lake where ice-off data were available (Pine River) and these ice-off data were used. Following the affirmation of off-reservation Ojibwe tribal fishing rights in the Ceded Territories of Wisconsin and the Upper Peninsula of Michigan in 1987, tribal spearfishers have targeted walleye during spring spawning (Mrnak et al. 2018). Nightly harvests are recorded as part of a compulsory creel survey (US Department of the Interior 1991). Using these records, we calculated the date of peak spawning activity in a given lake-year as the day of maximum tribal harvest. Although we were unable to account for varying effort in these data, a preliminary analysis comparing spawning dates estimated using tribal harvest to those determined from standardized agency surveys in the same lake and year showed that they were highly correlated (Pearson’s correlation: r = 0.91, P < 0.001). For lakes that had walleye spawning data from both agency surveys and tribal harvest, we used the data source with the greatest number of observation years. Ice-off phenology data was collected from two sources – either observed from the Global Lake and River Ice Phenology database (Benson et al. 2000)t, or modeled from a USGS region-wide machine-learning model which used North American Land Data Assimilation System (NLDAS) meteorological inputs combined with lake characteristics (lake position, clarity, size, depth, hypsography, etc.) to predict daily water column temperatures from 1979 - 2022, from which ice-off dates could be derived (https://www.sciencebase.gov/catalog/item/6206d3c2d34ec05caca53071; see Corson-Dosch et al. 2023 for details). Modeled data for our study lakes (see (Read et al. 2021) for modeling details), which performed well in reflecting ice phenology when compared to observed data (i.e., highly significant correlation between observed and modeled ice-off dates when both were available; r = 0.71, p < 0.001). Lake surface area (ha), latitude, and maximum depth (m) were acquired from agency databases and lake reports. Lake class was based on a WDNR lakes classification system (Rypel et al. 2019) that categorized lakes based on temperature, water clarity, depth, and fish community. Walleye stocking history was defined using the walleye stocking classification system developed by the Wisconsin Technical Working Group (see also Sass et al. 2021), which categorized lakes based on relative contributions of naturally-produced and stocked fish to adult recruitment by relying heavily on historic records of age-0 and age-1 catch rates and stocking histories. Wisconsin lakes were divided into three groups: natural recruitment (NR), a combination of stocking and natural recruitment (C-ST), and stocked only (ST). Walleye natural recruitment was indexed as age-0 walleye CPE (number of age-0 walleye captured per km of shoreline electrofished) from WDNR and GLIFWC fall electrofishing surveys (see Hansen et al. 2015 for details). We excluded lake-years where stocking of age-0 fish occurred before age-0 surveys to only include measurements of naturally-reproduced fish. 
    more » « less
  5. The Amundsen Sea sector of Antarctica has long been considered the most vulnerable part of the West Antarctic Ice Sheet (WAIS) because of the great water depth at the grounding line, a subglacial bed seafloor deepening toward the interior of the continent, and the absence of substantial ice shelves. Glaciers in this configuration are thought to be susceptible to rapid or runaway retreat. Ice flowing into the Amundsen Sea Embayment is undergoing the most rapid changes of any sector of the Antarctic ice sheets outside the Antarctic Peninsula, including substantial grounding-line retreat over recent decades, as observed from satellite data. Recent models suggest that a threshold leading to the collapse of WAIS in this sector may have been already crossed and that much of the ice sheet could be lost even under relatively moderate greenhouse gas emission scenarios. Drill cores from the Amundsen Sea provide tests of several key questions about controls on ice sheet stability. The cores offer a direct offshore record of glacial history in a sector that is exclusively influenced by ice draining the WAIS, which allows clear comparisons between the WAIS history and low-latitude climate records. Today, relatively warm (modified) Circumpolar Deep Water (CDW) is impinging onto the Amundsen Sea shelf and causing melting under ice shelves and at the grounding line of the WAIS in most places. Reconstructions of past CDW intrusions can assess the ties between warm water upwelling and large-scale changes in past grounding-line positions. Carrying out these reconstructions offshore from the drainage basin that currently has the most substantial negative mass balance of ice anywhere in Antarctica is thus of prime interest to future predictions. The scientific objectives for this expedition are built on hypotheses about WAIS dynamics and related paleoenvironmental and paleoclimatic conditions. The main objectives are: 1. To test the hypothesis that WAIS collapses occurred during the Neogene and Quaternary and, if so, when and under which environmental conditions; 2. To obtain ice-proximal records of ice sheet dynamics in the Amundsen Sea that correlate with global records of ice-volume changes and proxy records for atmospheric and ocean temperatures; 3. To study the stability of a marine-based WAIS margin and how warm deepwater incursions control its position on the shelf; 4. To find evidence for the earliest major grounded WAIS advances onto the middle and outer shelf; 5. To test the hypothesis that the first major WAIS growth was related to the uplift of the Marie Byrd Land dome. International Ocean Discovery Program (IODP) Expedition 379 completed two very successful drill sites on the continental rise of the Amundsen Sea. Site U1532 is located on a large sediment drift, now called the Resolution Drift, and it penetrated to 794 m with 90% recovery. We collected almost-continuous cores from recent age through the Pleistocene and Pliocene and into the upper Miocene. At Site U1533, we drilled 383 m (70% recovery) into the more condensed sequence at the lower flank of the same sediment drift. The cores of both sites contain unique records that will enable study of the cyclicity of ice sheet advance and retreat processes as well as ocean-bottom water circulation and water mass changes. In particular, Site U1532 revealed a sequence of Pliocene sediments with an excellent paleomagnetic record for high-resolution climate change studies of the previously sparsely sampled Pacific sector of the West Antarctic margin. Despite the drilling success at these sites, the overall expedition experienced three unexpected difficulties that affected many of the scientific objectives: 1. The extensive sea ice on the continental shelf prevented us from drilling any of the proposed shelf sites. 2. The drill sites on the continental rise were in the path of numerous icebergs of various sizes that frequently forced us to pause drilling or leave the hole entirely as they approached the ship. The overall downtime caused by approaching icebergs was 50% of our time spent on site. 3. A medical evacuation cut the expedition short by 1 week. Recovery of core on the continental rise at Sites U1532 and U1533 cannot be used to indicate the extent of grounded ice on the shelf or, thus, of its retreat directly. However, the sediments contained in these cores offer a range of clues about past WAIS extent and retreat. At Sites U1532 and U1533, coarse-grained sediments interpreted to be ice-rafted debris (IRD) were identified throughout all recovered time periods. A dominant feature of the cores is recorded by lithofacies cyclicity, which is interpreted to represent relatively warmer periods variably characterized by sediments with higher microfossil abundance, greater bioturbation, and higher IRD concentrations alternating with colder periods characterized by dominantly gray laminated terrigenous muds. Initial comparison of these cycles to published late Quaternary records from the region suggests that the units interpreted to be records of warmer time intervals in the core tie to global interglacial periods and the units interpreted to be deposits of colder periods tie to global glacial periods. Cores from the two drill sites recovered sediments of dominantly terrigenous origin intercalated or mixed with pelagic or hemipelagic deposits. In particular, Site U1533, which is located near a deep-sea channel originating from the continental slope, contains graded silts, sands, and gravels transported downslope from the shelf to the rise. The channel is likely the pathway of these sediments transported by turbidity currents and other gravitational downslope processes. The association of lithologic facies at both sites predominantly reflects the interplay of downslope and contouritic sediment supply with occasional input of more pelagic sediment. Despite the lack of cores from the shelf, our records from the continental rise reveal the timing of glacial advances across the shelf and thus the existence of a continent-wide ice sheet in West Antarctica during longer time periods since at least the late Miocene. Cores from both sites contain abundant coarse-grained sediments and clasts of plutonic origin transported either by downslope processes or by ice rafting. If detailed provenance studies confirm our preliminary assessment that the origin of these samples is from the plutonic bedrock of Marie Byrd Land, their thermochronological record will potentially reveal timing and rates of denudation and erosion linked to crustal uplift. The chronostratigraphy of both sites enables the generation of a seismic sequence stratigraphy for the entire Amundsen Sea continental rise, spanning the area offshore from the Amundsen Sea Embayment westward along the Marie Byrd Land margin to the easternmost Ross Sea through a connecting network of seismic lines. 
    more » « less