Profiling and Functional Analysis of Circular RNAs in Porcine Fast and Slow Muscles

Research Progress on the Regulating Factors of Muscle Fiber Heterogeneity in Livestock: A Review

The type of muscle fiber plays a crucial role in the growth, development, and dynamic plasticity of animals' skeletal muscle. Additionally, it is a primary determinant of the quality of both fresh and processed meat. Therefore, understanding the regulatory factors that contribute to muscle fibers' heterogeneity is of paramount importance. Recent advances in sequencing and omics technologies have enabled comprehensive cross-verification of research on the factors affecting the types of muscle fiber across multiple levels, including the genome, transcriptome, proteome, and metabolome. These advancements have facilitated deeper exploration into the related biological questions. This review focused on the impact of individual characteristics, feeding patterns, and genetic regulation on the proportion and interconversion of different muscle fibers. The findings indicated that individual characteristics and feeding patterns significantly influence the type of muscle fiber, which can effectively enhance the type and distribution of muscle fibers in livestock. Furthermore, non-coding RNA, genes and signaling pathways between complicated regulatory mechanisms and interactions have a certain degree of impact on muscle fibers' heterogeneity. This, in turn, changes muscle fiber profile in living animals through genetic selection or environmental factors, and has the potential to modulate the quality of fresh meat. Collectively, we briefly reviewed the structure of skeletal muscle tissue and then attempted to review the inevitable connection between the quality of fresh meat and the type of muscle fiber, with particular attention to potential events involved in regulating muscle fibers' heterogeneity.

A novel protein encoded by circKANSL1L regulates skeletal myogenesis via the Akt-FoxO3 signaling axis

Skeletal muscle is one the largest organs of the body and is involved in animal production and human health. Circular RNAs (circRNAs) have been implicated in skeletal myogenesis through largely unknown mechanisms. Herein, we report the phenotypic and metabolomic analysis of porcine longissimus dorsi muscles in Lantang and Landrace piglets, revealing a high-content of slow-oxidative fibers responsible for high-quality meat product in Lantang piglets. Using single-cell transcriptomics, we identified four myogenesis-related cell types, and the Akt-FoxO3 signaling axis was the most significantly enriched pathway in each subpopulation in the different pig breeds, as well as in fast-twitch glycolytic fibers. Using the multi-dimensional bioinformatic tools of circRNAome-seq and Ribo-seq, we identified a novel circRNA, circKANSL1L, with a protein-coding ability in porcine muscles, whose expression level correlated with myoblast proliferation and differentiation in vitro, as well as the transformation between distinct mature myofibers in vivo. The protein product of circKANSL1L could interact with Akt to decrease the phosphorylation level of FoxO3, which subsequently promoted FoxO3 transcriptional activity to regulate skeletal myogenesis. Our results established the existence of a protein encoded by circKANSL1L and demonstrated its potential functions in myogenesis.

Integrated analysis of circRNA, lncRNA, miRNA and mRNA to reveal the ceRNA regulatory network of postnatal skeletal muscle development in Ningxiang pig

Introduction: The development of skeletal muscle is regulated by regulatory factors of genes and non-coding RNAs (ncRNAs). Methods: The objective of this study was to understand the transformation of muscle fiber type in the longissimus dorsi muscle of male Ningxiang pigs at four different growth stages (30, 90, 150, and 210 days after birth, n = 3) by histological analysis and whole transcriptome sequencing. Additionally, the study investigated the expression patterns of various RNAs involved in muscle fiber transformation and constructed a regulatory network for competing endogenous RNA (ceRNA) that includes circular RNA (circRNA)/long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA). Results: Histomorphology analysis showed that the diameter of muscle fiber reached its maximum at 150 days after birth. The slow muscle fiber transformation showed a pattern of initial decrease followed by an increase. 29,963 circRNAs, 2,683 lncRNAs, 986 miRNAs and 22,411 mRNAs with expression level ≥0 were identified by whole transcriptome sequencing. Furthermore, 642 differentially expressed circRNAs (DEc), 505 differentially expressed lncRNAs (DEl), 316 differentially expressed miRNAs (DEmi) and 6,090 differentially expressed mRNAs (DEm) were identified by differential expression analysis. Functions of differentially expressed mRNA were identified by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). GO enrichment analysis indicates that 40 known genes and 6 new genes are associated with skeletal muscle development. Additionally, KEGG analysis shows that these genes regulate skeletal muscle development via MAPK, FoxO, Hedgehog, PI3K-Akt, Notch, VEGF and other signaling pathways. Through protein-protein interaction (PPI) and transcription factor prediction (TFP), the action mode of skeletal muscle-related genes was explored. PPI analysis showed that there were stable interactions among 19 proteins, meanwhile, TFP analysis predicted 22 transcription factors such as HMG20B, MYF6, MYOD1 and MYOG, and 12 of the 19 interacting proteins were transcription factors. The regulatory network of ceRNA related to skeletal muscle development was constructed based on the correlation of various RNA expression levels and the targeted binding characteristics with miRNA. The regulatory network included 31 DEms, 59 miRNAs, 667 circRNAs and 224 lncRNAs. conclusion: Overall, the study revealed the role of ceRNA regulatory network in the transformation of skeletal muscle fiber types in Ningxiang pigs, which contributes to the understanding of ceRNA regulatory network in Ningxiang pigs during the skeletal muscle development period.

Circular Intronic RNA circTTN Inhibits Host Gene Transcription and Myogenesis by Recruiting PURB Proteins to form Heterotypic Complexes

Muscle cell growth plays an important role in skeletal muscle development. Circular RNAs (circRNAs) have been proven to be involved in the regulation of skeletal muscle growth and development. In this study, we explored the effect of circTTN on myoblast growth and its possible molecular mechanism. Using C2C12 cells as a functional model, the authenticity of circTTN was confirmed by RNase R digestion and Sanger sequencing. Previous functional studies have showed that the overexpression of circTTN inhibits myoblast proliferation and differentiation. Mechanistically, circTTN recruits the PURB protein on the Titin (TTN) promoter to inhibit the expression of the TTN gene. Moreover, PURB inhibits myoblast proliferation and differentiation, which is consistent with circTTN function. In summary, our results indicate that circTTN inhibits the transcription and myogenesis of the host gene TTN by recruiting PURB proteins to form heterotypic complexes. This work may act as a reference for further research on the role of circRNA in skeletal muscle growth and development.

Melatonin promotes the development of the secondary hair follicles by regulating circMPP5

BACKGROUND

The quality and yield of cashmere fibre are closely related to the differentiation and development of secondary hair follicles in the skin of cashmere goats. The higher the density of secondary hair follicles, the higher the quality and yield of cashmere from the fleece. Development of secondary hair follicles commences in the embryonic stage of life and is completed 6 months after birth. Preliminary experimental results from our laboratory showed that melatonin (MT) treatment of goat kids after their birth could increase the density of secondary hair follicles and, thus, improve the subsequent yield and quality of cashmere. These changes in the secondary hair follicles resulted from increases in levels of antioxidant and expression of anti-apoptotic protein, and from a reduction in apoptosis. The present study was conducted to explore the molecular mechanism of MT-induced secondary hair follicle differentiation and development by using whole-genome analysis.

RESULTS

MT had no adverse effect on the growth performance of cashmere kids but significantly improved the character of the secondary hair follicles and the quality of cashmere, and this dominant effect continued to the second year. Melatonin promotes the proliferation of secondary hair follicle cells at an early age. The formation of secondary hair follicles in the MT group was earlier than that in the control group in the second year. The genome-wide data results involved KEGG analysis of 1044 DEmRNAs, 91 DElncRNAs, 1054 DEcircRNAs, and 61 DEmiRNAs which revealed that the mitogen-activated protein kinase (MAPK) signaling pathway is involved in the development of secondary hair follicles, with key genes (FGF2, FGF21, FGFR3, MAPK3 (ERK1)) being up-regulated and expressed. We also found that the circMPP5 could sponged miR-211 and regulate the expression of MAPK3.

CONCLUSIONS

We conclude that MT achieves its effects by regulating the MAPK pathway through the circMPP5 sponged the miR-211, regulating the expression of MAPK3, to induce the differentiation and proliferation of secondary hair follicle cells. In addition there is up-regulation of expression of the anti-apoptotic protein causing reduced apoptosis of hair follicle cells. Collectively, these events increase the numbers of secondary hair follicles, thus improving the production of cashmere from these goats.

Transcriptome-Wide Study Revealed That N6-Methyladenosine Participates in Regulation Meat Production in Goats

In mammals, skeletal muscle development is a complex biological process regulated by many factors. N6-methyladenosine (m6A) RNA modification plays an important role in many biological processes. However, the regulation of m6A on skeletal muscle growth and development in adult goats remains unclear. In this study, Duan goats (DA) and Nubia goats (NBY), both female and 12 months old, were selected as the research objects, and m6A-Seq and RNA-Seq were mainly used to detect the difference of m6A modification and gene expression during the development of the longissimus dorsi (LD) muscle in the two breeds. The results showed that compared with DA, the meat production performance of NBY was better than that of DA, and the modification level of m6A was higher than that of DA in LD. The m6A-Seq of LD indicated m6A peaks were mainly enriched in the coding sequence (CDS) and stop codon. A total of 161 differentially methylated genes (DMGs) and 1294 differentially expressed genes (DEGs) were identified in two breeds. GO and KEGG analysis showed that DMGs were closely related to cellular metabolism, and most of DMGs were enriched in pathways related to energy metabolism, muscle growth and development, mainly MAPK signaling pathway, Wnt signaling pathway and CGMP-PKG signaling pathway. The DEGs were significantly enriched in actin binding, calcium ion binding, angiogenesis, and other biological processes, and most of them were enriched in PI3K-Akt and CGMP-PKG signaling pathways. Combined analysis of m6A-Seq and RNA-Seq data revealed a negative correlation between differentially methylated m6A levels and mRNA abundance, and mRNA expression of the gene with m6A peak near 3′UTR will decrease. In addition, 11 DMGs regulating cell differentiation, muscle growth and development were identified. This study displayed the m6A profiles and distribution patterns in the goat transcriptome, determined the potential role of m6A modification in muscle growth and provided a new reference for the further study of goat skeletal muscle development.

Recent Progress on Circular RNAs in the Development of Skeletal Muscle and Adipose Tissues of Farm Animals

Circular RNAs (circRNAs) are a highly conserved and specifically expressed novel class of covalently closed non-coding RNAs. CircRNAs can function as miRNA sponges, protein scaffolds, and regulatory factors, and play various roles in development and other biological processes in mammals. With the rapid development of high-throughput sequencing technology, thousands of circRNAs have been discovered in farm animals; some reportedly play vital roles in skeletal muscle and adipose development. These are critical factors affecting meat yield and quality. In this review, we have highlighted the recent advances in circRNA-related studies of skeletal muscle and adipose in farm animals. We have also described the biogenesis, properties, and biological functions of circRNAs. Furthermore, we have comprehensively summarized the functions and regulatory mechanisms of circRNAs in skeletal muscle and adipose development in farm animals and their effects on economic traits such as meat yield and quality. Finally, we propose that circRNAs are putative novel targets to improve meat yield and quality traits during animal breeding.

Comparative Proteomic Analysis of Glycolytic and Oxidative Muscle in Pigs

The quality of meat is highly correlated with muscle fiber type. However, the mechanisms via which proteins regulate muscle fiber types in pigs are not entirely understood. In the current study, we have performed proteomic profiling of fast/glycolytic biceps femoris (BF) and slow/oxidative soleus (SOL) muscles and identified several candidate differential proteins among these. We performed proteomic analyses based on tandem mass tags (TMTs) and identified a total of 26,228 peptides corresponding to 2667 proteins among the BF and SOL muscle samples. Among these, we found 204 differentially expressed proteins (DEPs) between BF and SOL muscle, with 56 up-regulated and 148 down-regulated DEPs in SOL muscle samples. KEGG and GO enrichment analyses of the DEPs revealed that the DEPs are involved in some GO terms (e.g., actin cytoskeleton, myosin complex, and cytoskeletal parts) and signaling pathways (PI3K-Akt and NF-kappa B signaling pathways) that influence muscle fiber type. A regulatory network of protein-protein interaction (PPI) between these DEPs that regulates muscle fiber types was constructed, which demonstrates how three down-regulated DEPs, including PFKM, GAPDH, and PKM, interact with other proteins to potentially control the glycolytic process. This study offers a new understanding of the molecular mechanisms in glycolytic and oxidative muscles as well as a novel approach for enhancing meat quality by transforming the type of muscle fibers in pigs.

Phosphoproteomic analysis identifies differentially expressed phosphorylation sites that affect muscle fiber type in pigs

Skeletal muscle of livestock is composed of both fast- and slow-twitch muscle fibers, which are key factors in their meat quality. However, the role of protein phosphorylation in muscle fiber type is not completely understood. Here, a fast-twitch (biceps femoris, BF) and slow-twitch (soleus, SOL) muscle tissue sample was collected from three male offspring of Duroc and Meishan pigs. We demonstrate that the meat quality of SOL muscle is significantly better than that of BF muscle. We further used phosphoproteomic profiling of BF and SOL muscles to identify differences between these muscle types. A total of 2,327 phosphorylation sites from 770 phosphoproteins were identified. Among these sites, 287 differentially expressed phosphorylation sites (DEPSs) were identified between BF and SOL. GO and KEGG enrichment analysis of proteins containing DEPSs showed that these phosphorylated proteins were enriched in the glycolytic process GO term and the AMPK signaling pathway. A protein-protein interaction (PPI) analysis reveals that these phosphorylated proteins interact with each other to regulate the transformation of muscle fiber type. These analyses reveal that protein phosphorylation modifications are involved in porcine skeletal muscle fiber type transformation. This study provides new insights into the molecular mechanisms by which protein phosphorylation regulates muscle fiber type transformation and meat quality in pigs.

Whole transcriptome analyses and comparison reveal the metabolic differences between oxidative and glycolytic skeletal muscles of yak

Mammalian skeletal muscle is composed of various muscle fibers that exhibit different physiological and metabolic features. Muscle fiber type composition has significant influences on the meat quality of livestock. In this study, we comprehensively analyzed the whole transcriptome profiles of the oxidative muscle biceps femoris (BF) and the glycolytic muscle obliquus externus abdominis (OEA) of yak. A total of 1436 mRNAs, 1172 lncRNAs, and 218 circRNAs were differentially expressed in the oxidative muscles compared with the glycolytic muscles. KEGG annotation showed that differentially expressed mRNAs regulated by lncRNA and circRNA were mainly involved in PPAR signaling pathway, citrate cycle (TCA cycle), and PI3K-Akt signaling pathway, which reflect the different metabolic properties between oxidative and glycolytic muscles. In addition, regulatory networks associated with muscle fiber type conversion and mitochondria energy metabolism in muscles were constructed. Our study provides new evidence for a better understanding of the molecular mechanisms underlying skeletal muscle fiber determination and meat quality traits of yak.

Identification and characterization of circRNAs related to meat quality during embryonic development of the longissimus dorsi muscle in two pig breeds

Meat quality, an important economic trait, is regulated by many factors, especially by genetic factors, including coding genes, miRNAs, and lncRNAs. Recent studies have elucidated that circRNAs also play a key role in muscle development and lipid deposition. However, the functions and regulatory mechanisms of circRNAs in meat quality remain mostly unknown. The circRNA expression profiles between Huainan pigs (Chinese indigenous pigs, fat-type, Huainan HN) and Large White pigs (Western commercial pigs, lean-type, LW) in the longissimus dorsi (LD) muscle at 38, 58, and 78 days post conception (dpc) were compared by sequencing. In total, 39,887 circRNAs were identified in 18 samples, and 60, 78, and 86 differentially expressed circRNAs (DECs) were found at the three stages mentioned above between these two breeds. The parent genes of DECs were enriched in myogenesis, proliferation, adipogenesis and muscle fiber-type transition. The circRNA-miRNA interaction networks included 38 DECs and 47 miRNAs, and these miRNAs were involved in muscle development and lipid metabolism. Two shared DECs (circ_0030593 and circ_0032760) of these three stages were selected, their head-to-tail junction sites were validated by Sanger sequencing, and RT‒qPCR results suggested that these two DECs might be involved in intramuscular fat deposition. These findings provide a basis for understanding the role of circRNAs in meat quality.

A Combined Differential Proteome and Transcriptome Profiling of Fast- and Slow-Twitch Skeletal Muscle in Pigs

Skeletal muscle fiber types can contribute in part to affecting pork quality parameters. Biceps femoris (Bf) (fast muscle or white muscle) and Soleus (Sol) (slow muscle or red muscle) are two typical skeletal muscles characterized by obvious muscle fiber type differences in pigs. However, the critical proteins and potential regulatory mechanisms regulating porcine skeletal muscle fibers have yet to be clearly defined. In this study, the isobaric Tag for Relative and Absolute Quantification (iTRAQ)-based proteome was used to identify the key proteins affecting the skeletal muscle fiber types with Bf and Sol, by integrating the previous transcriptome data, while function enrichment analysis and a protein–protein interaction (PPI) network were utilized to explore the potential regulatory mechanisms of skeletal muscle fibers. A total of 126 differentially abundant proteins (DAPs) between the Bf and Sol were identified, and 12 genes were found to be overlapping between differentially expressed genes (DEGs) and DAPs, which are the critical proteins regulating the formation of skeletal muscle fibers. Functional enrichment and PPI analysis showed that the DAPs were mainly involved in the skeletal-muscle-associated structural proteins, mitochondria and energy metabolism, tricarboxylic acid cycle, fatty acid metabolism, and kinase activity, suggesting that PPI networks including DAPs are the main regulatory network affecting muscle fiber formation. Overall, these data provide valuable information for understanding the molecular mechanism underlying the formation and conversion of muscle fiber types, and provide potential markers for the evaluation of meat quality.

Circular RNA Profiling Identifies Novel circPPARA that Promotes Intramuscular Fat Deposition in Pigs

Intramuscular fat (IMF) content plays an important role in pork quality. Circular RNAs (circRNAs) implicate various biological processes; however, the regulatory mechanisms and functions of circRNAs in porcine IMF remains elusive. Hence, the study assessed the circRNA expression profiling in the longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to unravel their regulatory functions in improving meat quality. The RNA sequencing analysis identified 29,732 circRNAs from six sampled pigs, most of which were exon-derived. In the muscle, 336 were differentially expressed (DE) between the H and L IMF groups; 196 circRNAs were upregulated, and 140 were downregulated. Subsequent qRT-PCR validation of 10 DE circRNAs revealed expression patterns consistent with the RNA-seq data. Gene ontology and KEGG enrichment analysis revealed that most significantly enriched DE circRNAs' host genes were linked to lipid metabolism and adipogenesis processes. The circRNA-miRNA regulatory network analysis found several circRNAs targeting miRNAs associated with adipogenesis. Finally, a novel circRNA, circPPARA, was identified with the expression positively correlated with the IMF content. Detailed analysis revealed that circPPARA was formed via head-to-tail splicing and was more stable than the linear PPARA, predominantly located in the cytoplasm. Functional studies using overexpression and siRNA constructs demonstrated that circPPARA promotes differentiation and hinders the proliferation of porcine intramuscular preadipocytes. Moreover, the dual-luciferase assay revealed that circPPARA adsorbed miR-429 and miR-200b, thereby promoting intramuscular adipogenesis in pigs. Our results identified a candidate circRNA, circPPARA, that affects porcine IMF content. The study provides knowledge of the regulatory functions of circRNAs in intramuscular adipogenesis and abundant resource for future research on circRNAs in pigs.

circTAF8 Regulates Myoblast Development and Associated Carcass Traits in Chicken

Recent studies have shown that circular RNAs (circRNAs) play important roles in skeletal muscle development. CircRNA biogenesis is dependent on the genetic context. Single-nucleotide polymorphisms in the introns flanking circRNAs may be intermediate-inducible factors between circRNA expression and phenotypic traits. Our previous study showed that circTAF8 is an abundantly and differentially expressed circRNA in leg muscle during chicken embryonic development. Here, we aimed to investigate circTAF8 function in muscle development and the association of the SNPs in the circTAF8 flanking introns with carcass traits. In this study, we observed that overexpression of circTAF8 could promote the proliferation of chicken primary myoblasts and inhibit their differentiation. In addition, the SNPs in the introns flanking the circTAF8 locus and those associated with chicken carcass traits were analyzed in 335 partridge chickens. A total of eight SNPs were found associated with carcass traits such as leg muscle weight, live weight, and half and full-bore weight. The association analysis results of haplotype combinations were consistent with the association analysis of a single SNP. These results suggest that circTAF8 plays a regulatory role in muscle development. These identified SNPs were found correlated with traits to muscle development and carcass muscle weight in chickens.

Dynamic Expression and Regulatory Network of Circular RNA for Abdominal Preadipocytes Differentiation in Chicken (Gallus gallus)

Circular RNA (circRNA), as a novel endogenous biomolecule, has been emergingly demonstrated to play crucial roles in mammalian lipid metabolism and obesity. However, little is known about their genome-wide identification, expression profile, and function in chicken adipogenesis. In present study, the adipogenic differentiation of chicken abdominal preadipocyte was successfully induced, and the regulatory functional circRNAs in chicken adipogenesis were identified from abdominal adipocytes at different differentiation stages using Ribo-Zero RNA-seq. A total of 1,068 circRNA candidates were identified and mostly derived from exons. Of these, 111 differentially expressed circRNAs (DE-circRNAs) were detected, characterized by stage-specific expression, and enriched in several lipid-related pathways, such as Hippo signaling pathway, mTOR signaling pathway. Through weighted gene co-expression network analyses (WGCNA) and K-means clustering analyses, two DE-circRNAs, Z:35565770|35568133 and Z:54674624|54755962, were identified as candidate regulatory circRNAs in chicken adipogenic differentiation. Z:35565770|35568133 might compete splicing with its parental gene, ABHD17B, owing to its strictly negative co-expression. We also constructed competing endogenous RNA (ceRNA) network based on DE-circRNA, DE-miRNA, DE-mRNAs, revealing that Z:54674624|54755962 might function as a ceRNA to regulate chicken adipogenic differentiation through the gga-miR-1635-AHR2/IRF1/MGAT3/ABCA1/AADAC and/or the novel_miR_232-STAT5A axis. Translation activity analysis showed that Z:35565770|35568133 and Z:54674624|54755962 have no protein-coding potential. These findings provide valuable evidence for a better understanding of the specific functions and molecular mechanisms of circRNAs underlying avian adipogenesis.

miR-152 targets pyruvate kinase to regulate the glycolytic activity of pig skeletal muscles and affects pork quality

As a type of non-coding RNA, microRNAs are widely involved in the biological processes of animals. In the present study, the expression of miR-152 in glycolytic muscle fibers (Longissimus thoracis, LT) was lower than that of oxidative muscle fibers (Psoas major, PM). Using dual luciferase assay, miR-152 was shown to target muscle pyruvate kinase (PKM) to perform biological functions. Moreover, overexpression of miR-152 in primary porcine cells inhibited PKM gene expression and reduced lactic acid production in cells, whereas inhibition of miR-152 expression promoted PKM gene expression and increased lactic acid production. Correlation analysis showed that the expression of miR-152 was significantly positively correlated with the ultimate pH of LT after slaughter, while the expression of the PKM gene was significantly negatively correlated with the final pH of LT. In vivo and in vitro experiments discussed herein suggest that miR-152 may affect muscle pH by targeting the expression of the PKM gene. Our findings enrich the understanding of the genetic regulatory network that influences pork quality.

Expression Profiling and Functional Analysis of Circular RNAs in Inner Mongolian Cashmere Goat Hair Follicles

Background

Inner Mongolian cashmere goats have hair of excellent quality and high economic value, and the skin hair follicle traits of cashmere goats have a direct and important effect on cashmere yield and quality. Circular RNA has been studied in a variety of tissues and cells.

Result

In this study, high-throughput sequencing was used to obtain the expression profiles of circular RNA (circRNA) in the hair follicles of Inner Mongolian cashmere goats at different embryonic stages (45, 55, 65, and 75 days). A total of 21,784 circRNAs were identified. At the same time, the differentially expressed circRNA in the six comparison groups formed in the four stages were: d75vsd45, 59 upregulated and 33 downregulated DE circRNAs; d75vsd55, 61 upregulated and 102 downregulated DE circRNAs; d75vsd65, 32 upregulated and 33 downregulated DE circRNAs; d65vsd55, 67 upregulated and 169 downregulated DE circRNAs; d65vsd45, 96 upregulated and 63 downregulated DE circRNAs; and d55vsd45, 76 upregulated and 42 downregulated DE circRNAs. Six DE circRNA were randomly selected to verify the reliability of the sequencing results by quantitative RT-PCR. Subsequently, the circRNA corresponding host genes were analyzed by the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The results showed that the biological processes related to hair follicle growth and development enriched by GO mainly included hair follicle morphogenesis and cell development, and the signaling pathways related to hair follicle development included the Notch signaling pathway and NF-κB signaling pathway. We combined the DE circRNA of d75vsd45 with miRNA and mRNA databases (unpublished) to construct the regulatory network of circRNA–miRNA–mRNA, and formed a total of 102 pairs of circRNA–miRNA and 126 pairs of miRNA–mRNA interactions. The binding relationship of circRNA3236–chi-miR-27b-3p and circRNA3236–chi-miR-16b-3p was further verified by dual-luciferase reporter assays, and the results showed that circRNA3236 and chi-miR-27b-3p, and circRNA3236 and chi-miR-16b-3p have a targeted binding relationship.

Conclusion

To summarize, we established the expression profiling of circRNA in the fetal skin hair follicles of cashmere goats, and found that the host gene of circRNA may be involved in the development of hair follicles of cashmere goats. The regulatory network of circRNA–miRNA–mRNA was constructed and preliminarily verified using DE circRNAs.

Dietary Supplementation With Chlorogenic Acid Derived From Lonicera macranthoides Hand-Mazz Improves Meat Quality and Muscle Fiber Characteristics of Finishing Pigs via Enhancement of Antioxidant Capacity

Chlorogenic acid (CGA), one of the most abundant polyphenol compounds in nature, is regarded as a potential feed additive to promote animal health and enhance the meat products’ quality via its various biological properties. The current study aims: (1) to determine whether dietary CGA supplementation improves meat quality and muscle fiber characteristics, and (2) to ascertain whether the corresponding improvement is associated with enhancing the antioxidant capacity of the finishing pigs. Thirty-two (Large × White × Landrace) finishing pigs with an average initial body weight of 71.89 ± 0.92 kg were allotted to 4 groups, and each was fed diets supplemented with 0, 0.02, 0.04, or 0.08% (weight/weight) of CGA. The meat quality traits, muscle fiber characteristics, and the serum and muscle antioxidant capacity were assessed. Results suggested that, compared with the control group, dietary CGA supplementation at a level of 0.04% significantly decreased the b^∗ value and distinctly increased the inosinic acid content of longissimus dorsi (LD) and biceps femoris (BF) muscles (P < 0.01). Moreover, dietary supplementation with 0.04% of CGA markedly improved the amino acid composition of LD and BF muscles, as well as augmented the mRNA abundance of Nrf-2, GPX-1, MyoD, MyoG, and oxidative muscle fiber (I and IIa) in LD muscle (P < 0.05). This result indicates that a diet supplemented with 0.04% of CGA promotes myogenesis and induces a transformation toward more oxidative muscle fibers in LD muscle, subsequently improving meat quality. Besides, dietary supplementation with 0.02% and 0.04% of CGA notably enhanced the serum GSH-PX level (P < 0.01). Considering all these effects are closely related to the alteration of antioxidant activities of the finishing pigs, the underlying metabolism is likely connected to the boosting of their antioxidant capacity induced by dietary CGA.

The Roles of CircRNAs in Regulating Muscle Development of Livestock Animals

The muscle growth and development of livestock animals is a complex, multistage process, which is regulated by many factors, especially the genes related to muscle development. In recent years, it has been reported frequently that circular RNAs (circRNAs) are involved widely in cell proliferation, cell differentiation, and body development (including muscle development). However, the research on circRNAs in muscle growth and development of livestock animals is still in its infancy. In this paper, we briefly introduce the discovery, classification, biogenesis, biological function, and degradation of circRNAs and focus on the molecular mechanism and mode of action of circRNAs as competitive endogenous RNAs in the muscle development of livestock and poultry. In addition, we also discuss the regulatory mechanism of circRNAs on muscle development in livestock in terms of transcription, translation, and mRNAs. The purpose of this article is to discuss the multiple regulatory roles of circRNAs in the process of muscle development in livestock, to provide new ideas for the development of a new co-expression regulation network, and to lay a foundation for enriching livestock breeding and improving livestock economic traits.

Functional Role of circRNAs in the Regulation of Fetal Development, Muscle Development, and Lactation in Livestock

circRNAs are a class of endogenous noncoding RNA molecules with closed loop structures. They are mainly responsible for regulating gene expression in eukaryotic cells. With the emergence of high-throughput RNA sequencing (RNA-Seq) and new types of bioinformatics tools, thousands of circRNAs have been discovered, making circRNA one of the research hotspots. Recent studies have shown that circRNAs play an important regulatory role in the growth, reproduction, and formation of livestock products. They can not only regulate mammalian fetal growth and development but also have important regulatory effects on livestock muscle development and lactation. In this review, we briefly introduce the putative biogenic pathways and regulatory functions of circRNA and highlight our understanding of circRNA and its latest advances in fetal development, muscle development, and lactation biogenesis as well as expression in livestock. This review will provide a theoretical basis for the research and development of related industries.

Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.