A Comprehensive Gene Co-Expression Network Analysis Reveals a Role of GhWRKY46 in Responding to Drought and Salt Stresses

Abiotic stress, such as drought and salinity stress, seriously inhibit the growth and development of plants. Therefore, it is vital to understand the drought and salinity resistance mechanisms to enable cotton to provide more production under drought and salt conditions. In this study, we identified 8806 and 9108 differentially expressed genes (DEGs) through a comprehensive analysis of transcriptomic data related to the PEG-induced osmotic and salt stress in cotton. By performing weighted gene co-expression network analysis (WGCNA), we identified four co-expression modules in PEG treatment and five co-expression modules in salinity stress, which included 346 and 324 predicted transcription factors (TFs) in these modules, respectively. Correspondingly, whole genome duplication (WGD) events mainly contribute to the expansion of those TFs. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) analyses revealed those different modules were associated with stress resistance, including regulating macromolecule metabolic process, peptidase activity, transporter activity, lipid metabolic process, and responses to stimulus. Quantitative RT-PCR analysis was used to confirm the expression levels of 15 hub TFs in PEG6000 and salinity treatments. We found that the hub gene GhWRKY46 could alter salt and PEG-induced drought resistance in cotton through the virus-induced gene silencing (VIGS) method. Our results provide a preliminary framework for further investigation of the cotton response to salt and drought stress, which is significant to breeding salt- and drought-tolerant cotton varieties.