Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Investigating the intron recognition mechanism in eukaryotes

Recent studies indicate that many introns, as well as the complex spliceosomal mechanism to remove them, were present early in eukaryotic evolution. This study examines intron and exon characteristics from annotations of whole genomes to investigate the intron recognition mechanism. Exon definition uses the exon as the unit of recognition, placing length constraints on the exon but not on the intron (allowing it a greater range of lengths). In contrast, intron definition uses the intron itself as the unit of recognition and thus removes constraints on internal exon length forced by the use of an exon definition mechanism. Thus, intron and exon lengths within a genome can reflect the constraints imposed by its splicing. This study shows that it is possible firstly to recover valid intron and exon information from genome annotation. We then compare internal intron and exon information from a range of eukaryotic genomes and investigate possible evolutionary length constraints on introns and exons and how they can impact on the intron recognition mechanism. Results indicate that exon definition-based mechanisms may predominate in vertebrates although the exact system in fish is expected to show some differences with the better characterized system from mammals. We also raise the possibility that the last common ancestor of plants and animals contained some type of exon definition and that this mechanism was replaced in some genes and lineages by intron definition, possibly as a result of intron loss and/or intron shortening.     

Intron-exon structures of eukaryotic model organisms.

To investigate the distribution of intron-exon structures of eukaryotic genes, we have constructed a general exon database comprising all available intron-containing genes and exon databases from 10 eukaryotic model organisms: Homo sapiens, Mus musculus, Gallus gallus, Rattus norvegicus, Arabidopsis thaliana, Zea mays, Schizosaccharomyces pombe, Aspergillus, Caenorhabditis elegans and Drosophila. We purged redundant genes to avoid the possible bias brought about by redundancy in the databases. After discarding those questionable introns that do not contain correct splice sites, the final database contained 17 102 introns, 21 019 exons and 2903 independent or quasi-independent genes. On average, a eukaryotic gene contains 3.7 introns per kb protein coding region. The exon distribution peaks around 30-40 residues and most introns are 40-125 nt long. The variable intron-exon structures of the 10 model organisms reveal two interesting statistical phenomena, which cast light on some previous speculations. (i) Genome size seems to be correlated with total intron length per gene. For example, invertebrate introns are smaller than those of human genes, while yeast introns are shorter than invertebrate introns. However, this correlation is weak, suggesting that other factors besides genome size may also affect intron size. (ii) Introns smaller than 50 nt are significantly less frequent than longer introns, possibly resulting from a minimum intron size requirement for intron splicing.     

Statistical analysis of the exon-intron structure of higher eukaryote genes]

The exon-intron structure of human, insect (Drosophila sp.), and dicot plant (Arabidopsis thaliana) genes was considered. In each genome there exists a characteristic intron length. Anomalously long introns was usually the first introns in genes. In each sample there are correlations between the lengths of neighboring exons and between exon lengths and closeness to the consensus of the sites at exon boundaries. Exons and exon pairs containing an integer number of triplets are preferred. These results are relevant to the study of splicing mechanism and evolution of introns, as well as construction of gene recognition algorithms.     

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.     

Differential GC content between exons and introns establishes distinct strategies of splice-site recognition

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition.     

Changes of introns and exons length in genes of arabidopsis, rice, nematode and human

Genes of Arabidopsis thaliana, Orysa sativa, Caenorhabditis elegans, Homo sapiens have been studied by computer analysis. The average intron and exon lengths in genes of these organisms decreases with increase of intron number in genes. The length of introns and exons in A. thaliana and O. sativa genes is change with increase of intron number in genes by high coefficient of correlation. Linear dependence between the sum of exon lengths and intron number in genes increased proportionally to number of gene introns. The average length of introns and genes of human depend on density of genes in DNA.     

Evidence of splice signal migration from exon to intron during intron evolution

A comparison of the nucleotide sequences around the splice junctions that flank old (shared by two or more major lineages of eukaryotes) and new (lineage-specific) introns in eukaryotic genes reveals substantial differences in the distribution of information between introns and exons. Old introns have a lower information content in the exon regions adjacent to the splice sites than new introns but have a corresponding higher information content in the intron itself. This suggests that introns insert into nonrandom (proto-splice) sites but, during the evolution of an intron after insertion, the splice signal shifts from the flanking exon regions to the ends of the intron itself. Accumulation of information inside the intron during evolution suggests that new introns largely emerge de novo rather than through propagation and migration of old introns.     

Structure and expression of two porcine genomic clones encoding class I MHC antigens

Two nonallelic porcine class I MHC (SLA) genes have been isolated and characterized. Both genes are expressed in mouse L cells, directing the synthesis of class I SLA molecules that carry common monomorphic determinants but are serologically distinct. The corresponding DNA sequences have been determined. The organization of both of these genes is similar to that of other class I genes: a leader exon, three exons encoding extracellular domains, a transmembrane exon, and three intracytoplasmic exons. The two genes are highly homologous in both exon and intron segments, with average homologies of 88% and 80%, respectively. Nucleotide changes in exon 2 are clustered, whereas those in the other exons are dispersed throughout. Comparison of the swine DNA sequences with class I genes from other species reveals a generally high conservation of exons 2, 3, 4, and 6 with lower homology in the remaining protein-encoding domains. Introns are markedly less well conserved, although moderate homology is found between swine and human class I MHC genes in both introns and 3' flanking regions. Taken together with comparisons of the deduced protein sequences, these data indicate an order of swine greater than human greater than rabbit greater than mouse in the relationship of class I genes.     

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.     

Protein domains correlate strongly with exons in multiple eukaryotic genomes--evidence of exon shuffling?

We conducted a multi-genome analysis correlating protein domain organization with the exon-intron structure of genes in nine eukaryotic genomes. We observed a significant correlation between the borders of exons and domains on a genomic scale for both invertebrates and vertebrates. In addition, we found that the more complex organisms displayed consistently stronger exon-domain correlation, with substantially more significant correlations detected in vertebrates compared with invertebrates. Our observations concur with the principles of exon shuffling theory, including the prediction of predominantly symmetric phase of introns flanking the borders of correlating exons. These results suggest that extensive exon shuffling events during evolution significantly contributed to the shaping of eukaryotic proteomes.     

Isolation and characterization of the mouse CD8 beta-chain (Ly-3) genes. Absence of an intervening sequence between V- and J-like gene segments

We have isolated and determined the nucleotide sequence and genomic organization of the genes encoding Ly-3.1 and Ly-3.2. These genes span approximately 14 kb on chromosome 6 and consist of six exons and five introns. The exons correlate roughly with the putative functional domains, namely, a leader exon, a variable and joining region-like exon, a hinge region-like exon, a transmembrane exon, and two intracytoplasmic exons. There is no intervening sequence between V- and J-like gene segments, indicating that rearrangement is not necessary for the expression of the Ly-3 gene. In the 5'-flanking region there is no "TATA" box nor "CAAT" box; however, three "GC" boxes are located upstream of the ATG initiator codon. There are short stretches of sequence homologous to 5'-flanking sequences of the Ly-2 gene. In addition, the sequences CTCTGTGGCA at -748 exhibits homology to the enhancer core sequence of the human Ig H chain and TCR genes. Comparison of the nucleotide sequence corresponding to the extracellular portion between Ly-3.1 and Ly-3.2 revealed a single base difference which results in an amino acid substitution. Therefore it is likely that this amino acid difference is responsible for the previously defined Ly-3 allotypes.     

Distance analysis helps to establish characteristic motifs in intron sequences

Computer-assisted sequence analysis was applied to detect the most apparent nonrandom sequence motifs in eukaryotic introns. We describe in detail a method, which we call distance analysis, that we applied to the extensive study of 405 eukaryotic intron sequences. We observed very strong two-base periodicities for almost all tetranucleotides that are tandem repeats of nonhomopolymeric dinucleotides (the exception was GCGC and CGCG). We also observed, by using a fixed-point alignment method, that these periodic sequence motifs belong to large clusters of dinucleotides repeated tandemly as many as 15–35 times, which corresponds to the cluster lengths of 30–70 bases. We did not observe two-base periodicity of tetranucleotides in the collections of either 262 spliced eukaryotic exons or 107 bacterial genes. Instead, these sequences displayed strong three-base periodicity of some other tetranucleotides. These findings suggest that introns and exons display distinct sequence properties that can be used for mapping purposes.     

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