Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig.

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Associations of genome-wide structural variations with phenotypic differences in cross-bred Eurasian pigs.

BACKGROUND: During approximately 10,000 years of domestication and selection, a large number of structural variations (SVs) have emerged in the genome of pig breeds, profoundly influencing their phenotypes and the ability to adapt to the local environment. SVs (≥ 50 bp) are widely distributed in the genome, mainly in the form of insertion (INS), mobile element insertion (MEI), deletion (DEL), duplication (DUP), inversion (INV), and translocation (TRA). While studies have investigated the SVs in pig genomes, genome-wide association studies (GWAS)-based on SVs have been rarely conducted. RESULTS: Here, we obtained a high-quality SV map containing 123,151 SVs from 15 Large White and 15 Min pigs through integrating the power of several SV tools, with 53.95% of the SVs being reported for the first time. These high-quality SVs were used to recover the population genetic structure, confirming the accuracy of genotyping. Potential functional SV loci were then identified based on positional effects and breed stratification. Finally, GWAS were performed for 36 traits by genotyping the screened potential causal loci in the F2 population according to their corresponding genomic positions. We identified a large number of loci involved in 8 carcass traits and 6 skeletal traits on chromosome 7, with FKBP5 containing the most significant SV locus for almost all traits. In addition, we found several significant loci in intramuscular fat, abdominal circumference, heart weight, and liver weight, etc. CONCLUSIONS: We constructed a high-quality SV map using high-coverage sequencing data and then analyzed them by performing GWAS for 25 carcass traits, 7 skeletal traits, and 4 meat quality traits to determine that SVs may affect body size between European and Chinese pig breeds.

Characteristics of Transcriptome and Metabolome Concerning Intramuscular Fat Content in Beijing Black Pigs.

To study the characteristics of genes and metabolites related to intramuscular fat (IMF) content with less influence by breed background and individual differences, the skeletal muscle samples from 40 Beijing black pigs with either high or low IMF content were used to perform transcriptome and metabolome analyses. About 99 genes (twofold-change) were differentially expressed. Up-regulated genes in the high IMF pigs were mainly related to fat metabolism. The key genes in charge of IMF deposition are ADIPOQ, CIDEC, CYP4B1, DGAT2, LEP, OPRL1, PLIN1, SCD, and THRSP. KLHL40, TRAFD1, and HSPA6 were novel candidate genes for the IMF trait due to their high abundances. In the low IMF pigs, the differentially expressed genes involved in virus resistance were up-regulated. About 16 and 18 differential metabolites (1.5 fold-change) were obtained in the positive and negative modes, respectively. Pigs with low IMF had weaker fatty acid oxidation due to the down-regulation of various carnitines. Differentially expressed genes were more important in determining IMF deposition than differential metabolites because relatively few differential metabolites were obtained, and they were merely the products under the physiological status of diverged IMF content. This study provided valuable information for further studies on IMF deposition.

Metabolite and Proteomic Profiling of Serum Reveals the Differences in Molecular Immunity between Min and Large White Pig Breeds.

Pig diseases seriously threaten the health of pigs and the benefits of pig production. Previous research has indicated that Chinese native pigs, such as the Min (M) pig, has a better disease resistance ability than Large White (LW) pigs. However, the molecular mechanism of this resistance is still unclear. In our study, we used serum untargeted metabolomics and proteomics, interrogated to characterize differences in the molecular immunities between six resistant and six susceptible pigs raised in the same environment. A total of 62 metabolites were identified as being significantly exhibited in M and LW pigs. Ensemble feature selection (EFS) machine learning methods were used to predict biomarkers of metabolites and proteins, and the top 30 were selected and retained. Weighted gene co-expression network analysis (WGCNA) confirmed that four key metabolites, PC (18:1 (11 Z)/20:0), PC (14:0/P-18: 0), PC (18:3 (6 Z, 9 Z, 12 Z)/16:0), and PC (16:1 (9 Z)/22:2 (13 Z, 16 Z)), were significantly associated with phenotypes, such as cytokines, and different pig breeds. Correlation network analysis showed that 15 proteins were significantly correlated with the expression of both cytokines and unsaturated fatty acid metabolites. Quantitative trait locus (QTL) co-location analysis results showed that 13 of 15 proteins co-localized with immune or polyunsaturated fatty acid (PUFA)-related QTL. Moreover, seven of them co-localized with both immune and PUFA QTLs, including proteasome 20S subunit beta 8 (PSMB8), mannose binding lectin 1 (MBL1), and interleukin-1 receptor accessory protein (IL1RAP). These proteins may play important roles in regulating the production or metabolism of unsaturated fatty acids and immune factors. Most of the proteins could be validated with parallel reaction monitoring, which suggests that these proteins may play an essential role in producing or regulating unsaturated fatty acids and immune factors to cope with the adaptive immunity of different pig breeds. Our study provides a basis for further clarifying the disease resistance mechanism of pigs.

Metabolomics and lipidomics profiles related to intramuscular fat content and flavor precursors between Laiwu and Yorkshire pigs.

Chinese indigenous pig breeds have higher intramuscular fat content (IMF) and better meat quality than Western commercial pigs. The differential metabolites and lipids in the skeletal muscle associated with IMF contents and meat flavor in Laiwu and Yorkshire pigs were investigated in this study. As a result, 113 differential metabolites and 54 differential lipids were discovered. Lipidomics revealed that the Laiwu pig had a fast lipid droplet formation and contained more triglyceride than the Yorkshire pig, which was corresponded to its high IMF contents. Both the lipidomics and metabolomics results indicated that the Laiwu pig had a higher mitochondrial content and aerobic respiration, due to its larger percentage of oxidative fibers. In addition, differential metabolites, such as oxoglutaric acid, fumarate, and l-aspartate, were thought to be important flavor precursors contributing to the Laiwu pig's improved pork taste.

Integrating genome-wide association study with RNA-seq revealed DBI as a good candidate gene for intramuscular fat content in Beijing black pigs.

Increasing intramuscular fat (IMF) content can enhance the sensory quality of meat, including tenderness, juiciness, flavor, and color. Genome-wide association study and RNA-sequencing (RNA-seq) analysis were used to identify candidate IMF genes in Beijing Black pigs, a popular species among consumers in northern China. Two and three single nucleotide polymorphisms were significantly associated with IMF in SSC13 and SSC15 respectively. Solute carrier family 4 member 7 (SLC4A7) on SSC13 and insulin induced gene 2 (INSIG2), coiled-coil domain containing 93 (CCDC93), and diazepam binding inhibitor acyl-CoA binding protein (DBI) on SSC15 are good candidate genes in this population. Furthermore, RNA-seq analysis was performed between high and low IMF groups, and 534 differentially expressed genes were identified. In addition, based on differentially expressed genes, Kyoto Encyclopedia of Genes and Genomes analysis revealed that peroxisome proliferator-activated receptors and FoxO signaling pathway pathways might contribute to IMF. Moreover, the DBI gene was identified as a candidate for IMF both by genome-wide association study and RNA-seq analysis, suggesting that it might be a crucial candidate gene for influencing IMF in Beijing Black pigs.

Characterization and difference of lipids and metabolites from Jianhe White Xiang and Large White pork by high-performance liquid chromatography-tandem mass spectrometry.

The quality of pork, such as intramuscular fat (IMF) content and flavor, can affect the acceptance of the consumer. Many studies have reported on the pork quality of Chinese local and commercial pigs (for example Large White (LW) pigs). The Jianhe White Xiang (JWX) pig, one of the Chinese Xiang pigs, is known for its high IMF and pork flavor. However, studies investigating the characterization and difference of lipids and metabolites between the JWX and LW pork are limited. Herein, we performed metabolomic and lipidomic profiling of JWX and LW pork by high-performance liquid chromatography-tandem mass spectrometry. The IMF and meat redness (a*) of the JWX pork were significantly higher than those of the LW pork. Metabolomic profiling revealed that 118 out of 501 polar metabolites, such as carnitine, amino acids, sugar, and dipeptide, were significantly different between the two types of pork. Additionally, the screened metabolites were mainly related to carnitine synthesis, phospholipid metabolism, sugar metabolism, and amino acid metabolism. Lipidomic profiling identified 100 of 376 lipids, which contained carnitine, diglyceride, triglyceride (TG), sphingomyelin, cardiolipin, fatty acid, phosphatidylcholine (PC), and phosphatidylethanolamine (PE), which were significantly different between the two types of pork under the positive and negative ion modes (variable importance in projection (VIP) > 1, p < 0.05, and fold change (FC) > 2 or FC < 0.5). All of the different TG substances were up-regulated in the JWX pork, and their carbon chain length was longer than that of the residual TGs. In addition, the JWX pork had more double bonds of PC and PE substances than LW pork. Thus, our findings provide comprehensive metabolomic and lipidomic profiles between the JWX and LW pork and a basic understanding on increasing the pork quality.

Integrated Proteotranscriptomics Reveals Differences in Molecular Immunity between Min and Large White Pig Breeds.

Long-term selection or evolution is an important factor governing the development of disease resistance in pigs. To better clarify the molecular mechanisms underlying different levels of disease resistance, we used transcriptomics and proteomics analysis to characterize differences in the immunities between six resistant (Min pig) and six susceptible (Large White, LW) pigs which were raised in the same environment. A total of 135 proteins and 791 genes were identified as being differentially expressed between the Large White and Min pig groups. Protein expression clustering and functional analysis revealed that proteins related to immune system process, humoral immune response, the B cell receptor signaling pathway, lymphocyte-mediated immunity, and innate immune responses were more highly expressed in Min pigs. Transcriptome gene set enrichment analysis was used to reveal that pathways of cell adhesion molecules and antigen processing and presentation are significantly enriched in Min pigs. Integrated proteomics and transcriptomics data analysis identified 16 genes that are differentially expressed at both the mRNA and protein levels. In addition, 13 out of these 16 genes were related to the quantitative trait loci of immune diseases, including neural EGFL-like 2 (NELL2) and lactate dehydrogenase B (LDHB), which are involved in innate immunity. Correlation analysis between the genes/proteins and cytokines shows upregulated proteins in LW pigs in association with immunosuppressive/pro-inflammatory cytokines, such as interleukin (IL) 10, IL6, and tumor necrosis factor alpha. This was further validated using parallel reaction monitoring analysis. In summary, we discovered several potential candidate pathways and key genes/proteins involved in determining differences in disease resistance between the two studied pig breeds, which could provide new insights into the breeding of pigs for disease resistance.

Integrating Genome-Wide Association Study with RNA-Sequencing Reveals HDAC9 as a Candidate GeneInfluencing Loin Muscle Area in Beijing Black Pigs.

Loin muscle area (LMA) is an important meat production trait and plays a key role in determining carcass leanness. Genome-wide association study (GWAS) and RNA sequencing (RNA-seq) analysis were used to identify candidate LMA genes in Beijing Black pigs, a popular breed among consumers in northern China. Ten single nucleotide polymorphisms (SNPs) in sus scrofa chromosome (SSC) 9 were significantly associated with LMA. These SNPs were mapped to a 2.90 Mb (84.94-87.84 Mb) region. A total of 11 annotated genes were mapped on this region, namely MEOX2, CRPPA, SOSTDC1, LRRC72, ANKMY2, BZW2, TSPAN13, AGR2, AHR, SNX13, and HDAC9. In addition, RNA-seq analysis was performed between the high- and low-LMA groups, and 329 differentially expressed genes (DEGs) were identified. Further, Kyoto Encyclopedia of Genes and Genomes analysis based on DEGs revealed that the JAK/STAT signaling pathway and oxytocin signaling pathway may be responsible for LMA. Both GWAS and RNA-seq analysis identified the HDAC9 gene, indicating that it may be an important candidate gene affecting LMA in Beijing Black pigs. The findings provide valuable molecular insights into the mechanisms that influence LMA content in pigs, which can be utilized in targeted approaches to enhance meat quality and commercial profitability.

A strategy for accurately and sensitively quantifying free and esterified fatty acids using liquid chromatography mass spectrometry.

Fatty acid (FA) composition of foods dictates a diversity of aspects regarding food quality, ranging from product shelf life, sensory properties to nutrition. There is a challenge to quantitate FAs using liquid chromatography-mass spectrometry due to poor ionization efficiency and matrix effects. Here, an isotopic-tagged derivatization strategy was established to accurately and sensitively quantify free and esterified FAs. After derivatization reaction, the detection sensitivity of FAs was remarkably improved and the limit of quantitation was lower than 100 ng/L. The quantitative errors caused by matrix effects were diminished benefiting from isotope-derivatized internal standards. The established quantitation strategy was successfully applied to verify both free and esterified FA contents in meat after different post-harvest procedures, finding that free polyunsaturated FAs increased significantly during freezing process.

Genome-wide association study revealed ABCD4 on SSC7 and GREB1L and MIB1 on SSC6 as crucial candidate genes for rib number in Beijing Black pigs.

As one of the few animals with variation in the number of rib pairs (RIB), the pig is a good model to study the mechanism of RIB regulation. Quantitative trait loci (QTL) for porcine RIB are present on Sus scrofa chromosome 7 (SSC7). Although several candidate genes exist in this QTL region on SSC7, the causal gene has yet to be verified. Beijing Black pig with 14-17 RIB is a good population for candidate gene mining and 1104 individuals were genotyped using the Illumina Porcine 50K BeadChip. A total of 14 SNPs from 95.49 to 97.78 Mb on SSC7 showed genome-wide significant association with RIB. On SSC7, a locuszoom plot using pairwise linkage disequilibrium displayed the narrowest linkage region encompassing only two genes, ABCD4 and VRTN. In mice, a transcriptome expression profile was obtained using three embryos at E9.5 (the critical period for rib formation). ABCD4 was highly expressed, but no expression of VRTN was detected. On SSC6, there were four genome-wide significant SNPs from 106.42 to 106.92 Mb associated with RIB. GREB1L and MIB1, in this region, were regarded as novel candidate genes. These results revealed a crucial candidate causal gene on SSC7 and novel genes on SSC6 for rib number and provided interesting new insights into its genetic basis.

Copy Number Variations Contribute to Intramuscular Fat Content Differences by Affecting the Expression of PELP1 Alternative Splices in Pigs.

Intramuscular fat (IMF) is a key meat quality trait. Research on the genetic mechanisms of IMF decomposition is valuable for both pork quality improvement and the treatment of obesity and type 2 diabetes. Copy number variations (CNVs) are a type of variant that may influence meat quality. In this study, a total of 1185 CNV regions (CNVRs) including 393 duplicated CNVRs, 432 deleted CNVRs, and 361 CNVRs with both duplicated and deleted status were identified in a pig F2 resource population using next-generation sequencing data. A genome-wide association study (GWAS) was then performed between CNVs and IMF, and a total of 19 CNVRs were found to be significantly associated with IMF. QTL colocation analysis indicated that 3 of the 19 CNVRs overlapped with known QTLs. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3'UTR end of the proline, as well as glutamate and the leucine rich protein 1 (PELP1) gene may affect the expression of PELP1 alternative splices. Sequence alignment and Alphafold2 structure prediction results indicated that the two alternative splices of PELP1 have a 23 AA sequence variation and a helix-fold structure variation. This region is located in the region of interaction between PELP1 and other proteins which have been reported to be significantly associated with fat deposition or insulin resistance. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affects the structure of the PELP1 protein. In conclusion, we found some CNVRs, especially CNV150, located in PELP1 that affect IMF. These findings suggest a novel mechanistic approach for meat quality improvement in animals and the potential treatment of insulin resistance in human beings.

Optimizing the Construction and Update Strategies for the Genomic Selection of Pig Reference and Candidate Populations in China.

Optimizing the construction and update strategies for reference and candidate populations is the basis of the application of genomic selection (GS). In this study, we first simulated1200-purebred-pigs population that have been popular in China for 20 generations to study the effects of different population sizes and the relationship between individuals of the reference and candidate populations. The results showed that the accuracy was positively correlated with the size of the reference population within the same generation (r = 0.9366, p < 0.05), while was negatively correlated with the number of generation intervals between the reference and candidate populations (r = -0.9267, p < 0.01). When the reference population accumulated more than seven generations, the accuracy began to decline. We then simulated the population structure of 1200 purebred pigs for five generations and studied the effects of different heritabilities (0.1, 0.3, and 0.5), genotyping proportions (20, 30, and 50%), and sex ratios on the accuracy of the genomic estimate breeding value (GEBV) and genetic progress. The results showed that if the proportion of genotyping individuals accounts for 20% of the candidate population, the traits with different heritabilities can be genotyped according to the sex ratio of 1:1male to female. If the proportion is 30% and the traits are of low heritability (0.1), the sex ratio of 1:1 male to female is the best. If the traits are of medium or high heritability, the male-to-female ratio is 1:1, 1:2, or 2:1, which may achieve higher genetic progress. If the genotyping proportion is up to 50%, for low heritability traits (0.1), the proportion of sows from all genotyping individuals should not be less than 25%, and for the medium and high heritability traits, the optimal choice for the male-to-female ratio is 1:1, which may obtain the greatest genetic progress. This study provides a reference for determining a construction and update plan for the reference population of breeding pigs.

Integrative Analysis of Nanopore and Illumina Sequencing Reveals Alternative Splicing Complexity in Pig Longissimus Dorsi Muscle.

Alternative splicing (AS) is a key step in the post-transcriptional regulation of gene expression that can affect intramuscular fat (IMF). In this study, longissimus dorsi muscles from 30 pigs in high- and low- IMF groups were used to perform Oxford Nanopore Technologies (ONT) full-length sequencing and Illumina strand-specific RNA-seq. A total of 43,688 full-length transcripts were identified, with 4,322 novel genes and 30,795 novel transcripts. Using AStalavista, a total of 14,728 AS events were detected in the longissimus dorsi muscle. About 17.79% of the genes produced splicing isoforms, in which exon skipping was the most frequent AS event. By analyzing the expression differences of mRNAs and splicing isoforms, we found that differentially expressed mRNAs with splicing isoforms could participate in skeletal muscle development and fatty acid metabolism, which might determine muscle-related traits. SERBP1, MYL1, TNNT3, and TNNT1 were identified with multiple splicing isoforms, with significant differences in expression. AS events occurring in IFI6 and GADD45G may cause significant differences in gene expression. Other AS events, such as ONT.15153.3, may regulate the function of ART1 by regulating the expression of different transcripts. Moreover, co-expression and protein-protein interaction (PPI) analysis indicated that several genes (MRPL27, AAR2, PYGM, PSMD4, SCNM1, and HNRNPDL) may be related to intramuscular fat. The splicing isoforms investigated in our research provide a reference for the study of alternative splicing regulation of intramuscular fat deposition.

Identification of imprinted genes in the skeletal muscle of newborn piglets by high-throughput sequencing.

Imprinted genes - exhibiting parent-specific transcription - play essential roles in the process of mammalian development and growth. Skeletal muscle growth is crucial for meat production. To further understand the role of imprinted genes during the porcine skeletal muscle growth, DNA-seq and RNA-seq were used to explore the characteristics of imprinted genes from porcine reciprocal crosses. A total of 584 545 single-nucleotide variations were discovered in the DNA-seq data of F0 parents, heterozygous in two pig breeds (Yorkshire and Min pigs) but homozygous in each breed. These single-nucleotide variations were used to determine the allelic-specific expression in F1 individuals. Finally, eight paternal expression sites and three maternal expression sites were detected, whereas two paternally expressed imprinted genes (NDN and IGF2) and one maternally expressed imprinted gene (H1-3) were validated by Sanger sequencing. DNA methylation regulates the expression of imprinted genes, and all of the identified imprinted genes in this study were predicted to possess CpG islands. PBX1 and YY1 binding motifs were discovered in the promoter regions of all three imprinted genes, which were candidate elements regulating the transcription of imprinted genes. For these identified imprinted genes, IGF2 and NDN promoted muscle growth whereas H1-3 inhibited cell proliferation, corroborating the 'parental conflict' theory that paternally expressed imprinted genes assisted descendants' growth whereas maternally expressed imprinted genes inhibited it. This study discovered porcine imprinted genes in skeletal muscle and was the first to reveal that H1-3 was expressed by the maternal allele to our knowledge. Our findings provided valuable resources for the potential utilization of imprinted genes in pig breeding.

Identifying long non-coding RNAs and characterizing their functional roles in swine mammary gland from colostrogenesis to lactogenesis.

OBJECTIVE: This study was conducted to identify the functional long non-coding RNAs (lncRNAs) for swine lactation by RNA-seq data of mammary gland. METHODS: According to the RNA-seq data of swine mammary gland, we screened lncRNAs, performed differential expression analysis, and confirmed the functional lncRNAs for swine lactation by validation of genome wide association study (GWAS) signals, functional annotation and weighted gene co-expression network analysis (WGCNA). RESULTS: We totally identified 286 differentially expressed (DE) lncRNAs in mammary gland at different stages from 14 days prior to (-) parturition to day 1 after (+) parturition, and the expressions of most of lncRNAs were strongly changed from day -2 to day +1. Further, the GWAS signals of sow milk ability trait were significantly enriched in DE lncRNAs. Functional annotation revealed that these DE lncRNAs were mainly involved in mammary gland and lactation developing, milk composition metabolism and colostrum function. By performing weighted WGCNA, we identified 7 out of 12 lncRNA-mRNA modules that were highly associated with the mammary gland at day -14, day -2, and day +1, in which, 35 lncRNAs and 319 mRNAs were involved. CONCLUSION: This study suggested that 18 lncRNAs and their 20 target genes were promising candidates for swine parturition and colostrum occurrence processes. Our research provided new insights into lncRNA profiles and their regulating mechanisms from colostrogenesis to lactogenesis in swine.

Genome-Wide Expression Profiling of mRNAs, lncRNAs and circRNAs in Skeletal Muscle of Two Different Pig Breeds.

RNA-Seq technology is widely used to analyze global changes in the transcriptome and investigate the influence on relevant phenotypic traits. Beijing Black pigs show differences in growth rate and meat quality compared to western pig breeds. However, the molecular mechanisms responsible for such phenotypic differences remain unknown. In this study, longissimus dorsi muscles from Beijing Black and Yorkshire pigs were used to construct RNA libraries and perform RNA-seq. Significantly different expressions were observed in 1051 mRNAs, 322 lncRNAs, and 82 circRNAs. GO and KEGG pathway annotation showed that differentially expressed mRNAs participated in skeletal muscle development and fatty acid metabolism, which determined the muscle-related traits. To explore the regulatory role of lncRNAs, the cis and trans-target genes were predicted and these lncRNAswere involved in the biological processes related to skeletal muscle development and fatty acid metabolismvia their target genes. CircRNAs play a ceRNA role by binding to miRNAs. Therefore, the potential miRNAs of differentially expressed circRNAs were predicted and interaction networks among circRNAs, miRNAs, and key regulatory mRNAs were constructed to illustrate the function of circRNAs underlying skeletal muscle development and fatty acid metabolism. This study provides new clues for elucidating muscle phenotypic variation in pigs.

Identification of SNPs and Candidate Genes for Milk Production Ability in Yorkshire Pigs.

Sow milk production ability is an important limiting factor impacting suboptimal growth and the survival of piglets. Through pig genetic improvement, litter sizes have been increased. Larger litters need more suckling mammary glands, which results in increased milk from the lactating sow. Hence, there is much significance to exploring sow lactation performance. For milk production ability, it is not practical to directly measure the milk yield, we used litter weight gain (LWG) throughout sow lactation as an indicator. In this study, we estimated the heritability of LWG, namely, 0.18 ± 0.07. We then performed a GWAS, and detected seven significant SNPs, namely, Sus scrofa Chromosome (SSC) 2: ASGA0010040 (p = 7.73E-11); SSC2:MARC0029355 (p = 1.30E-08), SSC6: WU_10.2_6_65751151 (p = 1.32E-10), SSC7: MARC0058875 (p = 4.99E-09), SSC10: WU_10.2_10_49571394 (p = 6.79E-08), SSC11: M1GA0014659 (p = 1.19E-07), and SSC15: MARC0042106 (p = 1.16E-07). We performed the distribution of phenotypes corresponding to the genotypes of seven significant SNPs and showed that ASGA0010040, MARC0029355, MARC0058875, WU_10.2_10_49571394, M1GA0014659, and MARC0042106 had extreme phenotypic values that corresponded to the homozygous genotypes, while the intermediate values corresponded to the heterozygous genotypes. We screened for flanking regions ± 200 kb nearby the seven significant SNPs, and identified 38 genes in total. Among them, 28 of the candidates were involved in lactose metabolism, colostrum immunity, milk protein, and milk fat by functional enrichment analysis. Through the combined analysis between 28 candidate genes and transcriptome data of the sow mammary gland, we found nine commons (ANO3, MUC15, DISP3, FBXO6, CLCN6, HLA-DRA, SLA-DRB1, SLA-DQB1, and SLA-DQA1). Furthermore, by comparing the chromosome positions of the candidate genes with the quantitative trait locus (QTLs) as previously reported, a total of 17 genes were found to be within 0.86-94.02 Mb of the reported QTLs for sow milk production ability, in which, NAV2 was found to be located with 0.86 Mb of the QTL region ssc2: 40936355. In conclusion, we identified seven significant SNPs located on SSC2, 6, 7, 10, 11, and 15, and propose 28 candidate genes for the ability to produce milk in Yorkshire pigs, 10 of which were key candidates.

A Comparison of the Behavior, Physiology, and Offspring Resilience of Gestating Sows When Raised in a Group Housing System and Individual Stalls.

Being in a confined environment causes chronic stress in gestating sows, which is detrimental for sow health, welfare and, consequently, offspring physiology. This study assessed the health and welfare of gestating sows housed in a group housing system compared to individual gestation stalls. After pregnancy was confirmed, experimental sows were divided randomly into two groups: the group housing system (GS), with the electronic sow feeding (ESF) system; or individual stall (IS). The behavior of sows housed in the GS or IS was then compared; throughout pregnancy, GS sows displayed more exploratory behavior, less vacuum chewing, and less sitting behavior (p < 0.05). IS sows showed higher stress hormone levels than GS sows. In particular, at 41 days of gestation, the concentration of the adrenocorticotropic hormone (ACTH) and adrenaline (A) in IS sows was significantly higher than that of GS sows, and the A level of IS sows remained significantly higher at 71 days of gestation (p < 0.01). The lipopolysaccharide (LPS) test was carried out in the weaned piglets of the studied sows. Compared with the offspring of gestating sows housed in GS (PG) or IS (PS), PG experienced a shorter period of high temperature and showed a quicker return to the normal state (p < 0.05). Additionally, their lower levels of stress hormone (p < 0.01) suggest that PG did not suffer from as much stress as PS. These findings suggested that gestating sows housed in GS were more able to carry out their natural behaviors and, therefore, had lower levels of stress and improved welfare. In addition, PG also showed better disease resistance and resilience. These results will provide a research basis for the welfare and breeding of gestating sows.