Transcriptomic analyses identify albino-associated genes of a novel albino tea germplasm ‘Huabai 1’

Abstract

Albinism in shoots of tea plants is a common phenotypic expression which gives the tea infusion a pleasant umami taste. A novel natural albino mutant tea germplasm containing high amino acids content was found and named as ‘Huabai 1’. ‘Huabai 1’ has white jade tender shoots under low temperature and turns green with increased temperature. In order to understand the molecular mechanism of color change in leaf of ‘Huabai 1’, transcriptome analysis was performed to identify albino-associated differentially expressed genes (DEGs). A total of 483 DEGs were identified from white shoots of ‘Huabai 1’ compared to its green shoots. There were 15 DEGs identified to be involved in phenylpropanoid biosynthesis, which account for the majority of characterized DEGs. The metabolites related to phenylpropanoid biosynthesis revealed similar expression pattern of DEGs. Furthermore, metabolic pathways such as ubiquonone, tyrosine, and flavonoid biosynthesis associated with phenylpropanoid biosynthesis could also contribute to the color change in ‘Huabai 1’ tender shoots. Protein–protein interaction analysis revealed a hub protein NEDD8 (CSA009575) which interacted with many regulated genes in spliceosome, nitrogen metabolism, phenylpropanoid biosynthesis, and other pathways. In conclusion, the findings in this study indicate that the color change of ‘Huabai 1’ tender shoots is a combined effect of phenylpropanoid biosynthesis pathway and other metabolic pathways including flavonoid biosynthesis in tea plants. Chlorophyll biosynthesis-related genes LHCII and SGR may also play some roles in color change of ‘Huabai 1’.

References

[1] Wu, Z. et al. A chlorophyll-deficient rice mutant with impaired chlorophyllide esterification in chlorophyll biosynthesis. Plant Physiol. 145, 29-40 (2007).

[2] Campbell, B. W. et al. Identical substitutions in magnesium chelatase paralogs result in chlorophyll-deficient soybean mutants. G3 5, 123-131 (2015).

[3] Zhu, L. et al. Genetic characterisation and fine mapping of a chlorophyll-deficient mutant (BnaC.ygl) in Brassica napus. Mol. Breed. 34, 603-614 (2014).

[4] Choi, H. G. et al. Yield loss and quality degradation of strawberry fruits cultivated under the deficient insolation conditions by shading. Hortic. Environ. Biotech. 55, 263-270 (2014).

[5] Slattery, R. A. et al. Photosynthesis, light use efficiency, and yield of reduced-chlorophyll soybean mutants in field conditions. Front. Plant Sci. 8, 549 (2017).

[6] Feng, L. et al. Determination of quality constituents in the young leaves of albino tea cultivars. Food Chem. 155, 98-104 (2014).

[7] Du, Y. Y. et al. A study on the chemical composition of albino tea cultivars. J. Hortic. Sci. Biotech. 81, 809-812 (2006).

[8] Wei, K. et al. Comparison of catechins and purine alkaloids in albino and normal green tea cultivars (Camellia sinensis L.) by HPLC. Food Chem. 130, 720-724 (2012).

[9] Yuan, H., Zhang, J., Nageswaran, D. & Li, L. Carotenoid metabolism and regulation in horticultural crops. Hortic. Res. 2, 15036 (2015).

[10] Wang, M., Wang, G. & Ji, J. Suppression of the phytoene desaturase gene influence on the organization and function of photosystem II (PSII) and antioxidant enzyme activities in tobacco. Environ. Exp. Bot. 67, 460-466 (2010).

[11] Li, C.-F. et al. Differential Metabolic Profiles during the Albescent Stages of ‘Anji Baicha’ (Camellia sinensis). PLoS ONE 10, e0139996 (2015).

[12] Li, Q. et al. Proteomic analysis of young leaves at three developmental stages in an albino tea cultivar. Proteome Sci. 9, 44 (2011).

[13] Song, L. et al. Molecular link between leaf coloration and gene expression of flavonoid and carotenoid biosynthesis in Camellia sinensis cultivar ‘Huangjinya’. Front. Plant Sci. 8, 803 (2017).

[14] Xia, E.-H. et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol. Plant 10, 866-877 (2017).

[15] Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[16] Trapnell, C. et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511-515 (2010).

[17] Deng, Y. Y. et al. Integrated nr database in protein annotation system and Its localization. Computer Engineering 32, 71-73 (2006).

[18] Finn, R. D. et al. Pfam: the protein families database. Nucleic Acids Res. 42, D222-D230 (2014).

[19] Tatusov, R. L., Galperin, M. Y., Natale, D. A. & Koonin, E. V. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 28, 33-36 (2000).

[20] Koonin, E. V. et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 5, R7- R7 (2004).

[21] Apweiler, R. et al. UniProt: the Universal Protein knowledgebase. Nucleic Acids Res. 32, D115-D119 (2004).

[22] Kanehisa, M. et al. The KEGG resource for deciphering the genome. Nucleic Acids Res. 32, D277-D280 (2004).

[23] Ashburner, M. et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25, 25-29 (2000).

[24] Florea, L., Song, L. & Salzberg, S. L. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. F1000Res. 2, 188 (2013).

[25] Schulze, S. K. et al. SERE: single-parameter quality control and sample comparison for RNA-Seq. BMC Genomics 13, 524 (2012).

[26] Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol. 11, R106 (2010).

[27] Ma, Q. P., Hao, S., Chen, X. & Li, X. H. Validation of reliability for reference genes under various abiotic stresses in tea plant. Russ. J. Plant Physiol. 63, 423-432 (2016).

[28] Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402-408 (2001).

[29] Li, C.-F. et al. Biochemical and transcriptomic analyses reveal different metabolite biosynthesis profiles among three color and developmental stages in ‘Anji Baicha’ (Camellia sinensis). BMC Plant Biol. 16, 195 (2016).

[30] Vogt, T. Phenylpropanoid biosynthesis. Mol. Plant 3, 2-20 (2010).

[31] Boerjan, W., Ralph, J. & Baucher, M. Lignin Biosynthesis. Annu. Rev. Plant. Biol. 54, 519-546 (2003).

[32] Bouvier d’Yvoire, M. et al. Disrupting the cinnamyl alcohol dehydrogenase 1 gene (BdCAD1) leads to altered lignification and improved saccharification in Brachypodium distachyon. Plant J. 73, 496-508 (2013).

[33] Zhao, Q. et al. Loss of function of cinnamyl alcohol dehydrogenase 1 leads to unconventional lignin and a temperature-sensitive growth defect in Medicago truncatula. PNAS 110, 13660-13665 (2013).

[34] Takeda, Y. et al. Regulation of coniferaldehyde 5-hydroxylase expression to modulate cell wall lignin structure in rice. Planta 246, 337-349 (2017).

[35] Tian, X. et al. Sinapyl alcohol derivatives from the lipo-soluble part of Dichrocephala benthamii C. B. Clarke. Molecules 18, 1720 (2013).

[36] Pascal, A. A. et al. Molecular basis of photoprotection and control of photosynthetic light-harvesting. Nature 436, 134 (2005).

[37] Sakuraba, Y. et al. STAY-GREEN and chlorophyll catabolic enzymes interact at light-harvesting complex II for chlorophyll detoxification during leaf senescence in arabidopsis. Plant Cell 24, 507-518 (2012).

[38] Schwechheimer, C. & Mergner, J. The NEDD8 modification pathway in plants. Front. Plant Sci. 5, 103 (2014).

[39] Boh, B. K., Smith, P. G. & Hagen, T. Neddylation-induced conformational control regulates cullin RING ligase activity in vivo. J. Mol. Biol. 409, 136-145 (2011).

[40] Chen, S. et al. Quantitative proteomics analysis reveals the tolerance of Mirabilis jalapa L. to petroleum contamination. Environ. Sci. Pollut. Res. 24, 7375-7382 (2017).