Next Article in Journal
Intestinal Farnesoid X Receptor Modulates Duodenal Surface Area but Does Not Control Glucose Absorption in Mice
Next Article in Special Issue
Analysis of CircRNA Expression in Peripheral Blood of Holstein Cows in Response to Heat Stress
Previous Article in Journal
miR-103-3p Regulates the Differentiation and Autophagy of Myoblasts by Targeting MAP4
Previous Article in Special Issue
Convergent Genomic Signatures of Cashmere Traits: Evidence for Natural and Artificial Selection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 4
10.3390/ijms24044134
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Satellite-Free Centromere in Equus przewalskii Chromosome 10

by
Francesca M. Piras
1,†,
Eleonora Cappelletti
1,†,
Wasma A. Abdelgadir
1,
Giulio Salamon
2,
Simone Vignati
3,
Marco Santagostino
1,
Lorenzo Sola
1,
Solomon G. Nergadze
1 and
Elena Giulotto
1,*
1
Department of Biology and Biotechnology “Lazzaro Spallanzani”, University of Pavia, 27100 Pavia, Italy
2
Oasi di Sant’Alessio, Sant’Alessio con Vialone, 27016 Pavia, Italy
3
Independent Researcher, 27100 Pavia, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the work.
Int. J. Mol. Sci. 2023, 24(4), 4134; https://doi.org/10.3390/ijms24044134
Submission received: 20 January 2023 / Revised: 15 February 2023 / Accepted: 16 February 2023 / Published: 18 February 2023
(This article belongs to the Special Issue Animal Genomes and Epigenomes)
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

The centromere is the site of kinetochore assembly required for chromosome segregation during cell division. The centromeric function is epigenetically specified by the histone H3 variant CENP-A, which is highly conserved during evolution [1]. On the other hand, the underlying DNA sequences are extremely variable and are neither necessary nor sufficient to determine the centromeric function [2,3,4,5]. In mammals, centromeres are typically characterized by extended arrays of tandemly repeated sequences, called satellite DNA. These sequences represent the most rapidly evolving DNA sequences in eukaryotic genomes [6] and their presence has so far hampered comprehensive molecular analysis of these intriguing loci. In this context, our discovery that the species of the genus Equus are characterized by an extraordinarily high number of centromeres completely devoid of satellite DNA made these species a good model system for studying the epigenetic control and evolution of the centromere function [7,8,9,10,11,12,13,14,15,16,17].
The exceptional genome plasticity of equid species is reflected in their rapid karyotypic evolution and is marked by numerous events of centromere repositioning and chromosomal rearrangements [18,19,20]. Centromere repositioning, which is the shift of the centromere position on the chromosome without sequence rearrangement [21,22], gave rise to evolutionarily new centromeres [11,18,19,20,23]. A first cytogenetic approach revealed that this phenomenon played a key role during the evolution of equids [11,16,23]. Horse satellite repeats were cloned and analyzed by Southern blotting and FISH in early studies by different groups [24,25,26,27]. We then studied the chromosomal distribution of satellite sequences, namely 37cen and 2PI, in E. caballus (domestic horse), E. asinus (donkey), E. grevyi (Grevy’s zebra), and E. burchelli (Burchell’s zebra) which revealed a peculiar uncoupling between centromeric function and satellite repeats. In particular, several centromeres were found to be devoid of satellite DNA whereas blocks of satellite DNA were observed at some chromosomal termini, representing relics of ancestral inactivated centromeres [13,16].
Later on, using an anti-CENP-A antibody in ChIP experiments on chromatin extracted from primary fibroblast cell lines, we proved at the molecular level that several equid centromeres are completely devoid of satellite repeats [7,9,12,13,14,17]. While only one satellite-free centromere was found in the horse, on chromosome 11 (ECA11) [12], 16 were described in the donkey [7], 15 in the Burchell’s zebra, and 13 in the Grevy’s zebra [14]. Interestingly, it was observed that these CENP-A binding domains can slide within a genomic region of a few hundred kbs, probably limited by epigenetic boundaries. These domains were defined as epialleles and the phenomenon was called centromere sliding [7,9,13,14]. Studies on families composed of donkeys, horses, and their hybrid offspring (mule/hinny) revealed that the epialleles are inherited as Mendelian traits but their position can slide in one generation [7] while it is stable during propagation in cell culture, suggesting that the sliding can presumably take place during meiosis [7].
Among equid species, horses are considered the closest extant species to the equid ancestor, whose karyotype was mainly characterized by acrocentric chromosomes [19]. While the centromeres of the domestic horse were deeply studied both at the cytogenetic and at the molecular level [7,8,9,12,15,16,28,29], the centromeres of the other extant horse species, the Przewalski’s horse, are still lacking a detailed characterization.
Przewalski’s horses have never been successfully domesticated. It was suggested that the Przewalski’s horse species is a direct ancestor of the domestic horse [30,31]. However, this hypothesis is controversial with data indicating that Przewalski’s and domestic horses may represent two lineages that diverged between 38,000 and 72,000 years ago [32,33,34,35]. Recently, Gaunitz and colleagues demonstrated that Equus przewalskii is the feral descendant of horses first herded in the Central Asian steppes about 5500 years ago and their feralization likely involved multiple biological changes such as pigmentation as well as morpho-anatomical characteristics [36,37]. After the discovery of this species in the 1870s in the Asian steppes, Przewalski’s horses become extinct in the wild by 1969 but survived in captivity. Since the 1970s, horses were exchanged between zoos in order to reduce the level of inbreeding [38,39,40].
In the 1990s, reintroduction projects brought the species back in the wild but the current population is still endangered and entirely descending from only 13 founding individuals [32,41]. At present, the worldwide population of the living Przewalski’s horses has reached nearly 2500 individuals [42]. Among them, approximately 1360 horses live in the wild in China and Mongolia, 900 are distributed in zoos in Europe, and 120 in wildlife parks in the US [39,40,42,43,44].
Sequence analysis of the control regions from mitochondrial DNA allowed the identification of only four Przewalski’s horse matrilines [41,45]. In 2011, four complete mitogenomes, representing all four surviving mitochondrial lineages, were sequenced and assigned to haplogroups I, II, and III [32,46]. Later on, Achilli and colleagues grouped haplogroups I and II, which were similar, into haplogroup F, which is specific of Przewalski’s horse [47], while haplogroup III resulted closely related to the J and K haplogroups of domestic horses. This observation suggests that there was gene flow among ancestors of the present-day populations or common ancestry [31,47].
The karyotypes of Przewalski’s and domestic horses are very similar with 66 and 64 chromosomes, respectively. The two acrocentric chromosomes, 23 and 24, possibly representing the ancestral configuration, underwent a centric fusion to generate the metacentric chromosome 5 of the domestic horse [48,49,50]. Despite their different chromosomal numbers, the domestic and the Przewalski’s horses can produce viable and fertile offspring [51] and this is the reason why hybrids may be present in Przewalski’s populations both in the wild and in captivity.
Here, we investigated the chromosomal distribution of satellite DNA families in Equus przewalskii (EPR). Moreover, we demonstrated at the cytogenetic and the molecular level that the centromere of EPR10 is devoid of satellite DNA sequences. This centromere is located in the genomic region orthologous to the satellite-free centromeric region on ECA11.

2. Results

2.1. Characterization of a Fibroblasts Cell Line from Equus przewalskii

We established a primary fibroblast cell line from a skin biopsy of a 17-year-old male Przewalski’s horse. To confirm that this animal was a Przewalski’s horse (EPR) and not a hybrid with the domestic horse, we first analyzed its karyotype using G-banding obtained by DAPI staining (Figure 1). As expected for this species, the number of chromosomes was 66 with 12 metacentric pairs, 20 acrocentric pairs, including chromosomes 23 and 24 which are typical of this species, and the XY sex chromosomes. To better identify the EPR chromosomes 23 and 24, which are orthologous to ECA5p and ECA5q [49,50], we performed FISH experiments using the BAC probes CHORI241-17G8 and CHORI241-25H7 specific for ECA5p and ECA5q, respectively (Figure 2). The probes were used to hybridize metaphase spreads from a domestic horse (panel A), from a E. caballus/E. przewalskii hybrid (2n = 65) (panel B) and from the E. przewalskii cell line (panel C).
In the domestic horse, following hybridization with the CH241-17G8 probe, FISH signals were detected at the proximal region of the p-arm of chromosome 5 (Figure 2A left). As expected, in the hybrid individual, FISH signals were detected on the p-arm of the metacentric chromosome 5, derived from the domestic horse parent, and on the acrocentric Przewalski’s horse chromosome 23 (Figure 2B left). Finally, as expected for a Przewalski’s horse, in our individual, the BAC clone hybridized with two acrocentric chromosomes, which compose the pair of homologous chromosomes EPR23 (Figure 2C left). The BAC clone CH241-25H7 was previously mapped by FISH on the q-terminus of chromosome ECA5 [23]. This localization was confirmed in the present work (Figure 2A right). In the hybrid individual, FISH signals were observed on the metacentric horse-derived chromosome 5 and on the acrocentric chromosome 24 of Przewalski’s horse (Figure 2B right). As expected for a Przewalski’s horse, the BAC clone CH241-25H7 was located on the acrocentric EPR24 chromosome pair (Figure 2C right).
To further prove that the fibroblast cell line derives from a Przewalski’s horse, we assembled the sequence of its mitochondrial DNA (access number OQ285861) from whole-genome Illumina reads that were used as input in the ChIP-seq experiment described below. The sequence was then aligned with 104 mitogenomes available in GenBank from domestic horses, Przewalski’s horses, and their hybrids [47,52]. The mitogenome of our individual belongs to haplogroup F, which is typical of Przewalski’s horses.

2.2. Chromosomal Distribution of Satellite DNA Families in Equus przewalskii

To analyze the localization of the two major horse satellite DNA families in Przewalski’s horse, we carried out FISH experiments using the previously described plasmids containing 37cen and 2PI satellite sequences as probes [16,53]. Representative metaphase spreads hybridized with the two probes are shown in Figure 3 and the chromosomal distribution of the two satellite sequences is reported in Table 1.
Hybridization signals of the 37cen satellite are present at the primary constriction of all chromosomes with the exception of EPR10 (Table 1). Signal intensities are heterogeneous suggesting that the number of repeats is variable among different chromosomes (Figure 3A). The primary constrictions of chromosomes 2, 6, 18, and 29 show 37cen signals only on a subset of the analyzed chromosomes (Table 1), suggesting that the number of repeats at these loci is low.
Similarly to 37cen, the 2PI satellite was detected at the primary constriction of all chromosomes except EPR10 (Table 1) and the intensity of the signals was heterogeneous among different chromosomes, suggesting inter-chromosomal variability in the number of repeats (Figure 3B). The centromeres of chromosomes 1, 4, 11, 29, X, and Y were labeled only in a subset of the analyzed chromosomes (Table 1), suggesting a low number of repeats at these loci.
To test the presence of satellite DNA families other than 37cen and 2PI, we performed a FISH analysis using total genomic DNA as a probe. We previously used this procedure to identify regions containing highly repetitive sequences because of their different hybridization kinetics compared to single-copy sequences [16]. Even with this probe, chromosome 10 is lacking any hybridization signal while all other chromosomes are labeled at their primary constriction. These observations suggest that the centromere of chromosome 10 is the only one devoid of satellite DNA. EPRX, similarly to what was previously observed for the domestic horse [16], showed two hybridization signals, one at the primary constriction and one at an interstitial position of the q arm (Figure 3C and Table 1). In Figure 3D, an example of double-color FISH with 2PI satellite and total genomic DNA as probes is shown, confirming that chromosome 10 is the only one lacking any signal. The Przewalski’s specific acrocentric chromosomes 23 and 24, which are orthologous to horse chromosome 5p and 5q, respectively [49], are labeled at their centromeric end with the three probes.

2.3. ChIP-Seq Analysis of Centromeric Domains

To test, at the sequence level, whether the centromere EPR10 is satellite-free, a ChIP-seq experiment with an anti-CENP-A antibody was carried out on the chromatin extracted from the primary skin fibroblasts described above.
Since a chromosomal-level genome assembly for Przewalski’s horse is not available, given the high karyotypic and sequence identity with the domestic horse [35,49,50,54], we mapped the ChIP-seq reads on the horse reference genome. We used the EquCab2.0 reference where the satellite-free centromeric region of horse chromosome 11 is well assembled [7,12]. Using the pipeline that we developed for the identification of satellite-free centromeres in other equids [7,14], a CENP-A binding domain was identified in the region comprised between 27.59 and 27.79 Mb of chromosome 11 in the reference genome (Figure 4A). In the Przewalski’s horse, this region is located on chromosome 10, which is orthologous to horse chromosome 11 (Figure 4A). The regular shape of this peak indicates that the underlying genomic region is highly conserved between Przewalski’s and domestic horse.
We previously demonstrated that, in equids, homologous chromosomes can carry CENP-A binding domains in different positions, which correspond to different epialleles [7,8,9]. To analyze the organization of epialleles on EPR10 of the animal studied here, we used a single nucleotide variant (SNV)-based approach.
We identified heterozygous nucleotide positions within the centromeric domains using the ChIP reads. These positions are listed in Table 2 and marked as bars in Figure 4A. We detected eight heterozygous nucleotide positions which are distributed along the entire length of the peak, suggesting that the CENP-A binding domains on the two homologs are overlapping (Figure 4A, Table 2).
Since satellite-based centromeres contain megabases of unassembled satellite arrays, we could not obtain ChIP-seq peaks from these centromeres which are probably organized likewise the typical mammalian centromeres, as already shown for the domestic horse [28].
To define the satellite DNA family bearing the centromeric function, we aligned the ChIP-seq reads from CENP-A immunoprecipitated and from input chromatin with the consensus sequences of 37cen and 2PI. To quantify the enrichment of these sequences in CENP-A bound chromatin, we calculated the ratio between normalized read counts in immunoprecipitated and in input DNA (Figure 4B). As control, we used the ERE-1 retrotransposon, which is interspersed throughout the equid genomes and is not expected to be involved in the centromeric function [28]. A 3-fold enrichment was observed for the 37cen satellite while 2PI and ERE-1 were present at similar levels in the two fractions (Figure 4B). These results demonstrated that, as for the domestic horse [28], 37cen is the main functional centromeric satellite sequence.

3. Discussion

Since Przewalski’s and domestic horses can breed and give birth to fertile offspring, the presence of hybrid individuals, both in the wild and in captivity, represents one of the major problems in the conservation programs of Przewalski’s horses. Therefore, following the establishment of a new fibroblast cell line, we first tested whether the donor was a Przewalski’s horse not deriving from a recent hybridization with the domestic horse. Indeed, the karyotype of this individual showed 66 chromosomes comprising the acrocentric chromosome pairs 23 and 24, which are typical of this species, whereas the metacentric chromosome 5, which is typical of the domestic horse, was absent. Furthermore, sequence analysis of its mitogenome revealed that this individual belongs to haplogroup F, which is typical of the Przewalski’s species [47].
In previous work from our laboratory, the chromosomal distribution of the two major equid satellite DNA families, 37cen and 2PI, was investigated by FISH in E. caballus, E. asinus, E. grevyi, and E. burchelli [16]. This analysis demonstrated that, in the genus Equus, several centromeres are devoid of satellite DNA and that, in agreement with the library hypothesis for satellite DNA evolution [55,56], horses, donkeys, and zebras share a common set of satellite DNA families which underwent expansion or shrinkage in the different species [16]. In the horse, either one or both these satellites are present on all primary constrictions, except the one of chromosome 11, which contains the unique satellite-free centromere of this species [12,16].
Here, we showed that there is a good degree of conservation of satellite DNA in Przewalski’s and in domestic horse. However, some differences in the distribution of satellite DNA families between the two species were observed. For example, 37cen signals were detected at the primary constriction of EPR2 while the ortholog of the domestic horse (ECA2) lacks 37cen signals [16]. Moreover, the primary constrictions of chromosomes 1, 4, 11, and X of E. przewalskii are labeled by 2PI signals, while, in the orthologous chromosomes of the domestic horse (1, 4, 12, and X), such signals are undetectable [16]. These observations suggest that, in E. caballus, the above-mentioned centromeres may lack these satellite DNA families or carry a small number of repeated units which are undetectable at the FISH resolution level.
The high similarity between Przewalski’s and domestic horse was also confirmed following hybridization of metaphase chromosomes with total genomic DNA. This procedure was used to identify regions containing very abundant tandem repeats due to the different hybridization kinetics of highly reiterated sequences versus single copy DNA [16]. The results of the FISH experiments with the three probes (37cen, 2PI, and genomic DNA) suggest that the number of satellite repeats is variable at different primary constrictions (Figure 3 and Table 1). In addition, the appearance of a non-centromeric hybridization signal on the q arm of chromosome X, visible with the genomic DNA probe only, indicates that satellite repeats other than 37cen and 2PI are located at this site. We previously showed that, also in the domestic horse, the q arm of chromosome X contains satellite repeats other than 37cen and 2PI [16].
Interestingly, both major satellite families, 37cen and 2PI, were detected at the centromeres of the two acrocentric chromosomes EPR23 and EPR24 while the 37cen only is present at the primary constriction of the metacentric domestic horse ortholog (ECA5) [16]. Therefore, at the centromere of ECA5 the 2PI satellite is either absent or undetectable by FISH suggesting that this satellite may have been lost during fusion of two ancestral acrocentric chromosomes corresponding to EPR23 and EPR24. The generation of ECA5, in the domestic horse lineage, from a centric fusion was previously proposed [57] and the loss of satellite repeats as a consequence of Robertsonian fusion was recently demonstrated by us in the zebras [14]. However, since the ECA5 synteny is maintained in all Ceratomorpha and disrupted in all equids with the exception of the domestic horse, the phylogenetic history of this chromosome remains controversial. In particular, an ECA5 ortholog is acrocentric in rhinos, whose karyotypes remained quite stable during evolution and similar to its orthologous element in the hypothetical perissodactyl ancestral karyotype [18,19]. This situation suggests an alternative model: a fission event generated the two acrocentric chromosomes corresponding to EPR23 and EPR24. According to this hypothesis, the submetacentric chromosome of the domestic horse may have lost the 2PI satellite during evolution.
The FISH experiments demonstrated that EPR10, the ortholog of ECA11, is the only Przewalski’s chromosome devoid of satellite DNA signals at its primary constriction. Evidence of a possible lack of satellite repeats on Przewalski’s chromosome 10 was provided in an early pioneering study from Ryder and Hansen in 1979 [48] where a GC rich satellite band was purified by density centrifugation. Radioactively labeled cRNA was then hybridized to metaphase spreads from E. przewalskii. ChIP-seq results confirmed at the sequence level that the centromere of EPR10 is completely devoid of satellite DNA, similarly to what is observed for its ECA11 ortholog in the domestic horse. In agreement with our cytogenetic data, no other satellite-free centromere was detected in Przewalski’s horse.
The CENP-A enrichment peak detected on EPR10 lays in the same genomic region in which horse CENP-A binding domains are known to slide [7,9,12]. This result suggests that the event of centromere repositioning which gave rise to this centromere occurred in the common ancestor, before the separation of the two horse lineages, between 38,000 and 72,000 years ago [34]. In several domestic horses, two well separated CENP-A binding domains, corresponding to different epialleles on the two homologs, were detected [8,9]. In the Przewalski’s individual analyzed here, we observed a single Gaussian-like peak indicating that this individual is homozygous for centromere position, in agreement with the high level of inbreeding of this species [35,36,41,58,59,60,61].
We could not localize the CENP-A binding domains from the satellite-based centromeres since they contain megabases of satellite arrays which are not assembled in the reference genomes. However, as in the domestic horse [28] the 37cen satellite family is enriched in the CENP-A immunoprecipitated chromatin, thus bearing the centromeric function. The 2PI satellite is not enriched in the CENP-A bound chromatin and, as in the domestic horse, is preferentially distributed at pericentromeric regions.
The high degree of similarity of satellite DNA in terms of composition and distribution between Przewalski’s and domestic horses is in agreement with the recent divergence of the two lineages [32,36]. Interestingly, the localization of these satellite sequences in Burchell’s and Grevy’s zebras varies in the two lineages in a species-specific way [16]. This difference may be related to the fact that the radiation between the two horses is more recent if compared to the separation of the two zebras (>1MYA) [62].
In conclusion, the results of this work confirm that Przewalski’s and domestic horses are phylogenetically very close, supporting the hypothesis that their radiation from a common ancestor occurred in recent evolutionarily times.

4. Materials and Methods

4.1. Cell Lines

A primary fibroblast cell line was established from a skin biopsy from a male Przewalski’s horse living in captivity in the “Oasi di Sant’Alessio” zoological park. This animal was obtained from the Berlin Zoo (CITES number: De-BLN04031635 16.03.2004 Germany). The biopsy was taken post-mortem after the horse underwent euthanasia following an incurable fracture. The cell line was obtained from small fragments of dermis using a protocol well-established in our laboratory [7,8,9,11,12,15,16,23,28,29,63,64,65,66,67,68]. The horse primary fibroblast cell line was obtained from the skin of a slaughtered animal. The animal was being processed as part of the normal work of the abattoirs.
The cells were cultured in high-glucose DMEM medium (Euroclone, Milan, Italy) supplemented with 20% fetal bovine serum (Euroclone, Milan, Italy), 2 mM L-glutamine (Merck Life Science, Milan, Italy), 1% penicillin/streptomycin (Merck Life Science, Milan, Italy), and 2% non-essential amino acids (Euroclone, Milan, Italy). The cells were maintained at 37 °C in a humidified atmosphere of 5% CO2. Using this procedure, relatively pure cultures of skin fibroblasts are obtained. The fibroblasts show a typical elongated shape and grow aligned and in bundles when confluent. As shown also by other groups [69], the above-mentioned medium supports the growth of fibroblasts whereas other cell populations such as keratinocytes need additional supplements and growth factors or are mitotically less active.
Primary fibroblasts from a privately owned hybrid female individual derived from a cross between a domestic horse and a Przewalski’s horse were previously collected in connection with another project and kindly provided by Ernest Bailey from the University of Kentucky.

4.2. Metaphase Spread Preparation and FISH

For metaphase spread preparation, mitoses were mechanically collected by blowing the medium on the dish surface. Chromosome preparations were performed with the standard air-drying procedure and were used for karyotype analysis or fluorescence in situ hybridization (FISH).
Whole genomic DNA from Przewalski’s horse fibroblasts was extracted according to standard procedures. BAC clones (CHOR1241-17G8, coordinates in EquCab2.0: chr5:37,011,456–37,167,308; CHOR1241-25H7 coordinates in EquCab2.0: chr5: 98,031,303–98,224,666) and lambda phage 37cen and 2PI clones [16,53] were extracted from 10 mL bacterial cultures with the Quantum Prep Plasmid miniprep kit (BioRad, Milan, Italy), according to supplier instructions. All probes were labeled by nick translation with Cy3-dUTP (Enzo Life Sciences, Milan, Italy) or Alexa488-dUTP (Life Technologies, Monza, Italy) as previously described [64]. FISH was performed as previously described [68].
Chromosomes were counterstained by DAPI. Digital grayscale images for fluorescence signals were acquired with a fluorescence microscope (Zeiss Axio Scope.A1, Zeiss, Göttingen, Germany) equipped with a cooled CCD camera (Teledyne Photometrics, Birmingham, UK). Pseudocoloring and merging of images were performed using the IPLab 3.5.5 Imaging Software (Scanalytics Inc., Fairfax, VA, USA). Chromosomes were identified by DAPI banding according to the published karyotypes [49,50].

4.3. ChIP-Seq

Chromatin from primary fibroblasts was cross-linked with 1% formaldehyde, extracted, and sonicated to obtain DNA fragments ranging from 200 to 800 bp. Immunoprecipitation was performed as previously described [7] using a serum against horse CENP-A protein [8]. Paired-end sequencing was performed through an Illumina NovaSeq 6000 platform by IGA Technology Services (Udine, Italy).

4.4. Bioinformatic Analysis of ChIP-Seq Data

The identification of satellite-free centromeres was performed as previously described [14]. Reads were aligned with paired-end mode to the EquCab2.0 horse reference genomes with Bowtie2 aligner (version 2.4.2) using default parameters [70,71]. Normalization of read coverage of the ChIP datasets against the input datasets was performed using bamCompare (deepTools suite, 3.5.0 version) [72] using RPKM (Reads Per Kilobase per Million mapped reads) normalization in subtractive mode. The resulting coverage files were visualized using Integrative Genomics Viewer (IGV) software (2.9.2 version). Peaks were obtained with pyGenomeTracks (3.6 version) [73].
To quantify the number of reads corresponding to each repetitive element, the ChIP-seq reads from CENP-A immunoprecipitated chromatin and from input were aligned with the consensus sequences of 37cen (“SAT_EC”on Repbase), 2PI (“ES22”on Repbase and “AH010654.2” on NCBI Nucleotide), and ERE-1 retrotransposon (“ERE1”on Repbase) [74] using Bowtie2 (2.4.2 version) with single-end mode. RPKM normalized read counts were obtained using SAMtools (1.11 version) [75]. To measure enrichment due to immunoprecipitation with CENP-A, the ratio between normalized read counts in the immunoprecipitated and input samples was calculated.

4.5. SNV Analysis

To identify heterozygous single nucleotide positions in the CENP-A binding domain, we first removed duplicates from the ChIP BAM files using SAMtools [75,76]. SNVs in the centromeric region were identified using BCFtools (1.11 version) [76] and filtered for quality higher than 20 to exclude low quality calls, as suggested by the program. The minimum number of reads per position required for the analysis was 4, an arbitrary threshold that we previously used to identify SNVs in E. asinus [7]. The required frequency of the least frequent variant was equal or more than 0.3 [7]. We used only SNVs present within single copy sequences. As previously described [7], the presence of heterozygous nucleotide positions indicates the region of overlapping between the CENP-A binding domains of the two homologs.

4.6. Mitogenomic Analysis

The mitochondrial DNA sequence of our Przewalski’s horse (accession number OQ285861) was obtained from the Input reads aligned to chrM of the EquCab2.0 reference genome. The obtained sequence was aligned with mitochondrial genomes retrieved in GeneBank from 82 domestic horses (JN398377-JN398457, HQ439467, KT757764), 20 Przewalski’s horses (KT221844, KT221845, KT368742-KT368756, KT757761, HQ439484, JN398402, JN398403) and two Przewalski’s horse/domestic horse hybrids (KT368757 and KT368758). Sequence alignment and phylogenetic analysis was performed using MEGAX [77] as previously described by Kvist and collaborators [52].

Author Contributions

Conceptualization, E.G.; data curation, F.M.P., E.C., W.A.A. and S.G.N.; investigation, F.M.P., E.C., W.A.A., M.S., L.S., S.G.N. and E.G.; resources, G.S. and S.V.; writing—original draft preparation, F.M.P. and E.C.; writing—review and editing, E.G., F.M.P. and E.C.; supervision, E.G.; project administration, E.G.; funding acquisition, E.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Education, University and Research (MIUR) (Dipartimenti di Eccellenza Program (2018–2022)—Department of Biology and Biotechnology “L. Spallanzani”, University of Pavia) and by Animal Breeding and Functional Annotation of Genomes (A1201) Grant 2019-67015-29340/Project Accession 1018854 from the USDA National Institute of Food and Agriculture.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw sequencing data from this study have been submitted to the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/ (accessed on 1 February 2023)) under accession number PRJNA922937. The assembled mitochondrial DNA sequence of Przewalski’s horse from this has been submitted to the NCBI Nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide/ (accessed on 1 February 2023)) under accession number OQ285861.

Acknowledgments

We would like to thank Ernest Bailey and Udeni Balasuriya (University of Kentucky, USA) for providing to us the fibroblast cell line from the Przewalski’s horse hybrid, Tosso Leeb (University of Bern, Switzerland) for the horse BAC clones and Alessandro Achilli (University of Pavia, Italy) for helpful suggestions on mitochondrial DNA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allshire, R.C.; Karpen, G.H. Epigenetic regulation of centromeric chromatin: Old dogs, new tricks? Nat. Rev. Genet. 2008, 9, 923–937. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Choo, K.H. Centromerization. Trends Cell Biol. 2000, 10, 182–188. [Google Scholar] [CrossRef] [PubMed]
  3. Marshall, O.J.; Chueh, A.C.; Wong, L.H.; Choo, K.H. Neocentromeres: New insights into centromere structure, disease development, and karyotype evolution. Am. J. Hum. Genet. 2008, 82, 261–282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Earnshaw, W.C.; Migeon, B.R. Three related centromere proteins are absent from the inactive centromere of a stable isodicentric chromosome. Chromosoma 1985, 92, 290–296. [Google Scholar] [CrossRef] [PubMed]
  5. Voullaire, L.E.; Slater, H.R.; Petrovic, V.; Choo, K.H. A functional marker centromere with no detectable alpha-satellite, satellite III, or CENP-B protein: Activation of a latent centromere? Am. J. Hum. Genet. 1993, 52, 1153–1163. [Google Scholar]
  6. Plohl, M.; Luchetti, A.; Mestrović, N.; Mantovani, B. Satellite DNAs between selfishness and functionality: Structure, genomics and evolution of tandem repeats in centromeric (hetero)chromatin. Gene 2008, 409, 72–82. [Google Scholar] [CrossRef]
  7. Nergadze, S.G.; Piras, F.M.; Gamba, R.; Corbo, M.; Cerutti, F.; McCarter, J.G.W.; Cappelletti, E.; Gozzo, F.; Harman, R.M.; Antczak, D.F.; et al. Birth, evolution, and transmission of satellite-free mammalian centromeric domains. Genome Res. 2018, 28, 789–799. [Google Scholar] [CrossRef] [Green Version]
  8. Cappelletti, E.; Piras, F.M.; Badiale, C.; Bambi, M.; Santagostino, M.; Vara, C.; Masterson, T.A.; Sullivan, K.F.; Nergadze, S.G.; Ruiz-Herrera, A.; et al. CENP-A binding domains and recombination patterns in horse spermatocytes. Sci. Rep. 2019, 9, 15800. [Google Scholar] [CrossRef] [Green Version]
  9. Purgato, S.; Belloni, E.; Piras, F.M.; Zoli, M.; Badiale, C.; Cerutti, F.; Mazzagatti, A.; Perini, G.; Della Valle, G.; Nergadze, S.G.; et al. Centromere sliding on a mammalian chromosome. Chromosoma 2015, 124, 277–287. [Google Scholar] [CrossRef] [Green Version]
  10. Peng, S.; Petersen, J.L.; Bellone, R.R.; Kalbfleisch, T.; Kingsley, N.B.; Barber, A.M.; Cappelletti, E.; Giulotto, E.; Finno, C.J. Decoding the Equine Genome: Lessons from ENCODE. Genes 2021, 12, 1707. [Google Scholar] [CrossRef]
  11. Carbone, L.; Nergadze, S.G.; Magnani, E.; Misceo, D.; Francesca Cardone, M.; Roberto, R.; Bertoni, L.; Attolini, C.; Francesca Piras, M.; de Jong, P.; et al. Evolutionary movement of centromeres in horse, donkey, and zebra. Genomics 2006, 87, 777–782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wade, C.M.; Giulotto, E.; Sigurdsson, S.; Zoli, M.; Gnerre, S.; Imsland, F.; Lear, T.L.; Adelson, D.L.; Bailey, E.; Bellone, R.R.; et al. Genome sequence, comparative analysis, and population genetics of the domestic horse. Science 2009, 326, 865–867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Piras, F.M.; Cappelletti, E.; Santagostino, M.; Nergadze, S.G.; Giulotto, E.; Raimondi, E. Molecular Dynamics and Evolution of Centromeres in the Genus Equus. Int. J. Mol. Sci. 2022, 23, 4183. [Google Scholar] [CrossRef] [PubMed]
  14. Cappelletti, E.; Piras, F.M.; Sola, L.; Santagostino, M.; Abdelgadir, W.A.; Raimondi, E.; Lescai, F.; Nergadze, S.G.; Giulotto, E. Robertsonian Fusion and Centromere Repositioning Contributed to the Formation of Satellite-free Centromeres During the Evolution of Zebras. Mol. Biol. Evol. 2022, 39, msac162. [Google Scholar] [CrossRef]
  15. Roberti, A.; Bensi, M.; Mazzagatti, A.; Piras, F.M.; Nergadze, S.G.; Giulotto, E.; Raimondi, E. Satellite DNA at the Centromere is Dispensable for Segregation Fidelity. Genes 2019, 10, 469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Piras, F.M.; Nergadze, S.G.; Magnani, E.; Bertoni, L.; Attolini, C.; Khoriauli, L.; Raimondi, E.; Giulotto, E. Uncoupling of satellite DNA and centromeric function in the genus Equus. PLoS Genet. 2010, 6, e1000845. [Google Scholar] [CrossRef] [Green Version]
  17. Giulotto, E.; Raimondi, E.; Sullivan, K.F. The Unique DNA Sequences Underlying Equine Centromeres. Prog. Mol. Subcell Biol. 2017, 56, 337–354. [Google Scholar] [CrossRef] [PubMed]
  18. Trifonov, V.A.; Musilova, P.; Kulemsina, A.I. Chromosome evolution in Perissodactyla. Cytogenet. Genome Res. 2012, 137, 208–217. [Google Scholar] [CrossRef]
  19. Trifonov, V.A.; Stanyon, R.; Nesterenko, A.I.; Fu, B.; Perelman, P.L.; O’Brien, P.C.; Stone, G.; Rubtsova, N.V.; Houck, M.L.; Robinson, T.J.; et al. Multidirectional cross-species painting illuminates the history of karyotypic evolution in Perissodactyla. Chromosome Res. 2008, 16, 89–107. [Google Scholar] [CrossRef]
  20. Musilova, P.; Kubickova, S.; Vahala, J.; Rubes, J. Subchromosomal karyotype evolution in Equidae. Chromosome Res. 2013, 21, 175–187. [Google Scholar] [CrossRef]
  21. Montefalcone, G.; Tempesta, S.; Rocchi, M.; Archidiacono, N. Centromere repositioning. Genome Res. 1999, 9, 1184–1188. [Google Scholar] [CrossRef] [Green Version]
  22. Rocchi, M.; Archidiacono, N.; Schempp, W.; Capozzi, O.; Stanyon, R. Centromere repositioning in mammals. Heredity 2012, 108, 59–67. [Google Scholar] [CrossRef] [PubMed]
  23. Piras, F.M.; Nergadze, S.G.; Poletto, V.; Cerutti, F.; Ryder, O.A.; Leeb, T.; Raimondi, E.; Giulotto, E. Phylogeny of horse chromosome 5q in the genus Equus and centromere repositioning. Cytogenet. Genome Res. 2009, 126, 165–172. [Google Scholar] [CrossRef]
  24. Wijers, E.R.; Zijlstra, C.; Lenstra, J.A. Rapid evolution of horse satellite DNA. Genomics 1993, 18, 113–117. [Google Scholar] [CrossRef] [PubMed]
  25. Sakagami, M.; Hirota, K.; Awata, T.; Yasue, H. Molecular cloning of an equine satellite-type DNA sequence and its chromosomal localization. Cytogenet. Cell Genet. 1994, 66, 27–30. [Google Scholar] [CrossRef]
  26. Broad, T.E.; Forrest, J.W.; Lewis, P.E.; Pearce, P.D.; Phua, S.H.; Pugh, P.A.; Stewart-Scott, I.A. Cloning of a DNA repeat element from horse: DNA sequence and chromosomal localization. Genome 1995, 38, 1132–1138. [Google Scholar] [CrossRef] [PubMed]
  27. Broad, T.E.; Ede, A.J.; Forrest, J.W.; Lewis, P.E.; Phua, S.H.; Pugh, P.A. Families of tandemly repeated DNA elements from horse: Cloning, nucleotide sequence, and organization. Genome 1995, 38, 1285–1289. [Google Scholar] [CrossRef]
  28. Cerutti, F.; Gamba, R.; Mazzagatti, A.; Piras, F.M.; Cappelletti, E.; Belloni, E.; Nergadze, S.G.; Raimondi, E.; Giulotto, E. The major horse satellite DNA family is associated with centromere competence. Mol. Cytogenet. 2016, 9, 35. [Google Scholar] [CrossRef] [Green Version]
  29. Nergadze, S.G.; Belloni, E.; Piras, F.M.; Khoriauli, L.; Mazzagatti, A.; Vella, F.; Bensi, M.; Vitelli, V.; Giulotto, E.; Raimondi, E. Discovery and comparative analysis of a novel satellite, EC137, in horses and other equids. Cytogenet. Genome Res. 2014, 144, 114–123. [Google Scholar] [CrossRef]
  30. Ryder, O. Genetic studies of Przewalski’s horses and their impact on conservation. In Przewalski’s Horse: The History and Biology of an Endangered Species; Boyd, L., Houpt, K.A., Eds.; The State University of New York Press: Albany, NY, USA, 1994; pp. 75–92. [Google Scholar]
  31. Ransom, J.; Kaczensky, P. Wild Equids: Ecology, Management, and Conservation; JHU Press: Baltimore, MD, USA, 2016; pp. 1–229. [Google Scholar]
  32. Goto, H.; Ryder, O.A.; Fisher, A.R.; Schultz, B.; Kosakovsky Pond, S.L.; Nekrutenko, A.; Makova, K.D. A massively parallel sequencing approach uncovers ancient origins and high genetic variability of endangered Przewalski’s horses. Genome Biol. Evol. 2011, 3, 1096–1106. [Google Scholar] [CrossRef] [Green Version]
  33. Wallner, B.; Brem, G.; Müller, M.; Achmann, R. Fixed nucleotide differences on the Y chromosome indicate clear divergence between Equus przewalskii and Equus caballus. Anim. Genet. 2003, 34, 453–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Orlando, L.; Ginolhac, A.; Zhang, G.; Froese, D.; Albrechtsen, A.; Stiller, M.; Schubert, M.; Cappellini, E.; Petersen, B.; Moltke, I.; et al. Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse. Nature 2013, 499, 74–78. [Google Scholar] [CrossRef] [PubMed]
  35. Der Sarkissian, C.; Ermini, L.; Schubert, M.; Yang, M.A.; Librado, P.; Fumagalli, M.; Jónsson, H.; Bar-Gal, G.K.; Albrechtsen, A.; Vieira, F.G.; et al. Evolutionary Genomics and Conservation of the Endangered Przewalski’s Horse. Curr. Biol. 2015, 25, 2577–2583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gaunitz, C.; Fages, A.; Hanghøj, K.; Albrechtsen, A.; Khan, N.; Schubert, M.; Seguin-Orlando, A.; Owens, I.J.; Felkel, S.; Bignon-Lau, O.; et al. Ancient genomes revisit the ancestry of domestic and Przewalski’s horses. Science 2018, 360, 111–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Orlando, L. Ancient Genomes Reveal Unexpected Horse Domestication and Management Dynamics. Bioessays 2020, 42, e1900164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bouman, D.T.; Bouman, J.G. The history of Przewalski’s Horse. In Przewalski’s Horse: The History and Biology of an Endangered Species; Boyd, L., Houpt, K.A., Eds.; State University of New York Press: Albany, NY, USA, 1994; pp. 5–38. [Google Scholar]
  39. Jiang, Z.; Zong, H. Reintroduction of the Przewalski’s Horse in China: Status Quo and Outlook. Nat. Conserv. Res. 2019, 4, 15–22. [Google Scholar] [CrossRef]
  40. Turghan, M.A.; Jiang, Z.; Niu, Z. An Update on Status and Conservation of the Przewalski’s Horse (Equus ferus przewalskii): Captive Breeding and Reintroduction Projects. Animals 2022, 12, 3158. [Google Scholar] [CrossRef]
  41. Oakenfull, E.A.; Ryder, O.A. Mitochondrial control region and 12S rRNA variation in Przewalski’s horse (Equus przewalskii). Anim. Genet. 1998, 29, 456–459. [Google Scholar] [CrossRef] [Green Version]
  42. Kerekes, V.; Sándor, I.; Nagy, D.; Ozogány, K.; Göczi, L.; Ibler, B.; Széles, L.; Barta, Z. Trends in demography, genetics, and social structure of Przewalski’s horses in the Hortobagy National Park, Hungary over the last 22 years. Glob. Ecol. Conserv. 2021, 25, e01407. [Google Scholar] [CrossRef]
  43. Bakirova, R.T.; Zharkikh, T.L. Programme on Establishing a Semi-Free Population of Przewalski’s Horse in Orenburg State Nature Reserve: The First Successful Project on the Reintroduction of the Species in Russia. Nat. Conserv. Res. 2019, 4, 57–64. [Google Scholar] [CrossRef]
  44. Bernatkova, A.; Oyunsaikhan, G.; Šimek, J.; Komárková, M.; Bobek, M.; Ceacero, F. Influence of weather on the behaviour of reintroduced Przewalski’s horses in the Great Gobi B Strictly Protected Area (Mongolia): Implications for conservation. BMC Zool. 2022, 7. [Google Scholar] [CrossRef]
  45. Oakenfull, E.A.; Lim, H.N.; Ryder, O.A. A survey of equid mitochondrial DNA: Implications for the evolution, genetic diversity and conservation of Equus. Conserv. Genet. 2000, 1, 341–355. [Google Scholar] [CrossRef]
  46. Jansen, T.; Forster, P.; Levine, M.A.; Oelke, H.; Hurles, M.; Renfrew, C.; Weber, J.; Olek, K. Mitochondrial DNA and the origins of the domestic horse. Proc. Natl. Acad. Sci. USA 2002, 99, 10905–10910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Achilli, A.; Olivieri, A.; Soares, P.; Lancioni, H.; Hooshiar Kashani, B.; Perego, U.A.; Nergadze, S.G.; Carossa, V.; Santagostino, M.; Capomaccio, S.; et al. Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. Proc. Natl. Acad. Sci. USA 2012, 109, 2449–2454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ryder, O.A.; Hansen, S.K. Molecular cytogenetics of the Equidae. I. Purification and cytological localization of a (G + C)-rich satellite DNA from Equus przewalskii. Chromosoma 1979, 72, 115–129. [Google Scholar] [CrossRef]
  49. Myka, J.L.; Lear, T.L.; Houck, M.L.; Ryder, O.A.; Bailey, E. FISH analysis comparing genome organization in the domestic horse (Equus caballus) to that of the Mongolian wild horse (E. przewalskii). Cytogenet. Genome Res. 2003, 102, 222–225. [Google Scholar] [CrossRef]
  50. Yang, F.; Fu, B.; O’Brien, P.C.; Robinson, T.J.; Ryder, O.A.; Ferguson-Smith, M.A. Karyotypic relationships of horses and zebras: Results of cross-species chromosome painting. Cytogenet. Genome Res. 2003, 102, 235–243. [Google Scholar] [CrossRef]
  51. Ahrens, E.; Stranzinger, G. Comparative chromosomal studies of E. caballus (ECA) and E. przewalskii (EPR) in a female F1 hybrid. J. Anim. Breed Genet. 2005, 122 (Suppl. S1), 97–102. [Google Scholar] [CrossRef] [PubMed]
  52. Kvist, L.; Niskanen, M. Modern Northern Domestic Horses Carry Mitochondrial DNA Similar to Przewalski’s Horse. J. Mamm. Evol. 2021, 28, 371–376. [Google Scholar] [CrossRef]
  53. Anglana, M.; Bertoni, L.; Giulotto, E. Cloning of a polymorphic sequence from the nontranscribed spacer of horse rDNA. Mamm. Genome 1996, 7, 539–541. [Google Scholar] [CrossRef]
  54. Benirschke, K.; Malouf, N.; Low, R.J.; Heck, H. Chromosome complement: Differences between Equus caballus and Equus przewalskii, Poliakoff. Science 1965, 148, 382–383. [Google Scholar] [CrossRef] [PubMed]
  55. Fry, K.; Salser, W. Nucleotide sequences of HS-alpha satellite DNA from kangaroo rat Dipodomys ordii and characterization of similar sequences in other rodents. Cell 1977, 12, 1069–1084. [Google Scholar] [CrossRef]
  56. Salser, W.; Bowen, S.; Browne, D.; el-Adli, F.; Fedoroff, N.; Fry, K.; Heindell, H.; Paddock, G.; Poon, R.; Wallace, B.; et al. Investigation of the organization of mammalian chromosomes at the DNA sequence level. Fed. Proc. 1976, 35, 23–35. [Google Scholar] [PubMed]
  57. Ryder, O.A.; Epel, N.C.; Benirschke, K. Chromosome banding studies of the Equidae. Cytogenet Cell Genet. 1978, 20, 332–350. [Google Scholar] [CrossRef] [PubMed]
  58. Liu, G.; Xu, C.Q.; Cao, Q.; Zimmermann, W.; Songer, M.; Zhao, S.S.; Li, K.; Hu, D.F. Mitochondrial and pedigree analysis in Przewalski’s horse populations: Implications for genetic management and reintroductions. Mitochondrial DNA 2014, 25, 313–318. [Google Scholar] [CrossRef] [PubMed]
  59. Breen, M.; Downs, P.; Irvin, Z.; Bell, K. Intrageneric amplification of horse microsatellite markers with emphasis on the Przewalski’s horse (E. przewalskii). Anim. Genet. 1994, 25, 401–405. [Google Scholar] [CrossRef] [PubMed]
  60. Liu, G.; Shafer, A.B.A.; Zimmermann, W.; Hu, D.; Wang, W.; Chu, H.; Cao, J.; Zhao, C. Evaluating the reintroduction project of Przewalski’s horse in China using genetic and pedigree data. Biol. Conserv. 2014, 171, 288–298. [Google Scholar] [CrossRef]
  61. Tang, Y.; Liu, G.; Zhao, S.; Li, K.; Zhang, D.; Liu, S.; Hu, D. Major Histocompatibility Complex (MHC) Diversity of the Reintroduction Populations of Endangered Przewalski’s Horse. Genes 2022, 13, 928. [Google Scholar] [CrossRef]
  62. Jónsson, H.; Schubert, M.; Seguin-Orlando, A.; Ginolhac, A.; Petersen, L.; Fumagalli, M.; Albrechtsen, A.; Petersen, B.; Korneliussen, T.S.; Vilstrup, J.T.; et al. Speciation with gene flow in equids despite extensive chromosomal plasticity. Proc. Natl. Acad. Sci. USA 2014, 111, 18655–18660. [Google Scholar] [CrossRef] [Green Version]
  63. Vidale, P.; Magnani, E.; Nergadze, S.G.; Santagostino, M.; Cristofari, G.; Smirnova, A.; Mondello, C.; Giulotto, E. The catalytic and the RNA subunits of human telomerase are required to immortalize equid primary fibroblasts. Chromosoma 2012, 121, 475–488. [Google Scholar] [CrossRef] [Green Version]
  64. Vidale, P.; Piras, F.M.; Nergadze, S.G.; Bertoni, L.; Verini-Supplizi, A.; Adelson, D.; Guérin, G.; Giulotto, E. Chromosomal assignment of six genes (EIF4G3, HSP90, RBBP6, IL8, TERT, and TERC) in four species of the genus Equus. Anim. Biotechnol. 2011, 22, 119–123. [Google Scholar] [CrossRef]
  65. Santagostino, M.; Khoriauli, L.; Gamba, R.; Bonuglia, M.; Klipstein, O.; Piras, F.M.; Vella, F.; Russo, A.; Badiale, C.; Mazzagatti, A.; et al. Genome-wide evolutionary and functional analysis of the Equine Repetitive Element 1: An insertion in the myostatin promoter affects gene expression. BMC Genet. 2015, 16, 126. [Google Scholar] [CrossRef] [Green Version]
  66. Santagostino, M.; Piras, F.M.; Cappelletti, E.; Del Giudice, S.; Semino, O.; Nergadze, S.G.; Giulotto, E. Insertion of Telomeric Repeats in the Human and Horse Genomes: An Evolutionary Perspective. Int. J. Mol. Sci. 2020, 21, 2838. [Google Scholar] [CrossRef] [Green Version]
  67. Nergadze, S.G.; Lupotto, M.; Pellanda, P.; Santagostino, M.; Vitelli, V.; Giulotto, E. Mitochondrial DNA insertions in the nuclear horse genome. Anim. Genet. 2010, 41 (Suppl. S2), 176–185. [Google Scholar] [CrossRef]
  68. Raimondi, E.; Piras, F.M.; Nergadze, S.G.; Di Meo, G.P.; Ruiz-Herrera, A.; Ponsà, M.; Ianuzzi, L.; Giulotto, E. Polymorphic organization of constitutive heterochromatin in Equus asinus (2n = 62) chromosome 1. Hereditas 2011, 148, 110–113. [Google Scholar] [CrossRef]
  69. Vangipuram, M.; Ting, D.; Kim, S.; Diaz, R.; Schüle, B. Skin punch biopsy explant culture for derivation of primary human fibroblasts. J. Vis. Exp. 2013, 77, e3779. [Google Scholar] [CrossRef] [Green Version]
  70. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  71. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, R25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Ramírez, F.; Ryan, D.P.; Grüning, B.; Bhardwaj, V.; Kilpert, F.; Richter, A.S.; Heyne, S.; Dündar, F.; Manke, T. deepTools2: A next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016, 44, W160–W165. [Google Scholar] [CrossRef] [PubMed]
  73. Lopez-Delisle, L.; Rabbani, L.; Wolff, J.; Bhardwaj, V.; Backofen, R.; Grüning, B.; Ramírez, F.; Manke, T. pyGenomeTracks: Reproducible plots for multivariate genomic datasets. Bioinformatics 2021, 37, 422–423. [Google Scholar] [CrossRef] [PubMed]
  74. Jurka, J.; Kapitonov, V.V.; Pavlicek, A.; Klonowski, P.; Kohany, O.; Walichiewicz, J. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 2005, 110, 462–467. [Google Scholar] [CrossRef] [PubMed]
  75. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Subgroup, G.P.D.P. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve years of SAMtools and BCFtools. Gigascience 2021, 10, giab008. [Google Scholar] [CrossRef] [PubMed]
  77. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
Ijms 24 04134 g001 550
Figure 1. DAPI-banded karyotype of a male Equus przewalskii. The karyotype (2n = 66) consists of 24 biarmed autosomes, 40 acrocentric autosomes, a biarmed X chromosome, and the acrocentric Y chromosome. Chromosomes are arranged and numbered according to standard nomenclature.
Figure 1. DAPI-banded karyotype of a male Equus przewalskii. The karyotype (2n = 66) consists of 24 biarmed autosomes, 40 acrocentric autosomes, a biarmed X chromosome, and the acrocentric Y chromosome. Chromosomes are arranged and numbered according to standard nomenclature.
Ijms 24 04134 g001
Ijms 24 04134 g002 550
Figure 2. Chromosomal localization of two horse BAC clones specific for ECA5p and ECA5q. Metaphase spreads from E. caballus (A), E. caballus × E. przewalskii (B), and E. przewalskii (C) hybridized with BAC clone CHORI241-17G8 (ECA5p) on the left and BAC clone CHORI241-25H7 (ECA5q) on the right. Chromosomes were stained with DAPI (blue). The white arrows indicate BAC signals (red).
Figure 2. Chromosomal localization of two horse BAC clones specific for ECA5p and ECA5q. Metaphase spreads from E. caballus (A), E. caballus × E. przewalskii (B), and E. przewalskii (C) hybridized with BAC clone CHORI241-17G8 (ECA5p) on the left and BAC clone CHORI241-25H7 (ECA5q) on the right. Chromosomes were stained with DAPI (blue). The white arrows indicate BAC signals (red).
Ijms 24 04134 g002
Ijms 24 04134 g003 550
Figure 3. Hybridization of satellite DNA probes with metaphase chromosomes of E. przewalskii. Single-color FISH experiments with three satellite DNA probes: (A) 37cen (red signals), (B) 2PI (red signals), and (C) total genomic DNA (green signals). In panel (D), two-color FISH with 2PI probe (red signals) and total genomic DNA (green signals). The white arrows indicate chromosome 10, which is the only chromosome pair lacking any hybridization signal.
Figure 3. Hybridization of satellite DNA probes with metaphase chromosomes of E. przewalskii. Single-color FISH experiments with three satellite DNA probes: (A) 37cen (red signals), (B) 2PI (red signals), and (C) total genomic DNA (green signals). In panel (D), two-color FISH with 2PI probe (red signals) and total genomic DNA (green signals). The white arrows indicate chromosome 10, which is the only chromosome pair lacking any hybridization signal.
Ijms 24 04134 g003
Ijms 24 04134 g004 550
Figure 4. ChIP-seq characterization of centromeres in Equus przewalskii. (A) ChIP-seq identification of the EPR10 satellite-free centromere. ChIP-seq reads from primary fibroblasts of Przewalski’s horse were mapped on the EquCab2.0 horse reference genome. The CENP-A enriched domain is visualized as a peak. The y-axis reports the normalized read counts whereas the x-axis reports the coordinates on the reference genome. Under the peak, heterozygous nucleotide positions detected in the reads from immunoprecipitated DNA are indicated as rectangles. (B) Enrichment of 37cen and 2PI satellite families in the CENP-A immunoprecipitated chromatin. Values were measured as the ratio between normalized read counts (RPKM) in immunoprecipitated and input DNA. The ERE-1 retrotransposon is used as a control.
Figure 4. ChIP-seq characterization of centromeres in Equus przewalskii. (A) ChIP-seq identification of the EPR10 satellite-free centromere. ChIP-seq reads from primary fibroblasts of Przewalski’s horse were mapped on the EquCab2.0 horse reference genome. The CENP-A enriched domain is visualized as a peak. The y-axis reports the normalized read counts whereas the x-axis reports the coordinates on the reference genome. Under the peak, heterozygous nucleotide positions detected in the reads from immunoprecipitated DNA are indicated as rectangles. (B) Enrichment of 37cen and 2PI satellite families in the CENP-A immunoprecipitated chromatin. Values were measured as the ratio between normalized read counts (RPKM) in immunoprecipitated and input DNA. The ERE-1 retrotransposon is used as a control.
Ijms 24 04134 g004
Table
Table 1. Distribution of satellite DNA families on Equus przewalskii chromosomes.
Table 1. Distribution of satellite DNA families on Equus przewalskii chromosomes.
EPR
Chromosome		37cen2PIGenomic DNA
Number of
Labeled
Chromosomes (%)		Number of
Analyzed Chromosomes		Number of
Labeled
Chromosomes (%)		Number of Analyzed ChromosomesNumber of
Labeled
Chromosomes (%)		Number of Analyzed Chromosomes
1cen60 (100)6035 (85.4)4110 (100)10
2cen3 (5.1)5937 (100)3710 (100)10
3cen53 (100)5334 (100)3410 (100)10
4cen57 (100)5727 (71.1)3810 (100)10
5cen51 (100)5136 (100)3610 (100)10
6cen12 (20.3)5936 (100)361 (10)10
7cen52 (100)5238 (100)3810 (100)10
8cen56 (100)5638 (100)3810 (100)10
9cen49 (100)4939 (100)399 (100)9
10cen0 (0)580 (0)340 (0)10
11cen48 (100)4811 (40.7)2713 (100)13
12cen48 (100)4829 (100)299 (100)9
13cen47 (100)4734 (100)349 (100)9
14cen41 (100)4133 (100)338 (100)8
15cen45 (100)4532 (100)329 (100)9
16cen33 (100)3328 (100)288 (100)8
17cen30 (100)3027 (100)278 (100)8
18cen17 (44.7)3828 (100)288 (100)8
19cen45 (100)4525 (100)2510 (100)10
20cen24 (100)2423 (100)238 (100)8
21cen24 (100)2425 (100)258 (100)8
22cen20 (100)2024 (100)248 (100)8
23cen20 (100)2024 (100)248 (100)8
24cen20 (100)2024 (100)248 (100)8
25cen20 (100)2024 (100)248 (100)8
26cen18 (100)1824 (100)248 (100)8
27cen18 (100)1824 (100)248 (100)8
28cen18 (100)1823 (100)237 (100)7
29cen4 (10.5)388 (36.4)2211 (100)11
30cen18 (100)1824 (100)248 (100)8
31cen17 (100)1724 (100)248 (100)8
32cen17 (100)1724 (100)248 (100)8
Xcen29 (100)2915 (71.4)215 (100)5
Xqinter0 (0)290 (0)215 (100)5
Ycen17 (100)173 (25)126 (100)6
Table
Table 2. SNV analysis of the genomic region corresponding to the CENP-A binding domain of chromosome EPR10.
Table 2. SNV analysis of the genomic region corresponding to the CENP-A binding domain of chromosome EPR10.
Position
(chr11, EquCab2.0 Reference)		Reference AlleleAlternative AlleleNumber of Reads
Reference Allele		Number of Reads
Alternative Allele		Total Number of ReadsReference Allele
Frequency		Alternative
Allele
Frequency
27648739TC815230.30.7
27667251TC2015350.60.4
27671967CG1418320.40.6
27680741GC79160.40.6
27680744CT88160.50.5
27680746CT99180.50.5
27712271CT1119300.40.6
27728102TA816240.30.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Piras, F.M.; Cappelletti, E.; Abdelgadir, W.A.; Salamon, G.; Vignati, S.; Santagostino, M.; Sola, L.; Nergadze, S.G.; Giulotto, E. A Satellite-Free Centromere in Equus przewalskii Chromosome 10. Int. J. Mol. Sci. 2023, 24, 4134. https://doi.org/10.3390/ijms24044134

AMA Style

Piras FM, Cappelletti E, Abdelgadir WA, Salamon G, Vignati S, Santagostino M, Sola L, Nergadze SG, Giulotto E. A Satellite-Free Centromere in Equus przewalskii Chromosome 10. International Journal of Molecular Sciences. 2023; 24(4):4134. https://doi.org/10.3390/ijms24044134

Chicago/Turabian Style

Piras, Francesca M., Eleonora Cappelletti, Wasma A. Abdelgadir, Giulio Salamon, Simone Vignati, Marco Santagostino, Lorenzo Sola, Solomon G. Nergadze, and Elena Giulotto. 2023. "A Satellite-Free Centromere in Equus przewalskii Chromosome 10" International Journal of Molecular Sciences 24, no. 4: 4134. https://doi.org/10.3390/ijms24044134

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop