Phylogeography of European moose (Alces alces) based on contemporary mtDNA data and archaeological records

Spatial ecology of moose in Sweden: Combined Sr-O-C isotope analyses of bone and antler

The study of spatial (paleo)ecology in mammals is critical to understand how animals adapt to and exploit their environment. In this work we analysed the 87Sr/86Sr, δ18O and δ13C isotope composition of 65 moose bone and antler samples from Sweden from wild-shot individuals dated between 1800 and 1994 to study moose mobility and feeding behaviour for (paleo)ecological applications. Sr data were compared with isoscapes of the Scandinavian region, built ad-hoc during this study, to understand how moose utilise the landscape in Northern Europe. The 87Sr/86Sr isoscape was developed using a machine-learning approach with external geo-environmental predictors and literature data. Similarly, a δ18O isoscape, obtained from average annual precipitation δ18O values, was employed to highlight differences in the isotope composition of the local environment vs. bone/antler. Overall, 82% of the moose samples were compatible with the likely local isotope composition (n = 53), suggesting that they were shot not far from their year-round dwelling area. 'Local' samples were used to calibrate the two isoscapes, to improve the prediction of provenance for the presumably 'non-local' individuals. For the latter (n = 12, of which two are antlers and ten are bones), the probability of geographic origin was estimated using a Bayesian approach by combining the two isoscapes. Interestingly, two of these samples (one antler and one bone) seem to come from areas more than 250 km away from the place where the animals were hunted, indicating a possible remarkable intra-annual mobility. Finally, the δ13C data were compared with the forest cover of Sweden and ultimately used to understand the dietary preference of moose. We interpreted a difference in δ13C values of antlers (13C-enriched) and bones (13C-depleted) as a joint effect of seasonal variations in moose diet and, possibly, physiological stresses during winter-time, i.e., increased consumption of endogenous 13C-depleted lipids.

Opportunities and challenges for monitoring a recolonizing large herbivore using citizen science

Monitoring is a prerequisite for evidence-based wildlife management and conservation planning, yet conventional monitoring approaches are often ineffective for species occurring at low densities. However, some species such as large mammals are often observed by lay people and this information can be leveraged through citizen science monitoring schemes. To ensure that such wildlife monitoring efforts provide robust inferences, assessing the quantity, quality, and potential biases of citizen science data is crucial. For Eurasian moose (Alces alces), a species currently recolonizing north-eastern Germany and occurring in very low numbers, we applied three citizen science tools: a mail/email report system, a smartphone application, and a webpage. Among these monitoring tools, the mail/email report system yielded the greatest number of moose reports in absolute and in standardized (corrected for time effort) terms. The reported moose were predominantly identified as single, adult, male individuals, and reports occurred mostly during late summer. Overlaying citizen science data with independently generated habitat suitability and connectivity maps showed that members of the public detected moose in suitable habitats but not necessarily in movement corridors. Also, moose detections were often recorded near roads, suggestive of spatial bias in the sampling effort. Our results suggest that citizen science-based data collection can be facilitated by brief, intuitive digital reporting systems. However, inference from the resulting data can be limited due to unquantified and possibly biased sampling effort. To overcome these challenges, we offer specific recommendations such as more structured monitoring efforts involving the public in areas likely to be roamed by moose for improving quantity, quality, and analysis of citizen science-based data for making robust inferences.

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

Clear phylogeographic pattern and genetic structure of wild boar Sus scrofa population in Central and Eastern Europe

The wild boar Sus scrofa is one of the widely spread ungulate species in Europe, yet the origin and genetic structure of the population inhabiting Central and Eastern Europe are not well recognized. We analysed 101 newly obtained sequences of complete mtDNA genomes and 548 D-loop sequences of the species and combined them with previously published data. We identified five phylogenetic clades in Europe with clear phylogeographic pattern. Two of them occurred mainly in western and central part of the continent, while the range of the third clade covered North-Eastern, Central and South-Eastern Europe. The two other clades had rather restricted distribution. In Central Europe, we identified a contact zone of three mtDNA clades. Population genetic structure reflected clear phylogeographic pattern of wild boar in this part of Europe. The contribution of lineages originating from the southern (Dinaric-Balkan) and eastern (northern cost of the Black Sea) areas to the observed phylogeographic pattern of the species in Central and Eastern Europe was larger than those from the regions located in southern France, Iberian, and Italian Peninsulas. The present work was the first mitogenomic analysis conducted in Central and Eastern Europe to study genetic diversity and structure of wild boar population.

The declining occurrence of moose (Alces alces) at the southernmost edge of its range raise conservation concerns

The border region between Austria, the Czech Republic, and Germany harbors the most south-western occurrence of moose in continental Europe. The population originated in Poland, where moose survived, immigrated from former Soviet Union or were reintroduced after the Second World War expanded west- and southwards. In recent years, the distribution of the nonetheless small Central European population seems to have declined, necessitating an evaluation of its current status. In this study, existing datasets of moose observations from 1958 to 2019 collected in the three countries were combined to create a database totaling 771 records (observations and deaths). The database was then used to analyze the following: (a) changes in moose distribution, (b) the most important mortality factors, and (c) the availability of suitable habitat as determined using a maximum entropy approach. The results showed a progressive increase in the number of moose observations after 1958, with peaks in the 1990s and around 2010, followed by a relatively steep drop after 2013. Mortality within the moose population was mostly due to human interactions, including 13 deadly wildlife-vehicle collisions, particularly on minor roads, and four animals that were either legally culled or poached. Our habitat model suggested that higher altitudes (ca. 700-1,000 m a.s.l.), especially those offering wetlands, broad-leaved forests and natural grasslands, are the preferred habitats of moose whereas steep slopes and areas of human activity are avoided. The habitat model also revealed the availability of large core areas of suitable habitat beyond the current distribution, suggesting that habitat was not the limiting factor explaining the moose distribution in the study area. Our findings call for immediate transboundary conservation measures to sustain the moose population, such as those aimed at preventing wildlife-vehicle collisions and illegal killings. Infrastructure planning and development activities must take into account the habitat requirements of moose.

Moose genomes reveal past glacial demography and the origin of modern lineages

Background

Numerous megafauna species from northern latitudes went extinct during the Pleistocene/Holocene transition as a result of climate-induced habitat changes. However, several ungulate species managed to successfully track their habitats during this period to eventually flourish and recolonise the holarctic regions. So far, the genomic impacts of these climate fluctuations on ungulates from high latitudes have been little explored. Here, we assemble a de-novo genome for the European moose (Alces alces) and analyse it together with re-sequenced nuclear genomes and ancient and modern mitogenomes from across the moose range in Eurasia and North America.

Results

We found that moose demographic history was greatly influenced by glacial cycles, with demographic responses to the Pleistocene/Holocene transition similar to other temperate ungulates. Our results further support that modern moose lineages trace their origin back to populations that inhabited distinct glacial refugia during the Last Glacial Maximum (LGM). Finally, we found that present day moose in Europe and North America show low to moderate inbreeding levels resulting from post-glacial bottlenecks and founder effects, but no evidence for recent inbreeding resulting from human-induced population declines.

Conclusions

Taken together, our results highlight the dynamic recent evolutionary history of the moose and provide an important resource for further genomic studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07208-3.

Size of the lower carnassial in the arctic and the red fox from Late Pleistocene in Belgium compared to other ancient and extant populations

Lengths, widths, and size proportions (length to width) of the lower carnassial were measured in 45 teeth of the arctic fox and 35 teeth of the red fox from Belgium radiocarbon dated to 46 640–14 120 ka BP. Data the Late Pleistocene foxes from Belgium were compared to 20 ancient and extant populations form Europe, Asia, and North America. The Pleistocene arctic fox from Belgium showed larger carnassial than in all recent samples of this species, whereas the Belgian fossil red foxes were characterized by the carnassial size comparable to that of the recent Siberian red foxes. Both fox species from the Pleistocene of Belgium showed the highest index of the carnassials length to width, which means increase in carnivorous adaptation. We conclude that the higher level of carnivorous specialization reached by the Belgian arctic and red foxes at the end of the Late Pleistocene reflected their scavenging on kills of large carnivores and human hunters (remains of megafauna). Harsh environmental conditions of that period and specific composition of ecosystems led to adapting to a more carnivorous food niche in both foxes.

Phylogeography of moose in western North America

Subspecies designations within temperate species’ ranges often reflect populations that were isolated by past continental glaciation, and glacial vicariance is believed to be a primary mechanism behind the diversification of several subspecies of North American cervids. We used genetics and the fossil record to study the phylogeography of three moose subspecies (Alces alces andersoni, A. a. gigas, and A. a. shirasi) in western North America. We sequenced the complete mitochondrial genome (16,341 base pairs; n = 60 moose) and genotyped 13 nuclear microsatellites (n = 253) to evaluate genetic variation among moose samples. We also reviewed the fossil record for detections of all North American cervids to comparatively assess the evidence for the existence of a southern refugial population of moose corresponding to A. a. shirasi during the last glacial maximum of the Pleistocene. Analysis of mtDNA molecular variance did not support distinct clades of moose corresponding to currently recognized subspecies, and mitogenomic haplotype phylogenies did not consistently distinguish individuals according to subspecies groupings. Analysis of population structure using microsatellite loci showed support for two to five clusters of moose, including the consistent distinction of a southern group of moose within the range of A. a. shirasi. We hypothesize that these microsatellite results reflect recent, not deep, divergence and may be confounded by a significant effect of geographic distance on gene flow across the region. Review of the fossil record showed no evidence of moose south of the Wisconsin ice age glaciers ≥ 15,000 years ago. We encourage the integration of our results with complementary analyses of phenotype data, such as morphometrics, originally used to delineate moose subspecies, for further evaluation of subspecies designations for North American moose.

Winter temperature correlates with mtDNA genetic structure of yellow-necked mouse population in NE Poland

We analysed a fragment (247 bp) of cytochrome b of mitochondrial DNA sequenced using 353 samples of yellow-necked mice Apodemus flavicollis trapped in seven forests and along three woodlot transects in north-eastern Poland. Our aims were to identify the phylogeographic pattern and mtDNA structure of the population and to evaluate the role of environmental conditions in shaping the spatial pattern of mtDNA diversity. We found out that three European haplogroups occurred sympatrically in north-eastern Poland. Inferences based on mtDNA haplotype distribution and frequency defined five subpopulations. The mtDNA-based structure of mice significantly correlated with winter temperature: frequency of Haplogroup 1 was positively, and that of Haplogroup 3 negatively correlated to mean temperature of January in the year of trapping. Synthesis of the published pan-European data on the species phylogeography also showed that the possibly 'thermophilous' Haplogroup 1 has the westernmost occurrence, whereas the more 'cold-resistant' Haplogroup 3 occurs much further to north-east than the other haplogroups. The observed patter may be a byproduct of the tight coevolution with nuclear genes, as we have earlier found that - in mice population in NE Poland - the spatial pattern of nuclear DNA was best explained by January temperature. Alternatively, the observed association of mitochondrial genetic variation with temperature is possible to be adaptive as cytochrome b is involved in the process of ATP production via oxidative phosphorylation.

Genetic analyses of brown hare (Lepus europaeus) support limited migration and translocation of Greek populations

Numerous studies have shown that the phylogeography of many species, including European brown hare, has been affected by the climatic oscillations of the Pleistocene. During this period the Balkans acted as a major refugium offering habitable conditions for many species. However, few studies have focused on the specific role of the Greek peninsula in the phylogeographic history of species in this southernmost margin of Balkans. We, therefore analyzed a 528 bp fragment of the D-loop region of mitochondrial DNA in 154 wild brown hare individuals from unsampled areas from both mainland and island Greece and compared it to 310 available brown hare sequences (including 110 Greek samples). Newly identified haplotypes show characteristic distribution in specific Greek areas reinforcing the theory that Greece can be considered as a subrefuge within Balkans for a number of species, with several “refugia within refugia” spots, holding significant genetic diversity. No haplotypes from wild Greek individuals clustered with the Central and Western Europe group revealing a minimal contribution of this area to the colonization of central Europe. One hundred and ten reared brown hares were also analyzed to elucidate the impact of introductions on local populations. Most of these samples presented close genetic affinity with haplotypes from Central and Western Europe indicating that farms in Greece use breeders from those areas. Therefore, despite human translocation of individuals, the genetic structure of brown hare has mostly been influenced by paleoclimatic conditions and minimally by human actions.

Mitochondrial DNA of pre-last glacial maximum red deer from NW Spain suggests a more complex phylogeographical history for the species

The major climatic oscillations that characterized the Quaternary had a great influence on the evolution and distribution of several species. During cold periods, the distribution of temperate‐adapted species became fragmented with many surviving in southern refugia (Iberian, Italian, and Balkan Peninsulas). Red deer was one of the species that contracted its original range to southern refugia. Currently, two main lineages have been described for the species: western and eastern. We have analyzed fossils pre‐dating the last glacial maximum (LGM) from Liñares cave (NW Spain) that belongs to the peripheral range of the western clade, and fossils from the Danish Holocene belonging to the central part of the same clade. Phylogenetic analyses place our samples in the western clade. However, some specimens from Liñares represent an early split in the tree along with other pre‐LGM western samples from previous studies. Despite low bootstrap values in the Bayesian phylogenies, haplotype networks connect these foreign haplotypes to the eastern clade. We suggest a mixed phylogeographical model to explain this pattern with range expansions from the east during the expansion phase after the cold periods in marine isotope stage 3. We find slight isolation by distance in post‐LGM populations that could be a consequence of the recolonization from southern refugia after the LGM.