Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations.

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

A Comprehensive Analysis of the Myocardial Transcriptome in ZBED6-Knockout Bama Xiang Pigs.

The ZBED6 gene is a transcription factor that regulates the expression of IGF2 and affects muscle growth and development. However, its effect on the growth and development of the heart is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) can regulate genes at the epigenetic, transcriptional, and posttranscriptional levels and play an important role in the development of eukaryotes. To investigate the function of ZBED6 in the cardiac development of pigs, we constructed the expression profiles of mRNAs and lncRNAs in myocardial tissue obtained from Bama Xiang pigs in the ZBED6 knockout group (ZBED6-KO) and the wild-type group (ZBED6-WT). A total of 248 differentially expressed genes (DEGs) and 209 differentially expressed lncRNAs (DELs) were detected, and 105 potential cis target genes of DELs were identified. The functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) databases revealed two GO items related to muscle development by the cis target genes of DELs. Moreover, IGF2 was the direct target gene of ZBED6 by ChIP-PCR experiment. Our results explored the mechanism and expression profile of mRNAs and lncRNAs of ZBED6 gene knockout on myocardium tissue development, mining the key candidate genes in that process like IGF2.

Genome-Wide Association Studies, Runs of Homozygosity Analysis, and Copy Number Variation Detection to Identify Reproduction-Related Genes in Bama Xiang Pigs.

Litter size and teat number are economically important traits in the porcine industry. However, the genetic mechanisms influencing these traits remain unknown. In this study, we analyzed the genetic basis of litter size and teat number in Bama Xiang pigs and evaluated the genomic inbreeding coefficients of this breed. We conducted a genome-wide association study to identify runs of homozygosity (ROH), and copy number variation (CNV) using the novel Illumina PorcineSNP50 BeadChip array in Bama Xiang pigs and annotated the related genes in significant single nucleotide polymorphisms and common copy number variation region (CCNVR). We calculated the ROH-based genomic inbreeding coefficients (F ROH) and the Spearman coefficient between F ROH and reproduction traits. We completed a mixed linear model association analysis to identify the effect of high-frequency copy number variation (HCNVR; over 5%) on Bama Xiang pig reproductive traits using TASSEL software. Across eight chromosomes, we identified 29 significant single nucleotide polymorphisms, and 12 genes were considered important candidates for litter-size traits based on their vital roles in sperm structure, spermatogenesis, sperm function, ovarian or follicular function, and male/female infertility. We identified 9,322 ROHs; the litter-size traits had a significant negative correlation to F ROH. A total of 3,317 CNVs, 24 CCNVR, and 50 HCNVR were identified using cnvPartition and PennCNV. Eleven genes related to reproduction were identified in CCNVRs, including seven genes related to the testis and sperm function in CCNVR1 (chr1 from 311585283 to 315307620). Two candidate genes (NEURL1 and SH3PXD2A) related to reproduction traits were identified in HCNVR34. The result suggests that these genes may improve the litter size of Bama Xiang by marker-assisted selection. However, attention should be paid to deter inbreeding in Bama Xiang pigs to conserve their genetic diversity.

Non-targeted Metabolomics Reveals Metabolic Characteristics of Porcine Atretic Follicles.

Follicular atresia is one of the main factors limiting the reproductive power of domestic animals. At present, the molecular mechanisms involved in porcine follicular atresia at the metabolic level remain unclear. In this study, we divided the follicles of Bama Xiang pigs into healthy follicles (HFs) and atretic follicles (AFs) based on the follicle morphology. The expression of genes related to atresia in granulosa cells (GCs) and the concentration of hormones in the follicular fluid (FF) from HFs and AFs were detected. We then used liquid chromatography-mass spectrometry-based non-targeted metabolomic approach to analyze the metabolites in the FF from HFs and AFs. The results showed that the content of estradiol was significantly lower in AFs than in HFs, whereas that of progesterone was significantly higher in AFs than that in HFs. The expression of BCL2, VEGFA, and CYP19A1 was significantly higher in HFs than in AFs. In contrast, the expression of BAX and CASPASE3 was significantly lower in HFs. A total of 18 differential metabolites (DMs) were identified, including phospholipids, bioactive substances, and amino acids. The DMs were involved in 12 metabolic pathways, including arginine biosynthesis and primary bile acid biosynthesis. The levels of eight DMs were higher in the HF group than those in the AF group (p < 0.01), and those of 10 DMs were higher in the AF group than those in the HF group (p < 0.01). These findings indicate that the metabolic characteristics of porcine AFs are lower levels of lipids such as phospholipids and higher levels of amino acids and bile acids than those in HFs. Disorders of amino acid metabolism and cholic acid metabolism may contribute to porcine follicular atresia.

Analysis of long non-coding RNAs in skeletal muscle of Bama Xiang pigs in response to heat stress.

The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) molecules and predict their target genes related to muscle development and lipid metabolism in longissimus dorsi (LD) muscles of Bama Xiang pigs under constant heat stress. Ten male Bama Xiang pigs with an average initial body weight of 14 kg were randomly divided into control group (22°C) and heat stress (35 °C) group. The experiment lasted for 28 days. All the pigs were slaughtered at the end of the experiment, and LD muscles were collected for muscle quality analysis and transcriptome sequencing. Heat stress reduced meat quality of Bama Xiang pigs. lncRNAs in LD were identified systematically by deep RNA sequencing between the two groups. The results showed that 365 lncRNAs from the LD were identified, including 128 intergenic lncRNAs, 82 intronic lncRNAs, and 155 anti-sense lncRNAs. The differences lie in transcript of length, number of exons and wider size distribution, and expression level per KB fragment in three subtypes of lncRNAs. The three types of transposable elements coverage, including Line/L1, SINE/tRNA, and LTR/ERVL-MaLR, are the highest in mRNA and the three subtypes of lncRNAs in pigs. lncRNAs and mRNAs were different in comparison of features. The results predicted the target genes of the significant differentially expressed lncRNAs related to muscle development and lipid metabolism. This is the first study to expand the knowledge about muscle-related lncRNAs biology in Bama Xiang pigs under heat stress and will contribute to the development of alleviating the adverse effects of heat stress on pork quality targeting lncRNAs.

Comprehensive inbred variation discovery in Bama pigs using de novo assemblies.

The Bama mini-pig (BM pig) is an inbred strain of the Bama Xiang pig (BX pig) and an important animal model used for studying human diseases. The extremely long inbreeding period renders a clear distinction between the features of the BM and BX strains, such as in their metabolism and olfactory system. However, there is limited information about differences between BM and BX animals at the genomic level. In this study, we generated genome sequencing data and used the assembly-vs-assembly approach to evaluate the phenotypic variations caused by inbreeding in these strains. Moreover, we detected differential expression of mutant genes related to the phenotypes in BX and BM pigs. We sequenced the genome of the BX pig strain and performed a series of analyses to reveal the comprehensive inbred genetic variants between BX and BM pigs. Here, the 2.56-Gb draft genome assembly for the BX pig and an N50 contig length of approximately 11.87 kb is described, and an N50 scaffold length of approximately 99 kb and the variations in the BX pig genome were identified by comparison with the BM pig reference genome. There were 1,424,354 single nucleotide polymorphisms (SNPs), 2,961,891 insertions and deletions (indels), 13,772 structural variants (SVs), and 20,606 copy number variants (CNVs) identified in the BX genome. Functional annotation of SVs and CNVs showed that the genes (ADGRE2, GPR143, olfactory receptor 52B4-like, olfactory receptor 10H1-like and SHROOM2) with both SVs and CNVs were enriched in the most of all KEGG pathways and gene ontology (GO) terms of mutant genes. ADGRE2, GPR143 and SHROOM2 were both found to have significant higher expression levels in BX pigs than in BM pigs. In the contrary, the expressions of olfactory receptor 52B4-like and olfactory receptor 10H1-like were significant lower in BX pigs than in BM pigs. In conclusion, sequence analysis of the BX pig genome revealed that the genome structure of the two pig strains has considerable genomic variation that was caused by the long inbreeding period. Moreover, qRT-PCR analysis of the mutant genes displayed a significant distinction that may be associated with phenotypic differences between these pig strains.

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs.

The Bama Xiang pig (BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat. To uncover the genetic basis of the BMX phenotype, we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs, respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2NL, FRS2, KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.