Genome-wide association analysis of the primary feather growth traits of duck: identification of potential Loci for growth regulation

Identification of SNPs and INDELS associated with duck egg quality traits through a genome-wide association analysis

Egg quality traits are economically important in the poultry industry. To explore the genetic architecture and identify potential candidate genes, a genome-wide association study (GWAS) was performed for 13 egg quality traits using data from whole-genome sequencing of 299 Longyan Shan-ma female ducks, including 12 quantitative traits and one qualitative trait, eggshell color (ESC; white, light green, green). From estimation of pedigree genetic parameters, heritability (h2) ranged from 0.022 to 0.996 for the 12 quantitative traits, with the highest h2 (0.996) for eggshell color a* value (ESCA) and the lowest h2 (0.022) for egg yolk percentage relative to EW. A total of 8,874 single nucleotide polymorphism (SNP)-based significant associations (1.0 × 10-6) and 247 insertion-deletion (indel)-based significant associations (1.00 × 10-5) were identified, including 5,980 SNPs and 159 indel markers. From 5,924 SNPs and 143 indels associated with ESC traits, 181 potential candidate genes were identified, and most significant SNPs and indels (P < 1.0 × 10-20) were located at 1.86 Mb (44.29-46.15 Mb) on chromosome 4. The top SNP (chr4:45325309:C>A; P = 7.97 × 10-43) and the top indel (chr4:45299595:delTTCCACTCCAC; P = 4.20 × 10-36) for the ESC a* value were within two known ESC candidate genes; ATP-binding cassette subfamily G member 2 (ABCG2) and protein kinase cGMP-dependent 2 (PRKG2). Of 56 SNPs and 16 indels associated with other egg quality traits, 46 potential candidate genes were identified including synapse differentiation-inducing 1-like (SYNDIG1L) for EW, and core histone macro-H2A.1 (LOC101795967) and neurogenin 1 (NEUROG1) for egg shape index; and four genes including collagen type VI alpha 3 chain (COL6A3), lysine demethylase 7A (KDM7A), LOC101802169, and sperm-associated antigen 16 (SPAG16) for egg yolk weight and the percentage of yolk to total egg weight. Of the 46 genes, the molecular functions of 22 are related to protein binding, indicating important roles in the formation of egg quality traits. Our findings provide new insight into the genetic basis of egg quality traits in ducks.

Identification of candidate genes associated with primary feathers of tianfu nonghua ducks based on Genome-wide association studies

The primary feathers of ducks have important economic value in the poultry industry. This study quantified the primary feather phenotype of Nonghua ducks, including the primary feathers' length, area, distribution of black spots, and feather symmetry. And genome-wide association analysis was used to screen candidate genes that affect the primary feather traits. The genome-wide association study (GWAS) results identified the genetic region related to feather length (FL) on chromosome 2. Through Linkage disequilibrium (LD) analysis, candidate regions (chr2: 115,246,393-116,501,448 bp) were identified and were further annotated to 5 genes: MRS2, GPLD1, ALDH5A1, KIAA0319, and ATP9B. Secondly, candidate regions related to feather black spots were identified on chromosome 21. Through LD analysis, the candidate regions (chr21: 163,552-2,183,853 bp) were screened and further annotated to 47 genes. Among them, STK4, CCN5, and YWHAB genes were related to melanin-related pathways or pigment deposition, which may be key genes affecting the distribution of black spots on feathers. In addition, we also screened 125 genes on multiple chromosomes that may be related to feather symmetry. Among them, significant SNPs on chromosome 1 were further identified as candidate regions (chr1: 142,118,209-142,223,605 bp) through LD analysis and annotated into 2 genes, TGFBRAP1 and LOC113839965. These results reported the genetic basis of the primary feather from multiple phenotypes, and offered valuable insights into the genetic basis for the growth and development of duck feathers and feather color pattern.

Application of GWAS and mGWAS in Livestock and Poultry Breeding

In recent years, genome-wide association studies (GWAS) and metabolome genome-wide association studies (mGWAS) have emerged as crucial methods for investigating complex traits in animals and plants. These have played pivotal roles in research on livestock and poultry breeding, facilitating a deeper understanding of genetic diversity, the relationship between genes, and genetic bases in livestock and poultry. This article provides a review of the applications of GWAS and mGWAS in animal genetic breeding, aiming to offer reference and inspiration for relevant researchers, promote innovation in animal genetic improvement and breeding methods, and contribute to the sustainable development of animal husbandry.

Genome-wide association analysis explores the genetic loci of amino acid content in duck's breast muscle

BACKGROUND

Amino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.

RESULT

In the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55-56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2-21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.

CONCLUSION

In this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.

Identification of candidate genes affecting the tibia quality in Nonghua duck

The skeleton is a vital organ providing structural support in poultry. Weakness in bone structure can lead to deformities, osteoporosis, cage fatigue, and fractures, resulting in economic losses. Research has substantiated that genetic factors play a significant role in influencing bone quality. The discovery of genetic markers associated with bone quality holds paramount importance for enhancing genetic traits related to the skeletal system in poultry. This study analyzed nine phenotypic indicators of tibia quality in 120-day-old ducks. The phenotypic correlation revealed a high correlation among diameter, Perimeter, and weight (0.69-0.78), and a strong correlation was observed between toughness and breaking strength (0.62). Then, we conducted a genome-wide association analysis of the phenotypic indicators to elucidate the genetic basis of tibial quality in Nonghua ducks. Among the 11 candidate genes that were annotated, TAPT1, BST1, and STIM2 were related to the diameter indicator, ZNF652, IGF2BP1, CASK, and GREB1L were associated with the weight and toughness indicators. RFX8, GLP1R, and DNAAF5 were identified for ash, calcium, and phosphorus content, respectively. Finally, KEGG and GO analysis for annotated genes were performed. STIM2 and BST1 were enriched into the Calcium signalling pathway and Niacin and nicotinamide metabolic pathway, which may be key candidate genes affecting bone quality phenotypes. Gene expression analysis of the candidate genes, such as STIM2, BST1, TAPT1, and CASK showed higher expression levels in bones compared to other tissues. The obtained results can contribute to new insights into tibial quality and provide new genetic biomarkers that can be employed in duck breeding.

Genome-Wide Association Studies Reveal Candidate Genes Associated with Pigmentation Patterns of Single Feathers of Tianfu Nonghua Ducks

In modern advanced genetics and breeding programs, the study of genes related to pigmentation in ducks is gaining much attention and popularity. Genes and DNA mutation cause variations in the plumage color traits of ducks. Therefore, discovering related genes responsible for different color traits and pigment patterns on each side of the single feathers in Chinese ducks is important for genetic studies. In this study, we collected feather images from 340 ducks and transported them into Image Pro Plus (IPP) 6.0 software to quantify the melanin content in the feathers. Thereafter, a genome-wide association study was conducted to reveal the genes responsible for variations in the feather color trait. The results from this study revealed that the pigmented region was larger in the male ducks as compared to the female ducks. In addition, the pigmented region was larger on the right side of the feather vane than on the left side in both dorsal and ventral feathers, and a positive correlation was observed among the feather color traits. Further, among the annotated genes, WNT3A, DOCK1, RAB1A, and ALDH1A3 were identified to play important roles in the variation in pigmented regions of the various feathers. This study also revealed that five candidate genes, including DPP8, HACD3, INTS14, SLC24A1, and DENND4A, were associated with the color pigment on the dorsal feathers of the ducks. Genes such as PRKG1, SETD6, RALYL, and ZNF704 reportedly play important roles in ventral feather color traits. This study revealed that genes such as WNT3A, DOCK1, RAB1A, and ALDH1A3 were associated with different pigmentation patterns, thereby providing new insights into the genetic mechanisms of single-feather pigmentation patterns in ducks.

Whole-genome resequencing to unveil genetic characteristics and selection signatures of specific pathogen-free ducks

Specific pathogen-free ducks are important high-grade laboratory animals, with a key role in research related to poultry biosecurity, production, and breeding. However, the genetic characteristics of experimental duck varieties remain poorly explored. Herein we performed whole-genome resequencing to construct a single nucleotide polymorphism genetic map of the genomes of 3 experimental duck varieties [Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)] to determine their genetic characteristics and identify selection signatures. Subsequent analyses of population structure and genetic diversity revealed that each duck variety formed a monophyletic group, with SM showing richer genetic diversity than JD and SX. Further, on exploring shared selection signatures, we found 2 overlapping genomic regions on chromosome Z of all experimental ducks, which comprised immune response-related genes (IL7R and IL6ST). Moreover, growth and skeletal development (IGF1R and GDF5), meat quality (FoxO1), and stress resistance (HSP90B1 and Gpx8-b) candidate gene loci were identified in strongly selected signatures specific to JD, SM, and SX, respectively. Our results identified the population genetic basis of experimental ducks at the whole-genome level, providing a framework for future molecular investigations of genetic variations and phenotypic changes. We believe that such studies will eventually contribute to the management of experimental animal resources.