An Analysis of Skin Thickness in the Dezhou Donkey Population and Identification of Candidate Genes by RNA-seq.

The aim of the present study was to analyze the main factors that have a significant impact on skin thickness, and to further identify the genes and signaling pathways regulating skin growth by RNA-seq in Dezhou donkeys. Skin samples from different body regions of 15 slaughtered donkeys were obtained to study variations in skin thickness over the bodies. Skin thickness data for another 514 donkeys was obtained by minimally invasive skin sampling from the back, and measurements of the donkeys' body size traits and pedigree data were also collected. These data were used to analyze changes in skin thickness and estimate genetic parameters. In addition, transcriptomic analysis was conducted on the skin tissues of individuals from two groups with significant differences in skin thickness. Our results showed that skin thickness over the bodies ranged from 1.08 to 4.36 mm. The skin from the back was the thickest and had the highest correlation with that of other regions of the body. The skin thickness decreased from the back to the side of the ventral abdomen, and the skin thickness on the limbs increased from the proximal end to the distal end. The results also showed that the skin from the same body regions of jacks was thicker than that of jennies in the same age group. The skin thickness of jennies increased from birth to the age of 2 and then clearly decreased after 2 years of age. The estimated heritability of skin thickness was 0.15, and the genetic correlations between skin thickness and body size traits were negligible. Transcriptome analysis showed that the thick-skin group had 65 up-regulated genes and 38 down-regulated genes compared with the thin-skin group. The differentially expressed genes were highly enriched in epidermal development and cell adhesion molecule signaling pathways. We identified the candidate genes responsible for variations in skin thickness in the Dezhou donkey, including KRT10, KRT1, CLDN9, MHCII and MMP28. These results contribute to a better understanding of the growth and development of donkey skin, reveal the molecular mechanism responsible for donkey skin thickness and suggest directions for genetic selection in the Dezhou donkey population.