http://orcid.org/0000-0002-8312-0483
Figures
Abstract
The Cervidae family comprises more than fifty species divided into three subfamilies: Capreolinae, Cervinae and Hydropotinae. A characteristic attribute for the species included in this family is the great karyotype diversity, with the chromosomal numbers ranging from 2n = 6 observed in female Muntiacus muntjak vaginalis to 2n = 70 found in Mazama gouazoubira as a result of numerous Robertsonian and tandem fusions. This work reports chromosomal homologies between cattle (Bos taurus, 2n = 60) and nine cervid species using a combination of whole chromosome and region-specific paints and BAC clones derived from cattle. We show that despite the great diversity of karyotypes in the studied species, the number of conserved chromosomal segments detected by 29 cattle whole chromosome painting probes was 35 for all Cervidae samples. The detailed analysis of the X chromosomes revealed two different morphological types within Cervidae. The first one, present in the Capreolinae is a sub/metacentric X with the structure more similar to the bovine X. The second type found in Cervini and Muntiacini is an acrocentric X which shows rearrangements in the proximal part that have not yet been identified within Ruminantia. Moreover, we characterised four repetitive sequences organized in heterochromatic blocks on sex chromosomes of the reindeer (Rangifer tarandus). We show that these repeats gave no hybridization signals to the chromosomes of the closely related moose (Alces alces) and are therefore specific to the reindeer.
Citation: Frohlich J, Kubickova S, Musilova P, Cernohorska H, Muskova H, Vodicka R, et al. (2017) Karyotype relationships among selected deer species and cattle revealed by bovine FISH probes. PLoS ONE 12(11): e0187559. https://doi.org/10.1371/journal.pone.0187559
Editor: Roscoe Stanyon, University of Florence, ITALY
Received: August 23, 2017; Accepted: October 21, 2017; Published: November 7, 2017
Copyright: © 2017 Frohlich et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper.
Funding: This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic under the project CEITEC 2020 (LQ1601) and by the Ministry of Agriculture of the Czech Republic (RO 0516). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: AAL, moose (Alces alces); AAM, pronghorn (Antilocapra americana); BAC, Bacterial Artificial Chromosome; BTA, cattle (Bos taurus); CAX, Cetartiodactyla ancestral X; CCA, roe deer (Capreolus capreolus); CEL, red deer (Cervus elaphus); CLI, pygmy hippo (Choeropsis/Hexaprotodon liberiensis); CTR, rusa deer (Cervus timorensis russa); DDA, fallow deer (Dama dama); dist, distal in relation to the centromere; EDA, milu deer (Elaphurus davidianus); FISH, fluorescence in situ hybridization; GCA, giraffe (Giraffa camelopardalis); MRE, Chinese muntjac (Muntiacus reevesi); OJO, okapi (Okapia johnstoni); PAK, Pecoran ancestral karyotype; PAR, pseudoautosomal region; PAX, Pecoran ancestral X; prox, proximal in relation to the centromere; RAK, Ruminantia ancestral karyotype; REL, Eld's deer (Rucervus eldii); RTA, reindeer (Rangifer tarandus)
Introduction
The Cervidae family includes more than fifty different species which are divided into three subfamilies: Capreolinae, Cervinae and Hydropotinae [1]. Several of the Cervidae species have a growing economic potential as farm animals, but they are also hunted for trophies. However, most species are threatened and need protection. A distinctive trait of this family is the great diversity of karyotypes amongst respective species, with the chromosomal numbers (2n) ranging from 6 in female Muntiacus muntjak vaginalis to 70 in Mazama gouazoubira [2]. The chromosomal evolution in particular cervid subfamilies occurred by different types of rearrangements. For example, Robertsonian fusions are the most common type of chromosomal rearrangements in Cervini, whereas tandem fusions are the major chromosomal rearrangements underlying the karyotype diversification of muntjacs (Muntiacini) [2–5].
The karyotypes of most Cervidae species have been predominantly studied by standard cytogenetic methods (G-, R-banding, etc.) [3,6,7] and the evolutionary relationships have been studied by comparative chromosome painting in a limited number of cervid species. Cross-species chromosome painting is a variation of the fluorescence in-situ hybridization technique (FISH) and enables us to establish chromosome homology maps, define the sites of chromosome fusions and fissions and to investigate chromosome rearrangements which occurred during karyotype evolution [8,9].
The Capreolinae subfamily consists of nine genera (predominantly with 2n = 70) where chromosomes in only two species were characterized by comparative chromosome painting. Dementeyeva et al. have established the chromosome map of C. pygargus by cross-species chromosome painting using dromedary probes (Camelus dromedarius) [10] and homologies between Mazama gouazoubira and Chinese muntjac (Muntiacus reevesi) have been established by hybridization with set of Chinese muntjac painting probes [11].
The subfamily Cervinae comprises two tribes Cervini (seven genera) and Muntiacini (two genera). Whereas only Robertsonian translocations were found in karyotypes of the Cervini species, repeated tandem fusions took part in karyotype differentiation of the Muntiacini [2,12]. Generally as the tandem fusions in the karyotype evolution are scarce, the karyotypes of all Muntiacini species were intensively studied by molecular cytogenetic methods based on FISH with painting and BAC probes [4,13,14]. High resolution cross-species comparative mapping between all Muntiacini species were performed, while only two Cervini species were analyzed by chromosomal painting [12].
The subfamily Hydropotinae includes only one species (Hydropotes inermis), the most primitive deer, which has retained the karyotype (2n = 70) closely related to the ancestral karyotype of Cervidae. This fact was corroborated by hybridization of painting probes prepared from M. reevesi [11].
Although the family Bovidae is evolutionally closely related to the Cervidae, cattle painting probes have never been applied to deer chromosomes. The cattle karyotype (2n = 60) is composed of 58 acrocentric autosomes and 2 sex chromosomes and differs only by 1 fission from the Pecoran ancestral karyotype (PAK, 2n = 58)[15,16]. The cattle (Bos taurus) fluorescent painting probes have been exploited for comparative karyotype studies not only in many closely related wild animal species [17,18], but also in more phylogenetically distant Cetartiodactyla families [19,20].
Evolutionary rearrangements of ruminant autosomes are rare, except for centric and tandem fusions [18,21–23]. Only the use of high-density BAC (Bacterial Artificial Chromosome) mapping enables precise detection of intrachromosomal rearrangements and their breakpoints in different species. Changes in homologous chromosomes were detected for example in BTA3 homolog between Bovini and Antilopini [24] or in BTA1 homolog between Bovidae and Cervidae, Giraffidae and Antilocapridae [25,19]. An alternative to the high resolution BAC mapping is the in silico bioinformatic approach, which allows to detect cryptic chromosomal divergencies. However, this strategy is limited only to the species with sequenced and well assembled genomes [26].
Unlike autosomes, the sex chromosomes differ by more complex chromosomal rearrangements including inversions, centromere shift, heterochromatic variation, and autosomal translocations [22,27–30]. Interspecies variation of X chromosome provides a rich source of phylogenetic information but detailed structure of chromosome X has not been studied in more detail in most Ruminantia species for some considerable time. Just recently, a new investigation of the X chromosome evolution in different representatives of Cetartiodactyla employing high resolution BAC mapping was published [31]. In Cervidae, comparative FISH with whole X chromosome painting probes was employed in several studies [10,12,32]. Better insights into the organization of X chromosomes are enabled by using detailed mapping comparison of BAC or cosmid clones which were used for analysis of the X chromosome rearrangements in various species of the family Bovidae. For instance, Gallagher et al. confirmed the conservation of X chromosome homology regions among members of Bovini, Boselaphini, Tragelaphini, domestic sheep and two deer species using BAC probes [27]. Comparative mapping of three X chromosome-specific bovine cosmids was employed in Rangifer tarandus gonosomes [33]. The results of another detailed study concerning conserved syntenies of the X chromosome in four Cervidae species are currently available [31].
A substantial part of the X chromosome can be formed by heterochromatin which can be interspersed or present in blocks in the intercalary or centromeric/pericentromeric regions. The centromeric heterochromatin is mostly composed of satellite DNAs that are also present in centromeres of autosomes. Besides these, the repetitive DNA specific for sex chromosomes can exist, as described in Bovidae [18,24,34] and in Antilocapridae [19]. In Cervidae, an increased size of gonosomes caused by intercalary heterochromatin blocks has been observed in Rangifer tarandus and Elaphodus cephalophus [12,32].
In the present study, we report a comparative chromosome painting between cattle and various Cervidae species using chromosome painting probes prepared by flow cytometry and laser microdissection. Moreover, we used region specific and BAC probes for more detailed and accurate structural analysis of the autosomal rearrangements and the closer inspection of the structure of the cervid X chromosomes.
Material and methods
Ethical statement
All procedures performed in this study were in accordance with the ethical standards of the Veterinary Research Institute (Brno, Czech Republic), which complies with the Czech and European Union Legislation for the protection of animals used for scientific purposes. According to these regulations ethics approval was not required, as the biological material (blood/tissue) was obtained postmortem from animals upon animal slaughter in abattoir or which died during the hunting. The blood from living animals was collected by a ZOO veterinarian during other medical procedures. All collaborating ZOOs have license issued by the Ministry of the Environment of the Czech Republic (Act No 162/2003 Coll.).
Samples
Whole peripheral blood samples of the 3 specimens of red deer (Cervus elaphus, 2n = 68), milu deer (Elaphurus davidianus, 2n = 68), rusa deer (Cervus timorensis russa, 2n = 60), Eld's deer (Rucervus eldii, 2n = 58), moose (Alces alces, 2n = 68), reindeer (Rangifer tarandus, 2n = 70), roe deer (Capreolus capreolus, 2n = 70), fallow deer (Dama dama, 2n = 68), Chinese muntjac (Muntiacus reevesi, 2n = 46) and giraffe (Giraffa camelopardalis; 2n = 30) were obtained from captive born animals held in the Czech zoological gardens and a private deer farm (Bila Lhota, Czech Republic) (Table 1). Blood samples were cultured, harvested and fixed according to standard protocols as described previously [35]. Metaphase chromosome spreads for laser microdissection as well as slides for FISH analysis were prepared according to the procedures described previously [36]. For the assessment of deer karyotypes, GTG-banding was performed using the standard trypsin/Giemsa procedure [37].
Tissue cultures from cattle (Bos taurus, 2n = 60), which were used for flow sorting of chromosomes, were established in our laboratory according to [20].
Painting probes
Whole chromosomes or chromosomal regions for the construction of painting probes were isolated by flow sorting using MoFlo XDP Cell Sorter (Beckman Coulter, USA) [20] or microdissected by PALM Microlaser system (Carl Zeiss MicroImaging GmbH, Munich, Germany). The chromosomal DNA was then amplified by DOP-PCR (degenerate oligonucleotide primed polymerase chain reaction) [38]. Probe labelling was performed during the secondary PCR with Green-dUTP or Orange-dUTP (Abbott, IL, USA) [39].
BAC clones
BAC clones were selected from the CHORI-240 cattle library on the basis of NCBI Bos_taurus_UMD_3.1.1 Primary Assembly data and obtained from the Children’s Hospital Oakland Research Institute, BACPAC Resources, USA. Genomic BAC DNA was labelled with biotin-16-dUTP or digoxigenin-11-dUTP (Roche, Mannheim, Germany) using BioPrime Array CGH Genomic Labeling Module (Invitrogen, Carlsbad, CA, USA). Wherever it was possible, two to three BAC clones were selected from the same region and used together as a single probe to generate high quality hybridization signals in the phylogenetically more distant species. Table 2 lists the BAC clones used in the study.
BAC clones from the bovine CHORI-BAC library CH240 were selected for the detection of BTA1, BTA3 and BTAX chromosomes. Gray coloring marks two or three BAC clones, which were selected from the same region to generate more intense hybridization signal.
FISH
The analysis was performed following the protocols for FISH with painting or BAC probes. The posthybridization washing, image capture and processing were carried out as described in [36]. The BAC probes labelled with biotin-16-dUTP and digoxigenin-11-dUTP were detected with Avidin-CY3 (Amersham Pharmacia Biotech, Piscataway, NJ, USA) and antidigoxigenin-fluorescein (Roche). For the analysis, a Zeiss Axio Imager Z2 fluorescence microscope equipped with appropriate fluorescent filters was used. Images of well-spread metaphase chromosomes were captured and analysed using ISIS software (MetaSystems, Altlussheim, Germany).
Isolation of gonosome specific heterochromatin
R. tarandus chromosome X and Y were obtained by laser microdissection, amplified by DOP-PCR and subsequently cloned into a pDrive vector (Qiagen, Hilden, Germany). Species-specific clones were selected by dot-blot hybridization following the protocol of [34], fluorescently labeled and checked by FISH. Subsequently, plasmid DNA was isolated and sequenced. The resulting sequences were analyzed using BLASTN (Nucleotide collection (nr/nt) in Cetartiodactyla (taxid:91561)), BLAST2, RepeatMasker software and deposited in the GeneBank.
Results
Evolutionary rearrangements of Bovidae and Cervidae autosomes
To establish the chromosomal homologies between Bovidae and Cervidae families we hybridized painting probes derived from cattle (BTA) on metaphase chromosomes of nine species representing two main subfamilies of Cervidae (Cervinae and Capreolinae) (Fig 1 and Table 1). E. davidianus (EDA) was chosen as a model representative of family Cervidae whose karyotype has never been studied by molecular cytogenetic methods. Its karyotype is composed of 68 chromosomes, while only one autosomal pair is submetacentric. For the basic homology assessment, we hybridized whole chromosome painting probes (BTA1-BTA29) to chromosomal spreads of EDA. The painting probes from all 29 cattle autosomes revealed 35 conserved chromosomal segments on EDA genome. Only six probes for entire individual cattle chromosomes BTA1, 2, 5, 6, 8 and 9 each hybridized to two separate EDA chromosomes. Cattle chromosomes 26 and 28 are tandemly fused into one acrocentric chromosome, whereas BTA17 and BTA19 correspond to the arm of the sole submetacentric chromosome in Elaphurus genome. Each of all other cattle probes had homology only to one EDA chromosome. For detailed description of deer orthologs which represent parts of cattle chromosomes, distal region-specific probes were used. Distal region-specific probes also enabled us to distinguish between tandem and centric fusions in other studied species.
The numbers correspond to the bovine chromosome equivalents. The distances between species are not representative of the evolution time.
In another two Cervini species, C. elaphus (CEL) and D. dama (DDA), the hybridization of BTA1, 2, 5, 6, 8, 9, 17, 19, 26 and 28 painting probes on appropriate metaphase spreads revealed no obvious differences in their karyotypes in comparison with the EDA. Rusa deer (C. timorensis russa, CTR) and Eld's deer (R. eldii, REL) belong as well to the tribe Cervini and they share evolutionary chromosome changes as mentioned in other Cervinae (fusions of BTA17/19 and BTA26/28). However, their chromosome number is lowered due to the other centric fusions. Fusions of the chromosomes BTA2dist/7 and BTA5dist/8prox are shared in both species, whereas fusions BTA18/3, BTA5prox/22, and BTA18/1prox, BTA5prox/10, BTA22/1dist are exclusive to CTR and REL, respectively (Fig 1 and Table 1). G-banded karyotype of CTR with cattle chromosomal homologies is displayed in Fig 2.
Lines on the sides of CTR1-6 chromosome indicate boundaries between two cattle probes.
Three studied Capreolinae species differ in their karyotypes. The karyotype of C. capreolus (CCA) is composed of 68 acrocentric autosomes, whereas R. tarandus (RTA) has, besides acrocentric chromosomes, also one small submetacentric autosome pair and A. alces (AAL) has similar karyotype as RTA with another larger submetacentric pair. In all the above mentioned Capreolinae species we detected FISH signals of BTA1, 2, 5, 6, 8 and 9 on two separate chromosomes and tandem fusion of BTA28/26. The small submetacentric chromosome in RTA and AAL is homologous to the distal two-thirds part of BTA1, whereas the larger submetactric chromosome in AAL is a result of the centric fusion of BTA29/17. Examples of fluorescence in situ hybridization are displayed in Fig 3.
(A) Tandem fusion of BTA28/26 in milu deer (EDA) detected by region specific painting probes BTA26dist (red) and BTA28dist (green) (B) Centric fusion of BTA17 (red) and 19 (green) in rusa deer (CTR). (C) Centric fusion of BTA29/17 in moose (AAL) validated by probes BTA29dist (green) and BTA17dist (red). (D) Hybridization of BTA1 on reindeer (RTA) submetacentric and acrocentric orthologs. Centromeres are marked by lines.
The genus Muntiacus is famous for its high degree of interspecific karyotype diversity caused by numerous tandem fusions of ancestral cervid chromosomes. Using cattle painting probes we analyzed selected chromosomes of M. reevesi in which tandem fusions of cattle orthologs occurred during evolution. Hybridization results and order of cattle homologs on MRE1, 2, 3, 4, 5 and 11 are shown in Fig 4. Each of BTA 1, 5, 6, 8 and 9 probes always hybridized as two homologous blocks on different chromosomes in MRE genome. The BTA2 probe gave a fluorescence signal on MRE3 chromosome however, the signal was split into two blocks separated by BTA9 signal.
To find possible intrachromosomal rearrangements on autosomes we selected BAC probes derived from BTA1 and BTA3. Using twelve BAC clones for eleven different loci derived from BTA3 chromosome in EDA, CCA and MRE, we detected the same order of BAC probes as in the cattle. On the other hand, the ortholog of BTA1 in Cervidae is not only split into two separate chromosomes, but the larger of them is intrachromosomaly rearranged. BAC clones derived from the proximal parts of BTA1 (BAC1prox) hybridized to the proximal parts of the smaller ortholog, whereas BAC clones from distal (BAC1dist) and telomeric (BAC1tel) parts of BTA1 hybridized to the middle of the larger acrocentric ortholog (in EDA, DDA, CEL, CTR, REL, CCA and MRE) or to the distal (in relation to the centromere) part of the p-arm of submetacentric chromosome in AAL and RTA. Nonetheless, the position of BAC1dist and BAC1tel between acrocentric and submetacentric chromosomes is altered (Fig 5A).
(A) Hybridization of cattle probes BAC1qp (pink), BAC1qd (green), and BAC1qt (red) probes on chromosomes of cattle (BTA), rusa deer (CTR) and reindeer (RTA). (B) Schematic illustration demonstrating rearrangements of BTA1 orthologs in rusa deer (CTR), Chinese muntjac (MRE), reindeer (RTA), pronghorn (AAM), giraffe (GCA), okapi (OJO) and pygmy hippo (CLI). The dots approximate the positions of BAC1qp (pink), BAC1qd (green), and BAC1qt (red) probes.
Evolutionary rearrangements of the X chromosomes in Bovidae and Cervidae
For better depiction of X chromosome rearrangements, we used appropriate BAC clones derived from BTAX. The used BAC X clones were divided into four groups according to their location in cattle submetacentric X chromosome (BTAXp, BTAXqprox, BTAXqdist and BTAXPAR) (Table 2). In the studied Cervidae species, we observed morphologically two different types of X chromosomes, acrocentric type, which is present in Cervinae and a sub/metacentric one in Capreolinae. Using appropriate BAC clones, we also revealed two types of cervid X chromosome. The BTAXqdist and BTAXPAR regions always form one block and they are in the same orientation in acrocentric and submetacentric X variants and they are even conserved within Bovidae and Cervidae. In Capreolinae, the BTAXp and BTAXqprox regions form two blocks and they are separately inverted in comparison to cattle X. The acrocentric variant is much more complex, both BTAXp and BTAXqprox regions are rearranged (Fig 6). The BTAXp block is divided into two separate regions with a breakpoint site between BAC clones 67P21 and 442N4 (34–38 Mb on cattle X). Similarly, BTAXqprox is also separated in two regions with a breakpoint site between BACs 198N19 and 93K24 (55–57 Mb on cattle X).
The cattle BAC probes were divided into four groups marked with different colours. Positions of the BAC clones are on the side of the chromosome.
Isolation of gonosome specific heterochromatin
We prepared four different plasmid clones for heterochromatin regions present on sex chromosomes of RTA. Two of the obtained clones (6 and 56) gave hybridization signal on both X and Y chromosomes. Both clones gave the same overlapping signal on the Y chromosome, whereas at the end of q-arm of X chromosome clone 56 hybridized more proximally than clone 6. Clone 5 was localized in one large block under the centromeric region of the X chromosome. Clone 24 hybridized at the end of the q-arm of X chromosome and in the centromeric regions of most autosomes. Hybridization results are shown in Fig 7. No fluorescence signal was observed when these probes were hybridized on metaphase spreads of related species (ALL).
(X, Y—gonosomes; A-autosome).
The sequences of all clones were deposited in NCBI. No significant homologies between clone sequences were found using BLAST2 program. RepeatMasker revealed internal repeat only in clone 6. Three of the clones showed a similarity to different parts of the 135Mb long Ovis canadensis chromosome X sequence (CP011912) (Table 3).
Discussion
In this work, we provide the first genome-wide investigation of the karyotypes of various cervid species by cross-species FISH using cattle painting probes. Although the family Bovidae is in evolutionary terms a closely related group to the Cervidae, comparative evolutionary studies based on reciprocal chromosome painting are limited. Only Burkin et al. hybridized all sheep (Ovis aries) chromosome painting probes on chromosomes of Muntiacus muntjak vaginalis and Chi et al. applied whole chromosome probes from M. reevesi on metaphases of gayal (Bos frontalis) [40,41].
Our comparative chromosomal analysis confirmed that common karyotype differences between cattle and studied Cervidae species includes disruptions of orthologs of BTA1, 2, 5, 6, 8, 9 and tandem fusion BTA28/26. Whereas fissions of six cattle orthologs are characteristic evolutionary changes in Cervidae, the fusion BTA28/26 is typical for all pecoran species, excluding Bovidae. The same conclusions were deduced by zoo-FISH using dromedary probes (Camelus dromedarius) on Siberian roe deer (Capreolus pygargus) [10].
The main characteristic in all studied Cervini species is one common submetacentric chromosome which arose via Robertsonian fusion of segments homologous to the BTA17 and 19. This is in accordance with previous studies, based on G, R—banded chromosomes [3,7,25]. We confirmed centric fusions of cattle homologs in CTR and REL biarmed chromosomes, which were previously described by R-banding and BAC mapping [3].
From Capreolinae, we studied three species, where two of them (AAL and RTA) had a small metacentric chromosome in their karyotype. This submetacentric is not formed by centric fusion of small acrocentric ancestral chromosomes, but it originated from distal two-thirds of BTA1 homolog as in other cervid species, where it occurs in acrocentric form. In addition to this small submetacentric, another larger submetacentric chromosome is present in AAL genome. This chromosome arose from centric fusion of BTA29/17.
The tribe Muntiacini is well known for its high degree of interspecific karyotype diversity between closely related species, which is caused by numerous tandem fusions. For instance, M. muntjak vaginalis has the lowest chromosome number known in mammals (2n = 6/7), whereas M. reevesi (another morphologically similar species used in this study) has the highest diploid number (2n = 46) of all Muntiacus species [42]. More than one cattle painting probe hybridized on six pairs of MRE chromosomes. Number of synteny blocks on these six pairs of chromosomes ranged from two (MRE5, MRE11) to five (MRE2) (Fig 4). We confirmed that each of the orthologs of BTA1, 2, 5, 6, 8 and 9 always create two separate homologous blocks in MRE genome, although both BTA2 blocks are present on the MRE3 chromosome. Our obtained hybridization data are mostly in accordance with the results of Chi et al. who used probes from MRE on gayal (B. frontalis) metaphase chromosomes [41]. The only difference was found on MRE3 chromosome, where we distinguished two separated homologous blocks of BTA2, which cannot be revealed with the use of MRE probes. In MRE, the ancestral tandem fusion BTA28/26 is not present as a separate chromosome, but it forms one block on MRE2.
Evolutionary intrachromosomal rearrangements in ruminant autosomes are rare. However, the changes on homologs of BTA3 and BTA1 were detected in Antilopini, Giraffidae, Cervidae and Antilocapridae [19,24,25]. Therefore, we used BAC probes derived from BTA1 and BTA3 to detect possible intrachromosomal rearrangements in selected Cervidae species. Despite using high density BAC mapping on the BTA3 ortholog in studied Cervidae species, no changes in the order of BAC clones were revealed.
On the other hand, on the ortholog of BTA1 in all studied Cervidae species the complex rearrangements were detected. Those findings are not so surprising, because BTA1 orthologs are very diverse in many species due to the disruptions, inversions, and translocations. In all our nine Cervidae species, the BTA1 ortholog is split into two chromosomes of different sizes. The larger chromosome is acrocentric in most species with the exception of AAL and RTA where due to the centromere shift it has submetacentric form. The state of BTA1 orthologues in AAL and RTA represents simple disruption into two chromosomes with retained BAC order. The situation with BTA1 ortholog in EDA, DDA, CEL, CTR, REL, CCA and MRE is more complex. In these species, rearrangement of the larger ortholog representing two-thirds of BTA1 occurred (Fig 5A). This rearrangement of the BTA1 ortholog was previously described in Vietnamese sika deer (Cervus nippon) using a different set of BAC probes [25]. Another structure of the BTA1 ortholog is present in the genome of giraffe (Giraffa camelopardalis), pronghorn (Antilocapra americana) and okapi (Okapia johnstoni), where the BTA1 equivalent has been disrupted and inverted [19]. Outside the suborder Ruminantia, the BTA1 ortholog rearrangements were detected in pygmy hippo (Choeropsis liberiensis, CLI) from the family Hippopotamidae. The ortholog of BTA1 in pygmy hippo is split into two chromosomal arms (CLI6p, CLI11q) with other multiple rearrangements [20]. All mentioned BTA1 ortholog changes are summarized and illustrated in Fig 5B. Previous studies, which used dromedary and human painting probes, found that Ruminantia and Pecora ancestral characteristics of the BTA1 ortholog are retained in the giraffe and pronghorn chromosomes GCA2 and AAM1 [16,43]. From that ancestral state, one fission and one inversion occurred in the evolutionary line leading to Bovidae, Moschidae and Cervidae. In Cervidae, another fission of BTA1 ortholog formed two separate chromosomes (RTA and AAL) and subsequently in Cervinae and CCA another break and translocation occurred on the larger of them. Considering complexity of rearrangements of BTA1 ortholog in different species (Fig 5B), the employment of the high density BAC mapping could reveal progressive chromosomal changes which occurred during evolution. The increased rates of chromosomal rearrangements in some Pecoran chromosomes suggest the presence of the breakpoint hotspots that can explain such increased levels of rearrangements on these particular chromosomes [44,45].
Unlike autosomes, the sex chromosomes in Cetartiodactyla are variable between species due to the complex chromosomal changes including inversions, centromere shift and heterochromatic variation [22,30]. The X chromosome rearrangements have not been studied in more detail in most Cervidae species even though they proved to be a rich source of phylogenetic information [46]. For better depiction of X chromosome rearrangements, we used appropriate BAC clones (Table 2). Two different morphological types of X chromosomes were observed in studied Cervidae species, acrocentric type present in Cervinae and submetacentric one in Capreolinae. Both types of X differed not only by morphology, but also by the order of used BAC clones. The structure of the submetacentric variant is closer to the BTAX, whereas the acrocentric type is much more rearranged (Fig 6).
In Cervidae, the number of papers using Bovidae X chromosome BAC/cosmid clones is limited. Only a few cosmid clones were hybridized on the X chromosome of RTA and several BAC clones on two Cervidae species (C. nippon and Odocoileus virginianus), which are carriers of both X chromosome variants [27,33]. Recently, the structure of the X chromosome was described in four Cervidae species (Capreolus pygargus, AAL, DDA and EDA) by high resolution BAC mapping [31]. Our hybridization data are mostly in agreement with the above mentioned papers and substantially extend knowledge concerning the structure of cervid X chromosomes. In addition, our selection of X BAC clones enabled us to detected not only four, but five synteny blocks on the X chromosome of DDA that were not previously revealed by other BAC mapping [31]. We observed the same structure of X chromosome with five synteny blocks in all six studied Cervinae species (DDA, CEL, EDA, CTR, REL and MRE). On the other hand, Proskuryakova et al. studied only two Cervinae species (DDA and EDA) and they described two different types of X chromosome rearrangements in those two closely related animals. Studies that describe the detailed structure and rearrangements of the X chromosome outside the family Bovidae were scarce and limited only to a few species from the Giraffidae, Cervidae, Antilocapridae and Hippopotamidae [19,20,27]. Recently, new investigation of the X chromosome structure in eighteen different Cetartiodactyla species using comparative BAC mapping appeared [31]. Taking into account all known structures of X chromosomes in different species, the acrocentric (Cervinae) variant is unique, whereas the submetacentric (Capreolinae) type is similar to the X chromosome of the giraffe and hippo (without regard to the centromere position). For confirmation of the similarity between Capreolinae and giraffe X chromosome structure we hybridized the same set of BAC probes on giraffe chromosomes (Fig 6). The order of BAC clones on giraffe X was same as in Capreolus. With regard to those facts, the X chromosome structure present in Capreolinae, Giraffidae and Hippopotamidae appears to be the ancestral type. Those findings are in agreement with the most recently published data in which authors postulate that the Pecoran ancestral X chromosome (PAX) is preserved in the genomes of giraffe and moose and differs from the Cetartiodactyla ancestral X chromosome (CAX) only by centromere position [31].
Rearrangements on X chromosomes can be caused by the presence of heterochromatin blocks, which are composed of repetitive sequences. The occurrence of those sequences was described in Bovidae [18,24,34] and Antilocapridae [19]. Among Cervidae species, such heterochromatin blocks have only been observed in R. tarandus and Elaphodus cephalophus, which are responsible for the increased size of gonosomes [12,32]. In RTA, several heterochromatin motives were characterized by in situ hybridization of microdissected X chromosome regions [32]. In our study, we also used the microdissection method for generation of heterochromatin specific clones with the subsequent FISH localization and Sanger sequencing. Similarly to [32], we characterized several variants of heterochromatin repeats in RTA with different hybridization patterns. Although all examined clones were short (101–447 bp), they gave strong fluorescent signals confirming their repetitive character. Internal tandem repeat motif was not found in any of the studied clones, thus, the obtained sequence cannot be characterized as microsatellites. We confirmed that these repetitive sequences are characteristic for RTA as they do not hybridize to metaphase chromosomes of AAL.
Conclusion
In this study, we provide comparative chromosome painting between various Cervidae species and cattle. For the first time, we report in-depth examination of C. timorensis russa, E. davidianus, R. eldii, A. alces and R. tarandus karyotypes. We also carried out a detailed and accurate structural analysis of the cervid X chromosomes using bovine BAC clones. Heterochromatic blocks causing elongation of R. tarandus gonosomes were studied to the sequence level.
Acknowledgments
This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic under the project CEITEC 2020 (LQ1601) and by the Ministry of Agriculture of the Czech Republic (RO 0516). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors wish to thank Mr. Paul Veater (Bristol, UK) for proofreading the translated manuscript.
References
- 1.
Wilson DE, Reeder DM. Mammal Species of the World. Johns Hopkins University Press, Baltimore; 2005.
- 2. Fontana F, Rubini M. Chromosomal evolution in Cervidae. BioSystems. 1990;24: 157–174. pmid:2249009
- 3. Bonnet-Garnier A, Claro F, Thévenon S, Gautier M, Hayes H. Identification by R-banding and FISH of chromosome arms involved in Robertsonian translocations in several deer species. Chromosome Res. 2003;11: 649–663. pmid:14606627
- 4. Chi JX, Huang L, Nie W, Wang J, Su B, Yang F. Defining the orientation of the tandem fusions that occurred during the evolution of Indian muntjac chromosomes by BAC mapping. Chromosoma. 2005;114: 167–172. pmid:16010580
- 5. Yang F, O’Brien PC, Wienberg J, Ferguson-Smith MA. A reappraisal of the tandem fusion theory of karyotype evolution in Indian muntjac using chromosome painting. Chromosome Res. 1997;5: 109–117. pmid:9146914
- 6. Duarte JMB, Jorge W. Morphologic and cytogenetic description of the small red brocket (Mazama bororo Duarte,1996) in Brazil. Mammalia. 2009;67: 403–410.
- 7. Rubini M, Negri E, Fontana F. Standard karyotype and chromosomal evolution of the fallow deer (Dama dama L.). Cytobios. 1990;64: 155–161. pmid:2090390
- 8. Pinkel D, Landegent J, Collins C, Fuscoe J, Segraves R, Lucas J, et al. Fluorescence in situ hybridization with human chromosome-specific libraries: detection of trisomy 21 and translocations of chromosome 4. Proc Natl Acad Sci U S A. 1988;85: 9138–9142. pmid:2973607
- 9. Yang F, Fu B, O’Brien PCM, Nie W, Ryder OA, Ferguson-Smith MA. Refined genome-wide comparative map of the domestic horse, donkey and human based on cross-species chromosome painting: insight into the occasional fertility of mules. Chromosome Res. 2005;12: 65–76.
- 10. Dementyeva PV, Trifonov VA, Kulemzina AI, Graphodatsky AS. Reconstruction of the putative cervidae ancestral karyotype by chromosome painting of Siberian roe deer (Capreolus pygargus) with dromedary probes. Cytogenet Genome Res. 2010;128: 228–235. pmid:20413959
- 11. Yang F, O’Brien PC, Wienberg J, Neitzel H, Lin CC, Ferguson-Smith MA. Chromosomal evolution of the Chinese muntjac (Muntiacus reevesi). Chromosoma. 1997;106: 37–43. pmid:9169585
- 12. Huang L, Chi J, Nie W, Wang J, Yang F. Phylogenomics of several deer species revealed by comparative chromosome painting with Chinese muntjac paints. Genetica. 2006;127: 25–33. pmid:16850210
- 13. Huang L, Chi J, Wang J, Nie W, Su W, Yang F. High-density comparative BAC mapping in the black muntjac (Muntiacus crinifrons): Molecular cytogenetic dissection of the origin of MCR 1p+4 in the X1X2Y1Y2Y3 sex chromosome system. Genomics. 2006;87: 608–615. pmid:16443346
- 14. Yang F, Carter NP, Shi L, Ferguson-Smith MA. A comparative study of karyotypes of muntjacs by chromosome painting. Chromosoma. 1995;103: 642–652. pmid:7587587
- 15. Gallagher DS Jr, Womack JE. Chromosome conservation in the Bovidae. J Hered. 1992;83: 287–298. pmid:1401875
- 16. Kulemzina AI, Yang F, Trifonov VA, Ryder OA, Ferguson-Smith MA, Graphodatsky AS. Chromosome painting in Tragulidae facilitates the reconstruction of Ruminantia ancestral karyotype. Chromosome Res. 2011;19: 531–539. pmid:21445689
- 17. Ropiquet A, Hassanin A, Pagacova E, Gerbault-Seureau M, Cernohorska H, Kubickova S, et al. A paradox revealed: karyotype evolution in the four-horned antelope occurs by tandem fusion (Mammalia, Bovidae, Tetracerus quadricornis). Chromosome Res. 2010;18: 277–286. pmid:20204496
- 18. Rubes J, Kubickova S, Pagacova E, Cernohorska H, Di Berardino D, Antoninova M, et al. Phylogenomic study of spiral-horned antelope by cross-species chromosome painting. Chromosome Res. 2008;16: 935–947. pmid:18704723
- 19. Cernohorska H, Kubickova S, Kopecna O, Kulemzina AI, Perelman PL, Elder FFB, et al. Molecular cytogenetic insights to the phylogenetic affinities of the giraffe (Giraffa camelopardalis) and pronghorn (Antilocapra americana). Chromosome Res. 2013;21: 447–460. pmid:23896647
- 20. Frohlich J, Kubickova S, Cernohorska H, Musilová P, Muskova H, Rubes J. A comparative study of pygmy hippopotamus (Choeropsis liberiensis) karyotype by cross-species chromosome painting. J Mamm Evol. 2016;
- 21. Kulemzina AI, Trifonov VA, Perelman PL, Rubtsova NV, Volobuev V, Ferguson-Smith MA, et al. Cross-species chromosome painting in Cetartiodactyla: Reconstructing the karyotype evolution in key phylogenetic lineages. Chromosome Res. 2009;17: 419–436. pmid:19350402
- 22. Rubes J, Musilova P, Kopecna O, Kubickova S, Cernohorska H, Kulemsina AI. Comparative molecular cytogenetics in Cetartiodactyla. Cytogenet Genome Res. 2012;137: 194–207. pmid:22627059
- 23. Slate J, Van Stijn TC, Anderson RM, McEwan KM, Maqbool NJ, Mathias HC, et al. A deer (subfamily Cervinae) genetic linkage map and the evolution of ruminant genomes. Genetics. 2002;160: 1587–1597. pmid:11973312
- 24. Cernohorska H, Kubickova S, Kopecna O, Vozdova M, Matthee CA, Robinson TJ, et al. Nanger, Eudorcas, Gazella, and Antilope form a well-supported chromosomal clade within Antilopini (Bovidae, Cetartiodactyla). Chromosoma. 2014;124: 235–247. pmid:25416455
- 25. Bonnet A, Thévenon S, Claro F, Gautier M, Hayes H. Cytogenetic comparison between Vietnamese sika deer and cattle: R-banded karyotypes and FISH mapping. Chromosome Res. 2001;9: 673–687. pmid:11778690
- 26. Lorenzi LD, Planas J, Rossi E, Malagutti L, Parma P. New cryptic karyotypic differences between cattle (Bos taurus) and goat (Capra hircus). Chromosome Res. 2015;23: 225–235. pmid:25612562
- 27. Gallagher DS, Davis SK, De Donato M, Burzlaff JD, Womack JE, Taylor JF, et al. A molecular cytogenetic analysis of the tribe Bovini (Artiodactyla: Bovidae: Bovinae) with an emphasis on sex chromosome morphology and NOR distribution. Chromosome Res. 1999;7: 481–492. pmid:10560971
- 28. Iannuzzi L, King WA, Di Berardino D. Chromosome evolution in domestic bovids as revealed by chromosome banding and FISH-mapping techniques. Cytogenet Genome Res. 2009;126: 49–62. pmid:20016156
- 29. Piumi F, Schibler L, Vaiman D, Oustry A, Cribiu EP. Comparative cytogenetic mapping reveals chromosome rearrangements between the X chromosomes of two closely related mammalian species (cattle and goats). Cytogenet Cell Genet. 1998;81: 36–41. pmid:9691172
- 30. Robinson TJ, Harrison WR, Ponce de León FA, Davis SK, Elder FF. A molecular cytogenetic analysis of X chromosome repatterning in the Bovidae: transpositions, inversions, and phylogenetic inference. Cytogenet Cell Genet. 1998;80: 179–184. pmid:9678354
- 31. Proskuryakova AA, Kulemzina AI, Perelman PL, Makunin AI, Larkin DM, Farré M, et al. X Chromosome evolution in Cetartiodactyla. Genes. 2017;8: 216. pmid:28858207
- 32. Lee C, Griffin DK, O’Brien PC, Yang F, Lin CC, Ferguson-Smith MA. Defining the anatomy of the Rangifer tarandus sex chromosomes. Chromosoma. 1998;107: 61–69. pmid:9567201
- 33. Prakash B, Kuosku V, Olsaker I, Gustavsson I, Chowdhary BP. Comparative FISH mapping of bovine cosmids to reindeer chromosomes demonstrates conservation of the X-chromosome. Chromosome Res. 1996;4: 214–217. pmid:8793206
- 34. Pauciullo A, Kubickova S, Cernohorska H, Petrova K, Berardino D, Ramunno L, et al. Isolation and physical localization of new chromosome-specific centromeric repeats in farm animals. VetMed 2006;51:0–9
- 35. Cernohorska H, Kubickova S, Vahala J, Robinson TJ, Rubes J. Cytotypes of Kirk’s dik-dik (Madoqua kirkii,Bovidae) show multiple tandem fusions. Cytogenet Genome Res. 2011;132: 255–263. pmid:21124018
- 36. Cernohorska H, Kubickova S, Vahala J, Rubes J. Molecular insights into X;BTA5 chromosome rearrangements in the tribe Antilopini (Bovidae). Cytogenet Genome Res. 2012;136: 188–198. pmid:22327909
- 37. Seabright M. A rapid banding technique for human chromosomes. Lancet. 1971;2: 971–972.
- 38. Telenius H, Carter NP, Bebb CE, Nordenskjold M, Ponder BAJ, Tunnacliffe A. Degenerate oligonucleotide-primed PCR: General amplification of target DNA by a single degenerate primer. Genomics. 1992;13: 718–725. pmid:1639399
- 39. Kubickova S, Cernohorska H, Musilova P, Rubes J. The use of laser microdissection for the preparation of chromosome-specific painting probes in farm animals. Chromosome Res. 2002;10: 571–577. pmid:12498346
- 40. Burkin DJ, Yang F, Broad TE, Wienberg J, Hill DF, Ferguson-Smith MA. Use of the Indian muntjac idiogram to align conserved chromosomal segments in sheep and human genomes by chromosome painting. Genomics. 1997;46: 143–147. pmid:9403070
- 41. Chi J, Fu B, Nie W, Wang J, Graphodatsky AS, Yang F. New insights into the karyotypic relationships of Chinese muntjac (Muntiacus reevesi), forest musk deer (Moschus berezovskii) and gayal (Bos frontalis). CGR. 2005;108: 310–316. pmid:15627750
- 42. Wurster DH, Benirschke K. Indian muntjac, Muntiacus muntjak: a deer with a low diploid chromosome number. Science. 1970;168: 1364–1366. pmid:5444269
- 43. Kulemzina AI, Perelman PL, Grafodatskaya DA, Nguyen TT, Thompson M, Roelke-Parker ME, et al. Comparative chromosome painting of pronghorn (Antilocapra americana) and saola (Pseudoryx nghetinhensis) karyotypes with human and dromedary camel probes. BMC Genetics. 2014;15: 68. pmid:24923361
- 44. Bailey JA, Baertsch R, Kent WJ, Haussler D, Eichler EE. Hotspots of mammalian chromosomal evolution. Genome Biology. 2004;5: R23. pmid:15059256
- 45. Pevzner P, Tesler G. Human and mouse genomic sequences reveal extensive breakpoint reuse in mammalian evolution. Proc Natl Acad Sci USA. 2003;100: 7672–7677. pmid:12810957
- 46. Robinson TJ, Harrison WR, Ponce de León A, Elder FF. X chromosome evolution in the suni and eland antelope: detection of homologous regions by fluorescence in situ hybridization and G-banding. Cytogenet Cell Genet. 1997;77: 218–222. pmid:9284920