The development of hair follicles in yak shows significant seasonal cycles. In our previous research, transcriptome data including mRNAs and lncRNAs in five stages during the yak hair follicles (HFs) cycle were detected, but their regulation network and the hub genes in different periods are yet to be explored. This study aimed to screen and identify the hub genes during yak HFs cycle by constructing a mRNA-lncRNA co-expression network. A total of 5000 differently expressed mRNA (DEMs) and 729 differently expressed long noncoding RNA (DELs) were used to construct the co-expression network, based on weighted genes co-expression network analysis (WGCNA). Four temporally specific modules were considered to be significantly associated with the HFs cycle of yak. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the modules are enriched into Wnt, EMC-receptor interaction, PI3K-Akt, focal adhesion pathways, and so on. The hub genes, such as FER, ELMO1, PCOLCE, and HOXC13, were screened in different modules. Five hub genes (WNT5A, HOXC13, DLX3, FOXN1, and OVOL1) and part of key lncRNAs were identified for specific expression in skin tissue. Furthermore, immunofluorescence staining and Western blotting results showed that the expression location and abundance of DLX3 and OVOL1 are changed following the process of the HFs cycle, which further demonstrated that these two hub genes may play important roles in HFs development.

Investigating the critical genes related to milk synthesis is essential for the improvement of the milk yield of the yak. Real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable and widely used method to measure and evaluate gene expression levels. Selection of suitable reference genes is mandatory to acquire accurate normalization of gene expression results from RT-qPCR. To select the most stable reference genes for reliable normalization of mRNA expression by RT-qPCR in the mammary gland of the Ashidan yak, we selected 16 candidate reference genes and analyzed their expression stability at different physiological stages (lactation and dry period). The expression stability of the candidate reference genes was assessed using geNorm, NormFinder, BestKeeper, Delta Ct, and RefFinder methods. The results showed that the hydroxymethylbilane synthase gene (HMBS) and the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide gene (YWHAZ) were the most stable genes across all treatment samples. The reliability of selected reference genes was validated by normalizing relative expression of the lactation-related 60S ribosomal protein L35 gene (RPL35). The relative expression of RPL35 varied considerably according to the different reference genes. This work provides valuable information to further promote research in the molecular mechanisms involved in lactation and mammary gland development and provides a foundation for the improvement of the milk yield and quality of the Ashidan yak.

Single-cell sequencing technology can fully reflect the heterogeneity of cell populations at the single cell level, making it possible for us to re-recognize various tissues and organs. At present, the sequencing study of hair follicles is transiting from the traditional ordinary transcriptome level to the single cell level, which will provide diverse insights into the function of hair follicle cells. This review focuses on research advances in the hair follicle microenvironment obtained from scRNA-seq studies of major cell types in hair follicle development, with a special emphasis on the discovery of new subpopulations of hair follicles by single-cell techniques. We also discuss the problems and current solutions in scRNA-seq observation and look forward to its prospects.

Testis has an important function in male reproduction. Its development is regulated by a large number of genes. The real-time reserve transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a useful tool to evaluate the gene expression levels. However, unsuitable reference genes (RGs) can cause the misinterpretation of gene expression levels. Unfortunately, the ideal RGs for yak testis development are yet to be studied. In this study, 13 commonly used RGs were selected to identify the most stable RGs in yak testis at four different developmental stages, including two immature stages (6 months and 18 months) and two mature stages (30 months and 6 years). This study used GeNorm, NormFinder, BestKeeper, ∆Ct, and RefFinder programs to evaluate the stability of 13 candidate genes. The results of RefFinder showed that the stabilities of TATA box-binding protein (TBP) and ubiquitously expressed transcript protein (UXT) were ranked the top two across all developmental stages. TBP and hydroxymethylbilane synthase (HMBS) were stably expressed in immature stages, while mitochondrial ribosomal protein L39 (MRPL39) and TBP had higher stability than other candidate genes in mature stages. This study provided valuable information for gene expression studies to assist further investigation on the molecular mechanisms in underlying yak testis development.