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.     

Characterization and Gene Organization of Taiwan Banded Krait (Bungarus multicinctus) γ-Bungarotoxin

γ-Bungarotoxin was isolated from Bungarus multicinctus (Taiwan banded krait) venom using a combination of chromatography on a SP-Sephadex C-25 column and a reverse-phase high-performance liquid chromatography column. Circular dichroism (CD) measurement revealed that its secondary structure was dominant with β-sheet structure as is that of snake venom α-neurotoxins and cardiotoxins. γ-Bungarotoxin exhibits activity on inhibiting the binding of [³H]quinuclidinyl benzilate to the M2 muscarinic acetylcholine receptor subtype, and competes weakly with radioiodinated α-bungarotoxin for binding to the Torpedo nicotinic acetylcholine receptor. Moreover, the toxin inhibits collagen-induced platelet aggregation, with an IC₅₀ of approximately 200 nM. The genomic DNA encoding the γ-bungarotoxin precursor is amplified by polymerase chain reaction (PCR). The gene is organized with three exons separated by two introns, and shares virtually identical overall organization with those reported for α-neurotoxin and cardiotoxin genes, including similar intron insertions. The intron sequences of these genes share sequence identity up to 85%, but the exon sequences are highly variable. These observations suggest that γ-bungarotoxin, α-neurotoxins, and cardiotoxins originate from a common ancestor, and the evolution of these genes shows a tendency to diversify the functions of snake venom proteins.     

Hypervariable and Highly Divergent Intron–Exon Organizations in the Chordate Oikopleura dioica

Oikopleura dioica is a pelagic tunicate with a very small genome and a very short life cycle. In order to investigate the intron–exon organizations in Oikopleura, we have isolated and characterized ribosomal protein EF-1, Hox, and -tubulin genes. Their intron positions have been compared with those of the same genes from various invertebrates and vertebrates, including four species with entirely sequenced genomes. Oikopleura genes, like Caenorhabditis genes, have introns at a large number of nonconserved positions, which must originate from late insertions or intron sliding of ancient insertions. Both species exhibit hypervariable intron–exon organization within their -tubulin gene family. This is due to localization of most nonconserved intron positions in single members of this gene family. The hypervariability and divergence of intron positions in Oikopleura and Caenorhabditis may be related to the predominance of short introns, the processing of which is not very dependent upon the exonic environment compared to large introns. Also, both species have an undermethylated genome, and the control of methylation-induced point mutations imposes a control on exon size, at least in vertebrate genes. That introns placed at such variable positions in Oikopleura or C. elegans may serve a specific purpose is not easy to infer from our current knowledge and hypotheses on intron functions. We propose that new introns are retained in species with very short life cycles, because illegitimate exchanges including gene conversion are repressed. We also speculate that introns placed at gene-specific positions may contribute to suppressing these exchanges and thereby favor their own persistence.Supplementary material () is available for this article.     

Characterization and gene organization of Taiwan banded krait (Bungarus multicinctus) gamma-bungarotoxin

-Bungarotoxin was isolated from Bungarus multicinctus (Taiwan banded krait) venom using a combination of chromatography on a SP-Sephadex C-25 column and a reverse-phase high-performance liquid chromatography column. Circular dichroism (CD) measurement revealed that its secondary structure was dominant with -sheet structure as is that of snake venom -neurotoxins and cardiotoxins. -Bungarotoxin exhibits activity on inhibiting the binding of [³H]quinuclidinyl benzilate to the M2 muscarinic acetylcholine receptor subtype, and competes weakly with radioiodinated -bungarotoxin for binding to the Torpedo nicotinic acetylcholine receptor. Moreover, the toxin inhibits collagen-induced platelet aggregation, with an IC₅₀ of approximately 200 nM. The genomic DNA encoding the -bungarotoxin precursor is amplified by polymerase chain reaction (PCR). The gene is organized with three exons separated by two introns, and shares virtually identical overall organization with those reported for -neurotoxin and cardiotoxin genes, including similar intron insertions. The intron sequences of these genes share sequence identity up to 85%, but the exon sequences are highly variable. These observations suggest that -bungarotoxin, -neurotoxins, and cardiotoxins originate from a common ancestor, and the evolution of these genes shows a tendency to diversify the functions of snake venom proteins.     

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.     

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.     

Splicing-Related Features of Introns Serve to Propel Evolution

The role of spliceosomal intronic structures played in evolution has only begun to be elucidated. Comparative genomic analyses of fungal snoRNA sequences, which are often contained within introns and/or exons, revealed that about one-third of snoRNA-associated introns in three major snoRNA gene clusters manifested polymorphisms, likely resulting from intron loss and gain events during fungi evolution. Genomic deletions can clearly be observed as one mechanism underlying intron and exon loss, as well as generation of complex introns where several introns lie in juxtaposition without intercalating exons. Strikingly, by tracking conserved snoRNAs in introns, we found that some introns had moved from one position to another by excision from donor sites and insertion into target sties elsewhere in the genome without needing transposon structures. This study revealed the origin of many newly gained introns. Moreover, our analyses suggested that intron-containing sequences were more prone to sustainable structural changes than DNA sequences without introns due to intron''s ability to jump within the genome via unknown mechanisms. We propose that splicing-related structural features of introns serve as an additional motor to propel evolution.     

Exon junction sequences as cryptic splice sites: implications for intron origin

Introns are flanked by a partially conserved coding sequence that forms the immediate exon junction sequence following intron removal from pre-mRNA. Phylogenetic evidence indicates that these sequences have been targeted by numerous intron insertions during evolution, but little is known about this process. Here, we test the prediction that exon junction sequences were functional splice sites that existed in the coding sequence of genes prior to the insertion of introns. To do this, we experimentally identified nine cryptic splice sites within the coding sequence of actin genes from humans, Arabidopsis, and Physarum by inactivating their normal intron splice sites. We found that seven of these cryptic splice sites correspond exactly to the positions of exon junctions in actin genes from other species. Because actin genes are highly conserved, we could conclude that at least seven actin introns are flanked by cryptic splice sites, and from the phylogenetic evidence, we could also conclude that actin introns were inserted into these cryptic splice sites during evolution. Furthermore, our results indicate that these insertion events were dependent upon the splicing machinery. Because most introns are flanked by similar sequences, our results are likely to be of general relevance.     

Structure and expression of the mouse prealbumin gene

We cloned a genomic DNA fragment which covers the entire sequence of the mouse prealbumin gene and then studied the structure. The coding regions are separated into four exons by three introns, and these numbers, the sizes of the exons and the relative sites of the exon-intron junctions are all in complete agreement with those determined for the human gene. The sequences of four exons can be aligned perfectly with that of the previously determined mouse prealbumin cDNA. In addition to the exon regions, we found two highly conserved DNA regions between the mouse and human prealbumin genes, one in the 5'-flanking region of the gene and the other in the 3' end region of the first intron. These DNA regions contain several consensus glucocorticoid receptor-binding site sequences, and the latter also contains an enhancer sequence present in the immunoglobulin kappa-chain joining-constant kappa intron. RNA hybridizing to the mouse prealbumin cDNA was detected in the extracts from liver, brain, and kidney, but was not detected in testes, spleen, or heart. Little change was caused in the level of prealbumin mRNA in the liver by administration of dexamethasone to mice.     

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.     

昆虫JHE基因结构及与甲基化相关的CpG O/E值分析

本研究选择西方蜜蜂、黑腹果蝇、致倦库蚊和家蚕四种昆虫为例，以降解保幼激素的特异性酯酶（juvenile hormone esterase，JHE）为研究对象，首先运用NCBI的Spidey在线软件将各昆虫的JHE基因与其相应的基因组序列作比对，分别确定各昆虫JHE基因的上游2000bp的具体序列和外显子及内含子区域。发现西方蜜蜂JHE基因含有7个外显子和6个内含子，而其余三种昆虫的JHE基因均含有6个外显子和5个内含子。接着，分别统计基因各区域的CpG，G，C位点的占有量，并计算出CpG O/E值（CpG位点的实际值与期望值之间的比值），发现各区域CpG O/E平均值的大小依次为：基因上游2000bp区>内含子区>外显子区，基因上游2000bp区的CpG O/E平均值大于1.0，而第1到第4外显子的CpG O/E平均值均只有0.8左右。表明在昆虫JHE基因中，若有CpG甲基化位点发生，应主要发生在第1到第4外显子区域。     

Determinants of plant U12-dependent intron splicing efficiency

Factors affecting splicing of plant U12-dependent introns have been examined by extensive mutational analyses in an in vivo tobacco (Nicotiana tabacum) protoplast system using introns from three different Arabidopsis thaliana genes: CBP20, GSH2, and LD. The results provide evidence that splicing efficiency of plant U12 introns depends on a combination of factors, including UA content, exon bridging interactions between the U12 intron and flanking U2-dependent introns, and exon splicing enhancer sequences (ESEs). Unexpectedly, all three plant U12 introns required an adenosine at the upstream purine position in the branchpoint consensus UCCUURAUY. The exon upstream of the LD U12 intron is a major determinant of its higher level of splicing efficiency and potentially contains two ESE regions. These results suggest that in plants, U12 introns represent a level at which expression of their host genes can be regulated.     

Nucleotide diversity and molecular evolution of the WAG-2 gene in common wheat (Triticum aestivum L) and its relatives

In this work, we examined the genetic diversity and evolution of the WAG-2 gene based on new WAG-2 alleles isolated from wheat and its relatives. Only single nucleotide polymorphisms (SNP) and no insertions and deletions (indels) were found in exon sequences of WAG-2 from different species. More SNPs and indels occurred in introns than in exons. For exons, exons+introns and introns, the nucleotide polymorphism π decreased from diploid and tetraploid genotypes to hexaploid genotypes. This finding indicated that the diversity of WAG-2 in diploids was greater than in hexaploids because of the strong selection pressure on the latter. All dn/ds ratios were < 1.0, indicating that WAG-2 belongs to a conserved gene affected by negative selection. Thirty-nine of the 57 particular SNPs and eight of the 10 indels were detected in diploid species. The degree of divergence in intron length among WAG-2 clones and phylogenetic tree topology suggested the existence of three homoeologs in the A, B or D genome of common wheat. Wheat AG-like genes were divided into WAG-1 and WAG-2 clades. The latter clade contained WAG-2, OsMADS3 and ZMM2 genes, indicating functional homoeology among them.     

Nested Evolution of a tRNALeu(UAA) Group I Intron by both Horizontal Intron Transfer and Recombination of the Entire tRNA Locus

The origin and evolution of bacterial introns are still controversial issues. Here we present data on the distribution and evolution of a recently discovered divergent tRNA(Leu)(UAA) intron. The intron shows a higher sequence affiliation with introns in tRNA(Ile)(CAU) and tRNA(Arg)(CCU) genes in alpha- and beta-proteobacteria, respectively, than with other cyanobacterial tRNA(Leu)(UAA) group I introns. The divergent tRNA(Leu)(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNA(Leu)(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNA(Leu)(UAA) gene which is located within a unstable genome region.     

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.