RESUMEN
The skin microbiome of amphibians can influence host susceptibility towards the fungal pathogen Batrachochytrium dendrobatidis (Bd), while simultaneously having the potential to be altered by Bd. Severe Bd infections are known to alter the amphibian skin microbiome; however, little is known about microbiome interactions in amphibians with low infection intensity. In addition to disease dynamics, environmental factors may influence the microbiome. To test for patterns in bacterial diversity based on pathogen infection and environmental factors, 399 Columbia spotted frogs (Rana luteiventris) were sampled throughout northern Idaho and northeastern Washington across two years. Bd prevalence and intensity were measured in 376 frogs, revealing a prevalence of 69%, but generally low infection intensity (Mean = 127 Bd zoospore equivalents among infected frogs). Skin bacterial communities were characterized in 92 frogs using 16S rRNA gene amplicon sequencing. Our results indicated correlations of decreasing Shannon diversity and evenness as infection intensity increased. Latitude was correlated with bacterial richness and Faith's Phylogenetic Diversity measures, indicating increased diversity in northern locations. Beta diversity (UniFrac) analyses revealed that skin microbiomes were distinct between infected and uninfected frogs, and infection intensity had a significant effect on microbiome composition. Site explained the majority of microbiome variation (weighted UniFrac: 57.5%), suggesting a combination of local habitat conditions explain variation, as only small proportions of variation could be explained by year, month, temperature, elevation, and latitude individually. Bacterial genera with potential for Bd-inhibitory properties were found with differential relative abundance in infected and uninfected frogs, with higher Stenotrophomonas and lower Pseudomonas relative abundance observed in infected frogs. Further study may indicate if Bd inhibition by members of the skin microbiome is an influence behind the low infection intensities observed and whether low Bd infection intensities are capable of altering skin microbiome composition.
RESUMEN
The Northern Rocky Mountain ecosystem supports rich biological diversity with many endemic and rare species. Extant endemics display two biogeographic patterns: widespread species with fragmented populations, and narrow-range endemics. These distributions are shown by the congeneric snails Anguispira kochi occidentalis and Anguispira nimapuna. These two taxa are disjunct from the remaining species of the genus, which achieves its greatest diversity in eastern North America. Given the disjunct nature of A. k. occidentalis and A. nimapuna, we here present a mtDNA phylogeny of the genus that includes both eastern and western species to assess the phylogenetic position of A. k. occidentalis and A. nimapuna. We then reconstruct the demographic history of A. k. occidentalis and A. nimapuna by analysing current patterns of genetic variation and interpreting the results considering the historical biogeography of the region. Both A. k. occidentalis and A. nimapuna represent unique taxa that are genetically and geographically distinct from their congeners. The current distribution and genetic structure of A. k. occidentalis has been shaped by both historical isolation in refugia and more recent northward shifts, whereas A. nimapuna is represented by two populations with shallow divergence in an area of long-term habitat stability.
RESUMEN
The application of species distribution models (SDMs) to areas outside of where a model was created allows informed decisions across large spatial scales, yet transferability remains a challenge in ecological modeling. We examined how regional variation in animal-environment relationships influenced model transferability for Canada lynx (Lynx canadensis), with an additional conservation aim of modeling lynx habitat across the northwestern United States. Simultaneously, we explored the effect of sample size from GPS data on SDM model performance and transferability. We used data from three geographically distinct Canada lynx populations in Washington (n = 17 individuals), Montana (n = 66), and Wyoming (n = 10) from 1996 to 2015. We assessed regional variation in lynx-environment relationships between these three populations using principal components analysis (PCA). We used ensemble modeling to develop SDMs for each population and all populations combined and assessed model prediction and transferability for each model scenario using withheld data and an extensive independent dataset (n = 650). Finally, we examined GPS data efficiency by testing models created with sample sizes of 5%-100% of the original datasets. PCA results indicated some differences in environmental characteristics between populations; models created from individual populations showed differential transferability based on the populations' similarity in PCA space. Despite population differences, a single model created from all populations performed as well, or better, than each individual population. Model performance was mostly insensitive to GPS sample size, with a plateau in predictive ability reached at ~30% of the total GPS dataset when initial sample size was large. Based on these results, we generated well-validated spatial predictions of Canada lynx distribution across a large portion of the species' southern range, with precipitation and temperature the primary environmental predictors in the model. We also demonstrated substantial redundancy in our large GPS dataset, with predictive performance insensitive to sample sizes above 30% of the original.
RESUMEN
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 km2 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.
RESUMEN
Winter bait stations are becoming a commonly used technique for multispecies inventory and monitoring but a technical evaluation of their effectiveness is lacking. Bait stations have three components: carcass attractant, remote camera, and hair snare. Our 22,975 km2 mountainous study area was stratified with a 5 × 5 km sampling grid centered on northern Idaho and including portions of Washington, Montana, and British Columbia. From 2010-14, we conducted 563 sampling sessions at 497 bait stations in 453 5 × 5 km cells. We evaluated the effectiveness of cameras and hair snare DNA collection at stations to detect species and individual animals, factors affecting DNA viability, the effectiveness of re-visiting stations, and the influence of elevation, seasonality, and latency on species detections. Cameras were more effective at detecting multiple species than DNA hair snaring. Length of deployment time and elevation increased genetic species ID success but individual ID success rates were increased only by collecting hairs earlier in the season. Re-visiting stations did not change camera or genetic species detection results but did increase the number of individual genotypes identified. Marten and fisher were detected quickly while bobcat and coyote showed longer latency to detection. Seasonality significantly affected coyote and bobcat detections but not marten, fisher, or weasel. Multispecies bait station study design should incorporate mixed elevation sites with stratified seasonality. Priority should be given to including cameras as components of bait stations over hair snares, unless there is a specific genetic goal to the study. A hair snare component should be added, however, if individual ID or genetic data are necessary. Winter stations should be deployed a minimum of 45-60 days to allow for detection of low density species and species with long latency to detection times. Hair samples should be collected prior to DNA-degrading late season rain events. Re-visiting stations does not change which species are detected at stations; therefore, studies with objectives to delineate species presence or distribution will be more effective if they focus on deploying more stations across a broader landscape in lieu of surveying the same site multiple times.
