Genetic parameters and direct, maternal and heterosis effects on litter size in a diallel cross among three commercial varieties of Iberian pig.

The Iberian pig is one of the pig breeds that has the highest meat quality. Traditionally, producers have bred one of the available varieties, exclusively, and have not used crosses between them, which has contrasted sharply with other populations of commercial pigs for which crossbreeding has been a standard procedure. The objective of this study was to perform an experiment under full diallel design among three contemporary commercial varieties of Iberian pig and estimate the additive genetic variation and the crossbreeding effects (direct, maternal and heterosis) for prolificacy. The data set comprised 18 193 records for total number born and number born alive from 3800 sows of three varieties of the Iberian breed (Retinto, Torbiscal and Entrepelado) and their reciprocal crosses (Retinto × Torbiscal, Torbiscal × Retinto, Retinto × Entrepelado, Entrepelado × Retinto, Torbiscal × Entrepelado and Entrepelado × Torbiscal), and a pedigree of 4609 individuals. The analysis was based on a multiple population repeatability model, and we developed a model comparison test that indicated the presence of direct line, maternal and heterosis effects. The results indicated the superiorities of the direct line effect of the Retinto and the maternal effect of the Entrepelado populations. All of the potential crosses produced significant heterosis, and additive genetic variation was higher in the Entrepelado than it was in the other two populations. The recommended cross for the highest yield in prolificacy is a Retinto father and an Entrepelado mother to generate a hybrid commercial sow.

Five genomic regions have a major impact on fat composition in Iberian pigs.

The adipogenic nature of the Iberian pig defines many quality attributes of its fresh meat and dry-cured products. The distinct varieties of Iberian pig exhibit great variability in the genetic parameters for fat deposition and composition in muscle. The aim of this work is to identify common and distinct genomic regions related to fatty acid composition in Retinto, Torbiscal, and Entrepelado Iberian varieties and their reciprocal crosses through a diallelic experiment. In this study, we performed GWAS using a high density SNP array on 382 pigs with the multimarker regression Bayes B method implemented in GenSel. A number of genomic regions showed strong associations with the percentage of saturated and unsaturated fatty acid in intramuscular fat. In particular, five regions with Bayes Factor >100 (SSC2 and SSC7) or >50 (SSC2 and SSC12) explained an important fraction of the genetic variance for miristic, palmitoleic, monounsaturated (>14%), oleic (>10%) and polyunsaturated (>5%) fatty acids. Six genes (RXRB, PSMB8, CHGA, ACACA, PLIN4, PLIN5) located in these regions have been investigated in relation to intramuscular composition variability in Iberian pigs, with two SNPs at the RXRB gene giving the most consistent results on oleic and monounsaturated fatty acid content.

Genetic parameters and crossbreeding effects of fat deposition and fatty acid profiles in Iberian pig lines.

The aim of this study was to estimate the genetic and environmental parameters and crossbreeding effects on fatty acid and fat traits in the Iberian pig. Our final goal is to explore target selection traits and define crossbreeding strategies. The phenotypes were obtained under intensive management from 470 animals in a diallelic experiment involving Retinto, Torbiscal, and Entrepelado lines. The data set was composed of backfat thickness at the fourth rib (BFT), intramuscular fat (IMF) in the longissimus thoracis (LT), and the fatty acid profile for IMF and subcutaneous fat (SCF) traits. Data were analyzed through a Bayesian bivariate animal model by using a reparameterization of Dickerson's model. The results obtained showed an important genetic determinism for all traits analyzed with heritability ranging from 0.09 to 0.67. The common environment litter effect also had an important effect on IMF (0.34) and its fatty acid composition (0.06-0.53) at slaughter. The additive genetic correlation between BFT and IMF (additive genetic correlation [] = 0.31) suggested that it would be possible to improve lean growth independent of the IMF with an appropriate selection index. Furthermore, the high additive genetic correlation ( = 0.68) found between MUFA tissues would seem to indicate that either the LT or SCF could be used as the reference tissue for MUFA selection. The relevance of the crossbreeding parameters varied according to the traits analyzed. Backfat thickness at the fourth rib and the fatty acid profile of the IMF showed relevant differences between crosses, mostly due to line additive genetic effects associated with the Retinto line. On the contrary, those for IMF crosses were probably mainly attributable to heterosis effects. Particularly, heterosis effects were relevant for the Retinto and Entrepelado crosses (approximately 16% of the trait), which could be valuable for a crossbreeding system involving these lines.

Endometrial gene expression profile of pregnant sows with extreme phenotypes for reproductive efficiency.

Prolificacy can directly impact porcine profitability, but large genetic variation and low heritability have been found regarding litter size among porcine breeds. To identify key differences in gene expression associated to swine reproductive efficiency, we performed a transcriptome analysis of sows' endometrium from an Iberian x Meishan F2 population at day 30-32 of gestation, classified according to their estimated breeding value (EBV) as high (H, EBV > 0) and low (L, EBV < 0) prolificacy phenotypes. For each sample, mRNA and small RNA libraries were RNA-sequenced, identifying 141 genes and 10 miRNAs differentially expressed between H and L groups. We selected four miRNAs based on their role in reproduction, and five genes displaying the highest differences and a positive mapping into known reproductive QTLs for RT-qPCR validation on the whole extreme population. Significant differences were validated for genes: PTGS2 (p = 0.03; H/L ratio = 3.50), PTHLH (p = 0.03; H/L ratio = 3.69), MMP8 (p = 0.01; H/L ratio =4.41) and SCNN1G (p = 0.04; H/L ratio = 3.42). Although selected miRNAs showed similar expression levels between H and L groups, significant correlation was found between the expression level of ssc-miR-133a (p < 0.01) and ssc-miR-92a (p < 0.01) and validated genes. These results provide a better understanding of the genetic architecture of prolificacy-related traits and embryo implantation failure in pigs.

Associations between pig adiponectin (ADIPOQ) genotype and serum lipid levels are modulated by age-specific modifiers.

The adiponectin (ADIPOQ) locus is a positional and functional candidate gene for 2 porcine chromosome 13 (SSC13) QTL influencing cholesterol (CHOL) and low-density lipoprotein (LDL) concentrations in 190-d-old pigs. By sequencing 2.37 kb of the pig ADIPOQ cDNA, we have identified 1 c.*1512G>T 3' untranslated region polymorphism that has been genotyped in a Duroc pig commercial population with records for serum lipid levels at 45 and 190 d of age. Statistical analysis of the data have revealed significant associations between the ADIPOQ genotype and CHOL (P=0.0040) and LDL (P=0.0011) concentrations at 190 d but not at 45 d. In family 3, most of the SSC13 QTL effects on LDL levels at 190 d were explained by the ADIPOQ genotype. We also found an association with triglyceride levels at 45 d (P=0.0060) but not at 190 d. Measurement of allelic mRNA imbalance demonstrated that the G and T alleles are expressed at very similar levels in muscle and fat tissues, indicating that the c.*1512G>T polymorphism does not affect transcript abundance. As a whole, results obtained in the current work as well as previous data gathered in humans and pigs provide evidence that the magnitude of associations between blood lipid phenotypes and candidate loci genotypes may vary depending on the age of the individual, therefore suggesting the existence of dynamic genotype×environment interactions changing on a temporal scale.

Crossbreeding effects on pig growth and carcass traits from two Iberian strains.

An experiment of a 2×2 full diallelic cross between two contemporary Iberian pig strains (Retinto: RR, and Torbiscal: TT) was conducted to estimate the crossbreeding effects for growth and carcass traits. Phenotypic records were obtained under intensive management and consisted of two different data sets. The first set comprised growth traits until weaning and was collected at two different farms (6236 and 1208 records, respectively). Specific data included individual piglet weight at birth and at weaning at 28 days and average daily gain from birth to weaning at 28 days (ADG28) for both RR and TT and their reciprocal crosses. The second set comprised growth data from birth to slaughter (~340 days and ~160 kg) and carcass traits from 349 individuals (randomly) sampled at weaning from the first dataset. Data were analyzed through a Bayesian analysis by using a reparameterization of Dickerson's model that allowed estimation of the posterior distributions of the following crossbreeding effects: average maternal breed effect (gM), average paternal breed effect (gP) and individual heterosis (hI). Results showed that the relative magnitude of crossbreeding effects depends on the trait analyzed. Crosses where Torbiscal strain was used as mother (RT and TT) achieved the greatest performance for all growth traits at weaning, leading to remarkable gM effects. The most outstanding example is the case of ADG28 where the probability of relevance was one. In contrast, TR cross showed the greatest differences from RR cross for all growth at slaughter and carcass traits. These differences were mainly due to hI and gP crossbreeding parameters. In particular, the posterior mean of hI was more noticeable for live weight at slaughter, average daily gain at slaughter and carcass length, while gP was more relevant for hams (kg) and loins (kg) representing from 3% to10% of average performance of traits. Hence, growth traits at weaning did not reveal any notable advantage of the crossbreeding scheme because of the superiority of the Torbiscal strain with respect to its mothering ability and the small hI. However, results from growth and carcass traits at slaughter would support the implementation of a TR crossbred system. It would allow exploitation of both the gP of the Torbiscal strain and the hI between these two Iberian pig strains. Additionally, gP estimates and phenotypic differences between reciprocal crosses might suggest signs of the presence of paternal genetic imprinting in primal cuts traits.

Transcriptional analysis of intramuscular fatty acid composition in the longissimus thoracis muscle of Iberian × Landrace back-crossed pigs.

This study aimed at identifying differential gene expression conditional on the fatty acid profile of the longissimus thoracis (Lt) muscle, a prime cut of economic relevance for fresh and cured pork production. A population of 110 Iberian (25%) × Landrace (75%) back-crossed pigs was used, because these two breeds exhibit extreme profiles of intramuscular saturated fatty acid, monounsaturated fatty acid (MUFA) and polyunsaturated fatty acid (PUFA) contents. Total RNA from Lt muscle was individually hybridized to GeneChip Porcine Genome arrays (Affymetrix). A principal component analysis was performed with data from the 110 animals to select 40 extreme animals based on the total fatty acid profile and the MUFA composition (MAP). Comparison of global transcription levels between extreme fatty acid profile pigs (n = 40) resulted in 219 differentially expressed probes (false discovery rate <0.10). Gene ontology, pathway and network analysis indicated that animals with higher percentages of PUFA exhibit a shift toward a more oxidative muscular metabolism state, with a raise in mitochondria function (PPARGC1A, ATF2), fatty acid uptake and oxidation (FABP5, MGLL). On the other hand, 87 probes were differentially expressed between MUFA composition groups (n = 40; false discovery rate <0.10). In particular, muscles rich in n-7 MUFA expressed higher levels of genes involved in lipid metabolism (GLUL, CRAT, PLA2G15) and lower levels of fatty acid elongation genes (ELOVL5). Moreover, the chromosomal position of FABP5, PAQR3, MGLL, PPARGC1A, GLUL and ELOVL5 co-localized with very relevant QTL for fat deposition and composition described in the same resource population. This study represents a complementary approach to identifying genes underlying these QTL effects.

Variability-specific differential gene expression across reproductive stages in sows.

Differential gene expression analyses typically focus on departures across mathematical expectations (i.e. mean) from two or more groups of microarrays, without considering alternative patterns of departure. Nevertheless, recent studies in humans and great apes have suggested that differential gene expression could also be characterized in terms of heterogeneous dispersion patterns. This must be viewed as a very interesting genetic phenomenon clearly linked to the regulation mechanisms of gene transcription. Unfortunately, we completely lack information about the incidence and relevance of dispersion-specific differential gene expression in livestock species, although a specific Bayes factor (BF) for testing this kind of differential gene expression (i.e. within-probe heteroskedasticity) has been recently developed. Within this context, our main objective was to characterize the incidence of dispersion-specific differential gene expression in pigs and, if possible, providing the first evidence of this phenomenon in a livestock species. We evaluated dispersion-specific differential gene expression on ovary, uterus and hypophysis samples from 22 F2 Iberian × Meishan sows, where a total of 15,252 probes were interrogated. For each tissue, heteroskedasticity of probe-specific residual variances was evaluated by three pairwise comparisons involving three physiological stages, that is, heat, 15 days of pregnancy and 45 days of pregnancy. Between 2.9% and 37.4% of the analyzed probes provided statistical evidence of within-tissue across-physiological stages dispersion-specific differential gene expression (BF >1), and between 0.1% and 3.0% of them reported decisive evidence (BF >100). It is important to highlight that <8% of the heteroskedastic probes were also linked to differential gene expression in terms of departures among the probe-specific mathematical expectation of each physiological stage. This discarded the disturbance of scale effects in a high percentage of probes and suggested that probe-specific heteroskedasticity must be viewed as an independent phenomenon within the context of differential gene expression. As a whole, our results report a remarkable incidence of dispersion-specific differential gene expression across the whole genome of the pig, establishing a very interesting starting point for further studies focused on deciphering the genetic mechanisms underlying heteroskedasticity.

Genome-wide association study for intramuscular fatty acid composition in an Iberian × Landrace cross.

The lipid content and fatty acid (FA) profile have an important impact in human health as well as in the technological transformation and nutritional and organoleptic quality of meat. A genome-wide association study (GWAS) on 144 backcross pigs (25% Iberian × 75% Landrace) was performed for 32 traits associated with intramuscular FA composition and indices of FA metabolism. The GWAS was carried out using Qxpak 5.0 and the genotyping information obtained from the Porcine SNP60K BeadChip (Illumina Inc., San Diego, CA). Signals of significant association considering a false- discovery rate (q-value < 0.05) were observed in 15 of the 32 analyzed traits, and a total of 813 trait-associated SNP (TAS), distributed in 43 chromosomal intervals on almost all autosomes, were annotated. According to the clustering analysis based on functional classification, several of the annotated genes are related to FA composition and lipid metabolism. Some interesting positional concordances among TAS and previously reported QTL for FA compositions and/or other lipid traits were also found. These common genomic regions for different traits suggest pleiotropic effects for FA composition and were found primarily on SSC4, SSC8, and SSC16. These results contribute to our understanding of the complex genetic basis of FA composition and FA metabolism.

Quantitative trait loci analysis of a Duroc commercial population highlights differences in the genetic determination of meat quality traits at two different muscles.

We performed a whole-genome scan with 110 informative microsatellites in a commercial Duroc population for which growth, fatness, carcass and meat quality phenotypes were available. Importantly, meat quality traits were recorded in two different muscles, that is, gluteus medius (GM) and longissimus thoracis et lumborum (LTL), to find out whether these traits are determined by muscle-specific genetic factors. At the whole-population level, three genome-wide QTL were identified for carcass weight (SSC7, 60 cM), meat redness (SSC13, 84 cM) and yellowness (SSC15, 108 cM). Within-family analyses allowed us to detect genome-wide significant QTL for muscle loin depth between the 3rd and 4th ribs (SSC15, 54 cM), backfat thickness (BFT) in vivo (SSC10, 58 cM), ham weight (SSC9, 69 cM), carcass weight (SSC7, 60 cM; SSC9, 68 cM), BFT on the last rib (SSC11, 48 cM) and GM redness (SSC8, 85 cM; SSC13, 84 cM). Interestingly, there was low positional concordance between meat quality QTL maps obtained for GM and LTL. As a matter of fact, the three genome-wide significant QTL for colour traits (SSC8, SSC13 and SSC15) that we detected in our study were all GM specific. This result suggests that QTL effects might be modulated to a certain extent by genetic and environmental factors linked to muscle function and anatomical location.

Sequencing and gene expression of the porcine ITIH SSC13 cluster and its effect on litter size in an Iberian × Meishan F2 population.

The aim of the present study was to identify polymorphisms and to analyze endometrial gene expression of the porcine SSC13 ITIH cluster that could explain differences in prolificacy of 255 F(2) sows derived from an Iberian (Ib)×Meishan (Me) intercross in which QTL for the number of piglets born alive (NBA) and total number of piglets born (TNB) were previously detected on this chromosome. Sequencing of ITIH-1, -3, and -4 mRNAs was done and several polymorphisms segregating within the Ib×Me population were found in all three genes. Significant associations with NBA were found for two SNPs from ITIH-1, four from ITIH-3, and four SNPs from ITIH-4 (p<0.05). Haplotypes for the significant SNPs were calculated by segregation analysis and a marker assisted association test indicated that the alleles coming from the Meishan breed had a favorable effect on NBA for all three genes. Interestingly, some of the significant SNPs were located within the von Willebrand domain of the ITIH proteins, the binding site of molecules essential for the synthesis of the extracellular matrix during cumulus expansion. Gene expression analyses also revealed differences in the expression level of the ITIH-3 gene regarding the prolificacy performance (high or low) and the uterus sample (apical or basal).

Quantitative trait loci for fatness at growing and reproductive stages in Iberian × Meishan F(2) sows.

A considerable number of fatness QTL have been identified in growing pigs, but there is a lack of knowledge about the genetic architecture of this trait in gilts and sows. We have performed a genome scan, in 255 Iberian × Meishan F(2) sows, for backfat thickness (BF) at 150 (BF(150) ) and 210 (BF(210)) days of age, 30 days after conception (BF(30)) and 7-10 days before farrowing (BF(bf)). We have found one BF150 QTL in SSC6 (120 cM) that was highly significant (P < 0.001) at the chromosome-wide level and suggestive at the genome-wide level (P < 0.1). Ten additional chromosome-wide significant QTL were found for sow BF(150) (SSC1, SSC13), BF(210) (SSC6, SSC8, SSC15), BF(30) (SSC5, SSC6) and BF(bf) (SSC1, SSC6, SSC13). The location of several of the BF QTL varied depending on the growing and reproductive status of the sow, suggesting that part of these genetic effects may have a temporal pattern of phenotypic expression.

Building phenotype networks to improve QTL detection: a comparative analysis of fatty acid and fat traits in pigs.

Models in QTL mapping can be improved by considering all potential variables, i.e. we can use remaining traits other than the trait under study as potential predictors. QTL mapping is often conducted by correcting for a few fixed effects or covariates (e.g. sex, age), although many traits with potential causal relationships between them are recorded. In this work, we evaluate by simulation several procedures to identify optimum models in QTL scans: forward selection, undirected dependency graph and QTL-directed dependency graph (QDG). The latter, QDG, performed better in terms of power and false discovery rate and was applied to fatty acid (FA) composition and fat deposition traits in two pig F2 crosses from China and Spain. Compared with the typical QTL mapping, QDG approach revealed several new QTL. To the contrary, several FA QTL on chromosome 4 (e.g. Palmitic, C16:0; Stearic, C18:0) detected by typical mapping vanished after adjusting for phenotypic covariates in QDG mapping. This suggests that the QTL detected in typical mapping could be indirect. When a QTL is supported by both approaches, there is an increased confidence that the QTL have a primary effect on the corresponding trait. An example is a QTL for C16:1 on chromosome 8. In conclusion, mapping QTL based on causal phenotypic networks can increase power and help to make more biologically sound hypothesis on the genetic architecture of complex traits.

Expression study on the porcine PTHLH gene and its relationship with sow teat number.

Teat number is an important trait in sows that should accompany the increase in litter size that has been achieved in the last decades through selection. We have previously identified a genome-wide significant QTL for teat number in porcine chromosome SSC5 by means of an experimental Meishan by Iberian F(2) intercross population. In the present report, we have studied the porcine parathyroid hormone-like hormone (PTHLH) gene, which maps to SSC5, as a candidate gene for this trait, as PTHLH is involved in nipple formation during embryogenesis and nipple development during pregnancy and lactation. We have found that porcine PTHLH gene is transcribed into three mRNA species differing in the 5'UTR region. Two of these variants are reported in pigs here for the first time: one was similar to variant 1 described in humans while the other, which was generated by the retention of two small introns, has not been identified before in any other species. In addition, mRNA expression profile for two of the mRNA variants was assessed in 19 pig tissues. Porcine PTHLH showed a widespread expression as it was present in all tested tissues and relative expression of each variant was tissue dependent. Finally, we have performed an association study between a non-synonymous mutation in the coding region of this gene and sow teat number. The PTHLH polymorphism was segregating in our Meishan by Iberian F(2) population at intermediate allelic frequencies. We compared here six different statistical models to choose the one with a better fit and a lower degree of complexity. However, despite the potential negative effect of the PTHLH mutation in the signal peptide of this protein, we did not detect any association between the PTHLH genotype and the sow teat number phenotype, concluding that the causal mutation of the observed QTL is very likely not related to this gene.

Genetic evaluation combining purebred and crossbred data in a pig breeding scheme.

Genetic evaluations using purebred data alone and combined purebred and crossbred information were performed for lean meat percentage in a pig breeding scheme. One purebred (PB) model and 2 crossbred models (CCPS1 and CCPS2) were used in the analyses. Data were obtained from the Selección Batallé S.A. Company (Riudarenes, Spain) and spanned a period of 4 yr (2006 to 2009). The data corresponded to 3 nuclei of purebred populations, Landrace (LD), Duroc (DU), and Pietrain (PI); 1 multiplying farm with animals from a 2-way cross (TB1; DU × LD); and commercial farms with animals from a 3-way cross (TB2; TB1 × PI). Genetic parameters were similar across the models, with the exception of purebred PI. The DU and LD purebreds presented large heritabilities (0.5 to 0.6) for lean meat percentage, whereas the PI purebred showed a lower heritability (approximately 0.1) for the PB model and moderate heritability for the CCPS1 and CCPS2 models (0.2 to 0.3). The mean reliability of the predicted purebred breeding values was clearly increased when the CCPS1 and CCPS2 models were used. Moreover, a reranking of the animals with important changes in the selection decisions was observed in the PI purebred. In a simulation study, the CCPS1 model achieved a greater response to selection than the PB model for the PI purebred. On another hand, between the CCPS1 and CCPS2 models, CCPS1 was slightly superior in terms of predictive ability, exhibiting a greater robustness. These results illustrate the usefulness of using crossbred models to evaluate lean meat percentage in this pig breeding scheme.

Porcine intramuscular fat content and composition are regulated by quantitative trait loci with muscle-specific effects.

Intramuscular fat (IMF) storage is a biological process with a strong impact on nutritional and technological properties of meat and also with relevant consequences on human health. The genetic architecture of IMF content and composition phenotypes has been thoroughly studied in pigs through the identification of QTL and the estimation of genetic parameters. A question that has not been elucidated yet is if the genetic determinants of IMF-related phenotypes are muscle specific or, conversely, have broad effects on the whole skeletal muscle compartment. We have addressed this question by generating lipid QTL maps for 2 muscles with a high commercial value, gluteus medius (GM) and longissimus thoracis et lumborum (LTL), in a Duroc commercial population (n = 350). Our data support a lack of concordance between the GM and LTL QTL maps, suggesting that the effects of polymorphisms influencing IMF, cholesterol, and fatty acid contents are modulated to some extent by complex spatial factors related to muscle location, metabolism, and function. These results have important implications on the implementation of genomic selection schemes aimed to improve the lipid profile of swine meat.

An application of change-point recursive models to the relationship between litter size and number of stillborns in pigs.

We developed and implemented change-point recursive models and compared them with a linear recursive model and a standard mixed model (SMM), in the scope of the relationship between litter size (LS) and number of stillborns (NSB) in pigs. The proposed approach allows us to estimate the point of change in multiple-segment modeling of a nonlinear relationship between phenotypes. We applied the procedure to a data set provided by a commercial Large White selection nucleus. The data file consisted of LS and NSB records of 4,462 parities. The results of the analysis clearly identified the location of the change points between different structural regression coefficients. The magnitude of these coefficients increased with LS, indicating an increasing incidence of LS on the NSB ratio. However, posterior distributions of correlations were similar across subpopulations (defined by the change points on LS), except for those between residuals. The heritability estimates of NSB did not present differences between recursive models. Nevertheless, these heritabilities were greater than those obtained for SMM (0.05) with a posterior probability of 85%. These results suggest a nonlinear relationship between LS and NSB, which supports the adequacy of a change-point recursive model for its analysis. Furthermore, the results from model comparisons support the use of recursive models. However, the adequacy of the different recursive models depended on the criteria used: the linear recursive model was preferred on account of its smallest deviance value, whereas nonlinear recursive models provided a better fit and predictive ability based on the cross-validation approach.

Bayes factor analyses of heritability for serum and muscle lipid traits in Duroc pigs.

Concern about pork quality has increased during last decades. Given the influence of fat content and composition on sensorial, nutritional, and technological variables of pork meat, an accurate knowledge about genetic control of pig lipid metabolism is required. This study focused on providing a broad characterization for serum and meat lipid trait heritability estimates in pigs. Analyses were performed on a population of 370 Duroc barrows and measured the additive polygenic background for the serum concentrations of cholesterol, triglyceride, and low- and high-density lipoproteins at 45 and 190 d of age (at slaughter), as well as intramuscular fat, cholesterol content, and C:12 to C:22 fatty acid content in longissimus thoracis et lumborum and gluteus medius muscles at slaughter. These traits were analyzed under Bayesian univariate animal linear models, and the statistical relevance of heritability estimates was evaluated through Bayes factor (BF); the model with polygenic additive effects was favored when BF >1. All serum lipid traits showed relevant genetic determinism, but the BF reached greater values at 190 d of age. Serum lipid traits displayed moderate modal estimates for heritability that ranged from 0.18 to 0.30. On the other hand, the genetic determinism for meat quality traits showed a heterogeneous behavior with large and less-than-1 BF. In general, longissimus thoracis et lumborum and gluteus medius muscles showed a similar pattern, with strong evidence of polygenic additive effects for intramuscular fat and palmitic, stearic, and cis-vaccenic fatty acids content, whereas oleic and muscle cholesterol content showed moderate to weak BF with moderate heritabilities. Similarly, results regarding linoleic, arachidonic, n-3, and n-6 fatty acids suggested a moderate genetic determinism, but only in gluteus medius muscle. For the remaining traits (myristic and palmitoleic fatty acids in both muscles, along with linoleic, arachidonic, n-3, and n-6 fatty acids in the longissimus thoracis et lumborum muscle), no statistical evidence for genetic control was observed in this study. As a whole, these results confirm the complexity of lipid metabolism in pigs.

Analysis of candidate genes underlying two epistatic quantitative trait loci on SSC12 affecting litter size in pig.

The previous results from a genome scan for total number of piglets born and number of piglets born alive in a F(2) Iberian by Meishan intercross showed several single and epistatic QTL. One of the most interesting results was obtained for SSC12, where two QTL affecting both traits showed epistatic interaction. In this study, we proposed two genes (SLC9A3R1 and NOS2) as biological and potentially positional candidates underlying these QTL. Both cDNAs were characterized and 23 polymorphisms were detected. A chromosome scan was conducted with 12 markers, plus one SNP in SLC9A3R1 and one in NOS2, covering 110 cM of SSC12. The epistatic QTL (QTL1 at 15 cM and QTL2 at 97 cM) were confirmed, and SLC9A3R1 and NOS2 were mapped around the QTL1 and QTL2 regions respectively. Several SNPs in both genes were tested with standard animal model and marker assisted association tests. The most significant results were obtained with the NOS2 haplotype defined by one missense SNP c.2192C > T (Val to Ala) and a 15 bp duplication at the 3'UTR. This duplication seems to include AU-rich elements, and could be a target site for miRNA, therefore there are statistical and biological indications to consider this haplotype as the potential causal mutation underlying QTL2. SLC9A3R1 results were not conclusive. Although the interaction between the SNPs was not significant, we cannot reject the possibility of interaction of the NOS2 haplotype with other polymorphisms closely linked to the SL9A3R1 SNPs analysed.

Functional and association studies on the pig HMGCR gene, a cholesterol-synthesis limiting enzyme.

The 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) is the rate-limiting enzyme in the biosynthesis of cholesterol. We have studied the role of the HMGCR gene in pig lipid metabolism by means of expression and structural analysis. We describe here the complete coding region of this gene in pigs and report two synonymous single nucleotide polymorphisms in the coding region. We have, additionally, studied the association of one of these polymorphisms (HMGCR:c.807A>C) with several lipid deposition- and cholesterol-related traits in a half-sib population generated from a commercial Duroc line, showing in some families a positive relationship of HMGCR:c.807A allele with serum low-density lipoprotein (LDL)-bound cholesterol and triglyceride levels, and also with intramuscular fat (IMF) content of gluteus medius muscle. We have also assessed the expression levels in muscle and in liver from 68 Duroc individuals corresponding to the most extreme animals for the analysed traits. Liver HMGCR expression correlated negatively with the serum high-density lipoprotein (HDL) levels, carcass lean percentage and stearic acid content, while muscle expression correlated also negatively with the carcass lean percentage, stearic and linoleic acids content, but showed a positive correlation with the serum lipid cholesterol (HDL, LDL and total cholesterol), IMF and muscle oleic and palmitic fatty acid content. With this information, we have performed an association analysis of expression data with lipid metabolism phenotypic levels and the HMGCR genotype. The results indicate that HMGCR expression levels in muscle are different in the two groups of pigs with extreme values for fat deposition and total cholesterol levels, and also between animals with the different HMGCR genotypes.