A gene co‐expression network analysis of the candidate genes and molecular pathways associated with feather follicle traits of chicken skin

Abstract Back and thigh skin of chickens showed significant differences in the thickness and the feather follicle density and size, which are important traits for slaughtered chickens' appearance. In the present study, global gene expression profiling was conducted in the back and thigh skin of chickens using Microarray technology. The results showed that 676 genes were differentially expressed between back and thigh skin. The expression of the differentially expressed genes (DEGs), including PPP1R3C, IGF1, PTCHD1, HOXB6, FGF9, CAMK4, SHH, BMP8B, FOXN1 and PTGER2, was validated by real‐time quantitative polymerase chain reaction (RT‐qPCR), and the results were consistent with microarray results. Functional analysis revealed that the DEGs were significantly involved in cell proliferation, differentiation, apoptosis, adhesion and transport process, and the pathways were significantly mapped into the ECM‐receptor interaction, peroxisome, focal adhesion, Hedgehog and PPAR signalling pathways. Protein–protein interaction network analysis suggested that signalling pathways related to feathers morphogenesis and development, such as Wnt, FGF, MAPK, SHH and BMP signalling pathways, occupied important positions in the network. Genes involved in these signalling pathways and adhesion molecules might play a vital role in skin and feather follicle development. Further single nucleotide polymorphism (SNP) association analysis of Wnt3A showed that the AC genotype of SNP g.255361 C>A significantly increased the feather follicle density of thigh skin. Our findings may provide new insights on candidate genes and pathways related to skin and feather follicle formation of chickens.