The Mexican turkey population is considered to be the descendant of the original domesticated wild turkey and it is distinct from hybrid strains obtained by the intense artificial selection activity that has occurred during the last 40 years. In this study 30 Mexican turkeys were genomically compared to 38 commercial hybrids using 327,342 SNP markers in order to elucidate the differences in genome variability resulting from different types of selection, i.e., only adaptive for Mexican turkey, and strongly directional for hybrids. Runs of homozygosity (ROH) were detected and the two inbreeding coefficients (F and FROH) based on genomic information were calculated. Principal component and admixture analyses revealed two different clusters for Mexican turkeys (MEX_cl_1 and MEX_cl_2) showing genetic differentiation from hybrids (HYB) (FST equal 0.168 and 0.167, respectively). A total of 3602 ROH were found in the genome of the all turkeys populations. ROH resulted mainly short in length and the ROH_island identified in HYB (n = 9), MEX_cl_1 (n = 1), and MEX_cl_2 (n = 2) include annotated genes related to production traits: abdominal fat (percentage and weight) and egg characteristics (egg shell color and yolk weight). F and FROH resulted correlated to each other only for Mexican populations. Mexican turkey genomic variability allows us to separate the birds into two subgroups according to the geographical origin of samples, while the genomic homogeneity of hybrid birds reflected the strong directional selection occurring in this population.
Keywords: ROH; SNPs; autochthonous populations; genome variability; inbreeding; run of homozygosity; selected breeds; turkey.