A chromosome-scale reference genome of Aquilegia oxysepala var. kansuensis

Abstract

The genus Aquilegia (Ranunculaceae) has been cultivated as ornamental and medicinal plants for centuries. With petal spurs of strikingly diverse size and shape, Aquilegia has also been recognized as an excellent system for evolutionary studies. Pollinator‐mediated selection for longer spurs is believed to have shaped the evolution of this genus, especially the North American taxa. Recently, however, an opposite evolutionary trend was reported in an Asian lineage, where multiple origins of mini- or even nonspurred morphs have occurred. Interesting as it is, the lack of genomic resources has limited our ability to decipher the molecular and evolutionary mechanisms underlying spur reduction in this special lineage. Using long-read sequencing (PacBio Sequel), in combination with optical maps (BioNano DLS) and Hi–C, we assembled a high-quality reference genome of A. oxysepala var. kansuensis, a sister species to the nonspurred taxon. The final assembly is approximately 293.2 Mb, 94.6% (277.4 Mb) of which has been anchored to 7 pseudochromosomes. A total of 25,571 protein-coding genes were predicted, with 97.2% being functionally annotated. When comparing this genome with that of A. coerulea, we detected a large rearrangement between Chr1 and Chr4, which might have caused the Chr4 of A. oxysepala var. kansuensis to partly deviate from the “decaying” path that was taken before the split of Aquilegia and Semiaquilegia. This high-quality reference genome is fundamental to further investigations on the development and evolution of petal spurs and provides a strong foundation for the breeding of new horticultural Aquilegia cultivars.

References

[1] Whittall, J. B., Medina-Marino, A., Zimmer, E. A. & Hodges, S. A. Generating single-copy nuclear gene data for a recent adaptive radiation. Mol. Phylogenet. Evol. 39, 124-134 (2006).

[2] Kramer, E. M. & Hodges, S. A. Aquilegia as a model system for the evolution and ecology of petals. Philos. Trans. R. Soc. Lond. B 365, 477-490 (2010).

[3] Kramer, E. M. Aquilegia: a new model for plant development, ecology, and evolution. Annu. Rev. Plant. Biol. 60, 261-277 (2009).

[4] Hodges, S. A. & Derieg, N. J. Adaptive radiations: from field to genomic studies. Proc. Natl Acad. Sci. USA 106, 9947-9954 (2009).

[5] Whittall, J. B. & Hodges, S. A. Pollinator shifts drive increasingly long nectar spurs in columbine flowers. Nature 447, 706-709 (2007).

[6] Hodges, S. A., Whittall, J. B., Fulton, M. & Yang, J. Y. Genetics of floral traits influencing reproductive isolation between Aquilegia formosa and Aquilegia pubescens. Am. Nat. 159, S51-S60 (2002).

[7] Fulton, M. & Hodges, S. A. Floral isolation between Aquilegia formosa and Aquilegia pubescens. Philos. Trans. R. Soc. Lond. B 266, 2247-2252 (1999).

[8] Huang, L. et al. Genetic diversity and evolutionary history of four closely related Aquilegia species revealed by 10 nuclear gene fragments. J. Syst. Evol. 56, 129-138 (2018).

[9] Xue, C. et al. Morphological variation pattern of Aquilegia ecalcarata and its relatives. J. Syst. Evol. 58, 221-233 (2020).

[10] 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, e36426 (2018).

[11] Fior, S. et al. Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions. N. Phytol. 198, 579-592 (2013).

[12] Zhu, R. R., Gao, Y., Zhang, Q. X. & Debener, T. Quantitative trait locus mapping of floral and related traits using an F₂ population of Aquilegia . Plant Breed. 133, 153-161 (2014).

[13] Yant, L., Collani, S., Puzey, J., Levy, C. & Kramer, E. M. Molecular basis for three-dimensional elaboration of the Aquilegia petal spur. Proc. R. Soc. Lond. B 282, 20142778 (2015).

[14] Whittall, J. B., Voelckel, C., Kliebenstein, D. J. & Hodges, S. A. Convergence, constraint and the role of gene expression during adaptive radiation: floral anthocyanins in Aquilegia. Mol. Ecol. 15, 4645-4657 (2006).

[15] Voelckel, C., Borevitz, J. O., Kramer, E. M. & Hodges, S. A. Within and between whorls: comparative transcriptional profiling of Aquilegia and Arabidopsis. PLoS ONE 5, e9735 (2010).

[16] Sharma, B., Yant, L., Hodges, S. A. & Kramer, E. M. Understanding the development and evolution of novel floral form in Aquilegia. Curr. Opin. Plant Biol. 17, 22-27 (2014).

[17] Puzey, J. R., Gerbode, S. J., Hodges, S. A., Kramer, E. M. & Mahadevan, L. Evolution of spur-length diversity in Aquilegia petals is achieved solely through cell-shape anisotropy. Proc. R. Soc. Lond. B 279, 1640-1645 (2012).

[18] Ballerini, E. S., Kramer, E. M. & Hodges, S. A. Comparative transcriptomics of early petal development across four diverse species of Aquilegia reveal few genes consistently associated with nectar spur development. BMC Genomics 20, 668 (2019).

[19] Li, L., Luo, Y., Yang, C., Deng, J. & Erst, A. S. Aquilegia yangii (Ranunculaceae), a new species from western China. Phytotaxa 348, 289-296 (2018).

[20] Dolezel, J. & Bartos, J. Plant DNA flow cytometry and estimation of nuclear genome size. Ann. Bot. 95, 99-110 (2005).

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

[22] Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589-595 (2010).

[23] Simao, 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).

[24] Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078-2079 (2009).

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

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

[27] Emms, D. M. & Kelly, S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20, 238 (2019).

[28] Smaczniak, C., Immink, R. G. H., Angenent, G. C. & Kaufmann, K. Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies. Development 139, 3081-3098 (2012).

[29] Theissen, G., Melzer, R. & Rumpler, F. MADS-domain transcription factors and the floral quartet model of flower development: linking plant development and evolution. Development 143, 3259-3271 (2016).

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

[31] Akoz, G. & Nordborg, M. The Aquilegia genome reveals a hybrid origin of core eudicots. Genome Biol. 20, 256 (2019).

[32] Nevado, B., Ramos-Onsins, S. E. & Perez‐Enciso, M. Resequencing studies of nonmodel organisms using closely related reference genomes: optimal experimental designs and bioinformatics approaches for population genomics. Mol. Ecol. 23, 1764-1779 (2014).

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

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

[35] Boetzer, M. & Pirovano, W. SSPACE-LongRead: scaffolding bacterial draft genomes using long read sequence information. BMC Bioinform. 15, 211 (2014).

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

[37] Wingett, S. et al. HiCUP: pipeline for mapping and processing Hi-C data. F1000Res 4, 1310 (2015).

[38] Dudchenko, O. et al. The Juicebox Assembly Tools module facilitates de novo assembly of mammalian genomes with chromosome-length scaffolds for under $1000. Preprint at bioRxiv https://doi.org/10.1101/254797 (2018).

[39] Tarailo‐Graovac, M. & Chen, N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinform. Chapter 4, Unit 4.10 (2009).

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

[41] Kalvari, I. et al. Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families. Nucleic Acids Res. 46, D335-D342 (2018).

[42] Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933-2935 (2013).

[43] Birney, E., Clamp, M. & Durbin, R. GeneWise and genomewise. Genome Res. 14, 988-995 (2004).

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

[45] Alioto, T., Blanco, E., Parra, G. & Guigó, R. Using geneid to identify genes. Curr. Protoc. Bioinform. 64, e56 (2018).

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

[47] Majoros, W. H., Pertea, M. & Salzberg, S. L. TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics 20, 2878-2879 (2004).

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

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

[50] Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562-578 (2012).

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

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

[53] Marcais, G. et al. MUMmer4: a fast and versatile genome alignment system. PLoS Comput. Biol. 14, e1005944 (2018).

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

[55] Eddy, S. Profile hidden Markov models. Bioinformatics 14, 755-763 (1998).

[56] Suyama, M., Torrents, D. & Bork, P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 34, W609-W612 (2006).

[57] Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870-1874 (2016).

[58] Wickham, H. ggplot2: Elegant Graphics for Data Analysis. (Springer, 2016).

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

[60] Capella-Gutiérrez, S., Silla-Martínez, J. M. & Gabaldón, T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972-1973 (2009).