Effective genome editing and identification of a regiospecific gallic acid 4-O-glycosyltransferase in pomegranate (Punica granatum L.)

Abstract

Pomegranate (Punica granatum L.) trees are woody perennials that bear colorful and nutritious fruits rich in phenolic metabolites, e.g., hydrolyzable tannins (HTs) and flavonoids. We here report genome editing and gene discovery in pomegranate hairy roots using Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) (CRISPR/Cas9), coupled with transcriptome and biochemical analyses. Single guide RNAs (sgRNAs) were designed to target two UDP-dependent glycosyltransferases (UGTs), PgUGT84A23 and PgUGT84A24, which possess overlapping activities in β-glucogallin (a galloylglucose ester; biosynthetic precursor of HTs) biosynthesis. A unique accumulation of gallic acid 3- O- and 4- O-glucosides (galloylglucose ethers) was observed in the PgUGT84A23 and PgUGT84A24 dual CRISPR/Cas9-edited lines (i.e., ugt84a23 ugt84a24) but not the control (empty vector) or PgUGT84A23/ PgUGT84A24 single edited lines (ugt84a23 or ugt84a24). Transcriptome and real-time qPCR analyses identified 11 UGTs with increased expression in the ugt84a23 ugt84a24 hairy roots compared to the controls. Of the 11 candidate UGTs, only PgUGT72BD1 used gallic acid as substrate and produced a regiospecific product gallic acid 4- O-glucoside. This work demonstrates that the CRISPR/Cas9 method can facilitate functional genomics studies in pomegranate and shows promise for capitalizing on the metabolic potential of pomegranate for germplasm improvement.

References

[1] Wu, S. & Tian, L. Diverse phytochemicals and bioactivities in the ancient fruit and modern functional food pomegranate (Punica granatum). Molecules 22, 1606 (2017).

[2] Bar-Ya’akov, I., Tian, L., Amir, R. & Holland, D. Primary metabolites, anthocyanins, and hydrolyzable tannins in the pomegranate fruit. Front. Plant Sci. 10, 620 (2019).

[3] Holland, D., Hatib, K. & Bar-Ya’akov, I. in Horticultural Reviews (ed. Janick, J.) 127-191 (John Wiley & Sons, 2009).

[4] Wang, H., Russa, M. L. & Qi, L. S. CRISPR/Cas9 in genome editing and beyond. Annu. Rev. Biochem. 85, 227-264 (2016).

[5] Ma, X., Zhu, Q., Chen, Y. & Liu, Y.-G. CRISPR/Cas9 platforms for genome editing in plants: developments and applications. Mol. Plant 9, 961-974 (2016).

[6] Koonin, E. V., Makarova, K. S. & Zhang, F. Diversity, classification and evolution of CRISPR-Cas systems. Curr. Opin. Microbiol. 37, 67-78 (2017).

[7] Knott, G. J. & Doudna, J. A. CRISPR-Cas guides the future of genetic engineering. Science 361, 866-869 (2018).

[8] Chen, K., Wang, Y., Zhang, R., Zhang, H. & Gao, C. CRISPR/Cas genome editing and precision plant breeding in agriculture. Annu. Rev. Plant Biol. 70, 667-697 (2019).

[9] Ono, N., Bandaranayake, P. C. G. & Tian, L. Establishment of pomegranate (Punica granatum) hairy root cultures for genetic interrogation of the hydrolyzable tannin biosynthetic pathway. Planta 236, 931-941 (2012).

[10] Ono, N., Qin, X., Wilson, A., Li, G. & Tian, L. Two UGT84 family glycosyltransferases catalyze a critical reaction of hydrolyzable tannin biosynthesis in pomegranate (Punica granatum). PLoS ONE 11, e0156319 (2016).

[11] Lim, E. et al. The activity of Arabidopsis glycosyltransferases toward salicylic acid, 4-hydroxybenzoic acid, and other benzoates. J. Biol. Chem. 277, 586-592 (2002).

[12] Schuster, B. & Herrmann, K. Hydroxybenzoic and hydroxycinnamic acid derivatives in soft fruits. Phytochemistry 24, 2761-2764 (1985).

[13] Pawlowska, A. M., De Leo, M. & Braca, A. Phenolics of Arbutus unedo L. (Ericaceae) fruits: identification of anthocyanins and gallic acid derivatives. J. Agric. Food Chem. 54, 10234-10238 (2006).

[14] Lu, Y. & Foo, L. Y. The polyphenol constituents of grape pomace. Food Chem. 65, 1-8 (1999).

[15] Lim, E.-K. et al. Evolution of substrate recognition across a multigene family of glycosyltransferases in Arabidopsis. Glycobiology 13, 139-145 (2003).

[16] Qin, G. et al. The pomegranate (Punica granatum L.) genome and the genomics of punicalagin biosynthesis. Plant J. 91, 1108-1128 (2017).

[17] Yuan, Z. et al. The pomegranate (Punica granatum L.) genome provides insights into fruit quality and ovule developmental biology. Plant Biotechnol. J. 16, 1363-1374 (2018).

[18] Park, J., Bae, S. & Kim, J.-S. Cas-Designer: a web-based tool for choice of CRISPR-Cas9 target sites. Bioinformatics 31, 4014-4016 (2015).

[19] Zuker, M . Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 31, 3406-3415 (2003).

[20] Mao, Y. et al. Application of the CRISPR-Cas system for efficient genome engineering in plants. Mol. Plant 6, 2008-2011 (2013).

[21] Clark, J . Cetyltrimethyl ammonium bromide (CTAB) DNA miniprep for plant DNA isolation. Cold Spring Harb. Protoc. 2009, pdb.prot5177 (2009).

[22] Wilson, A. E., Wu, S. & Tian, L. PgUGT95B2 preferentially metabolizes flavones/flavonols and has evolved independently from flavone/flavonol UGTs identified in Arabidopsis thaliana. Phytochemistry 157, 184-193 (2019).

[23] Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[24] Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011).

[25] Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550- 550 (2014).

[26] Bradford, M. A rapid and sensitive for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248-254 (1976).

[27] Wilson, A. E. et al. Characterization of a UGT84 family glycosyltransferase provides new insights into substrate binding and reactivity of galloylglucose ester-forming UGTs. Biochemistry 56, 6389-6400 (2017).

[28] Edgar, R. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792-1797 (2004).

[29] Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870-1874 (2016).