Transcript-wide identification and expression pattern analysis to comprehend the roles of AP2/ERF genes under development and abiotic stress in Trichosanthes kirilowii

Background The APETALA 2/ ethylene-responsive element binding factors ( AP2/ERF ), are thought to be associated with plant abiotic stress response, and involved in some plant hormone signaling pathways. Trichosanthes kirilowii is an important edible and medicinal crop, so far no research has been conducted on the TkAP2/ERF genes. Result In this study, a total of 135 TkERFs were identified, these genes were divided into 4 subfamilies and clustered into 13 groups. Moreover, 37 paralogous pairs were identified, with only two having Ka/Ks values greater than 1, proving that most TkERF genes underwent purifying selection during evolution. Co-expression networks constructed using transcriptome data at various flowering stages revealed that 50, 64, and 67 AP2/ERF genes correlated with members of the ethylene, gibberellin, and abscisic acid signaling pathways, respectively. When tissue cultured seedlings were treated with ETH, GA3 and ABA, 11, 12 and 17 genes were found to be up-regulated, respectively, suggesting that some members of the TkERF gene family may be involved in plant hormone signaling pathways. And under 4 ℃, PEG and NaCl treatment, 15, 20 and 19 genes were up-regulated, respectively, this suggested that these selected genes might be involved in plant abiotic stresses. Conclusions Overall, we identified 135 AP2/ERF family members, a comprehensive analysis of AP2/ERF gene expression patterns by RNA-seq and qRT-PCR showed that they played important roles in flower development and abiotic stress. This study provided a theoretical basis for the functional study of TkAP2/ERF genes and the genetic improvement of T. kirilowii .