Runs of homozygosity and gene identification in Pelibuey sheep using genomic data

Abstract

The runs of homozygosity (ROHs), the inbreeding coefficient, and the effective population size (Ne) in Pelibuey sheep were analyzed in 24 Pelibuey ewes from two lambs at parturition and 24 ewes that gave birth to a single lamb using the Ilumina OvineSNP50 BeadChip. The Ne decreased from 535 to 192 in the first ten generations. A total of 2194 ROHs were identified on the basis of single nucleotide polymorphisms (SNPs), were identified in the prolific group and 2185 SNPs in ROH in the non-prolific group. The distribution of the lengths of the ROH, considering both groups, were found to be: 4065 less than 6 Mb, 213 between 6 and 12 Mb, 72 between 12 and 24 Mb, twenty between 24 and 48 Mb and 8 greater than 48 Mb. In prolific sheep, the ROH associated with prolificacy were identified near the LINGO2, FLRT2, ADGRB3 genes, related to “positive regulation of synapse assembly”; and the DGKG, DGKE, DGKB and DGKI genes, related to “protein kinase C-activating G-protein coupled receptor signaling pathway”. The present work present genes that can function as signal mediators or have activity in embryonic development, which is relevant to the economic activity of this species.