
==== Front
Poult Sci
Poult Sci
Poultry Science
0032-5791
1525-3171
Elsevier

S0032-5791(24)00790-9
10.1016/j.psj.2024.104211
104211
GENETICS AND MOLECULAR BIOLOGY
Genome-wide association studies reveal the genetic basis of growth and carcass traits in Sichuan Shelduck
Yang Zhao *
Xi Yang *
Qi Jingjing *
Li Liang *
Bai Lili *
Zhang Jun †
Lv Jia †
Li Bo ‡
Liu Hehe liuee1985@sicau.edu.cn
⁎1
⁎ Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 610000, China
† Rural Revitalization Development Service Centre, Zigong, China
‡ Farming Service Centre, Rong County, Zigong, China
1 Corresponding author: liuee1985@sicau.edu.cn
14 8 2024
11 2024
14 8 2024
103 11 10421125 4 2024
8 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
China has abundant local duck resource populations, and evaluating the characteristics of these breeds will help improve development and utilization. In this study, we conducted the first investigations of growth and slaughter performance on Sichuan Shelduck (n = 240), an endangered duck local breed. The average body weight is 1497.91 g at 90 d of age. According to the growth curve through data recorded every 2 wk, we observed a low relative growth rate (RGR) for the early growth stage. The RGR shows a decreasing trend with age increasing in the stage from 0 to 56 d of age. The SNP-based heritability estimation showed the growth rate has a relatively high heritability, indicating high genetic stability for this trait. In the correlation analysis, the percentage of leg muscle is positively correlated with the absolute growth rate (AGR) at 28 to 42 d of age, whereas it is negatively correlated with the earlier stages, exhibiting a time-specific correlation result. Additionally, genome-wide association studies (GWAS) identified PCSK6, TOX2, and TOMM7 as potential candidate genes influencing AGR (42–56) and AGR (56–90), while the candidate genes of slaughter traits were PTP4A2, FAM110B, TOX, UBXN2B, and FCHSD2. These results provide an important reference for further understanding the genetic basis of growth and meat production performance of Sichuan Shelduck.

Key words

GWAS
Sichuan Shelduck
growth
meat production performance
candidate gene
==== Body
pmcINTRODUCTION

Globally, particularly in Asia, ducks have become an important part of meat consumption and developed towards multi-type consumption ways (Tai and Tai, 2001; Magdelaine, et al., 2008; Whitton, et al., 2021). Consequently, breeders aim to develop diverse duck breeds to meet future market demands. The availability of abundant duck genetic resources is crucial for breeding improved varieties (Liu, et al., 2021). Investigation and screening of the genes underlying these excellent traits of local ducks can also provide a scientific basis for genetic improvement.

Recently, High-throughput sequencing and genome-wide association studies (GWAS) have been widely applied to farm animals to screen candidate genes and mutations related to duck growth and production traits (Ding et al., 2022; Li, et al., 2022). Notably, functional genes and loci have been discovered in various local duck breeds (Li, et al., 2023c). In terms of growth and development, whole-genome data elucidated the genetic pattern of Zhongshan ducks, highlighting the roles of SEMA5B and MIB1 in growth and development (Chen, et al., 2023). In Jinding ducks, Shanma ducks, and Shaoxing ducks, Whole genome resequencing was used to identify growth-related genes (IGF1R) and skeletal development genes (CDF5) (Li, et al., 2023b). The researchers identified EIF2AK3 gene was associated with body weight development in Shanma duck, Jinding duck, Gaoyou duck, Maple Leaf duck, Shaoxing duck (Zhu, et al., 2021). Studies on Gaoyou and Jinding ducks revealed that different expression patterns of the MYOD and MYF6 genes in breast muscle and leg muscle may be related to muscle development and differentiation (Li, et al., 2014). In addition, Pekin duck is renowned globally as a high-quality meat breed and has been extensively studied (Zhang, et al., 2018). The GWAS results indicated that IGF2BP1 gene was a key candidate gene affecting body size and feed efficiency in Peking ducks (Zhou, et al., 2018). And different genotypes of MAGI3 may influence breast muscle development of Peking ducks, and then affect the thickness of the perimysium (Tang, et al., 2023). Zhu et al. showed that solute carrier proteins (SLC39A10, SLC10A2) may affect the growth of meat-type Peking ducks (Zhu, et al., 2019). Deng et al. identified PLXDC2 associated with breast width and TNS3 associated with fossil bone length (Deng, et al., 2019). In terms of slaughtering performance, LOC101791418, TUBGCP3 and ATP11A play an important role in eviscerated weight and leg muscle weight percentage (Deng, et al., 2019). CTDSPL and PKP1 are new candidate genes for breast muscle thickness used by genotyping-by-sequencing (Deng, et al., 2020). Furthermore, genotyping and GWAS analyses have identified 2 SNPs (chr29:2,296,787 and chr29:2,296,832) in the AUTS2 gene that were related to high leanness in Pekin ducks (Liang, et al., 2023).

Despite the increasing demand for higher growth rates and meat production of ducks, the development and utilization of local duck resources have been neglected. The Sichuan Shelduck, well-known for its small body size, early maturity, strong adaptability, well-developed chest and leg muscles, and delicious meat, is a dual-purpose duck breed for meat and egg production. Combining its advantages with fast growth and high egg production from duck genetic resources can result in economically favorable hybridization combinations. However, there is currently a lack of comprehensive evaluation of growth and meat production traits, unclear genetic basis, and insufficient molecular markers, which limits its application in hybrid breeding. This study conducted a systematic evaluation of the growth performance and slaughter performance of both males and females. Additionally, we used GWAS to identify candidate genes and associated SNPs. The research results can enhance our understanding of the genetic basis of Sichuan Shelduck and promote its application in genetic improvement.

MATERIALS AND METHODS

Animals

All ducks utilized in this experiment were sourced from the Rong County Hemp Duck Breeding Farm in Zigong City, Sichuan Province. A total of 240 ducks were included for the assessment of growth performance, body size, and slaughter performance. These ducks were raised according to standard feeding procedures and nutritional requirements (Supplementary Table S1). The experimental ducks were initially kept in the brooding room for 0 to 14 d. From 15 to 90 d of age, the ducks were transferred to a ground flat breeding house for feeding. Throughout the experiment, the same feeding density and environment were maintained, and the ducks had access to food ad libitum. Body weight was measured every 2 wk from 0 to 90 d of age. At the age of 90 d, 240 ducks were randomly selected for slaughter, during which slaughter performance metrics, breast muscle samples, and fresh blood samples were collected for subsequent analyses. The Animal Care and Use Committee of the Institute of Animal Sciences, Sichuan Agricultural University, approved all animal experimentation and ducks. The guidelines approved the methods and protocols used in this study. The guidelines for the methods and protocols used in this study were also approved. All efforts were made to minimize any potential suffering experienced by the ducks.

Collection of Phenotypic Value

Size and weight. Body weight gain between 2 weighing times is defined as the absolute growth rate (AGR). The relative growth rate (RGR) is determined by the ratio of AGR to the previous weight. The calculation formulas are as follows (W0: measured weight at last time; W1: measured weight at next time):AGR(g)=W1−W0

RGR(%)=W1−W0W0×100%

at 90 d of age, body measurements were collected for all ducks. Half-Divina dept (H-Dd): Measured with a tape measure from the tip of the beak to the midpoint of the line connecting the hip bones. Neck length (NL): Distance from the first neck cone to the base of the neck. Breast width (BW): Measured with calipers as the distance between the 2 shoulder joints. Shank Length (SL): Measured with calipers as the straight-line distance from the upper joint of the tibia to the midpoint between the third and fourth toes. Fossil bone length (FBL): Measured with a tape measure as the distance between the 2 ends of the keel bone.

Slaughter trait. Slaughter Weight (SW): The weight of the duck after 6 h of fasting before slaughter. Dressed Weight (DW): The weight of the carcass after bloodletting and removal of feathers, foot cuticles, toe shells, and beak shells. Half-eviscerated weight (H-EW): The weight of the carcass after removal of the trachea, esophagus, crop, intestines, spleen, pancreas, gallbladder, reproductive organs, gizzard contents, and keratin membrane. Eviscerated Weight (EW): Half-eviscerated weight minus heart, liver, proventriculus, gizzard, and abdominal fat. Abdominal fat weight (AFW): Weight of abdominal fat and fat around the muscular stomach. Skin fat weight (SFW): Weight of skin and subcutaneous fat. The dressed percentage (DP) was calculated as the DW / SW × 100%. The percentage of half-eviscerated yield (H-EYP) was calculated as the H-EW/SW × 100%. The eviscerated yield percentage (EYP) was calculated as the EW/SW × 100%. The breast muscle percentage (BMP) was calculated as the breast muscle weight (BMW)/EW × 100%. The leg muscle percentage (LMP) was calculated as the leg muscle weight (LMW)/EW × 100%. The abdominal fat percentage (AFP) was calculated as the AFW/(EW+AFW) × 100%. The skin fat percentage (SFP) was calculated as (SFW + AFW) / EW × 100%.

DNA Extraction and Whole-Genome Re-Sequencing

The DNA samples were all extracted from duck whole blood using a standard Phenol-Chloroform extraction protocol. The quality of DNA was assessed using a NanoDrop2000 (Thermo. NanoDrop2000C) instrument and 0.8% agarose gel electrophoresis. After the examinations, standard procedures were used to generate paired-end libraries for each eligible sample. The average insert size was 500 bp, and the average read length was 150 bp. All libraries were sequenced on an Illumina HiSeq X-Ten platform or HiSeq 4000 platform to an average raw read sequence coverage of 5 × for the duck population. The raw reads were filtered using the NGS QC (v2.3.3) Toolkit with the default parameter. The depth ensured the accuracy of variant calling and genotyping and met the requirements for population genetic analyses.

Genomic SNP calling

The raw sequence data was aligned to the duck reference genome (ZJU 1.0, RefSeq assembly accession: GCF_015476345.1) using Burrows-Wheeler alignment (BWA aln) with default parameters (Li and Durbin, 2009). SNP calling was performed exclusively using GATK (version 3.5) (DePristo, et al., 2011), and the output was further filtered using VCFtools (version 0.1.15) (Danecek, et al., 2011). The SNPs were filtered based on the following criteria: (1) the SNPs were required to have a minor allele frequency > 0.03 (2) the maximum missing rate was < 0.1; and (3) the SNPs with only 2 alleles; (4) SNPs with heterozygosity greater than 80% were removed. After filtering, 17,967,246 SNPs were finally obtained and used for the subsequent analysis.

SNP-Based Heritability Calculation

Using genotypic and phenotypic data, a genetic relationship matrix (GRM) was constructed using Plink software. By combining GRM and phenotypic data, Genome-wide Complex Trait Analysis (GCTA) was used to partition the phenotypic variance into genetic and nongenetic parts to estimate heritability.

GWAS

GWAS was performed on the Sichuan Shelduck population to detect genomic regions that affect growth and meat production traits in ducks using the mixed linear model program EMMAX (Kang, et al., 2010). The random effect was the phylogenetic matrix estimated by all genome-wide SNPs. The whole-genome significance cutoff was defined as the Bonferroni threshold, 0.05 / Total SNPs (−log10 (P) = 8.56). The linear model is as follows:y=Xα+Zβ+Wμ+e

where y is the vector of phenotypic values of indicators, Xα is the fixed effects. To reduce the impact of gender on the results, “gender” was used as a covariate. Zβ represents the effect of SNP, and β represents the allele substitution effect. Wμ represents random animal effects with a variance–covariance structure based on the kinship matrix estimated using whole-genome SNP genotypes, and e is random residuals for data.

The Manhattan plots were generated using R (v3.5.1) and the ‘qqman’ package (Turner, 2018). The QQ plots were also generated to detect false positives from population stratification. The QQ plot shows that the ordinate represents the observed SNP p-value, while the abscissa represents the theoretical p-value generated using a chi-squared distribution. SNP variant information was combined to screen for potential molecular markers. Gene annotation was performed using SnpEff software, and individual genotypes for candidate regions were extracted using VCFtools in a Linux environment.

Candidate Gene Selection

The KEGG and GO enrichment was conducted using KOBAS 3.0 online tools (https://bioinfo.org/kobas/gene list/). Finally, we performed GO and KEGG analysis and visualization using R software packages “cluster Profiler”, “enrich plot”, and “ggplot2”.

Statistical Analysis

The study utilized statistical methods to calculate the mean and standard deviation for various phenotypic data. An inter-group difference analysis was conducted using independent sample T-tests with SPSS to assess differences in slaughter indicators among different genders in Sichuan Shelduck. Spearman was used to perform a correlation analysis between growth rate and meat production performance.

RESULTS

Descriptive Statistics

To better evaluate the growth rate of the research population, we collected the body weight data of ducks every 2 wk from hatching to 56 d of age (Supplementary Table S2). At 90 d of age, the body weight was 1556 ± 139.5 g for males and 1459 ± 136.29 g for females. The growth curve results indicate that ducks exhibit relatively rapid growth before 42 d, gradually slowing down after 56 d of age (Figure 1A, Supplementary Table S3). In addition, we calculated the RGR and AGR for all ducks at 4 stages before 56 d of age. The RGR and AGR patterns were slightly different during the peak growth period. The peak of AGR appeared at 28 to 42 d of age. Both male and female ducks showed decreasing growth for the first time at 42 to 56 d of age (Figures 1B and 1C). The RGRs of the population reached a maximum at the beginning and then decreased gradually (Figures 1D and 1E). The body size results of male and female ducks showed that all body size indicators showed significant differences between male and female ducks (P < 0.01), except SL (Figure 1F).Figure 1 The descriptive statistics of phenotypic traits of Sichuan Shelducks. (A) Male and female ducks' body weight growth curve from birth to 90 d old. The X-axis and Y-axis showed the ages and body weights, respectively. (B–E) Absolute growth rate of female and male ducks for each growth stage, respectively. Dots in the bar graphs represent the RGR or AGR of the individual sex group. x- and y-axes represent different growth stages and growth rate data. (F) Comparison of body size differences between males and females of Sichuan Shelduck. H-Dd, Half-Divina dept; NL, Neck length; BW, Breast width; SL, Shank length; FBL, Fossil bone length.

Figure 1

In this study, 16 slaughter performance indexes of Sichuan Shelducks were measured and calculated (Table 1). The results showed that the slaughter weight of male ducks was 1520.67 ± 152.95 g, and the female ones were 1429.74 g. The DP of male ducks was 78.029%, and the female ones were 79.005%. The DP and H-EYP reached more than 70%. The BMP and LMP were 10.84% ​​and 12.60%, respectively. The average LP of ducks was 23.43%. In addition, there were significant differences between males and females in some carcass indicators, including EW, LMW, AFW, BMP, SFP and LP (P < 0.05).Table 1 Comparison of carcass phenotypes between male and female ducks.

Table 1Trait type	Trait	Males	Female	P	
N	Mean ± SD	N	Mean ± SD	
Carcass yield	SW(g)	120	1540.92 ± 148.76	113	1451.42 ± 139.60	<0.01	
DW(g)	120	1202.83 ± 131.39	113	1147.52 ± 129.74	<0.01	
DP (%)	120	78.03 ± 0.03	113	79.01 ± 0.04	<0.05	
H-EW (g)	119	1104.71 ± 123.34	113	1053.01 ± 126.01	<0.01	
H-EYP (%)	119	71.61 ± 0.03	113	72.48 ± 0.04	0.05	
EW (g)	118	1004.24 ± 114.22	113	960.35 ± 111.73	<0.01	
EYP (%)	118	65.10 ± 0.03	113	66.12 ± 0.04	<0.05	
LMW (g)	120	126.135 ± 15.64	113	120.39 ± 15.70	<0.01	
LMP (%)	118	12.62 ± 0.01	113	12.59 ± 0.01	0.9	
BMW (g)	120	104.88 ± 22.02	113	109.52 ± 19.06	0.09	
BMP (%)	118	10.40 ± 0.01	113	11.38 ± 0.01	<0.01	
AFW (g)	86	6.28 ± 5.01	96	8.91 ± 5.62	<0.01	
AFP (%)	85	0.58 ± 0.00	96	0.88 ± 0.00	<0.01	
SFW (g)	119	187.32 ± 177.84	113	179.61 ± 43.98	0.65	
SFP (%)	85	20.65 ± 0.22	96	19.95 ± 0.04	0.76	
LP (%)	118	23.01 ± 0.02	113	23.97 ± 0.02	<0.01	
Internal organs (g)	GW	120	46.91 ± 7.64	113	39.12 ± 6.32	<0.01	
GSW	119	4.28 ± 0.78	113	4.08 ± 0.91	0.08	
HW	120	9.83 ± 1.49	113	9.41 ± 1.19	<0.05	
LW	120	34.33 ± 4.99	113	32.14 ± 5.01	<0.01	
SW	115	0.58 ± 0.24	106	0.58 ± 0.26	0.89	
Note: SW, Slaughter weight; DW, Dressed weight; DP, Dressed percentage; H-EW, Half-eviscerated weight; H-EYP, Percentage of half-eviscerated yield; EW, Eviscerated weight; EYP, Percentage of eviscerated yield; LMW, Leg muscle weight; LMP, Leg muscle percentage; BMW, Breast muscle weigh; BMP, Breast muscle percentage; AFW, Abdominal fat weight; AFP, Abdominal fat percentage; SFW, Skin fat weight; SFP, Skin fat percentage; LP, Percentage of lean meat BMT, Breast muscle thickness; GW, Gizzard weight; GSW, Glandular stomach weight; HW, Heart weight; LW, Liver weight; SW, Spleen weight. N is the number of samples; the mean is the average value; S.D. is the standard deviation. The P-values by T-test for the significance between groups of males and females were provided.

Phenotypic Correlation Analysis

The results of the phenotypic correlation analysis showed a significant correlation between meat production performance and body size development indicators (Figure 2, Supplementary Table S4). BW showed substantial correlations with DP (R = 0.218, P < 0.01) and H-EYP (R = 0.166, P < 0.05). In terms of growth and development, AGR (14–28), AGR (28–42), and AGR (56–90) were significantly correlated with the FBL in body size (P < 0.01). This suggests that FBL may be an important factor influencing growth rates during these periods. The AGR (28–42) was negatively correlated with most traits, except for a positive correlation with LMP. Additionally, LMP showed a significant negative correlation with AGR (14–28). It can be inferred that the growth in weight from 28 to 42 d of age primarily depends on the development of leg muscles. The correlation analyses further revealed significant correlations between different meat production phenotypes such as AFP and H-EYP (R = 0.344), EYP (R = 0.327) and LMP (R = −0.222) (P < 0.01).Figure 2 Correlation analysis chart among phenotypic indicators. H-EYP, percentage of half-eviscerated yield; EYP, Eviscerated yield percentage; LMP, Leg muscle percentage; BMP, Breast muscle percentage; AFP, Abdominal fat percentage; SFP, Skin fat percentage; LP, Percentage of lean meat; BMT, Breast muscle thickness; BMA, Breast muscle area. The depth of colors in the chart indicates the strength of the correlation, with positive and negative correlations represented by warm and cool colors, respectively. Asterisks serve as markers for significance, with more asterisks indicating higher significance. Deeper colors and more asterisks signify a stronger correlation. A single table is diagonally symmetrical, with diagonal values set to 1. ** P < 0.01 and * P < 0.05.

Figure 2

Genetic Parameter Estimation Based on SNPs

The heritability based on SNP is the extent to which it affects an individual trait or a specific genetic characteristic, reflecting the strength of association between the genotype and phenotype of the SNP. The SNP-based heritability estimation showed (Figure 3, Supplementary Table S5) that the heritability differs among stages. Among the growth traits, RGR (28–42) has the highest heritability, followed by 0-day-old weight, AGR (14–28), and 90-day-old, all exceeding a heritability of 0.5. Body size traits with higher heritability included H-Dd (0.37) and FBL (0.27). Regarding slaughter traits, the heritability for 8 indicators, SFW, H-EW, EW, DW, LMW, SW, GW, and BMP exceeded 0.5. This may indicate that genetic factors largely influence these traits.Figure 3 SNP-based genetic parameter estimation The Y-axis represents SNP heritability, and the X-axis represents phenotype. The higher the bar, the higher the heritability.

Figure 3

GWAS Screening for Candidate Genes Related to Production Performance

The GWAS results were derived from a dataset comprising 17,967,246 genomic SNP. We conducted GWAS analysis on 6 growth development and slaughter performance phenotypes: AGR (42–56), AGR (56–90), EW, BMW, AFW, and LP. The top 20 SNP in the GWAS results were selected for gene annotation.

Growth

Based on the GWAS analysis, the AGR (42-56) results showed clear signals on the chromosomes in the Manhattan map (Figure 4A), with 21 SNP sites surpassing the threshold line. Through SNP annotation, we discovered PCSK6, LOC101795342, TOX2, LOC101795089, JADE2, PCDH11X, LOC110352402, TOMM22, LOC106019054 and FSTL4 genes (Supplementary Table S6). The AGR (56–90) GWAS analysis results showed 22 significant SNP sites above the threshold line (Figure 4B). These loci were associated with TOMM7, SEMA5A, CDH22, LOC113842992 and LOC101795342 genes according to the SNP annotation (Supplementary Table S6).Figure 4 The Manhattan and quantile-quantile (QQ) plots of genome-wide association analysis of AGR (42-56), AGR (56-90), and GO functional enrichment analysis plots. A Manhattan and QQ plots of AGR (42-56). B, Manhattan and QQ plots of AGR (56-90). Each point represents one SNP, and abscissa numbers represent different chromosomes. The dotted line in the Manhattan plot represents the threshold level. C, GO enrichment of genes identified in AGR (42-56). D, GO enrichment of genes identified in AGR (56-90). Each row represents an enriched function, and the length of the bar represents the enrich ratio, which is calculated as “input gene number”/“background gene number”. The bar's color is the same as the color in the circular network above, representing different clusters.

Figure 4

In addition, we performed GO enrichment analysis (Supplementary Table S7) on other identified genes to explore potential candidate genes responsible for duck growth rate. The GO terms are divided into 3 categories: Molecular Function (MF), Cellular Component (CC), and Biological Process (BP). In the AGR (42–56) group, multiple potential growth-related pathways and clusters were identified, including negative regulation of collateral sprouting, the nerve growth factor binding, secretion by cell, nerve growth factor binding, the multicellular organism development, among others, in which the PCSK6, JADE2, and PCDH11X genes were enriched (Figure 4C). In the GO analysis results of AGR (56-90), TOMM7 was found to be enriched in protein import into the mitochondrial matrix, regulation of protein stability, and integral component of membrane. The SEMA5A gene was also enriched in 9 other pathways (Figure 4D).

Slaughter Performance

GWAS results showed potential signals for 4 traits (EW, BMW, AFW, LP); however, no significant results were found after the Bonferroni correction. Therefore, we selected the top 20 SNPs for each of the 4 traits for annotation to explore candidate genes affecting meat performance (Supplementary Table S8). Specifically, the GWAS analysis results of EW and BMW identified (Figures 5A, 5B) a total of 40 SNPs associated with these 2 phenotypes (Chr3: 36123219 C>G, chr21: 3347530 T>C, chr21: 2334214 G>T, etc.). These SNPs were annotated to 5 genes, including LOC110352184, ASIP, ZNF341, CBFA2T2, and LOC101802175. TMEM68, TOX, FAM110B, PTP4A2, and UBXN2B genes were annotated in AFW, while the IQSEC1 and KIF13B genes were annotated in the LP.Figure 5 The Manhattan and quantile-quantile (QQ) plots of genome-wide association analysis of EW, BMW, AFW, LP, and GO functional enrichment analysis plots. A, B, C, and D, Manhattan plot of EW, BMW, AFW, and LP respectively, in slaughter performance indicators. Each point represents one SNP, and abscissa numbers represent different chromosomes. The dotted line in the Manhattan plot represents the threshold level. E, Plot of GO functional enrichment analysis of candidate genes in EW and BMW. F, Plot of GO functional enrichment analysis of candidate genes in AFW. G, Plot of GO functional enrichment analysis of candidate genes in LP. Each row represents an enriched function, and the length of the bar represents the enrich ratio, which is calculated as “input gene number”/“background gene number”. The bar's color is the same as the color in the circular network above, representing different clusters.

Figure 5

We then performed GO enrichment analysis on these genes to explore potential genes that determine duck meat production performance (Figures 5E–G, Supplementary Table S9). The findings indicated that the gene CBFA2T2, identified through annotation of SNPs in EW, exhibited enrichment in 3 significant GO terms: intestinal epithelial cell differentiation, negative regulation of the Notch signaling pathway, and positive regulation of neuron projection development. These terms fall within the domain of biological functions and are likely associated with animal development and nerve growth (Figure 5E). In abdominal fat, we identified multiple pathways related to fat deposition, including white fat cell differentiation and positive regulation of brown fat cell differentiation, among which the PRDM16 gene was enriched. In addition, the FCHSD2 gene was enriched in the positive regulation of Arp2/3 complex-mediated actin nucleation and regulation of actin filament polymerization pathways, which are related to actin. In addition, the gene IQSEC1 annotated in the LP indicator was enriched in actin cytoskeleton organization.

DISCUSSION

Growth and carcass traits are crucial components for the profitability of the duck meat production industry to meet the increasing need for top-notch yet cost-effective poultry meat. Understanding the genetic foundations of these traits can enhance meat yield and quality, ultimately boosting financial gains and ensuring long-term sustainability in the industry. Sichuan Shelduck has an important genetic resource, but its phenotypic and genetic characteristics are still unclear. In this study, we evaluated its phenotypic and genetic parameters. In general, this breed is relatively small in body size and has a slow growth rate, which is possibly helpful for balancing body weight and bone development (Jain and Singhal, 2012; Duggan, et al., 2016; Zhang, et al., 2019). Regarding heritability, animal breeding is a complex process, and pedigree errors are common in commercial breeding practices, which can lower the accuracy of genetic parameter estimation (Zhu, et al., 2019). Utilizing a relationship matrix based on high-density genotyping effectively addresses these errors (Goddard, et al., 2011; Munoz, et al., 2014). Our results showed that RGR (28–42) had the highest heritability, followed by day-old weight, AGR (14–28) (0.66), SFW (0.67). It indicates that these traits are genetically stable and predictable, and the selection potential for this trait could be high, implying significant improvement may be obtained during crossing with other varieties. Currently, there are limited studies on the heritability of slaughter performance indexes in Sichuan Shelducks. The present study is the first evaluation work on this breed, which helps to deepen the understanding of the genetic basis of growth and meat production in this breed. In addition, the results of the correlation analysis revealed that different tissue sites exhibit distinct developmental timelines. Notably, indicators such as breast muscle, abdominal fat and skin fat mainly develop early, while leg muscles develop later. We speculate that there may be temporal specificity in regulating gene expression in different muscle tissues. This result suggests that appropriate time point sampling must be considered during future studies of gene expression regulation in the breast and leg muscles.

Animal growth is a complex physiological process influenced by various factors (Hu, et al., 2015). This study identified multiple SNPs associated with growth weight traits in 240 Sichuan Shelduck. Among all the identified SNPs, 41 were found in the GWAS results for AGR (42–56) and AGR (56–90), annotated to 6 genes, with PCSK6 being the most significant. PCSK6, identified to be highly expressed in muscle tissue, encodes a protease involved in muscle cell differentiation, proliferation, and apoptosis, which is crucial for regulating muscle growth (Rykaczewska, et al., 2020). PCSK6 gene knockout in mice results in abnormal bone development (Malfait, et al., 2012). Additionally, PCSK6 may affect energy metabolism regulation (Li, et al., 2023a). LH-regulated PCSK6 inhibits apoptosis of human granulosa cells through activin A and TGFb2 activation (Wang, et al., 2014). Overall, the impact of the PCSK6 gene on animal growth is multifaceted. The TOX2 gene primarily influences animal growth by promoting the differentiation of T follicular helper cells (Xu, et al., 2019), thus affecting the immune response through the proliferation and maturation of B lymphocytes. In the AGR (56–90) GWAS results, the TOMM7 gene was found. TOMM7's gene product, a subunit of the translocase of the outer mitochondrial membrane, regulates the assembly and stability of the translocase complex (Hönlinger, et al., 1996). TOMM7 was identified to be associated with severe growth retardation and premature aging traits (Garg, et al., 2022). These genes may serve as key candidate genes in the later growth stages of Sichuan Shelducks.

In terms of meat production performance, the presence of the CBFA2T2 gene was observed in the EW of Sichuan Shelduck. This gene is related to the weight of chickens at various growth stages (Gu, et al., 2011; Sharma, et al., 2015). The encoded product of CBFA2T2 is believed to regulate skeletal development. CBFA2T2 is essential for osteogenic differentiation of mesenchymal stem cells (Huang, et al., 2018). Moreover, CBFA2T2 facilitates adipogenic differentiation of mesenchymal stem cells by regulating CEBPA (Luo, et al., 2020). It influences cell function and tissue development, including being a key regulator of fat generation (Luo, et al., 2020), suggesting an association with muscle growth, fat metabolism, and potential impact on meat production performance. Abdominal fat, the primary deposition site for body fat, is highly correlated with overall body fat. In the GWAS results of this study, AFW displayed a distinct signal, with annotated genes including PTP4A2, FAM110B, TOX, UBXN2B, and FCHSD2. PTP4A2, identified as a protein tyrosine phosphatase, demonstrates higher expression in female elephant adipose tissue than males (Nilsson, et al., 2014). It has been suggested that PTP4A2/PRL2 may regulate leptin receptor signaling in adipose tissue (Huan, et al., 2003). FAM110B has been associated with the lean meat rate of pigs (Wang, et al., 2022). In addition, in animals (Machado, et al., 2022), UBXN2B, TOX, and FAM110B genes were related to carcass weight and quality (Alam, et al., 2023). For example, the FAM110B and TOX genes influence meat marbling (Li and Kim, 2015; Seabury, et al., 2017; Srikanth, et al., 2020). TOX and FAM110B impact reproductive traits in cattle(Fortes, et al., 2011), influencing maternal ability and calving ease in Nellore females (Silva, et al., 2020). UBXN2B's enrichment in various cellular functional pathways underscores its importance in cell development (Machado, et al., 2022). TOX, FAM110B, and UBXN2B have been identified to affect multiple traits (Naserkheil, et al., 2022), including growth, birth weight, carcass weight, average daily gain, feed intake, meat tenderness, height, and stature, in different beef cattle breeds (Magalhães, et al., 2016; Seabury, et al., 2017; Brunes, et al., 2021). FCHSD2, involved in cytoskeletal remodeling (Abouelezz and Almeida-Souza, 2022), may play a role in processes such as cell proliferation, differentiation, and migration, which impact animal fat production and metabolism. Studies have linked FCHSD2 to milk fat percentage (Marete, et al., 2018), although further research is needed for confirmation.

CONCLUSION

In summary, Sichuan Shelducks are generally small in body size, and slow-growing rapid, with a body weight of 1497.91 g at 90 d of age. However, the performance of the meat was better. The DP and H-EYP reached more than 70%. The BMP and LMP were 10.84% ​​and 12.60%, respectively. The average LP of ducks was 23.43%. GWAS highlights PCSK6, TOX2, and TOMM7 as candidate genes influencing later-stage growth (42–56, 56–90 d). PTP4A2, FAM110B, TOX, and UBXN2B also affect AFW, while CBFA2T2 is associated with EW. The identified SNPs serve as valuable molecular markers for selecting growth rate and meat production traits, enhancing production efficiency for Sichuan Shelducks.

DISCLOSURES

The authors declare no conflicts of interest.

Appendix Supplementary materials

Image, application 1

Image, application 2

Image, application 3

Image, application 4

Image, application 5

Image, application 6

Image, application 7

Image, application 8

Image, application 9

ACKNOWLEDGMENTS

This work was supported by grants from the National Key R&D Program of China (2022YFF1000100 , 2023YFD1300302), the China Agriculture Research System of Waterfowl (CARS-42), and Key Technology Support Program of Sichuan Province (2021YFYZ0014 ). Demonstration Extension of Breeding Ability Selection System of Nonghua Drakes (22ZDYFZF0005 ).

Author contributions: ZY is responsible for constructing article ideas and drafting manuscripts. HL and LL participated in the writing instruction and revision of the manuscript. YX and JQ are responsible for technical guidance. In addition, ZY, LB, JZ, JL, and BL conducted sample collection and data statistics. All listed authors have made substantial contributions to the research and publication. Thank you to all the authors for their contributions to the study.

Supplementary material associated with this article can be found in the online version at doi:10.1016/j.psj.2024.104211.
==== Refs
REFERENCES

Abouelezz A. Almeida-Souza L. The mammalian endocytic cytoskeleton Eur J Cell Biol 101 2022 151222
Alam M.Z. Haque M.A. Iqbal A. Lee Y.M. Ha J.J. Jin S. Park B. Kim N.Y. Won J.I. Kim J.J. Genome-wide association study to identify QTL for Carcass traits in Korean Hanwoo cattle Animals (Basel) 13 2023 2737 37685003
Brunes L.C. Baldi F. Lopes F.B. Lôbo R.B. Espigolan R. Costa M.F.O. Stafuzza N.B. Magnabosco C.U. Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle J. Anim. Breed. Genet. 138 2021 23 44 32654373
Chen H. Luo K. Wang C. Xuan R. Zheng S. Tang H. Li Y. Xiong Y. Wu Y. Wang L. Ouyang J. Yan X. Genomic characteristics and selection signals of Zhongshan ducks Animal 17 2023 100797
Danecek P. Auton A. Abecasis G. Albers C.A. Banks E. DePristo M.A. Handsaker R.E. Lunter G. Marth G.T. Sherry S.T. McVean G. Durbin R. The variant call format and VCFtools Bioinf. 27 2011 2156 2158
Deng M.T. Zhang F. Zhu F. Yang Y.Z. Yang F.X. Hao J.P. Hou Z.C. Genome-wide association study reveals novel loci associated with fat-deposition and meat-quality traits in Pekin ducks Anim. Genet. 51 2020 953 957 32844456
Deng M.T. Zhu F. Yang Y.Z. Yang F.X. Hao J.P. Chen S.R. Hou Z.C. Genome-wide association study reveals novel loci associated with body size and carcass yields in Pekin ducks BMC Genom. 20 2019 1
DePristo M.A. Banks E. Poplin R. Garimella K.V. Maguire J.R. Hartl C. Philippakis A.A. del Angel G. Rivas M.A. Hanna M. McKenna A. Fennell T.J. Kernytsky A.M. Sivachenko A.Y. Cibulskis K. Gabriel S.B. Altshuler D. Daly M.J. A framework for variation discovery and genotyping using next-generation DNA sequencing data Nat. Genet. 43 2011 491 498 21478889
Ding J. Ying F. Li Q. Zhang G. Zhang J. Liu R. Zheng M. Wen J. Zhao G. A significant quantitative trait locus on chromosome Z and its impact on egg production traits in seven maternal lines of meat-type chicken J. Anim. Sci. Biotechnol. 13 2022 96 35941697
Duggan B.M. Hocking P. Clements D. Gait in ducks (Anas platyrhynchos) and chickens (Gallus gallus)–similarities in adaptation to high growth rate Biol. Open 5 2016 1077 1085 27387535
Fortes M.R. Reverter A. Nagaraj S.H. Zhang Y. Jonsson N.N. Barris W. Lehnert S. Boe-Hansen G.B. Hawken R.J. A single nucleotide polymorphism-derived regulatory gene network underlying puberty in 2 tropical breeds of beef cattle J. Anim. Sci. 89 2011 1669 1683 21357453
Garg A. Keng W.T. Chen Z. Sathe A.A. Xing C. Kailasam P.D. Shao Y. Lesner N.P. Llamas C.B. Agarwal A.K. Mishra P. Autosomal recessive progeroid syndrome due to homozygosity for a TOMM7 variant J. Clin. Invest. 132 2022 e156864
Goddard M.E. Hayes B.J. Meuwissen T.H. Using the genomic relationship matrix to predict the accuracy of genomic selection J. Anim. Breed. Genet. 128 2011 409 421 22059574
Gu X. Feng C. Ma L. Song C. Wang Y. Da Y. Li H. Chen K. Ye S. Ge C. Hu X. Li N. Genome-wide association study of body weight in chicken F2 resource population PLoS One 6 2011 e21872 21779344
Hönlinger A. Bömer U. Alconada A. Eckerskorn C. Lottspeich F. Dietmeier K. Pfanner N. Tom7 modulates the dynamics of the mitochondrial outer membrane translocase and plays a pathway-related role in protein import EMBO J. 15 1996 2125 2137 8641278
Hu Y. Liu H. Song C. Xu W. Ji G. Zhu C. Shu J. Li H. Profiles of mRNA expression of related genes in the duck hypothalamus-pituitary growth axis during embryonic and early post-hatch development Gene. 559 2015 38 43 25577952
Huan J.N. Li J. Han Y. Chen K. Wu N. Zhao A.Z. Adipocyte-selective reduction of the leptin receptors induced by antisense RNA leads to increased adiposity, dyslipidemia, and insulin resistance J. Biol. Chem. 278 2003 45638 45650 12925533
Huang H. Dou L. Song J. Luo J. CBFA2T2 is required for BMP-2-induced osteogenic differentiation of mesenchymal stem cells Biochem. Biophys. Res. Commun. 496 2018 1095 1101 29378183
Jain V. Singhal A. Catch up growth in low birth weight infants: striking a healthy balance Rev. Endocr. Metab. Disord. 13 2012 141 147 22415299
Kang H.M. Sul J.H. Service S.K. Zaitlen N.A. Kong S.Y. Freimer N.B. Sabatti C. Eskin E. Variance component model to account for sample structure in genome-wide association studies Nat Genet 42 2010 348 354 20208533
Li C. Guo Z. Liu F. An P. Wang M. Yang D. Tang Q. PCSK6 attenuates cardiac dysfunction in doxorubicin-induced cardiotoxicity by regulating autophagy Free. Radic. Biol. Med. 203 2023 114 128 37061139
Li H. Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform Bioinf. 25 2009 1754 1760
Li H. Zhu C. Tao Z. Xu W. Song W. Hu Y. Zhu W. Song C. MyoD and Myf6 gene expression patterns in skeletal muscle during embryonic and posthatch development in the domestic duck (Anas platyrhynchos domestica) J. Anim. Breed. Genet. 131 2014 194 201 24180358
Li L. Quan J. Gao C. Liu H. Yu H. Chen H. Xia C. Zhao S. Whole-genome resequencing to unveil genetic characteristics and selection signatures of specific pathogen-free ducks Poult. Sci. 102 2023 102748
Li X. Xin A. Ma L. Gou X. Fang S. Dong X. Ni B. Tang L. Zhu L. Yan D. Kong X. Molecular genetic characterization and meat-use functional gene identification in Jianshui yellow-brown ducks through combined resequencing and transcriptome analysis Front. Vet. Sci. 10 2023 1269904
Li Y. Kim J.J. Multiple Linkage Disequilibrium Mapping Methods to Validate Additive Quantitative Trait Loci in Korean Native Cattle (Hanwoo) Asian. Australas. J. Anim. Sci. 28 2015 926 935 26104396
Li Y.D. Xue B. Xin L. Wang W.J. Li Z.W. Ning W. Fan X. Gao H.H. Guo H.S. Hui L. Wang S.Z. Integration of genome-wide association study and selection signatures reveals genetic determinants for skeletal muscle production traits in an F2 chicken population J. Integr. Agric. 21 2022 2065 2075
Liang S. Guo Z. Luo D. Tang J. Ji Z. Xie M. Hou S. Two variants of AUTS2 gene are associated with high lean meat percentage in Pekin ducks Gene. 848 2023 146864
Liu H. Wang L. Guo Z. Xu Q. Fan W. Xu Y. Hu J. Zhang Y. Tang J. Xie M. Zhou Z. Hou S. Genome-wide association and selective sweep analyses reveal genetic loci for FCR of egg production traits in ducks Genet. Sel. Evol. 53 2021 98 34930109
Luo J. Dou L. Yang Z. Zhou Z. Huang H. CBFA2T2 promotes adipogenic differentiation of mesenchymal stem cells by regulating CEBPA Biochem. Biophys. Res. Commun. 529 2020 133 139 32703401
Machado P.C. Brito L.F. Martins R. Pinto L.F.B. Silva M.R. Pedrosa V.B. Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems Animals (Basel) 12 2022 3526 36552446
Magalhães A.F. de Camargo G.M. Fernandes G.A.J. Gordo D.G. Tonussi R.L. Costa R.B. Espigolan R. Silva R.M. Bresolin T. de Andrade W.B. Takada L. Feitosa F.L. Baldi F. Carvalheiro R. Chardulo L.A. de Albuquerque L.G. Genome-wide association study of meat quality traits in Nellore cattle PLoS One 11 2016 e0157845
Magdelaine P. Spiess M. Valceschini E. Poultry meat consumption trends in Europe World's Poult. Sci. J. 64 2008 53 64
Malfait A.M. Seymour A.B. Gao F. Tortorella M.D. Le Graverand-Gastineau M.P. Wood L.S. Doherty M. Doherty S. Zhang W. Arden N.K. Vaughn F.L. Leaverton P.E. Spector T.D. Hart D.J. Maciewicz R.A. Muir K.R. Das R. Sorge R.E. Sotocinal S.G. Schorscher-Petcu A. Valdes A.M. Mogil J.S. A role for PACE4 in osteoarthritis pain: evidence from human genetic association and null mutant phenotype Ann. Rheum. Dis. 71 2012 1042 1048 22440827
Marete A.G. Guldbrandtsen B. Lund M.S. Fritz S. Sahana G. Boichard D. A meta-analysis including pre-selected sequence variants associated with seven traits in three French dairy cattle populations Front. Genet. 9 2018 522 30459810
Munoz P.R. Resende Jr M.F. Huber D.A. Quesada T. Resende M.D. Neale D.B. Wegrzyn J.L. Kirst M. Peter G.F. Genomic relationship matrix for correcting pedigree errors in breeding populations: impact on genetic parameters and genomic selection accuracy Crop. Sci. 54 2014 1115 1123
Naserkheil M. Manzari Z. Dang C.G. Lee S.S. Park M.N. Exploring and identifying candidate genes and genomic regions related to economically important traits in Hanwoo cattle Curr. Issues. Mol. Biol. 44 2022 6075 6092 36547075
Nilsson E.M. Fainberg H.P. Choong S.S. Giles T.C. Sells J. May S. Stansfield F.J. Allen W.R. Emes R.D. Mostyn A. Mongan N.P. Yon L. Molecular characterization of adipose tissue in the African elephant (Loxodonta Africana) PLoS One 9 2014 e91717 24633017
Rykaczewska U. Suur B.E. Röhl S. Razuvaev A. Lengquist M. Sabater-Lleal M. van der Laan S.W. Miller C.L. Wirka R.C. Kronqvist M. Gonzalez Diez M. Vesterlund M. Gillgren P. Odeberg J. Lindeman J.H. Veglia F. Humphries S.E. de Faire U. Baldassarre D. Tremoli E. Lehtiö J. Hansson G.K. Paulsson-Berne G. Pasterkamp G. Quertermous T. Hamsten A. Eriksson P. Hedin U. Matic L. PCSK6 is a key protease in the control of smooth muscle cell function in vascular remodeling Circ. Res. 126 2020 571 585 31893970
Seabury C.M. Oldeschulte D.L. Saatchi M. Beever J.E. Decker J.E. Halley Y.A. Bhattarai E.K. Molaei M. Freetly H.C. Hansen S.L. Yampara-Iquise H. Johnson K.A. Kerley M.S. Kim J. Loy D.D. Marques E. Neibergs H.L. Schnabel R.D. Shike D.W. Spangler M.L. Weaber R.L. Garrick D.J. Taylor J.F. Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle BMC Genom. 18 2017 386
Sharma A. Lee J.S. Dang C.G. Sudrajad P. Kim H.C. Yeon S.H. Kang H.S. Lee S.H. Stories and challenges of genome wide association studies in livestock - a review Asian. Australas. J. Anim. Sci. 28 2015 1371 1379 26194229
Silva D.B.S. Fonseca L.F.S. Pinheiro D.G. Magalhães A.F.B. Muniz M.M.M. Ferro J.A. Baldi F. Chardulo L.A.L. Schnabel R.D. Taylor J.F. Albuquerque L.G. Spliced genes in muscle from Nelore Cattle and their association with carcass and meat quality Sci. Rep. 10 2020 14701 32895448
Srikanth K. Lee S.H. Chung K.Y. Park J.E. Jang G.W. Park M.R. Kim N.Y. Kim T.H. Chai H.H. Park W.C. Lim D. A gene-set enrichment and protein-protein interaction network-based GWAS with regulatory SNPs identifies candidate genes and pathways associated with carcass traits in Hanwoo cattle Genes (Basel) 11 2020 316 32188084
Tai C. Tai J.J.L. Future prospects of duck production in Asia J. Poult. Sci. 38 2001 99 112
Tang H. Liu D. Zhang H. Fan W. Hu J. Xu Y. Guo Z. Huang W. Hou S. Zhou Z. Genome-wide association studies demonstrate the genes associated with perimysial thickness in ducks Anim. Genet. 54 2023 363 374 36697366
Turner S.D. qqman: An R package for visualizing GWAS results using Q-Q and Manhattan plots J. Open. Sour. Softw. 3 2018 731
Wang H. Wang X. Yan D. Sun H. Chen Q. Li M. Dong X. Pan Y. Lu S. Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology BMC Genom. 23 2022 594
Wang Y. Wang X.H. Fan D.X. Zhang Y. Li M.Q. Wu H.X. Jin L.P. PCSK6 regulated by LH inhibits the apoptosis of human granulosa cells via activin A and TGFβ2 J Endocrinol. 222 2014 151 160 24860148
Whitton C. Bogueva D. Marinova D. Phillips C.J.C. Are we approaching peak meat consumption? Analysis of meat consumption from 2000 to 2019 in 35 countries and its relationship to gross domestic product Animals (Basel) 11 2021 3466 34944243
Xu W. Zhao X. Wang X. Feng H. Gou M. Jin W. Wang X. Liu X. Dong C. The transcription factor tox2 drives t follicular helper cell development via regulating chromatin accessibility Immunit. 51 2019 826 839.e825
Zhang H.Y. Zeng Q.F. Bai S.P. Wang J.P. Ding X.M. Xuan Y. Su Z.W. Fraley G.S. Zhang K.Y. Study on the morphology and mineralization of the tibia in meat ducks from 1 to 56 d Poult. Sci. 98 2019 3355 3364 30916353
Zhang Y.R. Zhang L.S. Wang Z. Liu Y. Li F.H. Yuan J.M. Xia Z.F. Effects of stocking density on growth performance, meat quality and tibia development of Pekin ducks Anim. Sci. J. 89 2018 925 930 29682864
Zhou Z. Li M. Cheng H. Fan W. Yuan Z. Gao Q. Xu Y. Guo Z. Zhang Y. Hu J. Liu H. Liu D. Chen W. Zheng Z. Jiang Y. Wen Z. Liu Y. Chen H. Xie M. Zhang Q. Huang W. Wang W. Hou S. Jiang Y. An intercross population study reveals genes associated with body size and plumage color in ducks Nat. Commun. 9 2018 2648 30018292
Zhu F. Cheng S.R. Yang Y.Z. Hao J.P. Yang F.X. Hou Z.C. Genome-wide association study of growth and feeding traits in pekin ducks Front. Genet. 10 2019 702 31404312
Zhu T. Qi X. Chen Y. Wang L. Lv X. Yang W. Zhang J. Li K. Ning Z. Jiang Z. Qu L. Positive selection of skeleton-related genes during duck domestication revealed by whole genome sequencing BMC Ecol. Evol. 21 2021 165 34488647
