A QTL on pig chromosome 4 affects fatty acid metabolism: evidence from an Iberian by Landrace intercross

Genome-wide association analysis explores the genetic loci of amino acid content in duck's breast muscle

BACKGROUND

Amino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.

RESULT

In the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55-56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2-21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.

CONCLUSION

In this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Regulation of Swine Growth by Backfat Tissue during Growing and Finishing Stages

Simple Summary

Swine have a large influence on livestock animals. In particular, Korean native pigs (KNPs) have unique traits in their body composition including lipids and proteins. In this study, we performed RNA-sequencing analysis to identify porcine transcriptomic changes during growing and finishing stages in the backfat tissue of KNP and Yorkshire pig crossbreeds. Enrichment analysis revealed that differentially expressed genes (DEGs) were significantly influenced by lipid metabolism and hormonal changes. Network analysis showed that the LEP and ACTC1 genes were connected with significant terminologies which resulted from up- and down-regulated DEGs. The results of our analysis indicate that backfat tissue could regulate swine biology during stages of growth. Consequently, our analysis provided comprehensive understanding for transcriptomic changes during growth periods.

Abstract

Recently, interest in the function of pig backfat (BF) has increased in the field of livestock animals, and many transcriptome-based studies using commercial pig breeds have been conducted. However, there is a lack of comprehensive studies regarding the biological mechanisms of Korean native pigs (KNPs) and Yorkshire pig crossbreeds. In this study, therefore, BF samples of F1 crossbreeds of KNPs and Yorkshire pigs were investigated to identify differentially expressed genes (DEGs) and their related terms using RNA-sequencing analysis. DEG analysis identified 611 DEGs, of which 182 were up-regulated and 429 were down-regulated. Lipid metabolism was identified in the up-regulated genes, whereas growth and maturation-related terminologies were identified in the down-regulated genes. LEP and ACTC1 were identified as highly connected core genes during functional gene network analysis. Fat tissue was observed to affect lipid metabolism and organ development due to hormonal changes driven by transcriptional alteration. This study provides a comprehensive understanding of BF contribution to crossbreeds of KNPs and Yorkshire pigs during growth periods.

Expression analysis of porcine miR-33a/b in liver, adipose tissue and muscle and its potential role in fatty acid metabolism

mir-33a and mir-33b are co-transcribed with the SREBF2 and SREBF1 transcription factors, respectively. The main role of SREBF1 is the regulation of genes involved in fatty acid metabolism, while SREBF2 regulates genes participating in cholesterol biosynthesis and uptake. Our objective was to study the expression of both miR-33a and miR-33b, together with their host SREBF genes, in liver, adipose tissue and muscle to better understand the role of miR-33a/b in the lipid metabolism of pigs. In our study, the expression of miR-33a, miR-33b and SREBF2 in liver, adipose tissue, and muscle was studied in 42 BC1_LD (25% Iberian x 75% Landrace backcross) pigs by RT-qPCR. In addition, the expression of in-silico predicted target genes and fatty acid composition traits were correlated with the miR-33a/b expression. We observed different tissue expression patterns for both miRNAs. In adipose tissue and muscle a high correlation between miR-33a and miR-33b expression was found, whereas a lower correlation was observed in liver. The expression analysis of in-silico predicted target-lipid related genes showed negative correlations between miR-33b and CPT1A expression in liver. Conversely, positive correlations between miR-33a and PPARGC1A and USF1 gene expression in liver were observed. Lastly, positive and negative correlations between miR-33a/b expression and saturated fatty acid (SFA) and polyunsaturated fatty acid (PUFA) content, respectively, were identified. Overall, our results suggested that both miRNAs are differentially regulated and have distinct functions in liver, in contrast to muscle and adipose tissue. Furthermore, the correlations between miR-33a/b expression both with the expression of in-silico predicted target-lipid related genes and with fatty acid composition, opens new avenues to explore the role of miR33a/b in the regulation of lipid metabolism.

Identification of strong candidate genes for backfat and intramuscular fatty acid composition in three crosses based on the Iberian pig

Meat quality has an important genetic component and can be modified by the fatty acid (FA) composition and the amount of fat contained in adipose tissue and muscle. The present study aimed to find genomic regions associated with the FA composition in backfat and muscle (longissimus dorsi) in 439 pigs with three different genetic backgrounds but having the Iberian breed in common. Genome-wide association studies (GWAS) were performed between 38,424 single-nucleotide polymorphisms (SNPs) covering the pig genome and 60 phenotypic traits related to backfat and muscle FA composition. Nine significant associated regions were found in backfat on the Sus scrofa chromosomes (SSC): SSC1, SSC2, SSC4, SSC6, SSC8, SSC10, SSC12, and SSC16. For the intramuscular fat, six significant associated regions were identified on SSC4, SSC13, SSC14, and SSC17. A total of 52 candidate genes were proposed to explain the variation in backfat and muscle FA composition traits. GWAS were also reanalysed including SNPs on five candidate genes (ELOVL6, ELOVL7, FADS2, FASN, and SCD). Regions and molecular markers described in our study may be useful for meat quality selection of commercial pig breeds, although several polymorphisms were breed-specific, and further analysis would be needed to evaluate possible causal mutations.

Analysis of porcine IGF2 gene expression in adipose tissue and its effect on fatty acid composition

IGF2:g.3072G>A polymorphism has been described as the causal mutation of a maternally imprinted QTL for muscle growth and fat deposition in pigs. The objective of the current work was to study the association between the IGF2:g.3072G>A polymorphism and the IGF2 gene expression and its effect on fatty acid composition in adipose tissue in different pig genetic backgrounds. A cis-eQTL region associated with the IGF2 mRNA expression in adipose tissue was identified in an eGWAS with 355 animals. The IGF2 gene was located in this genomic interval and IGF2g.3072G>A was the most significant SNP, explaining a 25% of the gene expression variance. Significant associations between IGF2:g.3072G>A polymorphism and oleic (C18:1(n-9); p-value = 4.18x10^-07), hexadecanoic (C16:1(n-9); p-value = 4.04x10^-07), linoleic (C18:2(n-6); p-value = 6.44x10^-09), α-linoleic (C18:3(n-3); p-value = 3.30x10^-06), arachidonic (C20:4(n-6); p-value = 9.82x10^-08) FAs and the MUFA/PUFA ratio (p-value = 2.51x10^-9) measured in backfat were identified. Animals carrying the A allele showed an increase in IGF2 gene expression and higher PUFA and lower MUFA content. However, in additional studies was observed that there could be other proximal genetic variants affecting FA composition in adipose tissue. Finally, no differences in the IGF2 gene expression in adipose tissue were found between heterozygous animals classified according to the IGF2:g.3072G>A allele inherited from the father (A^PG^M or A^MG^P). However, pyrosequencing analysis revealed that there is imprinting of the IGF2 gene in muscle and adipose tissues, with stronger differences among the paternally and maternally inherited alleles in muscle. Our results suggested that IGF2:g.3072G>A polymorphism plays an important role in the regulation of IGF2 gene expression and can be involved in the fatty acid composition in adipose tissue. In both cases, further studies are still needed to deepen the mechanism of regulation of IGF2 gene expression in adipose tissue and the IGF2 role in FA composition.

Indel detection from Whole Genome Sequencing data and association with lipid metabolism in pigs

The selection in commercial swine breeds for meat-production efficiency has been increasing among the past decades, reducing the intramuscular fat content, which has changed the sensorial and technological properties of pork. Through processes of natural adaptation and selective breeding, the accumulation of mutations has driven the genetic divergence between pig breeds. The most common and well-studied mutations are single-nucleotide polymorphisms (SNPs). However, insertions and deletions (indels) usually represents a fifth part of the detected mutations and should also be considered for animal breeding. In the present study, three different programs (Dindel, SAMtools mpileup, and GATK) were used to detect indels from Whole Genome Sequencing data of Iberian boars and Landrace sows. A total of 1,928,746 indels were found in common with the three programs. The VEP tool predicted that 1,289 indels may have a high impact on protein sequence and function. Ten indels inside genes related with lipid metabolism were genotyped in pigs from three different backcrosses with Iberian origin, obtaining different allelic frequencies on each backcross. Genome-Wide Association Studies performed in the Longissimus dorsi muscle found an association between an indel located in the C1q and TNF related 12 (C1QTNF12) gene and the amount of eicosadienoic acid (C20:2(n-6)).

Variation in the FABP4 gene affects carcass and growth traits in sheep

Fatty acid-binding protein 4 (FABP4) participates in fatty-acid transportation and variation in the gene has been reported to affect fat deposition in mammals. To assess the effects of variation in FABP4 on carcass and growth traits in sheep, two regions of the gene were investigated in 749 New Zealand Romney lambs. Five variants (A₁ - E₁) were found in region-1 (exon 2 - intron 2) and three variants (A₂ - C₂) were found in region-2 (exon 3 - intron 3). A₁ was associated with a decrease in leg meat yield, loin meat yield and total meat yield, whereas A₂ was associated with a decrease in weaning-weight and pre-weaning growth-rate. Haplotype A₁-A₂ was found to be associated with a decrease in birth-weight, pre-weaning growth-rate, hot carcass weight, loin meat yield, shoulder meat yield and total meat yield, while haplotype A₁-B₂ was associated with increased fat depth at the 12th rib (V-GR). Taken together, this supports the contention that variation in FABP4 affects growth and meat production in sheep.

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

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

Expression analysis of candidate genes for fatty acid composition in adipose tissue and identification of regulatory regions

The aim of this work was to study the genetic basis of the backfat expression of lipid-related genes associated with meat quality traits in pigs. We performed a genome-wide association study with the backfat gene expression measured in 44 genes by qPCR and the PorcineSNP60 BeadChip genotypes in 115 Iberian x Landrace backcross animals. A total of 193 expression-associated SNPs located in 19 chromosomal regions were associated with expression levels of ACSM5, ELOVL6, FABP4, FADS2, and SLC27A4 genes. Three expression quantitative trail loci (eQTLs) corresponding to ACSM5, FABP4, and FADS2 were classified as cis-acting eQTLs, whereas the remaining 16 eQTLs have trans-regulatory effects. Remarkably, a SNP in the ACSM5 promoter region and a SNP in the 3′UTR region of FABP4 were the most associated polymorphisms with the ACSM5 and FABP4 expression levels, respectively. Moreover, relevant lipid-related genes mapped in the trans-eQTLs regions associated with the ACSM5, FABP4, FADS2, and SLC27A4 genes. Interestingly, a trans-eQTL hotspot on SSC13 regulating the gene expression of ELOVL6, ELOLV5, and SCD, three important genes implicated in the elongation and desaturation of fatty acids, was identified. These findings provide new data to further understand the functional regulatory mechanisms implicated in the variation of fatty acid composition in pigs.

Association of genetic variants and expression levels of porcine FABP4 and FABP5 genes

The FABP4 and FABP5 genes, coding for fatty acid transport proteins, have long been studied as positional candidate genes for SSC4 QTL affecting fat deposition and composition traits in pigs. Polymorphisms in these genes, FABP4:g.2634_2635insC and FABP5:g.3000T>G, have previously been associated with fatness traits in an Iberian by Landrace cross (IBMAP). The aim of the present work was to evaluate the functional implication of these genetic variants. For this purpose, FABP4 and FABP5 mRNA expression levels in 114 BC1_LD animals (25% Iberian × 75% Landrace) were analyzed using real-time quantitative PCR in backfat and muscle. FABP4 gene expression in backfat, but not in muscle, was associated with FABP4:g.2634_2635insC. In contrast, FABP5:g.3000T>G was not associated with gene expression levels. An expression-based genome-wide association study highlighted the FABP4:g.2634_2635insC polymorphism as the polymorphism most associated with FABP4 gene expression in backfat. Furthermore, other genomic regions associated in trans with the mRNA expression of FABP4 in backfat and FABP5 in muscle were also identified. Finally, two putative transcription binding sites for PPARG and NR4A2 may be affected by the FABP4:g.2634_2635insC polymorphism, modifying FABP4 gene expression. Our results reinforce FABP4 as a candidate gene for fatness traits on SSC4.

Genome-wide association study confirm major QTL for backfat fatty acid composition on SSC14 in Duroc pigs

Background

Fatty acid composition contributes importantly to meat quality and is essential to the nutritional value of the meat. Identification of genetic factors underlying levels of fatty acids can be used to breed for pigs with healthier meat. The aim of this study was to conduct genome-wide association studies (GWAS) to identify QTL regions affecting fatty acid composition in backfat from the pig breeds Duroc and Landrace.

Results

Using data from the Axiom porcine 660 K array, we performed GWAS on 454 Duroc and 659 Landrace boars for fatty acid phenotypes measured by near-infrared spectroscopy (NIRS) technology (C16:0, C16:1n-7, C18:0, C18:1n-9, C18:2n-6, C18:3n-3, total saturated fatty acids, monounsaturated fatty acids and polyunsaturated fatty acids). Two QTL regions on SSC4 and SSC14 were identified in Duroc for the de novo synthesized fatty acids traits, whereas one QTL on SSC8 was detected in Landrace for C16:1n-7. The QTL region on SSC14 has been reported in previous studies and a putative causative mutation has been suggested in the promoter region of the SCD gene. Whole genome re-sequencing data was used for genotype imputation and to fine map the SSC14 QTL region in Norwegian Duroc. This effort confirms the location of the QTL on this chromosome as well as suggesting other putative candidate genes in the region. The most significant single nucleotide polymorphisms (SNPs) located on SSC14 explain between 55 and 76% of the genetic variance and between 27 and 54% of the phenotypic variance for the de novo synthesized fatty acid traits in Norwegian Duroc. For the QTL region on SSC8 in Landrace, the most significant SNP explained 19% of the genetic variance and 5% of the phenotypic variance for C16:1n-7.

Conclusions

This study confirms a major QTL affecting fatty acid composition on SSC14 in Duroc, which can be used in genetic selection to increase the level of fatty acid desaturation. The SSC14 QTL was not segregating in the Landrace population, but another QTL on SSC8 affecting C16:1n-7 was identified and might be used to increase the level of desaturation in meat products from this breed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3752-0) contains supplementary material, which is available to authorized users.

A global analysis of CNVs in swine using whole genome sequence data and association analysis with fatty acid composition and growth traits

Copy number variations (CNVs) are important genetic variants complementary to SNPs, and can be considered as biomarkers for some economically important traits in domestic animals. In the present study, a genomic analysis of porcine CNVs based on next-generation sequencing data was carried out to identify CNVs segregating in an Iberian x Landrace backcross population and study their association with fatty acid composition and growth-related traits. A total of 1,279 CNVs, including duplications and deletions, were detected, ranging from 106 to 235 CNVs across samples, with an average of 183 CNVs per sample. Moreover, we detected 540 CNV regions (CNVRs) containing 245 genes. Functional annotation suggested that these genes possess a great variety of molecular functions and may play a role in production traits in commercial breeds. Some of the identified CNVRs contained relevant functional genes (e.g., CLCA4, CYP4X1, GPAT2, MOGAT2, PLA2G2A and PRKG1, among others). The variation in copy number of four of them (CLCA4, GPAT2, MOGAT2 and PRKG1) was validated in 150 BC1_LD (25% Iberian and 75% Landrace) animals by qPCR. Additionally, their contribution regarding backfat and intramuscular fatty acid composition and growth–related traits was analyzed. Statistically significant associations were obtained for CNVR112 (GPAT2) for the C18:2(n-6)/C18:3(n-3) ratio in backfat and carcass length, among others. Notably, GPATs are enzymes that catalyze the first step in the biosynthesis of both triglycerides and glycerophospholipids, suggesting that this CNVR may contribute to genetic variation in fatty acid composition and growth traits. These findings provide useful genomic information to facilitate the further identification of trait-related CNVRs affecting economically important traits in pigs.

Integration of liver gene co-expression networks and eGWAs analyses highlighted candidate regulators implicated in lipid metabolism in pigs

In the present study, liver co-expression networks and expression Genome Wide Association Study (eGWAS) were performed to identify DNA variants and molecular pathways implicated in the functional regulatory mechanisms of meat quality traits in pigs. With this purpose, the liver mRNA expression of 44 candidates genes related with lipid metabolism was analysed in 111 Iberian x Landrace backcross animals. The eGWAS identified 92 eSNPs located in seven chromosomal regions and associated with eight genes: CROT, CYP2U1, DGAT1, EGF, FABP1, FABP5, PLA2G12A, and PPARA. Remarkably, cis-eSNPs associated with FABP1 gene expression which may be determining the C18:2(n-6)/C18:3(n-3) ratio in backfat through the multiple interaction of DNA variants and genes were identified. Furthermore, a hotspot on SSC8 associated with the gene expression of eight genes was identified and the TBCK gene was pointed out as candidate gene regulating it. Our results also suggested that the PI3K-Akt-mTOR pathway plays an important role in the control of the analysed genes highlighting nuclear receptors as the NR3C1 or PPARA. Finally, sex-dimorphism associated with hepatic lipid metabolism was identified with over-representation of female-biased genes. These results increase our knowledge of the genetic architecture underlying fat composition traits.

Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat

The aim of this work is to better understand the genetic mechanisms determining two complex traits affecting porcine meat quality: intramuscular fat (IMF) content and its fatty acid (FA) composition. With this purpose, expression Genome-Wide Association Study (eGWAS) of 45 lipid-related genes associated with meat quality traits in swine muscle (Longissimus dorsi) of 114 Iberian × Landrace backcross animals was performed. The eGWAS identified 241 SNPs associated with 11 genes: ACSM5, CROT, FABP3, FOS, HIF1AN, IGF2, MGLL, NCOA1, PIK3R1, PLA2G12A and PPARA. Three expression Quantitative Trait Loci (eQTLs) for IGF2, ACSM5 and MGLL were identified, showing cis-acting effects, whereas 16 eQTLs had trans regulatory effects. A polymorphism in the ACSM5 promoter region associated with its expression was identified. In addition, strong candidate genes regulating ACSM5, FOS, PPARA, PIK3R1, PLA2G12A and HIF1AN gene expression were also seen. Notably, the analysis highlighted the NR3C1 transcription factor as a strong candidate gene involved in the regulation of the 45 genes analysed. Finally, the IGF2, MGLL, MC2R, ARHGAP6, and NR3C1 genes were identified as potential regulators co-localizing within QTLs for fatness and growth traits in the IBMAP population. The results obtained increase our knowledge in the functional regulatory mechanisms involved in these complex traits.

Analysis of the porcine APOA2 gene expression in liver, polymorphism identification and association with fatty acid composition traits

APOA2 is a protein implicated in triglyceride, fatty acid and glucose metabolism. In pigs, the APOA2 gene is located on pig chromosome 4 (SSC4) in a QTL region affecting fatty acid composition, fatness and growth traits. In this study, we evaluated APOA2 as a candidate gene for meat quality traits in an Iberian × Landrace backcross population. The APOA2:c.131T>A polymorphism, located in exon 3 of APOA2 and determining a missense mutation, was associated with the percentage of hexadecenoic acid [C16:1(n-9)], linoleic acid [C18:2(n-6)], α-linolenic acid [C18:3(n-3)], dihomo-gamma-linolenic acid [C20:3(n-6)] and polyunsaturated fatty acids (PUFAs) in backfat. Furthermore, this SNP was associated with the global mRNA expression levels of APOA2 in liver and was used as a marker to determine allelic expression imbalance by pyrosequencing. We determined an overexpression of the T allele in heterozygous samples with a mean ratio of 2.8 (T/A), observing a high variability in the allelic expression among individuals. This result suggests that complex regulatory mechanisms, beyond a single polymorphism (e.g. epigenetic effects or multiple cis-acting polymorphisms), may be regulating APOA2 gene expression.

Epigenetic regulation of the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Background

In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs.

Results

High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3’UTR length. Although several SNPs in the 3’UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1.

Conclusions

Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0111-y) contains supplementary material, which is available to authorized users.

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

BACKGROUND

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F2 intercross population and a Chinese Sutai half-sib population.

RESULTS

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F2 animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

CONCLUSIONS

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Recombination of the porcine X chromosome: a high density linkage map

Background

Linkage maps are essential tools for the study of several topics in genome biology. High density linkage maps for the porcine autosomes have been constructed exploiting the high density data provided by the PorcineSNP60 BeadChip. However, a high density SSCX linkage map has not been reported up to date. The aim of the current study was to build an accurate linkage map of SSCX to provide precise estimates of recombination rates along this chromosome and creating a new tool for QTL fine mapping.

Results

A female-specific high density linkage map was built for SSCX using Sscrofa10.2 annotation. The total length of this chromosome was 84.61 cM; although the average recombination rate was 0.60 cM/Mb, both cold and hot recombination regions were identified. A Bayesian probabilistic to genetic groups and revealed that the animals used in the current study for linkage map construction were likely to be carriers of X chromosomes of European origin. Finally, the newly generated linkage map was used to fine-map a QTL at 16 cM for intramuscular fat content (IMF) measured on longissimus dorsi. The sulfatase isozyme S gene constitutes a functional and positional candidate gene underlying the QTL effect.

Conclusions

The current study presents for the first time a high density linkage map for SSCX and supports the presence of cold and hot recombination intervals along this chromosome. The large cold recombination region in the central segment of the chromosome is not likely to be due to structural differences between X chromosomes of European and Asian origin. In addition, the newly generated linkage map has allowed us to fine-map a QTL on SSCX for fat deposition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0148-x) contains supplementary material, which is available to authorized users.

A co-association network analysis of the genetic determination of pig conformation, growth and fatness

BACKGROUND

Several QTLs have been identified for major economically relevant traits in livestock, such as growth and meat quality, revealing the complex genetic architecture of these traits. The use of network approaches considering the interactions of multiple molecules and traits provides useful insights into the molecular underpinnings of complex traits. Here, a network based methodology, named Association Weight Matrix, was applied to study gene interactions and pathways affecting pig conformation, growth and fatness traits.

RESULTS

The co-association network analysis underpinned three transcription factors, PPARγ, ELF1, and PRDM16 involved in mesoderm tissue differentiation. Fifty-four genes in the network belonged to growth-related ontologies and 46 of them were common with a similar study for growth in cattle supporting our results. The functional analysis uncovered the lipid metabolism and the corticotrophin and gonadotrophin release hormone pathways among the most important pathways influencing these traits. Our results suggest that the genes and pathways here identified are important determining either the total body weight of the animal and the fat content. For instance, a switch in the mesoderm tissue differentiation may determinate the age-related preferred pathways being in the puberty stage those related with the miogenic and osteogenic lineages; on the contrary, in the maturity stage cells may be more prone to the adipocyte fate. Hence, our results demonstrate that an integrative genomic co-association analysis is a powerful approach for identifying new connections and interactions among genes.

CONCLUSIONS

This work provides insights about pathways and key regulators which may be important determining the animal growth, conformation and body proportions and fatness traits. Molecular information concerning genes and pathways here described may be crucial for the improvement of genetic breeding programs applied to pork meat production.

Molecular characterization and different expression patterns of the FABP gene family during goat skeletal muscle development

The FABP (adipocyte fatty acid-binding protein) genes play an important role in intracellular fatty acid transport and considered to be candidate genes for fatness traits in domestic animal. In this study, we cloned the cDNA sequences of goat FABP family genes and their expression patterns were detected by semi-quantitative RT-PCR and quantitative real time RT-PCR. Expression analysis showed that goat FABP1 gene was predominantly expressed in liver, kidney and large intestine. While FABP4 was widely expressed in many tissues with a high expression level was observed in the fat, skeletal muscle, stomach and lung. Notably, FABP2 gene was expressed specifically in small intestine. Moreover, goat FABP3 was expressed at 60 day with the highest level, then significantly (p < 0.01) decreased at the 90 day. No significant expression differences were observed in longissimus dorsi muscles among 3 day, 30 day and 60 day. Goat FABP4 was expressed at 3 day with the lowest level, then significantly (p < 0.01) increased to a peak at the 60 day. In addition, a significant relationship between FABP3 mRNA expression levels and intramuscular fat (IMF) content was observed. These results suggest that the FABP3 and FABP4 may be important genes for meat quality and provides useful information for further studies on their roles in skeletal muscle IMF deposit.

Differences in muscle transcriptome among pigs phenotypically extreme for fatty acid composition

BACKGROUND

Besides having an impact on human health, the porcine muscle fatty acid profile determines meat quality and taste. The RNA-Seq technologies allowed us to explore the pig muscle transcriptome with an unprecedented detail. The aim of this study was to identify differentially-expressed genes between two groups of 6 sows belonging to an Iberian × Landrace backcross with extreme phenotypes according to FA profile.

RESULTS

We sequenced the muscle transcriptome acquiring 787.5 M of 75 bp paired-end reads. About 85.1% of reads were mapped to the reference genome. Of the total reads, 79.1% were located in exons, 6.0% in introns and 14.9% in intergenic regions, indicating expressed regions not annotated in the reference genome. We identified a 34.5% of the intergenic regions as interspersed repetitive regions. We predicted a total of 2,372 putative proteins. Pathway analysis with 131 differentially-expressed genes revealed that the most statistically-significant metabolic pathways were related with lipid metabolism. Moreover, 18 of the differentially-expressed genes were located in genomic regions associated with IMF composition in an independent GWAS study in the same genetic background. Thus, our results indicate that the lipid metabolism of FAs is differently modulated when the FA composition in muscle differs. For instance, a high content of PUFA may reduce FA and glucose uptake resulting in an inhibition of the lipogenesis. These results are consistent with previous studies of our group analysing the liver and the adipose tissue transcriptomes providing a view of each of the main organs involved in lipid metabolism.

CONCLUSIONS

The results obtained in the muscle transcriptome analysis increase the knowledge of the gene regulation of IMF deposition, FA profile and meat quality, in terms of taste and nutritional value. Besides, our results may be important in terms of human health.

New insight into the SSC8 genetic determination of fatty acid composition in pigs

Background

Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F₂ intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F₂ intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8.

Results

Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat.

Conclusions

Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.

Copy number variants in a highly inbred Iberian porcine strain

We carried out a comprehensive genomic analysis of porcine copy number variants (CNVs) based on whole-genome SNP genotyping data and provided new measures of genomic diversity (number, length and distribution of CNV events) for a highly inbred strain (the Guadyerbas strain). This strain represents one of the most ancient surviving populations of the Iberian breed, and it is currently in serious danger of extinction. CNV detection was conducted on the complete Guadyerbas population, adjusted for genomic waves, and used strict quality criteria, pedigree information and the latest porcine genome annotation. The analysis led to the detection of 65 CNV regions (CNVRs). These regions cover 0.33% of the autosomal genome of this particular strain. Twenty-nine of these CNVRs were identified here for the first time. The relatively low number of detected CNVRs is in line with the low variability and high inbreeding estimated previously for this Iberian strain using pedigree, microsatellite or SNP data. A comparison across different porcine studies has revealed that more than half of these regions overlap with previously identified CNVRs or multicopy regions. Also, a preliminary analysis of CNV detection using whole-genome sequence data for four Guadyerbas pigs showed overlapping for 16 of the CNVRs, supporting their reliability. Some of the identified CNVRs contain relevant functional genes (e.g., the SCD and USP15 genes), which are worth being further investigated because of their importance in determining the quality of Iberian pig products. The CNVR data generated could be useful for improving the porcine genome annotation.

Fine mapping and single nucleotide polymorphism effects estimation on pig chromosomes 1, 4, 7, 8, 17 and X

Fine mapping of quantitative trait loci (QTL) from previous linkage studies was performed on pig chromosomes 1, 4, 7, 8, 17, and X which were known to harbor QTL. Traits were divided into: growth performance, carcass, internal organs, cut yields, and meat quality. Fifty families were used of a F2 population produced by crossing local Brazilian Piau boars with commercial sows. The linkage map consisted of 237 SNP and 37 microsatellite markers covering 866 centimorgans. QTL were identified by regression interval mapping using GridQTL. Individual marker effects were estimated by Bayesian LASSO regression using R. In total, 32 QTL affecting the evaluated traits were detected along the chromosomes studied. Seven of the QTL were known from previous studies using our F2 population, and 25 novel QTL resulted from the increased marker coverage. Six of the seven QTL that were significant at the 5% genome-wide level had SNPs within their confidence interval whose effects were among the 5% largest effects. The combined use of microsatellites along with SNP markers increased the saturation of the genome map and led to smaller confidence intervals of the QTL. The results showed that the tested models yield similar improvements in QTL mapping accuracy.

Genome-wide analysis of porcine backfat and intramuscular fat fatty acid composition using high-density genotyping and expression data

Background

Porcine fatty acid composition is a key factor for quality and nutritive value of pork. Several QTLs for fatty acid composition have been reported in diverse fat tissues. The results obtained so far seem to point out different genetic control of fatty acid composition conditional on the fat deposits. Those studies have been conducted using simple approaches and most of them focused on one single tissue. The first objective of the present study was to identify tissue-specific and tissue-consistent QTLs for fatty acid composition in backfat and intramuscular fat, combining linkage mapping and GWAS approaches and conducted under single and multitrait models. A second aim was to identify powerful candidate genes for these tissue-consistent QTLs, using microarray gene expression data and following a targeted genetical genomics approach.

Results

The single model analyses, linkage and GWAS, revealed over 30 and 20 chromosomal regions, 24 of them identified here for the first time, specifically associated to the content of diverse fatty acids in BF and IMF, respectively. The analyses with multitrait models allowed identifying for the first time with a formal statistical approach seven different regions with pleiotropic effects on particular fatty acids in both fat deposits. The most relevant were found on SSC8 for C16:0 and C16:1(n-7) fatty acids, detected by both linkage and GWAS approaches. Other detected pleiotropic regions included one on SSC1 for C16:0, two on SSC4 for C16:0 and C18:2, one on SSC11 for C20:3 and the last one on SSC17 for C16:0. Finally, a targeted eQTL scan focused on regions showing tissue-consistent effects was conducted with Longissimus and fat gene expression data. Some powerful candidate genes and regions were identified such as the PBX1, RGS4, TRIB3 and a transcription regulatory element close to ELOVL6 gene to be further studied.

Conclusions

Complementary genome scans have confirmed several chromosome regions previously associated to fatty acid composition in backfat and intramuscular fat, but even more, to identify new ones. Although most of the detected regions were tissue-specific, supporting the hypothesis that the major part of genes affecting fatty acid composition differs among tissues, seven chromosomal regions showed tissue-consistent effects. Additional gene expression analyses have revealed powerful target regions to carry the mutation responsible for the pleiotropic effects.

Electronic supplementary material

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

Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition

Background

In pigs, adipose tissue is one of the principal organs involved in the regulation of lipid metabolism. It is particularly involved in the overall fatty acid synthesis with consequences in other lipid-target organs such as muscles and the liver. With this in mind, we have used massive, parallel high-throughput sequencing technologies to characterize the porcine adipose tissue transcriptome architecture in six Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition (three per group).

Results

High-throughput RNA sequencing was used to generate a whole characterization of adipose tissue (backfat) transcriptome. A total of 4,130 putative unannotated protein-coding sequences were identified in the 20% of reads which mapped in intergenic regions. Furthermore, 36% of the unmapped reads were represented by interspersed repeats, SINEs being the most abundant elements. Differential expression analyses identified 396 candidate genes among divergent animals for intramuscular fatty acid composition. Sixty-two percent of these genes (247/396) presented higher expression in the group of pigs with higher content of intramuscular SFA and MUFA, while the remaining 149 showed higher expression in the group with higher content of PUFA. Pathway analysis related these genes to biological functions and canonical pathways controlling lipid and fatty acid metabolisms. In concordance with the phenotypic classification of animals, the major metabolic pathway differentially modulated between groups was de novo lipogenesis, the group with more PUFA being the one that showed lower expression of lipogenic genes.

Conclusions

These results will help in the identification of genetic variants at loci that affect fatty acid composition traits. The implications of these results range from the improvement of porcine meat quality traits to the application of the pig as an animal model of human metabolic diseases.

Genome-wide association analyses for fatty acid composition in porcine muscle and abdominal fat tissues

Fatty acid composition is an important phenotypic trait in pigs as it affects nutritional, technical and sensory quality of pork. Here, we reported a genome-wide association study (GWAS) for fatty acid composition in the longissimus muscle and abdominal fat tissues of 591 White Duroc×Erhualian F2 animals and in muscle samples of 282 Chinese Sutai pigs. A total of 46 loci surpassing the suggestive significance level were identified on 15 pig chromosomes (SSC) for 12 fatty acids, revealing the complex genetic architecture of fatty acid composition in pigs. Of the 46 loci, 15 on SSC5, 7, 14 and 16 reached the genome-wide significance level. The two most significant SNPs were ss131535508 (P = 2.48×10(-25)) at 41.39 Mb on SSC16 for C20∶0 in abdominal fat and ss478935891 (P = 3.29×10(-13)) at 121.31 Mb on SSC14 for muscle C18∶0. A meta-analysis of GWAS identified 4 novel loci and enhanced the association strength at 6 loci compared to those evidenced in a single population, suggesting the presence of common underlying variants. The longissimus muscle and abdominal fat showed consistent association profiles at most of the identified loci and distinct association signals at several loci. All loci have specific effects on fatty acid composition, except for two loci on SSC4 and SSC7 affecting multiple fatness traits. Several promising candidate genes were found in the neighboring regions of the lead SNPs at the genome-wide significant loci, such as SCD for C18∶0 and C16∶1 on SSC14 and ELOVL7 for C20∶0 on SSC16. The findings provide insights into the molecular basis of fatty acid composition in pigs, and would benefit the final identification of the underlying mutations.

Molecular characterization, expression patterns, and association analysis with carcass traits of porcine USF1 gene

The upstream stimulatory factor 1 (USF1) is a transcription factor controlling expression of several genes involved in lipid and glucose homeostasis. In this study, two isoforms of the porcine USF1 gene were detected by reverse transcription polymerase chain reaction (RT-PCR), termed USF1 wild-type (wt) and USF1/CD, both of them contain a helix-loop-helix leucine zipper (HLH-LZ) conserved domain. Tissue distribution analysis showed that the two transcripts of porcine USF1 gene were ubiquitously expressed in all tested tissues, except for heart. Moreover, we found that a single nucleotide polymorphism (SNP, C/T) in intron 10 was significantly associated with ratio of lean to fat (P < 0.05), dress percentage (P < 0.05), average backfat thickness (P < 0.05), loin eye width (P < 0.05), lean meat percentage (P < 0.01), loin eye height (P < 0.01), and loin eye area (P < 0.01). This result suggests that porcine USF1 gene may be a candidate gene of meat production trait.

Detecção de locos de características quantitativas nos cromossomos 1, 2, 3, 12, 14, 15 e X de suínos: características de desempenho

Mapeou-se quantitative trait loci (QTL) associados a características de desempenho nos cromossomos 1, 2, 3, 12, 14, 15 e X de suínos pertencentes a uma população F2, formada a partir do cruzamento entre dois machos da raça naturalizada brasileira Piau e 18 fêmeas comerciais (Landrace x Large White x Pietrain). O mapa genético de ligação da população foi construído após a genotipagem dos animais para 35 marcadores microssatélites. As estimativas do conteúdo de informação polimórfica indicaram que os marcadores microssatélites foram adequados para as análises de QTL. Os dados foram analisados pelo mapeamento por intervalo usando-se o programa GridQTL. Encontraram-se seis QTL, sendo que o QTL genômico para idade ao abate atingiu a significância de 5% de probabilidade. As informações dos QTL detectados neste estudo são úteis para identificar genes que podem ser usados em conjunto com os métodos convencionais de seleção, aumentar a acurácia deles e prover uma compreensão dos fenótipos produtivos de suínos.

Polymorphism in the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Background

The ELOVL fatty acid elongase 6 (ELOVL6), the only elongase related to de novo lipogenesis, catalyzes the rate-limiting step in the elongation cycle by controlling the fatty acid balance in mammals. It is located on pig chromosome 8 (SSC8) in a region where a QTL affecting palmitic, and palmitoleic acid composition was previously detected, using an Iberian x Landrace intercross. The main goal of this work was to fine-map the QTL and to evaluate the ELOVL6 gene as a positional candidate gene affecting the percentages of palmitic and palmitoleic fatty acids in pigs.

Methodology and Principal Findings

The combination of a haplotype-based approach and single-marker analysis allowed us to identify the main, associated interval for the QTL, in which the ELOVL6 gene was identified and selected as a positional candidate gene. A polymorphism in the promoter region of ELOVL6, ELOVL6:c.-533C>T, was highly associated with the percentage of palmitic and palmitoleic acids in muscle and backfat. Significant differences in ELOVL6 gene expression were observed in backfat when animals were classified by the ELOVL6:c.-533C>T genotype. Accordingly, animals carrying the allele associated with a decrease in ELOVL6 gene expression presented an increase in C16:0 and C16:1(n-7) fatty acid content and a decrease of elongation activity ratios in muscle and backfat. Furthermore, a SNP genome-wide association study with ELOVL6 relative expression levels in backfat showed the strongest effect on the SSC8 region in which the ELOVL6 gene is located. Finally, different potential genomic regions associated with ELOVL6 gene expression were also identified by GWAS in liver and muscle, suggesting a differential tissue regulation of the ELOVL6 gene.

Conclusions and Significance

Our results suggest ELOVL6 as a potential causal gene for the QTL analyzed and, subsequently, for controlling the overall balance of fatty acid composition in pigs.

Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition

BACKGROUND

New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq.

RESULTS

The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism.

CONCLUSIONS

In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

Molecular characterization of differentially expressed TXNIP gene and its association with porcine carcass traits

Thioredoxin interacting protein (TXNIP), which plays a regulatory role in lipid metabolism and immune regulation, is down-regulated expressed in F(1) hybrids Landrace × Yorkshire skeletal muscle. Here we described the molecular characterization of porcine TXNIP gene. The full-length cDNA contains a coding sequence of 1,176 bp nucleotides with untranslated regions of 263 bp at 5'-end and 441 bp at 3'-end, respectively. The predicted molecular mass and isoelectric point of porcine TXNIP is 43.81 kDa and 7.385, respectively. The deduced 391 amino acids exhibit high identity with other mammalian TXNIP. The TXNIP gene contains eight coding exons and seven non coding introns, spans approximately 3,348 bp. The expression of porcine TXNIP mRNA is almost absent in Landrace × Yorkshire and lower level in 6-month-old pigs during skeletal muscle development. Other stages and breeds were high level expressed. Statistical analysis showed the TXNIP gene polymorphism (c.575-4T>C) was different between F(1) hybrids and their parents, was highly associated with dressing percentage (DP) and thorax-waist fat thickness (TFT) in the Yorkshire × Meishan F(2) population. The possible role of TXNIP was discussed.

The Human Obesity Gene Map: The 2003 Update

This is the tenth update of the human obesity gene map, incorporating published results up to the end of October 2003 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. Transgenic and knockout murine models relevant to obesity are also incorporated (N = 55). As of October 2003, 41 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. QTLs reported from animal models currently number 183. There are 208 human QTLs for obesity phenotypes from genome-wide scans and candidate regions in targeted studies. A total of 35 genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 272 studies reporting positive associations with 90 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, more than 430 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

The Human Obesity Gene Map: The 2004 Update

This paper presents the eleventh update of the human obesity gene map, which incorporates published results up to the end of October 2004. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTLs) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2004, 173 human obesity cases due to single-gene mutations in 10 different genes have been reported, and 49 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 166 genes which, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 221. The number of human obesity QTLs derived from genome scans continues to grow, and we have now 204 QTLs for obesity-related phenotypes from 50 genome-wide scans. A total of 38 genomic regions harbor QTLs replicated among two to four studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably with 358 findings of positive associations with 113 candidate genes. Among them, 18 genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, >600 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful publications and genomic and other relevant sites can be found at http://obesitygene.pbrc.edu.

Mapeamento de QTL nos cromossomos 1, 2, 3, 12, 14, 15 e X em suínos: características de carcaça e qualidade de carne

A realização do presente estudo teve como objetivo mapear Quantitative Trait Loci (QTL) de carcaça e qualidade de carne em uma população F2 de suínos desenvolvida pelo cruzamento de dois reprodutores da raça brasileira Piau com 18 fêmeas comerciais (Landrace x Large White x Pietrain). O mapa de ligação para essa população foi construído após a genotipagem de 684 animais para 35 marcadores microssatélites. Os dados foram analisados pelo mapeamento por intervalo usando-se sexo, lote e genótipo halotano como efeitos fixos e peso de carcaça ao abate, peso da carcaça direita e idade ao abate como covariáveis. Um total de 18 QTLs foi encontrado; os QTLs para maior espessura de toucinho na região da copa, na linha dorsolombar, e a perda por cozimento foram significativos em nível de 5% genômico. A característica espessura de toucinho foi essencialmente associada aos alelos da raça Piau, conhecido como porco tipo banha. As informações dos QTLs significativos encontrados servem para futuros estudos de mapeamento fino para identificação de genes a serem usados em conjunto com os métodos tradicionais de seleção, para melhorar a eficiência dos programas de melhoramento, assim como prover informação acerca da fisiologia envolvida no desenvolvimento das características quantitativas dos suínos.

Detection of quantitative trait loci for growth- and fatness-related traits in a large-scale White Duroc × Erhualian intercross pig population

Growth and fatness are economically important traits in pigs. In this study, a genome scan was performed to detect quantitative trait loci (QTL) for 14 growth and fatness traits related to body weight, backfat thickness and fat weight in a large-scale White Duroc × Erhualian F(2) intercross. A total of 76 genome-wide significant QTL were mapped to 16 chromosomes. The most significant QTL was found on pig chromosome (SSC) 7 for fatness with unexpectedly small confidence intervals of ∼2 cM, providing an excellent starting point to identify causal variants. Common QTL for both fatness and growth traits were found on SSC4, 5, 7 and 8, and shared QTL for fat deposition were detected on SSC1, 2 and X. Time-series analysis of QTL for body weight at six growth stages revealed the continuously significant effects of the QTL on SSC4 at the fattening period and the temporal-specific expression of the QTL on SSC7 at the foetus and fattening stages. For fatness traits, Chinese Erhualian alleles were associated with increased fat deposition except that at the major QTL on SSC7. For growth traits, most of White Duroc alleles enhanced growth rates except for those at three significant QTL on SSC6, 7 and 9. The results confirmed many previously reported QTL and also detected novel QTL, revealing the complexity of the genetic basis of growth and fatness in pigs.

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ∼ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ∼ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.

A genome scan for quantitative trait loci affecting milk somatic cell score in Israeli and Italian Holstein cows by means of selective DNA pooling with single- and multiple-marker mapping

Mastitis is an important and common dairy cattle disease affecting milk yield, quality, and consumer safety as well as cheese yields and quality. Animal welfare and residues of the antibiotics used to treat mastitis cause public concern. Considerable genetic variation may allow selection for increased resistance to mastitis. Because of high genetic correlation to milk somatic cell score (SCS), SCS can serve as a surrogate trait for mastitis resistance. The present study intended to identify quantitative trait loci (QTL) affecting SCS in Israeli and Italian Holstein dairy cattle (IsH and ItH, respectively), using selective DNA pooling with single and multiple marker mapping. Milk samples of 4,788 daughters of 6 IsH and 7 ItH sires were used to construct sire-family high- and low-tail pools, which were genotyped at 123 (IsH) and 133 (ItH) microsatellite markers. Shadow correction was used to obtain pool allele frequency estimates. Frequency difference between the tails and empirical standard error of D, SE(D), were used to obtain P-values. All markers significant by single marker mapping were also significant by multiple marker mapping, but not vice versa. Combining both populations, 22 QTL on 21 chromosomes were identified; all corresponded to previous reports in the literature. Confidence intervals were set by chi-squared drop method. Heterozygosity of QTL was estimated at 44.2%. Allele substitution effects ranged from 1,782 to 4,930 cells/mL in estimated breeding value somatic cell count units. Most (80%) of the observed variation in estimated breeding value somatic cell score could be explained by the QTL identified under the stringent criteria. The results found here can be used as a basis for further genome-wide association studies for the same trait.

Meta-analysis of results from quantitative trait loci mapping studies on pig chromosome 4

Meta-analysis of results from multiple studies could lead to more precise quantitative trait loci (QTL) position estimates compared to the individual experiments. As the raw data from many different studies are not readily available, the use of results from published articles may be helpful. In this study, we performed a meta-analysis of QTL on chromosome 4 in pig, using data from 25 separate experiments. First, a meta-analysis was performed for individual traits: average daily gain and backfat thickness. Second, a meta-analysis was performed for the QTL of three traits affecting loin yield: loin eye area, carcass length and loin meat weight. Third, 78 QTL were selected from 20 traits that could be assigned to one of three broad categories: carcass, fatness or growth traits. For each analysis, the number of identified meta-QTL was smaller than the number of initial QTL. The reduction in the number of QTL ranged from 71% to 86% compared to the total number before the meta-analysis. In addition, the meta-analysis reduced the QTL confidence intervals by as much as 85% compared to individual QTL estimates. The reduction in the confidence interval was greater when a large number of independent QTL was included in the meta-analysis. Meta-QTL related to growth and fatness were found in the same region as the FAT1 region. Results indicate that the meta-analysis is an efficient strategy to estimate the number and refine the positions of QTL when QTL estimates are available from multiple populations and experiments. This strategy can be used to better target further studies such as the selection of candidate genes related to trait variation.

Quantitative trait loci for meat quality traits in pigs considering imprinting and epistatic effects

The aim of the research was to gain a better understanding of the genomic regulation of meat quality by investigating individual and epistatic QTL in a three-generation full-sib population (Pietrain x crossbred dam line). In total, 386 animals were genotyped for 96 markers. Analysed traits included pH, reflectance value, conductivity, and meat colour. Thirteen significant individual QTL were identified. The most significant QTL were detected on SSC1 and SSC9 for pH, on SSC4 for meat colour, and on SSC8 for conductivity, accounting for 3.4% to 4.7% of the phenotypic variance. Nine significant epistatic QTL pairs were detected accounting for between 5.7% and 10.9% of the phenotypic variance. Epistatic QTL pairs showing the largest effects were for reflectance value between two locations of SSC4, and for pH between SSC10 and SSC13, explaining 9.5% and 10.9% of the phenotypic variance, respectively. This study indicates that meat quality traits are influenced by numerous QTL as well as a complex network of interactions.