| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
Jérôme Dejardin,2,
Jamal Tazi,1* and
Johann Soret1*
CNRS, UMR 5535, Institut de Génétique Moléculaire de Montpellier, Montpellier F-34293, France,1 CNRS, UPR 1142, Institut de Génétique Humaine, CNRS, Montpellier F-34293, France2
Received 5 October 2006/ Returned for modification 9 November 2006/ Accepted 18 January 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
SR proteins share a very similar domain organization consisting of one or two RNA recognition motifs at the N terminus and a region rich in arginine-serine dipeptides (RS domain) at the C terminus (18). However, they recognize distinct exonic splicing enhancer elements and are not functionally redundant in vitro or in vivo (5, 30). Evidence that SR proteins exert essential roles in vivo was first reported for B52/SRp55 in Drosophila melanogaster. Indeed, B52-null mutants do not develop past the second-instar larval stage even though the splicing of several endogenous transcripts tested was not significantly altered (37). Later, gene targeting of SF2/ASF in the DT40 chicken B-cell line and of SRp20 in mice was reported to cause lethality (21, 50). More recently, conditional deletion of SC35 in the thymus and in the heart was shown to cause a defect in T-cell maturation and to induce dilated cardiomyopathy in mice (13, 49). Interestingly, a similar deletion of SF2/ASF in the heart resulted in a significantly different phenotype (51), demonstrating that these two related splicing factors play fundamentally distinct roles in cardiac tissue.
While members of the SR protein family have been well characterized at the biochemical level, relatively little is known concerning their physiologically relevant target pre-mRNAs. Identifying the repertoire of genes specifically regulated by the different splicing factors is crucial not only to decipher their role in biological pathways but also to understand how they could be involved in human diseases which are increasingly associated with alterations of splicing processes (15).
As a way to address this question, we have taken advantage of Drosophila as a model system to investigate the consequences of tissue-specific expression of different members of the SR family. When expressed at similar elevated levels in the developing eye, SR proteins induce distinct developmental defects. Coimmunoprecipitation experiments performed with extracts from dASF/SF2 or B52 transgenic flies were used to isolate target transcripts, and these were detected by microarrays analyses. Each SR protein was found to interact with distinct sets of targets, including numerous genes whose functions are related to the phenotypes observed. In the case of B52 transgenic flies, which exhibit the more severe phenotype, we show that the high lethality is correlated with developmental defects in the optical lobe and with splicing alterations of B52-interacting mRNAs that are known to be involved in brain development.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Antibodies and immunofluorescence staining of larval tissues. Eye imaginal discs and brains from third-instar larvae were fixed for 20 min in PBT (1x phosphate-buffered saline, 0.1% Tween 20) with 4% formaldehyde. The tissues were extensively washed in PBT and incubated in PBT-1% bovine serum albumin (BSA) for 1 h at 4°C. Following incubation with the primary antibody diluted in PBT-1% BSA overnight at 4°C, the tissues were washed three times with PBT and incubated with the secondary antibody diluted in PBT-1% BSA for 2 h at 4°C. After three washes in PBT, the tissues were conserved in mounting agent (Vectashield mounting medium H1000; Vector) at 4°C. The polyclonal anti-cleaved-caspase-3 antibodies (no. 9661; Cell Signaling Technology) (kindly provided by P. Lassus), monoclonal anti-Cut 2B10 antibodies, and monoclonal anti-Fas2 1D4 antibodies from the Developmental Studies Hybridoma Bank (provided by F. Girard) were used at a 1/250 dilution. Cy3-conjugated anti-mouse (no. 315-166-047) and anti-rabbit (no. 111-166-047) secondary antibodies were diluted according to the manufacturer's instructions (Jackson Immunoresearch).
Microscopy and images analysis. Eye photographs were obtained using a Nikon Coolpix 990 charge-coupled device camera mounted on an MZFLIII binocular (Leica). All eye photographs represent the average eye pigmentation and phenotypes from each fly population, for both sexes, at 25°C.
Larval eye imaginal disc and brain images were taken with an LSM510 META confocal microscope (Carl Zeiss) using either a 25x plan Neofluar multi-immersion objective (numerical aperture, 0.8) or a 40x apochromat water immersion objective (numerical aperture, 1.2). Samples were sequentially illuminated at 488 and 543 nm, and subsequent GFP and Cy3 fluorescence was detected using 505-530 BP and 560LP filters, respectively.
Larval extract preparation and immunoprecipitation. Anterior quarters of 600 third-instar larvae were dissected in cold 1x phosphate-buffered saline and pooled in 1 ml of cold NETN buffer (50 mM Tris [pH 7.5], 150 mM NaCl, 1 mM EDTA, 0.5 mM MgCl2, 0.25% NP-40, 0.5% Triton X-100) with 400 µg tRNA from Escherichia coli (Roche), 20 µl RNasin (Promega), and 100 µl protease inhibitor cocktail (complete, EDTA-free; Roche). Tissue extracts were homogenized with a plastic pestle in a 1.5-ml microcentrifuge tube and sonicated twice for 5 seconds. The homogenates were centrifuged at 3,500 rpm at 4°C for 5 min. Protein G-Sepharose beads (Amersham Biosciences) were presaturated by incubating with NETN buffer with 1% BSA for 1 h at 4°C on a rotator and sequentially washed with coating buffer (50 mM Tris [pH 7.5], 150 mM NaCl). One hundred microliters of protein G-Sepharose beads was then coated with 15 µl of monoclonal anti-GFP antibody (Roche) in 500 µl of coating buffer with 1% BSA for 1 h at room temperature on a rotator and washed successively three times with coating buffer and NETN buffer. Coated beads were incubated for up to 2 h at 4°C with the supernatants supplemented with 5 µl of RNasin, 6 µl of 100 mM dithiothreitol, 18 µl of 0.5 M EDTA, and 60 µl of protease inhibitor cocktail (Roche). Following centrifugation for 5 min at 3,500 rpm, the supernatant was removed. The beads were washed six times in cold NETN buffer and resuspended in 100 µl of RNase-free water. The beads were then treated with proteinase K, and immunoprecipitated RNAs were isolated following phenol-chloroform extraction and precipitation. For Western blot analyses, total extracts from anterior quarters of 50 third-instar larvaes were prepared as described above.
Western blots. Both total extracts and protein samples collected during the immunoprecipitation protocol were run on an 11% sodium dodecyl sulfate (SDS)-polyacrylamide gel and transferred onto nitrocellulose by electroblotting for 90 min in 10 mM 3-(cyclohexylamino)propanesulfonic acid (pH 11.0) transfer buffer containing 10% ethanol. Blots were probed with monoclonal antibody 104, polyclonal anti-dASF/SF2 (Eurogentec), or monoclonal anti-GFP antibodies (Roche Diagnostics) and revealed with the LumiLight Western blot substrate (Roche) chemiluminescence kit.
Total RNA isolation, poly(A)+ RNA purification, and reverse transcription-PCR (RT-PCR). Anterior quarters of 600 third-instar larvae were homogenized with a plastic pestle in a 1.5-ml microcentrifuge tube, and total RNA extraction was performed with 2 ml of TriReagent (Sigma). First-strand cDNAs were synthesized from 5 µg of total RNA with a first-strand cDNA synthesis kit (Amersham). For PCR analyses, 1/15 of the reaction mixture was amplified with Taq polymerase (Invitrogen), and the cycle number was kept to a minimum to maintain linearity. Primer sequences are available upon request. PCR products were separated on 1.5 to 3% agarose gels containing ethidium bromide and visualized under UV light. Densitometric analyses of the amplified products were performed with the Gnome software.
For transcriptome microarray analyses, poly(A)+ RNAs were purified from 600 µg of total RNA, using the Oligotex mRNA midi kit (QIAGEN) according to the manufacturer's instructions.
Labeled cDNA and aRNA synthesis. One-tenth of immunoprecipitated RNA samples was amplified for 8 h, and corresponding Cy3- and Cy5-labeled antisense RNAs (aRNAs) were synthesized using the Amino Allyl Message Amp II aRNA amplification kit (Ambion) according to the manufacturer's instructions.
First-strand cDNA was synthesized with PowerScript reverse transcriptase (Clontech) from 1.5 µg of poly(A)+ RNA mixed with 1 µl mRNA spike control and 1.6 µl of anchored oligo(dT)20 (2.5 µg/µl; Invitrogen) in the presence of 40 U of RNase OUT (Invitrogen); 10 mM dithiothreitol; dATP, dCTP, and dGTP (0.5 mM each); 0.2 mM dTTP; and 0.3 mM aminoallyl-dUTP (Ambion) in Powerscript incubation buffer. The mixture was incubated at 42°C for 3 h, and the RNA template was degraded by alkaline hydrolysis (0.32 M NaOH at 65°C for 15 min). After neutralization with HEPES (free acid, 0.67 M), the aminoallyl-modified cDNAs were purified on QIAquick PCR purification columns (QIAGEN) according to the manufacturer's instructions, except that QIAquick wash buffer was replaced with 80% ethanol and cDNA was eluted in water. A second purification was done by ethanol precipitation.
The aminoallyl-modified cDNAs were chemically coupled to Cy3 and Cy5 N-hydroxysuccinimidylesters in 50 mM Na2CO3 solution (pH 9.0) for 1 h at room temperature in the dark. The coupling reaction was terminated by the addition of hydroxylamine (1.33 M). The fluorescent cDNAs were subsequently purified from unreacted CyDye using QIAquick PCR purification columns (QIAGEN).
Microarray hybridization. The microarrays designed with the Drosophila Operon 1.1 collection (Operon Biotechnologies), consisting of 14,593 oligonucleotides (70-mer) representing the 13,664 genes and 17,899 transcripts annotated in the Gadfly 3.1 database (Berkeley Drosophila Genome Project [http://www.bdgp.org]), are fully described in the Gene Expression Omnibus repository (http://www.ncbi.nlm.nih.gov/projects/geo/) under the reference entry GPL4275.
Prior to hybridization, excess oligonucleotides were removed from the arrays by shaking the arrays twice for 1 min in 0.2% SDS. Arrays were then washed two times in distilled water. Labeled cDNA was added to microarray hybridization buffer version 2 (Amersham Biosciences) with a final concentration of 50% formamide, denatured at 95°C for 3 min, and applied to the microarrays in individual chambers of an automated slide processor (Amersham Biosciences). Hybridization was carried out at 37°C for 16 h. Hybridized slides were washed at 37°C successively with 1x SSC (0.15 M NaCl plus 0.015 M sodium citrate)-0.2% SDS for 10 min, twice with 0.1x SSC-0.2% SDS for 10 min, with 0.1x SSC, for 1 min and with isopropanol before air drying.
Data acquisition. Microarrays were scanned at 10-µm resolution in both Cy3 and Cy5 channels with a GenePix 4200AL scanner with a variable photo multiplier tube voltage to obtain maximal signal intensities with <0.1% probe saturation. ArrayVision software was used for feature extraction. Spots with high local background or contaminating fluorescence were flagged manually. The local background, calculated for each spot as the median of the fluorescence intensities of four squares surrounding the spot, was subtracted from the foreground fluorescence intensity.
Statistical analysis of GO annotations. Gene ontology (GO) statistics were compared using the program GOstat2 (http://gostat.wehi.edu.au/) (3) in the pools of dASF/SF2 and B52 potential target genes and in the list of 13,664 genes present on the Drosophila microarray. Only overrepresented GO annotations with a maximal P value of 0.5 corrected by the method of Benjamini and Hochberg (3) were considered in the output lists, for subsets of GO hierarchies including development, RNA, apoptosis, and proteolysis.
Microarray accession number. Sample data have been deposited in the Gene Expression Omnibus repository (http://www.ncbi.nlm.nih.gov/projects/geo/) under reference entry GSE5872.
| RESULTS |
|---|
|
|
|---|
To drive expression in differentiating photoreceptor cells (R cells), transgenic flies carrying a single homozygous UAS-GFP-SR element were mated to flies from the GMR-Gal4 line (16). The latter express the yeast transcription factor in the posterior region of the eye imaginal disc under the control of GMR activating sequences derived from the glass gene promoter. As a control, we used transgenic flies carrying a GFP-NLS construct (expressing the GFP fused to an NLS). While no obvious phenotypes were noticed for the homozygous transgenic lines in the absence of Gal4 induction, expression of each of the GFP-SR proteins driven by GMR-Gal4 resulted in progeny flies that displayed defects in eye development of different degrees of severity (Fig. 1A). Mild phenotypes included pigmentation alterations and defects in ommatidia structures (GMR/GFP-d9G8#2 and GMR/GFP-dSC35#5), while stronger defects were characterized by reductions in eye size (GMR/GFP-dASF#1) and significant necrosis (GMR/GFP-dASF#1 and GMR/GFP-B52#5). Variations in the sex ratios and viability of the flies were also observed, again at various degrees, between the different Gal4-expressing progenies. As already reported (23, 24), viability was strongly impaired in B52 transgenic flies (Fig. 1B).
|
To assess whether the severity of the phenotypes depended upon the ratio between the exogenous and corresponding endogenous SR proteins, Western blot analyses were performed with dASF/SF2-specific antibodies and monoclonal antibody 104 on protein extracts obtained from anterior parts of third-instar larvae. As shown in Fig. 1E and F, the ratios of exogenous to endogenous SR protein expression were very similar in the GMR/GFP-dASF#1 and GMR/GFP-B52#6 larvae. Assuming that the GMR-expressing tissues represent up to 10% of the third-instar larva anterior quarters, these experiments indicated that the expression level of the GFP-SR proteins in these tissues is 5 to 10 times higher than that in their endogenous counterparts.
Taken together, these observations demonstrate that, when expressed at similar levels, distinct SR splicing factors induce different developmental defects.
B52, but not dASF/SF2, GMR-driven expression severely impairs brain development in Drosophila. To obtain a more detailed insight into the developmental defects giving rise to the observed phenotypes, direct and immunofluorescence analyses were performed on eye imaginal discs from third-instar larvae of the control GMR/GFP-NLS, GMR/GFP-dASF#1, and GMR/GFP-B52#6 transgenics. In all cases, GFP fluorescence was restricted to the differentiating R cells localized in the posterior part of the eye imaginal disc (Fig. 2) after the passage of the morphogenetic furrow (MF) (16). In GMR/GFP-dASF#1 larvae, the organization of the differentiating R cells was slightly modified compared to that in control GMR/GFP-NLS counterparts (Fig. 2A and B). As shown in Fig. 2C to F, this organization was dramatically altered in GMR/GFP-B52#6 larvae.
|
Since these phenotypes could result from increased apoptosis of the corresponding precursor cells, immunofluorescence staining was performed with an anti-active-caspase-3 antibody, which detects activated Drosophila effector caspases drICE and DCP-1 (52). As expected, apoptotic cells were detected on the posterior side of the MF in the control larvae (Fig. 2J), whereas significant apoptosis was observed throughout the posterior part of the eye imaginal disc in GMR/GFP-dASF#1 larvae (Fig. 2K). In GMR/GFP-B52#6 progeny, a high number of apoptotic cells were concentrated in the more posterior part of the eye disc (Fig. 2L), where R-cell organization and cone cell differentiation were primarily altered.
Taken together, these observations strongly suggest that the GMR-driven dASF/SF2 expression interferes with the differentiation program and/or triggers apoptosis of cone cells, whereas expression of GFP-B52 affects both R and cone cells. While this difference could account for eye phenotypes of distinct severity, it cannot explain the 100% lethality systematically observed at the third-instar larval stage in the GMR/GFP-B52#6 progeny.
Several studies previously revealed that development of the adult insect optic lobe, and more precisely the lamina ganglion, is strictly dependent on innervation from the compound eye (41). To determine whether the alterations of R-cell differentiation observed in GMR/GFP-B52#6 larvae impair the development of the optic lobe, we first analyzed the expression of GFP-NLS and GFP-SR proteins in this tissue. As shown in Fig. 3A and B, GFP-NLS expression was readily observed in axons of the R1 to R6 cells, which terminate their outgrowth in the lamina (32). GFP-NLS was also detected, although at a lower level, in R7 and R8 axons targeted to the distal medulla. These observations indicate that, in spite of its fusion to an NLS domain, the GFP protein is not restricted to the R-cell nucleus. Similar observations were obtained for the GFP-dASF protein (Fig. 3E and F) revealing on the one hand that GFP-SR protein localization is not restricted to the nuclei of R cells, a result that is similar to one recently published (17), and on the other hand that GFP-dASF/SF2 expression does not significantly affect the localization of the R-cell axons in the lamina. In sharp contrast, GFP-B52 induced a strong alteration in the projection of R-cell axons.
|
dASF/SF2 and B52 are associated with distinct sets of mRNA targets. To identify dASF- and B52-specific mRNA targets accounting for the distinct phenotypes induced in the transgenic flies, we analyzed transcripts coimmunoprecipitated with the GFP-SR proteins, assuming that their specificity is not significantly different from that of their endogenous counterparts. These experiments were performed using a monoclonal anti-GFP antibody in extracts from anterior quarters of third-instar larvae (Fig. 4A), and copurified mRNA species were used as probes to hybridize long oligonucleotide (QIAGEN Operon) microarrays specific for the 13,664 genes of the GadFly 3.2 Drosophila genome database (http://flybase.bio.indiana.edu). In this assay, cRNAs prepared from GFP-NLS-coimmunoprecipitated mRNAs were used as a nonspecific interaction control.
|
The results obtained from two independent immunoprecipitation experiments (see Table S2 in the supplemental material) indicate that mRNA species significantly enriched (enrichment factor of >1.5) in coimmunoprecipitates from GMR/GFP-dASF#1 and GMR/GFP-B52#6 larval extracts correspond to an equivalent but small number of genes (76 and 77, respectively). Among them, only 10 were consistently common to both GFP-dASF and GFP-B52 immunoprecipitates, confirming the lack of general redundancy between the dASF/SF2 and B52 activities. The majority of the genes identified by transcripts associated with GFP-dASF and GFP-B52 were found to belong to the same three functional classes (Fig. 4B, upper part) involved in RNA metabolism (dASF, 23.9%; B52, 30.9%), development (dASF, 32.2%; B52, 22.2%), and apoptosis or proteolysis (dASF, 7.8%; B52, 13.6%). Moreover, compared to the 13,664 genes represented in the array, both the dASF and B52 potential target genes revealed a significant enrichment in gene ontology terms associated with developmental processes, and more precisely in neural system development, RNA or mRNA metabolism, splicing, apoptosis, and proteolysis (see Table 3 in the supplemental material).
Around 60% of the genes corresponding to enriched target transcripts were not previously reported to be alternatively spliced (Fig. 4B, lower part), an observation possibly reflecting the involvement of SR proteins in constitutive splicing as well as in different steps of mRNA life such as export, stability, and translation (40). When only alternatively spliced genes were considered, 80 to 90% of them were found to be distributed in the functional classes concerning RNA metabolism, development, and apoptosis (Fig. 4C, upper part). Interestingly, the alternatively spliced genes associated with dASF and B52 are particularly well represented in the "development" class, where they correspond to 70 to 80% of the genes identified in this category (Fig. 4C, lower part). Taken together, these observations suggest that SR proteins may participate in developmental processes mainly through their activity in alternative splicing regulation.
Since the developmental defects induced by GFP-dASF and GFP-B52 in transgenic larvae differ essentially at the level of the optic lobe neurogenesis, we then searched for genes specifically involved in eye and nervous system development among the GFP-dASF and GFP-B52 targets (Table 1). As shown in Fig. 4D, approximately 70% of these genes are related to brain development in GFP-B52 targets, whereas the predominant fraction of these GFP-dASF targets are involved in eye development. Again, 70 to 80% of these targets correspond to alternatively spliced genes, raising the possibility that alterations of their splicing pattern are responsible for the observed phenotypes.
|
|
Similarly, the splicing profile of the potential B52 target Pyd (Polychaetoid), which is involved in eye development, as well as those of Nup154, Tm1 (Tropomyosin 1), and mnb (minibrain), which play a role in nervous system development, was essentially affected in GMR/GFP-B52#6 larvae (Fig. 6). To confirm that these genes indeed correspond to genuine B52 targets, their splicing patterns were also analyzed in transgenic larvae expressing an RNA aptamer which corresponds to tandemly repeated high-affinity BBS, which was previously reported to efficiently titrate this SR protein and inhibit B52-dependent splicing in vitro (42). The expression of this specific B52 inhibitor, under the control of the GMR driver, also induces eye developmental defects in transgenic flies (see Fig. S1A in the supplemental material) and is able to rescue the eye phenotype induced by B52 overexpression (see Fig. S1B and C in the supplemental material), thereby confirming its efficacy to counteract the activity of this SR protein in vivo. As shown in Fig. 6, the splicing profiles of the four genes were significantly altered in UAS-BBS transgenics, in a manner either similar to that observed in GMR/GFP-B52#6 larvae (Nup 154 and TM1) or clearly different (Pyd and mnb). Thus, both increased B52 expression and its in vivo depletion, targeted in the same tissue, resulted in modifications of alternative splicing patterns for multiple genes during Drosophila visual system development.
|
Taken together, these observations indicate that genes identified by coimmunoprecipitation experiments coupled to microarray analyses likely correspond to genuine dASF and B52 targets and confirm that these SR proteins regulate alternative splicing of specific subsets of pre-mRNAs that can be linked to altered phenotypes.
| DISCUSSION |
|---|
|
|
|---|
Our immunofluorescence experiments strongly suggest that the GMR-driven dASF/SF2 expression interferes with the differentiation program and/or triggers apoptosis of cone cells, whereas expression of GFP-B52 affects both R and cone cells. Interestingly, B52 overexpression in SL2 cells was recently shown to strongly repress the expression of dE2F1, an activator of cell proliferation, while the loss of expression of this SR protein is associated with splicing default of the dE2F2 gene, an antagonist of dE2F1 function (34). Consistently, the eye phenotypes of dE2F1 mutants resulting from cell cycle defects are suppressed by B52 mutations that probably affect the expression of dE2F2 mRNA. These results are in good agreement with our observation that the dE2F2 mRNA level is clearly increased in GMR/GFP-B52#6 larvae and with immunofluorescence experiments performed with antibodies raised against phosphorylated histone H3, which indicated that the number of dividing cells in the second mitotic wave prior to MF passage is significantly reduced in GMR/GFP-B52#6 eye imaginal disc compared to the corresponding control tissue (data not shown).
Besides the disorganization of ommatidia in the developing eye, a striking difference between the phenotypes induced by dASF/SF2 and B52 was the very high lethality provoked by the latter SR protein in third-instar larvae. Direct fluorescence analyses as well as immunostaining experiments have revealed that in B52 transgenic larvae, the R-cell axons project aberrantly within the optic lobe, an alteration likely responsible for the absence of normal lamina neurogenesis. Indeed, numerous studies have reported that retinal axon projections induce the organization of laminar and medullar neurons (41) and coordinate the distribution of glia to multiple target destinations, where they are required for axon guidance and neuronal survival (12). Strikingly, visual system mutants showing a reduced or absent retina also have a completely missing lamina and a greatly reduced medulla, indicating that retinal innervation is required for further development of these ganglia (33). Interestingly, the phenotype induced by B52 is very similar to that of eyeless mutants, which lack a lamina and display a significant size reduction of the other visual ganglia (10). In some mutants, eye and optic lobe defects are associated with an important lethality, suggesting that abnormal neurogenesis in the optic lobe impairs vital functions.
In the past few years, global assays aimed at identifying cellular RNA targets of splicing regulators have involved coimmunoprecipitation strategies followed by microarray analysis (25, 47), genomic systematic evolution of ligands by exponential enrichment (SELEX) (22), and, more recently, development of alternative splicing microarrays (5, 36, 48). In the present study, immunoprecipitation of RNA species bound to GFP-dASF and GFP-B52 proteins expressed in the Drosophila eye identified an equivalent number of candidate target genes (76 and 77, respectively). Among them, only 40% were already known to be alternatively spliced. Since many genes subjected to alternative splicing are not yet annotated (11, 43), a higher proportion of the candidate targets could actually correspond to alternatively spliced genes, as reflected by our RT-PCR experiments which revealed two neurotactin mRNA isoforms that had not been previously described (see Fig. S3 in the supplemental material). Moreover, 15 to 20% of the dASF/SF2 and B52 potential targets correspond to still-uncharacterized genes whose splicing status remains to be established.
Identification of constitutively spliced genes as candidate targets can also reflect the involvement of SR proteins in constitutive splicing. Indeed, it was recently reported that dASF/SF2 and B52 control 105 and 22 constitutive splicing events, respectively, in SL2 cells (5). On the other hand, such candidates could be subjected to other activities of SR proteins, such as regulation of mRNA stability and translation (27, 39). It is noteworthy that one of the B52 targets found in our screen (CG9080) has already been identified by genomic SELEX and shown to contain a BBS in the 3' untranslated region of the corresponding mRNA (22).
We also observed that only 6% of the identified genes correspond to targets of both dASF/SF2 and B52. Since splicing microarray analyses revealed that 5% of the splicing events are similarly affected by knockdown of either SR protein (5), these observations confirm the unique character of each SR protein and suggest that the functional overlap between these splicing factors is similar in the different tissues of a whole organism.
Results reported in this study indicate that our strategy efficiently identified specific target genes whose expression is modulated by dASF/SF2 and B52. To date, two different approaches have been used to identify B52 targets or alternative splicing events controlled by dASF/SF2 and B52. Genomic SELEX has identified 15 B52 candidate target genes, among which four exhibit splicing defects in a B52-null mutant (22). Splicing microarrays have revealed that dASF/SF2 and B52 control 319 and 107 splicing events, respectively, in SL2 cells (5), with 6 out of 6 tested being validated by RT-PCR experiments. In the present study, 36 dASF/SF2 and/or B52 potential target genes were analyzed by the same approach (Fig. 5 and 6; see Fig. S3 and S4 in the supplemental material), and 18 (50%) are significantly affected at the level of their splicing pattern in the corresponding transgenic larvae. Moreover, by using transgenic flies expressing an RNA aptamer capable of titrating out the endogenous B52 protein (42; this study), we have confirmed that B52 candidates validated by RT-PCR experiments (Fig. 6; see Fig. S4 in the supplemental material) represent genuine B52 targets. Of note, our RT-PCR experiments were designed to analyze sequence-defined splicing events reported in databases, and it is possible that candidates considered negative are nevertheless affected in other regions or are regulated by SR proteins at steps different from alternative splicing.
Comparison of the dASF/SF2 and B52 potential targets identified by the different approaches used so far indicates that three alternative splicing events regulated by dASF/SF2 in SL2 cells (5) affect the genes lola, flo, and garz, which we also identified as targets for this SR protein. On the other hand, three out of the four B52 targets (Rx, Mio, and RhoGap16F) validated by Kim et al. (22) were not identified in our screen. This is not surprising, since genomic SELEX identified all genomic sequences containing a high-affinity BBS, whereas our coimmunoprecipitation strategy was aimed at characterizing mRNAs bound by B52 and expressed in a particular cell type or tissue. Indeed, Rx is significantly expressed in the Drosophila larval brain, but not in R cells (14), and to our knowledge, expression of Mio and RhoGap16F has not been reported to occur in R cells. Among the four B52 targets found by both the genomic SELEX and splicing microarray analyses (fur1, lola, RhoGAp16F, and Syndecan), only lola (longitudinals lacking) has been identified in the present study. While these observations possibly reflect, as observed for fur2 (38), the lack of expression of fur1 in R cells, the situation is more surprising in the case of Syndecan, since expression of this gene has been detected in R-cell axons (35). This target, however, was not validated in B52 mutant larvae (22), suggesting that the event detected on splicing microarrays is restricted to specific tissues. On the other hand, such observations were not unexpected in the light of the many fundamental differences between the RNAi approach examining splicing changes upon protein depletion and our coimmunoprecipitation strategy aimed at identifying mRNAs bound to overexpressed proteins. Along this line, one can speculate that RNAi might primarily affect mRNAs with low-affinity sites since protein reduction is not complete, while overexpression might push occupancy of lower-affinity sites on a different set of mRNAs. In agreement with this hypothesis, a search for high-affinity binding sequences specifically recognized by the human ASF/SF2 homolog (44) and the B52/SRp55 protein (see reference 42 and references therein) did not reveal a specific enrichment of these motifs around the splice sites of the regulated exons in the validated dASF/SF2 and B52 targets (data not shown). However, the restricted number of splicing events considered might explain the lack of statistical significance and the only marginal enrichment of highly scoring putative binding sites for both SR proteins that we observed in our computational analysis.
Our analyses indicate that genes such as neurotactin, chromosome bows, and longitudinals lacking that are involved in axon guidance and neurogenesis constitute potential dASF/SF2 targets. Splicing alterations of the corresponding pre-mRNAs, however, do not appear to affect the correct projection of R-cell axons in the transgenic flies. On the other hand, one can envisage that splicing alterations of lola, which also corresponds to a B52 candidate target, are much less deleterious in dASF/SF2 than in B52 transgenics.
Another interesting B52 target involved in neurogenesis and whose alternative splicing is clearly altered in GMR/GFP-B52 transgenic larvae is minibrain (mnb). This gene encodes a serine/threonine protein kinase which appears to play an essential role during postembryonic neurogenesis (45). Mutations in the mnb gene induce an abnormal spacing of neuroblasts in the larval brain due to a reduced production of neuronal progeny. The human homolog (MNBH/DYRK1A) of the mnb gene is overexpressed in brains from fetuses with Down's syndrome (19), and the phenotypes of DYRK1A transgenic and knockout mice support the idea that DYRK1A contributes to some of the neuropathological traits described in Down's syndrome (2).
Although further experiments will be needed to precisely identify the genes whose splicing alterations account for the specific phenotypes induced by dASF/SF2 and B52, the present study indicates that the coimmunoprecipitation approach linked to microarray analyses represents a straightforward strategy to discover transcripts and associated pathways regulated by individual SR proteins and to understand how they participate in normal and pathological developmental processes.
| ACKNOWLEDGMENTS |
|---|
M. Gabut is a recipient of a fellowship from the Association pour la Recherche contre le Cancer (ARC). This work was supported by grants from the Agence Nationale de la Recherche (ANR-05-BLAN-0261-01) and European Alternative Splicing Network of Excellence (EURASNET, FP6 Life Sciences, Genomics and Biotechnology for Health).
| FOOTNOTES |
|---|
Published ahead of print on 5 February 2007. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
Present address: Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Canada M5S 3E1. ![]()
Present address: Harvard Medical School, Department of Genetics, and Massachusetts General Hospital, Department of Molecular Biology, Boston, MA 02114. ![]()
| REFERENCES |
|---|
|
|
|---|
2. Altafaj, X., M. Dierssen, C. Baamonde, E. Marti, J. Visa, J. Guimera, M. Oset, J. R. Gonzalez, J. Florez, C. Fillat, and X. Estivill. 2001. Neurodevelopmental delay, motor abnormalities and cognitive deficits in transgenic mice overexpressing Dyrk1A (minibrain), a murine model of Down's syndrome. Hum. Mol. Genet. 10:1915-1923.
3. Beissbarth, T., and T. P. Speed. 2004. GOstat: find statistically overrepresented gene ontologies within a group of genes. Bioinformatics 20:1464-1465.
4. Black, D. L. 2003. Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem. 72:291-336.[CrossRef][Medline]
5. Blanchette, M., R. E. Green, S. E. Brenner, and D. C. Rio. 2005. Global analysis of positive and negative pre-mRNA splicing regulators in Drosophila. Genes Dev. 19:1306-1314.
6. Blochlinger, K., L. Y. Jan, and Y. N. Jan. 1993. Postembryonic patterns of expression of cut, a locus regulating sensory organ identity in Drosophila. Development 117:441-450.[Abstract]
7. Brachmann, C. B., and R. L. Cagan. 2003. Patterning the fly eye: the role of apoptosis. Trends Genet. 19:91-96.[CrossRef][Medline]
8. Brand, A. H., A. S. Manoukian, and N. Perrimon. 1994. Ectopic expression in Drosophila. Methods Cell Biol. 44:635-654.[Medline]
9. Cagan, R. L., and D. F. Ready. 1989. The emergence of order in the Drosophila pupal retina. Dev. Biol. 136:346-362.[CrossRef][Medline]
10. Callaerts, P., S. Leng, J. Clements, C. Benassayag, D. Cribbs, Y. Y. Kang, U. Walldorf, K. F. Fischbach, and R. Strauss. 2001. Drosophila Pax-6/eyeless is essential for normal adult brain structure and function. J. Neurobiol. 46:73-88.[CrossRef][Medline]
11. Celniker, S. E., and G. M. Rubin. 2003. The Drosophila melanogaster genome. Annu. Rev. Genomics Hum. Genet. 4:89-117.[CrossRef][Medline]
12. Dearborn, R., Jr., and S. Kunes. 2004. An axon scaffold induced by retinal axons directs glia to destinations in the Drosophila optic lobe. Development 131:2291-2303.
13. Ding, J. H., X. Xu, D. Yang, P. H. Chu, N. D. Dalton, Z. Ye, J. M. Yeakley, H. Cheng, R. P. Xiao, J. Ross, J. Chen, and X. D. Fu. 2004. Dilated cardiomyopathy caused by tissue-specific ablation of SC35 in the heart. EMBO J. 23:885-896.[CrossRef][Medline]
14. Eggert, T., B. Hauck, N. Hildebrandt, W. J. Gehring, and U. Walldorf. 1998. Isolation of a Drosophila homolog of the vertebrate homeobox gene Rx and its possible role in brain and eye development. Proc. Natl. Acad. Sci. USA 95:2343-2348.
15. Faustino, N. A., and T. A. Cooper. 2003. Pre-mRNA splicing and human disease. Genes Dev. 17:419-437.
16. Freeman, M. 1997. Cell determination strategies in the Drosophila eye. Development 124:261-270.[Abstract]
17. Glanzer, J., K. Y. Miyashiro, J. Y. Sul, L. Barrett, B. Belt, P. Haydon, and J. Eberwine. 2005. RNA splicing capability of live neuronal dendrites. Proc. Natl. Acad. Sci. USA 102:16859-16864.
18. Graveley, B. R. 2000. Sorting out the complexity of SR protein functions. RNA. 6:1197-1211.[CrossRef][Medline]
19. Guimera, J., C. Casas, C. Pucharcos, A. Solans, A. Domenech, A. M. Planas, J. Ashley, M. Lovett, X. Estivill, and M. A. Pritchard. 1996. A human homologue of Drosophila minibrain (MNB) is expressed in the neuronal regions affected in Down syndrome and maps to the critical region. Hum. Mol. Genet. 5:1305-1310.
20. Huang, Y., and J. A. Steitz. 2001. Splicing factors SRp20 and 9G8 promote the nucleocytoplasmic export of mRNA. Mol. Cell 7:899-905.[CrossRef][Medline]
21. Jumaa, H., G. Wei, and P. J. Nielsen. 1999. Blastocyst formation is blocked in mouse embryos lacking the splicing factor SRp20. Curr. Biol. 9:899-902.[CrossRef][Medline]
22. Kim, S., H. Shi, D. K. Lee, and J. T. Lis. 2003. Specific SR protein-dependent splicing substrates identified through genomic SELEX. Nucleic Acids Res. 31:1955-1961.
23. Kraus, M. E., and J. T. Lis. 1994. The concentration of B52, an essential splicing factor and regulator of splice site choice in vitro, is critical for Drosophila development. Mol. Cell. Biol. 14:5360-5370.
24. Labourier, E., E. Allemand, S. Brand, M. Fostier, J. Tazi, and H. M. Bourbon. 1999. Recognition of exonic splicing enhancer sequences by the Drosophila splicing repressor RSF1. Nucleic Acids Res. 27:2377-2386.
25. Labourier, E., M. Blanchette, J. W. Feiger, M. D. Adams, and D. C. Rio. 2002. The KH-type RNA-binding protein PSI is required for Drosophila viability, male fertility, and cellular mRNA processing. Genes Dev. 16:72-84.
26. Lai, M. C., and W. Y. Tarn. 2004. Hypophosphorylated ASF/SF2 binds TAP and is present in messenger ribonucleoproteins. J. Biol. Chem. 279:31745-31749.
27. Lemaire, R., J. Prasad, T. Kashima, J. Gustafson, J. L. Manley, and R. Lafyatis. 2002. Stability of a PKCI-1-related mRNA is controlled by the splicing factor ASF/SF2: a novel function for SR proteins. Genes Dev. 16:594-607.
28. Li, X., and J. L. Manley. 2005. Inactivation of the SR protein splicing factor ASF/SF2 results in genomic instability. Cell 122:365-378.[CrossRef][Medline]
29. Lin, D. M., R. D. Fetter, C. Kopczynski, G. Grenningloh, and C. S. Goodman. 1994. Genetic analysis of Fasciclin II in Drosophila: defasciculation, refasciculation, and altered fasciculation. Neuron 13:1055-1069.[CrossRef][Medline]
30. Liu, H. X., M. Zhang, and A. R. Krainer. 1998. Identification of functional exonic splicing enhancer motifs recognized by individual SR proteins. Genes Dev. 12:1998-2012.
31. Maniatis, T., and B. Tasic. 2002. Alternative pre-mRNA splicing and proteome expansion in metazoans. Nature 418:236-243.[CrossRef][Medline]
32. Poeck, B., S. Fischer, D. Gunning, S. L. Zipursky, and I. Salecker. 2001. Glial cells mediate target layer selection of retinal axons in the developing visual system of Drosophila. Neuron 29:99-113.[CrossRef][Medline]
33. Rachidi, M., C. Lopes, and J. C. Benichou. 1997. Genetical analysis of visual system disorganizer (vid), a new gene involved in normal development of eye and optic lobe of the brain in Drosophila melanogaster. Genetica 99:31-45.[Medline]
34. Rasheva, V. I., D. Knight, P. Bozko, K. Marsh, and M. V. Frolov. 2006. Specific role of the SR protein splicing factor B52 in cell cycle control in Drosophila. Mol. Cell. Biol. 26:3468-3477.
35. Rawson, J. M., B. Dimitroff, K. G. Johnson, J. M. Rawson, X. Ge, D. Van Vactor, and S. B. Selleck. 2005. The heparan sulfate proteoglycans Dally-like and Syndecan have distinct functions in axon guidance and visual-system assembly in Drosophila. Curr. Biol. 15:833-838.[CrossRef][Medline]
36. Relogio, A., C. Ben Dov, M. Baum, M. Ruggiu, C. Gemund, V. Benes, R. B. Darnell, and J. Valcarcel. 2005. Alternative splicing microarrays reveal functional expression of neuron-specific regulators in Hodgkin lymphoma cells. J. Biol. Chem. 280:4779-4784.
37. Ring, H. Z., and J. T. Lis. 1994. The SR protein B52/SRp55 is essential for Drosophila development. Mol. Cell. Biol. 14:7499-7506.
38. Roebroek, A. J., T. A. Ayoubi, J. W. Creemers, I. G. Pauli, and W. J. Van de Ven. 1995. The Dfur2 gene of Drosophila melanogaster: genetic organization, expression during embryogenesis, and pro-protein processing activity of its translational product Dfurin2. DNA Cell Biol. 14:223-234.[Medline]
39. Sanford, J. R., N. K. Gray, K. Beckmann, and J. F. Caceres. 2004. A novel role for shuttling SR proteins in mRNA translation. Genes Dev. 18:755-768.
40. Sanford, J. R., J. Ellis, and J. F. Caceres. 2005. Multiple roles of arginine/serine-rich splicing factors in RNA processing. Biochem. Soc Trans. 33:443-446.[CrossRef][Medline]
41. Selleck, S. B., and H. Steller. 1991. The influence of retinal innervation on neurogenesis in the first optic ganglion of Drosophila. Neuron 6:83-99.[CrossRef][Medline]
42. Shi, H., B. E. Hoffman, and J. T. Lis. 1999. RNA aptamers as effective protein antagonists in a multicellular organism. Proc. Natl. Acad. Sci. USA 96:10033-10038.
43. Stolc, V., Z. Gauhar, C. Mason, G. Halasz, M. F. van Batenburg, S. A. Rifkin, S. Hua, T. Herreman, W. Tongprasit, P. E. Barbano, H. J. Bussemaker, and K. P. White. 2004. A gene expression map for the euchromatic genome of Drosophila melanogaster. Science 306:655-660.
44. Tacke, R., and J. L. Manley. 1995. The human splicing factors ASF/SF2 and SC35 possess distinct, functionally significant RNA binding specificities. EMBO J. 14:3540-3551.[Medline]
45. Tejedor, F., X. R. Zhu, E. Kaltenbach, A. Ackermann, A. Baumann, I. Canal, M. Heisenberg, K. F. Fischbach, and O. Pongs. 1995. minibrain: a new protein kinase family involved in postembryonic neurogenesis in Drosophila. Neuron 14:287-301.[CrossRef][Medline]
46. Tomlinson, A., and D. F. Ready. 1987. Cell fate in the Drosophila ommatidium. Dev. Biol. 123:264-275.
47. Ule, J., K. B. Jensen, M. Ruggiu, A. Mele, A. Ule, and R. B. Darnell. 2003. CLIP identifies Nova-regulated RNA networks in the brain. Science 302:1212-1215.
48. Ule, J., A. Ule, J. Spencer, A. Williams, J. S. Hu, M. Cline, H. Wang, T. Clark, C. Fraser, M. Ruggiu, B. R. Zeeberg, D. Kane, J. N. Weinstein, J. Blume, and R. B. Darnell. 2005. Nova regulates brain-specific splicing to shape the synapse. Nat. Genet. 37:844-852.[CrossRef][Medline]
49. Wang, H. Y., X. Xu, J. H. Ding, J. R. Bermingham, Jr., and X. D. Fu. 2001. SC35 plays a role in T cell development and alternative splicing of CD45. Mol. Cell 7:331-342.[CrossRef][Medline]
50. Wang, J., Y. Takagaki, and J. L. Manley. 1996. Targeted disruption of an essential vertebrate gene: ASF/SF2 is required for cell viability. Genes Dev. 10:2588-2599.
51. Xu, X., D. Yang, J. H. Ding, W. Wang, P. H. Chu, N. D. Dalton, H. Y. Wang, J. R. Bermingham, Jr., Z. Ye, F. Liu, M. G. Rosenfeld, J. L. Manley, J. Ross, Jr., J. Chen, R. P. Xiao, H. Cheng, and X. D. Fu. 2005. ASF/SF2-regulated CaMKIIdelta alternative splicing temporally reprograms excitation-contraction coupling in cardiac muscle. Cell 120:59-72.[CrossRef][Medline]
52. Yu, S. Y., S. J. Yoo, L. Yang, C. Zapata, A. Srinivasan, B. A. Hay, and N. E. Baker. 2002. A pathway of signals regulating effector and initiator caspases in the developing Drosophila eye. Development 129:3269-3278.
53. Zhang, Z., and A. R. Krainer. 2004. Involvement of SR proteins in mRNA surveillance. Mol. Cell 16:597-607.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||