Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize

Leveraging Multiomics Insights and Exploiting Wild Relatives' Potential for Drought and Heat Tolerance in Maize

Climate change, particularly drought and heat stress, may slash agricultural productivity by 25.7% by 2080, with maize being the hardest hit. Therefore, unraveling the molecular nature of plant responses to these stressors is vital for the development of climate-smart maize. This manuscript's primary objective was to examine how maize plants respond to these stresses, both individually and in combination. Additionally, the paper delved into harnessing the potential of maize wild relatives as a valuable genetic resource and leveraging AI-based technologies to boost maize resilience. The role of multiomics approaches particularly genomics and transcriptomics in dissecting the genetic basis of stress tolerance was also highlighted. The way forward was proposed to utilize a bunch of information obtained through omics technologies by an interdisciplinary state-of-the-art forward-looking big-data, cyberagriculture system, and AI-based approach to orchestrate the development of climate resilient maize genotypes.

Alternative splicing responses to salt stress in Glycyrrhiza uralensis revealed by global profiling of transcriptome RNA-seq datasets

Excessive reactive oxygen species stress due to salinity poses a significant threat to the growth of Glycyrrhiza uralensis Fisch. To adapt to salt stress, G. uralensis engages in alternative splicing (AS) to generate a variety of proteins that help it withstand the effects of salt stress. While several studies have investigated the impact of alternative splicing on plants stress responses, the mechanisms by which AS interacts with transcriptional regulation to modulate the salt stress response in G. uralensis remain poorly understood. In this study, we utilized high-throughput RNA sequencing data to perform a comprehensive analysis of AS events at various time points in G. uralensis under salt stress, with exon skipping (SE) being the predominant AS type. KEGG enrichment analysis was performed on the different splicing genes (DSG), and pathways associated with AS were significantly enriched, including RNA transport, mRNA surveillance, and spliceosome. This indicated splicing regulation of genes, resulting in AS events under salt stress conditions. Moreover, plant response to salt stress pathways were also enriched, such as mitogen-activated protein kinase signaling pathway - plant, flavonoid biosynthesis, and oxidative phosphorylation. We focused on four differentially significant genes in the MAPK pathway by AS and qRT-PCR analysis. The alternative splicing type of MPK4 and SnRK2 was skipped exon (SE). ETR2 and RbohD were retained intron (RI) and alternative 5'splice site (A5SS), respectively. The expression levels of isoform1 of these four genes displayed different but significant increases in different tissue sites and salt stress treatment times. These findings suggest that MPK4, SnRK2, ETR2, and RbohD in G. uralensis activate the expression of isoform1, leading to the production of more isoform1 protein and thereby enhancing resistance to salt stress. These findings suggest that salt-responsive AS directly and indirectly governs G. uralensis salt response. Further investigations into AS function and mechanism during abiotic stresses may offer novel references for bolstering plant stress tolerance.

Alternative Splicing Underpins the ALMT9 Transporter Function for Vacuolar Malic Acid Accumulation in Apple

Vacuolar malic acid accumulation largely determines fruit acidity, a key trait for the taste and flavor of apple and other fleshy fruits. Aluminum-activated malate transporter 9 (ALMT9/Ma1) underlies a major genetic locus, Ma, for fruit acidity in apple, but how the protein transports malate across the tonoplast is unclear. Here, it is shown that overexpression of the coding sequence of Ma1 (Ma1α) drastically decreases fruit acidity in "Royal Gala" apple, leading to uncovering alternative splicing underpins Ma1's function. Alternative splicing generates two isoforms: Ma1β is 68 amino acids shorter with much lower expression than the full-length protein Ma1α. Ma1β does not transport malate itself but interacts with the functional Ma1α to form heterodimers, creating synergy with Ma1α for malate transport in a threshold manner (When Ma1β/Ma1α ≥ 1/8). Overexpression of Ma1α triggers feedback inhibition on the native Ma1 expression via transcription factor MYB73, decreasing the Ma1β level well below the threshold that leads to significant reductions in Ma1 function and malic acid accumulation in fruit. Overexpression of Ma1α and Ma1β or genomic Ma1 increases both isoforms proportionally and enhances fruit malic acid accumulation. These findings reveal an essential role of alternative splicing in ALMT9-mediated malate transport underlying apple fruit acidity.

Functional Validation of Different Alternative Splicing Variants of the Chrysanthemum lavandulifolium ClNUM1 Gene in Tobacco

The Asteraceae are widely distributed throughout the world, with diverse functions and large genomes. Many of these genes remain undiscovered and unstudied. In this study, we discovered a new gene ClNUM1 in Chrysanthemum lavandulifolium and studied its function. In this study, bioinformatics, RT-qPCR, paraffin sectioning, and tobacco transgenics were utilized to bioinformatically analyze and functionally study the three variable splice variants of the unknown gene ClNUM1 cloned from C. lavandulifolium. The results showed that ClNUM1.1 and ClNUM1.2 had selective 3' splicing and selective 5' splicing, and ClNUM1.3 had selective 5' splicing. When the corresponding transgenic tobacco plants were subjected to abiotic stress treatment, in the tobacco seedlings, the ClNUM1.1 gene and the ClNUM1.2 gene enhanced salt and low-temperature tolerance and the ClNUM1.3 gene enhanced low-temperature tolerance; in mature tobacco plants, the ClNUM1.1 gene was able to enhance salt and low-temperature tolerance, and the ClNUM1.2 and ClNUM1.3 genes were able to enhance low-temperature tolerance. In summary, there are differences in the functions of the different splice variants and the different seedling stages of transgenic tobacco, but all of them enhanced the resistance of tobacco to a certain extent. The analysis and functional characterization of the ClNUM1 gene provided new potential genes and research directions for abiotic resistance breeding in Chrysanthemum.

Alternative Splicing Plays a Crucial Role in the Salt Tolerance of Foxtail Millet

Foxtail millet is an important cereal crop that is relatively sensitive to salt stress, with its yield significantly affected by such stress. Alternative splicing (AS) widely affects plant growth, development, and adaptability to stressful environments. Through RNA-seq analysis of foxtail millet under different salt treatment periods, 2078 AS events were identified, and analyses were conducted on differential gene (DEG), differential alternative splicing gene (DASG), and overlapping gene. To investigate the regulatory mechanism of AS in response to salt stress in foxtail millet, the foxtail millet AS genes SiCYP19, with two AS variants (SiCYP19-a and SiCYP19-b), were identified and cloned. Yeast overexpression experiments indicated that SiCYP19 may be linked to the response to salt stress. Subsequently, we conducted overexpression experiments of both alternative splicing variants in foxtail millet roots to validate them experimentally. The results showed that, under salt stress, both SiCYP19-a and SiCYP19-b jointly regulated the salt tolerance of foxtail millet. Specifically, overexpression of SiCYP19-b significantly increased the proline content and reduced the accumulation of reactive oxygen species (ROS) in foxtail millet, compared to that in SiCYP19-a. This shows that SiCYP19-b plays an important role in increasing the content of proline and promoting the clearance of ROS, thus improving the salt tolerance of foxtail millet.

The Emerging Role of Non-Coding RNAs (ncRNAs) in Plant Growth, Development, and Stress Response Signaling

Plant species utilize a variety of regulatory mechanisms to ensure sustainable productivity. Within this intricate framework, numerous non-coding RNAs (ncRNAs) play a crucial regulatory role in plant biology, surpassing the essential functions of RNA molecules as messengers, ribosomal, and transfer RNAs. ncRNAs represent an emerging class of regulators, operating directly in the form of small interfering RNAs (siRNAs), microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). These ncRNAs exert control at various levels, including transcription, post-transcription, translation, and epigenetic. Furthermore, they interact with each other, contributing to a variety of biological processes and mechanisms associated with stress resilience. This review primarily concentrates on the recent advancements in plant ncRNAs, delineating their functions in growth and development across various organs such as root, leaf, seed/endosperm, and seed nutrient development. Additionally, this review broadens its scope by examining the role of ncRNAs in response to environmental stresses such as drought, salt, flood, heat, and cold in plants. This compilation offers updated information and insights to guide the characterization of the potential functions of ncRNAs in plant growth, development, and stress resilience in future research.

Determination of Differential Alternative Splicing Under Stress Conditions

Alternative splicing (AS) is an important mechanism contributing to stress-induced regulation of gene expression and proteome diversity. Massive sequencing technologies allow the identification of transcripts generated via stress-responsive AS, potentially important for adaptation to stress conditions. Several bioinformatics tools have been developed to identify differentially expressed alternative splicing events/transcripts from RNA-sequencing results. This chapter describes a detailed protocol for differential alternative splicing analysis using the rMATS tool. In addition, we provide guidelines for validation of the detected splice variants by qRT-PCR based on the obtained output files.

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.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Effects of Climate Change and Drought Tolerance on Maize Growth

Climate change is affecting all regions of the world with different climates, and the scale of damage is increasing due to the occurrence of various natural disasters. In particular, maize production is highly affected by abnormal climate events such as heat waves and droughts. Increasing temperatures can accelerate growth and shorten the growing season, potentially reducing productivity. Additionally, enhanced temperatures during the ripening period can accelerate the process, reducing crop yields. In addition, drought stress due to water deficit can greatly affect seedling formation, early plant growth, photosynthesis, reproductive growth, and yield, so proper water management is critical to maize growth. Maize, in particular, is tall and broad-leaved, so extreme drought stress at planting can cause leaves to curl and stunt growth. It is important to understand that severe drought can have a detrimental effect on the growth and reproduction of maize. In addition, high temperatures caused by drought stress can inhibit the induction of flowering in male flowers and cause factors that interfere with pollen development. It is therefore important to increase the productivity of all food crops, including maize, while maintaining them in the face of persistent drought caused by climate change. This requires a strategy to develop genetically modified crops and drought-tolerant maize that can effectively respond to climate change. The aim of this paper is to investigate the effects of climate change and drought tolerance on maize growth. We also reviewed molecular breeding techniques to develop drought-tolerant maize varieties in response to climate change.

Characterization of the Isocitrate Dehydrogenase Gene Family and Their Response to Drought Stress in Maize

Isocitrate dehydrogenase (IDH) is a key rate-limiting enzyme in the tricarboxylic acid cycle and acts in glutamine synthesis. IDH also participates in plant growth and development and in response to abiotic stresses. We identified 11 maize IDH genes (ZmIDH) and classified these genes into ZmNAD-IDH and ZmNADP-IDH groups based on their different coenzymes (NAD+ or NADP+). The ZmNAD-IDH group was further divided into two subgroups according to their catalytic and non-catalytic subunits, as in Arabidopsis. The ZmIDHs significantly differed in physicochemical properties, gene structure, conserved motifs, and protein tertiary structure. Promoter prediction analysis revealed that the promoters of these ZmIDHs contain cis-acting elements associated with light response, abscisic acid, phytohormones, and abiotic stresses. ZmIDH is predicted to interact with proteins involved in development and stress resistance. Expression analysis of public data revealed that most ZmIDHs are specifically expressed in anthers. Different types of ZmIDHs responded to abiotic stresses with different expression patterns, but all exhibited responses to abiotic stresses to some extent. In addition, analysis of the public sequence from transcription data in an association panel suggested that natural variation in ZmIDH1.4 will be associated with drought tolerance in maize. These results suggested that ZmIDHs respond differently and/or redundantly to abiotic stresses during plant growth and development, and this analysis provides a foundation to understand how ZmIDHs respond to drought stress in maize.

Alternative splicing impacts the rice stripe virus response transcriptome

Alternative splicing (AS) is an important form of post transcriptional modification present in both animals and plants. However, little information was obtained about AS events in response to plant virus infection. In this study, we conducted a genome-wide transcriptome analysis on AS change in rice infected by a devastating virus, Rice stripe virus (RSV). KEGG analysis was performed on the differentially expressed (DE) genes and differentially alternative spliced (DAS) genes. The results showed that DE genes were significantly enriched in the pathway of interaction with plant pathogens. The DAS genes were mainly enriched in basal metabolism and RNA splicing pathways. The heat map clustering showed that DEGs clusters were mainly enriched in regulation of transcription and defense response while differential transcript usage (DTU) clusters were strongly enriched in mRNA splicing and calcium binding. Overall, our results provide a fundamental basis for gene-wide AS changes in rice after RSV infection.

Advances in alternative splicing identification: deep learning and pantranscriptome

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Transcriptome and WGCNA reveal hub genes in sugarcane tiller seedlings in response to drought stress

Drought stress can severely affect sugarcane growth and yield. The objective of this research was to identify candidate genes in sugarcane tillering seedlings in response to drought stress. We performed a comparative phenotypic, physiological and transcriptomic analysis of tiller seedlings of drought-stressed and well-watered "Guire 2" sugarcane, in a time-course experiment (5 days, 9 days and 15 days). Physiological examination reviewed that SOD, proline, soluble sugars, and soluble proteins accumulated in large amounts in tiller seedlings under different intensities of drought stress, while MDA levels remained at a stable level, indicating that the accumulation of osmoregulatory substances and the enhancement of antioxidant enzyme activities helped to limit further damage caused by drought stress. RNA-seq and weighted gene co-expression network analysis (WGCNA) were performed to identify genes and modules associated with sugarcane tillering seedlings in response to drought stress. Drought stress induced huge down-regulated in gene expression profiles, most of down-regulated genes were mainly associated with photosynthesis, sugar metabolism and fatty acid synthesis. We obtained four gene co-expression modules significantly associated with the physiological changes under drought stress (three modules positively correlated, one module negatively correlated), and found that LSG1-2, ERF1-2, SHKA, TIL, HSP18.1, HSP24.1, HSP16.1 and HSFA6A may play essential regulatory roles as hub genes in increasing SOD, Pro, soluble sugar or soluble protein contents. In addition, one module was found mostly involved in tiller stem diameter, among which members of the BHLH148 were important nodes. These results provide new insights into the mechanisms by which sugarcane tillering seedlings respond to drought stress.

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.

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.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses

BACKGROUND

Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions.

RESULTS

We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes.

CONCLUSION

Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

Two environmental signal-driven RNA metabolic processes: Alternative splicing and translation

Plants live in fixed locations and have evolved adaptation mechanisms that integrate multiple responses to various environmental signals. Among the different components of these response pathways, receptors/sensors represent nodes that recognise environmental signals. Additionally, RNA metabolism plays an essential role in the regulation of gene expression and protein synthesis. With the development of RNA biotechnology, recent advances have been made in determining the roles of RNA metabolism in response to different environmental signals-especially the roles of alternative splicing and translation. In this review, we discuss recent progress in research on how the environmental adaptation mechanisms in plants are affected at the posttranscriptional level. These findings improve our understanding of the mechanism through which plants adapt to environmental changes by regulating the posttranscriptional level and are conducive for breeding stress-tolerant plants to cope with dynamic and rapidly changing environments.

A bZIP transcription factor (CiFD) regulates drought- and low-temperature-induced flowering by alternative splicing in citrus

Drought and low temperature are two key environmental factors that induce adult citrus flowering. However, the underlying regulation mechanism is poorly understood. The bZIP transcription factor FD is a key component of the florigen activation complex (FAC) which is composed of FLOWERING LOCUS T (FT), FD, and 14-3-3 proteins. In this study, isolation and characterization of CiFD in citrus found that there was alternative splicing (AS) of CiFD, forming two different proteins (CiFDα and CiFDβ). Further investigation found that their expression patterns were similar in different tissues of citrus, but the subcellular localization and transcriptional activity were different. Overexpression of the CiFD DNA sequence (CiFD-DNA), CiFDα, or CiFDβ in tobacco and citrus showed early flowering, and CiFD-DNA transgenic plants were the earliest, followed by CiFDβ and CiFDα. Interestingly, CiFDα and CiFDβ were induced by low temperature and drought, respectively. Further analysis showed that CiFDα can form a FAC complex with CiFT, Ci14-3-3, and then bind to the citrus APETALA1 (CiAP1) promoter and promote its expression. However, CiFDβ can directly bind to the CiAP1 promoter independently of CiFT and Ci14-3-3. These results showed that CiFDβ can form a more direct and simplified pathway that is independent of the FAC complex to regulate drought-induced flowering through AS. In addition, a bHLH transcription factor (CibHLH96) binds to CiFD promoter and promotes the expression of CiFD under drought condition. Transgenic analysis found that CibHLH96 can promote flowering in transgenic tobacco. These results suggest that CiFD is involved in drought- and low-temperature-induced citrus flowering through different regulatory patterns.

Transcriptome diversity assessment of Gossypium arboreum (FDH228) leaves under control, drought and whitefly infestation using PacBio long reads

Alternative splicing (AS) and alternative polyadenylation (APA) are common mechanisms in eukaryotes to increase the complexity of transcriptomes and subsequently proteomes. Analysis of long reads transcriptomics data can result in the discovery of novel transcripts, splice sites, AS or APA events. Gossypium arboreum is an important cultivated cotton species and a putative contributor of the A sub-genome to the modern tetraploid cotton; and inherently tolerant to several biotic and abiotic stresses. Specifically, its variety 'FDH228' is considered to be an important resistance source. In this study, we sequenced the G. arboreum (var. FDH228) transcriptome using PacBio IsoSeq and illumina short read sequencing under three different conditions i.e. untreated/healthy, treated with biotic stress through whitefly infestation, and treated with abiotic stress via water deprivation, for the discovery and surveying of canonical and non-canonical AS, APA and transcript fusion events. We were able to obtain 15,419 unique transcripts from all samples representing 11,343 genes, out of which 10,832 were annotated and 520 were novel with respect to the published reference genome. These transcripts were grouped into different structural categories including 60 Antisense, 11,959 having a full-splice match, 999 with incomplete-splice match, 30 fusion transcripts, 177 genic, 479 intergenic, 771 novels in the catalog, and 944 Novel but not found in the catalog. Subsequently, randomly selected candidate transcripts were experimentally validated using qRT-PCR. Our comprehensive identification of canonical and non-canonical splicing events, and novel and fusion transcripts aids in the understanding of the resistance mechanisms for this specific germplasm.

The battle of crops against drought: Genetic dissection and improvement

With ongoing global climate change, water scarcity-induced drought stress remains a major threat to agricultural productivity. Plants undergo a series of physiological and morphological changes to cope with drought stress, including stomatal closure to reduce transpiration and changes in root architecture to optimize water uptake. Combined phenotypic and multi-omics studies have recently identified a number of drought-related genetic resources in different crop species. The functional dissection of these genes using molecular techniques has enriched our understanding of drought responses in crops and has provided genetic targets for enhancing resistance to drought. Here, we review recent advances in the cloning and functional analysis of drought resistance genes and the development of technologies to mitigate the threat of drought to crop production.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Global Landscape of Alternative Splicing in Maize Response to Low Temperature

Maize (Zea mays L.) is an important food crop planted across the world, and low-temperature stress can affect maize germination. Alternative splicing (AS) is widely present in plants under abiotic stress; however, the response of AS to low-temperature stress in maize remains unclear. In this study, a genome-wide analysis of AS during maize response to low temperatures was performed. AS events were distributed on each chromosome, approximately 2.05-2.09 AS events per gene. Seven genes only had AS in low-temperature-resistant inbred lines. A total of 278 KEGGs and 46 GOs were enriched based on overlapping AS genes, which were associated with hormone and oxidoreductase activity. The mutant was used to verify the function of AS gene ZmWRKY48, and the RGR, RSL, RRL, and RRSA of the mutant decreased by 15.16%-19.87% compared with the normal line. These results contribute to subsequent analysis of the regulatory mechanism of maize in response to low-temperature stress.

Transcriptome Analysis Reveals the Involvement of Alternative Splicing in the Nitrogen Starvation Response of Chlamydomonas reinhardtii

Alternative splicing (AS) is a regulatory mechanism of post-transcriptional regulation that plays an important role in plant response to abiotic stresses. However, corresponding research involving the mechanism of AS in the nitrogen starvation response of C. reinhardtii is rare. This study performed a comprehensive and systematic analysis of AS events in C. reinhardtii at nine time points (0 h, 10 m, 30 m, 1 h, 6 h, 8 h, 24 h, and 48 h) under nitrogen starvation. It used STAR and rMATS tools to identify and quantify the probability of the AS event happening through the transcriptome high-throughput sequencing data. A total of 5806 AS events in 3500 genes were identified, and the retained intron and skipped exon were considered the main AS types. The genes related to the AS event in nitrogen starvation were mainly involved in spliceosome and transporter and enriched in the citrate cycle and fatty acid degradation pathways. These results suggested that AS may play an important role in the nitrogen starvation response in C. reinhardtii, and provided insights into post-transcriptional regulation under nitrogen starvation.

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.

Time Course RNA-seq Reveals Soybean Responses against Root-Lesion Nematode and Resistance Players

Pratylenchus brachyurus causes serious damage to soybean production and other crops worldwide. Plant molecular responses to RLN infection remain largely unknown and no resistance genes have been identified in soybean. In this study, we analyzed molecular responses to RLN infection in moderately resistant BRSGO (Chapadões-BRS) and susceptible TMG115 RR (TMG) Glycine max genotypes. Differential expression analysis revealed two stages of response to RLN infection and a set of differentially expressed genes (DEGs) in the first stage suggested a pattern-triggered immunity (PTI) in both genotypes. The divergent time-point of DEGs between genotypes was observed four days post-infection, which included the activation of mitogen-activated protein kinase (MAPK) and plant-pathogen interaction genes in the BRS, suggesting the occurrence of an effector-triggered immunity response (ETI) in BRS. The co-expression analyses combined with single nucleotide polymorphism (SNP) uncovered a key element, a transcription factor phytochrome-interacting factor (PIF7) that is a potential regulator of moderate resistance to RLN infection. Two genes for resistance-related leucine-rich repeat (LRR) proteins were found as BRS-specific expressed genes. In addition, alternative splicing analysis revealed an intron retention in a myo-inositol oxygenase (MIOX) transcript, a gene related to susceptibility, may cause a loss of function in BRS.

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.

Paclobutrazol improves surface water use efficiency by regulating allometric trait behavior in maize

Paclobutrazol (PBZ) role in drought management of maize is least understood. In maize, root traits are linked with surface water management. Over three years, early and terminal deficit irrigation (EDI and TDI) with or without PBZ were imposed on DKC-9144 and PG-2475 maize varieties. Several allometric parameters viz. stem height, stem diameter, leaf area and root traits along with physiological processes were measured. Implication of these parameters in the management of soil surface irrigation in terms of water use efficiency (WUE) was demonstrated in maize. Increased number of lateral roots and root number density in DKC-9144 provided more surface area for water absorption for better management of EDI. Root growth rates showed a similar pattern with root length, root surface areas, and root numbers in EDI. Elevated expressions of ZmRTCL, ZmRTCS and ZmARF34 in EDI and EDI plus PBZ were associated with seminal roots and root laterals initiation. Under TDI alone or in combination with PBZ, root lengths (BRL, CRL, SRL) and root surface areas varied in DKC-9144 and PG-2475 over control. Furthermore, correlation analysis showed that decrease in WUE under TDI was significantly associated with a reduction in stem thickness and leaf surface area. For WUE_N in TDI and PBZ plus TDI, structural equation modelling proposed, brace root surface area (BRSA_N) as a positive contributor, while a negative contributor was seminal root surface area (SRSA_N). Present study explained the importance of specific root traits and their association with other allometric parameters for improving WUE in DKC-9144 variety of maize and the crop in general.

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.

A crosstalk of circadian clock and alternative splicing under abiotic stresses in the plants

The circadian clock is an internal time-keeping mechanism that synchronizes the physiological adaptation of an organism to its surroundings based on day and night transition in a period of 24 h, suggesting the circadian clock provides fitness by adjusting environmental constrains. The circadian clock is driven by positive and negative elements that regulate transcriptionally and post-transcriptionally. Alternative splicing (AS) is a crucial transcriptional regulator capable of generating large numbers of mRNA transcripts from limited numbers of genes, leading to proteome diversity, which is involved in circadian to deal with abiotic stresses. Over the past decade, AS and circadian control have been suggested to coordinately regulate plant performance under fluctuating environmental conditions. However, only a few reports have reported the regulatory mechanism of this complex crosstalk. Based on the emerging evidence, this review elaborates on the existing links between circadian and AS in response to abiotic stresses, suggesting an uncovered regulatory network among circadian, AS, and abiotic stresses. Therefore, the rhythmically expressed splicing factors and core clock oscillators fill the role of temporal regulators participating in improving plant growth, development, and increasing plant tolerance against abiotic stresses.

Recent developments in multi-omics and breeding strategies for abiotic stress tolerance in maize (Zea mays L.)

High-throughput sequencing technologies (HSTs) have revolutionized crop breeding. The advent of these technologies has enabled the identification of beneficial quantitative trait loci (QTL), genes, and alleles for crop improvement. Climate change have made a significant effect on the global maize yield. To date, the well-known omic approaches such as genomics, transcriptomics, proteomics, and metabolomics are being incorporated in maize breeding studies. These approaches have identified novel biological markers that are being utilized for maize improvement against various abiotic stresses. This review discusses the current information on the morpho-physiological and molecular mechanism of abiotic stress tolerance in maize. The utilization of omics approaches to improve abiotic stress tolerance in maize is highlighted. As compared to single approach, the integration of multi-omics offers a great potential in addressing the challenges of abiotic stresses of maize productivity.

Transcriptomic and alternative splicing analyses provide insights into the roles of exogenous salicylic acid ameliorating waxy maize seedling growth under heat stress

BACKGROUND

Salicylic acid (SA) is a phytohormone which works to regulate the abiotic stress response of plants. However, the molecular mechanism by which SA mediates heat tolerance in waxy maize (Zea mays L. sinsensis Kulesh) remains unknown.

RESULTS

Two varieties of waxy maize seedlings, heat-tolerant 'Yunuo7' (Y7) and heat-sensitive 'Suyunuo5' (S5), were pretreated with SA prior to heat stress (HTS). After treatment, physiological and transcriptomic changes were analyzed. Compared with HTS, the exogenous application of SA enhanced the shoot dry weight, the activities of antioxidant enzymes (e.g., SOD, POD, CAT and APX), and the concentration of endogenous phytohormones (e.g., SA, ABA, IAA, GA3), while decreased the MDA content. Transcriptome analysis showed that the number of differentially expressed genes (DEGs) identified in the control (CK) vs HTS and HTS vs HTS + SA comparisons were more in S5 than in Y7. HTS induced the downregulation of genes involved in photosynthesis and the upregulation of genes encoding heat shock transcription factors (HSFs) and heat shock proteins (HSPs). Compared with HTS, SA pretreatment reversed the expression of 5 photosynthesis-related genes, 26 phytohormone-related genes, and all genes encoding HSFs and HSPs in S5. Furthermore, the number of alternative splicing (AS) events increased under HTS treatment for both varieties, while decreased under SA pretreatment of S5. Differentially spliced genes (DSGs) showed little overlap with DEGs, and DEGs and DSGs differed significantly in functional enrichment.

CONCLUSIONS

Physiological and transcriptional together indicated that HTS and SA pretreatment had a greater effect on S5 than Y7. Additionally, it appears that transcriptional regulation and AS work synergistically to enhance thermotolerance in heat-sensitive waxy maize. Our study revealed the regulatory effects and underlying molecular mechanisms of SA on waxy maize seedling under HTS.

Integrative Proteome and Phosphoproteome Profiling of Early Cold Response in Maize Seedlings

Cold limits the growth and yield of maize in temperate regions, but the molecular mechanism of cold adaptation remains largely unexplored in maize. To identify early molecular events during cold shock, maize seedlings were treated under 4 °C for 30 min and 2 h, and analyzed at both the proteome and phosphoproteome levels. Over 8500 proteins and 19,300 phosphopeptides were quantified. About 660 and 620 proteins were cold responsive at protein abundance or site-specific phosphorylation levels, but only 65 proteins were shared between them. Functional enrichment analysis of cold-responsive proteins and phosphoproteins revealed that early cold response in maize is associated with photosynthesis light reaction, spliceosome, endocytosis, and defense response, consistent with similar studies in Arabidopsis. Thirty-two photosynthesis proteins were down-regulated at protein levels, and 48 spliceosome proteins were altered at site-specific phosphorylation levels. Thirty-one kinases and 33 transcriptional factors were cold responsive at protein, phosphopeptide, or site-specific phosphorylation levels. Our results showed that maize seedlings respond to cold shock rapidly, at both the proteome and phosphoproteome levels. This study provides a comprehensive landscape at the cold-responsive proteome and phosphoproteome in maize seedlings that can be a significant resource to understand how C₄ plants respond to a sudden temperature drop.

The Function of DNA Demethylase Gene ROS1a Null Mutant on Seed Development in Rice (Oryza Sativa) Using the CRISPR/CAS9 System

The endosperm is the main nutrient source in cereals for humans, as it is a highly specialized storage organ for starch, lipids, and proteins, and plays an essential role in seed growth and development. Active DNA demethylation regulates plant developmental processes and is ensured by cytosine methylation (5-meC) DNA glycosylase enzymes. To find out the role of OsROS1a in seed development, the null mutant of OsROS1a was generated using the CRISPR/Cas9 system. The null mutant of OsROS1a was stable and heritable, which affects the major agronomic traits, particularly in rice seeds. The null mutant of OsROS1a showed longer and narrower grains, and seeds were deformed containing an underdeveloped and less-starch-producing endosperm with slightly irregularly shaped embryos. In contrast to the transparent grains of the wild type, the grains of the null mutant of OsROS1a were slightly opaque and rounded starch granules, with uneven shapes, sizes, and surfaces. A total of 723 differential expression genes (DEGs) were detected in the null mutant of OsROS1a by RNA-Seq, of which 290 were downregulated and 433 were upregulated. The gene ontology (GO) terms with the top 20 enrichment factors were visualized for cellular components, biological processes, and molecular functions. The key genes that are enriched for these GO terms include starch synthesis genes (OsSSIIa and OsSSIIIa) and cellulose synthesis genes (CESA2, CESA3, CESA6, and CESA8). Genes encoding polysaccharides and glutelin were found to be downregulated in the mutant endosperm. The glutelins were further verified by SDS-PAGE, suggesting that glutelin genes could be involved in the null mutant of OsROS1a seed phenotype and OsROS1a could have the key role in the regulation of glutelins. Furthermore, 378 differentially alternative splicing (AS) genes were identified in the null mutant of OsROS1a, suggesting that the OsROS1a gene has an impact on AS events. Our findings indicated that the function on rice endosperm development in the null mutant of OsROS1a could be influenced through regulating gene expression and AS, which could provide the base to properly understand the molecular mechanism related to the OsROS1a gene in the regulation of rice seed development.

Analysis of Whole-Transcriptome RNA-Seq Data Reveals the Involvement of Alternative Splicing in the Drought Response of Glycyrrhiza uralensis

Alternative splicing (AS) is a post-transcriptional regulatory mechanism that increases protein diversity. There is growing evidence that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate the drought response in Glycyrrhiza uralensis remains unclear. In this study, we performed a genome-wide analysis of AS events in G. uralensis at different time points under drought stress using a high-throughput RNA sequencing approach. We detected 2,479 and 2,764 AS events in the aerial parts (AP) and underground parts (UP), respectively, of drought-stressed G. uralensis. Of these, last exon AS and exon skipping were the main types of AS. Overall, 2,653 genes undergoing significant AS regulation were identified from the AP and UP of G. uralensis exposed to drought for 2, 6, 12, and 24 h. Gene Ontology analyses indicated that AS plays an important role in the regulation of nitrogen and protein metabolism in the drought response of G. uralensis. Notably, the spliceosomal pathway and basal transcription factor pathway were significantly enriched with differentially spliced genes under drought stress. Genes related to splicing regulators in the AP and UP of G. uralensis responded to drought stress and underwent AS under drought conditions. In summary, our data suggest that drought-responsive AS directly and indirectly regulates the drought response of G. uralensis. Further in-depth studies on the functions and mechanisms of AS during abiotic stresses will provide new strategies for improving plant stress resistance.

Cloning and homologous characterization of geranylgeranyl pyrophosphate synthase (GGPPS) from Withania somnifera revealed alterations in metabolic flux towards gibberellic acid biosynthesis

Overexpression of a novel geranylgeranyl pyrophosphate synthase gene (WsGGPPS) in planta resulted in increased levels of gibberellic acid and decrease in withanolide content. Withania somnifera (L.) Dunal, the herb from family Solanaceae is one of the most treasured medicinal plant used in traditional medicinal systems owing to its unique stockpile of pharmaceutically active secondary metabolites. Phytochemical and pharmacological studies in this plant were well established, but the genes affecting the regulation of biosynthesis of major metabolites were not well elucidated. In this study cloning and functional characterization of a key enzyme in terpenoid biosynthetic pathway viz. geranylgeranyl pyrophosphate synthase (EC 2.5.1.29) gene from Withania somnifera was performed. The full length WsGGPPS gene contained 1,104 base pairs that encode a polypeptide of 365 amino acids. The quantitative expression analysis suggested that WsGGPPS transcripts were expressed maximally in flower tissues followed by berry tissues. The expression levels of WsGGPPS were found to be regulated by methyl jasmonate (MeJA) and salicylic acid (SA). Amino acid sequence alignment and phylogenetic studies suggested that WsGGPPS had close similarities with GGPPS of Solanum tuberosum and Solanum pennellii. The structural analysis provided basic information about three dimensional features and physicochemical parameters of WsGGPPS protein. Overexpression of WsGGPPS in planta for its functional characterization suggested that the WsGGPPS was involved in gibberellic acid biosynthesis.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the NAM pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

The Importance of a Genome-Wide Association Analysis in the Study of Alternative Splicing Mutations in Plants with a Special Focus on Maize

Alternative splicing is an important mechanism for regulating gene expressions at the post-transcriptional level. In eukaryotes, the genes are transcribed in the nucleus to produce pre-mRNAs and alternative splicing can splice a pre-mRNA to eventually form multiple different mature mRNAs, greatly increasing the number of genes and protein diversity. Alternative splicing is involved in the regulation of various plant life activities, especially the response of plants to abiotic stresses and is also an important process of plant growth and development. This review aims to clarify the usefulness of a genome-wide association analysis in the study of alternatively spliced variants by summarizing the application of alternative splicing, genome-wide association analyses and genome-wide association analyses in alternative splicing, as well as summarizing the related research progress.

High transcriptome plasticity drives phosphate starvation responses in tomato

Tomato is an important vegetable crop and fluctuating available soil phosphate (Pi) level elicits several morpho-physiological responses driven by underlying molecular responses. Therefore, understanding these molecular responses at the gene and isoform levels has become critical in the quest for developing crops with improved Pi use efficiency. A quantitative time-series RNA-seq analysis was performed to decipher the global transcriptomic changes that accompany Pi starvation in tomato. Apart from changes in the expression levels of genes, there were also alterations in the expression of alternatively-spliced transcripts. Physiological responses such as anthocyanin accumulation, reactive oxygen species generation and cell death are obvious 7 days after Pi deprivation accompanied with the maximum amount of transcriptional change in the genome making it an important stage for in-depth study while studying Pi stress responses (PSR). Our study demonstrates that transcriptomic changes under Pi deficiency are dynamic and complex in tomato. Overall, our study dwells on the dynamism of the transcriptome in eliciting a response to adapt to low Pi stress and lays it bare. Findings from this study will prove to be an invaluable resource for researchers using tomato as a model for understanding nutrient deficiency.

Rapid Regulation of Alternative Splicing in Response to Environmental Stresses

Plants overcome the changing environmental conditions through diverse strategies and complex regulations. In addition to direct regulation of gene transcription, alternative splicing (AS) also acts as a crucial regulatory mechanism to cope with various stresses. Generating from the same pre-mRNA, AS events allow rapid adjustment of the abundance and function of key stress-response components. Mounting evidence has indicated the close link between AS and plant stress response. However, the mechanisms on how environmental stresses trigger AS are far from understood. The advancing high-throughput sequencing technologies have been providing useful information, whereas genetic approaches have also yielded remarkable phenotypic evidence for AS control of stress responses. It is important to study how stresses trigger AS events for both fundamental science and applications. We review current understanding of stress-responsive AS in plants and discuss research challenges for the near future, including regulation of splicing factors, epigenetic modifications, the shared targets of splice isoforms, and the stress-adjusting ratios between splicing variants.

Genome-Wide Profiling of Alternative Splicing and Gene Fusion during Rice Black-Streaked Dwarf Virus Stress in Maize (Zea mays L.)

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease (MRDD), which is a viral disease that significantly affects maize yields worldwide. Plants tolerate stress through transcriptional reprogramming at the alternative splicing (AS), transcriptional, and fusion gene (FG) levels. However, it is unclear whether and how AS and FG interfere with transcriptional reprogramming in MRDD. In this study, we performed global profiling of AS and FG on maize response to RBSDV and compared it with transcriptional changes. There are approximately 1.43 to 2.25 AS events per gene in maize infected with RBSDV. GRMZM2G438622 was only detected in four AS modes (A3SS, A5SS, RI, and SE), whereas GRMZM2G059392 showed downregulated expression and four AS events. A total of 106 and 176 FGs were detected at two time points, respectively, including six differentially expressed genes and five differentially spliced genes. The gene GRMZM2G076798 was the only FG that occurred at two time points and was involved in two FG events. Among these, 104 GOs were enriched, indicating that nodulin-, disease resistance-, and chloroplastic-related genes respond to RBSDV stress in maize. These results provide new insights into the mechanisms underlying post-transcriptional and transcriptional regulation of maize response to RBSDV stress.

Analysis of gene expression in early seed germination of rice: landscape and genetic regulation

BACKGROUND

Seed germination is a crucial process, which determines the initiation of seed plant life cycle. The early events during this important life cycle transition that called early seed germination is defined as initially water uptake plus radicle growing out of the covering seed layers. However, a specific genome-wide analysis of early seed germination in rice is still obscure.

RESULTS

In this study, the physiological characteristics of rice seed during seed germination are determined to define key points of early seed germination. Transcriptome analyses of early phase of seed germination provided deeper insight into the genetic regulation landscape. Many genes involved in starch-to-sucrose transition were differentially expressed, especially alpha-amylase 1b and beta-amylase 2, which were predominantly expressed. Differential exon usage (DEU) genes were identified, which were significantly enriched in the pathway of starch and sucrose metabolism, indicating that DEU events were critical for starch-to-sucrose transition at early seed germination. Transcription factors (TFs) were also dramatic expressed, including the abscisic acid (ABA) responsive gene, OsABI5, and gibberellic acid (GA) responsive genes, GAI. Moreover, GAI transactivated GA responsive gene, GAMYB in vivo, indicating a potential pathway involved in early seed germination process. In addition, CBL-interacting protein kinase (CIPK) genes, such as CIPK13, CIPK14 and CIPK17 were potentially interacted with other proteins, indicating its pivotal role at early seed germination.

CONCLUSION

Taken together, gene regulation of early seed germination in rice was complex and protein-to-gene or protein-to-protein interactions were indispensable.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Drought Response in Wheat

RNA/DNA difference (RDD) is a post-transcriptional RNA modification to enrich genetic information, widely involved in regulating diverse biological processes in eukaryotes. RDDs in the wheat nuclear genome, especially those associated with drought response or tolerance, were not well studied up to now. In this study, we investigated the RDDs related to drought response based on the RNA-seq data of drought-stressed and control samples in wheat. In total, 21,782 unique RDDs were identified, of which 265 were found to be drought-induced, representing the first drought-responsive RDD landscape in the wheat nuclear genome. The drought-responsive RDDs were located in 69 genes, of which 35 were differentially expressed under drought stress. Furthermore, the effects of RNA/DNA differences were investigated, showing that they could result in changes of RNA secondary structure, miRNA-target binding as well as protein conserved domains in the RDD-containing genes. In particular, the A to C mutation in TraesCS2A02G053100 (orthology to OsRLCK) led to the loss of tae-miR9657b-5p targeting, indicating that RNA/DNA difference might mediate miRNA to regulate the drought-response process. This study reported the first drought-responsive RDDs in the wheat nuclear genome. It sheds light on the roles of RDD in drought tolerance, and may also contribute to wheat genetic improvement based on epi-transcriptome methods.

Alpha-Linolenic Acid Mediates Diverse Drought Responses in Maize (Zea mays L.) at Seedling and Flowering Stages

Water shortage caused by long-term drought is one of the most serious abiotic stress factors in maize. Different drought conditions lead to differences in growth, development, and metabolism of maize. In previous studies, proteomics and genomics methods have been widely used to explain the response mechanism of maize to long-term drought, but there are only a few articles related to metabolomics. In this study, we used transcriptome and metabolomics analysis to characterize the differential effects of drought stress imposed at seedling or flowering stages on maize. Through the association analysis of genes and metabolites, we found that maize leaves had 61 and 54 enriched pathways under seedling drought and flowering drought, respectively, of which 13 and 11 were significant key pathways, mostly related to the biosynthesis of flavonoids and phenylpropanes, glutathione metabolism and purine metabolism. Interestingly, we found that the α-linolenic acid metabolic pathway differed significantly between the two treatments, and a total of 10 differentially expressed genes and five differentially abundant metabolites have been identified in this pathway. Some differential accumulation of metabolites (DAMs) was related to synthesis of jasmonic acid, which may be one of the key pathways underpinning maize response to different types of long-term drought. In general, metabolomics provides a new method for the study of water stress in maize and lays a theoretical foundation for drought-resistant cultivation of silage maize.

Alternative Splicing Dynamics During the Lifecycle of Salvia miltiorrhiza Root Revealed the Fine Tuning in Root Development and Ingredients Biosynthesis

Alternative splicing (AS) is an essential post-transcriptional process that enhances the coding and regulatory potential of the genome, thereby strongly influencing multiple plant physiology processes, such as metabolic biosynthesis. To explore how AS affects the root development and synthesis of tanshinones and phenolic acid pathways in Salvia miltiorrhiza roots, we investigated the dynamic landscape of AS events in S. miltiorrhiza roots during an annual life history. Temporal profiling represented a distinct temporal variation of AS during the entire development stages, showing the most abundant AS events at the early seedling stage (ES stage) and troughs in 45 days after germination (DAG) and 120 DAG. Gene ontology (GO) analysis indicated that physiological and molecular events, such as lateral root formation, gravity response, RNA splicing regulation, and mitogen-activated protein kinase (MAPK) cascade, were greatly affected by AS at the ES stage. AS events were identified in the tanshinones and phenolic acids pathways as well, especially for the genes for the branch points of the pathways as SmRAS and SmKSL1. Fifteen Ser/Arg-rich (SR) proteins and eight phosphokinases (PKs) were identified with high transcription levels at the ES stage, showing their regulatory roles for the high frequency of AS in this stage. Simultaneously, a co-expression network that includes 521 highly expressed AS genes, SRs, and PKs, provides deeper insight into the mechanism for the variable programming of AS.

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.

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

The use of alternative transcription start or termination sites (aTSS or aTTS) as well as alternative splicing (AS) produce diverse transcript isoforms, playing indispensable roles in the plant development and environmental adaptations. Despite the advances in the finding of the genome-wide alternatively spliced genes in strawberry, it remains unexplored how AS responds to the developmental cues and what relevance do these outcomes have to the gene function. In this study, we have systematically investigated the transcriptome complexity using long-read Oxford Nanopore Technologies along the four successive developmental stages. The full-length cDNA sequencing results unraveled thousands of previously unexplored transcript isoforms raised from aTSS, aTTS, and AS. The relative contributions of these three processes to the complexity of strawberry fruit transcripts were compared. The aTSS and aTTS were more abundant than the AS. Differentially expressed transcripts unraveled the key transitional role of the white fruit stage. Isoform switches of transcripts from 757 genes were observed. They were associated with protein-coding potential change and domain gain or loss as the main consequences. Those genes with switched isoforms take part in the key processes of maturation in the late stages. A case study using yeast two hybrid analysis supported the functional divergence of the two isoforms of the B-box protein 22. Our results provided a new comprehensive overview of the dynamic transcriptomic landscape during strawberry fruit development and maturation.

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.

Increasing yield on dry fields: molecular pathways with growing potential

Drought stress constitutes one of the major constraints to agriculture all over the world, and its devastating effect is only expected to increase in the following years due to climate change. Concurrently, the increasing food demand in a steadily growing population requires a proportional increase in yield and crop production. In the past, research aimed to increase plant resilience to severe drought stress. However, this often resulted in stunted growth and reduced yield under favorable conditions or moderate drought. Nowadays, drought tolerance research aims to maintain plant growth and yield under drought conditions. Overall, recently deployed strategies to engineer drought tolerance in the lab can be classified into a 'growth-centered' strategy, which focuses on keeping growth unaffected by the drought stress, and a 'drought resilience without growth penalty' strategy, in which the main aim is still to boost drought resilience, while limiting the side effects on plant growth. In this review, we put the scope on these two strategies and some molecular players that were successfully engineered to generate drought-tolerant plants: abscisic acid, brassinosteroids, cytokinins, ethylene, ROS scavenging genes, strigolactones, and aquaporins. We discuss how these pathways participate in growth and stress response regulation under drought. Finally, we present an overview of the current insights and future perspectives in the development of new strategies to improve drought tolerance in the field.