- Open Access
Independent centromere formation in a capricious, gene-free domain of chromosome 13q21 in Old World monkeys and pigs
Genome Biology volume 7, Article number: R91 (2006)
Evolutionary centromere repositioning and human analphoid neocentromeres occurring in clinical cases are, very likely, two stages of the same phenomenon whose properties still remain substantially obscure. Chromosome 13 is the chromosome with the highest number of neocentromeres. We reconstructed the mammalian evolutionary history of this chromosome and characterized two human neocentromeres at 13q21, in search of information that could improve our understanding of the relationship between evolutionarily new centromeres, inactivated centromeres, and clinical neocentromeres.
Chromosome 13 evolution was studied, using FISH experiments, across several diverse superordinal phylogenetic clades spanning >100 million years of evolution. The analysis revealed exceptional conservation among primates (hominoids, Old World monkeys, and New World monkeys), Carnivora (cat), Perissodactyla (horse), and Cetartiodactyla (pig). In contrast, the centromeres in both Old World monkeys and pig have apparently repositioned independently to a central location (13q21). We compared these results to the positions of two human 13q21 neocentromeres using chromatin immunoprecipitation and genomic microarrays.
We show that a gene-desert region at 13q21 of approximately 3.9 Mb in size possesses an inherent potential to form evolutionarily new centromeres over, at least, approximately 95 million years of mammalian evolution. The striking absence of genes may represent an important property, making the region tolerant to the extensive pericentromeric reshuffling during subsequent evolution. Comparison of the pericentromeric organization of chromosome 13 in four Old World monkey species revealed many differences in sequence organization. The region contains clusters of duplicons showing peculiar features.
Human neocentromeres are functional, analphoid centromeres that emerge epigenetically in ectopic chromosomal regions [1–4]. In the majority of cases, neocentromeres appear to rescue the mitotic stability of acentric chromosomal fragments, often giving rise to aneuploidy . Studies on the evolutionary history of human chromosomes have shown that the centromere can reposition along the chromosome without marker order variation or leading to aneuploidy. This phenomenon is known as 'centromere repositioning' (CR) and has been reported in primates [6–11], in non-primate placental mammals [12, 13], marsupials , birds , and plants .
Recently, two cases of 'repositioned' centromeres in otherwise normal individuals were fortuitously discovered [11, 17]. These two cases can be considered as evolutionary centromere repositioning 'in progress' .
Two human neocentromeres were cytogenetically mapped to duplicons that flanked an ancestral centromere in band 15q24-26 that inactivated before hominoid divergence . Additionally, a neocentromere at 3q26.1 was found located in the same chromosomal domain where an evolutionarily new centromere appeared in the Old World monkey (OWM) ancestor . These studies, which used relatively low resolution cytogenetic techniques, suggested an intriguing relationship between human neocentromeres, ancestral centromeres and evolutionarily new centromeres (repositioned centromeres), which could represent the same phenomenon at different stages of fixation.
Sequences underlining some human neocentromeres have been identified using chromatin immunoprecipitation and genomic microarray (ChIP CHIP) analysis [18–20]. However, despite these recent advances, hypotheses to comprehensively explain the neocentromere phenomenon remain elusive. Phylogenomic studies could eventually provide information about the features, mechanisms and processes of evolutionarily new centromere seeding and development. We report here the evolutionary history of chromosome 13, which exhibits the most extensive clustering of neocentromeres of any human chromosome [3–5]. We found that this chromosome has been exceptionally conserved in evolution and we identified a locus, corresponding to human chromosome band 13q21, where CR events independently occurred in the OWM and pig lineages, whose ancestors diverged at least 95 million years ago (Mya). To further delineate the relationship between human neocentromeres and evolutionarily new centromeres, we used high resolution ChIP CHIP technology to determine the position of two human neocentromeres located in band 13q21. This analysis showed that the neocentromeres do not colocalize with each other or with the OWM and pig evolutionarily new centromeres at the sequence level, but instead map within a few Mb of them.
Chromosome 13 evolution
Chromosome 13 evolution was studied by co-hybridizing a panel of 12 single-copy human BAC clones distributed along chromosome 13 (Table 1, clones in bold). The probes were hybridized on metaphase spreads of 11 primate species (see Materials and methods), including great apes and representatives of OWMs and New World monkeys (NWMs). Examples of fluorescence in situ hybridization (FISH) experiments are reported in Figure 1a-f. Mapping comparison, not considering the centromere, clearly showed that humans and ancestors of OWMs and NWMs share an identical marker order arrangement (Figure 2a).
This analysis was also performed in selected mammals for which BAC libraries were available. The sequence encompassed by each human BAC of the basic panel (in bold in Table 1) was searched for conservation against mouse and rat genomes. 'Overgo' probes were designed on the most conserved region of each BAC. The only exception was marker G, for which a highly conserved region was identified 2 Mb apart from the corresponding human marker. The overgo probes were then used to screen BAC libraries from cattle (CHORI-240 library), horse (CHORI-241), pig (CHORI-242), and cat (RPCI-86). This approach facilitated comparative mapping by assembling a panel of mammalian probes (Additional data file 1) orthologous to the loci encompassed by the human BACs. The identification of additional cattle probes took advantage of the collection of BAC clones positioned on the human sequence by BAC ends, as reported in Larkin et al. . Results in non-primate mammals are also reported in Figure 2a, which shows the most parsimonious chromosomal changes necessary to reconstruct chromosome 13 evolution. In a comprehensive phylogenetic analysis of 64 species, Murphy et al.  defined four large superordinal clades of placental mammals, where Carnivora (cat (FCA)), Perissodactyla (horse (ECA)), and Cetartiodactyla (cattle (BTA), and pig (SUS)) belonged to clade IV and primates belonged to clade III. Remarkably, only a small inversion (markers E-F-G, Table 1, Figure 2a) distinguished the marker order of cat, horse, and pig (clade IV) with respect to primates (clade III).
The position of the centromere, operatively defined as the primary constriction of metaphase chromosomes, was found radically displaced, with respect to surrounding markers, in both OWM and pig (Figure 2a,b). We further refined the centromere position in OWMs using several human BACs spanning the interval H to I (BACs H1 to H9) that define a domain in band 13q21 of approximately 3.9 Mb (61,111,769 to 65,282,688 bp) (Table 1). FISH results (Figure 1a-f) showed that the centromeres in each of the four OWMs analyzed (Macaca mulatta (MMU), Papio hamadryas (PHA), Trachypithecus cristatus (TCR), and Cercopithecus aethiops (CAE)) were located in rather distinct chromosomal locations. In both MMU and PHA, several of these BACs were seen as duplicated signals on either side of the centromere. To better analyze the H1 to H9 region in pig (SUS), overgo probes were designed on the sequence of all H1 to H9 BACs. Only overgo probes from H1, H6, and H8 sequences provided positive results in library screening. Examples of co-hybridization FISH experiments using H6 and H8 probes on pig are reported in Figure 1e. These data showed that the pig centromere was also found in the region between probes H1 and H9 (Figure 2b). We can conclude that these CR events to the same chromosomal region were independent, because pig and OWMs diverged more than 90 Mya and belong to different mammalian clades. We extended our analysis to the elephant (Loxodonta africana, Afroteria, clade I), which is an outgroup with respect to clades III and IV . We found that the elephant homolog to chromosome 13 was highly rearranged and provided no information on the original position of the centromere of the ancestral chromosome 13 or on the E-F-G inversion (data not shown).
Cat chromosome A1 is the result of a fusion of two chromosomes corresponding to human chromosomes 13 and 5 (Figure 2a). Marker order of the portion corresponding to human chromosome 13 was substantially conserved with respect to the mammalian ancestor (Figure 2a). In cat, markers N (distal telomeric in most species) and H8 (pericentromeric in OWM and pig) yielded duplicated signals at the FCA-A1 centromere and at the H8 locus (Figure 1f). In cattle (BTA), chromosome BTA12 contained a large paracentric inversion relative to the ancestral chromosome. The breakpoints of this inversion were mapped at the ancestral centromere (marker A) and between markers H and I, where the OWM and pig evolutionarily new centromeres are located.
Chromosome 13 rearrangements in New World monkeys
Reiterative FISH experiments using additional human BAC clones were performed to more finely map the chromosome 13 fission and inversion breakpoints in NWMs (Table 1, Figure 2a). Human BAC probes that spanned breakpoints were identified by hybridization to both sides of the break on separate locations or chromosomes. In the dusky titi (CMO), the fission breakpoint was localized to probe C1 (BAC RP11-136G6, Table 1), which hybridized to the pericentromeric region of both CMO21 and CMO18 (Figure 2a). In the common marmoset (CJA, Callitrichinae), the fission breakpoint between markers F and G was localized to probe G2 (BAC RP11-939G7, Table 1), which hybridized to both the telomere of CJA5 and the centromere of CJA chromosome 1. The breakpoint of a subtelomeric inversion in CJA1 encompassing markers N-M-L was also mapped to probe K2 (RP11-351H1, Table 1), which hybridized to both the breakpoint and the telomere (Figure 2a).
In the wooly monkey (LLA), marker order was substantially conserved with respect to the NWM ancestor, except for the location of marker A between markers D and E, where an evolutionarily new centromere can be hypothesized (Figure 2a). Further analysis showed that a region of approximately 1.5 Mb delimited by markers B1 and A1 (Table 1) was found transposed from the pericentromeric region of the assumed primate ancestral chromosome to the new centromeric region of LLA8 (Figure 2a). Co-hybridization experiments suggested that the orientation of the segment toward the centromere is conserved with respect to humans. Marker B1 appeared to span the transposition breakpoint, while marker A1 was located adjacent to the heterochromatin/pericentromere boundary on the short-arm side. In this context it is worth noting that Cebidae (SSC and CJA) and Atelidae (LLA) diverged after the split of their common ancestor from Pitheciidae (CMO) .
At least four independent cases of human neocentromeres have been observed in band 13q21, and thus we explored whether the position of these neocentromeres corresponds to the evolutionarily new OWM and pig centromeres. Three 13q21 neocentromeres have been observed on inverted duplication (invdup) chromosomes [5, 24] and one on a small neocentric ring13q21 chromosome derived from a paracentric deletion of 13q21 (Figure 3a) [4, 25]. The size of the ring13q21 neocentric chromosome (and thus the region that contains the neocentromere) was determined by FISH mapping to be approximately 11 Mb, bounded by the absence of BACs RP11-468L10 (chr13: 64.0 Mb) and BAC RP11-332E3 (chr13: 75.3 Mb) . The neocentromere on an invdup13q14-qter chromosome (cell line 13a)  was confirmed to be located to this same approximately 11 Mb region by simultaneous FISH with BAC probes and immunofluorescence with antibodies to CENP-C (data not shown). This cytogenetic mapping showed that the 11 Mbp region containing the human 13q21 neocentromeres was overlapping with the locations of the CAE and SUS centromeres, between probes H7 and H9 (Table 1).
A genomic microarray (CHIP) was constructed containing 107 contiguous BACs (Figure 3d) spanning the entire ring13q21 chromosome, from BAC RP11-468L10 to RP11-332E3, inclusive. Neocentromere DNA from the invdup13q14 chromosome was obtained by chromatin immunoprecipitation (ChIP) from the cell line using antibodies to CENP-A. CENP-A ChIP DNA labeled with Cy-5 and input chromatin DNA labeled with Cy-3 were simultaneously hybridized to the genomic microarray, and positive BACs identified by the Cy-5/Cy-3 intensity ratios .
All BACs showed background ratios (log2 ≤ 1.19) except for three contiguous BACs (RP11-209P2 (log2 = 3.74 ± 1.16), -543G6 (4.23 ± 0.91), and -512J14 (3.07 ± 0.61)), which localized the CENP-A binding domain for this neocentromere (Figure 3b). The CENP-A ChIP protocol was technically not possible on the ring 13q21 cell line due to premature nuclear lysis and, therefore, an alternative ChIP protocol using antibodies to CENP-C was performed (see Materials and methods). All BACs showed background ratios (log2 ≤ 0.97) except for BAC RP11-23B16 (log2 = 3.87 ± 0.17), which localized the CENP-C binding domain for this neocentromere (Figure 3c). The primary data for the ChIP on CHIP analysis is provided in the Additional data files 7 to 17. Previous studies have shown that CENP-A and CENP-C co-immunoprecipitated onto the same chromatin at endogenous centromeres . ChIP on CHIP analysis at a 13q32 neocentromere showed that the CENP-A and CENP-C chromatin domains colocalized at the resolution of a BAC array (Additional data files 2, 3, 16, and 17). Therefore, the use of either CENP-A or CENP-C will accurately identify the 13q21 neocentromere DNA using the BAC array.
Thus, these two human 13q21 neocentromeres occupied distinct genomic regions in 13q21.33 separated by approximately 3 Mb (Table 1). Furthermore, this analysis demonstrated that these 13q21 neocentromeres were separated from the OWM centromeric region by approximately 4 to 7 Mb (the position of the centromere in OWMs was assumed to be located in the middle of the H1 to H8 region, at approximately 63 Mb). The size of the CENP binding domains at these neocentromeres was estimated by removing overlapping regions of neighboring BACs that were negative; for example, the overlap of BAC RP11-520F4 and -321F21 was removed from BAC RP11-23B16 (Figure 3d). This may either over- or underestimate the CENP binding domain, depending on the size and resolution of the BAC clones and/or the sensitivity of the ChIP on CHIP. Like other neocentromeres [18–20], the CENP binding regions of these 13q21 neocentromeres, as estimated from the BAC array, were somewhat enriched in percent AT and LINE elements, and reduced in SINE elements relative to the genome averages (invdup13q, 182 kb, 65.0% AT, 30.61% LINE, 9.69% SINE; Ring, 52 kb, 66.1% AT, 22.41% LINE, 6.15% SINE).
13q21 gene content
Chromosome 13 has one of the lowest gene densities for any human chromosome . We investigated the gene content of the approximately 3.9 Mb interval defined by the H1 to H9 clones (chr13: 61,178,154 to 65,079,597) harboring the OWM and pig evolutionarily new neocentromeres, by querying the specific tracks of the UCSC genome browser (UCSC, March 2006 release). RefSeq annotation does not report any gene in this region. AK127969 and AK098560, reported on the 'known genes' track appear as processed pseudogenes. Moreover, the region containing these RNAs is a portion of a duplicated sequence, with two palindromic copies in the region and a third copy at ch13:45,879,285-45,932,723. A copy of SATR2 satellite is localized in each of the three duplicated regions (chr13: 45,933,879 to 45,936,187, chr13: 63,306,505 to 63,309,036, chr13: 63,215,026 to 63,217,655). The genes PCDH20 (protocadherin 20, chr13: 60,881,821 to 60,887,282) and PCDH9 (protocadherin 9 isoform 1 precursor, chr13: 65,774,967 to 66,702,464) delimit this approximately 3.9 Mb gene desert.
Because of the association between ancestral centromeres and pericentromeric duplications [10, 28], we examined the duplication content of the region corresponding to the evolutionarily new centromeres of 13q21 of humans in more detail (61 to 72 Mb). No enrichment of segmental duplications (>1 kb in length, >90% sequence identity) was observed within this particular region of 13q21 (0.5%) when compared to the chromosome 13 average (2.83%). To identify more ancient segmental duplications, we implemented an alternative approach based on whole genome assembly comparisons using BLASTZ (see Materials and methods), which facilitates the detection of shorter and less homologous duplications (>250 bp, >80% sequence identity). Based on this analysis, we identified an excess of short, more divergent pairwise alignments: the number of older duplications is five times that of more recent duplications, while the ratio when compared to the whole chromosome 13 for older duplications is two-fold (Additional data file 4). Sequence similarity searches of these divergent, short duplications show that more than 64% (51/79) of these regions correspond to exonic portions of other genes (ovostatin, olfactory receptors, and so on) or spliced ESTs, although the genes are not annotated as such within 13q21. We identified the remnants of intron-exon structure for 25 of these regions (19.2 kb) consistent with unprocessed pseudogenes, which were duplicated early during primate evolution (Additional data file 4). Several of the alignments show sequence homology to extant pericentromeric regions on human chromosomes 1p, 2p, 6p and 9q (Additional data files 5 and 6). When the interval is refined further to 127 kb (chr13: 63,188,927 to 63,316,389), the number of pairwise alignments increases several orders of magnitude when compared to the chromosome 13 average (18,354 alignments/Mb versus 161 alignments/Mb). However, the actual number of non-redundant duplicated base-pairs increases only moderately, suggesting that there is a limited amount of sequence that has been the target of several independent duplications. Although not definitive, these sequence properties are potentially consistent with an ancient pericentromeric region.
We have investigated the evolutionary history of chromosome 13 in 11 primate species and in selected non-primate mammals by analysis of marker order arrangement and centromere position, using FISH co-hybridization experiments of appropriate panels of BAC clones (Table 1, Figure 2a). If the centromere position is not taken into account, the marker order of the human chromosome 13 is perfectly conserved in squirrel monkey (SSC, NWM), in OWMs, and in hominoids. This form, therefore, is considered ancestral to primates. Cat (Carnivora), horse (Perissodactyla), and pig (Cetartiodactyla), belonging to mammalian clade IV, share substantially the same marker order arrangement as the primate ancestor except for the inversion of the region encompassed by markers E-F-G (Figure 2a). A limited number of additional inversions and translocations accounted for the chromosome 13 arrangements from the other analyzed mammalian species (Figure 2a).
Analysis of the elephant from mammalian clade 1 (Afrotheria), which branched from placental mammals about 105 Mya , was not informative in resolving the origin of chromosome 13. In this respect, however, it is worth noting that Svartman et al.  have reported that the human chromosome 13 painting library yielded a single signal on the chromosome 2 short arm of the Afrotherian short-eared elephant shrew (Macroscelides proboscideus), indicating that the chromosome was a unique entity in the Afrotheria ancestor.
The radiation hybrid data reported by Murphy et al.  did not detect the E-F-G inversion that we detected in cat, perhaps due to the limits of the radiation hybrid data set they used. The detailed physical map of the cow genome  is completely consistent with our data. The analysis of the chicken genome revealed that the entire euchromatic portion of human chromosome 13 is a continuous segment constituting approximately one-third of the chicken chromosome 1 long arm, in inverted orientation (chr1: 132,956,175 to 171,200,311) ; UCSC, February 2004 assembly). Within this segment there is an inversion whose breakpoint limits are, in the human sequence, within the intervals 40,284,579 to 40,380,292 bp (first breakpoint) and 51,628,385 to 51,849,306 bp (second breakpoint), which corresponds to the E-F-G inversion detected in non-primate mammals. Amazingly, these findings strongly suggest that the chromosome 13 marker order is shared by birds and mammals, although they diverged more than 250 Mya .
Hominoids and squirrel monkeys (SSC, NWM; clade III) and horse (clade IV) all have centromeres adjacent to marker A, which may represent a centromere position shared by the ancestor of primates (clade III) and Cetartiodactyla (clade IV). However, in all the studied OWMs (CAE, MMU, TCR, PHA; clade III) and in the pig (Cetartiodactyla, clade IV), the centromere was found to have repositioned to a region corresponding to the human 13q21 chromosomal region. The CR events were investigated in detail using additional probes spanning the region encompassed by markers H and I. The results show an extensive reshuffling of sequences around the centromere among the four studied OWM species (Figure 2b). The locations of the centromeres appear dispersed in an area of approximately 3.9 Mb with respect to the human sequence representing the ancestral form of the region. These centromeres are, very likely, the result of a single CR event in the OWM ancestor. Cercopithecinae/Colobinae diverged approximately 14.4 to 17.9 Mya , while the genera Papio and Macaca branched about 8.6 to 10.9 Mya . The pericentromeric differences, therefore, provide evidence of an unprecedented degree of local centromeric/pericentromeric plasticity accompanied by extensive sequence remodeling. The position of the evolutionarily new centromere in pig falls within this region (Figure 2b). Mammalian clades III and IV diverged around 95 Mya , OWM and hominoids diverged approximately 25 Mya, and Cetartiodactyla (pig) and Perissodactyla (horse) diverged approximately 83 Mya . We conclude that the pig centromere represents an independent centromere repositioning event in the same chromosomal region. This suggests that there are some features, conserved in the mammalian lineages for at least 70 million years, predisposing this region to centromere formation.
The observation that the region where the OWM and pig evolutionarily new centromeres were seeded is completely devoid of genes is a key finding. Studies of human neocentromere cases have shown that the neocentromere does not influence gene expression per se [35–37]. However, the subsequent heterochromatization of the region that invariably follows the evolutionarily new centromere seeding could, in theory, negatively affect gene expression. The finding that the relatively large H1 to H9 interval is completely devoid of genes may permit extensive sequence reshuffling, an inherent property of eukaryotic centromeres. We have recently reported that the centromere of OWM chromosome 3 is an evolutionarily new centromere that was generated after OWM divergence from Hominoidea . Similarly to the present findings, the centromere seeding occurred within an approximately 430 kb region completely devoid of genes. It can be hypothesized, therefore, that the probability that an evolutionarily new centromere will become fixed during evolution is dependent upon the absence of nearby genes whose alteration due to the reshuffling could be selectively disadvantageous. Preliminary results of the characterization of the evolutionarily new centromere of macaque chromosome 6 are consistent with this model (our unpublished data).
The centromere of cat chromosome FCA-A1 was found located, with respect to the location of the ancestral mammalian centromere, at the opposite telomeric region, adjacent to marker N, where this chromosome fused with the homolog of human chromosome 5 (Figure 2a). This CR may represent a case of jumping centromeres, which is not a rare event in acrocentric chromosomes of primates . Murphy et al.  have noted that telomere-to-centromere conversions are not rare in non-primate mammals. The cross-hybridization signals of H8 and N markers in cat, and the fact that the two breakpoints of the paracentric inversion in the cattle chromosome BTA12 fall in the centromeric region and in the region corresponding to the 13q21 in humans, suggest a special connection in mammals between the chromosome 13 ancestral centromere and the 13q21 region.
Additional potential CR events were observed on mammalian chromosomes in this study. On chromosome LLA8, the best explanation for the emergence of the centromere between markers D and E is the transposition of a 1.5 Mb region that contained the original centromere, which remained active (Figure 2a). This event would not represent true CR as previously defined . However, it is possible that the transposition occurred after the CR event as a result of sequence movement from the inactivated, ancestral centromere to the newly formed centromere, as part of an accruing/degrading process. In this light, it is noteworthy that the transposed region is almost entirely composed of human pericentromeric segmental duplications.
To demonstrate a relationship between evolutionarily new centromeres and human neocentromeres, a chromosomal region must be found that contains both, as seen on human chromosome 3q26 . Endogenous centromeres are relatively easy to localize in mammalian chromosomes using cytogenetic techniques due to the primary constriction and, more importantly, the large region of repetitive satellite DNA that inevitably forms there over evolutionary time . However, human neocentromeres are more difficult to localize since they have not accumulated any repetitive DNA. Therefore, we used ChIP CHIP technology to precisely localize human neocentromeres and compared them to evolutionarily new centromeres.
Human chromosome 13 contains the highest number of reported neocentromeres of any human chromosome, which group into two major clusters at 13q21 and at 13q32 [4, 5]. Thus, we examined the correspondence of human 13q21 neocentromeres with the OWM and pig chromosome 13 evolutionarily new centromeres. This analysis showed that two independent 13q21 neocentromeres were located approximately 4 Mb and 7 Mb distal to the OWM/pig centromeres. The present study, which used high resolution 'ChIP on a CHIP' technology (Figure 3) , did not demonstrate a precise co-localization between the neo- and evolutionarily new centromeres on 13q21. The relatively long distance between human neocentromeres and newly emerged evolutionary centromeres could be taken as evidence against a relationship between them. The two neocentromeres at 15q24-26 reported by Ventura et al.  were found to map about 10 Mb apart. This entire area, however, was shown to represent the wide chromosomal region where pericentromeric duplications of an inactivated centromere were dispersed following the chromosomal fission that generated chromosomes 14 and 15 . As far as chromosome 13, the evolutionary history did not suggest any ancestral inactivated centromeres at 13q21. On the other hand, a significant enrichment of older (80% to 90%) segmental duplications corresponding to the exonic portions of genes was observed for a 127 kb portion of the region of 13q21. Although segmental duplications were seen throughout, the enrichment was most significant for this 127 kb segment.
This study adds two additional neocentromeres to the five that have been precisely mapped using chromatin immunoprecipitation . However, comparative sequence analysis of these seven neocentromere sequences revealed no specific features in common to which neocentromerization competence could be ascribed. The centromere forming potential of specific genomic regions may reflect some relatively long-range property of the chromosomal domain, as opposed to the presence of specific sequence elements.
The present study has tracked the extremely conserved evolution of human chromosome 13. The results defined important aspects of the complex scenario of centromere repositioning and human neocentromere emergence. The centromere of this chromosome repositioned in the same 13q21 region in OWMs and pigs, two species that diverged about 95 Mya. Fine-scale mapping of two clinical neocentromeres suggest that this propensity to form neocentromeres persists within the human population. The absence of genes in the region may be a critical component to progression/fixation of the novel centromere. Cross-species comparisons of chromosome 13 pericentromeric regions in OWM unveiled a striking reshuffling activity.
Materials and methods
Metaphase preparations were obtained from cell lines (lymphoblasts or fibroblasts) from the following species. Great apes: common chimpanzee (Pan troglodytes, PTR), gorilla (Gorilla gorilla, GGO), Borneo orangutan (Pongo pygmaeus pygmaeus, PPY-B). OWMs: rhesus monkey (Macaca mulatta, MMU, Cercopithecinae), sacred baboon (Papio hamadryas, PHA, Cercopithecinae), African green monkey (Cercopithecus aethiops, CAE, Cercopithecinae), silvered-leaf monkey (Trachypithecus cristatus, TCR, Colobinae). NWMs: wooly monkey (Lagothrix lagothricha, LLA, Atelinae), common marmoset (Callithrix jacchus, CJA, Callitrichinae), dusky titi (Callicebus moloch, CMO, Callicebinae), and squirrel monkey (Saimiri sciureus, SSC, Cebinae). Non-primate mammals: cattle (Bos taurus, BTA), horse (Equus caballus, ECA), pig (Sus scrofa, SUS), cat (Felix catus, FCA), and mouse (Mus musculus, MUS). EBV transformed human lymphoblasts containing the neocentric ring13q21 and invdup13q14 chromosomes were grown in standard RPMI media containing 10% fetal calf serum and antibiotics.
DNA extraction from BACs was reported previously . FISH experiments were performed essentially as described by Lichter et al. . Digital images were obtained using a Leica DMRXA2 epifluorescence microscope equipped with a cooled CCD camera (Princeton Instruments, Princenton, NJ, USA. Cy3-dCTP, FluorX-dCTP, DEAC, Cy5-dCTP and DAPI fluorescence signals, detected with specific filters, were recorded separately as gray scale images. Pseudocoloring and merging of images were performed using Adobe Photoshop™ software.
Sixteen overgo probes of 36 to 40 bp each were designed on sequences conserved between the human and mouse genomes according to the HomoloGene database  as described in . The probes were hybridized to high-density filters of mammalian BAC libraries (see Results) and the images were analyzed with ArrayVision Ver 6.0 (Imaging Research Inc., Linton, UK) Linton, UK. The sequence and location of overgo probes, along with clones they identified, are reported in Additional data file 1.
CENP-A immunoprecipitation was performed as described in Alonso et al. . CENP-C immunoprecipitation was performed using protocols modified from Oberley et al. . Approximately 5 × 107 growing cells were crosslinked in 0.5% formaldehyde at room temperature for 10 minutes, followed by addition of glycine to 0.125 M for 5 minutes. Cells were washed in cold phosphate-buffered saline and incubated in 2 ml of lysis buffer (5 mM Pipes pH 8.0, 85 mM KCl, 0.5% w/v NP40, 1 mM PMSF = Phenylmethylsulfonyl fluoride and protease inhibitor cocktail (Sigma, Saint Louis, MO, USA) for 10 minutes at 4°C. Cells were centrifuged and resuspended in 1 ml of cold nuclei lysis buffer (50 mM Tris-HCl pH 8.1, 10 mM EDTA, 1% SDS, 0.5 mM PMSF and protease inhibitor cocktail (Sigma), and sonicated to obtain a DNA ladder ranging from 500 to 1,000 bp. The lysate was centrifuged for 10 minutes at 12,000 g, 4°C. The supernatant was adjusted to 1% Triton, 2 mM EDTA, 20 mM Tris-HCl, pH 8.1, 150 mM NaCl, 0.1% SDS, and precleared for 20 minutes at 4°C with 5 μg of rabbit IgG and 2% of blocked Protein G (Amersham Pharmacia Biotech, Piscataway, NJ, USA), followed by centrifugation. Then, 10 μl of rabbit polyclonal anti-CENP-C (Bill Earnshaw, Edinburgh, Scotland) were added and incubated for 4 h at 4°C. The immunocomplexes were recovered by incubation with 6% blocked ProtG for 2 h at 4°C and centrifugation. They were washed consecutively with low salt buffer (0.1% SDS, 1% Triton, 2 mM EDTA, 20 mM Tris pH 8.1, 150 mM NaCl) and high salt buffer (0.1% SDS, 1% Triton, 2 mM EDTA, 20 mM Tris pH 8.1, 500 mM NaCl), LiCl buffer (0.25 M LiCl,1% NP40, 1% deoxycholate, 1 mM EDTA, 10 mM Tris pH 8.1, and TE (10 mM Tris, 1 mM EDTA), and resuspended in 10 mM Tris HCl, pH 7.5, 5 mM EDTA, 0.25% SDS. To reverse the crosslinks, both input chromatin and immunocomplexes were adjusted to 200 mM NaCl and incubated at 60°C for 8 h, in the presence of 150 μg of Proteinase K PCR grade (Roche Applied Science, Indianapolis, IN, USA) and 20 μg of RNAseA (Qiagen, Valencia, CA, USA). DNA was phenol-chloroform extracted, ethanol precipitated with 1 μg of glycogen (Roche Applied Science) and quantified. Sonicated DNA (25 ng) from CENP-C input and immunocomplexes were repaired with 2.4 units of Kinase, 8 units of Klenow, and 8 units of T4 polymerase (NEB, Beverly, MA, USA) for 1 h at 37°C and ligation-mediated PCR using 2 ng per PCR reaction was carried out as described in Alonso et al. . Ligation-mediated aminoallyl-dUTP PCR was carried out as described in .
13q21 BAC contig microarray construction and hybridization
The 11.4 Mb BAC contig containing BACs RP11-468L10 (AL356252) to RP11-332E3 (AL359392) was assembled by the UCSC Genome Bioinformatics group . Microarrays were prepared as described in  by spotting sonicated BAC DNA onto aminosilane-coated glass slides (GAPSII, Corning, NY, USA). A plasmid containing chromosome 17-specific alpha satellite DNA was included as a positive control for ChIP. For microarray hybridization, 2.5 μg of amplified aminoallyl-dUTP PCR product was conjugated with approximately 20 ng of Mono-Reactive-Cy3 (input chromatin) or Mono-Reactive-Cy5 (ChIP) Dye Pack in 100 mM NaHCO3 (Amersham Biosciences, Little Chalfont, UK) for an hour in the dark. For hybridization of microarrays, 2.5 μg of each -Cy3 and -Cy5 labeled probe (specific activity approximately 60 nt/dye) were denatured 10 minutes at 72°C, with 14 μg of Cot, 575 μg of yeast tRNA and 50 μg of sonicated Escherichia coli DNA, in 25 μl of UltrahybTM (Ambion, Austin, TX, USA) followed by 2 h of annealing at 42°C. Slides were warmed to 42°C for half an hour, the hybridization mix placed on the slide, covered with a 22 × 22 LifterSlip (Erie Scientific Company, Porthsmouth, NH, USA) and incubated overnight at 42°C in an ArrayIt hybridization chamber (Telechem International Inc., Sunnyvale, CA, USA). The slides were washed in 50% formamide, 2 × SSC (SSC = solution of trisodium citrate and sodium clorure), 0.1% Tween-20 for 10 minutes at 45°C, followed by 2 × SSC, 0.1% Tween-20 for 15 minutes at 45°C, and 1× phosphate-buffered saine, 0.1% Tween-20 10 minutes at room temperature. For each microarray, the mean and standard deviation of the triplicate normalized ratios (Lowess) were calculated. Spots with a greater than 25% standard deviation (SD) from the mean were rejected (less than 3% of total spots). For each microarray a scale normalization was performed (X-mean of experiment/SD of experiment) as described by Smyth et al. . Positive BACs were identified as those that were >3 SD from the mean and they represent 3 independent experiments.
To serve as negative controls, ChIP experiments were performed on a cell line containing a neocentromere in band 13q32  and hybridized to the 13q21 BAC array. For anti-CENP-A, all BACs showed background ratios of log2 ≤ 0.88, indicating no CENP-A binding in the 13q21 region, while the alpha satellite DNA showed a ratio of log2 = 6.67. For anti-CENP-C, all BACs showed background ratios of log2 < 1.48, indicating no CENP-C binding in the 13q21 region, while the alpha satellite DNA showed a ratio of log2 = 5. All the primary data for the ChIP on CHIP analysis are provided in Additional data files 7 to 17.
Segmental duplication analysis
We used a BLAST-based detection scheme  to identify all pairwise similarities representing duplicated regions (≥1 kb and ≥90% identity) within the finished sequence of chromosome 13 and compared to all other chromosomes in the NCBI genome assembly (May 2004, build 35). Divergence of duplication, and the number of substitutions per site between the two sequences, were calculated using Kimura's two-parameter method, which corrects for multiple events and transversion/transition mutational biases . In order to detect more divergent duplications, a second all-by-all genome BLASTZ discontinuous search was performed within the finished genome to recover more divergent (>80%) and shorter (>250 bp) alignments (XS, unpublished data) .
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a table listing non-primate mammalian BAC clones utilized in the study, and overgo probes used to screen them. Additional data file 2 is a figure showing ChIP on a CHIP analysis at a 13q32 neocentromere, indicating that the CENP-A and CENP-C chromatin domains colocalized at the resolution of a BAC array. Additional data file 3 contains the legend to the figure in Additional data file 2. Additional data file 4 is a table showing RefSeq genes with sequence similarity to 13q21duplicated regions. Additional data file 5 is a figure depicting segmental duplications (80% to 90% sequence identity and >1 kb or 250 bp) from a 11 Mb region (61 Mb to 72 Mb) of 13q21 to other regions of the genome. Additional data file 6 contains the legend to the figure in Additional data file 5. Additional data files 7 to 17 contain the primary data for the ChIP on CHIP analysis as follows. Additional data file 7 lists specifications of the custom made genomic BAC array (Array design). Additional data files 8, 9, 10, and 11 list primary data for three independent CenpA ChIP on a CHIP experiments and combined data for the three biological replicates for the cell line containing the invdup13q14 neocentric chromosome (cell line CHOP13). Additional data files 12, 13, 14, and 15 list primary data for three independent CenpC ChIP on a CHIP experiments and combined data for the three biological replicates for the cell line containing the ring13q21 neocentric chromosome (cell line Ring A). Additional data files 16 and 17 list data for control CENP-A and CENP-C ChIP performed on a cell line containing a 13q32 neocentromere (cell line BBB).
du Sart D, Cancilla MR, Earle E, Mao J, Saffery R, Tainton KM, Kalitsis P, Martin J, Barry AE, Choo KHA: A functional neo centromere formed through activation of a latent human centromere and consisting of non-alpha-satellite DNA. Nat Genet. 1997, 16: 144-153. 10.1038/ng0697-144.
Warburton PE: Epigenetic analysis of kinetochore assembly on variant human centromeres. Trends Genet. 2001, 17: 243-247. 10.1016/S0168-9525(01)02283-1.
Amor DJ, Choo KH: Neocentromeres: role in human disease, evolution, and centromere study. Am J Hum Genet. 2002, 71: 695-714. 10.1086/342730.
Warburton PE: Chromosomal dynamics of human neocentromere formation. Chromosome Res. 2004, 12: 617-626. 10.1023/B:CHRO.0000036585.44138.4b.
Warburton PE, Dolled M, Mahmood R, Alonso A, Li S, Naritomi K, Tohma T, Nagai T, Hasegawa T, Ohashi H, et al: Molecular cytogenetic analysis of eight inversion duplications of human chromosome 13q that each contain a neocentromere. Am J Hum Genet. 2000, 66: 1794-1806. 10.1086/302924.
Montefalcone G, Tempesta S, Rocchi M, Archidiacono N: Centromere repositioning. Genome Res. 1999, 9: 1184-1188. 10.1101/gr.9.12.1184.
Ventura M, Archidiacono N, Rocchi M: Centromere emergence in evolution. Genome Res. 2001, 11: 595-599. 10.1101/gr.152101.
Carbone L, Ventura M, Tempesta S, Rocchi M, Archidiacono N: Evolutionary history of chromosome 10 in primates. Chromosoma. 2002, 111: 267-272.
Eder V, Ventura M, Ianigro M, Teti M, Rocchi M, Archidiacono N: Chromosome 6 phylogeny in primates and centromere repositioning. Mol Biol Evol. 2003, 20: 1506-1512. 10.1093/molbev/msg165.
Ventura M, Mudge JM, Palumbo V, Burn S, Blennow E, Pierluigi M, Giorda R, Zuffardi O, Archidiacono N, Jackson MS, et al: Neocentromeres in 15q24-26 map to duplicons which flanked an ancestral centromere in 15q25. Genome Res. 2003, 13: 2059-2068. 10.1101/gr.1155103.
Ventura M, Weigl S, Carbone L, Cardone MF, Misceo D, Teti M, D'Addabbo P, Wandall A, Björck E, de Jong P, et al: Recurrent sites for new centromere seeding. Genome Res. 2004, 14: 1696-1703. 10.1101/gr.2608804.
Murphy WJ, Larkin DM, Everts-van der Wind A, Bourque G, Tesler G, Auvil L, Beever JE, Chowdhary BP, Galibert F, Gatzke L, et al: Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps. Science. 2005, 309: 613-617. 10.1126/science.1111387.
Carbone L, Nergadze SG, 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. 10.1016/j.ygeno.2005.11.012.
De Leo AA, Guedelha N, Toder R, Voullaire L, Ferguson-Smith MA, O'Brien PC, Graves JA: Comparative chromosome painting between marsupial orders: relationships with a 2n = 14 ancestral marsupial karyotype. Chromosome Res. 1999, 7: 509-517. 10.1023/A:1009233327176.
Kasai F, Garcia C, Arruga MV, Ferguson-Smith MA: Chromosome homology between chicken (Gallus gallus domesticus) and the red-legged partridge (Alectoris rufa); evidence of the occurrence of a neocentromere during evolution. Cytogenet Genome Res. 2003, 102: 326-330. 10.1159/000075770.
Nagaki K, Cheng Z, Ouyang S, Talbert PB, Kim M, Jones KM, Henikoff S, Buell CR, Jiang J: Sequencing of a rice centromere uncovers active genes. Nat Genet. 2004, 36: 138-145. 10.1038/ng1289.
Amor DJ, Bentley K, Ryan J, Perry J, Wong L, Slater H, Choo KH: Human centromere repositioning "in progress". Proc Natl Acad Sci USA. 2004, 101: 6542-6547. 10.1073/pnas.0308637101.
Lo AW, Magliano DJ, Sibson MC, Kalitsis P, Craig JM, Choo KH: A novel chromatin immunoprecipitation and array (cia) analysis identifies a 460-kb cenp-a-binding neocentromere DNA. Genome Res. 2001, 11: 448-457. 10.1101/gr.GR-1676R.
Alonso A, Mahmood R, Li S, Cheung F, Yoda K, Warburton PE: Genomic microarray analysis reveals distinct locations for the CENP-A binding domains in three human chromosome 13q32 neocentromeres. Hum Mol Genet. 2003, 12: 2711-2721. 10.1093/hmg/ddg282.
Lo AW, Craig JM, Saffery R, Kalitsis P, Irvine DV, Earle E, Magliano DJ, Choo KH: A 330 kb CENP-A binding domain and altered replication timing at a human neocentromere. EMBO J. 2001, 20: 2087-2096. 10.1093/emboj/20.8.2087.
Larkin DM, Everts-van der Wind A, Rebeiz M, Schweitzer PA, Bachman S, Green C, Wright CL, Campos EJ, Benson LD, Edwards J, et al: A cattle-human comparative map built with cattle BAC-ends and human genome sequence. Genome Res. 2003, 13: 1966-1972.
Murphy WJ, Eizirik E, Johnson WE, Zhang YP, Ryder OA, O'Brien SJ: Molecular phylogenetics and the origins of placental mammals. Nature. 2001, 409: 614-618. 10.1038/35054550.
Ray DA, Xing J, Hedges DJ, Hall MA, Laborde ME, Anders BA, White BR, Stoilova N, Fowlkes JD, Landry KE, et al: Alu insertion loci and platyrrhine primate phylogeny. Mol Phylogenet Evol. 2005, 35: 117-126. 10.1016/j.ympev.2004.10.023.
Li S, Malafiej P, Levy B, Mahmood R, Field M, Hughes T, Lockhart LH, Wu Z, Huang M, Hirschhorn K, et al: Chromosome 13q neocentromeres: Molecular cytogenetic characterization of three additional cases and clinical spectrum. Am J Med Genet. 2002, 110: 258-267. 10.1002/ajmg.10454.
Knegt AC, Li S, Engelen JJ, Bijlsma EK, Warburton PE: Prenatal diagnosis of a karyotypically normal pregnancy in a mother with a supernumerary neocentric 13q21→13q22 chromosome and balancing reciprocal deletion. Prenat Diagn. 2003, 23: 215-220. 10.1002/pd.559.
Ando S, Yang H, Nozaki N, Okazaki T, Yoda K: CENP-A, -B, and -C chromatin complex that contains the I-type alpha-satellite array constitutes the prekinetochore in HeLa cells. Mol Cell Biol. 2002, 22: 2229-2241. 10.1128/MCB.22.7.2229-2241.2002.
Dunham A, Matthews LH, Burton J, Ashurst JL, Howe KL, Ashcroft KJ, Beare DM, Burford DC, Hunt SE, Griffiths-Jones S, et al: The DNA sequence and analysis of human chromosome 13. Nature. 2004, 428: 522-528. 10.1038/nature02379.
She X, Horvath JE, Jiang Z, Liu G, Furey TS, Christ L, Clark R, Graves T, Gulden CL, Alkan C, et al: The structure and evolution of centromeric transition regions within the human genome. Nature. 2004, 430: 857-864. 10.1038/nature02806.
Springer MS, Murphy WJ, Eizirik E, O'Brien SJ: Placental mammal diversification and the Cretaceous-Tertiary boundary. Proc Natl Acad Sci USA. 2003, 100: 1056-1061. 10.1073/pnas.0334222100.
Svartman M, Stone G, Page JE, Stanyon R: A chromosome painting test of the basal Eutherian karyotype. Chromosome Res. 2004, 12: 45-53. 10.1023/B:CHRO.0000009294.18760.e4.
Everts-van der Wind A, Larkin DM, Green CA, Elliott JS, Olmstead CA, Chiu R, Schein JE, Marra MA, Womack JE, Lewin HA: A high-resolution whole-genome cattle-human comparative map reveals details of mammalian chromosome evolution. Proc Natl Acad Sci USA. 2005, 102: 18526-18531. 10.1073/pnas.0509285102.
Consortium ICGS: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.
Raaum RL, Sterner KN, Noviello CM, Stewart CB, Disotell TR: Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence. J Hum Evol. 2005, 48: 237-257. 10.1016/j.jhevol.2004.11.007.
Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392: 917-920. 10.1038/31927.
Saffery R, Sumer H, Hassan S, Wong LH, Craig JM, Todokoro K, Anderson M, Stafford A, Choo KH: Transcription within a functional human centromere. Mol Cell. 2003, 12: 509-516. 10.1016/S1097-2765(03)00279-X.
Wong NC, Wong LH, Quach JM, Canham P, Craig JM, Song JZ, Clark SJ, Choo KH: Permissive transcriptional activity at the centromere through pockets of DNA hypomethylation. PLoS Genet. 2006, 2: e17-10.1371/journal.pgen.0020017.
Lam AL, Boivin CD, Bonney CF, Rudd MK, Sullivan BA: Human centromeric chromatin is a dynamic chromosomal domain that can spread over noncentromeric DNA. Proc Natl Acad Sci USA. 2006, 103: 4186-4191. 10.1073/pnas.0507947103.
Lichter P, Tang Chang C-J, Call K, Hermanson G, Evans GA, Housman D, Ward DC: High resolution mapping of human chromosomes 11 by in situ hybridization with cosmid clones. Science. 1990, 247: 64-69. 10.1126/science.2294592.
McPherson JD, Marra M, Hillier L, Waterston RH, Chinwalla A, Wallis J, Sekhon M, Wylie K, Mardis ER, Wilson RK, et al: A physical map of the human genome. Nature. 2001, 409: 934-941. 10.1038/35057157.
Bourque G, Pevzner PA: Genome-scale evolution: reconstructing gene orders in the ancestral species. Genome Res. 2002, 12: 26-36.
Oberley MJ, Farnham PJ: Probing chromatin immunoprecipitates with CpG-island microarrays to identify genomic sites occupied by DNA-binding proteins. Methods Enzymol. 2003, 371: 577-596.
University of California Santa Cruz Genome Bioinformatics. [http://www.genome.ucsc.edu]
Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.
Bailey JA, Yavor AM, Massa HF, Trask BJ, Eichler EE: Segmental duplications: organization and impact within the current human genome project assembly. Genome Res. 2001, 11: 1005-1017. 10.1101/gr.GR-1871R.
Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Cambridge Resource Centre for Comparative Genomics. [http://www.vet.cam.ac.uk/genomics]
MIUR (Ministero Italiano della Universita' e della Ricerca; Cluster C03, Prog. L.488/92), and European Commission (INPRIMAT, QLRI-CT-2002-01325) are gratefully acknowledged for financial support. This work was also supported in part by the National Institutes of Health R01 GM 061150 (to PEW) and R01 GM058815 (to EEE). EEE is an investigator of the Howard Hughes Medical Institute. Mammalian cell lines were kindly provided by the Cambridge Resource Centre . Anti-CENP-A kindly provided by Kinya Yoda (Nagoya, Japan) and anti-CENP-C kindly provided by William Earnshaw (University of Edinburgh).