Combined line-cross and half-sib QTL analysis in Duroc-Pietrain population

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Weighted gene co-expression network analysis reveals potential candidate genes affecting drip loss in pork

Drip loss is an essential evaluation indicator for pork quality. It is closely related to other meat quality indicators, including water-holding capacity, water loss rate and pH value at 45 min (pH₁ ) and 24 h post-mortem (pH₂ ), and is influenced by environmental and genetic factors and their interactions. We previously conducted differentially expressed gene analysis to identify candidate genes affecting drip loss using eight individuals with extremely high- and low-drip loss selected from 28 purebred Duroc pigs. Using 28 identical samples, in the present study, we performed weighted gene co-expression network analysis with drip loss and drip loss-related traits, including water-holding capacity, water loss rate, pH₁ and pH₂ . A total of 25 modules were identified, and five of them correlated with at least two drip loss or drip loss-related traits. After functional enrichment analysis of genes in the five modules, three modules were found to be critical, as their genes were significantly involved in amino acid metabolism, immune response and apoptosis, which have potential relationships with drip loss. Furthermore, we identified five candidate genes affecting drip loss in one critical module, AASS, BCKDHB, ALDH6A1, MUT and MCCC1, as they overlapped with differentially expressed genes detected in our previous study, exhibited protein-protein interactions and had potential biological functions in affecting drip loss according to the literature. The outcomes of the present study enhance our understanding of the molecular mechanisms underlying drip loss and will aid in improving the pork quality.

Haplotype Block Analysis Reveals Candidate Genes and QTLs for Meat Quality and Disease Resistance in Chinese Jiangquhai Pig Breed

The Jiangquhai (JQ) pig breed is one of the most widely recognized pig populations in China due to its unique and dominant characteristics. In this study, we examined the extent of Linkage disequilibrium (LD) and haplotype block structure of the JQ pig breed, and scanned the blocks for possible genes underlying important QTLs that could either be responsible for some adaptive features in these pigs or might have undergone some selection pressure. We compared some of our results with other Chinese and Western pig breeds. The results show that the JQ breed had the highest total block length (349.73 Mb ≈ 15% of its genome), and the coverage rate of blocks in most of its chromosomes was larger than those of other breeds except for Sus scrofa chromosome 4 (SSC4), SSC6, SSC7, SSC8, SSC10, SSC12, SSC13, SSC14, SSC17, SSC18, and SSCX. Moreover, the JQ breed had more SNPs that were clustered into haplotype blocks than the other breeds examined in this study. Our shared and unique haplotype block analysis revealed that the Hongdenglong (HD) breed had the lowest percentage of shared haplotype blocks while the Shanzhu (SZ) breed had the highest. We found that the JQ breed had an average r² > 0.2 at SNPs distances 10–20 kb and concluded that about 120,000–240,000 SNPs would be needed for a successful GWAS in the breed. Finally, we detected a total of 88 genes harbored by selected haplotype blocks in the JQ breed, of which only 4 were significantly enriched (p-value ≤ 0.05). These genes were significantly enriched in 2 GO terms (p-value < 0.01), and 2 KEGG pathways (p-value < 0.02). Most of these enriched genes were related to health. Also, most of the overlapping QTLs detected in the haplotype blocks were related to meat and carcass quality, as well as health, with a few of them relating to reproduction and production. These results provide insights into the genetic architecture of some adaptive and meat quality traits observed in the JQ pig breed and also revealed the pattern of LD in the genome of the pig. Our result provides significant guidance for improving the statistical power of GWAS and optimizing the conservation strategy for this JQ pig breed.

Genome-wide association study for backfat thickness at 100 kg and loin muscle thickness in domestic pigs based on genotyping by sequencing

Both backfat thickness at 100 kg (B100) and loin muscle thickness (LMT) are economically important traits in pigs. In this study, a total of 1,200 pigs (600 Landrace and 600 Yorkshire pigs) were examined with genotyping by sequencing. A total of 345,570 single nucleotide polymorphisms (SNPs) were obtained from 1,200 pigs. Then, a single marker regression test was used to conduct a genome-wide association study for B100 and LMT. A total of 8 and 90 significant SNPs were detected for LMT and B100, respectively. Interestingly, two shared significant loci [located at Sus scrofa chromosome (SSC) 6: 149876694 and SSC12: 46226580] were detected in two breeds for B100. Furthermore, three potential candidate genes were found for LMT and B100. The positional candidate gene FAM3C (SSC18: 25573656, P = 2.48 × 10^-9), which controls the survival, growth, and differentiation of tissues and cells, was found for LMT in Landrace pigs. At SSC9: 6.78-6.82 Mb in Landrace pigs, the positional candidate gene, INPPL1, which has a negative regulatory effect on diet-induced obesity and is involved in the regulation of insulin function, was found for B100. The candidate gene, RAB35, which regulates the adipocyte glucose transporter SLC2A4/GLUT4, was identified at approximately SSC14: 40.09-40.13 Mb in Yorkshire pigs. The results of this GWAS will greatly advance our understanding of the genetic architecture of the LMT and B100 traits. However, these identified loci and genes need to be further verified in more pig populations, and their functions also need to be validated by more biological experiments in pigs.

Identification of porcine CTLA4 gene polymorphism and their association with piglet diarrhea and performance traits

The objective of this study was to evaluate the association between the cytotoxic T-lymphocyte-associated antigen 4 (CTLA4) gene and piglet diarrhea. In this study, the mRNA expression of the CTLA4 gene increased significantly in IPEC-J2 cells after Escherichia coli K88 infection. Single nucleotide polymorphisms (SNPs) located in the 5' flanking region (SNPs g.107281989C>T) and 3'-untranslated region (3'-UTR; SNPs g.107288753C>A) were identified, and they were in linkage disequilibrium in both Min pigs and the Landrace population. Association analysis showed that Landrace piglets with a TT or AA genotype had a lower diarrhea index, and AA animals had higher average daily gain when compared to CC pigs, respectively (p < 0.05). However, the relationship between SNPs and diarrhea and performance traits in the Min population was not significant. Haplotype analysis indicated that the TC haplotype had the lowest diarrhea index. The 5' flanking deletion assay suggested that SNP g.107281989C>T was a molecular marker instead of the functional marker. This research demonstrated that genetic variances in the CTLA4 gene had significant effects on Landrace piglet diarrhea resistance.

Linkage disequilibrium pattern and genome-wide association mapping for meat traits in multiple porcine F2 crosses

In the present study, data from four F2 crosses were analysed and used to study the linkage disequilibrium (LD) structure within and across the crosses. Genome-wide association analyses (GWASes) for conductivity and dressing out meat traits were conducted using single-marker and Bayesian multi-marker models using the pooled data from all F2 crosses. Porcine F2 crosses generated from the distantly related founder breeds Wild Boar, Piétrain and Meishan, as well as from a porcine F2 cross from the closely related founder breed Piétrain and an F1 Large White × Landrace cross were pooled. A total of 2572 F2 animals were genotyped using a 62K SNP chip. The positions of the SNPs were based on genome assembly Sscrofa11.1. After post-alignment and genotype filtering, approximately 50K SNPs were usable for LD studies and GWASes. The main findings of the present study are that the breakdown of LD was faster in crosses from closely related founder breeds compared to crosses from distantly related founders. The fastest breakdown of LD was observed by pooling the data. Based on the single-marker results and LD structure, clusters and windows were built for 1-Mb intervals. For conductivity and dressing out, 183 and 191 nominal significant associations respectively and six and five clusters respectively were found. Dominance was important for conductivity, and considering dominance in GWASes improved the mapping signals. Most clear signals were found for conductivity on SSC6, 8 and 15 and for dressing out on SSC2 and 7. Considering dominance might contribute to the accuracy of genomic selection and serve as a guide for choosing mating pairs with good combining abilities. However, further research is needed to investigate if dominance is also important in crossbreed pig breeding schemes.

Detection of Selection Signatures in Chinese Landrace and Yorkshire Pigs Based on Genotyping-by-Sequencing Data

The domestic pigs have been undergone intense selection pressures for these development of interested traits following domestication and modern breeding. This has altered many traits in most of pig breeds, such as growth rate, body weight, fertility, and immunity. Thus, the objectives of this study were to (1) detect these selection signatures and identify the candidate genes that show evidences of recent artificial selection at the level of whole genome, (2) be beneficial to understand the relationship between genomic structure and phenotypic diversity, and (3) highlight the key roles of these candidate genes in growth and development in the two breeds. The data consisted of total raw number of 345570 single nucleotide polymorphisms (SNPs) in 1200 individuals from the Chinese Landrace pigs (L, n = 600) and Yorkshire pigs (Y, n = 600). Based on these SNPs data, two complementary methods, population differentiation (Fst) and composite likelihood ratio test (CLR), were carried out to detect the selection signatures in this study. A total of 540 potential selection regions (50 kb) which contained 111 candidate genes were detected for Landrace-Yorkshire pair (L-Y) by Fst. In addition, 73 and 125 candidate genes were found for Landrace pigs and Yorkshire pigs by CLR test based on 321 and 628 potential selection regions, respectively. Some candidate genes are associated with important traits and signaling pathways including the ACACA, MECR, COL11A1, GHR, IGF1R, IGF2R, IFNG, and MTOR gene. The ACACA and MECR gene are related to fatty acid biosynthesis. The COL11A1 gene is essential for the development of the normal differentiation. The GHR, IGF1R, and IGF2R gene are significant candidate genes which play major roles in the growth and development in animals. The IFNG gene is associated with some aspects of immune response. The MTOR gene regulates many signaling pathways and signaling transduction pathway.

Genomic regions associated with host response to porcine reproductive and respiratory syndrome vaccination and co-infection in nursery pigs

BACKGROUND

The WUR1000125 (WUR) single nucleotide polymorphism (SNP) can be used as a genetic marker for host response to porcine reproductive and respiratory syndrome (PRRS), PRRS vaccination, and co-infection with porcine circovirus type 2b (PCV2b). Objectives of this study were to identify genomic regions other than WUR associated with host response to PRRS vaccination and PRRSV/PCV2b co-infection and regions with a different effect on host response to co-infection, depending on previous vaccination for PRRS.

METHODS

Commercial crossbred nursery pigs were pre-selected for WUR genotype (n = 171 AA and 198 AB pigs) where B is the dominant and favorable allele. Half of the pigs were vaccinated for PRRS and 4 weeks later, all pigs were co-infected with PRRS virus and PCV2b. Average daily gain (ADG) and viral load (VL) were quantified post vaccination (Post Vx) and post co-infection (Post Co-X). Single-SNP genome-wide association analyses were then conducted to identify genomic regions associated with response to vaccination and co-infection.

RESULTS

Multiple SNPs near the major histocompatibility complex were significantly associated with PCV2b VL (-log ₁₀ P ≥ 5.5), regardless of prior vaccination for PRRS. Several SNPs were also significantly associated with ADG Post Vx and Post Co-X. SNPs with a different effect on ADG, depending on prior vaccination for PRRS, were identified Post Vx (-log ₁₀ P = 5.6) and Post Co-X (-log ₁₀ P = 5.5). No SNPs were significantly associated with vaccination VL (-log₁₀ P ≤ 4.7) or PRRS VL (-log₁₀ P ≤ 4.3). Genes near SNPs associated with vaccination VL, PRRS VL, and PCV2b VL were enriched (P ≤ 0.01) for immune-related pathways and genes near SNPs associated with ADG were enriched for metabolism pathways (P ≤ 0.04). SNPs associated with vaccination VL, PRRS VL, and PCV2b VL showed overrepresentation of health QTL identified in previous studies and SNPs associated with ADG Post Vx of Non-Vx pigs showed overrepresentation of growth QTL.

CONCLUSIONS

Multiple genomic regions were associated with PCV2b VL and ADG Post Vx and Post Co-X. Different SNPs were associated with ADG, depending on previous vaccination for PRRS. Results of functional annotation analyses and novel approaches of using previously-reported QTL support the identified regions.

Genome-wide association study reveals genetic loci and candidate genes for average daily gain in Duroc pigs

Objective

Average daily gain (ADG) is an important target trait of pig breeding programs. We aimed to identify single nucleotide polymorphisms (SNPs) and genomic regions that are associated with ADG in the Duroc pig population.

Methods

We performed a genome-wide association study involving 390 Duroc boars and by using the PorcineSNP60K Beadchip and two linear models.

Results

After quality control, we detected 3,5971 SNPs, which included seven SNPs that are significantly associated with the ADG of pigs. We identified six quantitative trait loci (QTL) regions for ADG. These QTLs included four previously reported QTLs on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, and SSC13, as well as two novel QTLs on SSC6 and SSC16. In addition, we selected six candidate genes (general transcription factor 3C polypeptide 5, high mobility group AT-hook 2, nicotinamide phosphoribosyltransferase, oligodendrocyte transcription factor 1, pleckstrin homology and RhoGEF domain containing G4B, and ENSSSCG00000031548) associated with ADG on the basis of their physiological roles and positional information. These candidate genes are involved in skeletal muscle cell differentiation, diet-induced obesity, and nervous system development.

Conclusion

This study contributes to the identification of the casual mutation that underlies QTLs associated with ADG and to future pig breeding programs based on marker-assisted selection. Further studies are needed to elucidate the role of the identified candidate genes in the physiological processes involved in ADG regulation.

Genomewide association of piglet responses to infection with one of two porcine reproductive and respiratory syndrome virus isolates

Porcine reproductive and respiratory syndrome (PRRS) is a devastating disease in the swine industry. Identification of host genetic factors that enable selection for improved performance during PRRS virus (PRRSV) infection would reduce the impact of this disease on animal welfare and production efficiency. We conducted genomewide association study (GWAS) analyses of data from 13 trials of approximately 200 commercial crossbred nursery-age piglets that were experimentally infected with 1 of 2 type 2 isolates of PRRSV (NVSL 97-7985 [NVSL] and KS2006-72109 [KS06]). Phenotypes analyzed were viral load (VL) in blood during the first 21 d after infection (dpi) and weight gain (WG) from 0 to 42 dpi. We accounted for the previously identified QTL in the region on SSC4 in our models to increase power to identify additional regions. Many regions identified by single-SNP analyses were not identified using Bayes-B, but both analyses identified the same regions on SSC3 and SSC5 to be associated with VL in the KS06 trials and on SSC6 in the NVSL trials ( < 5 × 10); for WG, regions on SSC5 and SSC17 were associated in the NVSL trials ( < 3 × 10). No regions were identified with either method for WG in the KS06 trials. Except for the region on SSC4, which was associated with VL for both isolates (but only with WG for NVSL), identified regions did not overlap between the 2 PRRSV isolate data sets, despite high estimates of the genetic correlation between isolates for traits based on these data. We also identified genomic regions whose associations with VL or WG interacted with either PRRSV isolate or with genotype at the SSC4 QTL. Gene ontology (GO) annotation terms for genes located near moderately associated SNP ( < 0.003) were enriched for multiple immunologically (VL) and metabolism- (WG) related GO terms. The biological relevance of these regions suggests that, although it may increase the number of false positives, the use of single-SNP analyses and a relaxed threshold also increased the identification of true positives. In conclusion, although only the SSC4 QTL was associated with response to both PRRSV isolates, genes near associated SNP were enriched for the same GO terms across PRRSV isolates, suggesting that host responses to these 2 isolates are affected by the actions of many genes that function together in similar biological processes.

Association study between single nucleotide polymorphisms in porcine genes and pork quality traits for fresh consumption and processing into Italian dry-cured ham

Single nucleotide polymorphisms (SNPs) of six genes (TTN, PRKAG3, CAST, CTSB, CTSF, and MYPN), known for associations with carcass and meat quality traits, post mortem proteolysis, were screened in a commercial crossed population of 368 heavy pigs (Large White x Landrace)×Duroc, reared according to the rules of Italian Protected Designation of Origin (PDO) dry-cured ham. Carcass, longissimus thoracis et lumborum muscle (LTL), and green ham traits were obtained after slaughtering, main weight losses of dry-cured hams were collected during processing. The results showed the impact of CAST variants on carcass weight, of CTSF on LTL tenderness, ham weight and fatness, of PRKAG3 and TTN on ultimate pH, hamweight. This study, while confirming significant associations between SNPs of genes and qualitative traits of carcass, longissimus and ham, supports CTSF as candidate gene suitable for fresh consumption purpose (tenderness of longissimus at 24h post mortem), and for dry-cured ham processing (higher thickness of ham subcutaneous fat).

Identification of growth trait related genes in a Yorkshire purebred pig population by genome-wide association studies

Objective

The aim of this study is to identify genomic regions or genes controlling growth traits in pigs.

Methods

Using a panel of 54,148 single nucleotide polymorphisms (SNPs), we performed a genome-wide Association (GWA) study in 562 pure Yorshire pigs with four growth traits: average daily gain from 30 kg to 100 kg or 115 kg, and days to 100 kg or 115 kg. Fixed and random model Circulating Probability Unification method was used to identify the associations between 54,148 SNPs and these four traits. SNP annotations were performed through the Sus scrofa data set from Ensembl. Bioinformatics analysis, including gene ontology analysis, pathway analysis and network analysis, was used to identify the candidate genes.

Results

We detected 6 significant and 12 suggestive SNPs, and identified 9 candidate genes in close proximity to them (suppressor of glucose by autophagy [SOGA1], R-Spondin 2 [RSPO2], mitogen activated protein kinase kinase 6 [MAP2K6], phospholipase C beta 1 [PLCB1], rho GTPASE activating protein 24 [ARHGAP24], cytoplasmic polyadenylation element binding protein 4 [CPEB4], GLI family zinc finger 2 [GLI2], neuronal tyrosine-phosphorylated phosphoinositide-3-kinase adaptor 2 [NYAP2], and zinc finger protein multitype 2 [ZFPM2]). Gene ontology analysis and literature mining indicated that the candidate genes are involved in bone, muscle, fat, and lung development. Pathway analysis revealed that PLCB1 and MAP2K6 participate in the gonadotropin signaling pathway and suggests that these two genes contribute to growth at the onset of puberty.

Conclusion

Our results provide new clues for understanding the genetic mechanisms underlying growth traits, and may help improve these traits in future breeding programs.

Integrative Analysis of Metabolomic, Proteomic and Genomic Data to Reveal Functional Pathways and Candidate Genes for Drip Loss in Pigs

The aim of this study was to integrate multi omics data to characterize underlying functional pathways and candidate genes for drip loss in pigs. The consideration of different omics levels allows elucidating the black box of phenotype expression. Metabolite and protein profiling was applied in Musculus longissimus dorsi samples of 97 Duroc × Pietrain pigs. In total, 126 and 35 annotated metabolites and proteins were quantified, respectively. In addition, all animals were genotyped with the porcine 60 k Illumina beadchip. An enrichment analysis resulted in 10 pathways, amongst others, sphingolipid metabolism and glycolysis/gluconeogenesis, with significant influence on drip loss. Drip loss and 22 metabolic components were analyzed as intermediate phenotypes within a genome-wide association study (GWAS). We detected significantly associated genetic markers and candidate genes for drip loss and for most of the metabolic components. On chromosome 18, a region with promising candidate genes was identified based on SNPs associated with drip loss, the protein "phosphoglycerate mutase 2" and the metabolite glycine. We hypothesize that association studies based on intermediate phenotypes are able to provide comprehensive insights in the genetic variation of genes directly involved in the metabolism of performance traits. In this way, the analyses contribute to identify reliable candidate genes.

Genome Wide Association Studies (GWAS) Identify QTL on SSC2 and SSC17 Affecting Loin Peak Shear Force in Crossbred Commercial Pigs

Of all the meat quality traits, tenderness is considered the most important with regard to eating quality and market value. In this study we have utilised genome wide association studies (GWAS) for peak shear force (PSF) of loin muscle as a measure of tenderness for 1,976 crossbred commercial pigs, genotyped for 42,721 informative SNPs using the Illumina PorcineSNP60 Beadchip. Four 1 Mb genomic regions, three on SSC2 (at 4 Mb, 5 Mb and 109 Mb) and one on SSC17 (at 20 Mb), were detected which collectively explained about 15.30% and 3.07% of the total genetic and phenotypic variance for PSF respectively. Markers ASGA0008566, ASGA0008695, DRGA0003285 and ASGA0075615 in the four regions were strongly associated with the effects. Analysis of the reference genome sequence in the region with the most important SNPs for SSC2_5 identified FRMD8, SLC25A45 and LTBP3 as potential candidate genes for meat tenderness on the basis of functional annotation of these genes. The region SSC2_109 was close to a previously reported candidate gene CAST; however, the very weak LD between DRGA0003285 (the best marker representing region SSC2_109) and CAST indicated the potential for additional genes which are distinct from, or interact with, CAST to affect meat tenderness. Limited information of known genes in regions SSC2_109 and SSC17_20 restricts further analysis. Re-sequencing of these regions for informative animals may help to resolve the molecular architecture and identify new candidate genes and causative mutations affecting this trait. These findings contribute significantly to our knowledge of the genomic regions affecting pork shear force and will potentially lead to new insights into the molecular mechanisms regulating meat tenderness.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: II. Genomewide association

Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

BACKGROUND

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

METHODS

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64,432 SNPs on the chip, 44,412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

RESULTS

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

CONCLUSIONS

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.

Investigation of four candidate genes (IGF2, JHDM1A, COPB1 and TEF1) for growth rate and backfat thickness traits on SSC2q in Large White pigs

As important quantitative traits, the growth rate and backfat thickness are controlled by multiple genes. The aim of this investigation was to evaluate the effect of the single and multiple SNPs of four candidate genes (IGF2, JHDM1A, COPB1 and TEF-1) on growth rate and backfat thickness. The four candidate genes were mapped on the p arm of SSC 2, and there are several QTLs, such as average daily gain, backfat thickness, an imprinted QTLs affecting muscle mass and fat deposition have been reported in this region. The polymorphisms of these genes were detected using PCR-RFLP methods, mixed procedure was used to analyze the single marker association with the growth and backfat thickness traits, and the gene-gene combination was investigated using multiple-markers analysis. The single marker association analysis indicated that the IGF2 intron-3 g.3072G > A and the substitution g.93G > A of TEF-1 gene were significantly associated with the age at 100 kg (P < 0.05). The JHDM1A 3′UTR g.224C > G, the c.3096C > T polymorphism of COPB1 gene and the substitution g.93G > A of TEF-1 gene were all significantly associated with the backfat at the shoulder (P < 0.05), backfat at the last rib, backfat at the lumbar, and the average backfat thickness, respectively. The multiple-markers analysis indicated that IGF2 and TEF-1 integrated gene networks for the age at 100 kg. Therefore, we can suggest that the polymorphism of IGF2 and TEF-1 gene could be used in marker-assisted selection for the age at 100 kg in Large White pigs.

Genome-wide copy number variations inferred from SNP genotyping arrays using a Large White and Minzhu intercross population

Copy number variations (CNVs) are one of the main contributors to genetic diversity in animals and are broadly distributed in the genomes of swine. Investigating the performance and evolutionary impacts of pig CNVs requires comprehensive knowledge of their structure and function within and between breeds. In the current study, 4 different programs (i.e., GADA, PennCNV, QuantiSNP, and cnvPartition) were used to analyze Porcine SNP60 genotyping data of 585 pigs from one Large White × Minzhu intercross population to detect copy number variant regions (CNVRs). Overlapping CNVRs recalled by at least 2 programs were used to construct a powerful and comprehensive CNVR map, which contained 249 CNVRs (i.e., 70 gains, 43 losses, and 136 gains/losses) and covered 26.22% of the regions in the swine genome. Ten CNVRs, representing different predicted statuses, were selected for validation via quantitative real-time PCR (QPCR); 9/10 CNVRs (i.e., 90%) were validated. When being traced back to the F0 generation, 58 events were identified in only Minzhu F0 parents and 2 events were identified in only Large White F0 parents. A series of CNVR function analyses were performed. Some of the CNVRs functions were predicted, and several interesting CNVRs for meat quality traits and hematological parameters were obtained. A comprehensive and lower false rate genome-wide CNV map was constructed for Large White and Minzhu pig genomes in this study. Our results may provide an important basis for determining the relationship between CNVRs and important qualitative and quantitative traits. In addition, it can help to further understand genetic processes in pigs.

A genetical genomics approach reveals new candidates and confirms known candidate genes for drip loss in a porcine resource population

In this study lean meat water-holding capacity (WHC) of a Duroc × Pietrain (DuPi) resource population with corresponding genotypes and transcriptomes was investigated using genetical genomics. WHC was characterized by drip loss measured in M. longissimus dorsi. The 60K Illumina SNP chips identified genotypes of 169 F2 DuPi animals. Whole-genome transcriptomes of muscle samples were available for 132 F2 animals using the Affymetrix 24K GeneChip® Porcine Genome Array. Performing genome-wide association studies of transcriptional profiles, which are correlated with phenotypes, allows elucidation of cis- and trans-regulation. Expression levels of 1,228 genes were significantly correlated with drip loss and were further analyzed for enrichment of functional annotation groups as defined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. A hypergeometric gene set enrichment test was performed and revealed glycolysis/glyconeogenesis, pentose phosphate pathway, and pyruvate metabolism as the most promising pathways. For 267 selected transcripts, expression quantitative trait loci (eQTL) analysis was performed and revealed a total of 1,541 significant associations. Because of positional accordance of the gene underlying transcript and the eQTL location, it was possible to identify eight eQTL that can be assumed to be cis-regulated. Comparing the results of gene set enrichment and the eQTL detection tests, molecular networks and potential candidate genes, which seemed to play key roles in the expression of WHC, were detected. The α-1-microglobulin/bikunin precursor (AMBP) gene was assumed to be cis-regulated and was part of the glycolysis pathway. This approach supports the identification of trait-associated SNPs and the further biological understanding of complex traits.

Association mapping of quantitative trait loci for carcass and meat quality traits at the central part of chromosome 2 in Italian Large White pigs

Association mapping of the central part of porcine chromosome 2 harboring QTLs for carcass and meat quality traits was performed with 17 gene-tagged SNPs located between 44.0 and 77.5 Mb on a physical map (Sscrofa10.2) in Italian Large White pigs. For the analyzed animals records of estimated breeding values for average daily gain, back fat thickness, lean cuts, ham weight, feed conversion ratio, pH1, pHu, CIE L*, CIE a*, CIE b* and drip loss were available. A significant QTL for fat deposition (adjusted P=0.0081) and pH1 (adjusted P=0.0972) to MYOD1 at position 44.4 Mb and a QTL for growth and meatiness (adjusted P=0.0238-0.0601) to UBL5 at position 68.9 Mb were mapped. These results from association mapping are much more accurate than those from linkage mapping and facilitate further search for position candidate genes and causative mutations needed for application of markers through marker assisted selection.

Association and expression analysis of porcine HNF1A gene related to meat and carcass quality traits

The aim of this study was to investigate the association and expression of HNF1A gene as a candidate gene for meat and carcass quality traits in pigs. Statistical analysis revealed that the g.8260 A>G polymorphism significantly associated with pH 24(H), meat percentage and muscle area in the F₂ Duroc × Pietrain (DuPi, n=313) and with pH 24(L), fat area and backfat thickness in the Pietrain (Pi, n=110) population. HNF1A mRNA and protein expressions were higher (p<0.05) in animals with the low post-mortem muscle pH 24(L). The promoter methylation profiling suggested that methylation was not involved on HNF1A expression regulation (p>0.05) in animal with divergent muscle pH. In conclusion, polymorphism in porcine HNF1A gene could be used as a candidate marker to improve the meat and carcass quality traits, with the consideration of breed-specific effect.

A consensus linkage map of common carp (Cyprinus carpio L.) to compare the distribution and variation of QTLs associated with growth traits

The ability to detect and identify quantitative trait loci (QTLs) in a single population is often limited. Analyzing multiple populations in QTL analysis improves the power of detecting QTLs and provides a better understanding of their functional allelic variation and distribution. In this study, a consensus map of the common carp was constructed, based on four populations, to compare the distribution and variation of QTLs. The consensus map spans 2371.6 cM across the 42 linkage groups and comprises 257 microsatellites and 421 SNPs, with a mean marker interval of 3.7 cM/marker. Sixty-seven QTLs affecting four growth traits from the four populations were mapped to the consensus map. Only one QTL was common to three populations, and nine QTLs were detected in two populations. However, no QTL was common to all four populations. The results of the QTL comparison suggest that the QTLs are responsible for the phenotypic variability observed for these traits in a broad array of common carp germplasms. The study also reveals the different genetic performances between major and minor genes in different populations.

Mutation in porcine Zip4-like zinc transporter is associated with pancreatic zinc concentration and apparent zinc absorption

The aim of the present study was to analyse the sequence variability of the porcine Zip4-like Zn transporter gene and the association of identified sequence variants with average daily gain, apparent Zn absorption, plasma Zn concentration and Zn concentration in the liver and pancreas. For the purpose of the study, two different sample sets were used. Set one, which was used for sequencing and association analysis, included mRNA from intestinal tissue from thirty-five piglets of a feeding trial. Sample set two consisted of forty-six samples of genomic DNA from sperm or tissue of wild boars and several pig breeds and was used to genotype animals of different breeds. The sequence analysis of porcine Zip4-like complementary DNA in sample set one revealed the presence of seven nucleotide substitutions. Of these, six were synonymous, whereas a substitution of A with C in exon IX (XM_001925360 c.1430A>C) causes an amino acid exchange from glutamic acid to alanine (p.Glu477Ala). The association analysis revealed no influence of the six synonymous substitutions on Zn values, but the non-synonymous nucleotide exchange significantly increased Zn concentration in the pancreas and apparent Zn absorption of the piglets in week 2 of the feeding trial. The parentage of the piglets and the genotyping results in sample set two suggest a breed-specific presence of the A allele in Piétrain for this amino acid substitution. These results indicate that genotype influences the Zn absorption abilities of individual animals, which should be taken into consideration in animal breeding as well as for the selection of experimental animals.

Association of the porcine transforming growth factor beta type I receptor (TGFBR1) gene with growth and carcass traits

BACKGROUND

Growth and carcass traits are of great economic importance in livestock production. A large number of quantitative trait loci (QTL) have been identified for growth and carcass traits on porcine chromosome one (SSC1). A key positional candidate for this chromosomal region is TGFBR1 (transforming growth factor beta type I receptor). This gene plays a key role in inherited disorders at cardiovascular, craniofacial, neurocognitive, and skeletal development in mammals.

RESULTS

In this study, 27 polymorphic SNPs in the porcine TGFBR1 gene were identified on the University of Illinois Yorkshire × Meishan resource population. Three SNPs (SNP3, SNP43, SNP64) representing major polymorphic patterns of the 27 SNPs in F1 and F0 individuals of the Illinois population were selected for analyses of QTL association and genetic diversity. An association analysis for growth and carcass traits was completed using these three representative SNPs in the Illinois population with 298 F2 individuals and a large commercial population of 1008 animals. The results indicate that the TGFBR1 gene polymorphism (SNP64) is significantly associated (p < 0.05) with growth rates including average daily gains between birth and 56 kg (p = 0.049), between 5.5 and 56 kg (p = 0.024), between 35 and 56 kg (p = 0.021). Significant associations (p < 0.05) were also identified between TGFBR1 gene polymorphisms (SNP3/SNP43) and carcass traits including loin-eye-area (p = 0.022) in the Illinois population, and back-fat thickness (p = 0.0009), lean percentage (p = 0.0023) and muscle color (p = 0.021) in the commercial population. These three SNPs were also used to genotype a diverse panel of 130 animals representing 11 pig breeds. Alleles SNP3_T and SNP43_G were fixed in Pietrain and Sinclair pig breeds. SNP64_G allele was uniquely identified in Chinese Meishan pigs. Strong evidence of association (p < 0.01) between both SNP3 and SNP64 alleles and reproductive traits including gestation length and number of corpora lutea were also observed in the Illinois population.

CONCLUSION

This study gives the first evidence of association between the porcine TGFBR1 gene and traits of economic importance and provides support for using TGFBR1 markers for pig breeding and selection programs. The genetic diversities in different pig breeds would be helpful to understand the genetic background and migration of the porcine TGFBR1 gene.

Genome-wide linkage analysis of QTL for growth and body composition employing the PorcineSNP60 BeadChip

BACKGROUND

The traditional strategy to map QTL is to use linkage analysis employing a limited number of markers. These analyses report wide QTL confidence intervals, making very difficult to identify the gene and polymorphisms underlying the QTL effects. The arrival of genome-wide panels of SNPs makes available thousands of markers increasing the information content and therefore the likelihood of detecting and fine mapping QTL regions. The aims of the current study are to confirm previous QTL regions for growth and body composition traits in different generations of an Iberian x Landrace intercross (IBMAP) and especially identify new ones with narrow confidence intervals by employing the PorcineSNP60 BeadChip in linkage analyses.

RESULTS

Three generations (F3, Backcross 1 and Backcross 2) of the IBMAP and their related animals were genotyped with PorcineSNP60 BeadChip. A total of 8,417 SNPs equidistantly distributed across autosomes were selected after filtering by quality, position and frequency to perform the QTL scan. The joint and separate analyses of the different IBMAP generations allowed confirming QTL regions previously identified in chromosomes 4 and 6 as well as new ones mainly for backfat thickness in chromosomes 4, 5, 11, 14 and 17 and shoulder weight in chromosomes 1, 2, 9 and 13; and many other to the chromosome-wide signification level. In addition, most of the detected QTLs displayed narrow confidence intervals, making easier the selection of positional candidate genes.

CONCLUSIONS

The use of higher density of markers has allowed to confirm results obtained in previous QTL scans carried out with microsatellites. Moreover several new QTL regions have been now identified in regions probably not covered by markers in previous scans, most of these QTLs displayed narrow confidence intervals. Finally, prominent putative biological and positional candidate genes underlying those QTL effects are listed based on recent porcine genome annotation.

Genetic determinants for intramuscular fat content and water-holding capacity in mice selected for high muscle mass

Intramuscular fat content and water-holding capacity are important traits in livestock as they influence meat quality, nutritive value of the muscle, and animal health. As a model for livestock, two inbred lines of the Berlin Muscle Mouse population, which had been long-term selected for high muscle mass, were used to identify genomic regions affecting intramuscular fat content and water-holding capacity. The intramuscular fat content of the Musculus longissimus was on average 1.4 times higher in BMMI806 than in BMMI816 mice. This was accompanied by a 1.5 times lower water-holding capacity of the Musculus quadriceps in BMMI816 mice. Linkage analyses with 332 G(3) animals of reciprocal crosses between these two lines revealed quantitative trait loci for intramuscular fat content on chromosome 7 and for water-holding capacity on chromosome 2. In part, the identified loci coincide with syntenic regions in pigs in which genetic effects for the same traits were found. Therefore, these muscle-weight-selected mouse lines and the produced intercross populations are valuable genetic resources to identify genes that could also contribute to meat quality in other species.

Detection of quantitative trait loci affecting serum cholesterol, LDL, HDL, and triglyceride in pigs

Background

Serum lipids are associated with many serious cardiovascular diseases and obesity problems. Many quantitative trait loci (QTL) have been reported in the pig mostly for performance traits but very few for the serum lipid traits. In contrast, remarkable numbers of QTL are mapped for serum lipids in humans and mice. Therefore, the objective of this research was to investigate the chromosomal regions influencing the serum level of the total cholesterol (CT), triglyceride (TG), high density protein cholesterol (HDL) and low density protein cholesterol (LDL) in pigs. For this purpose, a total of 330 animals from a Duroc × Pietrain F2 resource population were phenotyped for serum lipids using ELISA and were genotyped by using 122 microsatellite markers covering all porcine autosomes for QTL study in QTL Express. Blood sampling was performed at approximately 175 days before slaughter of the pig.

Results

Most of the traits were correlated with each other and were influenced by average daily gain, slaughter date and age. A total of 18 QTL including three QTL with imprinting effect were identified on 11 different porcine autosomes. Most of the QTL reached to 5% chromosome-wide (CW) level significance including a QTL at 5% experiment-wide (GW) and a QTL at 1% GW level significance. Of these QTL four were identified for both the CT and LDL and two QTL were identified for both the TG and LDL. Moreover, three chromosomal regions were detected for the HDL/LDL ratio in this study. One QTL for HDL on SSC2 and two QTL for TG on SSC11 and 17 were detected with imprinting effect. The highly significant QTL (1% GW) was detected for LDL at 82 cM on SSC1, whereas significant QTL (5% GW) was identified for HDL/LDL on SSC1 at 87 cM. Chromosomal regions with pleiotropic effects were detected for correlated traits on SSC1, 7 and 12. Most of the QTL identified for serum lipid traits correspond with the previously reported QTL for similar traits in other mammals. Two novel QTL on SSC16 for HDL and HDL/LDL ratio and an imprinted QTL on SSS17 for TG were detected in the pig for the first time.

Conclusion

The newly identified QTL are potentially involved in lipid metabolism. The results of this work shed new light on the genetic background of serum lipid concentrations and these findings will be helpful to identify candidate genes in these QTL regions related to lipid metabolism and serum lipid concentrations in pigs.

Genes with expression levels correlating to drip loss prove association of their polymorphism with water holding capacity of pork

Six genes that were known to exhibit expression levels that are correlated to drip loss BVES, SLC3A2, ZDHHC5, CS, COQ9, and EGFR have been for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing, and used for genotyping. The SNPs were genotyped in about 1,800 animals from six pig populations including commercial herds of Pietrain (PI) and German Landrace (DL), different commercial herds of Pietrain×(German Large White×German Landrace) (PIF1(a/b/c)), and one experimental F2-population Duroc×Pietrain (DUPI). Comparative and genetic mapping established the location of BVES on SSC1, of SLC3A2 and ZDHHC5 on SSC2, of CS on SSC5, of COQ9 on SSC6 and of EGFR on SSC9, respectively, coinciding with QTL regions for carcass and meat quality traits. BVES, SLC3A2, and CS revealed association at least with drip loss and with several other measures of water holding capacity (WHC). Moreover, COQ9 and EGFR were associated with several meat quality traits such as meat color and/or thawing loss. This study reveals statistic evidence in addition to the functional relationship of these genes to WHC previously evidenced by expression analysis. This study reveals positional and genetic statistical evidence for a link of genetic variation at these loci or close to them and promotes those six candidate genes as functional and/or positional candidate genes for meat quality traits.

Identification of Carcass and Meat Quality QTL in an F(2) Duroc × Pietrain Pig Resource Population Using Different Least-Squares Analysis Models

A three-generation resource population was constructed by crossing pigs from the Duroc and Pietrain breeds. In this study, 954 F₂ animals were used to identify quantitative trait loci (QTL) affecting carcass and meat quality traits. Based on results of the first scan analyzed with a line-cross (LC) model using 124 microsatellite markers and 510 F₂ animals, 9 chromosomes were selected for genotyping of additional markers. Twenty additional markers were genotyped for 954 F₂ animals and 20 markers used in the first scan were genotyped for 444 additional F₂ animals. Three different Mendelian models using least-squares for QTL analysis were applied for the second scan: a LC model, a half-sib (HS) model, and a combined LC and HS model. Significance thresholds were determined by false discovery rate (FDR). In total, 50 QTL using the LC model, 38 QTL using the HS model, and 3 additional QTL using the combined LC and HS model were identified (q < 0.05). The LC and HS models revealed strong evidence for QTL regions on SSC6 for carcass traits (e.g., 10th-rib backfat; q < 0.0001) and on SSC15 for meat quality traits (e.g., tenderness, color, pH; q < 0.01), respectively. QTL for pH (SSC3), dressing percent (SSC7), marbling score and moisture percent (SSC12), CIE a* (SSC16), and carcass length and spareribs weight (SSC18) were also significant (q < 0.01). Additional marker and animal genotypes increased the statistical power for QTL detection, and applying different analysis models allowed confirmation of QTL and detection of new QTL.

Genome-wide mapping and identification of new quantitative trait loci affecting meat production, meat quality, and carcass traits within a Duroc purebred population

Most QTL detection studies in pigs have been carried out in experimental F(2) populations. However, segregation of a QTL must be confirmed within a purebred population for successful implementation of marker-assisted selection. Previously, QTL for meat quality and carcass traits were detected on SSC 7 in a Duroc purebred population. The objectives of the present study were to carry out a whole-genome QTL analysis (except for SSC 7) for meat production, meat quality, and carcass traits and to confirm the presence of segregating QTL in a Duroc purebred population. One thousand and four Duroc pigs were studied from base to seventh generation; the pigs comprised 1 closed population of a complex multigenerational pedigree such that all individuals were related. The pigs were evaluated for 6 growth traits, 7 body size traits, 8 carcass traits, 2 physiological traits, and 11 meat quality traits, and the number of pigs with phenotypes ranged from 421 to 953. A total of 119 markers were genotyped and then used for QTL analysis. We utilized a pedigree-based, multipoint variance components approach to test for linkage between QTL and the phenotypic values using a maximum likelihood method; the logarithm of odds score and QTL genotypic heritability were estimated. A total of 42 QTL with suggestive linkages and 3 QTL with significant linkages for 26 traits were detected. These included selection traits such as daily BW gain, backfat thickness, loin eye muscle area, and intramuscular fat content as well as correlated traits such as body size and meat quality traits. The present study disclosed QTL affecting growth, body size, and carcass, physiological, and meat quality traits in a Duroc purebred population.

Application of alternative models to identify QTL for growth traits in an F2 Duroc x Pietrain pig resource population

Background

A variety of analysis approaches have been applied to detect quantitative trait loci (QTL) in experimental populations. The initial genome scan of our Duroc x Pietrain F₂ resource population included 510 F₂ animals genotyped with 124 microsatellite markers and analyzed using a line-cross model. For the second scan, 20 additional markers on 9 chromosomes were genotyped for 954 F₂ animals and 20 markers used in the first scan were genotyped for 444 additional F₂ animals. Three least-squares Mendelian models for QTL analysis were applied for the second scan: a line-cross model, a half-sib model, and a combined line-cross and half-sib model.

Results

In total, 26 QTL using the line-cross model, 12 QTL using the half-sib model and 3 additional QTL using the combined line-cross and half-sib model were detected for growth traits with a 5% false discovery rate (FDR) significance level. In the line-cross analysis, highly significant QTL for fat deposition at 10-, 13-, 16-, 19-, and 22-wk of age were detected on SSC6. In the half-sib analysis, a QTL for loin muscle area at 19-wk of age was detected on SSC7 and QTL for 10th-rib backfat at 19- and 22-wk of age were detected on SSC15.

Conclusions

Additional markers and animals contributed to reduce the confidence intervals and increase the test statistics for QTL detection. Different models allowed detection of new QTL which indicated differing frequencies for alternative alleles in parental breeds.

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.

Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits

BACKGROUND

Genetic analysis of transcriptional profiles is a promising approach for identifying and dissecting the genetics of complex traits like meat performance. Accordingly, expression levels obtained by microarray analysis were taken as phenotypes in a linkage analysis to map eQTL. Moreover, expression levels were correlated with traits related to meat quality and principle components with high loadings of these traits. By using an up-to-date annotation and localization of the respective probe-sets, the integration of eQTL mapping data and information of trait correlated expression finally served to point to candidate genes for meat quality traits.

RESULTS

Genome-wide transcriptional profiles of M. longissimus dorsi RNAs samples of 74 F2 animals of a pig resource population revealed 11,457 probe-sets representing genes expressed in the muscle. Linkage analysis of expression levels of these probe-sets provided 9,180 eQTL at the suggestive significance threshold of LOD > 2. We mapped 653 eQTL on the same chromosome as the corresponding gene and these were designated as 'putative cis-eQTL'. In order to link eQTL to the traits of interest, probe-sets were addressed with relative transcript abundances that showed correlation with meat quality traits at p ≤ 0.05. Out of the 653 'putative cis-eQTL', 262 transcripts were correlated with at least one meat quality trait. Furthermore, association of expression levels with composite traits with high loadings for meat quality traits generated by principle component analysis were taken into account leading to a list of 85 genes exhibiting cis-eQTL and trait dependent expression.

CONCLUSION

Holistic expression profiling was integrated with QTL analysis for meat quality traits. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study.

Transcriptome analysis to identify differential gene expression affecting meat quality in heavy Italian pigs

Suppressive subtractive hybridization (SSH) was used to analyse the muscle transcriptome and identify genes affecting meat quality within an Italian pig population of Large White and Landrace purebred individuals. Seven phenotypes were recorded at slaughter: dorsal fat thickness, ham fat thickness, ham fat coverage, muscle compactness, marbling, meat colour and colour uniformity. Two subtractive libraries were created from longissimus dorsi tissue of selected pigs with extreme phenotypes for meat quality. Eighty-four differentially expressed ESTs were identified, which showed homology to expressed pig sequences and/or to genomic pig sequences produced within the pig genome project. Sixty-eight sequences were mapped on the pig genome, and most of these sequences co-localized with the same chromosomal positions as QTLs that have been previously identified for meat quality. Thirty sequences, including eight matching known genes previously related to muscle metabolic pathways, were selected to statistically validate their differential expression. Association analysis and t-test results indicated that 28 ESTs of the 30 analysed were associated with phenotypes investigated here and have significant differential expression levels (P≤ 0.05) between the two tails of the phenotypic distribution.

Mapping of quantitative trait loci for mycoplasma and tetanus antibodies and interferon-gamma in a porcine F(2) Duroc x Pietrain resource population

The aim of the present study was to detect quantitative trait loci (QTL) for innate and adaptive immunity in pigs. For this purpose, a Duroc x Pietrain F(2) resource population (DUPI) with 319 offspring was used to map QTL for the immune traits blood antibodies and interferon-gamma using 122 microsatellites covering all autosomes. Antibodies response to Mycoplasma hyopneumoniae and tetanus toxoid vaccine and the interferon-gamma (IFNG) serum concentration were measured at three different time points and were used as phenotypes. The differences of antibodies and interferon concentration between different time points were also used for the linkage mapping. Line-cross and imprinting QTL analysis, including two-QTL, were performed using QTL Express. A total of 30 QTL (12, 6, and 12 for mycoplasma, tetanus antibody, and IFNG, respectively) were identified at the 5% chromosome-wide-level significant, of which 28 were detected by line-cross and 2 by imprinting model. In addition, two QTL were identified on chromosome 5 using the two-QTL approach where both loci were in repulsion phase. Most QTL were detected on pig chromosomes 2, 5, 11, and 18. Antibodies were increased over time and immune traits were found to be affected by sex, litter size, parity, and month of birth. The results demonstrated that antibody and IFNG concentration are influenced by multiple chromosomal areas. The flanking markers of the QTL identified for IFNG on SSC5 did incorporate the position of the porcine IFNG gene. The detected QTL will allow further research in these QTL regions for candidate genes and their utilization in selection to improve the immune response and disease resistance in pig.

The insulin-like growth factor 2 (IGF2) gene intron3-g.3072G>A polymorphism is not the only Sus scrofa chromosome 2p mutation affecting meat production and carcass traits in pigs: evidence from the effects of a cathepsin D (CTSD) gene polymorphism

The objective of this study was to evaluate the effects of mutations in 2 genes [IGF2 and cathepsin D (CTSD)] that map on the telomeric end of the p arm of SSC2. In this region, an imprinted QTL affecting muscle mass and fat deposition was reported, and the IGF2 intron3-g.3072G>A substitution was identified as the causative mutation. In the same chromosome region, we assigned, by linkage mapping, the CTSD gene, a lysosomal proteinase, for which we previously identified an SNP in the 3'-untranslated region (AM933484, g.70G>A). We have already shown strong effects of this CTSD mutation on several production traits in Italian Large White pigs, suggesting a possible independent role of this marker in fatness and meat deposition in pigs. To evaluate this hypothesis, after having refined the map position of the CTSD gene by radiation hybrid mapping, we analyzed the IGF2 and the CTSD polymorphisms in 270 Italian Large White and 311 Italian Duroc pigs, for which EBV and random residuals from fixed models were calculated for several traits. Different association analyses were carried out to distinguish the effects of the 2 close markers. In the Italian Large White pigs, the results for IGF2 were highly significant for all traits when using either EBV or random residuals (e.g., using EBV: lean cuts, P = 2.2 x 10(-18); ADG, P = 2.6 x 10(-16); backfat thickness, P = 2.2 x 10(-9); feed:gain ratio, P = 2.3 x 10(-9); ham weight, P = 1.5 x 10(-6)). No effect was observed for meat quality traits. The IGF2 intron3-g.3072G>A mutation did not show any association in the Italian Duroc pigs, probably because of the small variability at this polymorphic site for this breed. However, a significant association was evident for the CTSD marker (P < 0.001) with EBV of all carcass and production traits in Italian Duroc pigs (lean content, ADG, backfat thickness, feed:gain ratio) after excluding possible confounding effects of the IGF2 mutation. The effects of the CTSD g.70G>A mutation were also confirmed in a subset of Italian Large White animals carrying the homozygous genotype IGF2 intron3-g.3072GG, and by haplotype analysis between the markers of the 2 considered genes in the complete data set. Overall, these results indicate that the IGF2 intron3-g.3072G>A mutation is not the only polymorphism affecting fatness and muscle deposition on SSC2p. Therefore, the CTSD g.70G>A polymorphism could be used to increase selection efficiency in marker-assisted selection programs that already use the IGF2 mutation. However, for practical applications, because the CTSD gene should not be imprinted (we obtained this information from expression analysis in adult skeletal muscle), the different modes of inheritance of the 2 genes have to be considered.

Functional genomics and genetical genomics approaches towards elucidating networks of genes affecting meat performance in pigs

The benefit of functional genomics is to identify key pathways and functional networks of genes and candidate genes underlying the genetic control of phenotypes. Genetical genomics, i.e. the integration of genetic analysis and expression phenotypes, has the potential to uncover regulatory networks controlling the coordinated expression of genes and to map variation on the level of DNA affecting the mRNA expression. Here we illustrate our own attempts to apply functional genomics and genetical genomics approaches in order to identify functional networks of genes relevant to traits related to meat performance. Expression data of 74 M longissimus dorsi samples obtained using Affymetrix GeneChips were correlated with drip loss and principal components (PCs) with high loadings of meat quality traits. Functional annotation analyses revealed that differences in water holding capacity, early pH decline and ultimate pH were related to the ubiquitin-proteasome system, mitochondrial metabolic pathways and muscle structural aspects. In particular, 1279 genes were correlated with drip loss (P <or= 0.001; q <or= 0.004). Negatively correlated transcripts were enriched in functional categories like extracellular matrix receptor interaction and Ca-signalling. Transcripts with a positive correlation represented oxidative phosphorylation, mitochondrial pathways and transporter activity. A linkage analysis revealed 897 expression QTL (eQTL) with 104 eQTL mapping in QTL regions for water holding capacity including 8 cis eQTL. The reduction of the multi-dimensional data sets of meat performance traits into lower dimensions of PC and the genetical genomics approach of eQTL analysis proved to be appropriate means to detect relevant biological pathways and to experimentally prioritize candidate genes.

Expression quantitative trait loci analysis of genes in porcine muscle by quantitative real-time RT-PCR compared to microarray data

Genetic analysis of transcriptional profiling is a promising approach for identifying biological pathways and dissecting the genetics of complex traits. Here, we report on expression quantitative trait loci (eQTL) that were estimated from the quantitative real-time RT-PCR data of 276 F(2) animals and compared with eQTL identified using 74 microarrays. In total, 13 genes were selected that showed trait-dependent expression in microarray experiments and exhibited 21 eQTL. Real-time RT-PCR and microarray data revealed seven cis eQTL in total, of which one was only detected by real-time RT-PCR, one was only detected by microarray analysis, three were consistently found in overlapping intervals and two were in neighbouring intervals on the same chromosome; whereas no trans eQTL was confirmed. We demonstrate that cis regulation is a stable characteristic of individual transcripts. Consequently, a global microarray eQTL analysis of a limited number of samples can be used for exploring functional and regulatory gene networks and scanning for cis eQTL, whereas the subsequent analysis of a subset of likely cis-regulated genes by real-time RT-PCR in a larger number of samples is relevant to narrow down a QTL region by targeting these positional candidate genes. In fact, when modelling SNPs of six genes as fixed effects in the eQTL analysis, eQTL peaks were shifted downwards, experimentally confirming the impact of the respective polymorphic genes, although these SNPs were not located in the regulatory sequence and these shifts occur as a result of linkage disequilibrium in the F(2) population.

Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2

Principal component analysis of traits related to carcass and meat properties were combined with microarray expression data for the identification of functional networks of genes and biological processes taking place during the conversion of muscle to meat. Principal components (PCs) with high loadings of meat quality traits were derived from phenotypic data of 572 animals of a porcine crossbreed population. Microarray data of 74 M. longissimus dorsi samples were correlated with PC datasets. Lists of significantly correlated genes were analyzed for enrichment of functional annotation groups as defined in the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases as well as the Ingenuity Pathways Analysis library. Ubiquitination, phosphorylation, mitochondrion dysfunction, actin, integrin, platelet-derived growth factor, epidermal growth factor, vascular endothelial growth factor, and Ca signaling pathways are correlated with meat quality. Among the significantly trait-associated genes, CAPZB, ANKRD1, and CTBP2 are promoted as candidate genes for meat quality that provide a link between the highlighted pathways. Knowledge of the relevant biological processes and the differential expression of members of the pathway will provide tools that are predictive for traits related to meat quality and that may also be diagnostic for many muscle defects or damages including muscle atrophy, dystrophy, and hypoxia.