Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

Lactone-ring-cleaving enzyme: Genetic analysis, novel RNA editing, and evolutionary implications

A lactonohydrolase from Fusarium oxysporum AKU 3702 is an enzyme catalyzing the hydrolysis of aldonate lactones to the corresponding aldonic acids. The amino acid sequences of the NH2 terminus and internal peptide fragments of the enzyme were determined to prepare synthetic oligonucleotides as primers for the PCR. An approximate 1,000-base genomic DNA fragment thus amplified was used as the probe to clone both genomic DNA and cDNA for the enzyme. The lactonohydrolase genomic gene consists of six exons separated by five short introns. A novel type of RNA editing, in which lactonohydrolase mRNA included the insertion of guanosine and cytidine residues, was observed. The predicted amino acid sequence of the cloned lactonohydrolase cDNA showed significant similarity to those of the gluconolactonase from Zymomonas mobilis, and paraoxonases from human and rabbit, forming a unique superfamily consisting of C-O cleaving enzymes and P-O cleaving enzymes. Lactonohydrolase was expressed under the control of the lac promoter in Escherichia coli.

Cloning of mDEAH9, a putative RNA helicase and mammalian homologue of Saccharomyces cerevisiae splicing factor Prp43

Yeast splicing factor Prp43, a DEAH box protein of the putative RNA helicase/RNA-dependent NTPase family, is a splicing factor that functions late in the pre-mRNA splicing pathway to facilitate spliceosome disassembly. In this paper we report cDNA cloning and characterization of mDEAH9, an apparent mammalian homologue of Prp43. Amino acid sequence comparison revealed that the two proteins are ≈65% identical over a 500-aa region spanning the central helicase domain and the C-terminal region. Expression of mDEAH9 in S. cerevisiae bearing a temperature-sensitive mutation in prp43 was sufficient to restore growth at the nonpermissive temperature. This functional complementation was specific, as mouse mDEAH9 failed to complement mutations in related splicing factor genes prp16 or prp22. Finally, double label immunofluorescence experiments performed with mammalian cells revealed colocalization of mDEAH9 and splicing factor SC35 in punctate nuclear speckles. Thus, the hypothesis that mDEAH9 represents the mammalian homologue of yeast Prp43 is supported by its high sequence homology, functional complementation, and colocalization with a known splicing factor in the nucleus. Our results provide additional support for the hypothesis that the spliceosomal machinery that mediates regulated, dynamic changes in conformation of pre-mRNA and snRNP RNAs has been highly conserved through evolution.

Suspensor-derived polyembryony caused by altered expression of valyl-tRNA synthetase in the twn2 mutant of Arabidopsis

The twn2 mutant of Arabidopsis exhibits a defect in early embryogenesis where, following one or two divisions of the zygote, the decendents of the apical cell arrest. The basal cells that normally give rise to the suspensor proliferate abnormally, giving rise to multiple embryos. A high proportion of the seeds fail to develop viable embryos, and those that do, contain a high proportion of partially or completely duplicated embryos. The adult plants are smaller and less vigorous than the wild type and have a severely stunted root. The twn2-1 mutation, which is the only known allele, was caused by a T-DNA insertion in the 5′ untranslated region of a putative valyl-tRNA synthetase gene, valRS. The insertion causes reduced transcription of the valRS gene in reproductive tissues and developing seeds but increased expression in leaves. Analysis of transcript initiation sites and the expression of promoter–reporter fusions in transgenic plants indicated that enhancer elements inside the first two introns interact with the border of the T-DNA to cause the altered pattern of expression of the valRS gene in the twn2 mutant. The phenotypic consequences of this unique mutation are interpreted in the context of a model, suggested by Vernon and Meinke [Vernon, D. M. & Meinke, D. W. (1994) Dev. Biol. 165, 566–573], in which the apical cell and its decendents normally suppress the embryogenic potential of the basal cell and its decendents during early embryo development.

A functional homolog of a yeast tRNA splicing enzyme is conserved in higher eukaryotes and in Escherichia coli

tRNA splicing in the yeast Saccharomyces cerevisiae requires an endonuclease to excise the intron, tRNA ligase to join the tRNA half-molecules, and 2′-phosphotransferase to transfer the splice junction 2′-phosphate from ligated tRNA to NAD, producing ADP ribose 1′′–2′′ cyclic phosphate (Appr>p). We show here that functional 2′-phosphotransferases are found throughout eukaryotes, occurring in two widely divergent yeasts (Candida albicans and Schizosaccharomyces pombe), a plant (Arabidopsis thaliana), and mammals (Mus musculus); this finding is consistent with a role for the enzyme, acting in concert with ligase, to splice tRNA or other RNA molecules. Surprisingly, functional 2′-phosphotransferase is found also in the bacterium Escherichia coli, which does not have any known introns of this class, and does not appear to have a ligase that generates junctions with a 2′-phosphate. Analysis of the database shows that likely members of the 2′-phosphotransferase family are found also in one other bacterium (Pseudomonas aeruginosa) and two archaeal species (Archaeoglobus fulgidus and Pyrococcus horikoshii). Phylogenetic analysis reveals no evidence for recent horizontal transfer of the 2′-phosphotransferase into Eubacteria, suggesting that the 2′-phosphotransferase has been present there since close to the time that the three kingdoms diverged. Although 2′-phosphotransferase is not present in all Eubacteria, and a gene disruption experiment demonstrates that the protein is not essential in E. coli, the continued presence of 2′-phosphotransferase in Eubacteria over large evolutionary times argues for an important role for the protein.

Explosive invasion of plant mitochondria by a group I intron

Group I introns are mobile, self-splicing genetic elements found principally in organellar genomes and nuclear rRNA genes. The only group I intron known from mitochondrial genomes of vascular plants is located in the cox1 gene of Peperomia, where it is thought to have been recently acquired by lateral transfer from a fungal donor. Southern-blot surveys of 335 diverse genera of land plants now show that this intron is in fact widespread among angiosperm cox1 genes, but with an exceptionally patchy phylogenetic distribution. Four lines of evidence—the intron’s highly disjunct distribution, many incongruencies between intron and organismal phylogenies, and two sources of evidence from exonic coconversion tracts—lead us to conclude that the 48 angiosperm genera found to contain this cox1 intron acquired it by 32 separate horizontal transfer events. Extrapolating to the over 13,500 genera of angiosperms, we estimate that this intron has invaded cox1 genes by cross-species horizontal transfer over 1,000 times during angiosperm evolution. This massive wave of lateral transfers is of entirely recent occurrence, perhaps triggered by some key shift in the intron’s invasiveness within angiosperms.

The role of alternative mRNA splicing in generating heterogeneity within the Anopheles gambiae class I glutathione S-transferase family

The class I glutathione S-transferases (GSTs) of Anopheles gambiae are encoded by a complex gene family. We describe the genomic organization of three members of this family, which are sequentially arranged on the chromosome in divergent orientations. One of these genes, aggst1-2, is intronless and has been described. In contrast, the two A. gambiae GST genes (aggst1α and aggst1β) reported within are interrupted by introns. The gene aggst1α contains five coding exons that are alternatively spliced to produce four mature GST transcripts, each of which contains a common 5′ exon encoding the N termini of the GST protein spliced to one of four distinct 3′ exons encoding the carboxyl termini. All four of the alternative transcripts of aggst1α are expressed in A. gambiae larvae, pupae, and adults. We report on the involvement of alternative RNA splicing in generating multiple functional GST transcripts. A cDNA from the aggst1β gene was detected in adult mosquitoes, demonstrating that this GST gene is actively transcribed. The percentage similarity of the six cDNAs transcribed from the three GST genes range from 49.5% to 83.1% at the nucleotide level.

Organization of the human myostatin gene and expression in healthy men and HIV-infected men with muscle wasting

Myostatin, a member of the transforming growth factor-β superfamily, is a genetic determinant of skeletal muscle growth. Mice and cattle with inactivating mutations of myostatin have marked muscle hypertrophy. However, it is not known whether myostatin regulates skeletal muscle growth in adult men and whether increased myostatin expression contributes to wasting in chronic illness. We examined the hypothesis that myostatin expression correlates inversely with fat-free mass in humans and that increased expression of the myostatin gene is associated with weight loss in men with AIDS wasting syndrome. We therefore cloned the human myostatin gene and cDNA and examined the gene’s expression in the skeletal muscle and serum of healthy and HIV-infected men. The myostatin gene comprises three exons and two introns, maps to chromosomal region 2q33.2, has three putative transcription initiation sites, and is transcribed as a 3.1-kb mRNA species that encodes a 375-aa precursor protein. Myostatin is expressed uniquely in the human skeletal muscle as a 26-kDa mature glycoprotein (myostatin-immunoreactive protein) and secreted into the plasma. Myostatin immunoreactivity is detectable in human skeletal muscle in both type 1 and 2 fibers. The serum and intramuscular concentrations of myostatin-immunoreactive protein are increased in HIV-infected men with weight loss compared with healthy men and correlate inversely with fat-free mass index. These data support the hypothesis that myostatin is an attenuator of skeletal muscle growth in adult men and contributes to muscle wasting in HIV-infected men.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5′ flanking sequence (−855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5′ flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.

RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs

On the basis of the sequence of the mitochondrial genome in the flowering plant Arabidopsis thaliana, RNA editing events were systematically investigated in the respective RNA population. A total of 456 C to U, but no U to C, conversions were identified exclusively in mRNAs, 441 in ORFs, 8 in introns, and 7 in leader and trailer sequences. No RNA editing was seen in any of the rRNAs or in several tRNAs investigated for potential mismatch corrections. RNA editing affects individual coding regions with frequencies varying between 0 and 18.9% of the codons. The predominance of RNA editing events in the first two codon positions is not related to translational decoding, because it is not correlated with codon usage. As a general effect, RNA editing increases the hydrophobicity of the coded mitochondrial proteins. Concerning the selection of RNA editing sites, little significant nucleotide preference is observed in their vicinity in comparison to unedited C residues. This sequence bias is, per se, not sufficient to specify individual C nucleotides in the total RNA population in Arabidopsis mitochondria.

Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

Development of an epithelium-specific expression cassette with human DNA regulatory elements for transgene expression in lung airways

The efficient expression of therapeutic genes in target cells or tissues is an important component of efficient and safe gene therapy. Utilizing regulatory elements from the human cytokeratin 18 (K18) gene, including 5′ genomic sequences and one of its introns, we have developed a novel expression cassette that can efficiently express reporter genes, as well as the human cystic fibrosis transmembrane conductance regulator (CFTR) gene, in cultured lung epithelial cells. CFTR transcripts expressed from the native K18 enhancer/promoter include two alternative splicing products, due to the activation of two cryptic splice sites in the CFTR coding region. Modification of the K18 intron and CFTR cDNA sequences eliminated the cryptic splice sites without changing the CFTR amino acid sequence, and led to enhanced CFTR mRNA and protein expression as well as biological function. Transgenic expression analysis in mice showed that the modified expression cassette can direct efficient and epithelium-specific expression of the Escherichia coli LacZ gene in the airways of fetal lungs, with no detectable expression in lung fibroblasts or endothelial cells. This is the first expression cassette which selectively directs lung transgene expression for CFTR gene therapy to airway epithelia.

tRNA Nuclear Export in Saccharomyces cerevisiae: In Situ Hybridization Analysis

To understand the factors specifically affecting tRNA nuclear export, we adapted in situ hybridization procedures to locate endogenous levels of individual tRNA families in wild-type and mutant yeast cells. Our studies of tRNAs encoded by genes lacking introns show that nucleoporin Nup116p affects both poly(A) RNA and tRNA export, whereas Nup159p affects only poly(A) RNA export. Los1p is similar to exportin-t, which facilitates vertebrate tRNA export. A los1 deletion mutation affects tRNA but not poly(A) RNA export. The data support the notion that Los1p and exportin-t are functional homologues. Because LOS1 is nonessential, tRNA export in vertebrate and yeast cells likely involves factors in addition to exportin-t. Mutation of RNA1, which encodes RanGAP, causes nuclear accumulation of tRNAs and poly(A) RNA. Many yeast mutants, including those with the rna1-1 mutation, affect both pre-tRNA splicing and RNA export. Our studies of the location of intron-containing pre-tRNAs in the rna1-1 mutant rule out the possibility that this results from tRNA export occurring before splicing. Our results also argue against inappropriate subnuclear compartmentalization causing defects in pre-tRNA splicing. Rather, the data support “feedback” of nucleus/cytosol exchange to the pre-tRNA splicing machinery.

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.”

Architectural limits on split genes

Exon/intron architecture varies across the eukaryotic kingdom with large introns and small exons the rule in vertebrates and the opposite in lower eukaryotes. To investigate the relationship between exon and intron size in pre-mRNA processing, internally expanded exons were placed in vertebrate genes with small and large introns. Both exon and intron size influenced splicing phenotype. Intron size dictated if large exons were efficiently recognized. When introns were large, large exons were skipped; when introns were small, the same large exons were included. Thus, large exons were incompatible for splicing if and only if they were flanked by large introns. Both intron and exon size became problematic at ≈500 nt, although both exon and intron sequence influenced the size at which exons and introns failed to be recognized. These results indicate that present-day gene architecture reflects at least in part limitations on exon recognition. Furthermore, these results strengthen models that invoke pairing of splice sites during recognition of pre-mRNAs, and suggest that vertebrate consensus sequences support pairing across either introns or exons.