Tandem repeats (TRs) are important components of eukaryotic genomes; they have both structural and functional roles: (i) they form essential chromosome structures such as centromeres and telomeres; (ii) they modify chromatin structure and affect transcription, resulting in altered gene expression and protein abundance. There are established links between variations in TRs and incompatibilities between species, evolutionary development, chromosome mis-segregation, aging, cancer outcomes and different diseases. Given the importance of TRs, it seemed essential to develop an efficient, sensitive and automated application for the identification of all kinds of TRs in various genomic sequences. Here, we present our new GRM application for identifying TRs, which is designed to overcome all the limitations of the currently existing algorithms. Our GRM algorithm provides a straightforward identification of TRs using the frequency domain but avoiding the mapping of the symbolic DNA sequence into numerical sequence, and using key string matching, but avoiding the statistical methods of locally optimizing individual key strings. Using the GRM application, we analyzed human, chimpanzee and mouse chromosome 19 genome sequences (RefSeqs), and showed that our application was very fast, efficient and simple, with a powerful graphical user interface. It can identify all types of TRs, from the smallest (2 bp) to the very large, as large as tens of kilobasepairs. It does not require any prior knowledge of sequence structure and does not require any user-defined parameters or thresholds. In this way, it ensures that a full spectrum of TRs can be detected in just one step. Furthermore, it is robust to all types of mutations in repeat copies and can identify TRs with various complexities in the sequence pattern. From this perspective, we can conclude that the GRM application is an efficient, sensitive and automated method for the identification of all kinds of TRs.

Background/Objectives: The ~1.6 kb NBPF repeat units in neuroblastoma breakpoint family (NBPF) genes are specific to humans and are associated with cognitive capacity in higher primates. While the number of NBPF monomers/Olduvai sequences in humans is approximately 2–3 times greater than in great apes, the difference in copy number values of canonical NBPF 3mer Higher-order repeats (HORs)/Olduvai triplets between humans and great apes is substantially larger. This study aims to analyze the organization and evolutionary significance of NBPF 3mer HORs/Olduvai triplets in fully sequenced primate genomes. Methods: We applied the global repeat map (GRM) algorithm to identify canonical and variant NBPF 3mer HORs/Olduvai triplets in the complete genomes of humans, chimpanzees, gorillas, and orangutans. The resulting monomer arrays were analyzed using the GRMhor algorithm to generate detailed schematic representations of NBPF HOR organization. Results: The analysis reveals a distinct difference in NBPF-related patterns among these primates, particularly in the number of tandemly organized canonical 3mer HORs/Olduvai triplets: 61 tandemly organized canonical NBPF 3mer HORs/Olduvai triplets in humans, compared to 0 in chimpanzees and orangutans, and 9 in gorillas. When considering only tandemly organized 3mer HORs/Olduvai triplets with more than three copies, the numbers adjust to 36 in humans and 0 in great apes. Furthermore, the divergence between individual NBPF monomers in humans and great apes is twice as high as that observed within great apes. Conclusions: These findings support the hypothesis that the tandem organization of NBPF 3mer HORs/Olduvai triplets plays a crucial role in enhancing cognitive capacity in humans compared to great apes, potentially providing a significant evolutionary advantage. This effect complements the impact of the increased number of individual NBPF monomers/Olduvai sequences, together contributing to a synergistic amplification effect.

From the recent genome assembly NHGRI_mPonAbe1-v2.0_NCBI (GCF_028885655.2) of orangutan chromosome 13, we computed the precise alpha satellite higher-order repeat (HOR) structure using the novel high-precision GRM2023 algorithm with Global Repeat Map (GRM) and Monomer Distance (MD) diagrams. This study rigorously identified alpha satellite HORs in the centromere of orangutan chromosome 13, discovering a novel 59mer HOR—the longest HOR unit identified in any primate to date. Additionally, it revealed the first intertwined sequence of three HORs, 18mer/27mer/45mer HORs, with a common aligned “backbone” across all HOR copies. The major 7mer HOR exhibits a Willard’s-type canonical copy, although some segments of the array display significant irregularities. In contrast, the 14mer HOR forms a regular Willard’s-type HOR array. Surprisingly, the GRM2023 high-precision analysis of chromosome 13 of human genome assembly T2T-CHM13v2.0 reveals the presence of only a 7mer HOR, despite both the orangutan and human genome assemblies being derived from whole genome shotgun sequences.

Unraveling the intricate centromere structure of human chromosomes holds profound implications, illuminating fundamental genetic mechanisms and potentially advancing our comprehension of genetic disorders and therapeutic interventions. This study rigorously identified and structurally analyzed alpha satellite higher-order repeats (HORs) within the centromere of human chromosome 15 in the complete T2T-CHM13 assembly using the high-precision GRM2023 algorithm. The most extensive alpha satellite HOR array in chromosome 15 reveals a novel cascading HOR, housing 429 15mer HOR copies, containing 4-, 7- and 11-monomer subfragments. Within each row of cascading HORs, all alpha satellite monomers are of distinct types, as in regular Willard’s HORs. However, different HOR copies within the same cascading 15mer HOR contain more than one monomer of the same type. Each canonical 15mer HOR copy comprises 15 monomers belonging to only 9 different monomer types. Notably, 65% of the 429 15mer cascading HOR copies exhibit canonical structures, while 35% display variant configurations. Identified as the second most extensive alpha satellite HOR, another novel cascading HOR within human chromosome 15 encompasses 164 20mer HOR copies, each featuring two subfragments. Moreover, a distinct pattern emerges as interspersed 25mer/26mer structures differing from regular Willard’s HORs and giving rise to a 34-monomer subfragment. Only a minor 18mer HOR array of 12 HOR copies is of the regular Willard’s type. These revelations highlight the complexity within the chromosome 15 centromeric region, accentuating deviations from anticipated highly regular patterns and hinting at profound information encoding and functional potential within the human centromere.

Background: Centromeric alpha satellite DNA is organized into higher-order repeats (HORs), whose precise structure is often difficult to resolve in standard genome assemblies. The recent telomere-to-telomere (T2T) assembly of the human genome enables complete analysis of centromeric regions, including the full structure of HOR arrays. Methods: We applied the novel high-precision GRMhor algorithm to the complete T2T-CHM13 assembly of human chromosome 21. GRMhor integrates global repeat map (GRM) and monomer distance (MD) diagrams to accurately identify, classify, and visualize HORs and their subfragments. Results: The analysis revealed a novel Cascading 11mer HOR array, in which each canonical HOR copy comprises 11 monomers belonging to 10 different monomer types. Subfragments with periodicities of 4, 7, 9, and 20 were identified within the array. A second, complex 23/25mer HOR array of mixed Willard’s/Cascading type was also detected. In contrast to the hg38 assembly, where a dominant 8mer and 33mer HOR were previously annotated, these structures were absent in the T2T-CHM13 assembly, highlighting the limitations of hg38. Notably, we discovered a novel 52mer HOR—the longest alpha satellite HOR unit reported in the human genome to date. Several subfragment repeats correspond to alphoid subfamilies previously identified using restriction enzyme digestion, but are here resolved with higher structural precision. Conclusions: Our findings demonstrate the power of GRMhor in resolving complex and previously undetected alpha satellite architectures, including the longest canonical HOR unit identified in the human genome. The precise delineation of superHORs, Cascading structures, and HOR subfragments provides unprecedented insight into the fine-scale organization of the centromeric region of chromosome 21. These results highlight both the inadequacy of earlier assemblies, such as hg38, and the critical importance of complete telomere-to-telomere assemblies for accurately characterizing centromeric DNA.