Circular RNAs (circRNAs) have been identified from various tissues and species, but their regulatory functions during developmental processes are not well understood. We examined circRNA expression profiles of two developmental stages of bovine skeletal muscle (embryonic and adult musculus longissimus) to provide first insights into their potential involvement in bovine myogenesis. We identified 12 981 circRNAs and annotated them to the Bos taurus reference genome, including 530 circular intronic RNAs (ciRNAs). One parental gene could generate multiple circRNA isoforms, with only one or two isoforms being expressed at higher expression levels. Also, several host genes produced different isoforms when comparing development stages. Most circRNA candidates contained two to seven exons, and genomic distances to back-splicing sites were usually less than 50 kb. The length of upstream or downstream flanking introns was usually less than 105 nt (mean≈11 000 nt). Several circRNAs differed in abundance between developmental stages, and real-time quantitative PCR (qPCR) analysis largely confirmed differential expression of the 17 circRNAs included in this analysis. The second part of our study characterized the role of circLMO7—one of the most down-regulated circRNAs when comparing adult to embryonic muscle tissue—in bovine muscle development. Overexpression of circLMO7 inhibited the differentiation of primary bovine myoblasts, and it appears to function as a competing endogenous RNA for miR-378a-3p, whose involvement in bovine muscle development has been characterized beforehand. Congruent with our interpretation, circLMO7 increased the number of myoblasts in the S-phase of the cell cycle and decreased the proportion of cells in the G0/G1 phase. Moreover, it promoted the proliferation of myoblasts and protected them from apoptosis. Our study provides novel insights into the regulatory mechanisms underlying skeletal muscle development and identifies a number of circRNAs whose regulatory potential will need to be explored in the future.

MicroRNAs (miRNAs), a class of single stranded, small (~22 nucleotides), non-coding RNAs, play an important role in muscle development. We focused on the role of the miR-30-5p family during bovine muscle development from previous high-throughput sequencing results and analyzed their expression profiles. MHC and MyoG mRNAs expression as well as their proteins were suppressed in differentiated C2C12 cells, suggesting the importance of miR-30-5p in muscle development. MBNL, the candidate target of miR-30-5p, is an alternative splicing regulation factor. MBNL1 and MBNL3 have opposite effects on muscle differentiation. Our results confirmed that miR-30a-5p and miR-30e-5p repress the expression of MBNL1, MBNL2 and MBNL3, whereas miR-30b-5p inhibits MBNL1 and MBNL2 expression. This provides direct evidence that MBNL expression can be flexibly regulated by miR-30-5p. Previous studies showed that MBNL1 promotes exon inclusion of two muscle-related genes (Trim55 and INSR). Through RNA splicing studies, we found that miR-30-5p had an effect on their alternative splicing, which means miR-30-5p via MBNL1 could be integrated into muscle signaling pathways in which INSR or Trim55 are located. In conclusion, miR-30-5p could inhibit muscle cell differentiation and regulate the alternative splicing of Trim55 and INSR by targeting MBNL. These results promote the understanding of the function of miRNAs in muscle development.

Circular RNAs (circRNAs) have been identified as a novel type of regulators involved in multiple biological processes. However, circRNAs with a potential function in fat deposition in buffalo are poorly understood. In this study, six RNA libraries of adipose tissue were constructed for three young and three adult Chinese buffaloes with paired-ends RNA sequencing using the Illumina HiSeq 3000 platform. A total of 5141 circRNAs were computationally identified. Among them, 252 circRNAs were differentially expressed (DE) between the young and adult buffaloes. Of these, 54 were upregulated and 198 were downregulated in the adult group. Eight DE circRNAs were further identified by quantitative real-time-PCR (qRT-PCR) and Sanger sequencing. Co-expression analysis revealed that 34 circRNAs demonstrated a strong correlation with fat deposition-associated genes (|r| > 0.980). Among these, expressional correlation between two circRNAs (19:45387150|45389986 and 21:6969877|69753491) and PR/SET domain 16 was further verified using qRT-PCR, and a strong correlation was revealed (1 > |r| > 0.8). These results strongly suggest that circRNAs 19:45387150|45389986 and 21:6969877|69753491 are potential regulators of buffalo fat deposition. In summary, this study characterized the circRNA profiles of adipose tissues at different stages for the first time and revealed two circRNAs strongly correlated with fat deposition-associated genes, which provided new candidate regulators for fat deposition in buffalo.

The buffalo (Bubalus bubalis L.) is prevalent in China and the increasing demand for meat production has changed its role from being a beast of burden to a meat source. The low fat deposition level has become one of the main barriers for its use in meat production. It is urgent to reveal factors involved in fat deposition in buffalo. This study performed RNA sequencing to investigate both long noncoding RNAs (lncRNAs) and mRNAs of adipose tissues in young and adult buffalos. A total of 124 lncRNAs and 2008 mRNAs showed differential expression patterns between young and adult samples. Coexpression analysis and functional enrichment revealed 585 mRNA–lncRNA pairs with potential function in fat deposition. After validation by qRT-PCR, we focused on a lncRNA transcribed from the ubiquinone oxidoreductase subunit C2 (NDUFC2) antisense (AS) strand which showed high correlation with thyroid hormone responsive protein (THRSP). NDUFC2-AS lncRNA is highly expressed in adipose tissue and maturation adipocytes and mainly exists in the nucleus. Functional assays demonstrated that NDUFC2-AS lncRNA promotes adipogenic differentiation by upregulating the expression levels of THRSP and CCAAT enhancer binding protein alpha (C/EBPα) in buffalo. These results indicate that NDUFC2-AS lncRNA promotes fat deposition in buffalo.

Adipogenesis is a complex and precisely orchestrated process mediated by a series of adipogenic regulatory factors. Recent studies have highlighted the importance of microRNAs (miRNAs) in diverse biological processes, most specifically in regulating cell differentiation and proliferation. However, the mechanisms of miRNAs in adipogenesis are largely unknown. In this study, we found that miR-107 expression was higher in bovine adipose tissue than that in other tissues, and there was a downregulation trend during adipocyte differentiation. To explore the function of miR-107 in adipocyte differentiation, agomiR-107 and antiagomiR-107 were transfected into bovine adipocytes, respectively. Oil Red O staining, CCK-8, EdU assays, RT-qPCR, and Western blotting were performed, and the results showed that overexpressed miR-107 significantly suppressed fat deposition and adipocyte differentiation, while knockdown of miR-107 promoted fat deposition and adipocytes differentiation. In addition, through bioinformatics analysis, luciferase reporter assays, RT-qPCR, and Western blotting, we identified apolipoprotein 2 (APOC2) as a target of miR-107. Transfection of siRNA-APOC2 into adipocytes led to suppression in adipocyte differentiation and proliferation, suggesting a positive role of APOC2 in bovine lipogenesis. In summary, our findings suggested that miR-107 regulates bovine adipocyte differentiation and lipogenesis by directly targeting APOC2, and these results. These theoretical and experimental basis for future clarification of the regulation mechanism of adipocyte differentiation and lipogenesis. Moreover, for the highly conserved among different species, miR-107 may be a potential molecular target to be used for the treatment of lipid-related diseases in the future.