Skip to main content

Genome-wide identification and expression analysis of the plant-specific PLATZ gene family in Tartary buckwheat (Fagopyrum tataricum)

Abstract

Background

Plant AT-rich sequence and zinc-binding (PLATZ) proteins belong to a novel class of plant-specific zinc-finger-dependent DNA-binding proteins that play essential roles in plant growth and development. Although the PLATZ gene family has been identified in several species, systematic identification and characterization of this gene family has not yet been carried out for Tartary buckwheat, which is an important medicinal and edible crop with high nutritional value. The recent completion of Tartary buckwheat genome sequencing has laid the foundation for this study.

Results

A total of 14 FtPLATZ proteins were identified in Tartary buckwheat and were classified into four phylogenetic groups. The gene structure and motif composition were similar within the same group, and evident distinctions among different groups were detected. Gene duplication, particularly segmental duplication, was the main driving force in the evolution of FtPLATZs. Synteny analysis revealed that Tartary buckwheat shares more orthologous PLATZ genes with dicotyledons, particularly soybean. In addition, the expression of FtPLATZs in different tissues and developmental stages of grains showed evident specificity and preference. FtPLATZ3 may be involved in the regulation of grain size, and FtPLATZ4 and FtPLATZ11 may participate in root development. Abundant and variable hormone-responsive cis-acting elements were distributed in the promoter regions of FtPLATZs, and almost all FtPLATZs were significantly regulated after exogenous hormone treatments, particularly methyl jasmonate treatment. Moreover, FtPLATZ6 was significantly upregulated under all exogenous hormone treatments, which may indicate that this gene plays a critical role in the hormone response of Tartary buckwheat.

Conclusions

This study lays a foundation for further exploration of the function of FtPLATZ proteins and their roles in the growth and development of Tartary buckwheat and contributes to the genetic improvement of Tartary buckwheat.

Peer Review reports

Background

Transcription factors (TFs) are sequence-specific binding proteins that can activate or inhibit the expression of target genes by recognizing and binding to cis-acting elements in their promoter regions of target genes to affect diverse biological processes at the transcriptional level [1]. Zinc finger proteins are an important class of TFs. Based on previous reports, more than 1500 TFs exist in Arabidopsis, accounting for approximately 5% of the Arabidopsis genome [2], of which approximately 15% are zinc finger proteins [3]. The zinc-finger protein consists of two cysteines and two histidines tetrahedrally coordinated with zinc atoms to form a compact finger-like structure. These proteins participate extensively in plant growth and development and actively respond to various stresses [3, 4]. Plant AT-rich sequence and zinc-binding (PLATZ) proteins are a novel class of plant-specific zinc-dependent DNA-binding proteins that preserve the unique structure of the zinc-finger protein family and contain two distantly conserved domains: C-x2-H-x11-C-x2-C-x(4–5)-C-x2-C-x(3–7)-H-x2-H, and C-x2-C-x(10–11)-C-x3-C [5]. Although the first PLATZ gene, PLATZ1, was isolated from peas in 2001 [5], it has attracted increasing attention from researchers. PLATZ1 can nonspecifically bind to A/T-rich sequences and inhibit transcription, as demonstrated by a transient assay [5]. Previous studies have shown that PLATZ proteins play essential roles in several biological processes in plants. For example, Li et al. reported that Floury3 (FL3) encodes a PLATZ protein in maize, which interacts with RNA polymerase III subunit 53 (RPC53) and transcription factor class C 1 (TFC1) to affect endosperm development and filling in seeds [6]. GL6, a PLATZ protein in rice, has been demonstrated to regulate grain length and spikelet number through the same interaction mechanism [7]. SHORT GRAIN6 (SG6) regulates the division of spikelet hull cells and determines seed size in rice by interacting with DP proteins and cell division regulators [8]. Plant regulation by PLATZ is not restricted to the seeds. In Arabidopsis, the PLATZ protein ORESARA15 (ORE15) could regulate leaf growth and senescence by promoting the rate and duration of early cell proliferation [9]. ABA-INDUCED expression 1 (AIN1) represses the elongation of the primary root of Arabidopsis upon ABA induction [10]. In addition, Chao et al. illustrated via transcriptome analysis that PLATZ TFs are important for the secondary growth of Populus stems [11]. Moreover, PLATZ proteins are extensively involved in the response to numerous abiotic stresses, including heat [12], drought [13, 14], salt and osmotic stresses [15, 16], and in response to hormones [17, 18].

Tartary buckwheat (Fagopyrum tataricum, 2n = 2x = 16) is primarily cultivated in Asia, Europe and North America [19]. As a traditional medicinal and edible crop, its grain has a balanced essential amino acid composition and is rich in phytochemicals and soluble fiber [20]. In particular, flavonoids, which have many important biomedical functions, are more abundant in Tartary buckwheat than in other main crops [21,22,23,24]. Tartary buckwheat has been recognized as a green food for humans in the twenty-first century, and has gained popularity among consumers. However, its low yield severely limits its industrial applications [25]. Therefore, identifying PLATZ proteins in Tartary buckwheat is necessary because of their functional potential, particularly their regulatory roles in the growth and development of plant seeds and their relationship to plant resistance, which could provide new insights into the yield improvement of Tartary buckwheat. The PLATZ family has been identified in several other plant species. To date, 12 members have been identified in Arabidopsis thaliana [26], 15 in Oryza sativa [26], 17 in Zea mays [26], 62 in Triticum aestivum [1] and 24 in Brassica rapa [27]. However, to the best of our knowledge, identification and functional characterization of the PLATZ gene family in Tartary buckwheat have not yet been reported. High-quality, chromosome-scale genome sequencing of Tartary buckwheat has recently been completed [28], laying the foundation for a systematic genome-wide study of the PLATZ gene family in Tartary buckwheat. In the present study, 14 PLATZ proteins were identified in Tartary buckwheat genome. We investigated the evolutionary relationships of FtPLATZs together with a comprehensive study of gene structures, conserved motif composition, and cis-acting elements in the promoter regions of FtPLATZs. Gene duplication events and their syntenic relationships with the six representative species were investigated. For functional characterization, we examined the expression profiles of FtPLATZs in different tissues of Tartary buckwheat and in grains at different developmental stages using real-time quantitative polymerase chain reaction (qRT-PCR). In addition, the responses of FtPLATZs to various exogenous hormones were investigated. This study aimed to form a foundation for further exploration of the functional mechanisms of FtPLATZs and contribute to the improvement of plant varieties and innovation of the germplasm in Tartary buckwheat.

Results

Identification of FtPLATZ proteins in Tartary buckwheat

Combining the results of the hidden Markov model (HMM) search and BLASTP operations and further examination of the conserved PLATZ domain, 14 putative FtPLATZ proteins were identified in Tartary buckwheat (Fig. S1 and Table S1). They were unevenly distributed on six chromosomes of Tartary buckwheat (Fig. 1). Chromosome Ft4 contained the largest number of FtPLATZ genes (four genes), followed by Ft1 and Ft8, both of which contained three genes. Ft3 contained two genes, whereas Ft2 and Ft6 contained only one. The FtPLATZ genes were not found on chromosomes Ft5 and Ft7. We designated these as FtPLATZ1 to FtPLATZ14 based on their location on the chromosomes. As shown in Table 1, the full-length cDNAs, predicted protein products and Mw of FtPLATZ genes varied greatly, ranging from 444 to 1590 bp, 148 to 530 aa, and 16.59 to 58.82 kDa, respectively. The average coding sequence (CDS) length, predicted protein products, and molecular weight (Mw) were 805 bp, 268 aa, and 30.20 kDa, respectively. The data clearly showed that FtPLATZ3 was the smallest, and FtPLATZ10 exhibited the largest size with the maximum level of CDS length, predicted protein products, and Mw among FtPLATZs. The difference in the theoretical isoelectric point (pI) values among FtPLATZ genes was relatively small, with an average of 8.80.

Fig. 1
figure 1

Schematic diagram of chromosomal distribution of FtPLATZ genes in Tartary buckwheat. The vertical bars represent the chromosomes of Tartary buckwheat, and the scale for chromosome length is shown on the left. The genes marked in red represent tandem duplication events

Table 1 PLATZ family genes in Tartary buckwheat

Phylogenetic analysis and classification of FtPLATZ proteins

To clarify the evolutionary relationship between the PLATZ proteins of Tartary buckwheat and the PLATZ proteins of two model plants, Arabidopsis and rice, we constructed a maximum likelihood (ML) tree with the 14 identified FtPLATZs, 12 AtPLATZs and 15 OsPLATZs (Fig. 2). The 41 PLATZ proteins were divided into five groups (I to V), and the FtPLATZ proteins were distributed in the four main groups (II to V). Group II contained the largest number of FtPLATZ members (6 of 14, 42.86%). Half of the PLATZ proteins in Group II were FtPLATZs. Group V contained four FtPLATZs, whereas the remaining 10 proteins were from Arabidopsis and rice. Group IV contained one FtPLATZ, one AtPLATZ member, and five OsPLATZs. In particular, group III was only composed of three FtPLATZ members, indicating no homology to AtPLATZs and OsPLATZs. Group I contained two AtPLATZs and three OsPLATZs but no FtPLATZ proteins. In addition, a phylogenetic tree for FtPLATZs was constructed and labeled based on the grouping in the overall phylogenetic tree to analyze the differences in gene structure and motif components among groups (Fig. 3a).

Fig. 2
figure 2

Unrooted phylogenetic tree constructed by the maximum likelihood method for PLATZ genes of Tartary buckwheat, Arabidopsis and rice

Fig. 3
figure 3

Phylogenetic relationship, gene structure and motif composition of FtPLATZ genes of Tartary buckwheat. a. Phylogenetic tree for PLATZ genes of Tartary buckwheat was constructed using the maximum likelihood method. b. Exon–intron structure of FtPLATZ genes. The legend is shown in the upper-right corner, and the introns are represented by grey lines. c. Motif composition of FtPLATZ genes. Different motifs are represented by different colors, as indicated in the legend on the right

Gene structure and conserved motifs analysis of FtPLATZ genes

The exon–intron structure of the FtPLATZ genes was investigated based on the genomic DNA sequence of Tartary buckwheat to understand the structural composition of the FtPLATZ genes (Fig. 3b). In general, the structures of FtPLATZ genes were distinguishable among the phylogenetic groups, and they showed similar characteristics within the groups. Most genes contained three introns (9 out of 14, 64.29%), and only five genes, FtPLATZ1/7/8/10/13, contained four introns. Group III was characterized by FtPLATZ genes with four introns, whereas groups IV and V contained only three-intron genes. In Group II, two genes, FtPLATZ1 and FtPLATZ13, had four introns, whereas the other genes had three introns.

Result similar to the exon–intron structure was also found in the motif composition of phylogenetically grouped FtPLATZ members (Table S2). As shown in Fig. 3c, motifs 1, 2, 3, and 6, which constituted the core domain of PLATZ, were universally present in the FtPLATZ members, except for one gene (FtPLATZ3) in group II, where motif 2 was not present, indicating a possible sequence loss during evolution. In addition, motifs 5 and 9 were uniquely present in group III, and motifs 4 and 8 were uniquely present in group V. Motif 10 appeared only in group II, and motifs 4 and 8 were present separately in the five FtPLATZs of group II. Only one member of group IV, FtPLATZ12, possessed motif 4 exclusively, in addition to the core domain of PLATZ.

Gene duplication events and synteny analysis of FtPLATZ genes

Possible gene duplication events among the FtPLATZs were investigated to explore the evolution of FtPLATZ genes. The results showed that tandem duplication and segmental duplication events were observed in FtPLATZs, where FtPLATZ7/FtPLATZ8 formed a tandem duplication event (Fig. 1) and FtPLATZ1/FtPLATZ13, FtPLATZ5/FtPLATZ9 and FtPLATZ7/FtPLATZ10 formed three segmental duplication events (Fig. 4). These results indicate that duplication events widely participated in the evolution of FtPLATZs.

Fig. 4
figure 4

Schematic diagram of the syntenic relationships of FtPLATZ genes in Tartary buckwheat. The grey ribbons represent syntenic blocks in the Tartary buckwheat genome, and the segmental duplication events are marked in red

Furthermore, we investigated the syntenic relationships between FtPLATZs and PLATZ genes from four representative dicotyledons (A. thaliana, G. max, V. vinifera, and S. lycopersicum) and two representative monocotyledons (O. sativa and Z. mays; Fig. 5). The number of orthologous gene pairs between Tartary buckwheat and the other six species was quite different: five pairs with Arabidopsis, ten with soybean, six with grape, five with tomato, one with rice, and one with maize (Table S3). In particular, among the 14 FtPLATZ genes, FtPLATZ5 (FtPinG0008634800.01.T01) was the only gene that was collinear with PLATZ proteins of the six representative plants. FtPLATZ5 was collinear with at least two PLATZ genes in dicotyledons and one in monocotyledons. The results indicated that these orthologous genes may exist before the differentiation of the ancestors.

Fig. 5
figure 5

Synteny analysis of PLATZ genes between Tartary buckwheat and the other six representative plants. The syntenic gene pairs were linked by red lines

Expression patterns of FtPLATZ genes in different tissues and grain developmental stages of Tartary buckwheat

The potential roles of the identified FtPLATZ genes in the growth and development of Tartary buckwheat were explored using qRT-PCR (Fig. 6a and Table S4). In general, the expression patterns of FtPLATZ genes varied greatly in different tissues, indicating their potential multiple functions in the growth and development of Tartary buckwheat. Two genes (FtPLATZ4 and FtPLATZ11) showed similar expression patterns, specifically expressed in the roots and slightly expressed in grains. Three genes (FtPLATZ6, FtPLATZ9, and FtPLATZ12) showed the highest expression levels in the stems. FtPLATZ5 showed the highest expression levels in the leaves and stems. Four genes (FtPLATZ1, FtPLATZ2, FtPLATZ7, and FtPLATZ13) were highly expressed in the flowers, whereas FtPLATZ2 and FtPLATZ13 were only slightly expressed in other tissues. In addition, four genes (FtPLATZ3, FtPLATZ8, FtPLATZ10, and FtPLATZ14) were highly expressed in the grains, reaching their highest expression levels successively in the S1, S2, S3, and S4 developmental stages of the grains.

Fig. 6
figure 6

Expression profiles of FtPLATZ genes in different tissues and different grain developmental stages of Tartary buckwheat and correlation analysis of the expression patterns of FtPLATZ genes. a. Expression profiles of 14 FtPLATZ genes in the root (R), stem (S), leaf (L), flower (F), and grain (G-S1, initial formation stage; G-S2, green grain stage; G-S3, discoloration stage; G-S4, initial maturity stage) of Tartary buckwheat. Error bars are obtained from three biological replicates. Lowercase letters above the bars represent significant differences between different treatments, as determined by Duncan’s multiple range test (p < 0.05). b. Pearson’s correlation of expression patterns among FtPLATZ genes. Red and blue represent positive and negative correlations, respectively. * and ** indicate significance at the levels of 0.05 and 0.01, respectively

The expression patterns of FtPLATZ genes in different developmental stages of Tartary buckwheat grains have drawn much attention. Six patterns are identified. The expression levels of FtPLATZ5, FtPLATZ11, and FtPLATZ14 increased with grain growth and development, whereas those of FtPLATZ1, FtPLATZ2, and FtPLATZ7 decreased with grain growth and development. In addition to the two monotonous expression patterns, some gene expression levels initially decreased and then increased (FtPLATZ6 and FtPLATZ9), and some initially increased and then decreased with grain development (FtPLATZ4, FtPLATZ8, and FtPLATZ12). In addition, the expression of FtPLATZ10 and FtPLATZ13 showed a wave-shaped trend, decreasing twice during stages S2 and S4 of the grains. In particular, FtPLATZ3 was highly expressed in the S1 stage, but not in the other stages. Collectively, FtPLATZ genes may play crucial roles during grain development in Tartary buckwheat.

Further correlation analysis indicated that the expression patterns of some FtPLATZ genes in different tissues of Tartary buckwheat and different developmental stages of the grain were significantly and positively correlated (Fig. 6b); that is, FtPLATZ1/FtPLATZ2 (p < 0.05), FtPLATZ1/FtPLATZ7 (p < 0.05), FtPLATZ2/FtPLATZ13 (p < 0.01), FtPLATZ4/FtPLATZ11 (p < 0.01), and FtPLATZ6/FtPLATZ12 (p < 0.01), which was consistent with the results shown in Fig. 6a, indicating that some FtPLATZ genes may act synergistically with one another during development.

Analysis of promoter cis-acting elements of FtPLATZ genes

The functional potential of the identified FtPLATZ genes was further explored by investigating cis-acting elements in the promoter regions of these genes. Various cis-acting elements were identified, as summarized in Table S5. Promoter-related elements (i.e., TATA-box and CAAT-box) and light-responsive elements (i.e., Box 4, G-Box, TCT-motif et al.) were most abundantly distributed in the promoter region of FtPLATZ genes. Notably, stress-related elements (i.e., ARE, LTR, and MBS) and hormone-responsive elements (i.e., ABRE, CGTCA-motif, and TCA-element) were also widely distributed in the promoter region of the FtPLATZ genes. In particular, the number of hormone-responsive elements varied considerably among the FtPLATZ genes (Fig. 7), suggesting that the 14 FtPLATZs may function specifically in response to different hormone stimulation. In addition, some development-related elements (i.e., O2-site, MSA-like and CAT-box) and site-binding-related elements (i.e., CCAAT-box, HD-Zip 3 and MBSI) were identified in the promoter region of the FtPLATZ genes, but not all FtPLATZs contained such elements.

Fig. 7
figure 7

Distribution of the cis-acting elements related to hormone response in the promoter region of FtPLATZ genes

Differential expression of FtPLATZ genes under different exogenous hormone treatments

The expression levels of the 14 identified FtPLATZ genes after treatment with five exogenous hormones and the control, were compared using qRT-PCR to investigate the response pattern of FtPLATZ genes to hormones (Fig. 8a and Table S6). The results showed that the expression levels of 11 of the 14 FtPLATZ genes were altered significantly after treatment with at least one type of exogenous hormone. MeJA treatment had the greatest impact on FtPLATZ genes among the five hormones, with significant upregulation of FtPLATZ2, FtPLATZ4, FtPLATZ6, and FtPLATZ9 and downregulation of FtPLATZ5, FtPLATZ7, FtPLATZ12, and FtPLATZ13. FtPLATZ genes, such as FtPLATZ3, FtPLATZ4, FtPLATZ5, FtPLATZ6, and FtPLATZ14, which responded significantly to SA treatment, were primarily upregulated. Only one gene, FtPLATZ12, was downregulated. Similar results were found for IAA and ABA treatments, in which the genes were primarily upregulated. In addition, under GA treatment, three genes were downregulated (FtPLATZ5, FtPLATZ9, and FtPLATZ14), and two genes were upregulated (FtPLATZ4 and FtPLATZ6). Notably, FtPLATZ6 was significantly upregulated by all five exogenous hormones, particularly ABA, GA, and SA. Moreover, FtPLATZ5 and FtPLATZ14 responded significantly to four hormones, and they showed similar response patterns to GA, IAA, and SA.

Fig. 8
figure 8

Expression profile of FtPLATZ genes under treatment with different exogenous hormones and their correlation analysis. a. Expression profiles of the 14 FtPLATZ genes after treatment with methyl jasmonate (MeJA), abscisic acid (ABA), gibberellin (GA), indole-3-acetic acid (IAA), and salicylic acid (SA) and the same amount of water as the control. Error bars are obtained from three biological replicates. The asterisks above the bars represent the level of significance of the expression differences under different exogenous hormone treatments compared with the control group, as determined by Student’s t-test. *, ** and *** indicate significance at the levels of 0.05, 0.01 and 0.001, respectively. b. Pearson’s correlation of response patterns to exogenous hormones among 14 FtPLATZ genes. Red and blue represent positive and negative correlations, respectively. * and ** indicate significance at the levels of 0.05 and 0.01, respectively

Furthermore, the results of the correlation analysis showed that the expression patterns of some genes were significantly correlated after treatment with exogenous hormones (Fig. 8b). The expression patterns of FtPLATZ2 and FtPLATZ4 were significantly and positively correlated (p < 0.01), whereas FtPLATZ3 and FtPLATZ7 showed a significant negative correlation (p < 0.05).

Subcellular localization of FtPLATZ proteins

The subcellular localization prediction results of CELLO and Plant-mPLoc consistently showed that most proteins were localized in the nucleus (11 out of 14), but the prediction results for the remaining three proteins (FtPLATZ1, FtPLATZ2, and FtPLATZ12) were inconsistent between the two prediction methods (Table 1). CELLO predicted that these proteins were localized extracellularly, whereas Plant-mPLoc predicted that they were located in the nucleus. Transient expression in Nicotiana benthamiana was examined to verify its subcellular localization (Fig. 9). These results indicated that the green fluorescent protein (GFP) fluorescent signals of the three fusion proteins were primarily localized in the nucleus. By contrast, the control 35S::GFP signal was detected in whole cells. The experimental results suggest that FtPLATZ proteins may function as conventional TFs.

Fig. 9
figure 9

Subcellular localization of FtPLATZ proteins. The control (35S::GFP), 35S::FtPLATZ1-GFP, 35S::FtPLATZ2-GFP, and 35S::FtPLATZ12-GFP fusion proteins were transiently expressed in Nicotiana benthamiana leaves, separately. GFP, green fluorescence of fusion proteins; NLS-mCherry, red fluorescence of the nucleus; Bright, bright field; Merged, merged microscopic images. Scale bars = 20 μm

Discussion

PLATZ TFs are a class of plant-specific zinc-dependent DNA-binding proteins that play important roles in the growth and development of plants and their response to stress [16]. In this study, we identified 14 PLATZ proteins in the Tartary buckwheat genome, all of which harbored conserved PLATZ domains. The amount of PLATZ proteins in Tartary buckwheat was similar to that identified in Arabidopsis (12) [26], rice (15) [26], and maize (17) [26]. However, the genome size of these species varied greatly (Tartary buckwheat, 489.3 Mb [28]; Arabidopsis, 125 Mb [29]; rice, 466 Mb [30], and maize, 2.3Gb [31]), implying that the amount of PLATZ proteins and the size of the genome were not closely related.

Gene duplication, including tandem duplication and segmental duplication, is regarded as a primary driving force in the evolution of genomes and genetic systems and is also a mechanism for organisms to adapt to changing environments [32, 33]. Fu et al. found that 21 of the 62 TaPLATZ genes identified in the wheat genome were from tandem duplications (33.9%) and two from segmental duplications and concluded that genomic duplication was the primary cause of the expansion of the TaPLATZ family [1]. Similarly, Azim et al. found a considerable number of gene duplication events in Brassica rapa, where 20 pairs of segmental duplication genes were detected among the 24 identified BrPLATZ genes, whereas no tandem duplication events were found [27]. In our study, a pair of tandem duplicated FtPLATZ genes (Fig. 1) and three pairs of segmental duplicated FtPLATZ genes (Fig. 4) were detected in Tartary buckwheat, accounting for 50% of the FtPLATZ genes (seven out of 14 genes), implying that gene duplication was the main driving force in the evolution of FtPLATZ genes. These duplicated genes had almost the same exon–intron structure and motif composition (Fig. 3b and c), but their expression preferences seemed to differ (Fig. 6a). Subfunctionalization of duplicated FtPLATZ genes may account for their different expression patterns [34]. Furthermore, synteny analysis showed that the PLATZ genes of Tartary buckwheat shared more orthologs with dicotyledons than with monocotyledons. Tartary buckwheat and soybean had the largest number of orthologous gene pairs (Table S3), implying that they could have a closer evolutionary relationship and may have evolved from a common ancestor, which conformed to previous findings [20, 35].

In the phylogenetic analysis, the 41 PLATZ proteins obtained from Tartary buckwheat, Arabidopsis and rice were classified into five groups based on their phylogenetic relationships, wherein 14 FtPLATZ proteins were distributed into four main groups (Groups II to V, Fig. 2). The exon–intron structures of FtPLATZ genes were similar, containing three or four introns (Fig. 3b), implying that FtPLATZ genes were relatively conserved during evolution [27]. FtPLATZ genes within the same group shared a similar gene structure and motif composition, whereas evident distinctions were found among different groups, particularly in motif composition, implying large functional differentiation of FtPLATZ genes. Ten motifs were detected in the FtPLATZ proteins, of which motifs 1, 2, 3, and 6 constituted the PLATZ domain (Fig. 3c). In our study, exon loss was observed in the FtPLATZ genes. In particular, FtPLATZ3 did not contain motif 2, part of the beginning of the PLATZ domain, which may be due to genetic variation that occurred during the evolution of FtPLATZ genes, thereby leading to the alteration of gene functions [36]. LOC_Os06g45540.1 (SG6; GL6) belonged to the same phylogenetic group as FtPLATZ3 (group II), which has been proven to regulate the grain size and spikelet number of rice [7, 8]. Meanwhile, time-course transcriptome analysis revealed that AT3G60670.1 in group II was involved in the development and maturation of Arabidopsis grains [37]. As shown in Fig. 6a, FtPLATZ3 was significantly expressed at the S1 stage of grains, which is considered a critical developmental period for grain size [35]. Collectively, we hypothesized that FtPLATZ3 may be involved in the regulation of grain size in Tartary buckwheat, and further experimental verification is necessary. In addition, FtPLATZ4 and FtPLATZ11 may play important roles in the development of Tartary buckwheat roots. Although they have been found to be specifically expressed in the roots through tissue expression profiles, AT2G12646.1 (RITF1), located in the closest phylogenetic branch with two FtPLATZ genes, has been demonstrated to play a central role in mediating root meristem growth factor 1 (RGF1) signalling and subsequently affecting the size of root meristems [38].

Plant hormones play important roles in numerous biological processes and contribute remarkably to the adaptability of plants to changing environments [39, 40]. Previous studies have shown that PLATZ genes are hormone responsive. GmPLATZ1 in soybeans [18] and PLATZ genes in Thellungiella salsuginea roots [17] could be induced by ABA. GhPLATZ1 from cotton is significantly upregulated in transgenic Arabidopsis under ABA and GA treatments [15]. Moreover, ABA can induce the expression of AIN1 in Arabidopsis, thereby affecting the elongation of the primary root [10]. PhePLATZ genes in moso bamboo were significantly regulated by GA, ABA, and MeJA treatments [41]. In the present study, hormone-responsive elements related to ABA, MeJA, SA, GA, and IAA were examined in the promoter region of FtPLATZ genes; however, their distribution across different FtPLATZ genes was diverse, which is similar to the findings of Fu et al. in the identification of TaPLATZs [1]. MeJA can activate the expression of defense genes, induce the synthesis of defensive compounds, and can also affect the antioxidant system [42]. Numerous studies have revealed that MeJA is involved in mediating defense responses against fungal pathogens [43], alleviating salt [44], drought [45], and chilling stresses [46]. The expression of nearly 80% of FtPLATZ genes (11 out of 14) changed significantly after treatment with exogenous hormones (Fig. 8a), among which MeJA treatment exhibited the widest effects in 8 of the 11 significantly disturbed genes, implying that FtPLATZ genes might be extensively involved in the stress response of Tartary buckwheat. FtPLATZ6 was the only gene that was significantly upregulated after all exogenous hormone treatments, which may indicate the critical role of FtPLATZ6 in biological processes involved in the hormone response of Tartary buckwheat. However, no abundant hormone-responsive elements were found in the promoter region of FtPLATZ6 (Fig. 7). Previous studies have reported that the distribution pattern of cis-acting elements is not directly related to the gene expression levels [47,48,49]. Therefore, the expression of FtPLATZ genes may involve complex regulatory mechanisms that require further experimental verification.

Conclusions

In this study, we systematically identified and characterized the PLATZ gene family in Tartary buckwheat. Fourteen FtPLATZ proteins were identified, which were unevenly distributed on six of the eight chromosomes in Tartary buckwheat. Based on phylogenetic analysis, the FtPLATZ proteins were classified into four groups, and each group shared a similar gene structure and motif composition. In addition, gene duplication, particularly segmental duplication, was the main driving force in the evolution of the FtPLATZ genes. We analyzed the expression levels of 14 FtPLATZ genes in different tissues and different grain developmental stages of Tartary buckwheat and their responses to five exogenous hormones. The results revealed, to a great extent, the important roles of FtPLATZ genes in the growth and development of Tartary buckwheat, such as FtPLATZ3, which might be involved in the regulation of grain size; FtPLATZ4 and FtPLATZ11, which played a role in root development; and FtPLATZ6, which was significantly upregulated after all the exogenous hormone treatments and may be critical for the hormone response of Tartary buckwheat. This study provides a foundation for further exploration of the functional characteristics of FtPLATZ genes and promotes targeted genetic breeding research for crop improvement in Tartary buckwheat.

Material and methods

Identification of FtPLATZ genes in Tartary buckwheat genome

The Tartary buckwheat genome was obtained from the Tartary Buckwheat Genome Project (TBGP; http://www. mbkbase.org/Pinku1/), and the gene anotation V2 version was used for subsequent analysis [28]. To identify PLATZ genes in the Tartary buckwheat genome, the HMM profile of PLATZ (PF04640) downloaded from the Pfam database (http://pfam.xfam.org/) was used to search against the Tartary buckwheat genome database via HMMER3.3 with default parameter settings [50]. The PLATZ genes of Arabidopsis and rice obtained from the TAIR database (https://www.arabidopsis.org/) and iTAK database (http://itak.feilab.net/cgi-bin/itak/index.cgi) [51], respectively, were used to perform a BLASTP operation to further retrieve possible FtPLATZ genes from the Tartary buckwheat genome with a score ≥ 100 and e-value ≤1 × 10 − 10 [52]. All putative FtPLATZ genes integrating the results of the HMM retrieval and BLASTP operations were submitted to the NCBI Conserved Domain Database (CDD, https://www.ncbi.nlm.nih.gov/cdd), SMART (http://smart.embl-heidelberg.de/), and Pfam to examine the existence of the conserved PLATZ domain.

Sequence characterization

We collected chromosomal location information for the identified FtPLATZ genes from the Tartary buckwheat genome database and visualized them using TBtools [53]. The properties of the identified FtPLATZ genes, including CDS length, protein length, Mw, and pI, were investigated using Expasy (https://web.expasy.org/compute_pi/). The exon–intron structure of the FtPLATZ members was investigated using TBtools based on Tartary buckwheat genome annotation information. Conserved motifs of FtPLATZ proteins were identified using the MEME Suite (https://meme-suite.org/meme/tools/meme) with default parameters, except that the maximum number of motifs was set to 10. Moreover, the cis-acting elements within the 2000 bp sequence upstream of FtPLATZ genes, usually regarded as the promoter region of a gene [1], were analyzed using PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [54].

Phylogenetic analysis

The identified genes, together with AtPLATZs and OsPLATZs, were used to construct a phylogenetic tree using the ML method with the MEGA X software [55]. The Jones–Taylor–Thornton (JTT) model combined with a discrete gamma distribution (+ G) was selected as the optimal model for constructing the phylogenetic tree. Sequences with more than 20% alignment gaps were removed, and a bootstrap test was conducted with 1000 replicates. A phylogenetic tree containing only FtPLATZs was constructed using these parameters. The classification of FtPLATZs based on the phylogenetic tree was referred to the method described by Wang et al. [26].

Gene duplication and synteny analysis

Possible gene duplication events among FtPLATZ genes were probed using multiple collinear scanning toolkits (MCScanX) [56]. Syntenic analyses were conducted using TBtools between the identified FtPLATZ proteins and PLATZ protein sequences of Glycine max, Vitis vinifera, Solanum lycopersicum, Oryza sativa, and Zea mays obtained from the iTAK database and the AtPLATZs obtained from the TAIR database.

Plant materials and treatments

Weining-14, a Tartary buckwheat variety provided by the Minor Grain Crops Research Centre of Northwest A & F University, was planted in the experimental field of Northwest A & F University, Yangling, Shaanxi, China in 2020. Different Tartary buckwheat tissues, including the roots, stems, leaves, flowers, and grains at different developmental stages (3, 10, 17, and 24 days after pollination, corresponding to the initial formation stage [G_S1], green grain stage [G_S2], discoloration stage [G_S3], and initial maturity stage [G_S4], respectively) were sampled.

To investigate the response of FtPLATZ genes to exogenous hormones, 21-day-old seedlings (Fig. S2) were treated with different exogenous hormones, including 100 μM methyl jasmonate (MeJA), abscisic acid (ABA), salicylic acid (SA), 10 μM indole-3-acetic acid (IAA), and gibberellin (GA) by foliar spraying. The control group was sprayed with equal amounts of water. After 6 h of treatment, the second leaves of the seedlings were collected separately [57]. All the samples were collected from at least three healthy plants, immediately frozen with liquid nitrogen, and then stored at − 80 °C for RNA extraction and subsequent qRT-PCR analysis.

Expression analyses of FtPLATZ genes by qRT-PCR

Total RNA was extracted from all samples using a MiniBEST Plant RNA Extraction Kit (TaKaRa). First-strand cDNA was synthesized using the PrimeScript™ II 1st Strand cDNA Synthesis Kit (TaKaRa). qRT-PCR was performed using TB Green™ Premix Ex Taq™ II (TaKaRa) on a Q7 Real-Time PCR System (Applied Biosystems™, Foster City, CA, USA) following the manufacturer’s instructions. The primers for qRT-PCR analysis were designed using Primer3 software (version 4.1.0, https://primer3.ut.ee/) based on the CDSs of the identified FtPLATZ genes obtained from TBGP, and the information of all primer sequences are listed in Table S7. FtH3 was selected as the internal reference gene, which has been proven to be stably expressed in Tartary buckwheat under any condition [58]. Expression data were analyzed using the 2−ΔΔCT method [59].

Subcellular localization of FtPLATZ proteins

The subcellular localization of the identified FtPLATZ proteins was predicted using CELLO (version 2.5, http://cello.life.nctu.edu.tw/) [60] and Plant-mPLoc (version 2.0, http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/) [61]. Three proteins, namely, FtPLATZ1, FtPLATZ2, and FtPLATZ12, with inconsistent results (CELLO predicted as extracellular, whereas Plant-mPLoc predicted as nuclear) were selected to verify the prediction of subcellular localization. The CDSs of FtPLATZs (excluding stop codons) were cloned from the cDNA for qRT-PCR using the primers listed in Table S7 and then inserted into the pCAMBIA2300-GFP vector driven by a 35S promoter. The recombinant plasmids 35S::FtPLATZ1-GFP, 35S::FtPLATZ2-GFP, and 35S::FtPLATZ12-GFP were constructed and transformed into Agrobacterium tumefaciens strain GV3101 (Shanghai Weidi Biotechnology Co., Ltd., Shanghai, China). Transient expression was performed in N. benthamiana leaves in accordance with the method of Fu et al. [1], and the GFP fluorescence signal was detected by confocal laser scanning microscopy (LSM880; Carl Zeiss, Germany).

Statistical analysis

Comparisons of the expression levels of the FtPLATZ genes in different tissues were statistically evaluated by one-way analysis of variance (ANOVA) using IBM SPSS Statistics 25 (IBM Corporation, Armonk, NY) [62]. Duncan’s multiple range test was used to determine significant differences between groups. Student’s t-test was carried out using R software (version 4.0.2) to examine whether the expression of FtPLATZ genes changed significantly after stimulation with exogenous hormones.

Availability of data and materials

The genome sequences of Tartary buckwheat used for identifying PLATZ genes in this study were located in the Tartary Buckwheat Genome Project (TBGP; http://www.mbkbase.org/Pinku1/). The Tartary buckwheat accession (Weining-14) used in the experiment was provided by the Minor Grain Crops Research Centre of Northwest A & F University. The datasets supporting the conclusions of this article are included in the article and its Supplementary files.

Abbreviations

aa:

Amino acid

ABA:

Abscisic acid

AIN1:

ABA-INDUCED expression 1

ANOVA:

One-way analysis of variance

At:

Arabidopsis thaliana

bp:

Base pair

Br:

Brassica rapa

CDD:

Conserved Domain Database

CDS:

Coding sequence

FL3:

Floury3

Ft:

Fagopyrum tataricum

GA:

gibberellin

GFP:

Green fluorescent proteins

Gh:

Gossypium hirtusum

Gm:

Glycine max

G_S1:

Initial formation stage

G_S2:

Green grain stage

G_S3:

Discoloration stage

G_S4:

Initial maturity stage

HMM:

Hidden Markov model

IAA:

Indole-3-acetic acid

JTT:

Jones-Taylor-Thornton

MeJA:

Methyl jasmonate

ML:

Maximum likelihood

Mw:

Molecular weight

ORE15:

ORESARA15

Os:

Oryza sativa

Phe:

Phyllostachys edulis

pI:

Isoelectric point

PLATZ:

Plant AT-rich sequence and zinc-binding

qRT-PCR:

Real-time quantitative polymerase chain reaction

RGF1:

Root meristem growth factor 1

RPC53:

RNA polymerase III subunit 53

SA:

Salicylic acid

SG6:

SHORT GRAIN6

Ta:

Triticum aestivum

TBGP:

Tartary Buckwheat Genome Project

TFC1:

Transcription factor class C 1

TFs:

Transcription factors

References

  1. Fu Y, Cheng M, Li M, Guo X, Wu Y, Wang J. Identification and characterization of PLATZ transcription factors in wheat. Int J Mol Sci. 2020;21(23):8934.

    Article  CAS  PubMed Central  Google Scholar 

  2. Riechmann JL, Heard J, Martin G, Reuber L, Jiang CZ, Keddie J, et al. Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000;290(5499):2105–10.

    Article  CAS  PubMed  Google Scholar 

  3. Kielbowicz-Matuk A. Involvement of plant C2H2-type zinc finger transcription factors in stress responses. Plant Sci. 2012;185-186:78–85.

    Article  CAS  PubMed  Google Scholar 

  4. Takatsuji H. Zinc-finger proteins: the classical zinc finger emerges in contemporary plant science. Plant Mol Biol. 1999;39(6):1073–8.

    Article  CAS  PubMed  Google Scholar 

  5. Nagano Y, Furuhashi H, Inaba T, Sasaki Y. A novel class of plant-specific zinc-dependent DNA-binding protein that binds to a/T-rich DNA sequences. Nucleic Acids Res. 2001;29(20):4097–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Li Q, Wang J, Ye J, Zheng X, Xiang X, Li C, et al. The maize imprinted gene Floury3 encodes a PLATZ protein required for tRNA and 5S rRNA transcription through interaction with RNA polymerase III. Plant Cell. 2017;29(10):2661–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Wang A, Hou Q, Si L, Huang X, Luo J, Lu D, et al. The PLATZ transcription factor GL6 affects grain length and number in Rice. Plant Physiol. 2019;180(4):2077–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhou SR, Xue HW. The rice PLATZ protein SHORT GRAIN6 determines GRAIN size by regulating spikelet hull cell division. J Integr Plant Biol. 2020;62(6):847–64.

    Article  CAS  PubMed  Google Scholar 

  9. Kim JH, Kim J, Jun SE, Park S, Timilsina R, Kwon DS, et al. ORESARA15, a PLATZ transcription factor, mediates leaf growth and senescence in Arabidopsis. New Phytol. 2018;220(2):609–23.

    Article  CAS  PubMed  Google Scholar 

  10. Dong T, Yin X, Wang H, Lu P, Liu X, Gong C, et al. ABA-INDUCED expression 1 is involved in ABA-inhibited primary root elongation via modulating ROS homeostasis in Arabidopsis. Plant Sci. 2021;304:110821.

    Article  CAS  PubMed  Google Scholar 

  11. Chao Q, Gao ZF, Zhang D, Zhao BG, Dong FQ, Fu CX, et al. The developmental dynamics of the Populus stem transcriptome. Plant Biotechnol J. 2019;17(1):206–19.

    Article  CAS  PubMed  Google Scholar 

  12. Blair EJ, Bonnot T, Hummel M, Hay E, Marzolino JM, Quijada IA, et al. Contribution of time of day and the circadian clock to the heat stress responsive transcriptome in Arabidopsis. Sci Rep. 2019;9(1):4814.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Zenda T, Liu S, Wang X, Liu G, Jin H, Dong A, et al. Key maize drought-responsive genes and pathways revealed by comparative transcriptome and physiological analyses of contrasting inbred lines. Int J Mol Sci. 2019;20(6):1268.

    Article  CAS  PubMed Central  Google Scholar 

  14. Gonzalez-Morales SI, Chavez-Montes RA, Hayano-Kanashiro C, Alejo-Jacuinde G, Rico-Cambron TY, de Folter S, et al. Regulatory network analysis reveals novel regulators of seed desiccation tolerance in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2016;113(35):E5232–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhang S, Yang R, Huo Y, Liu S, Yang G, Huang J, et al. Expression of cotton PLATZ1 in transgenic Arabidopsis reduces sensitivity to osmotic and salt stress for germination and seedling establishment associated with modification of the abscisic acid, gibberellin, and ethylene signalling pathways. BMC Plant Biol. 2018;18(1):218.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Liu S, Yang R, Liu M, Zhang S, Yan K, Yang G, et al. PLATZ2 negatively regulates salt tolerance in Arabidopsis seedlings by directly suppressing the expression of the CBL4/SOS3 and CBL10/SCaBP8 genes. J Exp Bot. 2020;71(18):5589–602.

    Article  CAS  PubMed  Google Scholar 

  17. Zhang Y, Shi SH, Li FL, Zhao CZ, Li AQ, Hou L, et al. Global transcriptome analysis provides new insights in Thellungiella salsuginea stress response. Plant Biol. 2019;21(5):796–804.

    Article  CAS  PubMed  Google Scholar 

  18. So HA, Choi SJ, Chung E, Lee JH. Molecular characterization of stress-inducible PLATZ gene from soybean (Glycine max L.). Plant Omics. 2015;8(6):479–84.

    CAS  Google Scholar 

  19. Li J, Feng S, Qu Y, Gong X, Luo Y, Yang Q, et al. Identifying the primary meteorological factors affecting the growth and development of Tartary buckwheat and a comprehensive landrace evaluation using a multi-environment phenotypic investigation. J Sci Food Agric. 2021;101(14):6104–16.

    Article  CAS  PubMed  Google Scholar 

  20. Liu M, Sun W, Ma Z, Huang L, Wu Q, Tang Z, et al. Genome-wide identification of the SPL gene family in Tartary buckwheat (Fagopyrum tataricum) and expression analysis during fruit development stages. BMC Plant Biol. 2019;19(1):299.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Li J, Yang P, Yang Q, Gong X, Ma H, Dang K, et al. Analysis of flavonoid metabolites in buckwheat leaves using UPLC-ESI-MS/MS. Molecules. 2019;24(7):1310.

    Article  CAS  PubMed Central  Google Scholar 

  22. Li J, Hossain MS, Ma HC, Yang QH, Gong XW, Yang P, et al. Comparative metabolomics reveals differences in flavonoid metabolites among different coloured buckwheat flowers. J Food Compos Anal. 2020;85:103335.

    Article  CAS  Google Scholar 

  23. Adom KK, Liu RH. Antioxidant activity of grains. J Agric Food Chem. 2002;50(21):6182–7.

    Article  CAS  PubMed  Google Scholar 

  24. Guo XD, Ma YJ, Parry J, Gao JM, Yu LL, Wang M. Phenolics content and antioxidant activity of tartary buckwheat from different locations. Molecules. 2011;16(12):9850–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Feng S, Li J, Qian GQ, Feng BL. Association between the yield and the main agronomic traits of Tartary buckwheat evaluated using the random forest model. Crop Sci. 2020;60(5):2394–407.

    Article  CAS  Google Scholar 

  26. Wang J, Ji C, Li Q, Zhou Y, Wu Y. Genome-wide analysis of the plant-specific PLATZ proteins in maize and identification of their general role in interaction with RNA polymerase III complex. BMC Plant Biol. 2018;18(1):221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Azim JB, Khan MFH, Hassan L, Robin AHK. Genome-wide characterization and expression profiling of plant-specific PLATZ transcription factor family genes in Brassica rapa L. Plant Breed Biotech. 2020;8(1):28–45.

    Article  Google Scholar 

  28. Zhang L, Li X, Ma B, Gao Q, Du H, Han Y, et al. The Tartary buckwheat genome provides insights into Rutin biosynthesis and abiotic stress tolerance. Mol Plant. 2017;10(9):1224–37.

    Article  CAS  PubMed  Google Scholar 

  29. Schneeberger K, Ossowski S, Ott F, Klein JD, Wang X, Lanz C, et al. Reference-guided assembly of four diverse Arabidopsis thaliana genomes. Proc Natl Acad Sci U S A. 2011;108(25):10249–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002;296(5565):79–92.

    Article  CAS  PubMed  Google Scholar 

  31. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: complexity, diversity, and dynamics. Science. 2009;326(5956):1112–5.

    Article  CAS  PubMed  Google Scholar 

  32. Moore RC, Purugganan MD. The early stages of duplicate gene evolution. Proc Natl Acad Sci U S A. 2003;100(26):15682–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Lopez-Maury L, Marguerat S, Bahler J. Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat Rev Genet. 2008;9(8):583–93.

    Article  CAS  PubMed  Google Scholar 

  34. Sun W, Jin X, Ma Z, Chen H, Liu M. Basic helix-loop-helix (bHLH) gene family in Tartary buckwheat (Fagopyrum tataricum): genome-wide identification, phylogeny, evolutionary expansion and expression analyses. Int J Biol Macromol. 2020;155:1478–90.

    Article  CAS  PubMed  Google Scholar 

  35. Liu M, Wang X, Sun W, Ma Z, Zheng T, Huang L, et al. Genome-wide investigation of the ZF-HD gene family in Tartary buckwheat (Fagopyrum tataricum). BMC Plant Biol. 2019;19(1):248.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Xu G, Guo C, Shan H, Kong H. Divergence of duplicate genes in exon-intron structure. Proc Natl Acad Sci U S A. 2012;109(4):1187–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Mizzotti C, Rotasperti L, Moretto M, Tadini L, Resentini F, Galliani BM, et al. Time-course transcriptome analysis of Arabidopsis siliques discloses genes essential for fruit development and maturation. Plant Physiol. 2018;178(3):1249–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yamada M, Han X, Benfey PN. RGF1 controls root meristem size through ROS signalling. Nature. 2020;577(7788):85–8.

    Article  CAS  PubMed  Google Scholar 

  39. Xia XJ, Zhou YH, Shi K, Zhou J, Foyer CH, Yu JQ. Interplay between reactive oxygen species and hormones in the control of plant development and stress tolerance. J Exp Bot. 2015;66(10):2839–56.

    Article  CAS  PubMed  Google Scholar 

  40. Peleg Z, Blumwald E. Hormone balance and abiotic stress tolerance in crop plants. Curr Opin Plant Biol. 2011;14(3):290–5.

    Article  CAS  PubMed  Google Scholar 

  41. Zhang K, Lan Y, Shi Y, Gao Y, Wu M, Xu Y, Xiang Y. Systematic analysis and functional characterization of the PLATZ transcription factors in Moso bamboo (Phyllostachys edulis). J Plant Growth Regul 2022. https://doi.org/https://doi.org/10.1007/s00344-021-10541-w.

  42. Yu X, Zhang W, Zhang Y, Zhang X, Lang D, Zhang X. The roles of methyl jasmonate to stress in plants. Funct Plant Biol. 2019;46(3):197–212.

    Article  CAS  PubMed  Google Scholar 

  43. Ramirez V, Coego A, Lopez A, Agorio A, Flors V, Vera P. Drought tolerance in Arabidopsis is controlled by the OCP3 disease resistance regulator. Plant J. 2009;58(4):578–91.

    Article  CAS  PubMed  Google Scholar 

  44. Lang D, Yu X, Jia X, Li Z, Zhang X. Methyl jasmonate improves metabolism and growth of NaCl-stressed Glycyrrhiza uralensis seedlings. Sci Hortic. 2020;266:109287.

    Article  CAS  Google Scholar 

  45. Wu H, Wu X, Li Z, Duan L, Zhang M. Physiological evaluation of drought stress tolerance and recovery in cauliflower (Brassica oleracea L.) seedlings treated with methyl Jasmonate and Coronatine. J Plant Growth Regul. 2012;31:113–23.

    Article  CAS  Google Scholar 

  46. Jin P, Zhu H, Wang J, Chen J, Wang X, Zheng Y. Effect of methyl jasmonate on energy metabolism in peach fruit during chilling stress. J Sci Food Agric. 2013;93(8):1827–32.

    Article  CAS  PubMed  Google Scholar 

  47. Yang C, Wang D, Zhang C, Ye M, Kong N, Ma H, et al. Comprehensive analysis and expression profiling of PIN, AUX/LAX, and ABCB auxin transporter gene families in Solanum tuberosum under Phytohormone stimuli and abiotic stresses. Biology. 2021;10(2):127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Weirauch MT, Hughes TR. Conserved expression without conserved regulatory sequence: the more things change, the more they stay the same. Trends Genet. 2010;26(2):66–74.

    Article  CAS  PubMed  Google Scholar 

  49. Li J, Zhang Y, Xu L, Wang C, Luo Y, Feng S, et al. Genome-wide identification of DNA binding with one finger (Dof) gene family in Tartary buckwheat (Fagopyrum tataricum) and analysis of its expression pattern after exogenous hormone stimulation. Biology. 2022;11(2):173.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7(10):e1002195.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zheng Y, Jiao C, Sun H, Rosli HG, Pombo MA, Zhang P, et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol Plant. 2016;9(12):1667–70.

    Article  CAS  PubMed  Google Scholar 

  52. Liu M, Ma Z, Wang A, Zheng T, Huang L, Sun W, et al. Genome-wide investigation of the auxin response factor gene family in Tartary buckwheat (Fagopyrum tataricum). Int J Mol Sci. 2018;19(11):3526.

    Article  PubMed Central  Google Scholar 

  53. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, 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 

  54. Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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 

  56. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Deng RY, Zhao HX, Xiao YH, Huang YJ, Yao PF, Lei YL, et al. Cloning, characterization, and expression analysis of eight stress-related NAC genes in Tartary buckwheat. Crop Sci. 2019;59(1):266–79.

    Article  CAS  Google Scholar 

  58. Li C, Zhao H, Li M, Yao P, Li Q, Zhao X, et al. Validation of reference genes for gene expression studies in tartary buckwheat (Fagopyrum tataricum Gaertn.) using quantitative real-time PCR. PeerJ. 2019;7:e6522.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  60. Yu CS, Lin CJ, Hwang JK. Predicting subcellular localization of proteins for gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004;13(5):1402–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Chou KC, Shen HB. Plant-mPLoc: a top-down strategy to augment the power for predicting plant protein subcellular localization. PLoS One. 2010;5(6):e11335.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Field A. Discovering statistics using IBM SPSS statistics. sage. 2013;

Download references

Acknowledgements

We thank all colleagues in our laboratory for providing useful discussions and technical assistance. We are very grateful to the editor and reviewers for critically evaluating the manuscript and providing constructive comments for its improvement.

Funding

This work was supported by the National Natural Science Foundation of China (31671631, 31501365); National Key R&D Program of China (2019YFD1000700, 2019YFD1000702); Shaanxi Province Modern Crops Seed Industry Project (20171010000004); and Technical System Project of Minor Grain Crops Industry in Shaanxi Province (NYKJ-2018-YL19).

Author information

Authors and Affiliations

Authors

Contributions

JL planned and designed the research, analyzed data, and wrote the original manuscript. SF analyzed data, and reviewed and edited the manuscript. YZ, LX, and YL performed the experiments. YY and QY reviewed and edited the manuscript. BF supervised the research. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Baili Feng.

Ethics declarations

Ethics approval and consent to participate

This article does not include any studies involving human participants or animals performed by the authors. These methods were carried out in accordance with the relevant guidelines and regulations. All experimental protocols were approved by Northwest A & F University.

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: Table S1.

CDS and protein sequences of FtPLATZs identified in this study.

Additional file 2: Table S2.

Putative motifs identified in FtPLATZ proteins by MEME.

Additional file 3: Table S3.

Orthologous gene pairs between Tartary buckwheat and other six representative species.

Additional file 4: Table S4.

Raw data of the expression profiles of FtPLATZ genes in different tissues and in different grain developmental stages of Tartary buckwheat analyzed by qRT-PCR.

Additional file 5: Table S5.

Cis-acting elements in the promoter regions of FtPLATZs.

Additional file 6: Table S6.

Raw data of the expression profiles of FtPLATZ genes in response to different exogenous hormone treatments analyzed by qRT-PCR.

Additional file 7: Table S7.

Primers of FtPLATZ genes used in this study.

Additional file 8: Fig. S1.

Multiple sequence alignment of PLATZ proteins in Tartary buckwheat.

Additional file 9: Fig. S2.

Picture of 21-day-old Tartary buckwheat seedlings treated with different exogenous hormones.

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

Li, J., Feng, S., Zhang, Y. et al. Genome-wide identification and expression analysis of the plant-specific PLATZ gene family in Tartary buckwheat (Fagopyrum tataricum). BMC Plant Biol 22, 160 (2022). https://doi.org/10.1186/s12870-022-03546-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-022-03546-4

Keywords