Integrated transcriptomics and metabolomics study of embryonic breast muscle of Jiaji ducks

doi:10.21203/rs.3.rs-3444280/v1

Download PDF

Research Article

Integrated transcriptomics and metabolomics study of embryonic breast muscle of Jiaji ducks

https://doi.org/10.21203/rs.3.rs-3444280/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Because number of matured muscle fibers in poultry does not increase after birth, the meat yield is mainly determined during embryogenesis. We previously indicated breast muscle grew rapidly from 18th day after hatching (E18) to E27, and almost stopped from E27 to E34 of Jiaji ducks, while the mechanism is unclear. This study utilized RNA-seq to explore the related genes of muscle development and their relationship with small molecule metabolites at E18, E27 and E34 of Jiaji ducks. Several thousand differentially expressed genes (DEGs) were detected. DEGs expression profiles included 8 trend maps, among which trend 1 was opposite to and trend 6 was consistent with breast muscle development trend of Jiaji ducks. Through joint analysis between trend 1 of DEGs and trend 1 of differential metabolites (DEMs), protein digestion and absorption pathway stood out. The decrease of COL8A2 gene expression will lead to the decrease of arginine content, which will inhibit the development of breast muscle in embryonic Jiaji duck. Similarly, joint analysis between trend 6 of DEGs and trend 6 of DEMs indicated the increase of GAMT gene expression will cause the increase of proline content, so as to promote the development of breast muscle of Jiaji duck in embryonic period.. These results will be helpful for further understanding the mechanism of muscle yields of Jiaji ducks.

Jiaji duck

transcriptome

metabolome

breast muscle

Joint analysis

Poultry is one of the most important meat resources in the world, providing a variety of animal products for consumers, such as high-quality meat, eggs, eiderdown and so on. To meet consumer's demand for animal products, the breeding, feeding and slaughtering technologies of traditional meat poultry such as broilers and meat ducks have been continuously innovated to improve the production efficiency and the quality of products of poultry. As an excellent meat type poultry, Muscovy duck (Cairina moschata) is very popular because of its fast growth, delicious meat and higher intramuscular fat [1, 3].

Poultry skeletal muscle development is a complex process, affected by genetic, nutritional and other factors [4, 5]. Skeletal muscle is made up of thousands of muscle cells, which have different sizes and shapes and are arranged in a specific way on the bone. During skeletal muscle development, mesodermal progenitor cells are selected to form myoblasts, a process known as muscle development or muscle formation. Myoblasts further fuse to form multinucleated myotubes, and myotubes further differentiate to form muscle fibers with contraction characteristics [6]. Since the number of muscle fibers in most animals does not increase after birth but is determined at the embryonic stage, muscle production in adult livestock and poultry is determined during embryogenesis [7].

Transcriptome sequencing is a basic method to study the genes and functions related to muscle development in livestock and poultry, and has been widely used. Feng et al found key genes related to porcine skeletal muscle development through transcriptome analysis [8]. Li et al revealed the changes in gene expression between wild chickens and domestic chickens induced by domestication through transcriptome analysis of Gastrocnemius [9]. Hu et al integrated transcriptomics and non-targeted metabolomics of Peking ducks, revealing the potential mechanism of skeletal muscle development in Peking ducks [10]. However, the proliferation and development of Jiaji duck breast muscle and its relationship with small molecule metabolites have not been studied.

Our previous study found that breast muscle cells of Jiaji ducks were mainly proliferated and the breast muscle weight increased rapidly E18 to E27. From E27 to E34, breast muscle cells almost stopped proliferating, and breast muscle weight basically stagnated (unpublished data). At the same time, we also carried out metabolomics analysis of E18, E27 and E34, and obtained a large number of small molecule metabolites related to the development and proliferation of Jiaji duck breast muscle (unpublished data). However, the genetic basis of muscle development and proliferation in Jiaji duck embryos and its relationship with metabolic small molecules are still unclear. Therefore, RNA-seq technology was used to explore the gene expression patterns, possible functions and their relationship with small molecule metabolites in the breast muscle of Jiaji ducks at E18, E27 and E34. The results of this paper will provide data support for further exploring the mechanism of proliferation and development of breast muscle for Jiaji duck and its relationship with small molecule metabolites.

2.1 Selection of samples

The duck eggs used in this paper were obtained from Hainan Chuanwei Muscovy Duck Breeding Co., Ltd., and a total of 200 eggs were selected with 81 g ± 5g. All eggs were incubated at 37℃±0.5℃ temperature and 86–87% humidity under the same conditions. Eight eggs were selected at E18, E27 and E34, respectively. The embryos were removed from the eggs, the entire embryo was weighed, one side of breast muscle were removed and weighed, and the other side of breast muscle were also removed and were snap-frozen in liquid nitrogen and stored at − 80℃.

2.2 RNA extraction, library construction and sequencing

Total RNA was extracted from 8 samples at each embryonic stage using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's protocol. The purity and concentration of RNA were verified by NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA) and agarose gel electrophoresis. The integrity of the sample RNA was evaluated using Bioanalyzer 2100 (Agilent, CA, USA) to ensure that the RIN value was greater than 0.7. About 1 ug total RNA was extracted and purified by Dynabeads Oligo (dT) 25-61005 (Thermo Fisher, CA, USA), and then the purified RNA was fragmented into small fragments by Magnesium RNA Fragmentation Module (NEB, cat.e6150, USA). Subsequently, the fragmented RNA was transcribed into complementary cDNA using reverse transcriptase. Then, with the help of E.coli DNA polymerase I and RNase H. Double-strand synthesis was performed to convert these RNA-DNA complexes into DNA double-strand structures. At the same time, dUTP solution was introduced into the synthesized DNA double strands to complete the end of the double-stranded DNA to the flat end. In order to allow it to be connected to an adaptor with a terminal T base, we added an A base at each end of it. Subsequently, we used magnetic beads to screen and purify its fragment size. The double strands were digested by UDG enzyme and amplified by PCR to control the size of the library fragment in the range of 300bp ± 50bp. RNA sequencing was performed in Hangzhou Lianchuan Biotechnology Co., Ltd., using Novaseq6000 sequencing platform, and the sequencing strategy was PE150 (150bp double-end sequencing).

2.3 Bioinformatics analysis of RNA-seq

We used Fastp software (https://github.com/OpenGene/fastp) to remove reads containing junction contamination, low-quality bases, and default parameter uncertain bases. Then use fastp to verify the sequence quality. We used HISAT2 ( https://ccb.jhu.edu/software/hisat2) ) to align sequencing data to the Muscovy duck reference genome (KizCaiMos1.0) [11]. StringTie software (https://ccb.jhu.edu/software/hisat2) was used to assemble genes or transcripts and FPKM (FPKM=[total_exon_fragments/mapped_reads(millions)×exon_length(kB)]) method was used to quantify the basic expression of all genes. The differentially expressed genes (DEGs) were selected with fold change > 2 or fold change < 0.5 and with parametric F-test comparing nested linear models (p value < 0.05) by R package edgeR (https://bioconductor.org/packages/release/bioc/html/edgeR.html). Finally, GO and KEGG enrichment analysis of genes was performed using DAVID software (https://david.ncifcrf.gov/).

2.4 Trend analysis

We used all DEGs obtained in this study to build a database according to the data requirements in the Lianchuan biological cloud platform through Excel software. After the database was completed, the data was analyzed through the STEM [12] cloud tool in the cloud platform to obtain the trend map of DEGs.

2.5 conjoint analysis

We first select the trend subgraphs that are consistent with (transcriptome trend 6) or opposite to (transcriptome trend 1) the trend of breast muscle development or myogenic cell proliferation from the trend diagram. Secondly, the required DEGs and related metabolites are used to construct the database through Excel, and then the veen cloud tool in the Lianchuan biological cloud platform is used to process the data to obtain the pathways enriched by DEGs and metabolites. The correlation between DEGs and related metabolites was analyzed by using the correlation clustering marker heat map cloud tool.

3.1 Overall evaluation of RNA-Seq data

To identify potential candidate genes affecting duck skeletal muscle during embryonic development, RNA-Seq was used to detect the gene expression profiles of breast muscle of Jiaji duck at E18, E27 and E34. The statistical information of RNA-seq data was presented in Table 1. From Table 1, it can be seen that the sequencing of breast muscle tissue produced an average of about 42.29 million raw reads per sample, with an average base of 6.34G. After filtering, the number of valid reads ranged from 31.90 million to 49.12 million, with an average base of 6.10 G, and the percentage of valid reads was above 94.22%. The Q20 and Q30 ranged from 99.96–99.97% and 97.89–98.38%, respectively. The results of valid data alignment to the duck reference genome were shown in Table S1 and Figure S1. The results showed that the average ratio of reads mapped to the reference genome of E18, E27 and E34 samples were 68.19%, 65.88% and 63.15%, respectively. On average, about 80% of the reads in each sample were mapped to exons, and a small portion was mapped to gaps and introns.

Table 1

Statistical information of RNA-seq data
Sample	Raw Data		Valid Data		Valid Ratio(%)	Q20%	Q30%	GC content%
	Read	Base	Read	Base
E18_1	51,318,794	7.70G	49,120,376	7.37G	95.72	99.96	98.02	48
E18_2	44,207,906	6.63G	42,637,844	6.40G	96.45	99.96	98.02	49
E18_3	41,625,394	6.24G	40,153,908	6.02G	96.46	99.96	98.00	50
E18_4	49,032,896	7.35G	47,004,982	7.05G	95.86	99.96	98.11	50
E18_5	34,816,218	5.22G	33,078,276	4.96G	95.01	99.97	98.38	49.50
E18_6	49,254,862	7.39G	47,081,622	7.06G	95.59	99.96	97.92	49
E18_7	37,277,092	5.59G	35,727,488	5.36G	95.84	99.96	97.91	50
E18_8	39,538,166	5.93G	38,019,306	5.70G	96.16	99.96	98.01	50
E27_1	44,807,248	6.72G	42,896,632	6.43G	95.74	99.96	97.89	50
E27_2	36,014,390	5.40G	33,931,090	5.09G	94.22	99.96	98.20	50
E27_3	50,835,772	7.63G	48,909,148	7.34G	96.21	99.96	97.89	50
E27_4	44,834,822	6.73G	42,701,490	6.41G	95.24	99.96	97.89	48.50
E27_5	39,579,520	5.94G	38,009,302	5.70G	96.03	99.97	98.27	49.50
E27_6	34,551,988	5.18G	33,392,038	5.01G	96.64	99.96	98.26	50
E27_7	42,407,282	6.36G	40,967,798	6.15G	96.61	99.97	98.28	50
E27_8	32,932,344	4.94G	31,900,926	4.79G	96.87	99.97	98.14	51
E34_1	38,982,814	5.85G	37,819,382	5.67G	97.02	99.97	98.19	50
E34_2	37,518,296	5.63G	36,292,284	5.44G	96.73	99.97	98.19	51
E34_3	45,707,548	6.86G	43,882,080	6.58G	96.01	99.97	98.23	50.50
E34_4	46,292,640	6.94G	44,547,798	6.68G	96.23	99.96	97.99	50.50
E34_5	48,190,280	7.23G	46,550,614	6.98G	96.60	99.97	98.07	50
E34_6	43,540,714	6.53G	42,164,096	6.32G	96.84	99.96	98.04	51
E34_7	45,208,206	6.78G	43,728,714	6.56G	96.73	99.96	98.06	51
E34_8	36,413,306	5.46G	35,162,300	5.27G	96.56	99.96	97.94	51

3.2 DEGs analysis

A total of 9,855 genes were detected among the three periods of Jiaji duck embryo breast muscle tissue sequencing samples, accounting for 62% of the total number of Muscovy duck annotated genes^[11]. Through differential expression analysis, 2,374 DEGs were detected for E18 VS E27 with 1,591 up-regulated and 783 down-regulated (Table S2), 5,033 DEGs were detected for E18 VS E34 with 3,999 up-regulated and 1,034 down-regulated (Table S3), 2,867 DEGs were detected for E27 VS E34 with 2,214 up-regulated and 653 down-regulated (Table S4, Fig. 1A). Through the heat map, we can find that the gene expression levels of each sample group were different (Fig. 1B), which showed that there were significant differences between the samples at different stages.

For DEGs between E18 VS E27, E18 VS E34, and E27 VS E34, GO enrichment and KEGG pathway analysis was performed and the top 20 GO biological processes and KEGG pathways were selected to construct rich factor figures. For DEGs in E18 VS E27, the GO terms of cell adhesion, homophilic cell adhesion via plasma membrane adhesion molecules, calcium ion binding, and chondrocyte differentiation were significantly enriched (P < 0.05). Cell adhesion, extracellular matrix and chondrocyte differentiation were related to cell proliferation and differentiation (Fig. 2A). KEGG pathway analysis showed that DEGs were mainly enriched in ECM-receptor interaction, Calcium signaling pathway, Focal adhesion, MAPK signaling pathway, and Wnt signaling pathway (P < 0.05, Fig. 2B). ECM-receptor interaction, MAPK signaling pathway, Calcium signaling pathway, Apelin signaling pathway, Wnt signaling pathway were related to cell proliferation and differentiation.

For DEGs in E27 VS E34, the GO terms of microtubule binding, centrosome, and structural constituent of cytoskeleton were significantly enriched (P < 0.05). In these 20 biological processes, mitotic cytokinesis, cholesterol biosynthetic process, cell adhesion, cellular response to retinoic acid were associated with cell proliferation and differentiation (Fig. 2C). KEGG pathway analysis showed that DEGs were mainly enriched in Cell cycle, ECM − receptor interaction, Adrenergic signaling in cardiomyocytes, Wnt signaling pathway, Focal adhesion and other signaling pathways (P < 0.05, Fig. 2D). Among them, ECM − receptor interaction, Focal adhesion, Cardiac muscle contraction, Wnt signaling pathway, GnRH signaling pathway, Progesterone − mediated oocyte maturation were the shared pathways enriched by DEGs in E18 VS E27. Most of these pathways were related to cell proliferation and differentiation.

For DEGs in E18 VS E34, the GO terms of DNA replication, embryonic limb morphogenesis, protein domain specific binding, and cell surface were significantly enriched (P < 0.05, Fig. 2E). In the enriched biological processes, G2 / M transition of mitotic cell cycle, positive regulation of cell population proliferation, chloride channel activity were related to proliferation and differentiation. KEGG pathway analysis showed that DEGs were mainly enriched in Cell cycle, Focal adhesion, ECM − receptor interaction, and ErbB signaling pathway (P < 0.05, Fig. 2F). ECM − receptor interaction, ErbB signaling pathway and Hedgehog signaling pathway were associated with cell proliferation and differentiation.

3.3 The expression profile analysis of DEGs

Through our previous results, we learned that the weight of breast muscle in Jiaji duck embryo grew rapidly from E18 to E27, and the weight of breast muscle was basically unchanged after E27. Through the expression profile analysis of DEGs, we obtained 8 trend maps (Fig. 3A). The expression of DEGs in trend 1 decreased from E18 to E27, and tended to be stable from E27 to E34 (Fig. 3B), which was opposite to the actual development trend of breast muscle. The expression of DEGs in trend 6 increased from E18 to E27, and tended to be stable from E27 to E34 (Fig. 3C), which was consistent with the actual breast muscle development trend. Therefore, we selected these two expression trends for further analysis.

There were 585 DEGs in trend 1 (Table S5). GO analysis of DEGs in trend 1 showed that most of the DEGs were significantly enriched in regulation of transcription, DNA templated, sequence specific DNA binding (P < 0.05, Fig. 4A and Table S6). KEGG analysis showed that DEGs were significantly enriched in TGF-beta signaling pathway, protein digestion and absorption, and vascular smooth muscle contraction. Among them, TGF-beta signaling pathway was extremely significant (P < 0.01, Fig. 4B and Table S7), and this pathway was closely related to cell proliferation and differentiation, containing ID4, GDF5, and ID2, the expression of these three genes showed a decreasing trend during muscle development (Fig. 6).

There were 332 DEGs in trend 6 (Table S8). GO analysis of DEGs in trend 6 showed that most of the DEGs were significantly enriched in plasma membrane, extracellular matrix structural constituent, sarcolemma (P < 0.05, Fig. 4C and Table S9). KEGG analysis showed that the significantly enriched DEGs were related to MAPK signaling pathway, Focal adhesion, ECM − receptor interaction and other pathways (P < 0.05, Fig. 4D and Table S10).

3.4 Joint analysis of transcriptomics and metabolomics

According to our preliminary experiments, we found a number of differential metabolites (DEMs) among E18, E27 and E34. Of them, the content of DEMs decreased from E18 to E27 and that of DEMs almost unchanged from E27 to E34 (Metabolome trend 1 and Table 11), which was particularly similar to the trend 1 of the DEGs analyzed above. There were also DEMs with their contents from E18 to E27, and were stable from E27 to E34 (Metabolome trend 6 and Table 12), which was basically consistent with the expression profiles of DEGs in transcriptome trend 6. Therefore, we conducted a joint analysis of the KEGG pathways for DEGs in transcriptome trend 1 and trend 6 and DEMs in the metabolome trend 1 and trend 6. By cross-comparing metabolomics and transcriptomics data, DEMs and DEGs were identified at the molecular and biochemical levels.

Between DEGs and DEMs in trend 1, 16 pathways were found to be shared by DEGs and DEMs, mainly including protein digestion and absorption, purine metabolism, and arginine biosynthesis (Fig. 5A). Protein is an important part of muscle. Genes that control protein synthesis may play important roles for muscle proliferation and differentiation. To this end, we further analyzed the protein digestion and absorption pathway, and explored DEGs and DEMs enriched in this pathway. The metabolites enriched in this pathway were arginine and L-arginine, DEGs included COL17A1, COL8A2, COL16A1, ENSAPLG00000010244.

The correlation between the gene and the metabolites was shown in Fig. 5B. From Fig. 5B, we can see that there was a significant positive correlation between arginine and COL17A1, COL8A2, ENSAPLG00000010244. A more in-depth analysis of these genes found that COL8A2 has a tight relationship with cell proliferation [13].

In trend 6, there were 11 pathways shared by DEMs and DEGs, including Sphingolipid metabolism, arginine and proline metabolism, pyrimidine metabolism, glutathione metabolism, etc (Fig. 5C). Arginine and proline metabolism pathway played an important role in cell proliferation and differentiation [14, 15] and the metabolite enriched in this pathway was Proline. DEGs enriched in this pathway were GAMT, NOS1, SAT1, and AGMAT. From the correlation heatmap of genes and metabolites, a significant relationship between proline and GAMT, SAT1, AGMAT can be found (Fig. 5D).

Skeletal muscle development is mainly regulated by genes and influenced by endogenous metabolites. The identification of key genes and metabolites of skeletal muscle development in embryonic stage is of great useful in explaining the molecular mechanism of muscle development in meat ducks. If only through single-omics analysis, it can not fully explain the relevant issues, so multi-omics analysis has been widely used [16, 17]. In this study, the transcriptome of Jiaji duck breast muscle during embryonic period was analyzed. The key genes and related metabolites of Jiaji duck breast muscle development during embryonic period were explained in this study. The results of this paper combined with our previous results of metabolomics analysis could provide a basis for the further study in exploring the mechanism of breast muscle development.

Through previous experiments, we found that the embryonic breast muscle cells proliferated rapidly and the weight of the breast muscle increased from E18 to E27. After E27, the cell proliferation almost stopped, and the weight of breast muscle did not change. The E34 was basically developed and ready for shelling. Therefore, we chose E18, E27, E34 to explore the key genes and metabolites that may be the regulators for skeletal muscle development of Jiaji ducks.

In this study, found many DEGs in E18 VS E27, E27 VS E34, E18 VS E34 and conducted GO and KEGG pathway analysis for these DEGs. We obtained two pathways shared by E18 VS E27, E27 VS E34, E18 VS E34, ECM-receptor interaction and focal adhesion. ECM consists of a complex mixture of structural and functional macromolecules, mainly including collagen, fibronectin, and laminin [18]. Studies have found that cell differentiation is affected by the interaction between cells and ECM [21]. For focal adhesion, studies have shown the site where integrin and proteoglycan mediated adhesion connects with the actin cytoskeleton is called focal adhesion (FA), which is dynamic multi-protein complexes that connect ECM with the intracellular cytoskeleton [22, 23]. The formation and maturation of FA is a key procedure during myoblast differentiation. Thus, the GO and KEGG pathway analysis of DEGs of this study indicate that DEGs tightly correlated with the proliferation of breast muscle cell and the weight of breast muscle for Jiajia ducks.

In this study, DEGs in transcriptome trend 1 was contrast to the trend of breast development and myoblast cell proliferation. Thus, DEGs in transcriptome trend 1 may inhibit the development of embryonic breast muscle of Jiaji duck. KEGG analysis showed that TGF-beta signaling pathway was the extremely significant pathway, which was related to cell proliferation and differentiation [24]. DEGs in this pathway included NOG, NBL1, CHRD, ID4, GREM1, GDF5, BMP2, BMP6, BMP7, SMAD6, LTBP1, ENSAPLG00000018058, ID2. Among them, ID4, ID2, and GDF5 were tightly related to skeletal muscle development. Study shows that ID4 gene can inhibit cell differentiation [25]. In stem cells, ID4 prevents other transcription factors from binding to DNA, thereby inhibiting the initiation of cell differentiation. This allows stem cells to continue to self-renew and proliferate, thereby maintaining their undifferentiated state [25]. Melnikova et al found that ID2 over-expression inhibited the differentiation of Sol8 myoblasts [26]. Sullivan et al also showed that ID2 inhibits cell differentiation but promotes the proliferation of different types of cells [27]. GDF5 belongs to the TGF-β family and follows the same Smad-dependent and Smad-independent cell signaling pathways [28–31]. Moore et al found that the addition of GDF5 to the medium can lead to higher proliferation and ECM deposition of various cell types in a dose-dependent manner [29].

We also found the expression profiles of DEGs in transcriptome trend 6 was consistent with the trend of breast muscle development and breast muscle cell proliferation of Jiaji ducks. Therefore, we supposed that DEGs in transcriptome trend 6 may promote the development of breast muscle and myoblast proliferation. In this study, we found DEGs were extremely significantly enriched in the MAPK signaling pathway, a pathway related to cell proliferation [32]. DEGs in this pathway included LRIT2 and DOK5. DOK5 protein can bind to a variety of growth factor receptors, such as EGFR and DGFR [33], and affect the transmission of growth factors, thereby affecting cell proliferation and differentiation. Xu et al showed that over-expression of DOK5 promoted cell proliferation and osteogenesis, and activated the canonical Wnt/β-catenin signaling pathway [34]. LRIT2 with its expression level consistent with skeletal muscle development and myoblast proliferation of Jiaji duck may be a promotor for skeletal muscle development and myoblast proliferation. However, there is no report about the role of LRIT2 up to date. Thus, the further study should be performed in the future. In addition, focal adhesion pathway, a pathway with extremely significant enrichment by DEGs in transcriptome trend 6, is closely related to cell proliferation and differentiation [35]. EGF is enriched in this pathway. Study shows EGF encodes epidermal growth factor [36], which is a signaling molecule that participates in the regulation of biological processes such as proliferation, differentiation by binding and activating EGF receptors on the cell surface [37].

Through the joint analysis, we found 16 overlapped pathways between transcriptome trend 1 and metabolome trend 1. Among them, the protein digestion and absorption pathway has been reported to be closely related to cell proliferation and differentiation [38, 39]. Further analysis found arginine and COL8A2 were in the pathway. Arginine can be used as a synthetic substance and involve in the synthesis of protein and other cellular active substances [40]. Some studies have found that arginine regulates signal transduction of cell proliferation and differentiation by participating in mechanisms such as nitric oxide synthesis and release [41]. In addition, arginine is also involved in the synthesis of DNA and protein [42], which plays an important role in cell proliferation and differentiation. COL8A2 involves in the construction and maintenance of extracellular matrix, and has a regulatory effect on cell proliferation and migration [13]. COL8A2 is co-expressed with activating enhancer binding protein 2 (TFAP2) during corneal development, which has been reported to be related to the proliferation of CEC [43]. The co-expression of COL8A2 and TFAP2 plays an important role in the migration and proliferation of vascular smooth muscle [44]. Hwang et al found that the cell survival rate, cell proliferation rate, cyclin D1 expression and the number of cells in S phase decreased in cells transfected with siCOL8A2 [45], indicating COL8A2 has a inhibitory role for cell proliferation. In this study, the expression trend of COL8A2 (Fig. 6) was contrast to the development trend of breast muscle of Jiaji duck, indicating that it had an inhibitory effect on the proliferation of muscle cells in the embryonic stage. In addition, a correlation analysis showed that arginine was significantly positively correlated with COL8A2, suggesting that the expression level of COL8A2 has a tight linkage with the content of arginine. At the same time, this also shows that the decrease of COL8A2 gene expression will lead to the decrease of arginine content, thus inhibiting the development of chest muscle of embryonic Jiaji ducks.

Through the joint analysis transcriptome trend 6 and metabolome trend 6, the arginine and proline metabolism pathway stood out because it was closely related to cell proliferation and differentiation. This pathway contains proline, which plays an important role in cell proliferation and differentiation [46] and can increase the rate of protein synthesis in muscle [47]. Proline can synthesize molecules related to protein synthesis, signal transduction and cell growth through the catalysis of a series of enzymes. Studies have shown that the regulation of the enzyme activity in the proline metabolic pathway can affect cell proliferation and differentiation [48]. DEGs enriched in this pathway are GAMT, SAT1, and AGMAT. The enzyme encoded by GAMT gene plays an important role in muscle cells [49], which is involved in the synthesis of creatine in muscle. Creatine directly affect the myogenic process (formation of muscle tissue), by altering secretions of myokines [50]. Arginine is the raw material for the synthesis of creatine [51]. There is a conversion relationship between proline and arginine[52–54], so proline has a certain effect on the synthesis of creatine. According to the correlation analysis, there is a significant positive correlation between GAMT and proline, so we believe that GAMT and proline jointly promote the proliferation of muscle cells.

According to the above analysis, we can draw the following conclusions. The decrease of COL8A2 gene expression will lead to the decrease of arginine content, which will inhibit the development of breast muscle in embryonic Jiaji duck. The increase of GAMT gene expression will cause the increase of proline content, so as to promote the development of breast muscle of Jiaji duck in embryonic period. In addition, ID2, ID4 and GDF5 genes may play a role in inhibiting muscle development during embryonic breast muscle development of Jiaji duck. DOK5 and EGF genes may promote muscle development in this process.These results will be helpful for further understanding the mechanism of muscle yields of Jiaji ducks.

E18:18th day after hatching

DEGs:differentially expressed genes

DEMs:differential metabolites

FA:focal adhesion

TFAP2:activating enhancer binding protein 2

Author Contributions

Conceptualization, T.S.X. and L.H.G; methodology, C.J.H; validation, D.F.W; formal analysis, S.Q.H; writing—original draft preparation, L.H.G and J.L.C; visualization, G.R. and R.L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China General Program, grant number 32272864 and 31972553; Key Research and Development Programs of Hainan Province, grant number ZDYF2022XDNY154; Chinese Modern Technology System of Agricultural Industry, grant number CARS-42-50.

Data Availability Statement

All data generated or analysed during this study are included in this article and its additional files. All RNA-seq data analysed during this study has been deposited in NCBI database (https://submit.ncbi.nlm.nih.gov/subs/sra/) under accession number: SUB13849747.

Acknowledgments

We thank Lianchuan Biotechnology Co., Ltd for the technical support provided, and at the same time thank Hainan Chuanwei Muscovy Duck Breeding Co., Ltd for the experimental samples required to provide material support.

Ethics approval and consent to participate

The animal experiment in this study were approved by the Animal Care and Use Committee of the Chinese Academy of Tropical Agricultural Sciences (No. CATAS-20220101). All experimental procedures were performed in agreement with relevant guidelines. This study was carried out in compliance with the ARRIVE guidelines.

Conflicts of Interest

No conflict of interest.

Guo L, Wei C, Yi L, Yang W, Geng Z, Chen X. Transcriptional Insights into Key Genes and Pathways Underlying Muscovy Duck Subcutaneous Fat Deposition at Different Developmental Stages. Animals (Basel). 2021;11(7):2099. Published 2021 Jul 15. doi:10.3390/ani11072099.
Kokoszyński D, Wilkanowska A, Arpášová H, Hrnčár C. Comparison of some meat quality and liver characteristics in Muscovy and mule ducks. Arch Anim Breed. 2020 May 5;63(1):137-144. doi: 10.5194/aab-63-137-2020. PMID: 32494586; PMCID: PMC7254918.
Kierończyk B, Rawski M, Mikołajczak Z, Wachowiak R, Homska N, Józefiak D. A Preliminary Study of Chemically Preserved and High-Moisture Whole Maize (Zea mays L.) Usage in Pekin Duck Nutrition: Effect on Growth Performance and Selected Internal Organ Traits. Animals (Basel). 2021 Apr 4;11(4):1018. doi: 10.3390/ani11041018. PMID: 33916534; PMCID: PMC8065844.
Wen C, Yan W, Mai C, Duan Z, Zheng J, Sun C, Yang N. Joint contributions of the gut microbiota and host genetics to feed efficiency in chickens. Microbiome. 2021 Jun 1;9(1):126. doi: 10.1186/s40168-021-01040-x. PMID: 34074340; PMCID: PMC8171024.
Ge C, Luo X, Wu L, Lv Y, Hu Z, Yu D, Liu B. Plant essential oils improve growth performance by increasing antioxidative capacity, enhancing intestinal barrier function, and modulating gut microbiota in Muscovy ducks. Poult Sci. 2023 Aug;102(8):102813. doi: 10.1016/j.psj.2023.102813. Epub 2023 May 25. PMID: 37343349; PMCID: PMC10404791.
Buckingham M. Gene regulatory networks and cell lineages that underlie the formation of skeletal muscle. Proc Natl Acad Sci U S A. 2017 Jun 6;114(23):5830-5837. doi: 10.1073/pnas.1610605114. PMID: 28584083; PMCID: PMC5468682.
Gu LH, Xu TS, Huang W, Xie M, Shi WB, Sun SD, Hou SS. Developmental characteristics of pectoralis muscle in Pekin duck embryos. Genet Mol Res. 2013 Dec 13;12(4):6733-42. doi: 10.4238/2013.December.13.6. PMID: 24391014.
Feng L, Si J, Yue J, Zhao M, Qi W, Zhu S, Mo J, Wang L, Lan G, Liang J. The Landscape of Accessible Chromatin and Developmental Transcriptome Maps Reveal a Genetic Mechanism of Skeletal Muscle Development in Pigs. Int J Mol Sci. 2023 Mar 29;24(7):6413. doi: 10.3390/ijms24076413. PMID: 37047386; PMCID: PMC10094211.
Li Q, Wang N, Du Z, Hu X, Chen L, Fei J, Wang Y, Li N. Gastrocnemius transcriptome analysis reveals domestication induced gene expression changes between wild and domestic chickens. Genomics. 2012 Nov;100(5):314-9. doi: 10.1016/j.ygeno.2012.07.008. Epub 2012 Jul 21. PMID: 22824655.
Hu Z, Liu X. Integration of Transcriptomics and Non-Targeted Metabolomics Reveals the Underlying Mechanism of Skeletal Muscle Development in Duck during Embryonic Stage. Int J Mol Sci. 2023 Mar 8;24(6):5214. doi: 10.3390/ijms24065214. PMID: 36982289; PMCID: PMC10049352.
Xu MM, Gu LH, Lv WY, Duan SC, Li LW, Du Y, Lu LZ, Zeng T, Hou ZC, Ma ZS, Chen W, Adeola AC, Han JL, Xu TS, Dong Y, Zhang YP, Peng MS. Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility. Genomics. 2022 Nov;114(6):110518. doi: 10.1016/j.ygeno.2022.110518. Epub 2022 Nov 5. PMID: 36347326.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006 Apr 5;7:191. doi: 10.1186/1471-2105-7-191. PMID: 16597342; PMCID: PMC1456994.
Hwang JS, Ma DJ, Choi J, Shin YJ. COL8A2 Regulates the Fate of Corneal Endothelial Cells. Invest Ophthalmol Vis Sci. 2020 Sep 1;61(11):26. doi: 10.1167/iovs.61.11.26. PMID: 32931574; PMCID: PMC7500139.
Chen X, Guo Y, Jia G, Liu G, Zhao H, Huang Z. Arginine promotes skeletal muscle fiber type transformation from fast-twitch to slow-twitch via Sirt1/AMPK pathway. J Nutr Biochem. 2018 Nov;61:155-162. doi: 10.1016/j.jnutbio.2018.08.007. Epub 2018 Sep 1. PMID: 30236872.
Brunton JA, Baldwin MP, Hanna RA, Bertolo RF. Proline supplementation to parenteral nutrition results in greater rates of protein synthesis in the muscle, skin, and small intestine in neonatal Yucatan miniature piglets. J Nutr. 2012 Jun;142(6):1004-8. doi: 10.3945/jn.111.154534. Epub 2012 Apr 25. PMID: 22535763.
Qiao N, Yang Y, Liao J, Zhang H, Yang F, Ma F, Han Q, Yu W, Li Y, Hu L, Pan J, Hussain R, Tang Z. Metabolomics and transcriptomics indicated the molecular targets of copper to the pig kidney. Ecotoxicol Environ Saf. 2021 May 1;218:112284. doi: 10.1016/j.ecoenv.2021.112284. Epub ahead of print. PMID: 33945902.
Sun X, Tu K, Li L, Wu B, Wu L, Liu Z, Zhou L, Tian J, Yang A. Integrated transcriptome and metabolome analysis reveals molecular responses of the clams to acute hypoxia. Mar Environ Res. 2021 Jun;168:105317. doi: 10.1016/j.marenvres.2021.105317. Epub 2021 Mar 26. PMID: 33819872.
Wang J, Zhang Q, Li S, Chen Z, Tan J, Yao J, Duan D. Low molecular weight fucoidan alleviates diabetic nephropathy by binding fibronectin and inhibiting ECM-receptor interaction in human renal mesangial cells. Int J Biol Macromol. 2020 May 1;150:304-314. doi: 10.1016/j.ijbiomac.2020.02.087. Epub 2020 Feb 11. PMID: 32057847.
Lock JG, Wehrle-Haller B, Strömblad S. Cell-matrix adhesion complexes: master control machinery of cell migration. Semin Cancer Biol. 2008 Feb;18(1):65-76. doi: 10.1016/j.semcancer.2007.10.001. Epub 2007 Oct 6. PMID: 18023204.
Hu C, Zuo Q, Jin K, Zhao Z, Wu Y, Gao J, Wang C, Wang Y, Zhan W, Zhou J, Cheng F, Sun H, Niu Y, Zhang Y. Retinoic acid promotes formation of chicken (Gallus gallus) spermatogonial stem cells by regulating the ECM-receptor interaction signaling pathway. Gene. 2022 Apr 30;820:146227. doi: 10.1016/j.gene.2022.146227. Epub 2022 Feb 3. PMID: 35124150.
Berman AE, Kozlova NI, Morozevich GE. Integrins: structure and signaling. Biochemistry (Mosc). 2003 Dec;68(12):1284-99. doi: 10.1023/b:biry.0000011649.03634.74. PMID: 14756624.
Rustad KC, Wong VW, Gurtner GC. The role of focal adhesion complexes in fibroblast mechanotransduction during scar formation. Differentiation. 2013 Oct;86(3):87-91. doi: 10.1016/j.diff.2013.02.003. Epub 2013 Apr 25. PMID: 23623400.
Mishra YG, Manavathi B. Focal adhesion dynamics in cellular function and disease. Cell Signal. 2021 Sep;85:110046. doi: 10.1016/j.cellsig.2021.110046. Epub 2021 May 15. PMID: 34004332.
Morikawa M, Derynck R, Miyazono K. TGF-β and the TGF-β Family: Context-Dependent Roles in Cell and Tissue Physiology. Cold Spring Harb Perspect Biol. 2016 May 2;8(5):a021873. doi: 10.1101/cshperspect.a021873. PMID: 27141051; PMCID: PMC4852809.
Holliday H, Roden D, Junankar S, Wu SZ, Baker LA, Krisp C, Chan CL, McFarland A, Skhinas JN, Cox TR, Pal B, Huntington ND, Ormandy CJ, Carroll JS, Visvader J, Molloy MP, Swarbrick A. Inhibitor of Differentiation 4 (ID4) represses mammary myoepithelial differentiation via inhibition of HEB. iScience. 2021 Jan 20;24(2):102072. doi: 10.1016/j.isci.2021.102072. PMID: 33554073; PMCID: PMC7851187.
Melnikova IN, Bounpheng M, Schatteman GC, Gilliam D, Christy BA. Differential biological activities of mammalian Id proteins in muscle cells. Exp Cell Res. 1999 Feb 25;247(1):94-104. doi: 10.1006/excr.1998.4330. PMID: 10047451.
Sullivan JM, Havrda MC, Kettenbach AN, Paolella BR, Zhang Z, Gerber SA, Israel MA. Phosphorylation Regulates Id2 Degradation and Mediates the Proliferation of Neural Precursor Cells. Stem Cells. 2016 May;34(5):1321-31. doi: 10.1002/stem.2291. Epub 2016 Feb 1. PMID: 26756672; PMCID: PMC5024556.
Docheva D, Müller SA, Majewski M, Evans CH. Biologics for tendon repair. Adv Drug Deliv Rev. 2015 Apr;84:222-39. doi: 10.1016/j.addr.2014.11.015. Epub 2014 Nov 21. PMID: 25446135; PMCID: PMC4519231.
Moore YR, Dickinson DP, Wikesjö UM. Growth/differentiation factor-5: a candidate therapeutic agent for periodontal regeneration? A review of pre-clinical data. J Clin Periodontol. 2010 Mar;37(3):288-98. doi: 10.1111/j.1600-051X.2009.01527.x. Epub 2010 Jan 19. PMID: 20088982.
Saiga K, Furumatsu T, Yoshida A, Masuda S, Takihira S, Abe N, Ozaki T. Combined use of bFGF and GDF-5 enhances the healing of medial collateral ligament injury. Biochem Biophys Res Commun. 2010 Nov 12;402(2):329-34. doi: 10.1016/j.bbrc.2010.10.026. Epub 2010 Oct 19. PMID: 20937261.
Wang RN, Green J, Wang Z, Deng Y, Qiao M, Peabody M, Zhang Q, Ye J, Yan Z, Denduluri S, Idowu O, Li M, Shen C, Hu A, Haydon RC, Kang R, Mok J, Lee MJ, Luu HL, Shi LL. Bone Morphogenetic Protein (BMP) signaling in development and human diseases. Genes Dis. 2014 Sep;1(1):87-105. doi: 10.1016/j.gendis.2014.07.005. PMID: 25401122; PMCID: PMC4232216.
Li Z, Liu S, Cai Y. EGFR/MAPK signaling regulates the proliferation of Drosophila renal and nephric stem cells. J Genet Genomics. 2015 Jan 20;42(1):9-20. doi: 10.1016/j.jgg.2014.11.007. Epub 2014 Dec 9. PMID: 25619598.
Zhang Y, Yan Z, Farooq A, Liu X, Lu C, Zhou MM, He C. Molecular basis of distinct interactions between Dok1 PTB domain and tyrosine-phosphorylated EGF receptor. J Mol Biol. 2004 Oct 29;343(4):1147-55. doi: 10.1016/j.jmb.2004.08.072. PMID: 15476828.
Xu L, Wu J, Yu Y, Li H, Sun S, Zhang T, Wang M. Dok5 regulates proliferation and differentiation of osteoblast via canonical Wnt/β-catenin signaling. J Musculoskelet Neuronal Interact. 2022 Mar 1;22(1):113-122. PMID: 35234166; PMCID: PMC8919664.
Mishra YG, Manavathi B. Focal adhesion dynamics in cellular function and disease. Cell Signal. 2021 Sep;85:110046. doi: 10.1016/j.cellsig.2021.110046. Epub 2021 May 15. PMID: 34004332.
Modamio-Hoybjor S, Mencia A, Goodyear R, del Castillo I, Richardson G, Moreno F, Moreno-Pelayo MA. A mutation in CCDC50, a gene encoding an effector of epidermal growth factor-mediated cell signaling, causes progressive hearing loss. Am J Hum Genet. 2007 Jun;80(6):1076-89. doi: 10.1086/518311. Epub 2007 Apr 24. PMID: 17503326; PMCID: PMC1867096.
Janmaat ML, Giaccone G. Small-molecule epidermal growth factor receptor tyrosine kinase inhibitors. Oncologist. 2003;8(6):576-86. doi: 10.1634/theoncologist.8-6-576. PMID: 14657536.
Lane KT, Beese LS. Thematic review series: lipid posttranslational modifications. Structural biology of protein farnesyltransferase and geranylgeranyltransferase type I. J Lipid Res. 2006 Apr;47(4):681-99. doi: 10.1194/jlr.R600002-JLR200. Epub 2006 Feb 13. PMID: 16477080..
Singh M, Bacolla A, Chaudhary S, Hunt CR, Pandita S, Chauhan R, Gupta A, Tainer JA, Pandita TK. Histone Acetyltransferase MOF Orchestrates Outcomes at the Crossroad of Oncogenesis, DNA Damage Response, Proliferation, and Stem Cell Development. Mol Cell Biol. 2020 Aug 28;40(18):e00232-20. doi: 10.1128/MCB.00232-20. PMID: 32661120; PMCID: PMC7459263.
Li P, Yin YL, Li D, Kim SW, Wu G. Amino acids and immune function. Br J Nutr. 2007 Aug;98(2):237-52. doi: 10.1017/S000711450769936X. Epub 2007 Apr 3. PMID: 17403271.
Wu G, Meininger CJ. Regulation of nitric oxide synthesis by dietary factors. Annu Rev Nutr. 2002;22:61-86. doi: 10.1146/annurev.nutr.22.110901.145329. Epub 2002 Jan 4. PMID: 12055338.
Miljkovic M, Seguin A, Jia X, Cox JE, Catrow JL, Bergonia H, Phillips JD, Stephens WZ, Ward DM. Loss of the mitochondrial protein Abcb10 results in altered arginine metabolism in MEL and K562 cells and nutrient stress signaling through ATF4. J Biol Chem. 2023 Jun 1;299(7):104877. doi: 10.1016/j.jbc.2023.104877. Epub ahead of print. PMID: 37269954; PMCID: PMC10316008.
Hara S, Kawasaki S, Yoshihara M, et al. Transcription factor TFAP2B up-regulates human corneal endothelial cell-specific genes during corneal development and mainte-nance. JBiolChem. 2019;294:2460–2469.
Marx SO, Totary-Jain H, Marks AR. Vascular smooth muscle cell proliferation in restenosis. Circ Cardiovasc Interv. 2011 Feb 1;4(1):104-11. doi: 10.1161/CIRCINTERVENTIONS.110.957332. PMID: 21325199; PMCID: PMC3816546.
Hwang JS, Ma DJ, Choi J, Shin YJ. COL8A2 Regulates the Fate of Corneal Endothelial Cells. Invest Ophthalmol Vis Sci. 2020 Sep 1;61(11):26. doi: 10.1167/iovs.61.11.26. PMID: 32931574; PMCID: PMC7500139.
Saibi W, Feki K, Yacoubi I, Brini F. Bridging Between Proline Structure, Functions, Metabolism, and Involvement in Organism Physiology. Appl Biochem Biotechnol. 2015 Aug;176(8):2107-19. doi: 10.1007/s12010-015-1713-0. Epub 2015 Jun 23. PMID: 26100388.
Brunton JA, Baldwin MP, Hanna RA, Bertolo RF. Proline supplementation to parenteral nutrition results in greater rates of protein synthesis in the muscle, skin, and small intestine in neonatal Yucatan miniature piglets. J Nutr. 2012 Jun;142(6):1004-8. doi: 10.3945/jn.111.154534. Epub 2012 Apr 25. PMID: 22535763.
Pandhare J, Cooper SK, Phang JM. Proline oxidase, a proapoptotic gene, is induced by troglitazone: evidence for both peroxisome proliferator-activated receptor gamma-dependent and -independent mechanisms. J Biol Chem. 2006 Jan 27;281(4):2044-52. doi: 10.1074/jbc.M507867200. Epub 2005 Nov 21. PMID: 16303758.
Baker SA, Gajera CR, Wawro AM, Corces MR, Montine TJ. GATM and GAMT synthesize creatine locally throughout the mammalian body and within oligodendrocytes of the brain. Brain Res. 2021 Nov 1;1770:147627. doi: 10.1016/j.brainres.2021.147627. Epub 2021 Aug 19. PMID: 34418357.
Farshidfar F, Pinder MA, Myrie SB. Creatine Supplementation and Skeletal Muscle Metabolism for Building Muscle Mass- Review of the Potential Mechanisms of Action. Curr Protein Pept Sci. 2017;18(12):1273-1287. doi: 10.2174/1389203718666170606105108. PMID: 28595527.
Wyss M, Kaddurah-Daouk R. Creatine and creatinine metabolism. Physiol Rev. 2000 Jul;80(3):1107-213. doi: 10.1152/physrev.2000.80.3.1107. PMID: 10893433.
Brunton JA, Bertolo RF, Pencharz PB, Ball RO. Proline ameliorates arginine deficiency during enteral but not parenteral feeding in neonatal piglets. Am J Physiol. 1999 Aug;277(2):E223-31. doi: 10.1152/ajpendo.1999.277.2.E223. PMID: 10444416.
Tomlinson C, Rafii M, Ball RO, Pencharz PB. Arginine can be synthesized from enteral proline in healthy adult humans. J Nutr. 2011 Aug;141(8):1432-6. doi: 10.3945/jn.110.137224. Epub 2011 Jun 15. PMID: 21677074.
Urschel KL, Rafii M, Pencharz PB, Ball RO. A multitracer stable isotope quantification of the effects of arginine intake on whole body arginine metabolism in neonatal piglets. Am J Physiol Endocrinol Metab. 2007;293:E811–8.

No competing interests reported.

Download PDF

Editorial decision: Revision requested
05 Apr, 2024
Reviews received at journal
29 Feb, 2024
Reviews received at journal
25 Feb, 2024
Reviewers agreed at journal
19 Feb, 2024
Reviewers agreed at journal
19 Feb, 2024
Reviewers invited by journal
19 Feb, 2024
Editor assigned by journal
19 Feb, 2024
Editor invited by journal
18 Oct, 2023
Submission checks completed at journal
18 Oct, 2023
First submitted to journal
14 Oct, 2023

You are reading this latest preprint version

Integrated transcriptomics and metabolomics study of embryonic breast muscle of Jiaji ducks

Status:

Version 1

Abstract

Figures

1. Introduction

2. Materials and methods

2.1 Selection of samples

2.2 RNA extraction, library construction and sequencing

2.3 Bioinformatics analysis of RNA-seq

2.4 Trend analysis

2.5 conjoint analysis

3. Results

3.1 Overall evaluation of RNA-Seq data

3.2 DEGs analysis

3.3 The expression profile analysis of DEGs

3.4 Joint analysis of transcriptomics and metabolomics

4. Discussion

Conclusions

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1