Elucidating the Role of Transcriptomic Networks and DNA Methylation in Collagen Deposition of Dezhou Donkey Skin

DNA methylation represents a predominant epigenetic modification with broad implications in various biological functions. Its role is particularly significant in the process of collagen deposition, a fundamental aspect of dermal development in donkeys. Despite its critical involvement, the mechanistic insights into how DNA methylation influences collagen deposition in donkey skin remain limited. In this study, we employed whole genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) to investigate the epigenetic landscape and gene expression profiles in the dorsal skin tissues of Dezhou donkeys across three developmental stages: embryonic (YD), juvenile (2-year-old, MD), and mature (8-year-old, OD). Our analysis identified numerous differentially methylated genes that play pivotal roles in skin collagen deposition and overall skin maturation, including but not limited to COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, GLUL, SFRP2, FOSL1, SERPINE1, MMP1, MMP2, MMP9, and MMP13. Notably, we observed an inverse relationship between gene expression and DNA methylation proximal to transcription start sites (TSSs), whereas a direct correlation was detected in regions close to transcription termination sites (TTSs). Detailed bisulfite sequencing analyses of the COL1A1 promoter region revealed a low methylation status during the embryonic stage, correlating with elevated transcriptional activity and gene expression levels. Collectively, our findings elucidate key genetic markers associated with collagen deposition in the skin of Dezhou donkeys, underscoring the significant regulatory role of DNA methylation. This research work contributes to the foundational knowledge necessary for the genetic improvement and selective breeding of Dezhou donkeys, aiming to enhance skin quality attributes.

Comprehensive transcriptomic analysis unveils the interplay of mRNA and LncRNA expression in shaping collagen organization and skin development in Dezhou donkeys

The primary focus of donkey hide gelatin processing lies in the dermal layer of donkey hide due to its abundant collagen content. However, the molecular mechanism involved in collagen organization and skin development in donkey skin tissue across various developmental stages remains incomplete. The current study aims to investigate the transcriptomic screening of lncRNAs and mRNA associated with skin development and collagen organization across different ages in Dezhou donkeys' skin. In the pursuit of this objective, we used nine skin tissue samples obtained from Dezhou donkeys at various ages including 8-month fetal stage, followed by 2 and 8 years. RNA-seq analysis was performed for the transcriptomic profiling of differentially expressed genes (DEGs) and lncRNAs associated with skin development in different age groups. Our investigation revealed the presence of 6,582, 6,455, and 405 differentially expressed genes and 654, 789, and 29 differentially expressed LncRNAs within the skin tissues of Dezhou donkeys when comparing young donkeys (YD) vs. middle-aged donkeys (MD), YD vs. old donkeys (OD), and MD vs. OD, respectively. Furthermore, we identified Collagen Type I Alpha 1 Chain (COL1A1), Collagen Type III Alpha 1 Chain (COL3A1), and Collagen Type VI Alpha 5 Chain (COL6A5) as key genes involved in collagen synthesis, with COL1A1 being subject to cis-regulation by several differentially expressed LncRNAs, including ENSEAST00005041187, ENSEAST00005038497, and MSTRG.17248.1, among others. Interestingly, collagen organizational and skin development linked pathways including Protein digestion and absorption, metabolic pathways, Phosphatidylinositol 3-Kinase-Protein Kinase B signaling pathway (PI3K-Akt signaling pathway), Extracellular Matrix-Receptor Interaction (ECM-receptor interaction), and Relaxin signaling were also reported across different age groups in Dezhou donkey skin. These findings enhance our comprehension of the molecular mechanisms underlying Dezhou donkey skin development and collagen biosynthesis and organization, thus furnishing a solid theoretical foundation for future research endeavors in this domain.

RNA-seq analysis identifies differentially expressed gene in different types of donkey skeletal muscles

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.

A novel SNP in NKX1-2 gene is associated with carcass traits in Dezhou donkey

BACKGROUND

At present, donkey meat in the market shows an imbalance between supply and demand, and there is an urgent need to cultivate a meat-type Dezhou donkey breed. On the one hand, it can improve the imbalance in the market, and on the other hand, it can promote the rapid development of the donkey industry. This study aimed to reveal significant genetic variation in the NK1 homeobox 2 gene (NKX1-2) of Dezhou donkeys and investigate the association between genotype and body size in Dezhou donkeys.

RESULTS

In this study, a SNP (g.54704925 A > G) was identified at the exon4 by high-depth resequencing of the Dezhou donkey NKX1-2 gene. The AA genotype is the dominant genotype. The g.54704925 A > G site was significantly associated with body length, thoracic girth, and hide weight (P < 0.05), while it was highly significantly associated with body height and carcass weight (P < 0.01) in Dezhou donkeys.

CONCLUSION

Overall, the results of this study showed that the NKX1-2 gene could be a candidate gene for breeding meat-type Dezhou donkeys, and the g.54704925 A > G locus could be used as a marker locus for selection and breeding.

Polymorphism detection of PRKG2 gene and its association with the number of thoracolumbar vertebrae and carcass traits in Dezhou donkey

BACKGROUND

Previous studies have shown that the protein kinase cGMP-dependent 2 (PRKG2) gene is associated with dwarfism in humans, dogo Argentines, and Angus cattle, as well as with height and osteoblastogenesis in humans. Therefore, the PRKG2 gene was used as the target gene to explore whether this gene is associated with several thoracolumbar vertebrae and carcass traits in Dezhou donkeys.

RESULTS

In this study, fifteen SNPs were identified by targeted sequencing, all of which were located in introns of the PRKG2 gene. Association analysis illustrated that the g.162153251 G > A, g.162156524 C > T, g.162158453 C > T and, g.162163775 T > G were significantly different from carcass weight. g.162166224 G > A, g.162166654 T > A, g.162167165 C > A, g.162167314 A > C and, g.162172653 G > C were significantly associated with the number of thoracic vertebrae. g.162140112 A > G was significantly associated with the number and the length of lumbar vertebrae, and g.162163775 T > G was significantly associated with the total number of thoracolumbar vertebrae.

CONCLUSION

Overall, the results of this study suggest that PRKG2 gene polymorphism can be used as a molecular marker to breed high-quality Dezhou donkeys.

Identification of LTBP2 gene polymorphisms and their association with thoracolumbar vertebrae number, body size, and carcass traits in Dezhou donkeys

The number of thoracolumbar vertebrae in Dezhou donkeys varies from 22 to 24 and is associated with body size and carcass traits. In mammals, the latent transforming growth factor beta binding protein 2 (LTBP2) has been found to have some functions in the development of thoracolumbar vertebrae. The relationship between LTBP2 and TLN (the number of thoracolumbar vertebrae) of Dezhou donkeys is yet to be reported. The purposes of this study are as follows: 1) to quantify the effect of thoracolumbar vertebrae number variation of Dezhou donkeys on body size and carcass trait; 2) to study the distribution of single nucleotide variants (SNVs) in the LTBP2 gene of Dezhou donkeys; and 3) to explore whether these SNVs can be used as candidate sites to study the mechanism of Dezhou donkey muti-thoracolumbar vertebrae development. The TLN, body size, and carcass traits of 392 individuals from a Dezhou donkey breed were recorded. All animals were sequenced for LTBP2 using GBTS liquid chip and 16 SNVs were used for further analysis. We then analyzed the relationship between these SNVs with TLN, body size, and carcass traits. The results showed that: 1) c.5547 + 860 C > T, c.5251 + 281 A > C, c.3769 + 40 C > T, and c.2782 + 3975 A > G were complete genetic linkages and significantly associated with thoracic vertebrae number (TN) (p < 0.05) (wild-type homozygotes had more TN than heterozygotes); 2) c.1381 + 768 T > G and c.1381 + 763 G > T were significantly associated with lumber vertebrae number (LN) (p < 0.05); 3) c.1003 + 704 C > T, c.1003 + 651 C > T, c.1003 + 626 A > G, and c.812 + 22526 T > G were significantly associated with chest circumference (CHC), front carcass weight (CWF), after carcass weight (CWA), and carcass weight (CW) (p < 0.05) (wild-type homozygotes were larger than other genotypes in CHC, CWF, CWA, and CW); and 4) the effect of variation is not consistent in c.565 + 11921 A > G, c.565 + 6840 A > G, c.565 + 3453 C > T, and c.494 + 5808 C > T. These results provide useful information that the polymorphism of LTBP2 is significantly associated with TLN, body size, and carcass traits in Dezhou donkeys, which can serve as a molecule marker to improve donkey production performance.

Association of HOXC8 Genetic Polymorphisms with Multi-Vertebral Number and Carcass Weight in Dezhou Donkey

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.

Identification of Candidate Genes for Twinning Births in Dezhou Donkeys by Detecting Signatures of Selection in Genomic Data

Twinning trait in donkeys is an important manifestation of high fecundity, but few reports are available elucidating its genetic mechanism. To explore the genetic mechanism underlying the twin colt trait in Dezhou donkeys, DNA from 21 female Dezhou donkeys that had birthed single or twin colts were collected for whole-genome resequencing. FST, θπ and Tajima’s D were used to detect the selective sweeps between single and twin colt fecundity in the Dezhou donkey groups. Another set of 20 female Dezhou donkeys with single or multiple follicles during estrus were selected to compare concentrations of reproductive hormone including follicle-stimulating hormone (FSH), luteinizing hormone (LH), estradiol (E2) and progesterone (P4). Four candidate genes including ENO2, PTPN11, SOD2 and CD44 were identified in the present study. The CD44 gene had the highest FST value, and ENO2, PTPN11 and SOD2 were screened by two joint analyses (FST and θπ, θπ and Tajima’s D). There was no significant difference in the LH, FSH and P4 levels between the two groups (p > 0.05); however, the serum E2 content in the multi-follicle group was significantly higher than that in the single-follicle group (p < 0.05). The identified candidate genes may provide new insights into the genetic mechanism of donkey prolificacy and may be useful targets for further research on high reproductive efficiency.

A Novel A > G Polymorphism in the Intron 1 of LCORL Gene Is Significantly Associated with Hide Weight and Body Size in Dezhou Donkey

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.

Multi-Thoracolumbar Variations and NR6A1 Gene Polymorphisms Potentially Associated with Body Size and Carcass Traits of Dezhou Donkey

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.

A novel 31bp deletion within the CDKL5 gene is significantly associated with growth traits in Dezhou donkey

The discovery of molecular markers which associate with livestock economic traits is of great significance for livestock breeding. Selective analysis has found a potential correlation between CDKL5 and growth traits, but there is still a lack of experimental proof. In this study, a 31-bp deletion (g.176595_176626delATGTCACATGTGGTACTGCCATGTGGAATTT) of CDKL5 gene was found by sequencing. The 31-bp indel was then genotyped in 380 individuals of Dezhou donkeys by polyacrylamide gel electrophoresis and there were three genotypes in this population. After the association analysis between growth traits and genotypes, it was found that this 31-bp indel polymorphism was significantly associated with the chest circumference of Dezhou donkeys (p < 0.05), and body length, chest depth and rump width (p < 0.01). In addition, all individuals with DD genotype were better than those with other genotypes in growth traits. This study revealed that a newly identified polymorphic locus in the CDKL5 gene is related to growth traits, which provides a molecular marker for genetic improvement of Dezhou donkey and may lay a solid foundation for the breeding of Dezhou donkey.

Estimation of Liveweight from Body Measurements through Best Fitted Regression Model in Dezhou Donkey Breed

Dezhou donkey is an excellent local hide, meat and milk breed in Shandong Province. To accelerate the speed of breeding, reproduction and conversation, correlation and multiple regression analysis between body weight (BW) and body size of Dezhou donkey rearing under intensive farms was made by SAS 9.4 software (Statistical Analysis for Windows). A total of 162 donkeys of both gender of age 2 ~ 10 years old were used to investigate the relationships between BW and body dimensional traits (cm) including height of withers, body length; thoracic depth, thoracic girth (TG), thoracic width (TW), circumference of cannon bone (CB), height of rump, rump length (RL) and rump width (RW). The results showed that BW and body measurements have positive and great correlations with R² value ranged from 0.58 to 0.88 (P <0.01). The R² values from the single-parameter equations showed that the TG was highly related to BW (0.72, P < 0.01). The stepwise regression equations were applied to obtain the best prediction equations, and the results indicated that the prediction accuracy for BW was improved with the inclusion of more body measurement variables. The "best fit" models were: BW (kg) = 1.88 × TG + 1.27 × BL + 2.55 × TW + 4.61 × CB + 2.18 × RW + 1.78 × RL - 422.8 (R² = 0.906, P < .01). The predicted BW from the present equations showed the nearest value to the real BW (R² > 0.94, P < .01). In addition, the equations derived to predict the BW of donkeys in Britain, Morocco were less satisfactory for use with the present Dezhou donkey breed because they overestimated or underestimated the BW due to the different donkey breed.

Daily Sperm Output, Spermatogenic Efficiency, and Sexual Behavior of Dezhou Donkey Jacks Mounting Jennies in Estrus

This study aimed to assess the sexual behavior of jacks mounting jennies in estrus and determine the daily sperm output (DSO) and spermatogenic efficiency using two equations to calculate testicular volume (TV). Eight sexually rested mature jacks, aging 5 to 10 years old, had semen collected once a day for 10 consecutive days using jennies in good standing estrus for mounting. Sexual behavior and semen parameters were assessed during each collection. Testicular measurements of height, width, and length were taken immediately before the first semen collection, and these measurements were used to calculate TV. After that, the TV was used to predict the DSO. The average total sperm number (TSN) obtained on days 8 to 10 was deemed the actual DSO. Differences in the predicted vs. the actual DSO were used to calculate the spermatogenic efficiency. In addition, the actual DSO was also used to calculate the number of inseminating doses a jack could produce for both on- and off-site breeding. Jack's sexual behavior and sperm motility did not vary across collection days. Sperm concentration and TSN reduced over time (P < .05). The actual DSO was 9.1 ± 4.1 billion, and the predicted DSO varied from 4.7 to 18 billion. Spermatogenic efficiency varied from 50 to 150% based on jack and the equation used to calculate TV. The number of inseminating doses ranged from 15 to 47 at 300-500 million progressively motile sperm (pms)/dose for on-site breeding. In contrast, the number of breeding doses with cooled-shipped semen (1 billion pms/dose) varied from 4 to 14 doses across donkeys. In conclusion, sexual behavior was not affected by daily semen collections. Sexual rest did not affect sperm motility. The predicted DSO varied with the equation used to determine TV. Clinically normal donkeys have high spermatogenic efficiency, which confirms previous histology reports.

Tissue expression profile, polymorphism of IGF1 gene and its effect on body size traits of Dezhou donkey

Insulin-like growth factor 1 (IGF1) is a multifunctional cell proliferation regulator that plays a critical role in regulating animal growth and development. In this study, the expression level of IGF1 gene in different tissues of Dezhou donkey in different periods was investigated by RT-qPCR. Meanwhile, two mutation sites were identified within the IGF1 gene and its effect on body size traits of Dezhou donkey was analysed. The results showed that the expression level of the adult donkey IGF1 gene in heart, liver, spleen, lung, renal and gastric tissues is higher than that of the young donkeys, but the young donkeys are significantly higher in muscle tissues than the adult donkeys. The IGF1-1 and IGF1-2 loci showed a trend that the GG mutant was larger than other genotypes in the growth traits of both male and female donkeys, among which the IGF1-1 loci had a significant association with the chest circumference and chest depth of male donkeys (P < 0.05), and the IGF1-2 loci had a significant association with the chest circumference of female donkeys. Haplotype combination Hap1Hap1 (GG-GG) showed a greater tendency than Hap2Hap2 (AA-GG) combination in terms of growth traits, reflecting that the results were consistent with the analysis results of genotypes, which also proved the analysis results of genotypes and growth traits had certain reliability. In summary, the IGF1 gene is a candidate gene for growth and development, and its polymorphisms can be used as the molecular markers for Dezhou donkey breeding.

Identification and Comparative Analysis of Long Non-Coding RNA in the Skeletal Muscle of Two Dezhou Donkey Strains

Long non-coding RNA (lncRNA) has been extensively studied in many livestock. However, compared with other livestock breeds, there is less research regarding donkey lncRNA function. It has been reported that lncRNA plays an important role in the timing control of development, aging, and death of livestock. Therefore, the study of donkey skeletal muscle lncRNA is of great significance for exploring donkey meat production performance. In this study, RNA-Seq was used to perform high-throughput sequencing of skeletal muscle (longissimus dorsi and gluteus) of two Dezhou donkey strains (SanFen and WuTou). The differentially expressed lncRNAs were screened between different strains and tissues. Then candidate genes for conjoint analysis were screened based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the accuracy of the sequencing data was verified by real-time quantitative polymerase chain reaction (RT-qPCR). Herein, we identified 3869 novel lncRNAs and 73 differentially expressed lncRNAs. Through the comparison between groups, the specific expression of lncRNAs were found in different strains and muscle tissues. Importantly, we constructed the lncRNA-miRNA-mRNA interaction network and found three important candidate lncRNAs (MSTRG.9787.1, MSTRG.3144.1, and MSTRG.9886.1) and four candidate genes (ACTN1, CDON, FMOD, and BMPR1B). Analysis of the KEGG pathway indicates that the TGF-β signaling pathway plays a pivotal role in the growth and development of skeletal muscle in Dezhou donkey strains.

Characterization of complete mitochondrial genome of Dezhou donkey (Equus asinus) and evolutionary analysis

Mitochondrial DNA (mtDNA) has been widely used in species identification and genetic diversification. Comparisons among mtDNAs of closely related species are valuable for phylogenetic studies. However, only the partial mtDNA Cytb gene and the D-loop sequences were used to analysis the phylogenetic relationship between donkey breeds due to lack of complete mitochondrial genome. Dezhou donkey, as a bigger somatotype ass, is one of Chinese domestic donkey breeds, and used by many places as breeding stock. To further investigate the phylogenetic relationship of Dezhou donkey with other breeds, the complete mtDNA was firstly sequenced and de novo assembled using next generation sequence data from total genomic DNA. The genome was 16,813 bp in length (NCBI submission number: KT182635) and contained 13 protein coding genes, 2 ribosomal RNA genes, 25 transfer RNA genes, and 1 control region. Based on the novel complete mtDNA sequence, the sequences of 13 protein coding genes and 2 ribosomal RNA genes were amplifying in other 2 Dezhou donkeys and 3 Yunnan donkeys, respectively. The pattern of genetic variation in horse, wild ass and domestic donkeys among these 15 genes indicated the sequence polymorphisms. The more accurate phylogenetic relationships of donkey species (Dezhou donkey, Yunnan donkey and previously published donkeys) were first obtained using the combined sequences of 12S rRNA+16S rRNA+13 protein-coding genes. Molecular-based phylogeny supported the hypothesis that Chinese domestic donkey breeds may have originated from Somali wild ass, not from Asian wild ass by analyzing mitochondrial genomes.