Assessing genetic diversity and differentiation in Portuguese coarse-wool sheep breeds with microsatellite markers

Genome-wide assessment of the population structure and genetic diversity of four Portuguese native sheep breeds

As the effects of global warming become increasingly complex and difficult to manage, the conservation and sustainable use of locally adapted sheep breeds are gaining ground. Portuguese native sheep breeds are important reservoirs of genetic diversity, highly adapted to harsh environments and reared in low input production systems. Genomic data that would describe the breeds in detail and accelerate the selection of more resilient animals to be able to cope with climatic challenges are still lacking. Here, we sequenced the genomes of 37 animals from four Portuguese native sheep breeds (Campaniça, Bordaleira Serra da Estrela, Merino Branco and Merino Preto) and 19 crossbred sheep to make inferences on their genomic diversity and population structure. Mean genomic diversities were very similar across these breeds (.30 ≤ H_o ≤ .34; .30 ≤ H_e ≤ .35; 1.7 × 10^-3 ≤ π ≤ 3.1 × 10^-3) and the levels of inbreeding were negligible (.005 ≤ F_IS ≤ .038). The Principal Components, Bayesian clustering and Treemix analyses split the Portuguese breeds in two main groups which are consistent with historical records: one comprising Campaniça and Serra da Estrela together with other European and transboundary dairy breeds; and another of the well-differentiated multi-purpose Merino and Merino-related breeds. Runs of homozygosity analyses yielded 1,690 ROH segments covering an average of 2.27 Gb across the genome in all individuals. The overall genome covered by ROH segments varied from 27,75 Mb in Serra da Estrela to 61,29 Mb in Campaniça. The phylogenetic analysis of sheep mitogenomes grouped the Portuguese native breeds within sub-haplogroup B1a along with two animals of the Akkaraman breed from Turkey. This result provides additional support to a direct influence of Southwest Asian sheep in local breeds from the Iberian Peninsula. Our study is a first step pertaining to the genomic characterization of Portuguese sheep breeds and the results emphasize the potential of genomic data as a valid tool to guide conservation efforts in locally adapted sheep breeds. In addition, the genomic data we generated can be used to identify markers for breed assignment and traceability of certified breed-products.

Genetic Diversity of Creole Sheep Managed by Indigenous Communities of the Central Region of Veracruz, Mexico

Simple Summary

The genetic diversity in three populations of creole sheep managed by indigenous communities of the central region of Veracruz, Mexico, is reported from blood samples of 90 sheep, taken from the herds of indigenous families in the municipalities of Tehuipango, Astacinga, and Tlaquilpa, Veracruz, Mexico. The genomic DNA of the sheep was evaluated using four microsatellites amplified by PCR and visualized on polyacrylamide gels. The four microsatellites were polymorphic, observed heterozygosity was lower than the expected level, and the indices of endogamy indicated a slight diminution of homozygotes and the variation was hosted at the individual level. The evaluated sheep present a genetic diversity that is conserved across endogamic crosses, for which reason the design of a plan of protection and use of these sheep populations would permit their sustainable management.

Abstract

In the indigenous communities of central Veracruz, herds of creole sheep have been established and managed through traditional practices of crossing, but their genetic characteristics have never been examined in order to evaluate their state of endogamy, and to help the management programs to protect this genetic resource. The objective of the present study was to characterize the genetic diversity of three populations of creole sheep managed by indigenous communities in the central region of Veracruz, Mexico. Indigenous family producers of creole sheep were located and blood samples taken from 90 individual sheep from the municipalities of Tehuipango, Astacinga and Tlaquilpa, Veracruz. In the laboratory, the genomic DNA was extracted and genetic diversity characterized using four microsatellites (ILSTS11, ILSTS5, SRCRSP9 and OarFCB128) amplified by PCR and visualized on polyacrylamide gels. The four microsatellites were highly informative (PIC = 85%) and presented values of 0.6 to 0.81 of heterozygosity, with an average number of 16 alleles. According to the Hardy–Weinberg equilibrium model, three of the loci were not significant (p < 0.05), presumably this means that they do not deviate significantly from H–W predictions and there was slight genetic differentiation (F_ST = 0.025), along with a slight decrease in homozygotes (F_IS = −0.021). According to the analysis of variance, 99% of the total variation was hosted at the individual level. It is concluded that the three creole sheep populations still present genetic diversity at the four loci and non-random pairings have occurred.

How have sheep breeds differentiated from each other in Morocco? Genetic structure and geographical distribution patterns

Background

Based on the relatively homogeneous origin of the sheep breeds in Morocco that originate mainly from Iberia, it is highly relevant to address the question of how these very diverse sheep populations differentiated from each other. The Mountains of the High Atlas and Middle Atlas are expected to constitute North–South and West–East geographical barriers, respectively, which could have shaped the history of the differentiation of sheep breeds. The aim of this study was to test this hypothesis by considering the genetic structure and the spatial distribution of five major breeds (Sardi, Timahdite, Beni Guil, Boujaad and D’man) and one minor breed (Blanche de Montagne), by analysing the mtDNA control region, using 30 individuals per breed.

Results

Phylogenetic and network analyses did not indicate any clear separation among the studied breeds and discriminant component principal analysis showed some overlap between them, which indicates a common genetic background. The calculated pairwise F_ST values and Nei’s genetic distances revealed that most breeds showed a moderate genetic differentiation. The lowest and highest degrees of differentiation were retrieved in the Beni Guil and Boujaad breeds, respectively. Analysis of molecular variance (AMOVA) indicated that more than 95% of the genetic diversity occurs within individuals, while between- and within-population variabilities represent only 1.332% and 2.881%, respectively. Isolation-by-distance, spatial Principal Component Analysis (sPCA), and spatial AMOVA analyses evidenced clear examples of geographical structuration among the breeds, both between and within breeds. However, several enigmatic relationships remain, which suggest the occurrence of complex events leading to breed differentiation.

Conclusions

The approaches used here resulted in a convergent view on the hypothetic events that could have led to the progressive differentiation between the Moroccan breeds. The major split seems to be linked to the West–East barrier of the Middle Atlas, whereas the influence of the High Atlas is less obvious and incompletely resolved. The study of additional breeds that have settled near the High Atlas should clarify the relationships between the breeds of the West part of the country, in spite of their small population size.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00679-2.

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.

Genetic diversity and differentiation of 12 eastern Adriatic and western Dinaric native sheep breeds using microsatellites

Nuclear genetic diversity and differentiation of 341 sheep belonging to 12 sheep breeds from Croatia and Bosnia and Herzegovina were examined. The aim of the study was to provide the understanding of the genetic structure and variability of the analysed pramenka sheep populations, and to give indications for conservation strategies based on the population diversity and structure information. The genetic variation of the sheep populations, examined at the nuclear level using 27 microsatellite loci, revealed considerable levels of genetic diversity, similar to the diversity found in other European indigenous low-production sheep breeds. Population-specific alleles were detected at most loci and in breeds analysed. The observed heterozygosity ranged from 0.643 (in Lika pramenka) to 0.743 (in Vlasic pramenka), and the expected heterozygosity ranged from 0.646 (in Lika pramenka) to 0.756 (in Dalmatian pramenka). Significant inbreeding coefficients were found for half of the populations studied and ranged from 0.040 (Pag island sheep) to 0.091 (Kupres pramenka). Moderate genetic differentiation was found between the studied sheep populations. The total genetic variability observed between different populations was 5.29%, whereas 94.71% of the variation was found within populations. Cres island sheep, Lika pramenka and Istrian sheep were identified as the most distinct populations, which was confirmed by the factorial analysis of correspondence and supported through a bootstrapping adjustment to correct for the difference in the sample sizes. The population structure analysis distinguished 12 clusters for the 12 sheep breeds analysed. However, the cluster differentiation was low for Dalmatian, Vlasic, Stolac and Krk pramenka. This systematic study identified Lika pramenka and Rab island sheep as those with the lowest diversity, whereas Istrian sheep and Pag island sheep had the highest. Conservation actions are proposed for Istrian, Rab and Cres island sheep, Lika and Kupres pramenka because of high estimated coefficients of inbreeding.

Mini review: breeding Awassi and Assaf sheep for diverse management conditions

The Local Awassi, a triple-purpose breed for meat, milk, and carpet-wool production, is a low-prolific, hardy breed that is well adapted to the unfavorable conditions of the Middle East, where it is managed under traditionally extensive to semi-extensive conditions. Breeding work with the Awassi has included within-breed selection, crossbreeding, and gene introgression. Those efforts resulted in a variety of Awassi-derived genotypes that successfully occupy semi-intensive as well as intensive production systems. Thus, within-breed selection resulted in development of the "Improved Awassi"-a dairy-type Awassi strain which, under intensive management, produces over 500 l milk/ewe annually; crossbreeding with the East Friesian breed led to the development of the Assaf dairy breed, which exceeds the Improved Awassi in prolificacy and in year-round breeding activity, and introgression of the B allele of the FecB locus into the Awassi and Assaf breeds resulted in the formation of the prolific Afec Awassi and Afec Assaf strains, with prolificacies of 1.9 and 2.5 lambs born per ewe lambing, respectively. Advanced molecular genetics tools have enabled a better understanding of how the Awassi breed was formed during domestication and have uncovered differences in its genetic structure compared to other breeds. Implementing large-scale selection schemes that implement emerging new information on the sheep genome, overcoming threats of inbreeding depression, and further breeding for high uterine capacity are the new breeding goals for the Awassi, Assaf, and their derivatives.

Genetic variation and relationship of six Indian sheep breeds adapted to the northwestern arid zone of Rajasthan

This study illustrates the genetic diversity and relationships within and among six Indian sheep breeds of the northwestern arid region of Rajasthan, based on microsatellite markers. The range of allele diversity was 7.72-9.56, and gene diversity was 0.686-0.766, revealing that these breeds possessed substantial amounts of genetic diversity. Positive F (IS) values suggested a deficit of heterozygotes in all six breeds. Despite the declining status of the Marwari, Chokla, Jaisalmeri, Magra and Pugal breeds, an absence of a recent genetic bottleneck was evident from the data. The genetic differentiation estimate (F (ST) = 6.1%) suggested low levels of differentiation between the breeds. Genetic distance estimates revealed a close relationship between the Magra-Pugal and Nali-Jaisalmeri breed pairs. This information forms a framework for designing genetic management and conservation programs for these valuable ovine breeds.

Diversity analysis of sheep breeds from Southern peninsular and Eastern regions of India

Information is presented on the genetic diversity and relationship among six Indian sheep breeds/populations belonging to the Southern peninsular and Eastern agroecological zones, based on microsatellite markers. Parameters of genetic variation, viz., allele diversity, observed heterozygosity, gene diversity and population inbreeding estimates, were calculated for the six breeds. The allele diversity ranged from 6.40 to 7.92, whereas the gene diversity varied from 0.617 to 0.727. The highest allele and gene diversity was observed for Nellore sheep, while the lowest was exhibited by Garole breed. Within population inbreeding estimate (F(IS)) revealed a significant deficit of heterozygotes in Deccani, Madgyal, Nellore and Garole, whereas Ganjam and Chhotanagpuri sheep showed an excess of heterozygotes. The contribution of each breed to the total diversity of the breeds was quantified by the Weitzman approach. The marginal loss of diversity incurred with removal of Nellore and Garole breeds was higher (>27%), whereas removal of Deccani breed resulted in lowest loss of diversity (3.84%) from the set. Estimation of the genetic differentiation (F(ST)) and genetic distance (D(A)) between the pairs of breeds revealed a close relationship between Deccani and Madgyal sheep (F(ST) = 0.017; D(A) = 0.080) and greatest demarcation between Madgyal and Garole breeds (F(ST) = 0.110; D(A) = 0.622). The information generated would help in shaping genetic management and conservation programs for the sheep breeds under consideration.

Genetic diversity in farm animals--a review

Domestication of livestock species and a long history of migrations, selection and adaptation have created an enormous variety of breeds. Conservation of these genetic resources relies on demographic characterization, recording of production environments and effective data management. In addition, molecular genetic studies allow a comparison of genetic diversity within and across breeds and a reconstruction of the history of breeds and ancestral populations. This has been summarized for cattle, yak, water buffalo, sheep, goats, camelids, pigs, horses, and chickens. Further progress is expected to benefit from advances in molecular technology.

Microsatellite analysis to estimate genetic relationships among five bulgarian sheep breeds

Herein, genetic relationships among five breeds of Bulgarian sheep were estimated using microsatellite markers. The total number of alleles identified was 226 at the 16 loci examined. D_A distance values were used for phylogenetic tree construction with the UPGMA algorithm. The two Tsigai and two Maritza populations were found to be geneticallvery closely related to each other y (0.198, and 0.258 respectively). The Pleven Black Head population was distinct from the other four. These results could be useful for preserving genes in these breeds, thereby ensuring their preservation in Bulgaria.

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 in mazandaranian native cattle: a comparison with Holstein cattle, using ISSR marker

This study was carried out to investigate genetic diversity in Mazandaranian native cattle population comparised to the Holstein breed, using Inter Simple Sequence Repeats (ISSR) marker. A total of 175 animals, including 71 native and 104 cattle of Holstein breed were screened. The extraction of DNA samples were carried out, using modified salting out method. A 19-mer oligonucleotide, (GA)9C, was used as primer in PCR reactions. The PCR products showed 15 different fragments with length ranged from 120 to 1600 bp in the two breeds.. Genetic variation indexes, including effective number of alleles, Shannon index, Nei's gene diversity and standard genetic distance were estimated, using POPGene software. Generally, the estimated genetic variation indexes showed low levels of diversity in the two breeds. However, Nei's gene diversity and Shannon index estimation was observed almost two folds in native cattle compared to Holstein breed. Less levels of diversity in Holstein cattle may be because of applying intensive selection programs. Conversely, native cattle have been less affected by selection. Therefore, it seems that Mazandaranian native cattle probably are better for breeding programs than Holstein cattle. Results showed that ISSR Markers are reliable and can be used in genetic diversity investigations.