Abstract
The donkey, like the horse, is a promising model for exploring karyotypic instability. We report the de novo whole-genome assemblies of the donkey and the Asiatic wild ass. Our results reflect the distinct characteristics of donkeys, including more effective energy metabolism and better immunity than horses. The donkey shows a steady demographic trajectory. We detected abundant satellite sequences in some inactive centromere regions but not in neocentromere regions, while ribosomal RNAs frequently emerged in neocentromere regions but not in the obsolete centromere regions. Expanded miRNA families and five newly discovered miRNA target genes involved in meiosis may be associated with fast karyotype evolution. APC/C, controlling sister chromatid segregation, cytokinesis and the establishment of the G1 cell cycle phase were identified by analysis of miRNA targets and rapidly evolving genes.
Similar content being viewed by others
Introduction
Donkeys and horses are globally important livestock, representing the Equus genus1,2. Compared with horses, donkeys have superior physiological characteristics, such as a better immune capacity and more effective energy metabolism3,4. The relationship between these species is complicated and confusing. For example, these animals can mate and produce mules or hinnies despite being different species. Karyotypic diversification is more prominent in Equus species than in other mammals5,6, suggesting that the Equus genus is a promising model for exploring the dynamics of chromosomal evolution7. A puzzling phenomenon is the relatively high frequency of centromere repositioning events8 in Equus, as at least seven cases have occurred between donkeys and horses, with at least six further in the donkey9. In our previous study, we investigated the mechanism of chromosomal rearrangement, including Robertsonian translocations and local rearrangements, using de novo assembled genome sequences from Przewalski’s wild horse (Equus przewalskii) and the Mongolian horse (Equus caballus)7. These results suggest that analysis based on whole genome sequences is a delicate and powerful method for studying chromosomal evolution. Here, we report the whole-genome sequence and de novo genome assembly of the donkey and Asiatic wild ass. Using these quality genome sequences, we addressed two questions: (1) what are the demographic and phylogenomic histories accompanying the speciation and genomic adaptive evolution in these representative Equus species and (2) what are the underlying genetic and epigenetic mechanisms of fast karyotype evolution and frequent centromere repositioning.
Results
Genome sequencing, assembly and annotation
The genome of one male donkey was sequenced and de novo assembled using a whole-genome shotgun strategy (Table 1). Eight paired-end libraries (a standard genomic library that was sequenced using paired-end reads with insert sizes of 400–1000 bp), one single-end library (insert size: 1.5–1.9 kbp) and eight mate-paired libraries (insert sizes: 3–15 kbp) were constructed for genome sequencing (Supplementary Table 1). Paired-end libraries were sequenced using the Illumina Miseq platform, the single-end library was sequenced using the Roche 454 FLX+ platform and mate-paired libraries were sequenced using the Illumina Hiseq2000 platform. The total sequence coverage was approximately 42.4-fold (genome size: ∼2.36 Gb) (Supplementary Table 2). For the Asiatic wild ass, one paired-end library (insert size: 500 bp) was constructed and sequenced using the Illumina Hiseq2000 platform (12.1-fold; Supplementary Tables 3 and 4). High contiguity genome sequences from the donkey were generated after their de novo assembly and they consisted of 2,166 scaffolds (>1k bp) with a total size of 2.36 Gb (Supplementary Table 5). The N50 lengths of the contigs and scaffolds were 66.7 kb and 3.8 Mb, respectively. Compared with other previously published genome sequences10,11,12,13 (Supplementary Figs 1 and 2), the contiguity of the contigs in the donkey assemblies was better. We also validated 248 core eukaryotic genes14 in the donkey genome assemblies and found considerable completeness (Supplementary Table 6). These improvements may be the result of longer lengths for the sequence reads, because the reads used in this study were mainly generated by the Illumina Miseq platform (2 × 251 bp) and longer than those generated by Hiseq2000 (2 × 100 bp)10,11,12,13.
To improve our gene prediction accuracy, eight types of tissue samples (heart, liver, spleen, lung, kidney, brain, spinal cord and muscle) from another female donkey were used to construct a normalized cDNA library. RNA-seq was performed using the Roche 454 FLX+ platform and 1,390,416 reads were generated with an average length of 522 bp (Supplementary Fig. 3 and Supplementary Table 7). Donkey genome annotation was performed using a dexterous genome annotation pipeline, including both ab initio predictions (Augustus and SNAP)15,16 and homology-based methods (RNA-seq of the female donkey and homologous proteins sequences of the Thoroughbred horse17). A total of 23,214 protein-coding genes were predicted in the donkey genome (Table 1, Supplementary Figs 4 and 5) averaging 1,281 bp coding sequences (CDSs) per gene. Among these genes, 15,648 could be confirmed with the RNA-seq sequences (Supplementary Fig. 6).
Demographic history and phylogenetic analysis
We identified 2,187,070 and 3,321,087 heterozygous SNPs (within each individual) in the donkey and the Asiatic wild ass genomes, respectively (Supplementary Table 8). The rate of heterozygosity was considerably higher in the Asiatic wild ass than in donkey. We also reconstructed the donkey, Asiatic wild ass and horse population demographics over the last one million years (Fig. 1a). Because Thoroughbred horse (Twilight) pedigrees show substantial levels of inbreeding17, we used heterozygous SNPs from the Mongolian horse7. Our demographic analysis revealed three horse population bottlenecks, which is consistent with the quaternary glaciations. Similar to those of the horse, Asiatic wild ass lineages show extremely dynamic demographic trajectories. Interestingly, the size of the donkey population was steady. We believe that this stability is because the donkey ancestors (African wild asses) living in northeast Africa2 may have been influenced by different climates during the quaternary glaciations, as climate changes could result in grassland contraction or expansion18.
The rich Equus fossil records have made this genus a model for evolutionary processes19. Previous research has shown that the donkey and the horse shares common ancestors approximately 6.4–12.7 million years ago20,21,22. In this paper, we constructed a phylogenetic tree using 5,665 single-copy orthologs from nine species17,23,24,25,26,27,28 (Fig. 1b). As shown in this tree, the Asiatic wild ass is most closely related to the donkey and together they form a sister group with the horse. Our results show that the donkey separated from the horse lineage approximately 7.7–15.4 million years ago, whereas the donkey and the Asiatic wild ass diverged approximately 1.5–3.3 million years ago. These estimates are comparable to the earliest divergence times reported20,21,22.
Genetic evolution
To obtain greater insight into the evolutionary dynamics of these genes, we calculated the expansion and contraction of orthologous gene clusters between the donkey and the horse. A total of 283 gene families in the donkey showed significant expansion (P < 0.05) compared with 206 in the horse (Fig. 1b). The functional categories that were enriched in significant donkey gene family expansions included olfactory transduction (KEGG:map04740, p = 5.355e–08) and protein digestion and absorption (KEGG:map04974, p = 0.01327) (Supplementary Table 9). The horse gene family expansions were primarily associated with defense responses (GO:0006952, p = 0.011693853) and responses to stress (GO:0006950, p = 0.028023107) (Supplementary Table 10).
Rapidly evolving genes are one of the primary contributors to such functional changes. We identified 1,292 genes evolving significantly (p < 0.05) faster in the donkey than in the horse and 706 genes evolving significantly (p < 0.05) faster in the horse than in the donkey. Rapidly evolving genes in domestic donkeys are significantly associated with aerobic respiration (GO:0009060, p = 0.027964968), forebrain development (GO:0030900, p = 0.006710136), regulation of lymphocyte differentiation (GO:0045619, p = 0.024669795), the tricarboxylic acid cycle (GO:0006099, p = 0.016761429) and the acetyl-CoA catabolic process (GO:0046356, p = 0.016761429) (Fig. 1c, Supplementary Table 11). These changes may be correlated with more effective energy metabolism4 and improved immune capacity in donkeys compared with horses. More specifically, twenty genes that are associated with forebrain development were found to be rapidly evolving in the donkey. Also, cell cycle arrest (GO:0007050, p = 0.020903) and telomere maintenance (GO:0000723, p = 0.002098) are rapidly evolving in the donkey genome, which may be associated with rapid karyotypic evolution. In contrast, rapidly evolving genes in the horse are significantly enriched in second-messenger-mediated signaling (GO:0019932, p = 0.001296105), heart looping (GO:0001947, p = 0.046821101), neural tube patterning (GO:0021532, p = 0.041550593), photoreceptor cell maintenance (GO:0045494, p = 0.041550593) and ribosome biogenesis (GO:0042254, p = 0.035248368) (Fig. 1c, Supplementary Table 12). These results may be associated with the animated disposition and greater athletic ability of the horse.
Synteny analysis and repetitive sequences
Dramatic chromosomal rearrangement in Equus individuals is a notable characteristic compared with other mammals5,6. However, genome-wide rearrangements between the donkey and horse have not been characterized given that donkeys have a different number of chromosomes (2n = 62) than horses (2n = 64)5. We performed whole-genome synteny analysis between the donkey and Thoroughbred horse genomes. A collinearity region between the donkey and Thoroughbred horse was approximately 1.89 Gb (Fig. 2a). Four types of rearrangements, BRK (insertion of unknown origin), DUP (inserted duplication), INV (inversion) and JMP (relocation) were identified. Rearrangement of the donkey genome was particularly evident when donkey-Thoroughbred horse genome alignments were compared with those of Thoroughbred horse-wild horse and Thoroughbred horse-Mongolian horse, as more large-scale chromosomal rearrangements can be found in the donkey genome.
Previous research has indicated that repetitive sequences are associated with syntenic breakpoints and chromosomal fragility. Seven types of common repetitive sequences in the donkey genome were identified: short interspersed repeated sequences (SINEs), long interspersed repeated sequences (LINEs), long terminal repeats (LTRs), DNA elements, satellites, simple repeats and low complexity. Overall, analyses of these sequences indicated that 42% of the donkey genome sequences are repetitive sequences (Supplementary Table 13), which is comparable to the horse (41.4%). Satellite sequences comprise 0.05% of the donkey genome, which is considerably lower than in the horse (1.59%). Satellite sequences are often associated with centromeres26, including some that are new in the donkey. The proportions of LINE_L1 and LTR_ERV1 increased, but those of LINE_L2 and several other repetitive sequences decreased in breakpoint regions (Fig. 2b). This phenomenon is more evident in the donkey genome, which is consistent with our previous findings7.
Chromosome rearrangements and sequence signatures in centromere regions
A striking phenomenon in Equus is the relatively high frequency of centromere repositioning events9. Although such events provide a potentially powerful evolutionary force for reproductive isolation and speciation, the underlying mechanisms remain unclear29. Comparative FISH studies have found that at least seven different centromere repositioning events occurred between the donkey and horse and at least six further occurred in the donkey alone9. Based on the quality of whole-genome donkey sequences, we were able to perform microscopic analyses across the normal centromere regions of the horse (Thoroughbred) as well as neocentromere regions and inactive centromere regions in donkeys. Using the same probes9 and two major Equus satellite sequences30 as in previous studies, we identified the centromere regions of donkey and horse chromosomes. Six types of regions were categorized into seven pairs of chromosomes in donkeys and horses (Fig. 3a, Supplementary Table 14) including the following: (1) centromere regions in horse chromosomes (“region #1” hereafter), (2) centromere regions in donkey chromosomes (at least six centromeres are neocentromeres, region #2), (3) homologous regions in horse chromosomes related to region #2 (region #3), (4) homologous regions in donkey chromosomes related to region #1 (region #4), (5) other regions in horse chromosomes (region #5), (6) other regions in donkey chromosomes (region #6).
Any shift in centromeric function without chromosomal rearrangement can be considered centromere repositioning30. However, some researchers have noted that centromere regions are hot spots for chromosomal changes in evolution and disease31. We explored chromosomal rearrangements using synteny analysis based on genomic sequences between homologous regions (Fig. 3a) (region #1 vs. region #4, region #3 vs. region #2, region #5 vs. region #6). Increased rearrangements were detected in several chromosomes (Fig. 3b,c) when region #1 and region #4 (representing inactive centromere regions) were compared and these were potentially caused by an accumulation or loss of satellite sequences. In contrast, no obvious increased rearrangements were detected (Fig. 3b,c) when region #3 and region #2 (representing neocentromere regions) were compared. Thus, we confirmed previous study results indicating that a DNA fragment can acquire centromere function without sequence alteration30,32.
To date, no prominent sequence characteristics have been confirmed to promote centromere repositioning, although it is widely accepted that neocentromeres can gradually accumulate satellite sequences accompanied by centromerization29. We noticed that the content of the satellite sequences in region #1 and region #4 in several chromosomes was increased (Fig. 3d,e, Supplementary Fig. 7). Twenty types of satellite sequences were examined in this study (Supplementary Fig. 8). Because region #1 of ECA14, ECA20, ECA22, ECA26 and region #4 of EAS8, EAS15 contained abundant SAT2p30 (Supplementary Fig. 8), we believe that the SAT2p were accumulated in the process of centromerization. Region #1 of ECA6, ECA11 and ECA17 does not contain abundant SAT2p, indicating these three centromeres may be novel. Unexpectedly, ribosomal RNAs were discovered in neocentromere regions and their homologous regions. In contrast, no ribosomal RNA could be detected in region #4 (inactive centromere regions) (Fig. 3d,f). A neighbor-joining tree constructed using conservative 5SrRNAs revealed that these 5SrRNAs are closely related (Supplementary Fig. 9). It is particularly worth mentioning that the ribosomal RNAs are located in the fibrillar centers of the nucleolus and play an important role in the organization of the nucleolus33. More research is needed to explain the genetic association between ribosomal RNA, nucleoli and centromere repositioning.
Small RNA-seq and prediction of novel miRNA targets
To understand the role of epigenetic regulation in karyotypic evolution29, we annotated the non-coding RNAs and analyzed differentially expanded miRNA families in the donkey vs. other mammalian species. A total of 1198 miRNAs, 512 snoRNAs, 530 snRNAs and 189 lncRNAs were identified in the donkey genome (Supplementary Table 15). The number of miRNAs in the donkey genome was comparable to humans (1215), mice (1497) and was higher than horses (881), dogs (647) and cattle (494)34. However, the distribution of miRNAs in the different miRNA families was quite different between donkeys and other mammals (Fig. 4a). We found that several miRNA families were expanded in the donkey genome (Fig. 4a, c). Targeted genes in donkey expanded miRNA families that were significantly enriched were related to the cell cycle, cancer and oocyte meiosis, which are probably associated with fast karyotype evolution (Supplementary Table 16). From the small RNA libraries that were constructed from the above-mentioned eight types of tissue samples (Supplementary Table 17), we identified 118 miRNA families matching those in the existing miRNA database and 40 novel miRNAs specific to the donkey (Supplementary Table 18). Five of the newly discovered donkey miRNAs target genes are involved in meiosis (Fig. 4b,c, Supplementary Table 19), suggesting fast karyotype evolution. In the meiosis pathway, another five genes are rapidly evolving (Fig. 4c). APC/C, which controls sister chromatid segregation, cytokinesis and establishment of the G1 phase of the cell cycle, was identified by analysis of miRNA target and rapidly evolving genes35,36.
Discussion
The donkey and Asiatic wild ass genomes supplement the reference genome for the Equus genus. Our comparative analysis based on these genomic sequences provides important insight into the demographic history and adaptive evolution of Equus. In addition, these results enhance our understanding of the chromosomal rearrangements and dynamics of characteristic sequences associated with centromere repositioning. These data will be beneficial to future research of the genomics of the Equus genus and mammalian chromosomal evolution.
Methods
Sampling and genome sequencing
All animal care and research procedures were carried out in accordance with the guiding principles for the care and use of laboratory animals and were approved by the Institutional Animal Care and Use Committee at Inner Mongolia Agricultural University. For donkey genome sequencing, a 7-year-old male donkey was selected from the Xilingol League of Inner Mongolia, China on 18 February 2010. For Asiatic wild ass genome sequencing, approximately 5 ml of blood from a male Asiatic wild ass was provided by the Bayan Nur Forestry Administration. The blood sample was collected during veterinary exams for several Asiatic wild asses on 1 March 2002. No Asiatic wild ass was hurt or captured as a result of these studies. DNA was extracted from peripheral blood cells. Eight paired-end libraries (insert sizes: 400, 450, 700 and 1000 bp), one single-end library (insert size 1.5–1.9 kb) and eight mate-paired libraries (insert sizes: 3, 5, 8, 12 and 15 kb) were constructed for donkey genome sequencing. Paired-end libraries were sequenced using the Illumina Miseq platform (2 × 251 bp), the single-end library was sequenced using the Roche 454 FLX+ platform (average: 510 bp) and Mate-paired libraries were sequenced using the Illumina Hiseq2000 platform (2 × 100 bp). For the Asiatic wild ass, one paired-end library (insert size 500 bp) was constructed and sequenced using the Illumina Hiseq2000 platform (2 × 100 bp). Library preparation and sequencing followed the manufacturer’s instructions.
Data filtering
Cutadapt1.2.1 (https://pypi.python.org/pypi/cutadapt/1.2.1) was used to trim adapter sequences from sequence reads generated by Illumina Miseq and Hiseq2000. Low-quality reads and reads with potential sequencing errors were also eliminated. For reads generated by Illumina Miseq, if the average phred quality scores for five consecutive bases were <Q20, we trimmed reads from the 3′-end. For reads generated by Illumina Hiseq2000, if the average phred quality scores of five consecutive bases were <Q20, we removed this read and its matching sequence.
Donkey genome assembly
We first assembled the sequence reads of the pair-end and single-end libraries into contigs and scaffolds using Newbler v2.8. Then, we used SSPACE software37 and information for the mate-pair libraries to construct longer scaffolds. Finally, Gapcloser (http://soap.genomics.org.cn/soapdenovo.html) was used to close gaps inside scaffolds.
Repetitive sequence and noncoding sequences analysis
RepeatMasker (http://www.repeatmasker.org/) was used to identify interspersed repeats and low complexity DNA sequences from the donkey and Thoroughbred horse genomes. Twenty types of satellite sequences were then plotted in a “heat map” using R software. 5SrRNA sequences were used to build the neighbor-joining tree using MEGA638. Genome noncoding sequence annotation was used in the Rfam database34. A small RNA library was constructed from eight types of tissue samples (heart, liver, spleen, lung, kidney, brain, spinal cord and muscle) from another female donkey. The library was sequenced using the Miseq platform. For novel miRNA identification, mireap (http://mireap.sourceforge.net/) was used. For target gene annotation, Miranda39 was used.
Genome annotation and RNA-seq
Donkey genome annotation was performed using the MAKER40 annotation pipeline, which included ab initio predictions and homology-based methods. Ab initio predictions were performed using Augustus15 and SNAP16. cDNA data were generated from multiple RNA sources. cDNA libraries were constructed from eight types of tissue samples (heart, liver, spleen, lung, kidney, brain, spinal cord and muscle) from another female donkey. The libraries were sequenced using the Roche 454 FLX+ platform. Homology-based prediction was performed by blasting against homologous protein sequences of Thoroughbred horse17 and cDNA sequences from donkey.
Heterozygosity rate and demographic history
Qualified sequence reads from pair-end libraries of the donkey and Asiatic wild ass were mapped to the scaffolds of the donkey. SNPs and InDels were called using the Genome Analysis Toolkit41 following its manual. We flagged a candidate SNP as a likely false-positive if it exhibited the following characteristics: (1) sequence coverage at that point is more than 200 or less than 4; (2) HaplotypeScore >13.0, MQ < 40, QD < 2; (3) ReadPosRankSum <−8.0, MQRankSum <−12.5. The demographic histories of the donkey, Asiatic wild ass and Mongolian horse were inferred using “pairwise sequentially Markovian coalescence” (PSMC)42 based on SNP distribution. Parameters were set as follows: −N30 −t15 −r5 −p 4+25*2+4+6. The Equus generation time (g) = 5 years and the neutral mutation rate per generation (μ) = 2.5 × 10−8 were set. Because low sequence coverage (below 20-fold) deeply impacted PSMC inference42, we performed a correction for Asiatic wild ass assuming a uniform False Negative Rate (uNFR = 26%) reported in previous research19.
Phylogeny analysis
Protein-coding genes from seven mammalian species (opossum, dog, pig, cattle, Thoroughbred horse, mouse and human) downloaded from Ensembl (http://www.ensembl.org) were used in addition to donkey genes to define gene families by OrthoMCL43. Thereafter, 5,665 single-copy families, which were generated from this analysis, were used to reconstruct phylogenies and estimate the time points of divergence. Protein-coding gene sequences from the Asiatic wild ass were generated by mapping reads from the Asiatic wild ass to the scaffolds of the donkey with samtools44 and genBlastG45. The protein sequences of orthologous gene sets were aligned by MUSCLE46 with its default settings. Poor alignment sites were eliminated using Gblock47. The phylogeny tree (including nine species) was drawn by PhyML48 using the JTT model. Based on the reconstructed phylogeny tree, we estimated the evolutionary time scales by PAML49. Calibration times were queried from the TimeTree database (http://www.timetree.org).
Gene family expansion and contraction
Gene families were defined by OrthoMCL43. Gene family expansion analysis was performed by CAFE50 based on a reconstructed phylogeny tree.
Rapidly evolving genes and dN/dS analysis
This analysis utilized 6,771 1:1 orthologous genes from seven species (donkey, Thoroughbred horse, dog, pig, cattle, mouse and human). The protein sequences from orthologous gene sets were aligned by MUSCLE46 using default settings. Gblock47 was used to eliminate poor alignment sites. Afterward, dN/dS ratios for each gene were estimated with the codeml function in the PAML package49. The maximum-likelihood method was used to estimate dN (the rate of non-synonymous substitution), dS (the rate of synonymous substitution) and dN/dS (the ratio of non-synonymous substitutions to the rate of synonymous substitutions). The likelihood ratio test was used to evaluate the p-value for each gene.
Synteny analysis and SV calling
We used Mauve Contig Mover51 to order donkey genome drafts relative to the Thoroughbred horse genome. Then, we used MUMmer52 to perform whole-genome synteny analysis. Genome rearrangements were identified using the nucmer module. The parameter was Options “-c 800 -g 300–l,100”.
Additional Information
Accession codes: The Whole Genome Shotgun project has been deposited in DDBJ/EMBL/GenBank as project accession PRJNA200657 and PRJNA200654 of donkey and Asiatic wild ass, respectively. The genome assembly of donkey has been deposited at DDBJ/EMBL/GenBank under the accession JREZ00000000 and this version described in this paper is version JREZ01000000. Transcript sequencing data have been deposited under Short Read Archive (SRA) accession SRR1562259. Asiatic wild ass genome sequencing data have been deposited under Short Read Archive (SRA) accession SRR1562345.
How to cite this article: Huang, J. et al. Donkey genome and insight into the imprinting of fast karyotype evolution. Sci. Rep. 5, 14106; doi: 10.1038/srep14106 (2015).
Change history
18 December 2015
A correction has been published and is appended to both the HTML and PDF versions of this paper. The error has not been fixed in the paper.
References
Rossel, S. et al. Domestication of the donkey: timing, processes and indicators. Proc. Natl. Acad. Sci. USA 105, 3715–3720 (2008).
Beja-Pereira, A. et al. African origins of the domestic donkey. Science 304, 1781 (2004).
Hintz, H. F. & Cymbaluk, N. F. Nutrition of the horse. Annu. Rev. Nutr. 14, 243–267 (1994).
Izraely, H., Choshniak, I., Stevens, C. E., Demment, M. W. & Shkolnik, A. Factors determining the digestive efficiency of the domesticated donkey (Equus asinus asinus). Q. J. Exp. Physiol. 74, 1–6 (1989).
Trifonov, V. A. et al. Multidirectional cross-species painting illuminates the history of karyotypic evolution in Perissodactyla. Chromosome Res. 16, 89–107 (2008).
Bush, G. L., Case, S. M., Wilson, A. C. & Patton, J. L. Rapid speciation and chromosomal evolution in mammals. Proc. Natl. Acad. Sci. USA 74, 3942–3946 (1977).
Huang, J. et al. Analysis of horse genomes provides insight into the diversification and adaptive evolution of karyotype. Sci. Rep. 4, 4958 (2014).
Montefalcone, G., Tempesta, S., Rocchi, M. & Archidiacono, N. Centromere repositioning. Genome Res. 9, 1184–1188 (1999).
Carbone, L. et al. Evolutionary movement of centromeres in horse, donkey and zebra. Genomics 87, 777–782 (2006).
Cho, Y. S. et al. The tiger genome and comparative analysis with lion and snow leopard genomes. Nat. Commun. 4, 2433 (2013).
Wan, Q. H. et al. Genome analysis and signature discovery for diving and sensory properties of the endangered Chinese alligator. Cell Res. 23, 1091–1105 (2013).
Wang, Z. et al. The draft genomes of soft-shell turtle and green sea turtle yield insights into the development and evolution of the turtle-specific body plan. Nat. Genet. 45, 701–706 (2013).
Zhou, X. et al. Baiji genomes reveal low genetic variability and new insights into secondary aquatic adaptations. Nat. Commun. 4, 2708 (2013).
Parra, G., Bradnam, K. & Korf, I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23, 1061–1067 (2007).
Stanke, M. & Waack, S. Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics 19 Suppl 2, ii215–225 (2003).
Korf, I. Gene finding in novel genomes. BMC Bioinformatics 5, 59 (2004).
Wade, C. M. et al. Genome sequence, comparative analysis and population genetics of the domestic horse. Science 326, 865–867 (2009).
Lorenzen, E. D. et al. Species-specific responses of Late Quaternary megafauna to climate and humans. Nature. 479, 359–364 (2011).
Orlando, L. et al. Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse. Nature 499, 74–78 (2013).
Huelsenbeck, J. P., Larget, B. & Swofford, D. A compound poisson process for relaxing the molecular clock. Genetics 154, 1879–1892 (2000).
Nikaido, M. et al. Maximum likelihood analysis of the complete mitochondrial genomes of eutherians and a reevaluation of the phylogeny of bats and insectivores. J. Mol. Evol. 53, 508–516 (2001).
Waddell, P. J., Cao, Y., Hasegawa, M. & Mindell, D. P. Assessing the Cretaceous superordinal divergence times within birds and placental mammals by using whole mitochondrial protein sequences and an extended statistical framework. Syst. Biol. 48, 119–137 (1999).
Elsik, C. G. et al. The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science 324, 522–528 (2009).
Groenen, M. A. et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature 491, 393–398 (2012).
Kirkness, E. F. et al. The dog genome: survey sequencing and comparative analysis. Science 301, 1898–1903 (2003).
Lander, E. S. et al. Initial sequencing and analysis of the human genome. Nature 409, 860–921 (2001).
Mikkelsen, T. S. et al. Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences. Nature 447, 167–177 (2007).
Waterston, R. H. et al. Initial sequencing and comparative analysis of the mouse genome. Nature 420, 520–562 (2002).
Amor, D. J. et al. Human centromere repositioning “in progress”. Proc. Natl. Acad. Sci. USA 101, 6542–6547 (2004).
Piras, F. M. et al. Uncoupling of satellite DNA and centromeric function in the genus Equus. PLoS Genet. 6, e1000845 (2010).
Villasante, A., Abad, J. P. & Mendez-Lago, M. Centromeres were derived from telomeres during the evolution of the eukaryotic chromosome. Proc. Natl. Acad. Sci. USA 104, 10542–10547 (2007).
Ventura, M. et al. Recurrent sites for new centromere seeding. Genome Res. 14, 1696–1703 (2004).
O’Sullivan, J. M., Pai, D. A., Cridge, A. G., Engelke, D. R. & Ganley, A. R. The nucleolus: a raft adrift in the nuclear sea or the keystone in nuclear structure? Biomol. Concepts 4, 277–286 (2013).
Burge, S. W. et al. Rfam 11.0: 10 years of RNA families. Nucleic. Acids. Res. 41, D226–232 (2013).
Chang, L., Zhang, Z., Yang, J., McLaughlin, S. H. & Barford, D. Molecular architecture and mechanism of the anaphase-promoting complex. Nature 513, 388–393 (2014).
Bharadwaj, R. & Yu, H. The spindle checkpoint, aneuploidy and cancer. Oncogene 23, 2016–2027 (2004).
Boetzer, M., Henkel, C. V., Jansen, H. J., Butler, D. & Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics 27, 578–579 (2011).
Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013).
Enright, A. J. et al. MicroRNA targets in Drosophila. Genome Biol. 5, R1 (2003).
Cantarel, B. L. et al. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 18, 188–196 (2008).
McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
Li, H. & Durbin, R. Inference of human population history from individual whole-genome sequences. Nature 475, 493–496 (2011).
Li, L., Stoeckert, C. J., Jr. & Roos, D. S. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189 (2003).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
She, R. et al. genBlastG: using BLAST searches to build homologous gene models. Bioinformatics 27, 2141–2143 (2011).
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic. Acids. Res. 32, 1792–1797 (2004).
Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540–552 (2000).
Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).
Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586–1591 (2007).
De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269–1271 (2006).
Darling, A. E., Mau, B. & Perna, N. T. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 5, e11147 (2010).
Delcher, A. L. et al. Alignment of whole genomes. Nucleic. Acids. Res. 27, 2369–2376 (1999).
Acknowledgements
This work was supported by National Natural Science Foundation of China (31472070), Ministry of Science and Technology of the People’s Republic of China specific scientific and technological cooperation with Russia (2011DFR30860), Inner Mongolia key laboratory project (20130902), National Natural Science Foundation of China (31160446) and National Natural Science Foundation of China (31360538). We thank Christa F. Honaker for valuable comments to this manuscript.
Author information
Authors and Affiliations
Contributions
M.D., H.M., Y.Z. and Z.S. designed and managed the project. D.D., W.B., Z.X., A.A., T.B. and J.H. collected samples and prepared the nucleic acid samples. B.L., L.Y., J.W., W.B., X.R., B.J., Q.Z., A.L., S.B., Y.G., R.W., W.Z., S.L., J.D., Y.Z., J.H. and F.D. performed the genomes sequencing. J.H., Y.Z., D.B., W.S., W.Y., Y.B., H.H., H.B. and Y.B. designed and performed the genome assembly and analyses. J.H., Y.S., H.M., Y.Z. and M.D. wrote and revised the paper.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Huang, J., Zhao, Y., Bai, D. et al. Donkey genome and insight into the imprinting of fast karyotype evolution. Sci Rep 5, 14106 (2015). https://doi.org/10.1038/srep14106
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep14106
This article is cited by
-
Genome collinearity analysis illuminates the evolution of donkey chromosome 1 and horse chromosome 5 in perissodactyls: A comparative study
BMC Genomics (2021)
-
Evaluation of the cell mediated immune responses using quantitative Real-time Polymerase Chain Reaction during the infestation of Rhinoestrus usbekistanicus (Diptera: Oestridae) in equine
International Journal of Tropical Insect Science (2021)
-
Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae
BMC Genomics (2020)