Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short‐read RNA sequencing, single molecule long‐read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron‐containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non‐conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.     

Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis

Background

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is an important gene regulation process that potentially regulates many physiological processes in plants, including the response to abiotic stresses such as salt stress.

Results

To analyze global changes in AS under salt stress, we obtained high-coverage (~200 times) RNA sequencing data from Arabidopsis thaliana seedlings that were treated with different concentrations of NaCl. We detected that ~49% of all intron-containing genes were alternatively spliced under salt stress, 10% of which experienced significant differential alternative splicing (DAS). Furthermore, AS increased significantly under salt stress compared with under unstressed conditions. We demonstrated that most DAS genes were not differentially regulated by salt stress, suggesting that AS may represent an independent layer of gene regulation in response to stress. Our analysis of functional categories suggested that DAS genes were associated with specific functional pathways, such as the pathways for the responses to stresses and RNA splicing. We revealed that serine/arginine-rich (SR) splicing factors were frequently and specifically regulated in AS under salt stresses, suggesting a complex loop in AS regulation for stress adaptation. We also showed that alternative splicing site selection (SS) occurred most frequently at 4 nucleotides upstream or downstream of the dominant sites and that exon skipping tended to link with alternative SS.

Conclusions

Our study provided a comprehensive view of AS under salt stress and revealed novel insights into the potential roles of AS in plant response to salt stress.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-431) contains supplementary material, which is available to authorized users.     

Evolutionary impact of limited splicing fidelity in mammalian genes

The functional significance of most alternative splicing (AS) events, especially frame-shifting ones, has been controversial. Using human-mouse comparison, we demonstrate that frame-preserving AS events adapt and get fixed more rapidly than frame-shifting AS events; selection for smaller exon size is stronger in frame-preserving exons than in frame-shifting ones. These results suggest AS events introducing mild changes are generally favored during evolution and explain the excess of shorter, frame-preserving cassette exons in present mammalian genomes.     

Autoregulation of Fox protein expression to produce dominant negative splicing factors

The Fox proteins are a family of regulators that control the alternative splicing of many exons in neurons, muscle, and other tissues. Each of the three mammalian paralogs, Fox-1 (A2BP1), Fox-2 (RBM9), and Fox-3 (HRNBP3), produces proteins with a single RNA-binding domain (RRM) flanked by N- and C-terminal domains that are highly diversified through the use of alternative promoters and alternative splicing patterns. These genes also express protein isoforms lacking the second half of the RRM (FoxΔRRM), due to the skipping of a highly conserved 93-nt exon. Fox binding elements overlap the splice sites of these exons in Fox-1 and Fox-2, and the Fox proteins themselves inhibit exon inclusion. Unlike other cases of splicing autoregulation by RNA-binding proteins, skipping the RRM exon creates an in-frame deletion in the mRNA to produce a stable protein. These FoxΔRRM isoforms expressed from cDNA exhibit highly reduced binding to RNA in vivo. However, we show that they can act as repressors of Fox-dependent splicing, presumably by competing with full-length Fox isoforms for interaction with other splicing factors. Interestingly, the Drosophila Fox homolog contains a nearly identical exon in its RRM domain that also has flanking Fox-binding sites. Thus, rather than autoregulation of splicing controlling the abundance of the regulator, the Fox proteins use a highly conserved mechanism of splicing autoregulation to control production of a dominant negative isoform.     

Differential Expression of the Peripheral Benzodiazepine Receptor and Gremlin During Adipogenesis

This study used the mRNA differential display technique to identify differentially expressed genes during the process of adipogenesis in the preadipocyte cell line, 3T3‐L1. 3T3‐L1 cells were treated with dexamethasone, isobutyl‐1‐methylxanthine, and insulin to induce differentiation into mature adipocytes. Cells were collected at three time‐points during differentiation: Day 0 (d0), or nondifferentiated; Day 3 (d3), during differentiation; and Day 10 (d10), >90% of the cells had differentiated into mature adipocytes. Initial studies yielded 18 potentially differentially regulated cDNA candidates (8 down‐regulated and 10 up‐regulated). Reverse Northern and Northern blots confirmed differential expression of six of the candidates. Four of the candidates up‐regulated on d3 and d10 were identified by sequence analysis to be lipoprotein lipase, a well‐known marker of adipocyte differentiation. A fifth candidate that was expressed in d0, but not d3 or d10, was identified as DRM/gremlin, a bone morphogenetic protein antagonist. Finally, a sixth candidate that was increased at d3 and d10 was identified as the peripheral benzodiazepine receptor, which has been implicated in proliferation, differentiation, and cholesterol transport in cells. This study is the first to show that peripheral benzodiazepine receptor and DRM/gremlin are expressed in preadipocyte cell lines and that they are differentially regulated during adipogenesis.     

Splicing factor and exon profiling across human tissues

It has been shown that alternative splicing is especially prevalent in brain and testis when compared to other tissues. To test whether there is a specific propensity of these tissues to generate splicing variants, we used a single source of high-density microarray data to perform both splicing factor and exon expression profiling across 11 normal human tissues. Paired comparisons between tissues and an original exon-based statistical group analysis demonstrated after extensive RT-PCR validation that the cerebellum, testis, and spleen had the largest proportion of differentially expressed alternative exons. Variations at the exon level correlated with a larger number of splicing factors being expressed at a high level in the cerebellum, testis and spleen than in other tissues. However, this splicing factor expression profile was similar to a more global gene expression pattern as a larger number of genes had a high expression level in the cerebellum, testis and spleen. In addition to providing a unique resource on expression profiling of alternative splicing variants and splicing factors across human tissues, this study demonstrates that the higher prevalence of alternative splicing in a subset of tissues originates from the larger number of genes, including splicing factors, being expressed than in other tissues.     

Differential glycosylation regulates processing of lipoprotein receptors by gamma-secretase

The low density lipoprotein (LDL) receptor-related protein 1 (LRP1) belongs to a growing number of cell surface proteins that undergo regulated proteolytic processing that culminates in the release of their intracellular domain (ICD) by the intramembranous protease gamma-secretase. Here we show that LRP1 is differentially glycosylated in a tissue-specific manner and that carbohydrate addition reduces proteolytic cleavage of the extracellular domain and, concomitantly, ICD release. The apolipoprotein E (apoE) receptor-2 (apoER2), another member of the LDL receptor family with functions in cellular signal transmission, also undergoes sequential proteolytic processing, resulting in intracellular domain release into the cytoplasm. The penultimate processing step also involves cleavage of the apoER2 extracellular domain. The rate at which this cleavage step occurs is determined by the glycosylation state of the receptor, which in turn is regulated by the alternative splicing of an exon encoding several O-linked sugar attachment sites. These findings suggest a role for differential and tissue-specific glycosylation as a physiological switch that modulates the diverse biological functions of these receptors in a cell-type specific manner.