Jennifer A. Moore

Effects of the landscape on boreal toad gene flow: does the pattern-process relationship hold true across distinct landscapes at the northern range margin?

Understanding the impact of natural and anthropogenic landscape features on population connectivity is a major goal in evolutionary ecology and conservation. Discovery of dispersal barriers is important for predicting population responses to landscape and environmental changes, particularly for populations at geographic range margins. We used a landscape genetics approach to quantify the effects of landscape features on gene flow and connectivity of boreal toad (Bufo boreas) populations from two distinct landscapes in south‐east Alaska (Admiralty Island, ANM, and the Chilkat River Valley, CRV). We used two common methodologies for calculating resistance distances in landscape genetics studies (resistance based on least‐cost paths and circuit theory). We found a strong effect of saltwater on genetic distance of CRV populations, but no landscape effects were found for the ANM populations. Our discordant results show the importance of examining multiple landscapes that differ in the variability of their features, to maximize detectability of underlying processes and allow results to be broadly applicable across regions. Saltwater serves as a physiological barrier to boreal toad gene flow and affects populations on a small geographic scale, yet there appear to be few other barriers to toad dispersal in this intact northern region.

Securing the Demographic and Genetic Future of Tuatara through Assisted Colonization

Climate change poses a particular threat to species with fragmented distributions and little or no capacity to migrate. Assisted colonization, moving species into regions where they have not previously occurred, aims to establish populations where they are expected to survive as climatic envelopes shift. However, adaptation to the source environment may affect whether species successfully establish in new regions. Assisted colonization has spurred debate among conservation biologists and ecologists over whether the potential benefits to the threatened species outweigh the potential disruption to recipient communities. In our opinion, the debate has been distracted by controversial examples, rather than cases where assisted colonization may be a viable strategy. We present a strategic plan for the assisted migration of tuatara (Sphenodon punctatus), an endemic New Zealand reptile. The plan includes use of extant populations as reference points for comparisons with assisted-colonization populations with respect to demography, phenotypic plasticity, and phenology; optimization of genetic variation; research to fill knowledge gaps; consideration of host and recipient communities; and inclusion of stakeholders in the planning stage. When strategically planned and monitored, assisted colonization could meet conservation and research goals and ultimately result in the establishment of long-term sustainable populations capable of persisting during rapid changes in climate.

Characterization of enteropathogenic and Shiga toxin-producing Escherichia coli in cattle and deer in a shared agroecosystem.

Shiga toxin-producing Escherichia coli (STEC) is an important foodborne pathogen. Cattle are suggested to be an important reservoir for STEC; however, these pathogens have also been isolated from other livestock and wildlife. In this study we sought to investigate transmission of STEC, enterohemorrhagic E. coli (EHEC) and enteropathogenic E. coli (EPEC) between cattle and white-tailed deer in a shared agroecosystem. Cattle feces were collected from 100 animals in a Michigan dairy farm in July 2012, while 163 deer fecal samples were collected during two sampling periods (March and June). The locations of deer fecal pellets were recorded via geographic information system mapping and microsatellite multi-locus genotyping was used to link the fecal samples to individual deer at both time points. Following subculture to sorbitol MacConkey agar and STEC CHROMagar, the pathogens were characterized by serotyping, stx profiling, and PCR-based fingerprinting; multilocus sequence typing (MLST) was performed on a subset. STEC and EHEC were cultured from 12 to 16% of cattle, respectively, and EPEC was found in 36%. Deer were significantly less likely to have a pathogen in March vs. June where the frequency of STEC, EHEC, and EPEC was 1, 6, and 22%, respectively. PCR fingerprinting and MLST clustered the cattle- and deer-derived strains together in a phylogenetic tree. Two STEC strains recovered from both animal species shared MLST and fingerprinting profiles, thereby providing evidence of interspecies transmission and highlighting the importance of wildlife species in pathogen shedding dynamics and persistence in the environment and cattle herds.

Application of large-scale parentage analysis for investigating natal dispersal in highly vagile vertebrates: a case study of American black bears (Ursus americanus).

Understanding the factors that affect dispersal is a fundamental question in ecology and conservation biology, particularly as populations are faced with increasing anthropogenic impacts. Here we collected georeferenced genetic samples (nu200a=u200a2,540) from three generations of black bears (Ursus americanus) harvested in a large (47,739 km2), geographically isolated population and used parentage analysis to identify mother-offspring dyads (nu200a=u200a337). We quantified the effects of sex, age, habitat type and suitability, and local harvest density at the natal and settlement sites on the probability of natal dispersal, and on dispersal distances. Dispersal was male-biased (76% of males dispersed) but a small proportion (21%) of females also dispersed, and female dispersal distances (mean ± SE u200a=u200a 48.9±7.7 km) were comparable to male dispersal distances (59.0±3.2 km). Dispersal probabilities and dispersal distances were greatest for bears in areas with high habitat suitability and low harvest density. The inverse relationship between dispersal and harvest density in black bears suggests that 1) intensive harvest promotes restricted dispersal, or 2) high black bear population density decreases the propensity to disperse. Multigenerational genetic data collected over large landscape scales can be a powerful means of characterizing dispersal patterns and causal associations with demographic and landscape features in wild populations of elusive and wide-ranging species.

Detection of Ophidiomyces, the Causative Agent of Snake Fungal Disease, in the Eastern Massasauga (Sistrurus catenatus) in Michigan, USA, 2014

Abstract Snake fungal disease (SFD), caused by Ophidiomyces ophiodiicola, threatens free-ranging snake populations across the US. We assayed 112 swabs from 102 individual eastern massasaugas (Sistrurus catenatus) at three locations in Michigan in 2014 for Ophidiomyces using quantitative PCR (qPCR). We observed a 12.7% qPCR prevalence of skin lesions. Individuals at each site had lesions, and occurrence of skin lesions was not significantly different between sites. We detected Ophidiomyces DNA at each of the three sites in five individuals (4.9%). We found no difference in detection probabilities between sites; however, snakes with dermatitis had higher Ophidiomyces DNA detection probabilities (P=0.15±0.08 SE) than snakes without dermatitis (P=0.02±0.01 SE, P=0.026). The emergence of SFD mortalities has potentially serious consequences for the viability of the eastern massasauga in Michigan. Future work should track temporal patterns in vital rates and health parameters, link health data to body condition indices for individual snakes, and conduct a “hotspot” analysis to examine health on a landscape scale.

Social network analysis of mating patterns in American black bears (Ursus americanus)

Nonrandom mating can structure populations and has important implications for population‐level processes. Investigating how and why mating deviates from random is important for understanding evolutionary processes as well as informing conservation and management. Prior to the implementation of parentage analyses, understanding mating patterns in solitary, elusive species like bears was virtually impossible. Here, we capitalize on a long‐term genetic data set collected from black bears (Ursus americanus) (N = 2422) in the Northern Lower Peninsula (NLP) of Michigan, USA. We identified mated pairs using parentage analysis and applied logistic regression (selection) models that controlled for features of the social network, to quantify the effects of individual characteristics, and spatial and population demographic factors on mating dynamics. Logistic regression models revealed that black bear mating was associated with spatial proximity of mates, male age, the time a pair had coexisted, local population density and relatedness. Mated pairs were more likely to contain older males. On average, bears tended to mate with nearby individuals to whom they were related, which does not support the existence of kin recognition in black bears. Pairwise relatedness was especially high for mated pairs containing young males. Restricted dispersal and high male turnover from intensive harvest mortality of NLP black bears are probably the underlying factors associated with younger male bears mating more often with female relatives. Our findings illustrate how harvest has the potential to disrupt the social structure of game species, which warrants further attention for conservation and management.

A Hybrid Zone between Canada Geese (Branta canadensis) and Cackling Geese (B. Hutchinsii)

ABSTRACT. n We studied patterns of geographic variation in structural size and genetic characteristics of white-cheeked geese inhabiting coastal areas of Hudson Bay, Canada, from northern Manitoba to southern Nunavut to determine the degree of morphological and spatial overlap, if any, between Cackling Geese (Branta hutchinsii) and Canada Geese (B. canadensis) in this region. Most Canada Geese occurred in sub-Arctic habitats south of 59°N latitude, and most Cackling Geese occurred in Arctic habitats north of 60°N, but the two species overlapped in a narrow zone between 59°N and 60°N latitude that coincided with the ecotone between sub-Arctic and Arctic ecozones. Mismatches between morphological and genetic characteristics of some individual females suggested that introgression had occurred in this area, and contrasting patterns in the nuclear and mitochondrial DNA (mtDNA) were consistent with female natal philopatry and male-biased dispersal. Evidence of introgression in the nuclear genome was geographically more widespread than evidence of introgression in the mtDNA genome. We suggest that the persistence of Canada Goose mtDNA in phenotypic Cackling Geese is a result of historical hybridization events that occurred when the Arctic—sub-Arctic ecotone was located farther north during a warmer climatic period. Despite evidence of introgression, most birds that we sampled appeared to belong to one or the other parental species, on the basis of their consistent identification using morphological, mtDNA, and nuclear DNA characteristics. We suggest that the area of overlap represents a tension zone between Canada Geese and Cackling Geese that is maintained by behavioral and ecological factors that limit effective dispersal.

Climatic and geographic predictors of life history variation in Eastern Massasauga (Sistrurus catenatus): A range-wide synthesis

Elucidating how life history traits vary geographically is important to understanding variation in population dynamics. Because many aspects of ectotherm life history are climate-dependent, geographic variation in climate is expected to have a large impact on population dynamics through effects on annual survival, body size, growth rate, age at first reproduction, size–fecundity relationship, and reproductive frequency. The Eastern Massasauga (Sistrurus catenatus) is a small, imperiled North American rattlesnake with a distribution centered on the Great Lakes region, where lake effects strongly influence local conditions. To address Eastern Massasauga life history data gaps, we compiled data from 47 study sites representing 38 counties across the range. We used multimodel inference and general linear models with geographic coordinates and annual climate normals as explanatory variables to clarify patterns of variation in life history traits. We found strong evidence for geographic variation in six of nine life history variables. Adult female snout-vent length and neonate mass increased with increasing mean annual precipitation. Litter size decreased with increasing mean temperature, and the size–fecundity relationship and growth prior to first hibernation both increased with increasing latitude. The proportion of gravid females also increased with increasing latitude, but this relationship may be the result of geographically varying detection bias. Our results provide insights into ectotherm life history variation and fill critical data gaps, which will inform Eastern Massasauga conservation efforts by improving biological realism for models of population viability and climate change.

Intestinal Microbial Community Dynamics of White-Tailed Deer (Odocoileus virginianus) in an Agroecosystem

The intestinal microbiota has important functions that contribute to host health. The compositional dynamics of microbial communities are affected by many factors, including diet and presence of pathogens. In contrast to humans and domestic mammals, the composition and seasonal dynamics of intestinal microbiota of wildlife species remain comparatively understudied. White-tailed deer (Odocoileus virginianus) is an ecologically and economically important wildlife species that inhabits agricultural ecosystems and is known to be a reservoir of enteric pathogens. Nevertheless, there is a lack of knowledge of white-tailed deer intestinal microbiota diversity and taxonomic composition. This study’s first objective was to characterize and compare the intestinal microbiota of 66 fecal samples from white-tailed deer collected during two sampling periods (March and June) using 16S rDNA pyrosequencing. Associations between community diversity and composition and factors including season, sex, host genetic relatedness, and spatial location were quantified. Results revealed that white-tailed deer intestinal microbiota was predominantly comprised of phyla Firmicutes and Proteobacteria, whose relative frequencies varied significantly between sampling periods. The second objective was to examine the associations between the presence of Escherichia coli and Salmonella, and microbiota composition and diversity. Results indicated that relative abundance of some microbial taxa varied when a pathogen was present. This study provides insights into microbial compositional dynamics of a wildlife species inhabiting coupled natural and agricultural landscapes. Data focus attention on the high prevalence of Proteobacteria particularly during the summer and highlight the need for future research regarding the role of white-tailed deer as a natural pathogen reservoir in agroecosystems.

Detecting black bear source-sink dynamics using individual-based genetic graphs.

Source–sink dynamics affects population connectivity, spatial genetic structure and population viability for many species. We introduce a novel approach that uses individual-based genetic graphs to identify source–sink areas within a continuously distributed population of black bears (Ursus americanus) in the northern lower peninsula (NLP) of Michigan, USA. Black bear harvest samples (n = 569, from 2002, 2006 and 2010) were genotyped at 12 microsatellite loci and locations were compared across years to identify areas of consistent occupancy over time. We compared graph metrics estimated for a genetic model with metrics from 10 ecological models to identify ecological factors that were associated with sources and sinks. We identified 62 source nodes, 16 of which represent important source areas (net flux > 0.7) and 79 sink nodes. Source strength was significantly correlated with bear local harvest density (a proxy for bear density) and habitat suitability. Additionally, resampling simulations showed our approach is robust to potential sampling bias from uneven sample dispersion. Findings demonstrate black bears in the NLP exhibit asymmetric gene flow, and individual-based genetic graphs can characterize source–sink dynamics in continuously distributed species in the absence of discrete habitat patches. Our findings warrant consideration of undetected source–sink dynamics and their implications on harvest management of game species.