Reconstructing range dynamics and range fragmentation of European bison for the last 8000 years

A reduced SNP panel optimised for non-invasive genetic assessment of a genetically impoverished conservation icon, the European bison

The European bison was saved from the brink of extinction due to considerable conservation efforts since the early twentieth century. The current global population of > 9500 individuals is the result of successful ex situ breeding based on a stock of only 12 founders, resulting in an extremely low level of genetic variability. Due to the low allelic diversity, traditional molecular tools, such as microsatellites, fail to provide sufficient resolution for accurate genetic assessments in European bison, let alone from non-invasive samples. Here, we present a SNP panel for accurate high-resolution genotyping of European bison, which is suitable for a wide variety of sample types. The panel accommodates 96 markers allowing for individual and parental assignment, sex determination, breeding line discrimination, and cross-species detection. Two applications were shown to be utilisable in further Bos species with potential conservation significance. The new SNP panel will allow to tackle crucial tasks in European bison conservation, including the genetic monitoring of reintroduced populations, and a molecular assessment of pedigree data documented in the world's first studbook of a threatened species.

Millennial processes of population decline, range contraction and near extinction of the European bison

European bison (Bison bonasus) were widespread throughout Europe during the late Pleistocene. However, the contributions of environmental change and humans to their near extinction have never been resolved. Using process-explicit models, fossils and ancient DNA, we disentangle the combinations of threatening processes that drove population declines and regional extinctions of European bison through space and across time. We show that the population size of European bison declined abruptly at the termination of the Pleistocene in response to rapid environmental change, hunting by humans and their interaction. Human activities prevented populations of European bison from rebounding in the Holocene, despite improved environmental conditions. Hunting caused range loss in the north and east of its distribution, while land use change was responsible for losses in the west and south. Advances in hunting technologies from 1500 CE were needed to simulate low abundances observed in 1870 CE. While our findings show that humans were an important driver of the extinction of the European bison in the wild, vast areas of its range vanished during the Pleistocene-Holocene transition because of post-glacial environmental change. These areas of its former range have been climatically unsuitable for millennia and should not be considered in reintroduction efforts.

Demographic Expansions and the Emergence of Host Specialization in Genetically Distinct Ecotypes of the Tick-Transmitted Bacterium Anaplasma phagocytophilum

In Europe, genetically distinct ecotypes of the tick-vectored bacterium Anaplasma phagocytophilum circulate among mammals in three discrete enzootic cycles. To date, potential ecological factors that contributed to the emergence of these divergent ecotypes have been poorly studied. Here, we show that the ecotype that predominantly infects roe deer (Capreolus capreolus) is evolutionarily derived. Its divergence from a host generalist ancestor occurred after the last glacial maximum as mammal populations, including roe deer, recolonized the European mainland from southern refugia. We also provide evidence that this host specialist ecotype's effective population size (N_e) has tracked changes in the population of its roe deer host. Specifically, both host and bacterium have undergone substantial increases in N_e over the past 1,500 years. In contrast, we show that while it appears to have undergone a major population expansion starting ~3,500 years ago, in the past 500 years, the contemporary host generalist ecotype has experienced a substantial reduction in genetic diversity levels, possibly as a result of reduced opportunities for transmission between competent hosts. IMPORTANCE The findings of this study reveal specific events important for the evolution of host specialization in a naturally occurring, obligately intracellular bacterial pathogen. Specifically, they show that host range shifts and the emergence of host specialization may occur during periods of population growth in a generalist ancestor. Our results also demonstrate the close correlation between demographic patterns in host and pathogen for a specialist system. These findings have important relevance for understanding the evolution of host range diversity. They may inform future work on host range dynamics, and they provide insights for understanding the emergence of pathogens that have human and veterinary health implications.

The Reintroduction Analysis of European Bison (Bison bonasus L., 1758) in the North of Romania and the Identification of the Most Favourable Locations

We analysed the possibility of reintroducing the European bison (Bison bonasus L.) in the north of Romania—in Suceava, Neamț, and Maramureș counties—as well as increasing the wild European bison population in Neamț county to improve the genetic quality of the existing population. Currently, there is a population of over 50 individuals in the wild in Vânători Neamț Natural Park, Romania. At the same time, an attempt was made to identify the connecting corridors between the free European bison in Neamț county and other populations through new nuclei of European bison released in Suceava and Maramureș counties. In this regard, the hunting grounds with the highest ecological potential for the analysed species were identified using GIS spatial analysis techniques. The aim was also to trace possible ecological corridors linking different reintroduction locations, taking into account the ecological claims of the species. The analysis also followed the size of the European bison groups to be released, the sex ratio of each group, and the periodicity of their releasing. In order to reach viable populations, scenarios and simulations were carried out depending on the age, number, and sex of the relocated specimens. In this regard, the dynamics and the minimum viable population that could survive without risk of extinction were highlighted. The analysis showed that the analysed area has a high potential for the reintroduction of European bison in the wild. The release and creation of new European bison nuclei in the wild creates the premises for natural contacts with the existing free populations in the wild, genetic improvement, and increasing fauna diversity with ecological, social, and economic implications.

Wild Steps in a semi-wild setting? Habitat selection and behavior of European bison reintroduced to an enclosure in an anthropogenic landscape

Recently, several wild or semi-wild herds of European bison have been reintroduced across Europe. It is essential for future successful bison reintroductions to know how the European bison use different habitats, which environmental parameters drive their habitat selection, and whether their habitat use and behavioural patterns in new reintroduction sites differ from habitats where European bison have been roaming freely for a long time. Here, we address these questions for a 40-ha enclosed site that has been inhabited by semi-free ranging European bison since 2012. The site, Vorup Meadows, is adjacent to the Gudenå river in Denmark and consists of human-modified riparian meadows. During 2013 we monitored the behavioural pattern and spatial use of the 11 bison present and in parallel carried out floristic analyses to assess habitat structure and food quality in the enclosure. We tested habitat use and selection against environmental parameters such as habitat characteristics, plant community traits, topography, and management area (release area vs. meadow area) using linear regression and spatial models. The bison herd had comparable diurnal activity patterns as observed in previous studies on free-roaming bison herds. Topography emerged as the main predictor of the frequency of occurrence in our spatial models, with high-lying drier areas being used more. Bison did not prefer open areas over areas with tree cover when accounting for habitat availability. However, they spent significantly more time in the release area, a former agricultural field with supplementary fodder, than expected from availability compared to the rest of the enclosure, a meadow with tree patches. We wish to increase awareness of possible long-term ethological effects of the release site and the management protocols accomplished here that might reduce the ecological impact by the bison in the target habitat, and thereby compromise or even oppose the conservation goals of the conservation efforts.

Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe’s largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

Millennial-scale faunal record reveals differential resilience of European large mammals to human impacts across the Holocene

The use of short-term indicators for understanding patterns and processes of biodiversity loss can mask longer-term faunal responses to human pressures. We use an extensive database of approximately 18 700 mammalian zooarchaeological records for the last 11 700 years across Europe to reconstruct spatio-temporal dynamics of Holocene range change for 15 large-bodied mammal species. European mammals experienced protracted, non-congruent range losses, with significant declines starting in some species approximately 3000 years ago and continuing to the present, and with the timing, duration and magnitude of declines varying individually between species. Some European mammals became globally extinct during the Holocene, whereas others experienced limited or no significant range change. These findings demonstrate the relatively early onset of prehistoric human impacts on postglacial biodiversity, and mirror species-specific patterns of mammalian extinction during the Late Pleistocene. Herbivores experienced significantly greater declines than carnivores, revealing an important historical extinction filter that informs our understanding of relative resilience and vulnerability to human pressures for different taxa. We highlight the importance of large-scale, long-term datasets for understanding complex protracted extinction processes, although the dynamic pattern of progressive faunal depletion of European mammal assemblages across the Holocene challenges easy identification of ‘static’ past baselines to inform current-day environmental management and restoration.

Deciphering the Wisent Demographic and Adaptive Histories from Individual Whole-Genome Sequences

As the largest European herbivore, the wisent (Bison bonasus) is emblematic of the continent wildlife but has unclear origins. Here, we infer its demographic and adaptive histories from two individual whole-genome sequences via a detailed comparative analysis with bovine genomes. We estimate that the wisent and bovine species diverged from 1.7 × 106 to 850,000 years before present (YBP) through a speciation process involving an extended period of limited gene flow. Our data further support the occurrence of more recent secondary contacts, posterior to the Bos taurus and Bos indicus divergence (∼150,000 YBP), between the wisent and (European) taurine cattle lineages. Although the wisent and bovine population sizes experienced a similar sharp decline since the Last Glacial Maximum, we find that the wisent demography remained more fluctuating during the Pleistocene. This is in agreement with a scenario in which wisents responded to successive glaciations by habitat fragmentation rather than southward and eastward migration as for the bovine ancestors. We finally detect 423 genes under positive selection between the wisent and bovine lineages, which shed a new light on the genome response to different living conditions (temperature, available food resource, and pathogen exposure) and on the key gene functions altered by the domestication process.

Land abandonment and changes in snow cover period accelerate range expansions of sika deer

Ongoing climate change and land‐use change have the potential to substantially alter the distribution of large herbivores. This may result in drastic changes in ecosystems by changing plant–herbivore interactions. Here, we developed a model explaining sika deer persistence and colonization between 25 years in terms of neighborhood occupancy and habitat suitability. We used climatic, land‐use, and topographic variables to calculate the habitat suitability and evaluated the contributions of the variables to past range changes of sika deer. We used this model to predict the changes in the range of sika deer over the next 100 years under four scenario groups with the combination of land‐use change and climate change. Our results showed that both climate change and land‐use change had affected the range of sika deer in the past 25 years. Habitat suitability increased in northern or mountainous regions, which account for 71.6% of Japan, in line with a decrease in the snow cover period. Habitat suitability decreased in suburban areas, which account for 28.4% of Japan, corresponding to land‐use changes related to urbanization. In the next 100 years, the decrease in snow cover period and the increase in land abandonment were predicted to accelerate the range expansion of sika deer. Comparison of these two driving factors revealed that climate change will contribute more to range expansion, particularly from the 2070s onward. In scenarios that assumed the influence of both climate change and land‐use change, the total sika deer range increased by between +4.6% and +11.9% from the baseline scenario. Climate change and land‐use change will require additional efforts for future management of sika deer, particularly in the long term.

Can Species Distribution Models Aid Bioassessment when Reference Sites are Lacking? Tests Based on Freshwater Fishes

Recent literature reviews of bioassessment methods raise questions about use of least-impacted reference sites to characterize natural conditions that no longer exist within contemporary landscapes. We explore an alternate approach for bioassessment that uses species site occupancy data from museum archives as input for species distribution models (SDMs) stacked to predict species assemblages of freshwater fishes in Texas. When data for estimating reference conditions are lacking, deviation between richness of contemporary versus modeled species assemblages could provide a means to infer relative biological integrity at appropriate spatial scales. We constructed SDMs for 100 freshwater fish species to compare predicted species assemblages to data on contemporary assemblages acquired by four independent surveys that sampled 269 sites. We then compared site-specific observed/predicted ratios of the number of species at sites to scores from a multimetric index of biotic integrity (IBI). Predicted numbers of species were moderately to strongly correlated with the numbers observed by the four surveys. We found significant, though weak, relationships between observed/predicted ratios and IBI scores. SDM-based assessments identified patterns of local assemblage change that were congruent with IBI inferences; however, modeling artifacts that likely contributed to over-prediction of species presence may restrict the stand-alone use of SDM-derived patterns for bioassessment and therefore warrant examination. Our results suggest that when extensive standardized survey data that include reference sites are lacking, as is commonly the case, SDMs derived from generally much more readily available species site occupancy data could be used to provide a complementary tool for bioassessment.

Next-generation dynamic global vegetation models: learning from community ecology

Dynamic global vegetation models (DGVMs) are powerful tools to project past, current and future vegetation patterns and associated biogeochemical cycles. However, most models are limited by how they define vegetation and by their simplistic representation of competition. We discuss how concepts from community assembly theory and coexistence theory can help to improve vegetation models. We further present a trait- and individual-based vegetation model (aDGVM2) that allows individual plants to adopt a unique combination of trait values. These traits define how individual plants grow and compete. A genetic optimization algorithm is used to simulate trait inheritance and reproductive isolation between individuals. These model properties allow the assembly of plant communities that are adapted to a site's biotic and abiotic conditions. The aDGVM2 simulates how environmental conditions influence the trait spectra of plant communities; that fire selects for traits that enhance fire protection and reduces trait diversity; and the emergence of life-history strategies that are suggestive of colonization-competition trade-offs. The aDGVM2 deals with functional diversity and competition fundamentally differently from current DGVMs. This approach may yield novel insights as to how vegetation may respond to climate change and we believe it could foster collaborations between functional plant biologists and vegetation modellers.