Combined haplotype blocks regression and multi-locus mixed model analysis reveals novel candidate genes associated with milk traits in dairy sheep

Insights into the genomic homogeneity of Moroccan indigenous sheep breeds though the lens of runs of homozygosity

Numerous studies have indicated that Morocco's indigenous sheep breeds are genetically homogenous, posing a risk to their survival in the challenging harsh climate conditions where they predominantly inhabit. To understand the genetic behind genetic homogeneity through the lens of runs of homozygosity (ROH), we analyzed the whole genome sequences of five indigenous sheep breeds (Beni Guil, Ouled Djellal, D'man, Sardi, Timahdite and Admixed).The results from principal component, admixture, Fst, and neighbour joining tree analyses consistently showed a homogenous genetic structure. This structure was characterized by an average length of 1.83 Mb for runs of homozygosity (ROH) segments, with a limited number of long ROH segments (24-48 Mb and > 48 Mb). The most common ROH segments were those ranging from 1-6 Mb. The most significant regions of homozygosity (ROH Islands) were mostly observed in two chromosomes, namely Chr1 and Chr5. Specifically, ROH Islands were exclusively discovered in the Ouled Djellal breed on Chr1, whereas Chr5 exhibited ROH Islands in all breeds. The analysis of ROH Island and iHS technique was employed to detect signatures of selection on Chr1 and Chr5. The results indicate that Chr5 had a high level of homogeneity, with the same genes being discovered across all breeds. In contrast, Chr1 displays some genetic variances between breeds. Genes identified on Chr5 included SLC39A1, IL23A, CAST, IL5, IL13, and IL4 which are responsible for immune response while genes identified on Chr1 include SOD1, SLAMF9, RTP4, CLDN1, and PRKAA2. ROH segment profile and effective population sizes patterns suggests that the genetic uniformity of studied breeds is the outcome of events that transpired between 250 and 300 generations ago. This research not only contributes to the understanding of ROH distribution across breeds but helps design and implement native sheep breeding and conservation strategies in Morocco. Future research, incorporating a broader sample size and utilizing the pangenome for reference, is recommended to further elucidate these breeds' genomic landscapes and adaptive mechanisms.

Genome-Wide SNP Analysis for Milk Performance Traits in Indigenous Sheep: A Case Study in the Egyptian Barki Sheep

Simple Summary

The Barki sheep is one of the three main breeds in Egypt, which is spread mainly throughout the northwestern coastal zone, which has harsh conditions. Considering the harsh, semi-arid habitat of this breed, milk performance traits such as milk yield and milk composition have a very important role in the feeding of newborn lambs and affect their growth during the early stage of life. In this study, rare milk performance data and genomic information of Barki sheep were used to uncover diversified genomic regions that could explain the variability of milk yield and milk quality traits in the studied population of Barki ewes. Genome-wide analysis identified genomic regions harboring interesting candidate genes such as SLC5A8, NUB1, TBC1D1, KLF3 and ABHD5 for milk yield and PPARA and FBLN1 genes for milk quality traits. The findings offer valuable information for obtaining a better understanding of the genetics of milk performance traits and contribute to the genetic improvement of these traits in Barki sheep.

Abstract

Sheep milk yield and milk composition traits play an important role in supplying newborn lambs with essential components such as amino acids, energy, vitamins and immune antibodies and are also of interest in terms of the nutritional value of the milk for human consumption. The aim of this study was to identify genomic regions and candidate genes for milk yield and milk composition traits through genome-wide SNP analyses between high and low performing ewes of the Egyptian Barki sheep breed, which is well adapted to the harsh conditions of North-East Africa. Therefore, out of a herd of 111 ewes of the Egyptian Barki sheep breed (IBD = 0.08), ewes representing extremes in milk yield and milk quality traits (n = 25 for each group of animals) were genotyped using the Illumina OvineSNP50 V2 BeadChip. The fixation index (F_ST) for each SNP was calculated between the diversified groups. F_ST values were Z-transformed and used to identify putative SNPs for further analysis (Z(F_ST) > 10). Genome-wide SNP analysis revealed genomic regions covering promising candidate genes related to milk performance traits such as SLC5A8, NUB1, TBC1D1, KLF3 and ABHD5 for milk yield and PPARA and FBLN1 genes for milk quality trait. The results of this study may contribute to the genetic improvement of milk performance traits in Barki sheep breed and to the general understanding of the genetic contribution to variability in milk yield and quality traits.