Skip to main content

Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis

Abstract

Background

The genetic locus responsible for duck body size has been fully explained before, but the growth trait-related genetic basis is still waiting to be explored. For example, the genetic site related to growth rate, an important economic trait affecting marketing weight and feeding cost, is still unclear. Here, we performed genome wide association study (GWAS) to identify growth rate-associated genes and mutations.

Result

In the current study, the body weight data of 358 ducks were recorded every 10 days from hatching to 120 days of age. According to the growth curve, we evaluated the relative and absolute growth rates (RGR and AGR) of 5 stages during the early rapid growth period. GWAS results for RGRs identified 31 significant SNPs on autosomes, and these SNPs were annotated by 24 protein-coding genes. Fourteen autosomal SNPs were significantly associated with AGRs. In addition, 4 shared significant SNPs were identified as having an association with both AGR and RGR, which were Chr2: 11483045 C>T, Chr2: 13750217 G>A, Chr2: 42508231 G>A and Chr2: 43644612 C>T. Among them, Chr2: 11483045 C>T, Chr2: 42508231 G>A, and Chr2: 43644612 C>T were annotated by ASAP1, LYN and CABYR, respectively. ASAP1 and LYN have already been proven to play roles in the growth and development of other species. In addition, we genotyped every duck using the most significant SNP (Chr2: 42508231 G>A) and compared the growth rate difference among each genotype population. The results showed that the growth rates of individuals carrying the Chr2: 42508231 A allele were significantly lower than those without this allele. Moreover, the results of the Mendelian randomization (MR) analysis supported the idea that the growth rate and birth weight had a causal effect on the adult body weight, with the growth rate having a greater effect size.

Conclusion

In this study, 41 SNPs significantly related to growth rate were identified. In addition, we considered that the ASAP1 and LYN genes are essential candidate genes affecting the duck growth rate. The growth rate also showed the potential to be used as a reliable predictor of adult weight, providing a theoretical reference for preselection.

Peer Review reports

Background

Meat is an important component of a healthy and well-balanced diet, providing high-quality protein and fat sources for the body [1]. However, with the rapid growth of the global population, determining ways to ensure an adequate food supply has become the focus of scientists [2]. Improving the growth performance of animals is perhaps one of the most direct methods. Growth performance can be affected by various factors, such as feed nutrition [3], feeding mode [4], and environmental factors [5]. In addition to those mentioned above, genetic factors are also one of the factors that cannot be ignored. The growth performance in different breeds of the same species often presents obvious differences [6, 7]. Therefore, improving animal growth performance at the genetic level has always been a major goal for breeders.

Originally, the candidate gene approach was a reliable way of identifying molecular markers and conducting genetic improvement [8,9,10]. However, with the present understanding of genetic research, the candidate gene approach has been shown to have some limitations in explaining complex traits. Fortunately, the maturity of high-throughput sequencing technology and the wide application of the genome-wide association study (GWAS) make it possible to comprehensively identify genetic loci affecting complex traits such as growth rate. To date, in farm animals such as pigs, cattle, and chickens, over 1,000 quantitative trait loci (QTLs) from the animal QTL database have been identified to be associated with growth-related traits such as body weight gain and average daily gain. A growth rate-related GWAS in a pig population including four breeds identified nine growth-related candidate genes [11]. A recent GWAS on growth traits in the Duroc population found that the most significant genetic region could explain 18.95% of the variation in average daily gain [12]. The above studies advanced our understanding of pig growth-related traits. In the Brahman cattle population, 167 and 262 significant SNPs were identified to be significantly associated with mature weight and maturity rate, respectively, and several candidate genes related to muscle and bone development were found [13]. In poultry, research on genetic mutations and candidate genes related to growth traits has also been ongoing. Recently, a study identified the TGF3 gene as a genetic marker for chicken growth traits by candidate gene approach [14]. However, a survey of the chicken pangenome has made great progress. Wang et al. successfully identified a deletion mutation in the promoter region of the IGF2BP1 gene as a causal mutation affecting growth traits such as chicken body size through GWAS and fine mapping methods [15]. Molecular marker-assisted selection using this mutation will greatly help to produce high-productivity chicken breeds. Interestingly, a few years ago, the same gene (IGF2BP1) in ducks was also identified as the gene responsible for affecting duck body size [16]. Subsequent studies also conducted a GWAS on Peking ducks and identified QTLs related to growth and feeding efficiency [2].

The duck is an important economic animal, especially in Asian countries. Although IGF2BP1 is responsible for body size in ducks, there are still many unknown growth-related genetic loci to be identified, such as the genetic site related to growth rate, which is one of the essential indices to evaluate growth performance and an important economic trait affecting the marketing age [17]. Therefore, in this study, we studied 358 ducks (Anas domesticus) and collected body weight data every 10 days from egg hatching (Day 1) to 120 days of age. A GWAS was conducted to identify genetic loci related to growth performance. The results deepen our understanding of the genetic basis of duck growth.

Results

Growth performance evaluation

To better evaluate the growth rate of the research population, we collected the body weight data of each duck every 10 days from hatching to 120 days of age and then constructed the growth curves of the male and female populations (Fig. 1A). As shown in the results, both male and female ducks were in a rapid growth stage before 40 days, but the growth rate gradually slowed down after 40 days, and the weight was stable after 50 days. Therefore, we calculated the RGRs and AGRs of 5 stages of all ducks before the age of 50 days, which were S1 (from birth weight to 10-day weight), S2 (from 10-day weight to 20-day weight), S3 (from 20-day weight to 30-day weight), S4 (from 30-day weight to 40-day weight) and S5 (from 40-day weight to 50-day weight). The RGR and AGR patterns showed a slight difference in the growth peak. The RGRs of the population reached a maximum at the beginning and then decreased gradually, while the peak of AGR appeared at S4. Both fell to the lowest level and showed a negative growth for the first time at S5 (Fig. 1B and C). Then, the relationship between growth rate and final adult body weight was investigated. With the body weight at 120 days as a reference, we compared the RGR and AGR differences between individuals with a heavier weight than average and a lower than average weight. At this stage, the weight of male and female ducks had already peaked and did not increase, indicating that the body weights of the ducks had attained somatic maturity. The results showed that the RGR and AGR were higher among heavier ducks at the age of 120 days (Fig. 1D and E). In particular, the AGR showed extremely significant differences (P < 0.01) at almost every growth stage (Fig. 1E). The correlation analysis also reflected similar results. AGRs and RGRs at each stage all showed varying degrees of positive correlation with the 50-day weight, while the RGRs could only explain the variation in 120-day body weight by 0.06% - 3.32%. For AGRs, this range was increased to 4.67% - 23.34% (Fig. 1F). The above results indicated that the AGR had a higher correlation with the final stable weight after the rapid growth stage of ducks and thus has the potential to be an important predictor for the early prediction of market weight.

Fig. 1
figure 1

Growth rate evaluation. A Male and female ducks' body weight growth curve from birth to 120-day old. The X-axis and Y-axis showed the ages and body weights, respectively. and C The RGRs B AGRs C of S1-S5 growth stages. The dots in the bar charts showed the RGR or AGR of one individual. The X-axis and Y-axis represent different growth stages and growth rates data. and E Comparison of RGRs and AGRs between the ducks heavier and lower than average body weight in 120-day old. The significance was shown as * (P < 0.05) and ** (P < 0.01). The contents represented by the horizontal and vertical coordinates were consistent with Fig. B and C. F The correlation analysis between growth rates and body weight data at 120-day old. The fitting curves and the Pearson correlation coefficient (R2) could reflect the correlation degree. The abscissa of each scatter plot represents the growth rates, and the ordinate represents 120-day body weight

GWAS based on RGR

A total of 4,237,004 filtered SNPs were obtained and used for the subsequent GWAS with a minor allele frequency (MAF) > 0.03 and a maximum missing rate < 0.1 among all 358 individuals from 5 synthetic lines (SL1-SL5) (Figure S1). Then, principal component analysis (PCA) was conducted to analyse the population structure. No obvious population stratification was found among the different synthetic lines, while 3 subpopulations were identified in SL1 ducks (Figure S2). Therefore, we used the first 3 principal components as the covariates to correct the influence of genetic backgrounds in the GWAS.

First, a GWAS based on the RGR was conducted. There were 31 SNPs identified on autosomes with a significant association with the RGRs of different growth stages by the criteria of –log10P > 8.59 (0.05/total SNP number) (Fig.  2A, Table S1). We carried out annotation analysis on these SNPs to investigate the genes related to these mutations (Table S1). After removing the pseudogenes and unannotated ncRNAs, a total of 24 protein-coding genes were identified, including a gene widely related to growth, growth-associated protein-43 (GAP43), which was related to the SNP significantly associated with the RGR at the S1 stage (Table S1), implying that this gene might participate in the early growth of ducklings. KEGG enrichment analysis (Table S2) was conducted for the other identified genes to explore the potential candidate genes that determine the duck growth rate. Only 6 known KEGG pathways in which the ERN1, GFPT2, and ASAP1 genes participate were identified (Table S2). In addition, the GO analysis provided more information. We clustered the KEGG and GO terms with similar functions (Table S2), and the top 5 clusters contained most pathways and GO functions (Fig. 2B). Multiple potential growth-related pathways and clusters were identified, including the regulation of growth and nervous system development in Cluster 1, in which the GAP43 gene was enriched (Fig. 2B, Table S2). In addition, we also found genes regulating microtubule development (Cluster 2, CAMSAP1), bone mineralization (Cluster 3, ANKH), cell proliferation and differentiation (Cluster 7: LYN), and adipose tissue development (Cluster 15, EBF2) (Table S2). Above all, these results indicated that the genes related to the significant SNP might play essential roles in duck growth and weight determination.

Fig. 2
figure 2

SNPs associated with the RGRs. A GWAS between genotypes and RGRs (left) and QQ plots (right). The red lines represent the Bonferroni corrected significance threshold (-log10P = 8.59). And the red dots were the significant SNPs that crossed the threshold. GWAS results of RGRs from S1 – S5 were shown from inside to outside. B The KEGG pathways (diamonds) and GO functions (circles) in which significant related SNPs-related genes are enriched. Each node represents an enriched term, and the node colour represents different clusters. The higher the enriched factors, the smaller the node sizes are. Only the first 5 terms with highest enrich ratios in each cluster were marked with text

GWAS based on AGR

Then, we continued to analyse the genetic loci associated with AGR. In the GWAS results of S1-S5, only three stages of AGR (S1, S2, and S5) were found to have significantly associated SNPs (Fig. 3A). Fourteen autosomal SNPs were found to be significantly associated with AGRs at the above mentioned stages (Table S3). We carried out gene annotation for these SNPs, and the related genes were enriched in various signalling pathways and gene functions by KEGG and GO analysis (Fig. 3B, Table S3, Table S4). It is worth noting that the pathways and functions enriched in Cluster 4 were highly coincident with Cluster 7 of the previous RGR-related gene enrichment result, which is mainly associated with cell proliferation and differentiation (Table S2 and Table S4). Therefore, we investigated the intersection of AGR- and RGR-associated SNPs. A total of 4 shared SNPs were found to be significantly associated with the AGR and RGR in the S5 growth stage and were all located on chromosome 2 (Fig. 3C). The linkage disequilibrium (LD) analysis identified two potential linkage intervals, Block1 (Chr2: 11483045 C>T & Chr2: 13750217 G>A) and Block2 (Chr2: 42508231 G>A & Chr2: 43644612 C>T), among these 4 SNPs (Figure S3). The annotated genes of these mutations were ASAP1 (Chr2: 11483045 C>T), LYN (Chr2: 42508231 G>A), and CABYR (Chr2: 43644612 C>T) (Table S1 and Table S3). Among these, Chr2: 42508231 G>A had the highest –log10P value in the GWAS based on both AGR and RGR (Table S1 and Table S3) and was located in the intron region of the LYN gene, which was reported to participate in epithelial and haematopoietic cell functions, including cell differentiation, growth, proliferation, and migration (Table S2 and Table S4). We then genotyped the Chr2: 42508231 G>A mutation and investigated the growth rates of individuals with different genotypes. The results showed that AGR and RGR at the S5 stage were significantly lower in the ducks carrying the Chr2: 42508231 A allele (Fig. 3D). According to the TF binding site prediction of the mutation, two transcription regulator binding sites (YY1, MF3) were found to be possibly lost after the mutation of Chr2: 42508231 (Fig. 3E).

Fig. 3
figure 3

SNPs associated with the AGRs. A GWAS between genotypes and AGRs (left) and QQ plots (right). The red lines represent the Bonferroni corrected significance threshold (-log10P = 8.59). And the red dots were the significant SNPs that crossed the threshold. GWAS results of AGRs from S1 – S5 were shown from inside to outside. B The KEGG pathways (diamonds) and GO functions (circles) in which significant related SNPs-related genes are enriched. Each node represents an enriched term, and the node colour represents different clusters. The higher the enriched factors, the smaller the node sizes are. Same as Fig. 2, only the first 5 terms with highest enrich ratios in each cluster were marked with text. C The shared significant SNPs between the AGRs and RGRs GWAS results. D Phenotypic comparisons between different genotypes based on the most significant SNP, Chr2: 42508231 G>A. The X-axis represents the different genotypes, and the Y-axis shows the AGR and RGR in the S5 growth stage. E Changes in the binding sites of Chr2: 42508231 G>A were demonstrated. The sequence used in the analysis was 5 bp of the left and right flanks of the mutation site

Growth rate is one of the causal factors determining adult body weight

In the first part of the results, we found correlations between AGRs and body weight at 120 days of age. Therefore, we intended to explore the causal relationship between growth rate and adult weight and whether individuals with a high growth rate will grow heavier. We conducted a GWAS on the body weight phenotypic data at each stage. From 70 days of age, a stable signal appeared on chromosome 28. By gene annotations, this signal corresponded to the gene IGF2BP1, which was proven to be associated with duck body size in a previous study [16]. However, it was interesting that this signal was not found in any GWAS results of body weight or growth rates at the earlier development stages until the age of 70 days (Figure S4 - S15, Fig. 4A). It seemed that there was no obvious relationship between the adult weight and the growth rate of ducks from the simple GWAS results. To examine this further, we took all SNP sites related to growth rate as instrumental variables after removing the SNPs with mutual LD relationships and conducted Mendelian randomization analysis. First, we analysed the effect of RGRs on 120-day body weight. A total of 29 filtered SNPs were involved in the subsequent analysis. Little indication for the existence of directional horizontal pleiotropy for the RGRs was observed by visual inspection of the funnel plot, which was consistent with the results of the directional horizontal pleiotropy test (P > 0.05) (Fig. 4B, Table S5). The TwoSampleMR results showed that as RGR increased, the 120-day body weight increased as well by using MR-egger (Beta = 0.21, P = 3.04E-02; OR = 1.23, 95% CI: [1.03, 1.47]), weighted median (Beta = 0.24, P = 3.29E-03; OR = 1.28, 95% CI: [1.08, 1.50]), IVW (Beta = 0.35, P = 1.71E-10; OR = 1.41, 95% CI: [1.27, 1.57]), simple mode (Beta = 0.43, P = 8.36E-03; OR = 1.54, 95% CI: [1.14, 2.07]), and weighted mode (Beta = 0.27, P = 1.46E-03; OR = 1.31, 95% CI: [1.13, 1.52]). In addition, the above five MR methods all showed that the significant promoting effect of a single SNP on the 120-day body weight increased with the single SNP effect on the RGRs (Fig. 4B, Table S6). The same analysis process was performed on AGRs as well. A potential horizontal pleiotropy was found. However, the IVW method (Beta = 0.45, P = 1.80E-03; OR = 1.57, 95% CI: [1.18, 2.08]) still showed a significant positive causal effect of AGRs on the 120-day body weight (Fig. 4C, Table S5, Table S6). The above results indicated that the growth rate was indeed one of the essential causes affecting adult weight. Then, we further explored the causal effect of birth weight on the 120-day body weight. The results also supported the idea that birth weight had a significant causal impact on 120-day body weight by all 5 analysis methods (Fig. 4D, Table S6). However, the effect value (Beta: about 0.01) of birth weight on adult weight was much lower than that of the growth rate, implying that the growth rate was more able to predict the adult body than the birth weight.

Fig. 4
figure 4

Causal effect of growth rate and birth weight on the adult body weight. A GWAS based on the 120-day body weight. The gray line represents the Bonferroni corrected significance threshold (-log10P = 8.59). The x-axis shows the physical positions of each marker along the chromosome, and the y-axis shows the - log10P values of SNPs. Band D Funnel plot displaying instrument strength (y-axis) plotted against causal effect estimates (x-axis) for SNPs associated with RGR B, AGR C, and birth weight D. Scatter plots displaying estimates of the association between each SNP and 120-day body weight (y-axis) against estimates of the association between each SNP and RGR B, AGR C, and birth weight D (x-axis). The error bars on each of the points indicated the 95% CI. The causal effect estimates of Inverse variance weighted, MR Egger, Simple mode, Weighted median, and Weighted mode were represented by solid lines in different colours

Discussion

Body growth is a physiological process affected by multiple complex factors. In the present study, we identified multiple SNPs associated with duck growth weight in 358 ducks. Among all the SNPs identified, 4 were shared in both the AGR and RGR GWAS results, and these SNPs were annotated by 3 genes: ASAP1 (Chr2: 11483045 C>T), CABYR (Chr2: 43644612 C>T) and LYN (Chr2: 42508231 G>A). ASAP1 is a GTPase-activating protein (GAP) that can negatively regulate ADP ribosylation factor (Arf). The Arf family is mainly involved in regulating membrane vesicle transport and the cytoskeleton [18]. Moreover, Arf was proven to be associated with early embryonic development in mice and regulates the growth, migration, and invasion of cancer cells [19, 20]. Studies on the ASAP1 gene have shown that this gene also plays an essential role in regulating growth and development. In mice, loss of the ASAP1 gene can lead to growth retardation, delayed ossification, adipocyte development, and fat depot formation reduction [21]. The SNPs related to the ASAP1 gene have also been associated with production and growth traits in cattle [22, 23]. The CABYR gene has been widely proven to be related to sperm morphology and sperm development [24, 25], but no research has yet shown that this gene is involved in the growth performance of animals or plants. However, the other gene located in the same LD block as CABYR, LYN, which belongs to the Src family of intracellular membrane-associated tyrosine kinases, is an important candidate gene for growth traits in many farm animals. Like most genes regulating growth traits, LYN plays an essential role in the occurrence and development of cancers [26,27,28] since it can function in the signal transduction of growth factor receptors [29]. LYN is a vital candidate gene affecting human height [30]. In animals, the LYN gene was identified to be associated with various growth traits of cattle, such as feed intake, growth, body weight, and birth weight. [31,32,33,34]. In chickens, the LYN gene was identified as being associated with comb development [35]. One study explained that the possible route by which the LYN gene affects growth traits is by regulating the intestinal microbial composition, 2 SNPs were identified in LYN and PLAG1 that were significantly associated with Methanobacterium, which had a higher abundance in high-weight female chickens [36]. In addition, the YY1 gene, which was predicted to lose a binding site at the regulatory region of the LYN gene, was previously proven to participate in fat metabolism. The loss of the YY1 protein in mice can protect against diet-induced obesity and increase energy consumption [37]. From the prediction result and similarity of biological function, a hypothesis of the potential regulatory relationship between the YY1 transcription factor and the LYN gene could be debated. However, the actual situation still needs further verification.

A previous study determined that IGF2BP1 is the major gene affecting the body shape of ducks [16], while in the present study, IGF2BP1 signals were only identified in the GWAS results of body weight after 60 days of age. This indicates that the IGF2BP1-related mutation may not be able to evaluate early growth traits very well. However, early growth trait indicators, such as birth weight, are usually used to predict adult body weight and conduct early selection. Meanwhile, it is unclear whether ducks' birth weight affects their subsequent growth performance. Studies in other species are still controversial as well. Some researchers have suggested that the birth weight of pigs has a greater impact on pig growth performance after weaning [38], while another study suggested that birth weight did not always determine subsequent growth performance [39]. The MR analysis showed that both birth weight and growth rate had causal effects on the 120-day slaughter weight, while the impact of birth weight on slaughter weight was slight. From the aspect of effect size, it could be even said that the duck birth weight is not so important to postnatal growth performance which is similar to the research result in other terrestrial vertebrates [40]. In contrast, the early growth rate can be used as a more reliable predictor of final adult body weight. In addition, these results indirectly indicated that SNPs that affect growth rate also play a role in duck adult body weight as minor effect sites to some extent.

Based on the above description, we considered the ASAP1 and LYN genes to be important candidate genes affecting duck growth rates. Moreover, the SNPs we identified could be used as potential molecular markers to select the growth rate phenotype, leading to improvements in the marketing age and growth performance.

Materials and methods

Animal samples and phenotype collection

A total of 358 ducks (Anas domesticus) (Male: n = 168; Female: n = 180; Record missing: n = 10) from 5 synthetic lines (SL1-SL5) bred by poultry breeding farm of Sichuan Agricultural University were included in the main research materials in the present study. The duck genetic background of each population is relatively independent, and there is no cross between them. The DNA samples were extracted from the blood cells of all 358 individuals by phenol-chloroform extraction method for genomic re-sequencing. The experimental ducks were reared in the brooding room for 0-14 days. The feeding density was 30 ducks / m2, and they were free-fed all day. During this time, the lighting time was 24 hours a day. At the age of 15~120 days, the ducks were fed in the ground flat breeding house, and the feeding density and lighting time was 1.5 ducks / m2 and 8~10 hours/day. At the age of 15-60 days, all ducks had 12 hours of free feeding. After 60 days of age, the feeding restriction was conducted, and the average daily intake was 0.13 Kg/duck. During the feeding period, each duck shall be weighed every ten days. And the gain of body weight between two weighing times is defined as AGR. The ratio of AGR to the previous weight is considered as the RGR. The calculation formulas are as follows (W0: measured weight at last time; W1: measured weight at next time):

$$\mathrm{AGR}={\mathrm{W}}_{1}- {\mathrm{W}}_{0}$$
$$\mathrm{RGR}=\frac{{\mathrm{W}}_{1}- {\mathrm{W}}_{0}}{{\mathrm{W}}_{0}}\times 100\%$$

All ducks were euthanized at 120 days old by inhaling carbon dioxide and cervical dislocation, which performed by competent personnel who experienced and correctly applied the technique. The Animal Welfare Committee approved the protocols for all Sichuan Agricultural University animal experiments, and all methods strictly obeyed the Guide for the Care and Use of Agricultural Animals in Research and Teaching . The metadata of all samples were collected in the Table S7.

Genomic SNP identification

The DNA samples were all extracted from duck whole blood using a standard Phenol-Chloroform extraction protocol. The quality and quantity of DNA were examined using a NanoDrop2000 device and by agarose gel electrophoresis. After the examinations, standard procedures generated 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 or HiSeq 4000 platform to an average raw read sequence coverage of ×5 for the duck population. The depth ensured the accuracy of variant calling and genotyping and met the requirements for population genetic analyses. Raw sequence data were mapped to the duck reference genome (ZJU 1.0, RefSeq assembly accession: GCF_015476345.1) with Burrows-Wheeler alignment (BWA aln) using the default parameters [41]. SNP calling was performed exclusively using GATK (version 3.5) [42], and the output was further filtered using VCFtools (version 0.1.15) [43]. The SNPs were filtered based on the following criteria: (1) the SNPs had to have a minor allele frequency > 0.03; (2) the maximum missing rate was < 0.1; and (3) the SNPs had only two alleles; (4) SNPs with heterozygosity greater than 80% will be removed. For all 358 individuals, over 4 million SNPs were finally obtained and used for the subsequent analysis.

PCA analysis

PCA was performed based on all SNPs using GCTA tool [44]. This software applies PCA to genetic data to analyse the population structure. The figures were plotted using the first, second and third principal components with R package ggplot2 (version: 3.3.6).

GWAS

Population structure and cryptic relatedness were used to minimize false positives and to increase statistical power. The mixed linear model program Emmax was used for the association analysis [45]. The first three principal component values (PCA eigenvectors) derived from the whole-genome SNPs and the forward/ backward crosses were fitted as fixed effects in the mixed model to correct population stratification [46]. The random effect was the kinship matrix estimated from all individual whole-genome SNPs. We defined the whole-genome significance cutoff as the Bonferroni threshold, which is 0.05/total SNPs (−log10P =8.59).

Two-sample mendelian randomization analysis

Two-sample MR analyses were performed using the GWAS results data. After removing the LD relationship, all SNPs associated with AGRs or RGRs and the top 20 SNPs associated with birth weight were chosen as the instruments for the MR analyses. To obtain estimates of the causal effect of adult weight, we performed Inverse variance weighted, MR Egger, Simple mode, Weighted median, and Weighted mode analysis to evaluate the causal effect of growth rate and birth weight on the adult body weight. Analyses were performed using the TwoSampleMR R package of MR-Base [47].

Data statistics and other bioinformatics tools

All correlation analysis was performed using the Pearson correlation method. The unpaired t-test was used to analyse whether there were significant differences (P < 0.05). The KEGG and GO enrichment was conducted using KOBAS 3.0 online tools (http://bioinfo.org/kobas/genelist/). The trans-factor binding site prediction was performed by PROMO (http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3).

Availability of data and materials

All Genomic re-sequencing raw data files were submitted to the SRA database (Accession: PRJNA907492 & PRJNA907501).

Abbreviations

GWAS:

Genome-wide association studies

QTLs:

Quantitative trait loci

LD:

Linkage disequilibrium

SNP:

Single nucleotide polymorphism

MAF:

Minor allele frequency

PCA:

Principal component analysis

MR:

Mendelian Randomization

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GO:

Gene ontology

AGR:

Absolute growth rate

RGR:

Relative growth rate

CABYR:

Calcium Binding Tyrosine Phosphorylation Regulated

ASAP1:

Development and Differentiation-Enhancing Factor 1

LYN:

LYN Proto-Oncogene, Src Family Tyrosine Kinase

IGF2BP1:

Insulin-Like Growth Factor 2 MRNA Binding Protein 1

References

  1. Pereira PMdCC, Vicente AFdRBJMs. Meat nutritional composition and nutritive role in the human diet. Meat sci. 2013;93(3):586–92.

    Article  CAS  PubMed  Google Scholar 

  2. Zhu F, Cheng S-R, Yang Y-z, Hao J-P, Yang F-X, Hou Z-CJFiG. Genome-wide association study of growth and feeding traits in pekin ducks. Front Genet. 2019;10:702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Gadde U, Kim W, Oh S, Lillehoj HSJAhrr. Alternatives to antibiotics for maximizing growth performance and feed efficiency in poultry: a review. Anim Health Res Rev. 2017;18(1):26–45.

    Article  CAS  PubMed  Google Scholar 

  4. Li Y, Luo C, Wang J, Guo FJJoAAR. Effects of different raising systems on growth performance, carcass, and meat quality of medium-growing chickens. J Appl Anim Res. 2017;45(1):326–30.

    Article  CAS  Google Scholar 

  5. Goo D, Kim JH, Park GH, Delos Reyes JB, Kil DYJA. Effect of heat stress and stocking density on growth performance, breast meat quality, and intestinal barrier function in broiler chickens. Animals. 2019;9(3):107.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Hassan FA, Roushdy EM, Zaglool AW, Ali MA, El-Araby IEJSVR. Growth performance, carcass traits and economic values of Pekin, Muscovy, and Mulard ducks. Slov Vet Res. 2018;55(20):357–65.

    Google Scholar 

  7. Nasr MA, Abd-Elhamid T, Hussein MAJMs. Growth performance, carcass characteristics, meat quality and muscle amino-acid profile of different rabbits breeds and their crosses. Meat sci. 2017;134:150–7.

    Article  CAS  PubMed  Google Scholar 

  8. Anh NTL, Kunhareang S, Duangjinda MJA-Ajoas. Association of chicken growth hormones and insulin-like growth factor gene polymorphisms with growth performance and carcass traits in Thai broilers. Asian-Australas J Anim Sci. 2015;28(12):1686.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jiang R, Li J, Qu L, Li H, Yang NJAG. A new single nucleotide polymorphism in the chicken pituitary-specific transcription factor (POU1F1) gene associated with growth rate. Anim Genet. 2004;35(4):344–6.

    Article  CAS  PubMed  Google Scholar 

  10. Vykoukalova Z, Knoll A, Dvořák J, Čepica SJJoAB. Genetics: New SNPs in the IGF2 gene and association between this gene and backfat thickness and lean meat content in Large White pigs. J Anim Breed Genet. 2006;123(3):204–7.

    Article  CAS  PubMed  Google Scholar 

  11. Tang Z, Xu J, Yin L, Yin D, Zhu M, Yu M, Li X, Zhao S, Liu XJFig. Genome-wide association study reveals candidate genes for growth relevant traits in pigs. Front genet. 2019;10:302.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ruan D, Zhuang Z, Ding R, Qiu Y, Zhou S, Wu J, Xu C, Hong L, Huang S, Zheng EJG. Weighted single-step GWAS identified candidate genes associated with growth traits in a Duroc pig population. Genes. 2021;12(1):117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Crispim AC, Kelly MJ, Guimarães SEF, e Silva FF, Fortes MRS, Wenceslau RR, Moore SJPO. Multi-trait GWAS and new candidate genes annotation for growth curve parameters in Brahman cattle. PLoS One. 2015;10(10):e0139906.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Hosnedlova B, Vernerova K, Kizek R, Bozzi R, Kadlec J, Curn V, Kouba F, Fernandez C, Machander V, Horna HJA. Associations between IGF1, IGFBP2 and tgfß3 genes polymorphisms and growth performance of broiler chicken lines. Animals. 2020;10(5):800.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Wang K, Hu H, Tian Y, Li J, Scheben A, Zhang C, Li Y, Wu J, Yang L, Fan XJMB, et al. The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size. Mol Biol Evol. 2021;38(11):5066–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhou Z, Li M, Cheng H, Fan W, Yuan Z, Gao Q, Xu Y, Guo Z, Zhang Y, Hu JJNC. An intercross population study reveals genes associated with body size and plumage color in ducks. Nat Commun. 2018;9(1):1–10.

    Google Scholar 

  17. Prakash A, Saxena VK, Singh MKJWsPSJ. Genetic analysis of residual feed intake, feed conversion ratio and related growth parameters in broiler chicken: a review. Worlds Poult Sci J. 2020;76(2):304–17.

    Article  Google Scholar 

  18. Inoue H, Randazzo PAJT. Arf GAPs and their interacting proteins. Traffic. 2007;8(11):1465–75.

    Article  CAS  PubMed  Google Scholar 

  19. Hayakawa N, Ogoh H, Sumiyoshi M, Matsui Y, Nishikawa S, Miyamoto K, Maede Y, Kiyonari H, Suzuki M, Watanabe TJB, et al. The ADP-ribosylation factor 1 gene is indispensable for mouse embryonic development after implantation. Biochem Biophys Res Commun. 2014;453(4):748–53.

    Article  CAS  PubMed  Google Scholar 

  20. Boulay P-L, Cotton M, Melançon P, Claing AJJoBC. ADP-ribosylation factor 1 controls the activation of the phosphatidylinositol 3-kinase pathway to regulate epidermal growth factor-dependent growth and migration of breast cancer cells. J Biol Chem. 2008;283(52):36425–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Schreiber C, Saraswati S, Harkins S, Gruber A, Cremers N, Thiele W, Rothley M, Plaumann D, Korn C. Armant OJPg: Loss of ASAP1 in mice impairs adipogenic and osteogenic differentiation of mesenchymal progenitor cells through dysregulation of FAK/Src and AKT signaling. PLoS genetics. 2019;15(6): e1008216.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Kour A, Deb SM, Nayee N, Niranjan SK, Raina VS, Mukherjee A, Gupta ID. Patil CSJAB: Novel insights into genome-wide associations in Bos indicus reveal genetic linkages between fertility and growth. Anim Biotechnol. 2022;34(1):39–55.

    Article  Google Scholar 

  23. Tizioto P, Meirelles S, Veneroni G, Tullio R, Rosa A, Alencar M, Medeiros S, Siqueira F, Feijó G, Silva LJMs. A SNP in ASAP1 gene is associated with meat quality and production traits in Nelore breed. Meat Sci. 2012;92(4):855–7.

    Article  CAS  PubMed  Google Scholar 

  24. Zhao Y, Yang L, Su G, Wei Z, Liu X, Song L, Hai C, Wu D, Hao Z, Wu YJL. Growth Traits and Sperm Proteomics Analyses of Myostatin Gene-Edited Chinese Yellow Cattle. Life. 2022;12(5):627.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Hu R, Jiang X, Liu G, Chi SJT. Polygenic co-expression changes the testis growth, hormone secretion and spermatogenesis to prompt puberty in Hu sheep. Theriogenology. 2022;194:116–25.

    Article  CAS  PubMed  Google Scholar 

  26. Ingley EJCC. Functions of the Lyn tyrosine kinase in health and disease. Cell Communication and signaling. 2012;10(1):1–11.

    Article  Google Scholar 

  27. Goldenberg-Furmanov M, Stein I, Pikarsky E, Rubin H, Kasem S, Wygoda M, Weinstein I, Reuveni H, Ben-Sasson SAJCr. Lyn is a target gene for prostate cancer: sequence-based inhibition induces regression of human tumor xenografts. Cancer res. 2004;64(3):1058–66.

    Article  CAS  PubMed  Google Scholar 

  28. Nisar A, Kayani MA, Nasir W, Mehmood A, Ahmed MW, Parvez A, Mahjabeen IJMG. Fyn and Lyn gene polymorphisms impact the risk of thyroid cancer. Mol Genet Genomics. 2022;297(6):1649–59.

    Article  CAS  PubMed  Google Scholar 

  29. Hibbs ML, Harder KWJGF. The duplicitous nature of the Lyn tyrosine kinase in growth factor signaling. Growth Factors. 2006;24(2):137–49.

    Article  CAS  PubMed  Google Scholar 

  30. Tripaldi R, Stuppia L, Alberti SJBeBA-RoC. Human height genes and cancer. Biochim Biophys Acta Reviews on Cancer. 2013;1836(1):27–41.

    Article  CAS  Google Scholar 

  31. Roberts AJLS. Genome-wide association study for carcass traits in a composite beef cattle breed. Livestock Science. 2018;213:35–43.

    Article  Google Scholar 

  32. Terakado A, Costa R, De Camargo G, Irano N, Bresolin T, Takada L, Carvalho C, Oliveira H, Carvalheiro R, Baldi FJA. Genome-wide association study for growth traits in Nelore cattle. Animal. 2018;12(7):1358–62.

    Article  CAS  PubMed  Google Scholar 

  33. Ghoreishifar SM, Eriksson S, Johansson AM, Khansefid M, Moghaddaszadeh-Ahrabi S, Parna N, Davoudi P, Javanmard AJGSE. Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds. Genet Sel Evol. 2020;52(1):1–15.

    Article  Google Scholar 

  34. Seabury CM, Oldeschulte DL, Saatchi M, Beever JE, Decker JE, Halley YA, Bhattarai EK, Molaei M, Freetly HC, Hansen SLJBg. Genome-wide association study for feed efficiency and growth traits in US beef cattle. BMC genomics. 2017;18(1):1–25.

    Article  Google Scholar 

  35. Tu Y, Liu Y, Zhang M, Shan Y, Ji G, Ju X, Zou J, Shu JJJoPS. Identifying Signatures of Selection Related to Comb Development. J Poult Sci. 2021;1:58.

    Google Scholar 

  36. Ji J, Luo C, Zou X, Lv X, Xu Y, Shu D, Qu HJPs. Association of host genetics with intestinal microbial relevant to body weight in a chicken F2 resource population. Poult sci. 2019;98(9):4084–93.

    Article  CAS  PubMed  Google Scholar 

  37. Verdeguer F, Soustek MS, Hatting M, Blättler SM, McDonald D, Barrow JJ, Puigserver PJM. biology c: Brown adipose YY1 deficiency activates expression of secreted proteins linked to energy expenditure and prevents diet-induced obesity. Mol cell biol. 2016;36(1):184–96.

    Article  CAS  Google Scholar 

  38. Wolter B, Ellis M, Corrigan B, DeDecker JJJoas. The effect of birth weight and feeding of supplemental milk replacer to piglets during lactation on preweaning and postweaning growth performance and carcass characteristics. J anim sci. 2002;80(2):301–8.

    Article  CAS  PubMed  Google Scholar 

  39. Camp Montoro J, Manzanilla EG, Solà-Oriol D, Muns R, Gasa J, Clear O, Calderón Díaz JAJA. Predicting productive performance in grow-finisher pigs using birth and weaning body weight. Animals. 2020;10(6):1017.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Case TJ. On the evolution and adaptive significance of postnatal growth rates in the terrestrial vertebrates. Q rev biol. 1978;53(3):243–82.

    Article  CAS  PubMed  Google Scholar 

  41. Li H, Durbin RJb. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, DelAngel G, Rivas MA, Hanna MJNg. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43(5):491.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry STJB. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Yang J, Lee SH, Goddard ME, Visscher PMJTAJoHG. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88(1):76–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Kang HM, Sul JH, Service SK, Zaitlen NA, Kong Sy, Freimer NB, Sabatti C, Eskin EJNg. Variance component model to account for sample structure in genome-wide association studies. Nat genet. 2010;42(4):348–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich DJNg. Principal components analysis corrects for stratification in genome-wide association studies. Nat genet. 2006;38(8):904–9.

    Article  CAS  PubMed  Google Scholar 

  47. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J, Langdon RJe. The MR-Base platform supports systematic causal inference across the human phenome. eLife. 2018;7:e34408.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Du Sert NP, Ahluwalia A, Alam S, Avey MT, Baker M, Browne WJ, Clark A, Cuthill IC, Dirnagl U, Emerson MJPb. Reporting animal research: Explanation and elaboration for the ARRIVE guidelines 2.0. PLoS biol. 2020;18(7):e3000411.

    Article  Google Scholar 

Download references

Acknowledgments

Thanks to all the authors for their contributions to the study.

Funding

The cost of experimental design, sequencing services, results analysis, main experimental reagents and manuscript writing is funded by the National Key R&D Program of China (SQ2022YFF1000051), Sichuan Science and Technology Program (2021JDJQ0008) and the China Agriculture Research System of Waterfowl (CARS-42).

Author information

Authors and Affiliations

Authors

Contributions

All experiments were designed by YX and QW. YX, JQ, JL, HH, and HX performed the data analysis of the genome. JH, XY, and LB participated in the sample collection and preparation work. YX, YZ, CH, SH, and HL completed the manuscript writing. JW and LL participated in the writing instruction and revision of the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Hehe Liu or Liang Li.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with relevant guidelines and regulations. The protocols for all animal experiments were approved by the Animal Welfare Committee of Sichuan Agricultural University, and all methods strictly obeyed the Guide for the ARRIVE (Animal Research: Reporting of In Vivo Experiments) guidelines 2.0 [48].

Consent for publication

Not applicable.

Conflict of Interest

The authors declared no potential conflicts of interest with respect to the research, author-ship and publication of this article.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table S1.

The significant SNPs associated with RGRs and the gene annotation. Supplementary Table S2. The KEGG and GO enrichment of RGR-associated genes. Supplementary Table S3. The significant SNPs associated with AGRs and the gene annotation. Supplementary Table S4. The KEGG and GO enrichment of AGR-associated genes. Supplementary Table S5. Heterogeneity and directional horizontal pleiotropy test of SNPs associated with AGRs, RGRs, and birth weight. Supplementary Table S6. The MR analysis results between exposures of AGR, RGR, birth weight and the outcome of the 120-day body weight. Supplementary Table S7. The metadata of all samples used in the research.

Additional file 2: Supplementary Figure S1.

SNP distribution. The number of SNPs within a 1 Mb window size distribute across all chromosomes. The color from green to red represents the gradual increasing SNP number in each window. Supplementary Figure S2. Principal Component Analysis. PCA was conducted on the SNP information of each sample using the GCTA tool. The blue, red, green, cyan and black dots represent SL1 - SL5 populations, respectively. Supplementary Figure S3. Linkage disequilibrium analysis. The numbers in each diamond-shaped box represent the correlation coefficients R2 between two SNPs, and the R2 > 0.2 was considered to indicate linkage. Supplementary Figure S4-S15. GWAS based on body weight from birth (Figure S4) to 110 days old ( Figure S15).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xi, Y., Wu, Q., Zeng, Y. et al. Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis. BMC Genomics 24, 285 (2023). https://doi.org/10.1186/s12864-023-09302-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09302-8

Keywords

  • Duck
  • GWAS
  • Growth rate