OpenNucleome for high-resolution nuclear structural and dynamical modeling.

The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome-an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of 'fixed points' within the nucleus-signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.

Circular Chromosome Conformation Capture (4C-Seq) Analysis of HBV-Host Chromosome DNA Interactome.

Hepatitis B virus (HBV) infects hepatocytes that are in the G0/G1 phase with intact nuclear membrane and organized chromosome architecture. In the nucleus of the infected cells, HBV covalently closed circular (ccc) DNA, an episomal minichromosome, serves as the template for all viral transcripts and the reservoir of persistent infection. Nuclear positioning of cccDNA can be assessed by the spatial distance between viral DNA and host chromosomal DNA through Circular Chromosome Conformation Capture (4C) combined with high-throughput sequencing (4C-seq). The 4C-seq analysis relies on proximity ligation and is commonly used for mapping genomic DNA regions that communicate within a host chromosome. The method has been tailored for studying nuclear localization of HBV episomal cccDNA in relation to the host chromosomes. In this study, we present a step-by-step protocol for 4C-seq analysis of HBV infection, including sample collection and fixation, 4C DNA library preparation, sequence library preparation, and data analysis. Although limited by proximity ligation of DNA fragments, 4C-seq analysis provides useful information of HBV localization in 3D genome, and aids the understanding of viral transcription in light of host chromatin conformation.

Homeotic DUX4 Genes Shape Dynamic Inter-Chromosomal Contacts with Nucleoli in Human Cells.

Nucleoli form interchromosomal contacts with genes controlling differentiation and carcinogenesis. DUX4 genes specify transcription factor possessing two homeodomains. Previously, using Circular Chromosome Conformation Capture (4С) approach on population of cells, it was demonstrated that DUX4 gene clusters form frequent contacts with nucleoli. It was found also that these contacts are almost completely abolished after heat shock treatment. 4C approach as all ligation-mediated methods is capable to detect rather close interactions between chromatin loops in nuclei. In order to independently confirm the formation and the frequency of the contacts in single cells we used FISH approach. Here, we show that DUX genes in single cells form stable contacts in all tested HEK293T cells. During heat shock, DUX4 genes reversibly move 1-3 µm away from the nuclei. We conclude that interchromosomal contacts formed by nucleoli are strong, dynamic, and reversible, providing both the initiation and maintenance of a differentiated state.

Validation of the Utility of the Genetically Shared Regions of Chromosomes (GD-ICS) Measuring Method in Identifying Complicated Genetic Relatedness.

BACKGROUND: Relatives share more genomic regions than unrelated individuals, with closer relatives sharing more regions. This concept, paired with the increased availability of high-throughput single nucleotide polymorphism (SNP) genotyping technologies, has made it feasible to measure the shared chromosomal regions between individuals to assess their level of relation to each other. However, such techniques have remained in the conceptual rather than practical stages in terms of applying measures or indices. Recently, we developed an index called "genetic distance-based index of chromosomal sharing (GD-ICS)" utilizing large-scale SNP data from Korean family samples and demonstrated its potential for practical applications in kinship determination. In the current study, we present validation results from various real cases demonstrating the utility of this method in resolving complex familial relationships where information obtained from traditional short tandem repeats (STRs) or lineage markers is inconclusive. METHODS: We obtained large-scale SNP data through microarray analysis from Korean individuals involving 13 kinship cases and calculated GD-ICS values using the method described in our previous study. Based on the GD-ICS reference constructed for Korean families, each disputed kinship was evaluated and validated using a combination of traditional STRs and lineage markers. RESULTS: The cases comprised those A) that were found to be inconclusive using the traditional approach, B) for which it was difficult to apply traditional testing methods, and C) that were more conclusively resolved using the GD-ICS method. This method has overcome the limitations faced by traditional STRs in kinship testing, particularly in a paternity case with STR mutational events and in confirming distant kinship where the individual of interest is unavailable for testing. It has also been demonstrated to be effective in identifying various relationships without specific presumptions and in confirming a lack of genetic relatedness between individuals. CONCLUSION: This method has been proven effective in identifying familial relationships across diverse complex and practical scenarios. It is not only useful when traditional testing methods fail to provide conclusive results, but it also enhances the resolution of challenging kinship cases, which suggests its applicability in various types of practical casework.

Chromosome Bandings and Recognition.

Chromosome banding can be defined as the lengthwise variation in staining properties along a chromosome stained with a dye. Chromosome banding became more practical in the early 1970s and is an essential technique used in karyotyping to identify human chromosomes for both clinical and research purposes. Most importantly, karyotyping is now considered a mandatory investigation of all newly diagnosed leukemias. Some banding methods, such as Giemsa (G)-, reverse (R)-, and centromere (C)-banding, still contribute greatly by being used as a routine procedure in clinical cytogenetic laboratory nowadays. Each chromosome has a unique sequence of bar code-like stripes, allowing the identification of individual homologues and the recognition of structural abnormalities through analyzing the disruption of the normal banding pattern at specific landmarks, regions, and bands as described in the ideogram. Since the quality of metaphases obtained from malignant cells is generally inferior to normal constitutional cells for karyotyping, a practical and accurate chromosome identification training guide is indispensable for a trainee or newly employed cytogenetic technologist in a cancer cytogenetic laboratory. The most common and currently used banding methods and chromosome recognition guide for distinguishable bands of each chromosome are described in detail in this chapter with an aim to facilitate quick and accurate karyotyping in cancer cells.

Profiling Chromosome Topological Features by Super-Resolution 3D Structured Illumination Microscopy.

Three-dimensional structured illumination microscopy (3D-SIM) and fluorescence in situ hybridization on three-dimensional preserved cells (3D-FISH) have proven to be robust and efficient methodologies for analyzing nuclear architecture and profiling the genome's topological features. These methods have allowed the simultaneous visualization and evaluation of several target structures at super-resolution. In this chapter, we focus on the application of 3D-SIM for the visualization of 3D-FISH preparations of chromosomes in interphase, known as Chromosome Territories (CTs). We provide a workflow and detailed guidelines for sample preparation, image acquisition, and image analysis to obtain quantitative measurements for profiling chromosome topological features. In parallel, we address a practical example of these protocols in the profiling of CTs 9 and 22 involved in the translocation t(9;22) in Chronic Myeloid Leukemia (CML). The profiling of chromosome topological features described in this chapter allowed us to characterize a large-scale topological disruption of CTs 9 and 22 that correlates directly with patients' response to treatment and as a possible potential change in the inheritance systems. These findings open new insights into how the genome structure is associated with the response to cancer treatments, highlighting the importance of microscopy in analyzing the topological features of the genome.

TERRA transcripts localize at long telomeres to regulate telomerase access to chromosome ends.

The function of TERRA in the regulation of telomerase in human cells is still debated. While TERRA interacts with telomerase, how it regulates telomerase function remains unknown. Here, we show that TERRA colocalizes with the telomerase RNA subunit hTR in the nucleoplasm and at telomeres during different phases of the cell cycle. We report that TERRA transcripts relocate away from chromosome ends during telomere lengthening, leading to a reduced number of telomeric TERRA-hTR molecules and consequent increase in "TERRA-free" telomerase molecules at telomeres. Using live-cell imaging and super-resolution microscopy, we show that upon transcription, TERRA relocates from its telomere of origin to long chromosome ends. Furthermore, TERRA depletion by antisense oligonucleotides promoted hTR localization to telomeres, leading to increased residence time and extended half-life of hTR molecules at telomeres. Overall, our findings indicate that telomeric TERRA transcripts inhibit telomere elongation by telomerase acting in trans, impairing telomerase access to telomeres that are different from their chromosome end of origin.

Amplification editing enables efficient and precise duplication of DNA from short sequence to megabase and chromosomal scale.

Duplication is a foundation of molecular evolution and a driver of genomic and complex diseases. Here, we develop a genome editing tool named Amplification Editing (AE) that enables programmable DNA duplication with precision at chromosomal scale. AE can duplicate human genomes ranging from 20 bp to 100 Mb, a size comparable to human chromosomes. AE exhibits activity across various cell types, encompassing diploid, haploid, and primary cells. AE exhibited up to 73.0% efficiency for 1 Mb and 3.4% for 100 Mb duplications, respectively. Whole-genome sequencing and deep sequencing of the junctions of edited sequences confirm the precision of duplication. AE can create chromosomal microduplications within disease-relevant regions in embryonic stem cells, indicating its potential for generating cellular and animal models. AE is a precise and efficient tool for chromosomal engineering and DNA duplication, broadening the landscape of precision genome editing from an individual genetic locus to the chromosomal scale.

Mosaicism for Autosomal Trisomies: A Comprehensive Analysis of 1266 Published Cases Focusing on Maternal Age and Reproductive History.

Mosaicism for autosomal trisomy is uncommon in clinical practice. However, despite its rarity among both prenatally and postnatally diagnoses, there are a large number of characterized and published cases. Surprisingly, in contrast to regular trisomies, no attempts at systematic analyses of mosaic carriers' demographics were undertaken. This is the first study aimed to address this gap. For that, we have screened more than eight hundred publications on mosaic trisomies, reviewing data including gender and clinical status of mosaic carriers, maternal age and reproductive history. In total, 596 publications were eligible for analysis, containing data on 948 prenatal diagnoses, including true fetal mosaicism (TFM) and confined placental mosaicism (CPM), and on 318 cases of postnatally detected mosaicism (PNM). No difference was found in maternal age between normal pregnancy outcomes with appropriate birth weight and those with intrauterine growth restriction. Unexpectedly, a higher proportion of advanced maternal ages (AMA) was found in normal outcomes compared to abnormal ones (abnormal fetus or newborn) and fetal losses, 73% vs. 56% and 50%, p = 0.0015 and p = 0.0011, correspondingly. Another intriguing finding was a higher AMA proportion in mosaic carriers with concomitant uniparental disomy (UPD) for chromosomes 7, 14, 15, and 16 compared to carriers with biparental disomy (BPD) (72% vs. 58%, 92% vs. 55%, 87% vs. 78%, and 65% vs. 24%, correspondingly); overall figures were 78% vs. 48%, p = 0.0026. Analysis of reproductive histories showed a very poor reporting but almost two-fold higher rate of mothers reporting a previous fetal loss from PNM cohort (in which almost all patients were clinically abnormal) compared to mothers from the TFM and CPM cohorts (with a large proportion of normal outcomes), 30% vs. 16%, p = 0.0072. The occurrence of a previous pregnancy with a chromosome abnormality was 1 in 13 in the prenatal cohort and 1 in 16 in the postnatal cohort, which are five-fold higher compared to published studies on non-mosaic trisomies. We consider the data obtained in this study to be preliminary despite the magnitude of the literature reviewed since reporting of detailed data was mostly poor, and therefore, the studied cohorts do not represent "big data". Nevertheless, the information obtained is useful both for clinical genetic counseling and for modeling further studies.

Reconstitution of nuclear envelope subdomain formation on mitotic chromosomes in semi-intact cells.

In metazoans, the nuclear envelope (NE) disassembles during the prophase and reassembles around segregated chromatids during the telophase. The process of NE formation has been extensively studied using live-cell imaging. At the early step of NE reassembly in human cells, specific pattern-like localization of inner nuclear membrane (INM) proteins, connected to the nuclear pore complex (NPC), was observed in the so-called "core" region and "noncore" region on telophase chromosomes, which corresponded to the "pore-free" region and the "pore-rich" region, respectively, in the early G1 interphase nucleus. We refer to these phenomena as NE subdomain formation. To biochemically investigate this process, we aimed to develop an in vitro NE reconstitution system using digitonin-permeabilized semi-intact mitotic human cells coexpressing two INM proteins, emerin and lamin B receptor, which were labeled with fluorescent proteins. The targeting and accumulation of INM proteins to chromosomes before and after anaphase onset in semi-intact cells were observed using time-lapse imaging. Our in vitro NE reconstitution system recapitulated the formation of the NE subdomain, as in living cells, although chromosome segregation and cytokinesis were not observed. This in vitro NE reconstitution required the addition of a mitotic cytosolic fraction supplemented with a cyclin-dependent kinase inhibitor and energy sources. The cytoplasmic soluble factor(s) dependency of INM protein targeting differed among the segregation states of chromosomes. Furthermore, the NE reconstituted on segregated chromosomes exhibited active nucleocytoplasmic transport competency. These results indicate that the chromosome status changes after anaphase onset for recruiting NPC components.

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs.

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

KaryoXpert: An accurate chromosome segmentation and classification framework for karyotyping analysis without training with manually labeled metaphase-image mask annotations.

Automated karyotyping is of great importance for cytogenetic research, as it speeds up the process for cytogeneticists through incorporating AI-driven automated segmentation and classification techniques. Existing frameworks confront two primary issues: Firstly the necessity for instance-level data annotation with either detection bounding boxes or semantic masks for training, and secondly, its poor robustness particularly when confronted with domain shifts. In this work, we first propose an accurate segmentation framework, namely KaryoXpert. This framework leverages the strengths of both morphology algorithms and deep learning models, allowing for efficient training that breaks the limit for the acquirement of manually labeled ground-truth mask annotations. Additionally, we present an accurate classification model based on metric learning, designed to overcome the challenges posed by inter-class similarity and batch effects. Our framework exhibits state-of-the-art performance with exceptional robustness in both chromosome segmentation and classification. The proposed KaryoXpert framework showcases its capacity for instance-level chromosome segmentation even in the absence of annotated data, offering novel insights into the research for automated chromosome segmentation. The proposed method has been successfully deployed to support clinical karyotype diagnosis.

Human telomere length is chromosome end-specific and conserved across individuals.

Short telomeres cause age-related disease, and long telomeres contribute to cancer; however, the mechanisms regulating telomere length are unclear. We developed a nanopore-based method, which we call Telomere Profiling, to determine telomere length at nearly single-nucleotide resolution. Mapping telomere reads to chromosome ends showed chromosome end-specific length distributions that could differ by more than six kilobases. Examination of telomere lengths in 147 individuals revealed that certain chromosome ends were consistently longer or shorter. The same rank order was found in newborn cord blood, suggesting that telomere length is determined at birth and that chromosome end-specific telomere length differences are maintained as telomeres shorten with age. Telomere Profiling makes precision investigation of telomere length widely accessible for laboratory, clinical, and drug discovery efforts and will allow deeper insights into telomere biology.

Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap.

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

Transcription-induced active forces suppress chromatin motion.

The organization of interphase chromosomes in a number of species is starting to emerge thanks to advances in a variety of experimental techniques. However, much less is known about the dynamics, especially in the functional states of chromatin. Some experiments have shown that the motility of individual loci in human interphase chromosome decreases during transcription and increases upon inhibiting transcription. This is a counterintuitive finding because it is thought that the active mechanical force (F) on the order of ten piconewtons, generated by RNA polymerase II (RNAPII) that is presumably transmitted to the gene-rich region of the chromatin, would render it more open, thus enhancing the mobility. We developed a minimal active copolymer model for interphase chromosomes to investigate how F affects the dynamical properties of chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of F and are enhanced at small F values, which has also been observed in experiments. In the intermediate F, the bond length between consecutive loci increases, becoming commensurate with the distance at the minimum of the attractive interaction between nonbonded loci. This results in a transient disorder-to-order transition, leading to a decreased mobility during transcription. Strikingly, the F-dependent change in the locus dynamics preserves the organization of the chromosome at [Formula: see text]. Transient ordering of the loci, which is not found in the polymers with random epigenetic profiles, in the gene-rich region might be a plausible mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

Lessons on harmonization of scoring criteria for dicentric chromosome assay in South Korea.

PURPOSE: Networking with other biodosimetry laboratories is necessary to assess the radiation exposure of many individuals in large-scale radiological accidents. The Korea biodosimetry network, K-BioDos, prepared harmonized scoring guidelines for dicentric chromosome assay to obtain homogeneous results within the network and investigated the efficiency of the guidelines. MATERIALS AND METHODS: Three laboratories in K-BioDos harmonized the scoring guidelines for dicentric chromosome assay. The results of scoring dicentric chromosomes using the harmonized scoring guidelines were compared with the laboratories' results using their own methods. Feedback was collected from the scorers following the three intercomparison exercises in 3 consecutive years. RESULTS: K-BioDos members showed comparable capacity to score dicentrics in the three exercises. However, the results of the K-BioDos guidelines showed no significant improvement over those of the scorers' own methods. According to the scorers, our harmonized guidelines led to more rejected metaphases and ultimately decreased the number of scorable metaphases compared with their own methods. Moreover, the scoring time was sometimes longer with the K-BioDos protocol because some scorers were not yet familiar with the guidelines, though most scorers reported that the time decreased or was unchanged. These challenges may cause low adherence to the guidelines. Most scorers expressed willingness to use the guidelines to select scorable metaphases or identify dicentrics for other biodosimetry works, whereas one did not want to use it due to the difference from their calibration curves. CONCLUSIONS: We identified potential resistance to following the harmonized guidelines and received requests for more detailed methods. Our findings suggest that the harmonized criteria should be continually updated, and education and training should be provided for all scorers. These changes could allow members within the biodosimetry network to successfully collaborate and support each other in large-scale radiological accidents.

Monoallelically expressed noncoding RNAs form nucleolar territories on NOR-containing chromosomes and regulate rRNA expression.

Out of the several hundred copies of rRNA genes arranged in the nucleolar organizing regions (NOR) of the five human acrocentric chromosomes, ~50% remain transcriptionally inactive. NOR-associated sequences and epigenetic modifications contribute to the differential expression of rRNAs. However, the mechanism(s) controlling the dosage of active versus inactive rRNA genes within each NOR in mammals is yet to be determined. We have discovered a family of ncRNAs, SNULs (Single NUcleolus Localized RNA), which form constrained sub-nucleolar territories on individual NORs and influence rRNA expression. Individual members of the SNULs monoallelically associate with specific NOR-containing chromosomes. SNULs share sequence similarity to pre-rRNA and localize in the sub-nucleolar compartment with pre-rRNA. Finally, SNULs control rRNA expression by influencing pre-rRNA sorting to the DFC compartment and pre-rRNA processing. Our study discovered a novel class of ncRNAs influencing rRNA expression by forming constrained nucleolar territories on individual NORs.

CBP and RAD21 Proteins Bind at the Termini of Forum Domains in Human Chromosomes.

Forum domains are 50-100-kb stretches of DNA delimited by the hotspots of double-strand breaks (DSBs). These domains possess coordinately expressed genes. However, molecular mechanisms of such regulation are not clear. It is assumed that the proteins specifically binding at the termini of domains can be involved in coordinated regulation of expression. In this study, we used the results of precise mapping of hotspots of DSBs and ChIP-Seq data for ten nuclear proteins in HEK293T cell line for a search of proteins specifically binding at forum-domain termini. We detected that two proteins, CBP and RAD24, which are known to be involved in epigenetic regulation of gene expression and formation of 3D chromosomal structures, bind at the termini. We assume that these proteins may be involved in coordinated expression of genes in forum domains.

PhaseDancer: a novel targeted assembler of segmental duplications unravels the complexity of the human chromosome 2 fusion going from 48 to 46 chromosomes in hominin evolution.

Resolving complex genomic regions rich in segmental duplications (SDs) is challenging due to the high error rate of long-read sequencing. Here, we describe a targeted approach with a novel genome assembler PhaseDancer that extends SD-rich regions of interest iteratively. We validate its robustness and efficiency using a golden-standard set of human BAC clones and in silico-generated SDs with predefined evolutionary scenarios. PhaseDancer enables extension of the incomplete complex SD-rich subtelomeric regions of Great Ape chromosomes orthologous to the human chromosome 2 (HSA2) fusion site, informing a model of HSA2 formation and unravelling the evolution of human and Great Ape genomes.

Tandem NBPF 3mer HORs (Olduvai triplets) in Neanderthal and two novel HOR tandem arrays in human chromosome 1 T2T-CHM13 assembly.

It is known that the ~ 1.6 kb Neuroblastoma BreakPoint Family (NBPF) repeats are human specific and contributing to cognitive capabilities, with increasing frequency in higher order repeat 3mer HORs (Olduvai triplets). From chimpanzee to modern human there is a discontinuous jump from 0 to ~ 50 tandemly organized 3mer HORs. Here we investigate the structure of NBPF 3mer HORs in the Neanderthal genome assembly of Pääbo et al., comparing it to the results obtained for human hg38.p14 chromosome 1. Our findings reveal corresponding NBPF 3mer HOR arrays in Neanderthals with slightly different monomer structures and numbers of HOR copies compared to humans. Additionally, we compute the NBPF 3mer HOR pattern for the complete telomere-to-telomere human genome assembly (T2T-CHM13) by Miga et al., identifying two novel tandem arrays of NBPF 3mer HOR repeats with 5 and 9 NBPF 3mer HOR copies. We hypothesize that these arrays correspond to novel NBPF genes (here referred to as NBPFA1 and NBPFA2). Further improving the quality of the Neanderthal genome using T2T-CHM13 as a reference would be of great interest in determining the presence of such distant novel NBPF genes in the Neanderthal genome and enhancing our understanding of human evolution.