Next Article in Journal
Effect of Date Palm Pollen Supplementation on the Egg Production, Ovarian Follicles Development, Hematological Variables and Hormonal Profile of Laying Hens
Previous Article in Journal
Plasma Carboxyl-Metabolome Is Associated with Average Daily Gain Divergence in Beef Steers
Previous Article in Special Issue
Characterization of Subcutaneous Fat of Toscano Dry-Cured Ham and Identification of Processing Stage by Multivariate Analysis Approach Based on Volatile Profile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs

1
Department of Agricultural and Food Sciences (DISTAL), University of Bologna, Viale Fanin 46, I-40127 Bologna, Italy
2
Stazione Sperimentale per l’Industria delle Conserve Alimentari (SSICA), Viale Faustino Tanara 31/A, I-43121 Parma, Italy
3
Department of Agrifood, Environmental and Animal Science, University of Udine, Via delle Scienze 208, I-33100 Udine, Italy
*
Authors to whom correspondence should be addressed.
Animals 2021, 11(1), 68; https://doi.org/10.3390/ani11010068
Submission received: 15 December 2020 / Revised: 27 December 2020 / Accepted: 28 December 2020 / Published: 1 January 2021

Abstract

:

Simple Summary

In recent years, the meat industry is looking with increased interest at the implementation of non-invasive new tools for predicting the final quality of dry-cured hams and monitoring ham aging. The selection of raw meat and the control of the salting procedure to predict the quality of dry-cured hams are of primary importance for meat processors. The identification of genetic markers associated with ham traits and related to different aptitudes of the thighs towards salting phases and weight losses is a primary goal for the Italian pig production chain. This paper addresses the need for investigating the associations between genomic markers and ham traits obtained through the application of non-invasive technologies for monitoring hams before and during the salting process. To our knowledge, this is the first study investigating the markers and genes associated with ham traits obtained through the use of Ham InspectorTM apparatus.

Abstract

Protected Designation of Origin (PDO) dry-cured hams production is greatly dependent on raw meat quality. This study was performed to identify genetic markers associated with the quality of dry-cured ham. Carcass traits of 229 heavy pigs belonging to three commercial genetic lines were registered (weight, EUROP classification). Phenotypic traits (Semimembranosus muscle ultimate pH, ham weight and lean meat content, adsorbed salt) of the corresponding thighs, undergone PDO ham process in three different plants, were measured, using a fast and non-invasive technology. Green ham weight and lean meat percentage influenced the estimated salt content and the weight loss during salting, even if the processing plant greatly affected the variability of the measured ham traits. The genomic data were obtained with the GeneSeek Genomic Profiler (GGP) 70k HD Porcine Array, using the slaughter day and the sex of the animals in the statistical analyses. The phenotypic traits were associated with the genotypes through GenAbel software. The results showed that 18 SNPs located on nine porcine chromosomes were found to be associated with nine phenotypic traits, mainly related to ham weight loss during salting. New associations were found between markers in the genes Neural Precursor Cell Expressed Developmentally Down-Regulated 9 (NEDD9, SSC7), T-Cell Lymphoma Invasion and Metastasis 2 (TIAM2, SSC1), and the ham quality traits. After validation, these SNPs may be useful to improve the quality of thighs for the production of PDO dry-cured hams.

1. Introduction

The production of Protected Designation of Origin (PDO) Parma and San Daniele hams plays an economic role of primary importance in the Italian pig production chain [1] and represents a point of excellence for the Italian pork chain. The quality of the raw meat and the carcass composition greatly influence both the suitability of the thighs to obtain high-quality PDO dry-cured hams such as Parma and San Daniele [2,3]. The fresh thighs used for the production of PDO Parma hams are obtained from heavy pigs slaughtered at a live weight of at least 140 to160 kg, with an age of at least nine months, and belonging to specific selected breeds defined by Parma ham Consortium [4]. These requirements are essential for obtaining hind legs with features suitable for long curing periods; additional required parameters of the thighs are the presence of subcutaneous fat of 15 to 30 mm of thickness to minimize processing losses [2]. These features affect both the aptitude of meat to adsorb salt during the ham salting periods and the weight loss of the ham after salting during the long maturation, which may last from 12 to 24 or more months [5]. However, over time, extensive analytical surveys were carried out to determine the contribution of this food to sodium intake, and medical and public health organizations recommended to reduce sodium dietary intake to prevent hypertension and other diseases [5,6]. On the other hand, an excessive reduction in the amount of adsorbed salt can worsen the technological performances, causing an increase in proteolysis, higher softness, and the production of off-flavors due to the release of free amino acids and peptides during ripening [5,6]. The research on environmental and genetic factors involved in salt adsorption and in the sodium amount of dry-cured ham has attracted great interest based on these contrasting nutritional and technological requirements behind the issue of sodium reduction. Several factors contribute to influencing the final salt amount in dry-cured hams; among them, a crucial role is represented by the know-how of the different ham processors [5]. The ham weight and size, the inter- and intramuscular fat content, the thickness of subcutaneous fat, and the lean meat content of the hind leg represent the main factors that can influence also the aptitude of the ham to adsorb salt [7,8,9,10]. The study and identification of genes and genetic markers associated with these traits could provide new and important tools for the ham industry for improving ham yield during processing. Also the physico-chemical modifications of thigh muscles post-mortem, such as pH decline, water activity, proteolytic and lipolytic reactions, produce changes in color, taste, flavor, texture [11], and salt absorption in the hams [5]. These traits are affected, at least in part, by the animal genetic background as reported by studies indicating genetic markers associated with Semimembranosus muscle (SM) pH and drip loss [12,13], and with enzymes activity in San Daniele hams [14]. Most of these association studies have been performed on the purebred pig lines used in the two-way or three-way crossbreeding schemes [12,15], but few association studies exist between genetic markers and ham ripening related traits measured directly in heavy pigs used for dry-cured ham production. Moreover, the genetic improvement of quality traits of dry-cured hams is difficult because these characteristics cannot be easily measured in live animals and are cost-effective phenotypes.
In recent years, different non-invasive technologies based on X-rays or near-infrared (NIR) spectroscopy were investigated as quality control techniques in the ham industry [16]. Magnetic induction technology is among the most promising to estimate lean content in green hams [17] and the salt absorbed by processed hams [18]. These new non-invasive technologies allow for the collecting of a huge number of phenotypic data. Thus, the identification of genetic markers and new tools to assess and help to predict these traits may be of primary importance for pig production chain and ham processors [10,14,19].
This work aims to investigate the variability of green ham traits in a sample of 229 commercial hybrid heavy pigs during salting and the genomic regions associated with the measured phenotypes. These regions will be singled out through a Genome-Wide Association Study (GWAS), which will allow identifying DNA markers and candidate genes associated with phenotypic characteristics of the fresh hams measured during the first stages of the ham processing. The identified SNPs might be of interest to improve ham yield after salting and the quality of the finished product. The identification of genetic markers associated with the quality of the thighs for PDO dry-cured hams can offer new helpful tools for the highly innovative pig production chain.

2. Materials and Methods

2.1. Animal Data Availability

The samples and hams used in the present study were obtained from the slaughtering of commercial heavy pigs intended for human consumption, and thus the present research did not need approval from a research ethics committee. The heavy pig hybrids were slaughtered in commercial abattoirs in four slaughter days between March and July 2018. The animals were slaughtered in compliance with the European rules [20,21] on the protection of animals during transport and related operations and slaughtering. All slaughter procedures were monitored by the veterinary team appointed by the Italian Ministry of Health. The carcasses came from crossbred pigs reared in three different farms and were obtained from the crossbreeding of the three main breeds reared in Italy for heavy pig production, namely Large White, Landrace, and Duroc. These commercial hybrids were kept in collective pens until they reached the slaughter weight. The pigs were fed diets complying with the dietary inclusion limits established by the Parma Ham Consortium [4] and were slaughtered at an average live weight of 160 kg. Each farm sent their heavy pigs to a specific abattoir, and the thighs were then processed by a specific processing plant. Thus, the farm of farrowing, abattoir, and ham processing plants are collinear variables. Similarly, also the slaughter days match with the farms, abattoirs, and plants, with the thighs coming from the same slaughter day for the hams of plant 1, from two slaughter days for hams processed in plant 2, and from one day of slaughtering for ham processed in plant 3.

2.2. Carcass and Ham Traits

A total of 229 fresh hams were collected from carcasses classified as “U”, “R”, and “O” in compliance with Commission Implementing Decision, 2014/38/EU [22]. Pig carcasses were graded according to the European EUROP carcass grading system [22] which considers the lean meat content estimated for the whole carcass. In particular, carcasses classified as “U” are characterized by a lean meat content between 50% and 55%, carcasses “R” are characterized by a lean meat content between 45% and 50%, and carcasses “O” are characterized by a lean meat content between 40% and 45%. Fresh hams were elaborated at three different dry-cured ham processing plants operating in accordance with tutelary regulations of Parma ham manufacturing [23,24]. A sample of SM was collected for the genomic analyses, and frozen at −20 °C until processed for DNA extraction.
The pH of fresh hams was measured in SM with a Hamilton glass electrode probe attached to a portable pH meter (WTW pH3110, Weilheim, Germany). Homogeneous sets of hams in terms of muscle pH at 24 h post-mortem (pHu) ranged between 5.50 and 5.90 were used for each plant.
Then, all fresh hams underwent traditional salting for Italian dry-cured ham, based on a two-step addition of salt [25], following the standard procedure of each plant.

2.3. Non-Invasive Magnetic Induction (MI) System Analysis

The lean amount of green hams and the salt content of the salted ones were determined by the Ham InspectorTM apparatus (Lenz, Barcelona, Spain). The system based on electromagnetic induction (MI) technology generates a signal with an amplitude depending on the ham lean amount and salt content, estimated using a proper calibration [17,18].
Two hundred and twenty-nine fresh hams were scanned with the MI system set in the “RAW” mode, installed at each plant. The lean amount expressed as a percentage of ham weight was estimated by using previously developed predictive models reported in Simoncini et al. [17]. This model was improved, including a greater number of dissected and analyzed hams, and estimating the new prediction accuracy (RMSE = 1.34%). After the first salting and at the end of the salting steps, the unabsorbed salt was brushed away, and the same salted hams were scanned with the MI system set in the “salted” mode. The salt content of the lean part, expressed as a percentage on a wet basis, was estimated by using a previously developed predictive model, in accordance with Schivazappa et al. [18]. Prediction accuracy of the model including a greater number of dissected and analyzed salted hams was estimated (RMSE = 0.14%).
All scanned hams were tempered at 3 ± 0.5 °C and temperature tested with the thermometer Ebro TFX 410 Pt1000 (Xylem Inc, Rye Brook, NY, USA) inserted into SM muscle (5 cm depth) to avoid possible drifts of MI signal, associated with variations in sample temperature.
The weight of the fresh hams as well as their weights during the salting process (after the first salting and at the end of the salting period) was recorded by the MI apparatus in order to calculate the corresponding weight losses, expressed as percentage loss of fresh ham weight.

2.4. Statistical Analyses of Ham Traits

All statistics of phenotypic traits were obtained by SPSS version 22.0 software platform (SPSS Inc., Chicago, IL, USA); normal distribution of data was investigated before statistical analyses. The boxplot procedure was applied to present the distribution of measured lean percentage of hams, differently labeled according to the EUROP grid (lean carcass grading). Data of ham traits were analyzed using Generalized Linear Model (GLM) procedure; the models included the processing plant and the sex of the pigs respectively, as fixed factors. These factors were included in the model since they are known to affect ham quality [3]. The Least Significant Difference (LSD) posthoc test was applied to compare the Estimated Marginal Means (EMMs). Finally, the Pearson’s correlation analysis was performed to investigate the relationships between all quality and technological traits of hams in each processing plant.

2.5. Genotyping and Association Study

The collected SM samples were used for the genomic analyses. DNA extractions were carried out using a standard protocol by an outsource laboratory (Agrotis S.r.L.-LGS, Cremona, Italy, http://www.lgscr.it/ENG/index.html) where the genotype analyses were also performed. For the genotyping, the GeneSeek® Genomic Profiler-GGP-70k HD Porcine chip (Illumina, San Diego, CA, USA; https://emea.illumina.com/products/by-type/microarray-kits/ggp-porcine.html) containing 68516 SNPs was used using the procedures indicated by the company. The SNPs were mapped using Sus scrofa Genome Assembly Build 11.1 (NCBI: https://www.ncbi.nlm.nih.gov/assembly/GCF_000003025.6; ENSEMBL: https://www.ensembl.org/Sus_scrofa/Info/Index). The sex of the animals (females, castrated males) was obtained on the basis of the genotypes for the sex chromosomes. The genotypic data were filtered using gPLINK (version 2.050, based on PLINK version 1.07) [26] and GenAbel (version 1.8-0, run on R version 3.4.4) [27] discarding the markers or samples not passing these thresholds: all the markers with call rate < 90%, with minor allele frequency < 5%, which are not in Hardy-Weinberg equilibrium [28] (p-value < 0.001); all the individuals with more than 10% missing genotypes and with Identity by State (IBS) > 90%. After filtering, 54,569 SNPs and 169 individuals were retained. The large number of excluded animals (n = 60) was because their IBS was above the threshold set at 90%. The association analyses were carried out using polygenic_hglm and qtscore functions implemented on the GenAbel package according to the procedure described in Nicolazzi et al. [29]. The statistical model included the day of slaughter and the sex as predictive variables, together with the effect of the SNP and the genomic kinship matrix. The genomic kinship matrix was obtained with GenAbel and was used to estimate the relatedness between the animals. Pedigree data were not available for these animals since commercial hybrids are obtained from heterospermic artificial insemination. Also, the information concerning litters and dams was not available. These data can be available in fully controlled experimental studies, but they are difficult to obtain in commercial farm conditions, where cross-fostering is a common practice. The p-values of the associations between ham traits and SNPs were corrected for the “deflation” factor, as reported by GenAbel. Markers were considered significant with corrected p-values below the chromosome-wide significant threshold (Table S1), and the trend towards the significance threshold was set at 5.00 × 10−5 [30]. The markers with corrected p-values below these thresholds were further considered for the estimation of the genotypic effects. Their genotypes were extracted from the PED file generated through the gPLINK tool using in-house implemented scripts in the R environment [31]. Linear mixed models were performed to obtain the Estimated Least Squares Means (LSM) for the genotypes of the significant markers. The associated traits were used as dependent variables; the genotype of each marker taken individually was modeled as a predictive variable together with the fixed effects of sex and slaughter day. The linear regression model did not include the random effect of the genomic kinship matrix. For each marker found with the GWAS, the additive and dominance genetic effects were also estimated. The additive effect was estimated as half of the difference between the two homozygous groups: a = 1/2(BB − AA), with A and B that indicate the first and the second allele of the analyzed markers, respectively. The dominance effect was estimated as the difference between the heterozygous group and the average of the two homozygous groups in each locus: d = AB − 1/2(AA + BB). These analyses were performed using functions in the packages nlme [32], lsmeans [33], lme4 [34], and car [35] in the R environment [31]. The genes located in the region flanking the identified markers (±500 kilobases from the associated marker) were further considered for the identification of candidate genes for the traits. The list of the flanking genes was obtained using the BioMart tool [36] and was submitted to David Bioinformatics Resources version 6.8 on-line tool (https://david.ncifcrf.gov/). Candidate genes were identified on the basis of their location (the nearest gene to the significant marker) and biological role. In order to find possible splice sites or motifs with biological relevance, the sequences flanking the intronic variants found associated with ham traits were submitted to the Tomtom tool in MEME Suite version 5.3.0 [37] (http://meme-suite.org/). The found motifs are ranked based on their Bonferroni significance (q-value) for the found match between the query sequence (or its complementary) and the motif in JASPAR CORE 2014 database.

3. Results

3.1. Ham Quality and Technological Traits

The variability noticed for ham lean % is reported in Figure 1. The boxplot displays the distribution of the lean % of hams estimated by Ham InspectorTM, in relation to the lean % of the corresponding carcasses according to EUROP classification [22]. The EUROP classification of the original carcass is printed on the skin of each fresh ham supplied to dry-cured ham producers, to provide information on the lean % expected in fresh ham.
The class R, corresponding to an estimated lean % of the carcass in the range 45 to 50% is the most abundant (n = 110) compared to classes O (lean % = 40 to 45%, n = 36) and U (lean %= 50 to 55%, n = 62). In the current study, the lean % of the hams scanned by Ham InspectorTM shows higher values (R = 62.3 ± 2.7, U = 63.0 ± 2.0, O = 61.6 ± 2.2) than those estimated from the EUROP classification of the corresponding whole carcass.
All fresh hams were processed in three different manufacturing plants in accordance with tutelary regulations of Parma ham [23,24].
The results of GLM analysis including “processing plant” and “sex” as fixed effects and their interaction are summarized in Table 1. For each trait, the EMMs (Estimated Marginal Means), the standard error, and the corresponding significance level are reported.
Differences in the measured ham traits between processing plants (p < 0.001) were detected. In particular, processing plant 1 handles the heaviest hams between the tested plants (15.1 kg), and hams leaner (62.7%) than plant 3 (60.6%). This condition affects salting weight losses, yielding higher values at the end of salting for plant 1 compared to other plants. The salt % predicted in ham lean content after both salting steps differed significantly between processing plants (p < 0.001), reporting greater values for plant 3 and 2, respectively (Table 1).
The sex of the pigs generated differences for the lean % of green hams (p < 0.001) and their weight loss after the first (p < 0.05) and at the end of the salting (p < 0.01; Table 1), and the predicted salt contents (p < 0.01). The higher lean % of green hams from females caused higher weight losses during salting, and higher values of salt content at the end of salting (2.60%) if compared to barrows (2.53%). As for the salt content at the end of salting, the interaction between processing plant and sex was significant (p < 0.05): the three plants differed regardless of the sex of the animal, but in the case of females the salt values were higher than in barrows.
The results of the Pearson’s correlation analysis performed between ham traits measured in each processing plant are reported in Tables S2–S4. The strongest and most significant correlations were found in plant 2 (Table S3): pH was negatively correlated with salt content measured at the first salting step (r = −0.241, p < 0.01) and the weight loss at the end of salting (r = −0.388, p < 0.001). The weight of green hams was negatively correlated with salt content after the first salting in all tested processing plants (Tables S2–S4), whereas, only for plant 2, the negative correlation coefficient with salt content at the end of the salting process was significant (r = −0.486, p < 0.001; Table S3). For all processing plants (Tables S2–S4), the lean % of green hams was positively related to ham weight losses and to the salt contents predicted after the first and at the end of the salting steps. In each processing plant, the highest positive correlation coefficients have been found between green ham weight and ham weight after the first salting (plant 1: r = 0.999, p < 0.001; plant 2 r = 0.999, p < 0.001; plant 3 r = 0.998, p < 0.001), green ham weight and ham weight at the end of salting (plant 1: r = 0.997, p < 0.001; plant 2 r = 0.999, p < 0.001; plant 3 r = 0.996, p < 0.001), and between ham weight after first and at the end of salting (plant 1: r = 0.997, p < 0.001; plant 2 r = 0.999, p < 0.001; plant 3 r = 0.999, p < 0.001).

3.2. Association Study Results

The performed GWAS allowed the identification of eight markers displaying a Bonferroni corrected p-value significant at the chromosome-wide level, reported in bold in Table 2. The most significant association was found for pHu, which is linked to WU_10.2_18_17949287 marker. This SNP is an intergenic variant located at 17 Mb on Sus scrofa chromosome 18 (SSC18). Green ham lean % showed an association with a region on SSC4 at 2 Mb, where the marker ASGA0016987 is located. This SNP is an exon variant of a non-coding transcript (ENSSSCT00000066959.1), and most of the genes comprised in the flanking region are non-coding genes. Green ham weight, ham weight after first salting, and ham weight at the end of the salting showed an association with the same two markers, namely CASI0010463 and WU_10.2_14_144250775. The first is located on SSC15 at 12 Mb, in an intron of the gene Neurexophilin 2 (NXPH2). WU_10.2_14_144250775 is an intergenic variant in a region characterized by high gene density, on SSC14 at 132 Mb. Two SNPs were significant for ham weight loss after first salting: ASGA0031014, an intron variant of the gene Neural Precursor Cell Expressed, Developmentally Down-Regulated 9 (NEDD9), and WU_10.2_10_74620421, a variant in a non-coding transcript on SSC10. Salt content at the end of salting showed two markers displaying a chromosome-wide significant association. The first (INRA0000796) is located on SSC1 at 13 Mb, and is an intron variant of the gene Regulator of G Protein Signaling 17 (RGS17); the second (ASGA0102337) is an intron variant of the gene GRAM Domain Containing 1B (GRAMD1B), located on SSC9 at 50 Mb. Eleven markers showed less significant associations with the ham traits, with adjusted p-values comprised between the chromosome-wide significance level and the threshold of 5.00 × 10−5 (Table 2). The complete lists of the top 20 markers most associated with the measured traits are reported in File S1. Table S5 reports the complete list of all the genes located in the regions flanking the markers reported in Table 2. The results of the functional association obtained for the genes in Table S5 are reported in Table S6. No significant terms were identified, and none of them seemed to indicate a direct involvement of the candidate genes in muscle or fat development. Thus, further discussion of the obtained associations was based on the review of the scientific literature.
For each of the found markers, the LSM of the genotypes were estimated and reported in Table 3. For pHu, both markers showed genotype distributions quite unbalanced, with GG genotype poorly represented (WU_10.2_18_17949287) or not represented at all (WU_10.2_4_91195648) in our sample. Other identified markers showed a genotypic class poorly represented (n < 10 animals, i.e., ASGA0016987, ALGA0002237, CASI0010463, WU_10.2_7_118557013) or completely lacking (ALGA0044906, INRA0000796). The weights of fresh thighs and hams after first and at the end of salting were associated with the same markers, namely CASI0010463 and WU_10.2_14_144250775. Furthermore, the marker ASGA0026341 showed an additive genetic effect on both ham weights after first salting and at the end of salting. In fact, the AA animals for this marker showed lower ham weights than GG pigs (Table 3). Additive genetic effects were also found for H3GA0000815, WU_10.2_7_118557013, WU_10.2_14_36295226, and ASGA0031014 (Table 3). In particular, the AA animals for the H3GA0000815 SNP were associated with a lower lean % in green hams; the A allele of WU_10.2_7_118557013 was related to higher weights at the end of salting; the G allele of ASGA0031014 showed lower weight losses, and the GG animals for the WU_10.2_14_36295226 SNP presented a higher % of salt adsorbed at first salting. Concerning the ASGA0031014 marker, the three genotypic classes were well represented, following the 1:2:1 ratio.
The results obtained with MEME Suite indicated that none of the identified intronic variants were located in splice sites, but the marker ASGA0031014 fell into an intronic region of the gene Neural Precursor Cell Expressed, Developmentally Down-Regulated 9 (NEDD9) harboring several binding sites recognized by Forkhead box O (FOXO) transcription factors (Table S7). The in-silico analysis indicated that the mutation ASGA0031014 changed the sequence recognized by FOXO3 (q-value = 0.050; Figure S1a), FOXO1 (q-value = 0.074; Figure S1b), and FOXO4 (q-value = 0.074; Figure S1c). The sequences recognized by these transcription factors are the reverse complement of the region flanking ASGA0031014, since the NEDD9 gene is located on the reverse strand (strand −1 in Table S5). For this reason, ASGA0031014 has T/C as alternate alleles in Figure S1.

4. Discussion

The development of new technologies enabling a fast and non-invasive measure of green ham potential to be processed into typical dry-cured ham is a goal for pig production chain and PDO ham producers. In the last decade, technologies such as those based on X-rays have been studied [16], with the aim of monitoring and gathering information in a non-invasive way on green hams intended for dry curing, focusing on their qualitative traits and the amount of salt absorbed during processing. The Ham InspectorTM system permits an on-line and non-invasive inspection of the hams, a classification of thighs based on the lean content (calculated as a percentage of green ham weight), and an estimation of salt amount adsorbed in the lean part during the salting stages. The use of Ham InspectorTM in three processing plants qualified for the production of PDO Parma ham provided the phenotypic data and information to perform the GWAS. At present, green hams allowed to be processed into PDO Parma hams are those labeled as U, R, and O according to the EUROP grid, with a few specific exceptions [38]. In the current study, a high variability in ham lean % within U, R, and O classes was detected. In particular, the thigh lean % estimated from Ham InspectorTM was higher than the lean % of the corresponding carcasses estimated from EUROP grading. Thus, the thigh lean % given by Ham InspectorTM and the carcass lean % obtained from EUROP grading provide different information. This is anyway expected, because ham is a lean cut, with a higher lean % than that of the corresponding EUROP carcass class. The measures estimated by Ham InspectorTM about the lean-to-fat ratio of each thigh proved to be crucial to homogenously group hams during processing [6]. Processing plants showed differences for all the measured phenotypes, except for ham weight loss after the first salting. This result must be carefully evaluated as the on-field conditions of the present study did not permit to distinguish the specific effect of the processing plant, which included also the effects of the pig genetic type, farm management, and abattoir. Thus, the effect of the processing plant may be overestimated in the present study. The effect of the lean amount of green hams on the processing weight losses is underlined by the positive and significant correlation coefficients displayed in the three plants. This is in line with previous findings: weight losses are indeed known to be greater in hams obtained from animals with greater lean %, whereas the presence of inter, intramuscular and covering fat, containing less water than muscular tissue, is associated with a reduction of weight losses [2]. As can be expected, the lean % of green hams was strongly related to salt content after the first salting period and at the end of salting: this behavior is related to the higher salt and water diffusion coefficients of Fick’s law in lean hams than in the fat ones [39]. Currently, the variability detected in the phenotypic traits of processed ham as a consequence of the differences in raw matter and in the production plants remains in full compliance with the tutelary guidelines of Parma hams [23].
With the aim of finding genomic regions and candidate genes associated with the green ham weight, lean %, salting losses, and contents of adsorbed salt, we used sex and slaughter day (and thus also the ham processing plant) as fixed effects in the GWAS model. GWAS indicated eight markers significant at the chromosome-wide level, and eleven with a Pc1df below the threshold of 5.00 × 10−5. The most significant association was found between pHu and the WU_10.2_18_17949287 marker, which is located in an intergenic region of SSC18, harboring the genes Plexin A4 (PLXNA4), and Muskelin 1 (MKLN1). This association signal with the PLXNA4 region could support the findings reported by Bordbar et al. [40] in Simmental beef cattle. These authors found that the bovine PLXNA4 gene harbored 18 SNPs significantly associated with muscle development in the Simmental breed. Despite this strong association, no clear biologic evidence supporting the PLXNA4 role in bovine muscle development was found [40]. The pHu variability was also found associated with the WU_10.2_4_91195648 marker, on SSC4. Anyway, both WU_10.2_18_17949287 and WU_10.2_4_91195648 markers showed unbalanced genotype distributions, which may have biased the association results, causing an overestimation of the marker effects.
Thigh lean % was associated with the marker H3GA0000815 located in an intronic region of the gene T-Cell Lymphoma Invasion And Metastasis 2 (TIAM2). Among the identified associations, the H3GA0000815 marker showed the clearest additive effect, with the GG animals showing a higher lean % than those with the AA genotype. The lower frequency observed for the GG genotype in our samples suggests that the favorable allele G may be used in selection schemes to increase lean mass deposition in the pig breeds used for the production of PDO hams. Interestingly, the TIAM2 gene was found in a previous study to be located in a region displaying a selective signature in a Duroc pig population [41]. This gene encodes a guanine nucleotide exchange factor with a role in intracellular signal transduction and in the regulation of cell migration and cell focal adhesions [42]. Cell migration and focal adhesions are essential steps in organogenesis, and also muscle development is known to involve a series of morphogenetic events including cell fusion, migration, and epidermal attachment [43]. Despite the fact that this result seems to agree with a previous study [41], to date, it is not possible to draw a clear hypothesis linking this gene with thigh lean % as the knowledge of the TIAM2 gene and its roles in pig muscle development is still mostly unknown.
The variability of green ham weight, and ham weight after first salting and at the end of the salting period showed to be associated with the same markers, indicating a pleiotropic effect of these mutations. The LSM estimated for these markers were concordant among the measured traits, with the same genotypic class displaying higher estimated LSM for the weights of the fresh thigh and the ham after first salting and at the end of salting. This pattern was expected given the degree of shared variability between these three traits, as demonstrated by their high and positive phenotypic correlations, indicating that the higher the weight of the fresh thigh, the heavier are the same hams after first salting and at the end of salting period. Ham weights at the first salting and at the end of salting were significantly associated with the two markers ASGA0026341 and WU_10.2_7_118557013. An additive effect was found for both these markers, which are located in intergenic regions on SSC5 and SSC7, respectively. This is not surprising, since the literature indicates that the majority of GWAS peaks are located in non-coding or intergenic regions [44], and may therefore be a signal indicating that the found association may be due to causal mutations in closely located genes or to mutations in a still unknown gene. Among the genes located near WU_10.2_7_118557013 is Calmodulin 1 (CALM1), which encodes for one of three calcium-binding calmodulin proteins. CALM1 mediates the control of a large number of enzymes, ion channels, aquaporins, and other proteins through calcium-binding, and controls also the transport of glucose and other sugars in hepatic cells [45]. Its gene expression was also found to be downregulated in the Longissimus muscle of pigs with a high deposition of intramuscular fat [46], suggesting that this gene may have a role also in the metabolic pathways influencing pig meat quality.
The ham weight loss after first salting was found to be associated with ASGA0031014, on SSC7. The alleles for this marker show a significant additive effect and balanced genotype distributions in the considered sample. The most favorable genotype for this locus was GG, which displayed the lowest weight losses during processing. This mutation is an intron variant of the gene NEDD9, which codes for a focal adhesion protein that acts as a scaffold to regulate signaling complexes important in cell attachment, migration, and invasion [47], as well as apoptosis and the cell cycle [48]. Some studies of different authors reported that NEDD9 is essential for the Transforming Growth Factor β (TGF-β) signaling pathway [49,50]. Moreover, members of the TGF-β superfamily can profoundly regulate mesenchymal stem cell differentiation, as well as adipogenesis and myogenesis [51,52]. Due to the involvement of NEDD9 in TGF-β signaling, it is, therefore, possible that changes in its gene sequence may affect also TGF-β signaling and the differentiation of mesenchymal stem cells during muscle development. Furthermore, NEDD9 is known to take part in the inhibition of primary cilia formation during embryogenesis [53]. Primary cilia are non-motile organelles that have been shown to play an important role as antennae to extracellular stimuli during embryogenesis [53], adipogenesis, and myofibrogenesis [54]. In our previous study, we found several genes related to primary cilia associated with intramuscular fat deposition [55], and NEDD9 may take part in this complex signaling pathway related to muscle and adipose tissue development in embryogenesis. Interestingly, the marker ASGA0031014 fell into an intronic region harboring several sequences recognized by FOXO transcriptional regulators. In particular, this intronic variant is located in a site matching the motifs recognized by FOXO1, FOXO3, and FOXO4 transcription regulators, suggesting that changes in this sequence may also affect the ability of those transcription regulators to recognize this binding site. This hypothesis seems to be supported in the literature, which indicated that NEDD9 gene expression was found to be regulated by another forkhead box member (FOXC1) [56] and that FOXO1 has a major role in the regulation of skeletal muscle differentiation and fiber type specification in mammals [57]. Taken together, these results suggest that ASGA0031014 may be a useful marker for the improvement of ham yield during processing. Our hypothesis needs anyway a validation with further specific studies.
In addition to these markers, other mutations were associated with the ham traits, even though they did not present significant additive effects. Among the markers associated with the ham weights at the first and at the end of the salting steps is CASI0010463, an intron variant of the Neurexophilin 2 (NXPH2) gene. A Run of Homozygosity (ROH) region was found in the NXPH2 gene in Valdostana Black Pied and Valdostana Chestnut bovine breeds [58]. This evidence seems to suggest a possible role of NXPH2 in livestock traits related to muscle development and production efficiency. However, this gene has never been studied before in pigs and the limited knowledge of the NXPH2 gene does not permit to reach a complete interpretation of the associations found for this gene in the present study. Another marker associated with green ham weight, ham weight after first salting, and ham weight at the end of salting is WU_10.2_14_144250775. This marker results to be unmapped in the latest release of the Ensembl database [59], while it was located at 132,664,262 bp on SSC14 in Illumina 11.1 mapping. Following the latter mapping, WU_10.2_14_144250775 is an intergenic variant. The closest genes to this marker are Acyl-CoA Dehydrogenase Short/Branched Chain (ACADSB), H6 Family Homeobox 2 (HMX2), H6 Family Homeobox 3 (HMX3), and BUB3 Mitotic Checkpoint Protein (BUB3). The sequence comprising ACADSB, HMX2, HMX3, and BUB3 is highly conserved between mammals [60], suggesting that this region could be of primary biological importance. Anyway, no clear association exists between muscle development and deposition, and these genes. Further investigation should be dedicated to testing the role of these genes located in this SSC14 region on production traits and porcine muscle development. Also of interest is the marker ASGA0102337, an intron variant of the GRAM Domain Containing 1B (GRAMD1B) gene. This gene codes for a cholesterol transporter that mediates non-vesicular transport of cholesterol from the plasma membrane to the endoplasmic reticulum [61,62]. GRAMD1B gene expression in muscle was found associated with obesity in the mouse model [63], and mutations in its sequence were associated with feed efficiency in a beef cattle population with individuals from various breeds [64]. As fat deposition is important in determining salting losses and salt uptake of the green hams, mutations in the GRAMD1B sequence may also influence salting and maturation traits. However, the possible effect of this gene on ham quality traits needs confirmation with studies in larger samples.

5. Conclusions

To our knowledge, this is the first study describing significant associations between porcine candidate genes and ham traits, measured in a non-invasive way using a Magnetic Induction System. This on-line and non-invasive technology allows the estimation of green ham lean % and the salt content in the salted ones, on a great number of hams during processing. The GWAS identified several markers associated with the measured traits, and SNPs in candidate genes related to the neuromuscular junction, and muscle development during embryonic stages were found associated with lean % and ham weights. Among the markers associated with ham traits, the marker located in NEDD9 could be of particular interest for further studies and for implementing selection schemes aimed at improving ham yields after the salting steps. Further studies including a greater number of hams and different processing conditions are needed to strengthen the investigated associations.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/11/1/68/s1, Table S1: The Chromosome-wide thresholds for significance calculated for each chromosome; Table S2: Pearson’s correlation coefficients (above the diagonal) and corresponding p-values (below the diagonal) between processing and analytical traits of hams (green and after 1st and at the end of salting steps) in Plant 1; Table S3: Pearson’s correlation coefficients (above the diagonal) and corresponding p-values (below the diagonal) between processing and analytical traits of hams (green and after 1st and at the end of salting steps) in Plant 2; Table S4: Pearson’s correlation coefficients (above the diagonal) and corresponding p-values (below the diagonal) between processing and analytical traits of hams (green and after 1st and at the end of salting steps) in Plant 3; Table S5: List of the genes in the regions flanking the markers associated with ham traits; Table S6: Results of the functional enrichment analysis performed on the list of the genes located in the regions flanking the markers associated with ham traits; Table S7: Results of the Tomtom tool with the binding sites recognized by different transcription regulators on the sequence flanking the marker ASGA0031014; Figure S1: The sequence flanking the marker ASGA0031014 (position 17) is recognized by the transcription factors (a) Forkhead box O 3 (FOXO3; upper sequence), with the alternate allele C (or G if we consider the complementary DNA sequence) possibly changing the binding site recognized by FOXO3 (q-value = 0.050); (b) Forkhead box O 1 (FOXO1; upper sequence),with the alternate allele C (or G if we consider the complementary DNA sequence) possibly changing the binding site recognized by FOXO1 (q-value = 0.074); (c) Forkhead box O 4 (FOXO4; upper sequence), with the alternate allele C (or G if we consider the complementary DNA sequence) possibly changing the binding site recognized by FOXO4 (q-value = 0.074); File S1: Tables with the top 20 markers associated with each ham trait and their complete output information obtained with GenAbel package.

Author Contributions

Conceptualization, M.Z. and R.D.; methodology, M.Z., N.S., and R.V.; formal analysis, P.Z., M.Z., N.S., and C.S.; investigation, M.Z., P.Z., N.S., B.S., C.S., R.V., and R.D.; resources, C.S., N.S., and M.Z.; data curation, P.Z., N.S., C.S., and R.D.; writing—original draft preparation, M.Z., N.S., R.D., and R.V.; writing—review and editing, M.Z., P.Z., N.S., R.D., and R.V.; supervision, R.D. and B.S.; project administration, R.V., R.D., and M.Z.; funding acquisition, R.D. and R.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the AGER CARIPLO Italian foundation, grant number 2017–2022, project “ProSuIT-Tecnologia (T) a favore (Pro) della suinicoltura italiana (Su-I)”, and by AGER-Agroalimentare e ricerca: Advanced research in genomics and processing technologies for the Italian heavy pig production-Hepiget project, grant number 2011-0279.

Institutional Review Board Statement

Ethical review and approval were waived for this study, as the samples and hams used in the present study were obtained from commercial heavy pigs intended for human consumption. This study did not involve live animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors. The data are not publicly available due to confidentiality agreements. Supporting data can be made available to bona fide researchers subject to a non-disclosure agreement.

Acknowledgments

The authors thank Raffaele Mazza from Agrotis-Laboratorio Genetica e Servizi-Associazione Italiana Allevatori (Cremona, Italy) for performing the genotyping and for his support and help. The authors thank Illumina Technical support for providing the updated manifest file with the SNPs annotated with version 11.1 of the porcine genome. The authors thank Tania Toscani from the Consorzio del Prosciutto di Parma for her support and help, and the Consorzio del Prosciutto di Parma for co-financing the trials with Ham InspectorTM in ham processing plants.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Carni, I.M. Carne Suina e Salumi. Available online: http://www.ismeamercati.it/carni/carne-suina-salumi (accessed on 22 October 2020).
  2. Bosi, P.; Russo, V. The production of the heavy pig for high quality processed products. Ital. J. Anim. Sci. 2004, 3, 309–321. [Google Scholar] [CrossRef] [Green Version]
  3. Candek-Potokar, M.; Skrlep, M. Factors in pig production that impact the quality of dry-cured ham: A review. Anim. Int. J. Anim. Biosci. 2012, 6, 327–338. [Google Scholar] [CrossRef] [PubMed]
  4. Prosciutto di Parma (Parma Ham) Protected Designation of Origin (Specifications and Dossier Pursuant to Article 4 of Council Regulation EEC no. 2081/92 Dated 14 July 1992). Offic. J. L 1992, 208, 9–14.
  5. Schivazappa, C.; Virgili, R. Impact of salt levels on the sensory profile and consumer acceptance of Italian dry-cured ham. J. Sci. Food Agric. 2020, 100, 3370–3377. [Google Scholar] [CrossRef] [PubMed]
  6. Pinna, A.; Saccani, G.; Schivazappa, C.; Simoncini, N.; Virgili, R. Revision of the cold processing phases to obtain a targeted salt reduction in typical Italian dry-cured ham. Meat Sci. 2020, 161, 107994. [Google Scholar] [CrossRef]
  7. Fernández, A.; de Pedro, E.; Núñez, N.; Silió, L.; García-Casco, J.; Rodríguez, C. Genetic parameters for meat and fat quality and carcass composition traits in Iberian pigs. Meat Sci. 2003, 64, 405–410. [Google Scholar] [CrossRef]
  8. Davoli, R.; Catillo, G.; Serra, A.; Zappaterra, M.; Zambonelli, P.; Zilio, D.M.; Steri, R.; Mele, M.; Buttazzoni, L.; Russo, V. Genetic parameters of backfat fatty acids and carcass traits in Large White pigs. Anim. Int. J. Anim. Biosci. 2019, 13, 924–932. [Google Scholar] [CrossRef]
  9. Davoli, R.; Luise, D.; Mingazzini, V.; Zambonelli, P.; Braglia, S.; Serra, A.; Russo, V. Genome-wide study on intramuscular fat in Italian Large White pig breed using the PorcineSNP60 BeadChip. J. Anim. Breed. Genet. 2016, 133, 277–282. [Google Scholar] [CrossRef]
  10. Davoli, R.; Schivazappa, C.; Zambonelli, P.; Braglia, S.; Rossi, A.; Virgili, R. Association study between single nucleotide polymorphisms in porcine genes and pork quality traits for fresh consumption and processing into Italian dry-cured ham. Meat Sci. 2017, 126, 73–81. [Google Scholar] [CrossRef]
  11. Toldrá, F.; Flores, M.; Sanz, Y. Dry-cured ham flavour: Enzymatic generation and process influence. Food Chem. 1997, 59, 523–530. [Google Scholar] [CrossRef]
  12. Davoli, R.; Zappaterra, M.; Zambonelli, P. Genome-wide association study identifies markers associated with meat ultimate pH in Duroc pigs. Anim. Genet. 2019, 50, 154–156. [Google Scholar] [CrossRef] [PubMed]
  13. Zappaterra, M.; Sami, D.; Davoli, R. Association between the splice mutation g.8283C>A of the PHKG1 gene and meat quality traits in Large White pigs. Meat Sci. 2019, 148, 38–40. [Google Scholar] [CrossRef]
  14. Renaville, B.; Piasentier, E.; Fan, B.; Vitale, M.; Prandi, A.; Rothschild, M.F. Candidate gene markers involved in San Daniele ham quality. Meat Sci. 2010, 85, 441–445. [Google Scholar] [CrossRef] [PubMed]
  15. Fontanesi, L.; Schiavo, G.; Gallo, M.; Baiocco, C.; Galimberti, G.; Bovo, S.; Russo, V.; Buttazzoni, L. Genome-wide association study for ham weight loss at first salting in Italian Large White pigs: Towards the genetic dissection of a key trait for dry-cured ham production. Anim. Genet. 2017, 48, 103–107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Fulladosa, E.; Muñoz, I.; Serra, X.; Arnau, J.; Gou, P. X-ray absorptiometry for non-destructive monitoring of the salt uptake in bone-in raw hams during salting. Food Control 2015, 47, 37–42. [Google Scholar] [CrossRef]
  17. Simoncini, N.; Virgili, R.; Schivazappa, C.; Pinna, A.; Rossi, A.; Álvarez, J. Assessment of fat and lean content in Italian heavy green hams by means of on-line non-invasive techniques. In Proceedings of the 58th International Congress of Meat Science and Technology, Montreal, QC, Canada, 12–17 August 2012; p. 4. [Google Scholar]
  18. Schivazappa, C.; Virgili, R.; Simoncini, N.; Tiso, S.; Álvarez, J.; Rodríguez, J.M. Application of the magnetic induction technique for the non-destructive assessment of salt gain after the salting process of Parma ham. Food Control 2017, 80, 92–98. [Google Scholar] [CrossRef]
  19. Otto, G.; Roehe, R.; Looft, H.; Thoelking, L.; Knap, P.W.; Rothschild, M.F.; Plastow, G.S.; Kalm, E. Associations of DNA markers with meat quality traits in pigs with emphasis on drip loss. Meat Sci. 2007, 75, 185–195. [Google Scholar] [CrossRef]
  20. EUR-Lex—32009R1099—EN—EUR-Lex. Available online: https://eur-lex.europa.eu/eli/reg/2009/1099/oj (accessed on 28 October 2020).
  21. EUR-Lex—32005R0001—EN—EUR-Lex. Available online: https://eur-lex.europa.eu/legal-content/EN/ALL/?uri=celex%3A32005R0001 (accessed on 28 October 2020).
  22. 2014/38/EU: Commission Implementing Decision of 24 January 2014 Authorising Methods for Grading Pig Carcases in Italy (Notified under Document C(2014) 279). Offic. J. Eur. Union 2014, 23, 35–40.
  23. Commission Implementing Regulation (EU) No 1208/2013 of 25 November 2013 Approving Minor Amendments to the Specification for a Name Entered in the Register of Protected Designations of Origin and Protected Geographical Indications (Prosciutto di Parma (PDO)); Publications Office of the European Union: Luxembourg, 2013; Volume 317.
  24. Council Regulation (EEC) No 2081/92 of 14 July 1992 on the Protection of Geographical Indications and Designations of Origin for Agricultural Products and Foodstuffs. Offic. J. L 1992, 208, 1–9.
  25. Parolari, G. Review: Achievements, needs and perspectives in dry-cured ham technology: The example of Parma ham/ Revisión: Avances, necesidades y perspectivas de la tecnología del jamón curado: El ejemplo del jamón de Parma. Food Sci. Technol. Int. 1996, 2, 69–78. [Google Scholar] [CrossRef]
  26. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  27. Aulchenko, Y.S.; Ripke, S.; Isaacs, A.; van Duijn, C.M. GenABEL: An R library for genome-wide association analysis. Bioinform. Oxf. Engl. 2007, 23, 1294–1296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Hardy, G.H. Mendelian proportions in a mixed population. Science 1908, 28, 49–50. [Google Scholar] [CrossRef] [PubMed]
  29. Nicolazzi, E.L.; Biffani, S.; Biscarini, F.; Wengel, P.O.t.; Caprera, A.; Nazzicari, N.; Stella, A. Software solutions for the livestock genomics SNP array revolution. Anim. Genet. 2015, 46, 343–353. [Google Scholar] [CrossRef] [PubMed]
  30. Burton, P.R.; Clayton, D.G.; Cardon, L.R.; Craddock, N.; Deloukas, P.; Duncanson, A.; Kwiatkowski, D.P.; McCarthy, M.I.; Ouwehand, W.H.; Samani, N.J.; et al. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 2007, 447, 661–678. [Google Scholar] [CrossRef] [Green Version]
  31. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  32. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; R Core Team. {nlme}: Linear and Nonlinear Mixed Effects Models; R Package Version 3.1-150. 2020. Available online: https://cran.r-project.org/web/packages/nlme/nlme.pdf (accessed on 28 October 2020).
  33. Lenth, R.V. Least-Squares Means: The R Package lsmeans. J. Stat. Softw. 2016, 69, 1–33. [Google Scholar] [CrossRef] [Green Version]
  34. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  35. Fox, J.; Weisberg, S. An R Companion to Applied Regression; SAGE Publications, Inc: Thousand Oaks, CA, USA, 2011; ISBN 978-1-4129-7514-8. [Google Scholar]
  36. Baker, M. Quantitative data: Learning to share. Nat. Methods 2012, 9, 39–41. [Google Scholar] [CrossRef]
  37. Gupta, S.; Stamatoyannopoulos, J.A.; Bailey, T.L.; Noble, W.S. Quantifying similarity between motifs. Genome Biol. 2007, 8, R24. [Google Scholar] [CrossRef] [Green Version]
  38. Istituto Parma Qualità. Istituto Nord Est Qualità Impiego delle Nuove Equazioni e Criteri di Accettazione delle Carcasse Suine- Circolare Unificata N. 11/2014; 2014. Available online: http://www.parmaqualita.it/IPQ/datiute/6/112014CRITERIACCETTABILITACARCASSE.pdf (accessed on 28 October 2020).
  39. Harkouss, R.; Chevarin, C.; Daudin, J.-D.; Sicard, J.; Mirade, P.-S. Development of a multi-physical finite element-based model that predicts water and salt transfers, proteolysis and water activity during the salting and post-salting stages of the dry-cured ham process. J. Food Eng. 2018, 218, 69–79. [Google Scholar] [CrossRef]
  40. Bordbar, F.; Jensen, J.; Zhu, B.; Wang, Z.; Xu, L.; Chang, T.; Xu, L.; Du, M.; Zhang, L.; Gao, H.; et al. Identification of muscle-specific candidate genes in Simmental beef cattle using imputed next generation sequencing. PLoS ONE 2019, 14, e0223671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Diao, S.; Luo, Y.; Ma, Y.; Deng, X.; He, Y.; Gao, N.; Zhang, H.; Li, J.; Chen, Z.; Zhang, Z. Genome-wide detection of selective signatures in a Duroc pig population. J. Integr. Agric. 2018, 17, 2528–2535. [Google Scholar] [CrossRef] [Green Version]
  42. Rooney, C.; White, G.; Nazgiewicz, A.; Woodcock, S.A.; Anderson, K.I.; Ballestrem, C.; Malliri, A. The Rac activator STEF (Tiam2) regulates cell migration by microtubule-mediated focal adhesion disassembly. EMBO Rep. 2010, 11, 292–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Maartens, A.P.; Brown, N.H. The many faces of cell adhesion during Drosophila muscle development. Dev. Biol. 2015, 401, 62–74. [Google Scholar] [CrossRef] [Green Version]
  44. Watanabe, K.; Taskesen, E.; van Bochoven, A.; Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 2017, 8, 1826. [Google Scholar] [CrossRef] [Green Version]
  45. Chen, Z.; Ding, L.; Yang, W.; Wang, J.; Chen, L.; Chang, Y.; Geng, B.; Cui, Q.; Guan, Y.; Yang, J. Hepatic Activation of the FAM3C-HSF1-CaM Pathway Attenuates Hyperglycemia of Obese Diabetic Mice. Diabetes 2017, 66, 1185–1197. [Google Scholar] [CrossRef] [Green Version]
  46. Hamill, R.M.; McBryan, J.; McGee, C.; Mullen, A.M.; Sweeney, T.; Talbot, A.; Cairns, M.T.; Davey, G.C. Functional analysis of muscle gene expression profiles associated with tenderness and intramuscular fat content in pork. Meat Sci. 2012, 92, 440–450. [Google Scholar] [CrossRef]
  47. Zhong, J.; Baquiran, J.B.; Bonakdar, N.; Lees, J.; Ching, Y.W.; Pugacheva, E.; Fabry, B.; O’Neill, G.M. NEDD9 Stabilizes Focal Adhesions, Increases Binding to the Extra-Cellular Matrix and Differentially Effects 2D versus 3D Cell Migration. PLoS ONE 2012, 7, e35058. [Google Scholar] [CrossRef] [Green Version]
  48. Semelakova, M.; Grauzam, S.; Betadthunga, P.; Tiedeken, J.; Coaxum, S.; Neskey, D.M.; Rosenzweig, S.A. Vimentin and Non-Muscle Myosin IIA are Members of the Neural Precursor Cell Expressed Developmentally Down-Regulated 9 (NEDD9) Interactome in Head and Neck Squamous Cell Carcinoma Cells. Transl. Oncol. 2018, 12, 49–61. [Google Scholar] [CrossRef]
  49. Vogel, T.; Ahrens, S.; Büttner, N.; Krieglstein, K. Transforming Growth Factor β Promotes Neuronal Cell Fate of Mouse Cortical and Hippocampal Progenitors In Vitro and In Vivo: Identification of Nedd9 as an Essential Signaling Component. Cereb. Cortex 2010, 20, 661–671. [Google Scholar] [CrossRef] [Green Version]
  50. Omata, Y.; Nakamura, S.; Koyama, T.; Yasui, T.; Hirose, J.; Izawa, N.; Matsumoto, T.; Imai, Y.; Seo, S.; Kurokawa, M.; et al. Identification of Nedd9 as a TGF-β-Smad2/3 Target Gene Involved in RANKL-Induced Osteoclastogenesis by Comprehensive Analysis. PLoS ONE 2016, 11, e0157992. [Google Scholar] [CrossRef] [PubMed]
  51. Ignotz, R.A.; Massagué, J. Type beta transforming growth factor controls the adipogenic differentiation of 3T3 fibroblasts. Proc. Natl. Acad. Sci. USA 1985, 82, 8530–8534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Rebbapragada, A.; Benchabane, H.; Wrana, J.L.; Celeste, A.J.; Attisano, L. Myostatin Signals through a Transforming Growth Factor β-Like Signaling Pathway To Block Adipogenesis. Mol. Cell. Biol. 2003, 23, 7230–7242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Bangs, F.K.; Schrode, N.; Hadjantonakis, A.-K.; Anderson, K.V. Lineage specificity of primary cilia in the mouse embryo. Nat. Cell Biol. 2015, 17, 113–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Arrighi, N.; Lypovetska, K.; Moratal, C.; Giorgetti-Peraldi, S.; Dechesne, C.A.; Dani, C.; Peraldi, P. The primary cilium is necessary for the differentiation and the maintenance of human adipose progenitors into myofibroblasts. Sci. Rep. 2017, 7, 15248. [Google Scholar] [CrossRef] [PubMed]
  55. Zappaterra, M.; Gioiosa, S.; Chillemi, G.; Zambonelli, P.; Davoli, R. Muscle transcriptome analysis identifies genes involved in ciliogenesis and the molecular cascade associated with intramuscular fat content in Large White heavy pigs. PLoS ONE 2020, 15, e0233372. [Google Scholar] [CrossRef]
  56. Shagisultanova, E.; Gaponova, A.V.; Gabbasov, R.; Nicolas, E.; Golemis, E.A. Preclinical and clinical studies of the NEDD9 scaffold protein in cancer and other diseases. Gene 2015, 567, 1–11. [Google Scholar] [CrossRef] [Green Version]
  57. Xu, M.; Chen, X.; Chen, D.; Yu, B.; Huang, Z. FoxO1: A novel insight into its molecular mechanisms in the regulation of skeletal muscle differentiation and fiber type specification. Oncotarget 2016, 8, 10662–10674. [Google Scholar] [CrossRef] [Green Version]
  58. Mastrangelo, S.; Jemaa, S.B.; Ciani, E.; Sottile, G.; Moscarelli, A.; Boussaha, M.; Montedoro, M.; Pilla, F.; Cassandro, M. Genome-wide detection of signatures of selection in three Valdostana cattle populations. J. Anim. Breed. Genet. 2020, 137, 609–621. [Google Scholar] [CrossRef]
  59. rs343048625 (SNP)—Explore This Variant—Sus_scrofa—Ensembl Genome Browser 101. Available online: http://uswest.ensembl.org/Sus_scrofa/Variation/Explore?v=rs343048625;vdb=variation (accessed on 12 October 2020).
  60. England, S.J.; Cerda, G.A.; Kowalchuk, A.; Sorice, T.; Grieb, G.; Lewis, K.E. Hmx3a does not require its homeodomain for its essential functions in spinal cord, ear and lateral line development. bioRxiv 2020, 68. [Google Scholar] [CrossRef] [Green Version]
  61. Sandhu, J.; Li, S.; Fairall, L.; Pfisterer, S.G.; Gurnett, J.E.; Xiao, X.; Weston, T.A.; Vashi, D.; Ferrari, A.; Orozco, J.L.; et al. Aster Proteins Facilitate Nonvesicular Plasma Membrane to ER Cholesterol Transport in Mammalian Cells. Cell 2018, 175, 514–529.e20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Kutyavin, V.I.; Chawla, A. Aster: A New Star in Cholesterol Trafficking. Cell 2018, 175, 307–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Farber, C.; Aten, J.; Farber, E.; Vera, V.; Gularte-Mérida, R.; Islas-Trejo, A.; Wen, P.; Horvath, S.; Lucero, M.; Lusis, A.; et al. Genetic dissection of a major mouse obesity QTL (Carfhg2): Integration of gene expression and causality modeling. Physiol. Genomics 2009, 37, 294–302. [Google Scholar] [CrossRef] [Green Version]
  64. Serão, N.V.; González-Peña, D.; Beever, J.E.; Faulkner, D.B.; Southey, B.R.; Rodriguez-Zas, S.L. Single nucleotide polymorphisms and haplotypes associated with feed efficiency in beef cattle. BMC Genet. 2013, 14, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Boxplots of the lean content (%) estimated by Ham InspectorTM for the green hams grouped according to the EUROP classes (U, R, and O) of the corresponding carcasses. For each EUROP class, the range of % estimated lean carcass is reported. Blank circles represent outliers.
Figure 1. Boxplots of the lean content (%) estimated by Ham InspectorTM for the green hams grouped according to the EUROP classes (U, R, and O) of the corresponding carcasses. For each EUROP class, the range of % estimated lean carcass is reported. Blank circles represent outliers.
Animals 11 00068 g001
Table 1. Effect of the processing plant (PP) and sex (S) of the animals on the traits (Estimated Marginal Means ± standard error) of hams (green and after salting steps).
Table 1. Effect of the processing plant (PP) and sex (S) of the animals on the traits (Estimated Marginal Means ± standard error) of hams (green and after salting steps).
Ham TraitsProcessing Plant (PP) Sex (S) PPxS
123p-ValueBarrowFemalep-Value
n. 60n.120n.50n. 105n. 125p-Value
pHu5.63c ± 0.025.69 b ± 0.015.79 a ± 0.02<0.0015.71 ± 0.015.70 ± 0.01n.s.n.s.
Weight GH, kg 115.10a ± 0.1013.90 b ± 0.1012.60 c ± 0.20<0.00113.90 ± 0.1113.80 ± 0.10n.s.n.s.
Lean GH, % 262.70a ± 0.3062.90 a ± 0.2060.60 b ± 0.30<0.00161.10 ± 0.2062.90 ± 0.20<0.001n.s.
Weight 1S, kg 314.90a ± 0.1013.70 b ± 0.1012.50 c ± 0.20<0.00113.70 ± 0.1013.60 ± 0.10n.s.n.s.
Weight Loss 1S, % 41.31 ± 0.031.26 ± 0.021.36 ± 0.04n.s.1.28 ± 0.031.35 ± 0.02<0.05n.s.
Salt 1S, % 51.06 c ± 0.011.21 b ± 0.011.46 a ± 0.02<0.0011.22 ± 0.011.27 ± 0.01<0.001n.s.
Weight ES, kg 614.60 a ± 0.1013.50 b ± 0.1012.20 c ± 0.20<0.00113.70 ± 0.1013.40 ± 0.10n.s.n.s.
Weight Loss ES, % 73.13 a ± 0.062.64 b ± 0.042.72 b ± 0.07<0.0012.73 ± 0.052.93 ± 0.04<0.01n.s.
Salt ES, % 82.37 c ± 0.022.80 a ± 0.022.52 b ± 0.03<0.0012.53 ± 0.022.60 ± 0.02<0.01<0.05
Salt 1S /Salt ES, % 944.70 b ± 0.3043.00 c ± 0.20 57.50 a ± 0.40<0.00148.50 ± 0.3048.30 ± 0.20n.s.n.s.
Estimated Marginal Means with different letters along rows are significantly different (LSD posthoc test); n.s: not significant. 1 Weight of green hams measured with Ham InspectorTM and expressed in kg. 2 Lean content of green hams estimated by Ham InspectorTM and expressed as a percentage of green ham weight (%). 3 Weight of hams at 1st salting measured with Ham InspectorTM and expressed in kg. 4 Weight loss measured after 1st salting and expressed as percentage loss of green ham weight (%). 5 Salt (as NaCl) content of the lean part of salted hams at 1st salting, estimated by Ham InspectorTM and expressed as a percentage of wet weight (%). 6 Weight of hams measured at the end of salting with Ham InspectorTM and expressed in kg. 7 Weight loss measured at the end of salting and expressed as percentage loss of green ham weight (%). 8 Salt (as NaCl) content of the lean part of salted hams at the end of salting, estimated by Ham InspectorTM and expressed as a percentage of wet weight (%). 9 Ratio between salt at 1st salting step and salt content at the end of salting, expressed as percentage (%).
Table 2. The list of the significant markers, with the associated traits, the marker position, the type of variant, and the list of the candidate genes in the regions flanking the significant markers. Chromosome-wide significant markers are in bold.
Table 2. The list of the significant markers, with the associated traits, the marker position, the type of variant, and the list of the candidate genes in the regions flanking the significant markers. Chromosome-wide significant markers are in bold.
TraitMarkerMarker Rs CodePosition 1MAF 2Pc1df 3Type of VariantCandidate Genes in the Region 4
pHuWU_10.2_18_17949287rs32131741418:17,103,7850.222.54 × 10−6intergenic variantPLXNA4, MKLN1
WU_10.2_4_91195648rs3429769524:83,519,8690.133.94 × 10−5intron variant of the gene CD247DCAF6, MPC2, ADCY10, MPZL1, RCSD1, CREG1, CD247, POU2F1, DUSP27, GPA33, MAEL
Lean GH, % 5ASGA0016987rs809945544:2,020,9900.251.42 × 10−5exon variant of a non-coding transcriptADGRB1, MROH5, PTP4A3, GPR20, SLC45A4
ALGA0002237rs808831861:30,193,3130.163.07 × 10−5intergenic variantSGK1, SLC2A12, TBPL1, TCF21, EYA4
H3GA0000815rs808489051:11,662,8200.354.72 × 10−5intron variant of the gene TIAM2NOX3, TFB1M, TIAM2, SCAF8
ASGA0001055rs809238301:11,684,4500.334.72 × 10−5intron variant of the gene TIAM2NOX3, TFB1M, TIAM2, SCAF8
Weight GH, kg 6CASI0010463rs33563591315:12,997,3620.135.48 × 10−6intron variant of the gene NXPH2NXPH2, SPOPL
WU_10.2_14_144250775rs34304862514:132,664,2620.401.33 × 10−5intergenic variantAWN, PSP-II, SPMI, PSTK, IKZF5, ACADSB, HMX2, HMX3, BUB3
Weight 1S, kg 7CASI0010463rs33563591315:12,997,3620.134.13 × 10−6intron variant of the gene NXPH2NXPH2, SPOPL
WU_10.2_14_144250775rs34304862514:132,664,2620.401.26 × 10−5intergenic variantAWN, PSP-II, SPMI, PSTK, IKZF5, ACADSB, HMX2, HMX3, BUB3
ASGA0026341rs808598295: 74,704,4570.334.57 × 10−5intergenic variantADAMTS20, PUS7L, IRAK4, TWF1, U6, TMEM117
Weight ES, kg 8CASI0010463rs33563591315:12,997,3620.132.15 × 10−6intron variant of the gene NXPH2NXPH2, SPOPL
WU_10.2_14_144250775rs34304862514:132,664,2620.404.27 × 10−6intergenic variantAWN, PSP-II, SPMI, PSTK, IKZF5, ACADSB, HMX2, HMX3, BUB3
ASGA0026341rs808598295: 74,704,4570.332.23 × 10−5intergenic variantADAMTS20, PUS7L, IRAK4, TWF1, U6, TMEM117
WU_10.2_7_118557013rs3258878617:111,991,7730.124.25 × 10−5intergenic variantEFCAB11, TDP1, KCNK13, PSMC1, NRDE2, CALM1
ALGA0044906rs808869097:112,364,4050.084.66 × 10−5intron variant of the gene TTC7BPSMC1, NRDE2, CALM1, TTC7B, RPS6K, A5TTC7B
Weight Loss 1S, % 9ASGA0031014rs809633187:8,081,7340.471.00 × 10−5intron variant of the gene NEDD9GCM2, ELOVL2, SMIM13, NEDD9, TMEM170B, ADTRP, HIVEP1
WU_10.2_10_74620421-10:67,919,7110.251.57 × 10−5sequence variant of a non-coding geneADARB2
Weight Loss ES, % 10ALGA0022599rs809171914:6,665,8560.272.90 × 10−5intergenic variantKHDRBS3
Salt 1S, % 11WU_10.2_14_36295226rs32387915414:34,218,1840.441.71 × 10−5intron variant of the genes SUDS3 and SRRM4SUDS3, TAOK3, VSIG10, WSB2, RFC5, KSR2SUDS3, SRRM4
ASGA0000817rs809212161:8,287,0080.422.22 × 10−5intergenic variantFNDC1, TAGAP, RSPH3, EZR, SYTL3, DYNLT1, TMEM181
Salt ES, % 12INRA0000796rs3324908621:13,370,6390.084.36 × 10−6intron variant of the gene RGS17RGS17,MTRFL1, FBXO5, VIP, MYCT1
ASGA0102337rs813236319:50,457,8990.351.46 × 10−5intron variant of the gene GRAMD1BHSPA8, CLMP, GRAMD1B, SCN3B, ZNF202, OR6X1, OR6M1, OR4D5, OR6T1,
1 The marker position in Sus scrofa Genome Assembly Build 11.1 is reported as Sus scrofa chromosome: nucleotide position in base pairs. 2 Minor Allele Frequency. 3 Significance estimates obtained with GenABEL with the correction for inflation factor (stratification effects). Bold values indicate markers significant at the chromosome-wide level. 4 The flanking region is the genomic region located between the marker position ± 500,000 nucleotides. 5 Lean content of green hams estimated by Ham InspectorTM and expressed as a percentage of green ham weight (%). 6 Weight of green hams measured with Ham InspectorTM and expressed in kg. 7 Weight of hams at 1st salting measured with Ham InspectorTM and expressed in kg. 8 Weight of hams at the end of salting measured with Ham InspectorTM and expressed in kg. 9 Weight loss measured after 1st salting and expressed as percentage loss of green ham weight (%). 10 Weight loss measured at the end of salting and expressed as percentage loss of green ham weight (%). 11 Salt (as NaCl) content of the lean part of salted hams at 1st salting, estimated by Ham InspectorTM and expressed as percentage of wet weight (%). 12 Salt (as NaCl) content of the lean part of salted hams at the end of salting, estimated by Ham InspectorTM and expressed as a percentage of wet weight (%).
Table 3. Estimated Least Squares Means (LSM) ± standard errors (S.E.), additive and dominance genetic effects of the markers associated with ham traits. Between brackets are reported the numbers of the observations used for each genotype.
Table 3. Estimated Least Squares Means (LSM) ± standard errors (S.E.), additive and dominance genetic effects of the markers associated with ham traits. Between brackets are reported the numbers of the observations used for each genotype.
TraitMarkerAlleleLSM ± S.E (N)Additive EffectDominance Effect
12111222
pHuWU_10.2_18_17949287GA5.71 a ± 0.045.73 a ± 0.025.70 a ± 0.01n.s.n.s.
(10)(55)(104)
WU_10.2_4_91195648GA-5.69 a ± 0.025.72 a ± 0.01--
(0)(43)(126)
Lean GH, % 1ASGA0016987GA63.91 a ± 0.8161.91 b ± 0.2762.61 ab ± 0.23n.s.0.042
(8)(64)(93)
ALGA0002237CA62.15 a ± 0.2163.04 a ± 0.3362.47 a ± 1.10n.s.n.s.
(115)(46)(4)
H3GA0000815GA63.57 a ± 0.4862.71 a ± 0.2561.75 b ± 0.260.0005n.s.
(20)(75)(70)
ASGA0001055GA61.72 a ± 0.5162.61 a ± 0.2662.40 a ± 0.25n.s.n.s.
(19)(70)(76)
Weight GH, kg 2CASI0010463CA13.58 ab ± 0.6814.05 a ± 0.1613.66 b ± 0.09n.s.n.s.
(2)(38)(125)
WU_10.2_14_144250775GA13.93 ab ± 0.2013.88 a ± 0.1013.48 b ± 0.13n.s.n.s.
(24)(85)(56)
Weight 1S, kg 3CASI0010463CA13.42 ab ± 0.6913.85 a ± 0.1613.48 b ± 0.09n.s.n.s.
(2)(37)(123)
WU_10.2_14_144250775GA13.74 ab ± 0.2013.69 a ± 0.1113.29 b ± 0.13n.s.n.s.
(24)(83)(55)
ASGA0026341GA14.01 a ± 0.2413.71 a ± 0.1113.29 b ± 0.120.009n.s.
(16)(78)(68)
Weight ES, kg 4CASI0010463CA13.12 ab ± 0.6913.63 a ± 0.1713.24 b ± 0.09n.s.n.s.
(2)(37)(116)
WU_10.2_14_144250775GA13.48 ab ± 0.2113.48 a ± 0.1113.04 b ± 0.13n.s.n.s.
(23)(79)(52)
ASGA0026341GA13.70 a ± 0.2513.50 a ± 0.1113.05 b ± 0.120.022n.s.
(16)(78)(68)
WU_10.2_7_118557013GA13.18 c ± 0.0913.74 b ± 0.1715.19 a ± 0.670.003n.s.
(120)(33)(2)
ALGA0044906GA13.24 b ± 0.0813.81 a ± 0.20---
(131)(24)(0)
Weight Loss 1S, % 5ASGA0031014GA 1.23 b ± 0.04 (42)n.s.
WU_10.2_10_74620421GA1.24 b ± 0.021.37 a ± 0.031.35 ab ± 0.08n.s.n.s.
(91)(61)(10)
Weight Loss ES, % 6ALGA0022599GA2.82 a ± 0.042.81 a ± 0.052.90 a ± 0.11n.s.n.s.
(86)(57)(12)
Salt 1S, % 7WU_10.2_14_36295226GA1.26 a ± 0.021.25 ab ± 0.011.21 b ± 0.010.036n.s.
(28)(85)(49)
ASGA0000817GA1.25 a ± 0.011.23 a ± 0.011.24 a ± 0.02n.s.n.s.
(51)(87)(24)
Salt ES, % 8INRA0000796CA-2.62 a ± 0.032.63 a ± 0.01--
(0)(21)(134)
ASGA0102337GA2.62 a ± 0.042.62 a ± 0.022.64 a ± 0.02n.s.n.s.
(17)(73)(65)
LSM with different letters along rows are significantly different for p < 0.10 adjusted with the Tukey test. “-” means that additive and dominance effect cannot be estimated for that marker as two genotypes were found in the samples. n.s. stands for not significant. 1 Lean content of green hams estimated by Ham InspectorTM and expressed as a percentage of green ham weight (%). 2 Weight of green hams measured with Ham InspectorTM and expressed in kg. 3 Weight of hams at 1st salting measured with Ham InspectorTM and expressed in kg. 4 Weight of hams at the end of salting measured with Ham InspectorTM and expressed in kg. 5 Weight loss measured after 1st salting and expressed as percentage loss of green ham weight (%). 6 Weight loss measured at the end of salting and expressed as percentage loss of green ham weight (%). 7 Salt (as NaCl) content of the lean part of salted hams at 1st salting, estimated by Ham InspectorTM and expressed as a percentage of wet weight (%). 8 Salt (as NaCl) content of the lean part of salted hams at the end of salting, estimated by Ham InspectorTM and expressed as a percentage of wet weight (%).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zappaterra, M.; Zambonelli, P.; Schivazappa, C.; Simoncini, N.; Virgili, R.; Stefanon, B.; Davoli, R. Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs. Animals 2021, 11, 68. https://doi.org/10.3390/ani11010068

AMA Style

Zappaterra M, Zambonelli P, Schivazappa C, Simoncini N, Virgili R, Stefanon B, Davoli R. Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs. Animals. 2021; 11(1):68. https://doi.org/10.3390/ani11010068

Chicago/Turabian Style

Zappaterra, Martina, Paolo Zambonelli, Cristina Schivazappa, Nicoletta Simoncini, Roberta Virgili, Bruno Stefanon, and Roberta Davoli. 2021. "Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs" Animals 11, no. 1: 68. https://doi.org/10.3390/ani11010068

APA Style

Zappaterra, M., Zambonelli, P., Schivazappa, C., Simoncini, N., Virgili, R., Stefanon, B., & Davoli, R. (2021). Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs. Animals, 11(1), 68. https://doi.org/10.3390/ani11010068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop