Spatial ecology of the stone marten in an Alpine area: combining camera-trapping and genetic surveys

A species’ potential distribution can be modelled adequately only if no factor other than habitat availability affects its occurrences. Space use by stone marten Martes foina is likely to be affected by interspecific competition with the strictly related pine marten Martes martes, the latter being able to outcompete the first species in forested habitats. Hence, to point out the environmental factors which determine the distribution and density of the stone marten, a relatively understudied mesocarnivore, we applied two non-invasive survey methods, camera-trapping and faecal-DNA based genetic analysis, in an Alpine area where the pine marten was deemed to be absent (Val Grande National Park N Italy). Camera trapping was conducted from October 2014 to November 2015, using up to 27 cameras. Marten scats were searched for between July and November 2015 and, to assess density, in spring 2017. Species identification was accomplished by a PCR-RFLP method, while 17 autosomal microsatellites were used for individual identification. The stone marten occurred in all available habitats (83% of trapping sites and 73.2% of scats); nonetheless, habitat suitability, as assessed using MaxEnt, depended on four major land cover variables—rocky grasslands, rocks and debris, beech forests and chestnut forests—, martens selecting forests and avoiding open rocky areas. Sixteen individuals were identified, of which 14 related to each other, possibly forming six different groups. Using capwire estimators, density was assessed as 0.95 (0.7–1.3) ind/km2. In the study area, the widespread stone marten selected forested areas, attaining density values like those reported for the pine marten in northern Europe and suggesting that patterns of habitat selection may depend on the relative abundance of the two competing martens.


Introduction
Affecting ecosystem function, structure and dynamics, mesocarnivores play important roles in natural communities and are considered sensitive indicators of environmental health and change in forested and aquatic habitats, particularly wherever large carnivores have been driven to extinction by human interference (Buskirk and Zielinski 2003;Roemer et al. 2009). Evaluating mesocarnivore distribution and abundance is thus essential for investigating trophic cascades and predator-prey density-dependent relationships and conservation-aimed management (Williams et al. 2002).
Nonetheless, with a few exceptions (red fox Vulpes vulpes, European badger Meles meles and Eurasian otter Lutra lutra), there is relatively little research on mesocarnivores (Brooke et al. 2014). Among the mustelid family, which largely contribute to the diversity of European mesocarnivores, the stone marten Martes foina has currently been understudied (Proulx et al. 2004), probably as an indirect consequence of its absence in the British Isles, which have long played a leader role in the investigation of mustelid ecology (e.g. Mcdonald 2002;O'Mahony et al. 2017;Mathews et al. 2018). Although the stone marten is widespread through much of continental Europe and central Asia-from Portugal in the west as far as north-western China in the east (Proulx et al. 2004)-, data on its spatial ecology are scarce (e.g. Abramov et al. 2006;Ruiz-Gonzalez et al. 2008), research having been focused on habitat selection (reviewed by Virgós et al. 2012). To the best of our knowledge, density data are reported by only three available studies, conducted in either rural (Switzerland 0.7-2 ind/km 2 , Lachat Feller 1993; Germany: 2 adult and 1.5 juvenile ind/ km 2 , Herrmann 2004) or urban areas (4.7-5.8 ind/km 2 , Herr et al. 2009) by means of radiotracking. In northern Italy, stone marten density has been assessed in an agricultural area using camera-trapping and the Random Encounter Model proposed by Rowcliffe et al. (2008). Although population density was similar to that recorded in rural Switzerland (0.96 ind/km 2 ; Ronchi 2016), the REM has been demonstrated to largely underestimate marten density (Balestrieri et al. 2016a).
Although being best adapted to warm climates (Proulx et al. 2004), the stone marten has been recorded from sea level up to 4200 m in Nepal (Oli 1994), while in Europe it occurs up to 2400 m a.s.l. on the Alps (Genovesi and De Marinis 2003). As other Martes species, the stone marten prefers forested habitats (Virgós et al. 2012); in southern Europe, it is often associated to mosaics of forest and field patches (Sacchi and Meriggi 1995;Werner 2012;Vergara et al. 2017); nonetheless, wherever available, forests are selected (Virgós and Casanovas 1998;Ruiz-Gonzalez et al. 2015;Zub et al. 2018).
Habitat use by the stone marten is considered to be driven by competition with the pine marten (Martes martes), which should be able to outcompete the stone marten in forested habitats (Delibes 1983). More recently, it has been suggested that the strictly nocturnal stone marten may be more tolerant of human disturbance than the cathemeral pine marten, and thus may manage better than the latter in rural habitats (Balestrieri et al. 2019). Nonetheless, in the Iberian Peninsula, southwards of the southern edge of pine marten distribution, the stone marten is a mainly forest-dwelling species, which selects low human density areas (Virgós and Casanovas 1998).
Currently, forests are more widespread in mountainous areas than in European lowlands: on the Alps, wood cover has progressively increased since the 1960s, following the abandonment of low-intensity farming and livestock rearing (Falcucci et al. 2007), with a positive effect on forest-dwelling species (MacDonald et al. 2000). Pine-and stone marten are sympatric throughout the Eastern Italian Alps, while, according to available data, the latter would be by far most widespread than the pine marten in the central part of the mountain range (Fonda 2019).
Aiming to assess stone marten distribution, habitat selection and density in a pine marten-free, forested habitat, we focused on the Val Grande National Park, a protected area in the Lepontine Alps, at the border between the western and central sector of the mountain range, where both available data for the whole Piedmont region (Sindaco and Carpegna 2010) and park rangers' records indicated that pine marten occurrence may be null or negligible (only two records were available: one in 1900 and one in 2006, both outside the Park). We applied two non-invasive methods, camera-trapping and faecal-DNA-based genetic sampling. Genotyping needed two steps, species identification, which was necessary to exclude the samples belonging to other sympatric mesocarnivores from further analyses, and microsatellite genotyping, to ascertain the minimum number of individuals occurring in the study area (Ruiz-Gonzalez et al. 2008, 2013.

Study area
The Val Grande National Park (Piedmont region, Verbano-Cusio-Ossola province; 46°01′ 45″ N 8°27′ 34″ E) is the largest wilderness area of the Alps (153.7 km 2 ). The abandonment of traditional land use practices since the end of World War II has led to the decline of cultivated lands (meadows, pastures, chestnut orchards, crops and vineyards) from 59% of the whole area at the end of the 19th century to 5% in 1999 (Höchtl et al. 2005). Most previously cultivated areas currently show various successional stages. Woods mainly consist of beech Fagus sylvatica and chestnut Castanea sativa and cover ca. 55% of the protected area. Mean yearly temperature and yearly rainfall are respectively 6.5°C and 2300 mm with wide variations depending on altitude above sea level, which ranges between 250 and 2300 m.

Materials and methods
Several studies have suggested that the simultaneous use of multiple survey methods may provide a more complete assessment of mammal diversity (Silveira et al. 2003;Li et al. 2012;Croose et al. 2019). Mustelids are elusive and their densities being typically low, large areas need to be sampled to assess their distribution and habitat preferences. Moreover, Martes spp. are very similar and their field signs cannot be distinguished by eye (Davison et al. 2002), making their monitoring even more challenging. To overcome these hindrances, we applied two non-invasive methods, which, based on the characteristics of both the study area and target species, were deemed to offer the best balance between cost-effectiveness and monitoring efficiency (Roberts 2011;Balestrieri et al. 2016aBalestrieri et al. , 2016b.

Camera-trapping
The study area was monitored by digital scouting cameratraps (Acorn II LTL 5210 with Passive Infra-Red motion sensor), tied to trees 30-50 cm above the ground level and set to record 30-s-long videoclips, with no interval between two successive recordings. Camera-trap sites were georeferenced and superimposed on digital maps. We used an active survey design, attracting animals into the detection zone of the camera trap by placing scent lures (cat food in carnivore-proof containers) in front of camera-traps (ca. 5 m away). The use of videoclips and lures aimed to improve the opportunity of observing the distinctive morphological traits (van Maanen 2013) of martens. In monochrome images, the three more conspicuous features are the shape and position of the ears, the paler colour of chest and thighs respect to forelimbs, hocks and tail in the stone marten and its chunkier overall silhouette.
All videos of Martes spp. were subjected to a blind identification procedure by three experienced researchers (AM, FF and AB) and discordant records were discarded. Capture independence was achieved by considering consecutive records of the same species at the same site within a 30-min interval as a single event (Kelly and Holub 2008;Monterroso et al. 2014).
From October 2014 to April 2015, the whole study area, between 450 and 1720 m a.s.l., was surveyed by deploying 21 camera-traps as regularly as possible, depending on accessibility (mean inter-trap distance ± SD = 2.6 ± 1.3 km; Fig. 1). In winter (Dec-Apr), only 10 trap-sites were activated, avoiding avalanche-prone sections. Between July and November 2015, we focused on the north-south corridor formed by the valleys "de il Fiume" and "Pogallo". Camera-traps (N = 27) were deployed within a 1 × 1 km grid superimposed on the kilometric grid of digitalized, 1:10,000 Regional Technical Maps, aiming to set each camera as much as possible in the centre of the grid mesh and sample the most representative habitat (mean inter-trap distance ± SD = 1.0 ± 0.3 km; Fig. 1).
Variation over sampling periods in each species' encounter rates was tested by the chi-squared (χ 2 ) test.

Faecal DNA-based species identification
To supplement the data on stone marten distribution collected by camera-trapping, non-invasive genetic sampling was conducted between July and November 2015. Fresh scats were searched for by two surveyors along linear transects coinciding with paths (N = 23; mean length = 7.2 km, min-max: 2.5-14.8 km). In the north-south corridor transects were surveyed monthly, while those in the rest of the protected area were sampled mostly once, and all samples were georeferenced by a GPS. A small portion (ca. 1 cm) of each scat suspected to belong to Martes spp. based on size and morphology (see Remonti et al. 2012) was picked up using sticks, stored in autoclaved tubes containing ethanol 96% and preserved at -20°C until processed.
DNA was isolated using the QIAamp DNA Stool Mini Kit (Qiagen) according to the manufacturer's instructions and species identification was accomplished by a PCR-RFLP method. Two primers amplify the mtDNA from Martes martes, M. foina and four Mustela species; then, the simultaneous digestion of amplified mtDNA by two restriction enzymes (RsaI and HaeIII) generates different restriction patterns for each mustelid species, providing for an effective genetic identification of sympatric marten species (Ruiz-Gonzalez et al. 2008, 2013.

Individual identification by microsatellite genotyping
To assess the minimum number of stone marten individuals, a specific sampling was carried out in the lower Pogallo valley, from April to June 2017, as to record only resident adult individuals (i.e. after juvenile dispersal and before kits-of-the year start marking; Libois and Waechter 1991). By superimposing a 1 × 1 km grid on 1:10,000 maps, we identified 24, 1-km 2 contiguous sub-areas with adequate paths and accessibility to the major habitats and altitude belts (500-1000, 1001-1500 and 1501-2000 m a.s.l.). Grid mesh size was based on available data on both stone-and pine marten density (0.1-3.16 ind/km 2 ; Marchesi 1989; Zalewski and Jedrzejewski 2006;  Balestrieri et al. 2016a), aiming to lower the chance of missing those individuals whose home ranges may fall in un-sampled areas (Kays and Slauson 2008). Each transect was surveyed twice (total length = 83.3 km), with a gap of 15-20 days between visits.
All samples were georeferenced and stored at -20°C in autoclaved tubes containing 96% ethanol until processed. Faecal DNA samples were genotyped using a multiplex panel of 17 autosomal microsatellite markers including 10 speciesspecific microsatellites (Mf 1.1, Mf 1.11, Mf 1.18, Mf 1.3, Mf 2.13, Mf 3.2, Mf 3.7, Mf 4.17, Mf 8.7, Mf 8.8; Basto et al. 2010) and 7 additional markers described in closely related mustelids (Ma1, Davis and Strobeck 1998;Mel1, Bijlsma et al. 2000;MLUT27, Cabria et al. 2007;Mvis072, Fleming et al. 1999;Mvi57, O'Connell et al. 1996;MP0059, Jordan et al. 2007; Lut453, Dallas et al. 2003), previously used in To lower the probability of retaining false homozygotes or false allele errors, a multitube-approach of 4 independent replicates was used (Taberlet et al. 1996), followed by stringent criteria to construct consensus genotypes (i.e. accepting heterozygotes if the two alleles are recorded in ≥ 2 replicates and homozygotes if a single allele is recorded in ≥ 3 replicates) (e.g. Frantz et al. 2003;Brzeski et al. 2013). Briefly, DNA quality was initially screened by PCR-amplifying each DNA sample four times at four loci (Mult-A) and only samples showing > 50% positive PCRs were further amplified four times at the remaining 12 loci. Samples with ambiguous results after four amplifications per locus or with < 50% successful amplifications across loci were not considered reliable genotypes and discarded (Ruiz-González et al. 2013). Multiplex PCR products were run on an ABI (Foster City, CA) 3130XL automated sequencer (Applied Biosystems), with the internal size standard GS500 LIZ™ (Applied Biosystems) and fragment analysis was conducted using the ABI software Genemapper 4.0. To test the discrimination power of our microsatellite set, we computed the probability of identity (PID) by GIMLET, using the unbiased equation for both small sample size and siblings. The more conservative PID for full-sibs (PID-Sib) was estimated as an upper limit to the probability that pairs of individuals would share the same genotype. Consensus genotypes from four replicates were reconstructed using GIMLET, which was also used to estimate genotyping errors: allelic dropout (ADO) and false alleles (FA) (Taberlet et al. 1996;Pompanon et al. 2005).

Estimation of population size and kinship from genetic data
Population size was assessed by capwire estimators (Miller et al. 2005), an urn model developed expressly for faecal DNA-based sampling which provides reliable estimates also for small populations and has been used to estimate population size in several species (e.g. Arrendal et al. 2007;Sugimoto et al. 2014). To obtain the maximum likelihood estimate (MLE) of population size, data were fitted to either the Equal Capture model, for which all individuals were assumed to have an equal probability of being sampled, or Two-Innate Rates model, assuming that the population contained a mixture of easy-to-capture and difficult-to-capture individuals. The fit of the two models was compared using a Likelihood Ratio Test (LRT) and the p-value was calculated by using a parametric bootstrap approach to estimate the distribution of the LRT for data simulated under the lessparameterized Equal Capture model (Pennell et al. 2013). Confidence intervals (CI) for population size were estimated using a parametric bootstrap approach (Miller et al. 2005).
Genetic relatedness and sibling analyses were calculated by ML-RELATE (Kalinowski et al. 2006) which uses a maximum likelihood method to compute pair-wise genetic relatedness (Rxy). Sibship analysis was conducted using COLONY 2.0.4 (Jones and Wang 2010), with the typing error rate set at 0.01. This approach considers the likelihood of the entire pedigree, as opposed to relatedness on a pair-wise basis.

Species distribution modelling
Land cover variables, extracted from available maps of forested areas of Piedmont region (Table 1; IPLA 2016), were resampled to a common resolution of 1 × 1 km cell size, using QuantumGIS 2.16.3. To test for multi-collinearity among variables, the Variance Inflation Factor (VIF) was calculated (Table 1), VIF values > 3 indicating highly correlated  (Phillips et al. 2006), a widely used method which applies the principle of maximum entropy to predict the potential distribution of species from presence-only data (Phillips and Dudík 2008;Elith et al. 2011), and has proved efficient for assessing habitat suitability for the stone marten (Vergara et al. 2016). Models were fitted to independent presence data, collected between October 2014 and June 2017, and an equal number of randomly selected background points (Drake 2014). To ensure more ecologically realistic response curves, MaxEnt was run using only linear and quadratic features (Bateman et al. 2012(Bateman et al. , 2016 and default values for all the other parameters (maximum number of iterations = 5000; convergence threshold = 10 -5 ; multiplier regularization = 1). To assess the relative contribution of each variable, a jackknife test was used (Phillips et al. 2006). Particularly, we estimated the arithmetic mean between the percent contribution and permutation importance, two measures which define the contribution of each variable to the final model Meyer et al. 2014).
To obtain the best model, all variables with importance < 5% were removed (Brambilla et al. 2013;Warren et al. 2014). Model accuracy was analysed by the area under the curve of the receiver operating characteristic (ROC) (Pearce and Ferrier 2000;Fawcett 2006). To assess the suitable area for the stone marten, we converted the continuous suitability map into a binary (suitable/unsuitable) classification, using the "equal training sensitivity and specificity" threshold (

Distribution
A total trapping effort of 4539 camera trap-days allowed recording of stone martens 362 times (one video-clip per 12.5 trap-days). The species was recorded at 83% of trapping sites (Table 2), occurring in all available habitats, up to 1850 m above sea level (Fig. 2). The pine marten was recorded for the first time in the north of the protected area in October 2014; throughout the study period, a total of 16 independent events occurred at six different sites (12.5%; Table 2). The use of baits and videos allowed the reliable identification of the two species for 79% of trapping events. For both martens, encounter-rates did not differ among sampling periods (χ 2 = 3.37, 2 df, P = 0.18 and χ 2 = 1.10, 2 df, P = 0.58, respectively). The mesocarnivore community also included the red fox (one video-clip per 28.1 trap-days) and European badger (one video-clip per 135.2 trap-days).
To maximise the cost effectiveness of genetic analyses, the apparently "less fresh" faecal samples-55 out of the 167 faecal samples collected between July and September 2015-were discarded. Of the remaining 112 samples, 96 (83.9%) could be assigned to a Martes species by the PCR-RFLP analysis-namely 82 to the stone marten and 12 to the pine marten-, while 16 samples did not amplify.

Density and kinship as assessed by the genetic sampling
To assess stone marten density, 99 out of 128 faecal samples collected in spring 2017 were selected for microsatellite genetic analyses based on their "freshness". The multiplex screening test was not passed by 55 samples, which were then discarded. Full multilocus microsatellite genotypes were obtained for the remaining 44 samples, of which 41 were assigned to the stone marten, two to the pine marten and one to Mustela sp.
The average proportion of positive PCRs (calculated from correctly and fully genotyped samples) was 82% and varied among loci from 70 to 100%. The observed average error rates across loci were ADO = 0.241 and FA = 0.021, while the number of alleles per locus ranged from 3 to 10 (mean 5.83). Average, non-biased observed and expected heterozygosities (Ho and He) were 0.6 and 0.62 respectively. PID analysis showed that the set of 17 loci would produce an identical genotype with a probability of 1.54 × 10 −9 , and with a probability of 3.43 × 10 −4 for a full-sib. After a regrouping procedure (i.e. pairwise comparison of the different genotypes obtained), we identified 16 stone marten individual genotypes with a complete multi-locus profile. The average number of detections (re-samplings) per individual was 2.6 (min-max 1-6), with six individuals detected only once and two individuals recorded six times. Using capwire estimators, population size was assessed as 21 (CI = 15 -28) individuals; including only the cells for which at least one sample was genotyped (N = 22), density was assessed at 0.95 (0.7-1.3) ind/km 2 .
Except for individuals 2 and 4, all stone martens were related to each other (14 half-sib pairs), individuals 3 and 5, 5 and 8, 6 and 13; and 11 and 16 being first order relatives (i.e. parent/offspring or full-sib dyad). Additionally, using ML-Relate, 13 full-siblings were acknowledged among the 16 stone martens. The two pine martens identified in the area did not result related to each other.
If we assume that mountain ridges that border the valleys coincide with the range limits of males, in the study area it was possible to identify six clusters consisting of 1-3 individuals (Fig. 3).

Habitat suitability
Habitat suitability, as assessed using 161 independent records of the stone marten in the protected area, depended on six environmental variables (Fig. 4), of which four provided the major contributions: rocky grasslands (31.67% importance), rocks and debris (24.63%), beech forests (16.64%) and chestnut forests (10.63%). The first had a non-linear quadratic effect, with the highest suitability at intermediate values of percent cover, although always > 0.4 (Fig. 4). Open rocky areas showed a negative trend, while suitability for the stone marten increased linearly for the other two variables (Fig. 4). The discriminatory ability of the MaxEnt model was sufficient (AUC = 0.74). Following the reclassification analysis based on ETSS = 0.47, 63.9% of the study area resulted suitable for the stone marten (Fig. 5), which, according to univariate analysis, selected beech-and other deciduous forests, while avoided rocky and scree areas (Table 3).

Discussion
Associative modelling approaches, such as SDMs, which consider the target species locations to be representative of ideal habitat conditions (in the multidimensional space described by the chosen variables), assess adequately the species potential distribution only if no other factor plays a major role in determining its occurrences (Gough and Rushton 2000). For the stone marten, interspecific competition with the strictly related pine marten has been  (Delibes 1983) and diet (Gazzola and Balestrieri 2020), suggesting that pine marten dominance may be a major biotic constraint of stone marten's niche.
In the VGNP, the results of both camera-traps and genetic surveys were consistent with the occasional record-based framework, providing straightforward evidence of the widespread occurrence of the stone marten and patchy distribution of the pine marten. If we assume that the frequency of occurrence of records is an index of the relative abundance of both species (Gese 2001;Carbone et al. 2001), the stone marten stood out also in terms of numbers. Low pine marten abundance may depend on the recent recolonisation of this sector of the Alpine range: in the last decade, the number of roadkills has progressively increased and some camera-trapping records have been collected in the western and northern parts of the province (Mosini and Balestrieri 2017), suggesting that the pine marten may be reinforcing its occurrences on the Alps as well as it is expanding in lowland areas of NW Italy (Balestrieri et al. 2015(Balestrieri et al. , 2016b. Whatever the reason, the large numerical difference recorded suggests that currently interspecific competition may play a negligible role in shaping habitat use by the stone marten, although, when sympatric, the pine marten usually dominates in forested habitat as those forming the bulk of our study area. Suitability models confirmed the preference of the stone marten for a mosaic of forested areas, particularly broad-leaved forests, which likely offer both food resources and cover from predators, and rocky grassland areas, which may offer cavities which are both safe and thermal regulated resting sites (Birks et al. 2005;Virgós et al. 2012). Pine martens, which usually nest on trees, in winter, in response to extreme cold, frequently use cavities at ground-level (Brainerd et al. 1995;Zalewski 1997). As the stone marten prefers warm climates (Vergara et al. 2016) and is also less arboreal than the pine marten (Goszczyński et al. 2007), in the study area may find suitable shelter sites in rocky crevices. Cold air temperatures and lack of cover would also explain the negative effect of "rock and debris" on habitat suitability for the stone marten, as this habitat mostly coincides with the top of mountain ridges. Detection probability may vary with the same covariates that affect occurrence probability, leading to biased estimates of their importance in determining occupancy (Yackulic et al. 2013). As, using camera-traps, encounterrates were constantly high throughout the study period and scat-based genetic sampling of stable marten populations has been demonstrated to provide detection probabilities next to 1 (Balestrieri et al. 2015), we are reasonably confident in the performance of our SDMs.
Besides competition with the pine marten, a further factor that has been reported to affect stone marten distribution and abundance is the availability of fleshy fruits (Mortelliti and Boitani 2008;Virgós et al. 2010), the latter being considered the most frugivorous mesocarnivore (Virgós et al. 2012). Nonetheless, recently Gazzola and Balestrieri (2020) reported that frugivory in the stone marten may depend on competition with the pine marten, the first preying mostly on rodents in areas of allopatry. In our study area, small mammals (Clethrionomys glareolus, Glis glis, Apodemus sp.) formed the bulk of stone marten diet (Balestrieri et al. 2018), suggesting that fruit availability did not affect the use of space.
In the last two decades, faecal DNA-based genotyping has proven an effective non-invasive method for estimating population size for several elusive species, including mustelids (e.g. pine marten, Ruiz-González et al. 2013;O'Mahony et al. 2012;Sheehy et al. 2014;Eurasian otter, Arrendal et al. 2007;Vergara et al. 2014). In our study area, mean stone marten density fell within the range recorded in rural Switzerland (Lachat Feller 1993) and was consistent with the density of closely related pine marten in areas on northern Europe showing similar climatic conditions (0.6-0.7 ind/km 2 ; Zalewski and Jedrzejewski 2006). Assuming that topographical constraints drive home range size and shape, boundaries tending to coincide with ridges (Powell and Mitchell 1998;Fig. 4 Response curves of the main land cover variables affecting stone marten distribution in the Val Grande National Park Monterroso et al. 2013), stone marten individuals may be split into six groups, each consisting of 2-3 individuals, in agreement with the species' intra-sexual spacing pattern (Powell 1978;Genovesi et al. 1997).
As recorded for both pine martens (Balestrieri et al. 2016a) and otters (Vergara et al. 2014), most individuals were related to each other, suggesting rapid population renewal and that dispersion occurs on relatively short distances. Heterozygosity matched with average values for the pine marten in continental Europe (min-max 0.56-0.7; Kyle et al. 2003;Pertoldi et al. 2008), while was slightly higher than that reported for the stone marten By applying two non-invasive methods, we pointed out that, as reported for feeding habits (Monterroso et al. 2016), in areas of sympatry with the pine marten patterns of habitat selection depend on the relative abundance of the two competing species. Being dominant, in the study area, the widespread stone marten  selected forested areas, attaining density values like those reported for the pine marten in northern Europe. Further studies are needed to confirm the intra-sexual territorial behaviour and determine how competition with the pine marten affects stone marten density in mountainous areas.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.