The Chrysanthemum lavandulifolium genome and the molecular mechanism underlying diverse capitulum types

Abstract

Cultivated chrysanthemum (Chrysanthemum × morifolium Ramat.) is a beloved ornamental crop due to the diverse capitula types among varieties, but the molecular mechanism of capitulum development remains unclear. Here, we report a 2.60 Gb chromosome-scale reference genome of C. lavandulifolium, a wild Chrysanthemum species found in China, Korea and Japan. The evolutionary analysis of the genome revealed that only recent tandem duplications occurred in the C. lavandulifolium genome after the shared whole genome triplication (WGT) in Asteraceae. Based on the transcriptomic profiling of six important developmental stages of the radiate capitulum in C. lavandulifolium, we found genes in the MADS-box, TCP, NAC and LOB gene families that were involved in disc and ray floret primordia differentiation. Notably, NAM and LOB30 homologs were specifically expressed in the radiate capitulum, suggesting their pivotal roles in the genetic network of disc and ray floret primordia differentiation in chrysanthemum. The present study not only provides a high-quality reference genome of chrysanthemum but also provides insight into the molecular mechanism underlying the diverse capitulum types in chrysanthemum.

References

[1] Song XB, Xu Y, Gao K et al. High-density genetic map construction and identification of loci controlling flower-type traits in chrysanthemum (chrysanthemum × morifolium Ramat.) . Horticulture Research. 2020; 7: 108.

[2] He SM, Dong X, Zhang G et al. High quality genome of Erigeron breviscapus provides a reference for herbal plants in Asteraceae. Mol Ecol Resour. 2020; 00: 1-17.

[3] Song C, Liu Y, Song A et al. The chrysanthemum nankingense genome provides insights into the evolution and diversification of chrysanthemum flowers and medicinal traits. Mol Plant. 2018; 11: 1482-91.

[4] Hirakawa H, Sumitomo K, Hisamatsu T et al. De novo whole-genome assembly in chrysanthemum seticuspe, a model species of chrysanthemums, and its application to genetic and gene discovery analysis. DNA Res. 2019; 26: 195-203.

[5] Wellmer F, Riechmann JL . Gene networks controlling the initiation of flower development. Trends Genet. 2010; 26: 519-527.

[6] Thomson B, Wellmer F . Molecular regulation of flower development. Curr Top Dev Biol. 2019; 131: 185-210.

[7] Krizek BA, Fletcher JC . Molecular mechanisms of flower development: an armchair guide. Nature Rev Genet. 2005; 6: 688-98.

[8] Chen F, Song Y, Li X et al. Genome sequences of horticultural plants: past, present, and future. Horticulture Research. 2019; 6: 112.

[9] Krizek BA . eLS . London: Wiley; 2020.

[10] Zhang QG, Liu KW, Li Z et al. The Apostasia genome and the evolution of orchids . Nature. 2017; 549: 379-83.

[11] Li MM, Zhang D, Gao Q et al. Genome structure and evolution of Antirrhinum majus L. Nature Plants. 2019; 5: 174-83.

[12] Zhang LS, Chen F, Zhang X et al. The water lily genome and the early evolution of flowering plants. Nature. 2020; 577: 79-84.

[13] Ding L, Song A, Zhang X et al. The core regulatory networks and hub genes regulating flower development in Chrysanthemum morifolium . Plant Mol Biol. 2020; 103: 669-88.

[14] Elomaa P, Zhao Y, Zhang T et al. Flower heads in Asteraceae - recruitment of conserved developmental regulators to control the flower-like inflorescence architecture. Horticulture research. 2018; 5: 1-10.

[15] Zoulias N, Duttke SHC, Garcês H et al. Auxin and pattern formation of the Asteraceae flower head (capitulum). Plant Physiol. 2019; 179: 391-401.

[16] Li F, Lan W, Zhou Q et al. Reduced expression of CbUFO is associated with the phenotype of a flower-defective Cosmos bipinnatus . Int J Mol Sci. 2019; 20: 2503.

[17] Dai SL, Zhang CJ, Chen J et al. Advances of researches on phylogeny of Dendranthema and origin of chrysanthemum. Journal of Beijing Forestry University. 2002; 24: 234-8.

[18] Yang LW, Wen XH, Fu JX et al. ClCRY2 facilitates floral transition in chrysanthemum lavandulifolium by affecting the transcription of circadian clock-related genes under short-day photoperiods. Horticulture Research. 2018; 5: 58.

[19] Wen XH, Qi S, Huang H et al. The expression and interactions of ABCE-class and CYC2-like genes in the capitulum development of chrysanthemum lavandulifolium and C. × morifolium . Plant Growth Regul. 2019; 88: 205.

[20] Qi S, Yang L, Wen X et al. Reference gene selection for RT-qPCR analysis of flower development in Chrysanthemum morifolium and chrysanthemum lavandulifolium . Front Plant Sci. 2016; 7: 651.

[21] Chakraborty M, Baldwin-Brown JG, Long AD et al. Contiguous and accurate de novo assembly of metazoan genomes with modest long read coverage. Nucleic Acids Res. 2016; 44: e147.

[22] Roach MJ, Schmidt SA, Borneman AR . Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinformatics. 2018; 19: 460.

[23] Qiao X, Li Q, Yin H et al. Gene duplication and evolution in recurring polyploidization - diploidization cycles in plants. Genome Biol. 2019; 20: 38.

[24] Badouin H, Gouzy J, Grassa CJ et al. The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature. 2017; 546: 7675.

[25] Liu B, Yan J, Li W et al. Mikania micrantha genome provides insights into the molecular mechanism of rapid growth . Nature Communications. 2020; 11: 1.

[26] Spencer V, Kim M . Re"CYC"ling molecular regulators in the evolution and development of flower symmetry. Semin Cell Dev Biol. 2018; 79: 16-26.

[27] Chen J, Shen CZ, Guo YP et al. Patterning the Asteraceae Capitulum: duplications and differential expression of the flower symmetry CYC2-like genes. Frontiers Plant Sci. 2018; 9: 1-14.

[28] Huang CH, Zhang C, Liu M et al. Multiple polyploidization events across Asteraceae with two nested events in the early history revealed by nuclear phylogenomics. Mol Biol Evol. 2016; 33: 2820-35.

[29] Zadnikova P, Simon R . How boundaries control plant development. Curr Opin Plant Biol. 2014; 17: 116-25.

[30] Mara C, Manrique S, Cuesta C et al. CUP-SHAPED COTYLEDON1 (CUC1) and CUC2 regulate cytokinin homeostasis to determine ovule number in Arabidopsis . J Exp Bot. 2018; 69: 5169-76.

[31] Rebocho AB, Kennaway JR, Bangham JA et al. Formation and shaping of the antirrhinum flower through modulation of the CUP boundary gene. Curr Biol. 2017; 27: 2610-2622.e3.

[32] Zhao Y, Zhang T, Broholm SK et al. Evolutionary co-option of floral meristem identity genes for patterning of the flower-like Asteraceae inflorescence. Plant Physiol. 2016; 172: 284-96.

[33] Wen XH, Pu Y, Liu Y . Effects of N, P and K nutrients on the growth and development of chrysanthemum lavandulifolium based on BBCH scale. Advanced in Ornamental Horticulture of China. 2019; 1: 76-84.

[34] Ramsköld D, Luo S, Wang YC et al. Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat Biotechnol. 2012; 30: 777-82.

[35] Koren S, Walenz BP, Berlin K et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017; 27: 722-36.

[36] Drake JP, Berlin K, Koren S et al. Assembling large genomes with single-molecule sequencing and locality-sensitive hashing. Nat Biotechnol. 2015; 33: 623-30.

[37] Chin CS, Peluso P, Sedlazeck FJ et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods. 2016; 13: 1050-4.

[38] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

[39] Vaser R, Sović I, Nagarajan N et al. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 2017; 27: 737-46.

[40] Marcais G, Delcher AL, Phillippy AM et al. MUMmer4: a fast and versatile genome alignment system. PLoS Comput Biol. 2018; 14: e1005944.

[41] Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34: 3094-100.

[42] Rao SS, Huntley MH, Durand NC et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014; 159: 1665-80.

[43] Li H, Durbin R . Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010; 26: 589-95.

[44] Burton JN, Adey A, Patwardhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[45] Singh KS, Wu Y, Ghosh JS et al. OPEN RNA-sequencing reveals global transcriptomic changes in Nicotiana tabacum responding to topping and treatment of axillary-shoot control chemicals. Sci Rep. 2016; 5: 18148.

[46] Kim D, Pertea G, Trapnell C et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013; 14: 4.

[47] Trapnell C, Roberts A, Goff L et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012; 7: 562-78.

[48] Burge C, Karlin S . Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997; 268: 78-94.

[49] Stanke M, Waack S . Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003; 19: ii215-25.

[50] Majoros WH, Pertea M, Salzberg SL . TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004; 20: 2878-9.

[51] Blanco E, Genís P, Roderic G . Using GENEID to identify genes. Curr Protoc Bioinformatics. 2007; 18: 1.

[52] Korf I. Gene finding in novel genomes. BMC bioinformatics. 2004; 5: 59.

[53] Birney E, Clamp M, Durbin R . GeneWise and Genomewise. Genome Res. 2004; 14: 988-95.

[54] Kim D, Langmead B, Salzberg SL . HISAT: a fast-spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[55] Pertea M, Pertea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Bio/technology (Nature Publishing Company). 2015; 33: 290-5.

[56] Tang SYY, Alexandre L, Mark B . Identification of protein coding regions in RNA transcripts. Nucleic Acids Symp Ser. 2015; 43: e78.

[57] Haas BJ, Salzberg SL, Zhu W et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008; 9: R7.

[58] Griffiths-Jones S, Moxon S, Marshall M et al. Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2004; 33: D121-4.

[59] Lowe TM, Eddy SR . tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997; 25: 955-64.

[60] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; 25: 1.

[61] Buchfink B, Xie C, Huson DH . Fast and sensitive protein alignment using diamond. Nat Methods. 2015; 12: 59-60.

[62] Emms DM, Kelly S . OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015; 16: 157.

[63] Lam-Tung N, Schmidt HA, von Haeseler A et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology & Evolution. 2015; 1: 268-74.

[64] Rannala YB . Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006; 23: 212-26.

[65] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[66] Katoh K, Standley DM . MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molec Biol Evol. 2013; 30: 772-80.

[67] Rice P, Longden I, Bleasby A . EMBOSS: the European molecular biology open software suite. Trends Genet. 2000; 16: 276-7.

[68] Langfelder P, Horvath S . WGCNA: an R package for weighted correlation network analysis. Bioinformatics. 2008; 9: 559.

[69] Langfelder P, Zhang B, Horvath S et al. Defining clusters from a hierarchical cluster tree: the dynamic tree cut package for R. Bioinformatics. 2008; 24: 719-20.

[70] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194.

[71] Morgan NP, Dehal PS, Arkin AP et al. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009; 7: 1641-50.