The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani.

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum.

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

Satellitome analyses in nematodes illuminate complex species history and show conserved features in satellite DNAs.

BACKGROUND: Satellite DNAs (satDNAs) are tandemly repeated non-coding DNA sequences that belong to the most abundant and the fastest evolving parts of the eukaryotic genome. A satellitome represents the collection of different satDNAs in a genome. Due to extreme diversity and methodological difficulties to characterize and compare satDNA collection in complex genomes, knowledge on their putative functional constraints and capacity to participate in genome evolution remains rather elusive. SatDNA transcripts have been detected in many species, however comparative studies of satDNA transcriptome between species are extremely rare. RESULTS: We conducted a genome-wide survey and comparative analyses of satellitomes among different closely related Meloidogyne spp. nematodes. The evolutionary trends of satDNAs suggest that each round of proposed polyploidization in the evolutionary history is concomitant with the addition of a new set of satDNAs in the satellitome of any particular Meloidogyne species. Successive incorporation of new sets of satDNAs in the genome along the process of polyploidization supports multiple hybridization events as the main factor responsible for the formation of these species. Through comparative analyses of 83 distinct satDNAs, we found a CENP-B box-like sequence motif conserved among 11 divergent satDNAs (similarity ranges from 36 to 74%). We also found satDNAs that harbor a splice leader (SL) sequence which, in spite of overall divergence, shows conservation across species in two putative functional regions, the 25-nt SL exon and the Sm binding site. Intra- and interspecific comparative expression analyses of the complete satDNA set in the analyzed Meloidogyne species revealed transcription profiles including a subset of 14 actively transcribed satDNAs. Among those, 9 show active transcription in every species where they are found in the genome and throughout developmental stages. CONCLUSIONS: Our results demonstrate the feasibility and power of comparative analysis of the non-coding repetitive genome for elucidation of the origin of species with a complex history. Although satDNAs generally evolve extremely quickly, the comparative analyses of 83 satDNAs detected in the analyzed Meloidogyne species revealed conserved sequence features in some satDNAs suggesting sequence evolution under selective pressure. SatDNAs that are actively transcribed in related genomes and throughout nematode development support the view that their expression is not stochastic.

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani.

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.

Isolation of High Molecular Weight DNA from the Model Beetle Tribolium for Nanopore Sequencing.

The long-read Nanopore sequencing has been recently applied for assembly of complex genomes and analysis of linear genome organization. The most critical factor for successful long-read sequencing is extraction of high molecular weight (HMW) DNA of sufficient purity and quantity. The challenges associated with input DNA quality are further amplified when working with extremely small insects with hard exoskeletons. Here, we optimized the isolation of HMW DNA from the model beetle Tribolium and tested for use in Nanopore sequencing. We succeeded in overcoming all the difficulties in HMW handling and library preparation that were encountered when using published protocols and commercial kits. Isolation of nuclei and subsequent purification of DNA on an anion-exchange chromatography column resulted in genomic HMW DNA that was efficiently relaxed, of optimal quality and in sufficient quantity for Nanopore MinION sequencing. DNA shearing increased average N50 read values up to 26 kb and allowed us to use a single flow cell in multiple library loads for a total output of more than 13 Gb. Although our focus was on T. castaneum and closely related species, we expect that this protocol, with appropriate modifications, could be extended to other insects, particularly beetles.

The Centromere Histone Is Conserved and Associated with Tandem Repeats Sharing a Conserved 19-bp Box in the Holocentromere of Meloidogyne Nematodes.

Although centromeres have conserved function, centromere-specific histone H3 (CenH3) and centromeric DNA evolve rapidly. The centromere drive model explains this phenomenon as a consequence of the conflict between fast-evolving DNA and CenH3, suggesting asymmetry in female meiosis as a crucial factor. We characterized evolution of the CenH3 protein in three closely related, polyploid mitotic parthenogenetic species of the Meloidogyne incognita group, and in the distantly related meiotic parthenogen Meloidogyne hapla. We identified duplication of the CenH3 gene in a putative sexual ancestral Meloidogyne. We found that one CenH3 (αCenH3) remained conserved in all extant species, including in distant Meloidogyne hapla, whereas the other evolved rapidly and under positive selection into four different CenH3 variants. This pattern of CenH3 evolution in Meloidogyne species suggests the subspecialization of CenH3s in ancestral sexual species. Immunofluorescence performed on mitotic Meloidogyne incognita revealed a dominant role of αCenH3 on its centromere, whereas the other CenH3s have lost their function in mitosis. The observed αCenH3 chromosome distribution disclosed cluster-like centromeric organization. The ChIP-Seq analysis revealed that in M. incognita αCenH3-associated DNA dominantly comprises tandem repeats, composed of divergent monomers which share a completely conserved 19-bp long box. Conserved αCenH3-associated DNA is also confirmed in the related mitotic Meloidogyne incognita group species suggesting preservation of both centromere protein and DNA constituents. We hypothesize that the absence of centromere drive in mitosis might allow for CenH3 and its associated DNA to achieve an equilibrium in which they can persist for long periods of time.

CenH3 distribution reveals extended centromeres in the model beetle Tribolium castaneum.

Centromeres are chromosomal domains essential for kinetochore assembly and correct chromosome segregation. Inconsistent in their underlying DNA sequences, centromeres are defined epigenetically by the presence of the centromere-specific histone H3 variant CenH3. Most of the analyzed eukaryotes have monocentric chromosomes in which CenH3 proteins deposit into a single, primary constriction visible at metaphase chromosomes. Contrary to monocentrics, evolutionary sporadic holocentric chromosomes lack a primary constriction and have kinetochore activity distributed along the entire chromosome length. In this work, we identified cCENH3 protein, the centromeric H3 histone of the coleopteran model beetle Tribolium castaneum. By ChIP-seq analysis we disclosed that cCENH3 chromatin assembles upon a repertoire of repetitive DNAs. cCENH3 in situ mapping revealed unusually elongated T. castaneum centromeres that comprise approximately 40% of the chromosome length. Being the longest insect regional centromeres evidenced so far, T. castaneum centromeres are characterized by metapolycentric structure composed of several individual cCENH3-containing domains. We suggest that the model beetle T. castaneum with its metapolycentromeres could represent an excellent model for further studies of non-canonical centromeres in insects.

Methylation profile of a satellite DNA constituting the intercalary G+C-rich heterochromatin of the cut trough shell Spisula subtruncata (Bivalvia, Mactridae).

Tandemly repeated DNAs usually constitute significant portions of eukaryotic genomes. In bivalves, however, repetitive DNAs are habitually not widespread. In our search for abundant repetitive DNAs in trough shells, we discovered a novel satellite DNA, SSUsat, which constitutes at least 1.3% of the genome of Spisula subtruncata. As foreseen by the satellite DNA library hypothesis, we confirmed that this satellite DNA is also present in two other Mactridae species, showing a highly conserved nucleotide sequence together with a dramatic diminution in the number of repeats. Predominantly located at the G + C-rich intercalary heterochromatin of S. subtruncata, SSUsat displays several DNA methylation peculiarities. The level of methylation of SSUsat is high (3.38%) in comparison with bivalve standards and triplicates the mean of the S. subtruncata genome (1.13%). Methylation affects not only the cytosines in CpG dinucleotides but also those in CHH and CHG trinucleotides, a feature common in plants but scarce and without any clear known relevance in animals. SSUsat segments enriched in methylated cytosines partly overlap those showing higher sequence conservation. The presence of a chromosome pair showing an accumulation of markedly under-methylated SSUsat monomers additionally indicates that the methylation processes that shape repetitive genome compartments are quite complex.

Structural and functional liaisons between transposable elements and satellite DNAs.

Transposable elements (TEs) and satellite DNAs (satDNAs) are typically identified as major repetitive DNA components in eukaryotic genomes. TEs are DNA segments able to move throughout a genome while satDNAs are tandemly repeated sequences organized in long arrays. Both classes of repetitive sequences are extremely diverse, and many TEs and satDNAs exist within a genome. Although they differ in structure, genomic organization, mechanisms of spread, and evolutionary dynamics, TEs and satDNAs can share sequence similarity and organizational patterns, thus indicating that complex mutual relationships can determine their evolution, and ultimately define roles they might have on genome architecture and function. Motivated by accumulating data about sequence elements that incorporate features of both TEs and satDNAs, here we present an overview of their structural and functional liaisons.

Centromere identity from the DNA point of view.

The centromere is a chromosomal locus responsible for the faithful segregation of genetic material during cell division. It has become evident that centromeres can be established literally on any DNA sequence, and the possible synergy between DNA sequences and the most prominent centromere identifiers, protein components, and epigenetic marks remains uncertain. However, some evolutionary preferences seem to exist, and long-term established centromeres are frequently formed on long arrays of satellite DNAs and/or transposable elements. Recent progress in understanding functional centromere sequences is based largely on the high-resolution DNA mapping of sequences that interact with the centromere-specific histone H3 variant, the most reliable marker of active centromeres. In addition, sequence assembly and mapping of large repetitive centromeric regions, as well as comparative genome analyses offer insight into their complex organization and evolution. The rapidly advancing field of transcription in centromere regions highlights the functional importance of centromeric transcripts. Here, we comprehensively review the current state of knowledge on the composition and functionality of DNA sequences underlying active centromeres and discuss their contribution to the functioning of different centromere types in higher eukaryotes.

Genomic size of CENP-A domain is proportional to total alpha satellite array size at human centromeres and expands in cancer cells.

Human centromeres contain multi-megabase-sized arrays of alpha satellite DNA, a family of satellite DNA repeats based on a tandemly arranged 171 bp monomer. The centromere-specific histone protein CENP-A is assembled on alpha satellite DNA within the primary constriction, but does not extend along its entire length. CENP-A domains have been estimated to extend over 2,500 kb of alpha satellite DNA. However, these estimates do not take into account inter-individual variation in alpha satellite array sizes on homologous chromosomes and among different chromosomes. We defined the genomic distance of CENP-A chromatin on human chromosomes X and Y from different individuals. CENP-A chromatin occupied different genomic intervals on different chromosomes, but despite inter-chromosomal and inter-individual array size variation, the ratio of CENP-A to total alpha satellite DNA size remained consistent. Changes in the ratio of alpha satellite array size to CENP-A domain size were observed when CENP-A was overexpressed and when primary cells were transformed by disrupting interactions between the tumor suppressor protein Rb and chromatin. Our data support a model for centromeric domain organization in which the genomic limits of CENP-A chromatin varies on different human chromosomes, and imply that alpha satellite array size may be a more prominent predictor of CENP-A incorporation than chromosome size. In addition, our results also suggest that cancer transformation and amounts of centromeric heterochromatin have notable effects on the amount of alpha satellite that is associated with CENP-A chromatin.

TCAGG, an alternative telomeric sequence in insects.

The TTAGG repeat, the only determined telomerase-dependent sequence in the Insecta, is generally reputed to be the canonical telomeric motif within the class. By studying the distribution of telomeric DNAs in 30 coleopteran beetles using Southern hybridization, BAL 31 DNA end-degradation assay and fluorescence in situ hybridization, we showed that arrays built of a TCAGG repeat substitute for (TTAGG)n sequences in all tested species within the superfamily Tenebrionoidea. We also provided the experimental evidence that (TCAGG)n repeats represent the terminal sequences on all chromosomes of the model species Tribolium castaneum. (TCAGG)n repeats are therefore promoted as the first sequence-motif alternative to TTAGG-type chromosome ends in insects. Detection of species negative for both TTAGG and TCAGG reveals that, although widespread, these motifs are not ubiquitous telomeric sequences within the order Coleoptera. In addition, Timarcha balearica proved to be a species that harbors (TTAGG)n repeats, but not at telomeric positions, thus further increasing the complexity of telomeric DNAs. Our experiments discarded CTAGG, CTGGG, TTGGG, and TTAGGG variants as potential replacements in TTAGG/TCAGG-negative species, indicating that chromosome termini of these beetles comprise other form(s) of telomeric sequences and telomere maintenance mechanisms.

Parallelism in evolution of highly repetitive DNAs in sibling species.

Characterization of heterochromatin in the flour beetle Tribolium audax revealed two highly repetitive DNA families, named TAUD1 and TAUD2, which together constitute almost 60% of the whole genome. Both families originated from a common ancestral approximately 110-bp repeating unit. Tandem arrangement of these elements in TAUD1 is typical for satellite DNAs, whereas TAUD2 represents a dispersed family based on 1412-bp complex higher-order repeats composed of inversely oriented approximately 110 bp units. Comparison with repetitive DNAs in the sibling species Tribolium madens showed similarities in nucleotide sequence and length of basic repeating units and also revealed structural and organizational parallelism in tandem and dispersed families assembled from these elements. In both Tribolium species, one tandem and one dispersed family build equivalent distribution patterns in the pericentromeric heterochromatin of all chromosomes including supernumeraries. Differences in the nucleotide sequence and in the complexity of higher-order structures between families of the same type suggest a scenario according to which rearranged variants of the corresponding ancestral families were formed and distributed in genomes during or after the speciation event, following the same principles independently in each descendant species. We assume that random effects of sequence dynamics should be constrained by organizational and structural features of repeating units and possible requirements for spatial distribution of particular sequence elements. An interspersed pattern of repetitive families also points to the intensive recombination events in heterochromatin. Synergy between the meiotic bouquet stage and satellite DNA sequence dynamics could make a positive feedback loop that promotes the observed genome-wide distribution. At the same time, considering the abundance of these DNAs in heterochromatin spanning the (peri)centromeric chromosomal segments, we speculate that diverged repetitive sequences might represent the DNA basis of reproductive barrier between the two sibling species.

Histone modifications within the human X centromere region.

Human centromeres are multi-megabase regions of highly ordered arrays of alpha satellite DNA that are separated from chromosome arms by unordered alpha satellite monomers and other repetitive elements. Complexities in assembling such large repetitive regions have limited detailed studies of centromeric chromatin organization. However, a genomic map of the human X centromere has provided new opportunities to explore genomic architecture of a complex locus. We used ChIP to examine the distribution of modified histones within centromere regions of multiple X chromosomes. Methylation of H3 at lysine 4 coincided with DXZ1 higher order alpha satellite, the site of CENP-A localization. Heterochromatic histone modifications were distributed across the 400-500 kb pericentromeric regions. The large arrays of alpha satellite and gamma satellite DNA were enriched for both euchromatic and heterochromatic modifications, implying that some pericentromeric repeats have multiple chromatin characteristics. Partial truncation of the X centromere resulted in reduction in the size of the CENP-A/Cenp-A domain and increased heterochromatic modifications in the flanking pericentromere. Although the deletion removed approximately 1/3 of centromeric DNA, the ratio of CENP-A to alpha satellite array size was maintained in the same proportion, suggesting that a limited, but defined linear region of the centromeric DNA is necessary for kinetochore assembly. Our results indicate that the human X centromere contains multiple types of chromatin, is organized similarly to smaller eukaryotic centromeres, and responds to structural changes by expanding or contracting domains.

Human gamma-satellite DNA maintains open chromatin structure and protects a transgene from epigenetic silencing.

The role of repetitive DNA sequences in pericentromeric regions with respect to kinetochore/heterochromatin structure and function is poorly understood. Here, we use a mouse erythroleukemia cell (MEL) system for studying how repetitive DNA assumes or is assembled into different chromatin structures. We show that human gamma-satellite DNA arrays allow a transcriptionally permissive chromatin conformation in an adjacent transgene and efficiently protect it from epigenetic silencing. These arrays contain CTCF and Ikaros binding sites. In MEL cells, this gamma-satellite DNA activity depends on binding of Ikaros proteins involved in differentiation along the hematopoietic pathway. Given our discovery of gamma-satellite DNA in pericentromeric regions of most human chromosomes and a dynamic chromatin state of gamma-satellite arrays in their natural location, we suggest that gamma-satellite DNA represents a unique region of the functional centromere with a possible role in preventing heterochromatin spreading beyond the pericentromeric region.

Satellite DNA junctions identify the potential origin of new repetitive elements in the beetle Tribolium madens.

Two related satellite DNA families (satellite I and satellite II) with complex higher-order repeat (HOR) monomers represent major DNA components equilocated in the pericentromeric heterochromatin of all Tribolium madens chromosomes. Fragments obtained upon genomic DNA restriction revealed two subfamilies of satellite II monomers, and also identified regions of transition between satellite I and satellite II sequences. The two subfamilies differ not only in diagnostic nucleotides, but also in flipped orientation of constituent subunits. Hybrid genomic fragments comprise directly linked satellite I and satellite II monomers that cannot be distinguished from randomly cloned monomers of corresponding families. An exception is the most proximal satellite I monomer in the hybrid fragment named TMADhinf, which shows sequence divergence typical for repeats evolving at array ends, in zones of low homogenization efficiency. This pattern points to the extensive rearrangement processes generating abrupt transitions between satellite arrays combined with array maintenance by unequal crossover. Switching points between adjacent satellites as well as the edges of flipped subunits are localized within a short sequence segment, indicating a preferential site of recombination within satellite subunits. Multiple copies of TMADhinf junction fragment support the hypothesis that sites of evolutionary origin of novel satellite repeat (sub)families can be localized at array ends, in regions of enhanced sequence divergence.

Preservation and high sequence conservation of satellite DNAs suggest functional constraints.

Due to a high evolutionary turnover many satellite DNAs are restricted to a group of closely related species. Here we demonstrate that the satellite DNA family PSUB, abundant in the beetle Palorus subdepressus, is distributed in a low number of copies among diverse taxa of Coleoptera (Insecta), some of them separated for an evolutionary period of up to 60 Myr. Comparison of PSUB cloned from the species Tribolium brevicornis with the PSUB family previously characterized in Palorus subdepressus revealed high sequence conservation and absence of fixed species-specific mutations. The most polymorphic sites are those with ancestral mutations shared among clones of both species. Since the ancestral mutations contribute significantly to overall diversity, it could be proposed that a similar mutational profile already existed in an ancestral species. The pattern of variability along the satellite monomer is characterized by the presence of conserved and variable regions. The nonrandom pattern of variability as well as the absence of sequence divergence is also discerned for PRAT satellite DNA, cloned previously from two Palorus species and a distantly related Pimelia elevata. Since PRAT and PSUB are present in parallel in diverse taxa of Coleoptera, we propose that their long evolutionary preservation suggests a possible functional significance. This indication is additionally supported not only by the high evolutionary conservation of the sequences, but also by the presence of significantly conserved and variable regions along the monomers.

Long inversely oriented subunits form a complex monomer of Tribolium brevicornis satellite DNA.

Highly abundant satellite DNA named TBREV is detected and characterized in the beetle Tribolium brevicornis (Insecta: Coleoptera). An outstanding peculiarity of the TBREV satellite monomer is its complex structure based on the two approximately 470-bp-long subunits, inversely oriented within a 1061-bp-long monomer sequence. The proposed evolutionary history demonstrates a clear trend toward increased complexity and length of the TBREV satellite monomer. This tendency has been observed on three levels: first as direct and inverted duplications of short sequence motifs, then by inverse duplication of the approximately 470-bp sequence segment, and, finally, by spread of inversely duplicated elements in a higher-order register and formation of extant monomers. Inversely oriented subunits share a similarity of 82% and have a high capacity to form a thermodynamically stable dyad structure that is, to our knowledge, the longest ever described in any satellite monomer. Analysis of divergences between inversely oriented subunits shows a tendency to a further reduction in similarity between them. Except in its centromeric localization, the TBREV satellite does not show similarity to other known Tribolium satellites, either in nucleotide sequence or in monomer length and complexity. However, TBREV shares common features of other Tribolium satellites that might be under functional constraints: nonconstant rate of evolution along the monomer sequence, short inverted repeats in the vicinity of an A+T tract, nonrandom distribution of A or T >/=3 tracts, and CENP-B box-like motifs. Although long inverted subunits might reinforce structural characteristics of the satellite monomer, their nucleotide sequence does not seem to be under constraints in order to preserve the dyad structure.

Conserved patterns in the evolution of Tribolium satellite DNAs.

Two satellite DNAs, TANAPH and TDEST, isolated from the beetle species Tribolium anaphe and Tribolium destructor, respectively, are characterized and compared with previously described Tribolium satellites, in order to deduce possible constraints on satellite sequence evolution between closely related species. Sequence diversity analysis of cloned monomers reveals the presence of variable and conserved segments in both satellites. In addition, non-random organization of As or Ts and their periodical distribution in the form of A or T >/=3 tracts, as well as CENP-B box-like motifs and dyad structures have been found in both satellites. Similar structural features are also present in satellites from other Tribolium species. We therefore propose that they, together with the observed non-constant rate of evolution along the satellite sequence, could be related to putative protein binding sites and suggest a possible selective pressure affecting these sequences. Tribolium satellites, including TANAPH and TDEST, are located in the pericentromeric heterochromatin of all chromosomes of the corresponding species. Since satellites from different species exhibit no significant sequence homology, we propose that they did not originate from a common ancestral sequence. More probably, they derive from simple sequence modules some of which could represent protein binding sites. Shuffling of simple sequence modules could generate different satellites, able to perform a similar role in different species.

Sequence of PRAT satellite DNA "frozen" in some Coleopteran species.

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera ( Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50-60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.