Identification and co-expression analysis of long noncoding RNAs and mRNAs involved in the deposition of intramuscular fat in Aohan fine-wool sheep

Abstract

Background

High-quality lamb meat is becoming increasingly popular as living standards improve and dietary patterns change. Currently, evaluations of the meat quality of livestock have revealed that the content of intramuscular fat (IMF) is lower in carcass fats, yet IMF has a critically important influence on the edibility and flavor of muscle meat [1]. Indeed, the quantity of IMF has become one of the most critical parameters of meat quality indicators, as it is considered to be positively related to meat quality and texture [2, 3]. When a certain amount of fat is deposited between the muscle bundles and muscle fibers, the marbled section of the meat has a high score, and the meat is fresh and juicy, which is often considered ideal [4, 5]. The selective deposition of fat can improve production efficiency and play a key role in improving meat quality. This practice is also a major focus and challenge of modern livestock breeding [6]. Therefore, ensuring the appropriate deposition of IMF in lean meat can enhance the future quality of sheep meat.

Studies have shown that intramuscular lipid deposition is affected by multiple genes and signaling pathways, such as the FAS, FAM134B, and HSL genes and the Wnt and AMPK signaling pathways [7–9]. Recently, long non-coding RNAs (lncRNAs) have received increased attention for their wide-ranging functions. LncRNAs refer to a class of non-coding RNAs longer than 200 nt in length [10]. Most lncRNAs have significant temporal and spatial expression specificity [11, 12] and have low sequence conservation among species [13–15]. LncRNAs can be divided into five types based on their positions relative to neighboring protein-coding genes: intronic lncRNAs, bidirectional lncRNAs, sense lncRNAs, intergenic lncRNAs and antisense lncRNAs [16].

LncRNAs can regulate various life activities of the body, including epigenetic regulation, transcriptional regulation and post-transcriptional regulation [17–19]. The most common regulation methods of lncRNAs include cis-regulation of the transcription of neighboring protein-coding genes and the trans-regulation of non-adjacent genes. In addition, lncRNAs can interact with miRNAs to affect the post-transcriptional translation of related mRNAs [20–22]. Studies have shown that lncRNAs can play direct or indirect roles in the process of lipid accumulation [23]. SRA (steroid receptor RNA activator) is one of the earliest discovered lncRNAs and plays an important role in lipid metabolism. SRA can bind to peroxisome proliferator-activated receptor gamma (PPAR γ) and enhance PPAR γ activity, thereby promoting the differentiation of pre-adipocytes [24]. A study of the expression levels of lncRNAs in the IMF of Jinhua and Landrace pigs revealed a total of 119 differentially expressed lncRNAs (DELs), six of which were involved in fat deposition and lipid metabolism-related pathways [25]. Furthermore, an analysis of transcriptome data from IMF in Inner Mongolia goats revealed that 1472 lncRNAs were involved in adipocyte growth regulation and morphological changes of adipocytes [26]. Another study has shown that lncRNAs can play a key regulatory role in fat deposition in sheep tails [27]. Overall, these findings demonstrate that lncRNAs can regulate lipid deposition through a variety of regulatory mechanisms. However, few studies have assessed the roles of lncRNAs in intramuscular lipid deposition in sheep.

Aohan fine-wool sheep (AFWS) is an important meat-hair, dual-purpose sheep breed in China that grows rapidly early in development. The elimination of male lambs for fat lamb production can increase both hair and meat gains as well as improve the overall benefits provided by fine wool sheep [28]. Exploring the developmental characteristics of IMF deposition and selecting candidate genes for AFWS provide references for future studies and applications in sheep breeding, improve the quality of mutton and accelerate the breeding process. The goal of our study was to systematically identify the profiles of differentially expressed mRNAs (DEMs) and DELs during intramuscular lipid deposition in sheep through high-throughput sequencing. We hoped that by studying the relationship between lncRNAs and lipid deposition, our findings would shed light on the mechanisms underlying selective muscle lipid deposition in sheep.

Results

Profiles of lncRNAs and mRNAs in sheep muscle

A total of six RNA expression profiles were generated in this experiment. The results are shown in Table 2. The average raw reading was 13.62 G (1G means 1*10⁹ base). After preprocessing the raw data, the average value of the filtered data obtained in each library was 12.82 G. The data obtained from the six expression profiles were relatively average with Q20 ≥ 99% and G/C contents ranging from 49 to 53%, indicating that the quality of the filtered data was reliable. The comparison rate between the filtered clean reads and the reference genome was greater than 88% in all six samples, indicating that the experiment was free of contamination and that the experimental results were robust.

Table 2

Statistical data derived from RNA Sequencing

	Mth-2-1	Mth-2-2	Mth-2-3	Mth-12–1	Mth-12–2	Mth-12–3
Raw reads	90,397,074	89,528,496	88,813,860	90,619,672	95,989,952	89,577,658
	(13.59G)	(14.40G)	(13.44G)	(13.56G)	(13.43G)	(13.32G)
Valid reads	85,098,620	84,675,502	81,473,556	85,252,798	91,240,846	85,201,776
	(12.79G)	(13.69G)	(12.78G)	(12.76G)	(12.70G)	(12.22G)
Valid Ratio(%)	94.14	94.58	91.74	94.08	95.05	95.11
Q20%	99.98	99.98	99.97	99.98	99.98	99.98
Q30%	98.02	97.98	97.83	98.03	98.15	98.17
GC content(%)	49.5	51	53	49	49	49
Mapped reads	76,713,303	74,525,753	69,966,338	77,129,925	83,264,100	77,731,307
	(90.15%)	(88.01%)	(85.88%)	(90.47%)	(91.26%)	(91.23%)
Expressed genes	24,537	24,229	24,295	24,605	24,563	24,077
Unique lncRNAs	6495	6492	6507	6532	6541	6425

Open in a separate window

An average of 24,384 expressed genes was identified in the six libraries, and a summary of the protein-coding genes identified is provided in Additional file 1 (Table S1). An average of 6499 unique lncRNAs was identified in the libraries. The information associated with all identified lncRNAs is shown in Additional file 2 (Table S2). We found that the number of reads was positively related to the length of the chromosome (Fig. 1a). Based on the locations of novel lncRNAs in the genome, we identified 525 antisense lncRNAs, 304 sense lncRNAs, 350 bidirectional lncRNAs, 1710 intronic lncRNAs and 4046 intergenic lncRNAs (Fig. ​(Fig.1b).1b). The sequence information for all identified lncRNAs is shown in Additional file 3 (Table S3).

Open in a separate window

Fig. 1

Characteristics of lncRNAs and mRNAs in the intramuscular fat of sheep. a: Density distribution of lncRNAs. b: Distribution of different types of lncRNAs. c: Length distribution of lncRNAs and mRNAs. d: Distribution of exon number for lncRNAs and mRNAs. e: Expression levels (log₁₀FPKM) and numbers of lncRNAs and mRNAs. f and g: Length distribution of ORFs of lncRNAs and coding genes

The structural characteristics and the expression levels of lncRNAs and mRNAs were different. The average length of lncRNAs was 868 nt, which was shorter than the average length of mRNAs (2131 nt) (Fig. ​(Fig.1c).1c). LncRNAs consisted of 1.7 exons on average, while mRNAs had 9.9 exons on average (Fig. ​(Fig.1d),1d), thus, lncRNAs had fewer exons than mRNAs. Meanwhile, lncRNAs had lower expression levels relative to mRNAs (Fig. ​(Fig.1e).1e). Moreover, the length of the open reading frame (ORF) of lncRNAs tended to be shorter than that of mRNAs (Fig. ​(Fig.1f,1f, g). Overall, lncRNAs were characterized by shorter lengths, fewer exons, lower expression levels and shorter ORF length distributions compared with mRNAs.

Identification of differentially expressed mRNAs and lncRNAs

A total of 199 DEMs were identified in muscle tissue (log₂ (fold change) ≥ 1 or log₂ (fold change) ≤ − 1 and q < 0.05). Of these differentially expressed genes (DEGs), 70 were up-regulated and 129 were down-regulated (Fig. 2a, c). A summary of DEGs is provided in Additional file 4 (Table S4A). We identified 61 lncRNAs that were differentially expressed, of which 25 lncRNAs were up-regulated and 36 lncRNAs were down-regulated (Fig. ​(Fig.2b,2b, d). Interestingly, 58 DELs were novel lncRNAs, we will focus on these novel lncRNAs in subsequent research. The list of DELs is provided in Additional file 4 (Table S4B). To illustrate the distribution of DEGs, we created clustering maps of top 100 DEMs and all the DELs (Fig. ​(Fig.2e,2e, f). Red indicates that the gene had a higher expression level, and blue indicates that the gene had a lower level of expression.

Open in a separate window

Fig. 2

The differentially expressed mRNAs and lncRNAs in the intramuscular fat of sheep. a and b: The number of up-regulated and down-regulated differentially expressed mRNAs and lncRNAs. The left red bars represent the number of genes up-regulated; the right blue bars represent the number of genes down-regulated. c and d: The volcano of expressed mRNAs and lncRNAs. The left blue points represent significantly decreased mRNAs and lncRNAs; gray points represent mRNAs and lncRNAs without significant changes. The right red points represent significantly increased mRNAs and lncRNAs. e and f: The clustering maps of top 100 differentially expressed mRNAs and lncRNAs. Red indicates that the gene had a higher expression level, and blue indicates that the gene had a lower level of expression

Among the DEMs, we found many genes that have been shown to be related to lipid deposition. IGF2, CAPN6, UCP2 and SOCS2 were highly expressed at Mth-2. IGF2 and CAPN6 were reported to affect the deposition of intramuscular fat and play an important role in meat efficiency [29, 30]. UCP2 and SOCS2 can regulate the proliferation and differentiation of preadipocytes [31, 32]. FOSB, SCD and CMYA5 were highly expressed at Mth-12. They have all been found to be potential candidate gene affecting meat quality [33–35]. Among the DELs, we found that MSTRG.13347.2, MSTRG.16157.3, MSTRG.11343.1, MSTRG.11343.4, MSTRG.10678.1 were highly expressed at Mth-2. We speculated that these novel lncRNAs might inhibit lipid deposition. Meanwhile, MSTRG.3004.3, MSTRG.21053.3, MSTRG.14887.1, MSTRG.790.1, MSTRG.10518.3 were highly expressed at Mth-12. We speculated that these novel lncRNAs might promote lipid deposition. However, the regulatory mechanisms underlying these lncRNAs require further study.

Enrichment analysis of differentially expressed mRNAs

GO functional enrichment analysis of DEGs revealed that these genes participated in a total of 346 significantly enriched functional classifications (P < 0.05), 235 of which were related to biological processes, 30 related to cellular components and 81 related to molecular functions (Additional file 5: Table S5A). The top 25 of biological processes, top 15 of cellular components, top 10 of molecular functions were shown in Fig. 3a. The most significantly enriched GO terms were: DNA binding (GO:0003677), extracellular region (GO:0005576), calcium ion binding (GO:0005509), extracellular space (GO:0005615), oxidation-reduction process (GO:0055114), oxidoreductase activity (GO:0016491).

Open in a separate window

Fig. 3

GO terms and pathways analysis of differentially expressed mRNAs and target genes. a: The top 25 of biological processes, top 15 of cellular components, top 10 of molecular functions of differentially expressed mRNAs. b: Top 20 pathways of differentially expressed mRNAs. The size of the points represents the number of significant differentially expressed mRNAs that matched to a GO term. The color of the points represents the significance of enrichment. c: The top 25 of biological processes, top 15 of cellular components, top 10 of molecular functions of differentially expressed target genes

In addition, results of the KEGG pathway analysis showed that these DEGs were involved in 126 biological pathways (Additional file 5: Table S5B), 18 pathways of which were significantly enriched, including arachidonic acid metabolism (ko00590), linoleic acid metabolism (ko00591), steroid hormone biosynthesis (ko00140), and retinol metabolism (ko00830), all of which were related to lipid metabolism. Moreover, the top 20 signaling pathways are shown in Fig. ​Fig.3b.3b. The results indicated that these pathways may have significantly contributed to the deposition of IMF.

Comprehensive analysis of candidate lncRNAs and mRNAs

To understand the potential function of novel lncRNAs, we performed cis-regulation and trans-regulation analyses on candidate lncRNAs. A total of 61 DELs regulated 49 DEMs, all the lncRNAs that acted on 49 mRNAs through trans-regulation (Additional file 6: Table S6). GO analysis of targets of lncRNAs revealed that these genes participated in a total of 422 GO terms, 180 of which were significantly enriched (P < 0.05) (Additional file 7: Table S7A). In GO annotation, these DEGs primarily played a role in biological processes. For example, glyceraldehyde-3-phosphate metabolic process (GO:0019682), positive regulation of phospholipid translocation (GO:0061092), glyceraldehyde-3-phosphate biosynthetic process (GO:0046166), cell-cell junction assembly (GO:0007043), muscle cell cellular homeostasis (GO:0046716) and actin filament polymerization (GO:0030041). The top 25 of biological processes, top 15 of cellular components, top 10 of molecular functions were shown in Fig. ​Fig.3c.3c. The KEGG pathway enrichment analysis of target genes revealed a total of 47 annotated pathways (Additional file 7: Table S7B), of these pathways, one was significantly enriched (P < 0.05), namely Tight junction (ko04530). Although some pathways were not significantly enriched, such as Wnt signaling pathway (ko04310), AMPK signaling pathway (ko04152), mTOR signaling pathway (ko04150), Cell adhesion molecules (CAMs)(ko04514) and MAPK signaling pathway (ko04010) and Jak-STAT signaling pathway, these pathways have been reported in the literature to play an important role in lipid deposition [36–41]. Overall, the significantly enriched pathways and GO terms involve 49 target genes. Of these, seven target genes are associated with lipid deposition. There were lncRNAs whose pearson correlation coefficients ≥0.9 or ≤ − 0.9 associated with these seven target genes (Table 3). Of these, seven have high expression levels: MSTRG.16157.3, MSTRG.21053.3, MSTRG.19941.2, MSTRG.2469.2, MSTRG.4051.3, MSTRG.21381.1, MSTRG.12864.1. A network describing the connections between the source genes and lncRNAs (whose pearson correlation coefficients ≥0.8 or ≤ − 0.8) was constructed (Fig. 4). However, the regulatory mechanisms underlying these lncRNAs require further study.

Table 3

Critical target mRNAs and candidate lncRNAs related to lipid deposition

Critical target mRNAs	Candidate lncRNAs
EXOC6	MSTRG.8047.1/MSTRG.3581.1/MSTRG.8215.2/MSTRG.2469.2
FZD4	MSTRG.21381.1/MSTRG.8573.1/MSTRG.4051.3/MSTRG.19788.2/MSTRG.1098.2/MSTRG.21622.1/MSTRG.2792.1/MSTRG.792.2/MSTRG.7213.2/MSTRG.8912.2
NCOA1	MSTRG.3592.2/MSTRG.21238.1
PAQR3	MSTRG.7951.1/MSTRG.792.1/MSTRG.21053.3/MSTRG.21238.1/MSTRG.2792.1/MSTRG.792.2
ULK1	MSTRG.10920.3/MSTRG.7213.2/MSTRG.20109.2/MSTRG.21622.1/MSTRG.1098.2/MSTRG.19788.2/MSTRG.16157.3
FHL1	MSTRG.6483.2/MSTRG.14888.1/MSTRG.6483.2/MSTRG.12097.3/MSTRG.12864.1/MSTRG.3013.3
TPI1	MSTRG.19941.2/MSTRG.21238.1/MSTRG.7951.1/MSTRG.792.1/MSTRG.21053.3

Open in a separate window

Open in a separate window

Fig. 4

LncRNA-mRNA co-expression network. The red points represent critical mRNAs involved in lipid deposition. The yellow points represent candidate lncRNAs, which can regulate target genes

Validation of lncRNA and mRNA expression by qRT-PCR

To validate the expression levels of DELs and DEMs, we randomly selected five DELs and five DEMs and detected their expression levels by qRT-PCR (Fig. 5a). The results of RNA-seq are shown in Fig. ​Fig.5b.5b. Comparison of the two sets of results above revealed consistent regulatory trends of genes detected by the two methods, indicating that the RNA-seq data were accurate.

Open in a separate window

Fig. 5

Five differentially expressed mRNAs and five differentially expressed lncRNAs, which were detected by qRT-PCR. a: Expression levels of genes by qRT-PCR. b: Expression levels of genes by RNA-seq

Discussion

IMF content increased gradually with growth, as significant differences were detected between Mth-2 and Mth-12. These findings were consistent with a previous study showing that the IMF content of sheep increased from 0 to 6 months but remained stable thereafter until 12 months of age [42]. Furthermore, these findings are consistent with the characteristics of muscle growth and the development of experimental sheep. The sheep switched to a fattening phase after weaning at 2 months. The weight of sheep increased rapidly between the ages of 4 to 6 months, after which weight gain stabilized. IMF is an important feature contributing to meat quality. Therefore, we selected the LTL samples at Mth-2 (less lipid deposition) and Mth-12 (more lipid deposition) for RNA-seq to provide a robust test of gene expression differences.

Overall, we identified a total of 26,247 genes and 6935 predicted novel lncRNAs in LTL samples of sheep by RNA-seq. Among these, 199 mRNAs (70 up-regulated and 129 down-regulated) and 61 lncRNAs (25 up-regulated and 36 down-regulated) were differentially expressed. We found many DEMs that have been shown to be related to lipid deposition. A marker-derived gene network reveals CAPN6 can regulate intramuscular fat deposition of beef cattle as modulators of carcass and meat quality traits [30]. Transcriptome analyses examining IMF content in the LTL in heavy Iberian Pigs identified FOSB as a candidate gene and other regulatory factors [33]. A study examining gene expression differences in metabolism and function between intramuscular and subcutaneous adipocytes in cattle found that SCD was highly expressed in adipocytes and closely associated with fat formation [34]. Similarly, we found many factors that affect the differentiation of pre-adipocytes, such as UCP2 and SOCS2. As the target gene of miR-132–3p, UCP2 can regulate the differentiation of sheep precursor fat cells [31]. SOCS2 can also act as a regulator of adipocyte size [32]. However, there are still many DEMs whose functions are unknown, and whether they are related to lipid deposition still needs further research.

To further characterize the mechanisms underlying DEGs, we performed GO and KEGG analysis of DEMs. In GO annotation, these DEGs primarily played a role in biological processes. These processes were closely related to fat formation and deposition, such as reverse cholesterol transport (GO: 0043691), positive regulation of cholesterol efflux (GO:0010875), glyceraldehyde-3-phosphate metabolic process (GO:0019682), positive regulation of phospholipid translocation (GO:0061092), glyceraldehyde-3-phosphate biosynthetic process (GO:0046166) and positive regulation of cholesterol esterification (GO:0010873) [43–48]. In addition, we found that many genes were enriched in biological processes, such as signal transmission, organ development, biosynthesis and cell proliferation, and these are also important processes in muscle development. KEGG pathway analysis revealed that the DEMs were significantly enriched in the immune system, inflammatory response and biological metabolism pathways, demonstrating that signal transmission between adipocytes and immune cells can greatly affect the function of adipose tissue [49]. This result was consistent with the fact that inflammatory cell infiltration has been documented to commonly occur in adipose tissue and stimulate the activation of the immune defense system [50]. We also found that many pathways related to lipid metabolism (Cholesterol metabolism and Arachidonic acid metabolism) were significantly enriched for many DEMs, such as HABP2, CST6, DLK1, PLA2G4E and FOSL1 that have been reported to participate in lipid metabolism [51–55]. Based on the GO and KEGG analysis, we obtained DEMs expression profiles that affected the IMF deposition of sheep.

In our study, 61 DELs participated in the regulation of target mRNAs, among them, 58 lncRNAs were novel lncRNAs and all of them acted on 49 target genes through trans-regulation. The functionality of lncRNA is reflected through the study of their target genes [56]. We performed GO and KEGG enrichment analysis on these target genes. We focused on the GO terms related to lipid deposition. These included terms under biological processes, such as glyceraldehyde-3-phosphate metabolic process (GO:0019682), positive regulation of phospholipid translocation (GO:0061092), glyceraldehyde-3-phosphate biosynthetic process (GO:0046166) and cell-cell junction assembly (GO:0007043), as well as molecular functions, such as Wnt-protein binding (GO:0017147), Atg1/ULK1 kinase complex (GO:1990316) and eukaryotic translation initiation factor 2B complex (GO:0005851). We found that the identified lncRNAs might be related to these GO terms. However, further research is required to identify the precise mechanisms.

The KEGG enrichment analysis revealed that target genes were significantly enriched in one pathway, Tight junction (ko04530). Tight junctions have been reported that it has a regulatory effect on the communication of various substances and conditions between cells [57]. Recently, two novel aspects were found that one was their involvement in signal transduction and the other was that tight junctions were considered to be a crucial component of innate immunity [58]. Autophagosome affects metabolism of various substances the by selectively degrading lysosomal targets [59]. Although some signaling pathways were not significantly enriched in our study, such as the “Wnt signaling pathway” “MAPK signaling pathway” and “AMPK signaling pathway”, they are critically important in the process of lipid deposition. Myriad studies have demonstrated the crucial role of canonical Wnt/β-catenin cascade in the development of organs physiological homeostasis and biological metabolism [60]. Wnt/β-catenin signaling pathway is an important developmental pathway that negatively regulates adipogenesis [61]. MAPK/ PI3K/Akt signaling pathway can improve glucose and lipid metabolism via regulation of the metabolic profiling [40]. The AMPK signaling pathway coordinates cell growth, autophagy and metabolism [62]. The above analysis and KEGG pathways might play important roles in lipid deposition and deserve further study.

To facilitate future studies of the mechanisms underlying lncRNAs, we constructed the lncRNA-mRNA co-expression network contained 41 novel lncRNAs and seven mRNAs based on analysis of critical target genes. The seven pairs of lncRNAs-mRNAs included MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1. Target genes of these lncRNAs have been reported to be involved in lipid deposition. FZD4 is highly expressed during fat production [63] and ULK1 participates in lipid metabolism [64]. PAQR3 has modulatory roles in obesity, energy metabolism, and leptin signaling [65]. TPI1 as a glycolytic enzyme, can catalyze the interconversion of dihydroxy-acetone phosphate and glyceraldehyde 3-phosphate in the glycolytic pathway [66]. It was proposed as a potential biomarker for IMF [67]. FHL1 regulates gene transcription, cell proliferation, metabolism and apoptosis, and plays a role in fat deposition [68]. Exoc6 participates in GLUT4 translocation in adipocytes thereby affecting glucose transport in fat and muscle cells [69]. Evidence has shown that the p160 co-activator family, consisting of NCOA1, NCOA2 and NCOA3, plays a critical role in adipogenesis and might have a concordant effect on lipid metabolism in mammals [70]. These results provide information for future studies examining how lncRNAs regulate IMF deposition in sheep. The specific regulatory mechanisms require further study and testing.

Conclusions

Our study systematically identified mRNA and lncRNA expression profiles during intramuscular lipid deposition in sheep. We obtained a total of 199 DEMs and 61 DELs and identified some important lncRNAs related to lipid deposition through GO and KEGG enrichment analysis. In addition, co-expression network analysis of lncRNAs and mRNAs involving 41 lncRNAs and seven mRNAs was conducted based on significant KEGG pathways. Seven pairs of lncRNA-mRNA, including MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1 were selected for further research. Our study provided a list of the lncRNAs and mRNAs related to intramuscular lipid deposition and laid a foundation for future research on the regulatory mechanisms of lncRNA on sheep muscle lipid deposition.

Methods

Supplementary Information

Acknowledgments

Abbreviations

Authors’ contributions

NL and JNH designed this study; FHH, RRZ, LRL, LLL, and QL participated in sample collection; FHH, RRZ, LLL and QL performed qRT-PCR validation; FHH, JL, NL, LRL and JNH analyzed the RNA-Seq data; FHH wrote the manuscript with contribution from JL, NL, LRL and JNH. All authors reviewed and approved the final manuscript.

Funding

The work was supported by the Earmarked Fund for Modern China Wool & Cashmere Technology Research System (CARS-39); the Project of Shandong Province Agricultural Variety Program (2019LZGC012); the National Natural Science Foundation of China (31402047), and the Project of Shandong Province Higher Educational Science and Technology Program(J18KA136). The funding body had no role in the design of the study, collection, analysis, interpretation of data, decision to publish or writing the manuscript.

Availability of data and materials

Additional data can be found in supplementary files. The RNA-Seq data was submitted to the SRA database under accession number SRR12247890.

Ethics approval and consent to participate

All the experimental procedures were in accordance with the Guidelines for Experimental Animals of the Ministry of Science and Technology (Beijing, China) and approved by the Experimental Animal Ethics Committee of Qingdao Agricultural University. The managements of laboratory animal has in keeping with 《Laboratory Animal-Requirements of Environment and Housing Facilities》(GB 14925–2001). The written informed consent to participate was obtained from the AFWS Stud Farm in Inner Mongolia Autonomous Region. All efforts were made to minimize suffering.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Footnotes

Contributor Information

Fuhui Han, Email: moc.qq@5446910292.

Jing Li, Email: moc.621@99ybabysys.

Ranran Zhao, Email: moc.qq@910838716.

Lirong Liu, Email: nc.cehac@gnoriluil.

Lanlan Li, Email: moc.qq@702330518.

Qian Li, Email: moc.qq@8789344041.

Jianning He, Email: moc.361@401gnixgnixeh.

Nan Liu, Email: moc.anis@uilnan.

References