Next Article in Journal
Aldehyde Dehydrogenase 2 (ALDH2) in Rat Fatty Liver Cold Ischemia Injury
Next Article in Special Issue
Genome-Wide Bioinformatics Analysis of MAPK Gene Family in Kiwifruit (Actinidia Chinensis)
Previous Article in Journal
Full-Length Transcriptome Survey and Expression Analysis of Cassia obtusifolia to Discover Putative Genes Related to Aurantio-Obtusin Biosynthesis, Seed Formation and Development, and Stress Response
Previous Article in Special Issue
Assessment of Genetic Diversity, Population Structure, and Evolutionary Relationship of Uncharacterized Genes in a Novel Germplasm Collection of Diploid and Allotetraploid Gossypium Accessions Using EST and Genomic SSR Markers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 9
10.3390/ijms19092478
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis and Cloning of the Apple Stress-Associated Protein Gene Family Reveals MdSAP15, Which Confers Tolerance to Drought and Osmotic Stresses in Transgenic Arabidopsis

by
Qinglong Dong
†,
Dingyue Duan
†,
Shuang Zhao
,
Bingyao Xu
,
Jiawei Luo
,
Qian Wang
,
Dong Huang
,
Changhai Liu
,
Chao Li
,
Xiaoqing Gong
,
Ke Mao
* and
Fengwang Ma
*
State Key Laboratory of Crop Stress Biology for Arid Areas/Shaanxi Key Laboratory of Apple, College of Horticulture, Northwest A & F University, Yangling 712100, China
*
Authors to whom correspondence should be addressed.
†
These two authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(9), 2478; https://doi.org/10.3390/ijms19092478
Submission received: 8 July 2018 / Revised: 13 August 2018 / Accepted: 13 August 2018 / Published: 21 August 2018
(This article belongs to the Special Issue Plant Genomics)
Browse Figures
Versions Notes

Abstract

:
Stress-associated proteins (SAPs) are novel A20/AN1 zinc finger domain-containing proteins that are now favorable targets to improve abiotic stress tolerance in plants. However, the SAP gene family and their biological functions have not been identified in the important fruit crop apple (Malus × domestica Borkh.). We conducted a genome-wide analysis and cloning of this gene family in apple and determined that the overexpression of MdSAP15 enhances drought tolerance in Arabidopsis plants. We identified 30 SAP genes in the apple genome. Phylogenetic analysis revealed two major groups within that family. Results from sequence alignments and analyses of 3D structures, phylogenetics, genomics structure, and conserved domains indicated that apple SAPs are highly and structurally conserved. Comprehensive qRT-PCR analysis found various expression patterns for MdSAPs in different tissues and in response to a water deficit. A transgenic analysis showed that the overexpression of MdSAP15 in transgenic Arabidopsis plants markedly enhanced their tolerance to osmotic and drought stresses. Our results demonstrate that the SAP genes are highly conserved in plant species, and that MdSAP15 can be used as a target gene in genetic engineering approaches to improve drought tolerance.

1. Introduction

The growth, development, and survival of plants is constantly challenged by a variety of biotic and abiotic environmental factors. Plants utilize complex molecular mechanisms that regulate patterns of gene expression to protect themselves against these stresses [1,2]. Some key modulators of stress responses have been characterized and have emerged as appropriate targets to enhance abiotic stress tolerance in many plants. They include NAC domain-containing transcription factors, DRE/CRT-binding transcription factors (DREBs/CBFs), mitogen-activated protein kinases (MAPKs), stress-associated proteins (SAPs), and heat shock factor/proteins (HSPs/HSF) [1,3,4,5,6]. Among these, the SAPs are a newly identified class of zinc finger proteins (ZFPs) that play crucial roles in various abiotic stress responses by numerous plants [1,2,7].
The SAP gene family members have two special ZF domains: the highly conserved A20 domain, which was first isolated in human umbilical vein endothelial cells with the characterization of a tumor necrosis factor (TNF)-α-inducible protein; and/or the AN1 domain, which is also highly conserved and first identified from Xenopus laevis animal hemisphere 1 (AN1) maternal RNA with the delineation of the ubiquitin-like protein [8,9]. The SAP proteins expressed in Arabidopsis thaliana (hereafter Arabidopsis), rice (Oryza sativa), tomato (Solanum lycopersicum), and cotton (Gossypium hirsutum) have been classified into five groups (I through V) based on results from their phylogenetic analyses [10,11]. One significant feature of plant SAPs is the very frequent occurrence of intronless genes [2]. For example, 11 rice SAP genes, 15 from desert poplar (Populus euphratica), and 30 from cotton lack introns and show a remarkably higher percentage of intronless genes [2,6,11].
The roles of SAP genes are increasingly being reported in plants. Transcriptional levels are induced by multiple stresses and provide a positive reinforcement of tolerance to abiotic stress. Rice A20/AN1 protein (OSISAP1/OsSAP1), the first identified plant SAP gene, is induced after different types of stress treatments are applied [1]. The overexpression of OSISAP1 confers tolerance to dehydration, cold, and salt in transgenic seedlings of tobacco (Nicotiana tabacum) [1]. Furthermore, ZFP177 (OsSAP9) and AtSAP5 are induced by numerous challenges and have significant roles in improving abiotic stress tolerance [12,13]. Similar results have been described for SAPs from maize (Zea mays) [14], medicago (Medicago truncatula) [15], banana (Musa sp.) [16], the halophyte grass Aeluropus littoralis [17], and poplar (Populus alba × P. glandulosa) [18]. These genes also function in biotic stress responses. For example, Tyagi et al. [19] have analyzed the expression patterns of rice SAP gene family members in response to pathogen elicitors, and discovered that OsSAP1, OsSAP8, and OsSAP11 are up-regulated. Transgenic tobacco overexpressing OsiSAP1 shows significantly enhanced basal resistance against infection by the bacterial pathogen Pseudomonas syringae pv. Tabaci [19].
The SAP genes also help regulate signal transduction and phytohormone synthesis. In rice, the overexpression of OsDOG (OsiSAP11) [20] and OsZFP185 (OsiSAP4) [21] results in dwarf phenotypes, a decrease in gibberellic acid (GA) contents, and deficient cell elongation. Furthermore, OsZFP185 negatively regulates the expression of several genes related to abscisic acid (ABA) biosynthesis, and interferes with ABA-mediated tolerance to salt, drought, and cold [21]. Various SAPs can function as E3 ubiquitin ligases, redox sensors, and/or regulators of gene expression under stress [7,13,22,23]. Other novel biological functions for SAP genes will continue to be reported.
Based on their highly conserved A20/AN1 domains, members of the SAP gene family have been identified and characterized in Arabidopsis [24], rice [24], maize [14], tomato [10], cotton [11], desert poplar [6], and medicago [25]. Although extensive genomic analyses have provided considerable details about this family in several species, members in apple (Malus × domestica Borkh.) have not been as thoroughly investigated. Nevertheless, recent completion of the draft genome sequence for apple has enabled genome-wide analyses of its SAP genes [26,27,28]. Here, we identified SAP members in apple and examined their A20/AN1 domain, protein and gene structures, conserved domains, phylogenetic relationships, chromosomal locations, cis-acting elements, and expression patterns for MdSAPs cloned in response to water deficits. We also overexpressed MdSAP10 in Arabidopsis and investigated its function. Our results will serve as a basis for exploring the molecular roles of SAPs. By facilitating further studies into their functions in abiotic stress responses, we can continue our efforts to introduce improved apple cultivars.

2. Results

2.1. Identification and Annotation of Apple SAP Genes

To identify the genes in the apple genome that encode SAP proteins, we conducted a BlastP of the apple genome database and identified 32 putative family members. We then used the Pfam and NCBI Conserved Domain Database (NCBI-CDD) databases for verification, searching for the A20/AN1 domain in the amino acid sequences encoded by all 32 genes. From this, we confirmed the identity of 30 typical apple SAP genes in the original dataset (Table 1).
We next cloned all of the full-length apple SAP genes based on predicted nucleotide sequences in the apple genome and in the NCBI nucleotide database. As shown in Table 1, this revealed that the full-length cDNAs of MdSAP7, -8, -10, -12, -14, -15, -16, -19, -21, -23, -25, -28, and -29 had been isolated and confirmed by RT-PCR (Supplementary Sequence A1). Their corresponding 5′- and 3′-UTRs were then amplified.

2.2. Structures and Conserved Domains of Apple SAP Genes

To gain insights into the structural diversity of SAP genes in apple, we analyzed the phylogenetic tree, exon–intron organization, and conserved domains in the coding sequences. The MdSAP proteins were classified as groups I and II based on their phylogenetic relationships (Figure 1A). Gene structure analysis indicated that MdSAP1 through MdSAP4, MdSAP6 through MdSAP21, MdSAP23, MdSAP24, MdSAP26, MdSAP29, and MdSAP30 contained no introns, whereas MdSAP5, MdSAP22, MdSAP25, MdSAP27, and MdSAP28 had one each (Figure 1B). Conserved domain analysis revealed that all of the MdSAP proteins included A20 and/or AN1 domain(s). MdSAP1 through MdSAP4, MdSAP7 through MdSAP20, MdSAP23, MdSAP24, MdSAP26, MdSAP29, and MdSAP30 contained an A20 domain and an AN1 domain; MdSAP5, MdSAP6, and MdSAP22 had single AN1 domains; and MdSAP21, MdSAP25, MdSAP27, and MdSAP28 each had two AN1 domains. In addition, MdSAP25 contained a C2H2 domain at the C terminal (Figure 1C and Figure 2).

2.3. Multiple Sequence Alignments and Three-Dimensional Structure of Apple A20/AN1 Domains

Multiple alignments demonstrated that the A20/AN1 domains are conserved among the MdSAP proteins (Figure 2). We then produced sequence logos that further showed that these domains were highly conserved at each residue position (Figure 3A,C; Supplementary Sequences A2 and A3). Afterward, the SWISS-MODEL web server was used for modeling and analysis of homology among protein structures. For this, we built the A20 domain and AN1 domain homology models and evaluated them using the homologous templates 2KZY.pdb and 1WFP, respectively (Figure 3B,D). The 3D models indicated that, respectively, the A20 domain and the AN1 domain in the MdSAP7 structure most closely matched the A20 domain of ubiquitin receptor ZNF216 and the zf-AN1 domain of the Arabidopsis F5O11.17 protein (PDB ID: 2KZY.1.A, 38% sequence identity for residues 18–50; 1WFP.1.A, 47% sequence identity for residues 96–148).

2.4. Phylogenetic Analysis of SAP Proteins

To examine the evolutionary relationships among plant SAP proteins, we used MEGA 6 and constructed unrooted phylogenetic trees from full-length protein sequences encoded by 453 SAP genes in 32 species (Supplementary File B1). Two major groups were revealed: I, containing an A20 domain and an AN1 domain; and II, containing two AN1 domains. Members in Group I were further classified into four subgroups (Ia-Id), while Group II members were assigned to two subgroups: IIa, containing two AN1 domains and one or two C2H2 domain(s); and IIb, containing only two AN1 domains. Among these 30 MdSAP proteins, 25 (MdSAP1-MdSAP20, MdSAP23, MdSAP24, MdSAP26, MdSAP29, and MdSAP30) could be unambiguously classified as Group I, while four (MdSAP21, MdSAP25, MdSAP27, and MdSAP28) were assigned to Group II based on their relationship with the other SAP proteins (Figure 4 and Supplementary File B2). Further analysis revealed that the 25 Group-I apple proteins belonged to subgroups Ia (seven genes), Ib (five), Ic (five), and Id (eight). Subgroup IIa contained MdSAP25, while IIb contained MdSAP21, MdSAP27, and MdSAP28.

2.5. Genome Distribution of Apple SAP Genes

We determined the genomic locations of these apple SAP genes based on their mapping coordinates. In all, 28 of the 30 MdSAP genes were assigned to chromosomes 1–4, 6–9, 11, 12, 14, and 17 (Table 1; Figure 5). However, we could not conclusively map two genes (MdSAP29 and MdSAP30) to any chromosome. The genes were unevenly distributed among the 12 chromosomes, with Chromosome 2 containing the most (eight genes), followed by Chromosome 7 (five genes), and one each for chromosomes 1, 3, 6, 12, and 14.
The apple gene family appears to have expanded during the process of genome evolution [26]. To uncover the mechanism underlying this expansion, we investigated gene duplication events, including tandem and segmental duplications, and found that many MdSAP genes (19/30, or 63.33%) were present in two or more copies (Figure 5). In all, 17 had undergone tandem duplication, while two were subjected to segment duplication. Those segment duplications produced many homologs of SAP genes on different chromosomes, while tandem duplications produced SAP gene clusters or hotspots (blue and red font in Figure 5). A relatively recent genome-wide duplication is that in the Pyreae tribe, which was thought to result in the transition of nine ancestral chromosomes to 17 chromosomes [26]. We noted here that multiple gene pairs were each linked to at least six potential chromosomal segmental duplications (Figure 5, pairs of bars in grey areas), e.g., large sections of chromosomes 9 and 17, 3 and 11, and 7 and 2.

2.6. Promoter Sequence Analysis of Apple SAP Genes

To investigate putative cis-acting elements in their promoter regions, we isolated approximately 1500-bp genomic sequences upstream of the start codon from our MdSAPs. Along with some cis elements involved in light-responsiveness (Table 2 and Supplementary File B3), we found that many were responsive to various stresses and correlative hormones. In total, 10 types of cis elements were discovered in the 13 promoters. They were associated with responses to hypoxia, heat, chilling, drought, pathogens, wounding, or hormones such as salicylic acid, methyl jasmonate, ABA, or ethylene. Therefore, we concluded that these cis elements play important roles in plant stress responses.
Sequences and functions for ABRE (ABA response element), ARE (anaerobic response element), CGTCA (MeJA-responsiveness), ERE (ethylene-responsive element), HSE (heat shock response element), LTR (low-temperature response element), MBS (MYB binding site involved in drought response), TCA (salicylic acid response element), TC-rich repeat (defense and stress responsiveness), and W-box (elicitation; wounding and pathogen responsiveness; binding site of WRKY type transcription factors) were obtained from the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). Digits represent the number of regulatory elements on plus/minus strand. Blank space indicates no corresponding cis-acting element in either strand of the promoter.

2.7. Expression Profiles of MdSAP Genes

Knowing the patterns of expression in various tissue types can help us understand gene functions. Firstly, we collected different tissues including young roots, stems, fully expanded leaves, flowers, and mature fruit (70 mm, red peel, 150 days after bloom), from apple plants that were five years old after bud grafting. The scion was Malus domestica “Golden Delicious”, and the rootstock was M. hupehensis. Then, we isolated the full-length cDNA, 5′-UTR, and 3′-UTR sequences for 13 MdSAP genes and used specific primers for our qRT-PCR assays. These MdSAPs were constitutively expressed in the five tissues examined here, albeit at different levels of transcription (Figure 6A). For example, MdSAP7, -8, -10, -15, -19, and -29 were most highly expressed in the leaves; MdSAP12, -16, -21, and -28 were most highly expressed in the fruits; and MdSAP14, -23, and -25 were most highly expressed in the roots. To induce a water deficit, irrigation was withheld up to 8 days, while the designated control plants continued to receive normally scheduled irrigation (“Golden Delicious” scions and M. hupehensis rootstocks). In response to drought stress, the expression of MdSAP15, -25, and -28 was significantly induced from that detected in the non-stressed control plants (Figure 6B), transcripts of MdSAP7 and -21 mRNAs were significantly reduced, and expression of the other MdSAP genes remained constant.

2.8. MdSAP15 Overexpression Enhances Osmotic Stress Tolerance by Arabidopsis Seedlings

Since the transcription of MdSAP15 mRNA was significantly accumulated under drought stress, we chose this gene for investigating its biological functions in Arabidopsis. After kanamycin-resistance screening and PCR detection using Arabidopsis genomic DNAs as templates, more than five transformants were identified and confirmed, with elevated levels of MdSAP15 transcripts (Figure 7A). From these, we selected three transgenic lines (L1, 4, and 5) with high MdSAP15 expression to evaluate its potential functioning in response to osmotic and drought stresses. For the osmotic stress assay, five-day-old seedlings grown on Murashige and Skoog (MS) agar plates were vertically plated on an MS agar medium supplemented with 0, 200, or 300 mM of mannitol. While the roots of “Col” and MdSAP15 OE seedlings displayed similar growth characteristics on MS agar medium plates (Figure 7B), their growth was affected when treated with different concentrations of mannitol. For example, the primary roots and fresh weights of the transgenics were longer and heavier than those of the wild-type (WT) “Col” upon exposure to 200 or 300 mM of mannitol (Figure 7C,D).

2.9. MdSAP15-Overexpressing Arabidopsis Seedlings Have Improved Physiological Traits Associated with Osmotic Stress Tolerance When Compared with “Col” Wild Type

For further investigation of this MdSAP15-mediated enhancement of tolerance to osmotic stress, we measured relative electrolyte leakage (REL) and concentrations of chlorophyll, proline, and malondialdehyde (MDA), all of which are important markers of such tolerance. Under normal growing conditions, we did not observe any obvious differences in REL. However, in response to mannitol exposure, REL values were significantly higher in “Col” seedlings than in the transgenics (Figure 7E). Levels of chlorophyll and proline were similar between “Col” and OE plants under control conditions, but were higher in the transgenics after they were treated with 200 mM of mannitol (Figure 7F,G). Furthermore, MDA concentrations did not differ among genotypes under normal conditions, but were lower in the overexpressed (OE) lines than in the WT following mannitol treatment (Figure 7H). These results suggested that the overexpression of MdSAP15 in Arabidopsis seedlings leads to enhanced osmotic stress tolerance.

2.10. MdSAP15 Overexpression Enhances Drought Tolerance in Transgenic Arabidopsis Plants

To evaluate drought tolerance, we simultaneously germinated seeds of the “Col” and MdSAP15-overexpressing lines and grew the seedlings on MS plates for one week before transplanting them into soil for another three weeks of culture. Drought conditions were then imposed by withholding water. After 20 days of stress, all of the plants exhibited symptoms related to severe water loss, although those symptoms were milder for the OE transgenics than for the WT. When rehydration began (2 days of refreshment), most of the “Col” plants did not recover (Figure 8A), whereas the survival rate was significantly higher for the three transgenic lines (Figure 8B). These results again indicated that the overexpression of MdSAP15 in Arabidopsis plants enhances their degree of drought tolerance.

3. Discussion

The proteins encoded by SAP genes comprise large families and are broadly distributed in higher plants [2]. Apple is an economically important woody plant and the most widely cultivated fruit crop in the world. Sequencing of its genome has provided a good platform for genome-wide analyses of all putative gene families in apple, including the DREB [29], MYB [30], MADS-box [31], and WRKY [32,33] families. However, genome-wide information about apple SAP genes has remained unknown, while members of that family have been identified in other plant species [6,10,11,14,24,25]. Moreover, the content of SAP genes varies substantially among species. For example, Brassica rapa, Glycine max, Solanum tuberosum, Salix purpurea, Populus trichocarpa, and cotton each have a relatively large number of SAP members, i.e., 28, 26, 19, 19, 19, and 19, respectively; while Chlamydomonas reinhardtii, Lotus japonicus, Carica papaya, and Amborella trichopoda have relatively few, i.e., 2, 6, 7, and 7, respectively (Table 3). Here, we determined that the apple genome contains 30 SAP genes, making this family much larger than in any other species.
Segmental, tandem, and whole-genome duplications are critical for both the diversification of gene functions and the rearrangement and expansion of genomes [31,32,33,34]. Whole-genome duplication events have occurred in apple [26], and tandem, segmental, and whole-genome duplications have caused some apple gene families to expand, including the MYB [30], MADS-box [31], and WRKY [32] families. We learned here that two MdSAP genes have undergone segmental duplication, and 17 have undergone tandem duplication. In addition, multiple gene pairs have each been linked to six potential chromosomal segmental duplications (Figure 5). Similar results have been reported for the Medicago SAP gene family. Our findings suggest that transposition events and the whole-genome and chromosomal segmental duplications have led to the expansion of the apple SAP gene family, and might partially explain why more SAP genes are present in apple than in any other species.
In several distinct species, the zinc finger types of some family members have either disappeared or increased in number. For example, the A20-A20-AN1 zinc finger occurs only in rice and Eucalyptus grandis; the A20 type is found in apple, rice, Amborella trichopoda, B. rapa, and grape (Vitis vinifera); the AN1-AN1-C2H2 zinc-finger exists in Arabidopsis thaliana, A. lyrata, B. rapa, Capsella rubella, Thellungiella parvula, and desert poplar (Table 3). We might speculate that the loss or the increase in zinc finger types of SAP genes in these genomes means that they are critical for the complicated enzymatic activity that is present in those species. In this study, we determined that apple SAPs are highly and structurally conserved based on analyses of gene structure, conserved domains, sequence alignments, 3D structures, and phylogenetics (Figure 1, Figure 2, Figure 3 and Figure 4). Similar results have been reported for Arabidopsis, rice, maize, tomato, cotton, desert poplar, and medicago [6,10,11,14,24,25].
Although members of the SAP gene family in Arabidopsis, rice, tomato, and cotton are arranged into five groups [10,11,24], the SAP proteins expressed in Arabidopsis, desert poplar, Populus trichocarpa, Salix purpurea, and S. suchowensis are classified into two major groups: I (Ia–If) and II (IIa and IIb) [6]. Due to the small number of plant species that have been examined, we cannot yet accurately analyze the evolutionary relationships within that gene family. In this study, we were able to divide those members into two major groups (Ia–Id/IIa and IIb) by comparing the apple genome with SAP proteins from 31 other species (Figure 4). We believe that this result from our evolution analysis is convincing. Usually, exon–intron structural diversity can provide important evidence for phylogenetic relationships and play a valuable role in the evolution of gene families [31]. An intronless structure is typical of SAP genes in various species, and is a key characteristic of that family. However, one exception is the grape genome, for which only two VvSAP members lack introns, while 10 members each contain one intron. Furthermore, we discovered here that the intronless gene structure of SAP genes is the dominant arrangement among Group I members, whereas most of the Group II members contain at least one intron (Table 4 and Supplementary File B2). Thus, the prevalence of an intronless gene structure reflects the ancient origin of SAP genes, and links well with their rapid accumulation of transcripts due to reduced post-transcriptional processing [2,35].
Plant SAPs are quickly induced by multiple abiotic stresses [1,12,13,14,36,37,38,39]. They include rice OsiSAP1/OsSAP1, which responds to drought, salt, cold, submergence, mechanical wounding, and ABA [1]; and ZFP177 (OsSAP9), which is also from rice, and shows enhanced expression in response to cold, heat, and PEG6000 [12]. The expression of OsiSAP8 in tobacco and rice is enhanced by salt, cold, heat, desiccation, wounding, submergence, heavy metals, and ABA [36]. Similarly, SAP genes in Aeluropus littoralis, banana, Arabidopsis, and maize respond to salt, cold, drought, and osmotic stresses in a tissue and stress-specific manner [13,14,16,37,38,39]. In this study, we comprehensively analyzed the expression patterns of 13 cloned SAP genes under drought stress. Whereas the expression of MdSAP15, -25, and -28 was significantly induced, transcripts levels for MdSAP7 and -21 mRNAs were significantly reduced (Figure 6B). Our results suggest that these genes have important roles in the response to water deficits.
The constitutive expression of SAP genes confers tolerance to multiple challenges. The overexpression of rice OsiSAP1/OsSAP1, OsiSAP8, and AlSAP in tobacco and rice increases their tolerance to numerous abiotic stresses [1,17,36,37,40]. Similar findings have been described for the overexpression of AtSAP5 in cotton and Arabidopsis [13,41]. The overexpression of AtSAP13 and MusaSAP1 in Arabidopsis and banana leads to greater drought and salt tolerances [16,39]. The downregulation of PagSAP1 improves salt tolerance in poplar and alters the regulation of genes involved in maintaining cellular ionic homeostasis [18]. We noted here that MdSAP15 overexpression conferred increased tolerance to osmotic stress by increasing the root lengths and fresh weights of transgenic Arabidopsis seedlings when compared with the WT. This overexpression also influenced a range of parameters associated with abiotic stress responses, including REL and the concentrations of chlorophyll, proline, and MDA, all of which are often used to evaluate the degree of plant tolerance under stress conditions [42,43,44]. Measured values for all of them were favorably affected in our transgenic lines (Figure 7). Finally, our experiments with induced water deficits demonstrated that the transgenic Arabidopsis plants showed milder stress symptoms when compared with the WT, and they also had higher survival rates during the period of rehydration and recovery from drought treatment (Figure 8). Taken together, these results indicate that the overexpression of MdSAP15 in Arabidopsis plants leads to enhanced drought tolerance. This work provides a basis for exploring the molecular roles of SAPs and facilitates further investigations into the functions of these genes in abiotic stress responses. Our data also lay a solid foundation for future efforts to introduce improved apple cultivars.

4. Materials and Methods

4.1. Identification of Apple SAP Genes

We downloaded the database of the Arabidopsis SAP family from the TAIR website (http://www.arabidopsis.org/) [24]. As query sequences for BlastP (http://www.rosaceae.org/tools/ncbi_blast) against predicted apple proteins, we used 14 Arabidopsis SAP proteins and the consensus protein sequences of the A20/AN1 domain hidden Markov model (HMM) profile (A20-like zinc finger, PF01754; AN1 zinc finger, PF01428) from the Pfam database (http://pfam.xfam.org/family/PF01754; http://pfam.xfam.org/family/PF01428). We then searched all of those SAP sequences against the apple genome database (https://www.rosaceae.org/gb/gbrowse/malus_x_domestica/) with HMMER v3.0 and BlastP [31,32]. Confirming the reliability of those protein sequences ensured that the A20 and/or AN1 domains were present in each candidate MdSAP protein. For this, we used the Pfam database (http://pfam.sanger.ac.uk/search) and NCBI Conserved Domain Database (NCBI-CDD; http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) [34].

4.2. Sequence Alignments and Phylogenetic Analysis

We performed multiple sequence alignments of 30 MdSAP protein sequences, using DNAMAN 6.0.3.99 with its default parameters [45]. A phylogenetic tree for the SAP gene family was constructed with MEGA 6.0 software (www.megasoftware.net) and the neighbor-joining (NJ) method, together with the full-length amino acid sequences of 453 SAPs from 32 plant species, including apple (Supplementary Sequence A2). Related sequences were downloaded from the resource Plaza 3.0 (http://bioinformatics.psb.ugent.be/plaza/). We used the following parameters in the NJ method: bootstrap (1000 replicates), complete deletion, and amino: p-distance [34].

4.3. Sequence Logos and Structure Model Analysis

Sequence logos for the A20 domain in 26 MdSAP genes and the AN1 domain in 28 MdSAP genes were generated by the application WebLogo (http://weblogo.threeplusone.com) (Supplementary Sequences A3 and A4). We used the web server SWISS-MODEL (http://swissmodel.expasy.org/) for modeling and predicting the homology of protein structures for those two domains [46]. The proposed 3D structure was modeled on the original NMR structure in PDB ID: 2KZY and 1WFP, and RasTop 2.2 software (http://www.geneinfinity.org/rastop/) was used to present that model [47].

4.4. Analyses of Intron–Exon Structure, Genome Distribution, and Gene Duplications

Genomic sequences (apple v1.0), gene distributions on chromosomes, and genome locations of SAPs in apple and 31 other species were downloaded from Plaza 3.0 (Supplementary Sequence A1) [48]. The structural features of these MdSAP genes, including exons/introns, numbers, and locations, were obtained and presented by using the gene structure display server (GSDS) web-based bioinformatics tool (http://gsds.cbi.pku.edu.cn/) [49]. The chromosomal positions of all of the MdSAP genes were located via MapInspect (www.plantbreeding.wur.nl/UK/software_mapinspect.html) [50]. Segmental and tandem-duplication events were investigated according to the method of Tian et al. [31].

4.5. Prediction of Cis-Acting Elements in Promoters

To examine the putative cis-acting elements in the promoters of apple SAP genes, we isolated sequences that were 1500 bp upstream of the translational start codon, using the contig sequences of that genome and PCR amplification. Details for the promoters used here are listed in Table 5 and Supplementary File B3. Possible cis-acting elements in those promoters were then predicted according to the Plant CARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [51].

4.6. Plant Materials, Growth Conditions, and Stress Treatments

Young roots, stems, and fully expanded leaves, as well as flowers and mature fruit (70 mm, red peel, 150 days after bloom), were collected from apple plants that were five years old after bud grafting. The scion was Malus domestica “Golden Delicious”, and the rootstock was M. hupehensis. Samples used for examining the effects of water deficits were harvested from plants three months after bud grafting was performed with “Golden Delicious” scions and M. hupehensis rootstocks. These grafted plants were grown in pots (height, 320 mm; diameter, 300 mm) in a greenhouse and treatments began when the plants were approximately 500-mm tall. To induce a water deficit, irrigation was withheld from certain plants for up to 8 days while the designated control plants continued to receive normally scheduled irrigation [52]. Our sampling schedule involved harvesting mature leaves at the middle nodes on days 0, 4, and 8 of the deficit period. All of the tissues were frozen immediately in liquid N2 and stored at −80 °C.
Seedlings of Arabidopsis thaliana L. (Heyn), cv. Columbia (“Col”), were used for genetic transformations and assays of osmotic and drought tolerance. They were cultured in a growth chamber under a 16-h photoperiod at 23 °C. For the drought tolerance assay, water was withheld from four-week-old plants for 20 days before they were rewatered. Survival rates were scored 2 days after rewatering began. Well-watered plants were used as the negative control. For the osmotic stress assay, five-day-old seedlings grown on MS agar plates were vertically plated on an MS agar medium supplemented with 0, 200, or 300 mM of mannitol. Their root lengths, fresh weights, relative electrolyte leakage (REL), and concentrations of chlorophyll, malondialdehyde (MDA), and proline were measured 11 days after that transfer. All of the experiments were repeated three times.

4.7. Cloning of MdSAPs and qRT-PCR Analysis

We extracted total RNA from previously frozen apple tissues according to the CTAB method and from Arabidopsis leaves, using Trizol reagent (Thermo Fisher Scientific-CN; Shanghai, China; https://www.thermofisher.com/cn/zh/home.html) [53]. Two micrograms of total RNA were collected for synthesizing first-strand cDNA. For cloning MdSAP, complete open reading frames were obtained via RT-PCR from fully expanded leaves of “Golden Delicious” apple, using specific primers listed in Table 5. The 5′- and 3′-untranslated regions (UTRs) were obtained with a Rapid Amplification for cDNA Ends kit (TaKaRa, Dalian, China). For the qRT-PCR assays, reverse transcription was performed with 1 μg of total RNA from each sample, followed by PCR-amplification of 1 μL of the product. We conducted the qRT-PCR assays in 20-µL reaction mixtures that contained 10 µL of SYBR® Premix Ex Taq™ (TaKaRa; Beijing, China; http://www.takarabiomed.com.cn), and used an iQ5 instrument (Bio-Rad, Hercules, CA, USA) as described before [34]. Thermal cycling included an initial 3 min at 95 °C; then 40 cycles of 10 s at 95 °C, 30 s at 58 °C, and 15 s at 72 °C; followed by 3 min at 72 °C and then 81 cycles of 7 s each, increasing by an increment of 0.5 °C from 55 °C to 95 °C. Three biological replicates were tested in each assay, and ∆Ct values were calculated by using MdMDH as our endogenous control [54]. Relative quantification was calculated according to the 2−∆∆Ct method [55], and dissociation curve analysis was performed for determining the specificity of the amplifications.

4.8. Vector Construction and Plant Transformation

To construct the MdSAP15 overexpression (OE) vectors, we performed RT-PCR to isolate the full-length cDNA of MdSAP15 from fully expanded leaves of the “Golden Delicious” apple. The cDNA was cloned into pRI 101-AN plant transformation vectors that were driven by the cauliflower mosaic virus (CaMV) 35S promoter. For Arabidopsis transformation, the recombinant plasmid described above was introduced into the “Col-0” ecotype via the Agrobacterium tumefaciens GV3101-mediated floral dip method. Seeds of the transgenic plants were individually harvested and screened with kanamycin monosulfate. Homozygous transgenic lines were used for further investigations.

4.9. Measurements of Physiological Indices

Relative electrolytic leakage was examined as described by Tan et al. [56]. Chlorophyll concentrations were determined using the protocol of Liang et al. [57], while MDA levels were obtained as described by Wei et al. [42]. Proline concentrations were calculated according to the method described by Dong et al. [43].

4.10. Statistical Analysis

All of the data were analyzed with IBM SPSS Statistics v.20 software (https://www.ibm.com/support/knowledgecenter/SSLVMB_20.0.0/com.ibm.spss.statistics_20.kc.doc/pv_welcome.html). One-way ANOVA and Duncan’s tests were used to compare the results. Differences between treatments were considered statistically significant at p < 0.05.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/9/2478/s1. Supplementary Sequence A1: Genomic information for cloned MdSAPs; Supplementary Sequence A2: Full-length amino acid sequences of 453 SAPs from 32 plant species; Supplementary Sequence A3: Sequence logo for the A20 domain in 26 MdSAP genes; Supplementary Sequence A4: Sequence logo for the AN1 domain in 28 MdSAP genes; Supplementary File B1: The IDs of SAP gene family members in various species; Supplementary File B2: The intron numbers and groups of SAP gene family members in various species; Supplementary File B3: Prediction of cis-acting elements in MdSAP promoters.

Author Contributions

F.M., K.M. and Q.D. collected the public dataset, performed bioinformatics analysis, and drafted the manuscript. D.D., S.Z., B.X., J.L., Q.W. and D.H. contributed to the bioinformatics analysis and preparation of all figures and tables. Q.D., D.D., C.L. (Changhai Liu), C.L. (Chao Li) and X.G. conducted the experiments. F.M. and Q.D. conceived this study and reviewed the manuscript. All authors have read and approved the final manuscript.

Funding

This work was supported by the State Key Program of the National Natural Science Foundation of China (31330068), the earmarked fund for the China Agriculture Research System (CARS-27), the National Natural Science Foundation of China (31401852 and 31701894), the China Postdoctoral Science Foundation (2017M620474), and the Fundamental Research Funds for the Central Universities (2452017064 and 2452016186).

Acknowledgments

The authors are grateful to Priscilla Licht for help in revising our English composition.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mukhopadhyay, A.; Vij, S.; Tyagi, A.K. Overexpression of a zinc-finger protein gene from rice confers tolerance to cold, dehydration, and salt stress in transgenic tobacco. Proc. Natl. Acad. Sci. USA 2004, 101, 6309–6314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Giri, J.; Dansana, P.K.; Kothari, K.S.; Sharma, G.; Vij, S.; Tyaqi, A.K. SAPs as novel regulators of abiotic stress response in plants. BioEssays 2013, 35, 639–648. [Google Scholar] [CrossRef] [PubMed]
  3. Mizoi, J.; Shinozaki, K.; Yamaguchi-Shinozaki, K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 2012, 1819, 86–96. [Google Scholar] [CrossRef] [PubMed]
  4. Nakashima, K.; Takasaki, H.; Mizoi, J.; Shinozaki, K.; Shinozaki, K.; Yamaguchi-Shinozaki, K. NAC transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 2012, 1819, 97–103. [Google Scholar] [CrossRef] [PubMed]
  5. Scharf, K.D.; Berberich, T.; Ebersberger, I.; Nover, L. The plant heat stress transcription factor (Hsf) family: Structure, function and evolution. Biochim. Biophys. Acta 2012, 1819, 104–109. [Google Scholar] [CrossRef] [PubMed]
  6. Jia, H.; Li, J.; Zhang, J.; Ren, Y.; Hu, J.; Lu, M. Genome-wide survey and expression analysis of the stress-associated protein gene family in desert poplar, Populus euphratica. Tree Genet. Genomes 2016, 12, 78. [Google Scholar] [CrossRef]
  7. Giri, J.; Vij, S.; Dansana, P.K.; Tyagi, A.K. Rice A20/AN1 zinc-finger containing stress-associated proteins (SAP1/11) and a receptor-like cytoplasmic kinase (OsRLCK253) interact via A20 zinc-finger and confer abiotic stress tolerance in transgenic Arabidopsis plants. New Phytol. 2011, 191, 721–732. [Google Scholar] [CrossRef] [PubMed]
  8. Dixit, V.M.; Green, S.; Sarma, V.; Holzman, L.B.; Wolf, F.W.; O’Rourke, K.; Ward, P.A.; Prochownik, E.V.; Marks, R.M. Tumor necrosis factor-alpha induction of novel gene products in human endothelial cells including a macrophage-specifc chemotaxin. J. Biol. Chem. 1990, 265, 2973–2978. [Google Scholar] [PubMed]
  9. Linnen, J.M.; Bailey, C.P.; Weeks, D.L. Two related localized mRNAs from Xenopus laevis encode ubiquitin-like fusion proteins. Gene 1993, 128, 181–188. [Google Scholar] [CrossRef]
  10. Solanke, A.U.; Sharma, M.K.; Tyagi, A.K.; Sharma, A.K. Characterization and phylogenetic analysis of environmental stress-responsive SAP gene family encoding A20/AN1 zinc fnger proteins in tomato. Mol. Genet. Genom. 2009, 282, 153–164. [Google Scholar] [CrossRef] [PubMed]
  11. Gao, W.; Long, L.; Tian, X.; Jin, J.; Liu, H.; Zhang, H.; Xu, F.; Song, C. Genome-wide identification and expression analysis of stress-associated proteins (SAPs) containing A20/AN1 zinc finger in cotton. Mol. Genet. Genom. 2016, 291, 2199–2213. [Google Scholar] [CrossRef] [PubMed]
  12. Huang, J.; Wang, M.M.; Jiang, Y.; Bao, Y.M.; Huang, X.; Sun, H.; Xu, D.Q.; Lan, H.X.; Zhang, H.S. Expression analysis of rice A20/AN1-type zinc fnger genes and characterization of ZFP177 that contributes to temperature stress tolerance. Gene 2008, 420, 135–144. [Google Scholar] [CrossRef] [PubMed]
  13. Kang, M.; Fokar, M.; Abdelmageed, H.; Allen, R.D. Arabidopsis SAP5 functions as a positive regulator of stress responses and exhibits E3 ubiquitin ligase activity. Plant Mol. Biol. 2011, 75, 451–466. [Google Scholar] [CrossRef] [PubMed]
  14. Xuan, N.; Jin, Y.; Zhang, H.; Xie, Y.; Liu, Y.; Wang, G. A putative maize zinc-finger protein gene, ZmAN13, participates in abiotic stress response. Plant Cell Tissue Organ Cult. 2011, 107, 101–112. [Google Scholar] [CrossRef]
  15. Gimeno-Gilles, C.; Gervais, M.L.; Planchet, E.; Satour, P.; Limami, A.M.; Lelievre, E. A stress-associated protein containing A20/AN1 zing-finger domains expressed in Medicago truncatula seeds. Plant Physiol. Biochem. 2011, 49, 303–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Sreedharan, S.; Shekhawat, U.K.; Ganapathi, T.R. MusaSAP1, a A20/AN1 zinc finger gene from banana functions as a positive regulator indifferent stress responses. Plant Mol. Biol. 2012, 80, 503–517. [Google Scholar] [CrossRef] [PubMed]
  17. Ben Saad, R.; Fabre, D.; Mieulet, D.; Meynard, D.; Dingkuhn, M.; Al-Doss, A.; Guiderdoni, E.; Hassairi, A. Expression of the Aeluropus littoralis AlSAP gene in rice confers broad tolerance to abiotic stresses through maintenance of photosynthesis. Plant Cell Environ. 2012, 35, 626–643. [Google Scholar] [CrossRef] [PubMed]
  18. Yoon, S.K.; Bae, E.K.; Lee, H.; Choi, Y.I.; Han, M.; Choi, H.; Kang, K.S.; Park, E.J. Downregulation of stress-associated protein 1 (PagSAP1) increases salt stress tolerance in poplar (Populus alba × P. glandulosa). Trees 2018, 32, 823–833. [Google Scholar] [CrossRef]
  19. Tyagi, H.; Jha, S.; Sharma, M.; Giri, J.; Tyagi, A.K. Rice SAPs are responsive to multiple biotic stresses and overexpression of OsSAP1, an A20/AN1 zinc-fnger protein, enhances the basal resistance against pathogen infection in tobacco. Plant Sci. 2014, 225, 68–76. [Google Scholar] [CrossRef] [PubMed]
  20. Liu, Y.; Xu, Y.; Xiao, J.; Ma, Q.; Li, D.; Xue, Z.; Chong, K. OsDOG, a gibberellin-induced A20/AN1 zinc-fnger protein, negatively regulates gibberellin-mediated cell elongation in rice. J. Plant Physiol. 2011, 168, 1098–1105. [Google Scholar] [CrossRef] [PubMed]
  21. Zhang, Y.; Lan, H.; Shao, Q.; Wang, R.; Chen, H.; Tang, H.; Zhang, H.; Huang, J. An A20/AN1-type zinc finger protein modulates gibberellins and abscisic acid contents and increases sensitivity to abiotic stress in rice (Oryza sativa). J. Exp. Bot. 2016, 67, 315–326. [Google Scholar] [CrossRef] [PubMed]
  22. Ströher, E.; Wang, X.J.; Roloff, N.; Klein, P.; Husemann, A.; Dietz, K.J. Redox-dependent regulation of the stress-induced zinc-finger protein SAP12 in Arabidopsis thaliana. Mol. Plant 2009, 2, 357–367. [Google Scholar] [CrossRef] [PubMed]
  23. Kothari, K.S.; Dansana, P.K.; Giri, J.; Tyagi, A.K. Rice stress associated protein 1 (OsSAP1) interacts with aminotransferase (OsAMTR1) and pathogenesis-related 1a protein (OsSCP) and regulates abiotic stress responses. Front. Plant Sci. 2016, 7, 1057. [Google Scholar] [CrossRef] [PubMed]
  24. Vij, S.; Tyagi, A.K. Genome-wide analysis of the stress associated protein (SAP) gene family containing A20/AN1 zinc-fnger(s) in rice and their phylogenetic relationship with Arabidopsis. Mol. Genet. Genom. 2006, 276, 565–575. [Google Scholar] [CrossRef] [PubMed]
  25. Zhou, Y.; Zeng, L.; Chen, R.; Wang, Y.; Song, J. Genome-wide identification and characterization of stress-associated protein (SAP) gene family encoding A20/AN1 zinc-finger proteins in Medicago truncatula. Arch. Biol. Sci. 2018, 70, 87–98. [Google Scholar] [CrossRef]
  26. Velasco, R.; Zharkikh, A.; Affourtit, J.; Dhingra, A.; Cestaro, A.; Kalyanaraman, A.; Fontana, P.; Bhatnagar, S.K.; Troggio, M.; Pruss, D.; et al. The genome of the domesticated apple (Malus × domestica Borkh.). Nat. Genet. 2010, 42, 833–839. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Li, X.; Ling, K.; Zhang, J.; Xie, Y.; Wang, L.; Yan, Y.; Wang, N.; Xu, J.; Li, J.; Li, C.; et al. Improved hybridde novogenome assembly of domesticated apple (Malus x domestica). Gigascience 2016, 5, 35. [Google Scholar] [CrossRef] [PubMed]
  28. Daccord, N.; Celton, J.M.; Linsmith, G.; Becker, C.; Choisne, N.; Schijlen, E.; van de Geest, H.; Bianco, L.; Micheletti, D.; Velasco, R.; et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat. Genet. 2017, 49, 1099–1106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Zhao, T.; Liang, D.; Wang, P.; Liu, J.; Ma, F. Genome-wide analysis and expression profiling of the DREB transcription factor gene family in Malus under abiotic stress. Mol. Genet. Genom. 2012, 287, 423–436. [Google Scholar] [CrossRef] [PubMed]
  30. Cao, Z.; Zhang, S.; Wang, R.; Zhang, R.; Hao, Y. Genome wide analysis of the apple MYB transcription factor family allows the identification of MdoMYB121 gene conferring abiotic stress tolerance in plants. PLoS ONE 2013, 8, e69955. [Google Scholar]
  31. Tian, Y.; Dong, Q.; Ji, Z.; Chi, F.; Cong, P.; Zhou, Z. Genome-wide identification and analysis of the MADS-box gene family in apple. Gene 2015, 555, 277–290. [Google Scholar] [CrossRef] [PubMed]
  32. Gu, Y.; Ji, Z.; Chi, F.; Qiao, Z.; Xu, C.; Zhang, J.; Dong, Q.; Zhou, Z. Bioinformatics and expression analysis of the WRKY gene family in apple. Sci. Agric. Sin. 2015, 48, 3221–3238. [Google Scholar]
  33. Meng, D.; Li, Y.; Bai, Y.; Li, M.; Cheng, L. Genome-wide identification and characterization of WRKY transcriptional factor family in apple and analysis of their responses to waterlogging and drought stress. Plant Physiol. Biochem. 2016, 103, 71–83. [Google Scholar] [CrossRef] [PubMed]
  34. Dong, Q.; Zhao, S.; Duan, D.; Tian, Y.; Wang, Y.; Mao, K.; Zhou, Z.; Ma, F. Structural and functional analyses of genes encoding VQ proteins in apple. Plant Sci. 2018, 272, 208–219. [Google Scholar] [CrossRef] [PubMed]
  35. Jeffares, D.C.; Penkett, C.J.; Bahler, J. Rapidly regulated genes are intron poor. Trends Genet. 2008, 24, 375–378. [Google Scholar] [CrossRef] [PubMed]
  36. Kanneganti, V.; Gupta, A.K. Overexpression of OsiSAP8, a member of stress associated protein (SAP) gene family of rice confers tolerance to salt, drought and cold stress in transgenic tobacco and rice. Plant Mol. Biol. 2008, 66, 445–462. [Google Scholar] [CrossRef] [PubMed]
  37. Ben Saad, R.; Zouari, N.; Ben Ramdhan, W.; Azaza, J.; Meynard, D.; Guiderdoni, E.; Hassairi, A. Improved drought and salt stress tolerance in transgenic tobacco overexpressing a novel A20/AN1 zinc-finger “AlSAP” gene isolated from the halophyte grass Aeluropus littoralis. Plant Mol. Biol. 2010, 72, 171–190. [Google Scholar] [CrossRef] [PubMed]
  38. Kang, M.; Lee, S.; Abdelmageed, H.; Reichert, A.; Lee, H.K.; Fokar, M.; Mysore, K.S.; Allen, R.D. Arabidopsis stress associated protein 9 mediates biotic and abiotic stress responsive ABA signaling via the proteasome pathway. Plant Cell Environ. 2017, 40, 702–716. [Google Scholar] [CrossRef] [PubMed]
  39. Dixit, A.; Tomar, P.; Vaine, E.; Abdullah, H.; Hazen, S.; Dhankher, O.P. A stress-associated protein, AtSAP13, from Arabidopsis thaliana provides tolerance to multiple abiotic stresses. Plant Cell Environ. 2018, 41, 1171–1185. [Google Scholar] [CrossRef] [PubMed]
  40. Dansana, P.K.; Kothari, K.S.; Vij, S.; Tyagi, A.K. OsiSAP1 overexpression improves water-deficit stress tolerance in transgenic rice by affecting expression of endogenous stress-related genes. Plant Cell Rep. 2014, 33, 1425–1440. [Google Scholar] [CrossRef] [PubMed]
  41. Hozain, M.; Abdelmageed, H.; Lee, J.; Kang, M.; Fokar, M.; Allen, R.D.; Holaday, A.S. Expression of AtSAP5 in cotton up-regulates putative stress-responsive genes and improves the tolerance to rapidly developing water deficit and moderate heat stress. J. Plant Physiol. 2012, 169, 1261–1270. [Google Scholar] [CrossRef] [PubMed]
  42. Wei, Z.; Gao, T.; Liang, B.; Zhao, Q.; Ma, F.; Li, C. Effects of exogenous melatonin on methyl viologen-mediated oxidative stress in apple leaf. Int. J. Mol. Sci. 2018, 19, 316. [Google Scholar]
  43. Dong, Q.L.; Liu, D.D.; An, X.H.; Hu, D.G.; Yao, Y.X.; Hao, Y.J. MdVHP1 encodes an apple vacuolar H+-PPase and enhances stress tolerance in transgenic apple callus and tomato. J. Plant Physiol. 2011, 168, 2124–2133. [Google Scholar] [CrossRef] [PubMed]
  44. Tu, M.; Wang, X.; Feng, T.; Sun, X.; Wang, Y.; Huang, L.; Gao, M.; Wang, Y.; Wang, X. Expression of a grape (Vitis vinifera) bZIP transcription factor, VlbZIP36, in Arabidopsis thaliana, confers tolerance of drought stress during seed germination and seedling establishment. Plant Sci. 2016, 252, 311–323. [Google Scholar] [CrossRef] [PubMed]
  45. Zhou, K.; Hu, L.; Li, P.; Gong, X.; Ma, F. Genome-wide identification of glycosyltransferases converting phloretin to phloridzin in malus species. Plant Sci. 2017, 265, 131–145. [Google Scholar] [CrossRef] [PubMed]
  46. Liu, Y.; Guan, X.; Liu, S.; Yang, M.; Ren, J.; Guo, M.; Huang, Z.; Zhang, Y. Genome-wide identification and analysis of tcp transcription factors involved in the formation of leafy head in Chinese cabbage. Int. J. Mol. Sci. 2018, 19, 847. [Google Scholar] [CrossRef] [PubMed]
  47. Gu, Y.B.; Ji, Z.R.; Chi, F.M.; Qiao, Z.; Xu, C.N.; Zhang, J.X.; Zhou, Z.S.; Dong, Q.L. Genome-wide identification and expression analysis of the WRKY gene family in peach. Hereditas 2016, 38, 254–270. [Google Scholar] [PubMed]
  48. Proost, S.; Van Bel, M.; Vaneechoutte, D.; Van de Peer, Y.; Inzé, D.; Mueller-Roeber, B.; Vandepoele, K. PLAZA 3.0: An access point for plant comparative genomics. Nucleic Acids Res. 2015, 43, 974–981. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, Q.; Dang, H.; Chen, Z.; Wu, J.; Chen, Y.; Chen, S.; Luo, L. Genome-wide identification, expression, and functional analysis of the sugar transporter gene family in cassava (manihot esculenta). Int. J. Mol. Sci. 2018, 19, 987. [Google Scholar] [CrossRef] [PubMed]
  50. Mao, K.; Dong, Q.; Li, C.; Liu, C.; Ma, F. Genome wide identification and characterization of apple bHLH transcription factors and expression analysis in response to drought and salt stress. Front. Plant Sci. 2017, 8, 480. [Google Scholar] [CrossRef] [PubMed]
  51. Wang, P.; Sun, X.; Jia, X.; Ma, F. Apple autophagy-related protein MdATG3s afford tolerance to multiple abiotic stresses. Plant Sci. 2017, 256, 53–64. [Google Scholar] [CrossRef] [PubMed]
  52. Shao, Y.; Qin, Y.; Zou, Y.; Ma, F. Genome-wide identification and expression profiling of the SnRK2 gene family in Malus prunifolia. Gene 2014, 552, 87–97. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, N.; Guo, T.; Sun, X.; Jia, X.; Wang, P.; Shao, Y.; Liang, B.; Gong, X.; Ma, F. Functions of two Malus hupehensis (Pamp.) Rehd. YTPs (MhYTP1 and MhYTP2) in biotic- and abiotic-stress responses. Plant Sci. 2017, 261, 18–27. [Google Scholar] [PubMed]
  54. Perini, P.; Pasquali, G.; Margis-Pinheiro, M.; de Oliviera, P.R.D.; Revers, L.F. Reference genes for transcriptional analysis of flowering and fruit ripening stages in apple (Malus × domestica Borkh.). Mol. Breed. 2014, 34, 829–842. [Google Scholar] [CrossRef]
  55. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  56. Tan, Y.; Li, M.; Yang, Y.; Sun, X.; Wang, N.; Liang, B.; Ma, F. Overexpression of MpCYS4, a phytocystatin gene from Malus prunifolia (Willd.) Borkh., enhances stomatal closure to confer drought tolerance in transgenic Arabidopsis and apple. Front. Plant Sci. 2017, 8, 33. [Google Scholar] [CrossRef] [PubMed]
  57. Liang, B.; Li, C.; Ma, C.; Wei, Z.; Wang, Q.; Huang, D.; Chen, Q.; Li, C.; Ma, F. Dopamine alleviates nutrient deficiency-induced stress in malus hupehensis. Plant Physiol. Biochem. 2017, 119, 346–359. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02478 g001 550
Figure 1. (A) Phylogenetic relationships; (B) Structures for 30 genes; and (C) Analysis of conserved domains for stress-associated protein (SAP) genes in apple. A phylogenetic tree for full-length amino acid sequences was constructed with MEGA software and the NJ method.
Figure 1. (A) Phylogenetic relationships; (B) Structures for 30 genes; and (C) Analysis of conserved domains for stress-associated protein (SAP) genes in apple. A phylogenetic tree for full-length amino acid sequences was constructed with MEGA software and the NJ method.
Ijms 19 02478 g001
Ijms 19 02478 g002 550
Figure 2. Multiple alignments of A20/AN1 domain and C2H2 amino acids in apple SAPs, using the DNAMAN program. Conserved domains are boxed, and identical amino acids are shown against a dark blue background (Similarity: dark blue = 100%; pink > 75%; cyan > 50%).
Figure 2. Multiple alignments of A20/AN1 domain and C2H2 amino acids in apple SAPs, using the DNAMAN program. Conserved domains are boxed, and identical amino acids are shown against a dark blue background (Similarity: dark blue = 100%; pink > 75%; cyan > 50%).
Ijms 19 02478 g002
Ijms 19 02478 g003 550
Figure 3. (A) Sequence logos for A20 domain in 26 MdSAP proteins, generated via WebLogo; (B) Three-dimensional tertiary structural model of A20 domain (PDB ID: 2KZY.1.A); (C) Sequence logos of AN1 domain in 28 MdSAP proteins; and (D) Three-dimensional tertiary structural model of AN1 domain (PDB ID: 1WFP.1.A). Within each stack, symbol height indicates the relative frequency of each amino acid at that position. Logos for the A20 domain and AN1 domain were obtained through multiple alignments of 26 and 28 MdSAP protein sequences, respectively. At the top of the corresponding amino acid sequences, arabic numbers (1–5) indicate β-sheets in A20 and AN1 domains. In (B,D), α-helices are red, β-sheets (arabic numbers 1–5) are yellow, and strands are blue/gray. Three-dimensional representations were generated with RasTop software.
Figure 3. (A) Sequence logos for A20 domain in 26 MdSAP proteins, generated via WebLogo; (B) Three-dimensional tertiary structural model of A20 domain (PDB ID: 2KZY.1.A); (C) Sequence logos of AN1 domain in 28 MdSAP proteins; and (D) Three-dimensional tertiary structural model of AN1 domain (PDB ID: 1WFP.1.A). Within each stack, symbol height indicates the relative frequency of each amino acid at that position. Logos for the A20 domain and AN1 domain were obtained through multiple alignments of 26 and 28 MdSAP protein sequences, respectively. At the top of the corresponding amino acid sequences, arabic numbers (1–5) indicate β-sheets in A20 and AN1 domains. In (B,D), α-helices are red, β-sheets (arabic numbers 1–5) are yellow, and strands are blue/gray. Three-dimensional representations were generated with RasTop software.
Ijms 19 02478 g003
Ijms 19 02478 g004 550
Figure 4. Phylogenetic analysis of 453 SAP proteins in 32 species. Unrooted NJ tree was constructed with MEGA 6 software, using full-length amino acid sequences. Tree comprises six subgroups (Ia–Id; IIa and IIb).
Figure 4. Phylogenetic analysis of 453 SAP proteins in 32 species. Unrooted NJ tree was constructed with MEGA 6 software, using full-length amino acid sequences. Tree comprises six subgroups (Ia–Id; IIa and IIb).
Ijms 19 02478 g004
Ijms 19 02478 g005 550
Figure 5. Chromosomal locations of 28 apple SAP genes. Scale is in megabases (Mb). Red font, tandem duplication; blue font, segmental duplication; grey area, genome-wide duplications.
Figure 5. Chromosomal locations of 28 apple SAP genes. Scale is in megabases (Mb). Red font, tandem duplication; blue font, segmental duplication; grey area, genome-wide duplications.
Ijms 19 02478 g005
Ijms 19 02478 g006a 550Ijms 19 02478 g006b 550
Figure 6. Tissue-specific expression and drought response of 13 MdSAP genes. (A) Expression patterns of 13 MdSAP genes in apple tissues. (B) Expression patterns of 13 MdSAP genes in response to drought stress, i.e., normally scheduled irrigation withheld for 0, 4, or 8 days. Three independent replicates were used for calculations. Error bars indicate standard deviation.
Figure 6. Tissue-specific expression and drought response of 13 MdSAP genes. (A) Expression patterns of 13 MdSAP genes in apple tissues. (B) Expression patterns of 13 MdSAP genes in response to drought stress, i.e., normally scheduled irrigation withheld for 0, 4, or 8 days. Three independent replicates were used for calculations. Error bars indicate standard deviation.
Ijms 19 02478 g006aIjms 19 02478 g006b
Ijms 19 02478 g007 550
Figure 7. Overexpression of MdSAP15 in Arabidopsis increased osmotic tolerance in response to mannitol treatment. (A) PCR identification and relative expression analysis of MdSAP15 in Col-0 wild-type (Col) and transgenic Arabidopsis lines. M, DNA marker; −, negative control (H2O); +, positive control (plasmid DNA of 35S::MdSAP15 pRI101AN vector). Specific primers for MdSAP15 were used to detect relative expression levels of Col and transgenic Arabidopsis lines. (B) Representative images of “Col” and transgenic seedlings at 5 days after seeds had been cultivated for 11 days on MS medium alone (0 mM mannitol) or MS medium supplemented with 200 or 300 mM of mannitol. Bars = 1 cm. (C) Primary root lengths and (D) fresh weights of “Col” and transgenic seedlings measured on Day 5 after plants had been exposed for 11 days to osmotic stress. (E) Relative electrolyte leakage, and levels of chlorophyll (F), proline (G), and malondialdehyde (MDA) (H) in “Col” and transgenic seedlings, measured on Day 5 following treatment on MS medium with 0 mM, 200 mM, or 300 mM of mannitol for 11 d. Error bars represent SD based on three independent replicates. For (C–H), bars not labeled with same letters in each panel indicate values are significantly different at p < 0.05, based on one-way ANOVA and Duncan’s tests.
Figure 7. Overexpression of MdSAP15 in Arabidopsis increased osmotic tolerance in response to mannitol treatment. (A) PCR identification and relative expression analysis of MdSAP15 in Col-0 wild-type (Col) and transgenic Arabidopsis lines. M, DNA marker; −, negative control (H2O); +, positive control (plasmid DNA of 35S::MdSAP15 pRI101AN vector). Specific primers for MdSAP15 were used to detect relative expression levels of Col and transgenic Arabidopsis lines. (B) Representative images of “Col” and transgenic seedlings at 5 days after seeds had been cultivated for 11 days on MS medium alone (0 mM mannitol) or MS medium supplemented with 200 or 300 mM of mannitol. Bars = 1 cm. (C) Primary root lengths and (D) fresh weights of “Col” and transgenic seedlings measured on Day 5 after plants had been exposed for 11 days to osmotic stress. (E) Relative electrolyte leakage, and levels of chlorophyll (F), proline (G), and malondialdehyde (MDA) (H) in “Col” and transgenic seedlings, measured on Day 5 following treatment on MS medium with 0 mM, 200 mM, or 300 mM of mannitol for 11 d. Error bars represent SD based on three independent replicates. For (C–H), bars not labeled with same letters in each panel indicate values are significantly different at p < 0.05, based on one-way ANOVA and Duncan’s tests.
Ijms 19 02478 g007
Ijms 19 02478 g008 550
Figure 8. Assessment of drought tolerance in “Col” and transgenic Arabidopsis plants. (A) Representative images of four-week-old plants; (B) Survival rates of WT and transgenic lines 2 days after rehydration began. Error bars represent SD based on three independent replicates. For (B), bars not labeled with same letters in each panel indicate values are significantly different at p < 0.05, based on one-way ANOVA and Duncan’s tests.
Figure 8. Assessment of drought tolerance in “Col” and transgenic Arabidopsis plants. (A) Representative images of four-week-old plants; (B) Survival rates of WT and transgenic lines 2 days after rehydration began. Error bars represent SD based on three independent replicates. For (B), bars not labeled with same letters in each panel indicate values are significantly different at p < 0.05, based on one-way ANOVA and Duncan’s tests.
Ijms 19 02478 g008
Table
Table 1. Properties of SAPs identified from apple genome.
Table 1. Properties of SAPs identified from apple genome.
Gene NameGene ID 1Zinc Finger DomainProtein Length (aa)Molecular Weight (kDa)Theoretical Isoeletrical PointChromosome Location
MdSAP1MDP0000494946A20-AN115116.538.41chr1:24227744..24228199
MdSAP2MDP0000588934A20-AN116217.718.62chr2:22718230..22718718
MdSAP3MDP0000122842A20-AN115617.078.27chr2:22745813..22746283
MdSAP4MDP0000237812A20-AN115617.078.27chr2:22753663..22754133
MdSAP5MDP0000316313A2016118.118.95chr2:24190887..24193156
MdSAP6MDP0000543745A2013615.186.37chr2:25865864..25866277
MdSAP7MDP0000362676A20-AN116918.429.05chr2:35848074..35848583
MdSAP8MDP0000874708A20-AN117218.248.12chr2:35871562..35872080
MdSAP9MDP0000362677A20-AN116918.429.05chr2:35878425..35878935
MdSAP10MDP0000164222A20-AN117218.417.52chr3:2250369..2250887
MdSAP11MDP0000516205A20-AN118920.096.78chr4:8801881..8802450
MdSAP12MDP0000506127A20-AN118920.096.78chr4:8807038..8807607
MdSAP13MDP0000286185A20-AN117618.977.46chr6:7080970..7081500
MdSAP14MDP0000263150A20-AN116818.519.16chr7:613109..613615
MdSAP15MDP0000292844A20-AN117118.178.12chr7:697510..698025
MdSAP16MDP0000294781A20-AN117618.977.46chr7:10610840..10611370
MdSAP17MDP0000133254A20-AN115016.418.41chr7:22721788..22722240
MdSAP18MDP0000139359A20-AN117018.598.67chr8:16795514..16796026
MdSAP19MDP0000707978A20-AN117018.598.67chr8:16801945..16802457
MdSAP20MDP0000296953A20-AN117018.598.67chr8:16804379..16804891
MdSAP21MDP0000211516AN1-AN119321.368.67chr9:2267267..2267848
MdSAP22MDP0000165407A2020923.098.29chr9:27693011..27693884
MdSAP23MDP0000231017A20-AN117218.187.51chr11:2865175..2865693
MdSAP24MDP0000683912A20-AN117218.187.51chr11:2875969..2876487
MdSAP25MDP0000652898AN1-AN1-C2H2-C2H229332.048.26chr12:7214047..7215925
MdSAP26MDP0000086327A20-AN117418.787.46chr14:7829264..7829785
MdSAP27MDP0000141121AN1-AN119221.398.72chr17:2716099..2716758
MdSAP28MDP0000853499AN1-AN119221.398.72chr17:2716210..2716869
MdSAP29MDP0000661416A20-AN120321.867.66unanchored:14381067..14381678
MdSAP30MDP0000284856A20-AN120321.867.66unanchored:14408647..14409258
1 Gene ID in apple genome (https://www.rosaceae.org/gb/gbrowse/malus_x_domestica/).
Table
Table 2. The cis-acting elements of 13 promoters in apple SAP genes.
Table 2. The cis-acting elements of 13 promoters in apple SAP genes.
Cis-Acting ElementsABREARECGTCAEREHSELTRMBSTCATC-Rich RepeatW-Box
Stress to ResponseABAHypoxiaMeJAEthyleneHeatChillingDroughtSADefensePathogen
MdSAP72/22/02/2 1/01/1 0/20/2
MdSAP8 0/1 0/1 3/00/3 0/1
MdSAP102/00/20/1 1/00/2
MdSAP12 1/0 1/0
MdSAP14 2/01/1 2/01/0 1/01/00/1
MdSAP150/10/2 0/1 0/10/11/01/0
MdSAP16 2/2 2/0 1/03/0
MdSAP19 0/2 0/1 1/1 1/12/0
MdSAP21 1/01/0 1/2 1/0
MdSAP231/01/40/2 1/0 1/0
MdSAP251/11/11/1 0/21/10/10/1
MdSAP28 0/41/2 1/11/11/1
MdSAP29 0/21/2 1/11/1 1/02/0
ABA: abscisic acid; ABRE: ABA response element; ARE: anaerobic response element; CGTCA: MeJA-responsiveness; ERE: ethylene-responsive element; HSE: heat shock response element; LTR: low-temperature response element; MBS: MYB binding site involved in drought response; TCA: salicylic acid response element; TC-Rich Repeat: defense and stress responsiveness; W-box: elicitation, wounding, and pathogen responsiveness/binding site of WRKY type transcription factors.
Table
Table 3. Numbers of SAP gene family members in various species.
Table 3. Numbers of SAP gene family members in various species.
Plant SpeciesA20-A20-A20AN1AN1-AN1-AN1-Total Number
AN1A20-AN1AN1-AN1-
AN1 C2H2C2H2-
C2H2
Malus domestica2303030130
Arabidopsis thaliana1000111114
Oryza sativa1111310118
Populus trichocarpa1500210119
Solanum lycopersicum900120113
Gossypium hirsutum1400220119
Populus euphratica1500020118
Arabidopsis lyrata1200011115
Amborella trichopoda30111017
Brassica rapa1801512128
Beta vulgaris60010018
Citrullus lanatus700210111
Cucumis melo1000010112
Carica papaya50001017
Chlamydomonas reinhardtii10010002
Capsella rubella1000011113
Citrus sinensis1000010112
Eucalyptus grandis810110011
Fragaria vesca1200110115
Glycine max1800220426
Gossypium raimondii1400210219
Lotus japonicus40011006
Manihot esculenta1400110117
Medicago truncatula1100210216
Physcomitrella patens600120110
Prunus persica800010211
Ricinus communis50021019
Solanum tuberosum1300310219
Theobroma cacao1000010112
Thellungiella parvula1100011114
Vitis vinifera401420011
Zea mays800110111
Table
Table 4. Statistics for numbers of intronless members within different groups of the SAP gene family.
Table 4. Statistics for numbers of intronless members within different groups of the SAP gene family.
GroupIaIbIcIdIIaIIb
Intronless number10543636115
Total number1235374783934
Percentage (%)85.3681.1385.1378.22.5614.7
Table
Table 5. Application of primers and sequences. ORF: open reading frames.
Table 5. Application of primers and sequences. ORF: open reading frames.
UsePrimer NameForward Primer (5′–3′)Reverse Primer (5′–3′)
CompleteMdSAP7ATGAAAAAAATGGCACAGAGAATCAAACCCGGACGATCTTTGCGG
ORFMdSAP8ATGGAGCACAATGAGACAGGATTCAGATTTTATCCAGCTTTTCT
amplificationMdSAP10ATGGAGCACGAGGAGACTGGATGTTAGATTTTATCGAGCTTCTCA
MdSAP12ATGGCGGAAGAGCACAGATGCGTCAAATCTTCTCGAGCTTCTCG
MdSAP14ATGAAAAAAATGGCACAGAGAATCAGAGCCGGACGATCTTCGCA
MdSAP15ATGGAGCACAATGAGACAGGATGTCAGATTTTATCCAGCTTGTCTG
MdSAP16ATGGAATCTCATGATGAAACTGCTAGATTTTGTCAAGTTTGTCTG
MdSAP19ATGGCGGAAGAGCATCGTTGCCATCAAATCTTATGCAGCTTCTCCG
MdSAP21ATGATGGGAGGAACAGAAGCTTTCAATACGCTCGAACAGATGGCC
MdSAP23ATGGAGCACGAGGAGACTGGATGTTAGATTTTACCAAGCTTGTCAG
MdSAP25ATGGGAACTCCGGAATTCCCAGACTATGCTCTTGAAGTACCGCCGT
MdSAP28ATGACGGGAGGAACAGAAGCTTTTCAATAAGCTCGAACAGAAGGC
MdSAP29ATGGCGGAAGAGCACAGATGCGATCACGGATTTGTACGTTGGCAA
PromoterMdSAP7CAGATTTTGTTCAAATGTAGGTGGGCGATGGAGGAGACAGAAAT
amplificationMdSAP8TGTTTCAATTGCGTTCTTGAGGCATTGTAATTCGCTTAAGTTCT
MdSAP10ACCTTTTCCAAAACCGTTATTAGTGCGAAAACCAACAATTAATGG
MdSAP12
MdSAP14GTAAAGAGGTTAGTGGCCCTGAACAAATTCTGATCGATCGATCGAT
MdSAP15ATGCGCTTTACTGTTTTTTCAGTCATTGTAATTCGCTAAGTCCTT
MdSAP16CACGAGGAGAGCACTAAAATGGACACCAAGAAAACCTCGCCGTTT
MdSAP19ACCTTTCTTTTGAGAAGTTTGTTGCAATTCCAAAACAAATTATTC
MdSAP21ATGGATTCTAGTTTGATTTGGGCGATTTTTCAGTTTGTTAAATTTT
MdSAP23ATATTTCCATCACATTGAATAACTACTCAGCTTACCTGCAAAGAG
MdSAP25GCAGGTAGAGTTTCAAAGTACGAAATTTTGTATGTACAACACTA
MdSAP28ACAGGTCACCGTGGTGACTCCGGGTCGGTCGGTCGGTCTGGGGTTG
MdSAP29GTGCTTTTTGTTGGAACACAAAGCGATCGAGAGGACAAAAATATTA
qRT-PCRMdSAP7TCGTCCGGGTTTGATGATTTTCCCCGGTCTCTGAATTTCG
MdSAP8GGGAAGCGGATAGGAACCATCTTGGGAGCTTCAGGAGGAG
MdSAP10GATTATCGCACTGCTGGACGAGTGCTAAGATACCGCTGCA
MdSAP12GTTGGTCATAGCCGAGAAGCATCAGCTTAATTCCCACGCG
MdSAP14GCTCTGACCGGTTTGACAATTTGCTGATGATCTCCGGGAG
MdSAP15ATGATTACCGGACTGCTGCTCCACATGGGTAGAAATGAGAGC
MdSAP16GCCAATCCTATCGTGAAGGCGAGACCTATGCAGACAAGAAGC
MdSAP19CGATTTCAGAGGGATGGGGACAACCATCCCCTACCCCAAT
MdSAP21AGGGAAAGAATGCGGGAAGACGAAGAAACATGAAACTGCGG
MdSAP23GCCAACCCTGTCGTAAAAGCTGCTAAGATACCGCTGCAGA
MdSAP25AATCCAATCCAAGCCTCGGATCCCATCCGAATTTTGCACG
MdSAP28TGCTTTGAGGGAAGGGAAAGAACATCGAATTGTGGAAGCAGA
MdSAP29TTCCTCCTCGCACAGATCAGTCCGCCATGTCTACAGTCAA
MdMDHCGTGATTGGGTACTTGGAACTGGCAAGTGACTGGGAATGA
pRI-101ANMdSAP15TTGATACATATGCCCGTCGACATGGAGCACAATAGAGTTGTTGATTCAGGATCCTCAGATTTTATC
35SCGCACAATCCCACTATCCTT
qRT-MdSAP15AGTCGTTGCAGCATCCATTGGGAAGCCTGTGTTGAGATAAGC
AtActin2GTGAAGGCTGGATTTGCAGGAAACCTCCGATCCAGACACTGT

Share and Cite

MDPI and ACS Style

Dong, Q.; Duan, D.; Zhao, S.; Xu, B.; Luo, J.; Wang, Q.; Huang, D.; Liu, C.; Li, C.; Gong, X.; et al. Genome-Wide Analysis and Cloning of the Apple Stress-Associated Protein Gene Family Reveals MdSAP15, Which Confers Tolerance to Drought and Osmotic Stresses in Transgenic Arabidopsis. Int. J. Mol. Sci. 2018, 19, 2478. https://doi.org/10.3390/ijms19092478

AMA Style

Dong Q, Duan D, Zhao S, Xu B, Luo J, Wang Q, Huang D, Liu C, Li C, Gong X, et al. Genome-Wide Analysis and Cloning of the Apple Stress-Associated Protein Gene Family Reveals MdSAP15, Which Confers Tolerance to Drought and Osmotic Stresses in Transgenic Arabidopsis. International Journal of Molecular Sciences. 2018; 19(9):2478. https://doi.org/10.3390/ijms19092478

Chicago/Turabian Style

Dong, Qinglong, Dingyue Duan, Shuang Zhao, Bingyao Xu, Jiawei Luo, Qian Wang, Dong Huang, Changhai Liu, Chao Li, Xiaoqing Gong, and et al. 2018. "Genome-Wide Analysis and Cloning of the Apple Stress-Associated Protein Gene Family Reveals MdSAP15, Which Confers Tolerance to Drought and Osmotic Stresses in Transgenic Arabidopsis" International Journal of Molecular Sciences 19, no. 9: 2478. https://doi.org/10.3390/ijms19092478

APA Style

Dong, Q., Duan, D., Zhao, S., Xu, B., Luo, J., Wang, Q., Huang, D., Liu, C., Li, C., Gong, X., Mao, K., & Ma, F. (2018). Genome-Wide Analysis and Cloning of the Apple Stress-Associated Protein Gene Family Reveals MdSAP15, Which Confers Tolerance to Drought and Osmotic Stresses in Transgenic Arabidopsis. International Journal of Molecular Sciences, 19(9), 2478. https://doi.org/10.3390/ijms19092478

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop