Ribosomal DNA Evolution at the Population Level: Nucleotide Variation in Intergenic Spacer Arrays of Daphnia Pulex

Nucleotide variation was surveyed in 21 subrepeat arrays from the ribosomal DNA intergenic spacer of three Daphnia pulex populations. Eighteen of these arrays contained four subrepeats. Contrary to expectations, each of the four positions within the array had a different consensus sequence. However, gene conversion, involving sequences less than the length of a subrepeat, had occurred between subrepeats in different positions. Three arrays had more than four subrepeats and were undoubtedly generated by unequal crossing over between standard-length arrays. The data strongly suggested that most unequal exchanges between arrays are intrachromosomal and that they occur much less frequently than unequal exchanges at the level of the entire rDNA repeat. Strong associations among variants at different positions allowed the recognition of five groups of arrays, two of which were found in more than one population. Five of the seven individuals surveyed had arrays from more than one group. Analysis of the distribution of nucleotide variation suggested that the populations were quite divergent, a result that is concordant with previous surveys of allozyme and mitochondrial DNA variation. It was suggested that some of the subrepeat array types are quite old, at least predating the recolonization of pond habitats in the midwestern United States after the last glaciation.     

Human rDNA: evolutionary patterns within the genes and tandem arrays derived from multiple chromosomes

Human rDNA forms arrays on five chromosome pairs and is homogenized by concerted evolution through recombination and gene conversion (loci RNR1, RNR2, RNR3, RNR4, RNR5, OMIM: 180450). Homogenization is not perfect, however, so that it becomes possible to study its efficiency within genes, within arrays, and between arrays by measuring and comparing DNA sequence variation. Previous studies with randomly cloned genomic DNA fragments showed that different parts of the gene evolve at different rates but did not allow comparison of rDNA sequences derived from specific chromosomes. We have now cloned and sequenced rDNA fragments from specific acrocentric chromosomes to (1) study homogenization along the rDNA and (2) compare homogenization within chromosomes and between homologous and nonhomologous chromosomes. Our results show high homogeneity among regulatory and coding regions of rDNA on all chromosomes, a surprising homogeneity among adjacent distal non-rDNA sequences, and the existence of one to three very divergent intergenic spacer classes within each array.     

Patterns of variation in the intergenic spacers of ribosomal DNA in Drosophila melanogaster support a model for genetic exchanges during X-Y pairing

Detailed analysis of variation in intergenic spacer (IGS) and internal transcribed spacer (ITS) regions of rDNA drawn from natural populations of Drosophila melanogaster has revealed contrasting patterns of homogenization although both spacers are located in the same rDNA unit. On the basis of the role of IGS regions in X-Y chromosome pairing, we proposed a mechanism of single-strand exchanges at the IGS regions, which can explain the different evolutionary trajectories followed by the IGS and the ITS regions. Here, we provide data from the chromosomal distribution of selected IGS length variants, as well as the detailed internal structure of a large number of IGS regions obtained from specific X and Y chromosomes. The variability found in the different internal subrepeat regions of IGS regions isolated from X and Y chromosomes supports the proposed mechanism of genetic exchanges and suggests that only the "240" subrepeats are involved. The presence of a putative site for topoisomerase I at the 5' end of the 18S rRNA gene would allow for the exchange between X and Y chromosomes of some 240 subrepeats, the promoter, and the ETS region, leaving the rest of the rDNA unit to evolve along separate chromosomal lineages. The phenomenon of localized units (modules) of homogenization has implications for multigene family evolution in general.     

Structure,molecular evolution,and phylogenetic utility of the 5(') region of the external transcribed spacer of 18S-26S rDNA in Lessingia (Compositae,Astereae)

The 18S-26S nuclear rDNA external transcribed spacer (ETS) has recently gained attention as a region that is valuable in phylogenetic analyses of angiosperms primarily because it can supplement nucleotide variation from the widely used and generally shorter internal transcribed spacers (ITS-1 and ITS-2) and thereby improve phylogenetic resolution and clade support in rDNA trees. Subrepeated ETS sequences (often occurring in the 5(') region) can, however, create a challenge for systematists interested in using ETS sequence data for phylogeny reconstruction. We sequenced the 5(')ETS for members of Lessingia (Compositae, Astereae) and close relatives (26 taxa total) to characterize the subrepeat variation across a group of closely related plant lineages and to gain improved understanding of the structure, molecular evolution, and phylogenetic utility of the region. The 5(')ETS region of Lessingia and relatives varied in length from approximately 245 to 1009 bp due to the presence of a variable number of subrepeats (one to eight). We assessed homology of the subrepeats using phylogenetic analysis and concluded that only two of the subrepeats and a portion of a third ( approximately 282 bp in total) were orthologous across Lessingia and could be aligned with confidence and included in further analyses. When the partial 5(')ETS data were combined with 3(')ETS and ITS data in phylogenetic analyses, no additional resolution of relationships among taxa was obtained beyond that found from analysis of 3(')ETS + ITS sequences. Inferred patterns of concerted evolution indicate that homogenization is occurring at a faster rate in the 3(')ETS and ITS regions than in the 5(')ETS region. Additionally, homogenization appears to be acting within but not among subrepeats of the same rDNA array. We conclude that challenges in assessing subrepeat orthology across taxa greatly limit the utility of the 5(')ETS region for phylogenetic analyses among species of Lessingia.     

Extraordinary ribosomal spacer length heterogeneity in a neotyphodium endophyte hybrid: implications for concerted evolution.

An extraordinary level of length heterogeneity was found in the ribosomal DNA (rDNA) of an asexual hybrid Neotyphodium grass endophyte, isolate Lp1. This hybrid Neotyphodium endophyte is an interspecific hybrid between two grass endophytes, Neotyphodium lolii, and a sexual form, Epichlöe typhina, and the length heterogeneity was not found in either of these progenitor species. The length heterogeneity in the hybrid is localized to the intergenic spacer (IGS) and is the result of copy-number variation of a tandemly repeated subrepeat class within the IGS, the 111-/119-bp subrepeats. Copy number variation of this subrepeat class appears to be a consequence of mitotic unequal crossing over that occurs between these subrepeats. This implies that unequal crossing over plays a role in the concerted evolution of the whole rDNA. Changes in the pattern of IGS length variants occurred in just two rounds of single-spore purification. Analysis of the IGS length heterogeneity revealed features that are unexpected in a simple model of unequal crossing over. Potential refinements of the molecular details of unequal crossing over are presented, and we also discuss evidence for a combination of homogenization mechanisms that drive the concerted evolution of the Lp1 rDNA.     

Aspects of nonrandom turnover involved in the concerted evolution of intergenic spacers within the ribosomal DNA of Drosophila melanogaster

Polymerase chain reaction (PCR)-amplified, sequenced, and digitally typed intergenic spacers (IGSs) of the ribosomal (r)DNA in D. melanogaster reveal unexpected features of the mechanisms of turnover involved with the concerted evolution of the gene family. Characterization of the structure of three isolated IGS length variants reveals breakage hot spots within the 330-base-pair (bp) subrepeat array found in the spacers. Internal mapping of variant repeats within the 240-bp subrepeat array using a novel digital DNA typing procedure (minisatellite variant repeat [MVR]-PCR) shows an unexpected pattern of clustering of variant repeats. Each 240-bp subrepeat array consists of essentially two halves with the repeats in each half identified by specific mutations. This bipartite structure, observed in a cloned IGS unit, in the majority of genomic DNA of laboratory and wild flies and in PCR-amplified products, has been widely homogenized yet is not predicted by a model of unequal crossing over with randomly placed recombination breakpoints. Furthermore, wild populations contain large numbers of length variants in contrast to uniformly shared length variants in laboratory stocks. High numbers of length variants coupled to the observation of a homogenized bipartite structure of the 240-bp subrepeat array suggest that the unit of turnover and homogenization is smaller than the IGS and might involve gene conversion. The use of PCR for the structural analysis of members of the rDNA gene family coupled to digital DNA typing provides powerful new inroads into the mechanisms of DNA turnover affecting the course of molecular evolution in this family. Correspondence to: G. A. Dover     

Subrepeats within the BR1β repeat unit inChironomus pallidivittatus can be classified into different types depending on codon usage

Summary A new type of repeat unit was isolated from Balbiani ring 1 ofChironomus pallidivittatus and designated, BR1 repeat. It consists of a constant and a subrepeated part, like previously described units belonging to the core blocks of the BR genes. The subrepeated part contains 10-codon subrepeats with an arrangement similar to the subrepeats of the previously described BR2 gene. The present unit differs from earlier reported core units firstly in a much lower number of copies (about 15) per genome, which are tandemly arranged. Secondly, the number of subrepeats per BR1 repeat unit can show great variations. On the basis of the pattern of codon usage, three types of subrepeats can be distinguished. One type lies 5-proximal in the subrepeat array and consists of variable numbers of subrepeats almost identical at the nucleotide level. The last complete subrepeat represents another type, with consistent differences in codon usage as compared to subrepeats of the proximal type. Finally, there is an intermediate type represented by the subrepeat preceding the distal one. Here, codon characteristics from proximal and distal subrepeats are mixed in a patchy and irregular way. The evolution of the arrays can be understood either as being the result of subrepeat formation in two steps (occurring before and after amplification of whole repeat units) or as the result of a continuous process in which there is evidence for participation of gene conversion.     

Sequence comparison of distal and proximal ribosomal DNA arrays in rice (Oryza sativa L.) chromosome 9S and analysis of their flanking regions

Rice (Oryza sativa ssp. japonica cv. Nipponbare) harbors a ribosomal RNA gene (rDNA) cluster in the nucleolar-organizing region at the telomeric end of the short arm of chromosome 9. We isolated and sequenced two genomic clones carrying rice rDNA fragments from this region. The rice rDNA repeat units could be classified into three types based on length, which ranged from 7,928 to 8,934 bp. This variation was due to polymorphism in the number of 254-bp subrepeats in the intergenic spacer (IGS). Polymerase chain reaction (PCR) analysis suggested that the rDNA units in rice vary widely in length and that the copy number of the subrepeats in the IGS ranges from 1 to 12 in the rice genome. PCR and Southern blot analyses showed that most rDNA units have three intact and one truncated copies of the subrepeats in the IGS, and distal (telomere-side) rDNA units have more subrepeats than do proximal (centromere-side) ones. Both genomic clones we studied contained rDNA-flanking DNA sequences of either telomeric repeats (5′-TTTAGGG-3′) or a chromosome-specific region, suggesting that they were derived from the distal or proximal end, respectively, of the rDNA cluster. A similarity search indicated that retrotransposons appeared more frequently in a 500-kb portion of the proximal rDNA-flanking region than in other subtelomeric regions or sequenced regions of the genome. This study reveals the repetitive nature of the telomeric end of the short arm of chromosome 9, which consists of telomeric repeats, an rDNA array, and a retrotransposon-rich chromosomal region.Sequence accession numbers in DDBJ assigned for OSJNOa063K24 and OSJNBb0013K10 are AP009051 and AP008245, respectively.     

Nuclear ribosomal DNA internal transcribed spacer 1 (ITS1) in Picea (Pinaceae): sequence divergence and structure

The nrDNA ITS1 of Picea is 2747-3271 bp, the longest known of all plants. We obtained 24 cloned ITS1 sequences from six individuals of Picea glehnii, Picea mariana, Picea orientalis, and Picea rubens. Mean sequence divergence within these individuals (0.018+/-0.009) is more than half that between the species (0.031+/-0.011) and may be maintained against concerted evolution by separation of Picea 18S-26S rDNA repeats on multiple chromosomes. Picea ITS1 contains three subrepeats with a motif (5'-GGCCACCCTAGTC) that is conserved across Pinaceae. Two subrepeats are tandem, remote from the third, and more closely related and significantly more similar to one another than either is to the third subrepeat. This correlation between similarity and proximity may be the result of subrepeat duplication or concerted evolution within rDNA repeats. In inferred secondary structures, subrepeats generally form long hairpins, with a portion of the Pinaceae conserved motif in the terminal loop, and tandem subrepeats pair with one another over most of their length. Coalescence of ITS1 sequences occurs in P. orientalis but not in the other species.     

Homogenization of Tandemly Repeated Nucleotide Sequences by Distance-Dependent Nucleotide Sequence Conversion

Previous work revealed that recurrent mutations (= mutation occurring more than once) in the tandemly repeated arrays present in nontranscribed spacers (NTS) of ribosomal RNA genes (rDNA) are clustered, i.e., they most frequently occur in repeats with adjacent or alternate distribution. A possible explanation is that the likelihood of heteroduplex formation, a prerequisite of gene conversion, decreases with the distance between repeats. To test this possibility, evolution of an array of 11 initially homogeneous repeats was computer simulated using three models, two assuming that the likelihood of heteroduplex formation decreases with increasing distance between the repeats and one assuming that it is constant. Patterns of mutation distribution obtained in computer simulations were compared with the distribution of mutations found in the repeated arrays in the NTS of seven rDNA clones. The patterns of mutations generated by the models assuming that the likelihood of heteroduplex formation decreases as distance between the repeats increases agreed with the patterns observed in rDNA; the patterns generated by the model assuming that the likelihood is independent of distance between repeats disagreed with the patterns observed in the rDNA clones. The topology of the heteroduplex formed between DNA in adjacent repeats predicts that the most frequently occurring conversions in the NTS repeated arrays will be shorter than the length of the repeat. The topology of the heteroduplex also predicts that if the heteroduplex leads to crossing over a circular repeat is excised. It is speculated that the circle can transpose or can be amplified via rolling circle replication and subsequently transpose. It is also shown that homogenization of the NTS repeated arrays proceeds at different rates in different species.     

Balbiani ring DNA: Sequence comparisons and evolutionary history of a family of hierarchically repetitive protein-coding genes

Summary All known types of Balbiani ring (BR) gene consist of multiple, tandemly arranged, ca. 180 to 300-bp repeat units that can be divided into a constant region and a subrepeat region. The latter region includes short tandem subrepeats (SRs). Comparison of all available BR sequences using computer methods has enabled us (a) to define more precisely the constant and subrepeat regions, (b) to infer the evolutionary relationships among the various types of BR repeats, (c) to derive a consensus approximation of an ancestral sequence from a small segment of which the highly diverse present-day SRs may have originated, and (d) to detect an underlying substructure in the constant region, evident in the consensus but not in the present-day sequences and possibly corresponding to an original 39-bp DNA segment from which the extant, giant BR sequences may have evolved. We discuss the processes of reduplication, diversification, and homogenization within the hierarchically repetitive BR sequences as examples of how a simple DNA element may evolve into a diverse family of large, protein-coding genes.     

Phylogenetic analysis of the rDNA intergenic spacer subrepeats and its implication for the domestication history of foxtail millet, Setaria italica

We sequenced ribosomal DNA intergenic spacer subrepeats and their flanking regions of foxtail millet landraces from various regions in Europe and Asia and its wild ancestor to elucidate phylogenetic differentiation within each of types I–III found in our previous work and to elucidate relationships among these three types. Type I was classified into seven subtypes designated as Ia–Ig based on subrepeat sequences; C repeats downstream of those subrepeats are also polymorphic. Of these, subtypes Ia–Id and Ig were found in foxtail millet landraces. Subtypes Ia and Ib were distributed broadly throughout Asia and Europe. Subtype Ic was distributed in China, Korea and Japan. Subtype Id has a 20-bp deletion in subrepeat 3 and has a unique C repeat sequence. This subtype was found in a morphologically primitive landrace group from Afghanistan and northwestern Pakistan and differed greatly from other type I subtypes, implying that these landraces were domesticated independently. Subtypes Ig was found in a landrace from Pakistan and Ia and Ie–Ig were in six wild ancestor accessions. Type II was also highly polymorphic and four subtypes were found and designated as subtypes IIa–IId, but sequence analyses indicated type III as monomorphic. The present work indicates that type III should be classified as a subtype of type II (subtype IIe). Sequence polymorphism of subrepeats of types I–III indicated that subrepeats of subtype IIa are greatly divergent from others. Relationships among types I–III are much more complicated than anticipated based on previous RFLP work.     

Nucleotide sequence of intergenic and external transcribed spacers of rDNA of diploid wheat Triticum urartu Thum. ex Candil

The primary structure of intergenic non-transcribed and external transcribed spacers of rDNA of diploid wheat Triticum urartu, cloned in pTu3 plasmid 2402 b.p. long was determined. The intergenic non-transcribed rDNA spacer of Tr. urartu was shown to consist of 8 subrepeats with an average of 133 b.p. long, heterogeneous in length and nucleotide sequence. A number of repeated sequences was revealed within each subrepeat. While comparing nucleotide sequences of rDNA subrepeats of Tr. urartu and Tr. aestivum a high homology was found (up to 82%). A high similarity between these plant species was also found in the promoter region and in the external transcribed rDNA spacer. Suppression of the nucleolar organizer of 1A chromosome in the presence of 1B and 6B chromosomes of Tr. aestivum is supposed to be connected with the existence of a great number of subrepeats in the intergenic non-transcribed rDNA spacer of B genome donors in polyploid wheat species of turgidum-aestivum row.     

Organization of rRNA genes in various tribes of plants in the cereal family (Poaceae barnh.)

The organization of 18S and 26S rRNA in 14 plant species, belonging to 8 tribes of the cereal family was studied. In rDNA of all the cereals studied, except maize and reed, the similar character of localization of nucleotide sequences, recognized by restrictases BamHI and EcoRI in 18S and 26S rRNA genes was revealed. The structural organization of rDNA of sainfoin (Papilionaceae) was shown to differ from genes, coding for high molecular rRNA in cereals. The primary structure of subrepeat of non-transcribed rDNA spacer of diploid wheat Tr. urartu, consisting of 132 base pairs was determined. The given subrepeat was hybridized with BamHI-fragments of DNA from cereals and sainfoin. It is shown to hybridize with rDNA of all the cereals studied, and it hardly hybridizes with rDNA of maize and sorghum, but doesn't hybridize with rDNA of sainfoin. The conclusion is made that the size polymorphism of restriction fragments in the coding rDNA region and the level of similarity of subrepeats of rDNA of the non-transcribed spacer may help to reveal the phylogenetic affinity of plants, belonging to different tribes within one family.     

Structure and variation in ribosomal RNA genes of pea

Summary A complete ribosomal DNA (rDNA) repeat unit has been cloned from the genome of Pisum sativum (garden pea) and used to construct a map containing a total of 58 cleavage sites for 23 different restriction enzymes. Regions encoding 18s and 25s ribosomal RNA (rRNA) were identified by R-loop analysis. A 180 bp sequence element is repeated eight times in the intergenic nontranscribed spacer (NTS) region, as defined by eight evenly spaced RsaI cleavage sites. Sequence heterogeneity among these elements (subrepeats) is indicated by the presence of an NcoI site within the five RsaI subrepeats distal to the 25s rRNA gene but not in the three subrepeats proximal to this gene, and also by the presence of an additional RsaI cleavage site in one subrepeat.The approximately 4000 copies of the rDNA repeat in the pea nuclear genome show considerable heterogeneity with respect to the length of the NTS region, and differences are also frequently observed between different genotypes. In both cases the length variation appears to be due primarily to differences in the number of subrepeat elements.Comparison of rDNA restriction maps for two pea genotypes separated for hundreds or perhaps thousands of generations reveals that they contain many rDNA identical repeat units. This data is consistent with the view that new rDNA variants are fixed only infrequently in the evolution of a species.Differences also exist between the rDNA repeats of a single genotype with respect to the degree of base modification at certain restriction sites. A large number of sites known to exist in the pea rDNA clone are not cleaved at all in genomic rDNA, or are cleaved in only some copies of the rDNA repeat. We believe these examples of incomplete cleavage results mostly from methylation, although it is difficult to rule out the possibility of sequence variation in all cases. Most putative modifications are best interpreted in terms of cytosine methylation in CG and CXG sequences, but at least one example is more consistent with adenine methylation.We also have constructed a more detailed restriction map of the wheat rDNA clone pTA71 and present a comparison of this map to our map of pea, pumpkin, and wheat in order to assess the amount of useful evolutionary information that can be obtained by comparison of such maps.     

Complex telomere-associated repeat units in members of the genus Chironomus evolve from sequences similar to simple telomeric repeats.

The dipteran Chironomus tentans has complex tandemly repeated 350-bp DNA sequences at or near the chromosome ends. As in Drosophila melanogaster, short simple repeats with cytosines and guanines in different strands have never been observed. We were therefore interested in learning whether the Chironomus repeats could have evolved from simple sequence telomeric DNA, which might suggest that they constitute a functional equivalent. We screened for repeat units with evolutionarily ancient features within the tandem arrays and recovered two clones with a less-evolved structure. Sequence analysis reveals that the present-day 350-bp unit probably evolved from a simpler 165-bp unit through the acquisition of transposed sequences. The 165-bp unit contains DNA with a highly biased distribution of cytosine and guanine between the two strands, although with the ratios inverted in two minor parts of the repeat. It is largely built up of short degenerate subrepeats for which most of the sequence can be reconstructed. The consensus for the subrepeat sequence is similar to the simple telomeric repeat sequences of several kinds of eukaryotes. We propose that the present-day unit has evolved from telomeric, simple sequence, asymmetric DNA from which it has retained some original sequence features and possibly functions.     

Sequence subfamilies of satellite repeats related to rDNA intergenic spacer are differentially amplified on<Emphasis Type="Italic"> Vicia sativa</Emphasis> chromosomes

We cloned and sequenced the Vicia sativa 25S-18S rDNA intergenic spacer (IGS) and the satellite repeat S12, thought to be related to the spacer sequence. The spacer was shown to contain three types of subrepeats (A, B, and C) with monomers of 173 bp (A), 10 bp (B), and 66 bp (C), separated by unique or partially duplicated sequences. Two spacer variants were detected in V. sativa that differed in length (2990 and 3168 bp) owing to an extra copy of the subrepeat A. The A subrepeats were also shown to be highly homologous to the satellite repeat S12, which is located in large clusters on chromosomes 4, 5, and 6, and is not associated with the rDNA loci. Sequencing of additional S12 clones retrieved from a shotgun genomic library allowed definition of three subfamilies of this repeat based on minor differences in their nucleotide sequences. Two of these subfamilies could be discriminated from the rest of the S12 sequences as well as from the IGS A subrepeats using specific oligonucleotide primers that labeled only a subset of the S12 loci when used in the primed in situ DNA labeling (PRINS) reaction on mitotic chromosomes. These experiments showed that, in spite of the high overall similarity of the IGS A subrepeats and the S12 satellite repeats, there are S12 subfamilies that are divergent from the common consensus and are present at only some of the chromosomes containing the S12 loci. Thus, the subfamilies may have evolved at these loci following the spreading of the A subrepeats from the IGS to genomic regions outside the rDNA clusters.Electronic Supplementary Material Supplementary material is available in the online version of this article at Accession numbers: GenBank AY234364–AY234374. The monomer sequences and additional information about the family of IGS-like repeat S12 will also appear in the PlantSat database (Macas et al. 2002, ) under Accession name Vicia_sativa_IGS-like