Comprehensive analysis of mutually exclusive alternative splicing in C. elegans.

Mutually exclusive selection of one exon in a cluster of exons is a rare form of alternative pre-mRNA splicing, yet suggests strict regulation. However, the repertoires of regulation mechanisms for the mutually exclusive (ME) splicing in vivo are still unknown. Here, we experimentally explore putative ME exons in C. elegans to demonstrate that 29 ME exon clusters in 27 genes are actually selected in a mutually exclusive manner. Twenty-two of the clusters consist of homologous ME exons. Five clusters have too short intervening introns to be excised between the ME exons. Fidelity of ME splicing relies at least in part on nonsense-mediated mRNA decay for 14 clusters. These results thus characterize all the repertoires of ME splicing in this organism.

Introduction

Alternative processing of precursor mRNAs (pre-mRNAs) is a major source of protein diversity and plays crucial roles in development, differentiation, and diseases in higher eukaryotes.- One of the most elaborately regulated forms of alternative pre-mRNA processing is mutually exclusive (ME) alternative splicing, through which only one exon is mutually exclusively selected from a cluster of exons at a time to determine critical aspects of the target genes such as ligand-binding specificity of receptors and properties of enzymes and channels.- The ME exons occur only as pairs in vertebrates, but the number of ME exons in a cluster in invertebrates can be more than two in some genes. The extreme example is the Dscam gene of a fruit fly Drosophila melanogaster, which has four ME exon clusters containing 12, 48, 33, and two variants.

Several mechanisms have been proposed for the mutually exclusive nature of the ME exon selection.- Steric interference between two splice sites of an intervening intron between two ME exons has been proposed to make the ME exons physically incapable of being spliced to each other in some mammalian genes.,, This is due to a shorter distance between the 5′ splice site and the branch point on the intron than minimal spacing required for a spliceosome to be productively assembled. Another mechanism proposed to prohibit double inclusion and double skipping of the ME exons is spliceosome incompatibility., This is proposed for a tandem exon pair flanked by U2- and U12-type introns on either side of the exons. Disposal of aberrantly spliced mRNAs by a surveillance system termed nonsense-mediated mRNA decay (NMD) is considered to play substantial roles for some genes.- Antagonism of repression by base-paring interaction of a docker with one of the selector sequences is proposed for the Dscam exon 6 cluster.,

In C. elegans, it is estimated by a recent genome-wide analysis that up to 25% of the protein-coding genes undergo alternative pre-mRNA processing and 55 events were assigned to ME alternative splicing. In previous studies, we have elucidated tissue-specific and/or developmental selection patterns and regulation mechanisms for some of the ME exon clusters in C. elegans by generating fluorescence alternative splicing reporters and isolating splicing factor mutants. In the case of exons 5B/5A of the egl-15 gene, encoding fibroblast growth factor receptors (FGFRs),- the RBFOX family and SUP-12 cooperatively bind to the upstream flanking intron of the upstream exon to repress the upstream exon in muscles. In the case of exons 9/10 of the let-2 gene, encoding α2 subunit of collagen type IV, ASD-2 binds to the downstream flanking intron of the downstream exon to promote inclusion of the downstream exon in muscle-specific and developmentally regulated manners., In the case of exons 7a/7b of the unc-32 gene, encoding subunit a of V0 domain of vacuolar proton-translocating ATPase (V-ATPase), UNC-75 binds to the intervening intron to repress the downstream exon and the RBFOX family binds to the downstream intron to promote inclusion of the upstream exon in the nervous system. Thus, the tissue-specificity, trans-acting factors, positions of the cis-elements, and functions of the factors for the regulation of ME exon clusters, vary from gene to gene in C. elegans. These findings raise questions about to what extent the repertoires and regulation mechanisms for the ME exon clusters have evolved in this organism.

Here we explore all the 55 putative ME splicing events in C. elegans listed in Ramani et al. that utilized high-throughput sequencing and microarray profiling of polyA+ RNAs isolated from four and five different developmental stages, respectively. We experimentally test whether the putative ME exons are actually mutually exclusively selected by reverse transcription-polymerase chain reaction (RT-PCR). For the verified ME exons, we analyze the nucleotide and amino acid sequence identities of the ME exons in each cluster. To also elucidate to what extent the mutually exclusive nature of the exon selection rely on NMD, we compare the RT-PCR patterns between a wild-type strain N2 and an NMD-deficient mutant smg-2.,

Results and Discussion

Table 1 summarizes the results of the comprehensive RT-PCR analyses at the L1 stage. The 55 events were assigned to 41 clusters in 37 genes. Eight of the clusters were considered to be tandem cassette exon pairs rather than ME exons because we detected in-frame double-inclusion and/or double-skipping isoforms (> 5% of the sum in molar concentration) in addition to the single-inclusion isoforms (Table 1; Fig. 1A; and data not shown). Notably, these cassette exons are multiple of three (3n) nucleotides (nt) in length except for those carrying natural termination codons (Table 1). Two exons in two genes were considered to be single cassette exons and two other exons in two genes appear to be constitutively included. The other 29 clusters in 27 genes were considered to be mutually exclusive (Table 1) since the single-inclusion isoforms were detected in our experiments and/or in the literature and other isoforms were almost undetectable or degraded by NMD in the wild-type background (see below). We confirmed that the single-inclusion isoforms were also almost exclusively expressed at the young adult stage (data not shown).

Table 1. Summary of experimental validation of the putative ME exon clusters

		Name and ID in WS239				Position in WS235				Supporting Reads*a	Sequence Identity		Intervening Intron(s) [nt]*b	NMD-Dependence*c
AS Type	Chr	Gene	WBGene ID	Str	Exon	Left end	Right end	Exon Length			Nucleotide	Amino Acid
Homologous ME	I	F30F8.9	WBGene00009276	+	4a	7837844	7837926	83	3n +2	9	41.0%	17.9%	683	Yes
					4b	7838610	7838692	83	3n +2	12
		lev-11	WBGene00002978	-	7b	14623010	14623148	139	3n +1	1615	71.2%	78.3%	422	No
					7a*d	14623571	14623709	139	3n +1	83
					5b	14624030	14624147	118	3n +1	478	81.4%	74.4%	155, 123	No
					5c	14624303	14624420	118	3n +1	94	47.5%	35.9%
					5a	14624544	14624661	118	3n +1	706	83.9%	71.8%
	II	snt-1	WBGene00004921	+	6B*d	6815048	6815187	140	3n +2	0	64.3%	55.3%	116	No
					6A	6815304	6815437	134	3n +2	63
		lat-1	WBGene00002251	+	3a	8903388	8903494	107	3n +2	24	55.8%	45.9%	112	Partially
					3b	8903607	8903719	113	3n +2	72
	III	let-805	WBGene00002915	+	19a	2676547	2676664	118	3n +1	44	40.4%	33.3%	97	No
					19b	2676762	2676870	109	3n +1	34
		gly-6	WBGene00001631	-	8b	4325419	4325523	105	3n	26	44.1%	28.9%	11	-
					8a	4325535	4325642	108	3n	49
		pdfr-1	WBGene00015735	-	8b	6636402	6636556	155	3n +2	77	53.4%	30.8%	336	Partially
					8a	6636893	6637017	125	3n +2	24
		coq-2	WBGene00000762	+	6a	6937600	6937733	134	3n +2	13	47.0%	40.0%	678	Yes
					6b	6938412	6938545	134*e	3n +2	12
		unc-32	WBGene00006768	+	7a	8909233	8909355	123	3n	46	47.2%	35.4%	725	-
					7b	8910081	8910221	141	3n	56
		gly-5	WBGene00001630	+	9b	13200771	13200881	111	3n	16	77.7%	54.1%	34, 1466	-
					9c	13200916	13201020	105	3n	13	66.1%	54.1%
					9a	13202487	13202588	102	3n	6	71.4%	48.6%
		W06F12.2	WBGene00012306	-	2b	13730072	13730243	172	3n +1	80	35.5%	17.5%	129	Yes
					2a	13730373	13730526	154	3n +1	28
	IV	gba-3	WBGene00008706	+	6a	17478400	17478647	248	3n +2	78	75.4%	80.5%	207	Yes
					6b	17478855	17479102	248	3n +2	82
	V	unc-62	WBGene00006796	+	7a*d	4504975	4505120	146	3n +2	4	49.1%	30.2%	238	No
					7b	4505359	4505516	158	3n +2	162
		gck-1	WBGene00001526	+	4a	8952725	8952879	155	3n +2	19	52.2%	48.1%	138	Partially
					4b	8953018	8953175	158	3n +2	31
		akt-1	WBGene00000102	+	6a	10250841	10251032	192	3n	136	62.8%	63.8%	294	-
					6b	10251327	10251533	207	3n	41
		atn-1	WBGene00000228	+	4a	12480728	12480888	161	3n +2	147	61.4%	57.1%	225	Partially
					4b	12481114	12481196	83	3n +2	26
		slo-1	WBGene00004830	+	10a	18499582	18499694	113	3n +2	21	65.5%	64.1%	856	Partially
					10b	18500551	18500663	113	3n +2	11
	X	mrp-1	WBGene00003407	-	13D	579498	579653	156	3n	31*f	68.6%	53.8%	334, 154, 435	-
					13C	579988	580143	156	3n	35	76.9%	59.6%
					13B	580298	580453	156	3n	37	74.4%	53.8%
					13A	580889	581044	156	3n	2	69.2%	50.0%
		cca-1	WBGene00000367	-	18b	7854051	7854201	151	3n +1	8	53.0%	44.0%	97	Partially
					18a	7854299	7854428	130	3n +1	20
		bet-2	WBGene00010199	+	3a	10571708	10571827	120	3n	35	45.8%	35.0%	37	-
					3b	10571865	10571984	120	3n	32
		let-2	WBGene00002280	-	10	16386613	16386723	111	3n	164	65.8%	56.8%	30	-
					9	16386754	16386861	108	3n	104
Non-Homologous ME	I	pqn-72	WBGene00004154	-	6	3274267	3274812	546	3n	53	37.0%	21.6%	117	-
					5	3274930	3275010	81	3n	3
		tom-1	WBGene00006594	+	14a	5552880	5553369	490	3n +1	172	50.0%	16.4%	1371	Partially
					14b	5554741	5554864	124	3n +1	72
		ttx-7	WBGene00008765	-	5b	7301663	7301774	112	3n +1	16	41.8%	17.1%	132	Partially
					5a	7301907	7302027	121	3n +1	54
	III	unc-32	WBGene00006768	+	4a	8906524	8906678	155	3n +2	34	60.0%	17.3%	288, 237	Partially
					4b	8906967	8907073	107	3n +2	49	51.6%	18.0%
					4c	8907311	8907432	122	3n +2	26	49.7%	23.1%
	IV	fbl-1	WBGene00001403	+	5D	9542035	9542292	258	3n	57	35.0%	22.1%	435	-
					5C	9542728	9542868	141	3n	10
	V	del-6	WBGene00011891	+	5	10574706	10574733	28	3n +1*g	4	53%*h	60%*h	200	Yes
					6	10574934	10575098	165	3n	230
	X	egl-15	WBGene00001184	+	5B	11017580	11017930	351	3n	38	40.1%	10.3%	14	-
					5A	11017945	11018139	195	3n	5
Tandem CE	I	tom-1	WBGene00006594	+	17	5557960	5557971	12	3n	-	-	-	-	-
Tandem CE					18	5558453	5558479	27	3n	-	-	-	-	-
Tandem CE	I	lev-11	WBGene00002978	-	9b	14622304	14622511	208	3n +1*g	-	-	-	-	-
Tandem CE					9a	14622658	14622743	86	3n +2*g	-	-	-	-	-
Single CE	II	etr-1	WBGene00001340	+	4	166071	166150	80	3n +2	-	-	-	-	-
Tandem CE	II	zyg-12	WBGene00006997	+	8	4952490	4952531	42	3n*g	-	-	-	-	-
Tandem CE					9	4952647	4952694	48	3n	-	-	-	-	-
Tandem CE	II	C34F11.3	WBGene00016415	+	10	5205040	5205075	36	3n	-	-	-	-	-
Tandem CE					11	5205444	5205632	189	3n	-	-	-	-	-
Tandem CE	III	clp-1	WBGene00000542	+	3	7981673	7981810	138	3n	-	-	-	-	-
Tandem CE					4	7982464	7982538	75	3n	-	-	-	-	-
Constitutive	IV	unc-44	WBGene00006780	+	13	5984796	5984980	185	3n +2	-	-	-	-	-
Single CE					14	5985093	5985170	78	3n	-	-	-	-	-
Undetected					15	5985577	5986355	779	3n +2	-	-	-	-	-
Alternative Acceptors	IV	unc-43	WBGene00006779	-	14S	10329097	10329195	99*i	3n	-	-	-	-	-
Tandem CE	V	K10D6.2	WBGene00010742	-	5L	11194203	11194342	140*j	3n +2*g	-	-	-	-	-
Tandem CE					4	11194574	11194627	54	3n	-	-	-	-	-
Tandem CE	V	Y69H2.3	WBGene00013481	-	7	18661622	18661711	90	3n	-	-	-	-	-
Tandem CE					6	18662119	18662313	195	3n	-	-	-	-	-
Constitutive	X	F49E2.5	WBGene00009888	+	5	9554299	9554547	249	3n	-	-	-	-	-
Tandem CE					6	9554692	9554823	132	3n	-	-	-	-	-
Constitutive	X	T23E7.2	WBGene00020732	+	16	17679922	17680035	114	3n	-	-	-	-	-

ME, mutually exclusive exons; CE, cassette exon. *aTotal number of sequence reads mapped to each ME exon and its junctions with the flanking exons out of 10.3 million mapped reads in RNA-seq analysis of polyA+ RNA from synchronized smg-2 L1 larvae. *bIntervening introns shorter than 40 nt are underlined. *cYes, NMD isoform(s) were significantly increased in the smg-2 mutant (P < 0.05 in a modified χ2 test of the RT-PCR products). Partially, NMD isoform(s) were slightly increased in the smg-2 mutant (difference in the amount of the splice variant > 0.4 [% of sum] in molar concentration). No, no apparent difference in the RT-PCR patterns between N2 and the smg-2 mutant. *dThese exons were expressed at very low levels at the L1 stage. At the young adult stage, lev-11 exon 7a and snt-1 exon 6B were still rare, while unc-62 exon 7a was readily detected. *eOut-of-frame tandem acceptor sites were frequently used for coq-2 exon 6b. *fmrp-1 exon 13D is aberrantly annotated in the RefSeq model used for mapping. *gThese exons carry in-frame termination codons. *hNucleotide and amino acid sequence identities of del-6 exons 5/6 appear high because of their short lengths. *iIn-frame tandem acceptor sites were frequently used for unc-43 exon 14. *jTandem donor sites were frequently used for K10D6.2 exon 5.

Figure 1. RT-PCR analyses of the putative ME exons in the wild-type (N2) and the smg-2 (yb979) mutant. (A) clp-1 exons 3/4. (B) del-6 exons 5/6. del-6 exon 5 is unique in that it carries a natural termination codon. (C) bet-2 exons 3a/3b. Note that the intervening intron is retained instead of double ME exon inclusion for this cluster. (D) unc-32 exons 7a/7b. (E) akt-1 exons 6a/6b. A non-productive exon 6a isoform utilizing an aberrant acceptor site is detected in the smg-2 mutant. (F) fbl-1 exons 5D/5C. (G) F30F8.9 exons 4a/4b. (H) coq-2 exons 6a/6b. (I) gck-1 exons 4a/4b. (J and K) lev-11 exons 5a/5c/5b. (L) let-805 exons 19a/19b. Splicing patterns are schematically indicated. Coding regions are in orange. Arrows indicate predicted positions of undetected isoforms indicated on the right. Asterisks indicate non-specific bands.

Features of the ME exon clusters

The 29 ME exon clusters can be divided into two groups according to sequence similarity of the ME exons. Homologous ME exon clusters include 19 pairs, two trios, and one quad of homologous ME exons, while non-homologous clusters include six pairs and one trio of non-homologous ME exons (Table 1). The homologous clusters may be originated from exon duplication., The lengths of the homologous ME exons are close to or exactly the same as the counterpart(s) except for atn-1 exons 4a/4b, while those of the non-homologous ME exons are often far different from the counterpart(s). Nevertheless, reading frames in the downstream common exons are preserved whichever exon in the clusters is selected in almost all cases. The only exception is exons 5/6 of the del-6 gene, encoding a degenerin-like ion channel protein, where exon 5 consists of 3n+1 nt and carries a natural termination codon while exon 6 consists of 3n nt and has no termination codon (Fig. 1B).

Four of the ME exon clusters consist of more than two ME exons. Exons 13A/13B/13C/13D of the mrp-1 gene, encoding an ATP-binding cassette (ABC) transporter, is the only cluster with four ME exons. All the four exons are 156 (3n) nt in length and are homologous to each other (Table 1). Exons 9b/9c/9a of the gly-5 gene, encoding a UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferase, is one of the three clusters with three ME exons. These exons are homologous to each other with almost the same size of 3n nt. The intervening intron between exons 9b and 9c is just 34 nt (discussed later). Exons 4a/4b/4c of the unc-32 gene are the only non-homologous trio of the ME exons. The unc-32 gene has another cluster of ME exons 7a/7b (Table 1) and we have recently reported that these two clusters in the single gene are independently regulated in tissue-specific manners. Exons 5a/5c/5b of the lev-11 gene, encoding multiple tropomyosin isoforms, are homologous and of exactly the same size (3n+1 nt). Pre-mRNA processing of the lev-11 gene is complex due to the combination of tissue-specific promoters, clusters of ME exons 4a/4b, 5a/5b/5c, and 7a/7b, and tandem cassette exons 9a/9b (Table 1). The complex structures and pre-mRNA processing patterns of the tropomyosin genes are evolutionarily conserved in metazoans, suggesting functional significance of multiple tropomyosin isoforms.

Steric interference to prohibit double inclusion of ME exons

Size distribution of the overall introns suggests that the minimal size of the introns is ~40 nt in C. elegans. Among the ME exon clusters we have already reported, the intervening introns for exons 5B/5A of the egl-15 gene- and exons 9/10 of the let-2 gene, are 14 nt and 30 nt, respectively, and we have never observed mRNA isoforms where these short introns are excised. These observations are consistent with the idea that the short introns of less than 40 nt cannot be excised because of the steric interference like in mammals, although no strong consensus are found for the branch point in C. elegans. According to this criterion, three more clusters are considered to be physically incapable of double exon inclusion: exons 3a/3b of bet-2 encoding a BET (two bromodomains) family protein (Fig. 1C), gly-5 exons 9b/9c, and exons 8a/8b of gly-6, a paralog of gly-5 (Table 1). Notably, the lengths of all the ME exons in these five clusters are 3n nt.

As U11 or U12 snRNA or an AT/AC splice-junction are not found in C. elegans, we need not consider the spliceosome incompatibility in the regulation of the ME exons here.

NMD-dependence of the mutually exclusive selection

If the lengths of the ME exons are 3n nt, inclusion or skipping of the ME exons does not cause a frame-shift or a premature termination codon (PTC) in the mRNA isoforms. Consistent with this idea, there was no apparent difference in the amounts of multiple-inclusion and all-skipping isoforms between the wild-type and the smg-2 mutant for such clusters (Table 1; Fig. 1D–F; and data not shown).

If the lengths of the ME exons are not 3n nt, multiple inclusion and all-skipping of the ME exons cause frame-shifts to create aberrant termination codons and such mRNA isoforms should be eliminated by NMD. We found that multiple-inclusion and/or all-skipping isoforms are evidently or slightly more abundant in the smg-2 mutant than in the wild-type for 14 out of the 19 clusters where the ME exons are not 3n nt (Table 1; Fig. 1G–I; and data not shown), indicating that these non-productive mRNA isoforms are actually eliminated by NMD in the wild-type. For lev-11 exons 7a/7b, snt-1 exons 6B/6A, and unc-62 exons 7a/7b, we confirmed predominant use of only one of the two ME exons (Table 1) as in the literature and this can be a reason why aberrantly spliced isoforms are rare and undetectable for these clusters. For the other two clusters, lev-11 exons 5a/5c/5b and let-805 exons 19a/19b, the RT-PCR patterns were indistinguishable between the wild-type and the smg-2 mutant (Table 1; Fig. 1J–L). All these results indicate that the fidelity of the splicing regulation varies among the ME exon clusters and some of them rely on the mRNA surveillance system.

Statistics of the ME exons and flanking introns

Figure 2 summarizes the statistics of the 63 experimentally verified ME exons and their flanking and intervening introns.

Figure 2. Statistics of the ME exons and their flanking introns. (A and B) Size distributions of the 63 verified ME exons (A) and their 92 flanking introns, including the intervening introns (B). The mean and median sizes are also indicated. (C) Sequence logos of the splice acceptor and donors sites of the 63 ME exons.

The median size of the ME exons (134 nt) (Fig. 2A) is similar to those of the entire unique exons in confirmed genes (144 nt). Most of the ME exons (60 of 63) are shorter than 260 nt and the average size (151 nt) is close to the median (Fig. 2A). In contrast, the size distribution of the entire unique exons has a fatter tail, making the average of 201 nt. The three ME exons longer than 350 nt belong to distinct non-homologous ME exon clusters (Table 1). The shortest ME exon (28 nt) exceptionally carries a natural termination codon (Table 1; Fig. 1B). Therefore, the size of the 48 homologous ME exons are in a relatively narrow range (83–248 nt) for C. elegans.

The mean size of the introns flanking the ME exons, including the five short intervening introns discussed above, is 378 nt (Fig. 2B), substantially longer than the overall average of the introns (267 nt). The median size of the introns flanking the ME exons is 145.5 nt (Fig. 2B), whereas more than half of all the C. elegans introns are 100 nt or less and most of them are near the minimal length, indicating that the introns flanking the ME exons tend to be longer than constitutive introns. This is consistent with a previous finding that many of cis-elements regulating alternative splicing in C. elegans are found in introns. Eleven out of the 29 ME exon clusters have UGCAUG stretch(es) in the flanking introns and/or in the ME exons (data not shown), suggesting tissue-specific splicing regulation by the RBFOX family splicing factors ASD-1 and FOX-1. Six out of the 29 clusters are affected in the unc-75 mutant, suggesting neuron-specific splicing regulation.

Figure 2C summarizes the sequences of the splice acceptor and donor sites for the verified ME exons. These are more diversified from the consensus sequences of the acceptor site (TTTTCAG/R) and the donor site (AG/GTAAGTT) in C. elegans, where R stands for A or G. Furthermore, two (2.2%) of the 92 flanking introns, let-2 intron 10 and del-6 intron 6, start with GC, a weaker donor than GT,, although GC-AG introns are rare (0.373%) in C. elegans like in other eukaryotes. Therefore, the splice sites of the ME exons are considered to be weaker than those of constitutive exons, consistent with previous findings on alternative splice sites in higher organisms.

Table 2 summarizes gene ontology (GO) analysis of 25 genes with GO terms out of the 27 genes with the verified ME exon clusters. It indicates enrichment of genes encoding membrane or extracellular matrix proteins (P < 0.001, Fisher’s exact test).

Table 2. Gene ontology analysis of 25 genes with the ME exons and GO terms

Ontology type	GO term ID	Fold enrichment	Count in 25 Genes with ME Exons and GO terms	Count in all genes with GO Terms (12,834)	P value (Fisher's Exact Test)	Term
Biological_process	GO:0046928	257	2	4	2.18E-05	regulation of neurotransmitter secretion
	GO:0030163	16	4	131	1.11E-04	protein catabolic process
	GO:0007166	79	2	13	2.80E-04	cell surface receptor linked signal transduction
	GO:0040011	3.5	9	1327	5.77E-04	locomotion
	GO:0034765	49	2	21	7.48E-04	regulation of ion transmembrane transport
	GO:0043050	45	2	23	8.99E-04	pharyngeal pumping
Cellular_component	GO:0005865	114	2	9	1.30E-04	striated muscle thin filament
	GO:0016021	2.9	12	2143	2.77E-04	integral to membrane
	GO:0016020	3.1	11	1847	3.37E-04	membrane
	GO:0005604	57	2	18	5.47E-04	basement membrane
	GO:0005578	47	2	22	8.22E-04	proteinaceous extracellular matrix
Molecular_function	GO:0005201	257	2	4	2.18E-05	extracellular matrix structural constituent
	GO:0005244	49	2	21	7.48E-04	voltage-gated ion channel activity

Conclusion

We demonstrated that the 29 ME exon clusters in the 27 genes are actually regulated in a mutually exclusive manner in C. elegans. Twenty-two of the 29 clusters consist of two to four homologous ME exons. Ten of the 29 clusters consist of ME exons with the lengths of 3n nt, five of which have too short intervening introns to be excised. Fourteen of the 19 clusters with the ME exons other than 3n nt in length rely at least in part on NMD. Nevertheless, many of the ME exon clusters appear to be strictly regulated. Further molecular and functional analyses of such clusters will elucidate novel mechanisms for mutually exclusive selection of the ME exons in vivo.

Materials and Methods

Total RNAs were extracted from synchronized L1 larvae of N2 and KH1668: smg-2 (yb979) I strains as described previously. RT-PCR was performed essentially as described previously. RT-PCR products were analyzed by using BioAnalyzer (Agilent) as described previously. The sequences of the RT-PCR products were confirmed by direct sequencing or by cloning and sequencing. Sequences of the primers used in the RT-PCR assays are available upon request to Kuroyanagi H.

A list of the GO terms was retrieved from the Gene Ontology website (http://www.geneontology.org/). Fisher’s exact test was performed by using Ekuseru-Toukei 2010 (Social Survey Research Information). Sequence logos were generated by using WebLogo3 at http://weblogo.threeplusone.com/create.cgi.