In this article, we show that mono/oligonucleotide quadruplets, as basic structures of DNA, along with our classification of trinucleotides, disclose an organization of genomes based on purine–pyrimidine symmetry. Moreover, the structure and stability of DNA are influenced by the Watson–Crick pairing and the natural law of DNA creation and conservation, according to which the same mono- or oligonucleotide insertion must be inserted simultaneously into both strands of DNA. Taken together, they lead to quadruplets with central mirror symmetry and bidirectional DNA strand orientation and are incorporated into Chargaff’s second parity rule (CSPR). Performing our quadruplet frequency analysis of all human chromosomes and of Neuroblastoma BreakPoint Family (NBPF) genes, which code Olduvai protein domains in the human genome, we show that the coding part of DNA violates CSPR. This may shed new light and give rise to a novel hypothesis on DNA creation and its evolution. In this framework, the logarithmic relationship between oligonucleotide order and minimal DNA sequence length, to establish the validity of CSPR, automatically follows from the quadruplet structure of the genomic sequence. The problem of the violation of CSPR in rare symbionts is discussed.

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.

Recent analyses have revealed that orangutan alpha satellite higher-order repeat (HOR) arrays in complete centromeres are composed of three to four distinct HOR blocks, each sharing only 80–90% sequence identity, thus forming a patchwork-quilt pattern of independent HOR expansions. In contrast, using our novel HOR-detection algorithm GRhor, we analyzed the complete Y chromosome centromere in orangutan and identified a highly ordered and complex alpha satellite 58mer superHOR array, comprising 67 HOR copies, including 46 highly identical canonical copies with a remarkably low divergence of only 0.25%. Given that the largest known human alpha satellite HOR is the 34mer on the Y chromosome, this novel 58mer structure qualifies as a superHOR. The canonical 58mer HOR contains only 44 distinct monomer types, with 14 types repeated within the unit, resulting in a unique five-row cascading organization. Such complexity is not detectable using standard HOR-searching tools employed in previous studies. Additionally, we identified a second, less pronounced 45mer cascading superHOR array with 0.81% divergence. For comparative purposes, we also detected a cascading 18mer HOR in gorilla and a Willard-type 28mer HOR in chimpanzee Y centromeres. Notably, preliminary genome-wide analysis in orangutan reveals other superHORs, including 84mer and 53mer arrays in chromosome 5; a 54mer in chromosome 10; a 51mer in chromosome 14; a 53mer in chromosome 15; and a 45mer in chromosome 22. These findings underscore the power of GRMhor in revealing highly structured and species-specific HOR architectures, with potential implications for centromere evolution and primate comparative genomics.