High-quality reference transcript datasets hold the key to transcript-specific RNA-sequencing analysis in plants

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

Comparison and improvement of RNA extraction methods in Sargassum (Phaeophyta)

Sargassum (Sargassaceae) is widely distributed globally and plays an important role in regulating climate change, but the landscape of genomes and transcripts is less known. High-quality nucleic acids are the basis for molecular biology experiments such as high-throughput sequencing. Although extensive studies have documented methods of RNA extraction, these methods are not very applicable to Sargassum, which contains high levels of polysaccharides and polyphenols. To find a suitable method to improve the quality of RNA extracted, we compared and modified several popular RNA extraction methods and screened one practical method with three specific Sargassum spp. The results showed that three CTAB methods (denoted as Methods 1, 2, and 3) and the RNAprep Pure Plant Kit (denoted as Method 4) could, with slight modifications, effectively isolate RNA from Sargassum species, except for Method 4 used with S. fusiforme. By performing further screening, we determined Method 4 was the best choice for S. hemiphyllum and S. henslowianum, as revealed by RNA yields, RNA Integrity Number (RIN), extraction time, and unigene mapped ratio. For S. fusiforme, Methods 1, 2, and 3 showed no obvious differences among the yields, quality, or time to perform. In addition, one other method was tested, but we found the quality of the RNA extracted by TRIzol reagent methods (denoted as Method 5) performed the worst when compared with the above four methods. Therefore, our study provides four suitable methods for RNA extraction in Sargassum and is essential for future genetic exploration of Sargassum.

High-quality reference transcriptome construction improves RNA-seq quantification in Oryza sativa indica

The Reference Transcriptomic Dataset (RTD) is an accurate and comprehensive collection of transcripts originating from a given organism. It holds the key to precise transcript quantification and downstream analysis of differential expressions and regulations. Currently, transcriptome annotations for most crop plants are far from complete. For example, Oryza sativa indica (O. sativa indica) is reported to have 40,759 transcripts in the Ensembl database without alternative transcript isoforms and alternative splicing (AS) events. To generate a high-quality RTD, we conducted RNA sequencing of rice leaf samples collected at various time points during Rhizoctonia solani infection. The obtained reads were analyzed by adopting the recently developed computational analysis pipeline to assemble the RTD with increased transcript and AS diversity for O. sativa indica (IndicaRTD). After stringent quality filtering, the newly constructed transcriptome annotation was comprised of 122,968 non-redundant transcripts from 53,695 genes. This study identified many novel transcripts compared to Ensembl deposited data that are important for regulating molecular and physiological processes in the plant system. Currently, the assembled IndicaRTD must allow fast quantification of transcript and gene expression with high precision.

Differential tolerance to heat stress of young leaves compared to mature leaves of whole plants relate to differential transcriptomes involved in metabolic adaptations to stress

Plants respond to heat shock by regulating gene expression. While transcriptomic changes in response to heat stress are well studied, it is not known whether young and old leaves reprogram transcription differently upon stress. When whole plants of Arabidopsis thaliana were subjected to heat shock, young leaves were affected significantly less than older leaves based on measurements of tissue damage. To identify quantitative changes to transcriptomes between young and old leaves upon heat stress, we used RNA sequencing on young and old leaves from plants exposed to control and heat stress at 42 °C for 1 h and 10 h. A total of 6472 differentially expressed genes between young and old leaf were identified under control condition, and 9126 and 6891 under 1 h and 10 h heat stress, respectively. Analyses of differentially expressed transcripts led to the identification of multiple functional clusters of genes that may have potential roles in the increased heat tolerance of young leaves including higher level of expression in young leaves of genes encoding chaperones, heat shock proteins and proteins known in oxidative stress resistance. Differential levels of transcripts for genes implicated in pectin metabolism, cutin and wax biosynthesis, pentose and glucuronate interconversions, cellulose degradation, indole glucosinolate metabolism and RNA splicing between young and old leaves under heat stress suggest that cell wall remodelling, cuticular wax synthesis and carbohydrate modifications impacted by alternative splicing may also have roles in the improved heat stress tolerance of young leaves.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables 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 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 SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction 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.

The value of genotype-specific reference for transcriptome analyses in barley

It is increasingly apparent that although different genotypes within a species share "core" genes, they also contain variable numbers of "specific" genes and different structures of "core" genes that are only present in a subset of individuals. Using a common reference genome may thus lead to a loss of genotype-specific information in the assembled Reference Transcript Dataset (RTD) and the generation of erroneous, incomplete or misleading transcriptomics analysis results. In this study, we assembled genotype-specific RTD (sRTD) and common reference-based RTD (cRTD) from RNA-seq data of cultivated Barke and Morex barley, respectively. Our quantitative evaluation showed that the sRTD has a significantly higher diversity of transcripts and alternative splicing events, whereas the cRTD missed 40% of transcripts present in the sRTD and it only has ∼70% accurate transcript assemblies. We found that the sRTD is more accurate for transcript quantification as well as differential expression analysis. However, gene-level quantification is less affected, which may be a reasonable compromise when a high-quality genotype-specific reference is not available.

Transcriptome Analyses of Barley Roots Inoculated with Novel Paenibacillus sp. and Erwinia gerundensis Strains Reveal Beneficial Early-Stage Plant-Bacteria Interactions

Plant growth-promoting bacteria can improve host plant traits including nutrient uptake and metabolism and tolerance to biotic and abiotic stresses. Understanding the molecular basis of plant–bacteria interactions using dual RNA-seq analyses provides key knowledge of both host and bacteria simultaneously, leading to future enhancements of beneficial interactions. In this study, dual RNA-seq analyses were performed to provide insights into the early-stage interactions between barley seedlings and three novel bacterial strains (two Paenibacillus sp. strains and one Erwinia gerundensis strain) isolated from the perennial ryegrass seed microbiome. Differentially expressed bacterial and barley genes/transcripts involved in plant–bacteria interactions were identified, with varying species- and strain-specific responses. Overall, transcriptome profiles suggested that all three strains improved stress response, signal transduction, and nutrient uptake and metabolism of barley seedlings. Results also suggested potential improvements in seedling root growth via repressing ethylene biosynthesis in roots. Bacterial secondary metabolite gene clusters producing compounds that are potentially associated with interactions with the barley endophytic microbiome and associated with stress tolerance of plants under nutrient limiting conditions were also identified. The results of this study provided the molecular basis of plant growth-promoting activities of three novel bacterial strains in barley, laid a solid foundation for the future development of these three bacterial strains as biofertilisers, and identified key differences between bacterial strains of the same species in their responses to plants.

Reading between the Lines: RNA-seq Data Mining Reveals the Alternative Message of the Rice Leaf Transcriptome in Response to Heat Stress

Rice (Oryza sativa L.) is a major food crop but heat stress affects its yield and grain quality. To identify mechanistic solutions to improve rice yield under rising temperatures, molecular responses of thermotolerance must be understood. Transcriptional and post-transcriptional controls are involved in a wide range of plant environmental responses. Alternative splicing (AS), in particular, is a widespread mechanism impacting the stress defence in plants but it has been completely overlooked in rice genome-wide heat stress studies. In this context, we carried out a robust data mining of publicly available RNA-seq datasets to investigate the extension of heat-induced AS in rice leaves. For this, datasets of interest were subjected to filtering and quality control, followed by accurate transcript-specific quantifications. Powerful differential gene expression (DE) and differential AS (DAS) identified 17,143 and 2162 heat response genes, respectively, many of which are novel. Detailed analysis of DAS genes coding for key regulators of gene expression suggests that AS helps shape transcriptome and proteome diversity in response to heat. The knowledge resulting from this study confirmed a widespread transcriptional and post-transcriptional response to heat stress in plants, and it provided novel candidates for rapidly advancing rice breeding in response to climate change.

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops have historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over-reliant on a small number of internationally traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (particularly where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies, we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 Mb genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realization of the benefits to global nutrition security and agrobiodiversity.

Root growth responses to mechanical impedance are regulated by a network of ROS, ethylene and auxin signalling in Arabidopsis

The growth and development of root systems is influenced by mechanical properties of the substrate in which the plants grow. Mechanical impedance, such as by compacted soil, can reduce root elongation and limit crop productivity. To understand better the mechanisms involved in plant root responses to mechanical impedance stress, we investigated changes in the root transcriptome and hormone signalling responses of Arabidopsis to artificial root barrier systems in vitro. We demonstrate that upon encountering a barrier, reduced Arabidopsis root growth and a characteristic 'step-like' growth pattern is due to a reduction in cell elongation associated with changes in signalling gene expression. Data from RNA-sequencing combined with reporter line and mutant studies identified essential roles for reactive oxygen species, ethylene and auxin signalling during the barrier response. We propose a model in which early responses to mechanical impedance include reactive oxygen signalling integrated with ethylene and auxin responses to mediate root growth changes. Inhibition of ethylene responses allows improved growth in response to root impedance, an observation that may inform future crop breeding programmes.

De novo transcriptome assembly, functional annotation, and expression profiling of rye (Secale cereale L.) hybrids inoculated with ergot (Claviceps purpurea)

Rye is used as food, feed, and for bioenergy production and remain an essential grain crop for cool temperate zones in marginal soils. Ergot is known to cause severe problems in cross-pollinated rye by contamination of harvested grains. The molecular response of the underlying mechanisms of this disease is still poorly understood due to the complex infection pattern. RNA sequencing can provide astonishing details about the transcriptional landscape, hence we employed a transcriptomic approach to identify genes in the underlying mechanism of ergot infection in rye. In this study, we generated de novo assemblies from twelve biological samples of two rye hybrids with identified contrasting phenotypic responses to ergot infection. The final transcriptome of ergot susceptible (DH372) and moderately ergot resistant (Helltop) hybrids contain 208,690 and 192,116 contigs, respectively. By applying the BUSCO pipeline, we confirmed that these transcriptome assemblies contain more than 90% of gene representation of the available orthologue groups at Virdiplantae odb10. We employed a de novo assembled and the draft reference genome of rye to count the differentially expressed genes (DEGs) between the two hybrids with and without inoculation. The gene expression comparisons revealed that 228 genes were linked to ergot infection in both hybrids. The genome ontology enrichment analysis of DEGs associated them with metabolic processes, hydrolase activity, pectinesterase activity, cell wall modification, pollen development and pollen wall assembly. In addition, gene set enrichment analysis of DEGs linked them to cell wall modification and pectinesterase activity. These results suggest that a combination of different pathways, particularly cell wall modification and pectinesterase activity contribute to the underlying mechanism that might lead to resistance against ergot in rye. Our results may pave the way to select genetic material to improve resistance against ergot through better understanding of the mechanism of ergot infection at molecular level. Furthermore, the sequence data and de novo assemblies are valuable as scientific resources for future studies in rye.

Exploring transcriptional switches from pairwise, temporal and population RNA-Seq data using deepTS

Transcriptional switch (TS) is a widely observed phenomenon caused by changes in the relative expression of transcripts from the same gene, in spatial, temporal or other dimensions. TS has been associated with human diseases, plant development and stress responses. Its investigation is often hampered by a lack of suitable tools allowing comprehensive and flexible TS analysis for high-throughput RNA sequencing (RNA-Seq) data. Here, we present deepTS, a user-friendly web-based implementation that enables a fully interactive, multifunctional identification, visualization and analysis of TS events for large-scale RNA-Seq datasets from pairwise, temporal and population experiments. deepTS offers rich functionality to streamline RNA-Seq-based TS analysis for both model and non-model organisms and for those with or without reference transcriptome. The presented case studies highlight the capabilities of deepTS and demonstrate its potential for the transcriptome-wide TS analysis of pairwise, temporal and population RNA-Seq data. We believe deepTS will help research groups, regardless of their informatics expertise, perform accessible, reproducible and collaborative TS analyses of large-scale RNA-Seq data.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

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.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq

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.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

High-Resolution RT-PCR Analysis of Alternative Barley Transcripts

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.

Nuclear Speckle RNA Binding Proteins Remodel Alternative Splicing and the Non-coding Arabidopsis Transcriptome to Regulate a Cross-Talk Between Auxin and Immune Responses

Nuclear speckle RNA binding proteins (NSRs) act as regulators of alternative splicing (AS) and auxin-regulated developmental processes such as lateral root formation in Arabidopsis thaliana. These proteins were shown to interact with specific alternatively spliced mRNA targets and at least with one structured lncRNA, named Alternative Splicing Competitor RNA. Here, we used genome-wide analysis of RNAseq to monitor the NSR global role on multiple tiers of gene expression, including RNA processing and AS. NSRs affect AS of 100s of genes as well as the abundance of lncRNAs particularly in response to auxin. Among them, the FPA floral regulator displayed alternative polyadenylation and differential expression of antisense COOLAIR lncRNAs in nsra/b mutants. This may explains the early flowering phenotype observed in nsra and nsra/b mutants. GO enrichment analysis of affected lines revealed a novel link of NSRs with the immune response pathway. A RIP-seq approach on an NSRa fusion protein in mutant background identified that lncRNAs are privileged direct targets of NSRs in addition to specific AS mRNAs. The interplay of lncRNAs and AS mRNAs in NSR-containing complexes may control the crosstalk between auxin and the immune response pathway.

SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions

Despite the many approaches to study differential splicing from RNA-seq, many challenges remain unsolved, including computing capacity and sequencing depth requirements. Here we present SUPPA2, a new method that addresses these challenges, and enables streamlined analysis across multiple conditions taking into account biological variability. Using experimental and simulated data, we show that SUPPA2 achieves higher accuracy compared to other methods, especially at low sequencing depth and short read length. We use SUPPA2 to identify novel Transformer2-regulated exons, novel microexons induced during differentiation of bipolar neurons, and novel intron retention events during erythroblast differentiation.

A high quality Arabidopsis transcriptome for accurate transcript-level analysis of alternative splicing

Alternative splicing generates multiple transcript and protein isoforms from the same gene and thus is important in gene expression regulation. To date, RNA-sequencing (RNA-seq) is the standard method for quantifying changes in alternative splicing on a genome-wide scale. Understanding the current limitations of RNA-seq is crucial for reliable analysis and the lack of high quality, comprehensive transcriptomes for most species, including model organisms such as Arabidopsis, is a major constraint in accurate quantification of transcript isoforms. To address this, we designed a novel pipeline with stringent filters and assembled a comprehensive Reference Transcript Dataset for Arabidopsis (AtRTD2) containing 82,190 non-redundant transcripts from 34 212 genes. Extensive experimental validation showed that AtRTD2 and its modified version, AtRTD2-QUASI, for use in Quantification of Alternatively Spliced Isoforms, outperform other available transcriptomes in RNA-seq analysis. This strategy can be implemented in other species to build a pipeline for transcript-level expression and alternative splicing analyses.

Herboxidiene triggers splicing repression and abiotic stress responses in plants

Background

Constitutive and alternative splicing of pre-mRNAs from multiexonic genes controls the diversity of the proteome; these precisely regulated processes also fine-tune responses to cues related to growth, development, and stresses. Small-molecule inhibitors that perturb splicing provide invaluable tools for use as chemical probes to uncover the molecular underpinnings of splicing regulation and as potential anticancer compounds.

Results

Here, we show that herboxidiene (GEX1A) inhibits both constitutive and alternative splicing. Moreover, GEX1A activates genome-wide transcriptional patterns involved in abiotic stress responses in plants. GEX1A treatment -activated ABA-inducible promoters, and led to stomatal closure. Interestingly, GEX1A and pladienolide B (PB) elicited similar cellular changes, including alterations in the patterns of transcription and splicing, suggesting that these compounds might target the same spliceosome complex in plant cells.

Conclusions

Our study establishes GEX1A as a potent splicing inhibitor in plants that can be used to probe the assembly, dynamics, and molecular functions of the spliceosome and to study the interplay between splicing stress and abiotic stresses, as well as having potential biotechnological applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3656-z) contains supplementary material, which is available to authorized users.