Search

Accurate characterization of splice junctions as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data and enable detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset (RTD) from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the splice junctions and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared to an earlier version, BaRTv1.0. The accuracy of transcript level quantification, splice junctions and transcript start and end sites has been validated extensively using parallel technologies and analysis, including high resolution RT PCR and 5’ RACE. BaRTv2.18 contains 39,434 genes and 148,260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

Background The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants. Results A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts – BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427–433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20–28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5′ and 3′ UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage. Conclusion A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley., Measurement(s) gene expression Technology Type(s) transcription profiling assay Factor Type(s) Genotype • Abiotic Stress • Developmental stage • Tissue • Biotic stress Sample Characteristic - Organism Hordeum vulgare Sample Characteristic - Environment sodium chloride salt • drought • increased temperature • decreased temperature Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13643387

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EoRNA* – Expression of RNA) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EoRNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley

Nonsense mediated RNA decay (NMD) is an evolutionary conserved RNA control mechanism that has also been implicated in the broader regulation of gene expression. Nevertheless, a role for NMD in genome regulation has not been fully assessed, partially because NMD inactivation is lethal in many organisms. Here, we performed in depth comparative analysis of Arabidopsis mutants lacking key proteins involved in different steps of NMD. We observed that UPF3, UPF1, and SMG7 have different impacts on NMD and the Arabidopsis transcriptome, with UPF1 having the biggest effect. Transcriptome assembly using stringent pipeline in UPF1-null plants revealed genome wide changes in alternative splicing, including switches in mRNA variants, suggesting a role for UPF1 in splicing. We further found that UPF1 inactivation leads to translational repression, manifested by a global shift in mRNAs from polysomes to monosomes and a downregulation of genes involved in translation and ribosome biogenesis. Despite this global change, NMD targets and low-expressed mRNAs with short half-lives were enriched in polysomes, indicating that UPF1 specifically suppresses the translation of aberrant RNAs. Particularly striking was an increase in the translation of TIR domain-containing, nucleotide-binding, leucine-rich repeat (TNL) immune receptors. The regulation of TNLs via UPF1/NMD-mediated mRNA stability and translational de-repression offers a dynamic mechanism for the rapid activation of TNLs in response to pathogen attack.

Nonsense-mediated RNA decay (NMD) is an RNA control mechanism that has also been implicated in the broader regulation of gene expression. Nevertheless, a role for NMD in genome regulation has not yet been fully assessed, partially because NMD inactivation is lethal in many organisms. Here, we performed an in-depth comparative analysis of Arabidopsis (Arabidopsis thaliana) mutants lacking the NMD-related proteins UPF3, UPF1, and SMG7. We found different impacts of these proteins on NMD and the Arabidopsis transcriptome, with UPF1 having the biggest effect. Transcriptome assembly in UPF1-null plants revealed genome-wide changes in alternative splicing, suggesting that UPF1 functions in splicing. The inactivation of UPF1 led to translational repression, as manifested by a global shift in mRNAs from polysomes to monosomes and the downregulation of genes involved in translation and ribosome biogenesis. Despite these global changes, NMD targets and mRNAs expressed at low levels with short half-lives were enriched in the polysomes of upf1 mutants, indicating that UPF1/NMD suppresses the translation of aberrant RNAs. Particularly striking was an increase in the translation of TIR domain-containing, nucleotide binding, leucine-rich repeat (TNL) immune receptors. The regulation of TNLs via UPF1/NMD-mediated mRNA stability and translational derepression offers a dynamic mechanism for the rapid activation of TNLs in response to pathogen attack.

Background The changing climate is altering timing of key fruit ripening processes and increasing the occurrence of fruit defects. To improve our understanding of the genetic control of raspberry fruit development an enhanced genetic linkage map was developed and used to examine ripening phenotypic data. Results In this study we developed an enhanced genetic linkage map for the raspberry cvs. Glen Moy x Latham reference mapping population using genotyping by sequencing (GbS). Alignment to a newly sequenced draft reference genome of red raspberry, cultivar (cv.) Glen Moy, identified 8019 single nucleotide polymorphisms (SNPs). After stringent filtering to take account of read coverage over all the progeny individuals, association with a single chromosome, heterozygosity and marker regression mapping, 2348 high confidence SNPs were retained and integrated with an existing raspberry genetic map. The linkage map contained many more SNPs segregating in Latham than in Glen Moy. This caused difficulties in quantitative trait loci (QTL) mapping with standard software and a novel analysis based on a hidden Markov model was used to improve the mapping. QTL mapping using the newly generated dense genetic map not only corroborated previously identified genetic locations but also provided additional genetic elements controlling fruit ripening in raspberry. Conclusion The high-density GbS map located the QTL peaks more precisely than in earlier studies, aligned the QTLs with Glen Moy genome scaffolds, narrowed the range of potential candidate genes to these regions that can be utilised in other populations or in gene expression studies to confirm their role and increased the repertoire of markers available to understand the genetic control of fruit ripening traits. Electronic supplementary material The online version of this article (10.1186/s12863-018-0666-z) contains supplementary material, which is available to authorized users.

How alternative splicing (AS) is regulated in plants has not yet been elucidated. Previously, we have shown that the nuclear cap-binding protein complex (AtCBC) is involved in AS in Arabidopsis thaliana. Here we show that both subunits of AtCBC (AtCBP20 and AtCBP80) interact with SERRATE (AtSE), a protein involved in the microRNA biogenesis pathway. Moreover, using a high-resolution reverse transcriptase-polymerase chain reaction AS system we have found that AtSE influences AS in a similar way to the cap-binding complex (CBC), preferentially affecting selection of 5' splice site of first introns. The AtSE protein acts in cooperation with AtCBC: many changes observed in the mutant lacking the correct SERRATE activity were common to those observed in the cbp mutants. Interestingly, significant changes in AS of some genes were also observed in other mutants of plant microRNA biogenesis pathway, hyl1-2 and dcl1-7, but a majority of them did not correspond to the changes observed in the se-1 mutant. Thus, the role of SERRATE in AS regulation is distinct from that of HYL1 and DCL1, and is similar to the regulation of AS in which CBC is involved.

Assembly of the barley genome and extensive use of RNA-seq has resulted in an abundance of gene expression data and the recognition of wide-scale production of alternatively spliced transcripts. Here, we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) that confirms the accuracy of alternatively spliced transcripts from RNA-seq and allows quantification of changes in the proportion of splice isoforms between different experimental conditions, time points, tissues, genotypes, ecotypes, and treatments. By validating a selection of barley genes, use of the panel gives confidence or otherwise to the genome-wide global changes in alternatively spliced transcripts reported by RNA-seq. This simple assay can readily be applied to perform detailed transcript isoform analysis for any gene in any species.

Assembly of the barley genome and extensive use of RNA-seq has resulted in an abundance of gene expression data and the recognition of wide-scale production of alternatively spliced transcripts. Here, we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) that confirms the accuracy of alternatively spliced transcripts from RNA-seq and allows quantification of changes in the proportion of splice isoforms between different experimental conditions, time points, tissues, genotypes, ecotypes, and treatments. By validating a selection of barley genes, use of the panel gives confidence or otherwise to the genome-wide global changes in alternatively spliced transcripts reported by RNA-seq. This simple assay can readily be applied to perform detailed transcript isoform analysis for any gene in any species.

SummaryFor plants, light is the source of energy and the most relevant regulator of growth and adaptations to the environment by inducing changes in gene expression at various levels, including alternative splicing. Chloroplasts trigger retrograde signals that control alternative splicing in leaves and roots in response to light. Here we provide evidence suggesting that sugars, derived from photosynthesis, act as mobile signals controlling alternative splicing in roots. The inhibition of TOR kinase activity diminishes the alternative splicing response to light and/or sugars in roots, showing the relevance of the TOR pathway in this signaling mechanism. Furthermore, disrupting the function of the mitochondria abolishes alternative splicing changes supporting a key role for these organelles in this signaling axis. We conclude that sugars can act as mobile signals coordinating alternative splicing responses to light throughout the whole plant, exerting this function in roots by activating the TOR pathway.Graphical AbstractArt by Dr. Luciana Giono.

Climate change is impacting soft fruit crops. In raspberry, uneven bud break, greater variability in time to fruit ripening and crumbly fruit are already in evidence. Understanding the developmental process and how the environment impacts will be crucial in sustaining the industry in this changing climate against a background of biotic stresses. This chapter reviews regulation of processes leading to flowering time and fertilisation, developing fruit, ripening, colour, flavour and size. Recent developments of genomic and transcriptome tools which will have a significant role in breeding of the next generation of raspberry fruit are considered.

In this chapter, we refer to the expressed portion of the barley genome as the relatively small fraction of the total cellular DNA that either contains the genes that ultimately produce proteins, or that directly/indirectly controls the level, location and/or timing of when these genes are expressed and proteins are produced. We start by describing the dynamics of tissue and time-dependent gene expression and how common patterns across multiple samples can provide clues about gene networks involved in common biological processes. We then describe some of the complexities of how a single mRNA template can be differentially processed by alternative splicing to generate multiple different proteins or provide a mechanism to regulate the amount of functional gene product in a cell at a given point in time. We extend our analysis, using a number of biological examples, to address how diverse families of small non-coding microRNAs specifically regulate gene expression, and complete our appraisal by looking at the physical/molecular environment around genes that can result in either the promotion or repression of gene expression. We conclude by assessing some of the issues that remain around our ability to fully exploit the depth and power of current approaches for analysing gene expression and propose improvements that could be made using new but available sequencing and bioinformatics technologies.

Recent developments in genetics and genomics have advanced research in all crops including soft fruit species. Molecular markers which detect genome-wide variability in both protein coding and non-coding regions have enabled genetic mapping studies to move beyond linkages between simple morphological traits (Jennings 1967a, 1988; Ourecky 1975; Crane and Lawrence 1931; Keep 1968) to linkage maps containing numerous genetic markers which can be utilised in marker assisted breeding. Until recently, mapping in blackberry and other Rubus species has lagged behind that of red raspberry due to their more complex genetic make-up and lesser economic importance.

Summary This paper examines the function of Arabidopsis thaliana AtPTB1 and AtPTB2 as plant splicing factors. The effect on splicing of overexpression of AtPTB1 and AtPTB2 was analysed in an in vivo protoplast transient expression system with a novel mini-exon splicing reporter. A range of mutations in pyrimidine-rich sequences were compared with and without AtPTB and NpU2AF65 overexpression. Splicing analyses of constructs in protoplasts and RNA from overexpression lines used high-resolution reverse transcription polymerase chain reaction (RT-PCR). AtPTB1 and AtPTB2 reduced inclusion/splicing of the potato invertase mini-exon splicing reporter, indicating that these proteins can repress plant intron splicing. Mutation of the polypyrimidine tract and closely associated Cytosine and Uracil-rich (CU-rich) sequences, upstream of the mini-exon, altered repression by AtPTB1 and AtPTB2. Coexpression of a plant orthologue of U2AF65 alleviated the splicing repression of AtPTB1. Mutation of a second CU-rich upstream of the mini-exon 3′ splice site led to a decline in mini-exon splicing, indicating the presence of a splicing enhancer sequence. Finally, RT-PCR of AtPTB overexpression lines with c. 90 known alternative splicing (AS) events showed that AtPTBs significantly altered AS of over half the events. AtPTB1 and AtPTB2 are splicing factors that influence alternative splicing. This occurs in the potato invertase mini-exon via the polypyrimidine tract and associated pyrimidine-rich sequence.

Light makes carbon fixation possible, allowing plant and animal life on Earth. We have previously shown that light regulates alternative splicing in plants. Light initiates a chloroplast retrograde signaling that regulates nuclear alternative splicing of a subset of Arabidopsis thaliana transcripts. Here, we show that light promotes RNA polymerase II (Pol II) elongation in the affected genes, whereas in darkness, elongation is lower. These changes in transcription are consistent with elongation causing the observed changes in alternative splicing, as revealed by different drug treatments and genetic evidence. The light control of splicing and elongation is abolished in an Arabidopsis mutant defective in the transcription factor IIS (TFIIS). We report that the chloroplast control of nuclear alternative splicing in plants responds to the kinetic coupling mechanism found in mammalian cells, providing unique evidence that coupling is important for a whole organism to respond to environmental cues.

Although significant work has been undertaken regarding the response of model and crop plants to heat shock during the acclimatory phase, few studies have examined the steady-state response to the mild heat stress encountered in temperate agriculture. In the present work, we therefore exposed tuberizing potato plants to mildly elevated temperatures (30/20 °C, day/night) for up to 5 weeks and compared tuber yield, physiological and biochemical responses, and leaf and tuber metabolomes and transcriptomes with plants grown under optimal conditions (22/16 °C). Growth at elevated temperature reduced tuber yield despite an increase in net foliar photosynthesis. This was associated with major shifts in leaf and tuber metabolite profiles, a significant decrease in leaf glutathione redox state and decreased starch synthesis in tubers. Furthermore, growth at elevated temperature had a profound impact on leaf and tuber transcript expression with large numbers of transcripts displaying a rhythmic oscillation at the higher growth temperature. RT-PCR revealed perturbation in the expression of circadian clock transcripts including StSP6A, previously identified as a tuberization signal. Our data indicate that potato plants grown at moderately elevated temperatures do not exhibit classic symptoms of abiotic stress but that tuber development responds via a diversity of biochemical and molecular signals.

Alternative splicing (AS) is a regulated mechanism that generates multiple transcripts from individual genes. It is widespread in eukaryotic genomes and provides an effective way to control gene expression. At low temperatures, AS regulates Arabidopsis clock genes through dynamic changes in the levels of productive mRNAs. We examined AS in barley clock genes to assess whether temperature-dependent AS responses also occur in a monocotyledonous crop species. We identify changes in AS of various barley core clock genes including the barley orthologues of Arabidopsis AtLHY and AtPRR7 which showed the most pronounced AS changes in response to low temperature. The AS events modulate the levels of functional and translatable mRNAs, and potentially protein levels, upon transition to cold. There is some conservation of AS events and/or splicing behaviour of clock genes between Arabidopsis and barley. In addition, novel temperature-dependent AS of the core clock gene HvPPD-H1 (a major determinant of photoperiod response and AtPRR7 orthologue) is conserved in monocots. HvPPD-H1 showed a rapid, temperature-sensitive isoform switch which resulted in changes in abundance of AS variants encoding different protein isoforms. This novel layer of low temperature control of clock gene expression, observed in two very different species, will help our understanding of plant adaptation to different environments and ultimately offer a new range of targets for plant improvement.

QTL mapping identifies a range of underlying and unrelated genes with apparent roles in raspberry fruit ripening and softening that show characteristic developing fruit expression profiles. Fruit softening is an important agronomical trait that involves a complex interaction of plant cell processes. We have used both qualitative and quantitative scoring of fruit firmness, length, mass, and resistance to applied force to identify QTL in a raspberry mapping population. QTLs were located primarily on linkage group (LG) 3 with other significant loci on LG 1 and LG 5 and showed mostly additive effects between the two parents. The expression of key genes that underlie these QTLs with roles in cell-wall solubility, water uptake, polyamine synthesis, transcription, and cell respiration was tested across five stages of fruit development, from immature green to red ripe fruit, using real-time RT-qPCR. Gene expression patterns showed variable expression patterns across fruit development with a highly significant positive and negative correlation between genes, supporting precise regulation of expression of different cell processes throughout raspberry fruit development. Variable timing in expression was also found in some genes at different fruit development stages between soft and firm cultivars. Multiple processes have a role to play in fruit softening and this will require development of multiple marker combinations to genes that characterise raspberry fruit softening.

Posttranscriptional control makes an important contribution to circadian regulation of gene expression. In higher plants, alternative splicing is particularly prevalent upon abiotic and biotic stress and in the circadian system. Here we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) to monitor alternative splicing events. The use of the panel allows the quantification of changes in the proportion of splice isoforms between different samples, e.g., different time points, different tissues, genotypes, ecotypes, or treatments.

Alternative splicing (AS) of pre-mRNAs is an important regulatory mechanism shaping the transcriptome. In plants, only few RNA-binding proteins are known to affect AS. Here, we show that the glycine-rich RNA-binding protein AtGRP7 influences AS in Arabidopsis thaliana. Using a high-resolution RT-PCR-based AS panel, we found significant changes in the ratios of AS isoforms for 59 of 288 analyzed AS events upon ectopic AtGRP7 expression. In particular, AtGRP7 affected the choice of alternative 5' splice sites preferentially. About half of the events are also influenced by the paralog AtGRP8, indicating that AtGRP7 and AtGRP8 share a network of downstream targets. For 10 events, the AS patterns were altered in opposite directions in plants with elevated AtGRP7 level or lacking AtGRP7. Importantly, RNA immunoprecipitation from plant extracts showed that several transcripts are bound by AtGRP7 in vivo and indeed represent direct targets. Furthermore, the effect of AtGRP7 on these AS events was abrogated by mutation of a single arginine that is required for its RNA-binding activity. This indicates that AtGRP7 impacts AS of these transcripts via direct interaction. As several of the AS events are also controlled by other splicing regulators, our data begin to provide insights into an AS network in Arabidopsis.

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT–PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5′ or 3′UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3′-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

AS (alternative splicing) is a post-transcriptional process which regulates gene expression through increasing protein complexity and modulating mRNA transcript levels. Regulation of AS depends on interactions between trans-acting protein factors and cis-acting signals in the pre-mRNA (precursor mRNA) transcripts, termed ‘combinatorial’ control. Dynamic changes in AS patterns reflect changes in abundance, composition and activity of splicing factors in different cell types and in response to cellular or environmental cues. Whereas the SR protein family of splicing factors is well-studied in plants, relatively little is known about other factors influencing the regulation of AS or the consequences of AS on mRNA levels and protein function. To address fundamental questions on AS in plants, we are exploiting a high-resolution RT (reverse transcription)–PCR system to analyse multiple AS events simultaneously. In the present paper, we describe the current applications and development of the AS RT–PCR panel in investigating the roles of splicing factors, cap-binding proteins and nonsense-mediated decay proteins on AS, and examining the extent of AS in genes involved in the same developmental pathway or process.

Here, we identify the Arabidopsis thaliana ortholog of the mammalian DEAD box helicase, eIF4A-III, the putative anchor protein of exon junction complex (EJC) on mRNA. Arabidopsis eIF4A-III interacts with an ortholog of the core EJC component, ALY/Ref, and colocalizes with other EJC components, such as Mago, Y14, and RNPS1, suggesting a similar function in EJC assembly to animal eIF4A-III. A green fluorescent protein (GFP)-eIF4A-III fusion protein showed localization to several subnuclear domains: to the nucleoplasm during normal growth and to the nucleolus and splicing speckles in response to hypoxia. Treatment with the respiratory inhibitor sodium azide produced an identical response to the hypoxia stress. Treatment with the proteasome inhibitor MG132 led to accumulation of GFP-eIF4A-III mainly in the nucleolus, suggesting that transition of eIF4A-III between subnuclear domains and/or accumulation in nuclear speckles is controlled by proteolysis-labile factors. As revealed by fluorescence recovery after photobleaching analysis, the nucleoplasmic fraction was highly mobile, while the speckles were the least mobile fractions, and the nucleolar fraction had an intermediate mobility. Sequestration of eIF4A-III into nuclear pools with different mobility is likely to reflect the transcriptional and mRNA processing state of the cell.

Alternative splicing (AS) increases the proteomic and functional capacity of genomes through the generation of alternative mRNA transcripts from the same gene. AS is now estimated to occur in a third of Arabidopsis and rice genes, and includes genes involved in the control of growth and development, responses to stress and signalling. Regulation of AS reflects the interactions between positive and negative cis sequences in the precursor messenger RNA and a range of trans-acting factors. The levels and activities of these factors differ in different cells and growth conditions. To identify changes in AS in multiple genes simultaneously, we have established a reproducible RT-PCR panel that can analyse 96 alternative splicing events and accurately measure the ratio of alternatively spliced products. This procedure detected statistically significant changes in AS in different plant organs, in plants grown under different light and day-length conditions, and in plants overexpressing splicing factors. The system provides a convenient, medium-throughput means of monitoring changes in AS in multiple genes. It can readily be applied to much larger or targeted sets of gene transcripts to generate information on the significance and regulation of AS in plant growth and development, specific processes and responses to external stimuli.

Transcript annotation in plant databases is incomplete and often inaccurate, leading to misinterpretation. As more and more RNA-seq data are generated, plant scientists need to be aware of potential pitfalls and understand the nature and impact of specific alternative splicing transcripts on protein production. A primary area of concern and the topic of this article is the (mis)annotation of open reading frames and premature termination codons. The basic message is that to adequately address expression and functions of transcript isoforms, it is necessary to be able to predict their fate in terms of whether protein isoforms are generated or specific transcripts are unproductive or degraded.

RNA-sequencing (RNA-seq) allows global gene expression analysis at the individual transcript level. Accurate quantification of transcript variants generated by alternative splicing (AS) remains a challenge. We have developed a comprehensive, nonredundant Arabidopsis reference transcript dataset (AtRTD) containing over 74 000 transcripts for use with algorithms to quantify AS transcript isoforms in RNA-seq. The AtRTD was formed by merging transcripts from TAIR10 and novel transcripts identified in an AS discovery project. We have estimated transcript abundance in RNA-seq data using the transcriptome-based alignment-free programmes Sailfish and Salmon and have validated quantification of splicing ratios from RNA-seq by high resolution reverse transcription polymerase chain reaction (HR RT-PCR). Good correlations between splicing ratios from RNA-seq and HR RT-PCR were obtained demonstrating the accuracy of abundances calculated for individual transcripts in RNA-seq. The AtRTD is a resource that will have immediate utility in analysing Arabidopsis RNA-seq data to quantify differential transcript abundance and expression. This research was supported by funding from the Biotechnology and Biological Sciences Research Council (BBSRC) (BB/K006568/1 to J.W.S.B.; BB/K006835/1 to H.G.N.), the Scottish Government Rural and Environment Science and Analytical Services division (RESAS) and by the Austrian Science Fund (FWF) (P26333) to M.K. and (DK W1207) to A.B. The authors acknowledge the European Alternative Splicing Network of Excellence (EURASNET), LSHG-CT-2005-518238 for catalysing important collaborations.

We describe a polymerase chain reaction (PCR) protocol suitable for use in secondary schools and colleges. This PCR protocol can be used to investigate genetic variation between plants. The protocol makes use of primers which are complementary to sequences of nucleotides that are highly conserved across different plant genera. The regions of chloroplast DNA amplified lie between these conserved sequences and are non-coding. These non-coding regions display a high frequency of mutations and show relatively high rates of evolutionary change. Thus it is possible to use the protocol to explore evolutionary relationships between plants. Results from Brassica oleracea can be used to highlight genetic similarity and differences within and across genera. The protocol is robust and is suitable for use either with a thermocycler or a series of water-baths, thus making it accessible for use in most schools and colleges.

In potato invertase genes, the constitutively included, 9-nucleotide (nt)-long mini-exon requires a strong branchpoint and U-rich polypyrimidine tract for inclusion. The strength of these splicing signals was demonstrated by greatly enhanced splicing of a poorly spliced intron and by their ability to support splicing of an artificial mini-exon, following their introduction. Plant introns also require a second splicing signal, UA-rich intronic elements, for efficient intron splicing. Mutation of the branchpoint caused loss of mini-exon inclusion without loss of splicing enhancement, showing that the same U-rich sequence can function as either a polypyrimidine tract or a UA-rich intronic element. The distinction between the splicing signals depended on intron context (the presence or absence of an upstream, adjacent and functional branchpoint), and on the sequence context of the U-rich elements. Polypyrimidine tracts tolerated C residues while UA-rich intronic elements tolerated As. Thus, in plant introns, U-rich splicing elements can have dual roles as either a general plant U-rich splicing signal or a polypyrimidine tract. Finally, overexpression of two different U-rich binding proteins enhanced intron recognition significantly. These results highlight the importance of co-operation between splicing signals, the importance of other nucleotides within U-rich elements for optimal binding of competing splicing factors and effects on splicing efficiency of U-rich binding proteins.

Constitutive splicing of the potato invertase miniexon 2 (9 nt long) requires a branchpoint sequence positioned around 50 nt upstream of the 5′ splice site of the adjacent intron and a U11 element found just downstream of the branchpoint in the upstream intron [Simpson, Hedley, Watters, Clark, McQuade, Machray and Brown (2000) RNA 6, 422–433]. The sensitivity of this in vivo plant splicing system has been used to demonstrate exon scanning in plants, and to characterize plant intronic elements, such as branchpoint and poly-pyrimidine tract sequences. Plant introns differ from their vertebrate and yeast couterparts in being UA- or U-rich (up to 85% UA). One of the key differences in splicing between plants and other eukaryotes lies in early intron recognition, which is thought to be mediated by UA-binding proteins. We are adopting three approaches to studying the RNA-protein interactions in plant splicing. First, overexpression of plant splicing factors and, in particular, UA-binding proteins, in conjunction with a range of mini-exon mutants. Secondly, the sequences of around 65% of vertebrate and yeast splicing factors have high-quality matches to Arabidopsis proteins, opening the door to identification and analysis of gene knockouts. Finally, to discover plant-specific proteins involved in splicing and in, for example, rRNA or small nuclear RNA processing, green fluorescent protein-cDNA fusion libraries in viral vectors are being screened.

The branchpoint sequence and associated polypyrimidine tract are firmly established splicing signals in vertebrates. In plants, however, these signals have not been characterized in detail. The potato invertase mini-exon 2 (9 nt) requires a branchpoint sequence positioned around 50 nt upstream of the 5' splice site of the neighboring intron and a U11 element found adjacent to the branchpoint in the upstream intron (Simpson et al., RNA, 2000, 6:422-433). Utilizing the sensitivity of this plant splicing system, these elements have been characterized by systematic mutation and analysis of the effect on inclusion of the mini-exon. Mutation of the branchpoint sequence in all possible positions demonstrated that branchpoints matching the consensus, CURAY, were most efficient at supporting splicing. Branchpoint sequences that differed from this consensus were still able to permit mini-exon inclusion but at greatly reduced levels. Mutation of the downstream U11 element suggested that it functioned as a polypyrimidine tract rather than a UA-rich element, common to plant introns. The minimum sequence requirement of the polypyrimidine tract for efficient splicing was two closely positioned groups of uridines 3-4 nt long (

International audience; Alternative splicing (AS) of pre-mRNA represents a major mechanism underlying increased transcriptome and proteome complexity. Here, we show that the nuclear speckle RNA-binding protein (NSR) and the AS competitor long noncoding RNA (or ASCO-IncRNA) constitute an AS regulatory module. AtNSR-GFP translational fusions are expressed in primary and lateral root (LR) meristems. Double Atnsr mutants and ASCO overexpressors exhibit an altered ability to form LRs after auxin treatment. Interestingly, auxin induces a major change in AS patterns of many genes, a response largely dependent on NSRs. RNA immunoprecipitation assays demonstrate that AtNSRs interact not only with their alternatively spliced mRNA targets but also with the ASCO-RNA in vivo. The ASCO-RNA displaces an AS target from an NSR-containing complex in vitro. Expression of ASCO-RNA in Arabidopsis affects the splicing patterns of several NSR-regulated mRNA targets. Hence, IncRNA can hijack nuclear AS regulators to modulate AS patterns during development.

Summary Exon definition is a mechanism whereby splice sites are selected initially via interactions between splicing factors across an exon, prior to spliceosome assembly and intron removal. It occurs in the splicing of vertebrate pre-mRNAs and, recently, evidence for exon definition and the role of exon sequences has been obtained in plant intron splicing. Here we demonstrate that interactions between plant introns influence splicing efficiency and that these interactions are consistent with an exon definition process. The splicing efficiency of a UA-poor, inefficiently spliced intron (wheat amylase) increases 3.5- to 4.4-fold when placed in tandem with a UA-rich, well spliced, intron (legumin). Enhanced splicing is also observed with partial pea legumin intron sequences. However, mutation of splice sites in the partial UA-rich intron sequences abolished the enhanced splicing effect such that intact splice sites at the 5′ and 3′ ends of the exon were required, thus pointing to exon definition. This was further supported by reducing the size of the intervening exon or replacing with a UA-rich sequence which leads to loss of splicing of the UA-poor intron. Finally, the results support UA-rich sequences functioning early in the splicing process in plants.

The impact of AS (alternative splicing) is well-recognized in animal systems as a key regulator of gene expression and proteome complexity. In plants, AS is of growing importance as more genes are found to undergo AS, but relatively little is known about the factors regulating AS or the consequences of AS on mRNA levels and protein function. We have established an accurate and reproducible RT (reverse transcription)–PCR system to analyse AS in multiple genes. Initial studies have identified new AS events confirming that current values for the frequency of AS in plants are likely to be underestimates.

The recent isolation of a number of plant box C/D small nucleolar (sno)RNAs demonstrates theconservation in plants of sequence and structural elements of processed box C/D snoRNAs. Boxes C andD, and terminal inverted repeats are known to be essential for accumulation and processing in vertebratesand yeast. Processing of vertebrate box C/D snoRNAs was examined by expression of various mousehsc70 intron 5-U14 constructs in tobacco protoplasts. Full-length U14 and internally deleted U14 accumu-lated in the plant cells. Human U3 and U8 fragments, consistent with processing to internal box C/C′sequences, also accumulated in the plant cells. The similarity of processing behaviour of the vertebratebox C/D constructs in tobacco protoplasts and Xenopus oocytes suggests the mechanism of processing,involving recognition and association of proteins, is conserved in plants.Keywords: small nucleolar RNA; box C/D small nucleolar RNA; small nucleolar RNA processing; U14.Eukaryotic nuclei contain two major classes of small nucleo- dependent manner (Kiss and Filipowicz, 1995; Cavaille´andlar RNAs (snoRNAs): box C/D and box H/ACA snoRNAs (re- Bachellerie,1996; Kiss et al., 1996). Processing is likely to in-viewed in Maxwell and Fournier, 1995; Tollervey and Kiss, volve exonucleolytic digestion of 5′and 3′sequences flanking1997; Smith and Steitz, 1997). Some box C/D snoRNAs, for the snoRNA, following debranching of the intron lariat (Kissexample, U3, U8, U13, U14 and U22 are involved in processing and Filipowicz, 1995; Cavaille´and Bachellerie, 1996; Kiss etsteps required in the production of mature 18S, 5.8S and 28S al., 1996). In Xenopus oocytes, where processing of U16 andribosomal RNA (rRNA; Maxwell and Fournier, 1995; Tollervey U18 occurs by initial endonucleolytic cleavage in the flankingand Kiss, 1997; Peculis and Steitz, 1994; Tycowski et al., 1994; intron sequences, processing to the mature snoRNA is againBeltrame and Tollervey, 1995; Liang and Fournier, 1995; Ca- likely to be exonucleolytic (Caffarelli et al., 1996).vaille´et al., 1996a). Other box C/D snoRNAs, including U14, Maturation of intronic box C/D and box H/ACA small nucle-contain sequences which are complementary to 18S and 28S olar ribonucleoprotein particles (snoRNPs) requires blockage ofrRNAs and act as guide RNAs to specify sites of 2′-O-ribose exonucleolytic digestion by snoRNP formation via binding ofmethylation of rRNAs (Bachellerie et al., 1995; Cavaille´et al., specific proteins to the snoRNA. The intronic box C/D snoRNAs1996b; Kiss-La´szlo´et al., 1996; Nicoloso et al., 1996; Tykow- contain the conserved box C and D sequences near their 5′andski et al., 1996a, b). The box H/ACA snoRNAs, E3 and U17, 3′termini, respectively, usually flanked by terminal inverted re-are involved in the 5′external transcribed spacer cleavage in peats. The formation of a stem between the inverted repeats and/precursor rRNA (pre-rRNA; Enright et al., 1996). E1 and E2 or the secondary structure of the snoRNA bring the boxes C andare needed for 18S rRNA production and E3 appears to be in- D into juxtaposition, and binding of proteins to the snoRNA andvolved in cleavage at or near the 5′end of 5.8S rRNA (Mishra more particularly to the stem-box C/D structure serves to blockand Eliceiri, 1997) and in yeast snR10 and snR30 are required exonucleolytic degradation (Tycowski et al., 1993; Cafarrelli etfor normal rRNA processing and cell growth (Tollervey, 1987; al., 1996; Watkins et al., 1996; Xia et al., 1997). These se-Morrisey and Tollervey, 1993). However, the majority of box quences/structure have been shown to be sufficient for efficientH/ACA snoRNAs are predicted to be involved in defining sites processing (Caffarelli et al., 1996; Cavaille´ and Bachellerie,of pseudouridylation in rRNAs again involving base-pairing in- 1996; Watkins et al., 1996; Xia et al., 1997). In addition, pro-teractions (Bousquet-Antonelli et al., 1997; Ganot et al., teins including fibrillarin have been shown to associate with the1997a, b; Ni et al., 1997). snoRNA prior to processing (Caffarelli et al., 1996; Watkins etAlthough U3, U8 and U13 are expressed from their own al., 1996). Box H/ACA snoRNAs are also protected from exo-promoters (Dahlberg and Lund, 1988; Tyc and Steitz, 1989), the nucleolytic digestion by virtue of stem-loop structures in the 5′majority of vertebrate box C/D and box H/ACA snoRNAs are and 3′halves of the molecules and by the association of proteins,encoded in introns of protein coding genes (Maxwell and Four- such as Gar1p. The protein components of box H/ACA snoRNPsnier, 1995; Tollervey and Kiss, 1997). The intronic snoRNAs probably bind the conserved box H and ACA sequences whichare produced by processing of the intron in a largely splicing- are essential for accumulation and processing (Kiss and Filipo-wicz, 1995; Kiss et al., 1996; Balakin et al., 1996; Bousquet

Plant polygalacturonase-inhibiting proteins (PGIPs) inhibit endo-polygalacturonases (endo-PGs) released by invasive fungi. Two potential PGIP cDNAs have been cloned from raspberry (Rubus idaeus). PGIP1 corresponds to a full-length PGIP cDNA with a high degree of identity with previously isolated genes and maintains all the characteristic features of PGIP peptides. A genomic fragment of PGIP1 identifies a single 243 bp intron which is efficiently spliced out of the PGIP1 pre-mRNA transcript. PGIP2 contains a frame shift mutation which would lead to the loss of 105 amino acids off the C-terminus. Contrary to PGIP activity levels in raspberry (Johnston et al., 1993), expression analysis detected PGIP1 at equal levels throughout fruit development from closed flower, through to ripe fruit.

Intron sequences from monocotyledonous and dicotyledonous origin were used to abolish marker gene expression in prokaryotes (Escherichia coli and Agrobacterium tumefaciens) but permit expression in selected eukaryotic systems using the eukaryotic specific splicing mechanism. A 1014 bp maize Shrunken-1 (Sh 1) intron 1 flanked by exon1 and exon2 sequences was cloned into the N-terminal of the NPT II-coding region. Transient gene expression analysis revealed that the modified neomycin phosphotransferase II (NPT II) gene, driven by the cauliflower mosaic virus (CaMV) 35S promoter, is expressed in barley protoplasts, but poorly expressed in tobacco protoplasts. In dicotyledonous cells AU-rich sequences are known to be important for efficient splicing and therefore an attempt was made to improve expression of the NPT II gene, containing the Sh 1 intron 1, in tobacco by increasing the AU content from 57% to 69%. Reverse transcriptase PCR analysis of RNA from transiently expressed NPT II transcripts from tobacco protoplasts revealed that despite the increase in AU-content, NPT II was still poorly expressed. Cryptic splice sites were identified as one possible cause for missplicing of the Sh1 intron 1 in dicots and poor levels of expression. Alternatively, cloning of the 198 bp intron 2 of the potato STLS 1 gene (81% AU) into the N-terminal part of the NPT II-coding region resulted in proper expression of NPT II in tobacco as well as in barley protoplasts and abolished marker gene expression in prokaryotes. The successful insertion of an intron into a selectable marker gene which completely abolishes gene expression in prokaryotes, without affecting expression of chimeric genes in monocotyledonous and dicotyledonous plant cells provides a suitable system to reduce the number of false-positives in transgenic plant production.

Although significant work has been undertaken regarding the response of model and crop plants to heat shock during the acclimatory phase, few studies have examined the steady-state response to the mild heat stress encountered in temperate agriculture. In the present work, we therefore exposed tuberizing potato plants to mildly elevated temperatures (30/20 °C, day/night) for up to 5 weeks and compared tuber yield, physiological and biochemical responses, and leaf and tuber metabolomes and transcriptomes with plants grown under optimal conditions (22/16 °C). Growth at elevated temperature reduced tuber yield despite an increase in net foliar photosynthesis. This was associated with major shifts in leaf and tuber metabolite profiles, a significant decrease in leaf glutathione redox state and decreased starch synthesis in tubers. Furthermore, growth at elevated temperature had a profound impact on leaf and tuber transcript expression with large numbers of transcripts displaying a rhythmic oscillation at the higher growth temperature. RT-PCR revealed perturbation in the expression of circadian clock transcripts including StSP6A, previously identified as a tuberization signal. Our data indicate that potato plants grown at moderately elevated temperatures do not exhibit classic symptoms of abiotic stress but that tuber development responds via a diversity of biochemical and molecular signals.

Intron lariat formation between the 5' end of an intron and a branchpoint adenosine is a fundamental aspect of the first step in animal and yeast nuclear pre-mRNA splicing. Despite similarities in intron sequence requirements and the components of splicing, differences exist between the splicing of plant and vertebrate introns. The identification of AU-rich sequences as major functional elements in plant introns and the demonstration that a branchpoint consensus sequence was not required for splicing have led to the suggestion that the transition from AU-rich intron to GC-rich exon is a major potential signal by which plant pre-mRNA splice sites are recognized. The role of putative branchpoint sequences as an internal signal in plant intron recognition/definition has been re-examined. Single nucleotide mutations in putative branchpoint adenosines contained within CUNAN sequences in four different plant introns all significantly reduced splicing efficiency. These results provide the most direct evidence to date for preferred branchpoint sequences being required for the efficient splicing of at least some plant introns in addition to the important role played by AU sequences in dicot intron recognition. The observed patterns of 3' splice site selection in the introns studied are consistent with the scanning model described for animal intron 3' splice site selection. It is suggested that, despite the clear importance of AU sequences for plant intron splicing, the fundamental processes of splice site selection and splicing in plants are similar to those in animals.

The circadian clock plays a crucial role in coordinating plant metabolic and physiological functions with predictable environmental variables, such as dusk and dawn, while also modulating responses to biotic and abiotic challenges. Much of the initial characterization of the circadian system has focused on transcriptional initiation, but it is now apparent that considerable regulation is exerted after this key regulatory step. Transcript processing, protein stability, and cofactor availability have all been reported to influence circadian rhythms in a variety of species. We used a genetic screen to identify a mutation within a putative RNA binding protein (SPLICEOSOMAL TIMEKEEPER LOCUS1 [STIPL1]) that induces a long circadian period phenotype under constant conditions. STIPL1 is a homolog of the spliceosomal proteins TFP11 (Homo sapiens) and Ntr1p (Saccharomyces cerevisiae) involved in spliceosome disassembly. Analysis of general and alternative splicing using a high-resolution RT-PCR system revealed that mutation of this protein causes less efficient splicing of most but not all of the introns analyzed. In particular, the altered accumulation of circadian-associated transcripts may contribute to the observed mutant phenotype. Interestingly, mutation of a close homolog of STIPL1, STIP-LIKE2, does not cause a circadian phenotype, which suggests divergence in function between these family members. Our work highlights the importance of posttranscriptional control within the clock mechanism.

A polymerase chain reaction (PCR) strategy designed to amplify DNA sequences between closely linked U2snRNA genes has generated extensive coding and 5' regulatory sequence information on the potato U2snRNA multigene family. Two of the U2snRNA coding sequences isolated differed substantially from normal U2snRNAs by containing both complementary deletions and regions of novel sequence. However, sequences such as Sm-binding sites and loops of stem-loops III and IV, which are some of the most highly conserved regions in U2snRNA, remain highly conserved in these genes. The complementary deletions would effectively remove stem-loop IIb which has been shown in yeast to be unnecessary for pre-mRNA splicing. Transcripts from one of the genes have been detected by reverse transcriptase-PCR (RT-PCR) in total RNA. These novel U2snRNA genes represent the first reported example of naturally occurring structural variants and provide support for the proposed non-essential role of U2snRNA stem-loop IIb.