Ramona N. Pena

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.

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.

Investigating reference genes for quantitative real-time PCR analysis across four chicken tissues

Accurate normalization of data is required to correct for different efficiencies and errors during the processing of samples in reverse transcription PCR analysis. The chicken is one of the main livestock species and its genome was one of the first reported and used in large scale transcriptomic analysis. Despite this, the chicken has not been investigated regarding the identification of reference genes suitable for the quantitative PCR analysis of growth and fattening genes. In this study, five candidate reference genes (B2M, RPL32, SDHA, TBP and YWHAZ) were evaluated to determine the most stable internal reference for quantitative PCR normalization in the two main commercial muscles (pectoralis major (breast) and biceps femoris (thigh)), liver and abdominal fat. Four statistical methods (geNorm, NormFinder, CV and BestKeeper) were used in the evaluation of the most suitable combination of reference genes. Additionally, a comprehensive ranking was established with the RefFinder tool. This analysis identified YWHAZ and TBP as the recommended combination for the analysis of biceps femoris and liver, YWHAZ and RPL32 for pectoralis major and RPL32 and B2M for abdominal fat and across-tissue studies. The final ranking for each tool changed slightly but overall the results, and most particularly the ability to discard the least robust candidates, were consistent between tools. The selection and number of reference genes were validated using SCD, a target gene related to fat metabolism. Overall, the results can be directly used to quantitate target gene expression in different tissues or in validation studies from larger transcriptomic experiments.

Selection of Internal Control Genes for Real-Time Quantitative PCR in Ovary and Uterus of Sows across Pregnancy

Background Reproductive traits play a key role in pig production in order to reduce costs and increase economic returns. Among others, gene expression analyses represent a useful approach to study genetic mechanisms underlying reproductive traits in pigs. The application of reverse-transcription quantitative PCR requires the selection of appropriate reference genes, whose expression levels should not be affected by the experimental conditions, especially when comparing gene expression across different physiological stages. Results The gene expression stability of ten potential reference genes was studied by three different methods (geNorm, NormFinder and BestKeeper) in ovary and uterus collected at five different physiological time points (heat, and 15, 30, 45 and 60 days of pregnancy). Although final ranking differed, the three algorithms gave very similar results. Thus, the most stable genes across time were TBP and UBC in uterus and TBP and HPRT1 in ovary, while HMBS and ACTB showed the less stable expression in uterus and ovary, respectively. When studied as a systematic effect, the reproductive stage did not significantly affect the expression of the candidate reference genes except at 30d and 60d of pregnancy, when a general drop in expression was observed in ovary. Conclusions Based in our results, we propose the use of TBP, UBC and SDHA in uterus and TBP, GNB2L1 and HPRT1 in ovary for normalization of longitudinal expression studies using quantitative PCR in sows.

Genetic Marker Discovery in Complex Traits: A Field Example on Fat Content and Composition in Pigs

Among the large number of attributes that define pork quality, fat content and composition have attracted the attention of breeders in the recent years due to their interaction with human health and technological and sensorial properties of meat. In livestock species, fat accumulates in different depots following a temporal pattern that is also recognized in humans. Intramuscular fat deposition rate and fatty acid composition change with life. Despite indication that it might be possible to select for intramuscular fat without affecting other fat depots, to date only one depot-specific genetic marker (PCK1 c.2456C>A) has been reported. In contrast, identification of polymorphisms related to fat composition has been more successful. For instance, our group has described a variant in the stearoyl-coA desaturase (SCD) gene that improves the desaturation index of fat without affecting overall fatness or growth. Identification of mutations in candidate genes can be a tedious and costly process. Genome-wide association studies can help in narrowing down the number of candidate genes by highlighting those which contribute most to the genetic variation of the trait. Results from our group and others indicate that fat content and composition are highly polygenic and that very few genes explain more than 5% of the variance of the trait. Moreover, as the complexity of the genome emerges, the role of non-coding genes and regulatory elements cannot be disregarded. Prediction of breeding values from genomic data is discussed in comparison with conventional best linear predictors of breeding values. An example based on real data is given, and the implications in phenotype prediction are discussed in detail. The benefits and limitations of using large SNP sets versus a few very informative markers as predictors of genetic merit of breeding candidates are evaluated using field data as an example.

A WUR SNP is associated with European Porcine Reproductive and Respiratory Virus Syndrome resistance and growth performance in pigs

Porcine reproductive and respiratory syndrome (PRRS) causes decreased reproductive performance and respiratory problems in pigs. The goals of the current study were 1) to examine whether individual variation applies to infection with PRRSV European strains and 2) to investigate the association of a single nucleotide polymorphism (SNP) WUR10000125 (WUR) at the interferon-inducible guanylate-binding protein 1 gene (GBP1) with average daily gain (ADG) in PRRSV infected and uninfected pigs. The experimental procedure consisted of two trials in which pigs from negative PRRSV farms were infected with a wild-type (n=80) or vaccinated with an attenuated European PRRS virus strain (n=40) and then monitored after infection or vaccination. Viral load and ADG were determined for each pig. In a third trial, the ADG for PRRSV-free pigs was monitored. All pigs were genotyped for WUR at the GBP1 gene (AA and AG genotype were defined). Results indicated that there was individual variation in the viral load from pigs challenged with a wild-type or low virulent European PRRSV strain. Secondly, our data showed that WUR SNP was associated to ADG in vaccinated pigs. Thus, ADG in AG pigs was significantly higher than in AA ones after vaccinating with an attenuated PRRSV strain. However, the reverse happened in a PRRSV-free environment where the AA pigs were those that grew faster. Based on these results, there is a scope for selecting pigs according to their responses to PRRS virus infection with European strains and that WUR SNP may play a role in causing PRRSV resistance.

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.

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 (pH24), 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 pH24 SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

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.

A polymorphism in the fatty acid desaturase-2 gene is associated with the arachidonic acid metabolism in pigs

Arachidonic acid (C20:4) is related to a wide range of biological effects including lipid homeostasis. The fatty acid desaturase-2 (FADS2) gene encodes for the delta-6-desaturase, which is involved in the biosynthesis of C20:4 from linoleic acid (C18:2). The purpose of this study was to characterise mutations in the promoter of the porcine FADS2, evaluating in particular the effect of one haplotype tagging polymorphism (rs321384923A > G) on the biosynthesis pathway of C20:4. A total of 1,192 Duroc barrows with records on fatty acid composition in muscle and subcutaneous fat were genotyped. Pigs carrying the A allele showed, irrespective of fat content, both enhanced FADS2 expression and higher C20:4 in muscle and exhibited increased ratios of C20:4 to C18:2 and of C20:4 to eicosadienoic acid (C20:2) in both muscle and adipose tissue. Despite the inverse relationship observed between C20:4 and fat content, the rs321384923 polymorphism had no impact on lean weight. It is concluded that the haplotype encompassing the rs321384923 polymorphism at the porcine FADS2 affects the n-6 fatty acid profile by specifically modifying the desaturation efficiency of C18:2 to C20:4 rather than by concomitant variations in C18:2 following changes in fat content.