Background: Reproductive performance of livestock is an economically important aspect of global food production. The Chinese Meishan pig is a prolific breed, with an average of three to five more piglets per litter than European breeds; however, the genetic basis for this difference is not well understood. Results: In this study, we investigated copy number variations (CNVs) of 32 Meishan pigs and 29 Duroc pigs by next-generation sequencing. A genome-wide analysis of 61 pigs revealed 12,668 copy number variable regions (CNVRs) that were further divided into three categories based on copy number (CN) of the whole population, i.e., gain (n = 7,638), and loss (n = 5,030) CNVRs. We then compared Meishan and Duroc pigs and identified 17.17 Mb of 6,387 CNVRs that only existing in Meishan pigs CNVRs that overlapped the reproduction-related gene encoding the aryl hydrocarbon receptor (AHR) gene. We found that normal AHR CN was more frequent than CN loss in four different pig breeds. An association analysis showed that AHR CN had a positive effect on litter size (P < 0.05) and that a higher CN was associated with higher total number born (P < 0.05), number born alive (P < 0.05), number of weaned piglets, and birth weight. Conclusions: The present study provides comprehensive CNVRs for Meishan and Duroc pigs through large-scale population resequencing. Our results provide a supplement for the high-resolution map of copy number variation in the porcine genome and valuable information for the investigation of genomic structural variation underlying traits of interest in pig. In addition, the association results provide evidence for AHR as a candidate gene associated with reproductive traits that can be used as a genetic marker in pig breeding programs.

Pseudo-phenotypes, such as 305-day yields, estimated breeding values or deregressed proofs, are usually used as response variables for genome-wide association studies (GWAS) of milk production traits in dairy cattle. Computational inefficiency challenges the direct use of test-day records for longitudinal GWAS with large datasets. We propose a rapid longitudinal GWAS method that is based on a random regression model. Our method uses Eigen decomposition of the phenotypic covariance matrix to rotate the data, thereby transforming the complex mixed linear model into weighted least squares analysis. We performed a simulation study that showed that our method can control type I errors well and has higher power than a longitudinal GWAS method that does not include time-varied additive genetic effects. We also applied our method to the analysis of milk production traits in the first three parities of 6711 Chinese Holstein cows. The analysis for each trait was completed within 1 day with known variances. In total, we located 84 significant single nucleotide polymorphisms (SNPs) of which 65 were within previously reported quantitative trait loci (QTL) regions. Our rapid method can control type I errors in the analysis of longitudinal data and can be applied to other longitudinal traits. We detected QTL that were for the most part similar to those reported in a previous study in Chinese Holstein. Moreover, six additional SNPs for fat percentage and 13 SNPs for protein percentage were identified by our method. These additional 19 SNPs could be new candidate quantitative trait nucleotides for milk production traits in Chinese Holstein.

(1) Background: Weaning is a challenging and stressful event in the pig’s life, which disrupts physiological balance and induces oxidative stress. Microbiota play a significant role during the weaning process in piglets. Therefore, this study aimed to investigate key gut microbiota and metabolites associated with weaning stress in piglets. (2) Methods: A total of ten newborn piglet littermates were randomly assigned to two groups: S (suckling normally) and W (weaned at 21 d; all euthanized at 23 d). Specimens of the cecum were dehydrated with ethanol, cleared with xylene, embedded in paraffin, and cut into 4 mm thick serial sections. After deparaffinization, the sections were stained with hematoxylin and eosin (H&E) for morphometric analysis. Cecal metagenomic and liver LC-MS-based metabolomics were employed in this study. Statistical comparisons were performed by a two-tailed Student’s t-test, and p < 0.05 indicated statistical significance. (3) Results: The results showed that weaning led to intestinal morphological damage in piglets. The intestinal villi of suckling piglets were intact, closely arranged in an orderly manner, and finger-shaped, with clear contours of columnar epithelial cells. In contrast, the intestines of weaned piglets showed villous atrophy and shedding, as well as mucosal bleeding. Metagenomics and metabolomics analyses showed significant differences in composition and function between suckling and weaned piglets. The W piglets showed a decrease and increase in the relative abundance of Bacteroidetes and Proteobacteria (p < 0.05), respectively. The core cecal flora in W piglets were Campylobacter and Clostridium, while those in S piglets were Prevotella and Lactobacillus. At the phylum level, the relative abundance of Bacteroidetes significantly decreased (p < 0.05) in weaned piglets, while Proteobacteria significantly increased (p < 0.05). Significant inter-group differences were observed in pathways and glycoside hydrolases in databases, such as the KEGG and CAZymes, including fructose and mannose metabolism, salmonella infection, antifolate resistance, GH135, GH16, GH32, and GH84. We identified 757 differential metabolites between the groups through metabolomic analyses—350 upregulated and 407 downregulated (screened in positive ion mode). In negative ion mode, 541 differential metabolites were identified, with 270 upregulated and 271 downregulated. Major differential metabolites included glycerophospholipids, histidine, nitrogen metabolism, glycine, serine, threonine, β-alanine, and primary bile acid biosynthesis. The significant differences in glycine, serine, and threonine metabolites may be potentially related to dysbiosis caused by weaning stress. Taken together, the identification of microbiome and metabolome signatures of suckling and weaned piglets has paved the way for developing health-promoting nutritional strategies, focusing on enhancing bacterial metabolite production in early life stages.

(1) Background: Piglet diarrhea is one of the most serious diseases in pigs and has brought great economic losses to the pig industry. Alteration of the gut microbiota is an important factor in the etiology of piglet diarrhea. Therefore, this study aimed to analyze the differences in the gut microbial structures and fecal metabolic profile between post-weaning diarrhea and healthy Chinese Wannan Black pigs. (2) Methods: An integrated approach of 16S rRNA gene sequencing combined with LC/MS-based metabolomics was employed in this study. (3) Results: We found an increase in the relative abundance of the bacterial genus Campylobacter and a decrease in phylum Bacteroidetes and the species Streptococcus gallolyticus subsp. macedonicus. (S. macedonicus) in piglet diarrhea. Meanwhile, obvious changes in the fecal metabolic profile of diarrheic piglets were also detected, particularly higher levels of polyamines (spermine and spermidine). Moreover, there were substantial associations between the disturbed gut microbiota and the altered fecal metabolites, especially a strong positive relationship between spermidine and Campylobacter. (4) Conclusions: These observations may provide novel insights into potential etiologies related to post-weaning diarrhea and further enhance our understanding of the role of gut microbiota in host homeostasis and in modulating gut microbial structure.

Intramuscular fat (IMF) content is vital for pork quality, serving an important role in economic performance in pig industry. Non-coding RNAs, with mRNAs, are involved in IMF deposition; however, their functions and regulatory mechanisms in porcine IMF remain elusive. This study assessed the whole transcriptome expression profiles of the Longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to identify genes implicated in porcine IMF adipogenesis and their regulatory functions. Hundreds of differentially expressed RNAs were found to be involved in fatty acid metabolic processes, lipid metabolism, and fat cell differentiation. Furthermore, combing co-differential expression analyses, we constructed competing endogenous RNAs (ceRNA) regulatory networks, showing crosstalk among 30 lncRNAs and 61 mRNAs through 20 miRNAs, five circRNAs and 11 mRNAs through four miRNAs, and potential IMF deposition-related ceRNA subnetworks. Functional lncRNAs and circRNAs (such as MSTRG.12440.1, ENSSSCT00000066779, novel_circ_011355, novel_circ_011355) were found to act as ceRNAs of important lipid metabolism-related mRNAs (LEP, IP6K1, FFAR4, CEBPA, etc.) by sponging functional miRNAs (such as ssc-miR-196a, ssc-miR-200b, ssc-miR10391, miR486-y). These findings provide potential regulators and molecular regulatory networks that can be utilized for research on IMF traits in pigs, which would aid in marker-assisted selection to improve pork quality.