A high-quality sponge gourd (Luffa cylindrica) genome

Abstract

Sponge gourd (Luffa cylindrica) is an important cultivated vegetable and medicinal plant in the family Cucurbitaceae. In this study, a draft genome sequence of the sponge gourd inbred line P93075 was analyzed. Using Illumina, PacBio, and 10× Genomics sequencing techniques as well as new assembly techniques such as FALCON and chromatin interaction mapping (Hi-C), a chromosome-scale genome of approximately 656.19 Mb, with an N50 scaffold length of 48.76 Mb, was generated. From this assembly, 25,508 protein-coding gene loci were identified, and 63.81% of the whole-genome consisted of transposable elements, which are major contributors to the expansion of the sponge gourd genome. According to a phylogenetic analysis of conserved genes, the sponge gourd lineage diverged from the bitter gourd lineage approximately 41.6 million years ago. Additionally, many genes that respond to biotic and abiotic stresses were found to be lineage specific or expanded in the sponge gourd genome, as demonstrated by the presence of 462 NBS-LRR genes, a much greater number than are found in the genomes of other cucurbit species; these results are consistent with the high stress resistance of sponge gourd. Collectively, our study provides insights into genome evolution and serves as a valuable reference for the genetic improvement of sponge gourd.

References

[1] Sohrab, S. S., Mandal, B., Pant, R. P. & Varma, A. First report of association of tomato leaf curl virus-New Delhi with yellow mosaic disease of Luffa cylindrica in India. Plant Dis. 87, 1148- 1148 (2003).

[2] Zhang, L., Shi, M., Ji, J., Hu, X. & Chen, F. Gut microbiota determines the prevention effects of Luffa cylindrica (L.) Roem supplementation against obesity and associated metabolic disorders induced by high-fat diet. FASEB J. 33, 10339- 10339 (2019).

[3] Azeez, M. A., Bello, O. S. & Adedeji, A. O. Traditional and medicinal uses of Luffa cylindrica: a review. J. Med. Plants Stud. 5, 102-111 (2013).

[4] Oboh, I. O. & Aluyor, E. O. Luffa cylindrica-an emerging cash crop. Afr. J. Agr. Res. 4, 684-688 (2009).

[5] Dairo, F. A. S., Aye, P. A. & Oluwasola, T. A. Some functional properties of loofah gourd (Luffa cylindrica L., M. J. Roem) seed. J. Food, Agric. Environ. 5, 97-101 (2007).

[6] Arise, R. O. et al. In vitro Angiotesin-1-converting enzyme, α-amylase and α-glucosidase inhibitory and antioxidant activities of Luffa cylindrical (L.)M. Roem seed protein hydrolysate. Heliyon 5, e01634 (2019).

[7] Abdel-Salam, I. M., Awadein, N. E. S. & Ashour, M. Cytotoxicity of Luffa cylindrica (L.) M.Roem. extract against circulating cancer stem cells in hepatocellular carcinoma. J. Ethnopharmacol. 229, 89-96 (2019).

[8] Abdel-Salam, I. M., Abou-Bakr, A. A. & Ashour, M. Cytotoxic effect of aqueous ethanolic extract of Luffa cylindrica leaves on cancer stem cells CD44+/24- in breast cancer patients with various molecular sub-types using tissue samples in vitro. J. Ethnopharmacol. 238, 111877 (2019).

[9] Khajuria, A., Gupta, A., Garai, S. & Wakhloo, B. P. Immunomodulatory effects of two sapogenins 1 and 2 isolated from Luffa cylindrica in Balb/C mice. Bioorg. Med. Chem. Lett. 17, 1608-1612 (2007).

[10] Wu, H. et al. Genetic linkage map construction and QTL analysis of two interspecific reproductive isolation traits in sponge gourd. Front. Plant Sci. 7, 980 (2016).

[11] Wu, H. et al. Large-scale development of EST-SSR markers in sponge gourd via transcriptome sequencing. Mol. Breeding 34, 1903-1915 (2014).

[12] Xie, D. et al. The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype. Nat. Commun. 10, 5158 (2019).

[13] Xu, Y. et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat. Genet. 45, 51-58 (2013).

[14] Garcia-Mas, J. et al. The genome of melon (Cucumis melo L.). Proc. Natl Acad. Sci. USA 109, 11872-11877 (2012).

[15] Huang, S. et al. The genome of the cucumber, Cucumis sativus L. Nat. Genet. 41, 1275-1281 (2009).

[16] Sun, H. et al. Karyotype stability and unbiased fractionation in the Paleo-Allotetraploid Cucurbita genomes. Mol. Plant 10, 1293-1306 (2017).

[17] Montero-Pau, J. et al. De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol. J. 16, 1161-1171 (2018).

[18] Peng, X. et al. Complete genome sequence of a new strain of Lagenaria siceraria endornavirus from China. Arch. Virol. 163, 805-808 (2018).

[19] Urasaki, N. et al. Draft genome sequence of bitter gourd (Momordica charantia), a vegetable and medicinal plant in tropical and subtropical regions. DNA Res. 24, 51-58 (2017).

[20] Li, H. & Wren, J. Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics 30, 2843-2851 (2014).

[21] Seppey, M., Manni, M. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness. Methods Mol. Biol. 1962, 227-245 (2019).

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

[23] Wang, J. et al. An overlooked paleotetraploidization in cucurbitaceae. Mol. Biol. Evol. 35, 16-26 (2018).

[24] Li, J. et al. Unique evolutionary pattern of numbers of gramineous NBS-LRR genes. Mol. Genet. Genomics 283, 427-438 (2010).

[25] Zhang, Y. M. et al. Uncovering the dynamic evolution of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes in Brassicaceae. J. Integr. Plant Biol. 58, 165-177 (2016).

[26] Zaafouri, K. et al. Enzymatic saccharification and liquid state fermentation of hydrothermal pretreated Tunisian Luffa cylindrica (L.) fibers for cellulosic bioethanol production. Renew. Energy 114, 1209-1213 (2017).

[27] Elisa, D., Heinz, E., Branislav, K., Petra, B. & Karel, R. Role of STN1 and DNA polymerase α in telomere stability and genome-wide replication in Arabidopsis. PLos Genet. 10, e1004682 (2014).

[28] Liao, C. T. & Lin, C. H. Photosynthetic responses of grafted bitter melon seedlings to flood stress. Environ. Exp. Bot. 2, 167-172 (1996).

[29] Lin, K. H. et al. Study of sponge gourd ascorbate peroxidase and winter squash superoxide dismutase under respective flooding and chilling stresses. Sci. Hortic. (Amst.) 162, 333-340 (2013).

[30] Parkash, A., Ng, T. B. & Tso, W. W. Isolation and characterization of luffacylin, a ribosome inactivating peptide with anti-fungal activity from sponge gourd (Luffa cylindrica) seeds. Peptides 23, 1019-1024 (2002).

[31] Fink, J. R. et al. Diffusion and uptake of phosphorus, and root development of corn seedlings, in three contrasting subtropical soils under conventional tillage or no-tillage. Biol. Fertil. Soils 52, 203-210 (2016).

[32] Zhao, Y., Wang, T., Zhang, W. & Li, X. SOS3 mediates lateral root development under low salt stress through regulation of auxin redistribution and maxima in Arabidopsis. N. Phytol. 189, 1122-1134 (2011).

[33] Xia, J., Yamaji, N. & Jian, F. M. An appropriate concentration of arginine is required for normal root growth in rice. Plant Signal. Behav. 9, e28717 (2014).

[34] Gad, M. Z. Anti-aging effects of l-arginine. J. Adv. Res. 1, 169-177 (2010).

[35] Kasinathan, V. & Wingler, A. Effect of reduced arginine decarboxylase activity on salt tolerance and on polyamine formation during salt stress in Arabidopsis thaliana. Physiol. Plant. 121, 101-107 (2010).

[36] Fernando, L. N. & Grün, I. U. Headspace-SPME analysis of volatiles of the ridge gourd (Luffa acutangula) and bitter gourd (Momordica charantia) flowers. Flavour Fragr. J. 16, 289-293 (2010).

[37] Peng, F., Meng, C. W., Zhou, Q. M., Chen, J. P. & Xiong, L. Cytotoxic evaluation against breast cancer cells of isoliquiritigenin analogues from spatholobus suberectus and their synthetic derivatives. J. Nat. Prod. 79, 248-251 (2016).

[38] Wang, H., Liu, Y., Zenc, Z. & He, W. Study on HPLC chromatographic fingerprint of anti-tumor active site SSCE of Caulis spatholobi. Zhongguo Zhongyao Zazhi 36, 2525-2529 (2011).

[39] Zhou, Z. Y. et al. Spatholobi Caulis extracts promote angiogenesis in HUVECs in vitro and in zebrafish embryos in vivo via up-regulation of VEGFRs. J. Ethnopharmacol. 200, 74-83 (2017).

[40] Bowerman, P. A., Ramirez, M. V., Price, M. B., Helm, R. F. & Winkel, B. S. J. Analysis of T-DNA alleles of flavonoid biosynthesis genes in Arabidopsis ecotype Columbia. BMC Res. Notes 5, 485 (2012).

[41] Saito, K. et al. The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol. Biochem. 72, 21-34 (2013).

[42] Aoki, T., Akashi, T. & Ayabe, S. Flavonoids of leguminous plants: structure, biological activity, and biosynthesis. J. Plant Res. 113, 475-488 (2000).

[43] Langfelder, P. & Horvath, S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9, 559 (2008).

[44] Sánchez-Rodríguez, C. et al. CHITINASE-LIKE1/POM-POM1 and its homolog CTL2 are glucan-lnteracting proteins important for cellulose biosynthesis in Arabidopsis. Plant Cell 24, 589-607 (2012).

[45] Vanholme, R., Demedts, B., Morreel, K., Ralph, J. & Boerjan, W. Lignin biosynthesis and structure. Plant Physiol. 153, 895-905 (2010).

[46] Eitan, Y. & Amos, T. Probabilistic modeling of Hi-C contact maps eliminates systematic biases to characterize global chromosomal architecture. Nat. Genet. 43, 1059 (2011).

[47] Liu, B. et al. Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects. http://arxiv.org/abs/1308.2012v1. (2013).

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

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

[50] Smit, A. & Hubley, R. RepeatModeler Open-1.0. GItHub (2018).

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

[52] Tarailo-Graovac, M. & Chen, N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Prot. Bioinform. 25, 4.10.1-4.10.14 (2009).

[53] Kim, D. et al. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

[54] Ghosh, S. & Chan, C. K. Analysis of RNA-Seq data using TopHat and Cufflinks. Methods Mol. Biol. 1374, 339-361 (2016).

[55] Haas, B. J. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 31, 5654-5666 (2003).

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

[57] Finn, R. D. et al. The Pfam protein families database. Nucleic Acids Res. 38, D211-D222 (2010).

[58] Finn, R. D. et al. HMMER web server: 2015 update. Nucleic Acids Res. 43, 30-38 (2015).

[59] Finn, R. D. et al. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res. 45, D190-D199 (2017).

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

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

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

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

[64] Wang, X. et al. Statistical inference of chromosomal homology based on gene colinearity and applications to Arabidopsis and rice. BMC Bioinformatics 7, 447 (2006).

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

[66] Edgar, R. C. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5, 113 (2004).