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.

CircRNAs play an important role in fat deposition, and testosterone-deficient boars exhibit significantly increased fat deposition; however, the mechanism by which testosterone regulates fat deposition through circRNAs remains unclear. In this study, circRNA-seq of backfat and abdominal fat from castrated and intact full-sib Yorkshire pigs was performed. The GO and KEGG enrichment analyses revealed that the host genes of the dorsal DE circRNAs were mainly involved in fatty acid transport, while in abdominal tissues, these genes were mainly involved in adipogenesis and inflammation. The interaction among sus_circPAPPA2, ssc-miR-2366 and GK was verified by dual fluorescence experiments and in porcine preadipocytes. The overexpression of sus_circPAPPA2 significantly inhibited the differentiation of preadipocytes. The expression of sus_circPAPPA2 was increased after adding 100 nM of testosterone, and preadipocyte differentiation was significantly inhibited. Testosterone can affect preadipocyte differentiation by upregulating the expression of sus_circPAPPA2, sponging miR-2366 and regulating the expression of genes, such as GK. These results indicate that testosterone can regulate the expression of adipocyte differentiation- and lipid metabolism-related genes by regulating the expression of circRNA, and ceRNA networks are different in the testosterone regulation of adipose deposition in different parts. This study provides basic data enhancing the understanding of the interaction between the hormone environment and mir-2366/GK to regulate trait performance in pigs.

Embryo loss is a major factor affecting profitability in the pig industry. Embryonic mortality occurs during peri-implantation and mid-late gestation in pigs. Previous investigations have shown that the embryo loss rate in Meishan pigs is significantly lower than in commercial breeds. Most studies have focused on embryonic mortality during early gestation, but little is known about losses during mid-late gestation. In this study, we performed a transcriptome analysis of endometrial tissue in mid-late gestation sows (gestation days 49 and 72) sampled from two breeds (Meishan (MS) and Duroc (DU)) that have different embryo loss rates. We identified 411, 1113, 697, and 327 differentially expressed genes, and 14, 36, 57, and 43 differentially expressed miRNAs in four comparisons (DU49 vs. DU72, DU49 vs. MS49, DU72 vs. MS72, and MS49 vs. MS72), respectively. Subsequently; seven differentially expressed mRNAs and miRNAs were validated using qPCR. Functional analysis suggested the differentially expressed genes and miRNAs target genes mainly involved in regulation of hormone levels, blood vessel development, developmental process involved in reproduction, embryonic placenta development, and the immune system. A network analysis of potential miRNA-gene interactions revealed that differentially expressed miRNAs in Meishan pigs are involved in the response to estradiol and oxygen levels, and affect angiogenesis and blood vessel development. The binding site on ssc-miR-503 for epidermal growth factor (EGF) and the binding site on ssc-miR-671-5p for estrogen receptor α (ESR1) were identified using a dual luciferase assay. The results of this study will enable further exploration of miRNA-mRNA interactions important in pig pregnancy and will help to uncover molecular mechanisms affecting embryonic mortality in pigs during mid-late gestation.

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.

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.