Application of a non-invasive tool and identification of genetic markers in swine to enhance ham quality: a preliminary study.

A genome-wide association study (GWAS) was performed as a preliminary step to identify regions potentially related to ham quality traits. In this research, genomic information was obtained from 238 commercial hybrid pigs utilising the GeneSeek® Genomic Profiler genome-wide porcine genotyping array. Carcasses were tested for hot weight, the thickness of backfat and loin, and lean meat percentage. The corresponding fresh hams were assayed for weight and ultimate pH; the activities of Cathepsin B and Ferrochelatase of Semimembranosus muscle were determined through fluorimetric methods. The lean meat percentage of fresh ham (LMPH), salt absorbed after first (SALT1) and overall salting stages (SALT) were estimated online by the Ham Inspector™ apparatus. Hams were processed in compliance with the procedures established for Protected Designation of Origin Parma ham, and ham weight losses were measured at the main processing stages. Hot carcass weights showed a significant negative correlation with their lean meat percentage and LMPH, while LMPH was correlated positively with carcass lean meat, SALT1, SALT, and weight losses. The GWAS detected genome-wide association for 12 single nucleotide polymorphisms with Ferrochelatase activity. The results obtained in this preliminary study were achieved by combining innovative and non-destructive technologies for screening hams under processing, measures of enzymatic muscle properties relevant to dry-cured ham quality, and genomic information obtained through a GWAS. Additional studies carried out in a larger number of pigs have been planned to investigate the effect of gene variants of Ferrochelatase activity in the dry-cured ham's quality with main reference to colour development and to confirm the GWAS results obtained in this study.

Genome-wide association study identifies markers associated with carcass and meat quality traits in Italian Large White pigs.

A GWAS was performed using the genotypes obtained by PorcineSNP60 v2 BeadChip and 11 phenotypic traits (carcass lean meat percentage; backfat thickness; Longissimus thoracis muscle thickness; lightness; backfat thickness measured with caliper at the midline; meat pH measured at about 1 h post mortem and 24 h post mortem; CIE L*, a* and b* color parameters; and water-holding capacity). Three markers were associated with three of the phenotypic traits considered: M1GA0009592 (SSC7) with backfat thickness and lean meat content, DIAS0002910 (SSC6) and ALGA0109856 (SSC6) with water-holding capacity. The marker M1GA0009592, associated with backfat thickness, lies in a QTL region near the gene JARID2, which is a transcription factor also involved in the regulation of adipose-derived stem cell pluripotency. The results seem to indicate a possible role of these genomic regions in the regulation of pig carcass fatness (i.e. backfat at last rib) and water-holding capacity.

Genome-wide association studies for the number of teats and teat asymmetry patterns in Large White pigs.

The number of teats is a morphological trait that influences the mothering ability of the sows and thus their reproduction performances. In this study, we carried out GWASs for the total number of teats and other 12 related parameters in 821 Italian Large White heavy pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip array. For four investigated parameters (total number of teats, the number of teats of the left line, the number of teats of the right line and the maximum number of teats comparing the two sides), significant markers were identified on SSC7, in the region of the vertnin (VRTN) gene. Significant markers for the numbers of posterior teats and the absolute difference between anterior and posterior teat numbers were consistently identified on SSC6. The most significant SNP for these parameters was an intron variant in the TOX high mobility group box family member 3 (TOX3) gene. For the other four parameters (absolute difference between the two sides; anterior teats; the ratio between the posterior and the anterior number of teats; and the absence or the presence of extra teats) only suggestively significant markers were identified on several other chromosomes. This study further supported the role of the VRTN gene region in affecting the recorded variability of the number of teats in the Italian Large White pig population and identified a genomic region potentially affecting the biological mechanisms controlling the developmental programme of morphological features in pigs.

Genome-wide association study identifies quantitative trait loci regions involved in muscle acidic profile in Large White heavy pigs.

The widespread use of genome-wide association studies resulted in the discovery of genomic regions associated with fatty acid (FA) composition in different porcine tissues, but little information exists about the genes involved in FA composition of meat obtained from heavy pigs selected for the production of Italian dry-cured hams. To this objective, we genotyped with a single nucleotide polymorphism (SNP) panel 795 Italian Large White heavy pigs to identify the markers and genomic regions associated with Semimembranosus muscle FA profile. Heritability estimates for intramuscular fat FA profile were of low-to-moderate magnitude, suggesting that these traits may be improved with genomic selection. On the whole, 45 SNPs were significantly associated with 14 FAs, and 4 of them (ALGA008109, ALGA0081097, CASI0010164 and SIRI0000267) were associated with more than 1 FA. The palmitoleic : palmitic and oleic : stearic ratios displayed the highest number of significant markers and the most significant associations (Bonferroni adjusted P < 5.00E-07). Of particular interest, the palmitoleic : palmitic ratio was strongly associated with markers located at 111 to 114 Mb on chromosome 14, in the same chromosomal region where Stearoyl-CoA desaturase Δ9 (SCD) gene is located. Several significant chromosomal regions were found; some of them harbour key genes playing pivotal roles in FA desaturation and elongation, such as SCD and some members of the Elongation of Very Long-Chain FA (ELOVL) gene family. The results suggest that the identification of causal mutations in these regions may provide a set of markers useful for selection schemes aimed at improving FA composition in pork products.

Genome-wide association study identifies markers associated with meat ultimate pH in Duroc pigs.

In this study we aimed to identify genomic regions associated with muscle pH, meat colour and water-holding capacity in a population of 280 Italian Duroc pigs genotyped by the Illumina PorcineSNP60 v2 Genotyping BeadChip. After quality control, the remaining 32 597 SNPs and 278 subjects were used to perform a genome-wide association study with the genabel package, using a kinship matrix in a model with the effects of sex, age and slaughter day. Bonferroni correction was applied, and the significant markers and regions were then further investigated to identify the nearest genes and the linkage disequilibrium (LD) between markers. Four markers (ASGA0082344, ASGA0095635, DBWU0000985 and CASI0005117) were significantly associated with ultimate pH (pHu ); no significant association was detected for the other traits. The four significant variants, located from 16.841 to 17.643 Mb on chromosome 3, were found within or close to the sequences of the sulfatase modifying factor 2 (SUMF2), lysine acetyltransferase 8 (KAT8), serine protease 8 (PRSS8) and phosphorylase kinase catalytic subunit gamma 2 (PHKG2) genes. The four associated markers lie in two LD blocks, suggesting that the observed effect is related to mutations located in two regions: the first one where SUMF2 is mapped and the second one where genes KAT8, PRSS8 and PHKG2 are located.

Genetic parameters of backfat fatty acids and carcass traits in Large White pigs.

Subcutaneous fat thickness and fatty acid composition (FAC) play an important role on seasoning loss and organoleptic characteristics of seasoned hams. Dry-cured ham industry prefers meats with low contents of polyunsaturated fatty acids (PUFA) because these negatively affect fat firmness and ham quality, whereas consumers require higher contents in those fatty acids (FA) for their positive effect on human health. A population of 950 Italian Large White pigs from the Italian National Sib Test Selection Programme was investigated with the aim to estimate heritabilities, genetic and phenotypic correlations of backfat FAC, Semimembranosus muscle intramuscular fat (IMF) content and other carcass traits. The pigs were reared in controlled environmental condition at the same central testing station and were slaughtered at reaching 150 kg live weight. Backfat samples were collected to analyze FAC by gas chromatography. Carcass traits showed heritability levels from 0.087 for estimated carcass lean percentage to 0.361 for hot carcass weight. Heritability values of FA classes were low-to-moderate, all in the range 0.245 for n-3 PUFA to 0.264 for monounsaturated FA (MUFA). Polyunsaturated fatty acids showed a significant genetic correlation with loin thickness (0.128), backfat thickness (-0.124 for backfat measured by Fat-O-Meat'er and -0.175 for backfat measured by calibre) and IMF (-0.102). Obviously, C18:2(n-6) shows similar genetic correlations with the same traits (0.211 with loin thickness, -0.206 with backfat measured by Fat-O-Meat'er, -0.291 with backfat measured by calibre and -0.171 with IMF). Monounsaturated FA, except with the backfat measured by calibre (0.068; P<0.01), do not show genetic correlations with carcass characteristics, whereas a negative genetic correlation was found between MUFA and saturated FA (SFA; -0.339; P<0.001). These results suggest that MUFA/SFA ratio could be increased without interfering with carcass traits. The level of genetic correlations between FA and carcass traits should be taken into account in dealing with the development of selection schemes addressed to modify carcass composition and/or backfat FAC.

Identification of differentially expressed small RNAs and prediction of target genes in Italian Large White pigs with divergent backfat deposition.

The identification of the molecular mechanisms regulating pathways associated with the potential for fat deposition in pigs can lead to the detection of key genes and markers for the genetic improvement of fat traits. Interactions of microRNAs (miRNAs) with target RNAs regulate gene expression and modulate pathway activation in cells and tissues. In pigs, miRNA discovery is far from saturation, and the knowledge of miRNA expression in backfat tissue and particularly of the impact of miRNA variations is still fragmentary. Using RNA-seq, we characterized the small RNA (sRNA) expression profiles in Italian Large White pig backfat tissue. Comparing two groups of pigs divergent for backfat deposition, we detected 31 significant differentially expressed (DE) sRNAs: 14 up-regulated (including ssc-miR-132, ssc-miR-146b, ssc-miR-221-5p, ssc-miR-365-5p and the moRNA ssc-moR-21-5p) and 17 down-regulated (including ssc-miR-136, ssc-miR-195, ssc-miR-199a-5p and ssc-miR-335). To understand the biological impact of the observed miRNA expression variations, we used the expression correlation of DE miRNA target transcripts expressed in the same samples to define a regulatory network of 193 interactions between DE miRNAs and 40 DE target transcripts showing opposite expression profiles and being involved in specific pathways. Several miRNAs and mRNAs in the network were found to be expressed from backfat-related pig QTL. These results are informative for the complex mechanisms influencing fat traits, shed light on a new aspect of the genetic regulation of fat deposition in pigs and facilitate the prospective implementation of innovative strategies of pig genetic improvement based on genomic markers.

Investigation of the Perilipin 5 gene expression and association study of its sequence polymorphism with meat and carcass quality traits in different pig breeds.

The growing request for healthier fatty acid composition of animal products is raising the necessity of a deeper knowledge of the main factors controlling fatty acids storage in muscle and backfat. Perilipin (PLIN) 5, and the whole Perilipin family, seem to play a crucial role in the regulation of lipids deposition as code for proteins coating intracellular lipid droplets surface. Nevertheless, the knowledge of these genes in pig is still incomplete. The present research was aimed at investigating in different pig breeds the PLIN5 gene, analysing its expression level and the associations of the variability in its downstream gene region with pork meat and carcass quality traits. Moreover, the PLIN5 protein localisation in porcine skeletal muscle was investigated through immunofluorescence, resulting to be widespread in Semimembranosus muscle (SM) myofibers. The identified single nucleotide polymorphism (SNP) rs327694326 (NC_010444.4:g.74314701T>C) located in PLIN5 downstream region was analysed in different pig populations, represented by 512 Italian Large White (ILW) pigs, 300 Italian Duroc (IDU) samples, 100 Italian Landrace, 100 Pietrain and 60 pigs belonging to three Italian native breeds (20 samples of Cinta Senese, 20 Calabrese and 20 Casertana pigs). The C allele was found in ILW, IDU and Pietrain pigs. In ILW pigs this SNP showed results indicating a possible association with oleic and cis-vaccenic fatty acid contents in backfat tissue. Furthermore, as PLINs are known to regulate lipases activity, we tested if the rs327694326 SNP was associated with differences in Hormone-sensitive lipase (LIPE) gene expression levels. In SM of ILW pigs, PLIN5 C allele was associated with significantly lower LIPE transcription levels than T allele (P=0.02 for Student's t test of TT v. CT samples, P<0.0001 for TT v. CC pigs), whereas in IDU breed no significant difference was noticed in LIPE transcription levels between TT and TC animals. The results may suggest that variations in the PLIN5 sequence may be linked to LIPE expression through a still poorly known regulative molecular process.

Transcriptional profiling of subcutaneous adipose tissue in Italian Large White pigs divergent for backfat thickness.

Fat deposition is a widely studied trait in pigs because of its implications with animal growth efficiency, technological and nutritional characteristics of meat products, but the global framework of the biological and molecular processes regulating fat deposition in pigs is still incomplete. This study describes the backfat tissue transcription profile in Italian Large White pigs and reports genes differentially expressed between fat and lean animals according to RNA-seq data. The backfat transcription profile was characterised by the expression of 23 483 genes, of which 54.1% were represented by known genes. Of 63 418 expressed transcripts, about 80% were non-previously annotated isoforms. By comparing the expression level of fat vs. lean pigs, we detected 86 robust differentially expressed transcripts, 72 more highly expressed (e.g. ACP5, BCL2A1, CCR1, CD163, CD1A, EGR2, ENPP1, GPNMB, INHBB, LYZ, MSR1, OLR1, PIK3AP1, PLIN2, SPP1, SLC11A1, STC1) and 14 lower expressed (e.g. ADSSL1, CDO1, DNAJB1, HSPA1A, HSPA1B, HSPA2, HSPB8, IGFBP5, OLFML3) in fat pigs. The main functional categories enriched in differentially expressed genes were immune system process, response to stimulus, cell activation and skeletal system development, for the overexpressed genes, and unfolded protein binding and stress response, for the underexpressed genes, which included five heat shock proteins. Adipose tissue alterations and impaired stress response are linked to inflammation and, in turn, to adipose tissue secretory activity, similar to what is observed in human obesity. Our results provide the opportunity to identify biomarkers of carcass fat traits to improve the pig production chain and to identify genetic factors that regulate the observed differential expression.

Genome-wide study on intramuscular fat in Italian Large White pig breed using the PorcineSNP60 BeadChip.

Genome-wide association study results are presented for intramuscular fat in Italian Large White pig breed. A total of 886 individuals were genotyped with PorcineSNP60 BeadChip. After quality control performed with plink software and in R environment, 49 208 markers remained for the association analysis. The genome-wide association studies was conducted using linear mixed model implemented in GenABEL. We detected seven new SNPs of genes till now not found associated to intramuscular fat (IMF). Three markers map in a wide intergenic region rich of QTL linked to fat traits, one map 388 kb upstream the gene SDK1, one map inside PPP3CA gene, one inside SCPEP1 gene and the last is not mapped in the porcine genome yet. Associations here presented indicate a moderate effect of these genes on IMF. In particular, PPP3CA, that is involved in the oxidative metabolism of skeletal muscle, could be considerated as an interesting candidate gene for IMF content in pigs. However, further studies are needed to clarify the role of these genes on the physiological processes involved in IMF regulation. These results may be useful to control this trait that is important in terms of nutritional, technological and organoleptic characteristics of fresh meat and processed products.

Relationship between perilipin genes polymorphisms and growth, carcass and meat quality traits in pigs.

The perilipins (PLIN) belong to a family of structural proteins that play a role regulating intracellular lipid storage and mobilization. Here, PLIN1 and PLIN2 have been evaluated as candidate genes for growth, carcass and meat quality traits in pigs. A sample of 607 Duroc pigs were genotyped for two single-nucleotide polymorphisms, one in intron 2 of the PLIN1 gene (JN860199:g.173G>A) and the other at the 3' untranslated region of the PLIN2 gene (GU461317:g.98G>A). Using a Bayesian approach, we have been able to find evidence of additive, dominant and epistatic associations of the PLIN1 and PLIN2 polymorphisms with early growth rate and carcass length. However, the major effects were produced by the dominant A allele at the PLIN2 polymorphism, which also affected the carcass lean weight. Thus, pigs carrying an additional copy of the A allele at the g.98G>A PLIN2 polymorphism had a probability of at least 98% of producing carcasses with heavier lean weight (+0.41 kg) and ham weight (+0.10 kg). The results obtained indicate that the PLIN2 polymorphism could be a useful marker for lean growth. In particular, it may help to reduce the undesired negative correlated response in lean weight to selection for increased intramuscular fat content, a common scenario in some Duroc lines involved in the production of high quality pork products.

Genetic analysis of seven Italian horse breeds based on mitochondrial DNA D-loop variation.

To understand the origin and genetic diversity of Italian horses, mitochondrial DNA D-loop sequences were generated for 163 horses from seven breeds. Sequence analysis of a 480-bp segment revealed a total of 84 haplotypes with 57 polymorphic sites, indicating multiple maternal origins and high genetic diversity. Comparison of the haplotypes with the equine mtDNA haplotype/haplogroup nomenclature showed a haplogroup distribution in the Italian breeds more similar to that found in the Middle East breeds than in the European breeds, probably due to the economic and cultural relationship with the Middle East in the past centuries.

Association mapping of quantitative trait loci for carcass and meat quality traits at the central part of chromosome 2 in Italian Large White pigs.

Association mapping of the central part of porcine chromosome 2 harboring QTLs for carcass and meat quality traits was performed with 17 gene-tagged SNPs located between 44.0 and 77.5 Mb on a physical map (Sscrofa10.2) in Italian Large White pigs. For the analyzed animals records of estimated breeding values for average daily gain, back fat thickness, lean cuts, ham weight, feed conversion ratio, pH1, pHu, CIE L*, CIE a*, CIE b* and drip loss were available. A significant QTL for fat deposition (adjusted P=0.0081) and pH1 (adjusted P=0.0972) to MYOD1 at position 44.4 Mb and a QTL for growth and meatiness (adjusted P=0.0238-0.0601) to UBL5 at position 68.9 Mb were mapped. These results from association mapping are much more accurate than those from linkage mapping and facilitate further search for position candidate genes and causative mutations needed for application of markers through marker assisted selection.

Detecting population structure and recent demographic history in endangered livestock breeds: the case of the Italian autochthonous donkeys.

Since its domestication, about 5000 years ago, the donkey (Equus asinus) has been extensively used as a work or draft animal in agricultural activities and for the transportation of people and goods. In the last century, technology improvement and growing mechanization strongly affected agriculture and the management and use of this livestock species in the industrialized countries. Nowadays, the use of donkeys for work or transport has almost disappeared, together with the need for mules or hinny breeding. During the last five decades, Italian autochthonous donkey populations suffered from a severe reduction in population size, which led to the extinction of several breeds. At present, eight breeds remain, all classified by FAO as critically endangered or endangered: Asinara, Pantesco, Grigio Siciliano, Romagnolo, Amiatino, Sardo Grigio, Martina Franca, and Ragusano. To evaluate the extant genetic variability of Italian donkeys, we typed 16 microsatellite loci in 258 individuals from these breeds. The results highlighted moderate levels of inbreeding ( F (IS) = 0.127) and a significant partition of genetic variation into breeds, as suggested by fixation index ( F (ST) = 0.109) and analysis of molecular variance (10.86% of total variation assigned to the between-breeds level) analyses. This was confirmed by a Bayesian clustering procedure that also highlighted a further partitioning at lower hierarchical levels corresponding to the farms of origin. This evidence suggests that an effective management strategy for Italian donkey populations should focus on breeds as conservation units. However, this requires a synergic management strategy at the farm level to maintain diversity and avoid inbreeding.

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design.

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.

Analysis of MC4R polymorphism in Italian Large White and Italian Duroc pigs: association with carcass traits.

The melanocortin-4 receptor (MC4R) gene codes for a G protein transmembrane receptor playing an important role in energy homeostasis control. In pig a single nucleotide polymorphism c.1426G>A has been identified and associated to average daily gain, feed intake and fatness traits but a lack of agreement on the effects of the gene on carcass traits in different breeds comes out from many studies. In the present study the c.1426G>A polymorphism is analysed in two Italian pig breeds, Large White and Duroc to study the association of the MC4R gene with some carcass traits. The results show that the c.1426G>A polymorphism affects daily gain, feed conversion ratio and ham weight in both breeds, lean cuts in the Italian Duroc and backfat thickness in the Italian Large White. The presence of MC4R mRNA transcript in different porcine tissues was analysed.

Perilipin 1 and perilipin 2 protein localization and gene expression study in skeletal muscles of European cross-breed pigs with different intramuscular fat contents.

This study investigated the lipid droplet coat proteins perilipin 1 (PLIN1) and perilipin 2 (PLIN2) localization in pig skeletal muscle and their relationship with intramuscular fat (IMF) content. PLIN1 and PLIN2 proteins were immunostained in semimembranosus muscle cross-sections from two groups of samples divergent for IMF and the gene expression was quantified. PLIN1 localized in the periphery of intramuscular adipocytes, whereas PLIN2 localized within myofibers with high lipid content. The high IMF group showed higher total cross-sectional area of PLIN1-stained adipocytes compared with the low IMF group (P<0.05), while the cross-sectional area and percentage of PLIN2-positive myofibers did not differ between IMF-divergent groups. This suggested that IMF content is mainly determined by extra-myocellular lipids. At mRNA level, PLIN2 expression was higher in high IMF muscles (P<0.05). The results indicate for the first time that in pig muscle PLIN1 and PLIN2 proteins are localized in correspondence with extra and intra-myocellular lipids, respectively.

Isolation of porcine expressed sequence tags for the construction of a first genomic transcript map of the skeletal muscle in pig.

To identify genes with effects on meat quality and production traits we developed an adult porcine skeletal muscle cDNA library. After pre-screening this library with seven genes highly expressed in skeletal muscle, 385 non-hybridizing clones were sequenced from both ends to yield 510 expressed sequence tags (ESTs). Together with those ESTs previously generated from this library, we have produced 701 porcine skeletal muscle ESTs. These ESTs were grouped into 306 different cDNA species and compared with the human skeletal muscle transcriptional profiles obtained from different databases. Furthermore we mapped 107 of these cDNAs using a somatic cell hybrid panel with genes mapping over all the autosomes (except on chromosome 11) and on chromosome X. The mapping of these cDNAs contributed to the construction of a first genomic transcript map of the skeletal muscle tissue in pig.