RNA-Seq transcriptome analysis of breast muscle in Pekin ducks supplemented with the dietary probiotic Clostridium butyricum
BMC Genomics volume 19, Article number: 844 (2018)
Increased attention is being paid to breast muscle yield and meat quality in the duck breeding industry. Our previous report has demonstrated that dietary Clostridium butyricum (C. butyricum) can improve meat quality of Pekin ducks. However, the potential biological processes and molecular mechanisms that are modulated by dietary C. butyricum in the breast muscle of Pekin ducks remain unknown.
Supplementation with C. butyricum increased growth performance and meat yield. Therefore, we utilized de novo assembly methods to analyze the RNA-Seq transcriptome profiles in breast muscle to explore the differentially expressed genes between C. butyricum-treated and control Pekin ducks. A total of 1119 differentially expressed candidate genes were found of which 403 genes were significantly up-regulated and 716 genes were significantly down-regulated significantly. qRT-PCR analysis was used to confirm the accuracy of the of RNA-Seq results. GO annotations revealed potential genes, processes and pathways that may participate in meat quality and muscle development. KEGG pathway analysis showed that the differentially expressed genes participated in numerous pathways related to muscle development, including ECM-receptor interaction, the MAPK signaling pathway and the TNF signaling pathway.
This study suggests that long-time dietary supplementation with C. butyricum can modulate muscle development and meat quality via altering the expression patterns of genes involved in crucial metabolic pathways. The findings presented here provide unique insights into the molecular mechanisms of muscle development in Pekin ducks in response to dietary C. butyricum.
Pekin duck is a popular poultry species worldwide because of its rapid growth and considerable breast muscle mass . Duck meat, in particular breast meat, is considered a delicacy worldwide, especially in China. Breast muscle development of Pekin ducks usually begins at 4.5 and ends at 7 weeks of age . However, to increase the economic gain, the market age of Pekin ducks can be advanced to 5 or 6 weeks by overfeeding . If the quality of the breast muscle has not fully developed by this age, huge economic losses may result. It is therefore crucial to explore new approaches to improve the duck breast muscle yield and meat quality in practical production.
Diet modulation using additives such as oil, antioxidants and probiotics can improve muscle production, meat quality, and fatty acid composition. Many studies have demonstrated that dietary addition of probiotics such as Lactobacillus plantarum and C. butyricum can improve growth performance and meat quality in animals kept for meat [4, 5]. C. butyricum, a Gram-positive anaerobe, can generate butyric acid as well as form spores to aid survival in the intestine of healthy organisms . Therefore, C. butyricum is a potential alternative to antibiotics to promote the growth rate and immune function of Cherry Valley ducks . C. butyricum has been extensively applied in the poultry industry, because it has been shown to enhance growth performance, meat quality, and muscular fatty acid composition in broilers [8, 9]. The production and quality of breast meat are of extreme importance for consumers. Moreover, the profitability of the duck industry is mainly dependent on increasing the proportion of main segments in the carcass, especially the breast meat . However, there has been insufficient research on the effects of C. butyricum in the breast muscle of Pekin ducks. Our previous work showed that appropriate supplementation (400 mg/kg) of C. butyricum can improve meat quality and lipid metabolism . However, the underlying molecular regulatory mechanisms are still unknown.
Systems biology approaches are rapidly developing, and high-throughput transcriptome sequencing has been applied to explore the differentially expressed candidate genes that may be involved in particular biological processes. Simultaneously, Illumina sequencing technology allows more reliable gene identification [12, 13]. The relatively short reads generated by Illumina sequencing can be effectively de novo assembled and used for gene exploration and comparison of gene expression patterns . Using a high-throughput RNA-Seq technique and secondary bioinformatics analysis, the differentially expressed genes (DEGs) in breast muscle were identified between Pekin ducks treated with C. butyricum-treated and a control group. This provided a basis for subsequent functional exploration of the DEGs in enhancing breast muscle development of Pekin ducks.
The aim of this study was to investigate the molecular mechanisms and the metabolic pathways involved in the regulation of the physiological changes induced by dietary C. butyricum which have the potential to enhance the efficiency of poultry production and meat quality, and to identify key genes involved in the promotion of duck muscle development of ducks. This work provides a theoretical basis for further study of the molecular mechanisms involved in Pekin duck muscle development. Elucidation of the underlying mechanisms could help us to seek novel approaches to optimize the use of probiotics in the poultry and livestock industries.
The present animal study protocol was approved by China Agricultural University Ethics Committee (Permit No. CAU20170505–3) and closely followed the recommendations of the Guidelines for Experimental Animals. Before tissue sampling, all animals were humanely euthanized with sodium pentobarbitone (30 mg/kg BW) anesthesia. All efforts were made to minimize the suffering of animals.
Experimental design, animals and management
Six hundred 1-day-old male Pekin ducks were obtained from a local commercial hatchery (Beijing Golden Star Duck Co., Ltd) and the experiment was carried out at the Experimental Center of China Agricultural University. The ducks were randomly allocated into two groups with five replicates of sixty ducks each, and raised in an air-conditioned room. The temperature was maintained at 32 to 34 °C from day 1 to 5, and then gradually reduced by 2 °C /week until it reached a final room temperature of between 23 and 25 °C. Ducks were allowed free access to feed and water with a 23 L: 1D lighting program during the experimental time. The groups were as follows: Control group, ducks were fed a corn-soybean basal diet; C. butyricum-treated group, ducks were given a basal diet supplemented with 400 mg/kg of C. butyricum (2.0 × 109 CFU/g). The basal diet was formulated to meet or exceed the nutrient requirements according to the National Research Council (NRC, 1994) (Additional file 1: Table S1). The probiotic strain C. butyricum (Batch No. 20170325003) was obtained from Beijing Shine Biology Technology Co., Ltd., China. Feed intake was recorded by each replicate at 42 d of age. Ducks were weighed after feed deprivation for 6 h to calculate average body weight (ABW), average daily gain (ADG), average daily feed intake (ADFI), and the feed conversion ratio (FCR) at 42 d of age.
At 42 d of age, twelve ducks per group were randomly selected and weighed. After dissection, the breast muscle was also weighed. Three muscle samples of each group were randomly selected and rapidly stored in liquid nitrogen for further analyses.
RNA extraction and library preparation
Total RNA of the breast muscle tissues was extracted from the six randomly selected ducks using RNAiso puls kit (Takara, Osaka, Japan) following the manufacturer’s guidelines. After extraction, the concentration and purity of RNA were checked by an NanoDrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE) at 260 nm; and the integrity of the total RNA was measured by 1.2% agarose gel electrophoresis and an Agilent 2100 Bioanalyzer system (Agilent Technologies, Santa Clara, CA). The high-quality RNA samples for library construction were selected according to 260/280 nm ratio (1.9–2.1), RNA concentration (≥ 500 ng/μL) and RNA integrity number (RIN ≥ 8.0), respectively (Additional file 1: Table S2). Samples that did not satisfy the above quality criteria were re-extracted or not used for subsequent analysis.
Sequencing libraries were prepared by the TruSeq RNA Sample Prep Kit (Illumina, San Diego, CA) following the manufacturer’s guidelines. Briefly, the mRNA was purified from the total RNA using poly-T oligo-attached magnetic beads (Invitrogen, USA) to pull down the poly-A mRNA and cleaved into short fragments of about 300 bases. These fragmented mRNAs were then used as templates to synthesize cDNA by reverse transcriptase and random primers. The PCR was performed to amplify cDNAs to construct the library. After PCR enrichment, the quantity and quality of cDNA were measured by a NanoDrop 2000 spectrophotometer and Agilent 2100 Bioanalyzer. According to the manufacturer’s suggestions, 10 moll/μl Tris buffer (10 mmol Tris-HCl, 0.1% Tween 20, pH 8.5) was used to normalize the six cDNA libraries after construction. The cDNA libraries were sequenced on the Illumina sequencing platform (HiSeq 4000).
Quality control and de novo transcriptome assembly
FastQC program (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to check the quality of raw sequence reads generated from the Illumina Hiseq 4000 platform before alignment. Low-quality reads below a threshold quality of 20; reads shorter than 50 bp in length as well as reads containing adapter sequences, ploy-N and the sequencing primer from the raw data were removed to produce clean reads . Meanwhile, Error %, Q30, GC-content % and sequence duplication level of the clean data were evaluated (Additional file 1: Table S3). All of the subsequent analyses and annotation were depended on high quality clean reads.
Raw reads were dealt with to filter redundant sequences, low-quality reads and artificial repeats to obtain clean reads by SeqPrep software (https://github.com/jstjohn/SeqPrep) and Sickle software (https://github.com/najoshi/sickle) . The first assembly was performed using clean data by SeqPrep software to generate a de novo dataset  . The contiguous nucleotide sequences (contigs) were further extended by combining the de novo dataset with all available duck mRNA sequences from Ensembl, and then, the sequences were assembled by Trinity software(http://trinityrnaseq.github.io) . At last only unique contigs ≥300 bp in length can be used for downstream study after assembly.
Annotation and classification of the de novo assembled transcriptome
BLASTX (https://blast.ncbi.nlm.nih.gov/Blast.cgi) search (E-value <1e− 5) was conducted to align the de novo assembled transcriptome against four databases (NR-Non-Redundant Protein Sequence Database, GO-Gene Ontology, KEGG-Kyoto Encyclopedia of Genes and Genomes, and COG-Clusters of Orthologous Groups) to obtain the annotation and classification information to predict the biological function.
GO annotation and KEGG pathway analysis of differentially expressed genes
The differentially expressed genes (DEGs) between C. butyricum-treated group and control group were identified using edgeR software (http://www.bioconductor.org/packages/2.12/bioc/html/edgeR.html). Genes with a P-value < 0.05 and fold-change ≥2 or ≤ 0.5 were considered significant. We used Goatools software(https://github.com/tanghaibao/GOatools)to evaluate DEGs between the control group and C. butyricum-treated group by GO analysis [19, 20]. DEGs at P-value < 0.05 and fold-change ≥2 or ≤ 0.5 were regarded as significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/) is a public database resource which was used to predict the potential metabolic pathways or signal transduction pathways of DEGs. KOBAS software 3.0 (http://kobas.cbi.pku.edu.cn/index.php) was used to evaluate the statistical enrichment of DEGs among KEGG pathways (Corrected P- value < 0.05).
qRT-PCR confirmation of differentially expressed genes
To confirm the reproducibility and repeatability of the differential expressed genes obtained from RNA-Seq data, we performed qRT-PCR to evaluate 10 randomly selected genes related to muscle development. Total RNA was extracted using TRIzol reagent (Invitrogen) and was reverse-transcribed using the TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China) according to the manufacturer’s guidelines. Specific primers of each gene were designed using primer premier 5.0 were synthesized by Sunbio Biotech Ltd. (Beijing, China) (Table 1). The Top Green qPCR SuperMix (TransGen Biotech, Beijing, China) was used to perform qRT-PCR procedure following the manufacturer’s instruction. The reaction mixtures were added in a 96-well plate at 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 55 °C for 34 s, and then followed by a melting curve using an ABI QuantStudio 7 Flex Sequence Detection System (Applied Biosystems, Foster City, CA). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene of Anas platyrhynchos was set as an endogenous reference gene. After each real-time experiment, dissociation curve analysis was conducted to make sure that only one product was amplified. All experiments were preformed in triplicates for each biological repeat.
Data are presented as mean ± S.D. Comparison between C. butyricum-treated group and control group of growth performance and meat quality data was performed using the univariate ANOVA procedure with SPSS19.0 software (SPSS Inc., Chicago, IL, USA). Relative fold changes of target gene expression was calculated by the 2-ΔΔCt method . The ΔCt value was confirmed by subtracting the target Ct of each sample from its GAPDH Ct value.
Effects of C. butyricum on performance and meat quality of Pekin ducks
During the experimental period (42 days), in comparison with the control group, increased ABW, average daily gain (ADG), feed conversion ratio (FCR), and percentage of breast muscle as well as decreased average daily feed intake (ADFI) were observed significantly (P < 0.05) in the C. butyricum-treated group (Table 2). Moreover, no significant differences were found for pH24h, the lightness and the yellowness of meat color, whereas increased pH45min, higher redness of meat color and lower drip loss and shear force of breast muscle in ducks were measured significantly (P < 0.05) between the C. butyricum-treated group and the control group (Additional file 1: Table S4).
De novo transcriptome assembly profiles evaluation
The summary information of the quality and alignment of transcriptome sequence of the six breast muscle samples were presented in Table 3. Approximately 6.0 Gb high-quality clean data per sample were obtained for further analysis after removing low-quality reads, cutting adapter, and filtering. The duck genome sequence had been released; however, the low quality of the assembled genome and annotations made the transcriptome analysis difficult. In this study, only approximately 55% of clean reads were mapped onto the reference duck genome (http://asia.ensembl.org/Anas_platyrhynchos/Info/Index?db=core) with varied alignments for our data (Additional file 1: Table S5). Thus, Trinity software was used for de novo assembly with the clean reads; more than 86.60% of clean reads per specimen can be mapped back to the reference transcript after alignment (Table 3, Additional file 1: Table S6). Almost 68.15–72.98% clean reads were aligned in a unique manner, while 15.53–18.46% as multiple-mapped reads. The CG% of each sample was all above 50% (Additional file 1: Table S6). A total of 73,740 unigenes sequences were de novo assembled with an N50 length of 2192 bp. After removing low-quality sequences and redundancies, the de novo transcript had 91,365 sequences with an N50 length of 2733 bp. The length of transcripts ranged from 201 to 61,746 bp. A comparison of the duck transcripts to the best candidate assembly which could be used as a reference was presented in Table 4. Our de novo assembled duck transcripts could be better performed assembly, increasing the mapping rate of over 86% per muscle sample. Thus, comparison among the three datasets revealed that our de novo assembled duck transcripts were better to be used to explore gene expression patterns instead of the reference duck genome.
Functional annotation and classification
Homology analysis of the obtained sequences after clustering was performed by BLASTX (https://blast.ncbi.nlm.nih.gov/Blast.cgi) against NR databases at NCBI. The distribution of top hit species were determined and the majority of the annotated sequences were consistent with known nucleotide sequences of bird species, with Anas platyrhynchos (4604) followed by Gallus gallus (1620) (Additional file 2: Fig. S1). Generally speaking, nearly all sequences of Pekin duck should be similar with Anas platyrhynchos, while the sequences information of Anas platyrhynchos is not complete in the NCBI.
A total of 73,740 genes were allocated to three primary GO categories (biological_process, cellular_component, and molecular_function), which were subsequent assigned to 65 functional terms. The total genes were mainly related to cellular process (7104), single-organism process (6452), metabolic process (5941), cell (6346), cell part (6345), binding (6083), and catalytic activity (3543) (Fig. 1, Additional file 2: Fig. S2). According to the KEGG analysis, the top 5 classifications of genes were: signal transduction, global and overview maps, immune system, endocrine system transport and catabolism, respectively (Fig. 2). The top 20 enriched pathways included “PI3K-Akt signaling pathway”, “Focal adhesion”, “Regulation of actin cytoskeleton”, “MAPK signaling pathway”, and “Protein processing in endoplasmic reticulum” which may be associated with muscle development (Additional file 2: Fig. S3). In addition to GO analysis, COG analysis was preformed to deeply identify the function of the assembled unigenes. Genes were annotated based on the Cluster of COG database (Fig. 3). The largest number of matched genes was associated with signal transduction mechanism which is corresponded with KEGG pathways annotation.
Analysis of differential expressed genes
Gene expression levels were defined using fragments per kilobase of transcript per million mapped reads (FPKM), which is proportional to the quantity of cDNA fragments derived from the gene transcripts. The lowest limit of gene expression level is set as 0.5 FPKM in at least one of the samples. Based on the limit, a total of 56,797 genes were confirmed to express in the breast muscle tissues. The correlation analysis according to the gene expression patterns revealed that the correlations of six samples were greater than 0.80 (Additional file 3).
A total of 1119 genes were found expressed differentially at a fold-change ≥2 or ≤ 0.5 at P- value < 0.05 in response to C. butyricum treatment. Cluster pattern analysis of DEGs between control group (A1, A2, and A3) and the C. butyricum-treated group (B1, B2, and B3) was listed in Fig. 4, indicating that DEGs of each three samples selected from the control group or C. butyricum-treated group presented the approximately similar expression pattern. Among these genes, 403 out of 26,733 genes were up-regulated significantly while 716 out of 30,064 genes were down-regulated significantly (Fig. 5). The fold changes induced by dietary supplementation with C. butyricum ranged from − 10.77 to 5.41 (Additional file 3). Compared to the control group, a total of 6 up-regulated genes and 6 down-regulated genes in C. butyricum-treated group which have orthologs in Anas platyrhynchos gene database and may be associated with muscle development were shown in Table 5.
GO annotations and KEGG pathway analysis of differential expressed genes
Gene Ontology (GO) enrichment analysis was done according to three categories including biological process, cellular component, and molecular function. Of these categories, most DEGs were enriched in the “biological process” category (Fig. 6, Additional file 2: Fig. S4, Additional file 4). Within the “biological process” category, “regulation of striated muscle contraction”, “regulation of muscle system process”, “regulation of actin filament-based process”, “muscle tissue morphogenesis”, and “triglyceride catabolic process” were the most primary subcategories. As for the “cellular component” category, most DEGs were assigned to “contractile fiber part”, and “myosin complex”. Within the “molecular function” category, the mainly enriched subcategories were “metal ion transmembrane transporter activity”, “cation channel activity”, “passive transmembrane transporter activity”, and “substrate-specific channel activity” (P < 0.05). The following genes participate in above GO terms: MYH6, SCN5A, TBX1, MYOZ2, AQP1, KCNC4, MYL3, CHRNA1, MGP, HSPB1, SCN4B, KCNQ5, ATP2A, CPS1, TPM3, TNNI1, FLNC, ASIC4, FBXO_25, and LAMB4 (Additional file 3). After that, the DEGs were annotated through KEGG to confirm enriched pathways. The significantly (Corrected P < 0.05) enriched pathways included “Tight junction”, “ECM-receptor interaction”, “MAPK signaling pathway”, “TNF signaling pathway”, “Nitrogen metabolism”, and “Alanine, aspartate and glutamate metabolism” in our study (Fig. 7, Additional file 5).
Confirmation of gene expression with quantitative RT-PCR
To confirm the differentially expressed genes of breast muscle between control and treatment groups obtained by RNA-Seq, ten differential expressed genes (MGP, FLNC, FBXO_25, LAMB4, MGAT3, MYL3, TPM3, PPP1R3B, HSPB1, and DGKD) were randomly selected and their expression patterns were quantified by qRT-PCR. MGP, FLNC, FBXO_25, LAMB4, and MGAT3 were up-regulated genes and MYL3, TPM3, PPP1R3B, HSPB1, and DGKD were down-regulated genes in high throughput RNA-Seq. The qRT-PCR results showed the similar down- or up-regulated trend in the expression of these genes (Fig. 8); therefore, the qRT-PCR expressions validate the findings by high throughput RNA-Seq.
There has recently been great interest in achieving higher growth performance and better meat quality in livestock . Increasing attention is being paid to the nutritional levels and health benefits of meat, especially Pekin duck meat . In this study, increased ABW and ADG and decreased FCR indicated that C. butyricum could feasibly be used to replace antibiotics to promote the growth of Pekin ducks. C. butyricum has beneficial effects on the promotion of growth, including enhanced immune function, improved meat quality, and balanced intestinal microflora [24,25,26]. The yield of breast muscle, the main source of Beijing roast duck, is very important . In the present study, the percentage of breast muscle increased significantly. This showed that long-term dietary supplementation of C. butyricum alters metabolic pathways in muscle, which leads to the increased breast meat yield and meat quality .
Muscle development is a complicated physiological process including cell regeneration, differentiation, and migration. The underlying molecular mechanisms have been elucidated by high-throughput transcriptome analysis; however, the reference duck genome is incomplete, influencing the accuracy and precision of further analysis. In one recent study, only approximately 50% of reads could be mapped to the reference duck genome , which is consistent with difficulties we experienced. Therefore de novo assembly of the duck transcriptome was performed. Based on KEGG and COG classifications, many unigenes were enriched in signal transduction, indicating the importance of signal regulation in muscle development. Meanwhile, many DEGs enriched in the regulation of muscle system processes, muscle tissue morphogenesis, regulation of actin filament-based process, myosin complex, triglyceride catabolic process, and signaling pathways related to signal transduction were identified, indicating that muscle development requires a great deal of protein synthesis and signal transduction . The DEGs identified by RNA-Seq were confirmed by qRT-PCR.
Myogenesis is an extremely complicated physiological process comprising myogenic progenitor proliferation, myoblast proliferation and differentiation, laeding to mature muscle . Myogenesis can be regulated by many genes and pathways. We observed the up-regulation of genes related to myogenesis in the breast muscle of Pekin ducks treated with dietary C. butyricum. MGP (Matrix Gla protein) expression was increased in the C. butyricum-treated group. MGP shows an up-regulated expression pattern early in muscle cell development, suggesting that this high expression could be associated with cell differentiation and proliferation [31, 32]. Filamins are a family of actin binding proteins that are specifically expressed in skeletal muscle and were first identified by their ability to crosslink actin filaments . Filamin C (FLNC), which is mainly expressed in muscle cells, interacts with myotilin at the Z-band edge  and is required for the maintenance of the structural integrity of skeletal muscles . The high expression level in breast muscle cells in the C. butyricum-treated group suggests that FLNC may be a potential regulator facilitating muscle development. FBXO_25 (F_box protein 25) also showed increased expression between the C. butyricum-treated and control group. This is a member of the MAPK pathway that regulates numerous biological processes including cell growth, differentiation and apoptosis in skeletal muscle development . FBXO_25 may facilitate muscle development in the C. butyricum-treated group by specifically activating the MAPK signaling pathway. The LAMB4 (laminin subunit beta 4) is differentially expressed during pig muscle development, and has been proposed to be a positive factor that can activate myofiber formation , hence the elevated expression of LAMB4 observed in the C. butyricum-treated group in our study is likely to contribute to the activation of myofiber formation in duck muscle. MGAT3 also showed increased expression in the C. butyricum-treated group. Triacylglycerol synthesis is a critical function in muscular physiological processes, and MGAT3 (mannosyl (beta-1,4-)-glycoprotein beta-1,4-N-acetylglucosaminyltransferase) can catalyze the synthesis of diacylglycerol from monoacylglycerol, a key step in improving meat quality and fatty acid composition . These up-regulated genes work together to promote the muscle development of C. butyricum-treated Pekin ducks.
The expression levels of some genes were downregulated in the C. butyricum-treated Pekin ducks. MYL3 (Myosin light chain 3) combines calcium ions, promotes muscle development and participates in the contraction of cardiac muscles. MYL3/2/11 and TNNI1/2/3 (troponin I type 1/2/3 (skeletal, slow)) are also candidate genes in chicken embryonic muscle development . However, MYL3 is not expressed during the myogenesis process in vitro or in mature skeletal muscles . Moreover, MYL3 was down-regulated in pork with high drip loss, which is the ability of the skeletal muscle to maintain water post-mortem . The down-regulated expression of MYL3 in the C. butyricum-treated group in our study confirmed these previous observations. TPM3 (tropomyosin 3) is associated with nemaline myopathy, which presents as typically type 1 fiber hypotrophy, and mutations in the TPM3 gene can induce congenital fiber-type disproportion, which is characterized by generalized muscle weakness [41, 42]. The down-regulated expression of TPM3 in the C. butyricum-treated group is therefore beneficial for muscle development. PPP1R3B (Protein phosphatase 1, regulatory (inhibitor) subunit 3B) is a promising functional candidate gene for muscle development or meat quality [43, 44]. Expression of the PPP1R3B gene is negatively associated with pH24h of muscle in pigs , and higher expression is seen in pork with low muscle pH in agreement with our results. This indicates that lower expression of PPP1R3B could contribute to the increased post-mortem muscle pH in the C. butyricum-treated Pekin ducks. DGKD (Diacylglycerol kinases delta) may be a common regulator of genes in co-expression networks affecting muscle and meat properties in pigs , and is also associated with muscle spasm in mice . The down-regulated expression of DGKD in the C. butyricum-treated group may therefore also be beneficial for muscle development. HSPB1(heat shock protein family B member 1) is not associated with myofibrils during normal muscle development . Its low expression level in the C. butyricum-treated group may be evidence for the increased antistress ability of the muscles in these ducks . In summary, the decreased expression levels of MYL3, TPM3, PPP1R3B, DGKD and HSPB1 in the C. butyricum-treated group indicate that Clostridium butyricum can modulate gene expression patterns to promote muscle development of Pekin ducks.
It showed that the ECM-receptor interaction, the MAPK signaling pathway and the TNF signaling pathway accounted for many of the DEGs enriched in KEGG pathway analysis. The increase in skeletal muscle yield is mainly caused by muscle hypertrophy . Components of the ECM-receptor interaction contribute to the formation of the muscle niche . ECM-receptor interactions are enriched in fat-tailed sheep breeds, and the level of fat deposition is a vital factor affecting meat quality. Interestingly, LAMB4, a gene related to the ECM-receptor interaction pathway, was also up-regulated in the tissues of fat-tailed sheep, concordantly with our results , suggesting that the ECM-receptor interaction may be a potent regulator of muscle cell proliferation and differentiation during skeletal muscle development. Recently, various studies have focused on the MAPK pathway in skeletal muscle development because this pathway directly affects muscle morphology by activating several myogenic transcription factors and cell differentiation [53, 54]. The DEGs between the C. butyricum-treated group and the control group were significantly enriched in the MAPK signaling pathway, an important regulator of skeletal muscle development. The MAPK signaling pathway can be activated during hypertrophy and facilitates increased muscle production through the increased muscle fiber size . Supplementation with C. butyricum may therefore activate the MAPK signaling pathway and ECM-receptor interactions in the breast muscle of Pekin ducks. Further investigations are needed to explore the roles of the candidate genes in the MAPK signaling pathway and the ECM-receptor interaction in breast muscle and in determining Pekin duck meat quality. Genes that were differentially expressed between control and C. butyricum-treated ducks were also related to the TNF signaling pathway. This pathway may be associated with the increased immune function of ducks in the C. butyricum-treated group [11, 56]. Interestingly, the KEGG pathway analysis also showed enrichment of “Alanine, aspartate and glutamate metabolism”. In our previous study, supplementation with C. butyricum has been proved to increase the concentrations of essential amino acids and flavor amino acids in the breast muscle of ducks . The enrichment of alanine, aspartate and glutamate metabolism in breast muscle provides supporting evidence for the increased meat quality of ducks in the C. butyricum-treated group.
To sum up, the transcriptome analysis revealed that the breast muscle yield and meat quality in C. butyricum-treated Pekin ducks were improved via altering the expression abundance of muscle system process genes (MGP, FLNC, FBXO_25, LAMB4, MGAT3, MYL3, TPM3, PPP1R3B, DGKD and HSPB1) and activating pathways associated with muscle development (mainly the ECM-receptor interaction and the MAPK signaling pathway). To obtain further information regarding the mechanism between C. butyricum and breast muscle yield, further investigation is demanded to optimize the effects of probiotic, especially in the livestock industry.
We characterized the transcriptome profiles of breast muscle tissue in Pekin ducks supplemented with dietary C. butyricum by de novo assembly RNA-Seq. These differentially expressed genes may be vital for understanding the molecular mechanisms of breast muscle development induced by dietary supplementation with C. butyricum. GO term enrichment of the DEGs showed that regulation of muscle system processes and muscle tissue morphogenesis was significantly enriched. KEGG pathway enrichment analysis revealed that the ECM-receptor interactions, and the MAPK signaling pathway, which can promote muscle development and improve meat quality, were activated after supplementation with C. butyricum in Pekin ducks. To the best of our knowledge, this is the first study on the relationship between breast muscle gene expression and C. butyricum supplementation in Pekin ducks. It provides novel insights into the responses to dietary supplementation with C. butyricum, and may contribute to the optimization of the use of probiotics as feed additives for breeding duck strains with high meat yields.
Acid sensing ion channel subunit family member 4
Sodium/potassium-transporting ATPase subunit alpha-2
Cholinergic receptor nicotinic alpha 1 subunit
Orthologous Groups of Proteins
Carbamoyl-phosphate synthase 1
Cysteine and glycine rich protein
- ECM-receptor interaction:
Extracellular matrix-receptor interaction
Potassium voltage-gated channel, Shaw-related subfamily, member 4
Potassium voltage-gated channel subfamily Q member 5
Mitogen-activated protein kinase
Myosin heavy chain, fast skeletal muscle
Myosin heavy chain 6
National Center for Biotechnology Information
National Research Council
Sodium voltage-gated channel beta subunit 4
Sodium voltage-gated channel alpha subunit 5
Tumor necrosis factor
Troponin I type 1 (skeletal, slow)
Yang QL, Lou XW, Wang Y, Pan DD, Sun YY, Cao JX. Effect of pH on the interaction of volatile compounds with the myofibrillar proteins of duck meat. Poult Sci. 2017;96(6):1963–9.
Xu T, Gu L, Schachtschneider KM, Liu X, Huang W, Xie M, Hou S. Identification of differentially expressed genes in breast muscle and skin fat of postnatal Pekin duck. PLoS One. 2014;9(9):e107574.
Chartrin P, Bernadet MD, Guy G, Mourot J, Hocquette JF, Rideau N, Duclos MJ, Baeza E. Does overfeeding enhance genotype effects on energy metabolism and lipid deposition in breast muscle of ducks? Comp Biochem Physiol A Mol Integr Physiol. 2006;145(4):413–8.
Suo C, Yin Y, Wang X, Lou X, Song D, Wang X, Gu Q. Effects of lactobacillus plantarum ZJ316 on pig growth and pork quality. BMC Vet Res. 2012;8:89.
Liao X, Wu R, Ma G, Zhao L, Zheng Z, Zhang R. Effects of Clostridium butyricum on antioxidant properties, meat quality and fatty acid composition of broiler birds. Lipids Health Dis. 2015;14:36.
Murayama TI, Mita N, Tanaka M, Kitajo T, Asano T, Mizuochi K, Kaneko KI. Effects of orally administered Clostridium butyricum MIYAIRI 588 on mucosal immunity in mice. Vet Immunol Immuno. 1995;48(3–4):333–42.
Zhuang S, Jiang FB, Jia ZX, Yan R. Clostridium butyricum can be used as a potential alternative for the antibiotic in Cherry Valley ducks. J Anim Plant Sci. 2015;25(5):1227–32.
Zhao X, Guo Y, Guo S, Tan J. Effects of Clostridium butyricum and Enterococcus faecium on growth performance, lipid metabolism, and cecal microbiota of broiler chickens. Appl Microbiol Biotechnol. 2013;97(14):6477–88.
Yang CM, Cao GT, Ferket PR, Liu TT, Zhou L, Zhang L, Xiao YP, Chen AG. Effects of probiotic, Clostridium butyricum, on growth performance, immune function, and cecal microflora in broiler chickens. Poult Sci. 2012;91(9):2121–9.
He J, Lu L, Tian Y, Tao Z, Wang D, Li J, Li G, Shen J, Fu Y, Niu D. Analysis of intramuscular fat and fatty acids of different duck breeds and their association with SNPs of duckA-FABPgene. Can J Anim Sci. 2011;91(4):593–6.
Liu YH, Li YY, Feng XC, Wang Z, Xia ZF. Dietary supplementation with Clostridium butyricum modulates serum lipid metabolism, meat quality, and the amino acid and fatty acid composition of Peking ducks. Poult Sci. 2018;97(9):3218–29.
Yi B, Chen L, Sa R, Zhong R, Xing H, Zhang H. High concentrations of atmospheric ammonia induce alterations of gene expression in the breast muscle of broilers (Gallus gallus) based on RNA-Seq. BMC Genomics. 2016;17(1):598.
Kong BW, Hudson N, Seo D, Lee S, Khatri B, Lassiter K, Cook D, Piekarski A, Dridi S, Anthony N, et al. RNA sequencing for global gene expression associated with muscle growth in a single male modern broiler line compared to a foundational barred Plymouth rock chicken line. BMC Genomics. 2017;18(1):82.
Moreton J, Dunham SP, Emes RD. A consensus approach to vertebrate de novo transcriptome assembly from RNA-seq data: assembly of the duck (Anas platyrhynchos) transcriptome. Front Genet. 2014;5(5):190.
Cock PJ, Fields CJ, Goto N, Heuer ML, Rice PM. The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids Res. 2010;38(6):1767–71.
Paszkiewicz K, Studholme DJ. De novo assembly of short sequence reads. Brief Bioinform. 2010;11(5):457–72.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12(1):323.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnol. 2011;29(7):644–52.
Gene Ontology C. The gene ontology project in 2008. Nucleic Acids Res. 2008;36:D440–4.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Cheng C, Liu Z, Zhou Y, Wei H, Zhang X, Xia M, Deng Z, Zou Y, Jiang S, Peng J. Effect of oregano essential oil supplementation. Meat Sci. 2017;133:103–9.
Chartrin P, Méteau K, Juin H, Bernadet MD, Guy G, Larzul C, Rémignon H, Mourot J, Duclos MJ, Baéza E. Effects of intramuscular fat levels on sensory characteristics of duck breast meat. Poult Sci. 2006;85(5):914–22.
Chen W, Wang JP, Yan L, Huang YQ. Evaluation of probiotics in diets with different nutrient densities on growth performance, blood characteristics, relative organ weight and breast meat characteristics in broilers. Bri Poult Sci. 2013;54(5):635–41.
Zhang B, Liu X, Li S, Wang Z. Effects of dietary lipids and Clostridium butyricum on chicken volatile flavour compounds. J Anim Res. 2017. https://doi.org/10.18805/ijar.v0iOF.6998.
Juan Z, Zhao-Ling S, Ming-Hua Z, Chun W, Hai-Xia W, Meng-Yun L, Jian-Qiong H, Yue-Jie Z, Xin S. Oral administration of Clostridium butyricum CGMCC0313-1 reduces ovalbumin-induced allergic airway inflammation in mice. Respirology. 2017;22(5):898–904.
Zhao X, Ding X, Yang Z, Shen Y, Zhang S, Shi S. Effects of Clostridium butyricum on breast muscle lipid metabolism of broilers. Ital J Anim Sci Published online 02 Apr. 2018. https://doi.org/10.1080/1828051X.2018.1453758.
Zhu F, Yuan JM, Zhang ZH, Hao JP, Yang YZ, Hu SQ, Yang FX, Qu LJ, Hou ZC. De novo transcriptome assembly and identification of genes associated with feed conversion ratio and breast muscle yield in domestic ducks. Anim Genet. 2015;46(6):636–45.
Zhang RP, Liu HH, Liu JY, Hu JW, Yan XP, Wang DM, Li L, Wang JW. Transcriptional profiling identifies location-specific and breed-specific differentially expressed genes in embryonic Myogenesis in Anas Platyrhynchos. PLoS One. 2015;10(12):e0143378.
Harfmann BD, Schroder EA, Esser KA. Circadian rhythms, the molecular clock, and skeletal muscle. J Bio Rhythms. 2015;30(2):84–94.
Correia E, Conceicao N, Cancela ML, Belo JA. Matrix Gla protein expression pattern in the early avian embryo. Int J Dev Biol. 2016;60(1–3):71–6.
Hur DJ, Raymond GV, Kahler SG, Riegert-Johnson DL, Cohen BA, Boyadjiev SA. A novel MGP mutation in a consanguineous family: review of the clinical and molecular characteristics of Keutel syndrome. Am J Med Genet A. 2005;135(1):36–40.
Feng Y, Chen MH, Moskowitz IP, Mendonza AM, Vidali L, Nakamura F, Kwiatkowski DJ, Walsh CA. Filamin a (FLNA) is required for cell-cell contact in vascular development and cardiac morphogenesis. Proc Nati Acad Sci USA. 2006;103(52):19836–41.
Raynaud F, Jond-Necand C, Marcilhac A, Furst D, Benyamin Y. Calpain 1-gamma filamin interaction in muscle cells: a possible in situ regulation by PKC-alpha. Int J Biochem Cell B. 2006;38(3):404–13.
Fujita M, Mitsuhashi H, Isogai S, Nakata T, Kawakami A, Nonaka I, Noguchi S, Hayashi YK, Nishino I, Kudo A. Filamin C plays an essential role in the maintenance of the structural integrity of cardiac and skeletal muscles, revealed by the medaka mutant zacro. Dev Biol. 2012;361(1):79–89.
Jeong HS, Jung ES, Sim YJ, Kim SJ, Jang JW, Hong KS, Lee WY, Chung HM, Park KT, Jung YS, et al. Fbxo25 controls Tbx5 and Nkx2-5 transcriptional activity to regulate cardiomyocyte development. Biochim Biophys Acta. 2015;1849(6):709–21.
Cao J, Cheng L, Shi Y. Catalytic properties of MGAT3, a putative triacylgycerol synthase. J Lipid Res. 2007;48(3):583–91.
Ouyang H, Wang Z, Chen X, Yu J, Li Z, Nie Q. Proteomic analysis of chicken skeletal muscle during embryonic development. Front Physiol. 2017;8:281.
Zhang C, Wang J, Wang G, Ji Z, Hou L, Liu Z, Chao T. Molecular cloning and mRNA expression analysis of sheep MYL3 and MYL4 genes. Gene. 2016;577(2):209–14.
Neuhoff C, Cinar MU, Tesfaye D, Tholen E, Schellander K, Looft C, Groe-Brinkhaus C. The Effect of Protein Abundance in Porcine Skeletal Muscle on the Expression of Water Holding Capacity. International Plant and Animal Genome XXII Conference, San Diego, CA, USA. January 10–15, 2014.
Lawlor MW, Dechene ET, Roumm E, Geggel AS, Moghadaszadeh B, Beggs AH. Mutations of tropomyosin 3 (TPM3) are common and associated with type 1 myofiber hypotrophy in congenital fiber type disproportion. Hum Mutat. 2010;31(2):176–83.
Clarke NF, Kolski H, Dye DE, Lim E, Smith RL, Patel R, Fahey MC, Bellance R, Romero NB, Johnson ES, et al. Mutations in TPM3 are a common cause of congenital fiber type disproportion. Ann Neurol. 2008;63(3):329–37.
Munro S, Cuthbertson DJ, Cunningham J, Sales M, Cohen PT. Human skeletal muscle expresses a glycogen-targeting subunit of PP1 that is identical to the insulin-sensitive glycogen-targeting subunit G(L) of liver. Diabetes. 2002;51(3):591–8.
Ponsuksili S, Murani E, Phatsara C, Jonas E, Walz C, Schwerin M, Schellander K, Wimmers K. Expression profiling of muscle reveals transcripts differentially expressed in muscle that affect water-holding capacity of pork. J Agric Food Chem. 2008;56(21):10311–7.
Cinar MU, Kayan A, Uddin MJ, Jonas E, Tesfaye D, Phatsara C, Ponsuksili S, Wimmers K, Tholen E, Looft C. Association and expression quantitative trait loci (eQTL) analysis of porcine AMBP, GC and PPP1R3B genes with meat quality traits. Mol Biol Rep. 2012;39(4):4809–21.
Ponsuksili S, Siengdee P, Du Y, Trakooljul N, Murani E, Schwerin M, Wimmers K. Identification of common regulators of genes in co-expression networks affecting muscle and meat properties. PLoS One. 2015;10(4):e0123678.
Leach NT, Sun Y, Michaud S, Zheng Y, Ligon KL, Ligon AH, Sander T, Korf BR, Lu W, Harris DJ, et al. Disruption of diacylglycerol kinase delta (DGKD) associated with seizures in humans and mice. Am J Hum Genet. 2007;80(4):792–9.
Middleton RC, Shelden EA. Small heat shock protein HSPB1 regulates growth of embryonic zebrafish craniofacial muscles. Exp Cell Res. 2013;319(6):860–74.
Ueda S, Kokaji Y, Simizu S, Honda K, Yoshino K, Kamisoyama H, Shirai Y, Yamanoue M. Chicken heat shock protein HSPB1 increases and interacts with αB-crystallin in aged skeletal muscle. Biosci Biotechnol Biochem. 2015;79(11):1867–75.
Yang L, Xu Z, Li H, Xiong Y, Bo Z. Differential transcriptional analysis between red and white skeletal muscle of Chinese Meishan pigs. Int J Biol Sci. 2010;6(4):350–60.
Jiang C, Shi P, Li S, Dong R, Tian J, Wei J, Luo S. Gene expression profiling of skeletal muscle of nursing piglets. Int J Biol Sci. 2010;6(7):627–38.
Li B, Qiao L, An L, Wang W, Liu J, Ren Y, Pan Y, Jing J, Liu W. Transcriptome analysis of adipose tissues from two fat-tailed sheep breeds reveals key genes involved in fat deposition. BMC Genomics. 2018;19(1):338.
Williamson D, Gallagher P, Harber M, Hollon C, Trappe S. Mitogen-activated protein kinase (MAPK) pathway activation: effects of age and acute exercise on human skeletal muscle. J Physio. 2010;547(3):977–87.
Akimoto T, Pohnert SC, Li P, Zhang M, Gumbs C, Rosenberq PB, Williams RS, Yan Z. Exercise stimulates Pgc-1α transcription in skeletal muscle through activation of the p38 MAPK pathway. J Biol Chem. 2005;280(20):19587–93.
Felicio AM, Cesar AM, Tavares Moura AS, Sartori JR, Coutinho LL. RNA-Seq and Pathway Analysis of Chicken Breast Muscle Revealed Important Genes Associated with Metabolism Under Heat Stress. San Diego: Plant Anim Genome XXIV Conference; 2016.
Chen G, Goeddel DV. TNF-R1 signaling: a beautiful pathway. Science. 2002;296(5573):1634–5.
We gratefully acknowledge Shanghai Majorbio Bio-pharm Technology Co., Ltd. for their technical support in performing the data analysis. We thank Catherine Perfect for editing the English text of a draft of this manuscript.
This research was supported by Beijing Innovation Consortium of Agriculture Research System (BAIC04–2018) and Beijing Natural Science Foundation (No. 6172016). Funders had no role in the study design, analysis, or writing of this article.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available at NCBI project PRJNA437203 (https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA437203) with accession number SRP134061, and SRA with accession number SRR6812634, SRR6812635, SRR6812636, SRR6812637, SRR6812638, and SRR6812639. Any reasonable requests are available from the corresponding author.
Ethics approval and consent to participate
The present animal study protocol was approved by China Agricultural University Ethics Committee (Permit No. CAU20170505–3) and closely followed the recommendations of the Guidelines for Experimental Animals. Before tissue sampling, all animals were humanely euthanized with sodium pentobarbitone (30 mg/kg BW) anesthesia.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Composition and nutrient level of diets (air-dry basis). Table S2. Quality control of each RNA sample for sequencing. Table S3. Evaluation of clean data of each sample in this study. Table S4. Effects of Clostridium butyricum on breast meat quality of Pekin ducks. Table S5. Mapping ratios of each sample to the reference duck genome in this study. Table S6. Mapping Statics and ratios of each sample to de novo assembled transcripts in this study. (DOCX 29 kb)
Figure S1. Species distribution of the top BLAST hits. Figure S2. Different levels of GO term function classification. Figure S3. The top 20 pathways according to numbers of unigene annotation. Figure S4. The numbers of up- or down-regulated genes annotated to each GO term function. (DOCX 1113 kb)
Detailed information of DEGs. (XLSX 135 kb)
Detailed information of GO enrichment of DEGs. (XLS 31 kb)
Detailed information of KEGG enrichment of DEGs. (XLSX 30 kb)
About this article
Cite this article
Liu, Y., Jia, Y., Liu, C. et al. RNA-Seq transcriptome analysis of breast muscle in Pekin ducks supplemented with the dietary probiotic Clostridium butyricum. BMC Genomics 19, 844 (2018). https://doi.org/10.1186/s12864-018-5261-1
- Breast muscle
- C. Butyricum
- De novo assembly
- Pekin duck