Genome-wide analysis of alternative splicing landscapes modulated during plant-virus interactions in Brachypodium distachyon

Evolutionarily Conserved Alternative Splicing Across Monocots

One difficulty when identifying alternative splicing (AS) events in plants is distinguishing functional AS from splicing noise. One way to add confidence to the validity of a splice isoform is to observe that it is conserved across evolutionarily related species. We use a high throughput method to identify junction-based conserved AS events from RNA-Seq data across nine plant species, including five grass monocots (maize, sorghum, rice, Brachpodium, and foxtail millet), plus two nongrass monocots (banana and African oil palm), the eudicot Arabidopsis, and the basal angiosperm Amborella In total, 9804 AS events were found to be conserved between two or more species studied. In grasses containing large regions of conserved synteny, the frequency of conserved AS events is twice that observed for genes outside of conserved synteny blocks. In plant-specific RS and RS2Z subfamilies of the serine/arginine (SR) splice-factor proteins, we observe both conservation and divergence of AS events after the whole genome duplication in maize. In addition, plant-specific RS and RS2Z splice-factor subfamilies are highly connected with R2R3-MYB in STRING functional protein association networks built using genes exhibiting conserved AS. Furthermore, we discovered that functional protein association networks constructed around genes harboring conserved AS events are enriched for phosphatases, kinases, and ubiquitylation genes, which suggests that AS may participate in regulating signaling pathways. These data lay the foundation for identifying and studying conserved AS events in the monocots, particularly across grass species, and this conserved AS resource identifies an additional layer between genotype to phenotype that may impact future crop improvement efforts.

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

Genome-wide analysis of specific alterations in transcript structure and accumulation caused by nutrient deficiencies in Arabidopsis thaliana

The alteration of transcript structure contributes to transcriptome plasticity. In this study, we analyzed the genome-wide response of exon combination patterns to deficiencies in 12 different nutrients in Arabidopsis thaliana roots. RNA sequencing analysis and bioinformatics using a simulation survey revealed more than 600 genes showing varying exon combinations. The overlap between genes showing differential expression (DE) and genes showing differential exon combination (DC) was notably low. Additionally, gene ontology analysis showed that gene functions were not shared between the DE and DC genes, suggesting that the genes showing DC had different roles than those showing DE. Most of the DC genes were nutrient specific. For example, two homologs of the MYB transcription factor genes MYB48 and MYB59 showed differential alternative splicing only in response to low levels of potassium. Alternative splicing of those MYB genes modulated DNA-binding motifs, and MYB59 is reportedly involved in the inhibition of root elongation. Therefore, the increased abundance of MYB isoforms with an intact DNA-binding motif under low potassium may be involved in the active inhibition of root elongation. Overall, we provide global and comprehensive data for DC genes affected by nutritional deficiencies, which contribute to elucidating an unknown mechanism involved in adaptation to nutrient deficiency.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)

Background

Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.

Results

In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.

Conclusions

These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Electronic supplementary material

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

Complete nucleotide sequences and virion particle association of two satellite RNAs of panicum mosaic virus

Over six decades ago, panicum mosaic virus (PMV) was identified as the first viral pathogen of cultivated switchgrass (Panicum virgatum). Subsequently, PMV was demonstrated to support the replication of both a satellite RNA virus (SPMV) and satellite RNA (satRNA) agents during natural infections of host grasses. In this study, we report the isolation and full-length sequences of two PMV satRNAs identified in 1988 from St. Augustinegrass (Stenotaphrum secundatum) and centipedegrass (Eremochloa ophiuroides) hosts. Each of these satellites have sequence relatedness at their 5'- and 3'-ends. In addition, satC has a region of ∼100 nt complementary to the 3'-end of the PMV genome. These agents are associated with purified virions of SPMV infections. Additionally, satS and satC RNAs contain conserved in-frame open reading frames in the complementary-sense sequences that could potentially generate 6.6- and 7.9-kDa proteins, respectively. In protoplasts and plants satS is infectious, when co-inoculated with the PMV RNA alone or PMV+SPMV RNAs, and negatively affects their accumulation.

Comprehensive Transcriptome Analyses Reveal that Potato Spindle Tuber Viroid Triggers Genome-Wide Changes in Alternative Splicing, Inducible trans-Acting Activity of Phased Secondary Small Interfering RNAs, and Immune Responses

Many pathogens express noncoding RNAs (ncRNAs) during infection processes. In the most extreme case, pathogenic ncRNAs alone (such as viroids) can infect eukaryotic organisms, leading to diseases. While a few pathogenic ncRNAs have been implicated in regulating gene expression, the functions of most pathogenic ncRNAs in host-pathogen interactions remain unclear. Here, we employ potato spindle tuber viroid (PSTVd) infecting tomato as a system to dissect host interactions with pathogenic ncRNAs, using comprehensive transcriptome analyses. We uncover various new activities in regulating gene expression during PSTVd infection, such as genome-wide alteration in alternative splicing of host protein-coding genes, enhanced guided-cleavage activities of a host microRNA, and induction of the trans-acting function of phased secondary small interfering RNAs. Furthermore, we reveal that PSTVd infection massively activates genes involved in plant immune responses, mainly those in the calcium-dependent protein kinase and mitogen-activated protein kinase cascades, as well as prominent genes involved in hypersensitive responses, cell wall fortification, and hormone signaling. Intriguingly, our data support a notion that plant immune systems can respond to pathogenic ncRNAs, which has broad implications for providing new opportunities for understanding the complexity of immune systems in differentiating "self" and "nonself," as well as lay the foundation for resolving the long-standing question regarding the pathogenesis mechanisms of viroids and perhaps other infectious RNAs.IMPORTANCE Numerous pathogens, including viruses, express pathogenic noncoding transcripts during infection. In the most extreme case, pathogenic noncoding RNAs alone (i.e., viroids) can cause disease in plants. While some work has demonstrated that pathogenic noncoding RNAs interact with host factors for function, the biological significance of pathogenic noncoding RNAs in host-pathogen interactions remains largely unclear. Here, we apply comprehensive genome-wide analyses of plant-viroid interactions and discover several novel molecular activities underlying nuclear-replicating viroid infection processes in plants, including effects on the expression and function of host noncoding transcripts, as well as the alternative splicing of host protein-coding genes. Importantly, we show that plant immunity is activated upon infection of a nuclear-replicating viroid, which is a new concept that helps to understand viroid-based pathogenesis. Our finding has broad implications for understanding the complexity of host immune systems and the diverse functions of noncoding RNAs.

Prediction and feature analysis of intron retention events in plant genome

Alternative splicing (AS) is a major contributor to increase the potential informational content of eukaryotic genomes by creating multiple mRNA species and proteins from a single gene. In plants, up to 60% genes are alternatively spliced and the most common type of AS is intron retention (IR). Genomic analyses of IR have illuminated its crucial role in shaping the evolution of genomes, in the control of developmental processes, and in the dynamic regulation of the transcriptome to influence phenotype. To explore the relationship between the sequence feature and the formation mechanism of IR, we statistically analyzed the retained introns and proposed an improved random forest-based hybrid method to predict intron retention events in plant genome. The results indicate that IR has significant relationship with individual introns which have weaker 5' splice sites, lower GC content and less termination codon occurrence. By the method we proposed, 93.48% retained introns can be correctly distinguished from constitutive introns. Strikingly, our study will facilitate a better understanding of underlying mechanisms involved in intron retention.

Global identification of alternative splicing via comparative analysis of SMRT- and Illumina-based RNA-seq in strawberry

Alternative splicing (AS) is a key post-transcriptional regulatory mechanism, yet little information is known about its roles in fruit crops. Here, AS was globally analyzed in the wild strawberry Fragaria vesca genome with RNA-seq data derived from different stages of fruit development. The AS landscape was characterized and compared between the single-molecule, real-time (SMRT) and Illumina RNA-seq platform. While SMRT has a lower sequencing depth, it identifies more genes undergoing AS (57.67% of detected multiexon genes) when it is compared with Illumina (33.48%), illustrating the efficacy of SMRT in AS identification. We investigated different modes of AS in the context of fruit development; the percentage of intron retention (IR) is markedly reduced whereas that of alternative acceptor sites (AA) is significantly increased post-fertilization when compared with pre-fertilization. When all the identified transcripts were combined, a total of 66.43% detected multiexon genes in strawberry undergo AS, some of which lead to a gain or loss of conserved domains in the gene products. The work demonstrates that SMRT sequencing is highly powerful in AS discovery and provides a rich data resource for later functional studies of different isoforms. Further, shifting AS modes may contribute to rapid changes of gene expression during fruit set.

Response of Gene Expression and Alternative Splicing to Distinct Growth Environments in Tomato

Phenotypic plasticity is the phenomenon that one particular genotype produces different phenotypes under different environmental conditions, but its underlying molecular and genetic mechanisms are poorly understood. Plastic traits may be under the control of genes whose expression is modulated by environmental cues. In this study, we investigated phenotypic plasticity in tomato (Solanum lycopersicum) and its ancestral species S. pimpinellifolium by comparing the global gene expression of young seedlings grown under two distinct growth conditions. Our results show that more than 7000 genes exhibited differential expression in response to environmental changes from phytotron to a plastic greenhouse, and 98 environmentally sensitive genes displayed the same patterns of expression response across the two tomato species. We also found that growth conditions had a remarkable impact on transcriptome complexity, attributable to alternative splicing (AS), in which 665 splice variants showed differential expression in response to the environmental changes. Moreover, more splice variants and AS events per gene were detected in plastic greenhouse-grown seedlings than their phytotron counterparts, and these seedlings also had higher percentages of intron retention events. The identification of the conserved environmentally-sensitive genes and the splice variants in this study will be useful for further analysis of gene regulation of environmental response in tomato and other crops.

Regulation of FT splicing by an endogenous cue in temperate grasses

Appropriate flowering timing is crucial for plant reproductive success. The florigen, FLOWERING LOCUS T (FT), interacts with 14-3-3 proteins and the bZIP transcription factor FD, functioning at core nodes in multiple flowering pathways. There are two FT homologues, FT1 and FT2, in Brachypodium distachyon. Here we show that FT2 undergoes age-dependent alternative splicing (AS), resulting in two splice variants (FT2α and FT2β). The FT2β-encoded protein cannot interact with FD or 14-3-3s but is able to form heterodimers with FT2α and FT1, thereby interfering with the florigen-mediated assembly of the flowering initiation complex. Notably, transgenic plants overproducing FT2β exhibit delayed flowering, while transgenic plants in which FT2β is silenced by an artificial microRNA display accelerated flowering, demonstrating a dominant-negative role of FT2β in flowering induction. Furthermore, we show that the AS splicing of FT2 is conserved in important cereal crops, such as barley and wheat. Collectively, these findings reveal a novel posttranscriptional mode of FT regulation in temperate grasses.

Alternative Splicing of Barley Clock Genes in Response to Low Temperature

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.

Dynamic transcriptome profiling of Bean Common Mosaic Virus (BCMV) infection in Common Bean (Phaseolus vulgaris L.)

BACKGROUND

Bean common mosaic virus (BCMV) is widespread, with Phaseolus species as the primary host plants. Numerous BCMV strains have been identified on the basis of a panel of bean varieties that distinguish the pathogenicity types with respect to the viral strains. The molecular responses in Phaseolus to BCMV infection have not yet been well characterized.

RESULTS

We report the transcriptional responses of a widely susceptible variety of common bean (Phaseolus vulgaris L., cultivar 'Stringless green refugee') to two BCMV strains, in a time-course experiment. We also report the genome sequence of a previously unreported BCMV strain. The interaction with the known strain NL1-Iowa causes moderate symptoms and large transcriptional responses, and the newly identified strain (Strain 2 or S2) causes severe symptoms and moderate transcriptional responses. The transcriptional profiles of host plants infected with the two isolates are distinct, and involve numerous differences in splice forms in particular genes, and pathway specific expression patterns.

CONCLUSIONS

We identified differential host transcriptome response after infection of two different strains of Bean common mosaic virus (BCMV) in common bean (Phaseolus vulgaris L.). Virus infection initiated a suite of changes in gene expression level and patterns in the host plants. Pathways related to defense, gene regulation, metabolic processes, photosynthesis were specifically altered after virus infection. Results presented in this study can increase the understanding of host-pathogen interactions and provide resources for further investigations of the biological mechanisms in BCMV infection and defense.

A survey of the sorghum transcriptome using single-molecule long reads

Alternative splicing and alternative polyadenylation (APA) of pre-mRNAs greatly contribute to transcriptome diversity, coding capacity of a genome and gene regulatory mechanisms in eukaryotes. Second-generation sequencing technologies have been extensively used to analyse transcriptomes. However, a major limitation of short-read data is that it is difficult to accurately predict full-length splice isoforms. Here we sequenced the sorghum transcriptome using Pacific Biosciences single-molecule real-time long-read isoform sequencing and developed a pipeline called TAPIS (Transcriptome Analysis Pipeline for Isoform Sequencing) to identify full-length splice isoforms and APA sites. Our analysis reveals transcriptome-wide full-length isoforms at an unprecedented scale with over 11,000 novel splice isoforms. Additionally, we uncover APA of ∼11,000 expressed genes and more than 2,100 novel genes. These results greatly enhance sorghum gene annotations and aid in studying gene regulation in this important bioenergy crop. The TAPIS pipeline will serve as a useful tool to analyse Iso-Seq data from any organism.

Alternative splicing and alternative polyadenylation (APA) contribute to mRNA diversity but are difficult to assess using short read RNA-seq data. Here, the authors use single molecule long-read isoform sequencing and develop a computational pipeline to identify full-length splice isoforms and APA sites in sorghum.

Genome-wide transcriptomic analysis reveals correlation between higher WRKY61 expression and reduced symptom severity in Turnip crinkle virus infected Arabidopsis thaliana

Turnip crinkle virus (TCV) is a carmovirus that infects many Arabidopsis ecotypes. Most studies mainly focused on discovery of resistance genes against TCV infection, and there is no Next Generation Sequencing based comparative genome wide transcriptome analysis reported. In this study, RNA-seq based transcriptome analysis revealed that 238 (155 up-regulated and 83 down-regulated) significant differentially expressed genes with at least 15-fold change were determined. Fifteen genes (including upregulated, unchanged and downregulated) were selected for RNA-seq data validation using quantitative real-time PCR, which showed consistencies between these two sets of data. GO enrichment analysis showed that numerous terms such as stress, immunity, defence and chemical stimulus were affected in TCV-infected plants. One putative plant defence related gene named WRKY61 was selected for further investigation. It showed that WRKY61 overexpression plants displayed reduced symptoms and less virus accumulation, as compared to wild type (WT) and WRKY61 deficient lines, suggesting that higher WRKY61 expression level reduced TCV viral accumulation. In conclusion, our transcriptome analysis showed that global gene expression was detected in TCV-infected Arabidopsis thaliana. WRKY61 gene was shown to be negatively correlated with TCV infection and viral symptoms, which may be connected to plant immunity pathways.

Satellite RNAs and Satellite Viruses

Satellite RNAs and satellite viruses are extraviral components that can affect either the pathogenicity, the accumulation, or both of their associated viruses while themselves being dependent on the associated viruses as helper viruses for their infection. Most of these satellite RNAs are noncoding RNAs, and in many cases, have been shown to alter the interaction of their helper viruses with their hosts. In only a few cases have the functions of these satellite RNAs in such interactions been studied in detail. In particular, work on the satellite RNAs of Cucumber mosaic virus and Turnip crinkle virus have provided novel insights into RNAs functioning as noncoding RNAs. These effects are described and potential roles for satellite RNAs in the processes involved in symptom intensification or attenuation are discussed. In most cases, models describing these roles involve some aspect of RNA silencing or its suppression, either directly or indirectly involving the particular satellite RNA.

Evolutionary Character of Alternative Splicing in Plants

Alternative splicing (AS) is one of the most important ways to enhance the functional diversity of genes. Huge amounts of data have been produced by microarray, expressed sequence tag, and RNA-seq, and plenty of methods have been developed specifically for this task. The most frequently asked questions in previous research were as follows. What is the content rate of AS genes among the whole gene set? How many AS types are presented in the genome, and which type is dominant? How about the conservation ability of AS among different species? Which kinds of isoforms from some genes have the environmental response to help individual adaptation? Based on this background, we collected analysis results from 17 species to try to map out the landscape of AS studies in plants. We have noted the shortages of previous results, and we appeal to all scientists working in the AS field to make a standard protocol so that analyses between different projects are comparable.

Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize

Alternative splicing plays a crucial role in plant development as well as stress responses. Although alternative splicing has been studied during development and in response to stress, the interplay between these two factors remains an open question. To assess the effects of drought stress on developmentally regulated splicing in maize (Zea mays), 94 RNA-seq libraries from ear, tassel, and leaf of the B73 public inbred line were constructed at four developmental stages under both well-watered and drought conditions. This analysis was supplemented with a publicly available series of 53 libraries from developing seed, embryo, and endosperm. More than 48,000 novel isoforms, often with stage- or condition-specific expression, were uncovered, suggesting that developmentally regulated alternative splicing occurs in thousands of genes. Drought induced large developmental splicing changes in leaf and ear but relatively few in tassel. Most developmental stage-specific splicing changes affected by drought were tissue dependent, whereas stage-independent changes frequently overlapped between leaf and ear. A linear relationship was found between gene expression changes in splicing factors and alternative spicing of other genes during development. Collectively, these results demonstrate that alternative splicing is strongly associated with tissue type, developmental stage, and stress condition.

Heat Stress Regulates the Expression of Genes at Transcriptional and Post-Transcriptional Levels, Revealed by RNA-seq in Brachypodium distachyon

Heat stress greatly affects plant growth/development and influences the output of crops. With the increased occurrence of extreme high temperature, the negative influence on cereal products from heat stress becomes severer and severer. It is urgent to reveal the molecular mechanism in response to heat stress in plants. In this research, we used RNA-seq technology to identify differentially expressed genes (DEGs) in leaves of seedlings, leaves and inflorescences at heading stage of Brachypodium distachyon, one model plant of grasses. Results showed many genes in responding to heat stress. Of them, the expression level of 656 DEGs were altered in three groups of samples treated with high temperature. Gene ontology (GO) analysis showed that the highly enriched DEGs were responsible for heat stress and protein folding. According to KEGG pathway analysis, the DEGs were related mainly to photosynthesis-antenna proteins, the endoplasmic reticulum, and the spliceosome. Additionally, the expression level of 454 transcription factors belonging to 49 gene families was altered, as well as 1,973 splicing events occurred after treatment with high temperature. This research lays a foundation for characterizing the molecular mechanism of heat stress response and identifying key genes for those responses in plants. These findings also clearly show that heat stress regulates the expression of genes not only at transcriptional level, but also at post-transcriptional level.

Comprehensive characterization of a time-course transcriptional response induced by autotoxins in Panax ginseng using RNA-Seq

BACKGROUND

As a valuable medicinal plant, the yield of Panax ginseng is seriously affected by autotoxicity, which is a common phenomenon due to continuous cropping. However, the mechanism of autotoxicity in P. ginseng is still unknown.

RESULTS

In total, high throughput sequencing of 18 RNA-Seq libraries produced 996,000,000 100-nt reads that were assembled into 72,732 contigs. Compared with control, 3697 and 2828 genes were significantly up- and down-regulated across different tissues and time points, respectively. Gene Ontology enrichment analysis showed that 'enzyme inhibitor activity', 'carboxylesterase activity', 'pectinesterase activity', 'centrosome cycle and duplication' and 'mitotic spindle elongation' were enriched for the up-regulated genes. Transcription factors including AP2s/ERFs, MYBs, and WRKYs were up-regulated in roots after benzoic acid treatment. Moreover, reactive oxygen species, peroxidases and superoxide dismutase contigs were up-regulated in roots after benzoic acid treatment. Physiological and biochemical indexes showed that the proline and malondialdehyde content were restored to lower levels at a later stage after benzoic acid treatment. Benzoic acid inhibited the root hair development in a dose-dependent manner, and several differential expressed genes potentially involved in hair development were identified. Several key contigs in the flavonoid and ginsenoside biosynthesis pathways were repressed. Finally, 58,518 alternative splicing (AS) events from 12,950 genes were found after benzoic acid treatment. Interestingly, contigs in the ginsenoside biosynthetic pathway underwent AS, providing useful information about post-transcriptional regulation in P. ginseng.

CONCLUSIONS

This study revealed the stress-response molecular mechanisms in P. ginseng induced by benzoic acid.

Identification of alternative splicing events by RNA sequencing in early growth tomato fruits

BACKGROUND

Alternative splicing (AS) regulates multiple biological processes including flowering, circadian and stress response in plant. Although accumulating evidences indicate that AS is developmentally regulated, how AS responds to developmental cues is not well understood. Early fruit growth mainly characterized by active cell division and cell expansion contributes to the formation of fruit morphology and quality traits. Transcriptome profiling has revealed the coordinated complex regulation of gene expression in the process. High throughput RNA sequencing (RNA-seq) technology is advancing the genome-wide analysis of AS events in plant species, but the landscape of AS in early growth fruit is still not available for tomato (Solanum lycopersicum), a model plant for fleshy fruit development study.

RESULTS

Using RNA-seq, we surveyed the AS patterns in tomato seedlings, flowers and young developing fruits and found that 59.3 % of expressed multi-exon genes underwent AS in these tissues. The predominant type of AS events is intron retention, followed by alternative splice donor and acceptor, whereas exon skipping has the lowest frequency. Although the frequencies of AS events are similar among seedlings, flowers and early growth fruits, the fruits generated more splice variants per gene. Further comparison of gene expression in early growth fruits at 2, 5 and 10 days post anthesis revealed that 5206 multi-exon genes had at least one splice variants differentially expressed during early fruit development, whereas only 1059 out of them showed differential expression at gene level. We also identified 27 multi-exon genes showing differential splicing during early fruit growth. In addition, the study discovered 2507 new transcription regions (NTRs) unlinked to the annotated chromosomal regions, from where 956 putative protein coding transcripts and 1690 putative long non-coding RNAs were identified.

CONCLUSIONS

Our genome-wide analysis of AS events reveals a distinctive AS pattern in early growth tomato fruits. The landscape of AS obtained in this study will facilitate future investigation on transcriptome complexity and AS regulation during early fruit growth in tomato. The newly found NTRs will also be useful for updating the tomato genome annotation.

AtRTD - a comprehensive reference transcript dataset resource for accurate quantification of transcript-specific expression in Arabidopsis thaliana

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.

Characterization of SCL33 splicing patterns during diverse virus infections in Brachypodium distachyon

In eukaryotes alternative splicing (AS) influences transcriptome and proteome diversity. The mechanism and the genetic components mediating AS during plant-virus interactions are not known. Using RNA sequencing approaches, we recently analyzed the global AS changes occurring in Brachypodium distachyon (Brachypodium) during infections of Panicum mosaic virus (PMV) and its satellite virus (SPMV). We reported AS of defense-related genes including receptor-like kinases, NB-LRR proteins and transcription factors. Strikingly, multiple spliceosome components are themselves alternatively spliced during PMV and SPMV infections. Here, we analyzed the temporal splicing patterns of a splicing factor, Bd-SCL33, following infection of Brachypodium with 6 additional viruses in diverse genera. Our results reveal both dynamic and conserved expression patterns of Bd-SCL33 splice variants during virus infection, and implicate Bd-SCL33 function in response to biotic stresses.

Genomic architecture and functional relationships of intronless, constitutively- and alternatively-spliced genes in Brachypodium distachyon

Splicing and alternative splicing (AS) are widespread co- and post-transcriptional regulatory processes in plants. Recently, we characterized genome-wide AS landscapes and virus-induced AS patterns in Brachypodium distachyon (Brachypodium), a C3 model grass. Brachypodium plants infected with Panicum mosaic virus (PMV) alone or in mixed infections with its satellite virus (SPMV) were used for high-throughput, paired-end RNA sequencing. Here, using gene attributes of ∼5,655 intronless genes, ∼13,302 constitutively spliced, and ∼7,564 alternatively spliced genes, we analyzed the influence of genomic features on splicing incidence and AS frequency. In Brachypodium, gene length, coding sequence length, and exon and intron number were positively correlated to splicing incidence and AS frequency. In contrast, exon length and the percentage composition of GC (%GC) content were inversely correlated with splicing incidence and AS frequency. Although gene expression status had little correlation with splicing occurrence per se, it negatively correlated to AS frequency: i.e., genes with ≥5 alternatively spliced transcripts were significantly less expressed compared to genes encoding <5 alternative transcripts. Further gene set enrichment analysis uncovered unique functional relationships among nonspliced, constitutively spliced and alternatively spliced genes.