The auxin response factor (ARF) gene family in Cyclocarya paliurus: genome-wide identification and their expression profiling under heat and drought stresses

Auxin response factors (ARFs), as the main components of auxin signaling, play a crucial role in various processes of plant growth and development, as well as in stress response. So far, there have been no reports on the genome-wide identification of the ARF transcription factor family in Cyclocarya paliurus, a deciduous tree plant in the family Juglaceae. In this study, a total of 34 CpARF genes were identified based on whole genome sequence, and they were unevenly distributed on 16 chromosomes, with the highest distribution on chromosome 6. Domain analysis of CpARF proteins displayed that 31 out of 34 CpARF proteins contain a typical B3 domain (DBD domain), except CpARF12/ CpARF14/CpARF31, which all belong to Class VI. And 20 CpARFs (58.8%) contain an auxin_IAA binding domain, and are mainly distributed in classes I, and VI. Phylogenetic analysis showed that CpARF was divided into six classes (I–VI), each containing 4, 4, 1, 8, 4, and 13 members, respectively. Gene duplication analysis showed that there are 14 segmental duplications and zero tandem repeats were identified in the CpARF gene family of the C. paliurus genome. The Ka/Ks ratio of duplicate gene pairs indicates that CpARF genes are subjected to strong purification selection pressure. Synteny analysis showed that C. paliurus shared the highest homology in 74 ARF gene pairs with Juglans regia, followed by 73, 51, 25, and 11 homologous gene pairs with Populus trichocarpa, Juglans cathayensis, Arabidopsis, and rice, respectively. Promoter analysis revealed that 34 CpARF genes had cis-elements related to hormones, stress, light, and growth and development except for CpARF12. The expression profile analysis showed that almost all CpARF genes were differentially expressed in at least one tissue, and several CpARF genes displayed tissue-specific expression. Furthermore, 24 out of the 34 CpARF genes have significantly response to drought stress (P < 0.05), and most of them (16) being significantly down-regulated under moderate drought treatment. Meanwhile, the majority of CpARF genes (28) have significantly response to drought stress (P < 0.05), and most of them (26) are significantly down-regulated under severe drought treatment. Furthermore, 32 out of the 34 CpARF genes have significantly response to high, middle, and low salt stress under salt treatment (P < 0.05). Additionally, subcellular localization analysis confirmed that CpARF16 and CpARF32 were all localized to nucleus. Thus, our findings expand the understanding of the function of CpARF genes and provide a basis for further functional studies on CpARF genes in C. paliurus.