The effects of habitat fragmentation on the genetic structure of wild boar (Sus scrofa) population in Lithuania
BMC Genomic Data volume 22, Article number: 53 (2021)
Wild boar (Sus scrofa) is a widely distributed ungulate whose success can be attributed to a variety of ecological features. The genetic variation and population structure of Lithuania’s wild boar population have not yet been thoroughly studied. The purposes of this study were to investigate the genetic diversity of S. scrofa and assess the effects of habitat fragmentation on the population structure of wild boar in Lithuania. A total of 96 S. scrofa individuals collected from different regions of Lithuania were genotyped using fifteen microsatellite loci.
The microsatellite analysis of the wild boars indicated high levels of genetic diversity within the population. Microsatellite markers showed evidence of a single panmictic wild boar population in Lithuania according to STRUCTURE’s highest average likelihood, which was K = 1. This was supported by pairwise Fst values and AMOVA, which indicated no differentiation between the four sampling areas. The results of the Mantel test revealed a weak isolation by distance and geographic diversity gradients that persisted between locations. Motorway fencing and heavy traffic were not an effective barrier to wild boar movement.
There was limited evidence of population genetic structure among the wild boar, supporting the presence of a single population across the study area and indicating that there may be no barriers hindering wild boar dispersal across the landscape. The widespread wild boar population in Lithuania, the high level of genetic variation observed within subpopulations, and the low level of variation identified between subpopulations suggest migration and gene flow between locations. The results of this study should provide valuable information in future for understanding and comparing the detailed structure of wild boar population in Lithuania following the outbreak of African swine fever.
The wild boar (Sus scrofa) is among the most widespread large mammals, with a natural range extending from western Europe and the Mediterranean Basin to the eastern Russian Federation and Japan, and throughout southeast Asia [1, 2]. Owing to this species’ remarkable adaptability, wild boar populations have expanded their geographical range and can be found in a variety of habitats and climates [3, 4]. Successful range expansion and the increasing abundance of wild boar populations are influenced by several factors such as a high ecological plasticity, high reproductive capacity, their ability to adapt to diverse foods , a lack of natural predators  and supplementary feeding . In light of these factors, the main regulatory mechanism for the rapid increase in the size of wild boar populations is wildlife management [7, 8].
Central European wild boar subspecies are also abundantly distributed in Lithuania . The size of the wild boar population has varied over time, and according to monitoring data in 2000 the estimated number of wild boars exceeded 23,000. Since 2001, the population continued to grow by 1500–2000 individuals a year (data of the Ministry of Environment of the Republic of Lithuania, https://am.lrv.lt). Hunting statistics reveal that 176,722 wild boars were hunted in the last 5 years (2014–2019). Hunting serves a population-regulation function and could affect the dispersal behaviour and population structure of this species.
Habitat fragmentation caused by transport infrastructures could have an influence on the wild boar population. Roads and traffic volumes act as barriers for individuals, can hinder migration, and strongly influence road mortality rates and gene flow between populations [10, 11]. There are two motorways that have a traffic volume of over 20,000 vehicles per day in Lithuania (data of the State Enterprise Lithuanian Road Administration, https://lakd.lrv.lt/en/). In Lithuania, Balčiauskas  recorded wildlife killed on roads and found that wild boar accounted for 9.8% of this figure between 2002 and 2007. The fencing of motorways started in 2004 in order to reduce number of wildlife-vehicle accidents in Lithuania. The use of exclusion fences is an effective method for reducing wildlife-vehicle collisions but increases the barrier effect that results in genetic subdivision [13, 14].
Numerous studies have focused on the ecological impact of wild boar, but little research has been performed on their genetic diversity and population structure in Lithuania. Molecular techniques can serve as valuable tools for improving understanding of genetic changes in populations, population structuring and genetic differentiation .
The main purposes of this study were to investigate the genetic diversity of S. scrofa and assess the effects of habitat fragmentation on the population structure of wild boar in Lithuania.
Genetic diversity analysis of wild boar in Lithuania
A total of 96 wild boars were successfully genotyped at 15 microsatellite loci, and total of 147 alleles were detected (Additional file 1). The number of alleles for each locus (NA) ranged from 2 to 13, with an average over all loci and all sample sites of 6.817 (Table 1, Additional file 1). Private alleles distinctive to a specific subpopulation were present in all sampling areas (Table 1). The greatest average number of unique alleles was detected in sampling areas II and III (Table 1). Overall, the heterozygosity values observed across all loci ranged from 0.614 (I) to 0.639 (IV), whereas expected heterozygosity values ranged from 0.651 (I) to 0.684 (IV) (Table 1).
After Bonferroni correction, exact tests for the Hardy-Weinberg equilibrium (HWE) revealed that 11 of the 15 markers were at equilibrium, while the remainder showed a significant departure from HWE in the wild population due to heterozygote deficiency (Additional file 1). At population level, three of the four sampling areas showed a significant deviation from the Hardy-Weinberg equilibrium after Bonferroni correction. According to the Hardy-Weinberg principle, the occurrence of homozygotes in the populations was higher than the expected values. Positive values of Fixation index (F), indicating a deficit of heterozygosity, were observed in all the sampling areas investigated (Table 1).
Genetic differentiation and population structure analysis
The genetic differentiation between subpopulations was established based on a priori grouping that corresponded to sampling areas in different geographical regions. These regions were divided according to their fragmentation by two motorways and calculated by pairwise Fst values (Fig. 1). Nei’s genetic distances (DNei) and Fst analysis indicated low or no genetic differentiation between all pairs of the subpopulations (Table 2). None of the wild boar subpopulation pairs differed significantly from one another.
The two dimensional PCoA plot showed that the first principal coordinate accounted for 6.92% of the total genetic variation, while the second coordinate produced 5.91% of the total genetic variation. The PcoA analysis indicated admixture between individuals from different locations and did not reveal distinct clustering (Fig. 2).
The result of the analysis of molecular variance (AMOVA) revealed that 86% of the variance was found to be between individuals, while 14% came from differences between individuals within population, and 0% was observed among subpopulations (Table 3). Statistical analysis of the fixation index (Fst = 0.004) and analysis of molecular variance revealed no significant genetic differentiation between the wild boar subpopulations (Table 3). Other F-statistics revealed significant values for Fis = 0.145 (p < 0.001) and Fit = 0.148 (p < 0.001). These data indicated that a greater genetic variability in S. scrofa was mainly distributed within individuals (Table 3).
In order to assess the correlation between the pairwise geographical and genetic distances observed across all collection localities, these were plotted as a linear regression and the Mantel permutation test was performed. As shown in Fig. 3, for all subpopulations, the positive correlation between both variables was weak but significant (R2 = 0.1658, P = 0.009). A relatively high regression coefficient and significant association between genetic and geographical distances were obtained among individuals from subpopulations II and III (R2 = 0.4848, P = 0.001), and a weak but significant (R2 = 0.209, P = 0.001) correlation was found among individuals from subpopulations I and III. However, there was no significant association between genetic and geographical distances for individuals in the combined subpopulations III and IV (R2 = 0.0032, P = 0.333), I and II (R2 = 0.0072, P = 0.351) (Additional file 2).
A Bayesian clustering approach was used to identify the hidden population structure of wild boars in Lithuania, and to establish whether the geographical grouping of samples corresponded with genetic groups. The Bayesian analysis of these data indicated an optimal value of ∆K = 4 for the clustering of the samples into four groups (as determined by the Evanno method ; Fig. 4a). However, these four clusters did not correspond to the four a priori designated subpopulations with geographical regions. According to the graphic visualisation of the population structure, no clear geographical clustering of wild boars was observed, and each of four subpopulations had all the genetic clusters with low to high confidence (Fig. 4c and d). The Evanno et al.  ΔK method does not evaluate K = 1 causing an overestimation . As is apparent in Fig. 4b, we also considered mean posterior probability [LnP(D)] to determine the number of clusters and the highest average log-likelihood value was associated with K = 1.
The genetic structure of the wild boar population in Lithuania has not previously been thoroughly studied. This present study showed evidence of a high gene flow between sampling areas with low or absent population differentiation and a lack of heterozygosity within certain loci. The high levels of genetic diversity within the population could be one of the determining factors that lead to species’ stronger environmental adaptability, survivability and numerous other features .
Genetic diversity and variation
The genetic variation patterns were similar to those detected in previously published genetic studies, although the markers used and sample numbers differ greatly between studies. Slightly lower (Ho = 0.622) mean observed heterozygosity values were detected in this study compared with those for wild boar populations in Bulgaria (Ho = 0.63) , Italy (Ho = 0.63), Hungary (Ho = 0.75)  and Portugal (Ho = 0.627) . Furthermore, mean Ho was slightly higher in the Lithuanian population than in the Polish population (Ho = 0.51) . A similar level of heterozygosity (Ho = 0.63) has been reported in wild boar populations inhabiting East Asia . The differences may be due to the use of different marker sets, or to differences in sample numbers. The observed heterozygosities (Ho) were lower than the expected heterozygosities (HE) in four subpopulations (I, II, III and IV) (Table 1). Heterozygote deficiency has also been found in wild boar populations in Bulgaria , Italy, Hungary , Portugal  and Poland . One possible cause of a significant deficit of heterozygosity in most populations, could be a recent demographic expansion [24, 25]. The significant deviations from the Hardy-Weinberg equilibrium observed in three sampling areas showed an excess of homozygotes, which can be attributed to several factors. The primary factors traditionally assumed to account for significant deviations from HWE are null alleles, inbreeding, the Wahlund effect, selection against heterozygotes, population admixture or a combination of these factors [26, 27]. On this basis it was hypothesised that the age-selective harvesting strategy implemented in Lithuania could have an impact on the genetics and sustainability of the wild boar population. In Lithuania, under the law on hunting of the Republic of Lithuania, (law no. 256 of 27 June 2000) traditional hunting practices have changed, with hunters being more motivated to harvest juveniles. The recommended structure for wild boar catches is 70–80% juveniles, 15–20% two-year old and middle-aged boars, and 5–10% adult wild boars . Selective harvesting of wild populations can induce changes in mating systems (such as mate choice), which in turn can induce changes in local gene pools .
Analysis of molecular variance (AMOVA) revealed a high intra-population genetic variation in the wild boar population in Lithuania (Table 3). A similar trend involving genetic variation mainly distributed within the population has also been reported for Bulgarian populations . The high intra-population variability and genetic homogeneity can be influenced by gene flow, which is impacted by the distribution and connectivity of populations [30, 31].
Population pairwise Fst values were effectively zero and none was significant, suggesting that the grouping of Lithuanian wild boar individuals due to traffic barriers had no tangible effect on population structure. Fst analyses found no evidence of genetic differentiation between subpopulations living on opposite sides of the motorway, but a significant positive correlation between genetic and geographical distance suggested isolation by distance. In contrast, the higher genetic differentiation (Fst = 0.0816) obtained from Bulgarian populations occurs due to geographical barriers such as mountain ridges and human impact . One possible explanation for the low FST values observed in the present study is that the wild boar is a migratory species and has a relatively large home range .
Analysis using the Mantel test did not show significant correlation between different sampling areas, except between western and eastern subpopulations (I and II; and III and IV). These results suggested that weak differentiation could occur due to habitat fragmentation by the main motorways: the E67 connecting Vilnius and Klaipėda and the E85 “Via Baltica” motorway connecting Lithuania and neighbouring Poland. The E67 and E85 are the busiest roads in Lithuania (6873 and 9523 vehicles per day respectively). Barriers created due to anthropogenic activity, such as the busiest motorways with high volumes of traffic, fencing and contiguous urban areas, could reduce gene flow and affect the formation of population structure. Conversely, the results of the Mantel test for all subpopulations together demonstrated weak isolation by distance, indicating that geographic distance weakly contributed to the genetic differentiation in the wild boar population.
The multi-locus microsatellite data presented here suggested that the most likely explanation for the lack of genetic structure was that wild boar breeding areas in Lithuania comprise a single panmictic unit. Dispersal of successfully reproducing animals among breeding areas exhibited high genetic connectivity between sampling areas. Bayesian analysis in STRUCTURE showed that all wild boar individuals could probably be assigned to one genetic cluster. This was supported by pairwise FST calculations that demonstrated little (ranging from - 0.000 to 0.007) or no significant differentiation between locations. Nikolov et al. 2009  identified two subgroups in their study of Bulgarian wild boar populations and detected that the Balkan mountain range acts as a natural migration barrier. The results of the present study differed from those of Ferreira et al. 2009 , which showed that Portuguese wild boar have formed three subpopulations (north, centre and south) due to a recent genetic bottleneck. Previous study  has found that strong population structuring exists in the western Iberian Peninsula, in other regions such as Central and Southern Iberia, Central and Eastern Europe or Continental Balkans, wild boar populations seem to be more admixed across wider areas. Similar pattern was revealed in our study. There is an absence of natural barriers in Lithuania, such as high mountains, that could separate subpopulations and explain the genetic difference.
There was limited evidence of population genetic structure among the wild boar, which therefore supported the presence of a single population across the study area. It indicated that there may be no barriers hindering wild boar dispersal across the landscape. The widespread wild boar population in Lithuania, the high level of genetic variation observed within subpopulations, and the low level of variation identified between subpopulations suggested migration and gene flow between locations. The results of this study should provide valuable information in future for understanding and comparing the detailed structure of wild boar population in Lithuania following the outbreak of African swine fever.
Study sites and sampling
Tissue samples of wild boar were collected over a five-year period (2009–2013) from 35 sampling sites across Lithuania (Fig.1). A total of 96 S. scrofa individuals legally harvested by licensed hunters in different parts of Lithuania were investigated. A decision was taken to focus on the single population of Lithuania and a sample size that would be sufficient to characterise population-level genetic diversity when using microsatellites.
The main habitats favoured by boars vary from semi-arid environments to marshes, forests and alpine grasslands . In Lithuania, wild boars mostly prefer habitats of deciduous with spruce and mixed spruce-deciduous forests . Wild boar samples were collected from different regions of Lithuania representing different landscapes. The samples were arranged by grouping individuals into four regional subpopulations (I, II, III, and IV) while also considering the country’s fragmentation by its major roads (E67, E85) with high volumes of traffic (Fig. 1). The first (I) and second (II) sampling areas covered mixed forests and grasslands, deciduous broad-leaved woods were dominant in the third (III) sampling area, and pine Pinus sylvestris forests were prevalent in the fourth (IV) sampling area.
Fresh muscle, spleen and blood were sampled from wild boars and either stored in plastic tubes (5–30 ml) filled with 96% alcohol or kept frozen at a temperature of − 20 °C. All the samples were legally collected and deposited at the State Food and Veterinary Service of the Republic of Lithuania (SFVS). The study did not involve the collection of samples from live animals. An ethics statement was not required. The hunters collected samples in accordance with national regulations on wild boar management.
Amplification and genotyping
In this research, samples were extracted using the “DNeasy Blood and Tissue Kit” (Qiagen, Catalogue No. 69506) following the manufacturer’s instructions. The concentration and purity of the isolated DNA were determined using Nanodrop 2000 Spectrophotometer (Thermo Scientific, DE, USA). Samples were used immediately for amplification or stored at − 20 °C for later use.
A set of 15 microsatellite markers were selected from the list of microsatellite markers recommended by the International Society of Animal Genetics (ISAG) – Food and Agriculture Organization (FAO) . The markers were grouped into two multiplex (SW24, S0386, S0355, SW353, SW936, SW72, S0070, S0107 and S0026, S0155, S0005, SW2410, SW830, SW632, SWR1941) reactions based on their size and annealing temperature.
PCR reactions were carried out in a total volume of 25 μL, containing 1 μL of DNA template, fluorescent forward primer (2 μM) and non-fluorescent reverse primer (2 μM), and 2x QuantiTect Multiplex PCR NoROX Master Mix (Ref. 204,743, QIAGEN GmbH). PCR reactions were carried out in the following steps: 10 min an initial denaturation at 95 °C, 30 or 35 cycles at 95 °C for 30 s depending on the primer set used, annealing at an optimal temperature ranging from 57 to 58 °C, extension at 72 °C for 1 min, then a final extension at 72 °C for 30 min. The ABI 3100 (Applied Biosystems, USA) DNA Analyzer was used to genotype alleles with a GeneScanTM-500 ROX size standard (Applied Biosystems). Gene Mapper 3.7 (Applied Biosystems) software was used to estimate the size of the alleles.
In order to estimate the population genetic structure of wild boars in Lithuania, the number of alleles per locus (NA), observed heterozygosity (Ho) and expected heterozygosity (HE) under Hardy-Weinberg assumptions were obtained in GenAlEx v6.1 . Deviations from the Hardy-Weinberg equilibrium (HWE) were tested with a Markov chain algorithm with 10,000 dememorisation steps, 100 batches and 1000 iterations using Genepop v.4.0 . The P values for HWE were corrected for multiple comparisons by applying a sequential Bonferroni correction, with an initial probability of p = 0.05 . To assess the genetic relationships between subpopulations, pairwise Nei’s genetic distances  were calculated between each pair of the sample sites using the same software. GenAlEx was further used to carry out principal coordinates analysis (PCoA) enabling the visualisation of genetic variation distribution across individuals, analysis of molecular variance (AMOVA) and F-statistics (Fst, Fis, and Fit). FST values were estimated according to Weir and Cockerham’s  version of Wright’s F-statistic using the FSTAT program package , followed by sequential Bonferroni correction for multiple tests .
The Mantel test  was used with 999 permutations in the GenAlEx software to test for evidence of isolation by distance (IBD).
The determination of the most probable number of clusters for Lithuania’s wild boar population (K value) was assessed by the STRUCTURE program version 2.3.4 . The probabilistic method was conducted with 200,000 replications in burn-in and 800,000 replications in the Markov Chain Monte Carlo (MCMC). Twenty clustering simulations (runs) were performed for each possible value of K (K = 1 to K = 10). The outputs of STRUCTURE were submitted to Structure Harvester online program version 0.6.94 (http://taylor0.biology.ucla.edu/structureHarvester/) to estimate the optimal value of K using the Delta (K) method  and allowing for different estimates of K in accordance with Janes et al. . Based on the resulting values of K, a clustering analysis of the studied sampling sites was performed and graphical output generated using CLUMPAK’s main pipeline (http://clumpak.tau.ac.il ).
Availability of data and materials
The datasets generated for this study were submitted to the European Variation Archive (EVA) public project PRJEB48125 (https://www.ebi.ac.uk/ena/browser/view/PRJEB48125). All other data generated or analysed during this study are included in this published article (and its supplementary information files).
African Swine Fever
Analysis of molecular variance
- AP :
Food and Agriculture Organization of the United Nations
- Fis :
- Fst :
Genetic differentiation between sub-populations
- He :
- Ho :
Isolation by distance
International Society of Animal Genetics
Markov chain Monte Carlo
Number of samples
- NA :
Number of alleles
Principal coordinate analysis
State Food and Veterinary Service of the Republic of Lithuania
Oliver WLR, Brisbin ILJ, Takahashi S. The Eurasian wild pig (Sus scrofa). In: WLR O, editor. Status survey and conservation action plan: pigs, Peccaries and Hippos. Gland: IUCN; 1993. p. 112–21.
Ballari SA, Barrios-García MN. A review of wild boar Sus scrofa diet and factors affecting food selection in native and introduced ranges. Mammal Rev. 2014;44(2):124–34. https://doi.org/10.1111/mam.12015.
Sjarmidi A, Gerard J. Autour de la systématique et la distribution des suidés. Monit Zool Ital. 1988;22:415–48.
Massei G, Genov PV. The environmental impact of wild boar. Galemys. 2004;16:135–45.
Johann F, Handschuh M, Linderoth P, Dormann CF, Arnold J. Adaptation of wild boar (Sus scrofa) activity in a human-dominated landscape. BMC Ecol. 2020;20(1):4. https://doi.org/10.1186/s12898-019-0271-7.
Oja R, Kaasik A, Valdmann H. Winter severity or supplementary feeding—which matters more for wild boar? Acta Theriol. 2014;59(4):553–9. https://doi.org/10.1007/s13364-014-0190-0.
Melis C, Szafrańska PA, Jędrzejewska B, Bartoń K. Biogeographical variation in the population density of wild boar (Sus scrofa) in western Eurasia. J Biogeogr. 2006;33(5):803–11. https://doi.org/10.1111/j.1365-2699.2006.01434.x.
Apollonio M, Belkin VV, Borkowski J, Borodin OI, Borowik T, Cagnacci F, et al. Challenges and science-based implications for modern management and conservation of European ungulate populations. Mamm Res. 2017;62(3):209–17. https://doi.org/10.1007/s13364-017-0321-5.
Baleišis R, Bluzma P, Balčiauskas L. Hoofed animals of Lithuania. Vilnius: Akstis; 2003.
van der Ree R, Jaeger JAG, van der Grift EA, Clevenger AP. Effects of roads and traffic on wildlife populations and landscape function: road ecology is moving towards larger scales. Ecol Soc. 2011;16(1):48. https://doi.org/10.5751/ES-03982-160148.
Bennett VJ. Effects of road density and pattern on the conservation of species and biodiversity. Curr Landscape Ecol Rep. 2017;2(1):1–11. https://doi.org/10.1007/s40823-017-0020-6.
Balčiauskas L. Distribution of species-specific wildlife–vehicle accidents on Lithuanian roads, 2002–2007. Est J Ecol. 2009;58(3):157–68. https://doi.org/10.3176/eco.2009.3.01.
Olsson MPO, Widen P. Effects of highway fencing and wildlife crossings on moose Alces alces movements and space use in southwestern Sweden. Wildlife Biol. 2008;14(1):111–7. https://doi.org/10.2981/0909-6396(2008)14[111:EOHFAW]2.0.CO;2.
Wierzchowski J, Kučas A, Balčiauskas L. Application of least-cost movement modeling in planning wildlife mitigation measures along transport corridors: case study of forests and moose in Lithuania. Forests. 2019;10(10):831. https://doi.org/10.3390/f10100831.
Shafer ABA, Wolf JBW, Alves PC, Bergström L, Bruford MW, Brännström I, et al. Genomics and the challenging translation into conservation practice. Trends Ecol Evol. 2015;30(2):78–87. https://doi.org/10.1016/j.tree.2014.11.009.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20. https://doi.org/10.1111/j.1365-294X.2005.02553.x.
Janes JK, Miller JM, Dupuis JR, Malenfant RM, Gorrell JC, Cullingham CI, et al. The K = 2 conundrum. Mol Ecol. 2017;26(14):3594–602. https://doi.org/10.1111/mec.14187.
Jump AS, Marchant R, Peñuelas J. Environmental change and the option value of genetic diversity. Trends Plant Sci. 2009;14(1):51–8. https://doi.org/10.1016/j.tplants.2008.10.002.
Nikolov IS, Gum B, Markov G, Kuehn R. Population genetic structure of wild boar Sus scrofa in Bulgaria as revealed by microsatellite analysis. Acta Theriol. 2009;54(3):193–205. https://doi.org/10.4098/j.at.0001-7051.049.2008.
Vernesi C, Crestanello B, Pecchioli E, Tartari D, Caramelli D, Hauffe H, et al. The genetic impact of demographic decline and reintroduction in the wild boar (Sus scrofa): a microsatellite analysis. Mol Ecol. 2003;12(3):585–95. https://doi.org/10.1046/j.1365-294X.2003.01763.x.
Ferreira E, Souto L, Soares AMVM, Fonseca C. Genetic structure of the wild boar population in Portugal: evidence of a recent bottleneck. Mamm Biol. 2009;74(4):274–85. https://doi.org/10.1016/j.mambio.2008.05.009.
Tajchman K, Drozd L, Karpiński M, Czyżowski P, Goleman M. Population genetic structure of wild boars in Poland. Russ J Genet. 2018;54(5):548–53. https://doi.org/10.1134/S1022795418050101.
Choi SK, Lee JE, Kim YJ, Min MS, Voloshina I, Myslenkov A, et al. Genetic structure of wild boar (Sus scrofa) populations from East Asia based on microsatellite loci analyses. BMC Genet. 2014;15(1):85. https://doi.org/10.1186/1471-2156-15-85.
Cornuet JM, Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996;144(4):2001–14. https://doi.org/10.1093/genetics/144.4.2001.
Leblois R, Estoup A, Streiff R. Genetics of recent habitat contraction and reduction in population size: does isolation by distance matter? Mol Ecol. 2006;15(12):3601–15. https://doi.org/10.1111/j.1365-294X.2006.03046.x.
Pemberton JM, Slate J, Bancroft DR, Barret JA. Nonamplifying alleles at microsatellite loci: a caution for parentage and population studies. Mol Ecol. 1995;4(2):249–52. https://doi.org/10.1111/j.1365-294X.1995.tb00214.x.
Robertson A, Hill WG. Deviations from hardy-Weinberg proportions: sampling variances and use in estimation of inbreeding coefficients. Genetics. 1984;107(4):703–18. https://doi.org/10.1093/genetics/107.4.703.
Janulaitis Z. Distribution, abundance and regulation of wild boar population in Lithuania. Acta Zool Litu. 2003;13(1):88. https://doi.org/10.1080/13921657.2003.10512549.
Allendorf FW, Hard JJ. Human-induced evolution caused by unnatural selection through harvest of wild animals. PNAS. 2009;106(1):9987–94. https://doi.org/10.1073/pnas.0901069106.
Sheng Y, Zheng W, Pei K, Ma K. Genetic variation within and among populations of a dominant desert tree Haloxylon ammodendron (Amaranthaceae) in China. Ann Bot. 2005;96(2):245–52. https://doi.org/10.1093/aob/mci171.
Hopley T, Byrne M. Gene flow and genetic variation explain signatures of selection across a climate gradient in two riparian species. Genes (Basel). 2019;10(8):579. https://doi.org/10.3390/genes10080579.
Veličković N, Ferreira E, Djan M, Ernst M, Obreht Vidaković D, Monaco A, et al. Demographic history, current expansion and future management challenges of wild boar populations in the Balkans and Europe. Heredity. 2016;117(5):348–57. https://doi.org/10.1038/hdy.2016.53.
FAO. Molecular genetic characterization of animal genetic resources. FAO Animal Production and Health Guidelines. 2011. Available at: http://www.fao.org/3/i2413e/i2413e00.htm. Accessed 20 May 2020.
Peakall R, Smouse PE. GENALEX 6: genetic analysis in excel population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95. https://doi.org/10.1111/j.1471-8286.2005.01155.x.
Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Heredity. 1995;86(3):248–9. https://doi.org/10.1093/oxfordjournals.jhered.a111573.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
Nei M. Genetic distance between populations. Am Nat. 1972;106(949):283–92. https://doi.org/10.1086/282771.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70. https://doi.org/10.1111/j.1558-5646.1984.tb05657.x.
Goudet J. FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995;86(6):485–6. https://doi.org/10.1093/oxfordjournals.jhered.a111627.
Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27(2):209–20.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. https://doi.org/10.1093/genetics/155.2.945.
Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. CLUMPAK: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15(5):1179–91. https://doi.org/10.1111/1755-0998.12387.
The authors would like to thank Gabija Kazakevičienė for her technical assistance.
Ethics approval and consent to participate
This study did not require official or institutional ethical approval. Animals were hunted in accordance with the order of the Minister of Environment of the Republic of Lithuania No. 258 (27.06.2000) concerning approval of hunting rules applicable on the territory of the Republic of Lithuania. Samples were collected as part of an infectious animal disease control programme in accordance with the order of the director of the State Food and Veterinary Service of the Republic of Lithuania on the control of contagious infectious diseases in animals, 30 March 2006, No. B1–265, Vilnius, Lithuania.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genetic diversity of wild boar in Lithuania estimated based on polymorphisms in 15 microsatellite loci.
Mantel tests of the relationships among genetic differentiation (Fst values) and geographical distance (km) in different sampling areas. A Mantel test between Fst values and the geographical distance in I and II sampling areas. B Mantel test between Fst values and the geographical distance in I and III sampling areas. C Mantel test between Fst values and the geographical distance in II and III sampling areas. D Mantel test between Fst values and the geographical distance in III and IV sampling areas.
About this article
Cite this article
Griciuvienė, L., Janeliūnas, Ž., Jurgelevičius, V. et al. The effects of habitat fragmentation on the genetic structure of wild boar (Sus scrofa) population in Lithuania. BMC Genom Data 22, 53 (2021). https://doi.org/10.1186/s12863-021-01008-8