Unraveling the regulatory role of MYC2 on ASMT gene expression in wheat: Implications for melatonin biosynthesis and drought tolerance.

Recognized for its multifaceted functions, melatonin is a hormone found in both animals and plants. In the plant kingdom, it plays diverse roles, regulating growth, development, and stress responses. Notably, melatonin demonstrates its significance by mitigating the effects of abiotic stresses like drought. However, understanding the precise regulatory mechanisms controlling melatonin biosynthesis genes, especially during monocots' response to stresses, requires further exploration. Seeking to understand the molecular basis of drought stress tolerance in wheat, we analyzed RNA-Seq libraries of wheat exposed to drought stress using bioinformatics methods. In light of our findings, we identified that the Myelocytomatosis oncogenes 2 (MYC2) transcription factor is a hub gene upstream of a main melatonin biosynthesis gene, N-acetylserotonin methyltransferase (ASMT), in the wheat drought response-gene network. Promoter analysis of the ASMT gene suggested that it might be a target gene of MYC2. We conducted a set of molecular and physiochemical assays along with robust machine learning approaches to elevate those findings further. MYC2 and ASMT were co-regulated under Jasmonate, drought, and a combination of them in the leaf tissues of wheat was detected. A meaningful correlation was observed among gene expression profiles, melatonin contents, photosynthetic activities, antioxidant enzyme activities, H2 O2 levels, and plasma membrane damage. The results indicated an evident relationship between jasmonic acid and the melatonin biosynthesis pathway. Moreover, it seems that the MYC2-ASMT module might contribute to wheat drought tolerance by regulating melatonin contents.

Potential role of the regulatory miR1119-MYC2 module in wheat (Triticum aestivum L.) drought tolerance.

MicroRNA (miRNA)-target gene modules are essential components of plants' abiotic stress signalling pathways Little is known about the drought-responsive miRNA-target modules in wheat, but systems biology approaches have enabled the prediction of these regulatory modules and systematic study of their roles in responses to abiotic stresses. Using such an approach, we sought miRNA-target module(s) that may be differentially expressed under drought and non-stressed conditions by mining Expressed Sequence Tag (EST) libraries of wheat roots and identified a strong candidate (miR1119-MYC2). We then assessed molecular and physiochemical differences between two wheat genotypes with contrasting drought tolerance in a controlled drought experiment and assessed possible relationships between their tolerance and evaluated traits. We found that the miR1119-MYC2 module significantly responds to drought stress in wheat roots. It is differentially expressed between the contrasting wheat genotypes and under drought versus non-stressed conditions. We also found significant associations between the module's expression profiles and ABA hormone content, water relations, photosynthetic activities, H2O2 levels, plasma membrane damage, and antioxidant enzyme activities in wheat. Collectively, our results suggest that a regulatory module consisting of miR1119 and MYC2 may play an important role in wheat's drought tolerance.

MicroRNA miR1118 contributes to wheat (Triticum aestivum L.) salinity tolerance by regulating the Plasma Membrane Intrinsic Proteins1;5 (PIP1;5) gene.

microRNAs (miRNAs) are important regulators of various adaptive stress responses in crops; however, many details about associations among miRNAs, their target genes and physiochemical responses of crops under salinity stress remain poorly understood. We designed this study in a systems biology context and used a collection of computational, experimental and statistical procedures to uncover some regulatory functions of miRNAs in the response of the important crop, wheat, to salinity stress. Accordingly, under salinity conditions, wheat roots' Expressed Sequence Tag (EST) libraries were computationally mined to identify the most reliable differentially expressed miRNA and its related target gene(s). Then, molecular and physiochemical evaluations were carried out in a separate salinity experiment using two contrasting wheat genotypes. Finally, the association between changes in measured characteristics and wheat salinity tolerance was determined. From the results, miR1118 was assigned as a reliable salinity-responsive miRNA in wheat roots. The expression profiles of miR1118 and its predicted target gene, Plasma Membrane Intrinsic Proteins1,5 (PIP1;5), significantly differed between wheat genotypes. Moreover, results revealed that expression profiles of miR1118 and PIP1;5 significantly correlate to Relative Water Content (RWC), root hydraulic conductance (Lp), photosynthetic activities, plasma membrane damages, osmolyte accumulation and ion homeostasis of wheat. Our results suggest a plausible regulatory node through miR1118 adjusting the wheat water status, maintaining ion homeostasis and mitigating membrane damages, mainly through the PIP1;5 gene, under salinity conditions. To our knowledge, this is the first report on the role of miR1118 and PIP1;5 in wheat salinity response.

Evidence that miR168a contributes to salinity tolerance of Brassica rapa L. via mediating melatonin biosynthesis.

Melatonin is a master regulator of diverse biological processes, including plant's abiotic stress responses and tolerance. Despite the extensive information on the role of melatonin in response to abiotic stress, how plants regulate endogenous melatonin content under stressful conditions remains largely unknown. In this study, we computationally mined Expressed Sequence Tag (EST) libraries of salinity-exposed Chinese cabbage (Brassica rapa) to identify the most reliable differentially expressed miRNA and its target gene(s). In light of these analyses, we found that miR168a potentially targets a key melatonin biosynthesis gene, namely O-METHYLTRANSFERASE 1 (OMT1). Accordingly, molecular and physiochemical evaluations were performed in a separate salinity experiment using contrasting B. rapa genotypes. Then, the association between B. rapa salinity tolerance and changes in measured molecular and physiochemical characteristics was determined. Results indicated that the expression profiles of miR168a and OMT1 significantly differed between B. rapa genotypes. Moreover, the expression profiles of miR168a and OMT1 significantly correlated with more melatonin content, robust antioxidant activities, and better ion homeostasis during salinity stress. Our results suggest that miR168a plausibly mediates melatonin biosynthesis, mainly through the OMT1 gene, under salinity conditions and thereby contributes to the salinity tolerance of B. rapa. To our knowledge, this is the first report on the role of miR168a and OMT1 in B. rapa salinity response.

A systems biology study unveils the association between a melatonin biosynthesis gene, O-methyl transferase 1 (OMT1) and wheat (Triticum aestivum L.) combined drought and salinity stress tolerance.

MAIN CONCLUSION: Enhanced levels of endogenous melatonin in the root of wheat, mainly through the OMT1 gene, augment the antioxidant system, reestablish redox homeostasis and are associated with combined stress tolerance. A systems biology approach, including a collection of computational analyses and experimental assays, led us to uncover some aspects of a poorly understood phenomenon, namely wheat (Triticum aestivum L.) combined drought and salinity stress tolerance. Accordingly, a cross-study comparison of stress experiments was performed via a meta-analysis of Expressed Sequence Tags (ESTs) data from wheat roots to uncover the overlapping gene network of drought and salinity stresses. Identified differentially expressed genes were functionally annotated by gene ontology enrichment analysis and gene network analysis. Among those genes, O-methyl transferase 1 (OMT1) was highlighted as a more important (hub) gene in the dual-stress response gene network. Afterwards, the potential roles of OMT1 in mediating physiochemical indicators of stress tolerance were investigated in two wheat genotypes differing in abiotic stress tolerance. Regression analysis and correspondence analysis (CA) confirmed that the expression profiles of the OMT1 gene and variations in melatonin content, antioxidant enzyme activities, proline accumulation, H2O2 and malondialdehyde (MDA) contents are significantly associated with combined stress tolerance. These results reveal that the OMT1 gene may contribute to wheat combined drought and salinity stress tolerance through augmenting the antioxidant system and re-establishing redox homeostasis, probably via the regulation of melatonin biosynthesis as a master regulator molecule. Our findings provide new insights into the roles of melatonin in wheat combined drought and salinity stress tolerance and suggest a novel plausible regulatory node through the OMT1 gene to improve multiple-stress tolerant crops.

Rubisco activase A (RcaA) is a central node in overlapping gene network of drought and salinity in Barley (Hordeum vulgare L.) and may contribute to combined stress tolerance.

Co-occurrence of abiotic stresses, especially drought and salinity, is a natural phenomenon in field conditions and is worse for crop production than any single stress. Nowadays, rigorous methods of meta-analysis and systems biology have made it possible to perform cross-study comparisons of single stress experiments, which can uncover main overlapping mechanisms underlying tolerance to combined stress. In this study, a meta-analysis of RNA-Seq data was conducted to obtain the overlapping gene network of drought and salinity stresses in barley (Hordeum vulgare L.), which identified Rubisco activase A (RcaA) as a hub gene in the dual-stress response. Thereafter, a greenhouse experiment was carried out using two barley genotypes with different abiotic stress tolerance and evaluated several physiochemical properties as well as the expression profile and protein activity of RcaA. Finally, machine learning analysis was applied to uncover relationships among combined stress tolerance and evaluated properties. We identified 441 genes which were differentially expressed under both drought and salinity stress. Results revealed that the photosynthesis pathway and, in particular, the RcaA gene are major components of the dual-stress responsive transcriptome. Comparative physiochemical and molecular evaluations further confirmed that enhanced photosynthesis capability, mainly through regulation of RcaA expression and activity as well as accumulation of proline content, have a significant association with combined drought and salinity stress tolerance in barley. Overall, our results clarify the importance of RcaA in combined stress tolerance and may provide new insights for future investigations.

Plausible association between drought stress tolerance of barley (Hordeum vulgare L.) and programmed cell death via MC1 and TSN1 genes.

Studying the drought-responsive transcriptome is of high interest as it can serve as a blueprint for stress adaptation strategies. Despite extensive studies in this area, there are still many details to be uncovered, such as the importance of each gene involved in the stress response as well as the relationship between these genes and the physiochemical processes governing stress tolerance. This study was designed to address such important details and to gain insights into molecular responses of barley (Hordeum vulgare L.) to drought stress. To that, we combined RNA-seq data analysis with field and greenhouse drought experiments in a systems biology approach. RNA-sequence analysis identified a total of 665 differentially expressed genes (DEGs) belonging to diverse functional categories. A gene network was derived from the DEGs, which comprised of a total of 131 nodes and 257 edges. Gene network topology analysis highlighted two programmed cell death (PCD) modulating genes, MC1 (metacaspase 1) and TSN1 (Tudor-SN 1), as important (hub) genes in the predicted network. Based on the field trial, a drought-tolerant and a drought-susceptible barley genotype was identified from eight tested cultivars. Identified genotypes exhibited different physiochemical characteristics, including proline content, chlorophyll concentration, percentage of electrolyte leakage and malondialdehyde content as well as expression profiles of MC1 and TSN1 genes. Machine learning and correspondence analysis revealed a significant relationship between drought tolerance and measured characteristics in the context of PCD. Our study provides new insights which bridge barley drought tolerance to PCD through MC1 and TSN1 pathway.

Molecular characterization of Brassica napus stress related transcription factors, BnMYB44 and BnVIP1, selected based on comparative analysis of Arabidopsis thaliana and Eutrema salsugineum transcriptomes.

Many studies have been performed to identify regulatory circuit underlying plant stress tolerance. However, the reliability of some findings has been criticized because of exclusive use of stress sensitive plant species such as Arabidopsis thaliana. Sensitive plant species often harbor narrow defensive mechanisms and have relatively low capacity for adaptive responses. Therefore, it is useful to employ tolerant model plants, such as Eutrema salsugineum, to provide comprehensive insights into various mechanisms involved in response to abiotic stresses. In this study, comparative transcriptome and regulatory network analysis of stress-sensitive (A. thaliana) and -tolerant (E. salsugineum) model plants uncovered regulatory hierarchies underlying response to abiotic stresses and suggested the transcription factor genes, MYB44 and VIP1 as the candidate hub genes to perform molecular analyses on their Brassica napus homologs, BnMYB44 and BnVIP1. The full-length coding sequence of BnMYB44 and BnVIP1 with 891 and 969 bp long were cloned and sequenced. They shared high similarity with their counterparts in other plants at nucleotide and amino acid levels. The expression patterns of BnMYB44 and BnVIP1 genes of the two B. napus cultivars under drought and salt stress conditions coupled with the data obtained from the physiological measurements as well as analysis of the BnMYB44 and BnVIP1 promoters suggested that BnMYB44 and BnVIP1 genes may contribute to responses to drought and salt stresses in B. napus.

In silico identification of miRNAs and their target genes and analysis of gene co-expression network in saffron (Crocus sativus L.) stigma.

As an aromatic and colorful plant of substantive taste, saffron (Crocus sativus L.) owes such properties of matter to growing class of the secondary metabolites derived from the carotenoids, apocarotenoids. Regarding the critical role of microRNAs in secondary metabolic synthesis and the limited number of identified miRNAs in C. sativus, on the other hand, one may see the point how the characterization of miRNAs along with the corresponding target genes in C. sativus might expand our perspectives on the roles of miRNAs in carotenoid/apocarotenoid biosynthetic pathway. A computational analysis was used to identify miRNAs and their targets using EST (Expressed Sequence Tag) library from mature saffron stigmas. Then, a gene co- expression network was constructed to identify genes which are potentially involved in carotenoid/apocarotenoid biosynthetic pathways. EST analysis led to the identification of two putative miRNAs (miR414 and miR837-5p) along with the corresponding stem- looped precursors. To our knowledge, this is the first report on miR414 and miR837-5p in C. sativus. Co-expression network analysis indicated that miR414 and miR837-5p may play roles in C. sativus metabolic pathways and led to identification of candidate genes including six transcription factors and one protein kinase probably involved in carotenoid/apocarotenoid biosynthetic pathway. Presence of transcription factors, miRNAs and protein kinase in the network indicated multiple layers of regulation in saffron stigma. The candidate genes from this study may help unraveling regulatory networks underlying the carotenoid/apocarotenoid biosynthesis in saffron and designing metabolic engineering for enhanced secondary metabolites.

Mining expressed sequence tags of rapeseed (Brassica napus L.) to predict the drought responsive regulatory network.

It is of great significance to understand the regulatory mechanisms by which plants deal with drought stress. Two EST libraries derived from rapeseed (Brassica napus) leaves in non-stressed and drought stress conditions were analyzed in order to obtain the transcriptomic landscape of drought-exposed B. napus plants, and also to identify and characterize significant drought responsive regulatory genes and microRNAs. The functional ontology analysis revealed a substantial shift in the B. napus transcriptome to govern cellular drought responsiveness via different stress-activated mechanisms. The activity of transcription factor and protein kinase modules generally increased in response to drought stress. The 26 regulatory genes consisting of 17 transcription factor genes, eight protein kinase genes and one protein phosphatase gene were identified showing significant alterations in their expressions in response to drought stress. We also found the six microRNAs which were differentially expressed during drought stress supporting the involvement of a post-transcriptional level of regulation for B. napus drought response. The drought responsive regulatory network shed light on the significance of some regulatory components involved in biosynthesis and signaling of various plant hormones (abscisic acid, auxin and brassinosteroids), ubiquitin proteasome system, and signaling through Reactive Oxygen Species (ROS). Our findings suggested a complex and multi-level regulatory system modulating response to drought stress in B. napus.

A novel pairwise comparison method for in silico discovery of statistically significant cis-regulatory elements in eukaryotic promoter regions: application to Arabidopsis.

Cis regulatory elements (CREs), located within promoter regions, play a significant role in the blueprint for transcriptional regulation of genes. There is a growing interest to study the combinatorial nature of CREs including presence or absence of CREs, the number of occurrences of each CRE, as well as of their order and location relative to their target genes. Comparative promoter analysis has been shown to be a reliable strategy to test the significance of each component of promoter architecture. However, it remains unclear what level of difference in the number of occurrences of each CRE is of statistical significance in order to explain different expression patterns of two genes. In this study, we present a novel statistical approach for pairwise comparison of promoters of Arabidopsis genes in the context of number of occurrences of each CRE within the promoters. First, using the sample of 1000 Arabidopsis promoters, the results of the goodness of fit test and non-parametric analysis revealed that the number of occurrences of CREs in a promoter sequence is Poisson distributed. As a promoter sequence contained functional and non-functional CREs, we addressed the issue of the statistical distribution of functional CREs by analyzing the ChIP-seq datasets. The results showed that the number of occurrences of functional CREs over the genomic regions was determined as being Poisson distributed. In accordance with the obtained distribution of CREs occurrences, we suggested the Audic and Claverie (AC) test to compare two promoters based on the number of occurrences for the CREs. Superiority of the AC test over Chi-square (2×2) and Fisher's exact tests was also shown, as the AC test was able to detect a higher number of significant CREs. The two case studies on the Arabidopsis genes were performed in order to biologically verify the pairwise test for promoter comparison. Consequently, a number of CREs with significantly different occurrences was identified between the promoters. The results of the pairwise comparative analysis together with the expression data for the studied genes revealed the biological significance of the identified CREs.