ColorHOR--novel graphical algorithm for fast scan of alpha satellite higher-order repeats and HOR annotation for GenBank sequence of human genome

MOTIVATION: GenBank data are at present lacking alpha satellite higher-order repeat (HOR) annotation. Furthermore, exact HOR consensus lengths have not been reported so far. Given the fast growth of sequence databases in the centromeric region, it is of increasing interest to have efficient tools for computational identification and analysis of HORs from known sequences. RESULTS: We develop a graphical user interface method, ColorHOR, for fast computational identification of HORs in a given genomic sequence, without requiring a priori information on the composition of the genomic sequence. ColorHOR is based on an extension of the key-string algorithm and provides a color representation of the order and orientation of HORs. For the key string, we use a robust 6 bp string from a consensus alpha satellite and its representative nature is tested. ColorHOR algorithm provides a direct visual identification of HORs (direct and/or reverse complement). In more detail, we first illustrate the ColorHOR results for human chromosome 1. Using ColorHOR we determine for the first time the HOR annotation of the GenBank sequence of the whole human genome. In addition to some HORs, corresponding to those determined previously biochemically, we find new HORs in chromosomes 4, 8, 9, 10, 11 and 19. For the first time, we determine exact consensus lengths of HORs in 10 chromosomes. We propose that the HOR assignment obtained by using ColorHOR be included into the GenBank database.     

A human chromosome 9-specific alphoid DNA repeat spatially resolvable from satellite 3 DNA by fluorescent in situ hybridization.

We have isolated a DNA clone (pMR9A) that identifies an alphoid DNA subset specific for chromosome 9. This alphoid subset is characterized by a dimeric organization as revealed by Southern blot analysis after digestion with HaeIII, HinfI, or StuI. Nonradioactive in situ hybridization demonstrated that pMR9A hybridizes only to the centromeric region of chromosome 9 and reveals chromosome 9 aneuploidies in interphase nuclei. In addition, the probe detects quantitative differences in alpha satellite DNA on chromosome 9, but these quantitative differences are not correlated with the size of the heterochromatic region. Double-labeling experiments, using a chromosome 9-specific satellite 3 clone and pMR9A, enabled us spatially to distinguish the alphoid and satellite 3 domains on metaphase chromosomes after treatment of the cultures with 5-azacytidine.     

Alpha Satellite Insertion Close to an Ancestral Centromeric Region

Human centromeres are mainly composed of alpha satellite DNA hierarchically organized as higher-order repeats (HORs). Alpha satellite dynamics is shown by sequence homogenization in centromeric arrays and by its transfer to other centromeric locations, for example, during the maturation of new centromeres. We identified during prenatal aneuploidy diagnosis by fluorescent in situ hybridization a de novo insertion of alpha satellite DNA from the centromere of chromosome 18 (D18Z1) into cytoband 15q26. Although bound by CENP-B, this locus did not acquire centromeric functionality as demonstrated by the lack of constriction and the absence of CENP-A binding. The insertion was associated with a 2.8-kbp deletion and likely occurred in the paternal germline. The site was enriched in long terminal repeats and located ∼10 Mbp from the location where a centromere was ancestrally seeded and became inactive in the common ancestor of humans and apes 20–25 million years ago. Long-read mapping to the T2T-CHM13 human genome assembly revealed that the insertion derives from a specific region of chromosome 18 centromeric 12-mer HOR array in which the monomer size follows a regular pattern. The rearrangement did not directly disrupt any gene or predicted regulatory element and did not alter the methylation status of the surrounding region, consistent with the absence of phenotypic consequences in the carrier. This case demonstrates a likely rare but new class of structural variation that we name “alpha satellite insertion.” It also expands our knowledge on alphoid DNA dynamics and conveys the possibility that alphoid arrays can relocate near vestigial centromeric sites.     

Intragene higher order repeats in neuroblastoma breakpoint family genes distinguish humans from chimpanzees

Much attention has been devoted to identifying genomic patterns underlying the evolution of the human brain and its emergent advanced cognitive capabilities, which lie at the heart of differences distinguishing humans from chimpanzees, our closest living relatives. Here, we identify two particular intragene repeat structures of noncoding human DNA, spanning as much as a hundred kilobases, that are present in human genome but are absent from the chimpanzee genome and other nonhuman primates. Using our novel computational method Global Repeat Map, we examine tandem repeat structure in human and chimpanzee chromosome 1. In human chromosome 1, we find three higher order repeats (HORs), two of them novel, not reported previously, whereas in chimpanzee chromosome 1, we find only one HOR, a 2mer alphoid HOR instead of human alphoid 11mer HOR. In human chromosome 1, we identify an HOR based on 39-bp primary repeat unit, with secondary, tertiary, and quartic repeat units, fully embedded in human hornerin gene, related to regenerating and psoriatric skin. Such an HOR is not found in chimpanzee chromosome 1. We find a remarkable human 3mer HOR organization based on the ~1.6-kb primary repeat unit, fully embedded within the neuroblastoma breakpoint family genes, which is related to the function of the human brain. Such HORs are not present in chimpanzees. In general, we find that human-chimpanzee differences are much larger for tandem repeats, in particularly for HORs, than for gene sequences. This may be of great significance in light of recent studies that are beginning to reveal the large-scale regulatory architecture of the human genome, in particular the role of noncoding sequences. We hypothesize about the possible importance of human accelerated HOR patterns as components in the gene expression multilayered regulatory network.     

Characterisation of a boundary between satellite III and alphoid sequences on human chromosome 10.

Alphoid and satellite III sequences are arranged as large tandem arrays in the centromeric regions of human chromosomes. Several recent studies using in situ hybridisation to investigate the relative positions of these sequences have shown that they occupy adjacent but non-overlapping domains in metaphase chromosomes. We have analysed the DNA sequence at the junction between alphoid and satellite III sequences in a cosmid previously mapped to chromosome 10. The alphoid sequence consists of tandemly arranged dimers which are distinct from the known chromosome 10-specific alphoid family. Polymerase chain reaction experiments confirm the integrity of the sequence data. These results, together with pulsed field gel electrophoresis data place the boundary between alphoid and satellite III sequences in the mapping interval 10 centromere-10q11.2. The sequence data shows that these repetitive sequences are separated by a partial L1 interspersed repeat sequence less than 500bp in length. The arrangement of the junction suggests that a recombination event has brought these sequences into close proximity.     

Structure, organization, and sequence of alpha satellite DNA from human chromosome 17: evidence for evolution by unequal crossing-over and an ancestral pentamer repeat shared with the human X chromosome.

The centromeric regions of all human chromosomes are characterized by distinct subsets of a diverse tandemly repeated DNA family, alpha satellite. On human chromosome 17, the predominant form of alpha satellite is a 2.7-kilobase-pair higher-order repeat unit consisting of 16 alphoid monomers. We present the complete nucleotide sequence of the 16-monomer repeat, which is present in 500 to 1,000 copies per chromosome 17, as well as that of a less abundant 15-monomer repeat, also from chromosome 17. These repeat units were approximately 98% identical in sequence, differing by the exclusion of precisely 1 monomer from the 15-monomer repeat. Homologous unequal crossing-over is suggested as a probable mechanism by which the different repeat lengths on chromosome 17 were generated, and the putative site of such a recombination event is identified. The monomer organization of the chromosome 17 higher-order repeat unit is based, in part, on tandemly repeated pentamers. A similar pentameric suborganization has been previously demonstrated for alpha satellite of the human X chromosome. Despite the organizational similarities, substantial sequence divergence distinguishes these subsets. Hybridization experiments indicate that the chromosome 17 and X subsets are more similar to each other than to the subsets found on several other human chromosomes. We suggest that the chromosome 17 and X alpha satellite subsets may be related components of a larger alphoid subfamily which have evolved from a common ancestral repeat into the contemporary chromosome-specific subsets.     

Organization and evolution of an alpha satellite DNA subset shared by human chromosomes 13 and 21

The structure of the alpha satellite DNA higher-order repeat (HOR) unit from a subset shared by human chromosomes 13 and 21 (D13Z1 and D21Z1) has been examined in detail. By using a panel of hybrids possessing either a chromosome 13 or a chromosome 21, different HOR unit genotypes on chromosomes 13 and 21 have been distinguished. We have also determined the basis for a variant HOR unit structure found on 8% of chromosomes 13 but not at all on chromosomes 21. Genomic restriction maps of the HOR units found on the two chromosome 13 genotypes and on the chromosome 21 genotype are constructed and compared. The nucleotide sequence of a predominant 1.9-kilobasepair HOR unit from the D13Z1/D21Z1 subset has been determined. The DNA sequences of different alpha satellite monomers comprising the HOR are compared, and the data are used to develop a model, based on unequal crossing-over, for the evolution of the current HOR unit found at the centromeres of both these chromosomes.Correspondence to: H.F. Willard     

Molecular analysis of a deletion polymorphism in alpha satellite of human chromosome 17: evidence for homologous unequal crossing-over and subsequent fixation.

The human alpha satellite DNA family is organized into chromosome-specific subsets characterized by distinct higher-order repeats based on a approximately 171 basepair monomer unit. On human chromosome 17, the predominant form of alpha satellite is a 16-monomer (16-mer) higher-order repeat present in 500-1000 copies per chromosome 17. In addition, less abundant 15-monomer and 14-monomer repeats are also found constitutively on chromosome 17. Polymorphisms in the form of different higher-order repeat lengths have been described for this subset, the most prominent polymorphism being a 13-monomer (13-mer) higher-order repeat present on approximately 35% of all chromosomes 17. To investigate the nature of this polymorphism, we have cloned, sequenced and compared the relevant regions of the 13-mer to the previously characterized 16-mer repeat. The results show that the repeats are virtually identical, with the principal difference being the exclusion of three monomers from the 13-mer repeat. We propose that the 13-mer is the product of an isolated homologous recombination event between two monomers of the 16-mer repeat. Sequence comparisons reveal the approximate site of recombination and flanking regions of homology. This recombination site corresponds to a position within the alphoid monomer which has been previously implicated in an independent homologous recombination event, suggesting that there may exist a preferred register for recombination in alphoid DNA. We suggest that these events are representative of an ongoing process capable of reorganizing the satellite subset of a given chromosome, thereby contributing to the establishment of chromosome-specific alpha satellite subsets.     

Definition of a second dimeric subfamily of human alpha satellite DNA

We describe a new human subfamily of alpha satellite DNA. The restriction endonuclease XbaI cleaves this subfamily into a collection of fragments which are heterogeneous with respect to size. We compared the sequences of 6 clones from four different XbaI size classes. Clones from a single size class were not necessarily more related than clones from different classes. Clones from different size classes were found to produce almost identical hybridization patterns with XbaI-digested human genomic DNA. All clones were found to share a common dimeric repeat organization, with dimers exhibiting about 84% sequence identities, indicating that the clones evolved from a common progenitor alphoid dimer. We show that this subfamily, and the EcoRI dimer subfamily originally described by Wu and Manuelidis, evolved from different progenitor alphoid dimers, and therefore represent distinct human alphoid subfamilies.     

Nucleotide sequence heterogeneity of alpha satellite repetitive DNA: a survey of alphoid sequences from different human chromosomes.

The human alpha satellite DNA family is composed of diverse, tandemly reiterated monomer units of approximately 171 basepairs localized to the centromeric region of each chromosome. These sequences are organized in a highly chromosome-specific manner with many, if not all human chromosomes being characterized by individually distinct alphoid subsets. Here, we compare the nucleotide sequences of 153 monomer units, representing alphoid components of at least 12 different human chromosomes. Based on the analysis of sequence variation at each position within the 171 basepair monomer, we have derived a consensus sequence for the monomer unit of human alpha satellite DNA which we suggest may reflect the monomer sequence from which different chromosomal subsets have evolved. Sequence heterogeneity is evident at each position within the consensus monomer unit and there are no positions of strict nucleotide sequence conservation, although some regions are more variable than others. A substantial proportion of the overall sequence variation may be accounted for by nucleotide changes which are characteristic of monomer components of individual chromosomal subsets or groups of subsets which have a common evolutionary history.     

The role of alphoid higher order repeats (HORs) in the centromere folding

Understanding the folding of centromere DNA in the maximally condensed methaphase chromosome remains a basic challenge in cell biology. We propose here a set of structural models with a graphical presentation of alphoid higher order repeat (HOR) distribution in the centromere folding, based on the assumption of encryption key for microtubule-centromere interaction which arises from chromosome-specific crystal-like structure of HORs. Specific HOR leads to a characteristic geometrical pattern which may be responsible for individual microtubule to recognize a specific structure of centromere in each chromosome.     

Consensus higher order repeats and frequency of string distributions in human genome

Key string algorithm (KSA) could be viewed as robust computational generalization of restriction enzyme method. KSA enables robust and effective identification and structural analyzes of any given genomic sequences, like in the case of NCBI assembly for human genome. We have developed a method, using total frequency distribution of all r-bp key strings in dependence on the fragment length l, to determine the exact size of all repeats within the given genomic sequence, both of monomeric and HOR type. Subsequently, for particular fragment lengths equal to each of these repeat sizes we compute the partial frequency distribution of r-bp key strings; the key string with highest frequency is a dominant key string, optimal for segmentation of a given genomic sequence into repeat units. We illustrate how a wide class of 3-bp key strings leads to a key-string-dependent periodic cell which enables a simple identification and consensus length determinations of HORs, or any other highly convergent repeat of monomeric or HOR type, both tandem or dispersed. We illustrated KSA application for HORs in human genome and determined consensus HORs in the Build 35.1 assembly. In the next step we compute suprachromosomal family classification and CENP-B box / pJalpha distributions for HORs. In the case of less convergent repeats, like for example monomeric alpha satellite (20-40% divergence), we searched for optimal compact key string using frequency method and developed a concept of composite key string (GAAAC--CTTTG) or flexible relaxation (28 bp key string) which provides both monomeric alpha satellites as well as alpha monomer segmentation of internal HOR structure. This method is convenient also for study of R-strand (direct) / S-strand (reverse complement) alpha monomer alternations. Using KSA we identified 16 alternating regions of R-strand and S-strand monomers in one contig in choromosome 7. Use of CENP-B box and/or pJalpha motif as key string is suitable both for identification of HORs and monomeric pattern as well as for studies of CENP-B box / pJalpha distribution. As an example of application of KSA to sequences outside of HOR regions we present our finding of a tandem with highly convergent 3434-bp Long monomer in chromosome 5 (divergence less then 0.3%).     

Chromosome-specific alpha satellite DNA: isolation and mapping of a polymorphic alphoid repeat from human chromosome 10

Distinct subsets of the human alpha satellite repetitive DNA family can be found in the centromeric region of each chromosome. Here we described the isolation and mapping of an alpha satellite repeat unit specific for human chromosome 10, using a somatic cell hybrid in which the only human centromere derives from chromosome 10. A hierarchical higher-order repeat unit, consisting of eight tandem approximately 171-bp alphoid monomer units, is defined by six restriction endonucleases. Under high-stringency conditions, a cloned representative of this 8-mer repeat family hybridizes to chromosome 10 only, both by Southern blot analysis of a somatic cell hybrid panel and by in situ hybridization. The probe furthermore detects a polymorphic restriction pattern of the alpha satellite array on chromosome 10. These features will make this probe a valuable genetic marker for studies of the centromeric region of chromosome 10.     

Study of alpha-satellite DNA in cosmid libraries, specific for chromosomes 13, 21, and 22, using fluorescence in situ hybridization

Fluorescent in situ hybridization (FISH) was employed in mapping the alpha-satellite DNA that was revealed in the cosmid libraries specific for human chromosomes 13, 21, and 22. In total, 131 clones were revealed. They contained various elements of centromeric alphoid DNA sequences of acrocentric chromosomes, including those located close to SINEs, LINEs, and classical satellite sequences. The heterochromatin of acrocentric chromosomes was shown to contain two different groups of alphoid sequences: (1) those immediately adjacent to the centromeric regions (alpha 13-1, alpha 21-1, and alpha 22-1 loci) and (2) those located in the short arm of acrocentric chromosomes (alpha 13-2, alpha 21-2, and alpha 22-2 loci). Alphoid DNA sequences from the alpha 13-2, alpha 21-2, and alpha 22-2 loci are apparently not involved in the formation of centromeres and are absent from mitotically stable marker chromosomes with a deleted short arm. Robertsonian translocations t(13q; 21q) and t(14q; 22q), and chromosome 21p-. The heterochromatic regions of chromosomes 13, 21, and 22 were also shown to contain relatively chromosome-specific repetitive sequences of various alphoid DNA families, whose numerous copies occur in other chromosomes. Pools of centromeric alphoid cosmids can be of use in further studies of the structural and functional properties of heterochromatic DNA and the identification of centromeric sequences. Moreover, these clones can be employed in high-resolution mapping and in sequencing the heterochromatic regions of the human genome. The detailed FISH analysis of numerous alphoid cosmid clones allowed the identification of several new, highly specific DNA probes of molecular cytogenetic studies--in particular, the interphase and metaphase analyses of chromosomes 2, 9, 11, 14, 15, 16, 18, 20, 21-13, 22-14, and X.     

Identification of DNA sequences flanking the breakpoint of human t(14q21q) Robertsonian translocations.

We have employed molecular probes and in situ hybridization to investigate the DNA sequences flanking the breakpoint of a group of t(14q21q) Robertsonian translocations. In all the families studied, the probands were patients with Down syndrome who carried a de novo t(14q21q) translocation. The DNA probes used were two alphoid sequences, alphaRI and alphaXT, which are specific for the centromeres of chromosomes 13 and 21 and of chromosomes 14 and 22, respectively; a satellite III sequence, pTRS-47, which is specific for the proximal p11 region of chromosomes 14 and 22; and a newly defined satellite III DNA, pTRS-63, which is specific for the distal p11 region of chromosome 14. The two alphoid probes detected approximately the same amount of autoradiographic signal on the translocated chromosomes as was expected for chromosomes 14 and 21 of the originating parent, suggesting that there has been no loss of these centromeric sequences during the translocation events. Results with the two satellite III probes indicated that the domain corresponding to pTRS-47 was retained in the translocated chromosomes, whereas the domain for pTRS-63 was lost. These results have allowed us to place the translocation breakpoint between the pTRS-47 and pTRS-63 domains within the p11 region of chromosome 14.     

Rapid generation of chromosome-specific alphoid DNA probes using the polymerase chain reaction

Summary Non-isotopic in situ hybridization of chromosome-specific alphoid DNA probes has become a potent tool in the study of numerical aberrations of specific human chromosomes at all stages of the cell cycle. In this paper, we describe approaches for the rapid generation of such probes using the polymerase chain reaction (PCR), and demonstrate their chromosome specificity by fluorescence in situ hybridization to normal human metaphase spreads and interphase nuclei. Oligonucleotide primers for conserved regions of the alpha satellite monomer were used to generate chromosome-specific DNA probes from somatic hybrid cells containing various human chromosomes, and from DNA libraries from sorted human chromosomes. Oligonucleotide primers for chromosome-specific regions of the alpha satellite monomer were used to generate specific DNA probes for the pericentromeric heterochromatin of human chromosomes 1, 6, 7, 17 and X directly from human genomic DNA.     

Labeling of the centromeric region on human chromosome 8 by in situ hybridization

Summary Probe DNA that binds preferentially to the centromeric region of human chromosomes 8 was synthesized. Alpha satellite probe DNA molecules were selectively amplified from sorter-purified human chromosomes 8 by in vitro DNA amplification using the polymerase chain reaction (PCR). Probe labeling was performed during PCR by incorporation of biotinylated deoxyuridine. In situ hybridization of unpurified probe DNA comprised of alpha satellite monomer and higher molecular weight DNA fragments with metaphase chromosome spreads showed binding to the centromeric regions of numerous chromosomes. However, blocking with unlabeled total human alphoid DNA dramatically reduced crosshybridization to chromosomes other than 8. Under these conditions, the degenerate probe DNA allowed unambiguous visualization of domains occupied by centromeric DNA of chromosome 8 in metaphase spreads and interphase cell nuclei, thus greatly facilitating the detection of numerical chromosome aberrations in tumor cells. In situ hybridization of size-fractionated alpha satellite DNA identified the monomeric fraction as the major cause of crosshybridization. Alpha satellite dimers and higher molecular weight DNA fragments showed relatively high specificity for human chromosomes 8.     

Comparative mapping of a gorilla-derived alpha satellite DNA clone on great ape and human chromosomes

We have isolated an alpha satellite DNA clone, pG3.9, from gorilla DNA. Fluorescence in situ hybridization on banded chromosomes under high stringency conditions revealed that pG3.9 identifies homologous sequences at the centromeric region of ten gorilla chromosomes, and, with few exceptions, also recognizes the homologous chromosomes in human. A pG3.9-like alphoid DNA is present on a larger number of orangutan chromosomes, but, in contrast, is present on only tow chromosomes in the chimpanzee. These results show that the chromosomal subsets of related alpha satellite DNA sequences may undergo different patterns of evolution.by J.B. Rattner     

Oligonucleotide-primed in situ DNA synthesis (PRINS): a method for chromosome mapping, banding, and investigation of sequence organization

Oligonucleotides were annealed to complementary sequences in fixed human metaphase chromosomes and extended with DNA polymerase. The newly synthesized fragments were labeled by incorporating bio-11-dUTP instead of TTP, and the sites of synthesis were detected by immunocytochemistry, using fluorochromes as the reporter molecules. We have obtained clear localization with oligonucleotides from alphoid (centromeric sequences), simple sequence (satellite) DNAs, a variety of Alu-dispersed repeated sequences, and oligonucleotides derived from the Tetrahymena and Trypanosoma telomere-specific sequences. The simple sequence and alphoid oligonucleotides gave results at least comparable to those obtained using the whole molecule as a probe for in situ hybridization, whereas the Alu oligonucleotides produced a diversity of results which depended on the absolute length and location of the oligonucleotide within the Alu sequence. The telomere-specific oligomers also produced a variety of results. The G-rich Trypanosoma oligomer and its complementary C-rich sequence produced strong telomeric signals and some interstitial signals on mouse chromosomes, but only weak telomeric signals on human chromosomes. The G-rich Tetrahymena oligomer produced detectable telomeric signals on human chromosomes. The technique appears to be a valuable extension of present tools for mapping and examining the organization of DNA sequences within chromosomes.