MicroRNAs (miRNAs) are endogenous small noncoding RNA molecules that negatively regulate gene expression. Herein, we investigated a selective number of miRNAs for their expression in skin tissue of Liaoning Cashmere goat during hair follicle cycles, and their intracellular regulatory networks were constructed based on bioinformatics analysis. The relative expression of six miRNAs (mir-103-3p, -15b-5p, 17-5p, -200b, -25-3p, and -30c-5p) at anagen phase is significantly higher than that at catagen and/or telogen phases. In comparison to anagen, the relative expression of seven miRNAs (mir-148a-3p, -199a-3p, -199a-5p, -24-3p, -30a-5p, -30e-5p, and -29a-3p) was revealed to be significantly up-regulated at catagen and/or telogen stages. The network analyses of miRNAs indicated those miRNAs investigated might be directly or indirectly involved in several signaling pathways through their target genes. These results provided a foundation for further insight into the roles of these miRNAs in skin tissue of Liaoning Cashmere goat during hair follicle cycles.

Long noncoding RNAs (lncRNAs) are a novel class of eukaryotic transcripts. They are thought to act as a critical regulator of protein-coding gene expression. Herein, we identified and characterized 13 putative lncRNAs from the expressed sequence tags from secondary hair follicle of Cashmere goat. Furthermore, we investigated their transcriptional pattern in secondary hair follicle of Liaoning Cashmere goat during telogen and anagen phases. Also, we generated intracellular regulatory networks of upregulated lncRNAs at anagen in Wnt signaling pathway based on bioinformatics analysis. The relative expression of six putative lncRNAs (lncRNA-599618, -599556, -599554, -599547, -599531, and -599509) at the anagen phase is significantly higher than that at telogen. Compared with anagen, the relative expression of four putative lncRNAs (lncRNA-599528, -599518, -599511, and -599497) was found to be significantly upregulated at telogen phase. The network generated showed that a rich and complex regulatory relationship of the putative lncRNAs and related miRNAs with their target genes in Wnt signaling pathway. Our results from the present study provided a foundation for further elucidating the functional and regulatory mechanisms of these putative lncRNAs in the development of secondary hair follicle and cashmere fiber growth of Cashmere goat.

Quantitative real-time PCR is the most sensitive technique for gene expression analysis. Data normalization is essential to correct for potential errors incurred in all steps from RNA isolation to PCR amplification. The commonly accepted approach for normalization is the use of reference gene. Until now, no suitable reference genes have been available for data normalization of gene expression in milk somatic cells of lactating yaks across lactation. In the present study, we evaluated the transcriptional stability of 10 candidate reference genes in milk somatic cells of lactating yak, including ACTB, B2M, GAPDH, GTP, MRPL39, PPP1R11, RPS9, RPS15, UXT, and RN18S1. Four genes, RPS9, PPP1R11, UXT, and MRPL39, were identified as being the most stable genes in milk somatic cells of lactating yak. Using the combination of RPS9, PPP1R11, UXT, and MRPL39 as reference genes, we further assessed the relative expression of 4 genes of interest in milk somatic cells of yak across lactation, including ELF5, ABCG2, SREBF2, and DGAT1. Compared with expression in colostrum, the overall transcription levels of ELF5, ABCG2, and SREBF2 in milk were found to be significantly upregulated in early, peak, and late lactation, and significantly downregulated thereafter, before the dry period. A similar pattern was observed in the relative expression of DGAT1, but no significant difference was revealed in its expression in milk from late lactation compared with colostrum. Based on these results, we suggest that the geometric mean of RPS9, PPP1R11, UXT, and MRPL39 can be used for normalization of real-time PCR data in milk somatic cells of lactating yak, if similar experiments are performed.

Circular RNAs (CircRNAs) are a type of non-coding RNAs, which contain a covalently closed loop structure without 5′ to 3′ free ends. CircRNAs play essential roles in the regeneration of secondary hair follicle (SHF) and cashmere growth in goats. CircRNA-1926 was previously identified in SHF of cashmere goats, but its potential roles are unclear. In this study, we confirmed the expression of circRNA-1926 in SHF bulge of nine cashmere goats with a significantly higher level at anagen than that of telogen. Through the use of both overexpression and siRNA interference, we showed that circRNA-1926 promoted the differentiation of SHF stem cell into hair follicle lineage in cashmere goats which was evaluated via indictor genes Keratin 7 and Keratin 17. Using RNA pull-down, we found that circRNA-1926 bound with miR-148a/b-3p. Additionally, our data indicated that circRNA-1926 promoted the expression of the CDK19 gene. Using dual-luciferase reporter assays, it was revealed that circRNA-1926 positively regulated the CDK19 expression through miR-148a/b-3p. The results from this study demonstrated that circRNA-1926 contributes the differentiation of SHF stem cells into hair follicle lineages in cashmere goats via sponging miR-148a/b-3p to enhance CDK19 expression.