Ningxiang pigs, a Chinese bred pig known for its tender meat and high quality unsaturated fatty acids. This study discovers the transcription profiles and functional networks in long non-coding RNA (lncRNA) and messenger RNA (mRNA) in subcutaneous adipose tissue. Subcutaneous adipose tissue was collected from piglet, nursery pig, early fattening, and late fattening stage of Ningxiang piglets, and lncRNA and mRNA transcription of each stage was profiled. A total of 339,204,926 (piglet), 315,609,246 (nursery), 266,798,202 (early fattening), and 343,740,308 (late fattening) clean reads were generated, and 2872 novel lncRNAs were identified. Additionally, 10,084 differential mRNAs (DEmRNAs) and 931 differential lncRNAs were determined. Most DEmRNAs were up-regulated in the piglet stage, while they were down-regulated in late fattening stage. A complicated interaction between mRNAs and lncRNAs was identified via STEM and WGCNA, demonstrated that lncRNAs are a significant regulatory component in mRNAs. The findings showed that modules 2 and 5 have a similar mode of transcription for both mRNA and lncRNA, and were mainly participated in steroid biosynthesis, glycosphingolipid biosynthesis, metabolic pathways, and glycerolipid metabolism. The mRNAs and lncRNAs transcription levels of both modules was higher in the early and late fattening stage, which may be due to the active activity of the metabolism in relation to fatty acids, sterols, steroids, and lipids in the subcutaneous adipose tissue during the early and late fattening stage. These findings could be expected to result in further research of the functional properties of lncRNA from subcutaneous adipose tissue at different stages of development in Ningxiang pigs.

Ningxiang pigs, a fat-type pig, are native to Ningxiang County in Hunan Province, with thousands of years of breeding history. This study aims to explore the expression profiles and functional networks on messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) in the liver. Liver tissue of Ningxiang piglets was collected at 30, 90, 150, and 210 days after birth (four development stages), and the mRNA and lncRNA expression was profiled. Compared to mRNA and lncRNA expression profiles, most differentially expressed mRNAs (DEmRNAs) were upregulated at 30 days; however, most DElncRNAs were downregulated at 210 days. Via Short Time-series Expression Miner (STEM) analysis and weighted gene co-expression network analysis (WGCNA), a complex interaction between mRNAs and lncRNAs was identified, indicating that lncRNAs may be a critical regulatory element for mRNAs. One module of genes in particular (module profile 4) was related to fibril organization, vasculogenesis, GTPase activator activity, and regulation of kinase activity. The mRNAs and lncRNAs in module profile 4 had a similar pattern of expression, indicating that they have functional and regulatory relationships. Only CAV1, PACSIN2, and CDC42 in the particular mRNA profile 4 were the target genes of lncRNAs in that profile, which shows the possible regulatory relationship between lncRNAs and mRNAs. The expression of these genes and lncRNAs in profile 4 was the highest at 30 days, and it is believed that these RNAs may play a critical role during the suckling period in order to meet the dietary requirements of piglets. In the lncRNA–mRNA co-expression network, the identified gene hubs and associated lncRNAs were shown to be involved in saccharide, lipid, and glucose metabolism, which may play an important role in the development and health of the liver. This result will lead to further investigation of liver lncRNA functions at various stages of development in Ningxiang pigs.

Recently, thousands of circular RNAs have been reported in different pig breeds. However, researches on the temporal and spatial expression patterns of circRNA over the period of animal growth are limited. Here, we systematically analyzed circRNAs in skeletal muscle and subcutaneous fat in four growth time points (30 days, 90 days, 150 days and 210 days after birth) of a Chinese native pig breed, Ningxiang pigs. A total of 1171 differentially expressed (DE) circRNAs between muscle and fat were identified, including 562 upregulated and 609 downregulated circRNAs. KEGG pathway enrichment analysis of these DE circRNAs revealed that host genes were mainly involved in glycogen metabolism signaling pathways, muscle development signaling pathways such as ErbB pathway and adipocytokine signaling pathways and AMPK signaling pathways and fatty acid biosynthesis. The circRNAs have striking spatiotemporal specificity in the form of dynamic expression at 90 d. Short Time-Series Expression Miner analysis showed multiple model spectra that are significantly enriched with time changes in muscle and fat. Our findings provide new ideas and perspectives about the role of circular RNAs and their targeting relations with mRNA and miRNA in skeletal muscle and fat tissue during pig growth.

Yaks are the main pillar of plateau animal husbandry and the material basis of local herdsmen’s survival. The level of mineral elements in the body is closely related to the production performance of yaks. In this study, we performed a comprehensive analysis of rumen epithelial morphology, transcriptomics and metagenomics to explore the dynamics of rumen functions, microbial colonization and functional interactions in yaks from birth to adulthood. Bacteria, eukaryotes, archaea and viruses colonized the rumen of yaks from birth to adulthood, with bacteria being the majority. Bacteroidetes and Firmicutes were the dominant phyla in five developmental stages, and the abundance of genus Lactobacillus and Fusobacterium significantly decreased with age. Glycoside hydrolase (GH) genes were the most highly represented in five different developmental stages, followed by glycosyltransferases (GTs) and carbohydrate-binding modules (CBMs), where the proportion of genes coding for CBMs increased with age. Integrating host transcriptome and microbial metagenome revealed 30 gene modules related to age, muscle layer thickness, nipple length and width of yaks. Among these, the MEmagenta and MEturquoise were positively correlated with these phenotypic traits. Twenty-two host genes involved in transcriptional regulation related to metal ion binding (including potassium, sodium, calcium, zinc, iron) were positively correlated with a rumen bacterial cluster 1 composed of Alloprevotella, Paludibacter, Arcobacter, Lactobacillus, Bilophila, etc. Therefore, these studies help us to understand the interaction between rumen host and microorganisms in yaks at different ages, and further provide a reliable theoretical basis for the development of feed and mineral element supplementation for yaks at different ages.