The Alcohol Dehydrogenase Gene Family in Melon (Cucumis melo L.): Bioinformatic Analysis and Expression Patterns

Jin, Yazhong; Zhang, Chong; Liu, Wei; Tang, Yufan; Qi, Hongyan; Chen, Hao; Cao, Songxiao

doi:10.3389/fpls.2016.00670

ORIGINAL RESEARCH article

Front. Plant Sci., 18 May 2016

Sec. Plant Breeding

Volume 7 - 2016 | https://doi.org/10.3389/fpls.2016.00670

The Alcohol Dehydrogenase Gene Family in Melon (Cucumis melo L.): Bioinformatic Analysis and Expression Patterns

$\r\nYazhong Jin,*&#x;$ Yazhong Jin^1,2*†

Chong Zhang^1†

Wei Liu¹

Yufan Tang¹

Hongyan Qi^1*

Hao Chen¹

Songxiao Cao¹

¹Key Laboratory of Protected Horticulture of Education Ministry and Liaoning Province, Department of Horticulture, Shenyang Agricultural University, Shenyang, China
²College of Agriculture, Heilongjiang Bayi Agricultural University, Daqing, China

Alcohol dehydrogenases (ADH), encoded by multigene family in plants, play a critical role in plant growth, development, adaptation, fruit ripening and aroma production. Thirteen ADH genes were identified in melon genome, including 12 ADHs and one formaldehyde dehydrogenease (FDH), designated CmADH1-12 and CmFDH1, in which CmADH1 and CmADH2 have been isolated in Cantaloupe. ADH genes shared a lower identity with each other at the protein level and had different intron-exon structure at nucleotide level. No typical signal peptides were found in all CmADHs, and CmADH proteins might locate in the cytoplasm. The phylogenetic tree revealed that 13 ADH genes were divided into three groups respectively, namely long-, medium-, and short-chain ADH subfamily, and CmADH1,3-11, which belongs to the medium-chain ADH subfamily, fell into six medium-chain ADH subgroups. CmADH12 may belong to the long-chain ADH subfamily, while CmFDH1 may be a Class III ADH and serve as an ancestral ADH in melon. Expression profiling revealed that CmADH1, CmADH2, CmADH10 and CmFDH1 were moderately or strongly expressed in different vegetative tissues and fruit at medium and late developmental stages, while CmADH8 and CmADH12 were highly expressed in fruit after 20 days. CmADH3 showed preferential expression in young tissues. CmADH4 only had slight expression in root. Promoter analysis revealed several motifs of CmADH genes involved in the gene expression modulated by various hormones, and the response pattern of CmADH genes to ABA, IAA and ethylene were different. These CmADHs were divided into ethylene-sensitive and –insensitive groups, and the functions of CmADHs were discussed.

Introduction

Alcohol dehydrogenases (ADH, EC 1.1.1.1) belong to dehydrogenase enzymes superfamily, and are widely distributed in all types of organisms (Chase, 1999; Jönvall et al., 2010; Strommer, 2011; Alka et al., 2013). ADH enzyme has long been the subject of molecular studies and encoded by multigene family in eukaryotic and prokaryotic kingdoms, catalyzing the interconversion between alcohols and aldehydes (Höög et al., 2003; Thompson et al., 2007). In humans, animals, yeast and bacterias, ADHs have been widely investigated (Khan et al., 2010; Kumar et al., 2012; Alka et al., 2013; Çelik and Aktas, 2013; Jönvall et al., 2013; Plapp et al., 2013; Quaglia et al., 2013), and ADH genes were involved in an astonishingly wide range of metabolic processes (Höög et al., 2003; Thompson et al., 2007; Strommer, 2011). Hence, these ADH genes were classed into several main superfamilies respectively, namly medium- (∼350 amino acid residues), short- (∼250 residues) and long-chain ADH or Iron-ADH genes superfamilies (600–750 residues or approximately 385 amino acid residues up to almost 900 residues) (Chase, 1999; Deng et al., 2002; Alka et al., 2013; Jönvall et al., 2013), and the medium-chain ADHs clustered into eight classes in vertebrates, based on sequence similarity, catalytic features and gene expression patterns (Thompson et al., 2007). But functional roles of these ADHs are not fully understood.

Sets of ADH genes or ADH-like genes have been identified in the genomes of poaceae, rosaceae, brassicaceae, fabaceae, and pinaceae plants (Thompson et al., 2010; Zheng et al., 2011). So far, most members of the ADHs in plants characterized at the gene level belonged to the medium-chain ADH protein superfamily (including ADH1, EC 1.1.1.1 and FDH, Class III ADH, EC 1.2.1.1; Chase, 1999), which usually contains zinc ligands in their active site (Tesnière and Verriès, 2000; Garabagi et al., 2005; ManrÍquez et al., 2006; Koutsompogeras et al., 2010; Singh et al., 2010; Komatsu et al., 2011; Pathuri et al., 2011; Bukh et al., 2012; Iaria et al., 2012; Min et al., 2012; Cheng F.F. et al., 2013), and the short-chain ADH protein superfamily, which lacks zinc-liganding cysteine residues in their coenzyme binding regions, and the molecular functions of only a few short-chain ADH genes were known (ManrÍquez et al., 2006; Gonzalez-Agüero et al., 2009; Kim et al., 2009; Strommer, 2011; Moummou et al., 2012). Other ADHs have been not noted in plants, such as long-chain ADHs or Iron-ADH genes. On the basis of phylogenetic analyses, the medium-chain ADH genes in plant were distributed in different subgroups (Perry and Furnier, 1996; Gaut et al., 1999; Small and Wendel, 2000; Komatsu et al., 2011; Zheng et al., 2011). Previous studies highlighted that the medium-chain ADHs in plant were involved in response to abiotic and biotic stress which induced the specific expression of these ADHs in different tissues of soybean, wheat and barley, implying that they may participate in different tissues development in stresses (Komatsu et al., 2011; Pathuri et al., 2011; Yamauchi et al., 2014). Moreover, regardless of medium- chain ADHs or short-chain ADHs, ADH genes have been shown to play a major role in fruit ripening and aroma synthesis (Tesnière and Verriès, 2000; ManrÍquez et al., 2006; Singh et al., 2010; Zhang et al., 2011; Zheng et al., 2011; Bukh et al., 2012; Iaria et al., 2012; Moummou et al., 2012). However, the actual participation of medium-chain ADHs in aroma volatile production in vivo has been only clearly demonstrated in tomato fruit by over-expressing or down-regulating the LeADH2 gene (Speirs et al., 1998). Other detailed studies related to role of ADH genes (including medium- and short-chain ADHs) in relation to aroma synthesis have been proved in tomato, melon and mango by recombinant protein (ManrÍquez et al., 2006; Singh et al., 2010; Moummou et al., 2012). In addition, persimmon DkADH, Artemisia annua AaADH2, Sedum sarmentosum SsADH (belonging to medium-chain ADH family) and other short-chain ADHs have been implicated in other primary and secondary metabolisms (Polichuk et al., 2010; Tonfack et al., 2011; Sung et al., 2012). In panax ginseng, the expression pattern of PgADH under abiotic stimuli suggested that a short-chain ADH, PgADH, was involved in responses to hormone-related environmental stresses (Kim et al., 2009). Meanwhile, the tissue-specific patterns of these ADH gene expressions among different tissue under various environmental stresses and during fruit ripening and aroma synthesis also provided a few evidences of functional specialization. We conclude further investigation is required to make functional annotation for the majority of predicted ADH in plant genomes.

In melon, two ADH genes, CmADH1 and CmADH2, have been identified and characterized (ManrÍquez et al., 2006). The release of melon genome provides an easy way to identify new members of ADH gene family. This study is a continuation of CAD genes (belonging to dehydrogenase enzymes superfamily) research in melon (Jin et al., 2014). In this study, we identified 12 CmADHs, including CmADH1 and CmADH2, and one FDH gene from the melon genome, and compared ADH sequences from a wide variety of plants, making full use of the available plant genome sequences (Arabidopsis, mango, tomato, grape et al.). Alignment and phylogenetic analysis of the ADH gene family with related ADH proteins from other species indicated that these ADH genes were classified into three groups and the medium-chain melon ADHs were clustered into six subgroups separately. We analyzed the structure and the promoter of melon ADH genes, and also investigated these ADH genes transcripts in response to various fruit development stages and various plant hormones, and monitored tissue-specific expression. We reported here the results of these analyses, suggesting that the CmADHs may be involved in melon fruit development induced by various hormones and during development and ripening. Information gained from this investigation will significantly advance the understanding of the function of predicted ADH genes in Melon.

Materials and Methods

Plant Materials

All the experiments were conducted on a commercial oriental melons (Cucumis melo var. makuwa Makino) cultivar, ‘CaiHong7.’ All melon plants were grown in a greenhouse under standard cultural practices at Shenyang Agricultural University, Shenyang, China, from March to June in 2012, and the melon fruit were harvested at different development stages (1, 5, 10, 15, 20, 25, 30, 33, 36, 39, 42, 45, and 48 days after anthesis) for expression analysis as previously described (Jin et al., 2014). Physiological maturity of this melon is about 36 days after anthesis. To allow analysis of tissue-specific gene expression, mature leaf, developing leaf, pistillate and staminate flower, young stems and root tissues were collected from melon plants grown in a greenhouse, then frozen with liquid nitrogen and kept in -80°C refrigerator for further use.

Hormone Treatments

Mature unripe oriental melons were harvested at 30 days after anthesis (pre-climacteric stage) with the same node of the plant at a mature green stage, before the onset of ripening. No blemished or diseased Fruits were chosen for treatment. Exogenous ethylene (100 μL/L), 1-methylcyclopropene (1-MCP, an ethylene perception inhibitor) (100 μL/L) and control (water) treatments were performed according to our previous report (Jin et al., 2014). Then melon fruit of three treatments were allowed to ripen for 12 days at 23°C in air only. Three replications were carried out for treatments, and each replicate consisted of 15 fruits unless indicated otherwise. Fresh tissue was sampled every 48 h, frozen in liquid nitrogen and stored at -80°C until further use. Furthermore, for abscisic acid (ABA;100 μM) and auxin (IAA;100 μM) treatment, the fruit disks were dipped in a solution containing ABA or IAA in 0.2% teepol (detergent) and vacuum infiltrated for 2 h as described previously (Jin et al., 2014). Infiltrated melon fruits with 0.2% teepol were used as control. After treatment, fruit flesh were frozen in liquid nitrogen and stored at -80°C. Three replications were carried out for treatments or control groups.

Identification of Melon ADH Genes

The keyword ‘alcohol dehydrogenase’ was used to search for the melon ADH sequences from the melon (Cucumis melon L.) genome¹ (Garcia-Mas et al., 2012). In addition, to confirm the accuracy of these genes, the predicted ADH-like gene sequences were compared to ADH proteins in other species by a BLASTp retrieve. Only those sequences with a high score (>200) were selected.

Sequence Analysis

The sequence analysis of CmADHs, including ADH protein prediction, functional domains and promoter analysis, signal peptides and disulfide bond prediction and CmADHs subcellular localization prediction, were performed according to the method described by Jin et al. (2014). Multiple alignments were performed with other known plant ADH proteins using Clustal Omega program and GENEDOC. Based on the neighbor-joining method (minimum evolution criterion, bootstrap values performed on 1000 replicates), phylogenetic analyses of putative melon ADH proteins were performed using MEGA5² program.

RNA Isolation and cDNA Synthesis

Total RNA from fruit samples, leaves, stem, root, and flower material, were extracted by the method described by Jin et al. (2014), quantified by the NanoDrop spectrophotometer ND-1000 and checked for integrity by electrophoresis (28S rRNA/18S rRNA ratios). RNA was treated by DNAse I (Promega, Madison, WI, USA) at 37°C for 50 min, re-precipitated and concentrated (40 μL) to remove any trace of genomic DNA. cDNA synthesis was initiated from DNase I-treated RNA using M-MLV RTase cDNA Synthesis Kit following the manufacturer’s instructions (Cat#D6130, TaKaRa, Tokyo, Japan).

Semi-quantitative PCR and Real-Time PCR

For semi-quantitative PCR and real-time PCR, gene-specific oligonucleotide primers for each ADH gene were designed by Primer3³ and were listed in Table 1. The specificity of each pair of primers was determined by agarose gel electrophoresis and PCR products resequencing, then semi-quantitative PCR and Real time PCR were performed for gene expression studies according to the method described by Jin et al. (2014). 18SrRNA DNA fragment (148 bp) of melon as an internal control.

TABLE 1

TABLE 1. Semi-quantitative PCR and Real-time PCR primers.

Statistical Analysis

Data are expressed as mean values ± standard deviation of three independent experiments (n = 3). The data were analyzed by the analysis of variance (ANOVA) using the SPSS 13.0 statistics program, and statistical significance of differences was were calculated by a one-way ANOVA following Duncan’s multiple range tests for each experiment at a P < 0.05 level. Origin 8.0 (OriginLab, Northampton, MA, USA) was used to draw the figures.

Results

Identification of Melon ADH Genes

All the predicted ADH genes in the melon genome were collected and compared with ADH genes in other species. Hence, the presence of functional domains was checked via NCBI’s Conserved Domain Database (CDD)⁴, and 13 genes (including 12ADHs and 1 FDH), encoding full-or nearly full-length functional proteins, were identified in the melon genome.

By ExPASy tools, we found that the longest ADHs protein consisted of 635 amino acid residues, and the shortest consisted of 266 amino acid residues. The ORF length ranged from 801 to 1908 nucleotides. The predicted molecular weight and isoelectric points of all ADHs proteins ranged from 29.0 kDa/5.19–68.5 kDa/9.06, respectively (Table 2). Of these ADH genes, two of them have been reported previously and were labeled with their published names (ManrÍquez et al., 2006) (Table 2). For the other 11 unreported genes, we gave them names by adding a number to their family name in the order which they were searched, in which one ADH was named as CmFDH1 on account of the presence of same functional domains with reported FDH genes in other species (Dolferus et al., 1997).

TABLE 2

TABLE 2. The information of ADH genes in melon.

Intron-Exon Structure of Melon ADH Genes

Structural analysis of the identified melon ADH genes revealed different intron-exon patterns both in relation to position and number of introns which ranged from 2 to 17 per gene (Figure 1). Furthermore, the substantial differences in the size between the exons were observed. The standard number of introns in plant ADH genes in the common ancestor is 9, located at equivalent positions in ADH genes throughout the plant kingdom (Strommer, 2011). However, plant ethanol-active ADH genes have evolved from class III ADH genes by gene duplication and acquisition of new substrate specificities (Dolferus et al., 1997). Based on Figure 1, CmFDH1 (Class III ADH;GSH-FDH) contained nine exons and eight introns, while the intron-exon structures of other 12 CmADHs were relatively complex. CmADH1, CmADH8, CmADH10 and CmADH11 had seven introns, and CmADH3, CmADH5, CmADH6 and CmADH7 had nine introns, and CmADH4 and CmADH9 contain eight introns; CmADH2 had two introns, moreover CmADH12 contained 17 introns (Figure 1). For ADH-plant genes, introns patterns contributed more evidences for gene acquisition through duplications, insertion or loss of introns (Charlesworth et al., 1998; Small and Wendel, 2000; Strommer, 2011). Additionally, both 76-bp and 83-bp sequence were found in six sequences, and a 76-bp sequence was found in these sequences, except CmADH2,8,12 (Figure 1). Nevertheless, it was unclear whether these sequence encoded the putative functional region of ADH genes.

FIGURE 1

FIGURE 1. Intron-exon structures of ADH genes from melon. Exons and introns were indicated by open boxes and lines respectively. Numbers above boxes indicated the exon sizes. The intron sizes were not to scale. The names of ADH genes and intron-exon structure were indicated at the left and right sides respectively.

Sequence Characterization of CmADH Genes

To understand the possible functional divergence of the individual members of melon ADH genes superfamily, the sequences of a set of ADH family members from other plant species, including cucumber, Arabidopsis, grape, potato, tomato, mango, apricot, peach, soybean, cocoa, were analyzed using a Neighbor-joining (NJ) phylogenetic tree by MAGA5. As shown in Figure 2, the phylogenetic tree divided these ADHs into three clades, namely medium-, short-, and long-chain ADHs, and the medium-chain ADHs in melon were clustered into six medium-chain ADHs subgroups, called as midum-chain groups I–VI, based on their amino acids sequence (Supplementary Figure S1). In Cantaloupe, CmADH1 and CmADH2, two highly divergent genes, have been isolated, and it has been characterized that CmADH1 belonged to the medium-chain zinc-binding type of ADHs and CmADH2 was positioned in the short-chain ADHs (Singh et al., 2010). Moreover, CmADH3-11 belonged to medium-chain zinc-binding type of ADHs and have not been characterized, in which CmADH1, Pa-ADH1 (apricot) and Le-ADHs (tomato), having been characterized, were involved in fruit-ripening and aroma volatiles biosynthesis (ManrÍquez et al., 2006; Gonzalez-Agüero et al., 2009; Moummou et al., 2012). CmADH12 belonged to long-chain ADHs which have not been characterized in plant; CmFDH1 clustered into medium-chain zinc-binding type of FDHs, known as Class III ADH, and AtFDH1 (Arabidopsi) and OsFDH1 (rice) have been cloned (Dolferus et al., 1997); In addition, it is now recognized that FDHs are involved in plant development and stress response (Feechan et al., 2005; Rustérucci et al., 2007; Lee et al., 2008; Espunya et al., 2012).

FIGURE 2

FIGURE 2. Phylogenetic relationship among CmADHs and ADHs from other plant species. The amino acid sequences were aligned by the Clustal Omega program, and the neighborjoining tree was drawn with MEGA5. The corresponding GenBank and the melon genome (https://melonomics.net/) were noted in the phylogenetic tree, and the accession number in the melon genome were CmADH1 (MELO3C023685P4), CmADH2 (MELO3C014897P1), CmADH3 (MELO3C026552P1), CmADH4 (MELO3 C027151P1), CmADH5 (MELO3C005792P1), CmADH6 (MELO3C026553P1), CmADH7 (MELO3C002189P1), CmADH8 (MELO3C003251P1), CmADH9 (MELO3C0110 43P1), CmADH10 (MELO3C026554P2), CmADH11 (MELO3C023687P1), CmADH12 (MELO3C019503P1), and CmFDH1 (MELO3C022399P1). The number for each interior branch was the percentage of bootstraps value (1000 replicates). Black circle denoted 13 CmADHs. ADHs from other plants in our paper were in supporting information (Supplementary Table S3).

Multiple alignments revealed that the CmADH genes showed lower identity in protein level with each other, as noted previously in Cantaloupe (ManrÍquez et al., 2006). Among the medium-chain zinc-binding type of CmADHs, there was lower identity between CmADH1 and other medium-chain zinc-binding CmADHs in the range of 30–88.51% (Figure 3), whereas CmADH1 showed the highest sequence identities to medium-chain ADH proteins from cucumber, mango, grape, tomato and Arabidopsis in range of 80–96.08% (Supplementary Figure S2), and other medium-chain CmADHs also showed higher sequence identities with other plants (Supplementary Figures S3–S9). All of the medium-chain zinc-binding CmADHs had the highly conserved Zn1 -binding signature [GHE(X)2G(X)5G(X)2V], the Zn2 structural motif [GD(X)9,10C(X)2C(X)2C(X)7C], and the NADPH-binding domain [GXG(X)2G] motif (so-called Rossmann fold) (Figure 3; Supplementary Figures S2–S9), which was the same with that of CmADH1 proved (Singh et al., 2010), suggesting that these proteins appear to be zinc-dependent ADHs and members of the plant ADH protein family and the medium-chain dehydrogenase/reductase (MDR) superfamily (Persson et al., 1991; McKie et al., 1993). Nevertheless, CmADH11, a medium-chain melon ADH, lacked the key determinant residues of NADPH-binding domain at protein level (Supplementary Figure S9).

FIGURE 3

FIGURE 3. Alignment of amino acid sequences of CmADHs. Conserved important regions identified previously were marked as follows: white arrows denoted catalytic zinc ion coordinating residue, black arrows denoted structural Zn ion coordinating residue, the black circle denoted key residues for substrate specificity. The black square denoted key Phe299/Gly(300) residues for substrate specificity. The white square denoted key Trp199 and Asp123 residues for substrate binding. Locations of the Zn1, Zn2, and NADPH binding domains were shown in boxes. The alignment was performed with the Clustal Omega program.

Furthermore, CmADH2, a short-chain ADH protein, and CmADH12, a long-chain ADH protein, had highly homology with corresponding ADHs from other plants (Supplementary Figures S10 and S11). The sequence alignment showed that the CmADH12 protein sequence was found most similar to corresponding ADH-like proteins from other plants, with identities close to 95.43% (Supplementary Figure S11), and contained an N-terminal cofactor binding pattern of typically TGXXXGXG and an active-site pattern of YXXXK which are two conserved regions around amino acids of short-chain dehydrogenases in plants (Joërnvall et al., 1999). Additionally, we used the protein database in NCBI as a source of amino acids sequence bioinformatics for CmADH12, and found that CmADH12 possessed a coenzyme (NADPH) binding domain AAAGGTG, suggesting that CmADH12 may belong to long-chain ADH subfamily, but may have the same catalytic function as shot-chain dehydrogenases. Interestingly, the analysis of the amino acids sequence of CmADH12 showed that it did not contain structurally and functionally important iron-binding motif (G-X-X-H-X-X-A-H-X-X-GX-X-X-X-X-P-H-G), and NAD(H)-or NADP(H)-binding fingerprint patterns G-X-G-X-X-G or G-X-G-X-X-A, GGGSXXD and A/GXXATGG, and Trp docking motif (A-x-[DN]-x-x-T-G-[DEK]-x-x-W), possessed in long-chain or iron-dependent ADH (group III ADH) protein of prokaryotes, human and other mammals (Inoue et al., 1989; Deng et al., 2002; Lyon et al., 2009; Elleuche et al., 2013, 2014). Additionally, Filling et al. (2002) described several conserved residues of short-chain dehydrogenase, involved in coenzyme binding and catalysis, such as Thr12, Asp60, Asn86, Asn87, Ala88, Gly89, Gly131, Gly183, Asn111, Ser151 and Tyz157. In melon, CmADH2 only contains five residues of these conserved residues, namely Thr12, Asp60, Gly183, Ser151 and Tyz157 (Supplementary Figure S10). These findings suggested that the biological functions of plant short-chain ADH proteins might be related to the diversity in amino acid sequence.

No typical signal peptides were found in all CmADHs after analyzing their N-terminals using the SignalP software. CmADHs subcellular localization prediction showed that these ADH genes might exist in the cytoplasm (Supplementary Table S2). Moreover, the transmembrane topology predictions of 13 CmADHs using TMHMM 2.0 software and ABTMpro showed that there was no internal transmembrane segment in CmADHs. We also found that there were 4 and 3 disulfide bonds in 12 melon ADH proteins by Scratch Protein Predictor, respectively, but there were 6 disulfide bonds in CmADH12 protein.

Promoter Sequence Analysis

Analysis of promoter sequences for the melon ADH genes allowed us to identify several motifs that were known to be involved in the regulation of gene expression in various developmental and physiological processes. A promoter motif search showed that CmADHs promoter contained putative regulatory elements corresponding to known cis-elements of eukaryotic genes (Barakat et al., 2009; Kim et al., 2010; Cheng H. et al., 2013). In melon ADH genes promoter, there were mainly two kinds of motifs, namely, cis-acting element involved in defense and stress responsiveness (such as hypoxia stress, heat stress) and cis-acting regulatory element involved in the response to various hormones, such as auxin (IAA), ethylene (ETH), ABA, salicylic acid (SA), and methyl jasmonate (MeJA) (Supplementary Table S1) (Barakat et al., 2009). Nevertheless, the most abundant motifs detected were the cis-elements involved in the response to biotic and abiotic stresses, including TC-rich repeats, LTR, W-BOX, ARE, HSE and MBS (Jin et al., 2014). At the same time, we still found that there were putative cis-acting regulatory elements involved in the MeJA-responsiveness, namely CGTA-motif, TGACG-motif and CCAAT-motif (Matton et al., 1990; Prashant et al., 2012; Cheng H. et al., 2013). TCA-element clustered in the promoters of CmFDH1, CmADH1, CmADH3, CmADH5, CmADH7-9 and CmADH12, and SARE were also present in CmADH3 promoter, both of which were cis-acting regulatory element involved in the SA-responsiveness (Prashant et al., 2012; Cheng H. et al., 2013). ERE, ethylene-responsive element, was found in the promoter of CmFDH1, CmADH1, CmADH2, CmADH3, CmADH5, CmADH6, CmADH8 and CmADH12. TGA (Auxin responsiveness), ABRE (ABA responsiveness), p-box and GARE (gibberellin responsiveness) were also present in several CmADHs promoter. Additionally, promoter sequence analysis revealed several conserved motifs, such as MBS-I and MBS-II (MYB binding site involved in flavonoid biosynthesis and drought-inducibility), HSE (heat stress responsiveness) and LTR (Low temperature responsiveness) (Prashant et al., 2012). CmFDH1, CmADH3 and CmADH9 possessed elicitor responsive element and enhancer. These results indicated that transcriptional regulation of these CmADH genes might be involved in fruit development and ripening and senescence and in the response to various stresses.

CmADHs Expression in Vegetative Tissues

In order to investigate the transcript levels of these 13 CmADH genes in different organs and tissues in melon plant, we collected samples of root, developing leaf, mature leaf, young stem, pistillate and staminate flower petal. Expression analysis using semi-quantitative PCR and real-time PCR showed that 12 CmADH and one CmFDH were constitutively expressed in these parts, but greatly varied in different tissues, and there were differences in same tissues between CmADHs (Figure 4; Supplementary Figure S12). Of the 12 ADH genes identified in melons, CmADH2 and CmADH10 were expressed most abundantly in root, developing leaf, mature leaf, young stem, pistillate and staminate flower. Because of the higher expression of CmADH2 and CmADH10 seleted as “1,” real-time PCR analysis of CmADH2 and CmADH10 showed a lower relative expression level (Figure 4). The expression of CmFDH1 was higher in root, young stem and pistillate flower. CmADH9 and CmADH11 were either not expressed or expressed at very low levels in these tissues. Of all analyzed tissues, the CmADH4 gene was highly expressed in the root, and only rarely expressed in other vegetative organs. CmADH1 expression was the greatest in root and young stem, and CmADH3 was expressed in vegetative organs, with the highest expression in the young leaf and stem. Furthermore, CmADH5, CmADH6 and CmADH7 showed similar expression patterns between vegetative organs, with lower transcript levels, and without significant differences except root (P < 0.05). CmADH8 and CmADH12 were also expressed in all tissues except mature leaf, showing lower expression in young leaf.

FIGURE 4

FIGURE 4. Transcript levels of CmADHs in different melon organs. Gene expressions of CmADHs in different organs in melon plants were determined by real-time PCR in root, developing leaves, mature leaves, young stems, pistillate flower petals and staminate flower petals in melon plants. 18s was used as internal control. The expression level of the genes in mature leaves was set as “1.0.” Data represented the means ± SD (n = 3) of three biological samples. The experiments were carried out in triplicate.

CmADHs Expressions during Melon Fruit Development

Transcript analysis indicated that these 12 CmADH genes and one CmFDH gene studied here were specifically expressed in fruit of different development stages, with remarkable differences in fruit (Figure 5; Supplementary Figure S12). The pattern of changes in transcript levels was similar for CmADH1, CmADH3, CmADH8, CmADH12 and CmFDH1 from 1 to 15 days after anthesis with gradual or sharp decrease at 15 days. While the expression of CmADH3, CmADH8 and CmADH12 showed an increase after 15 days, and the levels of these 3 ADH genes transcript were always maintained up to 48 days after anthesis. However, CmADH1 had higher expression in fruit during development, except from 30 to 32 days and 48 days after anthesis. In contrast, CmADH2 and CmADH10 were consistently expressed with an increase in transcript abundance, which subsequently remained at a relatively constant level through to harvest. Additionally, CmADH4, CmADH5, CmADH7, CmADH9 and CmADH11 were either not expressed or expressed at very low levels, and the transcript levels of CmADH6 were always weakly maintained during fruit development.

FIGURE 5

FIGURE 5. CmADHs relative expression in developing stages of melon fruit after anthesis were determined by real-time PCR. 18s was used as internal control. The expression level of CmADHs in melon fruit at 15 days after pollination was set as “1.0.” Data represented the means ± SD (n = 3) of three biological samples. The experiments were carried out in triplicate.

Effects of Hormones on CmADH Genes Expressions

The real-time PCR and semi-quantitative PCR results revealed that CmADHs were induced by IAA and ABA (though expression was lower in case of ABA), and showed different responses to IAA and ABA, except CmADH11 (Figure 6; Supplementary Figure S13). The CmADH1-4, CmADH10 and CmFDH1 were strongly induced by ABA, and the relative expression varied between these genes. Moreover, CmADH2 showed highest transcript levels. In contrast, IAA treatment strongly induced CmADH6, CmADH9 and CmADH12 expression, while other CmADHs were seemed to be insensitive to treatment with IAA and continued to be weakly expressed at a basal level.

FIGURE 6

FIGURE 6. The expression of CmADHs in melon fruit after treatments with IAA and ABA. IAA and ABA (100 μM) treatments were given for 2 h as described in Section “Materials and Methods.” Expression analysis was carried out by real-time PCR. For each gene, the relative abundance of mRNA was normalized against the 18S in the corresponding samples. The expression level of the genes in untreated melon fruit by IAA and ABA was set as “1.0.” Data represented the means ± SD (n = 3) of three biological samples. The experiments were carried out in triplicate.

Since oriental melons is a climacteric fruit which requires ethylene for initiation of ripening, we checked whether these CmADHs were regulated by ethylene. Transcript accumulation of different CmADHs in melon were investigated during the course of ethylene induced ripening of harvested melon for 12 days at 23°C. A basal level of expression of CmADAHs transcript was detected in all of 12 ADH genes and one FDH gene in ethylene untreated fruit on day 1. The transcript levels of CmADH1-3, CmADH8, CmADH12 and CmFDH1 shot up on day 1 in ethylene treated fruit, apart from CmADH4-7, CmADH9 and CmADH11 (Figure 7; Supplementary Figure S14). However, levels of CmADH1, CmADH2 and CmADH3 transcript gradually decreased from day 1 levels in subsequent days, but levels of CmADH8, CmADH12 and CmFDH1 gradually increased from day 1 levels, and reduced from day 6 thereafter. Of the 12 ADH genes, CmADH4-7, CmADH9 and CmADH11 transcripts showed the lowest steady state levels in all of treatments during the course of melon ripening. Since these 6 CmADHs clustered into the medium-chain ADH subfamily, we checked the expression of CmADH4-7, CmADH9 and CmADH11 transcripts by semi-quantitative PCR at high annealing temperature using gene-specific primers to ensure that the transcript patterns were not due to cross hybridization. We found that these 6 CmADHs transcript levels were still lower and did not give clear results when 20 cycles of PCR were carried out as in case of other ADH genes. So for all semi-quantitative PCR experiments related to CmADH4-7, CmADH9 and CmADH11, 35 cycles of PCR were carried out. Our results of semi-quantitative PCRs matched with real time qPCR analysis, and confirmed that all 12 ADHs and one FDH were expressed at different levels during ripening. The expressions of CmADHs were also studied in 1-MCP + ethylene treated fruits. All 12 ADHs and one FDH showed transcript accumulation in 1-MCP + ethylene treated fruits, but the levels were far lower compared to the levels in ethylene treated fruits except CmADH4-7, CmADH9, CmADH10 and CmADH11, especially on day 1 post treatment (Figure 7; Supplementary Figure S14). Furthermore, CmADH10 always showed stable and strong expression in all of treatments during the course of melon ripening. Therefore, these finding implied that CmADH4-7 and CmADH9-11 appeared to be not induced by ethylene. Since the patterns were same in real-time qPCR and semi-quantitative PCR experiments, further studies will be required.

FIGURE 7

FIGURE 7. The expression levels of CmADHs in melon fruit after treatments with ethylene and 1-MCP. The transcript levels of CmADHs were measured by real-time PCR in melon fruit treated, and 18S was used as internal control. The expression of the genes in untreated melon fruit after 1 day of storage was set to 1.0. Data represented the means ± SD (n = 3) of three biological samples. The experiments were carried out in triplicate.

Discussion

Alcohol dehydrogenases, belonging to a group of dehydrogenase/reductases superfamily, are broadly distributed in all types of organisms and represent one of the most abundant classes of enzymes. More importantly, the ADH genes are involved in stress responses, elicitors and ABA (Matton et al., 1990; Kim et al., 2009; Strommer, 2011), and also closely associated with seedling, pollen development, root tissues, fruit ripening and aroma synthesis (Dolferus et al., 1994; Chung and Ferl, 1999; ManrÍquez et al., 2006; Singh et al., 2010; Kubienova et al., 2014). The ADH gene family has been investigated in a number of plant species, such as grape (Tesniere et al., 2004), Mango (Singh et al., 2010), tomato (Speirs et al., 1998), peaches (Zhang et al., 2010), apple (Defilippi et al., 2005a), soybean (Komatsu et al., 2011) and so on. The aim of the present study was to identify and analyze genes coding for the ADH enzyme in melon.

Bioinformatic Analysis

In the melon genome database, we identified 13 ADH/ADH-like genes (including 12 ADHs and 1 FDH) coding for putative ADH enzymes, in which CmADH1 and CmADH2 have been cloned and characterized (ManrÍquez et al., 2006). Our research results also demonstrated that these melon ADH genes had very low sequence identity at protein level (ranging from 4.8 to 52.48%) and different intron-exon patterns at nucleotide level (Figure 3). Up to the present, research on ADH in plants have been focused on the medium- and short- chain ADH superfamilies (ManrÍquez et al., 2006; Kim et al., 2009; Singh et al., 2010; Iaria et al., 2012; Yamauchi et al., 2014; Chen et al., 2015). Furthermore, plant FDH, with hydroxylmethyl-glutathione or formaldehyde as preferred substrates, belonged to the medium-chain ADHs (Strommer, 2011), while it is weakly active or inactive on n-octanol or ethanol (Chase, 1999); Hence, in regard to substrate specifity, FDH was different from ADH1 in plant. In our study, CmFDH1 clustered in medium-chain ADH superfamily, but the amino acid sequence of CmFDH1 shared 28.33% identity with CmADH1. Based on the convergent evolution and structure of ADH gene-families in plant and animal, the plant ADH family was deemed to originate from a glutathione-dependent formaldehyde dehydrogenase gene (GSH-FDH) (Strommer, 2011). In maize, rice, pea and Arabidopsis, FDH genes were cloned and purified, but their functions in plants were not yet clear (Chase, 1999). As of present, there is still insufficient information on the function of FDHs in higher plants (Espunya et al., 2006; Strommer, 2011; Nian et al., 2013), hence it would be widely concerned by researchers in future.

The phylogenetic tree showed that ADH genes in melon clustered into 3 groups (Figure 2), indicating that ADH genes in melon exist in multiple isoforms. However, the phylogenetic tree obtained differed from previous analyses which grouped the plant ADH genes into two main classes (Chase, 1999; Duester et al., 1999; Thompson et al., 2007; Borràs et al., 2014), without the long-chain or Iron-activated ADH reported in plants. At the same time, we found that CmADHs sequences clustered close to that of other various species, suggesting that there are species- or lineage-specific ADH gene duplications (Komatsu et al., 2011) or a number of separate duplication events within melon which may be not homologous (Thompson et al., 2007). Additionally, these medium-chain CmADHs, divided into six subgroups (Supplementary Figure S1 and S4–S8), also suggested that it seems possible to be different substrate specificities among melon medium-chain ADH genes and to be the same with that of vertebrates in which the medium-chain ADHs clustered in different subgroups (Duester et al., 1999). Hitherto, melon CmADH1 (ManrÍquez et al., 2006), tomato Le-ADH2 and LeADH3 (Speirs et al., 1998), mango MiADH1,2 (Singh et al., 2010), grape Vv-ADH2 (Tesniere et al., 2004) and olive OeADH (Iaria et al., 2012), clustering a same class, belonged to the medium-chain zinc-binding type of ADHs (Supplementary Figure S1), and in vitro, the functional verification results revealed that this type of plant ADHs were implicated in plant development, fruit ripening and aroma production, and response to stresses.

Interestingly, contrary to medium-chain ADHs, CmADH11 only retained a conserved functional protein domain, Zn2-binding region [Cys(102), Cys(105), Cys(108) and Cys(116)] (Supplementary Figure S9). Although CmADH1, CmADH4 and CmADH11 were grouped into a class in the phylogenetic tree (Figure 2), CmADH11 showed a lower degree of homology to CmADH1 and CmADH4, implying that CmADH11 may be a member of zinc-medium-chain ADHs subfamily, and that these 3 CmADHs probably have similar functions, but more researches are required to examine this speculation. In addition, short-chain ADHs were found to be present in all living kingdoms. In tomato, a short-chain ADH, SlscADH1 protein contained motifs characteristic of classical short-chain dehydrogenases (Persson et al., 1991; Moummou et al., 2012). Nevertheless, in melon and ginseng, the protein sequence of short-chain ADH only encompassed some of these conserved regions, which were demonstrated to be related to substrate binding and active sites (ManrÍquez et al., 2006; Kim et al., 2009; Moummou et al., 2012). Furthermore, we found that the deduced protein of CmADH12, encoded by 635 amino acids, appeared to belong to long-chain ADH subfamily. In microorganisms, the sequences of the so-called long-chain or iron-dependent ADH genes were encoded by 600–750 amino acids, or up to over 900 amino acids (Williamson and Paquin, 1987; Goodlove et al., 1989; Chase, 1999). Whereas, CmADH12 did not has the typical protein domains of long-chain or iron-dependent ADH proteins, possessed in corresponding ADHs of prokaryotes, human and other mammals (Inoue et al., 1989; Deng et al., 2002; Lyon et al., 2009; Elleuche et al., 2013, 2014). Interestingly, CmADH12 protein sequence exhibited two motifs typical of classical SDR-ADHs, namely an N-terminal coenzyme-binding pattern of typically GXXXGXG and an active-site pattern of YXXXK (Supplementary Figure S11); Simultaneously, this sequence possessed a NADPH-binding domain AXXGGTG that appears to be motifs typical of classical medium-chain ADHs by NCBI search, implying that CmADH12 may be consisted of a short-chain dehydrogenase protein and a medium-chain dehydrogenase. Up to now, for the long-chain ADH genes or iron-activated ADHs, no member of this type ADHs has been described in plant. Here, we reported a novel melon ADH, and the majority of predicted CmADH12 in melon genome still awaits functional annotation. With the report of melon, CmADH1 and CmADH2, belonging to the medium-chain ADHs and short-chain ADHs respectively (ManrÍquez et al., 2006), the identification of a long-chain ADH, CmADH12, will make a bridge to connect the gap in the distribution of medium- and short-chain ADHs, and will represent significant improvements on the understanding of the evolution of ADH in higher plants.

Expression Analysis in Different Tissue

In plants, multiple homologous ADH genes, thought to have distinct roles, have been implicated in the response to biotic and abiotic stresses, particularly, in the tissues development under stresses (Abe et al., 2003; Delessert et al., 2004; Kumar et al., 2010; Frazier et al., 2011; Komatsu et al., 2011; Qi et al., 2012). At the same time, many previous studies indicated that plant ADHs have been involved in seedling, pollen development and fruit development (Kim et al., 2010; Thompson et al., 2010), suggesting that ADH genes play a important role in progress of the vegetative and reproductive organs development. In melon, CmADHs also showed different expression patterns in tissues and fruits during different development stages (Figures 4 and 5; Supplementary Figure S12). CmADH1 and CmADH2 were weakly expressed or not expressed in leaf, stem, seed, root, and flower of Cantaloupe (ManrÍquez et al., 2006), without the same as our results (Figure 4; Supplementary Figure S12). Hence, we speculated that it seems to be related to sample type or other, which are worthy of further study. In mango, all of three MiADHs were weakly expressed in flower, leaf, and stem, with different transcript levels between MiADHs (Singh et al., 2010). In petunia hybrid, three functional ADHs had different expressions in immature pollen grains, anther, stigma, petal, ovary and hypoxic root (Garabagi and Strommer, 2004); Moreover, particularly in the nectary, ADH3 was remarkably expressed, suggesting the association of ADH genes with production of aromatic compounds in the nectary of petunia hybrid (Garabagi et al., 2005). In the present study, CmADH1, CmADH2, CmADH10 and CmFDH1 were strongly expressed in petal, with a higher corresponding enzymatic activity in pistillate and staminate flower petal respectively (Supplementary Figure S15). Therefore, we infered that these four CmADHs participated in the flowers development and the production of volatile substances in melon flowers, and based on these findings, the previously assumed functional redundancy of these genes in the ADH network needs to be questioned (Strommer, 2011).

In addition, plant ADH genes were also implicated in fruit tissue development and fruit ripening, and were closely related to aroma volatile production. For example, mango flesh undergone softening during ripening, and MiADH1 and MiADH2 transcripts particularly accumulated at the onset of ripening in mango fruit, and gradually reduced thereafter; whereas, MiADH3 accumulated during early development of fruit (Singh et al., 2010). In tamato, SlscADH1, a fruit-ripening-associated short-chain ADH, has been related to the formation of aroma volatiles (Gaut et al., 1999). In the present study, the ADH enzyme activity increased at first and then decreased, and had a peak at 36 days after the anthesis during fruit development (Supplementary Figure S16). Nevertheless, only CmADH2 and CmADH10 were strongly expressed during the progress of fruit development except 15 days after the anthesis, suggesting that CmADH2 and CmADH10 may participate in fruit development of the whole development period of melon. In olea europaea, OeADH transcript guadually rose along with olea fruit development and peaked at fruit ripening, and the change of OeADH expression was coherent with the olea fruit growth curve (Iaria et al., 2012), and the similar pattern of ADH gene expression was found in apricot fruits (Gonzalez-Agüero et al., 2009). Additionally, in Cantaloupe, the transcript levels of CmADH1 and CmADH2 gradually increased at first and then decreased from 32 to 42 days with a peak at 39 days (ManrÍquez et al., 2006). However, our data showed CmADH1 was strongly expressed at the late stage of development; CmADH2 transcripts always maintained higher levels during fruit development (Supplementary Figure S12). In addition, VvADH1 and VvADH3 in grape were temporarily expressed in young and developing berry, whereas VvAdh2 transcript was strongly accumulated at later developmental stages, from the onset of ripening (Tesnière and Verriès, 2000). In melon, CmADH1, CmADH3, CmADH8 and CmFDH1 were highly expressed in fruit at 1 day after anthesis, suggesting that these 4 ADH genes likely participate in the ovary development of melon fruit. Of three functional petunia ADHs, PhADH2 and PhADH3 were expressed in ovary (Garabagi et al., 2005). As for melons, little information was available on the role of ADH in relation to the ovary development of fruit. Furthermore, expression analysis also indicated that CmADH8 and CmFDH1 were prominently expressed later in fruit development. However, little information was available on the role of another medium-chain ADH in relation to melon fruit development and ripening. These findings implied that CmADH8 and CmFDH1 might also be involved in fruit development or in other unknown functions at different developmental stages, and also suggested that medium-chain CmADHs might have evolved individual unique functions in fruit development, but it appears that a complex network, consisted of the members of ADH superfamily, plays a important role to modulate the fruit development and ripening. There is much to learn about medium-chain ADHs in melon.

There were significant expression differences between CmADH11 and other medium-chain CmADHs in different tissues and during development and ripening. These observed differences could partly be explained by the amino acids differences at the Zn1-binding sequences and the coenzyme binding position (Supplementary Figure S9). Based on our analysis of medium-chan CmADHs, CmADH11 had little key residues, and had Phe, Ala, and Tyr instead of three Gly in the coenzyme binding domain GXG(X)₂G of all functional medium-chain ADHs. Furthermore, some substituted residues in soybean ADH2 protein might retain the function of amino acids, and it is not clear in other substitutions (Komatsu et al., 2011). Accordingly, these findings implied difference in enzymatic properties between CmADH11 and CmADH1 and other plant ADHs. On the basis of the above results, it appears that CmADH11 is a pseudogene. But we used the expressed sequence tag (EST) database in NCBI as a source of mRNA sequence bioinformatics for CmADH11, and we found that CmADH11 showed the highest homology with a Cucumis melo cDNA clone (GenBank no.JG552191.1) (89% identity), obtained from callus, from EST database of melon (CM-DEa library). Hence, it appears that CmADH11 is specifically expressed in callus. Further research is required to confirm this observation.

Induced Expression Analysis

We also studied the effect of ABA and IAA on CmADHs expression. ABA is known to be extensively involved in plant developmental processes and the plant’s response to abiotic stresses, and can regulate the expression of relevant genes to increase plant adaptability and to regulate tissues development (Skriver and Mundy, 1990; Matsukura et al., 2010; Su et al., 2013). ABA has been shown to induce ADH gene expression in Arabidopsis (de Bruxelles et al., 1996; Su et al., 2013), barley (Macnicol and Jacobsen, 2001) and mango (Singh et al., 2010); In particular, maize (Hwang and Vantoai, 1991) and lettuce (Kato-Noguchi, 2000) ADH genes were strongly induced under abiotic stresses. Purportedly, ABA-mediated plant responses to drought stress might be related to the regulation of relevant genes by key regulator in ABA signaling, including MYB, MYC, NAC, PP2Cs, DREB family and so on (Abe et al., 2003; Fujita et al., 2011; Su et al., 2013). However, to the best of our information, there were little reports on effect of ABA or IAA on ripening specific ADHs in melon. Promoter analysis of the 12 melon ADH genes and one FDH gene suggested the presence of ABA responsive ABRE motif in the promoter of CmFDH1, CmADH2, CmADH3 and CmADH10 (Su et al., 2013) (Supplementary Table S1). Additionally, the promoter of mango ADHs has the ABA responsive ABRE motif (Gomez-Porras et al., 2007; Singh et al., 2010), while MiADH1 was induced by ABA, and MiADH2 was repressed. In the present study, ABA treatment boosted the expressions of CmADH1, CmADH2, CmADH3, CmADH4, CmADH10 and CmFDH1, but ABA did not obviously inhibit other ADHs transcripts (Figure 6; Supplementary Figure S13). Therefore, the ABA-induced these CmADHs expression observed in present study may be related to the upstream key transcription factors. However, this speculation has to be demonstrated by further cloning and function analyses of the promoter of CmADHs. Furthermore, co-treatment with ABA and ethylene enhanced MiADH1 expression which was higher than those by either of the hormone alone, and ethylene could relieve the inhibition of ABA on MiADH2 transcript, suggesting that there were likely mutual promotion between ABA and ethylene, or that ethylene was in some manner able to repress ABA action (Singh et al., 2010). All of these derived results also implied that it would be important to investigate whether ABA triggers different hormonal signal pahways, particularly the possible crosstalk between ABA and ethylene, because ABA and ethylene are key hormones in fruit ripening and senescence, as hormone crosstalk has been shown to be important for biotic and abiotic stresses (Bari and Jones, 2009; Su et al., 2013). More further work is needed to investigate the crosstalk in response to fruit ripening. Auxin (IAA) also plays a role in fruit development and ripening (Singh et al., 2010), and there was an active crosstalk between IAA and ethylene during fruit ripening (Trainotti et al., 2007). However, in the case of melon, there were not reports about the inducing of ripening in mature fruit by IAA as done by ethylene. In tamato, the mutual effect between IAA and ethylene boosted tomato fruit development at early developmental stages (Balbi and Lomax, 2003). Three mango ADH genes expressions were not obviously induced by IAA treatment, and IAA did not induced the mango fruit ripening, moreover there are not IAA responsive TGA motif in promoter in mango ADHs (Singh et al., 2010). CmADH6, CmADH9 and CmADH12 promoters contained IAA responsive TGA motif, and IAA treatment induced these three ADHs transcripts (Figure 6; Supplementary Figure S13). Other ADH genes were expressed at basal level or not expressed, but did not appear to be significantly regulated by IAA, which were consistent with observations in mango (Singh et al., 2010).

Our results also indicated that 12 CmADHs and 1 CmFDH showed different patterns of transcripts during ethylene induced ripening (Figure 7; Supplementary Figure S14), which could be explained by the presence of ethylene responsive ERE motifs in CmADH1, CmADH2, CmADH3, CmADH8, CmADH12 and CmFDH1 (Tesniere et al., 2004; Jin et al., 2014) (Supplementary Table S1). According to previous reports, there are ethylene responsive ERE motifs in Arabidopsis AtADH1 (Peng et al., 2001), mango MiADH1 (Singh et al., 2010) and tomato LeADH2 (Speirs et al., 1998) promoters, and these ADH genes were responsive to ethylene. So, in view of melon ADH genes responses to ethylene, CmADH1, CmADH2, CmADH3, CmADH8, CmADH12 and CmFDH1 were classed into ethylene-sensitive group, and other CmADHs clustered ethylene-insensitive group. Furthermore, Manråquez et al.(2006) also found that ethylene increased the levels of CmADH1 and CmADH2 transcripts in Cantaloupe during fruit ripening, while 1-MCP treatment inhibited the inducing effect of ethylene to CmADH1 and CmADH2. This was similar to the ripening related expression observed in mango ADHs (Singh et al., 2010), tomato LeADH2 (climacteric fruit) (Defilippi et al., 2005b; Barakat et al., 2009) and VvADH2 in grape (non-climacteric fruit) (Tesniere et al., 2004). Simultaneously, ethylene treatment enhanced the activity of ADH enzyme in mango and ethylene production and accelerated mango fruit ripening and senescence, but the increase effect of ethylene was weakened by 1-MCP (Singh et al., 2010). Whereas, in tomato, ADH enzymatic activity increased during fruit ripening, and LeADH2 was involved in fruit ripening and aroma compounds synthesis, but was not regulated by ethylene (Speirs et al., 1998, 2002); Ethylene induced the increases of ADH enzymatic activity and VvADH2 transcript in grape (a climacteric fruit) during fruit ripening (Tesniere et al., 2004), but did not clearly enhanced the activity of ADH enzyme in ACC oxidase antisense transgenic apple, without influence on the esters production in this kind of apple (Defilippi et al., 2005a). These findings showed that the ADH genes superfamily included ethylene non-dependent ADH gene members in climacteric kinds of fruit and ethylene dependent ADH gene members in non-climacteric kinds of fruit. In this paper, the changes of ADH enzyme activity in fruit treated by ethylene and 1-MCP (Supplementary Figure S17) and CmADH10 expression indicated that CmADH10 may be a non-dependent ADH gene member (Figure 7; Supplementary Figure S14). Whereas, more experimental works are needed to further elucidate the function of CmADH10 in aroma compounds synthesis and during fruit ripening and senescence.

Additionally, we also identified several hormone-responsive cis-regulatory elements in the CmADHs promoter region, such as GARE, TATC-box, P-box (gibberellin), TCA-element (SA), CCAAT-box, TGACG-motif and CGTCA-motif (MeJA) (Barakat et al., 2009; Kim et al., 2010; Jin et al., 2014) (Supplementary Table S3). Several previous studies showed a complex interplay of hormones was known to affect fruit development and ripening, with auxin and GA being important during fruit expansion and ABA and ethylene being important for ripening (Srivastava and Handa, 2005; Singh et al., 2010). In addition, oriental melon (Cucumis melo var. makuwa Makino) has distinct ripening patterns within the fruit (Jin et al., 2014). Therefore, based on our derived results in this study, we speculated that the different regulation of melon fruit ripening by ABA, IAA and ethylene might selectively affect the expression of one- or multi- ADHs, which prompted CmADHs to regulate fruit ripening and senescence. Maybe, it will be possible to retain a crosstalk between ABA/IAA/JA/GA and ethylene to regulate the ADH genes network involved in fruit development and ripening, as hormone crosstalk has been shown to be important for biotic and abiotic stresses (Bari and Jones, 2009; Su et al., 2013). Further research is required to confirm this hypothesis.

Conclusion

We identified 13 CmADHs (including 12 ADHs and 1 FDH) in melon genome, exhibiting lower identity in protein level with each other, and they were expressed differentially in different vegetative tissues and during fruit development and ripening. Phylogenetic analysis indicated that they belonged to three different ADH subfamilies and the medium-chain ADHs were clustered into six subgroups. On the transcript level, differential CmADHs expression suggested tight adaptation of the fruit to the developmental events. Promoter sequence analysis, intron-exon structure and the different responses of these ADH genes to ABA, IAA and ethylene implied that ADH genes had different functions or distinct metabolic roles in melon fruit development and ripening, and the different ADH genes expression at transcript induced by ethylene indicated that CmADH10 was an ethylene non-dependent ADH gene member in CmADH genes. Further and detailed research efforts should be made to find the detailed function of these ADH proteins in addition to CmADH1 and CmADH2 during melon fruit development and ripening. Therefore, selective silencing of the CmADH genes and more studies related to biochemical characterization of these CmADHs will be needed.

Author Contributions

Conceived and designed the experiments: YJ, HQ, CZ; Performed the experiments: YJ, CZ, WL, HC, SC; Analyzed the data: YJ, CZ, HQ, HC; Contributed reagents/materials/analysis tools: YJ, CZ, YT; Wrote the paper: YJ, CZ, HQ.

Funding

This work was supported by grants from the National Natural Science Foundation of China (31272154) and the Excellent Talent Project of Liaoning Provincial Department of Education (LR2014020).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgment

We thank Qiang Chen, Deqing Lv, and Tao Xu for excellent technical assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00670

FIGURE S1 | Phylogenetic tree of melon medium-chain ADH proteins. The amino acid sequences were aligned by the Clustal Omega program, and the neighborjoining tree was drawn with TreeView. The corresponding GenBank and the melon genome (https://melonomics.net/) were noted in the phylogenetic tree and the accession number in the melon genome were CmADH1 (MELO3C023685P4), CmADH3 (MELO3C026552P1), CmADH4 (MELO3C027151P1), CmADH5 (MELO3C00579 2P1), CmADH6 (MELO3C026553P1), CmADH7 (MELO3C002189P1), CmADH8 (ME LO3C003251P1), CmADH9 (MELO3C011043P1), CmADH10 (MELO3C026554P2), CmADH11 (MELO3C023687P1), and CmFDH1 (MELO3C022399P1). The number for each interior branch was the percentage of bootstraps value (1000 replicates). Black circle denoted 13 CmADHs. ADHs from other plants in our paper are in Supplementary Table S3.

FIGURE S2 | Amino acid sequence alignment of melon CmADH1 (MELO3C023685P4) and CmADH4 (MELO3C027151 P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S3 | Amino acid sequence alignment of melon CmFDH1 (MELO3C022399P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S4 | Amino acid sequence alignment of melon CmADH3 (MELO3C026552P1), CmADH6 (MELO3C026553P1) and CmADH10 (MELO3C026554P2) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S5 | Amino acid sequence alignment of melon CmADH5 (MELO3C005 792P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S6 | Amino acid sequence alignment of melon CmADH7 (MELO3C002 189P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S7 | Amino acid sequence alignment of melon CmADH8 (MELO3C003 251P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S8 | Amino acid sequence alignment of melon CmADH9 (MELO3C011043 P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S9 | Amino acid sequence alignment of melon CmADH11 (MELO3C023 687P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S10 | Amino acid sequence alignment of melon CmADH2 (MELO3C014 897P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S11 | Amino acid sequence alignment of melon CmADH12 (MELO3C019 503P1) with closely related sequences of other plants. GenBank accession numbers are in Supplementary Table S3.

FIGURE S12 | Transcript abundance of CmADHs in vegetative tissues and developing stages of melon fruit (1–48). R, root; ML, mature leaf; DL, developing leaf; YS, young stem; PF, pistillate flower petal; SF, staminate flower petal. RNA isolated from the samples. Expression analysis was carried out by semi-quantitative PCR using 18S as an internal control.

FIGURE S13 | mRNA abundance of CmADHs in mature green fruit after treatment with ABA and IAA. IAA (Auxin) and ABA (100μM) treatments were given for 2 h as described in Section “Materials and Methods.” Expression analysis was carried out by semi-quantitative PCR as described in section of Hormone treatments using actin as 18S internal control.

FIGURE S14 | CmADHs accumulation during different stages of ethylene and 1-MCP treatment induced ripening of melon by semi-quantitative PCR. 1–12 indicate days after ethylene and 1-MCP treatment. The accumulation level of the genes in untreated melon fruit was control. 18S was used as an internal control.

FIGURE S15 | Alcohol dehydrogenases activity in different melon organs. R, root; ML, mature leaf; DL, developing leaf; YS, young stem; PF, pistillate flower petal; SF, staminate flower petal.

FIGURE S16 | Alcohol dehydrogenases activity in fruits during melon development. 1, 5, 10, 15, 20, 25, 30, 33, 36, 39, 42, 45, and 48 indicate days after anthesis.

FIGURE S17 | Alcohol dehydrogenases activity in fruits during different stages of ethylene and 1-MCP treatment induced ripening of melon. 1, 3, 6, 9, and 12 indicate days after ethylene and 1-MCP treatment. All the treatments have been described in section of Hormone treatments.

Footnotes

References

Abe, H., Urao, T., Ito, T., Seki, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2003). Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 15, 63–78. doi: 10.1105/tpc.006130

PubMed Abstract | CrossRef Full Text | Google Scholar

Alka, K., Windle, H. J., Cornally, D., Ryan, B. J., and Henehan, G. T. M. (2013). A short chain NAD(H)-dependent alcohol dehydrogenase (HpSCADH) from Helicobacter pylori: a role in growth under neutral and acidic conditions. Int. J. Biochem. Cell Biol. 45, 1347–1355. doi: 10.1016/j.biocel.2013.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Balbi, V., and Lomax, T. L. (2003). Regulation of early tomato fruit development by the diageotropica gene. Plant Physiol. 131, 186–197. doi: 10.1104/pp.010132

PubMed Abstract | CrossRef Full Text | Google Scholar

Barakat, A., Bagniewska-Zadworna, A., Choi, A., Plakkat, U., DiLoreto, D., Yellanki, P., et al. (2009). The cinnamyl alcohol dehydrogenase gene family in Populus: phylogeny, organization, and expression. BMC Plant Biol. 9:26. doi: 10.1186/1471-2229-9-26

PubMed Abstract | CrossRef Full Text | Google Scholar

Bari, R., and Jones, J. D. (2009). Role of plant hormones in plant defence responses. Plant Mol. Biol. 69, 473–488. doi: 10.1007/s11103-008-9435-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Borràs, E., Albalat, R., Duester, G., Parés, X., and Farrés, J. (2014). The Xenopus alcohol dehydrogenase gene family: characterization and comparative analysis incorporating amphibian and reptilian genomes. BMC Genomics 15:216. doi: 10.1186/1471-2164-15-216

PubMed Abstract | CrossRef Full Text | Google Scholar

Bukh, C., Nord-Larsen, P. H., and Rasmussen, S. K. (2012). Phylogeny and structure of the cinnamyl alcohol dehydrogenase gene family in Brachypodium distachyon. J. Exp. Bot. 63, 6223–6236. doi: 10.1093/jxb/ers275

PubMed Abstract | CrossRef Full Text | Google Scholar

Çelik, A., and Aktas, F. A. (2013). A new NADH-dependent, zinc containing alcohol dehydrogenase from Bacillus thuringiensis serovar israelensis involved in oxidations of short to medium chain primary alcohols. J. Mol. Catal. B Enzym. 89, 114–121. doi: 10.1016/j.molcatb.2013.01.005