Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus

Abstract

Asparagus setaceus is a popular ornamental plant cultivated in tropical and subtropical regions globally. Here, we constructed a chromosome-scale reference genome of A. setaceus to facilitate the investigation of its genome characteristics and evolution. Using a combination of Nanopore long reads, Illumina short reads, 10× Genomics linked reads, and Hi-C data, we generated a high-quality genome assembly of A. setaceus covering 710.15 Mb, accounting for 98.63% of the estimated genome size. A total of 96.85% of the sequences were anchored to ten superscaffolds corresponding to the ten chromosomes. The genome of A. setaceus was predicted to contain 28,410 genes, 25,649 (90.28%) of which were functionally annotated. A total of 65.59% of the genome was occupied by repetitive sequences, among which long terminal repeats were predominant (42.51% of the whole genome). Evolutionary analysis revealed an estimated divergence time of A. setaceus from its close relative A. officinalis of ~9.66 million years ago, and A. setaceus underwent two rounds of whole-genome duplication. In addition, 762 specific gene families, 96 positively selected genes, and 76 resistance (R) genes were detected and functionally predicted in A. setaceus. These findings provide new knowledge about the characteristics and evolution of the A. setaceus genome, and will facilitate comparative genetic and genomic research on the genus Asparagus.

References

[1] Dahlgren, R. M. T. et al. (eds) The Families of the Monocotyledons (Springer-Verlag Press, 1985).

[2] Clifford, H. T. & Conran, J. G. Asparagaceae. in Flora of Australia (ed George, A. S.) (Australian Government Publishing Service, 1987).

[3] Kubituki, K. & Rudall, P. J. Asparagaceae. in The Families and Genera of Vascular Plants (ed Kubituki, K. ) (Springer-Verlag Press, 1998).

[4] Ellison, J. H. & Kinelski, J. J. Greenwich, a male asparagus hybrid. HortScience 21, 1249 (1986).

[5] Fukuda, T. et al. Molecular phylogeny of the genus Asparagus (Asparagaceae) inferred from plastid petB intron and petD-rpoA intergenic spacer sequences. Plant Spec. Biol. 20, 121-132 (2005).

[6] Li, S. F. et al. Analysis of transposable elements in the genome of Asparagus officinalis from high coverage sequence data. PLoS One 9, e97189 (2014).

[7] Li, S. F. et al. Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis). BMC Plant Biol. 17, 143 (2017).

[8] Harkess, A. et al. Sex-biased gene expression in dioecious garden asparagus (Asparagus officinalis). N. Phytol. 207, 883-892 (2015).

[9] Murase, K. et al. MYB transcription factor gene involved in sex determination in Asparagus officinalis. Genes Cells 22, 115-123 (2016).

[10] Dong, T. et al. Anthocyanins accumulation and molecular analysis of correlated genes by metabolome and transcriptome in green and purple asparagus (Asparagus officinalis L.). Food Chem. 271, 18-28 (2019).

[11] Harkess, A. et al. The asparagus genome sheds light on the origin and evolution of a young Y chromosome. Nat. Commun. 8, 1279 (2017).

[12] Pindel, A. Regeneration capacity of Asparagus setaceus (Kunth) Jessop ‘Pyramidalis’ in in vitro cultures. Acta Sci. Pol. Hortorum Cultus 16, 85-93 (2017).

[13] McGaw, L. J. & Eloff, J. N. Ethnoveterinary use of southern African plants and scientific evaluation of their medicinal properties. J. Ethnopharmacol. 119, 559-574 (2008).

[14] Bansal, R. K., Menzies, S. A. & Broadhurst, P. G. Screening of Asparagus species for resistance to Stemphylium leaf spot. NZ J. Agric. Res. 29, 539-545 (1986).

[15] Alberti, P. et al. Interspecific hybridization for Asparagus breeding. In Proc of the XLVIII Italian Society of Agricultural Genetics-SIFV-SIGA Joint Meeting (Lecce, Italy, 2004).

[16] Akter, S., Begum, K. N., Sultana, S. S. & Alam, S. S. Karyotype diversity in three Asparagus L. species. Cytologia 82, 551-557 (2017).

[17] Štajner, N., Bohanec, B. & Javornik, B. Genetic variability of economically important Asparagus species as revealed by genome size analysis and rDNA ITS polymorphisms. Plant Sci. 162, 931-937 (2002).

[18] Li, J. R. et al. Characterization of the complete chloroplast genome of Asparagus setaceus. Mitochondrial DNA B 4, 2639-2640 (2019).

[19] Kubota, S., Konno, I. & Kannos, A. Molecular phylogeny of the genus Asparagus (Asparagaceae) explains interspecific crossability between the garden asparagus (A. officinalis) and other Asparagus species. Theor. Appl. Genet. 124, 345-354 (2012).

[20] Norup, M. F. et al. Evolution of Asparagus L. (Asparagaceae): Out-of-South-Africa and multiple origins of sexual dimorphism. Mol. Phylogenet. Evol. 92, 25-44 (2015).

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

[22] Nowak, M. D. et al. The draft genome of Primula veris yields insights into the molecular basis of heterostyly. Genome Biol. 16, 12 (2015).

[23] Schatz, M. C., Witkowski, J. & McCombie, W. R. Current challenges in de novo plant genome sequencing and assembly. Genome Res. 13, 243 (2012).

[24] Yang, X. 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).

[25] Sun, D. et al. Draft genome sequence of cauliflower (Brassica oleracea L. var. botrytis) provides new insights into the C genome in Brassica species. Hortic. Res. 6, 82 (2019).

[26] Guan, R. et al. Draft genome of the living fossil Ginkgo biloba. Gigascience 5, 49 (2016).

[27] Wicker, T. et al. Impact of transposable elements on genome structure and evolution in bread wheat. Genome Biol. 19, 103 (2018).

[28] Samoluk, S. S. et al. Heterochromatin evolution in Arachis investigated through genome-wide analysis of repetitive DNA. Planta 249, 1405-1415 (2019).

[29] Li, S. F. et al. Chromosome evolution in connection with repetitive sequences and epigenetics in plants. Genes 8, 290 (2017).

[30] Uzunović, J., Josephs, E. B., Stinchcombe, J. R. & Wright, S. I. Transposable elements are important contributors to standing variation in gene expression in Capsella Grandiflora. Mol. Biol. Evol. 36, 1734-1745 (2019).

[31] Ming, R., Bendahmane, A. & Renner, S. S. Sex chromosomes in land plants. Annu. Rev. Plant Biol. 62, 485-514 (2011).

[32] Telgmann-Rauber, A., Jamsari, A., Kinney, M. S., Pires, J. C. & Jung, C. Genetic and physical maps around the sex-determining M-locus of the dioecious plant Asparagus. Mol. Genet. Genomics 278, 221-234 (2007).

[33] Wei, C. et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl Acad. Sci. USA 115, E4151-E4158 (2018).

[34] Ming, R. et al. The pineapple genome and the evolution of CAM photosynthesis. Nat. Genet. 47, 1345-1442 (2015).

[35] Meyers, B. C., Kaushik, S. & Nandety, R. S. Evolving disease resistance genes. Curr. Opin. Plant Biol. 8, 129-134 (2005).

[36] Die, J. V., Castro, P., Millán, T. & Gil, J. Segmental and tandem duplication driving the recent NBS-LRR gene expansion in the asparagus genome. Genes 9, 568 (2018).

[37] Li, S. F. et al. The landscape of transposable elements and satellite DNAs in the genome of a dioecious plant spinach (Spinacia oleracea L.). Mol. DNA 10, 3 (2019).

[38] Belton, J. M. et al. Hi-C: a comprehensive technique to capture the conformation of genomes. Methods 58, 268-276 (2012).

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

[40] Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722 (2017).

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

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

[43] Coombe, L. et al. ARKS: chromosome-scale scaffolding of human genome drafts with linked read kmers. BMC Bioinform. 19, 234 (2018).

[44] Warren, R. L. et al. LINKS: scalable, alignment-free scaffolding of draft genomes with long reads. Gigascience 4, 35 (2015).

[45] Pryszcz, L. P. & Gabaldon, T. Redundans: an assembly pipeline for highly heterozygous genomes. Nucleic Acids Res. 44, e113 (2016).

[46] Langmead, B. & Salzberg, S. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

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

[48] Keilwagen, J. et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 44, e89 (2016).

[49] Campbell, M. A., Hass, B. J., Hamilton, J. P., Mount, S. M. & Buell, C. R. Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics 7, 327 (2006).

[50] Stanke, M., Schöffmann, O., Morgenstern, B. & Waack, S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinform. 7, 62 (2006).

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

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

[54] Griffiths-Jones, S. et al. Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 33, D121-D124 (2005).

[55] Lagesen, K. et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100-3108 (2007).

[56] Bao, W., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mol. DNA 6, 11 (2015).

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

[58] Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-780 (2013).

[59] Castresana, J . Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540-552 (2000).

[60] Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688-2690 (2006).

[61] Yang, Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13, 555-556 (1997).

[62] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

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

[64] Marchler-Bauer, A. et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 43, D222-D226 (2015).

[65] Finn, R. D. et al. Pfam: the protein families database. Nucleic Acids Res. 42, D222-D230 (2014).