Genome-wide comparative analyses for selection signatures indicate candidate genes for between-breed variability in copper accretion in sheep

The problem of copper (Cu) intoxication and deficiency continues to impact economic gains and animal welfare in sheep husbandry. This study investigated the ovine genome for regions and potential genes under selection for Cu accretion between sheep breeds. For this, we compared ovine single nucleotide polymorphism (SNP) data of three Cu-susceptible breeds with three Cu-tolerant breeds. After merging SNP data of breeds and removal of related individuals, a total of 229 sheep and 45 640 autosomal SNPs were left. Then, we selected 14 individuals per breed into two datasets (datasets 1 and 2) for analysis of selection signatures using the Fixation index, cross-population extended haplotype homozygosity and haplotype-based FLK methods. Selection regions shared by both datasets detected by at least two methods revealed regions on OAR 4, 8 and 11 containing 54 candidate genes under selection for Cu accretion. Enrichment analysis revealed that 19 gene ontologies and 1 enriched Kyoto encyclopaedia of genes and genomes pathway terms were associated with the candidate genes under selection. Genes such as TP53, TNFSF13, TNFSF12, ALOX15, ALOX12, EIF5A and PREP are associated with the regulation of Cu homeostasis, programmed cell death or inflammatory response. We also found an enrichment of arachidonate 15-lipoxygenase activity, arachidonate 12-lipoxygenase activity and ferroptosis that influence cellular inflammation and cell death. These results shed light on ovine genomic regions under selection for Cu accretion and provide information on candidate genes for further studies on breed differences in ovine Cu accretion.