Identification, Conserved Domain and Sequence Feature Analysis of CsLAZYs
A total of six LAZY genes were identified in Arabidopsis thaliana [20]. Subsequently, the six AtLAZYs genes were used as queries through Basic Local Alignment Search Tool (BLAST) analysis against the tea plant genome (http://tpia.teaplant.org/Blast.html) [22]. Initially, a total of 15 candidate unique genes were obtained in tea plant, and multiple sequence alignments of all LAZY genes were performed among tea plant, Arabidopsis and rice (data not shown). As a result, only 3 unique genes contained regions of conserved sequence V that possessing an ethylene-responsive element-binding factor-associated amphiphilic repression (EAR) motif (LxLxL) (Figure 1A), and it is an indispensable conserved domain of LAZY [20, 23]. Thereafter, the obtained three genes are referred to as CsLAZY1 (CSS025254), CsLAZY2 (CSS049138) and CsLAZY3 (CSS020288), which were located in different scaffolds (Table 1). Their amino acid lengths are 399 aa (CsLAZY1), 367 aa (CsLAZY2) and 251 aa (CsLAZY3), respectively. Furthermore, the molecular weights (Mw) of CsLAZY1 to CsLAZY3 were 44.2, 41.2 and 29.0, respectively, and their isoelectric points (pI) were 6.55, 6.18 and 6.47, respectively (Table 1).
Subsequently, the exon/intron organization in the coding sequences of each CsLAZY gene was performed (Figure 1B), which contained 5 (CsLAZY1), 5 (CsLAZY2), and 4 (CsLAZY3) exons, respectively. In terms of intron and exon length, CsLAZY1 is the longest while CsLAZY3 is the shortest. The coding sequence of CsLAZY1 was cloned and sequenced, demonstrating that the cloned cDNA was totally consistent with the genomic reference sequence.
Evolution and Phylogenetic Analysis of LAZY Genes
Previous studies and numerous fully sequenced plant genomes make it possible to perform a comparative genomic analysis of LAZY genes across a broad range of plant species. The LAZY genes have been identified in many species that play a similar role, so we performed iterative BLAST searches to determine the phylogeny of LAZY1 genes. The LAZY1 genes were identified from 20 distinct plant species, and homology analysis of LAZY1 among algal, lowland species, monocots, and dicots provided further insight into the evolutionary processes of this gene family. Phylogenetic analyses showed these LAZY1 genes were highly conserved among algae, monocots and dicots, and LAZY1 gene evolved from primitive organisms despite their overall relatively low sequence similarities (Figure 2A). Furthermore, we analyzed them in evolutionary trees with other species that including Populus tomentosa, Vitis vinifera, Solanaceae lycopersicon, Oryza sativa, Arabidopsis thaliana, and Camellia sinensis. The complete LAZY gene families including 23 members from six plant species were used for phylogenetic analysis, displaying that the three CsLAZY genes were distinctly classified into two clades of class I and class II (Figure 2B). CsLAZY1 and CsLAZY2 were grouped into class I categories, and they both have high sequence similarity with the protein sequences of VvLAZY1, PtLAZY1 and PtLAZY2. CsLAZY3 was grouped into class II and showed the highest sequence similarity with VvLAZY3 gene.
Analysis of cis-elements in promoters and tissue expression patterns
To explore the potential difference in non-coding regions of CsLAZYs, a 2-kb flanking sequence upstream of the translation start codon was obtained, and several putative cis-regulatory elements in the promoter were identified using the PLACE and PlantCARE databases (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). The light-sensitive cis-elements comprised the largest part of all elements, including the Box 4, TCT-motif, ATC-motif, ATCT-motif, G-box, I-box, chs-CMA1a, MRE, and ACE (Figure 3A). The hormone-sensitive cis-elements have CGTCA-motif, GARE-motif, TCA-element, ERE, ABRE, TATC-box, and TGACG-motif, while only one type of cis-element contained in the promoter of each gene. It is noteworthy that the promoter region of CsLAZY1 contained many hormone cis-elements, including MeJA, GA, SA and ABA hormone responsive elements, implying that CsLAZY1 may plays an important role in tea plant by responds to hormones.
To understand the potential role of CsLAZYs in tea plant, we downloaded RNA-Seq data of eight tissues from the tea plant genome database, displaying that the expression levels of the three CsLAZY genes were obvious specificity in tissues (Figure 3B). For instance, CsLAZY1 was mainly expressed in stem, followed by in bud and leaf, while basically not expressed in fruit and root. In comparison, CsLAZY2 was mainly expressed in flower, and CsLAZY3 showed the highest expression level in the second leaf. Due to stem bending is one of the main causes of the branch angle, therefore CsLAZY1 was probably played a vital role in regulating branch angle of tea plant.
Expression pattern of CsLAZY1 in tissues among different tea varieties
To further verify the tissue expression pattern of CsLAZY1, we examined the tissue expression level of CsLAZY1 in different varieties. A total of eight varieties were selected, which possessed different branch angles, including four open type varieties (Benshan, Foshou, Yaoshanxiulv, and Tieguanyin) and four erect type varieties (Echa 5, Fuzao 2, Longjingchangye, and Zhenghedabaicha). The expression level of CsLAZY1 transcript varied significantly among the four tissues (leaf, bud, root and stem) (Figure 4). It was showed that CsLAZY1 transcript was not detected in root of the eight tea varieties, and it had the highest expression level in stem, followed by in leaf. Unexpectedly, CsLAZY1 had the highest expression level in leaf, followed by in stem in the cultivar of Tieguanyin.
Subcellular localization of CsLAZY1 protein
To obtain insight into the molecular function of CsLAZY1 protein, we constructed CsLAZY1-GFP and Pk7WGF2 35S-GFP fusion protein expression vectors to examine its subcellular localization. Transient expression in Arabidopsis protoplasts showed that CsLAZY1 protein was localized in the plasma membrane (Figure 5A). Besides, we transferred the CsLAZY1 protein to Agrobacterium to infect tobacco leaves, and obtained the identical result, namely CsLAZY1 protein was localized in the plasma membrane (Figure 5B).
Over-expression of CsLAZY1 in Arabidopsis
To further investigate the role of CsLAZY1 in shoot gravitropism, we transferred CsLAZY1 into Arabidopsis thaliana. The expression of CsLAZY1 was detected using real-time PCR assay in over-expressed (OE) plants but not detected in wild-type (WT) plants, and the three OE lines were named as OELAZY1-11, OELAZY1-20, and OELAZY1-24, respectively (Figure 6A). Subsequently, gravitropism assays through time-lapse imaging were implemented to survey the response to reorientation of the WT and three OE lines. All the seedlings with main stem of 5-10 cm were treated in 90° inverted gravity processing. In light, images were collected by computer-controlled cameras at 0, 30, 60, 90 and 120 minutes of inversion (Figure 6B), and the angle of the hypocotyls were measured from the images. Obviously, the OE plants bended upward slightly at 30 minutes, while no bend was seen from the WT plants. After 90 minutes of inversion, OE plants were reached the maximum bending angle, while the WT bended upward slightly at 90 minutes (Figure 6B). Compared to the WT plants, significant difference of bending angle was observed from OE plants after 30 minutes treatment (Figure 6C).
In dark, images were collected at 0, 30, 60, 90, 120, 150 and 180 minutes of inversion, and the angle of the hypocotyls were measured from the images. Both the WT and OE plants in dark bended upward later than their corresponding plants in light, demonstrating that the OE plants bended upward slightly after 60 minutes, while the WT plants bended upward slightly after 90 minutes (Figure 7A). In comparison, significant difference of bending angle was observed from OE plants after 90 minutes (Figure 7B). Consistent with the expression patterns, evidence indicate that CsLAZY1 may play a vital role in respond to gravitropism in stem of tea plant.