Investigation on muscle fiber types and meat quality and estimation of their heritability and correlation coefficients with each other in four pig populations.

The aim of this study was to investigate the muscle fiber types and meat quality in four populations and estimate the heritability and correlation coefficients of those traits in Shanxia long black pig (SX). In this study, a total of 318 pigs were recorded for 16 traits of the muscle fiber types and meat quality in four populations, including 256 individuals from the new breed SX. The population had a significant effect on all recorded traits, and the meat quality of the Lulai black pig was better than the remaining populations. The heritability (h2 ) of meat quality traits was from 0.06 (pH at 24 h) to 0.47 (shearing force), and the muscle fiber types belonged to the traits with low to medium heritability. The density of total fiber had the highest h2 (0.40), while the percentage of type IIA had the lowest h2 (0.04). Most traits are phenotypically correlated with each other, but only a small proportion of traits are genetically correlated with each other. None fiber type genetically correlated with meat quality significantly, because the genetic correlation coefficients had large standard errors. These results provided some insights into genetic improvements for the meat quality in pig breeds and also indicated that the parameters of muscle fiber characteristics can explain parts of the variation in meat quality.

Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations.

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

Association between MC1R gene and coat color segregation in Shanxia long black pig and Lulai black pig.

BACKGROUND: Coat color, as a distinct phenotypic characteristic of pigs, is often subject to preference and selection, such as in the breeding process of new breed. Shanxia long black pig was derived from an intercross between Berkshire boars and Licha black pig sows, and it was bred as a paternal strain with high-quality meat and black coat color. Although the coat color was black in the F1 generation of the intercross, it segregated in the subsequent generations. This study aims to decode the genetic basis of coat color segregation and develop a method to distinct black pigs from the spotted in Shanxia long black pig. RESULTS: Only a QTL was mapped at the proximal end of chromosome 6, and MC1R gene was picked out as functional candidate gene. A total of 11 polymorphic loci were identified in MC1R gene, and only the c.67_68insCC variant was co-segregating with coat color. This locus isn't recognized by any restriction endonuclease, so it can't be genotyped by PCR-RFLP. The c.370G > A polymorphic locus was also significantly associated with coat color, and has been in tightly linkage disequilibrium with the c.67_68insCC. Furthermore, it is recognized by BspHI. Therefore, a PCR-RFLP method was set up to genotype this locus. Besides the 175 sequenced individuals, another more 1,391 pigs were genotyped with PCR-RFLP, and all of pigs with GG (one band) were black. CONCLUSION: MC1R gene (c.67_68insCC) is the causative gene (mutation) for the coat color segregation, and the PCR-RFLP of c.370G > A could be used in the breeding program of Shanxia long black pig.

Patient-directed vs. fixed-volume PEG for colonoscopy preparation: a randomized controlled trial.

Background: Individualization using different volumes of polyethylene glycol is widely regarded as the optimal solution for bowel preparation, while the patient-directed regimen we propose may serve as a reliable individual solution. This study aimed to assess the efficacy, safety, and satisfaction of bowel preparation with a patient-directed regimen. Methods: Patients in the fixed-volume group ingested the same amount of PEG, while those in patient-directed group ingested different amount according to stool consistency or stool water content. Results: After filtering by exclusion criteria, 428 individuals in the fixed-volume group and 103 in the patient-directed group were successfully enrolled and analyzed. Eighty-three (80.6%) individuals in the patient-directed group had a reduced polyethylene glycol volume. There was no significant difference in the bowel preparation efficacy between the two groups (90.0% vs. 90.3%, χ² = 0.01; p = 0.918). Patients in the patient-directed group complained of fewer adverse effects (53.0% vs. 36.9%, χ² = 8.655; p = 0.003), especially vomiting (13.6% vs. 1.0%, χ² = 13.304; p < 0.001). Regarding comfort during bowel preparation, the degree of comfort was not significantly different between groups. Furthermore, the willingness rate for further colonoscopy in the patient-directed group was significantly higher than that in the fixed-volume group (90.3% vs. 77.1%, χ² = 8.912; p < 0.05). Multivariable logistic regression analysis showed that the body mass index served as an independent factor impacting quality of bowel preparation with the patient-directed regimen (OR 1.16, 95% CI 1.00-1.34; p = 0.043). Conclusions: Without decreasing the bowel preparation efficacy, the patient-directed regimen increased the safety and satisfaction of bowel preparation and is expected to be a regular and individual solution for bowel preparation. Individuals with a lower body mass index are more likely to undertake this new regimen. Trial registration number: ChiCTR1900022072 at ChiClinicalTrials.gov.

Fusobacterium nucleatum Promotes the Development of Ulcerative Colitis by Inducing the Autophagic Cell Death of Intestinal Epithelial.

There is a growing body of evidence which suggests that intestinal microbiota, especially Fusobacterium nucleatum (F. nucleatum), are associated with intestinal immune disease such as ulcerative colitis (UC). The mechanism by which F. nucleatum promotes intestinal epithelial cell (IEC) death remained undefined. Here, we investigated the potential mechanisms about how F. nucleatum aggravates IEC death in UC. We first detected the abundance of F. nucleatum in UC tissues and analyzed its relationship with the clinical characteristics of UC. Next, we explored whether F. nucleatum promotes intestinal epithelial cell death in vitro and in vivo. Furthermore, we extracted lipopolysaccharide (LPS) of the F. nucleatum and examined whether F. nucleatum exacerbates UC via LPS. Our results indicated that F. nucleatum was abundant in UC tissues and was correlated with clinical characteristics. In addition, we demonstrated that F. nucleatum and its LPS aggravated IEC death by promoted IEC autophagy. Furthermore, autophagy inhibitors, chloroquine (CQ), 3-methyladenine (3-MA) or Atg5 silencing prevented IEC death mediated by F. nucleatum, which suggests F. nucleatum may contribute to UC by activating autophagic cell death. All our results uncover a vital role of F. nucleatum in autophagic cell death and UC, giving rise to a new sight for UC therapy by inhibiting excessive IEC autophagy and autophagic cell death.

New insights into host adaptation to swine respiratory disease revealed by genetic differentiation and RNA sequencing analyses.

Swine respiratory disease (SRD) causes massive economic losses in the swine industry and is difficult to control and eradicate on pig farms. Here, we employed population genetics and transcriptomics approaches to decipher the molecular mechanism of host adaptation to swine respiratory disease. We recorded two SRD-related traits, the enzootic pneumonia-like (EPL) score and lung lesion (LL) levels, and performed four body weight measurements, at ages of 150, 180, 240, and 300 days, in a Chinese Bamaxiang pig herd (n = 314) raised under consistent indoor rearing conditions. We divided these animals into disease-resistant and disease-susceptible groups based on the most likely effects of both SRD-related traits on their weight gain, and performed genetic differentiation analyses in these two groups. Significant loci showing the top 1% of genetic differentiation values, exceeding the threshold of p = 0.005 set based on 1,000-times permutation tests, were defined as candidate regions related to host resistance or susceptibility to SRD. We identified 107 candidate genes within these regions, which are mainly involved in the biological processes of immune response, fatty acid metabolism, lipid metabolism, and growth factor signaling pathways. Among these candidate genes, TRAF6, CD44, CD22, TGFB1, CYP2B6, and SNRPA were highlighted due to their central regulatory roles in host immune response or fat metabolism and their differential expression between healthy lung tissues and lung lesions. These findings advance our understanding of the molecular mechanisms of host resistance or susceptibility to respiratory disease in pigs and are of significance for the breeding pigs resistant to respiratory disease in the swine industry.

Identification of copy number variations using high density whole-genome SNP markers in Chinese Dongxiang spotted pigs.

OBJECTIVE: Copy number variations (CNVs) are a major source of genetic diversity complementary to single nucleotide polymorphism (SNP) in animals. The aim of the study was to performed a comprehensive genomic analysis of CNVs based on high density whole-genome SNP markers in Chinese Dongxiang spotted pigs. METHODS: We used customized Affymetrix Axiom Pig1.4M array plates containing 1.4 million SNPs and the PennCNV algorithm to identify porcine CNVs on autosomes in Chinese Dongxiang spotted pigs. Then, the next generation sequence data was used to confirm the detected CNVs. Next, functional analysis was performed for gene contents in copy number variation regions (CNVRs). In addition, we compared the identified CNVRs with those reported ones and quantitative trait loci (QTL) in the pig QTL database. RESULTS: We identified 871 putative CNVs belonging to 2,221 CNVRs on 17 autosomes. We further discarded CNVRs that were detected only in one individual, leaving us 166 CNVRs in total. The 166 CNVRs ranged from 2.89 kb to 617.53 kb with a mean value of 93.65 kb and a genome coverage of 15.55 Mb, corresponding to 0.58% of the pig genome. A total of 119 (71.69%) of the identified CNVRs were confirmed by next generation sequence data. Moreover, functional annotation showed that these CNVRs are involved in a variety of molecular functions. More than half (56.63%) of the CNVRs (n = 94) have been reported in previous studies, while 72 CNVRs are reported for the first time. In addition, 162 (97.59%) CNVRs were found to overlap with 2,765 previously reported QTLs affecting 378 phenotypic traits. CONCLUSION: The findings improve the catalog of pig CNVs and provide insights and novel molecular markers for further genetic analyses of Chinese indigenous pigs.

Microbial communities in swine lungs and their association with lung lesions.

Under natural farming, environmental pathogenic microorganisms may invade and affect swine lungs, further resulting in lung lesions. However, few studies on swine lung microbiota and their potential relationship with lung lesions were reported. Here, we sampled 20 pigs from a hybrid herd raised under natural conditions; we recorded a lung-lesion phenotype and investigated lung microbial communities by sequencing the V3-V4 region of 16S rRNA gene for each individual. We found reduced microbial diversity but more biomass in the severe-lesion lungs. Methylotenera, Prevotella, Sphingobium and Lactobacillus were the prominent bacteria in the healthy lungs, while Mycoplasma, Ureaplasma, Sphingobium, Haemophilus and Phyllobacterium were the most abundant microbes in the severe-lesion lungs. Notably, we identified 64 lung-lesion-associated OTUs, of which two classified to Mycoplasma were positively associated with lung lesions and 62 showed negative association including thirteen classified to Prevotella and six to Ruminococcus. Cross-validation analysis showed that lung microbiota explained 23.7% phenotypic variance of lung lesions, suggesting that lung microbiota had large effects on promoting lung healthy. Furthermore, 22 KEGG pathways correlated with lung lesions were predicted. Altogether, our findings improve the knowledge about swine lung microbial communities and give insights into the relationship between lung microbiota and lung lesions.

Copy number variation in the MSRB3 gene enlarges porcine ear size through a mechanism involving miR-584-5p.

BACKGROUND: The size and type of ears are important conformation characteristics that distinguish pig breeds. A significant quantitative trait locus (QTL) for ear size has been identified on SSC5 (SSC for Sus scrofa chromosome) but the underlying causative gene and mutation remain unknown. Thus, our aim was to identify the gene responsible for enlarged ears in pig. RESULTS: First, we narrowed down the QTL region on SSC5 to a 137.85-kb interval that harbors only the methionine sulfoxide reductase B3 (MSRB3) gene. Then, we identified a 38.7-kb copy number variation (CNV) that affects the last two exons of MSRB3 and could be the candidate causative mutation for this QTL. This CNV showed complete concordance with genotype at the QTL of the founder animals in a white Duroc × Erhualian F2 intercross and was found only in pigs from six Chinese indigenous breeds with large ears and from the Landrace breed with half-floppy ears. Moreover, it accounted for the significant association with ear size on SSC5 across the five pig populations tested. eQTL mapping revealed that this CNV was significantly associated with the expression of the microRNA (miRNA) miR-584-5p, which interacts with MSRB3, one of its target genes. In vivo and in vitro experiments confirmed that miR-584-5p inhibits the translation of MSRB3 mRNA. Taken together, these results led us to conclude that presence of the 38.7-kb CNV in the genome of some pig breeds affects ear size by altering the expression of miR-584-5p, which consequently hinders the expression of one of its target genes (e.g. MSRB3). CONCLUSIONS: Our findings shed insight into the underlying mechanism of development of external ears in mammals and contribute to a better understanding of how the presence of CNV can regulate gene expression.

A genome-wide association study identifies genomic loci associated with backfat thickness, carcass weight, and body weight in two commercial pig populations.

Growth and fatness traits are economically important in the pig industry. To dissect the genetic architecture of these traits in commercial pigs, we conducted a genome-wide association study (GWAS) for carcass weight, backfat thickness, and body weight in two commercial populations: Duroc × (Landrace × Yorkshire) (DLY) and Duroc populations. To enhance the detection power, three GWAS approaches including single-trait GWAS, multi-trait GWAS and meta-analysis were used in this study. A total of 13 suggestive loci were identified on nine chromosomes. The most significant locus was detected at 272.05 Mb on SSC1, and it was associated with backfat thickness at the first rib in the DLY population. Three genes at the identified loci (TBC1D1, BAAT and PHLPP1) were highlighted as functionally plausible candidate genes for pig growth and fatness traits. Genome-wide significant locus was not evidenced in this study, indicating that large populations are required to identify QTL with minor effects on growth and fatness traits in commercial pig populations, in which intensively artificial selections have been imposed on these traits and small genetical variances usually retain in these traits.

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

BACKGROUND: Genome-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. RESULTS: We 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. CONCLUSIONS: We 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.

OBJECTIVE: Growth-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. METHODS: A 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. RESULTS: In 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. CONCLUSION: The 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.

Consequences of Epistasis on Growth in an Erhualian × White Duroc Pig Cross.

Epistasis describes an interaction between the effects of loci. We included epistasis in quantitative trait locus (QTL) mapping of growth at a series of ages in a cross of a Chinese pig breed, Erhualian, with a commercial line, White Duroc. Erhualian pigs have much lower growth rates than White Duroc. We improved a method for genomewide testing of epistasis and present a clear analysis workflow. We also suggest a new approach for interpreting epistasis results where significant additive and dominance effects of a locus in specific backgrounds are determined. In total, seventeen QTL were found and eleven showed epistasis. Loci on chromosomes 2, 3, 4 and 7 were highlighted as affecting growth at more than one age or forming an interaction network. Epistasis resulted in both the QTL on chromosomes 3 and 7 having effects in opposite directions. We believe it is the first time for the chromosome 7 locus that an allele from a Chinese breed has been found to decrease growth. The consequences of epistasis were diverse. Results were impacted by using growth rather than body weight as the phenotype and by correcting for an effect of mother. Epistasis made a considerable contribution to growth in this population and modelling epistasis was important for accurately determining QTL effects.

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

OBJECTIVE: Three genome-wide association studies (GWAS) and a meta-analysis of GWAS were conducted to explore the genetic mechanisms underlying variation in pig teat number. METHODS: We 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). RESULTS: We 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. CONCLUSION: The results verified the complex genetic architecture of pig teat number. The causative variants for teat number may be different in the three populations.

Multi-breed genome-wide association study reveals heterogeneous loci associated with loin eye area in pigs.

Numerous quantitative trait loci (QTL) for loin eye area had been identified by linkage mapping studies, but the lack of their precise position hinders their application in the pig breeding industry. To map QTL for loin eye area to a precise genomic region, we conducted a genome-wide association study (GWAS) using Illumina 60 K PorcineSNP60 Beadchip in four swine populations: 819 F2 pigs, 273 Laiwu pigs, 434 Sutai pigs, and 326 Erhualian pigs. In total, 26 single nucleotide polymorphisms (SNPs) deposited on seven chromosomes associated with loin eye area were identified, 11 of which surpassed the genome-wide significant threshold; of the 11 SNPs, seven located on SSC2 in F2 pigs and four located on SSC12 and SSC18 in Laiwu pigs. Of note, all of the identified QTL were breed specific and no common QTL was identified across the four populations in our study. These findings not only confirmed a previous QTL on SSC2 harboring the candidate gene insulin-like growth factor 2 (IGF2), but also identified some novel candidate genes, far upstream element binding protein 3 (FUBP3), myosin heavy chain (MYH) family, leucine-rich repeats and guanylate kinase domain containing (LRGUK). Our study will contribute to the further identification of the causal mutation underlying these QTL and improve our knowledge of the complex genetic architecture for loin eye area in pigs.

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

BACKGROUND: Limb 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. RESULTS: We 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. CONCLUSION: To 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.

A systems genetics study of swine illustrates mechanisms underlying human phenotypic traits.

BACKGROUND: The pig, which shares greater similarities with human than with mouse, is important for agriculture and for studying human diseases. However, similarities in the genetic architecture and molecular regulations underlying phenotypic variations in humans and swine have not been systematically assessed. RESULTS: We systematically surveyed ~500 F2 pigs genetically and phenotypically. By comparing candidates for anemia traits identified in swine genome-wide SNP association and human genome-wide association studies (GWAS), we showed that both sets of candidates are related to the biological process "cellular lipid metabolism" in liver. Human height is a complex heritable trait; by integrating genome-wide SNP data and human adipose Bayesian causal network, which closely represents bone transcriptional regulations, we identified PLAG1 as a causal gene for limb bone length. This finding is consistent with GWAS findings for human height and supports the common genetic architecture between swine and humans. By leveraging a human protein-protein interaction network, we identified two putative candidate causal genes TGFB3 and DAB2IP and the known regulators MESP1 and MESP2 as responsible for the variation in rib number and identified the potential underlying molecular mechanisms. In mice, knockout of Tgfb3 and Tgfb2 together decreases rib number. CONCLUSION: Our findings show that integrative network analyses reveal causal regulators underlying the genetic association of complex traits in swine and that these causal regulators have similar effects in humans. Thus, swine are a potentially good animal model for studying some complex human traits that are not under intense selection.

Sexually dimorphic genetic architecture of complex traits in a large-scale F2 cross in pigs.

BACKGROUND: It is common for humans and model organisms to exhibit sexual dimorphism in a variety of complex traits. However, this phenomenon has rarely been explored in pigs. RESULTS: To investigate the genetic contribution to sexual dimorphism in complex traits in pigs, we conducted a sex-stratified analysis on 213 traits measured in 921 individuals produced by a White Duroc × Erhualian F2 cross. Of the 213 traits examined, 102 differed significantly between the two sexes (q value <0.05), which indicates that sex is an important factor that influences a broad range of traits in pigs. We compared the estimated heritability of these 213 traits between males and females. In particular, we found that traits related to meat quality and fatty acid composition were significantly different between the two sexes, which shows that genetic factors contribute to variation in sexual dimorphic traits. Next, we performed a genome-wide association study (GWAS) in males and females separately; this approach allowed us to identify 13.6% more significant trait-SNP (single nucleotide polymorphism) associations compared to the number of associations identified in a GWAS that included both males and females. By comparing the allelic effects of SNPs in the two sexes, we identified 43 significant sexually dimorphic SNPs that were associated with 22 traits; 41 of these 43 loci were autosomal. The most significant sexually dimorphic loci were found to be associated with muscle hue angle and Minolta a* values (which are parameters that reflect the redness of meat) and were located between 9.3 and 10.7 Mb on chromosome 6. A nearby gene i.e. NUDT7 that plays an important role in heme synthesis is a strong candidate gene. CONCLUSIONS: This study illustrates that sex is an important factor that influences phenotypic values and modifies the effects of the genetic variants that underlie complex traits in pigs; it also emphasizes the importance of stratifying by sex when performing GWAS.

A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle.

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F2 intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F2 intercross and an F19 Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F2 animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F2 and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.

Three novel quantitative trait loci for skin thickness in swine identified by linkage and genome-wide association studies.

Skin is the largest organ in the pig body and plays a key role in protecting the body against pathogens and excessive water loss. Deciphering the genetic basis of swine skin thickness would enrich our knowledge about the skin. To identify the loci for porcine skin thickness, we first performed a genome scan with 194 microsatellite markers in a White Duroc × Erhualian F2 intercross. We identified three genome-wide significant QTL on pig chromosomes (SSC) 4, 7 and 15 using linkage analysis. The most significant QTL was found on SSC7 with a small confidence interval of ~5 cM, explaining 23.9 percent of phenotypic variance. Further, we conducted a genome-wide association study (GWAS) using Illumina PorcineSNP60 Beadchips for the F2 pedigree and a population of Chinese Sutai pigs. We confirmed significant QTL in the F2 pedigree and replicated QTL on SSC15 in Chinese Sutai pigs. A meta-analysis of GWASs on both populations detected a genomic region associated with skin thickness on SSC4. GWAS results were generally consistent with QTL mapping. Identical-by-descent analysis defined QTL on SSC7 in a 683-kb region harboring an interesting candidate gene: HMGA1. On SSC15, the linkage disequilibrium analysis showed a haplotype block of 2.20 Mb that likely harbors the gene responsible for skin thickness. Our findings provide novel insights into the genetic basis of swine skin thickness, which would benefit further understanding of porcine skin function.