TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Abstract

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N⁶-methyladenosine (m⁶A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

References

[1] Chen F, Song Y, Li X et al. Genome sequences of horticultural plants: past, present, and future. Hortic Res. 2019; 6: 112.

[2] Weng Y. Inaugural editorial: vegetable research. Veg Res. 2021; 1: 1.

[3] Pei Q, Yu T, Wu T et al. Comprehensive identification and analyses of the Hsf gene family in the whole-genome of three Apiaceae species. Hortic Plant J. 2021; 7: 457-68.

[4] Pei Q, Li N, Bai Y et al. Comparative analysis of the TCP gene family in celery, coriander and carrot (family Apiaceae). Vegetable Res. 2021; 1: 5.

[5] Mei Y, Jing D, Tang S et al. InsectBase 2.0: a comprehensive gene resource for insects. Nucleic Acids Res. 2021; 50: D1040-5.

[6] Huang S, Li R, Zhang Z et al. The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009; 41: 1275-81.

[7] Song X, Sun P, Yuan J et al. The celery genome sequence reveals sequential paleo-polyploidizations, karyotype evolution and resistance gene reduction in Apiales. Plant Biotechnol J. 2021; 19: 731-44.

[8] Song X, Wang J, Li N et al. Deciphering the high-quality genome sequence of coriander that causes controversial feelings. Plant Biotechnol J. 2020; 18: 1444-56.

[9] Sun X, Zhu S, Li N et al. A chromosome-level genome assembly of garlic (Allium sativum) provides insights into genome evolution and allicin biosynthesis . Mol Plant. 2020; 13: 1328-39.

[10] Finkers R, van Kaauwen M, Ament K et al. Insights from the first genome assembly of onion (Allium cepa) . G3 Genes|Genomes|Genetics. 2021; 11: jkab243.

[11] Golicz AA, Bayer PE, Barker GC et al. The pangenome of an agronomically important crop plant Brassica oleracea . Nat Commun. 2016; 7: 13390.

[12] He Z, Ji R, Havlickova L et al. Genome structural evolution in Brassica crops. Nat Plants. 2021; 7: 757-65.

[13] Cai X, Chang L, Zhang T et al. Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa . Genome Biol. 2021; 22: 166.

[14] Song JM, Guan Z, Hu J et al. Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus . Nat Plants. 2020; 6: 34-45.

[15] Bayer PE, Scheben A, Golicz AA et al. Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids. Plant Biotechnol J. 2021; 19: 2488-500.

[16] Barchi L, Rabanus-Wallace MT, Prohens J et al. Improved genome assembly and pan-genome provide key insights into eggplant domestication and breeding. Plant J. 2021; 107: 579-96.

[17] Gao L, Gonda I, Sun H et al. The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat Genet. 2019; 51: 1044- 51.

[18] Zhang X, Liu T, Wang J et al. Pan-genome of Raphanus highlights genetic variation and introgression among domesticated, wild, and weedy radishes. Mol Plant. 2021; 14: 2032-55.

[19] Sayers EW, Beck J, Bolton EE et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2021; 49: D10-7.

[20] Song X, Nie F, Chen W et al. Coriander Genomics Database: a genomic, transcriptomic, and metabolic database for coriander. Hortic Res. 2020; 7: 55.

[21] Wu T, Feng S, Yang Q et al. Integration of the metabolome and transcriptome reveals the metabolites and genes related to nutritional and medicinal value in Coriandrum sativum . J Integr Agric. 2021; 20: 1807-18.

[22] Yu HJ, Baek S, Lee YJ et al. The radish genome database (RadishGD): an integrated information resource for radish genomics. Database (Oxford). 2019; 2019: baz009.

[23] Wei Q, Wang J, Wang W et al. A high-quality chromosome-level genome assembly reveals genetics for important traits in eggplant. Hortic Res. 2020; 7: 153.

[24] Reyes-Chin-Wo S, Wang Z, Yang X et al. Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce. Nat Commun. 2017; 8: 14953.

[25] Rodriguez Del Rio A, Minoche AE, Zwickl NF et al. Genomes of the wild beets Beta patula and Beta vulgaris ssp. maritima . Plant J. 2019; 99: 1242-53.

[26] Lehner R, Blazek L, Minoche AE et al. Assembly and characterization of the genome of chard (Beta vulgaris ssp. vulgaris var. cicla) . J Biotechnol. 2021; 333: 67-76.

[27] Garcia-Mas J, Benjak A, Sansevarino W . The genome of melon (Cucumis melo L.) . Proc Natl Acad Sci USA. 2012; 109: 11872-7.

[28] Harkess A, Zhou J, Xu C et al. The asparagus genome sheds light on the origin and evolution of a young Y chromosome. Nat Commun. 2017; 8: 1279.

[29] Yang J, Moeinzadeh MH, Kuhl H et al. Haplotype-resolved sweet potato genome traces back its hexaploidization history. Nature Plants. 2017; 3: 696-703.

[30] Hufnagel B, Marques A, Soriano A et al. High-quality genome sequence of white lupin provides insight into soil exploration and seed quality. Nat Commun. 2020; 11: 492.

[31] Liu D, Yu L, Wei L et al. BnTIR: an online transcriptome platform for exploring RNA-seq libraries for oil crop Brassica napus . Plant Biotechnol J. 2021; 19: 1895-7.

[32] Chen H, Wang T, He X et al. BRAD V3.0: an upgraded Brassicaceae database. Nucleic Acids Res. 2022; 50: D1432-41.

[33] Liu Z, Li N, Yu T et al. The Brassicaceae genome resource (TBGR): a comprehensive genome platform for Brassicaceae plants. Plant Physiol. 2022; 190: 226-37.

[34] Fernandez-Pozo N, Menda N, Edwards JD et al. The Sol Genomics Network (SGN)-from genotype to phenotype to breeding. Nucleic Acids Res. 2015; 43: D1036-41.

[35] Zheng Y, Wu S, Bai Y et al. Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops. Nucleic Acids Res. 2019; 47: D1128-36.

[36] Yang Y, Sun P, Lv L et al. Prickly waterlily and rigid hornwort genomes shed light on early angiosperm evolution. Nat Plants. 2020; 6: 215-22.

[37] Yang J, Liu D, Wang X et al. The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat Genet. 2016; 48: 1225-32.

[38] Chalhoub B, Denoeud F, Liu S et al. Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014; 345: 950-3.

[39] Barrera-Redondo J, Ibarra-Laclette E, Vázquez-Lobo A et al. The genome of Cucurbita argyrosperma (silver-seed gourd) reveals faster rates of protein-coding gene and long noncoding RNA turnover and neofunctionalization within Cucurbita . Mol Plant. 2019; 12: 506-20.

[40] Sun H, Wu S, Zhang G et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid Cucurbita genomes. Mol Plant. 2017; 10: 1293-306.

[41] Montero-Pau J, Blanca J, Bombarely A 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. 2018; 16: 1161-71.

[42] Wang X, Wang H, Wang J et al. The genome of the mesopolyploid crop species Brassica rapa . Nat Genet. 2011; 43: 1035-9.

[43] Cheng SP, Jia KH, Liu H et al. Haplotype-resolved genome assembly and allele-specific gene expression in cultivated ginger. Hortic Res. 2021; 8: 188.

[44] Li HL, Wu L, Dong Z et al. Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway . Hortic Res. 2021; 8: 189.

[45] Song X, Li Y, Hou X . Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis) . BMC Genomics. 2013; 14: 573.

[46] Song X, Nie F, Chen W et al. Coriander genomics database: a genomic, transcriptomic, and metabolic database for coriander. Hortic Res. 2020; 7: 55.

[47] Moharana KC, Venancio TM . Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae). Plant J. 2020; 103: 726-41.

[48] Shirasawa K, Hirakawa H, Fukino N et al. Genome sequence and analysis of a Japanese radish (Raphanus sativus) cultivar named ’Sakurajima Daikon’ possessing giant root . DNA Res. 2020; 27: dsaa010.

[49] Song X, Wei Y, Xiao D et al. Brassica carinata genome characterization clarifies U’s triangle model of evolution and polyploidy in Brassica . Plant Physiol. 2021; 186: 388-406.

[50] Yue H, Nie X, Yan Z et al. N6-methyladenosine regulatory machinery in plants: composition, function and evolution. Plant Biotechnol J. 2019; 17: 1194-208.

[51] UniProt C . UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021; 49: D480-9.

[52] Mistry J, Chuguransky S, Williams L et al. Pfam: the protein families database in 2021. Nucleic Acids Res. 2021; 49: D412-9.

[53] Gene Ontology Consortium . The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021; 49: D325-34.

[54] Song X, Yang Q, Bai Y et al. Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants. Hortic Res. 2021; 8: 122.

[55] Song X, Ge T, Li Y et al. Genome-wide identification of SSR and SNP markers from the non-heading Chinese cabbage for comparative genomic analyses. BMC Genomics. 2015; 16: 328.

[56] Song X, Li N, Guo Y et al. Comprehensive identification and characterization of simple sequence repeats based on the whole-genome sequences of 14 forest and fruit trees. For Res. 2021; 1: 1-10.

[57] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 238.

[58] Wang Y, Tang H, Debarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49.

[59] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[60] Li H, Fan Y, Yu J et al. Genome-wide identification of flowering-time genes in Brassica species and reveals a correlation between selective pressure and expression patterns of vernalization-pathway genes in Brassica napus . Int J Mol Sci. 2018; 19: E3632.

[61] Bouche F, Lobet G, Tocquin P et al. FLOR-ID: an interactive database of flowering-time gene networks in Arabidopsis thaliana . Nucleic Acids Res. 2016; 44: D1167-71.

[62] Aach J, Mali P, Church GM . CasFinder: flexible algorithm for identifying specific Cas9 targets in genomes. bioRxiv. 2014; 005074.

[63] Wang H, Wu J, Sun S et al. Glucosinolate biosynthetic genes in Brassica rapa . Gene. 2011; 487: 135-42.

[64] Cheng F, Wu J, Wang X . Genome triplication drove the diversification of Brassica plants. Hortic Res. 2014; 1: 14024.

[65] Li P, Quan X, Jia G et al. RGAugury: a pipeline for genome-wide prediction of resistance gene analogs (RGAs) in plants. BMC Genomics. 2016; 17: 852.

[66] Sekhwal MK, Li P, Lam I et al. Disease resistance gene analogs (RGAs) in plants. Int J Mol Sci. 2015; 16: 19248-90.

[67] Kozomara A, Birgaoanu M, Griffiths-Jones S . miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019; 47: D155-62.

[68] Chen C, Li J, Feng J et al. sRNAanno-a database repository of uniformly annotated small RNAs in plants. Hortic Res. 2021; 8: 45.

[69] Lorenz R, Bernhart SH, Siederdissen HZ et al. ViennaRNA package 2.0. Algorithms Mol Biol. 2011; 6: 26.

[70] Kielbasa SM, Bluthgen N, Fahling M et al. Targetfinder.org: a resource for systematic discovery of transcription factor target genes. Nucleic Acids Res. 2010; 38: W233-8.

[71] Tarailo-Graovac M, Chen N . Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009; Chapter 4, Unit 4.10;25.

[72] Giannoulatou E, Park SH, Humphreys DT et al. Verification and validation of bioinformatics software without a gold standard: a case study of BWA and bowtie. Bioinformatics. 2014; 15: S15.

[73] Yu T, Bai Y, Liu Z et al. Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants. Hortic Res. 2022; 9: uhac035.