A single nucleotide polymorphism in the porcine cathepsin K (CTSK) gene is associated with back fat thickness and production traits in Italian Duroc pigs

A SNP in the ovine cathepsin K (CTSK) gene is associated with yearling growth performance in a crossbred sheep population

Cathepsin K (CTSK) is a lysosomal protease existent in the skeletal muscles which is involved in biochemical processes related to obesity. Several studies have reported the effects of CTSK gene on body weight and fat deposition in human, mice and pigs. However, information about its structure and functions in sheep is very limited. Thus, this study was performed to evaluate the association between CTSK gene variants and yearling growth performance in Afshari × Booroola-Merino crossbred sheep. A fragment of 500 bp in exon 6 and partial of intron 5 of CTSK gene was amplified with polymerase chain reaction (PCR). All animals were genotyped by single-stranded conformation polymorphism (SSCP) and further confirmed by sequencing. Association analysis using a fixed linear model indicated that g.106510225G > A SNP was significantly related to average daily weight gain (ADWG) per year, fat-tail weight to carcass weight ratio (FW/CW), muscle thickness (MT) and muscle cross-sectional area (MCSA) of animals (p < 0.05). Due to the low polymorphic information content (PIC <0.25) for targeted locus in studied population, more association studies are needed to confirm the CTSK gene effects on growth traits in sheep.

Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits

Maternal genetic effects (MGE) could affect meat quality traits such as intramuscular fat (IMF) and its fatty acid composition. However, it has been scarcely studied, especially in rabbits. The objectives of the present study were, first, to assess the importance of MGE on intramuscular fat and fatty acid composition by applying a Bayesian maternal animal model in two rabbit lines divergently selected for IMF. The second objective was to identify genomic regions and candidate genes of MGE that are associated with the traits of these offspring, using Bayesian methods in a Genome Wide Association Study (GWAS). Quantitative analyses were performed using data from 1982 rabbits, and 349 animals from the 9th generation and 76 dams of the 8th generation with 88,512 SNPs were used for the GWAS. The studied traits were IMF, saturated fatty acids (total SFA, C14:0; myristic acid, C16:0; palmitic acid and C18:0; stearic acid), monounsaturated fatty acids (total MUFA, C16:1n-7; palmitoleic acid and C18:1n-9; oleic acid), polyunsaturated fatty acids (total PUFA, C18:2n-6; linoleic acid, C18:3n-3; α-linolenic acid and C20:4n-6; arachidonic acid), MUFA/SFA and PUFA/SFA. The proportion of phenotypic variance explained by the maternal genetic effect ranged from 8 to 22% for IMF, depending on the model. For fatty acid composition, the proportion of phenotypic variance explained by maternal genetic effects varied from 10% (C18:0) to 46% (MUFA) in a model including both direct and additive maternal genetic effects, together with the common litter effect as a random variable. In particular, there were significant direct maternal genetic correlations for C16:0, C18:1n9, C18:2n6, SFA, MUFA, and PUFA with values ranging from -0.53 to -0.89. Relevant associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU5 and OCU19 containing some relevant candidates (TANC2, ACE, MAP3K3, TEX2, PRKCA, SH3GL2, CNTLN, RPGRIP1L and FTO) related to lipid metabolism, binding, and obesity. These regions explained about 1.2 to 13.9% of the total genomic variance of the traits studied. Our results showed an important maternal genetic effect on IMF and its fatty acid composition in rabbits and identified promising candidate genes associated with these traits.

Relationships of genomic estimated breeding values for age at puberty, birth weight, and growth during development in normal cyclic and acyclic gilts

Managing replacement gilts to reach optimal body weight and growth rate for boar stimulation and first breeding is a key component for sow reproductive longevity and producer profitability. Failure to display pubertal estrus remains a major reason that gilts are culled from the herd. Puberty is metabolically gated so evaluating phenotypic and genetic relationships between birth weight and growth traits with age at puberty and acyclicity can provide valuable insight for efficient gilt development. Data on a litter of origin of the gilt, average daily gain at different stages of development, and age at puberty were available for age-matched cyclic (n = 4,861) and acyclic gilts (prepubertal anestrus, n = 578; behavioral anestrus, n = 428). Genomic estimated breeding values were predicted for each trait using genomic best linear unbiased prediction. Primiparous sows produced more acyclic gilts than multiparous sows (P < 0.05). Accounting for effects of parity and litter size, prepubertal anestrus gilts were heavier at birth and behaviorally anestrus gilts grew faster during the finisher period compared to cyclic gilts (P < 0.05), reflecting possible prenatal programming that negatively affects optimal pubertal development and antagonistic effects between adolescent growth and expression of estrus of gilts from first parity sows. Regression of phenotypic age at puberty with lifetime growth rate (birth to selection) showed a negative linear relationship whereas genomic estimated breeding values showed a negative quadratic relationship indicating that gilts with the least and greatest growth are less optimal as replacements. The slopes of these relationships are small with low negative phenotypic (r = -0.06) and genetic correlations (r = -0.13). The addition of data from acyclic gilts did not substantially change the estimates for genetic relationships between growth and pubertal onset. Although this study identified differences in birth weight and growth rate between cyclic and acyclic gilts the genetic relationships are weak, suggesting that genetic selection for these traits can be achieved separately. Avoiding the smallest and largest gilts in a cohort born to first parity sows could result in gilts with optimal development and reduce the proportion of replacement gilts that are acyclic.

Genome-wide analyses reveal the regions involved in the phenotypic diversity in Sicilian pigs

Nero Siciliano (Sicilian Black, SB) is a local pig breed generally of uniform black color. In addition to this officially recognized breed, there are animals showing morphological characteristics resembling the SB but with gray hair (Sicilian Grey, SG). The SG, compared with the SB, also shows a more compact structure with greater transverse diameters, higher average daily gains and lower thickness of the back fat. In this study, using the Illumina PorcineSNP60 BeadChip, we run genome-wide analyses to identify regions that may explain the phenotypic differences between SB (n = 21) and SG (n = 27) individuals. Combining the results of the two case-control approaches (GWAS and F_ST ), we identified two significant regions, one on SSC5 (95 401 083 bp) and one on SSC15 (55 051 435 bp), which contains several candidate genes related to growth traits in pig. The results of the Bayesian population differentiation approach identified a marker near the MGAT4C, a gene associated with average daily gain in pigs. Finally, scanning the genome for runs of homozygosity islands, we found that the two groups have different runs of homozygosity islands, with several candidate genes involved in coat color (in SG) or related to different pig performance traits (in SB). In summary, the two analyzed groups differed for several phenotypic traits, and genes involved in these traits (growth, meat traits and coat color) were detected. This study provided another contribution to the identification of genomic regions involved in phenotypic variability in local pig populations.

Genome-wide association study for backfat thickness at 100 kg and loin muscle thickness in domestic pigs based on genotyping by sequencing

Both backfat thickness at 100 kg (B100) and loin muscle thickness (LMT) are economically important traits in pigs. In this study, a total of 1,200 pigs (600 Landrace and 600 Yorkshire pigs) were examined with genotyping by sequencing. A total of 345,570 single nucleotide polymorphisms (SNPs) were obtained from 1,200 pigs. Then, a single marker regression test was used to conduct a genome-wide association study for B100 and LMT. A total of 8 and 90 significant SNPs were detected for LMT and B100, respectively. Interestingly, two shared significant loci [located at Sus scrofa chromosome (SSC) 6: 149876694 and SSC12: 46226580] were detected in two breeds for B100. Furthermore, three potential candidate genes were found for LMT and B100. The positional candidate gene FAM3C (SSC18: 25573656, P = 2.48 × 10^-9), which controls the survival, growth, and differentiation of tissues and cells, was found for LMT in Landrace pigs. At SSC9: 6.78-6.82 Mb in Landrace pigs, the positional candidate gene, INPPL1, which has a negative regulatory effect on diet-induced obesity and is involved in the regulation of insulin function, was found for B100. The candidate gene, RAB35, which regulates the adipocyte glucose transporter SLC2A4/GLUT4, was identified at approximately SSC14: 40.09-40.13 Mb in Yorkshire pigs. The results of this GWAS will greatly advance our understanding of the genetic architecture of the LMT and B100 traits. However, these identified loci and genes need to be further verified in more pig populations, and their functions also need to be validated by more biological experiments in pigs.

Using genome wide association studies to identify common QTL regions in three different genetic backgrounds based on Iberian pig breed

One of the major limitation for the application of QTL results in pig breeding and QTN identification has been the limited number of QTL effects validated in different animal material. The aim of the current work was to validate QTL regions through joint and specific genome wide association and haplotype analyses for growth, fatness and premier cut weights in three different genetic backgrounds, backcrosses based on Iberian pigs, which has a major role in the analysis due to its high productive relevance. The results revealed nine common QTL regions, three segregating in all three backcrosses on SSC1, 0–3 Mb, for body weight, on SSC2, 3–9 Mb, for loin bone-in weight, and on SSC7, 3 Mb, for shoulder weight, and six segregating in two of the three backcrosses, on SSC2, SSC4, SSC6 and SSC10 for backfat thickness, shoulder and ham weights. Besides, 18 QTL regions were specifically identified in one of the three backcrosses, five identified only in BC_LD, seven in BC_DU and six in BC_PI. Beyond identifying and validating QTL, candidate genes and gene variants within the most interesting regions have been explored using functional annotation, gene expression data and SNP identification from RNA-Seq data. The results allowed us to propose a promising list of candidate mutations, those identified in PDE10A, DHCR7, MFN2 and CCNY genes located within the common QTL regions and those identified near ssc-mir-103-1 considered PANK3 regulators to be further analysed.

Polymorphisms of the porcine cathepsins, growth hormone-releasing hormone and leptin receptor genes and their association with meat quality traits in Ukrainian Large White breed

Cathepsins, growth hormone-releasing hormone (GHRH) and leptin receptor (LEPR) genes have been receiving increasing attention as potential markers for meat quality and pig performance traits. This study investigated the allele variants in four cathepsin genes (CTSB, CTSK, CTSL, CTSS), GHRH and LEPR in pure-bred Ukrainian Large White pigs and evaluated effects of the allele variants on meat quality characteristics. The study was conducted on 72 pigs. Genotyping was performed using PCR-RFLP technique. Meat quality characteristics analysed were intramuscular fat content, tenderness, total water content, ultimate pH, crude protein and ashes. A medium level of heterozygosity values was established for GHRH and LEPR genes which corresponded to very high levels of informativeness indexes. Cathepsins CTSL, CTSB and CTSK had a low level of heterozygosity, and CTSS did not segregate in this breed. Association studies established that intramuscular fat content and tenderness were affected by the allele variance in GHRH and LEPR but not by CTSB and CTSL genes. The GHRH results could be particularly relevant for the production of lean prime cuts as the A allele is associated with both, a lower meat fat content and better tenderness values, which are two attributes highly regarded by consumers. Results of this study suggest that selective breeding towards GHRH/AA genotype would be particularly useful for improving meat quality characteristics in the production systems involving lean Large White lines, which typically have less than 2 % intramuscular fat content.

A genome-wide association scan in pig identifies novel regions associated with feed efficiency trait

Feed conversion ratio (FCR) is an economically important trait in pigs, and feed accounts for a significant proportion of the costs involved in pig production. In this study we used a high-density SNP chip panel, Porcine SNP60 BeadChip, to identify the association between FCR and SNP markers and to study the genetic architecture of the trait. After quality control, a total of 30,847 SNP that could be mapped to the 18 porcine autosomes (SSC) using the pig genome assembly 10.2 were used in the analyses. Deregressed estimated breeding value was used as the response variable. A total of 3,071 Duroc pigs had both FCR data and genotype data. The linkage disequilibrium (r(2)) between adjacent markers was 0.56. Two association mapping approaches were used: a linear mixed model (LMM) based on single-locus regression analysis and a Bayesian variable selection approach (BVS). A total of 79 significant (P < 0.0001) SNP associations on 6 chromosomes were identified by LMM analyses. Out of these, 10 SNP crossed the genome-wide significance threshold. These 10 SNP were all located on SSC 4 and 14. In the BVS analysis, a total of 44 SNP located on 12 chromosomes had posterior probability more than or equal to 0.05 (i.e., Bayes factor ≥ 10). Thirteen SNP were identified by both LMM and BVS. These 13 SNP were located on 4 chromosomes: SSC 4, 7, 8, and 14. Hypoxia inducible factor 1, alpha subunit inhibitor (HIF1AN) and ladybird homeobox 1 (LBX1) are 2 possible candidate genes affecting FCR on SSC 4 and 14, respectively. The study provides a list of SNP associated with FCR and also offers valuable information on the genetic architecture and candidate genes for this trait.

Genome-wide linkage analysis of QTL for growth and body composition employing the PorcineSNP60 BeadChip

BACKGROUND

The traditional strategy to map QTL is to use linkage analysis employing a limited number of markers. These analyses report wide QTL confidence intervals, making very difficult to identify the gene and polymorphisms underlying the QTL effects. The arrival of genome-wide panels of SNPs makes available thousands of markers increasing the information content and therefore the likelihood of detecting and fine mapping QTL regions. The aims of the current study are to confirm previous QTL regions for growth and body composition traits in different generations of an Iberian x Landrace intercross (IBMAP) and especially identify new ones with narrow confidence intervals by employing the PorcineSNP60 BeadChip in linkage analyses.

RESULTS

Three generations (F3, Backcross 1 and Backcross 2) of the IBMAP and their related animals were genotyped with PorcineSNP60 BeadChip. A total of 8,417 SNPs equidistantly distributed across autosomes were selected after filtering by quality, position and frequency to perform the QTL scan. The joint and separate analyses of the different IBMAP generations allowed confirming QTL regions previously identified in chromosomes 4 and 6 as well as new ones mainly for backfat thickness in chromosomes 4, 5, 11, 14 and 17 and shoulder weight in chromosomes 1, 2, 9 and 13; and many other to the chromosome-wide signification level. In addition, most of the detected QTLs displayed narrow confidence intervals, making easier the selection of positional candidate genes.

CONCLUSIONS

The use of higher density of markers has allowed to confirm results obtained in previous QTL scans carried out with microsatellites. Moreover several new QTL regions have been now identified in regions probably not covered by markers in previous scans, most of these QTLs displayed narrow confidence intervals. Finally, prominent putative biological and positional candidate genes underlying those QTL effects are listed based on recent porcine genome annotation.

Association of A-FABP gene polymorphism in intron 1 with meat quality traits in Junmu No. 1 white swine

This study was designed to investigate the single nucleotide polymorphism (SNP) in intron 1 of the gene A-FABP in 127 Junmu No. 1 white swine using PCR-SSCP. The association between the polymorphism and meat quality traits was also studied. The cloning and sequencing results indicated that the polymorphism of intron 1 was due to a point mutation in position 3481bp of A-FABP, giving 3 genotypes (CC, CD and DD). Association analysis indicated that the polymorphism had a significant effect on marbling (P<0.05). Genotype DD had higher marbling than CD and CC, but the difference between CD and CC was no significant. Polymorphism had a highly significant effect on intramuscular fat (IMF) content (P<0.01). DD was higher than CD, which was higher than CC. No significant conclusions can be drawn regarding other traits. Immunoblot analysis of A-FABP levels was carried out on 3 different genotype individuals. Expression was markedly reduced in DD compared with genotype CC. Thus A-FABP may be a candidate gene or a quantitative trait locus-linked gene associated with meat quality traits.

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.

The porcine TBC1D1 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.

Identification of three novel SNPs and association with carcass traits in porcine TNNI1 and TNNI2

In this study, two novel SNPs (EU743939:g.5174T>C in intron 4 and EU743939:g.8350C>A in intron 7) in TNNI1 and one SNP (EU696779:g.1167C>T in intron 3) in TNNI2 were identified by PCR-RFLP (PCR restriction fragment length polymorphism) using XbaI, MspI and SmaI restriction enzyme, respectively. The allele frequencies of three novel SNPs were determined in the genetically diverse pig breeds including ten Chinese indigenous pigs and three Western commercial pig breeds. Association analysis of the SNPs with the carcass traits were conducted in a Large White × Meishan F(2) pig population. The linkage of two SNPs (g.5174T>C and g.8350C>A) in TNNI1 gene had significant effect on fat percentage. Besides these, the g.5174T>C polymorphism was also significantly associated with skin percentage (P < 0.05), shoulder fat thickness (P < 0.05) and backfat thickness between sixth and seventh ribs (P < 0.05). The significant effects of g.1167C>T polymorphism in TNNI2 gene on fat percentage (P < 0.01), lean meat percentage (P < 0.05), lion eye area (P < 0.05), thorax-waist backfat thickness (P < 0.01) and average backfat thickness (P < 0.05) were also found.

Confirmed association between a single nucleotide polymorphism in the FTO 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).