Abstract
Water deficit stress at the early stage of development is one of the main factors limiting pearl millet production. One practice to counteract this limitation would be to resort to the application of hormones to stimulate plant growth and development at critical stages. Exogenous methyl jasmonate (MeJA) can improve drought tolerance by modulating signaling, metabolism, and photosynthesis pathways, therefore, we assumed that can occur in pearl millet during the early stage of development. To decipher the molecular mechanisms controlling these pathways, RNAseq was conducted in two pearl millet genotypes, drought-sensitive SosatC88 and drought-tolerant Souna3, in response to 200 μM of MeJA. Pairwise comparison between the MeJA-treated and non-treated plants revealed 3270 differentially expressed genes (DEGs) among 20,783 transcripts in SosatC88 and 127 DEGs out of 20,496 transcripts in Souna3. Gene ontology (GO) classification assigned most regulated DEGs in SosatC88 to heme binding, oxidation–reduction process, response to oxidative stress and membrane, and in Souna3 to terpene synthase activity, lyase activity, magnesium ion binding, and thylakoid. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis reveals that DEGs in SosatC88 are related to the oxidation–reduction process, the biosynthesis of other secondary metabolites, the signal transduction, and the metabolism of terpenoids, while in Souna3, DEGs are related to the metabolism of terpenoids and the energy metabolism. Two genes encoding a diterpenoid biosynthesis-related (Pgl_GLEAN_10009413) and a Glutathione S transferase T3 (Pgl_GLEAN_10034098) were contra-regulated between SosatC88 and Souna3. Additionally, five random genes differentially expressed by RNAseq were validated using qPCR, therefore, they are potential targets for the development of novel strategies breeding schemes for plant growth under water deficit stress. These insights into the molecular mechanisms of pearl millet genotype tolerance at the early stage of development contribute to the understanding of the role of hormones in adaptation to drought-prone environments.
Similar content being viewed by others
Introduction
Worldwide, water deficit contributes to crop production losses, particularly in the Sub-Saharan Africa countries, where agriculture depends on rainfall. In addition, population growth increases food demand while arable land is continuously shrinking. Meanwhile, climate forecasts are predicted to be extremes by 20501 in Africa, therefore, there is a need to breed high-yielding drought-tolerant genotypes. One of the major staple food and drought-tolerant plant is pearl millet (Pennisetum glaucum L. Br), ranked fourth among the cultivated tropical cereals. Since its domestication from West Africa 4,900 years ago2, pearl millet constitutes a staple food for more than 100 million people, with 90% grown by smallholder farmers3, hence contributing to nutrition and food security. Indeed, its grain is gluten-free, rich in vitamins, micronutrients (iron and zinc), and protein (8–19%)4. The genome size of the crop is estimated to 1.79 gigabytes containing 38,579 genes and recent studies led to the identification of specific marker-traits linked to key agronomic traits, such as grain yield, stem and leaf biomass as well as tremendous adaptive traits for heat and drought5. As a cross-pollinated C4 plant (2n = 2x = 14), pearl millet endows with a large genetic diversity and grows in some of the hottest and driest areas across the arid and semi-arid regions. Its vulnerability to drought is genotype dependent6.
Drought-related stresses at early stage of pearl millet plant development are limiting factors to its production7. Water deficit for example induces a limitation of photosynthesis activities due to reduced leaf expansion, impaired photosynthetic machinery, premature leaf senescence8,9, stomatal closure10, and the activity of Calvin’s cycle enzymes such as Rubisco and phosphoenolpyruvate carboxylase11. Reactive oxygen species (ROS), produced under water deficit, target various organelles, including chloroplasts, mitochondria, and peroxisomes, resulting in cell membrane instability, senescence, or plant death12. Other physiological responses to water deficit include decreasing stem height and diameter, leaf number, leaf area ratio, dry matter, shoot/root weight ratio, net CO2 assimilation and chlorophyll fluorescence, photosynthetic activity, altered cell wall elasticity and generation of toxic metabolites13,14. Another important factor leading to plant tolerance to water deficit is hormonal regulation. The expression of several transcription factors and their target genes involved in mediating phytohormones abscisic acid (ABA) under water deficit conditions, key components of perception and signaling, are induced to attenuate the negative impacts due to water deficit13,15.
Exogenous treatment of abscisic acid16 and salicylic acid17 improves tolerance under water deficit in wheat (Triticum aestivum). The application of jasmonates induced improved drought tolerance in barley (Hordeum vulgare),18, soybean (Glycine max)19, and sugar beet (Beta vulgaris)20. A treatment of Methyl jasmonate (MeJA) can enhance tolerance to dehydration of plants subjected to osmotic stress caused by polyethylene glycol21 and water retention22. MeJA governs many aspects of plant development, including seed germination, root growth, flowering, fruit ripening, senescence23,24 and abiotic stresses tolerance25. This hormone mitigates the influence of water deficit by promoting stomatal closure induced by H2O2 generation26, acting principally on the expression of genes involved in signaling processes.
Transcriptomic studies in pearl millet have been used to reveal responsive genes under drought, heat, or salinity conditions27,28. However, the mechanism of pearl millet in response to MeJA treatment has a negligible focus. Proteome has also been analyzed in pearl millet that provided insights into the functional enzyme involves in drought responses. Comprehensive physiological and proteomic responses of the drought stressed pearl millet revealed significant tissue-specific signatures during the adaptation to drought conditions29. In addition, a stay-green phenotype signature in drought-tolerant pearl millet was defined that contributes to a significant effect on yield, grain filling and chlorophyll loss30. The present study investigates genes that are differentially expressed after an exogenous MeJA treatment and with a mimic of water deficit by withholding water supply to decipher the molecular basis modulating the response of pearl millet under water scarcity conditions. By spraying pearl millet plants at the early stage of development with MeJA, the assumption that a differential gene expression could activate or repress genes underlying key metabolic and biosynthesis pathways that lead to coping with water deficit-related stress was tested. Comparative transcriptome analysis of two pearl millet genotypes contrasting by their level of tolerance to water deficit at the early stage identified sets of differentially expressed genes (DEGs) between non-treated and MeJA-treated plants. The regulation of DEGs modulated by MeJA and new genetic insights of adaptive mechanisms of pearl millet in response to water deficit stress are discussed.
Results
MeJA effects at early-stage development of pearl millet
Plants of two genotypes contrasting for the response to drought (sensitive SosatC88 and tolerant Souna3) showed no significant developmental changes for leaves number and plant height, p > 0.05 at 17 days after sowing, i.e. T0. The aerial part of these plants was then daily sprayed with MeJA (200 µM) for 10 days over T0 and we observed that the non-treated plants of both SosatC88 and Souna3 appeared withered compared to those treated (Fig. 1a,b). Plant height, leaves number and chlorophyll content were measured in both genotypes. MeJA did not significantly change leaves number in SosatC88 at T1, i.e. five days after treatment (p = 0.72), and T2, i.e. 10 days after treatment (p = 0.052). However, in Souna3 at T1 MeJA-treated plants have significantly more leaves (p = 0.049) than non-treated plants (Fig. 1c,d). In both genotypes, the height of MeJA-treated plants stayed lower than that of non-treated plants until the end of the treatment with a significant difference observed (p = 0.03) in SosatC88 at T2 (Fig. 1e,f) in treated and non-treated plants. In both genotypes, the chlorophyll content decreased during the 10 days of water deficit and remained lower in MeJA-treated plants, however, this change was only significant in Souna3 at T1 (p = 0.03) (Fig. 1g,h).
MeJA treatment in modulating gene expression
To evaluate whether the MeJA treatment in both SosatC88 and Souna3 mimicked drought tolerance response at the gene expression level, we conducted qPCR using RNA extracted from the aerial part of plants sprayed or not for 10 days with MeJA. The relative expression levels of pearl millet lipoxygenase 2 (PgLox2) and of 9-cis-epoxycarotenoid (PgNCED) were determined. The two genes were previously used as endogenously regulated after a MeJA treatment or in response to drought-induced stress, respectively. The results showed that PgLox2 was over-expressed in MeJA-treated plants compared to the non-treated (Supplemental figures, Fig. S1), indicating that MeJA penetrated pearl millet plants and did change gene expression. A significant positive fold change PgNCED expression profile was observed in both drought-sensitive SosatC88 and drought-tolerant Souna3, indicating that the MeJA treatment modules the response to drought-induced stress in pearl millet by modulating gene expression (Supplemental figures, Fig. S1).
RNA-seq assembly and analysis after MeJA treatment
To further explore this contrasting behavior between genotypes, we analyzed the transcriptomic responses to MeJA treatment on drought-sensitive SosatC88 and drought-tolerant Souna3. A total of 408,935.888 raw reads were generated, with an average of 34,077.991 reads per sample. After filtration and quality check, 379,364.208 high-quality reads were obtained with an average of 30,624.486 reads per sample. On average, 92.77% of total data passed a Phred score of ≥ Q30 bases for each sample. A final number of 367,493.828 high-quality reads (passing filter and correctly assigned to a sample after demultiplexing) were then aligned onto the pearl millet reference genome and 89.29–90.81% reads from twelve samples were uniquely mapped (Supplemental file, Table S1).
Differentially expressed genes (DEGs) for all samples were identified using edgeR package with |log2FoldChange|> 2 and p-adj < 0.05. In the drought-sensitive SosatC88, 3270 of 20,783 transcripts were differentially expressed (717 down-regulated and 2553 up-regulated). In drought-tolerant Souna3, 127 of 20,496 transcripts were differentially expressed (56 down-regulated and 71 up-regulated). This indicates that there were more differential genes in the drought-sensitive SosatC88 than in the drought-tolerant Souna3. Pairwise comparison (Fig. 2a-c) between non-treated SosatC88 and Souna3 plants revealed a total of 517 up-regulated genes and 149 down-regulated genes. Between MeJA-treated SosatC88 and Souna3, 265 genes are up-regulated and 473 genes are down-regulated. Comparison from non-treated SosatC88 and Souna3 plants and MeJA-treated SosatC88 and Souna3 revealed that 70 and 40 genes are up-regulated and down-regulated, respectively. However, 12 genes were contra-regulated (Fig. 2d, Supplemental figures). Comparison from SosatC88 and Souna3 DEGs showed 49 and 25 genes are up-regulated and down-regulated, respectively. However, two genes encoding a diterpenoid biosynthesis-related (Pgl_GLEAN_10009413) and a Glutathione S transferase T3 (Pgl_GLEAN_10034098), respectively, were contra-regulated between SosatC88 and Souna3 (Fig. 2e, Supplemental figures).
(a, b and c) Number of DEGs in all combinations with fold change > 2 or < -2 and FDR-corrected p-value < 0.05) Red and blue bars indicate up-regulated and down-regulated respectively; (d) Venn diagram of number of differentially expressed transcripts in ‟SosatC88 vs Souna3 (non-treated)” (yellow) and ‟SosatC88 vs Souna3 (MeJA-treated)” ( blue); (e) Venn diagram of number of differentially expressed transcripts in ‟SosatC88_MeJA- treated vs non-treated” ( yellow) and ‟Souna3_MeJA- treated vs non-treated” (blue).
DEGs functional analysis after MeJA treatment
Gene ontology (GO) assignments of biological process (BP), cellular component (CC) and molecular function (MF) were used to classify the differentially expressed genes and their expected functions from the datasets identified above. GO classification of pearl millet gene lists from DEGs were also carried out. Comparing non-treated SosatC88 and non-treated Souna3 , a total of 45 GO terms was assigned (13 CC, 10 MF and 22 BP) (Fig. 3a), 40 GO terms between MeJA-treated SosatC88 and MeJA-treated Souna3 (13 CC, 9 MF and 18 BP) (Fig. 3b), 48 GO terms between non-treated and MeJA-treated SosatC88(14 CC, 10 MF and 24 BP) (Fig. 4a), and 30 GO terms between non-treated and MeJA-treated Souna3 (10 CC, 4 MF and 16 BP) (Fig. 4b). This classification reveals that in non-treated conditions, differences between the two genotypes are particularly in cell (141 genes) cell part (140 genes), membrane (122 genes), catalytic activity (273 genes), binding (235 genes), metabolic process (252 genes) and cellular process (176 genes). In MeJA-treated conditions, the most differences appeared in cell (179 genes), cell part (177 genes), membrane (191 genes), membrane part (149 genes), catalytic activity (305 genes), binding (294 genes), cellular process (223 genes), and metabolic process (273 genes) indicating that some important cellular processes and metabolic activities occurred in the leaves of water-stressed millet in response to MeJA treatment.
According to the GO ontology, MeJA treatment modulated several genes involved in cell part, cell, organelle, protein-containing complex, organelle part, membrane, membrane part, extracellular region, catalytic activity, binding, transporter activity, localization, biological regulation, response to stimulus, metabolic process, cellular process, regulation of the biological process, signaling, cellular component organization or biogenesis. In drought-sensitive SosatC88, the DEGs are greater than in drought-tolerant Souna3 (Fig. 4a,b). However, MeJA treatment during non-watered conditions has a particular impact on the genes involved in the structure of the membrane and the binding genes.
GO terms enrichment analysis
Goseq analysis was used to classify the functions of DEGs. Between both non-treated SosatC88 and Souna3, with related terms ‛extracellular region’ (cellular component category), ‛pathogenesis’ (biological process category) and ‛transferase activity, transferring acyl group’ (molecular function category) were significantly enriched (Fig. 5a). Between MeJA-treated SosatC88 and Souna3, none of the categories are significantly enriched; but molecular function category including ‛magnesium ion binding’, ‛lysase activity’, ‛terpene synthase activity’, ‛DNA binding’, ‛nucleotide binding’, ‛polysaccharide binding’, ‛calcium binding’, ‛oxidoreductase activity, acting on the CH-OH group of donors, NAD or NADP as acceptor’ and cellular component including ‛nucleus’ and ‛nucleosome’ were the top over-represented categories (Fig. 5b).
Between SosatC88 MeJA-treated and non-treated, molecular function category including ‛heme binding’, ‛iron ion binding’, ‛terpene synthase activity’, ‛lyase activity’ ‛oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen’, ‛electron transfer activity’, ‛peroxidase activity’ and biological process category including ‛oxidation–reduction process’, ‛photosynthesis, light harvesting’, ‛response to oxidative stress’ were the most highly enriched (Fig. 6a). On the other hand, molecular function category, ‛terpene synthase activity’, ‛lyase activity’, ‛magnesium ion binding’; cellular component category, ‛thylakoid’ and biological process ‛photosynthetic electron transport chain’ were significantly enriched in Souna3_MeJA-treated vs non-treated. The GO molecular function category, ‛terpene synthetase activity’ and ‛lysase activity’ were also highly enriched in both SosatC88_MeJA-treated vs non-treated and Souna3_MeJA-treated vs non-treated (Fig. 6b).
Pathways analysis
Based on the KEGG database, several signaling pathways and genes involved in the response of both millet varieties to MeJA treatment in water-stressed conditions were discovered. Comparing non-treated SosatC88 and Souna3 and MeJA-treated SosatC88 and Souna3, an increase in carbon metabolism, pyruvate metabolism, oxidative phosphorylation was observed (Fig. 7a,b and Supplementary file). Propanoate metabolism, inositol phosphate metabolism, carbon fixation in photosynthetic organisms, biosynthesis of unsaturated fatty acids are induced by MeJA treatment. The most representative pathways of the overexpressed genes in between non-treated and MeJA-treated SosatC88 and non-treated and MeJA-treated Souna3 were those associated with metabolism (KO 01100), biosynthesis of secondary metabolites (KO 01110), antibiotic biosynthesis (KO 01,130), carbon metabolism (KO 01,200), photosynthesis (KO 00,195) and plant hormone signal transduction (KO 04,075) (Fig. 7c,d and Supplementary file). Some enriched pathways such as oxidative phosphorylation (KO 00,190), carbon fixation in photosynthetic organisms (KO 00,710), phenylalanine metabolism (KO 00,360), phenylalanine tyrosine tryptophan metabolism (KO 00,400), phenylpropanoid (KO 00,940), ubiquitin-mediated proteolysis (KO 04,120), MAPK signaling (KO 04,016) and circadian rhythm (KO 04,712) are only represented in SosatC88, indicating that these pathways were recruited by the MeJA hormone to enhance SosatC88 adaptively under water deficit.
KEGG pathway enrichment analysis of the differentially expressed genes (DEGs)72.
Validation of gene expression using qPCR
Five genes differentially expressed under or not treatment of MeJA of SosatC88 and Souna3 were selected based on log2 fold change and adjusted p-value parameters (FDR < 0.05) and their expression was validated by qPCR. These DEGs encode for: a small Auxin Up-Regulated RNA 50 (SAUR50), a peroxidase A2 (PA2), a myeloblastosis 30 (MYB30), a receptor kinase-like protein Xa21 and a glutathione S transferase T3 (GstT3). These DEGs included plant hormone signal transduction related genes, stress-responsive genes, peroxidase activity-related genes. Differential expression of transcripts was observed in SosatC88 and Souna3 between non-treated and MeJA-treated plants in concordance with the RNAseq data (Fig. 8a,b). The results also validated the contra-regulation revealed in RNAseq of PgGstT3 by MeJA treatment in SosatC88 and Souna3. A Pearson correlation was performed to test whether the qPCR results correlated with the RNAseq results. The ratio of expression levels obtained between non-treated and treated plants by qPCR was compared to the ratio of expression measured by RNA-Seq. A significant correlation (r2 = 0.777, p = 0.02) was observed between qPCR and RNA-Seq data that validate the differential expression of genes (Fig. 8c).
Validation of RNASeq data with qPCR. (a) and (b) Expression of five randomly selected genes was examined by qPCR analysis. For each gene, fold changes were calculated by ΔΔCt method and log2Fold change were compared between qPCR and RNAseq. (c) Correlation between RNAseq and qPCR data based on log2fold change of the five selected genes: y = 0.31x – 0.8.
Discussion
Water deficit is considered one of the most threatening abiotic stresses in agriculture in the era of climate change31. Water deficit resulting from these drought periods affects different aspects of plant growth and development, leading in some cases to a considerable crop production loss. These negative impacts can be alleviated by a wide range of plant hormonal treatments. MeJA, which is known to be involved in stress resistance mechanisms in plants, may have a key role to play in providing useful solutions. MeJA significantly reduced growth and leaf expansion in soybeans (Glycine max), hampered seedling diameter in scots pine (Pinus sylvestris)32,33,34, and inhibited both root and shoot growth in rice (Oryva sativa)35. The treatment of MeJA 200 µM at the early stage of pearl millet development did not induce significant visible morphological change between genotypes. However, the MeJA-treated plants retained more vigor than the non-treated plants. Similar results were observed in rice (Oryza sativa L.) which the drought-sensitive genotype raised with primed seeds with 100 µM of MeJA retained its morphology better than the tolerant genotype under drought condition36. In wheat (Triticum aestivum L.), 0.25 µM of MeJA delayed senescence under drought condition37. In pearl millet, key characteristics of drought resistance are related to root morphology as well as an efficient photosynthetic machinery30. Our results revealed a reduction in chlorophyll content that suggests a loss of photosynthetic pigments or activity. These are in accordance with previous studies33,38,39 that observed a decline in photosynthesis of soybean (Glycine max), barley (Hordeum vulgare) leaves and Arabidopsis thaliana. The loss of photosynthetic pigments would decrease the amount of energy absorbed by the photosynthetic system, thus reducing the energy required for metabolism, which consequently impacts plant growth and development. The loss of chlorophyll content could also be attributed to leaf senescence40. Furthermore, it has been shown that jasmonates attenuate the influence of drought by favoring the closure of stomata26,41. Our findings suggest that treatment of MeJA could cause the stomatal closure, leading to the inhibition of CO2 absorption and the reduction of photosynthetic activity11,42,43.
Despite non-significant morphological changes between genotypes, MeJA treatment was effective in modulating endogenous stress-related genes expression. It is known that gene encoding Lox2 is highly induced by jasmonates (JA and MeJA) or in response to biotic stress44 and NCED is the key enzyme in the ABA biosynthetic pathway in plants45,46. The overexpression of pearl millet Lox2 (PgLox2) is induced by the jasmonates produced by the plants in response and in response to MeJA treatment. The expression profile of pearl millet NCED (PgNCED) in MeJA-treated and non-treated plants suggests that water deficiency and MeJA treatment induced an increase in ABA biosynthesis, which is most important with external MeJA input. These results suggest that the decrease in plant growth is caused by increased ABA and MeJA accumulation.
To gain insights in the molecular basis of pearl millet responses to MeJA treatment, a transcriptome analysis was carried out at the early stage of plant development. Transcriptomic data were generated to identify the molecular basis at the genome-wide level underlying the response of pearl millet to MeJA treatment at the early developmental stage. Under non-treated conditions, distribution of DEGs in drought-sensitive SosatC88 and drought-tolerant Souna3 in extracellular region term, pathogenesis term and transferase activity, transferring acyl group term according to GO analysis suggested that water-deficiency involves genes related to these terms. In contrast, when MeJA was applied there was no significant difference in the transcriptomic response between SosatC88 and Souna3. This suggests that MeJA stimulates in SosatC88 pathways and biological processes to better adapt to water deficiency as in Souna3. In addition, 12 genes were found to be differentially expressed in the opposite manner depending on whether we are in MeJA-treated or non-treated conditions. Among these genes, two are up-regulated by the MeJA treatment. These two genes are Pgl_GLEAN_10004431 encoding for an unknown protein and Pgl_GLEAN_10018402 encoding for an uncharacterized protein family UPF0136 transmembrane. The additional ten genes were down-regulated by the MeJA treatment. They are mostly Tyrosine-protein kinase (Pgl_GLEAN_10005256), Multicopper oxidase, type 2 (Pgl_GLEAN_10001778), Peptidase (Pgl_GLEAN_10016112), Phospholipid/glycerol acyltransferase (Pgl_GLEAN_10012582), Glycoside hydrolase (Pgl_GLEAN_10004142) and uncharacterized genes. About the uncharacterized genes, they could likely represent key regulators of the variation in flowering time or drought tolerance among different genotypes. Further investigation is required for understanding the role in pearl millet adaptation to stress.
Data reported here indicate in one hand that the MeJA triggers higher transcriptional responses in SosatC88 to better adapt and cope with water deficiency in opposite to Souna3. The striking difference in DEGs between SosatC88 and Souna3 provides information on the amplitude of the effect of the MeJA treatment on the transcriptome of SosatC88. Consistent with data reported here, other transcriptome studies report a higher number of DEGs under water stress-sensitive variety in wheat (Triticum aestivum)47 and sorghum (Sorghum bicolor)48. MeJA induced in both genotype enrichment of GO terms related to terpene synthetase activity and lyase activity. One way to protect against oxidative stress is the biosynthesis of volatile terpenoids, which is thought to quench ROS49. Lyases are known to play a role in antioxidant activity in cases of biotic or abiotic stress50. We theorize that these groups of genes probably constitute the characteristic basis for the adaptation to water deficit conferred by the treatment of MeJA in both genotypes SosatC88 and Souna3. On the other hand, the results of the GO enrichment analysis of up-regulated gene groups provide evidence of striking antioxidant enzyme activities. Several important functional gene groups related to water deficiency stress including oxidation–reduction process, heme binding, electron transfer activity, response to oxidative stress, membrane, peroxidase activity, among others, are significantly enriched in SosatC88 treated with MeJA. Peroxidases are involved in the defense against abiotic stress through their role in the scavenging of ROS51. They are also involved in cell wall formation, confer cell wall rigidity and reducing leaf expansion32. It has also been shown in several studies that MeJA treatment increases peroxidase activity32,37,39. A cross-species meta-analysis of progressive drought stress identifies the heme-binding family protein as key drought adaptive genes conserved52. Our results show that MeJA improves SosatC88 adaptation, drought-sensitive genotype, to water deficiency by increasing and diversifying its antioxidant response mechanisms to cope under these unfavorable conditions.
On the other hand, proteins such as peroxidase superfamily protein, thioredoxin-dependent peroxidase, prohibitin 3 and a subunit of the NAD(P)H complex located in the chloroplast thylakoid lumen (Pgl_GLEAN_10037666) identified with high level in pearl millet leaves under drought stress 29 have their genes expression up-regulated under drought stress conditions. These groups of proteins were also up-regulated in the leaves under drought in bean (Phaseolus vulgaris L.53, in wheat (Triticum aestivum L.)54 and in maize (Zea mays)55. Furthermore, ribulose bisphosphate carboxylase small chain (Pgl_GLEAN_10020566), a protein identified as associated with a degradation of chlorophyll30. Together, they are associated with the processes that underlie the senescence or the stay-green phenotype and overexpressed in SosatC88 by MeJA. This shows the importance of these proteins as they have passed post-transcriptional regulation to play a major role in the response to drought.
Conclusion
A successful treatment of exogenous MeJA in two genotypes with different response to drought did not alter plant morphology, but modulated gene expression. RNAseq analysis revealed a transcriptional process that mobilizes a large number of DEGs related to cell wall protection and detoxification, implying different mechanisms involved in response to drought at the early stage of the pearl millet developmental cycle. Five genes differentially expressed are genuine candidates as molecular markers in functional genomics or in the breeding of pearl millet. Overall, data give an indication that, in pearl millet, exogenous MeJA treatment may play a dual role by simulating abiotic stress defenses to protect plants under water deficit and contribute to crop production in drought-prone environments.
Materials and methods
Plant material and growth Conditions
Seeds of two pearl millet open-pollinated varieties, SosatC88 and Souna3, were used. SosatC88 was developed by recombining 19 S1 progeny selected at Cinzana (Mali) in 1988 from a cross between the local Souna and Sanio genotypes. It was considered a water deficit-sensitive. Souna3 is a synthetic variety resulting from the recombination of eight lines from the PC 28 and PC 32 populations (106–7, 108–4, 113–3, 115–4, 134–5, 142–4, 143–4, 148–3). Souna3 is a drought tolerant genotype which reaches maturity between 85–95 days and can be grown in areas with moderate rainfall ranging between 400 and 750 mm. All plant material used in this study are provided by the ISRA genebank and their uses comply with relevant institutional, national, and international guidelines and legislation.
The experiment was performed in the shadehouse of the Centre d'Etude Régional pour l'Amélioration de l'Adaptation à la Sécheresse (CERAAS) in Thiès (Sénégal; 14°45′57.17″N, 16°53′31″W) from March to April 2018 with a photoperiod of 13 h of light/day with a light intensity of 1158 µmol lux, a vapor pressure deficit (Vpd) of 2.3 Kpascal, an average air temperature of 35.3 ºC day and 23.5 ºC night.
Both genotypes were sown in 10 L plastic buckets containing a non-autoclaved sandy soil with low levels of clay and silt (12%) and organic matter (0.4%) from Bambey (Senegal)7. In each bucket, 4 plants of the same variety were grown. These buckets were completely randomized with 6 replicates per treatment.
The plants were watered three times a week by pouring an equivalent volume of water on each soil bucket until 17 days after sowing (DAS). Then, water was withdrawn and plant leaves were sprayed with a solution containing 200 μM MeJA (Methyl jasmonate 95%, CAS Number: 39924–52-2 (Methyl 3-oxo-2-(2-pentenyl) cyclopentaneacetate from Sigma-Aldrich) dissolved in 0.1% tween 20 and with 0.1% tween 20 for the MeJA-treated and the non-treated, respectively, every day for 10 days24. Three days after the last treatment, the aerial parts (shoot and leaf) of the MeJA-treated and non-treated plants were collected at late seedling growth stage (i.e. phase 17 according to the BBCH scale56, https://en.wikipedia.org/wiki/Cereal_growth_staging_scales), immediately frozen in liquid nitrogen and stored at – 80 °C for further RNA extraction.
Development and Chlorophyll density parameters measurement
Leaf number and plant height were assessed at T0, T1 and T2. T0 corresponds to the first day just before treatment (both non-treated and MeJA-treated plants), T1 and T2 correspond to 5 and 10 days after treatment, respectively. The chlorophyll density level was measured at T1 and T2 using a SPAD 502 PLUS CHLOROPHYLL-METER57. Four plants from the same pot were used for measurement of leaves number, plant height and chlorophyll density level, which resulted in six biological replicates for each genotype in each sampling date. A significant difference in mean between the samples as obtained by Student’s t-test at P < 0.05.
RNA extraction
Leaves from four plants of each genotype were pooled in one sample. Sampling was performed on three biological replicates (MeJA-treated and non-treated plants)58. The samples were ground in liquid nitrogen and then, 200 mg of powder was used for total RNA isolation with TRI Reagent (Sigma-Aldrich) according to the instructions of the manufacturer. RNA was then DNase treated (Qiagen, Germany) and purified using Qiagen RNeasy Minikit. The quality of the RNA was evaluated on 1.2% agarose gel (Supplemental figures, Fig. S2) and electrophoresis was carried out at 100 Volts for 30 min and the quantification was performed using a Nano Drop Lite and the Agilent 2100 bioanalyzer.
Gene markers expression level assay for MeJA treatment validation
To validate the results observed at the end of the MeJA treatment, we identified genes that are involved in response to plant biotic and abiotic stresses24. Among these, we used Lipoxygenase 2 (Lox2) involved in the biosynthesis of MeJA and 9-cis-epoxycarotenoid (NCED) involved in abscisic acid (ABA) biosynthesis. Both are associated in plant response to abiotic/biotic stresses. Lox2 encodes the enzyme that catalyzes an important step in jasmonic acid (JA) biosynthesis from linolenic acid derived from the membrane59,60 while NCED is the key enzyme in ABA biosynthesis involved in plant response to water deficit46.
Primer pairs for PgLox2 and PgNCED were designed using Primer3Plus online tool61 according to the target gene-specific sequence design, the option of qPCR setting Primer3Plus online tool platform was used. The primer sequences of transcripts are listed in Supplemental file, Table S2.
qPCR assays
Five hundred (500) ng of the total RNA was reverse transcribed using the Goscript Reverse Transcriptase system (Promega) according to the manufacturer’s instructions. The qPCR reaction was performed in a total volume of 10 µl containing 5 µl of GoTaq® qPCR Master Mix (2X), 0.25 μM of each primer, 1 μl (100 ng) of diluted cDNA and free nuclease water and then, loaded in a StepOnePlus™ Real-Time PCR Systems. The program used for qPCR contains three steps: the first step is the initial activation at 95 °C for 2 min; the next step is the PCR amplification performed up to 40 cycles with 95 °C for 15 s, 58 °C for 60 s and the third step is a melt curve analysis ramping from 58 °C to 95 °C.
The relative expression levels of transcripts under the different conditions were normalized based on the expression of an endogenous control gene, ubiquitin 5 (UBQ5)62 using the 2−∆∆Ct method63. The qbase + software Version: 3.264 was used for analysis.
RNAseq
Library construction and sequencing
For library construction, RNAs were purified on RNAClean XP beads and then analyzed by capillary electrophoresis (Fragment Analyser). Highly pure mRNA was isolated from the total RNA using oligo (dT) beads. The Illumina TruSeq RNA Library Prep Kit v2 was used to synthesize the second strand cDNAs library. Sequencing was performed on an Illumina HiSeq 2500 using the Sequence By Synthesis (SBS) technique using the TruSeq Rapid SBS kit (Illumina). Sequencing results were obtained as single-end reads (50 bp each) in the FASTQ format.
Differential gene expression analysis
Clean reads were aligned to the reference genome5 and counted with Bioconductor R package Rsubrea65. The read counts were used to performed the Differential Expression Genes analysis with edgeR package of Bioconductor66. Pairwise comparisons were performed between four datasets: non-treated SosatC88 vs non-treated Souna3, MeJA-treated SosatC88 vs MeJA-treated Souna3, MeJA-treated SosatC88 vs non-treated SosatC88, and MeJA-treated Souna3 vs non-treated Souna3. False Discovery Rate (FDR) < 0.05 and log|Fold change|> 2 were considered as the conditions to state the genes as differentially expressed.
Functional annotation of transcripts and genes enrichment analysis
The sequence similarity search was carried out against the National Center for Biotechnology Information (NCBI) non-redundant (nr) protein database and the Swiss-Prot protein database using the BLASTx algorithm from Blast2GO67 specifying E-values < 10–5. Gene Ontology (GO) categorization was done with Blast2GO. From the annotated GO file, we classified and plotted the DEGs according to the official classification (Molecular Function, Biological Process and Cellular Component), using WEGO 2.068,69 and functional enrichment analysis was also performed using Goseq70, a bioconductor R package implemented in galaxy70 using Wallenius method.
DEGs were then mapped on the biological pathways of the web-based Kyoto Encyclopedia of Genes and Genomes (KEGG) web-based annotation server (KASS) by running Blastx against KEGG GENES (Kyoto Encyclopedia) https://www.genome.jp/kegg/kaas/. We completed the KEGG annotation using the KEGG result files of millet genome information uploaded to GIGAdb71. The result contains KO assignments (KEGG Orthology) and KEGG pathways were automatically generated. KEGG enrichment analysis was also performed based on the same criteria described above for GO.
Data availability
All data provided in the article and Supplementary files are available from the corresponding author on reasonable request.
References
IPCC, 2014: Summary for policymakers. In: Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Field, C.B., V.R. Barros, D.J. Dokken, K.J. Mach, M.D. Mastrandrea, T.E. Bilir, M. Chatterjee, K.L. Ebi, Y.O. Estrada, R.C. Genova, B. Girma, E.S. Kissel, A.N. Levy, S. MacCracken, P.R. Mastrandrea, and L.L.White (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, pp. 1–32.
Burgarella, C. et al. A western Sahara centre of domestication inferred from pearl millet genomes. Nat. Ecol. Evol. 2, 1377–1380 (2018).
Taylor JRN. Millet Pearl: Overview. in Encyclopedia of Food Grains - 2nd Edition (ed. Press, O. A.) 190–198 (2016).
Nambiar, V. S., Dhaduk, J. J., Sareen, N., Shahu, T. & Desai, R. Potential functional implications of pearl millet (Pennisetum glaucum) in health and disease. J. Appl. Pharm. Sci. 1, 62–67 (2011).
Varshney, R. K. et al. Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments. Nat. Biotechnol. 35, 969–976 (2017).
Rhoné, B. et al. Pearl millet genomic vulnerability to climate change in West Africa highlights the need for regional collaboration. Nat. Commun. 11 (2020).
Debieu, M. et al. Response to early drought stress and identification of QTLs controlling biomass production under drought in pearl millet. PLoS ONE 13, 1–19 (2018).
Diouf, O. et al. Response of Pearl Millet to nitrogen as affected by water deficit. Agronomie 24, 77–84 (2004).
Kusaka, M., Lalusin, A. G. & Fujimura, T. The maintenance of growth and turgor in pearl millet (Pennisetum glaucum [L] Leeke) cultivars with different root structures and osmo-regulation under drought stress. Plant Sci. 168, 1–14 (2005).
Zheng, Y. X., Wu, J. C., Cao, F. L. & Zhang, Y. P. Effects of water stress on photosynthetic activity, dry mass partitioning and some associated metabolic changes in four provenances of neem (Azadirachta indica A. Juss). Photosynthetica 48, 361–369 (2010).
Bota, J., Medrano, H. & Flexas, J. Is photosynthesis limited by decreased Rubisco activity and RuBP content under progressive water stress?. New Phytol. 162, 671–681 (2004).
Ma, C., Wang, Z., Kong, B. & Lin, T. Exogenous trehalose differentially modulate antioxidant defense system in wheat callus during water deficit and subsequent recovery. Plant Growth Regul. 70, 275–285 (2013).
Ahuja, I., de Vos, R. C. H., Bones, A. M. & Hall, R. D. Plant molecular stress responses face climate change. Trends Plant Sci. 15, 664–674 (2010).
Ramírez, V. et al. Drought tolerance in Arabidopsis is controlled by the OCP3 disease resistance regulator. Plant J. 58, 578–591 (2009).
Sirichandra, C., Wasilewska, A., Vlad, F., Valon, C. & Leung, J. The guard cell as a single-cell model towards understanding drought tolerance and abscisic acid action. J. Exp. Bot. 60, 1439–1463 (2009).
Du, Y. L. et al. Exogenous abscisic acid reduces water loss and improves antioxidant defence, desiccation tolerance and transpiration efficiency in two spring wheat cultivars subjected to a soil water deficit. Funct. Plant Biol. 40, 494–506 (2013).
Kang, G. et al. Proteomics reveals the effects of salicylic acid on growth and tolerance to subsequent drought stress in wheat. J. Proteome Res. 11, 6066–6079 (2012).
Bandurska, H., Stroiński, A. & Kubiś, J. The effect of jasmonic acid on the accumulation of ABA, proline and spermidine and its influence on membrane injury under water deficit in two barley genotypes. Acta Physiol. Plant. 25, 279–285 (2003).
Mohamed, H. I. & Latif, H. H. Improvement of drought tolerance of soybean plants by using methyl jasmonate. Physiol. Mol. Biol. Plants 23, 545–556 (2017).
Fugate, K. K. et al. Methyl jasmonate alleviates drought stress in young sugar beet (Beta vulgaris L.) plants. J. Agron. Crop Sci. 204, 566-576. (2018).
Li, L., Van Staden, J. & Jäger, A. K. Effects of plant growth regulators on the antioxidant system in seedlings of two maize cultivars subjected to water stress. Plant Growth Regul. 25, 81–87 (1998).
Wang, S. Y. Methyl jasmonate reduces water stress in strawberry. J. Plant Growth Regul. 18, 127–134 (1999).
Wasternack, C. Jasmonates: An update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 100, 681–697 (2007).
Diallo, A. O. et al. Transcriptome analysis of an mvp mutant reveals important changes in global gene expression and a role for methyl jasmonate in vernalization and flowering in wheat. J. Exp. Bot. 65, 2271–2286 (2014).
Ahmad, P. et al. Jasmonates: Multifunctional roles in stress tolerance. Front. Plant Sci. 7, 813 (2016).
Suhita, D., Raghavendra, A. S., Kwak, J. M. & Vavasseur, A. Cytoplasmic alkalization precedes reactive oxygen species production during methyl jasmonate- and abscisic acid-induced stomatal closure. Plant Physiol. 134, 1536–1545 (2004).
Jaiswal, S. et al. Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources. Sci. Rep. 8, 3382 (2018).
Mukesh Sankar, S. et al. Differential modulation of heat-inducible genes across diverse genotypes and molecular cloning of a sHSP from Pearl Millet [Pennisetum glaucum (L.) R. Br.]. Front. Plant Sci. 12, 1–21 (2021).
Ghatak, A. et al. Comprehensive tissue-specific proteome analysis of drought stress responses in Pennisetum glaucum (L.) R. Br. (Pearl millet). J. Proteomics 143, 122–135 (2016).
Ghatak, A. et al. Physiological and Proteomic Signatures Reveal Mechanisms of Superior Drought Resilience in Pearl Millet Compared to Wheat. Front. Plant Sci. 11, 600278 (2021).
Sultan, B. & Gaetani, M. Agriculture in West Africa in the twenty-first century: Climate change and impacts scenarios, and potential for adaptation. Front. Plant Sci. 7, 1262 (2016).
Moore, J. P., Paul, N. D., Whittaker, J. B. & Taylor, J. E. Exogenous jasmonic acid mimics herbivore-induced systemic increase in cell wall bound peroxidase activity and reduction in leaf expansion. Funct. Ecol. 17, 549–554 (2003).
Anjum, S. A., Wang, L., Farooq, M., Khan, I. & Xue, L. Methyl Jasmonate-induced alteration in lipid peroxidation, antioxidative defence system and yield in soybean under drought. J. Agron. Crop Sci. 197, 296–301 (2011).
Heijari, J. et al. Application of methyl jasmonate reduces growth but increases chemical defence and resistance against Hylobius abietis in Scots pine seedlings. in Entomologia Experimentalis et Applicata vol. 115 117–124 (Wiley, 2005).
Tsai, F.-Y., Lin, C. C. & Kao, C. H. A comparative study of the effects of abscisic acid and methyl jasmonate on seedling growth of rice. Plant Growth Regul. 21, 37–42 (1997).
Sasi, M. et al. Plant growth regulator induced mitigation of oxidative burst helps in the management of drought stress in rice (Oryza sativa L.). Environ. Exp. Bot. 185, 104413 (2021).
Ma, C., Wang, Z. Q., Zhang, L. T., Sun, M. M. & Lin, T. B. Photosynthetic responses of wheat (Triticum aestivum L.) to combined effects of drought and exogenous methyl jasmonate. Photosynthetica 52, 377–385 (2014).
Hristova, V. A. & Popova, L. P. Treatment with methyl jasmonate alleviates the effects of paraquat on photosynthesis in barley plants. Photosynthetica 40, 567–574 (2002).
Jung, S. Effect of chlorophyll reduction in Arabidopsis thaliana by methyl jasmonate or norflurazon on antioxidant systems. Plant Physiol. Biochem. 42, 225–231 (2004).
He, Y., Fukushige, H., Hildebrand, D. F. & Gan, S. Evidence supporting a role of jasmonic acid in Arabidopsis leaf senescence. Plant Physiol. 128, 876–884 (2002).
Gehring, C. A., Irving, H. R., McConchie, R. & Parish, R. W. Jasmonates induce intracellular alkalinization and closure of Paphiopedilum guard cells. Ann. Bot. 80, 485–489 (1997).
Nayyar, H. & Gupta, D. Differential sensitivity of C3 and C4 plants to water deficit stress: Association with oxidative stress and antioxidants. Environ. Exp. Bot. 58, 106–113 (2006).
Chaves, M. M., Flexas, J. & Pinheiro, C. Photosynthesis under drought and salt stress: Regulation mechanisms from whole plant to cell. Ann. Bot. 103, 551–560 (2009).
Edreva, A. Pathogenesis-related proteins: research progress in the last 15 years. Gen Appl Plant Physiol 31, 105–124 (2005).
Zhang, M., Leng, P., Zhang, G. & Li, X. Cloning and functional analysis of 9-cis-epoxycarotenoid dioxygenase (NCED) genes encoding a key enzyme during abscisic acid biosynthesis from peach and grape fruits. J. Plant Physiol. 166, 1241–1252 (2009).
Iuchi, S. et al. Regulation of drought tolerance by gene manipulation of 9-cis- epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis. Plant J. 27, 325–333 (2001).
Mia, M. S., Liu, H., Wang, X., Zhang, C. & Yan, G. Root transcriptome profiling of contrasting wheat genotypes provides an insight to their adaptive strategies to water deficit. Sci. Rep. 10, 1–11 (2020).
Fracasso, A., Trindade, L. M. & Amaducci, S. Drought stress tolerance strategies revealed by RNA-Seq in two sorghum genotypes with contrasting WUE. BMC Plant Biol. 16, 1–18 (2016).
Possell, M. & Loreto, F. The role of volatile organic compounds in plant resistance to abiotic stresses: Responses and mechanisms. in 209–235 (Springer, Dordrecht, 2013). https://doi.org/10.1007/978-94-007-6606-8_8.
Gholizadeh, A. Anti-oxidation profile in the leaves of Maize Inbreds: Elevation in the activity of Phenylalanine Ammonia Lyase under drought Stress. J. Plant Sci. 5, 137–145 (2010).
Yazawa, T., Kawahigashi, H., Matsumoto, T. & Mizuno, H. Simultaneous transcriptome analysis of sorghum and bipolaris sorghicola by using RNA-seq in combination with De Novo Transcriptome Assembly. PLoS One 8, (2013).
Shaar-Moshe, L., Hübner, S. & Peleg, Z. Identification of conserved drought-adaptive genes using a cross-species meta-analysis approach. BMC Plant Biol. 15, 1–18 (2015).
Zadražnik, T., Hollung, K., Egge-Jacobsen, W., Meglič, V. & Šuštar-Vozlič, J. Differential proteomic analysis of drought stress response in leaves of common bean (Phaseolus vulgaris L.). J. Proteomics 78, 254–272 (2013).
Komatsu, S., Kamal, A. H. M. & Hossain, Z. Wheat proteomics: Proteome modulation and abiotic stress acclimation. Front. Plant Sci. 5, 1–19 (2014).
Riccardi, F., Gazeau, P., Jacquemot, M. P., Vincent, D. & Zivy, M. Deciphering genetic variations of proteome responses to water deficit in maize leaves. Plant Physiol. Biochem. 42, 1003–1011 (2004).
Hess, M. et al. Use of the extended BBCH scale - General for the descriptions of the growth stages of mono- and dicotyledonous weed species. Weed Res. 37, 433–441 (1997).
Dudhate, A., Shinde, H., Tsugama, D., Liu, S. & Takano, T. Transcriptomic analysis reveals the differentially expressed genes and pathways involved in drought tolerance in pearl millet [pennisetum glaucum (l.) r. Br.]. PLoS ONE 13, 1–14 (2018).
Auer, P. L. & Doerge, R. W. Statistical design and analysis of RNA sequencing data. Genetics 185, 405–416 (2010).
Beale, M. H. & Ward, J. L. Jasmonates: Key players in the plant defence. Nat. Prod. Rep. 15, 533–548 (1998).
Farmer, E. E., Alméras, E. & Krishnamurthy, V. Jasmonates and related oxylipins in plant responses to pathogenesis and herbivory. Curr. Opin. Plant Biol. 6, 372–378 (2003).
Untergasser, A. et al. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 35, W71 (2007).
Saha, P. & Blumwald, E. Assessing reference genes for accurate transcript normalization using quantitative real-time PCR in Pearl Millet [Pennisetum glaucum (L.) R. Br.]. PLoS ONE 9, e106308 (2014).
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408 (2001).
Hellemans, J., Mortier, G., De Paepe, A., Speleman, F. & Vandesompele, J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 8, R19 (2007).
Liao, Y., Smyth, G. K. & Shi, W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 47, e47 (2019).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2010).
Götz, S. et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435 (2008).
Ye, J. et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 34, 293–297 (2006).
Ye, J. et al. WEGO 2.0: A web tool for analyzing and plotting GO annotations, 2018 update. Nucleic Acids Res. 46, W71–W75 (2018).
Young, M. D., Wakefield, M. J., Smyth, G. K. & Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, R14 (2010).
GigaDB Dataset - DOI https://doi.org/10.5524/100192 - Genomic data of Pearl Millet (Pennisetum glaucum). http://gigadb.org/dataset/100192.
Kanehisa, M. & Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30 (2000).
Acknowledgements
We are grateful to Professor Fathey SARHAN (University of Quebec at Montreal (Qc, Canada) for providing the Methyl jasmonate 95% product. This research was supported by the United States Agency for International Development under Cooperative Agreement No. AOO-0868-S008-A04, through the FEED THE FUTURE Sorghum and Millet Innovation Lab agreement between Kansas State University and ISRA (GenMil Project S19182).
Author information
Authors and Affiliations
Contributions
A.O.D., D.D., and N.A.K. conceived and designed the study. A.N., N.C.F. and R.D.D. performed the experiments. A.N. performed RNA extraction, qPCR and transcriptomic analysis. A.N., A.O.D., D.D., and N.A.K. wrote the manuscript. N.A.K. performed the final editing of the manuscript. All authors reviewed and approved this submission.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ndiaye, A., Diallo, A.O., Fall, N.C. et al. Transcriptomic analysis of methyl jasmonate treatment reveals gene networks involved in drought tolerance in pearl millet. Sci Rep 12, 5158 (2022). https://doi.org/10.1038/s41598-022-09152-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-09152-6