Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 25 October 2022
Sec. Plant Abiotic Stress
This article is part of the Research Topic Environmental Extremes Threatening Food Crops View all 37 articles

Brassinosteroids induced drought resistance of contrasting drought-responsive genotypes of maize at physiological and transcriptomic levels

  • 1Gansu Provincial Key Lab of Arid Land Crop Science, College of Agronomy, Lanzhou, China
  • 2College of Agronomy, Hunan Agricultural University, Changsha, China
  • 3Crop Breeding Department, Jilin Changfa Modern Agricultural Science and Technology Group, co., Ltd., Changchun, China
  • 4College of Plant Protection, Gansu Agricultural University, Lanzhou, China
  • 5Department of Environmental Sciences, Kohsar University, Murree, Pakistan
  • 6Institute of Botany, University of the Punjab, Lahore, Pakistan
  • 7Department of Botany, Division of Science and Technology, University of Education, Lahore, Pakistan

The present study investigated the brassinosteroid-induced drought resistance of contrasting drought-responsive maize genotypes at physiological and transcriptomic levels. The brassinosteroid (BR) contents along with different morphology characteristics, viz., plant height (PH), shoot dry weight (SDW), root dry weight (RDW), number of leaves (NL), the specific mass of the fourth leaf, and antioxidant activities, were investigated in two maize lines that differed in their degree of drought tolerance. In response to either control, drought, or brassinosteroid treatments, the KEGG enrichment analysis showed that plant hormonal signal transduction and starch and sucrose metabolism were augmented in both lines. In contrast, the phenylpropanoid biosynthesis was augmented in lines H21L0R1 and 478. Our results demonstrate drought-responsive molecular mechanisms and provide valuable information regarding candidate gene resources for drought improvement in maize crop. The differences observed for BR content among the maize lines were correlated with their degree of drought tolerance, as the highly tolerant genotype showed higher BR content under drought stress.

Introduction

Maize (Zea mays L.), the most significant grain crop on the earth, is generally cultivated for food, feed, and biofuel production (Ansari et al., 2022). Climatic changes have significantly affected crop growth and yield in recent times (Fàbregas et al., 2018; Singhal et al., 2022). Moreover, different climate-induced abiotic stresses severely affect the growth and yield of major crops worldwide.

Among these abiotic constraints, drought is of prime importance because of its rapidly increasing severity in the coming future (Chen Y. et al., 2019). Drought stress inhibits maize seedling growth and crop development and reduces yield, particularly at the grain filling stage. However, yield losses caused by water scarcity in maize may vary depending on drought severity and the phenological stage of the crop (Hu et al., 2018). In this regard, plant growth-promoting hormones play a key role in safeguarding maize growth under drought stress, as evidenced by an earlier study, which showed that ABA levels speed up the kernel filling process (Jin et al., 2019).

The maize plant has a simple structure, and its unique germplasm makes it a more valuable crop than other cereals (Corso et al., 2015). Although significant progress has been made in investigating maize response to drought stress, previous studies have focused on physiological and metabolic or single-trait responses (Estrella-Maldonado et al., 2021). The high-throughput sequencing frameworks that rely on the RNA level and their analysis results have profoundly changed our viewpoint on the extent and heterogeneity of the maize transcriptome to a substantial degree over the last decade (Faghani et al., 2015). Whatever the case, most transcriptomic studies focused on the vegetative stage of the tested crop plants. Our perception of maize drought stress response structures and attributes highlights the need for further studies (Kosová et al., 2016).

Plants experience various environmental challenges during their life cycle that impact their growth and development (Anjum et al., 2017a). Drought is one of the most detrimental environmental conditions, eventually lowering agricultural yields (Anjum et al., 2016). During the vegetative phase, drought stress can limit growth rate, extend the vegetative development stage, and change the distribution of carbohydrates in maize (Aleem et al., 2021).

It has been demonstrated that short-duration water deficiencies result in 28%–32% losses in dry weight 66%–93% during the tasseling and ear formation stages of the maize crop, respectively. The number of kernel rows is established at the V9 stage when several ear shoots appear (Fracasso et al., 2016). The maize plant starts accumulating nutrients and dry weight quickly and steadily at the V10 stage and keeps doing so throughout the reproductive period (Hao et al., 2020). Long-term dryness during the pre-flowering stage has also been demonstrated to decrease the ultimate size of certain leaves and internodes, postpone the development of tassels and silk, and result in 15% to 20% yield losses (Haider et al., 2017).

Different agronomic and physiological techniques are used to reduce the negative effects of drought and to help plants develop a resistance to it (Janiak et al., 2018). The use of plant growth regulators is one of the potential methods to increase agricultural output under difficult conditions. Exogenously administered growth regulators have been shown to enhance plants’ ability to withstand a variety of abiotic stressors, including drought and heavy metal salt pressure (Tanveer et al., 2019). Brassinosteroids (BRs) serve various crucial functions in plant growth and belong to a novel family of phytohormones and development. On the steroidal side, each of these BRs has a 22R, 23R-diol structural characteristic. In light of this, molecules with 22R and 23R hydroxyls are much more active than their synthetic equivalents with 22S and 23S hydroxyls. In contrast, BR members with a 24S-methyl or -ethyl group exhibit 10 times more bioactivity in the native 22R, 23R configuration than do hormones with a 24R-alkyl function, highlighting the importance of the stereochemistry at this asymmetric center (Kang et al., 2019).

The steroidal phytohormones known as brassinosteroids (BRs) control the growth and development of plants. While the network-like pathways for BR production in Arabidopsis have been well characterized, it is still unclear what functions some biosynthetic enzymes play in the various sub-pathways (Kumar et al., 2018). Broad substrate specificities of the involved enzymes enable metabolite flow across a number of routes during BR production, which takes place along a metabolic grid. In the plant kingdom, three substances—cholesterol, sitosterol, and campesterol—are used to create BLs. The two sterols most prevalent in plant membranes are sitosterol and campesterol; however, there have been two methods for the biosynthesis of brassinolides that start with campesterol (Tanveer et al., 2018).

Similarly, around the time of pollination, 5 days of drought stress causes abnormal embryo development and a significant reduction in kernel number. The maize ear leaf makes a significant contribution to the buildup of biomass as a result of photosynthesis (Liu et al., 2019). Five or six leaves close to and above the ear generate most of the photosynthate for kernel production. The size of the “source” in plants is reduced due to drought stress, which causes a sharp reduction in the photosynthetic rate. Plants at the seedling stage were primarily used to survey and describe gene regulatory networks of the drought stress response (Wang et al., 2014).

Different plants face osmotic and drought stress, which leads to the expression of many genes. They may generally be separated into ABA-dependent and other signaling pathways. Since it affects stomatal closure and the production of stress-responsive genes, ABA is a key phytohormone in the drought stress response in plants. Hormones play a key role in plant response to abiotic stresses; osmotic stress signaling, the calcium-dependent route, mitogen-activated kinase-mediated signaling, phospholipid signaling, and reactive oxygen species (ROS) signaling are some other signaling mechanisms for drought stress tolerance (Anjum et al., 2017b).

In national plants, maize cross varieties are often preferred over homozygous regular lines since they are the most heterosis-rich crop in modern agricultural production. Separated from the normal line, the cross-variety plant is more rooted (Luo et al., 2018). The improvement steps are not performed in comparable field circumstances. There were not as many features passed on to the cultivars by the drought treatment as there would have been under normal conditions. Samples showed enormous rates of vegetation when exposed to the stress of the drought. Like all other samples, these exhibited the most pronounced leaves (Li et al., 2017).

While the more established leaves of different plants exposed to dry solid season areas showed signs of aging (the lower two leaves were dry or yellow, and the third leaf also started to senesce), this did not significantly affect other samples, which commonly developed under stress conditions (just toning down) (Liu et al., 2021). Our findings have improved the translation of maize drought response instruments and provided a theoretical foundation for breeding novel maize cultivars resistant to dry periods (Singh et al., 2017).

As normal biological changes occur, the significant financial difficulties caused by drought in agricultural productivity will get worse, and the financial difficulties caused by droughts in agricultural productivity will worsen (Dossa et al., 2017). Several cutting-edge technologies are being used and exploited to forecast these shifts (Mishra et al., 2022). It has been recommended to use polyethylene glycol (PEG), amino acids, cell protectors, phytohormones, minerals, typical eccentric treatments, and so on (Muthusamy et al., 2016). Brassinosteroids are a group of steroidal phytohormones included in combinations (BR). Due to their non-toxicity, non-mutagenicity, and eco-friendly nature, as well as their effectiveness at low concentrations, ease of use, and potential for deceptive mixing on a commercial scale, the BRs are a good choice for agricultural use (Zhao et al., 2020). Drought tolerance is a complex polygenic trait and is still not fully understood. This study aimed to investigate the drought tolerance of two contrasting maize lines. We performed the transcriptome analysis to identify the different genes controlling drought tolerance in maize and also checked the role of BRs in mitigating the drought-induced changes in plants.

Material and methods

Plant material, experimental design, and BR treatment

Seeds of two maize inbred lines (ILs), drought-sensitive (H21) and drought-tolerant (478), were provided by the Gansu Provincial Key Lab of Arid Land Crop Science, College of Agronomy, Gansu Agricultural University, Lanzhou, China, and were sown in pots. Two inbred lines were chosen to identify the tolerant genes and to check the Br role in drought tolerance. Br hormones also improve the plant growth in maize inbred lines; this showed that Br hormones could improve drought tolerance in sensitive lines, and tolerant genes identified in tolerant inbred lines can be used in molecular breeding. The pots were placed in a temperature-controlled environment in the laboratory (temperature 25 ± 5; relative humidity 65 ± 5; light/dark: 16/8 h). Up to 10 days following emergence, all the plants received enough water to maintain a steady moisture level; at that point, the application of treatments started. Treatments consisted of two levels of drought stress (20% and 0.3%) and three levels of BR (0, 25, and 50 ml). Fifty milliliters of 20% PEG 6000 was applied every 3 days to induce drought, and distilled water was applied in control. Seedlings at the four fully developed leaf stages were selected to determine morphological and physiological parameters. For transcriptome analysis, samples were stored at -80°C in liquid nitrogen. The experiment was carried out under a completely randomized design (CRD) with three replications. The BRs (including castasterone, 24-epicastasterone, brassinolide, and 24-epibrassinolide) have been extracted with ice-cold 80% aqueous methanol solution. The stock solutions of BRs were prepared at a concentration of 100 mg/l in methanol and stored at -18°C in the dark. Standards of lower concentrations were prepared weekly by the serial dilution of the stock solution with methanol and were stored at 4°C in the dark. During this analysis, a number of physiological and morphological parameters were analyzed, some of which involved growth parameters, for example, plant height, root dry weight, shoot dry weight, the number of leaves, leaf area, the specific mass of the fourth leaf, antioxidant activity, and malondialdehyde and chlorophyll contents.

Transcriptomic analysis

RNA extraction library construction and sequencing

The outright RNA was removed using a TRIzol reagent (Thermo Fisher, 15596018). The total amount of RNA and its quality were determined using the Bioanalyzer 2100 and the RNA 6000 Nano LabChip Kit (Agilent, California, United States of America, 5067-1511), and first-rate RNA tests with RIN numbers greater than 7.0 were used in the construction of the sequencing library. After the total RNA had been removed, the mRNA remaining in the pure RNA (5 μg) was cleaned using two rounds of filtering with Dynabeads Oligo (dT) (Thermo Fisher in California, USA). After being filtered, the mRNA was fragmented into shorter portions by utilizing divalent cations at a higher temperature (using the Magnesium RNA Fragmentation Module from NEB, cat. e6150, USA) for 5–7 min at 94°C. After that, the cut RNA pieces were inversely unraveled to make the cDNA using SuperScript™ II Reverse Transcriptase (Invitrogen, cat. 1896649, USA). The cDNA was then utilized to incorporate U-named second-deserted DNAs with E. coli DNA polymerase I (NEB, cat.m0209, USA) and RNase H (NEB, cat.m0297, USA) and dUTP Solution (Thermo Fisher, cat. R0133, USA). The unpolished terminations of each strand were given an A-base coating, which prepared the strands to be ligated to the recorded connectors.

Each connector had a T-base shade built into it to ligate the connection to the A-followed partitioned DNA. After affixing twofold document connections to the parts, AM Pure XP dabs were used to carry out size verification. After the power labile UD Genzyme (NEB, cat.m0280, USA) treatment of the U-named second-deserted DNAs, the ligated things were escalated with PCR using the following conditions: starting denaturation at 95°C for 3 min; eight examples of denaturation at 98°C for 15 s, treating at 60°C for 15 s, and increase at 72°C for 30 s; and a short-time later last extension at 72°C for 5 min. The preceding cDNA libraries had expansion sizes that ranged from 300 to 50 bp on average. Following the manufacturer’s instruction, we completed the 2 × 150 bp matched end sequencing (PE150) using an Illumina NovaSeq™ 6000 (LC-Bio Technology Co., Ltd., Hangzhou, China).

Sequence and filtering of clean reads

A cDNA library constructed by invention using the pooled RNA from the tests done with the Illumina NovaSeq™ 6000 succession stage was sequenced. We sequenced the transcriptome using the Illumina paired-end RNA-seq approach; reads obtained from the sequencing machines include raw reads containing adapters or low-quality bases, affecting the following assembly and analysis (Martin, 2011). Thus, to get high-quality clean reads, reads were further filtered by Cutadapt. The parameters were as follows: 1) removing reads containing adapters; 2) removing reads containing polyA and poly; 3) removing reads containing more than 5% of unknown nucleotides (N); 4) removing low-quality reads containing more than 20% of low-quality (Q-value ≤ 20) bases. After that, the quality of the sequence was checked using FastQC. It was observed that Q20, Q30, and GC content of the data had not been contaminated. Following that, a total of gigabase pairs worth of reads with clean ends was generated (Thompson et al., 2020). The raw sequencing data have been submitted to the NCBI Gene Expression Omnibus (GEO) databases or the NCBI Short Read Archive (SRA) with accession numbers. Both of these archives have been given accession numbers.

Pathway enrichment analysis (KEGG)

Qualities generally cooperate in assuming parts in specific organic capacities. Pathway-based investigation assists with understanding the rates of different physical processes. KEGG is a significant public pathway-related information base. Pathway improvement examination distinguished enhanced metabolic pathways or sign transduction pathways in DEGs in contrast with the entire genome foundation. Pathway advancement examination utilized Omic Share to determine extraordinarily upgraded metabolic pathways or sign transduction pathways in differentially communicated qualities. A hypergeometric test-taking FDR characterized extensively improved pathways in differentially communicated qualities (Aleem et al., 2021).

RT-PCR

The PCR assays were performed on the Rotor-Gene 3000 detection system (Corbett Research, Singapore) using One-Step PrimeScript™ RT-PCR Kit (Perfect Real Time) (Takara). All sets of reactions were carried out in a final volume of 20 l, each containing 10 l of 2× One-Step RT-PCR Buffer III, 0.4 l of MCMVf (10 M), 0.4 l of MCMVr (10 M), 0.8 l of probe (10 M), 0.4 l of Ex Taq™ (Takara) HS (5 U/l), 0.4 l of PrimeScript™ RT Enzyme Mix II, 1.0 l of total RNA or 1.0 l of RNA transcripts, and 6.6 l of RNase-free dH2O. Amplification reactions were performed as follows: 42°C for 5 min; 95°C for 10 s; 40 cycles of 95°C for 5 s, and 60°C for 20 s. The specificity of this TaqMan assay was evaluated using six different reactions, including water control. Using the TaqMan probe, strong fluorescent signals were detected only from reactions with samples, while the signals from four other samples along with the water control were superposed to the baseline under optimized reaction conditions. These samples can be differentiated from the four other maize samples by comparing the signals of different levels. The PCR products were analyzed further by agarose gel electrophoresis. The assay with the sample displayed the expected band of 67 bp and an unexpected faint band above the 67-bp band, whereas those of the four others did not.

Statistical analysis

Statistical significances between control and drought treatment were tested using the Newman–Keuls method at P < 0.05 by IBM® SPSS®, and data were generated by SigmaPlot 12.5. PCA was performed with the help of XLSTAT statistical software.

Results

Morphological and physiological traits of maize cultivars

During this analysis, a number of physiological and morphological parameters were analyzed, some of which involved growth parameters, for example, plant height, root dry weight, number of leaves, leaf area, and number of stems. Similarly, as far as functional parameters are concerned, they involve photosynthetic activity, leaf hydration, chlorophyll content, carbon dioxide diffusion, and CAT and pigment content. Phenotypic diversity for all traits was found under control and drought. Chlorophyll content and dry load were reduced by 90% during the drought, but a 40% reduction was seen in plant height. The two maize cultivars have been growing under identical circumstances, except for certain drought treatments. In both cultivars, as shown in the figure, there is a noticeable difference in plant height under two treatments; plants in drought treatment are shorter than those in water treatment. Both cultivars of drought-treated plants displayed withered leaves. However, under the two treatments, there was no visible variation in the shape of the ears between the two cultivars. When drought stress was applied to the samples under observation, they showed visible variations. According to the data obtained after analysis, root dry weight showed rough and symmetric results under the two conditions. Exceptions of plant height and any other variables’ phenotypic advantages were obtained in two situations. A couple of layouts look at each boundary where H21 displayed a more pronounced plant height. Additionally, while talking about root dry weight, 478 performs better than H21. A comparison of results for all boundaries is shown in Tables 1, 2.

TABLE 1
www.frontiersin.org

Table 1 Advanced test results of ANOVA.

TABLE 2
www.frontiersin.org

Table 2 ANOVA test result.

Pathway enrichment analysis (KEGG)

The results of BRs were significantly different between the two genotypes, and occasionally, H21 displayed more BRs in contrast to 478 in its leaves. The DEGs were subjected to KEGG pathway improvement investigation in the systems of 2 versus 0 h, 16 versus 2 h, and 16 versus 0 h. The paths for plant height and root dry weight, chlorophyll content, number of leaves, shoot dry weight, and many more were improved in the DEGs with wider articulation. Numerous DEGs were improved into signal transduction pathways following the treatment for drought stress, including the calcium flagging pathway-plant and plant chemical signal transduction. Contrarily, under the drought treatment, only a small number of DEGs independently enhanced the KEGG pathway in drought-sensitive maize (Figure 1) and drought-tolerant maize (Figure 2). The x-hub shows improvement in KEGG pathway examination of DEGs of drought-sensitive (Figure S1) and drought-tolerant lines (Figure 3) under the drought treatment, individually (Figure 4). The x-hub shows an improvement score, while the y-hub demonstrates the pathway name (Figure 5). The size of the speck addresses the quantity of DEGs. The various shades of specks address different p-values (Figure 6).

FIGURE 1
www.frontiersin.org

Figure 1 KEGG pathway enrichment scatter plot. The color of the point represents the size of the q-value. The smaller the q-value, the closer the point color to the red color. The point size expresses the number of differential genes in each pathway.

FIGURE 2
www.frontiersin.org

Figure 2 KEGG enrichment axial plot shows the DEG enrichment analysis results in the KEGG pathway.

FIGURE 3
www.frontiersin.org

Figure 3 Classification of differentially expressed genes. Total gene counts of KEEG terms are associated with all genes.

FIGURE 4
www.frontiersin.org

Figure 4 KEGG pathway enrichment scatter plot. The vertical axis represents the path name, and the horizontal axis represents the path factor corresponding to the rich factor.

FIGURE 5
www.frontiersin.org

Figure 5 KEGG pathway enrichment scatter plot. The vertical axis represents the path name, and the horizontal axis represents the path factor corresponding to the rich factor.

FIGURE 6
www.frontiersin.org

Figure 6 The differentially expressed genes or proteins are grouped hierarchically stretched GO terms, biological process, cellular components, and molecular functions.

Transcriptome analysis for drought tolerance

The results obtained for normal conditions and drought or comparatively different as the sample show different results under normal conditions and similarly different under drought stress. The point at which they are linked showed that 18,184 rates divided among every treatment went from 83.86% to 90.21%. Two hundred nine characteristics could only be found in the drought-resistant cultivar (H21) after receiving drought treatment. H21 displayed a better plant height under drought stress conditions than 478. This cultivar (H21) had more apparent leaves than 478 (Figure 7). Even though the more mature leaves of 478 plants that were subjected to dry-season areas of strength showed signs of aging (the first two leaves were dry or yellow, and the third leaf additionally began to show signs of aging), this did not have any significant impact on H21, which consistently grew (just slightly toning down) significantly even when subjected to stress conditions.

FIGURE 7
www.frontiersin.org

Figure 7 Transcriptomic profiling reveals shared signaling networks in different samples of maize.

Twelve cDNA libraries of the over two maize genotypes at the control and drought stress were built and sequenced with three organic duplicates to uncover the subatomic system of maize reaction to drought stress at the seedling stage. A sum of 20 G of clean information was obtained utilizing a reasonable stage. The complete base of each example was dispersed between 6.87 and 6.94 G, Q30 bases were circulated between 92.94% and 93.86%, and the typical GC content was 45.45%. By contrasting checks with the reference genome, the genomic arrangement of each example was obtained, and the arrangement rate was 96.25%–97.50% (Figure 7).

PCA

It has been seen in Figure 8 that PCA showed different values for H21 and different control values for 478, and the overall PCA value was found to be 13.98%, but on another note, they showed some values which are 4.87% for H21L0R1, 20.21% for control 478, and 10.67% for H21L2R3; similarly, values of 15.87%, 20.01%, 20.09%, 22.5%, and 17.5% for many other variables.

FIGURE 8
www.frontiersin.org

Figure 8 Categorization of different maize genotypes for drought tolerance and sensitivity. PCA with 78.95% values on the x-axis and 13.98% values on the y-axis.

The number of H21 plants in stressful circumstances was much higher than in the control plants. This improvement was less notable in H21 than in 478. The 478 line exhibited greater peroxidation of membrane lipids based on the PCA content than in genotype H21, while lipid peroxidation appeared to rise more in H21 than in 478 following drought exposure (Figure 9).

FIGURE 9
www.frontiersin.org

Figure 9 Categorization of different maize genotypes for drought tolerance and sensitivity. Principal component analysis of maize varieties based on different parameters visible in the form of colored spots with different values.

The vital part examination has been displayed in (Figure S2), which shows that drought additionally decreased the effectiveness of the essential photosynthetic cycles. There were obvious contrasts between the control and drought-pushed plants, especially for the boundaries portraying plant height, chlorophyll content, etc. In any case, no obvious contrasts were traced between the genotypes for the various limitations of the PCA examination.

With drought stress, the distribution of excess energy was effectively extended (boundaries, plant height, and dry root weight). According to the extent of the red spots, which were concentrated on plants of one or the other genotype, the specific mass of the fourth leaf did not seem to be negatively impacted during the drought. The drought negatively influenced the availability of the individual plant units, as inferred from the blue areas on the figure, but the two genotypes did not differ. Our samples remained unaffected by drought for dry root weight H21. In any event, 478 demonstrated a more noticeable decline in the rate of chlorophyll content due to drought at the end of the chain because of the observed differences in the location of the specific orange spots (Figure S3)

Similarly, although transcriptome analyses are widely used, they can only provide a limited picture of substantial changes in value clarity. There is no need to be concerned about changes in plant cell transcriptome brought by drought reflecting in the cell proteome (because of many advances and factors affecting quality verbalization among RNA and protein levels); at the same time, it was evident that there was no undeniable connection between these two strategies for substantial value verbalization information, as the drought’s administration of alterations in plant proteome was segregated from the transcriptome’s advancements. For further evaluation, it may be helpful to examine the relationship between the transcriptional, translational, and posttranslational degrees of articulation of BR-related qualities in plants familiar with the dry period and their effects on the components of intracellular BRs and the overall plant response to this strain component.

RT-PCR

RT-PCR gave us information to identify differentially expressed genes that are mapped; reads of the genes were calculated using the sequence alignment/map (SAM) files created by Tophat2. Then, differential expression analyses between two groups of samples were performed using the R package “DESeq2.” DEGs can be considered the best method for a graphical or visual representation of RT-PCR. Genes with |log2 (fold changes) |>1 and false discovery rate (FDR) <0.02 were recognized as differentially expressed. To detect the DEGs, four comparison groups, i.e., LMC_LMD (Lv28 under moderate drought control vs. Lv28 under moderate drought), LSC_LSD (Lv28 under severe drought control vs. Lv28 under severe drought), HMC_HMD (H21 under moderate drought control vs. 478 under moderate drought), and HSC_HSD (H21 under severe drought control vs. 478 under severe drought), were included as shown in Figure 10. Differential expression analysis was performed by comparing the gene expression profiles between different drought conditions to identify genes responding to drought in the two lines. In total, 100% and 51% of genes showed drought response under MD and SD in 478, respectively, among which 28% and 26% were upregulated under MD and SD, while 16% and 12% of genes were downregulated, respectively (Figure 11). In the drought-tolerant line H21, 68% and 33% of genes showed response to drought stress under MD and SD, respectively, among which 62% and 59% genes were upregulated, while 16% and 15% genes were downregulated, respectively (Figure 12). Remarkably, the number of the DEGs increased in 478 from MD to SD while it decreased in Lv28. In addition, the results of RT-PCR displayed as DEGs in H21 showed intensively upregulated under MD, with 82% genes upregulated but only 32% genes downregulated. Among these DEGs, 84% of genes responded to drought in H21 and 478 (Figure S4).

FIGURE 10
www.frontiersin.org

Figure 10 A variable peak with the second minor peak showing gene expression density; the gene density indicated the change in morphological parameters that were very clear during and after the whole analysis.

FIGURE 11
www.frontiersin.org

Figure 11 Results of different samples showing gene expression in the form of violin; the violin having different colors, indicating different samples analyzed at regular intervals.

FIGURE 12
www.frontiersin.org

Figure 12 Differently expressed genes displayed as a bar graph in the different groups; the two bars showed higher and lower gene expression values in maize samples.

Discussion

In the current work, we have presented comprehensive details of maize’s physiological and transcriptomic data under drought stress. Under simulated drought stress, it may be generalized that maize demonstrates stress tolerance mechanisms. The levels of activity arise following stress treatment, along with conductivity and soluble sugar content, suggesting that maize had improved scavenging abilities to withstand drought (Table 1). It was displayed that a rise in enzyme activity might remove peroxides produced by stress, preventing the plasma membrane’s oxidation and safeguarding cells against injury. According to our findings, maize controls the activity of defense enzymes, to reduce the accumulation of harmful substances (Zenda et al., 2018).

Additionally, we have also performed many analyses that gave us graphical results (Muthurajan et al., 2018). The results were incredibly significant and demonstrated that hormones could reduce growth under drought stress (Chen G. et al., 2019). Moreover, although transcriptome analyses are widely used, they only provide a very restricted result that causes changes (Table 3). According to the results of the KEEG improvement analysis, a total of the top 20 results were improved in the BRs of H21 and 478, including 15 in the cosmology of “organic cycle,” nine in the philosophy of “sub-atomic capacity,” and seven in the metaphysics of “cell component (Tiwari et al., 2021).” The DEGs of H21 and 478 had several advanced results, including eight organic cycle terms affecting the stress response (Chen G. et al., 2019).

TABLE 3
www.frontiersin.org

Table 3 Morphological and physiological trait analysis.

At the same time, if we talk about the results of our samples, i.e., H21 and 478, they did not show similar changes that were related to KEEG results (Xuan et al., 2022). Furthermore, all four sets of BRs showed that the “extracellular region” and “starch metabolic cycle,” showing two lines, included similar drought blockage components (Zenda et al., 2018). PCA revealed the identical characteristics of both plants, including the relevance of numerous limitations (Estrella-Maldonado et al., 2021). Each rate is shown with the eigenvalues in the PCA evaluation, indicating that these characteristics are accurate and connected to the other rates (Table 3).

Drought significantly reduced the growth of maize lines. We used seedlings of the maize line under controlled and drought conditions for a time-series transcriptomics study to identify DEGs and their patterns of response to increasing drought imposition (Zenda et al., 2019). In conclusion, drought stress inhibited enzyme activity linked to the production of cellular components, which is consistent with the suppression of plant growth (Yang et al., 2015). On the other hand, stress resistance mechanisms and regulatory activities are enhanced. We can evaluate dynamic gene responses to steadily rising drought stress using time-series transcriptomic data (Zeng et al., 2019). The samples, significance to drought tolerance, and interactions with sRNAs are other ways to define DEGs (Ghodke et al., 2020).

These results showed that the drought resistance of H21 may be increased by the root development, stability of the cytoskeleton, and performance of various cell cycles during drought. According to the KEGG enhancement analysis, 10 metabolic pathways were enhanced in the two lines (Wang et al., 2019). Starch and sucrose digestion, as well as plant chemical signal transduction, were fully developed in both lines. Still, their levels in H21L0R3 were typically greater than those in 478. Many researchers who worked with exogenous BR samples, transgenic plants, and earlier plants that had been kept with BRs before being exposed to lack of circumstances investigated the heights of numerous express records that were linked with drought reaction (Table 3).

Arabidopsis thaliana, Brassica napus, Cucumis sativus, Nicotiana tabacum, Solanum tuberosum, Brachypodium, or grains were the subjects of these tests (Ghodke et al., 2020). Even while the levels of CPD recordings dynamically increased over the first 24 h of PEG-induced osmotic stress in potatoes, it has not been thoroughly investigated whether drought-stressed plants exhibit increased BR biosynthetic properties or not (Wang et al., 2016). The early C-22 oxidation pathway and the late C-6 oxidation pathway both include several catalyzed stages (Qiu et al., 2019). In maize coding for H21, supplementary BR biosynthetic quality was dispersed.

This rate-limiting protein catalyzes different hydroxylation reactions at the beginning of early and late C-6 corrosion trials and the first H21. A vast phenotypic variability was discovered for all traits to handle the circumstances and start a dry phase (He et al., 2020). Chlorophyll content and the dry load decreased by 90% during the dry spell, while plant height decreased by 40%. Except for plant height, the variety ranges were smaller under drought stress compared to control circumstances, and the base was zero for all characteristics (Wang et al., 2015). Information obtained indicated that the properties were transported in a symmetrical manner (Tang et al., 2017).

In addition to identifying genes that may be involved in drought tolerance and other abiotic stresses, our research provides extensive transcriptome data enabling comparisons with gene responses to drought in other tissues under various drought circumstances (Zhang et al., 2017). Our results show that early and middle DEGs have greater genetic diversity and are more likely to be associated with drought resistance (Hu et al., 2018). The role of early and middle DEGs in drought tolerance involves huge levels of genetic variability. They include the gain and loss of gene function that may emerge from maize lines’ adaptability to changing growth circumstances (Jin et al., 2019).

A type of CNV called presence and absence variation (PAV) results from the lack of stress-responsive genes in the pathways that respond to environmental stressors under specific circumstances to reduce fitness costs. The discovery corroborates this that several crop species, including maize, soybean, and grape, have genes with CNV overrepresentation in the pathways of stress responses (Thirunavukkarasu et al., 2017).

Under drought circumstances, genetically speaking, the tolerant genotype H21 showed greater levels of metabolites, lower levels of lipid peroxidation, and better cell water retention than the sensitive genotype 478. In all, sample genomes analyzed with the help of RT-PCR in the form of DEGs were identified in our RNA-seq data as being expressed in response to drought (Zhang et al., 2019), with many DEGs being specifically discovered in H21. We discovered that genes linked to previously described pathways implicated in the drought stress response, such as those linked to the production of secondary metabolites, transcription factor regulation, detoxification, and stress defense, were changed upon exposure to drought stress.

The contribution of H21's characteristics associated with cell wall biosynthesis to the water maintenance of cell wall mix has regions of strength thanks to the input of a coordinated multienzyme complex in polysaccharide production. The raised articulation of glycosyltransferase added to the drought barrier of Arabidopsis thaliana cytokinin-inadequate freaks (Zenda et al., 2018). GO advancement examination uncovered that numerous H21 DEGs were connected with cell wall association and RT-PCR, particularly the expansion family prompted by different abiotic stresses and ABA. Stomatal thickness was displayed to diminish by overexpression. Furthermore, RT-PCR functional validation analysis supported the discovered genes’ differential expression. Our discoveries not only increase our understanding of the mechanisms that allow maize to withstand drought stress but also offer a priceless genetic resource or selection target for maize’s genetic advancement (Yang et al., 2015).

Additionally, GO and KEGG enrichment analyses on DEGs showed that several differentially expressed genes were most strongly enriched in biological processes, such as citrate cycle (TCA cycle), fatty acid metabolism, carbon fixation in photosynthetic organisms, and ribosome, enzymes distributed with a high proportion (Lei et al., 2022), indicating that a variety of metabolites were synthesized via active metabolic processes in leaves of maize. The RNA-seq technique can potentially identify novel genes from non-model organisms.

These observations have also provided a helpful resource for studying specific processes and pathways involved in growth and development of various plants and other plant species. The pathway-based analysis helps in comprehending the traits and subsequent real cycles. Maize has economic and conventional significance due to its distinct model and unique germplasm (Ereful et al., 2020). Additionally, it serves as a model plant for research on regulatory interactions. The investigation of maize drought tolerance has produced incredible advances in understanding this complex trait (Xia et al., 2020).

The results obtained after the analysis will help further research in the same area related to drought or any other similar plant that showed the same mechanism. Finding genes from maize that are connected to drought resistance and using them to increase the yield in various drought environments becomes a crucial breeding goal. It is theoretically and practically very important to analyze maize’s mechanism for drought tolerance (Vojta et al., 2016). Numerous regulatory mechanisms, including increased proline, soluble sugar, malondialdehyde, superoxide dismutase, potential operational delineations, and catchment area treatment content buildup were evolved in maize to deal with drought and to prevent or limit cellular damage induced by drought stress. Improved metabolic pathways were identified by pathway improvement assessment in distinct DEGs and the establishment of the entire genome. Omic Share was used for pathway headway evaluation to identify drastically altered metabolic or signal transduction pathways in several features (Tarun et al., 2020). Our study showed that applications of Br hormone could reduce the toxic effects in maize in inbred lines Br application increased the gene expressions and increased the drought tolerance in maize inbred lines.

Conclusions

The transcriptome analysis uncovered the differential expression of genes controlling drought tolerance in maize. Two contrasting maize genotypes performed differently under drought stress conditions. Drought stress leads to noticeable changes in maize genotypes. Results indicated that the drought-responsive maize genotype performed well under drought stress and as indicated by a differentially expressed gene network. DEGs are involved in different stress-responsive pathways and control several vital traits under drought stress. BRs played a key role in promoting plant growth under drought-stress conditions. The application of BRs improved the growth of maize genotypes and activated several stress-responsive genes. BRs compensated for the drought-induced toxic changes in maize genotypes and improved root growth, shoot growth, and chlorophyll contents. Using maize lines with reciprocal drought resistance, transcriptomes were examined under dry periods to identify several BRs. The BRs linked with plant chemical sign transduction, osmotic change, root advancement, and receptive oxygen rummaging were enhanced under drought lenient lines. Our outcomes showed drought-responsive components and gave new insights into understanding drought resistance in maize.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

SG and YP contributed to the study concept and design and statistical analysis. SG, AA, IU, and ZZ contributed to the analysis and interpretation of data. YP contributed to the investigation and resources. SG and SA contributed to the drafting of the manuscript. MK, AR, RS, and YW contributed to the review, editing, and proofreading of this manuscript. YP contributed to funding acquisition and study supervision. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the industrial support plan for colleges and universities of Gansu,China (No. 2022CYZC-46), the Fuxi Talent Project of Gansu Agricultural University, China (No. GAUFX-47102Y09), and the Lanzhou Science and Technology Bureau (No. 2020-RC-122).

Conflict of interest

The author Dr. Adnan Rasheed is employed by Jilin Changfa Modern Agricultural Science and Technology Group Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be constructed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.961680/full#supplementary-material

Supplementary Figure 1 | Venn diagram showing the relationship between the differentially expressed genes. The values in each circle represent the number of genes.

Supplementary Figure 2 | Principal Component Analysis (PCA) shows results using a scattered graph with 54.46% variance on the x-axis and 14.29% on the y-axis.

Supplementary Figure 3 | The expression numbers of transcriptomic responses in each treatment using a PCA, three peak lines showing variations in the analysis values.

Supplementary Figure 4 | Differently expressed genes displayed as Bars graph in different groups two bars showed gene expression one having a higher number of transcripts and other with lower.

Abbreviations

ANOVA, Analysis of variance; PCA, Principal component analysis; Q-PCR, Quantitative polymerase chain reaction; BRs, Brassinosteroids; KEGG, Kyoto Encyclopedia of Genes and Genomes; SOD, Superoxide dismutase; POD, Peroxidase; CAT, Catalase.

References

Aleem, M., Raza, M. M., Haider, M. S., Atif, R. M., Ali, Z., Bhat, J. A., et al. (2021). Comprehensive RNA-seq analysis revealed molecular pathways and genes associated with drought tolerance in wild soybean (Glycine soja sieb. and zucc.). Physiol. Plant 172, 707–732. doi: 10.1111/ppl.13219

PubMed Abstract | CrossRef Full Text | Google Scholar

Anjum, S. A., Ashraf, U., Tanveer, M., Khan, I., Hussain, S., Shahzad, B., et al. (2017a). Drought induced changes in growth, osmolyte accumulation and antioxidant metabolism of three maize hybrids. Front. Plant Sci. 8, 69. doi: 10.3389/fpls.2017.00069

PubMed Abstract | CrossRef Full Text | Google Scholar

Anjum, S. A., Ashraf, U., Zohaib, A., Tanveer, M., Naeem, M., Ali, I., et al. (2017b). Growth and developmental responses of crop plants under water deficit: a review. Zemdirbyste-Agriculture 104, 267–276. doi: 10.13080/z-a.2017.104.034

CrossRef Full Text | Google Scholar

Anjum, S. A., Tanveer, M., Ashraf, U., Hussain, S., Shahzad, B., Khan, I., et al. (2016). Effect of progressive drought stress on growth, leaf gas exchange, and antioxidant production in two maize cultivars. Environ. Sci. pollut. Res. 23 (17), 17132–17141. doi: 10.1007/s11356-016-6894-8

CrossRef Full Text | Google Scholar

Ansari, M. A., Bano, N., Kumar, A., Dubey, A. K., Asif, M. H., Sanyal, I., et al. (2022). Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) taub.) genotypes with contrasting drought tolerance. Funct. Integr. Genomics 1, 18. doi: 10.1007/s10142-022-00860-w

CrossRef Full Text | Google Scholar

Chen, Y., Chen, Y., Shi, Z., Jin, Y., Sun, H., Xie, F., et al. (2019). Biosynthesis and signal transduction of ABA, JA, and BRs in response to drought stress of Kentucky bluegrass. Int. J. Mol. Sci. 20, 1289. doi: 10.3390/ijms20061289

CrossRef Full Text | Google Scholar

Chen, G., Wang, Y., Wang, X., Yang, Q., Quan, X., Zeng, J., et al. (2019). Leaf epidermis transcriptome reveals drought-induced hormonal signaling for stomatal regulation in wild barley. Plant Growth Regul. 87, 39–54. doi: 10.1007/s10725-018-0450-0

CrossRef Full Text | Google Scholar

Corso, M., Vannozzi, A., Maza, E., Vitulo, N., Meggio, F., Pitacco, A., et al. (2015). Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance. J. Exp. Bot. 66, 5739–5752. doi: 10.1093/jxb/erv274

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, D., Bogeat-Triboulot, M. B., Tisserant, E., Balzergue, S., Martin-Magniette, M. L., Lelandais, G., et al. (2011). Comparative transcriptomics of drought responses in Populus: a meta-analysis of genome-wide expression profiling in mature leaves and root apices across two genotypes. BMC Genomics 11 (1), 1–21.

Google Scholar

Dossa, K., Li, D., Wang, L., Zheng, X., Liu, A., Yu, J., et al. (2017). Transcriptomic, biochemical and physio-anatomical investigations shed more light on responses to drought stress in two contrasting sesame genotypes. Sci. Rep. 7, 1–14. doi: 10.1038/s41598-017-09397-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ereful, N. C., Liu, L.-Y., Greenland, A., Powell, W., Mackay, I., Leung, H. (2020). RNA-Seq reveals differentially expressed genes between two indica inbred rice genotypes associated with drought-yield QTLs. Agronomy 10, 621. doi: 10.3390/agronomy10050621

CrossRef Full Text | Google Scholar

Estrella-Maldonado, H., Ramírez, A. G., Ortiz, G. F., Peraza-Echeverría, S., Martinez-De La Vega, O., Góngora-Castillo, E., et al. (2021). Transcriptomic analysis reveals key transcription factors associated to drought tolerance in a wild papaya (Carica papaya) genotype. PloS One 16, e0245855. doi: 10.1371/journal.pone.0245855

PubMed Abstract | CrossRef Full Text | Google Scholar

Fàbregas, N., Lozano-Elena, F., Blasco-Escámez, D., Tohge, T., Martínez-Andújar, C., Albacete, A., et al. (2018). Overexpression of the vascular brassinosteroid receptor BRL3 confers drought resistance without penalizing plant growth. Nat. Commun. 9, 1–13. doi: 10.1038/s41467-018-06861-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Faghani, E., Gharechahi, J., Komatsu, S., Mirzaei, M., Khavarinejad, R. A., Najafi, F., et al. (2015). Comparative physiology and proteomic analysis of two wheat genotypes contrasting in drought tolerance. J. Proteomics 114, 1–15. doi: 10.1016/j.jprot.2014.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Fracasso, A., Trindade, L. M., Amaducci, S. (2016). Drought stress tolerance strategies revealed by RNA-seq in two sorghum genotypes with contrasting WUE. BMC Plant Biol. 16, 1–18. doi: 10.1186/s12870-016-0800-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghodke, P., Khandagale, K., Thangasamy, A., Kulkarni, A., Narwade, N., Shirsat, D., et al. (2020). Comparative transcriptome analyses in contrasting onion (Allium cepa L.) genotypes for drought stress. PloS One 15, e0237457. doi: 10.1371/journal.pone.0237457

PubMed Abstract | CrossRef Full Text | Google Scholar

Haider, M. S., Zhang, C., Kurjogi, M. M., Pervaiz, T., Zheng, T., Zhang, C., et al. (2017). Insights into grapevine defense response against drought as revealed by biochemical, physiological and RNA-seq analysis. Sci. Rep. 7, 1–15. doi: 10.1038/s41598-017-13464-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, L.-Y., Liu, X.-Y., Zhang, X.-J., Sun, B.-C., Cheng, L., Zhang, D.-F., et al. (2020). Genome-wide identification and comparative analysis of drought related genes in roots of two maize inbred lines with contrasting drought tolerance by RNA sequencing. J. Integr. Agric. 19, 449–464. doi: 10.1016/S2095-3119(19)62660-2

CrossRef Full Text | Google Scholar

He, C., Du, Y., Fu, J., Zeng, E., Park, S., White, F., et al. (2020). Early drought-responsive genes are variable and relevant to drought tolerance. G3: Genes Genom. Genet. 10, 1657–1670. doi: 10.1534/g3.120.401199

CrossRef Full Text | Google Scholar

Hu, L., Xie, Y., Fan, S., Wang, Z., Wang, F., Zhang, B., et al. (2018). Comparative analysis of root transcriptome profiles between drought-tolerant and susceptible wheat genotypes in response to water stress. Plant Sci. 272, 276–293. doi: 10.1016/j.plantsci.2018.03.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Janiak, A., Kwasniewski, M., Sowa, M., Gajek, K., Żmuda, K., Kościelniak, J., et al. (2018). No time to waste: Transcriptome study reveals that drought tolerance in barley may be attributed to stressed-like expression patterns that exist before the occurrence of stress. Front. Plant Sci. 8, 2212. doi: 10.3389/fpls.2017.02212

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, H., Liu, S., Zenda, T., Wang, X., Liu, G., Duan, H. (2019). Maize leaves drought-responsive genes revealed by comparative transcriptome of two cultivars during the filling stage. PloS One 14, e0223786. doi: 10.1371/journal.pone.0223786

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Z., Babar, M. A., Khan, N., Guo, J., Khan, J., Islam, S., et al. (2019). Comparative metabolomic profiling in the roots and leaves in contrasting genotypes reveals complex mechanisms involved in post-anthesis drought tolerance in wheat. PloS One 14, e0213502. doi: 10.1371/journal.pone.0213502

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosová, K., Urban, M. O., Vítámvás, P., Prášil, I. T. (2016). “"Drought stress response in common wheat, durum wheat, and barley: transcriptomics, proteomics, metabolomics, physiology, and breeding for an enhanced drought tolerance,",” in Drought stress tolerance in plants (Cham: Springer), 277–314.

Google Scholar

Kumar, J., Gunapati, S., Kianian, S. F., Singh, S. P. (2018). Comparative analysis of transcriptome in two wheat genotypes with contrasting levels of drought tolerance. Protoplasma 255, 1487–1504. doi: 10.1007/s00709-018-1237-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lei, P., Liu, Z., Li, J., Jin, G., Xu, L., Ji, X., et al. (2022). Integration of the physiology, transcriptome and proteome reveals the molecular mechanism of drought tolerance in cupressus gigantea. Forests 13, 401. doi: 10.3390/f13030401

CrossRef Full Text | Google Scholar

Li, H., Li, M., Wei, X., Zhang, X., Xue, R., Zhao, Y., et al. (2017). Transcriptome analysis of drought-responsive genes regulated by hydrogen sulfide in wheat (Triticum aestivum l.) leaves. Mol. Genet. Genom. 292, 1091–1110. doi: 10.1007/s00438-017-1330-4

CrossRef Full Text | Google Scholar

Liu, S., Zenda, T., Dong, A., Yang, Y., Wang, N., Duan, H. (2021). Global transcriptome and weighted gene co-expression network analyses of growth-stage-specific drought stress responses in maize. Front. Genet. 12, 645443. doi: 10.3389/fgene.2021.645443

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Zhang, X., Sun, B., Hao, L., Liu, C., Zhang, D., et al. (2019). Genome-wide identification and comparative analysis of drought-related microRNAs in two maize inbred lines with contrasting drought tolerance by deep sequencing. PloS One 14, e0219176. doi: 10.1371/journal.pone.0219176

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, M., Zhao, Y., Wang, Y., Shi, Z., Zhang, P., Zhang, Y., et al. (2018). Comparative proteomics of contrasting maize genotypes provides insights into salt-stress tolerance mechanisms. J. Proteome Res. 17, 141–153. doi: 10.1021/acs.jproteome.7b00455

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17 (1), 10–12.

Google Scholar

Mishra, S., Sahu, G., Shaw, B. P. (2022). Integrative small RNA and transcriptome analysis provides insight into key role of miR408 towards drought tolerance response in cowpea. Plant Cell Rep. 41, 75–94. doi: 10.1007/s00299-021-02783-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Muthurajan, R., Rahman, H., Manoharan, M., Ramanathan, V., Nallathambi, J. (2018). Drought responsive transcriptome profiling in roots of contrasting rice genotypes. Indian J. Plant Physiol. 23, 393–407. doi: 10.1007/s40502-018-0381-9

CrossRef Full Text | Google Scholar

Muthusamy, M., Uma, S., Backiyarani, S., Saraswathi, M. S., Chandrasekar, A. (2016). Transcriptomic changes of drought-tolerant and sensitive banana cultivars exposed to drought stress. Front. Plant Sci. 7, 1609. doi: 10.3389/fpls.2016.01609

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, C.-W., Zhao, J., Chen, Q., Wu, F. (2019). Genome-wide characterization of drought stress responsive long non-coding RNAs in Tibetan wild barley. Environ. Exp. Bot. 164, 124–134. doi: 10.1016/j.envexpbot.2019.05.002

CrossRef Full Text | Google Scholar

Singh, D., Singh, C. K., Taunk, J., Tomar, R. S. S., Chaturvedi, A. K., Gaikwad, K., et al. (2017). Transcriptome analysis of lentil (Lens culinaris medikus) in response to seedling drought stress. BMC Genom. 18, 1–20. doi: 10.1186/s12864-017-3596-7

CrossRef Full Text | Google Scholar

Singhal, J., Nakamura, S., Ota, Y. (2022). Effect of epi-brassinolide on gram (Cicer arietinum) plants grown under water stress in juvenile stage. Indian J. Agri. Sci. 63, 395–397.

Google Scholar

Thompson, M. J. (2020). High-throughput SNP genotyping to accelerate crop improvement. PlantBreed. Biotech. 3, 195–212. doi: 10.9787/PBB.2014.2.3.195

CrossRef Full Text | Google Scholar

Tang, S., Li, L., Wang, Y., Chen, Q., Zhang, W., Jia, G., et al. (2017). Genotype-specific physiological and transcriptomic responses to drought stress in setaria italica (an emerging model for panicoideae grasses). Sci. Rep. 7, 1–15. doi: 10.1038/s41598-017-08854-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanveer, M., Shahzad, B., Sharma, A., Biju, S., Bhardwaj, R. (2018). 24-epibrassinolide; an active brassinolide and its role in salt stress tolerance in plants: a review. Plant Physiol. Biochem. 130, 69–79. doi: 10.1016/j.plaphy.2018.06.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanveer, M., Shahzad, B., Sharma, A., Khan, E. A. (2019). 24-epibrassinolide application in plants: An implication for improving drought stress tolerance in plants. Plant Physiol. Biochem. 135, 295–303. doi: 10.1016/j.plaphy.2018.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarun, J. A., Mauleon, R., Arbelaez, J. D., Catausan, S., Dixit, S., Kumar, A., et al. (2020). Comparative transcriptomics and co-expression networks reveal tissue-and genotype-specific responses of qDTYs to reproductive-stage drought stress in rice (Oryza sativa l.). Genes 11, 1124–1129. doi: 10.3390/genes11101124

CrossRef Full Text | Google Scholar

Thirunavukkarasu, N., Sharma, R., Singh, N., Shiriga, K., Mohan, S., Mittal, S., et al. (2017). Genomewide expression and functional interactions of genes under drought stress in maize. Int. J. Gen. Med. 20, 17. doi: 10.1155/2017/2568706

CrossRef Full Text | Google Scholar

Tiwari, P., Srivastava, D., Chauhan, A. S., Indoliya, Y., Singh, P. K., Tiwari, S., et al. (2021). Root system architecture, physiological analysis and dynamic transcriptomics unravel the drought-responsive traits in rice genotypes. Ecotoxicol. Environ. Saf. 207, 111252. doi: 10.1016/j.ecoenv.2020.111252

PubMed Abstract | CrossRef Full Text | Google Scholar

Vojta, P., Kokáš, F., Husičková, A., Grúz, J., Bergougnoux, V., Marchetti, C. F., et al. (2016). Whole transcriptome analysis of transgenic barley with altered cytokinin homeostasis and increased tolerance to drought stress. New Biotechnol. 33, 676–691. doi: 10.1016/j.nbt.2016.01.010

CrossRef Full Text | Google Scholar

Wang, L., Duan, C., Wu, D., Guan, Y. (2014). Quantification of endogenous brassinosteroids in sub-gram plant tissues by in-line matrix solid-phase dispersion-tandem solid phase extraction coupled with high performance liquid chromatography-tandem mass spectrometry. J. Chromatogr. A. 1359, 44–51. doi: 10.1016/j.chroma.2014.07.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Qin, Q., Sun, F., Wang, Y., Xu, D., Li, Z., et al. (2016). Genome-wide differences in DNA methylation changes in two contrasting rice genotypes in response to drought conditions. Front. Plant Sci. 7, 1675. doi: 10.3389/fpls.2016.01675

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Zenda, T., Liu, S., Liu, G., Jin, H., Dai, L., et al. (2019). Comparative proteomics and physiological analyses reveal important maize filling-kernel drought-responsive genes and metabolic pathways. Int. J. Mol. Sci. 20, 3743–3748. doi: 10.3390/ijms20153743

CrossRef Full Text | Google Scholar

Wang, N., Zhao, J., He, X., Sun, H., Zhang, G., Wu, F. (2015). Comparative proteomic analysis of drought tolerance in the two contrasting Tibetan wild genotypes and cultivated genotype. BMC Genomics 16, 1–19. doi: 10.1186/s12864-015-1657-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, H., Ma, X., Xu, K., Wang, L., Liu, H., Chen, L., et al. (2020). Temporal transcriptomic differences between tolerant and susceptible genotypes contribute to rice drought tolerance. BMC Genom. 21, 1–18. doi: 10.1186/s12864-020-07193-7

CrossRef Full Text | Google Scholar

Xuan, H., Huang, Y., Zhou, L., Deng, S., Wang, C., Xu, J., et al. (2022). Key soybean seedlings drought-responsive genes and pathways revealed by comparative transcriptome analyses of two cultivars. Int. J. Mol. Sci. 23, 2893. doi: 10.3390/ijms23052893

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Fountain, J. C., Wang, H., Ni, X., Ji, P., Lee, R. D., et al. (2015). Stress sensitivity is associated with differential accumulation of reactive oxygen and nitrogen species in maize genotypes with contrasting levels of drought tolerance. Int. J. Mol. Sci. 16, 24791–24819. doi: 10.3390/ijms161024791

PubMed Abstract | CrossRef Full Text | Google Scholar

Zenda, T., Liu, S., Wang, X., Jin, H., Liu, G., Duan, H. (2018). Comparative proteomic and physiological analyses of two divergent maize inbred lines provide more insights into drought-stress tolerance mechanisms. Int. J. Mol. Sci. 19, 3225. doi: 10.3390/ijms19103225

CrossRef Full Text | Google Scholar

Zenda, T., Liu, S., Wang, X., Liu, G., Jin, H., Dong, A., et al. (2019). Key maize drought-responsive genes and pathways revealed by comparative transcriptome and physiological analyses of contrasting inbred lines. Int. J. Mol. Sci. 20, , 1268. doi: 10.3390/ijms20061268

CrossRef Full Text | Google Scholar

Zeng, W., Peng, Y., Zhao, X., Wu, B., Chen, F., Ren, B., et al. (2019). Comparative proteomics analysis of the seedling root response of drought-sensitive and drought-tolerant maize varieties to drought stress. Int. J. Mol. Sci. 20, 2793. doi: 10.3390/ijms20112793

CrossRef Full Text | Google Scholar

Zhang, X., Liu, X., Zhang, D., Tang, H., Sun, B., Li, C., et al. (2017). Genome-wide identification of gene expression in contrasting maize inbred lines under field drought conditions reveals the significance of transcription factors in drought tolerance. PloS One 12, e0179477. doi: 10.1371/journal.pone.0179477

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Chen, B., Lu, G., Han, B. (2019). Overexpression of a NAC transcription factor enhances rice drought and salt tolerance. Biochem Bioph Res Co. 379, 985–989.

Google Scholar

Zhao, X., Bai, S., Li, L., Han, X., Li, J., Zhu, Y., et al. (2020). Comparative transcriptome analysis of two aegilops tauschii with contrasting drought tolerance by RNA-seq. Int. J. Mol. Sci. 21, 3595. doi: 10.3390/ijms21103595

CrossRef Full Text | Google Scholar

Keywords: GO analysis, KEGG, metabolism, transcriptome, zea mays

Citation: Gillani SFA, Zhuang Z, Rasheed A, Haq IU, Abbasi A, Ahmed S, Wang Y, Khan MT, Sardar R and Peng Y (2022) Brassinosteroids induced drought resistance of contrasting drought-responsive genotypes of maize at physiological and transcriptomic levels. Front. Plant Sci. 13:961680. doi: 10.3389/fpls.2022.961680

Received: 05 June 2022; Accepted: 20 September 2022;
Published: 25 October 2022.

Edited by:

Tanveer Alam Khan, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Germany

Reviewed by:

Mohsin Tanveer, University of Tasmania, Australia
Taiba Saeed, Aligarh Muslim University, India

Copyright © 2022 Gillani, Zhuang, Rasheed, Haq, Abbasi, Ahmed, Wang, Khan, Sardar and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yunling Peng, pengyunlingpyl@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.