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.

The Full-Length Transcriptome Provides New Insights Into the Transcript Complexity of Abdominal Adipose and Subcutaneous Adipose in Pekin Ducks.

Adipose tissues have a central role in organisms, and adipose content is a crucial economic trait of poultry. Pekin duck is an ideal model to study the mechanism of abdominal and subcutaneous adipose deposition for its high ability of adipose synthesis and deposition. Alternative splicing contributes to functional diversity in abdominal and subcutaneous adipose. However, there has been no systematic analysis of the dynamics of differential alternative splicing of abdominal and subcutaneous adipose in Pekin duck. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to explore the transcriptional complexity of abdominal and subcutaneous adipose in Pekin ducks. In total, 143,931 and 111,337 full-length non-chimeric transcriptome sequences of abdominal and subcutaneous adipocytes were obtained from 41.78 GB raw data, respectively. These data led us to identify 19,212 long non-coding RNAs (lncRNAs) and 74,571 alternative splicing events. In addition, combined with the next-generation sequencing technology, we correlated the structure and function annotation with the differential expression profiles of abdominal and subcutaneous adipose transcripts. This study identified lots of novel alternative splicing events and major transcripts of transcription factors related to adipose synthesis. STAT3 was reported as a vital gene for adipogenesis, and we found that its major transcript is STAT3-1, which may play a considerable role in the process of adipose synthesis in Pekin duck. This study greatly increases our understanding of the gene models, genome annotations, genome structures, and the complexity and diversity of abdominal and subcutaneous adipose in Pekin duck. These data provide insights into the regulation of alternative splicing events, which form an essential part of transcript diversity during adipogenesis in poultry. The results of this study provide an invaluable resource for studying alternative splicing and tissue-specific expression.

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.

MYOD1 inhibits avian adipocyte differentiation via miRNA-206/KLF4 axis.

BACKGROUND: A considerable number of muscle development-related genes were differentially expressed in the early stage of avian adipocyte differentiation. However, the functions of them in adipocyte differentiation remain largely known. In this study, the myoblast determination protein 1 (MYOD1) was selected as a representative of muscle development. We investigated its expression, function, and regulation in avian adipocyte differentiation. RESULTS: The expression of MYOD1 decreased significantly in the early stage of avian adipocyte differentiation. CRISPR/Cas9-mediated deletion of MYOD1 induced adipocyte differentiation, whereas over-expression of MYOD1 inhibited adipogenesis. The mRNA-seq data showed that MYOD1 could perturb the lipid biosynthetic process during differentiation. Our results showed that MYOD1 directly up-regulates the miR-206 expression by binding the upstream 1200 bp region of miR-206. Then, over-expression of miR-206 can inhibit the adipogenesis. Furthermore, MYOD1 affected the expression of endogenous miR-206 and its target gene Kruppel-like factor 4 (KLF4), which is an important activator of adipogenesis. Accordingly, the inhibition of miR-206 or over-expression of KLF4 could counteract the inhibitory effect of MYOD1 on adipocyte differentiation. CONCLUSIONS: Our results establish that MYOD1 inhibits adipocyte differentiation by up-regulating miR-206 to suppress the KLF4 expression. These findings identify a novel function of MYOD1 in adipocyte differentiation, suggesting a potential role in body-fat distribution regulation.

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.

The transcriptome landscapes of ovary and three oviduct segments during chicken (Gallus gallus) egg formation.

The avian embryo develops within a specialized biological container (eggshell) that contains crucial nutritional compartments (albumen, yolk). We analyzed the transcriptome of ovary and three segments of oviduct, including magnum, isthmus and uterus in the chicken during egg formation. RNA-Seq libraries (42 in total) for ovary and three different parts of the oviduct were sequenced for two different phases of egg formation. We obtained 8365 novel transcripts with an mRNA length longer than 200 bp; of these, 6832 were long intergenic non-coding RNA transcripts. We identified 547 differentially expressed genes in magnum (actively secreting albumen versus inactive) and 585 in uterus (active eggshell calcification versus quiescent). By combining QTL, transcriptome and proteome data, we obtained high quality gene lists for chicken egg formation. This is the first study to describe the ovary and oviduct transcriptomes by mRNA sequencing, and to elucidate the global repertoire of functional genes involved in egg formation.

Full-length transcriptome sequencing from multiple tissues of duck, Anas platyrhynchos.

Duck (Anas platyrhynchos), one of the most economically important waterfowl, is an ideal model for studying the immune protection mechanism of birds. An incomplete duck reference genome and very limited availability of full-length cDNAs has hindered the identification of alternatively spliced transcripts and slowed down many basic studies in ducks. We applied PacBio Iso-Seq technologies to multiple tissues from duck for use in transcriptome sequencing. We obtained 199,993 full-length transcripts and comprehensively annotated these transcripts. 23,755 lncRNAs were predicted from all identified transcripts and 35,031 alternative splicing events, which divided into 5 models, were accurately predicted from 3,346 genes. Our data constitute a large increase in the known number of both lncRNA, and alternatively spliced transcripts of duck and plays an important role in improving current genome annotation. In addition, the data will be extremely useful for functional studies in other birds.

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.