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.

Splicing of precursor messenger RNAs (pre-mRNAs) is an essential step in the expression of most eukaryotic genes. Both constitutive splicing and alternative splicing, which produces multiple messenger RNA (mRNA) isoforms from a single primary transcript, are modulated by reversible protein phosphorylation. Although the plant splicing machinery is known to be a target for phosphorylation, the protein kinases involved remain to be fully defined. We report here the identification of pre-mRNA processing 4 (PRP4) KINASE A (PRP4KA) in a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana (Arabidopsis). Prp4 kinase is the first spliceosome-associated kinase shown to regulate splicing in fungi and mammals but it has not yet been studied in plants. In the same screen we identified mutants defective in SAC3A, a putative mRNA export factor that is highly coexpressed with PRP4KA in Arabidopsis. Whereas the sac3a mutants appear normal, the prp4ka mutants display a pleiotropic phenotype featuring atypical rosettes, late flowering, tall final stature, reduced branching, and lowered seed set. Analysis of RNA-sequencing data from prp4ka and sac3a mutants identified widespread and partially overlapping perturbations in alternative splicing in the two mutants. Quantitative phosphoproteomic profiling of a prp4ka mutant detected phosphorylation changes in several serine/arginine-rich proteins, which regulate constitutive and alternative splicing, and other splicing-related factors. Tests of PRP4KB, the paralog of PRP4KA, indicated that the two genes are not functionally redundant. The results demonstrate the importance of PRP4KA for alternative splicing and plant phenotype, and suggest that PRP4KA may influence alternative splicing patterns by phosphorylating a subset of splicing regulators.

To investigate factors influencing pre-mRNA splicing in plants, we conducted a forward genetic screen using an alternatively-spliced GFP reporter gene in Arabidopsis thaliana. This effort generated a collection of sixteen mutants impaired in various splicing-related proteins, many of which had not been recovered in any prior genetic screen or implicated in splicing in plants. The factors are predicted to act at different steps of the spliceosomal cycle, snRNP biogenesis pathway, transcription, and mRNA transport. We have described eleven of the mutants in recent publications. Here we present the final five mutants, which are defective, respectively, in RNA-BINDING PROTEIN 45D (rbp45d), DIGEORGE SYNDROME CRITICAL REGION 14 (dgcr14), CYCLIN-DEPENDENT KINASE G2 (cdkg2), INTERACTS WITH SPT6 (iws1) and CAP BINDING PROTEIN 80 (cbp80). We provide RNA-sequencing data and analyses of differential gene expression and alternative splicing patterns for the cbp80 mutant and for several previously published mutants, including smfa and new alleles of cwc16a, for which such information was not yet available. Sequencing of small RNAs from the cbp80 mutant highlighted the necessity of wild-type CBP80 for processing of microRNA (miRNA) precursors into mature miRNAs. Redundancy tests of paralogs encoding several of the splicing factors revealed their functional non-equivalence in the GFP reporter gene system. We discuss the cumulative findings and their implications for the regulation of pre-mRNA splicing efficiency and alternative splicing in plants. The mutant collection provides a unique resource for further studies on a coherent set of splicing factors and their roles in gene expression, alternative splicing and plant development.

Summary: Alternative splicing is essential for plants, enabling a single gene to produce multiple transcript variants to boost functional diversity and fine‐tune responses to environmental and developmental cues. Arabidopsis thaliana At‐RS31, a plant‐specific splicing factor in the Serine/Arginine‐rich (SR) protein family, responds to light and the Target of Rapamycin (TOR) signalling pathway, yet its downstream targets and regulatory impact remain unknown. To identify At‐RS31 targets, we applied individual‐nucleotide resolution crosslinking and immunoprecipitation (iCLIP) and RNAcompete assays. Transcriptomic analyses of At‐RS31 mutant and overexpressing plants further revealed its effects on alternative splicing. iCLIP identified 4034 At‐RS31 binding sites across 1421 genes, enriched in CU‐rich and CAGA RNA motifs. Comparative iCLIP and RNAcompete data indicate that the arginine/serine (RS) domain of At‐RS31 may influence its binding specificity in planta , underscoring the value of combining in vivo and in vitro approaches. Transcriptomic analysis showed that At‐RS31 modulates diverse splicing events, particularly intron retention and exitron splicing, and influences other splicing modulators, acting as a hierarchical regulator. By regulating stress response genes and genes in both TOR and abscisic acid signalling pathways, At‐RS31 may help integrate these signals, balancing plant growth with environmental adaptability through alternative splicing.

Size control is a fundamental question in biology, showing incremental complexity in plants, whose cells possess a rigid cell wall. The phytohormone auxin is a vital growth regulator with central importance for differential growth control. Our results indicate that auxin-reliant growth programs affect the molecular complexity of xyloglucans, the major type of cell wall hemicellulose in eudicots. Auxin-dependent induction and repression of growth coincide with reduced and enhanced molecular complexity of xyloglucans, respectively. In agreement with a proposed function in growth control, genetic interference with xyloglucan side decorations distinctly modulates auxin-dependent differential growth rates. Our work proposes that auxin-dependent growth programs have a spatially defined effect on xyloglucan’s molecular structure, which in turn affects cell wall mechanics and specifies differential, gravitropic hypocotyl growth.

Plants adjust their architecture to a constantly changing environment, requiring adaptation of differential growth. Despite their importance, molecular switches, which define growth transitions, are largely unknown. Apical hook development in dark grown Arabidopsis thaliana (A. thaliana) seedlings serves as a suitable model for differential growth transition in plants. Here, we show that the phytohormone auxin counteracts the light-induced growth transition during apical hook opening. We, subsequently, identified genes which are inversely regulated by light and auxin. We used in silico analysis of the regulatory elements in this set of genes and subsequently used natural variation in gene expression to uncover correlations between underlying transcription factors and the in silico predicted target genes. This approach uncovered that MADS box transcription factor AGAMOUS-LIKE 8 (AGL8)/FRUITFULL (FUL) modulates apical hook opening. Our data shows that transient FUL expression represses the expression of growth stimulating genes during early phases of apical hook development and therewith guards the transition to growth promotion for apical hook opening. Here, we propose a role for FUL in setting tissue identity, thereby regulating differential growth during apical hook development.