Bahareldin Ali Abdalla

Integrated Analysis of Long Non-coding RNAs (LncRNAs) and mRNA Expression Profiles Reveals the Potential Role of LncRNAs in Skeletal Muscle Development of the Chicken

Long non-coding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation. However, little is currently known about the mechanisms by which they regulate skeletal muscle development in the chicken. In this study, we used RNA sequencing to profile the leg muscle transcriptome (lncRNA and mRNA) at three stages of skeletal muscle development in the chicken: embryonic day 11 (E11), embryonic day 16 (E16), and 1 day after hatching (D1). In total, 129, 132, and 45 differentially expressed lncRNAs, and 1798, 3072, and 1211 differentially expressed mRNAs were identified in comparisons of E11 vs. E16, E11 vs. D1, and E16 vs. D1, respectively. Moreover, we identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and constructed lncRNA-gene interaction networks. In total, 126 and 200 cis-targets, and two and three trans-targets were involved in lncRNA-gene interaction networks that were constructed based on the E11 vs. E16, and E11 vs. D1 comparisons, respectively. The comparison of the E16 vs. D1 lncRNA-gene network comprised 25 cis-targets. We determined that lncRNA target genes are potentially involved in cellular development, and cellular growth and proliferation using Ingenuity Pathway Analysis. The gene networks identified for the E11 vs. D1 comparison were involved in embryonic development, organismal development and tissue development. The present study provides an RNA sequencing based evaluation of lncRNA function during skeletal muscle development in the chicken. Comprehensive analysis facilitated the identification of lncRNAs and target genes that might contribute to the regulation of different stages of skeletal muscle development.

Genome-wide association study of aggressive behaviour in chicken

In the poultry industry, aggressive behaviour is a large animal welfare issue all over the world. To date, little is known about the underlying genetics of the aggressive behaviour. Here, we performed a genome-wide association study (GWAS) to explore the genetic mechanism associated with aggressive behaviour in chickens. The GWAS results showed that a total of 33 SNPs were associated with aggressive behaviour traits (P < 4.6E-6). rs312463697 on chromosome 4 was significantly associated with aggression (P = 2.10905E-07), and it was in the intron region of the sortilin-related VPS10 domain containing receptor 2 (SORCS2) gene. In addition, biological function analysis of the nearest 26 genes around the significant SNPs was performed with Ingenuity Pathway Analysis. An interaction network contained 17 genes was obtained and SORCS2 was involved in this network, interacted with nerve growth factor (NGF), nerve growth factor receptor (NGFR), dopa decarboxylase (L-dopa) and dopamine. After knockdown of SORCS2, the mRNA levels of NGF, L-dopa and dopamine receptor genes DRD1, DRD2, DRD3 and DRD4 were significantly decreased (P < 0.05). In summary, our data indicated that SORCS2 might play an important role in chicken aggressive behaviour through the regulation of dopaminergic pathways and NGF.

LncRNA-Six1 Encodes a Micropeptide to Activate Six1 in Cis and Is Involved in Cell Proliferation and Muscle Growth.

Long non-coding RNAs (lncRNAs) play important roles in epigenetic regulation of skeletal muscle development. In our previous RNA-seq study (accession number GSE58755), we found that lncRNA-Six1 is an lncRNA that is differentially expressed between White Recessive Rock (WRR) and Xinghua (XH) chicken. In this study, we have further demonstrated that lncRNA-Six1 is located 432 bp upstream of the gene encoding the protein Six homeobox 1 (Six1). A dual-luciferase reporter assay identified that lncRNA-Six1 overlaps the Six1 proximal promoter. In lncRNA-Six1, a micropeptide of about 7.26 kDa was found to play an important role in the lncRNA-Six1 in cis activity. Overexpression of lncRNA-Six1 promoted the mRNA and protein expression level of the Six1 gene, while knockdown of lncRNA-Six1 inhibited Six1 expression. Moreover, tissue expression profiles showed that both the lncRNA-Six1 and the Six1 mRNA were highly expressed in chicken breast tissue. LncRNA-Six1 overexpression promoted cell proliferation and induced cell division. Conversely, its loss of function inhibited cell proliferation and reduced cell viability. Similar effects were observed after overexpression or knockdown of the Six1 gene. In addition, overexpression or knockdown of Six1 promoted or inhibited, respectively, the expression levels of muscle-growth-related genes, such as MYOG, MYHC, MYOD, IGF1R, and INSR. Taken together, these data demonstrate that lncRNA-Six1 carries out cis-acting regulation of the protein-encoding Six1 gene, and encodes a micropeptide to activate Six1 gene, thus promoting cell proliferation and being involved in muscle growth.

miR-16 controls myoblast proliferation and apoptosis through directly suppressing Bcl2 and FOXO1 activities

Myogenesis mainly involves several steps including myoblast proliferation, differentiation, apoptosis and fusion. Except for muscle specific regulators, few miRNAs were proved to coordinate this complex process. Here, we reported that miR-16 inhibited myoblast proliferation and promoted myoblast apoptosis by directly targeting Bcl2 and FOXO1. The expression level of miR-16 was significantly decreased in the hypertrophic pectoral muscle compared to the normal pectoral muscle in chicken. In vitro, elevating miR-16 significantly inhibited myoblast proliferation and promoted myoblast apoptosis, resulting in about 11.2% cells arrested in G1 phase and 12.3% apoptotic cells in the early stage. Bioinformatic and biochemical analyses revealed Bcl2 and FOXO1 as direct targets of miR-16. Consist to the effect of miR-16 on myogenesis, specific inhibition of Bcl2 or FOXO1 significantly suppressed myoblast proliferation and induced myoblast apoptosis, indicating that both Bcl2 and FOXO1 contributed to miR-16 regulatory function in myogenesis. Interestingly, FOXO1, as the core target, mediated multiple growth-related pathways induced by miR-16 such as PI3K-AKT-MAPK and PI3K-AKT-mTOR. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) revealed that 234 annotated genes bound by FOXO1 in the early-differentiated myoblasts, which were significantly enriched in myogenic proliferation, death and hypotrophy. Altogether, we proposed that miR-16 acted as a coordinated mediator to suppress myogenesis in avian through the control of myoblast proliferation and apoptosis. These findings have provided a novel mechanism whereby miR-16 represses Bcl2 and FOXO1 expression to maintain myoblast growth and skeletal muscle mass.

Circular RNAs are abundant and dynamically expressed during embryonic muscle development in chickens

Abstract The growth and development of skeletal muscle is regulated by proteins as well as non-coding RNAs. Circular RNAs (circRNAs) are universally expressed in various tissues and cell types, and regulate gene expression in eukaryotes. To identify the circRNAs during chicken embryonic skeletal muscle development, leg muscles of female Xinghua (XH) chicken at three developmental time points 11 embryo age (E11), 16 embryo age (E16) and 1 day post hatch (P1) were performed RNA sequencing. We identified 13,377 circRNAs with 3,036 abundantly expressed and most were derived from coding exons. A total of 462 differentially expressed circRNAs were identified (fold change > 2; q-value < 0.05). Parental genes of differentially expressed circRNAs were related to muscle biological processes. There were 946 exonic circRNAs have been found that harbored one or more miRNA-binding site for 150 known miRNAs. We validated that circRBFOX2s promoted cell proliferation through interacted with miR-206. These data collectively indicate that circRNAs are abundant and dynamically expressed during embryonic muscle development and could play key roles through sequestering miRNAs as well as other functions.

MiR-34b-5p Suppresses Melanoma Differentiation-Associated Gene 5 (MDA5) Signaling Pathway to Promote Avian Leukosis Virus Subgroup J (ALV-J)-Infected Cells Proliferaction and ALV-J Replication

Avian leukosis virus subgroup J (ALV-J) is an oncogenic retrovirus that has a similar replication cycle to multiple viruses and therefore can be used as a model system for viral entry into host cells. However, there are few reports on the genes or microRNAs (miRNAs) that are responsible for the replication of ALV-J. Our previous miRNA and RNA sequencing data showed that the expression of miR-34b-5p was significantly upregulated in ALV-J-infected chicken spleens compared to non-infected chicken spleens, but melanoma differentiation-associated gene 5 (MDA5) had the opposite expression pattern. In this study, a dual-luciferase reporter assay showed that MDA5 is a direct target of miR-34b-5p. In vitro, overexpression of miR-34b-5p accelerated the proliferation of ALV-J-infected cells by inducing the progression from G2 to S phase and it promoted cell migration. Ectopic expression of MDA5 inhibited ALV-J-infected cell proliferation, the cell cycle and cell migration, and knockdown of MDA5 promoted proliferation, the cell cycle and migration. In addition, during ALV-J infections, MDA5 can detect virus invasion and it triggers the MDA5 signaling pathway. MDA5 overexpression can activate the MDA5 signaling pathway, and thus it can inhibit the mRNA and protein expression of the ALV-J env gene and it can suppress virion secretion. In contrast, in response to the knockdown of MDA5 by small interfering RNA (siRNA) or an miR-34b-5p mimic, genes in the MDA5 signaling pathway were significantly downregulated (P < 0.05), but the mRNA and protein expression of ALV-J env and the sample-to-positive ratio of virion in the supernatants were increased. This indicates that miR-34b-5p is able to trigger the MDA5 signaling pathway and affect ALV-J infections. Together, these results suggest that miR-34b-5p targets MDA5 to accelerate the proliferation and migration of ALV-J-infected cells, and it promotes ALV-J replication, via the MDA5 signaling pathway.

miRNA-223 upregulated by MYOD inhibits myoblast proliferation by repressing IGF2 and facilitates myoblast differentiation by inhibiting ZEB1

Skeletal muscle differentiation can be regulated by various transcription factors and non-coding RNAs. In our previous work, miR-223 is differentially expressed in the skeletal muscle of chicken with different growth rates, but its role, expression and action mechanism in muscle development still remains unknown. Here, we found that MYOD transcription factor can upregulate miR-223 expression by binding to an E-box region of the gga-miR-223 gene promoter during avian myoblast differentiation. IGF2 and ZEB1 are two target genes of miR-223. The target inhibition of miR-223 on IGF2 and ZEB1 are dynamic from proliferation to differentiation of myoblast. miR-223 inhibits IGF2 expression only in the proliferating myoblast, whereas it inhibits ZEB1 mainly in the differentiating myoblast. The inhibition of IGF2 by miR-223 resulted in the repression of myoblast proliferation. During myoblast differentiation, miR-223 would be upregulated owing to the promoting effect of MYOD, and the upregulation of miR-223 would inhibit ZEB1 to promote myoblast differentiation. These results not only demonstrated that the well-known muscle determination factor MYOD can promote myoblast differentiation by upregulate miR-223 transcription, but also identified that miR-223 can influence myoblast proliferation and differentiation by a dynamic manner regulates the expression of its target genes.

MiR-16-5p targets SESN1 to regulate the p53 signaling pathway, affecting myoblast proliferation and apoptosis, and is involved in myoblast differentiation.

SummaryThe proliferation, apoptosis, and differentiation of myoblasts are essential processes in skeletal muscle development. During this developmental process, microRNAs (miRNAs) play crucial roles. In our previous RNA-seq study (accession number GSE62971), we found that miR-16-5p was differentially expressed between fast and slow growth in chicken. In this study, we report that miR-16-5p could inhibit myoblast proliferation, promote myoblast apoptosis, and repress myoblast differentiation by directly binding to the 3′ UTR of SESN1, which is also differentially expressed. Overexpression of SESN1 significantly promoted the proliferation, inhibited apoptosis, and induced differentiation of myoblasts. Conversely, its loss of function hampered myoblast proliferation, facilitated myoblast apoptosis, and inhibited myoblast differentiation. Interestingly, we found SESN1 could regulate p53 by a feedback mechanism, thereby participating in the regulation of p53 signaling pathway, which suggests that this feedback is indispensable for myoblast proliferation and apoptosis. Altogether, these data demonstrated that miR-16-5p directly targets SESN1 to regulate the p53 signaling pathway, and therefore affecting myoblast proliferation and apoptosis. Additionally, SESN1 targets myogenic genes to control myoblast differentiation.

Characterization of miRNA and their target gene during chicken embryo skeletal muscle development

MicroRNAs (miRNAs) are non-coding RNAs that regulate mRNA expression by degradation or translational inhibition. We investigated the underlying molecular mechanisms of skeletal muscle development based on differentially expressed genes and miRNAs. We compared mRNA and miRNA from chicken skeletal muscle at embryonic day E11, E16 and one day post-hatch (P1). The interaction networks were constructed, according to target prediction results and integration analysis of up-regulated genes with down regulated miRNAs or down-regulated genes with up-regulated miRNAs with |log2fold change| ≥ 1.75, P < 0.005. The miRNA-mRNA integration analysis showed high number of mRNAs regulated by a few number of miRNAs. In the E11_VS_E16, comparison group we identified biological processes including muscle maintenance, myoblast proliferation and muscle thin filament formation. The E11_VS_P1 group comparison included negative regulation of axon extension, sarcomere organization, and cell redox homeostasis and kinase inhibitor activity. The E16_VS_P1 comparison group contained genes for the negative regulation of anti-apoptosis and axon extension as well as glomerular basement membrane development. Functional in vitro assays indicated that over expression of miR-222a and miR-126–5p in DF-1 cells significantly reduced the mRNA levels of the target genes CPEB3 and FGFR3, respectively. These integrated analyses provide several candidates for future studies concerning miRNAs-target function on regulation of embryonic muscle development and growth.

Molecular characterization, expression profile of the FSHR gene and its association with egg production traits in muscovy duck

Follicle-stimulating hormone (FSH) and its receptor play a key role in the follicular development and regulation of steroidogenesis in the ovary and spermatogenesis in the testis. The purpose of this study was to characterize the muscovy duck FSHR gene, identify SNPs and their association with egg production traits in muscovy ducks. Here, we cloned the complementary DNA (cDNA) sequence of FSHR, and examined the expression patterns of FSHR gene in adult female muscovy duck tissues. The cloned cDNA of the muscovy duck FSHR gene shared high similarity to those of pekin duck (Anas platyrhynchos) (95.7%) and chicken (93.2%). Three different muscovy duck FSHR transcripts were identified. Quantitative real-time PCR (RT-qPCR) results showed that the FSHR gene was expressed in all the 14 tested tissues, and the highest expression level was seen in the ovary. A total of 16 SNPs were identified, among which, four SNPs were located in the coding region of FSHR. The SNP C320T is significantly associated with egg production at 59 weeks of age (