Recent Duplications Dominate VQ and WRKY Gene Expansions in Six <i>Prunus</i> Species

Zhong, Yan; Wang, Ping; Zhang, Xiaohui; Cheng, Zong-Ming

doi:https://doi.org/10.1155/2021/4066394

International Journal of Genomics

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 4066394 | https://doi.org/10.1155/2021/4066394

Recent Duplications Dominate VQ and WRKY Gene Expansions in Six Prunus Species

Yan Zhong,¹Ping Wang,¹Xiaohui Zhang,²and Zong-Ming Cheng¹

Academic Editor: Olasz Ferenc György

Received09 Apr 2021

Revised09 Jul 2021

Accepted19 Nov 2021

Published17 Dec 2021

Abstract

Genes encoding VQ motif-containing (VQ) transcriptional regulators and WRKY transcription factors can participate separately or jointly in plant growth, development, and abiotic and biotic stress responses. In this study, 222 VQ and 645 WRKY genes were identified in six Prunus species. Based on phylogenetic tree topologies, the VQ and WRKY genes were classified into 13 and 32 clades, respectively. Therefore, at least 13 VQ gene copies and 32 WRKY gene copies were present in the genome of the common ancestor of the six Prunus species. Similar small Ks value peaks for the VQ and WRKY genes suggest that the two gene families underwent recent duplications in the six studied species. The majority of the Ka/Ks ratios were less than 1, implying that most of the VQ and WRKY genes had undergone purifying selection. Pi values were significantly higher in the VQ genes than in the WRKY genes, and the VQ genes therefore exhibited greater nucleotide diversity in the six species. Forty-one of the Prunus VQ genes were predicted to interact with 44 of the WRKY genes, and the expression levels of some predicted VQ-WRKY interacting pairs were significantly correlated. Differential expression patterns of the VQ and WRKY genes suggested that some might be involved in regulating aphid resistance in P. persica and fruit development in P. avium.

1. Introduction

During their growth and development, plants are constantly threatened by an array of adverse environmental conditions. These include abiotic stresses such as drought, heat, and salt, and biotic stresses in the form of various pathogens and insects [1–4]. Plants have numerous transcriptional regulators and transcription factors (TFs) that help them cope with such ecological stresses. These regulators and TFs also have important functions in different aspects of physiological metabolism [5, 6].

VQ motif-containing (VQ) genes have been widely identified in different plants. For example, 34 VQ family members have been identified in Arabidopsis thaliana [7], 75 in soybean (Glycine max) [8], 26 in tomato (Solanum lycopersicum) [9], 49 in apple [10], and 59 in tobacco (Nicotiana tabacum L.) [11]. VQ genes have been found to play significant roles in plant growth, development, seed germination, and resistance to abiotic and biotic stresses [12–16]. The first identified VQ gene, sigma factor binding protein1 (SIB1), was found to be involved in resistance to Botrytis cinerea in A. thaliana [12, 13], and AtVQ14 was shown to participate in the HAIKU (IKU) pathway, regulating endosperm development and seed size [14].

The WRKY gene family is one of the largest TF superfamilies in eukaryotes [17, 18], and WRKY genes are generally more numerous than VQ genes in various species. There are 56 known WRKY genes in rose (Rosa chinensis) [19], 94 in sorghum (Sorghum bicolor [L.] Moench) [20], 104 in poplar (Populus trichocarpa) [21], and 112 in cotton (Gossypium raimondii) [22]. WRKY genes are associated with seed dormancy, flowering, lignin biosynthesis, and leaf senescence in rice (Oryza sativa), Arabidopsis, Chinese flowering cabbage (Brassica rapa L. ssp. pekinensis), poplar (P. trichocarpa), and wheat (Triticum aestivum L.) [23–28]. WRKY genes also participate extensively in plant responses to abiotic and biotic stresses [29, 30]. VvWRKY30 and GmWRKY12 have been shown to enhance the salt tolerance of grapevine (Vitis vinifera) and soybean (G. max), respectively [31, 32], and rice and Arabidopsis WRKY genes have been shown to participate in the plant immune response to pathogens [13, 33].

Interaction models of VQ and WRKY proteins have revealed that their interactions play important roles in plant immunity, growth, and development [5, 34–41]. For example, the interaction between AtVQ10 and AtWRKY8 reduces the damage caused by B. cinerea infection in Arabidopsis [42], and the interaction between MaVQ5 and MaWRKY26 enhances the survival of banana (Musa acuminata) under cold stress [43]. The interaction of AtVQ20 with AtWRKY2 and AtWRKY34 mediates male gametophytic functions and pollen development in Arabidopsis [44, 45].

Prunus species are important members of the Rosaceae family and include various ornamental and economically important fruit trees such as P. yedoensis, P. persica, and P. dulcis [46–51]. Although the different species show a variety of karyotypes in nature [52], Prunus genomes are usually very similar in size and basic synteny [53, 54]. The chromosome numbers of the genus Prunus are , and Prunus species evolved from a common ancestor with [55, 56]. The karyotypes of Prunus plants remain highly conserved, implying close relationships among these species [54, 55]. However, no genome-wide analysis has been performed to investigate the evolutionary patterns of VQ and WRKY genes and their interactions in closely related Prunus species. In this study, we identified VQ and WRKY genes in the genomes of P. yedoensis, P. domestica, P. avium, P. dulcis, P. persica, and P. yedoensis var. nudiflora. We then analyzed the phylogenetic relationships, evolutionary patterns, predicted interaction networks, and expression profiles of the VQ and WRKY genes from these species.

2. Materials and Methods

2.1. The Identification of VQ and WRKY Genes

The whole genome sequences of P. yedoensis, P. domestica, P. avium, P. dulcis, P. persica, and P. yedoensis var. nudiflora were downloaded from the Genome Database for Rosaceae (GDR; https://www.rosaceae.org/). All proteins in the genomes of each species were predicted using InterProScan with default parameters. Genes encoding a VQ or WRKY domain were identified as VQ or WRKY genes based on the InterProScan predictions. For genes with multiple transcripts, only the first transcript was retained for further analysis. The amino acid sequences of the VQ domain (PF05678) and the WRKY domain (PF03106) were obtained from the Pfam database (http://pfam.xfam.org/) and used as query sequences in TBLASTN searches against all the nucleotide coding sequences (CDSs) in the six genomes ( value ≤ 10⁻⁴) [57]. The BLAST hit sequences were verified by Pfam analysis, and hits that encoded VQ or WRKY domains were considered to be members of the VQ or WRKY gene family. Finally, the VQ and WRKY candidate genes were cross-verified by InterProScan and BLAST analysis. The sequences of VQ and WRKY genes from Fragaria vesca were obtained from previous studies [58, 59].

2.2. Construction of Phylogenetic Trees and Classification of Clades

The CDSs of the two gene families from the six Prunus species and the outgroup F. vesca were converted into amino acid sequences for alignment, and then, the amino acid alignments were converted back to nucleotide alignments to construct the phylogenetic trees. The VQ and WRKY phylogenetic trees were constructed in MEGA 7 using the neighbor-joining (NJ) method with pairwise deletion, and bootstrap values were obtained from 1000 replicates [60]. The trees were divided into various clades using two criteria: (i) each clade should contain one or more genes from F. vesca and (ii) each clade should contain genes from five or six Prunus species.

2.3. Genetic Parameters of VQ and WRKY Genes

Alignments of the CDSs from each clade in the two phylogenetic trees were obtained using ClustalW2 based on the alignment of amino acid sequences from the six Prunus species [61]. MEGA 7 was then used to calculate the synonymous substitution rate (Ks), nonsynonymous substitution rate (Ka), Ka/Ks ratio, and genetic diversity (Pi) values in each clade of the two gene families [60]. Due to Ks saturation, Ks values larger than 1 were discarded in the subsequent analysis. The duplication times () of the VQ and WRKY genes were estimated based on a mutation rate of point mutations per site per generation and 3 years/generation in peach: [62]. At the same time, the frequency of sequence exchanges in VQ and WRKY genes was examined using the GENECONV program with default parameters [63].

2.4. Prediction of Interactions between VQ and WRKY Genes

The VQ and WRKY CDSs from A. thaliana were downloaded from TAIR10 (http://www.Arabidopsis.org) [64]. These Arabidopsis genes were used as reference genes to uncover potential interactions between VQ and WRKY genes in the six Prunus species. The VQ and WRKY genes were numbered based on their gene ID orders (Tables S1 and S2). The CDSs of the VQ and WRKY genes from A. thaliana and the six studied species were then aligned, and 12 phylogenetic trees were constructed using the methods described above. Prunus and A. thaliana genes located in the same clade with bootstrap were considered to be homologs based on the phylogenetic tree topologies. Interaction relationships between the VQ and WRKY genes in the six Prunus species could then be predicted based on those of their VQ-WRKY homologs in A. thaliana [6–9, 65, 66]. The predicted interaction networks were then visualized using Cytoscape 3.7.1 (https://cytoscape.org/) [67].

2.5. Expression Patterns of VQ and WRKY Genes in P. persica and P. avium

RNA sequencing (RNA-seq) data for two P. persica peach lines (susceptible S38 and resistant R36) infested with the green peach aphid (GPA) from 3 h to 72 h were obtained from a previous study [68]. We then analyzed these data further to examine the expression of VQ and WRKY genes. Differences in VQ and WRKY gene expression were identified based on thresholds of and false discovery rate (FDR) adjusted value < 0.05. Expression heatmaps were created for VQ and WRKY genes that were differentially expressed between the R36 and S38 lines using the pheatmap package in RStudio (https://www.rstudio.com/) based on FPKM values (Fragments per Kilobase of transcript per Million mapped reads). The FPKM values of VQ and WRKY genes in sweet cherry (P. avium) from 3 to 94 days after full bloom (DAFB) were also used to create heatmaps using the same methods [69]. Pearson’s correlation coefficients (PCCs) between FPKM values of VQ and WRKY genes with predicted interaction relationships were calculated using SPSS version 20.0 (IBM Corp., Armonk, NY, USA).

3. Results

3.1. VQ and WRKY Genes in the Six Prunus Species

A total of 222 VQ genes were identified in the six Prunus species: 55 in P. yedoensis, 70 in P. domestica, 25 in P. avium, 23 in P. dulcis, 26 in P. persica, and 23 in P. yedoensis var. nudiflora (Table 1). P. domestica contained the largest number of VQ genes, followed by P. yedoensis. The other four Prunus species had lower and similar VQ gene numbers. Because the P. yedoensis var. nudiflora genome assembly is approximately half the size of the P. yedoensis assembly, it contained approximately half as many VQ genes.

The total number of WRKY genes identified in the six Prunus species (645) was greater than that of VQ genes (Table 1). The number of WRKY genes in each species was also greater than the number of VQ genes. Unsurprisingly, as with the VQ genes, the largest number of WRKY genes (262) was found in P. domestica, and this number was twice that found in P. yedoensis (139). Very similar WRKY gene numbers were found in P. avium (53), P. dulcis (56), and P. persica (58). However, the number of WRKY genes in P. yedoensis var. nudiflora (77) was slightly higher than that in P. avium, P. dulcis, and P. persica.

The average CDS lengths of the VQ and WRKY genes in the six Prunus species were irregular (Table 1), and the average CDS lengths were greater for WRKY genes than for VQ genes in each species.

3.2. Evolutionary Events Associated with VQ and WRKY Genes in the Six Prunus Species

The phylogenetic tree of VQ genes could be divided into 13 clades based on the two criteria described in Materials and Methods (Figure 1), indicating that no fewer than 13 ancient VQ copies existed in the genome of the six species’ common ancestor. Gene duplication and loss events were detected in each VQ clade during the evolution of the six Prunus species. Among all 222 VQ genes in the six species, approximately 60.81% (135/222) and 2.25% (5/222) had undergone gene duplication or gene loss, respectively (Table 2). The two species with the most VQ genes, P. domestica and P. yedoensis, also displayed the highest frequencies of gene duplication. Sixty-six VQ genes in P. domestica and 51 VQ genes in P. yedoensis were associated with 55 and 39 duplication events, respectively. These numbers of genes and duplication events were more than two times greater than those in the other four species. In general, there were very few VQ gene loss events in the six species. Two loss events occurred in P. domestica and P. yedoensis var. nudiflora, and one occurred in P. yedoensis. No loss events occurred in the remaining three species.

A phylogenetic tree was constructed using 645 WRKY genes from the six Prunus species and 61 WRKY genes from F. vesca, and the genes were classified into 32 clades (Figure 2). At least 32 ancient WRKY genes were therefore present in the ancestral species before the differentiation of the six extant Prunus species. Among all the WRKY genes, 66.67% (430/645) and 0.47% (3/645) were associated with gene duplication and loss events, respectively (Table 2). P. domestica was again associated with the most WRKY gene duplication events (223) and the largest number of WRKY genes involved in duplication (254). P. domestica was followed by P. yedoensis, with 103 duplication events and 135 duplication-related genes; these values were much lower than those of P. domestica. In the other four species, the number of WRKY gene duplication events ranged from 17 to 39, and the numbers of associated genes ranged from 29 to 60. Only two WRKY gene loss events were observed in P. avium, and one was observed in P. yedoensis var. nudiflora. No loss events were observed among the WRKY genes of the other species.

3.3. Duplication Times of VQ and WRKY Genes in the Six Prunus Species

The distribution of paralog Ks values revealed the duplication times and scales of the VQ and WRKY gene families. For the VQ genes, the Ks values were not continuously distributed. Most Ks values ranged from 0 to 0.5, and almost none were found between 0.5 and 1.0 (Figure 3(a)). This result indicated that almost no VQ gene duplication events occurred at a relatively ancient stage (0.5–1.0); instead, their duplications occurred relatively recently (0–0.5). In particular, the distinct Ks peak from 0 to 0.1 clearly demonstrated that the vast majority of VQ genes in the six species arose from recent duplications. To further investigate VQ gene duplication times, Ks values in the range of 0 to 0.01 were refined into ten smaller units (Figure 3(b)). At this smaller scale, the Ks values were distributed continuously and showed peaks at 0.01–0.02 and 0–0.01. This result indicated that very recent duplications have shaped the VQ genes in the six Prunus species.

(a)

(b)

(c)

(d)

Figure 3

The frequency of Ks values of paralogs and orthologs of VQ genes and WRKY genes in the six Prunus species. The -axis represents the Ks range, and the -axis represents the number of gene pairs. (a) The distribution of Ks values of VQ genes within the range of 0–1.0. (b) The distribution of Ks values of VQ genes within the range of 0–0.1. (c) The distribution of Ks values of WRKY genes within the range of 0–1.0. (d) The distribution of Ks values of WRKY genes within the range of 0–0.1.

The Ks values of paralogous WRKY genes were distributed within the range of 0–0.9 (Figure 3(c)), indicating that duplication events have occurred continuously over the course of WRKY gene evolution. The Ks peak of the WRKY genes was also detected between 0 and 0.1, similar to that of the VQ genes. This indicated that large-scale WRKY gene duplications occurred during the same recent period as the VQ gene duplications. Within the more detailed range of 0 to 0.1, the Ks values peaked at 0–0.01 and then gradually declined from 0.01 to 0.1 (Figure 3(d)). These results showed that the evolution of the WRKY genes in the six species was also driven mainly by very recent duplications. The similar duplication times of the dominant gene expansions in the two families may be indicative of close evolutionary relationships between some VQ and WRKY genes in the six species.

3.4. Evolutionary Patterns of VQ and WRKY Genes in the Six Prunus Species

In general, the Ka/Ks ratios were greater for the VQ genes than for the WRKY genes. This suggests that fewer functional restrictions operate on the VQ genes than on the WRKY genes (Figure 4(a)). In the VQ gene family, the Ka/Ks ratios for both paralogs and orthologs were generally less than 1 (Figure 4(b)). Among all pairs of VQ genes in the six Prunus species, more than 90% (1201/1307) had undergone purifying selection. The remaining gene pairs had Ka/Ks ratio values greater than 1, indicating that a minority of the VQ genes were under positive selection. There was also a difference between the Ka/Ks ratios of paralogs and orthologs: orthologs had significantly higher Ka/Ks ratios (-test, ). This suggested that paralogs were subjected to more functional restrictions than orthologs among the six species. The Ka/Ks ratios of the WRKY gene pairs were similar to those of the VQ gene pairs; most were less than 1, but approximately 10% (673/6407) were greater than 1. However, the median, mean, and quartiles of the Ka/Ks ratios were slightly higher for the paralogs than for the orthologs (Figure 4(b)).

(a)

(b)

(c)

(d)

The Pi values were significantly greater for the VQ genes than for the WRKY genes (-test, ), and the VQ genes therefore displayed greater nucleotide diversity (Figure 4(c)). Among the WRKY genes, the Pi values were significantly higher for paralogs than for orthologs (-test, ; Figure 4(d)).

The number of sequence exchange events was more than four times greater between the WRKY genes (1039) than between the VQ genes (220), a result that reflects the larger numbers of WRKY genes in the six species (Table S3). Significantly more frequent sequence exchanges were detected between orthologs than between paralogs for both VQ (-test, ) and WRKY genes (-test, ).

3.5. Interactions between the VQ and WRKY Genes

Based on the topologies of the phylogenetic trees (Figures S1 and S2), AtVQ and AtWRKY homologs were identified in the six Prunus species (Figure S3). By analogy to the interactions of their Arabidopsis homologs, 41 Prunus VQ genes were predicted to interact with 44 WRKY genes: 8 PyVQ genes with 4 PyWRKY genes in P. yedoensis, 10 PgVQs with 9 PgWRKYs in P. domestica, 3 PvVQs with 6 PvWRKYs in P. avium, 4 PaVQs with 9 PaWRKYs in P. dulcis, 10 PpVQs with 7 PpWRKYs in P. persica, and 6 PcVQs with 9 PcWRKYs in P. yedoensis var. nudiflora. One-to-one, one-to-many, and many-to-many interaction relationships were all predicted between the VQ and WRKY genes of the studied Prunus species. The single predicted one-to-one interaction relationship was PvVQ5-PvWRKY38 in P. avium. PpVQ1-PpWRKY3/PpWRKY31/PpWRKY37 in P. persica and PcVQ10/PcVQ17/PcVQ19-PcWRKY14 in P. yedoensis var. nudiflora were good examples of one-to-many interaction relationships. All the predicted VQ and WRKY genes in P. yedoensis and P. domestica exhibited many-to-many interaction relationships.

3.6. Expression Profiles of PpVQ and PpWRKY Genes in Peach during Aphid Infestation

Eleven PpVQ and 24 PpWRKY genes were significantly differentially expressed between the susceptible peach line, S38, and the resistant line, R36, after GPA infestation for 0, 3, 6, 9, 12, 24, 48, and 72 h (Figure 5).

(a)

(b)

(c)

(d)

The expression levels of the differentially expressed VQ genes were significantly lower in S38 than in R36 (-test, ), demonstrating that these VQ genes responded more strongly to GPA infestation in the resistant line than in the susceptible line. In particular, two genes (PpVQ22 and PpVQ24) had higher expression levels than the other genes in the two peach lines. They had high expression levels at 0 h, then displayed fluctuating upregulation at subsequent time points in S38 and R36. These patterns indicated that PpVQ22 and PpVQ24 expression continuously fluctuated in both lines during GPA infestation. In the R36 line, PpVQ8 and PpVQ9 also had relatively high expression levels that displayed wave-like oscillations.

The expression levels of differentially expressed WRKY genes were likewise higher in R36 than in S38. PpWRKY19 had the highest expression level of the differentially expressed WRKY genes, and its expression also exhibited waves of upregulation from 0 to 72 h in both lines. Similar tendencies were observed for PpWRKY46 and PpWRKY53 in both lines and for PpWRKY7, PpWRKY8, PpWRKY9, PpWRKY16, PpWRKY41, PpWRKY44, and PpWRKY50 in R36 only.

Although some of the same genes exhibited similar expression patterns in both lines, more VQ and WRKY genes had higher expression levels in the resistant line and may have contributed to its ability to defend against GPA infestation.

3.7. Expression Patterns of PvVQ and PvWRKY Genes during Sweet Cherry Exocarp Development

Thirteen PvVQ genes and 32 PvWRKY genes were differentially expressed during the development of sweet cherry fruit exocarps from 3 to 94 DAFB. The experimental period encompassed three developmental stages: (i) 3–30 DAFB (cell division and expansion), (ii) 31–45 DAFB (seed development), and (iii) 52–94 DAFB (rapid cell expansion; Figure 6). The differentially expressed VQ genes exhibited four expression patterns from full bloom to fruit ripening. The first pattern, which was displayed by genes such as PvVQ17, was characterized by higher expression at 3 DAFB. The second pattern, shown by PvVQ4 and PvVQ20, was characterized by similar expression at all three stages. The third pattern involved relatively high expression at 31 and 52 DAFB but lower expression at other time points. PvVQ7, PvVQ16, and PvVQ21 displayed this pattern and may therefore be involved in the early phases of seed development and rapid cell expansion. The fourth pattern was shown by PvVQ12, PvVQ19, and PvVQ25 and was characterized by relatively low expression at stages I and II but higher expression in the later period of stage III.

(a)

(b)

The differentially expressed WRKY genes also displayed four expression patterns. PvWRKY8, PvWRKY18, and PvWRKY34 displayed the first pattern, with higher expression in the early phase of stage I. PvWRKY3, PvWRKY19, and PvWRKY20 represented the second pattern, showing little variation in expression among the three stages. WRKY genes with distinctly higher expression in the early phase of stage II (e.g., PvWRKY25 and PvWRKY35) or in the later phase of stage III (e.g., PvWRKY13 and PvWRKY16) were assigned to the third or fourth pattern, respectively. These results indicated that some PvVQ and PvWRKY genes were differentially expressed during fruit development and maturation in sweet cherry.

3.8. Correlations between Related VQ and WRKY Genes

Using FPKM values of the VQ and WRKY genes in the two peach lines, PCC values were calculated between predicted VQ-WRKY interaction partners (Table 3). Ten VQ-WRKY gene pairs had significant correlations () or very significant correlations () in expression in the resistant peach line, R36. Nine of the significantly correlated gene pairs had significant positive correlations, but the remaining pair (PpVQ20-PpWRKY37) had a highly significant negative correlation. Only one sweet cherry gene pair (PvVQ22-PvWRKY48) showed a significant correlation in expression. These significantly correlated gene pairs may participate in responses to pathogen infection or in other development-related physiological activities.

4. Discussion

4.1. Gene Duplications Have Shaped the VQ and WRKY Families in Six Prunus Species

Gene duplication is one of the mechanisms by which genetic materials are produced for plant evolution [70]. Numerous gene families have been shaped by gene duplication in plant genomes. Examples include the VQ gene family in soybean (G. max) and pear (Pyrus bretschneideri) [8, 71] and the WRKY gene family in sorghum (S. bicolor) and potato (Solanum tuberosum) [20, 72].

In this study, at least 60% of the VQ and WRKY family members in six Prunus genomes were generated by gene duplication events. Gene duplication was therefore one of the most important genetic events promoting gene expansions during the evolution of these two gene families in the Prunus species. Although the VQ and WRKY genes of the six species originated from common ancestral VQ and WRKY copies, the different scales of the duplication events they experienced led to different numbers of VQ and WRKY genes in the six genomes. More gene duplication events were detected for the VQ and WRKY genes of P. yedoensis and P. domestica than for those of the other four species; this was the principal cause of the greater gene numbers in P. yedoensis and P. domestica (Table 2).

4.2. Relatively Recent Duplications Were the Main Drivers of VQ and WRKY Gene Expansions in the Six Prunus Species

Recent duplications have driven the expansion and evolution of plant gene families [73], including the NBS-LRR gene families in grapevine (V. vinifera) and poplar (P. trichocarpa) [74]. Most soybean VQ genes have undergone recent segmental duplications [8]. In the present study, the majority of VQ genes in the six Prunus species were generated by recent duplication events, as evidenced by paralogous Ks values that peaked at 0.01–0.02 and 0–0.01 (Figure 3). The most recent duplication of the VQ genes appears to have occurred approximately <6.33 million years ago (MYA), based on the mutation rate in peach (Materials and Methods). This duplication therefore arose after the divergence of the Prunus species, which occurred between 36 and 44 MYA [47]. Furthermore, the Ks values of WRKY gene paralogs peaked at 0–0.01, indicating that the WRKY genes of the six species began to experience recent expansion events approximately <3.16 MYA. Thus, large-scale, recent gene duplications were the dominant force driving VQ and WRKY gene expansions among the six Prunus species. No whole-genome duplication (WGD) events have occurred recently in Prunus species such as P. persica [49, 75–77]. Therefore, expansions of the VQ and WRKY genes may be only loosely connected with WGDs in these Prunus species.

4.3. Different Evolutionary Patterns between VQ and WRKY Gene Families in the Six Prunus Species

The vast majority of VQ and WRKY genes have experienced purifying selection (Figures 4(a) and 4(b); ) to eliminate deleterious mutations and maintain gene functions during evolution [78]. These results are consistent with those for VQ genes from the genomes of three legumes and tobacco (N. tabacum L.) [8, 11, 79] and for the WRKY genes of Arabidopsis, cotton (G. raimondii), and rose (R. chinensis) [19, 22, 80]. Therefore, most members of the two gene families in the six species studied here had relatively low variation rates and conserved functions. A fraction of the VQ and WRKY genes were under positive selection (), which may have been related to their ability to help the plant cope with various environmental conditions [77]. Trends in Ka/Ks distribution were similar between the VQ and WRKY genes. However, Ka/Ks ratios, which are indicative of selection pressure, were higher in the VQ genes than in the WRKY genes. The VQ genes also had more significant nucleotide diversity than the WRKY genes (-test, ). VQ proteins interact with a variety of partners, and VQ genes therefore encode extremely diverse proteins with different primary structures [37]. Therefore, WRKY genes may show greater conservation and are more likely to be trapped by functional restrictions than VQ genes in the six Prunus genomes.

4.4. PpVQ and PpWRKY Genes Respond Together to Peach Aphid Infestation

VQ and WRKY genes have previously been reported to participate in plant disease response and defense [33, 36, 39, 42]. In the present study, 11 PpVQ and 24 PpWRKY genes were differentially expressed under GPA infestation in two peach lines (Figure 5). For example, PpVQ24 expression responded strongly to GPA infestation; PpVQ24 is a homolog of AtVQ16 and AtVQ23, which are associated with defense against B. cinerea in Arabidopsis [13]. Likewise, PpWRKY19, which displayed the highest expression levels after GPA infestation, is homologous to AtWRKY54 and AtWRKY70. These two Arabidopsis WRKY TFs negatively mediate responses to drought and necrotrophic pathogens [39, 81, 82]. These results demonstrate that PpVQ24 and PpWRKY19 may be involved in the GPA response in peach.

Overall, PpVQ and PpWRKY genes responded more strongly to GPA infestation in the resistant line, R36, than in the susceptible line, S38. This was consistent with the longer duration of defense in R36 than in S38 [68]. Stronger responses of VQ and WRKY genes to pathogen infection in resistant plants have been reported previously in other species. Examples include the GmVQ genes in soybean (G. max) infested with the common cutworm (Spodoptera litura Fabricius), OsVQ genes in rice (O. sativa) after infection with three pathogens, and WRKY genes in chrysanthemum (Chrysanthemum morifolium) after aphid infestation [15, 16, 35].

It is well-known that VQ proteins interact with WRKY TFs to participate jointly in plant defense against pathogens and insects. In peaches subjected to GPA infestation, the FPKM values of ten predicted interacting gene pairs (e.g., PpVQ8-PpWRKY37, PpVQ8-PpWRKY44, and PpVQ20-PpWRKY37) were significantly correlated. Interestingly, the Arabidopsis genes AtVQ16 (SIB2) and AtVQ23 (SIB1), which are homologs of PpVQ8, increase B. cinerea resistance by interacting with AtWRKY33, which is a PpWRKY44 homolog [13]. We speculate that PpVQ and PpWRKY genes with highly correlated expression patterns may interact to respond together to GPA infestation in peach, and this possibility awaits further research.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Yan Zhong and Ping Wang contributed equally to this work.

Acknowledgments

This study was supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions, the National Natural Science Foundation of China (31870204), and the high-performance computing platform of Bioinformatics Center, Nanjing Agricultural University.

Supplementary Materials

Supplementary 1. Table S1: the correspondence between gene IDs and VQ names in the six Prunus species.

Supplementary 2. Table S2: the correspondence between gene IDs and WRKY names in the six Prunus species.

Supplementary 3. Table S3: the sequence exchange events of VQ and WRKY genes in the six Prunus species.

Supplementary 4. Figure S1: phylogenetic trees of VQ genes from Arabidopsis and the six Prunus species. (a) P. yedoensis. (b) P. domestica. (c) P. avium. (d) P. dulcis. (e) P. persica. (f) P. yedoensis var. nudiflora. The square line indicates homologous VQ genes of A. thaliana and Prunus species in the same clade with .

Supplementary 5. Figure S2: phylogenetic trees of WRKY genes from Arabidopsis and the six Prunus species. (a) P. yedoensis. (b) P. domestica. (c) P. avium. (d) P. dulcis. (e) P. persica. (f) P. yedoensis var. nudiflora. The square line indicates homologous WRKY genes of A. thaliana and Prunus species in the same clade with .

Supplementary 6. Figure S3: predicted interaction networks of VQ and WRKY proteins from the six Prunus species. (a) P. yedoensis. (b) P. domestica. (c) P. avium. (d) P. dulcis. (e) P. persica. (f) P. yedoensis var. nudiflora.

References

A. Banerjee and A. Roychoudhury, “WRKY proteins: signaling and regulation of expression during abiotic stress responses,” ScientificWorldJournal, vol. 2015, article 807560, 2015.
View at: Google Scholar
Y. L. Bai, C. Kissoudis, Z. Yan, R. G. F. Visser, and G. van der Linden, “Plant behaviour under combined stress: tomato responses to combined salinity and pathogen stress,” Plant Journal, vol. 93, no. 4, pp. 781–793, 2018.
View at: Publisher Site | Google Scholar
J. J. Jiang, S. H. Ma, N. H. Ye, M. Jiang, J. Cao, and J. Zhang, “WRKY transcription factors in plant responses to stresses,” Journal of Integrative Plant Biology, vol. 59, no. 2, pp. 86–101, 2017.
View at: Publisher Site | Google Scholar
S. S. K. P. Vurukonda, S. Vardharajula, M. Shrivastava, and A. SkZ, “Enhancement of drought stress tolerance in crops by plant growth promoting rhizobacteria,” Microbiological Research, vol. 184, pp. 13–24, 2016.
View at: Publisher Site | Google Scholar
Y. Chi, Y. Yang, Y. Zhou et al., “Protein-protein interactions in the regulation of WRKY transcription factors,” Molecular Plant, vol. 6, no. 2, pp. 287–300, 2013.
View at: Publisher Site | Google Scholar
M. Wang, A. Vannozzi, G. Wang et al., “A comprehensive survey of the grapevine VQ gene family and its transcriptional correlation with WRKY proteins,” Frontiers in Plant Science, vol. 6, p. 417, 2015.
View at: Publisher Site | Google Scholar
Y. Cheng, Y. Zhou, Y. Yang et al., “Structural and functional analysis of VQ motif-containing proteins in Arabidopsis as interacting proteins of WRKY transcription factors,” Plant Physiology, vol. 159, no. 2, p. 810, 2012.
View at: Publisher Site | Google Scholar
Y. B. Wang, Z. F. Jiang, Z. X. Li et al., “Genome-wide identification and expression analysis of theVQgene family in soybean (glycine max),” PeerJ, vol. 7, article e7509, 2019.
View at: Publisher Site | Google Scholar
H. D. Ding, G. B. Yuan, S. R. Mo et al., “Genome-wide analysis of the plant-specific VQ motif-containing proteins in tomato (Solanum lycopersicum_ ) and characterization of SlVQ6 in thermotolerance,” Plant Physiology and Biochemistry, vol. 143, pp. 29–39, 2019.
View at: Publisher Site | Google Scholar
Q. L. Dong, S. Zhao, D. Y. Duan et al., “Structural and functional analyses of genes encoding VQ proteins in apple,” Plant Science, vol. 272, pp. 208–219, 2018.
View at: Publisher Site | Google Scholar
C. H. Liu, H. Liu, C. Y. Zhou, and M. P. Timko, “Genome-wide identification of the VQ protein gene family of tobacco (Nicotiana tabacum L.) and analysis of its expression in response to phytohormones and abiotic and biotic stresses,” Genes (Basel), vol. 11, no. 3, p. 284, 2020.
View at: Publisher Site | Google Scholar
K. Morikawa, T. Shiina, S. Murakami, and Y. Toyoshima, “Novel nuclear-encoded proteins interacting with a plastid sigma factor, Sig1, in Arabidopsis thaliana,” FEBS Letters, vol. 514, no. 2-3, pp. 300–304, 2002.
View at: Publisher Site | Google Scholar
Z. B. Lai, Y. Li, F. Wang et al., “Arabidopsis sigma factor binding proteins are activators of the WRKY33 transcription factor in plant defense,” Plant Cell, vol. 23, no. 10, pp. 3824–3841, 2011.
View at: Publisher Site | Google Scholar
A. H. Wang, D. Garcia, H. Y. Zhang et al., “The VQ motif protein IKU1 regulates endosperm growth and seed size in Arabidopsis,” Plant Journal, vol. 63, no. 4, pp. 670–679, 2010.
View at: Publisher Site | Google Scholar
X. Li, R. Qin, Q. Du et al., “Knockdown of GmVQ58 encoding a VQ motif-containing protein enhances soybean resistance to the common cutworm (Spodoptera litura Fabricius),” Journal of Experimental Botany, vol. 71, no. 10, pp. 3198–3210, 2020.
View at: Publisher Site | Google Scholar
N. Li, X. Li, J. Xiao, and S. Wang, “Comprehensive analysis of VQ motif-containing gene expression in rice defense responses to three pathogens,” Plant Cell Reports, vol. 33, no. 9, pp. 1493–1505, 2014.
View at: Publisher Site | Google Scholar
J. P. Jin, H. Zhang, L. Kong, G. Gao, and J. Luo, “PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors,” Nucleic Acids Research, vol. 42, no. D1, pp. D1182–D1187, 2014.
View at: Publisher Site | Google Scholar
Y. J. Zhang and L. J. Wang, “The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants,” BMC Evolutionary Biology, vol. 5, no. 1, p. 1, 2005.
View at: Publisher Site | Google Scholar
X. T. Liu, D. D. Li, S. Y. Zhang, Y. Xu, and Z. Zhang, “Genome-wide characterization of the rose (Rosa chinensis) WRKY family and role of RcWRKY41 in gray mold resistance,” BMC Plant Biology, vol. 19, no. 1, p. 522, 2019.
View at: Publisher Site | Google Scholar
E. H. Baillo, M. S. Hanif, Y. H. Guo, Z. Zhang, P. Xu, and S. A. Algam, “Genome-wide identification of WRKY transcription factor family members in sorghum (Sorghum bicolor(L.) moench),” PLoS One, vol. 15, no. 8, p. e0236651, 2020.
View at: Publisher Site | Google Scholar
H. S. He, Q. Dong, Y. H. Shao et al., “Genome-wide survey and characterization of the WRKY gene family in Populus trichocarpa,” Plant Cell Reports, vol. 31, no. 7, pp. 1199–1217, 2012.
View at: Publisher Site | Google Scholar
M. Q. Ding, J. D. Chen, Y. R. Jiang et al., “Genome-wide investigation and transcriptome analysis of the WRKY gene family in Gossypium,” Molecular Genetics and Genomics, vol. 290, no. 1, pp. 151–171, 2015.
View at: Publisher Site | Google Scholar
C. Zhou, Q. Lin, J. Lan et al., “WRKY transcription factor OsWRKY29 represses seed dormancy in rice by weakening abscisic acid response,” Frontiers in Plant Science, vol. 11, p. 691, 2020.
View at: Publisher Site | Google Scholar
W. Li, H. P. Wang, and D. Q. Yu, “_Arabidopsis_ WRKY Transcription Factors WRKY12 and WRKY13 Oppositely Regulate Flowering under Short-Day Conditions,” Molecular Plant, vol. 9, no. 11, pp. 1492–1503, 2016.
View at: Publisher Site | Google Scholar
L. Yang, X. Zhao, F. Yang, D. Fan, Y. Jiang, and K. Luo, “PtrWRKY19, a novel WRKY transcription factor, contributes to the regulation of pith secondary wall formation in _Populus trichocarpa,” Scientific Reports, vol. 6, no. 1, p. 18643, 2016.
View at: Publisher Site | Google Scholar
L. G. Chen, S. Y. Xiang, Y. L. Chen, D. Li, and D. Yu, “Arabidopsis WRKY45 Interacts with the DELLA Protein RGL1 to Positively Regulate Age-Triggered Leaf Senescence,” Molecular Plant, vol. 10, no. 9, pp. 1174–1189, 2017.
View at: Publisher Site | Google Scholar
Z. Q. Fan, X. L. Tan, W. Shan, J. F. Kuang, W. J. Lu, and J. Y. Chen, “BrWRKY65, a WRKY transcription factor, is involved in regulating three leaf senescence-associated genes in Chinese flowering cabbage,” International Journal of Molecular Sciences, vol. 18, no. 6, p. 1228, 2017.
View at: Publisher Site | Google Scholar
L. Zhao, W. Zhang, Q. Song et al., “A WRKY transcription factor, TaWRKY40-D, promotes leaf senescence associated with jasmonic acid and abscisic acid pathways in wheat,” Plant Biology, vol. 22, no. 6, pp. 1072–1085, 2020.
View at: Publisher Site | Google Scholar
R. Li, J. Zhang, J. C. Li et al., “Prioritizing plant defence over growth through WRKY regulation facilitates infestation by non-target herbivores,” eLife, vol. 4, article e04805, 2015.
View at: Publisher Site | Google Scholar
R. P. Birkenbihl, B. Kracher, A. Ross, K. Kramer, I. Finkemeier, and I. E. Somssich, “Principles and characteristics of the Arabidopsis WRKY regulatory network during early MAMP-triggered immunity,” Plant Journal, vol. 96, no. 3, pp. 487–502, 2018.
View at: Publisher Site | Google Scholar
D. Zhu, L. X. Hou, P. L. Xiao, Y. Guo, M. K. Deyholos, and X. Liu, “VvWRKY30_, a grape WRKY transcription factor, plays a positive regulatory role under salinity stress,” Plant Science, vol. 280, pp. 132–142, 2019.
View at: Publisher Site | Google Scholar
W. Y. Shi, Y. T. du, J. Ma et al., “The WRKY transcription factor GmWRKY12 confers drought and salt tolerance in soybean,” International Journal of Molecular Sciences, vol. 19, no. 12, p. 4087, 2018.
View at: Publisher Site | Google Scholar
H. Lee, J. Cha, C. Choi et al., “Rice WRKY11 plays a role in pathogen defense and drought tolerance,” Rice, vol. 11, no. 1, p. 5, 2018.
View at: Publisher Site | Google Scholar
D. M. Yao, C. Zou, Y. N. Shu, and S. S. Liu, “WRKY transcription factors in Nicotiana tabacum modulate plant immunity against whitefly via interacting with MAPK cascade pathways,” Insects, vol. 12, no. 1, p. 16, 2021.
View at: Publisher Site | Google Scholar
P. Li, A. Song, C. Gao et al., “The over-expression of a chrysanthemum _WRKY_ transcription factor enhances aphid resistance,” Plant Physiology and Biochemistry, vol. 95, pp. 26–34, 2015.
View at: Publisher Site | Google Scholar
H. Wang, Y. Hu, J. Pan, and D. Yu, “Arabidopsis VQ motif-containing proteins VQ12 and VQ29 negatively modulate basal defense against Botrytis cinerea,” Scientific Reports, vol. 5, no. 1, p. 14185, 2015.
View at: Publisher Site | Google Scholar
Y. Jing and R. Lin, “The VQ motif-containing protein family of plant-specific transcriptional regulators,” Plant Physiology, vol. 169, no. 1, pp. 371–378, 2015.
View at: Publisher Site | Google Scholar
P. Hu, W. Zhou, Z. W. Cheng, M. Fan, L. Wang, and D. Xie, “JAV1 Controls Jasmonate-Regulated Plant Defense,” Molecular Cell, vol. 50, no. 4, pp. 504–515, 2013.
View at: Publisher Site | Google Scholar
J. Li, R. Zhong, and E. T. Palva, “WRKY70 and its homolog WRKY54 negatively modulate the cell wall-associated defenses to necrotrophic pathogens in Arabidopsis,” PLoS One, vol. 12, no. 8, article e0183731, 2017.
View at: Publisher Site | Google Scholar
S. Besseau, J. Li, and E. T. Palva, “WRKY54 and WRKY70 co-operate as negative regulators of leaf senescence in Arabidopsis thaliana,” Journal of Experimental Botany, vol. 63, no. 7, pp. 2667–2679, 2012.
View at: Publisher Site | Google Scholar
Y. Y. Lei, Y. P. Sun, B. T. Wang et al., “Woodland strawberry WRKY71 acts as a promoter of flowering via a transcriptional regulatory cascade,” Horticulture Research, vol. 7, no. 1, p. 137, 2020.
View at: Publisher Site | Google Scholar
J. Q. Chen, H. P. Wang, Y. Li, J. Pan, Y. Hu, and D. Yu, “ArabidopsisVQ10 interacts with WRKY8 to modulate basal defense againstBotrytis cinerea,” Journal of Integrative Plant Biology, vol. 60, no. 10, pp. 956–969, 2018.
View at: Publisher Site | Google Scholar
Y. J. Ye, Y. Y. Xiao, Y. C. Han et al., “Banana fruit VQ motif-containing protein5 represses cold-responsive transcription factor MaWRKY26 involved in the regulation of JA biosynthetic genes,” Scientific Reports, vol. 6, no. 1, p. 23632, 2016.
View at: Publisher Site | Google Scholar
R. H. Lei, X. L. Li, Z. B. Ma, Y. Lv, Y. Hu, and D. Yu, “Arabidopsis WRKY2 and WRKY34 transcription factors interact with VQ20 protein to modulate pollen development and function,” Plant Journal, vol. 91, no. 6, pp. 962–976, 2017.
View at: Publisher Site | Google Scholar
R. Lei, Z. Ma, and D. Yu, “WRKY2/34-VQ20 modules in Arabidopsis thaliana negatively regulate expression of a trio of related MYB transcription factors during pollen development,” Frontiers in Plant Science, vol. 9, p. 331, 2018.
View at: Publisher Site | Google Scholar
K. Shirasawa, T. Esumi, H. Hirakawa et al., “Phased genome sequence of an interspecific hybrid flowering cherry, ‘Somei-Yoshino’ (Cerasus × yedoensis),” DNA Research, vol. 26, no. 5, pp. 379–389, 2019.
View at: Publisher Site | Google Scholar
S. Baek, K. Choi, G. B. Kim et al., “Draft genome sequence of wild Prunus yedoensis reveals massive inter-specific hybridization between sympatric flowering cherries,” Genome Biology, vol. 19, no. 1, p. 127, 2018.
View at: Publisher Site | Google Scholar
K. Shirasawa, K. Isuzugawa, M. Ikenaga et al., “The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding,” DNA Research, vol. 24, no. 5, pp. 499–508, 2017.
View at: Publisher Site | Google Scholar
The International Peach Genome Initiative, I. Verde, A. G. Abbott et al., “The high-quality draft genome of peach ( _Prunus persica_ ) identifies unique patterns of genetic diversity, domestication and genome evolution,” Nature Genetics, vol. 45, no. 5, pp. 487–494, 2013.
View at: Publisher Site | Google Scholar
T. Zhebentyayeva, V. Shankar, R. Scorza et al., “Genetic characterization of worldwide _Prunus domestica_ (plum) germplasm using sequence-based genotyping,” Horticulture Research, vol. 6, no. 1, p. 12, 2019.
View at: Publisher Site | Google Scholar
T. Alioto, K. G. Alexiou, A. Bardil et al., “Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence,” Plant Journal, vol. 101, no. 2, pp. 455–472, 2020.
View at: Publisher Site | Google Scholar
C. J. Ye, L. Stilgenbauer, A. Moy, G. Liu, and H. H. Heng, “What is karyotype coding and why is genomic topology important for cancer and evolution?” Frontiers in Genetics, vol. 10, p. 1082, 2019.
View at: Publisher Site | Google Scholar
S. Vilanova, D. J. Sargent, P. Arús, and A. Monfort, “Synteny conservation between two distantly-related Rosaceae genomes: Prunus (the stone fruits) and Fragaria (the strawberry),” BMC Plant Biology, vol. 8, no. 1, p. 67, 2008.
View at: Publisher Site | Google Scholar
S. Jung, D. Jiwan, I. Cho et al., “Synteny of Prunus and other model plant species,” BMC Genomics, vol. 10, no. 1, p. 76, 2009.
View at: Publisher Site | Google Scholar
E. Dirlewanger, E. Graziano, T. Joobeur et al., “Comparative mapping and marker-assisted selection in Rosaceae fruit crops,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 26, pp. 9891–9896, 2004.
View at: Publisher Site | Google Scholar
E. Illa, D. J. Sargent, E. L. Girona et al., “Comparative analysis of rosaceous genomes and the reconstruction of a putative ancestral genome for the family,” BMC Evolutionary Biology, vol. 11, no. 1, p. 9, 2011.
View at: Publisher Site | Google Scholar
Y. LI, Y. ZHONG, K. HUANG, and Z. M. CHENG, “Genomewide analysis of NBS-encoding genes in kiwi fruit (Actinidia chinensis),” Journal of Genetics, vol. 95, no. 4, pp. 997–1001, 2016.
View at: Publisher Site | Google Scholar
Y. Zhong, C. Guo, J. J. Chu, H. Liu, and Z. M. Cheng, “Microevolution of the VQ gene family in six species of Fragaria,” Genome, vol. 61, no. 1, pp. 49–57, 2018.
View at: Publisher Site | Google Scholar
W. Wei, Y. Hu, Y. T. Han, K. Zhang, F. L. Zhao, and J. Y. Feng, “The WRKY transcription factors in the diploid woodland strawberry _Fragaria vesca_ : Identification and expression analysis under biotic and abiotic stresses,” Plant Physiology and Biochemistry, vol. 105, pp. 129–144, 2016.
View at: Publisher Site | Google Scholar
S. Kumar, G. Stecher, and K. Tamura, “MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets,” Molecular Biology and Evolution, vol. 33, no. 7, pp. 1870–1874, 2016.
View at: Publisher Site | Google Scholar
M. A. Larkin, G. Blackshields, N. P. Brown et al., “Clustal W and clustal X version 2.0,” Bioinformatics, vol. 23, no. 21, pp. 2947-2948, 2007.
View at: Publisher Site | Google Scholar
Z. Xie, L. Wang, L. Wang et al., “Mutation rate analysis via parent-progeny sequencing of the perennial peach. I. A low rate in woody perennials and a higher mutagenicity in hybrids,” Proceedings of the Biological Sciences, vol. 283, no. 1841, p. 1016, 2016.
View at: Google Scholar
S. Sawyer, “Statistical tests for detecting gene conversion,” Molecular Biology and Evolution, vol. 6, no. 5, pp. 526–538, 1989.
View at: Google Scholar
P. Lamesch, T. Z. Berardini, D. H. Li et al., “The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools,” Nucleic Acids Research, vol. 40, no. D1, pp. D1202–D1210, 2012.
View at: Publisher Site | Google Scholar
J. Garrido-Gala, J. J. Higuera, J. Muñoz-Blanco, F. Amil-Ruiz, and J. L. Caballero, “The VQ motif-containing proteins in the diploid and octoploid strawberry,” Scientific Reports, vol. 9, no. 1, p. 4942, 2019.
View at: Publisher Site | Google Scholar
J. Guo, J. Chen, J. Yang, Y. Yu, Y. Yang, and W. Wang, “Identification, characterization and expression analysis of the VQ motif-containing gene family in tea plant (Camellia sinensis),” BMC Genomics, vol. 19, no. 1, p. 710, 2018.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
L. Niu, L. Pan, W. Zeng et al., “Dynamic transcriptomes of resistant and susceptible peach lines after infestation by green peach aphids (Myzus persicae Sülzer) reveal defence responses controlled by the Rm3 locus,” BMC Genomics, vol. 19, no. 1, p. 846, 2018.
View at: Publisher Site | Google Scholar
M. Alkio, U. Jonas, M. Declercq, S. van Nocker, and M. Knoche, “Transcriptional dynamics of the developing sweet cherry ( _Prunus avium_ L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes,” Horticulture Research, vol. 1, no. 1, p. 11, 2014.
View at: Publisher Site | Google Scholar
N. Panchy, M. Lehti-Shiu, and S. H. Shiu, “Evolution of gene duplication in plants,” Plant Physiology, vol. 171, no. 4, pp. 2294–2316, 2016.
View at: Publisher Site | Google Scholar
Y. Cao, D. Meng, M. Abdullah, Q. Jin, Y. Lin, and Y. Cai, “Genome wide identification, evolutionary, and expression analysis of VQ genes from two Pyrus species,” Genes (Basel), vol. 9, no. 4, p. 224, 2018.
View at: Publisher Site | Google Scholar
C. Zhang, D. D. Wang, C. H. Yang et al., “Genome-wide identification of the potato WRKY transcription factor family,” PLoS One, vol. 12, no. 7, article e0181573, 2017.
View at: Publisher Site | Google Scholar
C. S. Brandon, M. J. Greenwold, and J. L. Dudycha, “Ancient and recent duplications support functional diversity of daphnia opsins,” Journal of Molecular Evolution, vol. 84, no. 1, pp. 12–28, 2017.
View at: Publisher Site | Google Scholar
S. H. Yang, X. H. Zhang, J. X. Yue, D. Tian, and J. Q. Chen, “Recent duplications dominate NBS-encoding gene expansion in two woody species,” Molecular Genetics and Genomics, vol. 280, no. 3, pp. 187–198, 2008.
View at: Publisher Site | Google Scholar
Q. M. Chen, Q. H. Li, X. Qiao, H. Yin, and S. Zhang, “Genome-wide identification of lysin motif containing protein family genes in eight Rosaceae species, and expression analysis in response to pathogenic fungus Botryosphaeria dothidea in Chinese white pear,” BMC Genomics, vol. 21, no. 1, p. 612, 2020.
View at: Publisher Site | Google Scholar
Y. P. Cao, Y. H. Han, D. D. Meng et al., “Expansion and evolutionary patterns of GDSL-type esterases/lipases in Rosaceae genomes,” Functional & Integrative Genomics, vol. 18, no. 6, pp. 673–684, 2018.
View at: Publisher Site | Google Scholar
Y. Zhong, X. H. Zhang, Q. L. Shi, and Z. M. Cheng, “Adaptive evolution driving the young duplications in six Rosaceae species,” BMC Genomics, vol. 22, no. 1, p. 112, 2021.
View at: Publisher Site | Google Scholar
P. Soares, L. Ermini, N. Thomson et al., “Correcting for purifying selection: an improved human mitochondrial molecular clock,” American Journal of Human Genetics, vol. 84, no. 6, pp. 740–759, 2009.
View at: Publisher Site | Google Scholar
L. Ling, Y. Qu, J. T. Zhu, D. Wang, and C. Guo, “Genome-wide identification and expression analysis of theVQgene family inCicer arietinumandMedicago truncatula,” PeerJ, vol. 8, article e8471, 2020.
View at: Publisher Site | Google Scholar
Q. S. Wang, M. H. Wang, X. Z. Zhang, B. Hao, S. K. Kaushik, and Y. Pan, “WRKY gene family evolution in Arabidopsis thaliana,” Genetica, vol. 139, no. 8, pp. 973–983, 2011.
View at: Publisher Site | Google Scholar
J. N. Chen, T. M. Nolan, H. X. Ye et al., “Arabidopsis WRKY46, WRKY54, and WRKY70 transcription factors are involved in brassinosteroid-regulated plant growth and drought responses,” Plant Cell, vol. 29, no. 6, pp. 1425–1439, 2017.
View at: Publisher Site | Google Scholar
J. Li, S. Besseau, P. Törönen et al., “Defense-related transcription factorsWRKY70 andWRKY54 modulate osmotic stress tolerance by regulating stomatal aperture inArabidopsis,” New Phytologist, vol. 200, no. 2, pp. 457–472, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yan Zhong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

603

Downloads

843

Citations