Genome-Wide Detection of Selection Signatures in Duroc Revealed Candidate Genes Relating to Growth and Meat Quality

Whole-genome SNP allele frequency differences between Tibetan and Large white pigs reveal genes associated with skeletal muscle growth

BACKGROUND

The skeletal muscle growth rate and body size of Tibetan pigs (TIB) are lower than Large white pigs (LW). However, the underlying genetic basis attributing to these differences remains uncertain. To address this knowledge gap, the present study employed whole-genome sequencing of TIB (slow growth) and LW (fast growth) individuals, and integrated with existing NCBI sequencing datasets of TIB and LW individuals, enabling the identification of a comprehensive set of genetic variations for each breed. The specific and predominant SNPs in the TIB and LW populations were detected by using a cutoff value of 0.50 for SNP allele frequency and absolute allele frequency differences (△AF) between the TIB and LW populations.

RESULTS

A total of 21,767,938 SNPs were retrieved from 44 TIB and 29 LW genomes. The analysis detected 2,893,106 (13.29%) and 813,310 (3.74%) specific and predominant SNPs in the TIB and LW populations, and annotated to 24,560 genes. Further GO analysis revealed 291 genes involved in biological processes related to striated and/or skeletal muscle differentiation, proliferation, hypertrophy, regulation of striated muscle cell differentiation and proliferation, and myoblast differentiation and fusion. These 291 genes included crucial regulators of muscle cell determination, proliferation, differentiation, and hypertrophy, such as members of the Myogenic regulatory factors (MRF) (MYOD, MYF5, MYOG, MYF6) and Myocyte enhancer factor 2 (MEF2) (MEF2A, MEF2C, MEF2D) families, as well as muscle growth inhibitors (MSTN, ACVR1, and SMAD1); KEGG pathway analysis revealed 106 and 20 genes were found in muscle growth related positive and negative regulatory signaling pathways. Notably, genes critical for protein synthesis, such as MTOR, IGF1, IGF1R, IRS1, INSR, and RPS6KA6, were implicated in these pathways.

CONCLUSION

This study employed an effective methodology to rigorously identify the potential genes associated with skeletal muscle development. A substantial number of SNPs and genes that potentially play roles in the divergence observed in skeletal muscle growth between the TIB and LW breeds were identified. These findings offer valuable insights into the genetic underpinnings of skeletal muscle development and present opportunities for enhancing meat production through pig breeding.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Persistence of autozygosity in crossbreds between autochthonous and cosmopolitan breeds of swine: a simulation study

Crossbreeding might be a valid strategy to valorize local pig breeds. Crossbreeding should reduce homozygosity and, as a consequence, yield hybrid vigor for fitness and production traits. This study aimed to quantify the persistence of autozygosity in terminal crossbred pigs compared with purebreds and, in turn, identify genomic regions where autozygosity's persistence would not be found. The study was based on genotyping data from 20 European local pig breeds and three cosmopolitan pig breeds used to simulate crossbred offspring. This study consisted of two steps. First, one hundred matings were simulated for each pairwise combination of the 23 considered breeds (for a total of 276 combinations), ignoring the sex of the parent individuals in order to generate purebred and crossbred matings leveraging all the germplasm available. Second, a few preselected terminal-maternal breed pairs were used to mimic a realistic terminal crossbreeding system: (i) Mora Romagnola (boars) or Cinta Senese (boars) crossed with Large White (sows) or Landrace (sows); (ii) Duroc (boars) crossed with Mora Romagnola (sows) or Cinta Senese (sows). Runs of homozygosity was used to estimate genome-wide autozygosity (FROH). Observed FROH was higher in purebreds than in crossbreds, although some crossbred combinations showed higher FROH than other purebred combinations. Among the purebreds, the highest FROH values were observed in Mora Romagnola and Turopolje (0.50 and 0.46, respectively). FROH ranged from 0.04 to 0.16 in the crossbreds Alentejana × Large White and Alentejana × Iberian, respectively. Persistence of autozygosity was found in several genomic segments harboring regions where quantitative trait loci (QTLs) were found in the literature. The regions were enriched in QTLs involved in fatty acid metabolism and associated with performance traits. This simulation shows that autozygosity persists in most breed combinations of terminal crosses. Results suggest that a strategy for crossbreeding is implemented when leveraging autochthonous and cosmopolitan breeds to obtain most of the hybrid vigor.

Expressions of Insulin-like Growth Factor System among Different Breeds Impact Piglets' Growth during Weaning

The present study investigated the insulin-like growth factors (IGFs) and their receptors and binding proteins among three pig breeds during weaning. Sixty Duroc (DR), Taoyuan black (TYB), and Xiangcun black (XCB) piglets (20 piglets per breed) were selected at 21 and 24 (3 days of post-weaning) days of age to analyze organ indices, plasma concentrations of IGF and IGF-binding proteins (IGFBPs) using ELISA kits, and gene expression of IGF-system-related components in different tissues. The plasma IGFBP-3 concentration in TYB piglets was higher (p > 0.05) than in the XCB and DR piglets at 21 days of age. At 21 days of age, compared with the DR piglets, the IGF-1 expression was lower (p < 0.05) in the kidney, but it was higher (p < 0.05) in the spleen of XCB and TYB piglets. At 24 days of age, the IGF-1 expression was higher (p < 0.05) in the kidney of TYB piglets than in the XCB and DR piglets, while IGFBP-3 in the stomach and IGFBP-4 in the liver of XCB and TYB piglets were lower (p < 0.05) compared with the DR piglets. Weaning down-regulated (p < 0.05) IGF-1 expression in the jejunum, spleen, and liver of piglets, while it up-regulated (p < 0.05) IGFBP-3 expression in the stomach, IGFBP-4 in the liver, IGFBP-5 in the ileum, and IGFBP-6 in the jejunum of DR piglets. Spearman's correlation analysis showed a negative correlation (p < 0.05) between plasma IGFBP-2 and IGFBP-5 concentration and the organ indices of piglets. Collectively, there were significant differences in the IGF system components among the three pig breeds. The IGF system components were altered during weaning, which might be involved in weaning stress to decrease the growth of piglets.

Genome-wide detection of CNV regions between Anqing six-end-white and Duroc pigs

BACKGROUND

Anqing six-end-white pig is a native breed in Anhui Province. The pigs have the disadvantages of a slow growth rate, low proportion of lean meat, and thick back fat, but feature the advantages of strong stress resistance and excellent meat quality. Duroc pig is an introduced pig breed with a fast growth rate and high proportion of lean meat. With the latter breed featuring superior growth characteristics but inferior meat quality traits, the underlying molecular mechanism that causes these phenotypic differences between Chinese and foreign pigs is still unclear.

RESULTS

In this study, copy number variation (CNV) detection was performed using the re-sequencing data of Anqing Six-end-white pigs and Duroc pigs, A total of 65,701 CNVs were obtained. After merging the CNVs with overlapping genomic positions, 881 CNV regions (CNVRs) were obtained. Based on the obtained CNVR information combined with their positions on the 18 chromosomes, a whole-genome map of the pig CNVs was drawn. GO analysis of the genes in the CNVRs showed that they were primarily involved in the cellular processes of proliferation, differentiation, and adhesion, and primarily involved in the biological processes of fat metabolism, reproductive traits, and immune processes.

CONCLUSION

The difference analysis of the CNVs between the Chinese and foreign pig breeds showed that the CNV of the Anqing six-end-white pig genome was higher than that of the introduced pig breed Duroc. Six genes related to fat metabolism, reproductive performance, and stress resistance were found in genome-wide CNVRs (DPF3, LEPR, MAP2K6, PPARA, TRAF6, NLRP4).

Genome-wide scans for selection signatures in indigenous pigs revealed candidate genes relating to heat tolerance

Heat stress is a major problem that constrains pig productivity. Understanding and identifying adaptation to heat stress has been the focus of recent studies, and the identification of genome-wide selection signatures can provide insights into the mechanisms of environmental adaptation. Here, we generated whole-genome re-sequencing data from six Chinese indigenous pig populations to identify genomic regions with selection signatures related to heat tolerance using multiple methods: three methods for intra-population analyses (Integrated Haplotype Score, Runs of Homozygosity and Nucleotide diversity Analysis) and three methods for inter-population analyses (Fixation index (F_ST), Cross-population Composite Likelihood Ratio and Cross-population Extended Haplotype Homozygosity). In total, 1 966 796 single nucleotide polymorphisms were identified in this study. Genetic structure analyses and F_ST indicated differentiation among these breeds. Based on information on the location environment, the six breeds were divided into heat and cold groups. By combining two or more approaches for selection signatures, outlier signals in overlapping regions were identified as candidate selection regions. A total of 163 candidate genes were identified, of which, 29 were associated with heat stress injury and anti-inflammatory effects. These candidate genes were further associated with 78 Gene Ontology functional terms and 30 Kyoto Encyclopedia of Genes and Genomes pathways in enrichment analysis (P < 0.05). Some of these have clear relevance to heat resistance, such as the AMPK signalling pathway and the mTOR signalling pathway. The results improve our understanding of the selection mechanisms responsible for heat resistance in pigs and provide new insights of introgression in heat adaptation.

Colicins of Escherichia coli Lead to Resistance against the Diarrhea-Causing Pathogen Enterotoxigenic E. coli in Pigs

Gut microbes can affect host adaptation to various environment conditions. Escherichia coli is a common gut species, including pathogenic strains and nonpathogenic strains. This study was conducted to investigate the effects of different E. coli strains in the gut on the health of pigs. In this study, the complete genomes of two E. coli strains isolated from pigs were sequenced. The whole genomes of Y18J and the enterotoxigenic E. coli strain W25K were compared to determine their roles in pig adaptation to disease. Y18J was isolated from feces of healthy piglets and showed strong antimicrobial activity against W25K in vitro. Gene knockout experiments and complementation analysis followed by modeling the microbe-microbe interactions demonstrated that the antagonistic mechanism of Y18J against W25K relied on the bacteriocins colicin B and colicin M. Compared to W25K, Y18J is devoid of exotoxin-coding genes and has more secondary-metabolite-biosynthetic gene clusters. W25K carries more genes involved in genome replication, in accordance with a shorter cell cycle observed during a growth experiment. The analysis of gut metagenomes in different pig breeds showed that colicins B and M were enriched in Laiwu pigs, a Chinese local breed, but were scarce in boars and Duroc pigs. IMPORTANCE This study revealed the heterogeneity of E. coli strains from pigs, including two strains studied by both in silico and wet experiments in detail and 14 strains studied by bioinformatics analysis. E. coli Y18J may improve the adaptability of pigs toward disease resistance through the production of colicins B and M. Our findings could shed light on the pathogenic and harmless roles of E. coli in modern animal husbandry, leading to a better understanding of intestinal-microbe-pathogen interactions in the course of evolution.

Transcriptome Analysis of the Adipose Tissue of Luchuan and Duroc Pigs

Simple Summary

Fat is a vital body tissue of pigs and a crucial index that affects the production efficiency of pigs. In this study, Duroc pigs and Luchuan pigs were used as animal materials, transcriptome sequencing was used to compare the back adipose tissue of the two breeds, to explore the key reason of difference in fat deposition. The result provided new ideas and reference for further study of fat development.

Abstract

Fat deposition is a crucial element in pig production that affects production efficiency, quality and consumer choices. In this study, Duroc pigs, a Western, famous lean pig breed, and Luchuan pigs, a Chinese, native obese pig breed, were used as animal materials. Transcriptome sequencing was used to compare the back adipose tissue of Duroc and Luchuan pigs, to explore the key genes regulating fat deposition. The results showed that 418 genes were highly expressed in the Duroc pig, and 441 genes were highly expressed in the Luchuan pig. In addition, the function enrichment analysis disclosed that the DEGs had been primarily enriched in lipid metabolism, storage and transport pathways. Furthermore, significant differences in the metabolic pathways of alpha-linolenic acid, linoleic acid and arachidonic acid explained the differences in the flavor of the two kinds of pork. Finally, the gene set enrichment analysis (GSEA) exposed that the difference in fat deposition between Duroc and Luchuan pigs may be due to the differential regulation of the metabolism pathway of fatty acid. Therefore, this study described the differential expression transcriptional map of adipose tissue of Duroc pig and Luchuan pig, identified the functional genes regulating pig fat deposition, and provided new hypotheses and references for further study of fat development.

Genome-Wide Association Analysis and Genetic Parameters for Feed Efficiency and Related Traits in Yorkshire and Duroc Pigs

Simple Summary

Genetic improvements in feed efficiency (FE) and related traits could considerably reduce pig production costs and energy consumption. Thus, we performed a genetic parameter estimation and genome-wide association study of four FE and FE-related traits, namely, average daily feed intake, average daily gain, the feed conversion ratio, and residual feed intake, of two pig breeds, Yorkshire and Duroc. The results demonstrate the genetic relationships of FE and FE-related traits with two growth traits, age and backfat thickness at 100 kg. We also identified many single-nucleotide polymorphisms (SNPs) and novel candidate genes related to these traits. In addition, we found many pathways significantly associated with FE and FE-related traits, and they are generally involved in digestive and metabolic processes. The results of this study are expected to provide a valuable reference for the genomic selection of FE and FE-related traits in pigs.

Abstract

Feed efficiency (FE) traits are key factors that can influence the economic benefits of pig production. However, little is known about the genetic architecture of FE and FE-related traits. This study aimed to identify SNPs and candidate genes associated with FE and FE-related traits, namely, average daily feed intake (ADFI), average daily gain (ADG), the feed conversion ratio (FCR), and residual feed intake (RFI). The phenotypes of 5823 boars with genotyped data (50 K BeadChip) from 1365 boars from a nucleus farm were used to perform a genome-wide association study (GWAS) of two breeds, Duroc and Yorkshire. Moreover, we performed a genetic parameter estimation for four FE and FE-related traits. The heritabilities of the FE and FE-related traits ranged from 0.13 to 0.36, and there were significant genetic correlations (−0.69 to 0.52) of the FE and FE-related traits with two growth traits (age at 100 kg and backfat thickness at 100 kg). A total of 61 significant SNPs located on eight different chromosomes associated with the four FE and FE-related traits were identified. We further identified four regions associated with FE and FE-related traits that have not been previously reported, and they may be potential novel QTLs for FE. Considering their biological functions, we finally identified 35 candidate genes relevant for FE and FE-related traits, such as the widely reported MC4R and INSR genes. A gene enrichment analysis showed that FE and FE-related traits were highly enriched in the biosynthesis, digestion, and metabolism of biomolecules. This study deepens our understanding of the genetic mechanisms of FE in pigs and provides valuable information for using marker-assisted selection in pigs to improve FE.

Identification and functional prediction of long non-coding RNAs related to skeletal muscle development in Duroc pigs

Objective

The growth of pigs involves multiple regulatory mechanisms, and modern molecular breeding techniques can be used to understand the skeletal muscle growth and development to promote the selection process of pigs. This study aims to explore candidate lncRNAs and mRNAs related to skeletal muscle growth and development among Duroc pigs with different average daily gain (ADG).

Methods

A total of 8 pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG). And followed by whole transcriptome sequencing to identify differentially expressed (DE) lncRNAs and mRNAs.

Results

In RNA-seq, 703 DE mRNAs (263 up-regulated and 440 down-regulated) and 74 DE lncRNAs (45 up-regulated and 29 down-regulated) were identified. In addition, 1,418 Transcription factors (TFs) were found. Compared with mRNAs, lncRNAs had fewer exons, shorter transcript length and open reading frame length. DE mRNAs and DE lncRNAs can form 417 lncRNA-mRNA pairs (antisense, cis and trans). DE mRNAs and target genes of lncRNAs were enriched in cellular processes, biological regulation, and regulation of biological processes. In addition, quantitative trait locus (QTL) analysis was used to detect the functions of DE mRNAs and lncRNAs, the most of DE mRNAs and target genes of lncRNAs were enriched in QTLs related to growth traits and skeletal muscle development. In single-nucleotide polymorphism/insertion-deletion (SNP/INDEL) analysis, 1,081,182 SNP and 131,721 INDEL were found, and transition was more than transversion. Over 60% of percentage were skipped exon events among alternative splicing events.

Conclusion

The results showed that different ADG among Duroc pigs with the same diet maybe due to the DE mRNAs and DE lncRNAs related to skeletal muscle growth and development.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia

BACKGROUND

The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS).

RESULTS

The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance.

CONCLUSIONS

Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

Signatures of selection reveal candidate genes involved in production traits in Chinese crossbred buffaloes

Identification of selection signature is important for a better understanding of genetic mechanisms that affect phenotypic differentiation in livestock. However, the genome-wide selection responses have not been investigated for the production traits of Chinese crossbred buffaloes. In this study, an SNP data set of 133 buffaloes (Chinese crossbred buffalo, n = 45; Chinese local swamp buffalo, n = 88) was collected from the Dryad Digital Repository database (https://datadryad.org/stash/). Population genetics analysis showed that these buffaloes were divided into the following 2 groups: crossbred buffalo and swamp buffalo. The crossbred group had higher genetic diversity than the swamp group. Using 3 complementary statistical methods (integrated haplotype score, cross population extended haplotype homozygosity, and composite likelihood ratio), a total of 31 candidate selection regions were identified in the Chinese crossbred population. Here, within these candidate regions, 25 genes were under the putative selection. Among them, several candidate genes were reported to be associated with production traits. In addition, we identified 13 selection regions that overlapped with bovine QTLs that were mainly involved in milk production and composition traits. These results can provide useful insights regarding the selection response for production traits of Chinese crossbred buffalo, as identified candidate genes influence production performance.

Genome-wide analysis of genetic diversity and artificial selection in Large White pigs in Russia

Breeding practices adopted at different farms are aimed at maximizing the profitability of pig farming. In this work, we have analyzed the genetic diversity of Large White pigs in Russia. We compared genomes of historic and modern Large White Russian breeds using 271 pig samples. We have identified 120 candidate regions associated with the differentiation of modern and historic pigs and analyzed genomic differences between the modern farms. The identified genes were associated with height, fitness, conformation, reproductive performance, and meat quality.

Lipidomic and Transcriptomic Analysis of the Longissimus Muscle of Luchuan and Duroc Pigs

Meat is an essential food, and pork is the largest consumer meat product in China and the world. Intramuscular fat has always been the basis for people to select and judge meat products. Therefore, we selected the Duroc, a western lean pig breed, and the Luchuan, a Chinese obese pig breed, as models, and used the longissimus dorsi muscle for lipidomics testing and transcriptomics sequencing. The purpose of the study was to determine the differences in intramuscular fat between the two breeds and identify the reasons for the differences. We found that the intramuscular fat content of Luchuan pigs was significantly higher than that of Duroc pigs. The triglycerides and diglycerides related to flavor were higher in Luchuan pigs compared to Duroc pigs. This phenotype may be caused by the difference in the expression of key genes in the glycerolipid metabolism signaling pathway.