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.

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.

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.

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.

Jinhua pigs have excellent meat quality and intramuscular fat content (IMF). Crossbreeding of Jinhua with Duroc pigs with high productivity was conducted to develop a new composite breed that possesses the beneficial characteristics of both breeds. The objective of this study was to estimate the crossbreeding parameters (additive breed, dominance, and epistatic loss effects) using the Kinghorn model and genetic parameters (heritability and genetic correlation) for carcass and meat quality traits by analysing the phenotypic data of Jinhua, Duroc, and their crossbred pigs. Backfat thickness at the thinnest point above the 9th to 13th thoracic vertebrae (BF), longissimus muscle area between the 4th and 5th thoracic vertebrae (LMA), meat shear force value (SFV), and IMF were measured. The additive breed effects were significant for all traits: 1.59 cm, -8.30 cm2 , -6.38 lb/cm2 , and 1.76% for BF, LMA, SFV, and IMF, respectively. The dominance effect was significant for LMA (7.41 cm2 ) and IMF (-2.46%), whereas the epistatic loss effect was significant for only LMA (-15.18 cm2 ). The estimated heritability values were high, ranging from 0.58 for IMF to 0.76 for LMA. A negative but non-significant genetic correlation of -0.11 was estimated between BF and IMF; however, previous studies have reported that the genetic correlation between these traits is moderately positive in modern western pigs. Our results imply that, with the estimation of crossbreeding and genetic parameters, genetic improvement could be implemented to produce a new composite breed with good meat quality and productivity, to meet Japanese market requirements, by crossbreeding Jinhua 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.
CONCLUSIONS: 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).

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.

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.

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.

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.