Phylogeography and genetic variation in Western Jacob's ladder (Polemonium occidentale) provide insights into the origin and conservation of rare species in the Great Lakes region

The perennial herb Western Jacob's Ladder (Polemonium occidentale, Polemoniaceae) is widespread in the mountains of western North America but reappears as a disjunct in the Great Lakes Region in Minnesota and Wisconsin, USA as the narrow endemic P. occidentale subsp. lacustre. This distribution is shown by a diverse assemblage of angiosperms. It has been hypothesized that these species became isolated just after the Last Glacial Maximum, but this has not been tested. Additionally, the genetic diversity and population connectivity of the endemic Great Lakes flora has been understudied, with important conservation implications. Using genotyping-by-sequencing, we examined the relationship of P. occidentale subsp. lacustre to its closest relatives, relationships among all known populations, and genetic diversity within these populations. Polemonium occidentale subsp. lacustre represents an isolated, unique lineage that diverged from its closest relatives 1.3 Ma and arrived in the Great Lakes Region by at least 38 ka. Nearly all extant populations diverged prior to the Last Glacial Maximum, are genetically distinct, and show little within-population genetic diversity. Clonality may mitigate reduction in diversity due to drift. Mixed population signal between Wisconsin and some Minnesota populations may be due to gene flow during the Late Pleistocene. While populations of P. occidentale subsp. lacustre may be relictual from a now extinct western relative, it is best treated as a distinct species. Conservation efforts should focus more on ensuring that current populations remain rather than maintaining large populations sizes across a few populations. However, encouraging habitat heterogeneity may accomplish both simultaneously.

Spatial and temporal variability in the diet of Pacific marten (Martes caurina) on Haida Gwaii: an apex predator in a highly modified ecosystem

Knowledge of the diet ecology of apex predators in insular island ecosystems has direct applications to the conservation of endemic species at risk of extinction. We used stable isotopes of carbon and nitrogen to infer resource-use strategies of an indigenous predator, the Pacific marten (Martes caurina (Merriam, 1890)), in a highly modified ecosystem on Haida Gwaii, British Columbia, Canada. We used Bayesian isotopic mixing models to estimate the relative contributions of different food sources to marten diet across seasons and during overall lifetime, and to determine how diet varied with different levels of access to marine resources. Isotopes of carbon and nitrogen measured in hair and muscle tissue suggested that marten consumed salmon (3%–17%) and berries (25%–37%) seasonally; these diet groups were relatively minor components of the lifetime diet. Analysis of bone collagen suggested that terrestrial fauna — including birds, deer, small mammals, and invertebrates — contributed the most to diet (41%–55%), and marine invertebrates (38%–48%), not salmon (0%–3%), were the main allochthonous marine nutrient subsidy to lifetime diet. Plasticity in foraging ecology, combined with a broad dietary niche, introduced prey, notably the invasive Sitka black-tailed deer (Odocoileus hemionus sitkensis Merriam, 1898), as well as abundant marine resources, may allow marten to outcompete other indigenous and endemic carnivores on Haida Gwaii.

Appendicular skeletal morphology of North American Martes reflect independent modes of evolution in conjunction with Pleistocene glacial cycles

Pleistocene glacial cycles are thought to have driven ecological niche shifts, including novel niche formation. North American pine martens, Martes americana and M. caurina, are exemplar taxa thought to have diverged molecularly and morphologically during Pleistocene glaciation. Previous research found correlations between Martes limb morphology with biome and climate, suggesting that appendicular evolution may have occurred via adaptation to selective pressures imposed by novel and shifting habitats. Such variation can also be achieved through non-adaptive means such as genetic drift. Here, we evaluate whether regional genetic differences reflect limb morphology differences among populations of M. americana and M. caurina by analyzing evolutionary tempo and mode of six limb elements. Our comparative phylogenetic models indicate that genetic structure predicts limb shape better than size. Marten limb size has low phylogenetic signal, and the best supported model of evolution is punctuational (kappa), with morphological and genetic divergence occurring simultaneously. Disparity through time analysis suggests that the tempo of limb evolution in Martes tracks Pleistocene glacial cycles, such that limb size may be responding to shifting climates rather than population genetic structure. Contrarily, we find that limb shape is strongly tied to genetic relationships, with high phylogenetic signal and a lambda mode of evolution. Overall, this pattern of limb size and shape variation may be the result of geographic isolation during Pleistocene glacial advance, while declines in disparity suggest hybridization during interglacial periods. Future inclusion of extinct populations of Martes, which were more morphologically and ecologically diverse, may further clarify Martes phenotypic evolution.

Carnivore Contact: A Species Fracture Zone Delineated Amongst Genetically Structured North American Marten Populations (Martes americana and Martes caurina)

North American martens are forest dependent, influenced by human activity, and climate vulnerable. They have long been managed and harvested throughout their range as the American marten (Martes americana). Recent work has expanded evidence for the original description of two species in North America — M. americana and the Pacific Coast marten, M. caurina — but the geographic boundary between these groups has not been described in detail. From 2010 to 2016 we deployed 734 multi-taxa winter bait stations across a 53,474 km² study area spanning seven mountain ranges within the anticipated contact zone along the border of Canada and the United States. We collected marten hair samples and developed genotypes for 15 polymorphic microsatellite loci for 235 individuals, and 493 base-pair sequences of the mtDNA gene COI for 175 of those individuals. Both nuclear and mitochondrial genetic structure identified a sharp break across the Clark Fork Valley, United States with M. americana and M. caurina occurring north and south of the break, respectively. We estimated global effective population size (N_e) for each mountain range, clinal genetic neighborhood sizes (NS), calculated observed (H_o) and expected (H_e) heterozygosity, fixation index (F_ST), and clinal measures of allelic richness (Ar), H_o, and inbreeding coefficient (F_IS). Despite substantial genetic structure, we detected hybridization along the fracture zone with both contemporary (nuclear DNA) and historic (mtDNA) gene flow. Marten populations in our study area are highly structured and the break across the fracture zone being the largest documented in North America (F_ST range 0.21–0.34, mean = 0.27). With the exception of the Coeur d’Alene Mountains, marten were well distributed across higher elevation portions of our sampling area. Clinal NS values were variable suggesting substantial heterogeneity in marten density and movement. For both M. americana and M. caurina, elevationaly dependent gene flow and high genetic population structure suggest that connectivity corridors will be important to ensuring long-term population persistence. Our study is an example of how a combination of global and clinal molecular data analyses can provide important information for natural resource management.

Unraveling the mystery of the glacier bear: Genetic population structure of black bears (Ursus americanus) within the range of a rare pelage type

Glacier bears are a rare grey color morph of American black bear (Ursus americanus) found only in northern Southeast Alaska and a small portion of western Canada. We examine contemporary genetic population structure of black bears within the geographic extent of glacier bears and explore how this structure relates to pelage color and landscape features of a recently glaciated and highly fragmented landscape. We used existing radiocollar data to quantify black bear home-range size within the geographic range of glacier bears. The mean home-range size of female black bears in the study area was 13 km2 (n = 11), whereas the home range of a single male was 86.9 km2. We genotyped 284 bears using 21 microsatellites extracted from noninvasively collected hair as well as tissue samples from harvested bears. We found ten populations of black bears in the study area, including several new populations not previously identified, divided largely by geographic features such as glaciers and marine fjords. Glacier bears were assigned to four populations found on the north and east side of Lynn Canal and the north and west side of Glacier Bay with a curious absence in the nonglaciated peninsula between. Lack of genetic relatedness and geographic continuity between black bear populations containing glacier bears suggest a possible unsampled population or an association with ice fields. Further investigation is needed to determine the genetic basis and the adaptive and evolutionary significance of the glacier bear color morph to help focus black bear conservation management to maximize and preserve genetic diversity.

Ecology and Infection Dynamics of Multi-Host Amdoparvoviral and Protoparvoviral Carnivore Pathogens

Amdoparvovirus and Protoparvovirus are monophyletic viral genera that infect carnivores. We performed surveillance for and sequence analyses of parvoviruses in mustelids in insular British Columbia to investigate parvoviral maintenance and cross-species transmission among wildlife. Overall, 19.1% (49/256) of the tested animals were parvovirus-positive. Aleutian mink disease virus (AMDV) was more prevalent in mink (41.6%, 32/77) than martens (3.1%, 4/130), feline panleukopenia virus (FPV) was more prevalent in otters (27.3%, 6/22) than mink (5.2%, 4/77) or martens (2.3%, 3/130), and canine parvovirus 2 (CPV-2) was found in one mink, one otter, and zero ermines (N = 27). Viruses were endemic and bottleneck events, founder effects, and genetic drift generated regional lineages. We identified two local closely related AMDV lineages, one CPV-2 lineage, and five FPV lineages. Highly similar viruses were identified in different hosts, demonstrating cross-species transmission. The likelihood for cross-species transmission differed among viruses and some species likely represented dead-end spillover hosts. We suggest that there are principal maintenance hosts (otters for FPV, raccoons for CPV-2/FPV, mink for AMDV) that enable viral persistence and serve as sources for other susceptible species. In this multi-host system, viral and host factors affect viral persistence and distribution, shaping parvoviral ecology and evolution, with implications for insular carnivore conservation.