Whole-genome resequencing of Osmanthus fragrans provides insights into flower color evolution

Abstract

Osmanthus fragrans is a well-known ornamental plant that has been domesticated in China for 2500 years. More than 160 cultivars have been found during this long period of domestication, and they have subsequently been divided into four cultivar groups, including the Yingui, Jingui, Dangui, and Sijigui groups. These groups provide a set of materials to study genetic evolution and variability. Here, we constructed a reference genome of O. fragrans ‘Liuyejingui’ in the Jingui group and investigated its floral color traits and domestication history by resequencing a total of 122 samples, including 119 O. fragrans accessions and three other Osmanthus species, at an average sequencing depth of 15×. The population structure analysis showed that these 119 accessions formed an apparent regional cluster. The results of linkage disequilibrium (LD) decay analysis suggested that varieties with orange/red flower color in the Dangui group had undergone more artificial directional selection; these varieties had the highest LD values among the four groups, followed by the Sijigui, Jingui, and Yingui groups. Through a genome-wide association study, we further identified significant quantitative trait loci and genomic regions containing several genes, such as ethylene-responsive transcription factor 2 and Arabidopsis pseudoresponse regulator 2, that are positively associated with petal color. Moreover, we found a frameshift mutation with a 34-bp deletion in the first coding region of the carotenoid cleavage dioxygenase 4 gene. This frameshift mutation existed in at least one site on both alleles in all varieties of the Dangui group. The results from this study shed light on the genetic basis of domestication in woody plants, such as O. fragrans.

References

[1] Zang, D. K. & Xiang, Q. B. Studies on Osmanthus fragrans cultivars. J. Nanjing Forestry Univ. 28, 7-13 (2004).

[2] Xiang, Q. B. & Liu, Y. L. An Illustrated Monograph of the Sweet Osmanthus Variety in China 93-260. (Zhejiang Science & Technology Press, Hangzhou, Zhejiang Province, 2008).

[3] He, Y. X., Yuan, W. J., Dong, M. F., Han, Y. J. & Shang, F. D. The first genetic map in sweet osmanthus (Osmanthus fragrans Lour.) using specific locus amplified fragment sequencing. Front. Plant Sci. 8, 1621 (2017).

[4] Duan, Y. F. et al. Genetic diversity of androdioecious Osmanthus fragrans (Oleaceae) cultivars using microsatellite markers. Appl. Plant Sci. 1, 1200092 (2013).

[5] Cai, X. et al. Analysis of aroma-active compounds in three sweet osmanthus (Osmanthus fragrans) cultivars by gas-chromatogolfactometry and GC-mass spectrometry. J. Zhejiang. Univ. Sci. B. 15, 638-648 (2014).

[6] Fu, J. X. et al. The Emission of the Floral Scent of Four Osmanthus fragrans Cultivars in Response to Different Temperatures. Molecules 22, 430 (2017).

[7] Fu, J. X. et al. Identifcation of foral aromatic volatile compounds in 29 cultivars from four groups of Osmanthus fragrans by gas chromatography-mass spectrometry. Hortic. Environ. Biote. 60, 611-623 (2019).

[8] Han, Y. J., Liu, L. X., Dong, M. F. & Shang, F. D. cDNA cloning of the phytoene synthase (PSY) and expression analysis of PSY and carotenoid cleavage dioxygenase genes in Osmanthus fragrans. Biologia 68, 258-263 (2013).

[9] Han, Y. J. et al. Mechanism of floral scent production in Osmanthus fragrans and the production and regulation of its key floral constituents, β-ionone and linalool. Hortic. Res. 6, 106 (2019).

[10] Wang, Y. G. et al. Carotenoid accumulation and its contribution to flower coloration of Osmanthus fragrans. Front. Plant Sci. 9, 1499 (2018).

[11] Baldermann, S. et al. Functional characterization of a carotenoid cleavage dioxygenase 1 and its relation to the carotenoid accumulation and volatile emission during the floral development of Osmanthus fragrans Lour. J. Exp. Bot. 61, 2967-2977 (2010).

[12] Han, Y. J. et al. Characterization of OfWRKY3, a transcription factor that positively regulates the carotenoid cleavage dioxygenase gene OfCCD4 in Osmanthus fragrans. Plant Mol. Biol. 91, 485-496 (2016).

[13] Yang, L. et al. The chromosome-level quality genome provides insights into the evolution of the biosynthesis genes for aroma compounds of Osmanthus fragrans. Hortic. Res. 5, 72 (2018).

[14] Yang, K. M. Chinese Osmanthus 47-67 (China Forestry Publishing House, 2020).

[15] Zou, J. J., Zhou, Y., Cai, X. & Wang, C. Y. Increase in DNA fragmentation and the role of ethylene and reactive oxygen species in petal senescence of Osmanthus fragrans. Postharvest Biol. Tec. 93, 97-105 (2014).

[16] Rogers, S. O. & Bendich, A. J. Plant Molecular Biology Manual 73-83 (eds Gelvin, S. B., Schilperoort, R. A. & Verma, D. P. S.) (Springer, 1989).

[17] Wang, B. et al. Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat. Comms. 7, 11708 (2016).

[18] Vaser, R., Sović, I., Nagarajan, N. & Šikić, M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 27, 737-746 (2017).

[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] Yang, X. et al. HTQC: a fast quality control toolkit for Illumina sequencing data. BMC Bioinforma. 14, 33 (2013).

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

[22] Tarailo-Graovac, M. & Chen, N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinformatics 4, 10.1-10.14 (2009).

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

[24] Bao, Z. R. & Eddy, S. R. Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 12, 1269-1276 (2002).

[25] Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, I351-I358 (2005).

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

[27] Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573-580 (1999).

[28] Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, W435-W439 (2006).

[29] Burge, C. & Karlin, S. Prediction of complete gene structures in human genomic DNA. J. Mol. Biol. 268, 78-94 (1997).

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

[31] Birney, E. & Durbin, R. Using GeneWise in the Drosophila annotation experiment. Genome Res. 10, 547-548 (2000).

[32] Yu, X. J. et al. Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup. Genomics 88, 745-751 (2006).

[33] Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105-1111 (2009).

[34] Trapnell, C. et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511-515 (2010).

[35] Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658-1659 (2006).

[36] Wu, T. D. & Colin, K. W. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics 21, 1859-1875 (2005).

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

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

[39] Mitchell, A. et al. The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res. 43, D213-D221 (2015).

[40] Kanehisa, M. et al. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109-D114 (2012).

[41] Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402 (1997).

[42] Camacho, C. et al. BLAST plus: architecture and applications. BMC Bioinform. 10, 421 (2009).

[43] Mistry, J. et al. Predicting active site residue annotations in the Pfam database. BMC Bioinform. 8, 298 (2007).

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

[45] Finn, R. D. et al. The Pfam protein families database. Nucleic Acids Res. 36, D281-D288 (2008).

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

[47] Conesa, A. & Gotz, S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int. J. Plant Genom. 2008, 1-12 (2008).

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

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

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

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

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

[53] Han, M. V. et al. 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).

[54] Kiełbasa, S. M. et al. Adaptive seeds tame genomic sequence comparison. Genome Res. 21, 487-493 (2011).

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

[56] Poplin, R. et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv. 24, 1-14 (2018).

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

[58] Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256-W259 (2019).

[59] Yang, J., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76-82 (2011).

[60] Alexander, D. H. et al. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655-1664 (2009).

[61] Zhang, C. et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 10, 1093 (2018).

[62] Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156-2158 (2011).

[63] Tang, Y. et al. GAPIT version 2: an enhanced integrated tool for genomic association and prediction. Plant Genome 9, 1-9 (2016).

[64] Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Statis. Soc. B. 57, 289-300 (1995).

[65] Dolezel, J., Greilhuber, J. & Suda, J. Estimation of nuclear DNA content in plants using flow cytometry. Nat. Protoc. 2, 2233-2244 (2007).

[66] Duan, Y. F. et al. Analysis on the current reseach situation of ancient sweet osmanthus and discussion on the origin and evolution of the sweet osmanthus cultivars in China. J. Hubei Univ. Natl 28, 375-379 (2010).

[67] Cun, Y. Z. & Wang, X. Q. Plant recolonization in the Himalaya from the southeastern Qinghai-Tibetan Plateau: Geographical isolation contributed to high population differentiation. Mol. Phylogenet. Evol. 56, 972-982 (2010).

[68] Gong, W. et al. Low genetic diversity and high genetic divergence caused by inbreeding and geographical isolation in the populations of endangered species Loropetalum subcordatum (Hamamelidaceae) endemic to China. Conserv. Genet. 11, 2281-2288 (2010).

[69] Chen, J. H. et al. Genome-wide analysis of Cushion willow provides insights into alpine plant divergence in a biodiversity hotspot. Nat. commun. 10, 5230 (2019).

[70] Wang, X. F. et al. Identification and RElationships of Osmanthus fragrans cultivar groups. J. Northeast Forestry Univ. 41, 71-74 (2013).

[71] Pan, Y. et al. Network inference analysis identifies an APRR2-like gene linked to pigment accumulation in tomato and pepper fruits. Plant Physiol. 161, 1476-1485 (2013).

[72] Oren, E. et al. The multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon. J. Exp. Bot. 70, 15 (2019).

[73] Brandi, F. et al. Study of’Redhaven’ peach and its white-fleshed mutant suggests a key role of CCD4 carotenoid dioxygenase in carotenoid and nor-isoprenoid volatile metabolism. BMC Plant Biol. 11, 24 (2011).

[74] Zheng, X. J. et al. Natural variation in CCD4 promoter underpins species-specific evolution of red coloration in citrus peel. Mol. Plant 12, 1294-1307 (2019).

[75] Bodin, P. et al. Expression of CCD4 gene involved in carotenoid degradation in yellow-flowered Petunia × hybrida. Sci. Hortic. 261, 108916 (2020).