## 3 Results
### Genetic Polymorphism of NR6A1 Gene in Dezhou Donkey  
The _NR6A1_ gene of donkey was located in chromosome 10, consists of nine exons and 10 introns. Seven polymorphism sites including g.18093100G \(>\) T, g.18094587C \(>\) G, g.18106043A \(>\) C, g.18108764A \(>\) G, g.18110615G \(>\) C, g.18112000T \(>\) G and g.18114954C \(>\) T were first identified in the introns of the _NR6A1_ gene of Dezhou donkey by direct sequencing and comparing these samples with the reference sequence (GenBank accession no. NC_052186.1 Figure 1). The genetic information for the seven loci is shown in Table 3. In the Hardy-Weinberg equilibrium test, six loci were in the Hardy-Weinberg equilibrium (\(p>0.05\)). The genetic PIC of _NR6A1_ mutation sites ranged from 0.417 to 0.506, which reflects that the genetic diversity of _NR6A1_ gene in Dezhou donkeys is moderate. \begin{table}
\begin{tabular}{c c c c c c c} \hline
**Type** & **N** & **BL** & **CHC** & **CW** & **TL** & **LL** & **TLL** \\ \hline  & & **(Mean \(\pm\) SE, cm)** & **(Mean \(\pm\) SE, cm)** & **(Mean \(\pm\) SE, cm)** & **(Mean \(\pm\) SE, cm)** & **(Mean \(\pm\) SE, cm)** \\ \hline
17 & 65 & 131.33 \(\pm\) 0.6 & 141.74 \(\pm\) 2.29 & 148.9 \(\pm\) 1.71 & 68.86 \(\pm\) 0.35 & 26.27 \(\pm\) 0.24 & 95.13 \(\pm\) 0.47 \\
18 & 385 & 132.66 \(\pm\) 0.32 & 144.66 \(\pm\) 0.28 & 152.36 \(\pm\) 0.86 & 73.65 \(\pm\) 0.17 & 23.78 \(\pm\) 0.1 & 97.43 \(\pm\) 0.23 \\
19 & 5 & 132.2 \(\pm\) 3.85 & 144.4 \(\pm\) 2.5 & 151.5 \(\pm\) 6.26 & 75.2 \(\pm\) 1.59 & 23.2 \(\pm\) 0.37 & 98.4 \(\pm\) 1.75 \\ Pearson correlation & 0.070 & 0.113 ** & 0.049 & 0.099 ** & \(-\)0.171 *** & 0.012 \\ \(p\) value & 0.137 & 0.016 & 0.297 & 0.035 & 0.000 & 0.796 \\
5 & 361 & 132.28 \(\pm\) 0.33 & 144.54 \(\pm\) 0.29 & 151.89 \(\pm\) 0.91 & 73.53 \(\pm\) 0.18 & 23.31 \(\pm\) 0.08 & 96.86 \(\pm\) 0.23 \\
6 & 94 & 133.2 \(\pm\) 0.56 & 143.1 \(\pm\) 1.59 & 151.74 \(\pm\) 1.41 & 70.9 \(\pm\) 0.39 & 27.23 \(\pm\) 0.14 & 98.05 \(\pm\) 0.48 \\ Pearson correlation & 0.060 & \(-\)0.068 & 0.040 & \(-\)0.024 & 0.312 ** & 0.060 \\ \(p\) value & 0.199 & 0.148 & 0.390 & 0.603 & 0.000 & 0.200 \\
22 & 11 & 128.64 \(\pm\) 1.44 & 141.18 \(\pm\) 1.41 & 142.3 \(\pm\) 3.6 & 68.55 \(\pm\) 0.79 & 23.05 \(\pm\) 0.38 & 91.59 \(\pm\) 0.96 \\
23 & 399 & 132.32 \(\pm\) 0.31 & 144.27 \(\pm\) 0.45 & 151.92 \(\pm\) 0.85 & 73.03 \(\pm\) 0.18 & 23.82 \(\pm\) 0.1 & 96.86 \(\pm\) 0.21 \\
24 & 45 & 134.67 \(\pm\) 0.92 & 144.73 \(\pm\) 0.57 & 153.52 \(\pm\) 1.99 & 73.65 \(\pm\) 0.5 & 27.09 \(\pm\) 0.31 & 100.66 \(\pm\) 0.66 \\ Pearson correlation & 0.146 *** & 0.042 & 0.101 ** & 0.077 * & 0.184 *** & 0.084 * \\ \(p\) value & 0.002 & 0.376 & 0.032 & 0.099 & 0.000 & 0.073 \\ \hline \end{tabular} Note: *** means \(p<0.01\), Component differences reached a very significant level; ** means \(0.01<p<0.05\), Component differences reached significant level; ** means \(0.05<p<0.1\), There was a significant difference trend between the groups. BL means body length, CHC means chest circumference, CW means carcass weight, TL means the length of thoracic vertebrae, LL means the length of lumbar vertebrae, TLL means the total length of thoracolumbar vertebrae.  
\end{table} 
Table 2: The association analysis of TN, LN, and TLN with body size and carcass traits of Dezhou donkeys. Figure 1: Seven SNVs in the _NR6A1_ gene of the Dezhou donkey. Note: The red dotted line indicates the peak map of relevant SNV sites.