Landscape relatedness: detecting contemporary fine-scale spatial structure in wild populations

Methods for detecting contemporary, fine-scale population genetic structure in continuous populations are scarce. Yet such methods are vital for ecological and conservation studies, particularly under a changing landscape. Here we present a novel, spatially explicit method that we call landscape relatedness (LandRel). With this method, we aim to detect contemporary, fine-scale population structure that is sensitive to spatial and temporal changes in the landscape. We interpolate spatially determined relatedness values based on SNP genotypes across the landscape. Interpolations are calculated using the Bayesian inference approach integrated nested Laplace approximation. We empirically tested this method on a continuous population of brown bears (Ursus arctos) spanning two counties in Sweden. Two areas were identified as differentiated from the remaining population. Further analysis suggests that inbreeding has occurred in at least one of these areas. LandRel enabled us to identify previously unknown fine-scale structuring in the population. These results will help direct future research efforts, conservation action and aid in the management of the Scandinavian brown bear population. LandRel thus offers an approach for detecting subtle population structure with a focus on contemporary, fine-scale analysis of continuous populations.


Introduction
Knowledge of contemporary spatial structuring of populations is an important basis for ecological studies in addition to informing and facilitating conservation of a species (Bossart and Prowell 1998;Palsbøll 1999). With the recent advent of high-throughput technologies, studies of spatial structuring have been increasingly considering genetic structure enabling an understanding of how variations in the landscape affect gene flow. Many methods have been developed to study genetic structure and while each method is informative, there are limitations that make certain types of analysis difficult.
Most genetic structure methods are based on assignment tests where individuals are 'assigned' to a subpopulation that is most fitting to their genotype (Manel et al. 2005). In practice, most assignment methods apply a Bayesian clustering algorithm to identify population structure (Guillot and Foll 2009) such as the program Structure (Pritchard et al. 2000). Kinship-based methods are a type of assignment test that identifies and locates highly related individuals as inferred through molecular markers (Broquet et al. 2009;Palsbøll et al. 2010). Other methods include multidimensional scaling (e.g. PCA see Jombart et al. 2010) or spatial autocorrelation and are also based on allele frequencies (see Peakall et al. 2003). One common limitation of assignment methods is that they are based on the island population model (Latter 1973) and thus designed only for discrete populations. Furthermore, most methods readily identify subpopulations that are highly differentiated (i.e. with little gene flow between them), but often fail to detect structure in populations with high connectivity (Manel et al. 2005;Saenz-Agudelo et al. 2009;Lowe and Allendorf 2010). Here, we present a new method designed to detect contemporary, fine-scale spatial structure across the landscape in a continuous population, thereby offering a complementary method to the existing set of methods.

Limitations of assignment-based methods
While assignment-based methods often reveal structuring in a population, assumptions such as random sampling, Hardy-Weinberg equilibrium, and marker neutrality are common (Saenz-Agudelo et al. 2009). Violations of the assumption of marker neutrality are confounded by the variations in selective pressures on different loci, thus influencing estimations of population differentiation (Whitlock and McCauley 1999;Waples and Gaggiotti 2006). Even loci assumed to be selectively neutral can be subjected to the same selective pressures as a result of genetic linkage. Moreover, ascertainment of genetic markers, their representativeness of the whole genome, and total number of loci can have effects on resulting estimates (e.g. Helyar et al. 2010). Markers ascertained either outside the population of interest or using too few individuals to properly assess population-level allele frequencies can bias estimates of structure (i.e. ascertainment bias). Likewise, marker characteristics such as high polymorphism may mask true underlying structure (Putman and Carbone 2014). For one of the most commonly used methods, cluster-based modeling, Putman and Carbone (2014) recommend a minimum of 50 loci for reliable results. However, if the aim is to detect fine-scale structure, prohibitively more than 50 would be required (Peery et al. 2008), thus limiting studies to broad-scale structure. Furthermore, assignment methods often require a priori knowledge of source populations (Waples and Gaggiotti 2006). Finally, two confounding factors in many assignment-based methods are the presence of related individuals (Manel et al. 2005;Putman and Carbone 2014) and isolation by distance (IBD) (Blair et al. 2012;Meirmans 2012;Cushman et al. 2014;Ruiz-Gonzalez et al. 2015), which can result in family structure or allelic clines being confused with population structure.

Kinship-based methods
Several methods exist for estimating structure by using kinship to estimate migration rates (Palsbøll 1999;Peery et al. 2008;Saenz-Agudelo et al. 2009). Parentage methods provide accurate estimates of connectivity between two subpopulations when connectivity between subpopulations is high, which is the opposite of other assignment-based methods, making the two approaches complementary (Saenz-Agudelo et al. 2009). However, sufficient sampling coverage of candidate parents can be challenging. Several variations of kinship methods exist (see Vitalis 2002;Fontanillas et al. 2004;Broquet et al. 2009;Økland et al. 2010;Iacchei et al. 2013;Dharmarajan et al. 2014). In contrast to other assignment methods, kinship-based methods are useful in detecting finescale structuring (Saenz-Agudelo et al. 2009;Palsbøll et al. 2010). Another key advantage is that they estimate contemporary population structure based on the existing generations at the time of sampling (Palsbøll et al. 2010), whereas non-kinship assignment methods based on genetic divergence only provide estimates of structure that has accumulated over time. These estimates are, therefore, averages of processes occurring in the past tens to thousands of generations. Thus, kinship-based methods provide reliable estimates of contemporary structure in discrete populations.

The landscape relatedness (LandRel) method
Our landscape relatedness (LandRel) method is based on pairwise relatedness values interpolated over the landscape. This method is designed for detecting finescale structure in continuous and contemporary populations, but can theoretically be used for discrete populations. As LandRel relies on quantitative relatedness values, there is no distinction of kin categories (e.g. parent-offspring, full-siblings, etc.) thereby differentiating it from kinship-based methods. LandRel relies on the assumption that within a perfectly homogeneous landscape in which a population is randomly mating, the mean relatedness throughout the landscape will be uniform. Any divergence from these two assumptions will result in a heterogeneous distribution of relatedness. Consider, for example, that a fenced highway was built splitting an otherwise panmictic population in two. Individuals born near the road will no longer disperse in random directions, but will disperse alongside or away from the road on either side. This will lead to a sorting of individuals regarding relatedness. LandRel seeks to detect divergence from landscape homogeneity and random mating behaviour by testing for non-uniform distribution of relatedness throughout the landscape.
LandRel shares many of the advantages of kinshipbased methods including the ability to detect fine-scale structuring. However, unlike kinship-based methods, LandRel can be used for continuous populations and does not require a priori knowledge of source populations. Furthermore, since it is based on relatedness, there is no need to remove highly related individuals, which is necessary for some kinship-based methods (e.g. Saenz-Agudelo et al. 2009) and most other assignment methods. Additionally, kinship-based methods, particularly parentage methods, rely heavily on high quality genotyping to infer kinship. LandRel does not require knowledge of kinship making it more tolerant to genotyping errors. Furthermore, LandRel is spatially explicit thereby elucidating where in the landscape structure exists.

Study system
The Scandinavian brown bear (Ursus arctos) consists of three genetically distinct subpopulations with some connectivity between them (Manel et al. 2004;Norman et al. 2013). The two northernmost subpopulations originate from the eastern European lineage while the southernmost subpopulation is a relic of the western European Iberian lineage (Taberlet and Bouvet 1994). The population size within Scandinavia has been steadily increasing (Kindberg et al. 2011) since 1930 when a hunting-induced bottleneck reduced the size to approximately 130 individuals (Swenson et al. 1995). However, the latest population estimate revealed a decline in just a few years (Kindberg and Swenson 2014) prompting Artdatabanken (Uppsala, Sweden) to change the national status to Near Threatened. Furthermore, based on samples from up to 29 years ago, Tallmon et al. (2004) found that the southernmost population had a low effective population size and low immigration and recommended that this population be monitored for signs of inbreeding.
This study empirically tests LandRel on the southernmost subpopulation of the Scandinavian brown bear. This population is currently estimated to consist of approximately 791 individuals (Kindberg and Swenson 2014), which has not changed significantly since the previous census (Kindberg et al. 2009), and is subjected to quota-based hunting. To our knowledge, no structuring within this subpopulation has been previously identified and we therefore assume that it is a panmictic population. However, considering that the brown bear exhibits male-biased dispersal and female philopatry (Støen et al. 2006;Zedrosser et al. 2007;Norman and Spong 2015), we expect some structure among females. Additionally, we assume the landscape to be homogeneous for the purposes of testing the method. Our expectations are thus that males will show uniform levels of relatedness throughout the landscape while females will have areas with significantly high relatedness akin to core areas.

Sample collection and genetic processing
Samples were collected in the autumn of 2012 in the neighbouring counties of Dalarna and Gävleborg in central Sweden (Fig. 1). For the purposes of conducting a census, the Swedish County Administration Board (Länstyrelsen) organised citizen volunteers to collect faeces opportunistically, place them in vials and send them with the GPS coordinates to Bioforsk (now NIBIO) (Norway) for DNA extraction. Further details of the sampling procedure can be found in Kindberg et al. (2011) and the DNA extraction protocol and individual identification in Schregel et al. (2012). DNA extracts from uniquely identified individuals were then sent to our laboratory (SLU, Umeå, Sweden).
DNA extracts from all sampled individuals were SNP (Single Nucleotide Polymorphism) genotyped at 96 loci as described in Norman and Spong (2015). The panel of SNPs included four Y-chromosome markers, three X-chromosome SNPs, four diagnostic mitochondrial markers and 85 autosomal SNPs. SNPs were ascertained de novo in brown bears across the geographic range in Sweden and included the population being analysed in this study (Norman et al. 2013). SNPs were selected to be highly discriminatory with characteristics such as high minor allele frequency and low levels of linkage between SNPs (Norman et al. 2013).
Individuals with multiple sample locations were analysed using the median-centre of all locations to estimate the most probable home range core. Samples were collected in the autumn, when bears tend to remain within or close by their home range , thereby minimizing potential error from wandering bears. All other individuals were analysed using their single GPS point. Further details of the median-centre analysis can be found in Norman and Spong (2015). We first used our genotype data to perform a principal components analysis (PCA) to discern how a commonly used and statistically rigourous method (Patterson et al. 2006) for analysing population structure performed with our data. We used the ''adegenet'' package version 1.4-2 (Jombart and Ahmed 2011) implemented in R (R Development Core 2013) and chose a two-dimensional analysis. Here a matrix of genotype distances was created and subsequently analysed for genetic differentiation based on eigenvalues. If there is clear population structure, points separate into distinguishable clouds. To help determine if there are true clusters, we tested the resulting values in principal component 1 to identify if it was unimodal (i.e. normally distributed), suggesting one cluster, or multimodal, suggesting two or more clusters, using a quantile-to-quantile plot (qqplot) implemented in R (R Development Core 2013).

Relatedness estimates
Estimates of relatedness between pairs of individuals were derived using the Lynch-Ritland coefficient of relatedness (Lynch and Ritland 1999) with the R package ''related'' version 0.8 (R Development Core 2013; Pew et al. 2014). Relatedness values (r-values) reflect the proportion of the genome that is identical by descent between two individuals. First order relatives (i.e. parent-offspring and full siblings) share approximately 50 % of their genome and have an r-value around 0.50. Second-order relatives (i.e. half-siblings and grandparent-grandoffspring) have an average r-value around 0.25 and those that are unrelated have an r-value around 0.00. If a pair is more unrelated than expected under panmictic conditions, the r-value becomes negative. The Lynch-Ritland r-value was chosen as it has been show to outperform others in simulated and empirical studies (Thomas 2005;Csilléry et al. 2006).

Interpolations
Interpolations of relatedness across the study area were conducted using integrated nested Laplace approximations (Rue et al. 2009) with the package INLA (version 0.014) implemented in R (R Development Core 2013). INLA is an approximate Bayesian inference designed for structured latent Gaussian models; a type of additive regression model. INLA uses direct numerical integration to approximate marginal posterior densities as does the Markov Chain Monte Carlo method, but with much less computational time (Holand et al. 2013).

Procedure
Interpolations for mean relatedness were calculated based on one focal individual at a time. The geographic location of each non-focal individual was represented by the pairwise r-value with the focal individual, whereas the focal individual was not represented by any value. Interpolations were repeated in the same manner for each individual (N = 412) so that every individual was included in an interpolation N-1 times. The results of all interpolations were overlaid together, with the overall result being the sum of the individual-based mean relatedness values across the landscape. To determine if there existed any sexspecific patterns of relatedness, the same process was repeated using only males or only females as focal individuals. All individuals, regardless of sex, were still included as non-focal individuals.
Statistical significance was calculated using overall interpolated values at each grid point by dividing mean relatedness with root-mean-square of the standard deviation. We used an alpha level of 0.05 and, being a two-tailed test, normalised values greater than 0.975 or less than 0.025 were considered significant.
Areas of significance were in turn analysed following the same procedure as for the entire study area (global area), but only including the geographic area that contained areas of significance (local area). We use the term ''global'' in the context of ''all-encompassing'' as opposed to the global brown bear population and is used to distinguish the entire study area from the smaller ''local'' areas that were analysed, and that are contained within the global area. For the local areas, only individuals located within defined boundaries were included in the analysis; however, the same pairwise r-values based on the allele frequencies of the entire sample set were used. As with the global area, for each local area we first analysed all individuals and subsequently partitioned males and females. In addition to focusing on areas determined to be significant, we included an area with almost no global significance as a control.
Landscape Ecol (2017) 32:181-194 185 INLA Prior to running the interpolations, relatedness coefficients were first checked to ensure a Gaussian distribution. Interpolations were performed following the guidelines of Lindgren and Rue (2015) using the spatial SPDE-model (Lindgren et al. 2011). A twodimensional mesh was created for each area being analysed. INLA was run with family set to Gaussian and using the following model: where y is the pairwise r-value. We ran other models with different families and including Euclidean distance as a covariate and found that the model above performed best as it resulted in the lowest deviance information criterion (DIC). Mean relatedness and standard deviations were interpolated at every grid point on grids of 150 9 150 for the global area and 100 9 100 for the local areas. Maps were created using the R package lattice v. 0.20-30 (Sarkar 2008).

Inbreeding
Areas that showed statistically significant relatedness patterns in the interpolations (i.e. the local areas) were further investigated for evidence of inbreeding. For each area, we extracted all pairwise relatedness values equal to or exceeding 0.40, to capture all first-order relatives (i.e. parent-offspring and full-siblings). We then performed a Wilcoxon rank sum test between each pairwise area to detect if there was a statistical significance between areas based on an alpha level of 0.05.

Sample collection and genotyping
We obtained DNA extracts and GPS coordinates for 412 individuals: 180 males and 232 females (Fig. 2). All individuals were genotyped successfully at all 96 SNP loci. Mean minor allele frequency of autosomal loci was 0.37 (range 0.13-0.50). The call rate excluding Y-chromosome loci was 0.997 and error rate based on heterozygous loci was 3.8 9 10 -4 .

Population structure with PCA
The PCA resulted in one loose cloud with no apparent substructuring throughout all sampled individuals (Fig. 3a). The qqplot indicated a normal distribution with, perhaps, some outliers at the ends, which could represent migrants (Fig. 3b). The best interpretation of these results is thus that this group of individuals comprises one population with a few distantly related individuals and, therefore, that no substructuring exists.

Interpretation of maps
We generated maps showing interpolations as heat maps in which redder areas have high degrees of relatedness to all individuals in the study area (rvalue [ 0). Green areas indicate a level of unrelatedness that would be expected in panmictic populations (r-value * 0). Bluer areas indicate increasing degrees of unrelatedness, which would be reflective of nonpanmictic conditions (r-value \ 0), which could be caused either by structure in populations or isolationby-distance.
Global scale (entire study area) Figure 4a shows the overall interpolation for the entire study area. Most of the interpolation is blue, thereby indicating that the population deviates from expectations of panmixia, which is likely attributed to effects of isolation-by-distance. Figure 4b shows areas that are statistically significant (i.e. more unrelated to the population as a whole than expected by chance). The northern parts of the study area show large areas of significance. Partitioning males and females results in differing patterns (see Fig. 4c-f). Males show a significant degree of unrelatedness in northern Dalarna and northern Gävleborg relative to the population as a whole, whereas females show a significantly high degree of unrelatedness in northern Gävleborg only. This suggests that individuals in these areas are segregated from the rest of the population, but it does not reveal how individuals in these areas are related to each other. To explore this, we conducted the same analysis but with a focus specifically within these areas (the local scale) (Fig. 5).
Local scale (portions of the entire study area) Control (CA) We first analysed an area with only one small patch of significance at the global scale, which we used as a control (CA). The mean predicted relatedness across CA shows patterns of high relatedness as depicted in orange shades (Fig. 6a). This was in contrast to the results at the global scale and is likely explained by smaller geographic area and lack of isolation-by-distance. We checked for significance for each sex. Males showed no significant areas (Fig. 6b), whereas females showed large patches of significant relatedness (Fig. 6c). Since the brown bear exhibits male-biased dispersal and female philopatry, these results are consistent with our expectations: panmictic conditions for males and core areas with highly related females. Northern Dalarna (ND) As for CA, ND showed patterns of relatedness across the landscape (Fig. 6d). When testing for statistical significance, almost the entire area was significant. This held true for males (Fig. 6e). Females showed significance for a large proportion of the area (Fig. 6f). At the global scale, we found that males were significantly unrelated to the population as a whole (Fig. 5d). At the local scale, we found that they were significantly more related to each other than expected by chance. This indicates that the individuals in ND are segregated from the rest of the population in the study area.
Northern Gävleborg (NG) The predictions in NG follow the same pattern of relatedness across the landscape as other local areas (Fig. 6g). As for CA and unlike in ND, there were some areas of high relatedness that were statistically significant next to areas that were not significant. However, unlike CA, the significance existed for both sexes (Fig. 6h & i) suggesting that males are not panmictic and hence that there is population structuring occurring.

Discussion
We have developed a new spatially explicit method for detecting population structure in a continuous population that we call LandRel. Assuming random mating, a homogenous landscape and non-sex-biased dispersal, populations should be distributed with relatively equal levels of relatedness throughout the landscape.
Where levels of relatedness are higher or lower than expected by chance, it can be inferred that one or more of the above assumptions have been violated and that population structure exists. This is what LandRel seeks to discover. Where a priori knowledge of landscape structure, dispersal characteristics or mating behaviours exist and differ from the above assumptions, adjustments of expectations of relatedness across the landscape can be made accordingly. We empirically tested LandRel in the south-central population of the Scandinavian brown bear, which exhibits male-biased dispersal and female philopatry (Støen et al. 2006;Zedrosser et al. 2007). Our expectations were therefore that males would display equal levels of relatedness throughout the study area, while females would display areas of high relatedness, also known as core areas. Isolation-by-distance was previously tested (Mantel Test) using the same data and resulted in significance for females (p value \ 0.001), but not for males (p value = 0.080) (see Norman and Spong 2015), thus further supporting the assumption that males are panmictic throughout the study area.
In contrast to our expectations, our LandRel results showed evidence of population structure in two areas in the northern part of the study area. This was first identified at the global scale, where individuals in these northern areas were more unrelated to the population as a whole than expected by chance. As we focused our analyses on these areas, we determined that these individuals were also significantly more related to each other than expected by chance. Since females are philopatric, and are expected to form areas of high relatedness (Støen et al. 2005), it is difficult to draw conclusions regarding population structure of females. However, in areas where males appear to be structured, females show similar patterns thereby further suggesting a barrier to gene flow, and thus population structuring for all individuals.
Since the population is continuous beyond the study area, particularly to the north and northwest, processes occurring beyond the study boundary may influence these northern areas identified as segregated. The influence may be in the form of introgression from the distinct subpopulation that exists north of the study area. Still, this does little to explain the striking differences in population-based relatedness between c Fig. 4 The first column shows interpolations of the entire study area (''global'') for pairwise relatedness of a all individuals to all individuals, c males to all, and e females to all. The second column shows areas of statistical significance derived when the cumulated mean over the root-mean-square falls within the alpha level of 0.05. Areas that are significant indicate that individuals in these areas are significantly more (if red) or less (if blue) related to the population as a whole than expected by chance 123 neighbouring individuals in northern Gävleborg. Additionally, we lack information about individuals from the gap between Dalarna and Gävleborg in the north, making it difficult to determine whether the segregated individuals are a part of one continuous population, or if they are, indeed, two differentiated areas. At this time, it is unclear how biologically arbitrary sampling boundaries affect relatedness in the landscape when the population is continuous beyond the sampled boundary.

Implications for brown bear conservation and management
If the identified structure is a result of physical or behavioural isolation from the remaining population, there is a risk of inbreeding, especially if there are few individuals in the isolated area. We tested both areas that appeared to be segregated for signs of inbreeding and found that northern Dalarna had significantly higher relatedness values than the control, suggesting that inbreeding has occurred. Tallmon et al. (2004) expressed their concern for possible inbreeding in this subpopulation after detecting a low effective population size and low immigration rate. As inbreeding can have severe consequences on population fitness (Newman and Tallmon 2006), this evidence for inbreeding in northern Dalarna is of conservation concern. Though further investigation is needed to determine the cause(s) of the structuring, Nellemann et al. (2007) conducted a habitat suitability model for brown bear in an area just to the south of our local northern Dalarna area and suggested that due to the close proximity of human settlements to the west, south and east, the population may become confined. The segregated population in NW Gävleborg did not show a significant difference in relatedness values when compared with the control. However, in our analysis NW Gävleborg contains individuals that are not segregated from the main population thereby possibly balancing out the results. We therefore cannot make any conclusive statement that inbreeding has not occurred in NW Gävleborg. Our results pinpoint parts of the population that are at greater risk for reduced fitness due to isolation. Further investigation is needed to determine the cause of isolation, whether it is, for example, a barrier in the landscape or deviations from known mating or dispersal patterns.

Landscape relatedness (LandRel) method
Besides detecting population structuring, the LandRel method can be informative for other factors affecting populations. In addition to detecting areas with inbreeding, it can provide an indication of where barriers exist in the landscape irrespective of whether it is a natural barrier or a recently developed humaninduced alteration in the landscape. Additionally, if used in a monitoring programme, LandRel will be sensitive to changes within the population with no lag effect, making it an ideal monitoring tool.
One key advantage of LandRel is that it provides insights into contemporary population processes. Some of the most common methods used to analyse population structure are based on genetic differentiation between areas. This gives rise to historically derived population structure going back several to hundreds of generations. While this may be highly informative for understanding the demographic history of a species and how it has evolved, it falls short if more recent population processes need to be investigated. LandRel thus provides more immediate feedback from effects of recent anthropogenic alterations to the landscape for example.
Taking genomic data one step further from allelebased to relatedness-based analysis minimises some complications that arise with other assignment-based  Fig. 4 showing the three ''local'' areas analysed as follows: orange represents the control area (CA), blue represents northern Dalarna (ND), and green represents northern Gävleborg (NG). Males are represented as triangles and females as circles methods. While the r-value itself is defined using estimates of allele frequencies in the population, the resulting comparisons between individuals should not be affected due to its relative nature. This is also true for loci that are not in Hardy-Weinberg equilibrium. Additionally, r-values are not sensitive to loci under selection. Unlike most assignment methods, LandRel will indicate if there is no structure in the population. While many kinship-based methods share these advantages, LandRel is useful for continuous populations and requires no a priori knowledge of the number of source populations unlike kinship-based methods (Iacchei et al. 2013). LandRel has much in common with the spatial autocorrelation method as developed by Peakall et al. (2003). However, LandRel uses pairwise relatedness instead of allele frequencies, thus minimizing issues associated with incorrect allele frequency estimation. In addition, LandRel results in spatially explicit structure across the landscape whereas the spatial autocorrelation method is visualized through correlograms (Peakall et al. 2003;Smouse et al. 2008). Finally, LandRel enabled us to Fig. 6 The first column shows the overall interpolation for the three local areas: a control area (CA); d northern Dalarna (ND); and g northern Gävleborg (NG). The second and third column shows significant areas of relatedness for males and females respectively Landscape Ecol (2017) 32: 181-194 191 detect fine-scale structuring that a commonly used method, PCA, was not able to detect. There are some limitations with the LandRel method that should be considered. First, if there exists two panmictic subpopulations of equal size that are isolated from each other, the initial results from LandRel will not be apparent as there would appear to be an even spread of relatedness throughout. Even so, the overall interpolated mean relatedness values will end up being lower than if the entire area was panmictic due to the pairwise comparisons between individuals residing in opposite groups. However, the same result would occur if there exists isolation-bydistance. This can be investigated further by testing for isolation-by-distance and focusing the interpolation on smaller areas, as we did in this study, since any structure between areas will become more apparent. Sampling characteristics pose further limitations to LandRel, which is an issue for most landscape genetics studies (see Schwartz and McKelvey 2009). Sampling needs to be conducted at a scale large enough to encompass most natal dispersal distances in order for LandRel to be effective. For this, a priori knowledge of dispersal would be beneficial and can be accomplished using the same data needed for LandRel (see Norman and Spong 2015). While sampling intensity is not as important as scale per se, the higher the sampling intensity, the more accurate the interpolations will be. Finally, more research is needed to bring LandRel to its full potential. For example, testing how it works with various degrees of structuring and simulating different sampling intensities would identify important characteristics of the method.
Collecting samples noninvasively is one approach that enables data collection on elusive and rare species from individuals that are living without ever needing to disturb them (Taberlet et al. 1999;Waits and Paetkau 2005;Smith and Wang 2014). Since noninvasive samples collected from the environment are often degraded due to exposure to UV rays and other chemical processes, use of high quality molecular markers that are insensitive to fragmented DNA are recommended. SNPs are particularly useful for degraded DNA since they require only small DNA fragments (65-100 base pairs). Their qualitative nature also makes them less sensitive to false alleles. Additionally, it is possible to genotype using single copy detection, thereby dramatically reducing the occurrence of allelic dropouts and making it possible to accurately genotype samples with very low DNA quantities. In this study, we used SNPs with a call rate of almost 1.00 and a genotyping error rate of less than 0.0004 enabling us to maximise the amount of information extracted from the samples.

Conclusion
In this study, we show that the LandRel method can be used for identifying contemporary, fine-scale population structure within a continuous population of the southernmost Scandinavian brown bear population. It also helped identify areas with higher levels of inbreeding. As inbreeding can have negative consequences on a population, this population should be monitored regularly. LandRel is a complementary method to the many other methods used to identify population structure and has a focus on fine-scale structure across the landscape within a continuously distributed population. The primary advantage of LandRel is that it provides insight into contemporary processes within a landscape context, something that is highly sought after in ecological and conservation-oriented studies.