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.     

Identification of an intronic splicing regulatory element involved in auto‐regulation of alternative splicing of SCL33 pre‐mRNA

In Arabidopsis, pre‐mRNAs of serine/arginine‐rich (SR) proteins undergo extensive alternative splicing (AS). However, little is known about the cis‐elements and trans‐acting proteins involved in regulating AS. Using a splicing reporter (GFP–intron–GFP), consisting of the GFP coding sequence interrupted by an alternatively spliced intron of SCL33, we investigated whether cis‐elements within this intron are sufficient for AS, and which SR proteins are necessary for regulated AS. Expression of the splicing reporter in protoplasts faithfully produced all splice variants from the intron, suggesting that cis‐elements required for AS reside within the intron. To determine which SR proteins are responsible for AS, the splicing pattern of the GFP–intron–GFP reporter was investigated in protoplasts of three single and three double mutants of SR genes. These analyses revealed that SCL33 and a closely related paralog, SCL30a, are functionally redundant in generating specific splice variants from this intron. Furthermore, SCL33 protein bound to a conserved sequence in this intron, indicating auto‐regulation of AS. Mutations in four GAAG repeats within the conserved region impaired generation of the same splice variants that are affected in the scl33 scl30a double mutant. In conclusion, we have identified the first intronic cis‐element involved in AS of a plant SR gene, and elucidated a mechanism for auto‐regulation of AS of this intron.     

Alternative splicing coupled mRNA decay shapes the temperature‐dependent transcriptome

Mammalian body temperature oscillates with the time of the day and is altered in diverse pathological conditions. We recently identified a body temperature‐sensitive thermometer‐like kinase, which alters SR protein phosphorylation and thereby globally controls alternative splicing (AS). AS can generate unproductive variants which are recognized and degraded by diverse mRNA decay pathways—including nonsense‐mediated decay (NMD). Here we show extensive coupling of body temperature‐controlled AS to mRNA decay, leading to global control of temperature‐dependent gene expression (GE). Temperature‐controlled, decay‐inducing splicing events are evolutionarily conserved and pervasively found within RNA‐binding proteins, including most SR proteins. AS‐coupled poison exon inclusion is essential for rhythmic GE of SR proteins and has a global role in establishing temperature‐dependent rhythmic GE profiles, both in mammals under circadian body temperature cycles and in plants in response to ambient temperature changes. Together, these data identify body temperature‐driven AS‐coupled mRNA decay as an evolutionary ancient, core clock‐independent mechanism to generate rhythmic GE.     

Global genome splicing analysis reveals an increased number of alternatively spliced genes with aging

Alternative splicing (AS) is a key regulatory mechanism for the development of different tissues; however, not much is known about changes to alternative splicing during aging. Splicing events may become more frequent and widespread genome‐wide as tissues age and the splicing machinery stringency decreases. Using skin, skeletal muscle, bone, thymus, and white adipose tissue from wild‐type C57BL6/J male mice (4 and 18 months old), we examined the effect of age on splicing by AS analysis of the differential exon usage of the genome. The results identified a considerable number of AS genes in skeletal muscle, thymus, bone, and white adipose tissue between the different age groups (ranging from 27 to 246 AS genes corresponding to 0.3–3.2% of the total number of genes analyzed). For skin, skeletal muscle, and bone, we included a later age group (28 months old) that showed that the number of alternatively spliced genes increased with age in all three tissues (P < 0.01). Analysis of alternatively spliced genes across all tissues by gene ontology and pathway analysis identified 158 genes involved in RNA processing. Additional analysis of AS in a mouse model for the premature aging disease Hutchinson–Gilford progeria syndrome was performed. The results show that expression of the mutant protein, progerin, is associated with an impaired developmental splicing. As progerin accumulates, the number of genes with AS increases compared to in wild‐type skin. Our results indicate the existence of a mechanism for increased AS during aging in several tissues, emphasizing that AS has a more important role in the aging process than previously known.     

The mitochondrial complexome of Arabidopsis thaliana

Mitochondria are central to cellular metabolism and energy conversion. In plants they also enable photosynthesis through additional components and functional flexibility. A majority of those processes relies on the assembly of individual proteins to larger protein complexes, some of which operate as large molecular machines. There has been a strong interest in the makeup and function of mitochondrial protein complexes and protein–protein interactions in plants, but the experimental approaches used typically suffer from selectivity or bias. Here, we present a complexome profiling analysis for leaf mitochondria of the model plant Arabidopsis thaliana for the systematic characterization of protein assemblies. Purified organelle extracts were separated by 1D Blue native (BN) PAGE, a resulting gel lane was dissected into 70 slices (complexome fractions) and proteins in each slice were identified by label free quantitative shot‐gun proteomics. Overall, 1359 unique proteins were identified, which were, on average, present in 17 complexome fractions each. Quantitative profiles of proteins along the BN gel lane were aligned by similarity, allowing us to visualize protein assemblies. The data allow re‐annotating the subunit compositions of OXPHOS complexes, identifying assembly intermediates of OXPHOS complexes and assemblies of alternative respiratory oxidoreductases. Several protein complexes were discovered that have not yet been reported in plants, such as a 530 kDa Tat complex, 460 and 1000 kDa SAM complexes, a calcium ion uniporter complex (150 kDa) and several PPR protein complexes. We have set up a tailored online resource ( https://complexomemap.de/at_mito_leaves ) to deposit the data and to allow straightforward access and custom data analyses.     

AtBUD13 affects pre‐mRNA splicing and is essential for embryo development in Arabidopsis

Pre‐mRNA splicing is an important step for gene expression regulation. Yeast Bud13p (bud‐site selection protein 13) regulates the budding pattern and pre‐mRNA splicing in yeast cells; however, no Bud13p homologs have been identified in plants. Here, we isolated two mutants that carry T‐DNA insertions at the At1g31870 locus and shows early embryo lethality and seed abortion. At1g31870 encodes an Arabidopsis homolog of yeast Bud13p, AtBUD13. Although AtBUD13 homologs are widely distributed in eukaryotic organisms, phylogenetic analysis revealed that their protein domain organization is more complex in multicellular species. AtBUD13 is expressed throughout plant development including embryogenesis and AtBUD13 proteins is localized in the nucleus in Arabidopsis. RNA‐seq analysis revealed that AtBUD13 mutation predominantly results in the intron retention, especially for shorter introns (≤100 bases). Within this group of genes, we identified 52 genes involved in embryogenesis, out of which 22 are involved in nucleic acid metabolism. Our results demonstrate that AtBUD13 plays critical roles in early embryo development by effecting pre‐mRNA splicing.