key: cord-0077905-gl6q7zl8
authors: Bano, Nasreen; Fakhrah, Shafquat; Mohanty, Chandra Sekhar; Bag, Sumit Kumar
title: Transcriptome Meta-Analysis Associated Targeting Hub Genes and Pathways of Drought and Salt Stress Responses in Cotton (Gossypium hirsutum): A Network Biology Approach
date: 2022-04-25
journal: Front Plant Sci
DOI: 10.3389/fpls.2022.818472
sha: 00dd8e31f731fa09569ab032b3a00784758c7acf
doc_id: 77905
cord_uid: gl6q7zl8

Abiotic stress tolerance is an intricate feature controlled through several genes and networks in the plant system. In abiotic stress, salt, and drought are well known to limit cotton productivity. Transcriptomics meta-analysis has arisen as a robust method to unravel the stress-responsive molecular network in crops. In order to understand drought and salt stress tolerance mechanisms, a meta-analysis of transcriptome studies is crucial. To confront these issues, here, we have given details of genes and networks associated with significant differential expression in response to salt and drought stress. The key regulatory hub genes of drought and salt stress conditions have notable associations with functional drought and salt stress-responsive (DSSR) genes. In the network study, nodulation signaling pathways 2 (NSP2), Dehydration-responsive element1 D (DRE1D), ethylene response factor (ERF61), cycling DOF factor 1 (CDF1), and tubby like protein 3 (TLP3) genes in drought and tubby like protein 1 (TLP1), thaumatin-like proteins (TLP), ethylene-responsive transcription factor ERF109 (EF109), ETS-Related transcription Factor (ELF4), and Arabidopsis thaliana homeodomain leucine-zipper gene (ATHB7) genes in salt showed the significant putative functions and pathways related to providing tolerance against drought and salt stress conditions along with the significant expression values. These outcomes provide potential candidate genes for further in-depth functional studies in cotton, which could be useful for the selection of an improved genotype of Gossypium hirsutum against drought and salt stress conditions.

Gossypium (cotton) species is an excellent model plant system for polyploid plant study (Rong et al., 2005) . The Gossypium genus comprises six tetraploids species and 45 diploids (Hawkins et al., 2006; Grover et al., 2015) . Approximately two million years ago (MYA), allotetraploid species (Gossypium hirsutum and Gossypium barbadense) were formed from the interspecific hybridization events between Gossypium herbaceum (A1) or Gossypium arboreum (A2) and Gossypium raimondii (D5) or Gossypium gossypioides (D6) followed by polyploidization (Senchina et al., 2003) . However, all the polyploid genomes are conserved in terms of gene content but they diversified at the subgenomic level (Chen et al., 2020) . The most widely available commercial species of cotton is upland cotton (G. hirsutum) because it provides >90% production of world cotton while G. barbadense gives 3-4% and G. herbaceum and G. arboreum give 2% (Chen et al., 2007) . G. hirsutum is the most crucial fiber-producing crop in the world and is grown in more than 80 countries.

Globally, India is the major producer of cotton, although yield productions are low due to crops suffering from one or more abiotic stress factors, such as salt, drought, heat, and waterlogging (Manikandan et al., 2020) . Abiotic stresses are a prime factor for limiting plant growth and yield productivity due to constant climate change. Understanding the molecular aspects in plants with response to abiotic stresses is at primacy now (Ambrosino et al., 2020) . These responses of plants are crucial to deal with environmental stresses in order to survive. The stress signal passes through the external environment, which continues to move toward the cell to trigger adaptive responses. Among the abiotic stress, drought and salinity show the most acute abiotic stress problems in the production of cotton in subtropic and tropic areas (Abdelraheem et al., 2021) . Drought stress majorly affects plant growth and causes morphological, physiochemical, and physiological changes in plants, which hamper plant growth and productivity (Chandra et al., 2021) , and salinity of soil also causes significant loss in the production of cotton (Manikandan et al., 2020) . Both drought and salt stress have negative impacts in all the stages of growth and can reduce cotton yield by up to 50-67% (Dabbert and Gore, 2014; Ullah et al., 2017; Abdelraheem et al., 2019 Abdelraheem et al., , 2021 . Thus, an understanding of the genomic level of cotton's response to various environmental stresses would give details on how to improve its genetic makeup with increased resistance to salt and drought stress. The availability of a large amount of transcriptome dataset deposition in Sequence Read Archive (SRA) provides an opportunity to accomplish a comprehensive gene network study (Liao et al., 2003) . Understanding these gene networks in the biological insights regarding abiotic stress response in plants is quite challenging (Hartwell et al., 1999; Cramer et al., 2011; Todaka et al., 2012) .

For gene regulation and construction of gene networks, transcriptome data is important for functional annotation and discovery of putative genes (Reményi et al., 2004; Ma et al., 2009) . Using the genes that were expressed in transcriptome data, gene networks are built based on similarity (correlation) in gene expression (coexpression). The co-expression networks (CNs) of candidate genes help in the study of expression patterns across multiple stress conditions among a large number of genes (Merico et al., 2009; Serin et al., 2016) . The Pearson correlation coefficient (PCC) calculated for the similarity score of CNs from genes with significant expression levels, such as sets of genes that have immensely correlated scores tend to group a small module. These correlated genes could be involved in associated signaling pathways and similar functions. Therefore, unknown or poorly characterized co-expressed genes in a module can be predicted by using genes of known function in certain biological processes (Mochida et al., 2011; Rhee and Mutwil, 2014) . In CNs of genes, transcription factor hubs are possible contenders for gene regulatory function. Transcription factors help to suggest the co-regulation of subsets of genes that are present in the same co-functionality or signal transduction pathways during the biological processes. At the transcriptional level, co-expression of genes relies upon the promoter region with similar cis-regulatory elements. So, coexpression profiles of genes are integrated into cis-motifs in the respective promoter sequences of a module, which give insights into the regulation of coordinated genes in CNs (Ma and Bohnert, 2007; Todaka et al., 2012; Sharma et al., 2013; Wilkins et al., 2016) .

High-throughput comprehensive analysis of expressed genes has become the important approach for showing the putative genes, finding cis-regulatory motifs, predicting gene function, and describing gene regulatory networks (Goda et al., 2008; Brady and Provart, 2009; Mochida et al., 2011) . Due to the availability of large-scale transcriptome data in the public domain, we have been able to construct biological networks and pathways. Much effort has been done to make use of co-expression analysis of transcriptome data. These studies used several approaches to find genes and predict the function of genes (Aoki et al., 2007) . Later, these pathways and networks have been used to understand functional interaction amongst genes (Moreno-Risueno et al., 2010) .

The large-scale availability of transcriptome data in NCBI has given us the chance to build a co-expression network for Gossypium transcriptome in different abiotic stress conditions. These networks might be used to determine the genes that have a significant role in salt and drought stress conditions in specific gene regulatory networks. In this study, we analyzed all the available transcriptome data in NCBI related to salt and drought stress till January 2021 in different tissues and cultivars. For the construction of the co-expression network of cotton, PCC was calculated for drought and salt stress conditions in different tissues (root, leaf, and seed) along with biological processes. Moreover, hub genes were identified from co-expressed modules and also retrieved their functional enrichment and pathways and these identified hub genes were also characterized through expression profiles. These investigations showed that a coexpression module study enables to find the putative genes that could be suitable targets for crop improvement approaches aimed to protect the cotton from drought and salt stress conditions.

A total of 118 (SRA accession number = SRP242577, SRP043419, SRP228010, SRP157859, and SRP192537) and 33 (SRA accession number = SRP192537 and SRP166405) RNA-seq datasets of treated and untreated cotton samples were retrieved from the NCBI SRA database till January 2021 using specific keywords "G. hirsutum" and individual abiotic stress. i.e., "salinity and drought." These datasets were extracted from the leaf, root, and seed tissues, and 4 cultivars (TM1, TM2, H15, and ZM12) in salinity and drought stress datasets were extracted from the leaf and root tissue that includes TM1 culture.

The parameters for the dataset taken into consideration for further study were the experiments which were conducted for drought and salinity treatment of cotton with at least two replicates (Walia et al., 2005 (Walia et al., , 2007 (Walia et al., , 2009 Jain et al., 2007; Hu et al., 2009; Mustroph et al., 2010) . In total, our final meta-datasets comprised 151 RNA-seq representing 33 droughts and 118 salt-treated in 7 experiments as summarized in Supplementary Table 1 . The quality control of these raw fastq datasets was processed by the Trimmomatic tool (Bolger et al., 2014) with average quality of 20 per base and adapter trimming was also done where it was required. Filtered reads were mapped on the G. hirsutum genome by using a STAR aligner (Dobin et al., 2013) with the default parameter. The quantification of differentially expressed genes (DEGs) was done through cuffdiff in Cufflinks version 2.2.1 software (Trapnell et al., 2012) with p-values ≤ 0.05, and log2fold-change (log2FC) values ≥ 2 for upregulated and ≤ −2 for downregulated genes. The top-ranked significant DEGs were compared in different tissue under drought and salt stress conditions. Unbiased clustering of drought and salt stress-responsive (DSSR) genes from different tissues with their log2FC were displayed through heatmap using the "Pheatmap" package function in R with normalized scale and euclidean distance method. 

Drought and salt stress-responsive genes were selected based on their differential expression in several experiments. Each tissue with higher differential expression was considered a reference in DSSR-GCN (gene co-expression network). The DSSR-GCN was constructed by the FPKM along with < 0.05 p-values. The "expression correlation networks" module in Cytoscape version 3.8 1 (Smoot et al., 2011) was used to construct the network. This module calculated a positive PCC (r ≥ 0.95) among the interacting members of the network. The co-expressed genes were visualized through Cytoscape.

To identify functionally grouped GO and pathway, the ClueGo module in Cytoscape was used (Bindea et al., 2009) . Using ClueGo, a biological process enrichment analysis of co-expressed genes in up and downregulated genes was performed in root and leaf tissue of drought and root, leaf, and seed tissues of salt stress datasets. Moreover, common enriched pathways between drought and salt were identified and represented through the EnrichmentMap module in Cytoscape. Node color shows the normalized enrichment score (NES) from the core enrichment genes. A positive score (red) shows gene set enrichment at the top of the ranked list a negative score (blue) indicates gene set enrichment at the bottom of the ranked list. The color of the border represents the p-value of the pathway enrichment from gene set enrichment analysis (GSEA). The node size shows the number (setSize) of DEGs identified in the pathways.

In a co-expression network, each edge shows an interaction between genes and the node shows a gene. We have estimated closeness connectivity, betweenness centrality (Cb), eccentricity, and degree centrality (k) using the network analyzer (Doncheva et al., 2012) plugin in Cytoscape software. Identification of hub genes was done based on connection or degree or they were shared in DSSR-GCN. Highly connected genes that have a high correlation with other genes in the network are known as hub genes. Any alteration in hub genes expression could lead to affect the major part of the network. Rapid information can be transferred in the network by the genes that have high connectivity compared to the other genes.

Significantly highly expressed top-ranked genes that have a higher betweenness centrality value and degree of connectivity (k) were taken as hub genes. Each node is allotted by a specific value based on links held through it, known as degree centrality (k). The degree helps in estimating the significance of node to control the network and it demonstrates the number of interactions retrained through a node in the network with other nodes.

Where K u denotes the node-set comprising entire neighbors of node u, and w (u,v) indicates the edge weight which connects node v with node u.

The betweenness centrality describes how often a node falls on the shortest path with other nodes in its neighborhood.

Betweenness centrality (Cb) measures how many times a node falls on the shortest path with other neighboring nodes. It refers to the ability of a node to control information flow and signal processing in the network.

Where, p (k,u,f) represents the number of interactions from node k to node f that pass through node u, whereas p (k,f) shows the total number of shortest interactions between nodes k and f (Prasad et al., 2020) .

For extensively analyzing the functional annotations of DEGs, we used the GO analysis toolkit for the agricultural community (agriGO version 2) (Tian et al., 2017) . The gene ontology provided the annotations at the cellular level, molecular level, and biological level. The GO graph was generated in RStudio using the ggplot package. Further, the functions of DSSR hub genes were confirmed through the GeneMANIA server (Warde-Farley et al., 2010) , which predicts the functions of gene sets. The functional study was performed at a 0.05 p-value cutoff, with the Arabidopsis thaliana database as a reference. For KEGG pathways analysis was employed using clusterProfiler package in RStudio with the number of permutation (nPerm) set to 10,000, org.At.tair.db and pAjustMethod set to "none." P-value ≤ 0.05 was considered for KEGG pathways. The PlantCARE database (Lescot et al., 2002) was searched for cis-regulatory elements in the 2 kb sequences upstream of identified hub genes identified by using the "Signal Scan Search" program.

To check the amplification of the key hub genes of cotton (G. hirsutum), plants were grown in the growth chamber under a controlled temperature of 30-37 • C with normal photoperiodic conditions. The seedlings were grown in separate triplicate pots for 1 month, stress treatment was given with, 20% PEG 8,000 (approx 300 ml/kg of soil) in the pots for the drought stress (Patade et al., 2009) and NaCl (300-mM) for salt stress . One time stress treatment was provided to each pot and then after intervals of 12, 24, 48, and 72 h RNA was isolated following the protocol (spectrum TM Plant Total RNA Isolation Kit). Once the RNA was isolated, cDNA was prepared (1 µg/µl) with verso cDNA synthesis kit (Thermo scientific). The amplification of the four putative genes for drought and four genes for salt stress was designed with primer3 Plus and checked with Quantitative Real-Time PCR (qRT-PCR) (HiMedia Insta Q 48 M4). The reaction was prepared with 10 µl of fast the SYBER TM green master mix (Applied Biosystem) and 1 µl of 100 ng cDNA adding 5 picomoles (forward and reverse) primers in a total 20 µl reaction. Ubiquitin was taken as internal control and the program for qRT-PCR was set at 95?C for 2 min, followed by cycling for 40 cycles of denaturation at 95 • C for 15 s, annealing at 60 C for 1 min, and extension at 72 C for 30 s. Relative expression of the employed genes was calculated with mean ± SD of biological triplicate samples by the 2ˆ− ct method (Livak and Schmittgen, 2001) .

The statistical analysis of qRT-PCR was carried out, using GraphPad Prism version 8.4.3 software, with two-tailed Student's t-test in triplicate sample repeats.

A total of 2 SRA accessions (SRP166405 and SRP192537) of drought and 5 SRA accessions (SRP242577, SRP043419, SRP228010, SRP157859, and SRP192537) of salt containing 33 and 118 fastq files were processed through the Trimmomatic tool. These RNA-Seq data included root, leaf, and seed of different cultivars. As a result, 2,168,052,606 clean reads with an average length of 277 bp were produced in salt stress data and 776,795,545 clean reads with an average length of 293 bp were generated in drought stress after removing the adapter sequences and lowquality reads ( Table 1) . The clean reads were then mapped to the G. hirsutum (cotton) genome, 1,960,725,723 (90.43%) reads in salt (Supplementary Table 1A ) and 722984665 (93.06%) reads in drought (Supplementary Table 1B ) were mapped to the reference cotton genome.

Further, 5,962 and 3,510 differentially expressed genes (DEGs) were identified in drought and salt stress data. Total 3,132 and 2,830 up and downregulated genes in drought and 2,265 and 1,245 up and downregulated genes (P-value < 0.05 and log2fold change > 2) in salt stress data. To make the summary of DSSR genes, which were expressed exclusively in drought or salt, DEGs were compared and found that 3,841 were unique in drought and 1,508 were unique in salt stress data (Figure 1) . The unbiased clustering of all significant DEGs of drought (leaf and root) and salt (leaf, root, and seed) stress-responsive genes from different tissues with their expression values (log2FC) 

To build the cotton-DSSR gene network, we have used the analyzed transcriptome data of cotton in two abiotic stress (drought and salt) conditions in different tissues. The proteinprotein interaction network was constructed in a circular layout after co-expression calculation using the Expression Correlation module in the Cytoscape tool. Further, among the analyzed transcriptome data, the top 100 DEGs genes with FPKM values and significant p-value (< 0.05) were selected from each tissue for further co-expression network analysis using the PCC method. Based on correlation score calculation droughtresponsive gene co-expression network (GCN) was generated FIGURE 1 | Venn diagram showing the number of shared and unique genes in differentially expressed genes (DEGs) of drought and salt stress data. In order to make drought and salt stress-responsive (DSSR) genes set compendium, we focused on 3,841 and 1,508 DSSR uniquely differentially expressed genes in drought and salt. which comprised 100 nodes connected by 1,704 and 2,747 coexpressed edges with minimum network density (ND) 0.34 and 0.55 in up-and downregulated genes of leaf tissue at PCC ≥ 0.95. Similarly, drought-responsive GCN of root tissue was generated which comprised 100 nodes connected by 1,505 and 2,314 co-expressed edges with minimum ND 0.30 and 0.46 in up and downregulated genes. Moreover, the salt responsive GCN of leaf tissue was generated which comprised 98 and 9 nodes connected by 4,565 and 45 co-expressed edges in up and downregulated genes, whereas 100 nodes connected with 3,252 and 2,604 co-expressed edges with minimum ND 0.65 and 0.52 in up-and down-regulated genes of salt root tissue. Similarly, salt responsive GCN of seed tissue was generated which comprised 27 and 74 nodes connected with 351 and 2,701 co-expressed edges in up and downregulated genes. These co-expressed genes were further used for network analysis of significant DEGs for drought and salt stress conditions. For this we have estimated closeness connectivity, betweenness centrality (Cb), eccentricity, and degree centrality (k) using network analyzer plugin in Cytoscape software and circular layout obtained network for drought (Figures 2A,C) and salt stress dataset at different tissue levels (Figures 3A,C,E) . The genes are ranked based on betweenness centrality scores and degree connectivity in drought (Supplementary Table 2A ) and salt stress dataset (Supplementary Table 2B ). Among these, SBT, EXPB2, CHS, TLP3, and CML44 have a high degree and betweenness centrality value in the drought stress dataset while WRKY6, ELF4, TLP, EXLB1, AATP1, and PER53 comprises high degree and betweenness centrality values in the salt stress.

To acquire an in-depth understanding of the drought and salt stress-responsive gene network, biological processes analysis has been applied on top selected differentially expressed genes using the ClueGo module with a p-value of 0.05 plugin in Cytoscape. Supplementary Table 3A for drought  and Supplementary Table 3B for salt datasets. In drought stress, the interacting genes in biological process (BP) were mainly enriched in the response to water deprivation, cellular response to abiotic stimulus, flavonoid biosynthetic process, phenylpropanoid metabolic process, response to jasmonic acid, and trehalose biosynthetic process in upregulated ( Figure 2B ) genes of leaf tissue while, flavonoid biosynthetic process, response to decreased oxygen levels, response to hydrogen peroxide and response to water deprivation in upregulated ( Figure 2D ) genes of root tissue. Moreover, in salt stress condition, upregulated BP were enriched with cellular response to abiotic stimulus, flavonoid metabolic process, isoprenoid metabolic process, hyperosmotic salinity response, response to abscisic acid, response to decreased oxygen levels, response to extracellular stimulus, and response to jasmonic acid in root tissue ( Figure 3B ) while regulation of secondary metabolite biosynthetic process, induced systemic resistance, oxylipin biosynthetic process and cellular response to boron-containing substance levels were upregulated BP in leaf tissue ( Figure 2D) . Moreover, upregulated BP in seed tissue were malate metabolic process peptidylthreonine phosphorylation, basipetalauxin transport, peptidylthreonine phosphorylation, and maintenance of protein location in the cell (Figure 3F) . The co-expression network and their functional enrichment study highlighted the key role of these putative genes in drought and salt stress responses in both (leaf and root) tissues except seed. Therefore seed tissue was dropped from further in-depth studies.

The gene Ontology is used to assign putative functions to all identified DEGs into three categories: biological processes (BP), cellular components (CC), and molecular functions (MF). Significantly differentially expressed genes were taken for further functional and biological pathways investigation. In drought stress conditions, gene ontology analysis in both the tissues (leaf and root) showed a significant presence of peroxidase activity, carboxylesterase, and petinesterase activity in molecular functions (MF). Biological processes in drought were enriched with response to abiotic stimulus, response to stress, and defense response (Supplementary Table 4A) . Similarly, in salt stress conditions, both tissues (leaf, and root) showed the significant presence of lipoxygenase activity, oxidoreductase activity, and cytokinin dehydrogenase activity, and phosphoric ester hydrolase activity in MF. BP in salt was enriched with the carboxylic acid biosynthetic process, response to stress, response to abiotic stimulus, response to oxidative stress, and cell redox homeostasis (Supplementary Table 4B ). In cellular component, apoplast, cell wall, and extracellular region were enriched in drought and cell wall, extracellular matrix, apoplast, cell cortex, and oxygenevolving complex were enriched in salt stress data.

Overall, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis reveals the role of drought stress response in Ubiquinone and other terpenoid-quinone biosynthesis and MAPK signaling pathways Table 5A ) and salt stress responses in phenylpropanoid biosynthesis and Isoquinoline alkaloid biosynthesis (Supplementary Table 5B ). Based on significant log2fold change (≤ -2 or ≥ 2) and p-value ≤ 0.05, TLP3, TLP5, TLP6, TLP, NCED1, EXLB1, PBS1, PCL1, BT2, TAR4, ITN1, ORG2, WRK54, and CAX3 genes were highly upregulated in drought stress condition and TLP1, TLP, GSTXA, NCS2, AKR1, LAC14, PDR1, ERD10, INV1, HSP70, DOX1, AOC, GPAT3, ATHB7, PER1, ELF4, EXLB1, RAV1, PER30, and CIPK6 were highly upregulated in salt stress condition. Overall, this study highlights the key role of these pathways in drought and salt stress responses. Additionally, the common differentially expressed genes (DEGs) were identified in both drought and salt stress conditions, and details of significant DEGs are listed in Supplementary Table 6. By utilizing these genes the common enriched pathways were identified and represented through Cytoscape using the EnrichmentMap module (Figure 4) . Biosynthesis of secondary metabolites, metabolic pathways, tryptophan metabolism, MAPK signaling pathways, protein processing in the endoplasmic reticulum, and phenylpropanoid biosynthesis pathways were commonly found in both drought and salt datasets.

In the co-expression network, the genes that have a high correlation with other genes are known as hub genes. Therefore, a topological feature of a network can be used to know the specific drought and salt stress-responsive genes. Using Cytoscape's Network analyzer tool, candidate hub nodes were determined using two topological features, betweenness (the fraction of all shortest paths that include a node within a network) and degree (number of connections) calculation. Altogether, we have identified 86 and 61 nodes with a high degree of connectivity and betweenness value and were subsequently referred to as hub genes in the network in drought (Table 1A) and salt stress conditions (Table 1B) .

A functional enrichment study of 86 and 61 putative drought and salt stress-responsive hub genes was performed by using GeneMania (GM) webserver and GO analysis. Through the GM webserver, we also predicted interactions between hub genes in the network using Arabidopsis thaliana as an additional parameter. The GM outcomes showed their role in response to water, oxidoreductase activity, acting on single donors with incorporation of molecular oxygen, isoprenoid metabolic process, dioxygenase activity, terpenoid biosynthetic process, isoprenoid biosynthetic process, and abscisic acid metabolic process in drought stress conditions ( Figure 5A ) while response to hypoxia, cellular response to decreased oxygen levels, cellular response to oxygen levels, response to water and aging in salt stress condition ( Figure 5B) . Furthermore, to confirm the function of GM, gene ontology classification has also been done. The GO analysis showed their role in response to stress, response to stimulus, response to abiotic stimulus, response to abscisic acid, and oxidoreductase activity in drought ( Figure 6A ). The highest number of genes were involved in response to stimulus, then in response to salt stress, catalytic activity and response to endogenous stimulus in salt stress condition ( Figure 6B) .

The KEGG pathways analysis demonstrated the role of drought stress response in phenylpropanoid biosynthesis, biosynthesis of secondary metabolites and diterpenoid biosynthesis, and sulfur metabolism pathways in which the highest number of genes are involved in the biosynthesis of secondary metabolites ( Figure 7A ) and salt stress, the higher the number of genes is involved in metabolic pathways than in phenylpropanoid biosynthesis and pyruvate metabolism ( Figure 7B ).

Based on degree centrality cutoff ≥ 30 and significant differential values (log2FC ≥ 2 or ≤ −2 with p-value ≤ 0.05) the expression profiles were illustrated for identified hub genes in both drought and salt stress conditions. The expression pattern revealed that NSP2, CDF1, DRE1D, ERF61, P2C08, HFB2A, TLP3, MYB15, ORG2, and TLP5 genes have significantly higher expression in drought stress conditions ( Figure 8A ). All the drought-responsive genes were majorly clustered in four groups, three of them showing upregulated genes, in which the TLP5 gene does not cluster with another gene. Moreover, TLP1, TLP, P2C08, EF109, ATHB7, ELF4, bHLH25, DRE1D, WRKY6, bHLH126, and NAC83 genes were expressed significantly higher in salt stress (Figure 8B) , in which Tubby-like protein 1 (TLP1) showed the highest expression (>10-fold) among the identified genes also making the separate cluster from other. These outcomes suggested that TLP might have some specific role in drought and salt stress conditions. Further, to validate the putative hub genes of G. hirsutum in response to drought and salt stresses, the qRT-PCR validation of four key genes in drought (NSP2, DRE1D, TLP5, and TLP3) and four in salt (TLP1, WRKY6, ATHB7, and EF109) stress conditions was carried out. The corresponding primers are listed in Table 2 . The expression of four putative genes was significantly upregulated in salt-stressed conditions, where the majority of the genes showed higher responses at 12, 24, and 72 h (Figures 9A-D) in terms of fold change (FC). TLP1 showed the highest expression at 24 h (>23 FC) and 12 h (>14 FC) compared with control ( Figure 9A) . Moreover, WRKY6 (>5 FC) and EF109 (>10 FC) also showed the highest expression at 12 h and 24 h compared with control ( Figures 9B,D) whereas ATHB7 showed the highest expression at 12 h (>12 FC) in comparison to control ( Figure 9C) . Likewise, the expression of four putative genes was also upregulated in drought-stressed conditions, where most of the genes showed higher expression at 12 and 72 h (Figures 9E,F) , where NSP2 showed the highest expression at 24 h (>16 FC) and 72 h (>29 FC) ( Figure 9E) and DRE1D, TLP5, and TLP3 showed relatively higher expression at 12 and 72 h compared with control (Figures 9F-H) . The correlation analysis of key drought (NSP2, DRE1D, TLP5, and TLP3) and salt (TLP1, WRKY6, ATHB7, and EF109) stress-responsive hub genes also represented a higher positive correlation between transcriptome and qRT-PCR dataset (Figures 9I,J) . Altogether, the differential responses of these putative hub genes in salt and drought stresses suggest that these genes may confer to resist the abiotic stresses in cotton.

The functions of highly expressed hub genes in the expression profiles were further confirmed through cis-regulatory elements analysis. We identified the binding frequency of cis-regulatory elements in 2,000 kb upstream sequences. In drought stress, TLP5, NSP2, DRE1D, NAC29, TLP3, and HFA6B comprises the higher number of abiotic stress-responsive elements (Table 3A) , and EF109, MIP1B, ATHB7, and TLP contains the higher number of abiotic stress-responsive elements in salt stress ( Table 3B ). The upstream region of hub genes showed the abundance of MBS (MYB binding site involved in drought-inducibility) in drought stress conditions (Figure 10A) , whereas ABRE (cis-acting element involved in the abscisic acid responsiveness) proximal elements were dominant in salt stress conditions ( Figure 10B) .

Cotton is a very important fiber-producing crop that is being cultivated across the globe. It is grown in semi-arid and arid regions where abiotic stresses namely salt and drought are highly prevalent due to which yields get affected. Among these, scarcity of water is a major limitation for the production of the crop. Water deficit conditions can severely lessen the crop yield and adversely influence the biochemical and physiological processes of plants which further lead to the lint yield reduction. The severity of the drought depends on various factors such as moisture storing ability of soils, rainfall amount and distribution, and evaporative demands. In cotton, the growth of the cell is affected by the reduction in turgor pressure due to drought stress. Drought also affects carbohydrate metabolism and photosynthesis directly or indirectly, changes in photosynthesis lead to the reduction of the boll maintenance in cotton (Leakey et al., 2009) . The important physiological defense mechanisms to cope with drought stress conditions are root development, stomata closure, photosynthesis, cellular adaptations, JA and ABA hormone production and ROS scavenging have been determined in cotton (Ullah et al., 2017) . Apart from water scarcity, soil salinity also gives rise to significant loss (5-9%) of cotton production (Nirmala et al., 2012) . Several defense mechanisms such as hormone regulation, membrane stability maintenance, antioxidants generation, and stress proteins induction, and carbon fixation rate have been found to be involved in plant survival under moisture stress conditions. Drought and salt stress conditions enhance the expression level of stress-related transcription factors and genes which further induces several droughts and salt stress-related pathways to stimulate the tolerance in the plant. There are several transcription factors such as WRKY, NAC, MYB, bHLH, ERF, ORG2, DREB, and TLP which have an important role in DSSR. To know the highly putative genes in drought and stress responses in cotton, comprehensive cotton-DSSR hub genes, and their functional enrichment studies were performed. The overall flow chart and strategy for designing and analyzing drought and salt stress-responsive gene co-expression networks (Figure 11) . By using the transcriptome (RNA-Seq) dataset of different tissues of cotton (G. hirsutum) in drought and salt stress conditions, the cotton-DSSR network study was done. Firstly, the GCN with the PCC method was constructed. Further, these GCN were taken for estimating the closeness connectivity, betweenness centrality (Cb), eccentricity, and degree centrality (k) using a network analyzer plugin with a circular layout. Furthermore, the genes are ranked based on betweenness centrality scores and degree connectivity (Figures 3, 4) . In drought, SBT, EXPB2, CHS, TLP3, and CML44 genes and in salt stress conditions, WRKY6, ELF4, TLP, EXLB1, AATP1, and PER53 genes were found with a higher degree of connectivity and betweenness centrality scores. Previously, it has been reported that phytosulfokine (PSK) precursor processing by subtilase SBT signaling improves drought stress tolerance in Arabidopsis (Stuhrwohldt et al., 2021) , expression of expansinlike B2 (EXPB2) enhanced for cell wall extensibility during drought stress conditions in Solanum pennellii (Egea et al., 2018) . Overexpression of Chalcone synthase (CHS) increases the production of flavonoids, which served as antioxidants in response to different stresses, namely drought, salt, and heavy metal . In Arabidopsis, TLP3 was found to function in response to ABA (Bao et al., 2014) and overexpression of calmodulin-like (CML44) stress-responsive gene from Solanum habrochaites improves tolerance to various abiotic stresses (Munir et al., 2016) . These genes (SBT, EXPB2, CHS, TLP3, and CML44) were also detected in drought stress conditions. Therefore, putative drought-responsive genes might be crucial for drought stress responses. In salt stress conditions, the WRKY6 expression was found to be enhanced in response to ROS treatment and it also triggers other plant defense responses and senescence processes (Robatzek and Somssich, 2001; Li et al., 2020) . The member of EARLY FLOWERING3 (ELF) have an important role in providing tolerant to high NaCl (Sakuraba et al., 2017; Cheng et al., 2020) and Tubby-like protein (TLP) members have been identified to play a crucial role in several abiotic responses such as drought, salt, oxidative (Wardhan et al., 2012; Xu et al., 2016; . Furthermore, the expression of expansin-like B1 (EXPB1) was found to be upregulated by salt, oxidative, and osmotic conditions Brassica rapa (Muthusamy et al., 2020) . The involvement of these genes (WRKY6, ELF, TLP, and EXLB1) demonstrated an imperative role in salt stress, and these genes, which were also detected in salt stress data, could therefore play an important role in salt stress. For an in-depth understanding of the DSSR network, the biological processes analysis of putative co-expressed responsive genes (Figures 2A,C,E) was carried out with the CluoGo module with a significant p-value. These results demonstrated response to water deprivation, cellular response to abiotic stimulus, flavonoid biosynthetic process, phenylpropanoid metabolic process, response to jasmonic acid and trehalose biosynthesis ( Figure 2B ) biological process in upregulated genes of leaf tissue while, flavonoid biosynthetic process, response to decreased oxygen levels, response to hydrogen peroxide and response to water deprivation ( Figure 2D ) biological process in upregulated genes of root tissue in drought stress condition. It has been reported that by regulating the homeostasis of ROS, flavonoids improve drought tolerance of maize seedlings and regulation of phenylpropanoid biosynthesis contributes to drought resistance in apple roots by MYB88 and MYB124, responsible for lignin accumulation under stress conditions (Geng et al., 2020) . Moreover, the accumulation of jasmonic acid is required for abscisic acid biosynthesis under drought stress conditions in citrus roots (de Ollas et al., 2013) and enhanced trehalose biosynthesis improves yield production under drought stress so there is a positive correlation between high-yield parameters and trehalose overproduction (Joshi et al., 2020) . The outcomes of biological processes in salt stress conditions showed cellular response to abiotic stimulus, flavonoid metabolic process, isoprenoid metabolic process, hyperosmotic salinity response, response to abscisic acid, response to decreased oxygen levels, response to extracellular stimulus, and response to jasmonic acid in upregulated genes of root tissue (Figure 3B) whereas, regulation of secondary metabolite biosynthetic process, induced systemic resistance, oxylipin biosynthetic process and cellular response to boron-containing substance levels were upregulated in leaf tissue ( Figure 3D) . Seed tissue was enriched with malate metabolic process, peptidyl-threonine phosphorylation, basipetalauxin transport, and maintenance of protein location in cells in upregulated genes ( Figure 3F ). An earlier report stated that MYB111 is a positive regulator in salt stress response and bioflavonoids addition able to rescue the loss of salt tolerance in myb111 mutants , and salt stress condition also enhances the biosynthesis of isoprenoid (Zuo et al., 2019) , indicating their importance in salt stress condition. Furthermore, oxylipin biosynthesis provides adaptation to plants under salt stress conditions (Savchenko et al., 2014) and TCA cycle intermediates (malate, α-ketoglutarate, aconitate, isocitrate, citrate, and succinate) enhanced during salt stress exposure in the barley cultivar Sahara (salt-sensitive) (Bandehagh and Taylor, 2020) . These overall outcomes suggest that high-ranked genes and biological processes of co-expressed networks play a crucial role in drought and salt stress responses at both (root and leaf) tissue levels.

Significantly differentially expressed genes were taken for further biological pathways investigation. The KEGG pathways study revealed the role of interacting genes in Ubiquinone and other terpenoid-quinone biosynthesis and MAPK signaling pathways in drought stress response (Supplementary Table 5A ) and phenylpropanoid biosynthesis and Isoquinoline alkaloid biosynthesis in salt stress response (Supplementary Table 5B ). Ubiquinone (UQ) is an important prenylquinones that acts as an antioxidant, participates in a plant's response to stress, and regulates gene expression and cell signal transduction (Liu and Lu, 2016) . The expression of genes in terpenoids biosynthesis is affected by drought stress and modulated by transcriptional regulation in trichomes of Nicotiana tabacum . Moreover, enhanced phenylpropanoid concentrations lead to improving the nutritional quality of wheat sprouts under salt stress conditions (Cuong et al., 2020) and increased alkaloid concentration plays an important role in the adaptation to salt conditions (Shi and Gu, 2020) . Overall, the KEGG pathway showed that DSSR genes mainly target secondary metabolitesrelated pathways to provide tolerance against drought and salt stress conditions. Moreover, the common pathways in both drought and salt stress datasets showed that a higher number of genes significantly involved in the biosynthesis of secondary metabolites, metabolic pathways, phenylpropanoid biosynthesis, and tryptophan metabolism pathways (Figure 4) . It is reported that tryptophan is a precursor for melatonin and serotonin biosynthesis (Murch et al., 2000) , which have a regulatory role in abiotic stress (drought and salt) (Kaur et al., 2015) . From outcomes of common enriched pathways, we can hypothesize that the involved genes might have an important role in providing resistance against drought as well as salt stress conditions in cotton.

Based on two topological features, betweenness (the fraction of all shortest paths that include a node within a network) and degree (the number of connections), candidate hub genes were determined using the network analyzer tool. A total of 86 and 61 hub genes were identified with a high degree of connectivity and betweenness value in drought (Table 1A) and salt stress conditions (Table 1B) . These identified genes mostly belong to the family of secondary metabolites responsive genes (SMRGs). SMRGs including MYB15 (Marchev et al., 2020) , NAC29 Wang T. T. et al., 2019) , TLP3 (Ibrahim and Jaafar, 2013; Bao et al., 2014) , HFB2A (Paupiere et al., 2020) , WRKY27, OLE16 (Shao et al., 2019) , BBE8 (Daniel et al., 2017) , LEA3 (Zhao et al., 2011) , CHSY (Dao et al., 2011) , DRE1D (Agarwal et al., 2017) , Malate synthase (de Kraker and Gershenzon, 2011) , cycloartenol synthase (Jin et al., 2017) and CASL1 (Taura et al., 2007) in drought enhances secondary metabolite responsive pathways and exert to provide tolerance in the adverse condition in drought stress condition. In salt, several identified genes such as TLP (Zare-Hassani et al., 2019; de Jesus-Pires et al., 2020) , WRKY6 (Robatzek and Somssich, 2001; Afrin et al., 2015) , MMP-2 , peroxidase 53 (Novo-Uzal et al., 2014) , bHLH25 (Goossens et al., 2019; Meraj et al., 2020; Yang et al., 2020) , malate synthase (Brito et al., 2020) , carboxylesterase Cao et al., 2019) , 1-deoxy-D-xylulose-5-phosphate synthase 2 (Paetzold et al., 2010) , NAC (He et al., 2019) , DCS2 (George et al., 2017) , XTH (Choi et al., 2011; Guerriero et al., 2018) , AtHB-7 (Olsson et al., 2004) , NCED5 (Udomchalothorn et al., 2017) , CAS (Shirazi et al., 2019) , ERF (Meraj et al., 2020) , galactinol synthase (Patel et al., 2020) , EXPA7 (Jadamba et al., 2020) and glutathione S-transferase (Nazar et al., 2011; Czerniawski and Bednarek, 2018) are SMRGs exhibited crucial role in providing resistance under salt stress condition. Overall, these putative hub genes might have putative functions in drought and salt stress responses in cotton.

A functional enrichment study of 86 and 61 candidate drought and salt stress-responsive hub genes was performed by using GeneMania (GM) webserver, further their functions were confirmed through gene ontology (GO) classification. The functional enrichment analysis of 86 candidate hub genes in drought stress by GM demonstrated their role in response to water, oxidoreductase activity, acting on single donors with incorporation of molecular oxygen, isoprenoid metabolic process, dioxygenase activity, terpenoid biosynthetic process, and abscisic acid metabolic process ( Figure 5A ). In rice, oxidoreductase activity implies a positive role in drought tolerance (Huang et al., 2014) , in transgenic tobacco, isoprenoid metabolism has an important function during drought (Tattini et al., 2014) and drought stress also affects the gene expression of terpenoids biosynthesis and control its transcriptional regulation in tobacco . At the molecular level, several physiological processes such as stomatal closure, germination, seed dormancy, and change in gene expression are regulated through ABA, which is a critical hormone, that serves as the core regulator of several stresses (drought, salt, and low temperature) in plants . These GM functions were confirmed through GO classification, in which molecular functions and biological processes showed similar outcomes such as response to stress, response to stimulus, response to abiotic stimulus, response to abscisic acid, and oxidoreductase activity in drought ( Figure 6A) . Moreover, in the functional enrichment study of 61 putative hub genes in salt stress conditions, the GM outcomes showed response to hypoxia, cellular response to decreased oxygen levels, cellular response to oxygen levels, response to water, and aging ( Figure 5B ). GM results revealed the response to decreased oxygen levels and hypoxia, which were caused by several factors such as excess nutrients, phosphorus, water body stratification (layering), and nitrogen, one of the main reasons is saline gradients and also showing the response to water, indicating their important role during salt stress condition. Salt stress-responsive functions of GM were further affirmed through the GO category, the biological process (BP) of GO demonstrating the significant direct role of these candidate genes in salt stress responses also showing significant catalytic activity in molecular function (MF) (Figure 6B) . The catalytic activity has been reported to be involved in abiotic stress tolerances by modulating genes (Chin et al., 2019) . Overall, a functional enrichment study in drought and salt stress conditions showed a significant critical role in both the stress conditions.

In drought stress, the KEGG pathways are enriched with phenylpropanoid biosynthesis, biosynthesis of secondary metabolites, diterpenoid biosynthesis, and sulfur metabolism pathways ( Figure 7A) . The phenylpropanoid biosynthetic pathway is triggered under several abiotic stresses such as drought, heavy metals, salinity, high/low temperatures, and ultraviolet radiations, resulting in the accumulation of several phenolic compounds that have the potential to scavenge harmful reactive oxygen species (Sharma et al., 2019) . Terpenoid phytoalexins comprise two zealexins and kauralexins families, accumulate in roots in response to both abiotic and biotic stress such as drought, high salinity, Fusarium verticillioides infection, and diabrotica balteata herbivory. In maize roots, the accumulation of terpenoids is associated with drought tolerance (Vaughan et al., 2015) . The previous report stated that sulfur not only acts like other macronutrients but is also imperative during metabolic adaptation reactions to drought stress, indicating sulfur metabolism is required in response to drought stress (Abuelsoud et al., 2016; Ahmad et al., 2016) . In salt stress, the KEGG pathways are enriched with metabolic pathways, amino sugar metabolism, phenylpropanoid biosynthesis, and pyruvate metabolism (Figure 7B) . The expression profile of phenylpropanoid biosynthetic pathway-related genes and phenylpropanoid compounds levels were investigated and found that enhanced phenylpropanoid concentrations lead to improve nutritional quality of wheat sprouts in the presence of salt treatment (Cuong et al., 2020) and there is a positive correlation of pyruvate pathways under salinity stress condition, as pyruvate pathways get activated under salt stress condition (Che-Othman et al., 2017) . Previously it was reported that the amino acid and carbohydrate metabolic pathways play a crucial role during salt stress in S. lycopersicum . These results demonstrated important signaling pathways that have crucial functions under drought and salt stress conditions.

The transcriptome differential expression profiles of the identified drought and salt stress-responsive hub genes were constructed using significant log2FC and degree centrality values. In drought stress, NSP2, CDF1, DRE1D, ERF61, P2C08, HFB2A, TLP3, MYB15, ORG2, and TLP5 genes show significantly higher expression ( Figure 8A ). The DEGs in drought were mainly clustered into four groups, in which Group II, III, and IV comprised upregulated genes, while Group I contained downregulated genes. In the upregulated genes, TLP5 does not group with any other genes, as group II only includes the TLP5 gene. In salt stress, the genes TLP1, TLP, P2C08, EF109, ATHB7, ELF4, bHLH25, DRE1D, WRKY6, bHLH126, and NAC83 show a significantly higher expression ( Figure 8B) . The DEGs in salt were also clustered into four groups, in which Group II, III, and IV comprised upregulated genes while Group I contained downregulated genes. In the upregulated genes, TLP1 does not group with other genes, since Group II only includes the TLP1 gene and showed the highest expression (>10-fold) among the upregulated genes. The TLP gene members have been identified in some dicotyledonous and monocotyledonous plants such as maize , apple (Xu et al., 2016) , and Arabidopsis, poplar, rice (Yang et al., 2008) . In Arabidopsis, AtTLP3, and AtTLP9 were found to have a redundant effect in response to osmotic stress treatments and ABA (Lai et al., 2004) , while AtTLP9 also responded to salt and drought stress (Lai et al., 2004; Bao et al., 2014) . In Malus domestica, several TLP genes showed higher expression during abiotic stress treatments (Xu et al., 2016) . It was found that CaTLP1 was expressed in response to dehydration stress in Cicer arietinum, while its expression led to increased tolerance to salt, drought, and oxidative stress in tobacco (Wardhan et al., 2012) . Based on the higher expression of hub genes in the transcriptome dataset, four genes (NSP2, DRE1D, TLP5, and TLP3) in drought and four genes (TLP1, WRKY6, ATHB7, and EF109) in salt stress were taken for further qRT-PCR validation (Figure 9) . The results showed that tubby-like protein (TLP) showed higher expression in both drought and salt stress responses. Further, correlation analysis also supported the significant positive correlation between qRT-PCR and transcriptome data (Figures 9I,J) . As per the previous report, the WRKY6 gene helps in the induction of salicylic acid (Robatzek and Somssich, 2001) in A. thaliana which has reported an important role in mitigating the adverse effects of salt stress on plants (Ma et al., 2017) . The higher expression of TLP1 in Cicer arietinum led to increased tolerance to salt, drought, and oxidative stress in tobacco (Wardhan et al., 2012) . Moreover, AtTLP3 member has been found to respond to abiotic stresses . The overall outcomes suggested that NSP2, DRE1D, TLP5, and TLP3 in drought and TLP1, WRKY6, ATHB7, and EF109 in salt may be crucial in cotton improvement.

Altogether, these studies suggested that TLP gene family members might play an important role in drought as well as in salt stress tolerance and could be suitable targets to protect the cotton from drought and salt stress conditions.

In order to confirm TLP's functional role and other highly expressed genes in drought and stress conditions, the cisregulatory element analysis was performed. In drought stress, TLP5, NSP2, DRE1D, NAC29, TLP3, and HFA6B include a higher number of abiotic stress-responsive elements (Table 3A) , the upstream region of these putative genes show the abundance of MBS (MYB binding site involved in drought-inducibility) ( Figure 10A ). This MBS-cis-element is required by the MYB TF for the gene expression of drought-inducible genes (Kaur et al., 2017) . Moreover, the TGA-element (auxin-responsive element) shows dominance in the TLP3 gene. The plant can regulate stomatal opening, integrate growth and stomatal regulation in drought conditions by utilizing auxin in Arabidopsis (Salehin et al., 2019) . In salt stress, EF109, MIP1B, ATHB7, and TLP contain a higher number of elements that respond to abiotic stress in salt stress (Table 3B) , where ABRE elements (cisacting element involved in the abscisic acid responsiveness) were abundant ( Figure 10B) . Previously it was stated that ABA is a crucial plant hormone with an important role in regulating salt stress responses (Ha et al., 2014) , salinity exposure enhances the accumulation of ABA, resulting in increased tolerance (Zhu, 2002) . In addition, TLP1 and WRKY comprise the abundance of P-box (gibberellin-responsive element) in their upstream region. Gibberellic acid is one of the essential plant growth regulators that help improve salt tolerance and also reduce the effects of salt stress on plants (Chauhan et al., 2019) and also induces high salinity tolerance in Pisum sativum through upregulation of antiporter genes, antioxidants, and secondary metabolites (Ahmad et al., 2021) . Altogether, these studies suggested that NSP2, DRE1D, TLP5, and TLP3 genes in drought and TLP1, WRKY6, ATHB7, and EF109 genes in salt stress might play an important role in cotton and could be suitable targets to protect the cotton crop from drought and salt stress condition.

Overall, our integrative transcriptome meta-analysis and its network topology analysis revealed that NSP2, DRE1D, ERF61, CDF1, and TLP3 genes induced strong drought stress responses, accentuated by the increased expression levels and the response of TLP1, TLP, EF109, ELF4, and ATHB7 genes enhanced during salt stress condition. The functional enrichment study of candidate hub genes showed the role of oxidoreductase and dioxygenase activity, response to an abiotic stimulus along with several secondary metabolite responsive genes in drought stress conditions, whereas the response to stimulus and salt, response to hypoxia, and decreased oxygen levels with catalytic activity responsive genes in salt stress condition. Moreover pathways analysis also demonstrating the higher numbers of putative hub genes involved in secondary metabolite responsive pathways that have found an important role in drought and salt stress states. These drought and salt stress responsive functions were further confirmed through the higher expression values of putative genes in transcriptome data and qRT-PCR validation and proximal element studies. These droughts (NSP2, DRE1D, TLP5, and TLP3) and salt (TLP1, WRKY6, ATHB7, and EF109) stress-responsive hub genes exert tolerance against unfavorable conditions and protect the cotton crop from drought and stress conditions. These protective actions of hub genes might account to reduce the loss of crop productivity and provide resistivity in adverse states.

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

NB carried out the bioinformatics analysis, designed, and drafted the manuscript. SF performed quantitative expression analysis. CM and SB participated to supervise the study. All authors have read and approved the final manuscript.

This work was supported by a grant (MLP0035) from the Council of Scientific and Industrial Research (CSIR), India. NB acknowledges the financial support from CSIR-UGC SRF.

We acknowledge the HPC facility, CSIR-4PI, and Bangalore. We are also grateful to Mehar Hasan Asif for giving an idea of this work. The manuscript number is provided by the institute is CSIR-NBRI_MS/2021/11/09.

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

Progress and perspective on drought and salt stress tolerance in cotton

GWAS reveals consistent QTL for drought and salt tolerance in a MAGIC population of 550 lines derived from intermating of 11 upland cotton (Gossypium hirsutum) parents

Sulfur metabolism and drought stress tolerance in plants

JA-mediated transcriptional regulation of secondary metabolism in medicinal plants

Dehydration responsive element binding transcription factors and their applications for the engineering of stress tolerance

Seed priming with gibberellic acid induces high salinity tolerance in Pisum sativum through antioxidants, secondary metabolites and upregulation of antiporter genes

Drought stress in maize causes differential acclimation responses of glutathione and sulfur metabolism in leaves and roots

Implications of abscisic acid in the drought stress tolerance of plants

Bioinformatics resources for plant abiotic stress responses: state of the art and opportunities in the fast evolving -Omics era

Approaches for extracting practical information from gene co-expression networks in plant biology

Can alternative metabolic pathways and shunts overcome salinity induced inhibition of central carbon metabolism in crops? Front

Characterization of Arabidopsis tubby-like proteins and redundant function of AtTLP3 and AtTLP9 in plant response to ABA and osmotic stress

ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks

Trimmomatic: a flexible trimmer for Illumina sequence data

Web-queryable large-scale data sets for hypothesis generation in plant biology

Overexpression of Ricinus communis L. malate synthase enhances seed tolerance to abiotic stress during germination

Genome-wide identification and functional analysis of carboxylesterase and methylesterase gene families in peach (Prunus persica L. Batsch). Front

Strategies to mitigate the adverse effect of drought stress on crop plantsinfluences of soil bacteria: a review

Influence of gibberellic acid and different salt concentrations on germination percentage and physiological parameters of oat cultivars

Chalcone synthase EaCHS1 from Eupatorium adenophorum functions in salt stress tolerance in tobacco

Genome-wide identification and comparative analysis of the TUBBY-like protein gene family in maize

Toward sequencing cotton (Gossypium) genomes

Genomic diversifications of five Gossypium allopolyploid species and their impact on cotton improvement

The soybean gene j contributes to salt stress tolerance by up-regulating salt-responsive genes

Connecting salt stress signalling pathways with salinity-induced changes in mitochondrial metabolic processes in C3 plants

Plant cytosolic ascorbate peroxidase with dual catalytic activity modulates abiotic stress tolerances

Constitutive expression of CaXTH3, a hot pepper xyloglucan endotransglucosylase/hydrolase, enhanced tolerance to salt and drought stresses without phenotypic defects in tomato plants (Solanum lycopersicum cv. Dotaerang)

Effects of abiotic stress on plants: a systems biology perspective

Effect of salinity stress on phenylpropanoid genes expression and related gene expression in wheat sprout

Glutathione S-transferases in the biosynthesis of sulfur-containing secondary metabolites in brassicaceae plants

Challenges and perspectives on improving heat and drought stress resilience in cotton

The family of berberine bridge enzyme-like enzymes: a treasure-trove of oxidative reactions

Chalcone synthase and its functions in plant resistance

Plant thaumatin-like proteins: function, evolution and biotechnological applications

From amino acid to glucosinolate biosynthesis: protein sequence changes in the evolution of methylthioalkylmalate synthase in Arabidopsis

Jasmonic acid transient accumulation is needed for abscisic acid increase in citrus roots under drought stress conditions

STAR: ultrafast universal RNA-seq aligner

Topological analysis and interactive visualization of biological networks and protein structures

The drought-tolerant Solanum pennellii regulates leaf water loss and induces genes involved in amino acid and ethylene/jasmonate metabolism under dehydration

Regulation of phenylpropanoid biosynthesis by MdMYB88 and MdMYB124 contributes to pathogen and drought resistance in apple

Drought and salt stress in Chrysopogon zizanioides leads to common and specific transcriptomic responses and may affect essential oil composition and benzylisoquinoline alkaloids metabolism

The AtGenExpress hormone and chemical treatment data set: experimental design, data evaluation, model data analysis and data access

Means and methods for regulating secondary metabolite production in plants

Molecular confirmation of species status for the allopolyploid cotton species, Gossypium ekmanianum Wittmack

Production of plant secondary metabolites: examples, tips and suggestions for biotechnologists

Investigation on salt-response mechanisms in Arabidopsis thaliana from UniProt protein knowledgebase

Positive regulatory role of strigolactone in plant responses to drought and salt stress

From molecular to modular cell biology

Differential lineage-specific amplification of transposable elements is responsible for genome size variation in Gossypium

The NAC protein from Tamarix hispida, ThNAC7, confers salt and osmotic stress tolerance by increasing reactive oxygen species scavenging capability

Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice

Comparative transcriptome sequencing of tolerant rice introgression line and its parents in response to drought stress

TaNAC29, a NAC transcription factor from wheat, enhances salt and drought tolerance in transgenic Arabidopsis

Abscisic acid induced changes in production of primary and secondary metabolites, photosynthetic capacity, antioxidant capability, antioxidant enzymes and lipoxygenase inhibitory activity of Orthosiphon stamineus Benth

Overexpression of rice expansin7 (osexpa7) confers enhanced tolerance to salt stress in rice

Fbox proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress

Two cycloartenol synthases for phytosterol biosynthesis in Polygala tenuifolia willd

Enhancing trehalose biosynthesis improves yield potential in markerfree transgenic rice under drought, saline, and sodic conditions

In-silico analysis of cis-acting regulatory elements of pathogenesis-related proteins of Arabidopsis thaliana and Oryza sativa

Regulatory roles of serotonin and melatonin in abiotic stress tolerance in plants

Molecular analyses of the Arabidopsis TUBBY-like protein gene family

Elevated CO2 effects on plant carbon, nitrogen, and water relations: six important lessons from FACE

PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences

The Arabidopsis MYB transcription factor, MYB111 modulates salt responses by regulating flavonoid biosynthesis

Flavonoids improve drought tolerance of maize seedlings by regulating the homeostasis of reactive oxygen species

GhTULP34, a member of tubby-like proteins, interacts with GhSKP1A to negatively regulate plant osmotic stress

Exogenous melatonin confers salt stress tolerance to watermelon by improving photosynthesis and redox homeostasis

Function and mechanism of WRKY transcription factors in abiotic stress responses of plants

Network component analysis: reconstruction of regulatory signals in biological systems

Plastoquinone and ubiquinone in plants: biosynthesis, physiological function and metabolic engineering

Analysis of relative gene expression data using real-time quantitative PCR and the 2-CT method

Enhanced tolerance to chilling stress in OsMYB3R-2 transgenic rice is mediated by alteration in cell cycle and ectopic expression of stress genes

Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression

Salicylic acid alleviates the adverse effects of salt stress on Dianthus superbus (Caryophyllaceae) by activating photosynthesis, protecting morphological structure, and enhancing the antioxidant system

Identification of relationship among exogenous NaCl with cotton leaves on cation uptake, nutrient ratios and status in rhizosphere soil

Green (cell) factories for advanced production of plant secondary metabolites

Transcriptional factors regulate plant stress responses through mediating secondary metabolism

How to visually interpret biological data using networks

Global landscape of a co-expressed gene network in barley and its application to gene discovery in Triticeae crops

Omics meet networks -using systems approaches to infer regulatory networks in plants

Overexpression of calmodulin-like (ShCML44) stress-responsive gene from Solanum habrochaites enhances tolerance to multiple abiotic stresses

Tryptophan is a precursor for melatonin and serotonin biosynthesis in in vitro regenerated St. John's Wort (Hypericum perforatum L. cv. Anthos) plants

Cross-kingdom comparison of transcriptomic adjustments to lowoxygen stress highlights conserved and plant-specific responses

BrEXLB1, a Brassica rapa expansin-like B1 gene is associated with root development, drought stress response, and seed germination

Understanding the significance of sulfur in improving salinity tolerance in plants

In vitro screening for salinity stress at seedling stage of cotton

Molecular cloning of two novel peroxidases and their response to salt stress and salicylic acid in the living fossil Ginkgo biloba

The homeobox genes ATHB12 and ATHB7encode potential regulators of growth in response to water deficit in Arabidopsis

The isogene 1-Deoxy-D-Xylulose 5-Phosphate Synthase 2 controls isoprenoid profiles, precursor pathway allocation, and density of tomato trichomes

Halopriming imparts tolerance to salt and PEG induced drought stress in sugarcane

Enhancing salt tolerance of plants: from metabolic reprogramming to exogenous chemical treatments and molecular approaches

Reprogramming of tomato leaf metabolome by the activity of heat stress transcription factor HsfB1

Targeting hub genes and pathways of innate immune response in COVID-19: A network biology perspective

Combinatorial control of gene expression

Towards revealing the functions of all genes in plants

A new member of the Arabidopsis WRKY transcription factor family, AtWRKY6, is associated with both senescence-and defence-related processes

Comparative genomics of Gossypium and Arabidopsis: unraveling the consequences of both ancient and recent polyploidy

Arabidopsis EARLY FLOWERING3 increases salt tolerance by suppressing salt stress response pathways

Auxin-sensitive Aux/IAA proteins mediate drought tolerance in Arabidopsis by regulating glucosinolate levels

Oxylipins and plant abiotic stress resistance

Rate variation among nuclear genes and the age of polyploidy in Gossypium

Learning from Co-expression networks: possibilities and challenges

New insights into the role of seed oil body proteins in metabolism and plant development

Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress

Recent advances in dissecting stress-regulatory crosstalk in rice

Transcriptome analysis and differential gene expression profiling of two contrasting quinoa genotypes in response to salt stress

Triterpenoid gene expression and phytochemical content in Iranian licorice under salinity stress

Cytoscape 2.8: new features for data integration and network visualization

Phytosulfokine (PSK) precursor processing by subtilase SBT3.8 and PSK signaling improve drought stress tolerance in Arabidopsis

Isoprene production in transgenic tobacco alters isoprenoid, non-structural carbohydrate and phenylpropanoid metabolism, and protects photosynthesis from drought stress

Cannabidiolic-acid synthase, the chemotype-determining enzyme in the fiber-type Cannabis sativa

Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks

OsNucleolin1-L expression in Arabidopsis enhances photosynthesis via transcriptome modification under salt stress conditions

Drought coping strategies in cotton: increased crop per drop

Accumulation of terpenoid phytoalexins in maize roots is associated with drought tolerance

Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage

Comparing genomic expression patterns across plant species reveals highly diverged transcriptional dynamics in response to salt stress

Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage

Genetic variance and transcriptional regulation modulate terpenoid biosynthesis in trichomes of Nicotiana tabacum under drought

Reference genome sequences of two cultivated allotetraploid cottons, Gossypium hirsutum and Gossypium barbadense

Transcriptional regulation of Lonicera japonica Thunb. during flower development as revealed by comprehensive analysis of transcription factors

Salt stress effects on secondary metabolites of cotton in relation to gene expression responsible for aphid development

The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function

Overexpression of CaTLP1, a putative transcription factor in chickpea (Cicer arietinum L.), promotes stress tolerance

EGRINs (environmental gene regulatory influence networks) in rice that function in the response to water deficit, high temperature, and agricultural environments

Genomewide identification and expression analysis of the tubby-like protein family in the Malus domestica genome

Photosynthetic regulation under salt stress and salt-tolerance mechanism of sweet sorghum

Genomewide comparative phylogenetic and molecular evolutionary analysis of tubby-like protein family in Arabidopsis, rice, and poplar

The effects of methyl jasmonate and salicylic acid on the production of secondary metabolites in organ culture of Ziziphora persica

The amino acid metabolic and carbohydrate metabolic pathway play important roles during salt-stress response in tomato

Group 3 late embryogenesis abundant protein in Arabidopsis: structure, regulation, and function

Salt and drought stress signal transduction in plants

Isoprene acts as a signaling molecule in gene networks important for stress responses and plant growth

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