Intron positions correlate with module boundaries in ancient proteins

We analyze the three-dimensional structure of proteins by a computer program that finds regions of sequence that contain module boundaries, defining a module as a segment of polypeptide chain bounded in space by a specific given distance. The program defines a set of “linker regions” that have the property that if an intron were to be placed into each linker region, the protein would be dissected into a set of modules all less than the specified diameter. We test a set of 32 proteins, all of ancient origin, and a corresponding set of 570 intron positions, to ask if there is a statistically significant excess of intron positions within the linker regions. For 28-Å modules, a standard size used historically, we find such an excess, with P < 0.003. This correlation is neither due to a compositional or sequence bias in the linker regions nor to a surface bias in intron positions. Furthermore, a subset of 20 introns, which can be putatively identified as old, lies even more explicitly within the linker regions, with P < 0.0003. Thus, there is a strong correlation between intron positions and three-dimensional structural elements of ancient proteins as expected by the introns-early approach. We then study a range of module diameters and show that, as the diameter varies, significant peaks of correlation appear for module diameters centered at 21.7, 27.6, and 32.9 Å. These preferred module diameters roughly correspond to predicted exon sizes of 15, 22, and 30 residues. Thus, there are significant correlations between introns, modules, and a quantized pattern of the lengths of polypeptide chains, which is the prediction of the “Exon Theory of Genes.”

Chloroplast protein and centrosomal genes, a tRNA intron, and odd telomeres in an unusually compact eukaryotic genome, the cryptomonad nucleomorph

Cells of several major algal groups are evolutionary chimeras of two radically different eukaryotic cells. Most of these “cells within cells” lost the nucleus of the former algal endosymbiont. But after hundreds of millions of years cryptomonads still retain the nucleus of their former red algal endosymbiont as a tiny relict organelle, the nucleomorph, which has three minute linear chromosomes, but their function and the nature of their ends have been unclear. We report extensive cryptomonad nucleomorph sequences (68.5 kb), from one end of each of the three chromosomes of Guillardia theta. Telomeres of the nucleomorph chromosomes differ dramatically from those of other eukaryotes, being repeats of the 23-mer sequence (AG)7AAG6A, not a typical hexamer (commonly TTAGGG). The subterminal regions comprising the rRNA cistrons and one protein-coding gene are exactly repeated at all three chromosome ends. Gene density (one per 0.8 kb) is the highest for any cellular genome. None of the 38 protein-coding genes has spliceosomal introns, in marked contrast to the chlorarachniophyte nucleomorph. Most identified nucleomorph genes are for gene expression or protein degradation; histone, tubulin, and putatively centrosomal ranbpm genes are probably important for chromosome segregation. No genes for primary or secondary metabolism have been found. Two of the three tRNA genes have introns, one in a hitherto undescribed location. Intergenic regions are exceptionally short; three genes transcribed by two different RNA polymerases overlap their neighbors. The reported sequences encode two essential chloroplast proteins, FtsZ and rubredoxin, thus explaining why cryptomonad nucleomorphs persist.

Phylogenetic relationships among group II intron ORFs

Group II introns are widely believed to have been ancestors of spliceosomal introns, yet little is known about their own evolutionary history. In order to address the evolution of mobile group II introns, we have compiled 71 open reading frames (ORFs) related to group II intron reverse transcriptases and subjected their derived amino acid sequences to phylogenetic analysis. The phylogenetic tree was rooted with reverse transcriptases (RTs) of non-long terminal repeat retroelements, and the inferred phylogeny reveals two major clusters which we term the mitochondrial and chloroplast-like lineages. Bacterial ORFs are mainly positioned at the bases of the two lineages but with weak bootstrap support. The data give an overview of an apparently high degree of horizontal transfer of group II intron ORFs, mostly among related organisms but also between organelles and bacteria. The Zn domain (nuclease) and YADD motif (RT active site) were lost multiple times during evolution. Differences in domain structures suggest that the oldest ORFs were concise, while the ORF in the mitochondrial lineage subsequently expanded in three locations. The data are consistent with a bacterial origin for mobile group II introns.

Comparison of intron-containing and intron-lacking human genes elucidates putative exonic splicing enhancers

Of the rules used by the splicing machinery to precisely determine intron–exon boundaries only a fraction is known. Recent evidence suggests that specific short sequences within exons help in defining these boundaries. Such sequences are known as exonic splicing enhancers (ESE). A possible bioinformatical approach to studying ESE sequences is to compare genes that harbor introns with genes that do not. For this purpose two non-redundant samples of 719 intron-containing and 63 intron-lacking human genes were created. We performed a statistical analysis on these datasets of intron-containing and intron-lacking human coding sequences and found a statistically significant difference (P = 0.01) between these samples in terms of 5–6mer oligonucleotide distributions. The difference is not created by a few strong signals present in the majority of exons, but rather by the accumulation of multiple weak signals through small variations in codon frequencies, codon biases and context-dependent codon biases between the samples. A list of putative novel human splicing regulation sequences has been elucidated by our analysis.

Intron insertion facilitates amplification of cloned virus cDNA in Escherichia coli while biological activity is reestablished after transcription in vivo.

Insertion of introns into cloned cDNA of two isolates of the plant potyvirus pea seedborne mosaic virus facilitated plasmid amplification in Escherichia coli. Multiple stop codons in the inserted introns interrupted the open reading frame of the virus cDNA, thereby terminating undesired translation of virus proteins in E. coli. Plasmids containing the full-length virus sequences, placed under control of the cauliflower mosaic virus 35S promoter and the nopaline synthase termination signal, were stable and easy to amplify in E. coli if one or more introns were inserted into the virus sequence. These plasmids were infectious when inoculated mechanically onto Pisum sativum leaves. Examination of the cDNA-derived viruses confirmed that intron splicing of in vivo transcribed pre-mRNA had occurred as predicted, reestablishing the virus genome sequences. Symptom development and virus accumulation of the cDNA derived viruses and parental viruses were identical. It is proposed that intron insertion can be used to facilitate manipulation and amplification of cloned DNA fragments that are unstable in, or toxic to, E. coli. When transcribed in vivo in eukaryotic cells, the introns will be eliminated from the sequence and will not interfere with further analysis of protein expression or virus infection.

Alternatively spliced mRNAs predicted to yield frame-shift proteins and stable intron 1 RNAs of the herpes simplex virus 1 regulatory gene alpha 0 accumulate in the cytoplasm of infected cells.

The infected cell protein no. 0 (ICP0), the product of the alpha 0 gene, and an important herpes simplex virus 1 regulatory protein is encoded by three exons. We report that intron 1 forms a family of four stable nonpolyadenylylated cytoplasmic RNAs sharing a common 5' end but differing in 3' ends. The 5' and 3' ends correspond to the accepted splice donor and four splice acceptor sites within the mapped intron domain. The most distant splice acceptor site yields the mRNA encoding the 775-aa protein known as ICP0. The mRNAs resulting from the use of alternative splice acceptor sites were also present in the cytoplasm of infected cells and would be predicted to encode proteins of 152 (ICP0-B), 87 (ICP0-C), and 90 (ICP0-D) amino acids, respectively. Both the stability of the alpha 0 mRNA and the utilization of at least one splice acceptor site was regulated by ICP22 and or US1.5 protein inasmuch as cells infected with a mutant from which these genes had been deleted accumulated smaller amounts of alpha 0 mRNA than would be predicted from the amounts of accumulated intron RNAs. In addition, one splice acceptor site was at best underutilized. These results indicate that both the splicing pattern and longevity of alpha 0 mRNA are regulated. These and other recent examples indicate that herpes simplex virus 1 regulates its own gene expression and that of the infected cells through control of mRNA splicing and longevity.

Two mitochondrial group I introns in a metazoan, the sea anemone Metridium senile: one intron contains genes for subunits 1 and 3 of NADH dehydrogenase.

Mitochondrial genes for cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 5 (ND5) of the sea anemone Metridium senile (phylum Cnidaria) each contain a group I intron. This is in contrast to the reported absence of introns in all other metazoan mtDNAs so far examined. The ND5 intron is unusual in that it ends with A and contains two genes (ND1 and ND3) encoding additional subunits of NADH dehydrogenase. Correctly excised ND5 introns are not circularized but are precisely cleaved near their 3' ends and polyadenylylated to provide bicistronic transcripts of ND1 and ND3. COI introns, which encode a putative homing endonuclease, circularize, but in a way that retains the entire genome-encoded intron sequence (other group I introns are circularized with loss of a short segment of the intron 5' end). Introns were detected in the COI and ND5 genes of other sea anemones, but not in the COI and ND5 genes of other cnidarians. This suggests that the sea anemone mitochondrial introns may have been acquired relatively recently.

Intercellular mobility and homing of an archaeal rDNA intron confers a selective advantage over intron- cells of Sulfolobus acidocaldarius.

Some intron-containing rRNA genes of archaea encode homing-type endonucleases, which facilitate intron insertion at homologous sites in intron- alleles. These archaeal rRNA genes, in contrast to their eukaryotic counterparts, are present in single copies per cell, which precludes intron homing within one cell. However, given the highly conserved nature of the sequences flanking the intron, homing may occur in intron- rRNA genes of other archaeal cells. To test whether this occurs, the intron-containing 23S rRNA gene of the archaeal hyperthermophile Desulfurococcus mobilis, carried on nonreplicating bacterial vectors, was electroporated into an intron- culture of Sulfolobus acidocaldarius. PCR experiments demonstrated that the intron underwent homing and spread through the culture. By using a double drug-resistant mutant of S. acidocaldarius, it was shown that spreading resulted partly from a selective advantage of intron+ cells and partly from intercellular mobility of the intron and homing.

Intron phase correlations and the evolution of the intron/exon structure of genes.

Two issues in the evolution of the intron/exon structure of genes are the role of exon shuffling and the origin of introns. Using a large data base of eukaryotic intron-containing genes, we have found that there are correlations between intron phases leading to an excess of symmetric exons and symmetric exon sets. We interpret these excesses as manifestations of exon shuffling and make a conservative estimate that at least 19% of the exons in the data base were involved in exon shuffling, suggesting an important role for exon shuffling in evolution. Furthermore, these excesses of symmetric exons appear also in those regions of eukaryotic genes that are homologous to prokaryotic genes: the ancient conserved regions. This last fact cannot be explained in terms of the insertional theory of introns but rather supports the concept that some of the introns were ancient, the exon theory of genes.

Nonsense mutation in the phosphofructokinase muscle subunit gene associated with retention of intron 10 in one of the isolated transcripts in Ashkenazi Jewish patients with Tarui disease.

Mutations in the human phosphofructokinase muscle subunit gene (PFKM) are known to cause myopathy classified as glycogenosis type VII (Tarui disease). Previously described molecular defects include base substitutions altering encoded amino acids or resulting in abnormal splicing. We report a mutation resulting in phosphofructokinase deficiency in three patients from an Ashkenazi Jewish family. Using a reverse transcription PCR assay, PFKM subunit transcripts differing by length were detected in skeletal muscle tissue of all three affected subjects. In the longer transcript, an insertion of 252 nucleotides totally homologous to the structure of the 10th intron of the PFKM gene was found separating exon 10 from exon 11. In addition, two single base transitions were identified by direct sequencing: [exon 6; codon 95; CGA (Arg) to TGA (stop)] and [exon 7; codon 172; ACC (Thr) to ACT (Thr)] in either transcript. Single-stranded conformational polymorphism and restriction enzyme analyses confirmed the presence of these point substitutions in genomic DNA and strongly suggested homozygosity for the pathogenic allele. The nonsense mutation at codon 95 appeared solely responsible for the phenotype in these patients, further expanding genetic heterogeneity of Tarui disease. Transcripts with and without intron 10 arising from identical mutant alleles probably resulted from differential pre-mRNA processing and may represent a novel message from the PFKM gene.

Seven newly discovered intron positions in the triose-phosphate isomerase gene: evidence for the introns-late theory.

The gene encoding the glycolytic enzyme triose-phosphate isomerase (TPI; EC 5.3.1.1) has been central to the long-standing controversy on the origin and evolutionary significance of spliceosomal introns by virtue of its pivotal support for the introns-early view, or exon theory of genes. Putative correlations between intron positions and TPI protein structure have led to the conjecture that the gene was assembled by exon shuffling, and five TPI intron positions are old by the criterion of being conserved between animals and plants. We have sequenced TPI genes from three diverse eukaryotes--the basidiomycete Coprinus cinereus, the nematode Caenorhabditis elegans, and the insect Heliothis virescens--and have found introns at seven novel positions that disrupt previously recognized gene/protein structure correlations. The set of 21 TPI introns now known is consistent with a random model of intron insertion. Twelve of the 21 TPI introns appear to be of recent origin since each is present in but a single examined species. These results, together with their implication that as more TPI genes are sequenced more intron positions will be found, render TPI untenable as a paradigm for the introns-early theory and, instead, support the introns-late view that spliceosomal introns have been inserted into preexisting genes during eukaryotic evolution.

A mutated intron sequence codes for an antigenic peptide recognized by cytolytic T lymphocytes on a human melanoma.

We have identified an antigen recognized on a human melanoma by autologous cytolytic T lymphocytes. It is encoded by a gene that is expressed in many normal tissues. Remarkably, the sequence coding for the antigenic peptide is located across an exon-intron junction. A point mutation is present in the intron that generates an amino acid change that is essential for the recognition of the peptide by the anti-tumor cytotoxic T lymphocytes. This observation suggests that the T-cell-mediated surveillance of the integrity of the genome may extend to some intronic regions.

Catalytically critical nucleotide in domain 5 of a group II intron.

Domain 5 (D5) is a small hairpin structure within group II introns. A bimolecular assay system depends on binding by D5 to an intron substrate for self-splicing activity. In this study, mutations in D5 identify two among six nearly invariant nucleotides as being critical for 5' splice junction hydrolysis but unimportant for binding. A mutation at another site in D5 blocks binding. Thus, mutations can distinguish two D5 functions: substrate binding and catalysis. The secondary structure of D5 may resemble helix I formed by the U2 and U6 small nuclear RNAs in the eukaryotic spliceosome. Our results support a revision of the previously proposed correspondence between D5 and helix I on the basis of the critical trinucleotide 5'-AGC-3' present in both. We suggest that this trinucleotide plays a similar role in promoting the chemical reactions for both splicing systems.

Total chemical synthesis of a ribozyme derived from a group I intron.

We describe the complete chemical synthesis of a ribozyme that catalyzes template-directed oligonucleotide ligation. The specific activity of the synthetic ribozyme is nearly identical to that of the same enzyme generated by in vitro transcription with T7 RNA polymerase. The ribozyme is derived from a group I intron and consists of three RNA fragments of 36, 43, and 59 nt that self-assemble to form a catalytically active complex. We have site-specifically substituted ribonucleotide analogs into this enzyme and have identified two 2'-hydroxyl groups that are required for full catalytic activity. In contrast, neither the 2'-hydroxyl nor the exocyclic amino group of the conserved guanosine in the guanosine binding site is necessary for catalysis. By allowing the ribozyme to be modified as easily as its substrates, this synthetic ribozyme system should be useful for testing specific hypotheses concerning ribozyme-substrate interactions and tertiary interactions within the ribozyme.