Sus Scrofa Chromosome Research Articles

Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances.

Selection affects the patterns of linkage disequilibrium (LD) around the site of a beneficial allele with an increase in LD among the hitchhiking alleles. Comparing the differences in regional LD between pig populations could help to identify putative genomic regions with potential adaptations for economic traits. In this study, using Illumina Porcine SNP60K BeadChip genotyping data from 207 Chinese indigenous, 117 South American village and 408 Large White pigs, we estimated the variation of genome-wide LD between populations using the varld program. The top 0.1% standardized VarLD scores were used as a criterion for all comparisons, and compared with LD blocks, a total of four selection signatures on Sus scrofa chromosome (SSC) 7, 9, 13 and 14 were identified in all populations. These signatures overlapped with quantitative trait loci for linoleic acid content, age at puberty, number of muscle fibers per unit area, hip structure and body weight traits in pigs. Among them, one of the signatures (56.5-56.6Mb on SSC7) in Large White pigs harbored the ADAMTSL3 gene, which is known to affect body length. The findings of this study seem to point toward recent selection in different pig populations. Further investigations are encouraged to confirm the selection signatures detected by varld in the present study.

Genome-wide association study identifies candidate genes for piglet splay leg syndrome in different populations

BackgroundPiglet splay leg syndrome (PSL) is one of the most frequent genetic defects, and can cause considerable economic loss in pig production. The present understanding of etiology and pathogenesis of PSL is poor. The current study focused on identifying loci associated with PSL through a genome-wide association study (GWAS) performed with the Illumina Porcine60 SNP Beadchip v2.0. The study was a case/control design with four pig populations (Duroc, Landrace, Yorkshire and one crossbred of Landrace × Yorkshire).ResultAfter quality control of the genotyping data, 185 animals (73 cases, 112 controls) and 43,495 SNPs were retained for further analysis. Principal components (PCs) identified from the genomic kinship matrix were included in the statistical model for correcting the effect of population structure. Seven chromosome-wide significant SNPs were identified on Sus scrofa chromosome 1 (SSC1), SSC2 (2 SNPs), SSC7, SSC15 (2 SNPs) and SSC16 after strict Bonferroni correction. Four genes (HOMER1 and JMY on SSC2, ITGA1 on SSC16, and RAB32 on SSC1) related to muscle development, glycogen metabolism and mitochondrial dynamics were identified as potential candidate genes for PSL.ConclusionsWe identified seven chromosome-wide significant SNPs associated with PSL and four potential candidate genes for PSL. To our knowledge, this is the first pilot study aiming to identify the loci associated with PSL using GWAS. Further investigations and validations for those findings are encouraged.

Haplotypes on pig chromosome 3 distinguish metabolically healthy from unhealthy obese individuals.

We have established a pig resource population specifically designed to elucidate the genetics involved in development of obesity and obesity related co-morbidities by crossing the obesity prone Göttingen Minipig breed with two lean production pig breeds. In this study we have performed genome wide association (GWA) to identify loci with effect on blood lipid levels. The most significantly associated single nucleotide polymorphisms (SNPs) were used for linkage disequilibrium (LD) and haplotype analyses. Three separate haploblocks which influence the ratio between high density lipoprotein cholesterol and total cholesterol (HDL-C/CT), triglycerides (TG) and low density lipoprotein cholesterol (LDL-C) levels respectively were identified on Sus Scrofa chromosome 3 (SSC3). Large additive genetic effects were found for the HDL-C/CT and LDL-C haplotypes. Haplotypes segregating from Göttingen Minipigs were shown to impose a positive effect on blood lipid levels. Thus, the genetic profile of the Göttingen Minipig breed seems to support a phenotype comparable to the metabolic healthy obese (MHO) phenotype in humans.

Genome-wide association studies reveal additional related loci for fatty acid composition in a Duroc pig multigenerational population.

The aim of the present study was to detect quantitative trait loci affecting fatty acid composition in back fat and intramuscular fat in a Duroc pig population comprising seventh-generation pedigrees using genome-wide association studies (GWAS). In total, 305 animals were genotyped using single nucleotide polymorphisms (SNPs) array and five selected SNPs from regions containing known candidate genes related to fatty acid synthesis or metabolism. In total, 24 genome-wide significant SNP regions were detected in 12 traits, and 76 genome-wide suggestive SNP regions were detected in 33 traits. The Sus scrofa chromosome (SSC) 7 at 10.3Mb was significantly associated with C17:0 in intramuscular fat, while the SSC9 at 13.6Mb was significantly associated with C14:0 in intramuscular fat. The SSC12 at 1.0Mb was significantly associated with C14:0 in back fat and the SSC14 at 121.0Mb was significantly associated with C18:0 in intramuscular fat. These regions not only replicated previously reported loci containing some candidate genes involved in fatty acid composition (fatty acid synthase and stearoyl-CoA desaturase) but also included several additional related loci.

Genome-wide detection of genetic markers associated with growth and fatness in four pig populations using four approaches

BackgroundGenome-wide association studies (GWAS) have been extensively used to identify genomic regions associated with a variety of phenotypic traits in pigs. Until now, most GWAS have explored single-trait association models. Here, we conducted both single- and multi-trait GWAS and a meta-analysis for nine fatness and growth traits on 2004 pigs from four diverse populations, including a White Duroc × Erhualian F2 intercross population and Chinese Sutai, Laiwu and Erhualian populations.ResultsWe identified 44 chromosomal regions that were associated with the nine traits, including four genome-wide significant single nucleotide polymorphisms (SNPs) on SSC2 (SSC for Sus scrofa chromosome), 4, 7 and X. Compared to the single-population GWAS, the meta-analysis was less powerful for the identification of SNPs with population-specific effects but more powerful for the detection of SNPs with population-shared effects. Multiple-trait analysis reduced the power to detect trait-specific SNPs but significantly enhanced the power to identify common SNPs across traits. The SNP on SSC7 had pleiotropic effects on the nine traits in the F2 and Erhualian populations. Another pleiotropic SNP was observed on SSCX for these traits in the F2 and Sutai populations. Both population-specific and shared SNPs were identified in this study, thus reflecting the complex genetic architecture of pig growth and fatness traits.ConclusionsWe demonstrate that the multi-trait method and the meta-analysis on multiple populations can be used to increase the power of GWAS. The two significant SNPs on SSC7 and X had pleiotropic effects in the F2, Erhualian and Sutai populations.

Genome-wide association study identifies 22 new loci for body dimension and body weight traits in a White Duroc×Erhualian F2 intercross population.

ObjectiveGrowth-related traits are important economic traits in the swine industry. However, the genetic mechanism of growth-related traits is little known. The aim of this study was to screen the candidate genes and molecular markers associated with body dimension and body weight traits in pigs.MethodsA genome-wide association study (GWAS) on body dimension and body weight traits was performed in a White Duroc×Erhualian F2 intercross by the illumina PorcineSNP60K Beadchip. A mixed linear model was used to assess the association between single nucleotide polymorphisms (SNPs) and the phenotypes.ResultsIn total, 611 and 79 SNPs were identified significantly associated with body dimension traits and body weight respectively. All SNPs but 62 were located into 23 genomic regions (quantitative trait loci, QTLs) on 14 autosomal and X chromosomes in Sus scrofa Build 10.2 assembly. Out of the 23 QTLs with the suggestive significance level (5×10−4), three QTLs exceeded the genome-wide significance threshold (1.15×10−6). Except the one on Sus scrofa chromosome (SSC) 7 which was reported previously all the QTLs are novel. In addition, we identified 5 promising candidate genes, including cell division cycle 7 for abdominal circumference, pleiomorphic adenoma gene 1 and neuropeptides B/W receptor 1 for both body weight and cannon bone circumference on SSC4, phosphoenolpyruvate carboxykinase 1, and bone morphogenetic protein 7 for hip circumference on SSC17.ConclusionThe results have not only demonstrated a number of potential genes/loci associated with the growth-related traits in pigs, but also laid a foundation for studying the genes’ role and further identifying causative variants underlying these loci.

Genetic analysis of teat number in pigs reveals some developmental pathways independent of vertebra number and several loci which only affect a specific side

BackgroundNumber of functional teats is an important trait in commercial swine production. As litter size increases, the number of teats must also increase to supply nutrition to all piglets. Therefore, a genome-wide association analysis was conducted to identify genomic regions that affect this trait in a commercial swine population. Genotypic data from the Illumina Porcine SNP60v1 BeadChip were available for 2951 animals with total teat number (TTN) records. A subset of these animals (n = 1828) had number of teats on each side recorded. From this information, the following traits were derived: number of teats on the left (LTN) and right side (RTN), maximum number of teats on a side (MAX), difference between LTN and RTN (L − R) and absolute value of L − R (DIF). Bayes C option of GENSEL (version 4.61) and 1-Mb windows were implemented. Identified regions that explained more than 1.5% of the genomic variation were tested in a larger group of animals (n = 5453) to estimate additive genetic effects.ResultsMarker heritabilities were highest for TTN (0.233), intermediate for individual side counts (0.088 to 0.115) and virtually nil for difference traits (0.002 for L − R and 0.006 for DIF). Each copy of the VRTN mutant allele increased teat count by 0.35 (TTN), 0.16 (LTN and RTN) and 0.19 (MAX). 15, 18, 13 and 18 one-Mb windows were detected that explained more than 1.0% of the genomic variation for TTN, LTN, RTN, and MAX, respectively. These regions cumulatively accounted for over 50% of the genomic variation of LTN, RTN and MAX, but only 30% of that of TTN. Sus scrofa chromosome SSC10:52 Mb was associated with all four count traits, while SSC10:60 and SSC14:54 Mb were associated with three count traits. Thirty-three SNPs accounted for nearly 39% of the additive genetic variation in the validation dataset. No effect of piglet sex or percentage of males in litter was detected, but birth weight was positively correlated with TTN.ConclusionsTeat number is a heritable trait and use of genetic markers would expedite selection progress. Exploiting genetic variation associated with teat counts on each side would enhance selection focused on total teat counts. These results confirm QTL on SSC4, seven and ten and identify a novel QTL on SSC14.

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 F2 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 F2 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 F2 intercross, our approach has limited power to distinguish one particular positional candidate gene from a QTL region.

Genetic and genomic basis of antibody response to porcine reproductive and respiratory syndrome (PRRS) in gilts and sows.

BackgroundOur recent research showed that antibody response to porcine reproductive and respiratory syndrome (PRRS), measured as sample-to-positive (S/P) ratio, is highly heritable and has a high genetic correlation with reproductive performance during a PRRS outbreak. Two major quantitative trait loci (QTL) on Sus scrofa chromosome 7 (SSC7; QTLMHC and QTL130) accounted for ~40 % of the genetic variance for S/P. Objectives of this study were to estimate genetic parameters for PRRS S/P in gilts during acclimation, identify regions associated with S/P, and evaluate the accuracy of genomic prediction of S/P across populations with different prevalences of PRRS and using different single nucleotide polymorphism (SNP) sets.MethodsPhenotypes and high-density SNP genotypes of female pigs from two datasets were used. The outbreak dataset included 607 animals from one multiplier herd, whereas the gilt acclimation (GA) dataset included data on 2364 replacement gilts from seven breeding companies placed on health-challenged farms. Genomic prediction was evaluated using GA for training and validation, and using GA for training and outbreak for validation. Predictions were based on SNPs across the genome (SNPAll), SNPs in one (SNPMHC and SNP130) or both (SNPSSC7) QTL, or SNPs outside the QTL (SNPRest).ResultsHeritability of S/P in the GA dataset increased with the proportion of PRRS-positive animals in the herd (from 0.28 to 0.47). Genomic prediction accuracies ranged from low to moderate. Average accuracies were highest when using only the 269 SNPs in both QTL regions (SNPSSC7, with accuracies of 0.39 and 0.31 for outbreak and GA validation datasets, respectively. Average accuracies for SNPALL, SNPMHC, SNP130, and SNPRest were, respectively, 0.26, 0.39, 0.21, and 0.05 for the outbreak, and 0.28, 0.25, 0.22, and 0.12, for the GA validation datasets.ConclusionsModerate genomic prediction accuracies can be obtained for PRRS antibody response using SNPs located within two major QTL on SSC7, while the rest of the genome showed limited predictive ability. Results were obtained using data from multiple genetic sources and farms, which further strengthens these findings. Further research is needed to validate the use of S/P ratio as an indicator trait for reproductive performance during PRRS outbreaks.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-016-0230-0) contains supplementary material, which is available to authorized users.

Genome-wide study refines the quantitative trait locus for number of ribs in a Large White×Minzhu intercross pig population and reveals a new candidate gene.

In China, sparerib is one of the most valuable parts of the pork carcass. As a result, candidate gene mining for number of ribs has become an interesting study focus. To examine the genetic basis for this major trait, we genotyped 596 individuals from an F2 Large White×Minzhu intercross pig population using the PorcineSNP60 Genotyping BeadChip. The genome-wide association study identified a locus for number of ribs in a 2.38-Mb region on Sus scrofa chromosome 7 (SSC7 of Sus scrofa genome assembly, Sscrofa10.2). We identified the top significant SNP ASGA0035536, which explained 16.51% of the phenotypic variance. A previously reported candidate causal mutation (g.19034 A>C) in vertebrae development-associated gene VRTN explained 8.79% of the phenotypic variation on number of ribs and had a much lower effect than ASGA0035536. Haplotype sharing analysis in F1 boars localized the rib number QTL to a 951-kb interval on SSC7. This interval encompassed 17 annotated genes in Sscrofa10.2, including the previously reported VRTN candidate gene. Of the 17 candidate genes, LTBP2, which encodes a latent transforming growth factor beta binding protein, was previously reported to indirectly regulate the activity of growth differentiation factor Gdf11, which has been shown to increase the number of ribs in knock-out mice. Thus, we propose LTBP2 as a good new candidate gene for number of ribs in the pig population. This finding advances our understanding of the genetic architecture of rib number in pigs.

Identification of loci affecting teat number by genome-wide association studies on three pig populations.

ObjectiveThree genome-wide association studies (GWAS) and a meta-analysis of GWAS were conducted to explore the genetic mechanisms underlying variation in pig teat number.MethodsWe performed three GWAS and a meta-analysis for teat number on three pig populations, including a White Duroc×Erhualian F2 resource population (n = 1,743), a Chinese Erhualian pig population (n = 320) and a Chinese Sutai pig population (n = 383).ResultsWe detected 24 single nucleotide polymorphisms (SNPs) that surpassed the genome-wide significant level on Sus Scrofa chromosomes (SSC) 1, 7, and 12 in the F2 resource population, corresponding to four loci for pig teat number. We highlighted vertnin (VRTN) and lysine demethylase 6B (KDM6B) as two interesting candidate genes at the loci on SSC7 and SSC12. No significant associated SNPs were identified in the meta-analysis of GWAS.ConclusionThe results verified the complex genetic architecture of pig teat number. The causative variants for teat number may be different in the three populations

Bioinformatic analyses in early host response to Porcine Reproductive and Respiratory Syndrome virus (PRRSV) reveals pathway differences between pigs with alternate genotypes for a major host response QTL

BackgroundA region on Sus scrofa chromosome 4 (SSC4) surrounding single nucleotide polymorphism (SNP) marker WUR10000125 (WUR) has been reported to be strongly associated with both weight gain and serum viremia in pigs after infection with PRRS virus (PRRSV). A proposed causal mutation in the guanylate binding protein 5 gene (GBP5) is predicted to truncate the encoded protein. To investigate transcriptional differences between WUR genotypes in early host response to PRRSV infection, an RNA-seq experiment was performed on globin depleted whole blood RNA collected on 0, 4, 7, 10 and 14 days post-infection (dpi) from eight littermate pairs with one AB (favorable) and one AA (unfavorable) WUR genotype animal per litter.ResultsGene Ontology (GO) enrichment analysis of transcripts that were differentially expressed (DE) between dpi across both genotypes revealed an inflammatory response for all dpi when compared to day 0. However, at the early time points of 4 and 7dpi, several GO terms had higher enrichment scores compared to later dpi, including inflammatory response (p < 10-7), specifically regulation of NFkappaB (p < 0.01), cytokine, and chemokine activity (p < 0.01). At 10 and 14dpi, GO term enrichment indicated a switch to DNA damage response, cell cycle checkpoints, and DNA replication. Few transcripts were DE between WUR genotypes on individual dpi or averaged over all dpi, and little enrichment of any GO term was found. However, there were differences in expression patterns over time between AA and AB animals, which was confirmed by genotype-specific expression patterns of several modules that were identified in weighted gene co-expression network analyses (WGCNA). Minor differences between AA and AB animals were observed in immune response and DNA damage response (p = 0.64 and p = 0.11, respectively), but a significant effect between genotypes pointed to a difference in ion transport/homeostasis and the participation of G-coupled protein receptors (p = 8e-4), which was reinforced by results from regulatory and phenotypic impact factor analyses between genotypes.ConclusionWe propose these pathway differences between WUR genotypes are the result of the inability of the truncated GBP5 of the AA genotyped pigs to inhibit viral entry and replication as quickly as the intact GBP5 protein of the AB genotyped pigs.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2547-z) contains supplementary material, which is available to authorized users.

Genome-wide association study reveals regions associated with gestation length in two pig populations.

Reproduction traits, such as gestation length (GLE), play an important role in dam line breeding in pigs. The objective of our study was to identify single nucleotide polymorphisms (SNPs) that are associated with GLE in two pig populations. Genotypes and deregressed breeding values were available for 2081 Dutch Landrace-based (DL) and 2301 Large White-based (LW) pigs. We identified two QTL regions for GLE, one in each population. For DL, three associated SNPs were detected in one QTL region spanning 0.52Mbp on Sus scrofa chromosome (SSC) 2. For LW, four associated SNPs were detected in one region of 0.14Mbp on SSC5. The region on SSC2 contains the heparin-binding EGF-like growth factor (HBEGF) gene, which promotes embryo implantation and has been described to be involved in embryo survival throughout gestation. The associated SNP can be used for marker-assisted selection in the studied populations, and further studies of the HBEGF gene are warranted to investigate its role in GLE.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population.

BackgroundIn many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.ResultsIn total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.ConclusionsTo the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

Multi-breed genome-wide association study reveals novel loci associated with the weight of internal organs.

BackgroundRecently, many genome-wide association studies (GWAS) have been conducted to understand the genetic architecture of economic important traits in farm animals. Pig is widely used as a biomedical animal model for its similarity with humans in terms of organ formation and disease mechanisms. Moreover, understanding the mechanisms that underlie the development of internal organs will impact the productive potential of pigs. Our aim was to uncover new single nucleotide polymorphisms (SNPs) associated with the weight of internal organs and carcass and also potential candidate genes.MethodsWe performed GWAS for the weight of heart, liver, spleen, kidney and carcass on five pig populations (White Duroc × Erhualian F2 intercross, Sutai population, Laiwu population, Erhualian population and commercial population, for a total of 2650 individuals). Genotype data was produced using the PorcineSNP60 Beadchip array. After quality control, the data was used for association tests under a general linear mixed model. Population stratification was adjusted by including a random polygenic effect based on a matrix of genotypic relationships. A meta-analysis of our GWAS datasets was conducted by summing up the Chi square values across breeds, with the degrees of freedom of the Chi square distribution equal to the effective number of breeds.ResultsThirty-nine quantitative trait loci (QTL) located on 15 chromosomes were identified by the single-population GWAS at the suggestive level. Among these, nine QTL surpassed the 5 % genome-wide significance threshold, including four for heart weight on SSC (Sus scrofa chromosome) 2, 4, 7 and 10, two for liver weight on SSC7, two for spleen weight on SSC5 and SSC7 and one for carcass weight on SSC11. The QTL on SSC7 showed pleiotropic effects for heart, liver and spleen weights in the F2 population. In addition, two QTL were detected in several populations, including one on SSC2 for heart weight in the F2 and Sutai populations and one on SSC7 for liver weight in the F2 and Laiwu populations. The meta-analysis detected four novel QTL on SSC1, 3, 8 and 16 for carcass weight.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-015-0168-7) contains supplementary material, which is available to authorized users.

Quantitative trait loci for the number of vertebrae on Sus scrofa chromosomes 1 and 7 independently influence the numbers of thoracic and lumbar vertebrae in pigs

Although quantitative trait loci (QTLs) for number of thoracic-lumbar vertebrae have been identified on Sus scrofa chromosomes (SSCs) 1 and 7, the influence of these QTLs on the thoracic and lumbar vertebrae is not clear. The aim of this study was to identify single nucleotide polymorphisms (SNPs) associated with total number of thoracic-lumbar vertebrae and for each trait (number of thoracic and lumbar vertebrae) separately. A total of 581 individuals from an F2 Large White×Minzhu population were genotyped using an SNP60K chip. Performing a genome-wide association study (GWAS) for total number of thoracic-lumbar vertebrae, 38 significant SNPs were identified in two QTL regions located on SSC1 and SSC7. Performing a GWAS for number of thoracic vertebrae only, 72 significant SNPs were located on SSC7. While performing a GWAS for number of lumbar vertebrae only, 17 significant SNPs were identified on SSC1. Gene mining suggested that the gene encoding orphan nuclear receptor, germ cell nuclear factor (NR6A1) on SSC1 was a strong candidate affecting the number of lumbar vertebrae in pigs. Additionally, genes encoding vertnin (VRTN), prospero homeobox 2 (PROX2), Finkel-Biskis-Jinkins murine osteosarcoma viral oncogene homolog (FOS), and transforming growth factor beta 3 (TGFB3) may be important candidates affecting the number of thoracic vertebrae in pigs. QTLs on SSC1 and SSC7 independently influenced the numbers of thoracic and lumbar vertebrae. These results shed light on the complex genetic background of vertebrae development in pigs.

A meta analysis of genome-wide association studies for limb bone lengths in four pig populations.

BackgroundLimb bone length is an economically important trait in pigs, because it is negatively correlated with backfat thickness, and is also a determinant to the yield of hip and loin. Moreover, abnormal growth of the limb bone leads to leg structural weakness. Until now, the genetic architecture of the pig lime bone length remains poorly understood. The object of this study was to map genomic loci for limb bone length by genome-wide association study (GWAS) on 4 pig populations.ResultsWe measured the lengths of five limb bones including scapula, humerus, ulna, femur and tibia that were dissected from the right-side carcass of 925, 331, 314 and 434 animals from White Duroc × Erhualian F2 intercross, Erhualian, Laiwu and Sutai populations, respectively. We genotyped the 2004 pigs for 62,163 single nucleotide polymorphisms (SNPs) on the Porcine SNP60 BeadChip, and performed GWAS and a GWAS meta analysis in the 4 populations. In total, we identified 12 and 4 loci associated with the limb bone lengths at suggestive and genome-wide significant levels respectively, of which 4 loci were reported for the first time. The most prominent locus was identified in a 924-kb (kilo base pairs) linkage disequilibrium block on Sus Scrofa chromosome (SSC) 7, and High Mobility Group AT-hook 1 (HMGA1) appears to be a strong candidate gene in this region. Another promising locus is located in the middle of SSC4, and Pleiomorphic Adenoma Gene 1 (PLAG1) is a functionally plausible candidate gene underlying the locus. Because the lengths of the 5 limb bones are highly correlated to each other, most of significant loci were associated with all of the 5 traits; however, several loci showed specific effect on the length of one limb bone, such as the locus at the proximal end of SSC2 associated with only the scapula length.ConclusionTo our knowledge, this study was the first GWAS meta analysis for limb bone lengths in pigs. As expected, the meta analysis is more powerful to identify genomic loci. A total of 16 loci were identified in this study, including four novel loci. HMGA1 and PLAG1 are two appearing candidate genes for pig limb bone lengths, which warrant further investigations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-015-0257-1) contains supplementary material, which is available to authorized users.

Identification of a putative quantitative trait nucleotide in guanylate binding protein 5 for host response to PRRS virus infection.

BackgroundPreviously, we identified a major quantitative trait locus (QTL) for host response to Porcine Respiratory and Reproductive Syndrome virus (PRRSV) infection in high linkage disequilibrium (LD) with SNP rs80800372 on Sus scrofa chromosome 4 (SSC4).ResultsWithin this QTL, guanylate binding protein 5 (GBP5) was differentially expressed (DE) (p < 0.05) in blood from AA versus AB rs80800372 genotyped pigs at 7,11, and 14 days post PRRSV infection. All variants within the GBP5 transcript in LD with rs80800372 exhibited allele specific expression (ASE) in AB individuals (p < 0.0001). A transcript re-assembly revealed three alternatively spliced transcripts for GBP5. An intronic SNP in GBP5, rs340943904, introduces a splice acceptor site that inserts five nucleotides into the transcript. Individuals homozygous for the unfavorable AA genotype predominantly produced this transcript, with a shifted reading frame and early stop codon that truncates the 88 C-terminal amino acids of the protein. RNA-seq analysis confirmed this SNP was associated with differential splicing by QTL genotype (p < 0.0001) and this was validated by quantitative capillary electrophoresis (p < 0.0001). The wild-type transcript was expressed at a higher level in AB versus AA individuals, whereas the five-nucleotide insertion transcript was the dominant form in AA individuals. Splicing and ASE results are consistent with the observed dominant nature of the favorable QTL allele. The rs340943904 SNP was also 100 % concordant with rs80800372 in a validation population that possessed an alternate form of the favorable B QTL haplotype.ConclusionsGBP5 is known to play a role in inflammasome assembly during immune response. However, the role of GBP5 host genetic variation in viral immunity is novel. These findings demonstrate that rs340943904 is a strong candidate causal mutation for the SSC4 QTL that controls variation in host response to PRRSV.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1635-9) contains supplementary material, which is available to authorized users.

Genome-wide association studies (GWAS) identify a QTL close to PRKAG3 affecting meat pH and colour in crossbred commercial pigs.

BackgroundImproving meat quality is a high priority for the pork industry to satisfy consumers’ preferences. GWAS have become a state-of-the-art approach to genetically improve economically important traits. However, GWAS focused on pork quality are still relatively rare.ResultsSix genomic regions were shown to affect loin pH and Minolta colour a* and b* on both loin and ham through GWAS in 1943 crossbred commercial pigs. Five of them, located on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, SSC16 and SSCX, were associated with meat colour. However, the most promising region was detected on SSC15 spanning 133–134 Mb which explained 3.51% - 17.06% of genetic variance for five measurements of pH and colour. Three SNPs (ASGA0070625, MARC0083357 and MARC0039273) in very strong LD were considered most likely to account for the effects in this region. ASGA0070625 is located in intron 2 of ZNF142, and the other two markers are close to PRKAG3, STK36, TTLL7 and CDK5R2. After fitting MARC0083357 (the closest SNP to PRKAG3) as a fixed factor, six SNPs still remained significant for at least one trait. Four of them are intragenic with ARPC2, TMBIM1, NRAMP1 and VIL1, while the remaining two are close to RUFY4 and CDK5R2. The gene network constructed demonstrated strong connections of these genes with two major hubs of PRKAG3 and UBC in the super-pathways of cell-to-cell signaling and interaction, cellular function and maintenance. All these pathways play important roles in maintaining the integral architecture and functionality of muscle cells facing the dramatic changes that occur after exsanguination, which is in agreement with the GWAS results found in this study.ConclusionsThere may be other markers and/or genes in this region besides PRKAG3 that have an important effect on pH and colour. The potential markers and their interactions with PRKAG3 require further investigationElectronic supplementary materialThe online version of this article (doi:10.1186/s12863-015-0192-1) contains supplementary material, which is available to authorized users.

Genome-wide association study using deregressed breeding values for cryptorchidism and scrotal/inguinal hernia in two pig lines.

BackgroundCryptorchidism and scrotal/inguinal hernia are the most frequent congenital defects in pigs. Identification of genomic regions that control these congenital defects is of great interest to breeding programs, both from an animal welfare point of view as well as for economic reasons. The aim of this genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) that are strongly associated with these congenital defects. Genotypes were available for 2570 Large White (LW) and 2272 Landrace (LR) pigs. Breeding values were estimated based on 1 359 765 purebred and crossbred male offspring, using a binary trait animal model. Estimated breeding values were deregressed (DEBV) and taken as the response variable in the GWAS.ResultsHeritability estimates were equal to 0.26 ± 0.02 for cryptorchidism and to 0.31 ± 0.01 for scrotal/inguinal hernia. Seven and 31 distinct QTL regions were associated with cryptorchidism in the LW and LR datasets, respectively. The top SNP per region explained between 0.96% and 1.10% and between 0.48% and 2.77% of the total variance of cryptorchidism incidence in the LW and LR populations, respectively. Five distinct QTL regions associated with scrotal/inguinal hernia were detected in both LW and LR datasets. The top SNP per region explained between 1.22% and 1.60% and between 1.15% and 1.46% of the total variance of scrotal/inguinal hernia incidence in the LW and LR populations, respectively. For each trait, we identified one overlapping region between the LW and LR datasets, i.e. a region on SSC8 (Sus scrofa chromosome) between 65 and 73 Mb for cryptorchidism and a region on SSC13 between 34 and 37 Mb for scrotal/inguinal hernia.ConclusionsThe use of DEBV in combination with a binary trait model was a powerful approach to detect regions associated with difficult traits such as cryptorchidism and scrotal/inguinal hernia that have a low incidence and for which affected animals are generally not available for genotyping. Several novel QTL regions were detected for cryptorchidism and scrotal/inguinal hernia, and for several previously known QTL regions, the confidence interval was narrowed down.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-015-0096-6) contains supplementary material, which is available to authorized users.