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

Recombination between heterologous human acrocentric chromosomes

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

The complete sequence of a human genome

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

Dark Matter of Primate Genomes: Satellite DNA Repeats and Their Evolutionary Dynamics

A substantial portion of the primate genome is composed of non-coding regions, so-called "dark matter", which includes an abundance of tandemly repeated sequences called satellite DNA. Collectively known as the satellitome, this genomic component offers exciting evolutionary insights into aspects of primate genome biology that raise new questions and challenge existing paradigms. A complete human reference genome was recently reported with telomere-to-telomere human X chromosome assembly that resolved hundreds of dark regions, encompassing a 3.1 Mb centromeric satellite array that had not been identified previously. With the recent exponential increase in the availability of primate genomes, and the development of modern genomic and bioinformatics tools, extensive growth in our knowledge concerning the structure, function, and evolution of satellite elements is expected. The current state of knowledge on this topic is summarized, highlighting various types of primate-specific satellite repeats to compare their proportions across diverse lineages. Inter- and intraspecific variation of satellite repeats in the primate genome are reviewed. The functional significance of these sequences is discussed by describing how the transcriptional activity of satellite repeats can affect gene expression during different cellular processes. Sex-linked satellites are outlined, together with their respective genomic organization. Mechanisms are proposed whereby satellite repeats might have emerged as novel sequences during different evolutionary phases. Finally, the main challenges that hinder the detection of satellite DNA are outlined and an overview of the latest methodologies to address technological limitations is presented.

Centromere studies in the era of 'telomere-to-telomere' genomics

We are entering into an exciting era of genomics where truly complete, high-quality assemblies of human chromosomes are available end-to-end, or from 'telomere-to-telomere' (T2T). This technological advance offers a new opportunity to include endogenous human centromeric regions in high-resolution, sequence-based studies. These emerging reference maps are expected to reveal a new functional landscape in the human genome, where centromere proteins, transcriptional regulation, and spatial organization can be examined with base-level resolution across different stages of development and disease. Such studies will depend on innovative assembly methods of extremely long tandem repeats (ETRs), or satellite DNAs, paired with the development of new, orthogonal validation methods to ensure accuracy and completeness. This review reflects the progress in centromere genomics, credited by recent advancements in long-read sequencing and assembly methods. In doing so, I will discuss the challenges that remain and the promise for a new period of scientific discovery for satellite DNA biology and centromere function.

Centromere Repeats: Hidden Gems of the Genome

Satellite DNAs are now regarded as powerful and active contributors to genomic and chromosomal evolution. Paired with mobile transposable elements, these repetitive sequences provide a dynamic mechanism through which novel karyotypic modifications and chromosomal rearrangements may occur. In this review, we discuss the regulatory activity of satellite DNA and their neighboring transposable elements in a chromosomal context with a particular emphasis on the integral role of both in centromere function. In addition, we discuss the varied mechanisms by which centromeric repeats have endured evolutionary processes, producing a novel, species-specific centromeric landscape despite sharing a ubiquitously conserved function. Finally, we highlight the role these repetitive elements play in the establishment and functionality of de novo centromeres and chromosomal breakpoints that underpin karyotypic variation. By emphasizing these unique activities of satellite DNAs and transposable elements, we hope to disparage the conventional exemplification of repetitive DNA in the historically-associated context of ‘junk’.

Alpha satellite DNA biology: finding function in the recesses of the genome

Repetitive DNA, formerly referred to by the misnomer "junk DNA," comprises a majority of the human genome. One class of this DNA, alpha satellite, comprises up to 10% of the genome. Alpha satellite is enriched at all human centromere regions and is competent for de novo centromere assembly. Because of the highly repetitive nature of alpha satellite, it has been difficult to achieve genome assemblies at centromeres using traditional next-generation sequencing approaches, and thus, centromeres represent gaps in the current human genome assembly. Moreover, alpha satellite DNA is transcribed into repetitive noncoding RNA and contributes to a large portion of the transcriptome. Recent efforts to characterize these transcripts and their function have uncovered pivotal roles for satellite RNA in genome stability, including silencing "selfish" DNA elements and recruiting centromere and kinetochore proteins. This review will describe the genomic and epigenetic features of alpha satellite DNA, discuss recent findings of noncoding transcripts produced from distinct alpha satellite arrays, and address current progress in the functional understanding of this oft-neglected repetitive sequence. We will discuss unique challenges of studying human satellite DNAs and RNAs and point toward new technologies that will continue to advance our understanding of this largely untapped portion of the genome.

Higher-order repeat structure in alpha satellite DNA occurs in New World monkeys and is not confined to hominoids

Centromeres usually contain large amounts of tandem repeat DNA. Alpha satellite DNA (AS) is the most abundant tandem repeat DNA found in the centromeres of simian primates. The AS of humans contains sequences organized into higher-order repeat (HOR) structures, which are tandem arrays of larger repeat units consisting of multiple basic repeat units. HOR-carrying AS also occurs in other hominoids, but results reported to date for phylogenetically more remote taxa have been negative. Here we show direct evidence for clear HOR structures in AS of the owl monkey and common marmoset. These monkeys are New World monkey species that are located phylogenetically outside of hominoids. It is currently postulated that the presence of HOR structures in AS is unique to hominoids. Our results suggest that this view must be modified. A plausible explanation is that generation of HOR structures is a general event that occurs occasionally or frequently in primate centromeres, and that, in humans, HOR-carrying AS became predominant in the central region of the centromere. It is often difficult to assemble sequence reads of tandem repeat DNAs into accurate contig sequences; our careful sequencing strategy allowed us to overcome this problem.

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%).

In situ aneuploidy assessment in human sperm: the use of primed in situ and peptide nucleic acid-fluorescence in situ hybridization techniques

Both the primed in situ (PRINS) and the peptide nucleic acid-fluorescence in situ hybridization (PNA-FISH) techniques constitute alternatives to the conventional (fluorescence in situ hybridization, FISH) procedure for chromosomal investigations. The PRINS reaction is based on the use of a DNA polymerase and labeled nucleotide in an in situ primer extension reaction. Peptide nucleic acid probes are synthetic DNA analogs with uncharged polyamide backbones. The two procedures present several advantages (specificity, rapidity and discriminating ability) that make them very attractive for cytogenetic purposes. Their adaptation to human spermatozoa has allowed the development of new and fast procedures for the chromosomal screening of male gametes and has provided efficient complements to FISH for in situ assessment of aneuploidy in male gametes.

Localization of MMR proteins on meiotic chromosomes in mice indicates distinct functions during prophase I

Mammalian MutL homologues function in DNA mismatch repair (MMR) after replication errors and in meiotic recombination. Both functions are initiated by a heterodimer of MutS homologues specific to either MMR (MSH2-MSH3 or MSH2-MSH6) or crossing over (MSH4-MSH5). Mutations of three of the four MutL homologues (Mlh1, Mlh3, and Pms2) result in meiotic defects. We show herein that two distinct complexes involving MLH3 are formed during murine meiosis. The first is a stable association between MLH3 and MLH1 and is involved in promoting crossing over in conjunction with MSH4-MSH5. The second complex involves MLH3 together with MSH2-MSH3 and localizes to repetitive sequences at centromeres and the Y chromosome. This complex is up-regulated in Pms2-/- males, but not females, providing an explanation for the sexual dimorphism seen in Pms2-/- mice. The association of MLH3 with repetitive DNA sequences is coincident with MSH2-MSH3 and is decreased in Msh2-/- and Msh3-/- mice, suggesting a novel role for the MMR family in the maintenance of repeat unit integrity during mammalian meiosis.

Higher-order organization and compartmentalization of satellite DNA PIM357 in species of the coleopteran genus Pimelia

The PIM357 satellite DNA family is present in 26 Pimelia taxa (Tenebrionidae, Coleoptera) with endemic congeneric species from the Canary Islands showing higher interrepeat variability than continental ones. In this paper, we compare the repetitive DNA sequences of a Canarian species that has distinct subfamilies of repeat units, P. radula ascendens, with another without such subfamilies, P. sparsa sparsa. The chromosomal localization of the repeat units and the comparison of the variability of randomly cloned monomers to the one estimated by comparing repeat units from dimers and trimers suggest the absence of satellite subfamilies in P. sparsa sparsa. Hence, the repeat units of this species seem to be uniformly and randomly distributed throughout all chromosomes out of one chromosomal pair. On the contrary, P. radula ascendens shows four divergent subfamilies of repeat units supported by several diagnostic nucleotide substitutions. These subfamilies seem to form four distinct repeat units: monomer subfamily 1, monomer subfamily 4 and two higher-order units (dimer linking subfamily 1 and 4, and dimer linking subfamily 2 and 3). Moreover, monomers of subfamily 1 are present in three chromosomal pairs only. We discuss the effect of different potential factors acting in the concerted evolution and the genomic organization of stDNA sequences in these taxa.

Evidence for a fast, intrachromosomal conversion mechanism from mapping of nucleotide variants within a homogeneous alpha-satellite DNA array

Assuming that patterns of sequence variants within highly homogeneous centromeric tandem repeat arrays can tell us which molecular turnover mechanisms are presently at work, we analyzed the alpha-satellite tandem repeat array DXZ1 of one human X chromosome. Here we present accurate snapshots from this dark matter of the genome. We demonstrate stable and representative cloning of the array in a P1 artificial chromosome (PAC) library, use samples of higher-order repeats subcloned from five unmapped PACs (120-160 kb) to identify common variants, and show that such variants are presently in a fixed transition state. To characterize patterns of variant spread throughout homogeneous array segments, we use a novel partial restriction and pulsed-field gel electrophoresis mapping approach. We find an older large-scale (35-50 kb) duplication event supporting the evolutionarily important unequal crossing-over hypothesis, but generally find independent variant occurrence and a paucity of potential de novo mutations within segments of highest homogeneity (99.1%-99.3%). Within such segments, a highly nonrandom variant clustering within adjacent higher-order repeats was found in the absence of haplotypic repeats. Such variant clusters are hardly explained by interchromosomal, fixation-driving mechanisms and likely reflect a fast, localized, intrachromosomal sequence conversion mechanism.

The mosaic structure of human pericentromeric DNA: a strategy for characterizing complex regions of the human genome

The pericentromeric regions of human chromosomes pose particular problems for both mapping and sequencing. These difficulties are due, in large part, to the presence of duplicated genomic segments that are distributed among multiple human chromosomes. To ensure contiguity of genomic sequence in these regions, we designed a sequence-based strategy to characterize different pericentromeric regions using a single (162 kb) 2p11 seed sequence as a point of reference. Molecular and cytogenetic techniques were first used to construct a paralogy map that delineated the interchromosomal distribution of duplicated segments throughout the human genome. Monochromosomal hybrid DNAs were PCR amplified by primer pairs designed to the 2p11 reference sequence. The PCR products were directly sequenced and used to develop a catalog of sequence tags for each duplicon for each chromosome. A total of 685 paralogous sequence variants were generated by sequencing 34.7 kb of paralogous pericentromeric sequence. Using PCR products as hybridization probes, we were able to identify 702 human BAC clones, of which a subset, 107 clones, were analyzed at the sequence level. We used diagnostic paralogous sequence variants to assign 65 of these BACs to at least 9 chromosomal pericentromeric regions: 1q12, 2p11, 9p11/q12, 10p11, 14q11, 15q11, 16p11, 17p11, and 22q11. Comparisons with existing sequence and physical maps for the human genome suggest that many of these BACs map to regions of the genome with sequence gaps. Our analysis indicates that large portions of pericentromeric DNA are virtually devoid of unique sequences. Instead, they consist of a mosaic of different genomic segments that have had different propensities for duplication. These biologic properties may be exploited for the rapid characterization of, not only pericentromeric DNA, but also other complex paralogous regions of the human genome.

MFASAT: a new alphoid DNA sequence isolated from Macaca fascicularis (Cercopithecidae, Primates)

A new highly repeated DNA fragment isolated from Macaca fascicularis (MFASAT) is described. Our findings obtained by sequencing, Southern blot analysis, and fluorescent in situ hybridization (FISH) on metaphasic chromosomes strongly suggest that MFASAT can be considered as a member of the alphoid DNA family characteristic of Old World monkeys. The chromosomal localization of MFASAT, obtained by FISH, showed that this alphoid DNA is present in the peri-centromeric area of all the chromosomes. MFASAT showed a high degree of conservation when compared, by sequence alignment, to other Macaca species and Papio papio as expected for species with considerable genome conservation. A low degree of homology has been found comparing M. fascicularis alphoid DNA with a more distantly related Cercopithecidae species such as Cercopithecus aethiops.

A method for evaluating phylogenetic relationship of alpha-satellite DNA suprachromosomal family by nucleotide frequency calculation

The sequence similarity among chromosome-specific alpha-satellite DNA was quantitatively evaluated by a novel procedure: nucleotide frequency calculation. Tandem-arrayed repetitive DNA segments were aligned with unit length repeat, and the nucleotide frequency at each position was used to estimate the phylogenetic distance between repetitive DNA segments. The calculations for human and chimpanzee X chromosome alpha-satellites showed that the results were consistent with the known relationships of primates, indicating that the nucleotide frequency calculation worked effectively to estimate the distances between satellite arrays. Human chromosome-specific alpha-satellites had been grouped into three suprachromosomal families (I, II, and III), and in the current work the nucleotide frequency analysis has defined the quantitative distances between the chromosome-specific alpha-satellite DNA.

Structural bistability of repetitive DNA elements featuring CA/TG dinucleotide steps and mode of evolution of satellite DNA

Satellite DNA sequences are known to be important components required for the construction of centromeres and are common to all higher eukaryotes. Nevertheless, their nucleotide sequences vary significantly, even in evolutionarily related species. In order to elucidate how the nucleotide sequences define the conformational character of centromeric satellite DNA, an evolutionary path toward repetitive units has been hypothesized. In that context, the DNA conformation of fish satellite DNA was evaluated in two ways: the organization of subrepeats and sequence characteristics were compared, and the differences in stacking energies between A-helix and B-helix and the sequence-dependent bendability of the helices were evaluated. Our findings suggest that the monomeric units making up currently observed repetitive sequences have evolved through stepwise amplification of shorter, ancestral sequences by increasing the length of the units. In addition, we suggest that potentially key sequences required for DNA amplification comprise highly flexible structures. Thus flexibility of the DNA structure may be a primary prerequisite for DNA amplification.

Novel methodology for the detection of chromosome 21-specific alpha-satellite DNA sequences

We present a novel method, based on the hybridization of allele-specific oligonucleotide probes, that allows the specific detection of chromosome 21 alpha-satellite sequences. Absence of informative polymorphic markers from the centromeric region of chromosome 21 has constituted one of the difficulties in studying the centromere of this chromosome. The alpha-satellite subfamilies from chromosomes 21 and 13 are almost identical in sequence and thus cannot be distinguished using conventional hybridization techniques. Analysis using nuclear families showed that the centromeric polymorphism, detected using our specific probe and pulsed-field gel restriction analysis, segregates in a Mendelian fashion and exhibits a high degree of polymorphism among unrelated individuals. The alphoid DNA of chromosome 21 is highly polymorphic, useful not only as a definitive anchor for the genetic map, but also for studies of chromosome 21 nondisjunction, including the unequivocal assignment of meiotic origin.

A radiation hybrid map of 48 loci including the clouston hidrotic ectodermal dysplasia locus in the pericentromeric region of chromosome 13q

To facilitate the identification of the gene responsible for Clouston hidrotic ectodermal dysplasia (HED), we used a chromosome 13-specific radiation hybrid panel to map 54 loci in the HED candidate region. The marker retention data were analyzed using RHMAP version 3. The 54 markers have an average retention frequency of 31.6% with decreasing retention as a function of distance from the centromere. Two-point analysis identified three linkage groups with a threshold lod score of 4.00; one linkage group consisted of 49 loci including the centromeric marker D13Z1 and the telomeric flanking marker for the HED candidate region D13S143. Assuming a centromeric retention model, multipoint maximum likelihood analysis of these 49 loci except D13Z1 provided a 1000:1 framework map ordering 29 loci with 21 unique map positions and approximately 2000 times more likely than the next order. Loci that could not be ordered with this level of support were positioned within a range of adjacent intervals. This map spans 347 cR9000, has an average resolution of 17.3 cR9000, and includes 3 genes (TUBA2, GJbeta2, and FGF-9), 18 ESTs, 19 polymorphic loci, and 8 single-copy DNA segments. Comparison of our RH map to a YAC contig showed an inconsistency in order involving a reversed interval of 6 loci. Fiber-FISH and FISH on interphase nuclei analyses with PACs isolated from this region supported our order. We also describe the isolation of 8 new chromosome 13q polymorphic (CA)n markers that have an average PIC value of 0.67. These data and mapping reagents will facilitate the isolation of disease genes from this region.

Site-specific retrotransposition of L1 elements within human alphoid satellite sequences

In the course of a search for microsatellites as centromeric polymorphic markers at the 3' ends of Alu or L1 elements, we observed a much higher frequency of L1 than Alu elements embedded within alpha satellite DNA. By sequence analysis of the L1 elements at their alphoid locus of insertion, we found that the insertion site was specific, with the consensus being (Py)2-10/ (Pu)3-7. All potential sites within the consensus alphoid 171-bp repeat are occupied by such elements. This confirms the finding by Feng et al. (1996; Human retrotransposon encodes a conserved endonuclease required for retrotransposition, Cell 87:905-916) that the progenitor L1 elements encode a site-specific endonuclease and that they generate copies that are inserted at these specific sites. The analysis of retrotransposed L1 elements within the alphoid domains of the acrocentric chromosomes showed that a number of loci are shared among all five acrocentrics. This sheds light on the manner in which centromeric regions of these chromosomes are exchanging information during evolution.

Padlock probes reveal single-nucleotide differences, parent of origin and in situ distribution of centromeric sequences in human chromosomes 13 and 21

Chromosome centromeres, composed of repeated DNA sequences, orchestrate the correct segregation of chromatids in cell division. We have examined the centromeres of human chromosomes 13 and 21 by studying the distribution, in situ, of two alpha satellite sequences that differ in a single nucleotide position. This was possible using padlock probes, oligo-nucleotides that can be ligated into circles upon target recognition. The segregation of individual 13 and 21 homologues in a family was followed by monitoring of the signals from two differentially labelled probes, specific for either sequence variant. A characteristic arrangement of the repeat motifs in three separate spots, oriented transverse to the length axis of the metaphase chromosomes and bilaterally symmetric, indicates that only parts of the detected regions are involved in the centromeric region, joining the sister chromatids before anaphase.

Interhomologue sequence variation of alpha satellite DNA from human chromosome 17: evidence for concerted evolution along haplotypic lineages

Alpha satellite DNA is a family of tandemly repeated DNA found at the centromeres of all primate chromosomes. Different human chromosomes 17 in the population are characterized by distinct alpha satellite haplotypes, distinguished by the presence of variant repeat forms that have precise monomeric deletions. Pair-wise comparisons of sequence diversity between variant repeat units from each haplotype show that they are closely related in sequence. Direct sequencing of PCR-amplified alpha satellite reveals heterogeneous positions between the repeat units on a chromosome as two bands at the same position on a sequencing ladder. No variation was detected in the sequence and location of these heterogeneous positions between chromosomes 17 from the same haplotype, but distinct patterns of variation were detected between chromosomes from different haplotypes. Subsequent sequence analysis of individual repeats from each haplotype confirmed the presence of extensive haplotype-specific sequence variation. Phylogenetic inference yielded a tree that suggests these chromosome 17 repeat units evolve principally along haplotypic lineages. These studies allow insight into the relative rates and/or timing of genetic turnover processes that lead to the homogenization of tandem DNA families.