ABSTRACT: The crest trait is a specific and widely distributed phenotype in birds. However, the shape and function vary in different species of birds. To understand the mechanism of crest formation, the present study used RNA sequencing and weighted gene co-expression network analysis (WGCNA) to identify the crest-cushion-associated genes in the Chinese crested (CC) duck. As a result, 28, 40, 32, 33, and 126 differentially expressed genes (DEGs) were identified between CC and cherry valley (CV) ducks at the embryonic days (E)15, E22, E28, D7 (7 days old), and D42 stages, respectively. In addition, the results of WGCNA show that 3697 (turquoise module), 485 (green-yellow module), 687 (brown module), 205 (red module), and 1070 (yellow module) hub genes were identified in the E15, E22, E28, D7, and D42
