Optimized combination methods for exploring and verifying disease-resistant transcription factors in melon

A large amount of omics data and number of bioinformatics tools has been produced. However, the methods for further exploring omics data are simple, in particular, to mine key regulatory genes, which are a priority concern in biological systems, and most of the specific functions are still unknown. First, raw data of two genotypes of melon (susceptible and resistant) were obtained by transcriptome analysis. Second, 391 transcription factors (TFs) were identified from the plant transcription factor database and cucurbit genomics database. Then, functional enrichment analysis indicated that these genes were mainly annotated in the process of transcription regulation. Third, 243 and 230 module-specific TFs were screened by weighted gene coexpression network analysis and short time series expression miner, respectively. Several TF genes, such as WRKYs and bHLHs, were regarded as key regulatory genes according to the values of significantly different modules. The coexpression network showed that these TF genes were significant correlated with resistance (R) genes, such as DRP2, RGA3, DRP1 and NB-ARC. Fourth, cis-acting element analysis illustrated that these R genes may bind to WRKY and bHLH. Finally, the expression of WRKY genes was verified by quantitative reverse transcription PCR (RT-qPCR). Phylogenetic analysis was carried out to further confirm that these TFs may play a critical role in Curcurbitaceae disease resistance. This study provides a new optimized combination strategy to explore the functions of TFs in a wide spectrum of biological processes. This strategy may also effectively predict potential relationships in the interactions of essential genes.