Genomic diversity of the locally developed Latvian Darkheaded sheep breed

The Latvian Darkheaded is the only locally developed sheep breed. The breed was formed at the beginning of the 20th century by crossing local coarse-wooled sheep with the British Shropshire and Oxfordshire breeds. The breed was later improved by adding Ile-de-France, Texel, German blackheads, and Finnsheep to achieve higher prolificacy and better meat quality. Previous studies have reported the Latvian Darkheaded sheep to be closely related to Estonian and Lithuanian Blackface breeds, according to microsatellite data. To expand our knowledge of the genetic resources of the Latvian Darkheaded breed, we conducted a whole-genome resequencing analysis on 40 native sheep. The investigation showed that local sheep harbor genetic diversity levels similar to those observed among other improved breeds of European origin, including Charollais and Suffolk. Genome-wide nucleotide diversity (π) in Latvian Darkheaded sheep was 3.91 × 10-3, whereas the average observed heterozygosity among the 40 animals was 0.267 and 0.438 within the subsample of unrelated individuals. The Ne has rapidly decreased to 200 ten generations ago with a recent drop to Ne 73 four generations ago. However, inbreeding levels based on runs of homozygosity were, on average, low, with FROH ranging between 0.016 and 0.059. The analysis of the genomic composition of the breed confirmed shared ancestry with sheep of British origin, reflecting the history of the breed. Nevertheless, Latvian Darkheaded sheep were genetically separable. The contemporary Latvian Darkheaded sheep population is genetically diverse with a low inbreeding rate. However, further development of breed management programs is necessary to prevent an increase in inbreeding, loss of genetic diversity, and depletion of breed-specific genetic resources, ensuring the preservation of the native Latvian Darkheaded sheep.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

The genome-wide relationships of the critically endangered Quadricorna sheep in the Mediterranean region

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Nowadays most of the sheep European breeds derive from the second wave which is supposed to have largely replaced oldest genetic signatures, with the exception of several primitive breeds confined on the very edge of Northern Europe. Despite this, some populations also in the Mediterranean region are characterised by the presence of phenotypic traits considered ancestral such as the policeraty, large horns in the ram, short tail, and a moulting fleece. Italy is home of a large number of local breeds, albeit some are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterised by several traits considered as ancestral. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 41 breeds from the Mediterranean region and Middle-East, with the specific aim to reconstruct its origin. After retaining 32,862 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. The Quadricorna breed exhibit genetic diversity indices comparable with those of most of the other analysed breeds, however, the two populations showed opposing patterns of genetic diversity suggesting different levels of genomic inbreeding and drift (FIS and FROH). In general, all the performed genome-wide analyses returned complementary results, indicating a westward longitudinal cline compatible with human migrations from the Middle-East and several additional genetic footprints which might mirror more recent historical events. Interestingly, among the Italian breeds, the original Quadricorna (QUAD_SA) first separated showing its own ancestral component. In addition, the admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of this population. On the other hand, both the neighbor-joining tree and the treemix analysis seem to suggest a proximity of the Quadricorna populations to breeds of South-Eastern Mediterranean origin. Although our results do not support a robust link between the genetics of the first wave and the presence of primitive traits, the observed genetic uniqueness together with the inferred phylogeograpic reconstruction would suggest an ancient presence of the Quadricorna breed in the Italian Peninsula. Because of this singularity, urgent conservation actions are needed in order to keep the breed and all related cultural products alive.

Native Sheep Breeds in Poland—Importance and Outcomes of Genetic Resources Protection Programmes

Simple Summary

The aim of this paper was to present the current situation of native sheep breeding in Poland, in terms of the significance and the effects of genetic resources protection programmes. The sheep farming tradition in Poland is deeply rooted in culture, particularly in the mountain and foothill regions. Sheep are permanently linked to many areas of country, not only playing an indispensable part in shaping the landscape, but also providing unique products. The most effective method of conserving native breeds is their sustainable management. It is important to maintain and develop the promotion and certification of high-quality products from native breeds and to control the endangered status of local breeds, which must be monitored and updated on an ongoing basis.

Abstract

The sheep population of native breeds, despite their unique features and the ability to adapt to harsh environmental conditions, has significantly decreased in recent years. Due to the low profitability of breeding, many local breeds of sheep in Poland were exposed to the risk of extinction. Many years of crisis in sheep farming have exacerbated this situation. The aim of this paper was to present the current situation of native sheep breeding in Poland, in terms of significance and effects of genetic resources protection programmes. The conservation of genetic resources of sheep aims to maintain and increase the population size while striving to maintain the greatest possible genetic variability. There are 17 native breeds included in the Polish sheep genetic resources conservation programme. A positive element of the implementation of the conservation of genetic resources programme for sheep is the accompanying measures which are based on the use of the non-productive role of the species. Extensive sheep grazing, as a form of nature conservation, serves to preserve valuable natural landscapes and the culture of local communities associated with sheep farming. Production and promotion of quality products, especially using niche markets and short production chains, are essential to ensure the economic viability of farms. These activities must be accompanied by raising public awareness of indigenous breeds and their alternative use in environmental activities, as well as their role in preserving the cultural heritage of local communities, for example through mountain grazing and the production of traditional products.

SNP Genotyping Characterizes the Genome Composition of the New Baisary Fat-Tailed Sheep Breed

Lamb meat has become increasingly popular in several nations during the last few decades, especially in Kazakhstan. Due to the rising demand for lamb meat, our sheep breeders developed a new fat-tailed sheep and named the breed Baisary. Animals of the Baisary breed are characterized by a large physique, strong constitution, stretched body, deep and wide chest, medium or large-sized fat tail, long legs (height at the withers of adult rams 85–100 cm, sheep 75–90 cm), long lanceolate ears and strong hooves. Lambs of the Baisary breed surpass their peers of the original parent breeds by 15–20% in live weight at the weaning period. To characterize the genetic structure of Baisary sheep and compare it with the ancestral breeds, we genotyped 247 individuals from five sheep breeds with Ovine SNP50K. The estimated private allelic richness ranged from 0.0030 to 0.0047, with the minimum and maximum provided by the Gissar (Giss1) and Kazakh meat-wool breeds, respectively. The highest and lowest FIS values, meanwhile, were observed in the Afghan fat-tailed population and Baisary sheep, respectively. The calculated inbreeding coefficient showed that Edilbay and Baisary sheep have excess heterozygosity. According to principal components analysis, Baisary are close to Gissar populations, the Afghan fat-tailed breed and Edilbay sheep. These results were consistent with the Admixture and phylogenetic analysis. Overall, our results indicated that Baisary sheep differ genetically from their progenitors.

Genomic Population Structure of the Main Historical Genetic Lines of Spanish Merino Sheep

Simple Summary

Historical documentation shows that the Spanish Merino sheep was selected over many centuries due to the quality of wool, following which it was used to originate all other Merino breeds around the world, mainly by crossbreeding with local breeds. Today, the historical genetic lines that originated the Spanish Merino are still preserved in several closed herds in which they have been bred for nearly 200 years, maintaining their original genetic purity. Our study demonstrates, using a genomic approach, the exceptional genetic richness and variability of these lines, which are clearly differentiated from modern Merino breeds, and must therefore be protected to safeguard the large genetic pool they represent.

Abstract

According to historiographical documentation, the Romans first began to select Merino sheep in the Iberian Peninsula during the first century, with the aim of obtaining a breed appreciated for the quality of its wool. This process continued locally during the Middle Ages, when Spanish sheep were protected, and their export to foreign countries was banned. It was during the 16th century when individual Merino sheep were allowed to spread around the world to be used to improve the wool quality of local breeds. However, the wool crisis of the 1960s shifted the selection criteria of the Merino breed towards meat production at the expenses of wool. Consequently, individuals that display the genetic and phenotypic characteristics of those sheep originally bred in the kingdom of Spain in the Middle Ages are extremely difficult to find in commercial herds. In this study, we characterized the genetic basis of 403 individuals from the main historical Spanish Merino genetic lines (Granda, Hidalgo, Lopez-Montenegro, Maeso, Donoso and Egea), which were bred in isolation over the last 200 years, using a genomic approach based on genotyping data from the Axiom™ Ovine 50K SNP Genotyping Array. Our analysis included measuring population structure, genomic differentiation indexes, runs of homozygosity (ROH) patterns, and an analysis of molecular variance (AMOVA). The results showed large genetic differences between the historical lines, even though they belong to the same breed. In addition, ROH analysis showed differences due to increased inbreeding among the ancient generations compared with the modern Merino lines, confirming the breed’s ancestral and closed origin. However, our results also showed a high variability and richness within the Spanish historical Merino lines from a genetic viewpoint. This fact, together with their great ability to produce high-quality wool, suggests that ancestral Merino lines from Spain should be considered a valuable genetic population to be maintained as a resource for the improvement of wool-producing sheep breeds all around the world.

Exploring genetic diversity and population structure of Punjab goat breeds using Illumina 50 K SNP bead chip

Pakistan has 35 goat breeds. Moreover, the province of Punjab has highest goat population constituting 37% of country's total population with seven goat breeds including Beetal, Daira Deen Panah, Nachi, Barbari, Teddi, Pahari, and Pothwari. The diversity study of breeds warrants the documentation of breeds particularly using genome wide panel of markers, i.e., SNP chip. The objective of the current study was to fill this gap of information. Therefore, in current study we collected total of 879 unrelated goat blood samples along with data on body weight measurements; genomic DNA was extracted, and genotyping was carried out using 50 K SNP bead chip. Quality control measures were performed in Plink 1.07. Genetic diversity was observed among studied populations using heterozygosity and pairwise F_ST estimates, principal component analysis, admixture analysis in Plink software with visualization in Clumpak, and constructing phylogenetic tree in Mega 7 software. Moderate to high level of heterozygosity was observed among the studied populations. Coefficient of inbreeding varied from 0.0186 ± 0.0327 in Pahari to 0.183 ± 0.0715 in Barbari. Barbari and Daira Deen Panah had quite higher level of inbreeding coefficient as compared to all other breeds with value of 0.183 ± 0.0715 and 0.1378 ± 0.0741, respectively. PCA identified three steps of subdividing the seven goat breeds at various levels of K. All the seven breeds made independent clusters at various levels of PCA. Admixture analysis revealed the distinctness of Teddi and Barbari breeds. Genetic sub-structuring was observed in the admixture patterns of Beetal breed. Moreover, high level of genetic admixture was observed in Nachi, Pahari, Pothwari, and Daira Deen Panah breeds. Admixture results were further interpreted by calculating pairwise F_ST values. Our results provided first insights about genetic diversity of Pakistani goat breeds based on genomic data. To conclude, the enriched goat breed diversity in Pakistan could provide valuable genetic reservoir for national breeding schemes.

Runs of Homozygosity and NetView analyses provide new insight into the genome-wide diversity and admixture of three German cattle breeds

Angler (RVA) and Red-and-White dual-purpose (RDN) cattle were in the past decades crossed with influential Red Holstein (RH) sires. However, genome-wide diversity studies in these breeds are lacking. The objective of the present study was to elucidate the genome-wide diversity and population structure of the three German cattle breeds. Using 40,851 single nucleotide polymorphism markers scored in 337 individuals, runs of homozygosity (ROH) were analysed in each breed. Clustering and a high-resolution network visualisation analyses were performed on an extended dataset that included 11 additional (outgroup) breeds. Genetic diversity levels were high with observed heterozygosity above 0.35 in all three breeds. Only RVA had a recent past effective population size (Ne) estimate above 100 at 5 generations ago. ROH length distribution followed a similar pattern across breeds and the majority of ROH were found in the length class of >5 to 10 Mb. Estimates of average inbreeding calculated from ROH (FROH) were 0.021 (RVA), 0.045 (RDN) and 0.053 (RH). Moderate to high positive correlations were found between FROH and pedigree inbreeding (FPED) and between FROH and inbreeding derived from the excess of homozygosity (FHOM), while the intercept of the regression of FROH on FPED was above zero. The population structure analysis showed strong evidence of admixture between RVA and RH. Introgression of RDN with RH genes was minimally detected and for the first time, the study uncovered Norwegian Red Cattle ancestry in RVA. Highly heterogeneous genetic background was found for RVA and RH and as expected, the breeds of the extended dataset effectively differentiated mostly based on geographical origin, validating our findings. The results of this study confirm the impact of RH sires on RVA and RDN populations. Furthermore, a close monitoring is suggested to curb further reduction of Ne in the breeds.

Population structure and genetic diversity of 25 Russian sheep breeds based on whole-genome genotyping

Background

Russia has a diverse variety of native and locally developed sheep breeds with coarse, fine, and semi-fine wool, which inhabit different climate zones and landscapes that range from hot deserts to harsh northern areas. To date, no genome-wide information has been used to investigate the history and genetic characteristics of the extant local Russian sheep populations. To infer the population structure and genome-wide diversity of Russian sheep, 25 local breeds were genotyped with the OvineSNP50 BeadChip. Furthermore, to evaluate admixture contributions from foreign breeds in Russian sheep, a set of 58 worldwide breeds from publicly available genotypes was added to our data.

Results

We recorded similar observed heterozygosity (0.354–0.395) and allelic richness (1.890–1.955) levels across the analyzed breeds and they are comparable with those observed in the worldwide breeds. Recent effective population sizes estimated from linkage disequilibrium five generations ago ranged from 65 to 543. Multi-dimensional scaling, admixture, and neighbor-net analyses consistently identified a two-step subdivision of the Russian local sheep breeds. A first split clustered the Russian sheep populations according to their wool type (fine wool, semi-fine wool and coarse wool). The Dagestan Mountain and Baikal fine-fleeced breeds differ from the other Merino-derived local breeds. The semi-fine wool cluster combined a breed of Romanian origin, Tsigai, with its derivative Altai Mountain, the two Romney-introgressed breeds Kuibyshev and North Caucasian, and the Lincoln-introgressed Russian longhaired breed. The coarse-wool group comprised the Nordic short-tailed Romanov, the long-fat-tailed outlier Kuchugur and two clusters of fat-tailed sheep: the Caucasian Mountain breeds and the Buubei, Karakul, Edilbai, Kalmyk and Tuva breeds. The Russian fat-tailed breeds shared co-ancestry with sheep from China and Southwestern Asia (Iran).

Conclusions

In this study, we derived the genetic characteristics of the major Russian local sheep breeds, which are moderately diverse and have a strong population structure. Pooling our data with a worldwide genotyping set gave deeper insight into the history and origin of the Russian sheep populations.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0399-5) contains supplementary material, which is available to authorized users.