Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway

Abstract

Ginger (Zingiber officinale), the type species of Zingiberaceae, is one of the most widespread medicinal plants and spices. Here, we report a high-quality, chromosome-scale reference genome of ginger ‘Zhugen’, a traditionally cultivated ginger in Southwest China used as a fresh vegetable, assembled from PacBio long reads, Illumina short reads, and high-throughput chromosome conformation capture (Hi-C) reads. The ginger genome was phased into two haplotypes, haplotype 1 (1.53 Gb with a contig N50 of 4.68 M) and haplotype 0 (1.51 Gb with a contig N50 of 5.28 M). Homologous ginger chromosomes maintained excellent gene pair collinearity. In 17,226 pairs of allelic genes, 11.9% exhibited differential expression between alleles. Based on the results of ginger genome sequencing, transcriptome analysis, and metabolomic analysis, we proposed a backbone biosynthetic pathway of gingerol analogs, which consists of 12 enzymatic gene families, PAL, C4H, 4CL, CST, C3’H, C3OMT, CCOMT, CSE, PKS, AOR, DHN, and DHT. These analyses also identified the likely transcription factor networks that regulate the synthesis of gingerol analogs. Overall, this study serves as an excellent resource for further research on ginger biology and breeding, lays a foundation for a better understanding of ginger evolution, and presents an intact biosynthetic pathway for species-specific gingerol biosynthesis.

References

[1] Ravindran, P. N., Babu, K. N., Ravindran, P. N. & Babu, K. N. (eds) Ginger the genus Zingiber (CRC press, 2005).

[2] Li, H. et al. Ginger for health care: an overview of systematic reviews. Complement Ther. Med. 45, 114-23 (2019).

[3] Ahmad, B. et al. A Review on pharmacological properties of zingerone (4-(4-Hydroxy-3-methoxyphenyl)-2-butanone). Sci. World J. 2015, 816364 (2015).

[4] Mao, Q. Q. et al. Bioactive compounds and bioactivities of Ginger (Zingiber officinale Roscoe). Foods 8, 185 (2019).

[5] Khan, M., Ullah, S., Azhar, M., Komal, W., Muhammad, W. & Karbaschi, P. A mini-review on the therapeutic potential of Zingiber officinale (ginger). J. Nat. Prod. 15, 1 (2019).

[6] Semwal, R. B., Semwal, D. K., Combrinck, S. & Viljoen, A. M. Gingerols and shogaols: important nutraceutical principles from ginger. Phytochemistry 117, 554-68 (2015).

[7] Tanaka, Y., Sasaki, N. & Ohmiya, A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 54, 733-49 (2008).

[8] Bode, A. M., Ma, W. Y., Surh, Y. J. & Dong, Z. Inhibition of epidermal growth factor-induced cell transformation and activator protein 1 activation by [6]-gingerol. Cancer Res. 61, 850-3 (2001).

[9] Lee, H. S., Seo, E. Y., Kang, N. E. & Kim, W. K. [6]-Gingerol inhibits metastasis of MDA-MB-231 human breast cancer cells. J. Nutr. Biochem. 19, 313-9 (2008).

[10] Kress, W. J., Prince, L. M. & Williams, K. J. The phylogeny and a new classification of the gingers (Zingiberaceae): evidence from molecular data. Am. J. Bot. 89, 1682-96 (2002).

[11] Ramachandran, K. Chromosome numbers in zingiberaceae. Cytologia 34, 213-21 (1969).

[12] Ramachandran, K. Polyploidy induced in ginger by colchicine treatment. Curr. Sci. 51, 288-9 (1982).

[13] Adaniya, S. & Shoda, M. Meiotic irregularity in ginger (Zingiber officinale Roscoe). Chromosome Sci. 2, 141-4 (1998).

[14] Jha, N. N. et al. Effect of curcumin analogs on α-synuclein aggregation and cytotoxicity. Sci. Rep. 6, 28511 (2016).

[15] Paulraj, F., Abas, F., H Lajis, N., Othman, I. & Naidu, R. Molecular pathways modulated by curcumin analogue, diarylpentanoids in cancer. Biomolecules 9, 270 (2019).

[16] Ramirez-Ahumada Mdel, C., Timmermann, B. N. & Gang, D. R. Biosynthesis of curcuminoids and gingerols in turmeric (Curcuma longa) and ginger (Zingiber officinale): identification of curcuminoid synthase and hydroxycinnamoyl-CoA thioesterases. Phytochemistry 67, 2017-29 (2006).

[17] Katsuyama, Y., Miyazono, K., Tanokura, M., Ohnishi, Y. & Horinouchi, S. Structural and biochemical elucidation of mechanism for decarboxylative condensation of beta-keto acid by curcumin synthase. J. Biol. Chem. 286, 6659-68 (2011).

[18] Rodrigues, J. L., Prather, K. L., Kluskens, L. D. & Rodrigues, L. R. Heterologous production of curcuminoids. Microbiol. Mol. Biol. Rev. 79, 39-60 (2015).

[19] Wang, J. P., Liu, B., Sun, Y., Chiang, V. L. & Sederoff, R. R. Enzyme-enzyme interactions in monolignol biosynthesis. Front. Plant Sci. 9, 1942 (2018).

[20] Hassaninasab, A., Hashimoto, Y., Tomita-Yokotani, K. & Kobayashi, M. Discovery of the curcumin metabolic pathway involving a unique enzyme in an intestinal microorganism. Proc. Natl Acad. Sci. USA 108, 6615-20 (2011).

[21] Park, S. B. et al. Structural and biochemical characterization of the curcumin-reducing activity of CurA from vibrio vulnificus. J. Agric. Food Chem. 66, 10608-16 (2018).

[22] Alix, K., Gérard, P. R., Schwarzacher, T. & Heslop-Harrison, J. Polyploidy and interspecific hybridization: partners for adaptation, speciation and evolution in plants. Ann. Bot. 120, 183-94 (2017).

[23] Kyriakidou, M., Tai, H. H., Anglin, N. L., Ellis, D. & Strömvik, M. V. Current strategies of polyploid plant genome sequence assembly. Front. Plant Sci. 9, 1660 (2018).

[24] Yang, J. et al. Haplotype-resolved sweet potato genome traces back its hexaploidization history. Nat. Plants 3, 696-703 (2017).

[25] Zhang, J. et al. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nat. Genet. 50, 1565-73 (2018).

[26] Chen, H. et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat. Commun. 11, 2494 (2020).

[27] Wei, C. et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl Acad. Sci. USA 115, E4151-E4158 (2018).

[28] Chen, F. et al. Genome sequences of horticultural plants: past, present, and future. Hortic. Res. 6, 112 (2019).

[29] Hasing, T., Tang, H., Brym, M., Khazi, F., Huang, T. & Chambers, A. A phased Vanilla planifolia genome enables genetic improvement of flavour and production. Nat. Food 1, 811 (2020).

[30] Bonawitz, N. D. & Chapple, C. The genetics of lignin biosynthesis: connecting genotype to phenotype. Annu. Rev. Genet. 44, 337-63 (2010).

[31] Yang, F. et al. A maize gene regulatory network for phenolic Metabolism. Mol. Plant 10, 498-515 (2017).

[32] Sun, B. et al. Purple foliage coloration in tea (Camellia sinensis L.) arises from activation of the R2R3-MYB transcription factor CsAN1. Sci. Rep. 6, 32534 (2016).

[33] Akpomedaye, D. E. & Ejechi, B. O. The hurdle effect of mild heat and two tropical spice extracts on the growth of three fungi in fruit juices. Food Res. Int. 31, 339-41 (1998).

[34] Sahayaraj, K. Antifeedant effect of some plant extracts on the Asian armyworm, Spodoptera litura (Fabricius). Curr. Sci. 74, 523-6 (1998).

[35] Wang, C. et al. Genome-wide analysis of local chromatin packing in Arabidopsis thaliana. Genome Res. 25, 246-56 (2015).

[36] Marçais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27, 764-70 (2011).

[37] Vurture, G. W. et al. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics 33, 2202-4 (2017).

[38] Chin, C. S. et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat. Methods 13, 1050-4 (2016).

[39] Chin, C. S. et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 10, 563-9 (2013).

[40] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-60 (2009).

[41] Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963 (2014).

[42] Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 3, 95-98 (2016).

[43] Dudchenko, O. et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science 356, 92-95 (2017).

[44] Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-2 (2015).

[45] Ou, S., Chen, J. & Jiang, N. Assessing genome assembly quality using the LTR Assembly Index (LAI). Nucleic Acids Res. 46, e126 (2018).

[46] Tarazona, S., García-Alcalde, F., Dopazo, J., Ferrer, A. & Conesa, A. Differential expression in RNA-seq: a matter of depth. Genome Res. 21, 2213-23 (2011).

[47] Slater, G. S. & Birney, E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 6, 31 (2005).

[48] Stanke, M., Keller, O., Gunduz, I., Hayes, A., Waack, S. & Morgenstern, B. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, W435-439 (2006).

[49] Holt, C. & Yandell, M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinform. 12, 491 (2011).

[50] Gremme, G., Steinbiss, S. & Kurtz, S. GenomeTools: a comprehensive software library for efficient processing of structured genome annotations. IEEE/ACM Trans. Comput. Biol. Bioinform. 10, 645-56 (2013).

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

[52] Blanc, G. & Wolfe, K. H. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell 16, 1667-78 (2004).

[53] Li, L., Stoeckert, C. J. Jr & Roos, D. S. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-89 (2003).

[54] Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403-10 (1990).

[55] Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-80 (2013).

[56] Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W. & Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307-21 (2010).

[57] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-91 (2007).

[58] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-71 (2006).

[59] Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49 (2012).

[60] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 9, 559 (2008).

[61] Song, X. et al. Deciphering the high-quality genome sequence of coriander that causes controversial feelings. Plant Biotechnol. J. 18, 1444-56 (2020).