Study of the genetic origin of the Mexican creole donkey (Equus asinus) by means of the analysis of the D-loop region of mitochondrial DNA

Next-Generation Sequencing of the Complete Huaibei Grey Donkey Mitogenome and Mitogenomic Phylogeny of the Equidae Family

The Huaibei grey donkey (HGD) is an endangered species and a vital native breed in Anhui Province, China. However, its complete mitogenome, phylogeny, and maternal origin remain unclear. The objectives of this study were to detect the genetic diversity of the HGD and investigate its phylogenetic relationship with other breeds to inform conservation management. The complete mitogenome of the HGD was sequenced through next-generation sequencing, and the most variable region in the mitochondrial DNA displacement-loop (D-loop) was amplified via a polymerase chain reaction (PCR). Next, we used the median-joining network (MJN) to calculate the genetic relationships among populations and the neighbor-jointing method to build a phylogenetic tree and speculate as to its origin. The results showed that the mitogenome contains 22 tRNAs, 2 rRNAs, 13 PCGs, and 1 D-loop region. Analyzing the D-loop region of the HGDs, we identified 23 polymorphic sites and 11 haplotypes. The haplotype and nucleotide diversity were 0.87000 (Hd) and 0.02115 (Pi), respectively. The MJN analysis indicated that the HGD potentially has two maternal lineages, and phylogenetic analysis indicated that the Somali lineage could be the most probable domestication center for this breed. Therefore, our mitogenome analysis highlights the high genetic diversity of the HGD, which may have originated from the Somali wild ass, as opposed to the Asian wild ass. This study will provide a useful resource for HGD conservation and breeding.

Genetic origin of donkeys in Brazil

INTRODUCTION

The genetic groups of native donkeys in Brazil are characterized by adaptation to the local environment. However, the donkey population in the country is declining, mainly because of agricultural mechanization and transportation that has led to the abandonment and the consequent indiscriminate slaughter of these animals. There are three local genetic groups of distinct geographic and temporal formation. However, analyses of their origin, phylogenetic relationship, and population structure are scarce. Within this context, molecular markers such as the mitochondrial control region (D-loop) are useful for these analyses.

MATERIAL AND METHODS

This study aimed to evaluate the variation and origin of maternal lineages of groups of naturalized donkeys in Brazil (Brazilian, Nordestino, and Pêga).

RESULTS AND DISCUSSION

We detected five mitochondrial haplotypes with 19 polymorphic sites, two of them exclusively found in the Nordestino donkey; this group is in fact more distant from the others. Phylogenetic analysis indicates maternal contributions of two clades (Nubian and Somali) to the formation of the genetic groups of donkeys, a fact that explains the high diversity, structure, and distances of the groups, reported here for the first time. CONCLUSION : This analysis contributes production and conservation of native donkey breeds. It also gives clues about the formation of the Iberian breeds from which Brazilian donkeys originated.

Genetic diversity and population genetic structure in native Ethiopian donkeys (Equus asinus) inferred from equine microsatellite markers

We investigated the genetic diversity and population genetic structure of six morphologically distinct Ethiopian donkey populations using 12 equine microsatellite markers. The donkey populations were Abyssinian (AB), Afar (AF), Hararghe (HA), Ogaden (OG), Omo (OM) and Sinnar (SI). Blood samples were collected from 180 genetically unrelated donkeys (30 individuals per population). Population genetic diversity estimates showed that total number and mean number of observed alleles, average observed and expected heterozygosity were 94, 5.208 ± 0.0229, 0.555 ± 0.023 and 0.588 ± 0.022, respectively. Highly significant deficiency in heterozygote was detected within the overall samples (F_IS = 0.055 ± 0.021; P < 0.001). Though highly significant (P < 0.001), heterozygote deficiency within populations relative to total population was moderate (F_ST = 0.046 ± 0.016), suggesting a higher diversity within the populations (95.4%) than between populations. Various genetic distance estimation methods produced a similar topology of un-rooted dendrograms that grouped the overall Ethiopian donkeys into lowland (Ogaden, Omo and Sinnar) and highland (Abyssinian, Afar and Hararghe) genetic lineages. Likewise, Bayesian clustering analysis produced a similar pattern of clustering that was highly concordant with traditional donkey classification systems in Ethiopia.

Phylogenetic Relationships of Turkish Indigenous Donkey Populations Determined by Mitochondrial DNA D-loop Region

Simple Summary

This paper represents the first fundamental report of mtDNA diversity in Turkish indigenous donkey breeds and presents findings for the origin and genetic characterization of donkey populations dispersed in seven geographical regions in Turkey, and thus reveals insights into their genetic history. The median-joining network and phylogenetic tree exhibit two different maternal lineages of the 16 Turkish indigenous donkey populations.

Abstract

In this study, to analyze the mtDNA D-loop region and the origin of the maternal lineages of 16 different donkey populations, and to assess the domestication of Turkish indigenous donkeys in seven geographical regions, we investigated the DNA sequences of the D-loop region of 315 indigenous donkeys from Turkey. A total of 54 haplotypes, resulting from 35 polymorphic regions (27 parsimoniously informative and 6 singleton sites), were defined. Twenty-eight of these haplotypes are unique (51.85%), and 26 are shared among different Turkish indigenous donkey populations. The most frequent haplotype was Hap 1 (45.71%), followed by two haplotypes (Hap 4, 15.55% and Hap 7, 5.39%). The breed genetic diversity, evaluated by the haplotype diversity (H_D) and nucleotide diversity (π_D), for the Turkish donkey populations ranged from 0.533 ± 0.180 (Tekirdağ–Malkara, MAL) to 0.933 ± 0.122 (Aydin, AYD), and from 0.01196 ± 0.0026 (Antalya, ANT) to 0.02101 ± 0.0041 (Aydin, AYD), respectively. We observed moderate-to-high levels of haplotype diversity and moderate nucleotide diversity, indicating plentiful genetic diversity in all of the Turkish indigenous donkey populations. Phylogenetic analysis (NJT) and median-joining network analysis established that all haplotypes were distinctly grouped into two major haplogroups. The results of AMOVA analyses, based on geographic structuring of Turkish native donkey populations, highlighted that the majority of the observed variance is due to differences among samples within populations. The observed differences between groups were found to be statistically significant. Comparison among Turkish indigenous donkey mtDNA D-loop regions and haplotypes, and different countries’ donkey breeds and wild asses, identified two clades and which is named Somali (Clade IV) and Nubian (Clade V) lineages. The results can be used to understand the origin of Turkish donkey populations clearly, and to resolve the phylogenetic relationship among all of the different regions.

Genetic diversity and maternal origin of Northeast African and South American donkey populations

To investigate the mtDNA variation and origin of maternal lineages in South American donkeys and to reassess the domestication of donkeys in northeast Africa, we analyzed sequences (489 bp of the D-loop) from 323 domestic donkeys sampled from Peru, Brazil, Ethiopia and Egypt. Altogether, the 323 sequences displayed 53 different haplotypes (45 in Ethiopia, 14 in Egypt, eight in Peru and six in Brazil). Among the four populations, Egyptian donkeys possessed the highest haplotype diversity (0.910 ± 0.032), followed by Brazilian donkeys (0.879 ± 0.060). The Clade I haplotypes dominated in Peruvian donkeys (65%), whereas Clade II haplotypes dominated in Brazilian donkeys (67%). Estimates of F_ST values showed a high genetic differentiation between Peruvian and Brazilian donkey populations (F_ST  = 0.4066), which could be explained by the complex introduction history of South American donkeys. Phylogeographic analysis indicates that northeast Africa could be the most probable domestication center for Clade I donkeys. Analysis of molecular variance confirmed a weak genetic structure in domestic donkey populations among four continents (Europe, Asia, Africa and South America).

Mitochondrial DNA genetic diversity in six Italian donkey breeds (Equus asinus)

Donkeys have played an important role in agricultural land practices and in human historical periods of recent past and, still today, are used as a working power in several world areas. The objective of this study was to identify genetic variability in six Italian donkey breeds using mtDNA D-loop. Fifteen haplotypes, grouped in three haplogroups, were identified. The genetic indices were informative and showed a high population genetic variability. The results of AMOVA analyses based on geographic structuring of Italian populations highlighted that the majority of the observed variance is due to differences among samples within breeds. Comparison among Italian haplotypes and mtDNA D-loop sequences belonging to European domestic and Ethiopian donkeys and wild asses, clearly define two clades referred to Nubian lineage. The results can be useful to complement safeguard planes for donkey breeds that are considered to extinction endangered.

Genetic relationships among American donkey populations: insights into the process of colonization

This study presents the first insights into the genetic diversity and structure of the American donkey metapopulation. The primary objectives were to detect the main structural features underlying variability among American donkey populations, identify boundaries between differentiated gene pools, and draw the main colonization pathways since the introduction of donkeys into America in the 15th century. A panel of 14 microsatellite markers was applied for genotyping 350 American donkeys from 13 countries. The genetic structure of this metapopulation was analysed using descriptive statistics and Bayesian model-based methods. These populations were then compared to a database containing information on 476 individuals from 11 European breeds to identify the most likely ancestral donor populations. Results showed the presence of two distinct genetic pools, with confluence of the two in Colombia. The southern pool showed a unique genetic signature subsequent to an older founder event, but lacked any significant influence of modern gene flow from Europe. The northern pool, conversely, may have retained more ancestral polymorphisms and/or have experienced modern gene flow from Spanish breeds. The Andalusian and, to a lesser extent, the Catalan breeds have left a more pronounced footprint in some of the American donkey populations analysed.

Next Generation Semiconductor Based Sequencing of the Donkey (Equus asinus) Genome Provided Comparative Sequence Data against the Horse Genome and a Few Millions of Single Nucleotide Polymorphisms

Few studies investigated the donkey (Equus asinus) at the whole genome level so far. Here, we sequenced the genome of two male donkeys using a next generation semiconductor based sequencing platform (the Ion Proton sequencer) and compared obtained sequence information with the available donkey draft genome (and its Illumina reads from which it was originated) and with the EquCab2.0 assembly of the horse genome. Moreover, the Ion Torrent Personal Genome Analyzer was used to sequence reduced representation libraries (RRL) obtained from a DNA pool including donkeys of different breeds (Grigio Siciliano, Ragusano and Martina Franca). The number of next generation sequencing reads aligned with the EquCab2.0 horse genome was larger than those aligned with the draft donkey genome. This was due to the larger N50 for contigs and scaffolds of the horse genome. Nucleotide divergence between E. caballus and E. asinus was estimated to be ~ 0.52-0.57%. Regions with low nucleotide divergence were identified in several autosomal chromosomes and in the whole chromosome X. These regions might be evolutionally important in equids. Comparing Y-chromosome regions we identified variants that could be useful to track donkey paternal lineages. Moreover, about 4.8 million of single nucleotide polymorphisms (SNPs) in the donkey genome were identified and annotated combining sequencing data from Ion Proton (whole genome sequencing) and Ion Torrent (RRL) runs with Illumina reads. A higher density of SNPs was present in regions homologous to horse chromosome 12, in which several studies reported a high frequency of copy number variants. The SNPs we identified constitute a first resource useful to describe variability at the population genomic level in E. asinus and to establish monitoring systems for the conservation of donkey genetic resources.

Ancient DNA provides new insight into the maternal lineages and domestication of Chinese donkeys

Background

The donkey (Equus asinus) is an important domestic animal that provides a reliable source of protein and method of transportation for many human populations. However, the process of domestication and the dispersal routes of the Chinese donkey are still unclear, as donkey remains are sparse in the archaeological record and often confused with horse remains. To explore the maternal origins and dispersal route of Chinese donkeys, both mitochondrial DNA D-loop and cytochrome b gene fragments of 21 suspected donkey remains from four archaeological sites in China were amplified and sequenced.

Results

Molecular methods of species identification show that 17 specimens were donkeys and three samples had the maternal genetic signature of horses. One sample that dates to about 20,000 years before present failed to amplify. In this study, the phylogenetic analysis reveals that ancient Chinese donkeys have high mitochondrial DNA diversity and two distinct mitochondrial maternal lineages, known as the Somali and Nubian lineages. These results indicate that the maternal origin of Chinese domestic donkeys was probably related to the African wild ass, which includes the Nubian wild ass (Equus africanus africanus) and the Somali wild ass (Equus africanus somaliensis). Combined with historical records, the results of this study implied that domestic donkeys spread into west and north China before the emergence of the Han dynasty. The number of Chinese domestic donkeys had increased primarily to meet demand for the expansion of trade, and they were likely used as commodities or for shipping goods along the Silk Road during the Tang Dynasty, when the Silk Road reached its golden age.

Conclusions

This study is the first to provide valuable ancient animal DNA evidence for early trade between African and Asian populations. The ancient DNA analysis of Chinese donkeys also sheds light on the dynamic process of the maternal origin, domestication, and dispersal route of ancient Chinese donkeys.