Comparative analysis of alfalfa (Medicago sativa L.) seedling transcriptomes reveals genotype-specific drought tolerance mechanisms

Genome-Wide Association and Genomic Prediction of Alfalfa (Medicago sativa L.) Biomass Yield Under Drought Stress

Developing drought-resistant alfalfa (Medicago sativa L.) that maintains high biomass yield is a key breeding goal to enhance productivity in water-limited areas. In this study, 424 alfalfa breeding families were analyzed to identify molecular markers associated with biomass yield under drought stress and to predict high-merit plants. Biomass yield was measured from 18 harvests from 2020 to 2023 in a field trial with deficit irrigation. A total of 131 significant markers were associated with biomass yield, with 80 markers specifically linked to yield under drought stress; among these, 19 markers were associated with multiple harvests. Finally, genomic best linear unbiased prediction (GBLUP) was employed to obtain predictive accuracies (PAs) and genomic estimated breeding values (GEBVs). Removing low-informative SNPs [SNPs with p-values > 0.05 from the additive Genome-Wide Association (GWAS) model] for GBLUP increased PA by 47.3%. The high number of markers associated with yield under drought stress and the highest PA (0.9) represent a significant achievement in improving yield under drought stress in alfalfa.

Comparative transcriptomic and proteomic analyses of two salt-tolerant alfalfa (Medicago sativa L.) genotypes: investigation of the mechanisms underlying tolerance to salt

Abiotic stressors such as salt stress restrict plant development and output, which lowers agricultural profitability. In this study, alfalfa (Medicago sativa L.) varieties with different levels of salt tolerance were examined using high-throughput RNA sequencing (RNA-Seq) and Tandem Mass Tags (TMT) technologies to study the reactions of the root systems to salt stress, from transcriptomics and proteomics perspectives. The varieties Atlantic (AT) and Zhongmu-1 (ZM-1) were selected and evaluated after 2 h and 6 h of treatment with 150 mM NaCl. The results showed that under salt stress for 2 h, 1810 differentially expressed genes (DEGs) and 160 differentially expressed proteins (DEPs) in AT were screened, while 9341 DEGs and 193 DEPs were screened in ZM-1. Under salt stress for 6 h, 7536 DEGs and 118 DEPs were screened in AT, while 11,754 DEGs and 190 DEPs were screened in ZM-1. Functional annotation and pathway enrichment analyses indicated that the DEGS and DEPs were mainly involved in the glutathione metabolism, biosynthesis of secondary metabolites, glycolysis/gluconeogenesis, carbon fixation in photosynthetic organisms, and photosynthesis pathways. A series of genes related to salt tolerance were also identified, including GSTL3 and GSTU3 of the GST gene family, PER5 and PER10, of the PER gene family, and proteins such as APR and COMT, which are involved in biosynthesis of secondary metabolites. This study provides insights into salt resistance mechanisms in plants, and the related genes and metabolic pathways identified may be helpful for alfalfa breeding in the future.

The seed germination and seedling phytotoxicity of decabromodiphenyl ethane to tall fescue under citric acid amendment

The novel brominated flame retardant decabromodiphenyl ethane (DBDPE) has biological toxicity, persistence, long-range migration and bioaccumulation ability. However, there is currently little research on the phytotoxicity of DBDPE in plants. The perennial herbaceous plant tall fescue (Festuca elata Keng ex E. B. Alexeev) was selected as the model organism for use in seed germination experiments, and the phytotoxicity of DBDPE in the soil of tall fescue was studied. The results indicated that DBDPE had a significant effect on the germination and growth of tall fescue seedlings. Citric acid reduced the stress caused by DBDPE in plants, effectively alleviating the phytotoxicity of DBDPE in tall fescue. The root vitality and protein content significantly increased after the application of citric acid, increasing by 74.93-183.90%, 146.44-147.67%, respectively. The contents of proline and soluble sugars significantly decreased after the application of citric acid, decreasing by 45.18-59.69% and 23.03%, respectively (P < 0.05). There was no significant difference in superoxide dismutase (SOD) or peroxidase (POD) activity in tall fescue seedlings, and the catalase (CAT) activity and malondialdehyde (MDA) content were significantly lower after the application of citric acid, decreasing by 64.62-67.91% and 29.10-49.80%, respectively (P < 0.05). Tall fescue seedlings bioaccumulated DBDPE, with biological concentration factors (BCFs) ranging from 4.28 to 18.38 and transfer factors (TFs) ranging from 0.43 to 0.54. This study provides theoretical support for the study of the toxicity of DBDPE to plants and offers a research foundation for exploring the phytoremediation of DBDPE-contaminated soil by tall fescue.

Transcriptomic and Metabolomic Profiling of Root Tissue in Drought-Tolerant and Drought-Susceptible Wheat Genotypes in Response to Water Stress

Wheat is the most widely grown crop in the world; its production is severely disrupted by increasing water deficit. Plant roots play a crucial role in the uptake of water and perception and transduction of water deficit signals. In the past decade, the mechanisms of drought tolerance have been frequently reported; however, the transcriptome and metabolome regulatory network of root responses to water stress has not been fully understood in wheat. In this study, the global transcriptomic and metabolomics profiles were employed to investigate the mechanisms of roots responding to water stresses using the drought-tolerant (DT) and drought-susceptible (DS) wheat genotypes. The results showed that compared with the control group, wheat roots exposed to polyethylene glycol (PEG) had 25941 differentially expressed genes (DEGs) and more upregulated genes were found in DT (8610) than DS (7141). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the DEGs of the drought-tolerant genotype were preferably enriched in the flavonoid biosynthetic process, anthocyanin biosynthesis and suberin biosynthesis. The integrated analysis of the transcriptome and metabolome showed that in DT, the KEGG pathways, including flavonoid biosynthesis and arginine and proline metabolism, were shared by differentially accumulated metabolites (DAMs) and DEGs at 6 h after treatment (HAT) and pathways including alanine, aspartate, glutamate metabolism and carbon metabolism were shared at 48 HAT, while in DS, the KEGG pathways shared by DAMs and DEGs only included arginine and proline metabolism at 6 HAT and the biosynthesis of amino acids at 48 HAT. Our results suggest that the drought-tolerant genotype may relieve the drought stress by producing more ROS scavengers, osmoprotectants, energy and larger roots. Interestingly, hormone signaling plays an important role in promoting the development of larger roots and a higher capability to absorb and transport water in drought-tolerant genotypes.

Multi-omics integration analysis reveals the molecular mechanisms of drought adaptation in homologous tetraploid alfalfa(Medicago sativa 'Xinjiang-Daye')

Drought stress is a predominant abiotic factor leading to decreased alfalfa yield. Genomic ploidy differences contribute to varying adaptation mechanisms of different alfalfa cultivars to drought conditions. This study employed a multi-omics approach to characterize the molecular basis of drought tolerance in a tetraploid variant of alfalfa (Medicago sativa, Xinjiang-Daye). Under drought treatment, a total of 4446 genes, 859 proteins, and 524 metabolites showed significant differences in abundance. Integrative analysis of the multi-omics data revealed that regulatory modules involved in flavonoid biosynthesis, plant hormone signalling transduction, linoleic acid metabolism, and amino acid biosynthesis play crucial roles in alfalfa adaptation to drought stress. The severity of drought led to the substantial accumulation of flavonoids, plant hormones, free fatty acids, amino acids, and their derivatives in the leaves. Genes such as PAL, 4CL, CHI, CHS, PP2C, ARF_3, and AHP_4 play pivotal regulatory roles in flavonoid biosynthesis and hormone signalling pathways. Differential expression of the LOX gene emerged as a key factor in the elevated levels of free fatty acids. Upregulation of P5CS_1 and GOT1/2 contributed significantly to the accumulation of Pro and Phe contents. ERF19 emerged as a principal positive regulator governing the synthesis of the aforementioned compounds. Furthermore, observations suggest that Xinjiang-Daye alfalfa may exhibit widespread post-transcriptional regulatory mechanisms in adapting to drought stress. The study findings unveil the critical mechanisms by which Xinjiang-Daye alfalfa adapts to drought stress, offering novel insights for the improvement of alfalfa germplasm resources.

Tandem mass tag (TMT)-based quantitative proteomics analysis reveals the different responses of contrasting alfalfa varieties to drought stress

BACKGROUND

Drought stress restricts the growth, distribution and productivity of alfalfa (Medicago sativa L.). In order to study the response differences of alfalfa cultivars to drought stress, we previously carried out physiological and molecular comparative analysis on two alfalfa varieties with contrasting drought resistance (relatively drought-tolerant Longdong and drought-sensitive Algonquin). However, the differences in proteomic factors of the two varieties in response to drought stress still need to be further studied. Therefore, TMT-based quantitative proteomic analysis was performed using leaf tissues of the two alfalfa cultivars to identify and uncover differentially abundant proteins (DAPs).

RESULTS

In total, 677 DAPs were identified in Algonquin and 277 in Longdong under drought stress. Subsequently, we conducted various bioinformatics analysis on these DAPs, including subcellular location, functional classification and biological pathway enrichment. The first two main COG functional categories of DAPs in both alfalfa varieties after drought stress were 'Translation, ribosomal structure and biogenesis' and 'Posttranslational modification, protein turnover, chaperones'. According to KEGG database, the DAPs of the two alfalfa cultivars after drought treatment were differentially enriched in different biological pathways. The DAPs from Algonquin were enriched in 'photosynthesis' and 'ribosome'. The pathways of 'linoleic acid metabolism', 'protein processing in endoplasmic reticulum' and 'RNA transport' in Longdong were significantly enriched. Finally, we found significant differences in DAP enrichment and expression patterns between Longdong and Algonquin in glycolysis/glycogenesis, TCA cycle, photosynthesis, protein biosynthesis, flavonoid and isoflavonoid biosynthesis, and plant-pathogen interaction pathway after drought treatment.

CONCLUSIONS

The differences of DAPs involved in various metabolic pathways may explain the differences in the resistance of the two varieties to drought stress. These DAPs can be used as candidate proteins for molecular breeding of alfalfa to cultivate new germplasm with more drought tolerance to adapt to unfavorable environments.

MsCYP71 is a positive regulator for drought resistance in alfalfa

Cytochrome P450 (CYP450) family proteins play key roles in plant growth, development, stress responses, and other physiological processes. Here, we cloned the cytochrome P450 gene MsCYP71 in alfalfa and found that the expression of MsCYP71 was induced by drought stress. Silencing the MsCYP71 gene using virus-induced gene silencing technology significantly decreased the drought resistance of alfalfa, as indicated by their lower relative water content, net photosynthetic rate, and chlorophyll fluorescence maximum (Fm); further, the heterologous overexpression of MsCYP71 in tobacco significantly enhanced the drought resistance and Fm of transgenic tobacco. Furthermore, the expression of MsCYP71 across 45 alfalfa accessions under drought stress was investigated. A significant positive correlation between drought resistance and MsCYP71 expression was observed. The 45 alfalfa accessions were clustered into four groups, and drought resistance, Fm, and MsCYP71 were higher in group I than in the other groups, indicating that group I accessions can be used as candidate germplasm resources for the breeding of drought-resistant alfalfa varieties. Overall, our findings indicated that MsCYP71 is a positive regulator of drought resistance in alfalfa, and its expression can be used to evaluate the drought resistance of alfalfa.

Physiological response and transcriptome analyses of leguminous Indigofera bungeana Walp. to drought stress

Objective

Indigofera bungeana is a shrub with high quality protein that has been widely utilized for forage grass in the semi-arid regions of China. This study aimed to enrich the currently available knowledge and clarify the detailed drought stress regulatory mechanisms in I. bungeana, and provide a theoretical foundation for the cultivation and resistance breeding of forage crops.

Methods

This study evaluates the response mechanism to drought stress by exploiting multiple parameters and transcriptomic analyses of a 1-year-old seedlings of I. bungeana in a pot experiment.

Results

Drought stress significantly caused physiological changes in I. bungeana. The antioxidant enzyme activities and osmoregulation substance content of I. bungeana showed an increase under drought. Moreover, 3,978 and 6,923 differentially expressed genes were approved by transcriptome in leaves and roots. The transcription factors, hormone signal transduction, carbohydrate metabolism of regulatory network were observed to have increased. In both tissues, genes related to plant hormone signaling transduction pathway might play a more pivotal role in drought tolerance. Transcription factors families like basic helix-loop-helix (bHLH), vian myeloblastosis viral oncogene homolog (MYB), basic leucine zipper (bZIP) and the metabolic pathway related-genes like serine/threonine-phosphatase 2C (PP2C), SNF1-related protein kinase 2 (SnRK2), indole-3-acetic acid (IAA), auxin (AUX28), small auxin up-regulated rna (SAUR), sucrose synthase (SUS), sucrosecarriers (SUC) were highlighted for future research about drought stress resistance in Indigofera bungeana.

Conclusion

Our study posited I. bungeana mainly participate in various physiological and metabolic activities to response severe drought stress, by regulating the expression of the related genes in hormone signal transduction. These findings, which may be valuable for drought resistance breeding, and to clarify the drought stress regulatory mechanisms of I. bungeana and other plants.

Integrative analysis of transcriptome and metabolome revealed the mechanisms by which flavonoids and phytohormones regulated the adaptation of alfalfa roots to NaCl stress

INTRODUCTION

Salinity critically affects the growth and development of alfalfa (Medicago sativa), making it necessary to understand the molecular mechanism of alfalfa's adaptation to salt stress.

METHODS

In this study, alfalfa roots were subjected to salt stress and transcriptomics and metabolomics analyses were performed.

RESULTS

The results showed that flavonoid synthesis, hormone synthesis, and transduction pathways may be involved in the alfalfa salt stress adaptation reaction, and that they are related. Combined analysis of differential genes and differential metabolites found that dihydroquercetin and beta-ring hydroxylase (LUT5), ABA responsive element binding factor 2 (ABF2), protein phosphatase PP2C (PP2C) and abscisic acid (ABA) receptor PYL2 (PYL), luteolinidin was significantly correlated with PP2C and phytochrome-interacting factor 4 (PIF4) and (+)-7-isomethyl jasmonate were significantly correlated with flavonol synthase (FLS) gene. (+)-7-isomethyl jasmonate and homoeriodictyol chalcone were significantly correlated with peroxidase (POD). POD was significantly up-regulated under NaCl stress for 6 and 24 h. Moreover, flavonoids, gibberellin (GA), jasmonic acid (JA) and ABA were suggested to play an important role in alfalfa's response to salt stress. Further, GA,ABA, and JA may be involved in the regulation of flavonoids to improve alfalfa's salt tolerance, and JA may be a key signal to promote the synthesis of flavonoids.

DISCUSSION

This study revealed the possible molecular mechanism of alfalfa adaptation to salt stress, and identified a number of salt-tolerance candidate genes from the synthesis and signal transduction pathways of flavonoids and plant hormones, providing new insights into the regulatory network of alfalfa response to salt stress.

Regulation of Na+/H+ exchangers, Na+/K+ transporters, and lignin biosynthesis genes, along with lignin accumulation, sodium extrusion, and antioxidant defense, confers salt tolerance in alfalfa

Accumulation of high sodium (Na⁺) leads to disruption of metabolic processes and decline in plant growth and productivity. Therefore, this study was undertaken to clarify how Na⁺/H⁺ exchangers and Na⁺/K⁺ transporter genes contribute to Na⁺ homeostasis and the substantial involvement of lignin biosynthesis genes in salt tolerance in alfalfa (Medicago sativa L.), which is poorly understood. In this study, high Na⁺ exhibited a substantial reduction of morphophysiological indices and induced oxidative stress indicators in Xingjiang Daye (XJD; sensitive genotype), while Zhongmu (ZM; tolerant genotype) remained unaffected. The higher accumulation of Na⁺ and the lower accumulation of K⁺ and K⁺/(Na⁺ + K⁺) ratio were found in roots and shoots of XJD compared with ZM under salt stress. The ZM genotype showed a high expression of SOS1 (salt overly sensitive 1), NHX1 (sodium/hydrogen exchanger 1), and HKT1 (high-affinity potassium transporter 1), which were involved in K⁺ accumulation and excess Na⁺ extrusion from the cells compared with XJD. The lignin accumulation was higher in the salt-adapted ZM genotype than the sensitive XJD genotype. Consequently, several lignin biosynthesis-related genes including 4CL2, CCoAOMT, COMT, CCR, C4H, PAL1, and PRX1 exhibited higher mRNA expression in salt-tolerant ZM compared with XJD. Moreover, antioxidant enzyme (catalase, superoxide dismutase, ascorbate peroxidase, and glutathione reductase) activity was higher in ZM relative to XJD. This result suggests that high antioxidant provided the defense against oxidative damages in ZM, whereas low enzyme activity with high Na⁺ triggered the oxidative damage in XJD. These findings together illustrate the ion exchanger, antiporter, and lignin biosysthetic genes involving mechanistic insights into differential salt tolerance in alfalfa.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

Identification of QTL and candidate genes associated with biomass yield and Feed Quality in response to water deficit in alfalfa (Medicago sativa L.) using linkage mapping and RNA-Seq

Biomass yield and Feed Quality are the most important traits in alfalfa (Medicago sativa L.), which directly affect its economic value. Drought stress is one of the main limiting factors affecting alfalfa production worldwide. However, the genetic and especially the molecular mechanisms for drought tolerance in alfalfa are poorly understood. In this study, linkage mapping was performed in an F1 population by combining 12 phenotypic data (biomass yield, plant height, and 10 Feed Quality-related traits). A total of 48 significant QTLs were identified on the high-density genetic linkage maps that were constructed in our previous study. Among them, nine main QTLs, which explained more than 10% phenotypic variance, were detected for biomass yield (one), plant height (one), CP (two), ASH (one), P (two), K(one), and Mg (one). A total of 31 candidate genes were identified in the nine main QTL intervals based on the RNA-seq analysis under the drought condition. Blast-P was further performed to screen candidate genes controlling drought tolerance, and 22 functional protein candidates were finally identified. The results of the present study will be useful for improving drought tolerance of alfalfa varieties by marker-assisted selection (MAS), and provide promising candidates for further gene cloning and mechanism study.

Comparative physiological, transcriptomic, and WGCNA analyses reveal the key genes and regulatory pathways associated with drought tolerance in Tartary buckwheat

Drought stress is one of the major abiotic stress factors that affect plant growth and crop productivity. Tartary buckwheat is a nutritionally balanced and flavonoid-rich pseudocereal crop and also has strong adaptability to different adverse environments including drought. However, little is known about its drought tolerance mechanism. In this study, we performed comparative physiological and transcriptomic analyses of two contrasting drought-resistant Tartary buckwheat genotypes under nature drought treatment in the reproductive stage. Under drought stress, the drought-tolerant genotype XZSN had significantly higher contents of relative water, proline, and soluble sugar, as well as lower relative electrolyte leakage in the leaves than the drought-susceptible LK3. A total of 5,058 (2,165 upregulated and 2,893 downregulated) and 5,182 (2,358 upregulated and 2,824 downregulated) potential drought-responsive genes were identified in XZSN and LK3 by transcriptome sequencing analysis, respectively. Among the potential drought-responsive genes of XZSN, 1,206 and 1,274 genes were identified to be potential positive and negative contributors for XZSN having higher drought resistance ability than LK3. Furthermore, 851 out of 1,206 positive drought-resistant genes were further identified to be the core drought-resistant genes of XZSN based on WGCNA analysis, and most of them were induced earlier and quicker by drought stress than those in LK3. Functional annotation of the 851 core drought-resistant genes found that a large number of stress-responsive genes were involved in TFs, abscisic acid (ABA) biosynthesis, signal transduction and response, non-ABA signal molecule biosynthesis, water holding, oxygen species scavenging, osmotic adjustment, cell damage prevention, and so on. Transcriptional regulatory network analyses identified the potential regulators of these drought-resistant functional genes and found that the HD-ZIP and MYB TFs might be the key downstream TFs of drought resistance in Tartary buckwheat. Taken together, these results indicated that the XZSN genotype was more drought-tolerant than the LK3 genotype as evidenced by triggering the rapid and dramatic transcriptional reprogramming of drought-resistant genes to reduce water loss, prevent cell damage, and so on. This research expands our current understanding of the drought tolerance mechanisms of Tartary buckwheat and provides important information for its further drought resistance research and variety breeding.

Comparative physiological and root transcriptome analysis of two annual ryegrass cultivars under drought stress

Annual ryegrass is a widely cultivated forage grass with rapid growth and high productivity. However, drought is one of the abiotic stresses affecting ryegrass growth and quality. In this study, we compared the physiological and transcriptome responses of Chuansi No.1 (drought-tolerant, DT) and Double Barrel (drought-sensitive, DS) under drought stress simulated by PEG-6000 for 7 days. The results showed that Chuansi No. 1 had stronger physiological and biochemical parameters such as root properties, water content, osmotic adjustment ability and antioxidant ability. In addition, RNA-seq was used to elucidate the molecular mechanism of root drought resistance. We identified 8588 differentially expressed genes related to drought tolerance in root, which were mainly enriched in oxidation-reduction process, carbohydrate metabolic process, apoplast, arginine and proline metabolism, and phenylpropanoid biosynthesis pathways. The expression levels of DEGs were consistent with physiological changes of ryegrass under drought stress. We found that genes related to sucrose and starch synthesis, root development, osmotic adjustment, ABA signal regulation and specifically up-regulated transcription factors such as WRKY41, WRKY51, ERF7, ERF109, ERF110, NAC43, NAC68, bHLH162 and bHLH148 in Chuansi No. 1 may be the reason for its higher drought tolerance. This study revealed the underlying physiological and molecular mechanisms of root response to drought stress in ryegrass and provided some new candidate genes for breeding rye drought tolerant varieties.

Integrative small RNA and transcriptome analysis provides insight into key role of miR408 towards drought tolerance response in cowpea

Drought stress response studies and overexpression of vun-miR408 proved it to be essential for abiotic stress tolerance in cowpea. Small RNA and transcriptome sequencing of an elite high-yielding drought-tolerant Indian cowpea cultivar, Pusa Komal revealed a differential expression of 198 highly conserved, 21 legume-specific, 14 less-conserved, and 10 novel drought-responsive microRNAs (miRNAs) along with 3391 (up-regulated) and 3799 (down-regulated) genes, respectively, in the leaf and root libraries. Among the differentially expressed miRNAs, vun-miR408-3p, showed an up-regulation of 3.53-log₂-fold change under drought stress. Furthermore, laccase 12 (LAC 12) was identified as the potential target of vun-miR408-3p using 5' RNA ligase-mediated rapid amplification of cDNA ends. The stable transgenic cowpea lines overexpressing artificial vun-miR408-3p (OX-amiR408) displayed enhanced drought and salinity tolerance as compared to the wild-type plants. An average increase of 30.17% in chlorophyll, 26.57% in proline, and 27.62% in relative water content along with lesser cellular H₂O₂ level was observed in the transgenic lines in comparison with the wild-type plants under drought stress. Additionally, the scanning electron microscopic study revealed a decrease in the stomatal aperture and an increase in the trichome density in the transgenic lines. The expression levels of laccase 3 and laccase 12, the potential targets of miR408, related to lipid catabolic processes showed a significant reduction in the wild-type plants under drought stress and the transgenic lines, indicating the regulation of lignin content as a plausibly essential trait related to the drought tolerance in cowpea. Taken together, this study primarily focused on identification of drought-responsive miRNAs and genes in cowpea, and functional validation of role of miR408 towards drought stress response in cowpea.