Genome-wide identification of members in the YTH domain-containing RNA-binding protein family in apple and expression analysis of their responsiveness to senescence and abiotic stresses

Genome-Wide Identification and Expression Analysis of YTH Gene Family for Abiotic Stress Regulation in Camellia chekiangoleosa

N6-methyladenosine (m6A) is essential for RNA metabolism in cells. The YTH domain, conserved in the kingdom of Eukaryotes, acts as an m6A reader that binds m6A-containing RNA. In plants, the YTH domain is involved in plant hormone signaling, stress response regulation, RNA stability, translation, and differentiation. However, little is known about the YTH genes in tea-oil tree, which can produce edible oil with high nutritional value. This study aims to identify and characterize the YTH domains within the tea-oil tree (Camellia chekiangoleosa Hu) genome to predict their potential role in development and stress regulation. In this study, 10 members of the YTH family containing the YTH domain named CchYTH1-10 were identified from C. chekiangoleosa. Through analysis of their physical and chemical properties and prediction of subcellular localization, it is known that most family members are located in the nucleus and may have liquid-liquid phase separation. Analysis of cis-acting elements in the CchYTH promoter region revealed that these genes could be closely related to abiotic stress and hormones. The results of expression profiling show that the CchYTH genes were differentially expressed in different tissues, and their expression levels change under drought stress. Overall, these findings could provide a foundation for future research regarding CchYTHs in C. chekiangoleosa and enrich the world in terms of epigenetic mark m6A in forest trees.

The m6A reader MhYTP2 regulates the stability of its target mRNAs contributing to low nitrogen tolerance in apple (Malus domestica)

Studies have shown that the m⁶A reader primarily affects genes expression by participating in the regulation of mRNA localization, splicing, degradation, translation, and other metabolic processes. Previously, we discovered that the apple (Malus domestica) m⁶A reader MhYTP2 bound with and destabilized m⁶A-modified MdMLO19 mRNA. In addition, it enhanced the translation efficiency of m⁶A-modified mRNA of MdGDH1L, encoding a glutamate dehydrogenase, which confers resistance to powdery mildew. In this study, we report the function of MhYTP2 in the regulation of resistance to low nitrogen (N). The overexpression of MhYTP2 enhances the resistance of apple to low N. We show that MhYTP2 binds with and stabilizes the mRNAs of MdALN, which participates in the allantoin catabolic process and cellular response to N starvation in apple; MdPIDL, which participates in root hair elongation; MdTTG1, which is involved in the differentiation process of trichomes; and MdATG8A, which is a core participant in the regulation of autophagy. In addition, MhYTP2 accelerates the degradation of MdRHD3 mRNA, which regulates root development. RNA immunoprecipitation-seq and electrophoretic mobility shift assays show that the mRNAs of MdALN, MdATG8A, MdPIDL, MdTTG1, and MdRHD3 are the direct targets of MhYTP2. Overexpressing or knocking down the above genes in MhYTP2 overexpressing plants dismisses the function of MhYTP2 under low N, suggesting the role of MhYTP2 is dependent on those genes. Together, these results demonstrate that MhYTP2 enhances the resistance of apple to N deficiency by affecting the stability of the bound mRNAs.

Comprehensive Identification and Expression Analysis of the YTH Family of RNA-Binding Proteins in Strawberry

Plant growth and development processes are tightly regulated at multiple levels, including transcriptional and post-transcriptional levels, and the RNA-binding protein YTH regulates gene expression during growth and development at the post-transcriptional level by regulating RNA splicing, processing, stability, and translation. We performed a systematic characterization of YTH genes in diploid forest strawberry and identified a total of nine YTH genes. With the help of phylogenetic analysis, these nine genes were found to belong to two different groups, YTHDC and YTHDF, with YTHDF being further subdivided into three subfamilies. Replication analysis showed that YTH3 and YTH4 are a gene pair generated by tandem repeat replication. These two genes have similarities in gene structure, number of motifs, and distribution patterns. Promoter analysis revealed the presence of multiple developmental, stress response, and hormone-response-related cis-elements. Analysis of available transcriptome data showed that the expression levels of most of the YTH genes were stable with no dramatic changes during development in different tissues. However, YTH3 maintained high expression levels in all tissues and during fruit development, and YTH4 was expressed at higher levels in tissues such as flowers, leaves, and seedlings, while it was significantly lower than YTH3 in white fruits and ripening fruits with little fluctuation. Taken together, our study provides insightful and comprehensive basic information for the study of YTH genes in strawberry.

Genome-Wide Identification and Evolutionary Analysis of Gossypium YTH Domain-Containing RNA-Binding Protein Family and the Role of GhYTH8 in Response to Drought Stress

YTH domain-containing proteins are one kind of RNA-binding protein involved in post-transcriptional regulation and play multiple roles in regulating the growth, development, and abiotic stress responses of plants. However, the YTH domain-containing RNA-binding protein family has not been previously studied in cotton. In this study, a total of 10, 11, 22, and 21 YTH genes were identified in Gossypium arboreum, Gossypium raimondii, Gossypium barbadense, and Gossypium hirsutum, respectively. These Gossypium YTH genes were categorized into three subgroups by phylogenetic analysis. The chromosomal distribution, synteny analysis, structures of Gossypium YTH genes, and the motifs of YTH proteins were analyzed. Furthermore, the cis-element of GhYTH genes promoter, miRNA targets of GhYTH genes, and subcellular localization of GhYTH8 and GhYTH16 were characterized. Expression patterns of GhYTH genes in different tissues, organs, and in response to different stresses were also analyzed. Moreover, functional verifications revealed that silencing GhYTH8 attenuated the drought tolerance in the upland cotton TM-1 line. These findings provide useful clues for the functional and evolutionary analysis of YTH genes in cotton.

The HD-Zip I transcription factor MdHB7-like confers tolerance to salinity in transgenic apple (Malus domestica)

Salinity is a major environmental constraint that substantially limits global agricultural productivity. HD-Zip I transcription factors are involved in plant responses to salt stress, but little is known about the HD-Zip I genes in apple (Malus domestica). Here, we characterized the function of an apple HD-Zip I gene (MdHB7-like) and report that its expression is induced by salt stress. To further explore its role in salt stress, we created MdHB7-like overexpressing and RNAi transgenic apple plants. The overexpression of MdHB7-like improved the photosynthetic performance and reduced ROS and Na⁺ accumulation under salt stress. Plants that overexpressed MdHB7-like also showed increased accumulation of proline and soluble sugars, which may have played an important role in their salt stress tolerance. RNAi suppression of MdHB7-like had the opposite effects. Together, our results demonstrate that MdHB7-like is an important regulator of salt tolerance in apple. Our results provide new insights for future research on the mechanisms by which MdHB7-like promotes salt tolerance and provide a potential target for molecular breeding in apple.

Plant RNA Binding Proteins as Critical Modulators in Drought, High Salinity, Heat, and Cold Stress Responses: An Updated Overview

Plant abiotic stress responses are tightly regulated by different players at multiple levels. At transcriptional or post-transcriptional levels, several RNA binding proteins (RBPs) regulate stress response genes through RNA metabolism. They are increasingly recognized as critical modulators of a myriad of biological processes, including stress responses. Plant RBPs are heterogeneous with one or more conservative RNA motifs that constitute canonical/novel RNA binding domains (RBDs), which can bind to target RNAs to determine their regulation as per the plant requirements at given environmental conditions. Given its biological significance and possible consideration as a potential tool in genetic manipulation programs to improve key agronomic traits amidst frequent episodes of climate anomalies, studies concerning the identification and functional characterization of RBP candidate genes are steadily mounting. This paper presents a comprehensive overview of canonical and novel RBPs and their functions in major abiotic stresses including drought, heat, salt, and cold stress conditions. To some extent, we also briefly describe the basic motif structure of RBPs that would be useful in forthcoming studies. Additionally, we also collected RBP genes that were modulated by stress, but that lacked functional characterization, providing an impetus to conduct further research.

Linking the YTH domain to cancer: the importance of YTH family proteins in epigenetics

N6-methyladenosine (m6A), the most prevalent and reversible modification of mRNA in mammalian cells, has recently been extensively studied in epigenetic regulation. YTH family proteins, whose YTH domain can recognize and bind m6A-containing RNA, are the main "readers" of m6A modification. YTH family proteins perform different functions to determine the metabolic fate of m6A-modified RNA. The crystal structure of the YTH domain has been completely resolved, highlighting the important roles of several conserved residues of the YTH domain in the specific recognition of m6A-modified RNAs. Upstream and downstream targets have been successively revealed in different cancer types and the role of YTH family proteins has been emphasized in m6A research. This review describes the regulation of RNAs by YTH family proteins, the structural features of the YTH domain, and the connections of YTH family proteins with human cancers.

Common bean SNP alleles and candidate genes affecting photosynthesis under contrasting water regimes

Water deficit is a major worldwide constraint to common bean (Phaseolus vulgaris L.) production, being photosynthesis one of the most affected physiological processes. To gain insights into the genetic basis of the photosynthetic response of common bean under water-limited conditions, a collection of 158 Portuguese accessions was grown under both well-watered and water-deficit regimes. Leaf gas-exchange parameters were measured and photosynthetic pigments quantified. The same collection was genotyped using SNP arrays, and SNP-trait associations tested considering a linear mixed model accounting for the genetic relatedness among accessions. A total of 133 SNP-trait associations were identified for net CO2 assimilation rate, transpiration rate, stomatal conductance, and chlorophylls a and b, carotenes, and xanthophyll contents. Ninety of these associations were detected under water-deficit and 43 under well-watered conditions, with only two associations common to both treatments. Identified candidate genes revealed that stomatal regulation, protein translocation across membranes, redox mechanisms, hormone, and osmotic stress signaling were the most relevant processes involved in common bean response to water-limited conditions. These candidates are now preferential targets for common bean water-deficit-tolerance breeding. Additionally, new sources of water-deficit tolerance of Andean, Mesoamerican, and admixed origin were detected as accessions valuable for breeding, and not yet explored.

Genome-wide identification and expression analysis of YTH domain-containing RNA-binding protein family in common wheat

BACKGROUND

N6-Methyladenosine (m6A) is the most widespread RNA modification that plays roles in the regulation of genes and genome stability. YT521-B homology (YTH) domain-containing RNA-binding proteins are important RNA binding proteins that affect the fate of m6A-containing RNA by binding m6A. Little is known about the YTH genes in common wheat (Triticum aestivum L.), one of the most important crops for humans.

RESULTS

A total of 39 TaYTH genes were identified in common wheat, which are comprised of 13 homologous triads, and could be mapped in 18 out of the 21 chromosomes. A phylogenetic analysis revealed that the TaYTHs could be divided into two groups: YTHDF (TaDF) and YTHDC (TaDC). The TaYTHs in the same group share similar motif distributions and domain organizations, which indicates functional similarity between the closely related TaYTHs. The TaDF proteins share only one domain, which is the YTH domain. In contrast, the TaDCs possess three C3H1-type zinc finger repeats at their N-termini in addition to their central YTH domain. In TaDFs, the predicated aromatic cage pocket that binds the methylysine residue of m6A is composed of tryptophan, tryptophan, and tryptophan (WWW). In contrast, the aromatic cage pocket in the TaDCs is composed of tryptophan, tryptophan, and tyrosine (WWY). In addition to the general aspartic acid or asparagine residue used to form a hydrogen bond with N¹ of m6A, histidine might be utilized in some TaDFb proteins. An analysis of the expression using both online RNA-Seq data and quantitative real-time PCR verification revealed that the TaDFa and TaDFb genes are highly expressed in various tissues/organs compared with that of TaDFcs and TaDCs. In addition, the expression of the TaYTH genes is changed in response to various abiotic stresses.

CONCLUSIONS

In this study, we identified 39 TaYTH genes from common wheat. The phylogenetic structure, chromosome distribution, and patterns of expression of these genes and their protein structures were analyzed. Our results provide a foundation for the functional analysis of TaYTHs in the future.

Last-Generation Genome-Environment Associations Reveal the Genetic Basis of Heat Tolerance in Common Bean (Phaseolus vulgaris L.)

Genome-environment associations (GEAs) are a powerful strategy for the study of adaptive traits in wild plant populations, yet they still lack behind in the use of modern statistical methods as the ones suggested for genome-wide association studies (GWASs). In order to bridge this gap, we couple GEA with last-generation GWAS algorithms in common bean to identify novel sources of heat tolerance across naturally heterogeneous ecosystems. Common bean (Phaseolus vulgaris L.) is the most important legume for human consumption, and breeding it for resistance to heat stress is key because annual increases in atmospheric temperature are causing decreases in yield of up to 9% for every 1°C. A total of 78 geo-referenced wild accessions, spanning the two gene pools of common bean, were genotyped by sequencing (GBS), leading to the discovery of 23,373 single-nucleotide polymorphism (SNP) markers. Three indices of heat stress were developed for each accession and inputted in last-generation algorithms (i.e. SUPER, FarmCPU, and BLINK) to identify putative associated loci with the environmental heterogeneity in heat stress. Best-fit models revealed 120 significantly associated alleles distributed in all 11 common bean chromosomes. Flanking candidate genes were identified using 1-kb genomic windows centered in each associated SNP marker. Some of these genes were directly linked to heat-responsive pathways, such as the activation of heat shock proteins (MED23, MED25, HSFB1, HSP40, and HSP20). We also found protein domains related to thermostability in plants such as S1 and Zinc finger A20 and AN1. Other genes were related to biological processes that may correlate with plant tolerance to high temperature, such as time to flowering (MED25, MBD9, and PAP), germination and seedling development (Pkinase_Tyr, Ankyrin-B, and Family Glicosil-hydrolase), cell wall stability (GAE6), and signaling pathway of abiotic stress via abscisic acid (histone-like transcription factors NFYB and phospholipase C) and auxin (Auxin response factor and AUX_IAA). This work offers putative associated loci for marker-assisted and genomic selection for heat tolerance in common bean. It also demonstrates that it is feasible to identify genome-wide environmental associations with modest sample sizes by using a combination of various carefully chosen environmental indices and last-generation GWAS algorithms.

New candidate loci and marker genes on chromosome 7 for improved chilling tolerance in sorghum

Sorghum is often exposed to suboptimal low temperature stress under field conditions, particularly at the seedling establishment stage. Enhancing chilling tolerance will facilitate earlier planting and so minimize the negative impacts of other stresses experienced at later growth stages. Genome-wide association mapping was performed on a sorghum association panel grown under control (30/20 °C; day/night) and chilling (20/10 °C) conditions. Genomic regions on chromosome 7, controlling the emergence index and seedling (root and shoot) vigor, were associated with increased chilling tolerance but they did not co-localize with undesirable tannin content quantitative trait loci (QTLs). Shoot and root samples from highly contrasting haplotype pairs expressing differential responses to chilling stress were used to identify candidate genes. Three candidate genes (an alpha/beta hydrolase domain protein, a DnaJ/Hsp40 motif-containing protein, and a YTH domain-containing RNA-binding protein) were expressed at significantly higher levels under chilling stress in the tolerant haplotype compared with the sensitive haplotype and BTx623. Moreover, two CBF/DREB1A transcription factors on chromosome 2 showed a divergent response to chilling in the contrasting haplotypes. These studies identify haplotype differences on chromosome 7 that modulate chilling tolerance by either regulating CBF or feeding back into this signaling pathway. We have identified new candidate genes that will be useful markers in ongoing efforts to develop tannin-free chilling-tolerant sorghum hybrids.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

Overexpression of the RNA binding protein MhYTP1 in transgenic apple enhances drought tolerance and WUE by improving ABA level under drought condition

MhYTP1 is involved in post-transcriptional regulation as a member of YT521-homology (YTH) domain-containing RNA-binding proteins. We previously cloned MhYTP1 and found it participated in various biotic and abiotic stress responses. However, its function in long-term moderate drought has not been verified. Thus, we explored its biological role in response to drought. Under drought condition, the net photosynthesis rate (P_n) and water use efficiency (WUE) were significantly elevated in MhYTP1-overexpressing (OE) apple plants when compared with the non-transgenic (NT) controls. Further analysis indicated MhYTP1 expression was associated with elevated ABA content, increased stomatal density and reduced stomatal aperture. In addition, to gain insight into the function of stem-specific expression of MhYTP1, grafting experiments were performed. Interestingly, lower transpiration rate (T_r) and higher WUE were observed when transgenic plants were used as scions as opposed to rootstocks and when transgenic rather than NT plants were used as rootstocks, indicating MhYTP1 plays crucial roles in grafted plants. These results define a function for MhYTP1 in promoting tolerance to drought conditions, and suggest that MhYTP1 can serve as a candidate gene for future apple drought resistance breeding with the help of biotechnology.

Comprehensive genomic analysis of the TYROSINE AMINOTRANSFERASE (TAT) genes in apple (Malus domestica) allows the identification of MdTAT2 conferring tolerance to drought and osmotic stresses in plants

Tyrosine aminotransferase (TAT, EC 2.6.1.5) is the first key enzyme that catalyzes the reversible interconversion of tyrosine and 4-hydroxyphenylpyruvate in the tyrosine-derived pathway for syntheses of important secondary metabolites and compounds. Although plant TAT genes have been proposed to be important in response to abiotic stress, there is little information about TAT genes in woody perennial tree species, especially in economic fruit trees. Based on TAT domain searching, sequence homology screening and phylogenetic analysis, we identified four TATs in apple genome. Then, we carried out a detailed phylogenetic analysis of TAT genes from multi-species, focusing on apple (Malus domestica). The result showed that the TAT family comprises three major classes corresponding to genes from angiosperms, mammals, and bacteria. Angiosperm TAT genes could be further divided into six subclasses. Analysis of intron-exon structure revealed that the typical TAT gene contains six introns and seven exons, with exons of similar size at each exon location. Promoter analysis showed that the 5'-flanking region of apple MdTATs contain multiple cis-acting elements including those implicated in light, biotic stress, abiotic stress, and hormone response. MdTATs were expressed to various levels in all apple structures and organs evaluated, and showed distinct expression patterns under water deficit stress. Ectopic expression of MdTAT2 in Arabidopsis or over-expression of MdTAT2 in apple callus tissue conferred enhanced tolerance to drought and osmotic stress. Collectively, these results suggest a role for TAT genes in drought and osmotic stresses and provide valuable information for further research of TAT genes and their function in plants.

Emerging role of dynamic RNA modifications during animal development

The central dogma of molecular biology statically says that the information flows from DNA to messenger RNA to protein. But the recent advances in mass spectrometry and high throughput technology have helped the scientists to view RNA as little more than a courier of genetic information encoded in the DNA. The dynamics of RNA modifications in coding and non-coding RNAs are just emerging as a carrier of non-genetic information, uncovering a new layer of complexity in the regulation of gene expression and protein translation. In this review, we summarize about the current knowledge of N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C) and pseudouridine (Ψ) modifications in RNA, and described how these RNA modifications are implicated in early animal development and in several human diseases.

Genome-wide identification and expression analysis of YTH domain-containing RNA-binding protein family in cucumber (Cucumis sativus)

YTH domain-containing RNA-binding proteins are involved in post-transcriptional regulation and play important roles in the growth and development as well as abiotic stress responses of plants. However, YTH genes have not been previously studied in cucumber (Cucumis sativus). In this study, a total of five YTH genes (CsYTH1-CsYTH5) were identified in cucumber, which could be mapped on three out of the seven cucumber chromosomes. All CsYTH proteins had highly conserved C-terminal YTH domains, and two of them (CsYTH1 and CsYTH4) harbored extra CCCH and P/Q/N-rich domains. The phylogenesis, conserved motifs and exon-intron structure of YTH genes from cucumber, Arabidopsis and rice were also analyzed. The phylogenetically closely clustered YTHs shared similar gene structures and conserved motifs. An analysis of the cis-acting regulatory elements in the upstream region of these genes resulted in the identification of many cis-elements related to stress, hormone and development. Expression analysis based on the transcriptome data showed that some CsYTHs had development- or tissue-specific expression. In addition, their expression levels were altered under various stresses such as salt, drought, cold, and abscisic acid (ABA) treatments. These findings lay the foundation for the functional analysis of CsYTHs in the future.

MhYTP1 and MhYTP2 from Apple Confer Tolerance to Multiple Abiotic Stresses in Arabidopsis thaliana

The first YTH domain-containing RNA binding protein (YTP) was found in rat, where it was related to oxidative stress. Unlike characterizations in yeast and animals, functions of plant YTPs are less clear. Malus hupehensis (Pamp.) Rehd. YTP1 and YTP2 (MhYTP1 and MhYTP2) are known to be active in leaf senescence and fruit ripening. However, no research has been published about their roles in stress responses. Here, we investigate the stress-related functions of MhYTP1 and MhYTP2 in Arabidopsis thaliana. Both of the two genes participated in salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) signaling and play roles in plant responses to oxidative stress, chilling, high temperature, high salinity, and mannitol induced physiological drought stress. Moreover, MhYTP1 plays leading roles in SA and ABA signaling, and MhYTP2 plays leading roles in JA signaling and oxidative stress responses. These results will fill a gap in our knowledge about plant YTPs and stress responses and provide a foundation for future attempts to improve stress tolerance in apple.

Functions of two Malus hupehensis (Pamp.) Rehd. YTPs (MhYTP1 and MhYTP2) in biotic- and abiotic-stress responses

RNA binding proteins play important roles in plant responses to biotic and abiotic stresses. The YT521-B homology (YTH) domain-containing RNA binding protein (YTP) was first found in Rattus norvegicus and is related to oxygen-deficient stress. The Malus YTP gene family has 15 members. Results from their functional analysis will help researchers improve stress tolerance and fruit quality in apple. We cloned two homologous YTP family members in M. hupehensis - MhYTP1 and MhYTP2 - and identified their promoter regions that contain many cis-elements related to biotic and abiotic stresses. Both MhYTP1 and MhYTP2 can be induced by various treatments, e.g., methyl jasmonate (MeJA), salicylic acid (SA), abscisic acid (ABA), water-logging, water deficits, and high salinity. When compared with the wild type (WT), transgenic plants of 'GL-3' ('Royal Gala') apple that over-express MhYTP1 or MhYTP2 are more sensitive to D. mali infection, heat stress, and high salinity, more resistant to water-logging, chilling, drought and nutrition deficient conditions. All of these findings indicate that MhYTP1 and MhYTP2 participate in various biotic- and abiotic-stress responses.

Genome-wide identification and expression profiling of the SnRK2 gene family in Malus prunifolia

Sucrose non-fermenting-1-related protein kinase 2 (SnRK2) constitutes a small plant-specific serine/threonine kinase family with essential roles in the abscisic acid (ABA) signal pathway and in responses to osmotic stress. Although a genome-wide analysis of this family has been conducted in some species, little is known about SnRK2 genes in apple (Malus domestica). We identified 14 putative sequences encoding 12 deduced SnRK2 proteins within the apple genome. Gene chromosomal location and synteny analysis of the apple SnRK2 genes indicated that tandem and segmental duplications have likely contributed to the expansion and evolution of these genes. All 12 full-length coding sequences were confirmed by cloning from Malus prunifolia. The gene structure and motif compositions of the apple SnRK2 genes were analyzed. Phylogenetic analysis showed that MpSnRK2s could be classified into four groups. Profiling of these genes presented differential patterns of expression in various tissues. Under stress conditions, transcript levels for some family members were up-regulated in the leaves in response to drought, salinity, or ABA treatments. This suggested their possible roles in plant response to abiotic stress. Our findings provide essential information about SnRK2 genes in apple and will contribute to further functional dissection of this gene family.