Figures
Abstract
Background
Ripening affects the quality and nutritional contents of fleshy fruits and is a crucial process of fruit development. Although several studies have suggested that ubiquitin-conjugating enzyme (E2s or UBC enzymes) are involved in the regulation of fruit ripening, little is known about the function of E2s in papaya (Carica papaya).
Methodology/Principal findings
In the present study, we searched the papaya genome and identified 34 putative UBC genes, which were clustered into 17 phylogenetic subgroups. We also analyzed the nucleotide sequences of the papaya UBC (CpUBC) genes and found that both exon-intron junctions and sequence motifs were highly conserved among the phylogenetic subgroups. Using real-time PCR analysis, we also found that all the CpUBC genes were expressed in roots, stems, leaves, male and female flowers, and mature fruit, although the expression of some of the genes was increased or decreased in one or several specific organs. We also found that the expression of 13 and two CpUBC genes were incresesd or decreased during one and two ripening stages, respectively. Expression analyses indicates possible E2s playing a more significant role in fruit ripening for further studies.
Citation: Jue D, Sang X, Shu B, Liu L, Wang Y, Jia Z, et al. (2017) Characterization and expression analysis of genes encoding ubiquitin conjugating domain-containing enzymes in Carica papaya. PLoS ONE 12(2): e0171357. https://doi.org/10.1371/journal.pone.0171357
Editor: Sara Amancio, Universidade de Lisboa Instituto Superior de Agronomia, PORTUGAL
Received: August 23, 2016; Accepted: January 18, 2017; Published: February 23, 2017
Copyright: © 2017 Jue et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the Natural Science Foundation of Hainan Province (20163111), the Natural Science Foundation of China (31572087), the Natural Science Foundation of Guangdong Province of China (2014A030307031), the Basic Scientific Research Project of Non-profit Central Research Institutions (No. 1630062015004), the Germplasm Resource Protection of Ministry of Agriculture (2016NWB053), the Germplasm Resource Protection of Ministry of Agriculture Wild Plant (1251416305010), the Science and Technology Project of Guangdong Province (2016A020208002) and the Nature Science Foundation of Anhui Province (grant no. 1408085MKL63).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Papaya (Carica papaya) is an economically important fruit crop that is widely grown in tropical and sub-tropical regions, and its fleshy fruits play an important role in the human food supply, as well as in human nutrition and health [1, 2]. Gaining in popularity among tropical fruits worldwide, papaya is now ranked fourth in total tropical fruit production after bananas, oranges), and mango [3]. Fruit ripening is an important and complex process, and, therefore, the regulatory mechanisms of fruit ripening have been intensively studied [4]. Fruit with different ripening mechanisms can be divided into two groups: climacteric fruit (e.g., bananas, tomatoes, pears, mangos, and papayas), in which ripening is accompanied by a peak in respiration and a concomitant burst of ethylene, and non-climacteric fruit (e.g., pineapple, carambola, borojo, and watermelon), in which respiration exhibits no dramatic change and ethylene production remains at a very low level [5].
Ethylene, in particular, has been extensively studied, owing to its crucial role in the ripening of climacteric fruits [6–9], and the biosynthesis, perception, and signal transduction of ethylene, as well as downstream gene regulation, are well known [10]. Tomato (Solanum lycopersicum) has long served as an excellent model for studies of fleshy fruit development and ripening [11, 12], and to date, many of the underlying genes have been cloned, thus providing insight into the regulation of ripening. The ripening-associated genes that have been reported in tomato include: rin, which encodes a MADS-box transcription factor [13]; Nr, which encodes an ethylene receptor [14]; Cnr (Colorless nonripening), which encodes an SBP-box transcription factor [15]; Gr (Green ripe), which encodes a still poorly defined component of the ethylene transduction pathway [6]; SINAC4, a new tomato NAC transcription factor that positively regulates fruit ripening and carotenoid accumulation [16]; APETALA2a (AP2a), a transcription factor that regulates fruit ripening via regulation of ethylene biosynthesis and signaling [17]; LeNCED1, which initiates abscisic acid (ABA) biosynthesis at the onset of fruit ripening and might act as an original inducer in tomato [18]; and the tomato genes SHATTERPROOF, TAGL1 [19], TDR4/FUL1, and MBP7/FUL2 [20], which are also involved in regulating fruit ripening. However, in contrast to tomato, there have been only a few reports regarding fruit ripening in C. papaya.
Ubiquitination is an essential cellular process in all eukaryotes, from unicellular yeast to humans, and the ubiquitination-proteasomal pathway has been implicated in diverse aspects of eukaryotic cellular regulation, owing to its ability to degrade intracellular proteins [21, 22]. The process of protein ubiquitination is mediated through the action of three enzymes, which are known as ubiquitin-activating enzymes (E1s), ubiquitin-conjugating enzymes (E2s), and ubiquitin ligases (E3s) [23]. Ubiquitin-conjugating enzymes (E2s) function at an intermediate step in the protein ubiquitination pathway. The E2 family has expanded during evolution, and more ancestral eukaryotes possess fewer E2 enzymes (e.g., ≤20 in algae) than more derived ones (e.g., >40 in certain plants and animals) [24]. For example, 13 UBC domain proteins have been identified in Saccharomyces cerevisiae [25], 12, 18, and 19 E2 enzymes are encoded in the algae Ostreococcus tauri, Micromonas sp. RCC299, and Chlamydomonas reinhardtii, respectively (data not shown), 20 in Caenorhabditis elegans [26], 37 in both humans [24] and Arabidopsis thaliana [27], 48, 52, 74 and 75 UBC domain proteins in rice, tomato, banana and maize [28–31].
Previous studies have suggested that ubiquitination plays an important role in plant growth, development, and tolerance against various abiotic stresses, and several of the UBC genes in A. thaliana have been implicated specifically [32, 33]. For example, the A. thaliana ubiquitin-conjugating gene AtUBC13 has been implicated in epidermal cell differentiation and iron deficiency responses [34, 35], and AtUBC32 is an endoplasmic reticulum-associated degradation (ERAD) component that functions in brassinosteroid-mediated salt stress tolerance [36], whereas AtUBC21 (AtPEX4) is specialized for ubiquitination in peroxisome maintenance [37]. In rice, 14 OsUBCs were differentially expressed under drought and salt stress conditions [38], and in Zea mays, 48 and 16 ZmUBC genes were significantly upregulated under salt and drought stress conditions, respectively [30].
Several studies suggest that E2s also participate in the fruit ripening process. In tomato, for example, SlUBC32, which is upregulated during tomato fruit ripening and downregulated in the rin mutant, plays an important role in the regulation of fruit ripening [31]; and in banana, 32 MaUBC genes are reportedly up- and downregulated during different ripening stages [29]. However, functional studies of E2s involved in fruit ripening are limited and have not been conducted in papaya. Therefore, in the present study, we performed a genome-wide analysis of E2 domains, analyses of CpUBC gene structures and phylogenetic relationships, and expression profiling of CpUBC genes in various C. papaya plant tissues and at different stages of fruit ripening. To the best of our knowledge, this is the first reported genome-wide analysis of the papaya UBC gene family, and the results will facilitate further investigation of the roles of the UBC genes in fruit ripening and will aide in the functional validation of UBC genes in papaya.
Materials and Methods
2.1. Plant material and treatments
Papaya (C. papaya cv. ‘sun up’) plants were obtained from the South Subtropical Crops Research Institute of the Chinese Academy of Tropical Agricultural Science (Zhan-jiang, Guangdong Province, China), and root, stem, leaf, male flower, female flower, and fruit (mature green) organs were collected separately. For the fruit ripening experiment, we also collected fruit pulp from four development stages: immature green (IG), mature green (MG), breaker (Br), and mature fruit (MF). The experiments were performed in triplicate, and all samples were immediately frozen in liquid nitrogen and stored at −80°C for expression analysis.
2.2. Sequence retrieval and gene family member identification
To identify potential members of the papaya UBC protein family, published Arabidopsis and yeast (S. cerevisiae) UBC protein sequences were retrieved and used as queries in BLASTP searches against the C. papaya Genomics Database (CpGDB; http://www.plantgdb.org/CpGDB/) and PLAZA 3.0 (http://bioinformatics.psb.ugent.be/plaza/versions/plaza/) [39]. Amino acid sequences of all the UBC-encoding genes of Arabidopsis, S. cerevisiae, and S. lycopersicum were downloaded from the Arabidopsis Information Resource (http://www.arabidopsis.org/), Saccharomyces Genome Database (http://www.yeastgenome.org/), and Sol Genomics Network Tomato Database (http://solgenomics.net) websites, respectively. All putative papaya UBC protein sequences were then subject to domain analysis using Simple Modular Architecture Research Tool (SMART; http://smart.emblheidelberg.de/) [40], and information regarding the CpUBC genes, including location, open reading frame (ORF), amino acid (AA) length and intron number, were obtained from PLAZA 3.0 [39]. The molecular weight (MW) and theoretical isoelectric point (PI) of the candidate CpUBC proteins were investigated using ExPASy online tools (http://expasy.org/tools/) [41]. Finally, the structure of each gene was visualized using the Gene Structure Display Server (GSDS) (http://gsds.cbi.pku.edu.cn/) [42].
2.3. Sequence alignment and phylogenetic analysis
To evaluate the evolutionary relationship of the CpUBC proteins, full-length amino acid sequences of 34 CpUBC proteins, 15 UBC proteins from S. cerevisiae, 48 UBC proteins from Arabidopsis thaliana, and 52 UBC proteins from S. lycopersicum were subject to multiple sequence alignment, using MUSCLE [43], with default parameters, and manual optimization. We excluded ambiguously aligned sequences, in order to produce an alignment of 34 amino acid characters for subsequent phylogenetic analyses. An unrooted phylogenetic tree was constructed using the neighbor-joining (NJ) method in MEGA (Version 6.0; http://www.megasoftware.net/), and the significance of nodes was assessed using a bootstrap test with 1,000 iterations [44]. Representations of the calculated trees were constructed using TreeView (Version 1.6.6; http://taxonomy.zoology.gla.ac.uk/rod/treeview).
2.4. RNA isolation and expression analysis
Total RNA was extracted for cDNA synthesis using the Super fast new plants of RNA extraction kit (Huayueyang Bio Co., Ltd, Beijing, China), following the manufacturer's instructions. Reverse transcription reactions were performed using the PrimeScript RT reagent kit with gDNA Eraser (Takara Bio, Inc., Kusatsu, Japan), according to the supplier’s manual. Real-time PCR was performed using a Bio-Rad real-time thermal cycling system (LightCycler 480; Bio-Rad Laboratories, Inc., Hercules, CA, USA) and SYBR-green to assess the expression levels of 34 candidate CpUBC genes. The gene-specific primers were designed according to the CpUBC gene sequences using Primer 5 software and checked using Blast (S1 Table) In addition, the papaya actin 1 gene was used as an internal control for normalization. Each reaction consisted of 10 μL 2× SYBR Premix Ex Taq II (Takara Bio), 40 ng cDNA, and 250 nM of each primer, in a final volume of 20 μL. The following PCR program was used: 94°C for 10 min, followed by 40 cycles of 94°C for 10 s, 58–63°C for 20 s, and 72°C for 30 s. The relative mRNA levels of the genes were measured using the cycle threshold (Ct) 2(-ΔCt) method. The analysis included cDNA from the three biological samples for each tissue, and all the reactions were run in triplicate. In the comparative expression analysis of CpUBCs, genes that were up- or downregulated by at least two-fold were considered differentially expressed.
Results
Identification of UBC genes in papaya
In the present study, we used the published Arabidopsis and yeast (S. cerevisiae) UBC protein sequences as queries in BLASTP searches against the CpGDB and PLAZA 3.0[39], a total of 39 putative CpUBC genes were identified. After scanning of the 39 sequences for the UBC domain by motif scan using SMART search, we found only 34 sequences contain the UBC domain. According to their chromosome locations (Table 1), the 34 CpUBC genes were designated CpUBC1 to CpUBC34. In these 34 CpUBCs, two were identified as Related to Ubiqutin (RUB) conjugating enzymes (CpUBC18 and CpUBC29), and one was identified as a SUMO-conjugating enzyme (CpUBC11). Six other UBC proteins (CpUBC2/14/16/19/24/25) lacked the Cys active site, identifying them as ubiquitin-conjugating enzyme variants (UEVs), which are not active by themselves, thus leaving 25 potential ubiquitin E2s. To better understand the papaya E2s, we constructed a recombinant E2 protein library of all 34 CpUBCs. The predicted proteins ranged from 98 (CpUBC15) to 668 amino acids (CpUBC20) in length, with corresponding molecular masses of 10.74 kDa and 74.03 kDa, respectively, and predicted isoelectric points of between 4.29 (CpUBC33) and 9.64 (CpUBC27).
Structure and phylogenetic analysis of CpUBC genes
Using the GSDS website, we found that the number of introns in the 34 CpUBC genes ranged from zero (CpUBC8, CpUBC14 and CpUBC19) to seven (CpUBC20 and CpUBC23), with most of CpUBC genes containing three to five introns (n = 76.5%) (Table 1 and Fig 1). In addition, most CpUBC genes within the same subfamilies shared the same exon/intron structure. For example, in the UBC4/5 subfamily, CpUBC3, CpUBC4, CpUBC7, CpUBC9, and CpUBC26 contained three introns, whereas in the UBC11 subfamily, CpUBC17 and CpUBC21 harbored five introns.
Untranslated 5′ and 3′ regions, exons, and introns are indicated by red, blue, and gray, respectively. The scale bar represents 1000 bp.
Multiple sequence alignment of the predicted amino acid sequences revealed that the tryptophan (W) located at the C-terminal side of the active cysteine was conserved in most of the CpUBC genes (Fig 2A). The consensus active site motif “HPN” was found at six amino acids from the N-terminal site of the active cysteine, and the strongly conserved PxxPP motif was found at seven amino acids from the N-terminal side of the HPN motif (Fig 2A). In addition, the predicted amino acid sequences of the UBC domains of 34 CpUBC genes were analyzed using the MEME Suite website (http://meme-suite.org/index.html). The results indicated that the highly conserved sequence of the UBC domain is HPNINSNGSICLDILKEQWSP (Fig 2C).
(A) Multiple sequence alignment of CpUBC domains. Identical amino acids, conserved amino acids, and blocks of similar amino acid residues are shaded in black, charcoal gray, and gray, respectively. (B) Structures of the CpUBC proteins. The name of each corresponding protein is shown on the left. The position of the UBC domain is indicated in the figure. The length and order of the domains represent their actual location within each protein. The scale bar represents 100 amino acids. The different colors indicate the four E2 subtypes of UBC domains. (C) Conserved CpUBC domains. Logos of the protein alignment of the UBC domain is shown. The x-axis represents the conserved sequences, and the conservation of each residue is indicated by the height of each letter. The y-axis is a scale of the relative entropy, which reflects the conservation rate of each amino acid.
According to the UBC domain and the N- or C-terminal structure, E2 proteins are divided into four types. Class I E2s contain only the catalytic domain; Class II E2s contain N-terminal extensions; Class III E2s contain C-terminal extensions; and Class IV E2s have both N- and C-extensions [45, 46]. In the present study, we found that 23 of the CpUBCs belonged to Class I, and three, seven, and one belonged to Class II, Class III, and Class IV, respectively (Fig 2B).
Furthermore, phylogenetic analysis revealed that the UBC proteins could be divided into 14 E2 groups and three independent UEV groups (UBC17, UBC18, and UBC19, with UBC9 and UBC12 functioning in the SUMO and RUB1 conjugation pathways, based on >50% bootstrap support (Fig 3). Most of the CpUBC groups were named according to the identities of the S. cerevisiae proteins or proteins contained within the group (UBC1, UBC2, UBC3, UBC4/5, UBC6, UBC3/7, UBC8, UBC9, UBC10, UBC11, UBC12, UBC13, UBC17, and UBC18); however, there were no yeast proteins in groups UBC14, UBC15, UBC16, and UBC19, which indicated that the groups may be plant-specific or were lost in yeasts. Interestingly, the subfamilies UBC4/5 or UBC3/7 shared two highly identical paralogous yeast genes, respectively. For example, during the phylogenetic analysis of CpUBC genes using MEGA 5.0, UBC4 and UBC5 always clustered into the same subgroup, as did UBC3 and UBC7. According to previous reports [26, 27, 47], we designated the two clades as subfamilies UBC4/5 and UBC3/7.
The neighbor-joining tree includes 34 UBC proteins from papaya (squares), 15 from S. cerevisiae (rhombi), 48 from A. thaliana (circles), and 52 from S. lycopersicum (triangles). The different colors indicate different CpUBC subgroups.
Organ- and ripening stage-specific expression of CpUBC genes
All 34 of the candidate CpUBC genes were differentially expressed among the six papaya organs (Fig 4). In addition, we found that the expression of 13 (CpUBC1/4/5/6/9/10/13/17/18/24/25/26/34), 11 (CpUBC2/3/15/21/23/27/28/29/30/31/33), and seven (CpUBC7/8/11/12/14/16/20) genes were lowest in leaves, roots, and stems, respectively, and that the expression of only one gene was lowest in the male flowers (CpUBC32) and female flowers (CpUBC19), respectively. Meanwhile, 19 (CpUBC1/2/3/5/ 6/9/10/11/12/15/17/20/23/24/26/30/31/33/34) and two genes (CpUBC21 and CpUBC22) were highly expressed in male flowers and female flowers, respectively, which implicated their involvement in the development of floral sex organs. Four (CpUBC25/27/28/29), five (CpUBC7/8/14/16/19), and two (CpUBC4 and CpUBC32) genes were highly expressed in fruits, leaves, and stems, respectively, and a single gene (CpUBC13) was highly expressed in roots. Therefore, our results suggested that CpUBC genes play multiple roles in papaya development.
The x-axis indicates time course/treatment, and the y-axis indicates the relative expression level. Error bars indicate standard deviations of independent biological replicates (n ≥ 3). R, root; St, stem; L, leaf; T, tassel, Ys, young seed; and Si, silk.
In the present study, we also analyzed the expression of the 34 CpUBC genes at different fruit ripening stages. We found that 13 (CpUBC4/6/7/8/9/11/12/14/16/19/20/28/34) and two (CpUBC2 and CpUBC10) of the genes showed higher or lower expression during the progress of papaya fruit ripening, respectively, whereas the expression of the remaining genes was not significantly different (Fig 5). Of the 13 genes that were higher expressed in the fruit, the expression levels of eight genes (CpUBC7/8/9/11/12/14/19/20) were slightly increased at the IG, MG, and Br stages and most highly expressed at the MF stage. Meanwhile, the expression of four genes (CpUBC6/16/28/34) continuously increased during the first three stages of fruit ripening, before being decreased at the MF stage, and the expression of CpUBC4 was increased at MG stage and then reduced at the two late stages. These results suggest that the papaya E2s maybe playing a significant role in fruit development and ripening.
The x-axis indicates time course/treatment, and the y-axis indicates the relative expression level. Error bars indicate standard deviations of independent biological replicates (n ≥ 3). IG, immature green; MG, mature green; Br, breaker; and MF, mature fruit.
Discussion
Ripening affects the quality and nutritional contents of fleshy fruits and is a crucial process in the development of climacteric fruits [48]. In recent years, the study of genes involved in fruit ripening has rapidly progressed; however, most research has been conducted in tomato. In addition, E2 proteins are critically important in many aspects of plant growth and development; however, considering the potential functional significance of E2 proteins, only a few E2 family members have been described in higher plants[27, 31, 38], and only a few studies have investigated the role of the UBC gene family in fruit ripening [29, 31].
In the present study, we identified genes encoding 34 C. papaya UBC (CpUBC) proteins (Table 1), which includes two RUB-conjugating enzymes (CpUBC18 and CpUBC29), one SUMO-conjugating enzyme (CpUBC11), and six UEV proteins (CpUBC2/14/16/19/24/25). Similar results have also been reported by other studies. For example, 16 UBC proteins were identified in S. cerevisiae, including one RUB-conjugating enzyme (ScUBC9), one SUMO-conjugating enzyme (ScUBC12), and two UEV protein (ScMMS2, NP_011428 and ScSTP22, NP_009919) [49]. In addition, a total of 48 UBC domain-containing proteins have been identified in Arabidopsis, of which two are RUB-conjugating enzymes (RCE1, At4g36800 and RCE2, At2g18600), one is a SUMO-conjugating enzyme (AtSCE1, At3g57870), and eight others are UEVs, thus leaving 37 potential E2s [27]. Among 48 rice E2s, three proteins (OsUBC1/2/3) were identified as RUB E2s, three SUMO E2s (OsUBC4/5/6), and four UEVs (OsUBC28/29/30/31) [28]. Meanwhile, of the 52 SlUBC proteins in tomato, there are three SUMO E2s (SlUBC8/23/51) and five RUB E2s (SlUBC10/13/18/29/50) [31]; and of the ZmUBC proteins in maize, there are three RUB E2s (ZmUBC-37/40/43), eight SUMO E2s (ZmUBC-01/04/11/12/51/69/73/75), and six UEVs (ZmUBC-08/04/31/34/53/74) [30]. The number of CpUBC genes is far less than the 48, 48, 52, 75, and 74 UBC genes reported from the Arabidopsis [27], rice [28], tomato [31], maize [30], and banana genomes [29], respectively. The genome sizes of Arabidopsis, rice, tomato, maize, and banana are ~125, ~466, ~900, ~2300, and 523 Mb, respectively. However, although the genome sizes of rice and tomato are 3.7 and 7.2 times larger than that of Arabidopsis, respectively, the genomes of the three species contain a similar number of UBC genes (45, 48, and 52, respectively); and although the genome size of maize is 4.3 times larger than that of banana, the genomes of these two species also contain a similar number of UBC genes (75 and 74, respectively). Therefore, the difference in number of UBC genes among these species cannot be explained by genome size.
Alternatively, gene duplication events have been demonstrated to play a significant role in the expansion of gene family members in genomes [38], and an increasing number of studies has shown that segmental duplications are largely responsible for the expansion of gene families in maize, such as the CCCH, HD-Zip HSF, bZIP, and PRX gene families[50–53]. Research has also shown that segmental duplications are the main contributor to the expansion of UBC genes in maize and banana [29, 30]. However, owing to a lack of information regarding the chromosomes of papaya, we were unable to perform segmental and tandem duplication analysis in papaya.
Through systematic phylogenetic analyses, the present study provides a detailed classification and nomenclature of papaya E2s. The phylogenetic analysis categorized the 34 CpUBCs into 17 discrete groups, which was similar to the number of groups described in other species, such as tomato, maize, and banana [29–31]. However, there were still some minute differences between the topologies of the UBC genes in papaya and the other species. For example, there are 15 and 13 UBC groups in tomato and banana, respectively, not including the UEV subgroups [29, 31]. However, in Arabidopsis, rice, and maize, the UEV subgroups were taken into account [27]. In our study, to better understand the papaya E2s, the CpE2s were divided into 17 groups, which included three independent UEV groups. The corresponding homologs of AtUBC21 and AtUBC22 were grouped together in maize [30], banana [29], and papaya, whereas the corresponding homologs (and orthologs) were separated into two groups in tomato [31]. Additionally, homologs of AtUBC31, grouped together in the UBC4/5 group in maize and banana [29, 30], were not grouped in any groups in our study. These differences could have resulted from different parameter settings during the phylogenetic analyses.
In addition, the number of UBC genes was obviously different among the groups. The largest groups (UBC4/5) included eight members, whereas the UBC1, UBC2, UBC6, UBC9, UBC10, UBC13, UBC14, and UBC15 groups were only represented by one member each. The UBC4/5 groups was also the largest in other species, which indicates that, in plants, the group may possess more diverse functions than other groups. Phylogenetic data also suggested that the UBC9 and UBC12 groups were expanded in monocots but not in Arabidopsis, and previous research has shown that the genomes of maize and rice each contain eight and three members of the UBC9 and UBC12 groups, respectively [28, 30]. However, in papaya, we only identified one and two members in these groups. Furthermore, we also identified three UBC groups (UBC14, UBC15, and UBC16) in Arabidopsis and papaya for which no probable homologs exist in budding yeast; however, all three groups have potential homologs in animals [54], which indicates a possible gene loss during yeast evolution.
Although the genome sequence of papaya has been reported, the identification and functional studies of papaya genes have proceeded at a slow pace. Analysis of the temporal and spatial expression patterns of genes may provide useful information for characterizing their functions. Previous studies have shown that UBC genes have different expression patterns in different organs. In Arabidopsis, for example, AtUBC1 and AtUBC2 are ubiquitously expressed in roots, leaves, flowers, and seedlings, and the double mutant atubc1-1atubc2-1 exhibits a dramatically reduced number of rosette leaves and an early-flowering phenotype [55]. In banana, MaUBC10/11/33/34/61 are highly expressed in most organs, especially in roots, stems, leaves, and MaUBC6/11/34/35/45/61 were highly expressed in stems. In addition, 12 MaUBC genes (MaUBC13/18/29/33/34/36/43/46/48/53/67/70) were predominantly expressed in roots [29]. Meanwhile, in the present study, we found that 19 and two genes were highly expressed in male and female flowers, respectively, which suggests that the genes may be involved in the development of floral sex organs, and various other CpUBC genes were specifically higher expressed in other tissues, which that CpUBC genes play multiple roles in the development of papaya.
Previous studies have demonstrated that ubiquitin-conjugating enzymes (E2s) are critically important in many aspects of plant growth and development, as well as in physiological processes, such as stress responses [56]. For example, OgUBC1 from wild rice is involved in cellular responses to biotic and abiotic stresses [56], and three (OsUBC2/5/18) and five UBC genes (AtUBC13/17/20/26/31) from rice and Arabidopsis, respectively, are significantly reduced in response to salt and drought stress, whereas only three rice genes (OsUBC13/15/45) are significantly upregulated [38]. Although the function of E2s in plant development and stress responses have been well clarified, reports about the function of protein ubiquitination in fruit ripening remains extremely rare, except in tomato and banana [29, 31]. In tomato, a total of six UBC genes (SlUBC6/8/24/32/41/42) were directly regulated by the fruit-ripening regulator RIN, which suggests that specific UBC genes might be involved in fruit ripening, and in banana, five UBC genes (MaUBC1/9/70/68/71) exhibited approximately 10- to 40-fold higher levels of expression at the fifth stage than at other stages of fruit ripening, whereas seven other UBC genes (MaUBC8/16/17/33/34/56/61) exhibited continuously increasing expression throughout the fruit development. In the present study, 13 (CpUBC4/6/7/8/9/11/12/14/16/19/20/28/34) and two (CpUBC2 and CpUBC10) of the 34 CpUBC genes were up- or downregulated during the progression of papaya fruit ripening, respectively (Fig 5), and our results indicated that the expression patterns of some of the CpUBC genes, such as CpUBC17 and CpUBC56, was related to their phylogenetic relationships. By contrast, the expression patterns of some paralogs, such as CpUBC8 and CpUBC23, were quite different. Meanwhile, we also found that the genes with close phylogenetic relationships that exhibited similar expression patterns in banana, tomato, and papaya. For example, our data indicated that the CpUBC23 orthologs MaUBC3/8 and SlUBC32 were similarly expressed during fruit ripening. However, CpUBC23 was not considered differentially expressed genes in the present study since its expression was upregulated by less than two-fold. MaUBC1, MaUBC9, SlUBC41, SlUBC42, and CpUBC34, which belong to the UBC8 subgroup, also exhibited similar expression patterns during fruit ripening. Therefore, taken together, our results suggest that the papaya E2 family genes might be participated in the regulation of fruit development and ripening processes in papaya.
Conclusions
In the present study, we describe the genome-wide identification and analysis of UBC genes in papaya. A total of 34 putative CpUBC genes were identified, and phylogenetic analysis indicated that genes could be divided into 17 subfamilies. Analysis of exon-intron junctions and sequence motifs revealed high levels of conservation within and between phylogenetic groups. In addition, all the CpUBC genes were detected in roots, stems, leaves, male and female flowers, and mature fruit, although some genes showed higher or lower expression in one or several specific organs. Similar to the results in tomato and banana, we also found that the expression of 13 and two CpUBC genes were incresesd or decreased during one and two ripening stages, respectively. Therefore, the results of the present study suggest that CpUBC genes are involved in the regulation of fruit development and ripening processes. The results also provide novel insights into the function of plant UBC genes and will facilitate further investigation of the roles of the UBC genes in fruit ripening.
Acknowledgments
This work was supported by the Natural Science Foundation of Hainan Province (20163111), the Natural Science Foundation of China (31572087), the Natural Science Foundation of Guangdong Province of China (2014A030307031), the Basic Scientific Research Project of Non-profit Central Research Institutions (No. 1630062015004), the Germplasm Resource Protection of Ministry of Agriculture (2016NWB053), the Germplasm Resource Protection of Ministry of Agriculture Wild Plant (1251416305010), the Science and Technology Project of Guangdong Province (2016A020208002) and the Nature Science Foundation of Anhui Province (grant no. 1408085MKL63).
Author Contributions
- Conceptualization: DWJ SYS.
- Data curation: DWJ XLS.
- Formal analysis: DWJ XLS.
- Funding acquisition: DWJ SYS.
- Investigation: YCW YZ.
- Methodology: DWJ SYS.
- Project administration: DWJ SYS.
- Resources: DWJ ZWJ.
- Software: BS.
- Supervision: LQL.
- Validation: DWJ SYS.
- Visualization: DWJ.
- Writing – original draft: DWJ.
- Writing – review & editing: DWJ.
References
- 1. Devitt LC, Fanning K, Dietzgen RG, Holton TA. Isolation and functional characterization of a lycopene β-cyclase gene that controls fruit colour of papaya (Carica papaya L.). Journal of Experimental Botany. 2010;61(1):33–9. pmid:19887502
- 2. Ming R, Hou S, Feng Y, Yu Q, Dionnelaporte A, Saw JH, et al. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008;452(7190):991–6. pmid:18432245
- 3. Evans EA, Ballen FH. An Overview of Global Papaya Production, Trade, and Consumption. Food & Resource Economics. 2015.
- 4. Giovannoni JJ. Genetic regulation of fruit development and ripening. The plant cell. 2004;16(suppl 1):S170–S80.
- 5. White PJ. Recent advances in fruit development and ripening: an overview. Journal of Experimental Botany. 2002;53(377):1995–2000.
- 6. Barry CS, Giovannoni JJ. Ripening in the tomato Green-ripe mutant is inhibited by ectopic expression of a protein that disrupts ethylene signaling. Proceedings of the National Academy of Sciences. 2006;103(20):7923–8.
- 7. Hackett RM, Ho C-W, Lin Z, Foote HC, Fray RG, Grierson D. Antisense inhibition of the Nr gene restores normal ripening to the tomato Never-ripe mutant, consistent with the ethylene receptor-inhibition model. Plant Physiology. 2000;124(3):1079–86. pmid:11080285
- 8. Kevany BM, Tieman DM, Taylor MG, Cin VD, Klee HJ. Ethylene receptor degradation controls the timing of ripening in tomato fruit. The Plant Journal. 2007;51(3):458–67. pmid:17655616
- 9. Oeller PW, Min-Wong L. Reversible inhibition of tomato fruit senescence by antisense RNA. Science. 1991;254(5030):437. pmid:1925603
- 10. Lin Z, Zhong S, Grierson D. Recent advances in ethylene research. Journal of experimental botany. 2009;60(12):3311–36. pmid:19567479
- 11. Consortium TG. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485(7400):635–41. pmid:22660326
- 12. Moore S, Vrebalov J, Payton P, Giovannoni J. Use of genomics tools to isolate key ripening genes and analyse fruit maturation in tomato. Journal of Experimental Botany. 2002;53(377):2023–30. pmid:12324526
- 13. Vrebalov J, Pan IL, Arroyo AJM, McQuinn R, Chung M, Poole M, et al. Fleshy fruit expansion and ripening are regulated by the tomato SHATTERPROOF gene TAGL1. The Plant Cell. 2009;21(10):3041–62. pmid:19880793
- 14. Wilkinson JQ, Lanahan MB, Yen H-C, Giovannoni JJ, Klee HJ. An ethylene-inducible component of signal transduction encoded by Never-ripe. Science. 1995;270(5243):1807. pmid:8525371
- 15. Manning K, Tör M, Poole M, Hong Y, Thompson AJ, King GJ, et al. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nature genetics. 2006;38(8):948–52. pmid:16832354
- 16. Zhu M, Chen G, Zhou S, Tu Y, Wang Y, Dong T, et al. A new tomato NAC (NAM/ATAF1/2/CUC2) transcription factor, SlNAC4, functions as a positive regulator of fruit ripening and carotenoid accumulation. Plant and Cell Physiology. 2014;55(1):119–35. pmid:24265273
- 17. Karlova R, Rosin FM, Busscher-Lange J, Parapunova V, Do PT, Fernie AR, et al. Transcriptome and metabolite profiling show that APETALA2a is a major regulator of tomato fruit ripening. The Plant Cell. 2011;23(3):923–41. pmid:21398570
- 18. Zhang M, Yuan B, Leng P. The role of ABA in triggering ethylene biosynthesis and ripening of tomato fruit. Journal of Experimental Botany. 2009;60(6):1579–88.
- 19. Itkin M, Seybold H, Breitel D, Rogachev I, Meir S, Aharoni A. TOMATO AGAMOUS‐LIKE 1 is a component of the fruit ripening regulatory network. The Plant Journal. 2009;60(6):1081–95. pmid:19891701
- 20. Bemer M, Karlova R, Ballester AR, Tikunov YM, Bovy AG, Wolters-Arts M, et al. The tomato FRUITFULL homologs TDR4/FUL1 and MBP7/FUL2 regulate ethylene-independent aspects of fruit ripening. The Plant Cell. 2012;24(11):4437–51. pmid:23136376
- 21. Callis J, Vierstra RD. Protein degradation in signaling. Current opinion in plant biology. 2000;3(5):381–6. pmid:11019805
- 22.
Hershko A. The ubiquitin system. Ubiquitin and the Biology of the Cell: Springer; 1998. p. 1–17.
- 23. Mukhopadhyay D, Riezman H. Proteasome-independent functions of ubiquitin in endocytosis and signaling. Science. 2007;315(5809):201–5. pmid:17218518
- 24. van Wijk SJ, Timmers HM. The family of ubiquitin-conjugating enzymes (E2s): deciding between life and death of proteins. The FASEB Journal. 2010;24(4):981–93. pmid:19940261
- 25. Michelle C, Vourc’h P, Mignon L, Andres CR. What was the set of ubiquitin and ubiquitin-like conjugating enzymes in the eukaryote common ancestor? Journal of molecular evolution. 2009;68(6):616–28. pmid:19452197
- 26. Jones D, Crowe E, Stevens TA, Candido EPM. Functional and phylogenetic analysis of the ubiquitylation system in Caenorhabditis elegans: ubiquitin-conjugating enzymes, ubiquitin-activating enzymes, and ubiquitin-like proteins. Genome biology. 2001;3(1):1.
- 27. Kraft E, Stone SL, Ma L, Su N, Gao Y, Lau O-S, et al. Genome analysis and functional characterization of the E2 and RING-type E3 ligase ubiquitination enzymes of Arabidopsis. Plant physiology. 2005;139(4):1597–611. pmid:16339806
- 28. Bae H, Kim WT. Classification and interaction modes of 40 rice E2 ubiquitin-conjugating enzymes with 17 rice ARM-U-box E3 ubiquitin ligases. Biochemical and biophysical research communications. 2014;444(4):575–80. pmid:24486490
- 29. Dong C, Hu H, Jue D, Zhao Q, Chen H, Xie J, et al. The banana E2 gene family: Genomic identification, characterization, expression profiling analysis. Plant Science. 2016;245:11–24. pmid:26940488
- 30. Jue D, Sang X, Lu S, Dong C, Zhao Q, Chen H, et al. Genome-wide identification, phylogenetic and expression analyses of the ubiquitin-conjugating enzyme gene family in maize. PloS one. 2015;10(11):e0143488. pmid:26606743
- 31. Wang Y, Wang W, Cai J, Zhang Y, Qin G, Tian S. Tomato nuclear proteome reveals the involvement of specific E2 ubiquitin-conjugating enzymes in fruit ripening. Genome biology. 2014;15(12):1.
- 32. Sadanandom A, Bailey M, Ewan R, Lee J, Nelis S. The ubiquitin–proteasome system: central modifier of plant signalling. New Phytologist. 2012;196(1):13–28. pmid:22897362
- 33. Welchman RL, Gordon C, Mayer RJ. Ubiquitin and ubiquitin-like proteins as multifunctional signals. Nature reviews Molecular cell biology. 2005;6(8):599–609. pmid:16064136
- 34. Li W, Schmidt W. A lysine‐63‐linked ubiquitin chain‐forming conjugase, UBC13, promotes the developmental responses to iron deficiency in Arabidopsis roots. The Plant Journal. 2010;62(2):330–43. pmid:20113438
- 35. Wen R, Newton L, Li G, Wang H, Xiao W. Arabidopsis thaliana UBC13: implication of error-free DNA damage tolerance and Lys63-linked polyubiquitylation in plants. Plant molecular biology. 2006;61(1–2):241–53. pmid:16786304
- 36. Cui F, Liu L, Zhao Q, Zhang Z, Li Q, Lin B, et al. Arabidopsis ubiquitin conjugase UBC32 is an ERAD component that functions in brassinosteroid-mediated salt stress tolerance. The Plant Cell. 2012;24(1):233–44. pmid:22214659
- 37. Zolman BK, Monroe-Augustus M, Silva ID, Bartel B. Identification and functional characterization of Arabidopsis PEROXIN4 and the interacting protein PEROXIN22. The Plant Cell. 2005;17(12):3422–35. pmid:16272432
- 38. Zhiguo E, Zhang Y, Li T, Wang L, Zhao H. Characterization of the Ubiquitin-Conjugating Enzyme Gene Family in Rice and Evaluation of Expression Profiles under Abiotic Stresses and Hormone Treatments. PloS one. 2015;10(4):e0122621. pmid:25902049
- 39. Proost S, Van Bel M, Sterck L, Billiau K, Van Parys T, Van de Peer Y, et al. PLAZA: a comparative genomics resource to study gene and genome evolution in plants. The Plant Cell. 2009;21(12):3718–31. pmid:20040540
- 40. Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, et al. SMART 4.0: towards genomic data integration. Nucleic acids research. 2004;32(suppl 1):D142–D4.
- 41. Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic acids research. 2003;31(13):3784–8. pmid:12824418
- 42. Guo A, Zhu Q, Chen X, Luo J. [GSDS: a gene structure display server]. Yi chuan = Hereditas/Zhongguo yi chuan xue hui bian ji. 2007;29(8):1023–6.
- 43. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic acids research. 2004;32(5):1792–7. pmid:15034147
- 44. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular biology and evolution. 2013;30(12):2725–9. pmid:24132122
- 45. Papaleo E, Casiraghi N, Arrigoni A, Vanoni M, Coccetti P, De Gioia L. Loop 7 of E2 enzymes: an ancestral conserved functional motif involved in the E2-mediated steps of the ubiquitination cascade. PloS one. 2012;7(7):e40786. pmid:22815819
- 46. Schumacher F-R, Wilson G, Day CL. The N-terminal extension of UBE2E ubiquitin-conjugating enzymes limits chain assembly. Journal of molecular biology. 2013;425(22):4099–111. pmid:23871895
- 47. Winn PJ, Religa TL, Battey JN, Banerjee A, Wade RC. Determinants of functionality in the ubiquitin conjugating enzyme family. Structure. 2004;12(9):1563–74. pmid:15341722
- 48. Fabi JP, Cordenunsi BR, de Mattos Barreto GP, Mercadante AZ, Lajolo FM, Oliveira do Nascimento JR. Papaya fruit ripening: response to ethylene and 1-methylcyclopropene (1-MCP). Journal of agricultural and food chemistry. 2007;55(15):6118–23. pmid:17602654
- 49. Yi Sheng JHH, Doherty Ryan, Srikumar Tharan, Shloush Jonathan, Avvakumov George V., Walker John R., Xue Sheng, Neculai Dante, Wan Janet W., Kim Sung K., Arrowsmith Cheryl H., Raught Brian, Dhe-Paganon Sirano. A Human Ubiquitin Conjugating Enzyme (E2)-HECT E3 Ligase Structure-function Screen. Molecular & Cellular Proteomics Mcp. 2012;11(8):329–41.
- 50. Lin Y-X, Jiang H-Y, Chu Z-X, Tang X-L, Zhu S-W, Cheng B-J. Genome-wide identification, classification and analysis of heat shock transcription factor family in maize. BMC genomics. 2011;12(1):1.
- 51. Peng X, Zhao Y, Cao J, Zhang W, Jiang H, Li X, et al. CCCH-type zinc finger family in maize: genome-wide identification, classification and expression profiling under abscisic acid and drought treatments. PloS one. 2012;7(7):e40120. pmid:22792223
- 52. Wang Y, Wang Q, Zhao Y, Han G, Zhu S. Systematic analysis of maize class III peroxidase gene family reveals a conserved subfamily involved in abiotic stress response. Gene. 2015;566(1):95–108. pmid:25895479
- 53. Zhao Y, Zhou Y, Jiang H, Li X, Gan D, Peng X, et al. Systematic analysis of sequences and expression patterns of drought-responsive members of the HD-Zip gene family in maize. PloS one. 2011;6(12):e28488. pmid:22164299
- 54. Bachmair A, Novatchkova M, Potuschak T, Eisenhaber F. Ubiquitylation in plants: a post-genomic look at a post-translational modification. Trends in plant science. 2001;6(10):463–70. pmid:11590065
- 55. Xu L, Ménard R, Berr A, Fuchs J, Cognat V, Meyer D, et al. The E2 ubiquitin‐conjugating enzymes, AtUBC1 and AtUBC2, play redundant roles and are involved in activation of FLC expression and repression of flowering in Arabidopsis thaliana. The Plant Journal. 2009;57(2):279–88. pmid:18798874
- 56. Jeon EH, Pak JH, Kim MJ, Kim HJ, Shin SH, Lee JH, et al. Ectopic expression of ubiquitin-conjugating enzyme gene from wild rice, OgUBC1, confers resistance against UV-B radiation and Botrytis infection in Arabidopsis thaliana. Biochemical and biophysical research communications. 2012;427(2):309–14. pmid:23000158