Coat color is one of the most important characteristics for distinguishing Chinese indigenous pig breeds. In Wuzhishan pigs, the animals have black on the back and white on the abdomen. However, the molecular genetic basis of this phenotype is unclear. In this study, we used high-throughput RNA sequencing to compare expression profiles of coding and non-coding RNAs from white and black skin samples obtained from individual Wuzhishan pigs. The expression profiling revealed that 194 lncRNAs (long non-coding RNAs), 189 mRNAs (messenger RNAs), and 162 miRNAs (microRNAs) had significantly different levels of expression (|log₂ fold change| > 1, p-value < 0.05) in white and black skin. Compared to RNA levels in black skin, white skin had higher levels of expression of 185 lncRNAs, 181 mRNAs, and 23 miRNAs and lower levels of expression of 9 lncRNAs, 8 mRNAs, and 139 miRNAs. Functional analysis suggested that the differentially expressed transcripts are involved in biological processes such as melanin biosynthesis, pigmentation and tyrosine metabolism. Several key genes involved in melanogenesis, including MLANA, PMEL, TYR, TYRP1, DTC, TRPM1 and CAMK2A, had significantly different levels of expression in the two skin tissues. Potential lncRNA–miRNA–gene interactions were also examined. A total of 15 lncRNAs, 11 miRNAs and 7 genes formed 23 lncRNA–miRNA–gene pairs, suggesting that complex regulatory networks of coding and non-coding genes underlie the coat color trait in Wuzhishan pigs. Our study provides a foundation for understanding how lncRNA, miRNA and genes interact to regulate coat color in black-back/white-belly pigs. We also constructed lncRNA–miRNA–gene interaction networks to elucidate the complex molecular mechanisms underlying skin physiology and melanogenesis. The results extend our knowledge about the diversity of coat color among different domestic animals and provide a foundation for studying novel mechanisms that control coat color in Chinese indigenous pigs.

Castration is usually used to remove boar taint in commercial pork production, but the adipose accumulation was increased excessively, which affected the meat quality of pigs. Based on our previous study, secreted phosphoprotein 1 (SPP1) was significantly differentially expressed between castrated and intact male pigs. However, the role of SPP1 in regulating adipose growth and fat storage caused by castration is unknown. In this study, SPP1 was identified to inhibit adipogenesis by the expression of adipogenic markers PPARγ and FABP4 as well as Oil red staining assay during differentiation of porcine bone marrow mesenchymal stem cells (pBMSCs). Subsequently, testosterone was used to treat pBMSCs to simulate the androgen status of intact pigs. Compared with the control groups without testosterone, the SPP1 expression in the testosterone group was markedly increased in the late stage of pBMSCs differentiation. Furthermore, novel-miR-659 was predicted by TargetScan and miRDB to target SPP1 and verified through a dual-luciferase reporter assay. Oil Red O staining assay indicated that novel-miR-659 overexpression significantly promoted adipogenesis, whereas novel-miR-659 inhibition suppressed adipogenesis. The expressions of adipogenic markers PPARγ and FABP4 showed the same tendency. Taken together, our study found that the targeted interaction between novel-miR-659 and SPP1 is involved in regulation of fat deposition in castrated male pigs.

The porcine body length trait is an essential factor affecting meat production and reproductive performance. It is evident that the development/lengthening of individual vertebrae is one of the main reasons for increases in body length; however, the underlying molecular mechanism remains unclear. In this study, RNA-seq analysis was used to profile the transcriptome (lncRNA, mRNA, and miRNA) of the thoracic intervertebral cartilage (TIC) at two time points (1 and 4 months) during vertebral column development in Yorkshire (Y) and Wuzhishan pigs (W). There were four groups: 1- (Y1) and 4-month-old (Y4) Yorkshire pigs and 1- (W1) and 4-month-old (W4) Wuzhishan pigs. In total, 161, 275, 86, and 126 differentially expressed (DE) lncRNAs, 1478, 2643, 404, and 750 DE genes (DEGs), and 74,51, 34, and 23 DE miRNAs (DE miRNAs) were identified in the Y4 vs. Y1, W4 vs. W1, Y4 vs. W4, and Y1 vs. W1 comparisons, respectively. Functional analysis of these DE transcripts (DETs) demonstrated that they had participated in various biological processes, such as cellular component organization or biogenesis, the developmental process, the metabolic process, bone development, and cartilage development. The crucial bone development-related candidate genes NK3 Homeobox 2 (NKX3.2), Wnt ligand secretion mediator (WLS), gremlin 1 (GREM1), fibroblast growth factor receptor 3 (FGFR3), hematopoietically expressed homeobox (HHEX), (collagen type XI alpha 1 chain (COL11A1), and Wnt Family Member 16 (WNT16)) were further identified by functional analysis. Moreover, lncRNA, miRNA, and gene interaction networks were constructed; a total of 55 lncRNAs, 6 miRNAs, and 7 genes formed lncRNA–gene, miRNA–gene, and lncRNA–miRNA–gene pairs, respectively. The aim was to demonstrate that coding and non-coding genes may co-regulate porcine spine development through interaction networks. NKX3.2 was identified as being specifically expressed in cartilage tissues, and it delayed chondrocyte differentiation. miRNA-326 regulated chondrocyte differentiation by targeting NKX3.2. The present study provides the first non-coding RNA and gene expression profiles in the porcine TIC, constructs the lncRNA–miRNA–gene interaction networks, and confirms the function of NKX3.2 in vertebral column development. These findings contribute to the understanding of the potential molecular mechanisms regulating pig vertebral column development. They expand our knowledge about the differences in body length between different pig species and provide a foundation for future studies.

Despite identifying genes regulating the coat colour in Western pig breeds, the genetic basis of the coat colour in Chinese indigenous pigs is still not understood due to the diversity of indigenous breeds and their genetic differences from exotic pigs. In this study, 215 Wuzhishan pigs with three coat colour patterns (white, black, and black-back/white-belly) were used to conduct a genome-wide association analysis. We found that genes responsible for the coat colour in the Wuzhishan breed are located on chromosome 8. Ninety-seven genome-wide significant SNPs are related to the animal’s coat colour. Using a haplotype-sharing analysis, we narrowed the potential candidate region to a 10.1 Mb interval encompassing only one gene, RAPGEF2, which participates in the regulation of melanogenesis. Two additional candidate genes, PDGFRA and KIT, are located within 1 Mb of the genome-wide significant SNPs. Gene ontology analysis and literature mining suggest that these candidate genes are associated with the animal’s coat colour. mRNA expression results revealed that RAPGEF2 and PDGFRA had significantly higher expressions in black pigs than in white pigs and higher expressions in black skin than in white skin from the same black-back/white-belly pigs. These results suggest that RAPGEF2 and PDGFRA are potential candidate genes regulating the coat colour in Wuzhishan pigs. Interestingly, mutations of KIT (a gene duplication and a G to A substitution at the splicing site in intron 17) were detected in white Wuzhishan pigs but not in black-back/white-belly or black pigs, suggesting a close genetic relationship between white Wuzhishan pigs and Western white pig breeds. In summary, these results indicate that the expression of RAPGEF2 and PDGFRA may cause the coat colour variation by influencing the deposition of melanin, while the mutation of KIT causes the white coat colour. Our results may provide a theoretical basis for the breeding of white coat colour Wuzhishan pigs, and shed light on the complex genetic background of coat colour variations in indigenous Chinese pig breeds.

The screening of important candidate genes and the identification of genetic markers are important for molecular selection in the pig industry. The hematopoietically expressed homeobox (HHEX) gene plays an important role in embryonic development and organogenesis; however, the genetic variation and expression pattern of the porcine HHEX gene remains to be clarified. In this study, semiquantitative RT-PCR and immunohistochemistry results showed the specific expression of the HHEX gene in porcine cartilage tissues. A novel haplotype consisting of two SNPs rs80901185 (T > C) and rs80934526 (A > G) was detected in the promoter region of the HHEX gene. The expression of the HHEX gene was significantly higher in Yorkshire pigs (TA haplotype) than in Wuzhishan pigs (CG haplotype), and a population analysis showed that this haplotype was significantly associated with body length. An analysis subsequently revealed that the –586 to –1 bp region of the HHEX gene promoter showed the highest activity. Furthermore, we found that the activity of the TA haplotype was significantly higher than that of the CG haplotype by changing the potential binding of transcription factors YY1 and HDAC2. In summary, we conclude that the porcine HHEX gene may contribute to the breeding of pigs for body length traits.

MicroRNAs (miRNAs) constitute small regulatory molecules for a wide array of biological activities (18~24 nucleotides in length), including adipogenesis and adipose deposition. Their effect is, however, incompletely defined in inducing fat accumulation in castrated male pigs. Based on our study, four nine-times miRNAs were selected to examine their functions in adipose formation activities. In 3T3-L1 cells and backfat tissues of castrated and intact male pigs, miR-F4-C12 was identified as a factor in adipose development utilizing quantitative real-time PCR (qRT-PCR). Further, miR-F4-C12 was identified to promote fat development, suggesting that miR-F4-C12 was involved in adipogenesis. Moreover, phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) was proposed by the TargetScan, miRDB and starBase as a target of miR-F4-C12 and verified through a two-luciferase reporter assay. The over-expression of miR-F4-C12 dramatically decreases the PIK3R1 protein level in 3T3-L1 cells. The mRNA and protein levels of PIK3R1 in castrated pigs are reduced relative to intact pigs, providing further evidence that PIK3R1 is involved in regulating adipose accumulation. These results suggest that miR-F4-C12 involves adipose development and may regulate subcutaneous adipose tissue accumulation by targeting PIK3R1 in castrated male pigs.