Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys.

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

A mutation in POLR2A gene associated with body size traits in Dezhou donkeys revealed with GWAS.

The Dezhou donkey is a famous local donkey breed in China. The aim of the present study was to identify the genes associated with the body size traits of the Dezhou donkey and facilitate the breeding activities of the donkeys. A total of 349 donkeys from 2 generations (113 individuals in F0 and 236 in F1) were analyzed with restriction-site-associated DNA sequencing. A genome-wide association study revealed that the region between 13.7 and 15.6 Mb of chromosome 13 is significantly associated with body sizes. Candidate genes related to body size development, including POLR2A, CHRNB1, FGF11, and ZBTB4, were identified. The results of GO and KEGG analysis indicated that the genes involved in many GO terms were related to metabolic processes and developmental processes. Additionally, a T>C mutation (Chr13:14312485) was found at intron 10 of the POLR2A gene. The association analysis showed significant differences among genotypes for the size traits. The body size of the individuals with the TT genotype was significantly higher than that with the CC genotype. The results showed that the polymorphism of POLR2A has the potential to be used as a marker in the breeding programs of the Dezhou donkeys.

Body size is a crucial economic trait in donkeys, as it is closely related to meat and skin production. The aim of this study was to identify the genes and loci associated with body size traits, using the Dezhou donkey as an experimental population. The study findings make contributions to a better understanding on the molecular genetic mechanism of body size traits. The significant loci screened out in the present study may facilitate gene-assisted selection breeding and accelerate genetic selection, which is of great significance to the breeding of donkeys and the development of the donkey industry.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Donkeys assume relevant economic importance in several countries worldwide. However, the genetic structure of these populations is less investigated compared to other species. The aim of this study was to investigate the genetic background of nine different Italian donkey populations. A total of 220 animals were genotyped with about 27K markers extracted by the double-digest restriction site associated DNA sequencing. The consistency of the results across different approaches agreed with the demographic history, the origin, and previous results on the nine donkey populations, suggesting that our conclusions are robust. Moreover, the results of the present study highlighted low to moderate levels of admixture and, for some breeds, confirmed the success of proper management conservation plans.

Analysis of the Genetic Relationship and Inbreeding Coefficient of the Hetian Qing Donkey through a Simplified Genome Sequencing Technology.

The Hetian Qing donkey is an excellent local donkey breed in Xinjiang. It is of great significance to accelerate breeding and the speed of breeding and rejuvenation, as well as to understand the genetic basis of the strategies and population. This study collected a total of 4 male donkeys and 28 female donkeys. It then obtained genotype data through Simplified Genomic Sequencing (GBS) technology for data analysis. The results detected a total of 55,399 SNP loci, and the genotype detection rate of individuals was ≥90%. A total of 45,557 SNP loci were identified through quality control, of which 95.5% were polymorphic. The average minimum allele frequency was 0.250. The average observed heterozygosity was 0.347. The average expected heterozygosity was 0.340. The average IBS (state homologous) genetic distance was 0.268. ROH: 49 (homozygous fragments), with 73.47% of the length between 1 and 5 Mb. The average per-strip ROH length was 1.75 Mb. The mean inbreeding coefficient was 0.003. The 32 Hetian green donkeys could be divided into six families. The number of individuals in each family is significant. To sum up, the Hetian Qing donkey population has low heterozygosity, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the subsequent process of seed protection, seed selection should be conducted according to the divided pedigree to ensure the long-term protection of the genetic resources of Hetian green donkeys.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Double-Tube Multiplex TaqMan Real-Time PCR for the Detection of Eight Animal-Derived Dairy Ingredients.

Milk and dairy products represent important sources of nutrition in our daily lives. The identification of species within dairy products holds importance for monitoring food adulteration and ensuring traceability. This study presented a method that integrated double-tube and duplex real-time polymerase chain reaction (PCR) with multiplex TaqMan probes to enable the high-throughput detection of animal-derived ingredients in milk and dairy products. The detection system utilized one pair of universal primers, two pairs of specific primers, and eight animal-derived specific probes for cow, buffalo, goat, sheep, camel, yak, horse, and donkey. These components were optimized within a double-tube and four-probe PCR multiplex system. The developed double-tube detection system could simultaneously identify the above eight targets with a detection limit of 10-0.1 pg/µL. Validation using simulated adulterated milk samples demonstrated a detection limit of 0.1%. The primary advantage of this method lies in the simplification of the multiplex quantitative real-time PCR (qPCR) system through the use of universal primers. This method provides an efficient approach for detecting ingredients in dairy products, providing powerful technical support for market supervision.

Comparative molecular epidemiology, subtype distribution, and zoonotic potential of Blastocystis sp. in Equus animals (horses, donkeys, and mules) in northwestern Iran.

A total of 500 fecal samples were collected from Equus animals in six different cities (Ardabil, Namin, Nir, Meshginshahr, Germi, and Khalkhal) of Ardabil Province, northwestern Iran, with 200 samples from horses, 200 from donkeys, and 100 from mules. Of the horse samples, 100 were from racing horses under special monitoring and care, while the remaining 100 were from non-racing horses, including those used for herding or in rural areas. All fecal samples were examined for the presence of Blastocystis sp. using PCR amplification of the SSU rRNA gene's barcode region after DNA extraction. The molecular prevalence of Blastocystis infection in Equus animals was 7.6% (38/500). Blastocystis was more common in horses [11.5% (23/200)] than in donkeys [5.5% (11/200)] and mules [4% (4/100)] (P > 0.05). Compared to racing horses [3% (3/100)], non-racing/rural horses [20% (20/100)] exhibited a substantially higher prevalence of Blastocystis (P < 0.05). The prevalence of Blastocystis in diarrheal samples and younger animals was remarkably higher than in formed samples and older animals, respectively (P < 0.05). No significant difference in Blastocystis infection prevalence was found between the genders of examined animals (P > 0.05). In Equus animals, 38 Blastocystis isolates included eight STs: ST10 [31.6% (12/38)], ST1 [21.1% (8/38)], ST2 [15.8% (6/38)], ST3 [10.5% (4/38)], ST4 [7.9% (3/38)], ST7 [5.2% (2/38)], ST14 [5.2% (2/38)], and ST6 [2.6% (1/38)]. These results suggest that Equus animals act as a proper reservoir for numerous Blastocystis STs, consequently playing a crucial part in the spread of this protozoan infection to humans, animals, and water reservoirs.

Genome-wide analyses based on a novel donkey 40K liquid chip reveal the gene responsible for coat color diversity in Chinese Dezhou donkey.

Dezhou donkey is one of the representative local breeds in China, which is mainly divided into two strains: Sanfen and Wutou. There are obvious differences in coat color between the two strains. The former shows light points around the eyes, around the muzzle and under the belly, while the latter is completely solid black. In this study, genome-wide association analysis was performed for the differences in coat color traits between the Sanfen (n = 97) and Wutou (n = 108) strains using a novel donkey 40K liquid chip developed based on GenoBaits technology, to identify genomic regions and causal genes that could explain this variation. We also used FST and The cross-population composite likelihood ratio test (XPCLR) analyses to explore selected regions related to coat color differences. We identified one significant region on chromosome 15, with the most significant SNP located within the agouti signaling protein (ASIP) gene. At the same time, both FST and XPCLR methods detected the same selected region on chromosome 15, and ASIP was the gene with the strongest signal. ASIP and melanocortin 1 receptor (MC1R) control the ratio of eumelanin to pheomelanin through their protein activity. They are deeply involved in the process of melanosome organation and melanogenesis, thus affecting mammals' coat color variation. We used a range of genome-wide approach to identify the genetic basis of coat color variation in Dezhou donkeys. The results provide a supplement to the color variation study in Chinese donkeys at the genome-wide level, and preliminarily verified the reliability of the Molbreeding Donkey No. 1 40K liquid chip.

Exploring testicular miRNA profiles during developmental stages in Dezhou donkeys: A preliminary insight.

Testicular development and spermatogenesis are complex phenomena controlled by various genetic factors, including miRNA-based post-transcriptional gene expression regulation. Exploring the miRNA expression patterns during testicular development in Dezhou donkeys would enhance our understanding of equine fertility and spermatogenesis. In this investigation, we examined the testicular miRNA profiles at various stages of development. The experimental animals were divided into three groups based on their developmental stages: 2 months old (juvenile: n = 3), 12 months old (adolescent; n = 3) and 24 months old (adult; n = 3) donkeys. Total RNA was extracted from dissected testicles for miRNA sequencing and analysis. In total, 586 miRNAs, including 451 known miRNAs and 135 novel miRNAs, were identified. Among identified miRNAs, 315 displayed age-dependent expression differences. The levels of miRNA expression in the juvenile group were significantly higher than in the adolescent or adult groups. The MiR-483 exhibited the maximum fold change between juvenile and adolescent groups. Several screened genes, including SLC45A4 and TFCP2L1, have been linked to male reproductive pathways in donkeys. In addition, miR-744 was predicted to regulate SPIN2B, a gene implicated in spermatocyte cell cycle progression and genomic integrity of spermatozoa. These results contribute to our comprehension of microRNA regulation during testicular development and spermatogenesis in Dezhou donkeys. The identified microRNAs and their target genes have the potential to serve as biomarkers for evaluating the reproductive capacity of stud donkeys.

Transcriptome analysis reveals immune function-related mRNA expression in donkey mammary glands during four developmental stages.

The low susceptibility to mastitis of female donkey (jenny) mammary glands and the strong immune properties of donkey milk are acknowledged, but little is known about the genes involved in mammary gland immunity in jennies. Herein, we used RNA-sequencing and bioinformatics analyses to explore jenny mammary gland transcriptomes and detect potential functional differentially expressed (DE) mRNAs related to immunity during four specific developmental stages: foetal (F), pubertal (P), adult parous nonlactation (N) and lactation (L). A total of 2497, 583 and 1820 DE mRNAs were identified in jenny mammary glands at F vs. P, P vs. N, and N vs. L, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) analyses revealed numerous GO terms related to immune function, especially between F and P. Seven significantly enriched profiles were identified, among which 497 and 1261 DE mRNAs were upregulated in profiles 19 and 17. Eleven mRNAs were enriched in over 10 KEGG pathways. ß-2-microglobulin (B2M), immunoglobulin heavy constant mu (IGHM), toll like receptor 2 (TLR2), toll like receptor 4 (TLR4) and myeloid differentiation factor 88 (MYD88) were mainly involved in phosphoinositide 3-kinase (PI3K)-Akt signalling, phagosome and nuclear factor kappa-B (NF-kappa B) signalling pathways. The findings provide insight into the molecular features underpinning the low prevalence of intramammary infections (i.e., mastitis) in donkeys.

Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

Genomic and comparative analysis of the T cell receptor gamma locus in two Equus species.

The genus Equus is the only extant genus of the Equidae family, which belongs to Perissodactyla, an order of mammals characterized by an odd number of toes (odd-toes ungulates). Taking advantage of the latest release of the genome assembly, we studied, for the first time in two organisms belonging to the Equus genus, the horse (Equus caballus) and the donkey (Equus asinus), the T cell receptor gamma (TRG) locus encoding the gamma chain of the γδ T cell receptor. Forty-five Variable (TRGV) genes belonging to the seven IMGT-NC validated mammalian TRGV subgroups, 25 Joining (TRGJ) and 17 Constant (TRGC) genes organized in 17 V-J-(J)-C cassettes, in tandem on about 1100 Kb, characterize the horse TRG locus, making the horse TRG locus the one with the greatest extension and with a significantly higher number of genes than the orthologous loci of the other mammalian species. A clonotype analysis of an RNA-seq transcriptomic dataset derived from spleen of an adult healthy horse, using the complete set of the horse TRGJ germline gene sequences as a probe, revealed that, in addition to the most prominent V-J rearrangements within each cassette, there is a relevant proportion of trans-cassette V-J recombination, whereby the same TRGV genes can recombine with different TRGJ genes spliced to the corresponding TRGC genes. This recombinant event strongly contributes to the diversity of the γ chain repertoire. In the donkey TRG locus, 34 TRGV, 21 TRGJ and 14 TRGC genes distributed in 14 V-J-(J)-C cassettes were found in a region of approximately 860 kb. Although the donkey's TRG is smaller than that of the horse, in Equus genus, this is still the second largest locus so far found in any mammalian species. Finally, the comparative analysis highlighted differences in size and gene content between the horse and donkey TRG loci, despite belonging to the same genus, indicating a good level of diversification within Equus. These data is in agreement with the evolutionary idea of the existence of a Equus recent common ancestor in rapid evolution, for which a mutation rate between horses and donkeys is more comparable to that between species belonging to different genera rather than to species of the same genus.

Comparative Genomics Identifies the Evolutionarily Conserved Gene TPM3 as a Target of eca-miR-1 Involved in the Skeletal Muscle Development of Donkeys.

Species within the genus Equus are valued for their draft ability. Skeletal muscle forms the foundation of the draft ability of Equus species; however, skeletal muscle development-related conserved genes and their target miRNAs are rarely reported for Equus. In this study, a comparative genomics analysis was performed among five species (horse, donkey, zebra, cattle, and goat), and the results showed that a total of 15,262 (47.43%) genes formed the core gene set of the five species. Only nine chromosomes (Chr01, Chr02, Chr03, Chr06, Chr10, Chr18, Chr22, Chr27, Chr29, and Chr30) exhibited a good collinearity relationship among Equus species. The micro-synteny analysis results showed that TPM3 was evolutionarily conserved in chromosome 1 in Equus. Furthermore, donkeys were used as the model species for Equus to investigate the genetic role of TPM3 in muscle development. Interestingly, the results of comparative transcriptomics showed that the TPM3 gene was differentially expressed in donkey skeletal muscle S1 (2 months old) and S2 (24 months old), as verified via RT-PCR. Dual-luciferase test analysis showed that the TPM3 gene was targeted by differentially expressed miRNA (eca-miR-1). Furthermore, a total of 17 TPM3 gene family members were identified in the whole genome of donkey, and a heatmap analysis showed that EaTPM3-5 was a key member of the TPM3 gene family, which is involved in skeletal muscle development. In conclusion, the TPM3 gene was conserved in Equus, and EaTPM3-5 was targeted by eca-miR-1, which is involved in skeletal muscle development in donkeys.

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.

Mitogenomic phylogeny reveals the predominance of the Nubian lineage of African wild ass in Indian donkeys.

To contribute to the knowledge of maternal genetic diversity in domestic donkeys, this study investigated the mitochondrial DNA variations and analyzed the genetic structure in Indian donkeys based on 31 mitogenome sequences representing four breeds/populations (Agra, Halari, Kachchhi and Spiti). A total of 27 haplotypes with a haplotype diversity value of 0.989 were evident in the donkey genetic resources of India. The genetic differentiation between the investigated populations was evaluated using population pairwise FST values, which showed maximum differentiation between Kachchhi and Halari donkeys. The Neighbor-Joining (NJ) tree based on the whole mitogenome sequence and the Median-Joining (MJ) network for partial D-loop fragment showed clear demarcation of Indian donkeys into Nubian and Somali clades, substantiating African maternal origin of Indian domestic donkeys. The topology of the MJ network excluded the Asian wild asses as the possible progenitors of Indian donkeys. Halari and Agra donkeys showed conformity exclusively to the Nubian lineage of the African wild asses. However, representation of both the Nubian and Somali lineages was observed in Kachchhi and Spiti donkeys. Comprehensive analysis carried out by retrieving D-loop sequences from different countries representing Asia, Africa, Europe and South America revealed existence of shared haplotypes across geographically isolated regions of the globe. This observation is indicative of utility of donkeys as pack animals across inter-continental trading routes during development of human civilizations. Our results represent a valuable contribution to maternal genetic diversity of Indian donkeys and provide insights into the worldwide spread of the species following initial domestication in Africa.

Ancient segmentally duplicated LCORL retrocopies in equids.

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.

A Genome-Wide Association Study of the Chest Circumference Trait in Xinjiang Donkeys Based on Whole-Genome Sequencing Technology.

Animal genotyping by means of genome-wide association studies is important for connecting phenotypes of interest with their underlying genetics in livestock. However, the use of whole genome sequencing to investigate chest circumference (CC) in donkeys has rarely been reported. We aimed to use the genome-wide association study approach to detect significant single nucleotide polymorphisms (SNPs) and key genes associated with chest circumference traits in Xinjiang donkeys. We assessed 112 Xinjiang donkeys in this study. The chest circumference of each was measured 2 h before milking. We re-sequenced blood samples from the Xinjiang donkeys, and genome-wide association study analyses were performed using a mixed model with the PLINK, GEMMA, and REGENIE programs. We tested 38 donkeys for candidate SNPs for genome-wide association study using three software programs. Additionally, 18 SNP markers reached genome-wide significance (p < 1.61 × 10-9). On the basis of these, 41 genes were identified. Previously proposed candidate genes for CC traits were supported by this study, including NFATC2 (Nuclear Factor of Activated T Cells 2), PROP1 (PROP Paired-Like Homeobox 1), UBB (Ubiquitin B), and HAND2 (Heart and Neural Crest Derivatives Expressed 2). These promising candidates provide a valuable resource for validating potential meat production genes and will facilitate the development of high-yielding Xinjiang donkey breeds through marker-assisted selection or gene editing.

Effects of ex situ conservation on diversity and function of the gut microbiota of the Tibetan wild ass (Equus kiang).

Ex situ conservation is the main method for the protection of endangered wildlife. To explore the effect of ex situ conservation on the gut microbiota of the kiang (Equus kiang), metagenomic sequencing combined with bioinformatics analysis was used to investigate the composition and function of the gut microbiota of the kiang. The results showed that ex situ conservation not only protected wildlife, but also affected the composition and function of gut microbiota, as well as the health of animals. In the zoo, the ratio of the relative abundance of Firmicutes to that of Bacteroidetes (F/B) is higher, clusters of potentially pathogenic bacteria (such as Catonella, Catonella, and Mycoplasma) are more numerous, the abundance of resistance genes is higher, and the abundance of metabolic functions is increased. The dynamic changes of the gut microbiota also played an important role in the nutritional absorption, energy metabolism, and environmental adaptation of the kiang. Improving the rearing environment and increasing food diversity play important roles for increasing the diversity of gut microbiota, reducing the spread of potentially pathogenic bacteria, and reducing diseases. In the wild, especially in winter and in food-deficient areas, food supplementation can enhance the gut microbial homeostasis of wild animals and reduce the impact of crises. In depth studies of the gut microbial function of wildlife have important implications for improving ex situ conservation.

Identification and characterization of single nucleotide polymorphisms in DMRT3 gene in Indian horse (Equus caballus) and donkey (Equus asinus) populations.

Equines' ability in racing and riding as well as gaitedness have influenced the human civilization. Aim of this study was to identify and characterize the novel polymorphisms or SNPs in DMRT3 gene in Indian horse and donkey breeds. In this study, the DMRT3 gene was sequenced and characterized in 72 Indian horses' and 33 Indian donkeys' samples. One SNP (A > C) at 878 was found in studied horses while identical SNPs (A > C) at two different nucleotide positions i.e., 878 and 942 in DMRT3 gene (chromosome 23) were observed in studied Indian donkey breeds. Horses and donkeys both have a non-synonymous mutation (A > C) at nucleotide 878 (codon 61) that converts a Stop codon (TAG > TCG) to coding codon Serine, whereas donkeys have a synonymous mutation at nucleotide 942 (codon 82) that converts Serine (TCA > TCC) into Serine. A phylogenetic tree indicated that the DMRT3 gene was equally distributed among the equine breeds. Most of the donkey breeds have been shown high levels of genetic diversity while horse breeds and Halari donkey showed the least genetic diversity. Mutation in DMRT3 has a major impact on gaitedness in horses and is presented at a high frequency in gaited breeds and in horses breed for harness racing.