Genome-Wide Association Studies and Haplotype-Sharing Analysis Targeting the Egg Production Traits in Shaoxing Duck

Wenwu Xu; Zhenzhen Wang; Yuanqi Qu; Qingyi Li; Yong Tian; Li Chen; Jianhong Tang; Chengfeng Li; Guoqin Li; Junda Shen; Zhengrong Tao; Yongqing Cao; Tao Zeng; Lizhi Lu

doi:10.3389/fgene.2022.828884

Genome-Wide Association Studies and Haplotype-Sharing Analysis Targeting the Egg Production Traits in Shaoxing Duck

Front Genet. 2022 Mar 28:13:828884. doi: 10.3389/fgene.2022.828884. eCollection 2022.

Authors

Wenwu Xu¹, Zhenzhen Wang¹, Yuanqi Qu², Qingyi Li², Yong Tian¹, Li Chen¹, Jianhong Tang², Chengfeng Li², Guoqin Li¹, Junda Shen¹, Zhengrong Tao¹, Yongqing Cao¹, Tao Zeng¹, Lizhi Lu¹

Affiliations

¹ Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, China.
² Hubei Shendan Co., Ltd., Wuhan, China.

Abstract

Age at first egg (AFE) and egg number (EN) are economically important traits related to egg production, as they directly influence the benefits of the poultry industry, but the molecular genetic research that affects those traits in laying ducks is still sparse. Our objective was to identify the genomic regions and candidate genes associated with AFE, egg production at 43 weeks (EP43w), and egg production at 66 weeks (EP66w) in a Shaoxing duck population using genome-wide association studies (GWASs) and haplotype-sharing analysis. Single-nucleotide polymorphism (SNP)-based genetic parameter estimates showed that the heritability was 0.15, 0.20, and 0.22 for AFE, EP43w, and EP66w, respectively. Subsequently, three univariate GWASs for AFE, EP43w, and EP66w were carried out independently. Twenty-four SNPs located on chromosome 25 within a 0.01-Mb region that spans from 4.511 to 4.521 Mb were associated with AFE. There are two CIs that affect EP43w, i.e., twenty-five SNPs were in strong linkage disequilibrium region spanning from 3.186 to 3.247 Mb on chromosome 25, a region spanning from 4.442 to 4.446 Mb on chromosome 25, and two interesting genes, ACAD8 and THYN1, that may affect EP43w in laying ducks. There are also two CIs that affect EP66w, i.e., a 2.412-Mb region that spans from 127.497 to 129.910 Mb on chromosome 2 and a 0.355-Mb region that spans from 4.481 to 4.837 Mb on chromosome 29, and CA2 and GAMT may be the putative candidate genes. Our study also found some haplotypes significantly associated with these three traits based on haplotype-sharing analysis. Overall, this study was the first publication of GWAS on egg production in laying ducks, and our findings will be helpful to provide some candidate genes and haplotypes to improve egg production performance based on breeding in laying duck. Additionally, we learned from a method called bootstrap test to verify the reliability of a GWAS with small experimental samples that users can access at https://github.com/xuwenwu24/Bootstrap-test.

Keywords: egg production; gene; genome-wide association study; haplotype; laying duck.