The NAC gene family is one of the largest families of transcriptional regulators in plants, and it plays important roles in the regulation of growth and development as well as in stress responses. Genome-wide analyses have been performed in diverse plant species, but there is still no systematic analysis of the NAC genes of Coffea canephora Pierre ex A. Froehner. In this study, we identified 63 NAC genes from the genome of C. canephora. The basic features and comparison analysis indicated that the NAC gene members increased via duplication events during the evolution of the plant. Phylogenetic analysis divided the NAC proteins from C. canephora, Arabidopsis and rice into 16 subgroups. Analysis of the expression patterns of CocNACs under cold stress and coffee bean development indicated that 38 CocNACs were differentially expressed under cold stress; six genes may play important roles in the process of cold acclimation, and four genes among 54 CocNACs showing a variety of expression patterns during different developmental stages of coffee beans may be positively related to the bean development. This study can expand our understanding of the functions of the CocNAC gene family in cold responses and bean development, thereby potentially intensifying the molecular breeding programs of Coffea spp. plants.

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.

Coffee is one of the most popular beverages around the world, which is mainly produced from the allopolyploid Coffea arabica. The genomes of C. arabica and its two ancestors C. canephora and C. eugenioides have been released due to the development of next generation sequencing. However, few studies on C. arabica are related to the PIN-FORMED (PIN) auxin efflux transporter despite its importance in auxin-mediated plant growth and development. In the present study, we conducted a genome-wide analysis of the PIN gene family in the three coffee species. Totals of 17, 9 and 10 of the PIN members were characterized in C. Arabica, C. canephora and C. eugenioides, respectively. Phylogenetic analysis revealed gene loss of PIN1 and PIN2 homologs in C. arabica, as well as gene duplication of PIN5 homologs during the fractionation process after tetraploidy. Furthermore, we conducted expression analysis of PIN genes in C. arabica by in silico and qRT-PCR. The results revealed the existence of gene expression dominance in allopolyploid coffee and illustrated several PIN candidates in regulating auxin transport and homeostasis under leaf rust fungus inoculation and the tissue-specific expression pattern of C. arabica. Together, this study provides the basis and guideline for future functional characterization of the PIN gene family.

WUSCHEL-related homeobox (WOX) genes play significant roles in plant development and stress responses. Difficulties in somatic embryogenesis are a significant constraint on the uniform seedling production and genetic modification of Coffea arabica, hindering efforts to improve coffee production in Yunnan, China. This study comprehensively analyzed WOX genes in three Coffea species. A total of 23 CaWOXs, 12 CcWOXs, and 10 CeWOXs were identified. Transcriptomic profile analysis indicated that about half of the CaWOX genes were actively expressed during somatic embryogenesis. The most represented CaWOXs were CaWOX2a, CaWOX2b, CaWOX8a, and CaWOX8b, which are suggested to promote the induction and development of the embryogenic callus, whereas CaWOX13a and CaWOX13b are suggested to negatively impact these processes. Co-expression analysis revealed that somatic embryogenesis-related CaWOXs were co-expressed with genes involved in embryo development, post-embryonic development, DNA repair, DNA metabolism, phenylpropanoid metabolism, secondary metabolite biosynthesis, and several epigenetic pathways. In addition, qRT-PCR showed that four WOX genes responded to cold stress. Overall, this study offers valuable insights into the functions of CaWOX genes during somatic embryogenesis and under cold stress. The results suggest that certain WOX genes play distinct regulatory roles during somatic embryogenesis, meriting further functional investigation. Moreover, the cold-responsive genes identified here are promising candidates for further molecular analysis to assess their potential to enhance cold tolerance.

Coffee is one of the most popular beverages around the world. As one of the best-known coffee species, Liberian coffee (Coffea liberica Bull ex Hiern 1876) has a high resistance to leaf rust, a devasting disease caused by Hemileia vastatrix. However, there are few reports on the systematic position and phylogenetic relationship of C. liberica at the chloroplast (cp) genome level. Thus, we successfully assembled its cp genome. The full length is 154,799 bp with a GC content of 37.48%. We have further annotated the cp genome and predicted 85 protein-coding genes together with 8 rRNAs and 37 tRNAs. Furthermore, a large single copy region (LSC), a small single copy region (SSC), an inverted repeat region a (IRa) and an inverted repeat region b (IRb) are identified with lengths of 84,868 bp, 18,121 bp, 25,905 bp and 25,905 bp, respectively. The phylogenetic tree indicates that C. liberica is closely related to C. canephora, which is consistent with a previous result obtained from genotyping‐by‐sequencing.