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.
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, H2O2 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.
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.
The Arabidopsis thaliana accession C24 is a vernalization-responsive, moderately late flowering ecotype. We report that a mutation in AtMBD8, which encodes a protein with a putative Methyl-CpG-Binding Domain (MBD), in C24 background, results in a delay in flowering time during both long and short days. The atmbd8-1 mutant responded to vernalization as wild type (wt) plants. Consistent with a role in modulation of flowering time, an AtMBD8::GUS-reporter construct was expressed in the shoot meristem region and developing leaves. Full-genome transcriptional profiling revealed very few changes in gene expression between atmbd8-1 and wt plants. The expression level of FLC, the major repressor of transition to flowering, was unchanged in atmbd8-1, and in accordance with that, genes upstream of FLC were unaffected by the mutation. The expression level of CONSTANS, involved in photoperiodic control of flowering, was very similar in atmbd8-1 and wt plants. In contrast, the major promoters of flowering, FT, and SOC1, were both downregulated. As FT is a regulator of SOC1, we conclude that AtMBD8 is a novel promotor of flowering that acts upstream of FT in the C24 accession. In contrast to atmbd8-1, the Colombia (Col) SALK T-DNA insertion line, atmbd8-2, did not display a delayed transition to flowering. Transcriptional profiling revealed that a substantial number of genes were differentially expressed between C24 and Col wt seedlings. Several of these genes are also differentially expressed in late flowering mutants. We suggest that these differences contribute to the contrasting effect of a mutation in AtMBD8 in the two ecotypes.