
==== Front
Genome BiolGenome Biology1465-69061465-6914BioMed Central London gb-2004-5-12-r1021557596810.1186/gb-2004-5-12-r102SoftwareThe ASRG database: identification and survey of Arabidopsis thaliana genes involved in pre-mRNA splicing Wang Bing-Bing 1icewater@iastate.eduBrendel Volker 12vbrendel@iastate.edu1 Department of Genetics, Development and Cell Biology Iowa State University, Ames, IA 50011-3260, USA2 Department of Statistics, Iowa State University, Ames, IA 50011-3260, USA2004 29 11 2004 5 12 R102 R102 25 6 2004 6 9 2004 20 10 2004 Copyright © 2004 Wang and Brendel; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The database of Arabidopsis splicing related genes includes classification of genes encoding snRNAs and other splicing related proteins, together with information on gene structure, alternative splicing, gene duplications and phylogenetic relationships.

A total of 74 small nuclear RNA (snRNA) genes and 395 genes encoding splicing-related proteins were identified in the Arabidopsis genome by sequence comparison and motif searches, including the previously elusive U4atac snRNA gene. Most of the genes have not been studied experimentally. Classification of these genes and detailed information on gene structure, alternative splicing, gene duplications and phylogenetic relationships are made accessible as a comprehensive database of Arabidopsis Splicing Related Genes (ASRG) on our website.
==== Body
Rationale
Most eukaryotic genes contain introns that are spliced from the precursor mRNA (pre-mRNA). The correct interpretation of splicing signals is essential to generate authentic mature mRNAs that yield correct translation products. As an important post-transcriptional mechanism, gene function can be controlled at the level of splicing through the production of different mRNAs from a single pre-mRNA (reviewed in [1]). The general mechanism of splicing has been well studied in human and yeast systems and is largely conserved between these organisms. Plant RNA splicing mechanisms remain comparatively poorly understood, due in part to the lack of an in vitro plant splicing system. Although the splicing mechanisms in plants and animals appear to be similar overall, incorrect splicing of plant pre-mRNAs in mammalian systems (and vice versa) suggests that there are plant-specific characteristics, resulting from coevolution of splicing factors with the signals they recognize or from the requirement for additional splicing factors (reviewed in [2,3]).

Genome projects are accelerating research on splicing. For example, with the majority of splicing-related genes already known in human and budding yeast, these gene sequences were used to query the Drosophila and fission yeast genomes in an effort to identify potential homologs [4,5]. Most of the known genes were found to have homologs in both Drosophila and fission yeast. The availability of the near-complete genome of Arabidopsis thaliana [6] provides the foundation for the simultaneous study of all the genes involved in particular plant structures or physiological processes. For example, Barakat et al. [7] identified and mapped 249 genes encoding ribosomal proteins and analyzed gene number, chromosomal location, evolutionary history (including large-scale chromosomal duplications) and expression of those genes. Beisson et al. [8] catalogued all genes involved in acyl lipid metabolism. Wang et al. [9] surveyed more than 1,000 Arabidopsis protein kinases and computationally compared derived protein clusters with established gene families in budding yeast. Previous surveys of Arabidopsis gene families that contain some splicing-related genes include the DEAD box RNA helicase family [10] and RNA-recognition motif (RRM)-containing proteins [11]. At present, the Arabidopsis Information Resource (TAIR) links to more than 850 such expert-maintained collections of gene families [12].

Here we present the results of computational identification of potentially all or nearly all Arabidopsis genes involved in pre-mRNA splicing. Recent mass spectrometry analyses revealed more than 200 proteins associated with human spliceosomes ([13-17], reviewed in [18]). By extensive sequence comparisons using known plant and animal splicing-related proteins as queries, we have identified 74 small nuclear RNA (snRNA) genes and 395 protein-coding genes in the Arabidopsis genome that are likely to be homologs of animal splicing-related genes. About half of the genes occur in multiple copies in the genome and appear to have been derived both from chromosomal duplication events and from duplication of individual genes. All genes were classified into gene families, named and annotated with respect to their inferred gene structure, predicted protein domain structure and presumed function. The classification and analysis results are available as an integrated web resource, the database of Arabidopsis Splicing Related Genes (ASRG), which should facilitate genome-wide studies of pre-mRNA splicing in plants.

ASRG: a database of Arabidopsis splicing-related genes
Our up-to-date web-accessible database comprising the Arabidopsis splicing-related genes and associated information is available at [19]. The web pages display gene structure, alternative splicing patterns, protein domain structure and potential gene duplication origins in tabular format. Chromosomal locations and spliced alignment of cognate cDNAs and expressed sequence tags (ESTs) are viewable via links to the Arabidopsis genome database AtGDB [20], which also provides other associated information for these genes and links to other databases. Text-search functions are accessible from all the web pages. Sequence-analysis tools including BLAST [21] and CLUSTAL W [22] are integrated and facilitate comparison of splicing-related genes and proteins across various species.

Arabidopsis snRNA genes
A total of 15 major snRNA and two minor snRNA genes were previously identified experimentally in Arabidopsis [23-28]. These genes were used as queries to search the Arabidopsis genome for other snRNA genes. A total of 70 major snRNAs and three minor snRNAs were identified by this method. In addition, a single U4atac snRNA gene was identified by sequence motif search. We assigned tentative gene names and gene models as shown in Table 1, together with chromosome locations and similarity scores relative to a representative query sequence. The original names for known snRNAs were preserved, following the convention atUx.y, where x indicates the U snRNA type and y the gene number. Computationally identified snRNAs were named similarly, but with a hyphen instead of a period separating type from gene number (atUx-y). Putative pseudogenes were indicated with a 'p' following the gene name. Pseudogene status was assigned to gene models for which sequence similarity to known genes was low, otherwise conserved transcription signals are missing and the gene cannot fold into typical secondary structure. A recent experimental study of small non-messenger RNAs identified 14 tentative snRNAs in Arabidopsis by cDNA cloning ([29], GenBank accessions 22293580 to 22293592 and 22293600, Table 1). All these newly identified snRNAs were found in the set of our computationally predicted genes.

Conservation of major snRNA genes
As shown in Table 1, each of five major snRNA genes (U1, U2, U4, U5 and U6) exists in more than 10 copies in the Arabidopsis genome. U2 snRNA has the largest copy number, with a total of 18 putative homologs identified. Both U1 and U5 snRNAs have 14 copies, U6 snRNA has 13 copies, and U4 snRNA has only 11 copies. Sequence comparisons within Arabidopsis snRNA gene families showed that the U6 snRNA genes are the most similar, and the U1 snRNA genes are the most divergent. Eight active U6 snRNA copies are more than 93% identical to each other in the genic region, whereas active U1 snRNAs are on average only 87% identical. The U2 and U4 snRNAs are also highly conserved within each type, with more than 92% identity among the active genes. Details about the individual snRNAs and the respective sequence alignments are displayed at [30].

Previous studies identified two conserved transcription signals in most major snRNA gene promoters: USE (upstream sequence element, RTCCACATCG (where R is either A or G) and TATA box [24-27]. All 14 U5 snRNAs have the USE and TATA box. Furthermore, their predicted secondary structures are similar to the known structure of their counterparts in human, indicating that all these genes are active and functional (structure data not shown; for a review of the structures of human snRNAs, see [31]). Similarly, we identified 17 U2, 10 U1, nine U4, and nine U6 snRNA genes as likely active genes, with a few additional genes more likely to be pseudogenes because of various deletions. U4-10 and U6-7 do not have the conserved USE in the promoter region, but their U4-U6 interaction regions (stem I and stem II) are fairly well conserved. U2-16 is also missing the USE but has a secondary structure similar to other U2 snRNAs. These genes may be active, but differences in promoter motifs suggest that their expression may be under different control compared with other snRNAs homologs. The U2-17 snRNA has all conserved transcription signals, but 20 nucleotides are missing from its 3' end. The predicted secondary structure of U2-17 is similar to that of other U2 snRNAs, with a significantly shorter stem-loop in the 3' end as a result of the deletion. We are not sure if the U2-17 snRNA is functional, but the conserved transcription signals imply that it may be active.

Other conserved transcription signals were also identified in most active snRNAs, including the sequence element CAANTC (where N is either A, C, G or T) in U2 snRNAs (located at -6 to -1) [23], and the termination signal CAN3-10AGTNNAA in U snRNAs (U1, U2, U4 and U5) transcribed by RNA polymerase II (Pol II) [23,24,32]. The previously identified monocot-specific promoter element (MSP, RGCCCR, located upstream of USE) in U6.1 and U6.26 [33] is also found in five other U6 snRNA genes (U6.29, U6-2, U6-3, U6-4, U6-5). In all seven U6 snRNAs the consensus MSP sequence extends by two thymine nucleotides to RGCCCRTT. Although the MSP does not contribute significantly to U6 snRNA transcription initiation in Nicotiana plumbaginifolia protoplasts [33], the extended consensus may imply a role in gene expression regulation in Arabidopsis.

Low copy number of minor snRNA genes
The minor snRNAs are functional in the splicing of U12-type (AT-AC) introns. Four types of minor snRNAs, which correspond to four types of major snRNAs, exist in mammals. U11 is the analog of U1, U12 is the analog of U2, U4atac is the analog of U4, and U6atac is the analog of U6. The U5 snRNA seems to function in both the major and minor spliceosome [34]. Two minor snRNAs (atU12 and atU6atac) were experimentally identified in Arabidopsis [28]. Both have the conserved USE and TATA box in the promoter region. We identified another U6atac gene (atU6atac-2) by sequence mapping. This gene has a USE and a TATA box in the promoter region. The atU6atac-2 gene is more than 90% similar to atU6atac in both its 5' and 3' ends, with a 10-nucletotide deletion in the central region. The putative U4atac-U6atac interaction region in atU6atac-2 is 100% conserved with the interaction region previously identified in atU6atac [28,35].

U11 and U4atac have not been experimentally identified in Arabidopsis. BLAST searches using the human U11 and U4atac homologs as queries against the Arabidopsis genome failed to find any significant hits, indicating divergence of the minor snRNAs in plants and mammals. Using the strategy described below, we successfully identified a putative Arabidopsis U4atac gene. It is a single-copy gene containing all conserved functional domains. We also found a single candidate U11 snRNA gene (chromosome 5, from 17,492,101 to 17,492,600) that has the USE and TATA box in the promoter region. This gene also contains a putative binding site fr Sm protein and a region that could pair with the 5' splice site of the U12-type intron.

Identification of an Arabidopsis U4atac snRNA gene
Like U4 snRNA and U6 snRNA, human U4atac and U6atac snRNAs interact with each other through base pairing [36]. The same interaction is expected to exist between the Arabidopsis homologs. Therefore, we deduced the tentative AtU4atac stem II sequence (CCCGTCTCTGTCAGAGGAG) from AtU6atac snRNA and searched for matching sequences in the Arabidopsis genome. Hit regions together with flanking regions 500 base-pairs (bp) upstream and 500 bp downstream were retrieved and screened for transcription signals (USE and TATA box). One sequence was identified that contains both the USE and TATA box in appropriate positions, as shown in Figure 1.

The tentative U4atac snRNA gene contains not only the stem II sequence, but also the stem I sequence that presumably base-pairs with U6atac snRNA stem I. Furthermore, a highly conserved Sm-protein-binding region exists at the 3' end. The predicted secondary structure is nearly identical to hsU4atac, with a relative longer single-stranded region (data not shown). With the highly conserved transcriptional signals, functional domains and secondary structure, this candidate gene is likely to be a real U4atac snRNA homolog. We named it AtU4atac and assigned At4g16065 as its tentative gene model because it is located between gene models At4g16060 and At4g16070 on chromosome 4.

Tandem arrays of snRNAs genes
Some snRNAs genes exist as small groups on the Arabidopsis chromosomes [6]. We identified 10 snRNA gene clusters: seven U1-U4 snRNA clusters, one U2-U5 snRNA cluster, and a tandem duplication for both U2 snRNA (U2-10) and U5 snRNA (U5.1) (Figure 2). All seven Arabidopsis U1-U4 clusters have the U1 snRNA gene located upstream of the U4 snRNA gene, with a 180-300-nucleotide intergenic region. Five of the U1-U4 arrays are located on chromosome 5 (U1a/U4.1, U1-4/U4-5, U1-8/U4-7, U1-9/U4-8, and U1-13p/U4.3p), and the remaining two on chromosome 1 (U1-10/U4-6 and U1-14p/U4-10). The U2-17 and U5-10 occur in tandem array on chromosome 5, separated by fewer than 200 nucleotides.

Arabidopsis splicing-related protein-coding genes
Most of the proteins involved in splicing in mammals and Drosophila are known [4,37,38]. In addition, recent proteomics studies revealed many novel proteins associated with human spliceosomes (reviewed in [18]). Using all these animal proteins as query sequences, we identified a total of 395 tentative homologs in Arabidopsis. Sequence-similarity scores and comparison of gene structure and protein domain structure were used to assign the genes to families. Each gene was assigned a tentative name based on the name of its respective animal homolog. Different homologs within a gene family were labeled by adding an Arabic number (1, 2, and so on) to the name. Close family members with similar gene structure were indicated by adding -a, -b, and -c to the name. The 395 genes were classified into five different categories according to the presumed function of their products. Ninety-one encode small nuclear ribonucleoprotein particle (snRNP) proteins, 109 encode splicing factors, and 60 encode potential splicing regulators. Details of EST evidence, alternative splicing patterns, duplication sources and domain structure of these genes are listed in Table 2. We also identified 84 Arabidopsis proteins corresponding to 54 human spliceosome-associated proteins. The remaining 51 genes encode proteins with domains or sequences similar to known splicing factors, but without enough similarity to allow unambiguous classification. These two categories are not discussed in detail here, but information about these genes is available at our ASRG site [39].

The majority of snRNP proteins are conserved in Arabidopsis
There are five snRNPs (U1, U2, U4, U5 and U6) involved in the formation of the major spliceosome, corresponding to five snRNAs. Five snRNPs (U1 snRNP, U2 snRNP, U5 snRNP, U4/U6 snRNP and U4.U6/U5 tri-snRNP) have been isolated experimentally in yeast or human [40-45]. Each snRNP contains the snRNA, a group of core proteins, and some snRNP-specific proteins. Most of these proteins are conserved in Arabidopsis. All U snRNPs except U6 snRNP contain seven common core proteins bound to snRNAs. These core proteins all have an Sm domain and have been called Sm proteins. The U6 snRNP contains seven LSM proteins ('like Sm' proteins). Another LSM protein (LSM1) is not involved in binding snRNA (reviewed in [46]).

As shown in Table 2, all Sm and LSM proteins have homologs in Arabidopsis, and eight of them are duplicated. It is likely that these genes existed as single copies in the ancestor of animals and plants, but duplicated within the plant lineage. Only one of the 24 genes (LSM5, At5g48870) has been characterized experimentally in Arabidopsis. The LSM5 gene was cloned from a mutant supersensitive to ABA (abscisic acid) and drought (SAD1 [47]). LSM5 is expressed at low levels in all tissues and its transcription is not altered by drought stress [47]. cDNA and EST evidence exist for all other core protein genes, indicating that all 24 genes are active.

There are 63 Arabidopsis proteins corresponding to the 35 snRNP-specific proteins used as queries in our genome mapping. Very few of them have been characterized experimentally, including U1-70K, U1A and a tandem duplication pair of SAP130 [48-50]. U1-70K was reported as a single-copy essential gene. Expression of U1-70K antisense transcript under the APETALA3 promoter suppressed the development of sepals and petals [51]. We identified an additional homolog of U1-70K (At2g43370) and named it U1-70K2. The U1-70K2 proteins showed 48% similarity to the U1-70K protein according to Blast2 results. Both genes retain the sixth intron in some transcripts, a situation which would produce truncated proteins [48]. Interestingly, we found that five of the 10 Arabidopsis U1 snRNP proteins, including the U1-70K-coding genes, may undergo alternative splicing.

Several genes in U2, U5, U4/U6 and U4.U6/U5 snRNPs, but none in U1 snRNP, occur in more than three copies in the Arabidopsis genome. The atSAP114 family has five members, including two that occur in tandem (atSAP114-1a and atSAP114-1b). Three members have EST/cDNA evidence (Table 2). Interestingly, the predicted atSAP114p (At4g15580) protein contains a RNase H domain at the amino-terminal end, and thus atSAP114p shares similarity to At5g06805, a gene annotated as encoding a non-LTR retroelement reverse transcriptase-like protein. It is likely that the atSAP114p gene is a pseudogene that originated by retroelement insertion. There are three copies of the gene for the tri-snRNP 65 kilodalton (kDa) subunit, which are clustered on chromosome 4. Both the U4/U6 90 kDa protein and the U4/U6 15.5 kDa protein also have three gene copies, and the 116 kDa and 200 kDa subunits in U5 snRNP have four copies apiece.

The yeast U1 snRNP contains several specific proteins that are not present in mammalian U1 snRNPs [52]. As in mammals, Arabidopsis also lacks homologs of Prp42, a component of U1 snRNP in yeast [53]. However, Arabidopsis has two copies of the gene for Prp39, which are similar to Prp42. Furthermore, atPrp39a can produce a shorter protein isoform with a novel amino-terminal sequence by exon skipping. It is possible that the duplicates and alternative isoforms of plant U1 snRNP proteins are functional homologs of the yeast-specific proteins.

Several proteins specific to the minor spliceosome are also conserved in Arabidopsis. The human 18S U11/U12 snRNP contains several proteins found in U2 snRNP as well as seven novel proteins [14]. Four of the seven U11/U12-specific proteins (U11/U12-35K, 25K, 65K and 31K) are conserved in Arabidopsis, while the remaining three (59K, 48K and 20K) have no clear homologs. Interestingly, all four Arabidopsis genes are single copy in the genome, and three of them are apparently alternatively spliced (Table 2).

Splicing factors are slightly different in Arabidopsis than in other organisms
We divided the splicing factors into eight subgroups according to recent human spliceosome studies [13,14,16,18]: splice-site selection proteins; SR proteins; 17S U2 associated proteins; 35S U5 associated proteins; proteins specific to the BΔU1 complex; exon junction complex (EJC) proteins; second-step splicing factors and other known splicing factors. We focused our analysis on the first two subgroups because their functions in splicing are well established. A total of 109 proteins in Arabidopsis were identified, corresponding to 67 human queries from all eight subgroups. Most of the proteins are conserved among eukaryotes, but some human proteins have no obvious homologs in the Arabidopsis genome, and some novel splicing factors appear to exist in Arabidopsis. About 43% of genes encoding splicing factors are duplicated in the genome, whereas some proteins, such as SF1/BBP (branchpoint-binding protein, which facilitates U2 snRNP binding in fission yeast [54]) and cap-binding proteins (CBP20 and CBP80, possibly involved in cap proximal intron splicing [55]), derive from single-copy genes [56]. These single-copy gene products may work with all pre-mRNAs, including the ones with U12-type introns. Surprisingly, mutation of CBP80 (ABH1) is not lethal and is non-pleiotropic. The abh1 plants show ABA-hypersensitive closure of stomata and reduced wilting during drought [57].

Many splicing factors have been identified previously in Arabidopsis, including two U2AF65, two U2AF35, and 18 SR proteins [58-67]. The U2AF35-related protein atUrp, which could interact with U2AF65 and position RS-domain-containing splicing factors [68], is also present in the Arabidopsis genome. Although the Urp gene is expressed ubiquitously in human tissues, no ESTs from this gene were found in Arabidopsis. Three copies of PTB/hnRNP-I genes were identified in Arabidopsis. The PTB protein competes for the poly-pyrimidine tract with the U2AF large subunit, thus negatively regulating splicing [69].

We also identified a homolog related to atU2AF65 (At2g33440) and an additional SR protein (At2g46610). The U2AF65-related protein (atULrp, At2g33440) has 247 amino acids and shares over 40% similarity with the carboxy terminal region of the two atU2AF65 homologs. Only one RRM can be identified in atULrp, in contrast to three RRMs and one amino-terminal RS domain in atU2AF65 proteins, and there is no apparent RS domain in atULrp. No animal homolog of atULrp could be identified. The function of this one-RRM U2AF65-related protein is not clear. As it lacks other functional motifs, it might act as a competitor of U2AF65. A two-RRM U2AF65 protein can be produced through alternative splicing. The 11th intron of atU2AF65a can be retained (see RAFL full-length cDNA, gi:19310596) to produce a truncated protein with only the first two RRMs. Interestingly, the last RRM in atU2AF65a contains several amino-acid variations from the consensus pattern such that it could not be detected by InterPro and NCBI-CDD searches using default values, also suggesting that perhaps only the first two RRMs are essential.

The additional SR protein belongs to the atRSp31 family and was named atRSp32 (At2g46610). It shares 70% identity and 78% similarity with atRSp31. The protein is 250 amino acids in length and contains two RRMs and some RS dipeptides in the carboxy-terminal region. The gene structure of atRSp32 is similar to that of atRSp31. Two other genes (atRSp40 and atRSp41) are in the same family and also have similar exon and intron sizes (see gene structure information at [70]). Similarly to the previous classification of 18 SR proteins [61], the 19 SR proteins (including SR45) can be grouped into four large families of four to five members according to sequence similarity, gene structure and protein domain structure.

The atRSp31 family (atRSp31, atRSp32, atRSp40 and atRSp41) belongs to a novel plant SR family and has no clear animal ortholog. Other families include the SC35 (or SRrp/TASR2) family, SF2/ASF family, and the 9G8 family. Arabidopsis has a single copy of the SC35 ortholog and four SC35-like proteins (atSR33, atSCL30a, atSCL30 and atSCL28), which appear to have diverged significantly from SC35. It seems that this divergence predates the split of plants and animals because a similar SC35-like gene family exists in the human genome (SRrp35 and SRrp40). The SRrp35 and SRrp40 were found to antagonize other SR proteins in vitro and function in 5' splice-site selection [71]. SF2/ASF has four copies (atSR1/SRp34, atSRp30, atSRp34a and atSRp34b) with similar gene structures and domains. Human 9G8 protein has five homologs in Arabidopsis, with three (atRSZp21, atRSZp22 and atRSZp22a) containing one CCHC-type zinc finger and two (atRSZ33, atRSZ32) containing two CCHC-type zinc fingers in addition to an RRM and an RS domain. Interestingly, several SR proteins (atRSZp21, atRSZp22, SR45 and SCL33) were found to interact with atU1-70K, and some SR proteins can interact with each other, thus forming a complicated interaction network to facilitate splice-site selection and spliceosome assembly [3,61-63]. atSR45 was initially regarded as a novel plant SR protein [63], but by virtue of sequence-similarity scores it actually may be the ortholog of the human RNPS1 gene, which encodes an EJC protein. Other human SR proteins (SRp20, SRp30c, SRp40, SRp54, SRp55 and SRp75) lack clear orthologs in Arabidopsis. We conclude that SR protein families evolved differently in animals and plants from three to four common ancestors, including SC35, SF2/ASF and 9G8/RSZ. The SRrp (SC35-like in plants) family may not be classical SR proteins but they play important roles in splice-site selection.

Proteins in other subgroups, such as 17S U2 snRNP-associated proteins, 35S U5 snRNP-associated proteins, and protein specific to the BΔU1 complex, are also conserved in Arabidopsis. The BΔU1 complex is the spliceosome complex captured immediately before catalytic activation. Most proteins in the 35S U5 snRNP are absent in the BΔU1 complex but present in the active B complex, indicating the important roles of 35S U5 snRNP-associated proteins in spliceosome activation [13]. Conservation of these proteins in Arabidopsis revealed the same pathway of spliceosome activation in plants. A subcomplex named Prp19 complex in 35S U5 snRNP has a critical role in spliceosome activation [13,72]. All proteins in the human Prp19 complex have homologs in Arabidopsis, including a chromosomal duplication pair of Prp19 genes and a single copy of the CDC5 gene. For the BΔU1 complex, six human genes have homologs, and five of them are single copy in Arabidopsis. Two genes (NPW38BP/SNP70 and p220(NPAT)) in the human BΔU1 complex have no apparent Arabidopsis homologs.

Arabidopsis also lacks an SMN protein complex. In human, the SMN protein (survival of motor neurons) can interact with a series of proteins including Gemin2, Gemin3 (a helicase), Gemin4, Gemin5 and Gemin6 to form an SMN complex, which has important roles in the biogenesis of snRNPs and the assembly of the spliceosome through direct interactions with Sm proteins and snRNA [73]. Although the SMN protein exists in the fission yeast genome (GenBank accession CAA91173), no SMN complex members can be identified in the Arabidopsis genome.

Splicing regulators are expanded in Arabidopsis
Splicing regulators are proteins that can either modify splicing factors or compete with splicing factors for their binding site. Important splicing regulators are hnRNP proteins and SR protein kinases. The exact role of phosphorylation of SR proteins in splicing is not yet clear, but SR protein kinases are well conserved and exist as multiple copies in Arabidopsis. A total of eight SR protein kinases were identified in Arabidopsis, including three Lammer/CLK kinases (AFC1, AFC2 and AFC3), two SRPK1 homologs, and three SPRK2 homologs. The three Lammer/CLK kinases were identified previously, and AFC2 was shown to phosphorylate SR protein in vitro [63,74]. Overexpression of tobacco AFC2 homolog PK12 in Arabidopsis changed the alternative splice patterns of several genes, including atSRp30, atSR1/atSRp34 and U1-70K [75], indicating that these SR proteins may function to modulate splicing in plants.

The heterogeneous nuclear ribonucleoproteins (hnRNPs) bind to splice sites and to binding sites for splicing factors on nascent pre-mRNAs, thus competing with splicing factors to negatively control splicing (reviewed in [76]). Humans have about 20 hnRNP proteins, many of which function in splicing. A total of 35 potential hnRNP proteins possibly related to splicing was found in Arabidopsis by sequence-similarity searches, including a superfamily of glycine-rich RNA-binding proteins. This family contains 21 members similar to human hnRNP A1 and hnRNP A2/B1. It can be further divided into two subfamilies. One includes eight proteins containing one RRM, and another has 13 members with two RRMs. 12 of these proteins were identified previously, including AtGRP7, AtGRP8, UBA2a, UBA2b, UBA2c and AtRNPA/B1-6 [11,77,78]. AtGRP7 was found to be able to influence alternative splicing of its own transcripts as well as AtGRP8 transcripts [79]. UBA2 proteins can interact with UBP1 and UBA1 proteins, which have three RRMs and one RRM respectively, to recognize U-rich sequences in the 3' untranslated region (UTR) and stabilize mRNA [78]. Although the overexpression of UBA2 did not stimulate splicing of a reporter gene in tobacco protoplasts [78], we cannot rule out the possibility that it could be involved in splicing of other genes.

Other human hnRNPs related to splicing also have homologs in Arabidopsis. BLAST searches of the human (CUG)n triplet repeat RNA-binding protein (CUG-BP) against all Arabidopsis proteins revealed three putative homologs, including atFCA. atFCA and CUG-BP share similarity within the RRMs and a region approximately 40 amino acids in length. An additional protein (At2g47310) related to FCA was identified and named FCA2, as it shares about 50% similarity with FCA. The FCA proteins have two RRMs and a WW domain, which interact with the FY protein, a homolog of yeast polyadenylation factor Psf2p [80,81]. The FCA-FY complex negatively regulates the FCA protein by favoring a polyadenylation site from the third intron of FCA pre-mRNA [80,82]. FCA may be a multifunctional protein involved in mRNA processing, as human CUG-BP can function in both alternative splicing and deadenylation [83]. We also list 15 previously identified hnRNP-like proteins and two additional homologs as possible splicing regulators. The UBP1 proteins can strongly enhance splicing of some introns in protoplasts [84], whereas UBA1, RBP45 and RBP47 proteins have no similar function [78,85].

Unclassified splicing protein candidates
In addition to the 260 proteins in the above three categories, there are also 84 Arabidopsis proteins corresponding to human spliceosome-associated proteins identified in recent proteomic studies [15-18]. Some of these proteins function in other processes, such as transcription, polyadenylation and even translation. Their association with spliceosomes provides evidence for the coupling of splicing and other processes. Other proteins have no known functions. Only 35.8% of the proteins in this category are duplicated in Arabidopsis. We also identified a total of 51 Arabidopsis protein-coding genes similar to known splicing proteins. They have conserved domains and some level of sequence similarity to known splicing factors. We did not include these two categories in Table 2, but detailed information about them is available at ASRG [39].

Distribution and duplication of Arabidopsis splicing-related genes
The distribution of Arabidopsis snRNA and splicing-related proteins across the genome is shown in Figure 2 and at the ASRG website. Overall, the genes appear evenly distributed on the chromosomes, with several small gene clusters. Only four snRNA genes are located on chromosome 2, three of which are U2 snRNA genes. No U4 snRNA gene is located on chromosome 4. For the protein-coding genes, most functional categories have members located on each chromosome. The only exception is the SR protein kinase family, which has no member on chromosome 1. Interestingly, chromosome 1 contains the most snRNP proteins and splicing factors, but has the fewest splicing regulators. Several gene clusters encoding splicing-related proteins were also identified. Some clusters, such as tandemly duplicated gene pairs, include genes from the same category. One cluster located on chromosome 4 includes four genes encoding tri-snRNP proteins (atTri65a, atTri65b, atTri65c and atTri15.5c, homologs of tri-snRNP 65-KD protein and 15.5 KD protein). Two other clusters, atU2A-atCdc5 and atCUG-BP1-atU1C, include genes from different functional categories. No clear clusters of genes for snRNA-splicing-related proteins were identified. Although about one third of snRNA genes are located near other protein-coding genes, none of their neighboring genes is related to splicing. As a caveat, we should point out that our snRNA gene determination strongly suggests annotation errors in overlapping protein-coding gene models. Thus, atU2-1, atU2.3, atU4.2, atU4-11p, atU5-13 and atU6.26 overlap gene models At1g16820, At3g57770, At3g06895, At1g68390, At5g53740 and At3g13857, respectively, but none of these models is well supported by cDNA or EST evidence (see displays linked at ASRG [30]).

The 260 proteins in the first three categories could be grouped into 130 families, 66 of which consist of multiple members. The duplication rate is over 50%, which is higher than the 44% duplication rate of Arabidopsis transcription factors [86]. As shown in Table 3, about 50% of genes encoding snRNP proteins, 43% of splicing factors, and 78% of splicing regulators have duplications. The much higher duplication rate of splicing regulators may reflect diversification in splicing control.

At least 130 duplication events are required to yield the 260 proteins from 130 families given one single-copy ancestor per family. Thirty-three duplication events (about a quarter of the total) are likely to be the result of chromosome duplications. The chromosomal duplication ratio is 18.9-27.5% among the three groups (see Table 3). Some snRNA genes pairs, such as U2-14/U2-10, U5-3/U5-5 and (U6.1 U6.26)/(U6-8p U6-9p), may also have been produced by chromosome duplication. The C.D.2-3 region (chromosome duplication region between chromosomes 2 and 3, see [87]) has the most splicing-related gene pairs. Six genes in this region on chromosome 2 were duplicated in the same order on chromosome 3. EST evidence shows that all these genes are expressed. Three U5 snRNA genes (U5.1, U5.1b and U5-4) and four U2 snRNA genes (U2.2, U2.3, U2.4 and U2.6) also are located in the same region on chromosome 3. No U5 and U2 homologs exist in the corresponding region on chromosome 2, suggesting that the snRNA duplication events in that region may have happened after the chromosome duplication event, or that the snRNA duplicates were lost subsequent to chromosome duplication.

Chromosomal duplication rather than individual gene duplication appears to be the predominant mode of amplification for some types of genes. As shown in Table 2, the 24 genes encoding core proteins have nine duplication pairs, five of which can be attributed to chromosomal duplications. The 19 SR protein genes include eight duplication pairs, six of which are probably the results of chromosomal duplications. At least five chromosomal duplication events contributed to the superfamily of 21 hnRNP glycine-rich RBD and A/B genes. It is not clear why these functional categories have high chromosomal duplication ratios. It is possible that chromosomal duplication could create positive selection to maintain similar copy numbers of other genes encoding proteins that interact with the products of already duplicated genes.

Alternative splicing of Arabidopsis splicing-related genes
According to EST/cDNA alignments, 80 of the 260 protein coding genes show 66 alternative splicing events. This rate (30.8%) is much higher than the overall frequency of alternative splicing in Arabidopsis, which is about 13% using the same criteria (2,747 genes out of 20,446 genes with EST/cDNA evidence; B.-B.W. and V.B., unpublished work). As shown in Table 4, the snRNP protein-coding genes have the lowest alternative splicing ratio (24.2%), whereas the ratios for splicing factor and splicing regulator genes are both over 33%. More than half of the genes encoding EJC proteins, proteins specific for the BΔU1 complex, SR proteins, U11/U12 snRNP-specific proteins and U1 snRNP proteins undergo alternative splicing.

Among different types of alternative splicing, intron retention is the most abundant of the alternative transcripts identified for the 260 classified splicing-related genes. As shown in Table 4, 44 of the total 80 alternative splicing genes (about 55%) involve intron retention, 28 (35%) involve alternative acceptor-site selection and 15 (18.7%) are due to exon skipping. Compared with the corresponding ratio in all Arabidopsis alternative splicing events (55.3% intron retention, 23.4% alternative acceptor-site selection and 6.3% exon skipping; B.-B.W. and V.B., unpublished work), the ratio of intron retention in splicing-related genes is similar and the ratio of exon skipping is higher. Interestingly, only one of the 20 splicing regulator genes processed by alternative splicing (about 5%) shows exon skipping, indicating that exon skipping is an important post-transcriptional method for controlling the expression of splicing factor coding genes but not the splicing regulator genes.

Discussion
Previous studies had determined 30 snRNA genes and 46 protein-coding genes related to splicing in Arabidopsis (see Tables 1 and 2). In this study, we have computationally identified an additional 44 snRNA genes (Table 1) and 349 protein-coding genes (Table 2) that also may be involved in splicing. Among the five types of U snRNAs, U6 is the most conserved and U1 is the least conserved. We identified seven U1-U4 snRNA gene clusters. We were surprised to see so many U1-U4 clusters in Arabidopsis. In Drosophila, four snRNA clusters were reported [4], but none of them includes U1-U4 gene pairs. It is likely that a U1-U4 snRNA cluster existed in a progenitor of the current Arabidopsis genome, which was duplicated several times to form the extant seven clusters. The non-clustered U1 and U4 snRNA genes may have arisen by individual gene duplication or gene loss in duplicated clusters.

Among the proteins involved in splicing, most animal homologs are conserved in plants, indicating an ancient, monophylytic origin for the splicing mechanism. A striking feature of plant splicing-related genes is their duplication ratio. Fifty percent of the splicing genes are duplicated in Arabidopsis. The duplication ratio of the splicing-related genes increases from genes encoding snRNP proteins to genes encoding splicing regulators. These data strongly suggest that the general splicing mechanism is conserved, but that the control of splicing may be more diverse in plants.

The high duplication ratio of Arabidopsis splicing-related genes could be the result of evolutionary selection. Unlike animals, which can move around to maintain more homogeneous physiological conditions, plants are exposed to a larger range of stress conditions such as heat and cold. The duplicates will more probably be maintained in the genome as their functions become diversified, and potentially plant-specific, to ensure the fidelity of splicing under such varied conditions. Chromosome duplication has produced several Sm proteins, SR proteins and hnRNP proteins in Arabidopsis, which in turn could create positive selective pressures influencing the rate of duplication for functionally related genes. Because chromosome duplication occurred differentially within each plant lineage, we would expect different duplication patterns of these genes in, for example, monocots and dicots.

To confirm the above hypothesis, we searched the recently sequenced rice genome using the five Arabidopsis SC35 and SC35-like proteins as probes. Eight distinct genome loci were found to encode SC35 and SC35-like proteins, including three homologs of atSC35, two homologs of atSR33/SCL33 and atSCL30a, two homologs of atSCL30, and one homolog of atSCL28. Five of the eight rice genes are currently annotated in GenBank with accession numbers BAC79909 (osSC35a), BAD09319 (osSC35b), AAP46199 (osSR33-1), BAC799901 (osSCL30a/osSR33-2), and BAD19168 (osSCL30-1). As shown in the phylogenetic tree displayed in Figure 3, the two rice SC35 genes and atSC35 are likely to be orthologs of the animal SC35 gene. The other sequences cluster in SC35-like (SRrp/TASR) clades, indicating that the SC35 and SRrp/TASR genes diverged before the divergence of monocot and dicot plants (the divergence presumably happened even before the divergence of animals and plants, as described earlier). In addition, there are species-specific duplications. Thus, the Arabidopsis chromosomal duplication pair atSR33 - atSCL30a forms a clade, while their rice copies (osSR33-1 and osSCL30a) form another clade. Also there are additional duplications for the rice SC35 and SCL30 genes. We are currently working to identify all rice splicing related genes. The complete sets of these genes in two plant species should provide a good foundation for assessing similarities and differences in splicing mechanisms used by monocot and dicot plants.

As introns evolve rapidly, the mechanism to recognize and splice them should either evolve correspondingly or be flexible enough to accommodate the changes. It seems that plants deploy the most economic and practical way by keeping a largely conserved splicing mechanism and a very flexible recognition and control mechanism. Direct evidence comes from the presence of plant-specific splicing proteins, such as the novel SR protein family and the superfamily of hnRNP A/B. The absence of SMN complex and some yeast U1 snRNP proteins in Arabidopsis indicates that other organisms also have integrated new proteins or pathways into the splicing mechanism over the course of evolution relative to other eukaryotes. Other evidence supporting the conserved splicing but flexible regulating mechanism include differential conservation among U snRNAs (U1 snRNAs are less conserved than U6 snRNAs) and high alternative splicing frequency in U1 snRNP proteins, SR proteins and hnRNP proteins. The SR proteins and U1 snRNP proteins are involved in early steps of splicing and 5' and 3' splice-site selection; multiple isoforms of these proteins may be functionally significant in the control of splicing.

It is interesting to note that the overall alternative splicing frequency in splicing related genes is much higher than the frequency averaged over all Arabidopsis genes. More than half of SR proteins and U1 snRNP proteins show alternative splicing. Alternative splicing might increase protein diversity derived from splicing-related genes, which would further add flexibility to the splicing mechanism. The high frequency of alternative transcripts from splicing related genes raises another interesting question - how is splicing regulated in these splicing-related genes? One possible answer is that some splicing-related genes may be autoregulated. Accumulation of one transcript would feed back to inhibit/promote other isoforms. Several splicing-related genes have been reported to be regulated in this way. For example, AtGRP7 (hnRNP A/B superfamily) is a circadian clock-regulated protein which negatively autoregulates its expression [79]. When the AtGRP7 protein accumulates over the circadian cycle, it promotes production of alternative transcripts which use a cryptic 5' splice site. As a result of message instability, the alternative transcripts contain pre-mature stop codons and do not accumulate to high levels, thus decreasing the level of AtGRP7 protein [79]. atSRp30 has similar effects on its own transcripts [65]. Another possible answer is that some splicing-related genes might regulate the splicing of other splicing-related genes. For example, overexpression of AtGRP7 and atSRp30 is known to affect the splicing of AtGRP8 and atSR1, respectively [65,79]. A third possibility is that the environment could affect the alternative splicing pattern. A good example is the SR1 gene. The ratio of two transcripts from the SR1 gene (SR1B/SR1) increases in a temperature-dependent manner [67]. Generally, heat or cold stress could cause intron retention in some splicing regulators, which could further alter the splicing pattern of other genes. The fourth possible regulators are intronless genes. Combining all these possibilities, a pathway to regulate splicing could be inferred as follows: environmental changes → splicing pattern changes in some specific splicing-related genes and/or intronless genes → expression pattern changes (including splicing pattern changes) in general splicing related genes → changes in splicing patterns for specific genes.

Conclusions
A large number of Arabidopsis splicing-related genes were computationally identified in this study by means of sequence comparisons and motif searches, including a tentative U4atac snRNA gene containing all conserved motifs, a new SR protein-coding gene (atRSp32) belonging to the atRSp31 family, and several genes related to genes encoding known splicing-related proteins (atULrp and atFCA2). A web-accessible database containing all the Arabidopsis splicing related genes has been constructed and will be expanded to other organisms in the near future. This compilation should provide a good foundation to study the splicing process in more detail and to determine to what extent these genes are conserved across the entire plant kingdom. Our data show that about 50% of the splicing-related genes are duplicated in Arabidopsis. The duplication ratios for splicing regulators are even higher, indicating that the splicing mechanism is generally conserved among plants, but that the regulation of splicing may be more variable and flexible, thus enabling plants to respond to their specific environments.

Materials and methods
Search for Arabidopsis snRNAs
Sequences of the 15 experimentally identified major snRNAs were downloaded from GenBank. The two minor snRNAs sequences were compiled from the literature [28]. These genes were used to search against the Arabidopsis genome at the AtGDB BLAST server [88] and at the SALK T-DNA Express web server [89]. Our initial analysis was based on Release 3.0 of the Arabidopsis genome (GenBank accession numbers NC_003070.4, NC_003071.3, NC_003074.4, NC_003075.3, and NC_003076.4). Local BLAST [21] was used to derive the locations of the snRNA homologs from more recently sequenced regions of the genome. Criteria used for local BLAST were 'e 1 -F F -W 7' (cutoff eval is 1, dust filter on, with a minimum word size of 7). Human and maize snRNAs were also included as query sequences, and all hits with e-values less than 10-5 were regarded as possible homologs. A total of 70 major snRNAs and three minor snRNAs were identified by this method. Each major snRNA type has 10-18 copies in the genome. A tentative gene name and gene model were assigned to each snRNA gene after comparison with the snRNAs identified in MATDB [90]. Sequence-similarity values were based on BLAST alignments.

Search for Arabidopsis splicing-related proteins
A three-round BLAST search strategy was used to identify Arabidopsis splicing related protein-coding genes. First, sequences of splicing-related proteins from human and Drosophila were downloaded from GenBank according to several recent proteomic studies [15-18] and the website compilation of Stephen Mount's group available at [91]. Human hnRNP proteins identified in a recent review [76] were downloaded from GenBank. All these sequences were used as queries in a local BLAST search against Arabidopsis annotated proteins (obtained from TIGR at [92]). All hits with an e-value less than 10-10 were collected as candidates. Many of these candidates had highly significant e-values (usually 10-30 or below and much lower than other hits). These candidates were regarded as true homologs.

In the second step, all identified true homologs were used to query the Arabidopsis protein set again. An e-value of 10-20 was used as a cutoff value to find possible paralogs of the true homologs. Sequences identified in both rounds of BLAST hits were regarded as main candidates for splicing related proteins.

Finally, the main candidates were queried against GenPept and all annotated human proteins (obtained from Ensembl [93]). All candidates with significant similarity to proteins unrelated to splicing were removed from the main candidate list, and all candidates with significant similarity to proteins related to splicing were regarded as true splicing-related genes and were promoted to the status of true homologs. The remaining candidates were regarded as unclassified splicing-related proteins. BLAST results were initially analyzed by MuSeqBox [94]). Two custom scripts were written to read MuSeqBox output files, largely automating the search procedure.

Gene structure and chromosomal locations
The gene structure and chromosomal locations for the genes encoding splicing-related proteins were retrieved from AtGDB [95]. The chromosomal locations of the snRNA genes were inferred from the BLAST results. The location maps (Figure 1) were generated using the AtGDB advanced search function [96]. Spliced alignments of ESTs and cDNAs generated by GeneSeqer [97] were used to verify gene models. Gene structure information was used as an important criteria to group homologs into gene families.

Protein domains
InterProScan 3.3 was downloaded from [98] and was subsequently used to search protein domain databases using default parameters [99]. A Perl script was written to process the text results from InterProScan. Protein domain information was used in comparisons of homologs from different species. The search of the National Center for Biotechnology Information Conserved Domain Database (NCBI-CDD) [100] was conducted manually for certain genes to confirm the InterPro results.

Duplication source
The gene families with multiple copies were inspected to determine whether they were likely to have derived from chromosome-duplication events. Gene models of the duplicated gene were searched against the gene list of each chromosome redundancy region at MATDB [101]. If the gene and its duplicate were both in the list, they were regarded as a chromosome duplication pair. Otherwise, they were assumed to be produced by random gene duplication.

Identification of alternative splicing
All Arabidopsis ESTs and cDNAs were aligned against the genome using the spliced alignment program GeneSeqer as made available through AtGDB [102]. We retrieved the intron and exon coordinates of the reliable cognate alignments from the database. Scripts were written to identify introns that overlap with other introns or exons. We defined the alternative splicing cases as follows: alternative donor (AltD): an intron has the same 3'-end coordinate but different 5'-end coordinate as another overlapping intron; alternative acceptor (AltA): an intron has the same 5'-end coordinate but different 3'-end coordinate as another intron; alternative position (AltP): an intron has different 5'-end and 3'-end coordinates as another overlapping intron; exon skipping (ExonS): an annotated intron completely contains an alternatively identified exon in the same transcription direction; intron retention (IntronR): an annotated intron is completely contained by an alternatively identified exon.

Database and interface construction
Details about each splicing-related gene were saved in a MySQL database. PHP scripts were written to interact with the database and generate the interface web pages. Text and BLAST searches were implemented by Perl-cgi scripts.

Acknowledgements
We thank Shannon Schlueter for help with the web page and database design and implementation. We are also grateful to Shailesh Lal, Carolyn Lawrence and Michael Sparks for discussions and critical reading of the manuscript and to the anonymous reviewers for excellent suggestions. This work was supported in part by a grant from the ISU Plant Sciences Institute and NSF grants DBI-0110189 and DBI-0110254 to V.B.

Figures and Tables
Figure 1 Sequence alignments of U4atac and U6atac snRNAs. The tentative Arabidopsis U4atac snRNA was aligned against the human U4atac snRNA (U62822) using CLUSTAL W [22]. Possible sequence domains are indicated by different background colors, with cyan indicating transcription signals (USE, upstream sequence element; TATA, TATA box), green indicating the region involved in the stem-loop-stem structure, and pink indicating the domain that binds Sm proteins. The corresponding interaction region in U6atac snRNA is also marked in green. Red background indicates G-T base-pairs in the stem-loop structure. Grey letters indicate the genome sequence upstream and downstream of the putative U4atac gene. Asterisks (upper panel) and black shading (lower panel) show conserved positions in the alignment.

Figure 2 Chromosomal locations of Arabidopsis snRNAs. Chromosomes 1 to 5 are represented to scale by the long thick lines in dark green. The small bars above the chromosomes indicate the presence of an snRNA gene in that region. Different colors represent different snRNA types: red, U1 snRNA; magenta, U2 snRNA; blue, U4 snRNA; green, U5 snRNA; yellow, U6 snRNA; black, minor snRNA. The seven U1-U4 snRNA gene clusters (red-blue) and the single U2-U5 snRNA gene cluster (magenta-green) are indicated by red circles.

Figure 3 Phylogenetic tree of the SC35 protein family. The phylogenetic tree was constructed on the basis of protein sequence alignments of the SC35 homologs in human, Drosophila, Arabidopsis and rice. The GenBank accession numbers for the sequences are as follows: hsSC35, Q01130; hsSRrp40, AAL57514; hsSRrp35: AAL57515; dmSC35, AAF53192; atSC35, NP_851261; atSR33/SCL33, NP_564685; atSCL30a, NP_187966; atSCL30, NP_567021; atSCL28, NP_197382; osSC35a, BAC79909; osSC35b, BAD09319; osSR33-1, AAP46199; osSCL30a/SR33-2, BAC799901; osSCL30-2, BAD19168. The sequences were aligned using CLUSTALW [22] with default parameters, and the phylogenetic tree was produced according to the neighbor-joining method using PAM substitution model distances as implemented in the PHYLIP package [103].

Table 1 Arabidopsis snRNA genes

Gene	GeneID	Chromosome	Strand	From	To	Length (nucleotides)	e-value	Similarity	GenBank ID	
atU1a*	At5g49054	5	-	19903323	19903158	166	1E-89	1-166, 100%	gi17660	
atU1-2	At4g23415	4	+	12225621	12225786	166	1E-58	1-166, 92%	gi22293582	
atU1-3	At5g51675	5	+	21013986	21014149	164	4E-55	3-166, 91%		
atU1-4	At5g25774	5	-	8972971	8972807	165	2E-51	1-166, 90%	gi22293583	
atU1-5	At1g08115	1	-	2538238	2538073	166	1E-46	1-166, 89%	gi22293581	
atU1-6	At3g05695	3	+	1681815	1681977	163	4E-40	4-166, 87%		
atU1-7	At3g05672	3	+	1657766	1657928	163	4E-40	4-166, 87%	gi22293580	
atU1-8	At5g27764	5	+	9832576	9832740	165	1E-39	1-166, 87%		
atU1-9	At5g26694	5	-	9494594	9494430	165	1E-27	1-166, 84%		
atU1-10	At1g11884	1	-	4007396	4007236	161	1E-18	4-61, 93%; 80-166, 88%		
atU1-11p	At4g16645	4	+	9370786	9370841	56	7E-17	4-59, 94%		
atU1-12p	At4g23565	4	-	12298871	12298802	70	1E-15	94-163, 90%		
atU1-13p	At5g49524	5	-	20112431	20112275	157	2E-14	4-50, 91%; 91-166, 88%		
atU1-14p	At1g35354	1	+	12986822	12986908	87	1E-06	10-60, 88%; 84-118, 88%		
atU2-1	At1g16825	1	+	5758381	5758575	195	2E-88	1-196, 96%		
atU2.2*	At3g57645	3	+	21357718	21357913	196	1E-107	1-196, 100%	gi17661	
atU2.3	At3g57765	3	-	21408595	21408400	196	1E-95	1-196, 97%	gi17662	
atU2.4	At3g56825	3	-	21052994	21052800	195	5E-86	1-196, 95%	gi17663	
atU2.5	At5g09585	5	+	2975013	2975208	196	7E-79	1-196, 93%	gi17664	
atU2.6	At3g56705	3	+	21015472	21015667	196	1E-83	1-196, 94%	gi17665	
atU2.7	At5g61455	5	-	24730829	24730634	196	5E-86	1-196, 95%	gi17666	
atU2-8	At5g67555	5	+	26966884	26967079	196	5E-86	1-196, 95%		
atU2.9	At4g01885	4	+	815273	815466	194	2E-82	1-194, 94%	gi17667	
atU2-10	At2g02938	2	+	849777	849972	196	3E-93	1-196, 96%	gi22293586	
atU2-10b/12	At2g02940	2	+	852859	853054	196	3E-93	1-196, 96%		
atU2-11	At1g09805/09895	1	-	3180736	3180547	190	8E-85	1-190, 95%		
atU2-13	At2g20405	2	+	8809169	8809364	196	3E-81	1-196, 94%	gi22293584	
atU2-14	At1g14165	1	+	4842274	4842469	196	3E-81	1-196, 94%	gi22293585	
atU2-15	At5g62415	5	+	25083790	25083985	196	4E-74	1-196, 92%		
atU2-16	At5g57835	5	-	23448717	23448522	196	2E-67	1-196, 92%		
atU2-17	At5g14545	5	-	4690105	4690008	98	3E-44	1-98, 97%		
atU2-18p	At3g26815	3	+	9881236	9881303	68	2E-14	1-68, 89%		
atU4.1*	At5g49056	5	-	19902970	19902817	154	4E-80	1-154, 99%	gi17673	
atU4.2	At3g06900	3	-	2178343	2178190	154	2E-75	1-154, 98%	gi17674	
atU4.3p	At5g49526	5	-	20112072	20112030	43	2E-11	15-57, 95%	gi17675	
atU4-4	At1g49242/49235	1	-	18222354	18222201	154	2E-75	1-154, 98%	gi22293588	
atU4-5	At5g25776	5	-	8972618	8972465	154	1E-70	1-154, 96%		
atU4-6	At1g11886	1	-	4007020	4006867	154	1E-70	1-154, 96%	gi22293587	
atU4-7	At5g27766	5	+	9832934	9833083	150	7E-66	1-150, 96%		
atU4-8	At5g26996	5	-	9494230	9494081	150	7E-66	1-150, 96%		
atU4-9	At1g79965	1	+	30086031	30086168	138	9E-47	18-154, 92%		
atU4-10	At1g35356	1	+	12987189	12987313	125	3E-34	1-124, 90%		
atU4-11p	At1g68395	1	+	25647322	25647396	75	9E-07	18-37, 100%; 60-102, 90%		
atU5.1*	At3g55865	3	-	20740607	20740503	105	6E-35	1-105, 94%	gi17676	
atU5.1b	At3g55855	3	-	20736881	20736780	102	7E-38	1-102, 96%	gi22293592	
atU5-2	At1g65115	1	+	24194482	24194586	105	1E-39	1-105, 96%		
atU5-3	At1g70185	1	+	26433396	26433497	102	7E-38	1-102, 96%	gi22293590	
atU5-4	At3g55645	3	+	20653843	20653947	105	3E-37	1-105, 95%		
atU5-5	At1g24105/24095	1	-	8525204	8525103	102	2E-35	1-102, 95%	gi22293591	
atU5-6	At1g04475	1	-	1215831	1215730	102	2E-35	1-102, 95%	gi22293589	
atU5-7	At4g02535	4	-	1114629	1114528	102	1E-30	2-103, 93%		
atU5-8	At3g25445	3	-	9227212	9227116	97	1E-20	5-101, 89%		
atU5-9	At1g79545	1	-	29928543	29928447	97	1E-20	5-101, 89%		
atU5-10	At5g14547	5	-	4690412	4690370	43	3E-12	24-67, 97%		
atU5-11	At5g54065	5	-	21957066	21957023	44	2E-10	20-64, 95%		
atU5-12	At1g71355	1	+	26895255	26895298	44	2E-10	20-64, 95%		
atU5-13	At5g53745	5	-	21829988	21829943	46	3E-09	24-70, 93%		
atU6.1*	At3g14735	3	+	4951596	4951697	102	1E-51	1-102, 100%	gi16516	
atU6.26	At3g13855	3	+	4561111	4561212	102	2E-49	1-102, 99%	gi16517	
atU6.29	At5g46315	5	+	18804616	18804717	102	2E-49	1-102, 99%	gi16518	
atU6-2	At5g62995	5	+	25296825	25296926	102	1E-51	1-102, 100%		
atU6-3	At4g27595	4	+	13782215	13782316	102	1E-51	1-102, 100%		
atU6-4	At4g03375	4	-	1483121	1483020	102	1E-51	1-102, 100%		
atU6-5	At4g33085	4	-	15965258	15965158	101	8E-37	1-101, 94%		
atU6-6	At4g35225	4	+	16754836	16754931	96	1E-32	1-102, 93%		
atU6-7	At2g15532	2	+	6784793	6784869	77	7E-25	4-80, 93%		
atU6-8p	At1g52605	1	+	19596398	19596476	96	2E-19	4-99, 87%		
atU6-9p	At1g53465	1	-	19960538	19960485	54	9E-09	21-74, 88%		
atU6-10p	At3g45705	3	+	16792802	16792888	87	2E-06	1-46, 89%; 62-100, 89%		
atU6-11p	At5g11085	5	-	3522167	3522143	25	9E-06	1-25, 100%		
atU12*	At1g61275	1	+	22606785	22606960	176	1E-95	1-176, 100%	†gi22293600	
atU6atac*	At5g40395	5	-	16183534	16183413	122	1E-63	1-122, 100%	†	
atU6atac-2	At1g21395	1	-	7491489	7491378	112	5E-20	1-65, 95%; 81-110, 93%		
atU4atac	At4g16065	4	+	9096374	9096532	159	N/A	N/A		
Chromosomal locations were determined by conducting BLAST searches against the Arabidopsis genome (Release 5.0). *The gene used for query in the BLAST search; †atU12 and atU6atac sequences, which were experimentally identified [28]. Their sequences were compiled manually from the cited paper. The GenBank gi numbers for the chromosome sequences used are as follows: chromosome 1, 42592260; chromosome 2, 30698031; chromosome 3, 30698537; chromosome 4, 30698542; chromosome 5, 30698605.

Table 2 Arabidopsis splicing-related proteins

Human homologs	Saccharomyces cerevisiae	Gene name	GeneID	Chromosome	Tnb	AltS	Chromosomal duplication	Protein domain	Reference	
1.1 Sm core proteins										
SmB	SmB1	atSmB-a	At5g44500	5	7		>4-5a	Sm, 1		
		atSmB-b	At4g20440	4	21	IntronR (1);	>4-5a	Sm, 1		
			
SmD1	SmD1	atSmD1-a	At3g07590	3	7	IntronR (1);		Sm, 1		
		atSmD1-b	At4g02840	4	13			Sm, 1		
			
SmD2	SmD2	atSmD2-a	At2g47640	2	7	AltA (1); AltD (1);		Sm, 1		
		atSmD2-b	At3g62840	3	25	AltA (1);		Sm, 1		
			
SmD3	SmD3	atSmD3-a	At1g76300	1	9		>1-1c	Sm, 1		
		atSmD3-b	At1g20580	1	7		>1-1c	Sm, 1		
			
SmE	SmE	atSmE-a	At4g30330	4	2		>2-4b	Sm, 1		
		atSmE-b	At2g18740	2	10	AltA (1);	>2-4b	Sm, 1		
			
SmF	SmF	atSmF	At4g30220	4	6			Sm, 1		
			
SmG	SmG	atSmG-a	At2g23930	2	13			Sm, 1		
		atSmG-b	At3g11500	3	9			Sm, 1		
			
LSM2	LSm2	atLSM2	At1g03330	1	7			Sm, 1		
			
LSM3	LSm3	atLSM3a	At1g21190	1	6		>1-1c	Sm, 1		
		atLSM3b	At1g76860	1	16		>1-1c	Sm, 1		
			
LSM4	LSm4	atLSM4	At5g27720	5	13			Sm, 1		
			
LSM5	LSm5	atLSM5 /SAD1	At5g48870	5	7	AltA (1);		Sm, 1	[47]	
			
LSM6	LSm6	atLSM6a	At3g59810	3	7		>2-3	Sm, 1		
		atLSM6b	At2g43810	2	5		>2-3	Sm, 1		
			
LSM7	LSm7	atLSM7	At2g03870	2	6			Sm, 1		
LSM8	LSm8	atLSM8	At1g65700	1	9			Sm, 1		
			
LSM1	LSm1	atLSM1a	At1g19120	1	8			Sm, 1		
		atLSM1b	At3g14080	3	9	IntronR (1);		Sm, 1		
	
										
1.2 U1 snRNP specific proteins										
U1A Subunit	Mud1	atU1A	At2g47580	2	14	ExonS (1);		RRM, 2	[49]	
			
U1C Subunit	Yhc1	atU1C	At4g03120	4	5			C2H2, 1; mrCtermi, 3		
			
U1-70K	Snp1	atU1-70K	At3g50670	3	32	IntronR (1);		RRM, 1	[48]	
			
-	Prp39	atPrp39a	At1g04080	1	12	ExonS (6);		HAT, 7; TPR-like, 1		
		atPrp39b	At5g46400	5	1			HAT, 4;		
			
FBP11	Prp40	atPrp40a	At1g44910	1	10	IntronR (1);		WW, 2; FF, 5		
FBP11	Prp40	atPrp40b	At3g19670	3	5			WW, 2; FF, 5		
			
Luc7-like protein	Luc7	atLuc7a	At3g03340	3	6			DUF259, 1		
		atLuc7b	At5g17440	5	8			DUF259, 1		
			
Related to Luc7-like protein	Luc7	atLuc7-rl	At5g51410	5	7	IntronR (1);		DUF259, 1		
	
										
1.3 17S U2 snRNP specific proteins										
U2A' Subunit	Lea1p	atU2A	At1g09760	1	21			LRR 4;		
			
U2B" Subunit	Msl1p	atU2B"a	At1g06960	1	6	AltD (1);	>1-2a	RRM, 2		
		atU2B"b	At2g30260	2	13	AltA (1); IntronR (1);	>1-2a	RRM, 2;		
			
SF3a120/SAP114 Subunit	Prp21p	atSAP114-1a	At1g14650	1	17	AltB (1);		SWAP/Surp, 2; Ubiquitin, 1		
		atSAP114-1b	At1g14640	1				SWAP/Surp, 2		
		atSAP114-2	At5g06520	5				SWAP/Surp, 4		
		atSAP114-3	At4g16200	4	1			SWAP/Surp, 3		
		atSAP114p	At4g15580	4				SWAP/Surp, 3; Ubiquitin, 1		
			
SF3a60/SAP61 Subunit	Prp9p	atSAP61	At5g06160	5	10	AltD (1);		C2H2, 1		
			
SF3a66/SAP62 Subunit	Prp11p	atSAP62	At2g32600	2	13			C2H2, 1;		
			
SF3b120/SAP130 Subunit	Rse1p	atSAP130a	At3g55200	3	6			CPSF_A, 1; WD40-like, 1	[50]	
		atSAP130b	At3g55220	3	7			CPSF_A, 1; WD40-like, 1	[50]	
			
SF3b150/SAP145 Subunit	Cus1p	atSF3b150	At4g21660	4	16			PSP, 1; DUF382, 1		
		atSF3b150p	At1g11520	1						
			
SF3b160/SAP155 Subunit	Hsh155	atSAP155	At5g64270	5	11			HEAT, 1; ARM, 2; SAP_155, 1		
			
SF3b53/SAP49 Subunit	Hsh49p	atSAP49a	At2g18510	2	20			RRM, 2		
		atSAP49b	At2g14550	2				RRM, 2		
			
p14	Snu17p	atP14-1	At5g12190	5	7			RRM, 1;		
		atP14-2	At2g14870	2				RRM, 1;		
			
SF3b 14b /PHP5A	Rds3p	atSF3b_14b-a	At1g07170	1	10		>1-2a	UPF0123, 1;		
		atSF3b_14b -b	At2g30000	2	8		>1-2a	UPF0123, 1;		
			
SF3b 10		SF3b10a	At4g14342	4	11			SF3b10, 1;		
		SF3b10b	At3g23325	3	6			SF3b10, 1;		
	
										
1.4 U5 snRNP specific proteins										
15 kD Subunit	Dib1p	atU5-15	At5g08290	5	28			DIM1, 1; Thioredoxin_2; 1		
			
40 kD Subunit		atU5-40	At2g43770	2	21			WD-40, 7;		
			
100 kD Subunit	Prp28p	atU5-100KD	At2g33730	2	13			DEAD, 1; Helicase_C, 1		
			
102 KD/Prp6-like	Prp6p	atU5-102KD	At4g03430	4	18			Ubiquitin, 1; TPR, 3; HAT, 15; TPR-like, 2; Prp1_N, 1		
			
116 kD Subunit /elongation	Snu114p	atU5-116-1a	At1g06220	1	19	ExonS (1);		EFG_C, 1; GTP_EFTU, 1; GTP_EFTU_D2; 1; Small_GTP, 1; EFG_IV, 1;		
		atU5-116-1b	At5g25230	5				EFG_C, 1; GTP_EFTU, 1; GTP_EFTU_D2; 1; EFG_IV, 1;		
		atU5-116-2	At1g56070	1	214			EFG_C, 1; GTP_EFTU, 1; GTP_EFTU_D2; 1; EFG_IV, 1;		
		atU5-116-3	At3g22980	3	3			EFG_C, 1; GTP_EFTU, 1; Small_GTP, 1;		
			
200 kD Subunit/Helicase	Brr2p	atU5-200-1	At5g61140	5	11	IntronR (1);		DEAD, 2; Helicase_C, 2; Sec63, 2; ARM, 1		
		atU5-200-2a	At1g20960	1	23			DEAD, 2; Helicase_C, 2; Sec63, 2		
		atU5-200-2b	At2g42270	2	5			DEAD, 2; Helicase_C, 2; Sec63, 2		
		atU5-200-3	At3g27730	3				DEAD, 1; Sec63, 1; RuvA domain 2-like, 1		
			
220 kD Subunit	Prp8p	atU5-220/Prp8a	At1g80070	1	33			Mov34, 1		
		atU5-220/Prp8b	At4g38780	4	2			Mov34, 1		
	
										
1.5 U4/U6 snRNP specific proteins										
U4/U6-90K / SAP90	Prp3p	atSAP90-1	At1g28060	1	10					
		atSAP90-2	At3g55930	3						
		atSAP90-3	At3g56790	3						
			
U4/U6-60K / SAP60	Prp4p	atSAP60	At2g41500	2	8			WD-40, 7; SFM, 1; WD40-like, 1		
			
U4/U6-20K / CYP20		atTri-20	At2g38730	2	11			Pro_isomerase, 1		
			
U4/U6-61KD	Prp31	atU5-61/Prp31a	At1g60170	1	26			Nop, 1		
		atU5-61/Prp31b	At3g60610	3				Nop, 1		
			
U4/U6-15.5K	Snu13p	atU4/U6-15.5a	At5g20160	5	18	IntronR (2);		Ribosomal_L7Ae, 1		
		atU4/U6-15.5b	At4g12600	4	14			Ribosomal_L7Ae, 1		
		atU4/U6-15.5c	At4g22380	4	9			Ribosomal_L7Ae, 1		
	
										
1.6 Tri-snRNP specific proteins										
Tri-65 KD	Snu66p	atTri65a	At4g22350	4	7			UCH; 1; ZnF_UBP, 1		
		atTri65b	At4g22290	4	20			UCH; 1; ZnF_UBP, 1; Pentaxin, 1		
		atTri65c	At4g22410	4				UCH; 1; ZnF_UBP, 1		
			
Tri-110 KD	SAD1	atTri110	At5g16780	5	7			SART-1, 1		
			
Tri-27 kD/RY1		atTri-27 kD/RY1	At5g57370	5	14					
			
hSnu23/FLJ31121	Snu23p	atSnu23	At3g05760	3	7			ZnF_U1, 1;		
	
										
1.7 18S U11/U2 snRNP specific proteins										
U11/U12-35K		atU11/U12-35kD	At2g43370	2	7	IntronR (1);		RRM, 1		
			
U11/U12-25K (-99 protein)		atU11/U12-25K	At3g07860	3	6	IntronR (2);		C2H2, 1;		
			
U11/U12-65K		atU11/U12-65K	At1g09230	1	15	AltA (1);		RRM, 2;PHOSPHOPANTETHEINE, 2;		
			
U11/U12-31K (MADP1)		atU11/U12-31K	At3g10400	3	5			RRM, 1;CCHC, 1;		
	
										
2.1 Splice site selection										
U2AF35		atU2AF35a/AUSa	At1g27650	1	26			RRM, 1; CCCH, 2;		
		atU2AF35/AUSb	At5g42820	5	8			RRM, 1; CCCH, 2;	[58]	
			
U2AF65	Mud2	atU2AF65b/AULa	At1g60900	1	10			RRM, 3;	[58]	
		atU2AF65a/AULb	At4g36690	4	29	AltA (1); IntronR (2);		RRM, 2;	[58]	
		atULrp	At2g33440	2	2			RRM, 1		
		AUL3p	At1g60830	1						
			
U2AF35 related protein		atUrp	At1g10320	1				RRM, 1; CCCH, 2;		
			
SF1/BBP		atSF1/BBP	At5g51300	5	23	IntronR (1);		RRM, 1; CCHC, 2; KH, 1;		
			
CBP20	Cbc1	atCBP20	At5g44200	5	8			RRM, 1	[56]	
			
CBP80	Cbc2p	atCBP80	At2g13540	2	21			MIF4G, 1; ARM, 3	[56]	
			
PTB/hnRNP I		atPTB1	At1g43190	1	26			RRM, 4;		
		atPTB2a	At3g01150	3	21	AltD (1); ExonS (1);	RRM, 2			
		atPTB2b	At5g53180	5	17	ExonS (1);		RRM, 2		
	
										
2.2 SR proteins										
SC35		atSC35	At5g64200	5	32	AltD (1);		RRM, 1;	[61]	
			
SRrp40/TASR-2		atSR33/atSCL33	At1g55310	1	12	IntronR (1);	>1-3b	RRM, 1	[63]	
		atSCL30a	At3g13570	3	32	ExonS (2); IntronR (4);	>1-3b	RRM, 1	[61]	
		atSCL30	At3g55460	3	14	ExonS (1);		RRM, 1	[61]	
		atSCL28	At5g18810	5	5			RRM, 1	[61]	
			
SF2/ASF		atSR1/atSRp34	At1g02840	1	37	AltA (1); IntronR (1);	>1-4	RRM, 2	[64,67]	
		atSRp34a	At4g02430	4	13	AltA (1); ExonS (1); IntronR (4);	>1-4	RRM, 2		
		atSRp34b	At3g49430	3	3	ExonS (1); IntronR (1);	RRM, 2			
		atSRp30	At1g09140	1	15	AltA (1);		RRM, 2	[65]	
			
9G8		atRSZp22/atSRZ22	At4g31580	4	26		>2-4e	RRM, 1; CCHC, 1	[63,66]	
		atRSZp22a	At2g24590	2	7		>2-4e	RRM, 1; CCHC, 1	[63,66]	
		atRSzp21/atSRZ21	At1g23860	1	18			RRM, 1; CCHC, 1	[63,66]	
		atRSZ33	At2g37340	2	30	IntronR (1);	>2-3	RRM, 1; CCHC, 2	[61]	
		atRSZ34	At3g53500	3	36	AltA (1); IntronR (3);	>2-3	RRM, 1; CCHC, 2	[61]	
			
-		atRSp32	At2g46610	2	23	AltD (1); IntronR (1);	>2-3	RRM, 2		
		atRSp31	At3g61860	3	17	AltA (1);	>2-3	RRM, 2	[59]	
		atRSp41	At5g52040	5	34	AltA (1);	>4-5b	RRM, 2	[59]	
		atRSp40/atRSP35	At4g25500	4	15	ExonS (1); IntronR (1);	>4-5b	RRM, 2	[59]	
	
										
2.3 17S U2 associated proteins										
hPrp43	Prp43p	atPrp43-1	At5g14900	5				HA2, 1		
		atPrp43-2a	At3g62310	3	17	AltA (1);	>2-3	DEAD, 1; Helicase_C, 1; HA2, 1		
		atPrp43-2b	At2g47250	2	14		>2-3	DEAD, 1; Helicase_C, 1; HA2, 1		
			
SR140		atSR140-1	At5g25060	5	11			Surp, 1;RRM, 1;, 1;RPR, 1;		
		atSR140-2	At5g10800	5	2			Surp, 1;RRM, 1;RPR, 1;		
			
SPF45		atSPF45	At1g30480	1	9			D111/G-patch domain, 1; RRM, 1;		
			
SPF30		atSPF30	At2g02570	2	9	AltA (1);		Tudor, 1;		
	
										
2.4 35S U5 associated proteins										
hPrp19*	Prp19p	atPrp19a	At1g04510	1	18		>1-2a	WD-40, 7; Ubox, 1;		
		atPrp19b	At2g33340	2	27	IntronR (1);	>1-2a	WD-40, 7; Ubox, 1;		
			
CDC5*	Cef1	atCDC5	At1g09770	1	12			SANT, 2;	[104]	
			
PRL1*	Prp46p	atPRL1	At4g15900	4	14			WD-40, 2;WD40like, 1;		
		atPRL2	At3g16650	3	6			WD-40, 2;WD40like, 1;		
			
AD-002*	Cwc15p	atAD-002	At3g13200	3	22			Cwf_Cwc_15, 1;		
			
HSP73/HSPA8*		HSP73-1	At3g12580	3	35			Hsp70, 1;		
		HSP73-2	At5g42020	5	51	IntronR (1);		Hsp70, 1;		
		HSP73-3	At5g02500	5	553	IntronR (1);		Hsp70, 1;		
			
SPF27/BCAS2*		atSPF27	At3g18165	3	15			BCAS2, 1;		
			
beta catenin-like 1*		atCTNNNBL1	At3g02710	3	12			Armadillo, 1;ARM, 1;		
			
hSyf1	Syf1p	atSyf1	At5g28740	5	7			TPR, 1;HAT, 10;TPRlike, 3;		
			
hSyf3/CRN	Syf3	atCRN1a	At5g45990	5				TPR, 1; HAT, 14; TPR-like, 2		
		atCRN1b	At3g13210	3				TPR, 1; HAT, 12; TPR-like, 2		
		atCRN1c	At5g41770	5	13			TPR, 1; HAT, 14; TPR-like, 2		
		atCRN2	At3g51110	3	8			TPR, 1; HAT, 9; TPR-like, 1		
			
hIsy1	Isy1p	atlsy1	At3g18790	3	10			Isy1, 1;		
			
GCIP p29	Syf2	atGCIPp29	At2g16860	2	12					
			
SKIP	Prp45p	atSKIP	At1g77180	1	28			SKIP/SNW, 1;		
			
hECM2	Ecm2p	atECM2-1a	At1g07360	1	21		>1-2a	RRM, 1;CCCH, 1;		
		atECM2-1b	At2g29580	2	10		>1-2a	RRM, 1;CCCH, 1;		
		atECM2-2	At5g07060	5				CCCH, 1;		
			
KIAA0560		atAquarius	At2g38770	2	11					
			
MGC23918		atMGC23918	At3g05070	3	7					
			
G10	Cwc14p	atG10	At4g21110	4	12			G10, 1;		
			
Cyp E		atCypE1a/CYP2	At2g21130	2	4		>2-4c	Pro_isomerase, 1		
		atCypE1b	At4g38740	4	59		>2-4c	Pro_isomerase, 1;		
		atCypE2a/ROC3	At2g16600	2	39		>2-4a	Pro_isomerase, 1		
		atCypE2b	At4g34870	4	80		>2-4a	Pro_isomerase, 1;		
			
PPIase-like 1		atPPIase-like1	At2g36130	2	10			Pro_isomerase, 1;		
	
										
2.5 Proteins specific for BΔU1										
NPW38		atNPW38	At2g41020	2	16	AltD (1); IntronR (1);		WW, 2;		
			
N-CoR1		atN-CoR1	At3g52250	3	3			SANT, 2;Homeodomain_like, 2;		
			
hPrp4 kinase		atPRP4K-1	At3g25840	3	13	ExonS (1);		Pkinase, 1;TyrKc, 1;S_Tkc, 1;, 1;Kinase_like, 1;		
		atPRP4K-2	At1g13350	1	5	IntronR (1);		Pkinase, 1;TyrKc, 1;S_Tkc, 1;, 1;Kinase_like, 1;		
		atPRP4K-3	At3g53640	3				Pkinase, 1;TyrKc, 1;S_Tkc, 1;, 1;Kinase_like, 1;		
			
FBP-21		atFBP21	At1g49590	1	12	ExonS (1); IntronR (3);	C2H2, 1;			
			
TBL1-rp 1		atTBL1-rp1	At5g67320	5	14			WD-40, 5;Peptidase_S9A_N, 1;LisH, 1;WD40like, 1;		
			
Smc-1		atSmc1	At3g54670	3	12			ATP_GTP_A_BS, 1;SMC_N, 1;SMC_C, 1;ABC_transporter, 1;SMC_hinge, 1;		
	
										
2.6 Exon junction complex (EJC) proteins										
ALY	Yra1p	atALY-1a	At5g02530	5	19	IntronR (1);		RRM, 1;		
		atALY-1b	At5g59950	5	16	IntronR (1);		RRM, 1;		
		atALY-2a	At5g37720	5	17		>1-5b	RRM, 1;		
		atALY-2b	At1g66260	1	38	ExonS (1);	>1-5b	RRM, 1;		
			
Y14		atY14	At1g51510	1	10	IntronR (1);		RRM, 1;RBM8, 4;		
			
Srm160-like		atSRM102	At2g29210	2	18	AltA (1);		PWI, 1		
			
Magoh		atMagoh	At1g02140	1	19			Mago_nashi, 1;		
			
Nuk-34/eIF4A3/DDX48		atDDX48/eIF4A3-1	At3g19760	3	50		>1-3a	DEAD, 1;Helicase_C, 1;		
		atDDX48/eIF4A3-2	At1g51380	1	5		>1-3a	DEAD, 1;Helicase_C, 1;		
			
RNPS1		atSR45/atRNPS1	At1g16610	1	27	AltA (1);		RRM, 1	[63]	
			
UAP56		atUAP56a	At5g11200	5	21	AltA (1);		DEAD, 1; Helicase_C, 1		
		atUAP56b	At5g11170	5	25			DEAD, 1; Helicase_C, 1		
			
pinin		atPinin	At1g15200	1	9	AltA (1);		Pinin/SDK/memA, 1;		
	
										
2.7 Second step splicing factors										
Prp22	Prp22	atPrp22-1	At3g26560	3	11			DEAD, 1; Helicase_C, 1; S1, 1; HA2, 1;		
		atPrp22-2	At1g26370	1	5			DEAD, 1; Helicase_C, 1; HA2, 1		
		atPrp22-3	At1g27900	1	15			DEAD, 1; Helicase_C, 1; HA2, 1		
			
Prp17	Prp17p	atPrp17-1	At1g10580	1	10			WD-40, 7;		
		atPrp17-2	At5g54520	5	5	AltA (1);		WD-40, 6;		
			
Prp18	Prp18	atPrp18-1	At1g03140	1	16			Prp18, 1; SFM 1;		
		atPrp18-2	At1g54590	1				Prp18, 1		
			
Slu7	Slu7p	atSLU7-1a	At1g65660	1	6					
		atSLU7-1b	At4g37120	4	11					
		atSLU7-2	At3g45950	3						
			
Prp16	Prp16p	atPrp16	At5g13010	5	22			DEAD, 1; Helicase_C, 1; HA2, 1		
	
										
2.8 Other known splicing factors										
SRm300		atSRM300like	At3g23900	3	5	AltD (1);		RRM, 1; Filamin/ABP280 repeat, 1		
			
hTra-2/SFRS10		atTra/SFRS1	At1g07350	1	25	ExonS (1); IntronR (3);		RRM, 1		
			
Prp2		atPrp2-1a	At1g32490	1	9		>1-2c	DEAD, 1; Helicase_C, 1; HA2, 1		
		atPrp2-1b	At2g35340	2			>1-2c	DEAD, 1; Helicase_C, 1; HA2, 1		
		atPrp2-2	At4g16680	4				DEAD, 1; Helicase_C, 1; HA2, 1		
			
Prp5		atPrp5-1a	At3g09620	3				DEAD, 1; Helicase_C, 1		
		atPrp5-1b	At1g20920	1	11			DEAD, 1; Helicase_C, 1		
		atPrp5-2	At2g47330	2	9			DEAD, 1; Helicase_C, 1		
			
hDbr1	dbr1	atDbr1	At4g31770	4	12			Metallophos, 1; DBR1, 1		
	
										
3.1 SR protein kinase										
Lammer/CLK kinase		AFC1	At3g53570	3	11	AltA (1); IntronR (3);		PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1	[74]	
		AFC2	At4g24740	4	9	ExonS (1);		PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1	[74]	
		AFC3	At4g32660	4	9	AltD (1); IntronR (1);		PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1	[74]	
			
SRPK1		atSRPK1a	At2g17530	2	7		>2-4a	PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1		
		atSRPK1b	At4g35500	4	10		>2-4a	PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1		
			
SRPK2		atSRPK2a	At5g22840	5	2			PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1		
		atSRPK2b	At3g53030	3	7			PKinase, 1; TyrKc, 1; S_Tkc, 1; PKinase-like, 1		
		atSRPK2c	At3g44850	3	1			PKinase, 1; TyrKc, 1; PKinase-like, 1		
	
										
3.2 Glycine-rich RNA binding protein										
HnRNP A/B		atGRBP1a	At1g18630	1	5		>1-1c	RRM, 1		
		atGRBP1b	At1g74230	1	14		>1-1c	RRM, 1; Eggshell, 4		
		atGRBP1c	At4g13850	4	17		>3-4	RRM, 1		
		atGRBP1d	At3g23830	3	12	AltA (1);	>3-4	RRM, 1		
		atGRBP1e	At5g61030	5	8			RRM, 1; PfkB_Kinase, 1		
		atGRBP2	At2g16260	2				RRM, 1		
		AtGRP7/atGRBP3a	At2g21660	2	182	AltD (1); IntronR (3);	>2-4c	RRM, 1	[77]	
		AtGRP8/atGRBP3b	At4g39260	4	67	AltB (1); AltD (1); IntronR (5);	>2-4c	RRM, 1	[77]	
	
										
3.3 hnRNP A/B family										
hnRNP A/B		AtRNPA/B_1	At4g14300	4	4			RRM, 2	[11]	
		AtRNPA/B_2	At2g33410	2	13			RRM, 2	[11]	
		AtRNPA/B_3	At5g55550	5	13	IntronR (3);	>4-5c	RRM, 2	[11]	
		AtRNPA/B_4	At4g26650	4	21		>4-5c	RRM, 2	[11]	
		AtRNPA/B_5	At5g47620	5	12	AltD (2);		RRM, 2	[11]	
		AtRNPA/B_6	At3g07810	3	18	AltA (1);		RRM, 2; FKBP_PPIASE_2, 2	[11]	
		AtRNPA/B_7	At1g58470	1	6			RRM, 2		
		AtRNPA/B_8a	At5g40490	5	3			RRM, 2; Eggshell, 4		
		AtRNPA/B_8b	At1g17640	1				RRM, 2		
		AtRNP_N1	At3g13224	3	16	IntronR (1);		RRM, 2;HUDSXLRNA, 2;		
		UBA2a	At3g56860	3	23	IntronR (1);	>2-3	RRM, 2	[78]	
		UBA2b	At2g41060	2	9		>2-3	RRM, 2	[78]	
		UBA2c	At3g15010	3	10	IntronR (1);		RRM, 2	[78]	
	
										
3.4 Other hnRNPs (with animal homologs)										
hnRNP E1/E2		at-hnRNP-E	At3g04610	3	10			KH, 3;		
			
hnRNP F/ hnRNP H		at-hnRNP-F/ AtRNPH/F_1	At5g66010	5	9	AltA (1);		RRM, 2	[11]	
		at-hnRNP-H/ AtRNPH/F_2	At3g20890	3				RRM, 2	[11]	
			
hnRNP G		at-hnRNP-G1	At5g04280	5	6			RRM, 1; CCHC, 1		
		at-hnRNP-G2	At3g26420	3	35	AltA (1);		RRM, 1; CCHC, 1		
		at-hnRNP-G3	At1g60650	1	7			RRM, 1; CCHC, 1		
			
hnRNP P2		at-hnRNP-P	At1g50300	1	9			RRM, 1; ZnF_RBZ, 2		
			
hnRNP R/Q		hnRNP-R1	At4g00830	4	19			RRM, 3;		
		hnRNP-R2	At3g52660	3	1		>2-3	RRM, 3;		
		hnRNP-R3 / AtRNPA/B_9	At2g44710	2	13		>2-3	RRM, 3		
			
CUG-BP		AtCUG-BP1	At4g03110	4	4	AltA (1); IntronR (1);		RRM, 3; HUDSXLRNA, 4	[11]	
		AtCUG-BP2	At1g03457	1	9	AltA (1);		RRM, 3; HUDSXLRNA, 4	[11]	
			
(CUG-BP)		atFCA1	At4g16280	4	13	AltB (1); IntronR (1);		RRM, 2; WW, 1	[81]	
		atFCA2	At2g47310	2	6			RRM, 2; WW, 1		
	
										
3.5 Other plant hnRNPs										
		AtUBP1a	At1g54080	1	48	AltA (1);	>1-3b	RRM, 3	[84]	
		AtUBP1c	At3g14100	3	13		>1-3b	RRM, 3	[84]	
		AtUBP1b	At1g17370	1	17			RRM, 3	[84]	
			
		UBA1a	At2g22090	2	15		>2-4c	RRM, 1	[78]	
		UBA1b	At2g22100	2	2		>2-4c	RRM, 1	[78]	
		UBA1c	At2g19380	2	1			RRM, 1; C2H2, 3	[78]	
			
		atRBP45a	At5g54900	5	42		>4-5c	RRM, 3	[85]	
		atRBP45c	At4g27000	4	52		>4-5c	RRM, 3	[85]	
		AtRBP45b	At1g11650	1	53			RRM, 3	[85]	
		atRBP45d	At5g19350	5	10			RRM, 3		
			
		AtRBP47a	At1g49600	1	10		>1-3a	RRM, 3	[85]	
		AtRBP47b	At3g19130	3	21		>1-3a	RRM, 3	[85]	
		AtRBP47c	At1g47490	1	23	IntronR (1);		RRM, 3	[85]	
		AtRBP47c'	At1g47500	1	12			RRM, 3	[85]	
			
		Ath1	At4g16830	4	34			HANP4_PAI-RBP1, 1	[105]	
		Ath2	At4g17520	4	29		>4-5a	HANP4_PAI-RBP1, 1	[105]	
		Ath3	At5g47210	5	67	IntronR (1);	>4-5a	HANP4_PAI-RBP1, 1		
Gene names were kept consistent with names used in previous publications or derived from the names of the respective homologs (yeast names are given in the S.c. column, where available). The Tnb column gives the numbers of cognate cDNAs and ESTs supporting the gene structure. The AltS column indicates evidence for alternative splicing, including alternative donor site (AltD), alternative acceptor site (AltA), alternative position (AltP, both acceptor and donor sites are different), exon skipping (ExonS), and intron retention (IntronR). Chromosomal duplication indicates a known chromosome duplication region. Functional groups of proteins are separated by long lines spanning all columns. Different members in the group are separated by short lines starting at the Arabidopsis gene name. Genes duplicated in Arabidopsis are clustered together with no line between them. Dash line separate the Prp19 complex from other 35S U5 associated proteins and * indicates proteins in that complex. Abbreviations for domains are as follows: ABC_transporter: ABC transporter; Armadillo: Armadillo; ARM: ARM repeat fold; ATP_GTP_A_BS: ATP/GTP-binding site motif A (P-loop); BCAS2: Breast carcinoma amplified sequence 2; C2H2: Zn-finger, C2H2 matrin type; C2H2: Zn-finger, C2H2 type; CCCH: Zn-finger, C-x8-C-x5-C-x3-H type; CCHC: Zn-finger, CCHC type; CPSF_A: CPSF A subunit, C-terminal; Cwf_Cwc_15: Cwf15/Cwc15 cell cycle control protein; DBR1: Lariat debranching enzyme, C-terminal; DEAD: ATP-dependent helicase, DEAD-box; DEAD: DEAD/DEAH box helicase; DIM1: Pre-mRNA splicing protein; DUF259: Protein of unknown function DUF259; DUF382: Protein of unknown function DUF382; EFG_C: Elongation factor G, C-terminal; EFG_IV: Elongation factor G, domain IV; Eggshell: Eggshell protein; FF: FF domain; FKBP_PPIASE_2: Peptidylprolyl isomerase, FKBP-type; G10: G10 protein; GTP_EFTU_D2: Elongation factor Tu, domain 2; GTP_EFTU: Protein synthesis factor, GTP-binding; HA2: Helicase-associated region; HANP4_PAI-RBP1, 1: Hyaluronan/mRNA binding protein; HAT: RNA-processing protein, HAT helix; Helicase_C: Helicase, C-terminal; Homeodomain_like: Homeodomain-like; Hsp70: Heat shock protein Hsp70; HUDSXLRNA: Paraneoplastic encephalomyelitis antigen; Isy1: Isy1-like splicing; Kinase_like: Protein kinase-like; LisH: Lissencephaly type-1-like homology motif; LRR: Leucine-rich repeat,; Mago_nashi: Mago nashi protein; Metalloph: Metallophosphoesterase; MIF4G: Initiation factor eIF-4 gamma, middle; Mov34: Mov34/MPN/PAD-1; mrCtermi: Molluscan rhodopsin C-terminal tail; Nop: Pre-mRNA processing ribonucleoprotein, binding region; Peptidase_S9A_N: Peptidase S9A, prolyl oligopeptidase, N-terminal beta-propeller domain; PfkB_Kinase: Carbohydrate kinase, PfkB; PHOSPHOPANTETHEINE: Phosphopantetheine attachment site; Pinin/SDK/memA: Pinin/SDK/memA protein; Pkinase: Protein kinase; Pro_isomerase: Peptidyl-prolyl cis-trans isomerase, cyclophilin type; Prp18: Prp18 domain; Prp1_N: PRP1 splicing factor, N-terminal; PSP: PSP, proline-rich; PWI: Splicing factor PWI; RBM8: RNA binding motif protein 8; Ribosomal_L7Ae: Ribosomal protein L7Ae/L30e/S12e/Gadd45; RPR: Regulation of nuclear pre-mRNA protein; RRM: RNA-binding region RNP-1 (RNA recognition motif); S1: RNA binding S1; SANT: Myb DNA-binding domain; SAP_155: Splicing factor 3B subunit_1; SART-1: SART-1 protein; Sec63: Sec63 domain; SF3b10: Splicing factor 3B subunit 10; SFM: Splicing factor motif; SKIP/SNW: SKIP/SNW domain; Small_GTP: Small GTP-binding protein domain; SMC_C: Structural maintenance of chromosome protein SMC, C-terminal; SMC_hinge: SMCs flexible hinge; SMC_N: SMC protein, N-terminal; Sm_like_riboprot: Small nuclear-like ribonucleoprotein; Sm: Small nuclear ribonucleoprotein (Sm protein); S_Tkc: Serine/threonine protein kinase; Surp: SWAP/Surp; Thioredoxin_2: Thioredoxin domain 2; TPRlike: TPR-like; TPR: TPR repeat; Tudor: Tudor domain; TyrKc: Tyrosine protein kinase; Ubox: Zn-finger, modified RING; UCH: Peptidase C19, ubiquitin carboxyl-terminal hydrolase family 2; UPF0123: Protein of unknown function UPF0123; UPF0123: Protein of unknown function UPF0123; WD-40: G-protein beta WD-40 repeat; WD40like: WD40-like; WW: WW/Rsp5/WWP domain; ZnF_RBZ: Zn-finger, Ran-binding; ZnF_U1: Zn-finger, U1-like; ZnF_UBP: Zn-finger in ubiquitin thiolesterase.

Table 3 Duplication source involving Arabidopsis splicing-related proteins

	Genes	Family*	Single/multiple	Duplication ratio	Duplication events	Chromosomal duplications	Chromosomal duplication ratio	
snRNP proteins	91	54	27/27	50.0%	37	7	18.9%	
Splicing factors	109	58	33/25	43.1%	51	14	27.5%	
Splicing regulator	60	18	4/14	77.8%	42	11	26.2%	
Total	260	130	64/66	50.8%	130	32	24.6%	
*Family indicates both single copy gene and multiple-copy gene families. The Chromosomal duplication ratio column gives the fraction of all duplication events caused by chromosomal duplications.

Table 4 Alternative splicing in splicing-related genes

	Genes	AltA	AtlD	AltP	ExonS	IntronR	Overall	Ratio	
snRNP proteins	91	6	3	1	3	11	22	23.2%	
Splicing factors	109	14	5	0	11	21	38	34.9%	
Splicing regulator	60	8	4	2	1	12	20	33.3%	
Total	260	28	12	3	15	44	80	30.8%	
The column entries are the numbers of genes in which the respective alternative splicing events can occur. AltA, alternative acceptor site; AltD, alternative donor site; AltP, alternative intron position (both acceptor and donor sites are different); ExonS, exon skipping; IntronR, intron retention. The Overall and Ratio columns give the number and fraction of genes with any type of alternative splicing, respectively.
==== Refs
Kazan K  Alternative splicing and proteome diversity in plants: the tip of the iceberg has just emerged. Trends Plant Sci 2003 8 468 471 14557042 10.1016/j.tplants.2003.09.001 
Lorkovic ZJ Wieczorek Kirk DA Lambermon MH Filipowicz W  Pre-mRNA splicing in higher plants. Trends Plant Sci 2000 5 160 167 10740297 10.1016/S1360-1385(00)01595-8 
Reddy ASN  Nuclear pre-mRNA splicing in plants. Critical Rev Plant Sci 2001 20 523 571 10.1016/S0735-2689(01)80004-6 
Mount SM Salz HK  Pre-messenger RNA processing factors in the Drosophila genome. J Cell Biol 2000 150 F37 F44 10908584 10.1083/jcb.150.2.F37 
Käufer NF Potashkin J  Analysis of the splicing machinery in fission yeast: a comparison with budding yeast and mammals. Nucleic Acids Res 2000 28 3003 3010 10931913 10.1093/nar/28.16.3003 
Arabidopsis Genome Initiative  Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 2000 408 796 815 11130711 10.1038/35048692 
Barakat A Szick-Miranda K Chang IF Guyot R Blanc G Cooke R Delseny M Bailey-Serres J  The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome. Plant Physiol 2001 127 398 415 11598216 10.1104/pp.127.2.398 
Beisson F Koo AJ Ruuska S Schwender J Pollard M Thelen JJ Paddock T Salas JJ Savage L Milcamps A  Arabidopsis genes involved in acyl lipid metabolism. A 2003 census of the candidates, a study of the distribution of expressed sequence tags in organs, and a web-based database. Plant Physiol 2003 132 681 697 12805597 10.1104/pp.103.022988 
Wang D Harper JF Gribskov M  Systematic trans-genomic comparison of protein kinases between Arabidopsis and Saccharomyces cerevisiae. Plant Physiol 2003 132 2152 2165 12913170 10.1104/pp.103.021485 
Aubourg S Kreis M Lecharny A  The DEAD box RNA helicase family in Arabidopsis thaliana. Nucleic Acids Res 1999 27 628 636 9862990 10.1093/nar/27.2.628 
Lorkovic ZJ Barta A  Genome analysis: RNA recognition motif (RRM) and K homology (KH) domain RNA-binding proteins from the flowering plant Arabidopsis thaliana. Nucleic Acids Res 2002 30 623 635 11809873 10.1093/nar/30.3.623 
TAIR: gene family information 
Makarova OV Makarov EM Urlaub H Will CL Gentzel M Wilm M Lührmann R  A subset of human 35S U5 proteins, including Prp19, function prior to catalytic step 1 of splicing. EMBO J 2004 23 2381 2391 15175653 10.1038/sj.emboj.7600241 
Will CL Schneider C Hossbach M Urlaub H Rauhut R Elbashir S Tuschl T Lührmann R  The human 18S U11/U12 snRNP contains a set of novel proteins not found in the U2-dependent spliceosome. RNA 2004 10 929 941 15146077 10.1261/rna.7320604 
Zhou Z Sim J Griffith J Reed R  Purification and electron microscopic visualization of functional human spliceosomes. Proc Natl Acad Sci USA 2002 99 12203 12207 12215496 10.1073/pnas.182427099 
Makarov EM Makarova OV Urlaub H Gentzel M Will CL Wilm M Lührmann R  Small nuclear ribonucleoprotein remodeling during catalytic activation of the spliceosome. Science 2002 298 2205 2208 12411573 10.1126/science.1077783 
Rappsilber J Ryder U Lamond AI Mann M  Large-scale proteomic analysis of the human spliceosome. Genome Res 2002 12 1231 1245 12176931 10.1101/gr.473902 
Jurica MS Moore MJ  Pre-mRNA splicing: awash in a sea of proteins. Mol Cell 2003 12 5 14 12887888 10.1016/S1097-2765(03)00270-3 
Arabidopsis Splicing Related Genes Database 
Arabidopsis thaliana Genome Database 
Altschul SF Madden TL Schäffer AA Zhang J Zhang Z Miller W Lipman DJ  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997 25 3389 3402 9254694 10.1093/nar/25.17.3389 
Thompson JD Higgins DG Gibson TJ  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994 22 4673 4680 7984417 
Vankan P Filipowicz W  Structure of U2 snRNA genes of Arabidopsis thaliana and their expression in electroporated plant protoplasts. EMBO J 1988 7 791 799 
Vankan P Edoh D Filipowicz W  Structure and expression of the U5 snRNA gene of Arabidopsis thaliana. Conserved upstream sequence elements in plant U-RNA genes. Nucleic Acids Res 1988 16 10425 10440 3205713 
Vankan P Filipowicz W  A U-snRNA gene-specific upstream element and a -30 'TATA box' are required for transcription of the U2 snRNA gene of Arabidopsis thaliana. EMBO J 1989 8 3875 3882 2583119 
Waibel F Filipowicz W  U6 snRNA genes of Arabidopsis are transcribed by RNA polymerase III but contain the same two upstream promoter elements as RNA polymerase II-transcribed U-snRNA genes. Nucleic Acids Res 1990 18 3451 3458 2362802 
Hofmann CJ Marshallsay C Waibel F Filipowicz W  Characterization of the genes encoding U4 small nuclear RNAs in Arabidopsis thaliana. Mol Biol Rep 1992 17 21 28 1287476 
Shukla GC Padgett RA  Conservation of functional features of U6atac and U12 snRNAs between vertebrates and higher plants. RNA 1999 5 525 538 10199569 10.1017/S1355838299982213 
Marker C Zemann A Terhorst T Kiefmann M Kastenmayer JP Green P Bachellerie JP Brosius J Huttenhofer A  Experimental RNomics: identification of 140 candidates for small non-messenger RNAs in the plant Arabidopsis thaliana. Curr Biol 2002 12 2002 2013 12477388 10.1016/S0960-9822(02)01304-0 
ASRG snRNAs 
Patel AA Steitz JA  Splicing double: insights from the second spliceosome. Nat Rev Mol Cell Biol 2003 4 960 970 14685174 10.1038/nrm1259 
Connelly S Filipowicz W  Activity of chimeric U small nuclear RNA (snRNA)/mRNA genes in transfected protoplasts of Nicotiana plumbaginifolia: U snRNA 3'-end formation and transcription initiation can occur independently in plants. Mol Cell Biol 1993 13 6403 6415 8413239 
Connelly S Marshallsay C Leader D Brown JW Filipowicz W  Small nuclear RNA genes transcribed by either RNA polymerase II or RNA polymerase III in monocot plants share three promoter elements and use a strategy to regulate gene expression different from that used by their dicot plant counterparts. Mol Cell Biol 1994 14 5910 5919 8065324 
Tarn WY Steitz JA  Pre-mRNA splicing: the discovery of a new spliceosome doubles the challenge. Trends Biochem Sci 1997 22 132 137 9149533 10.1016/S0968-0004(97)01018-9 
Shukla GC Padgett RA  U4 small nuclear RNA can function in both the major and minor spliceosomes. Proc Natl Acad Sci USA 2004 101 93 98 14691257 10.1073/pnas.0304919101 
Shukla GC Cole AJ Dietrich RC Padgett RA  Domains of human U4atac snRNA required for U12-dependent splicing in vivo. Nucleic Acids Res 2002 30 4650 4657 12409455 10.1093/nar/gkf609 
Krämer A  The structure and function of proteins involved in mammalian pre-mRNA splicing. Annu Rev Biochem 1996 65 367 409 8811184 10.1146/annurev.bi.65.070196.002055 
Will CL Lührmann R  Protein functions in pre-mRNA splicing. Curr Opin Cell Biol 1997 9 320 328 9159080 10.1016/S0955-0674(97)80003-8 
ASRG proteins 
Stevens SW Abelson J  Purification of the yeast U4/U6.U5 small nuclear ribonucleoprotein particle and identification of its proteins. Proc Natl Acad Sci USA 1999 96 7226 7231 10377396 10.1073/pnas.96.13.7226 
Stevens SW Barta I Ge HY Moore RE Young MK Lee TD Abelson J  Biochemical and genetic analyses of the U5, U6, and U4/U6 × U5 small nuclear ribonucleoproteins from Saccharomyces cerevisiae. RNA 2001 7 1543 1553 11720284 
Gottschalk A Neubauer G Banroques J Mann M Lührmann R Fabrizio P  Identification by mass spectrometry and functional analysis of novel proteins of the yeast [U4/U6.U5] tri-snRNP. EMBO J 1999 18 4535 4548 10449419 10.1093/emboj/18.16.4535 
Caspary F Shevchenko A Wilm M Seraphin B  Partial purification of the yeast U2 snRNP reveals a novel yeast pre-mRNA splicing factor required for pre-spliceosome assembly. EMBO J 1999 18 3463 3474 10369685 10.1093/emboj/18.12.3463 
Krämer A Grüter P Gröning K Kastner B  Combined biochemical and electron microscopic analyses reveal the architecture of the mammalian U2 snRNP. J Cell Biol 1999 145 1355 1368 10385517 10.1083/jcb.145.7.1355 
Fabrizio P Esser S Kastner B Lührmann R  Isolation of S. cerevisiae snRNPs: comparison of U1 and U4/U6.U5 to their human counterparts. Science 1994 264 261 265 8146658 
Will CL Lührmann R  Spliceosomal UsnRNP biogenesis, structure and function. Curr Opin Cell Biol 2001 13 290 301 11343899 10.1016/S0955-0674(00)00211-8 
Xiong L Gong Z Rock CD Subramanian S Guo Y Xu W Galbraith D Zhu JK  Modulation of abscisic acid signal transduction and biosynthesis by an Sm-like protein in Arabidopsis. Dev Cell 2001 1 771 781 11740939 10.1016/S1534-5807(01)00087-9 
Golovkin M Reddy AS  Structure and expression of a plant U1 snRNP 70K gene: alternative splicing of U1 snRNP 70K pre-mRNAs produces two different transcripts. Plant Cell 1996 8 1421 1435 8776903 10.1105/tpc.8.8.1421 
Simpson GG Clark GP Rothnie HM Boelens W van Venrooij W Brown JW  Molecular characterization of the spliceosomal proteins U1A and U2B' from higher plants. EMBO J 1995 14 4540 4550 7556097 
Casacuberta E Puigdomenech P Monofort A  A genomic duplication in Arabidopsis thaliana contains a sequence similar to the human gene coding for SAP130. Plant Physiol Biochem 2001 39 565 573 10.1016/S0981-9428(01)01280-3 
Golovkin M Reddy AS  Expression of U1 small nuclear ribonucleoprotein 70K antisense transcript using APETALA3 promoter suppresses the development of sepals and petals. Plant Physiol 2003 132 1884 1891 12913145 10.1104/pp.103.023192 
Gottschalk A Tang J Puig O Salgado J Neubauer G Colot HV Mann M Seraphin B Rosbash M Lührmann R Fabrizio P  A comprehensive biochemical and genetic analysis of the yeast U1 snRNP reveals five novel proteins. RNA 1998 4 374 393 9630245 
McLean MR Rymond BC  Yeast pre-mRNA splicing requires a pair of U1 snRNP-associated tetratricopeptide repeat proteins. Mol Cell Biol 1998 18 353 360 9418882 
Huang T Vilardell J Query CC  Pre-spliceosome formation in S. pombe requires a stable complex of SF1-U2AF(59)-U2AF(23). EMBO J 2002 21 5516 5526 12374752 10.1093/emboj/cdf555 
Lewis JD Gorlich D Mattaj IW  A yeast cap binding protein complex (yCBC) acts at an early step in pre-mRNA splicing. Nucleic Acids Res 1996 24 3332 3336 8811086 10.1093/nar/24.17.3332 
Kmieciak M Simpson CG Lewandowska D Brown JW Jarmolowski A  Cloning and characterization of two subunits of Arabidopsis thaliana nuclear cap-binding complex. Gene 2002 283 171 183 11867224 
Hugouvieux V Kwak JM Schroeder JI  An mRNA cap binding protein, ABH1, modulates early abscisic acid signal transduction in Arabidopsis. Cell 2001 106 477 487 11525733 10.1016/S0092-8674(01)00460-3 
Domon C Lorkovic ZJ Valcarcel J Filipowicz W  Multiple forms of the U2 small nuclear ribonucleoprotein auxiliary factor U2AF subunits expressed in higher plants. J Biol Chem 1998 273 34603 34610 9852132 10.1074/jbc.273.51.34603 
Lopato S Waigmann E Barta A  Characterization of a novel arginine/serine-rich splicing factor in Arabidopsis. Plant Cell 1996 8 2255 2264 8989882 10.1105/tpc.8.12.2255 
Lopato S Mayeda A Krainer AR Barta A  Pre-mRNA splicing in plants: characterization of Ser/Arg splicing factors. Proc Natl Acad Sci USA 1996 93 3074 3079 8610170 10.1073/pnas.93.7.3074 
Lopato S Forstner C Kalyna M Hilscher J Langhammer U Indrapichate K Lorkovic ZJ Barta A  Network of interactions of a novel plant-specific Arg/Ser-rich protein, atRSZ33, with atSC35-like splicing factors. J Biol Chem 2002 277 39989 39998 12176998 10.1074/jbc.M206455200 
Golovkin M Reddy AS  The plant U1 small nuclear ribonucleoprotein particle 70K protein interacts with two novel serine/arginine-rich proteins. Plant Cell 1998 10 1637 1648 9761791 10.1105/tpc.10.10.1637 
Golovkin M Reddy AS  An SC35-like protein and a novel serine/arginine-rich protein interact with Arabidopsis U1-70K protein. J Biol Chem 1999 274 36428 36438 10593939 10.1074/jbc.274.51.36428 
Lazar G Schaal T Maniatis T Goodman HM  Identification of a plant serine-arginine-rich protein similar to the mammalian splicing factor SF2/ASF. Proc Natl Acad Sci USA 1995 92 7672 7676 7644475 
Lopato S Kalyna M Dorner S Kobayashi R Krainer AR Barta A  atSRp30, one of two SF2/ASF-like proteins from Arabidopsis thaliana, regulates splicing of specific plant genes. Genes Dev 1999 13 987 1001 10215626 
Lopato S Gattoni R Fabini G Stevenin J Barta A  A novel family of plant splicing factors with a Zn knuckle motif: examination of RNA binding and splicing activities. Plant Mol Biol 1999 39 761 773 10350090 10.1023/A:1006129615846 
Lazar G Goodman HM  The Arabidopsis splicing factor SR1 is regulated by alternative splicing. Plant Mol Biol 2000 42 571 581 10809003 10.1023/A:1006394207479 
Tronchere H Wang J Fu XD  A protein related to splicing factor U2AF35 that interacts with U2AF65 and SR proteins in splicing of pre-mRNA. Nature 1997 388 397 400 9237760 10.1038/41137 
Lin CH Patton JG  Regulation of alternative 3' splice site selection by constitutive splicing factors. RNA 1995 1 234 245 7489496 
ASRG SR protein gene structure 
Cowper AE Caceres JF Mayeda A Screaton GR  Serine-arginine (SR) protein-like factors that antagonize authentic SR proteins and regulate alternative splicing. J Biol Chem 2001 276 48908 48914 11684676 10.1074/jbc.M103967200 
Chan SP Kao DI Tsai WY Cheng SC  The Prp19p-associated complex in spliceosome activation. Science 2003 302 279 282 12970570 10.1126/science.1086602 
Yong J Pellizzoni L Dreyfuss G  Sequence-specific interaction of U1 snRNA with the SMN complex. EMBO J 2002 21 1188 1196 11867547 10.1093/emboj/21.5.1188 
Bender J Fink GR  AFC1, a LAMMER kinase from Arabidopsis thaliana, activates STE12-dependent processes in yeast. Proc Natl Acad Sci USA 1994 91 12105 12109 7991592 
Savaldi-Goldstein S Aviv D Davydov O Fluhr R  Alternative splicing modulation by a LAMMER kinase impinges on developmental and transcriptome expression. Plant Cell 2003 15 926 938 12671088 10.1105/tpc.011056 
Krecic AM Swanson MS  hnRNP complexes: composition, structure, and function. Curr Opin Cell Biol 1999 11 363 371 10395553 10.1016/S0955-0674(99)80051-9 
Heintzen C Melzer S Fischer R Kappeler S Apel K Staiger D  A light- and temperature-entrained circadian clock controls expression of transcripts encoding nuclear proteins with homology to RNA-binding proteins in meristematic tissue. Plant J 1994 5 799 813 8054987 10.1046/j.1365-313X.1994.5060799.x 
Lambermon MH Fu Y Wieczorek Kirk DA Dupasquier M Filipowicz W Lorkovic ZJ  UBA1 and UBA2, two proteins that interact with UBP1, a multifunctional effector of pre-mRNA maturation in plants. Mol Cell Biol 2002 22 4346 4357 12024044 10.1128/MCB.22.12.4346-4357.2002 
Staiger D Zecca L Wieczorek Kirk DA Apel K Eckstein L  The circadian clock regulated RNA-binding protein AtGRP7 autoregulates its expression by influencing alternative splicing of its own pre-mRNA. Plant J 2003 33 361 371 12535349 10.1046/j.1365-313X.2003.01629.x 
Simpson GG Dijkwel PP Quesada V Henderson I Dean C  FY is an RNA 3' end-processing factor that interacts with FCA to control the Arabidopsis floral transition. Cell 2003 113 777 787 12809608 10.1016/S0092-8674(03)00425-2 
Macknight R Bancroft I Page T Lister C Schmidt R Love K Westphal L Murphy G Sherson S Cobbett C Dean C  FCA, a gene controlling flowering time in Arabidopsis, encodes a protein containing RNA-binding domains. Cell 1997 89 737 745 9182761 10.1016/S0092-8674(00)80256-1 
Quesada V Macknight R Dean C Simpson GG  Autoregulation of FCA pre-mRNA processing controls Arabidopsis flowering time. EMBO J 2003 22 3142 3152 12805228 10.1093/emboj/cdg305 
Paillard L Legagneux V Osborne HB  A functional deadenylation assay identifies human CUG-BP as a deadenylation factor. Biol Cell 2003 95 107 113 12799066 10.1016/S0248-4900(03)00010-8 
Lambermon MH Simpson GG Wieczorek Kirk DA Hemmings-Mieszczak M Klahre U Filipowicz W  UBP1, a novel hnRNP-like protein that functions at multiple steps of higher plant nuclear pre-mRNA maturation. EMBO J 2000 19 1638 1649 10747031 10.1093/emboj/19.7.1638 
Lorkovic ZJ Wieczorek Kirk DA Klahre U Hemmings-Mieszczak M Filipowicz W  RBP45 and RBP47, two oligouridylate-specific hnRNP-like proteins interacting with poly(A)+ RNA in nuclei of plant cells. RNA 2000 6 1610 1624 11105760 10.1017/S1355838200001163 
Riechmann JL Heard J Martin G Reuber L Jiang C Keddie J Adam L Pineda O Ratcliffe OJ Samaha RR  Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science 2000 290 2105 2110 11118137 10.1126/science.290.5499.2105 
Vision TJ Brown DG Tanksley SD  The origins of genomic duplications in Arabidopsis. Science 2000 290 2114 2117 11118139 10.1126/science.290.5499.2114 
AtGDB BLAST 
T-DNAexpress: the SIGnAL Arabidopsis gene mapping tool 
MIPS: MATDB snRNAs 
Drosophila mRNA processing factors 
TIGR ftp site 
Ensembl 
Xing L Brendel V  Multi-query sequence BLAST output examination with MuSeqBox. Bioinformatics 2001 17 744 745 11524378 10.1093/bioinformatics/17.8.744 
AtGDB 
AtGDB advanced search 
Brendel V Xing L Zhu W  Gene structure prediction from consensus spliced alignment of multiple ESTs matching the same genomic locus. Bioinformatics 2004 20 1157 1169 14764557 10.1093/bioinformatics/bth058 
InterPro 
Zdobnov EM Apweiler R  InterProScan - an integration platform for the signature-recognition methods in InterPro. Bioinformatics 2001 17 847 848 11590104 10.1093/bioinformatics/17.9.847 
NCBI-CDD search 
Arabidopsis thaliana: MATDB Redundancy Viewer 
Zhu W Schlueter SD Brendel V  Refined annotation of the Arabidopsis genome by complete expressed sequence tag mapping. Plant Physiol 2003 132 469 484 12805580 10.1104/pp.102.018101 
PHYLIP 
Hirayama T Shinozaki K  A cdc5+ homolog of a higher plant, Arabidopsis thaliana. Proc Natl Acad Sci USA 1996 93 13371 13376 8917598 10.1073/pnas.93.23.13371 
Landsberger M Lorkovic ZJ Oelmuller R  Molecular characterization of nucleus-localized RNA-binding proteins from higher plants. Plant Mol Biol 2002 48 413 421 11905967 10.1023/A:1014089531125

