A novel QTL region for pH and meat color in Duroc pigs

One of the most important processes that occur during the transformation of muscle to meat is the pH decline as a consequence of the post-mortem metabolism of muscle tissue. Abnormal pH declines lead to pork defects such as pale, soft, and exudative meat. There is genetic variance for ultimate pH and the role of some genes on this phenotype is well established. After conducting a genome-wide association study on ultimate pH using 526 purebred Duroc pigs, we identified associated regions on Sus scrofa chromosomes (SSC) 3, 8, and 15. Functional candidate genes in these regions included PRKAG3 and PHKG1. The SSC8 region, at 71.6 Mb, was novel and, although no candidate causative gene could be identified, it may have regulatory effects. Subsequent analysis on 828 pigs from the same population confirmed the impact of the three associated regions on pH and meat color. We detected no interaction between the three regions. Further investigations are necessary to unravel the functional significance of the novel genomic region at SSC8. These variants could be used as markers in marker-assisted selection for improving meat quality.

Identification of breed-specific genomic variants in Colombian Creole pig breeds by whole-genome sequencing

Dissecting genetic variation of local breeds is important for the success of conservation. In this research, we investigated the genomic variation of Colombian Creole (CR) pigs, with a focus on the breed-specific variants in the exonic region of 34 genes with reported effects on adaptive and economic traits. Seven individuals of each of the three CR breeds (CM, Casco de Mula; SP, San Pedreño; and ZU, Zungo) were whole-genome sequenced along with 7 Iberian (IB) pigs and 7 pigs of each of the four most used cosmopolitan (CP) breeds (Duroc, Landrace × Large White, and Pietrain). Molecular variability in CR (6,451,218 variants; from 3,919,242, in SP, to 4,648,069, in CM) was comparable to that in CP, but higher than in IB. For the investigated genes, SP pigs displayed less exonic variants (178) than ZU (254), CM (263), IB (200), and the individual CP genetic types (201 to 335). Sequence variation in these genes confirmed the resemblance of CR to IB and indicates that CR pigs, particularly ZU and CM, are not exempt from selective introgression of other breeds. A total of 50 exonic variants were identified as being potentially specific to CR, including a high-impact deletion in the intron between exons 15 and 16 of the leptin receptor gene, which was only found in CM and ZU. The identification of breed-specific variants in genes related to adaptive and economical traits can bolster the understanding of the role of gene-environment interactions on local adaptation and points the way for effective breeding and conservation of CR pigs.

The leptin receptor gene affects piglet behavior and growth

Piglets with low birth weight present low vitality after farrowing, often leading to impaired weight gain during lactation. A recessive missense variant (C > T) for increased appetite and fatness in the porcine leptin receptor gene (rs709596309) causes a negative maternal effect on the weight of piglets at weaning. However, it is not known whether this variant already exerts an effect on the birth weight and vitality of newborn piglets and on their growing capacity during lactation. An experiment was conducted using 668 purebred Duroc piglets (131 CC, 311 CT, and 226 TT) from 74 multiparous sows (9 CC, 43 CT, and 22 TT) and 14 boars (1 CC, 10 CT, and 3 TT). All piglets were individually weighed at birth and tested for vitality, which was assessed on a scale from 1 (low vitality) to 3 (high vitality) based on behavioral observations, including the status of the piglet immediately before the test. Only non-adopted piglets were considered for piglet performance at weaning. Inferences on the effect of the genotype on birth and weaning traits were done on a Bayesian setting based on 2-trait bivariate models including the effects of the piglet and the litter, as well as the genotype of the sow and the piglet, the sex of the piglet, and the parity number. Vitality and the status of the piglet before the test were analyzed using a liability threshold (probit) model. As compared to other genotypes, TT newborn piglets were 28 g heavier, were more vital (the probability of being scored as highly vital was 6.5% higher) and were more often found suckling before the test (the probability of being suckling at test was 6.5% higher). As a result, they grew more during lactation (153 g) and were heavier at weaning (169 g) than littermates of the two other genotypes, thus partly compensating for the limited maternal capacity of TT sows. Our findings provide evidence that appetite-influencing genes, such as the leptin receptor gene, have developmental implications from very early life stages.

Identification of a missense variant in the porcine AGPAT gene family associated with intramuscular fat content through whole-genome sequencing

The 1-acylglycerol-3-phosphate O-acyltransferases (AGPATs) are enzymes that catalyze the conversion of lysophosphatidic acid to phosphatidic acid, which is a precursor of triacylglycerol, the main fat reservoir in mammals. We used whole-genome sequencing of 205 pigs to identify 6639 genetic variants in the porcine AGPAT gene family. Of these, 166 common variants in the AGPAT5 gene had significant associations with fat content and composition traits. We preselected a missense single nucleotide polymorphism in exon 6 of AGPAT5 (rs196952262, A>G) for validation of its associations in 1034 pigs from the same Duroc line. The A allele showed a positive additive effect for intramuscular fat content (+1.12% ± 0.21, p < 0.001, for gluteus medius and +0.89% ± 0.33, p < 0.01, for longissimus). We also observed significant associations with fatty acid composition that were, at least in part, independent of the increased intramuscular fat. The A allele resulted in more monounsaturated fatty acids (+0.34% ± 0.15, p < 0.05, for longissimus) and a greater monounsaturated/polyunsaturated fatty acids ratio (+0.11 ± 0.04, p < 0.01, for gluteus medius and +0.13 ± 0.05, p < 0.05, for longissimus). The effect of the AGPAT5 variant on intramuscular fat was more noticeable in fatter pigs, and AGPAT5 interacts with other genes that affect overall fatness such as LEPR. AGPAT5 was the most expressed gene of the AGPAT family in pig skeletal muscle. This variant can be used as a marker in assisted selection for modulating pig fat deposition and fatty acid content.

Age Evolution of Lipid Accretion Rate in Boars Selected for Lean Meat and Duroc Barrows

Fatty acid (FA) deposition in growing–fattening pigs is mainly based on endogenous lipid synthesis, but also direct FA incorporation from the diet. To evaluate the direct fat incorporation rates and the endogenous desaturation action of the stearoyl-CoA desaturase (SCD) enzyme, a deuterium (D)-labeled saturated FA (d35-C18:0) was added to the diet. Sixteen three-way (3W) crossbred boars, and thirty-two purebred Duroc barrows homozygous for the SCD single nucleotide polymorphism rs80912566 (16 CC/16 TT), were used. Half of the animals of each genotype belonged to the growing and fattening phases. The fractional incorporation rate (FIR) of dietary fat in growing pigs was generally higher in adipose tissues, whereas in fattening pigs it was higher in the liver. Duroc pigs exhibited lower FIRs than 3W pigs, suggesting lower rates of endogenous synthesis by 3W pigs. Real fractional unsaturation rates (FURs) increased with age by the higher FIRs in 3W pigs and the de novo synthesis pathway in Duroc genotypes. Moreover, pigs carrying the SCD_T allele showed more enhanced oleic acid biosynthesis than Duroc CC pigs. In conclusion, suitable feeding protocols should be designed for each pig type to optimize production traits, considering that the metabolic pathway of FA for its deposition may differ.

Transcriptome shifts triggered by vitamin A and SCD genotype interaction in Duroc pigs

BACKGROUND

The composition of intramuscular fat depends on genetic and environmental factors, including the diet. In pigs, we identified a haplotype of three SNP mutations in the stearoyl-coA desaturase (SCD) gene promoter associated with higher content of monounsaturated fatty acids in intramuscular fat. The second of these three SNPs (rs80912566, C > T) affected a putative retinol response element in the SCD promoter. The effect of dietary vitamin A restriction over intramuscular fat content is controversial as it depends on the pig genetic line and the duration of the restriction. This study aims to investigate changes in the muscle transcriptome in SCD rs80912566 TT and CC pigs fed with and without a vitamin A supplement during the fattening period.

RESULTS

Vitamin A did not affect carcass traits or intramuscular fat content and fatty acid composition, but we observed an interaction between vitamin A and SCD genotype on the desaturation of fatty acids in muscle. As reported before, the SCD-TT pigs had more monounsaturated fat than the SCD-CC animals. The diet lacking the vitamin A supplement enlarged fatty acid compositional differences between SCD genotypes, partly because vitamin A had a bigger effect on fatty acid desaturation in SCD-CC pigs (positive) than in SCD-TT and SCD-TC animals (negative). The interaction between diet and genotype was also evident at the transcriptome level; the highest number of differentially expressed genes were detected between SCD-TT pigs fed with the two diets. The genes modulated by the diet with the vitamin A supplement belonged to metabolic and signalling pathways related to immunity and inflammation, transport through membrane-bounded vesicles, fat metabolism and transport, reflecting the impact of retinol on a wide range of metabolic processes.

CONCLUSIONS

Restricting dietary vitamin A during the fattening period did not improve intramuscular fat content despite relevant changes in muscle gene expression, both in coding and non-coding genes. Vitamin A activated general pathways of retinol response in a SCD genotype-dependant manner, which affected the monounsaturated fatty acid content, particularly in SCD-CC pigs.

A Methodology to Quantify Resilience in Growing Pigs

Simple Summary

The concept of resilience can be defined as the ability of animals to maintain their productivity despite the stressors that might occur during the rearing period. Despite the growing interest in the genetic determinism of resilience and its possible inclusion in selection objectives, there is no straightforward method to measure this trait. Several quantifying methods estimate resilience based on productivity such as body weight or milk production, or non-productivity related traits such as immunity or stress responses. The objective of this study was to elaborate novel resilience indicators in growing pigs based on both productivity (body weight) and non-productivity (acute-phase proteins) related traits. We propose the deviation from the expected growth curve and the increment of the acute-phase protein haptoglobin, after applying a common vaccine, as resilience indicators in growing pigs under standard rearing conditions. We showed that the suggested resilience indicators are under a genetic control, show a substantial variability in the population, and may be improved through selection.

Abstract

There is a growing concern about the genetic determinism of resilience and its possible implementation in breeding programs. The objective of our study was to elaborate novel resilience indicators in growing pigs based on the deviation from the expected growth curve and the increment of the acute-phase protein haptoglobin (HP) after applying a common vaccine. A total of 445 pigs were vaccinated with an attenuated Aujeszky vaccine at 12 weeks of age. Deviation from the expected body weight (ΔBW) given the growth curve of unvaccinated pigs at 28 days post-vaccination (DPV) and the increment of HP at 4 DPV (ΔHP) were suggested as resilience indicators. Challenged pigs that maintained their productivity and had a minor activation of HP were deemed resilient, whereas pigs that had low ∆BW values and a high activation of HP were deemed susceptible. Pigs were also classified based on ∆BW and ∆HP relative to the expected BW at 28 DPV and to the basal level of HP, respectively. The concordance was high between both methods, indicating that ΔBW and ΔHP are not sensitive to the animal’s expected BW nor the basal level of HP. The heritability estimates were moderate for ∆BW (0.33) and low-to-moderate for ∆HP (0.16). Our study suggests ΔBW and ΔHP as novel resilience indicators in pigs. The suggested indicators capture different aspects of resilience, are easy to measure, and are genetically controlled. Thus, they may be improved through selective breeding. Further analyses are needed to validate our findings.

Antagonistic maternal and direct effects of the leptin receptor gene on body weight in pigs

Maternal effects on offspring growth can impact survival and evolution of natural and domesticated populations. Genetic correlation estimates often support a negative relationship between direct and maternal effects. However, the genetic underpinnings whereby this antagonism operates are unclear. In pigs, sow feeding status and body composition condition piglet development and growth. We hypothesized that variants in genes impacting these traits may be causative of maternal influences that could be antagonistic to the direct effects for piglet growth. A recessive missense mutation (C>T) in the porcine leptin receptor (LEPR) gene (rs709596309) has been identified as the possible causal polymorphism for increased feed intake and fatness. Using data from a Duroc line, we show that the TT sows exerted a negative impact on the body weight of their offspring at the end of the growing period of similar extent to the positive direct effect of the TT genotype over each individual. Thus, TT pigs from TT dams were about as heavy as CC and CT (C–) pigs from C–dams, but TT pigs from C–dams were around 5% heavier than C–pigs from TT dams. In contrast, body composition was only influenced by LEPR direct effects. This antagonism is due to a higher propensity of TT pigs for self-maintenance rather than for offspring investment. We show that TT pigs consumed more feed, favored fatty acid uptake over release, and produced lighter piglets at weaning than their C–counterparts. We conclude that LEPR underlies a transgenerational mechanism for energy distribution that allocates resources to the sow or the offspring according to whether selective pressure is exerted before or after weaning.

A probabilistic Poisson-based model to detect PRRSV recirculation using sow production records

Porcine reproductive and respiratory syndrome (PRRS) is a viral disease associated with a decrease in the number of born alive piglets (NBA) and an increase in the number of lost piglets (NLP) per farrowing. Under practical conditions, it is critical to assess whether a farm is suffering PRRSV recirculation in the sow herd as soon as possible. The aim of this research work was to develop a new method to detect potential PRRSV recirculation in sow production farms. Sow reproductive performance records from one farm (farm T) were used to set up the method and records from ten additional farms (farms V1 to V10) were used for validation. A conditional Poisson model of NLP on NBA was proposed to fit the data. A three-step procedure was implemented to detect potential PRRSV recirculation: (i) computation of the maximum-likelihood estimates of the expected values of NBA and NLP in a PRRSV non-recirculating scenario; (ii) calculation, for each farrowing, of the p-value associated with the probability of jointly observing deviations towards decreased NBA and increased NLP. The detection of a potential PRRSV recirculation was based on (iii) the combined p-value resulting from weighing the p-values of the last N farrowings by the chi-square-inverse method. In order to gain specificity, a displacement on the expected non-recirculating NBA and NLP values was used for tuning purposes. With this approach, two PRRSV circulating periods were detected in farm T, which were confirmed with standard laboratorial diagnostic techniques. The method was subsequently validated in farms V1 to V10, where ten PRRSV-recirculating time episodes had been diagnosed. The method proposed here was able to detect the ten PRRSV recirculations using a relatively small set of contiguous farrowings, with only two mismatched weeks, one as a false negative, in farm V1, and one as a false positive, in farm V4. It is concluded that a conditional Poisson-based model of NLP on NBA can be a useful tool for routinely detecting PRRSV recirculation in sow herds.

Identification of resilient sows in Porcine Reproductive and Respiratory Syndrome virus infected farms

The identification of resilient sows can improve reproductive performance in farms exposed to multiple challenges. A common challenge is the porcine reproductive and respiratory syndrome virus (PRRSV). A key issue to deal with disease resilience is to set up a feasible phenotyping strategy. Our aim was to develop a phenotyping criterion to discriminate susceptible from resilient sows in PRRSV-infected farms. A total of 517 Landrace x Large White gilts were classified as resilient (R) or susceptible (S) to PRRSV virus, following vaccination with MLV-PRRSV at 6-7 wk of age, in a PRRSV negative multiplication farm. Female piglets were phenotyped as R if their serum was negative to PRRSV at 7 and 21 d post-vaccination (DPV) or as S if their serum was positive at 7 and/or 21 DPV. Amongst them, 382 gilts were transferred to a PRRSV-positive production farm, where the number of piglets born alive (NBA), stillborn (NSB), mummified (NMU), lost (NLP=NSB+NMU) and total born (NTB = NBA+NLP) were recorded for almost three years. Data were collected during two periods according to the PRRSV farm health status, which were confirmed as either PRRSV-positive stable (endemic) or inestable (epidemic). Analyses were carried out under a Bayesian approach. The heritability for the resilience criterion was estimated using a threshold model. A linear (for NTB and NBA) and a binomial model (for NSB, NMU and NLP) on the resilience criterion by the farm health status were used to assess the difference between R and S sows. The heritability of the resilience criterion was 0.46 (SD 0.06). The probability of a piglet being lost was greater (≥0.97) in S than in R litters, regardless of whether the delivery occurred during a PRRSV outbreak (20.5% vs 17.0%) or not (15.8% vs 13.7%). The lower piglet mortality rate in R sows was due to NSB, in the endemic phase (13.0% vs 15.0% of NTB, with a posterior probability of 98% of S sows showing higher NSB than R sows), and to NMU, in the epidemic phase (4.0% vs 8.4% of NTB, with a posterior probability of >99% of S sows showing higher NMU than R sows). During a PRRSV outbreak, the S sows were twice as likely to give birth to a mummified piglet as compared to R sows. These findings provide evidence that the described phenotyping scheme has a potential use as a PRRSV resilience criterion.

Linoleic acid metabolic pathway allows for an efficient increase of intramuscular fat content in pigs

Background

Intramuscular fat (IMF) content is a relevant trait for high-quality meat products such as dry-cured ham, but increasing IMF has the undesirable correlated effect of decreasing lean growth. Thus, there is a need to find selection criteria for IMF independent from lean growth. In pigs, the proportion of linoleic (C18:2) and arachidonic (C20:4) acids decline with fat deposition and therefore they can be considered as indicators of fatness. The aim of this research was to estimate the genetic variation for C18:2 and C20:4 in IMF and their genetic correlations with IMF and lean growth traits, with the objective to assess their potential as specific biomarkers of IMF. The analysis was conducted using a full-pedigreed Duroc resource line with 91,448 records of body weight and backfat thickness (BT) at 180 days of age and 1371 records of fatty acid composition in the muscle gluteus medius.

Results

The heritability estimates for C18:2 and C20:4 in IMF, whether expressed in absolute (mg/g of muscle) or in relative (mg/g of fatty acid) terms, as well as for their ratio (C20:4/C18:2), were high (> 0.40), revealing that the C18:2 to C20:4 pathway is subjected to substantial genetic influence. Litter effects were not negligible, with values ranging from 8% to 15% of the phenotypic variance. The genetic correlations of C18:2 and C20:4 with IMF and BT were negative (- 0.75 to - 0.66, for IMF, and - 0.64 to - 0.36, for BT), if expressed in relative values, but almost null (- 0.04 to 0.07), if expressed in absolute values, except for C18:2 with IMF, which was highly positive (0.88). The ratio of C20:4 to C18:2 also displayed a stronger genetic correlation with IMF (- 0.59) than with BT (- 0.10).

Conclusions

The amount of C18:2 in muscle can be used as an IMF-specific biomarker. Selection for the absolute amount of C18:2 is expected to deliver a similar response outcome as selection for IMF at restrained BT. Further genetic analysis of the C18:2 metabolic pathway may provide new insights into differential fat deposition among adipose tissues and on candidate genes for molecular markers targeting specifically for one of them.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

Background

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

Results

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110–114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

Conclusions

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5557-9) contains supplementary material, which is available to authorized users.

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.

TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Genetics and breeding for intramuscular fat and oleic acid content in pigs

The intramuscular fat (IMF) and oleic acid (OL) content have been favorably related to pork quality and human health. This influences the purchasing behavior of consumers and, therefore, also shifts the attention of breeding companies toward whether these traits are included into the breeding goal of the lines producing for high-valued markets. Because IMF and OL are unfavorably associated with lean content, a key economic trait, the real challenge for the industry is not simply to increase IMF and OL, but rather to come up with the right trade-off between them and lean content. In this paper we review the efforts performed to genetically improve IMF and OL, with particular reference to the research we conducted in a Duroc line aimed at producing high quality fresh and dry-cured pork products. Based on this research, we conclude that there are selection strategies that lead to response scenarios where IMF, OL, and lean content can be simultaneously improved. Such scenarios involve regular recording of IMF and OL, so that developing a cost-efficient phenotyping system for these traits is paramount. With the economic benefits of genomic selection needing further assessment in pigs, selection on a combination of pedigree-connected phenotypes and genotypes from a panel of selected genetic markers is presented as a suitable alternative. Evidence is provided supporting that at least a polymorphism in the leptin receptor and another in the stearoyl-CoA desaturase genes should be in that panel. Selection for IMF and OL results in an opportunity cost on lean growth. The extent to which it is affordable relies on the consumers' willingness to pay for premium products and on the cost to benefit ratio of alternative management strategies, such as specific dietary manipulations. How the genotype can influence the effect of the diet on IMF and OL remains a topic for further research.

Carotenoid intake and SCD genotype exert complementary effects over fat content and fatty acid composition in Duroc pigs

Nutritional and genetic strategies are needed to enhance intramuscular fat (IMF) and MUFA content without altering carcass leanness. Dietary vitamin A restriction has been suggested to specifically promote IMF, whereas a polymorphism of the () gene has shown to specifically increase MUFA. The purpose of this study was to investigate the combined effects of provitamin A (PVA) carotenoid intake and genotype (>) on hepatic retinoid content and on the liver, muscle (LM and gluteus medius [GM]), and subcutaneous fat (SF) content and fatty acid composition. Following a split-plot design, 32 castrated Duroc pigs, half of each of the 2 homozygous genotypes (CC and TT), were subjected from 165 to 195 d of age to 2 finishing diets differing in the PVA carotenoid content (an enriched-carotene diet [C+] and a control diet [C-]). Both diets were identical except for the corn line used in the feed. The C+ was formulated with 20% of a carotenoid-fortified corn (M37W-Ph3) whereas the C- instead used 20% of its near isogenic M37W line, which did not contain PVA carotenoids. No vitamin A was added to the diets. The C- was estimated to provide, at most, 1,300 IU of vitamin A/kg and the C+ to supply an extra amount of at least 800 IU vitamin A/kg. Compared with the pigs fed the C-, pigs fed with C+ had 3-fold more retinoic acid ( < 0.01) and 4-fold more gene expression in the liver ( = 0.06). The diet did not affect performance traits and backfat thickness, but pigs fed the C+ had less fat (4.0 vs. 5.0%; = 0.07) and MUFA (18.3 vs. 22.5%; = 0.01) in the liver, less IMF (5.4 vs. 8.3%; = 0.04) in the GM, and more fat content (90.4 vs. 87.9%; = 0.09) and MUFA (48.0 vs. 46.6%; = 0.04) in SF. The TT genotype at the gene increased MUFA ( < 0.05) in all tissues (21.4 vs. 19.5% in the liver, 55.0 vs. 53.1% in the LM, 53.9 vs. 51.7% in the GM, and 48.0 vs. 46.7% in SF for TT and CC genotypes, respectively). Liver fat and MUFA content nonlinearly declined with liver all- retinoic acid, indicating a saturation point at relatively low all- retinoic acid content. The results obtained provide evidence for a complementary role between dietary PVA and genotype, in the sense that the TT pigs fed with a low-PVA diet are expected to show higher and more monounsaturated IMF without increasing total fat content.

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH₂₄), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH₂₄ SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Identification of signatures of selection for intramuscular fat and backfat thickness in two Duroc populations

Intramuscular fat (IMF) content is an important trait affecting the quality of pork. Two Duroc populations, one under positive selection for IMF and the other selected for decreased backfat but under stabilizing selection for IMF, were used to identify signatures of selection associated with IMF using 60,000 single-nucleotide polymorphism data. The effects of selection were analyzed between 2 lines or groups representing selected and control animals within each population using a discriminant analysis of principal components and Wright’s fixation index (FST). Moreover, extended haplotype homozygosity-based approaches were used to examine the changes in haplotype frequency due to recent selection. Each statistical method identified 10–20 selection signatures. A few haplotype-based signatures of selection agreed with results from a genome-wide association study (GWAS), while FST measures showed a better agreement with GWAS results. Agreement of marker-trait associations and signatures of selection was limited, and further examination will be necessary to understand the effect of selection on IMF and why some regions identified by GWAS did not appear to respond to the selection practiced. The genes in 21 consensus selection signatures were examined. Several genes with an effect on overall fatness were identified, but further research is needed to assess whether or not some of them could have a specific effect on IMF.

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.

The effect of SCD and LEPR genetic polymorphisms on fat content and composition is maintained throughout fattening in Duroc pigs

The effects of the stearoyl-CoA desaturase (SCD; AY487830:g.2228T>C) and leptin receptor (LEPR; NM_001024587:g.1987C>T) polymorphisms on fat content and fatty acid (FA) composition were investigated throughout fattening. Samples of Longissimus thoracis (LT) and subcutaneous fat (SF) from 214 Duroc barrows were collected from 160days to slaughter age (220days) using a longitudinal design. Results indicated that the positive effect of the T allele at the SCD gene on monounsaturated FA and of the T allele at the LEPR gene on saturated FA are maintained throughout the growing-finishing period, both in LT and SF. In LEPR, however, compositional changes, particularly in SF, are a result of increased fatness. There is very limited evidence of genotype by age interaction, and thus it is concluded that the combined selection for the SCD T and LEPR C alleles is a good strategy to increase the MUFA/SFA ratio regardless of the age at slaughter.

Genome-Wide Association Study Singles Out SCD and LEPR as the Two Main Loci Influencing Intramuscular Fat Content and Fatty Acid Composition in Duroc Pigs

Intramuscular fat (IMF) content and fatty acid composition affect the organoleptic quality and nutritional value of pork. A genome-wide association study was performed on 138 Duroc pigs genotyped with a 60k SNP chip to detect biologically relevant genomic variants influencing fat content and composition. Despite the limited sample size, the genome-wide association study was powerful enough to detect the association between fatty acid composition and a known haplotypic variant in SCD (SSC14) and to reveal an association of IMF and fatty acid composition in the LEPR region (SSC6). The association of LEPR was later validated with an independent set of 853 pigs using a candidate quantitative trait nucleotide. The SCD gene is responsible for the biosynthesis of oleic acid (C18:1) from stearic acid. This locus affected the stearic to oleic desaturation index (C18:1/C18:0), C18:1, and saturated (SFA) and monounsaturated (MUFA) fatty acids content. These effects were consistently detected in gluteus medius, longissimus dorsi, and subcutaneous fat. The association of LEPR with fatty acid composition was detected only in muscle and was, at least in part, a consequence of its effect on IMF content, with increased IMF resulting in more SFA, less polyunsaturated fatty acids (PUFA), and greater SFA/PUFA ratio. Marker substitution effects estimated with a subset of 65 animals were used to predict the genomic estimated breeding values of 70 animals born 7 years later. Although predictions with the whole SNP chip information were in relatively high correlation with observed SFA, MUFA, and C18:1/C18:0 (0.48–0.60), IMF content and composition were in general better predicted by using only SNPs at the SCD and LEPR loci, in which case the correlation between predicted and observed values was in the range of 0.36 to 0.54 for all traits. Results indicate that markers in the SCD and LEPR genes can be useful to select for optimum fatty acid profiles of pork.

Expression profiling of the GBP1 gene as a candidate gene for porcine reproductive and respiratory syndrome resistance

A genomic region in pig chromosome 4 has been previously associated with higher viraemia levels and lower weight gain following porcine reproduction and respiratory syndrome virus (PRRSV) infection. The region includes the marker WUR1000125, a G>A polymorphism next to a putative polyadenylation site in the 3'-untranslated region (3'-UTR) of the guanylate-binding protein 1, interferon-induced (GBP1) gene. The protein encoded by GBP1 is a negative regulator of T-cell responses. We show here that GBP1 expression is lower in liver and tonsils of pigs carrying the WUR1000125-G allele due to differential allele expression (allele A expression is 1.9-fold higher than for allele G). We also show that the GBP1 gene has two active polyadenylation signals 421 bp apart and that polyadenylation usage is dependent on the WUR1000125 genotype. The distal site is the most prevalently used in all samples, but the presence of the A allele favours the generation of shorter transcripts from the proximal site. This is confirmed by a differential allele expression study in AG genotype liver and tonsil samples. The interaction between WUR1000125 and other mutations identified in the 5'- and 3'-UTR regions of this gene needs to be studied. In conclusion, our study indicates that the WUR1000125 mutation is associated with changes in the expression of the negative T-cell regulator GBP1 gene. However, the chromosome 4 locus for PRRSV viraemia levels and weight gain contains a cluster of four other GBP genes that remain to be studied as candidate genes for this QTL.

[THE GENETIC DIVERSITY OF PIG BREEDS ON TEN PRODUCTION QUANTITATIVE TRAITS LOCI]

Genetic structure of 11 pig breeds and lines of different origins and production specialization for 10 quantitative traits loci: RYR1, GH, IGF2, CTSL, CTSS, CTSB, CTSK, MC4R, ESR1 & PRLR was analyzed. The results allow us to suppose that the uniqueness of genetic structure of each breed is largely determined by features of the allelic structure of a number of loci. Breeds of different production specialization differ significantly for such QTL, as RYR1, MC4R, ESR1 and IGF2. However, we have found that the loci CTSB, CTSL, CTSK, CTSS and PRLR in most of investigated breeds showed similar distribution of alleles, significant genetic differentiation between breeds and within breeds for these loci is missing. The exception is the Large Black breed for the CTSB locus and Ukrainian steppe spotted breed for the CTSL locus. Also, it can be assumed that allelic structure of some of investigated loci are more connected with the origins of pig breeds than to the selection on production traits. This refers, for example, to the growth hormone gene, the using of which in marker assisted selection could be effective only in some populations. The genetic distance between breeds and lines was defined and the dendrogram of the genetic relationship was built. Breeds are grouped in connection with the production specialization, as well as by their origins.

A polymorphism in the stearoyl-CoA desaturase gene promoter increases monounsaturated fatty acid content in dry-cured ham

Data on 125 dry-cured hams from purebred Duroc pigs were used to examine whether the favorable effect of the T allele in the promoter region of the stearoyl-CoA desaturase gene (AY487830:g.2228T>C) on monounsaturated fatty acid (MUFA) content in green ham is maintained after the curing process. It is shown that pigs carrying the T allele produced dry-cured hams with increased C16:1, C18:1n-9, C18:1n-7, and MUFA and decreased C18:0 and saturated fatty acid (SFA) content. The TT pigs had 1.39% more MUFA and 1.62% less SFA than the CC animals, while gilts had 0.74% more MUFA and 0.34% less SFA than barrows. The correlation between MUFA in green and dry-cured hams (n=53) was high (r=0.88), with TT pigs being more effective in retaining increased MUFA in green hams until the end of the curing period. It is concluded that increasing the presence of the T allele could have more impact than gender to produce hams with a high level of MUFA.

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.

A genome-wide association analysis for porcine serum lipid traits reveals the existence of age-specific genetic determinants

Background

The genetic determinism of blood lipid concentrations, the main risk factor for atherosclerosis, is practically unknown in species other than human and mouse. Even in model organisms, little is known about how the genetic determinants of lipid traits are modulated by age-specific factors. To gain new insights into this issue, we have carried out a genome-wide association study (GWAS) for cholesterol (CHOL), triglyceride (TRIG) and low (LDL) and high (HDL) density lipoprotein concentrations measured in Duroc pigs at two time points (45 and 190 days).

Results

Analysis of data with mixed-model methods (EMMAX, GEMMA, GenABEL) and PLINK showed a low positional concordance between trait-associated regions (TARs) for serum lipids at 45 and 190 days. Besides, the proportion of phenotypic variance explained by SNPs at these two time points was also substantially different. The four analyses consistently detected two regions on SSC3 (124 Mb, CHOL and LDL at 190 days) and SSC6 (135 Mb, CHOL and TRIG at 190 days) with highly significant effects on the porcine blood lipid profile. Moreover, we have found that SNP variation within SSC3, SSC6, SSC10, SSC13 and SSC16 TARs is associated with the expression of several genes mapping to other chromosomes and related to lipid metabolism.

Conclusions

Our data demonstrate that the effects of genomic determinants influencing lipid concentrations in pigs, as well as the amount of phenotypic variance they explain, are influenced by age-related factors.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-758) contains supplementary material, which is available to authorized users.

A functional variant in the stearoyl-CoA desaturase gene promoter enhances fatty acid desaturation in pork

There is growing public concern about reducing saturated fat intake. Stearoyl-CoA desaturase (SCD) is the lipogenic enzyme responsible for the biosynthesis of oleic acid (18 ∶ 1) by desaturating stearic acid (18 ∶ 0). Here we describe a total of 18 mutations in the promoter and 3' non-coding region of the pig SCD gene and provide evidence that allele T at AY487830:g.2228T>C in the promoter region enhances fat desaturation (the ratio 18 ∶ 1/18 ∶ 0 in muscle increases from 3.78 to 4.43 in opposite homozygotes) without affecting fat content (18 ∶ 0+18 ∶ 1, intramuscular fat content, and backfat thickness). No mutations that could affect the functionality of the protein were found in the coding region. First, we proved in a purebred Duroc line that the C-T-A haplotype of the 3 single nucleotide polymorphisms (SNPs) (g.2108C>T; g.2228T>C; g.2281A>G) of the promoter region was additively associated to enhanced 18 ∶ 1/18 ∶ 0 both in muscle and subcutaneous fat, but not in liver. We show that this association was consistent over a 10-year period of overlapping generations and, in line with these results, that the C-T-A haplotype displayed greater SCD mRNA expression in muscle. The effect of this haplotype was validated both internally, by comparing opposite homozygote siblings, and externally, by using experimental Duroc-based crossbreds. Second, the g.2281A>G and the g.2108C>T SNPs were excluded as causative mutations using new and previously published data, restricting the causality to g.2228T>C SNP, the last source of genetic variation within the haplotype. This mutation is positioned in the core sequence of several putative transcription factor binding sites, so that there are several plausible mechanisms by which allele T enhances 18 ∶ 1/18 ∶ 0 and, consequently, the proportion of monounsaturated to saturated fat.

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.

Association analysis with lipid traits of 2 candidate genes (LRP12 and TRIB1) mapping to a SSC4 QTL for serum triglyceride concentration in pigs

The performance of a genome scan for serum lipid traits at 45 and 190 d in 5 half-sib families of Duroc pigs allowed us to detect several pig chromosomal regions with significant effects on these phenotypes. In the current work, we aimed to refine the position of 1 chromosome 4 (SSC4) genome-wide significant QTL for serum triglyceride concentration at 190 d. Genotyping of 4 additional microsatellites allowed reduction of the 90% confidence interval of this QTL to the genomic interval between markers SW2409 and SW839. Sequencing experiments were performed to characterize the variability of 2 lipid-related genes, the lipoprotein receptor-related protein 12 (LRP12) and tribbles homolog 1 (TRIB1) loci, that map to this region. In this way, 2 (c.771A > G and c.1101A > G) and 1 (c.*156_157del) polymorphisms were identified at the LRP12 coding region and TRIB1 3' untranslated region, respectively. Association analyses between LRP12 and TRIB1 genotypes did not reveal any significant effect on serum lipid concentrations, suggesting that variation of these two loci does not explain the segregation of the SSC4 QTL. However, highly significant associations were observed for gluteus medius saturated fatty acid content (LRP12 c.1101A > G, P = 0.0006; TRIB1 c.*156_157del, P = 0.0003). In the light of these and other findings, the potential involvement of LRP12 and TRIB1 in muscle lipid metabolism deserves to be further explored.

Genomic architecture of heritability and genetic correlations for intramuscular and back fat contents in Duroc pigs

Genetic parameters such as heritability and correlations of fat traits in a Duroc population were dissected using molecular markers. The heritabilities of intramuscular fat in 2 muscles, the gluteus medius and LM, and back fat were 0.54, 0.48, and 0.23, respectively. The genetic correlations were well estimated with standardized SNP effects, being 0.65 between intramuscular fat traits and ∼0.37 between any intramuscular fat trait and back fat. Genetic correlations were overestimated when ignoring molecular information. Twelve chromosomes showed additive genetic variance for intramuscular fat compared with 8 for back fat. Population structure was accommodated using 4 different models. The number of significant, P < 5 × 10(-5) (suggestive, P < 2 × 10(-3)), SNP varied across models and ranged from 0 to 4 (2 to 261) for intramuscular fat in the gluteus medius, from 0 to 57 (9 to 564) for intramuscular fat in the LM, and from 3 to 4 (22 to 168) for back fat. Several SNP showed significant deviations from an additive mode of action. Only 2 SNP significantly affected 2 traits simultaneously.

Mapping and tissue mRNA expression analysis of the pig solute carrier 27A (SLC27A) multigene family

Solute-carrier family 27A molecules are integral transmembrane proteins that play a fundamental role in the uptake of long-chain fatty acids into mammalian cells. Our goal was to characterize this multigene family in pigs. Chromosomal location of the six porcine SLC27A genes was determined by radiation hybrid mapping and indicated that the six genes map to six different chromosomal locations. Moreover, we analyzed SLC27A mRNA expression in six pig tissues by quantitative RT-PCR. While SLC27A1, SLC27A3 and SLC27A4 were expressed in most, if not all, analyzed tissues, SLC27A2, SLC27A5 and SLC27A6 were predominantly expressed in the liver. In general, pig and human SLC27A mRNA expression profiles were remarkably concordant, although important differences were observed for SLC27A1 and SLC27A6 mRNAs. Discrepancies between mRNA expression profiles have been observed even in closely related primate species, and they might reflect the acquisition of regulatory changes promoting evolutionary adaptation.

Segregation of regulatory polymorphisms with effects on the gluteus medius transcriptome in a purebred pig population

Background

The main goal of the present study was to analyse the genetic architecture of mRNA expression in muscle, a tissue with an outmost economic importance for pig breeders. Previous studies have used F₂ crosses to detect porcine expression QTL (eQTL), so they contributed with data that mostly represents the between-breed component of eQTL variation. Herewith, we have analysed eQTL segregation in an outbred Duroc population using two groups of animals with divergent fatness profiles. This approach is particularly suitable to analyse the within-breed component of eQTL variation, with a special emphasis on loci involved in lipid metabolism.

Methodology/Principal Findings

GeneChip Porcine Genome arrays (Affymetrix) were used to determine the mRNA expression levels of gluteus medius samples from 105 Duroc barrows. A whole-genome eQTL scan was carried out with a panel of 116 microsatellites. Results allowed us to detect 613 genome-wide significant eQTL unevenly distributed across the pig genome. A clear predominance of trans- over cis-eQTL, was observed. Moreover, 11 trans-regulatory hotspots affecting the expression levels of four to 16 genes were identified. A Gene Ontology study showed that regulatory polymorphisms affected the expression of muscle development and lipid metabolism genes. A number of positional concordances between eQTL and lipid trait QTL were also found, whereas limited evidence of a linear relationship between muscle fat deposition and mRNA levels of eQTL regulated genes was obtained.

Conclusions/Significance

Our data provide substantial evidence that there is a remarkable amount of within-breed genetic variation affecting muscle mRNA expression. Most of this variation acts in trans and influences biological processes related with muscle development, lipid deposition and energy balance. The identification of the underlying causal mutations and the ascertainment of their effects on phenotypes would allow gaining a fundamental perspective about how complex traits are built at the molecular level.

Joint effects of porcine leptin and leptin receptor polymorphisms on productivity and quality traits

Leptin signalling plays a fundamental role in growth, fatness and body composition. The aim of this study was to investigate the porcine LEP gene sequence in an Iberian × Landrace experimental cross to identify polymorphisms associated with productivity and quality traits. Because of the documented effects on these traits of the LEPR c.1987C>T polymorphism, the LEP and LEPR c.1987C>T polymorphisms and their interactions have been jointly investigated. The LEP gene sequencing has allowed the identification of 39 polymorphisms, eight of which are novel. Three intronic SNPs, LEP g.1382C>T, LEP g.1387C>T and LEP g.1723A>G, have been genotyped, and association analyses have been carried out. Analyses of LEP g.1387C>T, fully linked to LEP g.1382C>T, have revealed additive effects on live and carcass weights and dominant effects on several backfat thickness measurements. Novel effects of both LEP and LEPR polymorphisms on fatty acid composition in subcutaneous fat have been detected, probably mediated through the effects on fatness. The results reported here suggest that the T alleles of both LEP g.1387C>T and LEPR c.1987C>T, which are fixed in the Iberian pigs, would lead to an increase in growth, fatness and saturated fatty acid content in fat, which could be explained by an increased feed intake.

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.

Gene encoding the prothoracicotropic hormone of a moth is expressed in the brain and gut

The molts of lepidopteran insects are typically controlled by the brain-derived prothoracicotropic hormone (PTTH) that stimulates ecdysteroidogenesis in the prothoracic glands (PGs). We report here that the larvae and pupae of the moth Sesamia nonagrioides can molt without brain (PGs must be present), suggesting that there might be a secondary source of PTTH. We addressed this issue by characterizing spatial and temporal expression patterns of the PTTH gene. To this end we identified a major part of the corresponding cDNA. Protein deduced from this cDNA fragment consisted of 128 amino acids and showed 48-85% homology with the matching regions of PTTHs known from other Lepidoptera. Quantification of PTTH expression in major body organs of the last instar larvae revealed high expression in the brain (fading in post-feeding larvae) and considerable expression in the gut (with a maximum in post-feeding larvae). The content of PTTH message in the gut was enhanced after decapitation. It is concluded that the molts of S. nonagrioides larvae are driven by PTTH gene expression in the gut.

Muscle transcriptomic profiles in pigs with divergent phenotypes for fatness traits

BACKGROUND

Selection for increasing intramuscular fat content would definitively improve the palatability and juiciness of pig meat as well as the sensorial and organoleptic properties of cured products. However, evidences obtained in human and model organisms suggest that high levels of intramuscular fat might alter muscle lipid and carbohydrate metabolism. We have analysed this issue by determining the transcriptomic profiles of Duroc pigs with divergent phenotypes for 13 fatness traits. The strong aptitude of Duroc pigs to have high levels of intramuscular fat makes them a valuable model to analyse the mechanisms that regulate muscle lipid metabolism, an issue with evident implications in the elucidation of the genetic basis of human metabolic diseases such as obesity and insulin resistance.

RESULTS

Muscle gene expression profiles of 68 Duroc pigs belonging to two groups (HIGH and LOW) with extreme phenotypes for lipid deposition and composition traits have been analysed. Microarray and quantitative PCR analysis showed that genes related to fatty acid uptake, lipogenesis and triacylglycerol synthesis were upregulated in the muscle tissue of HIGH pigs, which are fatter and have higher amounts of intramuscular fat than their LOW counterparts. Paradoxically, lipolytic genes also showed increased mRNA levels in the HIGH group suggesting the existence of a cycle where triacylglycerols are continuously synthesized and degraded. Several genes related to the insulin-signalling pathway, that is usually impaired in obese humans, were also upregulated. Finally, genes related to antigen-processing and presentation were downregulated in the HIGH group.

CONCLUSION

Our data suggest that selection for increasing intramuscular fat content in pigs would lead to a shift but not a disruption of the metabolic homeostasis of muscle cells. Future studies on the post-translational changes affecting protein activity or expression as well as information about protein location within the cell would be needed to to elucidate the effects of lipid deposition on muscle metabolism in pigs.

Evaluation of FABP2 as candidate gene for a fatty acid composition QTL in porcine chromosome 8

The objective of this work was to analyse the porcine Fatty acid binding protein 2, intestinal (FABP2) gene as a candidate gene for a fatty acid composition quantitative trait loci (QTL) previously described on porcine chromosome 8 in an Iberian by Landrace F(2) cross (IBMAP). Re-sequencing of the porcine FABP2 gene in three Iberian and eight Landrace parental animals resulted in the identification of three single-nucleotide polymorphisms, all of them localized in intron 1. The polymorphism FABP2:g.412T>C, localized in intron 1, and two additional microsatellites were genotyped in the IBMAP population in order to perform an association test of the FABP2 gene and to better define the QTL position previously described. Association analyses of the FABP2:g.412T>C with the fatty acid composition traits were not significant in simple association and marker-assisted association tests, suggesting that the FABP2 region sequenced is not responsible for the QTL. However, the addition of three new markers to the pedigree allowed us to define the S0144-SW61 marker interval as the most likely QTL position, facilitating the future study of other candidate genes for this QTL.

Mammary gland differentiation inversely correlates with GDF-8 expression

GDF-8 is recognised as an inhibitor of muscle cell growth and differentiation. Although initially thought to be restricted to muscle cells it is now accepted that GDF-8 expression has a broader tissue distribution. We demonstrate GDF-8 expression in the mouse mammary gland, which is predominantly associated with epithelial cells and displays an inverse correlation to the differentiated state of the gland. Specifically, the highest GDF-8 mRNA levels correlate with periods of maximal ductal growth, diminish as pregnancy progressed and are down-regulated to minimal levels by the onset of lactation as the epithelium differentiates. A similar profile is observed for both GDF-8 protein processing and reflects Smad2/3 phosphorylation profile. However, in contrast to muscle cells, GDF-8 neither reduces proliferation nor induces p21 expression levels in mammary epithelial cells. These data implicate a role for GDF-8 in mammary epithelial cell differentiation and demonstrate that GDF-8 has cell-type specific activities.

Mapping of quantitative trait loci for cholesterol, LDL, HDL, and triglyceride serum concentrations in pigs

The fine mapping of polymorphisms influencing cholesterol (CT), triglyceride (TG), and lipoprotein serum levels in human and mouse has provided a wealth of knowledge about the complex genetic architecture of these traits. The extension of these genetic analyses to pigs would be of utmost importance since they constitute a valuable biological and clinical model for the study of coronary artery disease and myocardial infarction. In the present work, we performed a whole genome scan for serum lipid traits in a half-sib Duroc pig population of 350 individuals. Phenotypic registers included total CT, TG, and low (LDL)- and high (HDL)-density lipoprotein serum concentrations at 45 and 190 days of age. This approach allowed us to identify two genomewide significant quantitative trait loci (QTL) for HDL-to-LDL ratio at 45 days (SSC6, 84 cM) and for TG at 190 days (SSC4, 23 cM) as well as a number of chromosomewide significant QTL. The comparison of QTL locations at 45 and 190 days revealed a notable lack of concordance at these two time points, suggesting that the effects of these QTL are age specific. Moreover, we have observed a considerable level of correspondence among the locations of the most significant porcine lipid QTL and those identified in humans. This finding might suggest that, in mammals, diverse polymorphisms located in a common set of genes are involved in the genetic variation of serum lipid levels.

Duplication of Stat5-binding sites within the β-lactoglobulin promoter compromises transcription in vivo

Achieving high and reliable levels of expression is essential for biotechnological application of transgene technology. In an attempt to maximise transgene expression, we analysed the effect of adding extra Stat5-binding sites, while maintaining their spatial arrangement, to the beta-lactoglobulin promoter. A 195 bp-fragment comprising the three Stat5-binding sites in the beta-lactoglobulin promoter was duplicated directly upstream, 88 bp or 177 bp upstream from the original site. We assessed the effect that addition of extra Stat5-binding sites had in transcriptional activation and chromatin structuring of the promoter, both in vitro in HC11 cells and in vivo in transgenic animals. Duplication of Stat5-binding sites in the beta-lactoglobulin promoter did favour chromatin remodelling in vivo but unexpectedly did not enhance transcriptional activation of the transgene, despite the higher basal and induced expression levels achieved in vitro in HC11 cells. Furthermore, the well-described positive relationship between beta-lactoglobulin transgene copy number and mammary expression was lost in the new lines. These results indicate that chromatin remodelling upon Stat5 binding is not sufficient to activate transgene expression, denote that strategies to achieve transgenic transcriptional activation are not obvious.

Transgene Methylation in Mice Reflects Copy Number But Not Expression Level

In mammals, CpG methylation is one of the mechanisms of epigenetic control over the linear sequence of bases of deoxyribonucleic acid (DNA); about 70% of CpG dinucleotides are methylated. The actual signal that triggers DNA methylation is not known, although repetitive DNA has been shown to be an attractive template for DNA methylases. To address methylation events associated with transgenic copy number, we have analyzed transgenes that are actively transcribed in a tissue-specific manner. We have compared gross transgene methylation by restriction-enzyme digestion in expressing and nonexpressing tissues. The observed pattern suggests that the DNA methylation machinery can recognize repeated genomic sequences independently of their transcriptional activity.

Characterization of genetic polymorphism in the goat beta-lactoglobulin gene

Two new variants have been detected and characterized for the goat beta-lactoglobulin gene at the cDNA level and confirmed at the genomic level. The two polymorphisms are located on exon 7 of the gene. One of the polymorphic sites is produced by a single nucleotide substitution in position +4601, allowing a polymerase chain reaction--restriction fragment length polymorphism (PCR-RFLP) genotyping procedure to be developed using SacII restriction enzyme. The other polymorphic position contains a 10 bp long insertion at position +4641 that can be detected by capillary electrophoresis of the PCR product amplified with a fluorescent primer. The association of these two polymorphisms was also investigated, resulting in the description of two new alleles. Both of these contained the point mutation at the SacII site, with or without the 10 bp insertion at position +4641. The distribution of these new polymorphisms was studied in a population of males of four different goat breeds. The gene frequencies for these variants were similar in Spanish and French breeds.

Isolation, sequencing and relative quantitation by fluorescent-ratio PCR of feline beta-lactoglobulin I, II, and III cDNAs

Beta-lactoglobulin (betaLG) is a whey protein found in the milk of most mammals, except those of humans and rodents. In ruminants, only one type of betaLG is expressed in milk, although the presence of pseudogenes has been reported. The milk of other mammals (dog, horse, dolphins) contains two types of betaLG: type I is related to ruminant betaLG, while type II has a strong similarity to the ruminant pseudogenes. The presence of three types of betaLG has been described only in cat milk. For the first time, we have cloned, sequenced, and characterized the three types of feline betaLG cDNAs from a mammary gland sample of a lactating cat. Since no specific probe could be easily used to differentiate them by Northern Blotting, we developed a new technique named Fluorescent-Ratio PCR (FR-PCR) in order to assess their level of expression during lactation. The relative amounts of each feline betaLG cDNA was quantified in two cats by capillary electrophoresis of the restricted RT-PCR product, labeled with fluorescent primers. We observed the same percentage of expression of the three betaLG genes in two cats. The differences of expression could be due to changes in the promoter and 3' non-coding region affecting the level of transcription and the mRNA stability. The FR-PCR technique was shown to be reproducible and accurate even for small percentage differences, being very useful when a small amount of sample is available.

Chromatin structures of goat and sheep beta-lactoglobulin gene differ

Different levels of the major milk protein beta-lactoglobulin are found in evolutionarily related ruminant species: with sheep milk containing as much as three times the concentration in goat milk. In an attempt to understand why these differences exist, we have characterised, using DNaseI as a probe of structure, the chromatin surrounding the goat beta-lactoglobulin promoter and compared it to that of the sheep homologue. The goat gene displays a mammary-specific chromatin pattern, which is reformed on expressing goat beta-lactoglobulin transgenes. This implies that this chromatin structure is sequence dependent and suggests that it plays a role in regulating beta-lactoglobulin gene expression. This pattern differs from that seen on the ovine beta-lactoglobulin gene in lactating sheep mammary chromatin. Thus, even between highly related species, the transcriptional mechanisms regulating activity of a gene can differ.