Testing the usefulness of the molecular coancestry information to assess genetic relationships in livestock using a set of Spanish sheep breeds1

Historical Westward Migration Phases of Ovis aries Inferred from the Population Structure and the Phylogeography of Occidental Mediterranean Native Sheep Breeds

In this study, the genetic relationship and the population structure of western Mediterranean basin native sheep breeds are investigated, analyzing Maghrebian, Central Italian, and Venetian sheep with a highly informative microsatellite markers panel. The phylogeographical analysis, between breeds’ differentiation level (Wright’s fixation index), gene flow, ancestral relatedness measured by molecular coancestry, genetic distances, divergence times estimates and structure analyses, were revealed based on the assessment of 975 genotyped animals. The results unveiled the past introduction and migration history of sheep in the occidental Mediterranean basin since the early Neolithic. Our findings provided a scenario of three westward sheep migration phases fitting properly to the westward Neolithic expansion argued by zooarcheological, historical and human genetic studies.

Genetic Differentiation among Livestock Breeds—Values for Fst

Simple Summary

The degree of relationship among livestock breeds can be quantified by the F_st statistic, which measures the extent of genetic differentiation between them. An F_st value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, F_st values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that F_st > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of F_st data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds.

Abstract

(1) Background: The F_st statistic is widely used to characterize between-breed relationships. F_st = 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 F_st 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 F_st 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 F_st data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed F_st of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated F_st data may have value for interdisciplinary research.

Inbreeding, effective population size, and coancestry in the Latxa dairy sheep breed

Traditionally, breeding programs have estimated and managed inbreeding based on pedigree information. The availability of genomic marker panels has made possible new alternatives to achieve more precise estimates, for example in case of missing pedigree. The objective of the present study was to assess and compare, different estimation methods (pedigree-based methodologies, single SNP-based approach (homozygosity) and runs of homozygosity-based method) to analyze the evolution of genetic diversity measured as inbreeding or as coancestry of 3 selected populations of Latxa dairy sheep (Latxa Cara Rubia and Latxa Cara Negra from Euskadi and Navarre). Genomic data came from 972 artificial insemination rams genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA) whose genealogy consisted of 4,484 animals. Inbreeding estimates based on molecular data were more similar between them than compared with those based on pedigree information. However, the SNP-based approach estimations of effective population size differed more, reflecting the sensitivity of effective population size to small changes in the evolution of inbreeding. The 2 Latxa Cara Negra populations showed increases of inbreeding rates with time and effective population sizes between 64 and 103 animals, depending on breed and methodology used. The Latxa Cara Rubia population did not show an increase in inbreeding rate, mainly due to semen importation from the related French population of Manech Tête Rousse. The effective size estimates based on coancestry increase show a higher variability and they are more sensitive to the source of information and the data structure considered. Realized effective population size based on individual increase in inbreeding were in agreement with the previous estimates. Coancestry evolution analysis based on DNA information showed an increase on coancestry during the last 10 yr in all breeds, as a consequence of the selection process. Moreover, the increase on coancestry between Latxa Cara Rubia and Manech Tête Rousse was more noticeable between than within each of those breeds.

Conservation genetics in Chinese sheep: diversity of fourteen indigenous sheep (Ovis aries) using microsatellite markers

The domestic sheep (Ovis aries) has been an economically and culturally important farm animal species since its domestication around the world. A wide array of sheep breeds with abundant phenotypic diversity exists including domestication and selection as well as the indigenous breeds may harbor specific features as a result of adaptation to their environment. The objective of this study was to investigate the population structure of indigenous sheep in a large geographic location of the Chinese mainland. Six microsatellites were genotyped for 611 individuals from 14 populations. The mean number of alleles (±SD) ranged from 7.00 ± 3.69 in Gangba sheep to 10.50 ± 4.23 in Tibetan sheep. The observed heterozygote frequency (±SD) within a population ranged from 0.58 ± 0.03 in Gangba sheep to 0.71 ± 0.03 in Zazakh sheep and Minxian black fur sheep. In addition, there was a low pairwise difference among the Minxian black fur sheep, Mongolian sheep, Gansu alpine merino, and Lanzhou fat-tailed sheep. Bayesian analysis with the program STRUCTURE showed support for 3 clusters, revealing a vague genetic clustering pattern with geographic location. The results of the current study inferred high genetic diversity within these native sheep in the Chinese mainland.

Evaluation of genetic variability within PrP genotyped sheep of endangered Italian Altamurana breed

In the last few decades, there has been increased awareness of preservation and exploitation of endangered breeds for the maintenance of biodiversity, as well as the concern for diseases in sheep breeding. This study was carried out in native endangered Altamurana dairy sheep breed from Southern Italy. The Altamurana breed was considered as two populations (Alt-Cav and Alt-Cra-Zoe), based on presumed cross-breed and remarkable differences in the PrP genotypes frequencies. The genetic diversity between the two Altamurana populations (Alt-Cav and Alt-Cra-Zoe) was evaluated in comparison to three Italian dairy breeds through fourteen microsatellite markers. Both measures of genetic distance and the population structure analysis highlighted that the Alt-Cav and Alt-Cra-Zoe sheep have a particular genetic background. The estimated fixation index (FST) and the genetic Nei's distances among the populations showed a higher level of genetic differentiation for Alt-Cav than Alt-Cra-Zoe. The Bayesian clustering analysis implemented in the STRUCTURE software showed clear and distinct clusters for the two Altamurana populations, confirming the hypothesis of Alt-Cav as a genetic group well differentiated from Alt-Cra-Zoe. Alt-Cav likely can be considered as belonging to the original strain of the Altamurana breed. This findings may be used to assist the programme for conservation and selection of scrapie resistance genotypes in endangered Altamurana sheep breed.

Recent advances in understanding the genetic resources of sheep breeds locally-adapted to the UK uplands: opportunities they offer for sustainable productivity

Locally adapted breeds of livestock are of considerable interest since they represent potential reservoirs of adaptive fitness traits that may contribute to the future of sustainable productivity in a changing climate. Recent research, involving three hill sheep breeds geographically concentrated in the northern uplands of the UK has revealed the extent of their genetic diversity from one another and from other breeds. Results from the use of SNPs, microsatellites, and retrovirus insertions are reviewed in the context of related studies on sheep breeds world-wide to highlight opportunities offered by the genetic resources of locally adapted hill breeds. One opportunity concerns reduced susceptibility to Maedi Visna, a lentivirus with massive impacts on sheep health and productivity globally. In contrast to many mainstream breeds used in farming, each of the hill breeds analyzed are likely to be far less susceptible to the disease threat. A different opportunity, relating specifically to the Herdwick breed, is the extent to which the genome of the breed has retained primitive features, no longer present in other mainland breeds of sheep in the UK and offering a new route for discovering unique genetic traits of use to agriculture.

Assessing introgression of Sahelian zebu genes into native Bos taurus breeds in Burkina Faso

A total of 350 samples were analyzed to estimate zebu gene proportions into two different taurine cattle breeds of Burkina Faso (Lobi and N'Dama) using 38 microsatellites and various statistical methodologies. West African and East African zebu samples were sequentially used as reference parental populations. Furthermore, N'Dama cattle from Congo, the composite South African Bonsmara cattle breed and a pool of European cattle were used successively as second parental populations. Independently of the methodology applied: (a) the use of West African zebu samples gave higher admixture coefficients than the East African zebu; (b) the higher zebu proportions were estimated when the European cattle was used as parental population 2; and (c) the use of the N'Dama population from Congo as parental population 2 gave the more consistent zebu proportion estimates for both the Lobi and the N'Dama breeds. In any case, the zebu admixture proportions estimated were not negligible and were always higher in the N'Dama cattle than in the Lobi cattle of Burkina Faso. This suggested that the introgression of Sahelian zebu genes into the taurine cattle of Southern West Africa can follow a complex pattern that can depend on local agro-ecological features. The current research pointed out that the estimation of admixture coefficients is highly dependent on both the assumptions underlying the methodologies applied and the selection of parental populations. Our analyses suggest that either too high or nil genetic identity between the parental and the expectedly derived populations must be avoided.

Genetic distinctiveness of the Herdwick sheep breed and two other locally adapted hill breeds of the UK

There is considerable interest in locally adapted breeds of livestock as reservoirs of genetic diversity that may provide important fitness traits for future use in agriculture. In marginal areas, these animals contribute to food security and extract value from land unsuitable for other systems of farming. In England, close to 50% of the national sheep flock is farmed on grassland designated as disadvantaged areas for agricultural production. Many of these areas are in the uplands, where some native breeds of sheep continue to be commercially farmed only in highly localised geographical regions to which they are adapted. This study focuses on three of these breeds, selected for their adaptation to near identical environments and their geographical concentration in regions close to one another. Our objective has been to use retrotyping, microsatellites and single nucleotide polymorphisms to explore the origins of the breeds and whether, despite their similar adaptations and proximity, they are genetically distinctive. We find the three breeds each have a surprisingly different pattern of retrovirus insertions into their genomes compared with one another and with other UK breeds. Uniquely, there is a high incidence of the R0 retrotype in the Herdwick population, characteristic of a primitive genome found previously in very few breeds worldwide and none in the UK mainland. The Herdwick and Rough Fells carry two rare retroviral insertion events, common only in Texels, suggesting sheep populations in the northern uplands have a historical association with the original pin-tail sheep of Texel Island. Microsatellite data and analyses of SNPs associated with RXFP2 (horn traits) and PRLR (reproductive performance traits) also distinguished the three breeds. Significantly, an SNP linked to TMEM154, a locus controlling susceptibility to infection by Maedi-Visna, indicated that all three native hill breeds have a lower than average risk of infection to the lentivirus.

Ascertaining gene flow patterns in livestock populations of developing countries: a case study in Burkina Faso goat

Background

Introgression of Sahel livestock genes southwards in West Africa may be favoured by human activity and the increase of the duration of the dry seasons since the 1970’s. The aim of this study is to assess the gene flow patterns in Burkina Faso goat and to ascertain the most likely factors influencing geographic patterns of genetic variation in the Burkina Faso goat population.

Results

A total of 520 goat were sampled in 23 different locations of Burkina Faso and genotyped for a set of 19 microsatellites. Data deposited in the Dryad repository: http://dx.doi.org/10.5061/dryad.41h46j37. Although overall differentiation is poor (F_ST = 0.067 ± 0.003), the goat population of Burkina Faso is far from being homogeneous. Barrier analysis pointed out the existence of: a) genetic discontinuities in the Central and Southeast Burkina Faso; and b) genetic differences within the goat sampled in the Sahel or the Sudan areas of Burkina Faso. Principal component analysis and admixture proportion scores were computed for each population sampled and used to construct interpolation maps. Furthermore, Population Graph analysis revealed that the Sahel and the Sudan environmental areas of Burkina Faso were connected through a significant number of extended edges, which would be consistent with the hypothesis of long-distance dispersal. Genetic variation of Burkina Faso goat followed a geographic-related pattern. This pattern of variation is likely to be related to the presence of vectors of African animal trypanosomosis. Partial Mantel test identified the present Northern limit of trypanosome vectors as the most significant landscape boundary influencing the genetic variability of Burkina Faso goat (p = 0.008). The contribution of Sahel goat genes to the goat populations in the Northern and Eastern parts of the Sudan-Sahel area of Burkina Faso was substantial. The presence of perennial streams explains the existence of trypanosome vectors. The South half of the Nakambé river (Southern Ouagadougou) and the Mouhoun river loop determined, respectively, the Eastern and Northern limits for the expansion of Sahelian goat genes. Furthermore, results from partial Mantel test suggest that the introgression of Sahelian goat genes into Djallonké goat using human-influenced genetic corridors has a limited influence when compared to the biological boundary defined by the northern limits for the distribution of the tsetse fly. However, the genetic differences found between the goat sampled in Bobo Dioulasso and the other populations located in the Sudan area of Burkina Faso may be explained by the broad goat trade favoured by the main road of the country.

Conclusions

The current analysis clearly suggests that genetic variation in Burkina Faso goat: a) follows a North to South clinal; and b) is affected by the distribution of the tsetse fly that imposes a limit to the Sahelian goat expansion due to their trypanosusceptibility. Here we show how extensive surveys on livestock populations can be useful to indirectly assess the consequences of climate change and human action in developing countries.

Assessing priorities for conservation in Tuscan cattle breeds using microsatellites

Preservation of rare genetic stocks requires assessment of within-population genetic diversity and between-population differentiation to make inferences on their degree of uniqueness. A total of 194 Tuscan cattle (44 Calvana, 35 Chianina, 25 Garfagnina, 31 Maremmana, 31 Mucca Pisana and 28 Pontremolese) individuals were genotyped for 34 microsatellite markers. Moreover, 56 samples belonging to Argentinean Creole and Asturiana de la Montaña cattle breeds were used as an outgroup. Genetic diversity was quantified in terms of molecular coancestry and allelic richness. STRUCTURE analyses showed that the Tuscan breeds have well-differentiated genetic backgrounds, except for the Calvana and Chianina breeds, which share the same genetic ancestry. The between-breed Nei's minimum distance (Dm) matrices showed that the pair Calvana-Chianina was less differentiated (0.049 ± 0.006). The endangered Tuscan breeds (Calvana, Garfagnina, Mucca Pisana and Pontremolese) made null or negative contributions to diversity, except for the Mucca Pisana contribution to allelic richness (CT = 1.8%). The Calvana breed made null or negative within-breed contributions (W = 0.0%; CW = -0.4%). The Garfagnina and Pontremolese breeds made positive contributions to between-breed diversity but negative and high within-breed contributions, thus suggesting population bottleneck with allelic losses and increase of homozygosity in the population. Exclusion of the four endangered Tuscan cattle breeds did not result in losses in genetic diversity (T = -0.7%; CT = -1.2%), whereas exclusion of the non-endangered breeds (Chianina and Maremmana) did (T = 2.1%; CT = 3.9%); the simple exclusion of the Calvana breed from the former group led to losses in genetic diversity (T = 0.47%; CT = 2.34%), indicating a diverse significance for this breed. We showed how quantifying both within-population diversity and between-population differentiation in terms of allelic frequencies and allelic richness provides different and complementary information on the genetic backgrounds assessed and may help to implement priorities and strategies for conservation in livestock.

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.

Genetic diversity of Swiss sheep breeds in the focus of conservation research

There is constant pressure to improve evaluation of animal genetic resources in order to prevent their erosion. Maintaining the integrity of livestock species as well as their genetic diversity is of paramount interest for long-term agricultural policies. One major use of DNA techniques in conservation is to reveal genetic diversity within and between populations. Forty-one microsatellites were analysed to assess genetic diversity in nine Swiss sheep breeds and to measure the loss of the overall diversity when one breed would become extinct. The expected heterozygosities varied from 0.65 to 0.74 and 10.8% of the total genetic diversity can be explained by the variation among breeds. Based on the proportion of shared alleles, each of the nine breeds were clearly defined in their own cluster in the neighbour-joining tree describing the relationships among the breeds. Bayesian clustering methods assign individuals to groups based on their genetic similarity and infer the number of populations. In STRUCTURE, this approach pooled the Valais Blacknose and the Valais Red. With BAPS method the two Valais sheep breeds could be separated. Caballero & Toro approach (2002) was used to calculate the loss or gain of genetic diversity when each of the breeds would be removed from the set. The changes in diversity based on between-breed variation ranged from -12.2% (Valais Blacknose) to 0% (Swiss Black Brown Mountain and Mirror Sheep); based on within-breed diversity the removal of a breed could also produce an increase in diversity (-0.6% to + 0.6%). Allelic richness ranged from 4.9 (Valais Red) to 6.7 (Brown Headed Meat sheep and Red Engadine Sheep). Breed conservation decisions cannot be limited to genetic diversity alone. In Switzerland, conservation goals are embedded in the desire to carry the cultural legacy over to future generations.

The genealogy and genetic viability of reintroduced Yellowstone grey wolves

The recovery of the grey wolf in Yellowstone National Park is an outstanding example of a successful reintroduction. A general question concerning reintroduction is the degree to which genetic variation has been preserved and the specific behavioural mechanisms that enhance the preservation of genetic diversity and reduce inbreeding. We have analysed 200 Yellowstone wolves, including all 31 founders, for variation in 26 microsatellite loci over the 10-year reintroduction period (1995-2004). The population maintained high levels of variation (1995 H(0) = 0.69; 2004 H(0) = 0.73) with low levels of inbreeding (1995 F(IS) = -0.063; 2004 F(IS) = -0.051) and throughout, the population expanded rapidly (N(1995) = 21; N(2004) = 169). Pedigree-based effective population size ratios did not vary appreciably over the duration of population expansion (1995 N(e)/N(g) = 0.29; 2000 N(e)/N(g) = 0.26; 2004 N(e)/N(g) = 0.33). We estimated kinship and found only two of 30 natural breeding pairs showed evidence of being related (average r = -0.026, SE = 0.03). We reconstructed the genealogy of 200 wolves based on genetic and field data and discovered that they avoid inbreeding through a wide variety of behavioural mechanisms including absolute avoidance of breeding with related pack members, male-biased dispersal to packs where they breed with nonrelatives, and female-biased subordinate breeding. We documented a greater diversity of such population assembly patterns in Yellowstone than previously observed in any other natural wolf population. Inbreeding avoidance is nearly absolute despite the high probability of within-pack inbreeding opportunities and extensive interpack kinship ties between adjacent packs. Simulations showed that the Yellowstone population has levels of genetic variation similar to that of a population managed for high variation and low inbreeding, and greater than that expected for random breeding within packs or across the entire breeding pool. Although short-term losses in variation seem minimal, future projections of the population at carrying capacity suggest significant inbreeding depression will occur without connectivity and migratory exchange with other populations.

Molecular analysis using mitochondrial DNA and microsatellites to infer the formation process of Japanese native horse populations

To assess the genetic diversity of Japanese native horse populations, we examined seven such populations using mitochondrial DNA (mtDNA) and microsatellite analyses. Four reference populations of Mongolian horses and European breeds were employed as other equids. In the mtDNA analysis, the control region (D-loop) of 411 bp was sequenced, and 12 haplotypes with 33 variable sites were identified in the Japanese native horses. The phylogenetic tree constructed by haplogrouping and using worldwide geographic references indicated that the haplotypes of the Japanese native horses were derived from six equid clusters. Compared with the foreign populations, the Japanese native populations showed lower within-population diversity and higher between-population differentiation. Microsatellite analysis, using 27 markers, found an average number of alleles per locus of 9.6 in 318 native and foreign horses. In most native populations, the within-population diversity was lower than that observed in foreign populations. The genetic distance matrix based on allelic frequency indicated that several native populations had notably high between-population differentiation. The molecular co-ancestry-based genetic distance matrix revealed that the European populations were differentiated from the Japanese and Mongolian populations, and no clear groups could be identified among the Japanese native horse populations. The genetic distance matrices had few correlations with the geographic distribution of the Japanese native populations. Based on the results of both mtDNA and microsatellite analyses, it could be speculated that each native population was formed by the founder populations derived from Mongolian horses. The genetic construction of each population appears to have been derived from independent breeding in each local area since the time of population fission, and this was accompanied by drastic genetic drift in recent times. This information will help to elucidate the ancestry of Japanese native horses.

Individual-based assessment of population structure and admixture in Austrian, Croatian and German draught horses

All over Europe, the number of draught horses has decreased drastically during the last 50 years. As a prerequisite for efficient management decisions, we analysed the conservation status in Austrian (Noriker Carinthia - NC, Noriker Salzburg - NS), Croatian (Croatian Coldblood - C, Posavina horse - P) and German (Altmaerkisch Coldblood - A, Black Forest horse - BF, Mecklenburg Coldblood - M, Rhenish German Draught horse - R, Saxon Thuringa Coldblood - ST, Schleswig Draught horse - Sch, South German Coldblood - SG) draught horses (434) using multilocus genotypic information from 30 (effectively 27) microsatellite loci. Populations located in areas with less intensive agricultural production (C, NC, NS, P and SG) had greater diversity within the population and estimated effective population size than A, BF, Sch, M, R and ST populations. The PCA plots revealed that populations form five separate groups. The 'Noriker' group (NC, NS and SG) and the 'Rhenish' group (A, M, R and ST) were the most distinctive (pairwise F(ST) values ranged from 0.078 to 0.094). The 'Croatian' group (C and P) was in the centre, while the BF and Sch populations formed two out-groups. A posterior Bayesian analysis detected further differentiation, mainly caused by political and geographical factors. Thus, it was possible to separate the South German Coldblood from the Austrian Noriker population where no subpopulation structure was detected. The admixture analysis revealed imprecise classification between C and P populations. A small but notable separation of R from A, M and ST populations was detected, while Sch and BF populations remained as out-groups. The information obtained should aid in making efficient conservation decisions.