Next Article in Journal
Diverse Regulation but Conserved Function: SOX9 in Vertebrate Sex Determination
Next Article in Special Issue
Genetic Variability and Population Structure of Polish Konik Horse Maternal Lines Based on Microsatellite Markers
Previous Article in Journal
Novel GRHL2 Gene Variant Associated with Hearing Loss: A Case Report and Review of the Literature
Previous Article in Special Issue
Waves Out of the Korean Peninsula and Inter- and Intra-Species Replacements in Freshwater Fishes in Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microsatellite DNA Analysis of Genetic Diversity and Parentage Testing in the Popular Dog Breeds in Poland

Department of Animal Molecular Biology, National Research Institute of Animal Production, Krakowska 1, 32-083 Balice, Poland
*
Author to whom correspondence should be addressed.
Genes 2021, 12(4), 485; https://doi.org/10.3390/genes12040485
Submission received: 25 February 2021 / Revised: 23 March 2021 / Accepted: 25 March 2021 / Published: 26 March 2021
(This article belongs to the Special Issue Genetic Structure of World Animal Populations)

Abstract

:
There is growing concern that extreme breed standardization contributes to a reduction of the effective population size and high levels of inbreeding, resulting in the loss of genetic diversity in many breeds. This study examined genetic diversity among eight popular dog breeds in Poland and evaluated the effectiveness of a 21-microsatellite (STR) panel recommended by the International Society for Animal Genetics (ISAG) for parent verification. The following breeds were characterized: German Shepherd, Maltese, Irish Wolfhound, Yorkshire Terrier, Biewer Yorkshire Terrier, Golden Retriever, Labrador Retriever, and French Bulldog. STRUCTURE analysis showed breed distinctiveness among all the dog breeds under study. Reynold’s distance ranged between θw = 0.634 and θw = 0.260. The studied breeds showed a medium level of genetic differentiation; the mean number of alleles per locus ranged from 3.4 to 6.6, and the effective number of alleles from 2.1 to 3.5. The mean degree of heterozygosity varied from 49% to 69% and from 47% to 68% for HO and HE, respectively. The population inbreeding coefficient (FIS) indicated an absence of inbreeding in the studied breeds. The average polymorphism information content (PIC) values for most of the breeds were higher than 0.5. The cumulative power of discrimination (PD) for all the markers in all breeds reached high values (close to 1.0), while the probability of identity (PID) was low, ranging between 10−11 and 10−19. The cumulative exclusion probability when the genotypes of one (PE1) and both parents (PE2) are known and showed that the parentage can be confirmed with a probability of 94.92% to 99.95% and 99.78% to 99.9999%, respectively.

1. Introduction

The dog (Canis familiaris) is one of the first animals domesticated by man and has accompanied humans in their day-to-day life for thousands of years. Historically, dogs were used as working animals to herd livestock, hunt, and guard the home, and now they are also treated as companion animals [1,2,3]. Their extensive use is associated with the human selection of dogs for certain phenotypes. Strong selection for desirable traits often causes an irreversible loss of alleles and a reduction of genetic diversity, which may produce undesirable characteristics and health problems [3,4,5]. Around the world, more than 400 breeds of dogs are recognized [1]. The American Kennel Club (AKC) currently registers 197 dog breeds, and the Fédération Cynologique Internationale (FCI) officially recognizes 353 breeds officially [6,7]. Breeding to achieve specific breed standards can lead to a reduction in effective population size and result in increased levels of inbreeding within breeds, resulting in a loss of genetic diversity in many breeds [5,8]. Therefore, it is important to carry out selective breeding while maintaining breed purity and high biodiversity. To this end, the genetic structure of different dog breeds needs to be studied and genetic changes occurring in breeds have to be monitored.
Microsatellite (STR) markers are a well-known effective and powerful tool widely used to investigate the genetic structure and diversity of dog breeds [9,10,11,12,13]. They are also the most important markers used for dog identification and parentage verification [9,14,15]. A panel of 21 markers recommended for parentage testing in domestic dogs by the International Society for Animal Genetics (ISAG) [16] is the panel standardized across multiple laboratories for canine genotyping and used for routine pedigree testing [14]. The analysis of polymorphism was carried out in the following breeds of dogs in Poland: Polish Hunting Dog [13] and Polish Greyhound [17] based on 21 STR recommended by ISAG, and in Borzoi and Tatra Shepherd Dog based on 18 STR [9,17]. No information is available regarding genetic variation among the most popular breeds of dogs in Poland, namely, German Shepherd, Yorkshire Terrier, Golden Retriever, Labrador Retriever, and French Bulldog [18]. In our study, the genetic analysis of these breeds was possible based on data collected as part of pedigree testing conducted at the National Research Institute of Animal Production (NRIAP). Additional analyses included Irish Wolfhound, Biewer Yorkshire Terrier, and Maltese, which are tested in fairly large numbers as part of parentage control. The objective of the study was to examine genetic variation within and among these dog breeds and to evaluate the effectiveness of an STR panel recommended for parentage testing.

2. Materials and Methods

The study was conducted based on the results of analysis of microsatellite polymorphism performed as part of canine pedigree testing at the National Research Institute of Animal Production in 2018–2020. A total of 903 samples were used from eight breeds, namely, German Shepherd (GS, n = 260), Maltese (M, n = 81), Irish Wolfhound (IW, n = 86), Biewer Yorkshire Terrier (BYT, n = 131), Yorkshire Terrier (YT, n = 77), Golden Retriever (GR, n = 48), Labrador Retriever (LR, n = 103), and French Bulldog (FB, n = 117). For all breeds, sampling was a selection of at least 70 unrelated animals (males and females) from different kennels, except for the GR breed for which 48 samples were collected.
This analysis uses the core panel of STR markers recommended by ISAG for individual identification and parentage analysis and a gender identification marker. The following microsatellite markers were used: AHTk211, CXX279, REN169O18, INU055, REN54P11, INRA21, AHT137, REN169D01, AHTh260, AHTk253, INU005, INU030, FH2848, AHT121, FH2054, REN162C04, AHTh171, REN247M23, AHTH130, REN105L03, REN64E19, and Amel locus. DNA was extracted from swabs and blood using the Sherlock AX Kit (A&A Biotechnology, Gdynia, Poland), following the manufacturer’s protocol. The extracts were quantified with a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The STR loci were amplified using Phusion U Hot Start DNA Polymerase (Thermo Scientific, Wilmington, DE, USA), and the PCR reaction was performed on Veriti® Thermal Cycler amplifier (Applied Biosystems, Foster City, CA, USA), using the following thermal profile: 5 min of initial DNA denaturation at 98 °C, followed by 30 cycles of denaturation at 98 °C for 15 s, annealing at 58 °C for 75 s, elongation of starters at 72 °C for 30 s, and final elongation of starters at 72 °C for 5 min. Analysis of the obtained PCR products was performed using an ABI 3130xl capillary sequencer (Applied Biosystems, Foster City, CA, USA). The amplified DNA fragments were subjected to electrophoresis in 7% denaturing POP-7 polyacrylamide gel in the presence of a standard length of 500 Liz and a reference sample. The results of the electrophoretic separation were analyzed automatically using the GeneMapper® Software 4.0 (Applied Biosystems, Foster City, CA, USA).

Data Analysis

Population structure was analyzed using STRUCTURE software version 2.3.3 [19] considering an admixture model with correlated allele frequencies between breeds. The length of the burn in and Markov chain Monte Carlo (MCMC) simulations was 50,000 and 100,000, respectively, in 10 runs for each number of clusters (K) ranging between 4 and 10. The results were exported to STRUCTURE HARVESTER [20] for plotting the likelihood membership coefficient (DeltaK) values so as to determine the most likely number of clusters. Genetic distance was analyzed using Reynolds’s distance—θw [21].
θ w = l u X u Y u 2 2 l 1 u X u Y u ,
where Xu, Yu are allele frequencies from the first and second populations.
Based on this genetic distance, an unweighted pair-group method with averages (UPGMA) dendrograms were constructed to illustrate the similarities between the populations. Observed heterozygosity—HO, expected heterozygosity—HE, and inbreeding coefficient—FIS, were calculated according to Nei and Roychoudhury [22], and Wright [23]. The Hardy–Weinberg equilibrium (HWE) of the 21-STR loci was tested by exact test using an algorithm based on Markov Chain Monte Carlo methods [24]. The genetics parameters were calculated as follows:
  • Polymorphic information content—PIC [25],
  • Power of discrimination—PD [26],
    P D = 1 j = k n k = 1 n p 2 jk
    CPD = 1 i = 1 n 1 P Di ,
    where pjk is the allele frequency j,k for i-locus; CPD is the cumulative power of discrimination.
  • Probability of identity—PID [27],
    P ID = p i 4 + 2 p i p j 2
    CP ID = P ID   ,
    where pi p,j is allele frequencies j,i; CPID is the cumulative probability of identity.
  • Probability of parentage exclusion for each locus, when the genotypes of one and both parents are known (PE1 and PE2) and the cumulative probability of parentage exclusion (CPE) [28],
    PE 1 = i = 1 n p i 2 1 p i 2 + i > j = 1 n 2 p i p j 1 p i + p j 2
    CPE 1 = 1 i = 1 n 1 PE 1 i
    PE 2 = i = 1 n p i 1 p 2 i > j = 1 n p i p j 2 4 3 p i + p j
    CPE 2 = 1 i = 1 n 1 PE 2 i   ,
    where pjk is allele frequency j,k for i-locus; CPE1, CPE2 are cumulative probabilities of identity.
The statistical analysis was performed using the IMGSTAT software, ver. 2.10.1 (2009), which supports the laboratory of the National Research Institute of Animal Production.

3. Results

A total of 185 alleles were detected at 21 STR loci across all breeds. The total number of alleles per locus ranged between 6 for AHTk211 and 12 for AHTk137 and REN169018.

3.1. Breed Relationships

The genetic population structure of each breed was determined based on the admixture level for each individual dog using the correlated allele frequencies model implemented within the STRUCTURE software. The results of Delta K indicated that the optimal number of genetic clusters representing most similar individuals for breeds was at K = 8 (Figure 1). The average proportion of assignment to the cluster above 90% was found for five breeds, i.e., GS (98%), GL (97%), IW and FB (96%), and M (93%). A lower assignment value was found in the BYT breed (88%) and the lowest in the LR and YT (85%).
Reynold’s (1983) genetic distance between eight breeds is summarized in Table 1. The genetic distance was the greatest between GS and IW (θw = 0.634), and the smallest between BYT and YT (θw = 0.260). The UPGMA dendrogram revealed that the YT and BYT were grouped together, and the Irish Wolfhound breed was farthest away from others (Figure 2).

3.2. Diversity Analysis

The greatest number of alleles was noted for the breeds BYT (138 alleles) and M (135 alleles), and the smallest for IW (72 alleles). The mean number of alleles per locus for breeds ranged from 3.4 for IW to 6.6 for BYT, whereas the effective number from 2.1 for WI to 3.5 for M and YT (Table 2).
Out of the 185 alleles found within these breeds, 20 were breed specific (Table 3). The greatest number of private alleles was observed at the locus REN169O18 (four alleles). There were two alleles each at REN247M23 and REN105L03, and one allele each at 12 loci. No breed-specific alleles occurred at the remaining six loci—AHTk211, INU055, REN169D01, INU030, AHTh171, FH2054. The greatest number of private alleles was identified in the GS breed. Specific alleles for this breed occurred at six loci, yet in most cases with low frequencies (<0.09), only two alleles occurred with frequencies of 0.198 and 0.217 for 212 bp allele in REN162C04 and 158 bp in REN169018, respectively. In FB, a higher frequency of 0.12 and 0.218 was noted for the 229 bp allele in REN105L03, and the 156 bp allele in REN169O18, respectively. In the YT and BYT breeds in CXX279, there was a 128 bp allele with a frequency exceeding 16% each, which did not occur in any of the other breeds. The other private alleles exhibited low frequency (less than 3%). No breed-specific alleles were found in the Irish wolfhound. The frequencies of the breed-specific alleles are presented in Table 3.
Estimates of within-breed genetic diversity are summarized in Table 4. The highest average heterozygosity was found for M (Ho = 0.685 and HE = 0.677), YT (Ho = 0.662 and HE = 0.698), and BYT (Ho = 0.661 and HE = 0.658). Among the eight dog breeds considered, the lowest values of heterozygosity were found for IW (Ho = 0.491 and HE = 0.474). The population inbreeding coefficient (FIS) ranged from −0.049 (GR) to 0.053 (YT). The average p-value of HWE for each breed was higher than 0.05 (Table 4).

3.3. Parentage Testing and Individual Identification

The parameters for determining the suitability of the analyzed STR for identification and parentage testing are presented in Table 4. Polymorphism exceeding 0.6 was observed only for five STR (AHT121, FH2054, AHTh171, REN54P11, and REN64E19), while a level over 0.5 was detected for 10 markers. The lowest polymorphism below 0.4 was noted for the rest of the loci. Mean PIC values for the studied breeds varied between 0.414 (for IW) and 0.655 (YT). The mean PD values for individual markers, calculated for all the eight dog breeds together, exceeded 0.8 for AHT121, FH2054, REN54P11, and REN169O18. For the other loci, PD had mean values exceeding 0.7 or close to 0.7 (REN247M23 and INU005) (Figure 3). The power of discrimination for the whole set of 21 STR, for each of the breeds, shows the high values of 0.999999999874911 (IW) to close to 1.0 (BYT, TY).
The cumulative probability of identity for 21 STR loci resulted in values as low as 4.5 × 10−19, 2.3 × 10−19 and 6.8 × 10−19 for M, YT, and BYT breeds, respectively. The higher values of 1.8 × 10−13 and 6.7 × 10−11 were obtained for GS and IW breeds, respectively (Table 4). The panel of 21 microsatellite markers was assessed for their power of exclusion to test parentage in dogs of eight breeds. The probabilities of exclusion were calculated for two hypothetical situations with one parental genotype available (PE1) and two parental genotypes available (PE2). The probability of exclusion for one parent available (PE1) ranged between 0.005 (REN169DO1 in IW) and 0.467 (CXX279 in YT) and when two parents were available (PE2) between 0.047 (REN169DO1 in IW) and 0.642 (REN169DO1 in YT) across different markers and breeds. The cumulative exclusion probability for PE1 and PE2 varied from 0.949242 (IW) to 0.999495 (YT) and from 0.997768 (IW) to 0.999999 (YT), respectively (Table 4).

4. Discussion

As a consequence of selection pressure, management in closed populations, and historical bottlenecks, many dog breeds are exposed to an increase in inbreeding and a loss of diversity [5,29]. It is therefore essential to analyze the genetic structure and evaluate the genetic variation of as many dog breeds as possible and keep track of changes in these breeds. It is also necessary to evaluate the usefulness of DNA markers used for identification and parentage testing of dogs.
With the advancement of science and technology, new genetic markers, such as single nucleotide polymorphisms, have become widely used; however, due to the cost and time of analysis, STR typing is still used in biodiversity and kinship analysis [8,9,10,11,12,13,14,15,17].
We used 21 STRs to determine the genetic population structure of eight canine breeds chosen for this study. Relationships of breeds were analyzed by two approaches—the Reynolds genetic distance, which provides the highest sensitivity for highly divergent populations [21,30], and model-based clustering [19]. The results of STRUCTURE confirmed that the eight dog breeds analyzed could be subdivided into eight genetic clusters. The results for K = 7, for Biewer Yorkshire Terrier and Yorkshire Terrier breeds, showed one cluster (Figure 4).
The two clusters formed by BYT and YT are suggestive of genetic sub-structuring that resulted from using animals belonging to divergent selection lines. The Biewer breed was founded in Germany in 1984 as a result of the selection for white hair genes to produce tricolor dogs. The Biewer Yorkshire Terriers were ultimately recognized as a distinct breed in 1989 by the Allgemeiner Club der Hundefreunde Deutschland (ACH). The Biewer Terrier became the first breed proven to be distinct and unique using genetic testing [31]. This breed became recognized in 2021 by the AKC as the organization’s 197th breed [5]. The estimated genetic distance confirmed the genetic similarity between YT and BYT (θw = 0.26). Both breeds showed genetic closeness to Maltese (θw < 0.4), which is considered, alongside the Black and Tan Terrier and Skye Terrier breeds, as the ancestor of the Yorkshire Terrier. Similar genetic proximity (θw = 0.41) was observed between Maltese and French Bulldog. The genetic distance between the other breeds was θw > 0.50. The differences between GR and LR, which belong to the same group of retrievers and flushing dogs (pointer type), may be due to their different origins. Golden Retriever is a breed developed in the 19th century in Great Britain, while Labrador Retriever dates back from the 18th century and comes from Newfoundland. Labrador Retriever shows similar genetic distance (θw = 0.55) with German Shepherd. Phylogenetically the most distant breed is Irish Wolfhound which originated in Ireland and was used to hunt and protect against wolves. The breed became extinct as wolf numbers decreased and post-19th century it was presumably recreated by dog fanciers [32].
The highest within-breed diversity was characteristic of the M, BYT, and YT breeds, in which the largest number of alleles was identified. In Maltese, 4 out of 135 identified alleles were specific for this breed, but their frequency was low (less than 3%). In contrast, a 128 bp allele at the CXX279 locus, common only to BYT and YT, was present in both breeds with a frequency of more than 16%. For these breeds, the highest degree of heterozygosity (over 65%) was also noted. For Yorkshire Terrier, Ho = 0.662, which is similar to the studies in the UK [8] and the US [33], where Ho for this breed was 0.66 and 0.789, respectively, indicating good levels of diversity in the YT breed. For the other breeds (GS, GR, LR, and FB), the degree of heterozygosity exceeded 50%. The genetic analysis of GS, which was performed based on the same STR panel and using a different 15-STR panel, showed similar values as in our study: Ho = 0.56 and Ho = 0.54, for dogs from the UK and Italian, respectively [8,9]. Tahir et al. [34], who used a 15-STR panel for GS breed in Pakistan, showed HO to be 0.742. In the same study, higher values were also observed for the Labrador Retriever breed (HO = 0.675) than in our study. A degree of heterozygosity below 50% was observed only for IW. This breed was characterized by the lowest total number of identified alleles. Furthermore, only two alleles each were determined in as many as four loci, among which 131 bp (AHTh130) and 216 bp alleles (REN169D01) had a frequency of 0.912 and 0.945, respectively. The mean effective number of alleles per locus was only 2.1. Limited genetic variation may be the result of the origin of this breed from a common founder population involving four founder lines [35]. Preliminary studies for Irish Wolfhound reported by UC Davis Veterinary Genetics Laboratory showed that HO = 0.502 for 33 STR, and when a 21-STR panel was used, HO = 0.483, and this value was similar to that obtained in our study (HO = 0.491) [36]. In the breeds under study, the mean HO and HE values were similar to one another. The inbreeding coefficient for five breeds was negative, but the mean FIS had low negative values (from −0.004 to −0.049), which suggests no inbreeding in these breeds. The lack of observed genetic bottlenecking in any breed, despite a breeding system, may be due to multiple pedigree lines used. The lack of inbreeding was also reported in Bracco Italiano (FIS = 0.061), Tatra Shepherd (FIS = −0.005), six livestock guard dog breeds (average FIS value = 0.024), and in Polish Hunting Dog (FIS = −0.01) [9,11,13,36].
The usefulness of the studied markers for individual identification of dogs was evaluated based on standard PIC, PD, and PID genetic parameters. As reported by Botstein et al. [25], markers with PIC values greater than 0.5 are considered to be informative, values between 0.25 and 0.50 are fairly informative, and values lower than 0.25 are not very informative. In our study, the mean PIC values were greater than 0.5 for most of the studied breeds, similar to many other breeds of dogs investigated [11,13,34,37,38,39,40]. Only in the GS and IW breeds were slightly lower PIC values (0.49 and 0.41, respectively) obtained (Table 4). An immediate measure of the usefulness of the analyzed STR panels for individual identification is the power of discrimination. The higher the power of discrimination of a given STR panel is, the greater the chance that it can be used for individual identification [41]. The PD for the panel of 16 microsatellite markers for Shiba Inu breed was more than 0.999999 [8], while cumulative PD for the 21 markers in all breeds in our study was close to 1.0. The 21-STR panel used for the human gave a similar value of 0.999999999999999999999999999967 [42]. The probability of identity was calculated to assess the suitability of tested panels for individual identification. PID shows the probability with which two unrelated, randomly selected individuals in the population will have the same genotype. It is accepted that CPID values ranging between 10−3 and 10−4 are sufficiently low for the identification of individuals in natural animal populations [43], whereas CPID value estimated only for 15 STR markers in canine amounted to 10−8 [44]. In our study, CPID ranged between 10−11 for IW and 10−19 for M and YT. For all breeds, the obtained low CPID values should be sufficient to distinguish individual dogs. The usefulness of the investigated panel of markers for parentage verification was determined by calculating the probability of exclusion (PE). For a previously used commercial panel of 10 STR, CPE was 0.994 [45]. The use of 17 or 18 STR markers gave a CPE of 0.99998%, and 0.99996%, respectively [11,46]. For all breeds, the recommended panel of 21 STR, used in this study, achieved a CPE1 above 0.99 and CPE2 above 0.999, except for GS (CPE1 = 0.986) and IW breeds (CPE1 = 0.949, CPE2 = 0.998). The values higher than CPE1 of 0.9999 and CPE2 of 0.99999 were obtained for YT, M, and BYT. The highest CPE2 = 0.999999 was observed for YT.

5. Conclusions

Analysis of the microsatellite DNA markers provides valuable information about canine diversity, and a 21-STR panel is an effective tool for individual identification and parentage testing. Our study showed the lowest PIC (0.414), PE1 (0.949), and PE2 (0.998) with the highest PID (6.71 × 10−11) for the Irish Wolfhound breed, illustrating the lower effectiveness of the STRs panel. In contrast, Yorkshire Terrier and Maltese breeds had the highest PIC (0.655 and 0.640), PE1 (0.9995 and 0.9994), and PE2 (0.999999 and 0.9999987) with the lowest PID (2.34 × 10−19 and 4.47 × 10−19). The results suggest the popular breeds in Poland have sufficient diversity relative to other populations that have been studied. The research here provides baseline data for monitoring and managing the breeds.

Author Contributions

Conceptualization, methodology, A.R. and A.P.; investigation, A.R. and A.P.; writing—review and Editing, A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wayne, R.K.; von Holdt, B.M. Evolutionary genomics of dog domestication. Mamm. Genome 2012, 23, 3–18. [Google Scholar] [CrossRef] [PubMed]
  2. Larson, G.; Karlsson, E.K.; Perri, A.; Webster, M.T.; Ho, S.Y.W.; Peters, J.; Stahl, P.W.; Piper, P.J.; Lingaas, F.; Fredholm, M.; et al. Rethinking dog domestication by integrating genetics, archeology, and biogeography. Proc. Natl. Acad. Sci. USA 2012, 109, 8878–8883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pedersen, N.C.; Liu, H.; Leonard, A.; Griffioen, L. A search for genetic diversity among Italian Greyhounds from Continental Europe and the USA and the effect of inbreeding on susceptibility to autoimmune disease. Canine Genet. Epidemiol. 2015, 2, 17. [Google Scholar] [CrossRef] [Green Version]
  4. Keijser, S.F.A.; Fieten, H.; Vos-Loohuis, M.; Piek, C.J.; Anderson, H.; Donner, J.; Scholten, I.; Nielen, M.; Hesselink, J.W.; Van Steenbeek, F.G. Heterozygosity testing and multiplex DNA panel screening as a potential tool to monitor health and inbreeding in a small, closed dog population. Canine Genet. Epidemiol. 2018, 5, 12. [Google Scholar] [CrossRef] [PubMed]
  5. Lampi, S.; Donner, J.; Anderson, H.; Pohjoismäki, J. Variation in breeding practices and geographic isolation drive subpopu-lation differentiation, contributing to the loss of genetic diversity within dog breed lineages. Canine Med. Genet. 2020, 7, 5. [Google Scholar] [CrossRef] [PubMed]
  6. American Kennel Club. Available online: https://www.akc.org/press-center/articles-resources/facts-and-stats/breeds-year-recognized (accessed on 1 January 2021).
  7. The Fédération Cynologique Internationale. Available online: http://www.fci.be/en/Presentation-of-our-organi-sation-4.html (accessed on 27 April 2020).
  8. Mellanby, R.J.; Ogden, R.; Clements, D.N.; French, A.T.; Gow, A.G.; Powell, R.; Corcoran, B.; Schoeman, J.P.; Summers, K.M. Population structure and genetic heterogeneity in popular dog breeds in the UK. Vet. J. 2013, 196, 92–97. [Google Scholar] [CrossRef] [Green Version]
  9. Bigi, D.; Marelli, S.P.; Randi, E.; Polli, M. Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers. Animal 2015, 9, 1921–1928. [Google Scholar] [CrossRef] [PubMed]
  10. Arata, S.; Asahi, A.; Takeuchi, Y.; Mori, Y. Microsatellite loci analysis for individual identification in Shiba Inu. J. Vet. Med. Sci. 2016, 78, 439–441. [Google Scholar] [CrossRef] [Green Version]
  11. Radko, A.; Rubiś, D.; Szumiec, A. Analysis of microsatellite DNA polymorphism in the Tatra Shepherd Dog. J. Appl. Anim. Res. 2017, 46, 254–256. [Google Scholar] [CrossRef]
  12. Turcsán, B.; Tátrai, K.; Petró, E.; Topál, J.; Balogh, L.; Egyed, B.; Kubinyi, E. Comparison of Behavior and Genetic Structure in Populations of Family and Kenneled Beagles. Front. Vet. Sci. 2020, 7, 183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Goleman, M.; Balicki, I.; Radko, A.; Jakubczak, A.; Fornal, A. Genetic diversity of the Polish Hunting Dog population based on pedigree analyses and molecular studies. Livest. Sci. 2019, 229, 114–117. [Google Scholar] [CrossRef]
  14. Van Asch, B.; Alves, C.; Gusmão, L.; Pereira, V.; Pereira, F.; Amorim, A. A new autosomal STR nineplex for canine iden-tification and parentage testing. Electrophoresis 2009, 30, 417–423. [Google Scholar] [CrossRef] [PubMed]
  15. Dimitrijevic, V.; Stevanovic, J.; Savic, M.; Petrujkic, B.; Simeunovic, P.; Milosevic, I.; Stanimirovic, Z. Validation of 10 mi-crosatellite loci for their use in parentage verification and individual identification in the Yugoslavian Shepherd Dog Sharplanina. Ann. Anim. Sci. 2013, 13, 715–722. [Google Scholar] [CrossRef] [Green Version]
  16. ISAG Panel DOG. 2005. Available online: www.isag.us/Docs/2005ISAGPanelDOG.pdf (accessed on 1 July 2005).
  17. Goleman, M.; Balicki, I.; Radko, A.; Rozempolska-Rucińska, I.; Zięba, G. Pedigree and Molecular Analyses in the Assessment of Genetic Variability of the Polish Greyhound. Animals 2021, 11, 353. [Google Scholar] [CrossRef] [PubMed]
  18. Dog Way. Available online: https://dogway.pl/blog/p/najpopularniejsze-rasy-psow-w-polsce/ (accessed on 7 May 2019).
  19. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2000, 155, 945–959. [Google Scholar]
  20. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Res. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  21. Reynolds, J.; Weir, B.S.; Cockerham, C.C. Eestimation of the coancestry coefficient: Basis for a short-term genetic distance. Genetics 1983, 105, 767–779. [Google Scholar] [CrossRef] [PubMed]
  22. Nei, M.; Roychoudhury, A.K. Sampling variances of heterozygosity and genetic distance. Genetics 1974, 76, 379–390. [Google Scholar] [CrossRef]
  23. Wright, S. Evolution and the Genetics of Populations; University of Chicago Press: Chicago, IL, USA, 1978. [Google Scholar]
  24. Guo, S.W.; Thompson, E.A. Performing the Exact Test of Hardy-Weinberg Proportion for Multiple Alleles. Biometrics 1992, 48, 361. [Google Scholar] [CrossRef] [PubMed]
  25. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 1980, 32, 314–331. [Google Scholar] [PubMed]
  26. Kimberly, A.H. Statistical analysis of STR data. Profiles DNA 1998, 1, 14–15. [Google Scholar]
  27. Paetkau, D.; Strobeck, C. Microsatellite analysis of genetic variation in black bear population. Mol. Ecol. 1994, 3, 189–195. [Google Scholar] [CrossRef]
  28. Jamieson, A.; Taylor, S.C.S. Comparisons of three probability formulae for parentage exclusion. Anim. Genet. 1997, 28, 397–400. [Google Scholar] [CrossRef] [PubMed]
  29. Leroy, G. Genetic diversity, inbreeding and breeding practices in dogs: Results from pedigree analyses. Veter. J. 2011, 189, 177–182. [Google Scholar] [CrossRef] [PubMed]
  30. Libiger, O.; Nievergelt, C.M.; Schork, N.J. Comparison of Genetic Distance Measures Using Human SNP Genotype Data. Hum. Biol. 2009, 81, 389–406. [Google Scholar] [CrossRef] [PubMed]
  31. Veterinary Genetics Laboratory, UC Davis, in Collaboration with Dr. Niels C. Pedersen and Staff. Genetic Diversity Testing for Biewer. Available online: https://vgl.ucdavis.edu/canine-genetic-diversity/biewer (accessed on 1 January 2019).
  32. Buldog Francuski. Available online: http://www.piesporadnik.pl/title,pid,45,oid,47,-%20cid,176.html (accessed on 1 January 2020).
  33. Wictum, E.; Kun, T.; Lindquist, C.; Malvick, J.; Vankan, D.; Sacks, B. Developmental validation of DogFiler, a novel multiplex for canine DNA profiling in forensic casework. Forensic Sci. Int. Genet. 2013, 7, 82–91. [Google Scholar] [CrossRef]
  34. Tahir, M.S.; Hussain, T.; Babar, M.E.; Nadeem, A.; Naseer, M.; Ullah, Z.; Intizar, M.; Hussain, S.M. A panel of microsatellite markers for genetic diversity and parentage analysis of dog breeds in Pakistan. J. Anim. Plant Sci. 2015, 25, 351–356. [Google Scholar]
  35. Veterinary Genetics Laboratory, UC Davis, in Collaboration with Dr. Niels C. Pedersen and Staff. Davice Genetic Diversity Testing for Irish Wolfhounds. Available online: https://vgl.ucdavis.edu/canine-genetic-diversity/irish-wolfhound (accessed on 19 August 2019).
  36. Veterinary Genetics Laboratory, UC Davis. Statistics and Breed-Wide Allele Frequency—Irish Wolfhound. Report Issued August 19. Available online: https://vgl.ucdavis.edu/canine-genetic-diversity/irish-wolfhound/stats (accessed on 19 August 2019).
  37. Radko, A.; Słota, E. Application of 19 microsatellite DNA markers for parentage control in Borzoi dogs. Pol. J. Vet. Sci. 2009, 12, 113–117. [Google Scholar] [PubMed]
  38. Ciampolini, R.; Cecchi, F.; Bramante, A.; Casetti, F.; Presciuttini, S. Genetic variability of the Bracco Italiano dog breed based on microsatellite polymorphism. Ital. J. Anim. Sci. 2011, 10, 267–270. [Google Scholar] [CrossRef]
  39. Ichikawa, Y.; Takagi, K.; Tsumagari, S.; Ishihama, K.; Morita, M. Test in based on microsatellite polymorphisms. J. Vet. Med. Sci. 2001, 63, 1209–1213. [Google Scholar] [CrossRef] [Green Version]
  40. Kang, B.-T.; Kim, K.-S.; Min, M.-S.; Chae, Y.-J.; Kang, J.-W.; Yoon, J.; Choi, J.; Seong, J.-K.; Park, H.-C.; An, J.; et al. Microsatellite loci analysis for the genetic variability and the parentage test of five dog breeds in South Korea. Genes Genet. Syst. 2009, 84, 245–251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Eichmann, C.; Berger, B.; Parson, W. Relevant aspects for forensic STR analysis of canine DNA: Repeat-based nomenclature and sensitive PCR multiplexes. Int. Congr. Ser. 2006, 1288, 813–815. [Google Scholar] [CrossRef]
  42. Boonderm, N.; Suriyanratakorn, D.; Sangpueng, S.; Onthong, N.; Nettakul, A.; Waiyawuth, W. Population genetic data of 21 STR markers in Thais of southern border provinces of Thailand. For. Sci. Int. Genet. 2017, 6, 523–525. [Google Scholar] [CrossRef]
  43. Waits, L.P.; Luikart, G.; Taberlet, P. Estimating the probability of identity among genotypes in natural populations: Cautions and guidelines. Mol. Ecol. 2001, 10, 249–256. [Google Scholar] [CrossRef] [PubMed]
  44. Eichmann, C.; Berger, B.; Steinlechner, M.; Parson, W. Estimating the probability of identity in a random dog population using 15 highly polymorphic canine STR markers. Forensic Sci. Int. 2005, 151, 37–44. [Google Scholar] [CrossRef] [PubMed]
  45. Radko, A.; Słota, E.; Kościelny, M. Polymorphism of 10 microsatellites and their usefulness for paternity control in dogs. In Biotechnology, Agriculture and the Food Industry; Zaikov, G.E., Ed.; Nova Scie. Publishers: New York, NY, USA, 2006; pp. 141–144. [Google Scholar]
  46. Dodd, J.N.; Morris, B.G.; Oliveira, D.A.A.; Bernoco, D. DNA testing for parentage verification and individual identification in seven breeds of dogs. Rev. Bras. Reprod. Anim. 2001, 25, 35–41. [Google Scholar]
Figure 1. Delta K values for STRUCTURE analysis of dog breeds obtained by the program Structure Harvester.
Figure 1. Delta K values for STRUCTURE analysis of dog breeds obtained by the program Structure Harvester.
Genes 12 00485 g001
Figure 2. Dendrogram of Reynolds genetic distance between dog breeds by unweighted pair-group method with averages (UPGMA) algorithm.
Figure 2. Dendrogram of Reynolds genetic distance between dog breeds by unweighted pair-group method with averages (UPGMA) algorithm.
Genes 12 00485 g002
Figure 3. Mean values of the polymorphic information content (PIC) and power of discrimination (PD) for the eight dog breeds under study.
Figure 3. Mean values of the polymorphic information content (PIC) and power of discrimination (PD) for the eight dog breeds under study.
Genes 12 00485 g003
Figure 4. STRUCTURE analysis of 21 STR genotypes from all study dogs (903 samples). The samples were grouped by the eight breeds (K = 8). In the case of K = 7, YT and BYT breeds were grouped in one cluster.
Figure 4. STRUCTURE analysis of 21 STR genotypes from all study dogs (903 samples). The samples were grouped by the eight breeds (K = 8). In the case of K = 7, YT and BYT breeds were grouped in one cluster.
Genes 12 00485 g004
Table 1. Reynolds genetic distance (θw) values of eight study dog breeds.
Table 1. Reynolds genetic distance (θw) values of eight study dog breeds.
BreedFBGRLRMGSIWYTBYT
BF0.000
GR0.5400.000
LR 0.4690.5770.000
M0.4090.5010.4590.000
GS0.5190.5960.5490.4590.000
IW0.5740.6060.5880.5390.6340.000
YT0.4300.5030.4650.3640.4840.5000.000
BYT0.4350.5110.4850.3720.5110.5100.2600.0000
Table 2. The number of alleles identified per locus (N), mean number of alleles per locus (A), and mean effective number of alleles per locus (Ae) for each breed.
Table 2. The number of alleles identified per locus (N), mean number of alleles per locus (A), and mean effective number of alleles per locus (Ae) for each breed.
LocusGSMIWBYTYTGRLRFBN
AAeAAeAAeAAeAAeAAeAAeAAe
AHT12171.894.953.194.585.262.472.484.611
AHT13772.064.521.294.173.454.593.274.012
AHTh17172.792.552.451.694.732.773.842.411
AHTh26082.664.031.683.563.741.493.162.510
AHTk21142.853.431.962.951.642.062.142.06
AHTk25361.552.531.962.853.331.852.043.27
CXX27952.843.452.283.583.541.383.152.68
FH205463.385.352.463.963.263.472.784.08
FH284862.262.232.273.153.331.843.952.88
INRA2163.053.842.963.863.542.653.151.48
INU00542.351.532.461.383.253.552.052.89
INU03052.162.621.762.573.232.741.831.87
INU05553.053.841.762.662.642.351.922.07
REN162C0452.461.932.173.473.441.652.131.78
REN169D0131.983.721.152.963.143.673.343.08
REN169O18104.352.442.162.453.263.352.353.112
REN247M2341.553.332.243.443.451.631.342.77
REN54P1152.185.431.974.463.873.762.963.09
AHTh13073.8105.121.383.553.032.663.983.610
REN105L0382.374.742.673.954.054.041.553.010
REN64E1971.272.242.863.273.742.753.254.29
N1251357213813192122106185
A5.96.43.46.66.24.45.85.1
Ae3.33.52.13.23.52.62.62.9
Table 3. Frequency of private alleles in the study breeds.
Table 3. Frequency of private alleles in the study breeds.
LocusAllele (bp)GSMBYTYTGRLRFB
ATH12192 0.031
ATH137127 0.0130
AHTH260256 0.0243
AHTK2532960.0354
CXX279128 0.1641 *0.1753 *
FH2848234 0.0208
INRA21109 0.0185
INU005134 0.0454
REN162C042120.1982
REN169O18156 0.1197
1580.2168
1760.0097
1780.0044
REN247M23274 0.0128
276 0.0146
REN54P11240 0.0062
AHTH130117 0.0123
REN105L03229 0.2179
2450.0079
REN64E19159 0.0115
* Alleles common to the Biewer Yorkshire Terrier (BYT) and Yorkshire Terrier (YT) breeds only.
Table 4. Mean values of genetic parameters assessed for 21 STR loci of the study breeds.
Table 4. Mean values of genetic parameters assessed for 21 STR loci of the study breeds.
BreedHOHEFISp-ValuePICCPDCPIDCPE1CPE2
GS0.54510.55410.01710.48400.49411 *1.80 × 10−130.9859910.9997326
M0.68550.6771−0.01270.49070.63981 *4.47 × 10−190.9994430.9999987
IW0.49110.4743−0.03730.39520.41391 *6.71 × 10−110.9492420.997768
BYT0.66080.6581−0.00410.34250.61661 *6.82 × 10−180.998750.999996
YT0.66230.69810.05330.31500.65451 *2.34 × 10−190.9994950.999999
GR0.59220.5620−0.04900.44140.51351 *3.37 × 10−140.990100.999875
LR0.59540.5886−0.00880.28130.54291 *4.38 × 10−150.993160.999941
FB0.60770.61770.01730.53950.56021 *8.65 × 10−160.996010.999964
* actual value < 1.0, equal to approximately > 0.9999999.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Radko, A.; Podbielska, A. Microsatellite DNA Analysis of Genetic Diversity and Parentage Testing in the Popular Dog Breeds in Poland. Genes 2021, 12, 485. https://doi.org/10.3390/genes12040485

AMA Style

Radko A, Podbielska A. Microsatellite DNA Analysis of Genetic Diversity and Parentage Testing in the Popular Dog Breeds in Poland. Genes. 2021; 12(4):485. https://doi.org/10.3390/genes12040485

Chicago/Turabian Style

Radko, Anna, and Angelika Podbielska. 2021. "Microsatellite DNA Analysis of Genetic Diversity and Parentage Testing in the Popular Dog Breeds in Poland" Genes 12, no. 4: 485. https://doi.org/10.3390/genes12040485

APA Style

Radko, A., & Podbielska, A. (2021). Microsatellite DNA Analysis of Genetic Diversity and Parentage Testing in the Popular Dog Breeds in Poland. Genes, 12(4), 485. https://doi.org/10.3390/genes12040485

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop