Genome-Wide Identification of WRKY Genes and Their Response to Cold Stress in Coffea canephora

WRKY transcription factors are known to play roles in diverse stress responses in plants. Low temperatures limit the geographic distribution of Coffea canephora Pierre ex A.Froehner. The WRKYs of C. canephora are still not well characterized, and the response of C. canephora WRKYs (CcWRKYs) under cold stress is still largely unknown. We identified 49 CcWRKYs from the C. canephora genome to gain insight into these mechanisms. These CcWRKYs were divided into three groups that were based on the conserved WRKY domains and zinc-finger structure. Gene expression analysis demonstrated that 14 CcWRKYs were induced during the cold acclimation stage, 17 CcWRKYs were preferentially upregulated by 4 °C treatment, and 12 CcWRKYs were downregulated by cold stress. Subsequently, we carried out a genome-wide analysis to predict 14,513 potential CcWRKY target genes in C. canephora. These isolated genes were involved in multiple biological processes, and most of them could be grouped by the response to stimulus. Among the putative CcWRKY target genes, 235 genes were categorized into response to the cold process, including carbohydrate metabolic, lipid metabolic, and photosynthesis process-related genes. Furthermore, the qRT-PCR and correlation analysis indicated that CcWRKY might control their putative targets that respond to cold stress. These results provide a basis for understanding the molecular mechanism for CcWRKY-mediated cold responses.