The chromosome-level reference genome of Coptis chinensis provides insights into genomic evolution and berberine biosynthesis

Abstract

Coptis chinensis Franch, a perennial herb, is mainly distributed in southeastern China. The rhizome of C. chinensis has been used as a traditional medicine for more than 2000 years in China and many other Asian countries. The pharmacological activities of C. chinensis have been validated by research. Here, we present a de novo high-quality genome of C. chinensis with a chromosome-level genome of ~958.20 Mb, a contig N50 of 1.58 Mb, and a scaffold N50 of 4.53 Mb. We found that the relatively large genome size of C. chinensis was caused by the amplification of long terminal repeat (LTR) retrotransposons. In addition, a whole-genome duplication event in ancestral Ranunculales was discovered. Comparative genomic analysis revealed that the tyrosine decarboxylase (TYDC) and (S)-norcoclaurine synthase (NCS) genes were expanded and that the aspartate aminotransferase gene (ASP5) was positively selected in the berberine metabolic pathway. Expression level and HPLC analyses showed that the berberine content was highest in the roots of C. chinensis in the third and fourth years. The chromosome-level reference genome of C. chinensis provides important genomic data for molecular-assisted breeding and active ingredient biosynthesis.

References

[1] Li, X.K., Motwani, M., Tong, W., Bornmann, W. & Schwartz, G.K. Huanglian, A chinese herbal extract, inhibits cell growth by suppressing the expression of cyclin B1 and inhibiting CDC2 kinase activity in human cancer cells. Mol. Pharmacol. 58, 1287 (2000).

[2] Franzblau, S.G. & Cross, C. Comparative in vitro antimicrobial activity of Chinese medicinal herbs. J. Ethnopharmacol. 15, 279-288 (1986).

[3] Liu, B., Li, W., Chang, Y., Dong, W. & Ni, L. Extraction of berberine from rhizome of Coptis chinensis Franch using supercritical fluid extraction. J. Pharm. Biomed. Anal. 41, 1056-1060 (2006).

[4] Lin, C.C., Ng, L.T., Hsu, F.F., Shieh, D.E. & Chiang, L.C. Cytotoxic effects of Coptis chinensis and Epimedium sagittatum extracts and their major constituents (berberine, coptisine and icariin) on hepatoma and leukaemia cell growth. Clin. Exp. Pharmacol. Physiol. 31, 65-69 (2004).

[5] Yan, D., Jin, C., Xiao, X.H. & Dong, X.P. Antimicrobial properties of berberines alkaloids in Coptis chinensis Franch by microcalorimetry. J. Biochem. Biophys. Methods 70, 845-849 (2008).

[6] Diogo, C.V. et al. Berberine as a promising safe anti-cancer agent-is there a role for mitochondria? Curr. Drug Targets 12, 850-859 (2011).

[7] Song, X., Pan, Y., Li, L., Wu, X. & Wang, Y. Composition and diversity of rhizosphere fungal community in Coptis chinensis Franch. continuous cropping fields. PLoS ONE 13, e0193811 (2018).

[8] Yan, H. et al. Anti-ultraviolet radiation effects of Coptis chinensis and Phellodendron amurense glycans by immunomodulating and inhibiting oxidative injury. Int. J. Biol. Macromol. 48, 720-725 (2011).

[9] Chang, W. Non-coding RNAs and berberine: a new mechanism of its antidiabetic activities. Eur. J. Pharmacol. 795, 8-12 (2017).

[10] Yang, T.C. et al. Alkaloids from Coptis chinensis root promote glucose uptake in C2C12 myotubes. Fitoterapia 93, 239-244 (2014).

[11] Liu, Y.M. et al. Berberine attenuates depressive-like behaviors by suppressing neuro-inflammation in stressed mice. Brain Res. Bull. 134, 220-227 (2017).

[12] Wang, L. et al. New enantiomeric isoquinoline alkaloids from coptis chinensis. Phytochemistry Lett. 7, 89-92 (2014).

[13] Chen, H.B. et al. Anti-inflammatory activity of coptisine free base in mice through inhibition of NF-κB and MAPK signaling pathways. Eur. J. Pharmacol. 811, 222-231 (2017).

[14] Wu, J. et al. Coptisine from Coptis chinensis inhibits production of inflammatory mediators in lipopolysaccharide-stimulated RAW 264.7 murine macrophage cells. Eur. J. Pharmacol. 780, 106-114 (2016).

[15] Wang, D., Liu, Z., Guo, M. & Liu, S. Structural elucidation and identification of alkaloids in Rhizoma Coptidis by electrospray ionization tandem mass spectrometry. J. Mass Spectrom. 39, 1356 (2004).

[16] Xie, W., Gu, D., Li, J., Cui, K. & Zhang, Y. Effects and action mechanisms of berberine and rhizoma coptidis on gut microbes and obesity in high-fat diet-fed C57BL/6J mice. Plos ONE 6, e24520 (2011).

[17] Dong, X. et al. De novo assembly of a wild pear (Pyrus betuleafolia) genome. Plant Biotechnol. J. 18, 581-595 (2020).

[18] Huang, X. et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 42, 961-967 (2010).

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

[20] Parra, G., Bradnam, K. & Korf, I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics 23, 1061-1067 (2007).

[21] Mochamad, A.F. et al. KNApSAcK family databases: integrated metabolite-plant species databases for multifaceted plant research. Plant Cell Physiol. 53, e1 (2012).

[22] Pei, L. et al. Genome and transcriptome of Papaver somniferum Chinese landrace CHM indicates that massive genome expansion contributes to high benzylisoquinoline alkaloid biosynthesis. Hortic. Res. 8, 5 (2021).

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

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

[25] Liu, X. et al. The genome of medicinal plant macleaya cordata provides new insights into benzylisoquinoline alkaloids metabolism. Mol. Plant 10, 975-989 (2017).

[26] Paterson, A.H. et al. Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 492, 423-427 (2012).

[27] Filiault, D.L. et al. The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history. Elife 7, https://doi.org/10.7554/eLife.36426 (2018).

[28] He, S.M. et al. Identification and characterization of genes involved in benzylisoquinoline alkaloid biosynthesis in Coptis species. Front. Plant Sci. 9, 731 (2018).

[29] Xie, J. et al. A chromosome-scale reference genome of Aquilegia oxysepala var. kansuensis. Hortic. Res. 7, 113 (2020).

[30] Sun, Q., He, X.H., Mu, J. & Cai, Y.F. Genome survey sequencing and analyzing of Coptis chinensis based on next-generation sequencing. J. Sichuan Univ. (Nat. Sci. Ed.) 51, 1330-1334 (2014).

[31] Meyers, B.C., Tingey, S.V. & Morgante M. Abundance, distribution, and transcriptional activity of repetitive elements in the maize genome. Genome Res. 11, 1660-1676 (2001).

[32] Ming, R., Vanburen, R., Wai, C.M., Tang, H. & Yu, Q. The pineapple genome and the evolution of CAM photosynthesis. Nat. Genet. 47, 1435-1442 (2015).

[33] 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).

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

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

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

[37] Mostovoy, Y. et al. A hybrid approach for de novo human genome sequence assembly and phasing. Nat. Methods 13, 587-590 (2016).

[38] Yaffe, E. & Tanay, A. Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture. Nat. Genet. 43, 1059-1065 (2011).

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

[40] 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-3212 (2015).

[41] Price, A.L., Jones, N.C. & Pevzner, P.A. De novo identification of repeat families in large genomes. Bioinformatics 21(Suppl 1), i351-358 (2005).

[42] Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265-268 (2007).

[43] Jurka, J. et al. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 110, 462-467 (2005).

[44] Haas, B.J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to Assemble Spliced Alignments. Genome Biol. 9, R7 (2008).

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

[46] Zdobnov, E.M. & Apweiler, R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17, 847-848 (2001).

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

[48] Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639-1645 (2009).

[49] Guo, L. et al. The opium poppy genome and morphinan production. Science 362, 343-347 (2018).

[50] Yang, Z. & Nielsen, R. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evol. 17, 32-43 (2000).

[51] He, S.M. et al. Identification of candidate genes involved in isoquinoline alkaloids biosynthesis in Dactylicapnos scandens by transcriptome analysis. Sci. Rep. 7, 9119 (2017).