Repeat Array Research Articles

ZBTB24 is a conserved multifaceted transcription factor at genes and centromeres that governs the DNA methylation state and expression of satellite repeats.

Since its discovery as a causative gene of the Immunodeficiency with Centromeric instability and Facial anomalies syndrome, ZBTB24 has emerged as a key player in DNA methylation, immunity and development. By extensively analyzing ZBTB24 genomic functions in ICF-relevant mouse and human cellular models, we document here its multiple facets as a transcription factor, with key roles in immune response-related genes expression and also in early embryonic development. Using a constitutive Zbtb24 ICF-like mutant and an auxin-inducible degron system in mouse embryonic stem cells, we showed that ZBTB24 is recruited to centromeric satellite DNA where it is required to establish and maintain the correct DNA methylation patterns through the recruitment of DNMT3B. The ability of ZBTB24 to occupy centromeric satellite DNA is conserved in human cells. Together, our results unveiled an essential and underappreciated role for ZBTB24 at mouse and human centromeric satellite repeat arrays by controlling their DNA methylation and transcription status.

The first chicken oocyte nucleus whole transcriptomic profile defines the spectrum of maternal mRNA and non-coding RNA genes transcribed by the lampbrush chromosomes.

Lampbrush chromosomes, with their unusually high rate of nascent RNA synthesis, provide a valuable model for studying mechanisms of global transcriptome up-regulation. Here, we obtained a whole-genomic profile of transcription along the entire length of all lampbrush chromosomes in the chicken karyotype. With nuclear RNA-seq, we obtained information about a wider set of transcripts, including long non-coding RNAs retained in the nucleus and stable intronic sequence RNAs. For a number of protein-coding genes, we visualized their nascent transcripts on the lateral loops of lampbrush chromosomes by RNA-FISH. The set of genes transcribed on the lampbrush chromosomes is required for basic cellular processes and is characterized by a broad expression pattern. We also present the first high-throughput transcriptome characterization of miRNAs and piRNAs in chicken oocytes at the lampbrush chromosome stage. Major targets of predicted piRNAs include CR1 and long terminal repeat (LTR) containing retrotransposable elements. Transcription of tandem repeat arrays was demonstrated by alignment against the whole telomere-to-telomere chromosome assemblies. We show that transcription of telomere-derived RNAs is initiated at adjacent LTR elements. We conclude that hypertranscription on the lateral loops of giant lampbrush chromosomes is required for synthesizing large amounts of transferred to the embryo maternal RNA for thousands of genes.

Accurate prenatal diagnosis of facioscapulohumeral muscular dystrophy 1 using nanopore sequencing

BackgroundFacioscapulohumeral muscular dystrophy 1 (FSHD1) is an autosomal dominant muscular disorder mainly caused by the contraction and hypomethylation of the D4Z4 repeat array in chromosome 4q35. Prenatal diagnosis of FSHD1...

Human Satellite 3 DNA encodes megabase-scale transcription factor binding platforms.

Eukaryotic genomes are frequently littered with large arrays of tandem repeats, called satellite DNA, which underlie the constitutive heterochromatin often found around centromeric regions. While some satellite DNA types have well-established roles in centromere biology, other abundant satellite DNAs have poorly characterized functions. For example, Human Satellite 3 (HSat3), which makes up roughly 2% of the human genome, forms enormous arrays up to tens of megabases, but these arrays play no known roles in centromere function and were almost fully excluded from genome assemblies until recently. As a result, these massive genomic regions have remained relatively understudied, and the potential functional roles of HSat3 have remained largely unknown. To address this, we performed a systematic screen for novel HSat3 binding factors. Our work revealed HSat3 arrays contain high densities of transcription factor (TF) motifs that are bound by factors related to multiple, highly conserved signaling pathways. Unexpectedly, the most enriched TFs in HSat3 belong to the Hippo pathway transcription effector family TEAD. We found that TEAD recruits the co-activator YAP to HSat3 regions in a cell-state specific manner. Using RNA polymerase-I reporter assays, targeted repression of HSat3, inducible degradation of YAP, and super-resolution microscopy, we show that HSat3 arrays can localize YAP/TEAD inside the nucleolus, where YAP regulates RNA Polymerase-I activity. Beyond revealing a direct relationship between the Hippo pathway and ribosomal DNA regulation, this work demonstrates that satellite DNA can encode multiple transcription factor binding motifs, defining a new role for these enormous genomic elements.

Phage-triggered reverse transcription assembles a toxic repetitive gene from a noncoding RNA.

Reverse transcription has frequently been co-opted for cellular functions and in prokaryotes is associated with protection against viral infection, but the underlying mechanisms of defense are generally unknown. Here, we show that in the DRT2 defense system, the reverse transcriptase binds a neighboring pseudoknotted noncoding RNA. Upon bacteriophage infection, a template region of this RNA is reverse transcribed into an array of tandem repeats that reconstitute a promoter and open reading frame, allowing expression of a toxic repetitive protein and an abortive infection response. Biochemical reconstitution of this activity and cryo-electron microscopy provide a molecular basis for repeat synthesis. Gene synthesis from a noncoding RNA is a previously unknown mode of genetic regulation in prokaryotes.

Complex genetic variation in nearly complete human genomes.

Diverse sets of complete human genomes are required to construct a pangenome reference and to understand the extent of complex structural variation. Here, we sequence 65 diverse human genomes and build 130 haplotype-resolved assemblies (130 Mbp median continuity), closing 92% of all previous assembly gaps1,2 and reaching telomere-to-telomere (T2T) status for 39% of the chromosomes. We highlight complete sequence continuity of complex loci, including the major histocompatibility complex (MHC), SMN1/SMN2, NBPF8, and AMY1/AMY2, and fully resolve 1,852 complex structural variants (SVs). In addition, we completely assemble and validate 1,246 human centromeres. We find up to 30-fold variation in α-satellite high-order repeat (HOR) array length and characterize the pattern of mobile element insertions into α-satellite HOR arrays. While most centromeres predict a single site of kinetochore attachment, epigenetic analysis suggests the presence of two hypomethylated regions for 7% of centromeres. Combining our data with the draft pangenome reference1 significantly enhances genotyping accuracy from short-read data, enabling whole-genome inference3 to a median quality value (QV) of 45. Using this approach, 26,115 SVs per sample are detected, substantially increasing the number of SVs now amenable to downstream disease association studies.

Atlas of telomeric repeat diversity in Arabidopsis thaliana

BackgroundTelomeric repeat arrays at the ends of chromosomes are highly dynamic in composition, but their repetitive nature and technological limitations have made it difficult to assess their true variation in genome diversity surveys.ResultsWe have comprehensively characterized the sequence variation immediately adjacent to the canonical telomeric repeat arrays at the very ends of chromosomes in 74 genetically diverse Arabidopsis thaliana accessions. We first describe several types of distinct telomeric repeat units and then identify evolutionary processes such as local homogenization and higher-order repeat formation that shape diversity of chromosome ends. By comparing largely isogenic samples, we also determine repeat number variation of the degenerate and variant telomeric repeat array at both the germline and somatic levels. Finally, our analysis of haplotype structure uncovers chromosome end-specific patterns in the distribution of variant telomeric repeats, and their linkage to the more proximal non-coding region.ConclusionsOur findings illustrate the spectrum of telomeric repeat variation at multiple levels in A. thaliana—in germline and soma, across all chromosome ends, and across genetic groups—thereby expanding our knowledge of the evolution of chromosome ends.

D4Z4 Hypomethylation in Human Germ Cells.

Expression of the double homeobox 4 (DUX4) transcription factor is highly regulated in early embryogenesis and is subsequently epigenetically silenced. Ectopic expression of DUX4 due to hypomethylation of the D4Z4 repeat array on permissive chromosome 4q35 alleles is associated with facioscapulohumeral muscular dystrophy (FSHD). In peripheral blood samples from 188 healthy individuals, D4Z4 methylation was highly variable, ranging from 19% to 76%, and was not affected by age. In 48 FSHD2 patients, D4Z4 methylation varied from 3% to 30%. Given that DUX4 is one of the earliest transcribed genes after fertilization, the D4Z4 array is expected to be unmethylated in mature germ cells. Deep bisulfite sequencing of 188 mainly normozoospermic sperm samples revealed an average methylation of 2.5% (range 0.3-22%). Overall, the vast majority (78%) of individual sperm cells displayed no methylation at all. In contrast, only 19 (17.5%) of 109 individual germinal vesicle (GV) oocytes displayed D4Z4 methylation <2.5%. However, it is not unexpected that immature GV oocytes which are not usable for assisted reproduction are endowed with D4Z4 (up to 74%) hypermethylation and/or abnormal (PEG3 and GTL2) imprints. Although not significant, it is interesting to note that the pregnancy rate after assisted reproduction was higher for donors of sperm samples and oocytes with <2.5% methylation.

Origin and maintenance of large ribosomal RNA gene repeat size in mammals.

The genes encoding ribosomal RNA are highly conserved across life and in almost all eukaryotes are present in large tandem repeat arrays called the rDNA. rDNA repeat unit size is conserved across most eukaryotes but has expanded dramatically in mammals, principally through the expansion of the intergenic spacer region that separates adjacent rRNA coding regions. Here, we used long-read sequence data from representatives of the major amniote lineages to determine where in amniote evolution rDNA unit size increased. We find that amniote rDNA unit sizes fall into two narrow size classes: "normal" (∼11-20 kb) in all amniotes except monotreme, marsupial, and eutherian mammals, which have "large" (∼35-45 kb) sizes. We confirm that increases in intergenic spacer length explain much of this mammalian size increase. However, in stark contrast to the uniformity of mammalian rDNA unit size, mammalian intergenic spacers differ greatly in sequence. These results suggest a large increase in intergenic spacer size occurred in a mammalian ancestor and has been maintained despite substantial sequence changes over the course of mammalian evolution. This points to a previously unrecognized constraint on the length of the intergenic spacer, a region that was thought to be largely neutral. We finish by speculating on possible causes of this constraint.

Clinical and genetic evaluation of hereditary myopathies in an adult Saudi cohort

BackgroundDiagnosis of hereditary myopathy is often challenging owing to overlapping clinical phenotypes and muscle histopathological findings. This retrospective study aimed to identify the phenotypic and genotypic spectra of hereditary myopathies at a tertiary hospital in Riyadh, Saudi Arabia.MethodsWe reviewed the medical records of patients with hereditary myopathy who were evaluated between January 2018 and December 2022.ResultsEighty-seven patients (78 families) were included, two-thirds were men with a mean age of 35 (SD 14.2) years. Limb-girdle muscular dystrophy (LGMD) was the most prevalent clinical diagnosis (25 cases; 29%), of whom, a genetic diagnosis was achieved in 15 of 22 patients tested (68%). In genetically confirmed LGMD, the most prevalent disorders were dysferlinopathy (27%) followed by fukutin-related protein (FKRP) - related limb girdle muscular dystrophy (20%), sarcoglycanopathy (20%), lamin A/C related myopathy (13%), and calpain-3 myopathy (13%). In 26 patients with pathogenic/likely pathogenic variants, the genetic testing method was whole exome sequencing (WES) (42%), Next generation sequencing (NGS) (31%), and targeted single gene analysis (27%). The sensitivity of each genetic testing method was as follows: 100% for targeted single-gene analysis, 100% for targeted analysis of D4Z4 repeat array units, 88% for myotonic dystrophy protein kinase (DMPK) repeat expansion analysis, 42% for NGS-neuromuscular panel, and 46% for WES.ConclusionThe prevalent types of hereditary myopathies were consistent with those reported locally and internationally. This study highlights the diagnostic yield of various molecular genetic tests for the diagnosis of hereditary myopathy in an adult cohort and the need for improved access to advanced molecular testing in cases suspected to have facioscapulohumeral muscular dystrophy (FSHD) or mitochondrial myopathies.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BackgroundCentromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.ResultsBased on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.ConclusionOur data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

ModDotPlot-rapid and interactive visualization of tandem repeats.

A common method for analyzing genomic repeats is to produce a sequence similarity matrix visualized via a dot plot. Innovative approaches such as StainedGlass have improved upon this classic visualization by rendering dot plots as a heatmap of sequence identity, enabling researchers to better visualize multi-megabase tandem repeat arrays within centromeres and other heterochromatic regions of the genome. However, computing the similarity estimates for heatmaps requires high computational overhead and can suffer from decreasing accuracy. In this work, we introduce ModDotPlot, an interactive and alignment-free dot plot viewer. By approximating average nucleotide identity via a k-mer-based containment index, ModDotPlot produces accurate plots orders of magnitude faster than StainedGlass. We accomplish this through the use of a hierarchical modimizer scheme that can visualize the full 128 Mb genome of Arabidopsis thaliana in under 5 min on a laptop. ModDotPlot is bundled with a graphical user interface supporting real-time interactive navigation of entire chromosomes. ModDotPlot is available at https://github.com/marbl/ModDotPlot.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

BackgroundThe honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses.ResultsPacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants.ConclusionsOur new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Analysis and occurrence of biallelic pathogenic repeat expansions in RFC1 in a German cohort of patients with a main clinical phenotype of motor neuron disease

Biallelic pathogenic repeat expansions in RFC1 were recently identified as molecular origin of cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) as well as of one of the most common causes of adult-onset ataxia. In the meantime, the phenotypic spectrum has expanded massively and now includes mimics of multiple system atrophy or parkinsonism. After identifying a patient with a clinical diagnosis of amyotrophic lateral sclerosis (ALS) as a carrier of biallelic pathogenic repeat expansions in RFC1, we studied a cohort of 106 additional patients with a clinical main phenotype of motor neuron disease (MND) to analyze whether such repeat expansions are more common in MND patients. Indeed, two additional MND patients (one also with ALS and one with primary lateral sclerosis/PLS) have been identified as carrier of biallelic pathogenic repeat expansions in RFC1 in the absence of another genetic alteration explaining the phenotype, suggesting motor neuron disease as another extreme phenotype of RFC1 spectrum disorder. Therefore, MND might belong to the expanding phenotypic spectrum of pathogenic RFC1 repeat expansions, particularly in those MND patients with additional features such as sensory and/or autonomic neuropathy, vestibular deficits, or cerebellar signs. By systematically analyzing the RFC1 repeat array using Oxford nanopore technology long-read sequencing, our study highlights the high intra- and interallelic heterogeneity of this locus and allows the identification of the novel repeat motif ‘ACAAG’.

Balancing act: BRCA2's elaborate management of telomere replication through control of G-quadruplex dynamicity.

In billion years of evolution, eukaryotes preserved the chromosome ends with arrays of guanine repeats surrounded by thymines and adenines, which can form stacks of four-stranded planar structure known as G-quadruplex (G4). The rationale behind the evolutionary conservation of the G4 structure at the telomere remained elusive. Our recent study has shed light on this matter by revealing that telomere G4 undergoes oscillation between at least two distinct folded conformations. Additionally, tumor suppressor BRCA2 exhibits a unique mode of interaction with telomere G4. To elaborate, BRCA2 directly interacts with G-triplex (G3)-derived intermediates that form during the interconversion of the two different G4 states. In doing so, BRCA2 remodels the G4, facilitating the restart of stalled replication forks. In this review, we succinctly summarize the findings regarding the dynamicity of telomeric G4, emphasize its importance in maintaining telomere replication homeostasis, and the physiological consequences of losing G4 dynamicity at the telomere.

Impact of viral telomeric repeat sequences on herpesvirus vector vaccine integration and persistence.

Marek's disease virus (MDV) vaccines were the first vaccines that protected against cancer. The avirulent turkey herpesvirus (HVT) was widely employed and protected billions of chickens from a deadly MDV infection. It is also among the most common vaccine vectors providing protection against a plethora of pathogens. HVT establishes latency in T-cells, allowing the vaccine virus to persist in the host for life. Intriguingly, the HVT genome contains telomeric repeat arrays (TMRs) at both ends; however, their role in the HVT life cycle remains elusive. We have previously shown that similar TMRs in the MDV genome facilitate its integration into host telomeres, which ensures efficient maintenance of the virus genome during latency and tumorigenesis. In this study, we investigated the role of the TMRs in HVT genome integration, latency, and reactivation in vitro and in vivo. Additionally, we examined HVT infection of feather follicles. We generated an HVT mutant lacking both TMRs (vΔTMR) that efficiently replicated in cell culture. We could demonstrate that wild type HVT integrates at the ends of chromosomes containing the telomeres in T-cells, while integration was severely impaired in the absence of the TMRs. To assess the role of TMRs in vivo, we infected one-day-old chickens with HVT or vΔTMR. vΔTMR loads were significantly reduced in the blood and hardly any virus was transported to the feather follicle epithelium where the virus is commonly shed. Strikingly, latency in the spleen and reactivation of the virus were severely impaired in the absence of the TMRs, indicating that the TMRs are crucial for the establishment of latency and reactivation of HVT. Our findings revealed that the TMRs facilitate integration of the HVT genome into host chromosomes, which ensures efficient persistence in the host, reactivation, and transport of the virus to the skin.

Behavioral signatures of the rapid recruitment of long-term memory to overcome working memory capacity limits.

Working- and long-term memory are often studied in isolation. To better understand the specific limitations of working memory, effort is made to reduce the potential influence of long-term memory on performance in working memory tasks (e.g., asking participants to remember artificial, abstract items rather than familiar real-world objects). However, in everyday life we use working- and long-term memory in tandem. Here, our goal was to characterize how long-term memory can be recruited to circumvent capacity limits in a typical visual working memory task (i.e., remembering colored squares). Prior work has shown that incidental repetitions of working memory arrays often do not improve visual working memory performance - even after dozens of incidental repetitions, working memory performance often shows no improvement for repeated arrays. Here, we used a whole-report working memory task with explicit rather than incidental repetitions of arrays. In contrast to prior work with incidental repetitions, in two behavioral experiments we found that explicit repetitions of arrays yielded robust improvement to working memory performance, even after a single repetition. Participants performed above chance at recognizing repeated arrays in a later long-term memory test, consistent with the idea that long-term memory was used to rapidly improve performance across array repetitions. Finally, we analyzed inter-item response times and we found a response time signature of chunk formation that only emerged after the array was repeated (inter-response time slowing after two to three items); thus, inter-item response times may be useful for examining the coordinated interaction of visual working and long-term memory in future work.

Fabrication Process Control to Realize High Yield, Uniform, Repeatable Low-Frequency Detector Arrays for the LiteBIRD CMB Experiment

Fabrication Process Control to Realize High Yield, Uniform, Repeatable Low-Frequency Detector Arrays for the LiteBIRD CMB Experiment

Exchange of subtelomeric regions between chromosomes 4q and 10q reverts the FSHD genotype and phenotype.

The most common form of facioscapulohumeral dystrophy (FSHD1) is caused by a partial loss of the D4Z4 macrosatellite repeat array in the subtelomeric region of chromosome 4. Patients with FSHD1 typically carry 1 to 10 D4Z4 repeats, whereas nonaffected individuals have 11 to 150 repeats. The ~150-kilobyte subtelomeric region of the chromosome 10q exhibits a ~99% sequence identity to the 4q, including the D4Z4 array. Nevertheless, contractions of the chr10 array do not cause FSHD or any known disease, as in most people D4Z4 array on chr10 is flanked by the nonfunctional polyadenylation signal, not permitting the DUX4 expression. Here, we attempted to correct the FSHD genotype by a CRISPR-Cas9-induced exchange of the chr4 and chr10 subtelomeric regions. We demonstrated that the induced t(4;10) translocation can generate recombinant genotypes translated into improved FSHD phenotype. FSHD myoblasts with the t(4;10) exhibited reduced expression of the DUX4 targets, restored PAX7 target expression, reduced sensitivity to oxidative stress, and improved differentiation capacity.

A bioinformatic approach to identify confirmed and probable CRISPR-Cas systems in the Acinetobacter calcoaceticus-Acinetobacter baumannii complex genomes.

The Acinetobacter calcoaceticus-Acinetobacter baumannii complex, or Acb complex, consists of six species: Acinetobacter baumannii, Acinetobacter calcoaceticus, Acinetobacter nosocomialis, Acinetobacter pittii, Acinetobacter seifertii, and Acinetobacter lactucae. A. baumannii is the most clinically significant of these species and is frequently related to healthcare-associated infections (HCAIs). Clustered regularly interspaced short palindromic repeat (CRISPR) arrays and associated genes (cas) constitute bacterial adaptive immune systems and function as variable genetic elements. This study aimed to conduct a genomic analysis of Acb complex genomes available in databases to describe and characterize CRISPR systems and cas genes. Acb complex genomes available in the NCBI and BV-BRC databases, the identification and characterization of CRISPR-Cas systems were performed using CRISPRCasFinder, CRISPRminer, and CRISPRDetect. Sequence types (STs) were determined using the Oxford scheme and ribosomal multilocus sequence typing (rMLST). Prophages were identified using PHASTER and Prophage Hunter. A total of 293 genomes representing six Acb species exhibited CRISPR-related sequences. These genomes originate from various sources, including clinical specimens, animals, medical devices, and environmental samples. Sequence typing identified 145 ribosomal multilocus sequence types (rSTs). CRISPR-Cas systems were confirmed in 26.3% of the genomes, classified as subtypes I-Fa, I-Fb and I-Fv. Probable CRISPR arrays and cas genes associated with CRISPR-Cas subtypes III-A, I-B, and III-B were also detected. Some of the CRISPR-Cas systems are associated with genomic regions related to Cap4 proteins, and toxin-antitoxin systems. Moreover, prophage sequences were prevalent in 68.9% of the genomes. Analysis revealed a connection between these prophages and CRISPR-Cas systems, indicating an ongoing arms race between the bacteria and their bacteriophages. Furthermore, proteins associated with anti-CRISPR systems, such as AcrF11 and AcrF7, were identified in the A. baumannii and A. pittii genomes. This study elucidates CRISPR-Cas systems and defense mechanisms within the Acb complex, highlighting their diverse distribution and interactions with prophages and other genetic elements. This study also provides valuable insights into the evolution and adaptation of these microorganisms in various environments and clinical settings.