Preadipocyte differentiation plays a critical role in subcutaneous fat deposition in pigs. However, the roles of different RNAs, such as messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) in the differentiation process of subcutaneous preadipocytes, are still largely unclear. In the present study, a transcriptome analysis, including the analysis of mRNAs, lncRNAs, and circRNAs, during different differentiation stages, namely, day 0 (D0), day 2 (D2), day 4 (D4), and day 8 (D8), of subcutaneous preadipocytes from Chinese Erhualian pigs was performed. A total of 1554, 470, 1344, 1777, and 676 differentially expressed (DE) mRNAs, 112, 58, 95, 136, and 93 DE lncRNAs, and 902, 787, 710, 932, and 850 DE circRNAs were identified between D2 and D0, between D4 and D2, between D8 and D4, between D4 and D0, and between D8 and D0, respectively. Furthermore, functional enrichment analysis showed that the common DE mRNAs during the entire differentiation process were mainly involved in lipid metabolic and cell differentiation processes. Additionally, co-expression network analysis identified the potential lncRNAs related to adipogenesis, e.g., MSTRG.131380 and MSTRG.62128. Our study provides new insights of the expression changes of RNAs during adipogenic differentiation, which might contribute to the phenotype of subcutaneous adipogenesis. These results greatly improve our understanding of the molecular mechanisms regulating subcutaneous fat deposition in pigs.

Kaiqing Liu, Wensai Yu, Wei Wei, Xinbao Zhang, Ye Tian, Melak Sherif, Xin Liu, Chao Dong, Wangjun Wu, Lifan Zhang, Jie Chen

Laying performance is an important economical trait of goose production. As laying performance is of low heritability, it is of significance to develop a marker-assisted selection (MAS) strategy for this trait. Definition of sequence variation related to the target trait is a prerequisite of quantitating MAS, but little is presently known about the goose genome, which greatly hinders the identification of genetic markers for the laying traits of geese. Recently developed restriction site-associated DNA (RAD) sequencing is a possible approach for discerning large-scale single nucleotide polymorphism (SNP) and reducing the complexity of a genome without having reference genomic information available. In the present study, we developed a pooled RAD sequencing strategy for detecting geese laying-related SNP. Two DNA pools were constructed, each consisting of equal amounts of genomic DNA from 10 individuals with either high estimated breeding value (HEBV) or low estimated breeding value (LEBV). A total of 139,013 SNP were obtained from 42,291,356 sequences, of which 18,771,943 were for LEBV and 23,519,413 were for HEBV cohorts. Fifty-five SNP which had different allelic frequencies in the two DNA pools were further validated by individual-based AS-PCR genotyping in the LEBV and HEBV cohorts. Ten out of 55 SNP exhibited distinct allele distributions in these two cohorts. These 10 SNP were further genotyped in a goose population of 492 geese to verify the association with egg numbers. The result showed that 8 of 10 SNP were associated with egg numbers. Additionally, liner regression analysis revealed that SNP Record-111407, 106975 and 112359 were involved in a multiplegene network affecting laying performance. We used IPCR to extend the unknown regions flanking the candidate RAD tags. The obtained sequences were subjected to BLAST to retrieve the orthologous genes in either ducks or chickens. Five novel genes were cloned for geese which harbored the candidate laying-related SNP, including membrane associated guanylate kinase (MAGI-1), KIAA1462, Rho GTPase activating protein 21 (ARHGAP21), acyl-CoA synthetase family member 2 (ACSF2), astrotactin 2 (ASTN2). Collectively, our data suggests that 8 SNP and 5 genes might be promising candidate markers or targets for marker-assisted selection of egg numbers in geese.