Hao Zhou1,2
Yu Yang1Lixia Wang1Shengqiang Ye1Jiajia Liu2Ping Gong1Yunguo Qian1Huijun Zeng3,4*
Xing Chen1*- 1Insitute of Animal Husbandry and Veterinary, Wuhan Academy of Agricultural Science, Wuhan, China
- 2Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
- 3Wuhan Institute for Food and Cosmetic Control, Wuhan, China
- 4Key Laboratory of Edible Oil Quality and Safety for State Market Regulation, Wuhan, China
The Liancheng white (LW) duck is one of the most valued Chinese indigenous poultry breeds. Its meat is rich in nutrients and has distinct flavors, but the molecular mechanisms behind them are unknown. To address this issue, we measured and compared multi-omic data (genome, transcriptome, and metabolome) of breast meat from LW ducks and the Mianyang Shelduck (MS) ducks. We found that the LW duck has distinct breed-specific genetic features, including numerous mutant genes with differential expressions associated with amino acid metabolism and transport activities. The metabolome driven by genetic materials was also seen to differ between the two breeds. For example, several amino acids that are beneficial for human health, such as L-Arginine, L-Ornithine, and L-lysine, were found in considerably higher concentrations in LW muscle than in MS duck muscle (p < 0.05). SLC7A6, a mutant gene, was substantially upregulated in the LW group (p < 0.05), which may lead to excessive L-arginine and L-ornithine accumulation in LW duck meat through transport regulation. Further, guanosine monophosphate (GMP), an umami-tasting molecule, was considerably higher in LW muscle (p < 0.05), while L-Aspartic acid was significantly abundant in MS duck meat (p < 0.05), showing that the LW duck has a different umami formation. Overall, this study contributed to our understanding of the molecular mechanisms driving the enriched nutrients and distinct umami of LW duck meat, which will provide a useful reference for duck breeding.
1 Introduction
The Liancheng white duck (LW, Anas platyrhynchos domesticus), also called the “White Dove Duck,” is one of the precious and rare waterfowl resources in China. Historical sources note that LW ducks have been raised in Lian city for over 700 years. Because of the long process of artificial and natural selection, LW ducks have genetic stability and have evolved various characteristics, such as white feathers, black beaks, and black feet. Its meat is considered to be enriched in nutrients and has unique flavor (Bai et al., 2021; Zhang et al., 2021b; Huo et al., 2021). However, the underlying molecular mechanism remains unknown.
Recently, multi-omic association analysis has become popular in the investigation of potential molecular mechanisms in complex trait and disease studies (Theron et al., 2013; Douglas et al., 2018; Chen et al., 2020; Arora et al., 2021; Cardoso et al., 2021; Doran et al., 2021; Lu et al., 2021). It has been proven a powerful tool for studying muscle regulation systems that drive meat quality traits (Ye et al., 2017; Wang et al., 2020; Dou et al., 2022). Dou et al. (2022) used transcriptome and metabolome data to show that considerably high levels of muscle-useful metabolites, such as 13-HODE, NMN, and LPL, in black-boned chickens revealed essential health and therapeutic benefits. The integration of transcriptome and metabolome analysis in Bamei pork and Gansu Black pork also found essential genes and fatty acids that contribute to pork quality (Wang et al., 2020). Pork quality is also markedly related to vast quantitative trait loci (QTL) on the pig genome in the pig QTL database (Zhan et al., 2022), suggesting that the genetic factors play an important role in defining the characteristics of meat quality (Ye et al., 2017). However, the lack of genomic, transcriptomic, and metabolomic resources has limited the investigation into the molecular mechanisms of the meat characteristics of the LW duck.
In this study, we performed genome, transcriptome, and metabolome analyses to understand better the molecular mechanisms of the nutrients and meat flavor in LW duck meat. Another common indigenous shelduck, Mianyang Shelduck (MS, A. platyrhynchos domesticus) that has the same matrilineal origin from mallard (A. platyrhynchos) with LW duck (Li et al., 2010), was also measured using multi-omics and used as a control. Genetic variations, gene expressions, and metabolite contents in the breast meat of LW ducks were compared to those in MS duck meat. The biological gene–metabolite regulation was investigated in the breast meat of LW ducks. Overall, this study adds to our understanding of the molecular mechanisms behind the nutrition usage and meat flavor of the LW duck and will provide reference for duck breeding.
2 Materials and methods
2.1 Sample collection
LW ducks and MS ducks were reared to 300 days of age on the same farm, under the same standard management conditions. They were kept together from duckling to adulthood, fed a commercial corn–soybean-based diet ad libitum with free access to water and fasted for 8 h before being slaughtered. All ducks were healthy and had not been treated with antibiotics. In total, 16 female ducks (eight from each group) were randomly selected and fasted for 6 h. We took breast muscle samples and immediately stored them in liquid nitrogen for the following experiments.
2.2 Genome sequencing and analysis
Tissue from sixteen breast muscles (eight from each breed) was sampled for DNA extraction. DN easy tissue extraction kits (Qiagen, Valencia, CA, United States) were used to extract genomic DNA following the manufacturer’s protocol. Next-generation sequencing (NGS) technology was used to perform paired-end (PE) sequencing of the libraries based on the Illumina sequencing platform. The sequencing depth was 40 × for each sample. Reads were filtered and aligned to the reference genome of A. platyrhynchos domesticus (CAUduck1.0, http://uswest.ensembl.org/Anas_platyrhynchos_platyrhynchos/Info/Index) using the bwa software (Li and Durbin, 2009). Duplicate reads were removed from individual alignments using the Picard MarkDuplicates tool (v1.115, https://github.com/broadinstitute/picard). The Genome Analysis Toolkit (GATK) Haplotype Caller protocol was used to call single-nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) via the local re-assembly of haplotypes (Mckenna et al., 2010). The SNPs and INDELs were filtered prior to analysis with the GATK Variant Filtration protocol (Mckenna et al., 2010). The filter settings were as follows: QD < 10.0, ReadPos RankSum < − 8.0, FS > 10, QUAL <30, DP < 4. Then, Fst values of single SNPs and INDELs were calculated using vcftools v.0.1.13 (Danecek et al., 2011) to evaluate genetic differentiation. The range of FST is between 0 and 1. An FST value of 0 indicated that the mutations have no genetic divergence between populations and a value of 1 indicated the mutations have complete extreme division (Pino Del Carpio et al., 2011). SNPs and INDELs with the top 1% Fst values were taken to indicate selected mutations, which referred to previous comparative genomic studies (Islam et al., 2019; Liu et al., 2019; Xiong et al., 2019). They then were annotated by the SnpEff tool, which offered putative variant impacts (high, moderate, low, modifier) to categorize and prioritize variants on genes (Cingolani et al., 2012). A functional enrichment analysis of genes with mutations was performed using the Functional Annotation Chart tool in Database for Annotation, Visualization and Integrated Discovery (DAVID) (https://david-d.ncifcrf.gov/), using duck as background. p < 0.05 was set as the threshold for significantly enriched GO terms.
2.3 RNA-seq and analysis
Six breast muscle samples (three from each breed) were collected at random. We used the Trizol protocol for total RNA extraction of the samples. RNA sequencing was carried out on the Illumina HiSeq platform using protocols for 2 × 150 bp paired-end sequencing with the aim of obtaining 6G raw data from mRNA. After sequencing, raw data were trimmed and filtered using Trimmomatic 0.33 (Bolger et al., 2014). We followed the HISAT2, StringTie, and R package Ballgown protocol for gene expression quantification (Kim et al., 2015; Pertea et al., 2016). Paired-end reads were aligned to the reference sequence of A. platyrhynchos domesticus (CAUduck1.0) using HISAT2, and BAM files were sorted and indexed by SAMtools (v.1.3) (Li et al., 2009). Then, StringTie (v.1.3.1) was used to assemble transcripts, estimate transcript abundances, and create table counts for each sample. Finally, Ballgown (v.2.8.4) was used to extract gene-level expression measurements and call DEGs (|Log2FoldChange|≥0.5, p-value <0.05) from the Ballgown objects obtained from StringTie. The statistical test in ballgown is a parametric F-test comparing nested linear models (Frazee et al., 2015). A functional enrichment analysis of DEGs was performed using the Functional Annotation Chart tool in DAVID (https://david-d.ncifcrf.gov/), using duck as background. p < 0.05 was set as the threshold for screening significantly enriched GO terms. The DEGs were overlapped with the mutated genes, these overlapped genes were input to the Functional Annotation Table tool in DAVID for functional analysis.
2.4 Metabolite-seq and analysis
Metabolites in the breast muscle tissue of twelve ducks (six from each breed) were quantitatively profiled using an LC-MS/MS-based targeted metabolomics platform, as described previously (Bao et al., 2019). In brief, metabolites in muscle homogenates were extracted with methanol. The samples were incubated on ice for 5 min and then were centrifuged at 15,000 g, 4°C for 20 min. Some of the supernatant was diluted to a final concentration containing 53% methanol by LC-MS grade water. The samples were subsequently transferred to a fresh Eppendorf tube and then were centrifuged at 15,000 g, 4°C for 20 min. Finally, the 2 μL supernatant was injected into a SCIEX QTRAP® 6,500 + system (SCIEX, Framingham, MA, United States). The eluents were eluent A (0.1% Formic acid-water) and eluent B (0.1%Formic acid-acetonitrile). The solvent gradient was set as follows: 2% B, 2 min; 2–100% B, 15.0 min; 100% B, 17.0 min; 100–2% B, 17.1 min; 2% B, 20 min.
Metabolomics data analyses were carried out using the SCIEX OS 1.4 software (SCIEX, Framingham, United States), which performed the integration and correction of chromatographic peaks. The peak area of each chromatogram peak reflects the relative content of the related metabolites, and all the chromatographic peak area-integral data are finally saved. We corrected the detected mass spectrum peaks of each metabolite in the samples to ensure the accuracy of the qualitative and quantitative analyses through the quality control sample. An in-house novoDB (novogene database) was used to identify metabolites. The variable importance in the projection (VIP) score is used to estimate the important metabolites in a PLS-DA model. These metabolites were subsequently assessed using Student’s t test. Significantly differential metabolites (SDMs) were identified according to the following criteria: VIP value >1, fold change > 1.2 or fold change < 0.833, and p-value < 0.05. SDMs were then used for the KEGG pathway enrichment analysis, with significant criteria considered when p < 0.05. Metabolite–metabolite correlations were analyzed using Pearson’s correlation in the R software. Statistical significance was determined using the correlation t-test with p < 0.05. The metabolite correlation heatmap was drawn based on the correlation coefficients.
2.5 Joint analysis of transcriptomic and metabolomic data
To investigate the potential regulatory network of SDMs and DEGs, the expression profile matrix of DEGs and corresponding metabolic profile matrix of SDMs of six individuals (three from each breed) were prepared for correlation analysis. Correlation analyses were carried out using the Pearson method provided by R packages psych (v.1.8.12, https://CRAN.R-project.org/package=psych). Test analyses for correlation coefficients were performed by corr. test using the FDR method for multiple tests. A strong correlation between DEGs and SDMs was identified with the threshold |r|>0.8 and corrected p < 0.05, and then the significantly correlated combinations were selected and used to construct the network.
3 Results
3.1 Genetic characteristics of liancheng white duck
We performed whole-genome resequencing of 16 samples from LW ducks and MS ducks. In total, 635 Gb of raw reads was generated, with an average depth of 40 × (Supplementary Table S1A). These data were mapped against the duck reference assembly and called variants. In total, 4,008,770 INDELs and 12,332,104 SNPs were identified and used for the subsequent analysis.
To determine the genetic features of LW ducks, we explored the genomic landscape of SNP Fst values to identify candidate genes that potentially control specific traits present in two breeds. We identified 133,448 SNPs and INDELs with Fst values that were in the top 1% (Fst >0.78) (Figure 1) and thus were potentially positively selected in LW ducks (Supplementary Table S2A). Annotation of these SNPs indicated them to be located in 4,891 genes (Supplementary Table S2B). Further, 120 high effect mutations and 1,212 moderate effect mutations were located in 65 genes and 368 genes, respectively (Supplementary Table S2B). Indeed, several candidate genes with genetic variations have been shown to be associated with functions related to the meat traits of LW ducks (Supplementary Table S2C). One high effect mutation caused a frameshift in the sphingosine-1-phosphate lyase 1 (SGPL1) gene that was involved in the fatty acid metabolic process. We identified 12 mutations, including one nonsynonymous and eight synonymous, and three splice region variants on TMLHE involved in Lysine degradation and Trimethyllysine dioxygenase activity. In addition, one missense variant and 10 downstream gene variants were identified in GLUL that play roles in arginine biosynthesis and aspartate and glutamate metabolism. SLC7A6 associated with Arginine transport has one missense variant and one upstream gene variant. PLA2G6 and ARRDC3, fat and lipid metabolism-related genes, have 13 and 45 modifier effect variants, respectively, which may affect the corresponding gene expression. Further, several intron variants were found in SIRT7 and NAMPT. Overall, these gene variants may be responsible for the meat characteristics of LW ducks.
FIGURE 1. Photograph of the Liancheng white duck (A) and Mianyang shelduck (B). (C) Manhattan plots illustrate the Fst value between Liancheng white duck and Mianyang Shelduck. The red line represents the Top 1% Fst value threshold (Fst >0.78).
The 390 genes with high and moderate effect mutations are functionally enriched in several related metabolic processes, such as the beta-alanine biosynthetic process, lipid metabolic process, catabolic process, and galactose catabolic process via UDP-galactose (Supplementary Table S2C). The PLCXD3, PLA2G6, PLCH1, PAFAH2, GAL3ST2, DPYD, and GALT genes were observed to be involved in most of these functions.
3.2 Gene expression profiles of breast meat differ between liancheng white ducks and mianyang shelduck
Transcriptome sequencing was used to detect the mRNA expression profiles of the breast muscles of LW ducks and MS ducks. In total, approximately 254 million clean reads were obtained and about Q30 of the sequence reads per sample were about 95, suggesting a good sequence quality (Supplementary Table S1B).
The transcriptome differences were analyzed by comparing LW ducks and MS ducks. The comparison revealed 859 significantly differentially expressed genes (DEGs), including 415 upregulated and 444 downregulated genes in the LW ducks (p < 0.05). Of them, ARRDC3,CASTOR2,CTH, GLUL,SLC7A6, and SIRT7 were significantly upregulated in LW duck compared to MS ducks, while NAMPT, CHRND, and PDGFRB were significantly upregulated in MS duck compared to LW ducks. The details of all DEGs found in the compared breeds are shown in Supplementary Table S3A, and DEGs are shown on the volcano plot (Figure 2A).
FIGURE 2. (A)Volcano plot of differentially expressed genes (DEGs) between Liancheng white (LW) duck and Mianyang Shelduck. (B) Bubble plot for visualizing several GO annotations of DEGs. The horizontal axis refers to gene ratio (count of input genes/count of genes in the GO term). The Venn plots show the overlaps between downregulated DEGs in LW duck and mutated genes (C) and upregulated DEGs in LW duck and mutated genes (D).
The DAVID website was used to identify the DEG functions. Overall, 81 significantly enriched entries were identified in the biological process category (p < 0.05, Supplementary Table S3B). Most biological process-enriched items in the comparison were significantly associated with peptidyl-aspartic acid hydroxylation, the D-gluconate catabolic process, the methylglyoxal catabolic process to D-lactate via S-lactoyl-glutathione, the glutathione metabolic process, and the lipid metabolic process (Figure 2B). In molecule-enriched functions, gluconokinase activity and glutathione transferase activity were also enriched.
3.3 Integrative analysis of the genome and transcriptome
To investigate the DEGs genetically regulated, one hundred downregulated DEGs and ninety-eight upregulated DEGs in LW duck were identified and may be regulated by genetic mutations (Figures 2C,D). Among them, we investigated 9 upregulated DEGs in LW duck with genomic mutations that are related to amino acid metabolic and lipid-related functions (Table 1), which are strongly associated with the formation and quality of meat, as a previous report demonstrated (Zhan et al., 2022). Several significantly upregulated genes in the LW group, such as CTH, CASTOR2, GLUL, SIRT7, and SLC7A6, participate in amino acid metabolism and transport. GLUL also play roles in aspartate and glutamate metabolism. ARRDC3, PLA2G6, and SGPL1, three important genes controlling lipid and fatty acid metabolic processes, were observed to be significantly upregulated expressed (p < 0.05). The NAMPT expression was significantly upregulated in the MS group when compared to LW duck (p < 0.05), which was related to functions of nicotinamide metabolism.
TABLE 1. The upregulated differentially expressed genes with genetic variations in LW duck shown to have functions involved in amino acid metabolic and lipid-related functions (p < 0.05).
3.4 Metabolic profiles of breast meat differ between liancheng white duck and mianyang shelduck duck
In total, 585 metabolites were identified. Among them, Creatine, D-Glucose 6-phosphate, D-Mannose 6-phosphate, Taurine, and Anserine are the top five abundant metabolites, supporting results from other studies on ducks (Liu et al., 2013; Wang et al., 2020). PLS-DA analysis demonstrates that the samples from different groups were separated and classified into two distinct clusters (Figure 3A). The first PLS component (PC 1) accounted for 20.5% of the whole variation (Figure 3A). Screening yielded 83 significantly differential metabolites (SDMs), comprised of 16 and 64 metabolites upregulated and downregulated in the muscles of LW duck when compared to MS duck, respectively (VIP value > 1, fold change>1.2 or fold change< 0.833, and p < 0.05, Supplementary Table S4). L-Arginine, L-Ornithine, Trimethyllysine, L-Lysine, Guanosine monophosphate (GMP), and Xanthosine were significantly upregulated in the LW group when compared to MS duck (Figure 3C). Of them, L-Ornithine was the most significantly upregulated metabolite (Fold change = 2.92, p = 0.001) (Table 2). Whiles L-Aspartic acid and Nicotinamide were significantly upregulated in MS ducks when compared to LW duck (Figure 3C). Nicotinamide metabolism was proven regulated by the NAMPT gene (Zhan et al., 2022), which was significantly upregulated in the MS group compared to LW group and was reflected by the significantly changed compound nicotinamide in MS duck meat. In addition, most of SDMs exhibited a high correlation with other different metabolites, and most SDMs were negatively associated with L-Ornithine and Guanosine monophosphate (Figure 3B).
FIGURE 3. (A) PLS-DA Score Scatter plot of metabolites in Liancheng white duck and Mianyang Shelduck. (B) Heatmaps showing correlation patterns of top 20 significantly differential metabolites (SDMs) between LW duck and MS duck. (C) Box plot showing six SDMs. *Indicate p-value < 0.05; **indicate p value < 0.01. (D) Bubble plot for visualizing GO Annotation of SDMs.
TABLE 2. The upregulated significantly differential metabolites (SDMs) in LW duck compared to MS duck (VIP value > 1, fold change>1.2 or fold change< 0.833, and p < 0.05).
By annotating their functions and clustering the differential metabolites based on KEGG analysis, the results indicate that SDMs were significantly enriched in 40 pathways (Supplementary Table S4B, Figure 3D). Among them, “D-Arginine and D-ornithine metabolism” (map00472), and “Aminoacyl-tRNA biosynthesis” (map00970) was significantly enriched (p < 0.05). The levels of six differential metabolites in these pathways differed significantly (p < 0.05), including L-Ornithine, L-Arginine, L-Leucine, L-Lysine, L-Aspartic acid, and L-Cysteine.
3.5 Interactive network analysis of metabolomes and transcriptomes
To measure the associations between metabolomes and transcriptomes in the breast muscle of LW ducks and MS ducks, 859 DEGs and 83 SDMs were analyzed together in a correlation analysis. The focus on a strong correction showed 702 and 2,857 connections between DEGs and upregulated metabolites in LW ducks and MS ducks, respectively (r > 0.8, corrected p < 0.01, Pearson method) (Supplementary Table S5). The network of DEGs and upregulated metabolites in LW ducks was focused, and the results indicated that L-Ornithine was a key metabolite significantly positively correlated with the expression of multiple genes (Figure 4; Supplementary Table S5), including SLC7A6 (r = 0.98, corrected p = 0.0006), SIRT7 (r = 0.91, corrected p = 0.01), GLUL (r = 0.96, corrected p = 0.003), CTH (r = 0.91, corrected p = 0.01), PLA2G6 (r = 0.95, corrected p = 0.003), and CASTOR2 (r = 0.82, corrected p = 0.05). L-Arginine also was positively connected to the SLC7A6 (r = 0.88, corrected p = 0.02), SIRT7 (r = 0.87, corrected p = 0.02), and GLUL (r = 0.96, corrected p = 0.002). In addition, L-Lysine positively correlated with the TMLHE expression (r = 0.87, corrected p = 0.02). These correlations suggested the possible regulation relationship between these genes and metabolites.
FIGURE 4. A Network integrating gene–metabolite associations. Orange triangles indicate metabolites, green circles indicate genes. Important connects are highlighted in orange.
4 Discussion
The LW duck experienced a long-term domestication process, which resulted in the peculiar breed genetic traits discovered in this study, including many mutant genes with differential expressions involved in amino acid and fatty acid metabolism or transport activities. For example, one high effect mutation led to a frameshift in the SGPL1 gene, which was involved in the fatty acid metabolic process that was previously reported (Bektas et al., 2010). The metabolome data, potentially driven by genetic material, were also discovered to differ from LW ducks to MS ducks, which is consistent with the finding that the composition of free amino acids in animal meat, which contributes to meat scent and flavor, differs among breeds (Dou et al., 2022; Zhan et al., 2022). All these findings may be the root cause of the exceptional taste of LW duck meat.
Several amino acids that are beneficial for human health, such as L-Arginine, L-Ornithine, and L-Lysine, showed significantly higher levels in LW breast meat compared to MS breast meat. L-arginine is an essential amino acid in the fetus and neonate, and it is conditionally an essential nutrient for adults, particularly in certain disease conditions. L-arginine administration is beneficial in improving reproductive, cardiovascular (Popolo et al., 2014), pulmonary (Mitani et al., 1997), renal (Ketteler et al., 1994), gastrointestinal (Lojo et al., 2016), liver (Lojo et al., 2016), and immune functions (Geiger et al., 2016), as well as in facilitating wound healing (Wu et al., 2000). In muscle cells, L-arginine can be hydrolyzed to L-ornithine, which exerts anabolic or anticatabolic effects through their involvement in protein metabolism, immune response, and cell proliferation (De Bandt and Cynober, 1998; Durante et al., 2007). Further investigation of the underlying mechanism in the present study found that because the content of metabolites, including L-arginine and L-ornithine, was upregulated, the expression level of SLC7A6 was increased. SLC7A6 is an important subunit of the transport system y + L for L-arginine transport (Rotmann et al., 2007; Barilli et al., 2014). A high correction between the SLC7A6 gene and L-arginine and L-ornithine further confirms the regulation relationship between SLC7A6 and these metabolites. A mutation is present in the five upstream regions of the SLC7A6 gene and hence could cause the altered expression levels in LW meat. Taken together, these findings elucidate the roles and underlying mechanisms of SLC7A6 in the genetic regulation of L-arginine and L-ornithine, which may contribute to high L-arginine and L-ornithine accumulation in the muscle of LW ducks. L-canavanine, a structural analogue of L-arginine (Rosenthal, 1977), was also highly accumulated in the muscle of LW ducks. L-canavanine showed demonstrative antineoplastic activity against many human cancers (Ding et al., 1999; Nurcahyanti and Wink, 2016), increasing the nutrition of LW duck meat.
L-lysine is an amino acid required by animals. It is known as the first limiting amino acid, and it plays a vital role in human physiologic development and fatty acid oxidation (Jin et al., 2020; Matthews, 2020). Although lysine is a common amino acid in body proteins, it is scarce in many key dietary sources (e.g., grains or meat of animals fed an exclusively grain-based diet) (Matthews, 2020). In the current study, a greater L-lysine level was discovered in the meat of LW ducks. Fürst (1993) discovered that adding L-lysine to one’s diet might help cure osteoporosis, while Shimomura et al. (2014) found that modest dietary supplementation with L-lysine might reduce vascular calcification in uremic rats. In the current work, the TMLHE gene, which encodes Trimethyllysine dioxygenase and is involved in the L-lysine degradation pathway, was shown to be significantly downregulated in LW ducks, and it might be responsible for the higher lysine content in LW duck meat. In addition, 12 affected variants in the TMLHE gene were found. These mutations, as intrinsic genetic factors, may be responsible for modifying TMLHE functions or expressions, further impacting L-lysine accumulation in the meat of LW ducks.
The breast meat of LW ducks and MS ducks is enriched differentially in Umami-tasting compounds. GMP was significantly upregulated in LW duck meat, and it is used widely as an umami taste stimulus and potent flavor enhancer, as it synergistically increases the umami taste elicited by glutamate, another important umami-tasting compound (Kuninaka, 1960; Ninomiya, 2002; Lin et al., 2003). In the present study, there is no difference in glutamate abundance between LW and MS ducks, suggesting GMP may be the major driver of umami taste in LW duck meat. However, another ubiquitous umami compound (Ney, 1971; Zhang et al., 2021a), L-Aspartic acid was higher abundant in the breast meat of MS ducks than that of LW duck, suggesting the different umami formations of LW duck meat. Aspartic Acid plays an important role in the synthesis of other amino acids and in the urea cycles. It is also a nitrogen donor for arginine synthesis (Reitzer, 2004). Lower abundant L-Aspartic acid in LW duck (higher abundant in MS duck) maybe due to consumption in the synthesis of highly abundant arginine in LW duck.
In this study, the regulation relationship of mutated DEGs and SDMs in LW duck meat was identified based on an integrative analysis of multi-omic data. For example, the mutated genes SLC7A6 and TMLHE, correlated with several beneficial metabolites, such as L-Arginine, L-Ornithine, and L-Lysine, may play roles in the high accumulation of these beneficial metabolites in LW duck meat. Differential content of umami metabolites, such as GMP and L-Aspartic acid in LW duck meat, may be an important factor in the formation of the unique flavor of LW duck meat. Overall, these results provide effective information and more evidence to support further insight into the biological mechanisms of LW duck meat.
Data availability statement
The multi-omics data analyzed during this study are available in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra) with the accession codes PRJNA864643.
Ethics statement
The animal study was reviewed and approved by The protocols in the present study were approved by the Laboratory Animal Research (ILAR) guide for the Care and Use of Laboratory Animals at Wuhan Academy of Agricultural Science, China.
Author contributions
HZ, HZ, and XC are the principal investigators and project managers in this research; HZ, XC, YY, LW, SY, PG, and YQ. collected the sample and conducted data sequencing for the study; HZ JL, and XC contributed to data presentation; HZ, JL, and XC performed the sequencing data analysis; HZ and XC. wrote and edited the manuscript, with input from all authors.
Funding
This work was supported by the National Waterfowl-Industry Technology Research System (CARS-42), Key Research and Development Program of Hubei Province (2020BBA034), Foundation of Hubei Hongshan Laboratory (2021hskf012), (Knowledge Innovation Project of Wuhan Science and Technology Bureau (2022020801020417)) and Innovation Project of Wuhan Academy of Agricultural Sciences (CXJSFW202102-1).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.939585/full#supplementary-material
References
Arora, D., Srikanth, K., Lee, J., Lee, D., Park, N., Wy, S., et al. (2021). Integration of multi-omics approaches for functional characterization of muscle related selective sweep genes in Nanchukmacdon. Sci. Rep. 11, 7219. doi:10.1038/s41598-021-86683-4
Bai, H., Li, X., Zhong, L., Song, Q., Jiang, Y., Zhang, Y., et al. (2021). Study on depositions of mineral elements and expression levels of key genes in different tissues of Liancheng white ducks. Acta Agric. Zhejiangensis 33, 2264–2274. doi:10.3969/j.issn.1004-1524.2021.12.06
Bao, X., Wu, J., Kim, S., Lorusso, P., and Li, J. (2019). Pharmacometabolomics reveals irinotecan mechanism of action in cancer patients. J. Clin. Pharmacol. 59, 20–34. doi:10.1002/jcph.1275
Barilli, A., Rotoli, B. M., Visigalli, R., and Dall'asta, V. (2014). Gliadin activates arginase pathway in RAW264.7 cells and in human monocytes. Biochim. Biophys. Acta 1842, 1364–1371. doi:10.1016/j.bbadis.2014.04.021
Bektas, M., Allende, M. L., Lee, B. G., Chen, W., Amar, M. J., Remaley, A. T., et al. (2010). Sphingosine 1-phosphate lyase deficiency disrupts lipid homeostasis in liver. J. Biol. Chem. 285, 10880–10889. doi:10.1074/jbc.M109.081489
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170
Cardoso, T. F., Coutinho, L. L., Bruscadin, J. J., Da Silva Diniz, W. J., Petrini, J., Andrade, B. G. N., et al. (2021). Multi-omics approach reveals miR-SNPs affecting muscle fatty acids profile in nelore cattle. Genes 12, 67. doi:10.3390/genes12010067
Chen, X., Cai, W., Xia, J., Yu, H., Wang, Q., Pang, F., et al. (2020). Metabolomic and transcriptomic analyses reveal that blue light promotes chlorogenic acid synthesis in strawberry. J. Agric. Food Chem. 68, 12485–12492. doi:10.1021/acs.jafc.0c05020
Cingolani, P., Platts, A., Wang Le, L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. (Austin) 6, 80–92. doi:10.4161/fly.19695
Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., Depristo, M. A., et al. (2011). The variant call format and VCFtools. Bioinformatics 27, 2156–2158. doi:10.1093/bioinformatics/btr330
De Bandt, J. P., and Cynober, L. A. (1998). Amino acids with anabolic properties. Curr. Opin. Clin. Nutr. Metab. Care 1, 263–272. doi:10.1097/00075197-199805000-00005
Ding, Y., Matsukawa, Y., Ohtanifujita, N., Kato, D., Dao, S., Fujii, T., et al. (1999). Growth inhibition of A549 human lung adenocarcinoma cells by L-canavanine is associated with p21/WAF1 induction. Jpn. J. Cancer Res. 90, 69–74. doi:10.1111/j.1349-7006.1999.tb00667.x
Doran, S., Arif, M., Lam, S., Bayraktar, A., Turkez, H., Uhlen, M., et al. (2021). Multi-omics approaches for revealing the complexity of cardiovascular disease. Brief. Bioinform. 22, bbab061. doi:10.1093/bib/bbab061
Dou, T., Yan, S., Liu, L., Wang, K., Jian, Z., Xu, Z., et al. (2022). Integrative analysis of transcriptomics and metabolomics to reveal the melanogenesis pathway of muscle and related meat characters in Wuliangshan black-boned chickens. BMC Genomics 23, 173. doi:10.1186/s12864-022-08388-w
Douglas, G. M., Hansen, R., Jones, C. M. A., Dunn, K. A., Comeau, A. M., Bielawski, J. P., et al. (2018). Multi-omics differentially classify disease state and treatment outcome in pediatric Crohn's disease. Microbiome 6, 13. doi:10.1186/s40168-018-0398-3
Durante, W., Johnson, F. K., and Johnson, R. A. (2007). Arginase: A critical regulator of nitric oxide synthesis and vascular function. Clin. Exp. Pharmacol. Physiol. 34, 906–911. doi:10.1111/j.1440-1681.2007.04638.x
Frazee, A. C., Pertea, G., Jaffe, A. E., Langmead, B., Salzberg, S. L., Leek, J. T., et al. (2015). Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 33, 243–246. doi:10.1038/nbt.3172
Fürst, P. (1993). Dietary L-lysine supplementation: A promising nutritional tool in the prophylaxis and treatment of osteoporosis. Nutrition 9, 71–72.
Geiger, R., Rieckmann, J. C., Wolf, T., Basso, C., Feng, Y., Fuhrer, T., et al. (2016). L-arginine modulates T cell metabolism and enhances survival and anti-tumor activity. Cell 167, 829–842. e813. doi:10.1016/j.cell.2016.09.031
Huo, W., Weng, K., Gu, T., Zhang, Y., Zhang, Y., Chen, G., et al. (2021). Effect of muscle fiber characteristics on meat quality in fast- and slow-growing ducks. Poult. Sci. 100, 101264. doi:10.1016/j.psj.2021.101264
Islam, R., Li, Y., Liu, X., Berihulay, H., Abied, A., Gebreselassie, G., et al. (2019). Genome-wide runs of homozygosity, effective population size, and detection of positive selection signatures in six Chinese goat breeds. Genes (Basel) 10, E938. doi:10.3390/genes10110938
Jin, C.-L., Zhang, Z.-M., Song, Z.-W., Gao, C.-Q., Yan, H.-C., Wang, X.-Q., et al. (2020). mTORC1-Mediated satellite cell differentiation is required for lysine-induced skeletal muscle growth. J. Agric. Food Chem. 68, 4884–4892. doi:10.1021/acs.jafc.0c01275
Ketteler, M., Border, W. A., and Noble, N. A. (1994). Cytokines and L-arginine in renal injury and repair. Am. J. Physiol. 267, F197–F207. doi:10.1152/ajprenal.1994.267.2.F197
Kim, D., Langmead, B., and Salzberg, S. L. (2015). Hisat: A fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kuninaka, A. (1960). Studies on taste of ribonucleic acid derivatives. J. Agric. Chem. Soc. Jpn. 34, 489–492. doi:10.1271/nogeikagaku1924.34.6_489
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Li, H. F., Zhu, W. Q., Song, W. T., Shu, J. T., Han, W., Chen, K. W., et al. (2010). Origin and genetic diversity of Chinese domestic ducks. Mol. Phylogenet. Evol. 57, 634–640. doi:10.1016/j.ympev.2010.07.011
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi:10.1093/bioinformatics/btp352
Lin, W., Ogura, T., and Kinnamon, S. C. (2003). Responses to di-sodium guanosine 5'-monophosphate and monosodium L-glutamate in taste receptor cells of rat fungiform papillae. J. Neurophysiol. 89, 1434–1439. doi:10.1152/jn.00994.2002
Liu, C., Pan, D., Ye, Y., and Cao, J. (2013). ¹H NMR and multivariate data analysis of the relationship between the age and quality of duck meat. Food Chem. 141, 1281–1286. doi:10.1016/j.foodchem.2013.03.102
Liu, L., Cui, H., Xing, S., Zhao, G., and Wen, J. (2019). Effect of divergent selection for intramuscular fat content on muscle lipid metabolism in chickens. Animals. 10, E4. doi:10.3390/ani10010004
Lojo, N., Rasic, Z., Zenko Sever, A., Kolenc, D., Vukusic, D., Drmic, D., et al. (2016). Effects of diclofenac, L-NAME, L-arginine, and pentadecapeptide BPC 157 on gastrointestinal, liver, and brain lesions, failed anastomosis, and intestinal adaptation deterioration in 24 hour-short-bowel rats. PLoS One 11, e0162590. doi:10.1371/journal.pone.0162590
Lu, S., Wang, J., Zhuge, Y., Zhang, M., Liu, C., Jia, H., et al. (2021). Integrative analyses of metabolomes and transcriptomes provide insights into flavonoid variation in grape berries. J. Agric. Food Chem. 69, 12354–12367. doi:10.1021/acs.jafc.1c02703
Matthews, D. E. (2020). Review of lysine metabolism with a focus on humans. J. Nutr. 150, 2548S–2555S. doi:10.1093/jn/nxaa224
Mckenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Mitani, Y., Maruyama, K., and Sakurai, M. (1997). Prolonged administration of L-arginine ameliorates chronic pulmonary hypertension and pulmonary vascular remodeling in rats. Circulation 96, 689–697. doi:10.1161/01.cir.96.2.689
Ney, K. (1971). Flavour enhancing effect of l-glutamate and similar compounds. Z. Leb. Unters. Forch. 146, 141–143. doi:10.1007/bf01882239
Nurcahyanti, A. D., and Wink, M. (2016). L-Canavanine potentiates the cytotoxicity of doxorubicin and cisplatin in arginine deprived human cancer cells. PeerJ 4, e1542. doi:10.7717/peerj.1542
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi:10.1038/nprot.2016.095
Pino Del Carpio, D., Basnet, R. K., De Vos, R. C., Maliepaard, C., Visser, R., Bonnema, G., et al. (2011). The patterns of population differentiation in a Brassica rapa core collection. Theor. Appl. Genet. 122, 1105–1118. doi:10.1007/s00122-010-1516-1
Popolo, A., Adesso, S., Pinto, A., Autore, G., and Marzocco, S. (2014). L-Arginine and its metabolites in kidney and cardiovascular disease. Amino Acids 46, 2271–2286. doi:10.1007/s00726-014-1825-9
Reitzer, L. (2004). Biosynthesis of glutamate, aspartate, asparagine, L-alanine, and D-alanine. EcoSal Plus 1. doi:10.1128/ecosalplus.3.6.1.3
Rosenthal, G. A. (1977). The biological effects and mode of action of L-canavanine, a structural analogue of L-arginine. Q. Rev. Biol. 52, 155–178. doi:10.1086/409853
Rotmann, A., Simon, A., Martiné, U., Habermeier, A., and Closs, E. I. (2007). Activation of classical protein kinase C decreases transport via systems y+ and y+L. Am. J. Physiol. Cell Physiol. 292, C2259–C2268. doi:10.1152/ajpcell.00323.2006
Shimomura, A., Matsui, I., Hamano, T., Ishimoto, T., Katou, Y., Takehana, K., et al. (2014). Dietary L-lysine prevents arterial calcification in adenine-induced uremic rats. J. Am. Soc. Nephrol. 25, 1954–1965. doi:10.1681/asn.2013090967
Theron, L., Fernandez, X., Marty-Gasset, N., Chambon, C., Viala, D., Pichereaux, C., et al. (2013). Proteomic analysis of duck fatty liver during post-mortem storage related to the variability of fat loss during cooking of “foie gras”. J. Agric. Food Chem. 61, 920–930. doi:10.1021/jf302979q
Wang, X., Jiang, G., Kebreab, E., Li, J., Feng, X., Li, C., et al. (2020). (1)H NMR-based metabolomics study of breast meat from Pekin and Linwu duck of different ages and relation to meat quality. Food Res. Int. 133, 109126. doi:10.1016/j.foodres.2020.109126
Wu, G., Meininger, C. J., Knabe, D. A., Bazer, F. W., and Rhoads, J. M. (2000). Arginine nutrition in development, health and disease. Curr. Opin. Clin. Nutr. Metab. Care 3, 59–66. doi:10.1097/00075197-200001000-00010
Xiong, Z., Dankova, G., Howe, L. J., Lee, M. K., Hysi, P. G., De Jong, M. A., et al. (2019). Novel genetic loci affecting facial shape variation in humans. Elife 8, e49898. doi:10.7554/eLife.49898
Ye, M., Ye, F., He, L., Liu, Y., Zhao, X., Yin, H., et al. (2017). Molecular cloning, expression profiling, and marker validation of the chicken Myoz3 gene. Biomed. Res. Int. 2017, 5930918. doi:10.1155/2017/5930918
Zhan, H., Xiong, Y., Wang, Z., Dong, W., Zhou, Q., Xie, S., et al. (2022). Integrative analysis of transcriptomic and metabolomic profiles reveal the complex molecular regulatory network of meat quality in Enshi black pigs. Meat Sci. 183, 108642. doi:10.1016/j.meatsci.2021.108642
Zhang, L., Hao, Z., Zhao, C., Zhang, Y., Li, J., Sun, B., et al. (2021a). Taste compounds, affecting factors, and methods used to evaluate chicken soup: A review. Food Sci. Nutr. 9, 5833–5853. doi:10.1002/fsn3.2501
Keywords: liancheng white ducks, multi-omics, umami-tasting compounds, poultry meat, meat nutrition
Citation: Zhou H, Yang Y, Wang L, Ye S, Liu J, Gong P, Qian Y, Zeng H and Chen X (2022) Integrated multi-omic data reveal the potential molecular mechanisms of the nutrition and flavor in Liancheng white duck meat. Front. Genet. 13:939585. doi: 10.3389/fgene.2022.939585
Received: 11 May 2022; Accepted: 14 July 2022;
Published: 15 August 2022.
Edited by:
Elena E. Terenina, Institut National de la Recherche Agronomique de Toulouse, FranceReviewed by:
Cécile M. D. Bonnefont, GenPhySE, Université de Toulouse, FranceYali Dang, Ningbo University, China
Shaohua Shi, Fujian Academy of Agricultural Sciences, China
Copyright © 2022 Zhou, Yang, Wang, Ye, Liu, Gong, Qian, Zeng and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huijun Zeng, zenghuijuner2017@163.com; Xing Chen, Ligang1997@126.com