QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Multivariate genome-wide associations for immune traits in two maternal pig lines

BACKGROUND

Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines.

RESULTS

In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified.

CONCLUSIONS

This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary.

Expression profiles of the CD274 and PLEKHH2 gene and association of its polymorphism with hematologic parameters in sheep

CD274 and PLEKHH2 genes have been identified as immune- and multiple diseases-related genes, and have recently garnered significant interest. However, their role in regulating immune functions in sheep remains largely unexplored. In this study, we aimed to investigate the effects of polymorphisms in CD274 and PLEKHH2 on hematologic parameters in 915 sheep. Our results showed that the CD274 and PLEKHH2 genes were most highly expressed in the spleen and tail fat, respectively, as determined by qRT-PCR. We also identified a G to A mutation (g 0.11858 G > A) in the exon 4 region of CD274, and a C to G mutation (g 0.38384 C > G) in the intron 8 region of PLEKH2. Association analysis revealed that CD274 g 0.11858 G > A was significantly associated with RBC, HCT, MCHC, and MCV (P < 0.05), while PLEKHH2 g 0.38384 C > G was significantly associated with HCT, MPV, MCHC, and MCV (P < 0.05). These results suggest that CD274 and PLEKHH2 genes may play a role in regulating blood physiological indicators and could be potential functional candidates for influencing immune traits in sheep breeding programs.

Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

High-quality metagenome-assembled genomes from proximal colonic microbiomes of synbiotic-treated korean native black pigs reveal changes in functional capacity

Synbiotics are feed supplements with the potential to promote health and productivity in pigs partly, through modulation of the intestinal microbiome. Our study used shotgun sequencing and 16S rRNA gene sequencing techniques to characterize the effect of a synbiotic containing three Lactobacillus species and a fructo-oligosaccharide on the proximal colonic microbiome of 4- to 7-month-old Korean native black gilts. With shotgun sequencing we constructed unique metagenome-assembled genomes of gut microbiota in Native Black Pig for the first time, which we then used for downstream analysis. Results showed that synbiotic treatment did not alter microbial diversity and evenness within the proximal colons, but altered composition of some members of the Lactobacillaceae, Enterococcaceae and Streptococcaceae families. Functional analysis of the shotgun sequence data revealed 8 clusters of orthologous groups (COGs) that were differentially represented in the proximal colonic microbiomes of synbiotic-treated Jeju black pigs relative to controls. In conclusion, our results show that administering this synbiotic causes changes in the functional capacity of the proximal colonic microbiome of the Korean native black pig. This study improves our understanding of the potential impact of synbiotics on the colonic microbiome of Korean native black pigs.

GWAS and Post-GWAS High-Resolution Mapping Analyses Identify Strong Novel Candidate Genes Influencing the Fatty Acid Composition of the Longissimus dorsi Muscle in Pigs

Fatty acid (FA) composition is one of the most important parameters for the assessment of meat quality in pigs. The FA composition in pork can also affect human health. Our aim was to identify quantitative trait loci (QTLs) and positional candidate genes affecting the FA profile of the longissimus dorsi muscle in a large F2 intercross between Landrace and Korean native pigs comprising 1105 F2 progeny by genome-wide association studies (GWAS) and post-GWAS high-resolution mapping analyses. We performed GWAS using the PorcineSNP60K BeadChip and a linear mixed model. Four genome-wide significant QTL regions in SSC8, SSC12, SSC14, and SSC16 were detected (p < 2.53 × 10-7). Several co-localizations of QTLs in SSC12 for oleic acid, linoleic acid, arachidonic acid, monounsaturated FAs, polyunsaturated FAs, and the polyunsaturated/saturated FA ratio were observed. To refine the QTL region in SSC12, a linkage and linkage disequilibrium analysis was applied and could narrow down the critical region to a 0.749 Mb region. Of the genes in this region, GAS7, MYH2, and MYH3 were identified as strong novel candidate genes based on further conditional association analyses. These findings provide a novel insight into the genetic basis of FA composition in pork and could contribute to the improvement of pork quality.

Myogenic regulatory factors are key players in determining muscle mass and meat quality in Jeju native and Berkshire pigs

BACKGROUND

Meat from Jeju native pigs (JNPs) is highly popular among Korean consumers; however, the production efficiency is limited due to the low adult body weight. In contrast, the Berkshire breed, which has a genetic background closely related to Asian native pigs, gains weight more efficiently.

OBJECTIVES

This study focused on the differential expression of genes related to muscle growth in postnatal myogenesis between Berkshire and JNPs, specifically the myogenic regulatory factor (MRF) genes (MyoD, Pax7, Myf5, Myf6 and MyBPH). The MRF family is primarily involved in the proliferation and development of muscle.

METHODS

Qualitative reverse transcription-polymerase chain reaction and western blot analyses revealed that expression of MyoD and Pax7 was significantly higher in Berkshire pigs than in JNPs. In addition, co-expression of MyoD and Pax7 was observed in myotubes formed in cultured C2C12 cells. ToppCluster was used to elucidate the relationship between biological processes of the MRFs and muscle-related signalling pathways.

RESULTS

MyoD and Pax7 are factors essential for the activation of satellite cell during myogenesis. However, the mRNA and protein levels of MyBPH (which is responsible for meat quality, e.g. water content, colour and tenderness) are significantly higher in both 1-day-old piglets and adult JNPs than in Berkshire pigs.

CONCLUSIONS

This study provides a genetic understanding of myogenesis in the postnatal and adult stages of Berkshire pigs and JNPs. Moreover, these results will help identify marker genes related to muscle mass, growth performance and meat quality in indigenous Korean pig breeds.

Genetic parameters and associated genomic regions for global immunocompetence and other health-related traits in pigs

The inclusion of health-related traits, or functionally associated genetic markers, in pig breeding programs could contribute to produce more robust and disease resistant animals. The aim of the present work was to study the genetic determinism and genomic regions associated to global immunocompetence and health in a Duroc pig population. For this purpose, a set of 30 health-related traits covering immune (mainly innate), haematological, and stress parameters were measured in 432 healthy Duroc piglets aged 8 weeks. Moderate to high heritabilities were obtained for most traits and significant genetic correlations among them were observed. A genome wide association study pointed out 31 significantly associated SNPs at whole-genome level, located in six chromosomal regions on pig chromosomes SSC4, SSC6, SSC17 and SSCX, for IgG, γδ T-cells, C-reactive protein, lymphocytes phagocytic capacity, total number of lymphocytes, mean corpuscular volume and mean corpuscular haemoglobin. A total of 16 promising functionally-related candidate genes, including CRP, NFATC2, PRDX1, SLA, ST3GAL1, and VPS4A, have been proposed to explain the variation of immune and haematological traits. Our results enhance the knowledge of the genetic control of traits related with immunity and support the possibility of applying effective selection programs to improve immunocompetence in pigs.

A functional regulatory variant of MYH3 influences muscle fiber-type composition and intramuscular fat content in pigs

Muscle development and lipid accumulation in muscle critically affect meat quality of livestock. However, the genetic factors underlying myofiber-type specification and intramuscular fat (IMF) accumulation remain to be elucidated. Using two independent intercrosses between Western commercial breeds and Korean native pigs (KNPs) and a joint linkage-linkage disequilibrium analysis, we identified a 488.1-kb region on porcine chromosome 12 that affects both reddish meat color (a*) and IMF. In this critical region, only the MYH3 gene, encoding myosin heavy chain 3, was found to be preferentially overexpressed in the skeletal muscle of KNPs. Subsequently, MYH3-transgenic mice demonstrated that this gene controls both myofiber-type specification and adipogenesis in skeletal muscle. We discovered a structural variant in the promotor/regulatory region of MYH3 for which Q allele carriers exhibited significantly higher values of a* and IMF than q allele carriers. Furthermore, chromatin immunoprecipitation and cotransfection assays showed that the structural variant in the 5'-flanking region of MYH3 abrogated the binding of the myogenic regulatory factors (MYF5, MYOD, MYOG, and MRF4). The allele distribution of MYH3 among pig populations worldwide indicated that the MYH3 Q allele is of Asian origin and likely predates domestication. In conclusion, we identified a functional regulatory sequence variant in porcine MYH3 that provides novel insights into the genetic basis of the regulation of myofiber type ratios and associated changes in IMF in pigs. The MYH3 variant can play an important role in improving pork quality in current breeding programs.

Genome-wide association studies for 30 haematological and blood clinical-biochemical traits in Large White pigs reveal genomic regions affecting intermediate phenotypes

Haematological and clinical-biochemical parameters are considered indicators of the physiological/health status of animals and might serve as intermediate phenotypes to link physiological aspects to production and disease resistance traits. The dissection of the genetic variability affecting these phenotypes might be useful to describe the resilience of the animals and to support the usefulness of the pig as animal model. Here, we analysed 15 haematological and 15 clinical-biochemical traits in 843 Italian Large White pigs, via three genome-wide association scan approaches (single-trait, multi-trait and Bayesian). We identified 52 quantitative trait loci (QTLs) associated with 29 out of 30 analysed blood parameters, with the most significant QTL identified on porcine chromosome 14 for basophil count. Some QTL regions harbour genes that may be the obvious candidates: QTLs for cholesterol parameters identified genes (ADCY8, APOB, ATG5, CDKAL1, PCSK5, PRL and SOX6) that are directly involved in cholesterol metabolism; other QTLs highlighted genes encoding the enzymes being measured [ALT (known also as GPT) and AST (known also as GOT)]. Moreover, the multivariate approach strengthened the association results for several candidate genes. The obtained results can contribute to define new measurable phenotypes that could be applied in breeding programs as proxies for more complex traits.

The first comprehensive description of the expression profile of genes involved in differential body growth and the immune system of the Jeju Native Pig and miniature pig

Sus scrofa provides a major source of animal protein for humans as well as being an excellent biomedical model. This study was carried out to understand, in detail, the genetic and functional variants of Jeju Native Pigs and miniature pigs through differential expression profiling of the genes controlling their immune response, growth performance, and meat quality. The Illumina HiSeq 2000 platform was used for generating 1.3 billion 90 bp paired-end reads, which were mapped to the S. scrofa genome using TopHat2. A total of 2481 and 2768 genes were differentially expressed with 8-log changes in muscle and liver samples, respectively. Five hundred forty-eight genes in muscle and 642 genes in liver samples had BLAST matches within the non-redundant database. GO process and pathway analyses showed enhanced biological processes related to the extracellular structural organization and skeletal muscle cell differentiation in muscle tissue, whereas the liver tissue shares functions related to the inflammatory response. Herein, we identify inflammatory regulatory genes in miniature pigs and growth response genes in Jeju Native Pigs, information which can provide a stronger base for the selection of breeding stock and facilitate further in vitro and in vivo studies for therapeutic purposes.

Rapid Communication: High-resolution quantitative trait loci analysis identifies LTBP2 encoding latent transforming growth factor beta binding protein 2 associated with thoracic vertebrae number in a large F2 intercross between Landrace and Korean native pigs

Number of vertebrae is associated with body size and meat productivity in pigs. The aim of this study was to identify QTL and associated positional candidate genes affecting the number of thoracic vertebrae (THO). A genomewide association study was conducted in a large resource population derived from an F intercross between Landrace and Korean native pigs using the Porcine SNP 60K BeadChip and the genomewide complex trait analysis (GCTA) program based on a linear mixed-effects model. A total of 38,385 SNP markers from 1,105 F progeny were analyzed for the THO trait after filtering for quality control. A total of 90 genomewide significant SNP markers ( < 1.30 × 10) on SSC 7 covering a 20-Mb region were identified for THO in this study. Several previous studies also mapped QTL for vertebral numbers in this region. The strongest association signals were detected at ASGA0035500 (-value = 4.46 × 10; 103,574,383 bp) and DIAS0000795 (-value = 4.46 × 10; 103,594,753 bp). The QTL region on SSC 7 for THO encompasses and , which are previously described candidate genes for vertebral number variation. To refine the QTL region, a haplotype-based linkage and linkage disequilibrium (LALD) analysis using the DualPHASE program was applied because subsequent conditional association and haplotype block analyses could not resolve the region that contains the 2 loci. The LALD analysis refined the critical region to a 533.9-kb region including ; was located outside the critical region. The gene encoding latent transforming growth factor beta binding protein 2 is involved in bone metabolisms. Based on these data, we propose as a positional candidate gene for THO in pigs. After further functional studies and verification of the association in other independent populations, these results could be useful for optimizing breeding programs that improve THO and other economically important traits in pigs.

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.

Whole genome scan reveals the genetic signature of African Ankole cattle breed and potential for higher quality beef

BACKGROUND

Africa is home to numerous cattle breeds whose diversity has been shaped by subtle combinations of human and natural selection. African Sanga cattle are an intermediate type of cattle resulting from interbreeding between Bos taurus and Bos indicus subspecies. Recently, research has asserted the potential of Sanga breeds for commercial beef production with better meat quality as compared to Bos indicus breeds. Here, we identified meat quality related gene regions that are positively selected in Ankole (Sanga) cattle breeds as compared to indicus (Boran, Ogaden, and Kenana) breeds using cross-population (XP-EHH and XP-CLR) statistical methods.

RESULTS

We identified 238 (XP-EHH) and 213 (XP-CLR) positively selected genes, of which 97 were detected from both statistics. Among the genes obtained, we primarily reported those involved in different biological process and pathways associated with meat quality traits. Genes (CAPZB, COL9A2, PDGFRA, MAP3K5, ZNF410, and PKM2) involved in muscle structure and metabolism affect meat tenderness. Genes (PLA2G2A, PARK2, ZNF410, MAP2K3, PLCD3, PLCD1, and ROCK1) related to intramuscular fat (IMF) are involved in adipose metabolism and adipogenesis. MB and SLC48A1 affect meat color. In addition, we identified genes (TIMP2, PKM2, PRKG1, MAP3K5, and ATP8A1) related to feeding efficiency. Among the enriched Gene Ontology Biological Process (GO BP) terms, actin cytoskeleton organization, actin filament-based process, and protein ubiquitination are associated with meat tenderness whereas cellular component organization, negative regulation of actin filament depolymerization and negative regulation of protein complex disassembly are involved in adipocyte regulation. The MAPK pathway is responsible for cell proliferation and plays an important role in hyperplastic growth, which has a positive effect on meat tenderness.

CONCLUSION

Results revealed several candidate genes positively selected in Ankole cattle in relation to meat quality characteristics. The genes identified are involved in muscle structure and metabolism, and adipose metabolism and adipogenesis. These genes help in the understanding of the biological mechanisms controlling beef quality characteristics in African Ankole cattle. These results provide a basis for further research on the genomic characteristics of Ankole and other Sanga cattle breeds for quality beef.

Genome scan linkage analysis identifies a major quantitative trait loci for fatty acid composition in longissimus dorsi muscle in an F2 intercross between Landrace and Korean native pigs

OBJECTIVE

This study was conducted to locate quantitative trait loci (QTL) influencing fatty acid (FA) composition in a large F₂ intercross between Landrace and Korean native pigs.

METHODS

Eighteen FA composition traits were measured in more than 960 F₂ progeny. All experimental animals were genotyped with 165 microsatellite markers located throughout the pig autosomes.

RESULTS

We detected 112 QTLs for the FA composition; Forty seven QTLs reached the genome-wide significant threshold. In particular, we identified a cluster of highly significant QTLs for FA composition on SSC12. QTL for polyunsaturated fatty acid on pig chromosome 12 (F-value = 97.2 under additive and dominance model, nominal p-value 3.6×10^-39) accounted for 16.9% of phenotypic variance. In addition, four more QTLs for C18:1, C18:2, C20:4, and monounsaturated fatty acids on the similar position explained more than 10% of phenotypic variance.

CONCLUSION

Our findings of a major QTL for FA composition presented here could provide helpful information to locate causative variants to improve meat quality traits in pigs.

Association of a missense mutation in the positional candidate gene glutamate receptor-interacting protein 1 with backfat thickness traits in pigs

OBJECTIVE

Previously, we reported quantitative trait loci (QTLs) affecting backfat thickness (BFT) traits on pig chromosome 5 (SW1482-SW963) in an F2 intercross population between Landrace and Korean native pigs. The aim of this study was to evaluate glutamate receptor-interacting protein 1 (GRIP1) as a positional candidate gene underlying the QTL affecting BFT traits.

METHODS

Genotype and phenotype analyses were performed using the 1,105 F₂ progeny. A mixed-effect linear model was used to access association between these single nucleotide polymorphism (SNP) markers and the BFT traits in the F₂ intercross population.

RESULTS

Highly significant associations of two informative SNPs (c.2442 T>C, c.3316 C>G [R1106G]) in GRIP1 with BFT traits were detected. In addition, the two SNPs were used to construct haplotypes that were also highly associated with the BFT traits.

CONCLUSION

The SNPs and haplotypes of the GRIP1 gene determined in this study can contribute to understand the genetic structure of BFT traits in pigs.

Association of the Single Nucleotide Polymorphisms in RUNX1, DYRK1A, and KCNJ15 with Blood Related Traits in Pigs

The aim of this study was to detect positional candidate genes located within the support interval (SI) regions based on the results of red blood cell, mean corpuscular volume (MCV), and mean corpuscular hemoglobin quantitative trait locus (QTL) in Sus scrofa chromosome 13, and to verify the correlation between specific single-nucleotide polymorphisms (SNPs) located in the exonic region of the positional candidate gene and the three genetic traits. The flanking markers of the three QTL SI regions are SW38 and S0215. Within the QTL SI regions, 44 genes were located, and runt-related transcription factor 1, dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 1A (DYRK1A), and potassium inwardly-rectifying channel, subfamily J, member 15 KCNJ15-which are reported to be related to the hematological traits and clinical features of Down syndrome-were selected as positional candidate genes. The ten SNPs located in the exonic region of the three genes were detected by next generation sequencing. A total of 1,232 pigs of an F₂ resource population between Landrace and Korean native pigs were genotyped. To investigate the effects of the three genes on each genotype, a mixed-effect model which is the considering family structure model was used to evaluate the associations between the SNPs and three genetic traits in the F₂ intercross population. Among them, the MCV level was highly significant (nominal p = 9.8×10^-9) in association with the DYRK1A-SNP1 (c.2989 G<A), and explained 4.9% of the phenotypic variance. However, since linkage disequilibrium is quite extensive in an F₂ intercross, our approach has limited power to distinguish one particular positional candidate gene from a QTL region.

Single- and Bayesian Multi-Marker Genome-Wide Association for Haematological Parameters in Pigs

Haematological traits are important traits that show associations with immune and metabolic status, as well as diseases in humans and animals. Mapping genome regions that affect the blood cell traits can contribute to the identification of genomic features useable as biomarkers for immune, disease and metabolic status. A genome-wide association study (GWAS) was conducted using PorcineSNP60 BeadChips. Single-marker and Bayesian multi-marker approaches were integrated to identify genomic regions and corresponding genes overlapping for both methods. GWAS was performed for haematological traits of 591 German Landrace pig. Heritability estimates for haematological traits were medium to high. In total 252 single SNPs associated with 12 haematological traits were identified (NegLog10 of p-value > 5). The Bayesian multi-marker approach revealed 102 QTL regions across the genome, indicated by 1-Mb windows with contribution to additive genetic variance above 0.5%. The integration of both methods resulted in 24 overlapping QTL regions. This study identified overlapping QTL regions from single- and multi-marker approaches for haematological traits. Identifying candidate genes that affect blood cell traits provides the first step towards the understanding of the molecular basis of haematological phenotypes.

An integrated in silico approach for functional and structural impact of non- synonymous SNPs in the MYH1 gene in Jeju Native Pigs

BACKGROUND

This study was performed to identify the non- synonymous polymorphisms in the myosin heavy chain 1 gene (MYH1) association with skeletal muscle development in economically important Jeju Native Pig (JNP) and Berkshire breeds. Herein, we present an in silico analysis, with a focus on (a) in silico approaches to predict the functional effect of non-synonymous SNP (nsSNP) in MYH1 on growth, and (b) molecular docking and dynamic simulation of MYH1 to predict the effects of those nsSNP on protein-protein association.

RESULTS

The NextGENe (V 2.3.4.) tool was used to identify the variants in MYH1 from JNP and Berkshire using RNA seq. Gene ontology analysis of MYH1 revealed significant association with muscle contraction and muscle organ development. The 95 % confidence intervals clearly indicate that the mRNA expression of MYH1 is significantly higher in the Berkshire longissimus dorsi muscle samples than JNP breed. Concordant in silico analysis of MYH1, the open-source software tools identified 4 potential nsSNP (L884T, K972C, N981G, and Q1285C) in JNP and 1 nsSNP (H973G) in Berkshire pigs. Moreover, protein-protein interactions were studied to investigate the effect of MYH1 mutations on association with hub proteins, and MYH1 was found to be closely associated with the protein myosin light chain, phosphorylatable, fast skeletal muscle MYLPF. The results of molecular docking studies on MYH1 (native and 4 mutants) and MYLFP demonstrated that the native complex showed higher electrostatic energy (-466.5 Kcal mol(-1)), van der Walls energy (-87.3 Kcal mol(-1)), and interaction energy (-835.7 Kcal mol(-1)) than the mutant complexes. Furthermore, the molecular dynamic simulation revealed that the native complex yielded a higher root-mean-square deviation (0.2-0.55 nm) and lower root-mean-square fluctuation (approximately 0.08-0.3 nm) as compared to the mutant complexes.

CONCLUSIONS

The results suggest that the variants at L884T, K972C, N981G, and Q1285C in MYH1 in JNP might represent a cause for the poor growth performance for this breed. This study is a pioneering in-depth in silico analysis of polymorphic MYH1 and will serve as a valuable resource for further targeted molecular diagnosis and population-based studies conducted for improving the growth performance of JNP.

Evaluation of body growth and immunity-related differentially expressed genes through deep RNA sequencing in the piglets of Jeju native pig and Berkshire

This study was carried out with the objective to investigate the differentially expressed genes (DEGs) between Jeju native pig (JNP) and Berkshire piglets. The RNA-Seq technique was used to investigate the transcriptomes in the fat, liver and longissimus dorsi muscle from these two breeds. Paired-end reads of the sequences that passed the quality filters were aligned to the Sus scrofa genome using tophat2 (v2.0.2). In this study, 65% of muscle, 20% of fat and 54% of liver genes showed higher expression in the piglets of JNP than in Berkshire. Gene Ontology and signaling pathways showed that immune response and lipid metabolisms were commonly enriched pathways in all three tissues. It was found that the genes pertaining to body growth and immune system are significantly (P < 0.01) more highly expressed in Berkshire piglets. DEGs explored between the piglets of the two breeds might influence the identification of the genetic markers for further breed improvement programs. Our findings provide a new perspective for understanding and identifying candidate genes that are involved in various biological functions. Moreover, transcriptome analysis makes it easier to understand the differences between genetic mechanisms of breeds.

Genome-wide QTL analysis of meat quality-related traits in a large F2 intercross between Landrace and Korean native pigs

Background

We conducted a genome-wide linkage analysis to identify quantitative trait loci (QTL) that influence meat quality-related traits in a large F₂ intercross between Landrace and Korean native pigs. Thirteen meat quality-related traits of the m. longissimus lumborum et thoracis were measured in more than 830 F₂ progeny. All these animals were genotyped with 173 microsatellite markers located throughout the pig genome, and the GridQTL program based on the least squares regression model was used to perform the QTL analysis.

Results

We identified 23 genome-wide significant QTL in eight chromosome regions (SSC1, 2, 6, 7, 9, 12, 13, and 16) (SSC for Sus Scrofa) and detected 51 suggestive QTL in the 17 chromosome regions. QTL that affect 10 meat quality traits were detected on SSC12 and were highly significant at the genome-wide level. In particular, the QTL with the largest effect affected crude fat percentage and explained 22.5% of the phenotypic variance (F-ratio = 278.0 under the additive model, nominal P = 5.5 × 10⁻⁵⁵). Interestingly, the QTL on SSC12 that influenced meat quality traits showed an obvious trend for co-localization.

Conclusions

Our results confirm several previously reported QTL. In addition, we identified novel QTL for meat quality traits, which together with the associated positional candidate genes improve the knowledge on the genetic structure that underlies genetic variation for meat quality traits in pigs.

Electronic supplementary material

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

Genome-wide detection and characterization of positive selection in Korean Native Black Pig from Jeju Island

Background

In the 1980s, Korean native black pigs from Jeju Island (Jeju black pigs) served as representative sample of Korean native black pigs, and efforts were made to help the species rebound from the brink of extinction, which occurred as a result of the introduction of Western pig breeds. Geographical separation of Jeju Island from the Korean peninsula has allowed Jeju black pigs not only to acquire unique characteristics but also to retain merits of rare Korean native black pigs.

Results

To further analyze the Jeju black pig genome, we performed whole-genome re-sequencing (average read depth of 14×) of 8 Jeju black pig and 6 Korean pigs (which live on the Korean peninsula) to compare and identify putative signatures of positive selection in Jeju black pig, the true and pure Korean native black pigs. The candidate genes potentially under positive selection in Jeju black pig support previous reports of high marbling score, rare occurrence of pale, soft, exudative (PSE) meat, but low growth rate and carcass weight compared to Western breeds.

Conclusions

Several candidate genes potentially under positive selection were involved in fatty acid transport and may have contributed to the unique characteristics of meat quality in JBP. Jeju black pigs can offer a unique opportunity to investigate the true genetic resource of once endangered Korean native black pigs. Further genome-wide analyses of Jeju black pigs on a larger population scale are required in order to define a conservation strategy and improvement of native pig resources.

Electronic supplementary material

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

Genome-wide association analysis to identify SNP markers affecting teat numbers in an F2 intercross population between Landrace and Korean native pigs

Most reproductive traits have low heritability and are greatly affected by environmental factors. Teat number and litter size are traits related to the reproduction ability of pigs. To identify quantitative trait loci (QTLs) for teat number traits, a genome-wide association study (GWAS) was conducted using an F2 intercross between Landrace and Korean native pigs. Genotype analysis was performed using the porcine SNP 60 K beadchip. The GWAS was performed using a mixed-effects model and linear regression approach. When a genome-wide threshold was determined using the Bonferroni method (P = 1.61 × 10(-6)), 38 single nucleotide polymorphism (SNP) markers in pig chromosome 7 (SSC7) were significantly associated with three teat number traits (total teat number, left teat number, and right teat number). Among these, SNPs in 5 genes (HDDC3, LOC100156276, LOC100155863, ANPEP, SCAMP2) were selected for further study based primarily on their statistical significance. A significant association was detected in SCAMP2 g.25280 G>A for total teat number (P = 2.0 × 10(-12)), HDDC3 g.1319 G>A SNP for left teat number (P = 2.3 × 10(-7)), and SCAMP2 g.14198 G>A for right teat number (P = 4.7 × 10(-12)). These results provide valuable information about the selective breeding for desirable teat numbers in pigs.

The Possibility of TBC1D21 as a Candidate Gene for Teat Numbers in Pigs

Based on a quantitative traits locus (QTL) study using a F₂ intercross between Landrace and Korean native pigs, a significant QTL affecting teat numbers in SSC7 was identified. The strong positional candidate gene, TBC1D21, was selected due to its biological function for epithelial mesenchymal cell development. Sequence analysis revealed six single nucleotide polymorphisms (SNPs) in the TBC1D21 gene. Among these, two SNP markers, one silent mutation (SNP01) for g.13,050A>G and one missense mutation (SNP04) for c.829A>T (S277C), were genotyped and they showed significant associations with teat number traits (p value = 6.38E-05 for SNP01 and p value = 1.06E-07 for SNP04 with total teat numbers). Further functional validation of these SNPs could give valuable information for understanding the teat number variation in pigs.

A Whole Genome Association Study to Detect Single Nucleotide Polymorphisms for Blood Components (Immunity) in a Cross between Korean Native Pig and Yorkshire

The purpose of this study was to detect significant SNPs for blood components that were related to immunity using high single nucleotide polymorphism (SNP) density panels in a Korean native pig (KNP)×Yorkshire (YK) cross population. A reciprocal design of KNP×YK produced 249 F₂ individuals that were genotyped for a total of 46,865 available SNPs in the Illumina porcine 60K beadchip. To perform whole genome association analysis (WGA), phenotypes were regressed on each SNP under a simple linear regression model after adjustment for sex and slaughter age. To set up a significance threshold, 0.1% point-wise p value from F distribution was used for each SNP test. Among the significant SNPs for a trait, the best set of SNP markers were determined using a stepwise regression procedure with the rates of inclusion and exclusion of each SNP out of the model at 0.001 level. A total of 54 SNPs were detected; 10, 6, 4, 4, 5, 4, 5, 10, and 6 SNPs for neutrophil, lymphocyte, monocyte, eosinophil, basophil, atypical lymph, immunoglobulin, insulin, and insulin-like growth factor-I, respectively. Each set of significant SNPs per trait explained 24 to 42% of phenotypic variance. Several pleiotropic SNPs were detected on SSCs 4, 13, 14 and 15.

Comparative transcriptomic analysis to identify differentially expressed genes in fat tissue of adult Berkshire and Jeju Native Pig using RNA-seq

Pork is a major source of animal protein for humans. The subcutaneous, intermuscular and the intramuscular fat are the factors responsible for meat quality. RNA-seq is rapidly adopted for the profiling of the transcriptomes in the studies related to gene regulation. The discovery of differentially expressed genes (DEGs) between adult animals of Jeju Native Pig (JNP) and Berkshire breeds are of particular interest for the current study. RNA-seq was used to investigate the transcriptome profiling in the fat tissue. Sequence reads were obtained from Ilumina HiSeq2000 and mapped to the pig genome using Tophat2. Total 153 DEGs were identified and 71 among the annotated genes, have BLAST matches in the non- redundant database. Metabolic, immune response and protein binding are enriched pathways in the fat tissue. In our study, biological adhesion, cellular, developmental and multicellular organismal processes in fat were up-regulated in JNP as compare to Berkshire. Multicellular organismal process, developmental process, embryonic morphogenesis and skeletal system development were the most significantly enriched terms in fat of JNP and Berkshire breeds (p = 1.17E-04, 0.044, 3.47E-04 and 4.48E-04 respectively). COL10A1, COL11A2, PDK4 and PNPLA3 genes responsible for skeletal system morphogenesis and body growth were down regulated in JNP. This study is the first statistical analysis for the detection of DEGs from RNA-seq data generated from fat tissue sample. This analysis can be used as stepping stone to understand the difference in the genetic mechanisms that might influence the identification of novel transcripts, sequence polymorphisms, isoforms and noncoding RNAs.

Comparative transcriptomic analysis by RNA-seq to discern differential expression of genes in liver and muscle tissues of adult Berkshire and Jeju Native Pig

RNA-seq is being rapidly adopted for the profiling of the transcriptomes in different areas of biology, especially in the studies related to gene regulation. The discovery of differentially expressed genes (DEGs) between adult animals of Jeju Native Pig (JNP) and Berkshire breeds of Sus scrofa, is of particular interest for the current study. For the better understanding of the gene expression profiles of the liver and longissimus dorsi muscle, DEGs were identified via RNA-seq. Sequence reads were obtained from Illumina HiSeq2000 and mapped to the pig reference genome (Sscrofa10.2) using Tophat2. We identified 169 and 39 DEGs in the liver and muscle of JNP respectively, by comparison with Berkshire breed. Out of all identified genes, 41 genes in the liver and 9 genes in the muscle have given significant expression. Gene ontology (GO) terms of developmental process and KEGG pathway analysis showed that metabolic, immune response and protein binding were commonly enriched pathways in the two tissues. Further the heat map analysis by ArrayStar has shown the different levels of expression in JNP with respect to the Berkshire breed. The validation through real time PCR and western blotting also confirmed the differential expression of genes in both breeds. Genes pertaining to metabolic process and inflammatory and immune system are more enriched in Berkshire breed. This comparative transcriptome analysis of two tissues suggests a subset of novel marker genes which expressed differently between the JNP and Berkshire.

Genome-wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome-wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F2 progeny were analyzed for each blood-related trait after filtering for quality control. Data were analyzed by the genome-wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10(-6) ) on SSC3, 6, 8, 13 and 17 were identified for blood-related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in 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.

Genome-wide association study reveals constant and specific loci for hematological traits at three time stages in a White Duroc × Erhualian F2 resource population

Hematological traits are important indicators of immune function and have been commonly examined as biomarkers of disease and disease severity in humans. Pig is an ideal biomedical model for human diseases due to its high degree of similarity with human physiological characteristics. Here, we conducted genome-wide association studies (GWAS) for 18 hematological traits at three growth stages (days 18, 46 and 240) in a White Duroc × Erhualian F₂ intercross. In total, we identified 38 genome-wide significant regions containing 185 genome-wide significant SNPs by single-marker GWAS or LONG-GWAS. The significant regions are distributed on pig chromosomes (SSC) 1, 4, 5, 7, 8, 10, 11, 12, 13, 17 and 18, and most of significant SNPs reside on SSC7 and SSC8. Of the 38 significant regions, 7 show constant effects on hematological traits across the whole life stages, and 6 regions have time-specific effects on the measured traits at early or late stages. The most prominent locus is the genomic region between 32.36 and 84.49 Mb on SSC8 that is associated with multiple erythroid traits. The KIT gene in this region appears to be a promising candidate gene. The findings improve our understanding of the genetic architecture of hematological traits in pigs. Further investigations are warranted to characterize the responsible gene(s) and causal variant(s) especially for the major loci on SSC7 and SSC8.

Association of the gene encoding stearoyl-CoA desaturase (SCD) with fatty acid composition in an intercross population between Landrace and Korean native pigs

Fatty acid (FA) composition is one of the most important parameters for the evaluation of meat quality. The stearoyl-CoA desaturase (SCD) gene is considered a positional candidate gene affecting FA composition in SSC14, based on previous quantitative trait loci studies. To evaluate the association of the SCD gene with FA composition in a Korean native pig × Landrace F(2) intercross population, we genotyped six single nucleotide polymorphisms (SNPs) of the SCD gene located in promoter region [2 SNPs (g.-353T>C, g.-233T>C)], exonic region [1 SNP (g.817C>T) in exon 2] and 3' UTR [3 SNPs (g.13311C>G, g.14384G>A, and g.14424C>T)] identified by massively parallel sequencing technology. Eighteen FA composition traits were measured in more than the 950 F(2) animals. A mixed-effect model was used to evaluate associations between these SNPs and FA composition traits in the F(2) intercross population. A detailed investigation detected that the five FA composition traits [palmitoleic acid (C16:1), stearic acid (C18:0), arachidic acid (C20:0), saturated FA, and unsaturated FA] were highly significant (P < 4.7 × 10(-5); C20:0) in association with the SNP g.-233T>C, SNP g.817C>T, SNP g.13311C>G and SNP g.14384G>A in the SCD gene, whereas SNP g.14424C>T was only significantly associated with palmitoleic acid (C16:1, P = 1.4 × 10(-3)). No significant association of FA composition traits with SNP g.-353T>C was detected. In particular, the SNP g.14384G>A accounted for 30.6 % of the additive genetic variance of palmitoleic acid (P = 1.9 × 10(-10)). These results suggest the SCD gene has a strong effect on FA composition in the crossbred pig population.

A missense mutation (c.1963A<G) of the complementary component 2 (C2) gene is associated with serum Ca⁺⁺ concentrations in pigs

Serum Ca(++) levels play important roles in the humoral immunity. The aim of this study was to detect quantitative trait loci and the associated positional candidate genes affecting baseline serum Ca(++) concentrations. A genome-wide association study was conducted in an F(2) intercross population between Landrace and Korean native pigs using the porcine single nucleotide polymorphism (SNP) 60 K beadchip and the PLINK program based on linear regression. Data used in the study included 410 F(2) pigs. All experimental animals were genotyped with 36,613 SNP markers located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 7 and serum Ca(++) levels (DIAS0002191, genomic control-corrected P = 7.7 × 10(-5)). The position of DIAS0002191 was closely located to SLA class III region containing the C2 gene encoding the complementary component 2 protein, a protein which is important in the humoral immune responses. De novo sequencing of the porcine C2 gene revealed a missense mutation [c.1963A<G (N655D)] and this missense mutation was also strongly associated with serum Ca(++) concentrations (genomic control-corrected P = 5.9 × 10(-5)). Further studies are necessary to investigate the effect of this missense mutation at a functional-molecular level. In conclusion, the missense mutation of the C2 gene identified in this study may help in elucidating the genetic factors underlying humoral immune reactions.

Genome-wide association study of porcine hematological parameters in a Large White × Minzhu F2 resource population

Hematological traits, which are important indicators of immune function in animals, have been commonly examined as biomarkers of disease and disease severity in humans and animals. Genome-wide significant quantitative trait loci (QTLs) provide important information for use in breeding programs of animals such as pigs. QTLs for hematological parameters (hematological traits) have been detected in pig chromosomes, although these are often mapped by linkage analysis to large intervals making identification of the underlying mutation problematic. Single nucleotide polymorphisms (SNPs) are the common form of genetic variation among individuals and are thought to account for the majority of inherited traits. In this study, a genome-wide association study (GWAS) was performed to detect regions of association with hematological traits in a three-generation resource population produced by intercrossing Large White boars and Minzhu sows during the period from 2007 to 2011. Illumina PorcineSNP60 BeadChip technology was used to genotype each animal and seven hematological parameters were measured (hematocrit (HCT), hemoglobin (HGB), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV), red blood cell count (RBC) and red blood cell volume distribution width (RDW)). Data were analyzed in a three step Genome-wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) method. A total of 62 genome-wide significant and three chromosome-wide significant SNPs associated with hematological parameters were detected in this GWAS. Seven and five SNPs were associated with HCT and HGB, respectively. These SNPs were all located within the region of 34.6-36.5 Mb on SSC7. Four SNPs within the region of 43.7-47.0 Mb and fifty-five SNPs within the region of 42.2-73.8 Mb on SSC8 showed significant association with MCH and MCV, respectively. At chromosome-wide significant level, one SNP at 29.2 Mb on SSC1 and two SNPs within the region of 26.0-26.2 Mb were found to be significantly associated with RBC and RDW, respectively. Many of the SNPs were located within previously reported QTL regions and appeared to narrow down the regions compared with previously described QTL intervals. In current research, a total of seven significant SNPs were found within six candidate genes SCUBE3, KDR, TDO, IGFBP7, ADAMTS3 and AFP. In addition, the KIT gene, which has been previously reported to relate to hematological parameters, was located within the region significantly associated with MCH and MCV and could be a candidate gene. These results of this study may lead to a better understanding of the molecular mechanisms of hematological parameters in pigs.

QTL analysis of clinical-chemical traits in an F₂ intercross between Landrace and Korean native pigs

Clinical-chemical traits are essential when examining the health status of individuals. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting clinical-chemical traits in a reciprocal F(2) intercross between Landrace and Korean native pigs. Following an overnight fast, 25 serum phenotypes related to clinical-chemical traits (e.g., hepatic function parameters, renal function parameters, electrolyte, lipids) were measured in >970 F(2) progeny. All experimental samples were subjected to genotyping analysis using 165 microsatellite markers located across the genome. We identified eleven genome-wide significant QTL in six chromosomal regions (SSC 2, 7, 8, 13, 14, and 15) and 59 suggestive QTL in 17 chromosomal regions (SSC 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, and 18). We also observed significant effects of reciprocal crosses on some of the traits, which would seem to result from maternal effect, QTL on sex chromosomes, imprinted genes, or genetic difference in mitochondrial DNA. The role of genomic imprinting in clinical-chemical traits also was investigated. Genome-wide analysis revealed a significant evidence for an imprinted QTL in SSC4 affecting serum amylase levels. Additionally, a series of bivariate linkage analysis provided strong evidence that QTL in SSC 2, 13, 15, and 18 have a pleiotropic effect on clinical-chemical traits. In conclusion, our study detected both novel and previously reported QTL influencing clinical-chemical traits in pigs. The identified QTL together with the positional candidate genes identified here could play an important role in elucidating the genetic structure of clinical-chemical phenotype variation in humans and swine.