LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development

Transcriptome analysis of long non-coding RNA associated with embryonic muscle development in chickens

1. Skeletal muscle is an important component of chicken carcass. In chickens, the number of muscle fibres is fixed during the embryonic period, and muscle development during the embryonic period determines the muscle development potential after hatching.2. Beijing-You (BY) and Cornish (CN) chickens show completely different growth rates and body types, and two breeds were used in this study to explore the role of lncRNAs in muscle development during different chicken embryonic periods. A systematic analysis of lncRNAs and mRNAs were conducted in the pectoral muscle tissues of BY and CN chickens at embryonic days 11 (ED11), 13 (ED13), 15 (ED15), 17 (ED17), and 1-day-old (D1) using RNA-seq. A total of 4,104 differentially expressed transcripts (DETs) were identified among the five stages, including 2,359 lncRNAs and 1,745 mRNAs.3. The number of DETs between the two breeds at ED17 (1,658 lncRNAs and 1,016 mRNAs) was much higher than the total number of DET at all the other stages (692 lncRNAs and 729 mRNAs), indicating that the two breeds show the largest difference in gene regulation at ED17.4. Correlation analysis was performed for all differentially expressed lncRNAs and mRNAs during the five periods. Forty-three, cis interaction pairs of lncRNA-mRNA related to chicken muscle development were predicted. The expression of four pairs was verified, and the results showed MSTRG.12395.2-FGFBP2 and MSTRG.18590.6-FMOD were significantly up-regulated in CN at ED11 compared to BY and might be important candidate genes for embryonic muscle development.

Whole transcriptome profiling reveals a lncMDP1 that regulates myogenesis by adsorbing miR-301a-5p targeting CHAC1

Myoblast proliferation and differentiation are essential for skeletal muscle development. In this study, we generated the expression profiles of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in different developmental stages of chicken primary myoblasts (CPMs) using RNA sequencing (RNA-seq) technology. The dual luciferase reporter system was performed using chicken embryonic fibroblast cells (DF-1), and functional studies quantitative real-time polymerase chain reaction (qPCR), cell counting kit-8 (CCK-8), 5-Ethynyl-2'-deoxyuridine (EdU), flow cytometry cycle, RNA fluorescence in situ hybridization (RNA-FISH), immunofluorescence, and western blotting assay. Our research demonstrated that miR-301a-5p had a targeted binding ability to lncMDP1 and ChaC glutathione-specific gamma-glutamylcyclotransferase 1 (CHAC1). The results revealed that lncMDP1 regulated the proliferation and differentiation of myoblasts via regulating the miR-301a-5p/CHAC1 axis, and CHAC1 promotes muscle regeneration. This study fulfilled the molecular regulatory network of skeletal muscle development and providing an important theoretical reference for the future improvement of chicken meat performance and meat quality.

Weighted gene co-expression network analysis identified hub genes critical to fatty acid composition in Gushi chicken breast muscle

BACKGROUND

The composition and content of fatty acids in the breast muscle are important factors influencing meat quality. In this study, we investigated the fatty acid composition and content in the breast muscle of Gushi chickens at different developmental stages (14 weeks, 22 weeks, and 30 weeks). Additionally, we utilized transcriptomic data from the same tissue and employed WGCNA and module identification methods to identify key genes associated with the fatty acid composition in Gushi chicken breast muscle and elucidate their regulatory networks.

RESULTS

Among them, six modules (blue, brown, green, light yellow, purple, and red modules) showed significant correlations with fatty acid content and metabolic characteristics. Enrichment analysis revealed that these modules were involved in multiple signaling pathways related to fatty acid metabolism, including fatty acid metabolism, PPAR signaling pathway, and fatty acid biosynthesis. Through analysis of key genes, we identified 136 genes significantly associated with fatty acid phenotypic traits. Protein-protein interaction network analysis revealed that nine of these genes were closely related to fatty acid metabolism. Additionally, through correlation analysis of transcriptome data, we identified 51 key ceRNA regulatory networks, including six central genes, 7 miRNAs, and 28 lncRNAs.

CONCLUSION

This study successfully identified key genes closely associated with the fatty acid composition in Gushi chicken breast muscle, as well as their post-transcriptional regulatory networks. These findings provide new insights into the molecular regulatory mechanisms underlying the flavor characteristics of chicken meat and the composition of fatty acids in the breast muscle.

Transcriptome Analysis Reveals the Profile of Long Non-Coding RNAs during Myogenic Differentiation in Goats

The long non-coding RNAs (lncRNAs) are emerging as essential regulators of the growth and development of skeletal muscles. However, little is known about the expression profiles of lncRNAs during the proliferation and differentiation of skeletal muscle satellite cells (MuSCs) in goats. In this study, we investigate potential regulatory lncRNAs that govern muscle development by performing lncRNA expression profiling analysis during the proliferation (cultured in the growth medium, GM) and differentiation (cultured in the differentiation medium, DM1/DM5) of MuSCs. In total, 1001 lncRNAs were identified in MuSC samples, and 314 differentially expressed (DE) (FDR < 0.05, |log2FC| > 1) lncRNAs were screened by pairwise comparisons from three comparison groups (GM-vs-DM1, GM-vs-DM5, DM1-vs-DM5). Moreover, we identified the cis-, trans-, and antisense-regulatory target genes of DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these target genes were significantly enriched in muscle development-related GO terms and KEGG pathways. In addition, the network of interactions between DE lncRNAs and their target genes was identified, which included well-known myogenesis regulators such as Myogenic differentiation 1 (MyoD), Myogenin (MyoG), and Myosin heavy chain (MyHC). Meanwhile, competing endogenous RNA (ceRNA) network analysis showed that 237 DE lncRNAs could bind to 329 microRNAs (miRNAs), while miRNAs could target 564 mRNAs. Together, our results provide a genome-wide resource of lncRNAs that may contribute to myogenic differentiation in goats and lay the groundwork for future investigation into their functions during skeletal muscle development.

CircRNAs Related to Breast Muscle Development and Their Interaction Regulatory Network in Gushi Chicken

Circular RNAs (circRNAs) play a significant regulatory role during skeletal muscle development. To identify circRNAs during postnatal skeletal muscle development in chickens, we constructed 12 cDNA libraries from breast muscle tissues of Chinese Gushi chickens at 6, 14, 22, and 30 weeks and performed RNA sequencing. In total, 2112 circRNAs were identified, and among them 79.92% were derived from exons. CircRNAs are distributed on all chromosomes of chickens, especially chromosomes 1-9 and Z. Bioinformatics analysis showed that each circRNA had an average of 38 miRNA binding sites, 61.32% of which have internal ribosomal entry site (IRES) elements. Furthermore, in total 543 differentially expressed circRNAs (DE-circRNAs) were identified. Functional enrichment analysis revealed that DE-circRNAs source genes are engaged in biological processes and muscle development-related pathways; for example, cell differentiation, sarcomere, and myofibril formation, mTOR signaling pathway, and TGF-β signaling pathway, etc. We also established a competitive endogenous RNA (ceRNA) regulatory network associated with skeletal muscle development. The results in this report indicate that circRNAs can mediate the development of chicken skeletal muscle by means of a complex ceRNA network among circRNAs, miRNAs, genes, and pathways. The findings of this study might help increase the number of known circRNAs in skeletal muscle tissue and offer a worthwhile resource to further investigate the function of circRNAs in chicken skeletal muscle development.

Identification and characterization of long non-coding RNAs in juvenile and adult skeletal muscle of largemouth bass (Micropterus salmoides)

Long non-coding RNAs (lncRNAs) are a class of transcriptional RNA molecules, which play critical roles in diverse biological processes. However, little is known about the overall expression pattern and roles of lncRNAs in skeletal muscle of largemouth bass (LMB). Here, we constructed two skeletal muscle RNA libraries to find lncRNAs that may involve in the regulation of skeletal muscle development between juvenile and adult LMB. A total of 16,147 lncRNAs and 4611 differentially expressed lncRNAs were identified. Among these identified lncRNAs, 10 lncRNAs were randomly selected to confirm their expression by real-time qPCR both in libraries, which were consistent with the RNA sequencing results. The target mRNAs of lncRNAs were predicted for GO enrichment analysis. Results showed that these targets associated with growth and development of muscle, such as skeletal muscle fiber development, myoblast proliferation and differentiation. Importantly, correlation analysis of lncRNA-miRNA-mRNA regulatory network revealed that several lncRNAs targeted miRNAs which are closely involved in the regulation of muscle development. It is the first time to identify a number of lncRNA that correlate with skeletal muscle development in LMB. Our results not only provide a comprehensive expression profile of muscle lncRNAs in this species, but also provide a theoretical basis for further elaborating genetic regulation mechanism of muscle growth and development, and pave the way for the future molecular assisted breeding in carnivorous fishes.

Regulation of Non-Coding RNA in the Growth and Development of Skeletal Muscle in Domestic Chickens

Chicken is the most widely consumed meat product worldwide and is a high-quality source of protein for humans. The skeletal muscle, which accounts for the majority of chicken products and contains the most valuable components, is tightly correlated to meat product yield and quality. In domestic chickens, skeletal muscle growth is regulated by a complex network of molecules that includes some non-coding RNAs (ncRNAs). As a regulator of muscle growth and development, ncRNAs play a significant function in the development of skeletal muscle in domestic chickens. Recent advances in sequencing technology have contributed to the identification and characterization of more ncRNAs (mainly microRNAs (miRNAs), long non-coding RNAs (LncRNAs), and circular RNAs (CircRNAs)) involved in the development of domestic chicken skeletal muscle, where they are widely involved in proliferation, differentiation, fusion, and apoptosis of myoblasts and satellite cells, and the specification of muscle fiber type. In this review, we summarize the ncRNAs involved in the skeletal muscle growth and development of domestic chickens and discuss the potential limitations and challenges. It will provide a theoretical foundation for future comprehensive studies on ncRNA participation in the regulation of skeletal muscle growth and development in domestic chickens.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Embryo Muscle of Chicken

Embryonic muscle development determines the state of muscle development and muscle morphological structure size. Recent studies have found that long non-coding RNAs (lncRNAs) could influence numerous cellular processes and regulated growth and development of flora and fauna. A total of 1056 differentially expressed lncRNAs were identified by comparing the different time points during embryonic muscle development, which included 874 new lncRNAs. Here, we found that there were different gene expression patterns on the 12th day of embryo development (E12). Herein, WGCNA and correlation analyses were used to predict lncRNA function on E12 through the screening and identification of lncRNAs related to muscle development in the embryo leg muscles of Chengkou mountain chickens at different times. GO and KEGG functional enrichment analysis was performed on target genes involved in cis-regulation and trans-regulation. An interaction network diagram was constructed based on the muscle development pathways, such as Wnt, FoxO, and PI3K-AKT signaling pathways, to determine the interaction between mRNAs and lncRNAs. This study preliminarily determined the lncRNA expression pattern of muscle development during the middle and late embryonic stages of Chengkou mountain chickens, and provided a basis to analyze the molecular mechanism of muscle development.

Weighted gene co-expression network indicates that the DYNLL2 is an important regulator of chicken breast muscle development and is regulated by miR-148a-3p

Background

The characteristics of muscle fibers determine the growth and meat quality of poultry. In this study, we performed a weighted gene co-expression network analysis (WGCNA) on the muscle fiber characteristics and transcriptome profile of the breast muscle tissue of Gushi chicken at 6, 14, 22, and 30 weeks.

Results

A total of 27 coexpressed biological functional modules were identified, of which the midnight blue module had the strongest correlation with muscle fiber and diameter. In addition, 7 hub genes were found from the midnight blue module, including LC8 dynein light chain 2 (DYNLL2). Combined with miRNA transcriptome data, miR-148a-3p was found to be a potential target miRNA of DYNLL2. Experiments on chicken primary myoblasts (CPMs) demonstrated that miR-148a-3p promotes the expression of myosin heavy chain (MYHC) protein by targeting DYNLL2, proving that it can promote differentiation of myoblasts.

Conclusions

This study proved that the hub gene DYNLL2 and its target miR-148-3p are important regulators in chicken myogenesis. These results provide novel insights for understanding the molecular regulation mechanisms related to the development of chicken breast muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08522-8.

miR-30a-3p can inhibit the proliferation and promote the differentiation of chicken primary myoblasts

1. Chicken muscle is an important factor in meat quality and its development is controlled by a complex regulatory network.2. The following study examined the expression of miR-30a-3p in Gushi chicken breast muscle tissue and found that it was differentially expressed at different embryonic stages, reaching a peak in the 14-day-old embryo (E14).3. The effect of miR-30a-3p on chicken primary myoblasts (CPMs) was explored. Results from both cell counting kit-8 (CCK-8) and 5-ethynyl-2'-deoxyuridine (EdU) showed that this can inhibit the proliferation of myoblasts, and through cell cycle experiments, the inhibition of myoblast proliferation was found, which may be due to G0/G1 arrest in the cell cycle.4. The effect of miR-30a-3p on the differentiation of myoblasts was studied. The results showed that miR-30a-3p can promote the expression of MYOD, myogenin (MYOG), and myosin heavy chain (MYHC) genes to promote the differentiation of myoblasts. Through MYHC protein immunofluorescence experiments, it was found that miR-30a-3p can effectively increase the area of myotubes.5. Finally, mRNA transcriptome data was analysed, which showed that miR-30a-3p has 51 potential target genes. Among them, forkhead box O3 (FOXO3), ankyrin repeat domain 1 (ANKRD1), and insulin-induced 1 (INSIG1) genes were differentially expressed at different developmental stages and were enriched in Gene Ontology (GO) terms, such as cell differentiation and cellular developmental process. The data showed that tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein gamma (YWHAG), BUB1 mitotic checkpoint serine/threonine kinase (BUB1), and growth arrest and DNA damage-inducible 45 (GADD45) genes were enriched in the cell cycle pathway.6. It can be speculated that miR-30a-3p plays roles through these genes in myoblast development. This research provides information for further improving knowledge of the chicken muscle development regulation network.

Integrative Analysis of lncRNA-mRNA Co-expression Provides Novel Insights Into the Regulation of Developmental Transitions in Female Varroa destructor

Varroa destructor is a major pathogenic driver of the Western honeybee colony losses globally. Understanding the developmental regulation of V. destructor is critical to develop effective control measures. Development is a complex biological process regulated by numerous genes and long non-coding RNAs (lncRNAs); however, the underlying regulation of lncRNAs in the development of V. destructor remains unknown. In this study, we analyzed the RNA sequencing (RNA-Seq) data derived from the four stages of female V. destructor in the reproductive phase (i.e., egg, protonymph, deutonymph, and adult). The identified differentially expressed mRNAs and lncRNAs exhibited a stage-specific pattern during developmental transitions. Further functional enrichment established that fat digestion and absorption, ATP-binding cassette (ABC) transporters, mitogen-activated protein kinase (MAPK) signaling pathway, and ubiquitin-proteasome pathway play key roles in the maturation of female V. destructor. Moreover, the lncRNAs and mRNAs of some pivotal genes were significantly upregulated at the deutonymph stage, such as cuticle protein 65/6.4/63/38 and mucin 5AC, suggesting that deutonymph is the key stage of metamorphosis development and pathogen resistance acquisition for female V. destructor. Our study provides novel insights into a foundational understanding of V. destructor biology.

Differentially Expressed lncRNAs Related to the Development of Abdominal Fat in Gushi Chickens and Their Interaction Regulatory Network

Chickens are one of the most important sources of meat worldwide, and the growth status of abdominal fat is closely related to production efficiency. Long noncoding RNAs (lncRNAs) play an important role in lipid metabolism and deposition regulation. However, research on the expression profile of lncRNAs related to the development of abdominal fat in chickens after hatching and their interaction regulatory networks is still lacking. To characterize the lncRNA expression profile during the development of chicken abdominal fat, abdominal adipose tissues from 6-, 14-, 22-, and 30-week-old Chinese Gushi chickens were herein used to construct 12 cDNA libraries, and a total of 3,827 new lncRNAs and 5,466 previously annotated lncRNAs were revealed. At the same time, based on the comparative analysis of five combinations, 276 differentially expressed lncRNAs (DE-lncRNAs) were screened. Functional enrichment analysis showed that the predicted target genes of these DE-lncRNAs were significantly enriched in pathways related to the posttranscriptional regulation of gene expression, negative regulation of cell proliferation, cell adhesion and other biological processes, glycosphingolipid biosynthesis, PPAR signaling, fatty acid degradation, fatty acid synthesis and others. In addition, association analysis of the lncRNA transcriptome profile was performed, and DE-lncRNA-related lncRNA-mRNA, lncRNA-miRNA and lncRNA-miRNA-mRNA interaction regulatory networks were constructed. The results showed that DE-lncRNA formed a complex network with PPAR pathway components, including PPARD, ACOX1, ADIPOQ, CPT1A, FABP5, ASBG2, LPL, PLIN2 and related miRNAs, including mir-200b-3p, mir-130b-3p, mir-215-5p, mir-122-5p, mir-223 and mir-125b-5p, and played an important regulatory role in biological processes such as lipid metabolism, adipocyte proliferation and differentiation. This study described the dynamic expression profile of lncRNAs in the abdominal fat of Gushi chickens for the first time and constructed the DE-lncRNA interaction regulatory network. The results expand the number of known lncRNAs in chicken abdominal fat and provide valuable resources for further elucidating the posttranscriptional regulatory mechanism of chicken abdominal fat development or deposition.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

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.