Genetic structure and admixture in sheep from terminal breeds in the United States

Analysing the genetic diversity of three sheep breeds in Turkey and nearby countries using 50 K SNPs data

This study analysed the genetic diversity and population structure of eight sheep breeds in Turkey and nearby countries. Moderate genetic diversity was observed, with the Sakiz (SKZ) exhibiting the highest diversity based on heterozygosity and allelic richness (AR) values. Genetic distances revealed differentiation between the populations, with the most significant divergence between the Cyprus Fat Tail (CFT) and SKZ breeds. PCA demonstrated SKZ and Chios (CHI) clustering together, indicating genetic similarity. Karakas (KRS), Norduz (NDZ), Afshari (AFS), Moghani (MOG) and others showed overlap, reflecting genetic relationships. Ancestry analysis found that KRS was predominantly inherited from the second ancestral population, while SKZ and NDZ were primarily derived from the first and second ancestral lineages. This illustrated the populations' diverse origins. Most genetic variation (96.84%) was within, not between, populations. The phi-statistic (PhiPT) indicated moderate differentiation overall. Phylogenetic analysis further demonstrated the genetic distinctiveness of the SKZ breed. ROH and FROH analyses showed that SKZ exhibited the highest homozygosity and inbreeding, while KRS displayed the lowest. This study elucidates these breeds' genetic diversity, structure and relationships. Key findings include moderate diversity, evidence of differentiation between breeds, diverse ancestral origins and distinct ROH patterns. This provides insights into the population's genetic characteristics and conservation requirements.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

DIA-Based Proteomic Analysis Reveals MYOZ2 as a Key Protein Affecting Muscle Growth and Development in Hybrid Sheep

Hybridization of livestock can be used to improve varieties, and different hybrid combinations produce unique breeding effects. In this study, male Southdown and Suffolk sheep were selected to hybridize with female Hu sheep to explore the effects of male parentage on muscle growth and the development of offspring. Using data-independent acquisition technology, we identified 119, 187, and 26 differentially abundant proteins (DAPs) between Hu × Hu (HH) versus Southdown × Hu (NH), HH versus Suffolk × Hu (SH), and NH versus SH crosses. Two DAPs, MYOZ2 and MYOM3, were common to the three hybrid groups and were mainly enriched in muscle growth and development-related pathways. At the myoblast proliferation stage, MYOZ2 expression decreased cell viability and inhibited proliferation. At the myoblast differentiation stage, MYOZ2 expression promoted myoblast fusion and enhanced the level of cell fusion. These findings provide new insights into the key proteins and metabolic pathways involved in the effect of male parentage on muscle growth and the development of hybrid offspring in sheep.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

Global DNA Methylation, miRNA, and mRNA Profiles in Sheep Skeletal Muscle Promoted by Hybridization

With the development of high-throughput sequencing technology, several nongenetic variations, including noncoding RNAs such as miRNAs, and DNA methylation, have been found to play an important role in animal muscle development and fat metabolism. In this study, Southdown and Suffolk were selected as male parents for hybridization with Hu sheep (Southdown × Hu (NH), Suffolk × Hu (SH), and Hu × Hu (HH)). RNA sequencing, bisulfite sequencing, and small-RNA sequencing were used to study the methylation patterns and differences in miRNA and mRNA expression in the F1 sheep longissimus dorsi muscle tissue. We identified 765 differentially expressed genes (DEGs), 10,161 differentially methylated regions (DMRs), and 164 differentially expressed miRNAs, which were significantly enriched in AMPK signaling, fatty acid degradation, metabolism, and other related pathways (P < 0.05). In addition, we constructed a DNA methylation-mRNA and miRNA-mRNA coexpression network. A total of 42 common genes were identified from DMRs and DEGs. Importantly, we predicted that 33 differentially expressed miRNAs directly or indirectly targeted the SLC27A6. The data obtained in this study provide useful information and evidence to support further understanding of the miRNA and DNA methylation of key genes regulating muscle growth and fat metabolism in hybrid sheep populations.

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Bioclimatic Zoning for Sheep Farming through Geostatistical Modeling in the State of Pernambuco, Brazil

The Intergovernmental Panel on Climate Change (IPCC) has pointed out the high vulnerability of developing countries to climate change, which is expected to impact food and income security. Sheep farming is one of the main animal productions among the families located in the most vulnerable regions of the semiarid region of Pernambuco state, a Brazilian territory known for its high temperatures, low relative humidity, and high net solar radiation. Therefore, the objective of this study was to identify different regions of Pernambuco that may be more suitable for different breeds of sheep, based on non-parametric statistics and kriging maps of the temperature and humidity index (THI). THI values were determined based on mean annual temperature and wind speed extracted from the TerraClimate remote sensing database. Pernambuco state presented THI values ranging from 66 to 79, with the hair breeds having a high potential for exploitation in almost all territories, including the main meat-producing breeds. The East Friesian breed, a high milk producer, would be well suited to the Agreste mesoregion, a territory that, like the Pajeú and Moxotó microregions, also proved favorable for the introduction of three wool breeds (Suffolk, Poll Dorset, and Texel) known as major meat producers. The kriging maps of the THI values successfully allowed the identification of strategic development regions of Pernambuco state with high potential for sheep breeding.

Genetic association of wool quality characteristics in United States Rambouillet sheep

Introduction: Fine wool production is an important source of revenue, accounting for up to 13% of total revenue in extensively managed wool sheep production systems of the United States. The Rambouillet are a predominant breed that excels in wool quality characteristics. Understanding the genetic basis of wool quality characteristics would aid in the development of genomic breeding strategies to facilitate genetic improvement. Methods: Wool characteristics and DNA were collected for rams enrolled in the North Dakota State University and University of Wyoming annual central performance ram tests over a three-year period (2019-2021, N = 313). The relationships of wool quality characteristics including grease fleece weight adjusted 365 days (wt. 365 adj.), clean fleece wt. 365 adj., staple length 365 adj., average fiber diameter, face wool cover, amount of skin wrinkles and belly wool were evaluated through genome-wide association studies (GWAS), Pearson correlation and ANOVA. Results: The GWAS identified four genome-wide significant genetic markers (p-value <1.19e-06) and five chromosome-wide significant markers (p-value <1.13e-05) on chromosomes 1, 2, 4, 15, and 19. Significant markers were associated with genes notable for relevant wool biological functions, including the gene ABCC8 which codes for SUR1, an ATP-sensitive potassium channel known to affect hair growth and 60S ribosomal protein L17-like, previously found to be expressed during follicle formation. The strongest Pearson correlation coefficients were identified between clean fleece wt. 365 adj. and grease fleece wt. 365 adj. (r = 0.83) and between clean fleece wt. 365 adj. and staple length 365 adj. (r = 0.53). Additionally, clean fleece wt. 365 adj. was correlated with final body weight (r = 0.35) and scrotal circumference (r = 0.16). Staple length 365 adj. (p-value = 5e-04), average fiber diameter (p-value = .0053) and clean fleece wt. 365 adj. (p-value = .014) were significantly associated with belly wool score. Discussion: The results of this study provide important insight into the relationships between wool quality characteristics and report specific markers that Rambouillet sheep producers may use to help inform selection and breeding decisions for improved wool quality.

Using whole genome sequence to compare variant callers and breed differences of US sheep

As whole genome sequence (WGS) data sets have become abundant and widely available, so has the need for variant detection and scoring. The aim of this study was to compare the accuracy of commonly used variant calling programs, Freebayes and GATK HaplotypeCaller (GATK-HC), and to use U.S. sheep WGS data sets to identify novel breed-associated SNPs. Sequence data from 145 sheep consisting of 14 U.S. breeds were filtered and biallelic single nucleotide polymorphisms (SNPs) were retained for genotyping analyses. Genotypes from both programs were compared to each other and to genotypes from bead arrays. The SNPs from WGS were compared to the bead array data with breed heterozygosity, principal component analysis and identifying breed associated SNPs to analyze genetic diversity. The average sequence read depth was 2.78 reads greater with 6.11% more SNPs being identified in Freebayes compared to GATK-HC. The genotype concordance of the variant callers to bead array data was 96.0% and 95.5% for Freebayes and GATK-HC, respectively. Genotyping with WGS identified 10.5 million SNPs from all 145 sheep. This resulted in an 8% increase in measured heterozygosity and greater breed separation in the principal component analysis compared to the bead array analysis. There were 1,849 SNPs identified in only the Romanov sheep where all 10 rams were homozygous for one allele and the remaining 135 sheep from 13 breeds were homozygous for the opposite allele. Both variant calling programs had greater than 95% concordance of SNPs with bead array data, and either was suitably accurate for ovine WGS data sets. The use of WGS SNPs improved the resolution of PCA analysis and was critical for identifying Romanov breed-associated SNPs. Subsets of such SNPs could be used to estimate germplasm composition in animals without pedigree information.

Genome-Wide Genetic Diversity and Population Structure of Local Sudanese Sheep Populations Revealed by Whole-Genome Sequencing

Local Sudanese sheep populations inhabiting diverse environmental conditions and holding opposing morphologies provide opportunities for molecular-genetic research. Characterizing their genome is crucial for sustainable breeding improvement and targeting favorable genes in breeding programs. However, the genome of these sheep populations, which comprises several subtypes, remains uncharacterized using whole-genome sequence data. This study aimed to elucidate genome-wide genetic diversity and population structure of 11 local Sudanese sheep populations, namely, Hammari, Kabbashi, Meidobe, Ashgar, Dubasi, Watish, Bega, Naili, Fulani, Zagawi, and Garag. Ninety whole blood samples were collected, and we extracted DNA using a Qiagen DNeasy® extraction kit. We used the Illumina HiSeq 2000 platform to sequence all the DNA samples. We included whole-genome sequence data of three Ethiopian sheep (Doyogena, Kefis, and Gafera) and one Libyan sheep (Libyan Barbary) in the study to infer the genetic relationships of local Sudanese sheep populations from a continental perspective. A total of 44.8 million bi-allelic autosomal SNPs were detected; 28.5% and 63.3% occur in introns and intergenic regions, respectively. The mean genetic diversity ranged from 0.276 for Garag to 0.324 for Kabbashi sheep populations. The lowest FST estimates were observed between Kabbashi and Ashgar and the highest between Bega and Fulani local Sudanese sheep populations. The principal component and population structure analyses of the 11 local Sudanese sheep populations indicated three separate genetic groups categorized following their tail morphotype, geographical distribution, and population subtype. The thin-tailed local Sudanese sheep populations exhibited independent clustering from the fat-tailed Ethiopian and Libyan sheep. We also observed distinct clustering between the fat-tailed Ethiopian and Libyan sheep. The present study’s findings demonstrated the population structure and principal components related to tail morphotype, geographical distribution, and population subtype of local Sudanese sheep populations. A clear signature of admixture was observed among the studied local Sudanese sheep populations.

Transcriptome-metabolome analysis reveals how sires affect meat quality in hybrid sheep populations

Crossbreeding improves and enhances meat quality and is widely used in sheep production; however, the molecular mechanisms underlying the meat quality of various crossbred sheep remain unknown. In this study, male Southdown, Suffolk and Hu sheep were crossbred with female Hu sheep, and the transcriptomes and metabolomes of the longissimus dorsi muscle of the F1 generation were sequenced to explore how different sire breeds affect meat quality. The results showed that 631 differentially expressed genes and 119 significantly altered metabolites contributed to muscle development characteristics and meat quality-related diversity (P < 0.05). These genes and metabolites were significantly enriched in lipid metabolism pathways, including arachidonic acid metabolism and PPAR signaling. Several candidate genes were associated with muscle growth, such as MYLK3, MYL10, FIGN, MYH8, MYOM3, LMCD1, and FLRT1. Among these, MYH8 and MYL10 participated in regulating muscle growth and development and were correlated with meat quality-related fatty acid levels (|r| > 0.5 and p < 0.05). We selected mRNA from four of these genes to verify the accuracy of the sequencing data via qRT-PCR. Our findings provide further insight into the key genes and metabolites involved in muscle growth and meat quality in hybrid sheep populations.

Baseline T-lymphocyte and cytokine indices in sheep peripheral blood

Background

Sheep are an important livestock species worldwide and an essential large-animal model for animal husbandry and veterinary research. Understanding fundamental immune indicators, especially T-lymphocyte parameters, is necessary for research on sheep diseases and vaccines, to better understand the immune response to bacteria and viruses for reducing the use of antibiotics and improving the welfare of sheep. We randomly selected 36 sheep of similar ages to analyze cell-related immune indicators in peripheral blood mononuclear cells (PBMCs). The proportions of CD4⁺ and CD8⁺ T cells in PBMCs were detected by flow cytometry. We used Concanavalin A (Con A) and Phorbol-12-myristate-13-acetate (PMA)/Ionomycin to stimulate PBMCs, and measured the expression of IFN-γ, IL-4, and IL-17A using enzyme-linked immunosorbent assay (ELISA) and enzyme-linked immunospot assay (ELISpot). Simultaneously, PMA/Ionomycin/brefeldin A (BFA) was added to PBMCs, then the expression of IFN-γ, IL-4, and IL-17A was detected by flow cytometry after 4 h of culturing. In addition, we observed the proliferation of PBMCs stimulated with Con A for 3, 4, and 5 days.

Results

The proportions of CD4⁺ T lymphocytes (18.70 ± 4.21%) and CD8⁺ T lymphocytes (8.70 ± 3.65%) were generally consistent among individuals, with a CD4/CD8 ratio of 2.40 ± 0.79. PBMCs produced high levels of IFN-γ, IL-4, and IL-17A after stimulation with PMA/Ionomycin and Con A. Furthermore, PMA/Ionomycin stimulation of PBMC yielded significantly higher cytokine levels than Con A stimulation. Flow cytometry showed that the level of IFN-γ (51.49 ± 11.54%) in CD8⁺ T lymphocytes was significantly (p < 0.001) higher than that in CD4⁺ T lymphocytes (14.29 ± 3.26%); IL-4 (16.13 ± 6.81%) in CD4⁺ T lymphocytes was significantly (p < 0.001) higher than that in CD8⁺ T lymphocytes (1.84 ± 1.33%), There was no difference in IL-17A between CD4⁺ (2.83 ± 0.98%) and CD8⁺ T lymphocytes (1.34 ± 0.67%). The proliferation of total lymphocytes, CD4⁺ T lymphocytes, and CD8⁺ T lymphocytes continued to increase between days 3 and 5; however, there were no significant differences in proliferation between the cell types during the stimulation period.

Conclusions

Evaluating primary sheep immune indicators, especially T lymphocytes, is significant for studying cellular immunity. This study provided valuable data and theoretical support for assessing the immune response of sheep to pathogens and improving sheep welfare.

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.

Natural Hybrid Origin of the Controversial "Species" Clematis × pinnata (Ranunculaceae) Based on Multidisciplinary Evidence

Interspecific hybridization is common and has often been viewed as a driving force of plant diversity. However, it raises taxonomic problems and thus impacts biodiversity estimation and biological conservation. Although previous molecular phylogenetic studies suggested that interspecific hybridization may be rather common in Clematis, and artificial hybridization has been widely applied to produce new Clematis cultivars for nearly two centuries, the issue of natural hybridization of Clematis has never been addressed in detail. In this study, we tested the hybrid origin of a mesophytic and cold-adapted vine species, Clematis pinnata, which is a rare and taxonomically controversial taxon endemic to northern China. Using field investigations, flow cytometry (FCM), phylogenomic analysis, morphological statistics, and niche modeling, we tested hybrid origin and species status of C. pinnata. The FCM results showed that all the tested species were homoploid (2n = 16). Phylonet and HyDe analyses based on transcriptome data showed the hybrid origins of C. × pinnata from either C. brevicaudata × C. heracleifolia or C. brevicaudata × C. tubulosa. The plastome phylogeny depicted that C. × pinnata in different sampling sites originated by different hybridization events. Morphological analysis showed intermediacy of C. × pinnata between its putative parental species in many qualitative and quantitative characters. Niche modeling results suggested that C. × pinnata had not been adapted to a novel ecological niche independent of its putative parents. These findings demonstrated that plants of C. × pinnata did not formed a self-evolved clade and should not be treated as a species. The present study also suggests that interspecific hybridization is a common mechanism in Clematis to generate diversity and variation, and it may play an important role in the evolution and diversification of this genus. Our study implies that morphological diversity caused by natural hybridization may overstate the real species diversity in Clematis.

Evolution of the sheep industry and genetic research in the United States: opportunities for convergence in the twenty-first century

The continuous development and application of technology for genetic improvement is a key element for advancing sheep production in the United States. The US sheep industry has contracted over time but appears to be at a juncture where a greater utilization of technology can facilitate industry expansion to new markets and address inefficiencies in traditional production practices. Significant transformations include the increased value of lamb in relation to wool, and a downtrend in large-scale operations but a simultaneous rise in small flocks. Additionally, popularity of hair breeds not requiring shearing has surged, particularly in semi-arid and subtropical US environments. A variety of domestically developed composite breeds and newly established technological approaches are now widely available for the sheep industry to use as it navigates these ongoing transformations. These genetic resources can also address long-targeted areas of improvement such as growth, reproduction and parasite resistance. Moderate progress in production efficiency has been achieved by producers who have employed estimated breeding values, but widespread adoption of this technology has been limited. Genomic marker panels have recently shown promise for reducing disease susceptibility, identifying parentage and providing a foundation for marker-assisted selection. As the ovine genome is further explored and genomic assemblies are improved, the sheep research community in the USA can capitalize on new-found information to develop and apply genetic technologies to improve the production efficiency and profitability of the sheep industry.

British Sheep Breeds as a Part of World Sheep Gene Pool Landscape: Looking into Genomic Applications

Sheep farming has been an important sector of the UK's economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.

Placental Characteristics Classification of Various Native Turkish Sheep Breeds

Simple Summary

The aim of this study was to classify placental characteristics of Akkaraman, Morkaraman, Karayaka, Awassi, Malya, and Bafra native sheep breeds using the hierarchical clustering method. As a result, six breeds were separated into three clusters: the first cluster consisted of Bafra, Karayaka, and Awassi breeds; the second consisted of Akkaraman and Malya breeds; and the third cluster included only the Morkaraman breed.

Abstract

The aim of this study was to classify placental characteristics of Akkaraman, Morkaraman, Karayaka, Awassi, Malya, and Bafra sheep breeds using the hierarchical clustering method. In total, 240 individual data records were used as experimental material. Placental characteristics such as total cotyledon surface area, small and large cotyledon length, small cotyledon depth, etc. were used as explanatory variables to classify the breeds’ characteristics. Hierarchical clustering was used with the nearest neighbour method with Euclidean distance in order to classify the sheep breeds’ variations. As a result, six breeds were separated into three clusters: the first cluster consisted of Bafra, Karayaka, and Awassi breeds; the second consisted of Akkaraman and Malya breeds; and the third cluster included only the Morkaraman breed. Bafra and Karayaka were pointed as the nearest breeds, with a similarity of 98.7% in terms of placental characteristics. The similarity rate of the Akkaraman and Malya breeds was at a level of 97.5%, whereas it was 96.8% for Bafra, Karayaka, and Awassi breeds. The similarity of Akkaraman, Karayaka, Awassi, Malya, and Bafra sheep breeds was estimated as 95.7%. The overall similarity was found to be at a level of 93.2% among sheep breeds. The outcomes of the study might be useful as a selection tool for reproductivity and can be used to select the breed to be reared.

Clustering patterns mirror the geographical distribution and genetic history of Lemnos and Lesvos sheep populations

Elucidating the genetic variation and structure of Lemnos and Lesvos sheep is critical for maintaining local genetic diversity, ecosystem integrity and resilience of local food production of the two North Aegean islands. In the present study, we explored genetic diversity and differentiation as well as population structure of the Lemnos and Lesvos sheep. Furthermore, we sought to identify a small panel of markers with the highest discriminatory power to assign animals across islands. A total number of n = 424 (n = 307, Lemnos and n = 117, Lesvos) ewes, sampled from n = 24 herds dispersed at different geographic regions on the two islands, were genotyped with the 50K SNP array. Mean observed heterozygosity was higher (but not statistically significantly different) in Lesvos than in Lemnos population (0.384 vs. 0.377) while inbreeding levels were higher in Lemnos than Lesvos herds (0.065 vs. 0.031). Results of principal components along with that of admixture analysis and estimated genetic distances revealed genetic clusters corresponding to Lesvos and Lemnos origin and the existence of infrastructure within islands that were associated with geographical isolation and genetic history of the studied populations. In particular, genetic analyses highlighted three geographically isolated herds in Lemnos that are located at mountainous areas of the island and are characterized as representatives of the local sheep by historic data and reports. Admixture analysis also showed a shared genetic background between Lemnos and Lesvos sheep attributable to past gene flow. Little overall genetic differentiation was detected between the two island sheep populations, while 150 discriminatory SNPs could accurately assign animals to their origin. Present results are comparable with those reported in the worldwide sheep breeds, suggesting geography related genetic patterns across and within islands and the existence of the local Lemnos sheep.