The main component of donkey meat is skeletal muscle, and different muscle fibers have been found to be associated with meat quality. However, RNA-seq technology has yet to be used to profile the transcriptomic changes of different muscles of the donkey. In this study, the characterizations of different muscles on the gene expression profiles of Dezhou donkey were obtained, the aim was to identify the important genes in donkey muscles, and aid in improving donkey meat quality via RNA-seq. In the donkey gluteus (DG) and donkey longissimus dorsi (DL) group, GO enrichment indicated that DEGs were mainly involved in the biological regulation and metabolic process, and KEGG analysis shows that a total of 427 DEGs were mapped to 216 KEGG pathways and 23 KEGG pathways were significantly enriched such as the ribosome, glycolysis/gluconeogenesis, glucagon signaling pathway and biosynthesis of amino acids pathways. Meanwhile, 504 DEGs were mapped to 223 KEGG pathways, in which 17 were significantly enriched including cardiac muscle contraction and oxytocin signaling pathway in donkey hamstring muscles (DH) and DL group. In addition, the tenderness in donkey meat might involve muscle fiber type and glucose metabolism, which might profit from the DEGs including MYH1, MYH7, TNNC1, TNNI3, TPM3, ALDOA, ENO3, and PGK1. The genes found in this study will provide some ideas for further understanding the molecular mechanism of donkey meat quality.

The number of thoracolumbar vertebrae is a quantitative trait positively correlated with the economic traits of livestock. More thoracolumbar vertebrae individuals could genetically be used to improve the livestock population, as more thoracolumbar vertebrae means a longer carcass, which could bring more meat production. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate gene for effecting the number of vertebrae in livestock. The purposes of this study are as follows: (a) Analyzing the effect of TLN variation on body size and carcass traits of Dezhou donkey; (b) Studying the distribution of seven single nucleotide variants (SNVs) in NR6A1 gene of Dezhou donkey; (c) Exploring the relationship between latent SNVs and TLN, the body size and carcass traits. We examined the thoracic and lumbar vertebrae number and seven SNVs in NR6A1 gene of 455 Dezhou donkeys, and analyzed the relationships between them. Five types of thoracolumbar combinations (T17L5 (individual with 17 thoracic and five lumbar vertebrae) 2.4%, T18L5 75.8%, T19L5 1.1%, T17L6 11.9%, and T18L6 8.8%) of Dezhou donkeys were detected in this study. For one thoracolumbar vertebra added, the body length of Dezhou donkey increases by 3 cm and the carcass weight increases by 6 kg. Seven SNVs (g.18093100G > T, g.18094587G > T, g.18106043G > T, g.18108764G > T, g.18110615T > G, g.18112000C > T and g.18114954T > G) of the NR6A1 gene were found to have a significant association with the TLN, body size and carcass traits of Dezhou donkey (p < 0.05), respectively. For instance, g.18114954C > T is significantly associated with lumber vertebrae number, the total number of thoracolumbar, and carcass weight, and individuals with TT genotype had significantly larger value than CC genotype (p < 0.05). Using these 7SNVs, 16 different haplotypes were estimated. Compared to Hap3Hap3, individuals homozygous for Hap2Hap2 showed significantly longer length in one thoracic spine (STL), the total thoracic vertebrae and one thoracolumbar spine. Our study will not only extend the understanding of genetic variation in the NR6A1 gene of Dezhou donkey, but also provide useful information for marker assisted selection in donkey breeding program.

Several studies have shown the association between the ligand-dependent nuclear receptor compression-like protein (LCORL) gene and body size in horses, pigs and donkeys. Based on previous studies, the LCORL gene was hypothesized to be associated with growth traits and hide weight in Dezhou donkeys. In this study, we aimed to reveal the variation of the LCORL gene in the Dezhou donkey and explore whether the gene is associated with hide weight and body size. In this study, genetic polymorphisms in the LCORL gene of the Dezhou donkey were studied using targeted sequencing technology, and single nucleotide polymorphisms (SNPs) of the LCORL gene were analyzed for association with hide weight and body size in Dezhou donkeys. The results showed that there was an SNP locus situated in intron 1 of the LCORL gene. Association analysis revealed that individuals with the GG genotype had significantly higher body height, body length, chest circumference and hide weight than those with the AA genotype (p < 0.05). Therefore, the g.112558859 A > G locus can be used as a potential candidate marker affecting body size and hide weight. This study provides the foundation for breeding high-quality donkeys with high hide yield.

An increase in the number of vertebrae can significantly affect the meat production performance of livestock, thus increasing carcass weight, which is of great importance for livestock production. The homeobox gene C8 (HOXC8) has been identified as an essential candidate gene for regulating vertebral development. However, it has not been researched on the Dezhou donkey. This study aimed to verify the Dezhou donkey HOXC8 gene’s polymorphisms and assess their effects on multiple vertebral numbers and carcass weight. In this study, the entire HOXC8 gene of the Dezhou donkey was sequenced, SNPs at the whole gene level were identified, and typing was accomplished utilizing a targeted sequencing genotype detection technique (GBTS). Then, a general linear model was used to perform an association study of HOXC8 gene polymorphism loci, multiple vertebral numbers, and carcass weight for screening candidate markers that can be used for molecular breeding of Dezhou donkeys. These findings revealed that HOXC8 included 12 SNPs, all unique mutant loci. The HOXC8 g.15179224C>T was significantly negatively associated with carcass weight (CW) and lumbar vertebrae length (LL) (p < 0.05). The g.15179674G>A locus was shown to be significantly positively associated with the number of lumbar vertebrae (LN) (p < 0.05). The phylogenetic tree constructed for the Dezhou donkey HOXC8 gene and seven other species revealed that the HOXC8 gene was highly conserved during animal evolution but differed markedly among distantly related animals. The results suggest that HOXC8 is a vital gene affecting multiple vertebral numbers and carcass weight in Dezhou donkeys, and the two loci g.15179224C>T and g.15179674G>A may be potential genetic markers for screening and breeding of new strains of high-quality and high-yielding Dezhou donkeys.