Identification of PbrPUB gene family members
In our study, we used a strictly pipeline to identify PUB genes in pear genome. First, the hidden Markov model (HMM) of the U-box domain (PF04564) was downloaded from the Pfam30.0 database, and used as a query to identify the candidate PUB members in Chinese white pear genomic database (http://peargenome.njau.edu.cn/) usingHMMER3.0. As a result, a total of 91 candidate PUB genes were identified in pear genome. Second, SMART tools were performed to verify the accuracy of 91 candidate PUB gens, and 29 PUB genes were removed because of the incomplete of PUB domain. At last, 62 PUB genes with complete U-box domain were obtained for further analysis (Table 1). We named these PbrPUB genes from PbrPUB1 to PbrPUB62 according to their location information on the chromosome. The ExPASy ProtParam (https://www.expasy.org/), an online proteomics and sequence analysis tool, was used to estimate molecular weight and isoelectric point (pI). The molecular weight for the PbrPUB gene family range from 39.33 kDa to 151.30 kDa (Kilodalton) and the pI range from 4.99 to 8.83, with an average of 6.78. Subcellular localization were also predicted by Cell-PLoc 2.0 (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/), and we found that most PUB gene were located in nucleus, except six pub located in cytoplasm and three located in cell membrane.
Phylogenetic analysis of PbrPUB gene family members
To investigate the evolutionary relationships of PbrPUB genes family, we constructed a neighbor joining(NJ) phylogenetic tree using the Mega-X program based on the full-length PUB proteins of pear (62 members), tomato (62 members) and Arabidopsis (61 members) (Fig. 1a). The protein sequences of PUB genes of tomato and Arabidopsis were obtained from previous study [10, 12]. Based the result of phylogenetic tree, 185 members of PUB genes from these three species were clustered into five groups, including Group I, Group II, Group III, Group IV and Group V. The member number of Group III was biggest in five subgroups, and it harboured 64 PUB genes. However, Group IV harboured least PUB genes, with 10 PUB genes. In generally, the PUB genes of pear and tomato were clustered into one subclade, suggesting that pear and tomato exhibited relatively closer relationship compared to Arabidopsis.
It is interesting to note that the member number of PUB gene family in these three species is similar. This result indicated that the number of PUB genes in these three species is conserved. To explore which group of pears had occurred expansion or lost during evolution process, we measure the number of PUB genes of each species in each group. In pear, group I, II, III, IV and V contain 11, 21, 21, 3 and 6 PbrPUB gene family members, respectively. In tomato, groups I, II, III, IV and V contain 12, 21, 21, 3, and 5 SIU-box genes, respectively. In Arabidopsis, groups I, II, III, IV and V contain 20, 12, 22, 4 and 3 AtPUB genes, respectively (Fig. 1b). The member number of each group in pear and tomato is almost equal, suggesting that pear had not undergone expansion or lost compared to tomato. However, compared to pear and tomato, the group I of Arabidopsis had undergone rapid expansion, while the group II of Arabidopsis had undergone rapid lost.
Analysis of PbrPUB gene family conserved motifs and gene structures
To further verify the classification results of phylogenetic tree, we investigated the conserved motif and gene structure of PbrPUB genes in pear. The MultipleEm for Motif Elicitation (MEME) motif search tool were used to predict the conserved motif of PbrPUB genes. A total of 20 motifs were estimated in our study, and we named as motif 1–20 (Fig. 2a, 2b, Additional file 1: Fig. S1). Among them, motif 1, 3 and 5 were found in all groups, indicating that were highly conserved in all PbrPUB proteins. Based on the SMATR website, we determined that the U-box was comprised of Motif1, Motif3 and Motif5 (Additional file 2: Fig.S2). This result provided evidence to support the accuracy of PUB gene set identified in our study. Based on the SMATR website, we also found the other conserved domain: ARM and Pkinase domin. The ARM is comprised of motif 2, 4 and 7; the Pkinase is comprised of motif 11, 13 and 20.Generally, most PbrPUB members in the same groups had similar motifs. For example, most of the members of group II contained motif 6, 10 and 8. This result indicated that these three motifs might be key functional domain of group II PUB genes, suggesting that these proteins might have conservative functions.
To explore the gene structure of PbrPUB genes in pear, we extracted the exon-intron information of 62 PbrPUB genes from pear database using in-house scripts. Based on the information, TBtools software were preformed to show the gene structure of PbrPUB (Fig. 2c). The number of exon in PUB genes was greatly divergent, ranging from 1 to 20. Among 62 PUB genes in pear, PbrPUB24 contained the greatest number of exon (20), while 16 PUB genes only contained one exon. Furthermore, the lengths of the exon and intron were differential. There are 30 PUB genes have been found contain UTR structure. Similarly, to the result of motif analysis, the PUB genes with similar gene structure were cluster into same subclade. For example, most members of class II only housed one exon. This result indicated that the members of same groups exhibited similar gene structure and conserved motifs. These results from conserved motifs and gene structure analysis provided strong evidence to support the accuracy of the classification result of phylogenetic tree.
Chromosomal localization and homologous gene analysis of PbrPUB genes
To further explore the distribution pattern of PUB gene in pear genome, the PbrPUB chromosomal distribution map was plotted using TBtools (Fig. 3a). The location information of PUB genes in pear were extracted by our in-house scripts. As a result, a total of 50 PbrPUB genes (82.26%) were unevenly mapped on the 17 pear chromosomes, and no member of PbrPUB gene family was mapped on chromosome 8. Therefore, we didn’t shown chromosome 8 in our Fig. 3. In addition, 12 genes were located on scaffold contigs, and we also didn’t show them in our Fig. 3. Chromosome 15 had the largest number of PbrPUB genes, with eight PbrPUB genes, followed by chromosome 5 with 6 genes. Chromosome 1, 2 and 12 each contained 4 PbrPUB genes. Two or three PbrPUB genes were mapped on chromosomes 3, 6, 7, 9, 10, 11, 13, 14, and 16. Chromosome 4 and chromosome 17 contained only one gene. We also identified the homologous genes of PUB gene family using MCscanX software. As result, 16 homologous gene pairs were identified in pear PUB gene family, which contained 26 homologous genes. Three homologous gene pairs were detected between chromosomes 5 and chromosomes 10 (Fig. 3b).
Cis-acting elements analysis in the putative promoter of PbrPUB genes
Cis-acting elements were important clues for the prediction of gene functions. To further investigate the function of PbrPUB genes, we predicted the cis-acting element of the putative promoter region of PbrPUB genes using PlantCARE databse. In this study, the 2000 bp upstream region of transcription initiation site of PbrPUB genes was defined as the putative promoter region. As a result, a tot al of 41 cis-acting elements were identified (Fig. 4), and we selected 15 interesting cis-acting elements for further analysis. These 15 cis-acting were associated with stress, hormone, plant growth and development. As shown in Fig. 4a, some diverse distribution patterns of cis-acting elements were observed in the promoter region of PbrPUB genes, indicating that the PUB gene family of pear particular in various different biology process. Meanwhile, we found that all PbrPUB genes contained the cis-acting related to hormone regulation, such as, abscisic acid (ABA), methyl jasmonate (MeJA), gibberellin (GA), salicylic acid, and auxin responsiveness elements. Previous study had reported that DSG1, which encodes a U-box domain, could regulate cell division and elongation by responding to multiple hormones, such as auxin, salicylic acid and ethylene [33]. In our study, 55 genes were identified as the responsiveness elements of ABA, suggesting that PUB gene family might particular in resistance under ABA treatment (Fig. 4b). In Arabidopsis, AtPUB9, AtPUB18, AtPUB19, and AtPUB44 were identified to involve in ABA response [34]. It is notable that the element related to MYB binding site involved in drought was predicted in 44 PbrPUB genes, suggesting that these 44 PbrPUB genes might mediated by MYB genes response to drought stress. Moreover, there were 30 PbrPUB genes have cis-acting elements related to cold, suggesting that these 30 PbrPUB genes might particular in resistance under low temp treatment. As we all know, flavonoid biosynthesis is one of important phenomenon during the process of response to stress in plant. In the study, we found PbrPUB10, PbrPUB24 and PbrPUB5 contained MYB binding site involved in flavonoid biosynthetic.
Expression profile of PbrPUB genes in different tissues of pear
To further explore the tissues-specific expression of PbrPUB genes, we conducted RNA-seq analysis using transcriptomic data from seven different tissues of ‘Dangshansuli’ cultivar pear, including stem, ovary, petal, sepal, bud, fruit and leaf [35]. We used RPKM (reads per kilobase per million) values to estimate the expression level of PbrPUB genes. The RPKM value equal to 0 means no expression in one library, and the RPKM value greater than 10 means highly expression in one library. Then, we investigated the expression level of 62 PUB genes. Pheatmap, an R package, was used to show the expression patterns of 62 PbrPUB genes (Fig. 5a). Based on the expression pattern of 62 PUB genes, they were clustered into four main classes. Gene in Class IV exhibited highly expression level in all of seven examined tissues, while gene in Class II exhibited almost no expression in all of seven tissues. Class I was specifically expressed in pear leaf, and a diversity of expression pattern were detected in Class III (Fig. 5a). Among the 62 PbrPUB genes, 52 genes (83.87%) were at least expressed in one tissue, even though the transcript abundance of several genes was relatively lower for certain tissues. Approximately 10 non-expressed PUB genes were identified in all of seven tissues, and they may be lost the function during the evolution process of PUB genes family in pear. 29 PbrPUB genes were expressed in all seven different tissues, indicating that they have various roles in the development of different tissues. Interesting, we found 28 PUB genes exhibited highest expression in leaf, suggesting that these 28 genes might involve the development of leaf. Due to leaf is an important plant organ involved resistance, we referred that these 28 PUB genes might particular in resistance in the process of pear growth and development.
To verify the transcriptome sequences analysis was reliable, 15 PbrPUB genes were randomly selected to conduct a quantitative real-time PCR (qRT-PCR) experiment to investigate the expression levels in seven different tissues of the ‘Dangshansuli’ pear, including stem, leaf, petal, fruit, sepal, bud, and ovary (Fig. 5b). We found that all of 15 PbrPUB genes exhibited a diversity of expression patterns in the seven different tissues, suggesting that PbrPUB genes family may function in different tissues and participate in diverse metabolic processes. Seven genes (PbrPUB1, PbrPUB3, PbrPUB7, PbrPUB9, PbrPUB18, PbrPUB36 and PbrPUB38) exhibited a similar expression pattern with a high expression level in leaf tissues, suggesting that PbrPUB genes play critical functions during leaf development. All of these seven genes exhibited highly expression level in leaf in transcriptome data. These result provided further evidence to support our transcriptome sequences analysis was reliable. Interestingly, most of 15 PUB genes were highly expressed in reproductive organ, suggesting that PbrPUB genes might associate with the development of reproductive organ.
The expression pattern of PbrPUB genes under abiotic stresses
Previous study had extensively reported PbrPUB gene family involved in various abiotic stresses [36]. To explore the function of PUB gene family in pear, Pyrus betulaefolia from grown 45-day-old seedlings were subjected to four different stress treatments including dehydration, low temperature, ABA and salt. To identify whether PUB genes in pear could respond to abiotic stresses, 11 PbrPUB genes were randomly selected from each group to conduct a qRT-PCR experiment to detect their expression level. 11 genes are comprise of 2 from group I (PbrPUB1 and PbrPUB14), 4 from group II (PbrPUB12, PbrPUB18, PbrPUB36 and PbrPUB38), 2 from group III (PbrPUB3 and PbrPUB25), 2 from group IV (PbrPUB7 and PbrPUB48) and 1 from group V (PbrPUB34).
To identify candidate PbrPUB genes related to dehydration treatment, the shoots of seedlings of pear were placed on dry filter papers for 0, 1, 6, 9, 12 and 24 h at ambient environment (Fig. 6a). Among the eleven PUB genes, 9 PUB genes were up-regulated expressed and one PUB gene (PbrPUB7) was down-regulated expressed under dehydration stress. However, PbrPUB14 was not significantly differential expressed under dehydration stress. Among the 9 up-regulated genes, PbrPUB18 exhibited highly increased expression level during the process of dehydration treatment, while PbrPUB12, PbrPUB3 and PbrPUB36 were up-regulated expressed during 12 h dehydration treatment and recovered to normal levels at 24 h. PbrPUB1, PbrPUB38 and PbrPUB25, exhibited highest expression level at 1 h, where PbrPUB12, PbrPUB14, PbrPUB3 and PbrPUB36 exhibited highest expression level at 12 h under dehydration treatment. These results suggested that PbrPUB1, PbrPUB38 and PbrPUB25 respond to dehydration treatment faster than that of PbrPUB12, PbrPUB14, PbrPUB3 and PbrPUB36. Therefore, PUB gene family in pear play vital role in the process of low temperature stress response.
To verify whether PUB gene in pear involve in low temperature stress, the seedlings were placed in the chamber set at 4°C for 0, 1, 6, 9, 12, 24, 48 and 96 h (Fig. 6b). We detected the expression level of those 11 PbrPUB genes by qRT-PCR experiment. 4 genes (PbrPUB12, PbrPUB3, PbrPUB36 and PbrPUB48) were up-regulated expressed under cold stress, suggesting that those PbrPUB genes might respond to low temperature. PbrPUB12, PbrPUB48 and PbrPUB36 were highly increased during the 48 h low temperature exposure. The expression level of PbrPUB3 was reached to double peak at 1 h and 48 h.
To explore the functions of PbrPUB genes under the salt stress, the seedlings were placed in solution containing 200 mM NaCl solution for 0, 2, 4, 6, 8, 12 and 36 h (Fig. 6c). We detected the expression level of those 11 PbrPUB genes by qRT-PCR experiment. On the whole, all of the selected 11 PbrPUB genes were significantly up-regulated expressed under the 200 mM salt stress treatment. The expression level of PbrPUB14, PbrPUB25, PbrPUB3, PbrPUB48 and PbrPUB7 were highly increased during the 12 h salt exposure. Moreover, PbrPUB1, PbrPUB12, PbrPUB18, PbrPUB34, PbrPUB36 and PbrPUB38 were highest expressed at 4 h under salt stress, suggesting that these 6 PbrPUB genes respond to salt treatment actively. We focus on the expression level of PbrPUB18. In the 4–8 h, the expression level of PbrPUB18 was significantly increased, and then it was down-regulated at 12 h, finally recovered normal level at 36 h.
Previous study had reported that PUB gene involved in ABA-mediated drought stress responses. To investigate the roles of PUB gene family in ABA stress, the seedlings were dipped in solution containing 100 µM ABA for 0, 1, 3, 6, 9, 12 and 36 h (Fig. 6d). As results, all of 11 PUB genes were respond to the ABA stress, and these gene were unregulated expressed at first, and then were down-regulated at 36 h after ABA treatment. These results indicated that PUB genes play important roles in ABA-regulated pathway. The expression levels of three genes (PbrPUB1, PbrPUB25, and PbrPUB36) were reached to peak at 1 h, suggesting that these three genes were actively responded to ABA stress. Interestingly, we found that PbrPUB18 was expressed in 6 h and 12 h after ABA treatment.
Subcellular localization of PbrPUB18
To further explore the function of PUB gene in pear, we constructed subcellular localization experiment to identify where the PUB genes act function. To further verify the biologic function of PbrPUB genes in pear under drought stress, PbrPUB18 was selected from 9 up-regulated expressed genes for further study. The full-length ORF of PbrPUB18 was fused to the N-terminal of GFP (Green Fluorescent Protein), under the control of CaMV 35S promoter to form a fusion construct 35S-PbrPUB18-GFP. Then the 35S-PbrPUB18-GFP fusion protein and 35S-GFP (as a control) was transient transformed into tobacco leaves and the fluorescence signal was observed with a confocal laser scanning microscope. The green fluorescence of GFP control was found in the membrane and the nucleus (Fig. 7a). In contrast, 35S-PbrPUB18-GFPfused GFP protein was only existed in the nucleus and integrated perfectly with DAPI (4′, 6-diamidino-2-phenylindole) regime (Fig. 7b), suggesting that PbrPUB18 was located in the nucleus, which was consistent with our prediction in Table 1.
Assessment of drought tolerance in transgenic lines of PbrPUB18
To further confirm the biologic function of PbrPUB18 gene under drought stress, Arabidopsis Col-0 plants (WT) were transformed by the floral dip method [37]. Two overexpression lines OE-4 and OE-5 were screened out by PCR identification at DNA level in T0 generation plants. The selected robust Arabidopsis thaliana seedlings were transplanted into the soil as T1 generation plants, and the T1 generation positive seedlings were identified again by semi-quantitative PCR at mRNA level. QRT-PCR also verified the expression of PbrPUB18 in OE-4 and OE-5 far above in WT (Additional file 3 Fig S3). To assess the function of overexpression PbrPUB18 in Arabidopsis on drought tolerance, 15-day-old WT and transgenic lines were used to expose to drought stress by withholding water for 12 days. There was no morphological difference between WTs and the transgenic lines in the normal condition. After 12 days without water, the two transgenic lines showed more tolerance to the drought stress, as manifested by lesser leaf-wilting symptoms compared with the WT plants (Fig. 8a). In addition, chlorophyll fluorescence measurements were recorded to further verify drought tolerance of WTs and the transgenic lines (Fig. 8b). The maximum quantum efficiency of the photochemistry (Fv/Fm) values was not affected by species and growth conditions, but under stress conditions, this parameter decreased significantly. After 12 days drought treatment, the Fv/Fm values of WT was significantly lower than of the two transgenic lines, suggesting WT showed more sensitivity to the drought stress (Fig. 8e). Electrolyte leakage (EL), an important indicator of cell injury, was measured after drought stress. The EL of two transgenic lines were only approximate 15 %–20 % compared to WT (37.3 %), suggesting that WT suffered more severe membrane damage than transgenic lines of Arabidopsis by overexpressing PbrPUB18 (Fig. 8c). The transgenic plants displayed significantly lower malondialdehyde (MDA) contents than WT exposure to drought condition (Fig. 8d).
In the drought tolerance stress assay, we found that two transgenic lines had lower values of EL and MDA, implying that they might be subjected to lighter extent oxidative stress than the wild type. Histochemical staining with 3, 3′-diaminobenzidine (DAB) and nitro-blue tetrazolium chloride (NBT) was used to analyze the in situ accumulation of H2O2 and O2-, two main reactive oxygen species (ROS), respectively. After drought stress the leaves of WT lines (Fig. 8f) were stained to a more serious extent compared with transgenic Arabidopsis, implying that more ROS was produced in the WT under the drought conditions. Similar to staining results, quantitative measurements further demonstrated that H2O2 contents in the two transgenic lines were remarkably lower than those of WT (Fig. 8g). And anti- O2-contents in the two transgenic lines were remarkably more than those of WT (Fig. 8h), implying that less ROS was produced in the two transgenic lines under the drought conditions.