MODMS: a multi-omics database for facilitating biological studies on alfalfa (Medicago sativa L.)

Abstract

Alfalfa (Medicago sativa L.) is a globally important forage crop. It also serves as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Multi-omics data on alfalfa continue to accumulate owing to recent advances in high-throughput techniques, and integrating this information holds great potential for expediting genetic research and facilitating advances in alfalfa agronomic traits. Therefore, we developed a comprehensive database named MODMS (multi-omics database of M. sativa) that incorporates multiple reference genomes, annotations, comparative genomics, transcriptomes, high-quality genomic variants, proteomics, and metabolomics. This report describes our continuously evolving database, which provides researchers with several convenient tools and extensive omics data resources, facilitating the expansion of alfalfa research. Further details regarding the MODMS database are available at https://modms.lzu.edu.cn/.

Acknowledgements

This work was partly supported by the National Key Research and Development Program (2022YFF1003200), and the National Natural Science Foundation of China (Grant No. 32201463 and 32201444). We received support from the Supercomputing Center of Lanzhou University.

Author contributions

L.F. and Z.L. designed the project, L.F., T.L., M.L., and X.D. per-formed data collection, data analysis, and database construction. J.Z., S.L., C.X., and X.H. performed the bioinformatics analysis. L.F., T.L., Y.H., and Z.L. wrote the manuscript. L.F. and Z.L. supervised the project. All authors read and approved the final manuscript. The authors declare no competing interests.

Data availability

All datasets have been made publicly available at https://modms.lzu.edu.cn/, which is accessible to all users free of charge and without requiring a login.

Conflict of interest

The authors declare that they have no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Feng Y, Shi Y, Zhao M. et al. Yield and quality properties of alfalfa (Medicago sativa L.) and their influencing factors in China. Eur J Agron. 2022;141:126637

[2] Dhakal M, West CP, Villalobos C. et al. Trade-off between nutri-tive value improvement and crop water use for an alfalfa-grass system. Crop Sci. 2020;60:1711-23

[3] Acharya JP, Lopez Y, Gouveia BT. et al. Breeding alfalfa (Med-icago sativa L.) adapted to subtropical agroecosystems. Agronomy. 2020;10:742

[4] Zhou Q, Jia C, Ma W. et al. MYB transcription factors in alfalfa (Medicago sativa): genome-wide identification and expression analysis under abiotic stresses. PeerJ. 2019;7:e7714

[5] He F, Zhang F, Jiang X. et al. A genome-wide association study coupled with a transcriptomic analysis reveals the genetic loci and candidate genes governing the flowering time in alfalfa (Medicago sativa L.). Front Plant Sci. 2022;13:913947

[6] Wang Q, Zhang Y, Hu H. et al. Comparative proteomics analysis of the responses to selenium in selenium-enriched alfalfa (Med-icago sativa L.) leaves. Plant Physiol Biochem. 2021;165:265-73

[7] Song T, Xu H, Sun N. et al. Metabolomic analysis of alfalfa (Med-icago sativa L.) root-symbiotic rhizobia responses under alkali stress. Front Plant Sci. 2017;8:1208

[8] Luo D, Zhou Q, Wu Y. et al. Full-length transcript sequenc-ing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.). BMC Plant Biol. 2019;19:1-20

[9] Chen H, Zeng Y, Yang Y. et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat Commun. 2020;11:2494

[10] Shen C, du H, Chen Z. et al. The chromosome-level genome sequence of the autotetraploid alfalfa and resequencing of core germplasms provide genomic resources for alfalfa research. Mol Plant. 2020;13:1250-61

[11] Chen L, He F, Long R. et al. A global alfalfa diversity panel reveals genomic selection signatures in Chinese varieties and genomic associations with root development. J Integr Plant Biol. 2021;63: 1937-51

[12] Li A, Liu A, Du X. et al. A chromosome-scale genome assembly of a diploid alfalfa, the progenitor of autotetraploid alfalfa. Hortic Res. 2020;7:194

[13] Zhang JJ, Lu YC, Zhang SH. et al. Identification of transcriptome involved in atrazine detoxification and degradation in alfalfa (Medicago sativa) exposed to realistic environmental contamina-tion. Ecotoxicol Environ Saf. 2016;130:103-12

[14] Liu W, Zhang Z, Chen S. et al. Global transcriptome profiling analysis reveals insight into saliva-responsive genes in alfalfa. Plant Cell Rep. 2016;35:561-71

[15] Wang J, Zhao Y, Ray I. et al. Transcriptome responses in alfalfa associated with tolerance to intensive animal grazing. Sci Rep. 2016;6:1-18

[16] Li Y, Nan Z, Duan T. Rhizophagus intraradices promotes alfalfa (Medicago sativa) defense against pea aphids (Acyrthosiphon pisum) revealed by RNA-seq analysis. Mycorrhiza. 2019;29: 623-35

[17] Dong W, Liu X, Li D. et al. Transcriptional profiling reveals that a MYB transcription factor MsMYB4 contributes to the salinity stress response of alfalfa. PLoS One. 2018;13:e0204033

[18] Postnikova OA, Shao J, Nemchinov LG. Analysis of the alfalfa root transcriptome in response to salinity stress. Plant Cell Physiol. 2013;54:1041-55

[19] Cui W, Yao P, Pan J. et al. Transcriptome analysis reveals insight into molecular hydrogen-induced cadmium tolerance in alfalfa: the prominent role of sulfur and (homo) glutathione metabolism. BMC Plant Biol. 2020;20:1-19

[20] Zhou Q, Luo D, Chai X. et al. Multiple regulatory networks are activated during cold stress in Medicago sativa L. Int J Mol Sci. 2018;19:3169

[21] Postnikova OA, Hult M, Shao J. et al. Transcriptome analy-sis of resistant and susceptible alfalfa cultivars infected with root-knot nematode Meloidogyne incognita. PLoS One. 2015;10: e0118269

[22] Cheng Q, Bai S, Ge G. et al. Study on differentially expressed genes related to defoliation traits in two alfalfa varieties based on RNA-Seq. BMC Genomics. 2018;19:1-8

[23] Zhang S, Shi Y, Cheng N. et al. De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa L.) leaf transcriptome and identification of candidate genes related to fall dormancy. PLoS One. 2015;10:e0122170

[24] Li X, Acharya A, Farmer AD. et al. Prevalence of single nucleotide polymorphism among 27 diverse alfalfa genotypes as assessed by transcriptome sequencing. BMC Genomics. 2012;13:1-11

[25] Gruber M, Xia J, Yu M. et al. Transcript analysis in two alfalfa salt tolerance selected breeding populations relative to a non-tolerant population. Genome. 2017;60:104-27

[26] O’Rourke JA, Fu F, Bucciarelli B. et al. The Medicago sativa gene index 1.2: a web-accessible gene expression atlas for investigat-ing expression differences between Medicago sativa subspecies. BMC Genomics. 2015;16:1-17

[27] Wang Y, Zhang P, Li L. et al. Proteomic analysis of alfalfa (Med-icago sativa L.) roots in response to rhizobium nodulation and salt stress. Genes. 2022;13:2004

[28] Chen L, Xia F, Wang M. et al. Metabolomic analyses of alfalfa (Medicago sativa L. cv.‘Aohan’) reproductive organs under boron deficiency and surplus conditions. Ecotoxicol Environ Saf. 2020;202:111011

[29] Yu T, Ma X, Liu Z. et al. TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies. Hortic Res. 2022;9:uhac213

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

[31] Li Q, Qi J, Qin X. et al. CitGVD: a comprehensive database of citrus genomic variations. Hortic Res. 2020;7:12

[32] Yang Z, Wang J, Huang Y. et al. CottonMD: a multi-omics database for cotton biological study. Nucleic Acids Res. 2023;51:D1446-56

[33] Zhang Z, Chai M, Yang Z. et al. GRAND: an integrated genome, transcriptome resources, and gene network database for gossyp-ium. Front Plant Sci. 2022;13:24

[34] Kataria R, Kaundal R. alfaNET: a database of alfalfa-bacterial stem blight protein-protein interactions revealing the molec-ular features of the disease-causing bacteria. Int J Mol Sci. 2021;22:8342

[35] Carrere S, Verdier J, Gamas P. MtExpress, a comprehensive and curated RNAseq-based gene expression atlas for the model legume Medicago truncatula. Plant Cell Physiol. 2021;62: 1494-500

[36] Roux B, Rodde N, Jardinaud MF. et al. An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing. Plant J. 2014;77:817-37

[37] Carrere S, Verdenaud M, Gough C. et al. LeGOO: an expertized knowledge database for the model legume Medicago truncatula. Plant Cell Physiol. 2020;61:203-11

[38] Pecrix Y, Staton SE, Sallet E. et al. Whole-genome landscape of Medicago truncatula symbiotic genes. Nat Plants. 2018;4:1017-25

[39] Buels R, Yao E, Diesh CM. et al. JBrowse: a dynamic web platform for genome visualization and analysis. Genome Biol. 2016;17: 1-12

[40] Camacho C, Coulouris G, Avagyan V. et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:1-9

[41] Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792-7

[42] Dewey CN. Whole-genome alignment. Methods Mol Biol. 2012;855:237-57

[43] Birney E, Clamp M, Durbin R. GeneWise and genomewise. Genome Res. 2004;14:988-95

[44] Kieber JJ, Schaller GE. Cytokinin signaling in plant development. Development. 2018;145:dev149344

[45] Heckmann AB, Sandal N, Bek AS. et al. Cytokinin induction of root nodule primordia in Lotus japonicus is regulated by a mech-anism operating in the root cortex. Mol Plant-Microbe Interact. 2011;24:1385-95

[46] Frébortová J, Novák O, Frébort I. et al. Degradation of cytokinins by maize cytokinin dehydrogenase is mediated by free radicals generated by enzymatic oxidation of natural benzoxazinones. Plant J. 2010;61:467-81

[47] Young ND, Debellé F, Oldroyd GED. et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011;480:520-4

[48] Bu D, Luo H, Huo P. et al. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021;49:W317-25

[49] Huerta-Cepas J, Szklarczyk D, Heller D. et al. eggNOG 5.0: a hier-archical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47:D309-14

[50] Liao Y, Smyth GK, Shi W. featureCounts: an efficient general pur-pose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923-30

[51] Pei G, Chen L, Zhang W. WGCNA application to proteomic and metabolomic data analysis. Methods Enzymol. 2017;585: 135-58

[52] Su G, Morris JH, Demchak B. et al. Biological network exploration with Cytoscape 3. Curr Protoc Bioinformatics. 2014;47:1-24

[53] Mahram A, Herbordt MC. NCBI BLASTP on high-performance reconfigurable computing systems. ACM Trans Reconfigurable Technol Syst. 2015;7:1-20

[54] 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-9

[55] Li H, Handsaker B, Wysoker A. et al. The sequence alignmen-t/map format and SAMtools. Bioinformatics. 2009;25:2078-9

[56] McKenna A, Hanna M, Banks E. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297-303

[57] Danecek P, Auton A, Abecasis G. et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156-8