Range-wide and temporal genomic analyses reveal the consequences of near-extinction in Swedish moose

Ungulate species have experienced severe declines over the past centuries through overharvesting and habitat loss. Even if many game species have recovered thanks to strict hunting regulation, the genome-wide impacts of overharvesting are still unclear. Here, we examine the temporal and geographical differences in genome-wide diversity in moose (Alces alces) over its whole range in Sweden by sequencing 87 modern and historical genomes. We found limited impact of the 1900s near-extinction event but local variation in inbreeding and load in modern populations, as well as suggestion of a risk of future reduction in genetic diversity and gene flow. Furthermore, we found candidate genes for local adaptation, and rapid temporal allele frequency shifts involving coding genes since the 1980s, possibly due to selective harvesting. Our results highlight that genomic changes potentially impacting fitness can occur over short time scales and underline the need to track both deleterious and selectively advantageous genomic variation.

Tradeoffs between resources and risks shape the responses of a large carnivore to human disturbance

Wide-ranging carnivores experience tradeoffs between dynamic resource availabilities and heterogeneous risks from humans, with consequences for their ecological function and conservation outcomes. Yet, research investigating these tradeoffs across large carnivore distributions is rare. We assessed how resource availability and anthropogenic risks influence the strength of lion (Panthera leo) responses to disturbance using data from 31 sites across lions' contemporary range. Lions avoided human disturbance at over two-thirds of sites, though their responses varied depending on site-level characteristics. Lions were more likely to exploit human-dominated landscapes where resources were limited, indicating that resource limitation can outweigh anthropogenic risks and might exacerbate human-carnivore conflict. Lions also avoided human impacts by increasing their nocturnal activity more often at sites with higher production of cattle. The combined effects of expanding human impacts and environmental change threaten to simultaneously downgrade the ecological function of carnivores and intensify human-carnivore conflicts, escalating extinction risks for many species.

Mountain- and brown hare genetic polymorphisms to survey local adaptations and conservation status of the heath hare (Lepus timidus sylvaticus, Nilsson 1831)

We provide the first whole genome sequences from three specimens of the mountain hare subspecies the heath hare (Lepus timidus sylvaticus), along with samples from two mountain hares (Lepus timidus timidus) and two brown hares (Lepus europaeus) from Sweden. The heath hare has a unique grey winter pelage as compared to other mountain hares (white) and brown hares (mostly brown), and face regional extinction, likely due to competitive exclusion from the non-native brown hare. Whole genome resequencing from the seven hare specimens were mapped to the Lepus timidus pseudoreference genome and used for detection of 11,363,883 polymorphic nucleotide positions. The data presented here could be useful for addressing local adaptations and conservation status of mountain hares and brown hares in Sweden, including unique subspecies.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina HiSeq X
Factor Type(s)	Species
Sample Characteristic - Organism	Lepus timidus • Lepus europaeus • Lepus timidus sylvaticus
Sample Characteristic - Location	Sweden

A century of genetic homogenization in Baltic salmon-evidence from archival DNA

Intra-species genetic homogenization arising from anthropogenic impacts is a major threat to biodiversity. However, few taxa have sufficient historical material to systematically quantify long-term genetic changes. Using archival DNA collected over approximately 100 years, we assessed spatio-temporal genetic change in Atlantic salmon populations across the Baltic Sea, an area heavily impacted by hydropower exploitation and associated with large-scale mitigation stocking. Analysis was carried out by screening 82 SNPs in 1680 individuals from 13 Swedish rivers. We found an overall decrease in genetic divergence and diminished isolation by distance among populations, strongly indicating genetic homogenization over the past century. We further observed an increase in genetic diversity within populations consistent with increased gene flow. The temporal genetic change was lower in larger wild populations than in smaller wild and hatchery-reared ones, indicating that larger populations have been able to support a high number of native spawners in relation to immigrants. Our results demonstrate that stocking practices of salmon in the Baltic Sea have led to the homogenization of populations over the last century, potentially compromising their ability to adapt to environmental change. Stocking of reared fish is common worldwide, and our study is a cautionary example of the potentially long-term negative effects of such activities.

African and Asian leopards are highly differentiated at the genomic level

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Moving far, staying close: red fox dispersal patterns revealed by SNP genotyping

The genetic structure of a population can provide important insights into animal movements at varying geographical scales. Individual and social behaviors, such as philopatry and dispersal, affect patterns of relatedness, age and sex structure, shaping the local genetic structure of populations. However, these fine scale patterns may not be detected within broader population genetic structure. Using SNP genotyping for pairwise relatedness estimates, we investigated the spatial and genetic structuring of 141 red foxes within south-central Sweden at two scales. First, we looked at broad scale population structuring among red foxes at the regional level. We then estimated pairwise relatedness values to evaluate the spatial and genetic structure of male, female and mixed sex pairs for patterns of philopatry and dispersal at a more localized scale. We found limited genetic differentiation at the regional scale. However, local investigations revealed patterns of female philopatry and male biased dispersal. There were significant differences in pairwise geographic distances between highly related same sex pairs with the average distance between related males, 37.8 km, being six times farther than that of related females, averaging 6.3 km. In summary, the low levels of genetic differentiation found in this study illustrates the mobility and dispersal ability of red foxes across scales. However, relatedness plays a strong role in the spatial organization of red foxes locally, ultimately contributing to male biased dispersal patterns.

De novo discovery of SNPs for genotyping endangered sun parakeets (Aratinga solstitialis) in Guyana

Parrots (Psittaciformes) are among the most endangered groups of birds today and remain threatened by habitat loss and exploitation for the live bird trade. Under such conditions, reliable and non-invasive monitoring techniques are crucial for successful conservation measures. In this study, we developed a panel of 86 high quality SNPs for genotyping endangered sun parakeets (Aratinga solstitialis) in Guyana, which form one of the last known breeding populations of this South American species in the wild. Genotyping was tested on different types of samples (blood, feathers, feces, beak and cloacal swabs). While blood performed best, feathers and feces also yielded reliable results and could thus be used as non-invasive sources of DNA for future population monitoring. Discriminant Analysis of Principal Components (DAPC) on genotypes revealed that Guyanese sun parakeets clustered separately from other psittacine species as well as conspecifics from a captive population. A priori known first-order kinships were also adequately detected by the SNP panel. Using a series of experimental contaminations, we found that contamination from other psittacine species and slight contamination (~ 10%) from conspecifics did not prevent successful genotyping and recognition of individuals. We show that instances of higher conspecific contamination (~ 50%) can be detected through an increased level of heterozygosity that falls outside the distribution of uncontaminated samples.

Large-scale spatial variation of chronic stress signals in moose

The physiological effects of short-term stress responses typically lead to increased individual survival as it prepares the body for fight or flight through catabolic reactions in the body. These physiological effects trade off against growth, immunocompetence, reproduction, and even long-term survival. Chronic stress may thus reduce individual and population performance, with direct implications for the management and conservation of wildlife populations. Yet, relatively little is known about how chronic stress levels vary across wild populations and factors contributing to increased chronic stress levels. One method to measure long-term stress in mammals is to quantify slowly incorporated stress hormone (cortisol) in hair, which most likely reflect a long-term average of the stress responses. In this study, we sampled 237 harvested moose Alces alces across Sweden to determine the relative effect of landscape variables and disturbances on moose hair cortisol levels. We used linear model combinations and Akaike’s Information Criterion (corrected for small sample sizes), and included variables related to human disturbance, ungulate competition, large carnivore density, and ambient temperature to estimate the covariates that best explained the variance in stress levels in moose. The most important variables explaining the variation in hair cortisol levels in moose were the long-term average temperature sum in the area moose lived and the distance to occupied wolf territory; higher hair cortisol levels were detected where temperatures were higher and closer to occupied wolf territories, respectively.

Droplet digital PCR assays for the quantification of brown trout (Salmo trutta) and Arctic char (Salvelinus alpinus) from environmental DNA collected in the water of mountain lakes

Classical methods for estimating the abundance of fish populations are often both expensive, time-consuming and destructive. Analyses of the environmental DNA (eDNA) present in water samples could alleviate such constraints. Here, we developed protocols to detect and quantify brown trout (Salmo trutta) and Arctic char (Salvelinus alpinus) populations by applying the droplet digital PCR (ddPCR) method to eDNA molecules extracted from water samples collected in 28 Swedish mountain lakes. Overall, contemporary fish CPUE (catch per unit effort) estimates from standardized survey gill nettings were not correlated to eDNA concentrations for either of the species. In addition, the measured environmental variables (e.g. dissolved organic carbon concentrations, temperature, and pH) appear to not influence water eDNA concentrations of the studied fish species. Detection probabilities via eDNA analysis showed moderate success (less than 70% for both species) while the presence of eDNA from Arctic char (in six lakes) and brown trout (in one lake) was also indicated in lakes where the species were not detected with the gillnetting method. Such findings highlight the limits of one or both methods to reliably detect fish species presence in natural systems. Additional analysis showed that the filtration of water samples through 1.2 μm glass fiber filters and 0.45 μm mixed cellulose ester filters was more efficient in recovering DNA than using 0.22 μm enclosed polyethersulfone filters, probably due to differential efficiencies of DNA extraction. Altogether, this work showed the potentials and limits of the approach for the detection and the quantification of fish abundance in natural systems while providing new insights in the application of the ddPCR method applied to environmental DNA.

Carnivores, competition and genetic connectivity in the Anthropocene

Current extinction rates are comparable to five prior mass extinctions in the earth’s history, and are strongly affected by human activities that have modified more than half of the earth’s terrestrial surface. Increasing human activity restricts animal movements and isolates formerly connected populations, a particular concern for the conservation of large carnivores, but no prior research has used high throughput sequencing in a standardized manner to examine genetic connectivity for multiple species of large carnivores and multiple ecosystems. Here, we used RAD SNP genotypes to test for differences in connectivity between multiple ecosystems for African wild dogs (Lycaon pictus) and lions (Panthera leo), and to test correlations between genetic distance, geographic distance and landscape resistance due to human activity. We found weaker connectivity, a stronger correlation between genetic distance and geographic distance, and a stronger correlation between genetic distance and landscape resistance for lions than for wild dogs, and propose a new hypothesis that adaptations to interspecific competition may help to explain differences in vulnerability to isolation by humans.

Widespread introgression of mountain hare genes into Fennoscandian brown hare populations

In Fennoscandia, mountain hare (Lepus timidus) and brown hare (Lepus europaeus) hybridize and produce fertile offspring, resulting in gene flow across the species barrier. Analyses of maternally inherited mitochondrial DNA (mtDNA) show that introgression occur frequently, but unavailability of appropriate nuclear DNA markers has made it difficult to evaluate the scale- and significance for the species. The extent of introgression has become important as the brown hare is continuously expanding its range northward, at the apparent expense of the mountain hare, raising concerns about possible competition. We report here, based on analysis of 6833 SNP markers, that the introgression is highly asymmetrical in the direction of gene flow from mountain hare to brown hare, and that the levels of nuclear gene introgression are independent of mtDNA introgression. While it is possible that brown hares obtain locally adapted alleles from the resident mountain hares, the low levels of mountain hare alleles among allopatric brown hares suggest that hybridization is driven by stochastic processes. Interspecific geneflow with the brown hare is unlikely to have major impacts on mountain hare in Fennoscandia, but direct competition may.

Estimating population size using single-nucleotide polymorphism-based pedigree data

Reliable population estimates are an important aspect of sustainable wildlife management and conservation but can be difficult to obtain for rare and elusive species. Here, we test a new census method based on pedigree reconstruction recently developed by Creel and Rosenblatt (2013). Using a panel of 96 single-nucleotide polymorphisms (SNPs), we genotyped fecal samples from two Swedish brown bear populations for pedigree reconstruction. Based on 433 genotypes from central Sweden (CS) and 265 from northern Sweden (NS), the population estimates (N = 630 for CS, N = 408 for NS) fell within the 95% CI of the official estimates. The precision and accuracy improved with increasing sampling intensity. Like genetic capture-mark-recapture methods, this method can be applied to data from a single sampling session. Pedigree reconstruction combined with noninvasive genetic sampling may thus augment population estimates, particularly for rare and elusive species for which sampling may be challenging.

Phantoms of the forest: legacy risk effects of a regionally extinct large carnivore

The increased abundance of large carnivores in Europe is a conservation success, but the impact on the behavior and population dynamics of prey species is generally unknown. In Europe, the recolonization of large carnivores often occurs in areas where humans have greatly modified the landscape through forestry or agriculture. Currently, we poorly understand the effects of recolonizing large carnivores on extant prey species in anthropogenic landscapes. Here, we investigated if ungulate prey species showed innate responses to the scent of a regionally exterminated but native large carnivore, and whether the responses were affected by human-induced habitat openness. We experimentally introduced brown bear Ursus arctos scent to artificial feeding sites and used camera traps to document the responses of three sympatric ungulate species. In addition to controls without scent, reindeer scent Rangifer tarandus was used as a noncarnivore, novel control scent. Fallow deer Dama dama strongly avoided areas with bear scent. In the presence of bear scent, all ungulate species generally used open sites more than closed sites, whereas the opposite was observed at sites with reindeer scent or without scent. The opening of forest habitat by human practices, such as forestry and agriculture, creates a larger gradient in habitat openness than available in relatively unaffected closed forest systems, which may create opportunities for prey to alter their habitat selection and reduce predation risk in human-modified systems that do not exist in more natural forest systems. Increased knowledge about antipredator responses in areas subjected to anthropogenic change is important because these responses may affect prey population dynamics, lower trophic levels, and attitudes toward large carnivores. These aspects may be of particular relevance in the light of the increasing wildlife populations across much of Europe.

Using eDNA to experimentally test ungulate browsing preferences

Large herbivores may affect ecosystem processes and states, but such effects can be difficult to quantify, especially within multispecies assemblages. To better understand such processes and improve our predictive ability of systems undergoing change, herbivore diets can be studied using controlled feeding trials (or cafeteria tests). With some wildlife, such as large herbivores, it is impractical to empirically verify these findings, because it requires visually observing animals in forested environments, which can disturb them from their natural behaviors. Yet, in field-based cafeteria trials it is nearly impossible to differentiate selection between herbivore species that forage on similar plants and make very similar bite marks. However, during browsing ungulates leave saliva residue which includes some buccal cells and DNA that can be extracted for species identification. Here we used a newly developed eDNA-based method (biteDNA) to test the browsing preferences of four sympatric ungulate species in the wild. Overall, food preferences varied between species, but all species strongly preferred deciduous over coniferous species. Our method allows the study of plant-animal interactions in multispecies assemblages at very fine detail.

Single nucleotide polymorphism-based dispersal estimates using noninvasive sampling

Quantifying dispersal within wild populations is an important but challenging task. Here we present a method to estimate contemporary, individual-based dispersal distance from noninvasively collected samples using a specialized panel of 96 SNPs (single nucleotide polymorphisms). One main issue in conducting dispersal studies is the requirement for a high sampling resolution at a geographic scale appropriate for capturing the majority of dispersal events. In this study, fecal samples of brown bear (Ursus arctos) were collected by volunteer citizens, resulting in a high sampling resolution spanning over 45,000 km² in Gävleborg and Dalarna counties in Sweden. SNP genotypes were obtained for unique individuals sampled (n = 433) and subsequently used to reconstruct pedigrees. A Mantel test for isolation by distance suggests that the sampling scale was appropriate for females but not for males, which are known to disperse long distances. Euclidean distance was estimated between mother and offspring pairs identified through the reconstructed pedigrees. The mean dispersal distance was 12.9 km (SE 3.2) and 33.8 km (SE 6.8) for females and males, respectively. These results were significantly different (Wilcoxon’s rank-sum test: P-value = 0.02) and are in agreement with the previously identified pattern of male-biased dispersal. Our results illustrate the potential of using a combination of noninvasively collected samples at high resolution and specialized SNPs for pedigree-based dispersal models.

DNA left on browsed twigs uncovers bite-scale resource use patterns in European ungulates

Fine-scale resource use by large herbivores is often difficult to quantify directly. This is particularly true for browsing ungulates due to the challenges in observing shy subjects in forested environments of low visibility. As a consequence we know relatively little about resource use by diverse browsing ungulates. When browsing, ungulates leave behind saliva on the browsed twig that includes their DNA, which can be used to identify the species that was responsible for browsing the twig. We used this method, which we term "biteDNA", to study bite-scale browsing patterns in a temperate ungulate community. This approach provides a level of detail in browsing patterns across species that was previously very hard to attain. We found that all deer species largely overlapped in terms of the tree species they used. Moose browsed larger diameters than red deer and roe deer, but these latter two species did not differ. Moose browsed at higher heights than red deer, and red deer higher than roe deer. Although the deer species differed in mean browsing height, species were comparable in terms of their minimum browsing height of ~20 cm. This means that height and diameter ranges of the smaller species were found to be completely inside the ranges of the larger species. Hence, while moose may access exclusive food resources in terms of browse height and diameter, red and roe deer cannot.

Browsed twig environmental DNA: diagnostic PCR to identify ungulate species

Ungulate browsing can have a strong effect on ecological processes by affecting plant community structure and composition, with cascading effects on nutrient cycling and animal communities. However, in the absence of direct observations of foraging, species-specific foraging behaviours are difficult to quantify. We therefore know relatively little about foraging competition and species-specific browsing patterns in systems with several browsers. However, during browsing, a small amount of saliva containing buccal cells is deposited at the bite site, providing a source of environmental DNA (eDNA) that can be used for species identification. Here, we describe extraction and PCR protocols for a browser species diagnostic kit. Species-specific primers for mitochondrial DNA were optimized and validated using twigs browsed by captive animals. A time series showed that about 50% of the samples will amplify up to 12 weeks after the browsing event and that some samples amplify up to 24 weeks after browsing (12.5%). Applied to samples of natural browsing from an area where moose (Alces alces), roe deer (Capreolus capreolus), fallow deer (Cervus dama) and red deer (Cervus elaphus) are sympatric, amplification success reached 75%. This method promises to greatly improve our understanding of multispecies browsing systems without the need for direct observations.

Helpers increase the reproductive potential of offspring in cooperative meerkats

In both animal and human societies, individuals may forego personal reproduction and provide care to the offspring of others. Studies aimed at investigating the adaptive nature of such cooperative breeding systems in vertebrates typically calculate helper 'fitness' from relationships of helper numbers and offspring survival to independence. The aim of this study is to use observations and supplemental feeding experiments in cooperatively breeding meerkats, Suricata suricatta, to investigate whether helpers influence the long-term reproductive potential of offspring during adulthood. We show that helpers have a significant and positive influence on the probability that offspring gain direct reproductive success in their lifetimes. This effect arises because helpers both reduce the age at which offspring begin to reproduce as subordinates and increase the probability that they will compete successfully for alpha rank. Supplemental feeding experiments confirm the causality of these results. Our results suggest that one can neither discount the significance of helper effects when none is found nor necessarily estimate accurately the fitness benefit that helpers accrue, unless their effects on offspring are considered in the long term.

Intrasexual competition and sexual selection in cooperative mammals

In most animals, the sex that invests least in its offspring competes more intensely for access to the opposite sex and shows greater development of secondary sexual characters than the sex that invests most. However, in some mammals where females are the primary care-givers, females compete more frequently or intensely with each other than males. A possible explanation is that, in these species, the resources necessary for successful female reproduction are heavily concentrated and intrasexual competition for breeding opportunities is more intense among females than among males. Intrasexual competition between females is likely to be particularly intense in cooperative breeders where a single female monopolizes reproduction in each group. Here, we use data from a twelve-year study of wild meerkats (Suricata suricatta), where females show high levels of reproductive skew, to show that females gain greater benefits from acquiring dominant status than males and traits that increase competitive ability exert a stronger influence on their breeding success. Females that acquire dominant status also develop a suite of morphological, physiological and behavioural characteristics that help them to control other group members. Our results show that sex differences in parental investment are not the only mechanism capable of generating sex differences in reproductive competition and emphasize the extent to which competition for breeding opportunities between females can affect the evolution of sex differences and the operation of sexual selection.

Deriving dispersal distances from genetic data

Dispersal is one of the most important factors determining the genetic structure of a population, but good data on dispersal distances are rare because it is difficult to observe a large sample of dispersal events. However, genetic data contain unbiased information about the average dispersal distances in species with a strong sex bias in their dispersal rates. By plotting the genetic similarity between members of the philopatric sex against some measure of the distance between them, the resulting regression line can be used for estimating how far dispersing individuals of the opposite sex have moved before settling. Dispersers showing low genetic similarity to members of the opposite sex will on average have originated from further away. Applying this method to a microsatellite dataset from lions (Panthera leo) shows that their average dispersal distance is 1.3 home ranges with a 95% confidence interval of 0.4-3.0 home ranges. These results are consistent with direct observations of dispersal from our study population and others. In this case, direct observations of dispersal distance were not detectably biased by a failure to detect long-range dispersal, which is thought to be a common problem in the estimation of dispersal distance.

High genetic variation in leopards indicates large and long-term stable effective population size

In this paper we employ recently developed statistical and molecular tools to analyse the population history of the Tanzanian leopard (Panthera pardus), a large solitary felid. Because of their solitary lifestyle little is known of their past or present population dynamics. Eighty-one individuals were scored at 18 microsatellite loci. Overall, levels of heterozygosity were high (0.77 +/- 0.03), with a small heterozygote deficiency (0.06 +/- 0.03). Effective population size (Ne) was calculated to be 38 000-48 000. A Ne:N ratio of 0.42 (average from four cat studies) gives a present population size of about 100 000 leopards in Tanzania. Four different bottleneck tests indicated that this population has been large and stable for a minimum of several thousand years. FST values were low and no significant genetic structuring of the population could be detected. This concurs well with the large migration values (Nm) obtained (>3.3 individuals/generation). Our analysis reveals that ecological factors (e.g. disease), which are known to have had major impact on other carnivore populations, are unlikely to have impacted strongly on the population dynamics of Tanzanian leopards. The explanation may be found in their solitary life-style, their often nonconfrontational behaviour toward interspecific competitors, or that any bottlenecks have been of limited size, localized, or too short to have affected genetic variation to any measurable degree. Since the genetic structuring is weak, gene flow is not restricted to within protected areas. Local loss of genetic variation is therefore not of immediate concern.