Transcriptome analysis in different rice cultivars provides novel insights into desiccation and salinity stress responses

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

OsRR26, a type-B response regulator, modulates salinity tolerance in rice via phytohormone-mediated ROS accumulation in roots and influencing reproductive development

MAIN CONCLUSION

OsRR26 is a cytokinin-responsive response regulator that promotes phytohormone-mediated ROS accumulation in rice roots, regulates seedling growth, spikelet fertility, awn development, represses NADPH oxidases, and negatively affects salinity tolerance. Plant two-component systems (TCS) play a pivotal role in phytohormone signaling, stress responses, and circadian rhythm. However, a significant knowledge gap exists regarding TCS in rice. In this study, we utilized a functional genomics approach to elucidate the role of OsRR26, a type-B response regulator in rice. Our results demonstrate that OsRR26 is responsive to cytokinin, ABA, and salinity stress, serving as the ortholog of Arabidopsis ARR11. OsRR26 primarily localizes to the nucleus and plays a crucial role in seedling growth, spikelet fertility, and the suppression of awn development. Exogenous application of cytokinin led to distinct patterns of reactive oxygen species (ROS) accumulation in the roots of both WT and transgenic plants (OsRR26OE and OsRR26KD), indicating the potential involvement of OsRR26 in cytokinin-mediated ROS signaling in roots. The application of exogenous ABA resulted in varied cellular compartmentalization of ROS between the WT and transgenic lines. Stress tolerance assays of these plants revealed that OsRR26 functions as a negative regulator of salinity stress tolerance across different developmental stages in rice. Physiological and biochemical analyses unveiled that the knockdown of OsRR26 enhances salinity tolerance, characterized by improved chlorophyll retention and the accumulation of soluble sugars, K+ content, and amino acids, particularly proline.

RNAi-Mediated Suppression of OsBBTI5 Promotes Salt Stress Tolerance in Rice

This study explores the impact of RNAi in terms of selectively inhibiting the expression of the OsBBTI5 gene, with the primary objective of uncovering its involvement in the molecular mechanisms associated with salt tolerance in rice. OsBBTI5, belonging to the Bowman-Birk inhibitor (BBI) family gene, is known for its involvement in plant stress responses. The gene was successfully cloned from rice, exhibiting transcriptional self-activation in yeast. A yeast two-hybrid assay confirmed its specific binding to OsAPX2 (an ascorbate peroxidase gene). Transgenic OsBBTI5-RNAi plants displayed insensitivity to varying concentrations of 24-epibrassinolide in the brassinosteroid sensitivity assay. However, they showed reduced root and plant height at high concentrations (10 and 100 µM) of GA3 immersion. Enzyme activity assays revealed increased peroxidase (POD) and superoxide dismutase (SOD) activities and decreased malondialdehyde (MDA) content under 40-60 mM NaCl. Transcriptomic analysis indicated a significant upregulation of photosynthesis-related genes in transgenic plants under salt stress compared to the wild type. Notably, this study provides novel insights, suggesting that the BBI gene is part of the BR signaling pathway, and that OsBBTI5 potentially enhances stress tolerance in transgenic plants through interaction with the salt stress-related gene OsAPX2.

Mechanisms of Plant Epigenetic Regulation in Response to Plant Stress: Recent Discoveries and Implications

Plant stress is a significant challenge that affects the development, growth, and productivity of plants and causes an adverse environmental condition that disrupts normal physiological processes and hampers plant survival. Epigenetic regulation is a crucial mechanism for plants to respond and adapt to stress. Several studies have investigated the role of DNA methylation (DM), non-coding RNAs, and histone modifications in plant stress responses. However, there are various limitations or challenges in translating the research findings into practical applications. Hence, this review delves into the recent recovery, implications, and applications of epigenetic regulation in response to plant stress. To better understand plant epigenetic regulation under stress, we reviewed recent studies published in the last 5-10 years that made significant contributions, and we analyzed the novel techniques and technologies that have advanced the field, such as next-generation sequencing and genome-wide profiling of epigenetic modifications. We emphasized the breakthrough findings that have uncovered specific genes or pathways and the potential implications of understanding plant epigenetic regulation in response to stress for agriculture, crop improvement, and environmental sustainability. Finally, we concluded that plant epigenetic regulation in response to stress holds immense significance in agriculture, and understanding its mechanisms in stress tolerance can revolutionize crop breeding and genetic engineering strategies, leading to the evolution of stress-tolerant crops and ensuring sustainable food production in the face of climate change and other environmental challenges. Future research in this field will continue to unveil the intricacies of epigenetic regulation and its potential applications in crop improvement.

Petunia as a model for MYB transcription factor action under salt stress

Salinity is a current and growing problem, affecting crops worldwide by reducing yields and product quality. Plants have different mechanisms to adapt to salinity; some crops are highly studied, and their salinity tolerance mechanisms are widely known. However, there are other crops with commercial importance that still need characterization of their molecular mechanisms. Usually, transcription factors are in charge of the regulation of complex processes such as the response to salinity. MYB-TFs are a family of transcription factors that regulate various processes in plant development, and both central and specialized metabolism. MYB-TFs have been studied extensively as mediators of specialized metabolism, and some are master regulators. The influence of MYB-TFs on highly orchestrated mechanisms, such as salinity tolerance, is an attractive research target. The versatility of petunia as a model species has allowed for advances to be made in multiple fields: metabolomic pathways, quality traits, stress resistance, and signal transduction. It has the potential to be the link between horticultural crops and lab models, making it useful in translating discoveries related to the MYB-TF pathways into other crops. We present a phylogenetic tree made with Petunia axillaris and Petunia inflata R2R3-MYB subfamily sequences, which could be used to find functional conservation between different species. This work could set the foundations to improve salinity resistance in other commercial crops in later studies.

Comparative transcriptomic data confirm the findings of dehydration stress-induced redox biology of indigenous aromatic rice cultivars

The present work compares the transcriptome data sets of post-imbibitional dehydration stress-raised seedlings of two contrasting indigenous aromatic rice cultivars (Oryza sativa L., Cultivars Jamainadu and Badshabhog) for unfolding genetic regulation of dehydration stress. The result of RNA-seq analysis in Illumina platform in general revealed significant cultivar-specific expression of genes under dehydration stress that substantiate the data of redox metabolic fingerprints (assessed in terms of differential efficacy of ascorbate-glutathione pathway, ROS-antioxidant interaction dynamics and sensitive biomarkers of oxidative stress). Both the cultivars showed a diverse global transcriptomic response under water-deficit condition. Transcripts selected for heatmap generation with proper annotation revealed genes that are significantly expressed and mainly involved in redox functions, signaling, membrane trafficking, replication, protein synthesis, etc. Gene ontology (GO) analysis proposed that dehydration stress in the drought-tolerant cultivar Badshabhog was attributable to the enhanced expression of genes associated with carbon dioxide-concentrating mechanism, peroxysomal biogenesis, protein modification, protein synthesis, mitochondrial electron transport chain functioning, intercellular protein transport, histone demethylation associated with developmental process, regulation of apoptosis, etc. The redox genes that got significantly over-expressed in the IARC Badshabhog vis-à-vis Jamainadu include l-ascorbate oxidase/peroxidase, monothiol glutaredoxin-S1, thioredoxin-like protein AAED1 (chloroplastic), thioredoxin-like protein CXXS1, NADH-dehydrogenase (ubiquinone)-1-beta subcomplex subunit 3-B, NADH-dehydrogenase subunit 6 and K, lipoxygenase 6 isoform-XI, etc. Overall, the results of the RNA-seq analysis led to the identification of cultivar-specific genes, with the cultivar Badshabhog exhibiting significantly greater molecular reprogramming for redox regulation and signaling necessary for combating dehydration stress.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03829-z.

Panicle transcriptome of high-yield mutant indica rice reveals physiological mechanisms and novel candidate regulatory genes for yield under reproductive stage drought stress

BACKGROUND

Reproductive stage drought stress (RDS) is a major global threat to rice production. Due to climate change, water scarcity is becoming an increasingly common phenomenon in major rice-growing areas worldwide. Understanding RDS mechanisms will allow candidate gene identification to generate novel rice genotypes tolerant to RDS.

RESULTS

To generate novel rice genotypes that can sustain yield under RDS, we performed gamma-irradiation mediated mutation breeding in the drought stress susceptible mega rice variety, MTU1010. One of the mutant MM11 (MTU1010 derived mutant11) shows consistently increased performance in yield-related traits under field conditions consecutively for four generations. In addition, compared to MTU1010, the yield of MM11 is sustained in prolonged drought imposed during the reproductive stage under field and in pot culture conditions. A comparative emerged panicle transcriptome analysis of the MTU1010 and MM11 suggested metabolic adjustment, enhanced photosynthetic ability, and hormone interplay in regulating yield under drought responses during emerged panicle development. Regulatory network analysis revealed few putative significant transcription factor (TF)-target interactions involved in integrated signalling between panicle development, yield and drought stress.

CONCLUSIONS

A gamma-irradiate rice mutant MM11 was identified by mutation breeding, and it showed higher potential to sustain yield under reproductive stage drought stress in field and pot culture conditions. Further, a comparative panicle transcriptome revealed significant biological processes and molecular regulators involved in emerged panicle development, yield and drought stress integration. The study extends our understanding of the physiological mechanisms and candidate genes involved in sustaining yield under drought stress.

Azospirillum brasilense improves rice growth under salt stress by regulating the expression of key genes involved in salt stress response, abscisic acid signaling, and nutrient transport, among others

Major food crops, such as rice and maize, display severe yield losses (30-50%) under salt stress. Furthermore, problems associated with soil salinity are anticipated to worsen due to climate change. Therefore, it is necessary to implement sustainable agricultural strategies, such as exploiting beneficial plant-microbe associations, for increased crop yields. Plants can develop associations with beneficial microbes, including arbuscular mycorrhiza and plant growth-promoting bacteria (PGPB). PGPB improve plant growth via multiple mechanisms, including protection against biotic and abiotic stresses. Azospirillum brasilense, one of the most studied PGPB, can mitigate salt stress in different crops. However, little is known about the molecular mechanisms by which A. brasilense mitigates salt stress. This study shows that total and root plant mass is improved in A. brasilense-inoculated rice plants compared to the uninoculated plants grown under high salt concentrations (100 mM and 200 mM NaCl). We observed this growth improvement at seven- and fourteen days post-treatment (dpt). Next, we used transcriptomic approaches and identified differentially expressed genes (DEGs) in rice roots when exposed to three treatments: 1) A. brasilense, 2) salt (200 mM NaCl), and 3) A. brasilense and salt (200 mM NaCl), at seven dpt. We identified 786 DEGs in the A. brasilense-treated plants, 4061 DEGs in the salt-stressed plants, and 1387 DEGs in the salt-stressed A. brasilense-treated plants. In the A. brasilense-treated plants, we identified DEGs involved in defense, hormone, and nutrient transport, among others. In the salt-stressed plants, we identified DEGs involved in abscisic acid and jasmonic acid signaling, antioxidant enzymes, sodium and potassium transport, and calcium signaling, among others. In the salt-stressed A. brasilense-treated plants, we identified some genes involved in salt stress response and tolerance (e.g., abscisic acid and jasmonic acid signaling, antioxidant enzymes, calcium signaling), and sodium and potassium transport differentially expressed, among others. We also identified some A. brasilense-specific plant DEGs, such as nitrate transporters and defense genes. Furthermore, our results suggest genes involved in auxin and ethylene signaling are likely to play an important role during these interactions. Overall, our transcriptomic data indicate that A. brasilense improves rice growth under salt stress by regulating the expression of key genes involved in defense and stress response, abscisic acid and jasmonic acid signaling, and ion and nutrient transport, among others. Our findings will provide essential insights into salt stress mitigation in rice by A. brasilense.

Metabolome plasticity in 241 Arabidopsis thaliana accessions reveals evolutionary cold adaptation processes

Acclimation and adaptation of metabolism to a changing environment are key processes for plant survival and reproductive success. In the present study, 241 natural accessions of Arabidopsis (Arabidopsis thaliana) were grown under two different temperature regimes, 16 °C and 6 °C, and growth parameters were recorded, together with metabolite profiles, to investigate the natural genome × environment effects on metabolome variation. The plasticity of metabolism, which was captured by metabolic distance measures, varied considerably between accessions. Both relative growth rates and metabolic distances were predictable by the underlying natural genetic variation of accessions. Applying machine learning methods, climatic variables of the original growth habitats were tested for their predictive power of natural metabolic variation among accessions. We found specifically habitat temperature during the first quarter of the year to be the best predictor of the plasticity of primary metabolism, indicating habitat temperature as the causal driver of evolutionary cold adaptation processes. Analyses of epigenome- and genome-wide associations revealed accession-specific differential DNA-methylation levels as potentially linked to the metabolome and identified FUMARASE2 as strongly associated with cold adaptation in Arabidopsis accessions. These findings were supported by calculations of the biochemical Jacobian matrix based on variance and covariance of metabolomics data, which revealed that growth under low temperatures most substantially affects the accession-specific plasticity of fumarate and sugar metabolism. Our findings indicate that the plasticity of metabolic regulation is predictable from the genome and epigenome and driven evolutionarily by Arabidopsis growth habitats.

Seed Halopriming: A Promising Strategy to Induce Salt Tolerance in Indonesian Pigmented Rice

Unfavorable environmental conditions and climate change impose stress on plants, causing yield losses worldwide. The Indonesian pigmented rice (Oryza sativa L.) cultivars Cempo Ireng Pendek (black rice) and Merah Kalimantan Selatan (red rice) are becoming popular functional foods due to their high anthocyanin contents and have great potential for widespread cultivation. However, their ability to grow on marginal, high-salinity lands is limited. In this study, we investigated whether seed halopriming enhances salt tolerance in the two pigmented rice cultivars. The non-pigmented cultivars IR64, a salt-stress-sensitive cultivar, and INPARI 35, a salt tolerant, were used as control. We pre-treated seeds with a halopriming solution before germination and then exposed the plants to a salt stress of 150 mM NaCl at 21 days after germination using a hydroponic system in a greenhouse. Halopriming was able to mitigate the negative effects of salinity on plant growth, including suppressing reactive oxygen species accumulation, increasing the membrane stability index (up to two-fold), and maintaining photosynthetic pigment contents. Halopriming had different effects on the accumulation of proline, in different rice varieties: the proline content increased in IR64 and Cempo Ireng Pendek but decreased in INPARI 35 and Merah Kalimantan Selatan. Halopriming also had disparate effects in the expression of stress-related genes: OsMYB91 expression was positively correlated with salt treatment, whereas OsWRKY42 and OsWRKY70 expression was negatively correlated with this treatment. These findings highlighted the potential benefits of halopriming in salt-affected agro-ecosystems.

Integrated transcriptome, proteome and metabolome analyses revealed secondary metabolites and auxiliary carbohydrate metabolism augmenting drought tolerance in rice

Drought is one of the major consequences of climate change and a serious threat to rice production. Drought stress activates interactions among genes, proteins and metabolites at the molecular level. A comparative multi-omics analysis of drought-tolerant and drought-sensitive rice cultivars can decipher the molecular mechanisms involved in drought tolerance/response. Here, we characterized the global-level transcriptome, proteome, and metabolome profiles, and performed integrated analyses thereof in a drought-sensitive (IR64) and a drought-tolerant (Nagina 22) rice cultivar under control and drought-stress conditions. The transcriptional dynamics and its integration with proteome analysis revealed the role of transporters in regulation of drought stress. The proteome response illustrated the contribution of translational machinery to drought tolerance in N22. The metabolite profiling revealed that aromatic amino acids and soluble sugars contribute majorly to drought tolerance in rice. The integrated transcriptome, proteome and metabolome analysis performed using statistical and knowledge-based methods revealed the preference for auxiliary carbohydrate metabolism through glycolysis and pentose phosphate pathway contributed to drought tolerance in N22. In addition, L-phenylalanine and the genes/proteins responsible for its biosynthesis were also found to contribute to drought tolerance in N22. In conclusion, our study provided mechanistic insights into the drought response/adaptation mechanism and is expected to facilitate engineering of drought tolerance in rice.

Comprehensive transcriptomic analysis of three varieties with different brown planthopper-resistance identifies leaf sheath lncRNAs in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) have been brought great attention for their crucial roles in diverse biological processes. However, systematic identification of lncRNAs associated with specialized rice pest, brown planthopper (BPH), defense in rice remains unexplored.

RESULTS

In this study, a genome-wide high throughput sequencing analysis was performed using leaf sheaths of susceptible rice Taichung Native 1 (TN1) and resistant rice IR36 and R476 with and without BPH feeding. A total of 2283 lncRNAs were identified, of which 649 lncRNAs were differentially expressed. During BPH infestation, 84 (120 in total), 52 (70 in total) and 63 (94 in total) of differentially expressed lncRNAs were found only in TN1, IR36 and R476, respectively. Through analyzing their cis-, trans-, and target mimic-activities, not only the lncRNAs targeting resistance genes (NBS-LRR and RLKs) and transcription factors, but also the lncRNAs acting as the targets of the well-studied stress-related miRNAs (miR2118, miR528, and miR1320) in each variety were identified. Before the BPH feeding, 238 and 312 lncRNAs were found to be differentially expressed in TN1 vs. IR36 and TN1 vs. R476, respectively. Among their putative targets, the plant-pathogen interaction pathway was significantly enriched. It is speculated that the resistant rice was in a priming state by the regulation of lncRNAs. Furthermore, the lncRNAs extensively involved in response to BPH feeding were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and the possible regulation networks of the key lncRNAs were constructed. These lncRNAs regulate different pathways that contribute to the basal defense and specific resistance of rice to the BPH.

CONCLUSION

In summary, we identified the specific lncRNAs targeting the well-studied stress-related miRNAs, resistance genes, and transcription factors in each variety during BPH infestation. Additionally, the possible regulating network of the lncRNAs extensively responding to BPH feeding revealed by WGCNA were constructed. These findings will provide further understanding of the regulatory roles of lncRNAs in BPH defense, and lay a foundation for functional research on the candidate lncRNAs.

Unearth of small open reading frames (sORFs) in drought stress transcriptome of Oryza sativa subsp. indica

Drought is a major abiotic stress that influences rice production. Although the transcriptomic data of rice against drought is widely available, the regulation of small open reading frames (sORFs) in response to drought stress in rice is yet to be investigated. Different levels of drought stress have different regulatory mechanisms in plants. In this study, drought stress was imposed on four-leaf stage rice, divided into two treatments, 40% and 30% soil moisture content (SMC). The RNAs of the samples were extracted, followed by the RNA sequencing analysis on their sORF expression changes under 40%_SMC and 30%_SMC, and lastly, the expression was validated through NanoString. A total of 122 and 143 sORFs were differentially expressed (DE) in 40%_SMC and 30%_SMC, respectively. In 40%_SMC, 69 sORFs out of 696 (9%) DEGs were found to be upregulated. On the other hand, 69 sORFs out of 449 DEGs (11%) were significantly downregulated. The trend seemed to be higher in 30%_SMC, where 112 (12%) sORFs were found to be upregulated from 928 significantly upregulated DEGs. However, only 8% (31 sORFs out of 385 DEGs) sORFs were downregulated in 30%_SMC. Among the identified sORFs, 110 sORFs with high similarity to rice proteome in the PsORF database were detected in 40%_SMC, while 126 were detected in 30%_SMC. The Gene Ontology (GO) enrichment analysis of DE sORFs revealed their involvement in defense-related biological processes, such as defense response, response to biotic stimulus, and cellular homeostasis, whereas enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that DE sORFs were associated with tryptophan and phenylalanine metabolisms. Several DE sORFs were identified, including the top five sORFs (OsisORF_3394, OsisORF_0050, OsisORF_3007, OsisORF_6407, and OsisORF_7805), which have yet to be characterised. Since these sORFs were responsive to drought stress, they might hold significant potential as targets for future climate-resilient rice development.

Transcriptomic Analysis Reveals Panicle Heterosis in an Elite Hybrid Rice ZZY10 and Its Parental Lines

Heterosis is the phenomenon in which some hybrid traits are superior to those of their parents. Most studies have analyzed the heterosis of agronomic traits of crops; however, heterosis of the panicles can improve yield and is important for crop breeding. Therefore, a systematic study of panicle heterosis is needed, especially during the reproductive stage. RNA sequencing (RNA Seq) and transcriptome analysis are suitable for further study of heterosis. Using the Illumina Nova Seq platform, the transcriptome of ZhongZheYou 10 (ZZY10), an elite rice hybrid, the maintainer line ZhongZhe B (ZZB), and the restorer line Z7-10 were analyzed at the heading date in Hangzhou, 2022. 581 million high-quality short reads were obtained by sequencing and were aligned against the Nipponbare reference genome. A total of 9000 differential expression genes were found between the hybrids and their parents (DG_HP). Of the DG_HP, 60.71% were up-regulated and 39.29% were down-regulated in the hybrid. Comparative transcriptome analysis revealed that 5235 and 3765 DG_HP were between ZZY10 and ZhongZhe B and between ZZY10 and Z7-10, respectively. This result is consistent with the transcriptome profile of ZZY10 and was similar to Z7-10. The expression patterns of DG_HP mainly exhibited over-dominance, under-dominance, and additivity. Among the DG_HP-involved GO terms, pathways such as photosynthesis, DNA integration, cell wall modification, thylakoid, and photosystem were significant. 21 DG_HP, which were involved in photosynthesis, and 17 random DG_HP were selected for qRT-PCR validation. The up-regulated PsbQ and down-regulated subunits of PSI and PSII and photosynthetic electron transport in the photosynthesis pathway were observed in our study. Extensive transcriptome data were obtained by RNA-Seq, providing a comprehensive overview of panicle transcriptomes at the heading stage in a heterotic hybrid.

Genome-wide discovery of genetic variations between rice cultivars with contrasting drought stress response and their potential functional relevance

Drought stress is a serious threat to rice productivity. Investigating genetic variations between drought-tolerant (DT) and drought-sensitive (DS) rice cultivars may decipher the candidate genes/regulatory regions involved in drought stress tolerance/response. In this study, whole-genome resequencing data of four DS and five DT rice cultivars were analyzed. We identified a total of approximately 4.8 million single nucleotide polymorphisms (SNPs) and 0.54 million insertions/deletions (InDels). The genetic variations (162,638 SNPs and 17,217 InDels) differentiating DS and DT rice cultivars were found to be unevenly distributed throughout the rice genome; however, they were more frequent near the transcription start and stop sites than in the genic regions. The cis-regulatory motifs representing the binding sites of stress-related transcription factors (MYB, HB, bZIP, ERF, ARR, and AREB) harboring the SNPs/InDels in the promoter regions of a few differentially expressed genes (DEGs) were identified. Importantly, many of these DEGs were located within the drought-associated quantitative trait loci. Overall, this study provides a valuable large-scale genotyping resource and facilitates the discovery of candidate genes associated with drought stress tolerance in rice.

New Insights into rice pyrimidine catabolic enzymes

INTRODUCTION

Rice is a primary global food source, and its production is affected by abiotic stress, caused by climate change and other factors. Recently, the pyrimidine reductive catabolic pathway, catalyzed by dihydropyrimidine dehydrogenase (DHPD), dihydropyrimidinase (DHP) and β-ureidopropionase (β-UP), has emerged as a potential participant in the abiotic stress response of rice.

METHODS

The rice enzymes were produced as recombinant proteins, and two were kinetically characterized. Rice dihydroorotate dehydrogenase (DHODH), an enzyme of pyrimidine biosynthesis often confused with DHPD, was also characterized. Salt-sensitive and salt-resistant rice seedlings were subjected to salt stress (24 h) and metabolites in leaves were determined by mass spectrometry.

RESULTS

The OsDHPD sequence was homologous to the C-terminal half of mammalian DHPD, conserving FMN and uracil binding sites, but lacked sites for Fe/S clusters, FAD, and NADPH. OsDHPD, truncated to eliminate the chloroplast targeting peptide, was soluble, but inactive. Database searches for polypeptides homologous to the N-terminal half of mammalian DHPD, that could act as co-reductants, were unsuccessful. OsDHODH exhibited kinetic parameters similar to those of other plant DHODHs. OsDHP, truncated to remove a signal sequence, exhibited a kcat/Km = 3.6 x 103 s-1M-1. Osb-UP exhibited a kcat/Km = 1.8 x 104 s-1M-1. Short-term salt exposure caused insignificant differences in the levels of the ureide intermediates dihydrouracil and ureidopropionate in leaves of salt-sensitive and salt-resistant plants. Allantoin, a ureide metabolite of purine catabolism, was found to be significantly higher in the resistant cultivar compared to one of the sensitive cultivars.

DISCUSSION

OsDHP, the first plant enzyme to be characterized, showed low kinetic efficiency, but its activity may have been affected by truncation. Osb-UP exhibited kinetic parameters in the range of enzymes of secondary metabolism. Levels of two pathway metabolites were similar in sensitive and resistant cultivars and appeared to be unaffected by short-term salt exposure."

Genome-wide identification of the class III peroxidase gene family of sugarcane and its expression profiles under stresses

INTRODUCTION

Plant-specific Class III peroxidases (PRXs) play a crucial role in lignification, cell elongation, seed germination, and biotic and abiotic stresses.

METHODS

The class III peroxidase gene family in sugarcane were identified by bioinformatics methods and realtime fluorescence quantitative PCR.

RESULTS

Eighty-two PRX proteins were characterized with a conserved PRX domain as members of the class III PRX gene family in R570 STP. The ShPRX family genes were divided into six groups by the phylogenetic analysis of sugarcane, Saccharum spontaneum, sorghum, rice, and Arabidopsis thaliana. The analysis of promoter cis-acting elements revealed that most ShPRX family genes contained cis-acting regulatory elements involved in ABA, MeJA, light responsiveness, anaerobic induction, and drought inducibility. An evolutionary analysis indicated that ShPRXs was formed after Poaceae and Bromeliaceae diverged, and tandem duplication events played a critical role in the expansion of ShPRX genes of sugarcane. Purifying selection maintained the function of ShPRX proteins. SsPRX genes were differentially expressed in stems and leaves at different growth stages in S. spontaneum. However, ShPRX genes were differentially expressed in the SCMV-inoculated sugarcane plants. A qRT-PCR analysis showed that SCMV, Cd, and salt could specifically induce the expression of PRX genes of sugarcane.

DISCUSSION

These results help elucidate the structure, evolution, and functions of the class III PRX gene family in sugarcane and provide ideas for the phytoremediation of Cd-contaminated soil and breeding new sugarcane varieties resistant to sugarcane mosaic disease, salt, and Cd stresses.

Transcriptome and Physio-Biochemical Profiling Reveals Differential Responses of Rice Cultivars at Reproductive-Stage Drought Stress

Drought stress severely affects the growth and development of rice, especially at the reproductive stage, which results in disturbed metabolic processes, reduced seed-set/grain filling, deteriorated grain quality, declined productivity, and lower yield. Despite the recent advances in understanding the responses of rice to drought stress, there is a need to comprehensively integrate the morpho-physio-biochemical studies with the molecular responses/differential expression of genes and decipher the underlying pathways that regulate the adaptability of rice at various drought-sensitive growth stages. Our comparative analysis of immature panicle from a drought-tolerant (Nagina 22) and a drought-sensitive (IR 64) rice cultivar grown under control (well-watered) and water-deficit/drought stress (treatment, imposed at the reproductive stage) conditions unraveled some novel stress-responsive genes/pathways responsible for reproductive-stage drought stress tolerance. The results revealed a more important role of upregulated (6706) genes in the panicle of N 22 at reproductive-stage drought stress compared to that (5590) in IR 64. Functional enrichment and MapMan analyses revealed that majority of the DEGs were associated with the phytohormone, redox signalling/homeostasis, secondary metabolite, and transcription factor-mediated mitigation of the adverse effects of drought stress in N 22. The upregulated expression of the genes associated with starch/sucrose metabolism, secondary metabolites synthesis, transcription factors, glutathione, linoleic acid, and phenylalanine metabolism in N 22 was significantly more than that in the panicle of IR 64. Compared to IR 64, 2743 genes were upregulated in N 22 under control conditions, which further increased (4666) under drought stress in panicle of the tolerant cultivar. Interestingly, we observed 6706 genes to be upregulated in the panicle of N 22 over IR 64 under drought and 5814 genes get downregulated in the panicle of N 22 over IR 64 under the stress. In addition, RT-qPCR analysis confirmed differential expression patterns of the DEGs. These genes/pathways associated with the reproductive-stage drought tolerance might provide an important source of molecular markers for genetic manipulation of rice for enhanced drought tolerance.

Long non-coding RNAs as the regulatory hubs in rice response to salt stress

Salinity seriously constrains growth and fertility of rice worldwide. Long non-coding RNAs (lncRNAs) play crucial roles in plant abiotic stress response. However, salt responsive lncRNAs are poorly understood in rice. Herein, salt responsive lncRNAs (DE-lncRNAs) were identified in FL478 (salt tolerant) compared to its susceptible parent (IR29) using RNA-seq in root tissues at seedling stage. In FL478 and IR29, 8724 and 9235 transcripts with length of > 200 bp were nominated as potential lncRNAs, respectively. Rigorous filtering left four (in FL478) and nine (in IR29) DE-lncRNAs with only 2 DE-lncRNAs in common. ATAC-seq data showed that the genomic regions of all four lncRNAs in FL478 and 6/9 in IR29 are significantly accessible for transcription. Weighted correlation network analysis (WGCNA) revealed that lncRNA.2-FL was highly correlated with 173 mRNAs as trans-targets and a gene encoding pentatricopeptide repeat (PPR) protein was predicted as cis-target of lncRNA.2-FL. In silico mutagenesis analysis proposed the same transcription factor binding sites (TFBSs) in vicinity of the trans- and cis-regulatory target genes of lncRNA.2-FL, which significantly affect their transcription start site (TSS). This study provides new insights into involvement of the DE-lncRNAs in rice response to salt stress. Among them, lncRNA.2-FL may play a significant regulatory role in the salt stress tolerance of FL478.

Comparative time-course transcriptome analysis of two contrasting alfalfa (Medicago sativa L.) genotypes reveals tolerance mechanisms to salt stress

Salt stress is a major abiotic stress affecting plant growth and crop yield. For the successful cultivation of alfalfa (Medicago sativa L.), a key legume forage, in saline-affected areas, it's essential to explore genetic modifications to improve salt-tolerance.Transcriptome assay of two comparative alfalfa genotypes, Adina and Zhaodong, following a 4 h and 8 h's 300 mM NaCl treatment was conducted in this study in order to investigate the molecular mechanism in alfalfa under salt stress conditions. Results showed that we obtained 875,023,571 transcripts and 662,765,594 unigenes were abtained from the sequenced libraries, and 520,091 assembled unigenes were annotated in at least one database. Among them, we identified 1,636 differentially expression genes (DEGs) in Adina, of which 1,426 were up-regulated and 210 down-regulated, and 1,295 DEGs in Zhaodong, of which 565 were up-regulated and 730 down-regulated. GO annotations and KEGG pathway enrichments of the DEGs based on RNA-seq data indicated that DEGs were involved in (1) ion and membrane homeostasis, including ABC transporter, CLC, NCX, and NHX; (2) Ca²⁺ sensing and transduction, including BK channel, EF-hand domain, and calmodulin binding protein; (3) phytohormone signaling and regulation, including TPR, FBP, LRR, and PP2C; (4) transcription factors, including zinc finger proteins, YABBY, and SBP-box; (5) antioxidation process, including GST, PYROX, and ALDH; (6) post-translational modification, including UCH, ubiquitin family, GT, MT and SOT. The functional roles of DEGs could explain the variations in salt tolerance performance observed between the two alfalfa genotypes Adina and Zhaodong. Our study widens the understanding of the sophisticated molecular response and tolerance mechanism to salt stress, providing novel insights on candidate genes and pathways for genetic modification involved in salt stress adaptation in alfalfa.

ERF subfamily transcription factors and their function in plant responses to abiotic stresses

Ethylene Responsive Factor (ERF) subfamily comprise the largest number of proteins in the plant AP2/ERF superfamily, and have been most extensively studied on the biological functions. Members of this subfamily have been proven to regulate plant resistances to various abiotic stresses, such as drought, salinity, chilling and some other adversities. Under these stresses, ERFs are usually activated by mitogen-activated protein kinase induced phosphorylation or escape from ubiquitin-ligase enzymes, and then form complex with nucleic proteins before binding to cis-element in promoter regions of stress responsive genes. In this review, we will discuss the phylogenetic relationships among the ERF subfamily proteins, summarize molecular mechanism how the transcriptional activity of ERFs been regulated and how ERFs of different subgroup regulate the transcription of stress responsive genes, such as high-affinity K⁺ transporter gene PalHKT1;2, reactive oxygen species related genes LcLTP, LcPrx, and LcRP, flavonoids synthesis related genes FtF3H and LhMYBSPLATTER, etc. Though increasing researches demonstrate that ERFs are involved in various abiotic stresses, very few interact proteins and target genes of them have been comprehensively annotated. Hence, future research prospects are described on the mechanisms of how stress signals been transited to ERFs and how ERFs regulate the transcriptional expression of stress responsive genes.

Variety-Specific Transcriptional and Alternative Splicing Regulations Modulate Salt Tolerance in Rice from Early Stage of Stress

Salt stress poses physiological drought, ionic toxicity and oxidative stress to plants, which causes premature senescence and death of the leaves if the stress sustained. Salt tolerance varied between different rice varieties, but how different rice varieties respond at the early stage of salt stress has been seldom studied comprehensively. By employing third generation sequencing technology, we compared gene expressional changes in leaves of three rice varieties that varied in their level of tolerance after salt stress treatment for 6 h. Commonly up-regulated genes in all rice varieties were related to water shortage response and carbon and amino acids metabolism at the early stage of salt stress, while reactive oxygen species cleavage genes were induced more in salt-tolerant rice. Unexpectedly, genes involved in chloroplast development and photosynthesis were more significantly down-regulated in the two salt tolerant rice varieties 'C34' and 'Nona Bokra'. At the same time, genes coding ribosomal protein were suppressed to a more severe extent in the salt-sensitive rice variety 'IR29'. Interestingly, not only variety-specific gene transcriptional regulation, but also variety-specific mRNA alternative splicing, on both coding and long-noncoding genes, were found at the early stage of salt stress. In summary, differential regulation in gene expression at both transcriptional and post-transcriptional levels, determine and fine-tune the observed response in level of damage in leaves of specific rice genotypes at early stage of salt stress.

Molecular basis of genetic plasticity to varying environmental conditions on growing rice by dry/direct-sowing and exposure to drought stress: Insights for DSR varietal development

Rice requires plenty of water for its cultivation by transplanting. This poses several challenges to its cultivation due to erratic rainfall resulting in drought, flood, and other abiotic stresses of varying intensity. Dry/direct-sown rice (DSR) has emerged as a water-saving/climate-smart alternative to transplanted rice (TPR). The performance of a rice cultivar on growing by different methods of planting under varying environmental conditions varies considerably. However, the molecular basis of the observed phenotypic plasticity of rice to varying environmental conditions is still elusive. Resilience to various environmental fluctuations is important to ensure sustainable rice production in the present era of global climate change. Our observations on exclusively up-regulated genes in leaf of Nagina 22 (N 22) grown by dry/direct-sowing and subjected to drought stress at panicle initiation stage (compared to that in leaf of IR 64), and another set of genes exclusively down-regulated in leaf of N 22 (compared to that in leaf of IR 64) indicate important roles of leaf in stress resilience. A large number of genes down-regulated exclusively in root of N 22 on dry/direct-sowing subjected to drought stress indicates a major contribution of roots in stress tolerance. The genes for redox-homeostasis, transcription factors, stress signaling, carbohydrate metabolism, and epigenetic modifications play important roles in making N 22 better adapted to DSR conditions. More importantly, the involvement of genes in rendering genetic plasticity to N 22 under changing environmental conditions was confirmed by reversal of the method of planting. To the best of our knowledge, this is the first report on decoding the molecular basis of genetic plasticity of rice grown by two different methods of planting subjected to drought stress at the reproductive stage of plant growth. This might help in DSR varietal development program to enhance water-productivity, conserve natural resources, and minimize the emission of greenhouse gases, thus achieving the objectives of negative-emission agriculture.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

Integrated analysis of transcriptome and metabolome reveals molecular mechanisms of salt tolerance in seedlings of upland rice landrace 17SM-19

Salt stress is a major abiotic stress that threatens global rice production. It is particularly important to improve salt tolerance in upland rice because of its growth environment. Upland rice landrace 17SM-19 with high salt tolerance was obtained from a previous study. In this study, an integrated analysis of transcriptome and metabolome was performed to determine the responses of the rice seedling to salt stress. When treated with 100 mm NaCl, the rice seedling growth was significantly inhibited at 5 d, with inhibition first observed in shoot dry weight (SDW). Changes in potassium (K⁺) content were associated with changes in SDW. In omics analyses, 1,900 differentially expressed genes (DEGs) and 659 differentially abundant metabolites (DAMs) were identified at 3 d after salt stress (DAS), and 1,738 DEGs and 657 DAMs were identified at 5 DAS. Correlation analyses between DEGs and DAMs were also conducted. The results collectively indicate that salt tolerance of upland rice landrace 17SM-19 seedlings involves many molecular mechanisms, such as those involved with osmotic regulation, ion balance, and scavenging of reactive oxygen species.

Artemisia argyi allelopathy: a generalist compromises hormone balance, element absorption, and photosynthesis of receptor plants

BACKGROUND

Allelopathy is expressed through the release of plant chemicals and is considered a natural alternative for sustainable weed management. Artemisia argyi (A. argyi) is widely distributed throughout Asia, and often dominates fields due to its strong allelopathy. However, the mechanism of A. argyi allelopathy is largely unknown and need to be elucidated at the physiological and molecular levels.

RESULTS

In this study, we used electron microscopy, ionomics analysis, phytohormone profiling, and transcriptome analysis to investigate the physiological and molecular mechanisms of A. argyi allelopathy using the model plant rice (Oryza sativa) as receptor plants. A. argyi water extract (AAWE)-treated rice plants grow poorly and display root morphological anomalies and leaf yellowing. We found that AAWE significantly inhibits rice growth by destroying the root and leaf system in multiple ways, including the integrity of ultrastructure, reactive oxygen species (ROS) homeostasis, and the accumulation of soluble sugar and chlorophyll synthesis. Further detection of the hormone contents suggests that AAWE leads to indole-3-acetic acid (IAA) accumulation in roots. Moreover, ionomics analysis shows that AAWE inhibits the absorption and transportation of photosynthesis-essential mineral elements, especially Mg, Fe, and Mn. In addition, the results of transcriptome analysis revealed that AAWE affects a series of crucial primary metabolic processes comprising photosynthesis in rice plants.

CONCLUSIONS

This study indicates that A. argyi realizes its strongly allelopathy through comprehensive effects on recipient plants including large-scale IAA synthesis and accumulation, ROS explosion, damaging the membrane system and organelles, and obstructing ion absorption and transport, photosynthesis and other pivotal primary metabolic processes of plants. Therefore, AAWE could potentially be developed as an environmentally friendly botanical herbicide due to its strong allelopathic effects.

Transcriptome Analysis Reveals Molecular Mechanisms under Salt Stress in Leaves of Foxtail Millet (Setaria italica L.)

Foxtail millet (Setaria italica L.) is an important cereal for managing future water scarcity and ensuring food security, due to its strong drought and salt stress resistance owing to its developed root system. However, the molecular responses of foxtail millet leaves to salt stress are largely unknown. In this study, seeds of 104 foxtail millet accessions were subjected to 0.17 mol·L^-1 NaCl stress during germination, and various germination-related parameters were analyzed to derive 5 salt-sensitive accessions and 13 salt-tolerant accessions. Hong Gu 2000 and Pu Huang Yu were the most salt-tolerant and salt-sensitive accessions, respectively. To determine the mechanism of the salt stress response, transcriptomic differences between the control and salt-treated groups were investigated. We obtained 2019 and 736 differentially expressed genes under salt stress in the salt-sensitive and salt-tolerant accessions, respectively. The transcription factor families bHLH, WRKY, AP2/ERF, and MYB-MYC were found to play critical roles in foxtail millet's response to salt stress. Additionally, the down-regulation of ribosomal protein-related genes causes stunted growth in the salt-sensitive accessions. The salt-tolerant accession alleviates salt stress by increasing energy production. Our findings provide novel insights into the molecular mechanism of foxtail millet's response to salt stress.

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.

A Review of Integrative Omic Approaches for Understanding Rice Salt Response Mechanisms

Soil salinity is one of the most serious environmental challenges, posing a growing threat to agriculture across the world. Soil salinity has a significant impact on rice growth, development, and production. Hence, improving rice varieties’ resistance to salt stress is a viable solution for meeting global food demand. Adaptation to salt stress is a multifaceted process that involves interacting physiological traits, biochemical or metabolic pathways, and molecular mechanisms. The integration of multi-omics approaches contributes to a better understanding of molecular mechanisms as well as the improvement of salt-resistant and tolerant rice varieties. Firstly, we present a thorough review of current knowledge about salt stress effects on rice and mechanisms behind rice salt tolerance and salt stress signalling. This review focuses on the use of multi-omics approaches to improve next-generation rice breeding for salinity resistance and tolerance, including genomics, transcriptomics, proteomics, metabolomics and phenomics. Integrating multi-omics data effectively is critical to gaining a more comprehensive and in-depth understanding of the molecular pathways, enzyme activity and interacting networks of genes controlling salinity tolerance in rice. The key data mining strategies within the artificial intelligence to analyse big and complex data sets that will allow more accurate prediction of outcomes and modernise traditional breeding programmes and also expedite precision rice breeding such as genetic engineering and genome editing.

Transcriptome Time-Course Analysis in the Whole Period of Cotton Fiber Development

Gossypium hirsutum and Gossypium barbadense are the widely cultivated tetraploid cottons around the world, which evolved great differences in the fiber yield and quality due to the independent domestication process. To reveal the genetic basis of the difference, we integrated 90 samples from ten time points during the fiber developmental period for investigating the dynamics of gene expression changes associated with fiber in G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124 and acc. 3-79. Globally, 44,484 genes expressed in all three cultivars account for 61.14% of the total genes. About 61.39% (N = 3,412) of the cotton transcription factors were involved in fiber development, which consisted of 58 cotton TF families. The differential analysis of intra- and interspecies showed that 3 DPA had more expression changes. To discover the genes with temporally changed expression profiles during the whole fiber development, 1,850 genes predominantly expressed in G. hirsutum and 1,050 in G. barbadense were identified, respectively. Based on the weighted gene co-expression network and time-course analysis, several candidate genes, mainly involved in the secondary cell wall synthesis and phytohormones, were identified in this study, underlying possibly the transcriptional regulation and molecular mechanisms of the fiber quality differences between G. barbadense and G. hirsutum. The quantitative real-time PCR validation of the candidate genes was consistent with the RNA-seq data. Our study provides a strong rationale for the analysis of gene function and breeding of high-quality cotton.

Identifying Novel Genes and Proteins Involved in Salt Stress Perception and Signaling of Rice Seedlings

Rice is one of the most important crops worldwide. Crop production is constrained markedly, however, by abiotic stresses such as salinity. To elucidate early stress response signaling networks involved in rice, we report in this study an original quantitative proteomic analysis of the rice seedlings subjected to short-term salt stress. We detected 570 differentially regulated proteins (DRPs) in the root sample. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis demonstrated that DRPs of the root were mainly involved in membrane trafficking, kinase activity, and ion toxicity responses. Interactome analysis revealed the central role of root proteins involved in membrane trafficking in the early response to salinity, such as cell surface receptor-like kinases (RLKs), phosphatidylinositols (PIs), calcium-dependent protein kinases 1 and 5, calcineurin B-like protein-interacting proteins, protein phosphatase 2C (PP2C) inhibitors, and abscisic acid receptors (PYL5/10), indicating activation of S-type anion channel. Furthermore, the proteogenomic analysis revealed 128 unique genome search-specific peptides with high-quality mass spectromety (MS/MS) spectra. We identified 38 novel protein-coding genes, refined the annotation of 17 existing gene models, and suggested several novel stress-responsive proteins, such as RLK5, peroxidase 27, and growth-regulating factor 2. Novel peptides had an ortholog match in the curated protein sequence set of other plant species. In conclusion, this study identifies novel stress-responsive proteins and genes of rice, thus warrant future consideration as candidates for molecular breeding of stress-tolerant crop varieties.

Seedling-stage salinity tolerance in rice: Decoding the role of transcription factors

Rice is an important staple food crop that feeds over half of the human population, particularly in developing countries. Increasing salinity is a major challenge for continuing rice production. Though rice is affected by salinity at all the developmental stages, it is most sensitive at the early seedling stage. The yield thus depends on how many seedlings can withstand saline water at the stage of transplantation, especially in coastal farms. The rapid development of "omics" approaches has assisted researchers in identifying biological molecules that are responsive to salt stress. Several salinity-responsive quantitative trait loci (QTL) contributing to salinity tolerance have been identified and validated, making it essential to narrow down the search for the key genes within QTLs. Owing to the impressive progress of molecular tools, it is now clear that the response of plants toward salinity is highly complex, involving multiple genes, with a specific role assigned to the repertoire of transcription factors (TF). Targeting the TFs for improving salinity tolerance can have an inbuilt advantage of influencing multiple downstream genes, which in turn can contribute toward tolerance to multiple stresses. This is the first comparative study for TF-driven salinity tolerance in contrasting rice cultivars at the seedling stage that shows how tolerant genotypes behave differently than sensitive ones in terms of stress tolerance. Understanding the complexity of salt-responsive TF networks at the seedling stage will be helpful to alleviate crop resilience and prevent crop damage at an early growth stage in rice.

RNA-seq analysis reveals the genes/pathways responsible for genetic plasticity of rice to varying environmental conditions on direct-sowing and transplanting

Rice cultivation by transplanting requires plenty of water. It might become a challenging task in future to grow rice by transplanting due to the climatic change, water and labor scarcities. Direct-sown rice (DSR) is emerging as a resource-conserving and climate-smart alternative to transplanted rice (TPR). However, no specific variety has been bred for dry/direct-sown conditions. The present study was undertaken to decipher the molecular basis of genetic plasticity of rice under different planting methods. Comparative RNA-seq analysis revealed a number (6133) of genes exclusively up-regulated in Nagina-22 (N-22) leaf under DSR conditions, compared to that (3538) in IR64 leaf. Several genes up-regulated in N-22 were down-regulated in IR64. Genes for growth-regulation and nutrient-reservoir activities, transcription factors, translational machinery, carbohydrate metabolism, cell cycle/division, and chromatin organization/epigenetic modifications were considerably up-regulated in the leaf of N-22 under DSR conditions. Complementary effects of these factors in rendering genetic plasticity were confirmed by the agronomic/physiological performance of rice cultivar. Thus, growth-regulation/nutrient-reservoir activities, transcription factors, and translational machinery are important molecular factors responsible for the observed genetic plasticity/adaptability of Nagina-22 to different planting methods. This might help to develop molecular markers for DSR breeding, replacing TPR with DSR for better water-productivity, and minimizing greenhouse-gas emission necessary for negative emission agriculture.

Molecular Insights into Salinity Responsiveness in Contrasting Genotypes of Rice at the Seedling Stage

Salinity is one of the most common unfavorable environmental conditions that limits plant growth and development, ultimately reducing crop productivity. To investigate the underlying molecular mechanism involved in the salinity response in rice, we initially screened 238 rice cultivars after salt treatment at the seedling stage and identified two highly salt-tolerant cultivars determined by the relative damage rate parameter. The majority of cultivars (94.1%) were ranked as salt-sensitive and highly salt-sensitive. Transcriptome profiling was completed in highly salt-tolerant, moderately salt-tolerant, and salt-sensitive under water and salinity treatments at the seedling stage. Principal component analysis displayed a clear distinction among the three cultivars under control and salinity stress conditions. Several starch and sucrose metabolism-related genes were induced after salt treatment in all genotypes at the seedling stage. The results from the present study enable the identification of the ascorbate glutathione pathway, potentially participating in the process of plant response to salinity in the early growth stage. Our findings also highlight the significance of high-affinity K⁺ uptake transporters (HAKs) and high-affinity K⁺ transporters (HKTs) during salt stress responses in rice seedlings. Collectively, the cultivar-specific stress-responsive genes and pathways identified in the present study act as a useful resource for researchers interested in plant responses to salinity at the seedling stage.

Transcriptome analysis of Auricularia fibrillifera fruit-body responses to drought stress and rehydration

BACKGROUND

Drought stress severely restricts edible fungus production. The genus Auricularia has a rare drought tolerance, a rehydration capability, and is nutrient rich.

RESULTS

The key genes and metabolic pathways involved in drought-stress and rehydration were investigated using a transcriptome analysis to clarify the relevant molecular mechanisms. In total, 173.93 Mb clean reads, 26.09 Gb of data bulk, and 52,954 unigenes were obtained. Under drought-stress and rehydration conditions, 14,235 and 8539 differentially expressed genes, respectively, were detected. 'Tyrosine metabolic', 'caffeine metabolism', 'ribosome', 'phagosome', and 'proline and arginine metabolism', as well as 'peroxisome' and 'mitogen-activated protein kinase signaling' pathways, had major roles in A. fibrillifera responses to drought stress. 'Tyrosine' and 'caffeine metabolism' might reveal unknown mechanisms for the antioxidation of A. fibrillifera under drought-stress conditions. During the rehydration process, 'diterpenoid biosynthesis', 'butanoate metabolism', 'C5-branched dibasic acid', and 'aflatoxin biosynthesis' pathways were significantly enriched. Gibberellins and γ-aminobutyric acid were important in the recovery of A. fibrillifera growth after rehydration. Many genes related to antibiotics, vitamins, and other health-related ingredients were found in A. fibrillifera.

CONCLUSION

These findings suggested that the candidate genes and metabolites involved in crucial biological pathways might regulate the drought tolerance or rehydration of Auricularia, shedding light on the corresponding mechanisms and providing new potential targets for the breeding and cultivation of drought-tolerant fungi.

Variety-specific transcript accumulation during reproductive stage in drought-stressed rice

The divergence of natural stress tolerance mechanisms between species is an intriguing phenomenon. To study it in rice, a comparative transcriptome analysis was carried out in 'heading' stage tissue (flag leaf, panicles and roots) of Nagina 22 (N22; drought-tolerant) and IR64 (drought-sensitive) plants subjected to field drought. Interestingly, N22 showed almost double the number of differentially expressed genes (DEGs) than IR64. Many DEGs colocalized within drought-related QTLs responsible for grain yield and drought tolerance and also associated with drought tolerance and critical drought-related plant traits such as leaf rolling, trehalose content, sucrose and cellulose content. Besides, co-expression analysis of the DEGs revealed several 'hub' genes known to actively regulate drought stress response. Strikingly, 1366 DEGs, including 21 'hub' genes, showed a distinct opposite regulation in the two rice varieties under similar drought conditions. Annotation of these variety-specific DEGs (VS-DEGs) revealed that they are distributed in various biological pathways. Furthermore, 103 VS-DEGs were found to physically interact with over 1300 genes, including 32 that physically interact with other VS-DEGs as well. The promoter region of these genes has sequence variations among the two rice varieties, which might be in part responsible for their unique expression pattern.

Germin like protein genes exhibit modular expression during salt and drought stress in elite rice cultivars

BACKGROUND

Germin-like proteins (GLPs) are ubiquitous plant proteins, which play significant role in plant responses against various abiotic stresses. However, the potential functions of GLPs in rice (Oryza sativa) against salt and drought stress are still unclear.

METHODS AND RESULTS

In this study, transcriptional variation of eight OsGLP genes (OsGLP3-6, OsGLP4-1, OsGLP8-4, OsGLP8-7, OsGLP8-10, OsGLP8-11 and OsGLP8-12) was analyzed in leaves and roots of two economically important Indica rice cultivars, KS282 and Super Basmati, under salt and drought stress at early seedling stage. The relative expression analysis from qRT-PCR indicated the highest increase in expression of OsGLP3-6 in leaves and roots of both rice varieties with a significantly higher expression in KS282. Moreover, relative change in expression of OsGLP8-7, OsGLP8-10 and OsGLP8-11 under salt stress and OsGLP8-7 under drought stress was also commonly higher in leaves and roots of KS282 as compared to Super Basmati. Whereas, OsGLP3-7 and OsGLP8-12 after salt stress and OsGLP8-4 and OsGLP8-12 after drought stress were observed with higher relative expression in roots of Super Basmati than KS282. Importantly, the OsGLP3-6 and OsGLP4-1 from chromosome 3 and 4 respectively showed higher expression in leaves whereas most of the OsGLP genes from chromosome 8 exhibited higher expression in roots.

CONCLUSION

Overall, as a result of this comparative analysis, OsGLP genes showed both general and specific expression profiles depending upon a specific rice variety, stress condition as well as tissue type. These results will increase our understanding of role of OsGLP genes in rice crop and provide useful information for the further in-depth research on their regulatory mechanisms in response to these stress conditions.

Transcriptome analysis of two Pogostemon cablin chemotypes reveals genes related to patchouli alcohol biosynthesis

Pogostemon cablin, a medicinally and economically important perennial herb, is cultivated around the world due to its medicinal and aromatic properties. Different P. cablin cultivars exhibit different morphological traits and patchouli oil components and contents (especially patchouli alcohol (PA) and pogostone (PO)). According to the signature constituent of the leaf, P. cablin was classified into two different chemotypes, including PA-type and PO-type. To better understand the molecular mechanisms of PA biosynthesis, the transcriptomes of Chinese-cultivated P. cablin cv. PA-type “Nanxiang” (NX) and PO-type “Paixiang” (PX) were analyzed and compared with ribonucleic acid sequencing (RNA-Seq) technology. We obtained a total of 36.83 G clean bases from the two chemotypes, compared them with seven databases and revealed 45,394 annotated unigenes. Thirty-six candidate unigenes participating in the biosynthesis of PA were found in the P. cablin transcriptomes. Overall, 8,390 differentially expressed unigenes were identified between the chemotypes, including 2,467 upregulated and 5,923 downregulated unigenes. Furthermore, six and nine differentially expressed genes (DEGs) were mapped to the terpenoid backbone biosynthetic and sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. One key sesquiterpene synthase gene involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, encoding patchoulol synthase variant 1, was significantly upregulated in NX. Additionally, GC-MS analysis of the two chemotypes in this study showed that the content of PA in NX was significantly higher than that of PX, while the content of PO showed the opposite phenotype. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that the DEG expression tendency was consistent with the transcriptome sequencing results. Overall, 23 AP2/ERF, 13 bHLH, 11 MYB, 11 NAC, three Trihelix, 10 WRKY and three bZIP genes that were differentially expressed may act as regulators of terpenoid biosynthesis. Altogether, 8,314 SSRs were recognized within 6,825 unigenes, with a distribution frequency of 18.32%, among which 1,202 unigenes contained more than one SSR. The transcriptomic characteristics of the two P. cablin chemotypes are comprehensively reported in this study, and these results will contribute to a better understanding of the molecular mechanism of PA biosynthesis. Our transcriptome data also provide a valuable genetic resource for further studies on P. cablin.

Evaluation of drought resistance and transcriptome analysis for the identification of drought-responsive genes in Iris germanica

Iris germanica, a species with very high ornamental value, exhibits the strongest drought resistance among the species in the genus Iris, but the molecular mechanism underlying its drought resistance has not been evaluated. To investigate the gene expression profile changes exhibited by high-drought-resistant I. germanica under drought stress, 10 cultivars with excellent characteristics were included in pot experiments under drought stress conditions, and the changes in the chlorophyll (Chl) content, plasma membrane relative permeability (RP), and superoxide dismutase (SOD), malondialdehyde (MDA), free proline (Pro), and soluble protein (SP) levels in leaves were compared among these cultivars. Based on their drought-resistance performance, the 10 cultivars were ordered as follows: 'Little Dream' > 'Music Box' > 'X'Brassie' > 'Blood Stone' > 'Cherry Garden' > 'Memory of Harvest' > 'Immortality' > 'White and Gold' > 'Tantara' > 'Clarence'. Using the high-drought-resistant cultivar 'Little Dream' as the experimental material, cDNA libraries from leaves and rhizomes treated for 0, 6, 12, 24, and 48 h with 20% polyethylene glycol (PEG)-6000 to simulate a drought environment were sequenced using the Illumina sequencing platform. We obtained 1, 976, 033 transcripts and 743, 982 unigenes (mean length of 716 bp) through a hierarchical clustering analysis of the resulting transcriptome data. The unigenes were compared against the Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG, and gene ontology (GO) databases for functional annotation, and the gene expression levels in leaves and rhizomes were compared between the 20% PEG-6000 stress treated (6, 12, 24, and 48 h) and control (0 h) groups using DESeq2. 7849 and 24,127 differentially expressed genes (DEGs) were obtained from leaves and rhizomes, respectively. GO and KEGG enrichment analyses of the DEGs revealed significantly enriched KEGG pathways, including ribosome, photosynthesis, hormone signal transduction, starch and sucrose metabolism, synthesis of secondary metabolites, and related genes, such as heat shock proteins (HSPs), transcription factors (TFs), and active oxygen scavengers. In conclusion, we conducted the first transcriptome sequencing analysis of the I. germanica cultivar 'Little Dream' under drought stress and generated a large amount of genetic information. This study lays the foundation for further exploration of the molecular mechanisms underlying the responses of I. germanica to drought stress and provides valuable genetic resources for the breeding of drought-resistant plants.

Transcriptome analysis of bread wheat leaves in response to salt stress

Salinity is one of the main abiotic stresses limiting crop productivity. In the current study, the transcriptome of wheat leaves in an Iranian salt-tolerant cultivar (Arg) was investigated in response to salinity stress to identify salinity stress-responsive genes and mechanisms. More than 114 million reads were generated from leaf tissues by the Illumina HiSeq 2500 platform. An amount of 81.9% to 85.7% of reads could be mapped to the wheat reference genome for different samples. The data analysis led to the identification of 98819 genes, including 26700 novel transcripts. A total of 4290 differentially expressed genes (DEGs) were recognized, comprising 2346 up-regulated genes and 1944 down-regulated genes. Clustering of the DEGs utilizing Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated that transcripts associated with phenylpropanoid biosynthesis, transporters, transcription factors, hormone signal transduction, glycosyltransferases, exosome, and MAPK signaling might be involved in salt tolerance. The expression patterns of nine DEGs were investigated by quantitative real-time PCR in Arg and Moghan3 as the salt-tolerant and susceptible cultivars, respectively. The obtained results were consistent with changes in transcript abundance found by RNA-sequencing in the tolerant cultivar. The results presented here could be utilized for salt tolerance enhancement in wheat through genetic engineering or molecular breeding.

Genome-wide association mapping of sodium and potassium concentration in rice grains and shoots under alternate wetting and drying and continuously flooded irrigation

Identification of a large number of QTL and candidate genes for sodium accumulation in a field grown population of rice derived from the aus subpopulation. Rice (Oryza sativa L.) is a globally important cereal crop. Sodium (Na⁺) and potassium (K⁺) are the major monovalent ions which affect rice growth, and exploring their uptake mechanisms will be useful for understanding rice biology. Since the balance of Na⁺ and K⁺ plays a significant role in adaptation of rice to salinity, that biology might inform the search for tolerance. In this study, the Na⁺ and K⁺ concentration and Na⁺/K⁺ ratio in grains and shoots were analyzed in the Bengal and Assam Aus Panel grown in field conditions under continuously flooded (CF) and alternate wetting and drying (AWD) irrigation. Overall, AWD irrigation significantly reduced the Na⁺ concentration and increased the K⁺ concentration in shoots and grains compared to the plants grown under CF. Genome-wide association mapping was conducted on Na⁺, K⁺ concentration and Na⁺/K⁺ ratio with 2 million SNPs using an efficient mixed model. Only QTLs which contained more than two significant SNPs (p < 0.0001) and where at least one of these significant SNPs passed a 10% false discovery rate were reported. A total of 106 QTLs were identified as being associated with Na⁺ concentration and Na⁺/K⁺ ratio across all traits and field conditions, with 48 QTLs found in multiple traits and/or water conditions. Four notable QTLs (one each on chromosomes 1 and 11, two on chromosome 2) and the haplotype variants of four candidate genes (OsHKT1;5, OsNHX2, LOC_Os02g32490 and OsFAD2_1) are discussed. The QTLs/candidate genes identified here could be useful for breeding rice that accumulates lower concentrations of sodium.

In Response to Abiotic Stress, DNA Methylation Confers EpiGenetic Changes in Plants

Epigenetics involves the heritable changes in patterns of gene expression determined by developmental and abiotic stresses, i.e., drought, cold, salinity, trace metals, and heat. Gene expression is driven by changes in DNA bases, histone proteins, the biogenesis of ncRNA, and changes in the nucleotide sequence. To cope with abiotic stresses, plants adopt certain changes driven by a sophisticated biological system. DNA methylation is a primary mechanism for epigenetic variation, which can induce phenotypic alterations in plants under stress. Some of the stress-driven changes in plants are temporary, while some modifications may be stable and inheritable to the next generations to allow them to cope with such extreme stress challenges in the future. In this review, we discuss the pivotal role of epigenetically developed phenotypic characteristics in plants as an evolutionary process participating in adaptation and tolerance responses to abiotic and biotic stresses that alter their growth and development. We emphasize the molecular process underlying changes in DNA methylation, differential variation for different species, the roles of non-coding RNAs in epigenetic modification, techniques for studying DNA methylation, and its role in crop improvement in tolerance to abiotic stress (drought, salinity, and heat). We summarize DNA methylation as a significant future research priority for tailoring crops according to various challenging environmental issues.

Genome-Wide Analysis of Potassium Channel Genes in Rice: Expression of the OsAKT and OsKAT Genes under Salt Stress

Potassium (K+), as a vital element, is involved in regulating important cellular processes such as enzyme activity, cell turgor, and nutrient movement in plant cells, which affects plant growth and production. Potassium channels are involved in the transport and release of potassium in plant cells. In the current study, three OsKAT genes and two OsAKT genes, along with 11 nonredundant putative potassium channel genes in the rice genome, were characterized based on their physiochemical properties, protein structure, evolution, duplication, in silico gene expression, and protein-protein interactions. In addition, the expression patterns of OsAKTs and OsKATs were studied in root and shoot tissues under salt stress using real-time PCR in three rice cultivars. K+ channel genes were found to have diverse functions and structures, and OsKATs showed high genetic divergence from other K+ channel genes. Furthermore, the Ka/Ks ratios of duplicated gene pairs from the K+ channel gene family in rice suggested that these genes underwent purifying selection. Among the studied K+ channel proteins, OsKAT1 and OsAKT1 were identified as proteins with high potential N-glycosylation and phosphorylation sites, and LEU, VAL, SER, PRO, HIS, GLY, LYS, TYR, CYC, and ARG amino acids were predicted as the binding residues in the ligand-binding sites of K+ channel proteins. Regarding the coexpression network and KEGG ontology results, several metabolic pathways, including sugar metabolism, purine metabolism, carbon metabolism, glycerophospholipid metabolism, monoterpenoid biosynthesis, and folate biosynthesis, were recognized in the coexpression network of K+ channel proteins. Based on the available RNA-seq data, the K+ channel genes showed differential expression levels in rice tissues in response to biotic and abiotic stresses. In addition, the real-time PCR results revealed that OsAKTs and OsKATs are induced by salt stress in root and shoot tissues of rice cultivars, and OsKAT1 was identified as a key gene involved in the rice response to salt stress. In the present study, we found that the repression of OsAKTs, OsKAT2, and OsKAT2 in roots was related to salinity tolerance in rice. Our findings provide valuable insights for further structural and functional assays of K+ channel genes in rice.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.

DNA hypomethylation in tetraploid rice potentiates stress-responsive gene expression for salt tolerance

Polyploidy is a prominent feature for genome evolution in many animals and all flowering plants. Plant polyploids often show enhanced fitness in diverse and extreme environments, but the molecular basis for this remains elusive. Soil salinity presents challenges for many plants including agricultural crops. Here we report that salt tolerance is enhanced in tetraploid rice through lower sodium uptake and correlates with epigenetic regulation of jasmonic acid (JA)-related genes. Polyploidy induces DNA hypomethylation and potentiates genomic loci coexistent with many stress-responsive genes, which are generally associated with proximal transposable elements (TEs). Under salt stress, the stress-responsive genes including those in the JA pathway are more rapidly induced and expressed at higher levels in tetraploid than in diploid rice, which is concurrent with increased jasmonoyl isoleucine (JA-Ile) content and JA signaling to confer stress tolerance. After stress, elevated expression of stress-responsive genes in tetraploid rice can induce hypermethylation and suppression of the TEs adjacent to stress-responsive genes. These induced responses are reproducible in a recurring round of salt stress and shared between two japonica tetraploid rice lines. The data collectively suggest a feedback relationship between polyploidy-induced hypomethylation in rapid and strong stress response and stress-induced hypermethylation to repress proximal TEs and/or TE-associated stress-responsive genes. This feedback regulation may provide a molecular basis for selection to enhance adaptation of polyploid plants and crops during evolution and domestication.

A Comparison of the Transcriptomes of Cowpeas in Response to Two Different Ionizing Radiations

In this study, gene expression changes in cowpea plants irradiated by two different types of radiation: proton-beams and gamma-rays were investigated. Seeds of the Okdang cultivar were exposed to 100, 200, and 300 Gy of gamma-rays and proton-beams. In transcriptome analysis, the 32, 75, and 69 differentially expressed genes (DEGs) at each dose of gamma-ray irradiation compared with that of the control were identified. A total of eight genes were commonly up-regulated for all gamma-ray doses. However, there were no down-regulated genes. In contrast, 168, 434, and 387 DEGs were identified for each dose of proton-beam irradiation compared with that of the control. A total of 61 DEGs were commonly up-regulated for all proton-beam doses. As a result of GO and KEGG analysis, the ranks of functional categories according to the number of DEGs were not the same in both treatments and were more diverse in terms of pathways in the proton-beam treatments than gamma-ray treatments. The number of genes related to defense, photosynthesis, reactive oxygen species (ROS), plant hormones, and transcription factors (TF) that were up-/down-regulated was higher in the proton beam treatment than that in gamma ray treatment. Proton-beam treatment had a distinct mutation spectrum and gene expression pattern compared to that of gamma-ray treatment. These results provide important information on the mechanism for gene regulation in response to two ionizing radiations in cowpeas.