The developmental transcriptome sequencing of bovine skeletal muscle reveals a long noncoding RNA, lncMD , promotes muscle differentiation by sponging miR-125b

Long non-coding RNA (CMR) involved in autoprotection in S. aureus mastitis in dairy cows by regulating miR-877/FOXM1

Long non-coding RNAs (LncRNAs) are dysregulated in a variety of human diseases and are highly involved in the development and progression of tumors. Studies on lncRNAs associated with cow mastitis have been lagging behind compared to humans or model animals, therefore, the aim of this study was to explore the mechanism of LncRNAs (CMR) involved in autoprotection against S. aureus mastitis in Bovine Mammary Epithelial Cells (BMECs). First, qRT-PCR was used to examine the relative expression of CMR in a S. aureus mastitis model of BMECs. Then, cell proliferation and apoptosis were detected by EdU and apoptosis assay. Finally, the targeting relationship between miRNAs and mRNA/LncRNAs was determined by dual luciferase reporter gene, qRT-PCR and western blotting techniques. The results showed that CMR was upregulated in the S. aureus mastitis model of BMECs and promoted the expression of inflammatory factors, and SiRNA-mediated CMR inhibited the proliferation of mammary epithelial cells and induced apoptosis. Mechanistically, CMR acts as a competitive endogenous RNA (ceRNA) sponge miR-877, leading to upregulation of FOXM1, a target of miR-877. Importantly, either miR-877 overexpression or FOXM1 inhibition abrogated CMR knockdown-induced apoptosis promoting cell proliferation and reducing inflammatory factor expression levels. In summary, CMR is involved in the regulation of autoprotection against S. aureus mastitis through the miR-877/FOXM1 axis in BMECs and induces immune responses in mammary tissues and cells of dairy cows, providing an important reference for subsequent prevention and control of cow mastitis and the development of targeted drugs.

Novel lncRNA regulatory elements in milk somatic cells of Holstein dairy cows associated with mastitis

Despite regulatory elements such as long non - coding RNAs representing most of the transcriptome, the functional understanding of long non - coding RNAs in relation to major health conditions including bovine mastitis is limited. This study examined the milk somatic cell transcriptome from udder quarters of 6 Holstein dairy cows to identify differentially expressed long non - coding RNAs using RNA - Sequencing. Ninety - four differentially expressed long non - coding RNAs are identified, 5 of which are previously annotated for gene name and length, 11 are annotated for gene name and 78 are novel, having no gene name or length previously annotated. Significant inflammatory response and regulation of immune response pathways (false discovery rate < 0.05) are associated with the differentially expressed long non - coding RNAs. QTL annotation analysis revealed 31 QTL previously annotated in the genomic regions of the 94 differentially expressed long non - coding RNAs, and the majority are associated with milk traits. This research provides a better understanding of long non - coding RNAs regulatory elements in milk somatic cells, which may enhance current breeding strategies for more adaptable or high mastitis resistant cattle.

Muscle fibre transition and transcriptional changes of horse skeletal muscles during traditional Mongolian endurance training

BACKGROUND

Traditional Mongolian endurance training is an effective way to improve the athletic ability of the horse for endurance events and is widely used. This incorporates aerobic exercise and intermittent fasting and these altered physiologic conditions are associated with switches between muscle fibre types.

OBJECTIVES

To better understand the adaption of horse skeletal muscle to traditional Mongolian endurance training from muscle fibre characteristics and transcriptional levels and to explore possible molecular mechanisms associated with the endurance performance of horses.

STUDY DESIGN

Before-after study.

METHODS

Muscle fibre type switches and muscle transcriptome changes in six Mongolian horses were assessed during 4 weeks of training. Transcriptomic and histochemical analyses were performed. The activities of oxidative and glycolytic metabolic enzymes were analysed and we generated deep RNA-sequencing data relating to skeletal muscles.

RESULTS

A fast-to-slow muscle fibre transition occurred in horse skeletal muscles, with a concomitant increase of oxidative enzyme activity and decreased glycolytic enzyme activity. Numerous differentially expressed genes were involved in the control of muscle protein balance and degradation. Differential alternative splicing events were also found during training which included exon-skipping events in Ttn that were associated with muscle atrophy. Differentially expressed noncoding RNAs showed connections with muscle protein balance-related pathways and fibre type specification via the post-transcriptional regulation of miRNA.

MAIN LIMITATIONS

The study focuses on horse athletic ability only from the aspect of muscular adaptation.

CONCLUSION

Traditional Mongolian endurance training-induced muscle fibre transition and metabolic and transcriptional changes. Muscle-specific non-coding RNAs could contribute to these transcriptomic changes during training.

MYOD induced lnc-MEG3 promotes porcine satellite cell differentiation via interacting with DLST

Long non-coding RNAs (lncRNAs) are involved in the process of muscle cell differentiation and play an important role. Previous studies have shown that lncRNA-MEG3 promotes the differentiation of porcine skeletal muscle satellite cells (PSCs), but the regulatory mechanism of MEG3 interaction with target protein has not been well studied. We demonstrated that MEG3 can bind dihydrolipoamide succinyltransferase (DLST) by RNA pull down and RIP-qPCR. Subsequently, knockdown and overexpression experiments showed that DLST promotes PSCs differentiation. Rescue experiments showed that the expression of DLST protein was significantly increased with MEG3 overexpression and decreased with MEG3 knockdown, while its mRNA expression was not changed. Furthermore, we have successfully predicted and validated that the transcription factor myogenic differentiation (MYOD) binds to the MEG3 core promoter though utilizing chromatin immunoprecipitation (CHIP) and luciferase reporter assays. The results indicated that MYOD acts as a transcription factor of MEG3 to promote MEG3 transcription. Knockdown of MEG3 in vivo indicated that MEG3 is involved in skeletal muscle regeneration. It is concluded that MYOD acts as a transcription factor to induce MEG3 expression. MEG3 acts as a molecular scaffold to bind and promote DLST protein expression. This paper provides a new molecular mechanism for MEG3 to promote the differentiation of PSCs.

Transcriptome Analysis of mRNA and lncRNA Related to Muscle Growth and Development in Gannan Yak and Jeryak

The production performance of Jeryak, resulting from the F1 generation of the cross between Gannan yak and Jersey cattle, exhibits a significantly superior outcome compared with that of Gannan yak. Therefore, we used an RNA-seq approach to identify differentially expressed mRNAs (DEMs) and differentially expressed lncRNAs (DELs) influencing muscle growth and development in Gannan yaks and Jeryaks. A total of 304 differentially expressed lncRNAs and 1819 differentially expressed mRNAs were identified based on the screening criteria of |log 2 FC| > 1 and FDR < 0.05. Among these, 132 lncRNAs and 1081 mRNAs were found to be down-regulated, while 172 lncRNAs and 738 mRNAs were up-regulated. GO and KEGG analyses showed that the identified DELs and DEMs were enriched in the entries of pathways associated with muscle growth and development. On this basis, we constructed an lncRNA-mRNA interaction network. Interestingly, two candidate DELs (MSTRG.16260.9 and MSTRG.22127.1) had targeting relationships with 16 (MYC, IGFBP5, IGFBP2, MYH4, FGF6, etc.) genes related to muscle growth and development. These results could provide a basis for further studies on the roles of lncRNAs and mRNAs in muscle growth in Gannan yaks and Jeryak breeds.

Bta-miR-206 and a Novel lncRNA-lncA2B1 Promote Myogenesis of Skeletal Muscle Satellite Cells via Common Binding Protein HNRNPA2B1

Skeletal muscle satellite cells (MuSCs) can proliferate, differentiate, and self-renew, and can also participate in muscle formation and muscle injury repair. Long noncoding RNAs (lncRNAs) can play an important role with the RNA binding protein and microRNAs (miRNAs) to regulate the myogenesis of bovine MuSCs, however, its molecular mechanism is still being explored. In this study, differentially expressed 301 lncRNAs were identified during the myogenic differentiation of cells based on an in vitro model of induced differentiation of bovine MuSCs using RNA sequencing (RNA-seq). Based on the ability of miR-206 to regulate myogenic cell differentiation, a new kind of lncRNA-lncA2B1 without protein-coding ability was found, which is expressed in the nucleus and cytoplasm. Subsequently, lncA2B1 inhibited cell proliferation by downregulating the expression of the proliferation marker Pax7 and promoted myogenic differentiation by upregulating the expression of the differentiation marker MyHC, whose regulatory function is closely related to miR-206. By RNA pulldown/LC-MS experiments, heterogeneous ribonucleoprotein A2/B1 (HNRNPA2B1), and DExH-Box Helicase 9 (DHX9) were identified as common binding proteins of lncA2B1 and miR-206. Overexpression of lncA2B1 and miR-206 significantly upregulated the expression level of HNRNPA2B1. Downregulation of HNRNPA2B1 expression significantly decreased the expression level of the differentiation marker MyHC, which indicates that miR-206 and lncA2B1 regulate myogenic differentiation of bovine MuSCs by acting on HNRNPA2B1. This study screened and identified a novel lncRNA-lncA2B1, which functions with miR-206 to regulate myogenesis via the common binding proteins HNRNPA2B1. The results of this study provide a new way to explore the molecular mechanisms by which lncRNAs and miRNAs regulate muscle growth and development.

Global Long Noncoding RNA Expression Profiling of MSTN and FGF5 Double-Knockout Sheep Reveals the Key Gatekeepers of Skeletal Muscle Development

Improving livestock and poultry growth rates and increasing meat production are urgently needed worldwide. Previously, we produced a myostatin (MSTN) and fibroblast growth factor 5 (FGF5) double-knockout (MF^-/-) sheep by CRISPR Cas9 system to improve meat production, and also wool production. Both MF^-/- sheep and the F1 generation (MF^+/-) sheep showed an obvious "double-muscle" phenotype. In this study, we identified the expression profiles of long noncoding RNAs (lncRNAs) in wild-type and MF^+/- sheep, then screened out the key candidate lncRNAs that can regulate myogenic differentiation and skeletal muscle development. These key candidate lncRNAs can serve as critical gatekeepers for muscle contraction, calcium ion transport and skeletal muscle cell differentiation, apoptosis, autophagy, and skeletal muscle inflammation, further revealing that lncRNAs play crucial roles in regulating muscle phenotype in MF^+/- sheep. In conclusion, our newly identified lncRNAs may emerge as novel molecules for muscle development or muscle disease and provide a new reference for MSTN-mediated regulation of skeletal muscle development.

lncRNA Sequencing of Antler Mesenchymal Tissue Revealed that the Regulatory Network of Antler Cell Proliferation and Differentiation

Antlers have been widely studied due to their unique physiological characteristics, such as rapid growth, periodic shedding and regeneration. However, little is known about how antler growth is regulated by long non-coding RNA (lncRNA). The aim of the present study was to identify the lncRNAs expression profile and explore the function of lncRNAs during the antler growth. Herein, RNA-sequencing technology (RNA-seq) was performed on the three growth periods (early developmental period: EP, middle developmental period: MP, later developmental period: LP) of male sika deer (Cervus nippon) antler, 16 differentially expressed lncRNAs (DE lncRNAs) and 11 DE lncRNAs were identified in EP vs MP and MP vs LP related to cell proliferation and cell differentiation, respectively. Finally, lncRNAs-mRNAs co-expression networks were constructed based on the identified DE lncRNAs and their potential trans-target genes. The result reveals that lncRNAs may play diverse roles in different periods of antler growth. It provides a novel perspective for revealing the molecular mechanism of antler growth.

CircCSDE1 Regulates Proliferation and Differentiation of C2C12 Myoblasts by Sponging miR-21-3p

The growth and development of skeletal muscle is regulated by many factors, and recent studies have shown that circular RNAs (circRNAs) can participate in this process. The model of porcine skeletal muscle injury was constructed to search for circRNAs that can regulate the growth and development of skeletal muscle in pigs. Using whole-transcriptome sequencing and bioinformatics analysis, a novel circRNA (circCSDE1) was screened out, which is highly expressed in skeletal muscle. Functional studies in C2C12 cells demonstrated that circCSDE1 could promote proliferation and inhibit myoblast differentiation, while opposing changes were observed by circCSDE1 knockdown. A dual-luciferase reporter assay revealed that circCSDE1 directly targeted miR-21-3p to regulate the expression of the downstream target gene (Cyclin-dependent kinase 16, CDK16). Moreover, miR-21-3p could inhibit proliferation and promote myoblast differentiation in C2C12 cells, opposite with the effects of circCSDE1. Additionally, the rescue experiments offered further evidence that circCSDE1 and its target, miR-21-3p, work together to regulate myoblast proliferation and differentiation. This study provides a theoretical basis for further understanding the regulatory mechanisms of circRNAs.

LncRNA-mRNA modules involved in goat rumen development: Insights from genome-wide transcriptome profiling

The rumen is an essential digestive and absorption organ of ruminants. During fetal life, lactation, and post-weaning period, goat rumen undergoes drastic morphological and metabolic-functional changes triggered by potential regulated genes and non-coding RNA molecules. As the essential regulatory factors, long non-coding RNAs (lncRNAs) have vital functions in various biological activities. However, their roles during rumen development are still poorly explored in ruminants. To explore the genome-wide expression profiles of lncRNAs and mRNAs in the goat rumens, we generated 5,007 lncRNAs and 19,738 mRNAs identified during the fetal and prepubertal stages by the high-throughput RNA sequencing. Notably, 365 lncRNAs and 2,877 mRNAs were considered to be differentially expressed. The weighted gene co-expression network analysis and functional analysis were performed to explore the regulatory roles of those differentially expressed molecules. The cis-and trans-target genes of differently expressed lncRNAs were enriched for pathways related to focal adhesion, cGMP-PKG signaling pathway, alpha-linolenic acid metabolism, arachidonic acid metabolism, and fat digestion and absorption. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes analyses showed that the differently expressed genes mainly participated in mitotic cytokinesis, desmosome, fatty acid degradation, cell adhesion molecules, and fatty acid metabolism. The prediction of lncRNA-mRNA interaction networks further revealed transcripts potentially involved in rumen development. The present study profiles a global overview of lncRNAs and mRNAs during rumen development. Our findings provide valuable resources for genetic regulation and molecular mechanisms of rumen development in ruminants.

Buffalo long non-coding RNA gene11007 promotes myoblasts proliferation

Buffalo meat is of good quality because it is lean and tender, and could bring significant cardiovascular benefits. The underlying difference in muscle development and meat quality is a complex and precisely orchestrated process which has been demonstrated to be regulated by long non-coding RNAs (lncRNAs). However, the regulatory role of lncRNAs in the growth and development of buffalo skeletal muscle is still unclear. In this study, the Ribo-Zero RNA-Seq method was used to explore the lncRNA expression profiles of buffalo myoblasts during the proliferation and differentiation phases. A specific set of 9,978 lncRNAs was found. By comparing the expression profiles of lncRNAs, it was found that there were 1,576 differentially expressed lncRNAs (DELs) during buffalo myoblast differentiation. Twelve DELs were chosen and subsequently verified in eight different buffalo tissues during fetal and adult stages by using qPCR. Gene11007 was found to be one of the most down-regulated lncRNAs during buffalo myoblasts differentiation and it was subsequently characterized. EdU, CCK-8, qPCR and western blotting assays showed that gene11007 promoted the proliferation of buffalo myoblasts but it had no effect on cell differentiation. Our research may enrich the genome annotations of buffalo and provide a new molecular target for the in-depth understanding of the regulation of lncRNAs in skeletal muscle.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Longissimus dorsi Skeletal Muscle of Shandong Black Cattle and Luxi Cattle

There is an increasing understanding of the possible regulatory role of long non-coding RNAs (LncRNA). Studies on livestock have mainly focused on the regulation of cell differentiation, fat synthesis, and embryonic development. However, there has been little study of skeletal muscle of domestic animals and the potential role of lncRNA. In this study, the transcriptome numbers of longissimus muscle of different beef cattle (Shandong black catle and Luxi catlle) were used to construct muscle related lncRNAs-miRNA-mRNA interaction network through bioinformatics analysis. This is helpful to clarify the molecular mechanism of bovine muscle development, and can be used to promote animal husbandry and improve animal husbandry production. According to the screening criteria of |FC|≧2 and q < 0.05, a total of 1,415 transcripts (of which 480 were LncRNAs) were differentially expressed (q < 0.05) in the different breeds. Further, we found that the most differentially expressed LncRNAs were found on chromosome 9, in which the differentially expressed LncRNAs targeted 1,164 protein coding genes (MYORG, Wnt4, PAK1, ADCY7,etc) (upstream and downstream<50 Kb). In addition, Pearson’s correlation coefficients of co-expression levels indicated a potential trans regulatory relationship between the differentially expressed LncRNAs and 43844 mRNAs (r > 0.9). The identified co-expressed mRNAs (MYORG, Dll1, EFNB2, SOX6, MYOCD, and MYLK3) are related to the formation of muscle structure, and enriched in muscle system process, strained muscle cell differentiation, muscle cell development, striated muscle tissue development, calcium signaling, and AMPK signaling. Additionally, we also found that some LncRNAs (LOC112444238, LOC101903367, LOC104975788, LOC112441863, LOC112449549, and LOC101907194) may interact with miRNAs related to cattle muscle growth and development. Based on this, we constructed a LncRNAs-miRNA-mRNA interaction network as the putative basis for biological regulation in cattle skeletal muscle. Interestingly, a candidate differential LncRNA (LOC104975788) and a protein-coding gene (Pax7) contain miR-133a binding sites and binding was confirmed by luciferase reporter assay. LOC104975788 may combined miR-133a competitively with Pax7, thus relieving the inhibitory effect of miR-133a on Pax7 to regulate skeletal muscle development. These results will provide the theoretical basis for further study of LncRNA regulation and activity in different cattle breeds.

Cellular Aquaculture: Prospects and Challenges

Aquaculture plays an important role as one of the fastest-growing food-producing sectors in global food and nutritional security. Demand for animal protein in the form of fish has been increasing tremendously. Aquaculture faces many challenges to produce quality fish for the burgeoning world population. Cellular aquaculture can provide an alternative, climate-resilient food production system to produce quality fish. Potential applications of fish muscle cell lines in cellular aquaculture have raised the importance of developing and characterizing these cell lines. In vitro models, such as the mouse C2C12 cell line, have been extremely useful for expanding knowledge about molecular mechanisms of muscle growth and differentiation in mammals. Such studies are in an infancy stage in teleost due to the unavailability of equivalent permanent muscle cell lines, except a few fish muscle cell lines that have not yet been used for cellular aquaculture. The Prospect of cell-based aquaculture relies on the development of appropriate muscle cells, optimization of cell conditions, and mass production of cells in bioreactors. Hence, it is required to develop and characterize fish muscle cell lines along with their cryopreservation in cell line repositories and production of ideal mass cells in suitably designed bioreactors to overcome current cellular aquaculture challenges.

Time-Series Clustering of lncRNA-mRNA Expression during the Adipogenic Transdifferentiation of Porcine Skeletal Muscle Satellite Cells

Skeletal muscle satellite cells (SMSCs), which are multifunctional muscle-derived stem cells, can differentiate into adipocytes. Long-chain non-coding RNA (lncRNA) has diverse biological functions, including the regulation of gene expression, chromosome silencing, and nuclear transport. However, the regulatory roles and mechanism of lncRNA during adipogenic transdifferentiation in muscle cells have not been thoroughly investigated. Here, porcine SMSCs were isolated, cultured, and induced for adipogenic differentiation. The expressions of lncRNA and mRNA at different time points during transdifferentiation were analysed using RNA-seq analysis. In total, 1005 lncRNAs and 7671 mRNAs showed significant changes in expression at differential differentiation stages. Time-series expression analysis showed that the differentially expressed (DE) lncRNAs and mRNAs were clustered into 5 and 11 different profiles with different changes, respectively. GO, KEGG, and REACTOME enrichment analyses revealed that DE mRNAs with increased expressions during the trans-differentiation were mainly enriched in the pathways for lipid metabolism and fat cell differentiation. The genes with decreased expressions were mainly enriched in the regulation of cell cycle and genetic information processing. In addition, 1883 DE mRNAs were regulated by 193 DE lncRNAs, and these genes were related to the controlling in cell cycle mainly. Notably, three genes in the fatty acid binding protein (FABP) family significantly and continuously increased during trans-differentiation, and 15, 13, and 11 lncRNAs may target FABP3, FABP4, and FABP5 genes by cis- or trans-regulation, respectively. In conclusion, these studies identify a set of new potential regulator for adipogenesis and cell fate and help us in better understanding the molecular mechanisms of trans-differentiation.

Characterization of novel lncRNA muscle expression profiles associated with meat quality in beef cattle

The aim of this study was to identify novel lncRNA differentially expressed (DE) between divergent animals for beef tenderness and marbling traits in Nellore cattle. Longissimus thoracis muscle samples from the 20 most extreme bulls (of 80 bulls set) for tenderness, tender (n = 10) and tough (n = 10) groups, and marbling trait, high (n = 10) and low (n = 10) groups were used to perform transcriptomic analysis using RNA‐Sequencing. For tenderness, 29 lncRNA were DE (p‐value ≤ 0.01) in tough beef animals in relation to tender beef animals. We observed that genic lncRNAs, for example, lncRNA_595.1, were overlapping exonic part of the PICK gene, while lncRNA_3097.2 and lncRNA_3129.5 overlapped intronic part of the genes GADL1 and PSMD6. The lncRNA associated with PICK1, GADL1, and PMD6 genes were enriched in the pathways associated with the ionotropic glutamate receptor, gamma‐aminobutyric acid synthesis, and the ubiquitin–proteasome pathway. For marbling, 50 lncRNA were DE (p‐value ≤ 0.01) in high marbling group compared with low marbling animals. The genic lncRNAs, such as lncRNA_3191.1, were overlapped exonic part of the ITGAL gene, and the lncRNA_512.1, lncRNA_3721.1, and lncRNA_41.4 overlapped intronic parts of the KRAS and MASP1 genes. The KRAS and ITGAL genes were enriched in pathways associated with integrin signaling, which is involved in intracellular signals in response to the extracellular matrix, including cell form, mobility, and mediates progression through the cell cycle. In addition, the lincRNAs identified to marbling trait were associated with several genes related to calcium binding, muscle hypertrophy, skeletal muscle, lipase, and oxidative stress response pathways that seem to play a role important in the physiological processes related to meat quality. These findings bring new insights to better understand the biology mechanisms involved in the gene regulation of these traits, which will be valuable for a further investigation of the interactions between lncRNA and mRNAs, and of how these interactions may affect meat quality traits.

Fast and slow myofiber-specific expression profiles are affected by noncoding RNAs in Mongolian horses

The heterogeneity and plasticity of muscle fibers are essential for the athletic performance of horses, mainly at the adaption of exercises and the effect on muscle diseases. Skeletal muscle fibers can be generally distinguished by their characteristics of contraction as slow and fast type myofibers. The diversity of contractile properties and metabolism enable skeletal muscles to respond to the variable functional requirements. We investigated the muscle fiber composition and metabolic enzyme activities of splenius muscle and gluteus medius muscle from Mongolian horses. The deep RNA-seq analysis of detecting differentially expressed mRNAs, lncRNAs, circRNAs and their correlation analysis from two muscles were performed. Splenius muscle and gluteus medius muscle from Mongolian horses showed a high divergence of myofiber compositions and metabolic enzyme activities. Corresponding to their phenotypic characteristics, 57 differentially expressed long noncoding RNAs and 12 differentially expressed circle RNAs were found between two muscles. The analysis results indicate multiple binding sites were detected in lncRNAs and circRNAs with myofiber-specific expressed miRNAs. Among which we found significant correlations between the above noncoding RNAs, miRNAs, their target genes, myofiber-specific developmental transcript factors, and sarcomere genes. We suggest that the ceRNA mechanism of differentially expressed noncoding RNAs by acting as miRNA sponges could be fine tuners in regulating skeletal muscle fiber composition and transition in horses, which will operate new protective measures of muscle disease and locomotor adaption for racehorses.

Whole transcriptome analysis of RNA expression profiles reveals the potential regulating action of long noncoding RNA in lactating cows fed a high concentrate diet

Subacute ruminal acidosis (SARA) is a common metabolic disease in the dairy farming industry which is usually caused by an excessive amount of high concentrate diet. SARA not only threatens animal welfare but also leads to economic losses in the farming industry. The liver plays an important role in the distribution of nutritional substances and metabolism; however, a high concentrate diet can cause hepatic metabolic disorders and liver injury. Recently, noncoding RNA has been considered as a critical regulator of hepatic disease, however, its role in the bovine liver is limited. In this study, 12 mid-lactating dairy cows were randomly assigned to a control (CON) group (40% concentrate of dry matter, n = 6) and a SARA group (60% concentrate of dry matter, n = 6). After 21 d of treatment, all cows were sacrificed, and liver tissue samples were collected. Three dairy cows were randomly selected from the CON and SARA groups respectively to perform whole transcriptome analysis. More than 20,000 messenger RNA (mRNA), 10,000 long noncoding RNA (lncRNA), 3,500 circular RNA (circRNA) and 1,000 micro RNA (miRNA) were identified. Furthermore, 43 mRNA, 121 lncRNA and 3 miRNA were differentially expressed, whereas no obvious differentially expressed circRNA were detected between the 2 groups. Gene Ontology (GO) annotation revealed that the differentially expressed genes were mainly enriched in oxidoreductase activity, stress, metabolism, the immune response, cell apoptosis, and cell proliferation. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the deferentially expressed genes were highly enriched in the phosphatidylinositol 3 kinase (PI3K)-serine/threonine kinase (AKT) signaling pathway (P < 0.05). According to KEGG pathway analysis, the differentially expressed lncRNA (DElncRNA) target genes were mainly related to proteasomes, peroxisomes, and the hypoxia-inducible factor-1 signaling pathway (P < 0.005). Further bioinformatics and integrative analyses revealed that the lncRNA were strongly correlated with mRNA; therefore, it is reasonable to speculate that lncRNA potentially play important roles in the liver dysfunction induced by SARA. Our study provides a valuable resource for future investigations on the mechanisms of SARA to facilitate an understanding of the importance of lncRNA, and offer functional RNA information.

Non-Coding RNAs as Regulators of Myogenesis and Postexercise Muscle Regeneration

miRNAs and lncRNAs do not encode proteins, but they play an important role in the regulation of gene expression. They differ in length, biogenesis, and mode of action. In this work, we focus on the selected miRNAs and lncRNAs involved in the regulation of myogenesis and muscle regeneration. We present selected miRNAs and lncRNAs that have been shown to control myogenic differentiation and show that manipulation of their levels could be used to improve myogenic differentiation of various types of stem and progenitor cells. Finally, we discuss how physical activity affects miRNA and lncRNA expression and how it affects muscle well-being.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

The circular RNA circCPE regulates myoblast development by sponging miR-138

Background

Skeletal muscle development, a long-term and complex process, is controlled by a set of the myogenic genes. Circular RNAs (circRNAs), a class of noncoding RNA, have been shown to regulate various biological processes. Recent studies indicate circRNAs may be involved in myogenesis, but the role and regulatory mechanism of circRNAs in myogenesis is largely unknown. In the present study, circCPE was firstly found to promote the bovine myoblast proliferation and inhibit cell apoptosis and differentiation by influencing the expression of FOXC1 in a miR138-mediated manner. And in vivo experiments revealed that overexpression of circCPE attenuates skeletal muscle regeneration.

Results

We identified a novel circular RNA circCPE by analyzing circRNAs sequencing data of bovine muscle tissue. Sequencing verification, RNase R treatment and Actinomycin D treatment confirmed the circular nature of circCPE in bovine muscle. Functional assays showed that overexpression of circCPE could inhibit bovine myoblast apoptosis and differentiation, as well as facilitate cell proliferation. Moreover, in vivo experiments revealed that overexpression of circCPE attenuates skeletal muscle regeneration. In consideration of circRNA action as miRNAs sponge, we found that circCPE harbors miR-138 binding sites and absorbed miR-138. Mechanistically, the rescue experiments showed that the overexpression of circCPE can counteract the inhibitory effect of miR-138 on the cell proliferation and the accelerated effects on the differentiation and apoptosis. Subsequently, we found that circCPE sequester the inhibitory effect of miR-138 on FOXC1 so as to involve in myogenesis.

Conclusions

Collectively, we constructed a novel circCPE/miR-138/FOXC1 regulatory network in bovine myogenesis, which further provide stronger evidence that circRNA involved in muscle development acting as miRNA sponge.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00618-7.

Transcriptome Analysis Reveals the Profile of Long Non-coding RNAs During Chicken Muscle Development

The developmental complexity of muscle arises from elaborate gene regulation. Long non-coding RNAs (lncRNAs) play critical roles in muscle development through the regulation of transcription and post-transcriptional gene expression. In chickens, previous studies have focused on the lncRNA profile during the embryonic periods, but there are no studies that explore the profile from the embryonic to post-hatching period. Here, we reconstructed 14,793 lncRNA transcripts and identified 2,858 differentially expressed lncRNA transcripts and 4,282 mRNAs from 12-day embryos (E12), 17-day embryos (E17), 1-day post-hatch chicks (D1), 14-day post-hatch chicks (D14), 56-day post-hatch chicks (D56), and 98-day post-hatch chicks (D98), based on our published RNA-seq datasets. We performed co-expression analysis for the differentially expressed lncRNAs and mRNAs, using STEM, and identified two profiles with opposite expression trends: profile 4 with a downregulated pattern and profile 21 with an upregulated pattern. The cis- and trans-regulatory interactions between the lncRNAs and mRNAs were predicted within each profile. Functional analysis of the lncRNA targets showed that lncRNAs in profile 4 contributed to the cell proliferation process, while lncRNAs in profile 21 were mainly involved in metabolism. Our work highlights the lncRNA profiles involved in the development of chicken breast muscle and provides a foundation for further experiments on the role of lncRNAs in the regulation of muscle development.

A novel lncRNA promotes myogenesis of bovine skeletal muscle satellite cells via PFN1-RhoA/Rac1

Myogenesis, the process of skeletal muscle formation, is a highly coordinated multistep biological process. Accumulating evidence suggests that long non-coding RNAs (lncRNAs) are emerging as a gatekeeper in myogenesis. Up to now, most studies on muscle development-related lncRNAs are mainly focussed on humans and mice. In this study, a novel muscle highly expressed lncRNA, named lnc23, localized in nucleus, was found differentially expressed in different stages of embryonic development and myogenic differentiation. The knockdown and over-expression experiments showed that lnc23 positively regulated the myogenic differentiation of bovine skeletal muscle satellite cells. Then, TMT 10-plex labelling quantitative proteomics was performed to screen the potentially regulatory proteins of lnc23. Results indicated that lnc23 was involved in the key processes of myogenic differentiation such as cell fusion, further demonstrated that down-regulation of lnc23 may inhibit myogenic differentiation by reducing signal transduction and cell fusion among cells. Furthermore, RNA pulldown/LC-MS and RIP experiment illustrated that PFN1 was a binding protein of lnc23. Further, we also found that lnc23 positively regulated the protein expression of RhoA and Rac1, and PFN1 may negatively regulate myogenic differentiation and the expression of its interacting proteins RhoA and Rac1. Hence, we support that lnc23 may reduce the inhibiting effect of PFN1 on RhoA and Rac1 by binding to PFN1, thereby promoting myogenic differentiation. In short, the novel identified lnc23 promotes myogenesis of bovine skeletal muscle satellite cells via PFN1-RhoA/Rac1.

Lnc-ORA interacts with microRNA-532-3p and IGF2BP2 to inhibit skeletal muscle myogenesis

Skeletal muscle is one of the most important organs of the animal body. Long noncoding RNAs play a crucial role in the regulation of skeletal muscle development via several mechanisms. We recently identified obesity-related lncRNA (lnc-ORA) in a search for long noncoding RNAs that influence adipogenesis, finding it impacted adipocyte differentiation by regulating the PI3K/protein kinase B/mammalian target of rapamycin pathway. However, whether lnc-ORA has additional roles, specifically in skeletal muscle myogenesis, is not known. Here, we found that lnc-ORA was significantly differentially expressed with age in mouse skeletal muscle tissue and predominantly located in the cytoplasm. Overexpression of lnc-ORA promoted C2C12 myoblast proliferation and inhibited myoblast differentiation. In contrast, lnc-ORA knockdown repressed myoblast proliferation and facilitated myoblast differentiation. Interestingly, silencing of lnc-ORA rescued dexamethasone-induced muscle atrophy in vitro. Furthermore, adeno-associated virus 9–mediated overexpression of lnc-ORA decreased muscle mass and the cross-sectional area of muscle fiber by upregulating the levels of muscle atrophy–related genes and downregulating the levels of myogenic differentiation–related genes in vivo. Mechanistically, lnc-ORA inhibited skeletal muscle myogenesis by acting as a sponge of miR-532-3p, which targets the phosphatase and tensin homolog gene; the resultant changes in phosphatase and tensin homolog suppressed the PI3K/protein kinase B signaling pathway. In addition, lnc-ORA interacted with insulin-like growth factor 2 mRNA-binding protein 2 and reduced the stability of myogenesis genes, such as myogenic differentiation 1 and myosin heavy chain. Collectively, these findings indicate that lnc-ORA could be a novel underlying regulator of skeletal muscle development.

Differential expression and prediction of function of lncRNAs in the ovaries of low and high fecundity Hanper sheep

Litter size is an important trait that determines the production efficiency of sheep bred for meat. Its detailed investigation can reveal the molecular mechanisms that control the fecundity of sheep and possibly accelerate the breeding process of new varieties of sheep that have high prolificacy. Long non-coding RNAs (lncRNAs) have proven to be an important factor in the regulation of follicular development. However, the mechanisms by which lncRNAs regulate litter size in sheep remain unclear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either monotocous (M) or polytocous (P; FM, FP, LM and LP groups) were collected and sequenced to identify differentially expressed lncRNAs and predict their function. The results indicate that the number of up- and down-regulated lncRNAs in the follicular phase (FM vs. FP) was 95 and 111 and 109 and 49, respectively, in the luteal phase (LM vs. LP). The functional enrichment of the different lncRNAs coexpressed with mRNA was analysed. The results demonstrated that the KISS1-GnRH-LH/FSH-E2 and EGF-EGFR-RAS-PI3K signalling pathways promoted the initiation of the primordial period, follicular development and ovulation in the follicular phase (FM vs. FP). During the luteal phase (LM vs. LP), the production and development of the corpus luteum in ewes was influenced by the KITLG-KIT/FGF-FGFR/HGF-MET-RAS-ERK signalling pathway. STEM clustering functional enrichment analysis of the differentially expressed lncRNAs indicated that profile11 was principally enriched in the Cytokine-Jak-STAT, PDGF-PDGFR-PI3K and KITLG-KIT-RAS-ERK signalling pathways. By analysis of the differential expression of the lncRNAs and their expression in each group, lncRNAs Xist (loc101112291) and Gtl2 (loc101123329) were found to be highly expressed, suggesting that regulation of follicular development was mediated through methylation processes.

Circular RNA circMYL1 Inhibit Proliferation and Promote Differentiation of Myoblasts by Sponging miR-2400

Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs (ncRNAs) involved in regulating skeletal muscle development by sponging miRNAs. In this study, we found that the circMYL1 expression was down-regulated during myoblast proliferation, while gradually up-regulated in myoblast differentiation. The potential role of circMYL1 was identified in the proliferation of bovine myoblast through mRNA and protein expression of proliferation marker genes (PCNA, CyclinD1, and CDK2), cell counting kit-8 assay, flow cytometry analysis, and 5-ethynyl 2′-deoxyuridine (EdU) assay. Analysis of the expression of differentiation marker genes (MyoD, MyoG, and MYH2) and immunofluorescence of Myosin heavy chain (MyHC) was used to assess cell differentiation. The proliferation analysis revealed that circMYL1 inhibited the proliferation of bovine primary myoblast. Furthermore, the differentiation analysis demonstrated that circMYL1 promoted the differentiation of bovine primary myoblast. The luciferase screening and RNA immunoprecipitation (RIP) assays found that circMYL1 could have interaction with miR-2400. Additionally, we demonstrated that miR-2400 promoted proliferation and inhibited differentiation of bovine primary myoblast, while circMYL1 may eliminate the effects of miR-2400, as showed by rescue experiments. Together, our results revealed that a novel circular RNA of circMYL1 could inhibit proliferation and promote differentiation of myoblast by sponging miR-2400.

MicroRNA bta-miR-365-3p inhibits proliferation but promotes differentiation of primary bovine myoblasts by targeting the activin A receptor type I

Background

MicroRNAs act as post-transcriptional regulators that repress translation or degrade mRNA transcripts. Each microRNA has many mRNA targets and each mRNA may be targeted by several microRNAs. Skeletal muscles express a plethora of microRNA genes that regulate muscle development and function by controlling the expression of protein-coding target genes. To expand our understanding of the role of microRNA, specifically bta-miR-365-3p, in muscle biology, we investigated its functions in regulating primary bovine myoblast proliferation and differentiation.

Results

Firstly, we found that bta-miR-365-3p was predominantly expressed in skeletal muscle and heart tissue in Chinese Qinchuan beef cattle. Quantitative PCR and western blotting results showed that overexpression of bta-miR-365-3p significantly reduced the expression levels of cyclin D1 (CCND1), cyclin dependent kinase 2 (CDK2) and proliferating cell nuclear antigen (PCNA) but stimulated the expression levels of muscle differentiation markers, i.e., MYOD1, MYOG at both mRNA and protein level. Moreover, downregulation of bta-miR-365-3p increased the expression of CCND1, CDK2 and PCNA but decreased the expression of MYOD1 and MYOG at both mRNA and protein levels. Furthermore, flow cytometry, EdU proliferation assays and immunostaining results showed that increased levels of bta-miR-365-3p suppressed cell proliferation but promoted myotube formation, whereas decreased levels of bta-miR-365-3p resulted in the opposite consequences. Finally, we identified that activin A receptor type I (ACVR1) could be a direct target of bta-miR-365-3p. It was demonstrated that bta-miR-365-3p can bind to the 3’UTR of ACVR1 gene to regulate its expression based on dual luciferase gene reporter assays. Consistently, knock-down of ACVR1 was associated with decreased expressions of CDK2, CCND1 and PCNA but increased expression of MYOG and MYOD1 both at mRNA and protein level.

Conclusion

Collectively, these data suggested that bta-miR-365-3p represses proliferation but promotes differentiation of bovine myoblasts through several biological mechanisms involving downregulation of ACVR1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-020-00528-0.

RIP-Seq of EZH2 Identifies TCONS-00036665 as a Regulator of Myogenesis in Pigs

Enhancer of zeste homolog 2 (EZH2) is the catalytic subunit of polycomb repressive complex 2 and contains a SET domain that catalyzes histone H3 trimethylation on lysine 27 (H3K27me3) to generate an epigenetic silencing mark. EZH2 interacts with transcription factors or RNA transcripts to perform its function. In this study, we applied RNA immunoprecipitation sequencing and long intergenic non-coding RNA (lincRNA) sequencing methods to identify EZH2-binding lincRNAs. A total of 356 novel EZH2-binding lincRNAs were identified by bioinformatics analysis and an EZH2-binding lincRNA TCONS-00036665 was characterized. TCONS-00036665 promoted pig skeletal satellite cell proliferation but inhibited cell differentiation, and this function was conserved between pigs and mice. Further mechanistic studies indicated that TCONS-00036665 can bind to EZH2 and recruits EZH2 to the promoters of the target genes p21, MyoG, and Myh4, which leads to the enrichment of H3K27me3 and the repression of target gene expression and pig myogenesis. In conclusion, the lincRNA TCONS-00036665 regulates pig myogenesis through its interaction with EZH2.

Effects of miR-125b-5p on Preadipocyte Proliferation and Differentiation in Chicken

Our previous studies have shown that miR-125b-5p was highly expressed and significantly upregulated during abdominal fat deposition in chickens. However, the role of miR-125b in the regulation of adipogenesis is not clear in chickens. Therefore, we evaluated the effects of miR-125b-5p on preadipocyte proliferation and differentiation and the interaction between miR-125b-5p and the acyl-CoA synthetase bubblegum family member 2 (ACSBG2) gene in adipogenesis in chicken abdominal adipose tissue. Here, transfection tests of miR-125b-5p mimic/inhibitor were performed in preadipocytes, and the effects of miR-125b-5p on preadipocytes proliferation and differentiation were analyzed. The target site of miR-125b-5p in the 3'UTR (untranslated region) of ACSBG2 were verified by a luciferase reporter assay. Our results showed that miR-125b-5p overexpression inhibited proliferation and reduced the number of cells in S phase and G2/M phase in preadipocytes; conversely, miR-125b-5p inhibition promoted the proliferation and increased the number of cells in S phase and G2/M phase. In adipocytes after induction, miR-125b-5p overexpression led to a notable increase in the accumulation of lipid droplets as well as in the concentration of triglycerides, while miR-125b-5p inhibition had the opposite effect. Furthermore, miR-125b-5p could directly bind to the 3'UTR of ACSBG2, and its overexpression could significantly repress the mRNA and protein expression of ACSBG2. These results indicate that miR-125b-5p can inhibit preadipocyte proliferation and can promote preadipocyte differentiation to affect adipogenesis in chicken abdominal adipose tissues, at least partially by downregulating ACSBG2.

Competing Endogenous RNA Networks as Biomarkers in Neurodegenerative Diseases

Protein aggregation is classically considered the main cause of neuronal death in neurodegenerative diseases (NDDs). However, increasing evidence suggests that alteration of RNA metabolism is a key factor in the etiopathogenesis of these complex disorders. Non-coding RNAs are the major contributor to the human transcriptome and are particularly abundant in the central nervous system, where they have been proposed to be involved in the onset and development of NDDs. Interestingly, some ncRNAs (such as lncRNAs, circRNAs and pseudogenes) share a common functionality in their ability to regulate gene expression by modulating miRNAs in a phenomenon known as the competing endogenous RNA mechanism. Moreover, ncRNAs are found in body fluids where their presence and concentration could serve as potential non-invasive biomarkers of NDDs. In this review, we summarize the ceRNA networks described in Alzheimer's disease, Parkinson's disease, multiple sclerosis, amyotrophic lateral sclerosis and spinocerebellar ataxia type 7, and discuss their potential as biomarkers of these NDDs. Although numerous studies have been carried out, further research is needed to validate these complex interactions between RNAs and the alterations in RNA editing that could provide specific ceRNET profiles for neurodegenerative disorders, paving the way to a better understanding of these diseases.

miR-125b prevent the progression of esophageal squamous cell carcinoma through the p38-MAPK signaling pathway

Background

To examine the clinical significance of miR-125b in esophageal squamous cell carcinoma (ESCC) and to research the effect of miR-125b on the biological function of ESCC cells and the relevant underlying mechanism.

Methods

The expression of miR-125b in ESCC tissues and cell lines were discovered by RT-PCR assay. The interrelation between miR-125b expression and clinicopathological parameters and the forecasting of ESCC patients were analyzed. CCK-8 method and Transwell methods were used to detect the increased growth, shifting, and irruption of ESCC cells. Bioinformatics analysis was applied to forecast the possible target genes of miR-125b and verified through dual-luciferase reporter gene assay. After that, the expression of p38-MAPK mRNA and protein were found out by RT-PCR and Western blot.

Results

The expression of miR-125b was down-regulated in ESCC tissues and cell lines (P<0.05). And the expression of miR-125b was closely about tumor differentiation, TNM level, and lymph node metastasis in ESCC patients. The low miR-125b formulation was closely related to rough forecasting in ESCC patients. Large scale expression of miR-125b can effectively decrease the acceleration, shifting, and irrupting strengths of ESCC cells. Bioinformatics analysis showed p38-MAPK was forecasted to be a potential mark of miR-125b, which was confirmed by dual luciferase assay, and extreme expression of miR-125b can stop the expression of p38-MAPK mRNA and protein.

Conclusions

miR-125b is down-regulated in ESCC. Moreover, its expression level is significant concerning tumor progression and prognosis in patients with ESCC. MiR-125b can stop the high growth and shifting of ESCC cells having p38-MAPK at target.

Long Non-coding RNA H19 Regulates Porcine Satellite Cell Differentiation Through miR-140-5p/SOX4 and DBN1

The H19 gene promotes skeletal muscle differentiation in mice, but the regulatory models and mechanisms of myogenesis regulated by H19 are largely unknown in pigs. Therefore, the regulatory modes of H19 in the differentiation of porcine skeletal muscle satellite cells (PSCs) need to be determined. We observed that H19 gene silencing could decrease the expressions of the myogenin (MYOG) gene, myogenic differentiation (MYOD), and myosin heavy chain (MYHC) in PSCs. Therefore, we constructed and sequenced 12 cDNA libraries of PSCs after knockdown of H19 at two differentiation time points to analyze the transcriptome differences. A total of 11,419 differentially expressed genes (DEGs) were identified. Among these DEGs, we found through bioinformatics analysis and protein interaction experiment that SRY-box transcription factor 4 (SOX4) and Drebrin 1 (DBN1) were the key genes in H19-regulated PSC differentiation. Functional analysis shows that SOX4 and DBN1 promote PSC differentiation. Mechanistically, H19 regulates PSC differentiation through two different pathways. On the one hand, H19 functions as a molecular sponge of miR-140-5p, which inhibits the differentiation of PSCs, thereby modulating the derepression of SOX4. On the other hand, H19 regulates PSC differentiation through directly binding with DBN1. Furthermore, MYOD binds to the promoters of H19 and DBN1. The knockdown of MYOD inhibits the expression of H19 and DBN1. We determined the function of H19 and provided a molecular model to elucidate H19’s role in regulating PSC differentiation.

Transcriptome analysis identified long non-coding RNAs involved in the adaption of yak to high-altitude environments

The mechanisms underlying yak adaptation to high-altitude environments have been investigated using various methods, but no report has focused on long non-coding RNA (lncRNA). In the present study, lncRNAs were screened from the gluteus transcriptomes of yak and their transcriptional levels were compared with those in Sanjiang cattle, Holstein cattle and Tibetan cattle. The potential target genes of the differentially expressed lncRNAs between species/strains were predicted using cis and trans models. Based on cis-regulated target genes, no KEGG pathway was significantly enriched. Based on trans-regulated target genes, 11 KEGG pathways in relation to energy metabolism and three KEGG pathways associated with muscle contraction were significantly enriched. Compared with cattle strains, transcriptional levels of acyl-CoA dehydrogenase, acyl-CoA-binding protein, 3-hydroxyacyl-CoA dehydrogenase were relatively higher and those of glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate mutase 1, pyruvate kinase and lactate/malate dehydrogenase were relatively lower in yak, suggesting that yaks activated fatty acid oxidation but inhibited glucose oxidation and glycolysis. Besides, NADH dehydrogenase and ATP synthase showed lower transcriptional levels in yak than in cattle, which might protect muscle tissues from deterioration caused by reactive oxygen species (ROS). Compared with cattle strains, the higher transcriptional level of glyoxalase in yak might contribute to dicarbonyl stress resistance. Voltage-dependent calcium channel/calcium release channel showed a lower level in yak than in cattle strains, which could reduce the Ca²⁺ influx and subsequently decrease the risk of hypertension. However, levels of EF-hand and myosin were higher in yak than in cattle strains, which might enhance the negative effects of reduced Ca²⁺ on muscle contraction. Overall, the present study identified lncRNAs and proposed their potential regulatory functions in yak.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

Comparative Transcriptomic Analysis of the Pituitary Gland between Cattle Breeds Differing in Growth: Yunling Cattle and Leiqiong Cattle

The hypothalamic-pituitary-thyroid (HPT) axis hormones regulate the growth and development of ruminants, and the pituitary gland plays a decisive role in this process. In order to identify pivotal genes in the pituitary gland that could affect the growth of cattle by regulating the secretion of hormones, we detected the content of six HPT hormones related to growth in the plasma of two cattle breeds (Yunling and Leiqiong cattle, both also known as the zebu cattle) with great differences in growth and compared the transcriptome data of their pituitary glands. Our study found that the contents of GH, IGF, TSH, thyroxine, triiodothyronine, and insulin were significantly different between the two breeds, which was the main cause of the difference in growth; 175 genes were identified as differentially expressed genes (DEGs). Functional association analyses revealed that DEGs were mainly involved in the process of transcription and signal transduction. Combining the enrichment analysis and protein interaction analysis, eight DEGs were predicted to control the growth of cattle by affecting the expression of growth-related hormones in the pituitary gland. In summary, our results suggested that SLC38A1, SLC38A3, DGKH, GNB4, GNAQ, ESR1, NPY, and GAL are candidates in the pituitary gland for regulating the growth of Yunling and Leiqiong cattle by regulating the secretion of growth-related hormones. This study may help researchers further understand the growth mechanisms and improve the artificial selection of zebu cattle.

lncRNA IGF2 AS Regulates Bovine Myogenesis through Different Pathways

The role of long non-coding RNA (lncRNA) in the regulation of bovine skeletal muscle development remains poorly understood. The present study investigated the function and regulatory mechanism of a novel lncRNA, insulin-like growth factor 2 antisense transcript (IGF2 AS), in bovine myoblast proliferation and differentiation. Gain or loss of IGF2 AS was performed using an expression plasmid or small interfering RNA (siRNA), respectively. Bovine myoblasts were used to investigate the biological function and mechanisms of IGF2 AS in vitro. Results were conjointly analyzed by celluar and molecular biology experiments as well as bioinformatics. Functionally, IGF2 AS could promote proliferation and differentiation of bovine myoblasts. The preliminary mechanism suggests, on the one hand, that IGF2 AS could complement the IGF2 gene intron region and affect the stability and expression of IGF2 mRNA. On the other hand, RNA pull-down and immunoprecipitation assays demonstrated that IGF2 AS could directly bind to the interleukin enhancer binding factor 3 (ILF3) protein and maybe partly though it to regulate myogenesis. In conclusion, the novel identified lncRNA IGF2 AS promoted proliferation and differentiation of bovine myoblasts through various pathways.

Identification of novel alternative splicing of bovine lncRNA lncFAM200B and its effects on preadipocyte proliferation

Adipogenesis is closely related to human health, livestock growth, and meat quality. A previous study identified that bovine lncFAM200B promoter has high activity in 3T3-L1 mice preadipocytes. Thus, lncFAM200B was a candidate gene for regulating adipogenesis. This study aimed to uncover the role of lncFAM200B in bovine adipogenesis and identify novel genetic variations within the bovine lncFAM200B gene. An expression analysis found that lncFAM200B was expressed higher in fat than that in muscle, but the difference was not related to the total methylation level of the promoter active region. Moreover, the expression of lncFAM200B exhibited a significant positive correlation with the expression of C/EBPa during bovine adipocyte differentiation. To uncover the function of lncFAM200B, the full-length lncFAM200B was cloned, and four kinds of transcript variants were found. Protein-coding potential prediction and prokaryotic expression system analysis showed that these four transcript variants were noncoding RNAs. The quantitative reverse-transcription polymerase chain reaction and 5-ethynyl-2'-deoxyuridine assay showed that the transcript variants decreased the messenger RNA expression of Cyclin D1 and inhibited the proliferation of bovine preadipocytes. Considering the important role of lncFAM200B in adipogenesis, we identified genetic variations in lncFAM200B. Three single-nucleotide polymorphisms (SNPs) were revealed, and two of them (SNP1 and SNP3) were associated with Nanyang cattle body measurement traits. In conclusion, this study found that bovine lncFAM200B inhibited preadipocyte proliferation, and two genetic variations of lncFAM200B could be used in cattle breeding.

A novel lncRNA, lnc403, involved in bovine skeletal muscle myogenesis by mediating KRAS/Myf6

Skeletal muscle, the most abundant and plasticity tissue in mammals, is essential for various functions such as movement, breathing, maintaining posture and metabolism. Myogenesis is a complex and precise process, which is regulated by the sequential expression of multiple transcription factors, and accumulating evidence have confirmed that multiple lncRNAs are involved in muscle development as the important transcriptional regulator. In this study, a novel lncRNA, named lnc403 was obtained, with a full-length 2689 bp, which had poor coding ability and was mainly expressed in the nucleus of myoblasts and myotubes. The expression of lnc403 was significantly different in the proliferation and differentiation stages of muscle cells. Then we successfully constructed lnc403 loss/gain-function cell models by transfecting silnc403 and pCDNA3.1-EGFP-lnc403 into satellite cells, respectively; and found that lnc403 inhibited skeletal muscle satellite cell differentiation but had no significant effect on cell proliferation, either in the case of lnc403 knockdown or overexpression. In order to further screen the target factors regulated by lncRNA in the process of myogenic differentiation, the RNA-pull down, mass spectrometry and bioinformatics analysis were performed. The results showed that lnc403 negatively regulated the expression of the adjacent gene Myf6 and positively regulated interaction proteins KRAS expression. The above results indicate that lnc403 affects skeletal muscle cell differentiation by affecting the expression of nearby genes and interacting proteins, implying lnc403 might participate in the bovine myoblasts differentiation through multi-pathway network regulation mode. This study provides a new perspective for further understanding of the regulation mechanism of lncRNAs on bovine myogenic process.

Genome-Wide Analysis Reveals Changes in Polled Yak Long Non-coding RNAs in Skeletal Muscle Development

Long non-coding RNAs (lncRNAs) have been extensively studied in recent years. Numerous lncRNAs have been identified in mice, rats, and humans, some of which play important roles in muscle formation and development. However, little is known about lncRNA regulators that affect muscle development in yak (Bos grunniens). LncRNA expression during skeletal muscle development in yak was analyzed by RNA sequencing at three development stages: 3 years (group A), 6 months (group M), and 90-day-old fetuses (group E). A total of 1180 lncRNAs were identified in the three development stages. Compared with group E, 154 were upregulated and 130 were downregulated in group A. Compared with group A, 31 were upregulated and 29 were downregulated in group M. Compared with group E, 147 were upregulated and 149 were downregulated in group M (padj < 0.001, |log2FC| > 1.2). In addition, functional annotation analysis based on gene ontology (GO) and the Kyoto protocol encyclopedia of genes and genomes (KEGG) database showed that differentially expressed lncRNAs (DElncRNAs) were cis–trans target genes. The results showed that DElncRNAs were mainly involved in PI3K-Akt signaling pathway, focal adhesion, MAPK signaling pathway, apoptosis, and p53 signaling pathway. Furthermore, RTL1, IGF2, MEF2C, Pax7, and other well-known muscle development regulators were included in a co-expression network of differentially expressed target genes and lncRNAs. These data will help to further clarify the function of lncRNAs in the different stages of skeletal muscle developmental in yak.

MEG3 Promotes Differentiation of Porcine Satellite Cells by Sponging miR-423-5p to Relieve Inhibiting Effect on SRF

Although thousands of long noncoding RNAs (lncRNAs) have been identified in porcine growth and development, the regulation mechanisms of functional lncRNAs have not been well explored. In this study, using 5′- and 3′-rapid amplification of cDNA ends (RACE) assays, we obtained two different variants of lncRNA maternally expressed gene 3 (MEG3), namely, MEG3 v1 and MEG3 v2, that were both highly expressed in porcine skeletal muscle and in the early stage of the differentiation of porcine satellite cells. Moreover, we identified the core transcript MEG3 v2. Functional analyses showed that MEG3 overexpression could effectively arrest myoblasts in the G1 phase, inhibit DNA replication, and promote myoblast differentiation, whereas MEG3 knockdown resulted in the opposite effects. Interestingly, the expression of serum response factor (SRF), a crucial transcription factor for myogenesis process, remarkably increased and decreased in mRNA and protein levels with the respective overexpression and knockdown of MEG3. Dual luciferase reporter assay showed that MEG3 could attenuate the decrease of luciferase activity of SRF induced by miR-423-5p in a dose-dependent manner. MEG3 overexpression could relieve the inhibitory effect on SRF and myoblast differentiation induced by miR-423-5p. In addition, results of RNA immunoprecipitation analysis suggested that MEG3 could act as a ceRNA for miR-423-5p. Our findings initially established a novel connection among MEG3, miR-423-5p, and SRF in porcine satellite cell differentiation. This novel role of MEG3 may shed new light on understanding of molecular regulation of lncRNA in porcine myogenesis.

Roles of lncRNAs and circRNAs in regulating skeletal muscle development

The multistep biological process of myogenesis is regulated by a variety of myoblast regulators, such as myogenic differentiation antigen, myogenin, myogenic regulatory factor, myocyte enhancer factor2A-D and myosin heavy chain. Proliferation and differentiation during skeletal muscle myogenesis contribute to the physiological function of muscles. Certain non-coding RNAs, including long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs), are involved in the regulation of muscle development, and the aberrant expressions of lncRNAs and circRNAs are associated with muscular diseases. In this review, we summarize the recent advances concerning the roles of lncRNAs and circRNAs in regulating the developmental aspects of myogenesis. These findings have remarkably broadened our understanding of the gene regulation mechanisms governing muscle proliferation and differentiation, which makes it more feasible to design novel preventive, diagnostic and therapeutic strategies for muscle disorders.

Bta-miR-885 promotes proliferation and inhibits differentiation of myoblasts by targeting MyoD1

The proliferation and differentiation of myoblasts are essential for the regeneration and development of skeletal muscles. However, the process of skeletal muscle development in cattle is complex and needs to be further investigated. The microRNAs (miRNAs) are endogenous, small noncoding RNAs that play a critical role during skeletal muscle development. In this study, we evaluated the function of miR-885 in muscle development in cattle. The results found that the expression of miR-885 was gradually upregulated during myoblast proliferation, whereas progressively downregulated during myoblast differentiation. The overexpression of miR-885 promoted cell proliferation of myoblast in cattle. Moreover, we further noted that the overexpression miR-885 triggered the expression level of various marker genes involved in cell proliferation, including proliferating cell nuclear antigen (PCNA), cyclin-dependent kinase 2 (CDK2), and cyclin B1 (CCNB1). Furthermore, it was observed that overexpression of miR-885 inhibited cell differentiation, and significantly decreased messenger RNA and protein expression levels of myogenic differentiation 1 (MyoD1) and myogenin (MyoG) in primary bovine myoblasts. Moreover, the miR-885 inhibitor revealed that miR-885 inhibited cell proliferation and promoted cell differentiation. In addition, the overexpression of miR-885 markedly decreased MyoD1 expression in primary bovine myoblasts. The luciferase reporter assay, quantitative real-time polymerase chain reaction, and western blot (WB) further indicated that miR-885 directly binding to 3' UTR of MyoD1 gene during transcriptional regulation. Conclusively, these results signified that miR-885 could be critical for the proliferation and differentiation in primary bovine myoblast cells by targeting the MyoD1 gene in cattle.

Genome-wide analysis of mRNAs and lncRNAs in Mycoplasma bovis infected and non-infected bovine mammary gland tissues

Mycoplasma bovis (M. bovis) causes diseases such as arthritis, pneumonia, abortion, and mastitis, leading to great losses in the bovine dairy industries. RNA types such as messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) play significant roles in regulating the immune responses triggered by bacteria. The expression profiles of mRNA and lncRNA as they occur in bovine mammary gland tissues infected with M. bovis are still not well understood. To illuminate this issue, transcription analysis of mRNA and LncRNAs were conducted on the mammary gland tissues belonging to Holstein cattle infected and not infected with M. bovis. The analysis revealed 1310 differentially expressed mRNAs and 57 differentially expressed lncRNAs in the bovine mammary gland tissues infected and not infected with M. bovis. In addition, 392 novel lncRNAs were detected, 19 of which were differentially expressed. Gene ontology analysis reveals that differentially expressed mRNAs and lncRNAs play significant roles in such vital biological pathways as metabolic pathways, T-cell receptor signaling, TGF-beta signaling, pathways in cancer, PI3K-Akt signaling, NF-kappa B signaling, mTOR signaling, and apoptosis, including in the immune response to cancer. Based on our literature review, this study is the first genome-wide lncRNA research conducted on bovine mammary gland tissues infected with M. bovis. Our results provide bovine mammary gland lncRNA and mRNA resources to understand their roles in the regulation of the immune response against the agent M. bovis in bovine mammary gland tissues.

MIR221HG Is a Novel Long Noncoding RNA that Inhibits Bovine Adipocyte Differentiation

Adipogenesis is a complicated but precisely orchestrated process mediated by a series of transcription factors. Our previous study has identified a novel long noncoding RNA (lncRNA) that was differentially expressed during bovine adipocyte differentiation. Because this lncRNA overlaps with miR-221 in the genome, it was named miR-221 host gene (MIR221HG). The purpose of this study was to clone the full length of MIR221HG, detect its subcellular localization, and determine the effects of MIR221HG on bovine adipocyte differentiation. The 5′ rapid amplification of cDNA ends (RACE) and 3′ RACE analyses demonstrated that MIR221HG is a transcript of 1064 nucleotides, is located on the bovine X chromosome, and contains a single exon. Bioinformatics analyses suggested that MIR221HG is an lncRNA and the promoter of MIR221HG includes the binding consensus sequences of the forkhead box C1 (FOXC1) and krüppel-like factor5 (KLF5). The semi-quantitative PCR and quantitative real-time PCR (qRT-PCR) of nuclear and cytoplasmic fractions revealed that MIR221HG mainly resides in the nucleus. Inhibition of MIR221HG significantly increased adipocyte differentiation, as indicated by a dramatic increment in the number of mature adipocytes and in the expression of the respective adipogenic markers, peroxisome proliferator-activated receptor γ (PPARγ), CCAAT/enhancer-binding protein α (C/EBPα), and fatty acid binding protein 4 (FABP4). Our results provide a basis for elucidating the mechanism by which MIR221HG regulates adipocyte differentiation.

A Novel Long Noncoding RNA, lncR-125b, Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Sponging miR-125b

Long noncoding RNAs (lncRNAs) have emerged as essential regulators of skeletal myogenesis, but few myogenesis-associated lncRNAs have been identified and our understanding of their regulatory mechanisms remains limited, particularly in goat. Here, we identified a novel lncRNA, TCONS_00006810 (named lncR-125b), from our previous lncRNA sequencing data on fetal (45, 60, and 105 days of gestation, three biological replicates for each point) and postnatal (3 days after birth, n = 3) goat skeletal muscle, and found that it is highly expressed in skeletal muscle and gradually upregulated during skeletal muscle satellite cell (SMSC) differentiation in goat. Notably, overexpression of lncR-125b accelerated the expression of myogenic differentiation 1 (MyoD 1) and myogenin (MyoG), and the formation of myotubes, and knockdown of lncR-125b showed opposite effects in SMSCs. Results of dual-luciferase assay and quantitative real-time polymerase chain reaction revealed that lncR-125b acts as a molecular sponge for miR-125b and that insulin-like growth factor 2 (IGF2), a critical regulator of skeletal myogenesis, is a direct target gene of miR-125b. Further analyses showed that lncR-125b negatively regulates miR-125b expression and positively regulates IGF2 expression in SMSCs. Mechanistically, lncR-125b promotes SMSC differentiation by functioning as a competing endogenous RNA (ceRNA) for miR-125b to control IGF2 expression. These findings identify lncR-125b as a novel noncoding regulator of muscle cell differentiation and skeletal muscle development in goat.

Novel lncRNA lncFAM200B: Molecular Characteristics and Effects of Genetic Variants on Promoter Activity and Cattle Body Measurement Traits

Skeletal muscle is one of the three major muscle types in an organism and has key roles in the motor system, metabolism, and homeostasis. RNA-Seq analysis showed that novel lncRNA, lncFAM200B, was differentially expressed in embryonic, neonatal, and adult cattle skeletal muscles. The main aim of this study was to investigate the molecular and expression characteristics of lncFAM200B along with its crucial genetic variations. Our results showed that bovine lncFAM200B was a 472 nucleotide (nt) non-coding RNA containing two exons. The transcription factor binding site prediction analysis found that lncFAM200B promoter region was enriched with SP1 transcription factor, which promotes the binding of myogenic regulatory factor MyoD and DNA sequence. The mRNA expression analysis showed that lncFAM200B was differentially expressed in embryonic, neonatal, adult bovine muscle tissues, and the lncFAM200B expression trend positively correlated with that of MyoG and Myf5 in myoblast proliferation and differential stages. To identify the promoter active region of lncFAM200B, we constructed promoter luciferase reporter gene vector pGL3-Basic plasmids containing lncFAM200B promoter sequences and transfected them into 293T, C2C12, and 3T3-L1 cells. Our results suggested that lncFAM200B promoter active region was from −403 to −139 (264 nt) of its transcription start site, covering 6 SP1 potential binding sites. Furthermore, we found a novel C-T variation, named as SNP2 (ERZ990081 in European Variation Archive) in the promoter active region, which was linked to the nearby SNP1 (rs456951291 in Ensembl database). The genotypes of SNP1 and combined genotypes of SNP1 and SNP2 were significantly associated with Jinnan cattle hip height. The luciferase activity analysis found that the SNP1-SNP2 haplotype CC had the highest luciferase activity, which was consistent with the association analysis result that the combined genotype CC-CC carriers had the highest hip height in Jinnan cattle. In conclusion, our data showed that lncFAM200B is a positive regulator of muscle development and that SNP1 and SNP2 could be used as genetic markers for marker-assisted selection (MAS) breeding of beef cattle.

lnc9141-a and -b Play a Different Role in Bovine Myoblast Proliferation, Apoptosis, and Differentiation

Previously, our transcriptome sequencing revealed that lnc9141 was differentially expressed in muscles of fetal bovine, calf, and adult bovine, which is considered to provide the basis for raising the muscle mass. In this study, we identified lnc9141 characters. lnc9141 has different transcription start sites and 3′ alternative splicing sites of exon 1, producing lnc9141-a and lnc9141-b transcripts that were highly expressed in the heart and lung. Moreover, neither lnc9141-a nor lnc9141-b had the ability to encode proteins. The functions of lnc9141-a and lnc9141-b were explored by cell cycle, 5-ethynyl-2'-deoxyuridine (EdU), and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The results showed that lnc9141-a or lnc9141-b overexpression decreased the number of myoblasts in the S phase and increased the proportion of cells in the G₀/G₁ phase. Furthermore, overexpressing lnc9141-a and lnc9141-b respectively downregulated the expression of Cyclin D1. However, lnc9141-a or lnc9141-b interference was found to increase the number of S-phase myoblasts, and upregulate Cyclin D1 and Cyclin E expression. Through Annexin V-FITC/propidium iodide (PI) double staining and the expression of apoptosis marker genes (Bax, Bcl2, and Caspase-3), it was found that lnc9141-b could regulate the expression of Bax gene. Meantime, high expression of lnc9141-b could decrease MyHC expression. In addition, the intergenic region between lnc9141 and IRX5 was 2.3 kb, with a head-to-head orientation. The study also revealed the core regions of the lnc9141 and IRX5 promoter. Our study demonstrated that both lnc9141-a and -b expression inhibited bovine myoblast proliferation. However, lnc9141-b regulated Bax and MyHC expression. The regulatory mechanism of lnc9141-a and lnc9141-b needs to be further explored.

Functions and Regulatory Mechanisms of lncRNAs in Skeletal Myogenesis, Muscle Disease and Meat Production

Myogenesis is a complex biological process, and understanding the regulatory network of skeletal myogenesis will contribute to the treatment of human muscle related diseases and improvement of agricultural animal meat production. Long noncoding RNAs (lncRNAs) serve as regulators in gene expression networks, and participate in various biological processes. Recent studies have identified functional lncRNAs involved in skeletal muscle development and disease. These lncRNAs regulate the proliferation, differentiation, and fusion of myoblasts through multiple mechanisms, such as chromatin modification, transcription regulation, and microRNA sponge activity. In this review, we presented the latest advances regarding the functions and regulatory activities of lncRNAs involved in muscle development, muscle disease, and meat production. Moreover, challenges and future perspectives related to the identification of functional lncRNAs were also discussed.