The development of hair follicle in cashmere goats shows significant periodic change, as with mice and humans. However, for cashmere goat with double-coat, the periodic change may be due to other regulatory molecules and signal pathways. To understand the mechanism of periodic development of hair follicle, we performed a weighted gene coexpression network analysis (WGCNA) to mine key genes and establish an interaction network by utilizing the NCBI public dataset. Ten coexpression modules, including 7689 protein-coding genes, were constructed by WGCNA, six of which are considered to be significantly related to the development of the hair follicle cycle. A functional enrichment analysis for each model showed that they are closely related to ECM- receptor interaction, focal adhesion, PI3K-Akt signaling pathway, estrogen signaling pathway, and so on. Combined with the analysis of differential expressed genes, 12 hub genes from coexpression modules were selected as candidate markers, i.e., COL1A1, C1QTNF6, COL1A2, AQP3, KRTAP3-1, KRTAP11-1, FA2H, NDUFS5, DERL2, MRPL14, ANTKMT and XAB2, which might be applied to improve cashmere production.

N⁶-methyladenosine (m⁶A) plays an important role in regulating gene expression. Previous studies found that m⁶A methylation affects skeletal muscle development. However, the effect of m⁶A methylases on bovine skeletal myogenesis is still unclear. Here, we found that the expression of m⁶A demethylases (FTO and ALKBH5) was significantly higher in the longissimus dorsi muscle of adult cattle than in newborn cattle. In contrast, the expression of m⁶A methyltransferases (METTL3, METTL14 and WTAP) was reduced. The mRNA expression of all five genes was found to be increased during the myogenesis of myoblasts in vitro. Knockdown of FTO or METTL3 promoted myoblast proliferation, inhibited myoblast apoptosis and suppressed myogenic differentiation, whereas ALKBH5 knockdown had the opposite effect. METTL14 knockdown enhanced myoblast proliferation and impaired myogenic differentiation. WTAP knockdown attenuated proliferation and contributed to apoptosis but did not affect differentiation. Furthermore, the functional domains of these five m⁶A methylases are conserved across species. Our results suggest that m⁶A methylases are involved in regulating skeletal muscle development and that there may be a complex network of m⁶A methylation regulating skeletal myogenesis.

The intramuscular fat (or marbling fat) content is an essential economic trait of beef cattle and improves the flavor and palatability of meat. Several studies have highlighted the correlation between long non-coding RNAs (lncRNAs) and intramuscular fat development; however, the precise molecular mechanism remains unknown. Previously, through a high-throughput sequencing analysis, we found a lncRNA and named it a long non-coding RNA BNIP3 (lncBNIP3). The 5′ RACE and 3′ RACE explored 1945 bp total length of lncBNIP3, including 1621 bp of 5′RACE, and 464 bp of 3′RACE. The nucleoplasmic separation and FISH results explored the nuclear localization of lncBNIP3. Moreover, the tissue expression of lncBNIP3 was higher in the longissimus dorsi muscle, followed by intramuscular fat. Furthermore, down-regulation of lncBNIP3 increased the 5-Ethynyl-2′- deoxyuridine (EdU)-EdU-positive cells. The flow cytometry results showed that the number of cells in the S phase was significantly higher in preadipocytes transfected with si-lncBNIP3 than in the control group (si-NC). Similarly, CCK8 results showed that the number of cells after transfection of si-lncBNIP3 was significantly higher than in the control group. In addition, the mRNA expressions of proliferative marker genes CyclinB1 (CCNB1) and Proliferating Cell Nuclear Antigen (PCNA) in the si-lncBNIP3 group were significantly higher than in the control group. The Western Blot (WB) results also showed that the protein expression level of PCNA transfection of si-lncBNIP3 was significantly higher than in the control group. Similarly, the enrichment of lncBNIP3 significantly decreased the EdU-positive cells in the bovine preadipocytes. The results of flow cytometry and CCK8 assay also showed that overexpression of lncBNIP3 inhibited the proliferation of bovine preadipocytes. In addition, the overexpression of lncBNIP3 significantly inhibited the mRNA expressions of CCNB1 and PCNA. The WB results showed that the overexpression of lncBNIP3 significantly inhibited the expression of the CCNB1 protein level. To further explore the mechanism of lncBNIP3 on the proliferation of intramuscular preadipocytes, RNA-seq was performed after interference with si-lncBNIP3, and 660 differentially expressed genes (DEGs) were found, including 417 up-regulated DEGs and 243 down-regulated DEGs. The KEGG pathway analysis showed that the cell cycle was the most significant pathway for the functional enrichment of DEGs, followed by the DNA replication pathway. The RT-qPCR quantified the expression of twenty DEGs in the cell cycle. Therefore, we speculated that lncBNIP3 regulated intramuscular preadipocyte proliferation through the cell cycle and DNA replication pathways. To further confirm this hypothesis, the cell cycle inhibitor Ara-C was used to inhibit DNA replication of the S phase in intramuscular preadipocytes. Herein, Ara-C and si-lncBNIP3 were simultaneously added to the preadipocytes, and the CCK8, flow cytometry, and EdU assays were performed. The results showed that the si-lncBNIP3 could rescue the inhibitory effect of Ara-C in the bovine preadipocyte proliferation. In addition, lncBNIP3 could bind to the promoter of cell division control protein 6 (CDC6), and down-regulation of lncBNIP3 promoted the transcription activity and the expression of CDC6. Therefore, the inhibitory effect of lncBNIP3 on cell proliferation might be understood through the cell cycle pathway and CDC6 expression. This study provided a valuable lncRNA with functional roles in intramuscular fat accumulation and revealed new strategies for improving beef quality.

Intramuscular fat (IMF) is a quality index associated with the taste and juiciness of meat. The deposition of IMF is affected by genetic and non-genetic factors, such as age, slaughter location, gender of the animal, and diet. Micro-ribonucleic acids (miRNA) are transcriptional regulators involved in adipogenesis, but the specific role of miR-376a in regulation of bovine adipocytes remains unknown. Our findings indicated that miR-376a was a potential negative regulator of bovine adipocyte differentiation. A bta-miR-376a mimic inhibited mRNA and protein expression of the marker genes, CDK1, CDK2, PCNA, C/EBPα, FAS, and PPAR γ, and significantly reduced ratios (%) of S-phase cells, the number of cells stained with 5-ethynyl-2′-deoxyuridine, and adipocyte proliferation. Oil red O staining and triglyceride content analysis also confirmed that bta-miR-376a was involved in adipocyte differentiation. Luciferase activities confirmed that Krüppel-like transcription factor 15 (KLF15) was a direct target gene of bta-miR-376a, and that KLF15 was a key transcription factor in adipogenesis. Therefore, bta-miR-376a might be a target for increasing beef IMF.

N⁶-methyladenosine (m⁶A) is the most prevalent methylation modification of eukaryotic mRNA, and it plays an important role in regulating gene expression. Previous studies have found that m6A methylation plays a role in mammalian skeletal muscle development. However, the effect of m⁶A on bovine skeletal myogenesis are still unclear. Here, we selected proliferating myoblasts (GM) and differentiated myotubes (on the 4th day of differentiation, DM) for m⁶A-seq and RNA-seq to explore the m⁶A methylation modification pattern during bovine skeletal myogenesis. m⁶A-seq analysis revealed that m⁶A methylation was an abundant modification of the mRNA in bovine myoblasts and myotubes. We scanned 5,691–8,094 m⁶A-modified transcripts, including 1,437 differentially methylated genes (DMGs). GO and KEGG analyses revealed that DMGs were primarily involved in transcriptional regulation and RNA metabolism, as well as insulin resistance and metabolic pathways related to muscle development. The combined analysis further identified 268 genes that had significant changes at both m⁶A and mRNA levels, suggesting that m⁶A modification may regulate myoblast differentiation by mediating the expression of these genes. Furthermore, we experimentally confirmed four genes related to myogenesis, including MYOZ2, TWIST1, KLF5 and MYOD1, with differential changes in both m⁶A and mRNA levels during bovine myoblast differentiation, indicating that they can be potential candidate targets for m⁶A regulation of skeletal myogenesis. Our results may provide new insight into molecular genetics and breeding of beef cattle, and provide a reference for investigating the mechanism of m⁶A regulating skeletal muscle development.