The genome sequence of star fruit (Averrhoa carambola)

Abstract

Oxalidaceae is one of the most important plant families in horticulture, and its key commercially relevant genus, Averrhoa, has diverse growth habits and fruit types. Here, we describe the assembly of a high-quality chromosome-scale genome sequence for Averrhoa carambola (star fruit). Ks distribution analysis showed that A. carambola underwent a whole-genome triplication event, i.e., the gamma event shared by most eudicots. Comparisons between A. carambola and other angiosperms also permitted the generation of Oxalidaceae gene annotations. We identified unique gene families and analyzed gene family expansion and contraction. This analysis revealed significant changes in MADS-box gene family content, which might be related to the cauliflory of A. carambola. In addition, we identified and analyzed a total of 204 nucleotide-binding site, leucine-rich repeat receptor (NLR) genes and 58 WRKY genes in the genome, which may be related to the defense response. Our results provide insights into the origin, evolution and diversification of star fruit.

References

[1] Veldkamp, J. F. in Flora Malesiana (ed. van Steenis, C. G. C. J.) 155-178 (Noordhoff, Leyden, 1971).

[2] Khare, C. P. Indian Medicinal Plants 74-75 (Springer, New York, 2007).

[3] Gunawardena, D. C., Jayasinghe, L. & Fujimoto, Y. Phytotoxic constituents of the fruits of Averrhoa carambola. Chem. Nat. Compd. 51, 532-533 (2015).

[4] The Angiosperm Phylogeny Group . An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG IV. Bot. J. Linn. Soc. 181, 1-20 (2016).

[5] Stevens, P. F. Angiosperm Phylogeny Website. Version 14, July 2017 [and more or less continuously updated since] (2001 onwards).

[6] Heibl, C. & Renner, S. S. Distribution models and a dated phylogeny for Chilean Oxalis species reveal occupation of new habitats by different lineages, not rapid adaptive radiation. Syst. Biol. 61, 823-834 (2012).

[7] Sun, M., Naeem, R., Su, J. X., Cao, Z. Y. & Chen, Z. D. Phylogeny of the Rosidae: a dense taxon sampling analysis. J. Syst. Evol. 54, 363-391 (2016).

[8] Fieischmann, P., Watanabe, N. & Winterhalter, P. Enzymatic carotenoid cleavage in star fruit (Averrhoa carambola). Phytochemistry 63, 131-137 (2003).

[9] Neto, M. M. et al. Intoxication by star fruit (Averrhoa carambola) in 32 uraemic patients: treatment and outcome. Nephrol. Dial. Transplant. 18, 120-125 (2003).

[10] Carolino, R. O. G. et al. Convulsant activity and neurochemical alterations induced by a fraction obtained from fruit Averrhoa carambola (Oxalidaceae: Geraniales). Neurochem. Int. 46, 523-531 (2005).

[11] Shui, G. G. & Leong, L. P. Residue from star fruit as valuable source for functional food ingredients and antioxidant nutraceuticals. Food Chem. 97, 277-284 (2005).

[12] Bin-Ngah, A.W., Ahmad, I. & Hassan, A. Carambola Production, Processing and Marketing in Malaysia (Inter-American Society for Tropical Horticulture (IASTH), 1992).

[13] Khoo, H. et al. A review on underutilized tropical fruits in Malaysia. Guangxi Agr. Sci. 41, 698-702 (2010).

[14] Fukushima, K. et al. Genome of the pitcher plant Cephalotus reveals genetic changes associated with carnivory. Nat. Ecol. Evol. 1, 0059 (2017).

[15] Wu, Q. Study on Technology of Breeding polyploid of Averrhoa Carambola by Biotechnique. Master dissertation. Southwest Agricultural University, Chongqing (2002).

[16] Wang, C. et al. Genome-wide analysis of local chromatin packing in Arabidopsis thaliana. Genome Res. 25, 246 (2015).

[17] Cui, L. et al. Widespread genome duplications throughout the history of flowering plants. Genome Res. 16, 738-749 (2006).

[18] Van de Peer, Y., Maere, S. & Meyer, A. The evolutionary significance of ancient genome duplications. Nat. Rev. Genet. 10, 725-732 (2009).

[19] Blanc, G. & Wolfe, K. H. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell 16, 1667-1678 (2004).

[20] Cai, J. et al. The genome sequence of the orchid Phalaenopsis equestris. Nat. Genet. 47, 65-72 (2015).

[21] Parenicova, L. et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: new openings to the MADS world. Plant Cell 15, 1538-1551 (2003).

[22] Argout, X. et al. The cacao Criollo genome v2.0: an improved version of the genome for genetic and functional genomic studies. BMC Genomics 18, 730 (2017).

[23] Argout, X. et al. The genome of Theobroma cacao. Nat. Genet. 43, 101-108 (2011).

[24] Willmann, M. R. & Poethig, R. S. The effect of the floral repressor FLC on the timing and progression of vegetative phase change in Arabidopsis. Development 138, 677-685 (2011).

[25] Gregis, V., Sessa, A., Dorca-Fornell, C. & Kater, M. M. The Arabidopsis floral meristem identity genes AP1, AGL24 and SVP directly repress class B and C floral homeotic genes. Plant J. 60, 626-637 (2009).

[26] Oliveira, G. P. et al. Origin and development of reproductive buds in jabuticaba cv. Sabará (Plinia jaboticaba Vell). Sci. Horticulturae 249, 432-438 (2019).

[27] Wu, R. Kiwifruit SVP2, a dormancy regulator. Eur. J. Hort. Sci. 83, 231-235 (2018).

[28] Li, D. et al. A repressor complex governs the integration of flowering signals in Arabidopsis. Dev. Cell. 15, 110-120 (2008).

[29] Shao, Z. et al. Large-scale analyses of angiosperm nucleotide-binding site-leucine-rich repeat genes reveal three anciently diverged classes with distinct evolutionary patterns. Plant Physiol. 170, 2095-2109 (2016).

[30] Maleck, K. et al. The transcriptome of Arabidopsis thaliana during systemic acquired resistance. Nat. Genet. 26, 403-410 (2000).

[31] Zhang, L. S. et al. The water lily genome and the early evolution of flowering plants. Nature 577, 79-84 (2020).

[32] Phukan, U. et al. WRKY transcription factors: molecular regulation and stress responses in plants. Front Plant Sci. 7, 760 (2016).

[33] Wang, Q., Wang, M., Zhang, X., Hao, B., Kaushik, S. K. & Pan, Y. WRKY gene family evolution in Arabidopsis thaliana. Genetica 139, 973 (2011).

[34] Schmidt, M. H. W. et al. De novo assembly of a new Solanum pennellii accession using nanopore sequencing. Plant Cell 29, 2336-2348 (2017).

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

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

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

[38] Stanke, M., Schoffmann, 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).

[39] Salamov, A. A. & Solovyev, V. V. Ab initio gene finding in Drosophila genomic DNA. Genome Res. 10, 516-522 (2000).

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

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

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

[43] Ogata, H. et al. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 27, 29-34 (1999).

[44] Hunter, S. et al. InterPro: the integrative protein signature database. Nucleic Acids Res. 37, D211-215 (2009).

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

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

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

[48] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

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

[50] Zhang, G. Q. et al. The Apostasia genome and the evolution of orchid. Nature 549, 379-383 (2017).

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

[52] Letunic, L., Doerks, T. & Bork, P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, D257-D260 (2014).

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

[54] Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. 2010 Gateway Computing Environments Workshop (GCE), (ed. Institute of Electrical and Electronic Engineers) 1-8 (New Orleans, LA, 2010).