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.

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.

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.

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).

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.

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.

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.

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.

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.

Bayesian threshold analysis of direct and maternal genetic parameters for piglet mortality at farrowing in Large White, Landrace, and Pietrain populations

A Bayesian threshold model was fitted to analyze the genetic parameters for farrowing mortality at the piglet level in Large White, Landrace, and Pietrain populations. Field data were collected between 1999 and 2006. They were provided by 3 pig selection nucleus farms of a commercial breeding company registered in the Spanish Pig Data Bank (BDporc). Analyses were performed on 3 data sets of Large White (60,535 piglets born from 4,551 litters), Landrace (57,987 piglets from 5,008 litters), and Pietrain (42,707 piglets from 4,328 litters) populations. In the analysis, farrowing mortality was considered as a binary trait at the piglet level and scored as 1 (alive piglet) or 0 (dead piglet) at farrowing or within the first 12 h of life. Each breed was analyzed separately, and operational models included systematic effects (year-season, sex, litter size, and order of parity), direct and maternal additive genetic effects, and common litter effects. Analyses were performed by Bayesian methods using Gibbs sampling. The posterior means of direct heritability were 0.02, 0.06, and 0.10, and the posterior means of maternal heritability were 0.05, 0.13, and 0.06 for Large White, Landrace, and Pietrain populations, respectively. The posterior means of genetic correlation between the direct and maternal genetic effects for Landrace and Pietrain populations were -0.56 and -0.53, and the highest posterior intervals at 95% did not include zero. In contrast, the posterior mean of the genetic correlation between direct and maternal effects was 0.15 in the Large White population, with the null correlation included in the highest posterior interval at 95%. These results suggest that the genetic model of evaluation for the Landrace and Pietrain populations should include direct and maternal genetic effects, whereas farrowing mortality could be considered as a sow trait in the Large White population.

A quantitative trait locus genome scan for porcine muscle fiber traits reveals overdominance and epistasis

Muscle histochemical characteristics are decisive determinants of meat quality. The relative percentage and diameters of the different muscular fiber types influence crucial aspects of meat such as color, tenderness, and ultimate pH. Despite its relevance, however, the information on muscle fiber genetic architecture is scant, because histochemical muscle characterization is a laborious task. Here we report a complete QTL scan of muscle fiber traits in 160 animals from a F(2) cross between Iberian and Landrace pigs using 139 markers. We identified 20 genome regions distributed along 15 porcine chromosomes (SSC1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, and X) with direct and(or) epistatic effects. Epistasis was frequent and some interactions were highly significant. Chromosomes 10 and 11 seemed to behave as hubs; they harbored 2 individual QTL, but also 6 epistatic regions. Numerous individual QTL effects had cryptic alleles, with opposite effects to phenotypic pure breed differences. Many of the QTL identified here coincided with previous reports for these traits in the literature, and there was overlapping with potential candidate genes and previously reported meat quality QTL.

Polymorphism of the pig 2,4-dienoyl CoA reductase 1 gene (DECR1) and its association with carcass and meat quality traits

We characterized the nearly complete coding sequence of the pig 2,4-dienoyl CoA reductase 1 (DECR1) gene, which encodes an enzyme involved in the beta-oxidation of polyunsaturated fatty enoyl-CoA esters and maps on a linoleic QTL located on Chromosome 4. Sequencing of a 937-bp fragment encompassing exons 2 and 10 revealed the existence of two missense SNP at exon 2 (C181 --> G181) and exon 5 (C458 -->G458). These two SNP are associated with Val (C) --> Leu (G) and Ser (C) --> Thr (G) conservative AA replacements at positions 61 and 153 of the DECR1 protein, respectively. Moreover, DECR1 genotyping in a representative sample of 184 pigs from the Large White, Pietrain, Iberian, Duroc, and Landrace breeds demonstrated the existence of disequilibrium linkage between these two SNP (Haplotype 1: C181C458; Haplotype 2: G181G458). An association analysis between DECR1 genotype and growth, carcass, and meat quality traits in a highly selected Landrace population (n = 470) revealed differences among genotypes for isocitrate dehydrogenase activity (highest posterior density [HPD] of 90%), longissimus thoracis pH (HPD of 95%), lightness (HPD of 90 to 95%), and redness (HPD of 95%). Because these associations were not consistently found in the three available genotype comparisons, we believe that exon 2 and 5 polymorphisms at the DECR1 gene might be in linkage disequilibrium with the true causal mutation influencing isocitrate dehydrogenase activity and muscle color and pH.

Adipocyte fatty-acid binding protein is closely associated to the porcine FAT1 locus on chromosome 41

We identified 22 polymorphisms in the adipocyte fatty-acid binding protein (FABP4) gene, a strong positional candidate gene for the FAT1 locus in porcine chromosome 4. The most informative polymorphism, an insertion/deletion in intron 1, together with a single nucleotide polymorphism in intron 3, was genotyped in a cross between Iberian and Landrace pigs. After performing QTL, single marker, and haplotype analyses, we showed that there were at least 2 quantitative trait genes in the FAT1 region and that the FABP4 polymorphism was tightly associated to fatness. A comparison of allelic frequencies in a panel of pig breeds suggested that the Del2634C polymorphism was under indirect selection. We also showed that FABP4 is tightly associated to fatness but not growth. Furthermore, a haplotype analysis suggests that there is genetic heterogeneity at the FAT1 locus within the Landrace breed.

Comparison of genetic models for analysing the effects of a PvuII polymorphism in the oestrogen receptor 1 (ESR1) gene on prolificacy in an Iberian�נMeishan pig population

The effect of a previously reported PvuII polymorphism in oestrogen receptor 1 (ESR1) was analysed in an F(2) population of Iberian x Meishan pigs. We tested three hypotheses: (1) that a causal mutation was fixed in the parental populations, (2) that a causal mutation existed that was in complete linkage disequilibrium with the alleles of the PvuII polymorphism and (3) that a causal mutation existed in linkage disequilibrium within each parental population. The third model was the most plausible based on the available data. ESR1 alleles displayed different patterns of linkage disequilibrium with the causal mutation in each of the parental populations and the PvuII polymorphism was clearly not the causal mutation. As a consequence, the use of the ESR1 mutation for selection must be evaluated for a particular pig population before it is applied.

Malic enzyme 1 genotype is associated with backfat thickness and meat quality traits in pigs

Malic enzyme 1 (ME1) is a part of the tricarboxylate shuttle that provides NADPH and acetyl-CoA required in fatty acid biosynthesis. The pig ME1 locus maps on the proximal end of chromosome 1, where a quantitative trait loci (QTL) affecting fat deposition has been previously described. We amplified fragments of 1457 and 1459 bp that corresponded to the complete coding region and the 3'-untranslated region (UTR), respectively, of the pig ME1 gene. The sequences of these two fragments in pigs from three breeds (Landrace, Large White and Piétrain) contained five single nucleotide polymorphisms (SNP) in the 3'-UTR: C1706T, G1762T, A1807C, C1857A and T1880A. Three haplotypes were found in two generations of a selected Landrace population: H1 (C1706 G1762 A1807 C1857 A1880), H2 (C1706 G1762 A1807 C1857 T1880) and H3 (T1706 T1762 C1807 A1857 T1880). Using Bayesian association analyses, significant associations (highest posterior density at 95%) between ME1 genotype and backfat (BF) thickness at 171 days and muscular pH were found in a Landrace population.

High amino acid variation in the intracellular domain of the pig prolactin receptor (PRLR) and its relation to ovulation rate and piglet survival traits1

Two polymorphisms of the porcine prolactin receptor (PRLR) gene were previously related to litter size by several authors; however, the magnitude and direction of such effects varied depending on the population analyzed. We have sequenced the complete coding region of the porcine PRLR gene and found 6 nonconservative SNP: C1217T (Leu/Pro406), C1283A (Asp/Ala428), G1439A (Lys/Arg480), T1528A (Met/Leu510), G1600A (Gly/Ser534), and G1789A (Gly/Ser597), within exon 10 of the gene, which encodes the entire intracytoplasmic domain of the protein. Eight haplotypes were found and were segregating at different frequencies in 6 porcine breeds. The effects of each individual SNP and haplotype were evaluated in an Iberian x Meishan F2 population using a univariate mixed-inheritance animal model. Significant effects on the number of corpora lutea were found for PRLR haplotypes (P < 0.012), confirming the previously reported associations of PRLR in this process and highlighting the importance of performing analysis of haplotypes rather than of individual SNP. Suggestive effects or tendencies were found for heart rate at birth (P < 0.10), rectal temperature (P < 0.05), and oxygen saturation (P < 0.10) 1 h after birth, and time to first suckle (P < 0.10). We found greater than expected levels of amino acid variability within the intracellular domain of the porcine PRLR, which have been associated with differences in the number of corpus lutea of sows and the preweaning survivability of piglets.

QTL mapping for teat number in an Iberian-by-Meishan pig intercross

The aim of this study was to investigate chromosomal regions affecting the number of teats in pigs and possible epistatic interactions between the identified quantitative trait loci (QTL). An experimental F2 cross between Iberian and Chinese Meishan lines was used for this purpose. A genomic scan was conducted with 117 markers covering the 18 porcine autosomes. Linkage analyses were performed by interval mapping using an animal model to estimate QTL and additive polygenic effects. Complementary analyses with models fitting two QTL were also carried out. The results showed three genomewide significant QTL mapping on chromosomes 5, 10 and 12, whose joint action control up to 30% of the phenotypic variance of the trait. Meishan alleles had a positive additive effect on teat number, and a positive-additive x additive-epistatic interaction was detected between QTL on chromosomes 10 and 12.

Polymorphisms of the porcine dopaminebeta-hydroxylase gene and their relation to reproduction and piglet survivability in an Iberian x Meishan F2 intercross

The goals of this study were to sequence and physically map the porcine dopaminebeta-hydroxylase (DBH) gene, as well as to perform an association study between polymorphisms of this gene and the reproductive performance and piglet survivability of F(2) pigs from an Iberian x Meishan cross. The porcine DBH gene was positioned by RH mapping near the telomere of chromosome 1q2.13, close to markers SSC10D08 and SW1301. Sequencing of DBH cDNAs from 10 pigs revealed the existence of six nucleotide polymorphisms, two of which led to non-synonymous amino acid substitutions within exon 3 at positions 463A>G and 616A>G that corresponded to Thr155Ala and Lys206Glu respectively. Three haplotypes segregated in an Iberian x Meishan population: DBH(X) (A(436)-A(616)), DBH(Y) (A(436)-G(616)) and DBH(Z) (G(436)-G(616)). The DBH haplotypes significantly affected rectal temperatures 1 h after birth (P = 0.002) and have a suggestive effect on the time to first colostrum intake (P = 0.019) and on birth weight (P = 0.019).

Characterization of the porcine acyl-CoA synthetase long-chain 4 gene and its association with growth and meat quality traits

Summary Long-chain acyl-CoA synthetase (ACSL) catalyses the formation of long-chain acyl-CoA from fatty acid, ATP and CoA, activating fatty acids for subsequent reactions. Long-chain acyl-CoA synthetase thus plays an essential role in both lipid biosynthesis and fatty acid degradation. The ACSL4 gene was evaluated as a positional candidate gene for the quantitative trait loci (QTL) located between SW2456 and SW1943 on chromosome X. We have sequenced 4906 bp of the pig ACSL4 mRNA. Sequence analysis allowed us to identify 10 polymorphisms located in the 3'-UTR region and to elucidate two ACSL4 haplotypes. Furthermore, a QTL and an association study between polymorphisms of the ACSL4 gene and traits of interest were carried out in an Iberian x Landrace cross. We report QTL that have not been previously identified, and we describe an association of the ACSL4 polymorphisms with growth and percentage of oleic fatty acid. Finally, we have determined allelic frequencies in 140 pigs belonging to the Iberian, Landrace, Large White, Meishan, Pietrain, Duroc, Vietnamese, Peccary and Babirusa populations.

Effect of the porcine IGF2-intron3-G3072A substitution in an outbred Large White population and in an Iberian × Landrace cross1

The IGF2-intron3-G3072A substitution has been recently described as the causal factor of the imprinted QTL for fat deposition and muscle growth detected within the porcine IGF2 region. The objective of this study was to investigate the IGF2 substitution effect in a Large White outbred population and in an Iberian x Landrace F2 cross. The results showed that the substitution has significant effects on fatness, growth, and shape traits with estimated effects in the expected direction. These results agree with those obtained in the F2 cross, where the IGF2-intron3-G3072A substitution is segregating only in a small family. In addition, a QTL scan has been performed in the F2 population for the traits used in the IGF2 substitution effect validation. Results of this study demonstrated that there are QTL segregating in swine chromosome 2 other than the IGF2 substitution for carcass weight, LM area, and pH measured at 24 h after slaughter. The results confirm the relevance of the IGF2 substitution, but they also show that there are still valuable mutations to be revealed in this chromosome.

The value of prior information for detection of QTL affecting longitudinal traits: an example using Von Bertalanffy growth function

A Bayesian procedure is presented for detecting quantitative trait loci (QTL) affecting longitudinal traits. The statistical model assumes a QTL affecting the prior distribution of the parameters of a given production function, under a hierarchical Bayesian scheme. Marginal posterior distributions for the effects associated with the QTL are calculated using Markov chain Monte Carlo methods. Furthermore, the Bayesian analysis allows the use of some available relevant information that can improve the detection of the QTL substantially. To illustrate the procedure, an example of QTL detection using the Von Bertalanffy growth function is presented with a F2 pig population bred from Iberian boars and Landrace sows. Animals of the F2 population were genotyped for seven markers in chromosome 2 (SSC2). Two prior distributions for the mean effect of the parameters related with birth and adult weight were compared. On the one hand, vague prior distributions were used, and, on the other, there were assumed univariate Gaussian distributions that ensure biologically meaningful adult and birth weights on the posterior growth curves. Results from the second prior distribution supported the presence of QTL, by showing that individuals with both alleles of Iberian origin had lower rates of maturation. On the contrary, when vague priors were used, the procedure was not able to detect QTL.

Fine mapping of porcine chromosome 6 QTL and LEPR effects on body composition in multiple generations of an Iberian by Landrace intercross

The leptin receptor gene (LEPR) is a candidate for traits related to growth and body composition, and is located on SSC6 in a region where fatness and meat composition quantitative trait loci (QTL) have previously been detected in several F2 experimental designs. The aims of this work were: (i) to fine map these QTL on a larger sample of animals and generations (F3 and backcross) of an Iberian×Landrace intercross and (ii) to examine the effects of LEPR alleles on body composition traits. Eleven single nucleotide polymorphisms (SNPs) were detected by sequencing LEPR coding regions in Iberian and Landrace pig samples. Three missense polymorphisms were genotyped by pyrosequencing in 33 F0, 70 F1, 418 F2, 86 F3 and 128 individuals coming from the backcross of four F2 males with 24 Landrace females. Thirteen microsatellites and one SNP were also genotyped. Traits analysed were: backfat thickness at different locations (BFT), intramuscular fat percentage (IMFP), eye muscle area (EMA), loin depth (LOD), weight of shoulder (SHW), weight of ribs (RIBW) and weight of belly bacon (BBW). Different statistical models were applied in order to evaluate the number and effects of QTL on chromosome 6 and the possible causality of the LEPR gene variants with respect to the QTL. The results support the presence of two QTL on SSC6. One, at position 60–100 cM, affects BFT and RIBW. The other and more significant maps in a narrow region (130–132 cM) and affects BFT, IMFP, EMA, LOD, SHW, RIBW and BBW. Results also support the association between LEPR alleles and BFT traits. The possible functional implications of the analysed polymorphisms are considered.

Viability of Iberian x Meishan F2 newborn pigs. I. Analysis of physiological and vitality variables

Eight physiological and vitality variables related to the first minutes of the life of Iberian x Meishan F2 piglets were analyzed. Variables included rectal temperature, arterial oxygen saturation, heart rate (all were monitored at birth and 60 min later), time to reach the udder, and time to the first colostrum intake. Litter effect within boar, a random source of variation, influenced all variables, whereas the random boar effect was not significant for heart rate and arterial oxygen saturation at birth and the time to reach the udder. Birth weight influenced rectal temperature at birth (P < 0.01) and 60 min later (P < 0.001), arterial oxygen saturation 1 h after birth (P < 0.05), heart rate 1 h after birth (P < 0.001), time to reach the udder (P < 0.05), and time to the first suckle (P < 0.01), whereas viability score influenced rectal temperatures (at birth and 60 min later; P < 0.001), and time to reach the udder (P < 0.001) and to suckle (P < 0.001). Finally, the order of birth showed significant effects for rectal temperature 1 h after birth (P < 0.001) and time to first colostrum intake (P < 0.001). Correlation coefficients between physiological variables were generally low, with the exception of the ones for rectal temperature 1 h after birth and arterial oxygen saturation 1 h after birth (0.38; P < 0.001), and for rectal temperature at birth and rectal temperature 1 h after birth (0.34; P < 0.001). Times to reach the udder and to suckle were highly correlated (0.67; P < 0.001) and also were moderately and negatively correlated with rectal temperature 1 h after birth (-0.36 and -0.38 respectively; P < 0.001). Heart rates at birth and 1 h later, as well as arterial oxygen saturation at birth and 1 h later, were not correlated, showing that the values at birth do not necessarily provide information about the physiological status of the pig 1 h after birth. The recording of physiological variables, birth weight, and viability score may be useful to identify weak piglets quickly and to establish palliative measures.

A simulation study on the detection of causal mutations from F2 experiments

A simulation study has been performed to evaluate the power and the rate of false positives for the detection of causal mutations under two different models of analysis. We used an F2 design generated from an F0 population of five sires of line 1 and 40 dams of line 2 to produce an F1 population of 10 sires and 80 dams. Two different locations of the causal mutation and several frequencies of the mutations in the parental populations were considered. The first model included only the genetic configuration of the mutation, while the second model also included the probability of line origin given the neutral markers. Both models performed well when the mutation at the candidate gene was the causal mutation, although a greater power was obtained using the first model, because of its relative simplicity compared to the second one. However, when the candidate gene mutation was a neutral mutation, the second model presented a lower rate of false positives than the first. Moreover, in some cases the second model allowed distinction between the neutral and the causal mutation. The F2 design has a great power to detect quantitative trait loci provided by linkage disequilibrium, but also makes it difficult to discriminate between causal and neutral mutations. Therefore a high percentage of false positives can be expected. The limitations of F2 designs for discriminating between neutral and causal mutations are discussed.

Derivation of a Bayes factor to distinguish between linked or pleiotropic quantitative trait loci

A simple procedure to calculate the Bayes factor between linked and pleiotropic QTL models is presented. The Bayes factor is calculated from the marginal prior and posterior densities of the locations of the QTL under a linkage and a pleiotropy model. The procedure is computed with a Gibbs sampler, and it can be easily applied to any model including the location of the QTL as a variable. The procedure was compared with a multivariate least-squares method. The proposed procedure showed better results in terms of power of detection of linkage when low information is available. As information increases, the performance of both procedures becomes similar. An example using data provided by an Iberian by Landrace pig intercross is presented. The results showed that three different QTL segregate in SSC6: a pleiotropic QTL affects myristic, palmitic, and eicosadienoic fatty acids; another pleiotropic QTL affects palmitoleic, stearic, and vaccenic fatty acids; and a third QTL affects the percentage of linoleic acid. In the example, the Bayes factor approach was more powerful than the multivariate least-squares approach.

Viability of Iberian × Meishan F2 newborn pigs. II. Survival analysis up to weaning1

Iberian x Meishan F2 piglet's preweaning survivability was analyzed using categorical data regression procedures within the proportional hazards assumption. A frailty sire model was assumed with the litter effect treated as an additional random source of variation. Moreover, the relative birth weight within litter and the litter effect were considered time-dependent covariates that changed their values in the second day of life due to cross fostering carried out to standardize litters. Six variables had a significant effect on survivability: birth weight (P < 0.01), relative birth weight within litter (P < 0.001), rectal temperature 60 min after birth (P < 0.01), type of presentation at birth (P < 0.05), presence of stillbirths (P < 0.001), and presence of mummified fetuses (P < 0.001). Small piglets (<0.98 kg) suffered a high hazard ratio (6.57; P < 0.001), with this variable being clearly lower for the rest of birth weight categories. Piglets that were small in relation to their siblings (relative birth weight within litter) also suffered an increased death risk, with a hazard ratio of 1.81 (P < 0.05), which was similar to animals with posterior presentations at birth (hazard ratio = 1.80; P < 0.05). Piglets with a rectal temperature lower than 35.4 degrees C 60 min after birth showed the highest hazard ratio (7.18; P < 0.01). Furthermore, the presence of mummified fetuses decreased the survivability of the remaining siblings, with a hazard ratio of 2.03 (P < 0.01), as did the presence of stillbirths (hazard ratio = 3.55; P < 0.001). The inclusion of the two random effects allowed us to estimate the mode of the joint posterior density of the sire variance (0.08) and the litter variance (1.98). The estimated heritability of preweaning survival reached a value of 0.03. We conclude that piglet survival involves several systematic influences related to birth weight, thermoregulatory ability, and injuries suffered during gestation and farrowing. The genetic variance was small compared with those generated by the common environment, for which the genetic improvement of piglet survival seems difficult.

Derivation of a Bayes Factor to Distinguish Between Linked or Pleiotropic Quantitative Trait Loci

A simple procedure to calculate the Bayes factor between linked and pleiotropic QTL models is presented. The Bayes factor is calculated from the marginal prior and posterior densities of the locations of the QTL under a linkage and a pleiotropy model. The procedure is computed with a Gibbs sampler, and it can be easily applied to any model including the location of the QTL as a variable. The procedure was compared with a multivariate least-squares method. The proposed procedure showed better results in terms of power of detection of linkage when low information is available. As information increases, the performance of both procedures becomes similar. An example using data provided by an Iberian by Landrace pig intercross is presented. The results showed that three different QTL segregate in SSC6: a pleiotropic QTL affects myristic, palmitic, and eicosadienoic fatty acids; another pleiotropic QTL affects palmitoleic, stearic, and vaccenic fatty acids; and a third QTL affects the percentage of linoleic acid. In the example, the Bayes factor approach was more powerful than the multivariate least-squares approach.

Full pedigree quantitative trait locus analysis in commercial pigs using variance components

In commercial livestock populations, QTL detection methods often use existing half-sib family structures and ignore additional relationships within and between families. We reanalyzed the data from a large QTL confirmation experiment with 10 pig lines and 10 chromosome regions using identity-by-descent (IBD) scores and variance component analyses. The IBD scores were obtained using a Monte Carlo Markov Chain method, as implemented in the LOKI software, and were used to model a putative QTL in a mixed animal model. The analyses revealed 61 QTL at a nominal 5% level (out of 650 tests). Twenty-seven QTL mapped to areas where QTL have been reported, and eight of these exceeded the threshold to claim confirmed linkage (P < 0.01). Forty-two of the putative QTL were detected previously using half-sib analyses, whereas 46 QTL previously identified by half-sib analyses could not be confirmed using the variance component approach. Some of the differences could be traced back to the underlying assumptions between the two methods. Using a deterministic approach to estimate IBD scores on a subset of the data gave very similar results to LOKI. We have demonstrated the feasibility of applying variance component QTL analysis to a large amount of data, equivalent to a genome scan. In many situations, the deterministic IBD approach offers a fast alternative to LOKI.

Identification of quantitative trait loci for production traits in commercial pig populations

The aim of this study was to investigate methods for detecting QTL in outbred commercial pig populations. Several QTL for back fat and growth rate, previously detected in experimental resource populations, were examined for segregation in 10 different populations. Two hundred trait-by-population-by-chromosome tests were performed, resulting in 20 tests being significant at the 5% level. In addition, 53 QTL tests for 11 meat quality traits were declared significant, using a subset of the populations. These results show that a considerable amount of phenotypic variance observed in these populations can be explained by major alleles segregating at several of the loci described. Thus, despite a relatively strong selection pressure for growth and back fat traits in these populations, these alleles have not yet reached fixation. The approaches used here demonstrate that it is possible to verify segregation of QTL in commercial populations by limited genotyping of a selection of informative animals. Such verified QTL may be directly exploited in marker-assisted selection (MAS) programs in commercial populations and their molecular basis may be revealed by positional candidate cloning.

Reallocation of body resources in lactating mice highly selected for litter size

The present study investigated differences in the allocation patterns of body stores in lactating female mice from a line selected for high litter size at birth (S-line, average litter size of 20) and dams from a nonselected control line (C-line, average litter size of 10). Body weight, litter size, litter weight, and absolute and relative lipid and protein mass were measured at peak lactation (2 wk in lactation) and at weaning (3 wk in lactation). Body size in S-line females has been increased as a correlated effect of selection for high litter size at birth, allowing for larger litters and higher absolute milk production. However, these dams produce larger litters relative to their own body weight. At peak lactation, lipid and protein percentage did not differ between lines. At weaning, S-line females had a higher protein percentage (P < 0.001) and lower lipid percentage (P < 0.05) than C-line females. Apparently, S-line females produce more offspring but at a greater cost to their own metabolism. This process was insufficient to supply the offspring with adequate resources, resulting in reduced (P < 0.0001) pup development and increased (P < 0.0001) preweaning mortality rates.

Quantitative trait locus mapping for meat quality traits in an Iberian x Landrace F2 pig population

An experimental F2 cross between Iberian and Landrace pig strains was performed to map quantitative trait loci (QTL) for diverse productive traits. Here we report results for meat quality traits from 369 F2 animals with records for pH 24 h postmortem (pH 24 h), muscle color Minolta measurements L* (lightness), a* (redness), and b* (yellowness), H* (hue angle), C* (chroma), intramuscular fat (IMF) and haematin pigment content measured in the longissimus thoracis. Pigs were genotyped for 92 markers covering the 18 porcine autosomes (SSC). Results of the genome scan show evidence for QTL for IMF (SSC6; F = 27.16), pH 24 h (SSC3; F = 7.73), haematin pigments (SSC4 and SSC7; F = 8.68 and 9.47 respectively) and Minolta color measurements L* (SSC4 and SSC7; F =16.42 and 7.17 respectively), and a* (SSC4 and SSC8; F = 8.05 and 7.36 respectively). No QTL were observed for the color measurements b*, H*, and C*. Alternative models fitting epistasis between QTL were also tested, but detected epistatic interactions were not significant at a genome-wise level. In this work we identify genomic regions related with meat quality traits. Improvement by traditional selection methods is complicated, and finer mapping would be required for their application in introgression programs.