A comprehensive investigation combining pooled whole-genome sequencing, extensive morphological evaluations, carcass composition assays, and targeted gene expression profiling has produced the first genomic atlas for five economically vital Egyptian sheep populations. Domestic sheep represent one of humanity\u2019s earliest pastoral domesticates, having emerged in the Fertile Crescent between Eastern Anatolia and the Zagros Mountains roughly eight to ten thousand years ago before expanding globally into a primary livestock species. In Egypt, where the national flock comprises approximately 5.5 million head and contributes significantly to the domestic red meat supply, indigenous breeds make up more than 65 percent of the sheep sector. To decipher the genetic architecture governing their unique productive attributes and local adaptations, an international research team analyzed Barki, Rahmani, Ossimi, and Awassi purebred lines alongside a first-generation crossbred population derived from Rahmani dams and Barki sires. The resulting dataset resolved millions of genomic variants, pinpointed prominent structural variant hotspots on specific autosomes, and illuminated how evolutionary pressures on metabolic and physiological pathways shape divergent breed profiles.<\/p>\n
The study cohort was established from foundation stock collected across five distinct northern Egyptian geographic zones, specifically Alexandria Governorate, the Salloum Plateau, Sakha City, and two separate locales within Matrouh Province including the El-Hammam district. To eliminate environmental confounding and nutritional variation, experimental animals were maintained under a standardized semi-intensive management system with shared barn housing, open range access, and identical feeding regimens at Elmabrok Farm in El-Hammam. For the primary whole-genome sequencing architecture, researchers sampled 120 adult male rams, selecting 24 unrelated individuals per breed cohort that were born within a narrow three-week window to strictly standardize age at 26.50 months. Venous blood was collected aseptically, and high-molecular-weight genomic DNA was extracted using silica-membrane columns, followed by rigorous quality verification through agarose gel electrophoresis and spectrophotometric absorbance ratios near 1.82. Equal quantities of purified DNA from the 24 animals in each group were combined into equimolar pools, providing an established, cost-effective framework for population-scale allele frequency estimation and variant discovery across broad pastoral cohorts.<\/p>\n
Genomic libraries were prepared from sheared 350-base-pair DNA fragments using enzymatic end-repair, unique dual-indexed adapter ligation, and bead-based size selection, followed by high-throughput paired-end sequencing on an Illumina NovaSeq 6000 platform. The sequencing campaign generated more than 470 gigabases of total sequence data, producing an average of 600 million paired-end reads per breed pool and achieving mean autosomal coverage depths ranging from 24.2x to 41.3x. Awassi achieved the highest coverage depth at 41.3x, with 78.3 percent of its reference genome covered at thirty-fold depth or greater, while Barki exhibited a notably elevated GC content of 50.7 percent and a mapping efficiency of 92.8 percent when aligned to the ARS-UI_Ramb_v3.0 sheep reference assembly. Bioinformatic filtering through fastp, Picard duplicate marking, and GATK HaplotypeCaller identified between 11.9 million single-nucleotide polymorphisms in Ossimi and 17.4 million in Rahmani. Across all examined breeds, the genomic distribution of polymorphisms remained highly conserved, with intronic segments harboring approximately 38.8 percent of all single-nucleotide variants and small insertions or deletions, followed closely by intergenic regions.<\/p>\n
Population-level heterozygosity metrics derived from pooled sequencing allele frequency spectra revealed pronounced contrasts in genetic variation that reflect distinct breeding histories. Purebred populations retained substantial genetic diversity, with Rahmani and Barki exhibiting population heterozygosity rates of 85.01 percent and 83.80 percent, whereas the first-generation Barki and Rahmani crossbred displayed a reduced heterozygosity estimate of 60.41 percent, representing a classic genetic signature of hybrid lineage structure. Genome-wide pairwise fixation index calculations across ten-kilobase windows revealed moderate to high genetic differentiation, ranging from 0.18 between Barki and Rahmani to 0.35 between Barki and Awassi. Nucleotide diversity mirrored these divergence patterns, spanning from 0.0021 in Ossimi to 0.0028 in Rahmani, with Barki showing an intermediate diversity value of 0.0025. Across all five cohorts, transition-to-transversion ratios held steady near 2.3, matching established mammalian baseline expectations and confirming the high technical fidelity of the sequence variant catalog.<\/p>\n
Functional classification of detected variants using Clusters of Orthologous Groups identified prominent, shared enrichment in categories governing lipid transport, lipid metabolism, and signal transduction mechanisms across every breed. Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathway mappings consolidated these functional observations, demonstrating highly significant enrichment with probability values below 0.0001 for the peroxisome proliferator-activated receptor signaling pathway and fatty acid degradation. Beyond lipid dynamics, functional annotations highlighted pathways responding to xenobiotic stimuli and chemical defense mechanisms, reflecting evolutionary pressures exerted by regional foraging environments and browse flora. Independent structural variant mapping executed through DELLY and copy number profiling revealed that structural variations were non-randomly clustered into distinct genomic hotspots on chromosomes 2, 6, and 18. These structural hotspots directly overlapped critical candidate loci, including fatty acid binding protein 4 and myostatin on chromosome 2, the keratin-associated protein gene cluster and bone morphogenetic protein receptor type 1B on chromosome 6, and the agouti signaling protein gene on chromosome 18.<\/p>\n
Phenotypic measurements systematically documented divergence across the broader experimental flock of 1,262 animals, confirming marked anatomical and physiological specialization among the breeds. Standardized body condition scores evaluated by three independent assessors demonstrated significant variation at probability values below 0.001, with Rahmani and the crossbred population scoring in the fat-to-obese range at 4.312 and 4.533, respectively. Conversely, Barki sheep displayed a lean body condition score of 2.534, whereas Ossimi and Awassi occupied intermediate positions with scores of 3.221 and 3.118. Detailed morphometric profiling of twenty-five linear body dimensions in one hundred adult animals per breed confirmed that Awassi and crossbred sheep possessed the largest skeletal frames, with Awassi recording the highest wither height at 67.59 centimeters and body length at 74.53 centimeters. In contrast, Barki animals ranked among the smallest across major skeletal dimensions, reflecting a compact conformation suited to arid pastoral mobility, while Rahmani presented the widest chest breadth at 21.42 centimeters and longest tail measurement at 20.53 centimeters.<\/p>\n
Detailed post-slaughter evaluations of twenty representative rams slaughtered at twenty-seven months of age demonstrated substantial, statistically significant disparities in meat yield and carcass fat distribution. Prior to slaughter, animals were subjected to an eighteen-hour fast, rendered unconscious via intravenous sodium thiopental at twenty milligrams per kilogram of live body weight, and humanely exsanguinated under surgical anesthesia in compliance with veterinary euthanasia standards. Carcass processing revealed that crossbred animals and purebred Rahmani yielded the highest trimmed meat weights, producing 23.95 kilograms and 18.40 kilograms of lean muscle tissue, respectively, alongside superior dressing percentages near 53.5 percent. Rahmani carcasses accumulated the largest volume of dissectible fat at 3.72 kilograms. Chemical proximate analysis of the longissimus dorsi muscle further demonstrated that Rahmani exhibited the highest intramuscular fat concentration at 4.20 percent, while crossbred muscle tissue showed the greatest crude protein density at 22.35 percent, and Awassi demonstrated an exceptional lean-to-fat ratio of 58.60.<\/p>\n
Quantitative reverse transcription polymerase chain reaction assays performed on skeletal muscle tissue from the twenty slaughtered rams functionally validated the transcriptional activity of candidate genes linked to production and muscularity. Normalizing target cycle thresholds against glyceraldehyde-3-phosphate dehydrogenase and beta-actin reference standards revealed that myostatin and growth differentiation factor 8 expression were elevated by 1.67-fold and 1.38-fold in Rahmani and crossbred muscle tissues. In contrast, insulin-like growth factors 1 and 2 exhibited higher relative expression in Barki and Ossimi, showing 1.51-fold and 1.74-fold variations across breeds. Upstream endocrine regulators of the somatotropic axis, specifically growth hormone and its receptor, demonstrated maximal expression in crossbred and Rahmani animals with 1.85-fold and 2.05-fold increases, while the proteolytic regulator calpain 3 was most active in Rahmani and Awassi, displaying a 1.34-fold elevation that corresponds with known differences in muscle protein turnover and meat tenderness characteristics.<\/p>\n
Transcriptional profiling of lipid metabolism genes corroborated the phenotypic fat deposition measurements observed across the fat-tailed and leaner populations. Key lipogenic enzymes governing de novo fatty acid synthesis, including acetyl-CoA carboxylase alpha and fatty acid synthase, exhibited significant upregulation in fat-tailed Rahmani and Awassi breeds, varying by 1.54-fold and 1.26-fold across cohorts. Stearoyl-CoA desaturase and lipoprotein lipase reached peak expression in Rahmani, displaying 1.58-fold and 1.28-fold differences, while fatty acid binding protein 4 showed a 1.31-fold elevation in Rahmani and Awassi. Circulating adipokines also mirrored these energetic divergences, as leptin expression varied 2.40-fold with highest levels in crossbred animals, and adiponectin showed a 1.39-fold elevation centered in Rahmani. Additional metabolic regulators, including cell death-inducing DFFA-like effector A, thyroid hormone responsive protein, and carnitine palmitoyltransferase 1, demonstrated coordinated, breed-specific expression shifts ranging from 1.36-fold to 1.64-fold, establishing a direct transcriptional correlation between enriched lipid pathways and physical carcass fat accretion.<\/p>\n
While this genomic atlas provides an unprecedented foundation for understanding small ruminant biodiversity in North Africa, the authors emphasized several methodological parameters and interpretative limitations. Because all experimental animals were maintained at a single commercial facility in Matrouh Governorate, the observed genetic distributions may not capture all localized micro-evolutionary variations present across wider Egyptian geographic zones. Furthermore, pooled sequencing designs yield population-level allele frequencies rather than individual genotypes, meaning that formal haplotype-based selection tests such as integrated haplotype scores or cross-population extended haplotype homozygosity could not be computed. Future investigations will require individual-level whole-genome resequencing across diverse regional flocks, broader crossbreeding designs, and high-density genome-wide association studies to definitively establish causal variant mechanisms. Nonetheless, these integrated findings provide essential molecular targets, including lipid and myogenic regulatory markers, to support marker-assisted selection programs and safeguard indigenous genetic resources against accelerating climate challenges.<\/p>\n
Subject of Research: Breed-pool whole-genome sequencing and phenotypic characterization of five Egyptian sheep populations<\/p>\n
Article Title: First breed-pool whole genome sequencing of egyptian sheep: a comprehensive genomic atlas revealing diversity and candidate genes for production and adaptation<\/p>\n
Article References: Hassanine, N. N. A. M., Amin, A. H., Hafez, E. E., Ahmed, M. M. M., Ismail, A. H., Mohamed, A. Z., El Halim, H. A., El-Sayed, G. M., Adss, I. A. A., Radwan, K. H., Abo-Zaid, M. A., Rashad, A. M. A., Husien, H. M., Saleh, A. A., & Wang, M. (2026). First breed-pool whole genome sequencing of egyptian sheep: a comprehensive genomic atlas revealing diversity and candidate genes for production and adaptation. BMC Biotechnology, 26(1), Article 96. https:\/\/doi.org\/10.1186\/s12896-026-01199-x<\/a><\/p>\n Image Credits: AI Generated<\/p>\n