Comparative transcriptome analysis reveals gene expression differences between two peach cultivars under saline-alkaline stress
Hereditas volume 157, Article number: 9 (2020)
Saline-alkaline stress is a major abiotic stress that is harmful to plant growth worldwide. Two peach cultivars (GF677 and Maotao) display distinct phenotypes under saline-alkaline stress. The molecular mechanism explaining the differences between the two cultivars is still unclear.
In the present study, we systematically analysed the changes in GF677 and Maotao leaves upon saline-alkaline stress by using cytological and biochemical technologies as well as comparative transcriptome analysis. Transmission electron microscopy (TEM) observations showed that the structure of granum was dispersive in Maotao chloroplasts. The biochemical analysis revealed that POD activity and the contents of chlorophyll a and chlorophyll b, as well as iron, were notably decreased in Maotao. Comparative transcriptome analysis detected 881 genes with differential expression (including 294 upregulated and 587 downregulated) under the criteria of |log2 Ratio| ≥ 1 and FDR ≤0.01. Gene ontology (GO) analysis showed that all differentially expressed genes (DEGs) were grouped into 30 groups. MapMan annotation of DEGs showed that photosynthesis, antioxidation, ion metabolism, and WRKY TF were activated in GF677, while cell wall degradation, secondary metabolism, starch degradation, MYB TF, and bHLH TF were activated in Maotao. Several iron and stress-related TFs (ppa024966m, ppa010295m, ppa0271826m, ppa002645m, ppa010846m, ppa009439m, ppa008846m, and ppa007708m) were further discussed from a functional perspective based on the phylogenetic tree integration of other species homologues.
According to the cytological and molecular differences between the two cultivars, we suggest that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes are crucial for saline-alkaline resistance in GF677. The results presented in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding to improve saline-alkaline tolerance in peach.
Saline-alkaline stress adversely affects plant growth and development. Based on the Food and Agriculture Organization of the United Nations (FAO) research, at least 242 million hectares of saline-alkaline soils are distributed within the Eurasian Region .
High salt concentrations and high pH levels of saline-alkaline soil generate reactive oxygen species (ROS), repress photosynthesis, energy production, and lipid metabolism, and damage the plant cell membrane and intracellular components . Upon saline-alkaline stress, a large cluster of gene expression is reprogrammed. Revealing the molecular mechanisms is helpful for saline-alkaline resistant plant breeding. In recent decades, several key genes such as OsLOL5, Gshdz4, and SsMT2, have been cloned and showed notable tolerance to saline-alkaline stress in transgenic plants [3,4,5]. Recently, with the rapid development of bioinformatics, transcriptome technology has become a powerful tool for elucidating the gene regulation of networks in many species [6, 7]. Using RNA-Seq technology, the effects of gene expression on saline-alkaline stress have been widely studied in alfalfa (Medicago sativa L.), flax (Linum usitatissinum L.), jujube (Ziziphus jujuba Mill.), black locust (Robinia pseudoacacia L.), and Chinese plum (Prunus salicina Lindl.), and thousands of differentially expressed genes have been detected [8,9,10,11].
Peach (Prunus persica L.) is an important deciduous fruit tree around the world. In China, the peach planting area is more than 700 thousand hectares, and the yield reaches more than one million tons per year . The North China Plain and Sichuan Basin, which partly suffer from saline-alkaline stress, are the main productive areas of peach planting [13, 14]. Saline-alkaline stress inhabits peach plant growth and affects fruit development, finally causing a reduction of in yield. GF677 (P. amygdalus x P. persica), bred by the French Institut National de la Recherche Agronomique (INRA) in the 1960s, has higher salt and alkali tolerance than Maotao (P. persica) . In 2017, Chen et al. identify more than one thousand genes differentially expressed between roots of GF677 and Maotao in response to saline-alkaline stresses . Considering leaf chlorosis usually acts as an early symptom of this stress , we carried out cytological observations and comparative transcriptome analysis of GF677 and Maotao leaves under salt-alkali stress in the present study, to elucidate the molecular mechanism governing the salt-alkali tolerance of GF677. These findings will provide a theoretical basis for genetic improvement and breeding of peach under salt-alkali stress.
The field experiment was performed in the Late White Peach Planting Base, Jianyang City, Sichuan Province (N30°30′50.97″, E104°26′35.21″). The pH of the test soil was 8.44 and the total salt content was 0.46% . Two-years-old GF677 and Maotao plantlets were planted in the field with an inter-row spacing of 1 m × 4 m. As a control, two varieties were planted in a greenhouse with normal soil (pH = 7.05). Nine GF677 and nine Maotao plantlets were selected for further experiments. The leaves from each plantlet were grouped into two parts: one for RNA isolation and another for cytological and biochemical analysis.
Each cultivar leaf sample was sliced into several sections, and soaked in 2% glutaraldehyde overnight at 4 °C. After rinsing with ddH2O, the samples were stained with 1% OsO4 for 2 h. The stained samples were further washed with ddH2O and dehydrated with acetone. The dehydrated samples were embedded in an epoxy resin and then sliced into ultrathin sections for ultrastructural observation as described by Jiang et al. .
Physiological parameters measurement
Chlorophyll a and b
Leaves (0.2 g) of GF677 and Maotao were homogenized with 80% acetone. After centrifugation at 12,000 rpm for 10 min, the supernatant was transferred separately into a clean microcentrifuge tube. The chlorophyll a and chlorophyll b contents were measured with a spectrophotometer (TU-1810, Beijing, China) at 663 nm and 645 nm, respectively .
Leaves (0.5 g) of GF677 and Maotao were washed with ddH2O. The total Fe concentrations were determined by inductively coupled plasma atomic emission spectroscopy (ICP-AES; Fisons ARL Accuris, Ecublens, Switzerland) .
Leaves (0.5 g) of GF677 and Maotao were homogenized with liquid nitrogen, and peroxidase (POD) activity was determined according to the protocol of the Peroxidase Assay Kit (A084–3-1, Nanjing Jiancheng Bioengineering Institute, China).
RNA extraction and transcriptome sequencing
For each cultivar, nine plantlets were randomly selected for RNA extraction. Total RNA of each plantlet was extracted according to the protocol of the RNAprep Pure Plant Plus Kit (Tiangen Biotech Co., Ltd., Beijing, China). RNA from three plantlets was equally pooled for cDNA library construction and transcriptome sequencing. For library preparation, mRNA was first extracted using oligo (dT) and then broken into 200 bp segments by adding fragmentation buffer. The first-stranded cDNA was synthesized by random hexamers, and the second-stranded cDNA was obtained by adding dNTPs and DNA polymerase I. The purified second strand cDNA was finally amplified using PCR. The libraries were sequenced on an Illumina HiSeq™ 4000 platform (Majorbio, Shanghai, China), and all raw reads were submitted to NCBI with the accession number GSE100180.
All analyses in the present research were carried out at least in triplicate. Based on the t-test, GraphPad Prism 5 was used for statistical analysis. Double stars showed a 0.01 significant difference.
Cytological and physiological features of two peach cultivars in tolerating saline-alkaline stress
There was no difference between GF677 and Maotao in normal soil (Fig. S1). However, the phenotypes of the two varieties differed upon saline alkaline field planting (Fig. 1a). In contrast to GF677, the leaves of Maotao were clear yellowish. The TEM analysis of leaves revealed a large structural change between the two cultivars. The chloroplast of GF677 had a compact granum structure, but the structure of the granum was dispersive in the Maotao chloroplast (Fig. 1b-c). Further chlorophyll determination showed that the contents of chlorophyll a or b in Maotao were only 22.3% or 15.9% in comparison with GF677. The POD enzyme activity and Fe content decreased by 50% in Maotao (Table 1).
Illumina sequence and DEG identification
A total of 29,006,496 and 31,874,666 clean reads were obtained in Maotao and GF677, respectively. The Q20 was more than 96.00% and the Q30 was more than 90.00%, suggesting the high quality of the clean reads. Of all clean reads, 83.68% ~ 90.43% were mapped to the reference genome (Table 2). To identify the gene response to saline-alkaline stress, the criteria of FDR ≤ 0.01 and |log2 ratio| ≥ 1 were used as thresholds. Among the DEGs, 294 and 587 were significantly up- and downregulated GF677 versus Maotao (Fig. 2a). In this study, the upregulation represented genes highly expressed in GF677 and downregulated showed genes highly expressed in Maotao. To validate the expression of DEGs, 11 genes involved in photosynthesis (ppa007547, ppa012123), carbohydrate metabolism (ppa007458), antioxidation (ppa027053, ppa011202), metal transporter (ppa003097), stress (ppa008441), polyamine biosynthesis (ppa007732), and transcription factor (ppa007708, ppa010846, ppa016095) were selected for the validation of DEGs, and the results showed that their expression pattern was highly related to the DEG data (Fig. 2b). To functionally characterize the DEGs, GO analysis was performed. A total of 30 GO terms were grouped, including several stress-related terms, such as “transporter activity”, “antioxidant activity”, “response to stimulus” (Fig. 2c). In the “molecular function” part, nine groups were obtained, and the “catalytic activity” was the most abundant group. In the “cell component” part, nine groups were obtained, and the “cell parts” was the most abundant group. In the “biological process” part, 12 groups were obtained, and the “metabolic process” was the most abundant group (Fig. 2c).
Overview of metabolism processes by MapMan
To provide an overview of the DEGs at the metabolic process level, the MapMan tool was used. It was clearly shown that photosynthesis was upregulated in GF677, while cell wall degradation, secondary metabolism, and starch degradation were activated in Maotao (Fig. 3; Table S3).
Category and verification of genes related to antioxidation, iron metabolism, transcription factor
Seven genes encoding glutathione S-transferase (GST) and six genes encoding POD were detected in DEGs (Fig. 4a). Of these genes, all seven GST genes and four POD genes were highly upregulated in GF677. ppa011202m encoding a GST and ppa027053m encoding a POD were further confirmed by qRT-PCR (Fig. 4b).
Eight iron-related genes were detected in the present research. Seven of them were upregulated, including two genes (ppa013552m, ppa019605m) encoding a metal ion transporter, one gene (ppa026143m) encoding calcium-regulated channel protein, one gene (ppa001921m) encoding an iron-sulfur cluster, one gene (ppa003816m) encoding an inorganic ion transport, one gene (ppa003097m) encoding an iron ion transmembrane transporter, and one gene (ppa009703m) encoding ferritin (Fig. 4a). Of these genes, ppa003097m was further confirmed by qRT-PCR (Fig. 4b).
Five genes encoding bHLH, three genes encoding MYB, and three genes encoding WRKY were detected in this study. Interestingly, the expression of bHLH and MYB downregulated, in contrast to the upregulation of WRKY (Fig. 4a). Furthermore, ppa016095m encoding a bHLH, ppa010846m encoding an MYB, and ppa007708m encoding a WRKY were confirmed by qRT-PCR (Fig. 4b).
We further constructed phylogenetic trees by integrating Arabidopsis homologues and other species functional verification genes (Fig. 5). The results showed that for each phylogenetic tree, the peach genes were scattered into different clades (Fig. 5a-c).
As a serious abiotic stress, saline-alkaline stress affects plant growth and reduces yield. In peach, comparative transcriptome analysis of roots of GF677 and Maotao suggests osmotic pressure increase and redox balance are crucial for GF677 response to saline-alkaline stresses . Leaf chlorosis usually acts as an early symptom of this stress [10, 11]. Our previous research showed that the photosynthetic rate was significantly higher in the leaves of GF677 than in Maotao . In the present research, TEM analysis revealed that the chloroplast granum structure was more intact and the content of chlorophyll a/b was higher, in GF677, suggesting that GF677 maintains higher photosynthesis efficiency under saline-alkaline stress (Fig. 1b-c, Table 1). Iron is an important ion that plays a key role in sustaining the activity of photosynthesis-related enzymes [22, 23]. In apple, Wang et al.  reported that Fe deficiency induces significant downregulation of genes involved in photosynthesis. In the present study, we show that iron content is higher in GF677, which is consistent with the upregulation of several iron and photosynthesis-related genes.
Saline-alkaline stress causes the accumulation of ROS, which severely damages such cell components as DNA, lipids, proteins, and sugars . Detoxification of ROS is the most important way in plants to weaken oxidative harmful stress . Of the ROS scavenging enzyme family, POD and GST are the key members [27, 28]. Previous studies have shown that these proteins usually play a synergetic role in scavenging ROS. In cotton, Li et al. identified that the POD and GST proteins are upregulated upon salt treatment, and both these upregulation effects are further verified by qRT-PCR . In the present research, notably, we found that POD activity was higher in GF677 (Table 1). This result is in accordance with the upregulation of most POD-related genes. Additionally, we also found that all seven GST genes were upregulated in GF677 (Fig. 4). Considering that the transcript level of a gene can largely reflect its translation state generally, we speculate that GST and POD have a synergistic role in scavenging ROS in GF677, which is consistent with the previous study .
TFs are key regulators of gene expression and have a variety of important functions in the plant response to abiotic stress [30, 31]. The identification of TFs involved in saline-alkaline stress is crucial to reveal the innate molecular mechanisms. In the present research, 11 genes related to bHLH, MYB, and WRKY TFs were detected (Fig. 4). As the most extensive TF class in eukaryotes, bHLH is not only universally involved in plant growth and metabolism but also plays an important role in plant response to stress [32, 33]. Fan et al.  revealed that Ib subgroup bHLH genes (AtbHLH38, AtbHLH39, AtbHLH100, and AtbHLH101) can bind the promoter of SKB1 and negatively regulate the accumulation of iron. In the present study, two-bHLH encoding genes (ppa024966m and ppa010295m), together with AtbHLH38, AtbHLH39, AtbHLH100, and AtbHLH101, were grouped into a subgroup by phylogenetic analysis (Fig. 5a). This result suggests ppa024966m and ppa010295m may have a similar role as Ib subgroup bHLHs, which is further confirmed by the higher content of Fe in GF677. Babitha et al.  showed that overexpression of AtbHLH17 enhances tolerance to NaCl, mannitol and oxidative stress in transgenic lines. Notably, ppa0271826m and ppa002645m show high homology to AtbHLH17, suggesting their putative function in the regulation of saline alkaline stress. MYB is one of the largest TF families in plants . The R2R3-MYB subfamily plays a key role in diverse biological processes, especially in response to various stresses [37, 38]. In soybean, overexpression of GmMYB3a negatively regulates saline-alkaline stress-related genes . In Arabidopsis, Cui et al.  showed that overexpression of AtMYB20 can enhance salt tolerance by negatively regulating type 2C serine/threonine protein phosphatases. In the present study, the phylogenetic tree shows that ppa010846m and ppa009439m are homologous to AtMYB20 and GmMYB3a, indicating a putative role in the regulation of saline-alkaline stress (Fig. 5b). WRKY is a class of TFs unique to plants that are mainly involved in development and stress responses [41, 42]. In Arabidopsis, overexpression of AtWRKY46, GmWRKY13 or VvWRKY11 can both positively regulate salt and drought stress tolerance [43,44,45]. In Nicotiana benthamiana, overexpression of GhWRKY41 confers transgenic plant salt and drought stress tolerance . In the present study, phylogenetic analysis shows that ppa008846m is highly similar to GmWRKY13 and VvWRKY11, while ppa007708m together with GhWRKY41 and AtWRKY46 are grouped into same subgroup (Fig. 5c), suggesting that it exhibits a similar function in the regulation of saline-alkaline stress.
In the present study, the biochemical, cytological and transcriptome differences between GF677 and Maotao were systematically analysed under saline-alkaline stress. The structure of chloroplast granum was intact in GF677, but dispersed in Maotao. Functional analysis of 881 DEGs showed that photosynthesis was activated, whereas cell wall degradation, secondary metabolism, and starch degradation were repressed in GF677. Based on the phylogenetic tree integration of other species homologues, several stress-related TFs were further functionally discussed. We speculate that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes, are crucial for saline-alkaline resistance in GF677. The results described in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding for improving saline-alkaline tolerance in peach.
Guan QJ, Ma HY, Wang ZJ, Wang ZY, Bu QY, Liu SK. A rice LSD1-like-type ZFP gene OsLOL5 enhances saline-alkaline tolerance in transgenic Arabidopsis thaliana, yeast and rice. BMC Genomics. 2016;17:142.
Cao L, Yu Y, DuanMu H, Chen C, Duan X, Zhu P, et al. A novel Glycine soja homeodomain-leucine zipper (HD-zip) I gene, Gshdz4, positively regulates bicarbonate tolerance and responds to osmotic stress in Arabidopsis. BMC Plant Biol. 2016;16:184.
Jin S, Xu C, Li G, Sun D, Li Y, Wang X, et al. Functional characterization of a type 2 metallothionein gene, SsMT2, from alkaline-tolerant Suaeda salsa. Sci Rep. 2017;7:17914.
Chen L, Qu H, Xia L, Liu Y, Jiang H, Sun Y, et al. Transcriptome profiling of the fertile parentand sterile hybrid in tea plant flower buds. Hereditas. 2019;156:12.
Jiang Y, Song HY, He JR, Wang Q, Liu J. Comparative transcriptome analysis provides global insight into gene expression differences between two orchid cultivars. PLoS One. 2018;13:e0200155.
An YM, Song LL, Liu YR, Shu YJ, Guo CH. De Novo transcriptional analysis of alfalfa in response to saline-alkaline stress. Front Plant Sci. 2016;7:931.
Guo M, Li S, Tian S, Wang B, Zhao X. Transcriptome analysis of genes involved in defense against alkaline stress in roots of wild jujube (Ziziphus acidojujuba). PLoS One. 2017;12:e0185732.
Liu J, Wang Y, Li Q. Analysis of differentially expressed genes and adaptive mechanisms of Prunus triloba Lindl. Under alkaline stress. Hereditas. 2017;154:10.
Meng F, Luo Q, Wang Q, Zhang X, Qi Z, Xu F, et al. Physiological and proteomic responses to salt stress in chloroplasts of diploid and tetraploid black locust (Robinia pseudoacacia L.). Sci Rep. 2016;6:23098.
Ji P. Research on international competitiveness and export influence factors of china’s peach industryMaster Thesis, North West Agriculture and Forestry University, College of Economics and Management (In English Abstract); 2011.
Yang Z, Wang BS. Present status of saline soil resources and countermeasures for improvement and utilization in China. Shandong Agric Sci. 2015;47:125–30 (In English Abstract).
Li Y, Cao K, Zhu G, Fang W, Chen C, Wang X, et al. Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. Genome Biol. 2019;20:36.
Gullo G, Motisi A, Zappia R, Dattola A, Diamanti J, Mezzetti B. Rootstock and fruit canopy position affect peach [Prunus persica (L.) Batsch] (cv. Rich may) plant productivity and fruit sensorial and nutritional quality. Food Chem. 2014;153:234–42.
Chen D, Song H, Sun S, Li J, Tu M, Jiang G. Scanning electron microscopy observation and comparative transcriptome analysis of roots of GF677 and Maotao in response to saline-alkali stresses. Mol Plant Breed. 2017;15:1580–6 (In English Abstract).
Song HY. Effects of chlorophyll synthesis and photosynthetic characteristics between GF677 and wide peach on alkaline soilMaster Thesis, Sichuan Agricultural University, College of Horticulture (In English Abstract); 2017.
Lei GJ, Zhu XF, Wang ZW, Dong F, Dong NY, Zheng SJ. Abscisic acid alleviates iron deficiency by promoting root iron reutilization and transport from root to shoot in Arabidopsis. Plant Cell Environ. 2014;37:852–63.
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39.
Tong Z, Gao Z, Wang F, Zhou J, Zhang Z. Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol Biol. 2009;10:71.
Tu M, Song H, Chen D, Li J, Sun S, Jiang G. Comparison of photosynthetic characteristics and chlorophyll fluorescence parameters between GF677 and wild peach on alkaline soil in Sichuan Hilly Basin. Mountain Res. 2018;1:153–62 (In English Abstract).
Varotto C, Maiwald D, Pesaresi P, Jahns P, Salamini F, Leister D. The metal ion transporter IRT1 is necessary for iron homeostasis and efficient photosynthesis in Arabidopsis thaliana. Plant J. 2002;31:589–99.
Pereira EG, Oliva MA, Rosado-Souza L, Mendes GC, Colares DS, Stopato CH, et al. Iron excess affects rice photosynthesis through stomatal and non-stomatal limitations. Plant Sci. 2013;201-202:81–92.
Wang YX, Hu Y, Zhu YF, Baloch AW, Jia XM, Guo AX. Transcriptional and physiological analyses of short-term iron deficiency response in apple seedlings provide insight into the regulation involved in photosynthesis. BMC Genomics. 2018;19:461.
Pandhair V, Sekhon BS. Reactive oxygen species and antioxidants in plants: an overview. J Plant Biochem Biotechnol. 2006;15:71–8.
Choudhury FK, Rivero RM, Blumwald E, Mittler R. Reactive oxygen species, abiotic stress and stress combination. Plant J. 2017;90:856–67.
Guan T, Shen J, Fa Y, Su Y, Wang X, Li H. Resistance-breaking population of Meloidogyne incognita utilizes plant peroxidase to scavenge reactive oxygen species, thereby promoting parasitism on tomato carrying Mi-1 gene. Biochem Biophys Res Commun. 2017;482:1–7.
Kumar S, Trivedi PK. Glutathione S-transferases: role in combating abiotic stresses including arsenic detoxification in plants. Front Plant Sci. 2018;9:751.
Li W, Zhao F, Fang W, Xie D, Hou J, Yang X, et al. Identification of early salt stress responsive proteins in seedling roots of upland cotton (Gossypium hirsutum L.) employing iTRAQ-based proteomic technique. Front Plant Sci. 2015;6:732.
Ohama N, Sato H, Shinozaki K, Yamaguchi-Shinozaki K. Transcriptional regulatory network of plant heat stress response. Trends Plant Sci. 2017;22:53–65.
Dong X, Yang Y, Zhang Z, Xiao Z, Bai X, Gao J, et al. Genome-wide identification of WRKY genes and their response to cold stress in Coffea canephora. Forests. 2019;10:335.
Ohashi-Ito K, Fukuda H. Functional mechanism of bHLH complexes during early vascular development. Curr Opin Plant Biol. 2016;33:42–7.
Sun X, Wang Y, Sui N. Transcriptional regulation of bHLH during plant response to stress. Biochem Biophys Res Commun. 2018;503:397–401.
Fan H, Zhang Z, Wang N, Cui Y, Sun H, Liu Y, et al. SKB1/PRMT5-mediated histone H4R3 dimethylation of Ib subgroup bHLH genes negatively regulates iron homeostasis in Arabidopsis thaliana. Plant J. 2014;77:209–21.
Babitha KC, Ramu SV, Pruthvi V, Mahesh P, Nataraja KN, Udayakumar M. Co-expression of AtbHLH17 and AtWRKY28 confers resistance to abiotic stress in Arabidopsis. Transgenic Res. 2013;22:327–41.
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15:573–81.
Roy S. Function of MYB domain transcription factors in abiotic stress and epigenetic control of stress response in plant genome. Plant Signal Behav. 2016;11:e1117723.
Hajiebrahimi A, Owji H, Hemmati S. Genome-wide identification, functional prediction, and evolutionary analysis of the R2R3-MYB superfamily in Brassica napus. Genome. 2017;60:797–814.
He Y, Yang X, Xu C, Guo D, Niu L, Wang Y, et al. Overexpression of a novel transcriptional repressor GmMYB3a negatively regulates salt-alkali tolerance and stress-related genes in soybean. Biochem Biophys Res Commun. 2018;498:586–91.
Cui MH, Yoo KS, Hyoung S, Nguyen HT, Kim YY, Kim HJ, et al. An Arabidopsis R2R3-MYB transcription factor, AtMYB20, negatively regulates type 2C serine/threonine protein phosphatases to enhance salt tolerance. FEBS Lett. 2013;587:1773–8.
Jiang J, Ma S, Ye N, Jiang M, Cao J, Zhang J. WRKY transcription factors in plant responses to stresses. J Integr Plant Biol. 2017;59:86–101.
Jing Z, Liu Z. Genome-wide identification of WRKY transcription factors in kiwifruit (Actinidia spp.) and analysis of WRKY expression in responses to biotic and abiotic stresses. Genes Genomics. 2018;40:429–46.
Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, et al. Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008;6:486–503.
Liu H, Yang W, Liu D, Han Y, Zhang A, Li S. Ectopic expression of a grapevine transcription factor VvWRKY11 contributes to osmotic stress tolerance in Arabidopsis. Mol Biol Rep. 2011;38:417–27.
Ding ZJ, Yan JY, Li CX, Li GX, Wu YR, Zheng SJ. Transcription factor WRKY46 modulates the development of Arabidopsis lateral roots in osmotic/salt stress conditions via regulation of ABA signaling and auxin homeostasis. Plant J. 2015;84:56–69.
Chu X, Wang C, Chen X, Lu W, Li H, Wang X, et al. The cotton WRKY gene GhWRKY41 positively regulates salt and drought stress tolerance in transgenic Nicotiana benthamiana. PLoS One. 2015;10:e0143022.