The chromosome-level Hemerocallis citrina Borani genome provides new insights into the rutin biosynthesis and the lack of colchicine

Abstract

Hemerocallis citrina Borani (huang hua cai in Chinese) is an important horticultural crop whose flower buds are widely consumed as a delicious vegetable in Asia. Here we assembled a high-quality reference genome of H. citrina using single-molecule sequencing and Hi-C technologies. The genome assembly was 3.77 Gb and consisted of 3183 contigs with a contig N50 of 2.09 Mb, which were further clustered into 11 pseudochromosomes. A larger portion (3.25 Gb or 86.20%) was annotated as a repetitive content and 54,295 protein-coding genes were annotated in the genome. Genome evolution analysis showed that H. citrina experienced a recent whole-genome duplication (WGD) event at ~15.73 million years ago (Mya), which was the main factor leading to many multiple copies of orthologous genes. We used this reference genome to predict 20 genes involved in the rutin biosynthesis pathway. Moreover, our metabolomics data revealed neither colchicine nor its precursors in H. citrina, challenging the long-standing belief that this alkaloid causes poisoning by the plant. The results of our disruptive research are further substantiated by our genomic finding that H. citrina does not contain any genes involved in colchicine biosynthesis. The high-quality genome lays a solid foundation for genetic research and molecular breeding of H. citrina.

References

[1] Liu, J. et al. Systematic identification metabolites of Hemerocallis citrina Borani by high-performance liquid chromatography/quadrupole-time-of-flight mass spectrometry combined with a screening method. J. Pharm. Biomed. Anal. 186, 113314 (2020).

[2] Ma, G. et al. iTRAQ-based quantitative proteomic analysis reveals dynamic changes during daylily flower senescence. Planta 248, 859-873 (2018).

[3] Li, C. F. et al. Evaluation of the toxicological properties and anti-inflammatory mechanism of Hemerocallis citrina in LPS-induced depressive-like mice. Biomed. Pharmacother. 91, 167-173 (2017).

[4] Yang, R. F., Geng, L. L., Lu, H. Q. & Fan, X. D. Ultrasound-synergized electrostatic field extraction of total flavonoids from Hemerocallis citrina baroni. Ultrason. Sonochem. 34, 571-579 (2017).

[5] Lin, S. H. et al. The antidepressant-like effect of ethanol extract of daylily flowers (Jīn Zhēn Huā) in rats. J. Tradit. Complement. Med. 3, 53-61 (2013).

[6] Wang, J. et al. Ethyl acetate fraction of Hemerocallis citrina Baroni decreases tert-butyl hydroperoxide-induced oxidative stress damage in BRL-3A cells. Oxid. Med. Cell Longev. 2018, 1-13 (2018).

[7] Tian, H. et al. Effects of phenolic constituents of daylily flowers on corticosterone-and glutamate-treated PC12 cells. BMC Complement. Alter. Med. 17, 69 (2017).

[8] Xu, P. et al. Antidepressant-like effects and cognitive enhancement of the total phenols extract of Hemerocallis citrina Baroni in chronic unpredictable mild stress rats and its related mechanism. J. Ethnopharmacol. 194, 819-826 (2016).

[9] Tang, M. N., Liu, X. B., Huang, J. L., Deng, F. M. & Zeng, J. G. Questioning and arguable research on edible Hemerocallis citrina containing colchicine. Chin. Tradit. Herb. Drugs 047, 3293-3300 (2016).

[10] Li, S. F. et al. Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus. Hortic. Res. 7, 48 (2020).

[11] Zhang, L. et al. The tartary buckwheat genome provides insights into rutin biosynthesis and abiotic stress tolerance. Mol. Plant. 10, 1224-1237 (2017).

[12] Nett, R. S., Lau, W. & Sattely, E. S. Discovery and engineering of colchicine alkaloid biosynthesis. Nature 584, 148-153 (2020).

[13] Klein, G. & Soos, G. Der mikrochemische Nachweis der Alkaloide in der Pflanze. Oesterr Bot. Z. 78, 157-163 (1929).

[14] Traub, H. P. Colchicine poisoning in relation to Hemerocallis and some other plants. Science 110, 686-687 (1949).

[15] Li, Z. H. Plant poisoning. Barefoot Dr. Mag. 8, 44-45 (1977).

[16] Zong, W., Zhang, L. & Wang, M. Z. (eds) Food Safety (Chemical Industry, 2016).

[17] Zhang, Z. J. et al. (eds) Introduction of Food Safety (Chemical Industry, 2015).

[18] Peng, W. X., Pan, T., Yuan, Y. Y. & Wang, L. (eds) Food Safety and Food Poisoning-First Aid Knowledge (GuiZhou, 2012).

[19] Zhou, C. Q. et al. (eds) Food Nutrition (China Metrology, 2006).

[20] Hong, Y. F., Cheng, Z. W., Li, J. H. & Hu, C. On different methods to treat the fresh Hemerocallis citrina and lead to the change of colchicine. J. Hunan Agric. Univ. 29, 500-502 (2003).

[21] Zhang, N. et al. Optimization of HPLC detection system for colchicine content in flower buds of Hemerocallis. J. Agric. Univ. Hebei. 9, 48-54 (2017).

[22] Doyle, J. J. & Doyle, J. L. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19, 11-15 (1987).

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

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

[25] Roach, M. J., Schmidt, S. & Borneman, A. R. Purge Haplotigs: synteny reduction for third-gen diploid genome assemblies. BMC Bioinformatics 19, 460 (2018).

[26] Yin, D. et al. Genome of an allotetraploid wild peanut Arachis monticola: a de novo assembly. Gigascience 7, giy066 (2018).

[27] Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at arXiv:1303.3997 (2013).

[28] Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119-1125 (2013).

[29] McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303 (2010).

[30] Simão, F. A. et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[31] Tarailo‐Graovac, M. & Chen, N. Using repeatMasker to identify repetitive elements in genomic gequences. Curr. Protoc. Bioinformatics 25, 4.10.1-4.10.14 (2004).

[32] Lowe, T. M. & Eddy, S. R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955-964 (1997).

[33] Nawrocki, E. P., Kolbe, D. L. & Eddy, S. R. Infernal 1.0: inference of RNA alignments. Bioinformatics 25, 1335-1337 (2009).

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

[35] Korf, I. Gene finding in novel genomes. BMC Bioinformatics 5, 59 (2004).

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

[37] Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644-652 (2011).

[38] Haas, B. J. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 31, 5654-5666 (2003).

[39] Cantarel, B. L. et al. MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 18, 188-196 (2008).

[40] Camacho, C. et al. BLAST plus: architecture and applications. BMC Bioinformatics.BioMed. Cent. 10, 1 (2009).

[41] Boeckmann, B. et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 31, 365-370 (2003).

[42] Mistry, J., Bateman, A. & Finn, R. D. Predicting active site residue annotations in the Pfam database. BMC Bioinformatics 8, 298 (2007).

[43] Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236-1240 (2014).

[44] Hotz, G. C., Forslund, K., Eddy, S. R., Sonnhammer, E. L. & Bateman, A. The Pfam protein families database. Nucleic Acids Res. 36, S281-S288 (2008).

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

[46] Conesa, A. & Götz, S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int. J. Plant Genomics. 2008, 619832 (2008).

[47] Kanehisa, M., Goto, S., Sato, Y., Furumichi, M. & Tanabe, M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109-D114 (2012).

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

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

[50] Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312-1313 (2014).

[51] Han, M. V., Thomas, G. W., Lugo-Martinez, J. & Hahn, M. W. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol. Biol. Evol. 30, 1987-1997 (2013).