Tracking diversity and differentiation in six sheep breeds from the North Iberian Peninsula through DNA variation

IBD sharing patterns as intra-breed admixture indicators in small ruminants

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

Genetic Differentiation among Livestock Breeds-Values for Fst

(1) Background: The Fst statistic is widely used to characterize between-breed relationships. Fst = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding Fst values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed Fst were as follows, cattle: microsatellite 0.06-0.12, SNP 0.08-0.15; sheep: microsatellite 0.06-0.10, SNP 0.06-0.17; horses: microsatellite 0.04-0.11, SNP 0.08-0.12; goats: microsatellite 0.04-0.14, SNP 0.08-0.16; pigs: microsatellite 0.06-0.27, SNP 0.15-0.22; chickens: microsatellite 0.05-0.28, SNP 0.08-0.26. (4) Conclusions: (1) Large amounts of Fst data are available for a substantial proportion of the world's livestock breeds, (2) the value for between-breed Fst of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated Fst data may have value for interdisciplinary research.

Genetic Characterization of a Sheep Population in Oaxaca, Mexico: The Chocholteca Creole

Creole sheep in México have undergone crossbreeding, provoking the loss of genetic variability. The objective of the present study is to determine the intra-racial genetic diversity, the genetic relationship with other genotypes, and the populational substructure of the Oaxacan Creole sheep. Twenty-nine blood samples were obtained of Creole sheep of the Oaxaca Mixteca region in México. A genetic analysis was made with 41 microsatellites recommended for studies of genetic diversity in sheep. An analysis was made of genetic diversity, populational structure, and genetic distance with 27 other sheep populations. The study found 205 alleles with a range of 2 to 9 by locus and an effective number of 3.33. The intra-racial analysis showed a moderate genetic diversity with values of expected heterozygosity of 0.686 and observed of 0.756, a mean polymorphic information content of 0.609, and a mean coefficient of consanguinity of -0.002. In interracial genetic diversity for the coefficients of consanguinity, the values were FIS = 0.0774, FIT = 0.16993, and FST = 0.10028, showing an elevated genetic distance with other creole breeds, but close to Argentine Creole, to another Creole of México and the Spanish Merino. Its genetic structure showed that it does not have any populational subdivision nor mixes with the others analyzed. It is concluded that it is a distinct and isolated population and is proposed as the creole breed "Chocholteca" for its conservation.

Genomic selection signatures in sheep from the Western Pyrenees

BACKGROUND

The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades.

RESULTS

Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index F_st. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program.

CONCLUSIONS

Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.

Analysis of population structure and genetic variability in Iranian buffaloes (Bubalus bubalis) using pedigree information

The objective of this study was to use pedigree analysis to evaluate the population structure, genetic variability and inbreeding in Iranian buffaloes. The analysis was based on the pedigree information of 42 285 buffaloes born from 549 sires and 6376 dams within 1697 herds. Pedigree information used in this study was collected during 1976 to 2012 by the Animal Breeding Centre of Iran. The CFC program was applied to calculate pedigree statistics and genetic structure analysis of the Iranian buffaloes. Also, the INBUPGF90 program was used for calculating regular inbreeding coefficients for individuals in the pedigree. The analysis of pedigree indicated that inbreeding coefficient ranged from 0% to 31% with an average of 3.42% and the trend of inbreeding was significantly positive over the years (P < 0.0001). Average coancestry was increased in recent years and overall generation interval was 6.62 years in Iranian buffaloes. Founder genome equivalent, founder equivalent, effective number of founders and effective number of non-founders were increased from 1976 to 2002, but their values decreased from 2002 onwards. A designed mating system to avoid inbreeding may be applied to this population of buffalo to maintain genetic diversity.

Population genetic structure analysis and effect of inbreeding on body weights at different ages in Iranian Mehraban sheep

The objective of this study was to describe the population structure and inbreeding, and to quantify their effects on weights at different ages of Mehraban sheep in Iran. The analysis was based on the pedigree information of 26990 animals and 10278 body weight records from birth to yearling age. Data and pedigree information were collected during 1994 to 2011 by the breeding station of Mehraban sheep. The population structure was analyzed using the CFC program. Inbreeding of all animals was calculated by INBUPGF90 program. All animals were grouped into three classes according to their inbreeding coefficients: the first class included non-inbred animals (F=0); and the second and third classes included inbred animals (0<F<0.05 and F≥0.05, respectively). The average inbreeding in Mehraban sheep was 1.69%. Founder equivalent (f_e) values were estimated to be 4244, 3116 and 2965 during 1994-1999, 2000-2005 and 2006-2011, respectively. The effective population sizes (N_e) were 363, 5080 and 5740 during 1994-1999, 2000-2005 and 2006-2011, respectively. Generation interval was 2.15 years for this breed of sheep. Regression coefficients of birth weight, 3-month weight, 6-month weight and yearling weight on lamb inbreeding were estimated to be -6.34±0.69, -14.68±5.33, 48.00±9.43 and 98.65±15.65, respectively. Both positive and negative inbreeding effects were found in the current study. The utilization of a program for designed mating system, in the present flock, could be a suitable approach to keep the level of inbreeding under control.

Conservation genetics of cattle, sheep, and goats

Cattle, sheep and goats were domesticated about 10,000 years ago, spread out of the domestication centers in Europe, Asia, and Africa during the next few thousands years, and gave many populations locally adapted. After a very long period of soft selection, the situation changed dramatically 200 years ago with the emergence of the breed concept. The selection pressure strongly increased, and the reproduction among breeds was seriously reduced, leading to the fragmentation of the initial gene pool. More recently, the selection pressure was increased again via the use of artificial insemination, leading to a few industrial breeds with very high performances, but with low effective population sizes. Beside this performance improvement of industrial breeds, genetic resources are being lost, because of the replacement of traditional breeds by high performance industrial breeds at the worldwide level, and because of the loss of genetic diversity in these industrial breeds. Many breeds are already extinct, and genetic resources in cattle, sheep, and goats are thus highly endangered, particularly in developed countries. The recent development of next generation sequencing technologies opens new avenues for properly characterizing the genetic resources, not only in the very diverse domestic breeds, but also in their wild relatives. Based on sound genetic characterization, urgent conservation measures must be taken to avoid an irremediable loss of farm animal genetic resources, integrating economical, sociological, and political parameters.

Genetic diversity of selected genes that are potentially economically important in feral sheep of New Zealand

Background

Feral sheep are considered to be a source of genetic variation that has been lost from their domestic counterparts through selection.

Methods

This study investigates variation in the genes KRTAP1-1, KRT33, ADRB3 and DQA2 in Merino-like feral sheep populations from New Zealand and its offshore islands. These genes have previously been shown to influence wool, lamb survival and animal health.

Results

All the genes were polymorphic, but no new allele was identified in the feral populations. In some of these populations, allele frequencies differed from those observed in commercial Merino sheep and other breeds found in New Zealand. Heterozygosity levels were comparable to those observed in other studies on feral sheep. Our results suggest that some of the feral populations may have been either inbred or outbred over the duration of their apparent isolation.

Conclusion

The variation described here allows us to draw some conclusions about the likely genetic origin of the populations and selective pressures that may have acted upon them, but they do not appear to be a source of new genetic material, at least for these four genes.

Microsatellite based ovine parentage testing to identify the source responsible for the killing of an endangered species

In this study, we present an animal forensic genetics case in which we applied ovine microsatellite based parentage testing to the resolution of an act against protected wildlife. The aim was to identify the origin of the animal responsible for the death of an Egyptian vulture (Neophron percnopterus) that had been poisoned through consumption of a Latxa Blond Face breed lamb. Thus, we carried out a 22 microsatellite based parentage test in order to identify the parents of the lamb in the flock which grazes in the same place where the vulture's remains were found. Multiple parentage analysis revealed two possible parents, one ewe and one ram, with a combined paternity/maternity index (PI) higher than 9.09E9 and a likelihood (W) value of 99.9999%, assuming 50% probability a priori. This result confirmed the flock of origin of the poisoned lamb and shows that the microsatellite panel described herein is a potentially useful tool for the resolution of animal forensic cases.

Genetic variation and population structure of Italian native sheep breeds undergoing in situ conservation

The genetic variability and presence of population substructures in 4 native Northern Italian sheep breeds, Alpagota, Brogna, Foza, and Lamon, undergoing in situ conservation, and 1 widespread Italian breed, Bergamasca, were studied by investigating 19 microsatellite markers. The breeds showed considerable genetic variability in terms of number of alleles and heterozygosity, with the exception of Alpagota, which was the least variable (0.607). Nevertheless, a significant deficit of heterozygotes was observed in each breed due to rather increased levels of inbreeding or to the presence of population substructures, probably caused by increased genetic variation in the founder populations. The analyses evidenced clear genetic differentiation (F(ST) = 0.085), reduced levels of admixture, and presence of private alleles among the breeds, confirming their genetic uniqueness. In particular, according to Reynolds genetic distances, Alpagota was the most differentiated, perhaps because it had been bred mostly in a rather isolated area. Loss of any of the investigated breeds would result in a loss of genetic diversity ranging from 0.5 to 1.6% of the total observed gene diversity. Results supported the decision to safeguard these breeds as important reservoirs of genetic diversity and suggested breeding and mating practices to maintain variability and to overcome within-breed substructures.

Population structure, genetic variation and morphological diversity in indigenous sheep of Ethiopia

We investigated genetic and morphological diversity and population structure of 14 traditional sheep populations originating from four ecological zones in Ethiopia (sub-alpine, wet highland, sub-humid lowland and arid lowland). All animals (n = 672) were genotyped for 17 microsatellite markers and scored for 12 morphological characters. The sheep were initially classified as fat-tailed (11 populations), thin-tailed (one population) and fat-rumped sheep (two populations). These classifications are thought to correspond to three consecutive introduction events of sheep from the Near-East into East Africa. For the 14 populations, allelic richness ranged from 5.87 to 7.51 and expected heterozygosity (H(E)) from 0.66 to 0.75. Genetic differentiations (F(ST) values) between all pairs of populations, except between sub-alpine populations, were significantly different from zero (P < 0.001). Cluster analysis of morphological characters and a dendrogram constructed from genetic distances were broadly consistent with the classification into fat-tailed, thin-tailed and fat-rumped sheep. Bayesian cluster analysis using microsatellite markers indicated that there has been further genetic differentiation after the initial introduction of sheep into Ethiopia. Investigation of factors associated with genetic variation showed that an isolation-by-distance model, independently of other factors, explained most of the observed genetic variation. We also obtained a strong indication of adaptive divergence in morphological characters, patterns of morphological variation being highly associated with ecology even when the effect of neutral genetic divergence (F(ST)) was parcelled out in partial Mantel tests. Using a combination of F(ST) values, Bayesian clustering analysis and morphological divergence, we propose a classification of Ethiopian sheep into six breed groups and nine breeds.

Genetic diversity and subdivision of 57 European and Middle-Eastern sheep breeds

The population structure and genetic diversity of 57 European and Middle Eastern marginal and cosmopolitan sheep breeds from 15 countries were analysed by typing 31 microsatellite markers. Mean unbiased expected heterozygosities ranged from 0.63 in British Exmoor Horn to 0.77 in Albanian Ruda. South-eastern European and Middle-Eastern sheep breeds were significantly more variable than northwestern and western European breeds. An overall heterozygote deficiency (f) across all loci was observed (P < 0.001), while genetic differentiation (theta) was 5.7%. Principal component analysis and Bayesian model-based clustering indicate a south-east to north-west cline, but also revealed distinct groups of Middle-Eastern fat-tailed sheep, south-eastern European sheep and north-western/western European sheep. Within the last group, two less-distinct clusters comprised the Merino-type and Alpine breeds respectively. The incomplete demarcations of most clusters probably reflects cross-breeding and/or upgrading.

Genetic Signatures of a Mediterranean Influence in Iberian Peninsula Sheep Husbandry

Highly adaptable and versatile populations of domestic sheep, the result of millennia of intense husbandry, are found in almost every corner of the world. Here we describe a genetic survey of sheep from the western fringe of its European distribution. We studied the mitochondrial DNA control region sequences from 161 individuals belonging to 7 Portuguese sheep breeds. Our study revealed a high level of genetic diversity, with an average breed haplotype diversity of 0.983, substantially above that observed in central European breeds, as well as the presence of maternal lineages until now only found in the Middle East and Asia. A broad north-south pattern describes the most important trend in the Portuguese sheep population with a southern population clearly distinct from most other breeds. A recurrent influx of new genetic diversity, probably via the Mediterranean Sea, may explain these patterns and appears to corroborate the importance of this maritime route in the history of both mankind and livestock. Zooarchaeological studies of sheep bones from southern Portugal indicate a marked size increase during the Moslem period that may reflect an improvement of this animal--perhaps part of the well known "Arab agricultural revolution" in Andalusia. This could have been a time when the gene pool of Iberian sheep was substantially enriched and may help to explain the history of modern sheep breeds in this peninsula.