Molecular relationships and phylogeny in a community of myxosporeans and actinosporeans based on their 18S rDNA sequences

The community of myxosporeans and actinosporeans inhabiting a typical Scottish highland stream and the outflow area of an adjacent salmon hatchery was analysed on the basis of their 18S rDNA sequences. Nine myxosporeans belonging to the genera Sphaerospora, Chloromyxum, Zschokkella, Myxidium, Hoferellus and Myxobilatus were identified from mature spores in different organs of the fish species present. Twelve actinosporean types belonging to the collective groups of neoactinomyxum, aurantiactinomyxon, raabeia, echinactinomyxon and synactinomyxon were found to be released from oligochaete worms collected from sediments. Twenty of the 21 sequences obtained from these myxozoans are new entries to the myxozoan database, and the genera Chloromyxum, Hoferellus and Myxobilatus were entered for the first time. Study of the molecular relationships between the different taxa and with other myxozoan sequences available showed that the myxosporeans inhabiting the urinary system clearly cluster together, independently of host species or spore morphology. However, the sequences of the two Sphaerospora species encountered show considerable differences from other members of this group and all other freshwater myxosporeans, and they were found to occupy an ancestral marine position. Three actinosporeans, i.e. Neoactinomyxum eiseniellae, Aurantiactinomyxon pavinsis and Raabeia 'type 3' were found to represent alternate life cycle stages of Chloromyxum sp., Chloromyxum truttae and Myxidium truttae, respectively (approximately 1400 identical base pairs each). Three other actinosporeans encountered (two echinactinomyxon and one raabeia type) showed over 92% sequence identity with myxosporeans from GenBank, whereas all other actinosporeans formed a closely related group devoid of any known myxosporeans.     

Sequencing and analysis of the internal transcribed spacers (ITSs) and coding regions in the EcoR I fragment of the ribosomal DNA of the Japanese pond frog Rana nigromaculata

The rDNA of eukaryotic organisms is transcribed as the 40S-45S rRNA precursor, and this precursor contains the following segments: 5' - ETS - 18S rRNA - ITS 1 - 5.8S rRNA - ITS 2 - 28S rRNA - 3'. In amphibians, the nucleotide sequences of the rRNA precursor have been completely determined in only two species of Xenopus. In the other amphibian species investigated so far, only the short nucleotide sequences of some rDNA fragments have been reported. We obtained a genomic clone containing the rDNA precursor from the Japanese pond frog Rana nigromaculata and analyzed its nucleotide sequence. The cloned genomic fragment was 4,806 bp long and included the 3'-terminus of 18S rRNA, ITS 1, 5.8S rRNA, ITS 2, and a long portion of 28S rRNA. A comparison of nucleotide sequences among Rana, the two species of Xenopus, and human revealed the following: (1) The 3'-terminus of 18S rRNA and the complete 5.8S rRNA were highly conserved among these four taxa. (2) The regions corresponding to the stem and loop of the secondary structure in 28S rRNA were conserved between Xenopus and Rana, but the rate of substitutions in the loop was higher than that in the stem. Many of the human loop regions had large insertions not seen in amphibians. (3) Two ITS regions had highly diverged sequences that made it difficult to compare the sequences not only between human and frogs, but also between Xenopus and Rana. (4) The short tracts in the ITS regions were strictly conserved between the two Xenopus species, and there was a corresponding sequence for Rana. Our data on the nucleotide sequence of the rRNA precursor from the Japanese pond frog Rana nigromaculata were used to examine the potential usefulness of the rRNA genes and ITS regions for evolutionary studies on frogs, because the rRNA precursor contains both highly conserved regions and rapidly evolving regions.     

Recent Advances in Our Knowledge of the Myxozoa

In the last few years two factors have helped to significantly advance our understanding of the Myxozoa. First, the phenomenal increase in fin fish aquaculture in the 1990s has lead to the increased importance of these parasites; in turn this has lead to intensified research efforts, which have increased knowledge of the development, diagnosis. and pathogenesis of myxozoans. The hallmark discovery in the 1980s that the life cycle of Myxobolus cerebralis requires development of an actinosporean stage in the oligochaete. Tubifex tubifex, led to the elucidation of the life cycles of several other myxozoans. Also, the life cycle and taxonomy of the enigmatic PKX myxozoan has been resolved: it is the alternate stage of the unusual myxozoan, Tetracapsula bryosalmonae, from bryozoans. The 18S rDNA gene of many species has been sequenced, and here we add 22 new sequences to the data set. Phylogenetic analyses using all these sequences indicate that: 1) the Myxozoa are closely related to Cnidaria (also supported by morphological data); 2) marine taxa at the genus level branch separately from genera that usually infect freshwater fishes; 3) taxa cluster more by development and tissue location than by spore morphology; 4) the tetracapsulids branched off early in myxozoan evolution, perhaps reflected by their having bryozoan, rather than annelid hosts; 5) the morphology of actinosporeans offers little information for determining their myxosporean counterparts (assuming that they exist); and 6) the marine actinosporeans from Australia appear to form a clade within the platysporinid myxosporeans. Ribosomal DNA sequences have also enabled development of diagnostic tests for myxozoans. PCR and in situ hybridisation tests based on rDNA sequences have been developed for Myxobolus cerebralis, Ceratomyxa shasta, Kudoa spp., and Tetracapsula bryosalmonae (PKX). Lectin-based and antibody tests have also been developed for certain myxozoans, such as PKX and C. shasta. We also review important diseases caused by myxozoans, which are emerging or re-emerging. Epizootics of whirling disease in wild rainbow trout (Oncorhynchus mykiss) have recently been reported throughout the Rocky Mountain states of the USA. With a dramatic increase in aquaculture of fishes using marine netpens, several marine myxozoans have been recognized or elevated in status as pathological agents. Kudoa thyrsites infections have caused severe post-harvest myoliquefaction in pen-reared Atlantic salmon (Salmo salar), and Ceratomyxa spp., Sphaerospora spp., and Myxidium leei cause disease in pen-reared sea bass (Dicentrarchus labrax) and sea bream species (family Sparidae) in Mediterranean countries.     

The 18S rRNA from Odontophrynus americanus 2n and 4n (Amphibia, Anura) reveals unusual extra sequences in the variable region V2.

The nucleotide sequence of the rDNA 18S region isolated from diploid and tetraploid species of the amphibian Odontophrynus americanus was determined and used to predict the secondary structure of the corresponding 18S rRNA molecules. Comparison of the primary and secondary structures for the 2n and 4n species confirmed that these species are very closely related. Only three nucleotide substitutions were observed, accounting for 99% identity between the 18S sequences, whereas several changes were detected by comparison with the Xenopus laevis 18S sequence (96% identity). Most changes were located in highly variable regions of the molecule. A noticeable feature of the Odontophrynus 18S rRNA was the presence of unusual extra sequences in the V2 region, between helices 9 and 11. These extra sequences do not fit the model for secondary structure predicted for vertebrate 18S rRNA.     

Structural constraints in expansion segments from a midge 26S rDNA

DNA sequences representing approximately 40% of the large-subunit rRNA gene from the lower dipteran Chironomus thummi were analyzed. Once aligned with their Drosophila counterparts, sequence and base content comparisons were carried out. Sequence identity was found to be high overall, except for six regions that displayed a local bias in nucleotide composition toward AT. These regions were identified as expansion segments D3, D4, D5, D6, D7a, and D12. Besides base sequence divergence, differences in length were observed between the respective variable domains of the two species, particularly for D7a. Prediction of secondary structure showed that the folding of the Chironomus expansion segments analyzed is in agreement with the general patterns proposed for eukaryotic LSU rRNA. The comparison with Drosophila revealed also that the Chironomus secondary structures of the variable domains are supported by multiple compensatory substitutions or even compensatory insertions. Chironomus D7a displayed an unusual structural feature with respect to the insect D7a models that have been inferred up to now. The structural constraint observed in the expansion segments of Diptera so distantly related as midges and Drosophila suggests that these regions contribute to some functional role. Concerning the D7a of insects so far analyzed, there can be, in addition to a conserved secondary structure, a nucleotide composition constraint that might be important for the process giving rise to the alpha and beta halves of the 26S rRNA. Correspondence to: E. Gorab     

Phylogenetic position of Sphaerospora testicularis and Latyspora scomberomori n. gen. n. sp. (Myxozoa) within the marine urinary clade

An amendment of the family Sinuolineidae (Myxosporea) is proposed in order to include a newly described genus Latyspora n. gen. The type species Latyspora scomberomori n. gen. n. sp. is a coelozoic parasite in the kidney tubules of Scomberomorus guttatus. In addition to the morphological and molecular characterization of L. scomberomori n. gen. n. sp., we also present novel SSU rDNA data on Sphaerospora testicularis, a serious parasite of Dicentrarchus labrax. Performed phylogenetic analyses revealed that both species cluster within the marine urinary clade encompassing the representatives with a shared insertion within their V4 SSU rRNA region and grouping according to the shape of their spores' sutural line and their similar tissue tropism in the host. Sphaerospora testicularis is the closest relative to Parvicapsula minibicornis within the Parvicapsula subclade and L. scomberomori n. gen. n. sp. is the basal species of the Zschokkella subclade. The phylogenetic position of S. testicularis, outwith the basal Sphaerospora sensu stricto clade, and its morphology suggest it being a non-typical Sphaerospora. The sequence data provided on S. testicularis can help in future revisions of the strongly polyphyletic genus Sphaerospora. We recommend re-sequencing of several sphaerosporids as an essential step before such taxonomic changes are accomplished.     

Complete sequences of the rRNA genes of Drosophila melanogaster

In this, the first of three papers, we present the sequence of the ribosomal RNA (rRNA) genes of Drosophila melanogaster. The gene regions of D. melanogaster rDNA encode four individual rRNAs: 18S (1,995 nt), 5.8S (123 nt), 2S (30 nt), and 28S (3,945 nt). The ribosomal DNA (rDNA) repeat of D. melanogaster is AT rich (65.9% overall), with the spacers being particularly AT rich. Analysis of DNA simplicity reveals that, in contrast to the intergenic spacer (IGS) and the external transcribed spacer (ETS), most of the rRNA gene regions have been refractory to the action of slippage-like events, with the exception of the 28S rRNA gene expansion segments. It would seem that the 28S rRNA can accommodate the products of slippage-like events without loss of activity. In the following two papers we analyze the effects of sequence divergence on the evolution of (1) the 28S gene "expansion segments" and (2) the 28S and 18S rRNA secondary structures among eukaryotic species, respectively. Our detailed analyses reveal, in addition to unequal crossing-over, (1) the involvement of slippage and biased mutation in the evolution of the rDNA multigene family and (2) the molecular coevolution of both expansion segments and the nucleotides involved with compensatory changes required to maintain secondary structures of RNA.      

Zschokkella hildae Auerbach, 1910: Phylogenetic position,morphology, and location in cultured Atlantic cod

The myxozoan Zschokkella hildae Auerbach, 1910, was detected with a prevalence of 100% in cultured Atlantic cod, Gadus morhua L. aged 1+ from a culture facility on the west coast of Scotland. Sporogonic stages of Z. hildae, plasmodia producing 2–5 mature spores, were located predominantly in the collecting ducts and ureters of the kidney, and spores were present in the urine collected from the bladder. Less frequently, plasmodia were detected in the interstitial tissue of the kidney. The parasite prevalence in cultured fish was considerably higher than reported in wild fish but no obvious signs of pathology were detected. SSU rDNA sequencing and phylogenetic analysis showed that Z. hildae is closely related to a Sinuolinea sp. from the urinary system of turbot, Psetta maxima (L.), and that these two species, together with other myxozoans from the urinary system of marine fish cluster together in a sub-clade of the recognised marine clade of myxozoans. This sub-clade is characterised by a specific linear expansion segment, helix E23_15 in the secondary structure of variable region V4 of the SSU rDNA. Z. hildae and Sinuolinea sp. show extraordinary large linear expansion segment in both V4 and V7 and an important number of complementary base changes in the conservative regions of the SSU rDNA, indicating considerable evolutionary changes in the SSU rDNA of these species when compared with other myxozoans from the marine environment.     

Hypervariable regions (V1–V9) of the dinoflagellate 18S rRNA using a large dataset for marker considerations

DNA sequencing methods have been used for the molecular taxonomic discrimination of dinoflagellate protists, particularly using partial 18S rRNA sequences. This study evaluated the taxonomic discrimination power of rRNA gene hypervariable regions (V1 to V9) in dinoflagellates from a large dataset. These included 77 dinoflagellate species (9 orders, 17 families, 40 genera). The complete 18S rRNA sequences of the dinoflagellates ranged from 1,787 to 1,813?bp in length, and consisted of eight V regions with a total combined length of 678 to 699?bp. Regions longer than 100?bp were recoded for V2, V4, and V8 regions; high nucleotide divergences were detected in V1, V2, and V4 regions. Statistic tests showed that the divergences of individual V regions were significantly different (t-test, P?<?0.05) compared with the complete 18S rRNA. The V2 region showed the highest score (83.5%) for PI sites. Moreover, intra-genus DNA similarities of the V2 were considerably low (<93%). Neighbor-joining analyses showed that phylogenetic resolution in the V2–V4 region was 1.32-fold higher than that of the complete 18S rRNA. These results demonstrate that V2 has the highest taxonomic resolving power within the 18S rRNA gene of dinoflagellates, suggesting the V2 and adjacent regions (e.g., V1 to V4) may be the best for marker considerations.     

Comparison of primary and secondary 26S rRNA structures in twoTetrahymena species: Evidence for a strong evolutionary and structural constraint in expansion segments

Summary We have determined the nucleotide sequence of the 26S large subunit (LSU) rRNA genes for twoTetrahymena species,T. thermophila andT. pyriformis. The inferred rRNA sequences are presented in their most probable secondary structures based on compensatory mutations, energy, and conservation criteria. The majority of the nucleotide changes between the twoTetrahymena LSU rRNAs and the positions of a relatively large deletion and of the processing cleavage sites resulting in the generation of the hidden break are all located within the so-called divergent domains or expansion segments. These are regions within the common core of secondary structure where expansions have taken place during the evolution of the rRNA of higher eukaryotes.The dispensable nature of some of the expansion segments has been taken as evidence of their non-functionality. However, our data show that a considerable selective constraint has operated to presesrve the secondary structure of these segments. Especially in the case of the D2 and D8 segments, the presence of a considerable number of compensatory base changes suggests that the secondary structure of these regions is of functional importance. Alternatively, these expansion segments may have maintained characteristic folding patterns because only such structures are being tolerated within otherwise functionally important regions.     

Mouse rDNA: sequences and evolutionary analysis of spacer and mature RNA regions.

Two regions of mouse rDNA were sequenced. One contained the last 323 nucleotides of the external transcribed spacer and the first 595 nucleotides of 18S rRNA; the other spanned the entire internal transcribed spacer and included the 3' end of 18S rRNA, 5.8S rRNA, and the 5' end of 28S rRNA. The mature rRNA sequences are very highly conserved from yeast to mouse (unit evolutionary period, the time required for a 1% divergence of sequence, was 30 X 10(6) to 100 X 10(6) years). In 18S rRNA, at least some of the evolutionary expansion and increase in G + C content is due to a progressive accretion of discrete G + C-rich insertions. Spacer sequence comparisons between mouse and rat rRNA reveal much more extensive and frequent insertions and substitutions of G + C-rich segments. As a result, spacers conserve overall G + C richness but not sequence (UEP, 0.3 X 10(6) years) or specific base-paired stems. Although no stems analogous to those bracketing 16S and 23S rRNA in Escherichia coli pre-rRNA are evident, certain features of the spacer regions flanking eucaryotic mature rRNAs are conserved and could be involved in rRNA processing or ribosome formation. These conserved regions include some short homologous sequence patterns and closely spaced direct repeats.     

Structural alignment of 18S and 28S rDNA sequences provides insights into phylogeny of Phytophaga (Coleoptera: Curculionoidea and Chrysomeloidea)

We performed a comparative study of partial rDNA sequences from a variety of Coleoptera taxa to construct an annotated alignment based on secondary structure information, which in turn, provides improved rRNA structure models useful for phylogenetic reconstruction. Subsequent phylogenetic analysis was performed to test monophyly and interfamilial relationships of the megadiverse plant feeding beetle group known as ‘Phytophaga’ (Curculionoidea and Chrysomeloidea), as well as to discover their closest relatives among the Cucujiformia. Parsimony and Bayesian analyses were performed based on the structural alignment of segments of 18S rRNA (variable regions V4‐V5, V7‐V9) and 28S rRNA (expansion segment D2). A total of 104 terminal taxa of Coleoptera were included: 96 species of Cucujiformia beetles, representing the families and most ‘subfamilies’ of weevils and chrysomeloids (Phytophaga), as well as several families of Cleroidea, Tenebrionoidea and Cucujoidea, and eight outgroups from three other polyphagan series: Scarabaeiformia, Elateriformia and Bostrichiformia. The results from the different methods of analysis agree — recovering the monophyly of the ‘Phytophaga’, including Curculionoidea and Chrysomeloidea as sister groups. The curculionoid and chrysomeloid phylogeny recovered from the aligned 18S and 28S rDNA segments, which is independent of morphological data, is in agreement with recent hypotheses or concepts based on morphological evidence, particularly with respect to familial relationships. Our results provide clues about the evolutionary origin of the phytophagan beetles within the megaclade Cucujiformia, suggesting that the sister group of ‘Curculionoidea + Chrysomeloidea’ is a clade of the ‘Cucujoidea’, represented in this study by species in Boganiidae, Erotylidae, Nitidulidae, Cucujidae and Silvanidae. The Coccinellidae and Endomychidae are not grouped with the latter, and the remaining terminal taxa are nested in Tenebrionoidea and Cleroidea. We propose that the combination of structurally aligned ribosomal RNA gene regions 18S (V4‐V5, V7‐V9) and 28S (D2) are useful in testing monophyly and resolving relationships among beetle superfamilies and families.     

Nucleotide sequence of Citrus limon 26S rRNA gene and secondary structure model of its RNA

The complete nucleotide sequence of Citrus limon 26S rDNA has been determined. The sequence has been aligned with large ribosomal RNA (L-rRNA) sequences of Escherichia coli, Saccharomyces cerevisiae and Oryza sativa. Nine extensive expansion segments in dicot 26S rRNA relative to E. coli 23S rRNA have been identified and compared with analogous segments of monocot, yeast, amphibian and human L-rRNAs. A secondary structure model for lemon 26S rRNA has been derived based on the refined model of E. coli 23S rRNA. It has been compared with other eukaryotic L-rRNAs models in terms of location of functionally important regions. Origin and evolution of L-rRNA expansion segments are discussed.     

A secondary structural model of the 28S rRNA expansion segments D2 and D3 for Chalcidoid wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.     

Analysis of the Primary Sequence and Secondary Structure of the Unusually Long SSU rRNA of the Soil Bug, Armadillidium vulgare

The complete nucleotide sequence of the SSU rRNA gene from the soil bug, Armadillidium vulgare (Crustacea, Isopoda), was determined. It is 3214 bp long, with a GC content of 56.3%. It is not only the longest SSU rRNA gene among Crustacea but also longer than any other SSU rRNA gene except that of the strepsipteran insect, Xenos vesparum (3316 bp). The unusually long sequence of this species is explained by the long sequences of variable regions V4 and V7, which make up more than half of the total length. RT-PCR analysis of these two regions showed that the long sequences also exist in the mature rRNA and sequence simplicity analysis revealed the presence of slippage motifs in these two regions. The putative secondary structure of the rRNA is typical for eukaryotes except for the length and shape variations of the V2, V4, V7, and V9 regions. Each of the V2, V4, and V7 regions was elongated, while the V9 region was shortened. In V2, two bulges, located between helix 8 and helix 9 and between helix 9 and helix 10, were elongated. In V4, stem E23-3 was dramatically expanded, with several small branched stems. In V7, stem 43 was branched and expanded. Comparisons with the unusually long SSU rRNAs of other organisms imply that the increase in total length of SSU rRNA is due mainly to expansion in the V4 and V7 regions. Received: 2 March 1999 / Accepted: 22 July 1999     

Probing the structure of mouse Ehrlich ascites cell 5.8S, 18S and 28S ribosomal RNA in situ.

The secondary structure of mouse Ehrlich ascites 18S, 5.8S and 28S ribosomal RNA in situ was investigated by chemical modification using dimethyl sulphate and 1-cyclohexyl-3-(morpholinoethyl) carbodiimide metho-p-toluene sulphonate. These reagents specifically modify unpaired bases in the RNA. The reactive bases were localized by primer extension followed by gel electrophoresis. The three rRNA species were equally accessible for modification i.e. approximately 10% of the nucleotides were reactive. The experimental data support the theoretical secondary structure models proposed for 18S and 5.8/28S rRNA as almost all modified bases were located in putative single-strand regions of the rRNAs or in helical regions that could be expected to undergo dynamic breathing. However, deviations from the suggested models were found in both 18S and 28S rRNA. In 18S rRNA some putative helices in the 5'-domain were extensively modified by the single-strand specific reagents as was one of the suggested helices in domain III of 28S rRNA. Of the four eukaryote specific expansion segments present in mouse Ehrlich ascites cell 28S rRNA, segments I and III were only partly available for modification while segments II and IV showed average to high modification.     

Secondary structure prediction for complete rDNA sequences (18S, 5.8S,and 28S rDNA) of Demodex folliculorum, and comparison of divergent domains structures across Acari

According to base pairing, the rRNA folds into corresponding secondary structures, which contain additional phylogenetic information. On the basis of sequencing for complete rDNA sequences (18S, ITS1, 5.8S, ITS2 and 28S rDNA) of Demodex, we predicted the secondary structure of the complete rDNA sequence (18S, 5.8S, and 28S rDNA) of Demodex folliculorum, which was in concordance with that of the main arthropod lineages in past studies. And together with the sequence data from GenBank, we also predicted the secondary structures of divergent domains in SSU rRNA of 51 species and in LSU rRNA of 43 species from four superfamilies in Acari (Cheyletoidea, Tetranychoidea, Analgoidea and Ixodoidea). The multiple alignment among the four superfamilies in Acari showed that, insertions from Tetranychoidea SSU rRNA formed two newly proposed helixes, and helix c3-2b of LSU rRNA was absent in Demodex (Cheyletoidea) taxa. Generally speaking, LSU rRNA presented more remarkable differences than SSU rRNA did, mainly in D2, D3, D5, D7a, D7b, D8 and D10.     

Molecular coevolution among cryptically simple expansion segments of eukaryotic 26S/28S rRNAs

The set of "expansion segments" of any eukaryotic 26S/28S ribosomal RNA (rRNA) gene is responsible for the bulk of the difference in length between the prokaryotic 23S rRNA gene and the eukaryotic 26S/28S rRNA gene. The expansion segments are also responsible for interspecific fluctuations in length during eukaryotic evolution. They show a consistent bias in base composition in any species; for example, they are AT rich in Drosophila melanogaster and GC rich in vertebrate species. Dot-matrix comparisons of sets of expansion segments reveal high similarities between members of a set within any 28S rRNA gene of a species, in contrast to the little or spurious similarity that exists between sets of expansion segments from distantly related species. Similarities among members of a set of expansion segments within any 28S rRNA gene cannot be accounted for by their base-compositional bias alone. In contrast, no significant similarity exists within a set of "core" segments (regions between expansion segments) of any 28S rRNA gene, although core segments are conserved between species. The set of expansion segments of a 26S/28S gene is coevolving as a unit in each species, at the same time as the family of 28S rRNA genes, as a whole, is undergoing continual homogenization, making all sets of expansion segments from all ribosomal DNA (rDNA) arrays in a species similar in sequence. Analysis of DNA simplicity of 26S/28S rRNA genes shows a direct correlation between significantly high relative simplicity factors (RSFs) and sequence similarity among a set of expansion segments. A similar correlation exists between RSF values, overall rDNA lengths, and the lengths of individual expansion segments. Such correlations suggest that most length fluctuations reflect the gain and loss of simple sequence motifs by slippage-like mechanisms. We discuss the molecular coevolution of expansion segments, which takes place against a background of slippage-like and unequal crossing-over mechanisms of turnover that are responsible for the accumulation of interspecific differences in rDNA sequences.