Alternative splicing in plants: directing traffic at the crossroads of adaptation and environmental stress

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

Divergence in regulatory mechanisms of GR-RBP genes in different plants under abiotic stress

The IVa subfamily of glycine-rich proteins (GRPs) comprises a group of glycine-rich RNA binding proteins referred to as GR-RBPa here. Previous studies have demonstrated functions of GR-RBPa proteins in regulating stress response in plants. However, the mechanisms responsible for the differential regulatory functions of GR-RBPa proteins in different plant species have not been fully elucidated. In this study, we identified and comprehensively studied a total of 34 GR-RBPa proteins from five plant species. Our analysis revealed that GR-RBPa proteins were further classified into two branches, with proteins in branch I being relatively more conserved than those in branch II. When subjected to identical stresses, these genes exhibited intensive and differential expression regulation in different plant species, corresponding to the enrichment of cis-acting regulatory elements involving in environmental and internal signaling in these genes. Unexpectedly, all GR-RBPa genes in branch I underwent intensive alternative splicing (AS) regulation, while almost all genes in branch II were only constitutively spliced, despite having more introns. This study highlights the complex and divergent regulations of a group of conserved RNA binding proteins in different plants when exposed to identical stress conditions. These species-specific regulations may have implications for stress responses and adaptations in different plant species.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

A dual-subcellular localized β-glucosidase confers pathogen and insect resistance without a yield penalty in maize

Maize is one of the most important crops for food, cattle feed and energy production. However, maize is frequently attacked by various pathogens and pests, which pose a significant threat to maize yield and quality. Identification of quantitative trait loci and genes for resistance to pests will provide the basis for resistance breeding in maize. Here, a β-glucosidase ZmBGLU17 was identified as a resistance gene against Pythium aphanidermatum, one of the causal agents of corn stalk rot, by genome-wide association analysis. Genetic analysis showed that both structural variations at the promoter and a single nucleotide polymorphism at the fifth intron distinguish the two ZmBGLU17 alleles. The causative polymorphism near the GT-AG splice site activates cryptic alternative splicing and intron retention of ZmBGLU17 mRNA, leading to the downregulation of functional ZmBGLU17 transcripts. ZmBGLU17 localizes in both the extracellular matrix and vacuole and contribute to the accumulation of two defence metabolites lignin and DIMBOA. Silencing of ZmBGLU17 reduces maize resistance against P. aphanidermatum, while overexpression significantly enhances resistance of maize against both the oomycete pathogen P. aphanidermatum and the Asian corn borer Ostrinia furnacalis. Notably, ZmBGLU17 overexpression lines exhibited normal growth and yield phenotype in the field. Taken together, our findings reveal that the apoplastic and vacuolar localized ZmBGLU17 confers resistance to both pathogens and insect pests in maize without a yield penalty, by fine-tuning the accumulation of lignin and DIMBOA.

SCR106 splicing factor modulates abiotic stress responses by maintaining RNA splicing in rice

Plants employ sophisticated molecular machinery to fine-tune their responses to growth, developmental, and stress cues. Gene expression influences plant cellular responses through regulatory processes such as transcription and splicing. Pre-mRNA is alternatively spliced to increase the genome coding potential and further regulate expression. Serine/arginine-rich (SR) proteins, a family of pre-mRNA splicing factors, recognize splicing cis-elements and regulate both constitutive and alternative splicing. Several studies have reported SR protein genes in the rice genome, subdivided into six subfamilies based on their domain structures. Here, we identified a new splicing factor in rice with an RNA recognition motif (RRM) and SR-dipeptides, which is related to the SR proteins, subfamily SC. OsSCR106 regulates pre-mRNA splicing under abiotic stress conditions. It localizes to the nuclear speckles, a major site for pre-mRNA splicing in the cell. The loss-of-function scr106 mutant is hypersensitive to salt, abscisic acid, and low-temperature stress, and harbors a developmental abnormality indicated by the shorter length of the shoot and root. The hypersensitivity to stress phenotype was rescued by complementation using OsSCR106 fused behind its endogenous promoter. Global gene expression and genome-wide splicing analysis in wild-type and scr106 seedlings revealed that OsSCR106 regulates its targets, presumably through regulating the alternative 3'-splice site. Under salt stress conditions, we identified multiple splice isoforms regulated by OsSCR106. Collectively, our results suggest that OsSCR106 is an important splicing factor that plays a crucial role in accurate pre-mRNA splicing and regulates abiotic stress responses in plants.

Phenylpropanoids Following Wounding and Infection of Sweet Sorghum Lines Differing in Responses to Stalk Pathogens

Sweet sorghum (Sorghum bicolor) lines M81-E and Colman were previously shown to differ in responses to Fusarium thapsinum and Macrophomina phaseolina, stalk rot pathogens that can reduce the yields and quality of biomass and extracted sugars. Inoculated tissues were compared for transcriptomic, phenolic metabolite, and enzymatic activity during disease development 3 and 13 days after inoculation (DAI). At 13 DAI, M81-E had shorter mean lesion lengths than Colman when inoculated with either pathogen. Transcripts encoding monolignol biosynthetic and modification enzymes were associated with transcriptional wound (control) responses of both lines at 3 DAI. Monolignol biosynthetic genes were differentially coexpressed with transcriptional activator SbMyb76 in all Colman inoculations, but only following M. phaseolina inoculation in M81-E, suggesting that SbMyb76 is associated with lignin biosynthesis during pathogen responses. In control inoculations, defense-related genes were expressed at higher levels in M81-E than Colman. Line, treatment, and timepoint differences observed in phenolic metabolite and enzyme activities did not account for observed differences in lesions. However, generalized additive models were able to relate metabolites, but not enzyme activities, to lesion length for quantitatively modeling disease progression: in M81-E, but not Colman, sinapic acid levels positively predicted lesion length at 3 DAI when cell wall-bound syringic acid was low, soluble caffeic acid was high, and lactic acid was high, suggesting that sinapic acid may contribute to responses at 3 DAI. These results provide potential gene targets for development of sweet sorghum varieties with increased stalk rot resistance to ensure biomass and sugar quality.

Alternative Splicing under Cold Stress in Paper Mulberry

The paper mulberry is a commonly found tree species with a long history of cultivation. It also serves as a crucial case study for understanding how woody plants adapt to low temperatures. Under cold treatment, we observed a substantial number of alternative splicing (AS) genes, showcasing the intricate landscape of AS events. We have detected all seven types of AS events, with the alternative 3' splice site (A3) having the most. We observed that many genes that underwent differential AS were significantly enriched in starch and sucrose metabolism and circadian rhythm pathways. Moreover, a considerable proportion of differentially spliced genes (DSGs) also showed differential expression, with 20.38% and 25.65% under 12 h and 24 h cold treatments, respectively. This suggests a coordinated regulation between gene AS and expression, playing a pivotal role in the paper mulberry's adaptation to cold stress. We further investigated the regulatory mechanisms of AS, identifying 41 serine/arginine-rich (SR) splicing factors, among which 11 showed differential expression under cold treatment, while 29 underwent alternative splicing. Additionally, genes undergoing AS displayed significantly higher DNA methylation levels under cold stress, while normal splicing (non-AS) genes exhibited relatively lower methylation levels. These findings suggest that methylation may play an important role in governing gene AS. Finally, our research will provide useful information on the role of AS in the cold acclimation tolerance of the paper mulberry.

Genome-Wide Analysis of Cation/Proton Antiporter Family in Soybean (Glycine max) and Functional Analysis of GmCHX20a on Salt Response

Monovalent cation proton antiporters (CPAs) play crucial roles in ion and pH homeostasis, which is essential for plant development and environmental adaptation, including salt tolerance. Here, 68 CPA genes were identified in soybean, phylogenetically dividing into 11 Na+/H+ exchangers (NHXs), 12 K+ efflux antiporters (KEAs), and 45 cation/H+ exchangers (CHXs). The GmCPA genes are unevenly distributed across the 20 chromosomes and might expand largely due to segmental duplication in soybean. The GmCPA family underwent purifying selection rather than neutral or positive selections. The cis-element analysis and the publicly available transcriptome data indicated that GmCPAs are involved in development and various environmental adaptations, especially for salt tolerance. Based on the RNA-seq data, twelve of the chosen GmCPA genes were confirmed for their differentially expression under salt or osmotic stresses using qRT-PCR. Among them, GmCHX20a was selected due to its high induction under salt stress for the exploration of its biological function on salt responses by ectopic expressing in Arabidopsis. The results suggest that the overexpression of GmCHX20a increases the sensitivity to salt stress by altering the redox system. Overall, this study provides comprehensive insights into the CPA family in soybean and has the potential to supply new candidate genes to develop salt-tolerant soybean varieties.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Alternative splicing: An efficient regulatory approach towards plant developmental plasticity

Alternative splicing (AS) is a gene regulatory mechanism that plants adapt to modulate gene expression (GE) in multiple ways. AS generates alternative isoforms of the same gene following various development and environmental stimuli, increasing transcriptome plasticity and proteome complexity. AS controls the expression levels of certain genes and regulates GE networks that shape plant adaptations through nonsense-mediated decay (NMD). This review intends to discuss AS modulation, from interaction with noncoding RNAs to the established roles of splicing factors (SFs) in response to endogenous and exogenous cues. We aim to gather such studies that highlight the magnitude and impact of AS, which are not always clear from individual articles, when AS is increasing in individual genes and at a global level. This work also anticipates making plant researchers know that AS is likely to occur in their investigations and that dynamic changes in AS and their effects must be frequently considered. We also review our understanding of AS-mediated posttranscriptional modulation of plant stress tolerance and discuss its potential application in crop improvement in the future. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > Splicing Mechanisms RNA-Based Catalysis > RNA Catalysis in Splicing and Translation.

Cation diffusion facilitator proteins of Beta vulgaris reveal diversity of metal handling in dicotyledons

Manganese (Mn), iron (Fe), and zinc (Zn) are essential for diverse processes in plants, but their availability is often limiting or excessive. Cation diffusion facilitator (CDF) proteins have been implicated in the allocation of those metals in plants, whereby most of our mechanistic understanding has been obtained in Arabidopsis. It is unclear to what extent this can be generalized to other dicots. We characterized all CDFs/metal tolerance proteins of sugar beet (Beta vulgaris spp. vulgaris), which is phylogenetically distant from Arabidopsis. Analysis of subcellular localization, substrate selectivities, and transcriptional regulation upon exposure to metal deficiencies and toxicities revealed unexpected deviations from their Arabidopsis counterparts. Localization and selectivity of some members were modulated by alternative splicing. Notably, unlike in Arabidopsis, Mn- and Zn-sequestrating members were not induced in Fe-deficient roots, pointing to differences in the Fe acquisition machinery. This was supported by low Zn and Mn accumulation under Fe deficiency and a strikingly increased Fe accumulation under Mn and Zn excess, coinciding with an induction of BvIRT1. High Zn load caused a massive upregulation of Zn-BvMTPs. The results suggest that the employment of the CDF toolbox is highly diverse amongst dicots, which questions the general applicability of metal homeostasis models derived from Arabidopsis.

De novo full-length transcriptome analysis of two ecotypes of Phragmites australis (swamp reed and dune reed) provides new insights into the transcriptomic complexity of dune reed and its long-term adaptation to desert environments

BACKGROUND

The extremely harsh environment of the desert is changing dramatically every moment, and the rapid adaptive stress response in the short term requires enormous energy expenditure to mobilize widespread regulatory networks, which is all the more detrimental to the survival of the desert plants themselves. The dune reed, which has adapted to desert environments with complex and variable ecological factors, is an ideal type of plant for studying the molecular mechanisms by which Gramineae plants respond to combinatorial stress of the desert in their natural state. But so far, the data on the genetic resources of reeds is still scarce, therefore most of their research has focused on ecological and physiological studies.

RESULTS

In this study, we obtained the first De novo non-redundant Full-Length Non-Chimeric (FLNC) transcriptome databases for swamp reeds (SR), dune reeds (DR) and the All of Phragmites australis (merged of iso-seq data from SR and DR), using PacBio Iso-Seq technology and combining tools such as Iso-Seq3 and Cogent. We then identified and described long non-coding RNAs (LncRNA), transcription factor (TF) and alternative splicing (AS) events in reeds based on a transcriptome database. Meanwhile, we have identified and developed for the first time a large number of candidates expressed sequence tag-SSR (EST-SSRs) markers in reeds based on UniTransModels. In addition, through differential gene expression analysis of wild-type and homogenous cultures, we found a large number of transcription factors that may be associated with desert stress tolerance in the dune reed, and revealed that members of the Lhc family have an important role in the long-term adaptation of dune reeds to desert environments.

CONCLUSIONS

Our results provide a positive and usable genetic resource for Phragmites australis with a widespread adaptability and resistance, and provide a genetic database for subsequent reeds genome annotation and functional genomic studies.

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.

The role of shoot-derived RNAs transported to plant root in response to abiotic stresses

A large number of RNA molecules are transported over long-distance between shoots and roots via phloem in higher plants. Mobile RNA signals are important for plants to tackle abiotic stresses. Shoot-derived mobile RNAs can be involved in the response to different developmental or environmental signals in the root. Some environmental conditions such as climate change, water deficit, nutrient deficiency challenge modern agriculture with more expeditious abiotic stress conditions. Root architecture determines the ability of water and nutrient uptake and further abiotic stress tolerance, and shoot tissue also determines the balance between shoot-root relationship in plant growth and adaptations. Thus, it is necessary to understand the roles of shoot-derived RNA signals and their potential function in roots upon abiotic stresses in the model plants (Arabidopsis thaliana and Nicotiana benthamiana) and agricultural crops. In this review, we summarize the so-far discovered shoot-derived mobile RNA transportation to the root under abiotic stress conditions, e.g. drought, cold stress and nutrient deficiencies. Furthermore, we will focus on the biological relevance and the potential roles of these RNAs in root development and stress responses which will be an asset for the future breeding strategies.

Dysregulated RNA processing and metabolism: a new hallmark of ageing and provocation for cellular senescence

The human genome is capable of producing hundreds of thousands of different proteins and non-coding RNAs from <20 000 genes, in a co-ordinated and regulated fashion. This is achieved by a collection of phenomena known as mRNA processing and metabolism, and encompasses events in the life cycle of an RNA from synthesis to degradation. These factors are critical determinants of cellular adaptability and plasticity, which allows the cell to adjust its transcriptomic output in response to its internal and external environment. Evidence is building that dysfunctional RNA processing and metabolism may be a key contributor to the development of cellular senescence. Senescent cells by definition have exited cell cycle, but have gained functional features such as the secretion of the senescence-associated secretory phenotype (SASP), a known driver of chronic disease and perhaps even ageing itself. In this review, I will outline the impact of dysregulated mRNA processing and metabolism on senescence and ageing at the level of genes, cells and systems, and describe the mechanisms by which progressive deterioration in these processes may impact senescence and organismal ageing. Finally, I will present the evidence implicating this important process as a new hallmark of ageing, which could be harnessed in the future to develop new senotherapeutic interventions for chronic disease.

Long non-coding RNAs reveal new regulatory mechanisms controlling gene expression

A plethora of non-coding RNAs have been found in eukaryotes, notably with the advent of modern sequencing technologies to analyze the transcriptome. Apart from the well-known housekeeping RNA genes (such as the ribosomal RNA or the transfer RNA), many thousands of transcripts detected are not evidently linked to a protein-coding gene. These, so called non-coding RNAs, may code for crucial regulators of gene expression, the small si/miRNAs, for small peptides (translated under specific conditions) or may act as long RNA molecules (antisense, intronic or intergenic long non-coding RNAs or lncRNAs). The lncRNAs interact with members of multiple machineries involved in gene regulation. In this review, we discussed about how plant lncRNAs permitted to discover new regulatory mechanisms acting in epigenetic control, chromatin 3D structure and alternative splicing. These novel regulations diversified the expression patterns and protein variants of target protein-coding genes and are an important element of the response of plants to environmental stresses and their adaptation to changing conditions.

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

Co-expression network of heat-response transcripts: A glimpse into how splicing factors impact rice basal thermotolerance

To identify novel solutions to improve rice yield under rising temperatures, molecular components of thermotolerance must be better understood. Alternative splicing (AS) is a major post-transcriptional mechanism impacting plant tolerance against stresses, including heat stress (HS). AS is largely regulated by splicing factors (SFs) and recent studies have shown their involvement in temperature response. However, little is known about the splicing networks between SFs and AS transcripts in the HS response. To expand this knowledge, we constructed a co-expression network based on a publicly available RNA-seq dataset that explored rice basal thermotolerance over a time-course. Our analyses suggest that the HS-dependent control of the abundance of specific transcripts coding for SFs might explain the widespread, coordinated, complex, and delicate AS regulation of critical genes during a plant's inherent response to extreme temperatures. AS changes in these critical genes might affect many aspects of plant biology, from organellar functions to cell death, providing relevant regulatory candidates for future functional studies of basal thermotolerance.

Post-transcriptional regulation of 2-acetyl-1-pyrroline (2-AP) biosynthesis pathway, silicon, and heavy metal transporters in response to Zn in fragrant rice

Fragrant rice (Oryza sativa L.) has a high economic and nutritional value, and the application of micronutrients regulates 2-acetyl-1-pyrroline (2-AP) production, which is responsible for aroma in fragrant rice. Alternative splicing (AS) is an important post-transcriptional regulatory mechanism to generate transcript variability and proteome diversity in plants. However, no systematic investigation of AS events in response to micronutrients (Zn) has been performed in fragrant rice. Furthermore, the post-transcriptional regulation of genes involved in 2-AP biosynthesis is also not known. In this study, a comprehensive analysis of AS events under two gradients of Zn treatment in two different fragrant rice cultivars (Meixiangzhan-2 and Xiangyaxiangzhan) was performed based on RNA-seq analysis. A total of 386 and 598 significant AS events were found in Meixiangzhan-2 treated with low and high doses of Zn, respectively. In Xiangyaxiangzhan, a total of 449 and 598 significant AS events were found in low and high doses of Zn, respectively. Go analysis indicated that these genes were highly enriched in physiological processes, metabolism, and cellular processes in both cultivars. However, genotype and dose-dependent AS events were also detected in both cultivars. By comparing differential AS (DAS) events with differentially expressed genes (DEGs), we found a weak overlap among DAS and DEGs in both fragrant rice cultivars indicating that only a few genes are post-transcriptionally regulated in response to Zn treatment. We further report that Zn differentially regulates the expression of 2-AP biosynthesis-related genes in both cultivars and Zn treatment altered the editing frequency of single nucleotide polymorphism (SNPs) in the genes involved in 2-AP biosynthesis. Finally, we showed that epigenetic modifications associated with active gene transcription are generally enriched over 2-AP biosynthesis-related genes. Similar to the 2-AP pathway, we found that heavy metal transporters (genes related to silicon, iron, Zn and other metal transport) are also regulated at transcriptional and post-transcriptional levels in response to Zn in fragrant rice. Taken together, our results provide evidence of the post-transcriptional gene regulation in fragrant rice in response to Zn treatment and highlight that the 2-AP biosynthesis pathway and heavy metal transporters may also be regulated through epigenetic modifications. These findings will serve as a cornerstone for further investigation to understand the molecular mechanisms of 2-AP biosynthesis and regulation of heavy metal transporters in fragrant rice.

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

Identifying Terpenoid Biosynthesis Genes in Euphorbia maculata via Full-Length cDNA Sequencing

The annual herb Euphorbia maculata L. produces anti-inflammatory and biologically active substances such as triterpenoids, tannins, and polyphenols, and it is used in traditional Chinese medicine. Of these bioactive compounds, terpenoids, also called isoprenoids, are major secondary metabolites in E. maculata. Full-length cDNA sequencing was carried out to characterize the transcripts of terpenoid biosynthesis reference genes and determine the copy numbers of their isoforms using PacBio SMRT sequencing technology. The Illumina short-read sequencing platform was also employed to identify differentially expressed genes (DEGs) in the secondary metabolite pathways from leaves, roots, and stems. PacBio generated 62 million polymerase reads, resulting in 81,433 high-quality reads. From these high-quality reads, we reconstructed a genome of 20,722 genes, in which 20,246 genes (97.8%) did not have paralogs. About 33% of the identified genes had two or more isoforms. DEG analysis revealed that the expression level differed among gene paralogs in the leaf, stem, and root. Whole sets of paralogs and isoforms were identified in the mevalonic acid (MVA), methylerythritol phosphate (MEP), and terpenoid biosynthesis pathways in the E. maculata L. The nucleotide information will be useful for identifying orthologous genes in other terpenoid-producing medicinal plants.

Function of Cajal Bodies in Nuclear RNA Retention in A. thaliana Leaves Subjected to Hypoxia

Retention of RNA in the nucleus precisely regulates the time and rate of translation and controls transcriptional bursts that can generate profound variability in mRNA levels among identical cells in tissues. In this study, we investigated the function of Cajal bodies (CBs) in RNA retention in A. thaliana leaf nuclei during hypoxia stress was investigated. It was observed that in ncb-1 mutants with a complete absence of CBs, the accumulation of poly(A⁺) RNA in the leaf nuclei was lower than that in wt under stress. Moreover, unlike in root cells, CBs store less RNA, and RNA retention in the nuclei is much less intense. Our results reveal that the function of CBs in the accumulation of RNA in nuclei under stress depends on the plant organ. Additionally, in ncb-1, retention of introns of mRNA RPB1 (largest subunit of RNA polymerase II) mRNA was observed. However, this isoform is highly accumulated in the nucleus. It thus follows that intron retention in transcripts is more important than CBs for the accumulation of RNA in nuclei. Accumulated mRNAs with introns in the nucleus could escape transcript degradation by NMD (nonsense-mediated mRNA decay). From non-fully spliced mRNAs in ncb-1 nuclei, whose levels increase during hypoxia, introns are removed during reoxygenation. Then, the mRNA is transferred to the cytoplasm, and the RPB1 protein is translated. Despite the accumulation of isoforms in nuclei with retention of introns in reoxygenation, ncb-1 coped much worse with long hypoxia, and manifested faster yellowing and shrinkage of leaves.

Full-Length Transcriptome Sequencing Reveals the Impact of Cold Stress on Alternative Splicing in Quinoa

Quinoa is a cold-resistant and nutrient-rich crop. To decipher the cold stress response of quinoa, the full-length transcriptomes of the cold-resistant quinoa variety CRQ64 and the cold-sensitive quinoa variety CSQ5 were compared. We identified 55,389 novel isoforms and 6432 novel genes in these transcriptomes. Under cold stress, CRQ64 had more differentially expressed genes (DEGs) and differentially alternative splicing events compared to non-stress conditions than CSQ5. DEGs that were specifically present only in CRQ64 were significantly enriched in processes which contribute to osmoregulation and ROS homeostasis in plants, such as sucrose metabolism and phenylpropanoid biosynthesis. More genes with differential alternative splicing under cold stress were enriched in peroxidase functions in CRQ64. In total, 5988 transcription factors and 2956 long non-coding RNAs (LncRNAs) were detected in this dataset. Many of these had altered expression patterns under cold stress compared to non-stress conditions. Our transcriptome results demonstrate that CRQ64 undergoes a wider stress response than CSQ5 under cold stress. Our results improved the annotation of the quinoa genome and provide new insight into the mechanisms of cold resistance in quinoa.

An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development

Background

Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies.

Result

In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses.

Conclusions

This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02620-2.

Analysis of Alternative Splicing During the Combinatorial Response to Simultaneous Copper and Iron Deficiency in Arabidopsis Reveals Differential Events in Genes Involved in Amino Acid Metabolism

Copper (Cu) and iron (Fe) constitute fundamental nutrients for plant biology but are often limited due to low bioavailability. Unlike responses to single Cu or Fe deprivation, the consequences of simultaneous Cu and Fe deficiency have not yet been fully deciphered. Previously, it was demonstrated that Cu and Fe deficiency applied in combination imposes transcriptome, proteome, and metabolome changes different from those triggered under each deficiency individually. Here, we evaluated the effect of alternative splicing (AS) on the transcriptome of rosette leaves under single and simultaneous Cu and Fe deficiency. Differentially spliced genes (DSGs) and differentially expressed genes (DEGs) coincided in number (2,600 approx.) although the overlapping fraction was minimal (15%). Functional annotation of changes exclusively detected under simultaneous Cu and Fe deficiency revealed that DEGs participated in general stress responses and translation, while DSGs were involved in metabolic reactions, especially amino acid biosynthesis. Interestingly, transcripts encoding central features for tryptophan (Trp) and asparagine (Asn) synthesis - two significantly altered metabolites under simultaneous Cu and Fe deficiency - underwent exclusive intron retention events under the double deficiency. However, transcript and protein amounts for these enzymes did not correlate with Trp and Asn concentration. In consequence, we propose that AS might act as a regulatory mechanism to modify the stability and/or functionality of the enzymes and therefore fine-tune amino acid production during the combinatorial response to simultaneous Cu and Fe deficiency.

Cold Response Transcriptome Analysis of the Alternative Splicing Events Induced by the Cold Stress in D. catenatum

Dendrobium catenatum Lindl is a valuable medicinal herb and gardening plant due to its ornamental value and special medical value. Low temperature is a major bottleneck restricting D. catenatum expansion towards the north, which influences the quality and yield of D. catenatum. In this study, we analysed the cold response of D. catenatum by RNA-Seq. A total of 4302 differentially expressed genes were detected under cold stress, which were mainly linked to protein kinase activity, membrane transport and the glycan biosynthesis and metabolism pathway. We also identified 4005 differential alternative events in 2319 genes significantly regulated by cold stress. Exon skipping and intron retention were the most common alternative splicing isoforms. Numerous genes were identified that differentially modulated under cold stress, including cold-induced transcription factors and splicing factors mediated by AS (alternative splicing). GO enrichment analysis found that differentially alternatively spliced genes without differential expression levels were related to RNA/mRNA processing and spliceosomes. DAS (differentially alternative splicing) genes with different expression levels were mainly enriched in protein kinase activity, plasma membrane and cellular response to stimulus. We further identified and cloned DcCBP20 in D. catenatum; we found that DcCBP20 promotes the generation of alternative splicing variants in cold-induced genes under cold stress via genetic experiments and RT–PCR. Taken together, our results identify the main cold-response pathways and alternative splicing events in D. catenatum responding to cold treatment and that DcCBP20 of D. catenatum get involved in regulating the AS and gene expression of cold-induced genes during this process. Our study will contribute to understanding the role of AS genes in regulating the cold stress response in D. catenatum.

How alternative splicing changes the properties of plant proteins

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

Essentiality for rice fertility and alternative splicing of OsSUT1

Sucrose-proton symporters play important roles in carbohydrate transport during plant growth and development. Their physiological functions have only been partly characterized and their regulation mechanism is largely unclear. Here we report that the knockout of a sucrose transporter gene, OsSUT1, by CRISPR-Cas9 mediated gene editing resulted in a slightly dwarf size and complete infertility of the gene's homozygous mutants. Observation of caryopsis development revealed that the endosperm of OsSUT1 mutants failed to cellularize and did not show any sign of seed-filling. Consistently, OsSUT1 was identified to express strongly in developing caryopsis of wild-type rice, particularly in the nucellar epidermis and aleurone which are critical for the uptake of nutrients into the endosperm. These results indicate that OsSUT1 is indispensable during the rice reproductive stage particularly for caryopsis development. Interestingly, OsSUT1 possesses at least 6 alternative splicing transcripts, including the 4 transcripts deposited previously and the other two identified by us. The differences among these transcripts primarily lie in their coding region of the 3' end and 3' UTR region. Real-time PCR showed that 4 of the 6 transcripts had different expressional patterns during rice vegetative and reproductive growth stages. Given the versatility of the gene, addressing its alternative splicing mechanism may expand our understanding of SUT's function substantially.

Temporal regulation of alternative splicing events in rice memory under drought stress

Plant adaptation to drought stress is essential for plant survival and crop yield. Recently, harnessing drought memory, which is induced by repeated stress and recovery cycles, was suggested as a means to improve drought resistance at the transcriptional level. However, the genetic mechanism underlying drought memory is unclear. Here, we carried out a quantitative analysis of alternative splicing (AS) events in rice memory under drought stress, generating 12 transcriptome datasets. Notably, we identified exon skipping (ES) as the predominant AS type (>80%) in differential alternative splicing (DAS) in response to drought stress. Applying our analysis pipeline to investigate DAS events following drought stress in six other plant species revealed variable ES frequencies ranging from 9.94% to 60.70% depending on the species, suggesting that the relative frequency of DAS types in plants is likely to be species-specific. The dinucleotide sequence at AS splice sites in rice following drought stress was preferentially GC-AG and AT-AC. Since U12-type splicing uses the AT-AC site, this suggests that drought stress may increase U12-type splicing, and thus increase ES frequency. We hypothesize that multiple isoforms derived from exon skipping may be induced by drought stress in rice. We also identified 20 transcription factors and three highly connected hub genes with potential roles in drought memory that may be good targets for plant breeding.

Salt responsive alternative splicing of a RING finger E3 ligase modulates the salt stress tolerance by fine-tuning the balance of COP9 signalosome subunit 5A

Increasing evidence points to the tight relationship between alternative splicing (AS) and the salt stress response in plants. However, the mechanisms linking these two phenomena remain unclear. In this study, we have found that Salt-Responsive Alternatively Spliced gene 1 (SRAS1), encoding a RING-Type E3 ligase, generates two splicing variants: SRAS1.1 and SRAS1.2, which exhibit opposing responses to salt stress. The salt stress-responsive AS event resulted in greater accumulation of SRAS1.1 and a lower level of SRAS1.2. Comprehensive phenotype analysis showed that overexpression of SRAS1.1 made the plants more tolerant to salt stress, whereas overexpression of SRAS1.2 made them more sensitive. In addition, we successfully identified the COP9 signalosome 5A (CSN5A) as the target of SRAS1. CSN5A is an essential player in the regulation of plant development and stress. The full-length SRAS1.1 promoted degradation of CSN5A by the 26S proteasome. By contrast, SRAS1.2 protected CSN5A by competing with SRAS1.1 on the same binding site. Thus, the salt stress-triggered AS controls the ratio of SRAS1.1/SRAS1.2 and switches on and off the degradation of CSN5A to balance the plant development and salt tolerance. Together, these results provide insights that salt-responsive AS acts as post-transcriptional regulation in mediating the function of E3 ligase.

High salinity severely affects plant growth and development, impairing crop production worldwide. E3 ligase is a stress-responsive regulator through ubiquitin-proteasome system for selective protein degradation. The E3s are regulated by transcriptional regulation and post-translational modifications. Here, we have discovered that stress-responsive AS acts as a post-transcriptional regulation modulating the function of E3 ligases. Intriguingly, the truncated proteins generated by salt-responsive AS play opposite roles compared with the full-length E3 ligase. The truncated isoform losing key domain could not degrade the target protein, instead, it interacts and competes with the E3 ligase through binding the same domain of the targets. This finding contributes significantly to a deeper mechanistic understanding of how AS regulates the function of E3 ligase in response to salt stress.

Regulation of alternative splicing in response to temperature variation in plants

Plants have evolved numerous molecular strategies to cope with perturbations in environmental temperature, and to adjust growth and physiology to limit the negative effects of extreme temperature. One of the strategies involves alternative splicing of primary transcripts to encode alternative protein products or transcript variants destined for degradation by nonsense-mediated decay. Here, we review how changes in environmental temperature-cold, heat, and moderate alterations in temperature-affect alternative splicing in plants, including crops. We present examples of the mode of action of various temperature-induced splice variants and discuss how these alternative splicing events enable favourable plant responses to altered temperatures. Finally, we point out unanswered questions that should be addressed to fully utilize the endogenous mechanisms in plants to adjust their growth to environmental temperature. We also indicate how this knowledge might be used to enhance crop productivity in the future.

Intron-retained radish (Raphanus sativus L.) RsMYB1 transcripts found in colored-taproot lines enhance anthocyanin accumulation in transgenic Arabidopsis plants

KEY MESSAGE

Overexpression of the naturally occurring intron-retained (IR) forms of radish RsMYB1 and RsTT8 transcripts in Arabidopsis causes a substantial increase in anthocyanin accumulation. The production of anthocyanins in plants is tightly controlled by the MYB-bHLH-WD40 (MBW) complex. In this study, analysis of four radish (Raphanus sativus L.) inbred lines with different colored taproots revealed that regulatory genes of anthocyanin biosynthesis, RsMYB1 and RsTT8, produce three transcripts, one completely spliced and two intron retention (IR1 and IR2) forms. Transcripts RsMYB1-IR1 and RsMYB1-IR2 retained the 1st (380 nt) and 2nd (149 nt) introns, respectively; RsTT8-IR1 retained the 4th intron (113 nt); RsTT8-IR2 retained both the 3rd (128 nt) and 4th introns. Levels of most IR forms were substantially low in radish samples, but the RsTT8-IR2 level was higher than RsTT8 in red skin/red flesh (RsRf) root. Since all IR forms contained a stop codon within the intron, they were predicted to encode truncated proteins with defective interaction domains, resulting in the inability to form the MBW complex in vivo. However, tobacco leaves transiently co-expressing RsMYB1-IRs and RsTT8-IRs showed substantially higher anthocyanin accumulation than those co-expressing their spliced forms. Consistently, co-expression of constructs encoding truncated proteins with spliced or IR forms of their interaction partner in tobacco leaves did not result in anthocyanin accumulation. Compared with RsMYB1, the overexpression of RsMYB1-IRs in Arabidopsis pap1 mutant increased anthocyanin accumulation by > sevenfold and upregulated the expression of Arabidopsis flavonoid biosynthesis genes including AtTT8. Our results suggest that the stable co-expression of RsMYB1-IRs in fruit trees and vegetable crops could be used to increase their anthocyanin contents.

Two FD homologs from London plane (Platanus acerifolia) are associated with floral initiation and flower morphology

The flowering-time gene FD encodes a bZIP transcription factor that interacts with FLOWERING LOCUS T (FT) to induce flowering in Arabidopsis. Previous research has identified two FT homologs of Platanus acerifolia, PaFT and PaFTL, which each have different expression patterns and are involved in diverse developmental processes. However, it is not known whether such FT/FD complexes participate in the flowering processes in P. acerifolia. Therefore, we isolated two closely related FD homologs, PaFDL1 and PaFDL2, and investigated their functions through the analysis of expression profiles, transgenic phenotypes, their interactions with different FT proteins, and potential cis-regulatory elements in their promoters. The PaFDL genes were found to display their maximal expression levels during the stage of floral transition, and subsequent expression patterns were also seen to be related to inflorescence developmental stage. In addition, both PaFDL1 and PaFDL2 were found to be subject to post-transcriptional alternative splicing, each gene producing two transcript forms. Transgenic tobacco overexpressing each of the four resulting transcript types displayed accelerated floral initiation and produced abnormal flowers. The results suggested that the complete PaFDL proteins may interact with different PaFT/PaFTL proteins in order to fulfill both conservative and diverse functions in floral initiation and floral development.

Reading between the Lines: RNA-seq Data Mining Reveals the Alternative Message of the Rice Leaf Transcriptome in Response to Heat Stress

Rice (Oryza sativa L.) is a major food crop but heat stress affects its yield and grain quality. To identify mechanistic solutions to improve rice yield under rising temperatures, molecular responses of thermotolerance must be understood. Transcriptional and post-transcriptional controls are involved in a wide range of plant environmental responses. Alternative splicing (AS), in particular, is a widespread mechanism impacting the stress defence in plants but it has been completely overlooked in rice genome-wide heat stress studies. In this context, we carried out a robust data mining of publicly available RNA-seq datasets to investigate the extension of heat-induced AS in rice leaves. For this, datasets of interest were subjected to filtering and quality control, followed by accurate transcript-specific quantifications. Powerful differential gene expression (DE) and differential AS (DAS) identified 17,143 and 2162 heat response genes, respectively, many of which are novel. Detailed analysis of DAS genes coding for key regulators of gene expression suggests that AS helps shape transcriptome and proteome diversity in response to heat. The knowledge resulting from this study confirmed a widespread transcriptional and post-transcriptional response to heat stress in plants, and it provided novel candidates for rapidly advancing rice breeding in response to climate change.

Targeting Alternative Splicing for Reversal of Cellular Senescence in the Context of Aesthetic Aging

SUMMARY

Cellular senescence is a state of stable cell cycle arrest that has increasingly been linked with cellular, tissue, and organismal aging; targeted removal of senescent cells brings healthspan and lifespan benefits in animal models. Newly emerging approaches to specifically ablate or rejuvenate senescent cells are now the subject of intense study to explore their utility to provide novel treatments for the aesthetic signs and diseases of aging in humans. Here, we discuss different strategies that are being trialed in vitro, and more recently in vivo, for the targeted removal or reversal of senescent cells. Finally, we describe the evidence for a newly emerging molecular mechanism that may underpin senescence; dysregulation of alternative splicing. We will explore the potential of restoring splicing regulation as a novel "senotherapeutic" approach and discuss strategies by which this could be integrated into the established portfolio of skin aging therapeutics.

Overlapping roles of spliceosomal components SF3B1 and PHF5A in rice splicing regulation

The SF3B complex, a multiprotein component of the U2 snRNP of the spliceosome, plays a crucial role in recognizing branch point sequence and facilitates spliceosome assembly and activation. Several chemicals that bind SF3B1 and PHF5A subunits of the SF3B complex inhibit splicing. We recently generated a splicing inhibitor-resistant SF3B1 mutant named SF3B1GEX1ARESISTANT 4 (SGR4) using CRISPR-mediated directed evolution, whereas splicing inhibitor-resistant mutant of PHF5A (Overexpression-PHF5A GEX1A Resistance, OGR) was generated by expressing an engineered version PHF5A-Y36C. Global analysis of splicing in wild type and these two mutants revealed the role of SF3B1 and PHF5A in splicing regulation. This analysis uncovered a set of genes whose intron retention is regulated by both proteins. Further analysis of these retained introns revealed that they are shorter, have a higher GC content, and contain shorter and weaker polypyrimidine tracts. Furthermore, splicing inhibition increased seedlings sensitivity to salt stress, consistent with emerging roles of splicing regulation in stress responses. In summary, we uncovered the functions of two members of the plant branch point recognition complex. The novel strategies described here should be broadly applicable in elucidating functions of splicing regulators, especially in studying the functions of redundant paralogs in plants.

Butt et al. used CRISPR-mediated directed evolution to generate rice mutants for the spliceosome components SF3B1 and PHF5A. They demonstrate that these mutants have different levels of sensitivity to salt treatments and suggest that the strategies they employed can be used in the future to study functions of redundant paralogs in plants.

Dual roles of the serine/arginine-rich splicing factor SR45a in promoting and interacting with nuclear cap-binding complex to modulate the salt-stress response in Arabidopsis

Alternative splicing (AS) is emerging as a critical co-transcriptional regulation for plants in response to environmental stresses. Although multiple splicing factors have been linked to the salt-sensitive signaling network, the molecular mechanism remains unclear. We discovered that a conserved serine/arginine-rich (SR)-like protein, SR45a, as a component of the spliceosome, was involved in post-transcriptional regulation of salinity tolerance in Arabidopsis thaliana. Furthermore, SR45a was required for the AS and messenger RNA (mRNA) maturation of several salt-tolerance genes. Two alternatively spliced variants of SR45a were induced by salt stress, full-length SR45a-1a and the truncated isoform SR45a-1b, respectively. Lines with overexpression of SR45a-1a and SR45a-1b exhibited hypersensitive to salt stress. Our data indicated that SR45a directly interacted with the cap-binding complex (CBC) subunit cap-binding protein 20 (CBP20) which mediated salt-stress responses. Instead of binding to other spliceosome components, SR45a-1b promoted the association of SR45a-1a with CBP20, therefore mediating salt-stress signal transduction pathways. Additionally, the mutations in SR45a and CBP20 led to different salt-stress phenotypes. Together, these results provide the evidence that SR45a-CBP20 acts as a regulatory complex to regulate the plant response to salt stress, through a regulatory mechanism to fine-tune the splicing factors, especially in stressful conditions.

Genome-Wide Investigation of N6-Methyladenosine Regulatory Genes and Their Roles in Tea (Camellia sinensis) Leaves During Withering Process

N6-methyladenosine (m6A), one of the internal modifications of RNA molecules, can directly influence RNA abundance and function without altering the nucleotide sequence, and plays a pivotal role in response to diverse environmental stresses. The precise m6A regulatory mechanism comprises three types of components, namely, m6A writers, erasers, and readers. To date, the research focusing on m6A regulatory genes in plant kingdom is still in its infancy. Here, a total of 34 m6A regulatory genes were identified from the chromosome-scale genome of tea plants. The expansion of m6A regulatory genes was driven mainly by whole-genome duplication (WGD) and segmental duplication, and the duplicated gene pairs evolved through purifying selection. Gene structure analysis revealed that the sequence variation contributed to the functional diversification of m6A regulatory genes. Expression pattern analysis showed that most m6A regulatory genes were differentially expressed under environmental stresses and tea-withering stage. These observations indicated that m6A regulatory genes play essential roles in response to environmental stresses and tea-withering stage. We also found that RNA methylation and DNA methylation formed a negative feedback by interacting with each other's methylation regulatory genes. This study provided a foundation for understanding the m6A-mediated regulatory mechanism in tea plants under environmental stresses and tea-withering stage.

A new mode of NPR1 action via an NB-ARC-NPR1 fusion protein negatively regulates the defence response in wheat to stem rust pathogen

NPR1 has been found to be a key transcriptional regulator in some plant defence responses. There are nine NPR1 homologues (TaNPR1) in wheat, but little research has been done to understand the function of those NPR1-like genes in the wheat defence response against stem rust (Puccinia graminis f. sp. tritici) pathogens. We used bioinformatics and reverse genetics approaches to study the expression and function of each TaNPR1. We found six members of TaNPR1 located on homoeologous group 3 chromosomes (designated as TaG3NPR1) and three on homoeologous group 7 chromosomes (designated as TaG7NPR1). The group 3 NPR1 proteins regulate transcription of SA-responsive PR genes. Downregulation of all the TaNPR1 homologues via virus-induced gene co-silencing resulted in enhanced resistance to stem rust. More specifically downregulating TaG7NPR1 homeologues or Ta7ANPR1 expression resulted in stem rust resistance phenotype. By contrast, knocking down TaG3NPR1 alone did not show visible phenotypic changes in response to the rust pathogen. Knocking out Ta7ANPR1 enhanced resistance to stem rust. The Ta7ANPR1 locus is alternatively spliced under pathogen inoculated conditions. We discovered a new mode of NPR1 action in wheat at the Ta7ANPR1 locus through an NB-ARC-NPR1 fusion protein negatively regulating the defence to stem rust infection.

Alternative splicing and its regulatory role in woody plants

Alternative splicing (AS) is an important post-transcriptional process to enhance proteome diversity in eukaryotic organisms. In plants, numerous reports have primarily focused on AS analysis in model plant species or herbaceous plants, leading to a notable lack of research on AS in woody plants. More importantly, emerging evidence indicates that many important traits, including wood formation and stress resistance, in woody plants are controlled by AS. In this review article, we summarize the current progress of all kinds of AS studies in different tree species at various stages of development and in response to various stresses, revealing the significant role played by AS in woody plants, as well as the similar properties and differential regulation within their herbaceous counterparts. Furthermore, we propose several potential strategies to facilitate the functional characterization of splicing factors in woody plants and evaluate a general pipeline for the systematic characterization of splicing isoforms in a complex AS regulatory network. The utilization of genetic studies and high-throughput omics integration approaches to analyze AS genes and splicing factors is likely to further advance our understanding of AS modulation in woody plants.

Investigation for a multi-silique trait in Brassica napus by alternative splicing analysis

Background

Flower and fruit development are vital stages of the angiosperm lifecycle. We previously investigated the multi-silique trait in the rapeseed (Brassica napus) line zws-ms on a genomic and transcriptomic level, leading to the identification of two genomic regions and several candidate genes associated with this trait. However, some events on the transcriptome level, like alternative splicing, were poorly understood.

Methods

Plants from zws-ms and its near-isogenic line (NIL) zws-217 were both grown in Xindu with normal conditions and a colder area Ma'erkang. Buds from the two lines were sampled and RNA was isolated to perform the transcriptomic sequencing. The numbers and types of alternative splicing (AS) events from the two lines were counted and classified. Genes with AS events and expressed differentially between the two lines, as well as genes with AS events which occurred in only one line were emphasized. Their annotations were further studied.

Results

From the plants in Xindu District, an average of 205,496 AS events, which could be sorted into five AS types, were identified. zws-ms and zws-217 shared highly similar ratios of each AS type: The alternative 5' and 3' splice site types were the most common, while the exon skipping type was observed least often. Eleven differentially expressed AS genes were identified, of which four were upregulated and seven were downregulated in zws-ms. Their annotations implied that five of these genes were directly associated with the multi-silique trait. While samples from colder area Ma'erkang generated generally reduced number of each type of AS events except for Intron Retention; but the number of differentially expressed AS genes increased significantly. Further analysis found that among the 11 differentially expressed AS genes from Xindu, three of them maintained the same expression models, while the other eight genes did not show significant difference between the two lines in expression level. Additionally, the 205 line-specific expressed AS genes were analyzed, of which 187 could be annotated, and two were considered to be important.

Discussion

This study provides new insights into the molecular mechanism of the agronomically important multi-silique trait in rapeseed on the transcriptome level and screens out some environment-responding candidate genes.

The CBL-CIPK Pathway in Plant Response to Stress Signals

Plants need to cope with multitudes of stimuli throughout their lifecycles in their complex environments. Calcium acts as a ubiquitous secondary messenger in response to numerous stresses and developmental processes in plants. The major Ca2+ sensors, calcineurin B-like proteins (CBLs), interact with CBL-interacting protein kinases (CIPKs) to form a CBL-CIPK signaling network, which functions as a key component in the regulation of multiple stimuli or signals in plants. In this review, we describe the conserved structure of CBLs and CIPKs, characterize the features of classification and localization, draw conclusions about the currently known mechanisms, with a focus on novel findings in response to multiple stresses, and summarize the physiological functions of the CBL-CIPK network. Moreover, based on the gradually clarified mechanisms of the CBL-CIPK complex, we discuss the present limitations and potential prospects for future research. These aspects may provide a deeper understanding and functional characterization of the CBL-CIPK pathway and other signaling pathways under different stresses, which could promote crop yield improvement via biotechnological intervention.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure

BACKGROUND

Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses.

RESULTS

By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing.

CONCLUSIONS

Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Gene editing: an instrument for practical application of gene biology to plant breeding

Plant breeding is well recognized as one of the most important means to meet food security challenges caused by the ever-increasing world population. During the past three decades, plant breeding has been empowered by both new knowledge on trait development and regulation (e.g., functional genomics) and new technologies (e.g., biotechnologies and phenomics). Gene editing, particularly by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) and its variants, has become a powerful technology in plant research and may become a game-changer in plant breeding. Traits are conferred by coding and non-coding genes. From this perspective, we propose different editing strategies for these two types of genes. The activity of an encoded enzyme and its quantity are regulated at transcriptional and post-transcriptional, as well as translational and post-translational, levels. Different strategies are proposed to intervene to generate gene functional variations and consequently phenotype changes. For non-coding genes, trait modification could be achieved by regulating transcription of their own or target genes via gene editing. Also included is a scheme of protoplast editing to make gene editing more applicable in plant breeding. In summary, this review provides breeders with a host of options to translate gene biology into practical breeding strategies, i.e., to use gene editing as a mechanism to commercialize gene biology in plant breeding.

A Collection of Pre-mRNA Splicing Mutants in Arabidopsis thaliana

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.

Role of MPK4 in pathogen-associated molecular pattern-triggered alternative splicing in Arabidopsis

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied alternative splicing induced by the PAMP flagellin in Arabidopsis. We identified 506 PAMP-induced differentially alternatively spliced (DAS) genes. Importantly, of the 506 PAMP-induced DAS genes, only 89 overlap with the set of 1950 PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DAS events. Global DAS analysis of mpk3, mpk4, and mpk6 mutants in the absence of PAMP treatment showed no major splicing changes. However, in contrast to MPK3 and MPK6, MPK4 was found to be a key regulator of PAMP-induced DAS events as the AS of a number of splicing factors and immunity-related protein kinases is affected, such as the calcium-dependent protein kinase CPK28, the cysteine-rich receptor like kinases CRK13 and CRK29 or the FLS2 co-receptor SERK4/BKK1. Although MPK4 is guarded by SUMM2 and consequently, the mpk4 dwarf and DEG phenotypes are suppressed in mpk4 summ2 mutants, MPK4-dependent DAS is not suppressed by SUMM2, supporting the notion that PAMP-triggered MPK4 activation mediates regulation of alternative splicing.

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied PAMP-induced alternative splicing in Arabidopsis and identified several hundred differentially alternatively spliced (DAS) genes. Importantly, of these PAMP-induced DAS genes, only 18% overlap with the set of PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DAS events. Global DAS analysis of MAPK mutants identified MPK4 as a key regulator of PAMP-induced DAS events. Although MPK4 is guarded by SUMM2 and consequently, the mpk4 dwarf and DEG phenotypes are suppressed in mpk4 summ2 mutants, MPK4-dependent DAS is not suppressed by SUMM2, showing that PAMP-triggered MPK4 activation mediates regulation of alternative splicing.

Distinct Evolutionary Origins of Intron Retention Splicing Events in NHX1 Antiporter Transcripts Relate to Sequence Specific Distinctions in Oryza Species

The genome of Asian cultivated rice (Oryza sativa L.) shows the presence of six organelle-specific and one plasma membrane (OsNHX1-7) NHX-type cation proton antiporters. Of these, vacuolar-localized OsNHX1 is extensively characterized. The genus Oryza consists of 27 species and 11 genome-types, with cultivated rice, diploid O. sativa, having an AA-type genome. Oryza NHX1 orthologous regions (gene organization, 5' upstream cis elements, amino acid residues/motifs) from closely related Oryza AA genomes cluster distinctly from NHX1 regions from more ancestral Oryza BB, FF and KKLL genomes. These sequence-specific distinctions also extend to two separate intron retention (IR) events involving Oryza NHX1 transcripts that occur at the 5' and 3' ends of the NHX1 transcripts. We demonstrate that the IR event involving the 5' UTR is present only in more recently evolved Oryza AA genomes while the IR event governing retention of the 13th intron of Oryza NHX1 (terminal intron) is more ancient in origin, also occurring in halophytic wild rice, Oryza coarctata (KKLL). We also report presence of a retro-copy of the OcNHX1 cDNA in the genome of O. coarctata (rOcNHX1). Preferential species and tissue specific up- or down-regulation of the correctly spliced NHX1 transcript/5' UTR/13th intron-retaining splice variants under salinity was observed. The implications of IR on NHX1 mRNA stability and ORF diversity in Oryza spp. is discussed.