Association of intron phases with conservation at splice site sequences and evolution of spliceosomal introns.

How exon-intron structures of eukaryotic genes evolved under various evolutionary forces remains unknown. The phases of spliceosomal introns (the placement of introns with respect to reading frame) provide an opportunity to approach this question. When a large number of nuclear introns in protein-coding genes were analyzed, it was found that most introns were of phase 0, which keeps codons intact. We found that the phase distribution of spliceosomal introns is strongly correlated with the sequence conservation of splice signals in exons; the relatively underrepresented phase 2 introns are associated with the lowest conservation, the relatively overrepresented phase 0 introns display the highest conservation, and phase 1 introns are intermediate. Given the detrimental effect of mutations in exon sequences near splice sites as found in molecular experiments, the underrepresentation of phase 2 introns may be the result of deleterious-mutation-driven intron loss, suggesting a possible genetic mechanism for the evolution of intron-exon structures.     

USING RDNA GENES TO UNDERSTAND THE ORIGIN AND EVOLUTION OF SPLICEOSOMAL AND GROUP I INTRONS

Ribosomal DNA genes in lichen algae and lichen fungi are astonishingly rich in spliceosomal and group I introns. We use phylogenetic, secondary structure, and biochemical analyses to understand the evolution of these introns. Despite the widespread distribution of spliceosomal introns in nuclear pre-mRNA genes, their general mechanism of origin remains an open question because few proven cases of recent and pervasive intron origin have been documented. The lichen introns are valuable in this respect because they are undoubtedly of a "recent" origin and limited to the Euascomycetes. Our analyses suggest that rDNA spliceosomal introns have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into r RNAs. We propose that the spliceosome itself (and not an external agent; e.g. transposable elements, group II introns) has given rise to the introns. The rDNA introns are found most often between the flanking sequence G (78%) - intron-G (72%), and their clustered positions on secondary structures suggest that particular r RNA regions are preferred sites (i.e., proto-splice sites) for insertion. Mapping of intron positions on the newly available tertiary structures show that they are found most often in exposed regions of the ribosomes. This again is consistent with an intron origin through reverse-splicing. Remarkably, the distribution and phylogenetic relationships of most group I introns in nuclear rDNA genes are also consistent with a reverse-splicing origin. These data underline the value of lichens as a model system for understanding intron origin and stress the importance of RNA-level processes in the spread of these sequences in nuclear coding regions.     

Group I introns in the liverwort mitochondrial genome: the gene coding for subunit 1 of cytochrome oxidase shares five intron positions with its fungal counterparts.

The complete nucleotide sequence of the mitochondrial DNA (mtDNA) from a liverwort, Marchantia polymorpha, contains thirty-two introns. Twenty-five of these introns possess the characteristic secondary structures and consensus sequences of group II introns. The remaining seven are group I introns, six of which happen to interrupt the gene coding for subunit 1 of cytochrome oxidase (cox1). Interestingly, the insertion sites of one group II and four group I introns in the cox1 gene coincide with those of the respective fungal mitochondrial interns. Moreover, comparison of the four group I introns with their fungal counterparts shows that group I introns inserted at identical genomic sites in different organisms are indeed related to one another, in terms of the peptide sequences generated from the complete or fragmental ORFs encoded by these introns. At the same time, the liverwort introns turned out to be more divergent from their fungal cognates than the latter are from one another. We therefore conclude that vertical transmission from a common ancestor organism is the simplest explanation for the presence of cognate introns in liverwort and fungal mitochondrial genomes.     

酵母基因内含子中二聚体寡核苷酸转录调控位点的统计分析

对高效和低效转录酵母基因内含子序列中寡核苷酸的出现频率进行对照分析, 结果显示高效和低效内含子序列的结构有差异, 而且高效转录内含子序列含有较多潜在的转录因子结合位点. 观察实验获得的转录调控位点, 发现许多调控位点不是相邻接的寡核苷酸,而是由一对保守寡核苷酸构成, 这对寡核苷酸被一段长度固定的非保守区域间隔开. 于是对此形式的二聚体寡核苷酸(dyad)在高效和低效内含子序列中出现的频率进行统计比较分析,抽提出在高效内含子组出现的频率显著高于在低效内含子组出现频率的二聚体寡核苷酸, 分析这些二聚体寡核苷酸在两组内含子序列中的分布特征, 并对照实验结果, 这些二聚体寡核苷酸可能与基因转录的正调控有关.     

Conservation of RNA secondary structures in two intron families including mitochondrial-, chloroplast- and nuclear-encoded members

Two families of fungal mitochondrial introns that include all known sequences have been recognized. These families are now extended to incorporate a plant mitochondrial intron and several introns in chloroplast- and nuclear-encoded rRNA and tRNA precursors. Members of the same family share distinctive sequence stretches and a number of potential RNA secondary structures that would bring these stretches and the intron-exon junctions into relatively close proximity. Using several of these introns which have been extensively studied by either biochemical or genetic means, an attempt is made to integrate the available data into a common picture.     

Sequence organization of the chloroplast ribosomal spacer ofSpinacia oleracea including the 3′ end of the 16S rRNA and the 5′ end of the 23S rRNA

The 2201-bp spacer between the chloroplast ribosomal 16S and 23S genes ofSpinacia oleracea was sequenced. It contains the genes of the tRNA^Ile (GAU) and tRNA^Ala (UGC) which are both interrupted by introns of respectively 728 and 816 bp. These introns belong to the class II according to the classfication of Michel and Dujon [17]. Comparison of the rDNA spacer sequence of maize, tobacco and spinach indicates that no conserved polypeptide is encoded within the introns of the two tRNA genes and that the two main insertions/deletions between the three plants are located within two loops of the class II introns secondary structure, which is therefore conserved. Based on the sequence complementarity observed between the upstream and downstream parts, of the 16S and 23S rRNA genes, RNase III-like secondary structures involved in the processing of the rRNA precursor are proposed.     

Birth of new spliceosomal introns in fungi by multiplication of introner-like elements

Spliceosomal introns are noncoding sequences that separate exons in eukaryotic genes and are removed from pre-messenger RNAs by the splicing machinery. Their origin has remained a mystery in biology since their discovery because intron gains seem to be infrequent in many eukaryotic lineages. Although a few recent intron gains have been reported, none of the proposed gain mechanisms can convincingly explain the high number of introns in present-day eukaryotic genomes. Here we report on particular spliceosomal introns that share high sequence similarity and are reminiscent of introner elements. These elements multiplied in unrelated genes of six fungal genomes and account for the vast majority of intron gains in these fungal species. Such introner-like elements (ILEs) contain all typical characteristics of regular spliceosomal introns (RSIs) but are longer and predicted to harbor more stable secondary structures. However, dating of multiplication events showed that they degenerate in sequence and length within 100,000 years to eventually become indistinguishable from RSIs. We suggest that ILEs not only account for intron gains in six fungi but also in ancestral eukaryotes to give rise to most RSIs by a yet unknown multiplication mechanism.     

Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses

Most eukaryotes have at least some genes interrupted by introns. While it is well accepted that introns were already present at moderate density in the last eukaryote common ancestor, the conspicuous diversity of intron density among genomes suggests a complex evolutionary history, with marked differences between phyla. The question of the rates of intron gains and loss in the course of evolution and factors influencing them remains controversial. We have investigated a single gene family, alpha-amylase, in 55 species covering a variety of animal phyla. Comparison of intron positions across phyla suggests a complex history, with a likely ancestral intronless gene undergoing frequent intron loss and gain, leading to extant intron/exon structures that are highly variable, even among species from the same phylum. Because introns are known to play no regulatory role in this gene and there is no alternative splicing, the structural differences may be interpreted more easily: intron positions, sizes, losses or gains may be more likely related to factors linked to splicing mechanisms and requirements, and to recognition of introns and exons, or to more extrinsic factors, such as life cycle and population size. We have shown that intron losses outnumbered gains in recent periods, but that "resets" of intron positions occurred at the origin of several phyla, including vertebrates. Rates of gain and loss appear to be positively correlated. No phase preference was found. We also found evidence for parallel gains and for intron sliding. Presence of introns at given positions was correlated to a strong protosplice consensus sequence AG/G, which was much weaker in the absence of intron. In contrast, recent intron insertions were not associated with a specific sequence. In animal Amy genes, population size and generation time seem to have played only minor roles in shaping gene structures.     

RNA splicing in the T-even bacteriophage

Group 1 introns, first demonstrated in the nuclear large rRNA of Tetrahymena thermophila and subsequently in many yeast, fungal mitochondrial, and chloroplast precursor RNAs, are capable of intron excision and exon ligation in vitro, although this process occurs much more rapidly in vivo. The discovery and characterization of a similar intron in the T4 phage thymidylate synthase gene (td) led to the finding of additional group 1 introns in other T4 genes and in genes of the related T2 and T6 phages. Because protein factors are not required in the splicing of group 1 introns in vitro, it has been postulated that the precursor RNA can assume a critical conformation enabling it to undergo site-specific autocatalytic cleavage and ligation (self-splicing). By means of site-directed mutation, it has been shown unequivocally that several sequence elements in the Tetrahymena rRNA intron are involved in the formation of base-paired stem structures that are essential for the self-splicing process. These sequence elements have been demonstrated in other eukaryotic group 1 introns, as well as in the td intron. In this brief review we shall describe the biochemical and structural properties of the td intron in relation to other newly found phage introns. The interesting implications arising from these revelations will also be discussed.     

Evolution of the intron-exon structure of eukaryotic genes

The origin and evolution of intron-exon structures continue to be controversial topics. Two alternative theories, the ‘exon theory of genes’ and the ‘insertional theory of introns’, debate the presence or absence of introns in primordial genes. Both sides of the argument have focused on the positions of introns with respect to protein and gene structures. A new approach has emerged in the study of the evolution of intron-exon structures: a population analysis of genes. One example is the statistical analysis of intron phases — the position of introns within or between codons. This analysis detected a significant signal of exon shuffling in the DNA sequence database containing both ancient and modern exon sequences: intron phase correlations, that is, the association together within genes of introns of the same phase. The results of this analysis suggest that exon shuffling played an important role in the origin of both ancient and modern genes.     

DNA sequence and secondary structures of the large subunit rRNA coding regions and its two class I introns of mitochondrial DNA fromPodospora anserina

Summary DNA sequence analysis has shown that the gene coding for the mitochondrial (mt) large subunit ribosomal RNA (rRNA) fromPodospora anserina is interrupted by two class I introns. The coding region for the large subunit rRNA itself is 3715 bp and the two introns are 1544 (r1) and 2404 (r2) bp in length. Secondary structure models for the large subunit rRNA were constructed and compared with the equivalent structure fromEscherichia coli 23S rRNA. The two structures were remarkably similar despite an 800-base difference in length. The additional bases in theP. anserina rRNA appear to be mostly in unstructured regions in the 3 part of the RNA. Secondary structure models for the two introns show striking similarities with each other as well as with the intron models from the equivalent introns inSaccharomyces cerevisiae, Neurospora crassa, andAspergillus nidulans. The long open reading frames in each intron are different from each other, however, and the nucleotide sequence similarity diverges as it proceeds away from the core structure. Each intron is located within regions of the large subunit rRNA gene that are highly conserved in both sequence and structure. Computer analysis showed that the open reading frame for intron r1 contained a common maturase-like polypeptide. The open reading frames of intron r2 apeared to be chimeric, displaying high sequence similarity with the open reading frames in the r1 and ATPase 6 introns ofN. crassa.     

Widespread occurrence of spliceosomal introns in the rDNA genes of ascomycetes

Spliceosomal (pre-mRNA) introns have previously been found in eukaryotic protein-coding genes, in the small nuclear RNAs of some fungi, and in the small- and large-subunit ribosomal DNA genes of a limited number of ascomycetes. How the majority of these introns originate remains an open question because few proven cases of recent and pervasive intron origin have been documented. We report here the widespread occurrence of spliceosomal introns (69 introns at 27 different sites) in the small- and large-subunit nuclear-encoded rDNA of lichen-forming and free-living members of the Ascomycota. Our analyses suggest that these spliceosomal introns are of relatively recent origin, i.e., within the Euascomycetes, and have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into rRNAs. The spliceosome itself, and not an external agent (e.g., transposable elements, group II introns), may have given rise to these introns. A nonrandom sequence pattern was found at sites flanking the rRNA spliceosomal introns. This pattern (AG-intron-G) closely resembles the proto-splice site (MAG-intron-R) postulated for intron insertions in pre-mRNA genes. The clustered positions of spliceosomal introns on secondary structures suggest that particular rRNA regions are preferred sites for insertion through reverse-splicing.     

The in vivo use of alternate 3'-splice sites in group I introns.

Alternative splicing of group I introns has been postulated as a possible mechanism that would ensure the translation of proteins encoded into intronic open reading frames, discontinuous with the upstream exon and lacking an initiation signal. Alternate splice sites were previously depicted according to secondary structures of several group I introns. We present here strong evidence that, in the case of Podospora anserina nad 1-i4 and cox1-i7 mitochondrial introns, alternative splicing events do occur in vivo. Indeed, by PCR experiments we have detected molecules whose sequence is precisely that expected if the predicted alternate 3'-splice sites were used.     

The unusual 23S rRNA gene of Coxiella burnetii: two self-splicing group I introns flank a 34-base-pair exon, and one element lacks the canonical omegaG

We describe the presence and characteristics of two self-splicing group I introns in the sole 23S rRNA gene of Coxiella burnetii. The two group I introns, Cbu.L1917 and Cbu.L1951, are inserted at sites 1917 and 1951 (Escherichia coli numbering), respectively, in the 23S rRNA gene of C. burnetii. Both introns were found to be self-splicing in vivo and in vitro even though the terminal nucleotide of Cbu.L1917 is adenine and not the canonical conserved guanine, termed OmegaG, found in Cbu.L1951 and all other group I introns described to date. Predicted secondary structures for both introns were constructed and revealed that Cbu.L1917 and Cbu.L1951 were group IB2 and group IA3 introns, respectively. We analyzed strains belonging to eight genomic groups of C. burnetii to determine sequence variation and the presence or absence of the elements and found both introns to be highly conserved (>/=99%) among them. Although phylogenetic analysis did not identify the specific identities of donors, it indicates that the introns were likely acquired independently; Cbu.L1917 was acquired from other bacteria like Thermotoga subterranea and Cbu.L1951 from lower eukaryotes like Acanthamoeba castellanii. We also confirmed the fragmented nature of mature 23S rRNA in C. burnetii due to the presence of an intervening sequence. The presence of three selfish elements in C. burnetii's 23S rRNA gene is very unusual for an obligate intracellular bacterium and suggests a recent shift to its current lifestyle from a previous niche with greater opportunities for lateral gene transfer.     

Genomic organization of human lactate dehydrogenase-A gene.

A human genomic clone containing the lactate dehydrogenase-A (LDH-A) gene of approx. 12 kilobases in length was isolated and characterized. The protein-coding sequence is interrupted by six introns, and the positions of these introns are at the random coil regions or near the ends of secondary structures located on the surface of the LDH-A molecule. An additional intron is present at 24 nucleotides 5' to the translation initiation codon ATG, while the 3' untranslated sequence of 565 nucleotides is not interrupted. The genomic blot analysis of human placenta DNA indicates the presence of multiple LDH-A gene-related sequences.     

Isolation and Characterization of Group II Introns from Pseudomonas alcaligenes and Pseudomonas putida

Group II introns isolated from Pseudomonas alcaligenes NCIB 9867, Pseudomonas putida NCIB 9869, and P. putida KT2440 were closely related with nucleotide sequence identities of between 87 and 96%. The genome of P. alcaligenes also harbored a truncated group II intron of 682 bp that lacks the gene for the intron-encoded protein (IEP). Unlike most bacterial group II introns, the Pseudomonas introns were found to lack the Zn domains in their IEPs, did not appear to interrupt any genes, and were located downstream of open reading frames which were adjacent to hairpin loop structures that resemble rho-independent terminators. These structures also contain the intron binding sites 1 and 2 (IBS1 and IBS2 sequences) that were required for intron target site recognition in transposition. One of the group II introns found in P. alcaligenes, Xln3, was shown to have transposed from the chromosome to the endogenous pRA2 plasmid at a site adjacent to IBS1- and IBS2-like sequences.     

The Euglena gracilis chloroplast ribulose-1,5-bisphosphate carboxylase gene. I. Complete DNA sequence and analysis of the nine intervening sequences

The nucleotide sequence of 6225 base pairs (bp) of Euglena gracilis chloroplast DNA including the complete DNA sequence of the chloroplast-encoded ribulose-1,5-bisphosphate carboxylase large subunit gene along with the flanking DNA sequences is presented. The gene is greater than 5.5 kilobase pairs in length and is organized as 10 exons coding for 475 amino acids, separated by 9 introns. The exons range in size from 45 to 438 bp, while the introns range in size from 382 to 568 bp. The introns have highly conserved boundary sequences with the consensus, 5'-N GTGTGGATTT...(intron)...TTAATTTTAT N-3'. The introns are 82-85 mol% AT, with a pronounced T greater than A greater than G greater than C base bias in the RNA-like strand. They do not appear to encode any polypeptides. In addition, the introns have a conserved sequence 30-50 bp from their 3'-ends with the consensus, 5'-TACAGTTTGAAAATGA-3'. The 5'-TACA sequence bears some homology to the 5'-end of the TACTAACA sequence found in a similar location in yeast nuclear mRNA introns. The conserved sequences of the Euglena rbcL introns may be indicative of a splicing mechanism similar to that of eucaryotic nuclear mRNA introns and group II mitochondrial introns.