A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes.

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

Selection on the promoter regions plays an important role in complex traits during duck domestication.

BACKGROUND: Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood. RESULTS: We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment. CONCLUSIONS: This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

Genome-wide association study identified the candidate genes associated with angel wing trait in Pekin duck.

Angel wing is a developmental wing deformity that can influence breeding and reproduction in the commercial duck industry. The nutrition foundation of angel wing trait was initially explored, but the genetic basic remains poorly understood. In this study, we identified candidate genes and single-nucleotide polymorphisms (SNPs) associated with angel wing trait in Pekin ducks using a genome-wide association study (GWAS) and selective sweep analysis. The GWAS results showed that nine SNPs across five chromosomes were significantly correlated with the angel wing trait. In total, 468 selection signals were shown between the angel wing ducks and normal ducks, and these signals harbored 154 genes, which were enriched in the nervous system and metabolism. This study provides the new insights into the genetic factors that may influence duck angel wing.

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.

Integrating transcriptome, proteome and QTL data to discover functionally important genes for duck eggshell and albumen formation.

Duck egg quality improvement is an essential target for Asian poultry breeding. In total, 15 RNA-Seq libraries (magnum, isthmus, and uterus at two different physiological states) were sequenced from 48 weeks old Pekin ducks. De novo assembly and annotation methods were utilized to generate new reference transcripts. Our results revealed that 1264 and 2517 genes were differentially expressed in magnum and uterus in the presence versus absence of an egg, respectively. We identified 1089 genes that were differentially expressed in isthmus compared to uterus (in both presence and absence of a calcifying egg). We observed that 11 common DEGs were detected in the egg white proteomes of 6 different bird species including domestic Chicken, Duck, Goose, Turkey, Quail, and Pigeon. On the other hand, only one of the top five most highly expressed genes in duck isthmus was in this category for the chicken isthmus (SPINK7). Among the large number of DEGs during eggshell formation in ducks, only 41 genes showed a similar differential expression pattern in both duck and chicken. By combining chicken QTL database, chicken oviduct transcriptome and egg proteome data for five bird species, we have obtained high-quality gene lists for egg formation. This is the first study to elucidate the transcriptomic changes in different duck oviduct segments during egg formation, and to integrate QTL, proteome and transcriptome data to probe the functional genes associated with albumen secretion and eggshell mineralization.

Dynamics of transcriptome changes during subcutaneous preadipocyte differentiation in ducks.

BACKGROUND: Pekin duck is an important animal model for its ability for fat synthesis and deposition. However, transcriptional dynamic regulation of adipose differentiation driven by complex signal cascades remains largely unexplored in this model. This study aimed to explore adipogenic transcriptional dynamics before (proliferation) and after (differentiation) initial preadipocyte differentiation in ducks. RESULTS: Exogenous oleic acid alone successfully induced duck subcutaneous preadipocyte differentiation. We explored 36 mRNA-seq libraries in order to study transcriptome dynamics during proliferation and differentiation processes at 6 time points. Using robust statistical analysis, we identified 845, 652, 359, 2401 and 1933 genes differentially expressed between -48 h and 0 h, 0 h and 12 h, 12 h and 24 h, 24 h and 48 h, 48 h and 72 h, respectively (FDR < 0.05, FC > 1.5). At the proliferation stage, proliferation related pathways and basic cellular and metabolic processes were inhibited, while regulatory factors that initiate differentiation enter the ready-to-activate state, which provides a precondition for initiating adipose differentiation. According to weighted gene co-expression network analysis, pathways positively related to adipogenic differentiation are significantly activated at the differentiation stage, while WNT, FOXO and other pathways that inhibit preadipocyte differentiation are negatively regulated. Moreover, we identified and classified more than 100 transcription factors that showed significant changes during differentiation, and found novel transcription factors that were not reported to be related to preadipoctye differentiation. Finally, we manually assembled a proposed regulation network model of subcutaneous preadipocyte differentiation base on the expression data, and suggested that E2F1 may serve as an important link between the processes of duck subcutaneous preadipocyte proliferation and differentiation. CONCLUSIONS: For the first time we comprehensively analyzed the transcriptome dynamics of duck subcutaneous preadipocyte proliferation and differentiation. The current study provides a solid basis for understanding the synthesis and deposition of subcutaneous fat in ducks. Furthermore, the information generated will allow future investigations of specific genes involved in particular stages of duck adipogenesis.

iTRAQ-Based Quantitative Proteomic Analysis of Duck Eggshell During Biomineralization.

The avian egg is a valuable model for the calcitic biomineralization process as it is the fastest calcification process occurring in nature and is a clear example of biomineralization. In this study, iTRAQ MS/MS is used to detect and study for the first time: 1) the overall duck eggshell proteome; 2) regional differences in the proteome between the inner and outer portions of the duck eggshell. The new reference protein datasets allow us to identify 179 more eggshell proteins than solely using the current release of Ensembl duck annotations. In total, 484 proteins are identified in the entire duck eggshell proteome. Twenty-eight novel proteins of unknown function that are involved in eggshell formation are also identified. Among the identified eggshell proteins, 54 proteins show differential abundances between the inner, partially mineralized eggshell (obtained 16 h after ovulation) compared to the overall complete eggshell (normally expulsed eggshell). At least 64 of the abundant matrix proteins are common to eggshell of 4 different domesticated bird species (chicken, duck, quail, turkey) and zebra finch. This study provides a new resource for avian eggshell proteomics, and augments the inventory of eggshell matrix proteins that will lead to a deeper understanding of calcitic biomineralization.

Revisiting avian 'missing' genes from de novo assembled transcripts.

BACKGROUND: Argument remains as to whether birds have lost genes compared with mammals and non-avian vertebrates during speciation. High quality-reference gene sets are necessary for precisely evaluating gene gain and loss. It is essential to explore new reference transcripts from large-scale de novo assembled transcriptomes to recover the potential hidden genes in avian genomes. RESULTS: We explored 196 high quality transcriptomic datasets from five bird species to reconstruct transcripts for the purpose of discovering potential hidden genes in the avian genomes. We constructed a relatively complete and high-quality bird transcript database (1,623,045 transcripts after quality control in five birds) from a large amount of avian transcriptomic data, and found most of the presumed missing genes (83.2%) could be recovered in at least one bird species. Most of these genes have been identified for the first time in birds. Our results demonstrate that 67.94% genes have GC content over 50%, while 2.91% genes are AT-rich (AT% > 60%). In our results, 239 (53.59%) genes had a tissue-specific expression index of more than 0.9 in chicken. The missing genes also have lower Ka/Ks values than average (genome-wide: Ka/Ks = 0.99; missing gene: Ka/Ks = 0.90; t-test = 1.25E-14). Among all presumed missing genes, there were 135 for which we did not find any meaningful orthologues in any of the 5 species studied. CONCLUSION: Insufficient reference genome quality is the major reason for wrongly inferring missing genes in birds. Those presumably missing genes often have a very strong tissue-specific expression pattern. We show multi-tissue transcriptomic data from various species are necessary for inferring gene family evolution for species with only draft reference genomes.