Intron Conservation Research Articles

Common introns within orthologous genes: software and application to plants

The residence of spliceosomal introns within protein-coding genes can fluctuate over time, with genes gaining, losing or conserving introns in a complex process that is not entirely understood. One approach for studying intron evolution is to compare introns with respect to position and type within closely related genes. Here, we describe new, freely available software called Common Introns Within Orthologous Genes (CIWOG), available at http://ciwog.gdcb.iastate.edu/, which detects common introns in protein-coding genes based on position and sequence conservation in the corresponding protein alignments. CIWOG provides dynamic web displays that facilitate detailed intron studies within orthologous genes. User-supplied options control how introns are clustered into sets of common introns. CIWOG also identifies special classes of introns, in particular those with GC- or U12-type donor sites, which enables analyses of these introns in relation to their counterparts in the other genes in orthologous groups. The software is demonstrated with application to a comprehensive study of eight plant transcriptomes. Three specific examples are discussed: intron class conversion from GT- to GC-donor-type introns in monocots, plant U12-type intron conservation and a global analysis of intron evolution across the eight plant species.

Wagner and Dollo: A Stochastic Duet by Composing Two Parsimonious Solos

New contributions toward generalizing evolutionary models expand greatly our ability to analyze complex evolutionary characters and advance phylogeny reconstruction. In this article, we extend the binary stochastic Dollo model to allow for multi-state characters. In doing so, we align previously incompatible Wagner and Dollo parsimony principles under a common probabilistic framework by embedding arbitrary continuous-time Markov chains into the binary stochastic Dollo model. This approach enables us to analyze character traits that exhibit both Dollo and Wagner characteristics throughout their evolutionary histories. Utilizing Bayesian inference, we apply our novel model to analyze intron conservation patterns and the evolution of alternatively spliced exons. The generalized framework we develop demonstrates potential in distinguishing between phylogenetic hypotheses and providing robust estimates of evolutionary rates. Moreover, for the two applications analyzed here, our framework is the first to provide an adequate stochastic process for the data. We discuss possible extensions to the framework from both theoretical and applied perspectives.

Lateral Transfer of a Lectin-Like Antifreeze Protein Gene in Fishes

Fishes living in icy seawater are usually protected from freezing by endogenous antifreeze proteins (AFPs) that bind to ice crystals and stop them from growing. The scattered distribution of five highly diverse AFP types across phylogenetically disparate fish species is puzzling. The appearance of radically different AFPs in closely related species has been attributed to the rapid, independent evolution of these proteins in response to natural selection caused by sea level glaciations within the last 20 million years. In at least one instance the same type of simple repetitive AFP has independently originated in two distant species by convergent evolution. But, the isolated occurrence of three very similar type II AFPs in three distantly related species (herring, smelt and sea raven) cannot be explained by this mechanism. These globular, lectin-like AFPs have a unique disulfide-bonding pattern, and share up to 85% identity in their amino acid sequences, with regions of even higher identity in their genes. A thorough search of current databases failed to find a homolog in any other species with greater than 40% amino acid sequence identity. Consistent with this result, genomic Southern blots showed the lectin-like AFP gene was absent from all other fish species tested. The remarkable conservation of both intron and exon sequences, the lack of correlation between evolutionary distance and mutation rate, and the pattern of silent vs non-silent codon changes make it unlikely that the gene for this AFP pre-existed but was lost from most branches of the teleost radiation. We propose instead that lateral gene transfer has resulted in the occurrence of the type II AFPs in herring, smelt and sea raven and allowed these species to survive in an otherwise lethal niche.

Rare Genomic Characters Do Not Support Coelomata: Intron Loss/Gain

Recently, a new phylogenetic method employing intron position sharing across species was proposed and support for a Coelomate clade reported (Zheng et al. 2007. A rigorous analysis of the pattern of intron conservation supports the Coelomata clade of animals. Mol Biol Evol. 24:2583-2592.). Here, we show that the previous analysis depends on: 1) an idiosyncratic definition of "conserved" introns, 2) exclusion of all phylogenetically informative introns present in outgroups, 3) incorrect inference of change along the critical branch, and 4) lack of variation in rates across branches. The method thus seems unlikely to give accurate results. In addition, we address differences in rates of loss across intron sites, which Zheng et al. claimed invalidates our previous analysis that supported Ecdysozoa (Roy and Gilbert. 2005a. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408.). Instead, we show that our conclusions are likely to be robust to such concerns.

Support for the Coelomata Clade of Animals from a Rigorous Analysis of the Pattern of Intron Conservation

Many intron positions are conserved in varying subsets of eukaryotic genomes and, consequently, comprise a potentially informative class of phylogenetic characters. Roy and Gilbert developed a method of phylogenetic reconstruction using the patterns of intron presence-absence in eukaryotic genes and, applying this method to the analysis of animal phylogeny, obtained support for an Ecdysozoa clade (Roy SW, Gilbert W. 2005. Resolution of a deep animal divergence by the pattern of intron conservation. Proc Natl Acad Sci USA. 102:4403-4408). The critical assumption in the method was the independence of intron loss in different branches of the phylogenetic tree. Here, this assumption is refuted by showing that the branch-specific intron loss rates are strongly correlated. We show that different tree topologies are obtained, in each case with a significant statistical support, when different subsets of intron positions are analyzed. The analysis of the conserved intron positions supports the Coelomata topology, that is, a clade comprised of arthropods and chordates, whereas the analysis of more variable intron positions favors the Ecdysozoa topology, that is, a clade of arthropods and nematodes. We show, however, that the support for Ecdysozoa is fully explained by parallel loss of introns in nematodes and arthropods, a factor that does not contribute to the analysis of the conserved introns. The developed procedure for the identification and analysis of conserved introns and other characters with minimal or no homoplasy is expected to be useful for resolving many hard phylogenetic problems.

Genomix: a method for combining gene-finders' predictions, which uses evolutionary conservation of sequence and intron–exon structure

Correct gene predictions are crucial for most analyses of genomes. However, in the absence of transcript data, gene prediction is still challenging. One way to improve gene-finding accuracy in such genomes is to combine the exons predicted by several gene-finders, so that gene-finders that make uncorrelated errors can correct each other. We present a method for combining gene-finders called Genomix. Genomix selects the predicted exons that are best conserved within and/or between species in terms of sequence and intron-exon structure, and combines them into a gene structure. Genomix was used to combine predictions from four gene-finders for Caenorhabditis elegans, by selecting the predicted exons that are best conserved with C.briggsae and C.remanei. On a set of approximately 1500 confirmed C.elegans genes, Genomix increased the exon-level specificity by 10.1% and sensitivity by 2.7% compared to the best input gene-finder. Scripts and Supplementary Material can be found at http://www.sanger.ac.uk/Software/analysis/genomix

The Emergence of Alternative 3′ and 5′ Splice Site Exons from Constitutive Exons

Alternative 3′ and 5′ splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3′ss and 5′ss exons. The results revealed that alternative 3′ss and 5′ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3′ss and 5′ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.

Origin of the genes for the isoforms of creatine kinase

Creatine kinase (CK) is a member of a family of phosphoryl transfer enzymes called phosphagen (guanidino) kinases which play a central role in cellular energy homeostasis. There are three CK isoform gene groups, each coding for proteins targeted to different intracellular compartments — cytoplasmic (CytCK), mitochondrial (MtCK) and flagellar (FlgCK). The former two CKs are either dimeric or octameric while FlgCKs are contiguous trimers consisting of three fused, complete CK domains. Conventional wisdom supports the view that CKs evolved from a cytoplasmic, monomeric ancestral protein closely related to a phosphagen kinase homologue, arginine kinase (AK). Recently, it has been shown that a demosponge (Phylum Porifera) expresses a true MtCK and two dimeric, protoflagellar CKs (protoflgCK) with great similarity to FlgCKs. To further probe the early evolution of CK, we have obtained additional sequences for Mt- and protoflgCKs from two more demosponges and from three hexactinellid (glass) sponges as well as an MtCK sequence from a basal metazoan cnidarian. Phylogenetic analyses using Maximum Likelihood (ML) of these new CK sequences with other CKs and phosphagen kinases yielded a consensus tree containing an assemblage of MtCKs and a supercluster consisting of protoflg-, Flg- and CytCKs. The MtCKs appear basal in the tree topology consistent with prior results. Within the protoflg-, Flg- and CytCK supercluster, the protoflgCKs appear to be allied to the domains of the FlgCKs, although the support is not robust. PCR amplification of genomic DNA and sequencing of the genes for Mt- and protoflgCK from the demosponge Suberites fuscus showed that the sponge MtCK shares four–five common intron:exon boundaries with invertebrate, protochordate and vertebrate MtCKs supporting a common ancestry and the extreme conservation of intron:exon organization in MtCK genes. The protoflgCK gene organization was highly divergent in relation to other CK genes but shares a common intron:exon boundary with domain 2 of the gene for the FlgCK from the tunicate Ciona intestinalis, providing support for the linkage of the protoflgCKs with the FlgCKs. Our results show that the two, major CK gene lineages are present in arguably the oldest, extant metazoan group, the hexactinellid sponges, indicating that these two genes are ancient and confirming prior work that the MtCK gene is likely basal and ancestral.

Large-scale intron conservation and order-of-magnitude variation in intron loss/gain rates in apicomplexan evolution

The age of modern introns and the evolutionary forces controlling intron loss and gain remain matters of much debate. In the case of the apicomplexan malaria parasite Plasmodium falciparum, previous studies have shown that while the positions of two thirds of P. falciparum introns are not shared with surveyed non-apicomplexans (leaving open the possibility that they were relatively recently gained), 99.1% are shared with Plasmodium yoelii, which diverged from P. falciparum at least 100 Mya. We show here that 60.6% of P. falciparum intron positions in conserved regions are shared with the distantly related apicomplexan Theileria parva, whereas only 18.2% of introns in the more intron-rich T. parva are shared with P. falciparum. Comparison of 3305 pairs of orthologous genes between T. parva and Theileria annulata showed that 7089/7111 (99.7%) introns in conserved regions are shared between species. These levels of conservation imply significant differences in rates of intron loss and gain through apicomplexan history. Because transposable elements (TEs) and/or (often TE-encoded) reverse transcriptase are implicated in models of intron loss and gain, the observed low rates of intron loss and gain in recent Plasmodium and Theileria evolution are consistent with the lack of known TE in those groups. We suggest that intron loss/gain in some eukaryotic lineages may be concentrated in relatively short episodes coincident with occasional TE invasions.

Intronic Alternative Splicing Regulators Identified by Comparative Genomics in Nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.

Evidence of mRNA-Mediated Intron Loss in the Human-Pathogenic Fungus Cryptococcus neoformans

Introns are a defining feature of eukaryotic genomes, though the mechanism of intron gain or loss is not well understood. Reverse transcription of mRNA followed by homologous recombination with the genome has been posited as a mechanism of intron loss, though little direct evidence of recent loss events has been described to support this model. We find supporting evidence for an mRNA-mediated mechanism of loss through comparative genome analyses that revealed a recent loss of 10 adjacent introns in a 22-exon gene in the human-pathogenic fungus Cryptococcus neoformans. We surveyed the gene structures of the entire genomes of Cryptococcus gattii, which diverged from the C. neoformans lineage 37 million years ago (Mya), and C. neoformans var. grubii and var. neoformans, which diverged 18 Mya. Our comparison revealed greater than 99.9% intron conservation, with evidence from 20 genes showing evidence of intron loss, but no convincing evidence of intron gain. Our findings confirm that Cryptococcus introns have been quite stable over recent evolutionary time, with occasional mRNA-mediated intron loss events.

Isolation and characterisation of six putative wheat cell wall-associated kinases

The plant cell wall-associated kinase (WAK) and WAK-like kinase (WAKL) make up a unique group in the receptor-like protein kinase (RLK) superfamily. Previous studies on Arabidopsis have revealed that the WAK gene family members play an important role in both cell elongation and stress response signalling. Here we show that four putative WAKs (TaWAK1, TaWAKL2, TaWAKL3, and TaWAK4) and two WAKLs (TaWAKL1 and TaWAKL2) were isolated from wheat based on the DNA sequence similarity and the protein structure conservation of Arabidopsis WAKs genes. TaWAK1, TaWAK2, TaWAK3 and TaWAKL1 each encode a putative intact protein with the characteristic of the WAK / WAKL gene family members, except for the abbreviated TaWAK4 and TaWAKL2 which were caused by nucleotide mutation and alternative splicing, respectively. Southern analysis revealed that TaWAKL1, TaWAK1, TaWAK2 and TaWAK3 are all multiple-copy members. Real-time PCR analysis revealed that the TaWAK1 and TaWAK3 displayed similar expression patterns, while expressions of TaWAKL1, TaWAKL2, and TaWAK2 were organ specific. Further, we analysed the conservation of introns and intron-exon structure and the putative protein structures between wheat and Arabidopsis, which showed the putative wheat WAKs are different from those of Arabidopsis and make up a new subgroup in the polygenetic tree.

Unusual Intron Conservation near Tissue-Regulated Exons Found by Splicing Microarrays

Alternative splicing contributes to both gene regulation and protein diversity. To discover broad relationships between regulation of alternative splicing and sequence conservation, we applied a systems approach, using oligonucleotide microarrays designed to capture splicing information across the mouse genome. In a set of 22 adult tissues, we observe differential expression of RNA containing at least two alternative splice junctions for about 40% of the 6,216 alternative events we could detect. Statistical comparisons identify 171 cassette exons whose inclusion or skipping is different in brain relative to other tissues and another 28 exons whose splicing is different in muscle. A subset of these exons is associated with unusual blocks of intron sequence whose conservation in vertebrates rivals that of protein-coding exons. By focusing on sets of exons with similar regulatory patterns, we have identified new sequence motifs implicated in brain and muscle splicing regulation. Of note is a motif that is strikingly similar to the branchpoint consensus but is located downstream of the 5′ splice site of exons included in muscle. Analysis of three paralogous membrane-associated guanylate kinase genes reveals that each contains a paralogous tissue-regulated exon with a similar tissue inclusion pattern. While the intron sequences flanking these exons remain highly conserved among mammalian orthologs, the paralogous flanking intron sequences have diverged considerably, suggesting unusually complex evolution of the regulation of alternative splicing in multigene families.

Conservation versus parallel gains in intron evolution

Orthologous genes from distant eukaryotic species, e.g. animals and plants, share up to 25–30% intron positions. However, the relative contributions of evolutionary conservation and parallel gain of new introns into this pattern remain unknown. Here, the extent of independent insertion of introns in the same sites (parallel gain) in orthologous genes from phylogenetically distant eukaryotes is assessed within the framework of the protosplice site model. It is shown that protosplice sites are no more conserved during evolution of eukaryotic gene sequences than random sites. Simulation of intron insertion into protosplice sites with the observed protosplice site frequencies and intron densities shows that parallel gain can account but for a small fraction (5–10%) of shared intron positions in distantly related species. Thus, the presence of numerous introns in the same positions in orthologous genes from distant eukaryotes, such as animals, fungi and plants, appears to reflect mostly bona fide evolutionary conservation.

Resolution of a deep animal divergence by the pattern of intron conservation

The relationship between three biologically important groups, arthropods, nematodes, and deuterostomes, remains unresolved. It is unknown whether arthropods are more closely related to nematodes (consistent with the "ecdysozoa" hypothesis) or to deuterostomes (consistent with "coelomata"). We present a method in which we use the pattern of spliceosomal intron conservation to develop a series of inequalities that characterize each possible relationship. We find that only the ecdysozoa grouping satisfies these predictions, with P < 10(-6). Simulations show that our method, unlike some previous methods, is largely insensitive to rate variation between branches.

Complex early genes.

We use the pattern of intron conservation in 684 groups of orthologs from seven fully sequenced eukaryotic genomes to provide maximum likelihood estimates of the number of introns present in the same orthologs in various eukaryotic ancestors. We find: (i) intron density in the plant-animal ancestor was high, perhaps two-thirds that of humans and three times that of Drosophila; and (ii) intron density in the ancestral bilateran was also high, equaling that of humans and four times that of Drosophila. We further find that modern introns are generally very old, with two-thirds of modern bilateran introns dating to the ancestral bilateran and two-fifths of modern plant, animal, and fungus introns dating to the plant-animal ancestor. Intron losses outnumber gains over a large range of eukaryotic lineages. These results show that early eukaryotic gene structures were very complex, and that simplification, not embellishment, has dominated subsequent evolution.

Unusual Intron Conservation Near Tissue-regulated Exons Found by Splicing Microarrays

Unusual Intron Conservation Near Tissue-regulated Exons Found by Splicing Microarrays

Patterns of Intron Gain and Loss in Fungi

Little is known about the patterns of intron gain and loss or the relative contributions of these two processes to gene evolution. To investigate the dynamics of intron evolution, we analyzed orthologous genes from four filamentous fungal genomes and determined the pattern of intron conservation. We developed a probabilistic model to estimate the most likely rates of intron gain and loss giving rise to these observed conservation patterns. Our data reveal the surprising importance of intron gain. Between about 150 and 250 gains and between 150 and 350 losses were inferred in each lineage. We discuss one gene in particular (encoding 1-phosphoribosyl-5-pyrophosphate synthetase) that displays an unusually high rate of intron gain in multiple lineages. It has been recognized that introns are biased towards the 5′ ends of genes in intron-poor genomes but are evenly distributed in intron-rich genomes. Current models attribute this bias to 3′ intron loss through a poly-adenosine-primed reverse transcription mechanism. Contrary to standard models, we find no increased frequency of intron loss toward the 3′ ends of genes. Thus, recent intron dynamics do not support a model whereby 5′ intron positional bias is generated solely by 3′-biased intron loss.

Structure and function of a cellulase gene in redclaw crayfish, Cherax quadricarinatus

The most abundant organic compound produced by plants is cellulose; however, it has long been accepted that most animals do not produce endogenous enzymes required for its degradation, but rely instead on symbiotic relationships with microbes that produce the necessary enzymes. Here, we present the genomic organisation of an endogenous glycosyl hydrolase family (GHF) 9 gene in redclaw crayfish ( Cherax quadricarinatus), consolidated from a cDNA sequence determined by Byrne et al. [Gene 239 (1999) 317–324.]. Comparison with several other invertebrate GHF9 genes reveals the conservation of both intron position/phase and splice sequence, which adds support to an argument for an ancestral animal cellulase gene. Furthermore, two introns in plant GHF9 genes are also identical in position, implying a more ancient origin for this class of animal cellulase. Protein purification from redclaw gastric fluid via fast performance liquid chromatography (FPLC) indicated the presence of two endoglucanase enzymes. The molecular weights of these components were determined by matrix-assisted laser desorption/ionisation—time-of-flight (MALDI-TOF) to be 47,887 Da (Cel1) and 50,295 Da (Cel2). Cel1 is possibly the functional product of the described cellulase gene, with N-terminal amino acid residues identical to the translated amino acid sequence from the corresponding gene region. Cel2 was identical to Cel1 for 7 of 11 N-terminal residues and likely to be the product of a paralogous endoglucanase gene. These results suggest that redclaw crayfish possess at least one and possibly two functional, endoglucanase enzymes, although further work is required to confirm their origin and attributes.

In NF1, CFTR, PER3, CARS and SYT7, alternatively included exons show higher conservation of surrounding intron sequences than constitutive exons.

It is still not fully understood to what extent intronic sequences contribute to the regulation of the different forms of alternative splicing. We are interested in the regulation of alternative cassette exon events, such as exon inclusion and exon skipping. We investigated these events by comparative genomic analysis of human and mouse in five experimentally well-characterized genes, neurofibromatosis 1 (NF1), cystic fibrosis transmembrane conductance regulator (CFTR), period 3 (PER3), cysteinyl-tRNA synthetase (CARS) and synaptotagmin 7 (SYT7). In NF1, high intron identity around the 52 constitutive and four alternatively skipped NF1 exons is restricted to the close vicinity of the exons. In contrast, we found on average high conservation of intron sequences over 300 base pairs up- and downstream of the five alternatively included NF1 exons. The investigation of alternatively included exons in CFTR, PER3, CARS and SYT7 supported this finding. In contrast, the mean intron identities around the alternatively skipped exons in CTFR and NF1 do not differ considerably from those around the constitutive exons. In these genes, the difference in intron conservation could point to a difference between the regulation of alternative exon inclusion and alternative exon skipping or constitutive exon splicing. Additional genome-wide investigations are necessary to elucidate to what extent our finding can be generalized.