Global profiling of alternative splicing landscape responsive to drought, heat and their combination in wheat (Triticum aestivum L.)

Alternative splicing analysis of stress tolerance to Al and flg22 in Vitis quinquangularis

MAIN CONCLUSION

Alternative splicing of transcriptomes after Al and flg22 treatment for 12 h in response to plant defense of Chinese wild Vitis quinquangularis: genes related to stress resistance and splicing factors were identified in response to Al and flg22 treatment. Alternative splicing (AS) is one of the major post-transcriptional regulation processes that potentially regulates the response to biotic and abiotic stresses in plants. So far, the insight into potential roles of AS in grapevine response to aluminium (Al) and flagellin 22 (flg22) stresses remains poorly understood. We performed transcriptome sequencing of grape leaves before and after Al treatment and flg22 treatment, respectively, to identify AS genes. In this study, a total of 11,805 AS events were identified in Al treatment, of which the skipped exon (SE; 88.72%) type was the most frequent. 9156 AS events were identified under flg22 treatment, of which the SE (88.52%) type was the most frequent. Compared with Al-treated and flg22-treated 0 h, there were 42 and 147 differential alternative splicing (DAS) genes differentially expressed (DASEGs) in Al-treated and flg22-treated 12 h, respectively. Functional analysis showed that DASEGs after Al treatment were mainly enriched in glutathione metabolism pathway; DASEGs after flg22 treatment were enriched in MAPK signaling and plant hormone signal transduction. We further verified seven resistance-related DASEGs with up-regulated expression in Al-treated 12 h, including beta-glucosidase, calcineurin B-like protein, synaptotagmin-3, cysteine synthase and glutathione reductase. Several genes function as leucine-rich repeats receptor-like serine/threonine protein kinase, BRI1 associated receptor kinase 1 and receptor-like protein kinase were also verified by RT-qPCR. We also verified four serine/arginine (SR)-rich proteins SCL30A, SCL28, RS2Z32 and SR45A, which were up-regulated in both Al and flg22 stresses. In conclusion, this study provides an in-depth analysis of the correlation between alternative splicing and grapevine stress tolerance, which helps to identify potential candidate genes for useful traits, provides a theoretical basis for grapevine breeding in plant stress tolerance, and offers new perspectives for understanding grapevine environmental adaptation strategies.

TaPP2C-a5 fine-tunes wheat seed dormancy and germination with a Triticeae-specific, alternatively spliced transcript

INTRODUCTION

The sessile plants often experience environmental conditions not ideal for growth, and therefore have evolved strategies to survive and adapt to stress conditions. Abscisic acid (ABA) regulates plant development and abiotic stress response. Clade A type 2C protein phosphatases (PP2Cs), act as co-receptors of ABA, negatively regulate ABA signalling. However, the biological function and detailed molecular mechanism of clade A PP2Cs in ABA signalling pathway remain to be elucidated in wheat.

OBJECTIVES

To analyze the mechanisms of stress response and development mediated by ABA signal precisely regulated by TaPP2C-a5 at the post-transcriptional level in wheat, providing candidate genes for wheat improvement.

METHODS

Based on our previous results of TaPP2Cs gene family analysis, the function and detailed regulation mechanisms of TaPP2C-a5 gene in seed dormancy and germination as well as drought response mediated by ABA signaling pathway were explored through reverse genetics technology.

RESULTS

We found that class A TaPP2C-a5 underwent alternative splicing (AS) to produce two transcripts encoding TaPP2C-a5.1 and TaPP2C-a5.2, respectively. Both TaPP2C-a5.1 and TaPP2C-a5.2 were highly expressed in mature seeds, and were upregulated by exogenous ABA in seedlings. Overexpression of TaPP2C-a5.1 and TaPP2C-a5.2 coordinately negatively regulated seed dormancy and ABA-mediated seed germination as well as post-germination developmental arrest in wheat. TaPP2C-a5.1 negatively regulated drought stress response, while TaPP2C-a5.2 did not participate in drought stress response. The homologous genes of TaPP2C-a5 underwent the same AS as TaPP2C-a5 in tetraploid wheat, but not in rice.

CONCLUSION

Our results revealed that TaPP2C-a5 gene underwent AS and was involved in the regulation of seed dormancy and germination, as well as drought stress response mediated by the ABA signaling at the post-transcriptional level. Our work not only provide a potential target gene to improve PHS resistance, but also emphasize alternative splicing as a strategy with evolution contexts to fine-tune ABA signaling and its involvement in certain biological process.

The Regulatory Roles of RNA-Binding Proteins in Plant Salt Stress Response

Salt stress is one of the most prominent abiotic stresses. Behind the intricate adaptive responses of plants to salt stress, the regulation of gene expression assumes a pivotal role. Complementing transcriptional mechanisms, post-transcriptional regulation performed by RNA-binding proteins provides an additional layer of control through sophisticated molecular machinery. RBPs interact with both RNA molecules and protein partners to coordinate RNA metabolism and, thus, fine-tune the expression of salt-responsive genes, enabling plants to rapidly adapt to ionic challenges. This review systematically evaluates the functional roles of RBPs localized in distinct subcellular compartments, including nuclear, cytoplasmic, chloroplastic, and mitochondrial systems, in mediating post-transcriptional regulatory networks under salinity challenges. Specific classes of RBPs are discussed in detail, including glycine-rich RNA-binding proteins (GR-RBPs), serine/arginine-rich splicing factors (SR proteins), zinc finger domain-containing proteins, DEAD-box RNA helicases (DBRHs), KH domain-containing proteins, Pumilio domain-containing proteins (PUMs), pentatricopeptide repeat proteins (PPRs), and RBPs involved in cytoplasmic RNA granule formation. By integrating their subcellular localization and current mechanistic insights, this review concludes by summarizing the current knowledge and highlighting potential future research directions, aiming to inspire further investigations into the complex network of RBPs in modulating plant responses to salt stress and facilitating the development of strategies to enhance plant salt tolerance.

The alternative splicing of HvLHCA4.2 enhances drought tolerance in barley by regulating ROS scavenging and stomatal closure

Alternative splicing (AS) provides flexible strategies for plants to cope with abiotic stress, but how AS regulates drought response in barley (Hordeum vulgare) remains poorly understood. This study investigates the genome-wide AS patterns in drought-sensitive (Baudin) and tolerant barley genotypes (Tadmor, EC_S1) during drought and recovery stages. Baudin showed maximal AS events during drought, while tolerant genotypes activated AS predominantly during recovery. Drought-tolerant EC_S1 exhibited differential AS events enriched in metabolic pathways, amino acid biosynthesis, and chromatin remodeling. WGCNA analysis identified eight modules linked to stomatal dynamics, biomass, and ROS scavenging, with the 'MEmidnightblue' module (related to stoma and ROS regulation) being drought-responsive. A drought-inducible splice variant, LHCA4.2b, was exclusively identified in tolerant genotypes. Silencing HvLHCA4.2a (lhca4.2a) and dual-silencing HvLHCA4.2a + b (lhca4.2a + b) drastically impaired drought tolerance, manifesting as wilting, accelerated water loss, reduced biomass, and elevated electrolyte leakage, coupled with chlorophyll degradation, ROS overaccumulation, and malondialdehyde overproduction. Enhanced stress sensitivity in lhca4.2a + b versus lhca4.2a underscores HvLHCA4.2b's critical role. Additionally, lhca4.2a + b displayed ABA-insensitive stomata with unaltered stomatal conductance under drought or exogenous ABA, implicating LHCA4.2b in ABA signaling. Mechanistically, AS-generated LHCA4.2b isoform enhances PSI-LHCI super-complex stability through strengthened interaction interfaces with LHCA1, thereby improving energy transfer efficiency and reducing ROS generation. This isoform simultaneously coordinates ABA signaling by elevating ABF1/ABF3 transcription factors while suppressing ABI3/ABI4/ABI5 repressors, collectively modulating stomatal closure and ROS homeostasis. Our findings elucidate AS-mediated drought adaptation mechanisms in barley and highlight HvLHCA4.2b as a potential target for breeding drought-tolerant cultivars.

Adaptation of High-Altitude Plants to Plateau Abiotic Stresses: A Case Study of the Qinghai-Tibet Plateau

High-altitude regions offer outstanding opportunities for investigating the impacts of combined abiotic stresses on plant physiological processes given their significant differences in terms of the ecological environment in high-elevation areas, low anthropogenic disturbance, and obvious distribution characteristics of plants along altitudinal gradients. Therefore, plants in high-altitude areas can be used as good targets for exploring plant adaptations to abiotic stress under extreme conditions. Plants that thrive in high-altitude environments such as the Qinghai-Tibet Plateau endure extreme abiotic stresses, including low temperatures, high UV radiation, and nutrient-poor soils. This study explores their adaptation mechanisms via phenotypic variation analyses and multiomics approaches. Key findings highlight traits such as increased photosynthetic efficiency, robust antioxidant systems, and morphological modifications tailored to high-altitude conditions. These insights advance our understanding of plant evolution in harsh environments and inform strategies to increase stress resistance in crops.

Physiological and transcriptional analysis provides insights into responses of a spring wheat variety to combination of salt and heat stresses

Enhancing the frequency and intensity of extreme high temperature conditions due to global warming largely deteriorates salt-induced harm to the crop plants living in saline lands, which leads to losses of agricultural production. In northern China, spring wheat is grown in many slightly saline areas and often subjected to a combination of salt and heat stresses. In this study, a spring wheat cultivar was selected as the experimental material and subjected to salt stress (S), heat stress (H) and their combination (S + H). Physiological analysis showed that the inhibitory effect of S + H stress on wheat growth was much stronger than that of individual salt stress due to aggravating Na+ toxicity caused by heat stress. We observed that many genes involved in plant hormones showed much higher expression under S + H stress than under salt stress and heat stress, including key ABA synthesis genes (NCEDs), core ABA signalling transduction genes, key ethylene synthesis genes, and core ethylene signalling transduction genes. Particularly, many ABA-responsive genes (HSFs, HSPs, DHNs and LEAs) were upregulated under S + H stress but not under salt stress and heat stress. DHNs and LEAs were identified to play an important role in preventing cytoplasmic dehydration, protein aggregation, and slowing Na+ migration, and ethylene was identified to contribute to Na+ detoxification. We propose that in response to S + H stress, wheat plants regulate the expression of DHNs, LEAs, HSPs and HSFs via the ABA pathway to prevent cell dehydration and protein aggregation and keep ion homeostasis via the ethylene pathway.

A gene expression atlas of Nicotiana tabacum across various tissues at transcript resolution

Alternative splicing (AS) expands the transcriptome diversity by selectively splicing exons and introns from pre-mRNAs to generate different protein isoforms. This mechanism is widespread in eukaryotes and plays a crucial role in development, environmental adaptation, and stress resistance. In this study, we collected 599 tobacco RNA-seq datasets from 35 projects. 207,689 transcripts were identified in this study, of which 35,519 were annotated in the reference genome, while 172,170 transcripts were newly annotated. Additionally, tissue-specific analysis revealed 4,585 transcripts that were uniquely expressed in different tissues, highlighting the complexity and specialization of tobacco gene expression. The analysis of AS events (ASEs) across different tissues showed significant variability in the expression levels of ASE-derived transcripts, with some of these transcripts being associated with stress resistance, such as the geranyl diphosphate synthase (GGPPS). Moreover, we identified 21,763 splicing quantitative trait locus (sQTLs), which were enriched in genes involved in biological processes such as histone acetylation. Furthermore, sQTLs involved genes related to plant hormone signal transduction, terpenoid backbone biosynthesis, and other resistance pathways. These findings not only reveal the diversity of gene expression in tobacco but also provide new insights and strategies for improving tobacco quality and resistance.

Genome-Wide Identification and Expression Analysis of TCP Transcription Factors Responding to Multiple Stresses in Arachis hypogaea L

The TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING-CELL-FACTOR (TCP) gene family, a plant-specific transcription factor family, plays pivotal roles in various processes such as plant growth and development regulation, hormone crosstalk, and stress responses. However, a comprehensive genome-wide identification and characterization of the TCP gene family in peanut has yet to be fully elucidated. In this study, we conducted a genome-wide search and identified 51 TCP genes (designated as AhTCPs) in peanut, unevenly distributed across 17 chromosomes. These AhTCPs were phylogenetically classified into three subclasses: PCF, CIN, and CYC/TB1. Gene structure analysis of the AhTCPs revealed that most AhTCPs within the same subclade exhibited conserved motifs and domains, as well as similar gene structures. Cis-acting element analysis demonstrated that the AhTCP genes harbored numerous cis-acting elements associated with stress response, plant growth and development, plant hormone response, and light response. Intraspecific collinearity analysis unveiled significant collinear relationships among 32 pairs of these genes. Further collinear evolutionary analysis found that peanuts share 30 pairs, 24 pairs, 33 pairs, and 100 pairs of homologous genes with A. duranensis, A. ipaensis, Arabidopsis thaliana, and Glycine max, respectively. Moreover, we conducted a thorough analysis of the transcriptome expression profiles in peanuts across various tissues, under different hormone treatment conditions, in response to low- and high-calcium treatments, and under low-temperature and drought stress scenarios. The qRT-PCR results were in accordance with the transcriptome expression data. Collectively, these studies have established a solid theoretical foundation for further exploring the biological functions of the TCP gene family in peanuts, providing valuable insights into the regulatory mechanisms of plant growth, development, and stress responses.

Enhancing Tolerance to Combined Heat and Drought Stress in Cool-Season Grain Legumes: Mechanisms, Genetic Insights, and Future Directions

The increasing frequency of concurrent heat and drought stress poses a significant challenge to agricultural productivity, particularly for cool-season grain legumes, including broad bean (Vicia Faba L.), lupin (Lupinus spp.), lentil (Lens culinaris Medik), chickpea (Cicer arietinum L.), grasspea (Lathyrus sativus L.), pea (Pisum sativum L.), and common vetch (Vicia sativa L.). These legumes play a vital role in sustainable agricultural systems due to their nitrogen-fixing ability and high nutritional value. This review synthesizes current knowledge of the impacts and tolerance mechanisms associated with combined heat and drought stresses in these crops. We evaluate physiological and biochemical responses to combined heat and drought stress, focusing on their detrimental effects on growth, development, and yield. Key genetic and molecular mechanisms, such as the roles of osmolytes, antioxidants, and stress-responsive genes, are explored. We also discuss the intricate interplay between heat and drought stress signaling pathways, including the involvement of Ca2+ ions, reactive oxygen species, transcription factor DREB2A, and the endoplasmic reticulum in mediating stress responses. This comprehensive analysis offers new insights into developing resilient legume varieties to enhance agricultural sustainability under climate change. Future research should prioritize integrating omics technologies to unravel plant responses to combined abiotic stresses.

Soybean gene GmMLP34 regulates Arabidopsis negative response to high temperature stress

The functions of major latex proteins (MLPs) in plant defense and stress responses have been widely documented; however, their roles in HT stress response in soybeans have not been elucidated. This study investigated the role of GmMLP34, a member of the major latex protein (MLP) family, in the response of soybeans to HT stress. Transcriptome analysis of HT-resistant (JD21) and HT-sensitive (HD14) soybean leaves under HT stress (43.40 ± 1.70 °C) and field conditions revealed differential expression of GmMLP34. Further examination across different HT-resistant varieties showed that GmMLP34 was down-regulated in the leaves of 6 HT-resistant varieties (85.7 %) and up-regulated in the leaves of 6 HT-sensitive varieties (85.7 %) under the HT treatment (45 °C for 3 h). The results of this study indicate that ectopic expression of the GmMLP34 gene in Arabidopsis led to a significant decrease in the survival rate of seedling when compared to the wild type (WT) under HT stress conditions of 37/28 °C (day/night) for 5 d, Moreover, the results indicated a significant decrease in primary root length and lateral root number under 45 °C/3 h HT stress followed by 12 h room temperature recovery. Additionally, the levels of abscisic acid (ABA), and flavonoids, and the activity of the peroxidase (POD) enzyme in the antioxidant system was decreased, while the activity of the superoxide dismutase (SOD) enzyme increased in GmMLP34-overexpressing transgenic Arabidopsis thaliana. The expression levels of the HT-response genes AtCHS1 and AtCHI2-A, were significantly down-regulated, whereas that of AtGBP1 was significantly up-regulated. These results suggest that GmMLP34 negatively regulates the response of Arabidopsis thaliana to HT stress by modulating flavonoid synthesis, hormone synthesis, and the antioxidant enzyme system. These findings provide theoretical information for the genetic improvement of HT tolerance in soybean and contribute to the understanding of the molecular mechanisms underlying plant responses to abiotic stress.

Population-level exploration of alternative splicing and its unique role in controlling agronomic traits of rice

Alternative splicing (AS) plays crucial roles in regulating various biological processes in plants. However, the genetic mechanisms underlying AS and its role in controlling important agronomic traits in rice (Oryza sativa) remain poorly understood. In this study, we explored AS in rice leaves and panicles using the rice minicore collection. Our analysis revealed a high level of transcript isoform diversity, with approximately one-fifth of the potential isoforms acting as major transcripts in both tissues. Regarding the genetic mechanism of AS, we found that the splicing of 833 genes in the leaf and 1,230 genes in the panicle was affected by cis-genetic variation. Twenty-one percent of these AS events could only be explained by large structural variations. Approximately 77.5% of genes with significant splicing quantitative trait loci (sGenes) exhibited tissue-specific regulation, and AS can cause 26.9% (leaf) and 23.6% (panicle) of sGenes to have altered, lost, or gained functional domains. Additionally, through splicing-phenotype association analysis, we identified phosphate-starvation-induced RING-type E3 ligase (OsPIE1; LOC_Os01g72480), whose splicing ratio was significantly associated with plant height. In summary, this study provides an understanding of AS in rice and its contribution to the regulation of important agronomic traits.

Phenotypic, Physiological and Biochemical Delineation of Wheat Genotypes Under Different Stress Conditions

Wheat is a vital crop, providing calories, nutrients and versatility in the food industry. However, the combination of heat and drought stress, exacerbated by climate change, poses a significant threat to wheat production, leading to potential yield losses. To ensure the sustainability of wheat production it is crucial to prioritize research on developing stress-tolerant wheat genotypes. The current study focused on identifying the traits that are important for developing stress-tolerant wheat varieties under timely sown irrigated, drought stress, heat stress, and combined stress conditions. It addresses the knowledge gap regarding the combined effects of heat and drought stress on wheat physiology and yield, aiming to shed light on the intricate interactions between these stresses. The experiment was conducted at CCS HAU, Hisar, during the Rabi seasons of 2019-2020 and 2020-2021. By evaluating variability parameters, conducting correlation analysis, and path coefficient analysis among 80 diverse wheat genotypes, this research identifies genetic factors contributing to stress tolerance and helps select plants with desirable characteristics. The results showed that traits i.e., malendialdehyde, wax covering on blade, wax covering on sheath and wax covering on spike had high potential for improvement through selection among genotypes for grain yield and its component traits. The study also highlighted the importance of selecting wheat varieties with early maturity to mitigate the risk of yield loss under combined stress conditions. Moreover, the interaction between drought and heat stress can increase oxidative stress, leading to elevated malondialdehyde levels. Selecting varieties with lower malondialdehyde and optimal canopy temperature is important. Understanding the complex response of wheat to heat, drought, and their combined stress is essential for improving crop quality and production potential. Overall, this research contributes to the field of plant breeding by facilitating the development of wheat varieties with high and stable yields in challenging environments.

Functional Characterization of Accessible Chromatin in Common Wheat

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Differential induction of defense genes in hexaploid wheat roots by the plant-parasitic nematodes Pratylenchus neglectus and P. thornei

Pratylenchus neglectus and P. thornei are among the most destructive root lesion nematodes of wheat in the Pacific Northwest, United States of America and throughout the world. The aim of this study was to determine whether both nematode species were similar in their ability to induce defense genes in roots of wheat genotype Scarlet, and whether a combination of both species induced a different pattern of gene induction than each species alone. The long-term aspect of the research was to identify nematode-inducible promoters for deploying defense genes in roots in breeding programs. The root transcriptomes of genotype Scarlet were obtained after a one-week infection period with each nematode species separately, or both species combined. Root defense gene expression was induced for all three treatments relative to the no-nematode control, but P. thornei affected expression to a greater extent compared to P. neglectus. The species combination induced the highest number of defense genes. This result was not predicted from nematode enumeration studies, in which P. thornei colonization was substantially lower than that of P. neglectus, and the nematode combination did not show a significant difference. Quantitative real time polymerase chain reaction (qRT-PCR) assays for Dehydrin2, Glucan endo-1,3-beta-glucosidase, 1-cys-Peroxiredoxin, Pathogenesis-related protein 1 and Late embryogenesis-abundant proteins 76 and group 3 authenticated the induction observed in the transcriptome data. In addition, a near-isogenic line of Scarlet harboring genetic resistance to fungal soilborne pathogens, called Scarlet-Rz1, showed similar or higher levels of defense gene expression compared to fungus-susceptible Scarlet in qRT-PCR assays. Finally, transcriptome expression patterns revealed nematode-inducible promoters that are responsive to both P. neglectus and P. thornei.

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.

RNA-seq analysis reveals transcriptome reprogramming and alternative splicing during early response to salt stress in tomato root

Salt stress is one of the dominant abiotic stress conditions that cause severe damage to plant growth and, in turn, limiting crop productivity. It is therefore crucial to understand the molecular mechanism underlying plant root responses to high salinity as such knowledge will aid in efforts to develop salt-tolerant crops. Alternative splicing (AS) of precursor RNA is one of the important RNA processing steps that regulate gene expression and proteome diversity, and, consequently, many physiological and biochemical processes in plants, including responses to abiotic stresses like salt stress. In the current study, we utilized high-throughput RNA-sequencing to analyze the changes in the transcriptome and characterize AS landscape during the early response of tomato root to salt stress. Under salt stress conditions, 10,588 genes were found to be differentially expressed, including those involved in hormone signaling transduction, amino acid metabolism, and cell cycle regulation. More than 700 transcription factors (TFs), including members of the MYB, bHLH, and WRKY families, potentially regulated tomato root response to salt stress. AS events were found to be greatly enhanced under salt stress, where exon skipping was the most prevalent event. There were 3709 genes identified as differentially alternatively spliced (DAS), the most prominent of which were serine/threonine protein kinase, pentatricopeptide repeat (PPR)-containing protein, E3 ubiquitin-protein ligase. More than 100 DEGs were implicated in splicing and spliceosome assembly, which may regulate salt-responsive AS events in tomato roots. This study uncovers the stimulation of AS during tomato root response to salt stress and provides a valuable resource of salt-responsive genes for future studies to improve tomato salt tolerance.

Genome-wide association study for seedling heat tolerance under two temperature conditions in bread wheat (Triticum aestivum L.)

BACKGROUND

As the greenhouse effect intensifies, global temperatures are steadily increasing, posing a challenge to bread wheat (Triticum aestivum L.) production. It is imperative to comprehend the mechanism of high temperature tolerance in wheat and implement breeding programs to identify and develop heat-tolerant wheat germplasm and cultivars.

RESULTS

To identify quantitative trait loci (QTL) related to heat stress tolerance (HST) at seedling stage in wheat, a panel of 253 wheat accessions which were re-sequenced used to conduct genome-wide association studies (GWAS) using the factored spectrally transformed linear mixed models (FaST-LMM). For most accessions, the growth of seedlings was found to be inhibited under heat stress. Analysis of the phenotypic data revealed that under heat stress conditions, the main root length, total root length, and shoot length of seedlings decreased by 47.46%, 49.29%, and 15.19%, respectively, compared to those in normal conditions. However, 17 varieties were identified as heat stress tolerant germplasm. Through GWAS analysis, a total of 115 QTLs were detected under both heat stress and normal conditions. Furthermore, 15 stable QTL-clusters associated with heat response were identified. By combining gene expression, haplotype analysis, and gene annotation information within the physical intervals of the 15 QTL-clusters, two novel candidate genes, TraesCS4B03G0152700/TaWRKY74-B and TraesCS4B03G0501400/TaSnRK3.15-B, were responsive to temperature and identified as potential regulators of HST in wheat at the seedling stage.

CONCLUSIONS

This study conducted a detailed genetic analysis and successfully identified two genes potentially associated with HST in wheat at the seedling stage, laying a foundation to further dissect the regulatory mechanism underlying HST in wheat under high temperature conditions. Our finding could serve as genomic landmarks for wheat breeding aimed at improving adaptation to heat stress in the face of climate change.

Exploring salt tolerance mechanisms using machine learning for transcriptomic insights: case study in Spartina alterniflora

Salt stress poses a significant threat to global cereal crop production, emphasizing the need for a comprehensive understanding of salt tolerance mechanisms. Accurate functional annotations of differentially expressed genes are crucial for gaining insights into the salt tolerance mechanism. The challenge of predicting gene functions in under-studied species, especially when excluding infrequent GO terms, persists. Therefore, we proposed the use of NetGO 3.0, a machine learning-based annotation method that does not rely on homology information between species, to predict the functions of differentially expressed genes under salt stress. Spartina alterniflora, a halophyte with salt glands, exhibits remarkable salt tolerance, making it an excellent candidate for in-depth transcriptomic analysis. However, current research on the S. alterniflora transcriptome under salt stress is limited. In this study we used S. alterniflora as an example to investigate its transcriptional responses to various salt concentrations, with a focus on understanding its salt tolerance mechanisms. Transcriptomic analysis revealed substantial changes impacting key pathways, such as gene transcription, ion transport, and ROS metabolism. Notably, we identified a member of the SWEET gene family in S. alterniflora, SA_12G129900.m1, showing convergent selection with the rice ortholog SWEET15. Additionally, our genome-wide analyses explored alternative splicing responses to salt stress, providing insights into the parallel functions of alternative splicing and transcriptional regulation in enhancing salt tolerance in S. alterniflora. Surprisingly, there was minimal overlap between differentially expressed and differentially spliced genes following salt exposure. This innovative approach, combining transcriptomic analysis with machine learning-based annotation, avoids the reliance on homology information and facilitates the discovery of unknown gene functions, and is applicable across all sequenced 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.

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

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

Polyploidization: A Biological Force That Enhances Stress Resistance

Organisms with three or more complete sets of chromosomes are designated as polyploids. Polyploidy serves as a crucial pathway in biological evolution and enriches species diversity, which is demonstrated to have significant advantages in coping with both biotic stressors (such as diseases and pests) and abiotic stressors (like extreme temperatures, drought, and salinity), particularly in the context of ongoing global climate deterioration, increased agrochemical use, and industrialization. Polyploid cultivars have been developed to achieve higher yields and improved product quality. Numerous studies have shown that polyploids exhibit substantial enhancements in cell size and structure, physiological and biochemical traits, gene expression, and epigenetic modifications compared to their diploid counterparts. However, some research also suggested that increased stress tolerance might not always be associated with polyploidy. Therefore, a more comprehensive and detailed investigation is essential to complete the underlying stress tolerance mechanisms of polyploids. Thus, this review summarizes the mechanism of polyploid formation, the polyploid biochemical tolerance mechanism of abiotic and biotic stressors, and molecular regulatory networks that confer polyploidy stress tolerance, which can shed light on the theoretical foundation for future research.

Responsive Alternative Splicing Events of Opisthopappus Species against Salt Stress

Salt stress profoundly affects plant growth, prompting intricate molecular responses, such as alternative splicing (AS), for environmental adaptation. However, the response of AS events to salt stress in Opisthopappus (Opisthopappus taihangensis and Opisthopappus longilobus) remains unclear, which is a Taihang Mountain cliff-dwelling species. Using RNA-seq data, differentially expressed genes (DEGs) were identified under time and concentration gradients of salt stress. Two types of AS, skipped exon (SE) and mutually exclusive exons (MXE), were found. Differentially alternative splicing (DAS) genes in both species were significantly enriched in "protein phosphorylation", "starch and sucrose metabolism", and "plant hormone signal transduction" pathways. Meanwhile, distinct GO terms and KEGG pathways of DAS occurred between two species. Only a small subset of DAS genes overlapped with DEGs under salt stress. Although both species likely adopted protein phosphorylation to enhance salt stress tolerance, they exhibited distinct responses. The results indicated that the salt stress mechanisms of both Opisthopappus species exhibited similarities and differences in response to salt stress, which suggested that adaptive divergence might have occurred between them. This study initially provides a comprehensive description of salt responsive AS events in Opisthopappus and conveys some insights into the molecular mechanisms behind species tolerance on the Taihang Mountains.

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.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

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.

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.

Alternative splicing of TaHSFA6e modulates heat shock protein-mediated translational regulation in response to heat stress in wheat

Heat stress greatly threatens crop production. Plants have evolved multiple adaptive mechanisms, including alternative splicing, that allow them to withstand this stress. However, how alternative splicing contributes to heat stress responses in wheat (Triticum aestivum) is unclear. We reveal that the heat shock transcription factor gene TaHSFA6e is alternatively spliced in response to heat stress. TaHSFA6e generates two major functional transcripts: TaHSFA6e-II and TaHSFA6e-III. TaHSFA6e-III enhances the transcriptional activity of three downstream heat shock protein 70 (TaHSP70) genes to a greater extent than does TaHSFA6e-II. Further investigation reveals that the enhanced transcriptional activity of TaHSFA6e-III is due to a 14-amino acid peptide at its C-terminus, which arises from alternative splicing and is predicted to form an amphipathic helix. Results show that knockout of TaHSFA6e or TaHSP70s increases heat sensitivity in wheat. Moreover, TaHSP70s are localized in stress granule following exposure to heat stress and are involved in regulating stress granule disassembly and translation re-initiation upon stress relief. Polysome profiling analysis confirms that the translational efficiency of stress granule stored mRNAs is lower at the recovery stage in Tahsp70s mutants than in the wild types. Our finding provides insight into the molecular mechanisms by which alternative splicing improves the thermotolerance in wheat.

When drought meets heat - a plant omics perspective

Changes in weather patterns with emerging drought risks and rising global temperature are widespread and negatively affect crop growth and productivity. In nature, plants are simultaneously exposed to multiple biotic and abiotic stresses, but most studies focus on individual stress conditions. However, the simultaneous occurrence of different stresses impacts plant growth and development differently than a single stress. Plants sense the different stress combinations in the same or in different tissues, which could induce specific systemic signalling and acclimation responses; impacting different stress-responsive transcripts, protein abundance and modifications, and metabolites. This mini-review focuses on the combination of drought and heat, two abiotic stress conditions that often occur together. Recent omics studies indicate common or independent regulators involved in heat or drought stress responses. Here, we summarize the current research results, highlight gaps in our knowledge, and flag potential future focus areas.

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.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

The Landscapes of Gluten Regulatory Network in Elite Wheat Cultivars Contrasting in Gluten Strength

Yangmai-13 (YM13) is a wheat cultivar with weak gluten fractions. In contrast, Zhenmai-168 (ZM168) is an elite wheat cultivar known for its strong gluten fractions and has been widely used in a number of breeding programs. However, the genetic mechanisms underlying the gluten signatures of ZM168 remain largely unclear. To address this, we combined RNA-seq and PacBio full-length sequencing technology to unveil the potential mechanisms of ZM168 grain quality. A total of 44,709 transcripts were identified in Y13N (YM13 treated with nitrogen) and 51,942 transcripts in Z168N (ZM168 treated with nitrogen), including 28,016 and 28,626 novel isoforms in Y13N and Z168N, respectively. Five hundred and eighty-four differential alternative splicing (AS) events and 491 long noncoding RNAs (lncRNAs) were discovered. Incorporating the sodium-dodecyl-sulfate (SDS) sedimentation volume (SSV) trait, both weighted gene coexpression network analysis (WGCNA) and multiscale embedded gene coexpression network analysis (MEGENA) were employed for network construction and prediction of key drivers. Fifteen new candidates have emerged in association with SSV, including 4 transcription factors (TFs) and 11 transcripts that partake in the post-translational modification pathway. The transcriptome atlas provides new perspectives on wheat grain quality and would be beneficial for developing promising strategies for breeding programs.

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.

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

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

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.

Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

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.

Transcriptomic analysis reveals mechanisms for the different drought tolerance of sweet potatoes

Drought is a common environmental stress with great negative impacts on plant growth, development and geographical distribution as well as agriculture and food production. Sweet potato is characterized by starchy, fresh and pigmented tuber, and is regarded as the seventh most important food crop. However, there has been no comprehensive study of the drought tolerance mechanism of different sweet potato cultivars to date. Here, we studied the mechanism for drought response of seven sweet potato drought-tolerant cultivars using the drought coefficients, physiological indicators and transcriptome sequencing. The seven sweet potato cultivars were classified into four groups of drought tolerance performance. A large number of new genes and transcripts were identified, with an average of about 8000 new genes per sample. Alternative splicing events in sweet potato, which were dominated by first exon and last exon alternative splicing, were not conserved among different cultivars and not significantly affected by drought stress. Furthermore, different drought-tolerance mechanisms were revealed through differentially expressed gene analysis and functional annotation. Two drought-sensitive cultivars, Shangshu-9 and Xushu-22, mainly resisted drought stress by up-regulating plant signal transduction. The other drought-sensitive cultivar Jishu-26 responded to drought stress by down-regulating isoquinoline alkaloid biosynthesis and nitrogen/carbohydrate metabolism. In addition, the drought-tolerant cultivar Chaoshu-1 and drought-preferred cultivar Z15-1 only shared 9% of differentially expressed genes, as well as many opposite metabolic pathways in response to drought. They mainly regulated flavonoid and carbohydrate biosynthesis/metabolism in response to drought, while Z15-1 increased photosynthesis and carbon fixation capacity. The other drought-tolerant cultivar Xushu-18 responded to drought stress by regulating the isoquinoline alkaloid biosynthesis and nitrogen/carbohydrate metabolism. The extremely drought-tolerant cultivar Xuzi-8 was almost unaffected by drought stress and responded to drought environment only by regulating the cell wall. These findings provide important information for the selection of sweet potatoes for specific purposes.

Temporal transcriptome analysis reveals several key pathways involve in cadmium stress response in Nicotiana tabacum L

Tobacco has a strong cadmium (Cd) enrichment capacity, meaning that it can absorb large quantities from the environment, but too much Cd will cause damage to the plant. It is not yet clear how the plant can dynamically respond to Cd stress. Here, we performed a temporal transcriptome analysis of tobacco roots under Cd treatment from 0 to 48 h. The number of differentially expressed genes (DEGs) was found to change significantly at 3 h of Cd treatment, which we used to define the early and middle stages of the Cd stress response. The gene ontology (GO) term analysis indicates that genes related to photosynthesis and fatty acid synthesis were enriched during the early phases of the stress response, and in the middle phase biological process related to metal ion transport, DNA damage repair, and metabolism were enriched. It was also found that plants use precursor mRNA (pre-mRNA) processes to first resist Cd stress, and with the increasing of Cd treatment time, the overlapped genes number of DEGs and DAS increased, suggesting the transcriptional levels and post-transcriptional level might influence each other. This study allowed us to better understand how plants dynamically respond to cadmium stress at the transcriptional and post-transcriptional levels and provided a reference for the screening of Cd-tolerant genes in the future.

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

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

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

Transcriptomic and metabolic regulatory network characterization of drought responses in tobacco

Drought stress usually causes huge economic losses for tobacco industries. Drought stress exhibits multifaceted impacts on tobacco systems through inducing changes at different levels, such as physiological and chemical changes, changes of gene transcription and metabolic changes. Understanding how plants respond and adapt to drought stress helps generate engineered plants with enhanced drought resistance. In this study, we conducted multiple time point-related physiological, biochemical,transcriptomic and metabolic assays using K326 and its derived mutant 28 (M28) with contrasting drought tolerance. Through integrative analyses of transcriptome and metabolome,we observed dramatic changes of gene expression and metabolic profiles between M28 and K326 before and after drought treatment. we found that some of DEGs function as key enzymes responsible for ABA biosynthesis and metabolic pathway, thereby mitigating impairment of drought stress through ABA signaling dependent pathways. Four DEGs were involved in nitrogen metabolism, leading to synthesis of glutamate (Glu) starting from NO-3 /NO-2 that serves as an indicator for stress responses. Importantly, through regulatory network analyses, we detected several drought induced TFs that regulate expression of genes responsible for ABA biosynthesis through network, indicating direct and indirect involvement of TFs in drought responses in tobacco. Thus, our study sheds some mechanistic insights into how plant responding to drought stress through transcriptomic and metabolic changes in tobacco. It also provides some key TF or non-TF gene candidates for engineering manipulation for breeding new tobacco varieties with enhanced drought tolerance.

Alternative Splicing of TaHsfA2-7 Is Involved in the Improvement of Thermotolerance in Wheat

High temperature has severely affected plant growth and development, resulting in reduced production of crops worldwide, especially wheat. Alternative splicing (AS), a crucial post-transcriptional regulatory mechanism, is involved in the growth and development of eukaryotes and the adaptation to environmental changes. Previous transcriptome data suggested that heat shock transcription factor (Hsf) TaHsfA2-7 may form different transcripts by AS. However, it remains unclear whether this post-transcriptional regulatory mechanism of TaHsfA2-7 is related to thermotolerance in wheat (Triticum aestivum). Here, we identified a novel splice variant, TaHsfA2-7-AS, which was induced by high temperature and played a positive role in thermotolerance regulation in wheat. Moreover, TaHsfA2-7-AS is predicted to encode a small truncated TaHsfA2-7 isoform, retaining only part of the DNA-binding domain (DBD). TaHsfA2-7-AS is constitutively expressed in various tissues of wheat. Notably, the expression level of TaHsfA2-7-AS is significantly up-regulated by heat shock (HS) during flowering and grain-filling stages in wheat. Further studies showed that TaHsfA2-7-AS was localized in the nucleus but lacked transcriptional activation activity. Ectopic expression of TaHsfA2-7-AS in yeast exhibited improved thermotolerance. Compared to non-transgenic plants, overexpression of TaHsfA2-7-AS in Arabidopsis results in enhanced tolerance to heat stress. Simultaneously, we also found that TaHsfA1 is directly involved in the transcriptional regulation of TaHsfA2-7 and TaHsfA2-7-AS. In summary, our findings demonstrate the function of TaHsfA2-7-AS splicing variant in response to heat stress and establish a link between regulatory mechanisms of AS and the improvement of thermotolerance in wheat.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

Two alternative splicing variants of a wheat gene TaNAK1, TaNAK1.1 and TaNAK1.2, differentially regulate flowering time and plant architecture leading to differences in seed yield of transgenic Arabidopsis

In wheat production, appropriate flowering time and ideal plant architecture are the prerequisites for high grain yield. Alternative splicing (AS) is a vital process that regulates gene expression at the post-transcriptional level, and AS events in wheat have been found to be closely related to grain-related traits and abiotic stress tolerance. However, AS events and their biological roles in regulating flowering time and plant architecture in wheat remain unclear. In this study, we report that TaNAK1 undergoes AS, producing three splicing variants. Molecular characterization of TaNAK1 and its splicing variants demonstrated that all three protein isoforms have a conserved NB-ARC domain and a protein kinase domain, but the positions of these two domains and the length of the protein kinase domains are different among them, implying that they may have different three-dimensional structures and therefore have different functions. Further investigations showed that the two splicing variants of TaNAK1, TaNAK1.1 and TaNAK1.2, exhibited different expression patterns during wheat growth and development, while the other one, TaNAK1.3, was not detected. Subcellular localization demonstrated that TaNAK1.1 was mainly localized in the cytoplasm, while TaNAK1.2 was localized in the nucleus and cytoplasm. Both TaNAK1.1 and TaNAK1.2 exhibit protein kinase activity in vitro. Ectopic expression of TaNAK1.1 and TaNAK1.2 in Arabidopsis demonstrated that these two splicing variants play opposite roles in regulating flowering time and plant architecture, resulting in different seed yields. TaNAK1.2 positive regulates the transition from vegetative to reproductive growth, plant height, branching number, seed size, and seed yield of Arabidopsis, while TaNAK1.1 negatively regulates these traits. Our findings provide new gene resource for regulating flowering time and plant architecture in crop breeding for high grain yield.

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

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

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

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

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 novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

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.

Comprehensive profiling of alternative splicing landscape during secondary dormancy in oilseed rape (Brassica napus L.)

Alternative splicing is a general mechanism that regulates gene expression at the post-transcriptional level, which increases the transcriptomic diversity. Oilseed rape (Brassica napus L.), one of the main oil crops worldwide, is prone to secondary dormancy. However, how alternative splicing landscape of oilseed rape seed changes in response to secondary dormancy is unknown. Here, we analyzed twelve RNA-seq libraries from varieties "Huaiyou-SSD-V1" and "Huaiyou-WSD-H2" which exhibited high (> 95%) and low (< 5%) secondary dormancy potential, respectively, and demonstrated that alternative splicing changes led to a significant increase with the diversity of the transcripts in response to secondary dormancy induction via polyethylene glycol 6000 (PEG6000) treatment. Among the four basic alternative splicing types, intron retention dominates, and exon skipping shows the rarest frequency. A total of 8% of expressed genes had two or more transcripts after PEG treatment. Further analysis revealed that global isoform expression percentage variations in alternative splicing in differently expressed genes (DEGs) is more than three times as much as those in non-DEGs, suggesting alternative splicing change is associated with the transcriptional activity change in response to secondary dormancy induction. Eventually, 342 differently spliced genes (DSGs) associated with secondary dormancy were identified, five of which were validated by RT-PCR. The number of the overlapped genes between DSGs and DEGs associated with secondary dormancy was much less than that of either DSGs or DEGs, suggesting that DSGs and DEGs may independently regulates secondary dormancy. Functional annotation analysis of DSGs revealed that spliceosome components are overrepresented among the DSGs, including small nuclear ribonucleoprotein particles (snRNPs), serine/arginine-rich (SR) proteins, and other splicing factors. Thus, it is proposed that the spliceosome components could be exploited to reduce secondary dormancy potential in oilseed rape.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01314-8.