LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development
Background: Chicken skeletal muscle is an important economic product. The late stages of chicken development constitute the main period that affects meat production. LncRNAs play important roles in controlling the epigenetic process of growth and development. However, studies on the role of lncRNAs in the late stages of chicken breast muscle development are still lacking. In this study, to investigate the expression characteristics of lncRNAs during chicken muscle development, 12 cDNA libraries were constructed from Gushi chicken breast muscle samples from 6-, 14-, 22-, and 30-week-old chickens.
Results: A total of 1,252 new lncRNAs and 1,376 annotated lncRNAs were identified. Furthermore, 53, 61, 50, 153, 117, and 78 DE-lncRNAs were found in the W14 vs. W6, W22 vs. W14, W22 vs. W6, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparison groups, respectively. After GO enrichment analysis of the DE-lncRNAs, several muscle development-related GO terms were found in the W22 vs. W14 comparison group. Moreover, it was found that the MAPK signaling pathway was one of the most frequently enriched pathways in the different comparison groups. In addition, 12 common target DE-miRNAs of DE-lncRNAs were found in different comparison groups, some of which were muscle-specific miRNAs, such as gga-miR-206, gga-miR-1a-3p, and miR-133a-3p. Interestingly, the precursors of four newly identified miRNAs were found to be homologous to lncRNAs. Additionally, we found some ceRNA networks associated with muscle development-related GO terms. For example, the ceRNA networks contained the DYNLL2 gene with 12 lncRNAs that targeted 2 miRNAs. We also constructed PPI networks, such as IGF-I-EGF and FZD6-WNT11.
Conclusions: This study revealed, for the first time, the dynamic changes in lncRNA expression in Gushi chicken breast muscle at different periods and revealed that the MAPK signaling pathway plays a vital role in muscle development. Furthermore, MEF2C and its target lncRNA may be involved in muscle regulation through the MAPK signaling pathway. This research provided valuable resources for elucidating posttranscriptional regulatory mechanisms to promote the development of chicken breast muscles after hatching.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
This is a list of supplementary files associated with this preprint. Click to download.
Figure S1. Heatmap showing DE-lncRNAs from different stages.
Figure S2. The chromosome distribution of DE-lncRNAs from different stages.
Figure S3. The enriched GO terms of the DE-lncRNA. (A-B) Cis-target genes in the W14 vs. W6 and W30 vs. W22 comparison groups. (C-D) Trans-target genes in the W14 vs. W6 and W30 vs. W22 comparison groups.
Figure S4. The enriched KEGG pathways of the DE-lncRNA. (A-C) Cis-target genes in the W14 vs. W6, W22 vs. W14, and W30 vs. W22 comparison groups. (D) Trans-target genes in the W14 vs. W6 comparison groups.
Figure S5. The lncRNA-miRNA-mRNA ceRNA networks. (A) W14 vs. W6 comparison group; (B) W22 vs. W14 comparison group; (C) W30 vs. W22 comparison group.
Table S1. Summary of draft reads of 12 cDNA libraries, determined by RNA sequencing. Abbreviations: W6_1, sample 1 of 6 weeks; W6_2, sample 2 of 6 weeks; W6_3, sample 3 of 6 weeks; W14_1, sample 1 of 14 weeks; W14_2, sample 2 of 14 weeks; W14_3, sample 3 of 14 weeks; W22_1, sample 1 of 22 weeks; W22_2, sample 2 of 22 weeks; W22_3, sample 3 of 22 weeks; W30_1, sample 1 of 30 weeks; W30_2, sample 2 of 30 weeks; W30_3, sample 3 of 30 weeks.
Table S2. Cis-regulatory interactions between DE-lncRNAs and DE-mRNAs.
Table S3. Trans-regulatory interactions between DE-lncRNAs and DE-mRNAs.
Table S4. The interaction between DE-lncRNAs and DE-miRNAs.
Table S5. Details on the novel miRNAs identified in this study. Abbreviations: W6, W14, W22, and W30 represent small RNA libraries obtained using samples from chickens aged 6, 14, 22, and 30 weeks, respectively.
Table S6. qRT-PCR primers. Abbreviation: AT refers to the annealing temperature; F and R refer to the forward and reverse primers, respectively.
Posted 09 Jan, 2021
On 15 Dec, 2020
Received 14 Dec, 2020
On 22 Nov, 2020
On 27 Jul, 2020
Invitations sent on 27 Jul, 2020
On 26 Jul, 2020
On 23 Jul, 2020
On 19 Jul, 2020
Received 13 Jul, 2020
On 26 May, 2020
Received 05 Apr, 2020
On 20 Mar, 2020
Invitations sent on 17 Mar, 2020
On 24 Feb, 2020
On 23 Feb, 2020
On 23 Feb, 2020
LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development
Posted 09 Jan, 2021
On 15 Dec, 2020
Received 14 Dec, 2020
On 22 Nov, 2020
On 27 Jul, 2020
Invitations sent on 27 Jul, 2020
On 26 Jul, 2020
On 23 Jul, 2020
On 19 Jul, 2020
Received 13 Jul, 2020
On 26 May, 2020
Received 05 Apr, 2020
On 20 Mar, 2020
Invitations sent on 17 Mar, 2020
On 24 Feb, 2020
On 23 Feb, 2020
On 23 Feb, 2020
Background: Chicken skeletal muscle is an important economic product. The late stages of chicken development constitute the main period that affects meat production. LncRNAs play important roles in controlling the epigenetic process of growth and development. However, studies on the role of lncRNAs in the late stages of chicken breast muscle development are still lacking. In this study, to investigate the expression characteristics of lncRNAs during chicken muscle development, 12 cDNA libraries were constructed from Gushi chicken breast muscle samples from 6-, 14-, 22-, and 30-week-old chickens.
Results: A total of 1,252 new lncRNAs and 1,376 annotated lncRNAs were identified. Furthermore, 53, 61, 50, 153, 117, and 78 DE-lncRNAs were found in the W14 vs. W6, W22 vs. W14, W22 vs. W6, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparison groups, respectively. After GO enrichment analysis of the DE-lncRNAs, several muscle development-related GO terms were found in the W22 vs. W14 comparison group. Moreover, it was found that the MAPK signaling pathway was one of the most frequently enriched pathways in the different comparison groups. In addition, 12 common target DE-miRNAs of DE-lncRNAs were found in different comparison groups, some of which were muscle-specific miRNAs, such as gga-miR-206, gga-miR-1a-3p, and miR-133a-3p. Interestingly, the precursors of four newly identified miRNAs were found to be homologous to lncRNAs. Additionally, we found some ceRNA networks associated with muscle development-related GO terms. For example, the ceRNA networks contained the DYNLL2 gene with 12 lncRNAs that targeted 2 miRNAs. We also constructed PPI networks, such as IGF-I-EGF and FZD6-WNT11.
Conclusions: This study revealed, for the first time, the dynamic changes in lncRNA expression in Gushi chicken breast muscle at different periods and revealed that the MAPK signaling pathway plays a vital role in muscle development. Furthermore, MEF2C and its target lncRNA may be involved in muscle regulation through the MAPK signaling pathway. This research provided valuable resources for elucidating posttranscriptional regulatory mechanisms to promote the development of chicken breast muscles after hatching.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10