Repetitive DNA Sequences Research Articles

ULTRA-Effective labeling of tandem repeats in genomic sequence

Abstract In the age of long read sequencing, genomics researchers now have access to accurate repetitive DNA sequence (including satellites) that, due to the limitations of short read sequencing, could previously be observed only as unmappable fragments. Tools that annotate repetitive sequence are now more important than ever, so that we can better understand newly uncovered repetitive sequences, and also so that we can mitigate errors in bioinformatic software caused by those repetitive sequences. To that end, we introduce the 1.0 release of our tool for identifying and annotating locally-repetitive sequence, ULTRA (ULTRA Locates Tandemly Repetitive Areas). ULTRA is fast enough to use as part of an efficient annotation pipeline, produces state-of-the-art reliable coverage of repetitive regions containing many mutations, and provides interpretable statistics and labels for repetitive regions. It is released under an open license, and is available for download at https://github.com/TravisWheelerLab/ULTRA

Novel centromeric repetitive DNA elements reveal karyotype dynamics in polyploid sainfoin (Onobrychis viciifolia).

Polyploidy is a common feature in eukaryotes with one of paramount consequences leading to better environmental adaptation. Heterochromatin is often located at telomeres and centromeres and contains repetitive DNA sequences. Sainfoin (Onobrychis viciifolia) is an important perennial forage legume for sustainable agriculture. However, there are only a few studies on the sainfoin genome and chromosomes. In this study, novel tandem repetitive DNA sequences of the sainfoin genome (OnVi180, OnVi169, OnVi176 and OnVidimer) were characterized using bioinformatics, molecular and cytogenetic approaches. The OnVi180 and OnVi169 elements colocalized within functional centromeres. The OnVi176 and OnVidimer elements were localized in centromeric, subtelomeric and interstitial regions. We constructed a sainfoin karyotype that distinguishes the seven basic chromosome groups. Our study provides the first detailed description of heterochromatin and chromosome structure of sainfoin and proposes an origin of heterozygous ancestral genomes, possibly from the same ancestral diploid species, not necessarily from different species, or for chromosome rearrangements after polyploidy. Overall, we discuss our novel and complementary findings in a polyploid crop with unique and complex chromosomal features.

Chromosome-level genome assembly of the morabine grasshopper Vandiemenella viatica19

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as “old” and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.

Centromere drive may propel the evolution of chromosome and genome size in plants.

Genome size is influenced by natural selection and genetic drift acting on variations from polyploidy and repetitive DNA sequences. We hypothesized that centromere drive, where centromeres compete for inclusion in the functional gamete during meiosis, may also affect genome and chromosome size. This competition occurs in asymmetric meiosis, where only one of the four meiotic products becomes a gamete. If centromere drive influences chromosome size evolution, it may also impact post-polyploid diploidization, where a polyploid genome is restructured to function more like a diploid through chromosomal rearrangements, including fusions. We tested if plant lineages with asymmetric meiosis exhibit faster chromosome size evolution compared to those with only symmetric meiosis, which lack centromere drive as all four meiotic products become gametes. We also examined if positive selection on centromeric histone H3 (CENH3), a protein that can suppress centromere drive, is more frequent in these asymmetric lineages. We analyzed plant groups with different meiotic modes: asymmetric in gymnosperms and angiosperms, and symmetric in bryophytes, lycophytes, and ferns. We selected species based on available CENH3 gene sequences and chromosome size data. Using Ornstein-Uhlenbeck evolutionary models and phylogenetic regressions, we assessed the rates of chromosome size evolution and the frequency of positive selection on CENH3 in these clades. Our analyses showed that clades with asymmetric meiosis have a higher frequency of positive selection on CENH3 and increased rates of chromosome size evolution compared to symmetric clades. Our findings support the hypothesis that centromere drive accelerates chromosome and genome size evolution, potenatially also influencing the process of post-polyploid diploidization. We propose a model which in a single famework helps explain the stability of chromosome size in symmetric lineages (bryophytes, lycophytes, and ferns) and its variability in asymmetric lineages (gymnosperms and angiosperms), providing a foundation for future research in plant genome evolution.

Evolutionary Dynamics of Satellite DNA Repeats across the Tettigoniidae Family: Insights from Genomic Analysis.

Satellite DNA repeats are repetitive DNA sequences found in eukaryotic genomes, typically consisting of short DNA motifs repeated in tandem arrays. Despite the vast body of literature on satellite DNA repeats in other taxa, investigations specifically targeting Tettigoniidae remain conspicuously absent. Our study aims to fill a critical gap in our understanding of satellitome evolutionary processes shaping Tettigoniidae genomes. Repeatome analysis revealed that the Meconema thalassinum genome comprises 92%, and Phryganogryllacris superangulata had the lowest value of 34%, with an average of 67% in other Tettigoniidae species. The analysis reveals significant variation in the number of satellite DNA repeats across species of the Tettigoniidae family, with M. thalassinum exhibiting the highest count, 246, reported in insects to date and the lowest count, 10, in Pholidoptera griseoptera. Ruspolia dubia and Ruspolia yunnana, which are congeneric species, showcase distinct counts of 104 and 84 families, respectively. Satellite DNA repeats in R. dubia exhibit the highest abundance, constituting 17.2% of the total genome, while the lowest abundance was reported in P. griseoptera, at 5.65%. The genome size correlates weakly with the satellite DNA family count (rs = 0.42, p = 0.29), but a strong correlation exists between satellite abundance and family number (rs = 0.73, p = 0.03). Moreover, the analysis of satellite DNA gain and loss patterns provides insights into the amplification and homogenization of satellite DNA families within the genome, with species-specific repeats exhibiting a positive trend toward amplification. The chromosomal distribution in M. thalassinum displayed that the highest accumulation was observed on Chr12, Chr01, and Chr04, constituting 17.79%, 17.4%, and 17.22% of the total chromosome size, respectively. The chromosome-specific propagation of satellite DNA families was evident, with MthSat01 solely on chromosome 1 and MthSat170 on chromosome 2, sharing 1.64% and 2.33%. The observed conservation and variations in satellite DNA number and abundances, along with distinct patterns of gain and loss, indicate the influence of potentially diverse evolutionary processes shaping the genomic landscape of these insects, which requires further investigation. Furthermore, the differential accumulation of satellite DNA on specific chromosomes implies that potential chromosome-specific functions or structural features influence the retention and proliferation of satellite sequences.

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles.

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Seamless site-directed mutagenesis in complex cloned DNA sequences using the RedEx method.

Seamless site-directed mutagenesis is an important technique for studying protein functions, tuning enzyme catalytic activities and modifying genetic elements in multiple rounds because it can insert, delete or substitute nucleotides, DNA segments or even entire genes at the target site without introducing any unwanted change. To facilitate seamless site-directed mutagenesis in large plasmids and bacterial artificial chromosomes (BACs) with repetitive sequences, we recently developed the RedEx strategy. Compared with previous methods, our approach achieves the recovery of correct recombinants with high accuracy by circumventing unwanted recombination between repetitive sequences. RedEx readily yields more than 80% accuracy in seamless DNA insertion and deletion in large multimodular polyketide synthase gene clusters, which are among the most difficult targets due to the large number of repetitive DNA sequences in modules encoding almost identical enzymes. Here we present the RedEx method by describing in detail the seamless site-directed mutagenesis in a BAC vector. Overall, the process includes three parts: (1) insertion of the RedEx cassette containing the desired mutation together with selection-counterselection markers flanked by unique restriction sites and 20-bp overlapping sequences into the target site by recombineering, (2) removal of the selection-counterselection markers in the BAC by restriction digestion and (3) circularization of the linear BAC by exonuclease-mediated in vitro DNA annealing. This protocol can be performed within 3 weeks and will enable researchers with DNA cloning experience to master seamless site-directed mutagenesis to accelerate their research.

GenRepAI: Utilizing Artificial Intelligence to Identify Repeats in Genomic Suffix Trees

Background: The human genome is densely populated with repetitive DNA sequences that play crucial roles in genomic functions and structures but are also implicated in over 40 human diseases. The computational challenge of identifying and characterizing these repeats is significant due to the complexity and size of the genome, which are overwhelming traditional algorithms. Methods: To address these challenges, we propose GenRepAI, a deep learning framework to navigate and analyze genomic suffix trees. GenRepAI employs supervised machine learning classifiers trained on labeled datasets of repeat annotations and unsupervised anomaly detection to identify novel repeat sequences. The models are trained using convolutional neural networks (CNNs), long short-term memory networks (LSTMs), and vision transformers to classify and annotate repeats within the human genome. Results: GenRepAI is designed to comprehensively profile repeats that underlie various neurological diseases, allowing researchers to identify pathogenic expansions. The framework will integrate into existing genomic analysis pipelines, with the capability to screen patient genomes and highlight potential causal variants for further validation. Conclusion: GenRepAI is set to become a foundational tool in genomics, leveraging artificial intelligence to enhance the characterization of repetitive sequences. It promises significant advancements in the molecular diagnosis of repeat expansion disorders and contributes to a deeper understanding of genomic structure and function, with broad applications in personalized medicine.

Short tandem repeat expansions in cortical layer-specific genes implicate in phenotypic severity and adaptability of autism spectrum disorder.

Short tandem repeats (STRs) are repetitive DNA sequences and highly mutable in various human disorders. While the involvement of STRs in various genetic disorders has been extensively studied, their role in autism spectrum disorder (ASD) remains largely unexplored. In this study, we aimed to investigate genetic association of STR expansions with ASD using whole genome sequencing (WGS) and identify risk loci associated with ASD phenotypes. We analyzed WGS data of 634 ASD families and performed genome-wide evaluation for 12,929 STR loci. We found rare STR expansions that exceeded normal repeat lengths in autism cases compared to unaffected controls. By integrating single cell RNA and ATAC sequencing datasets of human postmortem brains, we prioritized STR loci in genes specifically expressed in cortical development stages. A deep learning method was used to predict functionality of ASD-associated STR loci. In ASD cases, rare STR expansions predominantly occurred in early cortical layer-specific genes involved in neurodevelopment, highlighting the cellular specificity of STR-associated genes in ASD risk. Leveraging deep learning prediction models, we demonstrated that these STR expansions disrupted the regulatory activity of enhancers and promoters, suggesting a potential mechanism through which they contribute to ASD pathogenesis. We found that individuals with ASD-associated STR expansions exhibited more severe ASD phenotypes and diminished adaptability compared to non-carriers. Short tandem repeat expansions in cortical layer-specific genes are associated with ASD and could potentially be a risk genetic factor for ASD. Our study is the first to show evidence of STR expansion associated with ASD in an under-investigated population.

UBQLN1 links proteostasis and mitochondria function to telomere maintenance in human embryonic stem cells

BackgroundTelomeres consist of repetitive DNA sequences at the chromosome ends to protect chromosomal stability, and primarily maintained by telomerase or occasionally by alternative telomere lengthening of telomeres (ALT) through recombination-based mechanisms. Additional mechanisms that may regulate telomere maintenance remain to be explored. Simultaneous measurement of telomere length and transcriptome in the same human embryonic stem cell (hESC) revealed that mRNA expression levels of UBQLN1 exhibit linear relationship with telomere length.MethodsIn this study, we first generated UBQLN1-deficient hESCs and compared with the wild-type (WT) hESCs the telomere length and molecular change at RNA and protein level by RNA-seq and proteomics. Then we identified the potential interacting proteins with UBQLN1 using immunoprecipitation-mass spectrometry (IP-MS). Furthermore, the potential mechanisms underlying the shortened telomeres in UBQLN1-deficient hESCs were analyzed.ResultsWe show that Ubiquilin1 (UBQLN1) is critical for telomere maintenance in human embryonic stem cells (hESCs) via promoting mitochondrial function. UBQLN1 deficiency leads to oxidative stress, loss of proteostasis, mitochondria dysfunction, DNA damage, and telomere attrition. Reducing oxidative damage and promoting mitochondria function by culture under hypoxia condition or supplementation with N-acetylcysteine partly attenuate the telomere attrition induced by UBQLN1 deficiency. Moreover, UBQLN1 deficiency/telomere shortening downregulates genes for neuro-ectoderm lineage differentiation.ConclusionsAltogether, UBQLN1 functions to scavenge ubiquitinated proteins, preventing their overloading mitochondria and elevated mitophagy. UBQLN1 maintains mitochondria and telomeres by regulating proteostasis and plays critical role in neuro-ectoderm differentiation.

Y chromosome shredding in Anopheles gambiae: Insight into the cellular dynamics of a novel synthetic sex ratio distorter.

Despite efforts to explore the genome of the malaria vector Anopheles gambiae, the Y chromosome of this species remains enigmatic. The large number of repetitive and heterochromatic DNA sequences makes the Y chromosome exceptionally difficult to fully assemble, hampering the progress of gene editing techniques and functional studies for this chromosome. In this study, we made use of a bioinformatic platform to identify Y-specific repetitive DNA sequences that served as a target site for a CRISPR/Cas9 system. The activity of Cas9 in the reproductive organs of males caused damage to Y-bearing sperm without affecting their fertility, leading to a strong female bias in the progeny. Cytological investigation allowed us to identify meiotic defects and investigate sperm selection in this new synthetic sex ratio distorter system. In addition, alternative promoters enable us to target the Y chromosome in specific tissues and developmental stages of male mosquitoes, enabling studies that shed light on the role of this chromosome in male gametogenesis. This work paves the way for further insight into the poorly characterised Y chromosome of Anopheles gambiae. Moreover, the sex distorter strain we have generated promises to be a valuable tool for the advancement of studies in the field of developmental biology, with the potential to support the progress of genetic strategies aimed at controlling malaria mosquitoes and other pest species.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements.

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

ULTRA-Effective Labeling of Repetitive Genomic Sequence.

In the age of long read sequencing, genomics researchers now have access to accurate repetitive DNA sequence (including satellites) that, due to the limitations of short read sequencing, could previously be observed only as unmappable fragments. Tools that annotate repetitive sequence are now more important than ever, so that we can better understand newly uncovered repetitive sequences, and also so that we can mitigate errors in bioinformatic software caused by those repetitive sequences. To that end, we introduce the 1.0 release of our tool for identifying and annotating locally-repetitive sequence, ULTRA (ULTRA Locates Tandemly Repetitive Areas). ULTRA is fast enough to use as part of an efficient annotation pipeline, produces state-of-the-art reliable coverage of repetitive regions containing many mutations, and provides interpretable statistics and labels for repetitive regions. It released under an open license, and available for download at https://github.com/TravisWheelerLab/ULTRA.

Molecular motor in a box: Attractors, repellers, and ratchets in chromatin

The spatiotemporal organization of eukaryotic DNA is controlled by an intricate combination of passive (thermally activated) and active (ATP-consuming) processes. Based on recent experimental insights into the mechanochemical cycle of chromatin remodelers, molecular machines that actively control nucleosome positions, we introduce a model to study the competition between active and passive mechanisms at the most basic layer of DNA packaging, the wrapping of DNA into nucleosomes. Depending on the level of remodeler activity, the positions of nucleosomes are controlled by either the bending or the stretching energy of the wrapped DNA involved. Since these energies are highly sequence dependent, DNA guides its own packaging. However, since this dependence differs for the two deformation modes, active processes can drive nucleosomes from their equilibrium positions. Furthermore, for repetitive DNA sequences, such as telomeres, we find thermal ratchets that propel nucleosomes in a preferred direction. Published by the American Physical Society 2024

Immune function of a C-type lectin with long tandem repeats and abundant threonine in the ridgetail white prawn Exopalaemon carinicauda

C-type lectins (CTLs) are an important class of pattern recognition receptors (PRRs) that exhibit structural and functional diversity in invertebrates. Repetitive DNA sequences are ubiquitous in eukaryotic genomes, representing distinct modes of genome evolution and promoting new gene generation. Our study revealed a new CTL that is composed of two long tandem repeats, abundant threonine, and one carbohydrate recognition domain (CRD) in Exopalaemon carinicauda and has been designated EcTR-CTL. The full-length cDNA of EcTR-CTL was 1242 bp long and had an open reading frame (ORF) of 999 bp that encoded a protein of 332 amino acids. The genome structure of EcTR-CTL contains 4 exons and 3 introns. The length of each repeat unit in EcTR-CTL was 198 bp, which is different from the short tandem repeats reported previously in prawns and crayfish. EcTR-CTL was abundantly expressed in the intestine and hemocytes. After Vibrio parahaemolyticus and white spot syndrome virus (WSSV) challenge, the expression level of EcTR-CTL in the intestine was upregulated. Knockdown of EcTR-CTL downregulated the expression of anti-lipopolysaccharide factor, crustin, and lysozyme during Vibrio infection. The recombinant CRD of EcTR-CTL (rCRD) could bind to bacteria, lipopolysaccharides, and peptidoglycans. Additionally, rCRD can directly bind to WSSV. These findings indicate that 1) CTLs with tandem repeats may be ubiquitous in crustaceans, 2) EcTR-CTL may act as a PRR to participate in the innate immune defense against bacteria via nonself-recognition and antimicrobial peptide regulation, and 3) EcTR-CTL may play a positive or negative role in the process of WSSV infection by capturing virions.

Chromosomal mapping of repetitive DNA and retroelement sequences and its implications for the chromosomal evolution process in Ctenoluciidae (Characiformes)

Ctenoluciidae is a Neotropical freshwater fish family composed of two genera, Ctenolucius (C. beani and C. hujeta) and Boulengerella (B. cuvieri, B. lateristriga, B. lucius, B. maculata, and B. xyrekes), which present diploid number conservation of 36 chromosomes and a strong association of telomeric sequences with ribosomal DNAs. In the present study, we performed chromosomal mapping of microsatellites and transposable elements (TEs) in Boulengerella species and Ctenolucius hujeta. We aim to understand how those sequences are distributed in these organisms’ genomes and their influence on the chromosomal evolution of the group. Our results indicate that repetitive sequences may had an active role in the karyotypic diversification of this family, especially in the formation of chromosomal hotspots that are traceable in the diversification processes of Ctenoluciidae karyotypes. We demonstrate that (GATA)n sequences also accumulate in the secondary constriction formed by the 18 S rDNA site, which shows consistent size heteromorphism between males and females in all Boulengerella species, suggesting an initial process of sex chromosome differentiation.

Celine, a long interspersed nuclear element retrotransposon, colonizes in the centromeres of poplar chromosomes.

Centromeres in most multicellular eukaryotes are composed of long arrays of repetitive DNA sequences. Interestingly, several transposable elements, including the well-known long terminal repeat centromeric retrotransposon of maize (CRM), were found to be enriched in functional centromeres marked by the centromeric histone H3 (CENH3). Here, we report a centromeric long interspersed nuclear element (LINE), Celine, in Populus species. Celine has colonized preferentially in the CENH3-associated chromatin of every poplar chromosome, with 84% of the Celine elements localized in the CENH3-binding domains. In contrast, only 51% of the CRM elements were bound to CENH3 domains in Populus trichocarpa. These results suggest different centromere targeting mechanisms employed by Celine and CRM elements. Nevertheless, the high target specificity seems to be detrimental to further amplification of the Celine elements, leading to a shorter life span and patchy distribution among plant species compared with the CRM elements. Using a phylogenetically guided approach, we were able to identify Celine-like LINE elements in tea plant (Camellia sinensis) and green ash tree (Fraxinus pennsylvanica). The centromeric localization of these Celine-like LINEs was confirmed in both species. We demonstrate that the centromere targeting property of Celine-like LINEs is of primitive origin and has been conserved among distantly related plant species.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Influence of breed and environment on leukocyte telomere length in cattle.

High milk yield is associated with reduced longevity in high-producing dairy cattle breeds. Pre-term culling leads to high replacement heifer demand and economic losses for the dairy industry. Selection for this trait is limited because of low heritability and difficulties in phenotype measurement. Telomeres are elements found at the ends of chromosomes, consisting of repetitive DNA sequences, several thousand base pairs in length, coupled with nucleoprotein complexes. Eventually, in humans and most other animals, telomere length reduces with age. When telomeric DNA is truncated to a critical length, cell ageing, cell cycle arrest, and apoptosis are induced. As a result, telomere length can be considered as a predictor of health risks and an individual's lifespan. The leukocyte telomere length may be used as a proxy phenotype of productive lifespan to improve cattle selection. Our objectives were to assess the effects of breed and breed group (dairy vs. beef) on the leukocyte telomere length and to estimate the effect of cold climate on this trait in Kalmyk cattle populations from the South (Rostov Oblast) and Far North (Republic of Sakha) regions of Russia. The leukocyte telomere lengths were estimated computationally from whole-genome resequencing data. We leveraged data on leukocyte telomere length, sex, and age of 239 animals from 17 cattle breeds. The breed factor had a significant effect on leukocyte telomere length across our sample. There was no difference in leukocyte telomere length between dairy and beef groups. The population factor had a significant effect on leukocyte telomere length in Kalmyk animals. In conclusion, we found that breed, but not breed group (dairy vs. beef), was significantly associated with leukocyte telomere length in cattle. Residence in colder climates was associated with longer leukocyte telomere length in Kalmyk breed cattle.

Dominance of transposable element-related satDNAs results in great complexity of “satDNA library” and invokes the extension towards “repetitive DNA library”

Research on bivalves is fast-growing, including genome-wide analyses and genome sequencing. Several characteristics qualify oysters as a valuable model to explore repetitive DNA sequences and their genome organization. Here we characterize the satellitomes of five species in the family Ostreidae (Crassostrea angulata, C. virginica, C. hongkongensis, C. ariakensis, Ostrea edulis), revealing a substantial number of satellite DNAs (satDNAs) per genome (ranging between 33 and 61) and peculiarities in the composition of their satellitomes. Numerous satDNAs were either associated to or derived from transposable elements, displaying a scarcity of transposable element-unrelated satDNAs in these genomes. Due to the non-conventional satellitome constitution and dominance of Helitron-associated satDNAs, comparative satellitomics demanded more in-depth analyses than standardly employed. Comparative analyses (including C. gigas, the first bivalve species with a defined satellitome) revealed that 13 satDNAs occur in all six oyster genomes, with Cg170/HindIII satDNA being the most abundant in all of them. Evaluating the “satDNA library model” highlighted the necessity to adjust this term when studying tandem repeat evolution in organisms with such satellitomes. When repetitive sequences with potential variation in the organizational form and repeat-type affiliation are examined across related species, the introduction of the terms “TE library” and “repetitive DNA library” becomes essential.