Original Research ARTICLE
A Phylogenetically Informed Comparison of GH1 Hydrolases between Arabidopsis and Rice Response to Stressors
- 1State Key Laboratory of Crop Biology, College of Life Science, Shandong Agricultural University, Taian, China
- 2College of Life Sciences, Nantong University, Nantong, China
- 3Ministry of Agricultural Scientific Observing and Experimental Station of Maize in Plain Area of Southern Region, Nantong University, Nantong, China
- 4College of Chemistry, Central China Normal University, Wuhan, China
- 5School of Life Sciences and State Key Laboratory of Agrobiotechnology, The Chinese University of Hong Kong, Shatin, Hong Kong
- 6Shenzhen Research Institute, The Chinese University of Hong Kong, Shenzhen, China
- 7Jiangsu Entry-exit Inspection And Quarantine Bureau, Nanjing, China
Glycoside hydrolases Family 1 (GH1) comprises enzymes that can hydrolyze β-O-glycosidic bond from a carbohydrate moiety. The plant GH1 hydrolases participate in a number of developmental processes and stress responses, including cell wall modification, plant hormone activation or deactivation and herbivore resistance. A large number of members has been observed in this family, suggesting their potential redundant functions in various biological processes. In this study, we have used 304 sequences of plant GH1 hydrolases to study the evolution of this gene family in plant lineage. Gene duplication was found to be a common phenomenon in this gene family. Although many members of GH1 hydrolases showed a high degree of similarity in Arabidopsis and rice, they showed substantial tissue specificity and differential responses to various stress treatments. This differential regulation implies each enzyme may play a distinct role in plants. Furthermore, some of salt-responsive Arabidopsis GH1 hydrolases were selected to test their genetic involvement in salt responses. The knockout mutants of AtBGLU1 and AtBGLU19 were observed to be less-sensitive during NaCl treatment in comparison to the wild type seedlings, indicating their participation in salt stress response. In summary, Arabidopsis and rice GH1 glycoside hydrolases showed distinct features in their evolutionary path, transcriptional regulation and genetic functions.
As sessile organisms, plants encounter a variety of abiotic stresses during their life cycle. For example, cold or freezing stresses are major environmental factors that limit the distribution and productivity of plant species (Du et al., 2010). Low temperature may result in a series of detrimental phenomenon to plants such as pollen sterility and low membrane fluidity. Drought stress leads to water depletion of plant cells, which will destroy the bilayer of plant cell membrane and thereby increase the membrane permeability (Du et al., 2013). High salinity severely hinders crop yield. Its secondary effects, such as osmotic imbalance and ion toxicity, are adverse to plant development (Munns and Tester, 2008). In addition, hyper-osmolarity represents a generalized secondary stress response resulting from cold, drought and high salinity (Yamaguchi-Shinozaki and Shinozaki, 2005).
Glycoside hydrolases or glycosidase are widely distributed enzymes that can hydrolyze the glycosidic bond of carbohydrates (Chandrasekar et al., 2014). To date, 135 families of these enzymes have been classified according to the CAZy database (Xu et al., 2004; Cantarel et al., 2009; Lombard et al., 2014). Animal glycosidases have been shown to play tremendous roles in human disease such as Diabetes and Gaucher disease etc. (Körschen et al., 2013; Antu et al., 2014), indicating their crucial functions in human health. In comparison to animal glucosidases, higher plants contain a large number of genes encoding glycosidases (Xu et al., 2004; Opassiri et al., 2006), suggesting their extensive functional redundancy or differential responsiveness in plants.
In the last decades, although GH1 hydrolases from model plants Arabidopsis, rice and maize have been extensively studied (Xu et al., 2004; Opassiri et al., 2006; Gomez-Anduro et al., 2011; Zhao et al., 2012), their evolutionary origin from plant lineage remains unclear. Furthermore, the substrate specificities of many GH1 hydrolases have been tested (Wang et al., 2011; Luang et al., 2013; Rouyi et al., 2014) and the crystal structures of several representative members have been solved to further validate the molecular mechanisms of their enzymatic cleavage (Seshadri et al., 2009; Roston et al., 2014), which greatly promote the research in this field. However, relatively few functional studies in plants using genetic approach have been reported. Several plant GH1 hydrolases have been reported to play important roles in diverse biological processes including defense against herbivore, cell wall lignification, phytohormone activation and secondary metabolism (Wang et al., 2011; Chapelle et al., 2012; Xu et al., 2012; Zhou et al., 2012; Miyahara et al., 2013). For instance, AtBGLU42 has been reported to modulate iron deficiency response and regulate rhizobacteria-induced systemic resistance in Arabidopsis (Zamioudis et al., 2014). Four members AtBGLU34, AtBGLU35, AtBGLU37, and AtBGLU38, also known as β-S-glucosidases (myrosinases) have been demonstrated to hydrolyze glucosinolate and participate in plant defense responses (Barth and Jander, 2006; Wittstock and Burow, 2010). In addition, two BGLUs, AtBGLU18, and AtBGLU33 have been reported to hydrolyze the glucose ester of abscisic acid (ABA-GE), a storage form of plant hormone ABA (Lee et al., 2006; Xu et al., 2012). Thus, they play important regulatory role in ABA responses. Maize ZmBGLU1 has been characterized to release active cytokinins by hydrolyzing cytokinin-O-glucosides and was involved in defense responses (Kiran et al., 2006).
Given that the complex nature of plant GH1 hydrolases on deglycosylation, our group is thus interested in exploring the molecular mechanisms of this group of enzymes under abiotic stress. By locating the specific substrates (e.g., metabolites or glycosylated proteins), may help us to get a better understanding of their evolutionary divergence and biological significance. Therefore, as an initial step, we constructed a phylogenetic tree using 304 sequences of 14 species in order to reveal the evolution of BGLUs in plant lineages in this study. Subsequently, the spatio-temporal expressions of BGLUs from both Arabidopsis and rice have been studied. Following promoter analysis further demonstrated the transcriptional regulation of Arabidopsis BGLUs is tightly linked to their 5′-flanking regions. At last, the selected Arabidopsis BGLUs which are constitutively induced by salt stress were tested using their T-DNA insertion mutants. These knock-out mutants were less sensitive under sodium chloride (NaCl) treatment, demonstrating the potential roles of AtBGLUs in salt responses. To summarize, our analysis on plant BGLUs using a combination of bioinformatics, molecular and genetic approaches to investigate the evolutionary path, differential transcriptional regulation and potential protein functions of plant BGLUs.
Materials and Methods
Identification of GH1 Hydrolases and Sequence Alignment
Plant protein sequences consisting of a BGLU domain were obtained by searching nonredundant protein sequences from the NCBI database (National Center for Biotechnology Information, http://blast.ncbi.nlm.nih.gov/Blast.cgi) and Phytozome v11.0 (http://www.phytozome.net/). Fourteen available complete genome sequences selected for the phylogenetic analysis include those from model organisms including chlorophytes (Chlamydomonas reinhardtii), bryophyte (Physcomitrella patens), lycophyte (Selaginella moellendorffii), Arabidopsis (Arabidopsis thaliana), grape (Vitis vinifera), cucumber (Cucumis sativus), papaya (Carica papaya), castor bean (Ricinus communis), soybean (Glycine max), barrel medic (Medicago truncatula), sorghum (Sorghum bicolor), rice (Oryza sativa), maize (Zea mays), and cottonwood (Populus trichocarpa). There are 382 protein sequences with BGLU domains were obtained. The similarity amongst the sequences of every BGLUs were aligned using software ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/). The open gap penalty and extend gap penalty have been set to 10 and 0.2 in multiple alignments, respectively. The accession number or gene identifier, number of sequences and sample accession numbers used in this study are listed in the Table S1. Protein sequences (<30% identity) tested by the similarity matrix were removed before the phylogenetic analysis. A number of 304 sequences were selected and used to construct the evolutionary tree across plant lineages. In addition, 189 sequences of selected model plant and crop species were chosen to build the second phylogenetic tree.
Phylogenetic and Protein Structure Analysis
The phylogenetic trees were constructed based on the alignment of plant BGLU conserved domains (Figure S1) and generated using neighbor-joining (NJ) methods implemented in MEGA 7.0. The bootstrap values were calculated from 1,000 replicates.
The structural alignments and superimposition from 304 selected BGLU sequences were performed using the input from ClustalW2 based on crystal structure solved in experiments (i.e., PDB ID, 2JF7, 2JF6). The ConSurf Web (http://bental.tau.ac.il/new_ConSurfDB/) was used for the identification of functional protein regions (Ashkenazy et al., 2010). The conservation scores were calculated by using ML approach in ConSurf. Models and figures were drawn using Pymol (DeLano Scientific).
Multiple Em for Motif Elicitation (MEME)
The MEME analysis has been carried out using protein sequence of Arabidopsis and rice BGLUs. The MEME server (http://meme-suite.org/tools/meme) was used to calculate the conserved motifs within Arabidopsis BGLUs and rice BGLUs, respectively.
Plant Materials and Growth Conditions
The knockout mutant of AtBGLU1 (SALK_060948, GABI_341B12) and AtBGLU19 (SALK_007445, GABI_103F10) were obtained from Arabidopsis Biological Resource Center (ABRC) (https://www.arabidopsis.org/). Seeds used for seedling root length assays were surface-sterilized with 20% bleach and 0.05% Tween-20, and then grown on Murashige and Skoog (MS) medium (Murashige and Skoog, 1962) supplemented with 0.8% w/v agar and 1.5% w/v sucrose under the 16 h light (23°C)/8 h dark (21°C) cycles.
Stress Treatments, RNA Extraction and Quantitative Real-Time PCR
Arabidopsis (Col-0) were grown on MS plates supplemented with 0.8% w/v agar and 1.5% w/v sucrose for 12 days before treatments. The growth condition for rice seedlings (Oryza sativa, Nipponbare) was described previously (Meng et al., 2011). In brief, seeds were germinated in distilled water for 2–4 days at 28°C. Subsequently, 7-d-old seedlings were transferred to a tray and floated in nutrient solution (20 mg/L NH4NO3, 0.05 mg/L KH2PO4, 10 mg/L KCl, 1 mg/L MgCl2, 1 mg/L FeNaEDTA and trace elements in Hoagland solution) in greenhouse under 16 h light /8 h dark (28°C) cycles. 12-day-old Arabidopsis and 14-day-old rice seedlings were treated at 4°C (16 h light/8 h dark cycles) vs. the room temperature control, 15% (w/v) PEG and 125 and 150 mM NaCl vs. water control. Samples (whole seedlings of Arabidopsis and rice) for RNA extraction were collected at different time points after each treatment. Three plates with randomized position of each time point were pooled together. Each treatment was performed twice.
Plant total RNA was extracted by using the RNeasy Mini Kit (Qiagen, Germany) according to the manufacturer's instruction. Total RNA (5 μg) was reverse-transcribed into cDNA by using Superscript First-Strand Synthesis System (Invitrogen, USA) following the manufacturer's protocol. For both the evaluation of tissue-specific expression and transcript changes under stress treatment, real-time quantitative PCR was performed using FastStart Universal SYBR Green Master (Roche, Germany) in a total volume of 15 μl per reaction on StepOne Plus (Applied Biosystems) (Du et al., 2013). Transcript changes were normalized using Arabidopsis and rice ACTIN gene (At3g18780 and Os03g50885) and analyzed from three independent technical replicates as described previously (Czechowski et al., 2005). Each treatment and subsequent RNA extraction and qRT-PCR analysis were performed twice. The gene-specific primer sequences used in this study are shown in Table S2.
Measurement of Primary Root Length
Measurement of primary root length was carried out as described by Du et al. (2013). Four-day-old seedlings of Col-0, atbglu1, and atbglu19 grown on MS plates were transferred to MS control plates and MS medium supplemented with 125 and 150 mM NaCl. The seedlings were grown vertically under light/dark cycles of 16 h/8 h at 23°C. Photos were taken 10 d after transplantation.
Analysis of the 5′-Flanking Regions
The 1.5-kb 5′-flanking regions of Arabidopsis BGLUs was used to identify regulatory motifs. Transcriptional cis-elements were predicted using software server PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (Lescot et al., 2002).
Sequence data in this article are available in the Arabidopsis Genome Initiative and GenBank databases under the accession numbers listed in Figure 1 and Table S1.
Figure 1. Phylogenetic analysis of the β-glucosidases in plant lineage. Phylogenetic analysis of plant β-glucosidases using neighbor-joining (NJ) methods as implemented in MEGA 4.0. Three hundred and four sequences selected from different plant species were used to construct the evolutionary tree. The assigned gene identifiers or accession numbers for the genes are displayed. Bootstrap values are recorded at the nodes. Abbreviations of generic names are listed in the Materials and methods and Table S1.
Identification and Phylogenetic Analysis of Plant GH1 Hydrolases
In order to identify the plant GH1 hydrolases, BLASTp searches coupled with conserved domain analysis were conducted using NCBI database and Phytozome v11.0. A total of 382 GH1 hydrolases sequences from 14 plant species were collected and evaluated in terms of similarity. The 304 sequences with certain degree of similarities (greater than 30% amino acid identity) were then used to construct the phylogenetic tree (Table S1). The selected plant species represented different branches on the evolutionary tree for green plant GH1 hydrolases.
The topologies of phylogenetic tree were complicated due to the large number of candidate sequences (Figure 1). Twelve branches were constructed from the inner arch. The bootstrap value (BS) was low at some clades of the tree, implying that the sequences of this gene family were divergence and less reliable. Therefore, analyzing the whole phylogenetic tree may lose information which is significant in evolutionary path. However, the general trend of evolution, from lower plants to higher plants, was similar to other phylogenetic studies (Meng et al., 2011). For example, the green algae (Ostreococcus lucimarinus) in branch 1 (pointed by black arrow) is a representative member of lineage that diverged before the evolution of land plants. The GH1 hydrolase sequence of this species showed high possibility that was distantly related to the other plant GH1 hydrolases. Surprisingly, AtBGLU42 and OsBGLU40 was grouped in the same clade with the green algal GH1 hydrolase in branch 1. The long distance of these two GH1 hydrolases to the other members in Arabidopsis and rice may indicate they have different ancestors in comparison to the other Arabidopsis or rice GH1 members. Previous phylogenetic study indicates that OsBGLU40 has 70% identity with bacterial proteins (Xu et al., 2004; Opassiri et al., 2006; Gomez-Anduro et al., 2011; Zhao et al., 2012). They proposed that this gene may have been gene transcripts that were captured by retrotransposons and reincorporated into the rice genome, or by lateral gene transfer from a bacteria species. However, no such results mention AtBGLU42. Further analysis is needed to unravel this clustering. Besides, several other sequences from maize and soybean were also found in this branch (Figure 1), indicating that this gene transfer mechanism may be common amongst higher plants. Mosses, together with the lycophytes, are two early-divergent lineages of land plants. The two species selected in this study, Physcomitrella patens and Selaginella moellendorffii, were diversely distributed in several clades (Figure 1). For instance, the branches 5 and 8 were clustered by sequences from Physcomitrella patens, whereas branches 1, 2, 7, 9, and 10 consisted of sequences from Selaginella moellendorffii. They usually formed the basis of each sub-clade together with GH1 hydrolases from other plant species. This phenomenon is similar to the plant kingdom phylogeny, where the protein sequences from these two lower plants built up the basis to the higher plants. However, our analysis indicates that the majority of them were clustered in distinct branches. The GH1 hydrolases from higher plants were not always found to form a monophyletic group. Instead, sequences from the same species were observed to form sub-clades with relatively long evolutionary distance in branches 1, 3, 4, 6, 11, and 12. This phenomenon shown in the phylogenetic analysis may indicate their multi-ancestral origins and diverse functions in plant development and stress responses. Further evidence was obtained using Arabidopsis and rice as a representative dicot and monocot, respectively. The MEME analysis revealed that 2 of 3 motifs identified were not conserved between Arabidopsis BGLUs (Figure S2) and rice BGLUs (Figure S3), suggesting high divergence between the lineages of dicots and monocots.
Tissue Expression Pattern in Arabidopsis BGLUs (AtBGLUs) and Rice BGLUs (OsBGLUs)
Nowadays, transcript abundance is easily obtained from large-scale microarray and transcriptome analysis using online data mining tools or websites. However, the datasets generated from profiling experiments need to be validated using a second method, such as qRT-PCR. Furthermore, both microarray and RNA sequencing have uncertainties to differentiate transcripts with high similarity, which is the case in BGLU gene family. In addition, the steady state mRNA level will be altered due to a number of factors like growth condition, growth medium composition, sampling time and RNA extraction methods etc. Therefore, in order to accurately determine the expression of BGLU genes, primers that can differentiate homologous BGLU genes were designed for the targeted qRT-PCR gene expression analysis (Table S2).
To determine the tissue-specific expression pattern of different GH1 hydrolases, we further focused on Arabidopsis and rice sequences, which are representatives of dicots and monocots, respectively. Although some GH1 BGLUs have been demonstrated to be enzymes other than β-glucosidases (Xu et al., 2004), we examined the expression patterns of annotated 47 Arabidopsis (Col-0) BGLU genes and 37 rice (Nipponbare) BGLU genes in five tissues including seed, flower, leaf, stem and root. As shown in Figure 2, most AtBGLUs and OsBGLUs genes were expressed in different parts of tissues, among which AtBGLU6, AtBGLU22, OsBGLU24, OsBGLU28, OsBGLU29, and OsBGLU34 were highly expressed in roots than in leaves. Meanwhile, a goodly part of BGLU genes such as AtBGLU10, AtBGLU12, AtBGLU13, AtBGLU14, AtBGLU20, OsBGLU7, and OsBGLU26 were found to be highly expressed in flowers. Interestingly, some BGLUs specifically expressed in seeds (e.g., AtBGLU11, AtBGLU18, OsBGLU11, OsBGLU14, and OsBGLU25), suggesting their potential role in seed establishment or germination. The expression of OsBGLU14 was consistent with previous study, where it was highly expressed in embryo and endosperm (Baiya et al., 2014). In addition, the BGLUs classified in the same clades did not always show the co-expression pattern. For example, OsBGLU32 was classified in the same large clade as AtBGLU6. Considerably low expression was detected of OsBGLU32 in root tissue in comparison to the root transcript level of AtBGLU6. Taken together, our results showed that the expression pattern of the both Arabidopsis and rice BGLU gens was diverse and had tissue-specific preference.
Figure 2. Tissue-specific expression of Arabidopsis and rice BGLU genes. Hierarchical cluster analysis was applied to Arabidopsis and rice BGLU genes in different tissue types. The relative gene expression values (log2 scale of qRT-PCR, n = 3, technical replicates) were analyzed using the R language programmed heatmap format. Red color represents the level of expression.
AtBGLUs and OsBGLUs Are Differentially Regulated under Cold Stress
It has been mentioned in several publications that plant GH1 enzymes play an important role in the cold stress response (Fourrier et al., 2008). Therefore, we investigated the gene expression of both Arabidopsis and rice BGLUs in response to a time-course cold treatment. In comparison to rice BGLUs, the Arabidopsis BGLUs generally showed larger variation in their expression under cold treatment in comparison to the control (Figure 3). Some of BGLU genes in Arabidopsis were strongly induced upon the cold stress (Figure 3A). For example, the expression of AtBGLU7 was significantly up-regulated and the expressions of AtBGLU11 and AtBGLU46 were largely down-regulated when the cold treatment lasted for 24 h. Interestingly, very few of OsBGLUs were stimulated by the cold stress even at the different time points of cold treatment (4°C) (Figure 3B). This result suggested that OsBGLU family may not play major roles in cold stress responses in rice or the regulation is not at the transcription level. Since Oryza sativa L. (Nipponbare) used in this study is a tropical species, it maybe not generally adapted to cold stress.
Figure 3. Expression of Arabidopsis and rice BGLU genes under cold treatment. Hierarchical cluster analysis were used for the expression of Arabidopsis and rice BGLU genes after (A) 0.5, 3, 6, and 24 h or (B) 0.5, 3, 6, 12, and 24 h cold treatment. The relative gene expression values (log2 scale of qRT-PCR, n = 3, technical replicates) were analyzed using the R language programmed heatmap format. Red and green colors represent the up-regulation or down-regulation of gene expression, respectively.
Most of AtBGLUs and OsBGLUs Were Regulated under PEG Treatment
Osmotic stress is a representative abiotic stress worldwide. We used 15% PEG treatment to mimic the condition of osmotic stress. It appeared that a total of 16 of Arabidopsis BGLU genes were strongly induced under PEG treatment in comparison to the water control, whereas 11 genes were solely down-regulated (Figure 4A). The rest of the AtBGLUs were either up- or down-regulated in the different time point after PEG treatment. In contrast, most of OsBGLU genes exhibited low expression levels and were not induced or suppressed in early time points (i.e., <3 h) under PEG treatment. At the later time points (i.e., >3 h), the expression of a large number of OsBGLUs were altered (Figure 4B). Therefore, our results showed that the BGLU family genes in Arabidopsis and rice may play distinct roles in drought stress responses.
Figure 4. Expression of Arabidopsis and rice BGLU genes under PEG treatment. Hierarchical cluster analysis were used for the expression of Arabidopsis and rice BGLU genes after (A) 0.5, 3, 6, and 24 h or (B) 0.5, 3, 6, 12, and 24 h 15% (w/v) PEG treatment. The relative gene expression values (log2 scale of qRT-PCR, n = 3, technical replicates) were analyzed using the R language programmed heatmap format. Red and green colors represent the up-regulation or down-regulation of gene expression, respectively.
Gene Expression Patterns of AtBGLUs and OsBGLUs under Salt Treatment Are Similar To the PEG Treatment
High salinity is one of the most important stress conditions in many parts of the world, including northern China. To investigate the gene expression changes of AtBGLUs and OsBGLUs in response to the salt stress, the time course experiment was performed for both Arabidopsis and rice seedlings. Amongst them, large number of AtBGLU family members were up- or down-regulated by NaCl. For example, AtBGLU3, AtBGLU6, AtBGLU21, AtBGLU22 were significantly induced in the 6-h NaCl-treated samples, but was reduced at the 24 h (Figure 5A) in comparison to the water control. On the other hand, Some AtBGLUs like AtBGLU30, AtBGLU45, and AtBGLU46, were down-regulated under NaCl treatment (Figure 5A). The AtBGLU45 and 46 were described as monolignol glucosidases in lignification (Escamilla-Trevino et al., 2006), the relationship between these two enzymes and salt stress remains unclear. Differently, AtBGLU1, AtBGLU11, AtBGLU19, and AtBGLU42 were highly up-regulated upon the NaCl treatment, indicating that they may play crucial roles in response to high salinity (Figure 5A). In contrast, the expression level of OsBGLUs were not largely altered until a 6-h treatment by NaCl (Figure 5B), suggesting that they may participate in long term salt responses.
Figure 5. Expression of Arabidopsis and rice BGLU genes under NaCl treatment. Hierarchical cluster analysis were used for the expression of Arabidopsis and rice BGLU genes after (A) 0.5, 3, 6, and 24 h or (B) 0.5, 3, 6, 12, and 24 h 150 mM NaCl treatment. The relative gene expression values (log2 scale of qRT-PCR, n = 3, technical replicates) were analyzed using the R language programmed heatmap format. Red and green colors represent the up-regulation or down-regulation of gene expression, respectively.
The Arabidopsis T-DNA Mutants of Two AtBGLU Genes Were Less Sensitive to Salt Treatment
The constitutive induction of certain AtBGLUs under salt-treatments draw us attention to perform further analysis. Two knockout T-DNA insertion mutants of AtBGLU1 (SALK_060948, GABI_341B12) and AtBGLU19 (SALK_007445, GABI_103F10) were characterized and subjected to NaCl treatment, respectively (Figures 6A,B). Under 125 mM NaCl treatment, both the atbglu1 and atbglu19 mutants showed less sensitive phenotypes in comparison to the wild-type plants (Figures 6C,D). Longer primary root length was detected after 7-day growth on the NaCl-containing MS plates (Figures 6C,D), indicating that the constitutive induction of AtBGLU1 and AtBGLU19 under NaCl treatment successfully linked the genetic evidence. The potential roles of these two AtBGLUs are being further studied to reveal their molecular mechanisms in response to high salinity.
Figure 6. Phenotypic characterization of AtBGLU1 and AtBGLU19 under NaCl treatment. (A) Schematic view of T-DNA insertion for atbglu1 and atbglu19. Black and gray boxes are exons and 5′ or 3′-UTR regions, respectively. (B) Semi-quantitative reverse transcription polymerase chain reaction (RT-PCR) for identification of Arabidopsis knockout mutants. Plant RNA was extracted from rosette leaves of 4-wk-old WT, atbglu1, and atbglu19 plants. AtACTIN2 was used as control. (C) Relative primary root length amongst WT, atbglu1, and atbglu19 seedlings under NaCl treatment were shown. Values are means ±SE (n = 20; 20 seedlings per genotype were measured from 3 plates and the experiments were repeated twice). (D) WT, atbglu1, and atbglu19 plants seeds were germinated on MS medium for 4 days, and seedlings were subsequently transferred to MS supplemented with 125 mM NaCl for 10 days. Asterisks indicate significant differences compared with wild type by Student-T test (**P < 0.01).
Gene Structures and Genome Locations Imply the Duplication of GH1 Hydrolases in Higher Plants during Their Evolution
The output from BLASTp search for GH1 hydrolases revealed their complexity in gene structures and sequence homology (Figure 1). Eighty amino acid residues show over 90% conservation amongst 304 aligned sequences and were represented in a structure model based on known β-glucosidase structure information (Figure 7). Besides, the other parts of the protein sequences were diverse, even at the place of enzyme active site pockets, suggesting their potential ability to use a large spectrum of substrates. The selected proteins, ranging from 315 to 672 amino acids, were used to construct the evolutionary tree (Figure 1). Although the degree of similarity of some clades was low, some features were observed to reveal their evolutionary relationships. For instance, the presence of multiple homologous members in both early-divergent plant lineages and higher plants demonstrated that the function of GH1 hydrolases evolved before the evolution of land plants (Figure 1). Furthermore, gene duplication was a common feature across the lineages. The coexistence of multifold GH1 hydrolases amongst the plants may indicate their crucial roles in plant survival (Figure 1). Some researchers have studied the phylogenetic relationships within one species or several species belong to the higher plants (Xu et al., 2004; Opassiri et al., 2006; Gomez-Anduro et al., 2011; Zhao et al., 2012). No existing publications have mentioned the evolutionary relationships across the plant lineages (i.e., from lower plants to higher plants), maybe due to the sequence divergence of this gene family. From our results, the divergence of plant GH1 hydrolases started approximately at the stage between the evolution of mosses and lycophytes in comparison to the green algae (Figure 1), suggesting that the evolvement of neo-functionality in this family started as early in lower plants.
Figure 7. Structural alignments from selected BGLU homologs with the canonical BGLU structure. Ribbon representations of the BGLU structure are colored according to conservation in 304 selected BGLU sequences from phylogenetic analysis, the color representation blue to red indicates the increasing of conservation.
Promoter Analysis of AtBGLUs Correlates Their Expression Data in Response to Abiotic Stresses
Regulatory promoters have been considered to be one of the most important factors in regulating gene expression (Chen et al., 2014b). They often located in the few hundred base pairs before the transcription start site of a gene. To further explain the relationships between the gene expression and genetic phenotypes in this study, the 1.5-kb 5′-flanking regions of all the AtBGLUs were used to analyze the presence of stress-responsive cis-elements. The detailed results were listed in Table 1. Multiple cis-elements in response to temperature and drought stresses have been detected by software PlantCARE (Lescot et al., 2002), verifying that the regulation of the gene expression largely correlates to the motifs found in the promoter regions of AtBGLUs. For instance, LTR motifs have been reported to be regarded as responsive cis-element under low temperature (Hughes and Dunn, 1996). The expression of AtBGLU9 was up-regulated at 6 and 24 h after the cold treatment. The three LTR motifs located in the 5′-flanking region of this gene served as good candidates to this regulation (Table 1). Interestingly, the AtBGLU19, which was identified to be involved in salt stress response, had several ABA-responsive elements in their promoter regions such as G-box and ABA-responsive elements (ABRE) (Table 1). Abscisic acid constitutes one of the major pathways in response to salt stress (Chinnusamy et al., 2004; Hirayama and Shinozaki, 2010). Transcriptional regulation of ABA-responsive genes converges to transcription factors. To date, various transcription factors from different families in Arabidopsis have been reported to function in ABA-mediated gene activation or repression (Fujita et al., 2013). Both G-box and ABRE are candidate sequences for the binding of ABA-RESPONSIVE ELEMENT BINDING PROTEINS/ABA-RESPONSIVE ELEMENT BINDING FACTORS (AREB/ABFs) (Choi et al., 2000; Uno et al., 2000). This group of transcription factors belongs to the core module of ABA signal transduction. They belong to the basic leucine zipper (bZIP) transcription factor family (Jakoby et al., 2002) and use G-box type ABRE with a consensus ACGT core. Two independent groups isolated several members of this family through yeast-one hybrid screening (Choi et al., 2000; Uno et al., 2000). From phylogenetic analysis, all the ABA-responsive bZIP transcription factors belong to Group A bZIPs (Jakoby et al., 2002). Several transcription factors are well documented including AREB1/AtABF2, AREB2/AtABF4, ABF3, and ABI5 amongst the 13 members in this group (Brocard et al., 2002; Carles et al., 2002; Kim et al., 2004; Fujita et al., 2005; Yoshida et al., 2010). The high density of this type of motifs found in the promoters of AtBGLU19 indicates the possibility that the induction of these two genes under NaCl treatment results from the ABA-responsive core regulation. In contrast, only one G-box was found in the promoter of AtBGLU1 (Table 1), suggesting that other unknown motifs may participate in salt responses to drive its expression.
Table 1. List of stress-responsive elements located in the 1.5-kb 5′-flanking regions of Arabidopsis BGLUs.
Besides G-box and ABREs, some other cis-elements like P-box (Itoh et al., 1995), TATC-box (Gubler et al., 1995) and GARE (Sutliff et al., 1993) for gibberellins signal transduction, ERE for ethylene signaling pathways (Itzhaki and Woodson, 1993), HSE for heat-shock responses (Gurley and Key, 1991) and ARE for anaerobic induction (Walker et al., 1987) have been observed, indicating that multiple stress responses can alter the gene expression of AtBGLUs.
Importance of Plant GH1 Hydrolases in Plant Engineering to Tolerate Abiotic Stresses
In comparison to the biotic stress, abiotic stress showed more difficulties in the application of plant engineering (Wang et al., 2003; Bressan et al., 2009). Transformation of a single gene to target plants is difficult to enhance their abiotic resistance (Chen et al., 2014a). However, abiotic stresses such as low temperature, low water supply and high salinity, share some common secondary stress conditions like the adjustment of osmotic potential (Hirayama and Shinozaki, 2010; Mittler and Blumwald, 2010). Given these phenomenon, engineering effector proteins involved in cross-talk pathways may be beneficial.
Plant GH1 hydrolases are considered to be an ideal target for plant engineering (McManus and Osborne, 1991) since they are involved in different developmental processes and various stress responses by deglycosylation of targeted metabolites or proteins (Xu et al., 2004; Opassiri et al., 2006). In this study, experimental and bioinformatic analysis indicates that both rice and Arabidopsis BGLUs were involved in multiple abiotic stress responses. Thus, using GH1 hydrolases as engineering materials, we can achieve various goals such as manipulation of active metabolite homeostasis, enzyme subcellular localization or activity etc. However, the relatively few functional studies have been reported may due to several lines of difficulties in molecular biology and plant genetics. For example, the number of genes in this family is large and some of the members are highly conserved, giving high degree of redundancy. In addition, the spectrum of substrates is diverse and most of the potential substrates are not commercially available. Thus, the works to determine the substrate specificity for each GH1 hydrolase are time-consuming and technically difficult. Furthermore, some GH1 hydrolases locate in the vicinity of each other in plant chromosome, bringing difficulties for plant double mutant generation by conventional crossing. Therefore, the development of high-throughput technology to screen the substrates for GH1 hydrolases is necessary (Wei et al., 2015). The progress in plant mutant generation using CRISPR-Cas9 technologies (Miao et al., 2013) have served a possible solution on the generation of double or multiple mutants for the members of this gene family.
In conclusion, the evolution and expression of plant GH1 hydrolases have been studied in this research article. The phylogenetic analysis revealed gene duplication may be one of the most important factors in the evolution of GH1 hydrolases across plant lineages. Although some family members of AtBGLUs and OsBGLUs share high similarity, the tissue and stress-responsive expression patterns are diverse. Therefore, their distinct roles in plant development and stress responses need further investigation.
Y-YC, JZ, and Y-GL designed research scheme. Y-YC, J-FY, T-YL, ZF, F-YZ, Z-FS, TF, and L-JW performed experiments. M-XC, N-HY, and G-FH analyzed data. Y-YC, JZ, Y-GL, and M-XC wrote and revised the manuscript.
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.
This work was supported by University Natural Science Research Project of Jiangsu Province (13KJB210005), Science and technology program of Nantong (MS12016044), National Natural Science Foundation of China (NSFC31101099), Natural Science Foundation of Shandong Province (BS2015NY002), Science and technology program of Shandong Colleges and Universities (J10LC01) and Jiangsu Overseas Research and Training Program for University Prominent Young and Middle-aged Teachers and Presidents. We thank Ms. Zhang Y (The Chinese University of Hong Kong) for language editing.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00350/full#supplementary-material
Antu, K. A., Riya, M. P., Mishra, A., Anilkumar, K. S., Chandrakanth, C. K., Tamrakar, A. K., et al. (2014). Antidiabetic property of Symplocos cochinchinensis is mediated by inhibition of alpha glucosidase and enhanced insulin sensitivity. PLoS ONE 9:e105829. doi: 10.1371/journal.pone.0105829
Ashkenazy, H., Erez, E., Martz, E., Pupko, T., and Ben-Tal, N. (2010). ConSurf 2010, calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res. 38, W529–W533. doi: 10.1093/nar/gkq399
Baiya, S., Hua, Y., Ekkhara, W., and Ketudat Cairns, J. R. (2014). Expression and enzymatic properties of rice (Oryza sativa L.) monolignol beta-glucosidases. Plant Sci. 227, 101–109. doi: 10.1016/j.plantsci.2014.07.009
Barth, C., and Jander, G. (2006). Arabidopsis myrosinases TGG1 and TGG2 have redundant function in glucosinolate breakdown and insect defense. Plant J. 46, 549–562. doi: 10.1111/j.1365-313X.2006.02716.x
Brocard, I. M., Lynch, T. J., and Finkelstein, R. R. (2002). Regulation and role of the Arabidopsis abscisic acid-insensitive 5 gene in abscisic acid, sugar, and stress response. Plant Physiol. 129, 1533–1543. doi: 10.1104/pp.005793
Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The Carbohydrate-Active EnZymes database (CAZy), an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663
Carles, C., Bies-Etheve, N., Aspart, L., Leon-Kloosterziel, K. M., Koornneef, M., Echeverria, M., et al. (2002). Regulation of Arabidopsis thaliana Em genes, role of ABI5. Plant J. 30, 373–383. doi: 10.1046/j.1365-313X.2002.01295.x
Chandrasekar, B., Colby, T., Emran Khan Emon, A., Jiang, J., Hong, T. N., Villamor, J. G., et al. (2014). Broad-range glycosidase activity profiling. Mol. Cell. Proteomics 13, 2787–2800. doi: 10.1074/mcp.O114.041616
Chapelle, A., Morreel, K., Vanholme, R., Le-Bris, P., Morin, H., Lapierre, C., et al. (2012). Impact of the absence of stem-specific beta-glucosidases on lignin and monolignols. Plant Physiol. 160, 1204–1217. doi: 10.1104/pp.112.203364
Chen, M. X., Zheng, S. X., Yang, Y. N., Xu, C., Liu, J. S., Yang, W. D., et al. (2014b). Strong seed-specific protein expression from the Vigna radiata storage protein 8S Galpha promoter in transgenic Arabidopsis seeds. J. Biotechnol. 174, 49–56. doi: 10.1016/j.jbiotec.2014.01.027
Chinnusamy, V., Schumaker, K., and Zhu, J. K. (2004). Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J. Exp. Bot. 55, 225–236. doi: 10.1093/jxb/erh005
Czechowski, T., Stitt, M., Altmann, T., Udvardi, M. K., and Scheible, W.-R. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 139, 5–17. doi: 10.1104/pp.105.063743
Du, Z. Y., Chen, M. X., Chen, Q. F., Xiao, S., and Chye, M. L. (2013). Overexpression of Arabidopsis acyl-CoA-binding protein ACBP2 enhances drought tolerance. Plant Cell Environ. 36, 300–314. doi: 10.1111/j.1365-3040.2012.02574.x
Du, Z. Y., Xiao, S., Chen, Q. F., and Chye, M. L. (2010). Depletion of the membrane-associated acyl-coenzyme A-binding protein ACBP1 enhances the ability of cold acclimation in Arabidopsis. Plant Physiol. 152, 1585–1597. doi: 10.1104/pp.109.147066
Escamilla-Trevino, L. L., Chen, W., Card, M. L., Shih, M. C., Cheng, C. L., and Poulton, J. E. (2006). Arabidopsis thaliana beta-Glucosidases BGLU45 and BGLU46 hydrolyse monolignol glucosides. Phytochemistry 67, 1651–1660. doi: 10.1016/j.phytochem.2006.05.022
Fourrier, N., Bedard, J., Lopez-Juez, E., Barbrook, A., Bowyer, J., Jarvis, P., et al. (2008). A role for SENSITIVE TO FREEZING2 in protecting chloroplasts against freeze-induced damage in Arabidopsis. Plant J. 55, 734–745. doi: 10.1111/j.1365-313X.2008.03549.x
Fujita, Y., Fujita, M., Satoh, R., Maruyama, K., Parvez, M. M., Seki, M., et al. (2005). AREB1 is a transcription activator of novel ABRE-dependent ABA signaling that enhances drought stress tolerance in Arabidopsis. Plant Cell 17, 3470–3488. doi: 10.1105/tpc.105.035659
Fujita, Y., Yoshida, T., and Yamaguchi-Shinozaki, K. (2013). Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants. Physiol. Plant. 147, 15–27. doi: 10.1111/j.1399-3054.2012.01635.x
Gomez-Anduro, G., Ceniceros-Ojeda, E. A., Casados-Vazquez, L. E., Bencivenni, C., Sierra-Beltran, A., Murillo-Amador, B., et al. (2011). Genome-wide analysis of the beta-glucosidase gene family in maize (Zea mays L. var B73). Plant Mol. Biol. 77, 159–183. doi: 10.1007/s11103-011-9800-2
Gubler, F., Kalla, R., Roberts, J. K., and Jacobsen, J. V. (1995). Gibberellin-regulated expression of a myb gene in barley aleurone cells, evidence for Myb transactivation of a high-pI alpha-amylase gene promoter. Plant Cell 7, 1879–1891.
Itoh, K., Yamaguchi, J., Huang, N., Rodriguez, R. L., Akazawa, T., and Shimamoto, K. (1995). developmental and hormonal regulation of rice [alpha]-amylase(RAmy1A)-gusA fusion genes in transgenic rice seeds. Plant Physiol. 107, 25–31. doi: 10.1104/pp.107.1.25
Jakoby, M., Weisshaar, B., Dröge-Laser, W., Vicente-Carbajosa, J., Tiedemann, J., Kroj, T., et al. (2002). bZIP transcription factors in Arabidopsis. Trends Plant Sci. 7, 106–111. doi: 10.1016/S1360-1385(01)02223-3
Kim, S., Kang, J. Y., Cho, D. I., Park, J. H., and Kim, S. Y. (2004). ABF2, an ABRE-binding bZIP factor, is an essential component of glucose signaling and its overexpression affects multiple stress tolerance. Plant J. 40, 75–87. doi: 10.1111/j.1365-313X.2004.02192.x
Kiran, N. S., Polanska, L., Fohlerova, R., Mazura, P., Valkova, M., Smeral, M., et al. (2006). Ectopic over-expression of the maize beta-glucosidase Zm-p60.1 perturbs cytokinin homeostasis in transgenic tobacco. J. Exp. Bot. 57, 985–996. doi: 10.1093/jxb/erj084
Körschen, H. G., Yildiz, Y., Raju, D. N., Schonauer, S., Bönigk, W., Jansen, V., et al. (2013). The non-lysosomal β-glucosidase GBA2 is a non-integral membrane-associated protein at the endoplasmic reticulum (ER) and Golgi. J. Biol. Chem. 288, 3381–3393. doi: 10.1074/jbc.M112.414714
Lee, K. H., Piao, H. L., Kim, H. Y., Choi, S. M., Jiang, F., Hartung, W., et al. (2006). Activation of glucosidase via stress-induced polymerization rapidly increases active pools of abscisic acid. Cell 126, 1109–1120. doi: 10.1016/j.cell.2006.07.034
Lescot, M., Dehais, P., Thijs, G., Marchal, K., Moreau, Y., Van de Peer, Y., et al. (2002). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi: 10.1093/nar/30.1.325
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42, D490–495. doi: 10.1093/nar/gkt1178
Luang, S., Cho, J. I., Mahong, B., Opassiri, R., Akiyama, T., Phasai, K., et al. (2013). Rice Os9BGlu31 is a transglucosidase with the capacity to equilibrate phenylpropanoid, flavonoid, and phytohormone glycoconjugates. J. Biol. Chem. 288, 10111–10123. doi: 10.1074/jbc.M112.423533
McManus, M. T., and Osborne, D. J. (1991). Identification and characterization of an ionically-bound cell-wall glycoprotein expressed preferentially in the leaf rachis abscission zone of Sambucus-nigra L. J. Plant Physiol. 138, 63–67. doi: 10.1016/S0176-1617(11)80731-5
Meng, W., Su, Y. C., Saunders, R. M., and Chye, M. L. (2011). The rice acyl-CoA-binding protein gene family, phylogeny, expression and functional analysis. New Phytol. 189, 1170–1184. doi: 10.1111/j.1469-8137.2010.03546.x
Miyahara, T., Sakiyama, R., Ozeki, Y., and Sasaki, N. (2013). Acyl-glucose-dependent glucosyltransferase catalyzes the final step of anthocyanin formation in Arabidopsis. J. Plant Physiol. 170, 619–624. doi: 10.1016/j.jplph.2012.12.001
Opassiri, R., Pomthong, B., Onkoksoong, T., Akiyama, T., Esen, A., and Ketudat Cairns, J. R. (2006). Analysis of rice glycosyl hydrolase family 1 and expression of Os4bglu12 beta-glucosidase. BMC Plant Biol. 6:33. doi: 10.1186/1471-2229-6-33
Roston, R. L., Wang, K., Kuhn, L. A., and Benning, C. (2014). Structural determinants allowing transferase activity in SENSITIVE TO FREEZING 2, classified as a family I glycosyl hydrolase. J. Biol. Chem. 289, 26089–26106. doi: 10.1074/jbc.M114.576694
Rouyi, C., Baiya, S., Lee, S. K., Mahong, B., Jeon, J. S., Ketudat-Cairns, J. R., et al. (2014). Recombinant expression and characterization of the cytoplasmic rice beta-glucosidase Os1BGlu4. PLoS ONE 9:e96712. doi: 10.1371/journal.pone.0096712
Seshadri, S., Akiyama, T., Opassiri, R., Kuaprasert, B., and Cairns, J. K. (2009). Structural and enzymatic characterization of Os3BGlu6, a rice β-glucosidase hydrolyzing hydrophobic glycosides and (1, 3)-and (1, 2)-linked disaccharides. Plant Physiol. 151, 47–58. doi: 10.1104/pp.109.139436
Sutliff, T. D., Lanahan, M. B., and Ho, T. H. (1993). Gibberellin treatment stimulates nuclear factor binding to the gibberellin response complex in a barley alpha-amylase promoter. Plant Cell 5, 1681–1692. doi: 10.1105/tpc.5.11.1681
Uno, Y., Furihata, T., Abe, H., Yoshida, R., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2000). Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc. Natl. Acad. Sci. U.S.A. 97, 11632–11637. doi: 10.1073/pnas.190309197
Walker, J. C., Howard, E. A., Dennis, E. S., and Peacock, W. J. (1987). DNA sequences required for anaerobic expression of the maize alcohol dehydrogenase 1 gene. Proc. Natl. Acad. Sci. U.S.A. 84, 6624–6628. doi: 10.1073/pnas.84.19.6624
Wang, P., Liu, H., Hua, H., Wang, L., and Song, C. P. (2011). A vacuole localized β-glucosidase contributes to drought tolerance in Arabidopsis. Chinese Sci. Bull. 56, 3538–3546. doi: 10.1007/s11434-011-4802-7
Wang, W., Vinocur, B., and Altman, A. (2003). Plant responses to drought, salinity and extreme temperatures, towards genetic engineering for stress tolerance. Planta 218, 1–14. doi: 10.1007/s00425-003-1105-5
Wei, H., Brunecky, R., Donohoe, B. S., Ding, S. Y., Ciesielski, P. N., Yang, S., et al. (2015). Identifying the ionically bound cell wall and intracellular glycoside hydrolases in late growth stage Arabidopsis stems, implications for the genetic engineering of bioenergy crops. Front. Plant Sci. 6:315. doi: 10.3389/fpls.2015.00315
Xu, Z., Escamilla-Trevino, L., Zeng, L., Lalgondar, M., Bevan, D., Winkel, B., et al. (2004). Functional genomic analysis of Arabidopsis thaliana glycoside hydrolase family 1. Plant Mol. Biol. 55, 343–367. doi: 10.1007/s11103-004-0790-1
Xu, Z.-Y., Lee, K. H., Dong, T., Jeong, J. C., Jin, J. B., Kanno, Y., et al. (2012). A vacuolar β-glucosidase homolog that possesses glucose-conjugated abscisic acid hydrolyzing activity plays an important role in osmotic stress responses in Arabidopsis. Plant Cell 24, 2184–2199. doi: 10.1105/tpc.112.095935
Yamaguchi-Shinozaki, K., and Shinozaki, K. (2005). Organization of cis-acting regulatory elements in osmotic-and cold-stress-responsive promoters. Trends Plant Sci. 10, 88–94. doi: 10.1016/j.tplants.2004.12.012
Yoshida, T., Fujita, Y., Sayama, H., Kidokoro, S., Maruyama, K., Mizoi, J., et al. (2010). AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. Plant J. 61, 672–685. doi: 10.1111/j.1365-313X.2009.04092.x
Zamioudis, C., Hanson, J., and Pieterse, C. M. (2014). β-Glucosidase BGLU42 is a MYB72-dependent key regulator of rhizobacteria-induced systemic resistance and modulates iron deficiency responses in Arabidopsis roots. New Phytol. 204, 368–379. doi: 10.1111/nph.12980
Keywords: β-glucosidase, gene expression, multiple alignment, promoter, phylogenetic analysis, stress treatments
Citation: Cao Y-Y, Yang J-F, Liu T-Y, Su Z-F, Zhu F-Y, Chen M-X, Fan T, Ye N-H, Feng Z, Wang L-J, Hao G-F, Zhang J and Liu Y-G (2017) A Phylogenetically Informed Comparison of GH1 Hydrolases between Arabidopsis and Rice Response to Stressors. Front. Plant Sci. 8:350. doi: 10.3389/fpls.2017.00350
Received: 29 September 2016; Accepted: 28 February 2017;
Published: 24 March 2017.
Edited by:Nigel G. Halford, Rothamsted Research (BBSRC), UK
Reviewed by:Neil J. Willey, University of the West of England, UK
Gang Li, University of Nebraska-Lincoln, USA
Copyright © 2017 Cao, Yang, Liu, Su, Zhu, Chen, Fan, Ye, Feng, Wang, Hao, Zhang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work.