Comprehensive Analysis of lncRNAs and circRNAs Reveals the Metabolic Specialization in Oxidative and Glycolytic Skeletal Muscles

Unveiling the role of circRBBP7 in myoblast proliferation and differentiation: A novel regulator of muscle development

Muscle development is a multistep process regulated by diverse gene networks, and circRNAs are considered novel regulators mediating myogenesis. Here, we systematically analyzed the role and underlying regulatory mechanisms of circRBBP7 in myoblast proliferation and differentiation. Results showed that circRBBP7 has a typical circular structure and encodes a 13 -kDa protein. By performing circRBBP7 overexpression and RNA interference, we found that the function of circRBBP7 was positively correlated with the proliferation and differentiation of myoblasts. Using RNA sequencing, we identified 1633 and 532 differentially expressed genes (DEGs) during myoblast proliferation or differentiation, respectively. The DEGs were found mainly enriched in cell cycle- and skeletal muscle development-related pathways, such as the MDM2/p53 and PI3K-Akt signaling pathways. Further co-IP and IF co-localization analysis revealed that VEGFR-1 is a target of circRBBP7 in myoblasts. qRT-PCR and WB analysis further confirmed the positive correlation between VEGFR-1 and circRBBP7. Moreover, we found that in vivo transfection of circRBBP7 into injured muscle tissues significantly promoted the regeneration and repair of myofibers in mice. Therefore, we speculate that circRBBP7 may affect the activity of MDM2 by targeting VEGFR-1, altering the expression of muscle development-related genes by mediating p53 degradation, and ultimately promoting myoblast development and muscle regeneration. This study provides essential evidence that circRBBP7 can serve as a potential target for myogenesis regulation and a reference for the application of circRBBP7 in cattle genetic breeding and muscle injury treatment.

circUSP13 facilitates the fast-to-slow myofiber shift via the MAPK/ERK signaling pathway in goat skeletal muscles

Understanding how skeletal muscle fiber proportions are regulated is essential for understanding muscle function and improving the quality of mutton. While circular RNA (circRNA) has a critical function in myofiber type transformation, the specific mechanisms are not yet fully understood. Prior evidence indicates that circular ubiquitin-specific peptidase 13 (circUSP13) can promote myoblast differentiation by acting as a ceRNA, but its potential role in myofiber switching is still unknown. Herein, we found that circUSP13 enhanced slow myosin heavy chain (MyHC-slow) and suppressed MyHC-fast expression in goat primary myoblasts (GPMs). Meanwhile, circUSP13 evidently enhanced the remodeling of the mitochondrial network while inhibiting the autophagy of GPMs. We obtained fast-dominated myofibers, via treatment with rotenone, and further demonstrated the positive role of circUSP13 in the fast-to-slow transition. Mechanistically, activation of the mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) pathway significantly impaired the slow-to-fast shift in fully differentiated myotubes, which was restored by circUSP13 or IGF1 overexpression. In conclusion, circUSP13 promoted the fast-to-slow myofiber type transition through MAPK/ERK signaling in goat skeletal muscle. These findings provide novel insights into the role of circUSP13 in myofiber type transition and contribute to a better understanding of the genetic mechanisms underlying meat quality.

Circular RNAs Mediate the Effects of Dietary Tryptophan on the Transformation of Muscle Fiber Types in Weaned Piglets

The nutritional composition of the diet significantly impacts the overall growth and development of weaned piglets. The current study aimed to explore the effects and underlying mechanisms of dietary tryptophan consumption on muscle fiber type transformation during the weaning period. Thirty weaned piglets with an average body weight of 6.12 ± 0.16 kg were randomly divided into control (CON, 0.14% Trp diet) and high Trp (HT, 0.35% Trp) groups and maintained on the respective diet for 28 days. The HT group of weaned piglets exhibited highly significant improvements in growth performance and an increased proportion of fast muscle fibers. Transcriptome sequencing revealed the potential contribution of differentially expressed circular RNAs toward the transformation of myofiber types in piglets and toward the regulation of expression of related genes by targeting the microRNAs, miR-34c and miR-182, to further regulate myofiber transformation. In addition, 145 DE circRNAs were identified as potentially protein-encoding, with the encoded proteins associated with a myofiber type transformation. In conclusion, the current study greatly advances and refines our current understanding of the regulatory networks associated with piglet muscle development and myofiber type transformation and also contributes to the optimization of piglet diet formulation.

Tumor-derived lncRNAs: Behind-the-scenes mediators that modulate the immune system and play a role in cancer pathogenesis

Having been involved in complex cellular regulatory networks and cell-to-cell communications, non-coding RNAs (lncRNAs) have become functional carriers that transmit information between cells and tissues, modulate tumor microenvironments, encourage angiogenesis and invasion, and make tumor cells more resistant to drugs. Immune cells' exosomal lncRNAs may be introduced into tumor cells to influence the tumor's course and the treatment's effectiveness. Research has focused on determining if non-coding RNAs affect many target genes to mediate regulating recipient cells. The tumor microenvironment's immune and cancer cells are influenced by lncRNAs, which may impact a treatment's efficacy. The lncRNA-mediated interaction between cancer cells and immune cells invading the tumor microenvironment has been the subject of numerous recent studies. On the other hand, tumor-derived lncRNAs' control over the immune system has not gotten much attention and is still a relatively new area of study. Tumor-derived lncRNAs are recognized to contribute to tumor immunity, while the exact mechanism is unclear.

Advances in the discovery of genetic elements underlying longissimus dorsi muscle growth and development in the pig

As a major source of protein in human diets, pig meat plays a crucial role in ensuring global food security. Key determinants of meat production refer to the chemical and physical compositions or characteristics of muscle fibers, such as the number, hypertrophy potential, fiber-type conversion and intramuscular fat deposition. However, the growth and formation of muscle fibers comprises a complex process under spatio-temporal regulation, that is, the intermingled and concomitant proliferation, differentiation, migration and fusion of myoblasts. Recently, with the fast and continuous development of next-generation sequencing technology, the integration of quantitative trait loci mapping with genome-wide association studies (GWAS) has greatly helped animal geneticists to discover and explore thousands of functional or causal genetic elements underlying muscle growth and development. However, owing to the underlying complex molecular mechanisms, challenges to in-depth understanding and utilization remain, and the cost of large-scale sequencing, which requires integrated analyses of high-throughput omics data, is high. In this review, we mainly elaborate on research advances in integrative analyses (e.g. GWAS, omics) for identifying functional genes or genomic elements for longissimus dorsi muscle growth and development for different pig breeds, describing several successful transcriptome analyses and functional genomics cases, in an attempt to provide some perspective on the future functional annotation of genetic elements for muscle growth and development in pigs.

Identification and Functional Prediction of Circular RNAs Related to Growth Traits and Skeletal Muscle Development in Duroc pigs

Porcine skeletal muscle is a highly heterogeneous tissue type, and the Longissimus Dorsi muscle (LDM), as the most economical and physiologically metabolized skeletal muscle in pigs, has always been the focus of research and improvement in pig molecular breeding. Circular RNA, as an important new member of regulatory non-coding RNA after microRNA, has become a frontier hot spot in life science research. This study aims to explore candidate circRNAs related to growth, meat quality, and skeletal muscle development among Duroc pigs with different average daily gain (ADG). Eight pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG), followed by RNA-Seq analysis to identify circRNAs. The results showed that backfat at 6-7 rib (BF) and Intramuscular fat (IMF) content in the H group was significantly lower than L group, but ribeye area (REA) in the H group was higher than in the L group. In RNA-seq, 296 Differentially expressed (DE) circRNAs (157 upregulated and 139 downregulated) were identified and exons flanking long introns are easier to circularize to produce circRNAs. Most of the DE circRNAs were enriched in Quantitative trait locus (QTL) regions related to meat quality and growth traits. In addition, a gene can produce one or more circRNA transcripts. It was also found that the source genes of DE circRNAs were enriched in MAPK, FoXO, mTOR, PI3K-Akt, and Wnt signaling pathways. The results showed that different ADG, carcass, and meat quality traits among half-sibling Duroc pigs with the same diet may be due to the DE circRNAs related to skeletal muscle growth and development.

A Comprehensive Analysis of the Myocardial Transcriptome in ZBED6-Knockout Bama Xiang Pigs

The ZBED6 gene is a transcription factor that regulates the expression of IGF2 and affects muscle growth and development. However, its effect on the growth and development of the heart is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) can regulate genes at the epigenetic, transcriptional, and posttranscriptional levels and play an important role in the development of eukaryotes. To investigate the function of ZBED6 in the cardiac development of pigs, we constructed the expression profiles of mRNAs and lncRNAs in myocardial tissue obtained from Bama Xiang pigs in the ZBED6 knockout group (ZBED6-KO) and the wild-type group (ZBED6-WT). A total of 248 differentially expressed genes (DEGs) and 209 differentially expressed lncRNAs (DELs) were detected, and 105 potential cis target genes of DELs were identified. The functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) databases revealed two GO items related to muscle development by the cis target genes of DELs. Moreover, IGF2 was the direct target gene of ZBED6 by ChIP-PCR experiment. Our results explored the mechanism and expression profile of mRNAs and lncRNAs of ZBED6 gene knockout on myocardium tissue development, mining the key candidate genes in that process like IGF2.

Gut Microbiota Composition and Diversity in Different Commercial Swine Breeds in Early and Finishing Growth Stages

Simple Summary

The gut microbiota are involved in the metabolism of nutrients, and the growth and development of pig is strongly influenced by the gut microbiome. To maintain the integrity of the intestinal barrier and promote the digestion and absorption of nutrients and other physiological activities, it is beneficial if the host has a stable gut microbial community. The composition of the gut microbiota is influenced by many factors, such as genetic and environmental factors, and it changes with age. Throughout pig growth and development, the porcine gut microbiota constantly changes in composition. This study investigated the regulation of growth and development of body weight and body size index. We further examined changes in gut microbiota during early and finishing growth stages in Duroc, Landrace and Yorkshire pigs. Results showed that the microbiota of Landrace and Yorkshire pigs were more similar compared with Duroc pigs. There were significant differences in gut microbiota in the early and late growth stages. This study underlines the longitudinal variation in breed and lateral variation in age in gut microbiota.

Abstract

The gut microbiota affects the metabolism, health and growth rate of pigs. Understanding the characteristics of gut microbiota of different pig breeds at each growth stage will enable the design of individualized feeding strategies. The present study aimed to compare the growth curves and development patterns of pigs of three different breeds (Duroc, Landrace and Yorkshire) using the mathematical models Gompertz, Logistic, Von Bertalanffy and Richards. For Duroc pigs, the Gompertz model showed the highest prediction accuracy (R² = 0.9974). In contrast, the best models for Landrace and Yorkshire pigs were Richards (R² = 0.9986) and Von Bertalanffy (R² = 0.9977), respectively. Path analysis showed that body length (path coefficient  =  0.507) and chest circumference (path coefficient  =  0.532) contributed more significantly to the body weight of pigs at the early growth stage, while hip circumference (path coefficient  =  0.312) had a greater influence on pig body weight in the late growth stage. Moreover, the composition of the gut microbiota of pigs at two growth stages (60 kg of body weight in the early growth stage and 120 kg in the finishing stage) was studied using 16S rRNA sequencing technology. Variations in gut microbiota composition of pigs at different growth stages were observed. KEGG pathway enrichment analysis of annotated metagenomes revealed that protein synthesis and amino acid metabolism pathways were significantly enriched in pigs at the early growth stage, which may be related to nutritional requirements of pigs during this stage. This study confirmed longitudinal variation in the gut microbiota of pigs pertaining to age as well as lateral variation related to pig breed. The present findings expand the current understanding of the variations in swine gut microbiota during production stages.

Comprehensive Analysis of Long Noncoding RNA Modified by m6A Methylation in Oxidative and Glycolytic Skeletal Muscles

N⁶-methyladenosine (m⁶A) is the most common modification in eukaryotic RNAs. Accumulating evidence shows m⁶A methylation plays vital roles in various biological processes, including muscle and fat differentiation. However, there is a lack of research on lncRNAs’ m⁶A modification in regulating pig muscle-fiber-type conversion. In this study, we identified novel and differentially expressed lncRNAs in oxidative and glycolytic skeletal muscles through RNA-seq, and further reported the m⁶A-methylation patterns of lncRNAs via MeRIP-seq. We found that most lncRNAs have one m⁶A peak, and the m⁶A peaks were preferentially enriched in the last exon of the lncRNAs. Interestingly, we found that lncRNAs’ m⁶A levels were positively correlated with their expression homeostasis and levels. Furthermore, we performed conjoint analysis of MeRIP-seq and RNA-seq data and obtained 305 differentially expressed and differentially m⁶A-modified lncRNAs (dme-lncRNAs). Through QTL enrichment analysis of dme-lncRNAs and PPI analysis for their cis-genes, we finally identified seven key m⁶A-modified lncRNAs that may play a potential role in muscle-fiber-type conversion. Notably, inhibition of one of the key lncRNAs, MSTRG.14200.1, delayed satellite cell differentiation and stimulated fast-to-slow muscle-fiber conversion. Our study comprehensively analyzed m⁶A modifications on lncRNAs in oxidative and glycolytic skeletal muscles and provided new targets for the study of pig muscle-fiber-type conversion.

Characterization of the circRNA–miRNA–mRNA Network to Reveal the Potential Functional ceRNAs Associated With Dynamic Changes in the Meat Quality of the Longissimus Thoracis Muscle in Tibetan Sheep at Different Growth Stages

Circular RNAs (circRNAs) have a regulatory role in animal skeletal muscle development. In this study, RNA sequencing was performed to reveal the temporal regularity of circRNA expression and the effect of the circRNA–miRNA–mRNA ceRNA regulatory network on the meat quality of longissimus thoracis (LT) muscle in Tibetan sheep at different growth stages (4 months old, 4 m; 1.5 years old, 1.5 y; 3.5 years old, 3.5 y; 6 years old, 6 y). There were differences in the carcass performance and meat quality of Tibetan sheep at different ages. Especially, the meat tenderness significantly decreased (p < 0.05) with the increase of age. GO functional enrichment indicated that the source genes of the DE circRNAs were mainly involved in the protein binding, and myofibril and organelle assembly. Moreover, there was a significant KEGG enrichment in the adenosine 5′-monophosphate (AMP)-activated protein kinase (AMPK) signaling pathway, as well as the calcium signaling pathway, regulating the pluripotency of the stem cells. The circRNA–miRNA–mRNA ceRNA interaction network analysis indicated that circRNAs such as circ_000631, circ_000281, and circ_003400 combined with miR-29-3p and miR-185-5p regulate the expression of LEP, SCD, and FASN related to the transformation of muscle fiber types in the AMPK signaling pathway. The oxidized muscle fibers were transformed into the glycolytic muscle fibers with the increase of age, the content of intramuscular fat (IMF) was lowered, and the diameter of the muscle fiber was larger in the glycolytic muscle fibers, ultimately increasing the meat tenderness. The study revealed the role of the circRNAs in the transformation of skeletal muscle fiber types in Tibetan sheep and its influence on meat quality. It improves our understanding of the role of circRNAs in Tibetan sheep muscle development.

Circular RNA screening identifies circMYLK4 as a regulator of fast/slow myofibers in porcine skeletal muscles

The type of myofiber is related to the quality of meat. The slow oxidized myofiber helps to increase the tenderness and juiciness of muscle. Numerous studies have shown that circRNA plays a key role in skeletal muscle development. However, the role of circRNA in porcine skeletal myofiber types is unclear. In this study, we performed high-throughput RNA sequencing to study the differential expression of circRNA in the longissimus dorsi and the soleus muscle. A total of 40,757 circRNAs were identified, of which 181 were significantly different. Interestingly, some circRNAs were involved in metabolism pathways, AMPK, FoxO, and PI3K-Akt signaling pathways. Besides, we focused on a novel circRNA-circMYLK4. By injecting circMYLK4-AAV into piglets, we found that circMYLK4 significantly increased the mRNA and protein levels of the slow muscle marker genes. In summary, our study laid an essential foundation for further research of circRNA in myofiber type conversion and higher meat quality.

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.

Analysis of potential regulatory LncRNAs and CircRNAs in the oxidative myofiber and glycolytic myofiber of chickens

SART and PMM are mainly composed of oxidative myofibers and glycolytic myofibers, respectively, and myofiber types profoundly influence postnatal muscle growth and meat quality. SART and PMM are composed of lncRNAs and circRNAs that participate in myofiber type regulation. To elucidate the regulatory mechanism of myofiber type, lncRNA and circRNA sequencing was used to systematically compare the transcriptomes of the SART and PMM of Chinese female Qingyuan partridge chickens at their marketing age. The luminance value (L*), redness value (a*), average diameter, cross-sectional area, and density difference between the PMM and SART were significant (p < 0.05). ATPase staining results showed that PMMs were all darkly stained and belonged to the glycolytic type, and the proportion of oxidative myofibers in SART was 81.7%. A total of 5 420 lncRNAs were identified, of which 365 were differentially expressed in the SART compared with the PMM (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and KEGG pathways (p < 0.05), including striated muscle cell differentiation, regulation of cell proliferation, regulation of muscle cell differentiation, myoblast differentiation, regulation of myoblast differentiation, and MAPK signaling pathway. Pathways and coexpression network analyses suggested that XR_003077811.1, XR_003072304.1, XR_001465942.2, XR_001465741.2, XR_001470487.1, XR_003077673.1 and XR_003074785.1 played important roles in regulating oxidative myofibers by TBX3, QKI, MYBPC1, CALM2, and PPARGC1A expression. A total of 10 487 circRNAs were identified, of which 305 circRNAs were differentially expressed in the SART compared with the PMM (p < 0.05). Functional enrichment analysis showed that differentially expressed circRNAs were involved in host gene expression and were enriched in the AMPK, calcium signaling pathway, FoxO signaling pathway, p53 signaling pathway, and cellular senescence. Novel_circ_004282 and novel_circ_002121 played important roles in regulating oxidative myofibers by PPP3CA and NFATC1 expression. Using lncRNA-miRNA/circRNA-miRNA integrated analysis, we identified many candidate interaction networks that might affect muscle fiber performance. Important lncRNA-miRNA-mRNA networks, such as lncRNA-XR_003074785.1/miR-193-3p/PPARGC1A, regulate oxidative myofibers. This study reveals that lncXR_003077811.1, lncXR_003072304.1, lncXR_001465942.2, lncXR_001465741.2, lncXR_001470487.1, lncXR_003077673.1, XR_003074785.1, novel_circ_004282 and novel_circ_002121 might regulate oxidative myofibers. The lncRNA-XR_003074785.1/miR-193-3p/PPARGC1A pathway might regulate oxidative myofibers. All these findings provide rich resources for further in-depth research on the regulatory mechanism of lncRNAs and circRNAs in myofibers.

Comparison of long non-coding RNAs in adipose and muscle tissues between seven indigenous Chinese and the Yorkshire pig breeds

lncRNAs play crucial roles in fat metabolism in animals. Previously, we have compared the mRNA transcriptome profiles between seven fat-type Chinese pig breeds and one lean-type Western breed (Yorkshire, YY). The associations between differentially expressed (DE) genes and phenotypical traits were investigated. In the present study, to further explore the underlying regulatory mechanisms, lncRNAs were sequenced and compared between YY and Chinese indigenous breeds. The results showed 9114 and 7538 DE lncRNAs between at least one Chinese breed and the YY breed in the adipose and muscle tissue respectively. KEGG enrichment analysis revealed that the target genes of these DE lncRNAs mainly influenced the glucolipid metabolism, which is an important process affecting meat quality. Correlation analyses between the DE lncRNA and DE mRNA genes related to meat quality and growth traits were performed. The results showed that LTCONS_00073280 was associated with intramuscular fat content. Four lncRNAs (LTCONS_00101781, LTCONS_00037879, LTCONS_00088260 and LTCONS-00128343) might mediate backfat thickness. Overall, this study provides candidate lncRNAs that potentially affect meat quality, which might be useful for molecular breeding of pig breeds in future.

Function of Circular RNAs in Fish and Their Potential Application as Biomarkers

Circular RNAs (circRNAs) are an emerging class of regulatory RNAs with a covalently closed-loop structure formed during pre-mRNA splicing. Recent advances in high-throughput RNA sequencing and circRNA-specific computational tools have driven the development of novel approaches to their identification and functional characterization. CircRNAs are stable, developmentally regulated, and show tissue- and cell-type-specific expression across different taxonomic groups. They play a crucial role in regulating various biological processes at post-transcriptional and translational levels. However, the involvement of circRNAs in fish immunity has only recently been recognized. There is also broad evidence in mammals that the timely expression of circRNAs in muscle plays an essential role in growth regulation but our understanding of their expression and function in teleosts is still very limited. Here, we discuss the available knowledge about circRNAs and their role in growth and immunity in vertebrates from a comparative perspective, with emphasis on cultured teleost fish. We expect that the interest in teleost circRNAs will increase substantially soon, and we propose that they may be used as biomarkers for selective breeding of farmed fish, thus contributing to the sustainability of the aquaculture sector.

The Roles of LncRNAs in Osteogenesis, Adipogenesis and Osteoporosis

BACKGROUND

Osteoporosis (OP) is the most common bone disease, which is listed by the World Health Organization (WHO) as the third major threat to life and health among the elderly. The etiology of OP is multifactorial, and its potential regulatory mechanism remains unclear. Long non-coding RNAs (LncRNAs) are the non-coding RNAs that are over 200 bases in the chain length. Increasing evidence indicates that LncRNAs are the important regulators of osteogenic and adipogenic differentiation, and the occurrence of OP is greatly related to the dysregulation of the bone marrow mesenchymal stem cells (BMSCs) differentiation lineage. Meanwhile, LncRNAs affect the occurrence and development of OP by regulating OP-related biological processes.

METHODS

In the review, we summarized and analyzed the latest findings of LncRNAs in the pathogenesis, diagnosis and related biological processes of OP. Relevant studies published in the last five years were retrieved and selected from the PubMed database using the keywords of LncRNA and OP.

RESULTS/CONCLUSION

The present study aimed to examine the underlying mechanisms and biological roles of LncRNAs in OP, as well as osteogenic and adipogenic differentiation. Our results contributed to providing new clues for the epigenetic regulation of OP, making LncRNAs the new targets for OP therapy.

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.

Improvement, identification, and target prediction for miRNAs in the porcine genome by using massive, public high-throughput sequencing data

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.

An integration analysis based on genomic, transcriptomic and QTX information reveals credible candidate genes for fat-related traits in pigs

Meat quality improvement is of great interest to researchers in pig breeding and many researchers have identified abundant associated quantitative trait loci, genes and polymorphisms (QTXs) for fat-related traits. However, it is challenging to determine credible candidate genes from a mass of associations. The efficiency of identification of credible candidate genes in these QTXs is restricted by limited integration analyses of data from multiple omics. In this study, we constructed a 'candidate gene map' of fat-related traits in pigs based on published literature and the latest genome. In total, 6,861 QTXs, which covered 9,323 genes on the pig genome, were used. Combining the QTX hotspots and pathway analysis, we identified 180 candidate genes that may regulate the fat-related traits, and choose PNPLA2, PPARG, SREBF1, ACACA, PPARD and PPARA as credible candidate genes. In addition, we discussed the importance of incorporating transcriptome data and genomic data in causal gene identification, and the multi-omics information can effectively improve the credibility of identified candidate genes.

Comprehensive analysis of circular RNA profiles in skeletal muscles of aging mice and after aerobic exercise intervention

Aging induces gradual accumulation of damages in cells and tissues, which leads to physiological dysfunctions. Aging-associated muscle dysfunction is commonly seen in aged population and severely affects their physical activity and life quality, against which aerobic training has been shown to exert antagonizing or alleviating effects. Circular RNAs (circRNAs) play important roles in various physiological processes, yet their involvement in aging-associated muscle dysfunction is not well understood. In this study, we performed comprehensive analysis of circRNAs profiles in quadriceps muscles in sedentary young and aging mice, as well as aging mice with aerobic exercise using RNA sequencing. Our results identified circRNAs altered by factors of aging and aerobic exercise. Their host genes were then predicted and analyzed by gene ontology enrichment analysis. Importantly, we found that circBBS9 featured decreased levels in aging compared to young mice and elevated expression in exercise versus sedentary aging mice. Besides, we performed GO and KEGG analysis on circBBS9 target genes, as well as established the circBBS9-miRNA-mRNAs interaction network. Our results indicate that circBBS9 may play active roles in muscle aging by mediating the benefits of aerobic training intervention, thus may serve as a novel therapeutic target combating aging-associated muscle dysfunction.