Global Transcriptome Analysis of Combined Abiotic Stress Signaling Genes Unravels Key Players in Oryza sativa L.: An In silico Approach

Physiological and metabolic analyses reveal the proline-mediated flowering delay mechanism in Prunus persica

Peaches are susceptible to various environmental stresses. Particularly in late spring, freezing temperatures can damage peaches and consequently, affect their productivity. Therefore, flowering delay is a prominent strategy for avoiding spring frost damage. Our previous study confirmed that treatment with 5% sodium alginate and 100 mM CaCl2 (5AG) to avoid frost damage during the blooming stage delays flowering. To reveal the flowering delay mechanism of peaches, this study systematically analyzed the modification of amino acid profiles in control and 5AG-treated peach plants at different day intervals. Our findings indicate that arginine (Arg), glutamate (Glu), and proline (Pro) levels differed between the control and 5AG-treated peach shoots throughout the phenological development of flower buds. Furthermore, two amino acids (Arg and Glu) are involved in the Pro pathway. Thus, using a computational metabolomics method, Pro biosynthesis and its characteristics, gene ontology, gene synteny, cis-regulatory elements, and gene organizations were examined to decipher the involvement of Pro metabolism in peach flowering delay. In addition, qRT-PCR analysis revealed the transcriptional regulation of Pro-related and flowering-responsive genes and their role in flowering delay. Overall, this pilot study provides new insights into the role of Pro in the flowering delay mechanisms in Prunus persica through 5AG treatment.

Identification and virus-induced gene silencing (VIGS) analysis of methyltransferase affecting tomato (Solanum lycopersicum) fruit ripening

MAIN CONCLUSION

Six methyltransferase genes affecting tomato fruit ripening were identified through genome-wide screening, VIGS assay, and expression pattern analysis. The data provide the basis for understanding new mechanisms of methyltransferases. Fruit ripening is a critical stage for the formation of edible quality and seed maturation, which is finely modulated by kinds of factors, including genetic regulators, hormones, external signals, etc. Methyltransferases (MTases), important genetic regulators, play vital roles in plant development through epigenetic regulation, post-translational modification, or other mechanisms. However, the regulatory functions of numerous MTases except DNA methylation in fruit ripening remain limited so far. Here, six MTases, which act on different types of substrates, were identified to affect tomato fruit ripening. First, 35 MTase genes with relatively high expression at breaker (Br) stage of tomato fruit were screened from the tomato MTase gene database encompassing 421 genes totally. Thereafter, six MTase genes were identified as potential regulators of fruit ripening via virus-induced gene silencing (VIGS), including four genes with a positive regulatory role and two genes with a negative regulatory role, respectively. The expression of these six MTase genes exhibited diverse patterns during the fruit ripening process, and responded to various external ripening-related factors, including ethylene, 1-methylcyclopropene (1-MCP), temperature, and light exposure. These results help to further elaborate the biological mechanisms of MTase genes in tomato fruit ripening and enrich the understanding of the regulatory mechanisms of fruit ripening involving MTases, despite of DNA MTases.

Transcriptional survey of abiotic stress response in maize (Zea mays) in the level of gene co-expression network and differential gene correlation analysis

Abstract. Maize may be exposed to several abiotic stresses in the field. Therefore, identifying the tolerance mechanisms of natural field stress is mandatory. Gene expression data of maize upon abiotic stress were collected, and 560 differentially expressed genes (DEGs) were identified through meta-analysis. The most significant gene ontology terms in up-regulated genes were 'response to abiotic stress' and 'chitinase activity'. 'Phosphorelay signal transduction system' was the most significant enriched biological process in down-regulated DEGs. The co-expression analysis unveiled seven modules of DEGs, with a notable positive correlation between the modules and abiotic stress. Furthermore, the statistical significance was strikingly high for the turquoise, green and yellow modules. The turquoise group played a central role in orchestrating crucial adaptations in metabolic and stress response pathways in maize when exposed to abiotic stress. Within three up-regulated modules, Zm.7361.1.A1_at, Zm.10386.1.A1_a_at and Zm.10151.1.A1_at emerged as hub genes. These genes might introduce novel candidates implicated in stress tolerance mechanisms, warranting further comprehensive investigation and research. In parallel, the R package glmnet was applied to fit a logistic LASSO regression model on the DEGs profile to select candidate genes associated with abiotic responses in maize. The identified hub genes and LASSO regression genes were validated on an independent microarray dataset. Additionally, Differential Gene Correlation Analysis (DGCA) was performed on LASSO and hub genes to investigate the gene-gene regulatory relationship. The P value of DGCA of 16 pairwise gene comparisons was lower than 0.01, indicating a gene-gene significant change in correlation between control and abiotic stress. Integrated weighted gene correlation network analysis and logistic LASSO analysis revealed Zm.11185.1.S1_at, Zm.2331.1.S1_x_at and Zm.17003.1.S1_at. Notably, these 3 genes were identified in the 16 gene-pair comparisons. This finding highlights the notable significance of these genes in the abiotic stress response. Additional research into maize stress tolerance may focus on these three genes.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

Long-term cold, freezing and drought: overlapping and specific regulatory mechanisms and signal transduction in tea plant (Camellia sinensis (L.) Kuntze)

Introduction

Low temperatures and drought are two main environmental constraints reducing the yield and geographical distribution of horticultural crops worldwide. Understanding the genetic crosstalk between stress responses has potential importance for crop improvement.

Methods

In this study, Illumina RNA-seq and Pac-Bio genome resequencing were used to annotate genes and analyze transcriptome dynamics in tea plants under long-term cold, freezing, and drought.

Results

The highest number of differentially expressed genes (DEGs) was identified under long-term cold (7,896) and freezing (7,915), with 3,532 and 3,780 upregulated genes, respectively. The lowest number of DEGs was observed under 3-day drought (47) and 9-day drought (220), with five and 112 genes upregulated, respectively. The recovery after the cold had 6.5 times greater DEG numbers as compared to the drought recovery. Only 17.9% of cold-induced genes were upregulated by drought. In total, 1,492 transcription factor genes related to 57 families were identified. However, only 20 transcription factor genes were commonly upregulated by cold, freezing, and drought. Among the 232 common upregulated DEGs, most were related to signal transduction, cell wall remodeling, and lipid metabolism. Co-expression analysis and network reconstruction showed 19 genes with the highest co-expression connectivity: seven genes are related to cell wall remodeling (GATL7, UXS4, PRP-F1, 4CL, UEL-1, UDP-Arap, and TBL32), four genes are related to calcium-signaling (PXL1, Strap, CRT, and CIPK6), three genes are related to photo-perception (GIL1, CHUP1, and DnaJ11), two genes are related to hormone signaling (TTL3 and GID1C-like), two genes are involved in ROS signaling (ERO1 and CXE11), and one gene is related to the phenylpropanoid pathway (GALT6).

Discussion

Based on our results, several important overlapping mechanisms of long-term stress responses include cell wall remodeling through lignin biosynthesis, o-acetylation of polysaccharides, pectin biosynthesis and branching, and xyloglucan and arabinogalactan biosynthesis. This study provides new insight into long-term stress responses in woody crops, and a set of new target candidate genes were identified for molecular breeding aimed at tolerance to abiotic stresses.

Transcriptome meta-analysis of abiotic stresses-responsive genes and identification of candidate transcription factors for broad stress tolerance in wheat

Under field conditions, wheat is subjected to single or multiple stress conditions. The elucidation of the molecular mechanism of stress response is a key step to identify candidate genes for stress resistance in plants. In this study, RNA-seq data analysis identified 17.324, 10.562, 5.510, and 8.653 differentially expressed genes (DEGs) under salt, drought, heat, and cold stress, respectively. Moreover, the comparison of DEGs from each stress revealed 2374 shared genes from which 40% showed highly conserved expression patterns. Moreover, co-expression network analysis and GO enrichment revealed co-expression modules enriched with genes involved in transcription regulation, protein kinase pathway, and genes responding to phytohormones or modulating hormone levels. The expression of 15 selected transcription factor encoding genes was analyzed under abiotic stresses and ABA treatment in durum wheat. The identified transcription factor genes are excellent candidates for genetic engineering of stress tolerance in wheat.

Genome-Wide Identification of Potential mRNAs in Drought Response in Wheat (Triticum aestivum L.)

Plant cell metabolism inevitably forms an important drought-responsive mechanism, which halts crop productivity. Globally, more than 30% of the total harvested area was affected by dehydration. RNA-seq technology has enabled biologists to identify stress-responsive genes in relatively quick times. However, one shortcoming of this technology is the inconsistent data generation compared to other parts of the world. So, we have tried, here, to generate a consensus by analyzing meta-transcriptomic data available in the public microarray database GEO NCBI. In this way, the aim was set, here, to identify stress genes commonly identified as differentially expressed (p < 0.05) then followed by downstream analyses. The search term “Drought in wheat” resulted in 233 microarray experiments from the GEO NCBI database. After discarding empty datasets containing no expression data, the large-scale meta-transcriptome analytics and one sample proportional test were carried out (Bonferroni adjusted p < 0.05) to reveal a set of 11 drought-responsive genes on a global scale. The annotation of these genes revealed that the transcription factor activity of RNA polymerase II and sequence-specific DNA-binding mechanism had a significant role during the drought response in wheat. Similarly, the primary root differentiation zone annotations, controlled by TraesCS5A02G456300 and TraesCS7B02G243600 genes, were found as top-enriched terms (p < 0.05 and Q < 0.05). The resultant standard drought genes, glycosyltransferase; Arabidopsis thaliana KNOTTED-like; bHLH family protein; Probable helicase MAGATAMA 3; SBP family protein; Cytochrome c oxidase subunit 2; Trihelix family protein; Mic1 domain-containing protein; ERF family protein; HD-ZIP I protein; and ERF family protein, are important in terms of their worldwide proved link with stress. From a future perspective, this study could be important in a breeding program contributing to increased crop yield. Moreover, the wheat varieties could be identified as drought-resistant/sensitive based on the nature of gene expression levels.

Structural and Functional Characterization at the Molecular Level of the MATE Gene Family in Wheat in Silico

A series of multidrug extransporters known as the multidrug and potentially toxic extrusion (MATE) genes are found in all living things and are crucial for the removal of heavy metal ions, metalloids, exogenous xenobiotics, endogenous secondary metabolites, and other toxic substances from the cells. However, there has only been a small amount of them in silico analysis of the MATE family of genes in plant species. In the current study, the MATE gene family was characterized in silico where two families and seven subfamilies based on their evolutionary relationships were proposed. Plant breeders may use TraesCS1D02G030400, TraesCS4B02G244400, and TraesCS1A02G029900 genes for marker-assisted or transgenic breeding to develop novel cultivars since these genes have been hypothesized from protein-protein interaction study to play a critical role in the transport of toxic chemicals across cells. The exon number varies from 01 to 14. One exon has TraesCS1A02G188100, TraesCS5B02G562500, TraesCS6A02G256400, and TraesCS6D02G384300 genes, while 14 exons have only two genes that are TraesCS6A02G418800 and TraesCS6D02G407900. Biological stress (infestations of disease) affects the expression of most of the MATE genes, with the gene TraesCS5D02G355500 having the highest expression level in the wheat expression browser tool. Using the Grain interpretation search engine tool, it is found that the vast bulk of MATE genes are voiced throughout biotic environmental stresses caused by disease pests, with the genotype TraesCS5B02G326600.1 from family 1 exhibiting the greatest level of expression throughout Fusarium head blight infection by Fusarium graminearum after 4 days of infection. The researchers constructed 39 ternary plots, each with a distinct degree of expression under biotic and abiotic stress settings, and observed that 44% of the triplets have imbalanced outputs (extreme values) due to their higher tissue specificity and increased intensity.

Global Integrated Genomic and Transcriptomic Analyses of MYB Transcription Factor Superfamily in C3 Model Plant Oryza sativa (L.) Unravel Potential Candidates Involved in Abiotic Stress Signaling

Plant transcription factors (TFs) are significant players in transcriptional regulations, signal transduction, and constitute an integral part of signaling networks. MYB TFs are major TF superfamilies that play pivotal roles in regulation of transcriptional reprogramming, physiological processes, and abiotic stress (AbS) responses. To explore the understanding of MYB TFs, genome and transcriptome-wide identification was performed in the C3 model plant, Oryza sativa (OsMYB). This study retrieved 114 OsMYB TFs that were computationally analyzed for their expression profiling, gene organization, cis-acting elements, and physicochemical properties. Based on the microarray datasets, six OsMYB genes which were sorted out and identified by a differential expression pattern were noted in various tissues. Systematic expression profiling of OsMYB TFs showed their meta-differential expression of different AbS treatments, spatio-temporal gene expression of various tissues and their growth in the field, and gene expression profiling in responses to phytohormones. In addition, the circular ideogram of OsMYB genes in related C4 grass plants conferred the gene synteny. Protein–protein interactions of these genes revealed the molecular crosstalk of OsMYB TFs. Transcriptional analysis (qPCR) of six OsMYB players in response to drought and salinity stress suggested the involvement in individual and combined AbS responses. To decipher how these OsMYB play functional roles in AbS dynamics, further research is a prerequisite.

Multi-Omics and Integrative Approach towards Understanding Salinity Tolerance in Rice: A Review

Rice (Oryza sativa L.) plants are simultaneously encountered by environmental stressors, most importantly salinity stress. Salinity is the major hurdle that can negatively impact growth and crop yield. Understanding the salt stress and its associated complex trait mechanisms for enhancing salt tolerance in rice plants would ensure future food security. The main aim of this review is to provide insights and impacts of molecular-physiological responses, biochemical alterations, and plant hormonal signal transduction pathways in rice under saline stress. Furthermore, the review highlights the emerging breakthrough in multi-omics and computational biology in identifying the saline stress-responsive candidate genes and transcription factors (TFs). In addition, the review also summarizes the biotechnological tools, genetic engineering, breeding, and agricultural practicing factors that can be implemented to realize the bottlenecks and opportunities to enhance salt tolerance and develop salinity tolerant rice varieties. Future studies pinpointed the augmentation of powerful tools to dissect the salinity stress-related novel players, reveal in-depth mechanisms and ways to incorporate the available literature, and recent advancements to throw more light on salinity responsive transduction pathways in plants. Particularly, this review unravels the whole picture of salinity stress tolerance in rice by expanding knowledge that focuses on molecular aspects.

Meta-analysis of transcriptomic responses to cold stress in plants

Transcriptomic analyses are needful tools to gain insight into the molecular mechanisms underlying plant responses to abiotic stresses. The aim of this study was to identify key genes differentially regulated in response to chilling stress in various plant species with different levels of tolerance to low temperatures. A meta-analysis was performed using the RNA-Seq data of published studies whose experimental conditions were comparable. The results confirmed the importance of ethylene in the hormonal cross-talk modulating the defensive responses against chilling stress, especially in sensitive species. The transcriptomic activity of five Ethylene Response Factors genes and a REDOX Responsive Transcription Factor 1 involved in hormone-related pathways belonging to ethylene metabolism and signal transduction were induced. Transcription activity of two genes encoding for heat shock factors was enhanced, together with various genes associated with developmental processes. Several transcription factor families showed to be commonly induced between different plant species. Protein-protein interaction networks highlighted the role of the photosystems I and II, as well as genes encoding for HSF and WRKY transcription factors. A model of gene regulatory network underlying plant responses to chilling stress was developed, allowing the delivery of new candidate genes for genetic improvement of crops towards low temperatures tolerance.

Global analysis of switchgrass (Panicum virgatum L.) transcriptomes in response to interactive effects of drought and heat stresses

BACKGROUND

Sustainable production of high-quality feedstock has been of great interest in bioenergy research. Despite the economic importance, high temperatures and water deficit are limiting factors for the successful cultivation of switchgrass in semi-arid areas. There are limited reports on the molecular basis of combined abiotic stress tolerance in switchgrass, particularly the combination of drought and heat stress. We used transcriptomic approaches to elucidate the changes in the response of switchgrass to drought and high temperature simultaneously.

RESULTS

We conducted solely drought treatment in switchgrass plant Alamo AP13 by withholding water after 45 days of growing. For the combination of drought and heat effect, heat treatment (35 °C/25 °C day/night) was imposed after 72 h of the initiation of drought. Samples were collected at 0 h, 72 h, 96 h, 120 h, 144 h, and 168 h after treatment imposition, total RNA was extracted, and RNA-Seq conducted. Out of a total of 32,190 genes, we identified 3912, as drought (DT) responsive genes, 2339 and 4635 as, heat (HT) and drought and heat (DTHT) responsive genes, respectively. There were 209, 106, and 220 transcription factors (TFs) differentially expressed under DT, HT and DTHT respectively. Gene ontology annotation identified the metabolic process as the significant term enriched in DTHT genes. Other biological processes identified in DTHT responsive genes included: response to water, photosynthesis, oxidation-reduction processes, and response to stress. KEGG pathway enrichment analysis on DT and DTHT responsive genes revealed that TFs and genes controlling phenylpropanoid pathways were important for individual as well as combined stress response. For example, hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) from the phenylpropanoid pathway was induced by single DT and combinations of DTHT stress.

CONCLUSION

Through RNA-Seq analysis, we have identified unique and overlapping genes in response to DT and combined DTHT stress in switchgrass. The combination of DT and HT stress may affect the photosynthetic machinery and phenylpropanoid pathway of switchgrass which negatively impacts lignin synthesis and biomass production of switchgrass. The biological function of genes identified particularly in response to DTHT stress could further be confirmed by techniques such as single point mutation or RNAi.

Comparative transcriptome analysis of the cold resistance of the sterile rice line 33S

Rice (Oryza sativa L.) is one of the most important species for food production worldwide. Low temperature is a major abiotic factor that affects rice germination and reproduction. Here, the underlying regulatory mechanism in seedlings of a TGMS variety (33S) and a cold-sensitive variety (Nipponbare) was investigated by comparative transcriptome. There were 795 differentially expressed genes (DEGs) identified only in cold-treated 33S, suggesting that 33S had a unique cold-resistance system. Functional and enrichment analysis of these DEGs revealed that, in 33S, several metabolic pathways, such as photosynthesis, amino acid metabolism, secondary metabolite biosynthesis, were significantly repressed. Moreover, pathways related to growth and development, including starch and sucrose metabolism, and DNA biosynthesis and damage response/repair, were significantly enhanced. The expression of genes related to nutrient reserve activity were significantly up-regulated in 33S. Finally, three NAC and several ERF transcription factors were predicted to be important in this transcriptional reprogramming. This present work provides valuable information for future investigations of low-temperature response mechanisms and genetic improvement of cold-tolerant rice seedlings.

Inferring the genetic responses to acute drought stress across an ecological gradient

BACKGROUND

How do xerophytic species thrive in environments that experience extreme annual drought? Although critical to the survival of many species, the genetic responses to drought stress in many non-model organisms has yet to be explored. We investigated this question in Mentzelia section Bartonia (Loasaceae), which occurs throughout western North America, including arid lands. To better understand the genetic responses to drought stress among species that occur in different habitats, the gene expression levels of three species from Mentzelia were compared across a precipitation gradient. Two de novo reference transcriptomes were generated and annotated. Leaf and root tissues were collected from control and drought shocked plants and compared to one another for differential expression. A target-gene approach was also implemented to better understand how drought-related genes from model and crop species function in non-model systems.

RESULTS

When comparing the drought-shock treatment plants to their respective control plants, we identified 165 differentially expressed clusters across all three species. Differentially expressed genes including those associated with water movement, photosynthesis, and delayed senescence. The transcriptome profiling approach was coupled with a target genes approach that measured expression of 90 genes associated with drought tolerance in model organisms. Comparing differentially expressed genes with a ≥ 2 log-fold value between species and tissue types showed significant differences in drought response. In pairwise comparisons, species that occurred in drier environments differentially expressed greater genes in leaves when drought shocked than those from wetter environments, but expression in the roots mostly produced opposite results.

CONCLUSIONS

Arid-adapted species mount greater genetic responses compared to the mesophytic species, which has likely evolved in response to consistent annual drought exposure across generations. Drought responses also depended on organ type. Xerophytes, for example, mounted a larger response in leaves to downregulate photosynthesis and senescence, while mobilizing carbon and regulating water in the roots. The complexity of drought responses in Mentzelia suggest that whole organism responses need to be considered when studying drought and, in particular, the physiological mechanisms in which plants regulate water, carbon, cell death, metabolism, and secondary metabolites.

A Comparative Transcriptomic Meta-Analysis Revealed Conserved Key Genes and Regulatory Networks Involved in Drought Tolerance in Cereal Crops

Drought affects plant growth and development, causing severe yield losses, especially in cereal crops. The identification of genes involved in drought tolerance is crucial for the development of drought-tolerant crops. The aim of this study was to identify genes that are conserved key players for conferring drought tolerance in cereals. By comparing the transcriptomic changes between tolerant and susceptible genotypes in four Gramineae species, we identified 69 conserved drought tolerant-related (CDT) genes that are potentially involved in the drought tolerance of all of the analysed species. The CDT genes are principally involved in stress response, photosynthesis, chlorophyll biogenesis, secondary metabolism, jasmonic acid signalling, and cellular transport. Twenty CDT genes are not yet characterized and can be novel candidates for drought tolerance. The k-means clustering analysis of expression data highlighted the prominent roles of photosynthesis and leaf senescence-related mechanisms in differentiating the drought response between tolerant and sensitive genotypes. In addition, we identified specific transcription factors that could regulate the expression of photosynthesis and leaf senescence-related genes. Our analysis suggests that the balance between the induction of leaf senescence and maintenance of photosynthesis during drought plays a major role in tolerance. Fine-tuning of CDT gene expression modulation by specific transcription factors can be the key to improving drought tolerance in cereals.

Elucidating the Response of Crop Plants towards Individual, Combined and Sequentially Occurring Abiotic Stresses

In nature, plants are exposed to an ever-changing environment with increasing frequencies of multiple abiotic stresses. These abiotic stresses act either in combination or sequentially, thereby driving vegetation dynamics and limiting plant growth and productivity worldwide. Plants' responses against these combined and sequential stresses clearly differ from that triggered by an individual stress. Until now, experimental studies were mainly focused on plant responses to individual stress, but have overlooked the complex stress response generated in plants against combined or sequential abiotic stresses, as well as their interaction with each other. However, recent studies have demonstrated that the combined and sequential abiotic stresses overlap with respect to the central nodes of their interacting signaling pathways, and their impact cannot be modelled by swimming in an individual extreme event. Taken together, deciphering the regulatory networks operative between various abiotic stresses in agronomically important crops will contribute towards designing strategies for the development of plants with tolerance to multiple stress combinations. This review provides a brief overview of the recent developments in the interactive effects of combined and sequentially occurring stresses on crop plants. We believe that this study may improve our understanding of the molecular and physiological mechanisms in untangling the combined stress tolerance in plants, and may also provide a promising venue for agronomists, physiologists, as well as molecular biologists.

Elucidating the Response of Crop Plants towards Individual, Combined and Sequentially Occurring Abiotic Stresses

In nature, plants are exposed to an ever-changing environment with increasing frequencies of multiple abiotic stresses. These abiotic stresses act either in combination or sequentially, thereby driving vegetation dynamics and limiting plant growth and productivity worldwide. Plants' responses against these combined and sequential stresses clearly differ from that triggered by an individual stress. Until now, experimental studies were mainly focused on plant responses to individual stress, but have overlooked the complex stress response generated in plants against combined or sequential abiotic stresses, as well as their interaction with each other. However, recent studies have demonstrated that the combined and sequential abiotic stresses overlap with respect to the central nodes of their interacting signaling pathways, and their impact cannot be modelled by swimming in an individual extreme event. Taken together, deciphering the regulatory networks operative between various abiotic stresses in agronomically important crops will contribute towards designing strategies for the development of plants with tolerance to multiple stress combinations. This review provides a brief overview of the recent developments in the interactive effects of combined and sequentially occurring stresses on crop plants. We believe that this study may improve our understanding of the molecular and physiological mechanisms in untangling the combined stress tolerance in plants, and may also provide a promising venue for agronomists, physiologists, as well as molecular biologists.

Meta-Analysis of Yield-Related and N-Responsive Genes Reveals Chromosomal Hotspots, Key Processes and Candidate Genes for Nitrogen-Use Efficiency in Rice

Nitrogen-use efficiency (NUE) is a function of N-response and yield that is controlled by many genes and phenotypic parameters that are poorly characterized. This study compiled all known yield-related genes in rice and mined them from the N-responsive microarray data to find 1,064 NUE-related genes. Many of them are novel genes hitherto unreported as related to NUE, including 80 transporters, 235 transcription factors (TFs), 44 MicroRNAs (miRNAs), 91 kinases, and 8 phosphatases. They were further shortlisted to 62 NUE-candidate genes following hierarchical methods, including quantitative trait locus (QTL) co-localization, functional evaluation in the literature, and protein-protein interactions (PPIs). They were localized to chromosomes 1, 3, 5, and 9, of which chromosome 1 with 26 genes emerged as a hotspot for NUE spanning 81% of the chromosomes. Further, co-localization of the NUE genes on NUE-QTLs resolved differences in the earlier studies that relied mainly on N-responsive genes regardless of their role in yield. Functional annotations and PPIs for all the 1,064 NUE-related genes and also the shortlisted 62 candidates revealed transcription, redox, phosphorylation, transport, development, metabolism, photosynthesis, water deprivation, and hormonal and stomatal function among the prominent processes. In silico expression analysis confirmed differential expression of the 62 NUE-candidate genes in a tissue/stage-specific manner. Experimental validation in two contrasting genotypes revealed that high NUE rice shows better photosynthetic performance, transpiration efficiency and internal water-use efficiency in comparison to low NUE rice. Feature Selection Analysis independently identified one-third of the common genes at every stage of hierarchical shortlisting, offering 6 priority targets to validate for improving the crop NUE.

Comparative Expression Analysis of Stress-Inducible Candidate Genes in Response to Cold and Drought in Tea Plant [Camellia sinensis (L.) Kuntze]

Cold and drought are two of the most severe threats affecting the growth and productivity of the tea plant, limiting its global spread. Both stresses cause osmotic changes in the cells of the tea plant by decreasing their water potential. To develop cultivars that are tolerant to both stresses, it is essential to understand the genetic responses of tea plant to these two stresses, particularly in terms of the genes involved. In this study, we combined literature data with interspecific transcriptomic analyses (using Arabidopsis thaliana and Solanum lycopersicum) to choose genes related to cold tolerance. We identified 45 stress-inducible candidate genes associated with cold and drought responses in tea plants based on a comprehensive homologous detection method. Of these, nine were newly characterized by us, and 36 had previously been reported. The gene network analysis revealed upregulated expression in ICE1-related cluster of bHLH factors, HSP70/BAM5 connected genes (hexokinases, galactinol synthases, SnRK complex, etc.) indicating their possible co-expression. Using qRT-PCR we revealed that 10 genes were significantly upregulated in response to both cold and drought in tea plant: HSP70, GST, SUS1, DHN1, BMY5, bHLH102, GR-RBP3, ICE1, GOLS1, and GOLS3. SnRK1.2, HXK1/2, bHLH7/43/79/93 were specifically upregulated in cold, while RHL41, CAU1, Hydrolase22 were specifically upregulated in drought. Interestingly, the expression of CIP was higher in the recovery stage of both stresses, indicating its potentially important role in plant recovery after stress. In addition, some genes, such as DHN3, bHLH79, PEI54, SnRK1.2, SnRK1.3, and Hydrolase22, were significantly positively correlated between the cold and drought responses. CBF1, GOLS1, HXK2, and HXK3, by contrast, showed significantly negative correlations between the cold and drought responses. Our results provide valuable information and robust candidate genes for future functional analyses intended to improve the stress tolerance of the tea plant and other species.

Identification of conserved genes linked to responses to abiotic stresses in leaves among different plant species

As a consequence of global climate change, certain stress factors that have a negative impact on crop productivity such as heat, cold, drought and salinity are becoming increasingly prevalent. We conducted a meta-analysis to identify genes conserved across plant species involved in (1) general abiotic stress conditions, and (2) specific and unique abiotic stress factors (drought, salinity, extreme temperature) in leaf tissues. We collected raw data and re-analysed eight RNA-Seq studies using our previously published bioinformatic pipeline. A total of 68 samples were analysed. Gene set enrichment analysis was performed using MapMan and PageMan whereas DAVID (Database for Annotation, Visualisation and Integrated Discovery) was used for metabolic process enrichment analysis. We identified of a total of 5122 differentially expressed genes when considering all abiotic stresses (3895 were upregulated and 1227 were downregulated). Jasmonate-related genes were more commonly upregulated by drought, whereas gibberellin downregulation was a key signal for drought and heat. In contrast, cold stress clearly upregulated genes involved in ABA (abscisic acid), cytokinin and gibberellins. A gene (non-phototrophic hypocotyl) involved in IAA (indoleacetic acid) response was induced by heat. Regarding secondary metabolism, as expected, MVA pathway (mevalonate pathway), terpenoids and alkaloids were generally upregulated by all different stresses. However, flavonoids, lignin and lignans were more repressed by heat (cinnamoyl coA reductase 1 and isopentenyl pyrophosphatase). Cold stress drastically modulated genes involved in terpenoid and alkaloids. Relating to transcription factors, AP2-EREBP, MADS-box, WRKY22, MYB, homoebox genes members were significantly modulated by drought stress whereas cold stress enhanced AP2-EREBPs, bZIP members, MYB7, BELL 1 and one bHLH member. C2C2-CO-LIKE, MADS-box and a homeobox (HOMEOBOX3) were mostly repressed in response to heat. Gene set enrichment analysis showed that ubiquitin-mediated protein degradation was enhanced by heat, which unexpectedly repressed glutaredoxin genes. Cold stress mostly upregulated MAP kinases (mitogen-activated protein kinase). Findings of this work will allow the identification of new molecular markers conserved across crops linked to major genes involved in quantitative agronomic traits affected by different abiotic stress.

Computational solutions for modeling and controlling plant response to abiotic stresses: a review with focus on iron deficiency

Computational solutions enable plant scientists to model protein-mediated stress responses and characterize novel gene functions that coordinate responses to a variety of abiotic stress conditions. Recently, density functional theory was used to study proteins active sites and elucidate enzyme conversion mechanisms involved in iron deficiency responsive signaling pathways. Computational approaches for protein homology modeling and the kinetic modeling of signaling pathways have also resolved the identity and function in proteins involved in iron deficiency signaling pathways. Significant changes in gene relationships under other stress conditions, such as heat or drought stress, have been recently identified using differential network analysis, suggesting that stress tolerance is achieved through asynchronous control. Moreover, the increasing development and use of statistical modeling and systematic modeling of transcriptomic data have provided significant insight into the gene regulatory mechanisms associated with abiotic stress responses. These types of in silico approaches have facilitated the plant science community's future goals of developing multi-scale models of responses to iron deficiency stress and other abiotic stress conditions.

Integrated transcriptomic and metabolomic analyses of glutamine metabolism genes unveil key players in Oryza sativa (L.) to ameliorate the unique and combined abiotic stress tolerance

Plants can be considered to biosynthesize the specialized metabolites to adapt to various environmental stressors mainly on abiotic stresses (AbS). Among specialized metabolites, glutamine (Gln) is an essential plant metabolite to achieve sustainable plant growth, yield and food security. In this pilot study, swe employed computational metabolomics genome wide association survey (cmGWAS) of Gln metabolite profiling in Oryza sativa, targeting at the identification of abiotic stress responsible (AbSR) - Gln metabolite producing genes (GlnMPG). Identified 5 AbSR-GlnMPG alter the metabolite levels and play a predominant role in delineating the physiological significance of rice. These genes were systematically analysed for their biological features via OryzaCyc. Spatio-temporal and plant hormonal expression pattern of AbSR-GlnMPG was analysed and their differential expression profiling were noted in 48 different tissues and hormones, respectively. Furthermore, comparative ideogram of these genes revealed the chromosomal synteny with C4 grass genomes. Molecular crosstalks of these proteins, unravelled the various metabolic interaction. The systems expression profiling of AbSR-GlnMPG will lead to unravel the metabolite signaling and putative responses in multiple AbS. On the whole, this holistic study provides deeper insights on biomolecular features of AbSR-GlnMPG, which could be analysed further to decipher their functional metabolisms in AbS dynamism.

Genetic Basis of Maize Resistance to Multiple Insect Pests: Integrated Genome-Wide Comparative Mapping and Candidate Gene Prioritization

Several species of herbivores feed on maize in field and storage setups, making the development of multiple insect resistance a critical breeding target. In this study, an association mapping panel of 341 tropical maize lines was evaluated in three field environments for resistance to fall armyworm (FAW), whilst bulked grains were subjected to a maize weevil (MW) bioassay and genotyped with Diversity Array Technology's single nucleotide polymorphisms (SNPs) markers. A multi-locus genome-wide association study (GWAS) revealed 62 quantitative trait nucleotides (QTNs) associated with FAW and MW resistance traits on all 10 maize chromosomes, of which, 47 and 31 were discovered at stringent Bonferroni genome-wide significance levels of 0.05 and 0.01, respectively, and located within or close to multiple insect resistance genomic regions (MIRGRs) concerning FAW, SB, and MW. Sixteen QTNs influenced multiple traits, of which, six were associated with resistance to both FAW and MW, suggesting a pleiotropic genetic control. Functional prioritization of candidate genes (CGs) located within 10-30 kb of the QTNs revealed 64 putative GWAS-based CGs (GbCGs) showing evidence of involvement in plant defense mechanisms. Only one GbCG was associated with each of the five of the six combined resistance QTNs, thus reinforcing the pleiotropy hypothesis. In addition, through in silico co-functional network inferences, an additional 107 network-based CGs (NbCGs), biologically connected to the 64 GbCGs, and differentially expressed under biotic or abiotic stress, were revealed within MIRGRs. The provided multiple insect resistance physical map should contribute to the development of combined insect resistance in maize.

Global transcriptome analysis of novel stress associated protein (SAP) genes expression dynamism of combined abiotic stresses in Oryza sativa (L.)

Genes encoding proteins with A20/AN1 zinc-finger domains, belonging to the stress associated protein (SAP) gene family, are present in all eukaryotes and play a decisive role in plant response to diverse physiological and molecular activities particularly on biotic and abiotic stresses (AbS). In this first and foremost study, global transcriptome analysis of members of the SAP gene family was carried out in C3 model-Oryza sativa (OsSAP) aiming at the identification of OsSAP genes activated in response to unique or Combined AbS (CAbS). Based on the available spatio-temporal and phytohormonal RNA-Seq expression profile datasets, nine OsSAP genes were filtered out and identified by a differential expression signature noted in various tissues as well as plant hormones. Comparative genome ideogram of OsSAP genes confirmed the orthologous collinearity with C4 panicoid genomes. Interactome of these genes, revealed the molecular cross-talks of OsSAP. Thus, the computational expression signature of OsSAP genes led to a better understanding of gene dynamism in diverse developmental tissues/organs. Transcriptional regulation analysis of key OsSAP genes in response to stress (drought and salinity) suggested the novel role of OsSAP1, OsSAP2, OsSAP5, OsSAP7, OsSAP8 and OsSAP11 in AbS. Altogether, the study provides deeper insights on molecular characteristics of OsSAP genes, which could be deployed further to decipher their precise functional roles in AbS responses.Communicated by Ramaswamy H. Sarma.

Global integrated omics expression analyses of abiotic stress signaling HSF transcription factor genes in Oryza sativa L.: An in silico approach

Among the significant transcription factors (TFs), HSF proteins play pivotal roles in the regulation of hormonal signal transduction and different abiotic stress (AbS) responses. Hence considering its importance, global omics expression analysis of HSF candidates was performed in rice (OsHSF). The current study identified 25 HSF family members and physically plotted them against the rice genome. These proteins were systematically analyzed for their physicochemical features, organization and expression signatures. Further, heatmap of both spatio-temporal and global plant hormones revealed the developmental tissues and hormone specific expression profiling of these genes respectively. Comparative genome mapping between OsHSF players in interrelated C4 grass species revealed the chromosome level synteny. Signalome analysis revealed the protein - protein interactions of OsHSF. Expression profiling of key players in response to stresses exhibited the new involvement in combined AbS (CAbS) responses. Our results are significantly valuable to decipher their functional analysis of CAbS tolerant in rice.

Identification of key genes and its chromosome regions linked to drought responses in leaves across different crops through meta-analysis of RNA-Seq data

BACKGROUND

Our study is the first to provide RNA-Seq data analysis related to transcriptomic responses towards drought across different crops. The aim was to identify and map which genes play a key role in drought response on leaves across different crops. Forty-two RNA-seq samples were analyzed from 9 published studies in 7 plant species (Arabidopsis thaliana, Solanum lycopersicum, Zea mays, Vitis vinifera, Malus X domestica, Solanum tuberosum, Triticum aestivum).

RESULTS

Twenty-seven (16 up-regulated and 11 down-regulated) drought-regulated genes were commonly present in at least 7 of 9 studies, while 351 (147 up-regulated and 204 down-regulated) were commonly drought-regulated in 6 of 9 studies. Across all kind of leaves, the drought repressed gene-ontologies were related to the cell wall and membrane re-structuring such as wax biosynthesis, cell wall organization, fatty acid biosynthesis. On the other hand, drought-up-regulated biological processes were related to responses to osmotic stress, abscisic acid, water deprivation, abscisic-activated signalling pathway, salt stress, hydrogen peroxide treatment. A common metabolic feature linked to drought response in leaves is the repression of terpenoid pathways. There was an induction of AL1 (alfin-like), UGKYAH (trihelix), WRKY20, homeobox genes and members of the SET domain family in 6 of 9 studies. Several genes involved in detoxifying and antioxidant reactions, signalling pathways and cell protection were commonly modulated by drought across the 7 species. The chromosome (Chr) mapping of these key abiotic stress genes highlighted that Chr 4 in Arabidopsis thaliana, Chr 1 in Zea mays, Chr 2 and Chr 5 in Triticum aestivum contained a higher presence of drought-related genes compared to the other remaining chromosomes. In seedling studies, it is worth notice the up-regulation of ERF4 and ESE3 (ethylene), HVA22 (abscisic acid), TIR1 (auxin) and some transcription factors (MYB3, MYB94, MYB1, WRKY53 and WRKY20). In mature leaves, ERF1 and Alfin-like 1 were induced by drought while other transcription factors (YABBY5, ARR2, TRFL2) and genes involved phospholipid biosynthesis were repressed.

CONCLUSIONS

The identified and mapped genes might be potential targets of molecular breeding activities to develop cultivars with enhanced drought resistance and tolerance across different crops.

Abiotic stress induced miRNA-TF-gene regulatory network: A structural perspective

MicroRNAs (miRNAs) and transcription factors (TFs) are the largest families of trans-acting gene regulatory species, which are pivotal players in a complex regulatory network. Recently, extensive research on miRNAs and TFs in agriculture has identified these trans-acting regulatory species, as an effective tool for engineering new crop cultivars to increase yield and quality as well tolerance to environmental stresses but our knowledge of regulatory network is still not sufficient to decipher the exact mechanism. In the current work, stress-specific TF-miRNA-gene network was built for Arabidopsis under drought, cold, salt and waterlogging stress using data from reliable publically available databases; and transcriptome and degradome sequence data analysis by meta-analysis approach. Further network analysis elucidated significantly dense, scale-free, small world and hierarchical backbone of interactions. The various centrality measures highlighted several genes/TF/miRNAs as potential targets for tolerant variety cultivation. This comprehensive regulatory information will accelerate the advancement of current understanding on stress specific transcriptional and post-transcriptional regulatory mechanism and has promising utilizations for experimental biologist who are intended to improve plant crop performance under multiple Abiotic stress environments.

Comparative transcriptome meta-analysis of Arabidopsis thaliana under drought and cold stress

Multiple environmental stresses adversely affect plant growth and development. Plants under multiple stress condition trigger cascade of signals and show response unique to specific stress as well as shared responses, common to individual stresses. Here, we aim to identify common and unique genetic components during stress response mechanisms liable for cross-talk between stresses. Although drought and cold stress have been widely studied, insignificant information is available about how their combination affects plants. To that end, we performed meta-analysis and co-expression network comparison of drought and cold stress response in Arabidopsis thaliana by analyzing 390 microarray samples belonging to 29 microarray studies. We observed 6120 and 7079 DEGs (differentially expressed genes) under drought and cold stress respectively, using Rank Product methodology. Statistically, 28% (2890) DEGs were found to be common in both the stresses (i.e.; drought and cold stress) with most of them having similar expression pattern. Further, gene ontology-based enrichment analysis have identified shared biological processes and molecular mechanisms such as—‘photosynthesis’, ‘respiratory burst’, ‘response to hormone’, ‘signal transduction’, ‘metabolic process’, ‘response to water deprivation’, which were affected under cold and drought stress. Forty three transcription factor families were found to be expressed under both the stress conditions. Primarily, WRKY, NAC, MYB, AP2/ERF and bZIP transcription factor family genes were highly enriched in all genes sets and were found to regulate 56% of common genes expressed in drought and cold stress. Gene co-expression network analysis by WGCNA (weighted gene co-expression network analysis) revealed 21 and 16 highly inter-correlated gene modules with specific expression profiles under drought and cold stress respectively. Detection and analysis of gene modules shared between two stresses revealed the presence of four consensus gene modules.

Global analysis of threonine metabolism genes unravel key players in rice to improve the abiotic stress tolerance

The diversity in plant metabolites with improved phytonutrients is essential to achieve global food security and sustainable crop yield. Our study using computational metabolomics genome wide association study (cmGWAS) reports on a comprehensive profiling of threonine (Thr) metabolite in rice. Sixteen abiotic stress responsive (AbSR) – Thr metabolite producing genes (ThrMPG), modulate metabolite levels and play a significant role determining both physiological and nutritional importance of rice. These AbSR-ThrMPG were computationally analysed for their protein properties using OryzaCyc through plant metabolic network analyser. A total of 1373 and 1028 SNPs were involved in complex traits and genomic variations. Comparative mapping of AbSR-ThrMPG revealed the chromosomal colinearity with C4 grass species. Further, computational expression pattern of these genes predicted a differential expression profiling in diverse developmental tissues. Protein interaction of protein coding gene sequences revealed that the abiotic stresses (AbS) are multigenic in nature. In silico expression of AbSR-ThrMPG determined the putative involvement in response to individual AbS. This is the first comprehensive genome wide study reporting on AbSR –ThrMPG analysis in rice. The results of this study provide a pivotal resource for further functional investigation of these key genes in the vital areas of manipulating AbS signaling in rice improvement.

Molecular characterization of the 14-3-3 gene family in rice and its expression studies under abiotic stress

14-3-3 isoforms were relatively less conserved at the C-terminal region across plant groups. Both Os 14-3-3f and Os 14-3-3g were inducible with differential gene expression levels under different abiotic stress and developmental stages in sensitive and tolerant indica rice cultivars as confirmed both at transcript and protein level. Plant 14-3-3s has been well characterized to function in several signaling pathways, biotic as well as abiotic stress and nutrient metabolism. We attempted comprehensive analysis of 14-3-3 genes in different plant lineages such as green algae (Chlamydomonas reinhardtii), moss (Physcomitrella patens) and lycophyte (Selaginella moellendorffii), dicot Arabidopsis thaliana and monocot Oryza sativa sub sp. japonica at the gene and protein level. Sequence alignment results revealed that 14-3-3 isoforms were evolutionarily conserved across all taxa with variable C-terminal end. Phylogenetic analysis indicated that the majority of 14-3-3 isoforms in rice belong to the non-epsilon group that clustered separately from the dicot group. Segmental duplication event played a significant role in the expansion of both, Arabidopsis and rice, 14-3-3 isoforms as revealed by synteny studies. In silico gene expression using Massive Parallel Signature Sequencing and microarray analysis revealed that 14-3-3 isoforms have variable expression in different tissue types and under different abiotic stress regime in Arabidopsis and japonica rice. Both, semi-quantitative and qPCR results, confirmed that Os14-3-3f and Os14-3-3g were inducible under abiotic stress in lamina and roots of indica rice and relatively higher under salinity and cold stress in Nonabokra, under dehydration stress in N-22 and under exogenous ABA in IR-29 usually after 3-6 h of treatment. Both, 14-3-3f and 14-3-3g, were highly expressed in flag leaves, stems and panicles and mature roots. These results were further confirmed by immunoblot analysis of rice cultivars using Os14-3-3f antibody generated from recombinant Os14-3-3f protein. The results provide the first comprehensive report of Os14-3-3 gene expression in indica rice cultivars which differ in tolerance to abiotic stress that might be useful for further research.