Copy number variants in Italian Large White pigs detected using high‐density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance‐tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with penncnv software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by penncnv were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.     

Next Generation Semiconductor Based-Sequencing of a Nutrigenetics Target Gene (GPR120) and Association with Growth Rate in Italian Large White Pigs

The GPR120 gene (also known as FFAR4 or O3FAR1) encodes for a functional omega-3 fatty acid receptor/sensor that mediates potent insulin sensitizing effects by repressing macrophage-induced tissue inflammation. For its functional role, GPR120 could be considered a potential target gene in animal nutrigenetics. In this work we resequenced the porcine GPR120 gene by high throughput Ion Torrent semiconductor sequencing of amplified fragments obtained from 8 DNA pools derived, on the whole, from 153 pigs of different breeds/populations (two Italian Large White pools, Italian Duroc, Italian Landrace, Casertana, Pietrain, Meishan, and wild boars). Three single nucleotide polymorphisms (SNPs), two synonymous substitutions and one in the putative 3′-untranslated region (g.114765469C > T), were identified and their allele frequencies were estimated by sequencing reads count. The g.114765469C > T SNP was also genotyped by PCR-RFLP confirming estimated frequency in Italian Large White pools. Then, this SNP was analyzed in two Italian Large White cohorts using a selective genotyping approach based on extreme and divergent pigs for back fat thickness (BFT) estimated breeding value (EBV) and average daily gain (ADG) EBV. Significant differences of allele and genotype frequencies distribution was observed between the extreme ADG-EBV groups (P < 0.001) whereas this marker was not associated with BFT-EBV.     

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

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

A single nucleotide polymorphism in the porcine cathepsin K (<Emphasis Type="Italic">CTSK</Emphasis>) gene is associated with back fat thickness and production traits in Italian Duroc pigs

Cathepsin K (CTSK) was selected as a candidate gene for fat deposition in pigs because recently, in human and mouse, it was shown that this lysosomal proteinase is an obesity marker. A single nucleotide polymorphism (SNP) was identified in intron 4 of the porcine CTSK gene (g.15G>A; FM209043). Allele frequencies of this polymorphism were analysed in seven pig breeds. Radiation hybrid mapping confirmed the localization of CTSK to porcine chromosome 4, close to the FAT1 QTL region. Three populations of pigs (one Italian Large White and two Italian Duroc groups of pigs) were selected for association analysis. In the Italian Large White breed the g.15G>A SNP was not informative. Association analysis including all Italian Duroc pigs showed that the CTSK marker was associated with back fat thickness and lean cuts (P < 0.01), and average daily gain and feed:gain ratio (P < 0.05) estimated breeding values.     

The porcine <Emphasis Type="Italic">TBC1D1</Emphasis> gene: mapping,SNP identification,and association study with meat,carcass and production traits in Italian heavy pigs

TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) domain family, member 1] is a Rab-GTPase-activating related protein implicated in regulating the trafficking of glucose transporter 4 (GLUT4 or SLC2A4) storage vesicles to the cell surface in response to insulin and AMPK-activating stimuli in skeletal muscle. Mutations in the human and mouse TBC1D1 genes confer risk of obesity or leanness. We identified five single nucleotide polymorphisms (SNPs) in the porcine TBC1D1 gene. One of them (FN677935:g.219G>A) was genotyped either by high resolution melting and PCR-RFLP analyses to study allele frequencies in a few pig breeds and evaluate association with meat production and carcass traits in five groups of sib-tested pigs of Italian Large White and Italian Duroc breeds. The g.219G>A SNP was associated (P < 0.05) with ham weight, back fat thickness and lean cuts content in Italian Large White and with visible intermuscular fat in Italian Duroc pigs.     

SNPs of MYPN and TTN genes are associated to meat and carcass traits in Italian Large White and Italian Duroc pigs

Myopalladin (MYPN) and titin (TTN) play key roles in skeletal muscle structural organisation and for this reason they may be functional candidates for meat production traits in pigs. We analysed by PCR–RFLP two polymorphisms previously identified and described in the 3′UTR of MYPN and TTN genes in a group of Italian Large White (ILW) and Italian Duroc (ID) pigs with extreme divergent estimated breeding values for some production traits. Significant allele frequency differences were observed for lean cuts (LC) at the MYPN locus in ILW breed and for visible intermuscular fat (VIF) at MYPN and TTN loci in ID breed. Moreover, association analysis between the two genes and some production traits (average daily gain, LC, backfat thickness, ham weight, feed conversion ratio and VIF) was performed in a sample of 409 ILW and 381 ID pigs. On the whole, the results showed significant effects of both genes on some of the considered traits in ID breed, and for TTN gene in ILW animals. These results indicate that both genes can be considered as candidates for meat and carcass traits in ID pigs.     

Polymorphisms in an obesity-related gene (PCSK1) are associated with fat deposition and production traits in Italian heavy pigs

The proprotein convertase subtilisin/kexin type 1 (PCSK1) gene encodes the prohormone convertase 1/3 enzyme that processes prohormones into functional hormones that, in turn, regulate central and peripheral energy metabolism. Mutations in the human PCSK1 gene cause severe monogenic obesity or confer risk of obesity. We herein investigated the porcine PCSK1 gene with the aim of identifying polymorphisms associated with fat deposition and production traits in Italian heavy pigs. By re-sequencing about 5.1 kb of this gene in 21 pigs of different breeds, we discovered 14 polymorphisms that were organized in nine haplotypes, clearly distributed in two clades of putative European and Asian origin. Then we re-mapped this gene on porcine chromosome 2 and analysed its expression in several tissues including gastric oxyntic mucosa of weanling pigs in which PCSK1 processes the pre-pro-ghrelin into ghrelin, which in turn is involved in the control of feed intake and energy metabolism. Association analyses between PCSK1 single-nucleotide polymorphisms (SNPs) and production, carcass and several other traits were conducted on five groups of pigs from three different experimental designs, for a total of 1221 animals. Results indicated that the analysed SNPs were associated (P < 0.01 or P < 0.05) with several traits including backfat thickness and visible intermuscular fat in Italian Duroc (ID) and growth performances in Italian Large White (ILW) and in ILW × Italian Landrace pigs. However, the effects estimated in the ILW were opposite to the effects reported in the ID pigs. Suggestive association (P < 0.10) was observed with muscle cathepsin B activity, opening, if confirmed, potential applications to reduce the excessive softness defect of the green hams that is of particular concern for the processing industry. The results obtained supported the need to further investigate the PCSK1 gene to fully exploit the value of its variability and apply this information in pig breeding programmes.     

Mapping of quantitative trait loci on porcine chromosome 4

A F₂ population derived from a cross between European Large White and Chinese Meishan pigs was established in order to study the genetic basis of breed differences for growth and fat traits. Chromosome 4 was chosen for initial study as previous work had revealed quantitative trait loci (QTLs) on this chromosome affected growth and fat traits in a Wild Boar × Large White cross. Individuals in the F₂ population were typed for nine markers spanning a region of approximately 124 c m . We found evidence for QTLs affecting growth between weaning and the end of test (additive effect: 43·4 g/day) and fat depth measured in the mid-back position (additive effect: 1·82 mm). There was no evidence of interactions between the QTLs and sex, grandparents or F₁ sires, suggesting that the detected QTLs were fixed for alternative alleles in the Meishan and Large White breeds. Comparison of locations suggests that these QTLs could be the same as those found in the Wild Boar × Large White cross.     

Association between polymorphisms in cathepsin and cystatin genes with meat production and carcass traits in Italian Duroc pigs: confirmation of the effects of a cathepsin L (CTSL) gene marker

We genotyped single nucleotide polymorphisms (SNPs) in 5 cathepsin or cystatin genes (cathepsin F, CTSF; cathepsin L, CTSL; cathepsin S, CTSS; cathepsin Z, CTSZ; cystatin B, CSTB) in two groups of Italian Duroc pigs: the first group (n. 100) was chosen using a selective genotyping approach with extreme estimated breeding value (EBV) for visible intermuscular fat (VIF); the second group (n. 218) was made of performance-tested Duroc pigs not selected by any criteria. CTSL marker showed a tendency towards association (P < 0.10) with VIF (first group) and back fat thickness (BFT) and average daily gain (ADG; second group). In the second group, the CTSL polymorphism was associated with weight of lean cuts (LC; P < 0.05). Additive effects for all mentioned traits in the second group was significant (P < 0.05). The results we obtained in the Italian Duroc pigs confirmed the results and the direction of the effects already reported for the Italian Large White breed.     

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

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

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool‐seq approach

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool‐seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5–10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.     

Development of intra- and intermuscular adipose tissue in growing large white and Meishan pigs

Lipogenic enzyme activities of porcine intra- and intermuscular adipose tissues were determined in growing lean (Large White) and fat (Meishan) pigs. The activities of acetyl-CoA-carboxylase (ACX), malic enzyme (ME) and glucose-6-phosphate dehydrogenase (G6PDH) were compared in both breeds and at both adipose sites. All three enzyme activities were much lower in the intramuscular adipose tissue than in the intermuscular site. Although the lipogenic activity of the intramuscular adipose site was low, it appeared, however, to possess adequate levels of enzymes for in situ lipid synthesis. The highest differences in lipogenic enzyme activities between Meishan and Large White pigs were found in intramuscular adipose tissue, and essentially concerned the activity of malic enzyme which was much higher in Meishan pigs. A close relationship between ME activity and lipid content of intramuscular adipose tissue was observed in both breeds. It was concluded that ME appeared to be a major factor affecting the incidence of higher intramuscular fat in the pig.     

Metabolomics evidences plasma and serum biomarkers differentiating two heavy pig breeds

In pigs, many production traits are known to vary among breeds or lines. These traits can be considered end phenotypes or external traits as they are the final results of complex biological interactions and processes whose fine biological mechanisms are still largely unknown. This study was designed to compare plasma and serum metabolomic profiles between animals of two heavy pig breeds (12 Italian Large White and 12 Italian Duroc), testing indirectly the hypothesis that different genetic backgrounds might be the determining factors of differences observed on the level of metabolites in the analyzed biofluids between breeds. We used a targeted metabolomic approach based on mass spectrometric detection of about 180 metabolites and applied a statistical validation pipeline to identify differences in the metabolomic profiles of the two heavy pig breeds. Blood samples were collected after jugulation at the slaughterhouse and prepared for metabolomics analysis that was carried out using the Biocrates AbsoluteIDQ p180 Kit, covering five different biochemical classes: glycerophospholipids, amino acids, biogenic amines, hexoses and acylcarnitines. A statistical pipeline that included the selection of the most relevant metabolites differentiating the two breeds by sparse Partial Least Squares Discriminant Analysis (sPLS-DA) was coupled with a stability test and significance test determined with leave one out and permutation procedures. sPLS-DA plots clearly separated the pigs of the two investigated breeds. A few metabolites (a total of five metabolites considering the two biofluids) involved in key metabolic pathways largely contributed to these differences between breeds. In particular, a higher level of the sphingomyelins SM (OH) C14:1 (both in plasma and serum), SM (OH) C16:1 (in serum) and SM C16:0 (in serum) were observed in Italian Duroc than in Italian Large White pigs and the inverse was for the biogenic amine kynurenine (in plasma). The level of another biogenic amine (acetylornithine) was higher in Italian Large White than in Italian Duroc pigs in both analysed biofluids. These results provided biomarkers that could be important to understand the biological differences between these two heavy pig breeds. In particular, according to the functional role played by sphingomyelins in obesity-induced inflammatory responses, it could be possible to speculate that a higher level of sphingomyelins in Italian Duroc might be related to the higher interrmuscular fat deposition of this breed compared with the Italian Large White. Additional studies will be needed to evaluate the relevance of these biomarkers for practical applications in pig breeding and nutrition.     

Analysis of Association Between the MUC4 g.8227C>G Polymorphism and Production Traits in Italian Heavy Pigs Using a Selective Genotyping Approach

In pigs, susceptibility to enterotoxigenic Escherichia coli (ETEC) K88 strains (locus F4bcR) is determined by a dominant allele, with the recessive allele determining resistance. The susceptible allele also appeared to be associated with higher growth rate even with discordant results. A single nucleotide polymorphism (SNP) in exon 7 of the mucin 4 (MUC4) gene (DQ848681:g.8227C>G), shown to be in close linkage disequilibrium with the F4bcR locus, has been used as marker to identify susceptible pigs, substituting invasive villous adhesion tests. We herein analyzed this SNP in Italian local breeds and applied a selective genotyping approach in Italian Large White, Italian Landrace, and Italian Duroc comparing allele frequency distribution in groups of pigs with extreme estimated breeding values (EBV) for average daily gain (ADG) and backfat thickness (BFT) to evaluate if this marker is associated with these traits. Allele G (associated with susceptibility to ETEC) was associated with higher ADG and BFT in Italian Large White (P?=?6.66E-04 and P?=?0.012, respectively) and higher ADG in Italian Landrace (P?=?7.23E-12). This polymorphism was poorly informative in Italian Duroc. Antagonistic associations of the MUC4?g.8227C>G alleles on susceptibility to ETEC and growth performances evidence the complexity of applying marker assisted selection in pig breeding.     

Analysis of Association Between the MUC4 g.8227C>G Polymorphism and Production Traits in Italian Heavy Pigs Using a Selective Genotyping Approach

In pigs, susceptibility to enterotoxigenic Escherichia coli (ETEC) K88 strains (locus F4bcR) is determined by a dominant allele, with the recessive allele determining resistance. The susceptible allele also appeared to be associated with higher growth rate even with discordant results. A single nucleotide polymorphism (SNP) in exon 7 of the mucin 4 (MUC4) gene (DQ848681:g.8227C>G), shown to be in close linkage disequilibrium with the F4bcR locus, has been used as marker to identify susceptible pigs, substituting invasive villous adhesion tests. We herein analyzed this SNP in Italian local breeds and applied a selective genotyping approach in Italian Large White, Italian Landrace, and Italian Duroc comparing allele frequency distribution in groups of pigs with extreme estimated breeding values (EBV) for average daily gain (ADG) and backfat thickness (BFT) to evaluate if this marker is associated with these traits. Allele G (associated with susceptibility to ETEC) was associated with higher ADG and BFT in Italian Large White (P = 6.66E-04 and P = 0.012, respectively) and higher ADG in Italian Landrace (P = 7.23E-12). This polymorphism was poorly informative in Italian Duroc. Antagonistic associations of the MUC4 g.8227C>G alleles on susceptibility to ETEC and growth performances evidence the complexity of applying marker assisted selection in pig breeding.     

Confirmed association between a single nucleotide polymorphism in the <Emphasis Type="Italic">FTO</Emphasis> gene and obesity-related traits in heavy pigs

We recently showed that a polymorphism in the fat mass and obesity associated (FTO) gene (AM931150: g.276T > G) is associated with fat deposition traits in pigs. To confirm this result, we genotyped this polymorphism in an Italian Duroc population made up by 313 performance tested pigs with known estimated breeding values (EBVs) for average daily gain, back fat thickness (BFT), feed:gain ratio, lean cuts (LC), and visible intermuscular fat (VIF, a measure of intermuscular fat in the hams). In addition, we genotyped 148 commercial heavy pigs for which several fat deposition traits and lean meat percentage were measured. The results of the association analyses confirmed the effect of the FTO mutation on obesity-related traits (VIF, BFT and LC) in the Italian Duroc pigs (P < 0.01) and in the commercial pigs (intramuscular fat content of different muscles, P < 0.05 or P < 0.10; lean meat content, P < 0.05; BFT, P < 0.05; intermuscular fat content in the hams, P < 0.05).     

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

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

A study of associations of the H-FABP genotypes with fat and meat production of pigs

Although the content of intramuscular fat (IMF) influences significantly meat quality, it can be estimated only after the slaughter of animals. Variants of the H-FABP gene were suggested as candidate genes influencing the variability of IMF. The effect of H-FABP - HinfI polymorphism on the content of IMF, backfat thickness, the weight and percentages of major meat parts and of the leg in carcass weight was studied in a group of 97 pigs (46 gilts and 51 barrows) of Large White and Landrace breeds using the test of fattening capacity and carcass value. In the set of experimental animals, the frequencies of genes were H = 0.75 +/- 0.03 and h = 0.25 +/- 0.03. Biometric analyses did not corroborate differences among different H-FABP - HinfI genotypes and all the traits under study. Only in genotypes HH and Hh the differences between least-square means of phenotypic IMF values under study were close to the limit of significance (P = 0.06).     

Gene frequencies of blood groups in Dutch swine breeds

Gene frequencies of blood groups at 6 loci were estimated within two Dutch breeds of swine, Dutch Landrace and Dutch Large White. At all loci there were significant differences between breeds. The relationship of the two breeds with, respectively, the German Landrace breed and the German Large White breed could be confirmed by the gene frequencies.