|
|
||||||||
Systematics and Evolution |
Microbial Collection and Screening Laboratory, Institute of Molecular and Cell Biology, 30 Medical Drive, Singapore 1176091
Author for correspondence: Yue Wang. Tel: +65 7783207. Fax: +65 7791117. e-mail: mcbwangy{at}imcb.nus.edu.sg
| ABSTRACT |
|---|
|
|
|---|
Keywords: complexity hypothesis, lateral gene transfer, mosaic gene, phylogeny, evolution
Abbreviations: LGT; lateral gene transfer
The GenBank accession numbers for the 23S rRNA sequences determined in this study are AF192136AF192150.
| INTRODUCTION |
|---|
|
|
|---|
Although the complexity hypothesis is largely correct, it overlooks one important aspect of LGT: the gene is not the smallest unit of transfer. Sawyer (1989
) deduced that the average units of LGT are DNA patches of 100200 bp, suggesting that transfer of parts of genes is common in nature. When a gene is broken down to small pieces, the complexity will be greatly simplified. For example, a gene fragment encoding a domain, such as a hairpin of an rRNA molecule, may replace a corresponding region of the host gene via a single recombinational event. A new hairpin in an rRNA molecule may not be detrimental to the organism as long as the local secondary structure remains identical or similar (Gutell et al., 1994
; Van de Peer et al., 1996
; Asai et al., 1999
). This process may occur repeatedly, involving different parts of a gene, leading to gradual changes of the nucleotide sequence. This simplified complexity hypothesis predicts a mosaic nature for rRNA genes (as well as for other genes). Unlike the replacement of an entire gene, the step-by-step segmental replacement gradually obscures the evolutionary history documented in the nucleotide sequence. Although there have been previous reports of mosaic genes arising from local gene exchange between related bacterial species (Smith et al., 1991
; Groisman et al., 1992
), it is not clear whether and how frequently rRNA genes are affected by the lateral transfer of gene segments.
In this study, we investigated the distribution of nucleotide variations in the alignments of the 16S rDNA sequences of actinomycete species belonging to several well-defined genera. We frequently found the occurrence of two or three types of drastically different sequences in short and rather conserved regions of rRNA genes among species of the same genus. The nucleotide substitutions involve three to five base pairs of compensating covariations and are found within the stems of individual hairpins in the rRNA molecule. These observations strongly support our simplified complexity hypothesis to explain how LGT may affect the highly conserved genes. We also demonstrate how this segmental gene transfer may lead to confusion in rDNA sequence-based phylogenetic analysis.
| METHODS |
|---|
|
|
|---|
Preparation of genomic DNAs.
The genomic DNA of the actinomycetes was prepared as previously described (Wang et al., 1996
).
PCR amplification, cloning and sequencing.
The 5' one-third of the 23S rRNA genes was amplified by PCR using a pair of primers, one targeting a conserved region at the end of the 16S rRNA gene and the other a conserved block within the 23S rRNA gene. The PCR product also includes the spacer region. The sequences of the two oligonucleotides are as follows: 5'-GGTTGGATCCACCTCCTT-3', corresponding to nt 15251542 of E. coli 16S rRNA (Brosius et al., 1978
); and 5'-ACCAGTGAGCTATTAGCG-3' (nt 10901107). After the cloning of the PCR products, the M13 forward and reverse universal primers were used for sequencing the ends of each clone. The internal regions were sequenced in both orientations by using the following two sets of oligonucleotide primers targeting two conserved sequences within 23S rDNA. The first set of primers, targeting nt 4560 of E. coli 23S rRNA gene, are 23S-40f (5'-CCGATGAAGGACGTGGGA-3') and 23S-40r (5'-TCCCACGTCCTTCATCGG-3'); and the second set of primers, targeting nt 456472, are 23S-460f (5'-CCTTTCCCTCACGGTACT-3') and 23S-460r (5'-AGTACCGTGAGGGAAAGG-3').
Sequence alignment and phylogenetic analysis.
Multiple sequence alignment of the rRNA gene sequences was carried out by the CLUSTAL method of the DNASTAR program (DNASTAR, Inc.) and verified according to the consensus secondary structure model (Gutell et al., 1994
). Phylogenetic trees were constructed by using both a distance method (Saitou & Nei, 1987
) and a parsimony method (Swofford & Begle, 1993
). The two methods generated very similar trees and only the neighbour-joining trees are shown in this paper. The confidence level of phylogenetic tree topology was evaluated by the bootstrap method (Felsenstein, 1985
). The softwares for the tree construction and bootstrap analysis are contained in the CLUSTAL V (Higgins et al., 1992
) and the PAUP (Swofford & Begle, 1993
) phylogenetic analysis software packages.
| RESULTS |
|---|
|
|
|---|
|
1200 bases, nt 11108 of E. coli numbering) sequence analysis conducted in this study largely confirms the phylogenetic integrity of each genus, but demonstrates significant discrepancy at the subgenus level. For example, Fig. 2
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
We initiated this study by hypothesizing that LGT might affect complex genes, such as rRNA genes, through gradual substitution of gene segments encoding individual domains. Though the overall structure and its interaction with other parts of the translation machinery are complex, an rRNA molecule consists of many hairpin structures, each of which may interact with only one or a few other components of the translation apparatus. In addition, it is well known that the secondary but not the primary structure of many hairpins in rRNA is essential for function (Gutell et al., 1994
; Van de Peer et al., 1996
). These two properties of rRNA molecules underscore the high possibility of the exchange of individual hairpins without damaging effect on the proper function of a ribosome. This view is supported by reports that an organism may have distinct types of rRNA genes with high levels of base variations throughout the entire gene (Gunderson et al., 1987
; Mylvaganam & Dennis, 1992
; Carranza et al., 1996
; Wang et al., 1997
; Yap et al., 1999
). Niebel et al. (1987
) demonstrated that the rRNA cistrons from Proteus vulgaris were expressed and the products correctly processed and assembled into ribosomes when transformed into E. coli. Recently, Asai et al. (1999
) reported a complete exchange of rRNA genes between different bacterial species and replacement of a 23S rRNA gene segment of E. coli by the corresponding region of yeast. If a transgressed rRNA gene or gene segment brings to the recipient a beneficial property, such as antibiotic resistance (Green et al., 1997
; Mankin, 1997
), it will be kept and may eventually replace the endogenous gene(s) via gene conversion (Hillis et al., 1990
). High levels of localized intergenic base variations corresponding to the stem of a hairpin have also been reported (Ueda et al., 1999
; Yap & Wang, 1999
). Ueda et al. (1999
) found two types of base variations in the highly variable
region of the 16S rRNA gene and suggested that both random mutation and LGT contribute to rRNA gene heterogeneity of a bacterium.
In this study, we found high levels of base variations localized to the stem of hairpins of rRNA molecules that seldom experience random mutations between species of the same genus or closely related genera. We propose that the LGT of gene segments is the most plausible explanation of this phenomenon. We observed this phenomenon in the 16S rRNA genes of almost all genera we examined, indicating that such events probably happen with considerable frequency in nature. Although it is difficult, if not impossible, to determine the donor of a transferred DNA fragment, identical or nearly identical sequences are always found in species often belonging to distantly related taxa. Introduction of a piece of foreign DNA will certainly corrupt the evolutionary history recorded in the affected gene, though the extent of damage will depend on the size and the number of changed bases of the region involved. In most cases we have investigated, the presence of sequence variations in one hairpin does not affect the stable aggregation in one clade of all the species of a genus, reflecting the overall evolutionary stability of 16S rRNA genes. At subgenus level, the pattern of species clustering may be affected by the presence of these localized sequence variations, as shown in the analysis of the 16S rDNA of Nonomuraea species. However, it is not clear whether the sharing of a short DNA sequence is a result of common ancestry, or whether the sequences were introduced in independent events. The cluster pattern of Nonomuraea species shown in the 16S rDNA trees, constructed by including or excluding the short regions of sequence variations, was not reproducible in the 23S rRNA tree. In the study of Streptosporangium species, we speculate that the lateral transfer of a larger DNA fragment encoding two adjacent hairpins was the likely cause of the confusion at genus level of the positions of S. corrugatum and S. claviforme in the 16S rRNA tree. Regardless of whether the 16S or the 23S tree is correct, the dramatically different positions of the two species in the two trees are unlikely to be the result of random mutational events.
Conclusion
Complementing the complexity hypothesis, the simplified complexity hypothesis suggests that due to the overall functional constraints, the lateral transfer of parts of genes is likely to be the main mode of LGT occurring to the genes encoding components of complex systems. This mode of LGT results in gradual corruption of the evolutionary history written in the nucleotide sequences. The rate of sequence corruption of a gene should be a function of the complexity of interactions with other cellular components by the gene product. Even highly conserved genes like rRNA genes may only roughly record part of the evolutionary history of life, though how far back is not known. There is probably an increasingly high restriction for successful transfer of gene segments with widening evolutionary distance, which may explain why the three-domain separation of life largely holds when proteins of the information-processing machinery are used to derive phylogeny (Woese, 1998
; Jain et al., 1999
). The results of this study caution that in bacterial taxonomy the effect of LGT on the result of sequence-based phylogenetic analysis should be considered. Analysis of more than one gene set together with the examination of major phenotypic characteristics should be practised in future. Given that LGT is well established, it is now important to identify the principles that govern gene distribution patterns across prokaryotic genomes to determine the values of genes of different functional groups in bacterial taxonomy.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Brosius, J., Palmer, J. J., Kennedy, J. P. & Noller, H. F. (1978). Complete nucleotide sequence of a 16S ribosomal gene from Escherichia coli. Proc Natl Acad Sci USA 75, 4801-4805.
Carranza, S., Giribet, G., Riberat, C., Baguna, J. & Riutort, M. (1996). Evidence that two types of 18S rDNA coexist in the genome of Dugesia (Schmidtea) mediterranea (Platyhelminthes, Turbellaria, Tricladida). Mol Biol Evol 13, 824-832.[Abstract]
Chun, J. & Goodfellow, M. (1995). A phylogenetic analysis of the genus Nocardia with 16S rRNA gene sequences. Int J Syst Bacteriol 45, 240-245.
Doolittle, W. F. (1999). Phylogenetic classification and universal tree. Science 284, 2124-2128.
Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 46, 159-173.
Green, R., Samaha, R. R. & Noller, H. F. (1997). Mutations at nucleotides G2251 and U2585 of 23S rRNA perturb the peptidyl transferase center of the ribosome. J Mol Biol 266, 40-50.[Medline]
Groisman, E. A., Saier, M. H.Jr & Ochman, H. (1992). Horizontal transfer of a phosphatase gene as evidence for mosaic structure of the Salmonella genome. EMBO J 11, 1309-1316.[Medline]
Gunderson, J. H., Sogin, M. L., Wollet, G., Hollingdale, M., de la Cruz, V. F., Waters, A. P. & McCutchan, T. F. (1987). Structurally distinct, stage-specific ribosomes occur in Plasmodium. Science 238, 933-937.
Gutell, R. G., Larsen, N. & Woese, C. R. (1994). Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective. Microbiol Rev 58, 10-26.
Higgins, D. G., Bleasby, A. J. & Fuchs, R. (1992). CLUSTAL V: improved software for multiple sequence alignment. Comput Appl Biosci 8, 189-191.
Hillis, D. M., Moritz, C., Porter, C. A. & Baker, R. J. (1990). Evidence for biased gene conversion in concerted evolution of ribosomal DNA. Science 251, 308-310.
Jain, R., Rivera, M. C. & Lake, J. A. (1999). Horizontal gene transfer among genomes: the complexity hypothesis. Proc Natl Acad Sci USA 96, 3801-3806.
Koch, C., Kroppenstedt, R. M., Rainey, F. A. & Stackebrandt, E. (1996). 16S ribosomal DNA analysis of the genera Micromonospora, Actinoplanes, Catellatospora, Catenuloplanes, Couchioplanes, Dactylosporangium, and Pilimelia and emendation of the family Micromonosporaceae. Int J Syst Bacteriol 46, 765-768.
Koonin, E. V., Mushegian, A. R., Galperin, M. Y. & Walker, D. R. (1999). Comparison of archaeal and bacterial genomes: computer analysis of protein sequences suggests a chimeric origin for the archaea. Mol Microbiol 25, 619-637.
Lawrence, J. G. & Ochman, H. (1997). Amelioration of bacterial genome: rates of change and exchange. J Mol Evol 44, 383-397.[Medline]
Mankin, A. S. (1997). Pactamycin resistance mutations in functional sites of 16S rRNA. J Mol Biol 274, 8-15.[Medline]
Mylvaganam, S. & Dennis, P. P. (1992). Sequence heterogeneity between the two genes encoding 16S rRNA from the halophilic archaebacterium Halobacterium marismortui. Genetics 130, 399-410.[Abstract]
Nelson, K. E., Clayton, R. A., Gill, S. R. & 22 other authors (1999). Evidence for lateral gene transfer between Archaea and Bacteria from genome sequence of Thermotoga maritima. Nature 399, 323329.[Medline]
Niebel, H., Dorsch, M. & Stackebrandt, E. (1987). Cloning and expression in Escherichia coli of Proteus vulgaris genes for 16S ribosomal RNA. J Gen Microbiol 133, 2401-2409.[Medline]
Rheims, H., Schumann, P., Rohde, M. & Stackebrandt, E. (1998). Verrucosispora gifhornensis gen. nov., sp. nov., a new member of the actinobacterial family Micromonosporaceae. Int J Syst Bacteriol 48, 1119-1127.
Rivera, M. C., Jain, R., Moore, J. E. & Lake, J. A. (1998). Genomic evidence for two functionally distinct gene classes. Proc Natl Acad Sci USA 95, 6239-6244.
Saitou, N. & Nei, M. (1987). The neighbour-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4, 406-425.[Abstract]
Sawyer, S. (1989). Statistical tests for detecting gene conversion. Mol Biol Evol 6, 526-538.[Abstract]
Smith, J. M., Dowson, C. G. & Spratt, B. G. (1991). Localized sex in bacteria. Nature 349, 29-31.[Medline]
Swofford, D. L. & Begle, D. P. (1993). Phylogenetic Analysis Using Parsimony (Version 3.1), Users Manual. Champaign, IL: Smithsonian Institute Laboratory of Molecular Systematics.
Syvanen, M. (1994). Horizontal gene transfer: evidence and possible consequences. Annu Rev Genet 28, 237-261.[Medline]
Ueda, K., Seki, T., Kudo, T., Yoshida, T. & Kataoka, M. (1999). Two distinct mechanisms cause heterogeneity of 16S rRNA. J Bacteriol 181, 78-82.
Van de Peer, Y., Chapelle, S. & Wachter, R. D. (1996). A quantitative map of nucleotide substitution rates in bacterial rRNA. Nucleic Acids Res 24, 3381-3391.
Wang, Y., Zhang, Z. S. & Ruan, J. S. (1996). A proposal to transfer Microbispora bispora (Lechevalier 1965) to a new genus, Thermobispora gen. nov., as Thermobispora bispora comb. nov. Int J Syst Bacteriol 46, 933-938.
Wang, Y., Zhang, Z. S. & Ramanan, N. (1997). The actinomycete Thermobispora bispora contains two distinct types of transcriptionally active 16S rRNA genes. J Bacteriol 179, 3270-3276.
Ward-Rainey, N., Rainey, F. A. & Stackebrandt, E. (1996). The phylogenetic structure of the genus Streptosporangium. Syst Appl Microbiol 19, 50-55.
Woese, C. R. (1987). Bacterial evolution. Microbial Rev 51, 221-271.
Woese, C. R. (1998). The universal ancestor. Proc Natl Acad Sci USA 95, 6854-6859.
Woese, C. R., Kandler, O. & Wheelis, M. L. (1990). Towards a natural system of organisms: proposal for the domains Archaea, Bacteria and Eucarya. Proc Natl Acad Sci USA 87, 4576-4579.
Yap, W. H. & Wang, Y. (1999). Molecular cloning and comparative sequence analyses of rRNA operons in Streptomyces nodosus ATCC 14899. Gene 232, 77-85.[Medline]
Yap, W. H., Zhang, Z. S. & Wang, Y. (1999). Distinct types of rRNA operons exist in the genome of the actinomycete Thermomonospora chromogena and evidence for horizontal transfer of an entire rRNA operon. J Bacteriol 181, 5201-5209.
Zhang, Z. S., Wang, Y. & Ruan, J. S. (1998). Reclassification of Thermomonospora and Microtetraspora. Int J Syst Bacteriol 48, 411-422.
Received 8 February 2000;
revised 22 July 2000;
accepted 7 August 2000.
This article has been cited by other articles:
![]() |
D. Ghosh, K. Roy, K. E. Williamson, D. C. White, K. E. Wommack, K. L. Sublette, and M. Radosevich Prevalence of Lysogeny among Soil Bacteria and Presence of 16S rRNA and trzN Genes in Viral-Community DNA Appl. Envir. Microbiol., January 15, 2008; 74(2): 495 - 502. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Ueno, V. A. R. Huss, N. Urano, and S. Watabe Direct evidence for redundant segmental replacement between multiple 18S rRNA genes in a single Prototheca strain Microbiology, November 1, 2007; 153(11): 3879 - 3893. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. P. Andam and M. A. Parker Novel Alphaproteobacterial Root Nodule Symbiont Associated with Lupinus texensis Appl. Envir. Microbiol., September 1, 2007; 73(17): 5687 - 5691. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Mincer, W. Fenical, and P. R. Jensen Culture-Dependent and Culture-Independent Diversity within the Obligate Marine Actinomycete Genus Salinispora Appl. Envir. Microbiol., November 1, 2005; 71(11): 7019 - 7028. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. E. Dewhirst, Z. Shen, M. S. Scimeca, L. N. Stokes, T. Boumenna, T. Chen, B. J. Paster, and J. G. Fox Discordant 16S and 23S rRNA Gene Phylogenies for the Genus Helicobacter: Implications for Phylogenetic Inference and Systematics J. Bacteriol., September 1, 2005; 187(17): 6106 - 6118. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. T. Sacchi, D. Alber, P. Dull, E. A. Mothershed, A. M. Whitney, G. A. Barnett, T. Popovic, and L. W. Mayer High Level of Sequence Diversity in the 16S rRNA Genes of Haemophilus influenzae Isolates Is Useful for Molecular Subtyping J. Clin. Microbiol., August 1, 2005; 43(8): 3734 - 3742. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. Trujillo, A. Willems, A. Abril, A.-M. Planchuelo, R. Rivas, D. Ludena, P. F. Mateos, E. Martinez-Molina, and E. Velazquez Nodulation of Lupinus albus by Strains of Ochrobactrum lupini sp. nov. Appl. Envir. Microbiol., March 1, 2005; 71(3): 1318 - 1327. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Miller, S. Augustine, T. L. Olson, R. E. Blankenship, J. Selker, and A. M. Wood Discovery of a free-living chlorophyll d-producing cyanobacterium with a hybrid proteobacterial/cyanobacterial small-subunit rRNA gene PNAS, January 18, 2005; 102(3): 850 - 855. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. Schouls, C. S. Schot, and J. A. Jacobs Horizontal Transfer of Segments of the 16S rRNA Genes between Species of the Streptococcus anginosus Group J. Bacteriol., December 15, 2003; 185(24): 7241 - 7246. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. D. Crosbie, M. Pockl, and T. Weisse Dispersal and Phylogenetic Diversity of Nonmarine Picocyanobacteria, Inferred from 16S rRNA Gene and cpcBA-Intergenic Spacer Sequence Analyses Appl. Envir. Microbiol., September 1, 2003; 69(9): 5716 - 5721. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. van Berkum, Z. Terefework, L. Paulin, S. Suomalainen, K. Lindstrom, and B. D. Eardly Discordant Phylogenies within the rrn Loci of Rhizobia J. Bacteriol., May 15, 2003; 185(10): 2988 - 2998. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. P. Gogarten, W. F. Doolittle, and J. G. Lawrence Prokaryotic Evolution in Light of Gene Transfer Mol. Biol. Evol., December 1, 2002; 19(12): 2226 - 2238. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Parker, B. Lafay, J. J. Burdon, and P. van Berkum Conflicting phylogeographic patterns in rRNA and nifD indicate regionally restricted gene transfer in Bradyrhizobium Microbiology, August 1, 2002; 148(8): 2557 - 2565. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Parker Bradyrhizobia from Wild Phaseolus, Desmodium, and Macroptilium Species in Northern Mexico Appl. Envir. Microbiol., April 1, 2002; 68(4): 2044 - 2048. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Moreno, J. Romero, and R. T. Espejo Polymorphism in repeated 16S rRNA genes is a common property of type strains and environmental isolates of the genus Vibrio Microbiology, April 1, 2002; 148(4): 1233 - 1239. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. W. Friedrich Phylogenetic Analysis Reveals Multiple Lateral Transfers of Adenosine-5'-Phosphosulfate Reductase Genes among Sulfate-Reducing Microorganisms J. Bacteriol., January 1, 2002; 184(1): 278 - 289. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Keeling and J. D. Palmer Lateral transfer at the gene and subgenic levels in the evolution of eukaryotic enolase PNAS, August 23, 2001; (2001) 191337098. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Parker Case of Localized Recombination in 23S rRNA Genes from Divergent Bradyrhizobium Lineages Associated with Neotropical Legumes Appl. Envir. Microbiol., May 1, 2001; 67(5): 2076 - 2082. [Abstract] [Full Text] |
||||
![]() |
P. J. Keeling and J. D. Palmer Lateral transfer at the gene and subgenic levels in the evolution of eukaryotic enolase PNAS, September 11, 2001; 98(19): 10745 - 10750. [Abstract] [Full Text] [P |