Effects of hunting on mating, relatedness, and genetic diversity in a puma population.

Hunting mortality can affect population abundance, demography, patterns of dispersal and philopatry, breeding, and genetic diversity. We investigated the effects of hunting on the reproduction and genetic diversity in a puma population in western Colorado, USA. We genotyped over 11,000 single nucleotide polymorphisms (SNPs), using double-digest, restriction site-associated DNA sequencing (ddRADseq) in 291 tissue samples collected as part of a study on the effects of hunting on puma population abundance and demography in Colorado from 2004 to 2014. The study was designed with a reference period (years 1-5), during which hunting was suspended, followed by a treatment period (years 6-10), in which hunting was reinstated. Our objectives were to examine the effects of hunting on: (1) paternity and male reproductive success; (2) the relatedness between pumas within the population, and (3) genetic diversity. We found that hunting reduced the average age of male breeders. The number of unique fathers siring litters increased each year without hunting and decreased each year during the hunting period. Mated pairs were generally unrelated during both time periods, and females were more closely related than males. Hunting was also associated with increased relatedness among males and decreased relatedness among females in the population. Finally, genetic diversity increased during the period without hunting and decreased each year when hunting was present. This study demonstrates the utility of merging demographic data with large-scale genomic datasets in order to better understand the consequences of management actions. Specifically, we believe that this study highlights the need for long-term experimental research in which hunting mortality is manipulated, including at least one non-harvested control population, as part of a broader adaptive, zone management scheme.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species.

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250-km2 island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA), similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3000 single nucleotide polymorphisms in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing. Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Hunting alters viral transmission and evolution in a large carnivore.

Hunting can fundamentally alter wildlife population dynamics but the consequences of hunting on pathogen transmission and evolution remain poorly understood. Here, we present a study that leverages a unique landscape-scale quasi-experiment coupled with pathogen-transmission tracing, network simulation and phylodynamics to provide insights into how hunting shapes feline immunodeficiency virus (FIV) dynamics in puma (Puma concolor). We show that removing hunting pressure enhances the role of males in transmission, increases the viral population growth rate and increases the role of evolutionary forces on the pathogen compared to when hunting was reinstated. Changes in transmission observed with the removal of hunting could be linked to short-term social changes while the male puma population increased. These findings are supported through comparison with a region with stable hunting management over the same time period. This study shows that routine wildlife management can have impacts on pathogen transmission and evolution not previously considered.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi).

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

MrIML: Multi-response interpretable machine learning to model genomic landscapes.

We introduce a new R package "MrIML" ("Mister iml"; Multi-response Interpretable Machine Learning). MrIML provides a powerful and interpretable framework that enables users to harness recent advances in machine learning to quantify multilocus genomic relationships, to identify loci of interest for future landscape genetics studies, and to gain new insights into adaptation across environmental gradients. Relationships between genetic variation and environment are often nonlinear and interactive; these characteristics have been challenging to address using traditional landscape genetic approaches. Our package helps capture this complexity and offers functions that fit and interpret a wide range of highly flexible models that are routinely used for single-locus landscape genetics studies but are rarely extended to estimate response functions for multiple loci. To demonstrate the package's broad functionality, we test its ability to recover landscape relationships from simulated genomic data. We also apply the package to two empirical case studies. In the first, we model genetic variation of North American balsam poplar (Populus balsamifera, Salicaceae) populations across environmental gradients. In the second case study, we recover the landscape and host drivers of feline immunodeficiency virus genetic variation in bobcats (Lynx rufus). The ability to model thousands of loci collectively and compare models from linear regression to extreme gradient boosting, within the same analytical framework, has the potential to be transformative. The MrIML framework is also extendable and not limited to modelling genetic variation; for example, it can quantify the environmental drivers of microbiomes and coinfection dynamics.

Host relatedness and landscape connectivity shape pathogen spread in the puma, a large secretive carnivore.

Urban expansion can fundamentally alter wildlife movement and gene flow, but how urbanization alters pathogen spread is poorly understood. Here, we combine high resolution host and viral genomic data with landscape variables to examine the context of viral spread in puma (Puma concolor) from two contrasting regions: one bounded by the wildland urban interface (WUI) and one unbounded with minimal anthropogenic development (UB). We found landscape variables and host gene flow explained significant amounts of variation of feline immunodeficiency virus (FIV) spread in the WUI, but not in the unbounded region. The most important predictors of viral spread also differed; host spatial proximity, host relatedness, and mountain ranges played a role in FIV spread in the WUI, whereas roads might have facilitated viral spread in the unbounded region. Our research demonstrates how anthropogenic landscapes can alter pathogen spread, providing a more nuanced understanding of host-pathogen relationships to inform disease ecology in free-ranging species.

Urbanization impacts apex predator gene flow but not genetic diversity across an urban-rural divide.

Apex predators are important indicators of intact natural ecosystems. They are also sensitive to urbanization because they require broad home ranges and extensive contiguous habitat to support their prey base. Pumas (Puma concolor) can persist near human developed areas, but urbanization may be detrimental to their movement ecology, population structure, and genetic diversity. To investigate potential effects of urbanization in population connectivity of pumas, we performed a landscape genomics study of 130 pumas on the rural Western Slope and more urbanized Front Range of Colorado, USA. Over 12,000 single nucleotide polymorphisms (SNPs) were genotyped using double-digest, restriction site-associated DNA sequencing (ddRADseq). We investigated patterns of gene flow and genetic diversity, and tested for correlations between key landscape variables and genetic distance to assess the effects of urbanization and other landscape factors on gene flow. Levels of genetic diversity were similar for the Western Slope and Front Range, but effective population sizes were smaller, genetic distances were higher, and there was more admixture in the more urbanized Front Range. Forest cover was strongly positively associated with puma gene flow on the Western Slope, while impervious surfaces restricted gene flow and more open, natural habitats enhanced gene flow on the Front Range. Landscape genomic analyses revealed differences in puma movement and gene flow patterns in rural versus urban settings. Our results highlight the utility of dense, genome-scale markers to document subtle impacts of urbanization on a wide-ranging carnivore living near a large urban center.

Urbanization reduces genetic connectivity in bobcats (Lynx rufus) at both intra- and interpopulation spatial scales.

Urbanization is a major factor driving habitat fragmentation and connectivity loss in wildlife. However, the impacts of urbanization on connectivity can vary among species and even populations due to differences in local landscape characteristics, and our ability to detect these relationships may depend on the spatial scale at which they are measured. Bobcats (Lynx rufus) are relatively sensitive to urbanization and the status of bobcat populations is an important indicator of connectivity in urban coastal southern California. We genotyped 271 bobcats at 13,520 SNP loci to conduct a replicated landscape resistance analysis in five genetically distinct populations. We tested urban and natural factors potentially influencing individual connectivity in each population separately, as well as study-wide. Overall, landscape genomic effects were most frequently detected at the study-wide spatial scale, with urban land cover (measured as impervious surface) having negative effects and topographic roughness having positive effects on gene flow. The negative effect of urban land cover on connectivity was also evident when populations were analyzed separately despite varying substantially in spatial area and the proportion of urban development, confirming a pervasive impact of urbanization largely independent of spatial scale. The effect of urban development was strongest in one population where stream habitat had been lost to development, suggesting that riparian corridors may help mitigate reduced connectivity in urbanizing areas. Our results demonstrate the importance of replicating landscape genetic analyses across populations and considering how landscape genetic effects may vary with spatial scale and local landscape structure.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise.

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

Spinal arthritis in cane toads across the Australian landscape.

Loss of fitness can be a consequence of selection for rapid dispersal ability in invasive species. Increased prevalence of spinal arthritis may occur in cane toad populations at the invasion front as a cost of increased invasiveness, but our knowledge of the ecological drivers of this condition is lacking. We aimed to determine the factors explaining the prevalence of spinal arthritis in populations across the Australian landscape. We studied populations across a gradient of invasion histories. We collected 2415 toads over five years and determined the presence and size of spondylosis for each individual. We examined the effect of host size, leg length and invasion history on the prevalence of spondylosis. Host size was a significant predictor of spondylosis across populations. Contrary to our expectation, the overall prevalence of spondylosis was not positively related to invasion history and did not correlate with toad relative leg length. Rather than invasion age, the latitude at which populations were sampled provided an alternate explanation for the prevalence of spondylosis in cane toad populations and suggested that the incidence of this condition did not increase as a physiological cost of invasion, but is instead related to physical variables, such as climate.

Urban landscapes can change virus gene flow and evolution in a fragmentation-sensitive carnivore.

Urban expansion has widespread impacts on wildlife species globally, including the transmission and emergence of infectious diseases. However, there is almost no information about how urban landscapes shape transmission dynamics in wildlife. Using an innovative phylodynamic approach combining host and pathogen molecular data with landscape characteristics and host traits, we untangle the complex factors that drive transmission networks of feline immunodeficiency virus (FIV) in bobcats (Lynx rufus). We found that the urban landscape played a significant role in shaping FIV transmission. Even though bobcats were often trapped within the urban matrix, FIV transmission events were more likely to occur in areas with more natural habitat elements. Urban fragmentation also resulted in lower rates of pathogen evolution, possibly owing to a narrower range of host genotypes in the fragmented area. Combined, our findings show that urban landscapes can have impacts on a pathogen and its evolution in a carnivore living in one of the most fragmented and urban systems in North America. The analytical approach used here can be broadly applied to other host-pathogen systems, including humans.

Complementary Network-Based Approaches for Exploring Genetic Structure and Functional Connectivity in Two Vulnerable, Endemic Ground Squirrels.

The persistence of small populations is influenced by genetic structure and functional connectivity. We used two network-based approaches to understand the persistence of the northern Idaho ground squirrel (Urocitellus brunneus) and the southern Idaho ground squirrel (U. endemicus), two congeners of conservation concern. These graph theoretic approaches are conventionally applied to social or transportation networks, but here are used to study population persistence and connectivity. Population graph analyses revealed that local extinction rapidly reduced connectivity for the southern species, while connectivity for the northern species could be maintained following local extinction. Results from gravity models complemented those of population graph analyses, and indicated that potential vegetation productivity and topography drove connectivity in the northern species. For the southern species, development (roads) and small-scale topography reduced connectivity, while greater potential vegetation productivity increased connectivity. Taken together, the results of the two network-based methods (population graph analyses and gravity models) suggest the need for increased conservation action for the southern species, and that management efforts have been effective at maintaining habitat quality throughout the current range of the northern species. To prevent further declines, we encourage the continuation of management efforts for the northern species, whereas conservation of the southern species requires active management and additional measures to curtail habitat fragmentation. Our combination of population graph analyses and gravity models can inform conservation strategies of other species exhibiting patchy distributions.

Mixed population genomics support for the central marginal hypothesis across the invasive range of the cane toad (Rhinella marina) in Australia.

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235-33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad-scale, core-to-edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion-front region, northeast precipitation-limited region and southeast cold temperature-limited region. Core-to-edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.

Rangewide landscape genetics of an endemic Pacific northwestern salamander.

A species' genetic structure often varies in response to ecological and landscape processes that differ throughout the species' geographic range, yet landscape genetics studies are rarely spatially replicated. The Cope's giant salamander (Dicamptodon copei) is a neotenic, dispersal-limited amphibian with a restricted geographic range in the Pacific northwestern USA. We investigated which landscape factors affect D. copei gene flow in three regions spanning the species' range, which vary in climate, landcover and degree of anthropogenic disturbance. Least cost paths and Circuitscape resistance analyses revealed that gene flow patterns vary across the species' range, with unique combinations of landscape variables affecting gene flow in different regions. Populations in the northern coastal portions of the range had relatively high gene flow, largely facilitated by stream and river networks. Near the southeastern edge of the species' range, gene flow was more restricted overall, with relatively less facilitation by streams and more limitation by heat load index and fragmented forest cover. These results suggested that the landscape is more difficult for individuals to disperse through at the southeastern edge of the species' range, with terrestrial habitat desiccation factors becoming more limiting to gene flow. We suggest that caution be used when attempting to extrapolate landscape genetic models and conservation measures from one portion of a species' range to another.