Ovomucoid intervening sequences specify functional domains and generate protein polymorphism

Two thirds of the natural chicken ovomucoid gene has been sequenced, including all exons and the intron sequences surrounding all fourteen intron/ exon junctions. The junction sequences surrounding four of the introns are redundant; however, the sequences surrounding the other three introns contain no redundancies and thus the splicing sites at either end of these three introns are unambiguous. The splicing in all cases conforms to the GT-AG rule. The ovomucoid gene sequence around intron F can be used to predict the cause of an internal deletion polymorphism in the ovomucoid protein, which is an apparent error in the processing of the ovomucoid pre-mRNA. We also compare the structural organization of the ovomucoid gene with the ovomucoid protein sequence to examine theories of the evolution of ovomucoids as well as the origin of intervening sequences. This analysis suggests that the present ovomucoid gene evolved from a primordial ovomucoid gene by two separate intragenic duplications. Furthermore, sequence analyses suggest that introns were present in the primordial ovomucoid gene before birds and mammals diverged, about 300 million years ago. Finally, the positions of the introns within the ovomucoid gene support the theory that introns separate gene segments that code for functional domains of proteins and provide insight on the manner by which eucaryotic genes were constructed during the process of evolution.     

Complete genomic structure of human rpS3: identification of functional U15b snoRNA in the fifth intron

Analysis of the complete genomic structure of the human ribosomal protein S3 (rpS3) gene revealed the presence of a functional U15b snoRNA gene in its intron. Human ribosomal protein S3 (rpS3) gene of 6115 bp long has been identified to contain six introns and seven exons in this study. The first and fifth introns of human S3 gene contain functional U15 snoRNA genes. Although Xenopus and Fugu counterparts also have six introns and seven exons, S3 gene of Fugu contains two functional U15 snoRNAs in the fourth and sixth introns and two pseudo genes for U15 snoRNAs in the first and fifth introns. In Xenopus S1 gene encoding ribosomal protein S3, however, three of its six introns contain U15 snoRNA gene sequence. Sequence comparison of the U15 genes from Xenopus, Fugu and human revealed that the regions involved in binding to 28S rRNA and the consensus sequence (C, D and D' boxes) for snoRNAs are highly conserved among those genes from these three species. Human U15a and U15b RNAs which are derived from the first and the fifth introns, respectively, have been identified to be functional by microinjection of human U15a and U15b snoRNAs into Xenopus oocyte. Northern blot and primer extension analyses confirm that human U15b snoRNA is expressed in vivo.     

Computer identification of snoRNA genes using a Mammalian Orthologous Intron Database

Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.     

Evolution of collagen IV genes from a 54-base pair exon: A role for introns in gene evolution

The exon structure of the collagen IV gene provides a striking example for collagen evolution and the role of introns in gene evolution. Collagen IV, a major component of basement membranes, differs from the fibrillar collagens in that it contains numerous interruptions in the triple helical Gly-X-Y repeat domain. We have characterized all 47 exons in the mouse alpha 2(IV) collagen gene and find two 36-, two 45-, and one 54-bp exons as well as one 99- and three 108-bp exons encoding the Gly-X-Y repeat sequence. All these exons sizes are also found in the fibrillar collagen genes. Strikingly, of the 24 interruption sequences present in the alpha 2-chain of mouse collagen IV, 11 are encoded at the exon/intron borders of the gene, part of one interruption sequence is encoded by an exon of its own, and the remaining interruptions are encoded within the body of exons. In such "fusion exons" the Gly-X-Y encoding domain is also derived from 36-, 45-, or 54-bp sequence elements. These data support the idea that collagen IV genes evolved from a primordial 54-bp coding unit. We furthermore interpret these data to suggest that the interruption sequences in collagen IV may have evolved from introns, presumably by inactivation of splice site signals, following which intronic sequences could have been recruited into exons. We speculated that this mechanism could provide a role for introns in gene evolution in general.     

Structure and organization of the human transglutaminase 1 gene.

Membrane-associated transglutaminases (TGase1) have recently been found to be common in mammalian cells, but it is not clear whether these derive from the same or different genes. In order to determine the complexity of this system, we have isolated and characterized the human gene (TGM1). The gene of 14,133 base pairs was found to contain 15 exons spliced by 14 introns. Interestingly, the positions of these introns have been conserved in comparison with the genes of two other transglutaminase-like activities described in the literature, but the TGM1 gene is by far the smallest characterized to date because its introns are relatively smaller. On the other hand, the TGase1 enzyme is the largest known transglutaminase (about 90 kDa), apparently because its gene acquired tracts that encode additional sequences on its amino and carboxyl termini that confer its unique properties. Southern blot analyses of total human genomic DNA cut with several restriction enzymes reveal only one band. Use of human-rodent cell hybrid panels and chromosomal in situ hybridization with biotin-labeled probes revealed that the human TGM1 gene maps to chromosome position 14q11.2-13. Such data suggest there is a single gene copy per haploid human genome. Comparisons of sequence identities and homologies indicate that the transglutaminase family of genes arose by duplications and subsequent divergent evolution from a common ancestor but later became scattered in the human genome. Although our present Southern blot and chromosomal localization studies revealed no restriction fragment length polymorphisms, comparisons of published sequences and our genomic clone indicate there are two sequence variants for TGase1 within the human population. The rare smaller variant contains a two-nucleotide deletion near the 5'-end, uses an alternate initiation codon, and differs from the common larger variant only in the first 15 amino acids. Furthermore, the DNA sequences of intron 14 possess several tracts of dinucleotide repeats that by polymerase chain reaction analysis show wide size polymorphism within the human population. Accordingly, this gene system constitutes a useful polymorphic marker for genetic linkage analyses.     

Incipient mitochondrial evolution in yeasts. II. The complete sequence of the gene coding for cytochrome b in Saccharomyces douglasii reveals the presence of both new and conserved introns and discloses major differences in the fixation of mutations in evolution.

We have determined the complete sequence of the mitochondrial gene coding for cytochrome b in Saccharomyces douglasii. The gene is 6310 base-pairs long and is interrupted by four introns. The first one (1311 base-pairs) belongs to the group ID of secondary structure, contains a fragment open reading frame with a characteristic GIY ... YIG motif, is absent from Saccharomyces cerevisiae and is inserted in the same site in which introns 1 and 2 are inserted in Neurospora crassa and Podospora anserina, respectively. The next three S. douglasii introns are homologous to the first three introns of S. cerevisiae, are inserted at the same positions and display various degrees of similarity ranging from an almost complete identity (intron 2 and 4) to a moderate one (intron 3). We have compared secondary structures of intron RNAs, and nucleotide and amino acid sequences of cytochrome b exons and intron open reading frames in the two Saccharomyces species. The rules that govern fixation of mutations in exon and intron open reading frames are different: the relative proportion of mutations occurring in synonymous codons is low in some introns and high in exons. The overall frequency of mutations in cytochrome b exons is much smaller than in nuclear genes of yeasts, contrary to what has been found in vertebrates, where mitochondrial mutations are more frequent. The divergence of the cytochrome b gene is modular: various parts of the gene have changed with a different mode and tempo of evolution.     

Genetic organization of alpha-bungarotoxins from Bungarus multicinctus (Taiwan banded krait): evidence showing that the production of alpha-bungarotoxin isotoxins is not derived from edited mRNAs.

Two genomic DNAs with a size of approximately 2.8 kb, isolated from the liver of Bungarus multicinctus (Taiwan banded krait), encode the precursors of the long neurotoxins, alpha-Bgt(A31) and alpha-Bgt(V31), respectively. Both genes share virtually identical overall organization with three exons separated by two introns, which were inserted in the same positions in the coding regions of the genes. Moreover, their nucleotide sequences share approximately 98% identity. This result indicates that the two genes co-exist in the genome of B.multicinctus, and probably arose from gene duplication. The exon/intron structures of the alpha-Bgt genes were essentially the same as those reported for the short neurotoxins. This reflects that the long and short neurotoxins should share a common evolutionary origin. Comparative analyses on long neurotoxin and short neurotoxin genes showed that the protein coding regions of the exons were more diverse than the introns except for the signal peptide domain. This implies that the protein coding regions of the neurotoxins may have evolved via accelerated evolution. PCR amplification of venom gland cDNA mixtures revealed that only two amino acid sequences corresponding to alpha-Bgt(A31) and alpha-Bgt(V31) could be deduced from the cDNAs. The results of chromatographic analyses and protein sequencing again emphasized the view that, with the exception of alpha-Bgt(A31) and alpha-Bgt(V31), no other alpha-Bgt isotoxins with amino acid substitutions were present in B.multicinctus venom. In contrast to the proposition of Liu et al. ( Nucleic Acids Res., 1998,26, 5624-5629), our findings strongly suggest that each alpha-Bgt isotoxin is derived from the respective gene, and that alpha-Bgt RNA polymorphism does not originate from one single, intronless gene by the mechanism of RNA editing.     

Evolutionary implications of intron–exon distribution and the properties and sequences of the RPL10A gene in eukaryotes

The RPL10A gene encodes the RPL10 protein, required for joining 40S and 60S subunits into a functional 80S ribosome. This highly conserved gene, ubiquitous across all eukaryotic super-groups, is characterized by a variable number of spliceosomal introns, present in most organisms. These properties facilitate the recognition of orthologs among distant taxa and thus comparative studies of sequences as well as the distribution and properties of introns in taxonomically distant groups of eukaryotes. The present study examined the multiple ways in which RPL10A conservation vs. sequence changes in the gene over the course of evolution, including in exons, introns, and the encoded proteins, can be exploited for evolutionary analysis at different taxonomic levels. At least 25 different positions harboring introns within the RPL10A gene were determined in different taxa, including animals, plants, fungi, and alveolates. Generally, intron positions were found to be well conserved even across different kingdoms. However, certain introns seemed to be restricted to specific groups of organisms. Analyses of several properties of introns, including insertion site, phase, and length, along with exon and intron GC content and exon–intron boundaries, suggested biases within different groups of organisms. The use of a standard primer pair to analyze a portion of the intron-containing RPL10A gene in 12 genera of green algae within Chlorophyta is presented as a case study for evolutionary analyses of introns at intermediate and low taxonomic levels. Our study shows that phylogenetic reconstructions at different depths can be achieved using RPL10A nucleotide sequences from both exons and introns as well as the amino acid sequences of the encoded protein.     

Characterization of rbcL group IA introns from two colonial volvocalean species (Chlorophyceae)

Group I introns were reported for the first time in the large subunit of Rubisco (rbcL) genes, using two colonial green algae, Pleodorina californica and Gonium multicoccum (Volvocales). The rbcL gene of P. californica contained an intron (PlC intron) of 1320 bp harboring an open reading frame (ORF). The G. multicoccum rbcL gene had two ORF-lacking introns of 549 (GM1 intron) and 295 (GM2 intron) base pairs. Based on the conserved nucleotide sequences of the secondary structure, the PlC and GM1 introns were assigned to group IA2 whereas the GM2 intron belonged to group IA1. Southern hybridization analyses of nuclear and chloroplast DNAs indicated that such intron-containing rbcL genes are located in the chloroplast genome. Sequencing RNAs from the two algae revealed that these introns are spliced out during mRNA maturation. In addition, the PlC and GM1 introns were inserted in the same position of the rbcL exons, and phylogenetic analysis of group IA introns indicated a close phylogenetic relationship between the PlC and GM1 introns within the lineage of bacteriophage group IA2 introns. However, P. californica and G. multicoccum occupy distinct clades in the phylogenetic trees of the colonial Volvocales, and the majority of other colonial volvocalean species do not have such introns in the rbcL genes. Therefore, these introns might have been recently inserted in the rbcL genes independently by horizontal transmission by viruses or bacteriophage.     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

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.     

Intron Dynamics and the Evolution of Integrin β-Subunit Genes: Maintenance of an Ancestral Gene Structure in the Coral, Acropora millepora

We have determined the genomic structure of an integrin β-subunit gene from the coral, Acropora millepora. The coding region of the gene contains 26 introns, spaced relatively uniformly, and this is significantly more than have been found in any integrin β-subunit genes from higher animals. Twenty-five of the 26 coral introns are also found in a β-subunit gene from at least one other phylum, indicating that the coral introns are ancestral. While there are some suggestions of intron gain or sliding, the predominant theme seen in the homologues from higher animals is extensive intron loss. The coral baseline allows one to infer that a number of introns found in only one phylum of higher animals result from frequent intron loss, as opposed to the seemingly more parsimonious alternative of isolated intron gain. The patterns of intron loss confirm results from protein sequences that most of the vertebrate genes, with the exception of β4, belong to one of two β subunit families. The similarity of the patterns within each of the β1,2,7 and β3,5,6,8 groups indicates that these gene structures have been very stable since early vertebrate evolution. Intron loss has been more extensive in the invertebrate genes, and obvious patterns have yet to emerge in this more limited data set. Received: 5 March 2001 / Accepted: 17 May 2001     

Intron size, abundance, and distribution within untranslated regions of genes

Most research concerning the evolution of introns has largely considered introns within coding sequences (CDSs), without regard for introns located within untranslated regions (UTRs) of genes. Here, we directly determined intron size, abundance, and distribution in UTRs of genes using full-length cDNA libraries and complete genome sequences for four species, Arabidopsis thaliana, Drosophila melanogaster, human, and mouse. Overall intron occupancy (introns/exon kbp) is lower in 5' UTRs than CDSs, but intron density (intron occupancy in regions containing introns) tends to be higher in 5' UTRs than in CDSs. Introns in 5' UTRs are roughly twice as large as introns in CDSs, and there is a sharp drop in intron size at the 5' UTR-CDS boundary. We propose a mechanistic explanation for the existence of selection for larger intron size in 5' UTRs, and outline several implications of this hypothesis. We found introns to be randomly distributed within 5' UTRs, so long as a minimum required exon size was assumed. Introns in 3' UTRs were much less abundant than in 5' UTRs. Though this was expected for human and mouse that have intron-dependent nonsense-mediated decay (NMD) pathways that discourage the presence of introns within the 3' UTR, it was also true for A. thaliana and D. melanogaster, which may lack intron-dependent NMD. Our findings have several implications for theories of intron evolution and genome evolution in general.     

Comparative analysis of the structural organization of two closely related vitellogenin genes in X. laevis

The structural organization of the two closely related vitellogenin genes A1 and A2 has been determined and compared by electron microscopy. In both genes the mRNA-coding sequence of 6 kb is interrupted 33 times, leading to a total gene length of 21 kb for gene A1 and 16 kb for gene A2. Thus both genes have a mean exon length of 0.175 kb, while the mean intron length is 0.45 kb in gene A1 and 0.31 kb in gene A2. Because the introns interrupt the structural sequence at homologous positions in genes A1 and A2, we suggest that these two genes are the products of a duplication of an ancestral gene which had an intron-exon arrangement similar to that of the extant genes. Since the duplication event, the sequence and length of the analogous introns have changed rapidly, whereas homologous exons have diverged to an extent of only 5% of their sequences. The results suggest different mechanisms of evolution for exons and introns. While the exons evolved primarily by point mutations, such mutations, as well as deletion, insertion and duplication events, were important in the evolution of the introns.     

Different expression strategy: multiple intronic gene clusters of box H/ACA snoRNA in Drosophila melanogaster

The high degree of rRNA pseudouridylation in Drosophila melanogaster provides a good model for studying the genomic organization, structural and functional diversity of box H/ACA small nucleolar RNAs (snoRNAs). Accounting for both conserved sequence motifs and secondary structures, we have developed a computer-assisted method for box H/ACA snoRNA searching. Ten snoRNA clusters containing 42 box H/ACA snoRNAs were identified from D.melanogaster. Strikingly, they are located in the introns of eight protein-coding genes. In contrast to the mode of one snoRNA per intron so far observed in all animals, our results demonstrate for the first time a novel polycistronic organization that implies a different expression strategy for a box H/ACA snoRNA gene when compared to box C/D snoRNAs in D.melanogaster. Mutiple isoforms of the box H/ACA snoRNAs, from which most clusters are made up, were observed in D.melanogaster. The degree of sequence similarity between the isoforms varies from 99% to 70%, implying duplication events in different periods and a trend of enlarging the intronic snoRNA clusters. The variation in the functional elements of the isoforms could lead to partial alternation of snoRNA's function in loss or gain of rRNA complementary sequences and probably contributes to the great diversity of rRNA pseudouridylation in D.melanogaster.