Genome-wide association studies for seven production traits highlight genomic regions useful to dissect dry-cured ham quality and production traits in Duroc heavy pigs

Protected designation of origin dry-cured hams are obtained from heavy pigs (slaughtered at about 160 kg of live weight). A specific breeding program designed to improve meat quality for this production has included as key traits the level of intermuscular fat between the leg muscles and ham weight loss during the seasoning period together with a balance between fat and lean cuts. In this study we carried out genome-wide association studies for seven traits used in the genetic merit of Italian Duroc heavy pigs, five related to meat and carcass quality traits (visible intermuscular fat, ham weight loss at first salting, backfat thickness, ham weight and lean cuts), and two related to performance and efficiency traits (average daily gain and feed : gain ratio). A total of 573 performance-tested pigs were genotyped with the Illumina PorcineSNP60 BeadChip and genome-wide association analyses were carried out using the Bayes B approach with the 1 Mb window option of GenSel and random residuals for each of the seven traits. Detected windows were supported by independent single nucleotide polymorphism analyses with a linear mixed model (LMM) approach on the same animals for the same traits. A total of 30 windows identifying different quantitative trait loci (QTL) were detected and among those, 27 were confirmed by LMM in one of these traits. Among the confirmed windows, three QTL were reported for visible intermuscular fat, seven for ham weight loss at first salting and five and four for backfat thickness and lean cut, respectively. A total of eight QTL were detected for the other production traits. No overlapping QTL were reported except for one window on porcine chromosome 10 between lean cuts and ham weight that contained the CACNB2 gene that has been already associated with loin marbling score in other Duroc pigs. Several regions contained genes that have been already associated with production traits in other pig breeds, including Duroc lines, related to fat deposition or muscle structure. This work reports, for the first time, genome-wide association study results for several traits in Italian Duroc heavy pigs. These results will be useful to dissect the genetic basis for dry-cured ham production traits that determine the total genetic merit index of Italian Duroc pigs.     

A genome‐wide association study for a proxy of intermuscular fat level in the Italian Large White breed identifies genomic regions affecting an important quality parameter for dry‐cured hams

Intermuscular fat content in protected designations of origin dry‐cured hams is a very important meat quality trait that affects the acceptability of the product by the consumers. An excess in intermuscular fat (defined as the level of fat deposition between leg muscles) is a defect that depreciates the final product. In this study we carried out a genome‐wide association study for visible intermuscular fat (VIF) of hams in the Italian Large White pig breed. This trait was evaluated on the exposed muscles of green legs in 1122 performance‐tested gilts by trained personnel, according to a classification scale useful for routine and cheap evaluation. All animals were genotyped with the Illumina PorcineSNP60 BeadChip. The genome‐wide association study identified three QTL regions on porcine chromosome 1 (SSC1; accounting for ~79% of the SNPs below the 5.0E?04 threshold) and SSC2, two on SSC7 and one each on SSC3, SSC6, SSC9, SSC11, SSC13, SSC15, SSC16 and SSC17. The most significant SNP (ALGA0004143 on SSC1 at 77.3 Mb; P_FDR < 0.05), included in the largest QTL region which spanned about 6.8 Mb on SSC1, is located within the glutamate ionotropic receptor kainate type subunit 2 (GRIK2) gene. Functional annotation of all genes included in QTL regions for VIF suggested that intermuscular fat in the Italian Large White breed is a complex trait apparently influenced by complex biological mechanisms also involving obesity‐related processes. These QTL target mainly chromosome regions different from those affecting subcutaneous and intramuscular fat deposition.     

Factors in pig production that impact the quality of dry-cured ham: a review

This study reviews the factors of pig production that impact the quality of dry-cured ham. When processing is standardized, the quality of the final dry-cured product is primarily determined by the quality of the meat before curing (green ham). This has been defined as the aptitude for seasoning and is determined by the green ham weight, adipose tissue quantity and quality, meat physico-chemical properties and the absence of visual defects. Various ante-mortem factors including pig age and weight, genetic type, diet, feeding strategy and slaughter conditions determine green ham properties such as the dynamics of water loss, salt intake and, as a consequence, proteolysis and lipolysis. Muscle conditions (pH, salt concentration, water content and availability, temperature) influence enzymatic activity and development of characteristic texture and flavor. Generally, hams of older and heavier pigs present better seasoning aptitude because of higher adiposity. Adiposity is also positively correlated with fat saturation, which is desired to avoid rancidity and oiliness. The fatty acid profile of tissue lipids can be manipulated by diet composition. Feeding strategy affects tissue accretion and protein turnover, thus directly impacting proteolysis. With respect to the impact of pig genotype on dry-cured ham quality, local breeds are generally considered more suitable for producing quality dry hams; however, the majority of dry-cured hams on the market today are from modern pig breeds raised in conventional systems, providing lean hams. The importance of all these factors of pig production is discussed and synthesized, with an emphasis on the main difficulties encountered in dry-cured ham production.     

Genome‐wide association analysis for growth,muscularity and meat quality in Piétrain pigs

Improvement in growth and meat quality is one of the main objectives in sire line pig breeding programmes. Mapping quantitative trait loci for these traits using experimental crosses and a linkage‐based approach has been performed frequently in the past. The Piétrain breed often was involved as a founder breed to establish the experimental crosses. This breed was selected for muscularity and leanness but shows relatively poor meat quality. It is frequently used as a sire line breed. With the advent of genome‐wide and dense SNP chips in pig genomic research, it is possible to also conduct genome‐wide association studies within the Piétrain breed. In this study, around 500 progeny‐tested sires were genotyped with 60k SNPs. Data filtering showed that around 48k SNPs were useable in this sample. These SNPs were used to conduct a genome‐wide association study for growth, muscularity and meat quality traits. Because it is known that a mutation in the RYR1 gene located on chromosome 6 shows a major effect on meat quality, this mutation was included in the models. Single‐marker and multimarker association analyses were performed. The results revealed between zero and eight significant associations per trait with P < 5 × 10^?5. Of special interest are SNPs located on SSC6, SSC10 and SSC15.     

Prediction of dry-cured ham weight loss and prospects of use in a pig breeding program

Large ham weight losses (WL) in dry-curing are undesired as they lead to a loss of marketable product and penalise the quality of the dry-cured ham. The availability of early predictions of WL may ease the adaptation of the dry-curing process to the characteristics of the thighs and increase the effectiveness of selective breeding in enhancing WL. Aims of this study were (i) to develop Bayesian and Random Forests (RFs) regression models for the prediction of ham WL during dry-curing using on-site infrared spectra of raw ham subcutaneous fat, carcass and raw ham traits as predictors and (ii) to estimate genetic parameters for WL and their predictions (P-WL). Visible-near infrared spectra were collected on the transversal section of the subcutaneous fat of raw hams. Carcass traits were carcass weight, carcass backfat depth, lean meat content and weight of raw hams. Raw ham traits included measures of ham subcutaneous fat depth and linear scores for round shape, subcutaneous fat thickness and marbling of the visible muscles of the thigh. Measures of WL were available for 1672 hams. The best prediction accuracies were those of a Bayesian regression model including the average spectrum, carcass and raw ham traits, with R² values in validation of 0.46, 0.55 and 0.62, for WL at end of salting (23 days), resting (90 days) and curing (12 months), respectively. When WL at salting was used as an additional predictor of total WL, the R² in validation was 0.67. Bayesian regressions were more accurate than RFs models in predicting all the investigated traits. Restricted maximum likelihood (REML) estimates of genetic parameters for WL and P-WL at the end of curing were estimated through a bivariate animal model including 1672 measures of WL and 8819 P-WL records. Results evidenced that the traits are heritable (h² ± SE was 0.27 ± 0.04 for WL and 0.39 ± 0.04 for P-WL), and the additive genetic correlation is positive and high (r_a = 0.88 ± 0.03). Prediction accuracy of ham WL is high enough to envisage a future use of prediction models in identifying batches of hams requiring an adaptation of the processing conditions to optimise results of the manufacturing process. The positive and high genetic correlation detected between WL and P-WL at the end of dry-curing, as well as the estimated heritability for P-WL, suggests that P-WL can be successfully used as an indicator trait of the measured WL in pig breeding programs.     

New SNP of the porcine Perilipin 2 (<Emphasis Type="Italic">PLIN2</Emphasis>) gene,association with carcass traits and expression analysis in skeletal muscle

PLIN2 (perilipin 2) is a cytosolic protein that promotes the formation and stabilization of the intracellular lipid droplets, organelles involved in the storage of lipid depots. Porcine PLIN2 gene represents a biological and positional candidate for fat deposition, a polygenic trait that affects carcass and meat quality. The aim of the present study was to screen PLIN2 gene for polymorphisms, to evaluate the association with carcass quality traits, and to investigate the gene expression in skeletal muscle. Six new single nucleotide polymorphisms (SNP) were detected by sequencing 32 samples from five pig breeds (Italian Large White, Italian Duroc, Italian Landrace, Belgian Landrace, Pietrain). Two SNP localized in introns, two in the 3′-untranslated region (UTR), and two missense SNP were found in exons. A 3′-UTR mutation (GU461317:g.98G>A), genotyped in 290 Italian Duroc pigs by High Resolution Melting, resulted significantly associated (P < 0.01) with average daily gain, feed conversion ratio, lean cuts and hams weight estimated breeding values. PLIN2 gene expression analysis in skeletal muscle of Italian Large White and Italian Duroc pigs divergent for backfat thickness and visible intermuscular fat showed a trend of higher expression level in pigs with higher intermuscular fat. These results suggest that PLIN2 can be a marker for carcass quality in pigs. Further investigation at both gene and protein level could elucidate its role on fat deposition.     

Genome‐wide association study for reproductive traits in a Large White pig population

Using the PorcineSNP80 BeadChip, we performed a genome‐wide association study for seven reproductive traits, including total number born, number born alive, litter birth weight, average birth weight, gestation length, age at first service and age at first farrowing, in a population of 1207 Large White pigs. In total, we detected 12 genome‐wide significant and 41 suggestive significant SNPs associated with six reproductive traits. The proportion of phenotypic variance explained by all significant SNPs for each trait ranged from 4.46% (number born alive) to 11.49% (gestation length). Among them, 29 significant SNPs were located within known QTL regions for swine reproductive traits, such as corpus luteum number, stillborn number and litter size, of which one QTL region associated with litter size contained the ALGA0098819 SNP for total number born. Subsequently, we found that 376 functional genes contained or were near these significant SNPs. Of these, 14 genes—BHLHA15, OCM2, IL1B2, GCK, SMAD2, HABP2, PAQR5, GRB10, PRELID2, DMKN, GPI, GPIHBP1, ADCY2 and ACVR2B—were considered important candidates for swine reproductive traits based on their critical roles in embryonic development, energy metabolism and growth development. Our findings contribute to the understanding of the genetic mechanisms for reproductive traits and could have a positive effect on pig breeding programs.     

Associations between single nucleotide polymorphisms in 33 candidate genes and meat quality traits in commercial pigs

This study aimed to evaluate the effects of single nucleotide polymorphisms (SNPs) in candidate genes for meat quality using a custom 96‐SNP panel (Illumina Vera Code GoldenGate Assay) on 15 traits collected from 400 commercial pigs. Meat quality measurements included muscle pH, color (L*, a* and b*), drip loss, cooking loss, peak shear force and six sensory traits including appearance (outside and inside), tenderness, juiciness, flavor and overall liking as well as carcass weight and probe yield. Thirty‐five SNPs with minor allele frequencies > 0.10 remained for the multimarker association using the GLM procedure of sas 9.2. Results showed that 20 SNPs were significantly associated with at least one of the traits with either additive or dominance or both effects (P < 0.05). Among these significant SNPs, five of them in ADIPOQ, FTO, TNF, LEPR and AMPD1 had an effect on more than three traits simultaneously; those in MC4R, CAST, DGAT1 and MYF6 had an effect on two traits, while the others were associated with one trait. The results suggest that these markers could be incorporated into commercial pigs for marker‐assisted selection and breeding programs for carcass and meat quality trait improvement.     

Identification of four SNPs and association analysis with meat quality traits in the porcine <Emphasis Type="Italic">Pitx2c</Emphasis> gene

The association of the porcine Pitx2c gene with meat quality traits was investigated in the present study. A total of eight single nucleotide polymorphisms (SNPs) were found. Allele frequencies of four SNPs were further detected in four commercial breeds and eight Chinese indigenous breeds. Single SNP and meat quality associations were analyzed in a Yorkshire×Meishan F₂ population. The SNPs c.474C>T (P<0.01) and c.636C>T (P<0.05) showed a significant association with meat color (MCV1). The SNPs c.*37G>A and c.*47G>A were significantly associated with drip loss rate (DLR), water holding capacity (WHC) and meat color value (MCV1) consistently (P<0.05). Linkage disequilibrium (LD) analysis revealed that the adjacent SNPs were in LD. Two major haplotypes were identified, and association analysis between haplotype combinations and meat quality indicated that the presence of two copies of haplotype 1 -CCGG- may improve meat quality.     

RNA‐seq analysis reveals new candidate genes for drip loss in a Pietrain × Duroc × Landrace × Yorkshire population

Drip loss, one of the most important meat quality traits, is characterized by low heritability. To date, the genetic factors affecting the drip loss trait have not been clearly elucidated. The objective of this study was to identify critical candidate genes affecting drip loss. First, we generated a Pietrain × Duroc × Landrace × Yorkshire commercial pig population and obtained phenotypic values for the drip loss trait. Furthermore, we constructed two RNA libraries from pooled samples of longissimus dorsi muscles with the highest (H group) and lowest (L group) drip loss and identified the differentially expressed genes (DEGs) between these extreme phenotypes using RNA‐seq technology. In total, 25 883 genes were detected in the H and L group libraries, and none was specifically expressed in only one library. Comparative analysis of gene expression levels found that 150 genes were differentially expressed, of which 127 were upregulated and 23 were downregulated in the H group relative to the L group. In addition, 68 drip loss quantitative trait loci (QTL) overlapping with 63 DEGs were identified, and these QTL were distributed mainly on chromosomes 1, 2, 5 and 6. Interestingly, the triadin (TRDN) gene, which is involved in muscle contraction and fat deposition, and the myostatin (MSTN) gene, which has a role in muscle growth, were localized to more than two drip loss QTL, suggesting that both are critical candidate genes responsible for drip loss.     

Genome‐wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome‐wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F₂ progeny were analyzed for each blood‐related trait after filtering for quality control. Data were analyzed by the genome‐wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10^?6) on SSC3, 6, 8, 13 and 17 were identified for blood‐related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.     

Whole‐genome association study for fatty acid composition of oleic acid in Japanese Black cattle

Fatty acid composition, especially oleic acid (C18:1), plays an important role in the eating quality of meat in Japanese Black cattle. Therefore, the objective of this study was to identify loci associated with C18:1 in the intramuscular fat of the trapezius muscles in Japanese Black cattle using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. We also evaluated the relationship between C18:1 and three fatty acid synthesis genes, fatty acid synthase (FASN), stearoyl‐CoA desaturase and sterol regulatory element‐binding protein‐1. In this experiment, we applied a mixed model and Genomic Control approach using selective genotyping to perform a genome‐wide association study. A total of 160 animals (80 animals with higher values and 80 animals with lower values), selected from 3356 animals based on corrected phenotype, were genotyped using the Illumina BovineSNP50 BeadChip and three fatty acid synthesis genes, and the quality of these SNPs was assessed. In this study, a total of 38 955 SNPs, which included SNPs in the three fatty acid synthesis genes, were used, and the estimated inflation factor was 1.06. In the studied population, a total of 32 SNPs, including the FASN gene, had significant effects, and in particular 30 SNPs of all significant SNPs were located between 49 and 55 Mbp on chromosome 19. This study is one of the first genome‐wide association studies for fatty acid composition in a cattle population using the recently released Illumina BovineSNP50 BeadChip.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

Identification of genes related to intramuscular fat independent of backfat thickness in Duroc pigs using single‐step genome‐wide association

Intramuscular fat (IMF) is an important meat‐quality trait of pigs, which influences pork’s shearing force, hydraulics, tenderness and juicy flavor. However, to achieve a higher percentage of lean meat, pigs with lower backfat thickness (BF) are intensively selected for, which may lead to a reduction in pork quality. Therefore, the objective of this study was to locate loci that affect IMF without changing BF. A single‐step GWAS was performed on 950 Duroc pigs genotyped by a 50K SNP chip in order to detect genomic variants relevant to IMF and BF. The significant SNPs detected were afterwards divided into a BF subset (seven SNPs), an IMF subset (11 SNPs) and a subset of both traits (12 SNPs), according to their P‐value and LD. After SNP and QTL annotation, our results indicated that SSC1: 167938652, 166363826, 164829874 and 167171587 might be associated with IMF without changing BF. In the subset of both traits, we found that the combined effect of ALGA0006602 (SSC1: 159538854) and 12784636 (SSC1: 160773437) might improve the IMF without changing BF. Our gene annotation result showed that TLE3, ITGA11, SMAD6, PAQR5 and [RNF152^A/G × MC4R^A/A] genes might affect IMF independently of BF. We believe that the SNPs and genes identified in this study will be valuable for the future molecular breeding of IMF in Duroc pigs.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Molecular Characterization of Five Porcine Candidate Genes for Drip Loss in Pork

Drip loss is the loss of fluid from a piece of meat without mechanical force and represents an important meat quality trait. Previous work revealed a quantitative trait locus (QTL) for drip loss in pork in an experimental Duroc x Piétrain (DUPI) F2 family on SSC 5. Based on functional data indicating their possible involvement in water holding capacity and their expression in skeletal muscle, we selected five positional candidates (ACO2, ADSL, CBY1, KCNJ4, PLA2AG6) out of 130 predicted genes in the QTL interval for further analysis. We performed a mutation analysis of all coding exons and discovered 204 polymorphisms. We genotyped 39 single nucleotide polymorphisms (SNPs) in 192 Piétrain pigs with extreme drip loss phenotypes and detected a possible association with drip loss for one non-coding SNP in the ADSL gene (ss107793818, p_raw = 0.021). Correspondingly, ADSL diplotypes were associated with drip loss and pH1 of M. longissimus dorsi. However, after correction for multiple testing, none of the tested SNPs were significantly associated with drip loss. One possible explanation for these results is that one of the QTL-alleles from the experimental DUPI family may be fixed or nearly fixed in the tested Piétrain population.     

Sequence polymorphisms in porcine homologs of murine coat colour‐related genes

Herein, we report the variability among 57 porcine homologs of murine coat colour‐related genes. We identified single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within 44 expressed gene sequences by aligning eight pig complementary DNA (cDNA) samples. The sequence alignment revealed a total of 485 SNPs and 15 InDels. The polymorphisms were then validated by performing matrix‐assisted laser desorption/ionization time‐of‐flight mass spectrometry (MALDI‐TOF MS) with reference DNA samples obtained from 384 porcine individuals. Of the 384 individuals, three parents of the experimental F₂ family were included to detect polymorphisms between them for linkage mapping. We also genotyped previously reported polymorphisms of 12 genes, and one SNP each in three genes that were detected by performing a BLAST search of the Trace database. A total of 211 SNPs and three InDels were successfully genotyped from our porcine DNA panel. We detected SNPs in 33 of the 44 genes among the parents of an experimental F₂ family and then constructed a linkage map of the 33 genes for this family. The linkage assignment of each gene to the porcine chromosomes was consistent with the location of the BAC clone in the porcine genome and the corresponding gene sequence. We confirmed complete substitutions of EDNRB and MLPH in the Jinhua and Clawn miniature breeds, respectively. Furthermore, we identified polymorphic alleles exclusive to each pig group: 13 for Jinhua, two for Duroc, three for Meishan, four for the Japanese wild boar, one for the Clawn miniature pig and four for the Potbelly pig.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.