Subchromosomal karyotype evolution in Equidae

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

A Satellite-Free Centromere in Equus przewalskii Chromosome 10

In mammals, centromeres are epigenetically specified by the histone H3 variant CENP-A and are typically associated with satellite DNA. We previously described the first example of a natural satellite-free centromere on Equus caballus chromosome 11 (ECA11) and, subsequently, on several chromosomes in other species of the genus Equus. We discovered that these satellite-free neocentromeres arose recently during evolution through centromere repositioning and/or chromosomal fusion, after inactivation of the ancestral centromere, where, in many cases, blocks of satellite sequences were maintained. Here, we investigated by FISH the chromosomal distribution of satellite DNA families in Equus przewalskii (EPR), demonstrating a good degree of conservation of the localization of the major horse satellite families 37cen and 2PI with the domestic horse. Moreover, we demonstrated, by ChIP-seq, that 37cen is the satellite bound by CENP-A and that the centromere of EPR10, the ortholog of ECA11, is devoid of satellite sequences. Our results confirm that these two species are closely related and that the event of centromere repositioning which gave rise to EPR10/ECA11 centromeres occurred in the common ancestor, before the separation of the two horse lineages.

Robertsonian fusion and centromere repositioning contributed to the formation of satellite-free centromeres during the evolution of zebras

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated to highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and E. asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras E. burchelli (15 of 22) and E. grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At eleven chromosomes, whose primary constrictions seemed to be associated to satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres, therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.

Molecular Dynamics and Evolution of Centromeres in the Genus Equus

The centromere is the chromosomal locus essential for proper chromosome segregation. While the centromeric function is well conserved and epigenetically specified, centromeric DNA sequences are typically composed of satellite DNA and represent the most rapidly evolving sequences in eukaryotic genomes. The presence of satellite sequences at centromeres hampered the comprehensive molecular analysis of these enigmatic loci. The discovery of functional centromeres completely devoid of satellite repetitions and fixed in some animal and plant species represented a turning point in centromere biology, definitively proving the epigenetic nature of the centromere. The first satellite-free centromere, fixed in a vertebrate species, was discovered in the horse. Later, an extraordinary number of satellite-free neocentromeres had been discovered in other species of the genus Equus, which remains the only mammalian genus with numerous satellite-free centromeres described thus far. These neocentromeres arose recently during evolution and are caught in a stage of incomplete maturation. Their presence made the equids a unique model for investigating, at molecular level, the minimal requirements for centromere seeding and evolution. This model system provided new insights on how centromeres are established and transmitted to the progeny and on the role of satellite DNA in different aspects of centromere biology.

Genome collinearity analysis illuminates the evolution of donkey chromosome 1 and horse chromosome 5 in perissodactyls: A comparative study

BACKGROUND

It is important to resolve the evolutionary history of species genomes as it has affected both genome organization and chromosomal architecture. The rapid innovation in sequencing technologies and the improvement in assembly algorithms have enabled the creation of highly contiguous genomes. DNA Zoo, a global organization dedicated to animal conservation, offers more than 150 chromosome-length genome assemblies. This database has great potential in the comparative genomics field.

RESULTS

Using the donkey (Equus asinus asinus, EAS) genome provided by DNA Zoo as an example, the scaffold N50 length and Benchmarking Universal Single-Copy Ortholog score reached 95.5 Mb and 91.6%, respectively. We identified the cytogenetic nomenclature, corrected the direction of the chromosome-length sequence of the donkey genome, analyzed the genome-wide chromosomal rearrangements between the donkey and horse, and illustrated the evolution of the donkey chromosome 1 and horse chromosome 5 in perissodactyls.

CONCLUSIONS

The donkey genome provided by DNA Zoo has relatively good continuity and integrity. Sequence-based comparative genomic analyses are useful for chromosome evolution research. Several previously published chromosome painting results can be used to identify the cytogenetic nomenclature and correct the direction of the chromosome-length sequence of new assemblies. Compared with the horse genome, the donkey chromosomes 1, 4, 20, and X have several obvious inversions, consistent with the results of previous studies. A 4.8 Mb inverted structure was first discovered in the donkey chromosome 25 and plains zebra chromosome 11. We speculate that the inverted structure and the tandem fusion of horse chromosome 31 and 4 are common features of non-caballine equids, which supports the correctness of the existing Equus phylogeny to an extent.

Comparative chromosome painting in Chacoan peccary, Catagonus wagneri

Chacoan peccary (Catagonus wagneri, 2n=20) is the most endangered of three extant species of Tayassuidae. Its karyotype has been studied only by differential chromosome staining methods so far. To establish a comparative cytogenetic map of the peccary, we used cross-species hybridization with porcine (Sus scrofa, 2n=38) painting probes. Painting revealed 30 evolutionary conserved autosomal segments between pig and peccary. The q-arm of the submetacentric chromosome X is homologous to the porcine X chromosome, while the p-arm is composed of heterochromatin. Nucleolar organizer regions were detected on chromosomes 8 and 9 which are homologous to pig chromosomes 8 and 4/18, respectively. Fusions of chromosomes homologous to pig chromosomes 4/7 and 4/18 and fission of chromosome 7 are synapomorphic characters shared by Catagonus wagneri and Tayassu pecari but not by Pecari tajacu. Our results confirmed a high rate of karyotype evolution in Tayassuidae and a closer relationship of Catagonus wagneri with Tayassu pecari than with Pecari tajacu.

Biochemical evidence for diverse strategies in the inner kinetochore

The kinetochore is a complex structure whose function is absolutely essential. Unlike the centromere, the kinetochore at first appeared remarkably well conserved from yeast to humans, especially the microtubule-binding outer kinetochore. However, recent efforts towards biochemical reconstitution of diverse kinetochores challenge the notion of a similarly conserved architecture for the constitutively centromere-associated network of the inner kinetochore. This review briefly summarizes the evidence from comparative genomics for interspecific variability in inner kinetochore composition and focuses on novel biochemical evidence indicating that even homologous inner kinetochore protein complexes are put to different uses in different organisms.

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.

Birth, evolution, and transmission of satellite-free mammalian centromeric domains

Mammalian centromeres are associated with highly repetitive DNA (satellite DNA), which has so far hindered molecular analysis of this chromatin domain. Centromeres are epigenetically specified, and binding of the CENPA protein is their main determinant. In previous work, we described the first example of a natural satellite-free centromere on Equus caballus Chromosome 11. Here, we investigated the satellite-free centromeres of Equus asinus by using ChIP-seq with anti-CENPA antibodies. We identified an extraordinarily high number of centromeres lacking satellite DNA (16 of 31). All of them lay in LINE- and AT-rich regions. A subset of these centromeres is associated with DNA amplification. The location of CENPA binding domains can vary in different individuals, giving rise to epialleles. The analysis of epiallele transmission in hybrids (three mules and one hinny) showed that centromeric domains are inherited as Mendelian traits, but their position can slide in one generation. Conversely, centromere location is stable during mitotic propagation of cultured cells. Our results demonstrate that the presence of more than half of centromeres void of satellite DNA is compatible with genome stability and species survival. The presence of amplified DNA at some centromeres suggests that these arrays may represent an intermediate stage toward satellite DNA formation during evolution. The fact that CENPA binding domains can move within relatively restricted regions (a few hundred kilobases) suggests that the centromeric function is physically limited by epigenetic boundaries.

Improved de novo genomic assembly for the domestic donkey

Donkeys and horses share a common ancestor dating back to about 4 million years ago. Although a high-quality genome assembly at the chromosomal level is available for the horse, current assemblies available for the donkey are limited to moderately sized scaffolds. The absence of a better-quality assembly for the donkey has hampered studies involving the characterization of patterns of genetic variation at the genome-wide scale. These range from the application of genomic tools to selective breeding and conservation to the more fundamental characterization of the genomic loci underlying speciation and domestication. We present a new high-quality donkey genome assembly obtained using the Chicago HiRise assembly technology, providing scaffolds of subchromosomal size. We make use of this new assembly to obtain more accurate measures of heterozygosity for equine species other than the horse, both genome-wide and locally, and to detect runs of homozygosity potentially pertaining to positive selection in domestic donkeys. Finally, this new assembly allowed us to identify fine-scale chromosomal rearrangements between the horse and the donkey that likely played an active role in their divergence and, ultimately, speciation.

Intrachromosomal Rearrangements in Rodents from the Perspective of Comparative Region-Specific Painting

It has long been hypothesized that chromosomal rearrangements play a central role in different evolutionary processes, particularly in speciation and adaptation. Interchromosomal rearrangements have been extensively mapped using chromosome painting. However, intrachromosomal rearrangements have only been described using molecular cytogenetics in a limited number of mammals, including a few rodent species. This situation is unfortunate because intrachromosomal rearrangements are more abundant than interchromosomal rearrangements and probably contain essential phylogenomic information. Significant progress in the detection of intrachromosomal rearrangement is now possible, due to recent advances in molecular biology and bioinformatics. We investigated the level of intrachromosomal rearrangement in the Arvicolinae subfamily, a species-rich taxon characterized by very high rate of karyotype evolution. We made a set of region specific probes by microdissection for a single syntenic region represented by the p-arm of chromosome 1 of Alexandromys oeconomus, and hybridized the probes onto the chromosomes of four arvicolines (Microtus agrestis, Microtus arvalis, Myodes rutilus, and Dicrostonyx torquatus). These experiments allowed us to show the intrachromosomal rearrangements in the subfamily at a significantly higher level of resolution than previously described. We found a number of paracentric inversions in the karyotypes of M. agrestis and M. rutilus, as well as multiple inversions and a centromere shift in the karyotype of M. arvalis. We propose that during karyotype evolution, arvicolines underwent a significant number of complex intrachromosomal rearrangements that were not previously detected.

Pooling strategy and chromosome painting characterize a living zebroid for the first time

We have investigated the complex karyotype of a living zebra-donkey hybrid for the first time using chromosome-specific painting probes produced from flow-sorted chromosomes from a zebra (Equus burchelli) and horse (Equus caballus). As the chromosomes proved difficult to distinguish from one another, a successful new strategy was devised to resolve the difficulty and characterize each chromosome. This was based on selecting five panels of whole chromosome painting probes that could differentiate zebra and donkey chromosomes by labelling the probes with either FITC or Cy3 fluorochromes. Each panel was hybridized sequentially to the same G-Q-banded metaphases and the results combined so that every zebra and donkey chromosome in each suitable metaphase could be identified. A diploid number of 2n = 53, XY was found, containing haploid sets of 22 chromosomes from the zebra and 31 chromosomes from the donkey, without evidence of chromosome rearrangement. This new strategy, developed for the first time, may have several applications in the resolution of other complex hybrid karyotypes and chromosomal aberrations.

The Unique DNA Sequences Underlying Equine Centromeres

Centromeres are highly distinctive genetic loci whose function is specified largely by epigenetic mechanisms. Understanding the role of DNA sequences in centromere function has been a daunting task due to the highly repetitive nature of centromeres in animal chromosomes. The discovery of a centromere devoid of satellite DNA in the domestic horse consolidated observations on the epigenetic nature of centromere identity, showing that entirely natural chromosomes could function without satellite DNA cues. Horses belong to the genus Equus which exhibits a very high degree of evolutionary plasticity in centromere position and DNA sequence composition. Examination of horses has revealed that the position of the satellite-free centromere is variable among individuals. Analysis of centromere location and composition in other Equus species, including domestic donkey and zebras, confirms that the satellite-less configuration of centromeres is common in this group which has undergone particularly rapid karyotype evolution. These features have established the equids as a new mammalian system in which to investigate the molecular organization, dynamics and evolutionary behaviour of centromeres.

Equids

Alongside domestic horses and donkeys, the horse family, also known as equids, comprises six extant wild species of asses and zebras (Figure 1). Equids are extremely well represented in the fossil record, comprising a 55 million-year evolutionary history, punctuated by many episodes of innovation, extinction and migration. Limited to the single genus Equus today, in the Miocene (23.0–5.3 million years ago) the equid family flourished, comprising more than twenty genera. The group originated in Northern America, where the earliest fossil forms have been found, the so-called Hyracotheres, no larger than small dogs. These animals were soft-leaf browsers and in contrast to modern equids, which roam on a single toe with a solid keratin hoof, their hindlimbs were three-toed and their forelimbs four-toed. Equids thus form, together with rhinos and tapirs, the perissodactyls, an order of mammals characterized by an odd number of toes. Unlike ruminants, they are hindgut fermenters, which digest plant cellulose in their intestines and not in differentiated multiple stomach chambers.

Chromosomal polymorphism in mammals: an evolutionary perspective

Although chromosome rearrangements (CRs) are central to studies of genome evolution, our understanding of the evolutionary consequences of the early stages of karyotypic differentiation (i.e. polymorphism), especially the non-meiotic impacts, is surprisingly limited. We review the available data on chromosomal polymorphisms in mammals so as to identify taxa that hold promise for developing a more comprehensive understanding of chromosomal change. In doing so, we address several key questions: (i) to what extent are mammalian karyotypes polymorphic, and what types of rearrangements are principally involved? (ii) Are some mammalian lineages more prone to chromosomal polymorphism than others? More specifically, do (karyotypically) polymorphic mammalian species belong to lineages that are also characterized by past, extensive karyotype repatterning? (iii) How long can chromosomal polymorphisms persist in mammals? We discuss the evolutionary implications of these questions and propose several research avenues that may shed light on the role of chromosome change in the diversification of mammalian populations and species.

Impact of karyotype organization on interlocus recombination between T cell receptor genes in Equidae

The T cell receptor (TCR) genes (TRA, TRB, TRD and TRG) reside in 3 different chromosomal regions. During the maturation of T lymphocytes, the TCR genes are rearranged by site-specific recombination, a process that also predisposes T cells to aberrant rearrangements. Illegitimate recombination between the TCR genes occurs at a low level in healthy individuals, but this frequency may correlate with the risk of lymphoma. The aim of this work was to investigate interlocus recombination in equids. Illegitimate rearrangements were studied in peripheral blood lymphocytes by FISH with painting and BAC probes and by sequencing of PCR products, and the frequencies of recombination were assessed in horses and 4 other equids. The presence of several trans-rearrangement products between the TRA and TRG genes was verified by PCR in all investigated equids. Frequencies of trans-rearrangements in horses are higher than in humans, and colocalization of the TCR genes on the same chromosome increases the incidence of trans-rearrangements between them. The orientation of the TCR genes does not impact interlocus recombination itself but does affect the viability of cells carrying its products and consequently the number of trans-rearrangements observed in lymphocytes.

Speciation with gene flow in equids despite extensive chromosomal plasticity

Horses, asses, and zebras belong to a single genus, Equus, which emerged 4.0-4.5 Mya. Although the equine fossil record represents a textbook example of evolution, the succession of events that gave rise to the diversity of species existing today remains unclear. Here we present six genomes from each living species of asses and zebras. This completes the set of genomes available for all extant species in the genus, which was hitherto represented only by the horse and the domestic donkey. In addition, we used a museum specimen to characterize the genome of the quagga zebra, which was driven to extinction in the early 1900s. We scan the genomes for lineage-specific adaptations and identify 48 genes that have evolved under positive selection and are involved in olfaction, immune response, development, locomotion, and behavior. Our extensive genome dataset reveals a highly dynamic demographic history with synchronous expansions and collapses on different continents during the last 400 ky after major climatic events. We show that the earliest speciation occurred with gene flow in Northern America, and that the ancestor of present-day asses and zebras dispersed into the Old World 2.1-3.4 Mya. Strikingly, we also find evidence for gene flow involving three contemporary equine species despite chromosomal numbers varying from 16 pairs to 31 pairs. These findings challenge the claim that the accumulation of chromosomal rearrangements drive complete reproductive isolation, and promote equids as a fundamental model for understanding the interplay between chromosomal structure, gene flow, and, ultimately, speciation.

Two recessive mutations in FGF5 are associated with the long-hair phenotype in donkeys

Background

Seven donkey breeds are recognized by the French studbook. Individuals from the Pyrenean, Provence, Berry Black, Normand, Cotentin and Bourbonnais breeds are characterized by a short coat, while those from the Poitou breed (Baudet du Poitou) are characterized by a long-hair phenotype. We hypothesized that loss-of-function mutations in the FGF5 (fibroblast growth factor 5) gene, which are associated with a long-hair phenotype in several mammalian species, may account for the special coat feature of Poitou donkeys. To the best of our knowledge, mutations in FGF5 have never been described in Equidae.

Methods

We sequenced the FGF5 gene from 35 long-haired Poitou donkeys, as well as from a panel of 67 short-haired donkeys from the six other French breeds and 131 short-haired ponies and horses.

Results

We identified a recessive c.433_434delAT frameshift deletion in FGF5, present in Poitou and three other donkey breeds and a recessive nonsense c.245G > A substitution, present in Poitou and four other donkey breeds. The frameshift deletion was associated with the long-hair phenotype in Poitou donkeys when present in two copies (n = 31) or combined with the nonsense mutation (n = 4). The frameshift deletion led to a stop codon at position 159 whereas the nonsense mutation led to a stop codon at position 82 in the FGF5 protein. In silico, the two truncated FGF5 proteins were predicted to lack the critical β strands involved in the interaction between FGF5 and its receptor, a mandatory step to inhibit hair growth.

Conclusions

Our results highlight the allelic heterogeneity of the long-hair phenotype in donkeys and enlarge the panel of recessive FGF5 loss-of-function alleles described in mammals. Thanks to the DNA test developed in this study, breeders of non-Poitou breeds will have the opportunity to identify long-hair carriers in their breeding stocks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0065-5) contains supplementary material, which is available to authorized users.