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.

Geese feather production and the quality of downy feathers are additional economically important traits in the geese industry. However, little information is available about the molecular mechanisms fundamental to feather formation and the quality of feathers in geese. This study conducted de novo transcriptome sequencing analysis of two related geese species using the Illumina 4000 platform to determine the genes involved in embryonic skin feather follicle development. A total of 165,564,278 for Anser anser and 144,595,262 for Anser cygnoides clean reads were generated, which were further assembled into 77,134 unigenes with an average length of 906 base pairs in Anser anser and 66,041 unigenes with an average length of 922 base pairs in Anser cygnoides. To recognize the potential regulatory roles of differentially expressed genes (DEGs) during geese embryonic skin feather follicle development, the obtained unigenes were annotated to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis. In both species, GO and KOG had shown similar distribution patterns during functional annotation except for KEGG, which showed significant variation in signaling enrichment. Anser asnser was significantly enriched in the calcium signaling pathway, whereas Anser cygnoides was significantly enriched with glycerolipid metabolism. Further analysis indicated that 14,227 gene families were conserved between the species, among which a total of 20,715 specific gene families were identified. Comparative RNA-Seq data analysis may reveal inclusive knowledge to assist in the identification of genetic regulators at a molecular level to improve feather quality production in geese and other poultry species.