A genome-wide landscape of mRNAs, lncRNAs, circRNAs and miRNAs during intramuscular adipogenesis in cattle.

BACKGROUND: Intramuscular preadipocyte differentiation plays a critical role in bovine intramuscular fat (IMF) deposition. However, the roles of different RNAs, including mRNAs, circRNAs, lncRNAs and miRNAs, in regulating the adipogenic differentiation of intramuscular preadipocytes remain largely unclear. RESULTS: In the present study, a whole transcriptome sequencing and analysis, including the analysis of mRNAs, circRNAs, lncRNAs and miRNAs, during different differentiation stages (0, 3, 6, and 9 d) of intramuscular preadipocytes from Qinchuan cattle was performed. All samples were prepared with 3 biological replicates. Here, a total of 27,153 mRNAs, 14,070 circRNAs, 7035 lncRNAs, and 427 miRNAs were annotated. Among them, we identified 4848 differentially expressed mRNAs (DEMs), 181 DE circRNAs (DECs), 501 DE lncRNAs (DELs) and 77 DE miRNAs (DEmiRs) between 0 d and other differentiation days (3, 6, and 9 d). GO and KEGG functional enrichment analyses showed that these differentially expressed genes were mainly enriched in cell differentiation, fat metabolism and adipogenesis-related pathways. Furthermore, weighted gene coexpression network analysis (WGCNA) and co-expression network analysis screened out multiple important mRNAs, circRNAs and lncRNAs related to intramuscular adipogenesis. Based on the competing endogenous RNA (ceRNA) regulatory mechanism, we finally identified 24 potential ceRNA networks and 31 potential key genes, including FOXO1/miR-330/circRNA2018/MSTRG.20301, GPAM/miR-27b/ciRNA489 and SESN3/miR-433/circRNA2627MSTRG.20342. CONCLUSIONS: This study provides new insights into the differential expression patterns of different transcript types (i.e., mRNAs, circRNAs, lncRNAs and miRNAs) in intramuscular preadipocyte differentiation. Our findings provide data support for studying the molecular mechanism of key mRNAs and noncoding RNAs in IMF deposition, and provide new candidate markers for the molecular breeding of beef cattle.

SERPINE1AS2 regulates intramuscular adipogenesis by inhibiting PAI1 protein expression.

Antisense long non-coding RNAs (lncRNAs) played a crucial role in the precise regulation of essential biological processes and were abundantly present in animals. Many of these antisense lncRNAs have been identified as key roles in adipose tissue accumulation in livestock, underscoring their vital role in the regulation of animal physiology. Nonetheless, the functional roles of these antisense lncRNAs in regulating adipogenesis and the specific molecular mechanisms these processes were still unclear, which was a significant gap in current scientific research. In this study, we identified and characterized SERPINE1AS2, a novel natural antisense lncRNA, was highly expressed in the fat tissues of adult cattle and calves. Its expression gradually increased during the differentiation of intramuscular adipocytes. Through functional studies, we observed that knockdown of SERPINE1AS2 inhibited the proliferation and adipogenesis of intramuscular adipocytes, while overexpression of SERPINE1AS2 produced the opposite effect. RNA sequencing (RNA-seq) analysis following SERPINE1AS2 knockdown revealed that differential expression genes (DEGs) were significantly enriched in key signaling pathways, notably the MAPK, Wnt, and mTOR signaling pathways. Furthermore, SERPINE1AS2 interacted with Plasminogen Activator Inhibitor-1 (PAI1), forming RNA dimers through complementary base pairing and consequently influencing PAI1 expression. Interestingly, studies on PAI1 suggested that reduced expression facilitated adipogenesis and the downregulation of PAI1 alleviated the inhibitory effect of reduced SERPINE1AS2 on adipogenesis. In summary, this study suggested that SERPINE1AS2 played a crucial role in the adipogenesis of bovine intramuscular adipocytes by modulating the expression of PAI1. SERPINE1AS2 also regulated adipogenesis by engaging in the MAPK, Wnt, and mTOR signaling pathways. Our results suggested that SERPINE1AS2 had a complex regulatory mechanism on adipogenesis in intramuscular adipocytes.

Bta-miR-330 promotes bovine intramuscular pre-adipocytes adipogenesis via targeting SESN3 to activate the Akt-mTOR signaling pathway.

Consumers are more inclined to choose beef with a high intramuscular fat content (IMF), which regulated by lots of factors. It is very significant to find a miRNA that plays a key role in the accumulation of IMF. In our study, we found that bta-miR-330 was highly expressed in Japanese black cattle and differentially expressed at intramuscular pre-adipocytes differentiation processes. Furthermore, we transfected the bta-miR-330 mimic & inhibitor in intramuscular pre-adipocytes. The results showed that bta-miR-330 inhibits the proliferation but promotes the adipogenesis of intramuscular pre-adipocytes. Subsequently, our study showed that bta-miR-330 binds to SESN3, which inhibits the adipogenesis of intramuscular pre-adipocytes. Moreover, we established the mechanism that bta-miR-330 promotes the adipogenesis of intramuscular pre-adipocytes by targeting SESN3 to activate the Akt-mTOR signaling pathway. Overall, our results revealed that bta-miR-330-SESN3-Akt-mTOR axis plays an important role in adipogenesis of intramuscular pre-adipocytes, which provides a molecular basis for increasing IMF content in beef cattle.

RNA-seq analysis reveals the critical role of the novel lncRNA BIANCR in intramuscular adipogenesis through the ERK1/2 signaling pathway.

BACKGROUND: Long non-coding RNAs (lncRNAs) regulate numerous biological processes, including adipogenesis. Research on adipogenesis will assist in the treatment of human metabolic diseases and improve meat quality in livestock, such as the content of intramuscular fat (IMF). However, the significance of lncRNAs in intramuscular adipogenesis remains unclear. This research aimed to reveal the lncRNAs transcriptomic profiles in the process of bovine intramuscular adipogenesis and to identify the lncRNAs involved in the adipogenesis of bovine intramuscular adipocytes. RESULTS: In this research, a landscape of lncRNAs was identified with RNA-seq in bovine intramuscular adipocytes at four adipogenesis stages (0 d, 3 d, 6 d, and 9 d after differentiation). A total of 7035 lncRNAs were detected, including 3396 novel lncRNAs. Based on the results of differential analysis, co-expression analysis, and functional prediction, we focused on the bovine intramuscular adipogenesis-associated long non-coding RNA (BIANCR), a novel lncRNA that may have an important regulatory function. The knockdown of BIANCR inhibited proliferation and promoted apoptosis of intramuscular preadipocytes. Moreover, BIANCR knockdown inhibited intramuscular adipogenesis by regulating the ERK1/2 signaling pathway. CONCLUSION: This study obtained the landscape of lncRNAs during adipogenesis in bovine intramuscular adipocytes. BIANCR plays a crucial role in adipogenesis through the ERK1/2 signaling pathway. The results are noteworthy for improving beef meat quality, molecular breeding, and metabolic disease research.

Interactive regulation of DNA demethylase gene TET1 and m6A methyltransferase gene METTL3 in myoblast differentiation.

DNA methylation (5mC) and mRNA N6-methyladenosine (m6A) play an essential role in gene transcriptional regulation. DNA methylation has been well established to be involved in skeletal muscle development. Interacting regulatory mechanisms between DNA methylation and mRNA m6A modification have been identified in a variety of biological processes. However, the effect of m6A on skeletal muscle differentiation and the underlying mechanisms are still unclear. It is also unknown whether there is an interaction between DNA methylation and mRNA m6A modification in skeletal myogenesis. In the present study, we used m6A-IP-qPCR, LC-MS/MS and dot blot assays to determine that the DNA demethylase gene, TET1, exhibited increased m6A levels and decreased mRNA expression during bovine skeletal myoblast differentiation. Dual-luciferase reporter assays and RIP experiments demonstrated that METTL3 suppressed TET1 expression by regulating TET1 mRNA stability in a m6A-YTHDF2-dependent manner. Furthermore, TET1 mediated DNA demethylation of itself, MYOD1 and MYOG, thereby stimulating their expression to promote myogenic differentiation. Ectopic expression of TET1 rescued the effect of METTL3 knockdown on reduced myotubes. In contrast, TET1 knockdown impaired the myogenic differentiation promoted by METTL3 overexpression. Moreover, ChIP experiments found that TET1 could bind and demethylate METTL3 DNA, which enhanced METTL3 expression. In addition, TET1 knockdown decreased m6A levels. ChIP assays also showed that TET1 knockdown contributed to the binding of H3K4me3 and H3K27me3 to METTL3 DNA. Our results revealed a negative feedback regulatory loop between TET1 and METTL3 in myoblast differentiation, which unveiled the interplay among DNA methylation, RNA methylation and histone methylation in skeletal myogenesis.

m6A Methylases Regulate Myoblast Proliferation, Apoptosis and Differentiation.

N6-methyladenosine (m6A) plays an important role in regulating gene expression. Previous studies found that m6A methylation affects skeletal muscle development. However, the effect of m6A methylases on bovine skeletal myogenesis is still unclear. Here, we found that the expression of m6A demethylases (FTO and ALKBH5) was significantly higher in the longissimus dorsi muscle of adult cattle than in newborn cattle. In contrast, the expression of m6A 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 m6A methylases are conserved across species. Our results suggest that m6A methylases are involved in regulating skeletal muscle development and that there may be a complex network of m6A methylation regulating skeletal myogenesis.

Proliferation of bovine myoblast by LncPRRX1 via regulation of the miR-137/CDC42 axis.

Noncoding RNAs, such as long noncoding RNAs (lncRNAs), are abundant in livestock. Many lncRNAs that affect the growth rate of livestock have been identified in muscles. However, some of their physiological functions and regulatory mechanisms remain unclear. In this study, we identified a new lncRNA (lncPRRX1) and investigated its effect on the proliferation of bovine myoblasts. LncPRRX1 was highly expressed in muscle tissue, and interference with lncPRRX1 inhibited the proliferation of bovine myoblasts in vitro. The RNA molecules of lncPRRX1 act on miR-137 as competitive endogenous RNAs (ceRNAs). Overexpression of miR-137 suppressed the proliferation of myoblasts, while inhibition of miR-137 had the opposite effect. In addition, the predicted target genes of miR-137 were significantly enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, in which Cell Division Cycle 42 (CDC42) was shown to be the direct target gene of miR-137, and interference with CDC42 inhibited myoblast proliferation. Furthermore, interference with lncPRRX1 repaired the defects in CDC42 protein levels and cell proliferation caused by miR-137 inhibitors. Our results suggested that lncPRRX1 promoted bovine myoblast proliferation by regulating the miRNA-137/CDC42 axis.

Transcriptome-wide N 6-Methyladenosine Methylome Profiling Reveals m6A Regulation of Skeletal Myoblast Differentiation in Cattle (Bos taurus).

N 6 -methyladenosine (m6A) 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 m6A on bovine skeletal myogenesis are still unclear. Here, we selected proliferating myoblasts (GM) and differentiated myotubes (on the 4th day of differentiation, DM) for m6A-seq and RNA-seq to explore the m6A methylation modification pattern during bovine skeletal myogenesis. m6A-seq analysis revealed that m6A methylation was an abundant modification of the mRNA in bovine myoblasts and myotubes. We scanned 5,691-8,094 m6A-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 m6A and mRNA levels, suggesting that m6A 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 m6A and mRNA levels during bovine myoblast differentiation, indicating that they can be potential candidate targets for m6A 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 m6A regulating skeletal muscle development.

Corrigendum: Bovine Pre-adipocyte Adipogenesis Is Regulated by bta-miR-150 Through mTOR Signaling.

[This corrects the article DOI: 10.3389/fgene.2021.636550.].

Bovine Pre-adipocyte Adipogenesis Is Regulated by bta-miR-150 Through mTOR Signaling.

Micro RNA (miR) are recognized for their important roles in biological processes, particularly in regulatory componentization. Among the miR, miR-150 has been the focus of intense scrutiny, mostly due to its role in malignant tumors. A comparison between steer and bull adipose tissues identified bta-miR-150 as one of the nine downregulated miRNAs, although its function remains unknown (GEO:GSE75063). The present study aimed to further characterize the role of bta-miR-150 in cattle. bta-miR-150 has a negative regulatory effect on the differentiation of bovine adipocytes and promotes proliferation. Overexpression of bta-miR-150 can promote mRNA and protein expression of the marker genes CDK1, CDK2, and PCNA, increase the number of EdU-stained cells, promote adipocyte proliferation, inhibit adipocyte differentiation, and reduce lipid droplet formation. Results of RNA-seq and WGCNA analyses showed that the mammalian target of the rapamycin signaling pathway, which plays a major regulatory role, is dysregulated by the overexpression and inhibition of miR-150. We found that the target gene of bta-miR-150 is AKT1 and that bta-miR-150 affects AKT1 phosphorylation levels. These results showed that bta-miR-150 plays a role in adipogenic differentiation and might therefore have applications in the beef industry.

Bta-miR-376a Targeting KLF15 Interferes with Adipogenesis Signaling Pathway to Promote Differentiation of Qinchuan Beef Cattle Preadipocytes.

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.