Detection of satellite DNA in Palorus ratzeburgii: analysis of curvature profiles and comparison with Tenebrio molitor satellite DNA

Cytogenetic Analysis, Heterochromatin Characterization and Location of the rDNA Genes of Hycleus scutellatus (Coleoptera, Meloidae); A Species with an Unexpected High Number of rDNA Clusters

Simple Summary

The family Meloidae contains approximately 3000 species, commonly known as blister beetles for their ability to secrete a substance called cantharidin, which causes irritation and blistering in contact with animal or human skin. In recent years there have been numerous studies focused on the anticancer action of cantharidin and its derivatives. Despite the recent interest in blister beetles, cytogenetic and molecular studies in this group are scarce and most of them use only classical chromosome staining techniques. The main aim of our study was to provide new information in Meloidae. In this study, cytogenetic and molecular analyses were applied for the first time in the family Meloidae. We applied fluorescence staining with DAPI and the position of ribosomal DNA in Hycleus scutellatus was mapped by FISH. Hycleus is one of the most species-rich genera of Meloidae but no cytogenetic data have yet been published for this particular genus. Additionally, we isolated a satellite DNA family located within the pericentromeric regions of all chromosomes. The results obtained in this study may be a suitable starting point to initiate more extensive cytogenetic analyses in this important species-rich genus, and in the family Meloidae in general.

Abstract

Meloidae are commonly known as blister beetles, so called for the secretion of cantharidin, a toxic substance that causes irritation and blistering. There has been a recent increase in the interest of the cantharidin anticancer potential of this insect group. Cytogenetic and molecular data in this group are scarce. In this study, we performed a karyotype analysis of Hycleus scutellatus, an endemic species of the Iberian Peninsula. We determined its chromosome number, 2n = 20, as well as the presence of the X and Y sex chromosomes. In addition to a karyotype analysis, we carried out DAPI staining. By fluorescence in situ hybridization we mapped the rDNA clusters on 12 different chromosomes. Compared to others, this species shows an unusually high number of chromosomes carrying rDNA. This is one of the highest numbers of rDNA sites found in the Polyphaga suborder (Coleoptera). Additionally, we isolated a satellite DNA family (Hyscu-H), which was located within the pericentromeric regions of all chromosomes, including the sex chromosomes. The results suggest that Hyscu-H is likely to be one of the most abundant satellite DNA repeats in H. scutellatus.

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

Divergence in centromere structure distinguishes related genomes in Coix lacryma-jobi and its wild relative

Knowledge about the composition and structure of centromeres is critical for understanding how centromeres perform their functional roles. Here, we report the sequences of one centromere-associated bacterial artificial chromosome clone from a Coix lacryma-jobi library. Two Ty3/gypsy-class retrotransposons, centromeric retrotransposon of C. lacryma-jobi (CRC) and peri-centromeric retrotransposon of C. lacryma-jobi, and a (peri)centromere-specific tandem repeat with a unit length of 153 bp were identified. The CRC is highly homologous to centromere-specific retrotransposons reported in grass species. An 80-bp DNA region in the 153-bp satellite repeat was found to be conserved to centromeric satellite repeats from maize, rice, and pearl millet. Fluorescence in situ hybridization showed that the three repetitive sequences were located in (peri-)centromeric regions of both C. lacryma-jobi and Coix aquatica. However, the 153-bp satellite repeat was only detected on 20 out of the 30 chromosomes in C. aquatica. Immunostaining with an antibody against rice CENH3 indicates that the 153-bp satellite repeat and CRC might be both the major components for functional centromeres, but not all the 153-bp satellite repeats or CRC sequences are associated with CENH3. The evolution of centromeric repeats of C. lacryma-jobi during the polyploidization was discussed.

Centromere-competent DNA: structure and evolution

Although extant data favour centromere being an epigenetic structure, it is also clear that centromere formation is based on DNA, in particular, tandemly repeated satellite DNA and its transcripts. Presence of conserved structural motifs within satellite DNAs such as periodically distributed AT tracts, protein binding sites, or promoter elements indicate that despite sequence flexibility, there are structural determinants that are prerequisite for centromere function. In addition, existence of functional centromeric DNA transcripts indicates possible importance of structural elements at the level of RNA secondary or tertiary structure. Rapid centromere evolution is explained by homologous recombination followed by extrachromosomal rolling circle replication. This could lead to amplification of different satellite sequences within a genome. However, only those satellites that have inherent centromere-competence in the form of structural requirements necessary for centromere function are after amplification fixed in a population as a new centromere.

RNA Pol II promotes transcription of centromeric satellite DNA in beetles

Transcripts of centromeric satellite DNAs are known to play a role in heterochromatin formation as well as in establishment of the kinetochore. However, little is known about basic mechanisms of satellite DNA expression within constitutive heterochromatin and its regulation. Here we present comprehensive analysis of transcription of abundant centromeric satellite DNA, PRAT from beetle Palorus ratzeburgii (Coleoptera). This satellite is characterized by preservation and extreme sequence conservation among evolutionarily distant insect species. PRAT is expressed in all three developmental stages: larvae, pupae and adults at similar level. Transcripts are abundant comprising 0.033% of total RNA and are heterogeneous in size ranging from 0.5 kb up to more than 5 kb. Transcription proceeds from both strands but with 10 fold different expression intensity and transcripts are not processed into siRNAs. Most of the transcripts (80%) are not polyadenylated and remain in the nucleus while a small portion is exported to the cytoplasm. Multiple, irregularly distributed transcription initiation sites as well as termination sites have been mapped within the PRAT sequence using primer extension and RLM-RACE. The presence of cap structure as well as poly(A) tails in a portion of the transcripts indicate RNA polymerase II–dependent transcription and a putative polymerase II promoter site overlaps the most conserved part of the PRAT sequence. The treatment of larvae with alpha-amanitin decreases the level of PRAT transcripts at concentrations that selectively inhibit pol II activity. In conclusion, stable, RNA polymerase II dependant transcripts of abundant centromeric satellite DNA, not regulated by RNAi, have been identified and characterized. This study offers a basic understanding of expression of highly abundant heterochromatic DNA which in beetle species constitutes up to 50% of the genome.

Functional elements residing within satellite DNAs

Satellite DNAs represent a fast-evolving portion of the eukaryotic genome whose evolution is proposed to be driven by the stochastic process of molecular drive. Recent results indicate that satellite DNAs are subject to certain structural constraints, which are probably related to their interaction with proteins involved in the establishment of specific chromatin structures. The evolutionary persistence and high sequence conservation of some satellites, as well as the presence of stage- or tissue-specific, differentially expressed transcripts in several species, are consistent with the hypothesis that satellite DNA could have a regulatory role in eukaryotic organisms. Although the role of most transcripts is not known, some act as precursors of small interfering RNAs, which are now recognized as having an important role in chromatin modulation and the control of gene expression. Furthermore, some transcripts are involved in the cellular response to stress.

Molecular characterization of the first satellite DNA with CENP-B and CDEIII motifs in the bat Pipistrellus kuhli

The centromere is an essential structure in the chromosomes of all eukariotes and is central to the mechanism that ensures proper segregation during mitosis and meiosis. The comparison of DNA sequence motifs, organization and kinetocore components from yeast to man is beginning to indicate that, although centromeres are highly variable DNA elements, a conserved pattern of sequence arrangement and function is emerging. We have identified and characterized the first satellite DNA (P.k.SAT) from microbat species Pipistrellus kuhli. The presence of mammalian CENP-B box and yeast CDEIII box could indicate the participation of P.k.SAT in centromere organization.

Genomic dynamics of a low-copy-number satellite DNA family in Leptestheria dahalacensis (Crustacea, Branchiopoda, Conchostraca)

The LEP150 satellite DNA (satDNA) family found in Leptestheria dahalacensis (Ruppel, 1837) (Conchostraca) is a low-copy-number satellite with a canonical monomer of 150 bp. Nucleotide variation analyses suggest a 14-bp palindromic region as a possible protein binding site with constraints acting on the whole sequence but a 25-bp variable box. Besides the head-to-tail arrangement of 150 bp monomers, multimers analyses evidenced incomplete monomers, one duplication event, and three inversions. Both observed rearrangements and the higher values of sequence variability scored suggest that rearranged monomers reside in regions with a lower degree of homogenisation efficiency. Sixty-seven percent of the breakpoints occurs at kinkable dinucleotides, thus supporting their role in rearrangements as documented in alphoid satDNA recombination events. Monomers of different lengths may result from crossing over between repeats misaligned through the direct and inverted subrepeats of LEP150 monomers. ANOVA results indicate that the same range of sequence diversity is experienced at the individual and population ranks; therefore, the evolution of the L. dahalacensis satDNA is concerted.

Isolation and characterization of two families of satellite DNA with repetitive units of 135 bp and 2.5 kb in the ant Monomorium subopacum (Hymenoptera, Formicidae)

Analyzing the satellite DNA in the ant species Monomorium subopacum we found two unrelated families of satellite DNA. Because these satellite DNA families were isolated using the two enzymes HaeIII and EcoRI we called the two families HaeIII and EcoRI family, respectively. The HaeIII family proved to be organized in a 135-bp basic unit repeat, the EcoRI family in a 2.5-kb basic unit repeat. The latter represents perhaps the longest satellite DNA isolated up to now in insects. The HaeIII family apparently comprises about 10% of the total genomic DNA whereas the EcoRI family represents only about 1-2%. A comparative analysis of the two satellite DNA sequences showed no homology between the two families although both sequences possessed long A and T stretches. Eight of the 34 chromosomes showed hybridization with the HaeIII family and hybridization signals are visible in six chromosomes with the EcoRI family. Analysis of the electrophoretic mobility of satellite DNA on non-denaturing polyacrylamide showed that the HaeIII family is only slightly curved. However, the unit of the EcoRI satellite DNA family has curvature, especially the first 1000 bp of the monomeric repeat, in which this DNA is AT rich and has numerous A and T stretches. There are also internal inverted subrepeats in each family. The sequences of satellite DNA families found in Monomorium subopacum are different from the sequences of other satellite DNAs cloned in insects, including other species of ants.

Characterization of the heterochromatin of the darkling beetle Misolampus goudoti: cloning of two satellite DNA families and digestion of chromosomes with restriction enzymes

The darkling beetle Misolampus goudoti Er. has 58% of C-banded chromosome material. In this paper we deal with the study of the heterochromatin of this insect both by molecular and cytogenetical methods. Two different satellite DNA families have been characterized in Misolampus goudoti by agarose gel electrophoresis of EcoRI and PstI restriction fragments, respectively. The EcoRI family is composed of a monomeric unit of 196 bp (64.3% A-T rich) DNA sequence, representing about 120,000 copies per haploid genome. The presence of frequent intermediate-size satellite variants and an internal direct repetition of 61 bp in the EcoRI repetitive main monomer suggest that the evolution of this satellite proceeded by unequal crossing-over, occurring both within and between the 196 bp unit. Another highly repetitive sequence, defined by digestion of genomic DNA with PstI, has a more complex unit of 1.2 kb with about 70,000 copies per haploid genome. In situ digestion of M. goudoti chromosomes with restriction enzymes shows a non-specific chromosome DNA extraction from pericentromeric positions with EcoRI and chromosome specific extraction of DNA with PstI and HinfI. This is discussed in relation to the chromosomal location of both satellites.

Comparative analysis of different satellite DNAs in four Mytilus species

We report the characterization of three satellite DNAs in four species of mussel: Mytilus edulis, Mytilus galloprovincialis, Mytilus trossulus, and Mytilus californianus. The monomers of the Apa I satellite DNAs were 173, 161, and 166 bp long. These satellite monomers were used to construct phylogenetic trees to infer relationships among these species. The topologies obtained clearly indicate that M. californianus is the most divergent species with respect to the other three. Furthermore, localization of satellite DNAs on metaphase chromosomes was performed using fluorescent in situ hybridization (FISH). Fluorescent signals revealed a different organization and distribution of these three satellite DNAs.

Comparative study of satellite sequences and phylogeny of five species from the genusPalorus(Insecta, Coleoptera)

Major satellite sequences are analysed in the three tenebrionid beetles Palorus cerylonoides, P. genalis, and P. ficicola, and compared with the ones from P. ratzeburgii and P. subdepressus reported elsewhere. All of them are A+T rich, pericentromerically located, and with lengths of about 150 bp, either in the form of monomers or formed by more complex repeating units. A preliminary phylogenetic analysis of Palorus species using the 3' end of the mitochondrial Cytochrome Oxidase I gene shows that the five Palorus species have been diverging for a considerable amount of evolutionary time, with the pair P. ratzeburgii and P. genalis being the most closely related. Only these two taxa showed some similarity between their respective high-copy-number satellite sequences, while other satellites are mutually unrelated and might have originated independently. However, all the satellites have in common tertiary structure induced by intrinsic DNA curvature, a characteristic which is conserved within the genus. Palorus major satellites were previously detected in the genomes of congeneric species as low-copy-number clusters (Meštrovic et al., Mol. Biol. Evol. 15: 1062-1068. 1998). Given the divergences between the analysed species, the substitution rate deduced from high- and low-copy-number repeats is unexpectedly low. The presence of sequence-induced DNA curvature in all Palorus satellites and similar satellite DNAs in the species pair P. ratzeburgii and P. genalis suggest (i) that constraints are at the tertiary structure; and (ii) that the satellite DNA evolutionary turnover can be dependant on the history of the taxa under study, resulting in retention of similar satellites in related taxa.Key words: satellite DNA, evolution, mitochondrial cytochrome oxidase I, DNA curvature.

Satellite DNA in the ant Messor structor (Hymenoptera, Formicidae)

This paper is the first record of the satellite DNA of Formicidae. The satellite DNA of the ant Messor structor is organized in a tandem repeat of monomers of 79 bp. Like satellite DNAs of other insects, it is AT rich and presents direct and inverted internal repeats. Restriction analysis of the total DNA with methylation-sensitive enzymes strongly suggests that this DNA is undermethylated. The presence of this repetitive DNA in other species of the genus Messor is also tested.

A theoretical perusal of the satellite DNA curvature in tenebrionid beetles

The curvature patterns of seven satellite DNAs taken from beetles belonging to the family Tenebrionidae (Coleoptera) were modelled utilising a number of computer programs that describe and plot the curvature profiles of DNA. The theoretical analysis agreed with the experimentally observed curvature of most of these satellite DNAs, and its absence in Tribolium freemani and Tenebrio obscurus satellite I. In many cases, the tenebrionid satellite DNAs lack periodically repeated runs of phased-A-tracts, yet they represent a clear example of curved DNA. The macroscopic curvature of satellites from these closely related organisms confirmed that other sequence elements must be participating in the bending of these DNAs. Our modelling approaches are discussed, together with previous experimental results, in terms of the role played by DNA curvature in the organisation of satellite DNA and the tight compacting of heterochromatin.

Conservation of satellite DNA in species of the genus Pimelia (Tenebrionidae, Coleoptera)

Satellite DNA has been characterized in six allopatric species from the genus Pimelia: P. interjecta, P. integra, P. variolosa and P. baetica, inhabiting Iberian Peninsula, and P. elevata and P. criba, endemic to Balearic Islands Ibiza and Mallorca, respectively. All species show the presence of a single satellite DNA of a basic monomer length of 357 bp and A+T content of 69%, comprising a considerable amount of the genome (39%-45%, corresponding to about 4.5 x 10(5) copies per haploid genome). The sequence analysis of 22 cloned repeats reveals very high intra- and interspecific sequence similarity. Phylogenetic analysis separates the satellite sequences into two clusters, each comprising clones from three species exclusively. Within the clusters, satellite clones are not grouped species-specifically, except those of P. integra where species-diagnostic nt substitutions are detected with a pattern that could be produced by gene conversion. Such high sequence conservation could be related to preservation of satellite DNA curvature, resulting in a higher order helical structure, proposed to act as a specific protein binding domain.

Polymorphic curvature of satellite DNA in three subspecies of the beetle Pimelia sparsa

The curvature of the monomeric repeats of satellite DNAs from three subspecies of the beetle Pimelia sparsa (Coleoptera, Tenebrionidae) has been analysed. Evidence of curvature was inferred from their retarded migration in native polyacrylamide gels, which was confirmed by direct electron microscopy visualisation. Sequence-comparison analysis, which included sequence alignments and modelling studies, was used to reveal the patterns of local bending and curvature. The effects of the minor-groove-binding drugs distamycin and berenil on the curvature were assayed, and analysed with respect to their (A+T)-rich sequence preferences. Since our study deals with satellite DNAs from closely related organisms (three subspecies), we correlated the differences in sequence, and also the high similarity conserved in the (A+T)-rich regions, with the changes in the patterns of curvature.

Characterization of two abundant satellite DNAs from the mealworm Tenebrio obscurus

Two highly abundant satellite DNAs comprise 36% of the Tenebrio obscurus (Tenebrionidae, Coleoptera) genome. They are designated as satellite I and satellite II with the monomer length of 344 and 142 base pairs (bp), respectively. Both satellites differ in their nucleotide (nt) sequences, but the frequency of point mutations, well-conserved length of monomer variants, stretches of shared mutations characteristic for the process of gene conversion, and distribution of both satellites in regions of centromeric heterochromatin of all chromosomes indicate that the same evolutionary processes act on both of them with the same, or similar, rate. While satellite I shares no sequence similarity with any other known nt sequence, satellite II is 79.7% homologous with the highly abundant satellite from closely related Tenebrio molitor. Difference in the frequency of point mutations and absence of shared mutations indicating gene conversion strongly suggest that in these two closely related species mutational processes affecting satellite DNAs seem to be changed. Retarded electrophoretic mobility, due to sequence-induced curvature of DNA helix axis, was observed for T. obscurus satellite II, but not for satellite I. Although evolutionary processes act with different rates in T. obscurus and T. molitor satellites the monomer length and sequence-induced curvature are well preserved in both 142-bp satellites, as well as in, at the nt sequence level completely divergent, Palorus ratzeburgii (Tenebrionidae) satellite, indicating potential importance of these parameters in their evolution.

Analysis of divergence of Alphitobius diaperinus satellite DNA--roles of recombination, replication slippage and gene conversion

Satellite DNA is highly abundant in Alphitobius diaperinus (Tenebrionidae, Coleoptera), comprising 25% of the total genomic DNA. Sequence analysis reveals an average GC content of 50.8% and the presence of three different groups of satellite monomer variants, tH1, tH2 and tH3 with corresponding lengths of 123, 128 and 126 bp. Their mutual homologies range between 65 and 81%. Sequence comparison shows that the monomer variant tH2 has been formed by a recombination process between tH1 and tH3, which have a low average homology of only 65.15%. The longest stretch of 100% homology between the recombining units is 17 bp and is located 3' to the predicted recombination site. There is also an indication from sequence analysis that replication slippage and gene conversion play a part in the formation of satellite units and contribute to their divergence. The tH1, tH2 and tH3 monomer variants are organized in higher order repeating structures: a dimer, composed of tH1 and tH3, and a trimer containing tH1, tH2 and tH3 in series. The dimeric and trimeric repeat units furthermore create three higher order satellite subfamilies. Two of them contain either tandemly arranged dimers or trimers, while the third one is composed of both types of repeats, mutually interspersed.

Evolution of genetic redundancy for advanced players

An ever expanding database on the sequence organization and repetition of genic and non-genic components of nuclear and organelle genomes reveals that the vast majority of sequences are subject to one or other mechanism of DNA turnover (gene conversion, unequal crossing over, slippage, retrotransposition, transposition and others). Detailed studies, using novel methods of experimental detection and analytical procedures, show that such mechanisms can operate one on top of another and that wide variations in their unit lengths, biases, polarities and rates create bizarre and complex patterns of genetic redundancy. The ability of these mechanisms to operate both within and between chromosomes implies that realistic models of the evolutionary dynamics of redundancy, and of the potential interaction with natural selection in a sexual species, need to consider the diffusion of variant repeats across multiple chromosome lineages, in a population context. Recently, important advances in both experimental and analytical approaches have been made along these lines. There is increasing awareness that genetic redundancy and turnover induces a molecular co-evolution between functionally interacting genetic systems in order to maintain essential functions.

Localization of tandemly repeated DNA sequences in beetle chromosomes by fluorescent in situ hybridization

In situ hybridization to chromosomes and nuclei of Tenebrio molitor shows the massive presence of a species-specific satellite DNA in all chromosomes and six sites of rDNA in mitotic chromosomes. These sites are located in two autosomal pairs and in the X and Y chromosomes. In a related species, Misolampus goudoti, in which two different families of highly repetitive DNA have been previously characterized, one family is located in centromeric regions of all chromosomes with the exception of chromosome Y, while the other repeated DNA family is present both in centromeric and distal regions of all chromosomes. rRNA genes in this species are present in a medium-sized autosomal pair only. These results show that molecular cytogenetics can be applied to coleopteran chromosomes and open the way for a physical mapping of DNA sequences in these organisms. The results also provide insights into the type of meiotic association of the X and Y chromosomes in Coleoptera and the distribution of repeated DNAs within the genome of these insects.