Estimation of genetic diversity in Gute sheep: pedigree and microsatellite analyses of an ancient Swedish breed
© The Author(s) 2017
Received: 7 September 2016
Accepted: 20 January 2017
Published: 30 January 2017
Breeds with small population size are in danger of an increased inbreeding rate and loss of genetic diversity, which puts them at risk for extinction. In Sweden there are a number of local breeds, native breeds which have adapted to specific areas in Sweden, for which efforts are being made to keep them pure and healthy over time. One example of such a breed is the Swedish Gute sheep. The objective of this study was to estimate inbreeding and genetic diversity of Swedish Gute sheep.
Three datasets were analysed: pedigree information of the whole population, pedigree information for 100 animals of the population, and microsatellite genotypes for 94 of the 100 animals. The average inbreeding coefficient for lambs born during a six year time period (2007–2012) did not increase during that time period. The inbreeding calculated from the entire pedigree (0.038) and for a sample of the population (0.018) was very low. Sheep were more heterozygous at the microsatellite markers than expected (average multilocus heterozygosity and Ritland inbreeding estimates 1.01845 and -0.03931) and five of seven microsatellite markers were not in Hardy Weinberg equilibrium due to heterozygosity excess. The total effective population size estimated from the pedigree information was 155.4 and the average harmonic mean effective population size estimated from microsatellites was 88.3. Pedigree and microsatellite genotype estimations of inbreeding were consistent with a breeding program with the purpose of reducing inbreeding.
Our results showed that current breeding programs of the Swedish Gute sheep are consistent with efforts of keeping this breed viable and these breeding programs are an example for other small local breeds in conserving breeds for the future.
Gute sheep are an ancient breed from the Swedish island and province of Gotland, belonging to the North European short-tailed sheep group. Both Gute sheep and in general North European short-tailed sheep are characterized in part by their hardiness as well as coat colour and pattern variation . Many of the North European short-tailed sheep breeds are decreasing in population size [1, 2] and have low heterogeneity . However these breeds are unique and their conservation is relevant because they add to overall species diversity [3, 4]. The Gute sheep population experienced a severe bottleneck when polled sheep started to become more popular on Gotland and the horned sheep became rare. Beginning around 1920 the remaining horned sheep on Gotland were gathered by a few individuals interested in preserving the horned sheep. These sheep are the origins of the existing Gute sheep population today (personal communication with Gute sheep breed organizations in Sweden). Nevertheless, Gute sheep, unlike some other North European short-tailed sheep breeds, are not endangered at the moment, although the population is believed to be decreasing in size . Thanks to conservation efforts, there were an estimated 5200 animals in Sweden in 2012 . The first association working with the conservation of the Gute sheep, Föreningen Gutefåret, was formed in 1977, and an additional association, GutefårAkademin, was formed in 2007. Gute sheep breeding is focused on preserving the breed and not on improving production traits.
Gute sheep appear in the scientific literature, having been studied for both coat colour and population structure and genetic diversity. Classical genetic studies from the 1970s took advantage of the coat colour variation in Gute sheep, also then referred to as Goth sheep in English, to study coat colour inheritance in sheep . More recently, there have been population genetic studies which have included Gute sheep [3, 6], but these studies only had a small number of individuals; 20 and 12 respectively. There have been no studies of Gute sheep genetic diversity using pedigree information.
The objectives of this study were to estimate inbreeding and genetic diversity of Swedish Gute sheep using pedigree data and microsatellite genotypes to identify; if the population is at potential risk due to high or increasing inbreeding and decreased genetic diversity; and if strategies need to be implemented to conserve the Swedish Gute sheep in the future.
Pedigree information for all registered Gute sheep in Sweden born from 1960 to 2014 was provided by Elitlamm, the recording program for sheep in Sweden. For the early records only a small proportion of the population were registered in ElitLamm, but in recent years a very large proportion of the Gute sheep population in Sweden has been registered in ElitLamm because animals must be registered in order to be entitled to financial support for the keeping of endangered breeds. Pedigree information included animal identification numbers, year of birth and parent identification numbers. The animal identification numbers were anonymous.
Blood samples and phenotypic records (tail length and information about horns) were collected from different parts of Sweden by a trained technician in collaboration with the association “Föreningen Gutefåret”, one of the two Gute sheep breed organizations in Sweden. Ethical permission to collect blood samples from sheep and use for genetic studies was approved before the study (Dnr C102/13). Samples were stored in the -80 °C freezer of the SLU biobank. DNA was extracted from blood samples for 100 Gute sheep from a total of 13 flocks using the DNeasy mini kit for the QIAsymphony robot (Qiagen®, Hilden, Germany). Sheep were genotyped for a total of eight genetic markers; AME, a SNP marker in the amelogenin gene for sex determination; and seven microsatellite markers (INRA005, INRA023, INRA063, INRA172, MAF214, MAF65 and McM527). These markers are part of the ISAG kit for parentage testing. Genotyping was done at the certified lab at the Department of Animal Breeding and Genetics, SLU. The pedigree information for the 100 animals included flock identification, date of birth and official animal identification numbers for each animal, and their ancestors (which did not correspond with the anonymous animal identification numbers of the pedigree of the whole population).
Pedigree information of the whole population was used to estimate population parameters. Generation intervals, effective number of founders, effective number of ancestors, effective number of founder genomes and marginal contributions of ancestors were estimated with the software Pedig . Generation interval for sires of sires, dams of sires, sires of dams and dams of dams were calculated separately. The average was calculated from these four values and used as the overall generation interval. The coefficient of inbreeding by birth year , pedigree completeness indices, number of complete generation equivalents and average coancestry within birth cohorts were estimated with software, EVA . Effective population size was calculated as one divided by double the annual rate of inbreeding multiplied by the generation interval .
The pedigree information for each of the 100 sheep with blood samples available were combined into one pedigree. This pedigree was used for calculating average coefficient of inbreeding and pedigree completeness with the software EVA  for comparison with results for these parameters from the whole pedigree.
Diversity within populations (1-QINTRA) and between populations (1-QINTER), inbreeding in individuals compared with the population inbreeding (FIS), and Hardy Weinberg tests for heterozygosity deficiency and excess in microsatellite markers were estimated with the program Genepop v.4.3 . Inbreeding and relatedness among individuals  were estimated using the microsatellite genotypes using software Coancestry v.18.104.22.168 . Multilocus heterozygosity was calculated as the proportion of typed loci at which an individual was heterozygous divided by the population heterozygosity at those same loci . Past changes in effective population size (average harmonic mean effective population size) were estimated using microsatellite genotypes in the software VarEff v.1.2  with an assumed mutation rate of 0.00013  and generation interval 3 years. Population structure was estimated with the microsatellite genotypes using principal components analysis (PCA) using R package adegenet .
Generation interval (L) and effective population size (Ne) for Gute sheep in Sweden between 2007 and 2012 calculated using pedigree
Generation interval (L) in years
Effective population size (Ne)
Number of founders
Effective number of founders
Effective number of founder genomes
Effective number of ancestors
We present the first analyses of genetic diversity of the Gute sheep population in Sweden. Pedigree and microsatellite genotype estimations of inbreeding were consistent with a breeding program aiming at reducing inbreeding.
Previous studies have suggested that Gute sheep is a unique breed with low heterogeneity [3, 4], which we could not verify in our study. Tapio et al. estimated inbreeding in 32 north European sheep breeds, including Gute sheep, using microsatellite markers . All breeds in their study had greater FIS estimates than our estimate for Swedish Gute sheep and many of the FIS estimates indicated inbred populations (including for their Gute sheep population) . The heterogeneity in our study was higher than expected in the population, both seen in the multilocus heterozygosity (1.01845) and the fact that five out of seven microsatellite markers were not in Hardy Weinberg equilibrium because of excess heterozygosity. In addition, inbreeding calculated from the entire pedigree, from the pedigree of the 94 individuals and from their microsatellite markers was low. The following could have contributed to differences in results: our diversity and inbreeding estimates are from animals born after animals studied by Tapio et al. . Furthermore, we used pedigree data from the entire registered Swedish population and genotyped 94 animals from 13 different flocks with microsatellite genotypes while only 20 animals had been genotyped in the previous study . Finally, we only used seven microsatellite markers compared to the 25 used by Tapio et al.  and this is a weakness of our study.
Inbreeding estimated based on data from the entire pedigree was more than twice the estimated inbreeding from the sample population (0.038 versus 0.018) which could be due to the amount of pedigree information available for each animal. The entire pedigree is more complete and has more depth: the average depth of the entire population pedigree in generations was 6.06 compared in contrast with 5.41 for the sample population. The founders of a pedigree are assumed to be unrelated and when information from more generations is available, there could be more relationships accounted for. Nevertheless, both estimates of inbreeding were low and our hypothesis is, that this could be because neither of the datasets could account for bottlenecks occurring earlier in time (ie. Before pedigree information was collected). However, the estimation of inbreeding from microsatellite genotypes of the sample population, should be able to account for historical inbreeding, but the estimation of inbreeding from microsatellite data was also low (Ritland inbreeding = -0.03931). We are therefore confident, that our estimation of a low inbreeding reflects reality.
This was also verified by the further analysis, which indeed indicated the existence of bottlenecks in the population. The difference between the effective number of founders and the effective number of ancestors indicated that bottlenecks have occurred in this population and the difference of half the effective population size and the effective number of founder genomes indicated that drift has been accelerating since the founder population.
Previous studies found heterozygote deficiency in European sheep breeds (including North European short-tailed breeds), which was thought to be due to subdivision among flocks creating the Wahlund effect and due to non-random mating [3, 18–21]. In contrast we found that there was excess heterozygosity in Swedish Gute sheep. Excess heterozygosity in microsatellites in a population can be due to the structure of the breeding program, where unrelated animals are bred together as typically found in a population within a conservation program. For the sample population genotyped, 95% of animals had parents with different flock ID numbers indicating that they were born in different flocks, consistent with the breeding of animals that are unrelated. The use of pedigree information in a conservation breeding program has been shown to be a very powerful tool for maintaining genetic diversity and low inbreeding  and our study showed similar results.
Heterozygosity and inbreeding calculated from microsatellite markers were highly correlated with each other but were not with inbreeding calculated from the pedigree. Regardless, we can conclude from these results that there was low inbreeding. Multilocus heterozygosity is not a good indicator of inbreeding calculated from a pedigree because it is not highly correlated unless there is a large amount of variation in inbreeding in a population  and heterozygosity can be lost at a different rate than inbreeding is gained in a population . Additionally, correlation between these estimates can be lower because of amount and quality of information in the recorded pedigree. In a study of Finnsheep, estimates of inbreeding were correlated, however, estimates were not as correlated in one subpopulation of Finnsheep (grey individuals) and the authors believed this was in part because of their less complete pedigree .
The effective population size estimated from microsatellites was 88.3, which was lower than the effective population size calculated from the whole pedigree (155.4). However both of these estimates fall well within or even above the recommended effective population size for small populations, which should be at least between 50 and 100 [26, 27]. On the other hand it is still important to maintain the conservation breeding program and monitor progress because in order to maintain genetic variation in the long term, an effective population size of at least 500 is needed . Nevertheless, the comparison of our results with that of earlier studies suggests that the current conservation strategy in the Swedish Gute sheep population is successful and should be maintained.
The Swedish Gute sheep population can be an example of conservation for other small local breeds.
Despite evidence of historical bottlenecks, inbreeding estimates from the pedigree of all registered Gute sheep in Sweden and from the pedigree and microsatellite genotypes of a sample of the population indicated that the level of inbreeding was low and that inbreeding did not increase over the last six years of data analysed. While the Swedish Gute sheep population has an effective population size that is good for small populations, the current conservation program and monitoring of the average inbreeding should continue.
The project was performed in collaboration with the Föreningen Gutefåret association. The late professor Anna Näsholm was involved in the initial planning of the project.
The collection of blood samples and genotyping and part of the working time for analyses was financed by a grant from the Swedish Board of Agriculture to the Föreningen Gutefåret association. CMR benefited from a joint grant from the European Commission and the Swedish University of Agricultural Sciences, within the framework of the Erasmus-Mundus joint doctorate “European Graduate School in Animal Breeding and Genetics”.
Availability of data and materials
The microsatellite genotypes (Additional file 2) and the horn and tail phenotypes (Additional file 3) of the 94 genotyped Gute sheep are available as additional supporting files. Please contact authors for access to anonymized pedigree data.
CMR analysed the data and wrote the manuscript. AMJ conceived the project and assisted in manuscript writing. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethical permission to collect blood samples from sheep and use for genetic studies was approved before the study (Dnr C102/13).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Dýrmundsson ÓR, Niżnikowski R. North European short-tailed breeds of sheep: a review. Animal. 2010;4:1275–82.View ArticlePubMedGoogle Scholar
- Domestic Animal Diversity Information System (DAD-IS), Food and Agriculture Organization of the United Nations. 2016. http://dad.fao.org/. Accessed 29 Aug 2016.
- Tapio M, Tapio I, Grislis Z, Holm LE, Jeppsson S, Kantanen J, et al. Native breeds demonstrate high contributions to the molecular variation in northern European sheep. Mol Ecol. 2005;14:3951–63.View ArticlePubMedGoogle Scholar
- Tapio M, Ozerov M, Tapio I, Toro M, Marzanov N, Ćinkulov M, et al. Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia. BMC Genet. 2010;11:76.View ArticlePubMedPubMed CentralGoogle Scholar
- Adalsteinsson S, Lauvergne JJ, Boyazoglu JG, Ryder ML. A possible genetic interpretation of the colour variation in the fleece of the Gotland and Goth sheep. Ann Genet Sel Anim. 1978;10:329–42.PubMedPubMed CentralGoogle Scholar
- Mukiibi R, Rochus CM, Andersson G, Johansson AM. The use of endogenous retroviruses as markers to describe the genetic relationships among local Swedish sheep breeds. Anim Genet. 2015;46:220–3.View ArticlePubMedGoogle Scholar
- Boichard D. Pedig: A Fortran package for pedigree analysis suited for large populations. Montpellier: 7th World Cong Genet Appl Livest Prod; 2002. p. 28–9.Google Scholar
- Meuwissen THE, Luo Z. Computing inbreeding coefficients in large populations. Genet Sel Evol. 1992;24:305–13.View ArticlePubMed CentralGoogle Scholar
- Berg P, Nielsen J, Sørensen MK. EVA: Realized and predicted optimal genetic contributions. Belo Horizonte: 8th World Cong Genet Appl Livest Prod; 2006. p. 246p.Google Scholar
- Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Burt Mill, Harlow: Longman Scientific & Technical; 1996.Google Scholar
- Rousset F. Genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.View ArticlePubMedGoogle Scholar
- Ritland K. Estimators for pairwise relatedness and individual inbreeding coefficients. Genet Res. 1996;67:175–85.View ArticleGoogle Scholar
- Wang J. Coancestry: a program for simulating, estimating and analysing relatedness and inbreeding coefficients. Mol Ecol Resour. 2011;11:141–5.View ArticlePubMedGoogle Scholar
- Coltman DW, Pilkington JG, Smith JA, Pemberton JM. Parasite-mediated selection against inbred Soay sheep in a free-living, island population. Evolution. 1999;53:1259–67.View ArticleGoogle Scholar
- Nikolic N, Chevalet C. Detecting past changes of effective population size. Evol App. 2014;7:663–81.View ArticleGoogle Scholar
- Crawford AM, Cuthbertson RP. Mutations in sheep microsatellites. Genome Res. 1996;6:876–9.View ArticlePubMedGoogle Scholar
- Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5.View ArticlePubMedGoogle Scholar
- Lawson Handley LJ, Byrne K, Santucci F, Townsend S, Taylor M, Bruford MW, Hewitt GM. Genetic structure of European sheep breeds. Heredity. 2007;99:620–31.View ArticlePubMedGoogle Scholar
- Tapio M, Miceikiené I, Vikki J, Kantanen J. Comparison of microsatellite and blood protein diversity in sheep: inconsistencies in fragmented breeds. Mol Ecol. 2003;12:2045–56.View ArticlePubMedGoogle Scholar
- Peter C, Bruford M, Perez T, Dalamitra S, Hewitt G, Erhardt G, ECONOGENE Consortium. Genetic diversity and subdivision of 57 European and Middle-Eastern sheep breeds. Anim Genet. 2007;38:37–44.View ArticlePubMedGoogle Scholar
- Santucci F, Ibrahim KM, Bruzzone A, Hewit GM. Selection on MHC-linked microsatellite loci in sheep populations. Heredity. 2007;99:340–8.View ArticlePubMedGoogle Scholar
- Fernández J, Villanueva B, Pong-Wong R, Toro MA. Efficiency of the use of pedigree and molecular marker information in conservation programs. Genetics. 2005;170:1313–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Slate J, David P, Dodds KG, Veenvliet BA, Glass BC, Broad TE, McEwan JC. Understanding the relationship between the inbreeding coefficient and multilocus heterozygosity: theoretical expectations and empirical data. Heredity. 2004;93:255–65.View ArticlePubMedGoogle Scholar
- Bensch S, Andrén H, Hansson B, Pedersen HC, Sand H, Sejberg D, et al. Selection for heterozygosity gives hope to a wild population of inbred wolves. PLoS One. 2016; doi: 10.1371/journal.pone.000007.
- 25. Li MH, Strandén I, Tiirikka T, Sevón-Aimonen ML, Kantanen J. A comparison of approaches to estimate the inbreeding coefficient and pairwise relatedness using genomic and pedigree data in a sheep population. PLoS One. 2011; doi:10.1371/journal.pone.0026256.
- Wooliams JA, Gwaze GP, Meuwissen THE, Planchenault D, Renard JP, Thibier M, Wagner H, editors. Secondary guidelines for development of national farm animal genetic resources management plans: Management of small populations at risk. Food and Agriculture Organization of the United Nations. 1998.Google Scholar
- Bijma P. Long-term genetic contributions. Predictions of rates of inbreeding and genetic gain in selected populations, PhD Thesis. the Netherlands: Wageningen University; 2000.Google Scholar
- Franklin IR, Frankham R. How large must populations be to retain evolutionary potential? Animal Conserv. 1998;1:69–70.View ArticleGoogle Scholar