Integrated analysis of the DNA methylome and RNA transcriptome during the development of skeletal muscle in Duroc pigs.

BACKGROUND: Skeletal muscle development plays a crucial role in yield and quality of pork; however, this process is influenced by various factors. In this study, we employed whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing to comprehensively investigate the longissimus dorsi muscle (LDM), aiming to identify key genes that impact the growth and development of Duroc pigs with different average daily gains (ADGs). RESULTS: Eight pigs were selected and divided into two groups based on ADGs: H (774.89 g) group and L (658.77 g) group. Each pair of the H and L groups were half-siblings. The results of methylation sequencing revealed 2631 differentially methylated genes (DMGs) involved in metabolic processes, signalling, insulin secretion, and other biological activities. Furthermore, a joint analysis was conducted on these DMGs and the differentially expressed genes (DEGs) obtained from transcriptome sequencing of the same individual. This analysis identified 316 differentially methylated and differentially expressed genes (DMEGs), including 18 DMEGs in promoter regions and 294 DMEGs in gene body regions. Finally, LPAR1 and MEF2C were selected as candidate genes associated with muscle development. Bisulfite sequencing PCR (BSP) and quantitative real-time PCR (qRT-PCR) revealed that the promoter region of LPAR1 exhibited significantly lower methylation levels (P < 0.05) and greater expression levels (P < 0.05) in the H group than in the L group. Additionally, hypermethylation was observed in the gene body region of MEF2C, as was a low expression level, in the H group (P < 0.05). CONCLUSIONS: These results suggest that the differences in the ADGs of Duroc pigs fed the same diet may be influenced by the methylation levels and expression levels of genes related to skeletal muscle development.

Association between mitochondrial DNA copy number and production traits in pigs.

Mitochondria are essential organelles in the regulation of cellular energetic metabolism. Mitochondrial DNA copy number (mtDNA_CN) can be used as a proxy for mitochondria number, size, and activity. The aims of our study are to evaluate the effect of mtDNA_CN and mitochondrial haploblocks on production traits in pigs, and to identify the genetic background of this cellular phenotype. We collected performance data of 234 pigs and extracted DNA from skeletal muscle. Whole-genome sequencing data was used to determine mtDNA_CN. We found positive correlations of muscle mtDNA_CN with backfat thickness at 207 d (+0.14; p-value = 0.07) and negative correlations with carcase loin thickness (-0.14; p-value = 0.03). Pigs with mtDNA_CN values below the lower quartile had greater loin thickness (+4.1 mm; p-value = 0.01) and lower backfat thickness (-1.1 mm; p-value = 0.08), which resulted in greater carcase lean percentage (+2.4%; p-value = 0.04), than pigs with mtDNA_CN values above the upper quartile. These results support the hypothesis that a reduction of mitochondrial activity is associated with greater feed efficiency. Higher mtDNA_CN was also positively correlated with higher meat ultimate pH (+0.19; p-value <0.01) but we did not observe significant difference for meat ultimate pH between the two groups with extreme mtDNA_CN. We found no association of the most frequent mitochondrial haploblocks with mtDNA_CN or the production traits, but several genomic regions that harbour potential candidate genes with functions related to mitochondrial biogenesis and homeostasis were associated with mtDNA_CN. These regions provide new insights into the genetic background of this cellular phenotype but it is still uncertain if such associations translate into noticeable effects on the production traits.

Identification of postnatal development dependent genes and proteins in porcine epididymis.

BACKGROUND: The epididymis is a highly regionalized tubular organ possesses vectorial functions of sperm concentration, maturation, transport, and storage. The epididymis-expressed genes and proteins are characterized by regional and developmental dependent pattern. However, a systematic and comprehensive insight into the postnatal development dependent changes in gene and protein expressions of porcine epididymis is still lacking. Here, the RNA and protein of epididymis of Duroc pigs at different postnatal development stages were extracted by using commercial RNeasy Midi kit and extraction buffer (7 M Urea, 2 M thiourea, 3% CHAPS, and 1 mM PMSF) combined with sonication, respectively, which were further subjected to transcriptomic and proteomic profiling. RESULTS: Transcriptome analysis indicated that 198 and 163 differentially expressed genes (DEGs) were continuously up-regulated and down-regulated along with postnatal development stage changes, respectively. Most of the up-regulated DEGs linked to functions of endoplasmic reticulum and lysosome, while the down-regulated DEGs mainly related to molecular process of extracellular matrix. Moreover, the following key genes INSIG1, PGRMC1, NPC2, GBA, MMP2, MMP14, SFRP1, ELN, WNT-2, COL3A1, and SPARC were highlighted. A total of 49 differentially expressed proteins (DEPs) corresponding to postnatal development stages changes were uncovered by the proteome analysis. Several key proteins ACSL3 and ACADM, VDAC1 and VDAC2, and KNG1, SERPINB1, C3, and TF implicated in fatty acid metabolism, voltage-gated ion channel assembly, and apoptotic and immune processes were emphasized. In the integrative network, the key genes and proteins formed different clusters and showed strong interactions. Additionally, NPC2, COL3A1, C3, and VDAC1 are located at the hub position in each cluster. CONCLUSIONS: The identified postnatal development dependent genes and proteins in the present study will pave the way for shedding light on the molecular basis of porcine epididymis functions and are useful for further studies on the specific regulation mechanisms responsible for epididymal sperm maturation.

Genome-wide detection of CNV regions and their potential association with growth and fatness traits in Duroc pigs.

BACKGROUND: In the process of pig breeding, the average daily gain (ADG), days to 100 kg (AGE), and backfat thickness (BFT) are directly related to growth rate and fatness. However, the genetic mechanisms involved are not well understood. Copy number variation (CNV), an important source of genetic diversity, can affect a variety of complex traits and diseases and has gradually been thrust into the limelight. In this study, we reported the genome-wide CNVs of Duroc pigs using SNP genotyping data from 6627 animals. We also performed a copy number variation region (CNVR)-based genome-wide association studies (GWAS) for growth and fatness traits in two Duroc populations. RESULTS: Our study identified 953 nonredundant CNVRs in U.S. and Canadian Duroc pigs, covering 246.89 Mb (~ 10.90%) of the pig autosomal genome. Of these, 802 CNVRs were in U.S. Duroc pigs with 499 CNVRs were in Canadian Duroc pigs, indicating 348 CNVRs were shared by the two populations. Experimentally, 77.8% of nine randomly selected CNVRs were validated through quantitative PCR (qPCR). We also identified 35 CNVRs with significant association with growth and fatness traits using CNVR-based GWAS. Ten of these CNVRs were associated with both ADG and AGE traits in U.S. Duroc pigs. Notably, four CNVRs showed significant associations with ADG, AGE, and BFT, indicating that these CNVRs may play a pleiotropic role in regulating pig growth and fat deposition. In Canadian Duroc pigs, nine CNVRs were significantly associated with both ADG and AGE traits. Further bioinformatic analysis identified a subset of potential candidate genes, including PDGFA, GPER1, PNPLA2 and BSCL2. CONCLUSIONS: The present study provides a necessary supplement to the CNV map of the Duroc genome through large-scale population genotyping. In addition, the CNVR-based GWAS results provide a meaningful way to elucidate the genetic mechanisms underlying complex traits. The identified CNVRs can be used as molecular markers for genetic improvement in the molecular-guided breeding of modern commercial pigs.

Protein expression profiles in Meishan and Duroc sows during mid-gestation reveal differences affecting uterine capacity, endometrial receptivity, and the maternal-fetal Interface.

BACKGROUND: Embryonic mortality is a major concern in the commercial swine industry and primarily occurs early in gestation, but also during mid-gestation (~ days 50-70). Previous reports demonstrated that the embryonic loss rate was significant lower in Meishan than in commercial breeds (including Duroc). Most studies have focused on embryonic mortality in early gestation, but little is known about embryonic loss during mid-gestation. RESULTS: In this study, protein expression patterns in endometrial tissue from Meishan and Duroc sows were examined during mid-gestation. A total of 2170 proteins were identified in both breeds. After statistical analysis, 70 and 114 differentially expressed proteins (DEPs) were identified in Meishan and Duroc sows, respectively. Between Meishan and Duroc sows, 114 DEPs were detected at day 49, and 98 DEPs were detected at day 72. Functional enrichment analysis revealed differences in protein expression patterns in the two breeds. Around half of DEPs were more highly expressed in Duroc at day 49 (DUD49), relative to DUD72 and Meishan at day 49 (MSD49). Many DEPs appear to be involved in metabolic process such as arginine metabolism. Our results suggest that the differences in expression affect uterine capacity, endometrial matrix remodeling, and maternal-embryo cross-talk, and may be major factors influencing the differences in embryonic loss between Meishan and Duroc sows during mid-gestation. CONCLUSIONS: Our data showed differential protein expression pattern in endometrium between Meishan and Duroc sows and provides insight into the development process of endometrium. These findings could help us further uncover the molecular mechanism involved in prolificacy.

Weighted single-step GWAS identified candidate genes associated with semen traits in a Duroc boar population.

BACKGROUND: In the pig production industry, artificial insemination (AI) plays an important role in enlarging the beneficial impact of elite boars. Understanding the genetic architecture and detecting genetic markers associated with semen traits can help in improving genetic selection for such traits and accelerate genetic progress. In this study, we utilized a weighted single-step genome-wide association study (wssGWAS) procedure to detect genetic regions and further candidate genes associated with semen traits in a Duroc boar population. Overall, the full pedigree consists of 5284 pigs (12 generations), of which 2693 boars have semen data (143,113 ejaculations) and 1733 pigs were genotyped with 50 K single nucleotide polymorphism (SNP) array. RESULTS: Results show that the most significant genetic regions (0.4 Mb windows) explained approximately 2%~ 6% of the total genetic variances for the studied traits. Totally, the identified significant windows (windows explaining more than 1% of total genetic variances) explained 28.29, 35.31, 41.98, and 20.60% of genetic variances (not phenotypic variance) for number of sperm cells, sperm motility, sperm progressive motility, and total morphological abnormalities, respectively. Several genes that have been previously reported to be associated with mammal spermiogenesis, testes functioning, and male fertility were detected and treated as candidate genes for the traits of interest: Number of sperm cells, TDRD5, QSOX1, BLK, TIMP3, THRA, CSF3, and ZPBP1; Sperm motility, PPP2R2B, NEK2, NDRG, ADAM7, SKP2, and RNASET2; Sperm progressive motility, SH2B1, BLK, LAMB1, VPS4A, SPAG9, LCN2, and DNM1; Total morphological abnormalities, GHR, SELENOP, SLC16A5, SLC9A3R1, and DNAI2. CONCLUSIONS: In conclusion, candidate genes associated with Duroc boars' semen traits, including the number of sperm cells, sperm motility, sperm progressive motility, and total morphological abnormalities, were identified using wssGWAS. KEGG and GO enrichment analysis indicate that the identified candidate genes were enriched in biological processes and functional terms may be involved into spermiogenesis, testes functioning, and male fertility.

Genome-wide association study reveals genetic loci and candidate genes for average daily gain in Duroc pigs.

OBJECTIVE: Average daily gain (ADG) is an important target trait of pig breeding programs. We aimed to identify single nucleotide polymorphisms (SNPs) and genomic regions that are associated with ADG in the Duroc pig population. METHODS: We performed a genome-wide association study involving 390 Duroc boars and by using the PorcineSNP60K Beadchip and two linear models. RESULTS: After quality control, we detected 3,5971 SNPs, which included seven SNPs that are significantly associated with the ADG of pigs. We identified six quantitative trait loci (QTL) regions for ADG. These QTLs included four previously reported QTLs on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, and SSC13, as well as two novel QTLs on SSC6 and SSC16. In addition, we selected six candidate genes (general transcription factor 3C polypeptide 5, high mobility group AT-hook 2, nicotinamide phosphoribosyltransferase, oligodendrocyte transcription factor 1, pleckstrin homology and RhoGEF domain containing G4B, and ENSSSCG00000031548) associated with ADG on the basis of their physiological roles and positional information. These candidate genes are involved in skeletal muscle cell differentiation, diet-induced obesity, and nervous system development. CONCLUSION: This study contributes to the identification of the casual mutation that underlies QTLs associated with ADG and to future pig breeding programs based on marker-assisted selection. Further studies are needed to elucidate the role of the identified candidate genes in the physiological processes involved in ADG regulation.

Detecting selection signatures between Duroc and Duroc synthetic pig populations using high-density SNP chip.

The development of high throughput genotyping techniques has facilitated the identification of selection signatures of pigs. The detection of genomic selection signals in a population subjected to differential selection pressures may provide insights into the genes associated with economically and biologically important traits. To identify genomic regions under selection, we genotyped 488 Duroc (D) pigs and 155 D × Korean native pigs (DKNPs) using the Porcine SNP70K BeadChip. By applying the FST and extended haplotype homozygosity (EHH-Rsb) methods, we detected genes under directional selection associated with growth/stature (DOCK7, PLCB4, HS2ST1, FBP2 and TG), carcass and meat quality (TG, COL14A1, FBXO5, NR3C1, SNX7, ARHGAP26 and DPYD), number of teats (LOC100153159 and LRRC1), pigmentation (MME) and ear morphology (SOX5), which are all mostly near or at fixation. These results could be a basis for investigating the underlying mutations associated with observed phenotypic variation. Validation using genome-wide association analysis would also facilitate the inclusion of some of these markers in genetic evaluation programs.

SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population.

BACKGROUND: The aim of the present study was to compare the power of single nucleotide polymorphism (SNP)-based genome-wide association study (GWAS) and haplotype-based GWAS for quantitative trait loci (QTL) detection, and to detect novel candidate genes affecting economically important traits in a purebred Duroc population comprising seven-generation pedigree. First, we performed a simulation analysis using real genotype data of this population to compare the power (based on the null hypothesis) of the two methods. We then performed GWAS using both methods and real phenotype data comprising 52 traits, which included growth, carcass, and meat quality traits. RESULTS: In total, 836 animals were genotyped using the Illumina PorcineSNP60 BeadChip and 14 customized SNPs from regions of known candidate genes related to the traits of interest. The power of SNP-based GWAS was greater than that of haplotype-based GWAS in a simulation analysis. In real data analysis, a larger number of significant regions was obtained by SNP-based GWAS than by haplotype-based GWAS. For SNP-based GWAS, 23 genome-wide significant SNP regions were detected for 17 traits, and 120 genome-wide suggestive SNP regions were detected for 27 traits. For haplotype-based GWAS, 6 genome-wide significant SNP regions were detected for four traits, and 11 genome-wide suggestive SNP regions were detected for eight traits. All genome-wide significant SNP regions detected by haplotype-based GWAS were located in regions also detected by SNP-based GWAS. Four regions detected by SNP-based GWAS were significantly associated with multiple traits: on Sus scrofa chromosome (SSC) 1 at 304 Mb; and on SSC7 at 35-39 Mb, 41-42 Mb, and 103 Mb. The vertnin gene (VRTN) in particular, was located on SSC7 at 103 Mb and was significantly associated with vertebrae number and carcass lengths. Mapped QTL regions contain some candidate genes involved in skeletal formation (FUBP3; far upstream element binding protein 3) and fat deposition (METTL3; methyltransferase like 3). CONCLUSION: Our results show that a multigenerational pig population is useful for detecting QTL, which are typically segregated in a purebred population. In addition, a novel significant region could be detected by SNP-based GWAS as opposed to haplotype-based GWAS.

Relationship between perilipin genes polymorphisms and growth, carcass and meat quality traits in pigs.

The perilipins (PLIN) belong to a family of structural proteins that play a role regulating intracellular lipid storage and mobilization. Here, PLIN1 and PLIN2 have been evaluated as candidate genes for growth, carcass and meat quality traits in pigs. A sample of 607 Duroc pigs were genotyped for two single-nucleotide polymorphisms, one in intron 2 of the PLIN1 gene (JN860199:g.173G>A) and the other at the 3' untranslated region of the PLIN2 gene (GU461317:g.98G>A). Using a Bayesian approach, we have been able to find evidence of additive, dominant and epistatic associations of the PLIN1 and PLIN2 polymorphisms with early growth rate and carcass length. However, the major effects were produced by the dominant A allele at the PLIN2 polymorphism, which also affected the carcass lean weight. Thus, pigs carrying an additional copy of the A allele at the g.98G>A PLIN2 polymorphism had a probability of at least 98% of producing carcasses with heavier lean weight (+0.41 kg) and ham weight (+0.10 kg). The results obtained indicate that the PLIN2 polymorphism could be a useful marker for lean growth. In particular, it may help to reduce the undesired negative correlated response in lean weight to selection for increased intramuscular fat content, a common scenario in some Duroc lines involved in the production of high quality pork products.

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).

Genetic evaluation for production and body size traits using different animal models in purebred-Duroc pigs.

Duroc pigs are popular crossbred terminal sires, and accurate assessment of genetic parameters in the population can help to rationalize breeding programmes. The principle aim of this study were to evaluate the genetic parameters of production (birth weight, BW; age at 115 kg, AGE; feed conversion ratio, FCR) and body size (body length, BL; body height, BH; front cannon circumference, FCC) traits of Duroc pigs. The second objective was to analyze the fit of different genetic assessment models. The variance components and correlations of BW (28,348 records), AGE (28,335 records), FCR (11,135 records), BL (31,544 records), BH (21,862 records), and FCC (14,684 records) traits were calculated by using DMU and AIREMLF90 from BLUPF90 package. In the common environment model, the heritability of BW, AGE, FCR, BL, BH, and FCC traits were 0.17 ± 0.014, 0.30 ± 0.019, 0.28 ± 0.024, 0.16 ± 0.013, 0.14 ± 0.017, and 0.081 ± 0.016, with common litter effect values of 0.25, 0.20, 0.18, 0.23, 0.19, and 0.16, respectively. According to the results of the Akaike information criterion (AIC) calculations, models with smaller AIC values have a better fit. We found that the common environment model with litter effects as random effects for estimating genetic parameters had a better fit. In this Model, the estimated genetic correlations between AGE with BW, FCR, BL, BH, and FCC traits were -0.28 (0.040), 0.76 (0.038), -0.71 (0.036), -0.44 (0.060), and -0.60 (0.073), respectively, with phenotypic correlations of -0.17, 0.52, -0.22, -0.13 and -0.24, respectively. In our analysis of genetic trends for six traits in the Duroc population from 2012 to 2021, we observed significant genetic trends for AGE, BL, and BH. Particularly noteworthy is the rapid decline in the genetic trend for AGE, indicating an enhancement in the pig's growth rate through selective breeding. Therefore, we believe that some challenging-to-select traits can benefit from the genetic correlations between traits. By selecting easily measurable traits, they can gain from synergistic selection effects, leading to genetic progress. Conducting population genetic parameter analysis can assist us in devising breeding strategies.

Identification and Functional Prediction of Circular RNAs Related to Growth Traits and Skeletal Muscle Development in Duroc pigs.

Porcine skeletal muscle is a highly heterogeneous tissue type, and the Longissimus Dorsi muscle (LDM), as the most economical and physiologically metabolized skeletal muscle in pigs, has always been the focus of research and improvement in pig molecular breeding. Circular RNA, as an important new member of regulatory non-coding RNA after microRNA, has become a frontier hot spot in life science research. This study aims to explore candidate circRNAs related to growth, meat quality, and skeletal muscle development among Duroc pigs with different average daily gain (ADG). Eight pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG), followed by RNA-Seq analysis to identify circRNAs. The results showed that backfat at 6-7 rib (BF) and Intramuscular fat (IMF) content in the H group was significantly lower than L group, but ribeye area (REA) in the H group was higher than in the L group. In RNA-seq, 296 Differentially expressed (DE) circRNAs (157 upregulated and 139 downregulated) were identified and exons flanking long introns are easier to circularize to produce circRNAs. Most of the DE circRNAs were enriched in Quantitative trait locus (QTL) regions related to meat quality and growth traits. In addition, a gene can produce one or more circRNA transcripts. It was also found that the source genes of DE circRNAs were enriched in MAPK, FoXO, mTOR, PI3K-Akt, and Wnt signaling pathways. The results showed that different ADG, carcass, and meat quality traits among half-sibling Duroc pigs with the same diet may be due to the DE circRNAs related to skeletal muscle growth and development.

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.

Genome-Wide Association Study Reveals Additive and Non-Additive Effects on Growth Traits in Duroc Pigs.

Growth rate plays a critical role in the pig industry and is related to quantitative traits controlled by many genes. Here, we aimed to identify causative mutations and candidate genes responsible for pig growth traits. In this study, 2360 Duroc pigs were used to detect significant additive, dominance, and epistatic effects associated with growth traits. As a result, a total number of 32 significant SNPs for additive or dominance effects were found to be associated with various factors, including adjusted age at a specified weight (AGE), average daily gain (ADG), backfat thickness (BF), and loin muscle depth (LMD). In addition, the detected additive significant SNPs explained 2.49%, 3.02%, 3.18%, and 1.96% of the deregressed estimated breeding value (DEBV) variance for AGE, ADG, BF, and LMD, respectively, while significant dominance SNPs could explain 2.24%, 13.26%, and 4.08% of AGE, BF, and LMD, respectively. Meanwhile, a total of 805 significant epistatic effects SNPs were associated with one of ADG, AGE, and LMD, from which 11 sub-networks were constructed. In total, 46 potential genes involved in muscle development, fat deposition, and regulation of cell growth were considered as candidates for growth traits, including CD55 and NRIP1 for AGE and ADG, TRIP11 and MIS2 for BF, and VRTN and ZEB2 for LMD, respectively. Generally, in this study, we detected both new and reported variants and potential candidate genes for growth traits of Duroc pigs, which might to be taken into account in future molecular breeding programs to improve the growth performance of pigs.

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.

Identification of Homozygous Regions With Adverse Effects on the Five Economic Traits of Duroc Pigs.

Runs of homozygosity (ROH) are widely used to estimate genomic inbreeding, which is linked to inbreeding depression on phenotypes. However, the adverse effects of specific homozygous regions on phenotypic characteristics are rarely studied in livestock. In this study, the 50 K SNP data of 3,770 S21 Duroc (American origin) and 2,096 S22 Duroc (Canadian origin) pigs were used to investigate the harmful ROH regions on five economic traits. The results showed that the two Duroc lines had different numbers and distributions of unfavorable ROHs, which may be related to the different selection directions and intensities between the two lines. A total of 114 and 58 ROH segments were found with significant adverse effects on the economic traits of S21 and S22 pigs, respectively. Serval pleiotropic ROHs were detected to reduce two or multiple phenotypic performances in two Duroc populations. Candidate genes in these shared regions were mainly related to growth, fertility, immunity, and fat deposition. We also observed that some ROH genotypes may cause opposite effects on different traits. This study not only enhances our understanding of the adverse effects of ROH on phenotypes, but also indicates that ROH information could be incorporated into breeding programs to estimate and control the detrimental effects of homozygous regions.

Effect of ambient temperature on average daily gain of pigs evaluated using public weather data and a plateau-linear regression model.

We performed a plateau-linear regression model analysis of the average daily gain (ADG) of pigs on daily average temperature at the end of performance testing (T). Records for performance testing between 30 kg and 105 kg of 2268 purebred Duroc pigs raised at the National Livestock Breeding Center Miyazaki Station were used. Off-farm ambient temperatures were measured at the nearest Automated Meteorological Data Acquisition System station at Kobayashi, Miyazaki (Kobayashi station). A plateau-linear regression equation was obtained in which ADG decreased by 12.6 g for every 1°C when T > 21.1°C. We calculated the expected age in day at the end of testing (D105) using the regression equation obtained and T observed at the Kobayashi station in 2020. The number of days that D105 was prolonged due to higher T was 125 days, corresponding to approximately one third of the year. These results could contribute to planning and management of stable pork production in response to heat in Japan.

Heritability and genetic correlation estimates of semen production traits with litter traits and pork production traits in purebred Duroc pigs.

We estimated heritabilities of semen production traits and their genetic correlations with litter traits and pork production traits in purebred Duroc pigs. Semen production traits were semen volume, sperm concentration, proportion of morphologically normal sperms, total number of sperm, and total number of morphologically normal sperm. Litter traits at farrowing were total number born, number born alive, number stillborn, total litter weight at birth, mean litter weight at birth, and piglet survival rate at birth. Litter traits at weaning were litter size at weaning, total litter weight at weaning, mean litter weight at weaning, and piglet survival rate from birth to weaning. Pork production traits were average daily gain, backfat thickness, and loin muscle area. We analyzed 45,913 semen collection records of 896 boars, 6,950 farrowing performance records of 1,400 sows, 2,237 weaning performance records of 586 sows, and individual growth performance records of 9,550 animals measured at approximately 5 mo of age. Heritabilities were estimated using a single-trait animal model. Genetic correlations were estimated using a 2-trait animal model. Estimated heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of estimated genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Genetic parameters of 5 semen production traits, 10 litter traits, and 3 pork production traits in purebred Duroc pigs was estimated. Heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Identification of new semen trait-related candidate genes in Duroc boars through genome-wide association and weighted gene co-expression network analyses.

Semen traits are crucial in commercial pig production since semen from boars is widely used in artificial insemination for both purebred and crossbred pig production. Revealing the genetic architecture of semen traits potentially promotes the efficiencies of improving semen traits through artificial selection. This study is aimed to identify candidate genes related to the semen traits in Duroc boars. First, we identified the genes that were significantly associated with three semen traits, including sperm motility (MO), sperm concentration (CON), and semen volume (VOL) in a Duroc boar population through a genome-wide association study (GWAS). Second, we performed a weighted gene co-expression network analysis (WGCNA). A total of 2, 3, and 20 single-nucleotide polymorphisms were found to be significantly associated with MO, CON, and VOL, respectively. Based on the haplotype block analysis, we identified one genetic region associated with MO, which explained 6.15% of the genetic trait variance. ENSSSCG00000018823 located within this region was considered as the candidate gene for regulating MO. Another genetic region explaining 1.95% of CON genetic variance was identified, and, in this region, B9D2, PAFAH1B3, TMEM145, and CIC were detected as the CON-related candidate genes. Two genetic regions that accounted for 2.23% and 2.48% of VOL genetic variance were identified, and, in these two regions, WWC2, CDKN2AIP, ING2, TRAPPC11, STOX2, and PELO were identified as VOL-related candidate genes. WGCNA analysis showed that, among these candidate genes, B9D2, TMEM145, WWC2, CDKN2AIP, TRAPPC11, and PELO were located within the most significant module eigengenes, confirming these candidate genes' role in regulating semen traits in Duroc boars. The identification of these candidate genes can help to better understand the genetic architecture of semen traits in boars. Our findings can be applied for semen traits improvement in Duroc boars.