*Departamento de Biología Vegetal II, Facultad de Farmacia, Universidad Complutense de Madrid, Madrid, Spain;
and
School of Biological and Chemical Sciences, Birkbeck, University of London, London, England, and Mycology Section, Royal Botanic Garden, Kew, Richmond, England
Abstract
Twenty-four new insertions were obtained from seven different locations in the nuclear 18S rDNA for seven species of the lichen-forming fungal genus Physconia. They were analyzed allowing for terminal sequence conservation by adopting a flexible approach to exact insertion site position, and they were compared with 12 previously reported small insertion sequences from the 18S ribosomal RNA gene. Such insertions have previously been proposed to be degenerate self-splicing group I introns; however, the methodology used here identified consensus terminal sequences characteristic of spliceosomal introns. This finding is the first suggestion that multiple spliceosomal introns occur in ribosomal genes.
Introduction
Insertions in the nuclear small-subunit (SSU) rRNA gene have been described for many eukaryotic organisms. Such insertions occur in several discrete conserved positions (Gargas, DePriest, and Taylor 1995
) and can be optional even at the population level (DePriest 1993
; Crespo et al. 1997
). Most SSU rRNA introns so far reported are typically longer than 180 bp and have been identified as group I introns on the basis of their conserved structures (Cech 1988
; Gargas, DePriest, and Taylor 1995
). The existence of introns and their distribution among organisms and genes have given rise to a number of theories of intron origin and evolution. These theories can be grouped around two main theoretical branches: the exon theory of genes, or early-intron theory (Blake 1978
; Doolittle 1978
), and the insertional theory of introns, or late-intron theory (Rogers 1989
; Cavalier-Smith 1991
).
In recent years, there have been reports of small a number of insertions of less than 100 bp at several points in the nuclear SSU rDNA of different groups of filamentous ascomycete fungi (Rogers et al. 1993
; Gargas, DePriest, and Taylor 1995
; Grube, Gargas, and DePriest 1996
; Stenroos and DePriest 1998
; Winka, Alberg, and Eriksson 1998
). These short introns lack the necessary sequence and conserved structure to be considered group I introns, but in some cases sequences interpreted as internal guide sequences (IGSs) for splicing, similar to those in group I introns, have been described. These small insertions have been described at sites where group I introns have been found in other organisms (Gargas, DePriest, and Taylor 1995
). Site conservation, together with the IGSs, has been used to interpret the short insertions from the ascomycetes Arthonia lapidicola, Phialophora americana, and Porpidia crustulata as degenerate group I introns (Gargas, DePriest, and Taylor 1995
; Grube, Gargas, and DePriest 1996
). The small intron sequences from these fungi can be folded to give a secondary structure similar to that of the P1 and P2 stems of group I introns. Similar small insertions from Cladonia, Stereocaulon, Cladina, and Physcia are also considered to contain the two-stem structure, but without an IGS sequence (Stenroos and DePriest 1998
). One problem in accurately determining the origin of such small insertions is the short length of the sequences involved, along with the small number of such observations available.
The sequences of 24 new small insertions from the nuclear SSU rDNA of specimens of the lichen-forming ascomycete fungi Physconia were obtained. These were compared with the 12 previously reported similar short insertions (Rogers et al. 1993
; Gargas, DePriest, and Taylor 1995
; Grube, Gargas, and DePriest 1996
; Stenroos and DePriest 1998
; Winka, Ahlberg, and Eriksson 1998
). The aim of the program was to identify any conserved sequences or structures that could provide indications as to the origin and significance of small insertions and the surrounding rRNA gene block.
Materials and Methods
A cetyl-trimethyl ammonium bromide (CTAB)based extraction method (Cubero et al. 1999
) was used to extract the total DNA from samples of Physconia detersa, Physconia distorta, Physconia elegantula, Physconia enteroxantha, Physconia grisea, Physconia perisidiosa, Physconia subpulverulenta, and Physconia venusta (collection details are indicated in table 1
). Small portions of sample were ground in liquid nitrogen and incubated in CTAB extraction buffer (1% CTAB, 1 M NaCl, 100 mM Tris-HCl [pH 8.0], 20 mM EDTA, 1% polyvinyl polypyrrolidone) at 65°C. After two treatments with chloroform : iso-amylalcohol (CI, 24:1), DNA was precipitated by addition of two volumes of CTAB precipitation buffer (1% CTAB, 50 mM Tris-HCl [pH 8.0], 10 mM EDTA, 40 mM NaCl) and resuspended in 1.2 M NaCl. After a second treatment with CI, DNA was precipitated by addition of 0.6 volumes of isopropanol and resuspended in water.
|
PCR products were purified through commercial purification columns (Biotools) and sequenced using the PRISM Ready Reaction DyeDeoxy Terminator Cycle Sequencing kit and a 373 Applied Biosystems automated sequencer. Primers used for sequencing reaction were those used in the PCR and the internal primers nu-SSU-0497-3' and nu-SSU-1465-3' (Gargas and DePriest 1996
).
DNA sequences were assembled and manually aligned. Insertion sequences were detected by comparison with existing Physconia sequences known to lack introns and by comparison with the Saccharomyces cerevisiae nuclear SSU rDNA sequence (GenBank accession number J01353). Insertions were named according to the fungal sample and the insertion position with reference to the corresponding position in Escherichia coli (GenBank accession number J01695). Sequences obtained were submitted to the NCBI database under the accession numbers given in table 1 .
The specimens of P. distorta and P. perisidiosa were used for RNA isolation and amplification of the 18S rRNA to check for intron removal. Total RNA was isolated through RNA purification columns (Biotools). cDNA of the 18S rRNA was obtained by RT-PCR (Promega RT-PCR system). Product size was obtained by comparison with a 50-bp DNA ladder.
The hypothesis that insertions were degenerate group I introns was tested by obtaining base-pairing matrices from each alignment using the program DCSE (De Rijk and De Wachter 1993
). This method searches for compensatory segments within the alignment and identifies compensatory mutations which can indicate putative conserved secondary structures.
The sequences of 12 previously reported short insertions from fungal nuclear SSU rDNA (Rogers et al. 1993
; Gargas, DePriest, and Taylor 1995
; Grube, Gargas, and DePriest 1996
; Stenroos and DePriest 1998
) were compared with those obtained in this study in an attempt to identify common conserved sequences and structures. The additional fungal species and GenBank accession numbers for the small insertions were as follows: Lecanora dispersa, L37734; Porpidia crustulata, L37735; Cladonia bellidiflora, U64682; Cladonia diversa, U64683; Cladina portentosa, U64684; Stereocaulon taeniarum, U64687; Stereocaulo tomentosum, U64688; Pertusaria trachythallina, U64685; Graphis scripta, AF038878; and Physcia aipolia, U64686. The Phialophora americana insertion sequence was obtained from Rogers et al. (1993)
, and the Arthonia lapidicola insertion sequence was obtained from Grube, Gargas, and DePriest (1996)
.
Results and Discussion
Total length of nuclear SSU rDNA varied from approximately 1,720 to 3,200 bp. The product from P. subpulverulenta was 1,720 bp and did not contain any insertions. This sequence was therefore used as an initial guide to align and identify putative insertion elements. Fifty insertions, located at 15 different positions, were found in the remaining samples. Twenty-six of these insertions were typical group I introns (data not shown). Twenty-four small insertions (5070 bp) were too short to be identified as group I introns, and these were recovered from seven different positions in different species (table 1 ).
The elimination of the insertions from the mature 18S RNA was verified by reverse transcription and amplification of the 18S cDNA in P. distorta and P. perisidiosa. Six small insertions and five group I introns or three small insertions and six group I introns had previously been detected in these specimens. The RT-PCR products obtained in both cases were 1,720 bp. As this was the expected size of the amplified segment when introns were not present, it seems likely that both the group I introns and the small insertions were removed during the maturation process.
The exact site of small insertion could not always be determined solely by comparison with the insertion-lacking sequence of P. subpulverulenta. This uncertainty was due to the difficulty in identifying the 5' and 3' termini of the insertion within the sequence. This can be illustrated by considering a theoretical sequence containing an insertion CACGG(N)nAGCG and the corresponding intron-lacking sequence CACGCG. In this case, there are two different points at which the insertion could be located depending on the identity of the first and last base in the insertion sequence: CAC[GG(N)nA]GCG or CACG[G(N)nAG]CG. In order to accurately align the short insertions, they were all considered with both alternative termini. The main criterion used in this alignment was to conserve the insertion sequence wherever possible. Under this criterion, all short insertions were found to have GT at their 5' termini and AG at their 3' termini. The final intron positions are given in table 1 , and these correspond to the positions after nucleotides 296, 297, 330, 513, 673, 943, and 1129 in E. coli.
There is considerable sequence divergence between the different insertion sequences. Insertions from different locations in the 18S rDNA are not alignable, except for those at locations 296 and 297, which are similar in sequence (50% similar positions between all of the sequences and 60%83% similarity between pairs of sequences) and alignable. These insertions are interpreted as the result of an insertion sliding due to compensatory mutations at both intron-exon junctions (fig. 1 ), which is in agreement with the phylogenetic relationship between the species of the genus based on internal transcribed spacer sequences (unpublished data). For that reason, insertions at position 296 and insertions at position 297 are considered homologous in this work and are analyzed together.
|
|
Some other features observed in the Physconia insertions are also characteristic of spliceosomal introns. These include the position of the putative branch point from the 3' end of the sequence, the high frequency of G at position 5 of the intron sequence (Cech 1983
), and the presence of a pyrimidine immediately before the terminal AG motif (Jeffeys and Flavell 1977
). Although these features are not all universal for spliceosomal introns, they occur at high frequencies in messenger RNA introns and are usually included in consensus descriptions.
A general consensus description for the small introns in Physconia nuclear SSU rDNA can be given as GT(A/T)N(G/A) ... PurCTAAC ... PyrAG, where G/A represents a purine that is usually G but occasionally A. The sequence and motif conservation described here provides evidence to support the hypothesis that the small insertions found in the SSU rDNA gene of the Physconia samples are putative spliceosomal introns. They possess conserved functional features required for external splicing in messenger preRNAs. This is the first suggestion that multiple insertions in the ribosomal genes are putative spliceosomal introns.
Ribosomal RNA genes generally lack introns, but in some organisms these genes contain large insertions which have been identified as group I introns (for review, see Gargas et al. 1995
; Johansen, Muscarella, and Vogt 1996
). Prior to this study, only 12 nuclear SSU rDNA small insertions had been described, all of them in ascomycete fungi (fig. 3
). Of these, nine have been found in different species of the lichen-forming family lecanorales (Gargas, DePriest, and Taylor 1995
; Stenroos and DePriest 1998
), one has been found in the lichen-forming fungi A. lapidicola (arthoniales; Grube, Gargas, and DePriest 1996
) and G. scripta (ostropales; Winka, Ahlberg, and Eriksson 1998
), and one has been found in the parasitic fungus P. americana (dothideales; Rogers et al. 1993
). The latter is the only small insertion in the nuclear SSU rDNA described as being a putative spliceosomal type. However, the possibility that further small SSU rRNA insertions may be spliceosomal introns has not been generally discussed. In practice, most workers have placed prime importance on potential secondary structures and insertion site, rather than on sequence conservation. As a result, there have been few, if any, comparisons of possible alternative insertion sites and little critical examination of sequence features.
|
Putative Secondary Structures
In order to consider the possibility of conserved secondary structures among the Physconia small insertions, complementary segments were searched in the alignment. As the hypothesis of "degenerate group I introns" involves interactions between internal sequences of the intron- and exon-bordering sequences (Grube, Gargas, and DePriest 1996
), 10 nucleotides before and after the insertion were included in the analysis.
The results obtained did not show any evidence for the existence of conserved or similar structures in the different Physconia small introns. Attempts to generate two-stem secondary structures or detect an IGS sequence as proposed for the small introns of A. lapidicola and P. crustulata (Grube, Gargas, and DePriest 1996
) failed to give positive results, and although two-stem secondary structures were obtained by hand for some particular intron sequences, they were not in agreement with the conservation patterns shown in the alignments.
Distribution and Evolution of 18S rDNA Small Introns
Traditionally, ribosomal RNA genes are thought to lack spliceosomal introns of the type found in messenger RNAs. However, we have demonstrated that in some organisms there are introns in the 18S rDNA which clearly resemble those of protein spliceosomal introns. Currently, there is not enough information to establish a hypothesis about the origin of these ribosomal introns. However, any future hypothesis will need to consider several points:
The results of this study show that maximizing terminal sequence conservation provides sufficient evidence to identify the small insertion from the SSU rRNA genes of filamentous fungi as putative spliceosomal introns. However, further knowledge of the wider distribution of these ribosomal small insertions is required before a full hypothesis can be developed.
Acknowledgements
We thank F. García-Olmedo, D. L. Hawksworth, J. Valcárcel, P. Vargas, and two anonymous reviewers for their useful comments. This work was supported by grants CAM-07B-0014-1997-97 from the Comunidad Autónoma de Madrid and APC-960082 from the Spanish Ministerio de Educación y Cultura (MEC). O.F.C. was supported by an FPI fellowship from the MEC.
Footnotes
Thomas Eickbush, Reviewing Editor
2 Keywords: ribosomal insertion
spliceosomal intron
fungal insertion
group I intron
degenerate introns
Physconia.
1 Address for correspondence and reprints: Oscar F. Cubero, Departamento de Biología Vegetal II, Facultad de Farmacia, Universidad Complutense de Madrid, Madrid 28040, Spain.E-mail: oscarhfc{at}retemail.es
literature cited
Blake, C. C. F. 1978. Do genes-in-pieces imply proteins-in-pieces? Nature 273:267.
Cavalier-Smith, T. 1991. Intron phylogeny: a new hypothesis. TIG 7:145148.
Cech, T. R. 1983. RNA splicing: three themes with variations. Cell 34:713716.
. 1988. Conserved sequences and structures of group I introns: building an active site for RNA catalysisa review. Gene 73:258271.
Chen, W., L. E. Gray, and C. R. Grau. 1998. Phialophora gregata f. sp. adzukicola small subunit ribosomal RNA gene, partial sequence. GenBank accession number AF056487.
Crespo, A., P. D. Bridge, O. F. Cubero, and D. L. Hawksworth. 1997. Determination of genotypic variability in the lichen-forming fungus Parmelia sulcata. Bibl. Lichenol. 68:7379.
Cubero, O. F., A. Crespo, J. Fatehi, and P. D. Bridge. 1999. DNA extraction and PCR amplification method suitable for fresh, herbarium and lichenized fungi. Plant. Syst. Evol. 216:243249.[ISI]
DePriest, P. T. 1993. Small subunit rDNA variation in a population of lichen fungi due to optional group I introns. Gene 134:6774.
De Rijk, P., and R. De Wachter. 1993. DCSE, an interactive tool for sequence alignment and secondary structure research. Version 2.54.
Doolittle, W. F. 1978. Genes in pieces: were they ever together? Nature 272:581582.
Gargas, A., and P. T. DePriest. 1996. Polymerase chain reaction PCR primers for amplifying and sequencing nuclear 18S rDNA from lichenized fungi. Mycologia 885:745748.
Gargas, A., P. T. DePriest, M. Grube, and A. Tehler. 1995. Multiple origins of lichen symbioses in fungi suggested by SSU rDNA phylogeny. Science 268:14921495.
Gargas, A., P. T. DePriest, and J. W. Taylor. 1995. Position of multiple insertions in the SSU rDNA of lichen-forming fungi. Mol. Biol. Evol. 12:208218.[Abstract]
Grube, M., A. Gargas, and P. T. DePriest. 1996. A small insertion in the SSU rDNA of the lichen fungus Arthonia lapidicola is a degenerated group I intron. Curr. Genet. 29:582586.[ISI][Medline]
Jeffeys, A. J., and R. A. Flavell. 1977. The rabbit B-globin gene contains a large insert in the coding sequence. Cell 12:10971108.
Johansen, S., D. E. Muscarella, and V. M. Vogt. 1996. Insertion elements in ribosomal DNA. Pp. 89108 in R. A. Zimmermann and A. E. Dahlberg, eds. Ribosomal RNA. Structure, evolution, processing and function in protein biosynthesis. CRC Press, Boca Raton, Fla.
Langford, C. J., and D. Gallwitz. 1983. Evidence for an intron-contained sequence required for splicing of yeast RNA polymerase II transcripts. Cell 33:519527.
Lerner, M. R., J. A. Boyle, S. M. Mount, S. L. Wolin, and J. A. Steitz. 1980. Are snRNPs involved in splicing? Nature 238:220224.
Mount, S. M., I. Petterson, M. Hinterberger, A. Karmas, and J. Steitz. 1983. The U1 small nuclear RNA-protein complex selectively binds a 5' splice site in vitro. Cell 33:509518.
Raymond, B. C., and M. Rosbash. 1992. Yeast pre mRNA splicing. Pp. 143192 in J. Elizabeth, J. R. Pringle, and J. B. Broach, eds. The molecular and cellular biology of the yeast Saccharomyces: gene expression. Cold Spring Harbor Laboratory Press, New York.
Reich, C. I., R. W. Van Hoy, G. L. Porter, and J. A. Wise. 1992. Mutations at the 3' splice site can be suppressed by compensatory base changes in U1 snRNA in fission yeast. Cell 69:11591169.
Rogers, J. H. 1989. How were introns inserted into nuclear genes? TIG 5:213216.
Rogers, J., and R. Wall. 1980. A mechanism for RNA splicing. Proc. Natl. Acad. Sci. USA 77:18771879.
Rogers, S. O., Z. H. Yan, M. Shinohara, K. F. Lobuglio, and C. J. K. Wang. 1993. Messenger RNA intron in the nuclear 18S ribosomal RNA gene of deuteromycetes. Curr. Genet. 23:338342.[ISI][Medline]
Schneider, T. D., and R. M. Stephens. 1990. Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 18:60976100.[Abstract]
Stenroos, S., and P. T. DePriest. 1998. Small insertions in a shared position in the SSU rDNA of Lecanorales lichen-forming Ascomycetes. Curr. Genet. 33:124130.[ISI][Medline]
Stephens, R. M., and T. D. Schneider. 1992. Features of spliceosome evolution and function inferred from an analysis of the information at human splice sites. J. Mol. Biol. 226:11241136.
Winka, K., C. Ahlberg, and O. E. Eriksson. 1998. Are there lichenized Ostropales? Lichenologist 30:455462.
Wu, J., and J. L. Manley. 1989. Mammalian pre-mRNA branch site selection by U2 snRNP involves base pairings. Genes Dev. 3:15531561.[Abstract]
Zhuang, Y., and A. M. Weiner. 1989. A compensatory base change in human U2 snRNA can suppress a branch site mutation. Genes Dev. 3:15.