The objective of this study was to evaluate the changes in the goose embryo transcriptome during feather development. RNA-Sequencing (RNA-Seq) was used to find the transcriptome profiles of feather follicles from three stages of embryonic dorsal skin at embryonic day 13, 18, and 28 (E13, E18, E28). The results showed that 3001, 6634, and 13,780 genes were differently expressed in three stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed genes (DEGs) in E13 vs. E18 were significantly mapped into the GO term of extracellular structure organization and the pathway of extracellular matrix (ECM)-receptor interaction. In E18 vs. E28, the top significantly mapped into GO term was the single-organism developmental process; the pathway was also the ECM-receptor interaction. DEGs in E13 vs. E28 were significantly mapped into the GO term of the multicellular organismal process and the pathway of cell adhesion molecules. Subsequently, the union of DEGs was categorized by succession cluster into eight profiles, which were then grouped into four ideal profiles. Lastly, the seven genes spatio-temporal expression pattern was confirmed by real-time PCR. Our findings advocate that interleukin 20 receptor subunit alpha (IL20RA), interleukin 6 receptor (IL6R), interleukin 1 receptor type 1 (IL-1R1), Wnt family member 3A (WNT3A), insulin-like growth factor binding protein 3 (IGFBP3), bone morphogenetic protein 7 (BMP7), and secreted-frizzled related protein 2 (SFRP2) might possibly play vital roles in skin and feather follicle development and growth processes.
Keywords: cluster analysis; de novo transcriptome assembly; differentially expressed genes; feather follicle development.