Unproductive alternative splicing and nonsense mRNAs: A widespread phenomenon among plant circadian clock genes
© Filichkin and Mockler; licensee BioMed Central Ltd. 2012
Received: 14 December 2011
Accepted: 2 July 2012
Published: 2 July 2012
Recent mapping of eukaryotic transcriptomes and spliceomes using massively parallel RNA sequencing (RNA-seq) has revealed that the extent of alternative splicing has been considerably underestimated. Evidence also suggests that many pre-mRNAs undergo unproductive alternative splicing resulting in incorporation of in-frame premature termination codons (PTCs). The destinies and potential functions of the PTC-harboring mRNAs remain poorly understood. Unproductive alternative splicing in circadian clock genes presents a special case study because the daily oscillations of protein expression levels require rapid and steep adjustments in mRNA levels.
We conducted a systematic survey of alternative splicing of plant circadian clock genes using RNA-seq and found that many Arabidopsis thaliana circadian clock-associated genes are alternatively spliced. Results were confirmed using reverse transcription polymerase chain reaction (RT-PCR), quantitative RT-PCR (qRT-PCR), and/or Sanger sequencing. Intron retention events were frequently observed in mRNAs of the CCA1/LHY-like subfamily of MYB transcription factors. In contrast, the REVEILLE2 (RVE2) transcript was alternatively spliced via inclusion of a "poison cassette exon" (PCE). The PCE type events introducing in-frame PTCs are conserved in some mammalian and plant serine/arginine-rich splicing factors. For some circadian genes such as CCA1 the ratio of the productive isoform (i.e., a representative splice variant encoding the full-length protein) to its PTC counterpart shifted sharply under specific environmental stress conditions.
Our results demonstrate that unproductive alternative splicing is a widespread phenomenon among plant circadian clock genes that frequently generates mRNA isoforms harboring in-frame PTCs. Because LHY and CCA1 are core components of the plant central circadian oscillator, the conservation of alternatively spliced variants between CCA1 and LHY and for CCA1 across phyla  indicates a potential role of nonsense transcripts in regulation of circadian rhythms. Most of the alternatively spliced isoforms harbor in-frame PTCs that arise from full or partial intron retention events. However, a PTC in the RVE2 transcript is introduced through a PCE event. The conservation of AS events and modulation of the relative abundance of nonsense isoforms by environmental and diurnal conditions suggests possible regulatory roles for these alternatively spliced transcripts in circadian clock function. The temperature-dependent expression of the PTC transcripts among members of CCA1/LHY subfamily indicates that alternative splicing may be involved in regulation of the clock temperature compensation mechanism.
This article was reviewed by Dr. Eugene Koonin, Dr. Chungoo Park (nominated by Dr. Kateryna Makova), and Dr. Marcelo Yanovsky (nominated by Dr. Valerian Dolja).
KeywordsArabidopsis thaliana Alternative splicing Circadian clock RNA-seq Intron retention Cassette exon Nonsense mRNAs Premature termination codon CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) LATE ELONGATED HYPOCOTYL (LHY) REVEILLE 2 (RVE2).
Interrogation of the transcriptomes of several plant species using genome-scale expression microarrays revealed that many plant genes display diurnally driven rhythmic daily changes in mRNA abundance . The rhythmic expression of the majority of the core clock and circadian associated genes is phased to remarkably similar times of day in Arabidopsis thaliana Populus (poplar), Brachypodium, and Oryza (rice). Furthermore, the core circadian clock genes of Oryza subspecies japonica and indica displayed nearly indistinguishable oscillation profiles and phase calls . The microarray platform has been a method of choice for genome-scale profiling of daily cyclical changes in mRNA abundance. However, expression microarrays have a substantial limitation: in most cases probe sets cannot differentiate among alternatively spliced isoforms of an mRNA. Therefore, for many alternatively spliced genes the overall expression values represent an average of signals generated by the mixture of mRNA variants. In contrast, the massively parallel RNA sequencing (RNA-seq) method is a high-resolution alternative to the expression arrays. RNA-seq coupled with bioinformatic analyses allows direct inference of transcript structures relative to the reference genome and unbiased characterization of alternatively spliced events.
Our recent RNA-seq survey of the Arabidopsis transcriptome suggested that some circadian clock genes undergo extensive alternative splicing (AS) yielding single or multiple splice isoforms . Further, some specific splicing events were conserved across distantly related species of mono- and dicotyledonous plants . For example, RNA-seq data predicted that alternatively spliced mRNA of the CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) gene retains intron 4 (I4R) resulting in the introduction of several in-frame premature termination codons (PTCs). Validation of the I4R event in CCA1 gene using an isoform-specific set of primers and quantitative reverse transcription polymerase chain reaction (qRT-PCR) showed that nonsense I4R CCA1 transcripts can accumulate at high levels relative to the productive isoform in response to certain environmental stimuli such as high-intensity light stress . Furthermore, despite substantial differences in nucleotide sequence, intron length, and the context of splice signals in homologs, the I4R event in CCA1 mRNA was highly conserved among CCA1 clock genes in both eudicots (A. thaliana and P. trichocarpa) and monocots (O. sativa and B. distachyon) . This notable conservation suggests that the some AS events in circadian genes are highly conserved and therefore may play specific roles in the functioning of the plant circadian clock, possibly through the regulation of the abundance of circadian clock mRNAs. Here, we conducted a further survey and characterization of the AS events in plant circadian genes. Collectively, the RNA-seq, RT-PCR, qRT-PCR, and/or Sanger sequencing data suggested that several plant circadian genes of the central oscillator are extensively alternatively spliced.
Strategy for discovery and identification of novel alternative splicing events in circadian genes
Pre-mRNAs of the CCA1/LHY-like subfamily of MYB transcription factors undergo extensive and conserved alternative splicing
Validation of predicted alternative splicing events in Arabidopsis circadian genes
Predicted AS event
Fragment size, bp
Retained intron 4 (AT1G01060.2)
RT-PCR, qRT-PCR, Sanger sequencing
Retained intron 9 (AT1G01060.3)
Retained intron 7
Retained intron 8
Retained intron 4
RT-PCR, qRT-PCR, Sanger sequencing
Retained intron 3
RT-PCR, Sanger sequencing
Retained intron 3
Alternative Poison Cassette Exon (PCE) in intron 1
RT-PCR, Sanger sequencing
Retained intron 3
RT-PCR, Sanger sequencing
Retained intron 3
Retained intron 4
Retained intron 6
Retained intron 1
Retained intron 7
Retained intron 8
Retained portion of the intron 1 (alternative donor splice site)
RT-PCR, qRT-PCR, Sanger sequencing
Retained portion of the intron 1 (alternative donor splice site)
Unproductively spliced circadian transcripts can harbor PTCs at various locations in the mRNA
The location of a PTC within an mRNA can influence recognition of the transcript by the nonsense mediated mRNA decay (NMD) pathway in both mammals and plants and therefore may determine its NMD-eliciting properties. The distribution of AS events observed suggested that alternatively spliced circadian transcripts can harbor both early and late classes of PTCs (Table 1 and data not shown). The positions of retained introns along the length of the transcripts varied from early in RVE7 (I1R) to intermediate in CKB3 (I4R) to late in LCL1 (I7R) (Table 1). We monitored the relative levels of the LCL1 (I7R) transcripts during LCL1 peak expression using qRT-PCR. The isoforms with both retained (I7R) and spliced (I7S) intron 7 displayed similar temporal profiles with the peak expression at the dark/light transition period (ZT0, Figure 3A). In five out of twelve of the circadian genes inspected (LHY, LCL1, RVE2, RVE8, and CKB3) we identified at least two alternative splicing events (Table 1). The combinatorial nature of the alternative splicing events in these genes may yield transcripts that harbor both early and late PTCs.
Alternative splicing of 3’ untranslated regions
As initially suggested by RNA-seq read data (Additional file 1: Figure S1) and further validated by RT-PCR (Table 1), LHY, LCL1, and RVE8/LCL5 mRNAs had IR events in their 3’ untranslated regions (3’ UTRs). A time course of LHY transcript accumulation showed that isoforms with a retained intron in the 3’ UTR (I9R) and the fully spliced form (I9S) peaked at the dark/light transition period (ZT0, Figure 3B). The retention of introns in the 3’ UTRs of LCL1 (I8R), LHY (I9R), and RVE8 (I8R) results in transcript variants with abnormally long (500–800 nucleotides) 3’ UTRs that may be detected and degraded by NMD. Indeed, RNA-seq read coverage of the LHY 3' UTR also suggested that the levels of the LHY I9R isoform are increased in the NMD-impaired upf1-1 mutant (Additional file 1: Figure S1A).
An early PTC is introduced in the RVE2 pre-mRNA via a cassette exon
Alternative splicing of circadian clock-output genes
Interrogation of AS events using primers designed to hybridize to expected splice junctions confirmed that transcripts of CIRCADIAN RHYTHM AND RNA BINDING2 (CCR2, also known as GLYCINE RICH PROTEIN7 or GRP7) and GLYCINE RICH PROTEIN 8 (CCR1 or GRP8) are also alternatively spliced as suggested by previous studies [9–11]. Validation of RNA-seq data using RT-PCR and primers to splice junction regions showed that some GRP7 transcripts retained a portion of the first intron through the selection of an alternative donor splice site (Table 1 and Additional file 2). We confirmed that a similar partial retention event also occurred in the first intron of the GRP8 transcript through use of an alternative donor splice site. These two analogous events led to incorporation of early PTCs in both GRP7 and GRP8 transcripts (data not shown) and offer another example of conservation of specific AS events among homologs with overlapping but distinct functions.
Nonsense transcript isoforms can accumulate to substantial levels depending on environmental conditions
Differential RNA-seq coverage of introns suggests additional intron-retention events in key circadian genes
To identify additional intron retention events in circadian genes we subjected Arabidopsis seedlings to heat stress treatments followed by RNA-seq analysis as previously described . The RNA-seq dataset derived from the heat-treated plants was compared to the data from untreated controls. The statistical significance of differential expression of the introns of TAIR-annotated genes was evaluated using the GENE-Counter package  as described in the Methods. The differential coverage of introns by the 40-mer Illumina reads before and after heat stress treatment suggested that in addition to already validated transcripts (Table 1) additional IR events may occur in other circadian clock genes; these genes are listed in Additional file 3. The most notable of these potentially alternatively splice genes are the key circadian genes TIME FOR COFFEE (TIC), PHYTOCLOCK 1 (PCL1), and LOV KELCH PROTEIN 2 (LKP2). Therefore, the repertoire of unproductive alternative splicing events in circadian clock genes is likely to be even more diverse than those validated here. The relative levels of these putative novel isoforms with retained introns are likely to change in response to the thermal stress.
Our survey of the Arabidopsis thaliana transcriptome revealed that at least twelve genes associated with the circadian clock are subject to extensive alternative splicing. Three circadian genes harbored more than one AS event, suggesting a complex combinatorial pattern of resulting transcript isoforms. Notably, several AS events were conserved between homologous genes with overlapping but distinct functions. The specific alternative splicing events in some of the core clock genes such as CCA1 are conserved across distantly related plant species . Both CCA1 and LHY are morning-expressed transcription factors that act as negative regulators of TIMING OF CAB EXPRESSION 1 (TOC1) in the core oscillator feedback loop (reviewed in ). CCA1 and LHY repress TOC1 by binding to the EE motif within its promoter [4, 12, 13]). In turn, evening-expressed TOC1 activates the expression of CCA1 and LHY. The conservation of the CCA1 intron 4 retention event in different phyla suggests a selective evolutionary pressure to preserve analogous alternative splicing events in certain circadian genes . In this study we demonstrated that LHY, a CCA1 homolog with overlapping but distinct circadian clock functions, also was alternatively spliced through an analogous I4R event.
Further, IR events were observed in the 3’ UTRs of transcripts from LHY LCL1, and RVE8 genes. Because 3’ UTRs longer than 300 nt can elicit an NMD response in Arabidopsis it is likely that the LHY LCL1, and RVE8 transcripts harboring IR events in their 3' UTRs are targeted by NMD. Indeed, the distribution of RNA-seq read coverage of the 3' UTRs in the NMD-impaired mutant upf1-1 indicated that the LHY-I9R isoform was present at higher levels than those observed in the wild-type strain, and it is therefore likely that this isoform is degraded through the NMD pathway.
The GRP7/GRP8 pair of homologs further support the notion that that similar AS events introducing in-frame PTCs are conserved among homologous but functionally non-redundant circadian genes. The circadian slave oscillator GRP7 modulates clock output networks [4, 10], regulates mRNA export from the nucleus to the cytoplasm during cold stress , and mediates pathogen defense responses . GRP7 auto-regulates its abundance via a negative post-transcriptional feedback loop; it binds to its own pre-mRNA and influences splice site selection [9, 10]. Recently, it was proposed that the reciprocal regulation of GRP7 and GRP8 occurs via an interlocked feedback loop . Moreover, the alternative splicing of GRP7 and GRP8 pre-mRNAs influences circadian oscillations of both transcripts and could be coupled to the NMD pathway. The examples described above suggest certain unproductive AS events in circadian genes have been under stabilizing selection, and therefore it is likely that these transcripts have roles in regulation of the circadian clock.
The functional significance of CCA1 and LHY nonsense transcripts is not known. It is possible that PTC-harboring CCA1 and LHY mRNAs are rapidly degraded by NMD machinery. Therefore, the PTC-containing isoforms of CCA1 and LHY may contribute to regulation of the core feedback loop by fine-tuning the abundance of the full-length protein-coding transcripts. Another possibility is that transcripts harboring early PTCs in the 5’ region of the mRNA escape NMD detection and enter the steady-state translation pool; this occurs in mammalian cells . In this scenario truncated CCA1 and LHY proteins carrying DNA binding MYB domains may compete with the full-length proteins for binding to cis regulatory sites. Because CCA1 and LHY form homo- and heterodimers in vivo, truncated proteins may also disrupt the assembly of functional dimers and may act as dominant negative regulators.
We demonstrate that alternative splicing via poison cassette exon leads to incorporation of early PTCs in RVE2 transcripts. Conserved PCE events have been described in the mRNAs of several human SR splicing factors . The introduction of PTCs in mRNAs through PCE events was proposed to play a role in homeostatic regulation of abundance of the productive (i.e., encoding a full-length protein) transcript . Previously, we characterized an analogous PCE event in the mRNA of the Arabidopsis splicing factor SR34 . However, to our knowledge, the PTC-containing isoform of RVE2 mRNA is the first example of a characterized PCE event that introduces a nonsense codon into the mRNA of a key plant circadian clock gene. This finding suggests that the PCE-type AS events that occur in mammalian SR splicing factors genes  may represent a widespread strategy for introducing in-frame PTCs that is common among eukaryotic organisms including plants.
Cellular localization, recognition by NMD machinery, and translational activity of nonsense circadian transcripts remain unknown. According to the “stochastic noise” hypothesis many unproductive AS events may represent a "cellular noise" of an inherently error prone splicing machinery . However, the “noise” hypothesis cannot explain the fact that unproductive AS is widespread among eukaryotes: Up to 10 % of all mRNAs in yeast, fly, and human cells harbor PTCs . A “noisy” hypothesis is also in disagreement with the notion that many alternative splicing events (including those described here for several circadian genes) are highly conserved between homologs and/or across distantly related species [2, 8].
The final ratio of transcript isoforms is dictated by the choice of splicing signals and is determined by the selection and concentration of specific splicing factors with antagonistic or overlapping functions. The result of an accumulation of alternatively spliced transcripts could be a functional down regulation of the expression of a gene under specific conditions. Conserved alternative splicing patterns in plant SR splicing factors transcripts can be activated developmentally or by changes in external conditions such as environmental stresses [20–22]. This notion is in line with our finding that the level of the PTC-containing CCA1 transcript sharply decreased in cold-treated plants suggesting that unproductive splicing may be involved in the clock temperature compensation mechanism. Indeed, because the levels of CCA1 transcripts increased at lower temperatures, CCA1 was proposed to be responsible for low temperature compensation of the circadian clock .
Although the destiny and regulatory significance of alternatively spliced PTC-harboring transcripts in general and in circadian genes in particular remain unclear, the discovery of widespread unproductive AS in plant circadian genes raises a question of whether alternatively spliced mRNA serve as regulators of the circadian clock. It is possible that unproductively spliced mRNAs function compete with full-length mRNAs in a post-transcriptional mechanism for regulation of protein expression; these PTC-containing transcripts may eventually enter the NMD pathway.
Whether alternatively spliced isoforms can be translated normally or are recognized and targeted by NMD appears to depend on the location of the PTC within the mRNA. Indeed, recent evidence suggests that some transcripts with late PTCs (those in the 3’ region of the message) are translated to produce truncated proteins that act as dominant negative regulators. The examples of such regulation have been previously characterized biochemically for hnRNP A/B proteins  and SR splicing factors . In plants, early PTCs can either escape NMD  or elicit an NMD response less efficiently than the PTCs located in the middle of the coding sequence . We hypothesize that the transcripts harboring relatively early PTCs (such as CCA1-I4R and LHY-I4R) enter the translating pool and the resulting truncated proteins compete with and/or antagonize the normal functions of their full-length counterparts.
Increasing evidence suggests that the coupling of alternative splicing and NMD is a common phenomenon in mammalian cells and plays an important role in regulation of gene expression via regulated unproductive splicing and translation (RUST) [8, 28]. According to the RUST hypothesis, the abundance of a full-length protein-encoding transcript can be regulated by shifting the pre-mRNA splicing towards PTC-containing isoforms that are degraded by the NMD pathway. Because RUST occurs at the post-transcriptional level, it may provide a rapid response to the cell's changing requirements for particular mRNAs. Several instances of modulating productive transcript levels via unproductive alternative splicing and a negative feedback loop have already been described. For example, the mammalian splicing factor SC35 self-regulates its expression by modulating levels of unproductive splicing of the 3’ UTR of its own pre-mRNA .
We propose that RUST regulation of the PTC-containing Arabidopsis core circadian gene transcripts could be of particular importance because it allows rapid post-transcriptional adjustments in abundance of the oscillating productive transcript isoforms over the course of day. This hypothesis is in line with the recent discovery that AS of the core clock gene PRR9 is regulated by PROTEIN ARGININE METHYL TRANSFERASE 5 (PRMT5). Mutations in PRMT5 (which transfers methyl groups to arginine residues of certain spliceosomal proteins) impair circadian rhythms in Arabidopsis by causing severe changes in levels of PRR9 splice isoforms . PRMT5 expression is regulated by light/dark cycles, and it is a likely candidate for the factor that links environmental changes to levels of unproductive spliced products of plant circadian clock genes.
Many plant circadian clock genes undergo extensive alternative splicing that impacts the sequence of mRNA in both coding and untranslated regions. Most but not all of the alternatively spliced isoforms arise from full or partial intron retention events. In different circadian genes, the unproductive alternative splicing can introduce early, intermediate, or late in-frame premature termination codons. The relatively high abundance of transcripts that contain premature termination codons and the conservation of several characterized splicing events in homologous proteins of Arabidopsis and in multiple species suggest that the alternative splicing may play a role in regulation of the cyclical daily expression of plant circadian genes.
Plant growth conditions and sampling scheme
Wild-type control and upf1-1 mutant (ABRC stock CS870823, http://www.arabidopsis.org/) Arabidopsis thaliana (ecotype Columbia 0) seedlings were germinated and grown on plates containing Murashige and Skoog medium for two weeks. The seedlings were subjected to abiotic stress treatments such as high-intensity light, heat, cold, high salinity, and dehydration as previously described . For quantitative RT-PCR analysis, four-week-old Arabidopsis plants were entrained for ten days at constant temperature (20 °C) in a PGR15 growth chamber (Conviron) with day/night cycles of 12 hours light and 12 hours dark and a light intensity of 250 μmol m-2 s-1. For the diurnal profiling of splice variants (as shown in Figure 3) leaf tissues were collected every 3 hours. For the cold treatments (as shown in Figure 5) plants were first entrained as described above. The first time point (ZT0, 24 hours, Figure 5) was collected at 20 °C. Then, the temperature was decreased to 4 °C during the nights. The second time point (48 hours) was collected after exposing plants to 4 °C for 12 hours in the dark. The final time point (168 hours) was collected after the plants were entrained for an additional 96 hours at 4 °C during the night.
Seedlings or leaf tissues were frozen in liquid nitrogen and ground to a fine powder using mortar and pestle. Total RNA was isolated using modification of previously described protocol . Briefly, RNA was extracted using The Plant RNA Reagent (Invitrogen) and treated for 10 min at 65 °C with RNAsecure reagent (Ambion). To eliminate genomic DNA amplification, the RNA was treated for 15 min at 37 °C with RNase-free Turbo DNase (Ambion). Total RNA was additionally purified using RNAeasy Mini RNA kit (Qiagen) following the manufacturer’s RNA clean up protocol. Integrity and concentration of the RNA were monitored using Bioanalyzer 2100 (Agilent Technologies) and ND-1000 spectrophotometer (Thermo Fisher Scientific). Only RNA samples with RNA integrity numbers greater than 6 was used for cDNA synthesis.
Identification of potential alternative splicing events
The cycling profiles of Arabidopsis circadian genes were downloaded using Diurnal Search Tool through the Diurnal portal (http://diurnal.cgrb.oregonstate.edu/). The gene models and RNA-seq coverage among the key circadian clock genes were visually inspected using Arabidopsis RNA-seq GBrowse (http://athal.cgrb.oregonstate.edu/cgi-bin/gbrowse/arabidopsis-gbrowse/). The nucleotide sequences of candidate genes were obtained from Arabidopsis TAIR10 annotation (http://www.arabidopsis.org/) and splice junction specific primers were designed as outlined below.
Primer design and validation of isoform-specific AS events
A strategy for designing splice junction-specific primers was described previously ; examples are depicted in Figure 2. One of the primers was designed to hybridize to the splice junction of interest, the other to intron or exon sequence up- or downstream. The final optimization of primer pairs sequence was performed by using the Primer3 tool (http://primer3.sourceforge.net/). The sequences of oligonucleotide splice isoform-specific primers are listed in Additional file 2.
cDNA synthesis and RT-PCR
The first cDNA strand was synthesized using 100 ng of total mRNA, an anchored oligo d(T) primer, and Superscript III reverse transcriptase using the Invitrogen first-strand cDNA synthesis kit according to the manufacturer’s protocol. The first-strand cDNA reaction was diluted 10-fold and approximately 50 ng of cDNA was used as PCR template. The PCR amplification was carried out using DNA Engine Dyad (Bio-Rad) instrument and Phusion polymerase (New England BioLabs) in the presence of 2 % (v/v) dimethylsulfoxide. The generic “touch down” PCR protocol included 7 “touch down” cycles for the increased specificity (98 °C for 10 sec, 67 °C for 30 sec, 72 °C for 1 min, with annealing temperature decreasing 2 °C at each consequent cycle) followed by 20 to 30 cycles of 95 °C for 15 sec, 57–60 °C for 30 sec, and 72 °C for 1 min and a final extension at 72 °C for 10 min. The number of amplification cycles and annealing temperatures were optimized for each target cDNA depending on transcript abundance and melting temperatures of primer/target duplexes, respectively. Detailed descriptions of individual PCR protocols used for each target are summarized in Additional file 2. The PCR products were separated in 2.5 % agarose gels and stained with ethidium bromide. GOG (AT5G11980) and EF1α (AT5G60390) were used as internal housekeeping mRNA references to ensure an equal RT-PCR amplification and gel loading. Contamination by genomic DNA was monitored by amplification of the genomic ACTIN2 (AT3G18780) sequence.
qRT-PCR was performed using a Bio-Rad CFX96 Real-Time instrument and QuantiTect SYBR Green PCR master mix (Qiagen). The amplification conditions for all qRT-PCR were: 95 °C for 15 sec; then 94 °C for 15 sec, 56 °C for 30 sec (at a slow ramp rate 2 °C per sec), 70 °C for 20 sec (plate read), 75 °C for 15 sec (plate read) for 40 cycles; followed by melting curve step of 65 °C to 95 °C in increments of 0.5 °C (plate read). The normalized fold expression was calculated using CFX Manager software (Bio-Rad) and the -ΔΔC(t) method. GOG (AT5G11980) was used as a housekeeping reference gene for normalization of expression. The expected sizes of all qRT-PCR products were additionally confirmed by gel electrophoresis in 2 % agarose followed by staining with ethidium bromide. The sequences of oligonucleotide splice isoform-specific primers used in qRT-PCR are listed in Additional File 1.
Sanger sequencing was performed according to standard protocols using an ABI 3730 capillary sequencer (Applied Biosystems).
Statistical analysis of the intron coverage by RNA-seq reads
RNA-seq datasets derived from the heat stress treated and untreated control seedlings. The statistical significance of differential expression of the TAIR-annotated introns was evaluated using GENE-Counter tool . The introns were considered differentially expressed at the level of statistical significance if the P and Q values for normalized Illumina reads coverage were less or equal 0.05 (see Additional file 3). An additional filtering included for the fold change in Illumina read coverage (all the introns with fold coverage change less than two were removed from the list). Circadian genes were identified based on TAIR10 annotation (http://www.arabidopsis.org/).
Reviewer 1: Dr. Eugene V. Koonin
Reviewer's comments: Very interesting, solid work. The extensive alternative splicing yielding nonsense mRNAs of clock genes is demonstrated beyond doubt. The direct involvement of this phenomenon in clock regulation remains a hypothesis, and the authors are quite careful about. The only substantial question that I have is: would it be possible to compare the level of alternative splicing in clock genes to the overall background in Arabidopsis? Is it the case that the clock genes are particularly prone to alternative splicing in general and to accumulation of PTC in particular? An answer in the affirmative would extremely interesting and supportive of the regulation hypothesis but even a negative answer would be interesting.
Authors' response: Data from our previous transcriptome-wide RNA-seq study [REF  was incorporated into the Arabidopsis TAIR10 genome annotation and estimates that approximately 42 % of Arabidopsis genes are alternatively spliced. Currently, there are about 29 characterized core clock and circadian-associated genes in the TAIR10 annotation. Approximately half of them (15) are represented by two or more gene models. Therefore, at first glance the extent of alternative splicing in circadian genes appears similar to transcriptome-wide estimates. However, in this study we did not seek to explore statistics of alternative splicing among circadian genes per se. More importantly, we found that the patterns of alternative splicing in some key circadian genes can be very complex and in most cases generates nonsense isoforms. Moreover, circadian genes have a broad repertoire of alternative splicing including intron retention, alternative donor/acceptor sites, cassette exons, and alternatively spliced 5' and 3' UTRs. Thus, the central components of the clock oscillator CCA1 and LHY are represented by two and five gene models, respectively. Interestingly, we found that only one intron 4 retention event (I4R) introducing PTCs downstream of their DNA binding domains is conserved between CCA1 and LHY homologs. In addition, both the 5' and 3' UTRs of LHY are alternatively spliced suggesting that some of its isoforms may be regulated by NMD. Indeed, after this paper was submitted we found that the LHY isoform retaining an intron in its 3' UTR is up-regulated in NMD mutant (data not shown).
Concerning the question of whether the alternative splicing of the circadian genes generates more PTCs than average across the transcriptome, we have no statistically meaningful answer. However, the notion is that most of the alternative splicing events in circadian genes introduce PTCs and these events are likely to be at least as frequent as they are on average across the transcriptome. Not all PTCs will elicit mRNA degradation. Currently we are trying to determine which of the PTCs in circadian genes elicit NMD.
Reviewer 2: Dr. Chungoo Park
Reviewer's comments: In this manuscript, the authors showed that Arabidopsis thaliana circadian clock genes tend to be extensively alternatively spliced, and such isoforms especially retaining intron are subject to contain premature termination codons and to be conserved. However, these findings have been already observed by Filichkin et al. (2010) [Ref . In this study, which is different from previous study?
Authors' response: Even though the results reported in this paper are broadly related to our previous studies they are clearly unique. Further, this study addresses a different set of fundamental biological questions related to the role of unproductive alternative splicing (UAS) in circadian clock regulation for the following reasons.
First, unlike our previous study, here we interrogated novel splice junctions (SJs) in a subset of circadian gene transcripts at specific times of day. To pinpoint the peak expression phase for each transcript we used our recently obtained diurnal microarray datasets [REF . In contrast to the pooled sampling scheme this approach added an important advantage of enabling detection of minor splice variants at peak expression times, which can be missed if timing is not correct.
Second, using a diurnal time course sampling we showed that splice isoforms harboring in-frame premature termination codons (PTCs) can accumulate in a cyclical manner and in many (but not all) instances mirror the oscillations of the productive transcripts. Importantly, our time course data suggests a temporal shift in isoform ratios depending on the time of day or environmental conditions such as temperature changes.
Third, here we demonstrated that the conservation of similar AS events can be observed among circadian-regulated homologs within the same species. In the previous study [REF  we experimentally validated the majority of novel SJs predicted by RNA-seq for hundreds of randomly selected genes. Among these genes, CCA1, a central component of the circadian clock presented an interesting case because its intron 4 retention event was conserved across different plant phyla. Here we show that such conservation can be extended to the homologs within the same species with overlapping but distinct functions such as in the cases of CCA1/LHY and GRP7/GRP8. This notion further reinforces a hypothesis of the functional significance of UAS.
Finally, we showed that some circadian-associated genes such as RVE2 employ a strategy of introducing a PTC via a poison cassette exon (PCE) – a conserved UAS event previously shown to be involved in homeostatic regulation of some mammalian SR splicing factors [REF . To our knowledge this is the first example of such an AS event among circadian genes suggesting that the inclusion of a PTC via a PCE event may be a widespread mechanism, regulating in this instance homeostatic expression of a key circadian gene. Importantly, a time course analysis of the accumulation of RVE2 transcripts (Figure 3C) showed that the PCE isoform may not always mirror precisely the accumulation of its full-length protein-coding counterpart and suggests a possibility that the production of the PTC + transcript may undergo a phase shift under specific diurnal conditions. It is possible that similar to the case of some SR splicing factors, RVE2 expression is regulated by coupling of UAS with mRNA decay.
Combined, the approaches described above allowed us to detect novel AS events in circadian genes which otherwise could fall below statistical significance in RNA-seq analysis because of low isoform abundance at specific times of day and/or environmental conditions resulting in insufficient read coverage. A systematic interrogation of the putative novel splicing events allowed us to validate experimentally that AS of key circadian genes in most cases generates isoforms harboring PTCs. Future studies that explicitly demonstrate regulatory role(s) of UAS in clock function are likely to have major biological implications because the circadian clock orchestrates global gene expression in plants.
Reviewer's comments: One of main assertions in this study is that most circadian clock genes have alternative splicing forms. Whether this pattern is circadian gene-specific or most plant genes are extensively alternatively spliced is unclear.
Authors' response: Please see our response to the similar question from reviewer #1.
Reviewer's comments: To rule out noise hypothesis, abundance of unproductive AS among eukaryotes and conservation of alternative splicing events in relative species were mentioned. To this end, the authors should show direct evidences using plants and circadian genes tested.
Authors' response: The "stochastic splicing noise” versus regulated unproductive AS hypothesis is mentioned in the Discussion section of our manuscript. Rapidly accumulating evidence [for examples see references [8, 18, 19] indicates that unproductive alternative splicing (UAS) coupled with NMD is a regulatory gene expression mechanism that is widespread across eukaryotes. Even though our study suggests that UAS could be an important regulatory mechanism for several circadian genes it is not designed to directly resolve this overarching biological question. The references with the examples supporting "conservation of AS events between different species" have been added to the Discussion section.
Reviewer 3: Dr. Marcelo Yanovsky
Reviewer's comments: Recent work has provided evidence of widespread occurrence of alternative splicing (AS) in plants. Interestingly, genetic approaches have revealed an important role for AS in the proper regulation of circadian rhythms in plants. The manuscript by Filichkin and Mockler contributes to expand our knowledge of the interplay between AS and circadian networks through a thorough evaluation of AS of plant circadian clock genes using RNA-seq. The authors found that many circadian clock-associated genes in Arabidopsis thaliana are alternatively spliced. The results from RNA-seq data were confirmed using several approaches including RT-PCR, qRT-PCR and/or Sanger sequencing. Most AS events led to the incorporation of in-frame premature stop codons (PTC) through full or partial intron retention. This was a widespread phenomenon among mRNAs from the CCA1/LHY family of MYB transcription factors. Interestingly, an in-frame PTC also resulted from the inclusion of a "poison cassette exon" in one of the homologues of CCA1 known as REVEILLE 2 (RVE2). Intron retention events in the MYB family of transcription factors associated with the clock also took place in the 3´utr region, suggesting they may contribute to regulate mRNA stability. Interestingly, daily oscillations with different phases were observed for the different isoforms of the RVE2 gene. Finally, the authors confirmed and extended previous analysis of the effect of stress treatments on AS of CCA1, showing that the full length transcripts increases in response to cold treatments, whilst the isoform retaining intron 4 increases in response to heat, and these responses are sustained under prolonged stress conditions. The dependence of relative abundance of certain isoforms of clock associated genes on time of day or stress treatments suggest that AS of these genes may contribute to fine-tuning the regulation of physiological processes by the clock, as well as regulation of the clock in response to environmental signals.
Reviewer's comments: The manuscript is well written, clear and concise. The data is interesting and presented in a timely fashion for the field, reinforcing the importance of the interplay between AS and the regulation of circadian networks, through changes in AS of core-clock and clock-output genes.
Reviewer's comments: As far as data presentation is concerned, it would have been nice to see as part of the main figures, in addition to read coverage graphs that allow the identification of intron retention or exon skipping events, graphic data on exon/exon junctions, which should allow the identification of alternative donor and/or acceptor sites.
Authors' response: The graphical illustrations and criteria for coverage of exon/exon junctions by Illumina reads used for the identification of alternative donor and/or acceptor sites are shown in Supplementary Figure 7 of our previous publication [Ref .
Reviewer's comments: Finally, in the near future, it will be important to start testing the functional significance of the different AS isoforms of the clock genes, and the role of AS in the proper physiological regulation of the circadian network.
Authors' response: Indeed, this is very important point. Unproductive alternative splicing (UAS) in circadian clock genes presents a very special case study because the daily oscillations of transcripts levels require rapid and steep adjustments in mRNA levels. Currently, studies of the role of AS in regulation of the circadian network are under way in many labs worldwide using an array of mutant and transgenic lines.
Reviewer's comments: The authors mention at the beginning that “To maximize the sensitivity for minor spliced products and to avoid a possible light/dark-dependent shift in isoform ratios, the sampling was done at the phase of peak expression for each tested genes as shown in Figure 1”. Although this approach may simplify the analysis, I do not think it will maximize the sensitivity for minor spliced products. It is possible for instance that the relative abundance of different isoforms changes throughout the day, as shown by the authors for RVE2, and therefore this approach may lead to reduced rather than increased sensitivity.
Authors' response: We found that this approach works in most cases but we agree that for some genes under specific conditions the effect could be opposite: the unproductive isoform may exhibit different behavior than its productive counterpart. This is why we conducted a time course with limited points here (i.e. RVE2) and employ a full time course strategy in the follow up studies.
Reviewer's comments: Indeed, I believe the present description of AS events among plant clock associated genes underestimates the total number of events. For example, additional isoforms of PRR9 have been reported, besides those associated with retention of intron 3, which were not detected/reported here.
Authors' response: We agree and clearly state in the paper that our study is likely to underestimate the total number of AS and UAS events in circadian genes. The Supplemental Figure 3, for example, demonstrates that differential intron expression analysis suggests additional intron retention events under heat stress. Undoubtedly, future studies of abiotic or biotic stress treatments, or specific mutations in splicing-related genes will reveal more novel AS patterns.
Reviewer's comments: Finally, in Figure 2B, the authors should add the primer that allows amplification of the fully spliced isoform, in the 3´-5´direction.
Authors' response: In Figure 2B the primer pairs used for amplification of the fully spliced isoforms have been added.
Nonsense mediated mRNA decay
Premature termination codon
Poison cassette exon
High throughput RNA sequencing
Reverse transcription polymerase chain reaction, qRT-PCR, quantitative RT-PCR
Regulated unproductive splicing and translation
Serine/arginine splicing factors
This study was supported by the United States Department of Agriculture, National Institute of Food and Agriculture (National Research Initiative Competitive Grant 2008–01077 to TCM and SAF), and by the National Science Foundation (Plant Genome Research Program, Grant DBI-0605240 to TCM). We are grateful to Jeff Chang and Jason Cumbie for critical reading of the manuscript and assistance with the GENE-Counter computational pipeline. We thank Sanjuro Jogdeo for help in searching for potential alternative splicing events using GBrowse and Kate Peremyslova, Aurora Vangard, and Selina Liu for help with RNA isolation.
- Filichkin SA, Breton G, Priest HD, Dharmawardhana P, Jaiswal P, Fox SE, Chory J, Michael TP, Kay S, Mockler TC: Global profiling of rice and poplar transcriptomes highlights key conserved circadian-controlled pathways and cis-regulatory modules. PLoS ONE. 2011, 6: e16907-10.1371/journal.pone.0016907.PubMedPubMed CentralView ArticleGoogle Scholar
- Filichkin SA, Priest HD, Givan SA, Shen R, Bryant DW, Fox SE, Wong WK, Mockler TC: Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Res. 2010, 20: 45-58. 10.1101/gr.093302.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Cumbie JS, Kimbrel JA, Di Y, Schafer DW, Wilhelm L, Fox SE, Sullivan CM, Curzon AD, Carrington JC, Mockler TC, Chang JH: GENE-counter: a computational pipeline for the analysis of RNA-Seq data for gene expression differences. PLoS ONE. 2011, 6: e25279-10.1371/journal.pone.0025279.PubMedPubMed CentralView ArticleGoogle Scholar
- Pruneda-Paza JL, Kay SA: An expanding universe of circadian networks in higher plants. Trends Plant Sci. 2010, 15: 259-265. 10.1016/j.tplants.2010.03.003.View ArticleGoogle Scholar
- Wang BB, Brendel V: Genomewide comparative analysis of alternative splicing in plants. Proc Natl Acad Sci USA. 2006, 103: 7175-7180. 10.1073/pnas.0602039103.PubMedPubMed CentralView ArticleGoogle Scholar
- Yanhui C, Xiaoyuan Y, Kun H, Meihua L, Jigang L, Zhaofeng G, Zhiqiang L, Yunfei Z, Xiaoxiao W, Xiaoming Q, Yunping S, Li Z, Xiaohui D, Jingchu L, Xing-Wang D, Zhangliang C, Hongya G, Li-Jia Q: The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006, 60: 107-124. 10.1007/s11103-005-2910-y.PubMedView ArticleGoogle Scholar
- Zhang X, Chen Y, Wang ZY, Chen Z, Gu H, Qu LJ: Constitutive expression of CIR1 (RVE2) affects several circadian-regulated processes and seed germination in Arabidopsis. Plant J. 2007, 51: 512-25. 10.1111/j.1365-313X.2007.03156.x.PubMedView ArticleGoogle Scholar
- Lareau LF, Inada M, Green RE, Wengrod JC, Brenner SE: Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements. Nature. 2007, 446: 926-929. 10.1038/nature05676.PubMedView ArticleGoogle Scholar
- Staiger D, Zecca L, Wieczorek DA, Apel K, Eckstein L: The clock regulated RNA-binding protein AtGRP7 autoregulates its expression by influencing alternative splicing of its own pre-mRNA. Plant J. 2003, 33: 361-371. 10.1046/j.1365-313X.2003.01629.x.PubMedView ArticleGoogle Scholar
- Schöning JC, Streitner C, Page DR, Hennig S, Uchida K, Wolf E, Furuya M, Staiger D: Autoregulation of the circadian slave oscillator component ATGRP7 and regulation of its targets is impaired by a single RNA recognition motif point mutation. Plant J. 2007, 52: 1119-1130. 10.1111/j.1365-313X.2007.03302.x.PubMedView ArticleGoogle Scholar
- Schöning JC, Streitner C, Meyer IM, Gao Y, Staiger D: Reciprocal regulation of glycine-rich RNA-binding proteins via an interlocked feedback loop coupling alternative splicing to nonsense-mediated decay in Arabidopsis. Nucleic Acids Res. 2008, 36: 6977-6987. 10.1093/nar/gkn847.PubMedPubMed CentralView ArticleGoogle Scholar
- Alabadí D, Oyama T, Yanovsky MJ, Harmon FG, Más P, Kay SA: Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001, 293: 880-883. 10.1126/science.1061320.PubMedView ArticleGoogle Scholar
- Mizoguchi T, Wheatley K, Hanzawa Y, Wright L, Mizoguchi M, Song HR, Carré IA, Coupland G: LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis. Dev Cell. 2001, 2: 629-641.View ArticleGoogle Scholar
- Hori K, Watanabe Y: Context analysis of termination codons in mRNA that are recognized by plant NMD. Plant Cell Physiol. 2007, 48: 1072-1078. 10.1093/pcp/pcm075.PubMedView ArticleGoogle Scholar
- Fu ZQ, Guo M, Jeong BR, Tian F, Elthon TE, Cerny RL, Staiger D, Alfano JR: A type III effector ADP-ribosylates RNA-binding proteins and quells plant immunity. Nature. 2007, 447: 284-288. 10.1038/nature05737.PubMedView ArticleGoogle Scholar
- Inacio A, Silva AL, Pinto J, Ji X, Morgado A, Almeida F, Faustino P, Lavinha J, Liebhaber SA, Romao L: Nonsense mutations in close proximity to the initiation codon fail to trigger full nonsense-mediated mRNA decay. J Biol Chem. 2004, 279: 32170-32180. 10.1074/jbc.M405024200.PubMedView ArticleGoogle Scholar
- Lu SX, Knowles SM, Andronis C, Ong MS, Tobin EM: CIRCADIAN CLOCK ASSOCIATED1 and LATE ELONGATED HYPOCOTYL function synergistically in the circadian clock of Arabidopsis. Plant Physiol. 2009, 150: 834-843. 10.1104/pp.108.133272.PubMedPubMed CentralView ArticleGoogle Scholar
- Saltzman AL, Kim YK, Pan Q, Fagnani MM, Maquat LE, Blencowe BJ: Regulation of multiple core spliceosomal proteins by alternative splicing-coupled nonsense-mediated mRNA decay. Mol Cell Biol. 2008, 28: 4320-4330. 10.1128/MCB.00361-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Rehwinkel J, Raes J, Izaurralde E: Nonsense-mediated mRNA decay: Target genes and functional diversification of effectors. Trends Biochem Sci. 2006, 31: 639-646. 10.1016/j.tibs.2006.09.005.PubMedView ArticleGoogle Scholar
- Ali GS, Palusa SG, Golovkin M, Prasad J, Manley JL, Reddy AS: Regulation of plant developmental processes by a novel splicing factor. PLoS ONE. 2007, 2: e471-10.1371/journal.pone.0000471.PubMedPubMed CentralView ArticleGoogle Scholar
- Palusa SG, Ali GS, Reddy ASN: Alternative splicing of pre-mRNAs of Arabidopsis serine/arginine-rich proteins and its regulation by hormones and stresses. Plant J. 2007, 49: 1091-1107. 10.1111/j.1365-313X.2006.03020.x.PubMedView ArticleGoogle Scholar
- Kalyna M, Lopato S, Voronin V, Barta A: Evolutionary conservation and regulation of particular alternative splicing events in plant SR proteins. Nucleic Acids Res. 2006, 34: 4395-4405. 10.1093/nar/gkl570.PubMedPubMed CentralView ArticleGoogle Scholar
- Gould PD, Locke JC, Larue C, Southern MM, Davis SJ, Hanano S, Moyle R, Milich R, Putterill J, Millar AJ, Hall A: The molecular basis of temperature compensation in the Arabidopsis circadian clock. Plant Cell. 2006, 18: 1177-1187. 10.1105/tpc.105.039990.PubMedPubMed CentralView ArticleGoogle Scholar
- Mayeda A, Munroe SH, Cáceres JF, Krainer AR: Function of conserved domains of hnRNP A1 and other hnRNP A/B proteins. EMBO J. 1994, 15: 5483-5495.Google Scholar
- Zhu J, Krainer AR: Pre-mRNA splicing in the absence of an SR protein RS domain. Genes Dev. 2000, 14: 3166-3178. 10.1101/gad.189500.PubMedPubMed CentralView ArticleGoogle Scholar
- Hori K, Watanabe Y: Context analysis of termination codons in mRNA that are recognized by plant NMD. Plant Cell Physiol. 2007, 48: 1072-1078. 10.1093/pcp/pcm075.PubMedView ArticleGoogle Scholar
- Wu J, Kang JH, Hettenhausen C, Baldwin IT: Nonsense mediated mRNA decay (NMD) silences the accumulation of aberrant trypsin proteinase inhibitor mRNA in Nicotiana attenuata. Plant J. 2007, 51: 693-706. 10.1111/j.1365-313X.2007.03173.x.PubMedView ArticleGoogle Scholar
- Lewis BP, Green RE, Brenner SE: Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc Natl Acad Sci USA. 2003, 100: 189-192. 10.1073/pnas.0136770100.PubMedPubMed CentralView ArticleGoogle Scholar
- Dreumont N, Hardy S, Behm-Ansmant I, Kister L, Branlant C, Stévenin J, Bourgeois CF: Antagonistic factors control the unproductive splicing of SC35 terminal intron. Nucleic Acids Res. 2010, 38: 1353-1366. 10.1093/nar/gkp1086.PubMedPubMed CentralView ArticleGoogle Scholar
- Sanchez SE, Petrillo E, Beckwith EJ, Zhang X, Rugnone ML, Hernando CE, Cuevas JC, Godoy Herz MA, Depetris-Chauvin A, Simpson CG, Brown JW, Cerdán PD, Borevitz JO, Mas P, Ceriani MF, Kornblihtt AR, Yanovsky MJ: A methyl transferase links the circadian clock to the regulation of alternative splicing. Nature. 2010, 468: 112-116. 10.1038/nature09470.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.