Previous Article | Next Article ![]()
Journal of Bacteriology, February 2002, p. 666-671, Vol. 184, No. 3
0021-9193/01/$04.00+0 DOI: 10.1128/JB.184.3.666-671.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
MATFORSK, Norwegian Food Research Institute, 1430 Ås,1 Division of General Genetics, Department of Biology, University of Oslo, Blindern, 0315 Oslo, Norway2
Received 15 June 2001/ Accepted 29 October 2001
|
|
|---|
- and ß-proteobacteria, respectively, than with other cyanobacterial tRNALeu(UAA) group I introns. The divergent tRNALeu(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNALeu(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNALeu(UAA) gene which is located within a unstable genome region. |
|
|---|
A group I intron-encoded endonuclease specific for the anticodon region of tRNAfMet has been identified in Synechocystis sp. strain PCC 6803. The distribution pattern of this endonuclease suggests that it may be mobile (3). The presence of mobile endonucleases, or endonucleases located outside the intron, has previously been suggested based on the distribution pattern of the group I introns in cyanobacteria (19). The recent discovery of an intron-encoded endonuclease gives further support for a complex pattern of evolution of tRNA group I introns in bacteria.
The work presented here is a further step toward defining possible mechanisms for the evolution of tRNA introns in bacteria (5). We have addressed the evolution and origin of the divergent tRNALeu(UAA) intron within the Microcystis genus. The origin of the intron was investigated by comparative phylogenetic reconstructions of complete intron and the corresponding tRNALeu(UAA) exon sequences. These data support the intron being horizontally transferred between different tRNA genes. The next question we addressed was whether the sporadic intron distribution within the Microcystis genus (21) can be explained by intron mobility or, rather, by mechanisms involving instability of the whole intron-containing genomic region. Our results obtained by PCR amplification and Southern hybrization show that the whole region containing both the intron and the exon is absent in the intronless strains. Thus, the distribution of tRNALeu(UAA) group I introns within the Microcystis genus can be attributed to both intron transfer and recombination of the entire tRNA and intron-containing region.
|
|
|---|
50-mg (wet weight) cell pellets by a standard phenol-chloroform method (22). PCR amplification, cloning, and sequencing. The PCR primers used in this work are described in Table 1. Between 30 and 40 cycles were used in the amplification reactions, with a denaturation at 95°C for 30 s, annealing for 30 s at the temperatures described in Table 1, and extension at 72°C for 30 s to 1 min. All reactions were initiated by a 4-min denaturation at 94°C and terminated with a 7-min extension at 72°C. The amplified products were cloned by using the TA-cloning pGEM-T Vector system (Promega, Madison, Wis.). The sequencing was done with the Model 373 DNA Sequencer (Applied Biosystems, Foster City, Calif.). The TA-cloning and subsequent sequencing were done according to standard protocols (Promega and Applied Biosystems, respectively).
|
View this table: [in a new window] |
TABLE 1. Primer sequences with corresponding Tm values
|
Single-stranded 32P-labeled probes were generated from PCR-amplified DNA by using a random primer DNA labeling kit (Boehringer Mannheim, GmbH, Mannheim, Germany) as described in Espelund et al. (6). The regions used for probes are displayed in Fig. 1. The exon probe for tRNALeu(UAA) was made with a synthetic tRNALeu(UAA) type B gene (19) without intron as a template.
![]() View larger version (11K): [in a new window] |
FIG. 1. Cloned flanking region, included probe, and primer regions for the primers used to generate the probes for M. aeruginosa N-C 57 (EMBL accession no. AJ307004). Two partial ORFs (RF1 and RF2) were identified in the flanking region. The primer sites and orientation are indicated with arrows, while the probe regions are indicated with bars. The type B tRNALeu(UAA) exon probe is denoted with a bar with a stippled line to indicate that stippled region is not present in the probe.
|
Separate phylogenetic trees were constructed by using the maximum-parsimony method (9) and the distance method (25) provided in the software package PAUP*4.0 developed by D. L. Swoford (Florida State University, Tallahassee) and distributed by Sinauer Associates, Inc. (Sunderland, Mass.). In addition, maximum-likelihood methods were provided in the Puzzle 4.0.2 software (27). Search for the maximum-parsimony tree was done heuristically. Both maximum-likelihood distances with empirically determined base frequencies (7) and LogDet distances were used for the minimum evolution distance trees. We used the model developed by Hasegawa et al. (11) for the maximum-likelihood analysis. For the parsimony, and the distance methods consensus trees were constructed from 1,000 bootstrap replicates (8), while for the maximum-likelihood analysis 1,000 quarter puzzling steps were used (27). The tree length skewness, determined by the third momentum statistics (g1) for 10,000 randomly generated trees, was used to investigate the phylogenetic structure of the data. The results obtained were compared to critical values given by Hillis and Huelsenbeck (12).
|
|
|---|
Evolutionary origin of M. aeruginosa tRNALeu(UAA) exon sequences. The intron-containing Microcystis tRNALeu(UAA) gene has the typical cloverleaf secondary structure of tRNA genes (Fig. 2). To determine the evolutionary origin of the intron-containing Microcystis tRNALeu(UAA) gene, a representative set of known tRNA genes were selected from the tRNA database (26). The tRNA sequences were aligned both manually and by using the information obtained from the alignment provided in the tRNA database.
![]() View larger version (19K): [in a new window] |
FIG. 2. Secondary structure of the intron-containing tRNALeu(UAA) gene. The intron insertion site is indicated by an arrow. The structure was obtained by a combination of MFOLD computer prediction (GCG Package; Genetic Computer Group, Madison, Wis.) and manual editing.
|
![]() View larger version (38K): [in a new window] |
FIG. 3. Phylogenetic reconstruction of tRNA exon (A) and intron (B) sequences. The trees were built by using both distance and maximum-parsimony methods. The trees shown in panels A and B (based on 100 and 275 aligned positions, respectively) were built by the neighbor-joining method using maximum-likelihood distances. The distances are expressed as substitutions per nucleotide in the neighbor-joining tree. Numbers at the nodes (expressed as maximum likelihood/LogDet distances/maximum-likelihood distances/heuristic search for maximum parsimony trees) indicate the percentage of Puzzle (for maximum likelihood) and bootstrap trees for the others in which the cluster descending from the node was found. Numbers are only shown for nodes supported by 50% for all of the methods tested. Each tRNA gene is annotated in panel A, while the three clusters of introns identified (panel B) are indicated by I, II, and III.
|
The intron data set (10,000 randomly generated trees) also produced a significant (P = 0.01) g1 of 0.77 according to the criteria given by Hillis and Huelsenbeck (12). The tRNA introns identified within the Microcystis genus are closely related, with a homology of >99% (results not shown). The intron alignment supported a clustering of the Microcystis introns with the introns in tRNAIle(CAU) and tRNAArg(CCU) genes. The other cyanobacterial and chloroplast tRNALeu(UAA) introns, on the other hand, formed a separate cluster (Fig. 3B).
Flanking regions for tRNALeu(UAA) intron-containing gene. The generated sequences were deposited in the database with EMBL accession no. AJ307004. A partial open reading frame (ORF) of 218 amino acids was identified upstream of the intron-containing tRNA gene (Fig. 1). No homology to other amino acid sequences (ORFs) were identified by a TBLASTN search through all six reading frames in the NCBI database (www.ncib.nlh.gov; 1 May 2001). The most closely related sequence was a region on chromosome 1 in Deinococcus radiodurans R, sequence section 209, nucleotides 6238 to 6414 (accession no. AE002072). There was a 33% identity in this region. A 48-amino-acid partial ORF was identified on the complementary strand downstream of the tRNALeu(UAA) gene. No significant matches were found in the database with a TBLASTN search with this partial ORF.
Distribution of tRNALeu(UAA) introns, flanks, and exons. The distribution of the intron, flank, and exon sequences were accessed by both PCR amplification and Southern hybridization analyses. We constructed specific PCR primers (CAB and B-CBB) for selective amplification of the intron-containing tRNALeu(UAA) gene. This primer pair was then used in a PCR screening for isolates containing tRNALeu(UAA) genes of type B. All type B tRNALeu(UAA) genes identified contained introns, including the Nostoc strain N-C 308. The strains without identified introns (19, 21) did not yield any amplification products with this primer pair (including several strains belonging to the Nostoc group), suggesting the absence of type B tRNA genes (Fig. 4). The products were sequence verified, and the sequences obtained were confirmed through comparison to the published intron sequences sequences: accession no. y13474 for M. aeruginosa N-C 57, accession no. AJ228695 for M. aeruginosa N-C 228/1, and accession no. AJ228712 for Nostoc sp. strain N-C 308.
![]() View larger version (63K): [in a new window] |
FIG. 4. PCR amplification analyses with the primer pair CAB-B-CBB (A) and through control amplification of 16S rDNA with the primer pair CC-CD (B). All samples were electrophoresed in 1.5% agarose gels for 30 min at 100 V. Twenty percent of the amplification products was loaded in each lane. MW, molecular weight marker.
|
![]() View larger version (92K): [in a new window] |
FIG. 5. Southern hybridization analysis of selected Microcystis strains. Genomic DNA digested with HindIII was separated, blotted, and hybridized as described in Materials and Methods. The same membrane was used in each of the hybridizing experiments, providing an exact assignment of the hybridizing bands from different experiments. The hybridizing probes were as follows: flank 1 (A), tRNA type A (B), tRNA type B (C), intron (D), and 16S rDNA (E).
|
![]() View larger version (47K): [in a new window] |
FIG. 6. Southern hybridization analysis with AluI-digested DNA for the M. aeruginosa strains N-C 57 and 228/1. Digested genomic DNA was separated, blotted, and hybridized as described in Materials and Methods. The following probes were used: flank 2 (A) and intron (B).
|
|
|
|---|
Phylogenetic reconstruction support intron transfer. Previously, data have been presented that support horizontal intron transfer between tRNALeu(UAA), tRNAIle(CAU) and/or tRNAArg(CCA) genes (21). These results, however, were based on partial tRNA sequences (only 28 aligned positions) which do not yield sufficient information for rigid conclusions. Here we have cloned the complete tRNALeu(UAA) exon sequence from the intron-containing tRNA gene in the M. aeruginosa strain N-C 57 (20). The comparison of complete exon sequences gave further evidence for separate clustering of tRNALeu(UAA) apart from the tRNAIle(CAU) and tRNAArg(CCA) genes (see Fig. 3A).
The sequences for the tRNALeu(UAA) intron in Microcystis, on the other hand, have a clear affiliation with the introns in tRNAIle(CAU) and tRNAArg(CCA) (see Fig. 3B). Theoretically, it cannot be excluded that there has been concerted selection on certain primary sequence elements, resulting in substitutions that do not reflect the evolutionary history. Such an argument, however, questions the whole basis of modern phylogeny (for microorganisms, see the review by Woese [29]). The underlying assumption is that the majority of the genetic changes are due to neutral and not to selective substitutions (13). If we assume that the majority of genetic changes in our data set are neutral, then the most likely explanation for the different phylogenies for tRNALeu(UAA) exons and introns is horizontal intron transfer.
No homology to other known sequences for the identified reading frames. The inferred amino acid sequence of the partial ORFs upstream and downstream of the tRNALeu(UAA) gene did not show any significant similarities with previously characterized amino acid sequences. The potential functions can therefore not be identified from previously characterized genetic elements in the database. Research addressing whether or not the genetic elements in the flanking regions of the intron-containing tRNALeu(UAA) gene are important for the intron distribution is thus called for. The presence of potentially mobile endonucleases without homology to other known sequences has recently been identified for a cyanobacterial group I intron (3).
The distribution of tRNALeu(UAA) introns, exons, and flanking regions in the Microcystis genus suggests instability of the whole region. Both the flanks and the intron-containing tRNA gene are present in the strains N-C 57 and 228/1 but absent in the other investigated strains (see Fig. 5). These results strongly suggest that the observed intron distribution in the Microcystis genus is not due to the mobility of the intron itself but rather to instability of the whole intron-containing region. We do not know whether this distribution pattern is due to insertions or deletions or if the intron-containing tRNA gene is located on an episomal element. However, we conclude that relatively recent event(s) have caused the distribution pattern since the different Microcystis strains are closely related (23, 24).
Model for the evolution of the Microcystis tRNALeu(UAA) group I introns. We have presented the first statistically supported phylogenetic evidence for the horizontal transfer of group I introns been between different tRNA genes in bacteria. Furthermore, we have linked the location of the intron-containing tRNA gene to a region in the genome prone to recombinations. Our data support that the tRNA introns can evolve both by being localized in an unstable genetic element or region and via an intrinsic property of (infrequent) intron mobility (1, 15) (Fig. 7).
![]() View larger version (14K): [in a new window] |
FIG. 7. Model for the evolution of the Microcystis tRNALeu(UAA) group I intron. The model shows that the process of intron homing is relatively infrequent (A) and that the process of recombination of the entire tRNALeu(UAA) genome region is relatively frequent (B). Taken together, these two processes can explain the origin of the Microcystis introns and the current intron distribution within this genus.
|
|
|
|---|
The work was supported by grant 107622/420 from the Norwegian Research Council to K.S.J and in part by a research levy on certain agricultural products.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»