Gene Conversion Drives Within Genic Sequences: Concerted Evolution of Ribosomal RNA Genes in Bacteria and Archaea

Multiple copies of a given ribosomal RNA gene family undergo concerted evolution such that sequences of all gene copies are virtually identical within a species although they diverge normally between species. In eukaryotes, gene conversion and unequal crossing over are the proposed mechanisms for concerted evolution of tandemly repeated sequences, whereas dispersed genes are homogenized by gene conversion. However, the homogenization mechanisms for multiple-copy, normally dispersed, prokaryotic rRNA genes are not well understood. Here we compared the sequences of multiple paralogous rRNA genes within a genome in 12 prokaryotic organisms that have multiple copies of the rRNA genes. Within a genome, putative sequence conversion tracts were found throughout the entire length of each individual rRNA genes and their immediate flanks. Individual conversion events convert only a short sequence tract, and the conversion partners can be any paralogous genes within the genome. Interestingly, the genic sequences undergo much slower divergence than their flanking sequences. Moreover, genomic context and operon organization do not affect rRNA gene homogenization. Thus, gene conversion underlies concerted evolution of bacterial rRNA genes, which normally occurs within genic sequences, and homogenization of flanking regions may result from co-conversion with the genic sequence. Received: 31 March 2000 / Accepted: 15 June 2000     

Multiple rRNA variants in a single spore of the microsporidian Nosema bombi

To understand the source of the multiple DNA sequence variants of Nosema bombi ribosomal RNA (rRNA) found in a single bumble bee host, we PCR amplified, cloned, and sequenced the partial rRNA gene from 125 clones, which were derived from four out of 46 spores individually isolated from a single host by laser microdissection. At least two rRNA variants, characterized by either (GTTT)(2) or (GTTT)(3) repeat units within the internal transcribed spacer (ITS) region, were found per spore in approximately equal proportions, variants which were also found in approximately equal proportions in 55 clones of the two DNA extracts of multiple spores from the same host. Firstly, we demonstrate for the first time that DNA sequences can be obtained from single-binucleate microsporidia. Secondly, it appears that concerted evolution has not homogenized the sequences of all rRNA copies within a single N. bombi spore or even within a single nucleus. We thereby demonstrate unequivocally that two or more rRNA sequence variants exist per N. bombi spore, and urge caution in the use of multicopy rRNA genes for population genetic and phylogenetic analysis of this and other Microsporidia unless homologous copies can be reliably typed.     

rRNA gene organization in the Lyme disease spirochete, Borrelia burgdorferi.

Lyme disease is the most common vector-borne disease in the United States. The causative agent is the spirochete Borrelia burgdorferi. The copy number and organization of the genes encoding the rRNAs of this organism were determined. There is a single gene for 16S rRNA and two copies each of the 23S rRNA and 5S rRNA genes. All of the genes are located within a chromosomal fragment of approximately 9.5 to 10.0 kb. The 23S and 5S rRNA genes are tandemly duplicated in the order 23S-5S-23S-5S and are apparently not linked to the 16S rRNA gene, which is situated over 2 kb upstream from the 23S-5S duplication. The individual copies of the 23S-5S duplication are separated by a 182-bp spacer. Within each 23S-5S unit, an identical 22-bp spacer separates the 23S and 5S rRNA sequences from each other. The genome organization of the 23S-5S gene cluster in a number of different B. burgdorferi isolates obtained at a number of different geographical locations, as well as in several other species of Borrelia, was investigated. All isolates of B. burgdorferi tested displayed the tandem duplication, whereas the closely related species B. hermsii, B. anserina, and B. turicatae all contained a single copy of each of the genes. In addition, different geographical isolates of B. burgdorferi can be differentiated on the basis of a restriction fragment length polymorphism associated with the 23S-5S gene cluster. This polymorphism can be a useful tool for the determination of genetic relatedness between different isolates of B. burgdorferi.     

5S rDNA and U2 snRNA are linked in the genome of Crassostrea angulata and Crassostrea gigas oysters: does the (CT)n.(GA)n microsatellite stabilize this novel linkage of large tandem arrays?

The 5S rRNA genes from 2 species of the Ostreidae family, Crassostrea angulata and Crassostrea gigas, were molecularly characterized. The genes were amplified, cloned, and sequenced. The results revealed a 5S rDNA tandem array with a nucleotide sequence in an inverted position within the nontranscribed spacer region that corresponded to the U2 small nuclear RNA (snRNA) gene. The sequence analysis indicated that both genes could be functionally active. The presence of the microsatellite (CT)n x (GA)n at the 3' end of both genes and the possible involvement of concerted evolution are discussed.     

Complete rRNA gene sequences reveal that the microsporidium Nosema bombi infects diverse bumblebee (Bombus spp.) hosts and contains multiple polymorphic sites

Characterisation of microsporidian species and differentiation among genetic variants of the same species has typically relied on ribosomal RNA (rRNA) gene sequences. We characterised the entire rRNA gene of a microsporidium from 11 isolates representing eight different European bumblebee (Bombus) species. We demonstrate that the microsporidium Nosema bombi infected all hosts that originated from a wide geographic area. A total of 16 variable sites (all single nucleotid polymorphisms (SNPs)) was detected in the small subunit (SSU) rRNA gene and 42 (39 SNPs and 3 indels) in the large subunit (LSU) rRNA sequence. Direct sequencing of PCR-amplified DNA products of the internal transcribed spacer (ITS) region revealed identical sequences in all isolates. In contrast, ITS fragment length determined by PAGE and sequencing of cloned amplicons gave better resolution of sequences and revealed multiple SNPs across isolates and two fragment sizes in each isolate (six short and seven long amplicon variants). Genetic variants were not unique to individual host species. Moreover, two or more sequence variants were obtained from individual bumblebee hosts, suggesting the existence of multiple, variable copies of rRNA in the same microsporidium, and contrary to that expected for a class of multi-gene family under concerted evolution theory. Our data on within-genome rRNA variability call into question the usefulness of rRNA sequences to characterise intraspecific genetic variants in the Microsporidia and other groups of unicellular organisms.     

Organization of ribosomal RNA gene repeats of the mouse.

The organization of the ribosomal RNA (rRNA) genes of the mouse was determined by Southern blot hybridization using cloned rDNA fragments as probes, which could encompass the entire spacer region between two rRNA gene regions. The rRNA genes are organized into tandem repeats of nearly uniform length of about 44 kb. The heterogeneity detected in the nontranscribed spacer appears to be caused by its sequence rather than its length difference. At least three kinds of repetitive sequences are present in the non-transcribed spacer region; two of them are located 13 kb upstream from the 5'-end of 18S RNA gene and the other located 1 to 4 kb downstream from the 3'-end of 28S RNA gene.     

Contrasting nifD and ribosomal gene relationships among Mesorhizobium from Lotus oroboides in northern Mexico

PCR screens for length variation in a 5' portion of 23S ribosomal RNA and in the 3' end of the 16S rRNA-23S rRNA internal transcribed spacer (ITS) region indicated that nodule bacteria from a Mexican population of Lotus oroboides were diverse on a local scale. Three 23S rRNA length variants and five ITS length variants were detected among the 22 isolates. Sequencing of nearly full-length 16S rRNA genes in three isolates indicated that they fell into the genus Mesorhizobium, but comprised two distinct groups. Two isolates were closely related to M. loti LMG 6125T, while the other isolate clustered with an assemblage of Mesorhizobium taxa that included M. amorphae, M. plurifarium and M. huakuii. However, a phylogenetic tree based on 715 bp of the nitrogenase alpha-subunit (nifD) gene was significantly discordant with the relationships inferred from rRNA sequences. Two isolates that were nearly identical for 16S rRNA had nifD genes that varied at 2% of sites, and one of these nifD sequences was identical to that of another isolate with a strongly divergent 16S rRNA gene. A plasmid screen followed by Southern hybridization indicated that only one of these strains harbored a plasmid-borne nifD gene. These results imply that gene transfer events have altered the distribution of nifD sequences among lineages within this natural population of Mesorhizobium strains.     

Sequence arrangement of the rDNA of Drosophila melanogaster.

The sequence arrangement of genes coding for stable rRNA species and of the interspersed spacers on long single strands of rDNA purified from total chromosomal DNA of Drosophila melanogaster has been determined by a study of the structure of rRNA:DNA hybrids which were mounted for electron microscope observation by the gene 32-ethidium bromide technique. One repeat unit contains the following sequences in the order given. First, an 18 S gene of length 2.13 +/- 0.17 kb. Second, an internal transcribed spacer (Spl) of length 1.58 +/- 0.15 kb. A short sequence coding for the 5.8S and perhaps the 2S rRNA species is located within this spacer. Third, the 28S gene with a length of 4.36 +/- 0.23 kb. About 55% of the 28S genes are unbroken or continuous (C genes). However, about 45% of the 28S genes contain an insertion of an additional segment of DNA that is not complementary to rRNA (l genes). The insertion occurs at a reproducible point 2.99 +/- 0.26 kb from the junction with Spl. The insertions are heterogeneous in length and occur in three broad size classes: 1.42 +/- 0.47, 3.97 +/- 0.55, and 6.59 +/- 0.62 kb. Fourth, an external spacer between the 28S gene and the next 18S gene which is presumably mainly nontranscribed and which has a heterogeneous length distribution with a mean length and standard deviation of 5.67 +/- 1.92 kb. Short inverted repeat stems (100-400 nucleotide pairs) occur at the base of the insertion. It is known from other studies that I genes occur only on the X chromosome. The present study shows that the I and C genes on the X chromosomes are approximately randomly assorted. The sequence arrangement on the plasmid pDm103 containing one repeat of rDNA (Glover et al., 1975) has been determined by similar methods. The I gene on this plasmid contains an inverted repeat stem. The occurrence of inverted repeat sequences flanking the insertion supports the speculation that these sequences are translocatable elements similar to procaryotic translocons.     

Sequence arrangement of the rRNA genes of the dipteran Sarcophaga bullata

Velocity sedimentation studies of RNA of Sarcophaga bullata show that the major rRNA species have sedimentation values of 26S and 18S. Analysis of the rRNA under denaturing conditions indicates that there is a hidden break centrally located in the 26S rRNA species. Saturation hybridization studies using total genomic DNA and rRNA show that 0.08% of the nuclear DNA is occupied by rRNA coding sequences and that the average repetition frequency of these coding sequences is approximately 144. The arrangement of the rRNA genes and their spacer sequences on long strands of purified rDNA was determined by the examination of the structure of rRNa:DNA hybrids in the electron microscope. Long DNA strands contain several gene sets (18S + 26S) with one repeat unit containing the following sequences in order given: (a) An 18S gene of length 2.12 kb, (b) an internal transcribed spacer of length 2.01 kb, which contains a short sequence that may code for a 5.8S rRNA, (c) A 26S gene of length 4.06 kb which, in 20% of the cases, contains an intron with an average length of 5.62 kb, and (d) an external spacer of average length of 9.23 kb.     

Use of RFLPs larger than 100 kbp to map the position and internal organization of the nucleolus organizer region on chromosome 2 in Arabidopsis thaliana

In Arabidopsis thaliana ribosomal RNA genes (rRNA genes or rDNA) are grouped in two nucleolus organizer regions (NORs) that together comprise approximately 6% of the genome. The map positions of the NORs relative to other genetic markers are unknown. It was found that the restriction endonuclease HindIII cuts once in some but not all rRNA genes to yield strain-specific RFLPs of 100–700 kb that could be visualized by pulsed-field gel electrophoresis and Southern blotting. The HindIII RFLPs of the A. thaliana strains Columbia and Landsberg segregated among recombinant inbred lines derived from a cross between these two strains. Linkage analysis placed the NOR bearing the polymorphic HindIII sites to the top of the upper arm of chromosome 2. The name NOR2 is suggested for this locus. HindIII-bearing rRNA genes are interspersed with clusters of HindIII-less genes throughout NOR2. The observed clustering is most consistent with unequal crossing-over, or nearest-neighbor gene conversion, as the mechanism(s) that spread rRNA gene variants throughout an NOR. No meiotic cross-over events yielding a ‘hybrid’ NOR with multiple RFLPs from both parents were observed among the 47 recombinant inbred lines examined. However, the appearance of novel HindIII profiles in approximately 40% of the recombinant inbred lines demonstrates that fluctuations in the distribution of rRNA gene variants occur frequently and can be readily detected on pulsed-field gels.     

The 5S ribosomal RNA gene in Pythium species: two different genomic locations.

The 5S ribosomal RNA (rRNA) genes in eukaryotes may occur either interspersed with the other rRNA genes in the ribosomal DNA (rDNA) repeat, or in separate tandem arrays, or dispersed throughout the genome. In Pythium species and in several related Oomycetes, polymerase chain reaction (PCR) amplification of the nontranscribed spacer (NTS) region with one primer specific for the 5S gene revealed, with several exceptions, that the 5S rRNA gene was present in the rDNA repeat of those species with filamentous sporangia and was absent from the rDNA repeat of those with globose or unknown sporangia. When present, the gene was located approximately 1 kb downstream of the large-subunit rRNA gene and on the strand opposite that on which the other rRNA genes were located. This was confirmed in P. torulosum by sequencing of the gene and its flanking regions. The 5S rRNA genes of P. ultimum were found in tandem arrays outside the rDNA repeat. Evidence of dispersed 5S rRNA sequences was also observed. For many of the species of Pythium having the 5S gene in the NTS, the region between the large-subunit rRNA gene and the 5S gene was found to have length heterogeneity. Oomycetes related to Pythium were also found to have the 5S gene in the NTS, although sometimes in the opposite orientation. This may mean that the presence of the gene in the NTS is ancestral for the Oomycetes and that the absence of the gene in the NTS in those Pythiums with globose sporangia is due to loss of the gene from the rDNA repeat in an ancestor of this group of species. These results favor the view that 5S rRNA gene linkage to the rRNA cistron existed prior to the unlinked arrangement seen in most plants and animals.     

Evolution of dinoflagellate unigenic minicircles and the partially concerted divergence of their putative replicon origins

Dinoflagellate chloroplast genes are unique in that each gene is on a separate minicircular chromosome. To understand the origin and evolution of this exceptional genomic organization we completely sequenced chloroplast psbA and 23S rRNA gene minicircles from four dinoflagellates: three closely related Heterocapsa species (H. pygmaea, H. rotundata, and H. niei) and the very distantly related Amphidinium carterae. We also completely sequenced a Protoceratium reticulatum minicircle with a 23S rRNA gene of novel structure. Comparison of these minicircles with those previously sequenced from H. triquetra and A. operculatum shows that in addition to the single gene all have noncoding regions of approximately a kilobase, which are likely to include a replication origin, promoter, and perhaps segregation sequences. The noncoding regions always have a high potential for folding into hairpins and loops. In all six dinoflagellate strains for which multiple minicircles are fully sequenced, parts of the noncoding regions, designated cores, are almost identical between the psbA and 23S rRNA minicircles, but the remainder is very different. There are two, three, or four cores per circle, sometimes highly related in sequence, but no sequence identity is detectable between cores of different species, even within one genus. This contrast between very high core conservation within a species, but none among species, indicates that cores are diverging relatively rapidly in a concerted manner. This is the first well-established case of concerted evolution of noncoding regions on numerous separate chromosomes. It differs from concerted evolution among tandemly repeated spacers between rRNA genes, and that of inverted repeats in plant chloroplast genomes, in involving only the noncoding DNA cores. We present two models for the origin of chloroplast gene minicircles in dinoflagellates from a typical ancestral multigenic chloroplast genome. Both involve substantial genomic reduction and gene transfer to the nucleus. One assumes differential gene deletion within a multicopy population of the resulting oligogenic circles. The other postulates active transposition of putative replicon origins and formation of minicircles by homologous recombination between them.     

Chromosomal dynamics of nucleolar organizer regions (NORs) in the house mouse: micro-evolutionary insights

Variation in the number and chromosomal location of nucleolar organizer regions (NORs) was studied in the house mouse, Mus musculus (2n=40). From an origin in Western Asia, this species colonized the Middle East, Europe and Asia. This expansion was accompanied by diversification into five subspecies. NOR diversity was revealed by fluorescence in situ hybridization using 18S and 28S probes on specimens spanning Asia to Western Europe. The results showed that the house mouse genome possessed a large number of NOR-bearing autosomes and a surprisingly high rate of polymorphism for the presence/absence of rRNA genes on all these chromosomes. All NOR sites were adjacent to the centromere except for two that were telomeric. Subspecific differentiation established from the NOR frequency data was concordant with the overall pattern of radiation proposed from molecular studies, but highlighted several discrepancies that need to be further addressed. NOR diversity in M. musculus consisted of a large number of polymorphic NORs that were common to at least two subspecies, and a smaller number of NORs that were unique to one subspecies. The most parsimonious scenario argues in favor of a subspecific differentiation by lineage sorting of ancestral NOR polymorphisms; only the unique NORs would have appeared by inter-chromosomal transposition, except for the two telomeric ones that may have originated by hybridization with another species. Such a scenario provides an alternative view from the one prevailing in most systematic and phylogenetic analyses that NORs have a high transposition rate due to concerted evolution of rRNA genes.