In the current study, the genetic architecture of growth and linear type traits were investigated in Akkaraman sheep. Estimations of genomic heritability, genetic correlations, and phenotypic correlations were implemented for 17 growth and linear type traits of 473 Akkaraman lambs by the univariate and multivariate analysis of animal mixed models. Correspondingly, moderate heritability estimates, as well as high and positive genetic/phenotypic correlations were found between growth and type traits. On the other hand, 2 genome-wide and 19 chromosome-wide significant single nucleotide polymorphisms were found to be associated with the traits as a result of animal mixed model-based genome-wide association analyses. Accordingly, we propose several genes located on different chromosomes (e.g., PRDM2, PTGDR, PTPRG, KCND2, ZNF260, CPE, GRID2, SCD5, SPIDR, ZNF407, HCN3, TMEM50A, FKBP1A, TLE4, SP1, SLC44A1, and MYOM3) as putative quantitative trait loci for the 22 growth and linear type traits studied. In our study, specific genes (e.g., TLE4, PTGDR, and SCD5) were found common between the traits studied, suggesting an interplay between the genetic backgrounds of these traits. The fact that four of the proposed genes (TLE4, MYOM3, SLC44A1, and TMEM50A) are located on sheep chromosome 2 confirms the importance of these genomic regions for growth and morphological structure in sheep. The results of our study are therefore of great importance for the development of efficient selection indices and marker-assisted selection programs, as well as for the understanding of the genetic architecture of growth and linear traits in sheep.

The objective of this study was to uncover the genetic background of wool quality, a production trait, by estimating genomic heritability and implementing GWAS in Akkaraman sheep. The wool characteristics measured included fibre diameter (FD) and staple length (SL) at the age of 8 months and yearling fibre diameter (YFD), yearling staple length (YSL) and yearling greasy fleece weight (YGFW) at 18 months of age. Animals were genotyped using the Axiom 50 K Ovine Genotyping Array. Maximum likelihood estimations of a linear mixed model (LMM) were used to estimate genomic heritability, where GWAS was conducted following a score test of each trait. Genomic heritability estimates for the traits ranged between 0.22 and 0.63, indicating that phenotypes have a moderate range of heritability. One genome- and six chromosome-wide significant SNPs were associated with the wool traits in Akkaraman lambs. Accordingly, TRIM2, MND1, TLR2, RNF175, CEP290, TMTC3, RERE, SLC45A1, SOX2, MORN1, SKI, FAAP20, PRKCZ, GABRD, CFAP74, CALML6 and TMEM52 genes as well as nine uncharacterized regions (LOC101118971, LOC105609137, LOC105603067, LOC101122892, LOC106991694, LOC106991467, LOC106991455, LOC105616534 and LOC105609719) were defined as plausible candidates. The findings of this study shed light on the genetics of wool quality and yield for the Akkaraman breed and suggests targets for breeders during systematic breeding programmes.

Genome-wide association studies (GWAS) have been used as an effective tool to understand the genetics of complex traits such as gastrointestinal parasite (GIP) resistance. The aim of this study was to understand the genetics of gastrointestinal parasite (nematodes, Moniezia spp., Eimeria spp.) resistance in Akkaraman sheep by performing genomic heritability estimations and conducting GWAS to uncover responsible genomic regions. This is one of the first studies to examine the genetic resistance of Akkaraman sheep to the tapeworm parasite. The samples from 475 animals were genotyped using the Axiom 50K Ovine Genotyping Array. Genomic heritability estimates ranged from 0.00 to 0.34 for parasite resistance traits. This indicates that measured phenotypes have low to moderate heritability estimates. A total of two genome-wide significant SNP associated with TNEM3 and ATRNL1 genes and 10 chromosome-wide significant SNPs related with 10 genes namely NELL1, ST6GALNAC3, HIPK1, SYT1, ALK, ZNF596, TMCO5A, PTH2R, LARGE1, and SCG2 were suggested as candidates for parasite resistance traits. The majority of these candidate genes were involved in several basic biological processes that are essential and important for immune system functions and cellular growth; specifically, inflammatory responses, cellular transport, cell apoptosis, cell differentiation, histone de-acetylation, and endocytosis. These results have implications for animal breeding program studies due to the effect that the genetic background has on parasite resistance, which underlies many productive, health, and wellness-related traits.

The domestic sheep with over 1200 breeds descended from those early domesticated animals that are bred for a variety of resources such as meat, milk and wool. Akkaraman, a fat‐tailed indigenous sheep breed of Türkiye, is widespread throughout Central Anatolia, with the largest indigenous sheep population. Assessing the genetic diversity and genomic structure of animal breeds is among the key contributors to deciphering adaptation to environmental extremes and constructing efficient genetic improvement strategies. Therefore, this study aimed to characterize the genome of Akkaraman breed against various world‐renowned transboundary sheep and indigenous sheep with fat and thin tails. Genetic similarities and differences between those breeds have been displayed by estimating and comparing various genetic diversity indices, linkage disequilibrium (LD) estimates and fixation index (F_ST), runs of homozygosity (ROH) as well as PCA and neighbour‐joining tree analysis. Akkaraman sheep were observed to form a cluster alongside Moghani, Karakas, Tibetan and Cyprus Fat Tail sheep, which are primarily the sole representatives of fat‐tailed sheep in the study. This clustering was evident in both the PCA and neighbour‐joining tree analysis. The Akkaraman sheep was also observed to have the lowest genomic inbreeding and one of the lowest numbers of ROHs, which might also indicate that the breed has not been exposed to historical intensive selection pressure, inbred mating or a massive population bottleneck that might leave strong marks of genomic homozygosity. The results improve our understanding of the genetic diversity in Akkaraman sheep in comparison with certain mainstream sheep breeds as well as those indigenous breeds from around the world. Additionally, findings will also provide valuable insights to perform further GWAS effectively by considering population structure, diversity and LD patterns observed among the breeds while providing practical knowledge that will contribute to designing efficient and successful genome‐based selection programmes for worldwide sheep production systems.