Genetic diversity and variability in endangered Pantesco and two other Sicilian donkey breeds assessed by microsatellite markers

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.

An investigation of genetic diversity in three Dezhou donkey original breeding farms

Dezhou donkey is one of the excellent large donkey breeds in China. In our study, eight microsatellite markers were used to genotype from each of three populations of Chinese Dezhou donkeys: 67 individuals in Liaocheng (pop1), 103 individuals in Binzhou 1 (pop2), and 102 individuals in Binzhou 2 (pop3), in order to study the genetic diversity of these populations. A total of 213 alleles were detected, and the PIC results showed that eight loci were all highly polymorphic. The means of Ho and He of pop2 were the highest, which were 0.646 and 0.717, respectively. The PCoA analysis demonstrated that the samples from three conservation farms were mixed together. The phylogenetic tree showed that pop 2 and pop 3 were closely related. The phylogenetic tree results clustered that 272 donkeys were divided into six groups. AMOVA analysis showed that the genetic variation was mainly concentrated within the population and the genetic differentiation among the populations was low. Fst values between populations also indicated that genetic differentiation between populations was too small to be considered. This indicated a low probability of inbreeding within the population. And this also showed that the conservation and breeding of Dezhou donkeys in recent years have achieved excellent results. The investigation of genetic diversity in three Dezhou donkey original breeding farms can provide reference data for the selection and breeding of excellent breeds of Dezhou donkey.

Origin, Evolution, and Research Development of Donkeys

Lack of archaeological and whole-genome diversity data has restricted current knowledge of the evolutionary history of donkeys. With the advancement of science and technology, the discovery of archaeological evidence, the development of molecular genetics, and the improvement of whole-genome sequencing technology, the in-depth understanding of the origin and domestication of donkeys has been enhanced, however. Given the lack of systematic research, the present study carefully screened and collected multiple academic papers and books, journals, and literature on donkeys over the past 15 years. The origin and domestication of donkeys are reviewed in this paper from the aspects of basic information, cultural origin, bioarcheology, mitochondrial and chromosomal microsatellite sequences, and whole-genome sequence comparison. It also highlights and reviews genome assembly technology, by assembling the genome of an individual organism and comparing it with related sample genomes, which can be used to produce more accurate results through big data statistics, analysis, and computational correlation models. Background: The donkey industry in the world and especially in China is developing rapidly, and donkey farming is transforming gradually from the family farming model to large-scale, intensive, and integrated industrial operations, which could ensure the stability of product quality and quantity. However, theoretical research on donkey breeding and its technical development lags far behind that of other livestock, thereby limiting its industrial development. This review provides holistic information for the donkey industry and researchers, that could promote theoretical research, genomic selection (GS), and reproductive management of the donkey population.

Detection and Molecular Characterization of Two Gammaherpesviruses from Pantesco Breed Donkeys during an Outbreak of Mild Respiratory Disease

Equid and asinine gammaherpesviruses (GHVs; genus Percavirus) are members of the Herpesviridae family. Though GHVs have been reported in horse populations, less studies are available on gammaherpesviral infections in donkeys. This study reports the co-infection with two GHVs in Pantesco breed donkeys, an endangered Italian donkey breed. Samples (n = 124) were collected on a breeding farm in Southern Italy from 40 donkeys, some of which were healthy or presented erosive tongue lesions and/or mild respiratory signs. Samples were analysed by using a set of nested PCRs targeting the DNA polymerase, glycoprotein B, and DNA-packaging protein genes, and sequence and phylogenetic analyses were performed. Twenty-nine donkeys (72.5%) tested positive, and the presence of Equid gammaherpesvirus 7 and asinine herpesvirus 5 was evidenced. In 11 animals, we found evidence for co-infection with viruses from the two species. Virions with herpesvirus-like morphology were observed by electron microscopic examination, and viruses were successfully isolated in RK-13-KY cell monolayers. The histological evaluation of tongue lesions revealed moderate lympho-granulocytic infiltrates and rare eosinophilic inclusions. The detection of GHVs in this endangered asinine breed suggests the need long-life monitoring within conservation programs and reinforces the need for further investigations of GHV’s pathogenetic role in asinine species.

Maternal Phylogenetic Relationships and Genetic Variation among Rare, Phenotypically Similar Donkey Breeds

The mitochondrial DNA (mtDNA) D-loop of endangered and critically endangered breeds has been studied to identify maternal lineages, characterize genetic inheritance, reconstruct phylogenetic relations among breeds, and develop biodiversity conservation and breeding programs. The aim of the study was to determine the variability remaining and the phylogenetic relationship of Martina Franca (MF, with total population of 160 females and 36 males), Ragusano (RG, 344 females and 30 males), Pantesco (PT, 47 females and 15 males), and Catalonian (CT) donkeys by collecting genetic data from maternal lineages. Genetic material was collected from saliva, and a 350 bp fragment of D-loop mtDNA was amplified and sequenced. Sequences were aligned and evaluated using standard bioinformatics software. A total of 56 haplotypes including 33 polymorphic sites were found in 77 samples (27 MF, 22 RG, 8 PT, 19 CT, 1 crossbred). The breed nucleotide diversity value (π) for all the breeds was 0.128 (MF: 0.162, RG: 0.132, PT: 0.025, CT: 0.038). Principal components analysis grouped most of the haplogroups into two different clusters, I (including all haplotypes from PT and CT, together with haplotypes from MF and RG) and II (including haplotypes from MF and RG only). In conclusion, we found that the primeval haplotypes, haplogroup variability, and a large number of maternal lineages were preserved in MF and RG; thus, these breeds play putative pivotal roles in the phyletic relationships of donkey breeds. Maternal inheritance is indispensable genetic information required to evaluate inheritance, variability, and breeding programs.

Genomics and the Evolutionary History of Equids

The equid family contains only one single extant genus, Equus, including seven living species grouped into horses on the one hand and zebras and asses on the other. In contrast, the equine fossil record shows that an extraordinarily richer diversity existed in the past and provides multiple examples of a highly dynamic evolution punctuated by several waves of explosive radiations and extinctions, cross-continental migrations, and local adaptations. In recent years, genomic technologies have provided new analytical solutions that have enhanced our understanding of equine evolution, including the species radiation within Equus; the extinction dynamics of several lineages; and the domestication history of two individual species, the horse and the donkey. Here, we provide an overview of these recent developments and suggest areas for further research.

Genetic Characterization of Native Donkey (Equus asinus) Populations of Turkey Using Microsatellite Markers

Simple Summary

This study was conducted to evaluate the genetic variability of Turkish native donkey (Equus asinus) populations, using polymorphism of 17 microsatellite markers. The results revealed a highly mixed genotype of all the examined donkeys, suggesting that two different group of breeds can be distinguished from each other on the basis of microsatellite markers.

Abstract

This study presents the first insights to the genetic diversity and structure of the Turkish donkey populations. The primary objectives were to detect the main structural features of Turkish donkeys by microsatellite markers. A panel of 17 microsatellite markers was applied for genotyping 314 donkeys from 16 locations of Turkey. One hundred and forty-two alleles were identified and the number of alleles per locus ranged from 4 to 12. The highest number of alleles was observed in AHT05 (12) and the lowest in ASB02 and HTG06 (4), while ASB17 was monomorphic. The mean H_O in the Turkish donkey was estimated to be 0.677, while mean H_E was 0.675. The polymorphic information content (PIC) was calculated for each locus and ranged from 0.36 (locus ASB02) to 0.98 (locus AHT05), which has the highest number of alleles per locus in the present study. The average PIC in our populations was 0.696. The average coefficient of gene differentiation (G_ST) over the 17 loci was 0.020 ± 0.037 (p < 0.01). The G_ST values for single loci ranged from −0.004 for LEX54 to 0.162 for COR082. Nei’s gene diversity index (H_t) for loci ranged from 0.445 (ASB02) to 0.890 (AHT05), with an average of 0.696. A Bayesian clustering method, the Structure software, was used for clustering algorithms of multi-locus genotypes to identify the population structure and the pattern of admixture within the populations. When the number of ancestral populations varied from K = 1 to 20, the largest change in the log of the likelihood function (ΔK) was when K = 2. The results for K = 2 indicate a clear separation between Clade I (KIR, CAT, KAR, MAR, SAN) and Clade II (MAL, MER, TOK, KAS, KUT, KON, ISP, ANT, MUG, AYD and KAH) populations.

Banat donkey, a neglected donkey breed from the central Balkans (Serbia)

The dominant donkey breed in the Balkans is the mid-sized Balkan donkey with a grey to chocolate coat color. Local breeders from Serbia, however, still maintain a few larger individuals of a lighter coat color, named Banat donkey, and speculate that they are descendants of a Spanish donkey heard that had been transferred to the Banat region by the Hapsburg Queen Maria Theresa in the XVIII century for a specific purpose, to work in local vineyards. We have previously found a unique nuclear gene-pool and a prevalence of mitochondrial Clade 2 haplotypes in several such animals. In this study, we: (i) perform a comparative analysis of 18 morphological traits of the Banat donkey (seven individuals), Balkan donkey (53 individuals from two sub-populations of this breed) and the potential hybrids (eight individuals), and demonstrate the morphological distinctiveness of the Banat donkey, highlighting the diagnostic traits for distinguishing the breed: hip height, croup width, body length and chest depth; (ii) re-analyse published nuclear microsatellite data for these groups, and reveal that, although severely depopulated, the genetically distinct Banat donkey is not severely affected by the loss of genetic diversity and inbreeding; (iii) demonstrate that previously published Banat donkey mitochondrial haplotypes, analyzed genealogically together with those reported in ancient and modern individuals from Spain, Italy, Turkey, Cyprus and Africa, are shared with three Spanish breeds and individuals belonging to Amiata and some other Italian breeds. A unique morphological feature present in Banat and Somali wild donkeys, but also in Amiata donkeys, black stripes on legs, suggests that the origin of Clade 2 donkeys may be much more complex than previously thought. Actions to preserve the Banat donkey, a valuable but critically endangered genetic resource (<100 individuals), are urgent.

Genetics, Evolution, and Physiology of Donkeys and Mules

The genus Equus is made up of donkeys, horses, and zebras. Despite significant variation in chromosome number across these species, interspecies breeding results in healthy, although infertile, hybrid offspring. Most notable among these are the horse-donkey hybrids, the mule and hinny. Donkeys presently are used for everything from companion animals to beasts of burden. Although closely related from an evolutionary standpoint, differences in anatomy and physiology preclude the assumption that they can be treated identically to the domestic horse. Veterinarians should be aware of these differences and adjust their practice accordingly.

Genetic diversity and relationships of Chinese donkeys using microsatellite markers

Donkeys are one important livestock in China because of their nourishment and medical values. To investigate the genetic diversity and phylogenetic relationships of Chinese donkey breeds, a panel of 25 fluorescently labeled microsatellite markers was applied to genotype 504 animals from 12 Chinese donkey breeds. A total of 226 alleles were detected, and the expected heterozygosity ranged from 0.6315 (Guanzhong) to 0.6999 (Jiami). The mean value of the polymorphism information content, observed number of alleles, and expected number of alleles for all the tested Chinese donkeys were 0.6600, 6.890, and 3.700, respectively, suggesting that Chinese indigenous donkeys have relatively abundant genetic diversity. Although there were abundant genetic variations found, the genetic differentiation between the Chinese donkey breeds was relatively low, which displayed only 5.99 % of the total genetic variance among different breeds. The principal coordinates analysis clearly splits 12 donkey breeds into two major groups. The first group included Xiji, Xinjiang, Liangzhou, Kulun, and Guanzhong donkey breeds. In the other group, Gunsha, Dezhou, Biyang, Taihang, Jiami, Qingyang, and Qinghai donkeys were clustered together. This grouping pattern was further supported by structure analysis and neighbor-joining tree analysis. Furthermore, genetic relationships between different donkey breeds identified in this study were corresponded to their geographic distribution and breeding history. Our results provide comprehensive and precise baseline information for further research on preservation and utilization of Chinese domestic donkeys.

Simple sequence repeat (SSR) genotypic data reveal high genetic diversity in Rajasthan donkey of India

Spiti is the only breed of donkey in the list of 160 breeds registered so far for the indigenous germplasm. It points towards the fact that Indian donkey diversity is largely unexplored, in spite of its role in support system of poorest of the poor countrymen. Rajasthan state harbours more than eighty thousand donkey. Genetic diversity underpins population resilience and persistence so a need was felt for the systematic and scientific evaluation of existing diversity in this population. The diversity status of Rajasthan donkey population was investigated by using 20 microsatellite markers. Two of the loci, HMS5 and ASB17, depicted less than 4 alleles and were thus excluded from further analyses. Rajasthan donkey was found to host considerable diversity. This conclusion is evidenced by the number of alleles observed across loci (average 8.16±0.89), and by mean observed heterozygosity (0.688±0.06). Signatures of inbreeding could not be detected by the FIS indices in the Rajasthan donkey population. The population did not suffer any bottleneck in the recent past. The study provides the first scientific assessment of the genetic diversity status of Rajasthan donkey population. This information on the genetic analysis can be used in conjunction with the information on physical and phenotypic characteristics and management practices for registration of this population as a breed.

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.

Pedigree and herd characterization of a donkey breed vulnerable to extinction

Most donkey and local horse breeds are vulnerable to extinction as mechanization of agriculture progress throughout the world. The present study analyzed the pedigree and herd records of the donkey Asinina de Miranda breed (RAM), identifying genealogical and human factors that may affect the breed genetic diversity in the future and suggesting suitable strategies to breed preservation, early on the conservation program. The breeding rate was very low, with a ratio of foaling/live animals of 0.23 (178/760). The estimated number of founders and ancestors contributing to the reference population was 128 and 121. The number of founder herds in the reference population was 64, with an effective number of founder herds for the reference population of 7.6. The mean age of herd owners was 65.50 ± 0.884 years, with a negative association among the herd size and owner's age (P<0.001). In contrast, the size of the herd and the ownership of a male were both positively associated (P<0.001) with the herd number of in-born foals. Both the owners' age and the herd location (RAM home region v. dispersal region) were negatively associated with the foaling number (P<0.001). The main identified risk factors were: low breeding rates; low number of males and their unequal contribution to the genetic pool; unequal contribution of the herds to genetic pool; and advanced age of herd owners.

The mathematical modeling of the lactation curve for dairy traits of the donkey (Equus asinus)

In recent years, an increase in the number of donkeys farmed in Italy as a consequence of the growing demand for donkey milk for direct consumption has been observed. Some research has been carried out on jenny milk composition and on its nutritional properties, whereas milk production features are scarcely described for this species. In this work, the lactation curve shape of donkeys for milk yield and composition was investigated. A total of 453 test-day records for milk yield, fat and protein percentage, and somatic cell count of 62 lactations measured on 46 multiparous jennies of the Ragusano breed were considered. Effects of herd, age, and foaling season were assessed by using a mixed model analysis. Average and individual lactation curves were fitted using the Wood incomplete gamma function, the Cappio-Borlino modified gamma, and a third-order Legendre orthogonal polynomial model. Donkeys foaling between 6- and 10-yr-old had the highest test-day milk yield (about 1.85 kg/d). Donkeys foaling in winter and autumn had a higher daily milk yield compared with those foaling in summer and spring. Less defined results were obtained for composition traits. The general pattern of the donkey lactation curve is similar to the standard shape reported for the main dairy ruminant species, with a peak yield occurring at about 5 wk from parturition. Younger jennies tended to have lower production peaks and higher lactation persistency. Similarly to what is reported for dairy cattle, a large variability in individual patterns has been observed. No differences in goodness of fit have been observed between the models in the case of average lactation curves, whereas orthogonal polynomials were more efficient in fitting individual patterns.