Comparative study of the reproductive and growth performance of the Hamra and Rumbi ovine breeds.

Growth and reproductive performance traits are important economic indicators for analyzing the overall performance of breeding systems. This study aims to evaluate the comparative performance of two Algerian sheep (Rumbi and Hamra) in terms of growth and reproductive performance, and the effect of factors such as breed, season of birth, mode of birth and age of the mother on this performance in a semi-intensive breeding system. The reproductive performance of 577 Rumbi ewes and 1328 Hamra ewes bred at the Tiaret and Saïda experimental stations respectively, was analyzed using performance monitoring data. Fertility rates for the Rumbi and Hamra breeds of 87.14% and 78.8% respectively were practically similar (p > 0,05). Litter size at birth and weaning was significantly higher in the Hamra breed than in the Rumbi breed (p < 0,05). Weaning mortality was significantly higher in the Hamra breed than in the Rumbi breed, with an average of 22.60% versus 14.94% (p < 0,05). The effect of factors showed that there was a highly significant effect of the mother's age and season of birth on the reproductive performance of the Hamra and Rumbi breeds with a p < 0.0001 on fertility, litter size at birth, litter size at weaning and fertility. There was a significant effect of the year factor on reproductive performance with p = 0,013 for the Hamra breed and p = 0,031 for the Rumbi breed. The results of this study showed that Rumbi lambs were heavier at birth than Hamra lambs. The values observed were 4,86 kg versus 3,10 kg for the Hamra breed, with a highly significant difference (p < 0,0001), so that the average daily pre-weaning weight gains of Rumbi lambs were higher than those of Hamra lambs, at 0,195 kg/day versus 0,113 kg/day for Hamra lambs, with a high significance (p < 0,0001). The effect of factors showed that there was a significant effect of the mother's age on the ADGs (0-30), (30-70) and (70-90) of the Hamra and Rumbi breeds with a p = 0,034 and p = 0,02 respectively. There was also a highly significant effect of the birth mode effect on ADGs (0-30), (30-70) and (70-90) only for the Hamra breed with a p = 0,004. The effect of the birth weight on ADGs was not significant for both Hamra and Rumbi breeds with a p > 0,05. According to the findings of this study, the Hamra breed had superior reproductive potential and the Rumbi breed had superior growth. The Hamra breed showed better maternal skills in terms of fertility and prolificacy, while the Rumbi breed showed better lamb growth performance. Consequently, these results could be used for selective sheep breeding, taking into account the random effects of the environment and the potential of each breed.

Genome-wide discovery of selection signatures in four Anatolian sheep breeds revealed by ddRADseq.

High-density genomic data analyzed by accurate statistical methods are of potential to enlighten past breeding practices such as selection by unraveling fixed regions. In this study, four native Turkish sheep breeds (80 samples) were genotyped via 296.097 single nucleotide polymorphisms (SNPs) detected by double-digest restriction site-associated DNA (ddRADseq) library preparation combined with the Illumina HiSeq X Ten instrument in order to identify genes under selection pressure. A total of 32, 136, 133, and 119 protein-coding genes were detected under selection pressure by runs of homozygosity (ROH), integrated haplotype score (iHS), the ratio of extended haplotype homozygosity (Rsb), and fixation index (FST) approaches, respectively. Of these, a total of 129 genes were identified by at least two statistical models which overlapped with a total of 52 quantitative trait loci (QTL)-associated SNPs, known to be related to fiber diameter, milk content, body weight, carcass traits, some blood parameters, and entropion. A total of six genes under selection pressure were validated by three statistical approaches five of which are of potential to be integrated into animal breeding since they were associated with wool fiber diameter (ZNF208B), behaviors related to neurocognitive development (CBX1 and NFE2L1), adaptation to high-altitude (SDK1), and anxiety causing internal stress (GSG1L). The sixth gene (COPZ1) turned out to play an important role in coping with different types of cancer in mammals. In particular, ROH analysis uncovered significant findings that the Güney Karaman (GKR) had experienced different selection practices than the Akkaraman (AKR) breed. Moreover, some genes specifically under selection in the GKR breed turned out to be associated with olfaction (OR6K6, OR6N1, OR6N2, and OR4C16), survival during the gestation period (PRR15L), and heat stress (CDK5RAP9). The results of this study imply that GKR may become genetically different from the AKR breed at the genome level due to most probably experiencing different adaptation processes occurring in raised climatic conditions. These differences should be conserved to face future challenges, while other native Turkish sheep breeds could be monitored via genome-wide high-density SNP data to obtain deeper knowledge about the effects of natural selection.

Segregation GWAS to linearize a non-additive locus with incomplete penetrance: an example of horn status in sheep.

BACKGROUND: The objective of this study was to introduce a genome-wide association study (GWAS) in conjunction with segregation analysis on monogenic categorical traits. Genotype probabilities calculated from phenotypes, mode of inheritance and pedigree information, are expressed as the expected allele count (EAC) (range 0 to 2), and are inherited additively, by definition, unlike the original phenotypes, which are non-additive and could be of incomplete penetrance. The EAC are regressed on the single nucleotide polymorphism (SNP) genotypes, similar to an additive GWAS. In this study, horn phenotypes in Merino sheep are used to illustrate the advantages of using the segregation GWAS, a trait believed to be monogenic, affected by dominance, sex-dependent expression and likely affected by incomplete penetrance. We also used simulation to investigate whether incomplete penetrance can cause prediction errors in Merino sheep for horn status. RESULTS: Estimated penetrance values differed between the sexes, where males showed almost complete penetrance, especially for horned and polled phenotypes, while females had low penetrance values for the horned status. This suggests that females homozygous for the 'horned allele' have a horned phenotype in only 22% of the cases while 78% will be knobbed or have scurs. The GWAS using EAC on 4001 animals and 510,174 SNP genotypes from the Illumina Ovine high-density (600k) chip gave a stronger association compared to using actual phenotypes. The correlation between the EAC and the allele count of the SNP with the highest -log10(p-value) was 0.73 in males and 0.67 in females. Simulations using penetrance values found by the segregation analyses resulted in higher correlations between the EAC and the causative mutation (0.95 for males and 0.89 for females, respectively), suggesting that the most predictive SNP is not in full LD with the causative mutation. CONCLUSIONS: Our results show clear differences in penetrance values between males and female Merino sheep for horn status. Segregation analysis for a trait with mutually exclusive phenotypes, non-additive inheritance, and/or incomplete penetrance can lead to considerably more power in a GWAS because the linearized genotype probabilities are additive and can accommodate incomplete penetrance. This method can be extended to any monogenic controlled categorical trait of which the phenotypes are mutually exclusive.

Population structure and breed identification of Chinese indigenous sheep breeds using whole genome SNPs and InDels.

BACKGROUND: Accurate breed identification is essential for the conservation and sustainable use of indigenous farm animal genetic resources. In this study, we evaluated the phylogenetic relationships and genomic breed compositions of 13 sheep breeds using SNP and InDel data from whole genome sequencing. The breeds included 11 Chinese indigenous and 2 foreign commercial breeds. We compared different strategies for breed identification with respect to different marker types, i.e. SNPs, InDels, and a combination of SNPs and InDels (named SIs), different breed-informative marker detection methods, and different machine learning classification methods. RESULTS: Using WGS-based SNPs and InDels, we revealed the phylogenetic relationships between 11 Chinese indigenous and two foreign sheep breeds and quantified their purities through estimated genomic breed compositions. We found that the optimal strategy for identifying these breeds was the combination of DFI_union for breed-informative marker detection, which integrated the methods of Delta, Pairwise Wright's FST, and Informativeness for Assignment (namely DFI) by merging the breed-informative markers derived from the three methods, and KSR for breed assignment, which integrated the methods of K-Nearest Neighbor, Support Vector Machine, and Random Forest (namely KSR) by intersecting their results. Using SI markers improved the identification accuracy compared to using SNPs or InDels alone. We achieved accuracies over 97.5% when using at least the 1000 most breed-informative (MBI) SI markers and even 100% when using 5000 SI markers. CONCLUSIONS: Our results provide not only an important foundation for conservation of these Chinese local sheep breeds, but also general approaches for breed identification of indigenous farm animal breeds.

Transcriptomic and Metabolomic Analyses Reveal Molecular Regulatory Networks for Pigmentation Deposition in Sheep.

Domestic animals have multiple phenotypes of skin and coat color, which arise from different genes and their products, such as proteins and metabolites responsible with melanin deposition. However, the complex regulatory network of melanin synthesis remains to be fully unraveled. Here, the skin and tongue tissues of Liangshan black sheep (black group) and Liangshan semi-fine-wool sheep (pink group) were collected, stained with hematoxylin-eosin (HE) and Masson-Fontana, and the transcriptomic and metabolomic data were further analyzed. We found a large deposit of melanin granules in the epidermis of the black skin and tongue. Transcriptome and metabolome analysis identified 744 differentially expressed genes (DEGs) and 443 differentially expressed metabolites (DEMs) between the pink and black groups. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses revealed the DEGs and DEMs were mainly enriched in the pathways of secondary metabolic processes, melanin biosynthesis processes, melanin metabolism processes, melanosome membranes, pigment granule membranes, melanosome, tyrosine metabolism, and melanogenesis. Notably, we revealed the gene ENSARG00020006042 may be a family member of YWHAs and involved in regulating melanin deposition. Furthermore, several essential genes (TYR, TYRP1, DCT, PMEL, MLANA, SLC45A2) were significantly associated with metabolite prostaglandins and compounds involved in sheep pigmentation. These findings provide new evidence of the strong correlation between prostaglandins and related compounds and key genes that regulate sheep melanin synthesis, furthering our understanding of the regulatory mechanisms and molecular breeding of pigmentation in sheep.

ASReml-based estimation of early genetic parameters in hybrid mutton sheep populations.

The aim of this study was to analyze the effects of non-genetic factors on the estimation of genetic parameters of early growth traits in hybrid mutton sheep using ASReml software, in order to provide theoretical basis for screening the optimal hybriding combinations and accelerating the breeding process of new breeds of specialized housed-feeding mutton sheep. We selected the wellgrown hybrid Southhu (Southdown × Hu sheep) and Dorhu (Dorset × Hu sheep) sheep as the research objects, constructed weight correction formulae for SH and DH sheep at 60 and 180 days; and used ASReml software to investigate the effects of non-genetic factors on the estimation of genetic parameters of early growth traits in hybrid sheep. The results showed that the birth month and birth type were found significant for all traits (p < 0.001); the heritability of maternal effects ranged from 0.0709 to 0.1859. It was found that both SH and DH sheep emerged the potential for rapid early growth and development, early growth traits are significantly affected by maternal genetic effects, thereby the maternal effect should be taken into consideration for the purpose of improving accuracy in parameter estimations and therefore increasing the success of breeding programs.

Analysis of the genetic variance of fibre diameter measured along the wool staple for use as a potential indicator of resilience in sheep.

BACKGROUND: The effects of environmental disturbances on livestock are often observed indirectly through the variability patterns of repeated performance records over time. Sheep are frequently exposed to diverse extensive environments but currently lack appropriate measures of resilience (or sensitivity) towards environmental disturbance. In this study, random regression models were used to analyse repeated records of the fibre diameter of wool taken along the wool staple (bundle of wool fibres) to investigate how the genetic and environmental variance of fibre diameter changes with different growing environments. RESULTS: A model containing a fifth, fourth and second-order Legendre polynomial applied to the fixed, additive and permanent environmental effects, respectively, was optimal for modelling fibre diameter along the wool staple. The additive genetic and permanent environmental variance both showed variability across the staple length trajectory. The ranking of sire estimated breeding values (EBV) for fibre diameter was shown to change along the staple and the genetic correlations decreased as the distance between measurements along the staple increased. This result suggests that some genotypes were potentially more resilient towards the changes in the growing environment compared to others. In addition, the eigenfunctions of the random regression model implied the ability to change the fibre diameter trajectory to reduce its variability along the wool staple. CONCLUSIONS: These results show that genetic variation in fibre diameter measured along the wool staple exists and this could be used to provide greater insight into the ability to select for resilience in extensively raised sheep populations.

Integration analysis of pituitary proteome and transcriptome reveals fertility-related biomarkers in FecB mutant Small Tail Han sheep.

The Booroola fecundity mutation (FecB) in Small Tail Han sheep has been shown to enhance ovulation rates and litter sizes by affecting the hypothalamic-pituitary-gonadal (HPG) axis. Despite the pituitary's role in reproductive regulation, its involvement in FecB-induced ovulation remains understudied. Our study aimed to fill this gap by analyzing pituitary tissues from FecB homozygous (BB) and wild-type (WW) ewes during luteal and follicular phases using tandem mass tag-based protein quantification and the DIABLO framework for proteomic and transcriptomic data integration. Significant differences in 277 proteins were observed across estrus periods, with network analysis highlighting the voltage-dependent calcium channel L-type alpha-1C as a key convergence point in oxytocin signaling and GnRH secretion pathways. The DIABLO method revealed a strong correlation (0.98) between proteomic and transcriptomic datasets, indicating a coordinated response in FecB ewes. Notably, higher expression levels of Follicle Stimulating Hormone Subunit Beta (FSHB) and Luteinizing Hormone Subunit Beta (LHB) were found in BB ewes during the follicular phase, potentially due to elevated E2 concentrations. Furthermore, our analysis identified genes related to the Gamma-aminobutyric acid type A receptor family (GABRA2, GABRG1, GABRB1) in the pituitary, with GABRB1 showing higher expression in BB ewes. This suggests a role for GABA in modulating GnRH and gonadotropin feedback loops, potentially contributing to the FecB mutation's effect on ovulation. This study provides novel insights into the pituitary's role in fertility among FecB sheep, identifying GABA as a potential regulatory factor within the HPG axis. The findings also open avenues for discovering new biomarkers in pituitary endocrinology for sheep breeding purposes.

Impact of truncating historical data on prediction ability of dairy sheep selection candidates.

Along the last decades, the genetic evaluation methodology has evolved, improving breeding value estimates. Many breeding programmes have historical phenotypic records and large number of generations, but to make use of them could result in more inconveniences than benefits. In this study, the prediction ability of genotyped young animals was assessed by simultaneously evaluating the removal of historical data, two pedigree deepness and two methodologies (traditional BLUP and single-step genomic BLUP or ssGBLUP), using milk yield records of 40 years of three Latxa dairy sheep populations. The linear regression method was used to compare predictions of young rams before and after progeny testing, with six cut-off points, by intervals of 4 years (from 1992 to 2012), and statistics of ratio of accuracies, bias, and dispersion were calculated. The prediction accuracy of selection candidates, when genomic information was included, was the highest in all Latxa populations (between 0.54 and 0.69 with full data set). Nevertheless, the deletion of historical phenotypic data resulted on moderate accuracy gain in the bigger data size populations (mean gain 2.5%), and the smaller population took advantage of a moderate data deletion (2.7% gain by removing data until 2004), reducing accuracy when more records were removed. The bias of validation individuals was lower when the breeding value was predicted based on genomic information (between 2.1 and 13.9), being lower when the biggest amount of data was deleted in the bigger data size populations (5.2% reduction), and the smaller population was benefited from data deletion between 1996 and 2008 (3.8% bias reduction). Meanwhile, the slope of estimated genetic trend was lower when less data were included, and an overestimation of the unknown parent group estimates was observed. The results indicated that ssGBLUP evaluations were outstanding, compared with traditional BLUP evaluations, while the depth of pedigree had a very small influence, and deletion of historical phenotypic data was beneficial. Thus, Latxa routine genetic evaluations would benefit from truncating phenotypic records between 2000 and 2004, the use of two pedigree generations and the implementation of ssGBLUP methodology.

Milk transcriptome biomarker identification to enhance feed efficiency and reduce nutritional costs in dairy ewes.

In recent years, rising prices for high-quality protein-based feeds have significantly increased nutrition costs. Consequently, investigating strategies to reduce these expenses and improve feed efficiency (FE) have become increasingly important for the dairy sheep industry. This research investigates the impact of nutritional protein restriction (NPR) during prepuberty and FE on the milk transcriptome of dairy Assaf ewes (sampled during the first lactation). To this end, we first compared transcriptomic differences between NPR and control ewes. Subsequently, we evaluated gene expression differences between ewes with divergent FE, using feed conversion ratio (FCR), residual feed intake (RFI), and consensus classifications of high- and low-FE animals for both indices. Lastly, we assess milk gene expression as a predictor of FE phenotype using random forest. No effect was found for the prepubertal NPR on milk performance or FE. Moreover, at the milk transcriptome level, only one gene, HBB, was differentially expressed between the NPR (n = 14) and the control group (n = 14). Further, the transcriptomic analysis between divergent FE sheep revealed 114 differentially expressed genes (DEGs) for RFI index (high-FERFI = 10 vs low-FERFI = 10), 244 for FCR (high-FEFCR = 10 vs low-FEFCR = 10), and 1 016 DEGs between divergent consensus ewes for both indices (high-FEconsensus = 8 vs low-FEconsensus = 8). These results underscore the critical role of selected FE indices for RNA-Seq analyses, revealing that consensus divergent animals for both indices maximise differences in transcriptomic responses. Genes overexpressed in high-FEconsensus ewes were associated with milk production and mammary gland development, while low-FEconsensus genes were linked to higher metabolic expenditure for tissue organisation and repair. The best prediction accuracy for FE phenotype using random forest was obtained for a set of 44 genes consistently differentially expressed across lactations, with Spearman correlations of 0.37 and 0.22 for FCR and RFI, respectively. These findings provide insights into potential sustainability strategies for dairy sheep, highlighting the utility of transcriptomic markers as FE proxies.

Whole-genome scan for selection signature associated with temperature adaptation in Iranian sheep breeds.

The present study aimed to identify the selection signature associated with temperature adaptation in Iranian sheep breeds raised in cold and hot environments. The Illumina HD ovine SNP600K BeadChip genomic arrays were utilized to analyze 114 animals from eight Iranian sheep breeds, namely Ghezel, Afshari, Shall, Sanjabi, Lori-Bakhtiari, Karakul, Kermani, and Balochi. All animals were classified into two groups: cold-weather breeds and hot-weather breeds, based on the environments to which they are adapted and the regions where they have been raised for many years. The unbiased FST (Theta) and hapFLK tests were used to identify the selection signatures. The results revealed five genomic regions on chromosomes 2, 10, 11, 13, and 14 using the FST test, and three genomic regions on chromosomes 10, 14, and 15 using the hapFLK test to be under selection in cold and hot groups. Further exploration of these genomic regions revealed that most of these regions overlapped with genes previously identified to affect cold and heat stress, nervous system function, cell division and gene expression, skin growth and development, embryo and skeletal development, adaptation to hypoxia conditions, and the immune system. These regions overlapped with QTLs that had previously been identified as being associated with various important economic traits, such as body weight, skin color, and horn characteristics. The gene ontology and gene network analyses revealed significant pathways and networks that distinguished Iranian cold and hot climates sheep breeds from each other. We identified positively selected genomic regions in Iranian sheep associated with pathways related to cell division, biological processes, cellular responses to calcium ions, metal ions and inorganic substances. This study represents the initial effort to identify selective sweeps linked to temperature adaptation in Iranian indigenous sheep breeds. It may provide valuable insights into the genomic regions involved in climate adaptation in sheep.

Effect of DGAT1 gene polymorphisms in coarse-haired fat-tailed lambs of different genotypes.

The aim of this study was to investigate the DGAT1 gene polymorphism and its effects on lamb weight in kazakh and tajik sheep breeds. A total of 97 blood samples were collected from purebred (еdilbay х еdilbay) and crossbred lambs (еdilbay x gissar) breеd by the Baiserke Agro Scientific and Production Center in the Talgar District of the Almaty Region of Kazakhstan. Animals were genotyped for DGAT1-AluI polymorphism using the polymerase chain reaction-restriction length polymorphism (PCR-RFLP) method. The result of PCR-RFLP showed that purebred (еdilbay х еdilbay) sheep had three genotypes (CC, CT and TT) and crossbred sheep had two genotypes (CC and CT). The predominant genotype was CC with a frequency of 0.70 and 0.58 in purebred sheep and crossbred sheep breeds, respectively. The DGAT1 gene showed no significant association with live weight of lambs at different times in both breeds studied. However, the study showed that the CC genotype produced higher live weight at day 60 in purebred sheep (CC: 33,668 kg and CT: 32,444) and at day 120 (CC: 41,487 and CT: 40,929) in crossbred lambs. The present study was the first to investigate the polymorphism and relationships between genotypes and lamb live weights for DGAT1 gene in sheep breeds, purebred and crossbred. We conclude that further comprehensive investigations should be done for the exact evidence of the effects of DGAT1/Alui polymorphism on lamb live weights.

Genome-wide association studies for citric and lactic acids in dairy sheep milk in a New Zealand flock.

The objectives of this study were to estimate genetic parameters for citric acid content (CA) and lactic acid content (LA) in sheep milk and to identify the associated candidate genes in a New Zealand dairy sheep flock. Records from 165 ewes were used. Heritability estimates based on pedigree records for CA and LA were 0.65 and 0.33, respectively. The genetic and phenotypic correlations between CA and LA were strong-moderate and negative. Estimates of genomic heritability for CA and LA were also high (0.85, 0.51) and the genomic correlation between CA and LA was strongly negative (-0.96 ± 0.11). No significant associations were found at the Bonferroni level. However, one intragenic SNP in C1QTNF1 (chromosome 11) was associated with CA, at the chromosomal significance threshold. Another SNP associated with CA was intergenic (chromosome 15). For LA, the most notable SNP was intragenic in CYTH1 (chromosome 11), the other two SNPs were intragenic in MGAT5B and TIMP2 (chromosome 11), and four SNPs were intergenic (chromosomes 1 and 24). The functions of candidate genes indicate that CA and LA could potentially be used as biomarkers for energy balance and clinical mastitis. Further research is recommended to validate the present results.

Multi-omic Analyses Shed Light on The Genetic Control of High-altitude Adaptation in Sheep.

Sheep were domesticated in the Fertile Crescent and then spread globally, where they have been encountering various environmental conditions. The Tibetan sheep has adapted to high altitudes on the Qinghai-Tibet Plateau over the past 3000 years. To explore genomic variants associated with high-altitude adaptation in Tibetan sheep, we analyzed Illumina short-reads of 994 whole genomes representing ∼ 60 sheep breeds/populations at varied altitudes, PacBio High fidelity (HiFi) reads of 13 breeds, and 96 transcriptomes from 12 sheep organs. Association testing between the inhabited altitudes and 34,298,967 variants was conducted to investigate the genetic architecture of altitude adaptation. Highly accurate HiFi reads were used to complement the current ovine reference assembly at the most significantly associated ß-globin locus and to validate the presence of two haplotypes A and B among 13 sheep breeds. The haplotype A carried two homologous gene clusters: (1) HBE1, HBE2, HBB-like, and HBBC, and (2) HBE1-like, HBE2-like, HBB-like, and HBB; while the haplotype B lacked the first cluster. The high-altitude sheep showed highly frequent or nearly fixed haplotype A, while the low-altitude sheep dominated by haplotype B. We further demonstrated that sheep with haplotype A had an increased hemoglobin-O2 affinity compared with those carrying haplotype B. Another highly associated genomic region contained the EGLN1 gene which showed varied expression between high-altitude and low-altitude sheep. Our results provide evidence that the rapid adaptive evolution of advantageous alleles play an important role in facilitating the environmental adaptation of Tibetan sheep.

Unravelling the Genetic Architecture of Serum Biochemical Indicators in Sheep.

Serum biochemical indicators serve as vital proxies that reflect the physiological state and functions of different organs. The genetic parameters and molecular mechanisms underlying serum biochemical indicators of sheep (Ovis aries) have not been well understood. Therefore, the aim of the present study was to identify the genetic architecture and genomic loci underlying ten serum biochemical indicators in sheep, including alanine transaminase, aspartate transferase, lactate dehydrogenase, cholesterol, glucose, phosphorus, calcium, creatinine, urea and total protein levels. We implemented genetic parameter estimations and GWASs for each trait in 422 Akkaraman lambs. Overall, low to moderate heritability estimates were found in the range of 0.14-0.55. Additionally, low to high genetic correlations were observed among traits. In total, 23 SNP loci were associated with serum biochemical indicators leading to 19 genes. These were SPTA1, MGST2, CACUL1, IGFBP7, PARD3, PHB1, SLC15A5, TRIM35, RGS6, NUP93, CNTNAP2, SLC7A11, B3GALT5, DPP10, HST2ST1, NRP1, LRP1B, MAP3K9 and ENSOARG00020040484.1, as well as LOC101103187, LOC101117162, LOC105611309 and LOC101118029. To our knowledge, these data provide the first associations between SPTA1 and serum cholesterol and between ENSOARG00020040484.1 and serum glucose. The current findings provide a comprehensive inventory of the relationships between serum biochemical parameters, genetic variants and disease-relevant characteristics. This information may facilitate the identification of therapeutic targets and fluid biomarkers and establish a strong framework for comprehending the pathobiology of complex diseases as well as providing targets for sheep genetic improvement programs.

Effects of KRTAP20-1 Gene Variation on Wool Traits in Chinese Tan Sheep.

Chinese Tan sheep lambs are recognised for having tight 'spring-like' curly wool when young, but this phenotype disappears with age. This wool consists of shorter, fine wool fibres (which are usually unmedullated) and heterotypic hair fibres (which are frequently medullated), which are referred to as 'halo hair'. Both the wool and hair fibres consist of α-keratin proteins embedded in a keratin-associated protein (KAP) matrix. Of these KAPs, the KAP20-1 gene (designated KRTAP20-1) and its effect on four fibre traits (mean fibre curvature, mean fibre diameter, fibre diameter standard deviation, and coefficient of variation of fibre diameter) of Tan lambs was studied. Seven previously identified KRTAP20-1 variants (A, B, D, E, F, G, and H) of KRTAP20-1 were revealed, but the previously identified variant C was not present. Of the seven variants detected, only two (A and G) were common and present at frequencies greater than 5%, and the effect of these on the fibre traits of the finer wool fibres was assessed. It was found that variant G was associated with an increased mean fibre curvature in these wool fibres. This suggests that KRTAP20-1 might possibly be expressed differentially in the two fibre types, which may be of future value in breeding.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert.

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

A meta-analysis approach for annotation and identification of lncRNAs controlling perirenal fat deposition in suckling lambs.

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals.

Sheep litter size heredity basis using genome-wide selective analysis.

Sheep are important herbivorous domestic animal globally, and the Chinese indigenous sheep breed has a multitude of economically significant variations due to the diverse geographical and ecological conditions. In particular, certain native breeds exhibit a visible high litter size phenotype due to the selection pressure of natural and artificial for thousands of years, offering an ideal animal model for investigating sheep's fecundity. In this study, selective signal analysis was performed on public whole-genome sequencing data from 60 sheep across eight breeds to identify candidate genes related to litter size. Results revealed that a total of 34,065,017 single-nucleotide polymorphisms (SNPs) were identified from all sheep, and 65 candidate genes (CDGs) were pinpointed from the top 1% of interacted windows and SNPs between the pairwise fixation index (FST, >0.149543) and cross-population extended haplotype homozygosity (XP-EHH, >0.701551). A total of 41 CDGs (e.g. VRTN, EYA2 and MCPH1) were annotated to 576 GO terms, of which seven terms were directly linked to follicular and embryonic development (e.g. TBXT, BMPR1B, and BMP2). In addition, 73 KEGG pathways were enriched by 21 CDGs (e.g. ENTPD5, ABCD4 and RXFP2), mainly related to Hippo (TCF4, BMPR1B and BMP2), TGF-ß (BMPR1B and BMP2), PI3K-Akt (ITGB4, IL4R and PPP2R5A) and Jak-STAT signalling pathways (IL20RA and IL4R). Notably, a series of CDGs was under strong selection in sheep with high litter size traits. These findings result could improve the comprehension of the genetic underpinnings of sheep litter size. Furthermore, it provides valuable CDGS for future molecular breeding.

Chromosome-level genome assembly of Guide Black-Fur sheep (Ovis aries).

Guide Black-Fur sheep (GD) is a breed of Tibetan sheep (Ovis aries) that lives in the Qinghai-Tibetan plateau region at an altitude of over 4,000 m. However, a lack of genomic information has made it difficult to understand the high-altitude adaptation of these sheep. We sequenced and assembled the GD reference genome using PacBio, Hi-C, and Illumina sequencing technologies. The final assembled genome size was 2.73 Gb, with a contig N50 of 20.30 Mb and a scaffold N50 of 107.63 Mb. The genome is predicted to contain 20,759 protein-coding genes, of which 98.42 have functional annotations. Repeat elements account for approximately 52.2% of the genomic landscape. The completeness of the GD genome assembly is highlighted by a BUSCO score of 93.1%. This high-quality genome assembly provides a critical resource for future molecular breeding and genetic improvement of Tibetan sheep.