- Short Report
- Open Access
Microsatellite markers for identification and parentage analysis in the European wild boar (Sus scrofa)
BMC Research Notesvolume 5, Article number: 479 (2012)
The wild boar (Sus scrofa) is among the most widespread mammal species throughout the old world. Presently, studies concerning microsatellites in domestic pigs and wild boars have been carried out in order to investigate domestication, social behavior and general diversity patterns among either populations or breeds. The purpose of the current study is to develop a robust set of microsatellites markers for parentage analyses and individual identification.
A set of 14 previously reported microsatellites markers have been optimized and tested in three populations from Hungary, Portugal and Spain, in a total of 167 samples. The results indicate high probabilities of exclusion (0.99999), low probability of identity (2.0E-13 – 2.5E-9) and a parentage assignment of 100%.
Our results demonstrate that this set of markers is a useful and efficient tool for the individual identification and parentage assignment in wild boars.
The wild boar, Sus scrofa, is currently one of the most widespread wild mammal species, inhabiting an extensive range of environments . Its domestic form, the pig, is of economic importance, and the present day varieties are the result of multiple domestication events that occurred in different regions . Many studies have employed a wide array of genetic markers aimed at inferring domestication, migration processes, and examining patterns of genetic diversification for both domestic and wild forms. Nevertheless, the peculiar reproductive behaviour and mating system of the wild boar have been poorly characterized . While some authors state that litters are sired by a single father , others suggest the existence of multiple paternity [3, 5], thus pointing to either a polygynous (one male mating with several females) or polygynandrous mating system (both males and females mate with distinct individuals).
During the last decades, advances in genome sequencing and mapping studies have reported thousands of polymorphic neutral markers, such as microsatellites, which have proven to be powerful tools for parental and kinship analysis. In this study, we have chosen to use microsatellite loci given their higher variability (e.g., on average expected heterozygosity (He) is higher than 0.6) and consequently increased power for parentage assignment when compared to the same number of bi-allelic markers such as SNPs [6, 7]. Thus, we constructed a panel of polymorphic microsatellite based on two main criteria: 1) suitability to be combined and amplified in-group (plexes) and 2) an average He value higher than 0.6. We believe that the development and validation of such a panel of markers will provide a valuable tool to assess parental and kinship relationships, thus allowing considerable improvement for studies of mating behavior of both domestic and wild pigs.
Our results confirm the effectiveness of the microsatellite panel for establishment of parentage in wild boars. The probability of identity (PI, the probability of two independent samples having the same identical genotype), using all 14 microsatellites, results in values as low as 2.0E-13, 2.5E-9 and 1.2E-11 for the Hungarian, Portuguese and Spanish populations, respectively. Although the values obtained using all the 14 loci provide low values, similar PI values are obtained by combining only five loci, in the case of the Hungarian and Spanish populations, and six loci in the case of the Portuguese populations (Figure 1). The probability of identity when related individuals are included on the samples (PISibs) is also low, attaining the minimum probabilities for the 14 loci combination of 1.0E-5, 1.6E-4 and 3.1E-5 for Hungary, Portugal and Spain, respectively. The probability of exclusion when both parents are unknown (P1X), when one of the parents is known (P2X) and when the parents are putative (P3X), the maximum probability was different for each of the three populations (Figure 1).
Loci informativeness, heterozygosity levels and exclusion probability values are depicted in Table 1. The probability of finding null alleles is generally negligible, with the exception of three loci (Sw24, S0101, Sw857) in the Hungarian population and one locus on the Spanish population (S0226). Our analyses suggest that in both Hungarian and Spanish populations no microsatellite had significant deviations from Hardy-Weinberg proportions, following Bonferroni correction. In the Portuguese population, however, we found deviations on loci SW857 and SW72 (p < 0.003). Also, after Bonferroni corrections, only four out of 91 combinations exhibited significant deviations from a random association between alleles at different loci.
Assignment probabilities, using COLONY v.2.0, revealed a maternity probability of 1 for all mother/litter combinations. Nevertheless, none of the males captured within the same area of the pregnant females was found to sire any offspring. Nonetheless, the paternal genotypes were inferred and revealed that ten of the litters had one exclusive father while two litters from Hungary, two from Spain and one in Portugal presented genotypes consistent with the simultaneous occurrence of at least two different fathers (multi-paternity). The calculated full- and half-sibs probabilities were also consistent with these findings.
Our results support the usefulness of the described set of microsatellites as a valuable tool for parentage analysis in the wild boar, as the individual identification and power of exclusion levels reveal high power and accuracy (Figure 1). Indeed, due to the high information content (Table 1) of some of the microsatellite markers, it is even possible to obtain a high precision in individual identification and parentage assignment with a subset of only 6 markers.
Although many microsatellites have been described both for domestic pigs and wild boars, some of them are not suitable for assignments since 1) they are not sufficiently informative, thus requiring a higher number of markers to satisfactory results, or 2) present technical constrains, such as difficulty on amplification, scoring or poor performance in a multiplex setting. The microsatellite loci here reported were especially selected to overcome those technical limitations, providing diversity levels which are comparable to the molecular tools reported in previous works concerning parentage analysis [5, 8].
The selected microsatellites present high levels of informativeness, and null allele frequencies were found to be above 10%, even if only in a restricted number of loci, in Hungary and Spain (Table 1). Nevertheless, it should be taken into account that the estimation of null alleles is highly influenced by the sampling of substructured populations. Indeed, previous studies also detected similar deviations in other wild boar populations [4, 9], although in our case these deviations were only found in the Portuguese population. Since our data quality control rules out the possibility of genotyping errors as the possible source of these deviations, the most probable cause for the unbalance on the Portuguese wild boar population may be related to recent demographic fluctuation that certainly may have left a strong mark on the population genetic structure. Nonetheless, the similarity in the expected heterozygosity values (Table 1) obtained in these loci when compared with those reported for different species [9, 10], confirms the reliability of this set of markers and its power of resolution for parentage assignment.
Finally, we believe that the inability to detect any putative father within our samples might result from distinct, non-exclusive, causes ranging from 1) insufficient male sampling, 2) the specificity of the wild boar mating behavior where males do not guard females once copulation ends and 3) the indirect consequences of the hunting process which uses dogs and human noise to direct animals towards a group of hunters, thus unbalancing the opportunity to capture solitary adult males or groups of females [5, 11]. Even though the actual father was not retrieved from our own samples, we were nevertheless able to infer the potential fathers’ genotypes using the genotype reconstruction implemented in COLONY v. 2.0 [12, 13], an invaluable tool for the determination of the wild boar mating system.
A total of 167 tissue samples were collected from wild boars hunted in Portugal (36 males and 5 pregnant females bearing 31 offspring), Spain (17 males and 5 pregnant females bearing 24 offspring) and Hungary (15 males and 5 pregnant females bearing 29 offspring). All samples used in this work came from dead hunted animals. The dead of the animals did not result for this work, but from legal game hunting activities. Samples were collected, at the end of the day, from dead hunted animals after veterinary inspection. The tissue samples were extracted with JETQUICK Tissue DNA Spin Kit (Genomed, GmbH) according to the manufacturer’s protocol. A total of 27 fluorescent-labeled microsatellite markers were initially tested with a small panel of wild boar samples (n = 16), after a thorough selection based on available bibliographic data (ref. [14–19]; Table 2). Some of the analyzed markers were discarded either due to low levels of polymorphism (less than 3 alleles), lack of multiplex assay robustness, insufficient information content and difficulty to amplify or score. Finally, a total of 14 loci were chosen and optimized in two multiplex panels each containing 7 microsatellites (Table 2). PCR amplification were carried in two independent reactions with the same procedure – a total volume of 10 μl containing 10 ng of genomic DNA, 10 mM of primer mix (Table 2), Qiagen Multiplex PCR Master Mix (QIAGEN GmbH, Hilden) and water. The reaction conditions were as follow: (1) an initial denaturation at 95°C for 15 minutes; 2) 10 cycles of 95°C for 30 seconds, 60-56°C (ΔT −0.5°C) for 90 seconds and 72°C for 45 seconds; (3) 22 cycles of 95°C for 30 seconds, 56°C for 90 seconds and 72°C for 45 seconds; (4) 8 cycles of 95°C for 30 seconds, 53°C for 90 seconds and 72°C for 45 seconds, and (5) a final extension step on 72°C for 30 minutes. The samples were then tested in 2% agarose gel and their concentration normalized.
The multiplex products were added to a mixture of Hi-Di™ formamide and size standard (Gene Scan™ 500 LIZ size standard) and run in a 3130 XL Genetic Analyzer (Life Technologies) sequencer. GENE MAPPER v4.0 (Applied Biosystems, USA) software was used to analyze the resulting electropherograms in order to identify the obtained alleles. Data quality assessments were scattered along the process either in the form of negative controls (e.g., to exclude contamination problems), or in a final step consisting of a re-amplification and genotyping of randomly chosen samples (10%) to insure a perfect match in the obtained results .
All the data analyses were performed independently for each population and with the exception of maternity and paternity assignment and null alleles, were accomplished using only the adult individuals. The software GENALEX v. 6.41  was used to calculate the deviations from Hardy-Weinberg Equilibrium proportions (HWE), the probability of identity and the power of exclusion for the loci combinations when both parents are known, when one of the parents is known and for two putative parents. Cervus v. 3.0  software was used to analyze the number of alleles, observed and expected heterozygosity, combined non-exclusion probabilities (for the first parent and second parent, parent pair, identity and siblings identity) and the estimated null allele frequency (including mothers and offspring). Genepop v.1.2  software was utilized to calculate the gametic disequilibrium (1000 dememorization steps, 10 batches, 1000 interactions per batch) between alleles of different loci and for each population individually. A maximum-likelihood method implemented in COLONY v. 2.0 [12, 13] was then used to calculate the assignment probabilities in parentage and sib-ship analyses.
Chen K, Baxter T, Muir W, Groenen M, Schook L: Genetic Resources, Genome Mapping and Evolutionary Genomics of the Pig (Sus scrofa). Int J Biol Sci. 2007, 3: 153-165.
Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, Robins J, Lowden S, Finlayson H, Brand T, Willerslev E: Worldwide Phylogeography of Wild Boar Reveals Multiple Centers of Pig Domestication. Science. 2005, 307: 1618-1621. 10.1126/science.1106927.
Aguilera-Reyes U, Zavala-Páramo G, Valdez-Alarcón JJ, Cano-Camacho H, García-López GI, Pescador-Salas N: Multiple mating and paternity determinations in domestic swine (Sus scrofa). Anim Res. 2006, 55: 409-417. 10.1051/animres:2006024.
Poteaux C, Baubet E, Kaminski G, Brandt S, Dobson FS, Baudoin C: Socio-genetic structure and mating system of a wild boar population. J Zool. 2009, 278: 116-125. 10.1111/j.1469-7998.2009.00553.x.
Delgado R, Fernández-Llario P, Azevedo M, Beja-Pereira A, Santos P: Paternity assessment in free-ranging wild boar (Sus scrofa) - Are littermates full-sibs?. Mammalian Biology - Zeitschrift fur Saugetierkunde. 2008, 73: 169-176. 10.1016/j.mambio.2007.07.008.
Bruford MW, Bradley DG, Luikart G: DNA markers reveal the complexity of livestock domestication. Nat Rev Genet. 2003, 4: 900-910. 10.1038/nrg1203.
Schlotterer C: The evolution of molecular markers [mdash] just a matter of fashion?. Nat Rev Genet. 2004, 5: 63-69. 10.1038/nrg1249.
Hampton J, Pluske JR, Spencer PBS: A preliminary genetic study of the social biology of feral pigs in south-western Australia and the implications for management. Wildl Res. 2004, 31: 375-381. 10.1071/WR03099.
Chang WH, Chu HP, Jiang YN, Li SH, Wang Y, Chen CH, Chen KJ, Lin CY, Ju YT: Genetic variation and phylogenetics of Lanyu and exotic pig breeds in Taiwan analyzed by nineteen microsatellite markers. J Anim Sci. 2009, 87: 1-8. 10.2527/jas.2008-1738.
Li K, Geng J, Qu J, Zhang Y, Hu S: Effectiveness of 10 polymorphic microsatellite markers for parentage and pedigree analysis in plateau pika (Ochotona curzoniae). BMC Genet. 2010, 11: 101-
Fernández-Llario P, Mateos-Quesada P: Population structure of the wild boar (Sus scrofa) in two Mediterranean habitats in the western Iberian Peninsula. Folia Zoologica. 2003, 52: 143-148.
Wang J: Sibship Reconstruction From Genetic Data With Typing Errors. Genetics. 2004, 166: 1963-1979. 10.1534/genetics.166.4.1963.
Jones OR, Wang J: COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010, 10: 551-555. 10.1111/j.1755-0998.2009.02787.x.
Rohrer GA, Alexander LJ, Keele JW, Smith TP, Beattie CW: A Microsatellite Linkage Map of the Porcine Genome. Genetics. 1994, 136: 231-245.
Ellegren H, Chowdhary BP, Fredholm M, Hoyheim B, Johansson M, Nielsen PB, Thomsen PD, Andersson L: A physically anchored linkage map of pig chromosome 1 uncovers sex- and position-specific recombination rates. Genomics. 1994, 24: 342-350. 10.1006/geno.1994.1625.
Alexander L, Troyer D, Rohrer G, Smith T, Schook L, Beattie C: Physical assignments of 68 porcine cosmid and lambda clones containing polymorphic microsatellites. Mamm Genome. 1996, 7: 368-372. 10.1007/s003359900106.
Ernst CW, Robic A, Yerle M, Wang L, Rothschild MF: Mapping of calpastatin and three microsatellites to porcine chromosome 2q2·1-q2·4. Anim Genet. 1998, 29: 212-215. 10.1111/j.1365-2052.1998.00319.x.
Fredholm M, Wintero AK, Christensen K, Kristensen B, Nielsen PB, Davies W, Archibald A: Characterization of 24 porcine (dA-dC)n-(dT-dG)n microsatellites: genotyping of unrelated animals from four breeds and linkage studies. Mamm Genome. 1993, 4: 187-192. 10.1007/BF00417561.
Ellegren H, Chowdhary BP, Johansson M, Marklund L, Fredholm M, Gustawon I, Andersson L: A Primary Linkage Map of the Porcine Genome Reveals a Low Rate of Genetic Recombination. Genetics. 1994, 137: 1089-1110.
DeWoody J, Nason JD, Hipkins V: Mitigating scoring errors in microsatellite data from wild populations. Molecular Ecology Notes. 2006, 6: 951-957. 10.1111/j.1471-8286.2006.01449.x.
Peakall ROD, Smouse PE: Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.
Kalinowski ST, Taper ML, Marshall TC: Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007, 16: 1099-1106. 10.1111/j.1365-294X.2007.03089.x.
Raymond M, Rousset F: GENEPOP (Version 1.2): Population Genetics Software for Exact Tests and Ecumenicism. J Hered. 1995, 86: 248-249.
This work was funded by Fundação para a Ciência e Tecnologia (FCT) project PTDC/CVT/68907/2006. AB-P and NM were funded by FCT, Programa Ciência 2009. We thank to all hunters and game managers for assisting us with sampling.
The authors have no competing interests to declare.
VC carried out the microsatellite genotyping and drafted the manuscript. JP-G performed part of the statistical tests and sampling. PS, PF-L, JC, JB, GV, AI, AZ, organized the sampling campaigns. AB-P, PS, and NM participated in the design of the study and performed the statistical analysis. AB-P conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.