Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Splicing of constitutive upstream introns is essential for the recognition of intra-exonic suboptimal splice sites in the thrombopoietin gene

The human thrombopoietin (TPO) gene, which codes for the principal cytokine involved in platelet maturation, shows a peculiar alternative splicing of its last exon, where an intra-exonic 116 nt alternative intron is spliced out in a fraction of its mRNA. To characterize the molecular mechanism underlying this alternative splicing, minigenes of TPO genomic constructs with variable exon–intron configurations or carrying exclusively the TPO cDNA were generated and transiently transfected in the Hep3B cell line. We have found that the final rate of the alternative intron splicing is determined by three elements: the presence of upstream constitutive introns, the suboptimal splice sites of the alternative intron and the length of the alternative intron itself. Our results indicate that the recognition of suboptimal intra-exonic splice junctions in the TPO gene is influenced by the assembly of the spliceosome complex on constitutive introns and by a qualitative scanning of the sequence by the transcriptional/splicing machinery complex primed by upstream splicing signals.

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.

The miniaturized nuclear genome of eukaryotic endosymbiont contains genes that overlap, genes that are cotranscribed, and the smallest known spliceosomal introns.

Chlorarachniophyte algae contain a complex, multi-membraned chloroplast derived from the endosymbiosis of a eukaryotic alga. The vestigial nucleus of the endosymbiont, called the nucleomorph, contains only three small linear chromosomes with a haploid genome size of 380 kb and is the smallest known eukaryotic genome. Nucleotide sequence data from a subtelomeric fragment of chromosome III were analyzed as a preliminary investigation of the coding capacity of this vestigial genome. Several housekeeping genes including U6 small nuclear RNA (snRNA), ribosomal proteins S4 and S13, a core protein of the spliceosome [small nuclear ribonucleoprotein (snRNP) E], and a cip-like protease (clpP) were identified. Expression of these genes was confirmed by combinations of Northern blot analysis, in situ hybridization, immunocytochemistry, and cDNA analysis. The protein-encoding genes are typically eukaryotic in overall structure and their messenger RNAs are polyadenylylated. A novel feature is the abundance of 18-, 19-, or 20-nucleotide introns; the smallest spliceosomal introns known. Two of the genes, U6 and S13, overlap while another two genes, snRNP E and clpP, are cotranscribed in a single mRNA. The overall gene organization is extraordinarily compact, making the nucleomorph a unique model for eukaryotic genomics.

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.

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.

Rat growth hormone gene introns stimulate nucleosome alignment in vitro and in transgenic mice.

Average hepatic expression (mRNA per cell per gene) of a metallothionein-rat growth hormone (rGH) gene with its natural introns was about 15-fold higher than an intronless version when tested in transgenic mice. We examined the idea that intron removal leads to an alteration in chromatin structure that might be responsible for this effect. Using an in vitro chromatin assembly system, we observed that nucleosomes were aligned in a characteristic ordered array over the gene and promoter when all introns were present. Linker histones were necessary for this alignment to occur. In contrast, nucleosome alignment was perturbed in constructs lacking some or all of the introns. A similar disruption of nucleosome alignment was observed when comparing chromatin from livers of transgenic mice carrying rGH transgenes with or without introns. In vitro, sequences at the 3' end of the rGH gene position nucleosomes and facilitate nucleosome alignment upstream; however, nucleosome alignment does not occur on the approximately 3 kb of downstream flanking rat sequence. These observations suggest that signals present in genomic rGH DNA may serve to establish appropriate nucleosome alignment during development and, possibly, to restore nucleosome alignment to the transcribed region after disruption incurred by the passage of an RNA polymerase molecule, thereby facilitating subsequent rounds of transcription.

Rpp2, an essential protein subunit of nuclear RNase P, is required for processing of precursor tRNAs and 35S precursor rRNA in Saccharomyces cerevisiae

RPP2, an essential gene that encodes a 15.8-kDa protein subunit of nuclear RNase P, has been identified in the genome of Saccharomyces cerevisiae. Rpp2 was detected by sequence similarity with a human protein, Rpp20, which copurifies with human RNase P. Epitope-tagged Rpp2 can be found in association with both RNase P and RNase mitochondrial RNA processing in immunoprecipitates from crude extracts of cells. Depletion of Rpp2 protein in vivo causes accumulation of precursor tRNAs with unprocessed introns and 5′ and 3′ termini, and leads to defects in the processing of the 35S precursor rRNA. Rpp2-depleted cells are defective in processing of the 5.8S rRNA. Rpp2 immunoprecipitates cleave both yeast precursor tRNAs and precursor rRNAs accurately at the expected sites and contain the Rpp1 protein orthologue of the human scleroderma autoimmune antigen, Rpp30. These results demonstrate that Rpp2 is a protein subunit of nuclear RNase P that is functionally conserved in eukaryotes from yeast to humans.

A protein encoded by a group I intron in Aspergillus nidulans directly assists RNA splicing and is a DNA endonuclease

Some group I introns self-splice in vitro, but almost all are thought to be assisted by proteins in vivo. Mutational analysis has shown that the splicing of certain group I introns depends upon a maturase protein encoded by the intron itself. However the effect of a protein on splicing can be indirect. We now provide evidence that a mitochondrial intron-encoded protein from Aspergillus nidulans directly facilitates splicing in vitro. This demonstrates that a maturase is an RNA splicing protein. The protein-assisted reaction is as fast as that of any other known group I intron. Interestingly the protein is also a DNA endonuclease, an activity required for intron mobilization. Mobile elements frequently encode proteins that promote their propagation. Intron-encoded proteins that also assist RNA splicing would facilitate both the transposition and horizontal transmission of introns.

The product of ORF O located within the domain of herpes simplex virus 1 genome transcribed during latent infection binds to and inhibits in vitro binding of infected cell protein 4 to its cognate DNA site

The partially overlapping ORF P and ORF O are located within the domains of the herpes simplex virus 1 genome transcribed during latency. Earlier studies have shown that ORF P is repressed by infected cell protein 4 (ICP4), the major viral regulatory protein, binding to its cognate site at the transcription initiation site of ORF P. The ORF P protein binds to p32, a component of the ASF/SF2 alternate splicing factors; in cells infected with a recombinant virus in which ORF P was derepressed there was a significant decrease in the expression of products of key regulatory genes containing introns. We report that (i) the expression of ORF O is repressed during productive infection by the same mechanism as that determining the expression of ORF P; (ii) in cells infected at the nonpermissive temperature for ICP4, ORF O protein is made in significantly lower amounts than the ORF P protein; (iii) the results of insertion of a sequence encoding 20 amino acids between the putative initiator methionine codons of ORF O and ORF P suggest that ORF O initiates at the methionine codon of ORF P and that the synthesis of ORF O results from frameshift or editing of its RNA; and (iv) glutathione S-transferase–ORF O fusion protein bound specifically ICP4 and precluded its binding to its cognate site on DNA in vitro. These and earlier results indicate that ORF P and ORF O together have the capacity to reduce the synthesis or block the expression of regulatory proteins essential for viral replication in productive infection.

Intron “sliding” and the diversity of intron positions

Alignments of homologous genes typically reveal a great diversity of intron locations, far more than could fit comfortably in a single gene. Thus, a minority of these intron positions could be inherited from a single ancestral gene, but the larger share must be attributed to subsequent events of intron gain or intron “sliding” (movement from one position to another within a gene). Intron sliding has been argued from cases of discordant introns and from putative spatial clustering of intron positions. A list of 32 cases of discordant introns is presented here. Most of these cases are found to be artefactual. The spatial and phylogenetic distributions of intron positions from five published compilations of gene data, comprising 205 intron positions, have been examined systematically for evidence of intron sliding. The results suggest that sliding, if it occurs at all, has contributed little to the diversity of intron positions.

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.

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.