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

Exploring Alternative Splicing in Response to Salinity: A Tissue-Level Comparative Analysis Using Arabidopsis thaliana Public Transcriptomic Data

Increased soil salinity is a major threat to global agriculture and food security, caused mainly by anthropogenic activities and changing climatic cycles. Plants responses to salinity involve multiple regulatory layers, from transcriptome reprogramming to proteomic and metabolomic changes. Alternative splicing (AS) plays a role in coordinating the response to salinity, yet its extent, tissue, and condition specificity, remain poorly understood aspects. In this study, we used 52 publicly available RNA-seq datasets of salinity treatment to identify differential alternative splicing (DAS) events and genes participating in the response to this stimulus. Our findings reveal that either independently or coordinately, AS can regulate up to 20% of the transcriptome detected in Arabidopsis, with treatment intensity being the most determining factor. Moreover, we show that AS regulation was highly tissue-specific, with roots displaying strong AS-mediated stress responses. Furthermore, cross-stress comparisons showed that roots have a core set of AS-regulated genes associated with stress response and development, with functionally distinct sets of genes when comparing salt with other stresses, while also conserving a relevant condition-specific response. We demonstrate the need to integrate AS analysis to better understand plant adaptation mechanisms and highlight the key role of AS in salinity responses, revealing shared AS regulation between salt, heat, and drought responses.

Sporophyte-directed gametogenesis in Arabidopsis

Plants alternate between diploid sporophyte and haploid gametophyte generations1. In mosses, which retain features of ancestral land plants, the gametophyte is dominant and has an independent existence. However, in flowering plants the gametophyte has undergone evolutionary reduction to just a few cells enclosed within the sporophyte. The gametophyte is thought to retain genetic control of its development even after reduction2. Here we show that male gametophyte development in Arabidopsis, long considered to be autonomous, is also under genetic control of the sporophyte via a repressive mechanism that includes large-scale regulation of protein turnover. We identify an Arabidopsis gene SHUKR as an inhibitor of male gametic gene expression. SHUKR is unrelated to proteins of known function and acts sporophytically in meiosis to control gametophyte development by negatively regulating expression of a large set of genes specific to postmeiotic gametogenesis. This control emerged late in evolution as SHUKR homologues are found only in eudicots. We show that SHUKR is rapidly evolving under positive selection, suggesting that variation in control of protein turnover during male gametogenesis has played an important role in evolution within eudicots.

A barley pan-transcriptome reveals layers of genotype-dependent transcriptional complexity

A pan-transcriptome describes the transcriptional and post-transcriptional consequences of genome diversity from multiple individuals within a species. We developed a barley pan-transcriptome using 20 inbred genotypes representing domesticated barley diversity by generating and analyzing short- and long-read RNA-sequencing datasets from multiple tissues. To overcome single reference bias in transcript quantification, we constructed genotype-specific reference transcript datasets (RTDs) and integrated these into a linear pan-genome framework to create a pan-RTD, allowing transcript categorization as core, shell or cloud. Focusing on the core (expressed in all genotypes), we observed significant transcript abundance variation among tissues and between genotypes driven partly by RNA processing, gene copy number, structural rearrangements and conservation of promotor motifs. Network analyses revealed conserved co-expression module::tissue correlations and frequent functional diversification. To complement the pan-transcriptome, we constructed a comprehensive cultivar (cv.) Morex gene-expression atlas and illustrate how these combined datasets can be used to guide biological inquiry.

ZmHSFA2B self-regulatory loop is critical for heat tolerance in maize

The growth and development of maize (Zea mays L.) are significantly impeded by prolonged exposure to high temperatures. Heat stress transcription factors (HSFs) play crucial roles in enabling plants to detect and respond to elevated temperatures. However, the genetic mechanisms underlying the responses of HSFs to heat stress in maize remain unclear. Thus, we aimed to investigate the role of ZmHSFA2B in regulating heat tolerance in maize. Here, we report that ZmHSFA2B has two splicing variants, ZmHSFA2B-I and ZmHSFA2B-II. ZmHSFA2B-I encodes full-length ZmHSFA2B (ZmHSFA2B-I), whereas ZmHSFA2B-II encodes a truncated ZmHSFA2B (ZmHSFA2B-II). Overexpression of ZmHSFA2B-I improved heat tolerance in maize and Arabidopsis thaliana, but it also resulted in growth retardation as a side effect. RNA-sequencing and CUT&Tag analyses identified ZmMBR1 as a putative target of ZmHSFA2B-I. Overexpression of ZmMBR1 also enhanced heat tolerance in Arabidopsis. ZmHSFA2B-II was primarily synthesized in response to heat stress and competitively interacted with ZmHSFA2B-I. This interaction consequently reduced the DNA-binding activities of ZmHSFA2B-I homodimers to the promoter of ZmMBR1. Subsequent investigations indicate that ZmHSFA2B-II limits the transactivation and tempers the function of ZmHSFA2B-I, thereby reducing the adverse effects of excessive ZmHSFA2B-I accumulation. Based on these observations, we propose that the alternative splicing of ZmHSFA2B generates a self-regulatory loop that fine-tunes heat stress response in maize.

At-RS31 orchestrates hierarchical cross-regulation of splicing factors and integrates alternative splicing with TOR-ABA pathways

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. At-RS31, a plant-specific splicing factor in the Serine/Arginine (SR)-rich protein family, responds to light and the Target of Rapamycin (TOR) signaling 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 4,034 At-RS31 binding sites across 1,421 genes, enriched in CU-rich and CAGA RNA motifs. Comparative iCLIP and RNAcompete data indicate that the 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 (ABA) signaling pathways, At-RS31 may help integrate these signals, balancing plant growth with environmental adaptability through alternative splicing.

The nuclear exosome subunit HEN2 acts independently of the core exosome to assist transcription in Arabidopsis

Regulation of gene expression is at the frontier of plant responses to various external stimuli including stress. RNA polymerase-based transcription and post-transcriptional degradation of RNA play vital roles in this regulation. Here, we show that HUA ENHANCER 2 (HEN2), a co-factor of the nuclear exosome complex, influences RNAPII transcription elongation in Arabidopsis (Arabidopsis thaliana) under cold conditions. Our results demonstrate that a hen2 mutant is cold sensitive and undergoes substantial transcriptional changes compared to wild type when exposed to cold conditions. We found an accumulation of 5' fragments from a subset of genes (including C-repeat binding factors 1-3 [CBF1-3]) that do not carry over to their 3' ends. In fact, hen2 mutants have lower levels of full-length mRNA for a subset of genes. This distinct 5'-end accumulation and 3'-end depletion was not observed in other NEXT complex members or core exosome mutants, highlighting HEN2's distinctive role. We further used RNAPII-associated nascent RNA to confirm that the transcriptional phenotype is a result of lower active transcription specifically at the 3' end of these genes in a hen2 mutant. Taken together, our data point to the unique role of HEN2 in maintaining RNAPII transcription dynamics especially highlighted under cold stress.

The Arabidopsis splicing factor PORCUPINE/SmE1 orchestrates temperature-dependent root development via auxin homeostasis maintenance

Appropriate abiotic stress response is pivotal for plant survival and makes use of multiple signaling molecules and phytohormones to achieve specific and fast molecular adjustments. A multitude of studies has highlighted the role of alternative splicing in response to abiotic stress, including temperature, emphasizing the role of transcriptional regulation for stress response. Here we investigated the role of the core-splicing factor PORCUPINE (PCP) on temperature-dependent root development. We used marker lines and transcriptomic analyses to study the expression profiles of meristematic regulators and mitotic markers, and chemical treatments, as well as root hormone profiling to assess the effect of auxin signaling. The loss of PCP significantly alters RAM architecture in a temperature-dependent manner. Our results indicate that PCP modulates the expression of central meristematic regulators and is required to maintain appropriate levels of auxin in the RAM. We conclude that alternative pre-mRNA splicing is sensitive to moderate temperature fluctuations and contributes to root meristem maintenance, possibly through the regulation of phytohormone homeostasis and meristematic activity.

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.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Transcriptional and hormonal profiling uncovers the interactions between plant developmental stages and RNA virus infection

Arabidopsis thaliana is more susceptible to certain viruses during its later developmental stages. The differential responses and the mechanisms behind this development-dependent susceptibility to infection are still not fully understood. Here we explored the outcome of a viral infection at different host developmental stages by studying the response of A. thaliana to infection with turnip mosaic virus at three developmental stages: juvenile vegetative, bolting, and mature flowering plants. We found that infected plants at later stages downregulate cell wall biosynthetic genes and that this downregulation may be one factor facilitating viral spread and systemic infection. We also found that, despite being more susceptible to infection, infected mature flowering plants were more fertile (i.e. produce more viable seeds) than juvenile vegetative and bolting infected plants; that is, plants infected at the reproductive stage have greater fitness than plants infected at earlier developmental stages. Moreover, treatment of mature plants with salicylic acid increased resistance to infection at the cost of significantly reducing fertility. Together, these observations support a negative trade-off between viral susceptibility and plant fertility. Our findings point towards a development-dependent tolerance to infection.

Coordination of shoot apical meristem shape and identity by APETALA2 during floral transition in Arabidopsis

Plants flower in response to environmental signals. These signals change the shape and developmental identity of the shoot apical meristem (SAM), causing it to form flowers and inflorescences. We show that the increases in SAM width and height during floral transition correlate with changes in size of the central zone (CZ), defined by CLAVATA3 expression, and involve a transient increase in the height of the organizing center (OC), defined by WUSCHEL expression. The APETALA2 (AP2) transcription factor is required for the rapid increases in SAM height and width, by maintaining the width of the OC and increasing the height and width of the CZ. AP2 expression is repressed in the SAM at the end of floral transition, and extending the duration of its expression increases SAM width. Transcriptional repression by SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) represents one of the mechanisms reducing AP2 expression during floral transition. Moreover, AP2 represses SOC1 transcription, and we find that reciprocal repression of SOC1 and AP2 contributes to synchronizing precise changes in meristem shape with floral transition.

SMRT sequencing of a full-length transcriptome reveals cold induced alternative splicing in Vitis amurensis root

Alternative splicing enhances diversity at the transcriptional and protein levels that widely involved in plant response to biotic and abiotic stresses. V. amurensis is an extremely cold-tolerant wild grape variety, however, studies on alternative splicing (AS) in amur grape at low temperatures are currently poorly understood. In this study, we analyzed full-length transcriptome and RNA seq data at 0, 2, and 24 h after cold stress in V. amurensis roots. Following quality control and correction, 221,170 high-quality full-length non-concatemer (FLNC) reads were identified. A total of 16,181 loci and 30,733 isoforms were identified. These included 22,868 novel isoforms from annotated genes and 2815 isoforms from 2389 novel genes. Among the distinguished novel isoforms, 673 Long non-coding RNAs (LncRNAs) and 18,164 novel isoforms open reading frame (ORF) region were found. A total of 2958 genes produced 8797 AS events, of which 189 genes were involved in the low-temperature response. Twelve transcription factors show AS during cold treatment and VaMYB108 was selected for initial exploration. Two transcripts, Chr05.63.1 (VaMYB108short) and Chr05.63.2 (VaMYB108normal) of VaMYB108, display up-regulated expression after cold treatment in amur grape roots and are both localized in the nucleus. Only VaMYB108normal exhibits transcriptional activation activity. Overexpression of either VaMYB108short or VaMYB108normal in grape roots leads to increased expression of the other transcript and both increased chilling resistance of amur grape roots. The results improve and supplement the genome annotations and provide insights for further investigation into AS mechanisms during cold stress in V. amurensis.

AtSNU13 modulates pre-mRNA splicing of RBOHD and ALD1 to regulate plant immunity

Pre-mRNA splicing is a significant step for post-transcriptional modifications and functions in a wide range of physiological processes in plants. Human NHP2L binds to U4 snRNA during spliceosome assembly; it is involved in RNA splicing and mediates the development of human tumors. However, no ortholog has yet been identified in plants. Therefore, we report At4g12600 encoding the ortholog NHP2L protein, and AtSNU13 associates with the component of the spliceosome complex; the atsnu13 mutant showed compromised resistance in disease resistance, indicating that AtSNU13 is a positive regulator of plant immunity. Compared to wild-type plants, the atsnu13 mutation resulted in altered splicing patterns for defense-related genes and decreased expression of defense-related genes, such as RBOHD and ALD1. Further investigation shows that AtSNU13 promotes the interaction between U4/U6.U5 tri-snRNP-specific 27 K and the motif in target mRNAs to regulate the RNA splicing. Our study highlights the role of AtSNU13 in regulating plant immunity by affecting the pre-mRNA splicing of defense-related genes.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Identification of a highly drought-resistant pp7l hda6 mutant

Plants have developed efficient strategies to counteract drought stress, including stomata closure, significant changes in nuclear gene expression, and epigenetic mechanisms. Previously, we identified Arabidopsis thaliana PROTEIN PHOSPHATASE7-LIKE (PP7L) as an extrachloroplastic protein that promotes chloroplast development. In addition, it was shown that PP7L is involved in high light and salt tolerance. Here, we demonstrate that the pp7l mutant can withstand prolonged periods of drought stress. Interestingly, despite impaired growth under standard growth conditions, photosynthetic efficiency recovers in pp7l mutant plants experiencing drought conditions. To assess the (post)transcriptional changes occurring in the pp7l mutant under different durations of drought exposure, we used an RNA-sequencing technique that allows the simultaneous detection of organellar and nuclear transcripts. Compared with the previously reported drought-responsive changes in the wild type, the drought-responsive changes in organellar and nuclear transcripts detected in the pp7l mutant were negligible. Our analysis of the data generated in this study and review and analysis of previous literature motivated us to create a pp7l hda6 (histone deacetylase 6) mutant, which exhibits remarkable drought resistance. Notably, the growth penalty associated with pp7l was alleviated in the double mutant, ruling out a dwarf effect on the drought-tolerant trait of this genotype. Future studies may consider that multiple loci and factors are involved in stress resistance and explore combinations of these factors to create even more resilient plants.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

Mapping protein-RNA binding in plants with individual-nucleotide-resolution UV cross-linking and immunoprecipitation (plant iCLIP2)

Despite crucial roles of RNA-binding proteins (RBPs) in plant physiology and development, methods for determining their transcriptome-wide binding landscape are less developed than those used in other model organisms. Cross-linking and immunoprecipitation (CLIP) methods (based on UV-mediated generation of covalent bonds between RNAs and cognate RBPs in vivo, purification of the cross-linked complexes and identification of the co-purified RNAs by high-throughput sequencing) have been applied mainly in mammalian cells growing in monolayers or in translucent tissue. We have developed plant iCLIP2, an efficient protocol for performing individual-nucleotide-resolution CLIP (iCLIP) in plants, tailored to overcome the experimental hurdles posed by plant tissue. We optimized the UV dosage to efficiently cross-link RNA and proteins in plants and expressed epitope-tagged RBPs under the control of their native promoters in loss-of-function mutants. We select epitopes for which nanobodies are available, allowing stringent conditions for immunopurification of the RNA-protein complexes to be established. To overcome the inherently high RNase content of plant cells, RNase inhibitors are added and the limited RNA fragmentation step is modified. We combine the optimized isolation of RBP-bound RNAs with iCLIP2, a streamlined protocol that greatly enhances the efficiency of library preparation for high-throughput sequencing. Plant researchers with experience in molecular biology and handling of RNA can complete this iCLIP2 protocol in ~5 d. Finally, we describe a bioinformatics workflow to determine targets of Arabidopsis RBPs from iCLIP data, covering all steps from downloading sequencing reads to identifying cross-linking events ( https://github.com/malewins/Plant-iCLIPseq ), and present the R/Bioconductor package BindingSiteFinder to extract reproducible binding sites ( https://bioconductor.org/packages/release/bioc/html/BindingSiteFinder.html ).

Arabidopsis thaliana GLYCINE RICH RNA-BINDING PROTEIN 7 interaction with its iCLIP target LHCB1.1 correlates with changes in RNA stability and circadian oscillation

The importance of RNA-binding proteins (RBPs) for plant responses to environmental stimuli and development is well documented. Insights into the portfolio of RNAs they recognize, however, clearly lack behind the understanding gathered in non-plant model organisms. Here, we characterize binding of the circadian clock-regulated Arabidopsis thaliana GLYCINE-RICH RNA-BINDING PROTEIN 7 (AtGRP7) to its target transcripts. We identified novel RNA targets from individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) data using an improved bioinformatics pipeline that will be broadly applicable to plant RBP iCLIP data. 2705 transcripts with binding sites were identified in plants expressing AtGRP7-GFP that were not recovered in plants expressing an RNA-binding dead variant or GFP alone. A conserved RNA motif enriched in uridine residues was identified at the AtGRP7 binding sites. NMR titrations confirmed the preference of AtGRP7 for RNAs with a central U-rich motif. Among the bound RNAs, circadian clock-regulated transcripts were overrepresented. Peak abundance of the LHCB1.1 transcript encoding a chlorophyll-binding protein was reduced in plants overexpressing AtGRP7 whereas it was elevated in atgrp7 mutants, indicating that LHCB1.1 was regulated by AtGRP7 in a dose-dependent manner. In plants overexpressing AtGRP7, the LHCB1.1 half-life was shorter compared to wild-type plants whereas in atgrp7 mutant plants, the half-life was significantly longer. Thus, AtGRP7 modulates circadian oscillations of its in vivo binding target LHCB1.1 by affecting RNA stability.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Environmental regulation of male fertility is mediated through Arabidopsis transcription factors bHLH89, 91, and 10

Formation of functional pollen and successful fertilization rely on the spatial and temporal regulation of anther and pollen development. This process responds to environmental cues to maintain optimal fertility despite climatic changes. Arabidopsis transcription factors basic helix-loop-helix (bHLH) 10, 89, and 91 were previously thought to be functionally redundant in their control of male reproductive development, however here we show that they play distinct roles in the integration of light signals to maintain pollen development under different environmental conditions. Combinations of the double and triple bHLH10,89,91 mutants were analysed under normal (200 μmol m-2 s-1) and low (50 μmol m-2 s-1) light conditions to determine the impact on fertility. Transcriptomic analysis of a new conditionally sterile bhlh89,91 double mutant shows differential regulation of genes related to sexual reproduction, hormone signal transduction, and lipid storage and metabolism under low light. Here we have shown that bHLH89 and bHLH91 play a role in regulating fertility in response to light, suggesting that they function in mitigating environmental variation to ensure fertility is maintained under environmental stress.

Identification of epigenetically regulated genes involved in plant-virus interaction and their role in virus-triggered induced resistance

BACKGROUND

Plant responses to a wide range of stresses are known to be regulated by epigenetic mechanisms. Pathogen-related investigations, particularly against RNA viruses, are however scarce. It has been demonstrated that Arabidopsis thaliana plants defective in some members of the RNA-directed DNA methylation (RdDM) or histone modification pathways presented differential susceptibility to the turnip mosaic virus. In order to identify genes directly targeted by the RdDM-related RNA Polymerase V (POLV) complex and the histone demethylase protein JUMONJI14 (JMJ14) during infection, the transcriptomes of infected mutant and control plants were obtained and integrated with available chromatin occupancy data for various epigenetic proteins and marks.

RESULTS

A comprehensive list of virus-responsive gene candidates to be regulated by the two proteins was obtained. Twelve genes were selected for further characterization, confirming their dynamic regulation during the course of infection. Several epigenetic marks on their promoter sequences were found using in silico data, raising confidence that the identified genes are actually regulated by epigenetic mechanisms. The altered expression of six of these genes in mutants of the methyltransferase gene CURLY LEAF and the histone deacetylase gene HISTONE DEACETYLASE 19 suggests that some virus-responsive genes may be regulated by multiple coordinated epigenetic complexes. A temporally separated multiple plant virus infection experiment in which plants were transiently infected with one virus and then infected by a second one was designed to investigate the possible roles of the identified POLV- and JMJ14-regulated genes in wild-type (WT) plants. Plants that had previously been stimulated with viruses were found to be more resistant to subsequent virus challenge than control plants. Several POLV- and JMJ14-regulated genes were found to be regulated in virus induced resistance in WT plants, with some of them poisoned to be expressed in early infection stages.

CONCLUSIONS

A set of confident candidate genes directly regulated by the POLV and JMJ14 proteins during virus infection was identified, with indications that some of them may be regulated by multiple epigenetic modules. A subset of these genes may also play a role in the tolerance of WT plants to repeated, intermittent virus infections.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

BCL7A and BCL7B potentiate SWI/SNF-complex-mediated chromatin accessibility to regulate gene expression and vegetative phase transition in plants

Switch defective/sucrose non-fermentable (SWI/SNF) chromatin remodeling complexes are multi-subunit machineries that establish and maintain chromatin accessibility and gene expression by regulating chromatin structure. However, how the remodeling activities of SWI/SNF complexes are regulated in eukaryotes remains elusive. B-cell lymphoma/leukemia protein 7 A/B/C (BCL7A/B/C) have been reported as subunits of SWI/SNF complexes for decades in animals and recently in plants; however, the role of BCL7 subunits in SWI/SNF function remains undefined. Here, we identify a unique role for plant BCL7A and BCL7B homologous subunits in potentiating the genome-wide chromatin remodeling activities of SWI/SNF complexes in plants. BCL7A/B require the catalytic ATPase BRAHMA (BRM) to assemble with the signature subunits of the BRM-Associated SWI/SNF complexes (BAS) and for genomic binding at a subset of target genes. Loss of BCL7A and BCL7B diminishes BAS-mediated genome-wide chromatin accessibility without changing the stability and genomic targeting of the BAS complex, highlighting the specialized role of BCL7A/B in regulating remodeling activity. We further show that BCL7A/B fine-tune the remodeling activity of BAS complexes to generate accessible chromatin at the juvenility resetting region (JRR) of the microRNAs MIR156A/C for plant juvenile identity maintenance. In summary, our work uncovers the function of previously elusive SWI/SNF subunits in multicellular eukaryotes and provides insights into the mechanisms whereby plants memorize the juvenile identity through SWI/SNF-mediated control of chromatin accessibility.

REVEILLE2 thermosensitive splicing: a molecular basis for the integration of nocturnal temperature information by the Arabidopsis circadian clock

Cold stress is one of the major environmental factors that limit growth and yield of plants. However, it is still not fully understood how plants account for daily temperature fluctuations, nor how these temperature changes are integrated with other regulatory systems such as the circadian clock. We demonstrate that REVEILLE2 undergoes alternative splicing after chilling that increases accumulation of a transcript isoform encoding a MYB-like transcription factor. We explore the biological function of REVEILLE2 in Arabidopsis thaliana using a combination of molecular genetics, transcriptomics, and physiology. Disruption of REVEILLE2 alternative splicing alters regulatory gene expression, impairs circadian timing, and improves photosynthetic capacity. Changes in nuclear gene expression are particularly apparent in the initial hours following chilling, with chloroplast gene expression subsequently upregulated. The response of REVEILLE2 to chilling extends our understanding of plants immediate response to cooling. We propose that the circadian component REVEILLE2 restricts plants responses to nocturnal reductions in temperature, thereby enabling appropriate responses to daily environmental changes.

A GFP splicing reporter in a coilin mutant background reveals links between alternative splicing, siRNAs, and coilin function in Arabidopsis thaliana

Coilin is a scaffold protein essential for the structure of Cajal bodies, which are nucleolar-associated, nonmembranous organelles that coordinate the assembly of nuclear ribonucleoproteins (RNPs) including spliceosomal snRNPs. To study coilin function in plants, we conducted a genetic suppressor screen using a coilin (coi1) mutant in Arabidopsis thaliana and performed an immunoprecipitation-mass spectrometry analysis on coilin protein. The coi1 mutations modify alternative splicing of a GFP reporter gene, resulting in a hyper-GFP phenotype in young coi1 seedlings relative to the intermediate wild-type level. As shown here, this hyper-GFP phenotype is extinguished in older coi1 seedlings by posttranscriptional gene silencing triggered by siRNAs derived from aberrant splice variants of GFP pre-mRNA. In the coi1 suppressor screen, we identified suppressor mutations in WRAP53, a putative coilin-interacting protein; SMU2, a predicted splicing factor; and ZCH1, an incompletely characterized zinc finger protein. These suppressor mutations return the hyper-GFP fluorescence of young coi1 seedlings to the intermediate wild-type level. Additionally, coi1 zch1 mutants display more extensive GFP silencing and elevated levels of GFP siRNAs, suggesting the involvement of wild-type ZCH1 in siRNA biogenesis or stability. The immunoprecipitation-mass spectrometry analysis reinforced the roles of coilin in pre-mRNA splicing, nucleolar chromatin structure, and rRNA processing. The participation of coilin in these processes, at least some of which incorporate small RNAs, supports the hypothesis that coilin provides a chaperone for small RNA trafficking. Our study demonstrates the usefulness of the GFP splicing reporter for investigating alternative splicing, ribosome biogenesis, and siRNA-mediated silencing in the context of coilin function.

Response of the organellar and nuclear (post)transcriptomes of Arabidopsis to drought

Plants have evolved sophisticated mechanisms to cope with drought, which involve massive changes in nuclear gene expression. However, little is known about the roles of post-transcriptional processing of nuclear or organellar transcripts and how meaningful these changes are. To address these issues, we used RNA-sequencing after ribosomal RNA depletion to monitor (post)transcriptional changes during different times of drought exposure in Arabidopsis Col-0. Concerning the changes detected in the organellar transcriptomes, chloroplast transcript levels were globally reduced, editing efficiency dropped, but splicing was not affected. Mitochondrial transcripts were slightly elevated, while editing and splicing were unchanged. Conversely, alternative splicing (AS) affected nearly 1,500 genes (9% of expressed nuclear genes). Of these, 42% were regulated solely at the level of AS, representing transcripts that would have gone unnoticed in a microarray-based approach. Moreover, we identified 927 isoform switching events. We provide a table of the most interesting candidates, and as proof of principle, increased drought tolerance of the carbonic anhydrase ca1 and ca2 mutants is shown. In addition, altering the relative contributions of the spliced isoforms could increase drought resistance. For example, our data suggest that the accumulation of a nonfunctional FLM (FLOWERING LOCUS M) isoform and not the ratio of FLM-ß and -δ isoforms may be responsible for the phenotype of early flowering under long-day drought conditions. In sum, our data show that AS enhances proteome diversity to counteract drought stress and represent a valuable resource that will facilitate the development of new strategies to improve plant performance under drought.

Nuclear OsFKBP20-1b maintains SR34 stability and promotes the splicing of retained introns upon ABA exposure in rice

Alternative splicing (AS) is a critical means by which plants respond to changes in the environment, but few splicing factors contributing to AS have been reported and functionally characterized in rice (Oryza sativa L.). Here, we explored the function and molecular mechanism of the spliceosome-associated protein OsFKBP20-1b during AS. We determined the AS landscape of wild-type and osfkbp20-1b knockout plants upon abscisic acid (ABA) treatment by transcriptome deep sequencing. To capture the dynamics of translating intron-containing mRNAs, we blocked transcription with cordycepin and performed polysome profiling. We also analyzed whether OsFKBP20-1b and the splicing factors OsSR34 and OsSR45 function together in AS using protoplast transfection assays. We show that OsFKBP20-1b interacts with OsSR34 and regulates its stability, suggesting a role as a chaperone-like protein in the spliceosome. OsFKBP20-1b facilitates the splicing of mRNAs with retained introns after ABA treatment; some of these mRNAs are translatable and encode functional transcriptional regulators of stress-responsive genes. In addition, interacting proteins, OsSR34 and OsSR45, regulate the splicing of the same retained introns as OsFKBP20-1b after ABA treatment. Our findings reveal that spliceosome-associated immunophilin functions in alternative RNA splicing in rice by positively regulating the splicing of retained introns to limit ABA response.

Arabidopsis SRPKII family proteins regulate flowering via phosphorylation of SR proteins and effects on gene expression and alternative splicing

Alternative splicing of pre-mRNAs is crucial for plant growth and development. Serine/arginine-rich (SR) proteins are a conserved family of RNA-binding proteins that are critical for both constitutive and alternative splicing. However, how phosphorylation of SR proteins regulates gene transcription and alternative splicing during plant development is poorly understood. We found that the Arabidopsis thaliana L. SR protein-specific kinase II family proteins (SRPKIIs) play an important role in plant development, including flowering. SRPKIIs regulate the phosphorylation status of a subset of specific SR proteins, including SR45 and SC35, which subsequently mediates their subcellular localization. A phospho-dead SR45 mutant inhibits the assembly of the apoptosis-and splicing-associated protein complex and thereby upregulates the expression of FLOWERING LOCUS C (FLC) via epigenetic modification. The splicing efficiency of FLC introns was significantly increased in the shoot apex of the srpkii mutant. Transcriptomic analysis revealed that SRPKIIs regulate the alternative splicing of c. 400 genes, which largely overlap with those regulated by SR45 and SC35-SCL family proteins. In summary, we found that Arabidopsis SRPKIIs specifically affect the phosphorylation status of a subset SR proteins and regulate the expression and alternative splicing of FLC to control flowering time.

Synergistic action of the Arabidopsis spliceosome components PRP39a and SmD1b in promoting posttranscriptional transgene silencing

Besides regulating splicing, the conserved spliceosome component SmD1 (Small nuclear ribonucleoprotein D1)b promotes posttranscriptional silencing of sense transgenes (S-PTGS [post-transcriptional genesilencing]). Here, we show that the conserved spliceosome component PRP39 (Pre-mRNA-processing factor 39)a also plays a role in S-PTGS in Arabidopsis thaliana. However, PRP39a and SmD1b actions appear distinct in both splicing and S-PTGS. Indeed, RNAseq-based analysis of expression level and alternative splicing in prp39a and smd1b mutants identified different sets of deregulated transcripts and noncoding RNAs. Moreover, double mutant analyses involving prp39a or smd1b and RNA quality control (RQC) mutants revealed distinct genetic interactions for SmD1b and PRP39a with nuclear RQC machineries, suggesting nonredundant roles in the RQC/PTGS interplay. Supporting this hypothesis, a prp39a smd1b double mutant exhibited enhanced suppression of S-PTGS compared to the single mutants. Because the prp39a and smd1b mutants (i) showed no major changes in the expression of PTGS or RQC components or in small RNA production and (ii) do not alter PTGS triggered by inverted-repeat transgenes directly producing dsRNA (IR-PTGS), PRP39a, and SmD1b appear to synergistically promote a step specific to S-PTGS. We propose that, independently from their specific roles in splicing, PRP39a and SmD1b limit 3'-to-5' and/or 5'-to-3' degradation of transgene-derived aberrant RNAs in the nucleus, thus favoring the export of aberrant RNAs to the cytoplasm where their conversion into double-stranded RNA (dsRNA) initiates S-PTGS.

Do not panic: An intron-centric guide to alternative splicing

This review is an attempt to establish concepts of splicing and alternative splicing giving proper relevance to introns, the key actors in this mechanism. It might also work as a guide for those who found their favorite gene undergoes alternative splicing and could benefit from gaining a theoretical framework to understand the possible impacts of this process. This is not a thorough review of all the work in the field, but rather a critical review of some of the most relevant work done to understand the underlying mechanisms of splicing and the key questions that remain unanswered such as: What is the physiological relevance of alternative splicing? What are the functions of the different outcomes? To what extent do different alternative splicing types contribute to the proteome? Intron retention is the most frequent alternative splicing event in plants and, although scientifically neglected, it is also common in animals. This is a heterogeneous type of alternative splicing that includes different sub-types with features that have distinctive consequences in the resulting transcripts. Remarkably, intron retention can be a dead end for a transcript, but it could also be a stable intermediate whose processing is resumed upon a particular signal or change in the cell status. New sequencing technologies combined with the study of intron lariats in different conditions might help to answer key questions and could help us to understand the actual relevance of introns in gene expression regulation.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

The flowering time regulator FLK controls pathogen defense in Arabidopsis thaliana

Plant disease resistance is a complex process that is maintained in an intricate balance with development. Increasing evidence indicates the importance of posttranscriptional regulation of plant defense by RNA binding proteins. In a genetic screen for suppressors of Arabidopsis (Arabidopsis thaliana) accelerated cell death 6-1 (acd6-1), a small constitutive defense mutant whose defense level is grossly in a reverse proportion to plant size, we identified an allele of the canonical flowering regulatory gene FLOWERING LOCUS K HOMOLOGY DOMAIN (FLK) encoding a putative protein with triple K homology (KH) repeats. The KH repeat is an ancient RNA binding motif found in proteins from diverse organisms. The relevance of KH-domain proteins in pathogen resistance is largely unexplored. In addition to late flowering, the flk mutants exhibited decreased resistance to the bacterial pathogen Pseudomonas syringae and increased resistance to the necrotrophic fungal pathogen Botrytis cinerea. We further found that the flk mutations compromised basal defense and defense signaling mediated by salicylic acid (SA). Mutant analysis revealed complex genetic interactions between FLK and several major SA pathway genes. RNA-seq data showed that FLK regulates expression abundance of some major defense- and development-related genes as well as alternative splicing of a number of genes. Among the genes affected by FLK is ACD6, whose transcripts had increased intron retentions influenced by the flk mutations. Thus, this study provides mechanistic support for flk suppression of acd6-1 and establishes that FLK is a multifunctional gene involved in regulating pathogen defense and development of plants.

Transcription Factor Dynamics in Cross-Regulation of Plant Hormone Signaling Pathways

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

Mutation of Arabidopsis SME1 and Sm core assembly improves oxidative stress resilience

Alternative splicing is a key posttranscriptional gene regulatory process, acting in diverse adaptive and basal plant processes. Splicing of precursor-messenger RNA (pre-mRNA) is catalyzed by a dynamic ribonucleoprotein complex, designated the spliceosome. In a suppressor screen, we identified a nonsense mutation in the Smith (Sm) antigen protein SME1 to alleviate photorespiratory H₂O₂-dependent cell death in catalase deficient plants. Similar attenuation of cell death was observed upon chemical inhibition of the spliceosome, suggesting pre-mRNA splicing inhibition to be responsible for the observed cell death alleviation. Furthermore, the sme1-2 mutants showed increased tolerance to the reactive oxygen species inducing herbicide methyl viologen. Both an mRNA-seq and shotgun proteomic analysis in sme1-2 mutants displayed a constitutive molecular stress response, together with extensive alterations in pre-mRNA splicing of transcripts encoding metabolic enzymes and RNA binding proteins, even under unstressed conditions. Using SME1 as a bait to identify protein interactors, we provide experimental evidence for almost 50 homologs of the mammalian spliceosome-associated protein to reside in the Arabidopsis thaliana spliceosome complexes and propose roles in pre-mRNA splicing for four uncharacterized plant proteins. Furthermore, as for sme1-2, a mutant in the Sm core assembly protein ICLN resulted in a decreased sensitivity to methyl viologen. Taken together, these data show that both a perturbed Sm core composition and assembly results in the activation of a defense response and in enhanced resilience to oxidative stress.

Plant Intron-Splicing Efficiency Database (PISE): exploring splicing of ∼1,650,000 introns in Arabidopsis, maize, rice, and soybean from ∼57,000 public RNA-seq libraries

Intron retention is the most common alternative splicing event in plants and plays a crucial role in the responses of plants to environmental signals. Despite a large number of RNA-seq libraries from different treatments and genetic mutants stored in public domains, a resource for querying the intron-splicing ratio of individual intron is still required. Here, we established the first-ever large-scale splicing efficiency database in any organism. Our database includes over 57,000 plant public RNA-seq libraries, comprising 25,283 from Arabidopsis, 17,789 from maize, 10,710 from rice, and 3,974 from soybean, and covers a total of 1.6 million introns in these four species. In addition, we manually curated and annotated all the mutant- and treatment-related libraries as well as their matched controls included in our library collection, and added graphics to display intron-splicing efficiency across various tissues, developmental stages, and stress-related conditions. The result is a large collection of 3,313 treatment conditions and 3,594 genetic mutants for discovering differentially regulated splicing efficiency. Our online database can be accessed at https://plantintron.com/ .

A complex role of Arabidopsis CDKD;3 in meiotic progression and cytokinesis

Meiosis is a specialized cell division that halves the number of chromosomes in two consecutive rounds of chromosome segregation. In angiosperm plants is meiosis followed by mitotic divisions to form rudimentary haploid gametophytes. In Arabidopsis, termination of meiosis and transition to gametophytic development are governed by TDM1 and SMG7 that mediate inhibition of translation. Mutants deficient in this mechanism do not form tetrads but instead undergo multiple cycles of aberrant nuclear divisions that are likely caused by the failure to downregulate cyclin dependent kinases during meiotic exit. A suppressor screen to identify genes that contribute to meiotic exit uncovered a mutation in cyclin-dependent kinase D;3 (CDKD;3) that alleviates meiotic defects in smg7 deficient plants. The CDKD;3 deficiency prevents aberrant meiotic divisions observed in smg7 mutants or delays their onset after initiation of cytokinesis, which permits formation of functional microspores. Although CDKD;3 acts as an activator of cyclin-dependent kinase A;1 (CDKA;1), the main cyclin dependent kinase that regulates meiosis, cdkd;3 mutation appears to promote meiotic exit independently of CDKA;1. Furthermore, analysis of CDKD;3 interactome revealed enrichment for proteins implicated in cytokinesis, suggesting a more complex function of CDKD;3 in cell cycle regulation.

The Integration of Data from Different Long-Read Sequencing Platforms Enhances Proteoform Characterization in Arabidopsis

The increasing availability of massive omics data requires improving the quality of reference databases and their annotations. The combination of full-length isoform sequencing (Iso-Seq) with short-read transcriptomics and proteomics has been successfully used for increasing proteoform characterization, which is a main ongoing goal in biology. However, the potential of including Oxford Nanopore Technologies Direct RNA Sequencing (ONT-DRS) data has not been explored. In this paper, we analyzed the impact of combining Iso-Seq- and ONT-DRS-derived data on the identification of proteoforms in Arabidopsis MS proteomics data. To this end, we selected a proteomics dataset corresponding to senescent leaves and we performed protein searches using three different protein databases: AtRTD2 and AtRTD3, built from the homonymous transcriptomes, regarded as the most complete and up-to-date available for the species; and a custom hybrid database combining AtRTD3 with publicly available ONT-DRS transcriptomics data generated from Arabidopsis leaves. Our results show that the inclusion and combination of long-read sequencing data from Iso-Seq and ONT-DRS into a proteogenomic workflow enhances proteoform characterization and discovery in bottom-up proteomics studies. This represents a great opportunity to further investigate biological systems at an unprecedented scale, although it brings challenges to current protein searching algorithms.

Transcriptomic and splicing changes underlying tomato responses to combined water and nutrient stress

Tomato is a horticultural crop of high economic and nutritional value. Suboptimal environmental conditions, such as limited water and nutrient availability, cause severe yield reductions. Thus, selection of genotypes requiring lower inputs is a goal for the tomato breeding sector. We screened 10 tomato varieties exposed to water deficit, low nitrate or a combination of both. Biometric, physiological and molecular analyses revealed different stress responses among genotypes, identifying T270 as severely affected, and T250 as tolerant to the stresses applied. Investigation of transcriptome changes caused by combined stress in roots and leaves of these two genotypes yielded a low number of differentially expressed genes (DEGs) in T250 compared to T270, suggesting that T250 tailors changes in gene expression to efficiently respond to combined stress. By contrast, the susceptible tomato activated approximately one thousand and two thousand genes in leaves and roots respectively, indicating a more generalized stress response in this genotype. In particular, developmental and stress-related genes were differentially expressed, such as hormone responsive factors and transcription factors. Analysis of differential alternative splicing (DAS) events showed that combined stress greatly affects the splicing landscape in both genotypes, highlighting the important role of AS in stress response mechanisms. In particular, several stress and growth-related genes as well as transcription and splicing factors were differentially spliced in both tissues. Taken together, these results reveal important insights into the transcriptional and post-transcriptional mechanisms regulating tomato adaptation to growth under reduced water and nitrogen inputs.

Comparative Analysis of Alternative Splicing in Two Contrasting Apple Cultivars Defense against Alternaria alternata Apple Pathotype Infection

Alternaria blotch disease, caused by the Alternaria alternata apple pathotype (A. alternata AP), is one of the most serious fungal diseases in apples. Alternative splicing (AS), one of the pivotal post-transcriptional regulatory mechanisms, plays essential roles in various disease resistance responses. Here, we performed RNA-Seq for two apple cultivars (resistant cultivar 'Jonathan' (J) and susceptible cultivar 'Starking Delicious' (SD)) infected by A. alternata AP to further investigate their AS divergence. In total, 1454, 1780, 1367 and 1698 specifically regulated differential alternative splicing (DAS) events were detected in J36, J72, SD36 and SD72 groups, respectively. Retained intron (RI) was the dominant AS pattern. Conformably, 642, 764, 585 and 742 uniquely regulated differentially spliced genes (DSGs) were found during A. alternata AP infection. Comparative analysis of AS genes in differential splicing and expression levels suggested that only a small proportion of DSGs overlapped with differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis demonstrated that the DSGs were significantly enriched at multiple levels of gene expression regulation. Briefly, the specific AS was triggered in apple defense against A. alternata AP. Therefore, this study facilitates our understanding on the roles of AS regulation in response to A. alternata AP infection in apples.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

Genome-scale analysis of Arabidopsis splicing-related protein kinase families reveals roles in abiotic stress adaptation

Nearly 60 - 80 % of intron-containing plant genes undergo alternative splicing in response to either stress or plant developmental cues. RNA splicing is performed by a large ribonucleoprotein complex called the spliceosome in conjunction with associated subunits such as serine arginine (SR) proteins, all of which undergo extensive phosphorylation. In plants, there are three main protein kinase families suggested to phosphorylate core spliceosome subunits and related splicing factors based on orthology to human splicing-related kinases: the SERINE/ARGININE PROTEIN KINASES (SRPK), ARABIDOPSIS FUS3 COMPLEMENT (AFC), and Pre-mRNA PROCESSING FACTOR 4 (PRP4K) protein kinases. To better define the conservation and role(s) of these kinases in plants, we performed a genome-scale analysis of the three families across photosynthetic eukaryotes, followed by extensive transcriptomic and bioinformatic analysis of all Arabidopsis thaliana SRPK, AFC, and PRP4K protein kinases to elucidate their biological functions. Unexpectedly, this revealed the existence of SRPK and AFC phylogenetic groups with distinct promoter elements and patterns of transcriptional response to abiotic stress, while PRP4Ks possess no phylogenetic sub-divisions, suggestive of functional redundancy. We also reveal splicing-related kinase families are both diel and photoperiod regulated, implicating different orthologs as discrete time-of-day RNA splicing regulators. This foundational work establishes a number of new hypotheses regarding how reversible spliceosome phosphorylation contributes to both diel plant cell regulation and abiotic stress adaptation in plants.

Differential alternative splicing genes and isoform co-expression networks of Brassica napus under multiple abiotic stresses

Alternative splicing (AS) is an important regulatory process that affects plant development and stress responses by greatly increasing the complexity of transcriptome and proteome. To understand how the AS landscape of B. napus changes in response to abiotic stresses, we investigated 26 RNA-seq libraries, including control and treatments with cold, dehydration, salt, and abscisic acid (ABA) at two different time points, to perform comparative alternative splicing analysis. Apparently, AS events increased under all stresses except dehydration for 1 h, and intron retention was the most common AS mode. In addition, a total of 357 differential alternative splicing (DAS) genes were identified under four abiotic stresses, among which 81 DAS genes existed in at least two stresses, and 276 DAS genes were presented under only one stress. A weighted gene co-expression network analysis (WGCNA) based on the splicing isoforms, rather than the genes, pinpointed out 23 co-expression modules associated with different abiotic stresses. Among them, a number of significant hub genes were also found to be DAS genes, which encode key isoforms involved in responses to single stress or multiple stresses, including RNA-binding proteins, transcription factors, and other important genes, such as RBP45C, LHY, MYB59, SCL30A, RS40, MAJ23.10, and DWF4. The splicing isoforms of candidate genes identified in this study could be a valuable resource for improving tolerance of B. napus against multiple abiotic stresses.

Alternative splicing reprogramming in fungal pathogen Sclerotinia sclerotiorum at different infection stages on Brassica napus

Alternative splicing (AS) is an important post-transcriptional mechanism promoting the diversity of transcripts and proteins to regulate various life processes in eukaryotes. Sclerotinia stem rot is a major disease of Brassica napus caused by Sclerotinia sclerotiorum, which causes severe yield loss in B. napus production worldwide. Although many transcriptome studies have been carried out on the growth, development, and infection of S. sclerotiorum, the genome-wide AS events of S. sclerotiorum remain poorly understood, particularly at the infection stage. In this study, transcriptome sequencing was performed to systematically explore the genome-scale AS events of S. sclerotiorum at five important infection stages on a susceptible oilseed rape cultivar. A total of 130 genes were predicted to be involved in AS from the S. sclerotiorum genome, among which 98 genes were differentially expressed and may be responsible for AS reprogramming for its successful infection. In addition, 641 differential alternative splicing genes (DASGs) were identified during S. sclerotiorum infection, accounting for 5.76% of all annotated S. sclerotiorum genes, and 71 DASGs were commonly found at all the five infection stages. The most dominant AS type of S. sclerotiorum was found to be retained introns or alternative 3' splice sites. Furthermore, the resultant AS isoforms of 21 DASGs became pseudogenes, and 60 DASGs encoded different putative proteins with different domains. More importantly, 16 DASGs of S. sclerotiorum were found to have signal peptides and possibly encode putative effectors to facilitate the infection of S. sclerotiorum. Finally, about 69.27% of DASGs were found to be non-differentially expressed genes, indicating that AS serves as another important way to regulate the infection of S. sclerotiorum on plants besides the gene expression level. Taken together, this study provides a genome-wide landscape for the AS of S. sclerotiorum during infection as well as an important resource for further elucidating the pathogenic mechanisms of S. sclerotiorum.

Ecotype-specific blockage of tasiARF production by two different RNA viruses in Arabidopsis

Arabidopsis thaliana is one of the most studied model organisms of plant biology with hundreds of geographical variants called ecotypes. One might expect that this enormous genetic variety could result in differential response to pathogens. Indeed, we observed previously that the Bur ecotype develops much more severe symptoms (upward curling leaves and wavy leaf margins) upon infection with two positive-strand RNA viruses of different families (turnip vein-clearing virus, TVCV, and turnip mosaic virus, TuMV). To find the genes potentially responsible for the ecotype-specific response, we performed a differential expression analysis of the mRNA and sRNA pools of TVCV and TuMV-infected Bur and Col plants along with the corresponding mock controls. We focused on the genes and sRNAs that showed an induced or reduced expression selectively in the Bur virus samples in both virus series. We found that the two ecotypes respond to the viral infection differently, yet both viruses selectively block the production of the TAS3-derived small RNA specimen called tasiARF only in the virus-infected Bur plants. The tasiARF normally forms a gradient through the adaxial and abaxial parts of the leaf (being more abundant in the adaxial part) and post-transcriptionally regulates ARF4, a major leaf polarity determinant in plants. The lack of tasiARF-mediated silencing could lead to an ectopically expressed ARF4 in the adaxial part of the leaf where the misregulation of auxin-dependent signaling would result in an irregular growth of the leaf blade manifesting as upward curling leaf and wavy leaf margin. QTL mapping using Recombinant Inbred Lines (RILs) suggests that the observed symptoms are the result of a multigenic interaction that allows the symptoms to develop only in the Bur ecotype. The particular nature of genetic differences leading to the ecotype-specific symptoms remains obscure and needs further study.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

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.

Comparative analysis of wild-type accessions reveals novel determinants of Arabidopsis seed longevity

Understanding the genetic factors involved in seed longevity is of paramount importance in agricultural and ecological contexts. The polygenic nature of this trait suggests that many of them remain undiscovered. Here, we exploited the contrasting seed longevity found amongst Arabidopsis thaliana accessions to further understand this phenomenon. Concentrations of glutathione were higher in longer-lived than shorter-lived accessions, supporting that redox poise plays a prominent role in seed longevity. However, high seed permeability, normally associated with shorter longevity, is also present in long-lived accessions. Dry seed transcriptome analysis indicated that the contribution to longevity of stored messenger RNA (mRNAs) is complex, including mainly accession-specific mechanisms. The detrimental effect on longevity caused by other factors may be counterbalanced by higher levels of specific mRNAs stored in dry seeds, for instance those of heat-shock proteins. Indeed, loss-of-function mutant analysis demonstrated that heat-shock factors HSF1A and 1B contributed to longevity. Furthermore, mutants of the stress-granule zinc-finger protein TZF9 or the spliceosome subunits MOS4 or MAC3A/MAC3B, extended seed longevity, positioning RNA as a novel player in the regulation of seed viability. mRNAs of proteins with putative relevance to longevity were also abundant in shorter-lived accessions, reinforcing the idea that resistance to ageing is determined by multiple factors.

PIF4 enhances DNA binding of CDF2 to co-regulate target gene expression and promote Arabidopsis hypocotyl cell elongation

How specificity is conferred within gene regulatory networks is an important problem in biology. The basic helix-loop-helix PHYTOCHROME-INTERACTING FACTORs (PIFs) and single zinc-finger CYCLING DOF FACTORs (CDFs) mediate growth responses of Arabidopsis to light and temperature. We show that these two classes of transcription factor (TF) act cooperatively. CDF2 and PIF4 are temporally and spatially co-expressed, they interact to form a protein complex and act in the same genetic pathway to promote hypocotyl cell elongation. Furthermore, PIF4 substantially strengthens genome-wide occupancy of CDF2 at a subset of its target genes. One of these, YUCCA8, encodes an auxin biosynthesis enzyme whose transcription is increased by PIF4 and CDF2 to contribute to hypocotyl elongation. The binding sites of PIF4 and CDF2 in YUCCA8 are closely spaced, and in vitro PIF4 enhances binding of CDF2. We propose that this occurs by direct protein interaction and because PIF4 binding alters DNA conformation. Thus, we define mechanisms by which PIF and CDF TFs cooperate to achieve regulatory specificity and promote cell elongation in response to light.