Natural variation in Brachypodium distachyon responses to combined abiotic stresses.

The demand for agricultural production is becoming more challenging as climate change increases global temperature and the frequency of extreme weather events. This study examines the phenotypic variation of 149 accessions of Brachypodium distachyon under drought, heat, and the combination of stresses. Heat alone causes the largest amounts of tissue damage while the combination of stresses causes the largest decrease in biomass compared to other treatments. Notably, Bd21-0, the reference line for B. distachyon, did not have robust growth under stress conditions, especially the heat and combined drought and heat treatments. The climate of origin was significantly associated with B. distachyon responses to the assessed stress conditions. Additionally, a GWAS found loci associated with changes in plant height and the amount of damaged tissue under stress. Some of these SNPs were closely located to genes known to be involved in responses to abiotic stresses and point to potential causative loci in plant stress response. However, SNPs found to be significantly associated with a response to heat or drought individually are not also significantly associated with the combination of stresses. This, with the phenotypic data, suggests that the effects of these abiotic stresses are not simply additive, and the responses to the combined stresses differ from drought and heat alone.

Genome and time-of-day transcriptome of Wolffia australiana link morphological minimization with gene loss and less growth control.

Rootless plants in the genus Wolffia are some of the fastest growing known plants on Earth. Wolffia have a reduced body plan, primarily multiplying through a budding type of asexual reproduction. Here, we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas, which has the smallest genome size in the genus at 357 Mb and has a reduced set of predicted protein-coding genes at about 15,000. Comparison between multiple high-quality draft genome sequences from W. australiana clones confirmed loss of several hundred genes that are highly conserved among flowering plants, including genes involved in root developmental and light signaling pathways. Wolffia has also lost most of the conserved nucleotide-binding leucine-rich repeat (NLR) genes that are known to be involved in innate immunity, as well as those involved in terpene biosynthesis, while having a significant overrepresentation of genes in the sphingolipid pathways that may signify an alternative defense system. Diurnal expression analysis revealed that only 13% of Wolffia genes are expressed in a time-of-day (TOD) fashion, which is less than the typical ∼40% found in several model plants under the same condition. In contrast to the model plants Arabidopsis and rice, many of the pathways associated with multicellular and developmental processes are not under TOD control in W. australiana, where genes that cycle the conditions tested predominantly have carbon processing and chloroplast-related functions. The Wolffia genome and TOD expression data set thus provide insight into the interplay between a streamlined plant body plan and optimized growth.

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response.

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

Time of day and network reprogramming during drought induced CAM photosynthesis in Sedum album.

Plants with facultative crassulacean acid metabolism (CAM) maximize performance through utilizing C3 or C4 photosynthesis under ideal conditions while temporally switching to CAM under water stress (drought). While genome-scale analyses of constitutive CAM plants suggest that time of day networks are shifted, or phased to the evening compared to C3, little is known for how the shift from C3 to CAM networks is modulated in drought induced CAM. Here we generate a draft genome for the drought-induced CAM-cycling species Sedum album. Through parallel sampling in well-watered (C3) and drought (CAM) conditions, we uncover a massive rewiring of time of day expression and a CAM and stress-specific network. The core circadian genes are expanded in S. album and under CAM induction, core clock genes either change phase or amplitude. While the core clock cis-elements are conserved in S. album, we uncover a set of novel CAM and stress specific cis-elements consistent with our finding of rewired co-expression networks. We identified shared elements between constitutive CAM and CAM-cycling species and expression patterns unique to CAM-cycling S. album. Together these results demonstrate that drought induced CAM-cycling photosynthesis evolved through the mobilization of a stress-specific, time of day network, and not solely the phasing of existing C3 networks. These results will inform efforts to engineer water use efficiency into crop plants for growth on marginal land.

Single-molecule sequencing of the desiccation-tolerant grass Oropetium thomaeum.

Plant genomes, and eukaryotic genomes in general, are typically repetitive, polyploid and heterozygous, which complicates genome assembly. The short read lengths of early Sanger and current next-generation sequencing platforms hinder assembly through complex repeat regions, and many draft and reference genomes are fragmented, lacking skewed GC and repetitive intergenic sequences, which are gaining importance due to projects like the Encyclopedia of DNA Elements (ENCODE). Here we report the whole-genome sequencing and assembly of the desiccation-tolerant grass Oropetium thomaeum. Using only single-molecule real-time sequencing, which generates long (>16 kilobases) reads with random errors, we assembled 99% (244 megabases) of the Oropetium genome into 625 contigs with an N50 length of 2.4 megabases. Oropetium is an example of a 'near-complete' draft genome which includes gapless coverage over gene space as well as intergenic sequences such as centromeres, telomeres, transposable elements and rRNA clusters that are typically unassembled in draft genomes. Oropetium has 28,466 protein-coding genes and 43% repeat sequences, yet with 30% more compact euchromatic regions it is the smallest known grass genome. The Oropetium genome demonstrates the utility of single-molecule real-time sequencing for assembling high-quality plant and other eukaryotic genomes, and serves as a valuable resource for the plant comparative genomics community.

Architecture and evolution of a minute plant genome.

It has been argued that the evolution of plant genome size is principally unidirectional and increasing owing to the varied action of whole-genome duplications (WGDs) and mobile element proliferation. However, extreme genome size reductions have been reported in the angiosperm family tree. Here we report the sequence of the 82-megabase genome of the carnivorous bladderwort plant Utricularia gibba. Despite its tiny size, the U. gibba genome accommodates a typical number of genes for a plant, with the main difference from other plant genomes arising from a drastic reduction in non-genic DNA. Unexpectedly, we identified at least three rounds of WGD in U. gibba since common ancestry with tomato (Solanum) and grape (Vitis). The compressed architecture of the U. gibba genome indicates that a small fraction of intergenic DNA, with few or no active retrotransposons, is sufficient to regulate and integrate all the processes required for the development and reproduction of a complex organism.

Grasses suppress shoot-borne roots to conserve water during drought.

Many important crops are members of the Poaceae family, which develop root systems characterized by a high degree of root initiation from the belowground basal nodes of the shoot, termed the crown. Although this postembryonic shoot-borne root system represents the major conduit for water uptake, little is known about the effect of water availability on its development. Here we demonstrate that in the model C4 grass Setaria viridis, the crown locally senses water availability and suppresses postemergence crown root growth under a water deficit. This response was observed in field and growth room environments and in all grass species tested. Luminescence-based imaging of root systems grown in soil-like media revealed a shift in root growth from crown-derived to primary root-derived branches, suggesting that primary root-dominated architecture can be induced in S. viridis under certain stress conditions. Crown roots of Zea mays and Setaria italica, domesticated relatives of teosinte and S. viridis, respectively, show reduced sensitivity to water deficit, suggesting that this response might have been influenced by human selection. Enhanced water status of maize mutants lacking crown roots suggests that under a water deficit, stronger suppression of crown roots actually may benefit crop productivity.

Temporal and spatial transcriptomic and microRNA dynamics of CAM photosynthesis in pineapple.

The altered carbon assimilation pathway of crassulacean acid metabolism (CAM) photosynthesis results in an up to 80% higher water-use efficiency than C3 photosynthesis in plants making it a potentially useful pathway for engineering crop plants with improved drought tolerance. Here we surveyed detailed temporal (diel time course) and spatial (across a leaf gradient) gene and microRNA (miRNA) expression patterns in the obligate CAM plant pineapple [Ananas comosus (L.) Merr.]. The high-resolution transcriptome atlas allowed us to distinguish between CAM-related and non-CAM gene copies. A differential gene co-expression network across green and white leaf diel datasets identified genes with circadian oscillation, CAM-related functions, and source-sink relations. Gene co-expression clusters containing CAM pathway genes are enriched with clock-associated cis-elements, suggesting circadian regulation of CAM. About 20% of pineapple microRNAs have diel expression patterns, with several that target key CAM-related genes. Expression and physiology data provide a model for CAM-specific carbohydrate flux and long-distance hexose transport. Together these resources provide a list of candidate genes for targeted engineering of CAM into C3 photosynthesis crop species.

Comprehensive definition of genome features in Spirodela polyrhiza by high-depth physical mapping and short-read DNA sequencing strategies.

Spirodela polyrhiza is a fast-growing aquatic monocot with highly reduced morphology, genome size and number of protein-coding genes. Considering these biological features of Spirodela and its basal position in the monocot lineage, understanding its genome architecture could shed light on plant adaptation and genome evolution. Like many draft genomes, however, the 158-Mb Spirodela genome sequence has not been resolved to chromosomes, and important genome characteristics have not been defined. Here we deployed rapid genome-wide physical maps combined with high-coverage short-read sequencing to resolve the 20 chromosomes of Spirodela and to empirically delineate its genome features. Our data revealed a dramatic reduction in the number of the rDNA repeat units in Spirodela to fewer than 100, which is even fewer than that reported for yeast. Consistent with its unique phylogenetic position, small RNA sequencing revealed 29 Spirodela-specific microRNA, with only two being shared with Elaeis guineensis (oil palm) and Musa balbisiana (banana). Combining DNA methylation data and small RNA sequencing enabled the accurate prediction of 20.5% long terminal repeats (LTRs) that doubled the previous estimate, and revealed a high Solo:Intact LTR ratio of 8.2. Interestingly, we found that Spirodela has the lowest global DNA methylation levels (9%) of any plant species tested. Taken together our results reveal a genome that has undergone reduction, likely through eliminating non-essential protein coding genes, rDNA and LTRs. In addition to delineating the genome features of this unique plant, the methodologies described and large-scale genome resources from this work will enable future evolutionary and functional studies of this basal monocot family.

Climate-smart crops with enhanced photosynthesis.

The potential of enhanced photosynthetic efficiency to help achieve the sustainable yield increases required to meet future demands for food and energy has spurred intense research towards understanding, modeling, and engineering photosynthesis. These current efforts, largely focused on the C3 model Arabidopsis thaliana or crop plants (e.g. rice, sorghum, maize, and wheat), could be intensified and broadened using model systems closely related to our food, feed, and energy crops and that allow rapid design-build-test-learn cycles. In this outlooking Opinion, we advocate for a concerted effort to expand our understanding and improve our ability to redesign carbon uptake, allocation, and utilization. We propose two specific research directions that combine enhanced photosynthesis with climate-smart metabolic attributes: (i) engineering pathways for flexible (facultative) C3-C4 metabolism where plants will operate either C3 or C4 photosynthesis based on environmental conditions such as temperature, light, and atmospheric CO2 levels; and (ii) increasing rhizospheric sink strength for carbon utilization, including strategies that allow for augmented transport of carbon to the soil for improved soil properties and carbon storage without jeopardizing aboveground crop biomass. We argue that such ambitious undertakings be first approached and demonstrated by exploring the full genomic potential of two model grasses, the C3Brachypodium distachyon and the C4Setaria viridis. The development of climate-smart crops could provide novel and bold solutions to increase crop productivity while reducing atmospheric carbon and nitrogen emissions.

Translatome analyses capture of opposing tissue-specific brassinosteroid signals orchestrating root meristem differentiation.

The mechanisms ensuring balanced growth remain a critical question in developmental biology. In plants, this balance relies on spatiotemporal integration of hormonal signaling pathways, but the understanding of the precise contribution of each hormone is just beginning to take form. Brassinosteroid (BR) hormone is shown here to have opposing effects on root meristem size, depending on its site of action. BR is demonstrated to both delay and promote onset of stem cell daughter differentiation, when acting in the outer tissue of the root meristem, the epidermis, and the innermost tissue, the stele, respectively. To understand the molecular basis of this phenomenon, a comprehensive spatiotemporal translatome mapping of Arabidopsis roots was performed. Analyses of wild type and mutants featuring different distributions of BR revealed autonomous, tissue-specific gene responses to BR, implying its contrasting tissue-dependent impact on growth. BR-induced genes were primarily detected in epidermal cells of the basal meristem zone and were enriched by auxin-related genes. In contrast, repressed BR genes prevailed in the stele of the apical meristem zone. Furthermore, auxin was found to mediate the growth-promoting impact of BR signaling originating in the epidermis, whereas BR signaling in the stele buffered this effect. We propose that context-specific BR activity and responses are oppositely interpreted at the organ level, ensuring coherent growth.

The genome of black raspberry (Rubus occidentalis).

Black raspberry (Rubus occidentalis) is an important specialty fruit crop in the US Pacific Northwest that can hybridize with the globally commercialized red raspberry (R. idaeus). Here we report a 243 Mb draft genome of black raspberry that will serve as a useful reference for the Rosaceae and Rubus fruit crops (raspberry, blackberry, and their hybrids). The black raspberry genome is largely collinear to the diploid woodland strawberry (Fragaria vesca) with a conserved karyotype and few notable structural rearrangements. Centromeric satellite repeats are widely dispersed across the black raspberry genome, in contrast to the tight association with the centromere observed in most plants. Among the 28 005 predicted protein-coding genes, we identified 290 very recent small-scale gene duplicates enriched for sugar metabolism, fruit development, and anthocyanin related genes which may be related to key agronomic traits during black raspberry domestication. This contrasts patterns of recent duplications in the wild woodland strawberry F. vesca, which show no patterns of enrichment, suggesting gene duplications contributed to domestication traits. Expression profiles from a fruit ripening series and roots exposed to Verticillium dahliae shed insight into fruit development and disease response, respectively. The resources presented here will expedite the development of improved black and red raspberry, blackberry and other Rubus cultivars.

Seed desiccation mechanisms co-opted for vegetative desiccation in the resurrection grass Oropetium thomaeum.

Resurrection plants desiccate during periods of prolonged drought stress, then resume normal cellular metabolism upon water availability. Desiccation tolerance has multiple origins in flowering plants, and it likely evolved through rewiring seed desiccation pathways. Oropetium thomaeum is an emerging model for extreme drought tolerance, and its genome, which is the smallest among surveyed grasses, was recently sequenced. Combining RNA-seq, targeted metabolite analysis and comparative genomics, we show evidence for co-option of seed-specific pathways during vegetative desiccation. Desiccation-related gene co-expression clusters are enriched in functions related to seed development including several seed-specific transcription factors. Across the metabolic network, pathways involved in programmed cell death inhibition, ABA signalling and others are activated during dehydration. Oleosins and oil bodies that typically function in seed storage are highly abundant in desiccated leaves and may function for membrane stability and storage. Orthologs to seed-specific LEA proteins from rice and maize have neofunctionalized in Oropetium with high expression during desiccation. Accumulation of sucrose, raffinose and stachyose in drying leaves mirrors sugar accumulation patterns in maturing seeds. Together, these results connect vegetative desiccation with existing seed desiccation and drought responsive pathways and provide some key candidate genes for engineering improved drought tolerance in crop plants.

Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA precursors.

Artificial microRNAs (amiRNAs) are used for selective gene silencing in plants. However, current methods to produce amiRNA constructs for silencing transcripts in monocot species are not suitable for simple, cost-effective and large-scale synthesis. Here, a series of expression vectors based on Oryza sativa MIR390 (OsMIR390) precursor was developed for high-throughput cloning and high expression of amiRNAs in monocots. Four different amiRNA sequences designed to target specifically endogenous genes and expressed from OsMIR390-based vectors were validated in transgenic Brachypodium distachyon plants. Surprisingly, amiRNAs accumulated to higher levels and were processed more accurately when expressed from chimeric OsMIR390-based precursors that include distal stem-loop sequences from Arabidopsis thaliana MIR390a (AtMIR390a). In all cases, transgenic plants displayed the predicted phenotypes induced by target gene repression, and accumulated high levels of amiRNAs and low levels of the corresponding target transcripts. Genome-wide transcriptome profiling combined with 5'-RLM-RACE analysis in transgenic plants confirmed that amiRNAs were highly specific.

Subset of heat-shock transcription factors required for the early response of Arabidopsis to excess light.

Sunlight provides energy for photosynthesis and is essential for nearly all life on earth. However, too much or too little light or rapidly fluctuating light conditions cause stress to plants. Rapid changes in the amount of light are perceived as a change in the reduced/oxidized (redox) state of photosynthetic electron transport components in chloroplasts. However, how this generates a signal that is relayed to changes in nuclear gene expression is not well understood. We modified redox state in the reference plant, Arabidopsis thaliana, using either excess light or low light plus the herbicide DBMIB (2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone), a well-known inhibitor of photosynthetic electron transport. Modification of redox state caused a change in expression of a common set of about 750 genes, many of which are known stress-responsive genes. Among the most highly enriched promoter elements in the induced gene set were heat-shock elements (HSEs), known motifs that change gene expression in response to high temperature in many systems. We show that HSEs from the promoter of the ASCORBATE PEROXIDASE 2 (APX2) gene were necessary and sufficient for APX2 expression in conditions of excess light, or under low light plus the herbicide. We tested APX2 expression phenotypes in overexpression and loss-of-function mutants of 15 Arabidopsis A-type heat-shock transcription factors (HSFs), and identified HSFA1D, HSFA2, and HSFA3 as key factors regulating APX2 expression in diverse stress conditions. Excess light regulates both the subcellular location of HSFA1D and its biochemical properties, making it a key early component of the excess light stress network of plants.

Genome diversity in Brachypodium distachyon: deep sequencing of highly diverse inbred lines.

Brachypodium distachyon is small annual grass that has been adopted as a model for the grasses. Its small genome, high-quality reference genome, large germplasm collection, and selfing nature make it an excellent subject for studies of natural variation. We sequenced six divergent lines to identify a comprehensive set of polymorphisms and analyze their distribution and concordance with gene expression. Multiple methods and controls were utilized to identify polymorphisms and validate their quality. mRNA-Seq experiments under control and simulated drought-stress conditions, identified 300 genes with a genotype-dependent treatment response. We showed that large-scale sequence variants had extremely high concordance with altered expression of hundreds of genes, including many with genotype-dependent treatment responses. We generated a deep mRNA-Seq dataset for the most divergent line and created a de novo transcriptome assembly. This led to the discovery of >2400 previously unannotated transcripts and hundreds of genes not present in the reference genome. We built a public database for visualization and investigation of sequence variants among these widely used inbred lines.

A genetic linkage map of black raspberry (Rubus occidentalis) and the mapping of Ag(4) conferring resistance to the aphid Amphorophora agathonica.

KEY MESSAGE: We have constructed a densely populated, saturated genetic linkage map of black raspberry and successfully placed a locus for aphid resistance. Black raspberry (Rubus occidentalis L.) is a high-value crop in the Pacific Northwest of North America with an international marketplace. Few genetic resources are readily available and little improvement has been achieved through breeding efforts to address production challenges involved in growing this crop. Contributing to its lack of improvement is low genetic diversity in elite cultivars and an untapped reservoir of genetic diversity from wild germplasm. In the Pacific Northwest, where most production is centered, the current standard commercial cultivar is highly susceptible to the aphid Amphorophora agathonica Hottes, which is a vector for the Raspberry mosaic virus complex. Infection with the virus complex leads to a rapid decline in plant health resulting in field replacement after only 3-4 growing seasons. Sources of aphid resistance have been identified in wild germplasm and are used to develop mapping populations to study the inheritance of these valuable traits. We have constructed a genetic linkage map using single-nucleotide polymorphism and transferable (primarily simple sequence repeat) markers for F1 population ORUS 4305 consisting of 115 progeny that segregate for aphid resistance. Our linkage map of seven linkage groups representing the seven haploid chromosomes of black raspberry consists of 274 markers on the maternal map and 292 markers on the paternal map including a morphological locus for aphid resistance. This is the first linkage map of black raspberry and will aid in developing markers for marker-assisted breeding, comparative mapping with other Rubus species, and enhancing the black raspberry genome assembly.

Cell-autonomous circadian clock of hepatocytes drives rhythms in transcription and polyamine synthesis.

The circadian clock generates daily rhythms in mammalian liver processes, such as glucose and lipid homeostasis, xenobiotic metabolism, and regeneration. The mechanisms governing these rhythms are not well understood, particularly the distinct contributions of the cell-autonomous clock and central pacemaker to rhythmic liver physiology. Through microarray expression profiling in Met murine hepatocytes (MMH)-D3, we identified over 1,000 transcripts that exhibit circadian oscillations, demonstrating that the cell-autonomous clock can drive many rhythms, and that MMH-D3 is a valid circadian model system. The genes represented by these circadian transcripts displayed both cophasic and antiphasic organization within a protein-protein interaction network, suggesting the existence of competition for binding sites or partners by genes of disparate transcriptional phases. Multiple pathways displayed enrichment in MMH-D3 circadian transcripts, including the polyamine synthesis module of the glutathione metabolic pathway. The polyamine synthesis module, which is highly associated with cell proliferation and whose products are required for initiation of liver regeneration, includes enzymes whose transcripts exhibit circadian oscillations, such as ornithine decarboxylase and spermidine synthase. Metabolic profiling revealed that the enzymatic product of spermidine synthase, spermidine, cycles as well. Thus, the cell-autonomous hepatocyte clock can drive a significant amount of transcriptional rhythms and orchestrate physiologically relevant modules such as polyamine synthesis.

Longitudinal genome-wide association study reveals early QTL that predict biomass accumulation under cold stress in sorghum.

Introduction: Sorghum bicolor is a promising cellulosic feedstock crop for bioenergy due to its high biomass yields. However, early growth phases of sorghum are sensitive to cold stress, limiting its planting in temperate environments. Cold adaptability is crucial for cultivating bioenergy and grain sorghum at higher latitudes and elevations, or for extending the growing season. Identifying genes and alleles that enhance biomass accumulation under early cold stress can lead to improved sorghum varieties through breeding or genetic engineering. Methods: We conducted image-based phenotyping on 369 accessions from the sorghum Bioenergy Association Panel (BAP) in a controlled environment with early cold treatment. The BAP includes diverse accessions with dense genotyping and varied racial, geographical, and phenotypic backgrounds. Daily, non-destructive imaging allowed temporal analysis of growth-related traits and water use efficiency (WUE). A genome-wide association study (GWAS) was performed to identify genomic intervals and genes associated with cold stress response. Results: The GWAS identified transient quantitative trait loci (QTL) strongly associated with growth-related traits, enabling an exploration of the genetic basis of cold stress response at different developmental stages. This analysis of daily growth traits, rather than endpoint traits, revealed early transient QTL predictive of final phenotypes. The study identified both known and novel candidate genes associated with growth-related traits and temporal responses to cold stress. Discussion: The identified QTL and candidate genes contribute to understanding the genetic mechanisms underlying sorghum's response to cold stress. These findings can inform breeding and genetic engineering strategies to develop sorghum varieties with improved biomass yields and resilience to cold, facilitating earlier planting, extended growing seasons, and cultivation at higher latitudes and elevations.

Genome-wide mapping of alternative splicing in Arabidopsis thaliana.

Alternative splicing can enhance transcriptome plasticity and proteome diversity. In plants, alternative splicing can be manifested at different developmental stages, and is frequently associated with specific tissue types or environmental conditions such as abiotic stress. We mapped the Arabidopsis transcriptome at single-base resolution using the Illumina platform for ultrahigh-throughput RNA sequencing (RNA-seq). Deep transcriptome sequencing confirmed a majority of annotated introns and identified thousands of novel alternatively spliced mRNA isoforms. Our analysis suggests that at least approximately 42% of intron-containing genes in Arabidopsis are alternatively spliced; this is significantly higher than previous estimates based on cDNA/expressed sequence tag sequencing. Random validation confirmed that novel splice isoforms empirically predicted by RNA-seq can be detected in vivo. Novel introns detected by RNA-seq were substantially enriched in nonconsensus terminal dinucleotide splice signals. Alternative isoforms with premature termination codons (PTCs) comprised the majority of alternatively spliced transcripts. Using an example of an essential circadian clock gene, we show that intron retention can generate relatively abundant PTC(+) isoforms and that this specific event is highly conserved among diverse plant species. Alternatively spliced PTC(+) isoforms can be potentially targeted for degradation by the nonsense mediated mRNA decay (NMD) surveillance machinery or regulate the level of functional transcripts by the mechanism of regulated unproductive splicing and translation (RUST). We demonstrate that the relative ratios of the PTC(+) and reference isoforms for several key regulatory genes can be considerably shifted under abiotic stress treatments. Taken together, our results suggest that like in animals, NMD and RUST may be widespread in plants and may play important roles in regulating gene expression.