Skip to main content

Genome-wide identification and characterization of the TIFY gene family in kiwifruit

Abstract

Background

The TIFY gene family is a group of plant-specific transcription factors involved in regulation of plant growth and development and a variety of stress responses. However, the TIFY family has not yet been well characterized in kiwifruit, a popular fruit with important nutritional and economic value.

Results

A total of 27 and 21 TIFY genes were identified in the genomes of Actinidia eriantha and A. chinensis, respectively. Phylogenetic analyses showed that kiwifruit TIFY genes could be classified into four major groups, JAZ, ZML, TIFY and PPD, and the JAZ group could be further clustered into six subgroups (JAZ I to JAZ VI). Members within the same group or subgroup have similar exon-intron structures and conserved motif compositions. The kiwifruit TIFY genes are unevenly distributed on the chromosomes, and the segmental duplication events played a vital role in the expansion of the TIFY genes in kiwifruit. Syntenic analyses of TIFY genes between kiwifruit and other five plant species (including Arabidopsis thaliana, Camellia sinensis, Oryza sativa, Solanum lycopersicum and Vitis vinifera) and between the two kiwifruit species provided valuable clues for understanding the potential evolution of the kiwifruit TIFY family. Molecular evolutionary analysis showed that the evolution of kiwifruit TIFY genes was primarily constrained by intense purifying selection. Promoter cis-element analysis showed that most kiwifruit TIFY genes possess multiple cis-elements related to stress-response, phytohormone signal transduction and plant growth and development. The expression pattern analyses indicated that TIFY genes might play a role in different kiwifruit tissues, including fruit at specific development stages. In addition, several TIFY genes with high expression levels during Psa (Pseudomonas syringae pv. actinidiae) infection were identified, suggesting a role in the process of Pas infection.

Conclusions

In this study, the kiwifruit TIFY genes were identified from two assembled kiwifruit genomes. In addition, their basic physiochemical properties, chromosomal localization, phylogeny, gene structures and conserved motifs, synteny analyses, promoter cis-elements and expression patters were systematically examined. The results laid a foundation for further understanding the function of TIFY genes in kiwifruit, and provided a new potential approach for the prevention and treatment of Psa infection.

Peer Review reports

Background

The TIFY family is a plant-specific gene family coding for transcription factors. The AT4G24470 gene in Arabidopsis thaliana was the first member of the TIFY gene family to be characterized and was previously known as ZIM (Zinc-finger protein expressed in Inflorescence Meristem) because it contains a C2C2-GATA zinc-finger structure [1, 2]. Due to presence of the highly conserved TIF[F/Y]XG (X represents any amino acid) motif in the ZIM domain-containing protein sequences, the ZIM family was later renamed and catalogued as the TIFY family [1]. The TIFY gene family could be further classified into four major subfamilies based on the different domain architectures and phylogenetic analyses, including TIFY, JAZ (jasmonate-ZIM-domain), PPD (PEAPOD) and ZML (ZIM/ZIM-like) [3]. The TIFY subfamily proteins contain a conserved TIFY domain, which is also shared by the other three subfamilies. In addition to the TIFY domain, the other three subfamilies contain specific and conserved domains. For example, the JAZ subfamily contains another conserved JA-associated (Jas, also named CCT_2) functional domain with a special consensus sequence SLX2FX2KRX2RX5PY near the C-terminus [4, 5]. The PPD subfamily contains a unique PPD domain in the N-terminus and a diverged Jas domain that lacks the conserved P and Y amino acids at the C-terminus region [4]. The ZML subfamily proteins contain a C-terminus C2C2-GATA zinc-finger DNA-binding domain and a CCT (CONSTANS, CO-like and TOC1) domain involved in protein-protein interaction [4].

The TIFY gene family plays a critical role in plant growth and development as well as in stress responses. The TIFY family genes are widely involved in regulating the development processes of plant organs and tissues such as stem, leaf and flower. For instance, studies have shown that AtTIFY1 (ZIM), the first characterized member of the TIFY family in Arabidopsis thaliana, was not only related to the development of inflorescence and flowering, but also promoted the petioles and hypocotyl extension by mediating cell elongation [2, 6]. Moreover, AtTIFY4a (PPD1) and AtTIFY4b (PPD2) enhanced leaf growth by regulating lamina size and limiting curvature of the leaf blade [7], and the AtTIFY4b (PPD2) gene was also involved in regulating leaf flatness and lateral organ development [8, 9]. In rice, overexpression of OsTIFY11b/OsJAZ10 could increase grain-size by enhancing accumulation and translocation of carbohydrates in the stems and leaf sheaths [10]. Another rice JAZ gene, OsTIFY3/OsJAZ1, was involved in regulation of spikelet development [11]. Arabidopsis thaliana with overexpression of a Jas-domain deletion version of AtJAZ1 (AtJAZ1ΔJas) displayed an early flowering phenotype under short day conditions, while overexpression of the CmJAZ1-like gene resulted in a late flowering phenotype in Chrysanthemum morifolium [12, 13]. Tomato plants overexpressing SlJAZ2 exhibited accelerated vegetative growth and early flowering [14].

Apart from the critical functions in plant growth and development, the TIFY family genes also play important regulatory roles in defense to various abiotic and biotic stresses. In terms of abiotic stresses, the TIFY gene family participates in response to drought stress, salt stress, alkaline stress and other abiotic stresses. For instance, in Arabidopsis thaliana, overexpression of AtJAZ7 conferred drought tolerance by regulating plant photosynthesis, redox, amino acids, phytohormones and defense metabolites [15]. In addition, OsJAZ1 was also involved in drought tolerance in rice by interacting with OsbHLH148 in the jasmonate signaling pathway [16]. In rice, OsTIFY11a-overexpressing plants could significantly improve tolerance to salt and dehydration stresses [17]. Overexpression of GsJAZ2 was demonstrated to improve tolerance to alkaline stress in soybean [18]. Under high salinity conditions, the germination and growth rate of Arabidopsis seedlings overexpressing wheat TdTIFY11a were higher than wild type, showing higher salt stress tolerance [19]. In terms of biotic stresses, JAZ proteins in Arabidopsis played a role in response to wounding and herbivory [20]. Furthermore, transgenic bread wheat lines over-expressing TaJAZ1 improved the resistance to powdery mildew by promoting the accumulation of reactive oxygen species [21]. In addition, the TIFY gene family has several other functions. For example, AsJAZ1 in Astragalus sinicus was involved in nodule development and nitrogen fixation [22], and AhJAZ1 and AhTIFY8 in Arachis hypogaea also participated in the root nodule symbiosis process [23]. The multiple regulatory effects of the TIFY gene family in plant growth and stress resistance indicate that this family contains a large number of valuable gene resources related to plant life activities and stress resistance responses, and thus it is of great significance to mine, identify and characterize these genes.

In recent years, a large number of plant genomes have been sequenced and released, which laid the foundation for the identification and characterization of TIFY gene family at the whole genome level. The TIFY family genes have been identified and characterized in at least 25 plant species, including Arabidopsis [1], Oryza sativa (rice) [17], Zea mays (maize) [24, 25], Tricicum aestivum (wheat) [19, 26, 27], Vitis vinifera (grape) [28], Malus × domestica (apple) [29], Solanum lycopersicum (tomato) [25, 30], Pyrus pyrifolia (sand pear) [31], Citrullus lanatus (watermelon) [32], and Camellia sinensis (tea) [33] (Fig. 1). A variable number of TIFY genes content were identified in these plant species, of which up to 77 TIFY genes identified in Brassica napus [34], while only 15 TIFY genes identified in Salvia militiorrhiza [35] and Citrullus lanatus [32] (Fig. 1). Most of the identified TIFYs in these species, especially those in eudicots, could be clustered into four subfamilies (TIYF, JAZ, ZML, PPD). The TIFYs in monocots plants seemed to be only grouped into three subfamilies, and no PPD subfamily members were detected in monocot species (Fig. 1). With the identification of more TIFY genes in plants, the understanding of TIFY genes is gradually deepened, and the biological functions of TIFY genes are gradually clarified. However, the evolution and functional divergence of the kiwifruit TIFY genes remained unclear until now.

Fig. 1
figure 1

The number variations and compositions of TIFY gene families in different plant genomes. The plant species evolutionary relationships were redrawn according to APG (Angiosperm Phylogeny Group) IV

Kiwifruit is a popular fresh fruit consumed worldwide, with important nutritional and economic value. Because the fruits are remarkably rich in vitamin C, dietary fiber, mineral elements and other nutrients, kiwifruit is also well known as ‘the king of fruits’. The TIFY gene family has a variety of important biological functions and plays a very important regulatory role in plant growth and development and stress resistance. Therefore, it is of great significance to identify and clarify the biological functions of TIFY gene family in kiwifruit. Currently, several kiwifruit genomes have been released, such as Actinidia chinensis ‘Hongyang’ [36, 37], A. chinensis ‘red5’ [38] and A. eriantha ‘White’ [39], which will be helpful in identifying and characterizing kiwifruit TIFY gene family.

In this study, genome-wide identification and investigation of the TIFY family genes were carried out from two assembled kiwifruit genomes, including A. chinensis ‘Hongyang’ genomes (v3.0) and A. eriantha ‘White’ genome. Furthermore, gene structures, conserved domains, phylogenetic analysis, chromosomal locations, cis-regulatory elements, gene synteny analyses and expression characteristics were subsequently analyzed. The results of this study will pave a way to further understanding the evolution and biological functions of kiwifruit TIFY genes.

Results

Identification of TIFY genes in kiwifruit

To identify all the putative TIFY genes in kiwifruit, the seed profile of the TIFY domain (PF06200) was used to search against the annotated proteins of A. eriantha, and A. chinensis. All the putative members were verified for the presence of the conserved TIFY domain, and finally, a total of 27 and 21 TIFY genes were identified in the genome of A.eriantha and A. chinensis, respectively (Fig. 1). The CDS (coding sequence) and amino acid lengths of these identified TIFY sequences varied extensively. For example, the amino acid sequence lengths of AeTIFYs and AcTIFYs varied from 112 aa (AePPD1) to 718 aa (AeJAZ13) and from 107 aa (AcJAZ13) to 764 aa (AcZML1), respectively (Additional file 1). The predicted molecular weight of AeTIFYs and AcTIFYs also varied greatly, and ranged from 12.3 kDa (AePPD1) to 81.2 kDa (AeJAZ13) and from 12.2 kDa (AcJAZ13) to 84.1 kDa (AcZML1), respectively (Additional file 1). The theoretical pI values of 19 out of 27 TIFYs in A. eriantha and 18 out of 21 TIFYs in A. chinensis were higher than 7.0, indicating that most of the kiwifruit TIFYs were alkaline proteins. Subcellular localization analysis showed that, except for the AcZML1 in A. chinensis predicted to be located in the mitochondrion, all other kiwifruit TIFYs were predicted to be located in the nucleus (Additional file 1). The identified kiwifruit TIFY genes were named on the basis of their chromosomal positions and the phylogenetic relationships with the TIFY genes in Arabidopsis and rice described below. The detail information of these kiwifruit TIFY family genes is available in Additional file 1.

Phylogenetic analyses and classification of the kiwifruit TIFY gene family

To figure out the classification and evolutionary relationships of the identified kiwifruit TIFY gene family, a phylogenetic tree was constructed based on the multiple sequence alignment of 86 TIFY protein sequences, including 27 TIFYs in A. eriantha, 21 in A. chinensis, 18 in Arabidopsis thaliana and 20 in rice (Fig. 2). According to the constructed phylogenetic relationship, the 86 TIFYs were classified into four major phylogenetic groups: JAZ, ZML, TIFY and PPD. Among which the JAZ was the largest group with 54 TIFY family members, and could be further clustered into six subgroups (JAZ I to JAZ VI) (Fig. 2). Each of the six JAZ subgroups contained different numbers of JAZ proteins from A. eriantha, A. chinensis, Arabidopsis and rice. For example, the JAZ I, JAZ III, JAZ V and JAZ VI subgroups all contained members of JAZ proteins from the four species. However, the JAZ II subgroup only contained seven JAZ proteins from rice, and no JAZ protein from other species was clustered in this subgroup (Fig. 2). In addition, no JAZ protein from A. eriantha was found in the subgroup of JAZ IV. The group ZML was the second largest group and contained three ZML proteins from Arabidopsis, four ZML proteins from rice, nine ZML proteins from A. eriantha, four ZML proteins from A. chinensis. The TIFY subgroup consisted of six TIFY members, including two TIFY proteins from each of A. eriantha and A. chinensis, and one TIFY member from each of Arabidopsis and rice. The PPD subgroup contained two A. eriantha PPD proteins, two A. chinensis PPD proteins and two Arabidopsis PPD proteins, but did not contain rice PPD proteins (Fig. 2). In addition, phylogenetic analysis was also performed using only the alignment of the 48 kiwifruit TIFY protein sequences characterized herein (Fig. 3A). The topology of the resulting phylogenetic tree indicated that TIFY proteins from the same group tended to be clustered together, which was similar to that of the above phylogenetic tree constructed by TIFY sequences from the four plants (Fig. 2, Fig. 3A).

Fig. 2
figure 2

Neighbor-joining phylogenetic relationships of the TIFY proteins from kiwifruit, Arabidopsis and Oryza sativa. The phylogenetic tree was generated using MEGA-X software and the bootstrap method (1000 replicates) with full-length amino acid sequences from A. chinensis, A. eriantha, Arabidopsis and O. sativa

Fig. 3
figure 3

Phylogenetic tree, gene structure and conserved motif of kiwifruit TIFY gene family. A The neighbor-joining phylogenetic tree of kiwifruit TIFY family. B Exon-intron structure of kiwifruit TIFY genes. Black boxes indicate exons, and black lines indicate introns. The scale bar indicates 1 kb. C Distribution of the conserved motifs in kiwifruit TIFY proteins. Ten conserved motifs are marked with different colored boxes. The scale bar indicates 100 aa

Sequence analysis of kiwifruit TIFY family

The exon-intron structure of all the identified kiwifruit TIFY genes was investigated to better understand the structural diversity of these genes. As shown in Fig. 3B, there were some differences in the number of exons and the length of introns among members of the TIFY gene family. For example, the number of exons varied from 3 to 21, and the AcZML1 had 21 exons, which was much more than that of other TIFY genes. Furthermore, the second intron of AcJAZ12 was the longest among all the sequences, which resulted in the length of the genomic DNA sequence reaching 75.656 kb. The gene exon-intron organization also showed that the structural differences of genes in the same group/subgroup were small, and they often have similar gene structures and exon/intron numbers (Fig. 3B). For example, the genes in the subgroup JAZ IV and JAZ V had five and three exons, respectively. In the JAZ III subgroup, only AcJAZ4 had six exons, while the other nine genes in this subgroup had five exons. In the PPD group, except AePPD1 which only contained three exons, all other genes contained nine exons. In the TIFY group, the number of exons varied from four to six, and most genes had six exons (Fig. 3B).

To further reveal the structural and functional characteristics of kiwifruit TIFYs, the online website MEME was employed to analyze the conserved structural motifs of the TIFY proteins. A total of 10 conserved motifs were identified, and named motif one to motif 10. As showed in Fig. 3C, no member of the kiwifruit TIFY family had a complete set of 10 conserved motifs, the motif number of the TIFY family members ranged from one to six. Members of the ZML group generally contained more motifs, and eight out of 13 ZML members contained more than five motifs, while members of other groups or subgroups contained less than four motifs. Of the 10 conserved motifs, only motif one existed in all kiwifruit TIFY members, while motif two also existed in most members (Fig. 3C). In addition, some motifs existed only in specific groups, such as motif four only existed in JAZ VI subgroup, motif nine only existed in PPD group, motif 10 only existed in TIFY group. In general, members of the same group or subgroup usually had the same type and number of motifs, and the distribution of the motifs was often the same. For example, the members of TIFY group mainly contained four motifs, the members of PPD group contained three motifs, members of JAZ I subgroup mainly contained three motifs, and most members of JAZ III subgroup mainly contained two motifs (Fig. 3C).

In addition, the conserved domains in kiwifruit TIFYs were also examined using the Pfam web server, and seven putative conserved domains were identified, namely the TIFY domain, CCT domain, GATA domain, Jas_motif domain (also named CCT_2), Dynamin_M domain, GED domain and Transp_inhibit domain as showed in Fig. 4A. Among them, the Dynamin_M and GED domains only existed in AcZML1, and Transp_inhibit domain only existed in AeJAZ13. These three domains may be sequence specific functional domains. Other conserved domains of TIFY, CCT, GATA and Jas, however, tended to exist in multiple sequences and were group/subgroup specific. For example, most members of the ZML group contained four conservative domains, including TIFY, CCT, GATA and Jas. The JAZ group (including five subgroups) mainly contained TIFY and Jas domains, and the PPD group also mainly contained the two domains of TIFY and Jas, while the TIFY group only contained TIFY domain (Fig. 4A). The conserved TIFY domain was found in all the TIFY family members, and the Jas domain was also found in most TIFY members (Fig. 4A).

Fig. 4
figure 4

Conserved domain distributions and sequences logos of kiwifruit TIFY family. A The distribution of conserved domains in kiwifruit TIFY proteins. Different conserved domains are indicated by different colored boxes. The scale bar indicates 100 aa. B Sequence logo of the TIFY domain. C Sequence log of the Jas domain

To further display the conservation patterns of the conserved domains of TIFY and Jas, their sequence logos were generated as showed in Fig. 4B and C. The TIFY domain logo revealed that the TIFY domains were not well conserved, but most of them shared common motifs, such as TIF[Y/F]XG, TLXFXG, SLSFQG (Fig. 4B). Among these conserved motifs, the motif TIF[Y/F]XG was shared by the group members of PPD, TIFY and JAZ, while the motifs of TLXFXG and SLSFQG were only shared by the ZML members (Additional file 1). The TIF[Y/F]XG was the most dominant motif, and 33 out of 48 TIFY family members contained this motif (Additional file 1). Compared with TIFY domain, the Jas domain were more conserved, and shared more conserved residues at the motif of SLX2FX2KRX2R (Fig. 4C).

Analysis of the promoter cis-elements in the TIFY gene family

In order to explore the possible expression regulation patterns of the kiwifruit TIFY genes, the cis-elements in the promoter regions of the TIFY family members were predicted. The results showed that various putative cis-elements that are involved in stress-response, phytohormone, and plant growth and development were identified widely in the promoter sequence of each kiwifruit TIFY gene (Additional file 2). The promoter region of the TIFY gene family mainly contained six types of stress-response related cis-elements, including MYB binding site involved in drought-inducibility (MBS), low-temperature responsiveness element (LTR), anaerobic induction element (ARE), defense and responsiveness element (TC-rich repeats), wound-responsive element (WUN-motif) and light responsiveness element (G-box). Among them, the cis-elements ARE and G-box were the most distributed in the kiwifruit TIFY families. In the AeTIFY family, 26 members each contained these two cis-elements, and in the AcTIFY family, 19 members each contained these two cis-elements. The phytohormone related cis-acting elements primarily involved in the gibberellin-responsiveness (P-box, GARE-motif, TATC-box), MeJA-responsiveness (CGTCA-motif, TGACG-motif), auxin responsiveness (AuxRR-core, TGA-element), abscisic acid responsiveness (ABRE), salicylic acid responsiveness (TCA-element). Although three types of gibberellin-related cis-elements were identified in the TIFY promoter regions of kiwifruit, the number of TIFY gene family members containing gibberellin related cis-elements was less than that of MeJA, ABA and salicylic acid related cis-elements. As showed in Additional file 2, there were 17, 21, 23 and 16 members in the AeTIFY family, and 17, 17, 20 and 10 members in the AcTIFY family, containing MeJA (CGTCA motif), MeJA (TGACG motif), ABA and salicylic acid related cis-elements, respectively. In addition, the TIFY gene promoters also contained some cis-acting elements related to plant growth and development, including zein metabolism regulation (O2-site), meristem expression (CAT-box), meristem specific activation (CCGTCC-box), circadian control (circadian), endosperm expression (GCN4-motif) and seed-specific regulation (RY-element). Different numbers of regulatory elements were identified in the promoter regions of the TIFY genes, such as six to 14 regulatory elements were found in the promoter region of the AeTIFY members, and seven to 12 were found in AcTIFY members (Additional file 2).

Chromosomal distribution and gene duplication analyses of the TIFY genes

The results of chromosomal distribution analyses showed that the TIFY genes were distributed irregularly on the kiwifruit chromosomes (Fig. 5). The 27 genes of the AeTIFY family members were unevenly distributed on 19 of the 29 chromosomes of A. eriantha (Fig. 5 A). Among them, Chr06 had the most genes distributed with four AeTIFY members, followed by Chr29 with three AeTIFY members. Chr11, Chr20 and Chr23 with two AeTIFY members each, and only one AeTIFY member on each of the remaining 14 chromosomes (Fig. 5A). Of the 21 TIFY genes identified from the A. chinensis genome, 20 genes were unevenly distributed on 14 out of the 29 linkage groups (LGs) of A. chinensis, and one gene (AcJAZ13) was situated on Contig00986 (Additional file 1, Fig. 5B). Among the 14 LGs with TIFY gene distribution, LG6 contained four TIFYs and had the largest number of TIFY genes. LG7, LG17 and LG20 each had two TIFYs, and the remaining ten LGs (LG9, LG11, LG13, LG17, LG18, LG21, LG22, LG25, LG28 and LG29) each had one TIFY gene.

Fig. 5
figure 5

Chromosomal distribution and interchromosomal relationships of kiwifruit TIFY genes. A Chromosomes of A. eriantha; B Chromosomes of A. chinensis (Hongyang) v3.0. Gray lines indicate all synteny blocks in kiwifruit genomes, and the red lines indicate duplicated TIFY gene pairs

Several gene duplication events were detected throughout the TIFY genes in kiwifruit. The gene pair of AeZML8/AeZML9 in AeTIFY family was detected as tandem duplication genes according to previous descriptions of tandem duplication event [40, 41]. However, no tandem duplication event was detected in AcTIFY family (Fig. 5, Additional file 3). In addition to the tandem duplication event, 11 and 17 segmental duplication events were also identified with MCScanx method in AeTIFY family and AcTIFY-v3 family, respectively (Fig. 5, Additional file 3). Among the identified gene duplication events, most of them were segmental duplication events.

The Ka/Ks ratios were calculated to investigate potential selective pressure of the identified duplication gene pairs. The results showed that all the Ka/Ks values of the above detected tandemly and segmentally duplicated TIFY gene pairs were less than one (Additional file 3), suggesting that the repetitive TIFY genes in kiwifruit were primarily constrained by intense purification selection pressure.

Synteny and evolutionary analyses of kiwifruit TIFY genes and other plants TIFYs

To explore the potential evolutionary clues of the kiwifruit TIFY gene family, a serious of comparative syntenic graphs of kiwifruit associated with other five representative plant species were constructed. The five representative plant species contained four dicots plants (Arabidopsis thaliana, Camellia sinensis, Solanum lycopersicum and Vitis vinifera) and a monocots plant (Oryza sativa).

According to the comparative syntenic maps of A. eriantha associated with the five plant species (Fig. 6), a total of 19 AeTIFY genes showed syntenic relationships with those in Vitis vinifera, followed by Solanum lycopersicum (16), Camellia sinensis (13), Arabidopsis thaliana (12) and Oryza sativa (3). The numbers of AeTIFYs orthologous genes in Vitis vinifera, Solanum lycopersicum, Camellia sinensis, Arabidopsis thaliana and Oryza sativa were 25, 29, 19, 15 and 7, respectively (Additional file 4). Among the orthologous gene pairs, two AeTIFY genes of DTZ79_11g10370 in the syntenic analysis A.eriantha and Solanum lycopersicum and DTZ79_20g14030 in the syntenic analysis of A.eriantha and Oryza sativa were identified to be associated with at least three syntenic gene pairs. Although more syntenic gene pairs were identified between A. eriantha and dicots than those between A. eriantha and O. sativa, three collinear gene pairs (DTZ79_11g10370, DTZ79_13g04380, DTZ79_20g14030) were found both in dicots and monocots plants. However, some collinear gene pairs (DTZ79_04g05010, DTZ79_21g05530, DTZ79_28g03670) were just available in dicots plants, but not identified in the monocot plant of Oryza sativa (Additional file 4).

Fig. 6
figure 6

Synteny analysis of TIFY genes between A. eriantha and five representative plant species. Gray lines in the background indicate the collinear blocks within A. eriantha and other plant genomes, while the red lines highlight the syntenic TIFY gene pairs

As showed in the comparative syntenic maps of A. chinensis (Fig. 7), totally 19 AcTIFY family members showed a syntenic relationship with Vitis vinifera, followed by Camellia sinensis (17), Solanum lycopersicum (16), Arabidopsis thaliana (14) and Oryza sativa (7). The numbers of AcTIFY genes in Vitis vinifera, Camellia sinensis, Solanum lycopersicum, Arabidopsis thaliana and Oryza sativa were 24, 25, 24, 19 and 12, respectively (Additional file 5). In the identified orthologous gene pairs, three AcTIFY genes (Actinidia32203.t1 in A. chinensis and Solanum lycopersicum, Actinidia30242.t1 and Actinidia32203.t1 in A. chinensis and Oryza sativa) were also found to be associated with at least three corresponding gene pairs. For the syntenic analysis of the TIFY genes of AcTIFY and the five representative plants, five syntenic gene pairs (Actinidia09609.t1, Actinidia22128.t1, Actinidia30242.t1, Actinidia30601.t1, Actinidia32203.t1) were found both in dicots and monocots plants. However, another six syntenic gene pairs (Actinidia09108.t1, Actinidia10272.t1, Actinidia12384.t1, Actinidia13280.t1, Actinidia18522.t1 and Actinidia19959.t1) were only found in the four dicots plants, but not available in the monocot plant of O. sativa (Additional file 5).

Fig. 7
figure 7

Synteny analysis of TIFY genes between A. chinensis and five representative plant species. Gray lines in the background indicate the collinear blocks within A. chinensis (Hongyang) v3.0 and other plant genomes, while the red lines highlight the syntenic TIFY gene pairs

In addition, to further display the evolutionary relationships of the kiwifruit TIFY gene family, a multicollinearity plot of the TIFY genes among the two genomes of A. eriantha and A. chinensis was drawn (Fig. 8). The results showed that a total of 18 AeTIFY genes and 19 AcTIFY genes had a collinearity relationship, and 39 collinear gene pairs were identified between AeTIFY gene family and AcTIFY gene family (Fig. 8, Additional file 6).

Fig. 8
figure 8

Synteny analysis of TIFY genes between different kiwifruit taxa. Gray lines in the background represent the collinear blocks between A. eriantha and A. chinensis (Hongyang) v3.0 genomes, while the red lines highlight the syntenic TIFY gene pairs

Furthermore, the Ka/Ks ratios of all the above identified collinear gene pairs were also calculated, and the results showed that all collinear gene pairs displayed Ka/Ks < 1, indicating that the evolution of kiwifruit TIFY gene family might have suffered strong purifying selective pressure (Additional files 4, 5, 6).

Expression pattern analyses of kiwifruit TIFY genes

In this study, the expression patterns of kiwifruit TIFY genes in different tissue parts, different fruit development stages and Psa invasion stages were investigated. The expression patterns of the A. chinensis TIFY family genes in three main tissues (root, stem and leaf) were invested based on previous RNA-seq data. As showed in Fig. 9A, most of the TIFY genes were expressed at different levels in these three tissues, and the expression of some genes is highly tissue-specific. Interestingly, compared with other family members, AcZML1, AcZML4, AcTIFY2, AcJAZ11 and AcJAZ6 had relatively higher expression levels. Among them, AcJAZ11 had a high expression level in root, stem and leaf, and the expression level in leaf and stem was particularly prominent, indicating that this gene may be necessary for the development of tissues such as leaf and stem. AcJAZ6 was specifically highly expressed in leaf, suggesting that this gene may play an important role in leaf. Similar to AcJAZ6, the expression level of AcTIFY2 in leaf was also higher than that in root and stem. AcZML1 and AcZML4 were expressed in all expressed in these three tissues, and no obvious tissue-specific expression pattern was observed (Fig. 9A).

Fig. 9
figure 9

Expression profiles of kiwifruit TIFY genes in different tissues, different stages of fruit development and Psa invasion. A Expression profiles of AcTIFY genes in different tissues and different stages of fruit development. Leaf, root and stem indicate different tissues in A. chinensis. DAP20_immature, DAP120_mature green and DAP127_ripe represent different stages of fruit development after pollination. B Expression profiles of AcTIFY genes in response to Psa infection and ASM (acibenzolar-S-methy) treatment. HealtyControl represents samples without Psa inoculation or ASM treatment. Psa3, Psa24 and Psa48 represent hours post inoculation of Psa. C Expression profiles of AcTIFY genes in the process of Psa invasion. AH_0DPI, AH_2DPI, AH_14DPI represent days post inoculation with Psa in A. chinensis. D Expression profiles of AeTIFY genes in the process of Psa invasion. Ae_0DPI, Ae_2DPI, Ae_14DPI indicate days post inoculation with Psa in A. eriantha. The expression values of kiwifruit TIFYs were normalized to FPKM (fragments per kilobase of exon per million mapped fragments)

The expression patterns of A. chinensis TIFY gene members during fruit development and ripening (DAP20_immature, DAP120_mature green and DAP127_ripe) were further analyzed according to previous transcriptome data. As the results showed in Fig. 9A, the expression levels of AcZML2 and AcJAZ9 increased gradually along with fruit ripening, and the change trend of AcJAZ9 gene expression was more obvious. However, the expression levels of AcTIFY1 and AcTIFY2 decreased gradually with fruit ripening (Fig. 9A). In addition, there were also some genes that had higher expression levels at specific stages of fruit development. For example, the expression of AcJAZ11 was higher in immature and ripe stages, but lower in mature green stage. AcJAZ3 was highly expressed in ripe stage, while AcJAZ9 was highly expressed in mature green stage and ripe stage (Fig. 9A). The diversity of AcTIFYs expression patterns during fruit development indicated that these genes may played different roles along with fruit ripening.

The expression patterns of TIFY family members in A. chinensis and A. eriantha after the invasion of Psa were detected. The expression profile of AcTIFYs with or without ASM treatment in the process of Psa infection was shown in Fig. 9B. The expression levels of AcTIFY genes changed in varying degrees after Psa invasion. AcJAZ12, AcJAZ9, AcJAZ4 and AcJAZ5 had similar expression patterns. These genes all had high expression levels in the early stage of Psa infection, and the expression of these genes gradually decreased with the increase of infection. AcJAZ11 maintained a high expression level during the whole process of Psa infection. The expression of AcJAZ11 was decreased after ASM treatment compared with the samples without ASM treatment. On the contrary, the expression of AcJAZ9 after ASM treatment was higher than that without ASM treatment (Fig. 9B). Figure 9C was also the expression profile of AcTIFY genes to the invasion of Psa. Similar to the expression pattern of AcJAZ11 in Fig. 9B, AcJAZ11 in Fig. 9C also had the highest expression level in the early stage of Psa infection (AH_0DPI), and the expression level of this gene gradually decreased with the increase of infection time (Fig. 9C). Figure 9D was the expression profile of AeTIFY genes to the invasion of Psa. The expression levels of AeJAZ10 and AeJAZ11 genes remained at a high level at different periods after infection, while the expression level of AeJAZ11 was the highest in different periods (Fig. 9C), which also suggested AeJAZ11 might play a certain role in Psa infection.

Discussion

The TIFY gene family, as unique transcription factor family in plants, plays pivotal roles in regulating plant development, physiological processes and response to stresses. In the current study, the TIFY gene families in kiwifruit, a fruit crop with important economic values and popular worldwide, were systematically identified and analyzed.

The TIFY genes have been identified and characterized in several important horticultural crops, such as grape [28], apple [29], tomato [25, 30], pear [31], watermelon [32], tea [33]. The genus Actinidia contains 54 species and 74 taxa, but only a few species with great utilization values have been domesticated, such as A. chinensis, A. deliciosa, A. arguta and A. eriantha. At present, the commercial species widely cultivated in the world are A. chinensis and A. deliciosa, and A. eriantha has only a small amount of artificial cultivation. A. eriantha is easy to peel, enrich in vitamin C, and has strong resistance to Psa [39]. In this study, the genome-wide identification of the TIFY genes of A. chinensis and A. eriantha were performed. The results showed that there were at least 21 TIFY genes in the genome of A. chinensis, and at least 27 TIFY genes in A. eriantha (Fig. 1). The difference in the number of TIFY family genes between the two kiwifruit species was mainly due to the difference in the number of JAZ and TIFY group members. The JAZ and ZML groups in A. chinensis contained 13 and 4 members, respectively, while JAZ and ZML groups in A. eriantha contained 14 and 9 members, respectively (Fig. 1). The difference in the number of gene family members may be due to gene duplication or loss in the process of gene evolution. Gene duplication and loss were the main evolutionary driving forces for the expansion or contraction, and duplicated genes could lead to gene redundancy [42]. Duplication events in the critical sites such as the CDS, and the promoter sequence cause members of a gene family to receive new functions [43, 44].

The TIFY genes in kiwifruit genomes showed a high variation in their sequence structure. In terms of protein length, the variation range of amino acid sequence of AcTIFYs was 107 aa to 764 aa, while that of AeTIFYs was 112 aa to 718 aa (Additional file 1). In terms of gene structure, the variation range of kiwifruit TIFY gene exon was three to 21 (Fig. 3B). The variation range of the conserved motif of kiwifruit TIFY gene was one to six, and the variation range of the conserved domain of kiwifruit TIFY was also one to six. Some conserved motifs or domains were unique to specific sequences or subgroups (Figs. 3 C and 4A). Those high variation in the sequence structure revealed that TIFY family members have acquired changes in their genome during evolution events that affected their functions [43, 45].

Segmental and tandem gene duplications are the two major factors in the generation and maintenance of gene family, and the relative importance of segmental and tandem duplication in the evolution of gene family may correspond to functional differences of the gene family members [46]. In this study, only two tandem duplication events were detected in the kiwifruit TIFY gene families, while the rest were segmental duplications. This typical type of low tandem and high segmental duplications is consistent with the classification of gene duplication types in previous works, which mainly includes proteins involved in the roles of transcription factors, signaling, membrane transport and so on [46]. Previous studies have shown that the functions and expression patterns of segmentally duplicated genes were often similar [47, 48]. In this study, the expression pattern of segmentally duplicated gene pairs was not completely consistent, which was different from previous studies. For the duplicated gene pairs of AcJAZ9/AcJAZ4, the expression of AcJAZ9 was much higher than AcJAZ4 during fruit development. During Psa infection, the expression pattern of AcJAZ9 was also different from that of AcJAZ4. In A. eriantha, the expression patterns of the duplicated gene pairs AeJAZ10/AeJAZ11 were different, but the duplicated gene pairs AeZML1/AeZML6, AeZML3/AeZML4 had similar expression patterns. The different expression patterns between duplicated gene pairs indicated that the gene pairs may perform different functions [43]. In addition, the strong purifying selection signals detected in the duplication gene pairs also indicated the functional importance of kiwifruit TIFY genes. These results indicated that gene duplication events, especially segmental duplication, contributed to the evolution and expansion of the TIFY gene family in kiwifruit species.

The yield and quality of kiwifruit are easily affected by a variety of biotic and abiotic stresses. Abiotic stresses such as salt, temperature and waterlogging often have adverse effects on plant growth and development. Kiwifruit is a salt sensitive plant, and salt stress seriously affects the normal growth and physiological processes of kiwifruit plants [49]. Temperature influences shoot growth and maturation of fruit on kiwifruit. Under high temperature, the contents of carbohydrate and vitamin C in fruit reduced significantly, and the ripening rate and ethylene biosynthesis also decreased, which finally affected the growth and maturation of fruit [50, 51]. In addition, the biosynthesis and transportation of anthocyanins in red-fleshed kiwifruit were also inhibited by high temperature, which affected the nutrition and commercialization of the fruit [52]. Kiwifruit plants are very sensitive to waterlogging stress, which can lead to fruit yield reduction and even plant death in severe cases [53]. In this study, multiple putative cis-elements, which are mainly involved in stress-response, phytohormone, and plant growth and development, were identified from the promoter region of kiwifruit TIFY family genes. The diverse cis-regulatory elements in the promoter region of kiwifruit TIFY genes indicated that these genes may be involved in response to a variety of stresses and multiple plant hormones response processes, and played a certain role in the growth and development of kiwifruit.

Kiwifruit trees may encounter a variety of biotic stresses during growth, especially the threats of pests and diseases during the growth process. At present, the main threat to the development of kiwifruit industry is kiwifruit bacterial canker. The bacterial canker disease is a devastating disease incited by Psa, which seriously threatens the production and development of kiwifruit industry worldwide [54]. The disease has the characteristics of rapid transmission and strong pathogenicity, and can lead to the death of kiwifruit trees in a large area [55]. However, there is no effective method to prevent and treat the bacterial canker at present. The resistance of kiwifruit varieties and species to bacterial canker disease is very different. For example, A. chinensis ‘Hongyang’ and ‘Hort16A’ are highly susceptible to canker disease, while A. deliciosa ‘Jinkui’ and A. eriantha are highly resistant to canker disease [55, 56]. Therefore, breeding resistant varieties and exploring resistance genes and biological pathways may be an effective way to control the disease. This study analyzed the expression patterns TIFY genes in kiwifruit after Psa invasion. The gene expression of A. chinensis ‘Hongyang’, which is susceptible to Psa infection, showed that AcJAZ11 had a higher expression level during the process of Psa infection, but its expression level decreased with the increase of the infection time. A. eriantha has strong resistance to Psa, and the gene expression results showed that AeJAZ10 and AeJAZ11 had maintained high levels during the whole process of Psa infection. It was worth noting that the expression level of AeJAZ10 in the initial stage of infection is higher than that in other stages of infection, while the expression level of AeJAZ11 in the initial stage (Ae_0DPI) is lower than that of the Ae_2DPI stage of infection, indicating that the expression level of AeJAZ11 has a tendency to increase with the increase of infection time. These results indicated that AeJAZ11 may be a candidate gene in response to Psa infection and play a certain role in the process of Psa invasion. Previous studies have shown that many JAZ genes are key regulatory factors in the jasmonic acid (JA) signal pathway, and participate in response to the infection process of a variety of plant pathogens by mediating JA signal transduction [33, 34, 57,58,59,60,61]. The specific mechanism of the candidate gene AeJAZ11 in Psa infection, and whether JA pathway plays a role in the process of Psa infection still needs to be further verified in future studies.

Conclusions

The kiwifruit TIFY family genes were comprehensively and systematically characterized in the present study. A total of 27 and 21 TIFY family genes were genome-widely identified in the genomes of A. eriantha and A. chinensis, respectively. The phylogenetic analysis showed that the identified kiwifruit TIFY genes could be classified into four main groups of JAZ, ZML, TIFY and PPD, and members within the same group had similar gene structures and motif compositions. Collinearity analysis suggested that segmental duplication events played a major role in the expansion of kiwifruit TIFY family genes. The molecular evolution analysis indicated that the evolution of kiwifruit TIFY genes were dominated by purifying selection. Promoter cis-elements analysis, spatio-temporal expression pattern analysis and expression characteristics analysis of Psa invasion showed that a few TIFY genes might be involved in the growth and development of kiwifruit and the response to Psa invasion stress. The results presented in this study laid a foundation for further exploring and understanding the biological functions of TIFY family genes in kiwifruit.

Methods

Identification of TIFY family genes in kiwifruit genomes

In order to identify the TIFY family members from kiwifruit genomes, the whole-genome and proteome data of A. chinensis (Hong Yang) v3.0 and A. eriantha (White) were downloaded from the Kiwifruit Genome Database (KGD: http://kiwifruitgenome.org/) [62]. The hidden Markov model (HMM) profile of the conservative functional domain of TIFY (PF06200) was obtained from the Pfam database v34.0 (http://pfam.xfam.org/), and the HMM profile was further used to screen the kiwifruit proteomes using the hmmsearch software in the HMMER package v3.0 to obtain the potential gene family members of TIFY. After removing redundant and incomplete sequences, the conserved domain architectures of the acquired sequences were further confirmed by Pfam database and SMART website (http://smart.embl-heidelberg.de/). Sequences without the typical functional domain of TIFY were excluded from the dataset, and the amino acid sequences containing the conservative TIFY domain were regarded as potential members of the kiwifruit TIFY family and would be used for subsequent analysis. In addition, the TIFY protein sequences of model plants Arabidopsis thaliana and rice (Oryza sativa) were downloaded from TAIR (https://www.arabidopsis.org/) and TIGR (http://rice.plantbiology.msu.edu/) database, respectively.

The basic physicochemical properties, such as the molecular weight (MW) and isoelectric point (pI) of each kiwifruit TIFY protein, were predicted using the tool of ProtProm in ExPASy (https://web.expasy.org/protparam/). The subcellular localizations of kiwifruit TIFYs were predicted by the web-server of Cell-PLoc 2.0 (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/).

Sequence alignment and phylogenetic analysis

The multiple sequence alignment was performed using the program of MAFFT v7.409 [63] with auto strategy and default parameter settings. Then the resultant multiple alignment file was imported into MEGA-X software [64] to construct phylogenetic relationship using the Neighbor-Joining (NJ) algorithm, and the tree topology support was assessed by bootstrap analysis with 1000 replicates. Finally, the constructed phylogenetic tree was annotated and visualized using the online tool of EvolView (https://www.evolgenius.info/evolview/) [65].

Gene structure analysis and conserved motif discovery

The coding sequences (CDS) and their corresponding genomic sequences of all kiwifruit TIFY genes were obtained from the KGD, and then submitted to the online tool of Gene Structure Display Server (GSDS, http://gsds.cbi.pku.edu.cn/) to analyze and visualize the exon-intron organization of these genes. The MEME web server v5.33 (https://meme-suite.org/meme/tools/meme) was used to discover the conserved motif patterns of the identified TIFY proteins in kiwifruit. The maximum number of motifs was set at 10, and the default values were employed for the other parameters. In addition, the conserved domain composition of kiwifruit TIFY proteins were also analyzed using Pfam database and then visualized using TBtools [66]. Furthermore, the sequence logos of the conserved TIFY and Jas functional domains were generated using the web based application of WebLogo (http://weblogo.threeplusone.com/) [67].

Promoter cis-elements analysis of kiwifruit TIFY genes

To explore the putative cis-regulatory elements in the promoter regions of kiwifruit TIFY genes, the upstream 2000 bp sequences of the transcription initiation codon of all the TIFY genes were extracted from the three kiwifruit genomes, and then the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was used to predict the cis-elements in the promoter regions of kiwifruit TIFY genes.

Chromosomal locations, gene duplication and synteny analysis

To understand the chromosomal distributions of TIFY genes in the three kiwifruit genomes, the chromosomal position information of the TIFY genes were obtained from the genome sequence files and the corresponding gene structure annotation information files, which were downloaded from the KGD database. Then, the One Step MCSanX function in TBtools software was adopted to analyze the duplication types and intraspecific collinearity of TIFY family members. Finally, the Advanced Circles function in TBtools was used to draw the chromosomal location map and the collinearity relationship of the TIFY genes. In addition, the Multiple Synteny Plot function in TBtools was employed to constructed and exhibited the synteny relationship of the orthologous TIFY genes obtained from A. chinensis and A. eriantha. Furthermore, the homology of the TIFY genes between kiwifruit and the other five plants (including Arabidopsis thaliana, Camellia sinensis, Oryza sativa, Solanum lycopersicum and Vitis vinifera) were performed using the Dual Synteny Plot in TBtools. The genome sequences and general feature format files of the five selected plants were downloaded from NCBI (https://www.ncbi.nlm.nih.gov/genome). The Simple Ka/Ks Calculator (NG) of TBtools was used to calculate the nonsynonymous (Ka) and synonymous (Ks) substitution values and the Ka/Ks values of each duplicated TIFY gene pairs and the syntenic TIFY gene pairs, where Ka/Ks < 1 indicated purifying selection, Ka/Ks = 1 indicated neutral selection, and Ka/Ks > 1 indicated positive selection [68].

Expression pattern analysis of kiwifruit TIFY genes

To explore the expression patterns of kiwifruit TIFY genes, the RNA-seq expression profiles of the TIFY genes were mined from the KGD database under the project of PRJNA187369, PRJNA328414 and PRJNA436459. The PRJNA187369 project is the transcriptome analysis of A. chinensis ‘Hongyang’ leaves and fruits at different developmental stages [36]. The PRJNA328414 project is the transcriptome data of different kiwifruit taxa (A. chinensis ‘Hongyang’ and A.eriantha) infected by the pathogen Psa of kiwifruit canker disease [69, 70]. The PRJNA436459 project is the transcriptome analysis of A. chinensis with or without ASM (acibenzolar-S-methyl) treatment during the inoculation of Psa [71]. The heat maps of the expression levels of the TIFY genes were visualized using the Heatmap illustrator program in the toolkit of TBtools.

Availability of data and materials

The kiwifruit genome sequences of A. chinensis (Hong Yang) v3.0 and A. eriantha (White) were downloaded from the Kiwifruit Genome Database (KGD: http://kiwifruitgenome.org/), and the RNA-seq data were also downloaded from KGD. The TIFY protein sequences of Arabidopsis thaliana and rice were downloaded from TAIR (https://www.arabidopsis.org/) and TIGR (http://rice.plantbiology.msu.edu/) databases, respectively. The genome sequences and general feature format files of the five selected plants (Arabidopsis thaliana, Camellia sinensis, Oryza sativa, Solanum lycopersicum and Vitis vinifera) were downloaded from NCBI (https://www.ncbi.nlm.nih.gov/genome).

Abbreviations

MW:

Molecular weight

pI:

Isoelectric point

CDS:

Coding sequences

GSDS:

Gene Structure Display Server

KGD:

Kiwifruit Genome Database

Ka:

Nonsynonymous

Ks:

Synonymous

DPI:

Days post inoculation

Psa:

Pseudomonas syringae pv. actinidiae

JA:

Jasmonic acid

ASM:

Acibenzolar-S-methyl

DAP:

Days after pollination

HMM:

Hidden Markov model

References

  1. Vanholme B, Grunewald W, Bateman A, Kohchi T, Gheysen G. The tify family previously known as ZIM. Trends Plant Sci. 2007;12(6):239–44.

    Article  CAS  PubMed  Google Scholar 

  2. Nishii A, Takemura M, Fujita H, Shikata M, Yokota A, Kohchi T. Characterization of a novel gene encoding a putative single zinc-finger protein, ZIM, expressed during the reproductive phase in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2000;64(7):1402–9.

    Article  CAS  PubMed  Google Scholar 

  3. Bai Y, Meng Y, Huang D, Qi Y, Chen M. Origin and evolutionary analysis of the plant-specific TIFY transcription factor family. Genomics. 2011;98(2):128–36.

    Article  CAS  PubMed  Google Scholar 

  4. Chung HS, Niu Y, Browse J, Howe GA. Top hits in contemporary JAZ: an update on jasmonate signaling. Phytochemistry. 2009;70(13–14):1547–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Staswick PE. JAZing up jasmonate signaling. Trends Plant Sci. 2008;13(2):66–71.

    Article  CAS  PubMed  Google Scholar 

  6. Shikata M, Matsuda Y, Ando K, Nishii A, Takemura M, Yokota A, et al. Characterization of Arabidopsis ZIM, a member of a novel plant-specific GATA factor gene family. J Exp Bot. 2004;55(397):631–9.

    Article  CAS  PubMed  Google Scholar 

  7. White DW. PEAPOD regulates lamina size and curvature in Arabidopsis. Proc Natl Acad Sci U S A. 2006;103(35):13238–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Baekelandt A, Pauwels L, Wang Z, Li N, De Milde L, Natran A, et al. Arabidopsis leaf flatness is regulated by PPD2 and NINJA through repression of CYCLIN D3 genes. Plant Physiol. 2018;178(1):217–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zhu Y, Luo X, Liu X, Wu W, Cui X, He Y, et al. Arabidopsis PEAPODs function with LIKE HETEROCHROMATIN PROTEIN1 to regulate lateral organ growth. J Integr Plant Biol. 2020;62(6):812–31.

    Article  CAS  PubMed  Google Scholar 

  10. Hakata M, Kuroda M, Ohsumi A, Hirose T, Nakamura H, Muramatsu M, et al. Overexpression of a rice TIFY gene increases grain size through enhanced accumulation of carbohydrates in the stem. Biosci Biotechnol Biochem. 2012;76(11):2129–34.

    Article  CAS  PubMed  Google Scholar 

  11. Cai Q, Yuan Z, Chen M, Yin C, Luo Z, Zhao X, et al. Jasmonic acid regulates spikelet development in rice. Nat Commun. 2014;5:3476.

    Article  PubMed  Google Scholar 

  12. Zhai Q, Zhang X, Wu F, Feng H, Deng L, Xu L, et al. Transcriptional mechanism of jasmonate receptor COI1-mediated delay of flowering time in Arabidopsis. Plant Cell. 2015;27(10):2814–28.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Guan Y, Ding L, Jiang J, Shentu Y, Zhao W, Zhao K, et al. Overexpression of the CmJAZ1-like gene delays flowering in Chrysanthemum morifolium. Hortic Res. 2021;8(1):87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yu X, Chen G, Tang B, Zhang J, Zhou S, Hu Z. The Jasmonate ZIM-domain protein gene SlJAZ2 regulates plant morphology and accelerates flower initiation in Solanum lycopersicum plants. Plant Sci. 2018;267:65–73.

    Article  CAS  PubMed  Google Scholar 

  15. Meng L, Zhang T, Geng S, Scott PB, Li H, Chen S. Comparative proteomics and metabolomics of JAZ7-mediated drought tolerance in Arabidopsis. J Proteomics. 2019;196:81–91.

    Article  CAS  PubMed  Google Scholar 

  16. Seo JS, Joo J, Kim MJ, Kim YK, Nahm BH, Song SI, et al. OsbHLH148, a basic helix-loop-helix protein, interacts with OsJAZ proteins in a jasmonate signaling pathway leading to drought tolerance in rice. Plant J. 2011;65(6):907–21.

    Article  CAS  PubMed  Google Scholar 

  17. Ye H, Du H, Tang N, Li X, Xiong L. Identification and expression profiling analysis of TIFY family genes involved in stress and phytohormone responses in rice. Plant Mol Biol. 2009;71(3):291–305.

    Article  CAS  PubMed  Google Scholar 

  18. Zhao C, Pan X, Yu Y, Zhu Y, Kong F, Sun X, et al. Overexpression of a TIFY family gene, GsJAZ2, exhibits enhanced tolerance to alkaline stress in soybean. Mol Breeding. 2020;40(3):33.

    Article  CAS  Google Scholar 

  19. Ebel C, BenFeki A, Hanin M, Solano R, Chini A. Characterization of wheat (Triticum aestivum) TIFY family and role of Triticum Durum TdTIFY11a in salt stress tolerance. PLoS One. 2018;13(7):e0200566.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Chung HS, Koo AJ, Gao X, Jayanty S, Thines B, Jones AD, et al. Regulation and function of Arabidopsis JASMONATE ZIM-domain genes in response to wounding and herbivory. Plant Physiol. 2008;146(3):952–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Jing Y, Liu J, Liu P, Ming D, Sun J. Overexpression of TaJAZ1 increases powdery mildew resistance through promoting reactive oxygen species accumulation in bread wheat. Sci Rep. 2019;9(1):5691.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Li Y, Xu M, Wang N, Li Y. A JAZ protein in Astragalus sinicus interacts with a leghemoglobin through the TIFY domain and is involved in nodule development and nitrogen fixation. PLoS One. 2015;10(10):e0139964.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Sen S, DasGupta M. Involvement of Arachis hypogaea Jasmonate ZIM domain/TIFY proteins in root nodule symbiosis. J Plant Res. 2021;134(2):307–26.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang Z, Li X, Yu R, Han M, Wu Z. Isolation, structural analysis, and expression characteristics of the maize TIFY gene family. Mol Genet Genomics. 2015;290(5):1849–58.

    Article  CAS  PubMed  Google Scholar 

  25. Heidari P, Faraji S, Ahmadizadeh M, Ahmar S, Mora-Poblete F. New insights into structure and function of TIFY genes in Zea mays and Solanum lycopersicum: a genome-wide comprehensive analysis. Front Genet. 2021;12:657970.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Singh P, Mukhopadhyay K. Comprehensive molecular dissection of TIFY transcription factors reveal their dynamic responses to biotic and abiotic stress in wheat (Triticum aestivum L.). Sci Rep. 2021;11(1):9739.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xie S, Cui L, Lei X, Yang G, Li J, Nie X, et al. The TIFY gene family in wheat and its progenitors: genome-wide identification, evolution and expression analysis. Curr Genomics. 2019;20(5):371–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang Y, Gao M, Singer SD, Fei Z, Wang H, Wang X. Genome-wide identification and analysis of the TIFY gene family in grape. PLoS One. 2012;7(9):e44465.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Li X, Yin X, Wang H, Li J, Guo C, Gao H, et al. Genome-wide identification and analysis of the apple (Malus × domestica Borkh.) TIFY gene family. Tree Genet Genomes. 2014;11(1):808.

    Article  Google Scholar 

  30. Chini A, Ben-Romdhane W, Hassairi A, Aboul-Soud MAM. Identification of TIFY/JAZ family genes in Solanum lycopersicum and their regulation in response to abiotic stresses. PLoS One. 2017;12(6):e0177381.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Ma Y, Shu S, Bai S, Tao R, Qian M, Teng Y. Genome-wide survey and analysis of the TIFY gene family and its potential role in anthocyanin synthesis in Chinese sand pear (Pyrus pyrifolia). Tree Genet Genomes. 2018;14(2):25.

    Article  Google Scholar 

  32. Yang Y, Ahammed GJ, Wan C, Liu H, Chen R, Zhou Y. Comprehensive analysis of TIFY transcription factors and their expression profiles under jasmonic acid and abiotic stresses in watermelon. Int J Genomics. 2019;2019:6813086.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Zhang X, Ran W, Zhang J, Ye M, Lin S, Li X, et al. Genome-wide identification of the Tify gene family and their expression profiles in response to biotic and abiotic stresses in tea plants (Camellia sinensis). Int J Mol Sci. 2020;21(21):8316.

    Article  CAS  PubMed Central  Google Scholar 

  34. He X, Kang Y, Li W, Liu W, Xie P, Liao L, et al. Genome-wide identification and functional analysis of the TIFY gene family in the response to multiple stresses in Brassica napus L. BMC Genomics. 2020;21(1):736.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Li L, Liu Y, Huang Y, Li B, Ma W, Wang D, et al. Genome-wide identification of the TIFY family in Salvia miltiorrhiza reveals that SmJAZ3 interacts with SmWD40-170, a relevant protein that modulates secondary metabolism and development. Front Plant Sci. 2021;12:630424.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Huang S, Ding J, Deng D, Tang W, Sun H, Liu D, et al. Draft genome of the kiwifruit Actinidia chinensis. Nat Commun. 2013;4:2640.

    Article  PubMed  Google Scholar 

  37. Wu H, Ma T, Kang M, Ai F, Zhang J, Dong G, et al. A high-quality Actinidia chinensis (kiwifruit) genome. Hortic Res. 2019;6:117.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Pilkington SM, Crowhurst R, Hilario E, Nardozza S, Fraser L, Peng Y, et al. A manually annotated Actinidia chinensis var. chinensis (kiwifruit) genome highlights the challenges associated with draft genomes and gene prediction in plants. BMC Genomics. 2018;19(1):257.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Tang W, Sun X, Yue J, Tang X, Jiao C, Yang Y, et al. Chromosome-scale genome assembly of kiwifruit Actinidia eriantha with single-molecule sequencing and chromatin interaction mapping. Gigascience. 2019;8(4):giz027.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet. 2001;2(7):516–27.

    Article  CAS  PubMed  Google Scholar 

  41. Xie T, Chen C, Li C, Liu J, Liu C, He Y. Genome-wide investigation of WRKY gene family in pineapple: evolution and expression profiles during development and stress. BMC Genomics. 2018;19(1):490.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Tang C, Zhu X, Qiao X, Gao H, Li Q, Wang P, et al. Characterization of the pectin methyl-esterase gene family and its function in controlling pollen tube growth in pear (Pyrus bretschneideri). Genomics. 2020;112(3):2467–77.

    Article  CAS  PubMed  Google Scholar 

  43. Abdullah FS, Mehmood F, Malik HMT, Ahmed I, Heidari P, Poczai P. The GASA gene family in cacao (Theobroma cacao, Malvaceae): genome wide identification and expression analysis. Agronomy. 2021;11(7):1425.

    Article  CAS  Google Scholar 

  44. Musavizadeh Z, Najafi-Zarrini H, Kazemitabar SK, Hashemi SH, Faraji S, Barcaccia G, et al. Genome-wide analysis of potassium channel genes in rice: expression of the OsAKT and OsKAT genes under salt stress. Genes (Basel). 2021;12(5):784.

    Article  CAS  Google Scholar 

  45. Heidari P, Abdullah FS, Poczai P. Magnesium transporter gene family: genome-wide identification and characterization in Theobroma cacao, Corchorus capsularis, and Gossypium hirsutum of family Malvaceae. Agronomy. 2021;11(8):1651.

    Article  CAS  Google Scholar 

  46. Cannon SB, Mitra A, Baumgarten A, Young ND, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4:10.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Ahmad MZ, Sana A, Jamil A, Nasir JA, Ahmed S, Hameed MU, et al. A genome-wide approach to the comprehensive analysis of GASA gene family in Glycine max. Plant Mol Biol. 2019;100(6):607–20.

    Article  CAS  PubMed  Google Scholar 

  48. Faraji S, Filiz E, Kazemitabar SK, Vannozzi A, Palumbo F, Barcaccia G, et al. The AP2/ERF gene family in Triticum durum: genome-wide identification and expression analysis under drought and salinity stresses. Genes (Basel). 2020;11(12):1464.

    Article  CAS  Google Scholar 

  49. Abid M, Zhang Y, Li Z, Bai D, Zhong Y, Fang J. Effect of salt stress on growth, physiological and biochemical characters of four kiwifruit genotypes. Sci Hortic. 2020;271:109473.

    Article  CAS  Google Scholar 

  50. Antunes MDC, Sfakiotakis EM. Effect of high temperature stress on ethylene biosynthesis, respiration and ripening of ‘Hayward’ kiwifruit. Postharvest Biol Technol. 2000;20(3):251–9.

    Article  CAS  Google Scholar 

  51. Richardson AC, Marsh KB, Boldingh HL, Pickering AH, Bulley SM, Frearson NJ, et al. High growing temperatures reduce fruit carbohydrate and vitamin C in kiwifruit. Plant Cell Environ. 2004;27(4):423–35.

    Article  CAS  Google Scholar 

  52. Man YP, Wang YC, Li ZZ, Jiang ZW, Yang HL, Gong JJ, et al. High-temperature inhibition of biosynthesis and transportation of anthocyanins results in the poor red coloration in red-fleshed Actinidia chinensis. Physiol Plant. 2015;153(4):565–83.

    Article  CAS  PubMed  Google Scholar 

  53. Zhang J-Y, Huang S-N, Mo Z-H, Xuan J-P, Jia X-D, Wang G, et al. De novo transcriptome sequencing and comparative analysis of differentially expressed genes in kiwifruit under waterlogging stress. Mol Breeding. 2015;35(11):208.

    Article  Google Scholar 

  54. Tahir J, Hoyte S, Bassett H, Brendolise C, Chatterjee A, Templeton K, et al. Multiple quantitative trait loci contribute to resistance to bacterial canker incited by Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis). Hortic Res. 2019;6(1):101.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Song Y, Sun L, Lin M, Chen J, Qi X, Hu C, et al. Comparative transcriptome analysis of resistant and susceptible kiwifruits in response to Pseudomonas syringae pv. Actinidiae during early infection. PLoS One. 2019;14(2):e0211913.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Pei Y, Ma L, Sui L, Cui Y, Liu X, Gong G. Resistance evaluation of kiwifruit cultivars to Pseudomonas syringae pv. actinidiae and utilization. Journal of fruit. Science. 2021;38(7):1153–62.

    Google Scholar 

  57. Chini A, Fonseca S, Fernandez G, Adie B, Chico JM, Lorenzo O, et al. The JAZ family of repressors is the missing link in jasmonate signalling. Nature. 2007;448(7154):666–71.

    Article  CAS  PubMed  Google Scholar 

  58. Thines B, Katsir L, Melotto M, Niu Y, Mandaokar A, Liu G, et al. JAZ repressor proteins are targets of the SCFCOI1 complex during jasmonate signalling. Nature. 2007;448(7154):661–5.

    Article  CAS  PubMed  Google Scholar 

  59. Thatcher LF, Cevik V, Grant M, Zhai B, Jones JD, Manners JM, et al. Characterization of a JAZ7 activation-tagged Arabidopsis mutant with increased susceptibility to the fungal pathogen Fusarium oxysporum. J Exp Bot. 2016;67(8):2367–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Sun Q, Wang G, Zhang X, Zhang X, Qiao P, Long L, et al. Genome-wide identification of the TIFY gene family in three cultivated Gossypium species and the expression of JAZ genes. Sci Rep. 2017;7:42418.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Liu X, Zhao C, Yang L, Zhang Y, Wang Y, Fang Z, et al. Genome-wide identification, expression profile of the TIFY gene family in Brassica oleracea var. capitata, and their divergent response to various pathogen infections and phytohormone treatments. Genes (Basel). 2020;11(2):127.

    Article  CAS  Google Scholar 

  62. Yue J, Liu J, Tang W, Wu YQ, Tang X, Li W, et al. Kiwifruit genome database (KGD): a comprehensive resource for kiwifruit genomics. Hortic Res. 2020;7:117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Subramanian B, Gao S, Lercher MJ, Hu S, Chen WH. Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res. 2019;47(W1):W270–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202.

    Article  CAS  PubMed  Google Scholar 

  67. Crooks GE, Hon G, Chandonia J-M, Brenner SE. WebLogo: A sequence logo generator. Genome Res. 2004;14(6):1188–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Wang D, Zhang Y, Zhang Z, Zhu J, Yu J. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010;8(1):77–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Wang Z, Liu Y, Li L, Li D, Zhang Q, Guo Y, et al. Whole transcriptome sequencing of Pseudomonas syringae pv. actinidiae-infected kiwifruit plants reveals species-specific interaction between long non-coding RNA and coding genes. Sci Rep. 2017;7(1):4910.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Wang Z, Liu Y, Li D, Li L, Zhang Q, Wang S, et al. Identification of circular RNAs in kiwifruit and their species-specific response to bacterial canker pathogen invasion. Front Plant Sci. 2017;8:413.

    PubMed  PubMed Central  Google Scholar 

  71. Michelotti V, Lamontanara A, Buriani G, Orru L, Cellini A, Donati I, et al. Comparative transcriptome analysis of the interaction between Actinidia chinensis var. chinensis and Pseudomonas syringae pv. actinidiae in absence and presence of acibenzolar-S-methyl. BMC Genomics. 2018;19(1):585.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (grant no. 31960588, 31760567).

Author information

Authors and Affiliations

Authors

Contributions

ZW and CH conceived and designed this study, JT, HJ, MW, WZ and DJ analyzed the data, JT and ZW wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zupeng Wang or Chunhui Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. 

List of the identified TIFY family genes in kiwifruit.

Additional file 2. 

Analysis of cis-elements in promoter regions of kiwifruit TIFYs.

Additional file 3. 

Segmentally and tandemly duplicated kiwifruit TIFY gene pairs.

Additional file 4. 

One-to-one orthologous relationships between A. eriantha and other five plant species.

Additional file 5. 

One-to-one orthologous relationships between A. chinensis and other five plant species.

Additional file 6. 

The homologous relationships between A. eriantha and A. chinensis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tao, J., Jia, H., Wu, M. et al. Genome-wide identification and characterization of the TIFY gene family in kiwifruit. BMC Genomics 23, 179 (2022). https://doi.org/10.1186/s12864-022-08398-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08398-8

Keywords