Three chromosome-level duck genome assemblies provide insights into genomic variation during domestication.

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed "missing genes" in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Genome-wide association study reveals novel loci associated with feeding behavior in Pekin ducks.

BACKGROUND: Feeding behavior traits are an essential part of livestock production. However, the genetic base of feeding behavior traits remains unclear in Pekin ducks. This study aimed to determine novel loci related to feeding behavior in Pekin ducks. RESULTS: In this study, the feeding information of 540 Pekin ducks was recorded, and individual genotype was evaluated using genotyping-by-sequencing methods. Genome-wide association analysis (GWAS) was conducted for feeding behavior traits. Overall, thirty significant (P-value < 4.74E-06) SNPs for feeding behavior traits were discovered, and four of them reached the genome-wide significance level (P-value < 2.37E-07). One genome-wide significance locus associated with daily meal times was located in a 122.25 Mb region on chromosome 2, which was within the intron of gene ubiquitin-conjugating enzyme E2 E2 (UBE2E2), and could explain 2.64% of the phenotypic variation. This locus was also significantly associated with meal feed intake, and explained 2.72% of this phenotypic variation. CONCLUSIONS: This study is the first GWAS for feeding behavior traits in ducks. Our results provide a list of candidate genes associated with feeding behavior, and also help to better understand the genetic mechanisms of feeding behavior patterns in ducks.

Comparison of carcass and meat quality traits between lean and fat Pekin ducks.

OBJECTIVE: According to market demand, meat duck breeding mainly includes 2 breeding directions: lean Pekin duck (LPD) and fat Pekin duck (FPD). The aim of the present study was to compare carcass and meat quality traits between 2 strains, and to provide basic data for guidelines of processing and meat quality improvement. METHODS: A total of 62 female Pekin ducks (32 LPDs and 30 FPDs) were slaughtered at the age of 42 days. The live body weight and carcass traits were measured and calculated. Physical properties of breast muscle were determined by texture analyzer and muscle fibers were measured by paraffin sections. The content of inosine monophosphate (IMP), intramuscular fat (IMF) and fatty acids composition were measured by high-performance liquid chromatography, Soxhlet extraction method and automated gas chromatography respectively. RESULTS: The results showed that the bodyweight of LPDs was higher than that of FPDs. FPDs were significantly higher than LPDs in subcutaneous fat thickness, subcutaneous fat weight, subcutaneous fat percentage, abdominal fat percentage and abdominal fat shear force (p<0.01). LPDs were significantly higher than FPDs in breast muscle thickness, breast muscle weight, breast muscle rate and breast muscle shear force (p<0.01). The muscle fiber average area and fiber diameter of LPDs were significantly higher than those of FPDs (p<0.01). The muscle fiber density of LPDs was significantly lower than that of FPDs (p<0.01). The IMF of LPDs in the breast muscle was significantly higher than that in the FPDs (p<0.01). There was no significant difference between the 2 strains in IMP content (p>0.05). The polyunsaturated fatty acid content of LPDs was significantly higher than that of FPDs (p<0.01), and FPDs had higher saturated fatty acid and monounsaturated fatty acid levels (p<0.05). CONCLUSION: Long-term breeding work resulted in vast differences between the two strains Pekin ducks. This study provides a reference for differences between LPD and FPD that manifest as a result of long-term selection.

Genome-wide association study of bone quality and feed efficiency-related traits in Pekin ducks.

Feeding and bone traits are vital for breeding and reproduction in the commercial duck industry. In this study, we performed a genome-wide association study for feeding and bone traits in a population of 540 lean-type Pekin ducks, followed by genotyping-by-sequencing procedures. The genetic parameters of feeding and bone traits were also estimated using genomic information. In total, seventy-eight significant SNPs were determined, and eleven of them reached the genome-wide significant level for 7 traits except for body weight at 42-day old. A peak of genome-wide significant SNPs was detected on chromosome 2 for feed conversion ratio (P-value = 7.46E-11), and the top SNP (P-value = 2.23E-08) for bone-breaking strength was also obtained in the upstream of gene RAPGEF5. This study provided a list of novel markers and candidate genes associated with feeding and bone traits in Pekin ducks, which could contribute to the genetic selection in duck breeding.

Selection response and genetic parameter estimation of feeding behavior traits in Pekin ducks.

Body weight-related traits and feeding behavior traits are important in poultry breeding and production. To investigate the heritability of feeding behavior and their genetic correlation with body weight-related traits in Pekin ducks, 5,594 Pekin ducks were selected. The information about body weight-related traits and feeding behavior from 3 to 6 wk of age were recorded by automatic electronic feeders. The heritability estimates for body weight, residual feed intake, and feeding behavior traits are relatively high (ranging from 0.29 to 0.65). We observed that total feed intake, meal feed intake, body weight at the age of 3 wk, and daily body weight gain had strong positive genetic correlations with body weight at the age of 6 wk. Moreover, body weight at the age of 3 wk also showed a positive genetic correlation with the feed conversion ratio (0.33). Total feeding time, daily feed intake, and feeding rate had significant positive phenotypic correlations with feed efficiency. However, the average interval between meals, the number of daily visits, and the number of meals all had a low genetic or phenotypic relationship with body weight and feed efficiency. In conclusion, our study revealed that body weight, residual feed intake, and feeding behavior traits were all highly heritable traits, and the selection for certain feeding behaviors could improve feed efficiency in Pekin ducks. This study is the first report about genetic parameter estimates about feeding behaviors in ducks based on large datasets and provides solid data for genetic study in ducks.

Genome-wide association study of the level of blood components in Pekin ducks.

Blood components are considered to reflect nutrient metabolism and immune activity in both humans and animals. In this study, we measured 12 blood components in Pekin ducks and performed genome-wide association analysis to identify the QTLs (quantitative trait locus) using a genotyping-by-sequencing strategy. A total of 54 QTLs were identified for blood components. One genome-wide significant QTL for alkaline phosphatase was identified within the intron-region of the OTOG gene (P = 1.31E-07). Moreover, 21 genome-wide significant SNPs for the level of serum cholinesterase were identified on six different scaffolds. In addition, for serum calcium, one genome-wide significant QTL was identified in the upstream region of gene RAB11B. These results provide new markers for functional studies in Pekin ducks, and several candidate genes were identified, which may provide additional insights into specific mechanisms for blood metabolism in ducks and their potential application for duck breeding programs.

Genome-Wide Association Study of Growth and Feeding Traits in Pekin Ducks.

Growth rate and feeding efficiency are the most important economic traits for meat animals. Pekin duck is one of the major global breeds of meat-type duck. This study aims to identify QTL for duck growth and feeding efficiency traits in order to assist artificial selection. In this study, the growth and feeding related phenotypes of 639 Pekin ducks were recorded, and each individual genotype was evaluated using a genotyping-by-sequencing (GBS) protocol. The genetic parameters for growth and feeding efficiency related traits were estimated. Genome-wide association analysis (GWAS) was then performed for these traits. In total, 15 non-overlapping QTLs for the measured traits and 12 significant SNPs for feed efficiency traits were discovered using a mixed linear model. The most significant loci of feed intake (FI) is located in a 182Mb region on Chr1, which is downstream of gene RNF17, and can explain 2.3% of the phenotypic variation. This locus is also significantly associated with residual feed intake (RFI), and can explain 3% of this phenotypic variation. Among 12 SNPs associated with the feed conversion ratio (FCR), the most significant SNP (P-value = 1.65E-06), which was located in the region between the 3rd and 4th exon of the SORCS1 gene on Chr6, explained 3% of the phenotypic variance. Using gene-set analysis, a total of two significant genes were detected be associated with RFI on Chr1. This study is the first GWAS for growth and feeding efficiency related traits in ducks. Our results provide a list of candidate genes for marker assisted selection for growth and feeding efficiency, and also help to better understand the genetic mechanisms of feed efficiency and growth in ducks.

Genome-wide association study reveals novel loci associated with body size and carcass yields in Pekin ducks.

BACKGROUND: Pekin duck products have become popular in Asia over recent decades and account for an increasing market share. However, the genetic mechanisms affecting carcass growth in Pekin ducks remain unknown. This study aimed to identify quantitative trait loci affecting body size and carcass yields in Pekin ducks. RESULTS: We measured 18 carcass traits in 639 Pekin ducks and performed genotyping using genotyping-by-sequencing (GBS). Loci-based association analysis detected 37 significant loci for the 17 traits. Thirty-seven identified candidate genes were involved in many biological processes. One single nucleotide polymorphism (SNP) (Chr1_140105435 A > T) located in the intron of the ATPase phospholipid transporting 11A gene (ATP11A) attained genome-wide significance associated with five weight traits. Eight SNPs were significantly associated with three body size traits, including the candidate gene plexin domain containing 2 (PLXDC2) associated with breast width and tensin 3 (TNS3) associated with fossil bone length. Only two SNPs were significantly associated with foot weight and four SNPs were significantly associated with heart weight. In the gene-based analysis, three genes (LOC101791418, TUBGCP3 (encoding tubulin gamma complex-associated protein 3), and ATP11A) were associated with four traits (42-day body weight, eviscerated weight, half-eviscerated weight, and leg muscle weight percentage). However, no loci were significantly associated with leg muscle weight in this study. CONCLUSIONS: The novel results of this study improve our understanding of the genetic mechanisms regulating body growth in ducks and thus provide a genetic basis for breeding programs aimed at maximizing the economic potential of Pekin ducks.

In vivo prediction of the carcass fatness using live body measurements in Pekin ducks.

The purpose of this study was to evaluate the correlation between live body measurements and several fat traits in Pekin ducks, and ultimately to formulate multiple regression equations for the in vivo estimation of the carcass fatness of Pekin ducks. Several traits were measured in a total of 208 Pekin ducks aged 6 wk (107 males and 101 females). All ducks were weighed and measured for a set of body measurements including live body weight, body slope length, breast muscle thickness, skin fat thickness, chest width, keel length, and neck length. The breast muscle thickness and skin fat thickness was measured using B-scan sonography. Carcass information, including eviscerated weight, subcutaneous fat with skin weight, and abdominal fat weight, was collected after slaughter. Our results revealed that sex effects on most traits were significant (P < 0.05), and that the weight of subcutaneous fat with skin was significantly correlated with live body weight (r = 0.57 to 0.71, P < 0.01). Four additional traits of males were closely correlated with the weight of subcutaneous fat with skin, namely breast muscle thickness (r = 0.20, P < 0.01), skin fat thickness (r = 0.43, P < 0.01), chest width (r = 0.24, P < 0.01), and neck length (r = 0.20, P < 0.05). The abdominal fat weight, percentage of fat, and percentage of subcutaneous fat with skin of ducks were significantly correlated with live body weight (r = 0.38 to 0.43, P < 0.01), and skin fat thickness (r = 0.38 to 0.49, P < 0.01). These traits provided the basis for constructing regression equations to predict weight (or percentage) of subcutaneous fat and abdominal fat with high values of coefficients of multiple correlation (R) between the dependent variable and the independent variables. Two equations were verified to be applicable in other duck groups, with high accuracy, as more than 80% of estimated values were within the margin of error (<10%), compared with the actual values.

De novo transcriptome assembly and identification of genes associated with feed conversion ratio and breast muscle yield in domestic ducks.

Breast muscle yield and feed conversion efficiency are the major breeding aims in duck breeding. Understanding the role of specific transcripts in the muscle and small intestine might lead to the elucidation of interrelated biological processes. In this study, we obtained jejunum and breast muscle samples from two strains of Peking ducks that were sorted by feed conversion ratio (FCR) and breast muscle percentage into two-tailed populations. Ten RNA-Seq libraries were developed from the pooled samples and sequenced using the Hiseq2000 platform. We created a reference duck transcript database using de novo assembly methods, which included 16 663 irredundant contigs with an N50 length of 1530 bp. This new duck reference cDNA dataset significantly improved the mapping rate for RNA-Seq data, from 50% to 70%. Mapping and annotation were followed by Gene Ontology analysis, which showed that numerous genes were differentially expressed between the low and high FCR groups. The differentially expressed genes in the jejunum were enriched in biological processes related to immune response and immune response activation, whereas those in the breast muscle were significantly enriched in biological processes related to muscle cell differentiation and organ development. We identified new candidate genes, that is, PCK1, for improving the FCR and breast muscle yield of ducks and obtained much better reference duck transcripts. This study suggested that de novo assembly is essential when applying transcriptome analysis to a species with an incomplete genome.